 Open Access
 Total Downloads : 1872
 Authors : Majid Mehrabi, Drv.P.Singh
 Paper ID : IJERTV2IS90736
 Volume & Issue : Volume 02, Issue 09 (September 2013)
 Published (First Online): 25092013
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Vibration Analysis of Planetary Gear System
Majid Mehrabi (research scholar)
Mechanical Engg. DepartmentPEC University of Technology Chandigarh India Department of Mechanical engineering Islamic Azad UniversityTakestan Branch DrV.P.singh (Associate professor)
Mechanical Engg Department PEC University of Technology Chandigarh India
Abstract
In this work a dynamic model of a planetary gear transmission is developed to study the sensitivity of the natural frequencies and vibration modes to system parameters in perturbed situation. Parameters under consideration include component masses ,moment of inertia , mesh and support stiff nesses .The model admits three planar degree of freedom for planets ,sun, ring, and carrier. Vibration modes are classified into translational, rotational and planet modes .Welldefined modal properties of tuned (cyclically symmetric) planetary gears for linear ,timeinvariant case are used to calculate eigensensitivities and express them in simple formulae .These formulae provide efficient mean to determine the sensitivity to stiffness ,mass and inertia parameters in perturbed situation
1. Introduction
Planetary gearboxes are usually used in a wide variety of machinery such as automobiles, helicopters and aircraft engines .Their numerous advantages are high speed reductions in compact spaces, high torque/weight ratio , greater load sharing , diminished bearing loads and reduced noise and vibration. A typical simple planetary gear set consists of a sun gear, a ring gear and a number of identical planet gears (typically 36) meshing both with the sun and ring gears. They are well known for their symmetrical structure which allows an equal share of the total external torque applied between the planetary gears, the sun and the ring. However non stationary conditions of system such as overload
conditions, torque fluctuation may affect the dynamic behavior of a planetary gear transmission. The inequality of the load distribution, however arises on each planet gear because of random errors of manufacture , assembly and operating conditions. This results in noise and vibration which are key concerns in their applications and drop in efficiency of planetary gear system. In some helicopters planetary gear vibration is the primary source of cabin noise that can exceed 100dB.
Before 1990, the literature on analytical planetary gear dynamics is scarce. Cunliffe et al.[1] studied the eigenvalue problem for a thirteen degree of freedom system and identified the natural frequencies and vibration modes. Botman[2] investigated the vibration modes for a spur planetary gear with eighteen degrees of freedom. The effect of carrier rotation and variation of planet bearing stiffness on the natural frequencies are studied numerically. Frater et al.[3] extended Botman's natural frequency analysis by including unequal planet stiffness. Kahraman[4,5,6,7] presented a series of papers on planetary gear dynamics. He derived a nonlinear, timevarying planar dynamic model9 and subsequently extended it to threedimensions and examined the influence of planet phasing on dynamic response[5]. Kahraman[7] reduced his model to a purely torsional one to predict natural frequencies and vibration modes. Parker[8] rigorously proved the
effectiveness of using planet phasing schemes to suppress planetary gear vibration. Kahraman [6] showed the effects of mesh/support stiffnesses on the natural frequencies in his torsional model of planetary gears. Investigation of natural frequency
and vibration mode sensitivity to system parameters for tuned and mistuned planetary gears was done by J.Lin and R.G.Parker [15].They also investigated the natural frequency spectra and vibration modes of planetary gears to avoid power train resonances.[7 ]
Analytical investigation of the sensitivities of natural frequencies and vibration modes to stiffness and inertia parameters of general compound planetary gears for both tuned and mistuned are studied in[12]. How the highly structured freevibration properties of equally spaced planet systems change due to unequal planet spacing was studied in[5].Characterizing the effects of various errors on the dynamic properties of planetary gear system including tooth thickness
,runout and load sharing as the most common gear system manufacturing was done by Parker,R.G. and G.J.Cheon.[13].The effects of some important parameters ,such as the variation of mesh stiffness and static transmission errors on the nonlinear dynamics of a planetary gear system with multiple clearances was studied in[14].The objective of this paper is to analytically investigate the natural frequency and vibration modes of general planetary gears in perturbed situation and comparing them with unperturbed situation .This allows one to find the dominant parameter effecting the perturbation of the system .Investigating the natural frequency spectra one can tune the system to avoid resonance .
2 Modeling and equations
The analysis deals with planar vibration of single stage planetary gears .A Lumpedparameter model used in this work for dynamic analysis is shown in Fig 1. All gears are considered as rigid bodies and component supports are modeled by springs. Each of sun, ring, carrier and Z planets are treated as rigid bodies. Each component has three degrees of freedom: two translational and one rotational. The model is similar to that used by Parker [12], the planet deflections are described by radial and tangential coordinates and more naturally describe the vibration modes. The coordinates illustrated in Fig.1 are used. The carrier, ring and sun translations , , where g=c,r,s and planet translations , , L = 1,…Z are measured with
respect to a rotating frame fixed to the carrier with
origin O. The are directed towards the equilibrium position of planet 1, and , are the radial and tangential deflections of the zth planet. The rotating frame rotates with the constant carrier angular speed c. The rotational coordinates are =
,u = c,r,s,1,… Z, where is the component rotation; is the base circle radius for the sun, ring and planet, and the radius of the circle passing
through the planet centers for the carrier. Circumferential planet locations are specified by the fixed angles , where is measured relative to the
rotating basis vector i so that 1 = 0. Fig.1 shows a sun planet mesh with masses ms, mp and moment of inertia, . s is the pressure angle of sun/planet mesh and = s .
Figure 1. Lumped parameter model of planetary gears and system coordinates. (b) All translational coordinates ag , bg , g =c, r, s and , , L = 1,…z are with respect to the frame {i, j, k} rotating at constant carrier speed c .
Free vibration ,timeinvariant representation is considered and assumed identical. All planet bearing stiff nesses are equal kpn=kp, all sunplanet mesh stiffnesses are equal ksn= ksp and all ring planet mesh stiff nesses are equal krn=krp .The planets are equally spaced and cyclically symmetric structures, therefore can be divided into Z sector having a central angle = 2. Following the matrix derivation of Fox
Z
method [16, 17], the general eigenvalue equation for
perturbed free vibration is: [19] [M][X][] = 0 (1)
de = perturbation in design variable e
[I] = identity matrix [K1 ] = unperturbed stiffness matrix [dKe ] = perturbation in [K1 ] due to change in design variable e. [K] = perturbed stiffness matrix1 = unperturbed eigenvalue matrix (dagonal)
[de ] = perturbation in 1 due to change in design variable e [M1] = unperturbed mass matrix [dMe ] = perturbation in [M1] due to change in design variable e [M] = perturbed mass matrix[X1 ] = Unperturbed eigenvector matrix (mode shape)
Where is the spring deformation of the sunplanet z mesh in mode 1 given by
dXe = perturbation in [X1 ] due to change in design
= cos
sin
variable e
s
s
s
[X] = perturbed eigenvector matrixs + +
The rotational mode property implies that all sun
M is the inertia matrix and K
is the bearing stiffness
planet mesh deformations are equal i.e. = , so
b
matrix given in Appendix A.. To model the time varying stiffness associated with changing numbers of teeth in contact at each mesh, Km can be
decomposed into mean and timevarying
equation (3) becomes
= ( )2
components. Tooth separation nonlinearity is implicitly included in Km(t). Manipulating equation
1
1
(1) we obtain the first order equation
+ [ ] = [ ] (2)
The Eigen sensitivity to masses (ms , mr, mc , mp ) and
moment of inertia (Is, Ir, Ic, Ip )for sun ,ring ,carrier
1 1 and planets for a tuned perturbed system are
The diagonal terms give the eigenvalue perturbations
considered ,the eigenvalue derivatives for the three types of modes are
= [ ]
2 2 2
= + =
Or the derivatives
= 2 = 2
g = r ,s ,c
e
de e
2
i =
i = ii
Z Z
2
e de de
= (2 + 2) =
Considering the case of equal alteration of all sun planet mesh stiffness i.e., ksz = ksp , the perturbed system remains tuned. Eigenvector derivative can be
Z
L=1
L=1
written as [18]
= 2
dX = (In X1 eT ). 1 In M1 . (M1
2
L=1
j 1 1
T
= modal translational kinetic energy of ring, sun
X1 X1
)(dK . dMe)X
and carrier
1
1
XT M1X1 1 1
Considering the case of equal alteration of all sun planet mesh stiffness i.e., = ,the perturbed system remains tuned. The derivatives of the mass and stiffness matrices with respect to are:[5,12]
= modal rotational kinetic energy of ring, sun and carrier
Planet Mode
Characteristics of planet mode in perturbed situation are:
All sub matrices of
are zero except the four that

