# Dynamic Characteristics of Encased Differential Gear Train with Journal Bearing

Authors: Jie Yang; Yanjiong Yue; Rupeng Zhu; Weifang Chen; Miaomiao Li
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

Taking the marine encased differential gear train as an example, the relationship between the journal bearing parameters and the meshing force of the transmission system is analyzed. In this paper, the dynamic model of the encased differential gear train with journal bearing is established considering the factors of time-varying meshing stiffness and comprehensive transmission error. In this dynamic model, four stiffnesses and four damping coefficients are applied to characterize the asymmetry and interaction of the oil film stiffness and damping of planet bearing. The system responses are calculated by the Fourier series numerical algorithm. The results show that the introduction of journal bearing in encased differential gear train can contribute to gearbox vibration reduction. Moreover, the planet bearing parameters (e.g., clearance-to-radius ratio and eccentricity ratio) of the differential stage affect the meshing forces of both the differential and encased stages. In addition, the influence of the planet bearing parameters of the encased stage on the meshing force of the encased stage is more obvious than that of the differential stage. This work may develop a theoretical analysis framework for the design and manufacture of marine transmission systems in the future.


---

## Body

## 1. Introduction

The compound planetary gear train, which performs well in compactness, high power density, and large torque-to-weight ratio, is widely used in ship, helicopter, wind turbine, etc. As for heavy machinery transmission systems, journal bearings typically support the carrier and planet in compound planetary gear train, and their stiffness and damping are unequal in all directions [1]. Vibration noise is an important factor affecting the stability and reliability of the gearbox. Studying the dynamic characteristics of a compound planetary gear train may help to reduce the vibration noise.Some researchers [2–5] focused on the construction of the gear-rotor dynamics model. Using the finite element method, Baguet and Jacquenot [2] established a gear-rotor-bearing coupling model, in which the time-varying properties and nonlinearities of gears and bearings were introduced. Gao and Zhang [3] investigated the gear-bearing dynamic interactions using a dynamic model of a rotor-bearing system with confluence transmission. Theodossiades and Natsiavas [4] studied the dynamic response of the gear-rotor coupled system under the action of oil film force and meshing force and not only investigated the oil film stability but also analyzed the amplitude-frequency response. Based on nonlinear theory, Wan and Chang [5] carried out a systematic analysis of the dynamic behaviors of a gear-rotor-bearing with strongly nonlinear.Some scholars have thrown more lights to study the influence of backlash, time-varying meshing stiffness, and manufacturing error on the dynamic characteristics of planetary gear train. Based on elastodynamics, Guo and Parker [6] established a pure torsional dynamic model of the planetary gear train and analyzed the natural frequency and parametric sensitivity. Similarly, Zhang et al. [7] analyzed the characteristics of natural frequencies and mode shapes of an encased differential gear train. Wu et al. [8] studied the influence of backlash on the nonlinear characteristics of compound planetary gear set by analyzing the change of backlash with gear wear. Zhu et al. [9] explored the influence of multibacklash on bilateral shock of encased differential planetary gear train. Lin and Parker [10] used a multiscale method to investigate the influence of meshing stiffness waves on the stability of the planetary gear train, and they derived the analytical formula of the relationship between meshing stiffness variation and system dynamic stability. Al-Shyyab et al. [11] constructed a torsional dynamics model of a discrete two-stage planetary transmission system, and they studied the frequency response of the system under stiffness excitation. Considering time-varying meshing stiffness, Xiang et al. [12] developed a torsional nonlinear dynamic model of multistage planetary gear trains and studied the effects of bifurcation parameters such as excitation frequency, backlash, and damping.In addition, Iglesias et al. [13] investigated the effects of eccentricity and planet pin positioning errors on the quasi-static behavior of planetary gear train for different configurations of support parameters. In terms of multiple errors (e.g., eccentricity error, base pitch error, and assembly error), Mo et al. [14] proposed an improved mathematical model of the loading sharing coefficient for compound planetary gear train, which was verified by experiments. Using a lumped parameter method, Zhai et al. [15] developed the bending-torsional-pendular coupled vibration model of the wind turbine gearbox transmission system and studied the frequency response of carrier assembly errors (e.g., carrier horizontal assembly error, carrier axis assembly error, and carrier spatial assembly error) to the dynamic meshing force and dynamic transmission error.In most of the abovementioned studies, the planet bearing assumes a fixed stiffness spring and the bearing force is assumed to be linear with displacement. However, the asymmetry and interaction of oil film stiffness and damping of planet bearing is neglected. To fill this gap, we use four stiffnesses and four damping coefficients to explore the dynamic characteristics of encased differential gear train with journal bearing.In this study, the marine encased differential gear train with journal bearing is taken as the research object. First, the lumped parameter method is utilized to establish a dynamic model of the encased differential gear train considering time-varying meshing stiffness and comprehensive transmission errors. Second, the asymmetry and interaction of the oil film stiffness and damping of journal bearing is employed to the dynamic model, and the gears and carrier are treated according to the concentrated mass. Third, the time-varying meshing stiffness of herringbone gears is calculated based on the time-varying total contact line length [16]. Fourth, the dynamic equations of coupling journal bearings are constructed by Newton’s second law. The dynamic response of encased differential gear train with journal bearing are solved by the numerical algorithm of Fourier series. Finally, the effects of parameters of planet bearing on the meshing forces are explored in detail.


## 2. Dynamic Model of Compound Planetary Gear Train

