Next Article in Journal
Investigation on Unsteady Cavitation Flow and Excited Pressure Fluctuations in Regulating Valve
Previous Article in Journal
A CNN-Based Method for Counting Grains within a Panicle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Modeling and Analysis of Nonlinear Compound Planetary System

1
College of Mechanical and Vehicle Engineering, Chongqing University, Chongqing 400030, China
2
Structural Dynamics and Acoustic Systems Laboratory, University of Massachusetts Lowell, One University Avenue, Lowell, MA 01854, USA
*
Author to whom correspondence should be addressed.
Machines 2022, 10(1), 31; https://doi.org/10.3390/machines10010031
Submission received: 19 October 2021 / Revised: 10 November 2021 / Accepted: 12 November 2021 / Published: 1 January 2022
(This article belongs to the Section Vehicle Engineering)

Abstract

:
In the vehicle composite planetary gear transmission system, nonlinear excitations such as time-varying meshing stiffness, backlash and comprehensive error would lead to large vibration and noise, uneven load distribution, unstable operation and other problems. To address these issues, this work focuses on compound planetary gears and develops the bending-torsion coupling nonlinear dynamic model of the system based on the Lagrange equation. There are internal and external multi-source excitations applied to the system. This model is used to study the bending-torsion coupling meshing deformation relationship of each meshing pair along with the translational and torsional directions. The natural frequencies and vibration modal characteristics of the system are extracted from the model, and the influence of rotational inertia and meshing stiffness on the inherent characteristics of the system are studied. The coupling vibration characteristics of the system under operating condition are analyzed in terms of the inherent characteristics and time–frequency characteristics of the system. The simulation results exhibit that the planetary gear system has three modes. The change in natural frequency trajectory has two phenomena: modal transition and trajectory intersection. The main frequencies include engine rotating frequency, meshing frequency and its double frequency, and the rotation frequency and harmonic frequency of the engine have a great influence on the vibration response of the system. Finally, the virtual prototype of the composite planetary system is used to verify the accuracy of the established model from speed, inherent characteristics, meshing force and frequency composition.

1. Introduction

The planetary gear transmission system has, as its characteristics, a large transmission ratio, stable transmission, strong carrying capacity, high unit volume, mass power density, coaxial input and output and strong impact resistance and power shunt, which is widely used in various industrial fields. With the change in the number and phase of the planetary gear, as well as the diversification of the engagement with the gear ring sun gear, the design of the planetary gear system becomes complex and diverse. There are some nonlinear excitations in the gear transmission system, such as time-varying meshing stiffness, backlash and comprehensive transmission error, which cause the vibration and noise of the system, and directly affect the reliability and safety of the system.
Many scholars have made an in-depth study of the vibration characteristics of planetary systems. Ambarisha et al. [1] applied the lumped parameter and finite element model to investigate the complex and nonlinear dynamic behavior of spur planetary gears. Sun et al. [2,3,4] studied the strongly nonlinear dynamic behavior of planetary gear systems with multiple clearances. Kim et al. [5] offered the planetary gear dynamic model, assuming that the position of the contact line is determined by the average angular motion of the gear. Zhou et al. [6] employed numerical methods to analyze the planetary gear transmission system at high speed and light weight. Liu et al. [7,8,9] made a study of the planetary gear dynamic model considering time-varying stiffness, gear meshing error and gyroscopic effect. Kahraman et al. [10] established a planetary gear pure torsion dynamics model to obtain the planetary gear’s dynamics in which the natural frequency is compared with the dynamic model of the torsion translation planetary gear. Peng et al. [11] first used the meshing phase to study planetary gears with faults. Qiu et al. [12] reviewed the dynamics of the planetary gear transmission system of wind turbines from the aspects of research status, dynamic optimization design and development direction. Xiang et al. [13] discovered the non-linear dynamic characteristics of the system under the change in the planetary gear’s supporting stiffness and excitation frequency. Bozca et al. [14] presented a multi-body torsional vibration dynamic model of the rigid-flexible coupling planetary gear system by using the state-space model and considering the influence of box flexibility. Zhang et al. [15] raised a dynamic model of planetary transmission considering the flexibility of the ring gear and identified its inherent characteristics and relevant test verification. Qin et al. [16] advanced the coupling dynamics model of the wind turbine gear transmission system and examined the natural frequency, vibration response, dynamic meshing force and rolling bearing dynamic bearing force of the wind turbine gear transmission system. Dou et al. [17,18] studied the frequency coupling and coupling resonance of the composite planetary transmission system under complex excitations. Wang et al. [19] studied the time–frequency characteristics of planetary gear systems by establishing a flexible multi-body dynamics model. Li et al. [20] inspected the spectral characteristics of vibration signals of the compound planetary transmission system. Zhu et al. [21] explored the dynamic response characteristics of planetary transmission systems by gear modification. Inalpolat et al. [22] investigated the influence of system input speed on the dynamic response amplitude of multi-stage planetary gear systems and verified it through experiments. Chaari et al. [23] utilized the iterative method to explore the dynamic response characteristics of planetary transmission systems under nonlinear factors. Zhang et al. [24] established the dynamic model of planetary gear systems with sliding bearing and flexible structure, and the influence of different radial sliding bearings on load sharing characteristics was researched. Gu et al. [25] proposed a lumped parameter model of original planetary gear considering planetary position error and simulated their effects on quasi-static and dynamic load sharing among planets. Li et al. [26] observed the influence law of each excitation factor from the perspective of bifurcation characteristics, and researched whether the change in backlash would affect whether the system was periodic or chaotic. Ryali et al. [27] proposed a three-dimensional dynamic load distribution model of planetary gear sets. Mo et al. [28] inspected the load-sharing characteristics of flexible support when the sun gear is floating and normal. Sanchez-Espiga et al. [29] proposed a numerical approach for calculating the load distribution of planetary transmission. Iglesias et al. [30] found that the tooth load is much lower than the planetary non-uniform load in the transmission with and without defects.
These scholars have done a large amount of work on the dynamic modeling of single-stage planetary row, nonlinear dynamic analysis and load sharing characteristics. Few people have studied the dynamic characteristics of multi-stage planetary gear system under internal and external nonlinear excitation. However, these nonlinear factors bring about self-excited motion, which has a great influence on the system. In order to study the vibration problem caused by nonlinear excitation in the composite planetary system, this paper establishes the bending-torsion coupling dynamic model of the composite planetary transmission system, considering the external and internal excitations such as engine harmonic excitation, time-varying stiffness, gear time-varying phase, comprehensive error and dynamic backlash. The meshing deformation relationship of meshing pair in translational and torsional directions was analyzed in detail, and the correctness of the model is verified by the virtual prototype from the aspects of rotational speed, inherent characteristics, meshing force and frequency composition. The vibration characteristics of the nonlinear system are revealed by analyzing the natural frequency, which varies with stiffness and inertial trajectory, as well as the time–frequency characteristics of the system response.

