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

: This study establishes a vertical–horizontal coupling vibration model of hot rolling mill rolls under multi-piecewise nonlinear constraints considering the piecewise nonlinear spring force and piecewise nonlinear friction force constraints of the hydraulic cylinder in the vertical direction of the rolls, the piecewise stiffness constraints in the horizontal direction, and the influence of the nonlinear dynamic rolling force in the rolling process. Using the average method to solve the amplitude–frequency response equation of the coupled vibration system and taking the actual parameters of a 1780 mm hot rolling mill (Chengde Steel Co., Ltd., Chengde, China) as an example, we study the amplitude–frequency characteristics of the mill rolls under different parameter settings. The results show that the amplitude and resonance region can be reduced by appropriately reducing the external disturbance force and the nonlinear spring force of the hydraulic cylinder, appropriately increasing the nonlinear friction force, and eliminating the gap between the bearing seat and the mill housing, to avoid the amplitude jump phenomenon due to piecewise variation. Furthermore, using the singularity theory to study the static bifurcation characteristics of the coupled vibration system, we establish a relationship between the vibration parameters and the topological bifurcation solution of the coupled system. The transition sets and their corresponding bifurcation topological structure in three cases are given, and the steady and unsteady process parameter regions of the rolls are obtained. The dynamic behavior of the coupled vibration system can be controlled by varying the bifurcation parameter. This study provides a theoretical basis for restraining the vibration of hot rolling mill rolls and optimizing the process parameters.


Introduction
A rolling mill is an important piece of equipment in the field of modern heavy machinery, and its safe and stable operation is vital for ensuring a high-efficiency production of rolling products.However, rolling mill rolls are often unstable during the rolling process because of the high rolling speed and rolling strength of modern mills.This not only affects the quality and precision of the products but may also damage the rolling equipment [1][2][3][4].
The vibration problem in hot rolling mills has mainly been studied on the basis of the linear vibration theory in the early stage; however, practical results have shown that a linear analysis can sometimes lead to significant calculation errors.With the development of the nonlinear theory, the nonlinear vibration of hot rolling mill rolls has attracted the attention of researchers, who have studied it from various angles and achieved good results [5,6].Swiatoniowski et al. noticed that when the vertical system of a rolling mill is disturbed by an external disturbance, the stiffness of the workpiece varies periodically, leading to a nonlinear parametric vibration of the rolling mill and causing mill chatter [7].
To study the vibration characteristics of a rolling mill and to determine the vibration source, Fan et al. established a horizontal-direction friction chatter model of a roller [8].By analyzing the effect of varying the workpiece thickness and motor speed on the connection angle and roll gap friction force, Shi et al. established a nonlinear torsional vibration model of a rolling mill and showed that reducing the damping coefficient and nonlinear stiffness can help reduce the vibration intensity [9].The existing modeling and analysis of rolling mills have been conducted by considering a single-vibration direction, such as vertical vibration, horizontal vibration, or torsional vibration; this makes it impossible to accurately capture the characteristics of the rolling process and fully explain the complex phenomena that occur during vibration.In an actual rolling process, the rolling mill vibration has characteristics such as variability, time variation, nonlinearity, and multiconstraint, with most of these factors exhibiting different types of coupling vibrations.Therefore, the key to solving this problem is to analyze the rolling mill vibration from the perspective of coupling vibration [10].Liu et al. used the Hopf bifurcation theorem to analyze the stability and bifurcation characteristics of rolling mill rolls by subjecting the main drive system of a rolling mill to an electromechanical coupling vibration [11,12].Through testing and analyzing the vibration signal of a rolling mill in the rolling process, Peng et al. reported many different dynamic motion forms between the rolling mill rolls, strip steel, and rolling parameters; this interaction produces the coupling effect [13].Wang et al. established a multi-coupling vibration model comprising a rolling interface friction model, a rolling torque model, and a hydrodynamic lubrication model, and studied the interface dynamic characteristics of a high-speed rolling mill work interface based on an unsteady lubrication process [14].Xu et al. established a coupling dynamic model between a rolling mill gearbox and its various time-varying parameters and explained the dynamic vibration characteristics of the rolling mill under different rolling parameters [15].In addition, some scholars considered the interaction between the vibration of the rolling mill structure and the workpiece vibration and proposed a coupled vibration model for a rollrolled piece [16,17].Qian et al. studied the adaptive fuzzy control problem of a rolling mill vibration system and designed a controller to ensure the stability of the mechanical, hydraulic coupling rolling system [18].In summary, although some advances have been made in the behavior analysis of the coupled vibration of rolling mills, the complex vibration law of the hot rolling mill itself needs to be further explored and verified from different angles.This can help put forward more effective vibration suppression measures and thereby improve the production efficiency and product quality.In recent years, the dynamics of mechanical systems with piecewise mechanical factors, such as constraints, friction, and gap, have attracted extensive attention in the engineering field [19,20].Currently, there are few reports on the vibration of rolling mill rolls under the influence of piecewise nonlinear factors.Therefore, it is necessary to thoroughly study the influence of piecewise nonlinear factors on the vibration of rolling mill rolls.The hydraulic system used in a hot rolling mill is a typical nonlinear system.In the working process of the hydraulic cylinder, a nonlinear hydraulic stiffness is generated because of the compressibility of the oil.When a rolling mill roll vibrates, it exhibits a nonlinear spring force [21], and there is a nonlinear friction between the piston and the cylinder wall [22].Although the spring force and friction force in the hydraulic system are typically expressed in the form of weak nonlinearity, the nonlinear constraint effect cannot be ignored for high-precision position-controlled rolling mill systems.In the horizontal direction of the rolling mill rolls, there is a gap between the mill housing and the bearing seat, and because of the wear and eccentricity of the rolls, when the rolling mill rolls vibrate, the bearing seat will impact the mill housing, which will change the contact stiffness between the horizontal rolls, thus affecting the smooth operation of the rolling mill.
During rolling, horizontal vibration and vertical vibration often occur in a hot rolling mill.These interact and couple with each other [23].This study establishes a vertical-horizontal coupling vibration dynamic model of a hot rolling mill under multi-piecewise nonlinear constraints considering a dynamic rolling force model of the rolling process, a piecewise nonlinear spring force model of the hydraulic system, a piecewise nonlinear friction force model, and a piecewise stiffness model in the horizontal direction.The response of the coupled vibration system is solved using the nonlinear method.On this basis, taking the rolling parameters of an actual mill as an example, we study the effects of rolling and process parameters on the amplitude-frequency characteristics and static bifurcation characteristics of the coupled vibration system of the hot rolling mill.This research provides certain guidance for reducing and restraining the vibration of rolling mills.

