Bending-torsional-axial-pendular nonlinear dynamic modeling and frequency response analysis of a marine double-helical gear drive system considering backlash

Abstract

The double-helical gear system was widely used in ship transmission. In order to study the influence of backlash on the nonlinear frequency response characteristics of marine double-helical gear system, according to the structural characteristics of double-helical gear transmission, considering the time-varying meshing stiffness, backlash, damping, comprehensive transmission error, external load excitation, and other factors, a three-dimensional bending-torsional-axial-pendular coupling nonlinear dynamic modeling and dynamic differential equation of 24-DOF double-helical gear transmission system were established. The Runge–Kutta numerical method was used to analyze the influence of backlash, time-varying meshing stiffness, damping, error and external load excitation on the amplitude frequency characteristics. The results show that the backlash can cause the runout of the double-helical gear system, and the system has first harmonic and second harmonic response. With the increase of backlash, the amplitude of the system increases and the jumping phenomenon remains unchanged. The amplitude frequency response of the system is stimulated by time-varying meshing stiffness and comprehensive transmission error, and restrained by damping and external load excitation. The vibration displacement amplitude of the system increases with the increase of vibration displacement and has little effect on the state change of the system. The vibration test of double-helical gear is carried out. The frequency response components obtained by numerical simulation are basically consistent with the experimental results, which proves the correctness of the theoretical calculation. It provides a technical basis for the study of vibration and noise reduction performance of double-helical gear.

Introduction

Double-helical gear has the advantages of compact structure, large bearing capacity, stable operation, and can overcome large axial force. It is widely used in marine transmission device. The noise produced by gear vibration will not only affect the comfort of the cabin, but also threaten the safety of the ship. ¹

Many scholars at home and abroad have done a lot of research on gear vibration. Kahraman^2–4 used numerical method and analytical method to analyze the nonlinear dynamic frequency response characteristics of spur gear pair under internal excitation, and found that there are a lot of nonlinear phenomena such as jumping phenomenon and multi value solution in gear vibration response. Li ⁵ established the pure torsion dynamic model of spur gear pair, and studied the influence of system parameters on the nonlinear dynamic characteristics of spur gear system. Chen ⁶ established a gear-rotor-bearing coupling model and studied the circumferential and axial dynamic responses of double-helical gears with different helix angles. Liu ⁷ established the double-helical gear model based on the hybrid user-defined element method, and analyzed the vibration response of the system, which was in good agreement with the test. Yin ⁸ found that considering the time-varying oil film stiffness and time-varying damping of journal bearing can more accurately predict the dynamic characteristics of gear. Zhao ⁹ analyzed the influence of backlash, eccentricity, and tooth surface friction on the nonlinear vibration of high-speed gearbox system of wind turbine. Deng ¹⁰ studied the influence of different planet carriers and different support modes on the meshing performance of the double-helical planetary transmission system. Zhang ¹¹ established the dynamic model of double-helical gear split torque transmission system based on the finite element method, and simulated the response of the system in time domain and frequency domain. Zhang ¹² considered the phase relationship among multiple meshing of double-helical planetary gear transmission system, and studied the influence of meshing phase on the load sharing of double-helical planetary gear by integrating the support stiffness and error excitation of each component. Yuan ¹³ studied the influence of cumulative pitch error on the dynamic characteristics of double-helical gear system. Wang^14–15 studied the vibration characteristics of double-helical gear transmission system under different loads and different coincidence degrees. Mo ¹⁶ analyzed the nonlinear dynamic characteristics of double frequency of double-helical gear closed differential planetary transmission system under no-load or light load conditions. Zhou ¹⁷ established the bending-torsion-shaft model of double-helical gear, and analyzed the influence of the basic parameters such as the addendum height coefficient, the top clearance coefficient, and the tooth width on the vibration and noise of the reducer. Li ¹⁸ studied the influence of the fractal backlash on the dynamic characteristics of the gear bearing system. Zhi ¹⁹ proposed a gear nonlinear model of multi-degree of freedom considering the characteristics of gear tooth. Fractal theory was used to calculate the nonlinear backlash from the tribological aspect. Lian ²⁰ 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 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. Yi ²¹ 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. Ling ²² established a torsional nonlinear dynamic model of multi-stage planetary gear considering time-varying meshing stiffness, comprehensive gear error, and piece-wise backlash nonlinearities, and studied the effects of excitation frequency, backlash and damping on bifurcation characteristics. Siar ²³ proposed a nonlinear dynamic model of parallel shaft gear and intersecting shaft gear and solved the nonlinear algebraic Equations by multi-order harmonic balance method combined with Continuous Fourier transform. Ramadani^24–26 presents a new approach aiming to reduce gear vibration and weight by modifying its body structure. A new precise closed loop test rig was designed and produced to measure experimentally vibrations caused by rotating and lubricated gear pairs. Spitas ²⁷ analyzed the dynamic contact of gears and deduces the flexibility, contact geometry, backlash and torsional and transverse displacement and deflection in general three-dimensional multi degree of freedom systems. Sakaridis ²⁸ studied the influence of single tooth inertia on the dynamics of spur gears, and established a lumped element model, which decoupled the inertia characteristics of teeth and gears, and considered the meshing stiffness and backlash related to meshing position and load. Gkimisis ²⁹ proposed an accurate repositioning method for contact reversal of involute teeth, established a general approximate function of meshing stiffness of gear pair including load dependence, and successfully fitted it to the analysis data.

To sum up, the research on the nonlinear characteristics of double-helical gear transmission system at home and abroad mostly focuses on the establishment of bending torsion coupling model and the analysis of influence mechanism. The research on the establishment of nonlinear dynamic model considering BTAP coupling of double-helical gear system is insufficient, and the influence of time-varying meshing stiffness or backlash is often ignored. At the same time, there are few reports about the nonlinear frequency response of double-helical gear.

On this basis, the BTAP coupling dynamic model of double-helical gear pair system is established, which is more suitable for the actual working conditions, and the amplitude frequency response of double-helical gear pair system is solved. The effects of backlash, time-varying meshing stiffness, damping, comprehensive transmission error and external load excitation on frequency response characteristics are analyzed. Through the comparison with the vibration test data of double-helical gear, it provides a reference for the vibration reduction, noise reduction and structure optimization of double-helical gear.

BTAP coupling nonlinear dynamic modeling of double-helical gear pair considering backlash

Due to the deformation of elastic support and axial dynamic meshing force, there are not only torsional vibration, bending vibration and axial vibration, but also pendular vibration in double-helical gear transmission system. The BTAP coupling nonlinear dynamic modeling of double-helical gear pair is shown in in Figure 1.

Figure 1.

Bending-torsional-axial-pendular coupling nonlinear dynamic model of double-helical gear.

Here, T_pL, T_pR, T_gL, and T_gR represent the input torque and output torque acting on the driving helical gear and the driven helical gear, respectively. k_mj and c_mj represent the meshing stiffness and damping of the left and right helical gears. e _j represents the static transmission error of the left and right helical gears. The intermediate shaft section of the double-helical gear is regarded as a torsion bar unit with two nodes. k_it and c_it represent the torsional stiffness and damping of the intermediate shaft section of the double-helical gear. k _ijx , k _ijy , k _ijz , c _ijx , c _ijy , and c _ijz represent the support stiffness and support damping along the coordinate axis x, y, and z, respectively. k _ijθx , k _ijθy , c _ijθx , and c _ijθy represent the torsional support stiffness and torsional support damping along the coordinate axis x and y, respectively. Subscript i = p, g for driving and driven gears, subscript j= L, R for left and right equivalent gears.

Figure 2 shows the axial projection of double-helical gear pair.

The positive direction of rotation angular displacement of driving and driven wheels is shown in the Figure 2. The gear pressure angle is α. The included angle between center distance O_pO_g and x axis is φ. The included angle between meshing line MN and y axis is ψ.

Figure 2.

Axial projection of double-helical gear pair.

According to Figure 2, the angle between the meshing line MN and the y-axis ψ is

$$\psi =\alpha -\phi$$

(1)

If the driving gear rotates anticlockwise, the angle between the meshing line MN and the y-axis ψ is

$$\psi =-\left(\alpha +\phi \right)$$

(2)

According to Newton’s second law and the structural characteristics of double-helical gear, the BTAP coupling nonlinear dynamic differential equation is established. The differential equation of left driving gear of double-helical gear is as follows