2. Dynamic Model for Compound Planetary System

2.1. Model Description

The schematic of the compound planetary system is shown in Figure 1, which includes three pairs of planetary chains composed of sun gear s , planet carrier c , long planet gear a , short planet gear b , the small gear ring r 1 and the big gear ring r 2 . In addition, the brakes C 1 and C 2 , and the clutch C 3 , control different gears. The meshing relations in the system are equivalently simplified into a mass-spring-damping model, as shown in Figure 1. In this paper, the brake C 1 is engaged, and the brake C 2 and the clutch C 3 are disengaged. The small ring gear r1 sends power to the long planetary gear a, while a transmits the power to the sun gear s and the short planetary gear b; b transfers the power to the large ring gear r2 and the planet carrier c; finally, the power is output from both ends of the planet carrier.
To describe the meshing relationship between planetary gear systems more accurately, three coordinate systems are used in this paper: the first is the absolute coordinate system X _ O _ Y , which is used to describe the vibration displacement of all central rotating parts; the second is the involved coordinate system X i _ O i _ Y i , which describes the translational vibration displacement and which rotates synchronously with the planetary shelf; the third is the relative coordinate system x a j _ O a j _ y a j of the jth planetary wheel in the planetary row of a or b, whose origin is located in the center of each planetary wheel, and its coordinate axis is orthogonally along the radial and circumferential directions of the planetary carrier. In the dynamic model of the composite planetary gear system established in this paper, the vibration displacement of planetary gear is based on their respective relative coordinate systems x a j _ O a j _ y a j . The bending–torsion coupling mechanical model of the composite planetary gear system is shown in Figure 2. k m s a j , k m r a j , k m r b j   and   k m a b j are the meshing stiffness of the meshing pair in the torsional direction; θ c , θ r 1 , θ r 2 ,   and   θ s , respectively, represent the torsional displacement relative to the fixed reference frame X _ O _ Y ;   θ a j   and   θ b j indicate the relative torsional displacement of the planet wheel a j and the planet wheel b j in the relative coordinate system x a j _ O a j _ y a j , x b j _ O b j _ y b j , where j = 1 ,   2 ,   3 ; k b p x , k b q x , k b p y , and   k b q y express the supporting stiffness of p and q in the x direction and y direction in their respective coordinate systems, where p = s , c , r 1   and   r 2 , respectively, mean the sun wheel, planet rack, small ring gear and large ring gear in the composite planetary gear system; q = a j , b j represent the long planet wheel and the short planet wheel in the composite planetary gear.

2.2. The Deformation of Meshing Planetary Gears

The expression of meshing deformation of a pair of gear teeth is related to the meshing model, coordinate selection, positive rotation direction selection and meshing line selection. Meshing deformation is defined along the tangential direction of the base circle between meshing gears. There are four meshing relationships in Figure 3 in this paper; the meshing deformation relationship between planetary wheels is selected for research and the meshing compression is positive in the direction of deformation.
There are two different mesh lines, A B and A B , in Figure 3, though their meshing deformation calculation method is the same. This paper selects the A B as the mesh line. r b a and r b b are the base circle radius of the two gears. The mesh deformation caused by pure torsion of the gear can be described by the length variation between A and B points, which can be expressed as Equation (1):
δ a b = r b a θ a + r b b θ b
The projection relation of the translational displacement of the planetary wheel a j on the meshing line is shown in Figure 4a. The absolute translational displacement of the planetary rack L c is the vector sum of translational displacements x c and y c in the fixed coordinate, which is also the involvable displacement of the planetary wheel; x c and y c are projected to the orthogonal directions of x a j and y a j ; the relative displacement L a of the planet wheel a j is the vector sum of the translational displacements x a j and y a j in the relative coordinate x a j O a j y a j ; the absolute translational displacements L a can be regarded as the vector sum of its relative translational displacements L a and the involved translational displacements L c ; finally, the total deformation of x a j and y a j in two orthogonal directions is projected along the meshing line EF to obtain the change in the planetary wheel a j on the meshing line EF induced by the translation displacement σ a b b , as shown in Equation (8). Similar to the planetary gear a j , the projection relationship of the planetary gear b j ’s translation displacement on the meshing line is shown in Figure 4b.
In Figure 4, r c a and r c b , respectively, express the distance between the theoretical rotation center of the planet wheel a j and b j and the rotation center of the planet rack.   ψ a j and ψ b j represent the position angle of the planet wheels a j and b j in the fixed coordinate system X O Y , which are calculated as shown in Equations (2) and (3). β a and β b indicate the angle between the centerline of the planetary wheel a j and b j , meshing with each other and their respective coordinate axes; their relationship is shown in Equations (4) and (5). γ a and γ b indicate the angle between the meshing line of the a j and b j intermeshing planetary wheels and their respective coordinate axes X a j and X b j , which are calculated as shown in Equations (6) and (7); the other symbols have the same meaning as previously stated.
ψ a j = ω c t + 2 π ( j 1 ) / n + φ a
ψ b j = ω c t + 2 π ( j 1 ) / n + φ b
in which n is the number of planetary gears; φ a and φ b denote the initial circumferential installation angle of planetary gear a j and b j .
β a + β b = ψ b j ψ a j
s i n β a r c b = s i n β b r c a
γ a = β a + α π / 2
γ b = π / 2 ( β b α )
The total deformation σ a b a of planetary gear a j on the meshing line EF is:
σ a b a = ( y c s i n ψ a j + x c c o s ψ a j + x a j ) c o s γ a + ( y c c o s ψ a j x c s i n ψ a j + y a j ) s i n γ a
The total deformation σ a b a of planetary gear b j on the meshing line EF is:
σ a b b = ( y c s i n ψ b j + x c c o s ψ b j + x b j ) c o s γ b ( y c c o s ψ b j x c s i n ψ b j + y b j ) s i n γ b
The total meshing deformation σ a b caused by the translation displacement in the meshing line direction of planetary gear a j and planetary gear b j is:
σ a b = σ a b a + σ a b b
Similarly, in the direction of the outer meshing line between the sun wheel and the planetary gear, the meshing deformation δ s a and the total meshing deformation in external meshing line σ a s are shown in Equations (11) and (13). In the direction of the internal meshing line between the gear ring and the planetary gear, the meshing deformation δ r a and the total meshing deformation in external meshing line δ a r are shown in Equations (14) and (16).
δ s a = λ s a ( r b s ( θ s θ c ) + r b a θ a ) = σ s + σ a s
σ s = x s s i n ( ψ a j + α ) y s c o s ( ψ a j + α )
σ a s = x c s i n ( ψ a j + α ) + y c c o s ( ψ a j + α ) x a j s i n α + y a j c o s α
δ r a = λ r a ( r b a θ a + r b r ( θ r θ c ) ) = σ r + σ a r
σ r = x r s i n ( ψ a j α ) + y r c o s ( ψ a j α )
σ a r = x c s i n ( ψ a j α ) y c c o s ( ψ a j α ) x a j s i n α y a j c o s α

