Published: 31 December 2016

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

Bin Liu1
Jiahao Jiang2
Fei Liu3
Haoran Liu4
Peng Li5
1, 2, 5School of Electrical Engineering, Yanshan University, Qinhuangdao, China
3, 4School of Informational Science and Engineering, Yanshan University, Qinhuangdao, China
Corresponding Author:
Haoran Liu
Views 49
Reads 25
Downloads 1466

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=BlcQpKTK,

where:

lc=RΔh, Qp=1.08+1.79με1-εRH-1.02ε,
KT=1-0.7τb+0.3τfK, K=1.15a013ε0+23ε1a1,

where, B is the width of rolled piece; lc is the contact length between roll and rolled piece in deformation zone; Qp is the influential coefficient in stressed state; KT is the tension coefficient; K is the average deformation resistance of materials; a0, a1 is the regression coefficient of the model; R is the work roll radius; Δh is reduction quantity of rolled piece, Δh=H-h-2y, H is the entrance thickness of rolled piece, h is the exit thickness of rolled piece,yis the vibration displacement of rolls;εis the reduction rate of rolled piece, ε=Δh/H; ε0 is the reduction rate of frame ingress, ε0=(H--H)/H-; ε1 is the reduction rate of frame egress, ε1=(H--h)/H-, 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
μ=ΔhD0.5+K1-0.5e-K2v0+x˙,

where, K1 and K2 are the friction characteristic coefficients; D is the work roll diameter; v0 is work roll rotational speed; x˙ is the relative vibration speed of rolled piece in roll gap; considering that K2x˙1, so the Roberts formula can be simplified as:

3
μ=μ0,0+Δμ,

where:

μ0,0=H-hD0.5+K1-0.5e-K2v0,Δμ=b1x˙+b2y+b3x˙2+b4x˙y+b5y2+b6x˙3+b7x˙2y+b8x˙y2+b9y3,

where:

b0=μ0,0, b1=x˙μ0,0, b2=yμ0,0, b3=12!2x˙2μ0,0,
b4=2x˙yμ0,0, b5=12!2y2μ0,0, b6=13!3x˙3μ0,0,
b7=123x˙2yμ0,0, b8=123y2x˙μ0,0, b9=13!3y3μ0,0.

When rolling process is stable, 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=F0,0+ΔF,

where:

F0,0=1.15a0BL1L2RH-h,ΔF=I1x˙+I2y+I3x˙2+I4x˙y+I5y2+I6x˙3+I7x˙2y+I8x˙y2+I9y3,

where:

L1=13ε0+23ε1a1-0.7τb-0.3τf, L2=1.08+1.79b0e-K2v0hR-1.02HH2H-h,
I0=F0,0, I1=x˙F0,0, I2=yF0,0, I3=12!2x˙2F0,0,
I4=2x˙yF0,0, I5=12!2y2F0,0, I6=13!3x˙3F0,0,
I7=123x˙2yF0,0, I8=123y2x˙F0,0, I9=13!3y3F0,0.

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

5
Ff=2μ0+Δμ×F0,0+ΔF,

where, μ0 and F0,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
ΔFf=2μ0ΔF+2ΔμF0.

By substituting Eq. (3) and Eq. (4) into Eq. (6), the expression of ΔFf can be obtained:

7
ΔFf=2b0I1x˙+I2y+I3x˙2+I4x˙y+I5y2+I6x˙3+I7x˙2y+I8x˙y2+I9y3
+2I0b1x˙+b2y+b3x˙2+b4x˙y+b5y2+b6x˙3+b7x˙2y+b8x˙y2+b9y3.

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

Structure 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=BlcH+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: kf=τfSf/Lfd, kb=τbSb/Lbd; where τf and τb stand for forward and backward tension of rolled piece, Sf and Sb stand for forward and backward cross-sectional area of rolled piece, Lfd and Lbd 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. m1 and m2 are respectively the equivalent mass of upper and lower rolls; k1 and k2 are respectively the equivalent stiffness of upper and lower rolls; c1 and c2 are respectively the equivalent dumping of upper and lower rolls;y1andy2are 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
m1y¨1+c1y˙1+k1y1=ΔF+F*,m2y¨2+c2y˙2+k2y2=-ΔF-F*.

