Bifurcation and chaos for electromechanical coupling torsional vibration in rolling mill system driven by DC motor

Haoran Liu1, Jinjie Liu2, Fei Liu3, Zhanlong Zhu4, Kun Wang5, Shuang Liu6 1, 3School of Information Science and Engineering, Yanshan University, Qinhuangdao, China 1, 3The Key Laboratory for Special Fiber and Fiber Sensor of Hebei Province, Yanshan University, Qinhuangdao, China 2, 4, 6School of Electrical Engineering, Yanshan University, Qinhuangdao, China 5School of Science, Yanshan University, Qinhuangdao, China 2Corresponding author E-mail: 1liu.haoran@ysu.edu.cn, 2xiaoyuliujinjie@126.com, 3liufei@ysu.edu.cn, 4657590052@qq.com, 5791698058@qq.com, 6shliu@ysu.edu.cn


Introduction
With the rapid development of the industrial technology, the rolling mill equipment continues melting towards maximization, high speed, automation and high accuracy directions.Torsional vibration of the rolling mill main drive system will serious affect the normal work and cause the damage to the equipment [1][2][3].Torsional vibration is usually caused by nonlinear dynamic behaviors such as bifurcation and chaos.So there is great practice and immediate significance to study the nonlinear dynamic characteristics of the electromechanical coupling system.
The torsional vibration of rotational system was simulated and tested based on different work theory, which provides the experimental basis for torsional vibration analysis and control [4][5][6].The low frequency torsional gap in shafts provides a new idea for vibration control, the vibration suppression in a lumped torsional system using wave-absorption control with online computation of an imaginary wave-propagation system was discussed, and the theoretical analysis showed that smaller mass and larger spring stiffness of an imaginary system absorb vibration energy better [7].The propagation of torsional wave in the shaft with periodically attached local resonators was studied and the effect of shaft material on the vibration attenuation in band gap was investigated in Ref. [8].Considering the effect of the backlash nonlinear factor of the relative rotation nonlinear dynamic system, the high-dimensional torsional vibration global dynamical equation was established based on Lagrange equation, the bifurcation characteristic was analyzed and the effects of the parameters on the vibration were discussed [9].Abstracting the rolling mill main drive system into two masses or multi-mass relative rotation system, the bifurcation and the chaos which may cause torsional vibration of system were studied [10,11].A nonlinear feedback controller was also proposed to control the Hopf bifurcation point and the stability of periodic solutions.Gustavsson R. K. et al. [12] analyzed the effects of unbalance magnetic pull on the transverse vibration characteristics of motor rotor system.Taking high frequency converting current of inverter power switches into consideration, the electromechanical coupling dynamic model of high speed grinding system was established, and the physical mechanism of higher harmonic electromechanical coupling vibration was revealed in Ref. [13].Wang Xingyuan et al. [14,15] dealed with the existence of both Hopf bifurcation and topological horseshoe for a novel finance chaotic system.Kobayashi T. et al. [16] studied the effect of slot combination on acoustic noise from induction motors.Considering the high harmonic in electricity, Wu Huimin et al. [17] analyzed the effect of electromagnetic stiffness on the inherent frequency.Considering the influence of distribution of salient pole on air-gap permeance, model of electromechanical coupling torsional vibration for the rotor system under electromagnetic excitation was established by Xu Jinyou et al. [18].Taking the nonlinearity of the generator shaft and the interaction of mechanics and electrics in the generator sets into account, a transient model was obtained by combining Park equations and mechanics equations in et al. [19].The Hopf bifurcation, period-doubling bifurcation and chaos were investigated with nonlinear mode and Floquet theory.Considering the electromechanical devices with unsaturated magnetic circuit, the dynamic behavior of the system is analyzed by P. Woafo [20,21].
Up to now, there are few researches about the electrical parameters in rolling mill torsional vibration system driven by DC motor.In this Letter, a transmission system driven by DC motor is abstracted as two masses torsional vibration model.The nonlinear electromechanical coupling dynamical equation of transmission system is established considering the coupling effect of electrical parameters and mechanical parameters.The equivalent low dimensional bifurcation equation is obtained and the bifurcation characteristic is analyzed by using the singularity theory.Finally, numerical simulations were performed to study the influence of torsional rigidity and excitation amplitude on chaotic motion.These provide a theoretical basis for design and control of the electromechanical coupling transmission system which is widely used in practical engineering.

Electromechanical coupling torsional vibration system
The main drive system of rolling mill usually consists of motor, universal couplings and rollers.The method of lumped mass is adopted to abstract the rolling mill main drive system driven by DC synchronous motor into a two-mass relative rotation model.
Considering the rolling mill main drive system under the incentive of torsional rigidity and outside disturbance, the total kinetic and potential energy of the system can be described as: Damping force is: Generalized torque is: where ( = 1, 2) is generalized coordinate.Substituting Eq. (3), Eq. ( 4) into Eq.( 5), yields the generalized torque: Substituting Eq. ( 6), Eq. ( 7) into Lagrange equation: Yields: where and are inertia moments of the motor and mill roll, and ( = 1, 2) are rotational angle and the speed of the rotational angle of motor and mill roll, is linear torsional stiffness of the system, is structural damping of the system, is output torque, is the damping force.
During the functioning of the devices, nonlinearity appears in the components of the mechanical parts.Indeed, for hardening torsion or compression effects in mechanical problems, the stiffness is not constant but increases with the applied constraint.The coefficient of torsion or compression can be approximated by the relation: where is the stiffness for small torsion or compression, is the rotation angle displacement, is the coefficient of nonlinearity and ∂ is the normalization angle, the motor has two parts: a stator or inductor in a permanent magnet and a rotor, mobile around a revolution axis.To avoid the Foucault effect and to decrease the hysteresis loss, the following assumptions are considered: symmetrical magnetic circuit, negligible magnetic hysteresis, negligible magnetic flux leak [22].Using Kirchhoff's law, the equation of the electrical part is given by the following relation: where, in the right, the first term is the ohmic voltage, the second term is the voltage across the inductor, and the third term ⁄ represents the coupling term between the electrical and mechanical parts.Using the Newton second law of dynamics for rotational motions and taking into account the Laplace force, the equation of the mechanical part is obtained as: where is the moment of inertia of the rotor, = − ⁄ is the torque force due to viscous frictions, is the torque constant, is the electromagnetic torque due to the Laplace force, is the resistive torque, it is modeled by the relation: When = and = the system equations based on DC drive is: Eq. ( 14) is the nonlinear dynamic equation considering the torsional rigidity and the coupling term between the electrical and mechanical parts.In the process of motor running, is influenced by the motor load, the power supply voltage and the current.
is selected as a major parameter to analyze the system dynamic characteristics in the following.

Bifurcation analysis of high dimensional electromechanical coupling system
If we want to analyze system characteristics the first step is dimensionless the data processing.In this paper, we use linear dimensionless method for data processing, assuming that the original data and dimensionless data have a linear relationship: where is the dimensional scale, according to the practical parameters of electromechanical coupling torsional vibration system for 2150 rolling mill in a certain factory, a dimensionless form of linear proportional method is obtained by letting: where Ω ⁄ is system dimensionless ratio value, is standard reduced value for the system angle.Considering the excitation amplitude of is added at the load equipment, Then Eq. ( 14) can be reduced order for first-order equation as follows:

Lyapunov Schmidt reduction
Eq. ( 15) can be expressed as: where = ( , , , , ) is the state variables of dynamic system and ( ) = ( , , , , ) are the corresponding nonlinear vector functions.At the equilibrium point , system Eq.( 16) satisfies: Let ( ) be the Jacobian matrix of ( ) at the equilibrium: The Lyapunov stability is decided by the eigenvalues of ( ).The Jacobian matrix ( ) can be expressed as: This is a 5×5 square matrix and the value of the determinant is: There is a zero eigenvalue when + 3 − = 0, and is a singular matrix.It can be known that dim ( ) = 1, let be a basis vector of null space ( ), then we can get = (− ⁄ , 0,1,0,0).The adjoint operator of is the conjugate transpose matrix * .Because is a real matrix and the conjugate transpose matrix * is also a real matrix, and there is * = .We can see that dim ( * ) = 1 since is a Fredholm operator with zero index.One basis vector: of dim ( * ) = 1 is obtained by calculation, where = + 3 .So the five-dimensional Euclid space can be decomposed as follows: = ( ) ⊕ = ( ) ⊕ ( * ).
Define the projection operator: and the complementary projection operator: then the equation: equivalent to: where ∈ ( ) = span( ), ∈ .According to the definition of inner product in Euclid space, ( − ) ( + , ) = 0 is equivalent to: Thus, the following equation is obtained: Since + 3 − = 0, rewritten this equation and substitute it into Eq.( 24), then the expression of can be got: where: Assuming that a new state variable is = , then the bifurcation equation can be written as: where , as the open fold parameters of bifurcation equation are both for the physical parameters of the system.

Singularity analysis
The topology of Eq. ( 26 IV: Transition set: Transition set ∑ is shown in Fig. 1, the unfolding parametric space − is divided into four subregions (I, II, III, IV) by ∑ and Eqs. ( 5)-( 9) are bifurcation response curves of ∑ .Each subregion and its boundary bifurcation are shown in Fig. 2, Through the analysis of Fig. 2 we can draw the conclusion that: the roots of the Eq. ( 26) are smooth distribution in subregion I, the roots of the Eq. ( 26) have amplitude jump in subregion 2-9, This means that the system is not stable and such situations should be avoided in system designing and controlling.

Numerical simulations of chaotic processes
As harmonic excitation often occurs torsional vibration in the process of the motor operation, and the harmonic excitation is not a fixed value, but cyclical changes over time, by Fourier transform it generally expressed as a linear combination of trigonometric functions (sine/cosine) or their integration.This paper adopts = cos( ) for harmonic excitation, it is the general expression.In order to study the dynamic chaos of electromechanical coupling system, Runge-Kutta method is used for numerical simulation of torsion vibration system.Then Eq. ( 15) can be expressed as:  The system phase diagrams and Poincare maps diagram under different are shown in Fig. 5. From Fig. 3, we can see that the bifurcation diagram changed when = 1.4,= 1.5 and = 1.52, system response is a stable periodic motion.It can be seen from Fig. 5(a)-(c) that the phase trajectory is a closed curve and the Poincare map has finite fixed points.From Fig. 4, the maximum Lyapunov exponent is greater than zero when = 1.58, it can be convincing of occurrence of chaotic motion.We can also see that phase trajectory repeatedly winding in enclosed area but not closed, and Poincare map has the obvious fractal structure from Fig. 5(c1), (c2).When = 1.56 and keeps the other parameters constant, let change in a wide range.From the bifurcation diagram (Fig. 6.) and the maximum Lyapunov exponent (Fig. 7.) of , it can be seen that system jumps into single periodic motion from chaotic motion with the increase of .The system phase diagrams and Poincare maps with different are shown in Fig. 8 From Fig. 7, the maximum Lyapunov exponent is greater than zero when = 1.6, and the system appears chaotic motion.We can also see that phase trajectory repeatedly winding in enclosed area but not closed, and the Poincare map has the obvious fractal structure from Fig. 8(a).From Fig. 7, the maximum Lyapunov exponent is lower than zero when = 1.63, = 1.65 and = 1.7, system response appears a stable periodic motion.Then Fig. 8(b)-(d) show that phase trajectory is a closed curve, and Poincare map for a finite number of fixed points.
According to the numerical simulation, we come to a conclusion that both reducing torsional rigidity and increasing excitation amplitude can control the system from chaotic motion into stable periodic motion.
In order to control the chaotic motion, the external torque is introduced, assuming that kinetic equation with external torque is: The original chaotic motion will be changed and the new stable periodic motion will appear when we adjust the external torque .The external torque must be kept in non-zero value after the system movement is driven to the stable periodic motion, at this moment we can achieve the goal of chaos control.The value range of can be determined according to the Lyapunov stability theory and adjusting the Lyapunov index size.
In this paper the external torque is joined when the system appears chaotic motion, and we can select the appropriate torque force to control chaotic motion through analyzing the system bifurcation diagram and the maximum Lyapunov exponent when the torque changes.This method does not need to know the stable periodic orbit and state variables of the system.It only needs to load the system's motion to the periodic motion, by using the external torque, this method is simple and has little influence on the system's motor parameters.Contribution of each individual co-author to this article.Jinjie Liu -the dynamic model of the electromechanical coupling transmission system driven by DC motor is established, and the bifurcation characteristics of the system are analyzed by using the method of L-S reduction.Fei Liu -chaotic motion analysis of electromechanical coupling transmission system.Zhanlong Zhu -acquired data, and played an important role in interpreting the results.Kun Wang -helped perform the analysis with constructive discussions.Shuang Liu -revised the manuscript and approved the final version.

Conclusions
In this paper, considering the electromechanical devices with unsaturated magnetic circuit, a high dimensional dynamical equation of electromechanical coupling vibration in rolling mill drive system are deduced by using the dissipation Lagrange equation.The equivalent low-dimensional bifurcation equation which can reveal the influence of system parameters on the nonlinear dynamic characteristics is obtained by reducing the dimensionality of the system.The static bifurcation is studied by using the singularity theory.Considering outside incentive, choose and are chosen as bifurcation parameters, Numerical simulations are also given at last, which show that both reducing torsional rigidity and increasing excitation amplitude can control the system from chaotic motion into stable periodic motion, we can also add a constant torque direct feedback method to control chaos.These results have an important theoretical significance on reducing the vibration of electromechanical coupling transmission system.In addition, the research in our paper provides a reference for the real design of motor's parameters of electromechanical coupling transmission system.
) and the effects of the open fold parameters on the bifurcation behavior 2019.BIFURCATION AND CHAOS FOR ELECTROMECHANICAL COUPLING TORSIONAL VIBRATION IN ROLLING MILL SYSTEM DRIVEN BY DC MOTOR.HAORAN LIU, JINJIE LIU, FEI LIU, ZHANLONG ZHU, KUN WANG, SHUANG LIU will be studied by using the singularity theory.The following point sets are obtained based on the definition of the transition set: I: Bifurcation point set: