# A Novel Modeling Method in Metal Strip Leveling Based on a Roll-Strip Unit

Authors: Guodong Yi; Zili Wang; Zhouyu Hu
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

In the metal strip-multiroll leveling process, the action behavior of each roll is different. However, modeling each roll individually will result in redundancy which is not conducive to the modeling of the entire leveling process. To overcome this problem, a roll-strip unit (RSU) model is proposed to uniformly describe the behavior of each roll during the leveling process. The RSU and its equivalent model are defined on the basis of analyzing the force relationship between the roll and the strip. According to the linear distribution of the bending moment in the longitudinal direction of the strip, the position of the zero bending moment, that is, the virtual fulcrum, is obtained to determine the interval of the RSU. A plastic deformation function is established to describe the influence of plastic extension on the tension and velocity of the strip. The fitting of the deformation curve of the strip is optimized by the tension influence factor and the zero curvature moment. The static friction condition between the roll and the strip which ensures the normal operation of the RSU is given. The AMESim model of the RSU is established to lay the foundation for the dynamic modeling of the multiroll leveler.


---

## Body

## 1. Introduction

Precision metal strips are important materials with good mechanical properties and surface qualities in industrial production and are widely used in the manufacture of products in the aerospace, automotive, medical, electronics, and other industries.Cold rolling and hot rolling are the main processing methods for metal strips and often cause flatness defects due to uneven internal stresses. As an effective tool to eliminate the flatness defects of metal strips, the multiroll leveler uses a set of rolls to bend the strips with different loads to achieve continuous leveling.Currently, research on the leveling (straightening) methods for different types of metal profiles has focused on two aspects. One is the calculation of the leveling force parameters of the profile by the finite element method for designing the leveler. The other is the quality prediction and residual stress analysis of the leveled strips to obtain the optimum leveling parameters.The rectilinear structure has the characteristics of a large length-diameter aspect ratio, so it has similar characteristics of the straightening process. In bar straightening, Wang et al. established a mathematical displacement-force model of the tension straightening process based on the elastic-plastic mechanics’ theory and calculated the predicted displacement of tension straightening for various original deflections [1, 2]. Yu et al. established a quantitative analysis model and a finite element model of the two-roller straightening process and revealed the process mechanism of two-roller straightening [3]. Zhang et al. firstly proposed an asymmetrical hardening material model to describe the stress-strain curves considering the asymmetrical features of yield stress and plastic-hardening stage in tension and compression. The strain and stress distributions in elastic and plastic regions during the straightening process were discussed in their article [4]. Jin et al. researched the rotation blocking mechanism of the triple-roller equal curvature rotation blocking straightening system to solve the problem of bar rotation during the straightening process [5, 6]. In pipe straightening, Zhang et al. presented a simplified cantilever unit model and a mechanical model for calculating the straightening intermesh based on the structural features of the thin-walled tube and the numerical solution of the finite element analysis (FEA) method [7, 8]. Zhang et al. studied the cross-sectional ovalization of the thin-walled circular steel tube during the straightening process, deduced simplified normal strain component formulas based on the thin shell theory, and presented a rational model for predicting the maximum section flattening of the thin-walled circular steel tube under its straightening process by the principle of minimum potential energy [9, 10]. In beam straightening, Yin et al. analyzed the multiroller straightening process of section steel based on the springback theory of small curvature plane bending and redefined the large deformation straightening strategy and the small deformation straightening strategy to improve the straightening quality [11]. Kaiser et al. established a beam with a rectangular cross section to consider the development of curvature during the straightening process and designed a trend of curvature intending to tailor the final residual stress distribution to the desired optimum [12]. Wang et al. proposed the load-deflection pressure straightening model to study the variation of the guide rail bending formation in the contiguous steps of the multistep straightening process through FEM and regression analysis [13].In the plate leveling process, to obtain a flat product, the bending fibers or sections must be made straight from bending, and the residual stresses produced in the process of controlling material deformation must be eliminated [14, 15]. Liu et al. presented a new mechanics model based on the curvature integration method to predict the leveling process of plates and built an optimization model with equality and inequality constraints for the maximum yield stress search of each thickness of plates to judge whether incoming plates can be leveled quickly or not [16–18]. Cui et al. investigated the deformation characteristics and residual curvature during the longitudinal profile plate leveling process based on the curvature integral by elastic-plastic differences to solve the problem of changes in plate thickness [19]. Behrens et al. analyzed the edge and center waves of the sheet metal before and after deformation, investigated the effect of the leveling process on the sheet metal, and calculated the remaining shape defects after leveling [20]. Chen et al. investigated the essential deformation law of leveling for plates with transverse waves and proposed a 2.5-dimensional analytical model related to different initial unevenness to obtain the appropriate values of roller bending [21]. Brauneis et al. presented a fast and robust mathematical model for the leveling process of steel plates to calculate the leveling forces for a given plate and a given roll configuration of the leveler [22]. Among all the plate leveling processes, the heavy and composite plate leveling process is the most complex [23, 24]. Fan et al. established a stress neutral layer offset model based on three-point bending theory, steel plate bending characteristics, and layered theory in the straightening process to obtain the relationship between the neutral layer offset values and the reverse bending radius [25]. Baumgart et al. presented a detailed mathematical deflection model for heavy plate leveling to compensate for the effect of deflection on the roll intermesh and the plate flatness as well as to assess loads of critical parts [26].The high yield strength of the metal strip material is high, which results in high tensile strength. It is hard for the straightening process. Zhang et al. researched the elastic-plastic deformation of a wide strip in tension leveling with the FEM, discovered the rules of tension leveling and the relations between elongation and technique factors, and obtained conclusions about tension leveling technology and tension leveler behaviors [27]. Maksimov et al. presented the tractive force exerted on strips with different mechanical properties by a straightening machine and determined how the tractive forces are affected by the cross-sectional area of the strip and the unit back tension [28]. Stadler et al. investigated the curvature and the contact force at the strip-roll contact point of a metal strip in an experimental device and derived a material model and a steady-state strip deformation model to compute the feasible combinations of curvature and contact force [29]. Grüber and Hirt proposed a control strategy on a force measurement in the first load triangle of a leveling machine and investigated a feedforward control strategy based on control curves for varying yield strengths along the strip length [30, 31]. Nikula et al. analyzed the mechanical stress inflicted on a roller leveler processing cold steel strips based on a feature extracted from vibration measurements and introduced models specific to different steel grades to the prediction of the stress level [32, 33].The research mentioned above has made great progress in the leveling (or straightening) methods for different metal profiles. However, the action behavior of each roll on the strip is different in the multiroll leveling process, and modeling each roll individually will result in redundancy which is not conducive to the modeling of the entire leveling process. To describe the multiroll leveling process in a unified way, we propose an RSU model. To verify the effectiveness of the proposed model, we establish the RSU model in AMESim.The remainder of this paper is organized as follows: The definition of the RSU based on the force analysis is firstly proposed in Section2. Then, the virtual fulcrums of the RSU and the friction between the roll and the strip are analyzed in Section 3. In Section 4, the deformation curve is fitted and the deformation function of the strip is constructed. To verify the effectiveness of the proposed method, the RSU is modeled in AMESim and applied in the rolling production line in Section 5. At last, we draw conclusions in Section 6.


