Shock and Vibration
Volume 2021 (2021), Article ID 5581398, 10 pages
https://doi.org/10.1155/2021/5581398

Research Article

Coupled Vibration Characteristics Analysis of Hot Rolling Mill with Structural Gap

Guangxu Zhang,1,2 Jiahan Bao,1,2 Wenhao Li,1,2 Zhichong Wang,1,2 and Xiangshuai Meng1,2

¹Engineering Research Center of Hydraulic Vibration and Control, Ministry of Education, Anhui University of Technology, Maanshan, Anhui, China
²School of Mechanical Engineering, Anhui University of Technology, Maanshan, Anhui, China

Correspondence should be addressed to Jiahan Bao; baojiahan@126.com

Received 1 February 2021; Revised 6 June 2021; Accepted 9 June 2021; Published 19 June 2021

Academic Editor: Ömer Cívalek

Copyright © 2021 Guangxu Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

It is important to study the vibration of rolling mills to improve the stability of rolling production. A dynamic rolling process model is established by considering the elastic recovery of the exit strip and the influence of multiroll equilibrium, and the accuracy of the model is verified by experimental data. On this basis, based on the distribution of friction force in the deformation zone, the rolling force and rolling torque are nonlinearized. In addition, a rolling mill structure model is established by considering the structure gap and a piecewise nonlinear horizontal-vertical-torsional vibration model of the rolling mill is established by combining the structure model and dynamic rolling process model. Finally, the amplitude-frequency characteristics of the work roll under different external excitation amplitude and the dynamic bifurcation characteristics of the work roll under different gaps are analyzed. The study indicates that, by reducing excitation amplitude and structure gap, the system vibration can be reduced. The research results can provide a theoretical reference for further exploration of the coupling vibration of hot rolling mills.

1. Introduction

In the rolling process, rolling mill vibration is an inevitable problem, which causes vibration of steel strips and wears the parts of rolling mills [1–3]. Therefore, studying rolling mill vibration and put forward vibration suppression methods are beneficial to improve the service life of rolling mill parts and surface quality of rolled pieces.

The study of rolling mill vibration can be regarded as the study of rolling process and rolling mill structure. Most studies on rolling process focus on the roll gap [4, 5]. For instance, Li et al. [6], Chen et al. [7], and Heidari et al. [8] established the rolling process model by considering the friction on the roll gap and studied horizontal rolling force, vertical rolling force, and neutral angle. Cao et al. [9] and Ma et al. [10] calculated the rolling force and torque by using the functional relationship. Son et al. [11] established the rolling force model by using the neural network, and the relationship between rolling parameters and rolling force was studied without the modeling analysis of rolling deformation zone. Most studies on rolling mill structure focus on characteristics of horizontal vibration [12], vertical vibration [13], torsional vibration [14], and coupling vibration of rolling mills [15–17]. Wang et al. [18] established a coupled vibration model of rolling mill by considering the offset of the work roll. Peng et al. [19] established a horizontal-vertical vibration model of rolling mill by considering nonlinear hydraulic pressure, nonlinear friction force, and nonlinear rolling force. Zeng et al. [20, 21] established a horizontal-vertical-torsional vibration model of rolling mill by considering the relationship between dynamic rolling conditions and damping coefficient.

However, with the abovementioned study, elastic recovery of exit strip in the rolling process and the gap between work roll bearing seat and mill stand in the rolling mill structure were ignored. There are still some limitations or problems with existing models and studies. Therefore, a more complete coupled vibration model of rolling mill needs to be proposed and analyzed.

In this study, we consider the effects of the elastic recovery of exit strip, the force of backup roll on work roll, and the parameter coupling of horizontal, vertical, and torsional directions and establish a dynamic rolling process model. Then, we consider the gap between work roll bearing seat and mill stand and the dynamic rolling process and establish a dynamic nonlinear coupled vibration model of rolling mill. Finally, the bifurcation and amplitude-frequency characteristics of the coupled vibration system of the rolling mill under different conditions are studied.

2. Dynamic Rolling Process Model

