# Bifurcation and Chaos Analysis of the Spur Gear Transmission System for One-Way Clutch, Two-Shaft Assembly

Authors: Zhihui Liu; Hongzhi Yan; Yuming Cao; Yuqing Lai
Journal: Shock and Vibration
 (2017)

---

## Abstract

A four-degree-of-freedom nonlinear transverse and torsional vibration model of spur gear transmission system for one-way clutch, two-shaft assembly was developed, in which the one-way clutch was modeled as a piecewise nonlinear spring with discontinuous stiffness, considering the factors such as the time-varying gear mesh stiffness, static transmission error, and nonlinearity backlash. With the help of bifurcation diagrams, time domain response diagrams, phase plane diagrams, and Poincaré maps, the effects of the excitation frequency and the torsional stiffness of one-way clutch on the dynamic behavior of gear transmission system for one-way clutch, two-shaft assembly are investigated in detail by using Runge-Kutta method. Numerical results reveal that the system response involves period-1 motion, multiperiodic motion, bifurcation, and chaotic motion. Large torsional stiffness of one-way clutch can increase the impact and lead to instability in the system. The results can present a useful source of reference for technicians and engineers for dynamic design and vibration control of such system.


---

## Body

## 1. Introduction

The one-way clutch, two-shaft assembly is the key component of the two-speed transmission, which is widely used in loaders, helicopters, tanks, military trucks, and other gear transmission systems. The two-speed transmission is generally composed of the outer ring gear, intermediate input shaft, inner star wheel, and a plurality of cylinder rollers or sprags. The one-way clutch is made up of the outer ring gear, sprags, and the inner star wheel. The outer ring gear and the inner star wheel can be engaged or disengaged according to their relative angular speed, which can lead to the two-speed transmission having two kinds of output speed.It is well known that gear mechanisms are widely used owing to their rather large technical advantages. However, they are known as one of the most important sources of vibration and noise in mechanical transmission systems. The vibration, noise, and wear performances are mainly affected by the gear pair dynamic characteristics.There is a growing literature on nonlinear dynamics of gear pairs in the past decades. And to analyze the nonlinear dynamic behaviors of gear transmission system, the different excitation of static transmission error, dynamic meshing stiffness, bearing contact stiffness, prime mover, clearance nonlinearity, and load fluctuation are added [1–7]. Walha et al. [8] investigated the dynamics of a two-stage gear system with mesh stiffness and backlash. Parker et al. [9] used a finite element/contact mechanic model to investigate the dynamic response of a spur gear pair. And the dynamic characteristics of a gear pair are analyzed across a wide range of frequency and torques. Al-Shyyab and Kahraman [10] adopted lumped parameter method to establish a nonlinear time-varying dynamic model for multimesh gear, considering the nonlinear gear gap and gear meshing stiffness. The multiterm harmonic balance method (HBM) was used to gain the steady-state period-1 response, including the Root Mean Square (RMS) and the mean value of the vibration amplitude across a wide range of frequency, taking into account different gear meshing stiffness. He et al. [4] analyzed the dynamic characteristics of a spur gear pair, considering the dynamic meshing stiffness and sliding friction. The influence of friction on dynamic transmission error is discussed. In recent years, many scholars have studied the bifurcation and chaos characteristics of gear transmission system [11]. Wang et al. [12] made use of the bifurcation diagrams, phase plane diagrams, Poincaré maps, time domain response diagrams, and amplitude-frequency spectrums to investigate the effects of the pinion speed and stiffness on the dynamic behavior of gear transmission system for locomotive. Xiang et al. [13] built a six-degree-of-freedom nonlinear dynamic model of a spur gear-bearing system with time-varying stiffness, gear backlash, and surface friction. In [13], bifurcation diagrams of the mesh stiffness and the rotational speed ratio as control parameters were used to observe the dynamic orbits of the system. The chaotic motion was identified from the largest Lyapunov exponents, Poincaré maps, and the phase diagrams of the system. Chang-Jian and Chang [14] studied the dynamic orbits of a single-degree-of-freedom spur gear system by using bifurcation diagrams with the dimensionless damping coefficient and the dimensionless rotational speed ratio as control parameters. Gou et al. [15] used the simple cell mapping method and escape time algorithm to examine the bifurcation dynamic characteristics of the multiparameter coupling in torsion-vibration gear system. Although most scholars study the steady-state response analysis of gear transmission system, Khabou et al. [16] studied a spur gear dynamic behavior in the transient regime.The dynamic analysis of one-way clutch, two-shaft assembly is to study the coupling nonlinear characteristics of one-way clutch and gear transmission system. Though application of one-way clutch, two-shaft assembly is extensive, literature on its dynamics is limited. In the aspects of dynamic characteristics study of one-way clutch, Vernay et al. [17] pointed out and quantified the possible relative sliding effects between sprag-type races of one-way clutches by experiments. Gill-Jeong [18] verified the nonlinear behavior analysis of a paired spur gear system with a one-way clutch for the first time. He found that the dynamic transmission error of a gear system with a one-way clutch was reduced over the entire frequency range compared to a pure gear system with three-degree-of-freedom torsional vibration model. However, the effects of parameters such as bearing stiffness, bearing backlash, torsion stiffness, and damping of clutch on the nonlinear behavior of gear system coupled with a one-way clutch were not considered. In order to investigate the effect of the one-way clutch on the dynamic performance of the coupling transmission, a great deal of research effort has been devoted to the study of the nonlinear dynamics of one-way clutches in belt-pulley systems. In such system, the one-way clutch is usually modeled as a two-degree-of-freedom model, which is regarded as piecewise linear [19]. Zhu and Parker [20] used a two-degree-of-freedom model to examine the nonlinear dynamics of a two-pulley belt system with a one-way clutch through three methods of multiterm harmonic balance, numerical integration, and the bifurcation software AUTO separately. Zhu and Parker [21] also used the method of multiple scales to obtain the periodic solutions of a two-pulley belt system with clearance-type nonlinearity. The stable steady-state periodic responses of a belt-drive system with a one-way clutch were studied comprehensively [22–25]. Ding and Zu [22] studied the effects of the system parameters on the nonlinear steady-state response of a belt-drive system with a one-way clutch and had a conclusion that one-way clutch not only decreases the resonance amplitude of the driven pulley and shaft’s rotational vibration but also reduces the resonance region of the belt’s transverse vibration. Ding [23] also presented the stable steady-state periodic response of a two-pulley belt-drive system coupled with one-way clutch under double excitation for the first time. The Galerkin method was adapted to exam whether the one-way device can efficiently reduce the rotation vibration and transverse vibration of the pulley belt dynamic system either.In view of the lack of research for dynamic analysis for one-way clutch, two-shaft assembly, the present paper concentrates on building and analyzing a mathematical model for such system. For the first time, a four-degree-of-freedom nonlinear transverse and torsional vibration model of spur gear transmission system for one-way clutch, two-shaft assembly was modeled, considering the gear massesm1 and m2, a nonlinear mesh stiffness kh(t), and mesh damping ch. The one-way clutch is modeled as a piecewise nonlinear spring with discontinuous stiffness that separates the driven gear into two degrees of freedom. The main purpose is to inspect the influences of the gear speed and the torsional stiffness of one-way clutch on the system dynamic behavior and find out some useful reference information for technicians and engineers to dynamic design and vibration control of such system.