The encased differential gear train consisting of a single-degree-of-freedom encased planetary gear train (encased stage) and a two-degree-of-freedom differential planetary gear train (differential stage) is shown in Figure1. In the differential stage, No. 1, No. 2, No. 3, No. 4, No. 5, No. 6, and No. 7, respectively, represent the floating ring gear Zf1, the intermediate floating member Zg1, the composting ring gear Zr1, the planet Zpi (i = 1, 2, …, N, where N is the number of planets), the sun gear Zs1, the carrier H, and planet bearing. In the encased stage, No. 8, No. 9, No. 10, No. 11, No. 12, and No. 13, respectively, represent the floating ring gear Zf2, the intermediate floating memberZg2, the composting ring gear Zr2, the planet Zmj (j = 1, 2, …, M, where M is the number of planets), and the sun gear Zs2, and planet bearing. No. 14 and No. 15 represent journal bearings at input and output shafts, respectively.Figure 1
Diagram of the transmission system.Herringbone gears are employed in the sun gearsZs1/Zs2 (i.e., No. 5/No. 12) and planet Zpi/Zm (i.e., No. 4/No. 11). The composting ring gear Zr1/Zr2 (i.e., No. 3/No. 10) is composed of internal and external helical gears on the same gear tooth germ. The intermediate floating member Zg1/Zg2 (i.e., No. 2/No. 9) is composed of internal helical and external spur gears on the same gear tooth germ. The internal helical gear of Zr1/Zr2 meshes with the unilateral helical gear of planet Zpi/Zmj, and the external helical gear of Zr1/Zr2 meshes with the internal helical gear of Zg1/Zg2. The external spur gear of Zg1/Zg2 meshes with the internal spur gear of Zf1/Zf2 (i.e., No. 1/No. 8).Figure2 shows the dynamic model of the encased differential gear train. The rotational coordinate system based on the rotation angular velocity of carrier H is applied to the differential stage. The fixed coordinate system is used for the encased stage. The degree-of-freedom of the encased differential gear train is 19 + 3(N + M). Displacement parameter X is defined by(1)X=xs1,Hs1,Vs1,xpi,Hpi,Vpi,xr1,Hr1,Vr1,xg1,xf1,xs2,Hs2,Vs2,xmj,Hmj,Vmj,xr2,Hr2,Vr2,xg2,xf2,xH,HH,VH,where xs1, xpi, xr1, xs2, xmj, and xr2 are torsional line displacement parameters along base circle radius of Zs1, Zpi, Zr1, Zs2, Zmj, and Zr2, respectively. xg1, xf1, xg2, and xf2 are torsional line displacement parameters along the pitch circle radius of Zf1, Zg1, Zf2, and Zg2, respectively. xH, HH, and VH are torsional, lateral, and longitudinal displacements of carrier H, respectively. Hs1, Hpi, and Hr1 are lateral center displacements of Zs1, Zpi, and Zr1, respectively, in their rotational coordinates. Vsi, Vpi, and Vr1 are longitudinal center displacements of Zs1, Zpi, and Zr1, respectively, in their rotational coordinates. Hs2, Hmj, and Hr2 are lateral center displacements of Zs2, Zmj, and Zr2, respectively, in fixed coordinate. Vs2, Vmj, and Vr2 are longitudinal center displacements of Zs2, Zmj, and Zr2, respectively, in fixed coordinate.Figure 2
Dynamic model of the encased differential gear train.


### 2.1. Model for Calculation of Journal Bearing Stiffness and Damping

The dynamic model of the journal bearing is shown in Figure3. Assuming that displacement and velocity remain small disturbances in the vicinity of the static position O0 (x0, y0) [1], the first-order expansion of dynamic bearing forces (i.e., Fx and Fy) is defined bywhere ∂Fi/∂j∣0i=x,y;j=x,y,x˙,y˙ is the partial derivative at static position O0 (x0, y0).(2)
F
x
≈
F
x
0
+
∂
F
x
∂
x
0
x
−
x
0
+
∂
F
x
∂
y
0
y
−
y
0
+
∂
F
x
∂
x
0
x
˙
+
∂
F
x
∂
x
0
y
˙
,
F
y
≈
F
y
0
+
∂
F
y
∂
x
0
x
−
x
0
+
∂
F
y
∂
y
0
y
−
y
0
+
∂
F
y
∂
x
0
x
˙
+
∂
F
y
∂
y
0
y
˙
,Figure 3
Dynamic model of journal bearing.Here,(3)kij=∂Fi∂j0i=x,y;j=x,y,dij=∂Fi∂j0i=x,y;j=x˙,y˙,ΔFx=Fx−Fx0=kxxξ+kxyη+dxxξ˙+dxyη˙,ΔFy=Fy−Fy0=kyxξ+kyyη+dyxξ˙+dyyη˙,where Fx0 and Fy0 are static bearing forces in x and y directions, respectively. ΔFx and ΔFy are the increments of the oil film force. ξ and η are the increments of oil-film thickness in x and y directions, respectively. kij and dij (i, j = x, y) are the stiffness and damping coefficients of journal bearing, which are obtained by(4)kij=KijμωBψ3,dij=DijμBψ3i,j=x,y,Kxx=Bd24ε2π2+16−π2ε21−ε2216ε2+π21−ε2,Dxx=Bd22ππ2+2π2ε2−16ε21−ε2216ε2+π21−ε2,Kxy=Bd2π−π2+2π2ε2+16−π2ε41−ε25/216ε2+π21−ε2,Dxy=Bd28επ2+12π2ε2−16ε21−ε2216ε2+π21−ε2,Kyx=Bd2ππ2+π2+32ε2+216−π2ε41−ε25/216ε2+π21−ε2,Dyx=Dxy,Kyy=Bd24επ2+π2+32ε2+216−π2ε41−ε2216ε2+π21−ε2,Dyy=Bd22ππ2+12π2ε2−16ε21−ε2216ε2+π21−ε2,where B and d are the width and diameter of the journal bearing, μ is the oil viscosity, ε is the eccentricity ratio, ψ is the clearance-to-radius ratio, and Kij and Dij (i, j = x, y) are the corresponding dimensionless stiffness and damping coefficient of journal bearing.


### 2.2. Model for Calculation of Time-Varying Meshing Stiffness

The planets are evenly distributed among the sun gear in the planetary gear train. There are varying mesh phasing factors in meshing pairs (i.e., sun-planet and ring-planet) because of the influence of the number of teeth. The mesh phasing factor affects the phase of the time-varying meshing stiffness [17]. Therefore, the mesh phasing factor is incorporated into the encased differential gear train. Following Parker and Lin [18], the mesh phasing factor is defined by(5)γsm=Zsψm2π,γrm=−Zrψm2π,where γsm and γrm are the mesh phasing factor of sun-planet and ring-planet, respectively, and Zs and Zr are the teeth number of the sun and ring, respectively.Maatar and Velex [16] proposed the time-varying meshing stiffness of helical gear pairs based on the time-varying total contact line length. Following Maatar and Velex [16], the time-varying meshing stiffness model of the helical gear is applied to the herringbone gear by a parallel method in this study. The time-varying meshing stiffness model of the herringbone gear is constructed as follows:(6)kt=2k0Lτ,where k0 is the average value of the time-varying meshing stiffness, τ is dimensionless time, τ = t/Tm, t and Tm are meshing time and period, and L(τ) is the instantaneous total contact length of gear pair, which is obtained by(7)Lτ=1+∑k=1SAkcos2πkτ−γ+Bksin2πkτ−γbεαcosβb,where γ is the mesh phasing factor, εα is the end-face contact ratio of the herringbone gear, βb is the helical angle of the base circle of the herringbone gear, B is the unilateral tooth width of the herringbone gear, Ak and Bk are the kth-order harmonic amplitude of the time-varying meshing stiffness, S represents the order of the Fourier series, and the series is truncated at S = 5 in this article.The parametersAk and Bk are calculated by(8)Ak=12εαεβπ2k2cos2πkεβ+cos2πkεα−cos2πkεα+εβ−1,Bk=12εαεβπ2k2sin2πkεβ+sin2πkεα−sin2πkεα+εβ,where εβ is axial contact ratio of the herringbone gear.


### 2.3. The Equivalent Displacement and Meshing Forces

