Yaguo Lei^{1} , Xiao Yang^{4} , Zongyao Liu^{3} , Delong Wang^{2}
^{1, 4, 3, 2}Key Laboratory of Education Ministry for Modern Design and RotorBearing System, School of Mechanical Engineering, Xi’an Jiaotong University, Xi’an, 710049, China
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 14, 2017, p. 334339.
https://doi.org/10.21595/vp.2017.19139
Received 9 September 2017; accepted 17 September 2017; published 21 October 2017
Copyright © 2017  JVE International Ltd.
Timevarying gearmesh stiffness (TVGS) is the main cause of gear vibration, and its accuracy affects the responses of dynamic models. An exponential curve model based on the Saint Venant’s Principle is proposed to calculate the gearmesh stiffness of cracked spur gears in this paper. With the proposed model, the TVGS under the circumstances of healthy condition and four crack cases are computed, whose results have a good agreement with those of finite element method (FEM). Therefore, the exponential curve model can be used to estimate the TVGS and an alternative to FEM in gearmesh stiffness calculation is provided.
Keywords: timevarying gearmesh stiffness, tooth crack, exponential curve model.
Gearboxes are extensively employed in many industrial applications because of their accurate transmission ratio, high power density and large reduction. However, due to high service loads, harsh operating conditions and fatigue, faults may develop in gears, such as root crack and spalling [1]. Therefore, it is of vital importance to detect initial faults in gear transmission systems for the purpose of condition monitoring [2]. Generally, tooth damage often shows a reduction in gear mesh stiffness and causes the effective capability for bearing load decreasing. Therefore, the severity of tooth damage can be assessed by investigating the reduction of TVGS in gear meshing [3, 4]. In the field of damage mechanism analysis, dynamic modeling of a gear system with simulated faults has been a promising way.
The methods of TVGS calculation in previous studies can generally be categorized as finite element method (FEM) [5], experimental method [6] and analytical method [711], among which FEM and experimental method both require precise experiment facilities, though they have a higher accuracy. However, the analytical method provides an alternative to the FEM and the experimental method for its high efficiency to obtain gearmesh stiffness. The potential energy method [7], which belongs to analytical method, is widely adopted in gearmesh stiffness calculation. Tian [8] added the shear stiffness to polish the method, in which the concept of limiting line is proposed when calculating the effective area moment of inertia and effective section area using potential energy method. Liang [9] estimated the gearmesh stiffness of a planetary gear set using a straight limiting line. Chen [10] took nonuniformly distribution of crack along tooth width into consideration. Ma [11] evaluated the gearmesh stiffness of a cracked spur gear set by regarding the limiting line as a parabolic curve.
When a crack occurs and propagates, the effective area moment of inertia and effective section area can be calculated using potential energy method based on a limiting line. The straight line model [9] and parabolic curve model [11] are used to calculate the reduction of bend and shear stiffness in aforementioned works. However, it is obviously inaccurate because the existence of a crack in the beam memorably changes the stress and strain fields in the vicinity of crack, leading to stress concentration at the crack tip. Christides [12] pointed out that according to the Saint Venant’s Principle, the stress concentration at the crack tip decays exponentially with the increasing of distance from the crack location. Carneiro [13] extend this model to Timoshenko beam with shear stiffness considered. Hou [14] summarily reviewed crack beam models and indicated that the exponential curve model coincides best with the practical distribution of stiffness around the crack. Therefore, an exponential curve model for TVGS calculation of crack tooth is proposed in this paper.
In this study, tooth crack is modeled and equations are derived to estimate the gearmesh stiffness of cracked gear set using the potential energy method. Then, a case study is given to illustrate the crack effect on TVGS with an exponential curve model, in which the crack severity levels are modeled with 20 %, 40 %, 60 % and 80 % of the total length of the full through crack, respectively. Finally, the proposed method is compared with FEM to validate its accuracy.
The tooth shape between the base circle and root circle of an involute spur gear is determined by a cutter tip. In order to simulate the realistic tooth profile, the transition curve is assumed as an arc $AB$ in this paper, as shown in Fig. 1. To calculate the radius of arc $AB$, $A{O}_{1}$ is assumed to be perpendicular to $x$ axis, because of the thin tooth thickness. The radius $RR$ of arc $AB$ satisfies the following equation:
where ${\alpha}_{2}$ is the half of the base tooth angle shown in Fig. 1. Then ${h}_{{x}_{1}}$ can be calculated as:
In this paper, the potential energy method is adopted to calculate the TVGS of spur gears [7]. The gearmesh stiffness includes five components: Hertzian contact stiffness ${k}_{h}$, bend stiffness ${k}_{b}$, axial compressive stiffness ${k}_{a}$, shear stiffness ${k}_{s}$ and foundation stiffness ${k}_{f}$.Therefore, the total gearmesh stiffness $k$ can be expressed as Eq. (3):
where the subscripts 1 and 2 indicate the pinion and gear, respectively; $j=$1 and 2 represent the first and the second pair of meshing teeth, respectively.
The parameters of the gear set are as follows: the tooth number of pinion ${N}_{p}=$20; the tooth number of gear ${N}_{g}=$40; module $m=$1 mm; width of tooth $L=$10 mm; the hub radius of pinion ${R}_{int,p}=$ 4 mm; the hub radius of gear ${R}_{int,g}=$5 mm; Poisson’s ratio $\upsilon =$ 0.3, Young’s modulus $E=$2.068×10^{11} Pa, pressure angle on pitch circle $\alpha \text{'}=$ 20°. And the tooth is modeled as a nonuniform cantilever beam on the root circle, as described in Fig. 1. Then the gearmesh stiffness can be calculated as follows:
where $G$ and $E$ are the shear modulus and Young’s modulus, respectively; $L$ is the tooth width and $\upsilon $ is Poisson’s ratio; ${I}_{x}$ and ${A}_{x}$ represent the area moment of inertia and area of the section where the distance from the base circle is $x$, ${I}_{x}=(2{h}_{x}{)}^{3}L/12$, ${A}_{x}=2{h}_{x}L$; ${I}_{{x}_{1}}$ and ${A}_{{x}_{1}}$ represent the area moment of inertia and area of the section where the distance from the base circle is ${x}_{1}$, ${I}_{x1}=(2{h}_{x1}{)}^{3}L/12$, ${A}_{x1}=2{h}_{x1}L$. The details about ${\mu}_{f}$, ${S}_{f}$${L}^{*}$, ${M}^{*}$, ${P}^{*}$, and ${Q}^{*}$ can be found in Ref. [11]. The TVGS of a healthy spur gear set is shown in Fig. 4.
Fig. 1. Cantilever beam model for the spur gear tooth
The occurrence of tooth crack reduces the bend stiffness and shear stiffness, while the axial compressive stiffness and Hertzian contact stiffness remain the same [7]. An exponential curve model based on the Saint Venant’s Principle is proposed to describe the variation of the effective area moment of inertia ${I}_{x}$ and effective section area ${A}_{x}$. To illustrate the model, a rectangular section beam is selected to explain how the stress decays exponentially in the vicinity of the crack. And the differential equations of the beam can be presented as:
where ${C}_{1}={I}_{0}/{I}_{c}1$; ${I}_{0}$ is the area moment of inertia, ${I}_{0}={H}^{3}{L}_{1}/12$; ${L}_{1}$ and $H$ are the width and height of the beam section; ${I}_{c}$ represents the area moment of inertia of the cracked section, ${I}_{c}=(H{h}_{c}{)}^{3}{L}_{1}/12$; $M$ is the mass of a unit; ${x}_{c}$ and ${h}_{c}$ is the location and the depth of the crack, as shown in Fig. 2; $\lambda $ is an experimental coefficient, and $\lambda =$ 0.667 [12].
Based on Eq. (9) and Eq. (10), the effective area moment of inertia and section area of a rectangular beam, illustrated in Fig. 2, can be expressed as [12]:
where ${C}_{2}={A}_{0}/{A}_{c}1$; ${A}_{0}$ is the area of the undamaged beam, ${A}_{0}=H{L}_{1}$; ${A}_{c}$ is the area moment of inertia of the cracked section, ${A}_{c}=(H{h}_{c}){L}_{1}$.
Fig. 2. Effective loading area of cracked beam
Based on the aforementioned cracked beam model, an exponential curve model is proposed to compute the TVGS of cracked gears. The exponentially decayed curve is shown in Fig. 3, then $\stackrel{}{I}\left(x\right)$ and $\stackrel{}{A}\left(x\right)$ can be expressed as follows:
where ${x}_{Q}$ is the coordinate of crack tip $Q$ in $x$ direction; ${H}_{c}$ is the tooth thickness at the location of the crack tip; ${I}_{0}\text{,}$${I}_{c}\text{,}$${A}_{0}$ and ${A}_{c}$ can be expressed as ${I}_{0}={{H}_{c}}^{3}L/12\text{,}$${I}_{c}=[{H}_{c}({q}_{1}+{q}_{2})\mathrm{s}\mathrm{i}\mathrm{n}{\alpha}_{c}{]}^{3}L/12$, ${A}_{0}={H}_{c}L$, ${A}_{c}=[{H}_{c}({q}_{1}+{q}_{2}\left)\mathrm{s}\mathrm{i}\mathrm{n}{\alpha}_{c}\right]L$, respectively; ${q}_{1}$ and ${q}_{2}$ are the crack lengths when the crack is upon the central line and beyond the central line, respectively; ${\alpha}_{c}$ represents the direction of crack propagation, and there is ${\alpha}_{c}=$ 45°.
The bend stiffness and shear stiffness of the cracked tooth are presented as follows:
The location where crack starts is determined by 30° tangent method, shown as $\phi =$ 30° in Fig. 3. Based on the proposed method, four cases with different crack lengths are calculated, and the parameters are listed in Table 1. The TVGS of four cases is given in Fig. 4.
Table 1. Parameters of four crack cases
Crack level

20 %

40 %

60 %

80 %

Crack length

${q}_{1}=$0.6 mm,
${q}_{2}=$0

${q}_{1}=$1.1 mm,
${q}_{2}=$0

${q}_{1}={q}_{max}$,
${q}_{2}=$0.3 mm

${q}_{1}={q}_{max}$,
${q}_{2}=$0.8 mm

To validate the effectiveness of the TVGS calculated by potential energy method in the proposed model, a 3D finite element model using Solid185 elements with three teeth is adopted. A straight line with an angle of 45° is used to simulate the crack propagation on the pinion tooth, and four crack cases are modeled with the parameters in Table 1. In the finite element model, the load is applied on the hub bore of pinion and the hub bore of gear is fully constrained. Each node of the elements of the pinion is constrained to move only in the direction which is tangential to the hub bore. When the deformation in the direction of action line is acquired from FEM, the gearmesh stiffness can be obtained. Then the stiffness results of healthy condition and four crack cases by FEM are displayed in Fig. 5.
Fig. 3. Exponential curve model on cracked tooth: a) ${q}_{1}+{q}_{2}\le {q}_{\mathrm{m}\mathrm{a}\mathrm{x}}$; b) ${q}_{1}+{q}_{2}>{q}_{\mathrm{m}\mathrm{a}\mathrm{x}}$
a)
b)
Fig. 4. TVGS of four cases of cracked tooth
Fig. 5. The comparison between analytical method and FEM: a) Healthy condition; b) 20 % crack; c) 40 % crack; d) 60 % crack; e) 80 % crack
In the perfect condition, the gearmesh stiffness results obtained from the potential energy method and the FEM match significantly well, which validates the accuracy of cantilever beam model on the root circle. With the growth of crack, the gearmesh stiffness results of the proposed exponential curve model agree well with those obtained from FEM. Overall, the gearmesh stiffness equations derived by exponential curve model in this paper can be used to estimate the gearmesh stiffness with root crack if appropriate error of gearmesh stiffness is allowed.
In this paper, an exponential curve model based on the Saint Venant’s Principle is proposed for the gearmesh stiffness calculation of cracked gears, and the TVGS results obtained from this model are compared with those of FEM. Conclusions are summarized as follows.
The transition deflection and foundation deflection have a significant influence on the TVGS. Thus, a accurate description of transition curve is of vital importance. In addition, the shape of limiting line plays an important role in calculate bend stiffness and shear stiffness, and the exponentially decayed curve based on the Saint Venant’s Principle is able to explain the stress concentration at the crack tip reasonably. Furthermore, the gearmesh stiffness results of exponential curve model are consistent well with those of FEM well. Therefore, the exponential curve model can be used to estimate the gearmesh stiffness and functions as an alternative to FEM in TVGS calculation