Abstract
Due to the nonlinear factors of the single loop gear system and the intertooth system, is sometimes difficult to establish an accurate nonlinear dynamic model, which leads to the large deviation between the dynamic characteristics and the actual situation. According to the structural characteristics and dynamic mechanism of nonlinear factors, the bond graph power junction with switching characteristics is adopted. This method is used to establish the model of the single loop gear system and the intertooth system respectively. On this basis, the amplitudephasefrequency characteristics of PX single loop gear systems are obtained by numerical simulation analysis. Eventually, the dynamic stability of single loop gear system is judged based on this characteristic index. The study shows that the dynamic stability of PX type single loop gear system is unsteady in the nonlinear state. In addition, it is concluded that when the whole single loop gear system is unstable, its system module also has the same unstability.
1. Introduction
Under the rapid development of science and technology, the requirements of mechanical equipment for the transmission system are getting higher and higher, such as multipower source input and output, the transmission system with power diversion, confluence and other functions [13]. A certain gear system can make the structure of the transmission system more complex, such as a singleloop gear system, a doubleloop gear system, or a multiloop gear system. With the increase of loops number, the motion and dynamic characteristics have become extremely complex, which makes it difficult to control the dynamic characteristics of the loop gear system [4, 5]. In addition to the intricate structure of the transmission system, the diversity of transmission forms, such as mechanical, electromagnetic and fluid transmission modes, which leads to the modeling of the dynamic characteristics is particularly complicated. However, the traditional modeling method is extremely difficult to establish mechanical models of various transmission modes, and it cannot guarantee the presentation of real dynamic characteristics among various systems [68].
In the transmission system, the gear transmission is the most basic component unit. The nonlinear factors of the intertooth system have significant influence on the dynamic characteristics, such as timevarying mesh stiffness, tooth surface sliding friction, meshing damping and steadystate transmission error [912]. Considering the influence of the above nonlinear factors, the dynamic characteristics of the system are more similar to the actual transmission one. However, the dynamic modeling and solution of the transmission system has also brought many difficulties owing to the introduction of nonlinear factors. Usually, the springmassdamping element method and the lumpedparameter method are used to establish the mathematical model of the intertooth system [13, 14]. Because these methods adopt equivalent modeling, the overall influence of nonlinear factors in the system is ignored, so that the calculation results have a certain narrow sense. In these aspects, a lot of achievements have been accumulated. For example, the influence of flexible support added to a large fan tower on the dynamic characteristics of the wind motor gear transmission system. The nonlinear dynamic model including tower stiffness, timevarying meshing stiffness, meshing damping, backlash and other factors has been established using the vibration mechanics. Finally, the vibration characteristics of this system have been achieved through calculation and analysis [15, 16]. Some scholars have used the lumpedparameter method to establish a nonlinear dynamic model of wind motor gear transmission system including timevarying meshing stiffness, composite meshing error and backlash, and obtained the response of nonlinear factors to the system through calculation and analysis. Scholars have considered the nonlinear factors such as toothside clearance, timevarying engagement stiffness and comprehensive error of inclined gear, and established a coupled nonlinear dynamic model with three degrees of freedom. Based on the "piecewise linear" backlash function, the nonlinear dynamic response results under different backlash functions are calculated by numerical method [17, 18]. There are plenty of research results in this field, but these research methods all solve the macroscopic dynamic characteristics of the transmission system, and the authenticity of the intertooth system, the dynamic characteristics coupling between multiple energy domains is not deeply involved.
In view of the problems existing in the dynamic research of the complex multienergy transmission system, the bond graph method can solve them well [19, 20]. The bond graph theory is utilized to study the dynamic performance of systems in various engineering fields. Compared with other mechanical analysis methods, it has several unique features. Bond graph can be a natural development of block diagram. It can not only represent the signal flow direction of components in the system, but also indicate the power flow direction and the causal relationship of control signals. In addition, in the bond graph method, the causative relationship is determined by certain rules. The appearance and shape of the bonding diagram are basically the same for the same engineering system, so that it can display the composition characteristics of the transmission system very intuitively.
2. Structural model of the single loop gear system
As shown in Table 1, there is one typical structure of singleloop gear transmission system: PX type, which is composed of a group of singledegreeoffreedom fixed shaft gear train $P$ and twodegreeoffreedom turnover differential gear train $X$. ${J}_{0}$, ${J}_{1}$ refer to mechanical connection points; $I$ and $O$ are input and output terminals; Other characters represent the drive shaft.
Table 1Structural graph of single loop gear transmission system
Name  Structural drawing  Feature 
PX  Input shunt type 
3. Modeling principle of the intertooth nonlinear system
3.1. Timevarying mesh stiffness
3.1.1. Singletooth meshing stiffness
There are four kinds of energy stored in the gear teeth during meshing transmission: Hertz contact energy, radial compression deformation energy, bending potential energy and shear deformation energy. According to the calculation principle of energy method, combined with the force of the gear tooth shown in Figure 1, the calculation relationship of Hertz contact stiffness ${K}_{h}$, radial compression stiffness ${K}_{a}$, bending stiffness ${K}_{b}$ and shear stiffness ${K}_{s}$ of the gear tooth can be obtained respectively (the number of teeth is less than 41), as shown in Eq. (1). When the number of the gear teeth is greater than 41, remove the nonintegral term and change the upper limit from ${\alpha}_{2}$ to ${\alpha}_{4}$:
$\frac{1}{{K}_{a}}=\frac{({R}_{b}\mathrm{c}\mathrm{o}\mathrm{s}{\alpha}_{2}{R}_{f}\mathrm{c}\mathrm{o}\mathrm{s}{\alpha}_{3}){\mathrm{sin}}^{2}{\alpha}_{1}}{2EB{R}_{b}\mathrm{sin}{\alpha}_{2}}+{\int}_{\alpha 1}^{\alpha 2}\frac{\left(\alpha 2\alpha \right)\mathrm{c}\mathrm{o}\mathrm{s}\alpha {\mathrm{sin}}^{2}{\alpha}_{1}}{2EB[\mathrm{sin}\alpha +({\alpha}_{2}\alpha \left)\mathrm{cos}\alpha \right]}d\alpha ,$
$\begin{array}{l}\frac{1}{{K}_{b}}=\frac{{\left[1\frac{\left({R}_{b}\mathrm{c}\mathrm{o}\mathrm{s}{\alpha}_{2}{R}_{f}\mathrm{c}\mathrm{o}\mathrm{s}{\alpha}_{3}\right)\mathrm{c}\mathrm{o}\mathrm{s}{\alpha}_{1}\mathrm{c}\mathrm{o}\mathrm{s}{\alpha}_{3}}{{R}_{b}}\right]}^{3}(1\mathrm{c}\mathrm{o}\mathrm{s}{\alpha}_{1}\mathrm{c}\mathrm{o}\mathrm{s}{\alpha}_{2}{)}^{3}}{2EB\mathrm{c}\mathrm{o}\mathrm{s}{\alpha}_{1}\mathrm{s}\mathrm{i}{\mathrm{n}}^{3}{\alpha}_{2}}\\ +{\int}_{{\alpha}_{1}}^{{\alpha}_{2}}\frac{3(\alpha 2\alpha )\mathrm{c}\mathrm{o}\mathrm{s}\alpha \{1+\mathrm{c}\mathrm{o}\mathrm{s}{\alpha}_{1}[({\alpha}_{2}{\alpha}_{1})\mathrm{s}\mathrm{i}\mathrm{n}\alpha \mathrm{c}\mathrm{o}\mathrm{s}\alpha ]{\}}^{2}}{2EB[\mathrm{s}\mathrm{i}\mathrm{n}\alpha +({\alpha}_{2}\alpha )\mathrm{c}\mathrm{o}\mathrm{s}\alpha {]}^{3}}d\alpha ,\end{array}$
$\frac{1}{{K}_{s}}=\frac{1.2(1+\mu )({R}_{b}\mathrm{c}\mathrm{o}\mathrm{s}{\alpha}_{2}{R}_{f}\mathrm{c}\mathrm{o}\mathrm{s}{\alpha}_{3})\mathrm{c}\mathrm{o}{\mathrm{s}}^{2}{\alpha}_{1}}{EB{R}_{b}\mathrm{s}\mathrm{i}\mathrm{n}{\alpha}_{2}}+{\int}_{\alpha 1}^{\alpha 2}\frac{1.2(1+\mu )({\alpha}_{2}\alpha )\mathrm{c}\mathrm{o}\mathrm{s}\alpha \mathrm{c}\mathrm{o}{\mathrm{s}}^{2}{\alpha}_{1}}{EB[\mathrm{s}\mathrm{i}\mathrm{n}\alpha +({\alpha}_{2}\alpha \left)\mathrm{c}\mathrm{o}\mathrm{s}\alpha \right]}d\alpha ,$
where, $\mu $ is the Poisson ratio; $E$ is the elastic modulus; $B$ is the tooth width; $\alpha $ is the pressure angle; ${\alpha}_{2}=\left(\pi /2\mathrm{Z}\right)+\mathrm{t}\mathrm{a}\mathrm{n}\alpha \alpha $; ${\alpha}_{3}=\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{s}\mathrm{i}\mathrm{n}\left({R}_{b}\mathrm{s}\mathrm{i}\mathrm{n}{\alpha}_{2}/{R}_{f}\right)$; ${\alpha}_{4}=\left[\right({R}_{f1}/{R}_{b1}{)}^{2}1{]}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{1ex}{$2$}\right.}$.
Using the energy method, the Hertz contact stiffness, radial compression stiffness, bending stiffness, shear stiffness and the flexible deformation of the gear tooth are considered. The calculation relationship of the flexible deformation ${\delta}_{f}$ is shown in Eq. (2):
$\left.+{P}^{\mathrm{*}}\left({h}_{f},{\theta}_{f}\right)\left[1+{Q}^{\mathrm{*}}\left({h}_{f},{\theta}_{f}\right)\mathrm{t}\mathrm{a}{\mathrm{n}}^{2}{\alpha}_{m}\right]\right\},$
where, $B$ is tooth width, ${\alpha}_{m}$ is engagement angle, ${L}^{\mathrm{*}}$, ${M}^{\mathrm{*}}$, ${P}^{\mathrm{*}}$, ${Q}^{\mathrm{*}}$ are all represent polynomial coefficient, and satisfy the function relationship of Eq. (3):
where, ${X}^{\mathrm{*}}=\{{L}^{\mathrm{*}},{M}^{\mathrm{*}},{P}^{\mathrm{*}},{Q}^{\mathrm{*}}\}$, the flexible stiffness ${K}_{f}$ can be obtained from Eq. (4):
$+{P}^{\mathrm{*}}({h}_{f},{\theta}_{f})[1+{Q}^{\mathrm{*}}({h}_{f},{\theta}_{f}\left)\mathrm{t}\mathrm{a}{\mathrm{n}}^{2}{\alpha}_{m}\right]\}.$
When the driving and driven gears mesh in the single meshing area, according to the stiffness series principle of mechanical system, the comprehensive timevarying meshing stiffness ${K}_{sgle}$ can be calculated, as shown in Eq. (5):
When $i$ is 1, it means driving wheel. $i=$2 represents the driven wheel.
Fig. 1Force graph of single tooth engagement
3.1.2. Doubletooth meshing stiffness
Similarly, when the corresponding gear teeth of the driving and driven gears mesh in the double meshing area, they can be calculated according to the parallel stiffness principle. The calculation relationship of the comprehensive timevarying meshing stiffness ${K}_{ttal}$ of the gear teeth mesh is shown in Eq. (6):
where, when $i$ is 1, it means the first pair of meshing teeth. When $i$ is 2, it indicates the second pair of meshing teeth.
3.1.3. Timevarying meshing stiffness of single and double teeth
According to Eqs. (5) and (6), the comprehensive timevarying meshing stiffness ${K}_{\left(t\right)}$ of the gear transmission system is obtained, as shown in Eq. (7):
where, the period of gear teeth from the beginning to completely out of meshing is denoted $T$. $n$ is the sequence number of gear teeth. The time when the teeth of the second pair of driving and driven gears enter into mesh is ${t}_{1}$, ${t}_{1}=T/\epsilon $. The time when the teeth of the third pair of driving and driven gears enter the meshing stage is ${t}_{2}$, ${t}_{2}=11/\epsilon $. $\epsilon $ is the gear coincidence, and its value is within the range of (1, 2).
3.1.4. Bond graph model of timevarying meshing stiffness
Due to the nonlinear characteristics of gear tooth meshing stiffness, some signal channels that control nonlinear parameters to participate in system operation can be added to the static bond graph model, namely the switch type bonding graph element. The switch element is the characteristic mode that this system can open or close the information path on a key channel by control signals under specific conditions, so that the system has diverse states at different stages. Switch components, such as clutches in mechanical systems, electrical control switches and various control valves in hydraulic control systems, play a vital role in the precise control of system status. When the working state of these components changes in various systems, the working characteristics of this system will change from one state to another.
3.1.4.1. Modeling principle of switched power junctions (SPJ)
The method of switched power junctions (SPJ) is to add two mutually exclusive potential and flow ports on the basis of elementary 0sjunction and 1sjunction, as shown in Fig. 2. ${\mu}_{1}$ and ${\mu}_{2}$ represent two mutually exclusive Boolean parameters on key 1 and key 2, which control the working status of corresponding switches, such as on or off, as shown in Table 2. The signal flow in which the control parameters ${\mu}_{1}$ and ${\mu}_{2}$ are in the state type is also called the control key.
The control key is the key represented by the double arrow in the bond graph, which is also a power key. However, there is an essential difference between the control key and the other keys in the bond graph. The difference is that the control key controls the power variable of the corresponding key in the bond graph, while the control key itself has the function of the power variable, but the effect is small, so small that it can be ignored. In practical application, such as an ideal ammeter (voltmeter) only indicates the current value (voltage value) in the circuit, without causing the voltage value (current value) in the circuit to drop. In mechanical systems, for example, an ideal tachometer indicates only the value of the angular velocity on the rotating shaft, without consuming the power input from the spindle to the system, etc. These indicate that there is indeed a single signal connection in a local or subsystem of a system, and that the signal does not generate another power variable due to the power flow in the system (that is, the power response phenomenon). Therefore, the control key is indeed a power key in the various ports of the bond graph. However, the power consumption generated by the control key can be neglected when analyzing the system characteristics in general.
Table 2Boolean variable control mode
Component name  The way  Status type  
Put through / On  Break / Off  
Switch 1  ${\mu}_{1}$  1  0 
Switch 2  ${\mu}_{2}$ 
It can be seen from Fig. 2(a) that there are two input potential ports ${e}_{1}$ and ${e}_{2}$ and two output flow ports ${f}_{1}$ and ${f}_{2}$ on the copotential 0sjunction. If the signal flow on the 0sjunction activates a Boolean switch in the key path, the corresponding variable on the key path takes effect in the system. If key 1 is activated when ${\mu}_{1}=1$ and ${\mu}_{2}=0$, the flow quantity ${f}_{1}$ and potential variables ${e}_{1}$ on key 1 become effective. At this point, the bond path connected on the copotential 0sjunction is the bond 13, and there is a unique inflow flow quantity ${f}_{3}$ and potential variable ${e}_{1}$ at the junction. Similarly, It can be seen from Fig. 2(b) that when key 2 is activated at ${\mu}_{1}=0$ and ${\mu}_{2}=1$, flow quantity ${f}_{2}$ and potential variable ${e}_{2}$ on key 2 become effective. In addition, the bond path connected on the 1sjunction of the cocurrent junction is the bond 23, and there is a unique inflow flow quantity ${f}_{3}$ and potential variable ${e}_{2}$ at the junction. The calculation relationship between the copotential 0sjunction and the cocurrent 1sjunction is shown in Eq. (8) and Eq. (9) respectively:
Fig. 2The bond graph model of power junction
a) 0sjunction
b) 1sjunction
Copotential junction: 0sjunction:
Cocurrent junction: 1sjunction:
3.1.5. Bond graph model of timevarying stiffness
It can be seen from Eq. (7) that the timevarying stiffness $K\left(t\right)$ of gear tooth mesh is a piecewise function. According to the working principle of power junction and the conditions of state control distribution (Table 3), a bond graph model of gear mesh timevarying stiffness can be established, as shown in Fig. 3. ${\mu}_{1}$, ${\mu}_{2}$ and ${\mu}_{3}$ represent Boolean variables, which control the state of timevarying meshing stiffness system. At any time, only one boolean variable can be 1, and the other two boolean variables must be 0. Among them, the capacitive elements simulating timevarying meshing stiffness are ${C}_{1}$, ${C}_{2}$ and ${C}_{3}$ respectively. Where, ${C}_{1}={C}_{2}={C}_{3}=$1 / ${K}_{\left(t\right)}$. And the calculation relationship between each key path is shown in Eq. (10):
Table 3Control state table of timevarying meshing stiffness
Element name  Representation way  Status type  
$nT\le t<nT+{t}_{1}$  $nT+{t}_{1}\le t\le nT+{t}_{2}$  $nT+{t}_{2}\le t<(n+1)T$  
Switch 1  $\mu 1$  1  0  0 
Switch 2  $\mu 2$  0  1  0 
Switch 3  $\mu 3$  0  0  1 
Fig. 3The bond graph model of timevarying stiffness
3.2. Teeth surface sliding friction
3.2.1. Friction mechanism
When studying the friction mechanism of tooth surface, the influence of tooth side clearance and oil film on friction performance is not considered, and the tooth surface friction is regarded as Coulomb friction. In the process of gear transmission, the meshing teeth produce a sliding form of dry friction on the tooth surface. This form of friction has nothing to do with the contact area, but is proportional to the normal positive pressure at the contact surface, and the direction of friction is opposite to the relative sliding direction. When the meshing point of each pair of meshing gear teeth is at the upper and lower sections of the pitch circle, the direction of the sliding velocity of the tooth surface will change at the meshing point of the pitch circle, and the direction of the sliding friction force will also change. In addition, due to the fact that the geometric structure of the gear itself is not very ideal and is affected by the surrounding transmission environment, the position of each pair of gear teeth is different, and the meshing force will also change with the change of the meshing position, which results in the normal positive pressure at the meshing point. The size and direction of the tooth surface friction value will change periodically. According to the Coulomb sliding friction principle, the calculation relationship of basic friction is shown in Eq. (11):
where, ${F\text{'}}_{f}$ is the Coulomb friction force on the tooth surface, $\mu $ is the viscosity friction coefficient, ${F}_{N}$ is the positive pressure, $sign\left({V}_{\lambda}\right)$ is a symbolic function, ${V}_{\lambda}$ is the direction control variable of the friction force, as shown in Eq. (12):
where, ${V}_{\lambda}$ is the relative sliding speed of the engagement point, as shown in Eq. (13):
where, $\alpha $ is the gear pressure angle, ${r}_{a1}$ is the addendum circle radius of the driving gear, ${r}_{b1}$ is the addendum circle radius of the driving gear, ${r}_{a2}$ is the addendum circle radius of the driven gear, ${r}_{b2}$ is the addendum circle radius of the driven gear, ${\omega}_{1}$ is the rotational speed of the driving gear, and ${\omega}_{2}$ is the rotational speed of the driven gear.
3.2.2. Friction bond graph model of tooth surface
Friction is an energyconsuming system in gear transmission system, which can be simplified as an energyconsuming element $R$ to describe. The friction effect of tooth surface is simulated with 0sjunction and SPJ switch; The relationship between the direction of the friction force on the tooth surface and the relative meshing speed is simulated by using 1sjunction and SPJ switch. The friction coefficient $\mu $ is simulated by using the adjustable resistance conversion element (MR). ${F}_{N}$ represents control key of parameter signals. In addition, ${\mu}_{1}$ and ${\mu}_{2}$ are two mutually exclusive Boolean variables. Only one variable can be activated in the same time period to control the direction of friction, as shown in Table 4.
Table 4Boolean variable assignment table
$sign\left({V}_{\lambda}\right)$ function  Boolean variable  
${\mu}_{1}$  ${\mu}_{2}$  
1  1  0 
–1  0  1 
Combined with the modeling principle of switch type bond graph, the bond graph model of gear tooth surface friction can be obtained, as shown in Fig. 4.
Namely, the characteristic equation of the gear tooth surface friction system is shown in Eq. (14):
3.3. Engagement damping
The gear tooth meshing damping $\psi $ is approximately calculated by empirical formula, as shown in Eq. (15):
where, $\xi $ is damping ratio, ${K}_{m}$ is the average meshing stiffness, ${J}_{1}$ is the rotational inertia of the driving wheel, ${J}_{2}$ is the rotational inertia of driven wheel, ${r}_{b1}$ and ${r}_{b2}$ respectively represent the dividing circle radius of the driving and driven wheels.
Meshing damping is an energy dissipation system, which is simulated by resistive elements $R$. The bond graph model of gear tooth meshing damping is shown in Fig. 5.
Fig. 4Bond graph model of tooth surface friction
Fig. 5Bond graph model of meshing damping
3.4. Steadystate transmission error
3.4.1. Source of steadystate transmission error
There are many kinds of gear transmission error, mainly including manufacturing errors, installation errors and improper operating conditions, which often lead to different forms of damage to gears. These damage forms will cause vibration or noise in the system during gear operation. It not only worsens the system’s operating environment and reduces the transmission efficiency, but also decreases the service life of transmission components. Since steadystate transmission error generates a periodic excitation source for the gear system, it is described by a harmonic function ${e}_{\uff08t\uff09}$, shown in Eq. (16):
where, ${e}_{0}$ is the mean steadystate transmission error, ${e}_{r}$ is the amplitude of steadystate transmission error fluctuation, ${\omega}_{{\mathrm{}}_{h}}$ is the meshing frequency, ${\omega}_{{\mathrm{}}_{h}}={n}_{1}{z}_{1}/30$ ,${n}_{1}$ and ${z}_{1}$ are the rotational speed and teeth number of the driving wheel, respectively, ${\varphi}_{e}$ is the initial phase angle.
The derivative form of steadystate transmission error, whose derivative form is shown in Eq. (17):
where, ${e}_{r}$ is the fluctuation amplitude, ${\omega}_{{\mathrm{}}_{h}}$ is the meshing frequency, ${\varphi}_{e}$ is the initial phase angle.
3.4.2. Bond graph model of steadystate transmission error
The steadystate transmission error of the gear teeth is a displacement excitation, and its derivative form is the rotational speed characteristic of the system. Flow source $Sf$ is adopted for simulation, and its bond graph model is shown in Fig. 6.
Fig. 6The bond graph model of steadystate transmission error
4. Establishment of dynamic model of wheel system nonlinear system
The nonlinear dynamic bond graph model of PX singleloop gear train is shown in Fig. 7.
5. Study on dynamic characteristics of nonlinear single loop gear system
5.1. Dynamic characteristics of the singleloop ring system of PX type
Fig. 8(a) shows the phase trajectory diagram of PX singleloop gear system. The trajectory line circles point (0,0) clockwise, but does not surround point (–1, 0$j$), and passes through the negative real axis from the right side of the point. Therefore, the analysis shows that the stable state of PX single loop gear system is unstable. As shown in Fig. 8(b), the amplitudephase characteristic diagram of PX singleloop gear system shows that the amplitude margin of this system ${K}_{g}=$ 11.6 dB > 0 (0.136 rad/s).
5.2. Xtype subsystem
Fig. 9 shows the amplitudephase characteristic diagram. The subsystem is cutoff from the input end to type X, and the stable state of Xtype subsystem inside the singleloop system is studied. The analysis demonstrates that the dynamic characteristics of this subsystem in PX single loop gear system are unstable.
Fig. 7Nonlinear bond graph model for singleloop gear system of PX type
5.3. Ptype subsystem
Similarly, Fig. 10(a) shows the amplitudephase characteristics of the subsystem, which is cut off from the input to the Ptype. From Fig. 10(b), it can be concluded that the amplitude margin ${K}_{g}=$ 59.6 dB > 0 (2.13 rad/s).
Fig. 8Characteristic phase diagram of PXtype single loop gear system
a) Phase diagram
b) The amplitudephasefrequency characteristic
Fig. 9Characteristic phase diagram of the Xtype subsystem
a) Phase diagram
b) The amplitudephasefrequency characteristic
Fig. 10Characteristic phase diagram of the Ptype subsystem
a) Phase diagram
b) The amplitudephasefrequency characteristic
6. Conclusions
In this research, study on nonlinear dynamic characteristic of singleloop gear system based on bond graph method. Simulation has been conducted to prove the efficiency of this method. The conclusions are listed as follows:
1) Based on the mechanism study of nonlinear factors in the gear meshing system, such as timevarying meshing stiffness, tooth surface friction, gear tooth damping and steadystate transmission error, the bond graph model of these nonlinear factors is obtained.
2) According to the structural characteristics of the single loop gear system, the dynamic characteristics and modeling methods are studied. The nonlinear mathematical models of the PX are established by using the bond graph method. Then, the amplitudephasefrequency characteristics are obtained respectively to study the stability of the single loop gear system.
3) According to the dynamic characteristics of PX single loop gear system, it is concluded that, when the whole single loop gear system is unstable, its system module also has the same unstability.
References

