Shock and Vibration
Volume 2021 (2021), Article ID 6687467, 20 pages
https://doi.org/10.1155/2021/6687467

Research Article

Nonlinear Frequency Response Analysis of Double-helical Gear Pair Based on the Incremental Harmonic Balance Method

Hao Dong, Libang Wang, Haoqin Zhang, and Xiao-long Zhao

School of Mechatronic Engineering, Xi’an Technological University, Xi’an 710021, China

Correspondence should be addressed to Hao Dong; donghaofane@sina.com

Received 30 November 2020; Revised 9 March 2021; Accepted 25 March 2021; Published 16 April 2021

Academic Editor: Radoslaw Zimroz

Copyright © 2021 Hao Dong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The torsional dynamic model of double-helical gear pair considering time-varying meshing stiffness, constant backlash, dynamic backlash, static transmission error, and external dynamic excitation was established. The frequency response characteristics of the system under constant and dynamic backlashes were solved by the incremental harmonic balance method, and the results were further verified by the numerical integration method. At the same time, the influence of time-varying meshing stiffness, damping, static transmission error, and external load excitation on the amplitude frequency characteristics of the system was analyzed. The results show that there is not only main harmonic response but also superharmonic response in the system. The time-varying meshing stiffness and static transmission error can stimulate the amplitude frequency response of the system, while the damping can restrain the amplitude frequency response of the system. Changing the external load excitation has little effect on the amplitude frequency response state change of the system. Compared with the constant backlash, increasing the dynamic backlash amplitude can further control the nonlinear vibration of the gear system.

1. Introduction

Double-helical gear is widely used in vehicle, ship, and other heavy machinery transmission systems because of its strong bearing capacity, compact structure, and axial force cancellation. Its nonlinear vibration characteristics will directly affect the stability and reliability of the transmission system. Comparin and Singh [1] established a three-degree of freedom spur gear transmission system considering the random disturbance of low frequency external excitation caused by torque fluctuation and studied the influence of different backlash, gear damping ratio, and random disturbance on the dynamic behavior of the gear transmission system. Kahraman and Singh [2–4] analyzed the nonlinear dynamic frequency response characteristics of spur gear pair under internal excitation by the numerical method and analytical method and found that there are many nonlinear phenomena such as the jump phenomenon and multivalue solution in the vibration response of spur gear. Shen et al. [5, 6] established the torsional vibration model of spur gear pair and gave the unified form of nonlinear analytical solution of the gear system including time-varying meshing stiffness and backlash. Wang et al. [7] established the nonlinear dynamic model of spur gear pair with fault and obtained the phase diagram and bifurcation diagram of the single-stage gear system but did not study the nonlinear amplitude frequency characteristics of the gear system with or without faults. Li et al. [8] established the pure torsion dynamic model of spur gear pair and studied the influence of system parameters on the nonlinear dynamic characteristics of the spur gear system. Zhao et al. [9] analyzed the influence of gear side backlash, eccentricity, and tooth surface friction on the nonlinear vibration of the high-speed gearbox system of wind turbine. Sun et al. [10] developed the application of the single degree of freedom analytical harmonic balance method, derived the calculation formula of the analytic harmonic balance method for differential equations of the planetary gear system, and obtained the nonlinear characteristics of the planetary gear transmission system. Yuan et al. [11] studied the influence of cumulative pitch error on the dynamic characteristics of the double-helical gear system. Wang et al. [12, 13] studied the vibration characteristics of the double-helical gear transmission system under different loads and different coincidence degrees. In reference [9, 11–13], the vibration problem of the gear system is studied, but the strong nonlinear characteristics caused by backlash are not further explained. In fact, it still belongs to the category of linear dynamics. In addition, many scholars have performed a lot of research on the load sharing characteristics of double-helical gears [14–17]. Mo et al. [18] analyzed the nonlinear dynamic characteristics of double frequency problem of the double-helical gear closed differential planetary transmission system under no-load or light load conditions. Shin and Palazzolo [19] proposed a new method for modeling and analyzing the geared rotor-bearing system. Dong et al. [20] established a 6-DOF spur gear dynamic model including time-varying meshing stiffness, backlash, support backlash, and tooth surface friction. Yavuz et al. [21] first proposed a nonlinear dynamic model of parallel shaft gear and intersecting shaft gear and solved the nonlinear algebraic equations by the multiorder harmonic balance method combined with continuous Fourier transform, but the model is complex and not widely applicable. Xiang et al. [22] established a torsional nonlinear dynamic model of multistage planetary gear considering time-varying meshing stiffness, comprehensive gear error, and piecewise backlash nonlinearities and studied the effects of excitation frequency, backlash, and damping on bifurcation characteristics. Yang et al. [23] established a single degree of freedom spur gear pair model including time-varying mesh stiffness and static transmission error. The nonlinear frequency response characteristics of a spur gear system were studied. Lian et al. [24] established a multiple degrees of freedom nonlinear dynamic model of a gear pair with time-varying mesh stiffness, mesh damping, backlash, dynamic transmission errors, and radial clearance of ball bearing by using the mass centralized method. The effects of main parameters such as exciting frequency, radial clearance of ball bearing, and meshing stiffness ratio on chaos and bifurcation of the system were studied. Li and Peng [25] proposed a gear nonlinear model of multidegree of freedom considering the characteristics of gear tooth. Fractal theory was used to calculate the nonlinear backlash from the tribological aspect. Li et al. [26] studied the influence of the fractal backlash on the dynamic characteristics of the gear bearing system, and the fractal behavior of the backlash was described by W-M function. However, the influence of fractal on the nonlinear frequency response of gears is not studied in reference [25, 26].