The center displacement of each gear is projected into the meshing line direction. We ignored the change of the meshing angle caused by the movement of each gear center. Figure4 shows the positional relationship sun gear, planet, and composting ring gear in the differential stage.Figure 4
Position relationship of each component in the differential stage.The equivalent displacement of theith planet-composting ring and sun-planet along meshing line direction can be calculated by(9)Srpit=xpi−xr1+Hpisinαrp−Vpicosαrp,−Hr1sinαrp+ψi+Vr1cosαrp+ψi−erpit,Sspit=xs1−xpi+Hs1sinαsp−ψi+Vs1cosαsp−ψi,−Hpisinαsp−Vpicosαsp−espit,where erpi and espi are the equivalent error of the ith planet-composting ring and sun-planet along meshing line direction.According to the equivalent displacement calculation method of the differential stage, the equivalent displacement of thejth sun-planet and composting ring-planet along meshing line direction can be calculated by(10)Srmjt=xmj−xr2+Hmjsinαrm−Vmjcosαrm,−Hr2sinαrm+ψj+Vr2cosαrm+ψj−ermjt,Ssmjt=xs2−xmj+Hs2sinαsm−ψj+Vs2cosαsm−ψj,−Hmjsinαsm−Vmjcosαsm−esmjt.The planet is supported by planet bearing in the encased differential gear train, and thus the compression deformation between carrier and planet bearing can be obtained by(11)δHcpi=Hccosψi+Vcsinψi−Hpi,δVcpi=−Hcsinψi+Vccosψi−Vpi−xc,where δHcpi is the compression deformation of planet bearing along the lateral direction and δVcpi is the compression deformation along the longitudinal direction.Elastic and damping meshing forces in the differential and encased stages can be calculated by(12)Fit=KitSitDit=CiS˙it,i=spi,rpi,smj,andrmj,where S˙it is the first derivative of equivalent displacement. Ci is the meshing damping factor, which can be calculated by [19](13)Ci=2ξik0,i1/mpi+1/mgi,i=spi,rpi,smj,andrmj,where ζi denotes the constant damping ratios of the gear pair, ζi = 0.07; k0,i represents the average value of the meshing stiffness in the gear pair, and mpi and mgi are the masses of the driving and driven gears in the gear pair.


### 2.4. Differential Equations of the Compound Planetary Gear Trains

According to Figure2 of the dynamic model of the encased differential gear train, differential equations of 19 + 3 (N + M) degree-of-freedom are constructed for the encased differential gear train by the lumped parameter method. Differential equations of the differential stage can be deduced using Newton’s second law as follows.(1)
The sun gearZs1:(14)ms1x¨s1+∑i=1NFspit+Dspit=PD,Ms1H¨s1+∑i=1NFspit+Dspitsinαsp−ψi+Ks1Hs1=−Ms1gsinωHt,Ms1V¨s1+∑i=1NFspit+Dspitcosαsp−ψi+Ks1Vs1=−Ms1gcosωHt.(2)
Theith planet gear Zpi:(15)mpx¨pi−Fspit+Dspit+Frpit+Drpit=0,MpH¨pi−Fspit+Dspitsinαsp+Frpit+Drpitsinαrp+FHcpit=−MpigsinωHt+ψi+Fpi,MpV¨pi−Fspit+Dspitcosαsp+Frpit+Drpitcosαrp+FVcpit=−MpigcosωHt+ψi,whereFHcpi(t) is the oil film force of the ith planet bearing in the lateral direction, FVcpit is the oil film force of the ith planet bearing in the longitudinal direction, FHcpi(t) and FVcpit are calculated by equation (3), and Fpi is the centrifugal force of the ith planet Zpi.(3)
The composting ring gearZr1:(16)mr1x¨r1−∑i=1NFrpit+Drpit+kg1qxr1−xg1=0,Mr1H¨r1−∑i=1NFrpit+Drpitsinαrp+ψi+Kr1Hr1=−Mr1gsinωHt,Mr1V¨r1+∑i=1NFrpit+Drpitcosαrp+ψi+Kr1Vr1=−Mr1gcosωHt.(4)
The intermediate floating memberZg1 and floating ring Zf1:(17)mg1x¨g1−Kg1qxr1−xg1+Kf1qxg1−xf1=0,mf1x¨f1−Kf1qxg1−xf1+Kf1s2rf1bxf1rf1b−xs2rs2b=0.Differential equations of the encased stage can be deduced using Newton’s second law as follows:(1)
The sun gearZs2:(18)ms2x¨s2−Kf1s2rs2bxf1rf1b−xs2rs2b+∑j=1MFsmjt+Dsmjt=0,Ms2H¨s2+∑j=1MFsmjt+Dsmjtsinαsm−ψj+Ks2Hs2=0,Ms2V¨s2+∑j=1MFsmjt+Dsmjtcosαsm−ψj+Ks2Vs2=−Ms2g.(2)
Thejth planet gear Zmj:(19)mmx¨mj−Fsmjt+Dsmjt+Frmjt+Drmjt=0,MmH¨mj−Fsmjt+Dsmjtsinαsm+Frmjt+Drmjtsinαrm+FHmjt=−Mmgsinψj,MmV¨mj−Fsmjt+Dsmjtcosαsm+Frmjt+Drmjtcosαrm+FVmjt=−Mmgcosψj,whereFHmj(t) is the oil film force of the jth planet bearing in the lateral direction and FVmjt is the oil film force of the jth planet bearing in the longitudinal direction.(3)
The composting ring gearZr2:(20)mr2x¨r2−∑j=1MFrmjt+Drmjt+kg2qxr2−xg2=0,Mr2H¨r2−∑j=1MFrmjt+Drmjtsinαrm+ψj+Kr2Hr2=0,Mr2V¨r2+∑i=1MFrmjt+Drmjtcosαrm+ψj+Kr2Vr2=−Mr2g.(4)
The intermediate floating memberZg2 and floating ring Zf2:(21)mg2x¨g2−Kg2qxr2−xg2+Kf2qxg2−xf2=0,mf2x¨f2−Kf2qxg2−xf2+Kf2Lrf2bxf2rf2b−xcrc=0.The differential equation of carrierH can be deduced by Newton’s law:(22)mcx¨c+∑i=1NFVcpi+KcHrcxcrc−Kf2Hrcxf2rf2b−xcrc=−TLrc,McH¨c+∑i=1NFHcpicosψi−FVcpisinψi+FcH=0,McV¨c+∑i=1NFHcpisinψi+FVcpicosψi+FcV=−Mcg,where FcH is the oil film force of journal bearing in the lateral direction and FcVis the oil film force of journal bearing in the longitudinal direction.


## 2.1. Model for Calculation of Journal Bearing Stiffness and Damping