2.1. Dynamic Rolling Force

It is an ideal rolling process to keep the rolling parameters unchanged and the rolling mill has no vibration. The variable rolling parameters are the excitation source of rolling mill vibration. When the entrance strip thickness changes, the rolling mill will have coupling vibration and it can characterize the rolling continuity [22, 23]. Therefore, in this paper, the thickness change of the entry strip is an external excitation.

The vibration geometry model of the roll gap is presented by considering elastic recovery in Figure 1. is the unit rolling force. f₁, f₂, and f₃ are the frictions in the deformation zone. is the flattened work roll radius. Hin, Hn, and Hout are the strip thicknesses at the entry point, thinnest point, and exit point. l₁ and l₂ are the contact arc lengths in the deformation zone and elastic zone. , , and are the bite angle, neutral angle, and exit angle. and are the strip speeds at the entry and exit points. is the rotational speed of the work roll. and are the front tension and back tension. x₃ and y₃ are the dynamic vibration displacements of work roll in the horizontal direction and vertical direction.

Figure 1: Vibration geometry model of the roll gap.

Hitchcock indicated that the flattened work roll radius can be evaluated as [24]

where is the Poisson’s ratio; in this research, . E₁ is the Young’s modulus of the work roll; in this research, E₁ = 2.1 × 10⁵. Bm is the average width of the strip; in this research, Bm = 1178 mm. P is the total rolling force.

The contact arc length in the deformation zone can be calculated by

The contact arc length in the elastic zone is related to the variation of average exit strip thickness:

The strip thickness at the entry point can be calculated by

where is half of the strip thickness change at the entry point.

The strip thickness at the exit point can be calculated by

The strip thickness at the thinnest point can be calculated by

where Qhp is the strip elasticity coefficient; hp is the original roll gap spacing. Because of the rolling reaction force, the strip thickness at the thinnest point will be changed, .

The rolling speed at the neutral point can be calculated by

where is the angular velocity; is the vibration velocity in the horizontal direction; is the vibration velocity in the vertical direction.

The friction coefficient is related to roll roughness, rolling speed, roll wear, and lubrication, so the friction coefficient is calculated by [25]

where is the basic friction coefficient; is the speed influence coefficient of friction; is the rolling speed; is the work roll roughness coefficient; r and r₀ are actual work roll roughness and roughness of unused work roll; is the friction influence coefficient; L and L₀ are the actual rolling length and datum rolling length.

Friction is divided into three parts and the vertical force is balanced:

From (9), the neutral angle can be given as

where and can be calculated by

The rolling force can be calculated by [2]

where is the average value of tension:

T.-M. Zhi indicated that the unit rolling force p can be evaluated as follows [26]:

where

where t is the rolling temperature; c is the carbon content (%) in the strip; in this research, c = 0.0925%.

2.2. Dynamic Rolling Torque

The force of backup roll on work roll is presented in Figure 2. O₁ is the center of the work roll, and O₂ is the center of the backup roll. a, b, and m are the arms of force. FD is the force of backup roll on work roll. e is the horizontal distance of the work roll from the backup roll. RD is the backup roll radius. a and b are arms of force. is the rolling force-deflection angle. is the force-deflection angle of backup roll on work roll. is the work roll deflection angle.

Figure 2: Force of backup roll on work roll.

The Rolling torque is calculated by [26]

where

In this research, m = 0.1∼0.3 mm. is the friction coefficient of bearing; in this research, .

2.3. Data Validation

The rolling conditions in this research are shown in Table 1. The comparison between the actual and calculated parameters is presented in Figure 3.

Table 1: Rolling conditions.

Parameters
1 stand
2 stand
3 stand
4 stand
5 stand

Hin (mm)
34.89
17.27
8.63
5.41
3.59

Hn (mm)
15.86
8.32
5.29
3.56
2.52

Hout (mm)
17.27
8.63
5.41
3.59
2.53

(m/s)
1.05
2.50
4.05
6.15
8.79

(MPa)
0
9.05
9.97
12.39
14.57

