Next Article in Journal
Moving beyond the Content: 3D Scanning and Post-Processing Analysis of the Cuneiform Tablets of the Turin Collection
Previous Article in Journal
The Importance of Micronutrient Adequacy in Obesity and the Potential of Microbiota Interventions to Support It
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of the Dynamic Characteristics of Coaxial Counter-Rotating Planetary Transmission System

1
School of Mechanics Engineering, Harbin Institute of Technology, Harbin 150001, China
2
School of Mechanics and Mathematics, Belarusian State University, 220030 Minsk, Belarus
3
School of Mechanical Engineering, Guangxi University, Nanning 530004, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(11), 4491; https://doi.org/10.3390/app14114491
Submission received: 24 April 2024 / Revised: 16 May 2024 / Accepted: 22 May 2024 / Published: 24 May 2024

Abstract

:
This paper presents a coaxial counter-rotating planetary transmission system. The transmission system under study is a two-stage planetary gear train (PGT) comprising a fixed-axes PGT and a differential PGT. A dynamic model of the transmission system is established, considering both the excitations caused by the time-varying mesh stiffness (TMS) and the transmission errors, respectively. The Runge–Kutta algorithm is used to calculate and analyze the dynamic characteristics of the system. This includes studying dynamic meshing forces, planet gear displacements, and load-sharing coefficients (LSCs) under both internal and external excitations, as well as different input torques. The results indicate that when considering external excitations, the variations in the meshing force curves become more pronounced. The radial displacements of the planet gears in the differential PGT are greater than that in the fixed-axes PGT. With increasing input torque, the average displacements of the planet gears in all directions tend to increase. The differential PGT, transmitting a higher power, demonstrates a better load-sharing performance compared to the fixed-axes PGT.

1. Introduction

The coaxial counter-rotating planetary transmission system under investigation in this study is a two-stage PGT consisting of a fixed-axes PGT and a differential PGT. Planetary gear transmissions are favored for their high torque density, smooth output, compact design, and high efficiency. The coaxial counter-rotating planetary transmission system brings numerous advantages to the power transmission of helicopter rotor shafts and ship propellers. This type of transmission system enhances transmission efficiency, reduces energy loss, thereby extending equipment operating time, particularly suitable for equipment with long-duration operations, such as with helicopters and ships. Furthermore, the coaxial counter-rotating planetary transmission system, by concentrating power transmission along a single axis, saves space, allowing for more compact equipment design, which is adaptable to the limited space environment in aerospace and maritime sectors. The transmission system provides precise and smooth power transmission, improving the maneuverability and stability of helicopters and ships, particularly evident when operating in complex environments. Due to the typically fewer mechanical connections in coaxial counter-rotating planetary transmission systems, their maintenance costs are relatively low, contributing to reduced operational costs and prolonged equipment lifespan.
The field of PGTs has witnessed substantial enrichment through the noteworthy contributions of numerous scholars dedicated to exploring the dynamic characteristics of these systems. Bahk et al. [1] introduced a computational model to analyze the excitation of teeth profiles following modification. They conducted additional investigations to examine the effects of profile modification on the dynamic response of the system through the utilization of perturbation analysis. Shen et al. [2] introduced a torsional dynamic model specifically designed for planet gears. This model incorporates multiple elements, including gear wear, TMS, unloaded static transmission errors, and gear backlash. By incorporating these elements into the dynamic model, they aimed to provide a more comprehensive understanding of the behavior of planet gears under torsional loads. Xu et al. [3] employed a discretization technique to transform the continuous gear profile into a set of discrete points, which were subsequently utilized for planet gear modeling. Gu et al. [4,5] investigated the trajectories of components and the distribution of loads in planet gears, taking into account eccentricity errors and positional errors. Guo et al. [6] employed real-time excitations in harmonic balance to investigate the dynamics of wind turbine planetary gear sets. Huangfu et al. [7] proposed a continuous-discrete model for helical gears. Based on the dynamics model proposed, they investigated the transmission effects of planet gears, the influence of bolt constraints, and gear vibration. Han et al. [8] presented a three-dimensional load tooth contact analysis (3DLTCA) model considering changes in the center distance and root crack discontinuity of planet gears, and they calculated the dynamic response under the combined effects of assembly errors and crack discontinuities. Wang et al. [9], based on the translational–rotational dynamic model of planetary gear systems, studied the effects of multiple excitations on the vibration response of various planet gear components. Liu et al. [10] coupled the concentrated mass dynamic submodel of sun gears and planet gears with the reduced-order finite element submodel and used this model to analyze the coupled dynamics of flexible planet gears under high-speed and variable-speed conditions.
In the field of multistage compound planet gears and coaxial counter-rotating gears, Liu et al. [11] proposed a dynamic model for the Dual-Planetary Gear System (DPGS) that integrates the dynamic characteristics of the planet bearing rollers and retainers. Fyler et al. [12] proposed a two-dimensional steady-state discrete dynamic model for a dual-stag planetary gear set. They obtained the forced vibration response to gear mesh excitations by applying the modal summation technique. Guo et al. [13] investigated the vibration and acoustic characteristics of a power-split hybrid transmission utilizing a compound planetary gear set through numerical processes and refined dynamics modeling. Tung et al. [14] developed a time-varying reduced-order numerical model for a compound reducer, enabling the prediction of its meshing frequencies and harmonics. Zhang et al. [15] established a translational–torsional concentrated mass dynamic model for a compound planetary gear set, calculating eccentric load coefficients under the influence of excitations such as tooth wear accumulation and transmission errors. Yang et al. [16] proposed a multinode dynamic model for a marine compound gear transmission system, investigating its modal characteristics and dynamic responses. Li et al. [17] developed a dynamic model for a two-stage parallel-axis gear system with variable mesh stiffness, revealing the variations in system vibration under the interaction of coaxial gear ratios and meshing phases. Zhang et al. [18] employed a concentrated mass method to establish a translational–rotational dynamic model for a planetary gear system, analyzing the natural frequencies and mode characteristics of a coaxial counter-rotating helicopter main gearbox. Ryali et al. [19] developed a hybrid planet dynamic load distribution model, employing finite element substructuring techniques to ensure computational efficiency and investigated the influence of planet carrier flexibility on the quasi-static and dynamic responses of planetary gear sets. Lai et al. [20] proposed a flexible–rigid coupling dynamics (FRCD) model for Ravigneaux planetary gear sets (RPGS) with unloaded floating sun gears and studied the influence of the sun gear position on planet dynamic responses, revealing the mechanism of floating ring gear vibration. Wang et al. [21] proposed a novel split-type structure for first- and second-stage gears of planetary reducers. Through rigid body dynamics theory, dynamic simulations, noise tests, and dynamic stability tests were conducted to determine the dynamic characteristics of the entire system. Cui et al. [22], focusing on compound planetary gears, established a bending–torsional coupled dynamic model of the system based on Lagrange’s equations, investigating the bending–torsional coupled meshing deformations and translations and torsions of various engagements. Zhang et al. [23] studied the inherent vibration modes of Ravigneaux compound planetary gear trains (PGT), proposing a two degrees of freedom simplified dynamic model of planet gears considering nonlinearities, such as time-varying mesh stiffness and backlash, and obtaining numerical solutions for the bifurcation and dynamic instability of two groups of planetary gear systems.
In summary, the dynamics models of PGTs are mainly divided into concentrated mass analysis models, finite element analysis models, and rigid–flexible coupling analysis models. In the concentrated mass analysis model, the masses of all components (sun gear, planet gears, ring gear, carrier) in the PGTs are concentrated at single points, and these concentrated mass points are interconnected by springs with stiffness and damping. The finite element model utilizes finite element software to conduct dynamic analysis of gear transmission systems to obtain contact stresses, bending stresses, etc., while the rigid–flexible coupling model is a modeling method that treats components with large deformations as flexible based on the concentrated mass method.
The aforementioned studies on the dynamics of planet gears have mostly focused on single-stage transmissions or simple two-stage serial planetary transmissions with a single input and single output. There has been relatively limited research on coaxial counter-rotating planetary transmission systems composed of fixed-axes PGT and differential PGT. Although finite element analysis models and rigid–flexible coupling analysis models can extract more dynamic information, there are issues such as low efficiency when solving complex systems. Therefore, for complex mechanisms like coaxial counter-rotating planetary transmission systems, employing concentrated mass analysis models can most quickly extract dynamic information that is closer to reality. When using finite element analysis models and rigid–flexible coupling analysis models for calculations, it often takes several hours or even days. However, using the model proposed in this paper for numerical calculations only takes a few minutes. Moreover, when the parameters of the transmission system change, simply replacing the parameters in the program is sufficient, avoiding the complex process of repetitive modeling and grid partitioning required by the finite element method. From the results, which are provided in Section 5 of this paper, through the comparison between ADAMS dynamic simulation and the results of this paper, also proves the correctness of the model.
Therefore, this paper investigates the coaxial counter-rotating planetary transmission system, establishing a dynamic model of the transmission system that considers time-varying mesh stiffness and transmission error excitations. The dynamic characteristics of the system are computed and analyzed using the Runge–Kutta algorithm. This includes studying dynamic meshing forces, planet gear displacements, and load-sharing coefficients (LSCs) under both internal and external excitations, as well as different input torques. The main original findings in this study can be directly employed as preliminary principles for the design of coaxial counter-rotating planetary transmission systems.