The dynamic model of the journal bearing is shown in Figure3. Assuming that displacement and velocity remain small disturbances in the vicinity of the static position O0 (x0, y0) [1], the first-order expansion of dynamic bearing forces (i.e., Fx and Fy) is defined bywhere ∂Fi/∂j∣0i=x,y;j=x,y,x˙,y˙ is the partial derivative at static position O0 (x0, y0).(2)
F
x
≈
F
x
0
+
∂
F
x
∂
x
0
x
−
x
0
+
∂
F
x
∂
y
0
y
−
y
0
+
∂
F
x
∂
x
0
x
˙
+
∂
F
x
∂
x
0
y
˙
,
F
y
≈
F
y
0
+
∂
F
y
∂
x
0
x
−
x
0
+
∂
F
y
∂
y
0
y
−
y
0
+
∂
F
y
∂
x
0
x
˙
+
∂
F
y
∂
y
0
y
˙
,Figure 3
Dynamic model of journal bearing.Here,(3)kij=∂Fi∂j0i=x,y;j=x,y,dij=∂Fi∂j0i=x,y;j=x˙,y˙,ΔFx=Fx−Fx0=kxxξ+kxyη+dxxξ˙+dxyη˙,ΔFy=Fy−Fy0=kyxξ+kyyη+dyxξ˙+dyyη˙,where Fx0 and Fy0 are static bearing forces in x and y directions, respectively. ΔFx and ΔFy are the increments of the oil film force. ξ and η are the increments of oil-film thickness in x and y directions, respectively. kij and dij (i, j = x, y) are the stiffness and damping coefficients of journal bearing, which are obtained by(4)kij=KijμωBψ3,dij=DijμBψ3i,j=x,y,Kxx=Bd24ε2π2+16−π2ε21−ε2216ε2+π21−ε2,Dxx=Bd22ππ2+2π2ε2−16ε21−ε2216ε2+π21−ε2,Kxy=Bd2π−π2+2π2ε2+16−π2ε41−ε25/216ε2+π21−ε2,Dxy=Bd28επ2+12π2ε2−16ε21−ε2216ε2+π21−ε2,Kyx=Bd2ππ2+π2+32ε2+216−π2ε41−ε25/216ε2+π21−ε2,Dyx=Dxy,Kyy=Bd24επ2+π2+32ε2+216−π2ε41−ε2216ε2+π21−ε2,Dyy=Bd22ππ2+12π2ε2−16ε21−ε2216ε2+π21−ε2,where B and d are the width and diameter of the journal bearing, μ is the oil viscosity, ε is the eccentricity ratio, ψ is the clearance-to-radius ratio, and Kij and Dij (i, j = x, y) are the corresponding dimensionless stiffness and damping coefficient of journal bearing.


## 2.2. Model for Calculation of Time-Varying Meshing Stiffness

The planets are evenly distributed among the sun gear in the planetary gear train. There are varying mesh phasing factors in meshing pairs (i.e., sun-planet and ring-planet) because of the influence of the number of teeth. The mesh phasing factor affects the phase of the time-varying meshing stiffness [17]. Therefore, the mesh phasing factor is incorporated into the encased differential gear train. Following Parker and Lin [18], the mesh phasing factor is defined by(5)γsm=Zsψm2π,γrm=−Zrψm2π,where γsm and γrm are the mesh phasing factor of sun-planet and ring-planet, respectively, and Zs and Zr are the teeth number of the sun and ring, respectively.Maatar and Velex [16] proposed the time-varying meshing stiffness of helical gear pairs based on the time-varying total contact line length. Following Maatar and Velex [16], the time-varying meshing stiffness model of the helical gear is applied to the herringbone gear by a parallel method in this study. The time-varying meshing stiffness model of the herringbone gear is constructed as follows:(6)kt=2k0Lτ,where k0 is the average value of the time-varying meshing stiffness, τ is dimensionless time, τ = t/Tm, t and Tm are meshing time and period, and L(τ) is the instantaneous total contact length of gear pair, which is obtained by(7)Lτ=1+∑k=1SAkcos2πkτ−γ+Bksin2πkτ−γbεαcosβb,where γ is the mesh phasing factor, εα is the end-face contact ratio of the herringbone gear, βb is the helical angle of the base circle of the herringbone gear, B is the unilateral tooth width of the herringbone gear, Ak and Bk are the kth-order harmonic amplitude of the time-varying meshing stiffness, S represents the order of the Fourier series, and the series is truncated at S = 5 in this article.The parametersAk and Bk are calculated by(8)Ak=12εαεβπ2k2cos2πkεβ+cos2πkεα−cos2πkεα+εβ−1,Bk=12εαεβπ2k2sin2πkεβ+sin2πkεα−sin2πkεα+εβ,where εβ is axial contact ratio of the herringbone gear.


## 2.3. The Equivalent Displacement and Meshing Forces

The center displacement of each gear is projected into the meshing line direction. We ignored the change of the meshing angle caused by the movement of each gear center. Figure4 shows the positional relationship sun gear, planet, and composting ring gear in the differential stage.Figure 4
Position relationship of each component in the differential stage.The equivalent displacement of theith planet-composting ring and sun-planet along meshing line direction can be calculated by(9)Srpit=xpi−xr1+Hpisinαrp−Vpicosαrp,−Hr1sinαrp+ψi+Vr1cosαrp+ψi−erpit,Sspit=xs1−xpi+Hs1sinαsp−ψi+Vs1cosαsp−ψi,−Hpisinαsp−Vpicosαsp−espit,where erpi and espi are the equivalent error of the ith planet-composting ring and sun-planet along meshing line direction.According to the equivalent displacement calculation method of the differential stage, the equivalent displacement of thejth sun-planet and composting ring-planet along meshing line direction can be calculated by(10)Srmjt=xmj−xr2+Hmjsinαrm−Vmjcosαrm,−Hr2sinαrm+ψj+Vr2cosαrm+ψj−ermjt,Ssmjt=xs2−xmj+Hs2sinαsm−ψj+Vs2cosαsm−ψj,−Hmjsinαsm−Vmjcosαsm−esmjt.The planet is supported by planet bearing in the encased differential gear train, and thus the compression deformation between carrier and planet bearing can be obtained by(11)δHcpi=Hccosψi+Vcsinψi−Hpi,δVcpi=−Hcsinψi+Vccosψi−Vpi−xc,where δHcpi is the compression deformation of planet bearing along the lateral direction and δVcpi is the compression deformation along the longitudinal direction.Elastic and damping meshing forces in the differential and encased stages can be calculated by(12)Fit=KitSitDit=CiS˙it,i=spi,rpi,smj,andrmj,where S˙it is the first derivative of equivalent displacement. Ci is the meshing damping factor, which can be calculated by [19](13)Ci=2ξik0,i1/mpi+1/mgi,i=spi,rpi,smj,andrmj,where ζi denotes the constant damping ratios of the gear pair, ζi = 0.07; k0,i represents the average value of the meshing stiffness in the gear pair, and mpi and mgi are the masses of the driving and driven gears in the gear pair.


## 2.4. Differential Equations of the Compound Planetary Gear Trains