## 2. Definition of the RSU Based on Force Analysis

The multiroll leveling of strips can be abstracted as the interaction between the roll and the strip so that the roll and the strip section it acts on can be used as a basic component of the leveler, which is a roll-strip unit (RSU). To establish the RSU model, the force relationship between the roll and the strip during the leveling process is first analyzed, as shown in Figure1.Figure 1
Force relationship between the roll and the strip.In Figure1, the strip applies a pressure Pi to the roll i and produces a torque Ri due to the friction between them. The strip is mainly affected by tension and bending moment, without considering the effect of shear force. Within the length element, the force is shown in Figure 2.Figure 2
Force analysis of the strip element.Because the shear force is not considered, that is,∑Q=0, dQ=0 and Q′=dQ/dx=0. According to Figure 2, the bending moment of the strip element balanced, that is, ∑M=0, so Q=dM/dx. Take the derivatives on both sides of the equation, then M″x=0; that is to say, the bending moment is linearly distributed along the length in the strip. If Mi−1, Mi, and Mi+1 are the maximum bending moments acting on the strips corresponding to the rolls i−1, i, and i+1, respectively (the directions of Mi−1 and Mi−1 are opposite to that of Mi), the bending moments at the points A and B on the strip are zero, which are named virtual fulcrums. Ti−1 and Ti+1 are the forward and backward tensions of the strip segment, respectively. Fi is the static friction that the roll i acts on the strip.Based on the above-mentioned analysis, an equivalent model of the RSU is established, as shown in Figure3.Figure 3
Equivalent model of the RSU.In the leveling process, the roll is regarded as a rigid body with a fixed rotation center due to its extremely small deformation, which can be equivalent to a set of moments of inertia and damping. The strip can be divided into elastic deformation sections and plastic deformation sections, and each of them is equivalent to a spring-damping system. Seven nodes are set between these sections, and velocities and tensions at these nodes are defined asv1 to v7 and T1 to T7, respectively. The length change of the strip caused by plastic deformation is represented by ξ, which is a function of the internal stress, the elastic modulus, and the speed difference between the two ends of the plastic deformation section. If there is no plastic deformation in the leveling, the length of the plastic deformation section lP is 0. The drive unit of the roll is equivalent to a second-order system that applies a speed or torque to the roll based on a control signal. αi is the wrap angle of the strip on the roll. The key parameters and characteristics of the RSU model are discussed as follows.


## 3. Friction Analysis Based on RSU Virtual Fulcrums

### 3.1. Virtual Fulcrums of the RSU

According to the force analysis of the RSU, the bending moments are linearly distributed in the length direction of the strip with a maximum at the apex of the roll. The length direction of the strip is set to theX-axis, and xi and xi+1 are the x-coordinates of the center point of the roll i and the roll i+1, respectively, and xi=0 (i.e., the center of the roll i+1 is at the coordinate origin). Mi and Mi+1 are the maximum bending moments applied to the strip at the roll i and the roll i+1, respectively. The x-coordinate of the virtual fulcrum B on the right side of the roll i is expressed as(1)BX=xi+12⋅MiMi+Mi+1.The relationship between the bending moment and the deflectionfi of the strip is calculated by the Mohr integral method [34]:(2)fi=∫lKMidx,where Mi is the bending moment generated by the unit force acting on the desired position.According to the principle of material mechanics, in the elastic deformation range, the relationship between curvature and moment is(3)K=MEI.In the elastic-plastic deformation range, the relationship between curvature and moment is(4)K=MWEI3−2M/MW,where MW is the maximum elastic bending moment. When M is the maximum elastic bending moment, then equation (4) can be expressed the same as equation (3).Equations (3) and (4) are substituted into equation (2) and integrated to obtain(5)fi=fWΛ25−3+Λ3−2Λ,where fW is the maximum elastic bending deflection given as fW=AXBXMW/3EI, in which AX and BX are the coordinates of the virtual fulcrums A and B, respectively, and Λ=Mi/MWfW.Equation (5) determines the relationship between the total bending deflection fi and the maximum bending moment Mi of the strip at the roll i. fi=fi−1′+fi″, where fi−1′ is the residual deflection of the strip at the roll i−1 and fi″ is the reverse bending deflection of the strip at the roll i. fi″ is the sum of the residual curvature fi′ and the elastic recovery deflection fi'′' of the strip at the roll i. The deformation deflection of the strip is shown in Figure 4.Figure 4
Deformation deflection of the strip.During the leveling process, the original curvature of the strip and the relative position between rolls are known, so the total bending deflectionfi of the strip at each roll can be calculated. If the final curvature of the strip after leveling is zero, i.e., the strip is completely leveled, the positions of virtual fulcrums can be calculated by the iteration method, as shown in Figure 5.Figure 5
Solution process of the zero bending moment point of the strip.


### 3.2. Friction between the Roll and the Strip

In the RSU, the speed of the roll and the strip must be the same at the contact point; otherwise, there will be relative sliding between them, which will affect the surface quality of the strip and the service life of the roll. Therefore, the friction between the roll and the strip is static and depends on the speed relationship between them.A microelement of the strip at the point of contact with the roll is taken out for force analysis, as shown in Figure6.Figure 6
Forces acting on the microelement of the strip.
(a)
(b)When the difference between the tensions acting on the strip section corresponding to the wrap angle on the roll isΔT=T5−T3, the difference between the tensions acting on the microelement corresponding to the angle dα is(6)T′−T=ΔTrα1+α2⋅rdα=dαα1+α2ΔT.Since the strip deformation is continuous, the distribution of the leveling force on the strip is fairly smooth, so for the convenience of calculation, the leveling force applied by the roll to the strip can be considered a linear distribution, as shown in Figure7.Figure 7
Distribution of leveling force.If the strip width is a unit size, the area of the triangle in Figure7 is equal to the total leveling force P, and the maximum leveling force Pm is(7)Pm=2Pα1+α2.OncePm is determined, the distribution curve of leveling force is also determined. The force relationship of the microelement is expressed as(8)dF+T′−T+ma=dF+dαα1+α2⋅T+ρlrdα⋅a=0.Equation (8) is derived as follows:(9)dF=−ΔTα1+α2+ρlradα.Equation (9) is integrated to obtain the static friction F:(10)F=−∫0αΔTα1+α2+ρlradα=−ΔT−ρlrα1+α2a.The maximum static friction between the strip and the roll is(11)Fmax=−∫0αP′αμdα,where P′α is used to represent the positive pressure P′ acting on the microelement. The normal working condition of the RSU is F<Fmax.


