# Stability Analysis of a Nonlinear Coupled Vibration Model in a Tandem Cold Rolling Mill

Authors: Xing Lu; Jie Sun; Guangtao Li; Zhenhua Wang; Dianhua Zhang
Journal: Shock and Vibration
 (2019)

---

## Abstract

Mill chatter in tandem cold rolling mill is a major rejection to the quality and production of the strips. In most mill vibration models, either the roll mass is usually limited to vibrate in vertical direction and vertical-horizontal directions, or the multiple rolls system is simplified to a single mass system. However, the torsional chatter is also a typical type of mill chatter, and the presence of intermediate roll and backup roll will affect the overall vibration of the mill structure system. In this paper, a newly cold rolling mill vibration model coupled with the dynamic rolling processing model and nonlinear vibration model is proposed with the consideration of dynamic coupling and nonlinear characteristics of the rolling process, multiroll equilibrium, and roll movement in both vertical-horizontal-torsional directions. By using Hopf bifurcation theorem and Routh–Hurwitz determinant, the existence of the Hopf bifurcation point of the mill vibration system and bifurcation characteristics are analyzed. At last, the influence of different rolling conditions on the stability of the coupled mill system is investigated, and these results can also be used to design an optimum rolling schedule and determine the appearance of mill chatter under certain rolling conditions.


---

## Body

## 1. Introduction

With the optimization and upgrading in the steel industry, the urgent demands for high strip quality and strip yield have become more stringent. However, due to the abnormal mill vibration, the quality and production of the strip have been restricted. Abnormal vibration in the tandem cold rolling mill is a common and complicated phenomenon. Researchers tend to consider that mill chatter is a typical self-excited vibration [1–3], which consists of third-octave mode chatter (150–250 Hz), fifth-octave mode chatter (500–700 Hz), and torsional chatter.In the mill chatter analysis, the most typical models are finite element method (FEM) models, multiple degrees of freedom structure model, and mill drive system torsional structure models. Niroomand et al. [4–6] studied the third-octave mode chatter and fifth-octave mode chatter in the vertical direction by using FEM analysis models. In these papers, the backup roll was set as a mass point, and a connector was used to represent the displacement of the backup roll and mill stand housing. In FEM models, the structure mode coupling characteristics cannot be fully reflected. For most vertical vibration models, the mill stand was modeled with multiple degrees of freedom. Kapil et al. [7] used the further simplified structure model to analyze the frequency response in detail. Kimura et al. [8, 9] and Heidari et al. [10–13] illustrated the influence of lubrication characteristic on mill chatter. However, in these studies, only the vertical vibration phenomena were investigated. But, in the actual rolling process, the roll vibrated in more than one direction. So, it is necessary to establish a coupled chatter model which allows roll to vibrate in two directions. Yun et al. presented a dynamic coupled chatter model that allows roll vibrate in two directions and investigated the relationship between mill vibration and the rolling conditions [2]. Kim et al. [14] proposed a vertical and horizontal coupling chatter model and analyzed the dynamic response of roll displacement under different rolling pressures. In the research of torsional vibration on the rotating shaft, Krot [15] analyzed the nonlinear torsional vibration caused by angular and radial wear in gear. Liu et al. [16, 17] established a multiple degree of freedom of mill main drive system and studied the torsional chatter characteristics in the main drive system of the mill stand.Bifurcation was an important theory of nonlinear field, which implied the sudden change of the topological property of the dynamic system when the parameters change in the dynamic characteristic system. Hopf bifurcation was a typical type of bifurcation system. It was named after the name of the German mathematician Eberhard Hopf who proposed the concept. It referred to the bifurcation of the periodic oscillation and the appearance of a stable limit cycle, which is caused by the dynamic change of the parameter [18]. Hopf bifurcation had theoretical value not only in the dynamic bifurcation research but also closely related to the generation of self-excited vibration in engineering. Lee found that the Hopf bifurcation phenomenon in the wing flutter system was the main cause of the casing vibration [19], and the magnitude of the pitch angle played an important role in this vibration problem. The characteristics and localization of multiple pendulum-type centrifugal pendulum vibration absorbers were investigated to suppress torsional vibrations of shaft by using Hopf bifurcation theory [20]. Zhang and Dai used nonlinear creep theory to study bifurcation phenomenon and stability of China Railways High-speed vehicle (CRH-3) [21].Many researchers have found that the dynamic rolling mill system was often accompanied with the generation of limit cycle when unstable oscillation occurs during rolling. So, Hopf bifurcation theorem was also introduced into the mill stand vibration analysis. Gao et al. [22, 23] established a dynamic model of rolling mill system and analyzed the Hopf bifurcation and mill stability. Liu et al. [24] established a dynamic vibration considering the torsional vibration of gear and the horizontal vibration of the roll, and the stability of the mill drive system was analyzed. Peng et al. [25] proposed a swing chatter model coupling the vertical-horizontal-torsional vibration and studied the influence of meshing impact force on the roll system stability.Researchers have shown the existence of negative damping effects in mill structure. But the relationship between dynamic rolling conditions and damping coefficient is not clearly illustrated. Besides, the Hopf bifurcation phenomenon in the rolling mill drive vibration system is an important cause of system instability. The occurrence of the Hopf bifurcation limit cycle determines the stability and stable range of the mill. Then, vibration problem caused by negative damping effect can be transformed into the analysis of the Hopf bifurcation phenomenon in the vibration model. Scholars have done plenty scientific works about vibration in cold rolling. However, most of those theses were either concentrated on vertical vibration or torsional vibration, or the whole roll system was simplified by one single roll. Therefore, a further multiroll and nonlinear vertical-horizontal-torsional coupled vibration model of the universal crown control mill (UCM mill) needs to be proposed and analyzed. And based on this proposed rolling mill vibration model, the bifurcation phenomenon and stability of the mill system is analyzed to illustrate the dynamic characteristics of the UCM mill more comprehensively.


## 2. Dynamic Rolling Process Model