Carrier,ring and sun have zero translation and rotation i.e. = [0,0,0]

Natural frequencies has multiplicity Z =1 in perturbed situation.
Therefore planet mode will get the form
involve . For rotational mode ,the Eigen sensitivies are obtained as below
= 0,0,0, 1, 11 , 21, , 1
is a multiple of 1 = 1, then
= ( )2
(3)
z sinz = 0 z cosz = 0 z =0 (4)
1
In planet modes the sun ,carrier and ring have zero rotation and translation.The Eigen solution properties
=
is illustrated using the nominal system parameters
=1, =1
are listed in table 1.
Table1 Nominal system parameters of the planetary gear
Sun 
Ring 
Carrier 
Planet 

Mass (kg) 
0.4 
2.35 
5.43 
0.66 
I/r2 (kg) 
0.39 
3.00 
6.29 
0.61 
Base diameter (mm) 
77.42 
275.03 
177.8 
100.35 
Teeth number 
27 
99 
35 

Mesh stiffness (N/m) 
ksp = krp = km = 5×108 

Bearing stiffness (N/m) 
kp = ks = kr =kc = 108 

Torsional stiffness (N/m) 
kru = 109 ksu = kcu = 0 

Pressure angle (Â°) 
s = r = = 24.6 
The variation of eigenvalues in perturbed situation for two different conditions are shown below:
1 is assumed constant , K is assumed variable
.(Figure 2)
Figure – 2
2 M is assumed variable , K is assumed constant .(Figure 3)
Typical vibration modes for equally spaced planets are shown below. The movements of the ring and carrier are not shown in figure 4 .Solid lines are perturbed situation and dotted lines are unperturbed situation.
Conclusion and future work
The present study identifies the parameters affecting the natural frequencies of a general planetary gear system. The results can be applied to specific configurations (fixed sun, carrier or ring).For a tuned system a welldefined structure to the natural frequencies are obtained. The results are:
Sensitivity to mass, moment of inertia and stiffness in perturbed situation is shown in Figure2,3. It is observed that fluctuation in all the natural frequencies presented in figure 2 are on higher side than those in figure 3.It is concluded that stiffness component dominates prominently than mass on the natural frequency From figures 2 and 3 it is observed that in perturbed tuned situation with degenerate natural frequency multiplicity m is 1.Figure 4 shows typical vibration modes for equally spaced planets. For planet modes ,motion of ring ,carrier and sun is zero. In gear systems satisfying relation (4),translation and rotation modes have structured properties. Planet motions are a multiple of the motion of first planet. For sun, carrier and ring the rotational modes have pure rotation. All planets move in identical form. Their translational modes have pure translation.For vibration modes reducedorder Eigen value problems are achieved. Sensitivity to mass, moment of inertia and stiffness are investigated for perturbed situation .
The scope for future work is to apply the method used in this paper for various mass components and stiffnesses using relevant materials for manufacturing gears and propose the one having less influence on the vibration characteristic of the system.
Reference

Cunliffe, F., Smith, J.D. and Welboum, D.B., 1974, Dynamic Tooth Loads in Epicyclic Gears, Journal of Engineering for Industry, 94: 578584.

Botman, M., 1976, Epicyclic Gear Vibrations, Journal of Engineering for Industry, 96:811 815.

Frater, J.L., August, R. and Oswald, F.B., 1983, Vibration in Planetary Gear Systems with unequal Planet Stiffnesses, NASA Technical Memorandum 83428 .

Kahraman, A., 1994a, Load Sharing Characteristics of Planetary Transmissions, Mechanism and Machine Theory, 29:11511165.

Kahraman, A., 1994b, Planetary GearTrain Dynamics, ASME Journal of Mechanical Design, 116(3): 713720.
6 Kahraman, A., 1994c, Natural Modes of Planetary Gear Trains, Journal of Sound and Vibration, 173(1): 125130.

Kahraman, A., 1994d, Dynamic Analysis of a MultiMesh Helical Gear Train, ASME Journal of Mechanical Design, 116(3): 706
712.

Parker, R.G., 2000, A Physical Explanation for the Effectiveness of Planet Phasing to Suppress Planetary Gear Vibration, Journal of Sound and Vibration, 236(4): 561573.

Lin, J. and Parker, R.G., 2000a, Structured Vibration Characteristics of Planetary Gears with Unequally Spaced Planets, Journal of Sound and Vibration, 233(5): 921928.

Yichao,G and Parker,R.G.2010, Sensitivity of General Compound Planetary Gear Natural Frequencies and Vibration Modes to Modal Parameters, Journal of ASME
,vol.132/0110061

Kahraman, A., 1994a, Load Sharing Characteristics of Planetary Transmissions, Mechanism and Machine Theory, 29:1151 1165.

12 Lin, J. and Parker, R.G., 1999a, Analytical Characterization of the Unique Properties of Planetary Gear Free Vibration, Journal of Vibration and Acoustics, 121 : 316321.

GillJeon Cheon and Parker, R.G., 2004, Influence of Manufacturing Errors on the Dynamic Characteristic of Planetary Gear Systems , KSME International Journal,vol 18 No.4, pp.606621.

Tao Sun, HaiYan Hu ., 2003 Nonlinear Dynamics of a Planetary gear system with Multiple Clearances , Mechanism and Machine theory 38, pp. 13711390 .

Lin, J. and Parker, R.G., 1999b, Sensitivity of Planetary Gear Natural Frequencies andVibration Modes to Model Parameters, Journal of Sound and Vibration, 228(1 ): 109128.

Herting,D.N., Eigenvector Derivative with second order Terms,MSC Memo DNH43,july 6,1984.

Gerald D.High, 1990, An Iterative Method for Eigenvector Derivatives , The MacNeal Schwendler Corporation, Los Angeles ,CA.

Michele Bernasconi,Christine Choirat and Raffaello Seri, 2011, Differentials of Eigenvalues and Eigenvectors in undamped Discrete Systems under
Alternative Normalizations , proceedings of the World Congress on Engineering vol.

Majid Mehrabi , Raveesh Pratap , Analysis of Vibration Characteristics of Planetary GearSystem, 4th International Conference on Emerging Trends in Engineering and
sin sins sin coss sin
K z = k sin cos cos cos
K z = k sin cos cos cos
r2 rz cos s s
sins coss 1
sin2r cosrsinr sinr
Technology,Kurukshetra, India
APPENDIX A
M=diag(Mc,Mr,Ms,M1,…,MZ)
Kr z = krz
3
3
K z = k
cos2r cosr
symmetric 1
sin2 cos sin sin
cos2 cos
M = diag(m ,m ,I /r 2), L=c,r,s,1,,Z
s1 sz
g h h h h
symmetric 1
Kb=diag(Kcb,Krb,Ksb,0,,0)
Kgb = diag(kgx,kgy,kgu) g=c,r,s
sin sins sin coss sin
K z=k cos sins cos coss cos
Kc1Z 0 0 Kc21 Kc22 . Kc2z
s2 sz
sins coss 1
Kr1Z 0 Kr21 Kr22 . Kr2z
Ks1Z
Ks21
Ks22
. Ks2z
sin2s cosssins sins
K z=k cos2s coss
Km =
Kpp1 0 . 0
.
s3 sz
symmetric 1
symmetric
Kppz
= s = + s
z = K + K z + K
z = K + K z + K
Kpp c3z r3 s3z
1 0 sin
z 1 cos
Kc1 = kpz
symm. 1
cos sin 0
Kc2z = kpz
sin
cos 0
0 1 0
Kc3 pz pz
Kc3 pz pz
z = diag(k ,k ,0)
sin2 cos sin sin
K z = k cos2 cos
r1 rz
symmetric 1