Research Article  Open Access
Zizhen Qiao, Jianxing Zhou, Wenlei Sun, Xiangfeng Zhang, "A Coupling Dynamics Analysis Method for TwoStage Spur Gear under Multisource TimeVarying Excitation", Shock and Vibration, vol. 2019, Article ID 7350701, 18 pages, 2019. https://doi.org/10.1155/2019/7350701
A Coupling Dynamics Analysis Method for TwoStage Spur Gear under Multisource TimeVarying Excitation
Abstract
A new modeling method is proposed to simulate the dynamic response of a twostage gear transmission system using the finite element method (FEM). The continuous system is divided into four modules: shaftshaft element, shaftgear element, shaftbearing element, and geargear element. According to the FEM, the model is built with each element assembled. Meanwhile, the model considers the timevarying mesh stiffness (TVMS), bearing timevarying stiffness (BTVS), and the shaft flexibility. The Newmark integration method (NIM) is used to obtain the dynamic response of the spur gear system. Results show that vibration amplitude and the number of frequency components decrease after considering shaft flexibility through comparing the gear dynamic response under the condition of flexible shaft and rigid shaft. When the effect of bearing stiffness is considered, there will be a bearing passing frequency component in the frequency spectrum. In addition, the result shows that the simulation and experimental test of the frequency component are basically consistent. Furthermore, the theoretical model is validated against an experimental platform of the twostage gear transmission system and the dynamic responses are compared under the condition of increasing speed. Additionally, the increase of shaft stiffness not only changes some of the dominant mode shapes (torsional mode shapes) but also makes the number of primary resonance speeds added. The method can be used to guide the design of gear systems.
1. Introduction
The gear transmission system has many advantages, such as high efficiency, compact structure, and stable transmission ratio, so it is one of the most important mechanical transmission methods [1, 2]. Hence, it is extensively used in highspeed, largepower, and lowspeed, hightorque applications, such as automotive transmissions, gas turbines, marine main reducers, tanks, engineering machinery, and 3D printers [3]. Sometimes, the complexity of the structure of the gear system, the timevarying stiffness, and other nonlinear factors can become an excitation source leading to premature failure of the transmission system [4]. For this reason, the coupled dynamic response of the gear system has recently become an important research topic. Currently, the lumped parameter method and finite element methods (FEM) are two common modeling approaches applied to the analysis of coupling vibration in the gear system [5]. Presently, the lumped parameter method applied to model multistage gear dynamics assumes that the gear, shaft, and bearing are concentrated points. Additionally, the complex coupling effect between the elastic shaft elements is represented with simple bending, torsional, and axial stiffness. These assumptions often render the lumped parameter models less accurate because they are not able to reliably distinguish between the shaft section and bearing and are not able to account for the geometrical dimensions of the shaft [6–10]. In contrast, the FEM has good calculation accuracy and is of greater importance for predicting the vibration noise of the system. In the literature, use of the FEM for calculating dynamic response has been reported by many researchers. Neriya and Sankar [11] introduced the FEM in the study of the dynamics of gear rotor systems and established the linear modeling of meshing stiffness. Modak et al. [12] updated a finite element model of a structure so that an updated model predicts more accurately the dynamics of a structure. With the development of the research, the coupling between parts has been considered, and the coupling model of the gear system has been established. ÖzgÜven and Özkan [13] established the dynamic model of a shaftgearbearing system to study the response of a coupled system. Ouyang et al. [14] investigated dynamic characteristics of gearrollerbearing system by use of experimental test and theoretical modeling, and dynamic characteristics of gearrollerbearing system are investigated using of experimental test and theoretical modeling. Kubur et al. [15] established the dynamic model of a gear system with N shafts and N − 1 meshing pairs to study the relationship between the dynamic load of the bearing and the shaft length. Jiang et al. [16] established a complete mechanical model of rotorgear system using both a lumped parameter model and a finite element model (shaft, gear, and rotor elements are modeled using the lumped parameter method and assembled with the FEM). Guo et al. [17] established a gearbearinghousing system and investigated the vibration propagation of gear dynamics using mathematical modeling and finite element analysis. Ma et al. [18] studied the effects of varying amounts of profile relief on the mesh stiffness of a spur gear with FEM and substituted it into the dynamic model of the gear system. Jing et al. [4, 9] proposed a coupling dynamic method for a multistage planetary gear system and analyzed coupling vibration mechanism in the planetary gear system. Qian et al. [19] used finite element modeling to establish dynamic modeling of a singlestage gear transmission system considering the shaft, gear gyroscopic effects, and they calculated the natural characteristics of the system. Ankur et al. [20] established the shaft, bearing, and meshing gear pairs separately, and the viscous damping model is used to represent bearing and gear mesh damping. Wang [21] established a twodimensional (2D) spur gear with FEM to calculate the static transmission error, single and combined torsional mesh stiffness, tooth loadsharing ratio, maximum tooth root stress, and surface contact stress against various input torques over a complete mesh cycle.
The bulk of existing literatures established the coupling dynamic modeling for a gear system, and the dynamic characteristics were studied to some extent. Few studies have established the twostage spur gear system with FEM and conducted on the comprehensive influence of TVMS, BTVS and different shaft stiffness on the mode shape and the vibration transmission of the system simultaneously. This study considers the excitation of TVMS and BTVS as well as error excitation and also includes the flexibility of transmission shaft. Then, the FE modeling of the twostage spur gear transmission system is established with FEM. The influence of different shaft stiffness on mode shape and dynamic response are analyzed. The amplitude and frequencyresponse curves of the vibration acceleration of the bearing are obtained with simulation and compared with the experimental results.
2. Dynamic Model of Geared System
The general gear system includes input/output system, transmission system consisting of gear shaft system, and structural system made up of bearing seats. The twostage spur gear transmission system is shown in Figure 1.
The parameters of the system are shown in Tables 1 and 2.