Based on the assumption that the structure and vibration characteristic of rolling mill is symmetrical in relation to rolling line, there exist m1=m2; c1=c2; y1=-y2; k1=k2. 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
My¨+Cy˙+Ky=Δ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
mx¨+2cx˙+kf+kbx=ΔFf.

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

11
mx¨+2cx˙+(kf+kb)x=2b0I1x˙+I2y+I3x˙2+I4x˙y+I5y2+I6x˙3+I7x˙2y+I8x˙y2 +I9y3+2I0b1x˙+b2y+b3x˙2+b4x˙y+b5y2+b6x˙3+b7x˙2y+b8x˙y2+b9y3,My¨+Cy˙+Ky=I1x˙+I2y+I3x˙2+I4x˙y+I5y2+I6x˙3+I7x˙2y+I8x˙y2+I9y3+F*.

3. Coupling system solution of strip mill

Assuming that the system is subjected to periodic external disturbances, set F*=F0cosΩt. By transposition and replacement, Eq. (11) is transformed into a standard form:

12
x¨+ω102x=η11f1x˙,y+η12f2x˙,y-β1x˙,y¨+ω202y=η21f1x˙,y+η22cosΩt-β2y˙,

where:

f1x˙,y=I1x˙+I2y+I3x˙2+I4x˙y+I5y2+I6x˙3+I7x˙2y+I8x˙y2+I9y3,
f2x˙,y=b1x˙+b2y+b3x˙2+b4x˙y+b5y2+b6x˙3+b7x˙2y+b8x˙y2+b9y3,

where:

β1=2cm, ω10=(kf+kb)m, η11=2b0m, η12=2I0m,
β2=CM, ω20=KM, η21=1M, η22=1M.

Set, β1=εβ^1; β2=εβ^2; η11=εη^11; η12=εη^12; η13=εη^13; η21=εη^21; η22=εη^22.

By parameters replacement, Eq. (12) becomes:

13
x¨+ω102x=εη^11f1(x˙,y)+η^12f2(x˙,y)-β^1x˙,y¨+ω202y=εη^21f1(x˙,y)+η^22F0cosΩt-β^2y˙.

For finding a nonlinear approximate solution of Eq. (13), two time scales of T0=t and T1=εt are selected. Then, the time derivatives are defined as:

14
ddt=D0+εD1+,d2dt2=D02+2εD0D1+,

where, Dn=/Tn, ε is small parameter. By using multiple scales method, set the solution of Eq. (14) as:

15
x=x0T0,T1+εx1T0,T1,y=y0T0,T1+εy1T0,T1.

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

16
D02x0+ω102x0=0,D02y0+ω202y0=0,
17
D02x1+ω102x1=η^11I1D0x0+I2y0+I3(D0x0)2+I4D0x0y0+I5y02+I6(D0x0)3+I7(D0x0)2y0+I8D0x0y02+I9y03 +η^12b1D0x0+b2y0+b3(D0x0)2+b4D0x0y0+b5y02+b6(D0x0)3+b7(D0x0)2y0+b8D0x0y02+b9y03-β^1D0x0-2D0D1x0,D02y1+ω202y1=η^21I1D0x0+I2y0+I3(D0x0)2+I4D0x0y0+I5y02+I6(D0x0)3+I7(D0x0)2y0+I8D0x0y02+I9y03 -β^2D0y0+η^22F0cosΩT0-2D0D1y0.

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

18
x0=AT1eiω10T0+A-T1e-iω10T0,y0=BT1eiω20T0+B-T1e-iω20T0.

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
η^11I1Aiω10+I2Be-iσ1T1+3I6A2A-iω103+2I7AA-Bω102e-iσ1T1+2I8ABB-iω10+3I9B2B-e-iσ1T1 +η^12b1Aiω10+b2Be-iσ1T1+3b6A2A-iω103+2b7AA-Bω102e-iσ1T1+2b8ABB-iω10+3b9B2B-e-iσ1T1 -β^1Aiω10-2D1Aiω10=0,η^21I1Aiω10eiσ1T1+I2B+3I6A2A-iω103eiσ1T1+2I7AA-Bω102+2I8ABB-iω10eiσ1T1+3I9B2B- -β^ 2 Biω20+η^22F02eiσ1T1-2D1Biω20=0.

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

A=12aeiθ1, B=12beiθ2,

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

