Abstract

The nonlinear torsional vibration and instability oscillation caused by nonlinear damping in the shearer electromechanical coupling cutting transmission system in shearer driven by the permanent magnet synchronous motor (PMSM) is investigated in this paper. The electromechanical coupling transmission system in the shearer is equivalent to a concentrated mass model for the purpose of establishing the system dynamic model by the Lagrange–Maxwell equation. Then, the Routh–Hurwitz criterion is used to determine the torsional vibration critical point and stability region for the Hopf bifurcation for the cutting transmission system. According to the Routh–Hurwitz stability criterion, the Hopf bifurcation type and the effect of the supercritical Hopf bifurcation in the torsional vibration of the cutting transmission system are analyzed. Furthermore, based on the washout filter, the Hopf bifurcation controller is designed for suppressing the transmission system’s large vibration amplitude and unstable oscillation. In addition, the influences of the linear gain and nonlinear gain on the bifurcation point and the limit cycle amplitude are discussed. Finally, the numerical simulation results indicate the effectiveness of the designed controller. The research achievements can provide a theoretical basis for design or optimize the cutting transmission system of high-reliability shearer driven by PMSM.

1. Introduction

In recent years, with the continuous development of coal mining mechanization, it requires higher requirements for the fully mechanized mining equipment in high efficiency and high reliability [1, 2]. As an integral part of fully mechanized mining equipment, the shearers may easily break down due to the terrible working environment. Currently, the traditional main transmission mode of the cutting transmission system is “three-phase asynchronous motor + three-stage straight tooth retarder + planetary gear reducer + cutting drum,” which has a long transmission chain and many transmission components. The long-chain transmission mode is prone to failure [35]. Therefore, a semidirect driving system [6, 7] is proposed in this paper, which consists of PMSM, three-stage straight tooth retarder, and cutting drum, and is shown in Figure 1.

Comparing to the original system, the electromagnetic torque and load torque will greatly increase due to the system reduction ratio that is reduced by removing the planetary gear [8, 9], which is prone to torsional vibration of the cutting transmission system. Therefore, it has a great significance to research the torsional vibration of the cutting transmission system. Fortunately, this multimass elastic torsional system like the shearer cutting section has been researched by many domestic and foreign scholars who mainly study from the perspectives of natural frequency and disturbance. It is found that the nonlinear factors are the main cause of torsional vibration, and furthermore, bifurcation and chaos may occur in the system due to some nonlinear factors [1012]. On the coupled bending and torsional vibration model of rubbing rotor, with rotational speed as the bifurcation parameter, the bifurcation and chaotic behavior of the rotor system are studied by Al-Bedoor [13]. Liu et al. [14] studied the static bifurcation characteristics and rotor system stability with time-delay nonlinear parameters, and then the vibration was controlled with the delay feedback. In [15], the self-excited oscillations caused by nonlinear factors in the power system are analyzed. And the system Hopf bifurcation motion was researched with nonlinear mode and Floquet theory with different resistance parameters. The mechanical rotor system Hopf bifurcation was researched by Yuan [16], which was frequently encountered in the mechanical structure, and may cause the shafting instability and damage. Saigo et al. [17] studied the suppression of forced vibration of the concentrated mass torsion system by a wave absorption. The results demonstrate that smaller mass and larger spring stiffness of an imaginary system absorb vibration energy better.