## 3.1. Virtual Fulcrums of the RSU

According to the force analysis of the RSU, the bending moments are linearly distributed in the length direction of the strip with a maximum at the apex of the roll. The length direction of the strip is set to theX-axis, and xi and xi+1 are the x-coordinates of the center point of the roll i and the roll i+1, respectively, and xi=0 (i.e., the center of the roll i+1 is at the coordinate origin). Mi and Mi+1 are the maximum bending moments applied to the strip at the roll i and the roll i+1, respectively. The x-coordinate of the virtual fulcrum B on the right side of the roll i is expressed as(1)BX=xi+12⋅MiMi+Mi+1.The relationship between the bending moment and the deflectionfi of the strip is calculated by the Mohr integral method [34]:(2)fi=∫lKMidx,where Mi is the bending moment generated by the unit force acting on the desired position.According to the principle of material mechanics, in the elastic deformation range, the relationship between curvature and moment is(3)K=MEI.In the elastic-plastic deformation range, the relationship between curvature and moment is(4)K=MWEI3−2M/MW,where MW is the maximum elastic bending moment. When M is the maximum elastic bending moment, then equation (4) can be expressed the same as equation (3).Equations (3) and (4) are substituted into equation (2) and integrated to obtain(5)fi=fWΛ25−3+Λ3−2Λ,where fW is the maximum elastic bending deflection given as fW=AXBXMW/3EI, in which AX and BX are the coordinates of the virtual fulcrums A and B, respectively, and Λ=Mi/MWfW.Equation (5) determines the relationship between the total bending deflection fi and the maximum bending moment Mi of the strip at the roll i. fi=fi−1′+fi″, where fi−1′ is the residual deflection of the strip at the roll i−1 and fi″ is the reverse bending deflection of the strip at the roll i. fi″ is the sum of the residual curvature fi′ and the elastic recovery deflection fi'′' of the strip at the roll i. The deformation deflection of the strip is shown in Figure 4.Figure 4
Deformation deflection of the strip.During the leveling process, the original curvature of the strip and the relative position between rolls are known, so the total bending deflectionfi of the strip at each roll can be calculated. If the final curvature of the strip after leveling is zero, i.e., the strip is completely leveled, the positions of virtual fulcrums can be calculated by the iteration method, as shown in Figure 5.Figure 5
Solution process of the zero bending moment point of the strip.


## 3.2. Friction between the Roll and the Strip

In the RSU, the speed of the roll and the strip must be the same at the contact point; otherwise, there will be relative sliding between them, which will affect the surface quality of the strip and the service life of the roll. Therefore, the friction between the roll and the strip is static and depends on the speed relationship between them.A microelement of the strip at the point of contact with the roll is taken out for force analysis, as shown in Figure6.Figure 6
Forces acting on the microelement of the strip.
(a)
(b)When the difference between the tensions acting on the strip section corresponding to the wrap angle on the roll isΔT=T5−T3, the difference between the tensions acting on the microelement corresponding to the angle dα is(6)T′−T=ΔTrα1+α2⋅rdα=dαα1+α2ΔT.Since the strip deformation is continuous, the distribution of the leveling force on the strip is fairly smooth, so for the convenience of calculation, the leveling force applied by the roll to the strip can be considered a linear distribution, as shown in Figure7.Figure 7
Distribution of leveling force.If the strip width is a unit size, the area of the triangle in Figure7 is equal to the total leveling force P, and the maximum leveling force Pm is(7)Pm=2Pα1+α2.OncePm is determined, the distribution curve of leveling force is also determined. The force relationship of the microelement is expressed as(8)dF+T′−T+ma=dF+dαα1+α2⋅T+ρlrdα⋅a=0.Equation (8) is derived as follows:(9)dF=−ΔTα1+α2+ρlradα.Equation (9) is integrated to obtain the static friction F:(10)F=−∫0αΔTα1+α2+ρlradα=−ΔT−ρlrα1+α2a.The maximum static friction between the strip and the roll is(11)Fmax=−∫0αP′αμdα,where P′α is used to represent the positive pressure P′ acting on the microelement. The normal working condition of the RSU is F<Fmax.


## 4. Analysis of the Strip Deformation

### 4.1. Fitting of the Deformation Curve