(MPa)
9.05
9.97
12.39
14.57
16.55

t (°C)
985.60
967.74
951.80
931.37
916.92

P (N)
2×71×10⁷
2×44×10⁷
2×39×10⁷
2×24×10⁷
1×79×10⁷

MS (kN·m)
1914.96
1170.91
650.22
437.64
258.82

(a)

(b)

Figure 3: Comparison between the actual parameters and calculated parameters: (a) rolling force; (b) rolling torque.

As shown in Figure 3, the maximum relative error between the actual rolling force and the calculated rolling force with l₂ is less than 5%. The maximum relative error between the actual rolling torque and the calculated rolling torque with l₂ is less than 10%. The results show that this model has high accuracy. Without l₂, the maximum relative error between the actual rolling force and calculated rolling force is 15%, and the maximum relative error between the actual rolling torque and calculated rolling torque is 25%. The results show that the elastic recovery of the exit strip has more influence on rolling torque than rolling force.

2.4. Processing of Dynamic Rolling Parameters

The coupled vibration of rolling mills shows that the rolling parameters in horizontal, vertical, and torsional directions are coupled.

The horizontal force can be calculated by

Based on Taylor expansion, the strip thickness and rolling speed are separated, and the rolling force is nonlinear in the vertical direction.

where a₁–a₆ are the stiffness coefficients in the vertical direction; a₇–a₉ are the damping coefficients in the vertical direction.

ΔPx and ΔPy are much smaller than PSx and PSy and can be ignored. The deviation angle between the horizontal rolling force and rolling force can be calculated by

Based on (21) and Taylor expansion, the nonlinear horizontal rolling force can be expressed as

where b₁–b₆ are the stiffness coefficients in the horizontal direction; b₇–b₉ are the damping coefficients in the horizontal direction.

Based on Taylor expansion, the nonlinear rolling torque can be written as

where c₁–c₃ are the stiffness coefficients in the torsion direction; c₄–c₆ are the damping coefficients in the torsion direction.

In this study, P, PSy, PSx, and MS are variable mechanical parameters. Px, Py, and M are the mechanical parameters in the steady state of the vibration system of the hot rolling mill. ΔPx, ΔPy, and ΔM are the dynamic values of the mechanical parameters in the unsteady state. Therefore, equations (19), (21), and (22) are dynamic nonlinear equations and the rolling process model has a lot of nonlinear and coupling factors.

3. Nonlinear Vibration Model

Figure 4 is the nonlinear horizontal-vertical-torsional coupled vibration model considering the gap between work roll bearing seat and mill stand, friction on the stand, and nonlinearity of backup roll.

Figure 4: Geometry of the vertical-horizontal-torsional vibration model.

As shown in Figure 4, M₁, M₂, and M₃ are the equivalent mass of the upper part of mill stand, backup roll, and work roll. J₁ and J₂ are the rotational inertia of electrical machinery and upper rolls system of the rolling mill. K₁, K₂, K₃, and K₄ are the equivalent vertical stiffness. Based on Duffing theory, is the equivalent horizontal stiffness of backup roll. In this research, K1x = |b₆|. K3x is the equivalent horizontal stiffness of the work roll. K1z and Cpy are equivalent torsional stiffness and equivalent vertical damping. Ff2 and Ff3 are the frictions on the mill stand. e₁ is the gap between the work roll bearing seat and the mill stand. xi, yi, and are horizontal, vertical, and torsional displacements.

The equivalent horizontal stiffness K23x is calculated by

where is the roll flattening. Hertz indicated that the roll flattening can be evaluated as

The equivalent vertical damping between the work roll and strip can be calculated by

In this research, torsional vibration is considered to be caused by horizontal vibration.

According to the Lagrange theory, the horizontal-vertical-torsional coupled vibration equation can be written as

The gap between the work roll bearing seat and mill stand and the friction are expressed by piecewise functions:

Table 2 shows the vibration parameters of the 1580 hot rolling mill [27].

Table 2: Structural and technological parameters of the 1580 hot rolling mill.

