1. Introduction
The magnitude and variation of tooth pair compliance affect tooth loading and gear dynamics significantly [
1]. For an elastic gear, when a force
F is applied on a gear tooth, a tooth deflection
δ is generated along the line of force as illustrated in
Figure 1. The tooth deflection can be calculated using the finite element method (FEM) or analytical method [
2,
3]. With the development of computer science technology, the FEM can give very accurate results if the boundary conditions and the degree of discretization are properly defined [
4]. Meanwhile, many researchers have been developing analytical formulas to achieve the fast and accurate calculation of gear tooth deflections. Analytical formulas are easy to use and time-saving in computation [
5].
In Ref. [
1], gear teeth were assumed to be rigid in evaluating gear-body induced tooth deflections in their finite element model. This assumption will generate some errors. We will remove this assumption in our analysis and give a detailed explanation of how to develop a finite element model to evaluate gear-body induced tooth deflections.
In the existing analytical methods, tooth deflection is considered from two parts, the deflection from gear tooth (assuming rigid gear body) and the deflection caused by gear body. One popular analytical method to estimate the deflection from the first part is the potential energy method [
6,
7], in which the gear tooth is modeled as a cantilever beam of variable cross-section and the energy stored in the beam is considered as the summation of Hertzian contact energy, bending energy, shear energy, and axial compressive energy. Many researchers have applied this method to evaluate the time-varying meshing stiffness of gears [
8,
9,
10,
11]. The deflection induced by a gear body is also called tooth fillet/foundation deflection [
1]. Weber [
12], Attia [
13], Cornell [
14], and Sainsot and Velex [
1] developed analytical formulas to calculate the gear body-induced tooth deflections. Their analytical formulas all have the same format as Equation (1) but with different
L,
M,
P, and
Q values.
where
δ is tooth deflection under a force
F as shown in
Figure 1,
E represents Young’s modulus of gear body material,
denotes the gear pressure angle,
b is the tooth face width,
u represents the distance between the root circle and the intersection of the tooth center-line and the gear meshing action line (see
Figure 1),
h denotes the ratio of gear root circle radius (
Rf) and gear bore radius (
Ra),
θf represents the angle between tooth center-line and the junction with the root circle, and
Sf is the arc length corresponding to the angle
θf.
Table 1 gives a summary of the
L,
M,
P, and
Q values used in the four papers [
1,
12,
13,
14]. Authors in Refs. [
12,
13,
14] used constant
L,
M,
P, and
Q values in their formulas while Sainsot and Velex [
1] considered
L,
M,
P, and
Q as functions of
h and
θf. The formulas reported in Ref. [
1] were developed based on the theory of Muskhelishvili [
15] applied to circular elastic rings with the assumption of linear and constant stress variation at the root circle. It has been widely used in the gear meshing stiffness evaluation even for gears with a tooth crack [
16,
17,
18,
19]. Cirelli et al. [
20] studied non-linear dynamics of spur gears using a multi-body model. In their model, the gear body induced tooth deflection was calculated using the model proposed by Sainsot et al. [
1]. Vivet et al. [
21] gave one analytical formula to evaluate gear body induced tooth deflections. But, they did not give details how the formula was derived. Cappellini et al. [
22] developed a combined analytic-numerical contact model to simulate gear dynamics, but they did not specifically investigate gear body induced tooth deflections.
However, the analytical equations derived for
L,
M,
P, and
Q in Ref. [
1] are very complicated and hard to use. To address this problem, Sainsot and Velex [
1] numerically calculated the values of
L,
M,
P and
Q based on the plane strain assumption over a realistic range of values (
h between 1.4 and 7,
θf between 0.01 and 0.12) and fitted these values using polynomial functions of the form as shown in Equation (2)
where the constants
Ai,
Bi,
Ci,
Di,
Ei, and
Fi are given in
Table 2.
For a pair of spur gears with a contact ratio between one and two, single-tooth-pair meshing and double-tooth-pair meshing take place alternatively as shown in
Figure 2. The formula reported in Ref. [
1] can only be used in the single-tooth-pair meshing duration but not in the double-tooth-pair meshing duration as their model only considered a single loading force on a gear. To address this shortcoming, Xie et al. [
23] introduced a correction factor to consider the load dependence in the double-tooth-pair meshing duration and to account for the errors of using Equation (1) in the double-tooth-pair meshing duration. Meanwhile, Xie et al. [
24] derived new formulas for gear body induced deflections considering cubic and parabolic stress distribution around the tooth root and also load dependence between teeth. Chen et al. [
25] extended the formula reported in Ref. [
1] by considering the effect of tooth root crack on gear body-induced tooth deflections.
Based on our literature review, the formula reported in Ref. [
1] is most widely used by researchers to calculate gear body-induced deflections. The formula reported in Ref. [
24] is an improved version of the formula reported in Ref. [
1]. However, these methods still generate a large error in calculating gear body-induced deflections. There are mainly two reasons for the errors. First, the formulas were derived considering the gear body as an elastic ring without considering the effect of teeth on body deflections. Second, the stress distribution assumed in deriving these formulas is not proper. Linear and constant stress variation around the tooth root circle was assumed in Ref. [
1] while cubic and parabolic stress distributions were assumed in Ref. [
24]. Based on our finite element analysis, these assumptions are far from real situations. Detailed descriptions will be given in this study.
In this paper, we will discuss how to develop an accurate finite element model, and analyze the effect of different parameters such as gear number of teeth, location of loads, and gear geometry on gear body-induced deflections. Then, an improved analytical method will be proposed to deliver a more accurate estimation of gear body-induced deflections. The remaining part of this paper is organized as follows.
Section 2 describes our finite element model and investigates the effect of gear parameters on gear body-induced deflections.
Section 3 presents the proposed method to derive an improved formula to evaluate gear body-induced tooth deflections.
Section 4 gives case studies to validate the accuracy of the proposed method.
Section 5 draws a conclusion.
4. Validation
Two pairs of involute spur gears are used to validate the performance of our proposed method. The parameters of the gears are given in
Table 9. The torque per unit of tooth width applied on the driving gear is 10 N.m/m.
Figure 15 gives the comparisons between the proposed method and the FEM in evaluating the meshing stiffness of two gear pairs. The 2-D finite element model described in
Section 2 is used. The red lines in
Figure 15 are obtained from the 2-D finite element model. The blue lines are generated considering two parts. Part 1 is the gear body induced deflection that is evaluated using the analytical method proposed in
Section 2. Part 2 is the deflection from gear tooth (assuming rigid gear body) that is also obtained from the 2-D finite element model. We use the finite element model to evaluate the tooth induced deflections because our focus is to show the performance of the analytical method for evaluating gear body induced deflections.
From the comparison, we can see that the difference between the two methods in the single-tooth-pair meshing duration is very small. The difference in the double-tooth-pair meshing duration is larger but it is still below 10%. We think a 10% difference is still acceptable. A big reason for the difference may because we used the average value of the affiliated body stiffness in the meshing stiffness calculation. Actually, it is time-varying. A time-varying function may give results that are more accurate in gear meshing stiffness evaluation.
Our study also has limitations in that we did not give affiliated body stiffness values when gears are rotating fast with transmission errors. This will be part of our future work. This is also a challenging work. After this challenge is addressed, we will further investigate gear dynamics such as transmission errors in the future.