## 2. Working Principle of One-Way Clutch, Two-Shaft Assembly

Figure1 shows the structural schematic of one-way clutch, two-shaft assembly, which is composed primarily of bearing (1, 2), gear (1, 2), sprags, input shaft, output shaft, and inner ring. Gear 1 is fixed on the input shaft by the method of spline interference fit. Gear 2, sprags, and the inner ring, which is connected to the output shaft through key connection, are formed into a one-way clutch. Figure 2 is used as a reference to explain the work principle of the one-way clutch. When the one-way clutch is not working, the tension of the spring ensures that the sprags keep in touch with gear 2 and the inner ring. When gear 2 begins to rotate clockwise under the action of gear 1 or the rotation speed of gear 2 exceeds the rotation speed of the inner ring, the sprag will rotate around its center in clockwise direction with the spring tension and the tangential friction at the contact point. With the sprag rotating, the distance of the two contact points between the sprag and gear 2 and the inner ring would become longer, and the power can be transmitted from gear 2 to the inner ring. As gear 2 and the inner ring are fully wedged by the sprags, the rotation speed of gear 2 and the inner ring will be the same. Conversely, gear 2 will be overrunning and the power could not be transmitted to the inner ring when gear 2 rotates anticlockwise or the rotation speed of gear 2 is slower than the rotation speed of the inner ring. So it can be concluded that the rotation direction of gear 1 determines whether the load can be transferred for the one-way clutch, two-shaft assembly.Figure 1
Structural diagram of one-way clutch, two-shaft assembly.Figure 2
Structural schematic of one-way clutch.