2. Composition of the Transmission System

The diagrammatic representation and model of the coaxial counter-rotating planetary transmission system are shown in Figure 1. The entire system is composed of the central components (sun gears—a and s; left ring gears—dL; right ring gears—dR; left ring gears—rL; right ring gears—rR; carriers—c and h; and planet gears—mi and pi), as well as the input shaft and output shaft. This transmission design incorporates a combination of fixed-axes and differential PGTs. The input power is divided into two segments within the system. A specific portion is designated for the fixed-axes PGT, flowing from the sun gear, a, to the ring gear, d. Concurrently, the remaining section is allocated for the differential PGT, starting from the sun gear s and further distributed to the carrier, h, and ring gear, r. The transmission system is represented in the form of the Wolf–Arnaudov symbol [24], as shown in Figure 2. In the figure, t and t denote the transmission ratio of the fixed-axes PGT and the differential PGT, respectively, and t = t = t.

3. Dynamic Model of the Transmission System

3.1. Dynamic Model

Taking the example of the differential PGT, the dynamic model, illustrated in Figure 3, is established taking into account parameter excitations, such as TMS, dynamic transmission errors, and support conditions, during the gear meshing process. The fixed-axes PGT follows a similar model, with only different subscripts. Each component in this model has four degrees of freedom. The variables xj, yj, and zj (j = a, m, dL, dR, m1, …, mN, s, h, rL, rR, p1, …, pN) represent the translational displacements of the sun gear, carrier, ring gear, and each planet gear in the x, y, and z directions, respectively. The component torsional displacements around the z-axis are represented by θzj. The symbols k, c, and e in the figure stand for mesh stiffness, mesh damping, and transmission errors, respectively, with distinct components denoted by subscripts. The left and right gears are represented, respectively, by the subscripts L and R. Bearing stiffness is represented by subscripts with b. As an example, kbxs represents the support stiffness of the sun gear, s, in the x-direction. Similarly, kspiL represents the mesh stiffness of the left-side gear pair between the sun gear, s, and the planet gear, p. Likewise, the remaining symbols utilized in the model adhere to the identical notation.

3.2. Calculation of Component Relative Displacements

The generalized coordinate and mass matrix of each component are as follows:
X j = [ x j , y j , z j , θ z j ] T M j = d i a g [ m j , m j , m j , J j ] ,
where Mj represents the mass of the component, and Jj represents the moment of inertia with respect to the axis of rotation of the component.
For the purpose of modeling and solving the dynamic model, a fixed coordinate system, oxyz, is chosen for the fixed-axes PGT, while a follower coordinate system, ohxhyhzh, is selected for the differential PGT, with the carrier serving as the reference. The displacements and deformations of the components are projected onto meshing lines to obtain the relative displacements among the components.

3.2.1. Relative Displacements between Sun gear and Planet gears

As an illustration, the relationship between the i-th planet gear, p, and the sun gear, s, in the differential PGT, as shown in Figure 4. In the figure, ψpi represents the angular position of the planet gear.
Along the meshing line, the elastic compression deformation of sun gear, a, and planet gear, mi, can be expressed as positive for compression and negative for tension, as follows:
δ a m i j = ( x a sin ψ a m i + y a cos ψ a m i x m i sin α a m i y m i cos α a m i + r b a θ z a r b m θ z m i ) cos β b ± ( z a + z m i ) sin β b e a m i j ( t ) ,
where j can take on the values of either L or R; ψami represents the angle between the meshing line, ami, and the horizontal axis; αami represents the pressure angle; rba and rbm represent the base circle radius of the sun gear, a, and planet gear, m; βb represents the helix angle; eamij represents the transmission error between the sun gear, a, and planet gear, m; and the superscript sign (±) applies to the left-side gear, while the subscript sign (±) applies to the right-side gear.
The elastic compression deformation of the sun gear s and planet gear, pi, along the meshing line is as follows:
δ s p i j = ( x s sin ψ s p i + y s cos ψ s p i x p i sin α s p i y p i cos α s p i + r b s θ z s r b p θ z p i ) cos β b ± ( z s + z p i ) sin β b e s p i j ( t ) ,
where ψspi represents the angle between the meshing line, spi, and the horizontal axis; αspi represents the pressure angle; rbs and rbp represent the base circle radii of the sun gear, s, and planet gear, p; and espij represents the transmission error between the sun gear, s, and planet gear, p.

3.2.2. Relative Displacements between the Ring gear and Planet Gears

As an illustration, the relationship between the i-th planet gear, p, and the ring gear, r, in the differential PGT is shown in Figure 5.
The elastic compression deformation of the ring gear, d, and planet gear, mi, along the meshing line is as follows:
δ d m i j = ( x d j sin ψ d m i + y d j cos ψ d m i + x m i sin α d m i y m i cos α d m i + r b m θ z m i r b d θ z d j ) cos β b ± ( z d j z m i ) sin β b e d m i j ( t ) ,
where ψdmi represents the angle between the meshing line, dmi, and the horizontal axis; αdmi represents the pressure angle; rbd represents the base radius of the ring gear, d; and edmij represents the transmission error between the ring gear, d, and planet gear, m.
The elastic compression deformation of the ring gear r and planet gear pi along the meshing line is as follows:
δ r p i j = ( x r j sin ψ r p i + y r j cos ψ r p i + x p i sin α r p i y p i cos α r p i + r b p θ z p i j r b r θ z r j ) cos β b ± ( z r j z p i ) sin β b e r p i j ( t ) ,
where ψrpi represents the angle between the meshing line, rpi, and the horizontal axis; αrpi represents the pressure angle; rbr represents the base circle radius of the ring gear, r; and erpij represents the transmission error between the ring gear, r, and planet gear, p.

3.2.3. Relative Displacements between the Carrier and Planet Gears

The relationship between the planet gear, p, and the carrier, h, in the differential PGT is shown in Figure 6.
Neglecting the clearance of the bearings, the supporting bearings are approximated by a simplified model using stiffness and damping. Therefore, the relative displacements between the carrier, h, and the planet gear, pi, projected onto the carrier’s coordinate axes xh, yh, and zh are as follows:
δ x h p i = x h x p i cos ψ p i + y p i sin ψ p i r h θ z h sin ψ p i δ y h p i = y h x p i sin ψ p i y p i cos ψ p i + r h θ z h cos ψ p i δ z h p i = z h z p i θ z h p i = θ z h + θ z p i ,
where ψpi represents the angle between the xpi axis of the planet gear p and the xh axis of the carrier.
The relative displacements between the carrier, h, and the planet gear, pi, projected onto the planet gear’s coordinate axes xpi, ypi, and zpi are as follows:
δ x p i h = x p i x h cos ψ p i y h sin ψ p i δ y p i h = y p i + x h sin ψ p i y h cos ψ p i r h θ z h δ z p i h = z p i z h θ z p i h = θ z p i + θ z h ,

3.3. Derivation of the Motion Differential Equations

In the dynamic model, a load torque, Tout1, is applied to the output shaft of the ring gear, rR; a load torque, Tout2, is applied to the output shaft of the carrier, h; and an input torque, Tin, is applied to the sun gear, a. By analyzing the forces acting on each component and employing the Lagrange equations, the motion differential equations for each component, considering the translational motion along the xj, yj, and zj axes and the torsional motion around the zj axis (totaling four degrees of freedom), can be established.

3.3.1. Motion Differential Equation for the Sun Gear, a

The bearing support force acting on the sun gear, a, is as follows:
F b x a = c b x a x ˙ a + k b x a x a F b y a = c b y a y ˙ a + k b y a y a F b z a = c b z a z ˙ a + k b z a z a M b θ z a = c b θ z a θ ˙ z a + k b θ z a θ z a ,
where kxas, kyas, kzas, and kθzas represent the support stiffnesses of the bearing in the respective directions, and cbxa, cbya, cbza, and cbθza represent the support dampings of the bearing.
The support force of the shaft acting on the sun gear, a, is as follows:
F x a s = c x a s ( x ˙ a x ˙ s ) + k x a s ( x a x s ) F y a s = c y a s ( y ˙ a y ˙ s ) + k y a s ( y a y s ) F z a s = c z a s ( z ˙ a z ˙ s ) + k z a s ( z a z s ) M θ z a s = c θ z a s ( θ ˙ z a θ ˙ z s ) + k θ z a s ( θ z a θ z s ) ,
where kxas, kyas, kza, and kθzas represent the support stiffnesses of the shaft in the respective directions, and cxas, cyas, czas, and cθzas represent the support dampings.
The force acting on the sun gear, a, as a result of its mesh with the planet gear, mi, is as follows:
F a m i L = k a m i L δ a m i L + c a m i L δ ˙ a m i L F a m i R = k a m i R δ a m i R + c a m i R δ ˙ a m i R ,
where kamiL and kamiR represent the mesh stiffnesses of the left and right gear pairs between the sun gears, a, and the planet gears, m, respectively; and camiL and camiR represent the mesh dampings.
The motion differential equation for the sun gear, a, is as follows:
m a x ¨ a + i = 1 N ( F a m i L + F a m i R ) cos β b sin ψ a m i + F b x a + F x a s = 0 m a y ¨ a + i = 1 N ( F a m i L + F a m i R ) cos β b cos ψ a m i + F b y a + F y a s = 0 m a z ¨ a i = 1 N ( F a m i L F a m i R ) sin β b + F b z a + F z a s = 0 J a θ ¨ z a + i = 1 N ( F a m i L + F a m i R ) r b a cos β b + M b θ z a + M θ z a s = T i n ,

3.3.2. Motion Differential Equation for the Planet Gear, mi

The support force of the bearing acting on the planet gear, mi, is as follows:
F b x m i = k b x m x m i + c b x m x ˙ m i F b y m i = k b y m y m i + c b y m y ˙ m i F b z m i = k b z m z m i + c b z m z ˙ m i M b θ z m i = k b θ z m θ z m i + c b θ z m θ ˙ z m i ,
where kbxm, kbym, kbzm, and kbθzm represent the support stiffnesses of the bearing in the respective directions, and cbxm, cbym, cbzm, and cbθzm represent the support dampings of the bearing.
The meshing force acting on the planet gear, mi, as a result of its mesh with the ring gear, d, is as follows:
F d m i L = k d m i L δ d m i L + c d m i L δ ˙ d m i L F d m i R = k d m i R δ d m i R + c d m i R δ ˙ d m i R ,
where kdmiL and kdmiR represent the mesh stiffnesses of the gear pairs between the left and right ring gears, d, and planet gears, m; and cdmiL and cdmiR represent the mesh dampings.
The motion differential equation for the planet gear, mi, is as follows:
m m x ¨ m i ( F a m i L + F a m i R ) sin α a m i ( F d m i L + F d m i R ) sin α d m i cos β b + F b x m i = 0 m m y ¨ m i ( F a m i L + F a m i R ) cos α a m i + ( F d m i L + F d m i R ) cos α d m i cos β b + F b y m i = 0 m m z ¨ m i + ( F a m i L F a m i R F d m i L + F d m i R ) sin β b + F b z m i = 0 J m θ ¨ z m i + ( F d m i L + F d m i R F a m i L F a m i R ) r b m cos β b + M b θ z m i = 0 ,

3.3.3. Motion Differential Equation for the Ring Gear, d

The support force of the bearing acting on the left-side ring gear, dL, is as follows:
F b x d L = c b x d L x ˙ d L + k b x d L x d L F b y d L = c b y d L y ˙ d L + k b y d L y d L F b z d L = c b z d L z ˙ d L + k b z d L z d L M b θ z d L = c b θ z d L θ ˙ z d L + k b θ z d L θ z d L ,
where kbxdL, kbydL, kbzdL, and kbθzdL represent the support stiffnesses of the bearing in the respective directions, and cbxdL, cbydL, cbzdL, and cbθzdL represent the support dampings of the bearing.
The support force of the shaft acting on the left internal gear, dL, is as follows:
F x d L d R = c x d L d R ( x ˙ d L x ˙ d R ) + k x d L d R ( x d L x d R ) F y d L d R = c y d L d R ( y ˙ d L y ˙ d R ) + k y d L d R ( y d L y d R ) F z d L d R = c z d L d R ( z ˙ d L z ˙ d R ) + k z d L d R ( z d L z d R ) M θ z d L d R = c θ z d L d R ( θ ˙ z d L θ ˙ z d R ) + k θ z d L d R ( θ z d L θ z d R ) ,
where kxdLdR, kydLdR, kzdLdR, and kθzdLdR represent the support stiffnesses of the shaft in the respective directions, and cxdLdR, cydLdR, czdLdR, and cθzdLdR represent the support dampings.
The motion equation for the left-side ring gear, dL, is as follows:
m d x ¨ d L i = 1 N F d m i L sin ψ d m i cos β b + F b x d L + F x d L d R = 0 m d y ¨ d L + i = 1 N F d m i L cos ψ d m i cos β b + F b y d L + F y d L d R = 0 m d z ¨ d L + i = 1 N F d m i L sin β b + F b z d L + F z d L d R = 0 J d θ ¨ z d L i = 1 N F d m i L r b d cos β b + M b θ z d L + M θ z d L d R = 0 ,
The support force of the shaft acting on the right-side ring gear, dR, is as follows:
F x d L d R r L = c x d L d R ( x ˙ d R x ˙ d L ) + k x d L d R ( x d R x d L ) + c x d R r L ( x ˙ d R x ˙ r L ) + k x d R r L ( x d R x r L ) F y d L d R r L = c y d L d R ( y ˙ d R y ˙ d L ) + k y d L d R ( y d R y d L ) + c y d R r L ( x ˙ d R x ˙ r L ) + k y d R r L ( y d R y r L ) F z d L d R r L = c z d L d R ( z ˙ d R z ˙ d L ) + k z d L d R ( z d R z d L ) + c z d R r L ( z ˙ d R z ˙ r L ) + k z d R r L ( z d R z r L ) M θ z d L d R r L = c θ z d L d R ( θ ˙ z d R θ ˙ z d L ) + k θ z d L d R ( θ z d R θ z d L ) + c θ z d R r L ( θ ˙ z d R θ ˙ z r L ) + k θ z d R r L ( θ z d R θ z r L ) ,
where kxdRrL, kydRrL, kzdRrL, and kθzdRrL represent the support stiffnesses of the shaft in the respective directions, and cxdRrL, cydRrL, czdRrL, and cθzdRrL represent the support dampings.
The motion equation for the right-side ring gear, dR, is as follows:
m d x ¨ d R i = 1 N F d m i R sin ψ d m i cos β b + F x d L d R r L = 0 m d y ¨ d R + i = 1 N F d m i R cos ψ d m i cos β b + F y d L d R r L = 0 m d z ¨ d R i = 1 N F d m i R sin β b + F z d L d R r L = 0 J d θ ¨ z d R i = 1 N F d m i R r b d cos β b + M θ z d L d R r L = 0 ,

3.3.4. Motion Differential Equation for the Sun Gear s

The support force of the bearing acting on the sun gear, s, is as follows:
F b x s = c b x s x ˙ s ω h y s + k b x s x s F b y s = c b y s y ˙ s + ω h x s + k b y s y s F b z s = c b z s z ˙ s + k b z s z s M b θ z s = c b θ z s θ ˙ z s + k b θ z s θ z s ,
where kbxs, kbys, kbzs, and kbθzs represent the support stiffnesses of the bearing in the respective directions; cbxs, cbys, cbzs, and cbθzs represent the support dampings of the bearing; and wh represents the angular velocity of the carrier h.
The support force of the shaft acting on the sun gear, s, is as follows:
F x s a = c x a s ( x ˙ s x ˙ a ) + k x a s ( x s x a ) F y s a = c y a s ( y ˙ s y ˙ a ) + k y a s ( y s y a ) F z s a = c z a s ( z ˙ s z ˙ a ) + k z a s ( z s z a ) M θ z s a = c θ z a s ( θ ˙ z s θ ˙ z a ) + k θ z a s ( θ z s θ z a ) ,
The force acting on the sun gear, s, as a result of its mesh with the planet gear, pi, is as follows:
F s p i L = k s p i L δ s p i L + c s p i L δ ˙ s p i L F s p i R = k s p i R δ s p i R + c s p i R δ ˙ s p i R ,
where kspiL and kspiR represent the mesh stiffnesses of the gear pairs between the left and right sun gears, s, and planet gears, p; and cspiL and cspiR represent the mesh dampings.
The motion differential equation for the sun gear, s, is as follows:
m s x ¨ s 2 ω h y ˙ s ω h 2 x s + i = 1 N ( F s p i L + F s p i R ) cos β b sin ψ s p i + F b x s + F x s a = 0 m s y ¨ s + 2 ω h x ˙ s ω h 2 y s + i = 1 N ( F s p i L + F s p i R ) cos β b cos ψ s p i + F b y s + F y s a = 0 m s z ¨ s i = 1 N ( F s p i L F s p i R ) sin β b + F b z s + F z s a = 0 J s θ ¨ z s + i = 1 N ( F s p i L + F s p i R ) r b s cos β b + M b θ z s + M θ z s a = 0 ,

3.3.5. Motion Differential Equation for the Planet Gear pi

The support force of bearing acting on the planet gear pi is as follows:
F b x p i = k b x p δ x p i h + c b x p δ ˙ x p i h F b y p i = k b y p δ y p i h + c b y p δ ˙ y p i h F b z p i = k b z p δ z p i h + c b z p δ ˙ z p i h M b θ z p i = k b θ z p θ z p i h + c b θ z p θ ˙ z p i h ,
where kbxp, kbyp, kbzp and kbθzp represent the support stiffness of the bearing in the respective directions, and cbxp, cbyp, cbzp and cbθzp represent the support damping of the bearing.
The force acting on the planet gear pi as a result of its mesh with the ring gear r is as follows:
F r p i L = k r p i L δ r p i L + c r p i L δ ˙ r p i L F r p i R = k r p i R δ r p i R + c r p i R δ ˙ r p i R ,
where krpiL and krpiR represent the mesh stiffnesses of the gear pairs between the left and right ring gears r and planet gears p, and crpiL and crpiR represent the mesh dampings.
The motion differential equation for the planet gear pi is as follows:
m p ( x ¨ p i 2 ω h y ˙ p i ω h 2 x p i ) ( F s p i L + F s p i R ) sin α s p i ( F r p i L + F r p i R ) sin α r p i cos β b + F b x p i = 0 m p ( y ¨ p i + 2 ω h x ˙ p i ω h 2 y p i ) ( F s p i L + F s p i R ) cos α s p i + ( F r p i L + F r p i R ) cos α r p i cos β b + F b y p i = 0 m p z ¨ p i + ( F s p i L F s p i R F r p i L + F r p i R ) sin β b + F b z p i = 0 J p θ ¨ z p i + ( F r p i L + F r p i R F s p i L F s p i R ) r b p cos β b + M b θ z p i = 0 ,

3.3.6. Motion Differential Equation for the Ring Gear r

The support force of bearing acting on the left-side ring gear rL is as follows:
F x d R r L r R = c x d R r L ( x ˙ r L x ˙ d R ) + k x d R r L ( x r L x d R ) + c x r L r R ( x ˙ r L x ˙ r R ) + k x r L r R ( x r L x r R ) F y d R r L r R = c y d R r L ( y ˙ r L y ˙ d R ) + k y d R r L ( y r L y d R ) + c y r L r R ( y ˙ r L y ˙ r R ) + k y r L r R ( y r L y r R ) F z d R r L r R = c z d R r L ( z ˙ r L z ˙ d R ) + k z d R r L ( z r L z d R ) + c z r L r R ( z ˙ r L z ˙ r R ) + k z r L r R ( z r L z r R ) M θ z d R r L r R = c θ z d R r L ( θ ˙ z r L θ ˙ z d R ) + k θ z d R r L ( θ z r L θ z d R ) + c θ z r L r R ( θ ˙ z r L θ ˙ z r R ) + k θ z r L r R ( θ z r L θ z r R ) ,
where kxrLrR, kyrLrR, kzrLrR and kθzrLrR represent the support stiffnesses of the bearing in the respective directions, and cxrLrR, cyrLrR, czrLrR and cθzrLrR represent the support dampings of the bearing.
The motion equation for the left-side ring gear rL is as follows:
m r ( x ¨ r L 2 ω h y ˙ r L ω h 2 x r L ) i = 1 N F r p i L sin ψ r p i cos β b + F x d R r L r R = 0 m r ( y ¨ r L + 2 ω h x ˙ r L ω h 2 y r L ) + i = 1 N F r p i L cos ψ r p i cos β b + F y d R r L r R = 0 m r z ¨ r L + i = 1 N F r p i L sin β b + F z d R r L r R = 0 J r θ ¨ z r L i = 1 N F r p i L r b r cos β b + M θ z d R r L r R = 0 ,
The support force of bearing acting on the right-side ring gear rR is as follows:
F b x r R = c b x r R x ˙ r R ω h y r R + k b x r R x r R F b y r R = c b y r R y ˙ r R + ω h x r R + k b y r R y r R F b z r R = c b z r R z ˙ r R + k b z r R z r R M b θ z r R = c b θ z r R θ ˙ z r R + k b θ z r R θ z r R ,
where kbxrR, kbyrR, kbzrR, and kbθzrR represent the support stiffnesses of the bearing in the respective directions, and cbxrR, cbyrR, cbzrR, and cbθzrR represent the support dampings of the bearing.
The support force of the shaft acting on the right-side ring gear, rR, is as follows:
F x r R r L = c x r R r L ( x ˙ r R x ˙ r L ) + k x r R r L ( x r R x r L ) F y r R r L = c y r R r L ( y ˙ r R y ˙ r L ) + k y r R r L ( y r R y r L ) F z r R r L = c z r R r L ( z ˙ r R z ˙ r L ) + k z r R r L ( z r R z r L ) M θ z r R r L = c θ z r R r L ( θ ˙ z r R θ ˙ z r L ) + k θ z r R r L ( θ z r R θ z r L ) ,
The motion equation for the right-side ring gear, rR, is as follows:
m r ( x ¨ r R 2 ω h y ˙ r R ω h 2 x r R ) i = 1 N F r p i R sin ψ r p i cos β b + F b x r R + F x r R r L = 0 m r ( y ¨ r R + 2 ω h x ˙ r R ω h 2 y r R ) + i = 1 N F r p i R cos ψ r p i cos β b + F b y r R + F y r R r L = 0 m r z ¨ r R i = 1 N F r p i R sin β b + F b z r R + F z r R r L = 0 J d θ ¨ z r R i = 1 N F r p i R r b r cos β b + M b θ z r R + M θ z r R r L = T o u t 1 ,

3.3.7. Motion Differential Equation for the Carrier, h

The components of the reactive force acting on the carrier, h, from the planet gear, pi, along the carrier coordinate axes is as follows:
F b x h p i = k b x p δ x h p i + c b x p δ ˙ x h p i F b y h p i = k b y p δ y h p i + c b y p δ ˙ y h p i F b z h p i = k b z p δ z h p i + c b z p δ ˙ z h p i M b θ z h p i = k b θ z p θ z h p i + c b θ z h p i θ ˙ z h p i ,
The support force of the bearing acting on the carrier, h, is as follows:
F b x h = c b x h x ˙ h ω h y h + k b x h x h F b y h = c b y h y ˙ h + ω h x h + k b y h y h F z y h = c z y h z ˙ h + k z y h z h M b θ z h = c b θ z h θ ˙ z h + k b θ z h θ z h ,
where kbxh, kbyh, kbzh, and kbθzh represent the support stiffnesses of the bearing in the respective directions, and cbxh, cbyh, cbzh, and cbθzh represent the support dampings of the bearing.
The motion equation for the carrier, h, is as follows:
m h x ¨ h 2 ω h y ˙ h ω h 2 x h + i = 1 N F b x h p i + F b x h = 0 m h y ¨ h + 2 ω h x ˙ h ω h 2 y h + i = 1 N F b y h p i + F b y h = 0 m h z ¨ h + i = 1 N F b z h p i + F b z h = 0 J h θ z h + r h i = 1 N ( F b y h p i cos ψ p i F b x h p i sin ψ p i ) + M b θ z c = T o u t 2 ,
Taking the TMS, transmission errors, bearing stiffness, and damping into account, the motion differential equations of each subsystem can be assembled to obtain the overall dynamic motion differential equations of the entire system.
M X ¨ + C b + C m X ˙ + K b + K m ω h 2 K w + ω h K d X = F + T ,

3.4. The Definition of LSC

The LSC in gear systems refers to a coefficient employed for the computation of the forces acting on gears within a gear transmission. Owing to the influence of internal and external excitations, the meshing forces undergo continuous variations. Hence, the introduction of the LSC becomes essential to account for such fluctuations in gear transmissions. The LSCs at a specific temporal instance are to be delineated by bamij, bdmij, bspij, and brpij as follows:
b a m i j = 3 F a m i j / i = 1 3 F a m i j , b d m i j = 3 F d m i j / i = 1 3 F d m i j b s p i j = 3 F s p i j / i = 1 3 F s p i j , b r p i j = 3 F r p i j / i = 1 3 F r p i j ,
where i = 1, 2, 3; j = L, R.
Take the maximum values of the LSCs of each planet gear at each time point as the LSCs for the internal and external meshing of the planet gears at that moment; then, the dynamic LSCs can be defined as bamj, bdmj, bspj, and brpj as follows:
b a m j = max ( b a m i j ) , b d m j = max ( b d m i j ) b s p j = max ( b s p i j ) , b r p j = max ( b r p i j ) ,
Take the maximum values of the dynamic LSCs within a period of time as the LSCs as follows:
B a m j = max ( b a m j ) , B d m j = max ( b d m j ) B s p j = max ( b s p j ) , B r p j = max ( b r p j ) ,

4. Analysis of Dynamic Response Characteristics of the Transmission System

According to Reference [25], during the gear meshing process, the TMS and mesh damping are calculated. Considering the gear errors with a precision level of 5, transmission system dynamic response characteristics are analyzed using the Runge–Kutta algorithm. The basic parameters are listed in Table 1. The input power of the transmission system was 250 kW, and the input speed was 10,000 rpm.

4.1. Dynamic Meshing Forces

4.1.1. Dynamic Meshing Forces under Rated Condition

Calculations were performed to determine the time-domain response curves of the dynamic meshing forces for the fixed-axes PGT and the differential PGT on the left and right sides of the pairings. The respective response curves are shown in Figure 7 and Figure 8. As an example, the frequency-domain curves of the dynamic meshing forces for the external meshing pairs of on the left side are given, as shown in Figure 9.
On the basis of Figure 7 and Figure 8, it is evident that the dynamic meshing forces in the fixed-axes PGT and the differential PGT exhibited significant similarity. They followed the same temporal trend with only slight differences in magnitude. Therefore, the subsequent analysis exemplifies by the left-side meshing pair of the system. In the fixed-axes PGT, the average values of the meshing forces are approximately 425 N, while in the differential PGT, the average values are approximately 850 N. The meshing forces of each meshing pair fluctuate around their respective average values, and there are differences among them due to geometric relationships and meshing phases.
From Figure 9, it can be observed that due to the constant external load, the maximum amplitude of the meshing forces occurs at the 0 Hz position. In the fixed-axes PGT, the main frequency components are the rotation frequency, fa, of the sun gear, a, and the meshing frequency, f1, of the fixed-axes PGT and its second harmonic, 2f2. In the differential PGT, the main frequency components are the rotational frequency, f s h , of the sun gear and the meshing frequency, f2, of the differential PGT and its second harmonic, 2f2. This is caused by the periodic fluctuations in mesh stiffness and transmission errors.

4.1.2. Effect of External Excitation on Dynamic Meshing Forces

In practical conditions, the transmission system is subjected to complex and varying external loads rather than a single, steady torque. Studying the dynamic response characteristics of the transmission system under varying conditions and load fluctuations helps in comprehending its performance in practical applications when the system is subjected to alternating external loads, as represented by the Fourier series as follows:
T o u t = 596.875 + 50 cos f t 1 t + π 3 + 25 cos f t 2 t + π 2 ,
where ft1 is w2/3 and ft2 is w2/5.
The time-domain response curves of the meshing forces for the fixed-axes PGT and the differential PGT are shown in Figure 10 and Figure 11, respectively. As an example, the frequency-domain curves of the dynamic meshing forces for the external meshing pairs on the left side are given, as shown in Figure 12.
In Figure 10 and Figure 11, from the time-domain responses perspective, the overall trends of the dynamic meshing forces considering external loads are consistent with those under rated condition. This similarity is attributed to the variations in TMS and transmission errors. However, there are noticeable differences due to the alternating nature of the external loads. The fluctuation of the meshing forces becomes more prominent, resulting in sharper peaks compared to the meshing forces under the rated condition. The frequent variations in external loads may accelerate fatigue damage on the tooth surfaces. From the frequency-domain responses shown in Figure 12, compared to the frequency-domain response under the rated condition, there are additional frequency components in the spectrum due to the external loads, ft1 and ft2.

4.2. Analysis of Planet Gear Vibration Displacements

The vibration of the planet gears has the potential to change both the internal and external forces involved in the meshing process. Hence, it is imperative to study the vibrational characteristics of the planet gears.

4.2.1. Vibration Displacements of the Planet Gears under Rated Condition

The translational displacements and the torsional displacements of the planet gears are shown in Figure 13 and Figure 14.
According to Figure 2, in the fixed-axes PGT, the carrier, c, bears a torque of Tc = −1, while in the differential PGT, the carrier, h, bears a torque of T o u t 2 = t 2 1 + t . Thus, in Figure 13 and Figure 14, the radial displacements of the planet gears, m, are less than those of planet gears, p. In the x-direction, internal and external meshing forces cancel one another out, whereas in the y-direction, they superimpose. Additionally, in the differential PGT, there is a driving action on the planet gears, p, along the tangential direction of the carrier, resulting in larger displacements in the y-direction. Along the axial direction, the influence of transmission errors on the meshing forces is relatively small, as the errors primarily affect the radial direction. The primary cause of the axial vibrations is the TMS, with a frequency equal to the meshing frequency. In terms of the torsional direction, the planet gears oscillate around their equilibrium positions with an amplitude on the order of 10−4 rad. The frequency of the torsional vibration is a combination of transmission errors and time-varying meshing stiffness. Because of the geometric relationships and meshing phases differences, the amplitude and variation tendencies of the individual planet gears exhibit discrepancies.

4.2.2. Vibration Displacements of Planet gears under Different Input Torques

The vibration displacements of the planet gears are subject to variation when there are changes in the input torque. To investigate the vibration displacements of the planet gears in response to varying input torques, five operating conditions were set, as shown in Table 2. For each of the five operating conditions, 80 meshing cycles were selected, and the displacements of the planet gears were analyzed. The time-domain responses are shown in Figure 15 and Figure 16. The peak-to-peak values and variances of the vibration displacements are shown in Figure 17 and Figure 18, respectively.
From Figure 15 and Figure 16, it is observable that as the torque increases, except in the x-direction, the average displacements of the planet gears in various directions basically tended to increase. The displacements of the planet gears, m, in the x-direction are characterized by their distinctiveness due to their exceptionally small magnitude. Under the five operating conditions, the maximum magnitude remained at around 10−16 m. Therefore, even slight variations in parameters can lead to noticeable changes in its displacement. However, because of its small amplitude, this vibration can be practically negligible. Figure 17 and Figure 18 demonstrate that the peak-to-peak displacements and variances in the planet gears exhibited no discernible patterns in the x-direction. The peak-to-peak displacements and variances in the planet gear, m, in the y-direction exhibited a decrease as the torque increased, suggesting a reduction in the vibration amplitude. Conversely, the trend was opposite for the planet gear, p. In the axial and torsional directions, both the peak-to-peak displacements and variances in the planet gears, m and p, increased with increasing torque, indicating that the vibration became more pronounced in these directions as the torque increased.

4.3. Load-Sharing Coefficients

Using the previously described method in Section 3.4, the dynamic meshing forces of the transmission system were utilized to calculate the LSCs based on the analysis of the dynamic response.

4.3.1. LSCs under Rated Condition

Taking the dynamic LSCs on the left side as an example, the variation curves of the LSCs in the fixed-axes PGT and the differential PGT under the rated condition are shown in Figure 19 and Figure 20, respectively.
From Figure 19 and Figure 20, it can be observed that the variations in the LSCs and the meshing forces were generally consistent, fluctuating at around 1. The fluctuation range for the fixed-axes PGT was between 0.4 and 1.6, while for the differential PGT it is between 0.75 and 1.25. Because of the influence of the geometric relationships and mesh phases, there are differences in the LSCs among the meshing pairs. Since the fixed-axes PGT transmits lower torque compared to the differential PGT, the smaller meshing forces were more susceptible to fluctuations caused by errors. As a result, the fixed-axes PGT exhibited larger fluctuations in the LSCs compared to the differential PGT.

4.3.2. Dynamic LSCs under External Excitation

When subjected to external excitation, taking the dynamic LSCs on the left side as an example, the variations in the LSCs in the fixed-axes PGT and the differential PGT are shown in Figure 21 and Figure 22, respectively. The comparative analysis of the dynamic LSCs considering the impact of external excitation versus the absence of external excitation is shown in Figure 23.
From Figure 21 and Figure 22, it can be observed that under the influence of internal excitations, the overall trend of the LSCs was similar to that under the rated condition. However, the LSCs become sharper due to the alternating external loads. Figure 23 indicates that considering the external excitation results in higher LSCs compared to the case without external excitation, indicating an uneven load distribution caused by the external excitation.

4.3.3. Dynamic LSCs under Different Input Torques

Using the dynamic LSCs on the left side as an example, the variations in the LSCs in the fixed-axes PGT and the differential PGT under the different input torques are shown in Figure 24 and Figure 25, respectively. The comparison of the dynamic LSCs under different input torques are shown in Figure 26.
From Figure 24 and Figure 25, it can be observed that the LSCs exhibited a sinusoidal-like variation pattern, fluctuating approximately every 18 meshing cycles. As the torque increased, the average values and magnitudes of the fluctuations of the LSCs exhibited a decrease. This is consistent with earlier observation that the smaller meshing forces were more susceptible to error-induced fluctuations. Figure 26 shows that with an increasing torque, the mean value of the LSC decreased from around 1.74 to around 1.23 for the fixed-axes PGT and from around 1.26 to around 1.08 for the differential PGT.

5. Model Validation

After analyzing the operating principles of the transmission system and the motion relationships among its components, all parts were treated as rigid bodies in the ADAMS dynamic simulation. All rotating components were constrained as revolute joints, and contact pairs were added between meshing gears.
To validate the results of this study, a comparison was made between the meshing forces obtained from the ADAMS dynamic simulation and those obtained from the dynamic model. After analyzing the operating principles of the transmission system and the motion relationships among its components, all components were treated as rigid bodies in the ADAMS dynamic simulation. All rotating components were constrained as revolute joints, and contact pairs were added between meshing gears. A driving torque was applied to the sun gear, aL, while a load torque was applied to the ring gear, rR, and the carrier, h. The operating conditions were consistent with those in Section 4.1, and the constraints of the transmission system are illustrated in Figure 27.
The simulation duration was set to 0.0218 s. Taking the meshing forces between the sun gear and planet gears on the left side of both the fixed-axes PGT and the differential PGT as examples, the meshing force curves were obtained, which are shown in Figure 28.
By comparing Figure 28a with Figure 7a and Figure 28b with Figure 8a, it can be observed that the meshing force curves obtained from the ADAMS dynamic simulation are generally consistent with those calculated by the model in terms of the numerical values and trend variations. Within the allowable range of error, this validates the reliability of the theoretical calculations.

6. Conclusions

A comprehensive dynamic model of the coaxial counter-rotating planetary transmission system is constructed. The model incorporated TMS, transmission errors, and external excitations. The primary objective of the investigation was to conduct an analysis of the meshing forces, displacements of the planet gears, and the performance of the load distribution in the transmission system.
Firstly, considering the translational and torsional vibrations of each component, a dynamic model of the coaxial counter-rotating planetary transmission system was established using the concentrated mass method. The dynamic equations of each component were derived based on the relative displacements between gear teeth.
Secondly, the dynamic meshing forces, planet gear displacements, and load-sharing characteristics of the coaxial counter-rotating planetary transmission system were investigated. The meshing forces acting on both sides exhibited significant similarities. When considering the external excitations, the variations in the meshing force curves became more pronounced. The radial displacements of the planet gears in the differential PGT were greater than that in the fixed-axes PGT. With an increasing input torque, the average displacements of planet gears in all directions tended to increase. The differential PGT, transmitting higher power, demonstrated a better load-sharing performance compared to the fixed-axes PGT.
In summary, this study explored the dynamic responses of the coaxial counter-rotating planetary transmission system under different operating conditions. It is noted that while the quantitative results presented are specific to the system studied, the qualitative findings are generally applicable to coaxial counter-rotating planetary transmission systems. Although the model established by the concentrated mass method may lack the accuracy of a finite element analysis, its computational efficiency is high, and the results are also close to reality. When performing strength, fatigue, reliability analysis, or parameter optimization for coaxial counter-rotating planetary transmission systems, this method can be utilized to quickly obtain data, such as meshing forces, vibration displacements, and load-sharing coefficients, serving as the basis for further work. Mechanical system designers can also use the main findings of this study as general guidelines for designing coaxial counter-rotating planetary transmission systems.
Additionally, we recognize the importance of analytical methods in system analysis and optimization. In the future, we plan to combine numerical methods with analytical methods, based on simplified assumptions, to simplify and extend the model developed in this study and to obtain exact solutions and optimization results under certain conditions. This approach not only reduces computational complexity but also provides a theoretical foundation for further optimization. We will continue to explore these directions in future research to enhance and refine the study of the dynamic characteristics of transmission systems.

Author Contributions

Conceptualization, Z.Y. and Z.C.; methodology, Z.Y.; software, J.Q.; validation, Y.L., M.D., S.M. and G.Y.; formal analysis, Z.Y.; investigation, Z.Y.; resources, Y.L.; data curation, M.D.; writing—original draft preparation, Z.Y.; writing—review and editing, Z.Y.; visualization, S.M.; supervision, G.Y.; project administration, G.Y.; funding acquisition, G.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Heilongjiang Key R&D Project (grant number: GA21D004), Heilongjiang Major Scientific and Technological Achievements Transformation Project (grant number: CG21B010), and Harbin University of Technology Scientific Research Initiation Fund Project (grant number: FRFCU5710052921).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to laboratory regulations.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Boguski, B.; Kahraman, A. An experimental study on the motion transmission error of planetary gear sets. Comput. Methods. Appl. Mech. Engrg. 2015, 283, 956–970. [Google Scholar]
  2. Shen, Z.; Qiao, B.; Yang, L.; Luo, W.; Yang, Z.; Chen, X. Fault mechanism and dynamic modeling of planetary gear with gear wear. Mech. Mach. Theory 2021, 155, 104098. [Google Scholar] [CrossRef]
  3. Xu, H.; Ren, H.; Qin, D. Dynamic characteristics of the planetary gear system with rolling bearing. Multibody Syst. Dyn. 2023, 59, 171–191. [Google Scholar] [CrossRef]
  4. Gu, X.; Velex, P. A dynamic model to study the influence of planet position errors in planetary gears. J. Sound Vib. 2012, 331, 4554–4574. [Google Scholar] [CrossRef]
  5. Gu, X.; Velex, P. On the dynamic simulation of eccentricity errors in planetary gears. Mech. Mach. Theory 2013, 61, 14–29. [Google Scholar] [CrossRef]
  6. 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]
  7. Huangfu, Y.F.; Dong, X.J.; Chen, K.K.; Li, Z.; Peng, Z. An insight into the pass effect of the planet gear from an elastodynamics perspective. Sci. China Technol. Sci. 2023, 66, 2415–2431. [Google Scholar] [CrossRef]
  8. Han, H.; Ma, H.; Tian, H.; Peng, Z.; Zhu, J.; Li, Z. Sideband analysis of cracked planetary gear train considering output shaft radial assembly error. Mech. Syst. Signal Process 2023, 200, 110618. [Google Scholar] [CrossRef]
  9. Wang, Y.; Liao, Y.; Xu, H. Effects of multi-excitation on vibration characteristics of planetary gear system. Alex. Eng. J. 2022, 61, 10593–10602. [Google Scholar] [CrossRef]
  10. Liu, C.; Yang, C.; Zhao, Y.; Luo, J.; Chen, X. Dynamic modeling and analysis of high-speed flexible planetary gear transmission systems. Alex. Eng. J. 2023, 80, 444–464. [Google Scholar] [CrossRef]
  11. Liu, J.; Li, X.; Xia, M. A dynamic model for the planetary bearings in a double planetary gear set. Mech. Syst. Signal Process 2023, 194, 110257. [Google Scholar] [CrossRef]
  12. Fyler, D.C.; Inalpolat, M. A dynamic model for double-planet planetary gearsets. J. Vib. Acoust. 2016, 138, 021006. [Google Scholar] [CrossRef]
  13. Guo, H.; Zhang, J.; Yu, H. Robust optimisation of dynamic and NVH characteristics for compound power-split hybrid transmission. Proc. Inst. Mech. Eng. K J. Multi-Body Dyn. 2019, 233, 817–826. [Google Scholar] [CrossRef]
  14. Tung, L.C.; Chan, Y.J. A time-variant dynamic model for compound epicyclic–cycloidal reducers. Mech. Mach. Theory 2023, 179, 105095. [Google Scholar] [CrossRef]
  15. Zhang, H.; Shen, X. A dynamic tooth wear prediction model for reflecting “two-sides” coupling relation between tooth wear accumulation and load sharing behavior in compound planetary gear set. Proc. Inst. Mech. Eng. C J. Mech. Eng. Sci. 2020, 234, 1746–1763. [Google Scholar] [CrossRef]
  16. Yang, L.; Yuan, B.; Gong, J.; Qin, M.; Liu, G. Dynamic modelling and vibration characteristics of a marine compound gear transmission system. Proc. Inst. Mech. Eng. K J. Multi-Body Dyn. 2023, 237, 261–278. [Google Scholar] [CrossRef]
  17. Li, W.; Li, Z. Dynamic analysis of a multi-stage gear system considering the coupling between mesh phasing angle and coaxial teeth ratio. Nonlinear Dyn. 2023, 111, 19855–19878. [Google Scholar] [CrossRef]
  18. Zhang, D.; Zhu, R.; Fu, B.; Tan, W. Modal properties of contra-rotating encased differential gear train used in coaxial helicopter. J. Vib. Eng. Technol. 2020, 8, 799–814. [Google Scholar] [CrossRef]
  19. Ryali, L.; Talbot, D. A dynamic gear load distribution model for epicyclic gear sets with a structurally compliant planet carrier. Mech. Mach. Theory 2023, 181, 105225. [Google Scholar] [CrossRef]
  20. Lai, J.; Liu, Y.; Xu, X.; Li, H.; Xu, J.; Wang, S.; Guo, W. Dynamic modeling and analysis of Ravigneaux planetary gear set with unloaded floating ring gear. Mech. Mach. Theory 2022, 170, 104696. [Google Scholar] [CrossRef]
  21. Wang, D.; Zhao, Y.; Ren, J.; Yu, S. Dynamic characteristics of the ring gear structure of two-stage plastic planetary reducers. J. Braz. Soc. Mech. Sci. Eng. 2023, 45, 474. [Google Scholar] [CrossRef]
  22. Cui, T.; Li, Y.; Zan, C.; Chen, Y. Dynamic modeling and analysis of nonlinear compound planetary system. Machines 2022, 10, 31. [Google Scholar] [CrossRef]
  23. Zhang, J.; Guo, H.; Yu, H.; Zhang, T. Numerical and experimental investigation on nonlinear dynamic characteristics of planetary gear train. J. Theor. Appl. Mech. 2020, 58, 1009–1022. [Google Scholar] [CrossRef]
  24. Peneva, M.; Radkova, K.; Troha, S.; Karaivanov, D. On the using of two-carrier planetary gear trains with two compound and four external shafts as change-gears. MATEC Web Conf. 2022, 366, 01007. [Google Scholar] [CrossRef]
  25. Shuai, M.; Ting, Z.; Guo-Guang, J.; Xiao-Lin, C.; Han-Jun, G. Analytical investigation on load sharing characteristics of herringbone planetary gear train with flexible support and floating sun gear. Mech. Mach. Theory 2020, 144, 103670. [Google Scholar] [CrossRef]
Figure 1. Diagrammatic representation and model of the coaxial counter-rotating planetary transmission system.
Figure 1. Diagrammatic representation and model of the coaxial counter-rotating planetary transmission system.
Applsci 14 04491 g001
Figure 2. Diagram of the coaxial counter-rotating planetary transmission system represented by the Wolf–Arnaudov symbol.
Figure 2. Diagram of the coaxial counter-rotating planetary transmission system represented by the Wolf–Arnaudov symbol.
Applsci 14 04491 g002
Figure 3. Dynamic model of the system.
Figure 3. Dynamic model of the system.
Applsci 14 04491 g003
Figure 4. Relationship between the sun gear and planet gears in the differential PGT.
Figure 4. Relationship between the sun gear and planet gears in the differential PGT.
Applsci 14 04491 g004
Figure 5. Relationship between the ring gear and planet gears in the differential PGT.
Figure 5. Relationship between the ring gear and planet gears in the differential PGT.
Applsci 14 04491 g005
Figure 6. Relationship between the planet gear and the planet carrier in the differential PGT.
Figure 6. Relationship between the planet gear and the planet carrier in the differential PGT.
Applsci 14 04491 g006
Figure 7. Time-domain response curves of dynamic meshing forces in the fixed-axes PGT under rated condition: (a) left external meshing pairs; (b) left internal meshing pairs; (c) right external meshing pairs; (d) right internal meshing pairs.
Figure 7. Time-domain response curves of dynamic meshing forces in the fixed-axes PGT under rated condition: (a) left external meshing pairs; (b) left internal meshing pairs; (c) right external meshing pairs; (d) right internal meshing pairs.
Applsci 14 04491 g007
Figure 8. Time-domain response curves of dynamic meshing forces in the differential PGT under rated condition: (a) left external meshing pairs; (b) left internal meshing pairs; (c) right external meshing pairs; (d) right internal meshing pairs.
Figure 8. Time-domain response curves of dynamic meshing forces in the differential PGT under rated condition: (a) left external meshing pairs; (b) left internal meshing pairs; (c) right external meshing pairs; (d) right internal meshing pairs.
Applsci 14 04491 g008
Figure 9. Frequency-domain curves of the dynamic meshing forces for the transmission system under rated condition: (a) fixed-axes PGT; (b) differential PGT.
Figure 9. Frequency-domain curves of the dynamic meshing forces for the transmission system under rated condition: (a) fixed-axes PGT; (b) differential PGT.
Applsci 14 04491 g009
Figure 10. Time-domain response curves of the dynamic meshing forces in the fixed-axes PGT under external excitation: (a) left external meshing pairs; (b) left internal meshing pairs.
Figure 10. Time-domain response curves of the dynamic meshing forces in the fixed-axes PGT under external excitation: (a) left external meshing pairs; (b) left internal meshing pairs.
Applsci 14 04491 g010
Figure 11. Time-domain response curves of the dynamic meshing forces in the differential PGT under external excitation: (a) left external meshing pairs; (b) left internal meshing pairs.
Figure 11. Time-domain response curves of the dynamic meshing forces in the differential PGT under external excitation: (a) left external meshing pairs; (b) left internal meshing pairs.
Applsci 14 04491 g011
Figure 12. Frequency-domain curves of the dynamic meshing forces for the transmission system under external excitation: (a) fixed-axes PGT; (b) differential PGT.
Figure 12. Frequency-domain curves of the dynamic meshing forces for the transmission system under external excitation: (a) fixed-axes PGT; (b) differential PGT.
Applsci 14 04491 g012
Figure 13. Time-domain responses of the planet gear displacements in the fixed-axes PGT under the rated condition: (a) x-direction; (b) y-direction; (c) z-direction; (d) torsional direction around z-axis.
Figure 13. Time-domain responses of the planet gear displacements in the fixed-axes PGT under the rated condition: (a) x-direction; (b) y-direction; (c) z-direction; (d) torsional direction around z-axis.
Applsci 14 04491 g013
Figure 14. Time-domain responses of the planet gear displacements in the differential PGT under the rated condition: (a) x-direction; (b) y-direction; (c) z-direction; (d) torsional direction around z-axis.
Figure 14. Time-domain responses of the planet gear displacements in the differential PGT under the rated condition: (a) x-direction; (b) y-direction; (c) z-direction; (d) torsional direction around z-axis.
Applsci 14 04491 g014
Figure 15. Time-domain responses of the planet gear displacements in the fixed-axes PGT under the different input torques: (a) x-direction; (b) y-direction; (c) z-direction; (d) torsional direction around z-axis.
Figure 15. Time-domain responses of the planet gear displacements in the fixed-axes PGT under the different input torques: (a) x-direction; (b) y-direction; (c) z-direction; (d) torsional direction around z-axis.
Applsci 14 04491 g015
Figure 16. Time-domain responses of the planet gear displacements in the differential PGT under the different input torques: (a) x-direction; (b) y-direction; (c) z-direction; (d) torsional direction around z-axis.
Figure 16. Time-domain responses of the planet gear displacements in the differential PGT under the different input torques: (a) x-direction; (b) y-direction; (c) z-direction; (d) torsional direction around z-axis.
Applsci 14 04491 g016
Figure 17. Peak-to-peak displacements of the planet gears under the different input torques: (a) x-direction; (b) y-direction; (c) z-direction; (d) torsional direction around z-axis.
Figure 17. Peak-to-peak displacements of the planet gears under the different input torques: (a) x-direction; (b) y-direction; (c) z-direction; (d) torsional direction around z-axis.
Applsci 14 04491 g017
Figure 18. Variances of the planet gear displacements under the different input torques: (a) x-direction; (b) y-direction; (c) z-direction; (d) torsional direction around z-axis.
Figure 18. Variances of the planet gear displacements under the different input torques: (a) x-direction; (b) y-direction; (c) z-direction; (d) torsional direction around z-axis.
Applsci 14 04491 g018
Figure 19. Variations in the dynamic LSCs in the fixed-axes PGT under the rated condition: (a) left external meshing pairs; (b) left internal meshing pairs.
Figure 19. Variations in the dynamic LSCs in the fixed-axes PGT under the rated condition: (a) left external meshing pairs; (b) left internal meshing pairs.
Applsci 14 04491 g019
Figure 20. Variations in the dynamic LSCs in the differential PGT under the rated condition: (a) left external meshing pairs; (b) left internal meshing pairs.
Figure 20. Variations in the dynamic LSCs in the differential PGT under the rated condition: (a) left external meshing pairs; (b) left internal meshing pairs.
Applsci 14 04491 g020
Figure 21. Variations in dynamic LSCs in the fixed-axes PGT under external excitation. (a) Left external meshing pairs; (b) Left internal meshing pairs.
Figure 21. Variations in dynamic LSCs in the fixed-axes PGT under external excitation. (a) Left external meshing pairs; (b) Left internal meshing pairs.
Applsci 14 04491 g021
Figure 22. Variations in the dynamic LSCs in the differential PGT under external excitation: (a) left external meshing pairs; (b) left internal meshing pairs.
Figure 22. Variations in the dynamic LSCs in the differential PGT under external excitation: (a) left external meshing pairs; (b) left internal meshing pairs.
Applsci 14 04491 g022
Figure 23. Comparison of the dynamic LSCs of the transmission system with and without external excitation.
Figure 23. Comparison of the dynamic LSCs of the transmission system with and without external excitation.
Applsci 14 04491 g023
Figure 24. Variations in the dynamic LSCs in the fixed-axes PGT under the different input torques: (a) left external meshing pairs; (b) left internal meshing pairs.
Figure 24. Variations in the dynamic LSCs in the fixed-axes PGT under the different input torques: (a) left external meshing pairs; (b) left internal meshing pairs.
Applsci 14 04491 g024
Figure 25. Variations in the dynamic LSCs in the differential PGT under the different input torques: (a) left external meshing pairs; (b) left internal meshing pairs.
Figure 25. Variations in the dynamic LSCs in the differential PGT under the different input torques: (a) left external meshing pairs; (b) left internal meshing pairs.
Applsci 14 04491 g025
Figure 26. Comparison of the dynamic LSCs of the transmission system under the different input torques.
Figure 26. Comparison of the dynamic LSCs of the transmission system under the different input torques.
Applsci 14 04491 g026
Figure 27. Constraints of the transmission system in the ADAMS dynamic simulation.
Figure 27. Constraints of the transmission system in the ADAMS dynamic simulation.
Applsci 14 04491 g027
Figure 28. Meshing force curves obtained from the ADAMS dynamic simulation: (a) left external meshing forces of the fixed axes PGT; (b) left external meshing forces of the differential PGT.
Figure 28. Meshing force curves obtained from the ADAMS dynamic simulation: (a) left external meshing forces of the fixed axes PGT; (b) left external meshing forces of the differential PGT.
Applsci 14 04491 g028
Table 1. Parameters of the transmission system.
Table 1. Parameters of the transmission system.
ItemNumber of TeethAngle of Helix (°)Normal Pressure Angle (°)ModuleMass (kg)Moment of Inertia (kg·mm2)
Sun gear, a2220201.50.09217.664
Planet gear, m4420201.50.327268.521
Ring gear, dL11020201.50.1631325.860
Ring gear, dR11020201.50.1631325.860
Sun gear, s5520201.50.488525.752
Planet gear, p2720201.50.10835.531
Ring gear, rL11020201.50.1631325.860
Ring gear, rR11020201.50.1631325.860
Planet carrier, h 1.6484946.014
Table 2. Five different operating conditions.
Table 2. Five different operating conditions.
Torque (Nm)Speed (rpm)
Case 120011,938
Case 23007958
Case 34005968
Case 45004775
Case 56003979
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yue, Z.; Chen, Z.; Qu, J.; Li, Y.; Dzianis, M.; Mo, S.; Yu, G. Analysis of the Dynamic Characteristics of Coaxial Counter-Rotating Planetary Transmission System. Appl. Sci. 2024, 14, 4491. https://doi.org/10.3390/app14114491

AMA Style

Yue Z, Chen Z, Qu J, Li Y, Dzianis M, Mo S, Yu G. Analysis of the Dynamic Characteristics of Coaxial Counter-Rotating Planetary Transmission System. Applied Sciences. 2024; 14(11):4491. https://doi.org/10.3390/app14114491

Chicago/Turabian Style

Yue, Zongxiang, Zhaobo Chen, Jianjun Qu, Yang Li, Marmysh Dzianis, Shuai Mo, and Guangbin Yu. 2024. "Analysis of the Dynamic Characteristics of Coaxial Counter-Rotating Planetary Transmission System" Applied Sciences 14, no. 11: 4491. https://doi.org/10.3390/app14114491

APA Style

Yue, Z., Chen, Z., Qu, J., Li, Y., Dzianis, M., Mo, S., & Yu, G. (2024). Analysis of the Dynamic Characteristics of Coaxial Counter-Rotating Planetary Transmission System. Applied Sciences, 14(11), 4491. https://doi.org/10.3390/app14114491

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