# The Analytic Study on the Heavy Steel Plate Snake Rolling with the Same Roll Diameters

Authors: Lian-yun Jiang; Qing-cheng Meng; Chun-jiang Zhao; Shou-xin Wang; Yan-wei Liu
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

The deformation in the inner region along the thickness of the heavy steel plate can be improved by snake rolling method. Then the microstructure and property will be refined and the crack in the inner region may be avoided. Therefore, the in-depth research on snake rolling method mechanics parameter modeling should be conducted to guide production. A snake rolling process with the same roll diameters and different angular velocity was conducted in this paper. The rolling deformation zone will be divided into back slip zone, front slip zone, cross shear zone, and reverse deflection zone according to the direction of the friction during the snake rolling process. The four zones may not exist at the same time. The boundary conditions of existence of the back slip zone, front slip zone, and cross shear zone were established according to the relationship between threading angle and neutral angle. The calculating models which were used to calculate the snake rolling mechanical parameters including the rolling force and rolling torque were set up. The calculated models of unit compressive pressure in the four zones were set up by the slab method, and on this basis the accurate calculating models of the rolling force and rolling torque were set up according to the composition of the rolling deformation zone and the boundary condition. The mechanical parameters were calculated by the analytical method and the numerical method, and the relative deviation is less than 6% which can satisfy the industrial requirement. The present analytical model can predict the characteristics during snake rolling easily and quickly and it is also suitable for online control applications.


---

## Body

## 1. Introduction

The heavy steel plate with fine performance has been widely used in military equipment, warship, nuclear power, offshore platform, pressure vessel, heavy machinery, and so on. The production capacity of the cast ingot and the opening degree of the mill make it difficult to satisfy the required rolling reduction ratio [1]. Then the heavy steel plate with fine performance in the inner region cannot be easily obtained. Thus, one of the core problems that constrain the heavy steel plate rolling production process is the lack of sufficient deformation produced in the inner region, and then the problem of microstructure and properties or the crack will appear in the inner region.Asymmetrical rolling has been used in cold strip rolling production, and the application results show that the inner grains have been refined compared to the symmetrical rolling. Therefore, the inner region deformation can be improved by the asymmetrical rolling [2–4]. The heavy steel plate curvature problem will appear after asymmetrical rolling and the plate threading in the next pass, or the plate rotation will be difficult to achieve. The plate curvature problem that was produced during asymmetrical rolling for thin strip steel can be solved by the tension between the two stands, but it is not feasible to solve the plate curvature problem that is produced in the heavy steel plate asymmetrical rolling process by this method. The deformation in the inner region can also be improved by the snake rolling process, and the plate curvature after rolling can also be kept well [5, 6], so this problem can be solved by the snake rolling method which is a kind of asymmetrical rolling. The plate curvature and the strain along the thickness calculation or prediction must be conducted on the basis of the mechanical parameters modeling. Therefore, the in-depth research on snake rolling mechanical parameters modeling will be of vital importance for snake rolling industrial application.Currently, some scholars used the slab method, stream function method, and finite element method to study the rolling force and rolling torque modeling during asymmetrical rolling. Hwang Y. M. set up the mechanical parameters calculating models of asymmetrical hot-sheet rolling and asymmetrical clad sheet rolling, and the models assume that the normal stress was well distributed and the shear stress in the cross section was neglected [7, 8]. G.Y. Tzou and X. Liu constructed a model used to calculate the minimum thickness for asymmetrical hot and cold PV rolling which considered the effect of friction factor [9, 10]. H. Gao studied the influence of friction coefficient ratio on the cross shear deformation and rolling pressure and torque during asymmetrical rolling process by the slab method [4]. Tian Y. analyzed the asymmetrical rolling by slab method, and the deformation zone was divided into three zones according to the direction of the friction [11]. The contact arc was replaced by parabola, and the constant surface friction status was adopted during the analysis. M. Salimi provided an analytical model to calculate the strip curvature based on the slab method and modified slab method [12, 13]. F. Afrouz provided an analytical model of asymmetric rolling of unbounded clad sheets which considered the influence of the shear effects [14]. M. Kadkhodaei and S.H. Zhang set up the calculating models of mechanical parameters by the genetic algorithm and slab method [15, 16]. Su-wen Chen and A. Parvizi constructed the mechanical parameters calculating models by the slab method during large cylindrical shell rolling [17, 18].The mechanical parameters calculation models established by the above scholars are all based on the consideration of the deformation zone composed of back slip zone, cross shear zone, and front slip zone. The deformation zone of asymmetrical rolling has a variety of composition forms. The composition form of the deformation zone is different, and the boundary conditions corresponding to it are different. At present, no mechanical parameters calculation model considering the composition form of the deformation zone has been found. Therefore, it is necessary to establish a mechanical parameters calculation model that takes into account the composition form of the deformation zone.The snake rolling process has been used in the thick aluminum plate rolling process and achieved certain progress. S.Y. Li and G.F. Xu studied the microstructure and texture of aluminum alloy plate snake rolling and confirmed that snake rolling can refine the inner region grain of aluminum alloy plate and improve the mechanical properties of the inner region [19, 20]. T. Zhang and J.X. Yang used the numerical method and experimental method to study the bending behavior and strain heterogeneity distribution of thick aluminum alloy plate snake rolling [21, 22]. T. Zhang employed the numerical method and the experimental method to study the strain distribution and temperature distribution of thick aluminum alloy plate snake rolling [23–25].Three zones (back slip zone, front slip zone, and cross shear zone) will be divided in the asymmetrical rolling deformation zone. The distribution of the rolling deformation zone during snake rolling is different from asymmetrical rolling because of the offset distance of work rolls. Consequently, the calculating models that were used to calculate asymmetrical rolling mechanical parameters cannot be used to calculate snake rolling mechanical parameters. And the calculating models that were used to calculate snake rolling mechanical parameters taking into account the composition forms of the deformation zone should be set up in this paper.


## 2. Snake Rolling Basic Parameters Modeling

The necessary assumptions that need to be made in the process of establishing a mathematical model of snake rolling are as follows:(1) The width-to-thickness ratio of steel plate is very large for the snake rolling process. The width spread of the heavy steel plate is relatively small which can be neglected. Hence the heavy steel plate snake rolling model can be simplified as plane strain problem.(2) The deformation of the heavy steel plate is very large compared to the deformation of the work rolls, and the deformation of the work rolls can be neglected. As a result, the work rolls were assumed as a rigid body.(3) The diameter of the top and bottom work rolls are the same, and the angular velocity of the bottom work roll is larger than the angular velocity of the top work roll, which isR1=R2=R; n2>n1.


