# Coupled Vibration Behavior of Hot Rolling Mill Rolls under Multinonlinear Effects

Authors: Rongrong Peng; Xingzhong Zhang; Peiming Shi
Journal: Shock and Vibration
 (2020)

---

## Abstract

A coupled vibration model of hot rolling mill rolls under multiple nonlinear effects is established by considering the nonlinear spring force produced by the hydraulic cylinder, the nonlinear friction between the work rolls, the dynamic variation of rolling force, and the effect of external excitation as well as according to the structural constraints of a four-high hot rolling mill in the vertical and horizontal directions. The amplitude-frequency response equation of rolling mill rolls is determined by using a multiple-scale approximation method. Furthermore, use of actual data for simulation indicates that the internal resonance is the main cause of coupling vibration of the rolling mill rolls. In addition, changes in the movement displacement of the hydraulic cylinder and the coupling parameters strongly affect the coupling system of the rolling mill rolls. Finally, the study of the dynamic bifurcation characteristics of the rolling mill rolls indicates that, with varying external excitation amplitude, the vibration of rolls alternates between periodic motion, period-doubling motion, and chaotic motion in both vertical and horizontal directions. This is one of the reasons for the appearance of periodic light and dark stripes on the strip surface. Furthermore, the range of the external excitation amplitude (F0) at which the rolling mill roll system vibrates violently, that is, 5.68e5 N < F0 < 5.84e5 N and F0 > 6.12e5 N, must be avoided. The research results can provide a theoretical reference for further exploration of the coupling vibration mechanism of hot rolling mills.


---

## Body

## 1. Introduction

Hot rolling is an important process in the manufacturing and processing of steel strips. It is a high-temperature and large-reduction process, which imparts large fluidity in the workpiece. In addition, with the application of various new processing technologies, rolling mill vibration is an inevitable problem, which potentially loosens and wears out the mechanical parts, shortens the rolling mill life, and causes vibration of steel strips, thereby reducing the overall product quality and accuracy [1–3].Most studies on rolling mill vibration focus on a single vibration system, such as vertical, horizontal, or torsional vibration. For instance, Sun et al. considered the interaction between strip tension and thickness and analyzed the effect of multiple-source external disturbance on rolling mill vibration; they established a vertical vibration model of the rolling mill and further improved the control accuracy of tension and thickness of the hot rolling mill through optimization design [4]. Fan et al. analyzed and calculated the structural stability of hot rolling mill rolls and revealed that the gap between the framework and rolling mill rolls can cause the phenomenon of “vibration jump,” which destroys the balance of a rolling mill system [5]. Shi et al. studied the nonlinear torsional vibration characteristics of the main drive system of a rolling mill under different connection angles, damping coefficients, and friction force [6, 7]. Under continuous and high-speed rolling, different forms of coupled chatter are generated [8, 9]. This coupled vibration contributes to the low accuracy of a vibration model established on the basis of a single vibration system. Moreover, the research conclusions derived from such models cannot fully explain the complex nonlinear characteristics of rolling mill vibration. To address this problem, some research has been conducted on the chatter mechanism and behavior of coupled vibration [10–12]. Yun et al. established a 2-degree-of-freedom (DOF) linear coupling vibration model by considering the coupling effect of longitudinal and transverse movements of the roll and the effect of rolling speed on the rolling process [13, 14]. Yang et al. established a rolling mill model, including a rolling process model, a mill stand model, and a hydraulic servo system model, and they analyzed the effect of different working conditions on the coupling model [15]. Liu et al. considered the interaction between rolling mill structure vibration and workpiece vibration, proposed a model of workpiece-roll coupling vibration, studied the transition and bifurcation characteristics of vibration under the main resonance conditions, and introduced time-delay feedback control to suppress rolling mill vibration [16]. Wang et al. considered the friction force due to the roughness between rolls and the work roll movement on the basis of the rolling theory and lubrication theory. They established a coupling vibration model of the rolling interface under unsteady lubrication and the effect of multiple factors [17]. In addition, Lu et al. established a coupling vibration model of the dynamic rolling process and nonlinear vibration parameters and analyzed the bifurcation characteristics and stability of the vibration system using the Hopf bifurcation theorem [18]. In conclusion, although some achievements have been realized in the exploration of the coupling vibration mechanism of rolling mill rolls, there are still some limitations or problems with existing models and studies. Therefore, the complex vibration behavior of a rolling mill itself needs to be further explored and verified from different perspectives.In this study, we consider the effects of the nonlinear spring force of the hydraulic cylinder of a rolling mill, the nonlinear friction force between rolls, the nonlinear rolling force, and the structural constraints of the rolling mill and establish a coupling vibration model of a four-high hot rolling mill. We investigate the effects of tuning parameters, nonlinear stiffness, and coupling parameters on the amplitude-frequency characteristics of the vibration system. The simulation results show that a process of energy exchange occurs between the vertical and horizontal directions of the rolling mill vibration system. When the vibration parameters change, a jumping phenomenon occurs in both directions, causing instability of the coupling vibration system of the rolling mill. A change in the nonlinear spring force coefficient of the hydraulic cylinder considerably affects the vibration state of the rolling mill rolls, and changes in the coupling parameters lead to the complex nonlinear phenomenon of the vibration system. Thus, the coupling vibration model of a hot rolling mill established in this paper is confirmed to be effective. Finally, the bifurcation and chaos behavior of the coupled vibration system of the rolling mill under dynamic external excitation are studied. It is found that different periodic motions exist and the vibration alternates among different forms, which is one of the reasons for the appearance of periodic light and dark stripes on the strip surface. The results of this study can provide certain theoretical reference and technical support for reducing or restraining rolling mill vibration.


## 2. Vibration Modeling of Hot Rolling Mill Rolls