Parameter
Value

M₁ (kg)
94300

M₂ (kg)
50150

M₃ (kg)
28622

J₁ (kg·m²)
10500

J₂ (kg·m²)
1310

K₁ (N/m)
5.21 × 10¹⁰

K₂ (N/m)
2.89 × 10¹⁰

K₃ (N/m)
6.84 × 10¹⁰

K2x (N/m)
1.82 × 10⁹

K3x (N/m)
2.5 × 10⁹

K1z (N/m)
2.5 × 10⁷

e (mm)
10

4. Analysis of Vibration Stability

4.1. Analysis of Amplitude-Frequency Characteristics

In this research, , where . We use the average method to solve the dynamic response of the system [28]. Equation (29) has periodic solutions:

where

The phase angle and displacement of (29) change weakly.

From (31), the relationship between phase angle and displacement can be given as

Substituting equations (29)–(33) into (27), and the vibration is averaged, we obtain

where are the functions of horizontal, vertical, and torsional vibration.

We think , and amplitude-frequency characteristics equations of the rolling mill can be obtained.

The amplitude-frequency characteristics of the work roll under different thickness variations of the entry strip are shown in Figure 5. In this paper, H₁ could be considered an external excitation. The amplitudes of horizontal, vertical, and torsional vibrations of the system become bigger with the increase of the amplitude of external excitation. When H₁ = 0 mm, the rolling mill does not vibrate and this shows that the system is stable without excitation. When 38 Hz < f < 52 Hz, the resonant amplitudes of the system in three directions have multiple solutions under H₁ = 2 mm and this shows that the system is unstable. The resonant zones of the system in three directions are the same under H₁ = 2 mm and 70 Hz < f < 85 Hz, and this shows that the system has two horizontal-vertical-torsional coupling resonance zones under 38 Hz < f < 85 Hz and H₁ = 2 mm. In addition, the system has two separate resonance zones in the torsion direction under 0 < f < 38 Hz and H₁ = 2 mm, and the system has one separate resonance zone in the vertical vibration under 150 Hz < f < 180 Hz and H₁ = 2 mm. This is because the structure of the rolling mill in horizontal, vertical, and torsion directions is different, the system has single direction resonance. Furthermore, the coupling relationships of the rolling mill make the system have multidirectional coupling resonance. The violent resonance can be reduced by reducing the external excitation amplitude of the rolling mill.

(a)

(b)

(c)

Figure 5: Frequency domain vibration signal of the work roll: (a) vibration acceleration in the torsional direction; (b) vibration acceleration in the vertical direction; (c) vibration acceleration in the horizontal direction.

4.2. Analysis of Bifurcation Characteristics

Figure 6 shows the dynamic bifurcation diagram of the work roll vibration under different gaps between the work roll bearing seat and the mill stand. Figures 7 and 8 show the phase path and Poincare section under different gaps.

(a)

(b)

(c)

Figure 6: Bifurcation diagram under varying structure gap: (a) torsional direction; (b) vertical direction; (c) horizontal direction.

(a)

(b)

(c)

Figure 7: Bifurcation diagram under 0.3 mm structure gap: (a) torsional direction; (b) vertical direction; (c) horizontal direction.

(a)

(b)

(c)

Figure 8: Bifurcation diagram under 0.6 mm structure gap: (a) torsional direction; (b) vertical direction; (c) horizontal direction.

As Figures 6(a)–6(c) show, the dynamic bifurcation characteristics of horizontal, vertical, and torsional vibrations are the same, which also reflects the coupling. When e₁ < 0.4 mm, the system vibration is periodic motion and the corresponding phase path and Poincare section are single, which indicates that the system runs steadily. When e₁ > 0.4 mm, the system enters chaotic motion and the corresponding phase path and Poincare section are multiple, which indicates that the system runs unstable and affects the quality of strip. As Figure 6(b) shows, when e₁ > 0.7 mm, the horizontal vibration amplitude of the system is less than the gap. But in fact, the probability of this situation is very small, so we do not do research. Therefore, the unstable resonance can be avoided by reducing the gap between the work roll bearing seat and the mill stand.

