# Global Analysis of a Planetary Gear Train

Authors: Tongjie Li; Rupeng Zhu
Journal: Shock and Vibration
 (2014)

---

## Abstract

By using the Poincaré-like cell-to-cell mapping method and shooting method, the global characteristics of a planetary gear train are studied based on the torsional vibration model with errors of transmission, time-varying meshing stiffness, and multiple gear backlashes. The study results reveal that the planetary with a certain set of parameters has four coexisting periodic orbits, which are P-1, P-2, P-4, and P-8,  respectively. P-1 and P-2 motions are not of long-term stability, P-8 motion is of local stability, and P-4 motion is of global stability. Shooting method does not have the capacity of searching coexisting periodic orbits in a global scope, and it is easy to omit some periodic orbits which are far away from the main gropes of periodic orbits.


---

## Body

## 1. Introduction

Because of the existence of some nonlinear factors, such as multiple gear backlashes and time-varying mesh stiffness, planetary gear train is an inherent strongly nonlinear dynamical system. Generally speaking, a nonlinear dynamical system with a certain set of parameters may have several periodic orbits [1], and different orbits may have different domains of attraction [2]. Making certain of all coexisting periodic orbits and their domains of attraction, designers of a planetary gear train can forecast the motion state of the system and the scope of disturbance that the motion state can bear.Over the past three decades, the dynamics of planetary transmissions have drawn much attention, but almost all of the published studies on this field focused only on their linear vibration property [3–6]; the literature about nonlinear dynamics of a planetary transmission system is still limited [7, 8]. As to global analysis of a planetary gear train, we can not find a study on this aspect in the existing literature.Some scholars [9, 10] made development research about a method called cell-to-cell mapping method, which has demonstrated their great potential and efficiency to study the global characteristics of a general nonlinear dynamical system. In aspect of application of cell-to-cell mapping method to engineering, Liu and Chen [11] studied the global dynamic characteristics of a rotor-bearing nonlinear dynamical system by using Poincaré-like cell-to-cell mapping method. Polchai and Hsu [12] studied domains of stability of a synchronous generator by a cell-to-cell method.The objective of this paper is to study all of the coexisting periodic solutions and their domains of attraction of a planetary gear train based on a nonlinear torsional vibration model with multiple backlashes, time-varying mesh stiffness, and static transmission errors.


## 2. Dynamic Model of a Planetary Gear Train

In this study, the planetary gear system concerned is a single-stage 2K-H type planetary gear set as shown in Figure1. Some detailed modeling progress has been published in [8] by the author of this paper. The model consists of a sun gear S, a ring gear R, planet gears P, and a carrier C. Translational degrees of every components in this model are eliminated, and only rotational motions of every component bodies are considered. Rotational motions of the ring, carrier gear, sun gear, and planet gears are denoted by θh, h=r, c, s, p1, p2, …, pN, where N indicates the number of planet gears. In this paper, θr=0 since the ring gear is fixed on gearbox. The character rbh, h=s, r, p1, p2, …, pN, c, stands for the base circle radius for the sun, ring, planets, and the radius of the circle passing through the planet centers for the carrier, respectively.Figure 1
Nonlinear dynamical model of planetary gear train.k
s
p
i, cspi, bspi, and espi denote the backlash, time-varying mesh stiffness, viscous damping and the static transmission error of the sun gear-planet gear-i pair; krpi, crpi, brpi, and erpi denote the corresponding terms in the ring gear-planet gear-i pair. This nomenclature also applies to other physical parameters of the corresponding meshing pairs.


### 2.1. Time-Varying Mesh Stiffness

For spur gears, rectangular waves are often used to approximate mesh stiffness alternating betweend and d+1 pairs of teeth in contact [13]. Expanding rectangular wavesin Fourier series and taking the first harmonic term of the system’s mesh frequency, the time-varying mesh stiffness in corresponding mesh pair can be assumed as
(1)kspi(t)=kmspi+kaspisin⁡(ωt+φspi),krpi(t)=kmrpi+karpisin⁡(ωt+φrpi),
where kmspi and kmrpi are mean values of the stiffness in the corresponding meshing pair, kaspi and karpi are time-varying amplitude of the stiffness in the corresponding meshing pair, φspi and φrpi are initial phases of the stiffness in the corresponding meshing pair, and ω is the system’s mesh frequency.


### 2.2. Static Transmission Errors

The static transmission error is one of the primary excitations of planetary gear train’s vibration and noise. According to the study of Parker [14], the static transmission errors in the corresponding mesh pair espi(t) and erpi(t) can be expressed as
(2)espi(t)=easpisin⁡(ωt+βspi),erpi(t)=earpisin⁡(ωt+βrpi+γsr),
where easpi and erspi are time-varying amplitude of the static transmission error in the corresponding meshing pair, βspi and βrpi are initial phases of the static transmission error in the corresponding meshing pair, and γsr denotes the phase between the sun-planet and ring-planet mesh error for a given planet.


### 2.3. Nonlinear Displacement Function about Backlash

According to the study of Sun and Hu [7], the nonlinear displacement function about backlash is as follows:
(3)f(X,b)={X-bX>b0|X|≤bX+bX<-b,
where X denotes the relative gear mesh displacement and b denotes half backlash in the corresponding mesh pair.


### 2.4. Equations of Motion

In order to eliminate rigid body motion and simplify equations of the system,N1 (in this paper, N1=N+1) dimension equivalent displacements coordinate Xspi and Xsc are introduced according to the model shown in Figure 1:
(4)Xspi=θsrbs-θpirbpi-θcrccos⁡⁡α-espi(t),Xsc=θsrbs-2θcrccos⁡⁡α,
where α is the pressure angle; in this paper all the pressure angles of different gear are of the same value.The dimensionless torsional motion equations of the planetary gear system, which are shown as (5), can be established by using the Lagrange principle and equivalent displacements coordinate (4):
(5)X-¨spi+1ωn2meq1×∑i=13[kmspi+kaspisin⁡(τΩ+βsi)]f(X-spi,b-spi)+1ωnmeq1∑i=13cspiX-˙spi+1ωn2mpi×[kmspi+kaspisin⁡(τΩ+βsi)]f(X-spi,b-spi)+1ωnmpicspiX-˙spi-1ωn2mpi×[kmrpi+karpisin⁡(τΩ+βsi)]×f(X-sc-X-spi-e-spi-e-rpi,b-rpi)-1ωnmpicrpi(X-˙sc-X-˙spi-e-˙spi-e-˙rpi)+1ωn2mc∑i=13[kmrpi+karpisin⁡(τΩ+βsi)]×f(X-sc-X-spi-e-spi-e-rpi,b-rpi)+1ωnmc∑i=13crpi(X-˙sc-X-˙spi-e-˙spi-e-˙rpi)=FD(msbcωn2)+FL(mcbcωn2)-e-¨spi(τ),X-¨sc+1ωn2meq2∑i=13[kmspi+kaspisin⁡(τΩ+βsi)]×f(X-spi,b-spi)+1ωnmeq2×∑i=13cspiX-˙spi2ωnmc×∑i=13crpi(X-˙sc-X-˙spi-e-˙spi-e-˙rpi)+2ωn2mc∑i=13[kmrpi+karpisin⁡(τΩ+βsi)]×f(X-sc-X-spi-e-spi-e-rpi,b-rpi)=FD(msbcωn2)+2FL(mcbcωn2).
Some dimensionless parameters are introduced in the process of deducing (5):Ω
=
ω
/
ω
n, τ=ωnt, ωn=kmsp1/meq1, meq1=mcms/(mc+ms), meq2=mcms/(mc+2ms), ms=Js/rbs2, mpi=Jpi/rbpi2, mc=(Jc+∑i=13mpirc2)/rbc2, FD=TD/rbs, FL=TL/rbs, X=X-bc, X˙=X-˙bcωn, X¨=X-¨bcωn2, b=bcb-, e=bce-, where bc is a character length, Js, Jpi, Jc are the inertia of the sun gear, ith planet gear, and the carrier.


## 2.1. Time-Varying Mesh Stiffness

For spur gears, rectangular waves are often used to approximate mesh stiffness alternating betweend and d+1 pairs of teeth in contact [13]. Expanding rectangular wavesin Fourier series and taking the first harmonic term of the system’s mesh frequency, the time-varying mesh stiffness in corresponding mesh pair can be assumed as
(1)kspi(t)=kmspi+kaspisin⁡(ωt+φspi),krpi(t)=kmrpi+karpisin⁡(ωt+φrpi),
where kmspi and kmrpi are mean values of the stiffness in the corresponding meshing pair, kaspi and karpi are time-varying amplitude of the stiffness in the corresponding meshing pair, φspi and φrpi are initial phases of the stiffness in the corresponding meshing pair, and ω is the system’s mesh frequency.


## 2.2. Static Transmission Errors

The static transmission error is one of the primary excitations of planetary gear train’s vibration and noise. According to the study of Parker [14], the static transmission errors in the corresponding mesh pair espi(t) and erpi(t) can be expressed as
(2)espi(t)=easpisin⁡(ωt+βspi),erpi(t)=earpisin⁡(ωt+βrpi+γsr),
where easpi and erspi are time-varying amplitude of the static transmission error in the corresponding meshing pair, βspi and βrpi are initial phases of the static transmission error in the corresponding meshing pair, and γsr denotes the phase between the sun-planet and ring-planet mesh error for a given planet.


## 2.3. Nonlinear Displacement Function about Backlash

According to the study of Sun and Hu [7], the nonlinear displacement function about backlash is as follows:
(3)f(X,b)={X-bX>b0|X|≤bX+bX<-b,
where X denotes the relative gear mesh displacement and b denotes half backlash in the corresponding mesh pair.


## 2.4. Equations of Motion

In order to eliminate rigid body motion and simplify equations of the system,N1 (in this paper, N1=N+1) dimension equivalent displacements coordinate Xspi and Xsc are introduced according to the model shown in Figure 1:
(4)Xspi=θsrbs-θpirbpi-θcrccos⁡⁡α-espi(t),Xsc=θsrbs-2θcrccos⁡⁡α,
where α is the pressure angle; in this paper all the pressure angles of different gear are of the same value.The dimensionless torsional motion equations of the planetary gear system, which are shown as (5), can be established by using the Lagrange principle and equivalent displacements coordinate (4):
(5)X-¨spi+1ωn2meq1×∑i=13[kmspi+kaspisin⁡(τΩ+βsi)]f(X-spi,b-spi)+1ωnmeq1∑i=13cspiX-˙spi+1ωn2mpi×[kmspi+kaspisin⁡(τΩ+βsi)]f(X-spi,b-spi)+1ωnmpicspiX-˙spi-1ωn2mpi×[kmrpi+karpisin⁡(τΩ+βsi)]×f(X-sc-X-spi-e-spi-e-rpi,b-rpi)-1ωnmpicrpi(X-˙sc-X-˙spi-e-˙spi-e-˙rpi)+1ωn2mc∑i=13[kmrpi+karpisin⁡(τΩ+βsi)]×f(X-sc-X-spi-e-spi-e-rpi,b-rpi)+1ωnmc∑i=13crpi(X-˙sc-X-˙spi-e-˙spi-e-˙rpi)=FD(msbcωn2)+FL(mcbcωn2)-e-¨spi(τ),X-¨sc+1ωn2meq2∑i=13[kmspi+kaspisin⁡(τΩ+βsi)]×f(X-spi,b-spi)+1ωnmeq2×∑i=13cspiX-˙spi2ωnmc×∑i=13crpi(X-˙sc-X-˙spi-e-˙spi-e-˙rpi)+2ωn2mc∑i=13[kmrpi+karpisin⁡(τΩ+βsi)]×f(X-sc-X-spi-e-spi-e-rpi,b-rpi)=FD(msbcωn2)+2FL(mcbcωn2).
Some dimensionless parameters are introduced in the process of deducing (5):Ω
=
ω
/
ω
n, τ=ωnt, ωn=kmsp1/meq1, meq1=mcms/(mc+ms), meq2=mcms/(mc+2ms), ms=Js/rbs2, mpi=Jpi/rbpi2, mc=(Jc+∑i=13mpirc2)/rbc2, FD=TD/rbs, FL=TL/rbs, X=X-bc, X˙=X-˙bcωn, X¨=X-¨bcωn2, b=bcb-, e=bce-, where bc is a character length, Js, Jpi, Jc are the inertia of the sun gear, ith planet gear, and the carrier.


## 3. Introduction about Poincaré-Like Cell-to-Cell Mapping Method

The cell-to-cell mapping method, which is introduced by Hsu [9], is an efficient method that can be used to study the global characteristics of strong nonlinear systems. According to cell-to-cell mapping method, one can divide a large state space of a nonlinear system into finitely many small state spaces, which is called cell space, and study coexisting periodic orbits and their domains of attraction can be turned into the question about mapping relations of finitely many small cells in cell spaces.The main disadvantage of cell-to-cell mapping method is that the computation quantity is too large to study the global characteristics of high dimensional nonlinear dynamical system. Aiming at this defect, Levitas et al. [10] put forward to Poincaré-like cell-to-cell mapping method and bring down the computation quantity of global analysis for a nonlinear dynamical system.The state equations ofn dimension nonautonomous dynamical system with periodicity can be written as follows:
(6)x˙=f(t,x,μ)=f(t+T,x,μ),(t,x,μ)∈R×Rn×R,
where T is the period of the dynamical system and μ is a parameter of the system.Using Poincaré cell-to-cell mapping method, (6) can be changed into a discrete point-to-point mapping in a selected Poincaré section ∑, which is represented as
(7)x(k+1)=p(x(k),μ),k∈Z,
where p: ∑→∑.Searching the periodic solutionsx(t)=x(t+kT,μ), k∈Z is changed into calculating the fixed points x* in the point-to-point mapping (7), which satisfies the following formula
(8)x*=pk(x*,μ),k∈Z.Using point-to-point mapping (7), global analysis of the dynamical system (6) in the whole state space is simplified to calculating the fixed points and their domains of attraction in the selected Poincaré section ∑. Because the dimension of the Poincaré section ∑ is less than the whole state space, the computation quantity to study the global characteristics of a high dimensional dynamical system in ∑ can be reduced significantly.One can reduce the computation quantity significantly again by constructing a cell-to-cell mapping in the Poincaré section∑, which is called Poincaré-like cell-to-cell mapping:
(9)z(m+l)=pl(z(m),μ),l∈Z,
where Z(m) is a cell vector of integer elements which shows the position of the cell in state space and pl means executing the cell-to-cell mapping l times. For a given cell z(m), one can find its central point xd(m) as the start point and carry out point-to-point map (7) l times to get its image point xd(m+l). Judging which cell the map point xd(m+l) is located in, one can get the image cell z(m+1). It is easy to get the coexisting periodic cells and their domains of attraction by judging the relations of cells in cell apace.All of the cells wich the fixed pointsx* are located in make up the periodic cells which reflect the periodic solutions x(t) of dynamical system (6). In specific operation of Poincaré-like cell-to-cell mapping, let mapping times l be a large integer; for example, l=50; the calculated periodic cells will be of long-term stability.All of the cells which are attracted into a set of periodic cells constitute the domain of attraction of these periodic cells, and the magnitude of the domain of attraction reflects the degree of stability of the corresponding periodic solutionsx(t).


## 4. Global Analysis of Planetary Gear Train

### 4.1. Coexisting Periodic Orbits with Long-Term Stability

The main parameters of the planetary gear system showed in Figure1 are as follows m=3 mm, α=20°, easpi=earpi=10 μm, zs=15, zpi=24, zr=63, bc=10 μm, bspi=brpi=50 μm, N=3, kmspi=0.8256×109 N/m, kmrpi=1.06×109 N/m, βsi=βri=0,φspi=φrpi=0, kaspi=0.3×109 N/m, and karpi=0.3×109 N/m.Assuming all of the external meshing pairs are of the same mesh stiffness, static transmission error and backlash, and internal meshing pairs have similar assumptions, and (5) will become differential equation which is of 2 degrees of freedom about the relative displacements. Let sequence of solution vector x(t) of the state equations corresponding to (5) be x(t)=(x1,x2,x3,x4)=(X-spi(t),X-˙spi(t),X-sc(t),X-˙sc(t)).Becausex2 and x4 denote velocity of the relative displacements x1 and x3, respectively, all of the periodic orbits must pass through the following two hyperplane sections: x2=0 and x4=0. In this paper, hyperplane section ∑={(x1,x2,x3,x4)∣x4=0} is regarded as Poincaré section. The studied scope D in Poincaré section ∑ is a rectangular block showed as Figure 2, which is formed by the following 3 coordinate axes, x1, x2, x3. Global analysis in Poincaré section ∑ can bring down the work of calculation immensely.Figure 2
The studied scope D on Poincaré section∑={(x1,x2,x3,x4)∣x4=0}.Considering the position of periodic fixed points which can be solved by using shooting method [15], the studied scope D in Poincaré section ∑ can be determined as
(10)-25≤x1≤25,-10≤x2≤20,-10≤x3≤50.A cell space can be constructed by discretizing the studied scope D in∑ into 50×30×60 cells. In order to express periodic solution of the system that is of m times period of the system conveniently, the simplified expression P-m is introduced in this paper. By using shooting method, 3 coexisting periodic fixed points are searched in D, which are P-1 periodic fixed points p1=[1.2736,-2.7122,5.5372], P-2 periodic fixed points including P21=[0.3119,-2.9195,1.9834] and P22=[2.2329,-2.7524,8.5388], P-4 periodic fixed points including P41=[0.0584,-2.3862,2.5963], P42=[0.9718,-2.8886,7.4437], P43=[1.0232,-3.3769,2.3746], and P44=[2.9127,-2.6107,8.8167], respectively.Starting from P-1 periodic fixed points, integrate the state equation corresponding to (5) in a period of the system by using the fixed-step fourth order Runge-Kutta method; the phase trajectory of the P-1 periodic orbit which is projected into D can be obtained and shown as Figure 3. By using the same method, P-2 and P-4 periodic orbits projected into D can be obtained and shown as Figures 4 and 5, respectively. Asterisks in these figures denote the positions of corresponding periodic fixed points.Figure 3
P-1 fixed points and periodic orbits projected into D.Figure 4
P-2 fixed points and periodic orbits projected into D.Figure 5
P-4 fixed points and periodic orbits projected into D.Through the verification by using direct numerical integration method, one can find that only P-4 motion of these three coexisting periodic orbits is of the long-term stability. Namely, as to P-1 and P-2, there are no domains of attraction.By using Poincaré-like cell-to-cell mapping method, two kinds of coexisting periodic cells which are of long-term stability are found. One is P-4 periodic cell; another is P-4 periodic cell. By contrast, P-4 periodic cell found by Poincaré-like cell-to-cell mapping method and P-4 periodic fixed point found by shooting method stand for the same periodic orbit. P-8 periodic cells comprise the following points in cell space:P81=[-9.8296,3.4673,-8.4349], P82=[-10.0250,3.9632,-2.4987], P83=[-9.3758,3.9118,-6.4858], P84=[-10.1171,3.3448,-4.5905], P85=[-9.8296,3.4673,-8.4349], P86=[-10.0250,3.9632,-2.4987], P87=[-9.3758,3.9118,-6.4858], P88=[-10.1171,3.3448,-4.5905], respectively.The phase trajectory of the P-8 periodic orbit which is projected into D can be obtained and shown as Figure6, and the asterisks in Figure 6 denote the positions of corresponding periodic cells.Figure 6
P-8 fixed points and periodic orbits projected into D.Two conclusions can be gotten from the above study.(1)
With the given set of parameters, planetary gear set has two coexisting periodic orbits with long-term stability; one is P-4 periodic motion, and another is P-8 periodic motion.(2)
The shooting method does not have the capacity of searching coexisting periodic orbits in a global scope.Namely, the shooting method only finds some coexisting periodic orbits which are close to each other and easy to omit some periodic orbits which are far away from the main periodic orbit groups.


### 4.2. Domain of Attraction of Coexisting Periodic Orbits with Long-Term Stability

By using Poincaré-like cell-to-cell mapping method, domain of attraction of P-8 periodic orbits with long-term stability can be calculated and shown as Figure7. The marks “□” denote the position of P-8 periodic cell, the marks “·” denote cells of attraction of P-8 periodic orbits, and all of the cells constitute the domain of attraction of P-8 periodic orbits.Figure 7
Domain of attraction of P-8 periodic orbit.The degree of stability of P-8 periodic orbits can be measured by the scope of domain of attraction; namely, a bigger scope of domain of attraction shows that the corresponding periodic orbits are more stable. If the disturbance of P-8 periodic orbit is limited in its domain of attraction, the motion will keep stable.There are 5470 cells in the domain of attraction of P-8 periodic orbits showed in Figure6, and 6.08% of the whole cell space in D belonged to this domain of attraction. The domain of attraction showed in Figure 7 is composed of three subspaces. Periodic cells of P-8 are full of the three subspaces, but it is not continuous among each of the subspaces.Some conclusions about stability of P-8 periodic orbit can be gotten from the above discussions. Because the domain of attraction is full of the occupied cell spaces, P-8 periodic orbit is of local stability; namely, if disturbance of this orbit is limited in its domain of attraction, the orbit will keep stable, and if disturbance increases to be out of its domain of attraction, the orbit will lose its stability and turn into other coexisting periodic orbits.Figure8 is domain of attraction of P-4 periodic orbit. Except for the cell space which is occupied by domain of attraction of P-8, the domain of attraction of P-4 is almost full of the whole cell space. There are 84530 cells attracted into P-4 periodic orbit; namely, 93.92% of cells in the whole cell space belong to domain of attraction of P-4 periodic orbit. It is obvious that P-4 periodic orbit is of global stability.Figure 8
The domain of attraction of P-4 motion.


## 4.1. Coexisting Periodic Orbits with Long-Term Stability

The main parameters of the planetary gear system showed in Figure1 are as follows m=3 mm, α=20°, easpi=earpi=10 μm, zs=15, zpi=24, zr=63, bc=10 μm, bspi=brpi=50 μm, N=3, kmspi=0.8256×109 N/m, kmrpi=1.06×109 N/m, βsi=βri=0,φspi=φrpi=0, kaspi=0.3×109 N/m, and karpi=0.3×109 N/m.Assuming all of the external meshing pairs are of the same mesh stiffness, static transmission error and backlash, and internal meshing pairs have similar assumptions, and (5) will become differential equation which is of 2 degrees of freedom about the relative displacements. Let sequence of solution vector x(t) of the state equations corresponding to (5) be x(t)=(x1,x2,x3,x4)=(X-spi(t),X-˙spi(t),X-sc(t),X-˙sc(t)).Becausex2 and x4 denote velocity of the relative displacements x1 and x3, respectively, all of the periodic orbits must pass through the following two hyperplane sections: x2=0 and x4=0. In this paper, hyperplane section ∑={(x1,x2,x3,x4)∣x4=0} is regarded as Poincaré section. The studied scope D in Poincaré section ∑ is a rectangular block showed as Figure 2, which is formed by the following 3 coordinate axes, x1, x2, x3. Global analysis in Poincaré section ∑ can bring down the work of calculation immensely.Figure 2
The studied scope D on Poincaré section∑={(x1,x2,x3,x4)∣x4=0}.Considering the position of periodic fixed points which can be solved by using shooting method [15], the studied scope D in Poincaré section ∑ can be determined as
(10)-25≤x1≤25,-10≤x2≤20,-10≤x3≤50.A cell space can be constructed by discretizing the studied scope D in∑ into 50×30×60 cells. In order to express periodic solution of the system that is of m times period of the system conveniently, the simplified expression P-m is introduced in this paper. By using shooting method, 3 coexisting periodic fixed points are searched in D, which are P-1 periodic fixed points p1=[1.2736,-2.7122,5.5372], P-2 periodic fixed points including P21=[0.3119,-2.9195,1.9834] and P22=[2.2329,-2.7524,8.5388], P-4 periodic fixed points including P41=[0.0584,-2.3862,2.5963], P42=[0.9718,-2.8886,7.4437], P43=[1.0232,-3.3769,2.3746], and P44=[2.9127,-2.6107,8.8167], respectively.Starting from P-1 periodic fixed points, integrate the state equation corresponding to (5) in a period of the system by using the fixed-step fourth order Runge-Kutta method; the phase trajectory of the P-1 periodic orbit which is projected into D can be obtained and shown as Figure 3. By using the same method, P-2 and P-4 periodic orbits projected into D can be obtained and shown as Figures 4 and 5, respectively. Asterisks in these figures denote the positions of corresponding periodic fixed points.Figure 3
P-1 fixed points and periodic orbits projected into D.Figure 4
P-2 fixed points and periodic orbits projected into D.Figure 5
P-4 fixed points and periodic orbits projected into D.Through the verification by using direct numerical integration method, one can find that only P-4 motion of these three coexisting periodic orbits is of the long-term stability. Namely, as to P-1 and P-2, there are no domains of attraction.By using Poincaré-like cell-to-cell mapping method, two kinds of coexisting periodic cells which are of long-term stability are found. One is P-4 periodic cell; another is P-4 periodic cell. By contrast, P-4 periodic cell found by Poincaré-like cell-to-cell mapping method and P-4 periodic fixed point found by shooting method stand for the same periodic orbit. P-8 periodic cells comprise the following points in cell space:P81=[-9.8296,3.4673,-8.4349], P82=[-10.0250,3.9632,-2.4987], P83=[-9.3758,3.9118,-6.4858], P84=[-10.1171,3.3448,-4.5905], P85=[-9.8296,3.4673,-8.4349], P86=[-10.0250,3.9632,-2.4987], P87=[-9.3758,3.9118,-6.4858], P88=[-10.1171,3.3448,-4.5905], respectively.The phase trajectory of the P-8 periodic orbit which is projected into D can be obtained and shown as Figure6, and the asterisks in Figure 6 denote the positions of corresponding periodic cells.Figure 6
P-8 fixed points and periodic orbits projected into D.Two conclusions can be gotten from the above study.(1)
With the given set of parameters, planetary gear set has two coexisting periodic orbits with long-term stability; one is P-4 periodic motion, and another is P-8 periodic motion.(2)
The shooting method does not have the capacity of searching coexisting periodic orbits in a global scope.Namely, the shooting method only finds some coexisting periodic orbits which are close to each other and easy to omit some periodic orbits which are far away from the main periodic orbit groups.


## 4.2. Domain of Attraction of Coexisting Periodic Orbits with Long-Term Stability

By using Poincaré-like cell-to-cell mapping method, domain of attraction of P-8 periodic orbits with long-term stability can be calculated and shown as Figure7. The marks “□” denote the position of P-8 periodic cell, the marks “·” denote cells of attraction of P-8 periodic orbits, and all of the cells constitute the domain of attraction of P-8 periodic orbits.Figure 7
Domain of attraction of P-8 periodic orbit.The degree of stability of P-8 periodic orbits can be measured by the scope of domain of attraction; namely, a bigger scope of domain of attraction shows that the corresponding periodic orbits are more stable. If the disturbance of P-8 periodic orbit is limited in its domain of attraction, the motion will keep stable.There are 5470 cells in the domain of attraction of P-8 periodic orbits showed in Figure6, and 6.08% of the whole cell space in D belonged to this domain of attraction. The domain of attraction showed in Figure 7 is composed of three subspaces. Periodic cells of P-8 are full of the three subspaces, but it is not continuous among each of the subspaces.Some conclusions about stability of P-8 periodic orbit can be gotten from the above discussions. Because the domain of attraction is full of the occupied cell spaces, P-8 periodic orbit is of local stability; namely, if disturbance of this orbit is limited in its domain of attraction, the orbit will keep stable, and if disturbance increases to be out of its domain of attraction, the orbit will lose its stability and turn into other coexisting periodic orbits.Figure8 is domain of attraction of P-4 periodic orbit. Except for the cell space which is occupied by domain of attraction of P-8, the domain of attraction of P-4 is almost full of the whole cell space. There are 84530 cells attracted into P-4 periodic orbit; namely, 93.92% of cells in the whole cell space belong to domain of attraction of P-4 periodic orbit. It is obvious that P-4 periodic orbit is of global stability.Figure 8
The domain of attraction of P-4 motion.


## 5. Conclusions

(1)
By using Poincare-like cell-to-cell method and shooting method, the global characteristics of a planetary gear set are studied based on a nonlinear torsional dynamical model with gear backlashes and time-varying mesh stiffness.
(2)
With the given set of parameters, planetary gear set has two coexisting periodic orbits with long-term stability; one is P-4 periodic motion, and another is P-8 periodic motion. The P-4 periodic motion is of global stability and the P-8 periodic motion is of local stability.
(3)
The shooting method is not of the capacity of searching coexisting periodic orbits in a global scope. Namely, the shooting method can only find some coexisting periodic orbits which are close to each other and easy to omit some periodic orbits which are far away from main periodic orbit groups.



---

*Source: 930350-2014-02-16.xml*