## 3. Mathematical Model

Figure3 shows a schematic illustration of the dynamic model of one-way clutch, two-shaft assembly considered in the present analysis. In this model, r1 is the base circle radius of gear 1, r2 is the base circle radius of gear 2, I1 is the mass moments of inertia of gear 1, I2 is the mass moments of inertia of gear 2, Io is the mass moment of inertia of inner ring, m1 is the mass of gear 1, m2 is the mass of gear 2, ch is the damping coefficient of the gear mesh, kh(t) is the time-varying stiffness coefficient of the gear mesh, which can be obtained by assuming a rectangular wave as shown in (2) [8, 18], where kh is the average mesh stiffness value, khr and ϕr are the rth Fourier coefficient and phase angle of kh(t), respectively, and kh/ktp=εr,khr/ktp=2-2cos⁡(2πr(εr-1))/(πr),ϕr=atan⁡((1-cos(2πr(εr-1)))/(sin(2πr(εr-1)))), εr is the involute contact ratio, R=5, fh is the force-displacement on the mesh line function of gear mesh, c1 is the damping coefficient of bearing 1, c2 is the damping coefficient of bearing 2, ki, fi are the contact stiffness coefficient and the force-transverse displacement function of ith bearing (i=1,2), Kc is the torsional stiffness coefficient of one-way clutch, and Cc is the torsional damping coefficient of one-way clutch. From the lumped parameter mass method and Newton’s laws of motion, the five-degree-of-freedom nonlinear model of the geared clutch bearing system dynamic model can show that (1)m1y1′′+c1y1′+chr1θ1-r2θ2′+y1′-y2′-e′+k1f1y1+khtfhr1θ1-r2θ2+y1-y2-e=-F1,m2y2′′+c2y2′-chr1θ1-r2θ2′+y1′-y2′-e′+k2f2y2-khtfhr1θ1-r2θ2+y1-y2-e=F2,I1θ1′′+chr1θ1-r2θ2′+y1′-y2′-e′·r1+khtfhr1θ1-r2θ2+y1-y2-e·r1=T1t,I2θ2′′-chr1θ1-r2θ2′+y1′-y2′-e′·r2-khtfhr1θ1-r2θ2+y1-y2-e·r2=-T2t,Ioθo′′=T2t-To,where gear 1 and gear 2 are assumed perfect involute spur gears with no modifications and the torsional and transverse flexibility of the input shaft and output shaft are neglected. Here, yi and θi are the transverse and torsional displacements of the ith gear (i=1,2), θo is the torsional displacement of the inner ring, and ′ means derivative with respect to time t. The input torque T1(t) is assumed to be fluctuating as T1(t)=T1m+T1aT(t), where T1m is the average torque and T1aT(t) is the fluctuating torque, considered as T1aT(t)=T1aTsin⁡(ωat+ϕT). Forcing term T2(t), which is connected with the piecewise nonlinear force-displacement function fc(x) shown in (4), and force-velocity function fc′(x) [26], is the torque transferred from gear 2 to the inner ring, as shown in (3). And the output torque To(t) applied on the inner ring or output shaft is assumed to be constant: that is, To(t)=Tom. The external radial forces F1 and F2 are preloads which are applied to both rolling element bearings. The static transmission error e(t), connected with the manufacturing and assembling of the gears, is considered sinusoidal excitation as et=e~sin⁡ωet+ϕh(2)kht=kht+2πω=kh+∑r=1Rkhrcos⁡rωt-ϕr,(3)T2t=Kcfcθ2-θo3/2+Ccfc′θ2-θofcθ2-θo1/4,(4)fcx=xx>0,0x≤0.Figure 3
Schematic diagram of one-way clutch, two-shaft assemblies.It is assumed thatr1θ1=yg1, r2θ2=yg2, and r2θo=yco, where the parameters of yg1, yg2, and yco are defined as the displacement on the direction of mesh line of gear 1, gear 2 and the inner ring, respectively. Through defining new variables q(t) and qo(t) as q(t)=yg1-yg2+y1-y2-e(t) and qo(t)=yg2-yco, (1) can be simplified further as (5)m1y1′′+c1y1′+chq′+k1f1y1+khtfhq=-F1,m2y2′′+c2y2′-chq′+k2f2y2-khtfhq=F2,mcq′′-mcy1′′+mcy2′′+chq′+khtfhq=mcF1tr12I1+F2tr22I2-mce′′t,qo′′-chr22I2q′-khtr22I2fhq=-F2tr221I2+1Io+For22Io.Here,mc is the equivalent gear pair mass. The piecewise linear displacement functions fi(yi) and fh(q) include the half of backlash of bearings bi(i=1,2) and the half of backlash of the gears bh, which can be expressed as (7) and (8).(6)mc=1r12/I1+r22/I2,(7)fiyi=yi-biyi>bi0yi≤biyi+biyi<-bii=1,2,(8)fhq=q-bhq>bh0q≤bhq+bhq<-bh.A dimensionless form of (5) can be obtained by letting (9)zit=yitbc,γ1t=qtbc,γ2t=qotbc,ωi=kimii=1,2,ωn=khmc,τ=ωnt,ωh=ωeωn,ωT=ωaωn,ξii=ci2miωn,ξi3=ch2miωni=1,2,ξ33=ch2mcωn,ξ43=chr222I2ωn,kii=ωi2ωn2,ki3=mcmi+∑r=15khrmiωn2cos⁡rωhτ-ϕri=1,2,k33=1+∑r=15khrmcωn2cos⁡rωhτ-ϕr,k43=khr22I2ωn2+∑r=15khrr22I2ωn2cos⁡rωhτ-ϕr,where bc is the characteristic length and the terms ωa and ωe are the fundamental excitation frequencies of external torque fluctuations and the internal static transmission error, respectively. And the dimensionless governing equations of motion are shown as (10), where zi is the dimensionless transverse displacement of the ith gear (i=1,2), γ1 is dimensionless difference between the dynamic transmission error and the static transmission error of the two gears, and γ2 is the dimensionless dynamic transmission error between gear 2 and the inner ring.


