Published: 31 December 2017

# An analysis on strip vibration coupled with torsional vibration of main drive system of rolling mill

Chongyi Gao1
Guojun Du2
Rui Li3
Xiaoqiang Guo4
1College of Mechanical Engineering, Yanshan University, Qinhuangdao, China
1, 2, 3, 4Key Laboratory of Mechanics Reliability for Heavy Equipment and Large Structure of Hebei Province, Yanshan University, Qinhuangdao, China
Corresponding Author:
Guojun Du
Views 225

#### Abstract

According to the movement mechanism of strip and rollers during the continuous rolling process, the main drive system of each stand was simplified to a single degree of freedom discrete model, and the strip was simplified to an axially moving Euler beam. Then, a nonlinear continuous-discrete coupled vibration model between transverse and longitudinal vibrations of strip and torsional vibration of main drive system was established. According to Hamilton’s principle, the nonlinear differential equations were established. Moreover, modified iteration method and Kantorovich averaging method were used to solve the differential equations. Depending on numerical calculation, the amplitude-frequency responses of strip vibration coupled with torsional vibration of main drive system were obtained. Finally, the influences of the axial velocity, the strip tension, the torsional stiffness, and the rotational inertia on the vibration characteristics were discussed. The results would provide a theoretical reference for control and analysis of rolling mill vibration in engineering practice.

## 1. Introduction

The rolling mill vibration has brought serious troubles to iron and steel industry for many years, which is a significant technical problem. The vibration forms are various, and characteristics and causes are different. At present, the research on the different forms of rolling mill vibration mainly includes vertical vibration and horizontal vibration of rollers, transverse and longitudinal vibration of strip [1, 2], torsional vibration of main drive system [3] and the axial vibration, etc. As early as in 1967, the research of various vibrations of rolling mill was begun. Moller and Hoggart et al. [4] found the existence of torsional vibration of the rolling mill when they were employing the two rolls test machine, and the vibration behavior, which was considered as the self-excited vibration. Lawrence and Thomas [5] analyzed the chatter-mark on the strip surface, which caused torsional vibration by defect of transmission gear. Swiatoniowski and Bar [6] analyzed the chattering phenomenon during the rolling process, and the mathematical model of self-excited vibration was established. Wang et al. [7] analyzed the influence of rolling force on horizontal stiffness of rolling mill, and the functional relationship among the vertical vibration of rolling mill, the horizontal vibration of work roll and the torsional vibration of main drive system was established. Sun et al. [8] investigated the torsional vibration on the dynamic model of six rollers system of 1100 rolling mill, and the influence of torsional vibration on strip shape was acquired. Xu et al. [9, 10] established the hybrid system model of strip coupled with rollers, and the coupled vibration mechanism of roller and strip with tension fluctuation was obtained. Zou et al. [11] analyzed the different influences of cold rolling and hot rolling on axial force, and a mathematical model was established. Du et al. [12] established the nonlinear dynamic model of coupled model between rollers and strip, and the inertial boundary condition was proposed.

However, there are few researches studying on coupled vibration of multiple systems. In the earlier stage, the research of strip vibration coupled with torsional vibration of main drive system has been investigated by our group, such as the vibration analysis of main drive system of each stand by finite element software. Researchers found the energy was transferred from one to another main drive system through the strip when resonance occurred, which resulted to the great fluctuation of rollers, and the vibration amplitude may exceed self-resonance [13]. Therefore, it is significant to analyze the coupled vibration between strip and main drive system of rolling mill furtherly.

In this paper, the main drive system of rolling mill is simplified to a single degree of freedom model of discrete system, and the strip is simplified to an axially moving Euler beam [14, 15]. Then, a continuous-discrete coupled vibration model between transverse and longitudinal vibrations of strip and torsional vibration of rolling mill main drive system is established. The research results can provide important theoretical reference for control and analysis of strip vibration coupled with torsional vibration during the continuous rolling process.

## 2. Mechanical and mathematical models

The main drive system of each stand can be regarded as a mass-spring system in the tandem mill, which mainly composed of inertial components of electromotor, reduction gears and rollers etc., and elastic components of connecting shafts, a dimensional geometric model is shown in Fig. 1. Due to the traction behavior of motor drive, the motor of main drive system could be regarded as a fixed end [16]. And the rotational inertia of roller is much bigger than the other parts, so the main drive system can be simplified to a single degree of freedom discrete model of spring-mass system. The rollers can be equivalent to the rigid inertial components, $j$ is the rotational inertia and $\theta$ is the rotation angle; the connecting shafts can be equivalent to the elastic element, $k$ is the torsional stiffness, as shown in Fig. 2. When transverse and longitudinal vibrations of strip are considered, the mechanical geometry characteristic is similar to an axially moving beam [17, 18]. Also, the strip is an elastic continuum whose thickness is much less than the length, the strip is equivalent to the isotropic Euler beam based on the theory of axially moving beam. The equivalent mechanical model is shown in Fig. 3, and assumes that the strip with uniform motion is running, ${v}_{0}$is the axial velocity. It hypothesizes that there is no relative motion between rollers and strip, and the upper and lower rollers are symmetric discs of fixed axis rotation along the width of beam. The transverse displacement and the longitudinal displacement of Euler beam are $w\left({x}_{0},{y}_{0},t\right)$ and $u\left({x}_{0},{y}_{0},t\right)$ respectively; $l$ is the distance of Euler beam between two stands; the left tension and the right tension are ${P}_{1}$ and ${P}_{2}$ respectively, and ${P}_{1}={P}_{2}={P}_{0}$.

Fig. 1A dimensional geometric model of main drive system of one stand: 1 – electromotor, 2 – intermediate shaft, 3 – reduction gears, 4 – profile spindle, 5 – working roll, 6 – back-up roll, 7 – screwdown device, 8 – balancing device, 9 – unjamming gear, 10 – mill housing, 11 – platform

Based on Hamilton principle, the mathematical model of strip vibration coupled with torsional vibration of main drive system is established.

The kinetic energy ${T}_{1}$ of axially moving Euler beam can be written as:

1
$\begin{array}{l}{T}_{1}=\frac{1}{2}{\int }_{0}^{l}\rho A\left[{\left({v}_{0}+\frac{du}{dt}\right)}^{2}+{\left(\frac{dw}{dt}\right)}^{2}\right]d{x}_{0}\\ =\frac{1}{2}{\int }_{0}^{l}\rho A\left[\left({v}_{0}+{u}_{,t}+{v}_{0}{u}_{,{x}_{0}}{\right)}^{2}+\left({w}_{,t}+{v}_{0}{w}_{,{x}_{0}}{\right)}^{2}]d{x}_{0},\end{array}$

where $\rho$ is the density of strip, $A$ is the cross-sectional area of Euler beam, ${u}_{,t}$, ${w}_{,t}$ and ${u}_{,{x}_{0}}$, ${w}_{,{x}_{0}}$ are the $u\left({x}_{0},{y}_{0},t\right)$ and $w\left({x}_{0},{y}_{0},t\right)$ on variables $t$ and ${x}_{0}$ of the first partial derivatives respectively.

Fig. 2A simplified model of main drive system of rolling mill

Fig. 3A mechanical model of strip between main drive systems of two stands

The kinetic energy ${T}_{2}$ of roller can be written as:

2
${T}_{2}=\frac{1}{2}j{\left(\frac{du}{dt}\frac{1}{r}\right)}^{2}=\frac{j}{2{r}^{2}}\left({u}_{,t}+{v}_{0}{u}_{,{x}_{0}}{\right)}^{2},$

where $r$ is the radius of roller.

The potential energy ${U}_{1}$ of Euler beam can be written as:

3
$\begin{array}{l}{U}_{1}=\frac{1}{2}\iiint E{\epsilon }_{{x}_{0}}^{2}d{x}_{0}d{y}_{0}d{z}_{0}\\ \end{array}=\frac{1}{2}{\int }_{0}^{l}\left[EA{\left({u}_{,{x}_{0}}+\frac{1}{2}{w}_{,{x}_{0}}^{2}\right)}^{2}+EI{w}_{,{x}_{0}{x}_{0}}^{2}\right]d{x}_{0},$

where ${\epsilon }_{{x}_{0}}$ is the strain of Euler beam, $E$ is the Young’s modulus and $I$ is the moment of inertia.

The potential energy ${U}_{2}$ of tension can be written as:

4
${U}_{2}={\int }_{0}^{l}{P}_{0}\left({u}_{,{x}_{0}}+\frac{1}{2}{w}_{,{x}_{0}}^{2}\right)d{x}_{0}.$

Then according to Hamilton principle, the following equation can be obtained:

5
$\delta {\int }_{{t}_{2}}^{{t}_{1}}\left({T}_{1}+{T}_{2}-{U}_{1}-{U}_{2}\right)dt=0.$

Substitution of Eqs. (1-4) into Eq. (5), the motion equations can be expressed as:

6
$\rho A{u}_{,tt}+\left(\rho A{v}_{0}^{2}-EA\right){u}_{,{x}_{0}{x}_{0}}+2\rho A{v}_{0}{u}_{,{x}_{0}t}-EA{w}_{,{x}_{0}}{w}_{,{x}_{0}{x}_{0}}=0,$
7
$\begin{array}{l}\rho A{w}_{,tt}+EI{w}_{,{x}_{0}{x}_{0}{x}_{0}{x}_{0}}+\left(\rho A{v}_{0}^{2}-{P}_{0}\right){w}_{,{x}_{0}{x}_{0}}+2\rho A{v}_{0}{w}_{,{x}_{0}t}\\ -EA\left({w}_{,{x}_{0}{x}_{0}}{u}_{,{x}_{0}}+{u}_{,{x}_{0}{x}_{0}}{w}_{,{x}_{0}}+\frac{3}{2}{w}_{,{x}_{0}}^{2}{w}_{,{x}_{0}{x}_{0}}\right)=0.\end{array}$

In general, the axial kinetic energy caused by transverse vibration is relatively small. Let ${u}_{,t}={u}_{,tt}=0$, So, Eq. (6) is simplified as:

8
$\left(\rho A{v}_{0}^{2}EA\right){u}_{,{x}_{0}{x}_{0}}-EA{w}_{,{x}_{0}}{w}_{,{x}_{0}{x}_{0}}=0.$

The relationship between torque and torsional stiffness is expressed as: ${M}_{i}={k}_{i}\mathrm{\Delta }{\theta }_{i}$, where ${\theta }_{i}$ is relative torsional angle between two shaft sections, ${k}_{i}$ is the torsional stiffness of shaft segment. Therefore, boundary conditions can be written as:

When ${x}_{0}=0$:

9
$EAr{u}_{,{x}_{0}}-{P}_{0}r=-j\frac{{d}^{2}u}{d{t}^{2}}\frac{1}{r}-k\frac{u}{r}.$

When ${x}_{0}=l$:

10
$EAr{u}_{,{x}_{0}}+{P}_{0}r=-j\frac{{d}^{2}u}{d{t}^{2}}\frac{1}{r}-k\frac{u}{r},$
11
$w\left(0\right)=w\left(l\right)={w}_{,{x}_{0}{x}_{0}}\left(0\right)={w}_{,{x}_{0}{x}_{0}}\left(l\right)=0,$
12
${w}_{\mathrm{m}\mathrm{a}\mathrm{x}}={\phi }_{m}.$

The time variable of equation is separated from the space variable by substitution of $w={\varphi }_{0}\left({x}_{0}\right)\mathrm{c}\mathrm{o}\mathrm{s}{\omega }_{0}t$ into Eq. (8), the following equation can be obtained as: $u={\phi }_{0}\left({x}_{0}\right)\mathrm{c}\mathrm{o}{\mathrm{s}}^{2}{\omega }_{0}t$, where ${\omega }_{0}$ is vibration frequency.

Based on Kantorovich averaging method on the interval $\left[0,2\pi /{\omega }_{0}\right]$, the time variable can be eliminated. The motion equations can be simplified as:

13
$\left(\rho A{v}_{0}^{2}-EA\right){{\varphi }_{0}}_{,{x}_{0}{x}_{0}}-EA{{\phi }_{0}}_{,{x}_{0}}{{\phi }_{0}}_{,{x}_{0}{x}_{0}}=0,$
14
$EI{{\phi }_{0}}_{,{x}_{0}{x}_{0}{x}_{0}{x}_{0}}+\left(\rho A{v}_{0}^{2}-{P}_{0}\right){{\phi }_{0}}_{,{x}_{0}{x}_{0}}-{\omega }_{0}^{2}\rho A{\phi }_{0}$
$-\frac{3}{4}EA\left({{\phi }_{0}}_{,{x}_{0}{x}_{0}}{{\varphi }_{0}}_{,{x}_{0}}+{{\varphi }_{0}}_{,{x}_{0}{x}_{0}}{{\phi }_{0}}_{,{x}_{0}}+\frac{3}{2}{{\phi }_{0}}_{,{x}_{0}}^{2}{{\phi }_{0}}_{,{x}_{0}{x}_{0}}\right)=0.$

Boundary conditions can be written as, when ${x}_{0}=0$:

15
$EA{r}^{2}{{\varphi }_{0}}_{,{x}_{0}}+j{v}_{0}^{2}{{\varphi }_{0}}_{,{x}_{0}{x}_{0}}-\frac{4}{3}j{w}_{0}^{2}{\varphi }_{0}-\frac{4}{3}{P}_{0}r+k\frac{{\varphi }_{0}}{r}=0.$

When ${x}_{0}=l$:

16
$EA{r}^{2}{{\varphi }_{0}}_{,{x}_{0}}+j{v}_{0}^{2}{{\varphi }_{0}}_{,{x}_{0}{x}_{0}}-\frac{4}{3}j{w}_{0}^{2}{\varphi }_{0}+\frac{4}{3}{P}_{0}r+k\frac{{\varphi }_{0}}{r}=0.$
17
${\phi }_{0}\left(0\right)={\phi }_{0}\left(l\right)={{\phi }_{0}}_{,{x}_{0}{x}_{0}}\left(0\right)={{\phi }_{0}}_{,{x}_{0}{x}_{0}}\left(l\right)=0,$
18
${\phi }_{0}\left(\frac{1}{2}\right)={{\phi }_{0}}_{m}.$

The dimensionless quantities are given by:

$x=\frac{{x}_{0}}{l},\phi =\frac{{\phi }_{0}}{l},\varphi =\frac{{\varphi }_{0}}{l},v={v}_{0}\sqrt{\frac{\rho }{E}},\omega ={\omega }_{0}\sqrt{\frac{\rho A{l}^{4}}{EI}},$
$S=\frac{A{l}^{2}}{I},J=\frac{j}{\rho A{r}^{2}l},P=\frac{{P}_{0}{l}^{2}}{EI},K=\frac{kl}{EA{r}^{3}}.$

Then, the dimensionless form of Eqs. (13-18) can be obtained respectively:

19
$\left({v}_{0}^{2}-1\right){\varphi }_{,xx}-{\phi }_{,x}{\phi }_{,xx}=0,$
20
${\phi }_{,xxxx}+\left(S{v}^{2}-P\right){\phi }_{,xx}-{\omega }^{2}\phi -\frac{3}{4}S\left({\phi }_{,xx}{\varphi }_{,x}+{\varphi }_{,xx}{\phi }_{,x}+\frac{3}{2}{\phi }_{,x}^{2}{\phi }_{,xx}\right)=0.$

When $x=0$:

21
${\varphi }_{,x}+J{v}^{2}{\varphi }_{,xx}-\frac{4J}{3S}{\omega }^{2}\varphi -\frac{4P}{3S}+K\varphi =0.$

When $x=1$:

22
${\varphi }_{,x}+J{v}^{2}{\varphi }_{,xx}-\frac{4J}{3S}{\omega }^{2}\varphi +\frac{4P}{3S}+K\varphi =0,$
23
$\phi \left(0\right)=\phi \left(1\right)={\phi }_{,xx}\left(0\right)={\phi }_{,xx}\left(1\right)=0,$
24
$\phi \left(\frac{1}{2}\right)={\phi }_{m}.$

## 3. Analytical solution of vibration equations

Due to more complicated solving process of the Eqs. (19-24), the modified iteration method is used to solve the equations in this section.

### 3.1. The first-order approximate solution

Firstly, all the nonlinear terms of the Eq. (20) are omitted, the simplified motion equation can be obtained as follow:

25
${{\phi }_{1}}_{,xxxx}-{\omega }^{2}{\phi }_{1}=0.$

The series solution of Eq. (25) is:

26
${\phi }_{1}\left(x\right)={a}_{0}{M}_{0}\left(x\right)+{a}_{1}{N}_{0}\left(x\right)+{a}_{2}{I}_{0}\left(x\right)+{a}_{3}{K}_{0}\left(x\right),$

where:

${M}_{0}=\sum _{n=0}^{\infty }\frac{\left({\omega }^{2}{\right)}^{n}{x}^{4n}}{\left(4n\right)!},{N}_{0}=\sum _{n=0}^{\infty }\frac{\left({\omega }^{2}{\right)}^{n}{x}^{4n+1}}{\left(4n+1\right)!},{I}_{0}=\sum _{n=0}^{\infty }\frac{\left({\omega }^{2}{\right)}^{n}{x}^{4n+2}}{\left(4n+2\right)!},{K}_{0}=\sum _{n=0}^{\infty }\frac{\left({\omega }^{2}{\right)}^{n}{x}^{4n+3}}{\left(4n+3\right)!}.$

By substitution of Eq. (26) into Eqs. (23-24), one has ${\omega }_{1}=$ 18.71. Then, the coefficients of Eq. (26) have:

${a}_{0}=0,{a}_{1}={\mu }_{1}{\phi }_{m},{a}_{2}=0,{a}_{3}={\mu }_{2}{\phi }_{m}.$

Then:

27
${\phi }_{1}\left(x\right)={\phi }_{m}\left[{\mu }_{1}{N}_{0}\left(x\right)+{\mu }_{2}{K}_{0}\left(x\right)\right],$

where:

${\mu }_{1}=\frac{-{K}_{0}\left(1\right)}{{N}_{0}\left(1\right){K}_{0}\left(\frac{1}{2}\right)-{N}_{0}\left(\frac{1}{2}\right){K}_{0}\left(1\right)},{\mu }_{2}=\frac{{N}_{0}\left(1\right)}{{N}_{0}\left(1\right){K}_{0}\left(\frac{1}{2}\right)-{N}_{0}\left(\frac{1}{2}\right){K}_{0}\left(1\right)}.$

And substituting of Eq. (27) into Eq. (19), one can be got:

28
${\varphi }_{1}\left(x\right)=\frac{1}{2\left({v}^{2}-1\right)}\int {\phi }_{1,x}^{2}dx+{c}_{1}x+{c}_{2}.$

The coefficients ${c}_{1}$ and ${c}_{2}$ can be determined by Eqs. (21-22), and the results are as follows:

29
${c}_{1}=\frac{3S}{4J{\omega }^{2}-3SK}\left\{\frac{1}{2\left({v}^{2}-1\right)}{\left[2J{v}^{2}{\phi }_{1,x}{\phi }_{1,xx}-\left(\frac{4J}{3S}{\omega }^{2}-K\right)\int {\phi }_{1,x}^{2}dx+{\phi }_{1,x}^{2}\right]}_{x=1}+\frac{8P}{3S}\right\},$
30
$\begin{array}{l}{c}_{2}=\frac{3S}{4J{\omega }^{2}-3SK}\left\{\frac{1}{2\left({v}^{2}-1\right)}{\left[2J{v}^{2}{\phi }_{1,x}{\phi }_{1,xx}-\left(\frac{4J}{3S}{\omega }^{2}-K\right)\int {\phi }_{1,x}^{2}dx+{\phi }_{1,x}^{2}\right]}_{x=1}\right\\\ -\frac{4P}{3S}+\frac{8P}{4J{\omega }^{2}-3SK}}.\end{array}$

Then, by substitution of ${c}_{1}$ and ${c}_{2}$ into ${\varphi }_{1}\left(x\right)$, the first-order approximate solution can be obtained.

### 3.2. The second-order modified-iterative solution

In this subsection, the second-order modified-iterative solution can be given.

By substitution of Eqs. (27-28) into Eq. (20), one has:

31
${\phi }_{2,xxxx}+\left(S{v}^{2}-P\right){\phi }_{1,xx}-{\omega }^{2}{\phi }_{2}-\frac{3}{4}S\left({{\phi }_{1}}_{,xx}{{\varphi }_{1}}_{,x}+{{\phi }_{1}}_{,x}{{\varphi }_{1}}_{,xx}+\frac{3}{2}{{\phi }_{1}}_{,x}^{2}{{\phi }_{1}}_{,xx}\right)=0.$

That is:

32
${\phi }_{2,xxxx}-{\omega }^{2}{\phi }_{2}=\alpha {{\phi }_{1}}_{,xx}+\beta {{\phi }_{1}}_{,x}^{2}{{\phi }_{1}}_{,xx},$

where:

$\alpha =P+\frac{3}{4}S{c}_{1}-S{v}^{2},\beta =\frac{9S{v}^{2}}{8\left({v}^{2}-1\right)},{\phi }_{{1}_{,x}}^{2}\left(x\right)={\phi }_{m}^{2}\left(\sum _{n=0}^{\infty }{A}_{n}^{\left(1\right)}{x}^{4n}+\sum _{n=0}^{\infty }{B}_{n}^{\left(1\right)}{x}^{4n+2}\right),$
${\phi }_{{1}_{,xx}}\left(x\right)={\phi }_{m}\left(\sum _{n=1}^{\infty }{C}_{n}^{\left(1\right)}{x}^{4n-1}+\sum _{n=0}^{\infty }{D}_{n}^{\left(1\right)}{x}^{4n+1}\right),$

where:

${A}_{0}^{\left(1\right)}={\mu }_{1}^{2},$
${A}_{n}^{\left(1\right)}=\sum _{{n}_{1}=0}^{n}\frac{{\mu }_{1}^{2}{\omega }_{1}^{2n}}{\left(4{n}_{1}\right)!\left(4n-4{n}_{1}\right)!}+\sum _{{n}_{1}=0}^{n-1}\frac{{\mu }_{2}^{2}{\omega }_{1}^{2n-2}}{\left(4{n}_{1}+2\right)!\left(4n-4{n}_{1}-2\right)!},\left(n=1,2\dots \right),$
${B}_{n}^{\left(1\right)}=\sum _{{n}_{1}=0}^{n}\frac{2{\mu }_{1}{\mu }_{2}{\omega }_{1}^{2n}}{\left(4{n}_{1}\right)!\left(4n-4{n}_{1}+2\right)!},\left(n=0,1\dots \right),$
${C}_{n}^{\left(1\right)}=\frac{{\mu }_{1}{\omega }_{1}^{2n}}{\left(4n-1\right)!}\begin{array}{ll},& \left(n=1,2\dots \right)\end{array},{D}_{n}^{\left(1\right)}=\frac{{\mu }_{2}{\omega }_{1}^{2n}}{\left(4n+1\right)!}\begin{array}{ll},& \left(n=0,1...\right)\end{array}.$

With the above equations, that is obtained as follow:

33
${\phi }_{{1}_{,x}}^{2}{\phi }_{{1}_{,xx}}={\phi }_{m}^{3}\left(\sum _{n=1}^{\infty }{A}_{n}^{\left(2\right)}{x}^{4n-1}+\sum _{n=0}^{\infty }{B}_{n}^{\left(2\right)}{x}^{4n+1}+\sum _{n=0}^{\infty }{C}_{n}^{\left(2\right)}{x}^{4n+3}\right),$

where:

${A}_{n}^{\left(2\right)}=\sum _{m=0}^{n}\left({A}_{n}^{\left(1\right)}{C}_{n-m}^{\left(1\right)}\right),\mathrm{}\mathrm{}\mathrm{}\left(n=1,2\dots \right),{B}_{0}^{\left(2\right)}={\mu }_{1}^{2}{\mu }_{2},$
${B}_{n}^{\left(2\right)}=\sum _{m=0}^{n}\left({A}_{n}^{\left(1\right)}{D}_{n-m}^{\left(1\right)}+{B}_{n}^{\left(1\right)}{C}_{n-m}^{\left(1\right)}\right),\left(n=1,2\dots \right),$
${C}_{n}^{\left(2\right)}=\sum _{m=0}^{n}{B}_{n}^{\left(1\right)}{D}_{n-m}^{\left(1\right)},\left(n=0,1...\right).$

Due to the property of series solution, the solution of Eq. (32) can be written as:

34
$\begin{array}{l}{\phi }_{2}\left(x\right)={\phi }_{m}\left[{\zeta }_{1}\sum _{n=0}^{\infty }\frac{{\omega }^{2n}{x}^{4n+1}}{\left(4n+1\right)!}+{\zeta }_{2}\sum _{n=0}^{\infty }\frac{{\omega }^{2n}{x}^{4n+3}}{\left(4n+3\right)!}]\right\\ +{{\phi }_{m}}^{3}\left(\sum _{n=1}^{\infty }{C}_{n}{x}^{4n-1}+\sum _{n=0}^{\infty }{D}_{n}{x}^{4n+1}+\sum _{n=0}^{\infty }{E}_{n}{x}^{4n+3}\right)\\ +{\phi }_{m}\left(\sum _{n=1}^{\infty }{A}_{n}{x}^{4n-1}+\sum _{n=0}^{\infty }{B}_{n}{x}^{4n+1}\right),\\ \end{array}$

where ${\zeta }_{1}$ and ${\zeta }_{2}$ denote the undetermined coefficients:

${A}_{1}={B}_{0}={C}_{1}={D}_{0}={E}_{0}=0,{A}_{n+1}=\frac{{\omega }^{2}{A}_{n}+\alpha {C}_{n}^{\left(1\right)}}{\left(4n+3\right)\left(4n+2\right)\left(4n+1\right)4n},$
${B}_{n+1}=\frac{{\omega }^{2}{B}_{n}+\alpha {D}_{n}^{\left(1\right)}}{\left(4n+5\right)\left(4n+4\right)\left(4n+3\right)\left(4n+2\right)},{C}_{n+1}=\frac{{\omega }^{2}{C}_{n}+\beta {A}_{n}^{\left(2\right)}}{\left(4n+3\right)\left(4n+2\right)\left(4n+1\right)4n},$
${D}_{n+1}=\frac{{\omega }^{2}{D}_{n}+\beta {B}_{n}^{\left(2\right)}}{\left(4n+5\right)\left(4n+4\right)\left(4n+3\right)\left(4n+2\right)},{E}_{n+1}=\frac{{\omega }^{2}{E}_{n}+\beta {C}_{n}^{\left(2\right)}}{\left(4n+7\right)\left(4n+6\right)\left(4n+5\right)\left(4n+4\right)}.$

By substitution of Eq. (34) into Eqs. (23-24), one can be get:

35
$D\zeta =0,$

where:

$\zeta ={\left[\begin{array}{lll}{\zeta }_{1}& {\zeta }_{2}& 1\end{array}\right]}^{T},D=\left[\begin{array}{lll}{d}_{11}& {d}_{12}& {d}_{13}\\ {d}_{21}& {d}_{22}& {d}_{23}\\ {d}_{31}& {d}_{32}& {d}_{33}\end{array}\right],$

where:

${d}_{11}=\sum _{n=0}^{\infty }\frac{{\omega }^{2n}}{\left(4n+1\right)!},{d}_{12}=\sum _{n=0}^{\infty }\frac{{\omega }^{2n}}{\left(4n+3\right)!}$
${d}_{13}=\left(\sum _{n=1}^{\infty }{A}_{n}+\sum _{n=0}^{\infty }{B}_{n}\right)+{\phi }_{m}^{2}\left(\sum _{n=1}^{\infty }{C}_{n}+\sum _{n=0}^{\infty }{D}_{n}+\sum _{n=0}^{\infty }{E}_{n}\right),$
${d}_{21}=\sum _{n=0}^{\infty }\frac{{\omega }^{2n}}{\left(4n+1\right)!}{\left(\frac{1}{2}\right)}^{4n+1},{d}_{22}=\sum _{n=0}^{\infty }\frac{{\omega }^{2n}}{\left(4n+3\right)!}{\left(\frac{1}{2}\right)}^{4n+3},$
$\begin{array}{l}{d}_{23}=\left(\sum _{n=1}^{\infty }{A}_{n}{\left(\frac{1}{2}\right)}^{4n-1}+\sum _{n=0}^{\infty }{B}_{n}{\left(\frac{1}{2}\right)}^{4n+1}\right)\\ +{{\phi }_{m}}^{2}\left[\sum _{n=1}^{\infty }{C}_{n}{\left(\frac{1}{2}\right)}^{4n-1}+\sum _{n=0}^{\infty }{D}_{n}{\left(\frac{1}{2}\right)}^{4n+1}\right+\sum _{n=0}^{\infty }{E}_{n}{\left(\frac{1}{2}\right)}^{4n+3}]-1,\end{array}$
${d}_{31}=\sum _{n=1}^{\infty }\frac{{\omega }^{2n}}{\left(4n-1\right)!},{d}_{32}=\sum _{n=0}^{\infty }\frac{{\omega }^{2n}}{\left(4n+1\right)!},$
$\begin{array}{l}{d}_{33}=\left(\sum _{n=1}^{\infty }\left(4n-1\right)\left(4n-2\right){A}_{n}+\sum _{n=0}^{\infty }\left(4n+1\right)4n{B}_{n}\right)\\ +{{\phi }_{m}}^{2}\left[\sum _{n=1}^{\infty }\left(4n-1\right)\left(4n-2\right){C}_{n}\right+\sum _{n=0}^{\infty }\left(4n+1\right)4n{D}_{n}+\sum _{n=0}^{\infty }\left(4n+3\right)\left(4n+2\right){E}_{n}].\end{array}$

The frequency ${\omega }_{2}$ can be obtained by $\mathrm{d}\mathrm{e}\mathrm{t}D=$ 0. The coefficients ${\zeta }_{1}$ and ${\zeta }_{2}$ can be determined by Eq. (35), and then, the second-order modified-iterative solution can be obtained.

## 4. Numerical simulation and analysis

To verify the proposed method, some numerical simulations are carried out by GL-E36 strip. The main parameters of tandem mill and strip are listed as follows: the distance between stand F2 and stand F3 $l=$ 2.5 m, the thickness of strip $h=$ 0.018 m, the torsional stiffness $k=$ 9×103 N·m/rad, and the strip tension ${P}_{0}=$ 8×103 N. When the rotational inertia of roller $j$ is 200 kg·m2, 400 kg·m2, 600 kg·m2 and 800 kg·m2, respectively, the relationship curves between axial velocity and frequency of strip are shown in Fig. 4. It can be seen that the vibration frequency decreases with the increasing axial velocity gradually, and the greater the rotational inertia of roller, the faster the frequency decreases. In particular, the vibration frequency becomes more obvious, when ${v}_{0}<$ 0.01 m/s, greater rotational inertia and larger vibration frequency, and vice versa. That is to say, the axial velocity has strong influence on vibration characteristic of strip while the rotational inertia of roller is larger than the others.

Fig. 4Relationship between axial velocity and frequency

Fig. 5Relationship between tension and frequency

Fig. 5 shows the relationship curves between strip tension and frequency of strip with different values of rotational inertias, which are 200 kg·m2, 400 kg·m2 and 600 kg·m2, respectively. The torsional stiffness $k=$ 9×103 N·m/rad, the axial velocity ${v}_{0}=$ 0.008 m/s. From Fig. 5, the vibration frequency decreases with the increase of the strip tension. When ${P}_{0}\le$ 5×103 N, as long as the rotational inertia becomes larger, the frequency of strip gets stronger. Whereas, when ${P}_{0}>$ 5×103 N, the ${P}_{0}-{\omega }_{0}$ curves present fluctuation, then stabilize gradually. Thus, it follows that a smaller strip tension change results in a greater effect of frequency. Additionally, it is shown that the bigger the rotational inertia is, the faster the trend of ${P}_{0}-{\omega }_{0}$ curve declines. In other words, a further proof is given that a bigger variation of frequency caused by a greater rotational inertia.

Fig. 6 depicts the relationship curves between torsional stiffness of drive system and frequency of strip with different values of rotational inertias, which are 200 kg·m2, 400 kg·m2, 600 kg·m2, and 800 kg·m2, respectively. The strip tension ${P}_{0}=$ 8×103 N, the axial velocity ${v}_{0}=$ 0.008 m/s. As can be seen from the Fig. 6, with the increase of the torsional stiffness, the $k-{\omega }_{0}$ curves show a uniform change, and the vibration frequency decreases gradually. And a larger rotational inertia leads to a larger downtrend of $k-{\omega }_{0}$ curve. The result shows that the boundary of strip is close to the clamped condition with the infinite increase of rotational inertia. Compared with a model of moving strip without main drive system, the overall frequency curves tend to be more stable. Moreover, the greater the torsional stiffness is, the smaller the rolling forces exerts, which reduces the vibration frequency of strip to some extent.

From the frequency-response curves in Figs. 4-6, a common feature can be found that a larger rotational inertia of roller gives rise to a greater variation tendency of frequency, which causes a more obvious vibration. In other words, a larger rotational inertia has a greater effect on the relationships from the influencing factors that is the axial velocity, strip tension and torsional stiffness to the frequency. Therefore, it is very important to select the appropriate rotational inertia of roller in the research and practice.

Fig. 7 illustrates the amplitude-frequency responses under the condition of different axial velocities: 2×10–3 m/s, 2×10–3 m/s, 10×10–3 m/s, 14×10–3 m/s and 18×10–3 m/s. The torsional stiffness $k=$ 4×104 N·m/rad, the strip tension ${P}_{0}=$ 8×103 N, and the rotational inertia $j=$ 800 kg·m2. From Fig. 7, it is obtained that with the raise of the axial velocity, the vibration performance is transformed from sclerotic type into intenerate type gradually. When ${v}_{0}\le$ 10×10–3 m/s, with the increase of the amplitude, the vibration frequencies increase, and the hardening degree increases with the decreasing axial velocity. When ${v}_{0}\ge$ 14×10–3 m/s, the curve shows a trending of decrease, namely, the vibration performance is intenerated type. That is, being oversized or excessively small of the axial velocity will cause the great influence of the amplitude on frequency.

Fig. 6Relationship between torsional stiffness and frequency

Fig. 7Amplitude-frequency response with different axial velocities

Fig. 8 demonstrates the amplitude-frequency responses under the condition of different strip tensions: 4×103 N, 6×103 N, 8×103 N and 10×103 N. The torsional stiffness $k=$ 4×104 N·m/rad, the rotational inertia $j=$ 800 kg·m2, and the axial velocity ${v}_{0}=$ 0.01 m/s. In Fig. 8, with the increase of the amplitude, the vibration frequency increases, and the vibration performance keeps consistent. It also can be known that the changing of vibration frequency is small and stable under the condition of lower amplitude. However, at the point of amplitude ${\phi }_{0}=$ 0.006, there is an intersection. After that, with the raise of the amplitude, the vibration frequency increases fast. That is to say, a smaller strip tension results in a larger variation range of the amplitude-frequency response. Based on the mechanics of vibration, the smaller the tensile force, the more unstable the system, thereby, it makes strip vibration quite strong.

Fig. 8Amplitude-frequency response with different strip tensions

Fig. 9Amplitude-frequency response with different torsional stiffness

Fig. 9 displays the amplitude-frequency responses under the condition of different torsional stiffness: 2×104 N·m/rad, 4×104 N·m/rad, 6×104 N·m/rad and 8×104 N·m/rad. The rotational inertia $j=$ 800 kg·m2, the strip tension ${P}_{0}=$ 8×103 N, and the axial velocity ${v}_{0}=$ 0.01 m/s. From Fig. 9, as the torsional stiffness increases, the vibration performance changes from positive correlation to negative correlation gradually. When $k\le$ 4×104 N·m/rad, the frequencies increase with the increasing amplitude, and the vibration performance shows positive correlation. Keep on increasing the size of torsional stiffness, the amplitude can be decreased, and the vibration performance shows negative correlation. Due to the above analysis, unsuitable size of the torsional stiffness will lead to inconsistent movement between rollers and strip, therefore, it leads to the torsional vibration of main drive system. Consequently, the torsional stiffness should be selected appropriately in practical design to reduce the rolling instability.

## 5. Conclusions

In this paper, the strip vibration coupled with torsional vibration of main drive system of rolling mill is investigated. It is given that the influences of the axial velocity, strip tension and torsional stiffness of drive system on vibration frequency of strip, in which all of the influencing factors are related to the rotational inertia of roller. A larger rotational inertia of roller gives rise to a greater changing trend of frequency, and leads to a more obvious vibration of strip. Moreover, the axial velocity, the strip tension and the torsional stiffness of drive system have also great influences on the amplitude-frequency characteristic of strip vibration. The changing trends of the amplitude-frequency response curves are accordantly increasing with the different strip tension, and a larger tension results in a smaller increasing gradient. Whereas, the amplitude-frequency characteristic of strip vibration will transform from sclerotic type into intenerate type with the increasing axial velocity or with the increasing torsional stiffness. Therefore, the analysis results of this paper show that the importance of choosing appropriate axial velocity, strip tension, torsional stiffness and rotational inertia on control and optimization of strip vibration of rolling mill.

#### References

• Hou F. X., Zhang J., Cao J. G. Review of chatter studies in cold rolling. Journal of Iron and Steel Research, Vol. 19, Issue 10, 2007, p. 6-11, (in Chinese).
• Wright J. Mill drive system to minimize torque amplification. Iron and Steel Engineer, Vol. 53, Issue 7, 1976, p. 56-60.
• Fan X. B., Zang Y., Wang F., et al. Hot strip mill nonlinear torsional vibration with multi-stand coupling. Journal of Vibroengineering, Vol. 17, Issue 4, 2015, p. 1623-1633.
• Moller R. H., Hoggart J. S. Periodic surface finish and torque effects during cold strip rolling. Journal of the Australian Institute of Metals, Vol. 12, Issue 2, 1967, p. 155-164.
• Lawrence A. B., Thomas S. R. Winding reel involvement in temper mill chatter. Iron and Steel Engineer, Vol. 71, Issue 11, 1994, p. 27-29.
• Swiatoniowski A., Bar A. Parametrical excitement vibration in tandem mills-mathematical model and its analysis. Journal of Materials Processing Technology, Vol. 134, Issue 2, 2003, p. 214-224.
• Wang R. P., Peng Y., Zhang Y. Mechanism research of rolling mill coupled vibration. Journal of Mechanical Engineering, Vol. 49, Issue 12, 2013, p. 66-71, (in Chinese).
• Sun J. L., Peng Y., Liu H. M. Dynamic characteristics of cold rolling mill and strip based on flatness and thickness control in rolling process. Journal of Central South University, Vol. 21, Issue 2, 2014, p. 567-576.
• Wang Y., Xu F., Li Y. Dynamic modeling and coupling vibration analysis of hybrid systems consisting of strip, rolls and flexible supports. Journal of Vibration Engineering, Vol. 26, Issue 4, 2013, p. 599-607, (in Chinese).
• Xu F., Li Y., Li L. F. Parametric vibration analysis for sheet metal and its supporting system. Mechanical Research and Application, Vol. 26, Issue 3, 2013, p. 27-30, (in Chinese).
• Zou J. X., Xu L. J. Tandem Mill Vibration Control. Metallurgical Industry Press, Beijing, 1998, (in Chinese).
• Gao C. Y., Du G. J., Feng Y., et al. Nonlinear vibration analysis of moving strip with inertial boundary condition. Mathematical Problems in Engineering. Vol. 8, Issue 2015, 2015, p. 1-9.
• Gao C. Y., Du G. J., Li J. X. Numerical analysis on torsional vibration of main driving system of rolling mill with strip steel. Journal of Yanshan University, Vol. 40, Issue 1, 2016, p. 58-65, (in Chinese).
• Chen L. Q., Yang X. D. Steady-state response of axially moving viscoelastic beam with pulsating speed: comparison of two nonliner models. International Journal of Solids and Structures, Vol. 42, 2005, p. 10-15.
• Du G. J., Ma J. Q. Nonlinear vibration and buckling of circular sandwich plate under complex load. Applied Mathematics and Mechanics, Vol. 28, Issue 8, 2007, p. 1081-1091.
• Shen Y. Z., Liu H. M., Xiong J., and et al. Analysis of chaotic behavior of the main drive system with clearance of a heavey plate mill. Engineering Mechanics, Vol. 27, Issue 7, 2010, p. 232-236, (in Chinese).
• Han S. M., Benaroya H., Wei T. Dynamics of transversely vibrating beam using four engineering theories. Journal of Sound and Vibration, Vol. 225, Issue 5, 1999, p. 935-988.
• Huang J. L., Chen S. H. Study on nonlinear vibration of an axially moving beam with coupled transverse and longitudinal motions. Journal of Vibration and Shock, Vol. 30, Issue 8, 2011, p. 24-27, 50, (in Chinese).

#### Cited by

Influence Factors of Inherent Characteristics of F5 Cold-Rolling Mill’s Main Drive System
Impact of shafting angle of F5 cold rolling mill on the dynamics of main drive system
Torsional Vibration for Rolling Mill with the Drive System Shaft Axis Deviations
Advances in Asset Management and Condition Monitoring

31 January 2017
Accepted
26 July 2017
Published
31 December 2017
SUBJECTS
Mechanical vibrations and applications
Keywords
strip vibration
main drive system
torsional vibration
continuous-discrete coupled vibration model
modified iteration method
Acknowledgements

This research was supported in part by Natural Science Foundation of Hebei Province (E2017203115), in part by Key Project of Science and Technology of Hebei Higher School (No. ZD2015077), and in part by Doctor Foundation of Yanshan University (No. B992).