The length of the strip section in the RSU and its wrap angle on the roll are determined by the deformation curve of the strip in leveling.Traditional curve fitting methods consider that the strip and the roll are in line contact, and the influence of the tension is not taken into account in the calculation. However, the deformation curve is greatly affected by the tension in practice, and the wrap angle is not zero. In this paper, an improved curve fitting method based on traditional methods is proposed to make the fitting curve closer to the actual deformation curve, as shown in Figure8.Figure 8
Improved fitting of the deformation curve.Suppose the strip tension is small enough, the contact between the roll and the strip is a point. In this case, the strip deformation is continuous and the deformation curve is continuous and smooth everywhere with the first- and second-order continuity at the discontinuity point. A quartic polynomialy=ax4+bx3+cx2+dx+e with to-be-determined coefficients is used to fit the deformation curve of the strip with pure bending [35], as shown by the red line in Figure 8. The first-order derivatives at points Ci and Ci+1 and the second-order derivatives at point B are all zero. The following equation is derived based on the illustrated geometric relationships:(12)y=ax4+bx3+cx2+dx+ex=0=ϕi2,y=ax4+bx3+cx2+dx+ex=xi+1=yi+1−ϕi+12,y′=4ax3+3bx2+2cx+dx=0=0,y′=4ax3+3bx2+2cx+dx=xi+1=0,y″=12ax2+6bx+2cx=BX=0.Equation (12) is solved as follows:(13)a=3xi+1−2BXϕi+ϕi+1−2yi+12xi+136BX2−6BXxi+1+xi+12,b=23Bx2−xi+12ϕi+ϕi+1−2yi+1xi+136BX2−6BXxi+1+xi+12,c=3BX2xi+1−3BXϕi+ϕi+1−2yi+1xi+126BX2−6BXxi+1+xi+12,d=0,e=ϕi2.When there is tension in the leveling and the tension is large enough, the deformation curve of the strip is changed, as shown by the blue line in Figure8, which is fitted by the arc-tangent-arc, and the tangent points are Di and Di+1, respectively. The theoretical half wrap angles of the strip on the rolls are αi and αi′, respectively. Therefore, the deformation curve of the strip is completely determined by the coordinates of the two tangent points. The coordinates of Di and Di+1 are x0,y0 and x1,y1, respectively, and the following equation is obtained according to the geometric relationship mentioned above:(14)x02+y02=ϕi24,x1−xi+12+y1−yi+12=ϕi+124,y0x0⋅y1−y0x1−x0=−1,y1−y0x1−x0⋅y1−yi+1x1−xi+1=−1.Theoretical half wrap angles are calculated according to equation (14) as follows:(15)αi=αi′=arctanx0y0.Strip deformation with pure bending and that with large tension are two extreme cases of actual deformation. Therefore, the actual deformation curve of the strip is located in the space surrounded by the two curves formed by the two extreme cases and cross intersection pointQ of these two curves.If the distance between two rolls is small and the intermesh of the roll is large, a part of the strip will also completely fit the roll to form a wrap angle even if the strip is purely bent. Therefore, on the basis of the traditional curve fitting method, this paper introduces a tension factorλ to optimize the fitting curve. The improved algorithm flow is shown in Figure 9.Figure 9
Improved curve fitting algorithm.The tension factorλ is mainly affected by the strip stiffness, the tension, and the roll diameter and can be obtained by experiment and interpolation with the value range 0,1. λ increases with the increase of tension, decreases with the increase of stiffness, and increases with the increase of roll diameter. In leveling practice, the strip stiffness and the roll diameter are constant, so the tension plays a major role.The strip is assumed to be tangent to the rolls at both the entrance and the exit of the leveler; then, the deformation curve of the strip is obtained and denoted asyx via the curve fitting algorithm proposed above. The curvature of any point on the strip deformation curve y=yx can be calculated as(16)1ρ=y″1+y23/2.By fitting the deformation curve of the strip, the deformation curvature at any point of the strip and the wrap angle on each roll are obtained, which are the necessary parameters for dynamic modeling of the RSU.


### 4.2. Deformation Function of the Strip

When the deformation stress is less than the yield stress, the strip is elastically deformed; otherwise, the strip is plastically deformed. Elastic deformation is recoverable, while plastic deformation is unrecoverable and has an effect on the movement speed of the strip. The deformation of the strip in the RSU may be only elastic or both elastic and plastic depending on the stress distribution. Therefore, as long as the length of the plastic deformation section of the strip is determined, the length of the elastic deformation section is also determined.A microelement is taken from the strip for analysis. The main stress acting on the microelement is tensionT and bending moment M. The stress distribution of the microelement obtained by superimposing these two stresses is shown in Figure 10, where εp denotes the tension elongation rate and εs denotes the bending elongation rate [34].Figure 10
Stress distribution of a strip microelement. (a) Stress distribution with bending moment. (b) Stress distribution with tension. (c) Stress distribution with bending moment and tension.
(a)
(b)
(c)In Figure10, when the tension and the bending moment act on the microelement, the deformation neutral layer deviates from the thick symmetric layer with a distance ah. According to the force balance condition in the horizontal direction, it can be concluded as equation (17.1) and can be further simplified as equation (17.2):(17.1)ah+h/2ah−h/2=εs+εpεs−εp,(17.2)ah=h2⋅εsεp.According to the fitting curve of the deformed strip, the bending stress distribution inside the strip can be obtained by the flow shown in Figure9. Since the initial curvature of the strip is substantially uniform in mass production, the total bending deflection of the strip at the roll can be considered constant during leveling; that is, the bending moment inside the strip is constant.The tension is calculated by the velocity difference between adjacent velocity nodes. The change of the tension acting on the strip between velocity nodesvi and vi+1 from time t0 to t1 is expressed as(18)ΔT=Ti+1−Ti=sEl∫t0t1vi+1−vidt,where l=∫xixi+11+y′2dx. Since the initial value of the tension acting on the strip is 0, the tension acting on each segment at any time is determined by the velocity-time function vit of each node.Based on the known bending moment and tension, the stress at any position of the strip at any time can be determined by superimposing the bending stress and the tensile stress, and the elastic and plastic deformation sections can be determined depending on whether the stress exceeds the material yield limit. The plastic deformation of the strip consumes energy and may cause an extension in the length direction. The work done by external force in leveling is equal to the work done by internal force. If a strip with a lengthds is bent from the curvature 1/R1 to 1/R2, the work dA done by the bending moment M acting on the strip is [34](19)dA=MdsR2−dsR1=Mds1R2−1R1.For a strip with the abscissa fromxi to xi+1, when the total bending curvature is 1/R and the springback curvature is 1/ρ, the work done by the bending moment M is(20)A=Ml1R−12ρ=M∫xixi+11+y′2dx⋅1R−12ρ.Because part of the energy inA is returned to the system due to the springback of the strip, the actual consumed power Ap is(21)Ap=M∫xixi+11+y′2dx⋅1R−1ρ.In the range of the abscissaxi,xi+1, the work Aξ required for the length extension of the strip is(22)Aξ=σsbhlξ=σsbh∫xixi+11+y′2ξxdx.Because the tension and the bending moment act together on the strip, the deformation neutral layer deviates from the thick symmetric layer with a distancea, and the thick symmetric layer is stretched. If the stretch is plastic, the strip is extended. The tensile deformation of the thick symmetric layer plays an important role in improving the strip shape, but it is difficult to quantitatively calculate the deformation. At present, it is considered that the strip is plastically extended only when the thickness of the one-sided elastic deformation z0 is less than a, and the plastic deformation εr of the thick symmetric layer of the strip is as follows [34]:(23)εr=σsEσpσsKz−11−ηKz−11+σp/σsKz−1.The elongation rateξx at a point x,yx on the fitting curve is(24)ξx=σsEaxz0x−11−ηh/2z0x−11+2ax/hh/2z0x−1.


## 4.1. Fitting of the Deformation Curve