## 4. Parametric Studies on Characteristics of Bifurcation and Chaos

In this section, a case study is analyzed to demonstrate the effect of the rotation speed of gear 2 and the torsional stiffness of one-way clutch on its dynamic characteristics.For the structure of one-way clutch, two-shaft assemblies, gear 2 and the inner ring can be coupled or decoupled as the one-way clutch is active.The main parameters of the one-way clutch, two-shaft assemblies are listed in Table1. Equation (10) describes a strongly nonlinear dynamic system with time-varying gear mesh stiffness, nonlinear backlash, and nonlinear force-displacement function of one-way clutch. And it is solved by using the fourth-order Runge-Kutta method. In order to obtain a better understanding of the dynamic behavior of the one-way clutch, two-shaft assembly, the time series data were used to generate the bifurcation diagram, phase diagram, and Poincaré maps. The periodic responses are shown by varying the excitation frequency and the torsional stiffness of one-way clutch. On the other hand, since the new introduced coordinatesz1 andz2 have the similar nonlinear behavior, the dynamic response of the transverse vibrationz2 of gear 2 is selected to analyze the bifurcation and chaos characteristics. In order to find out the dynamic response of the windup angle of one-way clutch, the character of dynamic transmission error between gear 2 and the inner ring is chosen to be analyzed. Furthermore, in order to avoid the risk of obtaining transient solutions and minimize computing, the initial conditions of z2, z2′, γ1, γ1′, γ2, and γ2′ are all set to be zero.(10)10000100-11100001z1′′τz2′′τγ1′′τγ2′′τ+2ξ110ξ1300ξ22-ξ23000ξ33000-ξ430z1′τz2′τγ1′τγ2′τ+k110k1300k22-k23000k33000-k430f1z1τf2z2τfhγ1τ0=Fτ,Fτ=-F^1F^2mcr12I1F^1m+mcr22I2F^21-F^22r221I2+1Io+F^or22Io+00Fahωh20sin⁡ωhτ+ϕh+00mcr12I1F^1aT0sin⁡ωTτ+ϕT,F^i=Fimibcωn2,i=1,2,F^1m=F1mmcbcωn2,Fah=e~bc,F^1aT=F1aTmcbcωn2,F^o=Fobcωn2,F^21=Kcfcγ23/2bc3/2r25/2mcbcωn2+Ccfc′γ2fcγ21/4bc5/4r29/4mcbcωn,F^22=Kcfcγ23/2bc3/2r25/2bcωn2+Ccfc′γ2fcγ21/4bc5/4r29/4bcωn.Table 1
Properties of one-way clutch, two-shaft assembly.
Item
Notation
Value
Gears
Involute contact ratio
εr
1.5
Radius of gears
r
1, r2
0.04 m
Rotational inertia of gears
I
1, I2
0.0008 kg⋅m2
Mesh stiffness
k
h
2 × 108 N/m
Gear mesh damping
c
h
707
Bearing contact damping
c
1, c2
566
Quality of gears
m
1, m2
1 kg
One-way clutch
Rotational inertia of clutch
I
o
0.0008 kg⋅m2
Torsional stiffness
K
c
=
K
c
o
2467 N⋅m/rad
Torsional damping
C
c
0.217