$$\{\begin{array}{l}{m}_{\text{pL}}{\ddot{x}}_{\text{pL}}+c_{\text{pLx}}{\dot{x}}_{\text{pL}}+k_{\text{pLx}}{x}_{\text{pL}}+c_{\text{ptx}}({\dot{x}}_{\text{pL}}-{\dot{x}}_{\text{pR}})+k_{\text{ptx}}(x_{\text{pL}}-x_{\text{pR}})=-(k_{\text{mL}}{\delta }_{\text{L}}+c_{\text{mL}}{\dot{\delta }}_{\text{L}})\cos {\beta }_{\text{L}}\sin \psi \\ m_{\text{pL}}{\ddot{y}}_{\text{pL}}+c_{\text{pLy}}{\dot{y}}_{\text{pL}}+k_{\text{pLy}}{y}_{\text{pL}}+c_{\text{pty}}({\dot{y}}_{\text{pL}}-{\dot{y}}_{\text{pR}})+k_{\text{pty}}(y_{\text{pL}}-y_{\text{pR}})=-(k_{\text{mL}}{\delta }_{\text{L}}+c_{\text{mL}}{\dot{\delta }}_{\text{L}})\cos {\beta }_{\text{L}}\text{cos}\psi \\ m_{\text{pL}}{\ddot{z}}_{\text{pL}}+c_{\text{ptz}}({\dot{z}}_{\text{pL}}-{\dot{z}}_{\text{pR}})+k_{\text{ptz}}(z_{\text{pL}}-z_{\text{pR}})=(k_{\text{mL}}{\delta }_{\text{L}}+c_{\text{mL}}{\dot{\delta }}_{\text{L}})\text{sin}{\beta }_{\text{L}}\\ I_{\text{pL}}{\ddot{\theta }}_{\text{xpL}}+k_{\text{pL}\theta x}{\theta }_{\text{xpL}}+c_{\text{pL}\theta x}{\dot{\theta }}_{\text{xpL}}+(k_{\text{mL}}{\delta }_{\text{L}}+c_{\text{mL}}{\dot{\delta }}_{\text{L}})R_{\text{pL}}\sin {\beta }_{\text{L}}\sin \psi =0\\ I_{\text{pL}}{\ddot{\theta }}_{\text{ypL}}+k_{\text{pL}\theta y}{\theta }_{\text{ypL}}+c_{\text{pL}\theta y}{\dot{\theta }}_{\text{ypL}}+(k_{\text{mL}}{\delta }_{\text{L}}+c_{\text{mL}}{\dot{\delta }}_{\text{L}})R_{\text{pL}}\sin {\beta }_{\text{L}}\cos \psi =0\\ I_{\text{pL}}{\ddot{\theta }}_{\text{pL}}=T_{\text{pL}}-\cos {\beta }_{\text{L}}{R}_{\text{pL}}(k_{\text{mL}}{\delta }_{\text{L}}+c_{\text{mL}}{\dot{\delta }}_{\text{L}})-[k_{\text{pt}}({\theta }_{\text{pL}}-{\theta }_{\text{pR}})+c_{\text{pt}}({\dot{\theta }}_{\text{pL}}-{\dot{\theta }}_{\text{pR}})]\end{array}$$

(3)

The differential equation of left driven gear of double-helical gear is as follows

$$\{\begin{array}{l}{m}_{\text{gL}}{\ddot{x}}_{\text{gL}}+c_{\text{gLx}}{\dot{x}}_{\text{gL}}+k_{\text{gLx}}{x}_{\text{gL}}+c_{\text{gtx}}\left({\dot{x}}_{\text{gL}}-{\dot{x}}_{\text{gR}}\right)+k_{\text{gtx}}\left(x_{\text{gL}}-x_{\text{gR}}\right)=\left(k_{\text{mL}}{\delta }_{\text{L}}+c_{\text{mL}}{\dot{\delta }}_{\text{L}}\right)\cos \hspace{0.17em}{\beta }_{\text{L}}\sin \hspace{0.17em}\psi \\ m_{\text{gL}}{\ddot{y}}_{\text{gL}}+c_{\text{gLy}}{\dot{y}}_{\text{gL}}+k_{\text{gLy}}{y}_{\text{gL}}+c_{\text{gty}}\left({\dot{y}}_{\text{gL}}-{\dot{y}}_{\text{gR}}\right)+k_{\text{gty}}\left(y_{\text{gL}}-y_{\text{gR}}\right)=\left(k_{\text{mL}}{\delta }_{\text{L}}+c_{\text{mL}}{\dot{\delta }}_{\text{L}}\right)\cos \hspace{0.17em}{\beta }_{\text{L}}\text{cos}\hspace{0.17em}\psi \\ m_{\text{gL}}{\ddot{z}}_{\text{gL}}+c_{\text{gLz}}{\dot{z}}_{\text{gL}}+k_{\text{gLz}}{z}_{\text{gL}}+c_{\text{gtz}}\left({\dot{z}}_{\text{gL}}-{\dot{z}}_{\text{gR}}\right)+k_{\text{gtz}}\left(z_{\text{gL}}-z_{\text{gR}}\right)=-\left(k_{\text{mL}}{\delta }_{\text{L}}+c_{\text{mL}}{\dot{\delta }}_{\text{L}}\right)\text{sin}\hspace{0.17em}{\beta }_{\text{L}}\\ I_{\text{gL}}{\ddot{\theta }}_{\text{xgL}}+k_{\text{gL}\theta x}{\theta }_{\text{xgL}}+c_{\text{gL}\theta x}{\dot{\theta }}_{\text{xgL}}-\left(k_{\text{mL}}{\delta }_{\text{L}}+c_{\text{mL}}{\dot{\delta }}_{\text{L}}\right)R_{\text{gL}}\sin \hspace{0.17em}{\beta }_{\text{L}}\sin \hspace{0.17em}\psi =0\\ I_{\text{gL}}{\ddot{\theta }}_{\text{ygL}}+k_{\text{gL}\theta y}{\theta }_{\text{ygL}}+c_{\text{gL}\theta y}{\dot{\theta }}_{\text{ygL}}-\left(k_{\text{mL}}{\delta }_{\text{L}}+c_{\text{mL}}{\dot{\delta }}_{\text{L}}\right)R_{\text{gL}}\sin \hspace{0.17em}{\beta }_{\text{L}}\cos \hspace{0.17em}\psi =0\\ I_{\text{gL}}{\ddot{\theta }}_{\text{gL}}=-T_{\text{gL}}+\cos \hspace{0.17em}{\beta }_{\text{L}}{R}_{\text{gL}}\left(k_{\text{mL}}{\delta }_{\text{L}}+c_{\text{mL}}{\dot{\delta }}_{\text{L}}\right)-\left[k_{\text{gt}}\left({\theta }_{\text{gL}}-{\theta }_{\text{gR}}\right)+c_{\text{gt}}\left({\dot{\theta }}_{\text{gL}}-{\dot{\theta }}_{\text{gR}}\right)\right]\end{array}$$

(4)

The differential equation of right driving gear of double-helical gear is as follows

$$\{\begin{array}{l}{m}_{\text{pR}}{\ddot{x}}_{\text{pR}}+c_{\text{pRx}}{\dot{x}}_{\text{pR}}+k_{\text{pRx}}{x}_{\text{pR}}+c_{\text{ptx}}\left({\dot{x}}_{\text{pR}}-{\dot{x}}_{\text{pL}}\right)+k_{\text{ptx}}\left(x_{\text{pR}}-x_{\text{pL}}\right)=-\left(k_{\text{mR}}{\delta }_{\text{R}}+c_{\text{mR}}{\dot{\delta }}_{\text{R}}\right)\cos \hspace{0.17em}{\beta }_{\text{R}}\sin \hspace{0.17em}\psi \\ m_{\text{pR}}{\ddot{y}}_{\text{pR}}+c_{\text{pRy}}{\dot{y}}_{\text{pR}}+k_{\text{pRy}}{y}_{\text{pR}}+c_{\text{pty}}\left({\dot{y}}_{\text{pR}}-{\dot{y}}_{\text{pL}}\right)+k_{\text{pty}}\left(y_{\text{pR}}-y_{\text{pL}}\right)=-\left(k_{\text{mR}}{\delta }_{\text{R}}+c_{\text{mR}}{\dot{\delta }}_{\text{R}}\right)\cos \hspace{0.17em}{\beta }_{\text{R}}\text{cos}\hspace{0.17em}\psi \\ m_{\text{pR}}{\ddot{z}}_{\text{pR}}+c_{\text{ptz}}\left({\dot{z}}_{\text{pR}}-{\dot{z}}_{\text{pL}}\right)+k_{\text{ptz}}\left(z_{\text{pR}}-z_{\text{pL}}\right)=\left(k_{\text{mR}}{\delta }_{\text{R}}+c_{\text{mR}}{\dot{\delta }}_{\text{R}}\right)\text{sin}\hspace{0.17em}{\beta }_{\text{R}}\\ I_{\text{pR}}{\ddot{\theta }}_{\text{xpR}}+k_{\text{pR}\theta x}{\theta }_{\text{xpR}}+c_{\text{pR}\theta x}{\dot{\theta }}_{\text{xpR}}+\left(k_{\text{mR}}{\delta }_{\text{R}}+c_{\text{mR}}{\dot{\delta }}_{\text{R}}\right)R_{\text{pR}}\sin \hspace{0.17em}{\beta }_{\text{R}}\sin \hspace{0.17em}\psi =0\\ I_{\text{pR}}{\ddot{\theta }}_{\text{ypR}}+k_{\text{pR}\theta y}{\theta }_{\text{ypR}}+c_{\text{pR}\theta y}{\dot{\theta }}_{\text{ypR}}+\left(k_{\text{mR}}{\delta }_{\text{R}}+c_{\text{mR}}{\dot{\delta }}_{\text{R}}\right)R_{\text{pR}}\sin \hspace{0.17em}{\beta }_{\text{R}}\cos \hspace{0.17em}\psi =0\\ I_{\text{pR}}{\ddot{\theta }}_{\text{pR}}=T_{\text{pR}}-\cos \hspace{0.17em}{\beta }_{\text{R}}{R}_{\text{pR}}\left(k_{\text{mR}}{\delta }_{\text{R}}+c_{\text{mR}}{\dot{\delta }}_{\text{R}}\right)-\left[k_{\text{pt}}\left({\theta }_{\text{pR}}-{\theta }_{\text{pL}}\right)+c_{\text{pt}}\left({\dot{\theta }}_{\text{pR}}-{\dot{\theta }}_{\text{pL}}\right)\right]\end{array}$$

(5)

The differential equation of right driven gear of double-helical gear is as follows

$$\{\begin{array}{l}{m}_{\text{gR}}{\ddot{x}}_{\text{gR}}+c_{\text{gRx}}{\dot{x}}_{\text{gR}}+k_{\text{gRx}}{x}_{\text{gR}}+c_{\text{gtx}}\left({\dot{x}}_{\text{gR}}-{\dot{x}}_{\text{gL}}\right)+k_{\text{gtx}}\left(x_{\text{gR}}-x_{\text{gL}}\right)=\left(k_{\text{mR}}{\delta }_{\text{R}}+c_{\text{mR}}{\dot{\delta }}_{\text{R}}\right)\cos \hspace{0.17em}{\beta }_{\text{R}}\sin \hspace{0.17em}\psi \\ m_{\text{gR}}{\ddot{y}}_{\text{gR}}+c_{\text{gRy}}{\dot{y}}_{\text{gR}}+k_{\text{gRy}}{y}_{\text{gR}}+c_{\text{gty}}\left({\dot{y}}_{\text{gR}}-{\dot{y}}_{\text{gL}}\right)+k_{\text{gty}}\left(y_{\text{gR}}-y_{\text{gL}}\right)=\left(k_{\text{mR}}{\delta }_{\text{R}}+c_{\text{mR}}{\dot{\delta }}_{\text{R}}\right)\cos \hspace{0.17em}{\beta }_{\text{R}}\text{cos}\hspace{0.17em}\psi \\ m_{\text{gR}}{\ddot{z}}_{\text{gR}}+c_{\text{gRz}}{\dot{z}}_{\text{gR}}+k_{\text{gRz}}{z}_{\text{gR}}+c_{\text{gtz}}\left({\dot{z}}_{\text{gR}}-{\dot{z}}_{\text{gL}}\right)+k_{\text{gtz}}\left(z_{\text{gR}}-z_{\text{gL}}\right)=-\left(k_{\text{mR}}{\delta }_{\text{R}}+c_{\text{mR}}{\dot{\delta }}_{\text{R}}\right)\text{sin}\hspace{0.17em}{\beta }_{\text{R}}\\ I_{\text{gR}}{\ddot{\theta }}_{\text{xgR}}+k_{\text{gR}\theta x}{\theta }_{\text{xgR}}+c_{\text{gR}\theta x}{\dot{\theta }}_{\text{xgR}}-\left(k_{\text{mR}}{\delta }_{\text{R}}+c_{\text{mR}}{\dot{\delta }}_{\text{R}}\right)R_{\text{gR}}\sin \hspace{0.17em}{\beta }_{\text{R}}\sin \hspace{0.17em}\psi =0\\ I_{\text{gR}}{\ddot{\theta }}_{\text{ygR}}+k_{\text{gR}\theta y}{\theta }_{\text{ygR}}+c_{\text{gR}\theta y}{\dot{\theta }}_{\text{ygR}}-\left(k_{\text{mR}}{\delta }_{\text{R}}+c_{\text{mR}}{\dot{\delta }}_{\text{R}}\right)R_{\text{gR}}\sin \hspace{0.17em}{\beta }_{\text{R}}\cos \hspace{0.17em}\psi =0\\ I_{\text{gR}}{\ddot{\theta }}_{\text{gR}}=-T_{\text{gR}}+\cos \hspace{0.17em}{\beta }_{\text{R}}{R}_{\text{gR}}\left(k_{\text{mR}}{\delta }_{\text{R}}+c_{\text{mR}}{\dot{\delta }}_{\text{R}}\right)-\left[k_{\text{gt}}\left({\theta }_{\text{gR}}-{\theta }_{\text{gL}}\right)+c_{\text{gt}}\left({\dot{\theta }}_{\text{gR}}-{\dot{\theta }}_{\text{gL}}\right)\right]\end{array}$$

(6)

Here, I _ij represents the polar moment of inertia on the equivalent nodes O_pL, O_pR, O_gL, and O_gR at the center of each helical gear. θ _ij represents the rotational angular displacement of the equivalent nodes at the center of each helical gear. m _ij represents the equivalent node mass at the center of each helical gear. r _ij represents the base circle radius of each helical gear. β _j represents the helical angle of the left and right helical gear pairs. δ _j represents the relative meshing displacement of the left and right helical gear pairs. x _ij , y _ij , z _ij , θ _xij , θ _yij , and θ _ij respectively represent the translational vibration displacement of each equivalent node along the x, y, z directions of the coordinate axes and the rotational angular displacement around the x, y, z directions of the coordinate axes, respectively.

The time-varying meshing stiffness function and the comprehensive transmission error function are expressed in the form of Fourier series

$$k_{\text{m}j}\left(t\right)=k_{\text{mean}}\left[1-{\displaystyle \sum _{l=1}^L{\epsilon }_l\cos \left(l\omega t\right)}\right]$$

(7)

$$e_j\left(t\right)=e_{\text{m}}+{\displaystyle \sum _{l=1}^L{e}_l\sin \left(l\omega t\right)}$$

(8)

Here, kmean is the average stiffness. ɛl is the fluctuation coefficient reflecting the change of stiffness. em is the average error. el is the fluctuation coefficient reflecting the change of comprehensive transmission error. L is the expansion order.

m_c is the equivalent mass of left and right helical gear pair, which is expressed as follows

$$m_{\text{c}}=\frac{I_{\text{pL}}{I}_{\text{gL}}}{I_{\text{pL}}{r}_{\text{gL}}^2+I_{\text{gL}}{r}_{\text{pL}}^2}=\frac{I_{\text{pR}}{I}_{\text{gR}}}{I_{\text{pR}}{r}_{\text{gR}}^2+I_{\text{gR}}{r}_{\text{pR}}^2}$$

(9)

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

$$c_{\text{m}j}=2\xi \sqrt{m_{\text{c}}{k}_{\text{mean}}}\hspace{1em}j=\text{L},\hspace{0.17em}\text{R}$$

(10)

$$c_{i\text{t}}=2\xi \sqrt{m_{\text{c}}{k}_{i\text{t}}}\hspace{1em}i=\text{p},\hspace{0.17em}\text{g}$$

(11)

Where, ζ is the damping ratio of the system.

According to the knowledge of material mechanics, the calculation formulas of bending stiffness, tensile stiffness and torsional stiffness of intermediate shaft section are as follows

$$\{\begin{array}{l}{k}_{i\text{tx}}=\frac{12E{I}_x}{l^3};\hspace{1em}{k}_{i\text{ty}}=\frac{12E{I}_y}{l^3}\\ k_{i\text{tz}}=\frac{EA}{l};\hspace{1em}{k}_{i\text{t}}=\frac{GJ}{l}\end{array}\hspace{1em}\left(i=\text{p},\text{g}\right)$$

(12)

Here, E and G are young’s modulus of elasticity and shear modulus of gear material, respectively. I _x , I _y and J are moment of inertia and polar moment of inertia of cross section about x and y axis, respectively, which can be calculated by formula (13). l and A are length and cross-sectional area of undercut respectively. D and d are outer diameter and inner diameter of shaft, respectively

$$\{\begin{array}{l}A=\frac{\pi \left(D^2-d^2\right)}{4}\\ I_x=I_y=\frac{\pi \left(D^4-d^4\right)}{64}\\ J=\frac{\pi \left(D^4-d^4\right)}{32}\end{array}$$

(13)

The dynamic Equations shown in formulas (3) to (6) are expressed in matrix form as follows

$$M\ddot{q}+C\dot{q}+Kq=F$$

(14)

M, C, K and F are the mass matrix, damping matrix, stiffness matrix, and torque vector of the system, respectively. The generalized displacement array of the system with 24 degrees of freedom is expressed as follows

$$q=\left\{\begin{array}{l}{x}_{\text{pL}},y_{\text{pL}},z_{\text{pL}},{\theta }_{\text{xpL}},{\theta }_{\text{ypL}},{\theta }_{\text{pL}},x_{\text{gL}},y_{\text{gL}},z_{\text{gL}},{\theta }_{\text{xgL}},{\theta }_{\text{ygL}},{\theta }_{\text{gL}},x_{\text{pR}},y_{\text{pR}},z_{\text{pR}},{\theta }_{\text{xpR}},\\ {\theta }_{\text{ypR}},{\theta }_{\text{pR}},x_{\text{gR}},y_{\text{gR}},z_{\text{gR}},{\theta }_{\text{xgR}},{\theta }_{\text{ygR}},{\theta }_{\text{gR}}\end{array}\right\}$$

(15)

Before solving the BTAP coupling dynamic Equation of double-helical gear, because of the torsional displacement of each concentrated mass in the system is not constrained, it is necessary to use the relative displacement method to eliminate the rigid body displacement of the system to solve the differential equation. In addition to the relative meshing displacements δ_L and δ_R of the left and right helical gears along the meshing line, the relative torsional linear displacements γ_P and γ_G between the adjacent helical gears in the same axis should be introduced.

The expression of relative meshing displacement in the direction of meshing line of helical gear pair at left and right ends of double-helical gear is as follows

$$\{\begin{array}{l}{\delta }_{\text{L}}=\left[\left(x_{\text{pL}}-x_{\text{gL}}\right)\sin \hspace{0.17em}\psi +\left(y_{\text{pL}}-y_{\text{gL}}\right)\cos \hspace{0.17em}\psi +r_{\text{pL}}{\theta }_{\text{pL}}-r_{\text{gL}}{\theta }_{\text{gL}}\right]\cos \hspace{0.17em}{\beta }_{\text{L}}+\\ \left[\left(r_{\text{pL}}{\theta }_{\text{xpL}}-r_{\text{gL}}{\theta }_{\text{xgL}}\right)\sin \hspace{0.17em}\psi +\left(r_{\text{pL}}{\theta }_{\text{ypL}}-r_{\text{gL}}{\theta }_{\text{ygL}}\right)\cos \hspace{0.17em}\psi +\left(z_{\text{pL}}-z_{\text{gL}}\right)\right]\sin \hspace{0.17em}{\beta }_{\text{L}}-e_{\text{L}}\left(t\right)\\ {\delta }_{\text{R}}=\left[\left(x_{\text{pR}}-x_{\text{gR}}\right)\sin \hspace{0.17em}\psi +\left(y_{\text{pR}}-y_{\text{gR}}\right)\cos \hspace{0.17em}\psi +r_{\text{pR}}{\theta }_{\text{pR}}-r_{\text{gR}}{\theta }_{\text{gR}}\right]\cos \hspace{0.17em}{\beta }_{\text{R}}+\\ \left[\left(r_{\text{pR}}{\theta }_{\text{xpR}}-r_{\text{gR}}{\theta }_{\text{xgR}}\right)\sin \hspace{0.17em}\psi +\left(r_{\text{pR}}{\theta }_{\text{ypR}}-r_{\text{gR}}{\theta }_{\text{ygR}}\right)\cos \hspace{0.17em}\psi +\left(z_{\text{pR}}-z_{\text{gR}}\right)\right]\sin \hspace{0.17em}{\beta }_{\text{R}}-e_{\text{R}}\left(t\right)\end{array}$$

(16)

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

$$\{\begin{array}{l}{\gamma }_{\text{P}}=\left({\theta }_{\text{pL}}-{\theta }_{\text{pR}}\right)r_{\text{pL}}\\ {\gamma }_{\text{G}}=\left({\theta }_{\text{gL}}-{\theta }_{\text{gR}}\right)r_{\text{gL}}\end{array}$$

(17)

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

$$f\left({\delta }_j\right)=\{\begin{array}{l}{\delta }_j-b\hspace{1em}{\delta }_j>b\\ 0-b\hspace{1em}\le {\delta }_j\le b\\ {\delta }_j+b\hspace{1em}{\delta }_j<-b\end{array}$$

(18)

Here, 2b is the length of backlash, and the expression is as follows

$$b=b_{\text{0}}\left(1+W\right)$$

(19)

Here, b₀ is the average backlash. W is the change of dynamic backlash to describe tooth wear.

Other dimensionless parameters can be expressed as follows

$$\{\begin{array}{l}{\overline{x}}_{ij}=\frac{x_{ij}}{b_{\text{c}}},\hspace{1em}{\dot{\overline{x}}}_{ij}=\frac{{\dot{x}}_{ij}}{b_{\text{c}}{\omega }_n},\hspace{1em}{\ddot{\overline{x}}}_{ij}=\frac{{\ddot{x}}_{ij}}{b_{\text{c}}{\omega }_n^2}\\ {\Delta }_j=\frac{{\delta }_j}{b_c},\hspace{1em}{\Delta }_i=\frac{{\gamma }_i}{b_c}\\ {\omega }_n=\sqrt{k_{\text{mean}}/m_{\text{c}}},\hspace{1em}\Omega =\frac{\omega }{{\omega }_n}\\ \overline{t}={\omega }_{n}t,\hspace{1em}\xi =\frac{c}{2m_{\text{c}}{\omega }_n}\end{array}$$

(20)

where b _c is the nominal displacement scale. ${\overline{x}}_{ij}$ , ${\Delta }_i$ , ${\Delta }_j$ are the normalized vibration displacement. ${\omega }_n$ is the natural frequency. $\overline{t}$ is the normalized time. Ω is the normalized meshing frequency. The displacement $y_{ij}$ , $z_{ij}$ can also be replaced by dimensionless ${\overline{y}}_{ij}$ and ${\overline{z}}_{ij}$ in the dynamic Equation.

By substituting the above Equations into Equations (3) to (6), the system dynamic differential equations for eliminating rigid body displacement and dimensionalization are obtained as follows. Here, the differential equation of left driving gear of double-helical gear is as follows

$$\{\begin{array}{l}{\ddot{\overline{x}}}_{\text{pL}}+\frac{c_{\text{pLx}}{\dot{\overline{x}}}_{\text{pL}}}{m_{\text{pL}}{\omega }_n}+\frac{k_{\text{pLx}}{\overline{x}}_{\text{pL}}}{m_{\text{pL}}{\omega }_n^2}+\frac{c_{\text{ptx}}\left({\dot{\overline{x}}}_{\text{pL}}-{\dot{\overline{x}}}_{\text{pR}}\right)}{m_{\text{pL}}{\omega }_n}+\frac{k_{\text{ptx}}\left({\overline{x}}_{\text{pL}}-{\overline{x}}_{\text{pR}}\right)}{m_{\text{pL}}{\omega }_n^2}=\frac{-k_{\text{mL}}{\mathrm{\Delta }}_{\text{L}}\cos \hspace{0.17em}{\beta }_{\text{L}}\sin \hspace{0.17em}\psi }{m_{\text{pL}}{\omega }_n^2}-\frac{c_{\text{mL}}{\dot{\mathrm{\Delta }}}_{\text{L}}\cos \hspace{0.17em}{\beta }_{\text{L}}\sin \hspace{0.17em}\psi }{m_{\text{pL}}{\omega }_n}\\ {\ddot{\overline{y}}}_{\text{pL}}+\frac{c_{\text{pLy}}{\dot{\overline{y}}}_{\text{pL}}}{m_{\text{pL}}{\omega }_n}+\frac{k_{\text{pLy}}{\overline{y}}_{\text{pL}}}{m_{\text{pL}}{\omega }_n^2}+\frac{c_{\text{pty}}\left({\dot{\overline{y}}}_{\text{pL}}-{\dot{\overline{y}}}_{\text{pR}}\right)}{m_{\text{pL}}{\omega }_n}+\frac{k_{\text{pty}}\left({\overline{y}}_{\text{pL}}-{\overline{y}}_{\text{pR}}\right)}{m_{\text{pL}}{\omega }_n^2}=\frac{-k_{\text{mL}}{\mathrm{\Delta }}_{\text{L}}\cos \hspace{0.17em}{\beta }_{\text{L}}\text{cos}\hspace{0.17em}\psi }{m_{\text{pL}}{\omega }_n^2}-\frac{c_{\text{mL}}{\dot{\mathrm{\Delta }}}_{\text{L}}\cos \hspace{0.17em}{\beta }_{\text{L}}\text{cos}\hspace{0.17em}\psi }{m_{\text{pL}}{\omega }_n}\\ {\ddot{\overline{z}}}_{\text{pL}}+\frac{c_{\text{ptz}}\left({\dot{\overline{z}}}_{\text{pL}}-{\dot{\overline{z}}}_{\text{pR}}\right)}{m_{\text{pL}}{\omega }_n}+\frac{k_{\text{ptz}}\left({\overline{z}}_{\text{pL}}-{\overline{z}}_{\text{pR}}\right)}{m_{\text{pL}}{\omega }_n^2}=\frac{k_{\text{mL}}{\mathrm{\Delta }}_{\text{L}}\text{sin}\hspace{0.17em}{\beta }_{\text{L}}}{m_{\text{pL}}{\omega }_n^2}+\frac{c_{\text{mL}}{\dot{\mathrm{\Delta }}}_{\text{L}}\text{sin}\hspace{0.17em}{\beta }_{\text{L}}}{m_{\text{pL}}{\omega }_n}\end{array}$$

(21)

The differential equation of left driven gear of double-helical gear is as follow

$$\{\begin{array}{l}{\ddot{\overline{x}}}_{\text{gL}}+\frac{c_{\text{gLx}}{\dot{\overline{x}}}_{\text{gL}}}{m_{\text{gL}}{\omega }_n}+\frac{k_{\text{gLx}}{\overline{x}}_{\text{gL}}}{m_{\text{gL}}{\omega }_n^2}+\frac{c_{\text{gtx}}\left({\dot{\overline{x}}}_{\text{gL}}-{\dot{\overline{x}}}_{\text{gR}}\right)}{m_{\text{gL}}{\omega }_n}+\frac{k_{\text{gtx}}\left({\overline{x}}_{\text{gL}}-{\overline{x}}_{\text{gR}}\right)}{m_{\text{gL}}{\omega }_n^2}=\frac{k_{\text{mL}}{\mathrm{\Delta }}_{\text{L}}\cos \hspace{0.17em}{\beta }_{\text{L}}\sin \hspace{0.17em}\psi }{m_{\text{gL}}{\omega }_n^2}+\frac{c_{\text{mL}}{\dot{\mathrm{\Delta }}}_{\text{L}}\cos \hspace{0.17em}{\beta }_{\text{L}}\sin \hspace{0.17em}\psi }{m_{\text{gL}}{\omega }_n}\\ {\ddot{\overline{y}}}_{\text{gL}}+\frac{c_{\text{gLy}}{\dot{\overline{y}}}_{\text{gL}}}{m_{\text{gL}}{\omega }_n}+\frac{k_{\text{gLy}}{\overline{y}}_{\text{gL}}}{m_{\text{gL}}{\omega }_n^2}+\frac{c_{\text{gty}}\left({\dot{\overline{y}}}_{\text{gL}}-{\dot{\overline{y}}}_{\text{gR}}\right)}{m_{\text{gL}}{\omega }_n}+\frac{k_{\text{gty}}\left({\overline{y}}_{\text{gL}}-{\overline{y}}_{\text{gR}}\right)}{m_{\text{gL}}{\omega }_n^2}=\frac{k_{\text{mL}}{\mathrm{\Delta }}_{\text{L}}\cos \hspace{0.17em}{\beta }_{\text{L}}\text{cos}\hspace{0.17em}\psi }{m_{\text{gL}}{\omega }_n^2}+\frac{c_{\text{mL}}{\dot{\mathrm{\Delta }}}_{\text{L}}\cos \hspace{0.17em}{\beta }_{\text{L}}\text{cos}\hspace{0.17em}\psi }{m_{\text{gL}}{\omega }_n}\\ {\ddot{\overline{z}}}_{\text{gL}}+\frac{c_{\text{gLz}}{\dot{\overline{z}}}_{\text{gL}}}{m_{\text{gL}}{\omega }_n}+\frac{k_{\text{gLz}}{\overline{z}}_{\text{gL}}}{m_{\text{gL}}{\omega }_n^2}+\frac{c_{\text{gtz}}\left({\dot{\overline{z}}}_{\text{gL}}-{\dot{\overline{z}}}_{\text{gR}}\right)}{m_{\text{gL}}{\omega }_n}+\frac{k_{\text{gtz}}\left({\overline{z}}_{\text{gL}}-{\overline{z}}_{\text{gR}}\right)}{m_{\text{gL}}{\omega }_n^2}=\frac{-k_{\text{mL}}{\mathrm{\Delta }}_{\text{L}}\text{sin}\hspace{0.17em}{\beta }_{\text{L}}}{m_{\text{gL}}{\omega }_n^2}-\frac{c_{\text{mL}}{\dot{\mathrm{\Delta }}}_{\text{L}}\text{sin}\hspace{0.17em}{\beta }_{\text{L}}}{m_{\text{gL}}{\omega }_n}\end{array}$$

(22)

The differential equation of right driving gear of double-helical gear is as follows

$$\{\begin{array}{l}{\ddot{\overline{x}}}_{\text{pR}}+\frac{c_{\text{pRx}}{\dot{\overline{x}}}_{\text{pR}}}{m_{\text{pR}}{\omega }_n}+\frac{k_{\text{pRx}}{\overline{x}}_{\text{pR}}}{m_{\text{pR}}{\omega }_n^2}+\frac{c_{\text{ptx}}\left({\dot{\overline{x}}}_{\text{pR}}-{\dot{\overline{x}}}_{\text{pL}}\right)}{m_{\text{pR}}{\omega }_n}+\frac{k_{\text{ptx}}\left({\overline{x}}_{\text{pR}}-{\overline{x}}_{\text{pL}}\right)}{m_{\text{pR}}{\omega }_n^2}=\frac{-k_{\text{mR}}{\mathrm{\Delta }}_{\text{R}}\cos \hspace{0.17em}{\beta }_{\text{R}}\sin \hspace{0.17em}\psi }{m_{\text{pR}}{\omega }_n^2}-\frac{c_{\text{mR}}{\dot{\mathrm{\Delta }}}_{\text{R}}\cos \hspace{0.17em}{\beta }_{\text{R}}\sin \hspace{0.17em}\psi }{m_{\text{pR}}{\omega }_n}\\ {\ddot{\overline{y}}}_{\text{pR}}+\frac{c_{\text{pRy}}{\dot{\overline{y}}}_{\text{pR}}}{m_{\text{pR}}{\omega }_n}+\frac{k_{\text{pRy}}{\overline{y}}_{\text{pR}}}{m_{\text{pR}}{\omega }_n^2}+\frac{c_{\text{pty}}\left({\dot{\overline{y}}}_{\text{pR}}-{\dot{\overline{y}}}_{\text{pL}}\right)}{m_{\text{pR}}{\omega }_n}+\frac{k_{\text{pty}}\left({\overline{y}}_{\text{pR}}-{\overline{y}}_{\text{pL}}\right)}{m_{\text{pR}}{\omega }_n^2}=\frac{-k_{\text{mR}}{\mathrm{\Delta }}_{\text{R}}\cos \hspace{0.17em}{\beta }_{\text{R}}\text{cos}\hspace{0.17em}\psi }{m_{\text{pR}}{\omega }_n^2}-\frac{c_{\text{mR}}{\dot{\mathrm{\Delta }}}_{\text{R}}\cos \hspace{0.17em}{\beta }_{\text{R}}\text{cos}\hspace{0.17em}\psi }{m_{\text{pR}}{\omega }_n}\\ {\ddot{\overline{z}}}_{\text{pR}}+\frac{c_{\text{ptz}}\left({\dot{\overline{z}}}_{\text{pR}}-{\dot{\overline{z}}}_{\text{pL}}\right)}{m_{\text{pR}}{\omega }_n}+\frac{k_{\text{ptz}}\left({\overline{z}}_{\text{pR}}-{\overline{z}}_{\text{pL}}\right)}{m_{\text{pR}}{\omega }_n^2}=\frac{k_{\text{mR}}{\mathrm{\Delta }}_{\text{R}}\text{sin}\hspace{0.17em}{\beta }_{\text{R}}}{m_{\text{pR}}{\omega }_n^2}+\frac{c_{\text{mR}}{\dot{\mathrm{\Delta }}}_{\text{R}}\text{sin}\hspace{0.17em}{\beta }_{\text{R}}}{m_{\text{pR}}{\omega }_n}\end{array}$$

(23)

The differential equation of right driven gear of double-helical gear is as follows

$$\{\begin{array}{l}{\ddot{\overline{x}}}_{\text{gR}}+\frac{c_{\text{gRx}}{\dot{\overline{x}}}_{\text{gR}}}{m_{\text{gR}}{\omega }_n}+\frac{k_{\text{gRx}}{\overline{x}}_{\text{gR}}}{m_{\text{gR}}{\omega }_n^2}+\frac{c_{\text{gtx}}\left({\dot{\overline{x}}}_{\text{gR}}-{\dot{\overline{x}}}_{\text{gL}}\right)}{m_{\text{gR}}{\omega }_n}+\frac{k_{\text{gtx}}\left({\overline{x}}_{\text{gR}}-{\overline{x}}_{\text{gL}}\right)}{m_{\text{gR}}{\omega }_n^2}=\frac{k_{\text{mR}}{\mathrm{\Delta }}_{\text{R}}\cos \hspace{0.17em}{\beta }_{\text{R}}\sin \hspace{0.17em}\psi }{m_{\text{gR}}{\omega }_n^2}+\frac{c_{\text{mR}}{\dot{\mathrm{\Delta }}}_{\text{R}}\cos \hspace{0.17em}{\beta }_{\text{R}}\sin \hspace{0.17em}\psi }{m_{\text{gR}}{\omega }_n}\\ {\ddot{\overline{y}}}_{\text{gR}}+\frac{c_{\text{gRy}}{\dot{\overline{y}}}_{\text{gR}}}{m_{\text{gR}}{\omega }_n}+\frac{k_{\text{gRy}}{\overline{y}}_{\text{gR}}}{m_{\text{gR}}{\omega }_n^2}+\frac{c_{\text{gty}}\left({\dot{\overline{y}}}_{\text{gR}}-{\dot{\overline{y}}}_{\text{gL}}\right)}{m_{\text{gR}}{\omega }_n}+\frac{k_{\text{gty}}\left({\overline{y}}_{\text{gR}}-{\overline{y}}_{\text{gL}}\right)}{m_{\text{gR}}{\omega }_n^2}=\frac{k_{\text{mR}}{\mathrm{\Delta }}_{\text{R}}\cos \hspace{0.17em}{\beta }_{\text{R}}\text{cos}\hspace{0.17em}\psi }{m_{\text{gR}}{\omega }_n^2}+\frac{c_{\text{mR}}{\dot{\mathrm{\Delta }}}_{\text{R}}\cos \hspace{0.17em}{\beta }_{\text{R}}\text{cos}\hspace{0.17em}\psi }{m_{\text{gR}}{\omega }_n}\\ {\ddot{\overline{z}}}_{\text{gR}}+\frac{c_{\text{gRz}}{\dot{\overline{z}}}_{\text{gR}}}{m_{\text{gR}}{\omega }_n}+\frac{k_{\text{gRz}}{\overline{z}}_{\text{gR}}}{m_{\text{gR}}{\omega }_n^2}+\frac{c_{\text{gtz}}\left({\dot{\overline{z}}}_{\text{gR}}-{\dot{\overline{z}}}_{\text{gL}}\right)}{m_{\text{gR}}{\omega }_n}+\frac{k_{\text{gtz}}\left({\overline{z}}_{\text{gR}}-{\overline{z}}_{\text{gL}}\right)}{m_{\text{gR}}{\omega }_n^2}=\frac{-k_{\text{mR}}{\mathrm{\Delta }}_{\text{R}}\text{sin}\hspace{0.17em}{\beta }_{\text{R}}}{m_{\text{gR}}{\omega }_n^2}-\frac{c_{\text{mR}}{\dot{\mathrm{\Delta }}}_{\text{R}}\text{sin}\hspace{0.17em}{\beta }_{\text{R}}}{m_{\text{gR}}{\omega }_n}\end{array}$$

(24)

The differential equation on the meshing line of helical gear pair at the left end of double-helical gear is as follows

$$\{\begin{array}{l}{\ddot{\Delta }}_{\text{L}}=\left[\left({\ddot{\overline{x}}}_{\text{pL}}-{\ddot{\overline{x}}}_{\text{gL}}\right)\sin \hspace{0.17em}\psi +\left({\ddot{\overline{y}}}_{\text{pL}}-{\ddot{\overline{y}}}_{\text{gL}}\right)\cos \hspace{0.17em}\psi \right]\cos \hspace{0.17em}{\beta }_{\text{L}}+\left[\left({\ddot{\overline{z}}}_{\text{pL}}-{\ddot{\overline{z}}}_{\text{gL}}\right)\right]\sin \hspace{0.17em}{\beta }_{\text{L}}-\frac{c_{\text{mL}}{c}^2{\beta }_{\text{L}}{\dot{\mathrm{\Delta }}}_{\text{L}}}{m_{\text{c}}{\omega }_n}-\frac{k_{\text{mL}}{c}^2{\beta }_{\text{L}}f\left({\mathrm{\Delta }}_{\text{L}}\right)}{m_{\text{c}}{\omega }_n^2}\\ -\frac{c_{\text{pt}}{r}_{\text{pL}}^2\cos \hspace{0.17em}{\beta }_{\text{L}}}{I_{\text{pL}}{\omega }_n}{\dot{\mathrm{\Delta }}}_{\text{P}}-\frac{k_{\text{pt}}{r}_{\text{pL}}^2\cos \hspace{0.17em}{\beta }_{\text{L}}}{I_{\text{pL}}{\omega }_n^2}{\mathrm{\Delta }}_{\text{P}}+\frac{c_{\text{gt}}{r}_{\text{gL}}^2\cos \hspace{0.17em}{\beta }_{\text{L}}}{I_{\text{gL}}{\omega }_n}{\dot{\mathrm{\Delta }}}_{\text{G}}+\frac{k_{\text{gt}}{r}_{\text{gL}}^2\cos \hspace{0.17em}{\beta }_{\text{L}}}{I_{\text{gL}}{\omega }_n^2}{\mathrm{\Delta }}_{\text{G}}-\frac{c_{\text{mL}}{r}_{\text{pL}}^2{s}^2{\beta }_{\text{L}}{s}^2\psi {\dot{\mathrm{\Delta }}}_{\text{L}}}{I_{\text{pL}}{\omega }_n}\\ -\frac{k_{\text{mL}}{r}_{\text{pL}}^2{s}^2{\beta }_{\text{L}}{s}^2\psi f\left({\mathrm{\Delta }}_{\text{L}}\right)}{I_{\text{pL}}{\omega }_n^2}-\frac{c_{\text{pL}\theta x}{r}_{\text{pL}}\sin \hspace{0.17em}{\beta }_{\text{L}}\sin \hspace{0.17em}\psi {\dot{\theta }}_{\text{xpL}}}{I_{\text{pL}}{\omega }_n}-\frac{k_{\text{pL}\theta x}{r}_{\text{pL}}\sin \hspace{0.17em}{\beta }_{\text{L}}\sin \hspace{0.17em}\psi {\theta }_{\text{xpL}}}{I_{\text{pL}}{\omega }_n^2}-\frac{c_{\text{mL}}{r}_{\text{gL}}^2{s}^2{\beta }_{\text{L}}{s}^2\psi {\dot{\mathrm{\Delta }}}_{\text{L}}}{I_{\text{gL}}{\omega }_n}\\ -\frac{k_{\text{mL}}{r}_{\text{gL}}^2{s}^2{\beta }_{\text{L}}{s}^2\psi f\left({\mathrm{\Delta }}_{\text{L}}\right)}{I_{\text{gL}}{\omega }_n^2}+\frac{c_{\text{gL}\theta x}{r}_{\text{gL}}\sin \hspace{0.17em}{\beta }_{\text{L}}\sin \hspace{0.17em}\psi {\dot{\theta }}_{\text{xgL}}}{I_{\text{gL}}{\omega }_n}+\frac{k_{\text{gL}\theta x}{r}_{\text{gL}}\sin \hspace{0.17em}{\beta }_{\text{L}}\sin \hspace{0.17em}\psi {\theta }_{\text{xgL}}}{I_{\text{gL}}{\omega }_n^2}-\frac{c_{\text{mL}}{r}_{\text{pL}}^2{s}^2{\beta }_{\text{L}}{c}^2\psi {\dot{\mathrm{\Delta }}}_{\text{L}}}{I_{\text{pL}}{\omega }_n}\\ -\frac{k_{\text{mL}}{r}_{\text{pL}}^2{s}^2{\beta }_{\text{L}}{c}^2\psi f\left({\mathrm{\Delta }}_{\text{L}}\right)}{I_{\text{pL}}{\omega }_n^2}-\frac{c_{\text{pL}\theta y}{r}_{\text{pL}}\sin \hspace{0.17em}{\beta }_{\text{L}}\cos \hspace{0.17em}\psi {\dot{\theta }}_{\text{ypL}}}{I_{\text{pL}}{\omega }_n}-\frac{k_{\text{pL}\theta y}{r}_{\text{pL}}\sin \hspace{0.17em}{\beta }_{\text{L}}\cos \hspace{0.17em}\psi {\theta }_{\text{ypL}}}{I_{\text{pL}}{\omega }_n^2}-\frac{c_{\text{mL}}{r}_{\text{gL}}^2{s}^2{\beta }_{\text{L}}{c}^2\psi {\dot{\mathrm{\Delta }}}_{\text{L}}}{I_{\text{gL}}{\omega }_n}\\ -\frac{k_{\text{mL}}{r}_{\text{gL}}^2{s}^2{\beta }_{\text{L}}{c}^2\psi f\left({\mathrm{\Delta }}_{\text{L}}\right)}{I_{\text{gL}}{\omega }_n^2}+\frac{c_{\text{gL}\theta y}{r}_{\text{gL}}\sin \hspace{0.17em}{\beta }_{\text{L}}\cos \hspace{0.17em}\psi {\dot{\theta }}_{\text{ygL}}}{I_{\text{gL}}{\omega }_n}+\frac{k_{\text{gL}\theta y}{r}_{\text{gL}}\sin \hspace{0.17em}{\beta }_{\text{L}}\cos \hspace{0.17em}\psi {\theta }_{\text{ygL}}}{I_{\text{gL}}{\omega }_n^2}\\ +\frac{T_{\text{pL}}{r}_{\text{pL}}\cos \hspace{0.17em}{\beta }_{\text{L}}}{I_{\text{pL}}{b}_c{\omega }_n^2}+\frac{T_{\text{gL}}{r}_{\text{gL}}\cos \hspace{0.17em}{\beta }_{\text{L}}}{I_{\text{gL}}{b}_c{\omega }_n^2}-\frac{{\ddot{e}}_{\text{L}}\left(t\right)}{b_c{\omega }_n^2}\end{array}$$

(25)

The differential equation on the meshing line of helical gear pair at the right end of double-helical gear is as follows

$$\{\begin{array}{l}{\ddot{\mathrm{\Delta }}}_{\text{R}}=\left[\left({\ddot{\overline{x}}}_{\text{pR}}-{\ddot{\overline{x}}}_{\text{gR}}\right)\sin \hspace{0.17em}\psi +\left({\ddot{\overline{y}}}_{\text{pR}}-{\ddot{\overline{y}}}_{\text{gR}}\right)\cos \hspace{0.17em}\psi \right]\cos \hspace{0.17em}{\beta }_{\text{R}}+\left[\left({\ddot{\overline{z}}}_{\text{pR}}-{\ddot{\overline{z}}}_{\text{gR}}\right)\right]\sin \hspace{0.17em}{\beta }_{\text{R}}-\frac{c_{\text{mR}}{c}^2{\beta }_{\text{R}}{\dot{\mathrm{\Delta }}}_{\text{R}}}{m_{\text{c}}{\omega }_n}-\frac{k_{\text{mR}}{c}^2{\beta }_{\text{R}}f\left({\mathrm{\Delta }}_{\text{R}}\right)}{m_{\text{c}}{\omega }_n^2}\\ +\frac{c_{\text{pt}}{r}_{\text{pR}}^2\cos {\beta }_{\text{R}}}{I_{\text{pR}}{\omega }_n}{\dot{\mathrm{\Delta }}}_{\text{P}}+\frac{k_{\text{pt}}{r}_{\text{pR}}^2\cos {\beta }_{\text{R}}}{I_{\text{pR}}{\omega }_n^2}{\mathrm{\Delta }}_{\text{P}}-\frac{c_{\text{gt}}{r}_{\text{gR}}^2\cos \hspace{0.17em}{\beta }_{\text{R}}}{I_{\text{gR}}{\omega }_n}{\dot{\mathrm{\Delta }}}_{\text{G}}-\frac{k_{\text{gt}}{r}_{\text{gR}}^2\cos \hspace{0.17em}{\beta }_{\text{R}}}{I_{\text{gR}}{\omega }_n^2}{\mathrm{\Delta }}_{\text{G}}-\frac{c_{\text{mR}}{r}_{\text{pR}}^2{s}^2{\beta }_{\text{R}}{s}^2\psi {\dot{\mathrm{\Delta }}}_{\text{R}}}{I_{\text{pR}}{\omega }_n}\\ -\frac{k_{\text{mR}}{r}_{\text{pR}}^2{s}^2{\beta }_{\text{R}}{s}^2\psi f\left({\mathrm{\Delta }}_{\text{R}}\right)}{I_{\text{pR}}{\omega }_n^2}-\frac{c_{\text{pR}\theta x}{r}_{\text{pR}}\sin \hspace{0.17em}{\beta }_{\text{R}}\sin \hspace{0.17em}\psi {\dot{\theta }}_{\text{xpR}}}{I_{\text{pR}}{\omega }_n}-\frac{k_{\text{pR}\theta x}{r}_{\text{pR}}\sin \hspace{0.17em}{\beta }_{\text{R}}\sin \hspace{0.17em}\psi {\theta }_{\text{xpR}}}{I_{\text{pR}}{\omega }_n^2}-\frac{c_{\text{mR}}{r}_{\text{gR}}^2{s}^2{\beta }_{\text{R}}{s}^2\psi {\dot{\mathrm{\Delta }}}_{\text{R}}}{I_{\text{gR}}{\omega }_n}\\ -\frac{k_{\text{mR}}{r}_{\text{gR}}^2{s}^2{\beta }_{\text{R}}{s}^2\psi f\left({\mathrm{\Delta }}_{\text{R}}\right)}{I_{\text{gR}}{\omega }_n^2}+\frac{c_{\text{gR}\theta x}{r}_{\text{gR}}\sin \hspace{0.17em}{\beta }_{\text{R}}\sin \hspace{0.17em}\psi {\dot{\theta }}_{\text{xgR}}}{I_{\text{gR}}{\omega }_n}+\frac{k_{\text{gR}\theta x}{r}_{\text{gR}}\sin \hspace{0.17em}{\beta }_{\text{R}}\sin \hspace{0.17em}\psi {\theta }_{\text{xgR}}}{I_{\text{gR}}{\omega }_n^2}-\frac{c_{\text{mR}}{r}_{\text{pR}}^2{s}^2{\beta }_{\text{R}}{c}^2\psi {\dot{\mathrm{\Delta }}}_{\text{R}}}{I_{\text{pR}}{\omega }_n}\\ -\frac{k_{\text{mR}}{r}_{\text{pR}}^2{s}^2{\beta }_{\text{R}}{c}^2\psi f\left({\mathrm{\Delta }}_{\text{R}}\right)}{I_{\text{pR}}{\omega }_n^2}-\frac{c_{\text{pR}\theta y}{r}_{\text{pR}}\sin \hspace{0.17em}{\beta }_{\text{R}}\cos \hspace{0.17em}\psi {\dot{\theta }}_{\text{ypR}}}{I_{\text{pR}}{\omega }_n}-\frac{k_{\text{pR}\theta y}{r}_{\text{pR}}\sin \hspace{0.17em}{\beta }_{\text{R}}\cos \hspace{0.17em}\psi {\theta }_{\text{ypR}}}{I_{\text{pR}}{\omega }_n^2}-\frac{c_{\text{mR}}{r}_{\text{gR}}^2{s}^2{\beta }_{\text{R}}{c}^2\psi {\dot{\mathrm{\Delta }}}_{\text{R}}}{I_{\text{gR}}{\omega }_n}\\ -\frac{k_{\text{mR}}{r}_{\text{gR}}^2{s}^2{\beta }_{\text{R}}{c}^2\psi f\left({\mathrm{\Delta }}_{\text{R}}\right)}{I_{\text{gR}}{\omega }_n^2}+\frac{c_{\text{gR}\theta y}{r}_{\text{gR}}\sin \hspace{0.17em}{\beta }_{\text{R}}\cos \hspace{0.17em}\psi {\dot{\theta }}_{\text{ygR}}}{I_{\text{gR}}{\omega }_n}+\frac{k_{\text{gR}\theta y}{r}_{\text{gR}}\sin \hspace{0.17em}{\beta }_{\text{R}}\cos \hspace{0.17em}\psi {\theta }_{\text{ygR}}}{I_{\text{gR}}{\omega }_n^2}\\ +\frac{T_{\text{pR}}{r}_{\text{pR}}\cos \hspace{0.17em}{\beta }_{\text{R}}}{I_{\text{pR}}{b}_c{\omega }_n^2}+\frac{T_{\text{gR}}{r}_{\text{gR}}\cos \hspace{0.17em}{\beta }_{\text{R}}}{I_{\text{gR}}{b}_c{\omega }_n^2}-\frac{{\ddot{e}}_{\text{R}}\left(t\right)}{b_c{\omega }_n^2}\end{array}$$

(26)

The relative torsional linear displacement differential equation of adjacent coaxial helical gears is as follows

$$\{\begin{array}{l}{\ddot{\mathrm{\Delta }}}_{\text{P}}=\frac{T_{\text{pL}}}{I_{\text{pL}}{b}_{\text{c}}{\omega }_n^2}-\frac{T_{\text{pR}}}{I_{\text{pL}}{b}_{\text{c}}{\omega }_n^2}-\frac{k_{\text{mL}}f\left({\mathrm{\Delta }}_{\text{L}}\right)\cos \hspace{0.17em}{\beta }_{\text{L}}{R}_{\text{pL}}}{I_{\text{pL}}{\omega }_n^2}-\frac{c_{\text{mL}}{\dot{\mathrm{\Delta }}}_{\text{L}}\cos \hspace{0.17em}{\beta }_{\text{L}}{R}_{\text{pL}}}{I_{\text{pL}}{\omega }_n}+\frac{k_{\text{mR}}f\left({\mathrm{\Delta }}_{\text{R}}\right)\cos \hspace{0.17em}{\beta }_{\text{R}}{R}_{\text{pR}}}{I_{\text{pL}}{\omega }_n^2}\\ +\frac{c_{\text{mR}}{\dot{\mathrm{\Delta }}}_{\text{R}}\cos \hspace{0.17em}{\beta }_{\text{R}}{R}_{\text{pR}}}{I_{\text{pL}}{\omega }_n}-\frac{k_{\text{pt}}{\mathrm{\Delta }}_{\text{p}}{R}_{\text{pL}}}{I_{\text{pL}}{\omega }_n^2}-\frac{c_{\text{pt}}{\dot{\mathrm{\Delta }}}_{\text{p}}{R}_{\text{pL}}}{I_{\text{pL}}{\omega }_n}-\frac{k_{\text{pt}}{\mathrm{\Delta }}_{\text{p}}{R}_{\text{pR}}}{I_{\text{pL}}{\omega }_n^2}-\frac{c_{\text{pt}}{\dot{\mathrm{\Delta }}}_{\text{p}}{R}_{\text{pR}}}{I_{\text{pL}}{\omega }_n}\\ {\ddot{\mathrm{\Delta }}}_{\text{G}}=\frac{-T_{\text{gL}}}{I_{\text{gL}}{b}_{\text{c}}{\omega }_n^2}+\frac{T_{\text{gR}}}{I_{\text{gL}}{b}_{\text{c}}{\omega }_n^2}+\frac{k_{\text{mL}}f\left({\mathrm{\Delta }}_{\text{L}}\right)\cos \hspace{0.17em}{\beta }_{\text{L}}{R}_{\text{gL}}}{I_{\text{gL}}{\omega }_n^2}+\frac{c_{\text{mL}}{\dot{\mathrm{\Delta }}}_{\text{L}}\cos \hspace{0.17em}{\beta }_{\text{L}}{R}_{\text{gL}}}{I_{\text{gL}}{\omega }_n}-\frac{k_{\text{mR}}f\left({\mathrm{\Delta }}_{\text{R}}\right)\cos \hspace{0.17em}{\beta }_{\text{R}}{R}_{\text{gR}}}{I_{\text{gL}}{\omega }_n^2}\\ -\frac{c_{\text{mR}}{\dot{\mathrm{\Delta }}}_{\text{R}}\cos \hspace{0.17em}{\beta }_{\text{R}}{R}_{\text{gR}}}{I_{\text{gL}}{\omega }_n}-\frac{k_{\text{gt}}{\mathrm{\Delta }}_{\text{g}}{R}_{\text{gL}}}{I_{\text{gL}}{\omega }_n^2}-\frac{c_{\text{gt}}{\dot{\mathrm{\Delta }}}_{\text{g}}{R}_{\text{gL}}}{I_{\text{gL}}{\omega }_n}-\frac{k_{\text{gt}}{\mathrm{\Delta }}_{\text{g}}{R}_{\text{gR}}}{I_{\text{gL}}{\omega }_n^2}-\frac{c_{\text{gt}}{\dot{\mathrm{\Delta }}}_{\text{g}}{R}_{\text{gR}}}{I_{\text{gL}}{\omega }_n}\end{array}$$

(27)

The stiffness terms are as follows

$$k_{\text{m}j}\left(t\right)=1-{\displaystyle \sum _{l=1}^L{\epsilon }_l\cos \left(l\Omega \overline{t}\right)}\hspace{1em}j=\text{L},\hspace{0.17em}\text{R}$$

(28)

The comprehensive transmission error term in the excitation vector is as follows

$${\overline{e}}_j^{″}\left(t\right)=1-{\displaystyle \sum _{l=1}^L{\epsilon }_l{\left(l\Omega \right)}^2{e}_l\sin \left(l\Omega \overline{t}\right)}\hspace{1em}j=\text{L},\hspace{0.17em}\text{R}$$

(29)

For the piecewise backlash function, taking b _c =b₀ in the dimensionless process, the backlash model is transformed into

$$f\left({\Delta }_j\right)=\{\begin{array}{l}{\Delta }_j-\left(1+W\right)\hspace{1em}{\Delta }_j>1+W\\ 0-\left(1+W\right)\hspace{1em}\le {\Delta }_j\le 1+W\\ {\Delta }_j+\left(1+W\right)\hspace{1em}{\Delta }_j<-\left(1+W\right)\end{array}$$

(30)

After eliminating the rigid body displacement and dimensionalizing, the displacement vector of the system becomes

$$\overline{q}=\left\{{\overline{x}}_{\text{pL}},{\overline{y}}_{\text{pL}},{\overline{z}}_{\text{pL}},{\overline{x}}_{\text{gL}},{\overline{y}}_{\text{gL}},{\overline{z}}_{\text{gL}},{\overline{x}}_{\text{pR}},{\overline{y}}_{\text{pR}},{\overline{z}}_{\text{pR}},{\overline{x}}_{\text{gR}},{\overline{y}}_{\text{gR}},{\overline{z}}_{\text{gR}},{\Delta }_{\text{L}},{\Delta }_{\text{R}},{\Delta }_{\text{P}},{\Delta }_{\text{G}}\right\}$$

(31)

Model solving and influence analysis

The Runge–Kutta method is used to solve the BTAP coupling nonlinear dynamic Equations. The input torque is 1400 Nm. The other parameters are shown in Tables 1 and 2

Table 1.

System parameters.

Parameter	Driving gear	Driven gear
Number of teeth z	25	75
Helix angle β/(°)	24.43
Pressure angle/(°)	20
Modulus/(mm)	6
Tooth width b/(mm)	55
Moment of inertia/(Kg·m2)	0.012	0.7279
Quality/(Kg)	7.11	64.75
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

Table 2.

System stiffness and damping parameters.

Parameter item	Parameter value	Parameter item	Parameter value
kmean/N⋅m−1	1.8456×109	kgLx, kgRx/N⋅m−1	1.6741×109
kpLxkpRx/N⋅m−1	1.5424×109	kgLy, kgRy/N⋅m−1	4.5124×108
kpLy, kpRy/N⋅m−1	1.5746×109	kgLz, kgRz/N⋅m−1	1.3454×108
cpLx, cpRx/N⋅(m⋅Kg)−1	1.1745×105	cgLx, cgRx/N⋅(m⋅kg)−1	1.1500×105
cpLy, cpRy/N⋅(m⋅Kg)−1	3.7145×106	cgLy, cgRy/N⋅(m⋅kg) −1	3.6500×106
kpLθx, kpRθx/N⋅m⋅rad−1	2.3809×109	cgLz, cgRz/N⋅(m⋅Kg)−1	0.2416×105
kpLθy, kpRθy/N⋅m⋅rad−1	3.3607×108	kgLθx, kgRθx/N⋅m⋅rad−1	2.3988×109
		kgLθy, kgRθy/N⋅m⋅rad−1	3.3607×108

.

The basic system parameters are ζ=0.063, ɛ₁=0.16, ɛ₂=0.07, ɛ₃=0.03, e₁=0.03, e₂=0.01, e₃=0.01, W=0. Taking the left helical gear of double-helical gear system as an example, the phase diagram and time history diagram of primary resonance response and sub-harmonic response are obtained, as shown in Figure 3.

Figure 3.

Phase diagram and time history diagram of the system at different frequencies. (a) Ω=1 (b) Ω=1/2 (c) Ω=1/3.

Figure 3 shows that when Ω=1, the phase diagram is approximate to ellipse, the system has a single period harmonic response, and the gear is in a unilateral impact state. When Ω=1/2, Ω=1/3, the phase diagram is non-circular closed curve, the system response is subharmonic response, and the gear is in shock free state. The superharmonic response of the system is caused by the time-varying meshing stiffness and the high-order harmonic component under dynamic excitation.

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 curve of the system is obtained, and the results are shown in Figure 4.

Figure 4.

Amplitude frequency characteristic curves of the system.

It can be seen from Figure 4, when the dimensionless excitation frequency is 0.89, the amplitude frequency characteristic curve of double-helical gear system jumps, 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. The double-helical gear transmission system also has strong nonlinear dynamic characteristics caused by backlash.

The basic system parameters are ζ=0.063, ɛ₁=0.16, ɛ₂=0.07, ɛ₃=0.03, e₁=0.03, e₂=0.01, e₃=0.01. The influence of different constant backlash on the amplitude frequency characteristics of the system is shown in Figure 5.

Figure 5.

Amplitude frequency characteristic curves of the system with different constant backlash.

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. (5) W=0.6. Figure 5 shows that with the increase of the constant backlash, the vibration displacement amplitude of the gear system increases and the jumping phenomenon still exists. And the step frequency of the curve jumping has a weak right tilt phenomenon, At the same time, the amplitude frequency characteristic curve has smaller peak value at Ω=1/2 and Ω=1/3 of the dimensionless frequency, which corresponds to the sub-harmonic response of the system.

Influence of system parameters on amplitude frequency characteristics

Influence of time-varying meshing stiffness on amplitude frequency curve of system

The basic system parameters are ζ=0.063, e₁=0.03, e₂=0.01, e₃=0.01. The influence of time-varying meshing stiffness on the amplitude frequency curve of the system under different constant backlash W=0 and W=0.6 is shown in Figure 6. 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. (3) ɛ₁=0, ɛ₂=0, ɛ₃=0.

Figure 6.

Influence of time-varying meshing stiffness on amplitude frequency characteristic curve of system. (a) W=0 (b) W=0.6.

Figure 6 shows that the amplitude of amplitude frequency characteristic curve decreases with the decrease of time-varying meshing stiffness fluctuation coefficient, which is more obvious at lower frequencies. At the same time, with the decrease of the fluctuation coefficient of time-varying meshing stiffness, the jumping phenomenon of the system will gradually decrease or even disappear. Compared with W=0, when the constant clearance is 0.6, the system has shown linear vibration characteristics when the meshing stiffness fluctuation coefficient is 0.13, and the jumping phenomenon does not exist. When the meshing stiffness fluctuation coefficient is 0.16 and 0, the system state does not change, but the vibration displacement increases.

Influence of damping on amplitude frequency curve of system

The basic system parameters are ɛ₁=0.16, ɛ₂=0.07, ɛ₃=0.03, e₁=0.03, e₂=0.01, e₃=0.01. Under different constant backlash W=0 and W=0.6, the amplitude frequency characteristic curves of the system under different damping ratios are shown in Figure 7. Figure 7 shows that the damping can inhibit the step of the amplitude frequency characteristic curve of the system and the amplitude of nonlinear vibration displacement. Increasing the system damping can effectively control the nonlinear dynamic response of the gear system. Compared with W=0, when the constant backlash is W=0.6, the system shows linear vibration performance when the damping ratio is 0.07, and the jump phenomenon does not exist. When the damping ratio is 0.063 and 0.1, the system state does not change, but the vibration amplitude increases.

Figure 7.

Influence of damping on amplitude frequency characteristic curve of system. (a) W=0 (b) W=0.6.

Influence of comprehensive transmission error on amplitude frequency curve of system

The basic system parameters are ζ=0.063, ɛ₁=0.16, ɛ₂=0.07, ɛ₃=0.03. Under different constant backlash W=0 and W=0.6, the influence of comprehensive transmission error on the amplitude frequency curve of the system is shown in Figure 8. Three kinds of comprehensive transmission error amplitudes are considered: (1) e₁=0.1, e₂=0.05, e₃=0.03. (2) e₁=0.03, e₂=0.01, e₃=0.01. (3) e₁=0, e₂=0, e₃=0.

Figure 8.

Influence of comprehensive transmission error on amplitude frequency characteristic curve of system. (a) W=0 (b) W=0.6.

Figure 8 shows that the amplitude frequency characteristic curve of the system changes greatly with the change of the amplitude of the comprehensive transmission error. When e₁=0, e₂=0, e₃=0, the amplitude frequency curve still has jump phenomenon, indicating that the comprehensive transmission error 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 amplitude frequency curve of the system has a greater jump, and the larger the amplitude of static transmission error, the greater the amplitude of amplitude frequency characteristic curve. With the increase of backlash, the vibration amplitude increases, and the system state does not change.

Influence of external load excitation on amplitude frequency curve of system

The basic system parameters are ζ=0.063, ɛ₁=0.16, ɛ₂=0.07, ɛ₃=0.03, e₁=0.03, e₂=0.01, e₃=0.01. Under different constant backlash W=0 and W=0.6, the influence of external load excitation on the amplitude frequency curve of the system is shown in Figure 9.

Figure 9.

Influence of external load excitation on amplitude frequency characteristic curve of system. (a) W=0 (b) W=0.6.

Figure 9 shows that with the increase of external load excitation, the amplitude of the amplitude frequency characteristic curve of the system increases, and the step frequency of the curve jump has a weak right inclination phenomenon. With the increase of external load excitation, the jumping phenomenon will disappear, and the system shows linear vibration characteristics. Increasing external load excitation can effectively control the nonlinear dynamic response of gear system. With the increase of backlash, the vibration displacement increases, and the system state does not change.

Test verification

The physical diagram of double helical gear vibration test device is shown in Figure 10, including control box, double helical gear vibration test bench, data collector, analysis and post-processing software, etc.

Figure 10.

Vibration testing equipment for double-helical gear.

The structure diagram of double-helical gear vibration test equipment is shown in Figure 11.

Figure 11.

Structural diagram of double-helical gear vibration testing equipment. 1- Control box. 2- AC servo motor. 3, 6, 8- Elastic coupling. 4- Double-helical gear box. 5- Acceleration sensor. 7- Torque sensor. 9- Magnetic powder brake. 10- Data collector. 11- Computer. 12- Emergency stop switch.

The structure of double-helical gear vibration test-bed is shown in Figure 12. The test-bed is driven by a motor which can adjust the speed through a frequency converter, and its maximum speed is 1500 r/min. Four acceleration sensors are respectively installed on the box close to the rolling bearing. There is a dial on the end cover of the gearbox, which can adjust the center distance between the driving and driven gear.

Figure 12.

Vibration test bench for double-helical gear.

As shown in Figure 13, adjust the average backlash to 50 μm by adjusting the knob, and apply solid lubricating oil on the gears to reduce the friction between the tooth surfaces. For the acceleration signal measured by the experiment, the noise is reduced and the signal is amplified. At the same time, the displacement signal is obtained by quadratic integration.

Figure 13.

Backlash adjusting disc.

The FFT spectrum of double-helical gear system at Ω = 1/3 of the dimensionless frequency is obtained by theoretical calculation, as shown in Figure 14(a). In order to verify the theoretical reliability, the spectrum diagram of the test data at the corresponding speed of 400 r/min is measured, as shown in Figure 14(b).

Figure 14.

FFT Spectrum at the speed of 400 r/min. (a) Theoretical data (b) Experimental data.

As shown in Figure 14(b), when the rotating speed is 400r/min, the vibration displacement spectrum measured in the test mainly includes 166.7 Hz of meshing frequency, 333.4 Hz twice of meshing frequency and 506.73 Hz three times of meshing frequency. Comparing the theoretical calculation results with the experimental measurement results, the theoretical calculation results are basically the same as the experimental measurement results.

The FFT spectrum of double-helical gear system at Ω=1/2 of the dimensionless frequency is obtained by theoretical calculation, as shown in Figure 15(a). In order to verify the theoretical reliability, the spectrum of test data at the corresponding speed of 600r/min is measured, as shown in Figure 15(b).

Figure 15.

FFT Spectrum at the speed of 600 r/min. (a) Theoretical data (b) Experimental data.

Figure 15(b) shows that when the rotating speed is 600r/min, the vibration displacement spectrum measured in the test mainly includes 250 Hz of meshing frequency, 500 Hz twice of meshing frequency and 750 Hz three times of meshing frequency. Comparing the theoretical calculation results with the experimental measurement results, the theoretical calculation results are basically the same as the experimental measurement results, and the difference between them is that there is no third harmonic generation of meshing frequency in the theoretical calculation.

The FFT spectrum of double-helical gear system at dimensionless frequency Ω=1 is obtained by theoretical calculation, as shown in Figure 16(a). In order to verify the theoretical reliability, the spectrum of test data at corresponding speed of 1200r/min is measured, as shown in Figure 16(b).

Figure 16.

FFT Spectrum at the speed of 1200 r/min. (a) Theoretical data (b) Experimental data.

As shown in Figure 16(b), when the rotating speed is 1200r/min, the meshing frequency is 480.0 Hz, and an obvious sideband is formed near the meshing frequency, and the interval of sideband is twice of the rotating frequency. Comparing the theoretical calculation results with the experimental measurement results, the theoretical calculation results are basically the same as the experimental measurement results, and the difference between them is that there is no sideband near the meshing frequency in the theoretical calculation.

The test results show that the frequency components of the theoretical calculation results and the test results are the same, but there are manufacturing errors, installation errors and measurement errors in the test, which will produce certain errors with the numerical simulation. These factors will lead to different degrees of frequency component response.

When the test speed is 1200r/min, the displacement spectrum under different backlash b is shown in Figure 17.

Figure 17.

Displacement spectrum under different backlash at the speed of 1200 r/min. The variation of displacement amplitude with rotation speed under different backlash is shown in Figure 18. (a) b=50.77 μm (b) b=51.02 μm (c) b=51.20 μm (d) b=51.27 μm.

Figure 18.

Variation of displacement amplitude with rotating speed under different backlash.

The results show that the displacement amplitude increases with the increase of backlash at the same speed. Under the same backlash, the displacement amplitude increases gradually with the increase of rotational speed from 400 r/min to 1400 r/min. The experimental results are the same as the theoretical results, which further verify the reliability of the theory.

Conclusions

A nonlinear BTAP coupling dynamic modeling of double-helical gear is established considering backlash, time-varying meshing stiffness, damping, comprehensive transmission error and external load excitation, which improves the accuracy of the traditional model and is closer to the actual working state of double-helical gear. At the same time, the influence of each parameter on the frequency response characteristics of the system is analyzed, the results are as follows

1. The double helical gear system has the strong nonlinear characteristics of jump caused by backlash, as well as the main resonance response and subharmonic response. The generation of super harmonic response solution is caused by time-varying mesh stiffness and higher harmonic component under dynamic excitation.

2. With the increase of backlash, the vibration displacement amplitude of the gear system increases and the jumping phenomenon remains unchanged.

3. The amplitude frequency response of the system is excited by time-varying meshing stiffness and comprehensive transmission error, while the amplitude frequency response is constrained by damping and external load excitation. The increase of constant backlash will increase the amplitude of vibration displacement and has little effect on the change of system state.

4. The frequency components of the theoretical results are basically consistent with the experimental results. At the same speed, the system amplitude increases with the increase of backlash. Under the condition of the same backlash, the amplitude of the system increases with the increase of the speed in a certain range. The experimental results verify the correctness of the theory.