2.1. Coupling Model of Gear System
The gear transmission system is a continuous elastomer with countless degrees of freedom (DOFs). In order to accurately model the dynamic performance and efficiently understand the coupling vibration mechanism in the twostage gear system, the FEM is put forward in this study. In the proposed approach, the system is divided into many elements.
Then, the FEM dynamic model of the twostage gear transmission system is established, including the shaft, gear, and bearing elements, as shown in Figure 2. In this, the first shaft is divided into 12 elements, the second one 8 elements, and the third one 9 elements. Each element node is numbered. There are two rolling ball bearings on each shaft forming shaftbearing coupling element. Furthermore, on the No. 7 node of the first shaft and the No. 16 node of the second shaft are the firststage pinion and gear, respectively; on the No. 20 node of the second shaft and the No. 29 node of the third shaft are the secondstage pinion and gear, respectively. In this study, different nodes can realize physical connection and mechanical coupling of shaftbearing, shaftgear, and geargear through coupling. In the specific application, we can properly add the number of nodes at different shafting sections and divide more shafting elements according to the assembly requirements of shafts, the geometrical shape, practical accuracy requirements, and so on.
2.2. Shaft Element
The shaft is mainly for supporting rotation components and transmitting power, thus it is one of important parts in the gear transmission system [22]. The element models of elastic shafts are built in order to effectively consider shaft’s deformation, as shown in Figure 3. Space beam element theory is used for force analysis, which is mainly Euler–Bernoulli beam theory or Timoshenko beam theory [23]. Based on Table 1, the width diameter of the shaft in this gear system modeling is quite small, and the shearing deformation should be considered; thus, Timoshenko beam theory can be used to analyze the shafting section elements.
As shown in Figure 3, every simple shaft can be divided into many elastic shafting elements. The length of the element is l, the elastic modulus is E, and the rotational inertia is J. The two nodes of the shaft element are i and j. Since this model is a parallel shaft spur gear transmission system, there is almost no excitation in the axial direction. Therefore, each node has been given three DOFs, and each node of the shaft has two translational DOFs and 1 torsional ones. The generalized coordinates of the shaft element arewhere the subscripts indicate the shafting section, and are the translational DOFs, and is the torsional DOF. The stiffness matrix and the mass matrix of shafting element 6 × 6 are obtained as follows:
Damping can also be added to the shaft model via Rayleigh damping aswhere denote mass matrix and stiffness matrix of shafting element and denote the scale factors. The specific calculating method can be referred to [24].
2.3. Mesh Element and Mesh Stiffness
The meshing action of the gears is equivalent to that of a spring and damper. For DOF analysis, besides the torsional vibration, other kinds of vibration should be taken into account, because there is almost no excitation in the axial direction for spur gear. It can be thought that there are transverse and torsional vibrations without axial DOF of the meshing elements. As shown in Figure 4, the dynamics model of gear meshing pair is established, where and are pinion and gear, is the TVMS of meshing elements, denote the transverse vibration displacement and torsional angle of pinion and gear respectively, and α denotes the pressure angle.
The generalized coordinates of the meshing element can be expressed as
The vibration of the gear along the torsional DOF and translational DOF can change the meshing of gear pairs. The projection vector of the displacement of freedom of the gear in the meshing line direction [25] can be expressed in the following equation.where r_{p} and r_{g} denote the radius of pinion and gear.
The elastic meshing force of gear pair can be expressed in accordance with the following equation [25]:where is the comprehensive error and , in this, e_{1} denotes the error amplitude and is the firststage meshing frequency. The stiffness matrix of meshing pair [25] can be obtained as follows:
The finite element method was adopted to calculate the mesh stiffness of several meshing positions within a meshing period. In application of boundary condition, the gear is subjected to full constraint, the radial and axial constraint is applied to the pinion, and the unit force is applied to the circumferential node to solve, as shown in Figure 5(a). The gear pair appears deformed to a certain extent under the action of external load, and the pinion appears a small rotation angle. The microangle θ of the pinion can be expressed aswhere is the mean value of gear tooth deformation and is the radius of the inner ring gear of the pinion.
The small deformation of the inner ring of the pinion can be converted into the deformation of the contact gear tooth with the relationship between the small rotation angle and the radius of the base circle. The mesh deformation δ_{s} can be defined aswhere is the radius of the base circle of the pinion.
Through the relationship between the unit force applied by the inner ring node of the pinion and the radius of the inner ring and the radius of the base circle of the pinion, the node load applied to the inner ring of the pinion is converted into the mesh force of the gear teeth. The meshing force can be expressed from the following equation:where n is the number of inner ring nodes.
Finally, the stiffness can be expressed as follows [26]:
The mesh stiffness curve of the gear teeth from mesh to mesh (one mesh period) is obtained with FEM as shown in Figure 5.
2.4. Bearing Element and Bearing Stiffness
Bearing is mainly to support the shaft system and connect to shaft system and box. The stiffness of rolling bearing is timevarying under complicated working conditions of high rotating speed and varying loads. Since the model is a rolling ball bearing and there is almost no excitation in the axial direction, the model of deep groove ball bearing model and two bearing loading forms are shown in Figure 6 [27].
The stiffness matrix of bearing can be obtained as follows:
Usually, the form of motion of the bearing has two situations. One is shown in Figure 6(b)—radial load falling on the lowest roller bearing (odd pressure); the other is shown in Figure 6(c)—radial load falling between the two lowest roller bearings (even pressure). The curvature sum and subtraction of bearings are obtained from the previous work [25]. Table 3 gives the basic parameters of the bearing and the calculated numerical number.

For inner race contact, the contact displacement is δ_{1} = 0.59, and for outer race contact, it is δ_{0} = 0.66 [24]. The loaddisplacement coefficient during point contact can be expressed as
Then, we can obtain the loaddisplacement coefficients of inner and outer race contacts, respectively.
The loaddisplacement coefficient of normal load direction is as follows:
Radius displacement of bearing can be calculated in the following equation:where Q_{max} denotes the maximum normal load, D denotes the diameter of bearing race, and denotes the contact angle of deep groove ball bearing.
Bearing stiffness for odd pressure can be expressed aswhere F_{r} is the radial load.
Under even pressure, there is no direct force on the rollers with the influence of the radial load. Take the load of ψ = ±22.5° as the maximum load of rolling element to get the bearing stiffness under even pressure. Thereby, the stiffness under even pressure and odd pressure can be obtained. They can be expressed as follows: and . Hence, the timevarying stiffness is expressed aswhere denote the static bearing stiffness and the waving amplitude of timevarying stiffness respectively and denotes the bearing passing frequency and the initial phase angle of bearing. Meanwhile, .
3. System Model
3.1. TwoStage Phase Relationship
In order to establish twostage gear pair modeling and determine the twostage phase relationship, firstly a coordinate system is established as shown in Figure 7. It is assume that is theoretical central line of the firststage gears and is vertical to the central line of those two gears.
In the twostage parallel shaft spur gear system, there is a translational vibration with no vibration in the axis, as shown in Figure 7, where are the translational vibration displacement and torsional angle of the firststage pinion and gear, are the translational vibration displacement and torsional angle of the second stage pinion and gear, is the gear pressure angle, β is the shaft intersection angle among three shafts. According to the geometrical relationship, the intersection angle between the second stage gear meshing line and is reformulated as follows:
Therefore, in the second stage gear system, the projection displacement vector in the direction of meshing line can be expressed in accordance with the following equation [28]:
Meshing force of gear pair can be written as follows [28]:where are the radii of base circle of pinion and gear in the secondstage transmission system, respectively, is the TVMS of second stage gear, is the comprehensive error, further explanation, , where denotes the error amplitude, and in this, is the secondstage gear meshing frequency.
3.2. Overall System Model
According to the physical quantities of stiffness matrix, mass matrix, load vector, and so on, the elements are assembled with the FEM. More specifically, the stiffness matrix of each element is sent into the overall stiffness matrix according to the assembly rules of the overall system. The generalized displacement vector of the system is as follows:
where n is the number of nodes, n = 1, 2, …, 32. The twostage gear transmission system with three shafts, six bearings, and two gear pairs is divided into 29 elements and 95 DOFs with the FEM under the boundary conditions (the No. 1 node connects with input end and it is subjected to a torsional DOF constraint, and the last node bears load torque).
The overall stiffness matrix assembly of the twostage gear transmission system is shown in Figure 8, and the dynamic equations of the spur gear system can be expressed in accordance with equation (21).where denotes the node displacement vector and are the total mass, Rayleigh damping, stiffness matrix, and force vectors, respectively.
Dynamic calculation is done at last. Because stiffness matrix and displacement vector are periodic functions of time, they need to be calculated with the numerical integration. NIM has the format of unconditional stability with the adjusted parameters constant. It is verified that the calculation error and calculation are both less than the Cotes integration method, and the calculation efficiency is better than the traditional numerical method through tests [27]. Therefore, NIM can be used to analyze the dynamic response.
4. Static Characteristic Analysis
4.1. Transmission Shaft Deformation Analysis
Under loads, all the transmission shafts will bear the bending moment and torque produced by the gear’s peripheral force and radial force, so deformation is inevitable. If the inertial force and damping force of the shafting element are not considered, the force of gear shaft can be analyzed based on statics. The deformation of shaft can be defined through load and transmission shaft’s overall stiffness matrix. Therefore, the deformation can be solved with the following equation:where k_{s} is the overall stiffness matrix of the shaft and F is the force of shaft.
The deformation of gear shaft is shown in Figure 9 (compared with the actual results, it is amplified by 10^{5} times), under a constant input speed of 500 rpm and a constant load of 100 N·m with simulation. The shaft couples bearing support and meshing gear pair. It is supposed that the maximum deformation of each shaft is and the first shaft has the maximum force at the No. 7 coupling node of the gear, with the maximum deformation of 0.55 μm; the second stage gear has small rotated speed and large load, so compared with the firststage transmission shaft, it has larger deformation at the gear coupling position. The second transmission shaft has the maximum deformation of 1.98 μm at the No. 20 node; the third transmission shaft has the maximum deformation of 1.70 μm at the No. 29 node.
4.2. Static Transmission Error Analysis
The concept of transmission error was put forward by Harris [29] in 1958. It is defined as the active rotation over a certain angle of the pinion, the difference between the actual rotation angle and the theoretical rotation angle of the gear. If the gear has small rotated speed, the transmission error can be calculated based on static questions [30, 31]. In this study, the static transmission error (STE) is simulated. The STEs of the gears are, respectively, expressed aswhere (n = 18, 19, 20, 45, 46, 47, 57, 58, 59, 85, 86) is the displacement vector at the nth DOF.
The calculation result is shown in Figure 10. It can be found that in the process of single and double tooth alternation, because of the gear change of meshing stiffness, the STE will have great fluctuating changes. Compared with other positions, the STE has stable changes.
5. Modal Analysis
Without outer excitation of the system, the equation of free vibration can be obtained from formula (21) as follows:
During modal analysis, the damping exerts little influence on the natural frequency and mode shape of the gear system, so the damping can be neglected. Without considering the damping effect, the vibration equation of the system can be expressed as
The eigenvalues of the system can be expressed as follows:where is the average value of system stiffness and denote the natural frequency and mode shape of the system at i^{th} order, respectively. The characteristic equation can be written as
Finally, the natural frequency and mode shapes of the system with the influence of the flexible shaft can be obtained. The system has natural frequency of 95 orders in total, and the first 20 orders natural frequencies are shown in Table 4.

The continuous system is divided into shaftshaft, shaftgear, and other basic elements. The influence of geometric dimension and transmission shaft flexibility is considered to the model as described in Section 2.2. This section studies the effect of FE modeling with shaft flexibility on the coupling mode shapes of twostage gear transmission system consisting of the torsion modes of gear pair, translational vibration, and bending modes of the transmission shaft through simulation. Owing to the effect of transmission shaft’s flexibility, the transmission shafts show highorder bending mode shapes. The typical mode shapes are shown in Figure 11. The dominant vibration of the first order is torsional vibration of the first gear pair. The dominant vibration of third order is torsional vibration of the second stage gear. It has different vibration directions from torsional vibration of the first order, as shown in Figures 11(a) and 11(b) which show that the dominant vibration of the second order and fourthorder are translational vibrations of the gear; the dominant mode of sixth order and ninth order is loworder bending mode shape of the shaft, as shown in Figure 11(c). For the sixthorder mode shape, the shaft 1 and shaft 3 have concave bending vibration, and shaft 2 has convex ones; the three shafts at ninth order have convex vibration shape; because of the flexibility of transmission shaft, there are highorder bending mode shapes, as shown in Figure 11(d).
(a)
(b)
(c)
(d)
Consequently, the traditional lumpedmass model simplifies the transmission shaft into elastic elements of torsional and bending deformation. In modal analysis, the torsional mode shape and translational mode shape of the gear pair can only be analyzed, but the complex mode changes of the transmission shaft cannot be shown. Therefore, the FEM is used to optimize the modeling of twostage gear transmission and modify the gear structure to reanalyze its natural frequencies and the mode shapes. So, the method can be used for further studies.
6. Dynamic Characteristic Simulation Analysis
6.1. The Effect of Shaft Flexibility
Based on the twostage gear system proposed in Section 2, the system dynamic responses can be calculated with Newmark integration method. The curves of the dynamic meshing force and frequency spectrum are shown in Figures 12 and 13 under a constant input speed of 500 rpm and a constant load of 100 N·m. It can be clearly found that the dynamic response of the gear has periodic stable signals with no amplitude modulation. Because of small rotating speed but large load of the second stage gear, its dynamic load vibration amplitude is rather larger than that of the firststage gear. The major composition of the frequency spectrum is meshing frequency and its frequency multiplication (, ). Meanwhile, there is a natural frequency () at lowfrequency stage.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
Although the coefficient of the Fourier series reduced with the increase of harmonics order, because of the influence of frequencyresponse function and nonlinear feedback which makes the vibration energy of the system not strictly obey the change law [32]. In the firststage gear system, the maximum amplitude frequency component is f_{m1}. In the secondstage gear system, the peak amplitude is 2f_{m2}. This is because the frequency component of 2f_{m2} is similar to the natural frequency of the second order and intensifies the vibration.
Next, the effect of the shaft flexibility on the system dynamic response is described. If the flexibility of shaft (elastic modulus: E = 2.1 × 10^{11} Pa) is considered, the maximum dynamic load decreases, the shock vibration of the dynamic load of the twostage gears will greatly reduce, and the response signals are more stable compared with the timedomain process. The analysis of the frequency domain shows that, if the transmission shaft is considered a stiff body (elastic modulus: E = 2.1 × 10^{12} Pa), it shows a greater frequency signal, such as the natural amplitudes of the first order (), and as well as its harmonics (2f_{m1}, 3f_{m1}, etc) and 2f_{m2}, 5f_{m2}, 7f_{m2}, and so on. In this, the peak amplitude of f_{m1} is the largest. In summary, the dynamic load amplitude decreases, and the vibration response frequency composition increases under the influence of the transmission shaft flexibility. This is because the shaft’s flexibility can absorb energy and decreases the transmission of vibration of the gear system. The result is consistent with the analysis in the previous work [15].
6.2. The Effect of Bearing TimeVarying Stiffness
The curves in Figures 14 and 15 compare the radial force response of the shaft 2bearing coupled node of the system with constant bearing stiffness and TVMS in the time domain and frequency domain. It is found that considering BTVS, there is a bearing passing frequency () component at the low frequency. On the other hand, the introduction of BTVS has little effect on the gear dynamic response. This is because gear meshing frequency and its harmonics play a dominant role in response. f_{b2} has a very small influence compared with that. The static stiffness of the bearing has a significant influence on the dynamic response of the gear system, which is consistent with the conclusions of Cao et al. [33].
(a)
(b)
(a)
(b)
7. Experimental Study on Dynamic Characteristics
The experiment uses the test bench of the twostage gear transmission system made by SQI to analyze the dynamic characteristics and collect the timedomain signals of radial acceleration at the endcap position of the bearing, as shown in Figure 16(a). Meanwhile, the twostage gear transmission system has consistent parameters with the built model. The acceleration sensor 608A11 is placed in the gearbox. The measuring points are shown in Figure 16(b). Then, collection card DT9837 is used to output and analyze the vibration signals measured.
(a)
(b)
The experiment was done to collect the timedomain signals of radial acceleration at the endcap position of the bearing under an input speed of 960 rpm and constant load of 10 N·m with sampling frequency of 10 kHz. Then, the frequencydomain was obtained using Fourier transform, as shown in Figure 17. The figure shows that the major frequency composition includes f_{m1} and f_{m2} and their harmonics. In addition, at lowfrequency stage, there was a passing frequency composition (f_{b}) of the bearing.
(a)
(b)
The vibration response of radial acceleration at the endcap position of the bearing is simulated as shown in Figure 18, where f_{m1} is 576 Hz, f_{m2} is 186 Hz (, and ).
(a)
(b)
It is clearly seen from Figure 18 that only at a few moments there is little shock and vibration. The vibration response signals are relatively stable. The node displacement acceleration has periodic changes. Moreover, the major frequency composition includes the passing frequency of bearing () at the low frequency because of the BTVS of bearing. In addition, the peak amplitudes of and its harmonics and the amplitude of 2f_{m1} (1152 Hz) are large. This is because 2f_{m1} is basically the same as the fifthorder natural frequency of the system, which increases the vibration energy of the system. The and its doubling frequency also show vibration peaks. The harmonic (1299 Hz) is basically similar to the natural frequency of the sixth order of the system, so the vibration is the strongest here. There are obvious sidebands modulated by f_{b} and its harmonics. The loworder natural frequency of the system is listed in Table 4.
Furthermore, was also an excitation component in the experimental results, as shown in Figure 17(b). There is no such response for the simulation result because there are no basic dynamic characteristics of the testbench. In contrast, the gear system is at the test bench in the test. In the meantime, the natural frequency of the test bench is found to be through beat test.
When the torque is constant, the accelerations of the measuring point in the radial direction vary with the growth of speed, as shown in Figure 19. It is obvious that the two curves both have two resonance peaks of A and B. Peak A is for an input speed of 1700 rpm, where the 2f_{m2} is near the natural frequency of the third order (1055 Hz), so that the system has a strong vibration. Peak B is for an input speed of 2200 rpm, where f_{m1} and 2f_{m1} are near the natural frequency of orders six and eight (1230/2875 Hz), and the tripling frequency (1270 Hz) of the secondary transmission system is basically the same as the natural frequency of the sixth order, so that the system shows resonance peaks. It should be noted that the amplitude of the acceleration increases significantly after passing the resonance rotational speed area, and the acceleration has a stable change along with the rotating speed. The acceleration of the simulation calculation has no big difference from that of experiment. The two curves are basically identical.
Figure 19 shows that there remains a difference between the theoretical and test results. This could because the model applies Rayleigh damping. However, the vibration mode of the system is different with the changing speed. There are different mode shapes. For the further study, modal damping can be used in future research. Meanwhile, the coupling of air and lubricating oil also influence the system’s acceleration and cause error between the simulation and test [34].
8. Effect of Shaft Stiffness on System
8.1. Shaft Static Deformation Analysis
In order to study the influence of the shaft stiffness on the static deformation of the shaft, the static deformation of the transmission shaft was calculated with the transmission shaft stiffness k_{sh} being 0.25 × k_{sh}, 0.5 × k_{sh}, 1.0 × k_{sh}, 2.0 × k_{sh}, 4.0 × k_{sh}, and 10.0 × k_{sh} (initial stiffness is k_{sh}). Figure 20(a) shows that the changing law of the shaft static deformation is similar under the different shaft stiffness. At the shaftgear coupled node, the deformation is the most obvious. Moreover, along with the increase of the stiffness of the shaft, the deformation of the shaft at the same node becomes smaller and smaller. It is assumed that the deformation of the transmission shaft is D_{max}, showing that the value of deformation of each shaft also gradually decreases along with the increase of the stiffness of the transmission shaft through simulation. This is shown in Figure 20(b).
(a)
(b)
8.2. Effect of the Flexibility of Transmission Shaft on Natural Frequency and Mode Shape
The stiffness of the transmission shaft can affect the bearing loads, natural characteristics, and vibration noises of the reducer. The natural frequency and mode shape of the system are obtained with the stiffness of transmission shaft k_{sh} being 0.25 × k_{sh}, 0.5 × k_{sh}, 1.0 × k_{sh}, 2.0 × k_{sh}, 4.0 × k_{sh}, and 10.0 × k_{sh}. Along with the increase of transmission stiffness, the natural frequencies of the system at each order also increase from equation (15). For instance, with a stiffness of 0.25 × k, the natural frequency of the first order is 64.19 Hz, and with a stiffness of 0.5 × k, the natural frequency of the second order is 90.68 Hz, increased by 41.23%.
If the effect of the shaft stiffness on the system mode shape is considered, the curve of the first mode shapes of shaft 2 with different stiffness (k_{sh} = 0.25 × k_{sh}, k = 10 × k_{sh}) is shown in Figure 21. It is shown that the transmission shaft has the mode shape of translational and torsional coupling vibration. The vector of torsional mode shape is far more than the translational ones. Thus, torsional mode shape is a dominant vibration, which is the main cause of system vibration. So, it necessary to analyze the dominant vibration under different mode orders.
(a)
(b)
Meanwhile, it is found that the change of the shaft stiffness affects the dominant vibration of the system at orders 4, 5, 15, 19, and 20 with simulation. The changes are shown in Table 5.