The length of the strip section in the RSU and its wrap angle on the roll are determined by the deformation curve of the strip in leveling.Traditional curve fitting methods consider that the strip and the roll are in line contact, and the influence of the tension is not taken into account in the calculation. However, the deformation curve is greatly affected by the tension in practice, and the wrap angle is not zero. In this paper, an improved curve fitting method based on traditional methods is proposed to make the fitting curve closer to the actual deformation curve, as shown in Figure8.Figure 8
Improved fitting of the deformation curve.Suppose the strip tension is small enough, the contact between the roll and the strip is a point. In this case, the strip deformation is continuous and the deformation curve is continuous and smooth everywhere with the first- and second-order continuity at the discontinuity point. A quartic polynomialy=ax4+bx3+cx2+dx+e with to-be-determined coefficients is used to fit the deformation curve of the strip with pure bending [35], as shown by the red line in Figure 8. The first-order derivatives at points Ci and Ci+1 and the second-order derivatives at point B are all zero. The following equation is derived based on the illustrated geometric relationships:(12)y=ax4+bx3+cx2+dx+ex=0=ϕi2,y=ax4+bx3+cx2+dx+ex=xi+1=yi+1−ϕi+12,y′=4ax3+3bx2+2cx+dx=0=0,y′=4ax3+3bx2+2cx+dx=xi+1=0,y″=12ax2+6bx+2cx=BX=0.Equation (12) is solved as follows:(13)a=3xi+1−2BXϕi+ϕi+1−2yi+12xi+136BX2−6BXxi+1+xi+12,b=23Bx2−xi+12ϕi+ϕi+1−2yi+1xi+136BX2−6BXxi+1+xi+12,c=3BX2xi+1−3BXϕi+ϕi+1−2yi+1xi+126BX2−6BXxi+1+xi+12,d=0,e=ϕi2.When there is tension in the leveling and the tension is large enough, the deformation curve of the strip is changed, as shown by the blue line in Figure8, which is fitted by the arc-tangent-arc, and the tangent points are Di and Di+1, respectively. The theoretical half wrap angles of the strip on the rolls are αi and αi′, respectively. Therefore, the deformation curve of the strip is completely determined by the coordinates of the two tangent points. The coordinates of Di and Di+1 are x0,y0 and x1,y1, respectively, and the following equation is obtained according to the geometric relationship mentioned above:(14)x02+y02=ϕi24,x1−xi+12+y1−yi+12=ϕi+124,y0x0⋅y1−y0x1−x0=−1,y1−y0x1−x0⋅y1−yi+1x1−xi+1=−1.Theoretical half wrap angles are calculated according to equation (14) as follows:(15)αi=αi′=arctanx0y0.Strip deformation with pure bending and that with large tension are two extreme cases of actual deformation. Therefore, the actual deformation curve of the strip is located in the space surrounded by the two curves formed by the two extreme cases and cross intersection pointQ of these two curves.If the distance between two rolls is small and the intermesh of the roll is large, a part of the strip will also completely fit the roll to form a wrap angle even if the strip is purely bent. Therefore, on the basis of the traditional curve fitting method, this paper introduces a tension factorλ to optimize the fitting curve. The improved algorithm flow is shown in Figure 9.Figure 9
Improved curve fitting algorithm.The tension factorλ is mainly affected by the strip stiffness, the tension, and the roll diameter and can be obtained by experiment and interpolation with the value range 0,1. λ increases with the increase of tension, decreases with the increase of stiffness, and increases with the increase of roll diameter. In leveling practice, the strip stiffness and the roll diameter are constant, so the tension plays a major role.The strip is assumed to be tangent to the rolls at both the entrance and the exit of the leveler; then, the deformation curve of the strip is obtained and denoted asyx via the curve fitting algorithm proposed above. The curvature of any point on the strip deformation curve y=yx can be calculated as(16)1ρ=y″1+y23/2.By fitting the deformation curve of the strip, the deformation curvature at any point of the strip and the wrap angle on each roll are obtained, which are the necessary parameters for dynamic modeling of the RSU.


## 4.2. Deformation Function of the Strip

When the deformation stress is less than the yield stress, the strip is elastically deformed; otherwise, the strip is plastically deformed. Elastic deformation is recoverable, while plastic deformation is unrecoverable and has an effect on the movement speed of the strip. The deformation of the strip in the RSU may be only elastic or both elastic and plastic depending on the stress distribution. Therefore, as long as the length of the plastic deformation section of the strip is determined, the length of the elastic deformation section is also determined.A microelement is taken from the strip for analysis. The main stress acting on the microelement is tensionT and bending moment M. The stress distribution of the microelement obtained by superimposing these two stresses is shown in Figure 10, where εp denotes the tension elongation rate and εs denotes the bending elongation rate [34].Figure 10
Stress distribution of a strip microelement. (a) Stress distribution with bending moment. (b) Stress distribution with tension. (c) Stress distribution with bending moment and tension.
(a)
(b)
(c)In Figure10, when the tension and the bending moment act on the microelement, the deformation neutral layer deviates from the thick symmetric layer with a distance ah. According to the force balance condition in the horizontal direction, it can be concluded as equation (17.1) and can be further simplified as equation (17.2):(17.1)ah+h/2ah−h/2=εs+εpεs−εp,(17.2)ah=h2⋅εsεp.According to the fitting curve of the deformed strip, the bending stress distribution inside the strip can be obtained by the flow shown in Figure9. Since the initial curvature of the strip is substantially uniform in mass production, the total bending deflection of the strip at the roll can be considered constant during leveling; that is, the bending moment inside the strip is constant.The tension is calculated by the velocity difference between adjacent velocity nodes. The change of the tension acting on the strip between velocity nodesvi and vi+1 from time t0 to t1 is expressed as(18)ΔT=Ti+1−Ti=sEl∫t0t1vi+1−vidt,where l=∫xixi+11+y′2dx. Since the initial value of the tension acting on the strip is 0, the tension acting on each segment at any time is determined by the velocity-time function vit of each node.Based on the known bending moment and tension, the stress at any position of the strip at any time can be determined by superimposing the bending stress and the tensile stress, and the elastic and plastic deformation sections can be determined depending on whether the stress exceeds the material yield limit. The plastic deformation of the strip consumes energy and may cause an extension in the length direction. The work done by external force in leveling is equal to the work done by internal force. If a strip with a lengthds is bent from the curvature 1/R1 to 1/R2, the work dA done by the bending moment M acting on the strip is [34](19)dA=MdsR2−dsR1=Mds1R2−1R1.For a strip with the abscissa fromxi to xi+1, when the total bending curvature is 1/R and the springback curvature is 1/ρ, the work done by the bending moment M is(20)A=Ml1R−12ρ=M∫xixi+11+y′2dx⋅1R−12ρ.Because part of the energy inA is returned to the system due to the springback of the strip, the actual consumed power Ap is(21)Ap=M∫xixi+11+y′2dx⋅1R−1ρ.In the range of the abscissaxi,xi+1, the work Aξ required for the length extension of the strip is(22)Aξ=σsbhlξ=σsbh∫xixi+11+y′2ξxdx.Because the tension and the bending moment act together on the strip, the deformation neutral layer deviates from the thick symmetric layer with a distancea, and the thick symmetric layer is stretched. If the stretch is plastic, the strip is extended. The tensile deformation of the thick symmetric layer plays an important role in improving the strip shape, but it is difficult to quantitatively calculate the deformation. At present, it is considered that the strip is plastically extended only when the thickness of the one-sided elastic deformation z0 is less than a, and the plastic deformation εr of the thick symmetric layer of the strip is as follows [34]:(23)εr=σsEσpσsKz−11−ηKz−11+σp/σsKz−1.The elongation rateξx at a point x,yx on the fitting curve is(24)ξx=σsEaxz0x−11−ηh/2z0x−11+2ax/hh/2z0x−1.


## 5. AMESim Modeling of the RSU and Its Application in Rolling Production Line

### 5.1. AMESim Modeling of the RSU

According to the RSU model defined above, AMESim is used to implement the dynamic modeling, as shown in Figure11. The “roll/strip contact calculation unit” is a custom-developed component, which will be described below. The strip is divided into four mass-spring-damping segments, corresponding to four deformation sections. The motor is driven by a control signal, acts on a roll system with an equivalent rotational inertia of J via a spring-damping system, and finally acts on the strip by a roll. Adjacent RSUs are coupled by speed.Figure 11
Dynamic model of the RSU implemented in AMESim.The force relationship between the roll and the strip includes two aspects: one is the bending moment and the force at the contact point of the strip and roll in the tangential direction applied to the strip by the roll and the other is the additional torque applied to the roll by the strip. According to the analysis of the force relationship between the roll and the strip, the implementation of the “roll/strip contact calculation unit” is shown in Figure12.Figure 12
Roll/strip contact calculation unit.AMECustom is used to customize and package the RSU to establish a new basic element and its submodel, as shown in Figure13, which is also one of the basic components for building the AMESim model of a leveler.Figure 13
Customization of the RSU.


### 5.2. RSU Application in Rolling Production Line AMESim Modeling

The rolling production line is a highly integrated set of equipment which is designed and customized according to the process requirements of customers and the layout of the specific production environment [36]. Different rolling production lines contain different equipment and different layout of equipment; nevertheless, the rolling production line must contain the following basic equipment: uncoiler, pinch roller, rolling machine, tension roller, and winder. The basic equipment in the rolling production line is shown in Figure 14.Figure 14
Typical equipment in the rolling production line.


#### 5.2.1. Uncoiler and Winder

Due to the working stability, the uncoiler and the winder can be, respectively, defined as a single-input and a single-output roll/strip basic unit, as shown in Figure15. This process is also considered to simplify the modeling and improve the simulation speed.Figure 15
Uncoiler: (a) AMESim model and (b) customized components; winder: (c) AMESim model and (d) customized components.
(a)
(b)
(c)
(d)


#### 5.2.2. Pinch Roller

The pinch roller can be simplified as a basic roll/strip unit. According to the working principle of the pinch roller, the positive pressure is applied by the lower roller to the upper roller, which is irrelevant to the bending deformation of the strip, and is constant in the working process. When the strip passes through the pinch roller, the wrap angle is very small, approximating to zero. Considering the working state of the whole production line, there is no plastic deformation area in the equivalent roll/belt basic unit where the pinch roll is located. The AMESim model and its customized components of the pinch roller are shown in Figure16.Figure 16
Pinch roller: (a) AMESim model and (b) customized components.
(a)
(b)


#### 5.2.3. Tension Roller

The principle of the tension roller is that when the strip is wrapped around the tension roll, the friction force is produced in the contact arc, which enlarges or reduces the tension of the strip at the entrance and exit of the tension roll. In the production process, the practical wrap angle is always smaller than its theoretical value, and with the increase of plate stiffness, the practical wrap angle is smaller. The AMESim model and its customized components of the tension roller are shown in Figure17.Figure 17
Tension roller: (a) AMESim model and (b) customized components.
(a)
(b)


#### 5.2.4. Multiroller Leveling Machine

At present, most of the multiroller leveling machines are second-generation or second-and-a-half-generation products. In the mechanical configuration, the diameter of the working roll is uniform. The centralized driving mode is adopted in the traditional systems. The working rollers are driven by a high-power motor through coupling, gearbox, and universal coupling. The reduction of the roll system is adjusted by the inclination and reduction of the upper roll system. We take a 23-roller leveling machine designed and manufactured by Boya Precision Industrial Equipments Co., Ltd., as an example to illustrate the multiroller leveling machine modeling process using the RSU. The AMESim model and its customized components of the multiroller leveling machine are shown in Figure18.Figure 18
Multiroller leveling machine: (a) AMESim model and (b) customized components.
(a)
(b)To synthesize the above modeling, we can obtain the whole rolling production line AMESim modeling, as shown in Figure19.Figure 19
Whole rolling production line AMESim modeling.


### 5.3. Rolling Production Line Association Modeling System

The rolling production line AMESim modeling based on the RSU is utilized in developing the rolling production line association modeling system. The modeling system integrates MATLAB into the dynamic model of the rolling production line in AMESim. It helps to assign parameters, to read variables, and in simulation control. The specific framework of the rolling production line fast association modeling system is shown in Figure20.Figure 20
Framework of the rolling production line fast association modeling system.The initial parameter setting module calculates the parameters of each basic functional unit in the model according to the input information, including the inherent parameters of the equipment (roll diameter, roll distance, etc.) and the processing parameters (downpressure, strip material, and design tension). It is also responsible for transmitting the specific values of the corresponding parameters to the AMESim model after the calculation is completed. The parameters and variables of each submodel in the AMESim model can be found through a unique handle and read and write. When this module is initially executed, it will require basic information of the straightening production line model, including the model file name and device corresponding to the model instance name. Then, it can automatically update model parameters when it is executed again. The wizard parameter setting method is provided in the module to ensure the accuracy of the parameter input. The basic information of the rolling production line can be obtained in the function module shown in Figure21. The process parameters of the rolling production line can be obtained in the function module shown in Figure 22. After obtaining the basic parameters, the simulation process can be implemented. In the process of simulation, the parameters of the model are not a constant but a function of other parameters because of the elastic-plastic deformation of the strip. It may change in the process of simulation; for example, the length of the elastic interval of the strip may change when the tension changes. The function of the simulation execution control module is to recalculate the parameters and set the model after each step of simulation within a given integration time interval so as to achieve the purpose of making the simulation results more accurate. It should be pointed out that the efficiency of reading AMESim results in the MATLAB environment is low, so the module also provides the reading frequency control options for the simulation results. Figure 23 shows the setting function of the simulation execution module.Figure 21
Basic information of the rolling production line.Figure 22
Process parameters of the rolling production line.Figure 23
Simulation execution module.AMESim returns the simulation results to MATLAB in the form of two matrices at the end of each simulation stage, one of which is the result matrix, which saves the final values of all variables in the model at the end of simulation, and the other matrix saves the name information of all variables. When the information is transmitted to the platform, it needs to be decoded by the recognition module read by the simulation results to identify the required information and provide it to other modules for use, as shown in Figure24. Take tensor calculation at each node as an example. According to equation (12), the tensor difference between node 1 and node 2 can be calculated as(25)ΔT=T2−T1=sEl∫t1t2v2−v1dt=20.1N,where s and E are obtained in the module shown in Figure 21. The strip length between nodes can be calculated as l=∫x1x21+y′2dx, where y′ is expressed based on equations (6) and (7). So on the premise that T1 is known, we can obtain T2. Analogically, the tensions at other nodes are calculated and shown in the module in Figure 24. With the help of the rolling production line association modeling system based on the RSU, the parameter information of the complex multiroll leveling manufacturing system can be obtained more conveniently and the related performance parameters can be calculated quickly, which is helpful to realize the digitalization of the complex multiroll leveling manufacturing.Figure 24
Simulation result reading module.AMESim is a graphic modeling and simulation software program based on power flow data. For the convenience of calling, AMECustom is used to customize the RSU AMESim modeling process. It not only facilitated the centralized setting of submodel parameters but also reduced the complexity of subsequent processes. By analyzing the working principle of each equipment in the rolling production line, the dynamic model can be built in AMESim. Finally, the dynamic model database of the rolling equipment is established, and it helps to create the dynamic model of the rolling production line more efficiently.


## 5.1. AMESim Modeling of the RSU

According to the RSU model defined above, AMESim is used to implement the dynamic modeling, as shown in Figure11. The “roll/strip contact calculation unit” is a custom-developed component, which will be described below. The strip is divided into four mass-spring-damping segments, corresponding to four deformation sections. The motor is driven by a control signal, acts on a roll system with an equivalent rotational inertia of J via a spring-damping system, and finally acts on the strip by a roll. Adjacent RSUs are coupled by speed.Figure 11
Dynamic model of the RSU implemented in AMESim.The force relationship between the roll and the strip includes two aspects: one is the bending moment and the force at the contact point of the strip and roll in the tangential direction applied to the strip by the roll and the other is the additional torque applied to the roll by the strip. According to the analysis of the force relationship between the roll and the strip, the implementation of the “roll/strip contact calculation unit” is shown in Figure12.Figure 12
Roll/strip contact calculation unit.AMECustom is used to customize and package the RSU to establish a new basic element and its submodel, as shown in Figure13, which is also one of the basic components for building the AMESim model of a leveler.Figure 13
Customization of the RSU.


## 5.2. RSU Application in Rolling Production Line AMESim Modeling

The rolling production line is a highly integrated set of equipment which is designed and customized according to the process requirements of customers and the layout of the specific production environment [36]. Different rolling production lines contain different equipment and different layout of equipment; nevertheless, the rolling production line must contain the following basic equipment: uncoiler, pinch roller, rolling machine, tension roller, and winder. The basic equipment in the rolling production line is shown in Figure 14.Figure 14
Typical equipment in the rolling production line.


### 5.2.1. Uncoiler and Winder

Due to the working stability, the uncoiler and the winder can be, respectively, defined as a single-input and a single-output roll/strip basic unit, as shown in Figure15. This process is also considered to simplify the modeling and improve the simulation speed.Figure 15
Uncoiler: (a) AMESim model and (b) customized components; winder: (c) AMESim model and (d) customized components.
(a)
(b)
(c)
(d)


### 5.2.2. Pinch Roller

The pinch roller can be simplified as a basic roll/strip unit. According to the working principle of the pinch roller, the positive pressure is applied by the lower roller to the upper roller, which is irrelevant to the bending deformation of the strip, and is constant in the working process. When the strip passes through the pinch roller, the wrap angle is very small, approximating to zero. Considering the working state of the whole production line, there is no plastic deformation area in the equivalent roll/belt basic unit where the pinch roll is located. The AMESim model and its customized components of the pinch roller are shown in Figure16.Figure 16
Pinch roller: (a) AMESim model and (b) customized components.
(a)
(b)


### 5.2.3. Tension Roller

The principle of the tension roller is that when the strip is wrapped around the tension roll, the friction force is produced in the contact arc, which enlarges or reduces the tension of the strip at the entrance and exit of the tension roll. In the production process, the practical wrap angle is always smaller than its theoretical value, and with the increase of plate stiffness, the practical wrap angle is smaller. The AMESim model and its customized components of the tension roller are shown in Figure17.Figure 17
Tension roller: (a) AMESim model and (b) customized components.
(a)
(b)


### 5.2.4. Multiroller Leveling Machine

At present, most of the multiroller leveling machines are second-generation or second-and-a-half-generation products. In the mechanical configuration, the diameter of the working roll is uniform. The centralized driving mode is adopted in the traditional systems. The working rollers are driven by a high-power motor through coupling, gearbox, and universal coupling. The reduction of the roll system is adjusted by the inclination and reduction of the upper roll system. We take a 23-roller leveling machine designed and manufactured by Boya Precision Industrial Equipments Co., Ltd., as an example to illustrate the multiroller leveling machine modeling process using the RSU. The AMESim model and its customized components of the multiroller leveling machine are shown in Figure18.Figure 18
Multiroller leveling machine: (a) AMESim model and (b) customized components.
(a)
(b)To synthesize the above modeling, we can obtain the whole rolling production line AMESim modeling, as shown in Figure19.Figure 19
Whole rolling production line AMESim modeling.


## 5.2.1. Uncoiler and Winder

Due to the working stability, the uncoiler and the winder can be, respectively, defined as a single-input and a single-output roll/strip basic unit, as shown in Figure15. This process is also considered to simplify the modeling and improve the simulation speed.Figure 15
Uncoiler: (a) AMESim model and (b) customized components; winder: (c) AMESim model and (d) customized components.
(a)
(b)
(c)
(d)


## 5.2.2. Pinch Roller