5. Conclusion

(1)In the dynamic nonlinear rolling process model, elastic recovery, multiroll equilibrium, and parameter coupling were considered. Furthermore, the elastic recovery of the exit strip has more influence on rolling torque than rolling force. In the coupled structure model of the rolling mill, the structure gap, roll offset, and structure coupling were considered. Finally, the piecewise nonlinear horizontal-vertical-torsional coupled vibration model of the rolling mill was established as the entrance strip thickness changes, which reflects the continuity of rolling.(2)From the study of amplitude-frequency characteristics of the hot rolling mill under varying external excitation amplitude, it was found that, in the system, there exist single direction resonance and multidirectional coupling resonance. In addition, the resonant amplitude of the hot rolling mill can be reduced by reducing the external excitation amplitude.(3)From the study of the bifurcation characteristics of the coupling vibration of the hot rolling mill under varying structural gaps, it was found that period motion and chaos motion exist in horizontal, vertical, and torsional directions. Further, the coupling vibration of the hot rolling mill can be stable if the gap is maintained below the critical value.

To sum up, we can reduce the amplitude of external excitation and structure gap to suppress the vibration of the rolling mill and ensure the stable operation of the system.

Data Availability

The structural and conditions data of rolling mill used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the Key Program of Natural Science Research of Higher Education in Anhui Province of China (No. KJ2018ZD005).

References

V. Panjković, R. Gloss, J. Steward, S. Dilks, R. Steward, and G. Fraser, “Causes of chatter in a hot strip mill: observations, qualitative analyses and mathematical modelling,” Journal of Materials Processing Technology, vol. 212, no. 4, pp. 954–961, 2012.

Z. Gao, Y. Zang, and L.-Q. Zeng, “Review of modelling and theoretical studies on chatter in the rolling mills,” Journal of Mechanical Engineering, vol. 51, no. 16, pp. 87–112, 2015.

L. Hadji and M. Avcar, “Nonlocal free vibration analysis of porous FG nanobeams using hyperbolic shear deformation beam theory,” Advances in Nano Research, vol. 10, no. 3, pp. 281–293, 2021.

F.-S. Du, G.-G. Wang, X.-L. Zang, and X.-T. Li, “Friction model for strip rolling,” Journal of Iron and Steel Research International, vol. 17, no. 7, pp. 19–23, 2010.

Z. Y. Jiang and A. K. Tieu, “Elastic-plastic finite element method simulation of thin strip with tension in cold rolling,” Journal of Materials Processing Technology, vol. 130-131, no. 20, pp. 511–515, 2002.

W.-G. Li, C. Liu, N. Feng, X. Chen, and X.-H. Liu, “Friction estimation and roll force prediction during hot strip rolling,” Journal of Iron and Steel Research International, vol. 23, no. 12, pp. 1268–1276, 2016.

S. Chen, W. Li, and X. Liu, “Calculation of rolling pressure distribution and force based on improved karman equation for hot strip mill,” International Journal of Mechanical Sciences, vol. 89, no. 12, pp. 256–263, 2014.

A. Heidari, M. R. Forouzan, and M. R. Niroomand, “Development and evaluation of friction models for chatter simulation in cold strip rolling,” International Journal of Advanced Manufacturing Technology, vol. 96, no. 2, pp. 2055–2075, 2018.

J.-Z. Cao, D.-W. Zhao, S.-H. Zhang, W. Peng, S.-Z. Chen, and D.-H. Zhang, “Analysis of hot tandem rolling force with logarithmic velocity field and EA yield criterion,” Journal of Iron and Steel Research International, vol. 21, no. 3, pp. 295–299, 2014.

G.-S. Ma, Y.-M. Liu, W. Peng et al., “A new model for thermo-mechanical coupled analysis of hot rolling,” Journal of the Brazilian Society of Mechanical Sciences and Engineering, vol. 39, no. 2, pp. 1–8, 2015.

J.-S. Son, D.-M. Lee, I.-S. Kim, and S.-G. Choi, “A study on on-line learning neural network for prediction for rolling force in hot-rolling mill,” Journal of Materials Processing Technology, vol. 164-165, pp. 1612–1617, 2005.

H. Li, B. Wen, and J. Zhang, “Asymptotic method and numerical analysis for self-excited vibration in rolling mill with clearance,” Shock and Vibration, vol. 8, Article ID 262743, 2001.

M. Mosayebi, F. Zarrinkolah, and K. Farmanesh, “Calculation of stiffness parameters and vibration analysis of a cold rolling mill stand,” International Journal of Advanced Manufacturing Technology, vol. 91, no. 9, pp. 4359–4369, 2017.

D. Han, P. Shi, and K. Xia, “Nonlinear torsional vibration dynamics behaviors of rolling mill’s multi-DOF main drive system under parametric excitation,” Journal of Applied Mathematics, vol. 2014, Article ID 202686, 7 pages, 2014.

X.-B. Fan, Y. Zang, and K. Jin, “Rolling process and its influence analysis on hot continuous rolling mill vibration,” Applied Physics, vol. 122, 1008 pages, 2016.

H. Alwan and M. Avcar, “Analytical solution of free vibration of FG beam utilizing different types of beam theories: a comparative study,” Computers and Concrete, vol. 26, no. 3, pp. 285–292, 2020.

M. Avcar, “Effects of material non-homogeneity and two parameter elastic foundation on fundamental frequency parameters of timoshenko beams,” Acta Physica Polonica A, vol. 130, no. 1, pp. 375–378, 2016.

R. Wang, Y. Peng, Y. Zhang, and J.-L. Sun, “Mechanism research of rolling mill coupled vibration,” Journal of Mechanical Engineering, vol. 49, no. 12, pp. 66–71, 2013.

R.-R. Peng, X.-Z. Zhang, and P.-M. Shi, “Coupled vibration behavior of hot rolling mill rolls under multinonlinear effects,” Shock and Vibration, vol. 2020, Article ID 6104028, 14 pages, 2020.

L.-Q. Zeng, Y. Zang, and Z.-Y. Gao, “Multiple-modal-coupling modeling and stability analysis of cold rolling mill vibration,” Shock and Vibration, vol. 2016, Article ID 2347386, 26 pages, 2016.

L. Zeng, Y. Zang, Z.-Y. Gao, and X.-C. Liu, “Study on overall coupled modeling of the rolling mill,” Journal of Mechanical Engineering, vol. 51, no. 14, pp. 46–53, 2015.

X. Wang and X. Yan, “Dynamic model of the hot strip rolling mill vibration resulting from entry thickness deviation and its dynamic characteristics,” Mathematical Problems in Engineering, vol. 2019, Article ID 5868740, 11 pages, 2019.

X. Lu, J. Sun, G. Li, Q. Wang, and D. Zhang, “Dynamic analysis of vibration stability in tandem cold rolling mill,” Journal of Materials Processing Technology, vol. 272, pp. 47–57, 2019.

X. Lu, J. Sun, G. Li, Z. Wang, and D. Zhang, “Stability analysis of a nonlinear coupled vibration model in a tandem cold rolling mill,” Shock and Vibration, vol. 2019, Article ID 4358631, 14 pages, 2019.

S.-Z. Chen, D.-H. Zhang, J. Sun, J.-S. Wang, and J. Sun, “Online calculation model of rolling force for cold rolling mill based on numerical integration,” in Proceedings of the 2012 24th Chinese Control And Decision Conference (CCDC), Taiyuan, China, May 2012.

Y. K. Sun, Model and Control of Cold and Hot Rolling Strip Mill, Metallurgical Industry Press, Beijing, China, 2010.

W. Guo, Vertical-horizontal Vibration Research on Simplified Model of Hot Continuous Rolling Mill during Dynamic Rolling Process, Anhui University of Technology, Ma’anshan, China, 2019.

Y. Z. Liu, Nonlinear Vibrations, Higher Education Press, Beijing, China, 2001.