### 2.1. The Threading Angle Calculation

There is a roll offset between the top and bottom work rolls during snake rolling, and reductions of the top and bottom work rolls are different. Figure1 shows the geometric relationship of the snake rolling. A and B are threading points, which are, respectively, regarded asAC andBD parallel to the rolling direction, and the bottom work roll center which is in the vertical direction intersects at C and D points. From the geometric relationship in Figure 1, it can be obtained as follows:Figure 1
Schematic diagram of snake rolling plastic deformation.for the top work roll:(1)l′=R12-R1-Δh12-dfor the bottom work roll:(2)l′=R22-R2-Δh22along the vertical direction:(3)Δh1+Δh2=H-h0combining (1)~(3):(4)Δh1=ab+R+2abR+R2-a2b2+1(5)Δh2=H-h0-Δh1=H-h0-ab+R+2abR+R2-a2b2+1wherea=2R(H-h0)-(H-h0)2-d2/2d,b=2R-H+h0/d.Total reduction:(6)Δh=Δh1+Δh2=H-h0and, as shown in Figure1, Δh1=R1-R1·cos⁡α1. Hence it is concluded that (7)cos⁡α1=1-Δh1R.In the same way:(8)cos⁡α2=1-Δh2R.Because the threading angleα1 and α2 are very small, α=max⁡α1,α2, α<10o~15o, taking sin⁡α/2≈α/2 and then (9)α1≈2Δh1R(10)α2≈2Δh2R.Length of the deformation zone is as follows.(11)l=R12-R1-Δh12=2R·Δh1-Δh12


### 2.2. The Neutral Angle Calculation

The neutral angles of the top and bottom work rolls are also different because of the existence of the roll offset between the top and bottom work rolls and the difference of the angular velocity between the top and bottom work rolls.The forward slips of the top and bottom work rolls are, respectively, as follows.(12)S1=v0-v1v1S2=v0-v2v2According to the principle of metal equal flow during snake rolling,(13)hγ1vγ1=hγ2vγ2=h0v0.Equations (14) and (15) will be obtained from Figure 1.(14)hγ1=h0+R-Rcos⁡γ1+R-R2-R2sin2γ1+2Rdsin⁡γ1-d2=h0+2R-Rcos⁡γ1-R2cos2γ1+2Rdsin⁡γ1-d2(15)hγ2=h0+R-Rcos⁡γ2+R-R2-R2sin2γ2-2Rdsin⁡γ2-d2=h0+2R-Rcos⁡γ2-R2cos2γ2-2Rdsin⁡γ2-d2Sincevγ1=v1cos⁡γ1, vγ2=v2cos⁡γ2, substituting them into the (13) and (16), we obtain(16)v0v1=hγ1cos⁡γ1h0,v0v2=hγ2cos⁡γ2h0so(17)S1=v0v1-1=h0+2R-Rcos⁡γ1-R2cos2γ1+2Rd·sin⁡γ1-d2cos⁡γ1h0-1.Since2Rd·sin⁡γ1-d2 is much smaller than R2cos2γ1, it can be ignored, so (17) can be simplified as follows.(18)S1=h0+2R-2Rcos⁡γ1cos⁡γ1h0-1Asγ1 is very small in general, so the equation can be further simplified as follows.(19)S1=2R1-cos⁡γ1h0Thus the neutral angle of the top work roll can be obtained as follows.(20)γ1=arccos⁡1-S1h02RSimilarly, the neutral angle of the bottom work roll can be obtained as follows.(21)γ2=arccos⁡1-S2h02R


## 2.1. The Threading Angle Calculation

There is a roll offset between the top and bottom work rolls during snake rolling, and reductions of the top and bottom work rolls are different. Figure1 shows the geometric relationship of the snake rolling. A and B are threading points, which are, respectively, regarded asAC andBD parallel to the rolling direction, and the bottom work roll center which is in the vertical direction intersects at C and D points. From the geometric relationship in Figure 1, it can be obtained as follows:Figure 1
Schematic diagram of snake rolling plastic deformation.for the top work roll:(1)l′=R12-R1-Δh12-dfor the bottom work roll:(2)l′=R22-R2-Δh22along the vertical direction:(3)Δh1+Δh2=H-h0combining (1)~(3):(4)Δh1=ab+R+2abR+R2-a2b2+1(5)Δh2=H-h0-Δh1=H-h0-ab+R+2abR+R2-a2b2+1wherea=2R(H-h0)-(H-h0)2-d2/2d,b=2R-H+h0/d.Total reduction:(6)Δh=Δh1+Δh2=H-h0and, as shown in Figure1, Δh1=R1-R1·cos⁡α1. Hence it is concluded that (7)cos⁡α1=1-Δh1R.In the same way:(8)cos⁡α2=1-Δh2R.Because the threading angleα1 and α2 are very small, α=max⁡α1,α2, α<10o~15o, taking sin⁡α/2≈α/2 and then (9)α1≈2Δh1R(10)α2≈2Δh2R.Length of the deformation zone is as follows.(11)l=R12-R1-Δh12=2R·Δh1-Δh12


## 2.2. The Neutral Angle Calculation

The neutral angles of the top and bottom work rolls are also different because of the existence of the roll offset between the top and bottom work rolls and the difference of the angular velocity between the top and bottom work rolls.The forward slips of the top and bottom work rolls are, respectively, as follows.(12)S1=v0-v1v1S2=v0-v2v2According to the principle of metal equal flow during snake rolling,(13)hγ1vγ1=hγ2vγ2=h0v0.Equations (14) and (15) will be obtained from Figure 1.(14)hγ1=h0+R-Rcos⁡γ1+R-R2-R2sin2γ1+2Rdsin⁡γ1-d2=h0+2R-Rcos⁡γ1-R2cos2γ1+2Rdsin⁡γ1-d2(15)hγ2=h0+R-Rcos⁡γ2+R-R2-R2sin2γ2-2Rdsin⁡γ2-d2=h0+2R-Rcos⁡γ2-R2cos2γ2-2Rdsin⁡γ2-d2Sincevγ1=v1cos⁡γ1, vγ2=v2cos⁡γ2, substituting them into the (13) and (16), we obtain(16)v0v1=hγ1cos⁡γ1h0,v0v2=hγ2cos⁡γ2h0so(17)S1=v0v1-1=h0+2R-Rcos⁡γ1-R2cos2γ1+2Rd·sin⁡γ1-d2cos⁡γ1h0-1.Since2Rd·sin⁡γ1-d2 is much smaller than R2cos2γ1, it can be ignored, so (17) can be simplified as follows.(18)S1=h0+2R-2Rcos⁡γ1cos⁡γ1h0-1Asγ1 is very small in general, so the equation can be further simplified as follows.(19)S1=2R1-cos⁡γ1h0Thus the neutral angle of the top work roll can be obtained as follows.(20)γ1=arccos⁡1-S1h02RSimilarly, the neutral angle of the bottom work roll can be obtained as follows.(21)γ2=arccos⁡1-S2h02R