To sum up, the research on the nonlinear characteristics of the gear transmission system at home and abroad mostly focused on the spur gear and spur planetary gear system. In the reports on the nonlinear dynamics of double-helical gear, there are few reports on the nonlinear frequency response characteristics of double-helical gear. Based on the above research, a torsional dynamic model of the double-helical gear pair system is established in this study. The time-varying meshing stiffness, static backlash, dynamic backlash, static transmission error, external dynamic excitation, and other factors are considered in the model. The influence of parameter changes on the nonlinear dynamic characteristics of the double-helical gear pair system is studied by using the incremental harmonic balance method. The research results can provide reference for vibration reduction analysis and structure optimization of double-helical gear.

2. Torsional Nonlinear Dynamic Model of Double-Helical Gear Transmission

As shown in Figure 1, the torsional nonlinear dynamic model of double-helical gear pair is shown. The time-varying meshing stiffness, viscoelastic damping, backlash, static transmission error, and external dynamic excitation are considered in the model.

Figure 1: Torsional nonlinear dynamic model of double-helical gear pair.

According to Newton’s second law, the torsional nonlinear dynamic equation of double-helical gear pair can be written as

In the above formula, Iij(i = p, .j = L, R) represents the polar moment of inertia on the equivalent nodes OpL, OpR, OgL, and OgR at the center of each helical gear, θij(i = p,; j = L, R) represents the angular rotation displacement of the equivalent node at the center of each helical gear, rij(i = p, ; j = L, R) represents the radius of the base circle of each helical gear, βj(j = L, R) represents the helix angle of the left and right helical gear pairs, kmj, cmj(j = L, R) represents the meshing stiffness and damping of the left and right helical gear pairs, δj(j = L, R) represents the relative meshing displacement of the left and right helical gear pairs, ej(j = L, R) represents the static transmission error of the left and right helical gear pairs, kit,cit(i = p, ) represents the torsional stiffness and torsional damping of the intermediate shaft of double-helical gear, and TpL(t), TpR(t), TgL(t), TgR(t) represent the input torque and output torque acting on the driving helical gear and the driven helical gear, respectively.

Each moment can be expressed as follows:

where TpL, TpR, TgL, and TgR are the external nominal moments, and TpLa(t), TpRa(t), TgLa(t), and TgRa(t) are the external dynamic moments.

The four degrees of freedom θpL, θpR, θgL, and θgR in the nonlinear torsional dynamic equation of double-helical gear transmission are not independent of each other. The relative displacement method is used to eliminate the rigid body displacement of the system to solve the differential equation. The absolute torsional angular displacement θpL, θpR, θgL, and θgR on the original mass nodes is replaced by the relative meshing displacement δL, δR along the meshing line direction of the left and right helical gears and the relative torsional line displacement γP, γG between the coaxial adjacent helical gears.

The relative meshing displacement of the left and right helical gear pairs along the direction of the meshing line is expressed as follows:

The expression of relative torsional line displacement of adjacent coaxial helical gears is as follows:

The actual meshing displacement of double-helical gear along the meshing line is expressed as a piecewise function as follows:

where 2b is the length of backlash. For the actual double-helical gear pair transmission, with the difference between the relative meshing displacement along the meshing line direction and the backlash, that is, the actual meshing displacement along the meshing line direction, the double-helical gear pair transmission may be in three different working states, namely, no impact, single-side impact, and double-side impact [22, 23, 27], as shown in Figure 2.