Table 5 shows that different transmission shaft stiffness values can result in changes of dominant mode shapes of the gear system. For instance, if the shaft stiffnesses are 0.5k_{sh} and 0.25k_{sh}, the dominant mode shape of the coupling system will change from torsional vibration to translation vibration compared with the shaft stiffness of others. In a general way, the torsional vibration occurs more easily in the work time. So, the decrease of stiffness of the transmission shaft can effectively reduce the probability of resonance. The analytical results can provide some guidance for engineering application.
8.3. Bearing Dynamic Load Analysis
In the spur gear system, the distance between two bearings is different, leading to different dynamic loads on bearing and vibration response of gear shaftbearing system under imbalanced excitation. In addition to this, the flexibility of shaft can absorb energy and make an influence on vibration transmission, which also affects the load distribution of both ends of the bearing.
The effects of shaft stiffness on bearing load distribution are studied in this part. NEM is used to solve the model under the input speed of 500 rpm and the load of 100 N·m. The average radial dynamic loads at a and b ends of bearing of each shaft are obtained through simulation. The difference between the force of two bearing ends is defined aswhere and denote the average radial forces at a and b ends of bearing, respectively.
It shows that the change of stiffness of transmission shaft can affect the distribution of bearing force of a and b shafts at input end, middle end, and output end from Figure 22. Along with the increase of transmission shaft’s stiffness, the force of a and b ends tends to be balanced, so if the transmission shaft has greater flexibility, the force difference at two ends will be greater and deflection and warping of shaft can occur in an easier way. To some extent, increasing the stiffness of the shaft can add the fatigue life of the bearing.
8.4. Resonance Characteristics Analysis
Though the parameter of twostage gear transmission system is a linear problem, the dynamic response frequency composition is complicated due to the impact of TVMS and natural frequency of the system. Furthermore, different input speeds can directly affect the frequency composition of the excitation. In order to explore the effect of the input speed on the dynamic load fluctuation of the system, the change of fluctuation amplitude of the dynamic load of gear is calculated at major input speed of the reducer (500–10000 rpm) with different shaft stiffnesses, as shown in Figure 23.
(a)
(b)
(c)
(d)
The figure shows that, along with the increase of input speed, the dynamic load amplitude has great changes. At the initial shaft stiffness (k_{sh}), if the input speed is 668 rpm, 1671 rpm, and 2244 rpm, the twostage gear dynamic load has vibration peaks. This is because, at an input speed of 668 rpm, the firststage fourtime harmonic frequency component (4f_{m1}) is near the natural frequency of the fourthorder, and the system has greater vibration. At an input speed of 1671 rpm, the secondstage 10time harmonic frequency component (10f_{m2}) is near to the natural frequency of eighth order, so that the system has greater vibration energy. If the input speed increases and the input speed is 6493 rpm, the 2f_{m2} (1250 Hz) is not very different from the natural frequency of the third order and the torsional vibration of order three appears, leading to resonance of the system. However, the dynamic load amplitude value greatly decreases after passing the resonance area. In the same time, at an input speed of 7400 and 8498 rpm, the firststage gear has a great response, because the 2f_{m2} at such an input speed is consistent with the natural frequency of orders 14 and 16, making the first shaft vibrate greatly, as shown in Figure 23(a). The frequency component of the peak is listed in Table 6.

With the increase of shaft stiffness, the peak value of the dynamic load of the twostage gear increases greatly as well. For instance, if the input speed increases to 4.0 × k_{sh}, the system has resonance at the torsional direction at speeds of 5252, 6302, and 6350 rpm, and the gear’s dynamic load has sharp changes at the peak amplitudes, as shown in Figure 23(b). It was shown in Section 8.2 that, with the stiffness of transmission shaft k_{sh} increasing to 2k_{sh}, 4k_{sh}, and 10k_{sh}, the fourthorder vibration mode also changes. The system transforms from torsional vibration to bending vibration and has fewer vibration peak points. When the stiffness of the transmission shaft increases, the system’s natural frequency also changes and the torsional mode shape of the system is triggered to form many resonance areas, increasing the possibility of system resonance. Therefore, reasonable design of the stiffness of the transmission shaft can effectively avoid the production of the torsional resonance of the system. Thus, the design and analysis of the gear shaft cannot neglect the effect of the flexibility of the shaft.
9. Conclusion
In this study, a novel modeling method for the coupling vibration analysis of a twostage spur gear transmission system is presented with FEM. The model considered the TVMS, BTVS, and twostage phase relationship and solved the system dynamic characteristics with the NIM. Meanwhile, the influence of shaft stiffness change was investigated. The primary conclusions are as follows:(1)With the influence of bearing timevarying stiffness, there is a bearing passing frequency component f_{b} in lowfrequency location. Moreover, compared with a rigid shaft model, because the flexibility of the transmission shaft can absorb and isolate part of the energy, the dynamic load decreased, and the number of frequency components of the vibration response also decreased.(2)Considering the flexibility of the transmission shaft and the TVMS of the bearing and gear effectively, the frequency component’s distribution in the simulation is basically consistent with the test. Additionally, the effectiveness of the calculation method was verified comparing the simulation data with the experimental ones under the condition of increasing speed.(3)At the initial shaft stiffness (k_{sh}), if the input speed is 668 rpm, 1671 rpm, and 2244 rpm, the twostage gear dynamic loads show amplitude peaks. With the increase of the stiffness of the transmission shaft, the dominant mode shape (torsional mode shapes) of the system increases. Meanwhile, the dominant vibration areas appear divided, and the total amplitude value of the system also increases.(4)With the increase of the stiffness of transmission shaft, the forces at the two ends of the bearing become more consistent. This can provide a theoretical basis for the noise prediction of the gearbox.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors thank Xinjiang University for providing good facilities. The study was supported by the Chinese National Natural Science Foundation (51665054) and Opening Foundation of State Key Laboratory of Mechanical Transmissions (SKLMTKFKT201714).
References
 T. Ouyang, G. Huang, J. Chen, B. Gao, and N. Chen, “Investigation of lubricating and dynamic performances for highspeed spur gear based on tribodynamic theory,” Tribology International, vol. 136, no. 8, pp. 421–431, 2019. View at: Publisher Site  Google Scholar
 H. Liu, C. Zhang, and C. Wang, “Transverseswingtorsional coupling of twostage involute gear transmission system. Nonlinear dynamics modeling and experimental verification,” Vibration and Shock, vol. 36, no. 15, pp. 124–132, 2017. View at: Google Scholar
 Y. Pandya and A. Parey, “Simulation of crack propagation in spur gear tooth for different gear parameter and its influence on mesh stiffness,” Engineering Failure Analysis, vol. 30, pp. 124–137, 2013. View at: Publisher Site  Google Scholar
 W. Jing, A.Q. Zhang, D. Qin et al., “A coupling dynamics analysis method for a multistage planetary gear system,” Mechanism and Machine Theory, vol. 110, pp. 27–49, 2017. View at: Publisher Site  Google Scholar
 K. Umezawa, H. Houjoh, and H. Yoshimura, “The effect of shaft stiffness on the gear vibration,” Bulletin of the JSME, vol. 54, no. 499, pp. 699–705, 1987. View at: Google Scholar
 R. Guilbault, M. S. Lalonde, and M. Thomas, “Nonlinear damping calculation in cylindrical gear dynamic modeling,” Journal of Sound and Vibration, vol. 331, no. 9, pp. 2110–2128, 2012. View at: Publisher Site  Google Scholar
 F. Ren, D.T. Qin, and X.L. Wu, “Research on the dynamics of the double helical gear transmission,” Applied Mechanics and Materials, vol. 456, no. 15, pp. 256–259, 2013. View at: Publisher Site  Google Scholar
 S. Wu, M. J. Zuo, and A. Parey, “Simulation of spur gear dynamics and estimation of fault growth,” Journal of Sound and Vibration, vol. 317, no. 3–5, pp. 608–624, 2008. View at: Publisher Site  Google Scholar
 W. Jing, A.Q. Zhang, G.Q. Wang et al., “A study of nonlinear excitation modeling of helical gears with modification: theoretical analysis and experiments,” Mechanism and Machine Theory, vol. 128, pp. 314–335, 2018. View at: Publisher Site  Google Scholar
 Z. Wang, L. Zhang, Y.Q. Luo, and C.Z. Chen, “Dynamic model of spur gear pair with modulation internal excitation,” International Journal of Rotating Machinery, vol. 2017, Article ID 1264904, 8 pages, 2017. View at: Publisher Site  Google Scholar
 B. Neriya and T. S. Sankar, “The coupled torsional–flexural vibration of a geared shaft system using finite element method,” The Shock and Vibration Bulletin, vol. 55, no. 3, pp. 13–25, 1985. View at: Google Scholar
 S. V. Modak, T. K. Kundra, and B. C. Nakra, “Use of an updated finite element model for dynamic design,” Mechanical Systems and Signal Processing, vol. 16, no. 23, pp. 303–322, 2002. View at: Publisher Site  Google Scholar
 H. N. ÖzgÜven and Z. L. Özkan, “Whirl speeds and unbalance response of multi bearing rotors using finite elements,” Journal of Vibration, Acoustics, Stress and Reliability in Design, vol. 106, no. 1, pp. 72–79, 1984. View at: Publisher Site  Google Scholar
 T. Ouyang, Z. Su, S. Li, G. Huang, and N. Chen, “Experimental and numerical investigations on dynamic characteristics of gearrollerbearing system,” Mechanism and Machine Theory, vol. 140, no. 10, pp. 730–746, 2019. View at: Publisher Site  Google Scholar
 M. Kubur, A. Kahraman, D. M. Zini, and K. Kienzle, “Dynamic analysisof a multishaft helical gear transmission by finite elements: model and experiment,” Journal of Vibration and Acoustics, vol. 126, no. 3, pp. 398–406, 2004. View at: Publisher Site  Google Scholar
 Q.L. Jiang, D.Z. Wu, and S.G. Tan, “Development and application of a model for coupling gearedrotors system,” Journal of Vibration Engineering, vol. 23, no. 3, pp. 254–259, 2010. View at: Google Scholar
 Y. Guo, T. Eritenel, T. M. Ericson, and R. G. Parker, “Vibroacoustic propagation of gear dynamics in a gearbearinghousing system,” Journal of Sound and Vibration, vol. 333, no. 22, pp. 5762–5785, 2014. View at: Publisher Site  Google Scholar
 H. Ma, X. Pang, R. Feng, and B. Wen, “Evaluation of optimum profile modification curves of profile shifted spur gears based on vibration responses,” Mechanical Systems and Signal Processing, vol. 7071, pp. 1131–1149, 2016. View at: Publisher Site  Google Scholar
 L. Qian, J.Y. Tang, S.Y. Chen, and Y. Liu, “Dynamic modeling of a onestage gear system by finite element method and the dynamic analysis in high speed,” Journal of Mechanical Engineering, vol. 52, no. 17, pp. 155–161, 2016. View at: Publisher Site  Google Scholar
 S. Ankur, C. Manoj, and P. Anand, “Measurement of FRFs of coupled geared rotor system and the development of an accurate finite element model,” Mechanism and Machine Theory, vol. 123, no. 1, pp. 66–75, 2018. View at: Publisher Site  Google Scholar
 J. Wang, Numerical and Experimental Analysis of Spur Gears in Mesh, Curtain University of Technology, Bentley, Australia, 2003.
 A. Tatar, C. W. Schwingshackl, and M. I. Friswell, “Dynamic behaviour of threedimensional planetary geared rotor systems,” Mechanism and Machine Theory, vol. 134, pp. 39–56, 2019. View at: Publisher Site  Google Scholar
 J. Li, Theoretical Study of Space Beam Element and its Application in Bridges, Changsha University of Science and Technology, Changsha, China, 2014.
 C.B. Hu, Y.G. Wang, and D.S. Ling, “Effect of Rayleigh damping physical essence and parameters on dynamic response,” Journal of Zhe Jiang University (Engineering Science), vol. 51, no. 6, pp. 1196–1203, 2017. View at: Google Scholar
 Z.Z. Qiao, J.X. Zhou, and X.F. Zhang, “Research on finite element modeling method of twostage spur gear transmission system under multisource timevarying excitation,” Journal of Vibration and Shock, vol. 38, no. 15, pp. 182–189, 2019. View at: Google Scholar
 R. Chen, J. Zhou, and W. Sun, “Dynamic characteristics of a planetary gear system based on contact status of the tooth surface,” Journal of Mechanical Science and Technology, vol. 32, no. 1, pp. 69–80, 2018. View at: Publisher Site  Google Scholar
 L.H. Chang, Z.X. He, and G. Liu, “Finite element method for dynamic modeling of parallel shaft gear transmission system,” Journal of Vibration and Shock, vol. 35, no. 20, pp. 47–53, 2016. View at: Google Scholar
 H. Ma, R.Z. Song, X. Pang, and B. Wen, “Fault feature analysis of a cracked gear coupled rotor system,” Mathematical Problems in Engineering, vol. 2014, Article ID 832192, 22 pages, 2014. View at: Publisher Site  Google Scholar
 S. L. Harris, “Dynamic load on the teeth of spur gears,” Proceedings of the Institution of Mechanical Engineers, vol. 172, no. 1, pp. 87–112, 1958. View at: Publisher Site  Google Scholar
 Y. Yang, L.Y. Cao, H. Li, and Y. Dai, “Nonlinear dynamic response of a spur gear pair based on the modeling of periodic mesh stiffness and static transmission error,” Applied Mathematical Modelling, vol. 72, no. 6, pp. 444–469, 2019. View at: Publisher Site  Google Scholar
 B. Yuan, S. Chang, G. Liu, L.H. Chang, and L. Liu, “Quasistatic analysis based on generalized loaded static transmission error and dynamic investigation of widefaced cylindrical geared rotor systems,” Mechanism and Machine Theory, vol. 134, no. 8, pp. 289–296, 2019. View at: Publisher Site  Google Scholar
 Y. Li, K. Ding, G.L. He, and H. Lin, “Modulation sidebands of the vibration signal of gearbox,” Journal of Mechanical Engineering, vol. 54, no. 5, pp. 105–112, 2018. View at: Publisher Site  Google Scholar
 H. Cao, Y.M. Li, and Z.J. He, “Time varying bearing stiffness and vibration response analysis of high speed rolling bearingrotor systems,” Journal of Mechanical Engineering, vol. 50, no. 15, pp. 73–81, 2014. View at: Publisher Site  Google Scholar
 Xiaowei Cheng, Xiaohui Shi, and Quan Shi, “Study on the error wheel and experimental research of secondary gear transmission,” Journal of Chongqing University of Technology (Natural Science), vol. 24, no. 11, pp. 19–24, 2010. View at: Google Scholar
Copyright
Copyright © 2019 Zizhen Qiao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.