Hot strip mill nonlinear torsional vibration with multistand coupling
Xiaobin Fan^{1} , Yong Zang^{2} , Feng Wang^{3} , Ke Jin^{4}
^{1, 3, 4}School of Mechanical and Power Engineering, Henan Polytechnic University, Jiaozuo 454003, China
^{1, 2}Mechanical Engineering Institute, University of Science and Technology Beijing, Beijing 100083, China
^{1}Corresponding author
Journal of Vibroengineering, Vol. 17, Issue 4, 2015, p. 16231633.
Received 18 February 2015; received in revised form 22 April 2015; accepted 4 June 2015; published 30 June 2015
JVE Conferences
Strip mill vibration greatly affects the quality of products, rolling productivity, and even cause strip broken or equipment damage. In order to explore the serious vibration phenomena occured in CSP (Compact Strip Production) hot strip mill when rolling thin container strips, in this paper we established a mill main drive system nonlinear torsional vibration model, and emphatically studied the influence of rolling load, multistand coupling (tension) and interface friction on the torsional vibration of the system; Based on numerical simulation, field measurement and signal timefrequency analysis (WignerVille distribution), we analyzed the mills torsional vibration and horizontal vibration characteristics, the influence of multistand coupling (tension) on the rolls torsional vibration.
Keywords: rolling mill, torsional vibration, coupling, timefrequency analysis, nonlinear.
1. Introduction
Hot strip mill is a large complex integration equipment, there are a variety of complex coupling between systems and the process of rolling mill, in certain conditions or disturbances of rolling process, it is easy to induce varying degrees vibration and these vibration problems severely limit the quality and productivity of rolling strip products, and severe vibration may cause strip broken or equipment damage, which has serious threat to safety and causes huge economic losses. Despite there have extensive research, but due to the complexity of the issue, which has not been solved perfectly.
A. Bar presented the numerical analysis of vibrations with moderate frequencies excited during the milling process taking into account the transportation motion of the strip and its inertial properties [1]. Yongjiang Zheng carried out a complete spatial vibration characteristic analysis based on the modified Riccatioverall transfer matrix method and coupled vibration characteristics of cross and sway motions along with the complex mode shape analysis of the mills [2]. Adam Bar built the nonlinear mathematical model of oscillated system: the continuous group of the rolling stands coupled by transferring strip [3]. PeiHua studied threestand tandem mill stability in terms of the regenerative effect by a multistand chatter model that incorporates a homogeneous process model and a mill structural dynamics model [4]. Y. X Wu found the mill vibration is caused by amplitude and frequency modulation by the analysis of the mill roller bearings, the periodic change of stiffness of the bearings at the fixed load orientation, together with the periodic change of the effective errors of the bearings [5]. Vladimir Panjković finded that the frictional conditions in the roll gap are the principal cause of chatter in this mill and the frictional conditions appear to be associated with the thickness and properties of oxide formed on rolls [6]. Lipo Wang apply the dynamically generated fuzzy neural network (DGFNN) controller to a realworld application of controlling the torsional vibration of tandem coldrolling mill spindles [7]. Hideaki Furumoto designed a chamber in Mill Stabilizing Device and optimized its size [8].
The ${F}_{3}$ of the CSP hot strip mills group occurs strong vibration when rolling thin container srip (SPAH), ${F}_{2}$ and ${F}_{4}$ will also occur vibration when increasing its load correspondingly. The CSP hot rolling production line has six fourroll finishing mill and the production line layout is shown in Fig. 1. The work rolls of ${F}_{1}$${F}_{3}$ have large roller diameter (720800 mm, high chromium cast iron), ${F}_{4}$${F}_{6}$ are small diameter (540600 mm, unlimited chilled cast iron). The slab thickness is generally 5070 mm, the first three mills use large reduction and the reduction rate is up to 50 %.
Fig. 1. production line layout of CSP hot strip mill
2. Nonlinear torsional vibration characteristics of mill main drive system
2.1. Nonlinear torsional vibration model of main driving system
The dynamics model of mill main driving system is established with two inertia and nonlinear damping and nonlinear stiffness under alternating load, which is shown in Fig. 2. The mill rolls are dragged by AC synchronous motor through flexible shaft of the main driving system. It is simplified as twoinertia system. Let ${J}_{M}$, ${J}_{L}$ be the inertia moment of the motor and the load respectively. Its rotational speed are ${\omega}_{M}$ and ${\omega}_{L}\text{.}$ Angle are ${\theta}_{M}$ and ${\theta}_{L}$ respectively. The electromagnetic motor torque is ${T}_{M}$. The load torque is ${T}_{L}$. Shafting torsional damping function is ${C}_{ML}$ and expressed as a Van der Pol oscillator forms, namely ${C}_{ML}=C\left(1+\gamma {\left({\dot{\theta}}_{M}{\dot{\theta}}_{L}\right)}^{2}\right)$, where $\gamma $ is variable coefficient affected by rolling conditions; Shafting torsional stiffness function is denoted by ${K}_{ML}$ and expressed in the form of duffing oscillators, namely ${K}_{ML}=K\left(1\left({\left({\theta}_{M}{\theta}_{L}\right)}^{2}/6\right)\right)$; ${C}_{M}$ is the motor damping coefficient, ${C}_{L}$ is the roll (roll gap) damping coefficient.
Fig. 2. Nonlinear dynamic model of mill torsional vibration
Its dynamic model can be written as:
and:
where, $R$ is the radius of the work roll, $\mu $ is the coefficient of friction between the roll strip (detailed calculation is as Section 2.3), $P$ is the rolling force. The rolling force $P$ is composed of steady state rolling force ${P}_{0}$ and dynamic rolling force, and can be expressed as:
where, $\kappa $ is coefficient of rolling power fluctuations, which can be drawn by Section 2.4; ${\omega}_{r}$ is the angular velocity of the rotating roll, ${z}_{a}$ is the number of spindle teeth; The calculation of the steadystate rolling force ${P}_{0}$ is shown in Section 2.2.
Electromagnetic motor torque is ${T}_{M}={T}_{M0}(1+\xi \mathrm{c}\mathrm{o}\mathrm{s}{\omega}_{m}t)$, where, ${T}_{M0}$ is steadystate value of the motor electromagnetic torque, $\xi $ is the fluctuation coefficient, ${\omega}_{m}$ is the fluctuation frequency.
${T}_{B}$ and ${T}_{F}$ are the back and front tension respectively. And there are:
${T}_{B}={T}_{B0}\left(1+{\zeta}_{B}\mathrm{c}\mathrm{o}\mathrm{s}{\omega}_{TB}t\right),$
where, ${T}_{F0}$ and ${T}_{B0}$ are steadystate value of back and front tension. ${\zeta}_{F}$ and ${\zeta}_{B}$ are fluctuation coefficient of back and front tension; ${\omega}_{TF}$ and ${\omega}_{TB}$ are fluctuation frequency of back and front tension and its value can be passed around the mill measured vibration signal.
2.2. Steady rolling force calculation
Applying OROWAN deformation zone equilibrium theory, the hot rolled steady rolling force ${P}_{0}$ is calculated as follows:
where, $B$ is the strip width, ${l}_{c}^{\text{'}}$ is the horizontal projection length (mm) of the contact arc between flattening rolls and steel strip. Using the Hitchikok elastic flattening equation, we have:
where, $r$ is the absolute reduction and $m=8(1{\nu}^{2})/\pi {E}_{R}$ (${E}_{R}$ is roll elastic modulus, $\nu $ is Poisson coefficient).
${Q}_{P}$ is impact coefficient of stress state affected by the friction on the arc of contact. And which is calculated as follows by the regression method:
where, $h{}_{m}{}^{\mathrm{}}\mathrm{}=({h}_{1}+{h}_{2})/2$, and $\epsilon $ is the reduction rate, ${h}_{1}$ and ${h}_{2}$ are strip thickness of entrance and exit.
${K}_{q}$ is the impact factor of the rolling force caused by the back and front tension. Because of the hot strip mill rolling conditions at small tension and the fact that back tension has more influence than front tension, ${K}_{q}$ can be calculated as the following simplified formula:
where, ${q}_{f}$ and ${q}_{b}$ are front and back tension as before. $K$ is deformation resistance (Mpa) and its regression formulas of container strip deformation resistance is as follows:
where, strain rate:
and $N$ is roll speed, r/min; $T$ is rolling temperature.
2.3. Analysis of stand coupling (tension) effect
Mills coupled vibration is an important reason for the mill vibration instability. The mills distance and rolling speed and other rolling parameters has an important impact on the stability of the system. According to Hookeps law, the ${F}_{i}$ entrance tensile stress is proportional to the integral value of the difference between ${F}_{i}$ entrance velocity and ${F}_{i1}$ exit velocity. The tensile stress of ${F}_{i}$ mill are as follows:
where, ${E}_{s}$ is Young’s modulus of strip; ${L}_{i1}$ is the distance between ${F}_{i}$ and ${F}_{i1}$ mill; ${v}_{1}$ and ${v}_{2}$ are entrance velocity and exit velocity respectively. Relations between the stand are shown in Fig. 3.
Fig. 3. Multistand relation diagrams
According to the principle of equal mass flow, entrance velocity of rolling is as follows:
where, ${h}_{2m}$ is the steady state value of exit strip thickness, ${A}_{0}$ and ${\omega}_{h}$are the amplitude and frequency of rolls vertical vibration (all obtained from measured data). And the change of the back tensile stress is as follows:
where, $L$ is the distance between the mills. And the rolling force variation caused by $\mathrm{\Delta}{q}_{b}$ can be written as follows:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\frac{1.4{E}_{S}{v}_{2}}{L{h}_{1}{\omega}_{h}^{2}}B\sqrt{Rr}\left(0.8205+0.2376\frac{{l}_{c}}{{h}_{m}}+0.1006\epsilon \frac{{l}_{c}}{{h}_{m}}0.3768\epsilon \right)\dot{A}\left(t\right).$
From Eq. (14), we see that the phase of rolling force change $\mathrm{\Delta}{P}_{qb}$ caused by back tension in vibration is ahead the roll vertical vibration displacement of 90°phase angle.
2.4. Characteristics of rolling interface friction
Let us consider lubrication rolling. The interface is full of fluid lubrication and the strip is broadband steel ($l/\stackrel{}{h}>\text{5}$. Let $l$ is contact arc length, $\stackrel{}{h}$ be the average thickness of strip entrance and exit). Then, the onedimensional Reynolds equation can be established as follows for steady and unsteady rolling process:
where $x$ is the distance along the rolling direction, $h$ is the film thickness, $p$ is oil film pressure on the length of the lubricating wedge before the deformation zone entrance, $\stackrel{}{U}$ is average of work roll line speed and the strip velocity, $t$ is the time, ${\eta}_{0}$ is the dynamic viscosity of the lubricant, and it is not a constant for unsteady hot rolling process, which can be expressed as ${\eta}_{0}(T,p)={\eta}_{A}{e}^{\theta p\delta (T{T}_{s})}$, where ${\eta}_{A}$ is viscosity at room temperature ${T}_{s}$ and atmospheric pressure, θ is the viscosity pressure coefficient, $\delta $ is viscosity temperature coefficient.
From Eq. (15), we see that film pressure variation is composed by the squeeze effects $\partial h/\partial t$ and hydrodynamic effects $\partial /\partial x\left(\stackrel{}{U}h\right)$. In the steady rolling process there are no film thickness fluctuations. Thus, the impact of the former can be ignored and only the latter be considered. Under this condition, the film pressure $p\left(x\right)$ at any position x can be deduced as follows [9]:
$\left.\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\frac{1}{x}\left(1+\frac{3R{\xi}_{0}}{2{A}^{2}}\right)\left(1+\frac{{A}^{2}}{3{x}^{2}}\right)\right\},$
where, ${v}_{r}$ is line speed of work roll, ${v}_{1}$ is strip velocity, ${A}^{2}=R(\mathrm{\Delta}h2{\xi}_{0})$, $R$ is roll radius, lubricating layer thickness at the entrance is ${\xi}_{0}=3\theta {\eta}_{0}({v}_{r}{v}_{1})/\alpha [1{e}^{\theta (K{\sigma}_{0})}]$ (where $\alpha $ is the bite angle, $K$ is yield strength of the strip, ${\sigma}_{0}$ is back tensile stress of strip); The film shear stress $\tau \left(x\right)$ at an arbitrary position $x$ in the deformation zone is as follows [9]:
where, fluid resistance factor is $J=6{v}_{r}{\eta}_{0}/{\xi}_{0}\text{,}$ fluid lubrication coefficient is $I=\left({\eta}_{0}{l}_{d}^{2}{v}_{r}\right)/2{p}_{m}{\xi}_{0}^{2}$, $z={x}_{\varphi}/{l}_{d}$, ${x}_{\varphi}$ is neutral point coordinates, ${l}_{d}$ is contact arc length, $\epsilon $ is pass reduction rate. By the definition of friction factor $\mu \left(x\right)=p\left(x\right)/\tau \left(x\right)$, and according to the Eqs. (16) and (17), the friction coefficient average of the deformation zone can be approximated as follows [9]:
where, ${h}_{1}$ is entry thickness of strip, $D$ is the roll diameter. For the convenience of studying, it usually is taken as $\mu =a{e}^{bR{\omega}_{L}+c}$ or polynomial form reexpressed by the Taylor series expansion.
2.5. Spindle unbalanced excitation
The length is about 5 m and the weight is about 7 t of the curved tooth spindles. There are no support; if there exists eccentricity of the inner and outer gear and the imbalance caused by the shaft gravity center offset due to its weight or other factors, then great unbalanced force will appear. This force will damage the dynamic stability. When there is curved tooth headspace, the centrifugal force ${F}_{c}$ is as follows:
where, ${m}_{1}$ is the weight of the outer tooth sleeve, kg; $n$ is the operating speed of couplings, r/min; ${\mathrm{\Delta}}_{y}$ is relative radial displacement, mm. Due to the centrifugal force, in order to make the inner and outer tooth couplings can automatically make centering, the minimum torque ${T}_{min}$ is needed as follows:
where, ${m}_{2}$ is the weight of the intermediate shaft or the inner ring gear, kg; $e$ is relative eccentricity of inner and outer teeth, mm; $n$ is coupling rotational speed, r/min; $d$ is tooth pitch diameter, mm. Substituting these parameters into Eq. (20), there is ${T}_{min}=$4.52×10^{7} N∙m, but the measured torque is about 5×10^{5} N∙m, which cannot meet the requirements for automatic centering, that is there must have unbalanced centrifugal force.
2.6. Numerical simulation of torsional vibration mill
The ${F}_{3}$ mill structure parameters and process parameters values of the CSP hot strip line are as follows:
$C=$ 1.2×10^{4} N∙m/(rad∙s^{1}), ${\omega}_{TB}=$14∙2$\pi $ Hz, ${\omega}_{TF}=$76∙2$\pi $ Hz, ${\omega}_{r}=$ 0.8∙2$\pi $ Hz,
${T}_{B0}=$ 7.3×10^{5} N, ${T}_{F0}=$ 3.8×10^{5} N.
The influence of front and back tension (ie, adjacent mill) and strip resistance force on main driving system torsional vibration is simulated emphatically in this paper. The numerical simulation model was set up in Matlab/Simulink using 4order and 5order RungeKutta method. The influence of rolling material or specification were analyzed by Eq. (5).
The simulation results of the influence of the front and back tension on torsional vibration (main driveline torque) is shown in Fig. 4. We see that mill torsional vibration becomes serious when there is no back tension; the torque mean has increased when there is no front tension. Under these conditions the vibration may increase.
The influence of strip resistance force, ie, different strip materials or specifications, on the main driveline torque is shown in Fig. 5. It can be seen that if rolling containers strip (SPAH), the torque increases significantly. This is an important factor of torsional vibration startup.
3. Field experiment analysis
For ${F}_{2}$, ${F}_{3}$ and ${F}_{4}$, rolling mill main driving torque was tested using strain measurement and the signal is transmitted by the slip ring; The roll vibration signal was got by piezoelectric accelerometer installed on roller bearing housing; From the site, ${F}_{2}$${F}_{4}$ produce different vibration levels when rolling container strip (SPAH, finishing temperature is 900°C) with finishing thickness 1.6 mm. The parameters of the processes and force are shown in Table 1 for ${F}_{2}$, ${F}_{3}$ and ${F}_{4}$ mill. The working roll diameter is shown in Table 2. The chatter marks obtained by emergency stop sampling of strip and rollers are shown in Fig. 6.
Fig. 4. The influence of the front and back tension on torsional vibration
Fig. 5. The influence of strip materials or specifications on the main driveline torque
Table 1. Rolling processes and force parameters
${F}_{2}$

${F}_{3}$

${F}_{4}$


Entrance width (mm)

1142.54

1142.54

1142.54

Entrance thickness (mm)

23.911

10.546

5.265

Entrance temperature (°C)

991.53

971.53

951.89

Exit temperature (°C)

993.16

972.67

954.26

Contact arc length (mm)

74.9

46.7

26.1

Roller temperature (°C)

37.83

38.66

39.46

Strip hardness (kN)

46965

47298

37011

Mean rolling force (kN)

29503

25165

16080

Mean rolling torque (kN·m)

1663.7

861.8

296.8

Maximum rolling torque (kN·m)

19171.8

11544.2

5009.5

Rolling power (kW)

6069

6567

5241

Back tension (N/mm^{2})

109.28

60.24

60.15

Rolling velocity (m/s)

1.45

2.87

4.81

Bending force (kN)

1400.1

1400

1099.1

Roll axial movement (mm)

67.87

35.01

29.42

Roll speed (rps)

0.588

0.8

Table 2. Roll diameters
${F}_{2}$

${F}_{3}$

${F}_{4}$


Upper work roll diameter (mm)

791.55

769.232

549.1302

Lower work roll diameter (mm)

791.492

769.154

549.0517

It is found from the strip and roll chatter marks that the strip chatter marks space behind ${F}_{3}$ and under the strip surface is 5 cm. The chatter space is 4.67 cm above the strip surface.The chatter marks have same phase, i.e., the under and above are aligned. The frequency of them is about 51 Hz; the space of strip chatter marks behind ${F}_{2}$ is 13.89 cm. The phase of the lower part and upper part are inverted.
The time domain and timefrequency distribution (WignerVille distribution) of ${F}_{2}$ working roll torque is shown in Fig. 7. ${F}_{3}$ started biting strip at 3.22 seconds. It can be seen that the ${F}_{2}$ work rolls has slight torsional vibration. The dominant frequency is about 17 Hz, ${F}_{2}$ torque increased after ${F}_{2}$ biting strip. ${F}_{2}$ working roll horizontal vibration acceleration is shown in Fig. 8. Its dominant frequency is about 15 Hz.The frequency spectral energy of horizontal vibration is more concentrated than driveline torque, namely ${F}_{2}$ main vibration modes is roll horizontal vibration. The ${F}_{2}$ roll horizontal acceleration amplitude increases a little after ${F}_{3}$ biting strip, and the frequency spectrum diverse slightly.
Fig. 6. Strip and roll chatter marks
Fig. 7. ${F}_{2}$ torque time series and WignerVille distribution (${F}_{3}$ biting strip at 3.22 s)
Fig. 8. ${F}_{2}$ roll horizontal acceleration time series and WignerVille distribution (${F}_{3}$ biting strip at 3.22 s)
${F}_{3}$ work roll torque time series and WignerVille distribution is shown in Fig. 9. Strip left from ${F}_{2}$ at 4.78 seconds. It can be seen that the spectrum is dispersed, the firstorder relative advantage frequency is about 20 Hz. ${F}_{3}$ torque amplitude increases significantly after strip left from ${F}_{2}$, namely the absence of back tension will worsen the vibration. This coincides with the aforementioned simulation results; ${F}_{3}$ work rolls horizontal acceleration time series and WignerVille distribution are shown in Fig. 10. We see that it has “gourd” shape. So there exists obvious vibration. The spectrum has absolute maximum peak at 51 Hz, and coincides with the roll and the strip chatter marks frequencies. ${F}_{3}$ horizontal vibration amplitude significantly increases after strip left from ${F}_{2}$. Therefore, that is, the back tension has an important role in suppressing the vibration.
Fig. 9. ${F}_{3}$ torque time series and WignerVille distribution (strip left from ${F}_{2}$ at 4.78 s)
Fig. 10. ${F}_{3}$ horizontal acceleration time series and WignerVille distribution (strip left from ${F}_{2}$ at 4.78 s)
${F}_{3}$ roll torque time series and WignerVille distribution are shown in Fig. 11, and ${F}_{4}$ started biting strip at 1.64 seconds, we see that ${F}_{3}$ torque amplitude decreased slightly after ${F}_{4}$ biting strip and coincided with the aforementioned simulation results. The spectrum is dispersed, that is the front tension has suppression effect on ${F}_{3}$ torsional vibration. The ${F}_{3}$ work roll horizontal acceleration time series and WignerVille distribution are shown in Fig. 12. We can see that ${F}_{3}$ acceleration amplitude has been reduced after ${F}_{4}$ biting strip, that is the front tension has suppression effect on ${F}_{3}$ horizontal vibration and dominant frequency has no obvious change.
Fig. 11. ${F}_{3}$ roll torque time series and WignerVille distribution (${F}_{4}$ started biting strip at 1.64 s)
Fig. 12. ${F}_{3}$ roll horizontal acceleration time series and WignerVille distribution (${F}_{4}$ started biting at 1.64 s)
4. Conclusions
1) The two degrees of freedom mill drive system nonlinear torsional vibration model with nonlinear stiffness and nonlinear friction damping was established. The calculation method and its effect on the torsional vibration system of the rolling force, multistand coupling (tension), the interfacial friction characteristics and spindle unbalanced incentives and other factors are discussesd in this paper.
2) By the numerical simulation results analysis of the model, it is shown that the front and back tension have suppression effect on mill vibration. When there is no back tension, the mill torsional vibration is obvious. At the same time rolling resistance force is an important factor responsible for mill torsional vibration.
3) Based on longterm tracking test and signals timefrequency analysis, the experimental results show that the tension plays a very good role on mill vibration suppression.
References
 Bar A., Bar O. Types of midfrequency vibrations appearing during the rolling mill operation. Journal of Materials Processing Technology, Vol. 162163, Issue 5, 2005, p. 461464. [Search CrossRef]
 Zheng Yongjiang, Xie Zhaohui, Li Yigeng, Shen Guangxian Spatial vibration of rolling mills. Journal of Materials Processing Technology, Vol. 213, Issue 4, 2013, p. 581588. [Search CrossRef]
 Bar Adam, Świątoniowski Andrzej Interdependence between the rolling speed and nonlinear vibrations of the mill system. Journal of Materials Processing Technology, Vol. 155156, Issue 11, 2004, p. 21162121. [Search CrossRef]
 Hu PeiHua, Ehmann Kornel F. Regenerative effect in rolling chatter. Journal of Manufacturing Processes, Vol. 3, Issue 2, 2001, p. 8293. [Search CrossRef]
 Wu Y. X., Duan J. A. Frequency modulation of highspeed mill chatter. Journal of Materials Processing Technology, Vol. 129, Issue 13, 2002, p. 148151. [Search CrossRef]
 Panjković Vladimir, Gloss Ronald, Steward John, Dilks Stephen, Steward Robert Causes of chatter in a hot strip mill: observations, qualitative analyses and mathematical modelling. Journal of Materials Processing Technology, Vol. 212, Issue 4, 2012, p. 954961. [Search CrossRef]
 Wang Lipo, Frayman Yakov A dynamically generated fuzzy neural network and its application to torsional vibration control of tandem cold rolling mill spindles. Engineering Applications of Artificial Intelligence, Vol. 15, Issue 6, 2002, p. 541550. [Search CrossRef]
 Furumoto Hideaki, Kanemori Shinya, Hayashi Kanji, Sako Akira Enhancing technologies of stabilization of mill vibration by mill stabilizing device in hot rolling. Procedia Engineering, Vol. 81, 2014, p. 102107. [Search CrossRef]
 Xiaobin Fan Vibration Problem Research for CSP Mill Stand. University of Science and Technology Beijing, Beijing, 2007. [Search CrossRef]