## 3. The Mechanical Parameters Modeling

Figure1 is a schematic diagram of the snake rolling plastic deformation, in which the origin of the coordinate system is located at the throwing point (point “O” in Figure 1) of the top work roll, and the x axis direction is opposite to the rolling direction. The rolled pieces contact the bottom work roll firstly when the snake rolling process starts. Then the rolled pieces will be lifted, be moved forward, and contact with the top work roll with the action of the friction force of the bottom work roll. The rolled pieces can be threaded into the roll gap under the combined action of top and bottom work rolls. The positions of neutral points of the top and bottom work rolls are also different because of the existence of the roll offset between the top and bottom work rolls and the difference of the angular velocity between the top and bottom work rolls. Thus, the snake rolling plastic deformation zone can be divided into four different zones, as shown in Figure 1: I is back slip zone, II is cross shear zone, III is front slip zone, and IV is reverse deflection zone.As shown in Figure2, take out an element from the four different zones, respectively, in order to calculate the unit compressive pressure. The rolled pieces will contact the top and bottom work rolls at the same time in zone I and zone III. In zone I, as is shown in Figure 2(a), the directions of the friction stresses of rolled pieces top and bottom surfaces are the same with rolling direction, and rolled pieces can be threaded into the roll gap. The direction of friction stresses in zone III is opposite to zone I as is shown in Figure 2(c). In zone II, as is shown in Figure 2(b), the speed of rolled pieces is lower than the surface linear speed of the bottom work roll, but higher than the surface linear speed of the top work roll. Therefore, the direction of friction stress on the top contact surface of the rolled pieces is opposite to the rolling direction, but the direction of friction stress on the bottom contact surface of rolled pieces is the same as the rolling direction. The directions of friction stresses of the top and bottom contact surfaces are opposite in zone II, which makes the metal in the zone have a strong shearing effect and be favorable for the metal plastic deformation to be penetrated into the thickness center. The rolled pieces are affected only by the top work roll in zone IV, as is shown in Figure 2(d), and the speed of rolled pieces is higher than the surface linear speed of the top work roll, so the direction of the friction stress on the top contact surface of rolled pieces is opposite to the rolling direction.Figure 2
Stress states of a vertical element in four zones forv2>v1.
(a) (b) (c) (d)The contact arc between the work rolls and rolled pieces was simplified as parabolic curve during the snake rolling process, so the rolling reduction of the top and bottom work rolls can be obtained at this basis:Δh1(x)=x2/2R, Δh2(x)=(x-d)2/2R. Hence the rolled pieces thickness in the x cross section will be obtained: h=h0+[x2+(x-d)2]/(2R).Equation (22) was obtained according to the equilibrium equation of the force along the horizontal direction for every element in the deformation:(22)dσxh+p1tan⁡θ1+p2tan⁡θ2-τedx=0wheretan⁡θ1≈x/R, tan⁡θ2≈(x-d)/R.


### 3.1. Yield Criterion

The Von-Mises yield criterion at any point is written as below.(23)σx-σy2+σy-σz2+σz-σx2+6τxy2+τyz2+τzx2=2σs2Theτyz=τzx=0 for the plane strain condition, and σz=(σx+σy)/2 according to the flow rule. Substitute the two equations into (23), and (24) will be obtained.(24)32σx-σy2+6τxy2=2σs2The shear stressτxy produced in the surface of the deformation zone will reach the maximum value τs=mk.m is friction factor;k is mean yield shear stress and k=σs/3. Thus the shear yield strength of the materials of rolled pieces can be obtained.(25)τxy=mk=mσs3Substituting (25) into (24), the yield criterion for rolled pieces materials can be obtained.(26)σx-σy=2σs31-m2Consideringσx=q and σy=-p, the following relation was obtained.(27)σx-σy=p+q=M=2σs31-m2


### 3.2. Unit Compressive Pressure Modeling

#### 3.2.1. Zone I (xn1≤x≤l)

The forces in this zone have the following relationships:p 1 = p - τ 1 tan ⁡ θ 1, p2=p-τ2tan⁡θ2, τe=τ1+τ2, and substituting them into (22), the following can be obtained. (28)h0+x2+x-d22R·dpdx=2x-dRM-τ1x2R2-τ2x-d2R2-τ1-τ2=0Through the integration ofx in the above equation, the unit compressive pressure in zone I can be obtained. (29)pI=-2mσs3Rx+Mln⁡2x2-2dx+d2+2Rh0-4mσsR-h03d2+12Rh0arctan⁡2x-dd2+4Rh0+CI


#### 3.2.2. Zone II (xn2≤x≤xn1)

The forces in this zone have the following relationships:p 1 = p + τ 1 tan ⁡ θ 1, p2=p-τ2tan⁡θ2, τe=τ2-τ1, and substituting them into (22), the following can be obtained. (30)h0+x2+x-d22R·dpdx=2x-dRM+τ1x2R2-τ2x-d2R2+τ1-τ2Through the integration ofx in the above equation, the unit compressive pressure in zone II can be obtained. (31)pII=M+σsmd3Rln⁡2Rh0+2x2-2dx+d2+CII


#### 3.2.3. Zone III (d≤x≤xn2)

The forces in this zone have the following relationships:p 1 = p + τ 1 tan ⁡ θ 1, p2=p+τ2tan⁡θ2, τe=-τ1-τ2, and substituting them into (22), the following can be obtained.(32)h0+x2+x-d22R·dpdx=2x-dRM+τ1x2R2+τ2x-d2R2+τ1+τ2=0Through the integration ofx in the above equation, the unit compressive pressure in zone III can be obtained.(33)pIII=2mσs3Rx+Mln⁡2x2-2dx+d2+2Rh0+4mσsR-h03d2+12Rh0arctan⁡2x-dd2+4Rh0+CIII


#### 3.2.4. Zone IV (0≤x≤d)