### 4.1. Effect of Excitation Frequency

The ratio of the bearing stiffness to average gear mesh stiffnessk^ is introduced to represent the change of support stiffness. The frequency parameter ωh is assigned as control parameters. The bifurcation diagram of the dimensionless transverse displacement of gear 2 with excitation frequency ωh change is presented in Figure 4.Figure 4
z
2 versus ωh, plots for Fah=0.5, FaT=0, bh=0, b1=b2=bc=0.01, k^=5, Kc=5Kco, and Fm=1.It is well known that the jump phenomenon usually exists in the dynamic response of nonlinear system. It could be seen clearly in Figure4 that there is a phenomenon of jumping at value of ωh=0.61. The one-way clutch, two-shaft assemblies system mainly exhibits period-one motion at low values or at high excitation frequency: that is, ωh<0.69 or ωh>1.66. Just as shown in Figure 5(a), the time domain response diagram is a sine wave, the phase plane diagram is only one closed circle, and the Poincaré map is a single point. It is observed that the system goes into chaotic region within the range of 0.7≤ωh≤0.87,0.94≤ωh≤1.13, and 1.3≤ωh≤1.59. It can be seen in Figure 5(b) that the time domain response diagram shows a nonperiodic motion, the phase plane diagram is disorder, and the Poincaré map shows many discrete points. As ωh is increased from 0.88 to 0.93 and from 1.14 to 1.19, the system exhibits unstable period-one motion, and the period-2 motion will be observed over the range of 1.2≤ωh≤1.29. And as it is shown in Figure 5(c), the time domain response diagram shows a period motion, the phase plane diagram shows two closed circles, and the Poincaré map shows two unique points. Figure 5 presents time series, phase plane diagram, and Poincaré maps of Hopf bifurcation points at ωh=0.4, 0.7, 1.2, 1.3, 1.6, and 1.63. It can be seen that the motions of the one-way clutch, two-shaft assemblies system undergo period-1 motion, chaotic motion, period-2 motion, quasi-periodic motion, and period-7 motion.Figure 5
The time series ofz2, phase plane diagram, and Poincaré map under Fah=0.5, FaT=0, bh=0, b1=b2=bc=0.01, k^=5, Kc=5Kco, and Fm=1 and different excitation frequency with them; (a) ωh=0.4, (b) ωh=0.7, (c) ωh=1.2, (d) ωh=1.3, (e) ωh=1.6, and (f) ωh=1.63.
(a)
(b)
(c)
(d)
(e)
(f)Figure6 presents the bifurcation diagram of the system response of dimensionless dynamic transmission error between gear 2 and the inner ring, which is essentially used to study the dynamic windup angle of one-way clutch, using excitation frequency ωh as control parameter. It can be observed that the relative reverse motion between gear 2 and the inner ring exhibits rather complex phenomenon with the increase of the excitation frequency ωh. There exists some “frequency hopping” at ωh=0.26, 0.29, and 0.61, where the system undergoes the transient transition between period-one motion. Only period-one motion can be observed at low values of the excitation frequency, that is, ωh<0.69, or at high values of the excitation frequency, that is, ωh>1.68. The system performs period-2 motion within the range of 1.19≤ωh≤1.29. The system exhibits chaotic motions within the range of 0.7≤ωh≤0.88, 0.93≤ωh≤1.14, and 1.3≤ωh≤1.68. Unlike the transverse motion of gear 2, the relative reverse motion between gear 2 and the inner ring exhibits stable period-1 motion within the range of 0.89≤ωh≤0.92, 1.16≤ωh≤1.18. The system goes into chaotic motions through the quasi-periodic route around excitation frequency points ωh=0.71, 0.931, and 1.31 and goes into chaotic motions through the quasi-periodic route around excitation frequency points ωh=0.881, 1.141, and 1.681.Figure 6
γ
2 versus ωh, plots for Fah=0.5, FaT=0, bh=0, b1=b2=bc=0.01, k^=5, Kc=5Kco, and Fm=1.Figure7 shows the time series, phase plane diagram, and Poincaré maps of Hopf bifurcation points at ωh=0.4, 0.7, 1.2, 1.3, 1.6, and 1.63. It can be founded that the system behavior undergoes period-1 motion, chaotic motion, period-2 motion, quasi-periodic motion, and period-7 motion with the increase of the excitation frequency ωh.Figure 7
The time series ofγ2, phase plane diagram, and Poincaré map under Fah=0.5, FaT=0, bh=0, b1=b2=bc=0.01, k^=5, Kc=5Kco, and Fm=1 and different excitation frequency with them; (a) ωh=0.4, (b) ωh=0.7, (c) ωh=1.2, (d) ωh=1.3, (e) ωh=1.6, and (f) ωh=1.63.
(a)
(b)
(c)
(d)
(e)
(f)


### 4.2. Effect of Torsional Stiffness of One-Way Clutch

Figure8 presents the bifurcation diagram of the system response of dimensionless dynamic transmission error between gear 1 and gear 2, using the torsional stiffness of one-way clutch Kc as control parameter. It can be observed from the figure that, when the torsional stiffness is changed, there exist kinds of motion forms, such as periodic-1 motion, periodic-2 motion, periodic-4 motion, periodic-8 motion, and chaotic motion. In order to accurately distinguish the motion form, the time domain response diagrams, phase plane diagrams, and Poincaré maps are listed to illustrate the dynamic characteristics of the system in Figure 9. When Kc is less than 276000 N·m/rad, only period-one motion can be observed. The system motion state will be changed from periodic-1 motion to periodic-2 motion, when the torsional stiffness of one-way clutch is increased to 276000 N·m/rad, and the state is continued until 331000 N·m/rad. The system performs the transient periodic-4 motion within the range of 331000 N·m/rad ≤ Kc ≤ 341000 N·m/rad and exhibits chaotic motions within the range of 341000 N·m/rad ≤ Kc ≤ 676000 N·m/rad, just as shown in Figure 9(c). It can be seen that the time domain response diagram shows a nonperiodic motion, the Poincaré map shows many discrete points, and the phase plane diagram is disorder. With the increasing of the torsional stiffness, the system leaves chaos and goes into periodic-8 motion from 676000 N·m/rad to 696000 N·m/rad. When the torsional stiffness is greater than 696000 N·m/rad, the system goes into another region of chaotic motion again until 811000 N·m/rad. After undergoing the chaotic region, the system leaves chaotic motion and goes into periodic-1 motion again within the range of 811000 N·m/rad ≤ Kc ≤ 951000 N·m/rad. When Kc is higher than 951000 N·m/rad, the system leaves periodic-1 motion and returns into chaotic motion again.Figure 8
γ
1 versus Kc, plots for Fah=0.5, FaT=0, bh=0, b1=b2=bc=0.01, k^=5, ωh=1.25, and Fm=1.Figure 9
The time series ofγ1, phase plane diagram, and Poincaré map under Fah=0.5, FaT=0, bh=0, b1=b2=bc=0.01, k^=5, ωh=1.25, and Fm=1 and different torsional stiffness of one-way clutch with them; (a) Kc = 300000 N⋅m/rad, (b) Kc = 341000 N⋅m/rad, (c) Kc = 600000 N⋅m/rad, and (d) Kc = 688000 N⋅m/rad.
(a)
(b)
(c)
(d)