Multi-Piecewise Nonlinear Constraint Modeling of Hot Rolling Mill Rolls
Figure 2 shows the structural diagram of a four-high hot rolling mill, which is mainly composed of a mill stand, backup rolls, work rolls, a hydraulic cylinder, a balance cylinder, and a bending cylinder.The hydraulic cylinder is distributed between the frame and the upper backup roll and acts on the bearing seat of the upper backup roll.The upper backup roll is balanced under the joint restraint of the hydraulic and balance cylinders.
The bending cylinder acts on the bearing seat of the upper and lower work rolls.The bending cylinder maintains contact between the work roll and the backup roll, to realize the balance between the work roll and the backup roll.As the main parameters, the work roll body length of the rolling mill is 2080 mm, and the diameter is 840 mm.The backup roll body length is 1780 mm, and the diameter is 1600 mm.The hydraulic cylinder used is a double-acting single-piston servo hydraulic cylinder, the diameter and stroke of which are 1050 mm and 100 mm, respectively.

Piecewise Nonlinear Spring Force Model
The hydraulic elastic stiffness refers to the nonlinear elastic stiffness of the liquid when the working chamber of the hydraulic cylinder is completely closed, and the volume of the oil is compressed under the action of external forces [21,23].The spring force of a hydraulic cylinder can be expressed as: where α0 and β0 are, respectively, the linear and nonlinear equivalent stiffnesses of the hydraulic cylinder in motion.The hydraulic cylinder, shown in Figure 2, is a double-acting single-piston servo hydraulic cylinder.With regard to its operation, a change in the piston displacement changes the pressure and oil volume in the two cavities of the hydraulic cylinder, thereby changing the oil stiffness.Therefore, the elastic force of the hydraulic cylinder is a nonlinear function of the vertical vibration displacement.Because of the size of the hydraulic and bending cylinders, the oil bulk elastic modulus, piston stroke, and other factors are different, and so are the equivalent dynamic stiffnesses.To simplify the research, only the upper rolls of the mill are analyzed and simplified into a lumped mass block [7].
In Figure 3, m is the equivalent mass of the upper backup roll, upper work roll, and its bearing seat, and y is the vertical vibration displacement of the rolling mill rolls.e1 and e2 are the initial elastic deformations of the hydraulic screwdown cylinder and bending cylinder, respectively, under the action of the preload during steady-state rolling (e1 > 0, e2 < 0).ks1 and kw1 are the equivalent linear stiffnesses of the hydraulic screwdown cylinder and bending cylinder, respectively, ks2 and kw2 are the equivalent nonlinear stiffnesses of the hydraulic screwdown cylinder and bending cylinder, respectively.When the rolling mill is in operation, the downward movement of the rolls exceeds the initial elastic deformation of the bending cylinder, and the elastic stiffness between the hydraulic cylinder and the upper roll system is zero; in this case, the mill is mainly affected by the elastic force of the bending cylinder.Similarly, when the upward movement of the rolls exceeds the initial elastic deformation of the hydraulic cylinder, the elastic stiffness between the balance cylinder and the upper rolls is zero; in this case, the mill is mainly affected by the elastic force of the hydraulic cylinder.When the vertical vibration displacement of the roll system is between the two, the mill is affected by the elastic forces of the hydraulic and bending cylinders.The constraint is piecewise nonlinear, so the piecewise nonlinear spring force in the vertical direction can be expressed as: Figure 4 shows the variation law of the piecewise nonlinear spring force with the vertical vibration displacement.Evidently, the stiffness changes at e1 and e2.

Piecewise Nonlinear Friction Force Model
While the hydraulic cylinder converts hydraulic energy into mechanical energy, the friction force plays an important role.The main source of the friction force is the friction between the piston rod and the cylinder wall.In the rolling process, when the upper rolls vibrate, the friction force of the hydraulic and bending cylinders is not constant but varies with the periodic change in the vertical vibration speed ẏ, so the friction force is nonlinear [24].The law of the friction coefficient between the piston and the cylinder wall of the hydraulic cylinder with the vibration velocity can be expressed as follows: where μf is the friction coefficient between the piston and the cylinder wall, μs is the static friction coefficient, μm is the maximum dynamic friction coefficient, and vm is the vertical vibration speed when the friction coefficient is μm, so the piecewise nonlinear friction force in the vertical direction can be expressed as: where pf is the pressure between the piston rod and the cylinder wall, and its value is related to the material of the hydraulic cylinder, the tightness of the seal assembly, the hardness of the seal material, and the radial component of the load.
Figure 5 shows the variation law of the friction force with the vibration velocity.The friction force is symmetrical at the origin; when the vertical vibration velocity ẏ is positive, the friction force is greater than zero, and when the vertical vibration velocity ẏ is negative, the friction force is less than zero.

Piecewise Stiffness Model
Based on the stress characteristics of the roll, when the external disturbance force F1 in the horizontal direction is high, and there is a gap between the roller bearing seat and the housing, the work roll bearing seat will collide with the mill housing, thus destabilizing the roll system [8].According to the elastic collision theory and the structural characteristics of the rolling mill, a horizontal collision model of the hot rolling mill rolls is established, as shown in Figure 6.As shown in Figure 6, x is the horizontal vibration displacement of the mill rolls, k1 and c1 are respectively the linear stiffness and linear damping from the roll to the left housing, k2 is the lateral equivalent stiffness of the housing column, and Δx is the gap between the roller bearing seat and the mill housing.
There is a gap between the roller bearing seat and the liner plate of the mill housing.Therefore, the horizontal stiffness can be modeled as a piecewise function using a force function Fk(x), as expressed in Equation ( 5); Figure 7 shows the corresponding function curve.