According to Figure2 of the dynamic model of the encased differential gear train, differential equations of 19 + 3 (N + M) degree-of-freedom are constructed for the encased differential gear train by the lumped parameter method. Differential equations of the differential stage can be deduced using Newton’s second law as follows.(1)
The sun gearZs1:(14)ms1x¨s1+∑i=1NFspit+Dspit=PD,Ms1H¨s1+∑i=1NFspit+Dspitsinαsp−ψi+Ks1Hs1=−Ms1gsinωHt,Ms1V¨s1+∑i=1NFspit+Dspitcosαsp−ψi+Ks1Vs1=−Ms1gcosωHt.(2)
Theith planet gear Zpi:(15)mpx¨pi−Fspit+Dspit+Frpit+Drpit=0,MpH¨pi−Fspit+Dspitsinαsp+Frpit+Drpitsinαrp+FHcpit=−MpigsinωHt+ψi+Fpi,MpV¨pi−Fspit+Dspitcosαsp+Frpit+Drpitcosαrp+FVcpit=−MpigcosωHt+ψi,whereFHcpi(t) is the oil film force of the ith planet bearing in the lateral direction, FVcpit is the oil film force of the ith planet bearing in the longitudinal direction, FHcpi(t) and FVcpit are calculated by equation (3), and Fpi is the centrifugal force of the ith planet Zpi.(3)
The composting ring gearZr1:(16)mr1x¨r1−∑i=1NFrpit+Drpit+kg1qxr1−xg1=0,Mr1H¨r1−∑i=1NFrpit+Drpitsinαrp+ψi+Kr1Hr1=−Mr1gsinωHt,Mr1V¨r1+∑i=1NFrpit+Drpitcosαrp+ψi+Kr1Vr1=−Mr1gcosωHt.(4)
The intermediate floating memberZg1 and floating ring Zf1:(17)mg1x¨g1−Kg1qxr1−xg1+Kf1qxg1−xf1=0,mf1x¨f1−Kf1qxg1−xf1+Kf1s2rf1bxf1rf1b−xs2rs2b=0.Differential equations of the encased stage can be deduced using Newton’s second law as follows:(1)
The sun gearZs2:(18)ms2x¨s2−Kf1s2rs2bxf1rf1b−xs2rs2b+∑j=1MFsmjt+Dsmjt=0,Ms2H¨s2+∑j=1MFsmjt+Dsmjtsinαsm−ψj+Ks2Hs2=0,Ms2V¨s2+∑j=1MFsmjt+Dsmjtcosαsm−ψj+Ks2Vs2=−Ms2g.(2)
Thejth planet gear Zmj:(19)mmx¨mj−Fsmjt+Dsmjt+Frmjt+Drmjt=0,MmH¨mj−Fsmjt+Dsmjtsinαsm+Frmjt+Drmjtsinαrm+FHmjt=−Mmgsinψj,MmV¨mj−Fsmjt+Dsmjtcosαsm+Frmjt+Drmjtcosαrm+FVmjt=−Mmgcosψj,whereFHmj(t) is the oil film force of the jth planet bearing in the lateral direction and FVmjt is the oil film force of the jth planet bearing in the longitudinal direction.(3)
The composting ring gearZr2:(20)mr2x¨r2−∑j=1MFrmjt+Drmjt+kg2qxr2−xg2=0,Mr2H¨r2−∑j=1MFrmjt+Drmjtsinαrm+ψj+Kr2Hr2=0,Mr2V¨r2+∑i=1MFrmjt+Drmjtcosαrm+ψj+Kr2Vr2=−Mr2g.(4)
The intermediate floating memberZg2 and floating ring Zf2:(21)mg2x¨g2−Kg2qxr2−xg2+Kf2qxg2−xf2=0,mf2x¨f2−Kf2qxg2−xf2+Kf2Lrf2bxf2rf2b−xcrc=0.The differential equation of carrierH can be deduced by Newton’s law:(22)mcx¨c+∑i=1NFVcpi+KcHrcxcrc−Kf2Hrcxf2rf2b−xcrc=−TLrc,McH¨c+∑i=1NFHcpicosψi−FVcpisinψi+FcH=0,McV¨c+∑i=1NFHcpisinψi+FVcpicosψi+FcV=−Mcg,where FcH is the oil film force of journal bearing in the lateral direction and FcVis the oil film force of journal bearing in the longitudinal direction.


## 3. Numerical Studies and Analysis

The dynamics simulation is performed with a constant power of 1000 kW. Different rotational speeds are applied to the differential stage sun gear. Table1 lists the parameters of the encased differential gear train. The system response is solved by the numerical calculation method of the Fourier series [20]. The eccentricity ratio of the journal bearing is taken from 0.01 to 0.91. The clearance-to-radius ratio ranges from 0.1‰ to 2.5‰. The oil viscosity is 0.04 N·s/m2.Table 1
The parameters of the encased differential gear train.
Z
s1/Zs2
Z
pi/Zmj
Z
r1/Zr2
Z
g
1
/
Z
g
2
Z
f1/Zf2
H
Teeth number
38/98
96/53
232/202
251/231
375/348
—
Module
2.5
2.5
2.5
2.5
2
—
Helix angle (°)
25
25
25
25
0
—
Pressure angel (°)
20
20
20
20
20
—
Torsional stiffness (N·μm/rad)
—
—
—
9.5e5/6.9e5
1.2e9/4.3e8
9238
Support stiffness (N/μm)
6.2/2241
—
158.8/115.4
—
—
—


### 3.1. Comparative Analysis of Spring Bearing and Journal Bearing

Before performing planet bearing parameter analysis, we calculate the dynamic deviation of the meshing force to investigate whether the encased differential gear train with journal bearing is more stable than the encased differential gear train with spring (i.e., planet bearings are replaced by fixed stiffness springs). The dynamic model of the encased differential gear train with spring is constructed by equations (14)–(22), and the stiffness of spring is 6 × 106 N/m. According to the stiffness and damping coefficient of the journal bearing in Section 2.1, we establish the dynamic model of the encased differential planetary gear train with journal bearing. The dynamic deviation of meshing force (i.e., the degree of dynamic meshing force deviates from the theoretical static meshing force) is defined as follows [21]:(23)Se=1N∑j=1N1/Q∑i=1QFjie−Fs2,Si=1N∑j=1N1/Q∑i=1QFjii−Fs2,where Se and Si represent the dynamic deviation of external and internal meshing forces, respectively, N is the number of planets, Q is the sampling times of the dynamic meshing force, Fije and Fijiare external and internal instantaneous meshing force, respectively, and Fs is the theoretical static meshing force.Figure5(a) shows that, in the differential stage, the dynamic deviation curves of the external meshing force of the bearing model and the spring model are basically the same at a low rotational speed. However, the dynamic deviation curve of the external meshing force of the spring model has a peak when the rotational speed is 4800 r/min. It is worth noting that the dynamic deviation curve of the internal meshing force is approximately the same as that of the external meshing force (see Figures 5(a) and 5(b)). Figures 5(c) and 5(d) detail that, in the encased stage, the dynamic deviation curves of the meshing forces of the bearing model and the spring model are substantially the same. These findings indicate that at high rotational speeds, the oil film force of journal bearing can curb the vibration of the encased differential gear train. Therefore, we investigate the influence of the parameters of planet bearings on dynamic meshing forces in the next step.Figure 5
The dynamic deviation of meshing force: (a) external dynamic deviation of the differential stage, (b) internal dynamic deviation of the differential stage, (c) external dynamic deviation of the encased stage, and (d) internal dynamic deviation of the encased stage.
(a)
(b)
(c)
(d)


### 3.2. Influence of the Differential Planet Bearing Parameters on Meshing Forces