The forces in this zone have the following relationships:p 1 = p + τ 1 tan ⁡ θ 1, p2=0, τe=-τ1, and substituting them into (22), the following can be obtained.(34)h0+x22Rdpdx=xRM+τ1x2R2+τ1Through the integration ofx in the above equation, the unit compressive pressure in zone IV can be obtained.(35)pIV=2mσs3Rx+Mln⁡x2+2Rh0+2mσsR-2h03Rh0arctan⁡x2Rh0+CIV


### 3.3. The Rolling Force and Rolling Torque Modeling

There is a roll offset between the top and bottom work rolls during snake rolling, and angular velocities of the top and bottom work rolls are different. Consequently, the flow rule and the stress distribution of the mass point of the metal in the deformation zone are different from those of the symmetrical rolling and asymmetrical rolling.The reverse deflection zone (zone IV) always exists because of the existence of the roll offset between the top and bottom work rolls. During snake rolling, the neutral point of the top work roll moves toward the inlet direction, and the neutral point of the bottom work roll moves toward the outlet direction. The directions of the friction force of the top and bottom surfaces of the rolled pieces are opposite when the neutral point of the top work roll moves to the inlet and the neutral point of the bottom work roll moves to the outlet. In this condition, the deformation zone is composed of cross shear zone (zone II) and reverse deflection zone (zone IV). The neutral point cannot move to the inlet or the outlet under some restrictions, and the back slip zone (zone I) or front slip zone (zone III) may appear in the deformation zone. At this time, the deformation zone is composed of back slip zone (zone I), cross shear zone (zone II), and reverse deflection zone (zone IV) or composed of back slip zone (zone I), cross shear zone (zone II), front slip zone (zone III), and reverse deflection zone (zone IV).At the throwing point (point “O” in Figure1) of the top work roll the boundary conditions arex=0,q=0, so pIV=2σs/31-m2 can be obtained. CIV can be obtained by substituting them into (35).(1) The deformation zone is composed of cross shear zone (zone II) and reverse deflection zone (zone IV) whenγ1≥α1, γ2=0. In this case, the boundary conditions in the inlet position arex=l,q=0, and then pII=2σs/31-m2 can be obtained. CII can be obtained by substituting these equations into (31).As a rolled piece is in equilibrium along the vertical direction, the rolling forces of the top and bottom work rolls are equal; namely,F1=F2=F. Therefore, the unit compressive pressure of the top work roll should be used to calculate the rolling force. The rolling force can be obtained by integrating the unit compressive pressure along the contact arc.(36)F=B∫0dpIVdx+∫dlpIIdxThe rolling torque can be obtained by calculating the torque of friction force along the contact arc of top and bottom work rolls.(37)T1=-RB∫0lmkdx=-RBmkl(38)T2=RB∫dlmkdx=RBmkl-d(2) The deformation zone is composed of back slip zone (zone I), cross shear zone (zone II), and reverse deflection zone (zone IV) whenγ1<α1, γ2=0. In this case, the boundary conditions in the inlet position arex=l,q=0, so pI=2σs/31-m2 can be calculated. CI can be obtained by substituting these equations into (29).BecausepII(x=d)=pIV  (x=d) atx=d, CII can be calculated. Atx=xn1,  pI(x=xn1)=pII(x=xn1), soxn1 can be calculated by (29) and (31). The rolling force can be obtained by integrating the unit compressive pressure along the contact arc.(39)F=B∫0dpIVdx+∫dxn1pIIdx+∫xn1lpIdxThe rolling torque can be obtained by calculating the torque of friction force along the contact arc of top and bottom work rolls.(40)T1=RB∫xn1lmkdx-∫0xn1mkdx=RBmkl-2xn1(41)T2=RB∫xn1lmkdx+∫dxn1mkdx=RBmkl-d(3) The deformation zone is composed of back slip zone (zone I), cross shear zone (zone II), front slip zone (zone III), and reverse deflection zone (zone IV) whenγ1<α1,γ2<α2. In this case, the boundary conditions in the inlet position arex=l,q=0, so pI=2σs/31-m2 can be calculated. CI can be obtained by substituting these equations into (29). Because pIII(x=d)=pIV(x=d) atx=d, CIII can be calculated. Atx=xn1,  pI(x=xn1)=pII(x=xn1), so CII(x=xn1) can be calculated. Atx=xn2,  pII(x=xn2)=pIII(x=xn2), so CII(x=xn2) can be calculated. Since the unit compressive pressure is continuous in the neutral points of the top and bottom work rolls, the following equations was obtained as below.(42)CIIx=xn1=CIIx=xn2The volume remains constant during the rolling process, so the following relationship exists.(43)v1cos⁡γ1hxn1=v2cos⁡γ2hxn2Therefore the equation can be obtained as below.(44)v1·1-xn12R2·h0+xn12+xn1-d22R=v2·1-xn2-d2R2·h0+xn22+xn2-d22RBy combining (42) and (44), xn1, xn2, and CII can be obtained.The rolling force can be obtained by integrating the unit compressive pressure along the contact arc.(45)F=B∫0dpIVdx+∫dxn2pIIIdx+∫xn2xn1pIIdx+∫xn1lpIdxThe rolling torque can be obtained by calculating the torque of friction force along the contact arc of top and bottom work rolls.(46)T1=RB∫xn1lmkdx-∫xn2xn1mkdx-∫dxn2mkdx-∫0dmkdx=RBmkl-2xn1(47)T2=RB∫xn1lmkdx+∫xn2xn1mkdx-∫dxn2mkdx=RBmkl+d-2xn2


## 3.1. Yield Criterion

The Von-Mises yield criterion at any point is written as below.(23)σx-σy2+σy-σz2+σz-σx2+6τxy2+τyz2+τzx2=2σs2Theτyz=τzx=0 for the plane strain condition, and σz=(σx+σy)/2 according to the flow rule. Substitute the two equations into (23), and (24) will be obtained.(24)32σx-σy2+6τxy2=2σs2The shear stressτxy produced in the surface of the deformation zone will reach the maximum value τs=mk.m is friction factor;k is mean yield shear stress and k=σs/3. Thus the shear yield strength of the materials of rolled pieces can be obtained.(25)τxy=mk=mσs3Substituting (25) into (24), the yield criterion for rolled pieces materials can be obtained.(26)σx-σy=2σs31-m2Consideringσx=q and σy=-p, the following relation was obtained.(27)σx-σy=p+q=M=2σs31-m2


## 3.2. Unit Compressive Pressure Modeling

### 3.2.1. Zone I (xn1≤x≤l)

The forces in this zone have the following relationships:p 1 = p - τ 1 tan ⁡ θ 1, p2=p-τ2tan⁡θ2, τe=τ1+τ2, and substituting them into (22), the following can be obtained. (28)h0+x2+x-d22R·dpdx=2x-dRM-τ1x2R2-τ2x-d2R2-τ1-τ2=0Through the integration ofx in the above equation, the unit compressive pressure in zone I can be obtained. (29)pI=-2mσs3Rx+Mln⁡2x2-2dx+d2+2Rh0-4mσsR-h03d2+12Rh0arctan⁡2x-dd2+4Rh0+CI


### 3.2.2. Zone II (xn2≤x≤xn1)

The forces in this zone have the following relationships:p 1 = p + τ 1 tan ⁡ θ 1, p2=p-τ2tan⁡θ2, τe=τ2-τ1, and substituting them into (22), the following can be obtained. (30)h0+x2+x-d22R·dpdx=2x-dRM+τ1x2R2-τ2x-d2R2+τ1-τ2Through the integration ofx in the above equation, the unit compressive pressure in zone II can be obtained. (31)pII=M+σsmd3Rln⁡2Rh0+2x2-2dx+d2+CII


### 3.2.3. Zone III (d≤x≤xn2)

The forces in this zone have the following relationships:p 1 = p + τ 1 tan ⁡ θ 1, p2=p+τ2tan⁡θ2, τe=-τ1-τ2, and substituting them into (22), the following can be obtained.(32)h0+x2+x-d22R·dpdx=2x-dRM+τ1x2R2+τ2x-d2R2+τ1+τ2=0Through the integration ofx in the above equation, the unit compressive pressure in zone III can be obtained.(33)pIII=2mσs3Rx+Mln⁡2x2-2dx+d2+2Rh0+4mσsR-h03d2+12Rh0arctan⁡2x-dd2+4Rh0+CIII


### 3.2.4. Zone IV (0≤x≤d)

The forces in this zone have the following relationships:p 1 = p + τ 1 tan ⁡ θ 1, p2=0, τe=-τ1, and substituting them into (22), the following can be obtained.(34)h0+x22Rdpdx=xRM+τ1x2R2+τ1Through the integration ofx in the above equation, the unit compressive pressure in zone IV can be obtained.(35)pIV=2mσs3Rx+Mln⁡x2+2Rh0+2mσsR-2h03Rh0arctan⁡x2Rh0+CIV


## 3.2.1. Zone I (xn1≤x≤l)

The forces in this zone have the following relationships:p 1 = p - τ 1 tan ⁡ θ 1, p2=p-τ2tan⁡θ2, τe=τ1+τ2, and substituting them into (22), the following can be obtained. (28)h0+x2+x-d22R·dpdx=2x-dRM-τ1x2R2-τ2x-d2R2-τ1-τ2=0Through the integration ofx in the above equation, the unit compressive pressure in zone I can be obtained. (29)pI=-2mσs3Rx+Mln⁡2x2-2dx+d2+2Rh0-4mσsR-h03d2+12Rh0arctan⁡2x-dd2+4Rh0+CI


## 3.2.2. Zone II (xn2≤x≤xn1)

The forces in this zone have the following relationships:p 1 = p + τ 1 tan ⁡ θ 1, p2=p-τ2tan⁡θ2, τe=τ2-τ1, and substituting them into (22), the following can be obtained. (30)h0+x2+x-d22R·dpdx=2x-dRM+τ1x2R2-τ2x-d2R2+τ1-τ2Through the integration ofx in the above equation, the unit compressive pressure in zone II can be obtained. (31)pII=M+σsmd3Rln⁡2Rh0+2x2-2dx+d2+CII


## 3.2.3. Zone III (d≤x≤xn2)

The forces in this zone have the following relationships:p 1 = p + τ 1 tan ⁡ θ 1, p2=p+τ2tan⁡θ2, τe=-τ1-τ2, and substituting them into (22), the following can be obtained.(32)h0+x2+x-d22R·dpdx=2x-dRM+τ1x2R2+τ2x-d2R2+τ1+τ2=0Through the integration ofx in the above equation, the unit compressive pressure in zone III can be obtained.(33)pIII=2mσs3Rx+Mln⁡2x2-2dx+d2+2Rh0+4mσsR-h03d2+12Rh0arctan⁡2x-dd2+4Rh0+CIII


## 3.2.4. Zone IV (0≤x≤d)

The forces in this zone have the following relationships:p 1 = p + τ 1 tan ⁡ θ 1, p2=0, τe=-τ1, and substituting them into (22), the following can be obtained.(34)h0+x22Rdpdx=xRM+τ1x2R2+τ1Through the integration ofx in the above equation, the unit compressive pressure in zone IV can be obtained.(35)pIV=2mσs3Rx+Mln⁡x2+2Rh0+2mσsR-2h03Rh0arctan⁡x2Rh0+CIV


## 3.3. The Rolling Force and Rolling Torque Modeling

There is a roll offset between the top and bottom work rolls during snake rolling, and angular velocities of the top and bottom work rolls are different. Consequently, the flow rule and the stress distribution of the mass point of the metal in the deformation zone are different from those of the symmetrical rolling and asymmetrical rolling.The reverse deflection zone (zone IV) always exists because of the existence of the roll offset between the top and bottom work rolls. During snake rolling, the neutral point of the top work roll moves toward the inlet direction, and the neutral point of the bottom work roll moves toward the outlet direction. The directions of the friction force of the top and bottom surfaces of the rolled pieces are opposite when the neutral point of the top work roll moves to the inlet and the neutral point of the bottom work roll moves to the outlet. In this condition, the deformation zone is composed of cross shear zone (zone II) and reverse deflection zone (zone IV). The neutral point cannot move to the inlet or the outlet under some restrictions, and the back slip zone (zone I) or front slip zone (zone III) may appear in the deformation zone. At this time, the deformation zone is composed of back slip zone (zone I), cross shear zone (zone II), and reverse deflection zone (zone IV) or composed of back slip zone (zone I), cross shear zone (zone II), front slip zone (zone III), and reverse deflection zone (zone IV).At the throwing point (point “O” in Figure1) of the top work roll the boundary conditions arex=0,q=0, so pIV=2σs/31-m2 can be obtained. CIV can be obtained by substituting them into (35).(1) The deformation zone is composed of cross shear zone (zone II) and reverse deflection zone (zone IV) whenγ1≥α1, γ2=0. In this case, the boundary conditions in the inlet position arex=l,q=0, and then pII=2σs/31-m2 can be obtained. CII can be obtained by substituting these equations into (31).As a rolled piece is in equilibrium along the vertical direction, the rolling forces of the top and bottom work rolls are equal; namely,F1=F2=F. Therefore, the unit compressive pressure of the top work roll should be used to calculate the rolling force. The rolling force can be obtained by integrating the unit compressive pressure along the contact arc.(36)F=B∫0dpIVdx+∫dlpIIdxThe rolling torque can be obtained by calculating the torque of friction force along the contact arc of top and bottom work rolls.(37)T1=-RB∫0lmkdx=-RBmkl(38)T2=RB∫dlmkdx=RBmkl-d(2) The deformation zone is composed of back slip zone (zone I), cross shear zone (zone II), and reverse deflection zone (zone IV) whenγ1<α1, γ2=0. In this case, the boundary conditions in the inlet position arex=l,q=0, so pI=2σs/31-m2 can be calculated. CI can be obtained by substituting these equations into (29).BecausepII(x=d)=pIV  (x=d) atx=d, CII can be calculated. Atx=xn1,  pI(x=xn1)=pII(x=xn1), soxn1 can be calculated by (29) and (31). The rolling force can be obtained by integrating the unit compressive pressure along the contact arc.(39)F=B∫0dpIVdx+∫dxn1pIIdx+∫xn1lpIdxThe rolling torque can be obtained by calculating the torque of friction force along the contact arc of top and bottom work rolls.(40)T1=RB∫xn1lmkdx-∫0xn1mkdx=RBmkl-2xn1(41)T2=RB∫xn1lmkdx+∫dxn1mkdx=RBmkl-d(3) The deformation zone is composed of back slip zone (zone I), cross shear zone (zone II), front slip zone (zone III), and reverse deflection zone (zone IV) whenγ1<α1,γ2<α2. In this case, the boundary conditions in the inlet position arex=l,q=0, so pI=2σs/31-m2 can be calculated. CI can be obtained by substituting these equations into (29). Because pIII(x=d)=pIV(x=d) atx=d, CIII can be calculated. Atx=xn1,  pI(x=xn1)=pII(x=xn1), so CII(x=xn1) can be calculated. Atx=xn2,  pII(x=xn2)=pIII(x=xn2), so CII(x=xn2) can be calculated. Since the unit compressive pressure is continuous in the neutral points of the top and bottom work rolls, the following equations was obtained as below.(42)CIIx=xn1=CIIx=xn2The volume remains constant during the rolling process, so the following relationship exists.(43)v1cos⁡γ1hxn1=v2cos⁡γ2hxn2Therefore the equation can be obtained as below.(44)v1·1-xn12R2·h0+xn12+xn1-d22R=v2·1-xn2-d2R2·h0+xn22+xn2-d22RBy combining (42) and (44), xn1, xn2, and CII can be obtained.The rolling force can be obtained by integrating the unit compressive pressure along the contact arc.(45)F=B∫0dpIVdx+∫dxn2pIIIdx+∫xn2xn1pIIdx+∫xn1lpIdxThe rolling torque can be obtained by calculating the torque of friction force along the contact arc of top and bottom work rolls.(46)T1=RB∫xn1lmkdx-∫xn2xn1mkdx-∫dxn2mkdx-∫0dmkdx=RBmkl-2xn1(47)T2=RB∫xn1lmkdx+∫xn2xn1mkdx-∫dxn2mkdx=RBmkl+d-2xn2


## 4. Results and Discussions

The numerical method has been used in the rolling process simulation by many scholars, and the calculating precision has been verified by many scholars. As a new rolling technology, the heavy steel plate snake rolling experiment is difficult to be conducted in the laboratory because of the equipment casts. In this case, the numerical method will be used to simulate the snake rolling process, and the rolling force can be obtained by the numerical method. The calculation precision of the analytical model can also be verified.


### 4.1. Numerical Modeling

A finite element method (Arbitrary Lagrangian-Eulerian, ALE) can be used to simulate the rolling process. The snake rolling process can be simulated by Ansys LS-DYNA (explicit dynamic analysis), and the element type selected is a 3D solid element (ID 164). The top and bottom work rolls are modeled as rigid body. The rolls can only rotate along their axis and the other rotational and translational directions are constrained. The heavy steel plate is modeled as a bilinear isotropic material and constrained by the contact between the work rolls and the heavy steel plate. Table1 shows the initial data used in simulation and Figure 3 shows the meshing model used in the FEM (Finite Element Method) simulation. There are 253510 nodes and 196000 elements in the meshing model.Table 1
Data used in simulation.
Items Name Value Work roll diameter/mm 1051 Mill parameters Elastic modulus/GPa 206 Poisson ratio 0.25 Width/mm 3000 Thickness/mm 250 Length/mm 2000 Rolled piece parameters (performance) Rolling reduction/mm 32 Yield stress/MPa 100 Tangent modulus/MPa 10 Elastic modulus/GPa 110 Poisson ratio 0.29 Roll speed/m·s−1 1.3 Other setting parameters Dynamic friction coefficient 0.35 Static friction coefficient 0.4Figure 3
Meshing model used in the snake rolling process simulation.


### 4.2. Analytical and Numerical Results Verifying and Discussions

Thirteen kinds of parameters as shown in Table2 were selected for analytical and numerical calculation in order to verify the analytical results of the theoretical models setup in this paper. One of the time-snake rolling force curves was shown in Figure 4. Then the average rolling force can be obtained as listed in Table 2.Table 2
The relative error of analytic and numerical calculations of rolling force.
No. v 1  (m/s) v 2  (m/s) Reduction(mm) Roll offset (mm) Roll radius(mm) plate thickness(mm) plate width(mm) Rolling force(KN) relative deviation composition of the deformation zone Analytical results Numerical results 1 1.3 1.300 32 10 525.5 250 3000 48152 48727 1.18% I II III IV 2 1.3 1.305 32 10 525.5 250 3000 48120 48547 0.88% I II III IV 3 1.3 1.310 32 10 525.5 250 3000 48045 48389 0.71% I II III IV 4 1.3 1.325 32 10 525.5 250 3000 47949 48065 0.24% I II IV 5 1.3 1.325 32 15 525.5 250 3000 48863 48236 1.30% I II IV 6 1.3 1.325 32 18 525.5 250 3000 49401 48382 2.10% I II IV 7 1.3 1.325 32 20 525.5 250 3000 49756 48592 2.40% I II IV 8 1.3 1.325 32 22 525.5 250 3000 50107 49011 2.24% I II IV 9 1.3 1.45 32 10 525.5 250 3000 43857 46575 5.84% II IV 10 1.3 1.45 32 15 525.5 250 3000 44849 47680 5.94% II IV 11 1.3 1.45 32 18 525.5 250 3000 45442 48309 5.93% II IV 12 1.3 1.45 32 20 525.5 250 3000 45836 48369 5.24% II IV 13 1.3 1.45 32 22 525.5 250 3000 46229 48516 4.71% II IVFigure 4
The time-rolling force curve whenv2=1.325m/s and roll offset is 10 mm.The comparison of the relative deviation between the analytical results and the numerical results was shown in Table2. The composition of the deformation zone with different process parameters was also given in Table 2. It can be seen that the roll speed ratio is the most critical process parameter that affects the composition of the deformation zone.For the rolling parameters shown in Table1, the deformation zone is composed of back slip zone (zone I), cross shear zone (zone II), front slip zone (zone III), and reverse deflection zone (zone IV) when the roll speed ratio is less than 1.008. And the deformation zone is composed of back slip zone (zone I), cross shear zone (zone II), and reverse deflection zone (zone IV) when the roll speed ratio is greater than 1.008. Besides, the deformation zone is composed of cross shear zone (zone II) and reverse deflection zone (zone IV) when the roll speed ratio reaches 1.115.The longer the length of the cross shear zone is, the more beneficial it is to the inner region deformation of the heavy steel plate. Moreover, it can refine inner region grains and reduce the rolling force. And the roll speed ratio can also affect the plate curvature after rolling. Therefore, appropriate roll speed ratio should be confirmed according to the plate curvature requirement and the plate deformation permeability requirement.The curve of roll speed ratios-the rolling force is shown in Figure5. Among them, figure (a) is a full-scale curve, and figure (b) is a partial enlarged view of figure (a). A good agreement is observed between analytic results and numerical results. The figure shows that the rolling force gradually decreases as the roll speed ratio increases. This is because as the roll speed ratio increases, the neutral pointxn1 of the top work roll moves toward the inlet side of the deformation zone, and the neutral pointxn2 of the bottom work roll moves toward the outlet side of the deformation zone. As a result, the length of the cross shear zone increases, resulting in the shearing effect on the rolled pieces strengthen, so the rolling force is reduced.Figure 5
Curve of roll speed ratio-the rolling force.
(a) (b)The curve of roll offset-the rolling force is shown in Figure6. Among them, figure (a) is a full-scale curve, and figure (b) is a partial enlarged view of figure (a). A good agreement is observed between analytic results and numerical results. The figure shows that the rolling force gradually increases as the roll offset increases. This is because increasing the roll offset is equivalent to increasing the length of the reverse deflection zone; thereby the length of the deformation zone is increased.Figure 6
Curve of roll offset-the rolling force.
(a) (b)The curve of roll speed ratio-the rolling torque is shown in Figure7 when the deformation zone is composed of back slip zone (zone I), cross shear zone (zone II), front slip zone (zone III), and reverse deflection zone (zone IV). As can be seen from Figure 7, the rolling torque of the top work roll gradually decreases and the rolling torque of the bottom work roll gradually increases when the roll speed ratio increases. This is because the top work roll neutral pointxn1 moves toward the inlet of the deformation zone and the bottom work roll neutral pointxn2 moves toward the outlet of the deformation zone when the roll speed ratio increases. It can be seen from (46) and (47) that the rolling torque of the top work roll decreases as the roll speed ratio increases, while the rolling torque of the bottom work roll increases as the roll speed ratio increases.Figure 7
Curve of roll speed ratio-the rolling torque.The curve of roll offset-the rolling torque is shown in Figure8 when the deformation zone is composed of cross shear zone (zone II) and reverse deflection zone (zone IV). The figure shows that the rolling torques of the top and bottom work rolls are in the opposite direction. The rolling torque of the top work roll increases with the increase of the roll offset, while the rolling torque of the bottom work roll decreases with the increase of the roll offset. This is because increasing the roll offset is equivalent to increasing the length of the reverse deflection zone and thereby increasing the length of the deformation zone. It can be seen from (37) that the rolling torque of the top work roll increases as the roll offset increases. Although the increase of the roll offset d is equivalent to increasing the length of the reverse deflection zone, the difference between the length of the deformation zone and the length of the reverse deflection zonel-d is relatively reduced. It can be seen from (38) that the rolling torque of the bottom work roll decreases with the increase of the roll offset.Figure 8
Curve of roll offset-the rolling torque.


## 4.1. Numerical Modeling

A finite element method (Arbitrary Lagrangian-Eulerian, ALE) can be used to simulate the rolling process. The snake rolling process can be simulated by Ansys LS-DYNA (explicit dynamic analysis), and the element type selected is a 3D solid element (ID 164). The top and bottom work rolls are modeled as rigid body. The rolls can only rotate along their axis and the other rotational and translational directions are constrained. The heavy steel plate is modeled as a bilinear isotropic material and constrained by the contact between the work rolls and the heavy steel plate. Table1 shows the initial data used in simulation and Figure 3 shows the meshing model used in the FEM (Finite Element Method) simulation. There are 253510 nodes and 196000 elements in the meshing model.Table 1
Data used in simulation.
Items Name Value Work roll diameter/mm 1051 Mill parameters Elastic modulus/GPa 206 Poisson ratio 0.25 Width/mm 3000 Thickness/mm 250 Length/mm 2000 Rolled piece parameters (performance) Rolling reduction/mm 32 Yield stress/MPa 100 Tangent modulus/MPa 10 Elastic modulus/GPa 110 Poisson ratio 0.29 Roll speed/m·s−1 1.3 Other setting parameters Dynamic friction coefficient 0.35 Static friction coefficient 0.4Figure 3
Meshing model used in the snake rolling process simulation.


## 4.2. Analytical and Numerical Results Verifying and Discussions

Thirteen kinds of parameters as shown in Table2 were selected for analytical and numerical calculation in order to verify the analytical results of the theoretical models setup in this paper. One of the time-snake rolling force curves was shown in Figure 4. Then the average rolling force can be obtained as listed in Table 2.Table 2
The relative error of analytic and numerical calculations of rolling force.
No. v 1  (m/s) v 2  (m/s) Reduction(mm) Roll offset (mm) Roll radius(mm) plate thickness(mm) plate width(mm) Rolling force(KN) relative deviation composition of the deformation zone Analytical results Numerical results 1 1.3 1.300 32 10 525.5 250 3000 48152 48727 1.18% I II III IV 2 1.3 1.305 32 10 525.5 250 3000 48120 48547 0.88% I II III IV 3 1.3 1.310 32 10 525.5 250 3000 48045 48389 0.71% I II III IV 4 1.3 1.325 32 10 525.5 250 3000 47949 48065 0.24% I II IV 5 1.3 1.325 32 15 525.5 250 3000 48863 48236 1.30% I II IV 6 1.3 1.325 32 18 525.5 250 3000 49401 48382 2.10% I II IV 7 1.3 1.325 32 20 525.5 250 3000 49756 48592 2.40% I II IV 8 1.3 1.325 32 22 525.5 250 3000 50107 49011 2.24% I II IV 9 1.3 1.45 32 10 525.5 250 3000 43857 46575 5.84% II IV 10 1.3 1.45 32 15 525.5 250 3000 44849 47680 5.94% II IV 11 1.3 1.45 32 18 525.5 250 3000 45442 48309 5.93% II IV 12 1.3 1.45 32 20 525.5 250 3000 45836 48369 5.24% II IV 13 1.3 1.45 32 22 525.5 250 3000 46229 48516 4.71% II IVFigure 4
The time-rolling force curve whenv2=1.325m/s and roll offset is 10 mm.The comparison of the relative deviation between the analytical results and the numerical results was shown in Table2. The composition of the deformation zone with different process parameters was also given in Table 2. It can be seen that the roll speed ratio is the most critical process parameter that affects the composition of the deformation zone.For the rolling parameters shown in Table1, the deformation zone is composed of back slip zone (zone I), cross shear zone (zone II), front slip zone (zone III), and reverse deflection zone (zone IV) when the roll speed ratio is less than 1.008. And the deformation zone is composed of back slip zone (zone I), cross shear zone (zone II), and reverse deflection zone (zone IV) when the roll speed ratio is greater than 1.008. Besides, the deformation zone is composed of cross shear zone (zone II) and reverse deflection zone (zone IV) when the roll speed ratio reaches 1.115.The longer the length of the cross shear zone is, the more beneficial it is to the inner region deformation of the heavy steel plate. Moreover, it can refine inner region grains and reduce the rolling force. And the roll speed ratio can also affect the plate curvature after rolling. Therefore, appropriate roll speed ratio should be confirmed according to the plate curvature requirement and the plate deformation permeability requirement.The curve of roll speed ratios-the rolling force is shown in Figure5. Among them, figure (a) is a full-scale curve, and figure (b) is a partial enlarged view of figure (a). A good agreement is observed between analytic results and numerical results. The figure shows that the rolling force gradually decreases as the roll speed ratio increases. This is because as the roll speed ratio increases, the neutral pointxn1 of the top work roll moves toward the inlet side of the deformation zone, and the neutral pointxn2 of the bottom work roll moves toward the outlet side of the deformation zone. As a result, the length of the cross shear zone increases, resulting in the shearing effect on the rolled pieces strengthen, so the rolling force is reduced.Figure 5
Curve of roll speed ratio-the rolling force.
(a) (b)The curve of roll offset-the rolling force is shown in Figure6. Among them, figure (a) is a full-scale curve, and figure (b) is a partial enlarged view of figure (a). A good agreement is observed between analytic results and numerical results. The figure shows that the rolling force gradually increases as the roll offset increases. This is because increasing the roll offset is equivalent to increasing the length of the reverse deflection zone; thereby the length of the deformation zone is increased.Figure 6
Curve of roll offset-the rolling force.
(a) (b)The curve of roll speed ratio-the rolling torque is shown in Figure7 when the deformation zone is composed of back slip zone (zone I), cross shear zone (zone II), front slip zone (zone III), and reverse deflection zone (zone IV). As can be seen from Figure 7, the rolling torque of the top work roll gradually decreases and the rolling torque of the bottom work roll gradually increases when the roll speed ratio increases. This is because the top work roll neutral pointxn1 moves toward the inlet of the deformation zone and the bottom work roll neutral pointxn2 moves toward the outlet of the deformation zone when the roll speed ratio increases. It can be seen from (46) and (47) that the rolling torque of the top work roll decreases as the roll speed ratio increases, while the rolling torque of the bottom work roll increases as the roll speed ratio increases.Figure 7
Curve of roll speed ratio-the rolling torque.The curve of roll offset-the rolling torque is shown in Figure8 when the deformation zone is composed of cross shear zone (zone II) and reverse deflection zone (zone IV). The figure shows that the rolling torques of the top and bottom work rolls are in the opposite direction. The rolling torque of the top work roll increases with the increase of the roll offset, while the rolling torque of the bottom work roll decreases with the increase of the roll offset. This is because increasing the roll offset is equivalent to increasing the length of the reverse deflection zone and thereby increasing the length of the deformation zone. It can be seen from (37) that the rolling torque of the top work roll increases as the roll offset increases. Although the increase of the roll offset d is equivalent to increasing the length of the reverse deflection zone, the difference between the length of the deformation zone and the length of the reverse deflection zonel-d is relatively reduced. It can be seen from (38) that the rolling torque of the bottom work roll decreases with the increase of the roll offset.Figure 8
Curve of roll offset-the rolling torque.


## 5. Conclusions

(1) The rolling deformation zone will be divided into back slip zone, front slip zone, cross shear zone, and reverse deflection zone according to the direction of the friction during snake rolling process. The calculating models that used to calculate mechanical parameters were set up on this basis: the calculated models of unit compressive pressure in every zone were set up by the slab method; the accurate calculating models of the rolling force and rolling torque were set up according to the composition of the rolling deformation zone and the boundary condition. The mechanical parameters were calculated by the analytical method and the numerical method, and the relative deviation between the analytical results and numerical results was less than 6% which can satisfy the industrial requirement.(2) The reverse deflection zone always exists because of the existence of the roll offset between the top and bottom work rolls. The existence boundary conditions of the back slip zone, front slip zone, and cross shear zone was established according to the relationship between threading angle and neutral angle.(3) The rolling force gradually decreases as the roll speed ratio increases when the deformation zone is composed of a back slip zone (zone I), cross shear zone (zone II), front slip zone (zone III), and reverse deflection zone (zone IV). The rolling torque of the top work roll gradually decreases and the rolling torque of the bottom work roll gradually increases as the roll speed ratio increases.(4) The rolling torques of the top and bottom work rolls are in opposite directions when the deformation zone is composed of cross shear zone (zone II) and a reverse deflection zone (zone IV). The rolling torque of the top work roll increases, and the rolling torque of the bottom work roll decreases when the roll offset is increased.



---

*Source: 7807580-2018-11-07.xml*