Determination of Rolling Parameters in Deformation Zone
During rolling, the rolling force is an important factor affecting the vibration of the rolling mill [25,26].It is important to determine the rolling force correctly to improve the accuracy of the rolling mill vibration model.Figure 8 shows the dynamic deformation process of the workpiece in the rolling process.x / mm . Dynamic deformation process of the workpiece.
In Figure 8, the solid line indicates the roll position in the steady state, and the dotted line indicates the roll position in the vibration state.R is the roll radius, vR is the roll speed, and α is the bite angle.v0 and v1 are the entry and exit velocities of the workpiece in the steady state, respectively, h0 and h1 are the entry and exit thicknesses of the workpiece in the steady state, respectively, h2 is the exit thickness of the workpiece in the vibration state, τb and τf are the entry and exit tensions, respectively.l1 and l2 are the distances from the solid line roller center line to the vibration state at the entrance and exit, respectively.
Assuming that the shape of the roll gap is a parabola [27], when vibration occurs, the section thickness of the workpiece at any position l of the roll gap is as follows: Due to the vibration, the position l2 at the exit of the workpiece is the vibration displacement x of the roller center in the horizontal direction.The exit thickness of the workpiece can be obtained from Equation (6).
When the entry thickness h0 of the workpiece is stable, from Equation ( 6), we can express the position l1 at the entrance of the workpiece under vibration as follows: When the rolling mill rolls vibrate, the entrance speed of the workpiece will be composed of two parts: the velocity v0 at the entrance of the workpiece in the steady state and the change rate of the position ̇ at the entrance of the workpiece under vibration.This can be expressed as follows: According to Hu's law of the invariance of dynamic volume flows [28]: where ẋ is the horizontal vibration velocity of the rolling mill rolls, vl is the velocity of the workpiece at any section l, and (l1 − l)ẏ and (h0 − hl)ẋ are the volume flow rates of the workpiece in the vertical and horizontal directions, respectively.
There is a section in the deformation zone, where the horizontal velocity of the workpiece is equal to that of the roller at a particular point; this section is called the neutral surface, and the point on the roller is called the neutral point.Therefore, we have vγ = vRcosγ at the neutral point lγ, where vγ is the horizontal velocity of the workpiece at the neutral point, and γ is the neutral angle.From the geometric relationship shown in Figure 7, we find that the cross-sectional thickness of the workpiece is hγ = h2 + 2R(1 − cosγ) at the neutral point.Since the bite angle α is small, and γ < α, it can be approximately considered that cosγ ≈ 1, (1 -cosγ) ≈ γ 2 /2, and hγ = h2 + Rγ 2 .By substituting hγ and vγ into Equation (10), we can express the neutral point position lγ as follows: ℎ ln ℎ ℎ

Determination of Vertical-Horizontal Coupling Dynamic Rolling Force
Assuming that the workpiece is uniformly deformed in the deformation zone and considering that the sign of the friction force between the rolling interfaces changes at the neutral point lγ, any microelement with the same thickness is selected for analysis in the front and backward sliding areas of the rolling deformation area.Its thickness is dl, under the action of the tensile stress τx, rolling compressive stress p, and shear stress τs; Figure 9 shows the stress relation diagram.Based on the force balance theory provided by Kármán for the deformation zone [29], the force balance distribution equation along the x direction can be obtained as follows: where the positive and negative signs in "±" indicate the forward and backward sliding areas, respectively.The Coulomb friction model is used to describe the friction between the rolling interfaces, i.e., τs = μxδ, where μx is the friction coefficient between the rolling interfaces, and δ is the shear yield strength.During rolling, the tensile stress at any point l can be written as: According to von Mises yield criterion p = 2δ − τx, the unit rolling force can be obtained as follows: At this time, the force Fx of the workpiece in the horizontal direction and Fy in the vertical direction can be expressed as − − ( − ) Since the direction of friction between the rolling interfaces changes at the neutral point, the integral Equation ( 15) can be divided into two parts, namely l2~lγ and lγ~l1, and can be rewritten as: Equation ( 16) describes the dynamic rolling force of the hot rolling mill rolls in the horizontal and vertical directions.The expression is complex.To analyze the influence of roll vibration on the rolling mill rolls, the rolling force formula can be expanded under a steady rolling state, i.e., Taylor expansion near x = y = 0, ẋ = ẏ = 0. Thus, Equation ( 16) can be rewritten as: ( , , , ) = (0,0,0,0) + Δ ( , , , ) ( , , , ̇) = (0,0,0,0) + Δ ( , , , ) where F(0,0,0,0) is the rolling force in the steady state, and Δ ( , , , ) is the dynamic variation in the rolling force under vibration.
(0,0,0,0) = Because different rolling mills and steel grades have different process parameters, some of the parameters of the Taylor expansion are sensitive to mill vibration, while some can be ignored.The parameters and their specific gravity values are shown in Table A1 and Table A2 (Appendix A).Taking the actual parameters of a 1780 mm hot rolling mill (listed in Table 1) as an example, the specific gravity values of the various parameters in the dynamic rolling force variation are calculated and compared, so as to select the appropriate parameters.

Coupling Vibration Modeling and System Response Solution of Hot Rolling Mill Rolls
The piecewise nonlinear spring force constraint and the piecewise nonlinear friction force constraint of the hydraulic system must be considered in the production process of high-precision rolling mills.Moreover, the piecewise stiffness constraint in the horizontal direction and the influence of the dynamic rolling force on the rolling mill rolls should be considered.In this study, a vertical-horizontal coupling dynamic model of the hot rolling mill under multi-piecewise nonlinear constraints is established as shown in Figure 10.In Figure 10, only the upper part of the mill rolls is shown; this is because many researchers assume that the mechanical structure of rolling mill rolls is symmetrical on and off the rolling line when exploring the vibration of four-high hot rolling mills.Therefore, when constructing the vibration model, they only build the upper part of the rolling mill rolls, which helps reduce the calculation workload and does not affect the accuracy of the vibration model [30].
In Equation ( 19), the piecewise nonlinear spring force and the piecewise nonlinear friction force are given coefficients α and β, respectively, to study the vibration behavior of the hot rolling mill rolls under different constraint conditions.After simplifying Equation (19), we can obtain the vibration equation of the hot rolling mill rolls as follows: In Equation ( 20), the expression of ( ) is shown in Equation ( 21), and the expression of ( ) is shown in Equation (22).
The horizontal vibration equation contains the vertical vibration displacement y and velocity ẏ, and the vertical vibration equation also contains the horizontal vibration displacement x and velocity ẋ, which indicates a coupling between the horizontal and vertical vibrations in the hot rolling mill rolls.
Since the coupled vibration system is weakly nonlinear, the nonlinear term in the system is given a small parameter ε.Moreover, it is assumed that the external excitation frequencies in the horizontal and vertical directions are close to the natural frequencies of the system, i.e., = − , and = − , where σ1 and σ2 are the tuning parameters.When ε is sufficiently small, the motion of the system is closer to a periodic motion; then, Equation ( 20) can be expressed as: The average method is used to solve the system response [31].Let Equation ( 24) have the following solution form: Considering that the frequency of the external excitation in the horizontal and vertical directions is the same, there exists ω1 = ω2 = ω.Moreover, when ̇= ̇= ̇ = ̇ = 0 in Equation ( 26), the coupling system has a stable vibration amplitude and frequency.From Equation ( 27), we find that the amplitude and frequency of the external disturbance, the nonlinear stiffness of the hydraulic cylinder, the nonlinear friction force, and the clearance between the bearing pedestal and the archway affect the vibration of the rolling mill.

Analysis of Vibration Characteristics of Coupling System of Hot Rolling Mill Rolls
Table 1 lists the parameters of the 1780 mm hot rolling mill of Chengde Steel Co., Ltd.(Chengde, China) as an example.With the relevant data listed in Table 1, the nonlinear parameters of the dynamic rolling force can be obtained using Equation (18), as listed in Table 2.

Analysis of Amplitude-Frequency Characteristics
In the following, we mainly study the influences of different parameters on the amplitude-frequency characteristics of the coupled vibration system of the hot rolling mill.When the horizontal and vertical directions are subjected to the same external disturbance, i.e., F1 = F2 = F, the amplitude-frequency curve of the vertical-horizontal nonlinear coupled vibration system of the hot rolling mill is as shown in Figure 11.There are two resonance peak points in the amplitude-frequency curve in the horizontal and vertical directions.When the external disturbance frequency ω is close to the natural frequency in the horizontal or vertical directions, the coupling system resonates, causing the vibration amplitude to increase and the rolling mill to vibrate violently.The vibration amplitude in the horizontal direction is greater than that in the vertical direction, and the resonance frequency range is 65-155 Hz.In actual production, we can adjust the relevant vibration and process parameters to reduce the amplitude, avoid the resonance region, and suppress and avoid the occurrence of resonance phenomenon.With further increase in the external disturbance force, when F = 4 × 10 5 N, the amplitude-frequency curve of the coupled vibration system of the hot rolling mill is as shown in Figure 12.Compared with that shown in Figure 10 (F = 2 × 10 5 N), the vibration amplitudes in both the horizontal and vertical directions increase, and the resonance becomes more intense.A jump phenomenon is observed in one of the branches of the resonance peak, making system unstable.Taking the vertical direction as an example (Figure 12b), when the external disturbance frequency increases, the amplitude changes along the AB direction of the curve to point C. When the frequency exceeds point C, the amplitude jumps from C to point D, and the amplitude jumps.When the external disturbance frequency decreases, the amplitude changes along the ED direction to point F, jumps to point B, and then changes along the BA direction.Therefore, the vibration of the coupling system corresponding to the CF section of the amplitude-frequency curve is unstable.The horizontal direction (Figure 12a) also has similar vibration characteristics.Compared with that shown in Figures 11 and 12, it is found that appropriately reducing the amplitude of the external disturbance force can help optimize the control parameters of the vibration system, thus suppressing the vibration of the rolling mill and improving the stability of the system.When Δx = 0 mm, the amplitude-frequency characteristic curve of the coupling system of the hot rolling mill is as shown in Figure 13 (gap Δx = 1 mm in Figure 11) with the decrease in the gap Δx between the roller bearing seat and the mill housing in the horizontal direction.The vibration amplitudes in the horizontal and vertical directions are significantly reduced, the resonance region is narrowed, and there is only one resonance peak point.This shows that eliminating the gap between the bearing seat and the mill housing can help significantly reduce the vibration of the system and improve the stability of the vibration system.The main reason is that when the gap becomes smaller, the stiffness of the hot rolling rolls varies in the horizontal direction, which changes the natural frequency of the hot rolling mill rolls.In practice, to eliminate the clearance, a hydraulic liner is often added between the bearing seat and the housing of the hot rolling mill rolls.Figure 14 shows the amplitude-frequency characteristic curve of the coupled vibration system of the hot rolling mill with piecewise nonlinear spring force constraint coefficient α = 1.5.Compared with that shown in Figure 11 (α = 1 in Figure 11), when the nonlinear spring force increases slightly, the amplitudes in the horizontal and vertical directions increase significantly, one branch of the vibration curve bends and shifts to the right, and the resonance region evidently widens (65-190 Hz).In the vertical vibration curve (Figure 14b), we find evident jumping phenomena at piecewise nodes e1 and e2, the horizontal direction (Figure 14a) also has similar jumping characteristics.Which indicates that the roll coupling system of the hot rolling mill rolls is very sensitive to the piecewise nonlinear elastic constraint of the hydraulic cylinder, because of which the hot rolling mill is often unstable in the rolling process.Figure 15 shows the amplitude-frequency characteristic curve of the coupled vibration system of the hot rolling mill with piecewise nonlinear friction force constraint coefficient β = 5.Compared with that shown in Figure 11 (β = 1 in Figure 11), increasing the nonlinear friction force constraint is equivalent to increasing the damping of the rolling mill rolls, reducing the amplitude of the vibration system, and reducing the resonance area (54-110 Hz).Therefore, appropriately increasing the nonlinear friction force can help effectively suppress the vibration of the system.Another reason for rolling mill instability is the influence of the signum function in the cubic model, which has continuous and discontinuous characteristics.Figures [11][12][13][14][15] show that the nonlinear stiffness, nonlinear friction, gap between the bearing seat and the mill housing, and external disturbing force affect the amplitude-frequency curve variation of the coupled vibration system of the hot rolling mill.Therefore, the resonance behavior of the coupling system can be tailored by selecting appropriate system parameters to reduce the resonance damage on the rolling mill rolls.

Analysis of Bifurcation Characteristics
The singularity theory can be applied to identify the parameters that reflect the structural stability from a large number of parameters of a vibration system, so as to grasp the overall vibration characteristics of the system and predict the stability.Therefore, we used the singularity theory to analyze the static bifurcation characteristics of the coupled vibration system of the hot rolling mill under different parameters.
After the calculation, the amplitude a in Equation ( 27) is eliminated, and a Taylor series expansion is carried out at b = b0, with u = b 2 .After the linear transformation, the static bifurcation equation of the vertical horizontal coupling vibration system of the hot rolling mill rolls can be obtained.
This bifurcation equation is a universal unfolding of GS normal form u 4 +λ, and the other dimensions are set to 3. Here, λ is the bifurcation parameter, and d1, d2, and d3 are the unfolding parameters.Since it is difficult to express the transition set of the coupled vibration system of the hot rolling mill directly in a 3D space, the static bifurcation characteristics on the 2D projection plane in different parameter spaces are discussed in the following three cases.
(1) When d1 = 0, the bifurcation equation is u 4 +λ+d2u 2 +d3u 3 = 0, bifurcation point set B = ϕ (ϕ is the empty set), hysteretic point set = = , double limit point set D = {d2 > 0, d3 = 0}, and transition set Σ = B∪H∪D.Figure 16 shows the transition set of the coupled vibration system of the hot rolling mill rolls.Figure 16 shows the corresponding topological structure in the different regions.In Figure 16, the transition set divides the plane composed of unfolding parameters d3-d2 into six regions (I-VI).In the same region, the topological structure is equivalent, i.e., the bifurcation morphology is similar; however, in different regions, the topological structure is not equivalent.Figure 17 shows nine bifurcation modes, which reflect all the information of the coupled vibration system when the unfolding parameter d1 = 0. On the curves of zones I, II, IV, and double limit point set D, the system is stable, and there is no jump phenomenon.On the curves of regions III, V, VI, and set of hysteretic points H, a bifurcation parameter λ corresponds to multiple u, and there is a jump phenomenon.The main reason for the jump phenomenon is the nonlinear elastic force in the vertical direction, nonlinear friction force, and sectional stiffness in the horizontal direction.set Σ = B∪H∪D.Figure 18 shows the transition set of the coupled vibration system of the hot rolling mill rolls.Figure 19 shows the corresponding topological structure.The transition set divides the plane composed of the unfolding parameters d3-d1 into six regions; on the curves of regions II, VI, and set of hysteretic points H, the system jumps, and other regions are in the steady state.show that the vertical horizontal coupled vibration system of the hot rolling mill under multi-piecewise nonlinear constraints has a rich static bifurcation behavior, with different bifurcation characteristics under the influence of different parameters.Through the analysis of the bifurcation behavior, the stable region and unsteady parameter region of the coupling system of the hot rolling mill can be determined.Moreover, changing the bifurcation parameter λ can make the coupling vibration system to exhibit different bifurcation states, thus changing the vibration behavior of the system.This provides a theoretical reference for restraining the vibration of the hot rolling mill rolls.

Conclusions
1. To improve the accuracy of the vibration model of a hot rolling mill, the piecewise nonlinear spring force and piecewise nonlinear friction force were used to describe the constraint of the hydraulic screwdown system in the vertical direction of the rolling mill.A piecewise stiffness model was used to describe the impact of the horizontal bearing seat of the rolling mill rolls on the housing.Moreover, we derived an expression for the nonlinear coupling dynamic rolling force, which was mainly affected by the vertical vibration displacement and velocity and horizontal vibration displacement and speed.On this basis, a vertical-horizontal coupling vibration model of the hot rolling mill rolls under multisegment nonlinear constraints was established.
2. By comparing the effects of different parameters on the amplitude-frequency characteristics of the coupled vibration system of the hot rolling mill, we found an evident jump phenomenon at the piecewise nodes of the nonlinear elastic force, and the greater the nonlinear elastic force, the greater the vibration amplitude and the wider the resonance region.An increase in the external disturbance force also caused the jump phenomenon, thus destabilizing the vibration system of the hot rolling mill.Reducing the gap between the roller bearing seat and the mill housing and appropriately increasing the nonlinear friction force could help better reduce the amplitude and weaken the influence of piecewise characteristics to suppress the vibration of the hot rolling mill rolls.When rolling a sheet metal with a width of 1500 mm and a thickness of 14.1 mm, through certain measures, we could decrease the horizontal and vertical external disturbance forces to less than 3.5 × 10 5 N, reduce the gap between the horizontal roller bearing seat and the mill housing to less than 0.5 mm, and increase the friction coefficient between the piston and the cylinder wall of the hydraulic cylinder to greater than 0.6, thus ensuring an efficient and stable operation of the rolling mill.
3. A static bifurcation equation for the nonlinear coupled vibration system of the hot rolling mill rolls was obtained using the singularity theory.By analyzing three different transition sets and their corresponding bifurcation topological structures, we found that the coupled vibration system exhibits varying static bifurcation behaviors, indicating a relatively complex vibration behavior of the coupled system under the influence of multipiecewise nonlinear factors.Finally, the steady and unsteady process parameters of the hot rolling mill rolls were given, and the stability of the coupled vibration system of the hot rolling mill was predicted from a macroscopic viewpoint.
Author Contributions: Conceptualization, data curation, writing-original draft, R.P.; Supervision, writing-review and editing, X.Z.; Software, project administration, writing-review and editing, P.S.All authors have read and agreed to the published version of the manuscript.

Figure 1
Figure1shows a coupling dynamics diagram of a rolling mill.As shown, in the production process, a change in the dynamic parameters of the rolling mill structure model will lead to changes in the process parameters, force, and energy parameters of the rolling process model; the force and energy parameters in turn cause vibrations in the rolling mill structure, which interact with each other and lead to a coupling vibration.

Figure 1 .
Figure 1.Coupling dynamics diagram of a rolling mill structure and rolling process.

Figure 2 .
Figure 2. Structure diagram of a four-high hot rolling mill.

Figure 3 .
Figure 3. Piecewise nonlinear spring force model of hot rolling mill rolls in the vertical vibration direction.

Figure 4 .
Figure 4. Piecewise nonlinear spring force of hot rolling mill rolls in the vertical direction.

Figure 5 .
Figure 5. Piecewise nonlinear friction force of hot rolling mill rolls in the vertical direction.

Figure 6 .
Figure 6.Collision model between the roller bearing seat and the mill housing in the horizontal direction of hot rolling mill rolls.(a) The roller bearing seat is pressed against the right mill housing; (b) The roller bearing seat is between the mill housing; (c) The roller bearing seat is pressed against the left mill housing.

Figure 7 .
Figure 7. Piecewise stiffness force function of hot rolling mill rolls in the horizontal direction.

Figure 9 .
Figure 9. Stress relation diagram at any position in the deformation zone.

Figure 10 .
Figure 10.Vertical-horizontal coupling vibration model of hot rolling mill rolls under multi-piecewise nonlinear constraints.
where a, b, θ1, and θ2 are slowly varying functions of time t.Through derivation, we can obtain the differential equations of a, θ1, b, and θ2, then integrate them in the period (−π, π), and finally calculate the average value.The average equation is as follows:

Figure 11 .
Figure 11.Amplitude-frequency curve of the coupled vibration system of a hot rolling mill.(a) Amplitude-frequency curve of horizontal direction; (b) Amplitude-frequency curve of vertical direction.

Figure 12 .
Figure 12.Amplitude-frequency curve of the coupled vibration system of the hot rolling mill with external disturbing force F = 4 × 10 5 N. (a) Amplitude-frequency curve of horizontal direction; (b) Amplitude-frequency curve of vertical direction.

Figure 13 .
Figure 13.Amplitude-frequency curve of the coupled vibration system of hot rolling mill with gap Δx = 0 mm.(a) Amplitude-frequency curve of horizontal direction; (b) Amplitude-frequency curve of vertical direction.

Figure 14 .
Figure 14.Amplitude-frequency curve of the coupled vibration system of hot rolling mill with piecewise nonlinear spring force constraint coefficient α = 1.5.(a) Amplitude-frequency curve of horizontal direction; (b) Amplitude-frequency curve of vertical direction.

Figure 15 .
Figure 15.Amplitude-frequency curve of the coupled vibration system of hot rolling mill rolls with piecewise nonlinear friction force coefficient β = 5.(a) Amplitude-frequency curve of horizontal direction; (b) Amplitude-frequency curve of vertical direction.

Figure 16 .
Figure 16.Transition set of nonlinear coupled vibration system of hot rolling mill rolls with d1 = 0.

Figure 17 .
Figure 17.Topological structure of different regions in the transition set when d1 = 0.

Figure 18 .Figure 19 .
Figure 18.Transition set of a nonlinear coupled vibration system of hot rolling mill rolls with d2 = 0.

( 3 )
When d3 = 0, the bifurcation equation is u4 + λ + d1u + d2u 2 = 0, bifurcation point set B = ϕ, hysteretic point set = = − , double limit point set D = {d1 = 0, d2 ≤ 0},and transition set Σ = B∪H∪D.Figure20shows the transition set of the coupled vibration system.Figure21shows the corresponding bifurcation topology.The unfolding plane is divided into four regions, and eight bifurcation modes are given.The rolling mill is in a stable state on the curve of origin O, I, IV, and double limit point set D; in the other regions, the jump phenomenon appears.

Figure 20 .
Figure 20.Transition set of nonlinear coupled vibration system of hot rolling mill rolls with d3 = 0.

Figure 21 .
Figure 21.Topological structure diagram of different regions in the transfer set when d3 = 0.
Figures[16][17][18][19][20][21] show that the vertical horizontal coupled vibration system of the hot rolling mill under multi-piecewise nonlinear constraints has a rich static bifurcation behavior, with different bifurcation characteristics under the influence of different parameters.Through the analysis of the bifurcation behavior, the stable region and unsteady parameter region of the coupling system of the hot rolling mill can be determined.Moreover, changing the bifurcation parameter λ can make the coupling vibration system to exhibit different bifurcation states, thus changing the vibration behavior of the system.This provides a theoretical reference for restraining the vibration of the hot rolling mill rolls.

Table 1 .
Parameters of 1780 mm hot rolling mill.

Table 2 .
Nonlinear parameters of the dynamic rolling force model.

Table A1 .
Board Statement: Not applicable.This research was funded by the National Natural Science Foundation of China, grant number 61973262, and Natural Science Foundation of Hebei Province, grant number E2019203146, and the Graduate Innovation Funding Project of Hebei Province, grant number CXZZBS2021136.Conflicts of Interest:The authors declare no conflict of interest.Taylor series expansion of the dynamic rolling force in the horizontal direction.
Informed Consent Statement: Not applicable.Data Availability Statement: Not applicable.Funding:

Table A2 .
Taylor series expansion of the dynamic rolling force in the vertical direction.