The pinch roller can be simplified as a basic roll/strip unit. According to the working principle of the pinch roller, the positive pressure is applied by the lower roller to the upper roller, which is irrelevant to the bending deformation of the strip, and is constant in the working process. When the strip passes through the pinch roller, the wrap angle is very small, approximating to zero. Considering the working state of the whole production line, there is no plastic deformation area in the equivalent roll/belt basic unit where the pinch roll is located. The AMESim model and its customized components of the pinch roller are shown in Figure16.Figure 16
Pinch roller: (a) AMESim model and (b) customized components.
(a)
(b)


## 5.2.3. Tension Roller

The principle of the tension roller is that when the strip is wrapped around the tension roll, the friction force is produced in the contact arc, which enlarges or reduces the tension of the strip at the entrance and exit of the tension roll. In the production process, the practical wrap angle is always smaller than its theoretical value, and with the increase of plate stiffness, the practical wrap angle is smaller. The AMESim model and its customized components of the tension roller are shown in Figure17.Figure 17
Tension roller: (a) AMESim model and (b) customized components.
(a)
(b)


## 5.2.4. Multiroller Leveling Machine

At present, most of the multiroller leveling machines are second-generation or second-and-a-half-generation products. In the mechanical configuration, the diameter of the working roll is uniform. The centralized driving mode is adopted in the traditional systems. The working rollers are driven by a high-power motor through coupling, gearbox, and universal coupling. The reduction of the roll system is adjusted by the inclination and reduction of the upper roll system. We take a 23-roller leveling machine designed and manufactured by Boya Precision Industrial Equipments Co., Ltd., as an example to illustrate the multiroller leveling machine modeling process using the RSU. The AMESim model and its customized components of the multiroller leveling machine are shown in Figure18.Figure 18
Multiroller leveling machine: (a) AMESim model and (b) customized components.
(a)
(b)To synthesize the above modeling, we can obtain the whole rolling production line AMESim modeling, as shown in Figure19.Figure 19
Whole rolling production line AMESim modeling.


## 5.3. Rolling Production Line Association Modeling System

The rolling production line AMESim modeling based on the RSU is utilized in developing the rolling production line association modeling system. The modeling system integrates MATLAB into the dynamic model of the rolling production line in AMESim. It helps to assign parameters, to read variables, and in simulation control. The specific framework of the rolling production line fast association modeling system is shown in Figure20.Figure 20
Framework of the rolling production line fast association modeling system.The initial parameter setting module calculates the parameters of each basic functional unit in the model according to the input information, including the inherent parameters of the equipment (roll diameter, roll distance, etc.) and the processing parameters (downpressure, strip material, and design tension). It is also responsible for transmitting the specific values of the corresponding parameters to the AMESim model after the calculation is completed. The parameters and variables of each submodel in the AMESim model can be found through a unique handle and read and write. When this module is initially executed, it will require basic information of the straightening production line model, including the model file name and device corresponding to the model instance name. Then, it can automatically update model parameters when it is executed again. The wizard parameter setting method is provided in the module to ensure the accuracy of the parameter input. The basic information of the rolling production line can be obtained in the function module shown in Figure21. The process parameters of the rolling production line can be obtained in the function module shown in Figure 22. After obtaining the basic parameters, the simulation process can be implemented. In the process of simulation, the parameters of the model are not a constant but a function of other parameters because of the elastic-plastic deformation of the strip. It may change in the process of simulation; for example, the length of the elastic interval of the strip may change when the tension changes. The function of the simulation execution control module is to recalculate the parameters and set the model after each step of simulation within a given integration time interval so as to achieve the purpose of making the simulation results more accurate. It should be pointed out that the efficiency of reading AMESim results in the MATLAB environment is low, so the module also provides the reading frequency control options for the simulation results. Figure 23 shows the setting function of the simulation execution module.Figure 21
Basic information of the rolling production line.Figure 22
Process parameters of the rolling production line.Figure 23
Simulation execution module.AMESim returns the simulation results to MATLAB in the form of two matrices at the end of each simulation stage, one of which is the result matrix, which saves the final values of all variables in the model at the end of simulation, and the other matrix saves the name information of all variables. When the information is transmitted to the platform, it needs to be decoded by the recognition module read by the simulation results to identify the required information and provide it to other modules for use, as shown in Figure24. Take tensor calculation at each node as an example. According to equation (12), the tensor difference between node 1 and node 2 can be calculated as(25)ΔT=T2−T1=sEl∫t1t2v2−v1dt=20.1N,where s and E are obtained in the module shown in Figure 21. The strip length between nodes can be calculated as l=∫x1x21+y′2dx, where y′ is expressed based on equations (6) and (7). So on the premise that T1 is known, we can obtain T2. Analogically, the tensions at other nodes are calculated and shown in the module in Figure 24. With the help of the rolling production line association modeling system based on the RSU, the parameter information of the complex multiroll leveling manufacturing system can be obtained more conveniently and the related performance parameters can be calculated quickly, which is helpful to realize the digitalization of the complex multiroll leveling manufacturing.Figure 24
Simulation result reading module.AMESim is a graphic modeling and simulation software program based on power flow data. For the convenience of calling, AMECustom is used to customize the RSU AMESim modeling process. It not only facilitated the centralized setting of submodel parameters but also reduced the complexity of subsequent processes. By analyzing the working principle of each equipment in the rolling production line, the dynamic model can be built in AMESim. Finally, the dynamic model database of the rolling equipment is established, and it helps to create the dynamic model of the rolling production line more efficiently.


## 6. Conclusions

In this paper, a roll-strip unit is proposed for the unified modeling of multiroll leveling for metal strips. The virtual fulcrums of the unit are defined based on the force analysis of the RSU. To analyze the friction between the roll and the strip, the deformation curve of the strip is fitted and the deformation function is constructed. The strip deformation with pure bending and that with large tension are two extreme cases of actual deformation, and the actual deformation curve of the strip is located in the space surrounded by the two curves formed in the two extreme cases. The elastic deformation of the strip is recoverable, while plastic deformation is unrecoverable and has an effect on the movement speed of the strip.The plastic deformation of the strip consumes energy and may cause an extension in the length direction, which plays an important role in improving the strip shape. The friction between the roll and the strip is static, and it depends on the speed relationship between them. The RSU is modeled in AMESim. It is applied in a 23-roller rolling production modeling. With the help of the rolling production line association modeling system based on the RSU, the parameter information of the complex multiroll leveling manufacturing system can be obtained more conveniently and the related performance parameters can be calculated quickly, which is helpful to realize the digitalization of the complex multiroll leveling manufacturing [36].



---

*Source: 1486864-2020-01-31.xml*