After modeling differential equations of the encased differential gear train with journal bearing, the dynamic meshing force of differential equations can be calculated by equation (12). Besides, we use the root-mean-square to describe dynamic meshing forces.


#### 3.2.1. Clearance-to-Radius Ratio

To investigate the effects of the differential planet bearing clearance-to-radius ratios on meshing forces, we set the clearance-to-radius ratio as 0.1‰, 0.5‰, 1‰, 1.5‰, 2‰, and 2.5‰, respectively (see Figure6). The amplitude-frequency response curves of the meshing forces within the rotation speed range of 500–6000 r/min are calculated.Figure 6
Amplitude-frequency response curves of meshing forces under different clearance-to-radius ratios of the differential planet bearing: (a) external meshing force of the differential stage, (b) internal meshing force of the differential stage, (c) external meshing force of the encased stage, and (d) internal meshing force of the encased stage.
(a)
(b)
(c)
(d)Figures6(a) and 6(b) report that the amplitude-frequency response curves of meshing forces (i.e., external and internal) in the differential stage overlap at low rotational speed and diverge at a high rotational speed. This finding indicates that the impact of the clearance-to-radius ratio of the differential planet bearing on the meshing force is small when the rotational speed is low; the impact of clearance-to-radius ratio on the meshing force is obvious when the rotational speed is high. Figures 6(c) and 6(d) show that the amplitude-frequency response curves of meshing forces in the encased stage tend to diverge at low rotational speed and tend to overlap at a high rotational speed. Therefore, when the rotational speed in the encased planetary gear train is low, the clearance-to-radius ratio of the differential planet bearing should be reduced to diminish its influence on the meshing force to ensure the stability and reliability of the system.


#### 3.2.2. Eccentricity Ratio

To investigate the influence of the eccentricity ratio of the differential planet bearing on meshing forces, the eccentricity ratios amount on the planet bearing are taken as 0.01, 0.19, 0.37, 0.55, 0.73, and 0.91, respectively. We calculate the amplitude-frequency response curves of meshing forces within the rotation speed range of 500–6000 r/min, as shown in Figure7.Figure 7
Amplitude-frequency response curves of meshing forces under different eccentricity ratios of the differential planet bearing: (a) external meshing force of the differential stage, (b) internal meshing force of the differential stage, (c) external meshing force of the encased stage, and (d) internal meshing force of the encased stage.
(a)
(b)
(c)
(d)It can be seen from Figure7(a) that, for different eccentricity ratios of the differential planet bearing, amplitude-frequency response curves are basically consistent. This finding indicates that the eccentricity ratio of the differential planet bearing has little effect on the external meshing force in the differential planetary gear train. Figure 7(b) reveals that amplitude-frequency response curves tend to diverge when the rotational speed is between 2000 r/min and 5000 r/min. This result indicates that when the rotational speed is between 2000 r/min and 5000 r/min, the eccentricity ratio of the differential planet bearing has a significant influence on the internal meshing force in the differential stage. Amplitude-frequency response curves tend to diverge at low rotational speed and tend to overlap at high rotational speed (see Figures 7(c) and 7(d)), suggesting that the impact of eccentricity ratio of the differential planet bearing on the meshing forces in the encased stage is obvious when the rotational speed is low.


### 3.3. Influence of the Encased Planet Bearing Parameters on Meshing Forces

#### 3.3.1. Clearance-to-Radius Ratio

In this section, the impact of clearance-to-radius ratios of the encased planet bearing on meshing forces is explored. Similar to the differential planet bearing, we set the clearance-to-radius ratio of the encased planet bearing as 0.1‰, 0.5‰, 1‰, 1.5‰, 2‰, and 2.5‰, respectively, as shown in Figure8.Figure 8
Amplitude-frequency response curves of meshing forces under different clearance-to-radius ratios of the encased planet bearing: (a) external meshing force of the differential stage, (b) internal meshing force of the differential stage, (c) external meshing force of the encased stage, and (d) internal meshing force of the encased stage.
(a)
(b)
(c)
(d)Figures8(a) and 8(b) uncover that the amplitude-frequency response curves of meshing forces (i.e., external and internal) in the differential stage are substantially uniform for different clearance-to-radius ratios of the encased planet bearing. Figures 8(c) and 8(d) show that the amplitude-frequency response curves of meshing forces in the encased stage tend to overlap at low rotational speed and tend to diverge at a high rotational speed. It is worth noting that the amplitude-frequency response curves of meshing forces in the encased stage have a peak when the rotational speed is between 1500 r/min and 2000 r/min. These findings indicate that when the rotational speed is between 1500 r/min and 2000 r/min, the clearance-to-radius ratio of the encased planet bearing has an obvious impact on the meshing forces in the encased stage, especially 1‰ and 0.5‰ clearance-to-radius ratios.


#### 3.3.2. Eccentricity Ratio

Using the same eccentricity parameters as in Section3.2.2, the impact of the eccentricity ratio of the encased planet bearing on dynamic meshing forces is investigated, as shown in Figure 9.Figure 9
Amplitude-frequency response curves of meshing forces under different eccentricity ratios of the encased planet bearing: (a) external meshing force of the differential stage, (b) internal meshing force of the differential stage, (c) external meshing force of the encased stage, and (d) internal meshing force of the encased stage.
(a)
(b)
(c)
(d)Figures9(a) and 9(b) show that the corresponding amplitude-frequency response curves are the same regardless of which eccentricity ratio of the encased planet bearing is incorporated into our model. Similar to the clearance-to-radius ratio of the differential planet bearing, the eccentricity ratio of the encased planet bearing has a significant impact on the meshing forces in the encased stage when the rotational speed is between 1500 r/min and 2000 r/min (see Figures 9(c) and 9(d)). Besides, the eccentricity ratio of the encased planet bearing has more influence on the internal meshing force than on the external meshing force. Compared with Section 3.2.2, it can be found that the influence of the eccentricity ratio of the encased planet bearing on the meshing force of the encased stage is more significant than that of the differential stage.


## 3.1. Comparative Analysis of Spring Bearing and Journal Bearing