E. L. Esmail, E. Pennestrì, and M. Cirelli, “Powerflow and mechanical efficiency computation in twodegreesoffreedom planetary gear units: new compact formulas,” Applied Sciences, Vol. 11, No. 13, p. 5991, Jun. 2021, https://doi.org/10.3390/app11135991

W. Wang et al., “A multiobjective power flow optimization control strategy for a power split plugin hybrid electric vehicle using game theory,” Science China Technological Sciences, Vol. 64, No. 12, pp. 2718–2728, Dec. 2021, https://doi.org/10.1007/s1143102017703

X. Zhou, Y. Cui, L. Li, L. Wang, X. Liu, and B. Zhang, “Signal denoising in gear pitting fault identification by an improved singular value decomposition method,” Forschung im Ingenieurwesen, Vol. 84, No. 2, pp. 79–90, Jun. 2020, https://doi.org/10.1007/s10010020004007

Z. Jahanbin, A. Selk Ghafari, A. Ebrahimi, and A. Meghdari, “Multibody simulation of a flappingwing robot using an efficient dynamical model,” Journal of the Brazilian Society of Mechanical Sciences and Engineering, Vol. 38, No. 1, pp. 133–149, Jan. 2016, https://doi.org/10.1007/s4043001503504

Y. S. Hamed, K. M. Albogamy, and M. Sayed, “Nonlinear vibrations control of a contactmode AFM model via a timedelayed positive position feedback,” Alexandria Engineering Journal, Vol. 60, No. 1, pp. 963–977, Feb. 2021, https://doi.org/10.1016/j.aej.2020.10.024

Jedliński, A. Syta, J. Gajewski, and J. Jonak, “Nonlinear analysis of cylindrical gear dynamics under varying tooth breakage,” Measurement, Vol. 190, p. 110721, Feb. 2022, https://doi.org/10.1016/j.measurement.2022.110721

S. Liu, A. Hu, Y. Zhang, and L. Xiang, “Nonlinear dynamics analysis of a multistage planetary gear transmission system,” International Journal of Bifurcation and Chaos, Vol. 32, No. 7, Jun. 2022, https://doi.org/10.1142/s0218127422500961

Ouyang T. et al., “Effect of structural design parameters on nonlinear dynamic characteristics of the gear transmission,” Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, Vol. 236, No. 4, pp. 522–539, 2022.

H. Dong, H. Q. Zhang, X. L. Zhao, and L. L. Duan, “Study on dynamic loadsharing characteristics of face gear dualpower split transmission system with backlash, support and spline clearance,” Mechanical Sciences, Vol. 12, No. 1, pp. 573–587, May 2021, https://doi.org/10.5194/ms125732021

W. Zhou, R. Zhu, W. Liu, and Y. Shang, “An improved dynamic transmission error model applied on coupling analysis of gear dynamics and elastohydrodynamic lubrication,” Journal of Tribology, Vol. 144, No. 5, pp. 51601–51610, May 2022, https://doi.org/10.1115/1.4051813

XX.F. Gou, G.Y. Li, and L.Y. Zhu, “Dynamic characteristics of a straight bevel gear drive system considering multistate meshing and timevarying parameters,” Mechanism and Machine Theory, Vol. 171, p. 104779, May 2022, https://doi.org/10.1016/j.mechmachtheory.2022.104779

H. Li, S. Chen, J. Tang, Z. Sun, and Y. Hu, “Nonlinear dynamic modeling and analysis of spur gear based on gear compatibility conditions,” Mechanism and Machine Theory, Vol. 171, p. 104767, May 2022, https://doi.org/10.1016/j.mechmachtheory.2022.104767

Hou X. et al., “A modified damping model of vector form intrinsic finite element method for highspeed spiral bevel gear dynamic characteristics analysis,” Journal of Strain Analysis for Engineering Design, Vol. 57, No. 2, pp. 144–154, 2022.

Y. Nie et al., “Phenomenological vibration models of planetary gearboxes for gear local fault diagnosis,” Mechanism and Machine Theory, Vol. 170, p. 104698, Apr. 2022, https://doi.org/10.1016/j.mechmachtheory.2021.104698

H. T. Chen, X. L. Wu, D. T. Qin, J. Yang, and Z. G. Zhou, “Effects of gear manufacturing error on the dynamic characteristics of planetary gear transmission system of wind turbine,” Applied Mechanics and Materials, Vol. 86, pp. 518–522, Aug. 2011, https://doi.org/10.4028/www.scientific.net/amm.86.518