20
a˙=η^112I1a+34I6a3ω102+12I8ab2+η^122b1a+34b6a3ω102+12b8ab2 -12β^1a+η^11ω1012I2b+14I7a2bω102+38I9b3+η^12ω1012b2b+14b7a2bω102+38b9b3sinφ,aφ˙+ϕ˙=η^11ω1012I2b+14I7a2bω102+38I9b3+η^12ω1012b2b+14b7a2bω102+38b9b3cosφ+aσ-σ1,b˙=η^21ω2012I1aω10+38I6a3ω103+14I8ab2ω10cosφ-12β^ 2 b+η^22F02ω20sinϕ,bϕ˙=η^21ω2012I1aω10+38I6a3ω103+14I8ab2ω10sinφ+12I2b+14I7a2bω102+38I9b3+η^22F02ω20cosϕ+bσ.

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

21
η^112I1a+34I6a3ω102+12I8ab2+η^122b1a+34b6a3ω102+12b8ab2-12β^1a2
+aσ-σ12
=η^11ω1012I2b+14I7a2bω102+38I9b3+η^12ω1012b2b+14b7a2bω102+38b9b32,
22
η^21ω2012I1aω10+38I6a3ω103+14I8ab2ω10cosϕ-12β^2b2
+η^21ω2012I1aω10+38I6a3ω103+14I8ab2ω10sinϕ+12I2b+14I7a2bω102+38I9b3+bσ2=η^22F02ω202,

where:

cosϕ=-aσ-σ1η^11ω1012I2b+14I7a2bω102+38I9b3+η^12ω1012b2b+14b7a2bω102+38b9b3,
sinϕ=-η^112I1a+34I6a3ω102+12I8ab2+η^122b1a+34b6a3ω102+12b8ab2-12β^1aη^11ω1012I2b+14I7a2bω102+38I9b3+η^12ω1012b2b+14b7a2bω102+38b9b3.

4. Bifurcation characteristic analysis of strip mill

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

23
z6+α1z4+α2z3+α3z2+α4z+λ=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 z6+λ=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: H0=40969375α35+α44=0;

Transition set of system is: Σ=BHD.

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

Transition set and bifurcation diagram of system when α1=α2=0
Transition 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: H0=7293125α25+α43=0;

Transition set of system is: Σ=BHD.

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

Transition set and bifurcation diagram of system when α1=α3=0
Transition 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: H0=7294096α24-α33=0;

Transition set of system is: Σ=BHD.

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

Transition set and bifurcation diagram of system when α1=α4=0
Transition 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: H0=5123125α15+α42=0;

Transition set of system is: Σ=BHD.

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

Transition set and bifurcation diagram of system when α2=α3=0
Transition 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: H0=-13α12+α3=0;

Transition set of system is: Σ=BHD.

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

Transition set and bifurcation diagram of system when α2=α4=0
Transition 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: H0=128729α13+α22=0;

Transition set of system is: Σ=BHD.

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

Transition set and bifurcation diagram of system when α3=α4=0
Transition 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×105 kg
Stiffness of up rolls (K)
2.08×1010 N/m
Damping of up rolls (C)
1.04×106 N·s/m
Mass of rolled piece (m)
0.6318 kg
Forward stiffness of rolled piece (kf)
7.58×107 N/m
Backward stiffness of rolled piece (kb)
1.10×108 N/m
Damping of rolled piece (c)
5.20×103 N·s/m
Width of rolled piece (B)
1.5 m
Amplitude of external excitation (F0)
0.5 MN
Thickness of pre-rolling rolled piece (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 (v0)
2.5 m/s
Roll diameter (D)
0.84 m
Density of rolled piece (ρ)
7.8×103 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

Frequency responses of system with different amplitudes of external excitation
Frequency responses of system with different amplitudes of external excitation

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

Frequency responses of system with different values of internal nonlinear parameters
Frequency responses of system with different values of internal nonlinear parameters
Frequency 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 I1, the system becomes unstable, and frequency response in higher frequency disappeared; increasing the absolute value of the nonlinear parameter I2, 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 I9, 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • Johnson R., Qi Q. Chatter dynamics in sheet rolling. International Journal of Mechanical Sciences, Vol. 36, Issue 7, 1994, p. 617-630.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • Roberts W. Friction in the hot rolling of steel strip. Iron and Steel Engineer, Vol. 51, Issue 7, 1974, p. 56-62.
  • 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.

About this article

Received
08 April 2016
Accepted
22 August 2016
Published
31 December 2016
SUBJECTS
Chaos, nonlinear dynamics and applications
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.