The geometry of roll gap in the cold rolling is presented in Figure1. O is the original center of the work roll, and O′ is the new center of the flattened roll. xin, xn, and xout are the entry position, neutral position, and exit position. α and γ are the bite angle and neutral angle. vin and vout are the strip speed at the entry and exit point. hin, hn, and hout are the strip thickness at the entry, neutral, and exit point. xw is the dynamic vibration displacement of work roll in the horizontal direction. yw is the dynamic vibration displacement of work roll in the vertical direction. Tf and Tb are the front and back tension.Figure 1
Vibration geometry model of the roll gap.In cold rolling, Hitchcock indicated that the flattened work roll radiusR′ can be evaluated as(1)R′=R1+161−υ2πE1Pw⋅hin−hout,where υ is Poisson’s ratio; in this research υ=0.3. w is the width of the strip. E1 is Young’s modulus of the work roll. P is the total rolling force. R and R′ are the work roll original radius and roll flattened radius.Considering work roll elastic deformation, the contact arc length in the roll gap is(2)l′=Rhin−hout+81−υ2πE1Rp¯2+81−υ2πE1Rp¯,where p¯ is the average rolling force.Assuming the contour of the roll gap can be written as a parabolic curve, the strip thickness at any point in the roll gap can be determined as [11](3)hx=hc+x2R′,where hc is the roll gap spacing distance along the centerline of the work roll and x is the position of the strip along the rolling direction in the roll gap.From equation (3), the entry position of the strip can be given as(4)xin=xw−R′hin−hc.So, the dynamic roll gap spacing of the roll centerline can be related as [11, 22](5)hout=hc+2yw.Based on equation (5) and considering mill elastic deformation, the variation of strip thickness at exit is expressed as(6)dhout=2yw+dPKw,where dP is the rolling force variation and kw is the equivalent spring constant of the work roll.And strip thickness at the neutral point can be calculated by [8](7)hn=hin+hout2−12μloghinhout×1−Tf/khout1−Tb/khin.The neutral angle is calculated by(8)γ=houtR′×tanhoutR′×hn2,where khin and khout are the deformation resistance at strip entry and exit position and μ is the friction coefficient in the roll gap.Based on equation (7), strip thickness fluctuation at the neutral point can be expressed as(9)dhn=∂hn∂hindhin+∂hn∂houtdhout+∂hn∂TfdTf+∂hn∂TbdTb=12−12μhindhin+12+12μhoutdhout+12μkhout−TfdTf−12μkhin−TbdTb.Considering the mass conservation law to the mass flow continuity equation,(10)Vrhn=vinhin=vouthout,where Vr is the work roll circumference velocity.The strip speed variation at the entry and exit can be expressed as [8](11)dvin=Vrhindhn,dvout=Vrhout2houtdhn−hndhout.According to Hooke’s law, the variation of tension is proportional to integral of speed difference between entry ofi stand and exit of i−1 stand [26]. Assuming the exit speed of i−1 stand and entry speed of i+1 stand are always kept constant, the tension variation at the entry and exit can be expressed as(12)dTf,i=−E25.5025∫dvout,idt,dTb,i=E25.5025∫dvin,idt,where subscript i means i stand. And in this research, the strip length between i and i+1 stand is 5.5025 m.In the most rolling mill vibration models, the friction coefficient was usually set as a constant value. However, in the cold rolling, the lubrication characteristics in the roll gap were significantly different under the high rolling speed and low speed. Sims and Arthur [27] found that when the rolling speed was higher than 0.25 m/s, the friction coefficient showed a nonlinear characteristic which is caused by unsteady rolling conditions. As a result, it is necessary to investigate the dynamic vibration stability in the UCM mill with the consideration of nonlinear friction characteristics. Considering the roll wear and other conditions into consideration, the friction coefficient in this research is obtained by [28](13)μ=μ0+μV⋅e−Vr/V0⋅1+cR⋅Ra−Ra0⋅1+cW1+L/L0,where μ0 and μV are the basic friction coefficient that related to lubrication and lubrication coefficient related to rolling speed, V0 is the basic rolling speed, Ra and Ra0 are actual work roll roughness and roughness of unused work roll, CR is the work roll roughness coefficient, cW is the coefficient related to cumulative rolling strip length, and L and L0 are the cumulative rolling strip length of work roll and datum rolling length.The rolling force is calculated by Bland-Ford-Hill equation [29, 30]:(14)P=wl′k¯Qpnt,where(15)Qp=1.08+1.79εμR′hin−1.02ε,nt=1−1−μtTf+μtTbk¯,μt=0.7,ε=hin−houthin,and k¯ is the deformation resistance of the rolling strip. In this research, deformation resistance of the strip is expressed by the following equation:(16)k¯=A+B⋅εΣ1−C⋅e−D⋅εΣ,where εΣ is average strain, which can be given by εΣ=2/3⋅lnH0/h¯, h¯=hin+2hout/3. A, B, C, and D are references in deformation resistance model. In this research, A = 500.00 MPa, B = 179.00 MPa, C = 0.20, and D = 5.00. H0 is the thickness of the raw material.The investigation of UCM mill showed that No. 4 and No. 5 stand chatter frequently, so a comprehensive mill vibration monitoring test was implemented in a 1450 mm UCM plant by using CTC-type accelerometer sensors. The accelerometer sensors were mounted on the work roll bearings of No. 5 stand. When rolling the most common products of 0.35 mm strip in this plant, the vibration signals during the stable rolling stage are shown in Figure2. It illustrated that although the horizontal rolling force was relatively small compared with the vertical rolling force, the horizontal rolling force was still big enough. Therefore, it was sufficiently significant to establish a comprehensive model which allows the rolls to vibrate in both horizontal and vertical directions.Figure 2
Time domain vibration signal of the work roll: (a) vibration acceleration in the vertical direction; (b) vibration acceleration in the horizontal direction.
(a)
(b)Considering the dynamic equilibrium of strip inx direction, the horizontal rolling force can be calculated as follows:(17)Px=∫0aPsinθdθ+∫γaPμcosθdθ−∫0γPμcosθdθ=w2Tfhout−Tbhin.The vertical rolling force can be written as(18)Py=∫0αPcosθdθ+∫γαPμsinθdθ−∫0γPμsinθdθ=wl′k¯Qpnt∫0αcosθdθ+∫γαμsinθdθ−∫0γμsinθdθ=wl′k¯Qpnt∫0αcosθdθ−μcosα+2μcosγ−μ.Considering that the bite angleα and neutral angle γ are extremely small, it assumed that the vertical rolling force is equal to the total rolling force [14, 22], so equation (18) can be expressed as(19)Py≈wl′k¯Qpnt∫0αcosθdθ≈wl′k¯Qpnt=P.Assuming the thickness fluctuation and tension variation are caused by roll gap vibration as separated disturbances, the deviation of the rolling force in thex direction can be written as(20)dPx=∂Px∂hindhin+∂Px∂houtdhout+∂Px∂TbdTb+∂Px∂TfdTf.Based on equation (17), the deviation of the rolling force in the x direction can be also expressed as(21)dPx=−wTb2dhin+wTfdyw−whin2dTb+whout2dTf.The variation of the rolling force in they direction can be defined as(22)dPy=∂P∂hindhin+∂P∂houtdhout+∂P∂TfdTf+∂P∂TbdTb,where(23)∂P∂hin=wl′k¯nt−1.792μR′hin3+32×1.79μhoutR′hin5−1.02houthin2,∂P∂hout=wl′k¯nt−1.79μhinR′hin+1.021hin,∂P∂Tb=−0.7wl′Qp,∂P∂Tf=−0.3wl′Qp.The rolling torque of the upper roll system can be deduced as [23, 25](24)M=w∫0γμp¯Rdθ−w∫γαμp¯Rdθ=wμp¯R2γ−α=PμR2γ−α.Considering the bite angleα and neutral angle γ are extremely small, α≈hin−hout/R′, γ≈sinγ=xw−xn/R′, and xn=R′hn−hc.Also, the variation of the rolling torque can be calculated by(25)dM=∂M∂hindhin+∂M∂houtdhout+∂M∂TfdTf+∂M∂TbdTb+∂M∂xwdxw,where(26)∂M∂hin=wl′k¯nt2γ−α−1.792μR′hin3+32×1.79μhoutR′hin5−1.02houthin2−12wl′k¯QpntμR1hin−houtR′,∂M∂hout=wl′k¯nt2γ−α−1.79μhinR′hin+1.021hin+12wl′k¯QpntμR1hin−houtR′,∂M∂Tb=−0.7wl′QpμR2γ−α,∂M∂Tf=−0.3wl′QpμR2γ−α,∂M∂xw=wl′k¯QpntμR21R′−α.The calculated rolling conditions in this research are shown in Table1. The grade of rolling strip is MRT5, and the exit strip width is 900 mm.Table 1
Rolling conditions.
Parameters
#1 stand
#2 stand
#3 stand
#4 stand
#5 stand
Entry thickness,hin (mm)
2.000
1.269
0.762
0.460
0.300
Exit thickness,hout (mm)
1.269
0.762
0.460
0.300
0.200
Rolling force,P (kN)
7838.863
8455.166
7781.567
7500.996
7722.711
Exit tension,Tf (MPa)
140.554
165.252
188.160
205.575
68.000
Entry tension,Tb (MPa)
55.000
140.554
165.252
188.160
205.575
Friction coefficient,μ
0.065
0.060
0.040
0.030
0.020
Roll velocity,Vr (m s−1)
3.381
5.652
9.425
14.441
22.121
Strip entry speed,vin (m s−1)
2.250
3.546
5.903
9.774
14.977Figure3 shows the comparison of calculated rolling force in this rolling process model and actual rolling force in the experiment. The rolling conditions are shown in Table 1. As shown in Figure 3, the maximum relative error between actual rolling force and calculated rolling force is 7.23%. The comparative results of the simulation and experimental data show that this chatter model has a high accuracy.Figure 3
Comparison between the actual rolling force and calculated rolling force.