2.3. Internal Nonlinear Excitation

The internal nonlinear excitation would be caused by the comprehensive error, the dynamic backlash and the time-varying stiffness.
The eccentricity error (EE) between the sun and the planets is:
{ e s a i = E s s i n ( ω s t ε s + ψ s a i + α ) e a s i = E a i s i n ( ω a i t ε a i + α )
The EE between gear ring r 1 and planetary gear a is:
{ e r a i = E r s i n ( ω r t + ε r ψ r a i + α ) e a r i = E a i s i n ( ω a i t + ε a i + α )
The EE between gear ring r 2 and planetary gear b is:
{ e r b i = E r s i n ( ω r t ε r + ψ r b i + α ) e b r i = E b i s i n ( ω b i t ε b i + α )
The EE between planetary wheel a and planetary wheel b is:
{ e a b i = E a i s i n ( ω a i t + ε a i + α + γ a i ) e b a i = E b i s i n ( ω b i t + ε b i + α γ b i )
The EE caused by planetary carrier to internal and external meshing pairs is:
{ e c _ s a i = E c s i n ( ω c t + ψ j i + α ) e c _ r a i = E c s i n ( ω c t ψ j i + α ) e c _ a b i = E c s i n ( ω c t + ψ j i + α ) e c _ r b i = E c s i n ( ω c t ψ j i + α )
where   E l and E q i represent the amplitude of manufacturing eccentricity error, ε l   and   ε q i stand for the initial phase of eccentricity error corresponding, e l q i ,   e q l i   and   e q q i are the eccentricity error of sun gear, gear ring and planetary gear, l = s , r   and   c   are, respectively, sun gear, gear ring and planetary carrier; q = a , b are the two-stage planetary gear;   and ,   γ b i are the angle between the theoretical installation position of planetary wheel 1 p i in row a and 2 p i in row b ; ψ j i is the position angle between the central component j and the ith planetary wheel.
The assembly eccentricity error (AEE) between the sun gear and planetary wheel a is:
{ a s a i = A s s i n ( φ s + ψ s a i + α ) a a s i = A a i s i n ( φ a i + α )
The AEE between gear ring 1 and planetary gear a is:
{ a r a i = A r s i n ( φ r ψ r a i + α ) a a r i = A a i s i n ( φ a i + α )
The AEE between gear ring 2 and planetary wheel b is:
{ a r b i = A r s i n ( φ r + ψ r b i + α ) a b r i = A b i s i n ( φ b i + α )
The AEE caused by planetary carrier to internal and external meshing pairs is:
{ a c _ s a i = A c s i n ( φ c + ψ j i + α ) a c _ r a i = A c s i n ( φ c ψ j i + α ) a c _ a b i = A c s i n ( φ c + ψ j i + α ) a c _ r b i = A c s i n ( φ c ψ j i + α )
The AEE between planetary wheel a and planetary wheel b is:
{ a a b i = A a i s i n ( φ a i + α + γ a i ) a b a i = A b i s i n ( φ b i + α γ b i )
where A l   and   A q i stand for the assembly eccentricity error of amplitude, φ l   and   φ q i represent the eccentric error of the initial phase, a l q i ,   a q l i   and   a q q i are, respectively, projection to the corresponding meshing line wheel, gear ring, the sun and the planets round eccentric error;   l = s , r   and   c   are, respectively, sun gear, gear ring and planetary carrier; q = a , b   are the two-stage planetary gear.
The tooth profile equivalent meshing error is expressed as a harmonic function, and its excitation frequency is the meshing frequency of the tooth pair, as shown in Equation (27):
{ b s a i = B s a i s i n ( ω m s a t + 2 π γ s a i ) b r a i = B r a i s i n ( ω m r a t + 2 π γ r a i ) b r b i = B r b i s i n ( ω m r b t + 2 π γ r b i ) b a b i = B a b i s i n ( ω m a b t + 2 π γ a b i )
Here, b j i is the tooth profile error along the direction of the meshing line, B j i is the manufacturing error amplitude, ω m j is the meshing frequency, γ j i is the initial phase, i = 1 , 2 , 3 when j equals s a , r a , r b , a b , which are, respectively, sun and planetary wheel p a i , small gear and planetary wheel p a i , large gear and planetary wheel p b i and planetary wheel p a i .
Based on the mathematical model of the above, the comprehensive error along the engagement line in the composite planetary system is shown in Equation (28):
{ e s a i = e s a i + a s a i + e c _ s a i + a c _ s a i + e a b i + a a b i + b s a i e r 1 a i = e r a i + a r a i + e c _ r a i + a c _ r a i + e b a i + a b a i + b r a i e r 2 b i = e r b i + a r b i + e c _ r b i + a c _ r b i + e b a i + a b a i + b r b i
In summary, when considering the bending–torsion coupling mechanical model of the composite planetary gear system, the total deformation on the meshing line is the superposition of the torsional vibration displacement deformation of each gear meshing pair, the deformation induced by the translational vibration displacement and the comprehensive meshing error.
The total meshing deformation (TMD) of the solar wheel and the planetary wheel a j is:
ε s a j = δ s a j + σ s a j + e s a a j
The TMD between the pinion ring and the a j the planetary gear is:
ε r a j = δ r a j + σ r a j + e r a a j
The TMD between the large gear ring and the b j the planetary gear is:
ε r b j = δ r b j + σ r b j + e r b b j
The TMD ε a b j of the a j planetary gear and the b j planetary gear is:
ε a b j = δ a b j + σ a b j + e a b j
The nonlinear meshing forces of each meshing gear in the x and y directions on the meshing line and in the respective coordinate systems.
{ F m p j = k m p j ( t ) f ( ε p j , b ) + c m ε ˙ p j F m p j x = F m p j s i n ( ψ p j + α ) F m p j y = F m p j c o s ( ψ p j + α ) ( j = 1 , 2 , 3 )
when p equals s a , r   and   r b are, respectively, sun and planetary wheel, small gear and planetary wheel a j and large gear and planetary wheel b j .
The nonlinear meshing force of the planetary wheel a j and planetary wheel b j in the direction is:
{ F m a b j = k m a b j ( t ) f ( ε a b j , b ) + c m ε ˙ a b j F m a b j x = F m a b j s i n α F m a b j y = F m a b j c o s α ( j = 1 , 2 , 3 )
Bearing support in the translation direction is considered in the composite planetary gear system. Each component is subjected to bearing support reaction in x and y directions.
{ F b q x = k b q x x q + c b q x x ˙ q F b q y = k b q y y q + c b p y y ˙ q
when q equals s , r 1 , r 2 , a j   and   c are, respectively, sun gear, small gear, large gear, planetary wheel a j , planetary wheel b j and planetary frames.
Based on the above deformation relationship and elastic mechanics theory, considering the time-varying stiffness k m calculated by the Shichuan formula and tooth side clearance function f ( δ , b ) , the meshing force F m can be obtained by Equation (36).
F m = k m ( t ) f ( δ , b ) + c m δ ˙
The time-varying stiffness of teeth is:
k m ( t ) = F l o a d δ ( t ) b
in which F l o a d is the normal load on the tooth surface; b is the tooth width; δ ( t ) is the total tooth deformation.
Comprehensive deformation δ i in the direction of the single tooth meshing line can be expressed as the linear superposition of deformation of rectangular bending part δ B r , local deformation δ G , shear deformation δ S and trapezoidal bending deformation δ B t .
δ i = δ B r + δ B t + δ S + δ G
when i = p , n , respectively, are the single tooth deformation of driving wheel and driven gear.
The bending deformation of the rectangular part can be expressed as:
δ B r = 12 F N c o s 2 ω x / ( E b s F 3 [ h x h r ( h x h r ) + h r 3 / 3 ]
The bending deformation of trapezoidal part can be expressed as:
δ B t = 6 F N c o s 2 ω x / ( E b s F 3 ) [ h i h x h i h r ( 4 h i h x h i h r ) 2 ln h i h x h i h r 3 ] ( h i h r ) 3
The meshing shear deformation can be expressed as:
δ S = 2 ( 1 + ν ) F N c o s 2 ω 2 / E b s F [ h r + ( h i h r ln h i h r h i h r ) ]
The local deformation caused by the inclination of the substrate can be expressed as:
  δ G = 24 F N h x 2 c o s 2 ω 2 / π E b s F 2
here
h i = ( h s F h r s k ) / ( s F s k )
h = r k 2 ( s k 2 ) 2 r r 2 ( s F 2 ) 2
h x = r x c o s ( α x ω x ) r r 2 ( s F 2 ) 2
when r g r F , z 2 ( 1 x ) ( 1 c o s α 0 ) :
s F = 2 r F s i n { ( ( π + 4 x t g α 0 ) / 2 z ) + i n v α 0 i n v α F }
α F = arccos ( r g r F )
h r = r F 2 ( s F 2 ) 2 r r 2 ( s F 2 ) 2
when r g > r F , z < 2 ( 1 x ) ( 1 c o s α 0 ) :
s F = 2 r g s i n { ( π + 4 x t g α 0 ) / 2 z + i n v α 0 }
h r = r g 2 ( s F 2 ) 2 r r 2 ( s F 2 ) 2
where c represents Poisson’s ratio; z means number of teeth; x is displacement coefficient; r g   stands for the base circle;   r r denotes the root circle;   r x expresses the pitch circle; r k is the addendum circle;   r F means radius of effective tooth root circle participating in meshing; α 0 is the standard pressure angle; α x represents meshing angle of the gear.
When a pair of gears mesh, in addition to the deformation caused by each single tooth, there is also the local deformation δ p v caused by the meshing contact stress F n of the tooth pair, which can be expressed as
δ p v = 4 ( 1 ν 2 ) F n / π E b
Combined with the above formulas, the comprehensive deformation along the meshing direction at the meshing point is
δ ( t ) = δ p + δ n + δ p v

2.4. External Nonlinear Excitation

This article uses the four-stroke 12-cylinder V-type engine with rated power of 1000 kW and rated speed of 2300 r/min. Under the rated condition, the excitation torque characteristics of the single cylinder and the harmonic analysis of the output torque of the engine are shown in Figure 5a,b. It can be seen that the amplitude characteristics of the output torque of the engine are mainly composed of six-order harmonics, except for the average torque. In the nonlinear dynamic model, the torque can be equivalent to the superposition of the average amplitude and the harmonic function of the six-order harmonics. The expression is shown in Equation (53):
M i n = M 0 + μ M 0 s i n ( 2 π f e t + φ e )
where M 0 is the average output torque of the engine, f e is the torque fluctuation frequency, φ e is the initial phase and μ is the ratio of the torque amplitude to the torque M 0 .

2.5. Nonlinear Flexural and Torsional Coupling Dynamics Equation of the System

For multi-degree-of-freedom complex dynamical systems, the Lagrange equation is usually used to derive the system vibration differential equation.
The bending–torsional coupling vibration differential equation of sun gear is:
{ m s x ¨ s + j = 1 3 F m s a j x + F b s x = 0 m s y ¨ s j = 1 3 F m s a j y + F b s y = 0 J s θ ¨ s j = 1 3 r b s F m s a j = 0  
The bending–torsional coupling vibration differential equation of a small gear ring is:
{ m r 1 x ¨ r 1 + j = 1 3 F m r a j x + F b r 1 x = 0 m r 1 y ¨ r 1 j = 1 3 F m r a j y + F b r 1 y = 0 J r 1 θ ¨ r 1 + j = 1 3 r b r 1 F m r a j = T i n
The bending–torsional coupling vibration differential equation of the long planetary wheel is:
{ m a j ( x ¨ a j + x ¨ c c o s ψ a j + y ¨ c s i n ψ a j ) m a j ( 2 ω c y ˙ a j + ω c 2 x a j ) F m s a j s i n α + F m r a j s i n α + F m a b j x + F b a j x = 0 m a j ( y ¨ a j x ¨ c s i n ψ a j + y ¨ c c o s ψ a j ) + m a j ( 2 ω c x ˙ a j ω c 2 y a j ) F m s a j c o s α ( j = 1 , 2 , 3 ) + F m r a j c o s α + F m a b j y + F b a j y = 0 J a j ( θ ¨ c + θ ¨ a j ) + r b a j F m s a j r b a j F m r a j + r b a j F m a b j = 0
The bending–torsional coupling vibration differential equation of planetary carrier is:
{ m c x ¨ c + i = a , b j = 1 3 m i j x ¨ c + i = a , b j = 1 3 m i j ( x ¨ i j c o s ψ i j y ¨ i j s i n ψ i j ) + ω c 2 i = a , b j = 1 3 m i j ( x i j c o s ψ i j + y i j s i n ψ i j ) + 2 ω c i = a , b j = 1 3 m i j ( x ˙ i j s i n ψ i j y ˙ i j c o s ψ i j ) j = 1 3 F m s a j x j = 1 3 F m r a j x j = 1 3 F m r b j x + F b c x = 0 m c y ¨ c + i = a , b j = 1 3 m i j y ¨ c + i = a , b j = 1 3 m i j ( x ¨ i j s i n ψ i j + y ¨ i j c o s ψ i j ) + ω c 2 i = a , b j = 1 3 m i j ( x i j s i n ψ i j y i j c o s ψ i j ) + 2 ω c i = a , b j = 1 3 m i j ( x ˙ i j c o s ψ i j y ˙ i j s i n ψ i j ) j = 1 3 F m s a j y j = 1 3 F m r a j y j = 1 3 F m r b j y + F b c y = 0 J c θ ¨ c + i = a , b j = 1 3 m i j r c 2 θ ¨ c + j = 1 3 J a j ( θ ¨ c + θ ¨ a j ) j = 1 3 J b j ( θ ¨ c + θ ¨ b j ) ± j = 1 3 r b s F m s a j j = 1 3 r b r 1 F m r a j + j = 1 3 r b r 2 F m r b j = T o u t
The bending–torsional coupling vibration differential equation of the big gear ring is:
{ m r 2 x ¨ r 2 + j = 1 3 F m r b j x + F b r 2 x = 0 m r 2 y ¨ r 2 j = 1 3 F m r b j y + F b r 2 y = 0 J r 2 θ ¨ r 2 j = 1 3 r b r 2 F m r b j = T b r a k e
The bending-torsional coupling vibration differential equation of short planetary wheel is:
{ m b j ( x ¨ b j + x ¨ c c o s ψ b j + y ¨ c s i n ψ b j ) m b j ( 2 ω c y ˙ b j + ω c 2 x b j ) + F m a b j x + F m r b j s i n α + F b b j x = 0 m b j ( y ¨ b j x ¨ c s i n ψ b j + y ¨ c c o s ψ b j ) + m b j ( 2 ω c x ˙ b j ω c 2 y b j ) ( j = 1 , 2 , 3 ) F m a b j y + F m r b j y c o s α + F b b j y = 0 J b j ( θ ¨ c + θ ¨ b j ) r b b j F m a b j r b b j F m r b j = 0
In which m p is the equivalent mass of each component; J p is the equivalent moment of inertia of each inertial component; p = s ,   c ,   a j ,   b j ,   r 1   and   r 2 , respectively, represent each rotational component in the composite planetary gear; T i n is the nonlinear excitation torque, considering the engine fluctuation torque M i n ; T o u t is the load torque; T b r a k e is the braking torque.

3. Dynamic Characteristic Analysis

3.1. The Inherent Characteristics

The modulus of all gears is 4 and the displacement coefficient is 0; the pressure angle is 20 rad; the meshing clearance is 0.06 mm; the base circular tooth thicknesses of planet gear a and planet gear b are 7.7592 mm and 7.8679 mm, which are used to calculate the phase difference between gears. The transmission ratio can be calculated according to the number of teeth of each gear; the main parameters as shown in Table 1.
Based on the established dynamic model of the composite planetary gear system, the natural frequency and vibration mode of the system can be obtained. The analysis of its inherent characteristics can obtain the resonant frequency, vibration mode characteristics and vibration mode. The influence of the rotational inertia of each component of the system and the meshing stiffness between components on the natural frequency is further analyzed. The characteristics of the natural frequency and vibration mode of the composite planetary row under the working condition are analyzed, which can help to better understand the vibration law of the planetary gear system.
Figure 6 shows that the vibration modes of the composite planetary gear system can be divided into three types of modes: planetary gear vibration, global vibration and coupled vibration. When the system displays the planetary gear vibration, namely that only the planetary gear has torsional vibration, other rotating components do not vibrate, as shown in Figure 7a. When the system displays the global vibration, the vibration generated by the central rotating parts and planetary gear at all levels as well as the amplitude of the same type of planetary gear, are the same, as shown in Figure 7b. When the system displays the coupled vibration, except for the obvious characteristics of planetary wheel vibration and global vibration, the rest can be regarded as coupled vibration, as shown in Figure 7c.

3.2. The Natural Frequency Varies with the Trajectory of Stiffness and Inertia

The natural frequency is mainly related to the structural parameters of the system. It was found in the study that, with the change in parameters, the trajectory will appear as modal transition and trajectory intersection.
The trajectory of the natural frequencies of each order of the system changing with the meshing stiffness k m r a j is shown in Figure 8. It can be seen that ω 2 and ω 3 are the double roots of the same frequency before point A, and their trajectories are the same. At point A, there is a modal transition, ω 2 and ω 3 alter rapidly with different frequencies, respectively, then ω 3 changes with single frequency in the AB and ω 4 changes with a single root frequency before point B. ω 3 and ω 4 modify in the same trajectory after the trajectory intersection at B until the mode transition occurs again at C, and ω 3 and ω 4 transform rapidly with different frequencies, respectively. In this case, the complex trajectory of natural frequencies is induced by the third-order natural frequency ω 3 . The trajectories of ω 5 , ω 6 , ω 8   and   ω 9 always coincide, which is caused by the double roots under the same natural frequency. Trajectories of the meshing stiffness k m s a j , k m r b j and k m r b j with natural frequency are similar to that of k m r a j , as shown in Figure 9, Figure 10 and Figure 11.
Shown in Figure 12 are the trajectory changes in natural frequencies of each order when the rotational inertia J r 1 of the gear ring increases from 0.01   k g · m 2 to 0.5   k g · m 2 . It can be seen that there are more complex modal transitions and trajectory intersection phenomena when the natural frequency changes. The trajectory intersection of high-order frequencies ω 9 and ω 10 occur at point A, the two-order frequency trajectory remains consistent until the mode transition at point B separates rapidly; after that, these two phenomena occur simultaneously at point C. ω 7 and ω 8 change from the coincidence of trajectories before point C to the rapid separation, and then keep consistent with ω 9 , ω 6 and ω 7 , which modify at a single frequency to point D and intersect their trajectories; until the modal transition phenomenon occurs at point E, the two frequency trajectories quickly separate, and the same phenomenon as point C occurs again at point F. Before point F, the trajectories of ω 4 and ω 5 are the same. After point F, the trajectories of ω5 and ω4 are separated due to the modal transition and trajectory intersection, and ω 5 and ω 6 keep the same trajectory. This phenomenon is consistent with the torsional vibration, planetary wheel vibration and coupled vibration modes in the system. Trajectories of the inertia with natural J c , J s   ,   J r 2   , J a j   and   J b j frequency are similar to that of J r 1 , as shown in Figure 13, Figure 14, Figure 15, Figure 16 and Figure 17.

3.3. Time-Frequency Characteristic Analysis of System Coupling Response

Simulation condition: the rated speed is 2300   r / min , the mean torque is 2000   N . m and the fluctuation frequency is six times the rotation frequency of the crankshaft.
As seen in Figure 18a,b, the planet carrier and the gear ring are completely overlapped in frequency; namely, their vibration is caused by the same excitation and their vibration frequency coincides completely. It can be seen that the low-frequency components of the vibration are mainly composed of engine rotation frequency f e , modulation frequency f c , f s f c , the high-frequency components of the vibration are mainly composed of the planetary gear rotary frequency f m and the engine rotation frequency f m ± f e . It is also obvious that there are modulation components f m ± 6 f e induced by the fluctuation frequency 6 f e and the meshing frequency f m .
It can be seen from Figure 19 that the main peak frequencies include the engine rotating frequency f e , meshing frequency f m and its frequency doubling n f m . In addition, there are sidebands with planetary carrier rotating frequency f e , sun gear rotating frequency f s , planetary wheel rotating frequencies f a and f b , modulation frequencies f e , 6 f e , f b f e and n f m ± f c ,   n f m ± f s , n f m ± f a (n = 1, 2, 3…). The largest amplitude is the meshing frequency and double frequency of the system, while the smallest amplitude is the rotational frequency of the system and modulation frequency with each meshing frequency. In addition, there are special frequencies, which are mainly caused by nonlinear factors. The dominant frequency of meshing force is the same from Figure 20, the maximum amplitude is F m r b and the smallest amplitude is F m s a .

4. Verification of Model Accuracy

The virtual prototype model of the composite planetary gear system is established by combining CATIA and ADAMS for dynamic virtual simulation and vibration response analysis, the response results are compared and verified with the dynamic model.

4.1. Comparison of Natural Frequency

By comparing the mean value of the stabilized rotational speed with the calculated theoretical rotational speed, as shown in Table 2, it can be seen that the simulation errors of the rotational speed of the main components are all within 0.2%. It can be considered that the established virtual prototype model of the composite planetary wheel system has good accuracy and can be used for comparative analysis of numerical analysis results.
The undamped natural frequency of the composite planetary gear system simulation model is obtained by modal analysis and the results are compared with the numerical results as shown in Table 3. Due to some assumptions of lumped mass method equivalence and inaccurate measurement of parameters such as inertia and stiffness being adopted in dynamic modeling, there is a certain error between simulation results and numerical calculation values; however, the overall error is relatively small.

4.2. Comparison of Meshing Force

Simulation condition: the rated speed is 2300   r / min , the mean value is 2000   N . m and the fluctuation frequency is six times of the crankshaft frequency.
The numerical calculation values of meshing forces   F m r a , F m r b   and   F m a b are not completely consistent with the simulation results as shown in Table 4; the errors of all meshing forces are within 5%.
It can be clearly seen from Figure 21, Figure 22 and Figure 23 that, in each meshing force spectrum, the vibration response of wave frequency 6 f e is dominant, followed by the system meshing frequency and double frequency 2 f m and   3 f m , which are the same as that of numerical calculation. In the spectrum of the meshing force of the simulation results, there is a small value of the edge frequency band, mainly because the impact function is used in ADAMS to simulate the contact force of the gear meshing.
The comparative analysis of the virtual simulation and the numerical calculation results show that the coupling vibration response of the composite planetary gear system has a good correspondence in each vibration frequency component, which can form a good mutual verification with the established dynamic model and the numerical response results.

5. Conclusions

Taking time-varying stiffness, dynamic clearance, comprehensive error and engine harmonic excitation induced by the time-varying phase angle of the planetary gear into consideration, a time-varying dynamic model for composite planetary gear system is established in this paper. From this work, we can conclude that:
(1) The vibration of the composite planetary gear system has three modes: planetary gear vibration, global vibration and coupling vibration;
(2) The frequency trajectories of each order increase with the change in stiffness, although some frequency trajectories will have the phenomenon of modal transition and trajectory intersection. With the increase in inertia, the natural frequency of the system shows a decreasing trend and most orders change slightly. In addition, there are complex modal transitions and trajectory intersections in the natural frequency change trajectory, which are consistent with the torsional vibration, planetary wheel vibration and coupled vibration modes in the system;
(3) In bending–torsional direction, the vibration frequency components of different parts in the planetary gear are the same; however, their amplitudes are different. The low frequency is induced by the single frequency fluctuation torque converter input, which is independent of the system state. The high frequency is composed of the meshing frequency f m , six times the fluctuation frequency of engine rotating frequency 6 f e and its modulation frequencies f m ± 6 f e . The rotation frequency and harmonic frequency of the engine have a great influence on the vibration response of the system.;
(4) The main frequencies in the meshing force spectrum include the engine rotation frequency and the meshing frequency and its frequency multiplication. In addition, there are sideband frequencies, such as the rotation frequencies of planetary carrier, sun gear and planetary gear as well as their modulation frequency. The larger amplitudes are the meshing frequency and its frequency doubling, while the smaller amplitudes are the rotation frequencies and the modulation frequency of meshing frequency; there are also some special frequencies, which are mainly caused by the nonlinear factors of the system.

Author Contributions

Conceptualization, T.C. and C.Z.; validation, Y.L.; writing—original draft preparation, T.C. and C.Z.; writing—review and editing, Y.L. and Y.C.; supervision, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ambarisha, V.K.; Parker, R.G. Nonlinear dynamics of planetary gears using analytical and finite element models. J. Sound Vib. 2007, 302, 577–595. [Google Scholar] [CrossRef]
  2. Sun, T.; Hu, H. Nonlinear dynamics of planetary gear transmission by harmonic balance method based on DFT. Chin. J. Mech. Eng. 2002, 38, 58–61. [Google Scholar] [CrossRef]
  3. Yi, G.; Parker, G.R. Dynamic modeling and analysis of a spur planetary gear involving tooth wedging and bearing clearance nonlinearity. Eur. J. Mech.-A/Solids 2010, 29, 1022–1033. [Google Scholar]
  4. Guo, Y.; Keller, J.; Parker, R.G. Nonlinear dynamics and stability of wind turbine planetary gear sets under gravity effects. Eur. J. Mech.-A/Solids 2014, 47, 45–57. [Google Scholar] [CrossRef]
  5. Woohyung Kim, J.Y.L.J. Dynamic analysis for a planetary gear with time-varying pressure angles and contact ratios. J. Sound Vib. 2012, 331, 883–901. [Google Scholar]
  6. Zhou, J.; Sun, W.; Cao, L.; Wen, G. Vibro-Impact characteristics of planetary gear transmission. J. Xi’an Jiaotong Univ. 2016, 50, 16–21. [Google Scholar]
  7. Liu, H.; Cai, Z.; Xiang, C. Frequency coupling and dynamic characteristics of nonlinear meshing force for two-stage planetary gears. J. Vib. Shock 2015, 34, 13–23. [Google Scholar]
  8. Liu, C.; Qin, D.; Lim, T.C.; Liao, Y. Dynamic characteristics of the herringbone planetary gear set during the variable speed process. J. Sound Vib. 2014, 333, 6498–6515. [Google Scholar] [CrossRef]
  9. Lin, J.; Parker, R.G. Analytical Characterization of the Unique Properties of Planetary Gear Free Vibration. J. Vib. Acoust. 1999, 121, 316–321. [Google Scholar] [CrossRef]
  10. Kahraman, A. Natural Modes of Planetary Gear Trains. J. Sound Vib. 1994, 173, 125–130. [Google Scholar] [CrossRef]
  11. Peng, D.; Smith, W.A.; Randall, R.B.; Peng, Z. Use of mesh phasing to locate faulty planet gears. Mech. Syst. Signal Process. 2019, 116, 12–24. [Google Scholar] [CrossRef]
  12. Qiu, X.; Han, Q.; Chu, F. Review on Dynamic Analysis of Wind Turbine Geared Transmission Systems. J. Mech. Eng. 2014, 50, 23–36. [Google Scholar] [CrossRef]
  13. Xiang, L.; Gao, N.; Tang, L.; Guo, P. Nonlinear dynamic features of wind turbine’s gear systems subjected to internal and external excitations. Zhendong Yu Chongji/J. Vib. Shock 2018, 37, 126–132. [Google Scholar]
  14. Bozca, M. Kinematic analysis of a semi-rigid-elastic multi-body torsional vibration model of gears system. In Proceedings of the 12th International Carpathian Control Conference (ICCC), Velké Karlovice, Czech Republic, 25–28 May 2011; pp. 37–41. [Google Scholar]
  15. Zhang, J.; Liu, X.; Jiao, Y.; Song, Y. Vibration Analysis of Planetary Gear Trains Based on a Discrete-continuum Dynamic Model. J. Mech. Eng. 2014, 50, 104–112. [Google Scholar] [CrossRef]
  16. Qin, D.; Long, W.; Yang, J.; Zhou, H. Dynamic characteristics of wind turbine transmission system under varying wind speed and operation control conditions. J. Mech. Eng. 2012, 48, 1–8. [Google Scholar] [CrossRef]
  17. Dou, Z.; Li, Y.; Zeng, Z.; Du, M.; Yang, Y. Frequency coupling and coupling resonance of a composite planetary transmission system under complex excitations. J. Vib. Shock 2019, 3, 16–23. [Google Scholar]
  18. Wu, S.; Ren, H.; Zhu, E.; Wang, X.; Xu, Q.; Qian, B. Research advances for dynamics of planetary gear trains. Eng. J. Wuhan Univ. 2010, 43, 398–403. [Google Scholar]
  19. Wang, J.; Qin, D.; Ding, Y. Dynamic behavior of wind turbine by a mixed flexible-rigid multi-body model. J. Syst. Des. Dyn. 2009, 3, 403–419. [Google Scholar] [CrossRef] [Green Version]
  20. Li, G.; Li, F.; Zhu, Z. Simulation of dynamic characteristics of compound planetary gear trains. J. Cent. South Univ. (Sci. Technol.) 2015, 46, 1635–1640. [Google Scholar]
  21. Zhu, C.; Xu, X.; Lim, T.C.; Du, X.; Liu, M. Effect of flexible pin on the dynamic behaviors of wind turbine planetary gear drives. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2013, 227, 74–86. [Google Scholar] [CrossRef]
  22. Inalpolat, M.; Kahraman, A. A dynamic model to predict modulation sidebands of a planetary gear set having manufacturing errors. J. Sound Vib. 2010, 329, 371–393. [Google Scholar] [CrossRef]
  23. Chaari, F.; Fakhfakh, T.; Hbaieb, R.; Louati, J.; Haddar, M. Influence of manufacturing errors on the dynamic behavior of planetary gears. Int. J. Adv. Manuf. Technol. 2006, 27, 738–746. [Google Scholar] [CrossRef]
  24. Zhang, C.; Wei, J.; Wang, F.; Hou, S. Dynamic model and load sharing performance of planetary gear system with journal bearing. Mech. Mach. Theory 2020, 151, 1–22. [Google Scholar] [CrossRef]
  25. Gu, X.; Velex, P. A dynamic model to study the influence of plant position errors in planetary gears. J. Sound Vib. 2012, 331, 4554–4574. [Google Scholar] [CrossRef]
  26. Li, T.; Zhu, R.; Bao, H.; Xiang, C.; Liu, H. Nonlinear torsional vibration modeling and bifurcation characteristic study of a planetary gear train. J. Mech. Eng. 2011, 47, 76–83. [Google Scholar] [CrossRef]
  27. Ryali, L.; Talbot, D. A dynamic load distribution model of planetary gear sets. Mech. Mach. Theory 2020, 158, 104229. [Google Scholar] [CrossRef]
  28. Mo, S.; Zhang, T.; Jin, G.; Cao, X.; Gao, H. Analytical investigation on load sharing characteristics of herringbone planetary gear train with flexible support and floating sun gear. Mech. Mach. Theory 2020, 144, 1–27. [Google Scholar]
  29. Sanchez-Espiga, J.; Fernandez-Del-Rincon, A.; Iglesias, M.; Viadero, F. Planetary gear transmissions load sharing measurement from tooth root strains: Numerical evaluation of mesh phasing influence. Mech. Mach. Theory 2021, 163, 104370. [Google Scholar] [CrossRef]
  30. Iglesias, M.; Fernandez, D.; De-Juan, A.; Garcia, P.; Diez-Ibarbia, A.; Viadero, F. Planetary transmission load sharing: Manufacturing errors and system configuration study. Mech. Mach. Theory 2017, 111, 21–38. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic of the compound planetary transmission system.
Figure 1. Schematic of the compound planetary transmission system.
Machines 10 00031 g001
Figure 2. Bending–torsional coupling mechanical model of composite planetary gear system.
Figure 2. Bending–torsional coupling mechanical model of composite planetary gear system.
Machines 10 00031 g002
Figure 3. Torsional meshing deformation relationship between the planetary gear a j and the planetary gear b j .
Figure 3. Torsional meshing deformation relationship between the planetary gear a j and the planetary gear b j .
Machines 10 00031 g003
Figure 4. Translational meshing deformation relationship between the planetary gears a j and b j : (a) Projection relation of aj translation displacement; (b) Projection relation of bj translation displacement.
Figure 4. Translational meshing deformation relationship between the planetary gears a j and b j : (a) Projection relation of aj translation displacement; (b) Projection relation of bj translation displacement.
Machines 10 00031 g004
Figure 5. Engine harmonic excitation: (a) The excitation torque of engine single cylinder; (b) Harmonic relation of engine output torque.
Figure 5. Engine harmonic excitation: (a) The excitation torque of engine single cylinder; (b) Harmonic relation of engine output torque.
Machines 10 00031 g005
Figure 6. The vibration modes of each order of the composite planetary gear system.
Figure 6. The vibration modes of each order of the composite planetary gear system.
Machines 10 00031 g006
Figure 7. The vibration modes of each order of the composite planetary gear system: (a) Planetary wheel vibration; (b) Global vibration; (c) Coupled vibration.
Figure 7. The vibration modes of each order of the composite planetary gear system: (a) Planetary wheel vibration; (b) Global vibration; (c) Coupled vibration.
Machines 10 00031 g007
Figure 8. The trajectory of the natural frequency with k m r a j .
Figure 8. The trajectory of the natural frequency with k m r a j .
Machines 10 00031 g008
Figure 9. The trajectory of the natural frequency with k m s a j .
Figure 9. The trajectory of the natural frequency with k m s a j .
Machines 10 00031 g009
Figure 10. The trajectory of the natural frequency with k m r b j .
Figure 10. The trajectory of the natural frequency with k m r b j .
Machines 10 00031 g010
Figure 11. The trajectory of the natural frequency with k m a b j .
Figure 11. The trajectory of the natural frequency with k m a b j .
Machines 10 00031 g011
Figure 12. The trajectory of the natural frequency with J r 1 .
Figure 12. The trajectory of the natural frequency with J r 1 .
Machines 10 00031 g012
Figure 13. The trajectory of the natural frequency with J c .
Figure 13. The trajectory of the natural frequency with J c .
Machines 10 00031 g013
Figure 14. The trajectory of the natural frequency with J s .
Figure 14. The trajectory of the natural frequency with J s .
Machines 10 00031 g014
Figure 15. The trajectory of the natural frequency with J r 2 .
Figure 15. The trajectory of the natural frequency with J r 2 .
Machines 10 00031 g015
Figure 16. The trajectory of the natural frequency with J a j .
Figure 16. The trajectory of the natural frequency with J a j .
Machines 10 00031 g016
Figure 17. The trajectory of the natural frequency with J b j .
Figure 17. The trajectory of the natural frequency with J b j .
Machines 10 00031 g017
Figure 18. Frequency of vibration displacement of planetary carrier and small gear ring: (a) Spectrum of translational displacement; (b) Spectrum of torsional displacement.
Figure 18. Frequency of vibration displacement of planetary carrier and small gear ring: (a) Spectrum of translational displacement; (b) Spectrum of torsional displacement.
Machines 10 00031 g018
Figure 19. The spectrum of each meshing force.
Figure 19. The spectrum of each meshing force.
Machines 10 00031 g019
Figure 20. Time-frequency of each meshing force: (a) The meshing force F m r b between the large ring gear and planetary gear b; (b) The meshing force F m s a between the sun gear and the planet gear; (c) The meshing force F m r a between the small gear ring and the planetary gear a; (d) The meshing force F m a b between planetary gears.
Figure 20. Time-frequency of each meshing force: (a) The meshing force F m r b between the large ring gear and planetary gear b; (b) The meshing force F m s a between the sun gear and the planet gear; (c) The meshing force F m r a between the small gear ring and the planetary gear a; (d) The meshing force F m a b between planetary gears.
Machines 10 00031 g020
Figure 21. Frequency spectrum and time–frequency spectrum of meshing force F m a b .
Figure 21. Frequency spectrum and time–frequency spectrum of meshing force F m a b .
Machines 10 00031 g021
Figure 22. Frequency spectrum and time–frequency spectrum of meshing force F m r a .
Figure 22. Frequency spectrum and time–frequency spectrum of meshing force F m r a .
Machines 10 00031 g022
Figure 23. Frequency spectrum and time–frequency spectrum of meshing force F m r b .
Figure 23. Frequency spectrum and time–frequency spectrum of meshing force F m r b .
Machines 10 00031 g023
Table 1. Main parameters of the transmission system.
Table 1. Main parameters of the transmission system.
Number of TeethMoment of Inertia (kg*m2)Mass (kg)
The small ring r1770.254611.275
The big ring r2820.497714.434
sun gear340.01406.323
planet gear a210.00282.163
planet gear b220.00171.094
Table 2. Comparison between simulation value and theoretical values of rotational speed.
Table 2. Comparison between simulation value and theoretical values of rotational speed.
Ring Gear r1 (Input)Sun Gear sPlanet Gear ajPlanet Gear bjPlanet Carrier c (Output)
Theoretical values (r/min)23001572.55463.13235.41113.8
Simulation values (r/min)23001574.25452.73234.01113.9
Errors (%)00.110.190.0450.0048
Table 3. Comparison of natural frequency simulation value and theoretical value.
Table 3. Comparison of natural frequency simulation value and theoretical value.
Degree123456
Theoretical values(Hz)01339.12624.23281.35075.07420.5
Simulation values (Hz)01324.92508.03206.65092.76970.2
Errors(%)01.064.432.280.356.07
Table 4. Comparison between simulation and theoretical values of meshing forces.
Table 4. Comparison between simulation and theoretical values of meshing forces.
Mean Meshing Force F m r a F m r b F m a b
Theoretical values (N)4360.46204.66177.3
Simulation values (N)4391.95916.46295.9
Errors (%)0.724.641.92
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cui, T.; Li, Y.; Zan, C.; Chen, Y. Dynamic Modeling and Analysis of Nonlinear Compound Planetary System. Machines 2022, 10, 31. https://doi.org/10.3390/machines10010031

AMA Style

Cui T, Li Y, Zan C, Chen Y. Dynamic Modeling and Analysis of Nonlinear Compound Planetary System. Machines. 2022; 10(1):31. https://doi.org/10.3390/machines10010031

Chicago/Turabian Style

Cui, Tingqiong, Yinong Li, Chenglin Zan, and Yuanchang Chen. 2022. "Dynamic Modeling and Analysis of Nonlinear Compound Planetary System" Machines 10, no. 1: 31. https://doi.org/10.3390/machines10010031

APA Style

Cui, T., Li, Y., Zan, C., & Chen, Y. (2022). Dynamic Modeling and Analysis of Nonlinear Compound Planetary System. Machines, 10(1), 31. https://doi.org/10.3390/machines10010031

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop