Nonlinear vibration characteristic of strip mill under the coupling effect of roll-rolled piece

Liu, Bin; Jiang, Jiahao; Liu, Fei; Liu, Haoran; Li, Peng

doi:10.21595/jve.2016.17056

Journal of Vibroengineering

Browse Journal

Submit article

Published: 31 December 2016

Check for updates

Nonlinear vibration characteristic of strip mill under the coupling effect of roll-rolled piece

Bin Liu¹

Jiahao Jiang²

Fei Liu³

Haoran Liu⁴

Peng Li⁵

^{1, 2, 5}School of Electrical Engineering, Yanshan University, Qinhuangdao, China

^{3, 4}School of Informational Science and Engineering, Yanshan University, Qinhuangdao, China

Corresponding Author:

Haoran Liu

Cite the article Download PDF

Downloads 1445

WoS Core Citations 7

CrossRef Citations 7

Abstract

In the existing friction coefficient models, the influence of relatively horizontal vibration of rolled piece is always overlooked, this makes the friction coefficient model error and ultimately affects the accuracy of vibration model of rolling mill. In order to be more close to the actual situation, the vertical vibration model of strip mill is established by considering the effect of horizontal vibration of rolled piece. The approximate analytic solutions are obtained by multiple scales method. Then, bifurcation characteristics of system are analyzed by means of Singularity theory. Finally, frequency responses of system are investigated with the change of external excitation amplitudes and internal nonlinear parameters. Results show that vibration behavior of system changes with the external excitation amplitudes and frequencies, internal nonlinear parameters and coupling parameters. Moreover, the unstable frequency region of system can be obtained by amplitude-frequency curve.

1. Introduction

Plate and strip is the most popular rolling product, widely used in automobile manufacturing, shipping, aerospace and military field processing [1]. With the development of society, the size precision and surface quality of the plate and strip has been put forward higher requirements [2]. However, the vertical vibration of rolling mill has been a problem for enhancing the efficiency of cold rolling products all over the world. The vertical vibration of rolling mill will produce vibration mark in strip and roll surface, so that the quality of the products will be greatly discounted. Furthermore, the severe vibration of rolling mill even cause broken belt and equipment damage accidents, resulting in serious economic losses [3-5].

For the vertical vibration problem of strip mill in the rolling process, experts and scholars have established a variety vertical vibration models of rolling mill from different angles, and some suggests are given to restrain the vertical vibration of rolling mill. Yarita et al. presented a linear spring-dumping model to represent the interaction between roll system and upper beam of frame. He built a 4-DOF vibration model to study the vertical vibration of rolling mill [6]. Tamiya et al. considered that the mass of work roll is far more less than the backup roll, so he omitted the mass of work roll and built a 2-DOF vibration model based on the research of Yarita [7]; Soon after, Roberts et al. assumed that the rolling mill structure is equivalent along the rolling line, so he presented a single freedom vibration model of rolling mill [8]. With the development of nonlinear science, scholars begin to study the vertical vibration characteristics of rolling mill under the influence of nonlinear stiffness and nonlinear damping. Liu et al. set up a piecewise nonlinear vibration model of the hydraulic cylinder, and the vibration behavior of rolling mill with the nonlinear constraint of hydraulic cylinder were studied [9].

In recent years, with the further research on the vibration of rolling mill, some of vibration phenomena of the mill cannot be explained by simple rolling process theory. Scholars begin to study the mill vibration from coupling angles. Yun et al. formulated a 2-DOF coupling model by coupling the horizontal vibration and vertical vibration of rolls, then the expression of dynamic component of the rolling force was derived, and they applied it in the mill structure [10]. Zhang et al. researched on the rolling mill vibration caused by flexural-vibration of the strip, and presented an electromechanical coupling vibration model of rolling mill; their work explained the unsteady vibration phenomenon, and attached importance to the parametric vibration of the rolling mill [11]. Yang et al. proposed a coupling vibration model by coupling the rolling process model, the mill roll stand structure model and the hydraulic servo system model; meanwhile the effects of different working conditions on the stability of cold rolling mill system were discussed; It contributes to the further study and suppression of coupling vibration [12].

Although the numerous researches on strip mill vibration provide theoretical support and reference for the production and the design of rolling mill, the unsteady vibration of strip mill caused by horizontal vibration of rolled piece is always overlooked. However, in the actual rolling process of strip mill, there are two kinds of friction states between rolled piece and roll, static friction and sliding friction. These two different friction sates will bring about different effects on the vibration of the rolling mill system. The vibration of rolled piece will affect the rolling mill system by changing the roll gap friction coefficient. Furthermore, the coupling vibration of rolled piece and roll will be a direct threat to the quality of strip. While in the state of static friction, the effect is weaker. So the vibration of strip mill caused by horizontal vibration of rolled piece must be taken seriously and its further study is urgent [13].

The purpose of this paper is to analyze the dynamic behavior and vertical vibration characteristic of strip mill under the coupling effect of roll-rolled piece, especially considering the influence of horizontal vibration of rolled piece on the friction coefficient, which is always neglected in former researches. In order to couple the rolling mill vertical vibration and the horizontal vibration of rolled piece, a coupling vibration model of rolling mill is proposed. Then by using multiple scales method, the analytic solutions are obtained. Finally, the static bifurcation characteristics and amplitude-frequency characteristics of coupling system are analyzed. The research results may provide a theoretical reference for the vibration suppression of the rolling mill.

2. Mathematical modeling of the strip mill structure

2.1. Modeling of friction force and rolling force

Due to the influence factors of rolling force varied, there are various forms of rolling force formula. All in that formulas, the Hill formula is most widely used and has good accuracy [14]. Therefore, the rolling force can be expressed with Hill formula as follows:

1

F = B l_{c} Q_{p} K_{T} K,

where:

l_{c} = \sqrt{R Δ h}, Q_{p} = 1.08 + 1.79 μ ε \sqrt{1 - ε} \sqrt{\frac{R}{H}} - 1.02 ε,

K_{T} = 1 - \frac{0.7 τ_{b} + 0.3 τ_{f}}{K}, K = 1.15 a_{0} {(\frac{1}{3} ε_{0} + \frac{2}{3} ε_{1})}^{a_{1}},

where, $B$ is the width of rolled piece; $l_{c}$ is the contact length between roll and rolled piece in deformation zone; $Q_{p}$ is the influential coefficient in stressed state; $K_{T}$ is the tension coefficient; $K$ is the average deformation resistance of materials; $a_{0}$ , $a_{1}$ is the regression coefficient of the model; $R$ is the work roll radius; $Δ h$ is reduction quantity of rolled piece, $Δ h = H - h - 2 y$ , $H$ is the entrance thickness of rolled piece, $h$ is the exit thickness of rolled piece, $y$ is the vibration displacement of rolls; $ε$ is the reduction rate of rolled piece, $ε = Δ h / H$ ; $ε_{0}$ is the reduction rate of frame ingress, $ε_{0} = (\bar{H} - H) / \bar{H}$ ; $ε_{1}$ is the reduction rate of frame egress, $ε_{1} = (\bar{H} - h) / \bar{H}$ , $\bar{H}$ is the thickness of pre-rolling rolled piece; $τ_{f}$ and $τ_{b}$ respectively represent the forward and backward tensile stress of rolled piece; $μ$ is the friction coefficient of roll gap, it may be expressed as follows [15]:

2

μ = \sqrt{\frac{Δ h}{D}} [0.5 + (K_{1} - 0.5) e^{- K_{2} (v_{0} + \dot{x})}],

where, $K_{1}$ and $K_{2}$ are the friction characteristic coefficients; $D$ is the work roll diameter; $v_{0}$ is work roll rotational speed; $\dot{x}$ is the relative vibration speed of rolled piece in roll gap; considering that $K_{2} \dot{x} ≪ 1$ , so the Roberts formula can be simplified as:

3

μ = μ (0,0) + Δ μ,

where:

\{\begin{array}{l} μ (0,0) = \sqrt{\frac{H - h}{D}} [0.5 + (K_{1} - 0.5) e^{- K_{2} v_{0}}], \\ Δ μ = b_{1} \dot{x} + b_{2} y + b_{3} {\dot{x}}^{2} + b_{4} \dot{x} y + b_{5} y^{2} + b_{6} {\dot{x}}^{3} + b_{7} {\dot{x}}^{2} y + b_{8} \dot{x} y^{2} + b_{9} y^{3}, \end{array}

where:

b_{0} = μ (0,0), b_{1} = \frac{\partial}{\partial \dot{x}} μ (0,0), b_{2} = \frac{\partial}{\partial y} μ (0,0), b_{3} = \frac{1}{2!} \frac{\partial^{2}}{\partial {\dot{x}}^{2}} μ (0,0),

b_{4} = \frac{\partial^{2}}{\partial \dot{x} \partial y} μ (0,0), b_{5} = \frac{1}{2!} \frac{\partial^{2}}{\partial y^{2}} μ (0,0), b_{6} = \frac{1}{3!} \frac{\partial^{3}}{\partial {\dot{x}}^{3}} μ (0,0),

b_{7} = \frac{1}{2} \frac{\partial^{3}}{\partial {\dot{x}}^{2} \partial y} μ (0,0), b_{8} = \frac{1}{2} \frac{\partial^{3}}{\partial y^{2} \partial \dot{x}} μ (0,0), b_{9} = \frac{1}{3!} \frac{\partial^{3}}{\partial y^{3}} μ (0,0) .

When rolling process is stable, $\dot{x} = 0$ , $y = 0$ . Substituting Eq. (2) into Eq. (1), then expanding Eq. (1) at equilibrium point by Taylor formula, the rolling force can be represented as:

4

F = F (0,0) + Δ F,

where:

\{\begin{array}{l} F (0,0) = 1.15 a_{0} B L_{1} L_{2} \sqrt{R H - h}, \\ Δ F = I_{1} \dot{x} + I_{2} y + I_{3} {\dot{x}}^{2} + I_{4} \dot{x} y + I_{5} y^{2} + I_{6} {\dot{x}}^{3} + I_{7} {\dot{x}}^{2} y + I_{8} \dot{x} y^{2} + I_{9} y^{3}, \end{array}

where:

L_{1} = {(\frac{1}{3} ε_{0} + \frac{2}{3} ε_{1})}^{a_{1}} - 0.7 τ_{b} - 0.3 τ_{f}, L_{2} = 1.08 + \frac{1.79 b_{0} e^{- K_{2} v_{0}} \sqrt{h R} - 1.02 H}{H^{2}} (H - h),

I_{0} = F (0,0), I_{1} = \frac{\partial}{\partial \dot{x}} F (0,0), I_{2} = \frac{\partial}{\partial y} F (0,0), I_{3} = \frac{1}{2!} \frac{\partial^{2}}{\partial {\dot{x}}^{2}} F (0,0),

I_{4} = \frac{\partial^{2}}{\partial \dot{x} \partial y} F (0,0), I_{5} = \frac{1}{2!} \frac{\partial^{2}}{\partial y^{2}} F (0,0), I_{6} = \frac{1}{3!} \frac{\partial^{3}}{\partial {\dot{x}}^{3}} F (0,0),

I_{7} = \frac{1}{2} \frac{\partial^{3}}{\partial {\dot{x}}^{2} \partial y} F (0,0), I_{8} = \frac{1}{2} \frac{\partial^{3}}{\partial y^{2} \partial \dot{x}} F (0,0), I_{9} = \frac{1}{3!} \frac{\partial^{3}}{\partial y^{3}} F (0,0) .

According to the Friction theory, the friction formula can be written as:

5

F_{f} = 2 (μ_{0} + Δ μ) \times (F (0,0) + Δ F),

where, $μ_{0}$ and $F (0,0)$ are the steady-state values, $Δ μ$ and $Δ F$ are the dynamic variable quantities; Take the main part of the friction force, Eq. (5) can be written as:

6

Δ F_{f} = 2 μ_{0} Δ F + 2 Δ μ F_{0} .

By substituting Eq. (3) and Eq. (4) into Eq. (6), the expression of $Δ F_{f}$ can be obtained:

7

Δ F_{f} = 2 b_{0} [I_{1} \dot{x} + I_{2} y + I_{3} {\dot{x}}^{2} + I_{4} \dot{x} y + I_{5} y^{2} + I_{6} {\dot{x}}^{3} + I_{7} {\dot{x}}^{2} y + I_{8} \dot{x} y^{2} + I_{9} y^{3}]

+ 2 I_{0} [b_{1} \dot{x} + b_{2} y + b_{3} {\dot{x}}^{2} + b_{4} \dot{x} y + b_{5} y^{2} + b_{6} {\dot{x}}^{3} + b_{7} {\dot{x}}^{2} y + b_{8} \dot{x} y^{2} + b_{9} y^{3}] .

2.2. Dynamic model of strip mill mechanical structure

The rolling system structure is very complex, in order to facilitate the research, scholars often simplify the rolling system by using the lumped mass method, the method is simple and effective, it has been recognized by industry experts [16]. Because the mass of work roll is far more less than the backup roll, So, according to the method, the work roll and backup roll can be viewed as a mass block. Based on the mass-spring-damping classic model, a structural model of plate and strip rolling mill is established as shown in Fig. 1. The model especially consider the effect of horizontal vibration of rolled piece.

Fig. 1Structure model of strip mill

In Fig. 1, $m$ is the mass of rolled piece in roll gap, it can be obtained by the equation $m = ρ V$ ; where $ρ$ is the rolled piece density, $V$ is the volume of rolled piece in roll gap, and it may be approximately expressed as this: $V = B l_{c} (H + h) / 2$ ; The wave force of rolled piece in the forward and backward sliding zones is equivalent to the force function. Its stiffness can be written as: $k_{f} = {τ_{f} S_{f} / L}_{f d}$ , $k_{b} = {τ_{b} S_{b} / L}_{b d}$ ; where $τ_{f}$ and $τ_{b}$ stand for forward and backward tension of rolled piece, $S_{f}$ and $S_{b}$ stand for forward and backward cross-sectional area of rolled piece, $L_{f d}$ and $L_{b d}$ stand for forward and backward deformation length of rolled piece caused by external tension. The damping effect between the roll and rolled piece is equal to two equivalent damping $c$ . $m_{1}$ and $m_{2}$ are respectively the equivalent mass of upper and lower rolls; $k_{1}$ and $k_{2}$ are respectively the equivalent stiffness of upper and lower rolls; $c_{1}$ and $c_{2}$ are respectively the equivalent dumping of upper and lower rolls; $y_{1}$ and $y_{2}$ are the vibration displacement of upper and lower rolls; $R$ is the roll radium; $F^{*}$ is the external disturbance force of rolls.

According to the Lagrange principle of generalized dissipation, the dynamic equilibrium equation of the roll in the vertical direction can be expressed as:

8

\{\begin{array}{l} m_{1} {\ddot{y}}_{1} + c_{1} {\dot{y}}_{1} + k_{1} y_{1} = Δ F + F^{*}, \\ m_{2} {\ddot{y}}_{2} + c_{2} {\dot{y}}_{2} + k_{2} y_{2} = - Δ F - F^{*} . \end{array}

Based on the assumption that the structure and vibration characteristic of rolling mill is symmetrical in relation to rolling line, there exist $m_{1} = m_{2}$ ; $c_{1} = c_{2}$ ; $y_{1} = - y_{2}$ ; $k_{1} = k_{2}$ . So, the two equations in Eq. (8) have same expression form. In order to simplify the analysis procedure, only research on the upper rolls. The equivalent mass of upper rolls is denoted by $M$ , the equivalent stiffness of upper rolls is denoted by $K$ , the equivalent dumping of upper rolls is denoted by $C$ . Based on the simplification and equivalence, Eq. (8) can be written as:

9

M \ddot{y} + C \dot{y} + K y = Δ F + F^{*} .

Based on the truth that in cold rolling process the roll is close to completely elastic flattening, namely, the contact surface between roll and rolled piece is close to a plane. On this basis, horizontal vibration dynamic equation can be expressed as:

10

m \ddot{x} + 2 c \dot{x} + (k_{f} + k_{b}) x = Δ F_{f} .

By combining Eqs. (4, 7, 9, 10), coupling vibration equation of roll-rolled piece is obtained:

11

\{\begin{array}{l} m \ddot{x} + 2 c \dot{x} + (k_{f} + k_{b}) x = 2 b_{0} [I_{1} \dot{x} + I_{2} y + I_{3} {\dot{x}}^{2} + I_{4} \dot{x} y + I_{5} y^{2} + I_{6} {\dot{x}}^{3} + I_{7} {\dot{x}}^{2} y + I_{8} \dot{x} y^{2} \\ + I_{9} y^{3}] + 2 I_{0} [b_{1} \dot{x} + b_{2} y + b_{3} {\dot{x}}^{2} + b_{4} \dot{x} y + b_{5} y^{2} + b_{6} {\dot{x}}^{3} + b_{7} {\dot{x}}^{2} y + b_{8} \dot{x} y^{2} + b_{9} y^{3}], \\ M \ddot{y} + C \dot{y} + K y = I_{1} \dot{x} + I_{2} y + I_{3} {\dot{x}}^{2} + I_{4} \dot{x} y + I_{5} y^{2} + I_{6} {\dot{x}}^{3} + I_{7} {\dot{x}}^{2} y + I_{8} \dot{x} y^{2} + I_{9} y^{3} + F^{*} . \end{array}

3. Coupling system solution of strip mill

Assuming that the system is subjected to periodic external disturbances, set $F^{*} = F_{0} c o s Ω t$ . By transposition and replacement, Eq. (11) is transformed into a standard form:

12

\{\begin{array}{l} \ddot{x} + {ω_{10}}^{2} x = η_{11} f_{1} (\dot{x}, y) + η_{12} f_{2} (\dot{x}, y) - β_{1} \dot{x}, \\ \ddot{y} + {ω_{20}}^{2} y = η_{21} f_{1} (\dot{x}, y) + η_{22} c o s Ω t - β_{2} \dot{y}, \end{array}

where:

f_{1} (\dot{x}, y) = I_{1} \dot{x} + I_{2} y + I_{3} {\dot{x}}^{2} + I_{4} \dot{x} y + I_{5} y^{2} + I_{6} {\dot{x}}^{3} + I_{7} {\dot{x}}^{2} y + I_{8} \dot{x} y^{2} + I_{9} y^{3},

f_{2} (\dot{x}, y) = b_{1} \dot{x} + b_{2} y + b_{3} {\dot{x}}^{2} + b_{4} \dot{x} y + b_{5} y^{2} + b_{6} {\dot{x}}^{3} + b_{7} {\dot{x}}^{2} y + b_{8} \dot{x} y^{2} + b_{9} y^{3},

where:

β_{1} = \frac{2 c}{m}, ω_{10} = \sqrt{\frac{(k_{f} + k_{b})}{m}}, η_{11} = \frac{2 b_{0}}{m}, η_{12} = \frac{2 I_{0}}{m},

β_{2} = \frac{C}{M}, ω_{20} = \sqrt{\frac{K}{M}}, η_{21} = \frac{1}{M}, η_{22} = \frac{1}{M} .

Set, $β_{1} = ε {\hat{β}}_{1}$ ; $β_{2} = ε {\hat{β}}_{2}$ ; $η_{11} = ε {\hat{η}}_{11}$ ; $η_{12} = ε {\hat{η}}_{12}$ ; $η_{13} = ε {\hat{η}}_{13}$ ; $η_{21} = ε {\hat{η}}_{21}$ ; $η_{22} = ε {\hat{η}}_{22}$ .

By parameters replacement, Eq. (12) becomes:

13

\{\begin{array}{l} \ddot{x} + {ω_{10}}^{2} x = ε ({\hat{η}}_{11} f_{1} (\dot{x}, y) + {\hat{η}}_{12} f_{2} (\dot{x}, y) - {\hat{β}}_{1} \dot{x}), \\ \ddot{y} + {ω_{20}}^{2} y = ε ({\hat{η}}_{21} f_{1} (\dot{x}, y) + {\hat{η}}_{22} F_{0} c o s Ω t - {\hat{β}}_{2} \dot{y}) . \end{array}

For finding a nonlinear approximate solution of Eq. (13), two time scales of $T_{0} = t$ and $T_{1} = ε t$ are selected. Then, the time derivatives are defined as:

14

\begin{array}{l} \frac{d}{d t} = D_{0} + ε D_{1} + \dots, \\ \frac{d^{2}}{d t^{2}} = {D_{0}}^{2} + 2 ε D_{0} D_{1} + \dots, \end{array}

where, $D_{n} = \partial / \partial T_{n}$ , $ε$ is small parameter. By using multiple scales method, set the solution of Eq. (14) as:

15

\{\begin{array}{l} x = x_{0} (T_{0}, T_{1}) + ε x_{1} (T_{0}, T_{1}), \\ y = y_{0} (T_{0}, T_{1}) + ε y_{1} (T_{0}, T_{1}) . \end{array}

Substituting Eq. (14) and (15) into (13) and separating terms of each order of, one has:

16

\{\begin{array}{l} D_{0}^{2} x_{0} + ω_{10}^{2} x_{0} = 0, \\ D_{0}^{2} y_{0} + ω_{20}^{2} y_{0} = 0, \end{array}

17

\{\begin{array}{l} D_{0}^{2} x_{1} + ω_{10}^{2} x_{1} = {\hat{η}}_{11} [\begin{array}{l} I_{1} D_{0} x_{0} + I_{2} y_{0} + I_{3} (D_{0} x_{0})^{2} + I_{4} D_{0} x_{0} y_{0} + I_{5} {y_{0}}^{2} \\ + I_{6} (D_{0} x_{0})^{3} + I_{7} (D_{0} x_{0})^{2} y_{0} + I_{8} D_{0} x_{0} {y_{0}}^{2} + I_{9} {y_{0}}^{3} \end{array}] \\ + {\hat{η}}_{12} [\begin{array}{l} b_{1} D_{0} x_{0} + b_{2} y_{0} + b_{3} (D_{0} x_{0})^{2} + b_{4} D_{0} x_{0} y_{0} + b_{5} y_{0}^{2} \\ + b_{6} (D_{0} x_{0})^{3} + b_{7} (D_{0} x_{0})^{2} y_{0} + b_{8} D_{0} x_{0} {y_{0}}^{2} + b_{9} y_{0}^{3} \end{array}] - {\hat{β}}_{1} D_{0} x_{0} - 2 D_{0} D_{1} x_{0}, \\ {D_{0}}^{2} y_{1} + {ω_{20}}^{2} y_{1} = {\hat{η}}_{21} [\begin{array}{l} I_{1} D_{0} x_{0} + I_{2} y_{0} + I_{3} (D_{0} x_{0})^{2} + I_{4} D_{0} x_{0} y_{0} + I_{5} {y_{0}}^{2} \\ + I_{6} (D_{0} x_{0})^{3} + I_{7} (D_{0} x_{0})^{2} y_{0} + I_{8} D_{0} x_{0} {y_{0}}^{2} + I_{9} {y_{0}}^{3} \end{array}] \\ - {\hat{β}}_{2} D_{0} y_{0} + {\hat{η}}_{22} F_{0} c o s Ω T_{0} - 2 D_{0} D_{1} y_{0} . \end{array}

The solution of Eq. (16) is setting as:

18

\{\begin{array}{l} x_{0} = A (T_{1}) e^{i ω_{10} T_{0}} + \bar{A} (T_{1}) e^{- i ω_{10} T_{0}}, \\ y_{0} = B (T_{1}) e^{i ω_{20} T_{0}} + \bar{B} (T_{1}) e^{- i ω_{20} T_{0}} . \end{array}

Substituting Eqs. (18) into (13), and system internal resonance is taken into account. By using the small scale detuning parameters, the frequencies are redefined as: $Ω = ω {}_{20} + ε σ;$ $ω_{10} = ω_{20} + ε σ_{1}$ . Where, $σ$ and $σ_{1}$ are detuning parameters. In order to avoid the secular terms of equations, Eq. (17) must satisfy the conditions as follow:

19

\{\begin{array}{l} {\hat{η}}_{11} (\begin{array}{l} I_{1} A i ω_{10} + I_{2} B e^{- i σ_{1} T_{1}} + 3 I_{6} A^{2} \bar{A} i ω_{10}^{3} \\ + 2 I_{7} A \bar{A} B ω_{10}^{2} e^{- i σ_{1} T_{1}} + 2 I_{8} A B \bar{B} i ω_{10} + 3 I_{9} B^{2} \bar{B} e^{- i σ_{1} T_{1}} \end{array}) \\ {+ \hat{η}}_{12} (\begin{array}{l} b_{1} A i ω_{10} + b_{2} B e^{- i σ_{1} T_{1}} + 3 b_{6} A^{2} \bar{A} i ω_{10}^{3} \\ + 2 b_{7} A \bar{A} B ω_{10}^{2} e^{- i σ_{1} T_{1}} + 2 b_{8} A B \bar{B} i ω_{10} + 3 b_{9} B^{2} \bar{B} e^{- i σ_{1} T_{1}} \end{array}) \\ - {\hat{β}}_{1} A i ω_{10} - 2 D_{1} A i ω_{10} = 0, \\ {\hat{η}}_{21} (\begin{array}{l} I_{1} A i ω_{10} e^{i σ_{1} T_{1}} + I_{2} B + 3 I_{6} A^{2} \bar{A} i ω_{10}^{3} e^{i σ_{1} T_{1}} \\ + 2 I_{7} A \bar{A} B ω_{10}^{2} + 2 I_{8} A B \bar{B} i ω_{10} e^{i σ_{1} T_{1}} + 3 I_{9} B^{2} \bar{B} \end{array}) \\ - \hat{β} {}_{2} B i ω_{20} + \frac{{\hat{η}}_{22} F_{0}}{2} e^{i σ_{1} T_{1}} - 2 D_{1} B i ω_{20} = 0 . \end{array}

To solve Eq. (19), it is convenient to express the solution in polar form:

A = \frac{1}{2} a e^{i θ_{1}}, B = \frac{1}{2} b e^{i θ_{2}},

where, $a$ , $b$ , $θ_{1}$ , $θ_{2}$ both are the functions of $T_{1}$ . In order to obtain the solution of equation set, introducing intermediate variables $ϕ$ , $φ$ ; define that: $ϕ = θ_{2} - θ_{1} - σ_{1} T_{1}$ ; $φ = σ T_{1} - θ_{2}$ ; Substituting $A$ , $B$ , $ϕ$ , $φ$ into Eq. (19), the modulation equations are expressed as:

20

\{\begin{array}{l} \dot{a} = \frac{{\hat{η}}_{11}}{2} (I_{1} a + \frac{3}{4} I_{6} a^{3} ω_{10}^{2} + \frac{1}{2} I_{8} a b^{2}) + \frac{{\hat{η}}_{12}}{2} (b_{1} a + \frac{3}{4} b_{6} a^{3} ω_{10}^{2} + \frac{1}{2} b_{8} a b^{2}) \\ - \frac{1}{2} {\hat{β}}_{1} a + [\begin{array}{l} \frac{{\hat{η}}_{11}}{ω_{10}} (\frac{1}{2} I_{2} b + \frac{1}{4} I_{7} a^{2} b ω_{10}^{2} + \frac{3}{8} I_{9} b^{3}) \\ + \frac{{\hat{η}}_{12}}{ω_{10}} (\frac{1}{2} b_{2} b + \frac{1}{4} b_{7} a^{2} b ω_{10}^{2} + \frac{3}{8} b_{9} b^{3}) \end{array}] s i n φ, \\ a (\dot{φ} + \dot{ϕ}) = [\begin{array}{l} \frac{{\hat{η}}_{11}}{ω_{10}} (\frac{1}{2} I_{2} b + \frac{1}{4} I_{7} a^{2} b ω_{10}^{2} + \frac{3}{8} I_{9} b^{3}) \\ + \frac{{\hat{η}}_{12}}{ω_{10}} (\frac{1}{2} b_{2} b + \frac{1}{4} b_{7} a^{2} b ω_{10}^{2} + \frac{3}{8} b_{9} b^{3}) \end{array}] c o s φ + a (σ - σ_{1}), \\ \dot{b} = \frac{{\hat{η}}_{21}}{ω_{20}} (\frac{1}{2} I_{1} a ω_{10} + \frac{3}{8} I_{6} a^{3} ω_{10}^{3} + \frac{1}{4} I_{8} a b^{2} ω_{10}) c o s φ - \frac{1}{2} \hat{β} {}_{2} b + \frac{{\hat{η}}_{22} F_{0}}{2 ω_{20}} s i n ϕ, \\ b \dot{ϕ} = \frac{{\hat{η}}_{21}}{ω_{20}} (\begin{array}{l} (\frac{1}{2} I_{1} a ω_{10} + \frac{3}{8} I_{6} a^{3} ω_{10}^{3} + \frac{1}{4} I_{8} a b^{2} ω_{10}) s i n φ \\ + \frac{1}{2} I_{2} b + \frac{1}{4} I_{7} a^{2} b ω_{10}^{2} + \frac{3}{8} I_{9} b^{3} \end{array}) + \frac{{\hat{η}}_{22} F_{0}}{2 ω_{20}} c o s ϕ + b σ . \end{array}

Eliminating $ϕ$ , $φ$ , the frequency response is obtained in terms of two coupled equations as:

21

{[\frac{{\hat{η}}_{11}}{2} (I_{1} a + \frac{3}{4} I_{6} a^{3} ω_{10}^{2} + \frac{1}{2} I_{8} a b^{2}) + \frac{{\hat{η}}_{12}}{2} (b_{1} a + \frac{3}{4} b_{6} a^{3} ω_{10}^{2} + \frac{1}{2} b_{8} a b^{2}) - \frac{1}{2} {\hat{β}}_{1} a]}^{2}

+ {[a (σ - σ_{1})]}^{2}

= {[\frac{{\hat{η}}_{11}}{ω_{10}} (\frac{1}{2} I_{2} b + \frac{1}{4} I_{7} a^{2} b ω_{10}^{2} + \frac{3}{8} I_{9} b^{3}) + \frac{{\hat{η}}_{12}}{ω_{10}} (\frac{1}{2} b_{2} b + \frac{1}{4} b_{7} a^{2} b ω_{10}^{2} + \frac{3}{8} b_{9} b^{3})]}^{2},

22

{[\frac{{\hat{η}}_{21}}{ω_{20}} (\frac{1}{2} I_{1} a ω_{10} + \frac{3}{8} I_{6} a^{3} ω_{10}^{3} + \frac{1}{4} I_{8} a b^{2} ω_{10}) c o s ϕ - \frac{1}{2} {\hat{β}}_{2} b]}^{2}

{+ [\frac{{\hat{η}}_{21}}{ω_{20}} (\begin{array}{l} (\frac{1}{2} I_{1} a ω_{10} + \frac{3}{8} I_{6} a^{3} ω_{10}^{3} + \frac{1}{4} I_{8} a b^{2} ω_{10}) s i n ϕ \\ + \frac{1}{2} I_{2} b + \frac{1}{4} I_{7} a^{2} b ω_{10}^{2} + \frac{3}{8} I_{9} b^{3} \end{array}) + b σ]}^{2} = {(\frac{{\hat{η}}_{22} F_{0}}{2 ω_{20}})}^{2},

where:

c o s ϕ = \frac{- a (σ - σ_{1})}{\frac{{\hat{η}}_{11}}{ω_{10}} (\frac{1}{2} I_{2} b + \frac{1}{4} I_{7} a^{2} b {ω_{10}}^{2} + \frac{3}{8} I_{9} b^{3}) + \frac{{\hat{η}}_{12}}{ω_{10}} (\frac{1}{2} b_{2} b + \frac{1}{4} b_{7} a^{2} b {ω_{10}}^{2} + \frac{3}{8} b_{9} b^{3})},

s i n ϕ = - \frac{\frac{{\hat{η}}_{11}}{2} (I_{1} a + \frac{3}{4} I_{6} a^{3} {ω_{10}}^{2} + \frac{1}{2} I_{8} a b^{2}) + \frac{{\hat{η}}_{12}}{2} (b_{1} a + \frac{3}{4} b_{6} a^{3} {ω_{10}}^{2} + \frac{1}{2} b_{8} a b^{2}) - \frac{1}{2} {\hat{β}}_{1} a}{\frac{{\hat{η}}_{11}}{ω_{10}} (\frac{1}{2} I_{2} b + \frac{1}{4} I_{7} a^{2} b {ω_{10}}^{2} + \frac{3}{8} I_{9} b^{3}) + \frac{{\hat{η}}_{12}}{ω_{10}} (\frac{1}{2} b_{2} b + \frac{1}{4} b_{7} a^{2} b {ω_{10}}^{2} + \frac{3}{8} b_{9} b^{3})} .

4. Bifurcation characteristic analysis of strip mill

Substituting Eq. (21) into (22), and eliminating $a$ . Then the bifurcation equation is obtained:

23

z^{6} + α_{1} z^{4} + α_{2} z^{3} + α_{3} z^{2} + α_{4} z + λ = 0,

where, $α_{1}$ , $α_{2}$ , $α_{3}$ , $α_{4}$ are the unfolding parameters related to the internal parameters of system; $λ$ is the bifurcation parameter mainly related to external excitation. According to the Singularity theory, the Eq. (23) is the universal unfolding of the paradigm $z^{6} + λ = 0$ . Because there are too much unfolding parameters, the bifurcation characteristics cannot be directly displayed on the plane, so it is discussed in six cases.

Case 1: $α_{1} = α_{2} = 0$ .

Bifurcation set: $B = ϕ$ : Double limit point set: $D = ϕ$ ;

Lag point set: $H_{0} = \{\frac{4096}{9375} α_{3}^{5} + α_{4}^{4} = 0\}$ ;

Transition set of system is: $Σ = B \cup H \cup D$ .

Fig. 2Transition set and bifurcation diagram of system when α1=α2=0

Case 2: $α_{1} = α_{3} = 0$ .

Bifurcation set: $B = ϕ$ ; Double limit point set: $D = ϕ$ ;

Lag point set: $H_{0} = \{\frac{729}{3125} α_{2}^{5} + α_{4}^{3} = 0\}$ ;

Transition set of system is: $Σ = B \cup H \cup D$ .

Fig. 3Transition set and bifurcation diagram of system when α1=α3=0

Case 3: $α_{1} = α_{4} = 0$ .

Bifurcation set: $B = ϕ$ ; Double limit point set: $D = ϕ$ ;

Lag point set: $H_{0} = \{\frac{729}{4096} α_{2}^{4} - α_{3}^{3} = 0\}$ ;

Transition set of system is: $Σ = B \cup H \cup D$ .

Fig. 4Transition set and bifurcation diagram of system when α1=α4=0

Case 4: $α_{2} = α_{3} = 0$ .

Bifurcation set: $B = ϕ$ ; Double limit point set: $D = ϕ$ ;

Lag point set: $H_{0} = \{\frac{512}{3125} α_{1}^{5} + α_{4}^{2} = 0\}$ ;

Transition set of system is: $Σ = B \cup H \cup D$ .

Fig. 5Transition set and bifurcation diagram of system when α2=α3=0

Case 5: $α_{2} = α_{4} = 0$ .

Bifurcation set: $B = ϕ$ ; Double limit point set: $D = ϕ$ ;

Lag point set: $H_{0} = \{- \frac{1}{3} α_{1}^{2} + α_{3} = 0\}$ ;

Transition set of system is: $Σ = B \cup H \cup D$ .

Fig. 6Transition set and bifurcation diagram of system when α2=α4=0

Case 6: $α_{3} = α_{4} = 0$ .

Bifurcation set: $B = ϕ$ ; Double limit point set: $D = ϕ$ ;

Lag point set: $H_{0} = \{\frac{128}{729} α_{1}^{3} + α_{2}^{2} = 0\}$ ;

Transition set of system is: $Σ = B \cup H \cup D$ .

Letting two unfolding parameters of $α_{1}$ , $α_{2}$ , $α_{3}$ , $α_{4}$ equal to zero, six groups of projection of system transition sets are obtained accordingly. It can be seen from above six figures: on the same projection plane, it is divided into several different subregions by transition sets. In the same subregion, the bifurcation diagrams are all topological equivalent, and any points in different subregions are not equivalent. The system transition sets are divided into four subregions I, II, III, IV and four critical lines 1, 2, 3, 4. Selecting a point respectively in each subregions and the critical lines, then simulating their bifurcation characteristic in given system conditions, finally eight bifurcation diagrams are obtained. Simulation diagrams show that the system exhibits different bifurcation characteristics in different parameter regions. The eight groups of bifurcation diagrams indicate the change process of the system, which reflect the actual state of motion of the system. The analysis of bifurcation behavior can be used to determine the parameter region which leads to the instability of the system. Therefore, appropriate selection of system parameters can suppress the vibration of rolling mill.

Fig. 7Transition set and bifurcation diagram of system when α3=α4=0

5. Numerical simulation and analysis

Selecting the actual parameters of 1780 mm rolling mill of Chengde Iron and Steel Co., Ltd as an example. The characteristic of coupling vibration model is numerically analyzed. The corresponding system parameters in the model are shown in Table 1.

Table 1Parameters of rolling mill system

Parameters	Values
Mass of up rolls ( $M$ )	1.44×10⁵ kg
Stiffness of up rolls ( $K$ )	2.08×10¹⁰ N/m
Damping of up rolls ( $C$ )	1.04×10⁶ N·s/m
Mass of rolled piece ( $m$ )	0.6318 kg
Forward stiffness of rolled piece ( $k_{f}$ )	7.58×10⁷ N/m
Backward stiffness of rolled piece ( $k_{b}$ )	1.10×10⁸ N/m
Damping of rolled piece ( $c$ )	5.20×10³ N·s/m
Width of rolled piece ( $B$ )	1.5 m
Amplitude of external excitation ( $F_{0}$ )	0.5 MN
Thickness of pre-rolling rolled piece ( $\bar{H}$ )	0.0145 m
Entrance thickness of rolled piece ( $H$ )	0.0141 m
Exit thickness of rolled piece ( $h$ )	0.0082 m
Rotational speed of work roll ( $v_{0}$ )	2.5 m/s
Roll diameter ( $D$ )	0.84 m
Density of rolled piece ( $ρ$ )	7.8×10³ kg·m^-3
Forward tension of rolled piece ( $τ_{f}$ )	3.8 MPa
Backward tension of rolled piece ( $τ_{b}$ )	5.5 MPa

In this section, frequency responses of rolling system are analyzed. In order to investigate the efficiency of appropriate rolling parameters, comparisons between different nonlinear internal parameters and external excitation amplitudes are presented. Based on the coupling vibration Eq. (21) and (22), the frequency responses of rolling system are obtained, the results are shown as Fig. 8-9.

Fig. 8Frequency responses of system with different amplitudes of external excitation

Fig. 9Frequency responses of system with different values of internal nonlinear parameters

Fig. 8(a) and Fig. 8(b) is respectively the frequency response diagram of rolled piece and rolls. As is shown in Fig. 8, with the amplitudes changing of external excitation, the jump phenomenon is observed. In order to get the vibration rules of coupling system, three different external excitation amplitudes are compared in Fig. 8(a), 8(b), the results indicate that increasing the amplitude of the external excitation, the amplitude of the vertical vibration of the roll increases and the frequency domain of the jump phenomenon decreases. Vibration rules of the roll with the change of internal nonlinear parameters are shown in Fig. 9(a)-(c). By comparison and analysis, conclusion can be drawn: increasing the absolute value of the nonlinear parameter $I_{1}$ , the system becomes unstable, and frequency response in higher frequency disappeared; increasing the absolute value of the nonlinear parameter $I_{2}$ , the resonance frequency of the vertical vibration of the roll translates to the higher frequency, resonance frequency is far away from the natural frequency of the system; increasing the absolute value of the nonlinear parameter $I_{9}$ , the amplitude-frequency curve of roll vertical vibration migrates to the right, as the result, the unstable frequency region of the system increases.

6. Conclusions

In this paper, considering the relatively horizontal vibration of rolled piece, the friction coefficient is expressed as related with the vertical vibration displacement of roll and the horizontal vibration velocity of rolled piece. Furthermore, considering the dynamic interaction between roll and rolled piece, the coupling vibration model of strip mill is proposed. In order to analyze the characteristics of coupling system, the analytic solution of the coupling system is obtained by the multiple scales method. By numerical simulation and analysis, conclusions are obtained as follows:

1) According to the Singularity theory, six groups of transition sets and bifurcation diagrams of system are obtained. The results show that: with the bifurcation parameters changing, the system shows different stability. This paper proposed a method to avoid the unstable vibration of system by limiting the selection of rolling process parameters.

2) By using MATLAB simulation, the amplitude-frequency characteristics of system with different rolling parameters are analyzed. Simulation results indicate that: with the parameters changing, the “jump phenomenon” can be weakened or even avoided; The appearance of resonance in the system will be difficult. The results can be used as a theoretical basis for selecting parameters in engineering application, and also are of significance for design of strip mill systems.

References

Yildiz S. K., Forbes J. F., Huang B., et al. Dynamic modeling and simulation of a hot strip finishing mill. Applied Mathematical Modeling, Vol. 33, Issue 7, 2009, p. 3208-3225.

Publisher
Fujita N., Kimura Y., Kobayashi K., et al. Dynamic control of lubrication characteristics in high speed tandem cold rolling. Journal of Materials Processing Technology, Vol. 229, 2015, p. 407-416.

Publisher
Kim Y., Kim C. W., Lee S., et al. Dynamic modeling and numerical analysis of a cold rolling mill. International Journal of Precision Engineering and Manufacturing, Vol. 14, Issue 3, 2013, p. 407-413.

Publisher
Brusa E, Lemma L. Numerical and experimental analysis of the dynamic effects in compact cluster mills for cold rolling. Journal of Materials Processing Technology, Vol. 209, Issue 5, 2009, p. 2436-2445.

Publisher
Heidari A., Forouzan M. 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. 27-34.

Publisher
Yarita I., Furukawa K., Seino Y. Analysis of chattering in cold rolling for ultrathin gauge steel strip. Transactions of the Iron and Steel Institute of Japan, Vol. 18, Issue 1, 1978, p. 1-10.

Publisher
Tamiya T., Furui K., Lida H., et al. Analysis of chattering phenomenon in cold rolling. Proceedings of the International Conference on Steel Rolling, Tokyo, Vol. 2, 1980, p. 1191-1202.

Search CrossRef
Johnson R., Qi Q. Chatter dynamics in sheet rolling. International Journal of Mechanical Sciences, Vol. 36, Issue 7, 1994, p. 617-630.

Publisher
Liu F., Liu B., Liu H., et al. Vertical vibration of strip mill with the piecewise nonlinear constraint arising from hydraulic cylinder. International Journal of Precision Engineering and Manufacturing, Vol. 16, Issue 9, 2015, p. 1891-1898.

Publisher
Yun I., Ehmann K., Wilson W., et al. Chatter in the strip rolling process. Part 3: chatter model. Journal of Manufacturing Science and Engineering, Vol. 120, Issue 2, 1998, p. 343-348.

Publisher
Zhang Y., Yan X., Ling Q. Electromechanical coupling vibration of rolling mill excited by variable frequency harmonic. Advanced Materials Research, Vol. 912, Issue 914, 2014, p. 662-665.

Publisher
Yang X., Tong C. N. Coupling dynamic model and control of chatter in cold rolling. Journal of Dynamic Systems Measurement and Control, Vol. 134, Issue 4, 2012, p. 271-277.

Publisher
Peng Y, Zhang Y., Sun J., et al. Tandem strip mill’s multi-parameter coupling dynamic modeling based on the thickness control. Chinese Journal of Mechanical Engineering (English Edition), Vol. 28, Issue 2, 2015, p. 353-362.

Publisher
Hu P., Ehmann K. A dynamic model of the rolling process. part I: homogeneous model. International Journal of Machine Tools and Manufacture, Vol. 40, Issue 1, 2000, p. 1-19.

Publisher
Roberts W. Friction in the hot rolling of steel strip. Iron and Steel Engineer, Vol. 51, Issue 7, 1974, p. 56-62.

Search CrossRef
Yun I., Wilson W., Ehmann K. Review of chatter studies in cold rolling. International Journal of Machine Tools and Manufacture, Vol. 38, Issue 12, 1998, p. 1499-1530.

Publisher

Cited by

Resonance Analysis of Horizontal Nonlinear Vibrations of Roll Systems for Cold Rolling Mills under Double-Frequency Excitations

Li Jiang | Tao Wang | Qing-Xue Huang

(2023)

Stability of nonlinear vibrations induced by rolling force in a precise cold mill system

Chaofan Sun | Wu Zhao | Dan Huang | Hongbin Zhang

(2023)

Research on Dynamic Characteristics of the RBBH System Based on Dynamics Model and Vibration Data Fusion

Shuilin Lin | Jianliang Sun | Chen Zhao | Yan Peng

(2022)

Vertical–Horizontal Coupling Vibration of Hot Rolling Mill Rolls under Multi-Piecewise Nonlinear Constraints

Rongrong Peng | Xingzhong Zhang | Peiming Shi

(2021)

Nonlinear vibration characteristics and time-delayed displacement control of rolling mill under dynamic rolling force

Rongrong Peng

(2021)

Vibration characteristics of multi-parametric excitations and multi-frequency external excitations of rolling mill under entry thickness fluctuation of strip

Teng Xu | Dong-xiao Hou | Zheng-nan Sun | Da-wu Guo

(2020)

Bifurcation Characteristic Research on the Load Vertical Vibration of a Hydraulic Automatic Gauge Control System

(2019)

About this article

Received

08 April 2016

Accepted

22 August 2016

Published

31 December 2016

SUBJECTS

Chaos, nonlinear dynamics and applications

DOI

https://doi.org/10.21595/jve.2016.17056

Keywords

strip mill

nonlinear vibration

friction force

multiple scales method

bifurcation

Acknowledgements

This work is supported by the Natural Science Foundation of Hebei Province, China (Grant No. E2015203349).

Author Contributions

Liu Bin contributed to the conception of the study and revised the manuscript; Jiang Jiahao contributed significantly to date analyses and wrote the manuscript; Liu Fei contributed to the manuscript preparation and programming; Liu Haoran revised the manuscript and provided fund support; Li Peng helped perform the analysis with constructive discussions.

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

Previous article in issue Previous Next article in issue Next

Research article

2021 08 06

Nonlinear vibration characteristics and time-delayed displacement control of rolling mill under dynamic rolling force

Rongrong Peng

Research article

2020 12 21

Analysis of coupling fault correlation and nonlinear vibration of multi-stage gear transmission system

Qiaorong Zhao, Xin Wang, Taotao Li, Hongwei Zhang

Research article

2017 12 31

Dynamic analysis and control of strip mill vibration under the coupling effect of roll and rolled piece

Zhaolun Liu, Jiahao Jiang, Peng Li, Guixiang Pan, Bin Liu

Research article

2017 11 15

Influence of asymmetric structure parameters on rolling mill stability

Jinlei Huang, Yong Zang, Zhiying Gao, Lingqiang Zeng

B. Liu, J. Jiang, F. Liu, H. Liu, and P. Li, “Nonlinear vibration characteristic of strip mill under the coupling effect of roll-rolled piece,” Journal of Vibroengineering, Vol. 18, No. 8, pp. 5492–5504, Dec. 2016, https://doi.org/10.21595/jve.2016.17056

Copy Extrica

Copied to clipboard!

TY  - JOUR
DO  - 10.21595/jve.2016.17056
UR  - https://doi.org/10.21595/jve.2016.17056
TI  - Nonlinear vibration characteristic of strip mill under the coupling effect of roll-rolled piece
T2  - Journal of Vibroengineering
AU  - Liu, Haoran
AU  - Liu, Bin
AU  - Jiang, Jiahao
AU  - Liu, Fei
AU  - Li, Peng
PY  - 2016
DA  - 2016/12/31
PB  - JVE International Ltd.
SP  - 5492-5504
IS  - 8
VL  - 18
SN  - 1392-8716
ER  - 

Copy Ris

Copied to clipboard!

@article{Liu_2016,
	doi = {10.21595/jve.2016.17056},
	url = {https://doi.org/10.21595/jve.2016.17056},
	year = 2016,
	month = {dec},
	publisher = {{JVE} International Ltd.},
	volume = {18},
	number = {8},
	pages = {5492--5504},
	author = {Haoran Liu and Bin Liu and Jiahao Jiang and Fei Liu and Peng Li},
	title = {Nonlinear vibration characteristic of strip mill under the coupling effect of roll-rolled piece},
	journal = {Journal of Vibroengineering}
}

Copy Bibtex

Copied to clipboard!

[1]H. Liu, B. Liu, J. Jiang, F. Liu, and P. Li, “Nonlinear vibration characteristic of strip mill under the coupling effect of roll-rolled piece,” Journal of Vibroengineering, vol. 18, no. 8, pp. 5492–5504, Dec. 2016, doi: 10.21595/jve.2016.17056.

Copy IEEE

Copied to clipboard!

Liu, Haoran, Bin Liu, Jiahao Jiang, Fei Liu, and Peng Li. “Nonlinear Vibration Characteristic of Strip Mill under the Coupling Effect of Roll-Rolled Piece.” Journal of Vibroengineering 18, no. 8 (December 31, 2016): 5492–5504. https://doi.org/10.21595/jve.2016.17056.

Copy Chicago

Copied to clipboard!

Nonlinear vibration characteristic of strip mill under the coupling effect of roll-rolled piece

Abstract

1. Introduction

2. Mathematical modeling of the strip mill structure

2.1. Modeling of friction force and rolling force

2.2. Dynamic model of strip mill mechanical structure

3. Coupling system solution of strip mill

4. Bifurcation characteristic analysis of strip mill

5. Numerical simulation and analysis

6. Conclusions

References

Cited by

About this article

Related Articles