However, the cutting transmission system in shearer not only includes the motor but also mechanical parts, which belongs to the electromechanical coupling system. Due to the effect of electromechanical coupling, the cutting transmission system will exhibit complex dynamic behavior [1820]. Therefore, the torsional vibrations of the system cannot be interpreted rationally, neither from electrical nor from mechanical aspect. It is important to establish an accurate model of the electromechanical coupling system to further analyze the torsional vibration dynamic characteristics in the semidirect transmission system of the shearer and control the torsional vibration from the aspect of motor [21]. Gustavsson and Aidanpää [22] established the model of electromagnetic excitation produced by the uniform magnetic field, and the influence of nonlinear magnetic pull on rotor vibration of the hydrogenerator unit was discussed. Szolc et al. [23] researched the dynamic electromechanical coupling effects in machine drive systems driven by asynchronous motors, and the electromagnetic stiffness and damping coefficients of the interaction between motor and mechanical system were determined by means of the analytical method. Ju et al. [24] constructed the global electromechanical-coupling dynamic models for the main transmission system of the scraper conveyor, and the Hopf bifurcation characteristic is analyzed, which was caused by sliding friction. In view of the above mentioned, Hopf bifurcation is a significant factor affecting the stability of the electromechanical coupling transmission system. However, the dynamic analysis of the electromechanical coupling system mainly focus on the induction motor, and the PMSM are rarely considered. It is a remarkable fact that the damping of the semidirect transmission shafting would change with the shafting vibration, the shaft fatigue, the wear, and lubrication of the bearing [25, 26]. Torsional vibrations of the cutting transmission system caused by shaft damping variations would further lead to the shafting destabilization [27, 28].

Different from the above references, this paper analyses the effect of Hopf bifurcation produced by variable damping in shearer semidirect transmission system driven by the PMSM. Furthermore, based on the washout filter, the Hopf bifurcation controller is designed for suppressing the transmission system’s unstable oscillation. Finally, this paper is organized as follows. The electromechanical-coupling dynamics modelling for the main transmission system of the shearer driven by the high-torque and low-speed PMSM is provided in Section 2. The Hopf bifurcation characteristics analysis for the system torsional vibration and numerical simulation verification is given in Section 3. The washout filter controller for the system torsional vibration is designed, and numerical simulation verification is given in Section 4. Finally, the conclusions are drawn in Section 5.

2. Electromechanical Coupling Torsional Vibration Modelling

The semidirect cutting transmission system of the shearer driven by PMSM mainly consists of PMSM, three-stage gear, shafting, cutting drum, and so on. There are both resilient elements with small quality and inertia with large mass. Therefore, it is a multimass spring damping system with electromechanical coupling. To analyze the coupling dynamic characteristic of the electrical parameters and mechanical parameters, the PMSM and the cutting drum are regarded as the power source and load system, respectively, and ignore the effect of the reducer gear clearance. As a protective shaft for PMSM, the elastic torque shaft is more vulnerable to damage due to the relief groove and the low torsional strength. However, all other components have the feature of large mass and the small elasticity. Hence, the effect of the elastic torque shaft is only considered, and the other shaft and gear are equivalent to the PMSM and cutting drum. Then, the transmission system is equivalent to the two-mass spring torsional vibration system with resilient elements (PMSM and cutting drum) and inertia element (elastic torque shaft), as shown in Figure 2. To simplify the analysis, take the following assumptions: (1) neglecting the iron core saturation effect in the PMSM; (2) no consideration of eddy current and hysteresis loss in the PMSM; and (3) the motor stator winding current is a three-phase symmetric sine current.

In Figure 2, is the PMSM inertia; is the sum of the moment of the cutting drum and the transformation of the gear and other elements inertia; , , and are the three-phase stator current of the PMSM, respectively; and are the torsional stiffness and damping coefficient of the elastic torque shaft; and are the electromagnetic torque of the PMSM and load torque of cutting drum; and and are the mechanical angle of the PMSM shaft and cutting drum shaft.

Two mass systems are modelled by the Lagrange equation [29]. The kinetic energy of the shearer cutting transmission shafting mainly includes the PMSM rotational kinetic energy and the cutting drum rotational kinetic energy, which can be expressed as follows:

The system potential energy can be represented as

As a result of the downhole environment, water is sprayed to reduce the dust concentration, so the cutting environment is wet. In the process of cutting coal, system will produce nonlinear sliding friction forces due to relative displacement of cutting drum and coal seam. In view of the linear damping and nonlinear stick, slip frictional force is studied [26], which can be shown as follows:

The nonconservative generalized force of the electromagnetic system is electromagnetic torque and damping force. And the nonconservative generalized force of the mechanical part is the external load and damping force, and it are expressed aswhere is the generalized coordinate. Substituting equation (3) into equation (4) yields the generalized torque:

The system Lagrange–Maxwell equation is

For the mechanical system, the dynamic equation of the PMSM can be obtained by taking the rotation angle of the PMSM as the generalized coordinates () and substituting equations (1)∼(5) into equation (6):

Similarly, the dynamic equation of the cutting drum is

The mathematical model of PMSM in abc three-phase stationary coordinate system can be obtained [30] as follows:

Based on the electromagnetic field theory, electromagnetic torque is equal to the deflection of magnetic energy storage on the mechanical angular displacement . are the magnetic chain of three-phase winding. are the magnetic chain produced by a permanent magnet on the three-phase winding, and p is the number of pole pairs:

By introducing the Clark transform and Park transform, the stator current of the PMSM can be transformed from the three-phase stationary coordinate system (abc) to two-phase rotating coordinate system (dq). Equation (9) can be converted towhere p is the number of pole pairs, , id and iq are the direct and quadrature axis current in the coordinate system, and are the direct and quadrature axis inductor in the dq coordinate system and in stick-type PMSM, is the magnetic chain of the permanent magnets, and is the electric mechanical angular displacement.

Based on the electromagnetic field theory, in the two-phase rotating coordinate system (dq), the electromagnetic torque is

When the load torque and the electromagnetic torque are balanced, the mechanical angular velocity of the rotor is constant; therefore, , , , and [31], so

Finally, substituting equation (13) into equation (11) and combing with equations (7), (8), and (12), the electromechanical coupling torsional vibration model of the cutting transmission shafting is

In equation (13), are variable items which will change with time. This system will exhibit complex dynamical behavior due to the variation of the damping C.

3. Cutting Transmission System Bifurcation Analysis

For the convenience of the calculation, the difference value of two angles () is used to describe the strength of the vibration. Therefore, let the third type minus the fourth type in equation (14) merge nonvariable items in the system, defined as follows:where x1 is the quadrature axis current in the dq coordinate system, x2 is the direct axis current in the dq coordinate system, x3 is the difference between the mechanical angle of the PMSM and the angle of the cylinder, and x4 is x3’s derivative of time.

Therefore, equation (14) is simplified to

After the cutting transmission system of the shearer is changed to semidirect drive, the damping will change in the following cases: (1) the intensity of the vibration in the shaft system and (2) the wearing and the lubrication condition of the shaft. In addition, there may be oscillation vibration and Hopf bifurcation with the damping changes in the shafting, so the damping is selected as the bifurcation parameter.

The equilibrium point of equation (16) can be transferred to the coordinate origin through simple linear transformation. Therefore, the dynamic characteristics of the system balance at the origin are of universal significance without loss of generality. Hence, let , and then equilibrium point can be transferred to the coordinate origin for equation (16). The Jacobian matrix and nonlinear terms of the system at , respectively, are as follows:

According to the Hopf bifurcation theory, the necessary and sufficient conditions for Hopf bifurcation are expressed as follows: there are only one pair of conjugate imaginary roots in the characteristic root of and the other characteristic roots have negative real parts. The system is unstable at the origin when positive real part exists. While the eigenvalues have all negative real parts, the system is asymptotically stable. In order to simplify the calculation, the Hurwitz criterion is used to determine the existence of the Hopf bifurcation in this paper.

The characteristic polynomial of can be expressed aswhere , and .

Hurwitz criterion can be shown as

In equation (19), is the order principal minor determinant of equation (18):where the possible threshold is derived from . And if the resulting point satisfies the other requirements of equation (19) at the same time, the Hopf bifurcation occurs at this point. The parameters in the model are selected according to the data in the actual system shearer. Therefore, the main parameters of semidirect cutting transmission system in the form are as follows:where is the rated speed of PMSM, is the moment of inertia for PMSM, is the sum of the moment of the cutting drum and the transformation of the gear and other elements inertia, and is the stiffness coefficient of the elastic torque shaft. , , , and can be obtained, respectively, by off-line identification of PMSM.

By calculation, the couple transmission system of the shearer semidirect drive met equation (19) at . In this case, the characteristic root of has one pair of conjugate imaginary roots. The pure virtual root is , and vector and are right and left feature vectors corresponding to the characteristic root , respectively. Their relationship can be expressed as

Through equation (22), the stability index b of the Hopf bifurcation is defined as

In equation (23), is the conjugate plural of , and :

According to the value of , the stability of periodic solutions of Hopf bifurcation in the original system can be determined. If , the bifurcation type of the system is subcritical Hopf bifurcation, and the periodic solution is unstable. On the contrary, , the bifurcation type of the system is supercritical Hopf bifurcation, and the bifurcated periodic solution is stable. The stability index of Hopf bifurcation is through calculation, and therefore, the type of Hopf bifurcation is supercritical Hopf bifurcation.

The Runge–Kutta method is used to verify the dynamic characteristics generated by the system Hopf bifurcation point , and the parameters are brought into equation (16). In Figure 3 of torsional vibration amplitude figure variation with the damping , the system exhibits different dynamics on both sides of the bifurcation point because of supercritical Hopf bifurcation. When is greater than the bifurcation point , the running state of the main transmission system is stable and no limit cycle appears, as shown in Figure 4.

When is smaller than the bifurcation point , the system exhibits a stable limit cycle motion. Figure 3 shows that when the is less than the critical value , the smaller the value of the damping , the bigger the limit cycle. At the same time, when the damping is greater than the critical damping, the amplitude of the vibration is close to 0. It can be available from the comparison of time-history response of identical initial value in Figure 5, and the amplitude of the limit cycles is obviously larger in the case of smaller damping. Furthermore, when is fixed, the amplitude of the limit cycle has nothing to do with the initial value, as shown in Figure 6, which shows that the amplitudes are basically the same under two different initial values. From the above analysis, the conclusion of the simulation results is consistent with the theoretical analysis.

4. Cutting Transmission System Vibration Control

The numerical simulation results in Figures 3 and 5 show that, the system oscillates and the limit cycle appears when the damping C is smaller than . When the system responds to a smaller limit cycle, it has little harm to the elastic torque shaft. However, when the external conditions cause damping C to be smaller, the amplitude of the limit cycle will increase. At this time, the cutting transmission shaft system of the shearer will have strong vibration, which will cause great harm to the fatigue damage of the elastic torque shaft. In order to control the unstable oscillation of shafting, the washout filter controller is introduced, which is a kind of the extended linear or nonlinear state feedback method.

The state equation of the washout filter can be represented aswhere is the state variables of the washout filter, is the input variables, is the time constant, and is the output variable. In this paper, from the perspective of implementation, is used as the input variable, and is used as the state variables. In nonlinear systems with Hopf bifurcation, the squared term and the cubic term can affect the stability of Hopf bifurcation. So, nonlinear feedback controller can be shown as ; and are the linear gains and nonlinear control gain, respectively. is added to the second term in equation (16) as the output feedback.

So, the controlled system can be shown as

The Jacobian matrix and nonlinear terms of the control system can be represented as

If , the washout filter will not change the position of the equilibrium point. And the characteristic equation of in equation (26) iswhere .

If all eigenvalues of matrix have the negative real part, the controlled system will not occur Hopf bifurcation at , which means . When , , and , the roots of have negative real parts at this time, and the system is asymptotically stable, as shown in Figure 7(b). At this moment, the bifurcation characteristics are controlled at the Hopf bifurcation point , but the other dynamics of the system have not been determined.

According to the Hopf bifurcation theory, when Hopf bifurcation occurs at the equilibrium point, system (26) must satisfy the condition of equation (29). The relationship between the linear control gain and the bifurcation critical point can be drawn, which is shown in Figure 8. When , the controlled system has Hopf bifurcation same as the original system at . And if , the controlled system has the Hopf bifurcation point at . Comparing to the original system, the system stability domain increases greatly. Therefore, when the linear gain part of the controller selects the appropriate value, the bifurcation characteristics of the nonlinear system can be effectively controlled and the bifurcation point is transferred to the safe working area:

Take , and the Hopf bifurcation of the controlled system occurs at through calculation. The Jacobian matrix has a pair of pure imaginary roots , whose corresponding eigenvectors can be expressed as