## 4.1. Effect of Excitation Frequency

The ratio of the bearing stiffness to average gear mesh stiffnessk^ is introduced to represent the change of support stiffness. The frequency parameter ωh is assigned as control parameters. The bifurcation diagram of the dimensionless transverse displacement of gear 2 with excitation frequency ωh change is presented in Figure 4.Figure 4
z
2 versus ωh, plots for Fah=0.5, FaT=0, bh=0, b1=b2=bc=0.01, k^=5, Kc=5Kco, and Fm=1.It is well known that the jump phenomenon usually exists in the dynamic response of nonlinear system. It could be seen clearly in Figure4 that there is a phenomenon of jumping at value of ωh=0.61. The one-way clutch, two-shaft assemblies system mainly exhibits period-one motion at low values or at high excitation frequency: that is, ωh<0.69 or ωh>1.66. Just as shown in Figure 5(a), the time domain response diagram is a sine wave, the phase plane diagram is only one closed circle, and the Poincaré map is a single point. It is observed that the system goes into chaotic region within the range of 0.7≤ωh≤0.87,0.94≤ωh≤1.13, and 1.3≤ωh≤1.59. It can be seen in Figure 5(b) that the time domain response diagram shows a nonperiodic motion, the phase plane diagram is disorder, and the Poincaré map shows many discrete points. As ωh is increased from 0.88 to 0.93 and from 1.14 to 1.19, the system exhibits unstable period-one motion, and the period-2 motion will be observed over the range of 1.2≤ωh≤1.29. And as it is shown in Figure 5(c), the time domain response diagram shows a period motion, the phase plane diagram shows two closed circles, and the Poincaré map shows two unique points. Figure 5 presents time series, phase plane diagram, and Poincaré maps of Hopf bifurcation points at ωh=0.4, 0.7, 1.2, 1.3, 1.6, and 1.63. It can be seen that the motions of the one-way clutch, two-shaft assemblies system undergo period-1 motion, chaotic motion, period-2 motion, quasi-periodic motion, and period-7 motion.Figure 5
The time series ofz2, phase plane diagram, and Poincaré map under Fah=0.5, FaT=0, bh=0, b1=b2=bc=0.01, k^=5, Kc=5Kco, and Fm=1 and different excitation frequency with them; (a) ωh=0.4, (b) ωh=0.7, (c) ωh=1.2, (d) ωh=1.3, (e) ωh=1.6, and (f) ωh=1.63.
(a)
(b)
(c)
(d)
(e)
(f)Figure6 presents the bifurcation diagram of the system response of dimensionless dynamic transmission error between gear 2 and the inner ring, which is essentially used to study the dynamic windup angle of one-way clutch, using excitation frequency ωh as control parameter. It can be observed that the relative reverse motion between gear 2 and the inner ring exhibits rather complex phenomenon with the increase of the excitation frequency ωh. There exists some “frequency hopping” at ωh=0.26, 0.29, and 0.61, where the system undergoes the transient transition between period-one motion. Only period-one motion can be observed at low values of the excitation frequency, that is, ωh<0.69, or at high values of the excitation frequency, that is, ωh>1.68. The system performs period-2 motion within the range of 1.19≤ωh≤1.29. The system exhibits chaotic motions within the range of 0.7≤ωh≤0.88, 0.93≤ωh≤1.14, and 1.3≤ωh≤1.68. Unlike the transverse motion of gear 2, the relative reverse motion between gear 2 and the inner ring exhibits stable period-1 motion within the range of 0.89≤ωh≤0.92, 1.16≤ωh≤1.18. The system goes into chaotic motions through the quasi-periodic route around excitation frequency points ωh=0.71, 0.931, and 1.31 and goes into chaotic motions through the quasi-periodic route around excitation frequency points ωh=0.881, 1.141, and 1.681.Figure 6
γ
2 versus ωh, plots for Fah=0.5, FaT=0, bh=0, b1=b2=bc=0.01, k^=5, Kc=5Kco, and Fm=1.Figure7 shows the time series, phase plane diagram, and Poincaré maps of Hopf bifurcation points at ωh=0.4, 0.7, 1.2, 1.3, 1.6, and 1.63. It can be founded that the system behavior undergoes period-1 motion, chaotic motion, period-2 motion, quasi-periodic motion, and period-7 motion with the increase of the excitation frequency ωh.Figure 7
The time series ofγ2, phase plane diagram, and Poincaré map under Fah=0.5, FaT=0, bh=0, b1=b2=bc=0.01, k^=5, Kc=5Kco, and Fm=1 and different excitation frequency with them; (a) ωh=0.4, (b) ωh=0.7, (c) ωh=1.2, (d) ωh=1.3, (e) ωh=1.6, and (f) ωh=1.63.
(a)
(b)
(c)
(d)
(e)
(f)


## 4.2. Effect of Torsional Stiffness of One-Way Clutch

Figure8 presents the bifurcation diagram of the system response of dimensionless dynamic transmission error between gear 1 and gear 2, using the torsional stiffness of one-way clutch Kc as control parameter. It can be observed from the figure that, when the torsional stiffness is changed, there exist kinds of motion forms, such as periodic-1 motion, periodic-2 motion, periodic-4 motion, periodic-8 motion, and chaotic motion. In order to accurately distinguish the motion form, the time domain response diagrams, phase plane diagrams, and Poincaré maps are listed to illustrate the dynamic characteristics of the system in Figure 9. When Kc is less than 276000 N·m/rad, only period-one motion can be observed. The system motion state will be changed from periodic-1 motion to periodic-2 motion, when the torsional stiffness of one-way clutch is increased to 276000 N·m/rad, and the state is continued until 331000 N·m/rad. The system performs the transient periodic-4 motion within the range of 331000 N·m/rad ≤ Kc ≤ 341000 N·m/rad and exhibits chaotic motions within the range of 341000 N·m/rad ≤ Kc ≤ 676000 N·m/rad, just as shown in Figure 9(c). It can be seen that the time domain response diagram shows a nonperiodic motion, the Poincaré map shows many discrete points, and the phase plane diagram is disorder. With the increasing of the torsional stiffness, the system leaves chaos and goes into periodic-8 motion from 676000 N·m/rad to 696000 N·m/rad. When the torsional stiffness is greater than 696000 N·m/rad, the system goes into another region of chaotic motion again until 811000 N·m/rad. After undergoing the chaotic region, the system leaves chaotic motion and goes into periodic-1 motion again within the range of 811000 N·m/rad ≤ Kc ≤ 951000 N·m/rad. When Kc is higher than 951000 N·m/rad, the system leaves periodic-1 motion and returns into chaotic motion again.Figure 8
γ
1 versus Kc, plots for Fah=0.5, FaT=0, bh=0, b1=b2=bc=0.01, k^=5, ωh=1.25, and Fm=1.Figure 9
The time series ofγ1, phase plane diagram, and Poincaré map under Fah=0.5, FaT=0, bh=0, b1=b2=bc=0.01, k^=5, ωh=1.25, and Fm=1 and different torsional stiffness of one-way clutch with them; (a) Kc = 300000 N⋅m/rad, (b) Kc = 341000 N⋅m/rad, (c) Kc = 600000 N⋅m/rad, and (d) Kc = 688000 N⋅m/rad.
(a)
(b)
(c)
(d)


## 5. Summary and Conclusions

A general multi-degree-of-freedom nonlinear dynamic model of one-way clutch, two-shaft assembly based on lumped parameter method (LPM) is proposed in this paper. The gear backlash, bearing backlash, a nonlinear mesh stiffnesskh(t), and mesh damping ch are considered in the model. The one-way clutch is modeled as a piecewise nonlinear spring with discontinuous stiffness that separates gear 2 into two degrees of freedom. The dynamic response of one-way clutch, two-shaft assemblies is studied over a wide range of frequencies. Speed sweeping is used to check the effect of one-way clutch on the nonlinear behavior of gear transmission system, by applying variable step-size Runge-Kutta method.As mentioned above, with the increasing of excitation frequency and torsional stiffness of one-way clutch, the system of one-way clutch, two-shaft assembly is prone to involve periodic-1 motion, multiperiodic motion, and chaotic motion. Large torsional stiffness of one-way clutch can increase the impact and lead to instability in the system. To reduce the impact and increase stability, the excitation frequency of the engine and the torsional stiffness of the one-way clutch are taken into account to be based on matching and optimization principles for the system of one-way clutch, two-shaft assembly.



---

*Source: 8621514-2017-08-20.xml*