The mechanical structure of a hot rolling mill can be simply illustrated as shown in Figure1. In a rolling mill, the cylindrical block of the hydraulic cylinder is placed on the bearing bracket of the upper backup roll, and the piston is placed against the heel block under the mill stand. Therefore, the force is transmitted to the strip through the backup rolls and the upper and lower work rolls; then, the gap between the two work rolls is varied to enable rolling of strips to different thickness and specifications. Hence, the effect of nonlinear forces of the hydraulic system in the vertical direction cannot be ignored in the vibration model of rolling mill rolls.Figure 1
Structure of 1780 hot rolling mill.Furthermore, because the diameter and rotation speed of the backup rolls and work rolls differ, when the rolling speed changes, the friction force on the strip changes, and the rolling force of the roll system is not in the vertical direction. In other words, during tension rolling, the impact load of the hot rolling mill is extremely large when it bites into the strip, and the bearing bracket will impact the mill housing, which will degrade the stability of the rolling mill system, and the rolling force of the rolling mill rolls does not act in the vertical direction. Therefore, to improve the accuracy of a vibration model of the hot rolling mill, the nonlinear vibration factors of the rolling mill rolls in the horizontal direction must be considered. Consequently, by considering the effect of the nonlinear force of the rolling mill in the vertical and horizontal directions and combining it with rolling mill vibration, we can better explore the chatter mechanism and vibration behavior of a hot rolling mill.


### 2.1. Nonlinear Spring Force of the Hydraulic Cylinder

A hydraulic cylinder has many advantages in practical application and is widely used in the hydraulic screwdown device of a rolling mill. In this study, we consider a double-acting single-piston servo hydraulic cylinder in a hot rolling mill as an example for the analysis; its structure diagram is presented in Figure2.Figure 2
Structure diagram of a double-acting single-piston servo hydraulic cylinder.As Figure2 shows, A1 and A2 are the effective areas of piston rodless cavity and rod cavity of the hydraulic cylinder, respectively. VL1 and VL2 are the volumes of oil in the pipeline between the valve and the rodless cavity and between the valve and the rod cavity, respectively. L is the total stroke of the hydraulic cylinder, L0 is the initial position of the piston, y is the piston displacement, which is also the vibration displacement in the vertical direction of the rolling mill rolls, and k (y) is the total hydraulic spring stiffness of the hydraulic cylinder.During the working of the hydraulic cylinder, the piston can be regarded as a rigid body. The change in the piston displacement changes the pressure and oil volume in the two oil cavities, changing the oil stiffness accordingly. Therefore,k (y) is equivalent to the parallel stiffness of the two chamber’s hydraulic spring [19]:(1)ky=βeA12V1+βeA22V2=βeA12A1L0+y+VL1+βeA22A2L−L0−y+VL2,where βe is the elastic modulus of the oil volume and V1 and V2 are the volumes of the rodless cavity and rod cavity, respectively.BecauseVL1 and VL2 are very small compared with the volumes of the two cavities and may be ignored, the spring stiffness of the hydraulic cylinder can be simplified as(2)ky=βeA1L0+y+βeA2L−L0−y.Therefore, the spring stiffness of the hydraulic cylinder can be regarded as a function of vibration displacement, and the change law between them is shown in Figure3.Figure 3
Curve of the spring stiffness with respect to the vibration displacement of the hydraulic cylinder.By applying Taylor expansion to equation (2) at vibration displacement y = 0, we obtain(3)ky=βeA11L0−yL02+y2L03+βeA21L−L0−yL−L02+y2L−L03+Oy2.Due to the symmetry of spring potential energyU (y) and by neglecting the odd power term, we obtain(4)Uy=12αy2+14βy4,where α=βeA1/L0+βeA2/L−L0 and β=βeA1/L03+βeA2/L−L03.From the derivation of equation (4), the spring force of the hydraulic cylinder is(5)Fky=αy+βy3.Equation (5) indicates that the spring force of the hydraulic cylinder is a nonlinear function of vibration displacement.


### 2.2. Nonlinear Friction and Dynamic Rolling Force

The calculation accuracy of rolling force greatly affects the distribution of friction force and strip quality in the deformation zone. Many expressions for rolling force under different working conditions have been proposed by analyzing, simplifying, and regressing the characteristics of nonlinear coupling vibration of the rolling mill and dynamic components of rolling force based on actual tests [18, 20]. Of them, we adopt the Alexander–Ford formula for hot rolling [21]. The dynamic rolling process in the deformation zone is shown in Figure 4.Figure 4
Dynamic rolling process in the deformation zone.As Figure4 shows, τb and τf are the tension at the entrance and exit of the rolling, respectively. h0 and h1 are the thickness at the entrance and exit of the workpiece in steady state, respectively. h2 is the thickness at the exit of the workpiece in dynamic rolling, h2 = h1 + y, v0 is the rotation speed of the roll, and Fμ is the friction force on the work rolls:(6)P=2K−τb+τf2BQplc,where P is the rolling force, K is the strip deformation resistance, B is the average width of the workpiece, Qp is the influence coefficient of the stress state, and lc is the horizontal projection length of the roll contact arc in the deformation area:(7)lc=RΔh,Qp=1.08+1.79μδ⋅1−δRh2−1.02δ,K=1.15σ0expa1T+a2v010a3T+a4a6e0.4a5−a6−1e0.4,where Δh is the reduction, Δh = h0 − h1 − y, R is the roll radius, and μ is the friction coefficient between rolls. δ is the reduction rate, δ = Δh/h0, T is the deformation temperature, T = (t0 + 273)/1000, and e is the deformation degree, e = ln[1/(1 − δ)]. Regression coefficients σ0 and a1–a6 have a set of coefficients for different steel grades. In this study, the regression coefficient of the deformation resistance of ordinary carbon steel is taken [22] as follows: (8)a1=−2.878,a2=3.665,a3=0.1861,a4=−0.1216,a5=0.3795,a6=1.402,σ0=150.6.Considering that the friction coefficient between roll gaps varies with the fluctuation of rolling speed during rolling, Robert’s formula for friction coefficient is adopted [23]:(9)μ=Δh2R0.5+G1−0.5eG2v0e−G2x˙,where G1 and G2 are the friction characteristic coefficients, whose values are determined by specific friction model parameters; generally, G1 = 0.51 and G2 = 0.001. x˙ is the vibration speed in the horizontal direction. Because G2x˙≪1, equation (9) can be changed to μ=Δh/2R0.5+G1−0.5eG2v01−G2x˙+G2x˙2/2 and y ≤ h0 − h1; hence, equation (9) can be further expressed as(10)μ=μ0+Δμx˙,y,(11)Δμx˙,y=−b1x˙+b2x˙2−b3−b4x˙+b5x˙2y,μ0=g20.5+g1,b1=g1g2G2,b2=g1g2G222,b3=0.5+g1g3,b4=g1G2g3,b5=g1G222g3,g1=G1−0.5eG2v0,g2=h0−h12R,g3=2Rh0−h1,where μ0 is the friction coefficient of the roll gap in the steady state and Δμ (x˙, y) is the dynamic values of the friction coefficient.In steady-state rolling, the Taylor expansion of equation (6) at x˙ = 0 and y = 0 can be obtained as follows:(12)P=P0,0+ΔPx˙,y,(13)P0,0=d0=2K−τb+τf2BRh0−h1Qp′,ΔPx˙,y=d1x˙+d2y+d3x˙2+d4x˙y+d5y2,Qp′=1.08+1.79μ0Rh0−1.02h0−h1h0,d1=∂∂x˙P0,0,d2=∂∂yP0,0,d3=12!∂2∂x˙2P0,0,d4=∂2∂x∂yP0,0,d5=12!∂2∂2yP0,0.HereP0,0 is the rolling force in the steady state of the vibration system of the hot rolling mill; ΔPx,y is the dynamic value of the rolling force in the unsteady state. Therefore, equation (12) is the nonlinear rolling force, which is affected by the coupling of vertical vibration displacement and horizontal vibration speed, as shown in Figure 5.Figure 5
Curve of rolling force affected by the coupling of vertical vibration displacement and horizontal vibration velocity.


### 2.3. Nonlinear Coupled Vibration Model of Hot Rolling Mill Rolls

According to the structure diagram of the four-high hot rolling mill shown in Figure1, we consider the structural constraints on the roll and the effect of dynamic rolling force, the nonlinear spring force exerted by the hydraulic cylinder in the vertical direction, and the friction force change caused by the rolling speed fluctuation in the horizontal direction. Then, the mass concentration method is adopted, wherein the upper backup and work rolls are equivalent to a mass block, while the lower backup and work rolls are equivalent to a mass block. At the same time, in the experiment, the vibration displacement of the upper and lower rolls of the rolling mill is the same, but in the opposite directions. For the convenience of calculation and analysis, owing to the symmetry of rolling mill vibration, the rolling mill can be simplified as an effective and reliable 2-DOF vibration system [24]. The coupling vibration model of the upper rolls of the hot rolling mill under multiple nonlinear effects is established, as shown in Figure 6.Figure 6
Nonlinear coupled vibration model of hot rolling mill rolls.As Figure6 shows, m is the equivalent mass of the upper roll system of the rolling mill, v0 is the rotation speed of the work roll, F0 cosωt is the external excitation of the roll system, F0 is the external excitation amplitude, and ω is the external excitation frequency. k1 and c1 are the equivalent stiffness and damping between the upper roll system and the frame and the archway column, respectively, and k2 and c2 are the equivalent stiffness and damping between the upper roll system and the steel strip, respectively.The friction on the roll is small in the vertical direction and can be neglected. Therefore, from equations (10) and (12), the nonlinear friction on the roll in the horizontal direction can be obtained as follows:(14)Fμ=μP=μ0P0,0+ΔFμx˙,y,where ΔFμ (x˙, y) is the dynamic change of friction.According to the Lagrange theory, the coupled vibration dynamic equation of the hot rolling mill rolls can be obtained as follows:(15)mx¨+2c1x˙+2k1x=ΔFμx˙,y,my¨+c2y˙+k2y+γFky=ΔPx˙,y+F0cosωt,where x¨ and ÿ are the vibration accelerations in the horizontal and vertical directions, respectively. x˙ and y˙ are the vibration velocities in the horizontal and vertical directions, respectively. x and y are the vibration displacements in the horizontal and vertical directions, respectively, and γ is the constraint coefficient of the nonlinear spring force. Therefore, equation (15) can be simplified as(16)x¨+ω12x=−c1′x˙+κ1y+η1x˙2−ξ1x˙y−ς1y2+ϑ1x˙2y+ρ1x˙y2+ψ1x˙3,y¨+ω22y=−c2′y˙+κ2x˙+η2x˙2+ξ2x˙y+ς2y2−γβ′y3+F0∗cosωt,(17)ω12=2k1m,c1′=2c1−μ0d1+b1d1m,κ1=μ0d2−b1d2m,η1=μ0d3−b1d3m,ξ1=b3d1m,ς1=b3d2m,ϑ1=b4d1−b3d3+b2d2m,ρ1=b4d2m,ψ1=b2d1m,ω22=k2+γα−d2m,c2′=c2m,κ2=d1m,η2=d3m,ξ2=d4m,ς2=d5m,β′=γβm,F0∗=F0m.


## 2.1. Nonlinear Spring Force of the Hydraulic Cylinder

A hydraulic cylinder has many advantages in practical application and is widely used in the hydraulic screwdown device of a rolling mill. In this study, we consider a double-acting single-piston servo hydraulic cylinder in a hot rolling mill as an example for the analysis; its structure diagram is presented in Figure2.Figure 2
Structure diagram of a double-acting single-piston servo hydraulic cylinder.As Figure2 shows, A1 and A2 are the effective areas of piston rodless cavity and rod cavity of the hydraulic cylinder, respectively. VL1 and VL2 are the volumes of oil in the pipeline between the valve and the rodless cavity and between the valve and the rod cavity, respectively. L is the total stroke of the hydraulic cylinder, L0 is the initial position of the piston, y is the piston displacement, which is also the vibration displacement in the vertical direction of the rolling mill rolls, and k (y) is the total hydraulic spring stiffness of the hydraulic cylinder.During the working of the hydraulic cylinder, the piston can be regarded as a rigid body. The change in the piston displacement changes the pressure and oil volume in the two oil cavities, changing the oil stiffness accordingly. Therefore,k (y) is equivalent to the parallel stiffness of the two chamber’s hydraulic spring [19]:(1)ky=βeA12V1+βeA22V2=βeA12A1L0+y+VL1+βeA22A2L−L0−y+VL2,where βe is the elastic modulus of the oil volume and V1 and V2 are the volumes of the rodless cavity and rod cavity, respectively.BecauseVL1 and VL2 are very small compared with the volumes of the two cavities and may be ignored, the spring stiffness of the hydraulic cylinder can be simplified as(2)ky=βeA1L0+y+βeA2L−L0−y.Therefore, the spring stiffness of the hydraulic cylinder can be regarded as a function of vibration displacement, and the change law between them is shown in Figure3.Figure 3
Curve of the spring stiffness with respect to the vibration displacement of the hydraulic cylinder.By applying Taylor expansion to equation (2) at vibration displacement y = 0, we obtain(3)ky=βeA11L0−yL02+y2L03+βeA21L−L0−yL−L02+y2L−L03+Oy2.Due to the symmetry of spring potential energyU (y) and by neglecting the odd power term, we obtain(4)Uy=12αy2+14βy4,where α=βeA1/L0+βeA2/L−L0 and β=βeA1/L03+βeA2/L−L03.From the derivation of equation (4), the spring force of the hydraulic cylinder is(5)Fky=αy+βy3.Equation (5) indicates that the spring force of the hydraulic cylinder is a nonlinear function of vibration displacement.


## 2.2. Nonlinear Friction and Dynamic Rolling Force

The calculation accuracy of rolling force greatly affects the distribution of friction force and strip quality in the deformation zone. Many expressions for rolling force under different working conditions have been proposed by analyzing, simplifying, and regressing the characteristics of nonlinear coupling vibration of the rolling mill and dynamic components of rolling force based on actual tests [18, 20]. Of them, we adopt the Alexander–Ford formula for hot rolling [21]. The dynamic rolling process in the deformation zone is shown in Figure 4.Figure 4
Dynamic rolling process in the deformation zone.As Figure4 shows, τb and τf are the tension at the entrance and exit of the rolling, respectively. h0 and h1 are the thickness at the entrance and exit of the workpiece in steady state, respectively. h2 is the thickness at the exit of the workpiece in dynamic rolling, h2 = h1 + y, v0 is the rotation speed of the roll, and Fμ is the friction force on the work rolls:(6)P=2K−τb+τf2BQplc,where P is the rolling force, K is the strip deformation resistance, B is the average width of the workpiece, Qp is the influence coefficient of the stress state, and lc is the horizontal projection length of the roll contact arc in the deformation area:(7)lc=RΔh,Qp=1.08+1.79μδ⋅1−δRh2−1.02δ,K=1.15σ0expa1T+a2v010a3T+a4a6e0.4a5−a6−1e0.4,where Δh is the reduction, Δh = h0 − h1 − y, R is the roll radius, and μ is the friction coefficient between rolls. δ is the reduction rate, δ = Δh/h0, T is the deformation temperature, T = (t0 + 273)/1000, and e is the deformation degree, e = ln[1/(1 − δ)]. Regression coefficients σ0 and a1–a6 have a set of coefficients for different steel grades. In this study, the regression coefficient of the deformation resistance of ordinary carbon steel is taken [22] as follows: (8)a1=−2.878,a2=3.665,a3=0.1861,a4=−0.1216,a5=0.3795,a6=1.402,σ0=150.6.Considering that the friction coefficient between roll gaps varies with the fluctuation of rolling speed during rolling, Robert’s formula for friction coefficient is adopted [23]:(9)μ=Δh2R0.5+G1−0.5eG2v0e−G2x˙,where G1 and G2 are the friction characteristic coefficients, whose values are determined by specific friction model parameters; generally, G1 = 0.51 and G2 = 0.001. x˙ is the vibration speed in the horizontal direction. Because G2x˙≪1, equation (9) can be changed to μ=Δh/2R0.5+G1−0.5eG2v01−G2x˙+G2x˙2/2 and y ≤ h0 − h1; hence, equation (9) can be further expressed as(10)μ=μ0+Δμx˙,y,(11)Δμx˙,y=−b1x˙+b2x˙2−b3−b4x˙+b5x˙2y,μ0=g20.5+g1,b1=g1g2G2,b2=g1g2G222,b3=0.5+g1g3,b4=g1G2g3,b5=g1G222g3,g1=G1−0.5eG2v0,g2=h0−h12R,g3=2Rh0−h1,where μ0 is the friction coefficient of the roll gap in the steady state and Δμ (x˙, y) is the dynamic values of the friction coefficient.In steady-state rolling, the Taylor expansion of equation (6) at x˙ = 0 and y = 0 can be obtained as follows:(12)P=P0,0+ΔPx˙,y,(13)P0,0=d0=2K−τb+τf2BRh0−h1Qp′,ΔPx˙,y=d1x˙+d2y+d3x˙2+d4x˙y+d5y2,Qp′=1.08+1.79μ0Rh0−1.02h0−h1h0,d1=∂∂x˙P0,0,d2=∂∂yP0,0,d3=12!∂2∂x˙2P0,0,d4=∂2∂x∂yP0,0,d5=12!∂2∂2yP0,0.HereP0,0 is the rolling force in the steady state of the vibration system of the hot rolling mill; ΔPx,y is the dynamic value of the rolling force in the unsteady state. Therefore, equation (12) is the nonlinear rolling force, which is affected by the coupling of vertical vibration displacement and horizontal vibration speed, as shown in Figure 5.Figure 5
Curve of rolling force affected by the coupling of vertical vibration displacement and horizontal vibration velocity.


## 2.3. Nonlinear Coupled Vibration Model of Hot Rolling Mill Rolls

According to the structure diagram of the four-high hot rolling mill shown in Figure1, we consider the structural constraints on the roll and the effect of dynamic rolling force, the nonlinear spring force exerted by the hydraulic cylinder in the vertical direction, and the friction force change caused by the rolling speed fluctuation in the horizontal direction. Then, the mass concentration method is adopted, wherein the upper backup and work rolls are equivalent to a mass block, while the lower backup and work rolls are equivalent to a mass block. At the same time, in the experiment, the vibration displacement of the upper and lower rolls of the rolling mill is the same, but in the opposite directions. For the convenience of calculation and analysis, owing to the symmetry of rolling mill vibration, the rolling mill can be simplified as an effective and reliable 2-DOF vibration system [24]. The coupling vibration model of the upper rolls of the hot rolling mill under multiple nonlinear effects is established, as shown in Figure 6.Figure 6
Nonlinear coupled vibration model of hot rolling mill rolls.As Figure6 shows, m is the equivalent mass of the upper roll system of the rolling mill, v0 is the rotation speed of the work roll, F0 cosωt is the external excitation of the roll system, F0 is the external excitation amplitude, and ω is the external excitation frequency. k1 and c1 are the equivalent stiffness and damping between the upper roll system and the frame and the archway column, respectively, and k2 and c2 are the equivalent stiffness and damping between the upper roll system and the steel strip, respectively.The friction on the roll is small in the vertical direction and can be neglected. Therefore, from equations (10) and (12), the nonlinear friction on the roll in the horizontal direction can be obtained as follows:(14)Fμ=μP=μ0P0,0+ΔFμx˙,y,where ΔFμ (x˙, y) is the dynamic change of friction.According to the Lagrange theory, the coupled vibration dynamic equation of the hot rolling mill rolls can be obtained as follows:(15)mx¨+2c1x˙+2k1x=ΔFμx˙,y,my¨+c2y˙+k2y+γFky=ΔPx˙,y+F0cosωt,where x¨ and ÿ are the vibration accelerations in the horizontal and vertical directions, respectively. x˙ and y˙ are the vibration velocities in the horizontal and vertical directions, respectively. x and y are the vibration displacements in the horizontal and vertical directions, respectively, and γ is the constraint coefficient of the nonlinear spring force. Therefore, equation (15) can be simplified as(16)x¨+ω12x=−c1′x˙+κ1y+η1x˙2−ξ1x˙y−ς1y2+ϑ1x˙2y+ρ1x˙y2+ψ1x˙3,y¨+ω22y=−c2′y˙+κ2x˙+η2x˙2+ξ2x˙y+ς2y2−γβ′y3+F0∗cosωt,(17)ω12=2k1m,c1′=2c1−μ0d1+b1d1m,κ1=μ0d2−b1d2m,η1=μ0d3−b1d3m,ξ1=b3d1m,ς1=b3d2m,ϑ1=b4d1−b3d3+b2d2m,ρ1=b4d2m,ψ1=b2d1m,ω22=k2+γα−d2m,c2′=c2m,κ2=d1m,η2=d3m,ξ2=d4m,ς2=d5m,β′=γβm,F0∗=F0m.


## 3. Solution of the Coupling Internal Resonance Response of Hot Rolling Mill Rolls

The multiscale method has strong problem-solving ability and good computing ability, and considering that it is widely used in solving coupled vibration system problems, we use the multiscale method to solve the dynamic response of the system [25]. The nonlinear term in equation (16) is assigned a small parameter ε, and the fast and slow time scales T0 = t and T1 = εt, respectively, are introduced to obtain(18)ddt=D0+εD1+⋯,d2dt2=D02+2εD0D1+⋯,where Dn (n = 0, 1) is the partial differential sign, and Dn = ∂/∂Tn. The solution for equation (16) is set as follows:(19)x=x0T0,T1+εx1T0,T1,y=y0T0,T1+εy1T0,T1.Substituting equations (18) and (19) into equation (16), and separating the terms by the order of ε, we obtain(20)D02x0+ω12x0=0,D02y0+ω12y0=0,(21)D02x1+ω12x1=−2D0D1x0−c1′D0x0+κ1y0+η1D0x02−ξ1y0D0x0−ς1y02+ϑ1y0D0x02+ρ1y02D0x0+ψ1D0x03,D02y1+ω22y1=−2D0D1y0−c2′D0y0+κ2D0x0+η2D0x02+ξ2y0D0x0+ς2y02−γβ′y03+F0∗cosωT0.The solution of equation (20) is set as follows:(22)x0=AT1eiω1T0+cc,y0=BT1eiω2T0+cc,where A and B are the undetermined complex functions and cc is the complex conjugate of the left side items. By substituting equation (22) in equation (21), we obtain(23)D02x1+ω12x1=−2iω1D1A−iω1c1′A+2iρ1ω1ABB¯+3iψ1ω13A2A¯eiω1T0+κ1B+2ϑ1AA¯Beiω2T0−iω1ξ1ABeiω1+ω2T0+iω1ξ1A¯Beiω1−ω2T0−ω12ϑ1A2Beiω2+2ω1T0−ω12ϑ1A¯2Beiω2−2ω1T0+iω1ρ1AB2ei2ω2+ω1T0−iω1ρ1A¯B2ei2ω2−ω1T0−iω13ψ1A3e3iω1T0−η1ω12A2e2iω1T0−ς1B2e2iω2T0+2ω12η1AA¯−2ς1BB¯+cc,D02y1+ω22y1=−2iω2D1B−iω2c2′B−3γβ2B2B¯eiω2T0+iω1κ2Aeiω1T0+iω1ξ2ABeiω1+ω2T0−iω1ξ2A¯Beiω2−ω1T0−η2ω12A2e2iω1T0+ς2B2e2iω2T0−γβ′B3e3iω2T0+2ω12η2AA¯+2ς2BB¯+F0∗eiωT02+cc.Considering the internal resonance of the system,ω = ω2 + εδ and ω1 = ω2 + εδ1 are set, where δ and δ1 are the tuning parameters; δ can better describe the change between ω and ω2 and represent the value range between them. To avoid the duration term, A and B must satisfy the following:(24)−2iω1D1A−iω1c1′A+2iρ1ω1ABB¯+3iψ1ω13A2A¯+κ1B+2ϑ1AA¯Be−iδ1T1=0,−2iω2D1B−iω2c2′B−3γβ′B2B¯+iω1κ2Aeiδ1T1+F0∗eiδT12=0.We setA=1/2aT1iϕ1T1 and B=1/2bT1eiϕ2T1, where a and b are the amplitudes of rolling mill rolls in the horizontal and vertical directions, respectively. Substituting A and B into equation (24) and separating the real and imaging parts, the average equation of the coupling system can be obtained as follows:(25)a˙=−12c1′a+14ρ1ab2+38ψ1ω12a3+1ω112κ1+14ϑ1a2bsinθ1,aθ˙1+θ˙2=δ−δ1a−1ω112κ1+14ϑ1a2bcosθ1,b˙=−12c2′b+ω1ω2κ2acosθ1+F0∗2sinθ2,bθ˙2=δb−38ω2γβ′b3−ω1ω2κ2asinθ1−F0∗2cosθ2,where θ1 = φ2 − φ1 − δ1T1 and θ2 = δT1 − φ2. Considering the occurrence of a periodic movement in the coupling system of the hot rolling mill rolls, a˙=b˙=θ˙1=θ˙2=0 holds. Then, substituting this value into equation (25) and eliminating φ1 and φ2, the amplitude-frequency response equation of the coupling vibration system can be written as(26)−12c1′+14ρ1b2+38ψ1ω12a22a2+δ−δ12a2=1ω1212κ1+14ϑ1a22b2,−12c2′b+Mδ−δ12+δb−38ω2γβ′b3−M12c1′−14ρ1b2−38ψ1ω12a22=F0∗22,where M=ω12κ2a2/ω21/2κ1+1/4ϑ1a2b.


## 4. Amplitude Frequency Characteristics of Coupled Vibration of Hot Rolling Mill Rolls

Table1 shows the relevant vibration parameters of the 1780 hot rolling mill to conduct simulation analysis.Table 1
Structural and technological parameters of the 1780 hot rolling mill.
ParametersValueβe (GPa)1.6A1 (m2)0.6361A2 (m2)0.3243L (m)0.110L0 (m)0.064B (m)1.5R (m)0.42V0 (m·s−1)2.5τb (MPa)5.5τf (MPa)3.8H0 (m)0.0141H1 (m)0.0082t0 (°C)996m (Kg)1.44 × 105c1 (N·s·m−1)5.20 × 103k1 (N·m−1)7.31 × 109c2 (N·s·m−1)8.85 × 105k2 (N·m−1)2.08 × 1010F0 (MN)0.55Ε0.01Figure7 shows the amplitude-frequency diagram of the coupling vibration rolls of the hot rolling mill as a function of tuning parameter δ. Set initial value γβ = 1.05e12 (N/m−3) and ρ1 = −3.3. An energy exchange process is observed between the vertical and the horizontal direction, which is a special phenomenon that occurs in response to internal resonance. When tuning parameter δ > −10 Hz, amplitude a in the horizontal direction basically remains constant, while amplitude b in the vertical direction increases rapidly and then decreases after reaching the maximum. When 16.5 Hz < δ < 21.5 Hz, a increases and b decreases gradually. When −8.40 Hz < δ < −6.87 Hz, there are multiple solutions in the two directions and a jump phenomenon exists, causing instability of the coupling vibration system of the rolling mill in this range. Due to external disturbance, ω1≈ω2 and ω≈ω2, which imply that the external excitation frequency is close to the natural frequency in both directions and that resonance occurs. In addition, Figure 7 indicates a large linear frequency in the horizontal direction. This is because of the energy exchange to the horizontal direction when the vibration in the vertical direction is weakened. Moreover, the external excitation frequency is similar to the natural frequency in the horizontal direction; therefore, the amplitude in the horizontal direction increases rapidly. The resonance can be avoided by changing the external excitation amplitude or the linear stiffness between the rolls to maintain each frequency separated from each other.Figure 7
Amplitude-frequency curve of coupling vibration of the hot rolling mill roll system.Figure8 shows the amplitude-frequency curve of the coupling vibration of the hot rolling mill rolls under varying nonlinear spring force coefficient γβ of the hydraulic cylinder for tuning parameter δ = 14.5 Hz. In the figure, stable branches are plotted with solid lines, while unstable branches are plotted with dotted lines. When γβ = 0, that is, when the hydraulic cylinder piston is at the initial position, the vertical vibration of the rolling mill is only affected by the linear stiffness between the stands, and the amplitude of the vibration system is small. At this point, the rolling mill is in a steady state. With the increase in the piston movement displacement of the hydraulic cylinder, the corresponding nonlinear stiffness increases. For γβ ≥ 2.605e12 (N/m−3), there are multiple solutions, the amplitude increases, and the resonance region expands. This indicates that the change in the piston position of the hydraulic cylinder greatly affects the vibration state of the rolling mill rolls. In actual rolling, hydraulic oil with good temperature characteristics and large bulk modulus of elasticity should be used as much as possible. Moreover, seals with good performance should be selected to prevent an increase in the level of impurities due to the long-term use of oil. However, the change in nonlinear stiffness is caused by the change in viscosity, leading to abnormal vibration.Figure 8
Amplitude-frequency curve of the rolling mill rolls as a function of parameterγβ.Figure9 shows the amplitude-frequency curve of the coupling vibration of the hot rolling mill rolls under varying coupling term parameter ρ1, which is composed of horizontal vibration speed x˙ and vertical vibration displacement y. In the plot, the stable branches are indicated with solid lines, and the unstable branches are indicated with dotted lines. For tuning parameter δ = 14.5 Hz, when −12 < ρ1 < −3, there are multiple solutions and the vibration amplitude increases gradually, causing severe vibration of rolling mill rolls. The vibration form of the rolling mill in this range of coupling parameter is complex and variable. As ρ1 increases, the vibration state of the rolling mill gradually stabilizes. This further indicates that the study of the coupling vibration of the hot rolling mill rolls is meaningful and would provide some guiding significance to explain the vibration phenomenon of rolling mills.Figure 9
Amplitude-frequency curve of the rolling mill rolls as a function of coupling term parameterρ1.


## 5. Bifurcation Characteristics of Coupling Vibration of Hot Rolling Mill Rolls

In this study, using the bifurcation and chaos theory, we determine the critical value and parameter range from the bifurcation solution of the coupling vibration system of the hot rolling mill rolls, so as to restrain and avoid rolling mill vibration.


### 5.1. Bifurcation Characteristics of Vibration in the Vertical Direction

Figure10 shows the dynamic bifurcation diagram of the vertical coupling vibration of the hot rolling mill rolls. Figures 11 and 12 show the phase path and the Poincare section under different external excitation amplitude F0.Figure 10
Bifurcation diagram of the vertical vibration under varying external excitation amplitude.Figure 11
Phase path of the vertical vibration under varying external excitation amplitude: (a)F0 = 5.54e5 N, (b) F0 = 5.60e5 N, (c) F0 = 5.90e5 N, (d) F0 = 6.05e5 N, and (e) F0 = 6.55e5 N.
(a)(b)(c)(d)(e)Figure 12
Poincare section of vertical vibration under varying external excitation amplitude: (a)F0 = 5.54e5 N, (b) F0 = 5.60e5 N, (c) F0 = 5.90e N, (d) F0 = 6.05e5 N, and (e) F0 = 6.55e5 N.
(a)(b)(c)(d)(e)Figure10 indicates that when the external excitation amplitude F0 < 5.54e5 N, the vibration is period-1 motion, the corresponding phase path is a closed curve (Figure 11(a)), and the corresponding Poincare section is a fixed point (Figure 12(a)). This implies that the vertical vibration of the rolling mill has a stable solution and the vibration displacement is small. With increasing F0, the steady-state vibration becomes period-2 motion, the corresponding phase path is a closed curve formed after two circles (Figure 11(b)), and the Poincare section is two fixed points (Figure 12(b)). In this status, the vibration displacement of the rolling mill in the vertical direction increases gradually, resulting in an increase in the roll gap and the increase in the rolling workpiece reduction, causing uneven oscillation marks on the strip surface. With further increase in F0, the vibration displacement will also increase. As shown in Figure 10, the system will enter into paroxysmal chaos, resulting in up and down vibration disorder of the rolling mill rolls, seriously affecting strip quality and aggravating oscillation marks. Therefore, the range 5.68e5 N < F0 < 5.84e5 N should be avoided in an actual rolling process. Further, the vibration system degenerates into period-3 motion; the corresponding phase path is shown in Figure 11(c) and the Poincare section is shown in Figure 12(c). As F0 continues to increase, there is a period-doubling bifurcation, that is, the vibration becomes period-6 motion (Figures 11(d) and 12(d)), and when F0 > 6.12e5 N, it will transform to chaos motion from period-doubling motion; the phase path is an open circular curve (Figure 11(e)), while the Poincare section is a set of many points (Figure 12(e)). This is because the vibration frequency of the external excitation gradually becomes close to or equal to the derived frequency of the rolling mill rolls, thus achieving resonance, which causes the amplitude of vertical vibration to attain the maximum value. Timely measures to restrain this resonance must be taken to prevent the occurrence of breakage of the steel strip.


### 5.2. Bifurcation Characteristics of Vibration in the Horizontal Direction

Figure13 shows the dynamic bifurcation diagram of the horizontal vibration of the hot rolling mill rolls. Its specific vibration form is periodic-1 motion ⟶ periodic-2 motion ⟶ chaos motion ⟶ degeneration to period-2 motion ⟶ paroxysmal chaos ⟶ chaos motion. The phase path and Poincare sections corresponding to different external excitation amplitudes F0 are shown in Figures 14 and 15.Figure 13
Bifurcation diagram of horizontal vibration of the hot rolling mill rolls under varying external excitation amplitude.Figure 14
Phase paths of horizontal vibration of the rolling mill rolls under varying external excitation amplitude: (a)F0 = 5.31e5 N, (b) F0 = 5.95e5 N, and (c) F0 = 6.20e5 N.
(a)(b)(c)Figure 15
Poincare section of horizontal vibration of the rolling mill rolls under varying external excitation amplitude, (a)F0 = 5.31e5 N, (b) F0 = 5.95e5 N, and (c) F0 = 6.20e5 N.
(a)(b)(c)WhenF0 < 5.31e5 N, the horizontal vibration of the hot rolling mill is stable; the corresponding phase path is shown in Figure 14(a), and the Poincare section is shown in Figure 15(a). As F0 gradually increases, the horizontal vibration system directly enters chaotic motion after period-2 motion, and when F0 > 5.95e5 N, the system gradually degenerates into period-2 motion. The corresponding Poincare sections of the phase path are, respectively, shown in Figures 14(b) and 15(b). Then, the system passes through paroxysmal chaos when F0 > 6.20e5 N, and it finally enters chaotic motion. The corresponding phase trajectory is shown in Figure 14(c). The trajectory is a set of unclosed curves, indicating that the system does not have periodicity at this time; the corresponding Poincare section is shown in Figure 15(c). Note that the system undergoes a periodic motion, period-doubling motion, and chaotic motion in the horizontal direction, and the movement alternates among these forms. Thus, different forms of chatter marks are caused on the steel strip.Figures10 and 13 indicate that changes in the external excitation amplitude will cause changes in the vibration displacement in both vertical and horizontal directions. Therefore, fluctuation of the speed of the rolls is caused, which in turn leads to changes in nonlinear friction force and dynamic rolling force, resulting in abnormal vibration of the coupling system of the rolling mill rolls. Comparison of Figures 10 and 13 indicates that the two ranges of F0, namely, 5.68e5 N < F0 < 5.84e5 N and F0 > 6.12e5 N, can cause severe vibration of rolling mill rolls in the vertical and horizontal directions.


## 5.1. Bifurcation Characteristics of Vibration in the Vertical Direction

Figure10 shows the dynamic bifurcation diagram of the vertical coupling vibration of the hot rolling mill rolls. Figures 11 and 12 show the phase path and the Poincare section under different external excitation amplitude F0.Figure 10
Bifurcation diagram of the vertical vibration under varying external excitation amplitude.Figure 11
Phase path of the vertical vibration under varying external excitation amplitude: (a)F0 = 5.54e5 N, (b) F0 = 5.60e5 N, (c) F0 = 5.90e5 N, (d) F0 = 6.05e5 N, and (e) F0 = 6.55e5 N.
(a)(b)(c)(d)(e)Figure 12
Poincare section of vertical vibration under varying external excitation amplitude: (a)F0 = 5.54e5 N, (b) F0 = 5.60e5 N, (c) F0 = 5.90e N, (d) F0 = 6.05e5 N, and (e) F0 = 6.55e5 N.
(a)(b)(c)(d)(e)Figure10 indicates that when the external excitation amplitude F0 < 5.54e5 N, the vibration is period-1 motion, the corresponding phase path is a closed curve (Figure 11(a)), and the corresponding Poincare section is a fixed point (Figure 12(a)). This implies that the vertical vibration of the rolling mill has a stable solution and the vibration displacement is small. With increasing F0, the steady-state vibration becomes period-2 motion, the corresponding phase path is a closed curve formed after two circles (Figure 11(b)), and the Poincare section is two fixed points (Figure 12(b)). In this status, the vibration displacement of the rolling mill in the vertical direction increases gradually, resulting in an increase in the roll gap and the increase in the rolling workpiece reduction, causing uneven oscillation marks on the strip surface. With further increase in F0, the vibration displacement will also increase. As shown in Figure 10, the system will enter into paroxysmal chaos, resulting in up and down vibration disorder of the rolling mill rolls, seriously affecting strip quality and aggravating oscillation marks. Therefore, the range 5.68e5 N < F0 < 5.84e5 N should be avoided in an actual rolling process. Further, the vibration system degenerates into period-3 motion; the corresponding phase path is shown in Figure 11(c) and the Poincare section is shown in Figure 12(c). As F0 continues to increase, there is a period-doubling bifurcation, that is, the vibration becomes period-6 motion (Figures 11(d) and 12(d)), and when F0 > 6.12e5 N, it will transform to chaos motion from period-doubling motion; the phase path is an open circular curve (Figure 11(e)), while the Poincare section is a set of many points (Figure 12(e)). This is because the vibration frequency of the external excitation gradually becomes close to or equal to the derived frequency of the rolling mill rolls, thus achieving resonance, which causes the amplitude of vertical vibration to attain the maximum value. Timely measures to restrain this resonance must be taken to prevent the occurrence of breakage of the steel strip.


## 5.2. Bifurcation Characteristics of Vibration in the Horizontal Direction

Figure13 shows the dynamic bifurcation diagram of the horizontal vibration of the hot rolling mill rolls. Its specific vibration form is periodic-1 motion ⟶ periodic-2 motion ⟶ chaos motion ⟶ degeneration to period-2 motion ⟶ paroxysmal chaos ⟶ chaos motion. The phase path and Poincare sections corresponding to different external excitation amplitudes F0 are shown in Figures 14 and 15.Figure 13
Bifurcation diagram of horizontal vibration of the hot rolling mill rolls under varying external excitation amplitude.Figure 14
Phase paths of horizontal vibration of the rolling mill rolls under varying external excitation amplitude: (a)F0 = 5.31e5 N, (b) F0 = 5.95e5 N, and (c) F0 = 6.20e5 N.
(a)(b)(c)Figure 15
Poincare section of horizontal vibration of the rolling mill rolls under varying external excitation amplitude, (a)F0 = 5.31e5 N, (b) F0 = 5.95e5 N, and (c) F0 = 6.20e5 N.
(a)(b)(c)WhenF0 < 5.31e5 N, the horizontal vibration of the hot rolling mill is stable; the corresponding phase path is shown in Figure 14(a), and the Poincare section is shown in Figure 15(a). As F0 gradually increases, the horizontal vibration system directly enters chaotic motion after period-2 motion, and when F0 > 5.95e5 N, the system gradually degenerates into period-2 motion. The corresponding Poincare sections of the phase path are, respectively, shown in Figures 14(b) and 15(b). Then, the system passes through paroxysmal chaos when F0 > 6.20e5 N, and it finally enters chaotic motion. The corresponding phase trajectory is shown in Figure 14(c). The trajectory is a set of unclosed curves, indicating that the system does not have periodicity at this time; the corresponding Poincare section is shown in Figure 15(c). Note that the system undergoes a periodic motion, period-doubling motion, and chaotic motion in the horizontal direction, and the movement alternates among these forms. Thus, different forms of chatter marks are caused on the steel strip.Figures10 and 13 indicate that changes in the external excitation amplitude will cause changes in the vibration displacement in both vertical and horizontal directions. Therefore, fluctuation of the speed of the rolls is caused, which in turn leads to changes in nonlinear friction force and dynamic rolling force, resulting in abnormal vibration of the coupling system of the rolling mill rolls. Comparison of Figures 10 and 13 indicates that the two ranges of F0, namely, 5.68e5 N < F0 < 5.84e5 N and F0 > 6.12e5 N, can cause severe vibration of rolling mill rolls in the vertical and horizontal directions.


## 6. Conclusions

In this study, a coupled vibration model of hot rolling mill rolls under the effect of multiple nonlinear forces was established. Below is the summary of the analysis and steps involved in establishing the model.(1)
By analyzing the actual structure and working principle of a double-acting single-piston servo hydraulic cylinder, the nonlinear spring force produced by it was obtained. Considering the velocity fluctuation at the entry of the strip workpiece in the horizontal direction and the changes in the vibration displacement in the vertical direction, the nonlinear friction force between the rolls was obtained. Finally, by considering the dynamic variation of the rolling force, the coupling vibration model of a four-high hot rolling mill under the effect of multiple nonlinearity was established.(2)
Based on the amplitude-frequency response equation, and by using the actual rolling parameters of the 1780 four-high hot rolling mill, the main reason for the severe resonance of the rolling mill was found to be the occurrence of the internal resonance when the external excitation frequency is close to the derived frequency in the vertical and horizontal directions. This results in the instability of the system and occurrence of the jump phenomenon. Further, changes in the movement displacement of the hydraulic cylinder and the coupling term parameters considerably contribute to the changes in the amplitude and resonance range of the coupling vibration system of the hot rolling mill rolls.(3)
From the study of the bifurcation characteristics of the coupled vibration system of the hot rolling mill under varying external excitation amplitude, it was found that period and period-doubling motions exist in both vertical and horizontal directions, and the vibration alternates between different forms. Therefore, periodic light and dark stripes appear on the strip. The results indicate that the abnormal vibration of rolling mill rolls can be mitigated if the external excitation amplitude is maintained below the critical value.



---

*Source: 6104028-2020-06-10.xml*