Figure 2: Meshing impact state of double-helical gear pair.

In the existing research, the backlash function is expressed as piecewise linear. Generally, the backlash of gear teeth is obtained by static measurement, and the change of backlash in the process of gear impact is not considered, so the simulation results cannot truly reflect the transmission characteristics of the gear system. In order to reveal the influence of the actual tooth surface characteristics on the vibration and noise of the gear transmission system, the constant backlash model and dynamic backlash model will be considered simultaneously in this study.

The constant backlash model takes the backlash as a constant according to the traditional modeling method. The expression is as follows:

The dynamic backlash model assumes that the backlash is in the form of harmonics. The expression is as follows:

where b₀ is the equivalent backlash, W is the backlash variation describing the wear of gear tooth, and bl is the amplitude of dynamic backlash representing the deformation of the supporting bearing and axle. By substituting equations (3), (4), and (9) into (2), the dynamic equation (8) is obtained:

where mc is the equivalent mass of the left and right helical gear pairs, and fa is the fluctuation of external load excitation, which are, respectively, expressed as follows:

where fl is the fluctuation coefficient. The time-varying meshing stiffness function and static transmission error function are expressed in the form of Fourier series [2–6, 8, 28, 29]:

where kmean is the average stiffness, εl is the fluctuation coefficient reflecting the change of stiffness, and el is the fluctuation coefficient reflecting the change of static transmission error. When the excitation of static transmission error is maximum, the time-varying meshing stiffness should be the minimum, that is, the two have inverse phase relationship. The time-varying meshing stiffness curve is shown in Figure 3.

Figure 3: Time-varying mesh stiffness curve.

The meshing damping cmj and torsional damping cit are calculated by the following formula:

where ζ is the damping ratio of the system.

According to the knowledge of material mechanics, the calculation formula of torsional stiffness of intermediate shaft section is as follows:

where G is the shear elastic modulus of the gear material, J is the polar moment of inertia of the cross-section, l is the length of the intermediate shaft segment, and D and d are the outer diameter and inner diameter of the shaft, respectively.

Equation (8) is dimensionless, , , and , where bc is the nominal displacement scale, Δi, Δj is the normalized vibration displacement, ωn is the natural frequency, t is the normalized time, and Ω is the normalized meshing frequency. In this study, the derivative of time t in equations (8) is transformed into the derivative of dimensionless time t, and finally, the dimensional unified dynamic equation is obtained:

The stiffness term in the above formula is as follows:

The static transmission error term in the excitation vector is as follows:

The dynamic excitation of external torque is as follows:

For the piecewise backlash function, if bc = b₀ is taken in the dimensionless process, the two kinds of backlash models are transformed into

Equation (20) is written in the matrix form as follows:

where is the dimensionless mass matrix, is the dimensionless damping matrix, is the dimensionless stiffness matrix, is the dimensionless excitation vector, and is the dimensionless displacement vector. Limited to space, only the stiffness matrix is listed as follows:

3. Approximate Solution Based on the Incremental Harmonic Balance Method

By letting a new time scale , equation (22) can be represented as

The Newton–Raphson algorithm is used to solve the piecewise linear difference system represented by equation (24), and the solution of the differential equation is obtained as follows:

where is the approximate solution of equation (24), and is the incremental equation. Let any term in the n-order approximate solution of formula (25) and any term in the incremental equation are, respectively,

where are the Fourier coefficients, and N is the number of terms of the Fourier series. The nonlinear piecewise linear function of backlash may be expressed by a first order Taylor expansion as

By substituting equations (26) and equation (29) into equation (24), Taylor series expansion is carried out, and the high-order term is omitted, and equation (24) becomes linearized as

where

Order

By using the Galerkin process [5, 8, 30–32], the left and right ends of equation (29) are multiplied by cosiτ, sin(iτ) (i = 0, 1, 2, 3, …, N). The linear equations are obtained by integrating between 0 and 2.

Equation (36) is the iterative calculation formula of the steady-state periodic response of the nonlinear system with as unknown quantity, which is derived by using the incremental harmonic balance method.

According to equation (36), the approximate solution with arbitrary precision can be obtained, that is, given the initial value of iteration, is obtained; thus, the initial value A of the next iteration can be obtained. When the accuracy requirement is met, the iteration will be terminated.

4. Nonlinear Frequency Response Characteristics of the Gear System

4.1. Comparison between Analytical Solution and Numerical Solution

The design parameters of double-helical gear pair are given in Table 1.

Table 1: Basic parameters of the double-helical gear pair.

Parameter
Pinion
Gear

Number of teeth
25
75

Helix angle (°)
24.43

Normal pressure angle (°)
20

Normal module (mm)
6

Tooth width (mm)
55 × 2

Moment of inertia (kg·m²)
0.012
0.7279

Mass (kg)
7.11
64.75

Input torque (Nm)
1400

Meshing stiffness (N/m)
1.8456 × 10⁹

Length of intermediate shaft section (mm)
58

Outer diameter of intermediate shaft section (mm)
92
260

Inner diameter of intermediate shaft section (mm)
55
70

In order to fully study the nonlinear characteristics of double-helical gear under the coupling effect of time-varying meshing stiffness and backlash, the amplitude frequency characteristic curves of the system under four different models are obtained, and the results are shown in Figure 4. The basic system parameters are L = 3, N = 11, F₁ = 1.3534, ζ = 0.045, ɛ₁ = 0.16, ɛ₂ = 0.07, ɛ₃ = 0.03, e₁ = 0.03, e₂ = 0.01, e₃ = 0.01, f₁ = 0.06, f₂ = 0.03, f₃ = 0.01, W = 0, b₁ = 0.05, b₂ = 0.02, and b₃ = 0.01.

(a)

(b)

(c)

(d)

Figure 4: Amplitude frequency characteristic curve of the system under different models. (a) Constant backlash model under combined internal and external excitation, (b) constant backlash model under external excitation, (c) constant backlash model under internal excitation, and (d) dynamic backlash model under combined internal and external excitation.

In Figure 4, the ordinate represents the maximum ΔLmax and minimum ΔLmin nonlinear vibration displacement of the double-helical gear system along the meshing line, and the abscissa represents the dimensionalized frequency Ω. Four different models are, respectively, the constant backlash model under combined internal and external dynamic excitations as shown in Figure 4(a), the constant backlash model under external dynamic excitation as shown in Figure 4(b), the constant backlash model under internal dynamic excitation as shown in Figure 4(c), and the dynamic backlash model under combined internal and external dynamic excitations as shown in Figure 4(d). The solid line represents the result obtained by incremental harmonic balance method, and the circle represents the result obtained by the Runge–Kutta method with variable step size of fourth and fifth order. Four different models are as follows:

Figure 4 shows that there are multivalue solutions and jumping phenomena in the double-helical gear system under four different models. The multivalue solution interval of the constant backlash model is 0.78–0.86 and that of the dynamic backlash model is 0.83–0.84; the system includes no impact state and single-side impact state, no double-side impact state is found, and tooth back separation does not occur. In the high frequency range of dimensionless frequency from 1.794 to 1.803, transient instability occurs in the numerical method but not in the analytical method. At the same time, it can be found that the nonlinear dynamic characteristics under combined internal and external excitations are not simply obtained by linear superposition of internal and external excitations, respectively.

The multivalued properties and jump phenomena are corresponding with no impact motion and single-side impact motion. Figure 5 shows the steady-state response of the gear pair system when Ω = 0.83.

(a)

(b)

Figure 5: Steady-state response of the gear pair system with Ω = 0.83. (a) No impact and (b) single-side impact.

The basic system parameters are L = 3, N = 11, F₁ = 1.3534, ζ = 0.045, ɛ₁ = 0.16, ɛ₂ = 0.07, ɛ₃ = 0.03, e₁ = 0.03, e₂ = 0.01, e₃ = 0.01, f₁ = 0.06, f₂ = 0.03, f₃ = 0.01, and W = 0. The phase diagram and time history diagram of main resonance response and superharmonic response of the double-helical gear system are obtained by using the incremental harmonic balance method. Meanwhile, the Runge–Kutta numerical method is used for comparison and verification. The results are shown in Figure 6.

(a)

(b)

(c)

Figure 6: Phase diagram and time history diagram of the system at different frequencies. (a) Ω = 1, (b) Ω = ½, and (c) Ω = 1/3.

Figure 6 shows that when Ω = 1, the phase diagram is approximate ellipse, the system is a single period harmonic response, and the gear is in the state of no impact; when Ω = 1/2 and Ω = 1/3, the phase diagram is a noncircular closed curve, the system response is superharmonic response, and the meshing state of gear teeth is not changed. The generation of the superharmonic response solution of the system is caused by the high-order harmonic component in the time-varying meshing stiffness and dynamic excitation. The results of the incremental harmonic balance method are in good agreement with the numerical simulation results.

4.2. Influence of Backlash on Frequency Response Characteristics of the System

The basic system parameters are L = 3, N = 11, F₁ = 1.3534, ζ = 0.045, ɛ₁ = 0.16, ɛ₂ = 0.07, ɛ₃ = 0.03, e₁ = 0.03, e₂ = 0.01, e₃ = 0.01, f₁ = 0.06, f₂ = 0.03, f₃ = 0.01, W = 0, b₁ = 0.05, b₂ = 0.02, and b₃ = 0.01. The amplitude frequency characteristic curve of double-helical gear pair under the linear time-varying model (zero backlash) and nonlinear time-varying model (with backlash) is shown in Figure 7.

(a)

(b)

Figure 7: Influence of backlash on vibration response of the gear system. (a) Constant backlash models and (b) dynamic backlash models.

From the comparison between the linear time-varying model (zero backlash) and the nonlinear time-varying model (with backlash) in Figure 7, it can be seen that the existence of backlash nonlinearity leads to the nonlinear dynamic characteristics such as multivalued solution and amplitude jump of the double-helical gear pair system.

The basic system parameters are L = 3, N = 11, F₁ = 1.3534, ζ = 0.045, ɛ₁ = 0.16, ɛ₂ = 0.07, ɛ₃ = 0.03, e₁ = 0.03, e₂ = 0.01, e₃ = 0.01, f₁ = 0.06, f₂ = 0.03, and f₃ = 0.01. The influence of different constant backlash and dynamic backlash amplitude on the amplitude frequency characteristics of the system is shown in Figure 8.

(a)

(b)

Figure 8: Amplitude frequency characteristic curve of the system under two kinds of backlash models. (a) Constant backlash models and (b) dynamic backlash models.

For constant backlash, five models of constant backlash are considered: (1) W = −0.6; (2) W = −0.3; (3) W = 0; (4) W = 0.3; and (5) W = 0.6. It can be seen from Figure 8(a), with the increase of the constant backlash, the vibration displacement amplitude of the gear system increases, the multivalue solution region is wider, and the jumping phenomenon still exists. For dynamic backlash, four kinds of dynamic backlash models are considered: (1) b₁ = 0, b₂ = 0, b₃ = 0; (2) b₁ = 0.05, b₂ = 0.02, b₃ = 0.01; (3) b₁ = 0.10, b₂ = 0.07, b₃ = 0.03; and (4) b₁ = 0.15, b₂ = 0.08, b₃ = 0.04.

It can be seen from Figure 8(b) that with the increase of the dynamic backlash amplitude, the vibration displacement amplitude of the gear system slightly increases, the jumping phenomenon still exists, and the multivalue solution interval gradually decreases until it disappears. It can be seen that increasing the dynamic backlash can control the collision behavior of the gear system.

5. Influence of System Parameters on Amplitude Frequency Curve

5.1. Influence of Time-Varying Meshing Stiffness on Amplitude Frequency Curve of the System

The basic system parameters are L = 3, N = 11, F₁ = 1.3534, ζ = 0.045, e₁ = 0.03, e₂ = 0.01, e₃ = 0.01, f₁ = 0.06, f₂ = 0.03, f₃ = 0.01, W = 0, b₁ = 0.05, b₂ = 0.02, and b₃ = 0.01. The influence of time-varying meshing stiffness on the amplitude frequency curve of the system under the constant backlash model and dynamic backlash model is shown in Figure 9.

(a)

(b)

Figure 9: Influence of time-varying meshing stiffness on amplitude frequency characteristic curve of the system. (a) Constant backlash models and (b) dynamic backlash models.

In this study, three kinds of time-varying meshing stiffness fluctuation coefficients are considered: (1) ɛ₁ = 0.16, ɛ₂ = 0.07, ɛ₃ = 0.03; (2) ɛ₁ = 0.13 ɛ₂ = 0.04, ɛ₃ = 0.01; and (3) ɛ₁ = 0, ɛ₂ = 0, ɛ₃ = 0. It can be seen from Figures 9(a) and 9(b) that with the decrease of time-varying meshing stiffness fluctuation coefficient, the amplitude of amplitude frequency characteristic curve decreases, which is more obvious at lower frequencies; meanwhile, with the decrease of time-varying meshing stiffness fluctuation coefficient, the multivalue solution phenomenon and jump phenomenon of the system will gradually decrease or even disappear.

In addition, compared with the constant backlash model, the dynamic backlash model can control the nonlinear vibration of the system more effectively by reducing the fluctuation coefficient of the meshing stiffness.

5.2. Influence of Damping on Amplitude Frequency Curve of the System

The basic system parameters are L = 3, N = 11, F₁ = 1.3534, ɛ₁ = 0.16, ɛ₂ = 0.07, ɛ₃ = 0.03, e₁ = 0.03, e₂ = 0.01, e₃ = 0.01, f₁ = 0.06, f₂ = 0.03, f₃ = 0.01, W = 0, b₁ = 0.05, b₂ = 0.02, and b₃ = 0.01. The amplitude frequency characteristic curves of the system under different damping ratios are shown in Figure 10.

(a)

(b)

Figure 10: Influence of damping on amplitude frequency characteristic curve of the system. (a) Constant backlash models and (b) dynamic backlash models.

It can be seen from Figure 10 that damping can restrain the step of amplitude frequency characteristic curve, the emergence of multivalue solution interval, and the amplitude of nonlinear vibration displacement. Increasing the system damping can effectively control the nonlinear dynamic response of the gear system; compared with the constant backlash model, increasing the damping ratio under the dynamic backlash model can more effectively control the nonlinear vibration of the system.

5.3. Influence of Static Transmission Error on Amplitude Frequency Curve of the System

The basic system parameters are L = 3, N = 11, F₁ = 1.3534, ζ = 0.045, ɛ₁ = 0.16, ɛ₂ = 0.07, ɛ₃ = 0.03, f₁ = 0.06, f₂ = 0.03, f₃ = 0.01, W = 0, b₁ = 0.05, b₂ = 0.02, and b₃ = 0.01. The influence of static transmission error on the amplitude frequency curve of the system is shown in Figure 11.

(a)

(b)

Figure 11: Influence of static transmission error on amplitude frequency characteristic curve of the system. (a) Constant backlash models and (b) dynamic backlash models.

This study considers three kinds of static transmission error amplitudes: (1) e₁ = 0.1, e₂ = 0.05, e₃ = 0.03; (2) e₁ = 0.03, e₂ = 0.01, e₃ = 0.01; and (3) e₁ = 0, e₂ = 0, e₃ = 0. It can be seen from Figure 11 that the amplitude frequency characteristic curve of the system changes greatly with the change of the amplitude of the static transmission error, when e₁ = 0, e₂ = 0, and e₃ = 0, and the amplitude frequency curve still has a jump phenomenon and multivalue solution phenomenon, which indicates that the static transmission error in double-helical gear is not the main reason for the nonlinearity of the system. When e₁ = 0.1, e₂ = 0.05, e₃ = 0.03 and e₁ = 0.03, e₂ = 0.01, e₃ = 0.01, the larger the amplitude frequency curve of the system appears, and the larger the value of static transmission error amplitude is, the greater the amplitude frequency characteristic curve amplitude is.

5.4. Influence of Excitation Force Amplitude on Amplitude Frequency Curve of the System

The basic system parameters are L = 3, N = 11, ζ = 0.045, ɛ₁ = 0.16, ɛ₂ = 0.07, ɛ₃ = 0.03, e₁ = 0.03, e₂ = 0.01, e₃ = 0.01, f₁ = 0.06, f₂ = 0.03, f₃ = 0.01, W = 0, b₁ = 0.05, b₂ = 0.02, and b₃ = 0.01. In the constant backlash model and dynamic backlash model, the influence of excitation amplitude force on the amplitude frequency curve of the system is shown in Figure 12.

(a)

(b)

Figure 12: Influence of excitation force amplitude on amplitude frequency characteristic curve of the system. (a) Constant backlash models and (b) dynamic backlash models.

It can be seen from Figure 12 that with the increase of the excitation force amplitude, the amplitude frequency characteristic curve amplitude of the system increases, the jump phenomenon still exists, and the multivalue solution phenomenon disappears slowly. When F₁ = 3.3534, the amplitude frequency characteristic curve amplitude is the largest, while there is no multivalue solution interval. When F₁ = 1.3534, the amplitude of amplitude frequency characteristic curve is the smallest, and the range of unstable solution interval is the largest. The transition frequency of curve step basically does not change, and the results of the two backlash models have little difference. With the increase of the excitation amplitude, the amplitude frequency characteristic curve of the system still contains the no impact state and the single-side impact state, and the increase of the excitation amplitude has little effect on the change of the meshing impact state of the double-helical gear.

6. Conclusions

In this study, the torsional dynamic model of double-helical gear transmission is established. The factors such as time-varying meshing stiffness, static transmission error, constant backlash, dynamic backlash, and external dynamic excitation are considered in the model. The frequency response characteristic curve of the system is better solved by using the incremental harmonic balance method. It is found that the system has multivalue solution and jumping nonlinear characteristics, and the results are in good agreement with the numerical integration results.

The conclusions are as follows:

(1)There are main resonance response and superharmonic response in the double-helical gear system. The generation of the superharmonic response solution is caused by the high-order harmonic component in time-varying meshing stiffness and dynamic excitation.(2)With the increase of the constant backlash, the vibration displacement amplitude of the gear system increases, and the multivalue solution interval becomes larger; with the increase of the dynamic backlash amplitude, the vibration displacement amplitude of the gear system increases slightly, and the multivalue solution range decreases or even disappears.(3)The results show that the time-varying meshing stiffness and static transmission error can stimulate the amplitude frequency curve of the system, and the damping can restrain the amplitude frequency curve of the system. Compared with the constant backlash model, the vibration of the double-helical gear system can be further controlled by changing the system parameters in the dynamic backlash model.(4)The excitation amplitude has little effect on the change of the meshing impact state of the double-helical gear. Under the influence of excitation amplitude, the transition frequency of curve step has no change, and the results of the constant backlash model and dynamic backlash model have little difference.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was funded by the National Natural Science Foundation of China (51705390), Innovation Capability Support Program of Shaanxi (2020KJXX-016), Scientific Research Program Funded by Shaanxi Provincial Education Department Program (20JC015), Principal foundation project of Xi’an Technological University (XGPY200201), and the Natural Science Foundation of Shaanxi Province (2021JM-428).

References

R. J. Comparin and R. Singh, “Non-linear frequency response characteristics of an impact pair,” Journal of Sound and Vibration, vol. 134, no. 2, pp. 259–290, 2015.

A. Kahraman and R. Singh, “Non-linear dynamics of a spur gear pair,” Journal of Sound and Vibration, vol. 142, no. 1, pp. 49–75, 1990.

A. Kahraman and R. Singh, “Interactions between time-varying mesh stiffness and clearance non-linearities in a geared system,” Journal of Sound and Vibration, vol. 146, no. 1, pp. 135–156, 1991.

A. Al-Shyyab and A. Kahraman, “Non-linear dynamic analysis of a multi-mesh gear train using multi-term harmonic balance method, period-one motions,” Journal of Sound and Vibration, vol. 284, no. 1-2, pp. 151–172, 2005.

Y. J. Shen, S. P. Yang, C. Z. Peng et al., “Nonlinear dynamics of a spur gear pair with time-varying stiffness and backlash,” Journal of Low Frequency Noise Vibration and Active Control, vol. 23, no. 3, pp. 179–187, 2004.

S. P. Yang, Y. J. Shen, and X. D. Liu, “Nonlinear dynamics of gear system based on incremental harmonic balance method,” Journal of Sound and Vibration, vol. 3, no. 40–42, pp. 95–132, 2005, in Chinese.

Y. L. Wang, Y. J. Wan, F. W. Yin et al., “Analysis of nonlinear dynamics characteristic of spur gear pair system with fault,” Journal of Mechanical Transmission, vol. 43, no. 10, pp. 136–140, 2019, in Chinese.

Y. Li, T. Chen, and X. Wang, “Non-linear dynamics of gear pair with dynamic backlash subjected to combined internal and external periodic excitations,” Journal of Vibration and Control, vol. 22, no. 6, pp. 1693–1703, 2016.

X. Zhao, C. Z. Chen, J. Liu et al., “Nonlinear dynamic response analysis of high level of wind turbine gearbox transmission system considering eccentricity,” Journal of Solar Energy, vol. 41, no. 3, pp. 98–108, 2020, in Chinese.

T. Sun, Y. W. Shen, Z. M. Sun et al., “Study on nonlinear dynamic behavior of planetary gear train solution and dynamic behavior analysis,” Chinese Journal of Mechanical Engineering, vol. 38, no. 3, pp. 11–15, 2002, in Chinese.

B. Yuan, S. Chang, G. Liu et al., “Dynamic characteristics of a double helical gear system considering accumulation pitch error,” Journal of Sound and Vibration, vol. 39, no. 3, pp. 120–126, 2020, in Chinese.

F. Wang, Z. D. Fang, and S. J. Li, “Dynamic characteristics of a double helical gear under multi-load,” Journal of Sound and Vibration, vol. 32, no. 1, pp. 49–52, 2013, in Chinese.

F. Wang, Z. D. Fang, and S. J. Li, “Effect of contact ratio on dynamic behavior of a double-helical gear nonlinear system,” Journal of Sound and Vibration, vol. 33, no. 3, pp. 18–22, 2014, in Chinese.

