Bifurcation characteristics of torsionalhorizontal coupled vibration of rolling mill system
Shuang Liu^{1} , Yatao Shi^{2} , Yunpeng Zhang^{3} , Meng Zong^{4}
^{1, 2, 3, 4}Key Laboratory of Industrial Computer Control Engineering of Hebei Province, Yanshan University, Qinhuangdao, China
^{1, 4}National Engineering Research Center for Equipment and Technology of Cold Strip Rolling, Yanshan University, Qinhuangdao, China
^{2}Corresponding author
Journal of Vibroengineering, Vol. 19, Issue 3, 2017, p. 21882201.
https://doi.org/10.21595/jve.2016.17493
Received 1 August 2016; received in revised form 27 December 2016; accepted 29 December 2016; published 15 May 2017
JVE Conferences
In this paper, the static bifurcation and dynamic bifurcation of rolling mill system with torsionalhorizontal coupling vibration are studied. Firstly, the dynamic equation of torsionalhorizontal coupled vibration of rolling mill main drive system is established. Rolling mill is driven by the gear pair, so the torsional vibration of the gear pair, the horizontal vibration and the friction factor are considered in the dynamic equation. Then the equivalent lowdimensional bifurcation equation, which can reveal the system nonlinear characteristics, is obtained using LyapunovSchmidt method, and the system static bifurcation characteristics are studied using singularity theory. Lastly, the bifurcation condition and stability of the system with dynamic Hopf bifurcation are studied using Hopf bifurcation theorem. Numerical simulation of the actual parameter values confirms the analytical results.
Keywords: rolling mill, lubrication friction, LyapunovSchmidt method, singularity theory, bifurcation.
1. Introduction
With the development of modern steel rolling industry, rolling mill is required to have high precision and high dynamic performance. However, unstable vibration has been recognized as a major restriction in rolling productivity [1], so many scholars have been devoted to the study of rolling mill vibration. A large number of researchers have studied rolling mill system from the perspective of the torsional vibration. Liu H. [2] established a dynamic torsional vibration equation of a rotary nonlinear dynamic system and discussed the dynamic stability by the Hopf bifurcation theory. Scholars have made great achievements in the study of the mechanical systems dynamics. Shen Y. J. [3, 4] used singularity theory to analyze the dynamic response of a semiactive control system. Yang S. P. [5, 6] analyzed a vehicle suspension system’s singularity completely and studied the effects of system parameters on the bifurcation parameters and unfold parameters in detail. By theoretical analysis, it is shown that the design of parameters has a close relation with the system’s stability. According to the situation of bifurcation parameters, the reasonable parameters can achieve system stability [7, 8]. Li X. [9] studied the bursting phenomenon in a piecewise mechanical system with different time scales and provided important theoretical basis on the mechanical manufacturing and engineering practice. Shi P. [10], Liu G. [11] and Liu S. [12] established a nonlinear torsional vibration model of rolling mill system, respectively.
At the same time, some scholars have studied the dynamic behavior of rolling mill from the perspective of vertical and horizontal vibration. Yang X. [13] presented a vertical vibration model based on lubrication and friction theory. Hou D. [14, 15] established a verticalhorizontal coupling vibration dynamic model of rolling mill. It was found that the reasonable parameters can effectively restrain the bifurcation of system. Xu Y. [16] established a dynamic coupling model of the fourroll mill in horizontal and vertical. The stability of cold rolling mill system is discussed under different working conditions. Considering the friction between the roller and the workpiece, Zhang R. [17] derived a horizontal motion nonlinear equation of the singleroll driving mill and analyzed the effect of the stiffness fluctuation on the horizontal nonlinear resonance of rolling mill.
The purpose of this paper is to study horizontaltorsional coupled vibration characteristics of rolling mill system. As far as we know, rolling mill vibration is influenced by several coupling factors, which should be taken into account as much as possible. It was found that the rolling speed, the nonlinear friction and the thickness of the strip have strong influences on the stability of rolling mill [1820]. The research also found that there are a lot of horizontaltorsional coupling vibration phenomena in rolling mill [21]. In the actual rolling process, the horizontal vibration and torsional vibration are closely related to the change of friction force. Under certain conditions, these vibrations are coupled with each other. Therefore, it is necessary to study these two kinds of coupling vibration characteristics of rolling mill. In this paper, the influence of variable friction factors is considered and a horizontaltorsional coupled vibration model of rolling mill is established. Based on this model, dynamic behavior of system is researched and simulated. Besides, it reveals the characteristics and rules in instability of vibration system, and promotes the development of the research on rolling mill coupled vibration system.
The paper is organized as follows. Section 2 establishes a torsionalhorizontal coupling vibration dynamic model of rolling mill system. The static bifurcation of the system is studied, and the topological structure of the static bifurcation is obtained in Sections 3. Section 4 studies the system dynamic Hopf bifurcation, and verifies dynamic behavior by numerical simulation. Section 5 concludes this paper.
2. The coupling vibration model with friction of rolling mill
In the study of rolling mill driven by the motor, in order to simplify the model and conform to actual physical meaning, rolling mill drive system is usually abstracted into two parts, which are the motor terminal and the roller terminal according to lumped mass method. In rolling mill drive system, the motor and the roller are connected through the driving system, which consist of gearwheel box, universal couplings, arcgear synchronous belt and so on. For reflecting main influencing factors of the actual process, the connecting part can be regarded as a gear pair system, which is composed of the active gear and the passive gear. Because the motor and the driving gear, the roller and the passive gear are connected through a rigid shaft, the motor and the driving gear, the roller and the passive gear can be regarded as a concentrated mass, respectively. The model of rolling mill main drive system with the gear pair is as follows.
Fig. 1. The model of rolling mill main drive system with the gear pair
${I}_{p}$, ${I}_{g}$ are the rotational inertia of the input and the load after centralized mass, respectively. ${T}_{p}$, ${T}_{g}$ are the input torque and the load torque, respectively. Since the input torque is volatile, its general form is ${T}_{p}={T}_{a}+{T}_{b}\mathrm{c}\mathrm{o}\mathrm{s}\left(\omega \stackrel{}{t}\right)$. ${\theta}_{p}$, ${\theta}_{g}$ are the vibratory angular displacements of the input and the load, respectively. ${c}_{m}$, ${k}_{m}$ are the meshing damping coefficient, the meshing stiffness, respectively. ${R}_{p}$, ${R}_{g}$ are the base radius of the input and the load, respectively. ${c}_{b}$, ${k}_{b}$ are the bearing damping and the bearing stiffness on the driven shaft, respectively. ${F}_{b}$ is the force between the bearing and the driven gear on the drive shaft ${m}_{g}$, ${x}_{g}$ are the mass of driven gear, the horizontal displacement of the load center, respectively ${v}_{r}$, ${v}_{0}$ are the roller linear speed, the stable workpiece speed, respectively. ${F}_{f}$ is the friction between the roller and the workpiece.
In the actual working process of rolling mill, the roller is not often in contact with the workpiece directly, and it is necessary to do some lubrication. So, there is a lubricant film between the roller and the workpiece, and the friction force is a kind of lubrication friction force. Studies have shown that there is a correspondence between the friction coefficient and the thickness of the oil film. It can be affected by a variety of factors, such as the rolling speed, cold reduction ratio, the rolling out entry tension and so on. Base on the theory of fluid mechanics and Hill's formula, simplified the dynamics process and emphasized the influence of the lubricant viscosity coefficient $\eta $ and the speed of the roll line ${v}_{r}$, the model of oil film thickness under certain conditions can be simplified as:
Among them, ${b}_{0}$, ${c}_{0}$ are the corresponding coefficients.
In the rolling process, the corresponding data of film thickness and friction factor can be obtained through appropriate measuring and calculating, and the expression of the friction factor can be acquired by fitting the data. Expressing the friction coefficient in the form of index, the friction coefficient can be described as:
where, ${a}_{0}$, $a$, $b$ are the liquid friction coefficient, the dry friction coefficient, the corresponding coefficient of the friction, respectively. When the roller is in horizontal vibration, the linear speed of the roller ${v}_{r}$ can be expressed as:
where ${v}_{0}$ is the stable workpiece speed; ${\dot{x}}_{g}$ is the horizontal vibration speed. Substituting Eq. (3) into Eq. (2) and performing Taylor expansion, expanding retention in first several items have:
where $d=b\eta $, and ${\mu}_{0}=a{e}^{d{v}_{0}+c}$ is the coefficient of friction when the rolling speed is ${v}_{0}$. At this point the friction between the roller and the workpiece is:
where ${F}_{l}$ represents the rolling force.
According to the Newton’s law, the dynamic equation of the system is:
The transmission error of the linear displacement is defined as:
Substituting Eq. (6) with the linear displacement transfer error $z$, and considering the friction force ${F}_{f}$ in the formula, further simplified results can be expressed as:
where ${m}_{e}={I}_{p}{I}_{g}/\left({I}_{p}{R}_{g}^{2}+{I}_{g}{R}_{p}^{2}\right)$, ${F}_{1}={T}_{a}/{R}_{p}={T}_{g}/{R}_{g}$, ${F}_{2}={m}_{e}{R}_{p}{T}_{b}/{I}_{p}$.
Eq. (8) is the dynamics equation of rolling mill main drive system, which takes horizontaltorsional coupled vibration, lubrication friction and other factors into account. The dynamic behavior of rolling mill is very complex. The research shows that there is obvious horizontal chatter in the rolling process, and source of horizontal vibration is a tribological behavior between the coexistence of fluid lubrication and dry friction motion in rolling interface. Thus, taking the roller horizontal vibration and torsion vibration coupling effect into account, the model is built. Then the friction factor between roller and workpiece is analyzed in detail. Regarding the friction as the main factor, the coupling vibration of rolling mill is analyzed and the influence of parameter variation on vibration buckling is given. It can provide relative guidance for the design of actual system. In the following, the friction coefficient d would be taken as the main research parameter to analyze dynamic behavior of rolling mill.
3. Static bifurcation analysis of high dimensional nonlinear rolling mill system
Selecting ${z}_{0}$ as the nominal scale, and letting:
${k}_{1}=1,{P}_{1}=\frac{{F}_{1}}{\left({z}_{0}{m}_{e}{\omega}_{1}^{2}\right)},{P}_{2}=\frac{{F}_{2}}{\left({z}_{0}{m}_{e}{\omega}_{1}^{2}\right)},{\zeta}_{2}=\frac{{c}_{b}}{\left(2{m}_{g}{\omega}_{1}\right)},{\xi}_{3}=\frac{{c}_{m}}{\left(2{m}_{g}{\omega}_{1}\right)},$
${k}_{2}=\frac{{\omega}_{2}^{2}}{{\omega}_{1}^{2}},{k}_{3}=\frac{{k}_{m}}{\left({m}_{g}{\omega}_{1}^{2}\right)},{P}_{3}=\frac{{F}_{b}}{\left({m}_{g}{z}_{0}{\omega}_{1}^{2}\right),},{P}_{4}=\frac{{F}_{l}\left({a}_{0}+{\mu}_{0}\right)}{\left({m}_{g}{z}_{0}{\omega}_{1}^{2}\right)},$
${\alpha}_{1}=d{\omega}_{1},{\alpha}_{2}=\frac{1}{2}{d}^{2}{\omega}_{1}^{2},{\alpha}_{3}=\frac{1}{6}{d}^{3}{{\omega}_{1}}^{3}$.
Thus, the dimensionless system model can be expressed as:
Letting${x}_{1}={y}_{1}$, ${x}_{2}={\dot{y}}_{1}$, ${x}_{3}={y}_{2}$, ${x}_{4}={\dot{y}}_{2}$, and the state equation of the system can be expressed as:
$=\left(\begin{array}{c}{x}_{2}\\ \left({k}_{1}+{k}_{3}\right){x}_{1}2\left({\zeta}_{1}+{\zeta}_{3}\right){x}_{2}{k}_{2}{x}_{3}\left(2{\zeta}_{2}{\alpha}_{1}{P}_{4}\right){x}_{4}\\ {\alpha}_{2}{P}_{4}{x}_{4}^{2}+{\alpha}_{3}{P}_{4}{x}_{4}^{3}+{P}_{1}+{P}_{2}\mathrm{cos}\left(\mathrm{\Omega}t\right)+{P}_{3}{P}_{4}\\ {x}_{4}\\ {k}_{3}{x}_{1}2{\zeta}_{3}{x}_{2}{k}_{2}{x}_{3}\left(2{\zeta}_{2}{\alpha}_{1}{P}_{4}\right){x}_{4}{\alpha}_{2}{P}_{4}{x}_{4}^{2}+{\alpha}_{3}{P}_{4}{x}_{4}^{3}+{P}_{3}{P}_{4}\end{array}\right),$
where $x=({x}_{1},{x}_{2},{x}_{3},{x}_{4})$.
3.1. LyapunovSchmidt reduction
Ignoring the influence of the external excitation, Eq. (10) becomes autonomous system.
Let:
${l}_{1}={\alpha}_{2}{P}_{4},{l}_{2}={\alpha}_{3}{P}_{4},{l}_{3}={P}_{3}{P}_{4},{a}_{2}={k}_{3},{b}_{2}=2{\zeta}_{3}.$
The equation of motion can be expressed as follows:
Then Eq. (11) can be expressed as:
where $X={\left({x}_{1},{x}_{2},{x}_{3},{x}_{4}\right)}^{T}$ are the state variables of the dynamic system and $f\left(X\right)={\left({f}_{1},{f}_{2},{f}_{3},{f}_{4}\right)}^{T}$ are the corresponding nonlinear vector functions.
At the equilibrium point $\stackrel{}{X}$, Eq. (12) satisfies:
Letting $Df\left(\stackrel{}{X}\right)$ represent the Jacobian matrix of $f\left(X\right)$ at the equilibrium point, the Lyapunov stability is dependent on the eigenvalues of the Jacobian matrix. The Jacobian matrix $Df\left(\stackrel{}{X}\right)$ can be expressed as:
The value of the determinant is:
There is a zero eigenvalue when ${a}_{1}{d}_{1}+2{a}_{1}{l}_{1}{x}_{4}+3{a}_{1}{l}_{2}{x}_{4}^{2}+{a}_{1}{c}_{1}{a}_{2}{c}_{1}=0$, and $L$ is a singular matrix.
It can be known that $\mathrm{d}\mathrm{i}\mathrm{m}N\left(L\right)=\text{1}$. Letting ${v}_{0}$ represent a basis vector of null space $N\left(L\right)$, through calculating, ${v}_{0}=\left({c}_{1}/{a}_{1},\mathrm{0,1},0\right)$. The adjoint operator of $L$ is the conjugate transpose matrix ${L}^{*}$. Since $L$ is a real matrix, the conjugate transpose matrix ${L}^{*}$ is also a real matrix, and satisfies ${L}^{*}={L}^{T}$, we can obtain that $\mathrm{d}\mathrm{i}\mathrm{m}N\left({L}^{*}\right)=\text{1}$ since $L$ is a Fredholm operator with zero index. One basis vector ${v}_{1}=\left({a}_{1}{b}_{2}/{a}_{2}{b}_{1},1,{a}_{1}e/{a}_{2}{a}_{2}e/{a}_{2},{a}_{1}/{a}_{2}\right)$ of $N\left({L}^{*}\right)$ is obtained by calculating, where $e={d}_{1}+2{l}_{1}{x}_{4}+3{l}_{2}{x}_{4}^{2}$. So, the fourdimensional Euclid space ${R}^{4}$ can be decomposed as follows:
Defining the projection operator and the complementary projection operator as:
Then the equation can be expressed as:
Eq. (19) is equivalent to:
where $v\in N\left(L\right)=span\left({v}_{0}\right)$, $w\in M$. According to the definition of inner product in Euclid space, $\left(IP\right)f\left(v+w,{l}_{2}\right)=0$ is equivalent to:
Then, the following equation is obtained:
Since${a}_{1}{d}_{1}+2{a}_{1}{l}_{1}{x}_{4}+3{a}_{1}{l}_{2}{x}_{4}^{2}+{a}_{1}{c}_{1}{a}_{2}{c}_{1}=0\text{,}$ rewriting this equation and substituting it into Eq. (22), we have:
where
$C=\frac{8{a}_{1}^{2}{d}_{1}{l}_{1}}{{a}_{2}\left({a}_{1}{a}_{2}\right)}\frac{4{a}_{1}{l}_{1}}{{a}_{2}}+8{c}_{1}{l}_{1},$
$\begin{array}{l}D=F+{b}_{1}+{b}_{2}+6{l}_{3}{a}_{1}{b}_{1}+{a}_{1}{b}_{2}{a}_{2}{b}_{1}\frac{{a}_{1}}{{a}_{2}}+\frac{{a}_{1}^{2}{b}_{1}}{{a}_{2}}\\ \frac{2{a}_{1}{d}_{1}^{2}}{{a}_{1}}{a}_{2}\frac{2{a}_{1}{d}_{1}}{{a}_{2}}+\frac{2{a}_{1}^{2}{d}_{1}^{2}}{{a}_{2}\left({a}_{1}{a}_{2}\right)}+1.\end{array}$
Eq. (23) is a quartic algebraic equation and any univariate equation with a higher order can eliminate the second highest item by a proper transformation.
3.2. Singularity analysis
Define a new state variable $m={x}_{4}$, then the bifurcation equation can be written as:
where $\alpha $, $\beta $ are the open fold parameters of the bifurcation equation, $m$, $n$ are the physical parameters of the system, respectively.
Selecting ${g}_{0}\left(m,n\right)={m}^{4}+n$ as the GolubitskySchaeffer normal form, according to recognition condition, the necessary and sufficient condition for $G\left(m,n,\alpha ,\beta \right)$ to be a versal unfolding of ${g}_{0}\left(m,n\right)={m}^{4}+n$ is:
where $A$ is, a matrix constructed by some partial derivatives of $g$ and $G$, it can be expressed as:
Thus:
so ${g}_{0}\left(m,n\right)={m}^{4}+n$ is a 2parameter unfolding and $\alpha $, $\beta $ are the unfolding parameter.
The topology of Eq. (24) and the effects of the open fold parameters on the bifurcation behavior will be studied by singularity theory.
The following point sets are obtained based on the definition of the transition set:
Ⅰ: Bifurcation point set:
Ⅱ: Lag point sets:
${H}_{2}=\left\{\left(\alpha ,\beta \right)G={G}_{m}={G}_{mm}=0\right\}=\left\{(\alpha ,\beta )\alpha +\frac{4}{3}\beta \sqrt{\frac{\beta}{6}}=0,\beta \le 0\right\}.$
Ⅲ: Double limits point set:
Ⅳ: Transition set:
The transition set $\mathrm{\Sigma}$ is shown in Fig. 2.
Fig. 2. Transition set
It can be seen that the unfolding parametric space $\alpha \beta $ is divided into four sub regions by $\mathrm{\Sigma}$. The bifurcation diagrams of each subregion and its boundary are shown in Fig. 3, a total of nine parts with Fig. 3(a)(i) including the boundary. Among them, Fig. 3(a)(d) are four different bifurcation regions divided by the transition set and Fig. 3(e)(i) are bifurcation curves in the transition set. The topology of each sub region and its boundary are shown in Fig. 3. It can be seen from Fig. 3 that each sub region and its boundary have different bifurcation forms, and the system vibration mode will change when the unfolding parameters pass by the transition set.
Fig. 3. Bifurcation diagram
a)
b)
c)
d)
e)
f)
g)
h)
i)
4. Hopf bifurcation analysis of high dimensional nonlinear rolling mill system
4.1. The equilibrium point analysis
In last section, the main drive system of rolling mill can be expressed as:
The equilibrium point of the Eq. (27) can be transferred to the coordinate origin by simple linear transformation. Coordinate translation does not change the dynamic characteristics, so it is of significance to study the dynamic characteristics of the system equilibrium point at the origin point. When the equilibrium point is transformed to the coordinate origin, the Eq. (13) can be turned into the following equation:
After linearizing, the Jacobian matrix at the origin can be expressed as
The characteristic equation of the $J({x}_{0},{\alpha}_{2})$ is:
where:
$+{a}_{1}{d}_{1}{a}_{2}{d}_{1}{p}_{4}={a}_{1}{c}_{1}{a}_{2}{c}_{1}.$
It is easy to verify that the characteristic polynomial of Eq. (30) is equivalent to:
where:
Therefore, the Eq. (31) can be written as the following two equations:
Further, letting:
Fixing other system parameters as constants, when the friction coefficient $d$ changes, the characteristic root of Eq. (30) has two important distributions:
Situation 1: when ${\mathrm{\Delta}}_{1}<0$, ${r}_{1}>0$, ${r}_{2}=0$, ${s}_{2}>0$, Eq. (36) has a pair of conjugate complex roots with negative real parts, and Eq. (37) has a pair of pure imaginary roots.
Situation 2: when ${\mathrm{\Delta}}_{2}<0$, ${r}_{2}>0$, ${r}_{1}=0$, ${s}_{1}>0$, Eq. (36) has a pair of pure imaginary roots, and Eq. (37) has a pair of conjugate complex roots with negative real parts.
4.2. Hopf bifurcation analysis
The friction condition has great influence on the dynamic behavior of rolling mill system. So, friction coefficient $d$ is selected as a major parameter to analyze dynamic characteristics. The Situation 1 and Situation 2 all show that characteristic Eq. (30) has a pair of pure imaginary roots and a pair of conjugate complex roots with negative real parts. When the roots satisfy the two situations mentioned above, Hopf bifurcation is very likely to occur.
In order to find the Hopf bifurcation point and prove the existence of limit cycles, in general, we can calculate all the eigenvalues of Jacobian matrix $J({x}_{0},d)\text{,}$ and then judge whether characteristic roots pass through the imaginary axis with the change of parameters. But the amount of calculations is very large and it is difficult to write the analytic expression of characteristic roots. Thus, the Hurwitz determinant is used to determine the existence of the Hopf bifurcation.
The characteristic equation of matrix $J({x}_{0},d)$ is:
where:
where ${\mathrm{\Delta}}_{{}_{i}}^{\text{'}}$ are the Hurwitz determinants. Select $d={d}_{}^{*}$ as a Hopf bifurcation point, the conditions of Hopf bifurcation are:
Take ${k}_{1}=1$, ${k}_{2}=1$, ${k}_{3}=\text{0.6}$, ${\zeta}_{1}=\text{0.38}$, ${\zeta}_{2}=\text{0.38}$, ${\zeta}_{3}=\text{0.09}$, ${P}_{3}=1$, ${P}_{4}=\text{1}$. Hopf bifurcation occurs at ${d}^{*}=\text{1.2324}$.
Supposing that the pure imaginary roots are $\pm {\omega}_{0}i$, $U$ and $W$ are left and right normalized characteristic vectors of ${\omega}_{0}i$, respectively:
Letting:
$+U{f}_{xx}{W}^{*}\left[A\left(0\right)2i{\omega}_{0}I{]}^{1}{f}_{xx}WW\right),$
where ${W}^{*}$ is conjugate complex of $W$, $A\left(0\right)={\left.A\left({x}_{0},d\right)\right}_{d={d}^{*}}$:
The stability of bifurcating periodic solutions is decided by the sign of $\rho $. With ${d}^{*}=\text{1.2324}$, $\rho =\text{0.0900}<\text{0}$ can be obtained by calculation. Bifurcating periodic solutions aren’t asymptotical stable when $\rho <\text{0}$. The results show that it is a subcritical bifurcation point.
4.3. Numerical simulation
Numerical Simulation of the dynamical behavior is made near the bifurcation point. For subcritical bifurcation, when $d=\text{1.2100}<{d}^{*}$. It can be known that the equilibrium point is asymptotically stable from Fig. 4. The unstable periodic motion may occur in the region of $d=\text{1.2280}$.
Fig. 4. Asymptotic stability of the equilibrium point
a) The phase diagram
b) The time history
Fig. 5. Time history responses and phase diagram with $d=\text{1.2280}$ and initial values are a) [0.01 0.01 0.01 0.01], b) [0.12 0.12 0.12 0.12]
a1) The phase diagram
a2) The time history
b1) The phase diagram
b2) The time history
Initial values have a certain influence on the trajectory. When initial point is [0.01 0.01 0.01 0.01], which is close to the equilibrium point, Fig. 5(a1), (a2) show that time history responses are stable and phase diagram is stable focus. When initial point is [0.12 0.12 0.12 0.12], which is far away from the equilibrium point, Fig. 5(b1), (b2) show that time history responses are unstable and phase diagram gradually expands away from a certain limit cycle. Region of attraction on initial point can be obtained by numerical simulation. In the current parameter system, when initial point exceeds [0.09 0.09 0.09 0.09], the time history responses will be unstable.
When $d=\text{1.2432}>{d}^{*}$, the equilibrium point of the system is unstable as shown in Fig. 6, in which case the roller loses its stability and selfoscillation can be generated.
Fig. 6. Asymptotic stability of the equilibrium point
a) The phase diagram
b) The time history
In summary, the system shows a complex dynamic behavior. It can be known that the friction coefficient $d$ and the initial condition of the system are directly related to the stability of the system. When the value of friction coefficient $d$ is much smaller than the bifurcation point, the system is asymptotically stable; when the friction coefficient $d$ near the bifurcation point ${d}^{*}=\text{1.2324}$, the stability of the system is mainly affected by the initial values. When the initial values are close to the equilibrium point, the system is stable; and the stability becomes worse when the initial value is far away from the equilibrium point. When the value of friction coefficient $d$ is larger than the value of bifurcation parameter, the system is divergent, and the phenomenon of the selfexcitation oscillation may occur. In the vicinity of the bifurcation point, the system may have a violent shock and even lead to complex dynamic phenomena. These phenomena should be avoided in the actual production process, so as to ensure the smooth running of the whole rolling mill system.
5. Conclusions
In this paper, a torsionalhorizontal coupling vibration dynamic equation of rolling mill system is established to analyze the bifurcation characteristics. For the fourdimensional dynamic equation, firstly, the dimension is reduced by LyapunovSchmidt reduction method, the topological structure of the system static bifurcation is obtained by singularity analysis, and the parameter condition, which makes the system stable, is given. Then the friction coefficient $d$ is selected as the bifurcation parameter, the dynamic Hopf bifurcation of the system is studied, and the dynamic behavior of the system is verified by numerical simulation. The results show that the system may be unstable when the friction coefficient $d$ changes around the bifurcation point. These phenomena should be avoided in the actual production process to ensure the safe operation of the whole system. It provides a theoretical guide for the study of the horizontalrotational coupling mechanical vibration.
Acknowledgements
Project supported by the National Natural Science Foundation of China (Grant No. 61673334, 51575472). and supported by the National Natural Science Foundation of Hebei Province (Grant No. E2017203144).
References
 Ali Heidari, Mohammad Forouzan R. Optimization of cold rolling process parameters in order to increasing rolling speed limited by chatter vibrations. Journal of Advanced Research, Vol. 4, Issue 1, 2013, p. 2734. [Publisher]
 Liu H., Zhu Z., Shi P. Stability analysis of a relative rotation timedelay nonlinear dynamic system. Acta Physica Sinica, Vol. 59, Issue 10, 2010, p. 67706777. [Search CrossRef]
 Shen Y., Yang S., Li C., et al. Dynamic analysis of a nonlinear system under semiactive control. Journal of Vibration Engineering, Vol. 18, Issue 2, 2005, p. 219222. [Search CrossRef]
 Dong J., Shen Y., Yang S. Main resonance of duffingvan der pol system to multifrequency excitations. Journal of Shijiazhuang Railway Institute, Vol. 19, Issue 1, 2006, p. 6266. [Search CrossRef]
 Li S., Yang S. Subharmonic resonance and singularity of a nonlinear vehicle suspension system. Journal of Vibration Engineering, Vol. 20, Issue 2, 2007, p. 168173. [Search CrossRef]
 Gao G., Yang S., Chen E., et al. One local bifurcation of nonlinear system based on magnetorheological damper. Acta Mechanica Sinica, Vol. 35, Issue 5, 2004, p. 564568. [Search CrossRef]
 Shen Y., Yang S., Liu X. Nonlinear dynamics of a spur gear pair with timevarying stiffness and backlash based on incremental harmonic balance method. International Journal of Mechanical Sciences, Vol. 48, Issue 11, 2006, p. 12561263. [Publisher]
 Shen Y., Yang S., Pan C., et al. Nonlinear dynamics of a spur gear pair with timevarying stiffness and backlash. Journal of Low Frequency Noise Vibration and Active Control, Vol. 23, Issue 3, 2004, p. 179187. [Search CrossRef]
 Li X., Hou J. Bursting phenomenon in a piecewise mechanical system with parameter perturbation in stiffness. International Journal of NonLinear Mechanics, Vol. 81, 2016, p. 165176. [Publisher]
 Shi P., Ji Z., Liu B., et al. Nonlinear dynamics of torsional vibration for rolling mill’s main drive system under parametric excitation. Journal of Iron and Steel Research International, Vol. 20, Issue 1, 2013, p. 712. [Publisher]
 Liu G., Di H., Zhou C., et al. Tension and thickness control strategy analysis of two stands reversible cold rolling mill. Journal of Iron and Steel Research International, Vol. 19, Issue 10, 2012, p. 2025. [Publisher]
 Liu S., Liu B., Ye K., et al. Hopf bifurcation and stability of periodic solutions in a nonlinear relative rotation dynamical system with time delay. Acta Physica Sinica, Vol. 59, Issue 1, 2010, p. 3843. [Search CrossRef]
 Yang X., Li Q., Tong C., et al. Vertical vibration model for unsteady lubrication in rollsstrip interface of cold rolling mills. Advances in Mechanical Engineering, Vol. 4, Issue 12, 2012, p. 14321437. [Publisher]
 Hou D., Peng R., Liu H. Analysis of verticalhorizontal coupling vibration characteristics of rolling mill rolls based on strip dynamic deformation process. Shock and Vibration, Vol. 2014, Issue 2, 2014, p. 111. [Publisher]
 Hou D., Zhu Y., Liu H., et al. Research on nonlinear vibration characteristics of cold rolling mill based on dynamic rolling force. Journal of Mechanical Engineering, Vol. 49, Issue 14, 2013, p. 4550. [Publisher]
 Xu Y., Tong C., Yue G., et al. Coupling Dynamic model of chatter for cold rolling. Journal of Iron and Steel Research International, Vol. 17, Issue 12, 2010, p. 3034. [Search CrossRef]
 Zhang R., Chen Z., Wang F. Study on parametrically excited horizontal nonlinear vibration in singleroll driving mill system. Journal of Vibration and Shock, Vol. 29, Issue 6, 2010, p. 105108. [Search CrossRef]
 Sun J., Yan P., Liu H. Nonlinear vibration and stability of moving strip with timedependent tension in rolling process. Journal of Iron and Steel Research International, Vol. 17, Issue 6, 2010, p. 1115. [Publisher]
 Zeng L., Zang Y., Gao Z. Effect of rolling process parameters on stability of rolling mill vibration with nonlinear friction. Journal of Vibroengineering, Vol. 18 Issue 2, 2016, p. 12881296. [Search CrossRef]
 Liu H., Zhang J., Mi K., et al. Simulation on hydraulicmechanical coupling vibration of cold strip rolling mill vertical system. Advanced Materials Research, Vol. 694697, 2013, p. 407414. [Publisher]
 Senjanović I., Ćatipović I., Tomašević Coupled horizontal and torsional vibrations of a flexible barge. Engineering Structures, Vol. 30, Issue 1, 2008, p. 93109. [Publisher]
Articles Citing this One
Journal of Iron and Steel Research International
Dongxiao Hou, Liang Xu, Peiming Shi

2021

Mechanism and Machine Theory
Tiancheng Ouyang, Zixiang Su, Shuoyu Li, Guicong Huang, Nan Chen

2019