Since it is more complicated to calculate the three-order normal form of the controlled system through the central manifold theorem, the direct calculation method of normal form in [32] can be used. This method only needs to calculate the pure virtual root of the controlled system and its corresponding eigenvector, and then the normal form of the controlled system can be obtained through the nonlinear transformation. The linear transformation is introduced here:

The Hopf bifurcation specification of the controlled system can be shown as

The nonlinear transformation is introduced to the nonlinear terms, and then the polynomial form of and as follows:where is coefficient vector corresponding to .

In the above equation,

Finally, the parameters E is

So normal form of the controlled system can be expressed as

According to the Hopf bifurcation theory, the controlled system will have a supercritical bifurcation at when . The bifurcation type of the system stays the same by calculating. At the same time, subcritical Hopf bifurcation will occur at , as shown in Figure 9. It can be seen from Figure 9(b) that the time-history response is stable when the initial values are close to the equilibrium point. When initial values are far from the equilibrium point, the time-history response increases with time in Figure 9(a). It is necessary to avoid the generation of subcritical bifurcation in the actual situation.

The amplitude of the limit cycle can be effectively controlled through the selection of the value of and reasonably. Figure 10 shows that the amplitude of the controlled system is obviously lower than that of the original system without subcritical Hopf bifurcation (Figure 10(b) is the time-history diagram of the controlled system at ), and the significant effect is achieved after addition of the control system. In contrast to the two controlled systems in Figure 11 ( in Figure 11(a); in Figure 11(b)), if the case of is fixed, the smaller the is, the smaller the amplitude of the limit cycle is.

The bifurcation point will vary with , and the bifurcation type will change with . At the same time, it is necessary to calculate the eigenvalue and eigenvector of the system to determine the type of bifurcation, but the eigenvalue and eigenvector are difficult to calculate due to the change of . So it is hard to calculate the type of bifurcation for the system through the theoretical method when and change simultaneously. However according to the simulation, the subcritical Hopf bifurcation only occurs when is near 0. Figure 12 shows the amplitude variation of and change at the same time when is far less than 0. It can be seen from Figure 12 that the vibration amplitude decreases with the reduction of linear gain and nonlinear gain .

5. Conclusion

In this paper, the Hopf bifurcation caused by the change of the system damping was studied concerning the shearer transmission shafting semidirect driven by the PMSM. The electromechanical coupling dynamic model was established by the Lagrange–Maxwell equation. The dynamic bifurcation behavior of the system transmission shafting was analyzed through theoretical analysis and numerical simulation. Based on the washout filter, the torsional vibration bifurcation of the shafting was effectively controlled. The research work of this article mainly includes the following:(1)Due to the harsh environment of coal cutting, there is the supercritical and subcritical Hopf bifurcation with the diversification of the damping in the cutting transmission system. When the is less than the critical value , the smaller the value of the damping , the bigger the limit cycle in supercritical Hopf bifurcation. Simultaneously, the amplitude of the limit cycle has nothing to do with the initial value.(2)Based on the normative theory, the Hopf bifurcation controller of the shearer cutting transmission shaft was designed by the washout filter. The results show that the system stability domain increases with the linear gain of the controller and the generation of the subcritical Hopf bifurcation can be inhibited by the cube nonlinear parts.(3)For supercritical Hopf bifurcation, selecting the linear and nonlinear gain reasonably can effectively control the torsional vibration and the shafting torsional vibration. And the limit cycle amplitude decreases with the nonlinear gain of the controller when the linear gain is zero. At the same time, when kn and km are varied simultaneously, the limit cycle amplitude decreases with the reduction of linear gain and nonlinear gain in supercritical Hopf bifurcation.

All these results are of important theoretical significance to reduce the torsional vibration caused by nonlinear factors in the shearer cutting transmission system and provide a guarantee for the safe and reliable application of the low-speed and large-torque PMSM on the shearer.

Data Availability

All data included in this study are available upon request by contact with the corresponding author.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 51775543), the Key Research and Development Project of Xuzhou (Grant no. KC17014), and the project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).