S. Mo, T. Zhang, G. G. Jin et al., “Analytical investigation on load sharing characteristics of herringbone planetary gear train with flexible support and floating sun gear,” Mechanism and Machine Theory, vol. 144, no. 2, pp. 1–27, 2020.

M. Shuai, X. Y. Zong, Y. F. Zhi et al., “Analytical investigation on load-sharing characteristics for multi-power face gear split flow system,” Journal of Mechanical Engineering Science, vol. 234, no. 2, pp. 676–692, 2020.

S. Mo, J. Gong, G. Jin et al., “Precise modeling of complex tooth surface microtopography and multi-degree-of-freedom nonlinear friction dynamics for high-performance face gear,” Science Progress, vol. 103, no. 1, pp. 1–28, 2020.

Z. G. Wang and Y. C. Chen, “Design of a helical gear set with adequate linear tip-relief leading to improved static and dynamic characteristics,” Mechanism and Machine Theory, vol. 147, 2019.

W. C. Mo, Y. H. Jiao, Z. B. Chen et al., “Nonlinear dynamic analysis of a herringbone planetary gear system,” Journal of Harbin Engineering University, vol. 40, no. 10, pp. 1760–1766, 2019, in Chinese.

D. Shin and A. Palazzolo, “Nonlinear analysis of a geared rotor system supported by fluid film journal bearings,” Journal of Sound and Vibration, vol. 475, Article ID 115269, 2020.

H. Dong, J. Zhang, and L. Wang, “Study on bifurcation characteristics of multi-clearance bending torsional coupling gear transmission based on EHL,” Advances in Mechanical Engineering, vol. 12, no. 6, Article ID 168781402093750, 2020.

S. D. Yavuz, Z. B. Saribay, and E. Cigeroglu, “Nonlinear dynamic analysis of a drivetrain composed of spur, helical and spiral bevel gears,” Nonlinear Dynamics, vol. 100, no. 4, pp. 3145–3170, 2020.

L. Xiang, N. Gao, and A. Hu, “Dynamic analysis of a planetary gear system with multiple nonlinear parameters,” Journal of Computational and Applied Mathematics, vol. 327, pp. 325–340, 2017.

Y. Yang, L. Cao, H. Li, and Y. Dai, “Nonlinear dynamic response of a spur gear pair based on the modeling of periodic mesh stiffness and static transmission error,” Applied Mathematical Modelling, vol. 72, no. 8, pp. 444–469, 2019.

C. S. Lian, L. Wei, Q. W. Yu et al., “Nonlinear dynamic analysis and chaos control of multi-freedom semi-direct gear drive system in coal cutters,” Mechanical Systems and Signal Processing, vol. 116, pp. 62–77, 2018.

Z. Li and Z. Peng, “Nonlinear dynamic response of a multi-degree of freedom gear system dynamic model coupled with tooth surface characters, a case study on coal cutters,” Nonlinear Dynamics, vol. 84, no. 1, pp. 1–16, 2015.

X. P. Li, J. X. Mou, W. J. Pan et al., “Influence of fractal backlash on dynamic behavior of gear-bearing system,” Journal of Mechanical Engineering, vol. 54, no. 9, pp. 153–160, 2018, in Chinese.

L. Sheng, W. Li, Y. Wang, M. Fan, and X. Yang, “Nonlinear dynamic analysis and chaos control of multi-freedom semi-direct gear drive system in coal cutters,” Mechanical Systems and Signal Processing, vol. 116, pp. 62–77, 2019.

M. Yin, G. Chen, and H. Su, “Theoretical and experimental studies on dynamics of double-helical gear system supported by journal bearings,” Advances in Mechanical Engineering, vol. 8, no. 5, 2016.

J. Dong, S. Wang, H. Lin et al., “Dynamic modeling of double-helical gear with Timoshenko beam theory and experiment verification,” Advances in Mechanical Engineering, vol. 8, no. 5, 2016.

T. C. Kim, T. E. Rook, and R. Singh, “Super- and sub-harmonic response calculations for a torsional system with clearance nonlinearity using the harmonic balance method,” Journal of Sound & Vibration, vol. 281, no. 3–5, pp. 965–993, 2005.

W. Zhang and Q. Ding, “Torsion vibration and parametric instability analysis of a spur gear system with time-varying and square nonlinearities,” International Journal of Applied Mechanics, vol. 6, no. 1, pp. 26–77, 2014.

D. Y. Siar, B. S. Zihni, and C. Ender, “Nonlinear time-varying dynamic analysis of a spiral bevel geared system,” Nonlinear Dynamics, vol. 92, no. 4, pp. 1901–1919, 2018.