Before performing planet bearing parameter analysis, we calculate the dynamic deviation of the meshing force to investigate whether the encased differential gear train with journal bearing is more stable than the encased differential gear train with spring (i.e., planet bearings are replaced by fixed stiffness springs). The dynamic model of the encased differential gear train with spring is constructed by equations (14)–(22), and the stiffness of spring is 6 × 106 N/m. According to the stiffness and damping coefficient of the journal bearing in Section 2.1, we establish the dynamic model of the encased differential planetary gear train with journal bearing. The dynamic deviation of meshing force (i.e., the degree of dynamic meshing force deviates from the theoretical static meshing force) is defined as follows [21]:(23)Se=1N∑j=1N1/Q∑i=1QFjie−Fs2,Si=1N∑j=1N1/Q∑i=1QFjii−Fs2,where Se and Si represent the dynamic deviation of external and internal meshing forces, respectively, N is the number of planets, Q is the sampling times of the dynamic meshing force, Fije and Fijiare external and internal instantaneous meshing force, respectively, and Fs is the theoretical static meshing force.Figure5(a) shows that, in the differential stage, the dynamic deviation curves of the external meshing force of the bearing model and the spring model are basically the same at a low rotational speed. However, the dynamic deviation curve of the external meshing force of the spring model has a peak when the rotational speed is 4800 r/min. It is worth noting that the dynamic deviation curve of the internal meshing force is approximately the same as that of the external meshing force (see Figures 5(a) and 5(b)). Figures 5(c) and 5(d) detail that, in the encased stage, the dynamic deviation curves of the meshing forces of the bearing model and the spring model are substantially the same. These findings indicate that at high rotational speeds, the oil film force of journal bearing can curb the vibration of the encased differential gear train. Therefore, we investigate the influence of the parameters of planet bearings on dynamic meshing forces in the next step.Figure 5
The dynamic deviation of meshing force: (a) external dynamic deviation of the differential stage, (b) internal dynamic deviation of the differential stage, (c) external dynamic deviation of the encased stage, and (d) internal dynamic deviation of the encased stage.
(a)
(b)
(c)
(d)


## 3.2. Influence of the Differential Planet Bearing Parameters on Meshing Forces

After modeling differential equations of the encased differential gear train with journal bearing, the dynamic meshing force of differential equations can be calculated by equation (12). Besides, we use the root-mean-square to describe dynamic meshing forces.


### 3.2.1. Clearance-to-Radius Ratio

To investigate the effects of the differential planet bearing clearance-to-radius ratios on meshing forces, we set the clearance-to-radius ratio as 0.1‰, 0.5‰, 1‰, 1.5‰, 2‰, and 2.5‰, respectively (see Figure6). The amplitude-frequency response curves of the meshing forces within the rotation speed range of 500–6000 r/min are calculated.Figure 6
Amplitude-frequency response curves of meshing forces under different clearance-to-radius ratios of the differential planet bearing: (a) external meshing force of the differential stage, (b) internal meshing force of the differential stage, (c) external meshing force of the encased stage, and (d) internal meshing force of the encased stage.
(a)
(b)
(c)
(d)Figures6(a) and 6(b) report that the amplitude-frequency response curves of meshing forces (i.e., external and internal) in the differential stage overlap at low rotational speed and diverge at a high rotational speed. This finding indicates that the impact of the clearance-to-radius ratio of the differential planet bearing on the meshing force is small when the rotational speed is low; the impact of clearance-to-radius ratio on the meshing force is obvious when the rotational speed is high. Figures 6(c) and 6(d) show that the amplitude-frequency response curves of meshing forces in the encased stage tend to diverge at low rotational speed and tend to overlap at a high rotational speed. Therefore, when the rotational speed in the encased planetary gear train is low, the clearance-to-radius ratio of the differential planet bearing should be reduced to diminish its influence on the meshing force to ensure the stability and reliability of the system.


### 3.2.2. Eccentricity Ratio

To investigate the influence of the eccentricity ratio of the differential planet bearing on meshing forces, the eccentricity ratios amount on the planet bearing are taken as 0.01, 0.19, 0.37, 0.55, 0.73, and 0.91, respectively. We calculate the amplitude-frequency response curves of meshing forces within the rotation speed range of 500–6000 r/min, as shown in Figure7.Figure 7
Amplitude-frequency response curves of meshing forces under different eccentricity ratios of the differential planet bearing: (a) external meshing force of the differential stage, (b) internal meshing force of the differential stage, (c) external meshing force of the encased stage, and (d) internal meshing force of the encased stage.
(a)
(b)
(c)
(d)It can be seen from Figure7(a) that, for different eccentricity ratios of the differential planet bearing, amplitude-frequency response curves are basically consistent. This finding indicates that the eccentricity ratio of the differential planet bearing has little effect on the external meshing force in the differential planetary gear train. Figure 7(b) reveals that amplitude-frequency response curves tend to diverge when the rotational speed is between 2000 r/min and 5000 r/min. This result indicates that when the rotational speed is between 2000 r/min and 5000 r/min, the eccentricity ratio of the differential planet bearing has a significant influence on the internal meshing force in the differential stage. Amplitude-frequency response curves tend to diverge at low rotational speed and tend to overlap at high rotational speed (see Figures 7(c) and 7(d)), suggesting that the impact of eccentricity ratio of the differential planet bearing on the meshing forces in the encased stage is obvious when the rotational speed is low.


## 3.2.1. Clearance-to-Radius Ratio

To investigate the effects of the differential planet bearing clearance-to-radius ratios on meshing forces, we set the clearance-to-radius ratio as 0.1‰, 0.5‰, 1‰, 1.5‰, 2‰, and 2.5‰, respectively (see Figure6). The amplitude-frequency response curves of the meshing forces within the rotation speed range of 500–6000 r/min are calculated.Figure 6
Amplitude-frequency response curves of meshing forces under different clearance-to-radius ratios of the differential planet bearing: (a) external meshing force of the differential stage, (b) internal meshing force of the differential stage, (c) external meshing force of the encased stage, and (d) internal meshing force of the encased stage.
(a)
(b)
(c)
(d)Figures6(a) and 6(b) report that the amplitude-frequency response curves of meshing forces (i.e., external and internal) in the differential stage overlap at low rotational speed and diverge at a high rotational speed. This finding indicates that the impact of the clearance-to-radius ratio of the differential planet bearing on the meshing force is small when the rotational speed is low; the impact of clearance-to-radius ratio on the meshing force is obvious when the rotational speed is high. Figures 6(c) and 6(d) show that the amplitude-frequency response curves of meshing forces in the encased stage tend to diverge at low rotational speed and tend to overlap at a high rotational speed. Therefore, when the rotational speed in the encased planetary gear train is low, the clearance-to-radius ratio of the differential planet bearing should be reduced to diminish its influence on the meshing force to ensure the stability and reliability of the system.


## 3.2.2. Eccentricity Ratio

To investigate the influence of the eccentricity ratio of the differential planet bearing on meshing forces, the eccentricity ratios amount on the planet bearing are taken as 0.01, 0.19, 0.37, 0.55, 0.73, and 0.91, respectively. We calculate the amplitude-frequency response curves of meshing forces within the rotation speed range of 500–6000 r/min, as shown in Figure7.Figure 7
Amplitude-frequency response curves of meshing forces under different eccentricity ratios of the differential planet bearing: (a) external meshing force of the differential stage, (b) internal meshing force of the differential stage, (c) external meshing force of the encased stage, and (d) internal meshing force of the encased stage.
(a)
(b)
(c)
(d)It can be seen from Figure7(a) that, for different eccentricity ratios of the differential planet bearing, amplitude-frequency response curves are basically consistent. This finding indicates that the eccentricity ratio of the differential planet bearing has little effect on the external meshing force in the differential planetary gear train. Figure 7(b) reveals that amplitude-frequency response curves tend to diverge when the rotational speed is between 2000 r/min and 5000 r/min. This result indicates that when the rotational speed is between 2000 r/min and 5000 r/min, the eccentricity ratio of the differential planet bearing has a significant influence on the internal meshing force in the differential stage. Amplitude-frequency response curves tend to diverge at low rotational speed and tend to overlap at high rotational speed (see Figures 7(c) and 7(d)), suggesting that the impact of eccentricity ratio of the differential planet bearing on the meshing forces in the encased stage is obvious when the rotational speed is low.


## 3.3. Influence of the Encased Planet Bearing Parameters on Meshing Forces

### 3.3.1. Clearance-to-Radius Ratio

In this section, the impact of clearance-to-radius ratios of the encased planet bearing on meshing forces is explored. Similar to the differential planet bearing, we set the clearance-to-radius ratio of the encased planet bearing as 0.1‰, 0.5‰, 1‰, 1.5‰, 2‰, and 2.5‰, respectively, as shown in Figure8.Figure 8
Amplitude-frequency response curves of meshing forces under different clearance-to-radius ratios of the encased planet bearing: (a) external meshing force of the differential stage, (b) internal meshing force of the differential stage, (c) external meshing force of the encased stage, and (d) internal meshing force of the encased stage.
(a)
(b)
(c)
(d)Figures8(a) and 8(b) uncover that the amplitude-frequency response curves of meshing forces (i.e., external and internal) in the differential stage are substantially uniform for different clearance-to-radius ratios of the encased planet bearing. Figures 8(c) and 8(d) show that the amplitude-frequency response curves of meshing forces in the encased stage tend to overlap at low rotational speed and tend to diverge at a high rotational speed. It is worth noting that the amplitude-frequency response curves of meshing forces in the encased stage have a peak when the rotational speed is between 1500 r/min and 2000 r/min. These findings indicate that when the rotational speed is between 1500 r/min and 2000 r/min, the clearance-to-radius ratio of the encased planet bearing has an obvious impact on the meshing forces in the encased stage, especially 1‰ and 0.5‰ clearance-to-radius ratios.


### 3.3.2. Eccentricity Ratio

Using the same eccentricity parameters as in Section3.2.2, the impact of the eccentricity ratio of the encased planet bearing on dynamic meshing forces is investigated, as shown in Figure 9.Figure 9
Amplitude-frequency response curves of meshing forces under different eccentricity ratios of the encased planet bearing: (a) external meshing force of the differential stage, (b) internal meshing force of the differential stage, (c) external meshing force of the encased stage, and (d) internal meshing force of the encased stage.
(a)
(b)
(c)
(d)Figures9(a) and 9(b) show that the corresponding amplitude-frequency response curves are the same regardless of which eccentricity ratio of the encased planet bearing is incorporated into our model. Similar to the clearance-to-radius ratio of the differential planet bearing, the eccentricity ratio of the encased planet bearing has a significant impact on the meshing forces in the encased stage when the rotational speed is between 1500 r/min and 2000 r/min (see Figures 9(c) and 9(d)). Besides, the eccentricity ratio of the encased planet bearing has more influence on the internal meshing force than on the external meshing force. Compared with Section 3.2.2, it can be found that the influence of the eccentricity ratio of the encased planet bearing on the meshing force of the encased stage is more significant than that of the differential stage.


## 3.3.1. Clearance-to-Radius Ratio

In this section, the impact of clearance-to-radius ratios of the encased planet bearing on meshing forces is explored. Similar to the differential planet bearing, we set the clearance-to-radius ratio of the encased planet bearing as 0.1‰, 0.5‰, 1‰, 1.5‰, 2‰, and 2.5‰, respectively, as shown in Figure8.Figure 8
Amplitude-frequency response curves of meshing forces under different clearance-to-radius ratios of the encased planet bearing: (a) external meshing force of the differential stage, (b) internal meshing force of the differential stage, (c) external meshing force of the encased stage, and (d) internal meshing force of the encased stage.
(a)
(b)
(c)
(d)Figures8(a) and 8(b) uncover that the amplitude-frequency response curves of meshing forces (i.e., external and internal) in the differential stage are substantially uniform for different clearance-to-radius ratios of the encased planet bearing. Figures 8(c) and 8(d) show that the amplitude-frequency response curves of meshing forces in the encased stage tend to overlap at low rotational speed and tend to diverge at a high rotational speed. It is worth noting that the amplitude-frequency response curves of meshing forces in the encased stage have a peak when the rotational speed is between 1500 r/min and 2000 r/min. These findings indicate that when the rotational speed is between 1500 r/min and 2000 r/min, the clearance-to-radius ratio of the encased planet bearing has an obvious impact on the meshing forces in the encased stage, especially 1‰ and 0.5‰ clearance-to-radius ratios.


## 3.3.2. Eccentricity Ratio

Using the same eccentricity parameters as in Section3.2.2, the impact of the eccentricity ratio of the encased planet bearing on dynamic meshing forces is investigated, as shown in Figure 9.Figure 9
Amplitude-frequency response curves of meshing forces under different eccentricity ratios of the encased planet bearing: (a) external meshing force of the differential stage, (b) internal meshing force of the differential stage, (c) external meshing force of the encased stage, and (d) internal meshing force of the encased stage.
(a)
(b)
(c)
(d)Figures9(a) and 9(b) show that the corresponding amplitude-frequency response curves are the same regardless of which eccentricity ratio of the encased planet bearing is incorporated into our model. Similar to the clearance-to-radius ratio of the differential planet bearing, the eccentricity ratio of the encased planet bearing has a significant impact on the meshing forces in the encased stage when the rotational speed is between 1500 r/min and 2000 r/min (see Figures 9(c) and 9(d)). Besides, the eccentricity ratio of the encased planet bearing has more influence on the internal meshing force than on the external meshing force. Compared with Section 3.2.2, it can be found that the influence of the eccentricity ratio of the encased planet bearing on the meshing force of the encased stage is more significant than that of the differential stage.


## 4. Conclusions

Considering the factors of time-varying meshing stiffness and comprehensive transmission error, the dynamic model of encased differential gear train with journal bearing is established. The asymmetry and interaction of the oil film stiffness and damping of the planet bearings are employed in this dynamic model. The effects of planet bearing parameters on the meshing forces are explored, and the main conclusions are drawn as follows.When the rotational speed is low, the stability of the encased differential gear train with journal bearing is basically the same as that of the encased differential gear train with spring. However, when the rotational speed is high, the encased differential gear train with journal bearing is more stable than the encased differential gear train with spring.The clearance-to-radius and eccentricity ratio of the differential planet bearings affect the meshing forces of both the differential stage and encased stage. When the rotational speed is low, the influence of the clearance ratio and eccentricity ratio of the differential planet bearings on the meshing forces of the differential stage is less than that of the encased stage. However, this phenomenon is the opposite when the rotational speed is high.The clearance-to-radius and eccentricity ratio of the encased planet bearing obviously affect the forces meshing of the encased stage, while they slightly affect the meshing forces of the differential stage. When the rotational speed is between 1400 r/min and 2000 r/min, the clearance-to-radius and eccentricity ratio of the encased planet bearing have a significant impact on the meshing forces of the encased stage.



---

*Source: 2436191-2020-04-20.xml*