C. Dong and Y. Liu, “Experimental study on transmission error and dynamic backlash of elliptic gear transmission system,” Proceedings of the Institution of Mechanical Engineers, Vol. 236, No. 1, pp. 130–139, 2022.

Z. Cao et al., “Dynamic simulation and experimental study of electric vehicle motorgear system based on state space method,” Machines, Vol. 10, No. 7, p. 589, Jul. 2022, https://doi.org/10.3390/machines10070589

X. S. Wang, S. J. Wu, J. C. Hu, and J. Chen, “Chaos and bifurcation analysis of gear pair with wear fault,” in International Conference on Material Sciences and Manufacturing Technology, pp. 516–520, 2012.

J. RodríguezGuillén, R. SalasCabrera, and P. M. GarcíaVite, “Bond Graph as a formal methodology for obtaining a wind turbine drive train model in the perunit system,” International Journal of Electrical Power and Energy Systems, Vol. 124, p. 106382, Jan. 2021, https://doi.org/10.1016/j.ijepes.2020.106382

M. G. Tsegay, P. M. Pathak, A. K. Samantaray, and R. Merzouki, “Bond graph modeling of a spatial multisection soft bionic robot,” Mechanism and Machine Theory, Vol. 174, p. 104902, Aug. 2022, https://doi.org/10.1016/j.mechmachtheory.2022.104902
About this article
This work was supported by grants from the Introduce Highlevel Talents to Start Scientific Research Funds by Shaanxi Polytechnic Institute of China (No. BSJ 202111) and (No. 20209), which were highly appreciated by the authors.
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
Wenting Zhang: writing papers, mechanical modeling, simulation analysis, etc. Jiamin Liu: studied nonlinear dynamics of gear clearance system and worked on thesis translation, etc. Xintao Zhou: writing papers, providing research methods, investigation, validation, and experimental verification, etc. Yahui Cui: formal analysis, and investigation, etc.
The authors declare that they have no conflict of interest.