## 3. Dynamic Nonlinear Vibration Model of the UCM Mill

In most vibration research studies, the rolling mill stand was simplified as a mass-spring system. To investigate the displacements of each roll and interaction between rolls, a model with multiple degrees of freedom is necessary. However, with the consideration of the symmetry of the mill structure along the rolling pass line, the multiple degrees of freedom chatter model can be simplified. Figure4 is the simplified nonlinear vertical-horizontal-torsional coupled vibration model considering the roll movement in more than one direction.Figure 4
Geometry of the vertical-horizontal-torsional vibration model: (a) simplified vertical-horizontal-torsional coupled vibration model; (b) geometry of coupled vibration model inx-y direction.
(a)
(b)As shown in Figure4, the top work roll (WR) would vibrate in both vertical and horizontal directions. While the intermediate roll (IMR) and backup roll (BUR) were not allowed to vibrate in horizontal direction due to the assembly and restriction of the roll bearings [25].A horizontal roll vibration displacement of WR is generated due to the installation offset ofe between WR and IMR, and the equivalent spring constant in horizontal direction can be calculated by(27)Kwx=eRw+RimKw,where Rw and Rim are the roll radius of work roll and intermediate roll and Kw is the contact spring constant between upper work roll and intermediate roll.With the consideration of the coupling movement of the rolls in both vertical, horizontal, and torsional directions, the dynamic equilibrium equation of the upper roll system was given by(28)Mwx¨w+Cwxx˙w+Kwxxw=ΔPx,Mwy¨w+Cwyy˙w−y˙im+Kwyyw−yim=ΔPy,Mimy¨im−Cwyy˙w−y˙im+Cimy˙im−y˙b−Kwyyw−yim+Kimyim−yb=0,Mby¨b−Cimy˙im−y˙b+Cby˙b−Kimyim−yb+Kbyb=0,Jwθ¨w+Cwjθ˙w+Kwjθw=ΔM,where xw is the WR vibration displacement in horizontal direction, yw, yim, and yb are the vibration displacement of WR, IMR, and BUR in vertical direction, θw is the WR torsional vibration angular displacement, ΔPx and ΔPy are the fluctuation of rolling force in horizontal and vertical directions, ΔM is the variation of rolling torque, and Mw, Mim, and Mb are the mass of work roll, intermediate roll, and backup roll, respectively.The mill stand spring constant is the total stiffness of the whole rolling mill consisting of all the rolls. The mill spring constant (K) is measured by a stiffness calibration test. The spring constants of work roll (Kw), intermediate roll (Kim), and backup roll (Kb) are calculated by an FEM model of the UCM mill [5, 6, 12]. The roll vibration displacement under different rolling forces was obtained by a simulation experiment. The spring constants of the work roll, intermediate roll, and backup roll were obtained by analyzing the dynamic roll vibration displacement under different rolling forces.Damping coefficients of the intermediate and backup rolls are calibrated by a simulation program [31, 32]. First, the dynamic differential equation of the system is regularized to obtain the equivalent damping matrix. Then, the damping regular matrix is transformed into the original coordinate system, and let the damping ratio be 0.03. Finally, the equivalent damping value of the roll system is calculated. Work roll rotational inertia mass, work roll rotational spring constant, and work roll rotational damping coefficient are obtained from the work of Zeng et al. [23], as the design and mill structure of rolling stands are similar to one another. The roll diameters and work roll installation offset are obtained from the UCM rolling mill installation drawings. Table 2 shows the structure specifications of the UCM mill in equation (28).Table 2
Structure specifications of the UCM mill.
Parameters
Value
Work roll mass,Mw (Kg)
3914
Intermediate roll mass,Mim (Kg)
5410
Backup roll mass,Mb (Kg)
31330
Work roll contact spring constant,Kw (N m−1)
0.91 × 1010
Intermediate roll spring constant,Kim (N m−1)
1.03 × 1010
Backup roll spring constant,Kb (N m−1)
5.85 × 1010
Work roll damping coefficientx direction, Cwx (N·s m−1)
5.00 × 104
Work roll damping coefficient iny direction, Cwy (N·s m−1)
1.64 × 105
Intermediate roll damping coefficient,Cim (N·s m−1)
5.70 × 105
Backup roll damping coefficient,Cb (N·s m−1)
5.10 × 106
Work roll rotational inertia mass,Jw (Kg·m2)
1381
Work roll rotational spring constant,Kwj (N·m·rad−1)
7.90 × 106
Work roll rotational damping coefficient,Cwj (N·m·s rad−1)
4178
Work roll diameter,Rw (mm)
212.5
Intermediate roll diameter,Rim (mm)
245
Backup roll diameter,Rb (mm)
650
Work roll installation offset,e (mm)
5Equation (28) is the nonlinear dynamic equation of the vertical-horizontal-torsional coupled vibration model. Convert the above differential equation into the following state matrix form:(29)X˙=FX,μ+Hdhin,where(30)X=xw,x˙w,yw,y˙w,yim,y˙im,yb,y˙b,θw,θ˙w,dTf,dTbT=x1,x2,x3,…,x12T,H=0,−wTb2Mw,0,∂P∂hinMw,0,0,0,0,0,∂M∂hinJw,0,0T,F=010000000000F2,1F2,2F2,300000000000010000000000F4,3F4,4F4,5F4,600000000000100000000F6,3F6,4F6,5F6,6F6,7F6,800000000000100000000F8,5F8,6F8,7F8,800000000000001000F10,2F10,300000F10,9F10,10F10,11F10,1200F11,30000000F11,11F11,1200F12,30000000F12,11F12,12,F2,1=−KwxMw,F2,2=−CwxMw,F2,3=wTfMw,F4,3=2∂P/∂hout−KwyMw,F4,4=−CwyMw,F4,5=KwyMw,F4,6=CwyMw,F6,3=KwyMim,F6,4=CwyMim,F6,5=−Kwy+KimMim,F6,6=−Cwy+CimMim,F6,7=KimMim,F6,8=CimMim,F8,5=KimMb,F8,6=CimMb,F8,7=−Kim+KbMb,F8,8=−Cim+CbMb,F10,2=∂MJw∂xw,F10,3=2∂MJw∂hout,F10,9=−KwjJw,F10,10=−CwjJw,F10,11=∂MJw∂Tf,F10,12=∂MJw∂Tb,F11,3=−2E2Vr∂hn/∂houthout−hn5.5025hout2,F11,11=−E2Vr∂hn/∂Tf5.5025hout,F11,12=−E2Vr∂hn/∂Tb5.5025hout,F12,3=2E2Vr∂hn/∂hout5.5025hin,F12,11=E2Vr∂hn/∂Tf5.5025hin,F12,12=E2Vr∂hn/∂Tf5.5025hin.Then, equation (29) is the nonlinear dynamic state equation of the vertical-horizontal-torsional coupled vibration system considering the nonlinear characteristic during the rolling process, the movement of rolls in both vertical-horizontal directions, and the torsional vibration.


## 4. Analysis of Vibration Stability of the Roll System

### 4.1. Hopf Bifurcation Parameter Calculation

The friction characteristic would cause the instability of the dynamic state equation, so in this researchμ has been chosen as the Hopf bifurcation parameter. Assuming the entry thickness of each stand is constant, the equilibrium point is X¯, where FX¯=0. The system equilibrium point can be transformed to the origin by coordinate transformation, X0=0,0,0,0,0,0,0,0,0,0,0,0T. Let JFX be the Jacobian matrix of FX. So, the Hopf bifurcation can be obtained by judging all the eigenvalues of JFX. However, it is difficult to calculate the eigenvalues of JFX since the calculation amount is quite large. Then, the Routh–Hurwitz determinant is used to calculate the Hopf bifurcation point [33].The characteristic equation of Jacobian matrixJFX of the dynamic vertical-horizontal-torsional coupling vibration model is assumed to be(31)a12λ12+a11λ11+⋯+a2λ2+a1λ+a0=0.The main determinant(32)D1=a1,D2=a1a0a3a2,⋯Di=a1a0⋯0a3a2⋯0⋯⋯⋯⋯a2i−1a2i−2⋯ai,where ai=0 if i<0 or i>n, and a12=1.Assumingμ∗ is the Hopf bifurcation, the Hopf bifurcation conditions are [23, 24, 34](33)a0μ∗>0,D1μ∗>0,⋯,Dn−2μ∗>0,Dn−1μ∗=0,dΔi−1μ∗dμ≠0.According to these assumptions, the Hopf bifurcation point ofμ∗ occurs at μ∗=0.0048 and μ∗=0.0401.


### 4.2. Analysis of Stable Range of Friction

When the value ofμ∗ is smaller than the first bifurcation point of μ∗=0.0048, such as μ∗=0.0040, the response of work roll displacement and Hopf bifurcation diagram are shown in Figure 5. The roll displacement will diverge rapidly, and the equilibrium point of the equilibrium equation is unstable. Besides, as shown in Figure 5(b), the phase trajectory will disperse outward of the origin in a ring shape and the radius curvature of the phase trajectory increases gradually. It is concluded that the roll vibration system is unstable and mill chatter appears.Figure 5
Dynamic response of the work roll displacement and the phase diagram underμ∗=0.0040: (a) time history of work roll displacement; (b) phase diagram.
(a)
(b)As shown in Figure6(a), when μ∗=0.0048, the work roll displacement exhibits as an equal amplitude vibration, and the roll vibrates at its balanced position. This steady region refers to the maximum and minimum fluctuation range of work roll. Meanwhile, Figure 6(b) shows that the phase diagram will begin to converge toward a stable limit cycle. The strip will produce a slight thickness fluctuation within the quality requirements, and the mill system is stable under these rolling conditions.Figure 6
Dynamic response of the work roll displacement and the phase diagram underμ∗=0.0048: (a) time history of work roll displacement; (b) phase diagram.
(a)
(b)Whenμ∗=0.0300, the roll system’s time history response and the phase diagram are shown in Figure 7. It can be seen from these figures that the phase trajectory will converges toward a smaller and smaller circle near the origin point. The equilibrium point of work roll is asymptotically stable.Figure 7
Dynamic response of the work roll displacement and the phase diagram underμ∗=0.0300: (a) time history of work roll displacement; (b) phase diagram.
(a)
(b)The roll system stability about the critical parameter ofμ∗=0.0401 is shown in Figure 8. When at the critical Hopf bifurcation parameter, displacement of work roll firstly becomes larger and finally converges to an equal amplitude oscillation. The equilibrium point is unstable, and the Hopf bifurcation diagram presents a ring shape. However, dynamic range of the fixed limit cycle in Figure 8 is a little larger compared with Figure 6. This is mainly because the friction is bigger, and as a result, the fluctuation of the response rolling force ΔPi under the friction coefficient of 0.0401 is much bigger.Figure 8
Dynamic response of the work roll displacement and the phase diagram underμ∗=0.0401: (a) time history of work roll displacement; (b) phase diagram.
(a)
(b)The numerical solution of work roll displacement under the Hopf bifurcation point ofμ∗=0.0420 is shown in Figure 9. The roll displacement will disperse outward all the time. Besides, the phase diagram of the work roll will disperse outward of the cycle, and the radius curvature of the phase trajectory increases gradually. It illustrates that the mill vibration system is unstable and mill chatter is produced.Figure 9
Dynamic response of the work roll displacement and the phase diagram underμ∗=0.0420: (a) time history of work roll displacement; (b) phase diagram.
(a)
(b)Whenμ∗<0.0048, the roll system is unstable, and when μ∗=0.0048, the system is stable. So the Hopf bifurcation of μ∗=0.0048 is the subcritical bifurcation point, and μ∗=0.0401 is the supercritical bifurcation point.Many researchers attribute the change of stability in rolling mill system to the sudden change of the rolling force fluctuation in dynamic equilibrium equation [11, 25, 35]. The optimal setting range of friction coefficient is 0.0048–0.0401, and mill chatter appears when the friction is either too high or too low. This is mainly due to the fluctuation of rolling force. When the friction is high, with the increase of friction coefficient, the neutral angle becomes bigger, and then the neutral point shifts to the entry side, so the distance between the entry and neutral point becomes shorter, and the variation of the entry speed decreases [11]. However, the exit forward slip and exit speed fluctuation will increase. Based on equation (12), the tension variation at the entry position (dTb) and exit position (dTf) is smaller compared with low friction coefficient rolling condition, so the fluctuation of the response rolling force ΔPi increases, and the stability of the stand decreases. Meanwhile, when the friction is quite small, ∂P/∂hout will increase with the decrease of friction coefficient; F4,3 in equation (29) decreases, so the instantaneous dynamic structure spring constant decreases, and the fluctuation of the response rolling force ΔPi caused by the strip exit thickness variation will increase, so the stability of the stand is decreased.


### 4.3. Stable Range Analysis of Other Rolling Parameters

#### 4.3.1. Calculation of Limited Rolling Speed

Based on this coupled vertical-horizontal-torsional vibration analysis model, rolling speed has been chosen as the Hopf bifurcation parameter to obtain the limited stability rolling speed of the UCM tandem cold rolling mill, and other rolling conditions are shown in Table1. The actual bifurcation point of rolling speed is 25.1 m/s. Figure 10 shows the time history of work roll displacement under the rolling speed of 25.2 m/s. According to Figure 10, the displacement amplitude of the work roll shows a divergence trend, and the self-excited vibration happens. The critical rolling speed calculated by the proposed chatter model is 25.1 m/s.Figure 10
Work roll time history displacement under the rolling speed of 25.2 m/s.Figure11 shows the tension and roll gap variation in No. 4 and No. 5 stands obtained from the Production Data Acquisition (PDA) system. As is shown in Figure 11(b), when the rolling speed is 22.1 m/s, the back and front tension of the No. 5 stand fluctuate within a small range and the mill stand is stable. However, when the rolling speed is raised to 23.0 m/s, the back and front tension of the No. 5 stand fluctuate more violently, and the front tension of the No. 5 stand begins to diverge as shown in Figure 11(a). It is found that the roll gap of No. 4 and No. 5 stands fluctuates more violently by comparing Figures 11(c) and 11(d). Meanwhile, a sudden and harsh vibration noise will be heard in the No. 5 stand, and finally the strip breaks. It shows mill chatter appears under the speed of 23.0 m/s, and the actual critical rolling speed is 23.0 m/s.Figure 11
Tension and roll gap time history: (a) back and front tension of No. 5 stand under the speed of 23.0 m/s; (b) back and front tension of No. 5 stand under the speed of 22.1 m/s; (c) roll gap of No. 4 and No. 5 stands under the speed of 23.0 m/s; (d) roll gap of No. 4 and No. 5 stands under the speed of 22.1 m/s.
(a)
(b)
(c)
(d)The calculated critical rolling speed is 25.1 m/s, while the actual critical rolling speed is 23.0 m/s. There is a certain relative error between the calculated critical speed and actual critical rolling speed. This is due to the small coefficient difference in Jacobian matrix. With the change of rolling speed, the actual value of the tension and friction is changed. And the calculated value of tension factornt and friction coefficient in ∂P/∂h and ∂M/∂h is a little different with the actual value. So, the calculated limited rolling speed is a little bigger than the actual value. However, by comparing the critical rolling speed calculated by the proposed model with the actual critical rolling speed, it can be seen that the proposed model is still accurate.The vibration is largely influenced by the friction coefficient and rolling speed. Based on equation (13), Figure 12 shows the friction coefficient response of No. 5 stand versus speed. Despite the changes of speed, the friction coefficient is still in the stability range of 0.0048–0.0401. It is divided into two regions. In region 1, both the friction coefficient and rolling speed are in the stable range. In region 2, the friction coefficient is still in the stable range, while the speed exceeds the limited speed range. When the rolling speed is lower than 25.1 m/s, the mill system is stable without changing the lubrication conditions. And when the rolling speed is higher than 25.1 m/s, the stand is unstable. However, the friction coefficient is in the middle value of stable range under the speed is 25.1 m/s, so mill chatter might be eliminated by reducing the basic friction coefficient (by improving the lubrication characteristics, such as increasing the concentration of emulsion, increasing spray flow, and improving the property of the base oil).Figure 12
Friction coefficient response of No. 5 stand versus speed.


#### 4.3.2. Calculation of Maximum Reduction

Assuming the entry thickness is constant, the exit thickness is variable. And the exit thickness has been chosen as the Hopf bifurcation parameter. Other calculated rolling conditions are shown in Table1. Then, the calculated bifurcation point of exit thickness is 0.168 mm. The optimal setting range of reduction for No. 5 stand is 0–44%. When the reduction of No. 5 stand is 44%, the increase of reduction leads to a lower mill stand stability. The reason is also related to the increase of rolling force variation under the larger reduction rolling conditions.As can be seen in Table1, in normal rolling schedule, the reduction rate is 33.3%. And the max reduction is up to 44%. In reality, considering there are other factors that influence the stand stability, such as rolling speed, and lubrication conditions, the max reduction is lower than 44%. In general, the effect of the reduction in the mill stability is lower than the improvement of friction coefficient.


#### 4.3.3. Front Tension Stable Range

When the work roll velocity is 22.1 m/s, the calculated bifurcation point of exit tension is 1.5–133 MPa. Small front tension means that the rolling force will be quite large. So, under the same excitation interferences, the variation of rolling force will be much bigger, and the stand stability is decreased. When the front tension is big enough, tension factornt decreases, leading to a larger value in ∂P/∂h and ∂M/∂h. As a result, the fluctuation of rolling force and roll torque increases, and the stability of the system is reduced as well. So, the influence of tension on the mill stability is smaller than other parameters.These results can be used to estimate mill chatter in a certain stand. And the stable range of friction, reduction, and tension can be used to adjust the rolling schedule in the UCM mill to avoid mill chatter under high rolling speed.


## 4.1. Hopf Bifurcation Parameter Calculation

The friction characteristic would cause the instability of the dynamic state equation, so in this researchμ has been chosen as the Hopf bifurcation parameter. Assuming the entry thickness of each stand is constant, the equilibrium point is X¯, where FX¯=0. The system equilibrium point can be transformed to the origin by coordinate transformation, X0=0,0,0,0,0,0,0,0,0,0,0,0T. Let JFX be the Jacobian matrix of FX. So, the Hopf bifurcation can be obtained by judging all the eigenvalues of JFX. However, it is difficult to calculate the eigenvalues of JFX since the calculation amount is quite large. Then, the Routh–Hurwitz determinant is used to calculate the Hopf bifurcation point [33].The characteristic equation of Jacobian matrixJFX of the dynamic vertical-horizontal-torsional coupling vibration model is assumed to be(31)a12λ12+a11λ11+⋯+a2λ2+a1λ+a0=0.The main determinant(32)D1=a1,D2=a1a0a3a2,⋯Di=a1a0⋯0a3a2⋯0⋯⋯⋯⋯a2i−1a2i−2⋯ai,where ai=0 if i<0 or i>n, and a12=1.Assumingμ∗ is the Hopf bifurcation, the Hopf bifurcation conditions are [23, 24, 34](33)a0μ∗>0,D1μ∗>0,⋯,Dn−2μ∗>0,Dn−1μ∗=0,dΔi−1μ∗dμ≠0.According to these assumptions, the Hopf bifurcation point ofμ∗ occurs at μ∗=0.0048 and μ∗=0.0401.


## 4.2. Analysis of Stable Range of Friction

When the value ofμ∗ is smaller than the first bifurcation point of μ∗=0.0048, such as μ∗=0.0040, the response of work roll displacement and Hopf bifurcation diagram are shown in Figure 5. The roll displacement will diverge rapidly, and the equilibrium point of the equilibrium equation is unstable. Besides, as shown in Figure 5(b), the phase trajectory will disperse outward of the origin in a ring shape and the radius curvature of the phase trajectory increases gradually. It is concluded that the roll vibration system is unstable and mill chatter appears.Figure 5
Dynamic response of the work roll displacement and the phase diagram underμ∗=0.0040: (a) time history of work roll displacement; (b) phase diagram.
(a)
(b)As shown in Figure6(a), when μ∗=0.0048, the work roll displacement exhibits as an equal amplitude vibration, and the roll vibrates at its balanced position. This steady region refers to the maximum and minimum fluctuation range of work roll. Meanwhile, Figure 6(b) shows that the phase diagram will begin to converge toward a stable limit cycle. The strip will produce a slight thickness fluctuation within the quality requirements, and the mill system is stable under these rolling conditions.Figure 6
Dynamic response of the work roll displacement and the phase diagram underμ∗=0.0048: (a) time history of work roll displacement; (b) phase diagram.
(a)
(b)Whenμ∗=0.0300, the roll system’s time history response and the phase diagram are shown in Figure 7. It can be seen from these figures that the phase trajectory will converges toward a smaller and smaller circle near the origin point. The equilibrium point of work roll is asymptotically stable.Figure 7
Dynamic response of the work roll displacement and the phase diagram underμ∗=0.0300: (a) time history of work roll displacement; (b) phase diagram.
(a)
(b)The roll system stability about the critical parameter ofμ∗=0.0401 is shown in Figure 8. When at the critical Hopf bifurcation parameter, displacement of work roll firstly becomes larger and finally converges to an equal amplitude oscillation. The equilibrium point is unstable, and the Hopf bifurcation diagram presents a ring shape. However, dynamic range of the fixed limit cycle in Figure 8 is a little larger compared with Figure 6. This is mainly because the friction is bigger, and as a result, the fluctuation of the response rolling force ΔPi under the friction coefficient of 0.0401 is much bigger.Figure 8
Dynamic response of the work roll displacement and the phase diagram underμ∗=0.0401: (a) time history of work roll displacement; (b) phase diagram.
(a)
(b)The numerical solution of work roll displacement under the Hopf bifurcation point ofμ∗=0.0420 is shown in Figure 9. The roll displacement will disperse outward all the time. Besides, the phase diagram of the work roll will disperse outward of the cycle, and the radius curvature of the phase trajectory increases gradually. It illustrates that the mill vibration system is unstable and mill chatter is produced.Figure 9
Dynamic response of the work roll displacement and the phase diagram underμ∗=0.0420: (a) time history of work roll displacement; (b) phase diagram.
(a)
(b)Whenμ∗<0.0048, the roll system is unstable, and when μ∗=0.0048, the system is stable. So the Hopf bifurcation of μ∗=0.0048 is the subcritical bifurcation point, and μ∗=0.0401 is the supercritical bifurcation point.Many researchers attribute the change of stability in rolling mill system to the sudden change of the rolling force fluctuation in dynamic equilibrium equation [11, 25, 35]. The optimal setting range of friction coefficient is 0.0048–0.0401, and mill chatter appears when the friction is either too high or too low. This is mainly due to the fluctuation of rolling force. When the friction is high, with the increase of friction coefficient, the neutral angle becomes bigger, and then the neutral point shifts to the entry side, so the distance between the entry and neutral point becomes shorter, and the variation of the entry speed decreases [11]. However, the exit forward slip and exit speed fluctuation will increase. Based on equation (12), the tension variation at the entry position (dTb) and exit position (dTf) is smaller compared with low friction coefficient rolling condition, so the fluctuation of the response rolling force ΔPi increases, and the stability of the stand decreases. Meanwhile, when the friction is quite small, ∂P/∂hout will increase with the decrease of friction coefficient; F4,3 in equation (29) decreases, so the instantaneous dynamic structure spring constant decreases, and the fluctuation of the response rolling force ΔPi caused by the strip exit thickness variation will increase, so the stability of the stand is decreased.


## 4.3. Stable Range Analysis of Other Rolling Parameters

### 4.3.1. Calculation of Limited Rolling Speed

Based on this coupled vertical-horizontal-torsional vibration analysis model, rolling speed has been chosen as the Hopf bifurcation parameter to obtain the limited stability rolling speed of the UCM tandem cold rolling mill, and other rolling conditions are shown in Table1. The actual bifurcation point of rolling speed is 25.1 m/s. Figure 10 shows the time history of work roll displacement under the rolling speed of 25.2 m/s. According to Figure 10, the displacement amplitude of the work roll shows a divergence trend, and the self-excited vibration happens. The critical rolling speed calculated by the proposed chatter model is 25.1 m/s.Figure 10
Work roll time history displacement under the rolling speed of 25.2 m/s.Figure11 shows the tension and roll gap variation in No. 4 and No. 5 stands obtained from the Production Data Acquisition (PDA) system. As is shown in Figure 11(b), when the rolling speed is 22.1 m/s, the back and front tension of the No. 5 stand fluctuate within a small range and the mill stand is stable. However, when the rolling speed is raised to 23.0 m/s, the back and front tension of the No. 5 stand fluctuate more violently, and the front tension of the No. 5 stand begins to diverge as shown in Figure 11(a). It is found that the roll gap of No. 4 and No. 5 stands fluctuates more violently by comparing Figures 11(c) and 11(d). Meanwhile, a sudden and harsh vibration noise will be heard in the No. 5 stand, and finally the strip breaks. It shows mill chatter appears under the speed of 23.0 m/s, and the actual critical rolling speed is 23.0 m/s.Figure 11
Tension and roll gap time history: (a) back and front tension of No. 5 stand under the speed of 23.0 m/s; (b) back and front tension of No. 5 stand under the speed of 22.1 m/s; (c) roll gap of No. 4 and No. 5 stands under the speed of 23.0 m/s; (d) roll gap of No. 4 and No. 5 stands under the speed of 22.1 m/s.
(a)
(b)
(c)
(d)The calculated critical rolling speed is 25.1 m/s, while the actual critical rolling speed is 23.0 m/s. There is a certain relative error between the calculated critical speed and actual critical rolling speed. This is due to the small coefficient difference in Jacobian matrix. With the change of rolling speed, the actual value of the tension and friction is changed. And the calculated value of tension factornt and friction coefficient in ∂P/∂h and ∂M/∂h is a little different with the actual value. So, the calculated limited rolling speed is a little bigger than the actual value. However, by comparing the critical rolling speed calculated by the proposed model with the actual critical rolling speed, it can be seen that the proposed model is still accurate.The vibration is largely influenced by the friction coefficient and rolling speed. Based on equation (13), Figure 12 shows the friction coefficient response of No. 5 stand versus speed. Despite the changes of speed, the friction coefficient is still in the stability range of 0.0048–0.0401. It is divided into two regions. In region 1, both the friction coefficient and rolling speed are in the stable range. In region 2, the friction coefficient is still in the stable range, while the speed exceeds the limited speed range. When the rolling speed is lower than 25.1 m/s, the mill system is stable without changing the lubrication conditions. And when the rolling speed is higher than 25.1 m/s, the stand is unstable. However, the friction coefficient is in the middle value of stable range under the speed is 25.1 m/s, so mill chatter might be eliminated by reducing the basic friction coefficient (by improving the lubrication characteristics, such as increasing the concentration of emulsion, increasing spray flow, and improving the property of the base oil).Figure 12
Friction coefficient response of No. 5 stand versus speed.


### 4.3.2. Calculation of Maximum Reduction

Assuming the entry thickness is constant, the exit thickness is variable. And the exit thickness has been chosen as the Hopf bifurcation parameter. Other calculated rolling conditions are shown in Table1. Then, the calculated bifurcation point of exit thickness is 0.168 mm. The optimal setting range of reduction for No. 5 stand is 0–44%. When the reduction of No. 5 stand is 44%, the increase of reduction leads to a lower mill stand stability. The reason is also related to the increase of rolling force variation under the larger reduction rolling conditions.As can be seen in Table1, in normal rolling schedule, the reduction rate is 33.3%. And the max reduction is up to 44%. In reality, considering there are other factors that influence the stand stability, such as rolling speed, and lubrication conditions, the max reduction is lower than 44%. In general, the effect of the reduction in the mill stability is lower than the improvement of friction coefficient.


### 4.3.3. Front Tension Stable Range

When the work roll velocity is 22.1 m/s, the calculated bifurcation point of exit tension is 1.5–133 MPa. Small front tension means that the rolling force will be quite large. So, under the same excitation interferences, the variation of rolling force will be much bigger, and the stand stability is decreased. When the front tension is big enough, tension factornt decreases, leading to a larger value in ∂P/∂h and ∂M/∂h. As a result, the fluctuation of rolling force and roll torque increases, and the stability of the system is reduced as well. So, the influence of tension on the mill stability is smaller than other parameters.These results can be used to estimate mill chatter in a certain stand. And the stable range of friction, reduction, and tension can be used to adjust the rolling schedule in the UCM mill to avoid mill chatter under high rolling speed.


## 4.3.1. Calculation of Limited Rolling Speed

Based on this coupled vertical-horizontal-torsional vibration analysis model, rolling speed has been chosen as the Hopf bifurcation parameter to obtain the limited stability rolling speed of the UCM tandem cold rolling mill, and other rolling conditions are shown in Table1. The actual bifurcation point of rolling speed is 25.1 m/s. Figure 10 shows the time history of work roll displacement under the rolling speed of 25.2 m/s. According to Figure 10, the displacement amplitude of the work roll shows a divergence trend, and the self-excited vibration happens. The critical rolling speed calculated by the proposed chatter model is 25.1 m/s.Figure 10
Work roll time history displacement under the rolling speed of 25.2 m/s.Figure11 shows the tension and roll gap variation in No. 4 and No. 5 stands obtained from the Production Data Acquisition (PDA) system. As is shown in Figure 11(b), when the rolling speed is 22.1 m/s, the back and front tension of the No. 5 stand fluctuate within a small range and the mill stand is stable. However, when the rolling speed is raised to 23.0 m/s, the back and front tension of the No. 5 stand fluctuate more violently, and the front tension of the No. 5 stand begins to diverge as shown in Figure 11(a). It is found that the roll gap of No. 4 and No. 5 stands fluctuates more violently by comparing Figures 11(c) and 11(d). Meanwhile, a sudden and harsh vibration noise will be heard in the No. 5 stand, and finally the strip breaks. It shows mill chatter appears under the speed of 23.0 m/s, and the actual critical rolling speed is 23.0 m/s.Figure 11
Tension and roll gap time history: (a) back and front tension of No. 5 stand under the speed of 23.0 m/s; (b) back and front tension of No. 5 stand under the speed of 22.1 m/s; (c) roll gap of No. 4 and No. 5 stands under the speed of 23.0 m/s; (d) roll gap of No. 4 and No. 5 stands under the speed of 22.1 m/s.
(a)
(b)
(c)
(d)The calculated critical rolling speed is 25.1 m/s, while the actual critical rolling speed is 23.0 m/s. There is a certain relative error between the calculated critical speed and actual critical rolling speed. This is due to the small coefficient difference in Jacobian matrix. With the change of rolling speed, the actual value of the tension and friction is changed. And the calculated value of tension factornt and friction coefficient in ∂P/∂h and ∂M/∂h is a little different with the actual value. So, the calculated limited rolling speed is a little bigger than the actual value. However, by comparing the critical rolling speed calculated by the proposed model with the actual critical rolling speed, it can be seen that the proposed model is still accurate.The vibration is largely influenced by the friction coefficient and rolling speed. Based on equation (13), Figure 12 shows the friction coefficient response of No. 5 stand versus speed. Despite the changes of speed, the friction coefficient is still in the stability range of 0.0048–0.0401. It is divided into two regions. In region 1, both the friction coefficient and rolling speed are in the stable range. In region 2, the friction coefficient is still in the stable range, while the speed exceeds the limited speed range. When the rolling speed is lower than 25.1 m/s, the mill system is stable without changing the lubrication conditions. And when the rolling speed is higher than 25.1 m/s, the stand is unstable. However, the friction coefficient is in the middle value of stable range under the speed is 25.1 m/s, so mill chatter might be eliminated by reducing the basic friction coefficient (by improving the lubrication characteristics, such as increasing the concentration of emulsion, increasing spray flow, and improving the property of the base oil).Figure 12
Friction coefficient response of No. 5 stand versus speed.


## 4.3.2. Calculation of Maximum Reduction

Assuming the entry thickness is constant, the exit thickness is variable. And the exit thickness has been chosen as the Hopf bifurcation parameter. Other calculated rolling conditions are shown in Table1. Then, the calculated bifurcation point of exit thickness is 0.168 mm. The optimal setting range of reduction for No. 5 stand is 0–44%. When the reduction of No. 5 stand is 44%, the increase of reduction leads to a lower mill stand stability. The reason is also related to the increase of rolling force variation under the larger reduction rolling conditions.As can be seen in Table1, in normal rolling schedule, the reduction rate is 33.3%. And the max reduction is up to 44%. In reality, considering there are other factors that influence the stand stability, such as rolling speed, and lubrication conditions, the max reduction is lower than 44%. In general, the effect of the reduction in the mill stability is lower than the improvement of friction coefficient.


## 4.3.3. Front Tension Stable Range

When the work roll velocity is 22.1 m/s, the calculated bifurcation point of exit tension is 1.5–133 MPa. Small front tension means that the rolling force will be quite large. So, under the same excitation interferences, the variation of rolling force will be much bigger, and the stand stability is decreased. When the front tension is big enough, tension factornt decreases, leading to a larger value in ∂P/∂h and ∂M/∂h. As a result, the fluctuation of rolling force and roll torque increases, and the stability of the system is reduced as well. So, the influence of tension on the mill stability is smaller than other parameters.These results can be used to estimate mill chatter in a certain stand. And the stable range of friction, reduction, and tension can be used to adjust the rolling schedule in the UCM mill to avoid mill chatter under high rolling speed.


## 5. Conclusion

(1)
A new multiple degrees of freedom and nonlinear dynamic vertical-horizontal-torsional coupled mill chatter model for UCM cold rolling mill is established. In the dynamic rolling process model, the time-varying coupling effect between rolling force, roll torque variation, and fluctuation of strip thickness and tension is formulated. In the vibration model, the movement of work roll in both vertical-horizontal-torsional directions and the multiroll characteristics in the UCM rolling mill are considered. Finally, the dynamic relationship of roll movement responses and dynamic changes in the roll gap is formulated.
(2)
Selecting the friction coefficient as the Hopf bifurcation parameter and by using Routh–Hurwitz determinant, the Hopf bifurcation point of the mill vibration system is calculated, and the bifurcation characteristics are analyzed. Research shows that either too low or too high friction coefficient increases the probability of mill chatter. When the friction coefficient is extremely low, the mill system produces a stable limit cycle. With the friction continuing to increase, the phase trajectory will converge toward the origin point. However, the limit cycle of phase diagram disappears and disperses outward under high friction rolling conditions.
(3)
Based on the coupled mill chatter model, the mill system’s stable limited rolling speed, stable range of reduction, and front tension are investigated. Without changing the rolling speed, the most effective way to eliminate mill chatter is to change the lubrication conditions and adjust the reduction schedule. And the effect of tension on the vibration is quite limited. These methods can be used to improve the chatter phenomena in rolling mill under high rolling speed. The proposed vibration model can also be used to design an optimum rolling schedule in the UCM tandem cold rolling mill and to determine the appearance of mill chatter under certain rolling conditions.



---

*Source: 4358631-2019-02-17.xml*