Next Article in Journal
Effect of Plate Curvature on Heat Source Distribution in Induction Line Heating for Plate Forming
Previous Article in Journal
Water-Covered Roof Versus Inverted Flat Roof on the Mediterranean Coast: A Comparative Study of Thermal and Energy Behavior
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Detailed Investigation of Gear Body-Induced Tooth Deflections and Development of an Improved Analytical Solution

1
School of Mechatronics Engineering, Harbin Institute of Technology, Harbin 150001, China
2
National and Local Joint Engineering Research Center of Reliability Analysis and Testing for Mechanical and Electrical Products, Zhejiang Sci-Tech University, Hangzhou 310018, China
3
Department of Mechanical Engineering, University of Alberta, Edmonton, AB T6G 1H9, Canada
4
Department of Mechanical Engineering, University of Manitoba, Winnipeg, MB R3T 5V6, Canada
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(7), 2292; https://doi.org/10.3390/app10072292
Submission received: 8 February 2020 / Revised: 17 March 2020 / Accepted: 19 March 2020 / Published: 27 March 2020
(This article belongs to the Section Mechanical Engineering)

Abstract

:
Many researchers have developed analytical methods to evaluate gear meshing stiffness. Some of these methods ignored the effect of the gear body while others used a simplified model to consider its effect. Until now, a detailed investigation of gear body-induced tooth deflections has been rare, especially for the double-tooth-pair meshing period. In this study, we present a detailed investigation of gear body-induced tooth deflections. To be specific, we will discuss how to accurately evaluate gear body-induced tooth deflections using the finite element analysis, and what are the effect of parameters such as loading and gear parameters on gear body-induced tooth deflections. Then, an improved solution is developed for evaluating the body-induced tooth deflection. In the single-tooth-pair meshing period, the improved formula is developed based on a popular formula proposed by Sainsot and Velex. This is achieved by optimizing the coefficients used in their formula to make the formula more accurate to evaluate gear body-induced tooth deflections. Meanwhile, we introduce a new term called affiliated body stiffness to evaluate the body-induced tooth deflections in the double-tooth-pair meshing period. The improved method can give higher accuracy in evaluating gear body-induced tooth deflections of spur gears with a pressure angle of 20°.

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.
δ = 1 E F b cos 2 α [ L ( u S f ) 2 + M ( u S f ) + P × ( 1 + Q tan 2 α ) ]
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)
X i ( h , θ f ) = A i / θ f 2 + B i h 2 + C i h / θ f + D i / θ f + E i h + F i
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.

2. Finite Element Analysis of Gear Body-Induced Tooth Deflections

In this section, we will investigate the effect of important factors such as loading location and gear geometry such as the number of teeth and gear bore size on gear body-induced tooth deflections. A finite element model is developed to conduct the analysis. The knowledge generated by this analysis can also be used to develop an analytical or numerical solution to evaluate gear body-induced tooth deflections.

2.1. Comparisons between Rigid Teeth and Elastic Teeth Finite Element Models

In [1,24], finite element models were developed for gears with the assumption of the elastic body and rigid teeth to evaluate gear body-induced deflections. In these models, a contact force F was applied along the action line and the tooth deflection δ was considered as the deflection caused by gear body along the action line. The rigid teeth in these finite element models generate large errors because these rigid teeth increase gear body stiffness. In our analysis, we will use elastic teeth and body to demonstrate the difference.
Through our analysis, we find there are two ways to accurately obtain gear body-induced deflections through FEM. We will describe them one by one in the following two paragraphs.
Method 1: Build an elastic gear (elastic teeth and body) using FEM and apply a contact force F on a tooth (see Figure 1). Extract acting forces generated by the force F on the tooth root (see the acting forces in Figure 3). Then, apply the same acting forces on the root of the gear rather than apply the force F on a rigid tooth, and record gear tooth deflection δb along the action line. In this way, δb is the gear body-induced tooth deflection (the deflection induced by the tooth is not included) corresponding to the contact force F.
Method 2: Build an elastic gear (elastic teeth and body) using FEM and apply a contact force F on a tooth (see Figure 1). Record the tooth deflection δ along the action line that is the deflection induced by both gear teeth and gear body. Then, in the finite element model, change the gear body to be rigid. Apply the same contact force F on the same location of the tooth and record the tooth deflection δt along the action line that is the deflection caused by the gear tooth. The gear body-induced tooth deflection corresponding to the contact force F can be calculated as δb = δδt.
A finite element model reported in Ref. [11] is modified to test the performance of our two methods and the method used in Refs. [1,24]. The gears are modeled using the element type SOLID185 and the element shape of mapped hexahedral. Six teeth are refined when the middle one or two teeth are in meshing. In our model, meshing nodes of the pinion and gear are defined theoretically and then coupled. Gear parameters are given in Table 3 and Figure 3 shows the finite element model (a two-dimensional plane strain model) developed in Ansys.
Figure 4 gives a comparison of three methods: the method reported in [1,24] and our methods 1 and 2. A constant force F (contact force per unit of tooth width F/b = 10 N/m) is applied on a single gear tooth from the tooth tip to the tooth root. The force F mimics the meshing force between a pair of gears with the transmission ration 1:1. The x-axis is the gear angular displacement from tooth tip to tooth root. The tooth tip corresponds to the angle 0 at the very left of the x-axis while the tooth root corresponds to the angle at the very right of the x-axis. The y-axis is gear-body induced deflection along the action line. From the comparison, we can see that the rigid teeth in the finite element model will lead to less tooth deflection. The rigid teeth make the gear more rigid, which generates errors. Our two methods give almost the same result and they are more accurate than the model reported in [1,24].
Figure 5 gives deformation and stress information (obtained from our FEM model) along a tooth root when a meshing force (10 N/m) is applied on the tooth tip. In Figure 5, the x-axis represents the range of tooth angle on the root circle with the angle 0 corresponding to the tooth central line. The very left and right angles in the x-axis of Figure 5 correspond to tooth fillet. When the teeth are modeled as rigid, radial deformation changes linearly along the circumferential direction of the tooth root, and the circumferential deformation is almost constant. When the teeth are elastic, near a sinusoidal radial deformation is observed along the circumferential direction of the tooth root, and the circumferential deformation is almost a constant except for linear variation near the tooth fillet. In the middle part, the normal stress in the radial direction varies linearly and the shearing stress along the circumferential direction is almost constant, but, large stress variation occurs nears tooth fillets. In addition, the change trend is quite different between a rigid tooth and an elastic tooth. The big differences in the deformation and stress demonstrate that it is not proper to use rigid teeth in FEM to obtain gear body-induced deformations.

2.2. Effect of Loading and Gear Parameters on Tooth Root Stress

Linear and constant stress variation around tooth root circle was assumed in Ref. [1] while cubic and parabolic stress distribution was used in Ref. [24]. These assumptions may not be accurate. This sub-section will reveal more realistic stress distribution around tooth root circle and evaluate the effect of loading location and gear parameters on stress distribution.

2.2.1. Effect of Loading Location

A force F (contact force per unit of tooth width F/b = 10 N/m) is applied on four locations of a tooth along the tooth action line (see Figure 6) to test the effect of loading locations on tooth root stress. Locations p1 and p2 lie in the double-meshing period while locations p3 and p4 lie in the single-meshing period. These four locations are specified using the values given in Table 4.
Figure 7 presents normal stress distribution of the four loading locations while. Figure 8 gives shearing stress distribution. As the loading location moves from tooth tip to tooth root, the normal stress and shearing stress both decrease. The normal stress distribution is like a sinusoidal curve while the shearing stress distribution looks like a mountain valley. The tooth root stress distribution is not linear or cubic or parabolic.

2.2.2. Effect of Number of Teeth on Stress Distribution

This section investigates the effect of number of teeth on gear tooth root stress distribution. Gear parameters are given in Table 3 except we test 3 different teeth. A force F (contact force per unit of tooth width F/b = 10 N/m) is applied on the tooth tip.
With the increase of the number of teeth, the tooth root normal stress and shearing stress decrease slightly, and the circumferential angle range covered by a single tooth will become smaller, as shown in Figure 9.

2.2.3. Effect of Gear Bore Size on Tooth Root Stress Distribution

In this section, we only change the gear bore size to achieve different h values: 1.5, 2.1, and 3.5. A larger h value corresponds to a smaller gear bore size. From our analysis results, we find the effect of gear bore size on tooth root stress distribution is very small, as shown in Figure 10.
Similarly, we also analyzed the effect of gear modulus on tooth root stress distribution. We find its effect is negligible.
Overall, the gear tooth root stress is largely affected by the loading location, and slightly affected by the number of teeth and gear bore size. The effect of gear modulus on tooth root stress is negligible.

3. An Improved Formula to Evaluate Gear-Body Induced Deflections

In this section, we develop an improved formula to evaluate the gear body-induced tooth deflections. For the single-tooth-pair meshing duration, Equations (1) and (2) are still used following Ref. [1] but we will update the coefficients Ai, Bi, Ci, Di, Ei, and Fi of Equation (2) by optimizing them according to our finite element results. For the double-tooth-pair meshing duration, Equation (1) is not proper to use directly. We will give a new formula.

3.1. Single-Tooth-Pair Meshing Duration

This sub-section aims to optimize the coefficients Ai, Bi, Ci, Di, Ei, and Fi of Equation (2) in order to get higher accuracy when we evaluate gear body induced tooth deflection in the single-tooth-pair meshing duration. The objective function of this optimization problem is to minimize the difference between the results generated by Equation (1) and a well-established finite element model.
A finite element model reported in Ref. [11] is modified to evaluate the gear body-induced tooth deflections of 72 standard involute spur gears. The gears are modeled using the element type SOLID185 and the element shape of mapped hexahedral. Six teeth are refined when the middle one or two teeth are in meshing. In our model, meshing nodes of the pinion and gear are coupled. ANSYS is used to implement the finite element analysis. These gears (72 = 9 × 8) are a combination of of 9 teeth and 8 different h values, as shown in Table 5. These gears cover a realistic range of θf and h values: θf between 0.03 and 0.15 (radian), and h between 2.1 and 7 (remember that L, M, P and Q are functions of θf and h). All the gears have a module of 3.175 and a pressure angle of α = 20 ° . Gear body-induced tooth deflections of these 72 gears are evaluated using the finite element model to optimize the coefficients in Equation (2).
Based on our finite element analysis, gear body-induced tooth deflections are independent of the module of gears. This conclusion can also be confirmed by the analytical formulas reported in Ref. [1]. We can see that the module of gears is not a variable in Equation (1). Therefore, our equations derived in this paper can be used for a spur gear with any modulus. However, since we only used the pressure angle of α = 20 ° in our finite element analysis, the optimized coefficients may not give accurate results for other pressure angles of a spur gear.
In the reported finite element models [1,23], gear teeth were assumed as rigid in evaluating the gear body-induced tooth deflections. However, this assumption causes errors of about 10% in the gear body-induced tooth deflections based on our analysis. This assumption is released in our finite element analysis. When a force F is applied on the tooth profile (see Figure 3), we first extract acting forces generated by the force F on the tooth root. Then, we apply the same acting forces on the root of the gear and record gear tooth deflection δ along the force action line. In this way, δ is the gear body-induced tooth deflection and it is obtained from an elastic gear (elastic gear teeth and body).
A two-dimensional (2-D) plane strain model is used in our finite element analysis with element type PLANE183 as it can give similar accuracy in meshing stiffness evaluation as a 3-D finite element model. We use the 2-D model because it requires much lower computational cost. To demonstrate that it is acceptable to use the plane strain conditions, we did tests on gears with the same parameters (number of teeth: 19, gear module: 3.175, bore radius: 17.5 mm, meshing force per unit of tooth width F/b = 10 N/m, transmission ratio 1:1, pressure angle α = 20 ° ) except for the tooth face width. Figure 11 shows the comparison results. The x-axis of Figure 11 is the gear angular displacement corresponding to the meshing force moving from the tooth tip to root. The y-axis is the gear body stiffness (some researchers call it body stiffness) per unit of tooth width which is expressed as kb = F / b δ , where F is the meshing force, b represents tooth width, and δ denotes the gear body-induced tooth deflection. In the whole process, F is applied on a single gear tooth since we focus on the single-tooth-pair meshing duration in this sub-section. The force F mimics the meshing force between a pair of gears with the transmission ration 1:1. From Figure 11, we can see that when the tooth width is above 20 mm, the error generated by the 2-D plane strain model is negligible. Even when the tooth width is 10 mm, the error is less than 10%. It is acceptable to use the 2-D plane strain model in our finite element analysis. When the tooth width is much smaller than 10 mm, we do not recommend using the tooth deflection equations derived in this paper.
Using the above described 2-D plane strain finite element model, we evaluated the gear body stiffness for 72 standard involute spur gears as given in Table 5. Then the Genetic Algorithm was used to optimize the 24 coefficients that are required in Equation (2). The objective of the optimization is to let the gear body-induced tooth deflections evaluated using the Equation (1) closest to the results obtained from the finite element analysis. Table 6 gives the updated coefficients obtained from our optimization.
A comparison is conducted between our proposed coefficients and the coefficients reported in Ref. [1]. We use these coefficients to obtain the gear body-induced tooth deflections analytically, respectively. The y-axis of Figure 12 gives the maximum relative difference of the gear body stiffness per unit of tooth width between the analytical method and the 2-D plane strain finite element model. The maximum error generated using our proposed coefficients is below 10% whereas the coefficients reported in Ref. [1] can give the error as high as 20%. In 50 out of the total 72 cases (see Table 5), the errors generated using our coefficients are lower than those in Ref. [1]. In the other 22 cases, the method with our coefficients performs a little worse (the average difference between the two methods is smaller than 1.5%). Overall, our proposed coefficients give more accurate results in evaluating the gear body-induced tooth deflections in the single-tooth-pair meshing duration.

3.2. Double-Tooth-Pair Meshing Duration

In Section 2.1, we optimized the 24 coefficients used in Equation (1) that was derived by Sainsot et al. [1]. Our coefficients give better performance in the evaluation of gear body-induced tooth deflections during the single-tooth-pair meshing period. However, this analytical formula is not proper to use directly for the double-tooth-pair meshing duration. This section aims to derive an improved formula for the double-tooth-pair meshing duration.
In the double-tooth-pair meshing duration, there are two pairs of teeth meshing simultaneously. Figure 13 is a simple illustration of gear meshing forces and their resulting defections. In this illustration, we assume gear teeth do not have tooth deflection since our focus is only the gear body deflections. In the figure, F1 and F2 represent the gear meshing forces of the first and the second tooth pairs, respectively. And P1 and P2 denote the gear meshing points of the first and second gear tooth pairs, respectively. Imagine if the first tooth pair starts to mesh first and, correspondingly, a deflection ∆11 is generated in the mesh point P1 and a deflection ∆21 is generated in the meshing point P2. Later, when the second tooth pair starts to mesh (the first tooth pair is still in meshing), further deflections ∆12 and ∆22 are generated in the meshing points P1 and P2, respectively. Therefore, the total deflection is 11 + 12 (δ1 = ∆11 + ∆12) at the meshing point P1 and ∆21 + ∆22 (δ2 = ∆21 + ∆22) at the meshing point P2. The values of ∆11 and ∆22 can be obtained from the analytical equation derived for the single-tooth-pair meshing duration. Then, our remaining problem is to get the values of ∆12 and ∆21.
As we described in the previous paragraph, if we apply the force F1 followed by F2, the energy stored in the system can be obtained as follows:
W = 1 2 F 1 Δ 11 + 1 2 F 2 Δ 22 + F 1 Δ 12
Similarly, if we apply the force F2 followed by F1, the energy stored in the system can be expressed as Equation (4) based on one assumption that the deflection is proportional to the load. This assumption is generally valid for elastic materials with small deflection.
W = 1 2 F 1 Δ 11 + 1 2 F 2 Δ 22 + F 2 Δ 21
From Equations (3) and (4), we can conclude:
F 1 Δ 12 = F 2 Δ 21
Further, we can get:
Δ 21 = F 1 F 2 Δ 12 = λ 1 λ Δ 12
where λ denotes the load sharing ratio of a pair of spur gears with the expression λ = F 1 F 1 + F 2 . Equation (6) gives the relationship between Δ 12 and Δ 21 (see Figure 13 for details about these two deflections).
For better illustration, Δ 21 can be expressed as follows
Δ 21 = F 1 F 1 / Δ 21 = F 1 k 21 A
where k 21 A is the ratio of F1 and Δ 21 . We call k 21 A affiliated body stiffness in this paper. Similarly, we can define the other affiliated body stiffness k 12 A with the expression k 12 A = F 2 Δ 12 . From Equation (5), we can find that k 21 A is equal to k 12 A . From Equation (7), we can see that the deflection Δ 21 can be obtained if the affiliated body stiffness is known. Then, from Equation (6), the deflection of ∆12 is available if Δ 21 and the load sharing ratio λ are known.
We used the above-mentioned 2-D plane strain finite element model to evaluate the affiliated body stiffness for the 72 gears with the parameters given in Table 5. Figure 14 gives an example to illustrate the affiliated body stiffness of two gear pairs. Parameters of the two gear pairs are given in Table 7. The affiliated body stiffness is not constant while it is almost linearly proportional to the gear rotation angle (we know gear meshing point moves as the gear rotates). But, we find the deviation of the affiliated body stiffness values from their average value is small. The first gear pair has a maximum deviation of −3.18% while the second gear pair owns a maximum deviation of 3.23% (see Figure 14). For all other gears, the maximum deviation is around or below 3%. Therefore, we believe it is reasonable to use the average value of the affiliated body stiffness values of gear to represent this gear’s affiliated body stiffness. Some errors will be generated, but in most cases, the error is around or below 5%. Table 8 gives the average of the affiliated body stiffness values per unit of tooth width for 72 gears.
Notably, the affiliated body stiffness is independent of one gear’s module, the same as the gear body stiffness.

3.3. Meshing Stiffness of a Pair of Gears

According to the potential energy method [2], the meshing stiffness of a pair of spur gears can be evaluated as follows:
k l = { F 1 δ 1 + δ 1 + F 2 δ 2 + δ 2 , double   tooth   pairs   in   meshing F 1 δ 1 + δ 1 , sin gle   tooth   pair   in   meshing
where F1 and F2 are gear meshing forces, and δ 1 , δ 1 , δ 2 , and δ 2 are tooth deflections along the action line corresponding to the driving gear tooth of the first tooth pair, the driven gear tooth of the first tooth pair, the driving gear tooth of the second tooth pair, and the driven gear tooth of the second tooth pair, respectively.
Each of the tooth deflections ( δ 1 , δ 1 , δ 2 , and δ 2 ) can be considered as the summation of two parts: the deflection from gear tooth (assuming rigid gear body) and the deflection caused by gear body. The deflection from gear tooth (assuming rigid gear body) can be evaluated using the potential energy method [3] or the FEM [11]. The deflections from the gear body can be evaluated using the method proposed in this paper. In the single-tooth-pair meshing duration, gear body induced tooth deflection can be evaluated using Equation (1) with the parameters given in Table 6. In the double-tooth-pair meshing duration, the gear body induced tooth deflection needs to be considered as the summation of two parts. Part 1 can be evaluated using Equation (1) with the parameters given in Table 6. Part 2 needs to be evaluated using Equations (6) and (7) with the affiliated body stiffness given Table 8. In this process, the load sharing ratio of a pair of gears needs to be known.

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.

5. Conclusions

This study gives a detailed investigation of gear body-induced tooth deflections and presents an improved analytical solution to evaluate gear body-induced tooth deflections. First, we find it is not proper to assume rigid teeth in finite element models to evaluate gear body-induced tooth deflections. The assumption generates a large error. We also proposed two ways to accurately evaluate gear body-induced tooth deflections using the finite element method. Then, through finite element analysis, we find that loading location has a large effect on gear body induced tooth root stress distribution, the effect from the number of teeth and the bore size is small, and the effect from gear modulus is negligible. Gear tooth root stress distribution is not linear or cubic or parabolic. The gear tooth normal stress caused by a gear body is kind of a sinusoidal curve while the shearing stress looks like a mountain valley. Last but not least, an improved analytical solution is developed to accurately evaluate gear body-induced deflections. In the single-tooth-pair meshing duration, we optimized the coefficients of the formula reported in Ref. [1]. In the double-tooth-pair meshing duration, we introduced the concept of affiliated body stiffness and provided a table to obtain the affiliated body stiffness. The improved analytical solution is more accurate than existing analytical methods to evaluate the gear body-induced tooth deflections, which can also be used to incorporate the effect of gear body in gear meshing stiffness evaluation. However, our study is limited to spur gears with a pressure angle of 20°. In the future, we will investigate gear body-induced tooth deflections for spur gear with other pressure angles and other types of gears such as helical gear and bevel gear.

Author Contributions

H.Z. and X.L. together proposed this methodology. H.Z. did finite element analysis and calculations. X.L. wrote most of the manuscript. J.P. helped to develop the finite element models. M.R. helped to optimize formula paramters. Z.L. helped to write part of this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Harbin Institute of Technology Key Laboratory Open Fund, China [Grant No. 30620170208], the Natural Science and Engineering Research Council of Canada [Grant No. RGPIN-2019-05361], the International S&T Cooperation Program of China, China [Grant No.2015DFA71400], and the University of Manitoba University Research Grants Program.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sainsot, P.; Velex, P.; Duverger, O. Contribution of gear body to tooth deflections—A new bidimensional analytical formula. J. Mech. Des. 2004, 126, 748–752. [Google Scholar] [CrossRef]
  2. Liang, X.; Zhang, H.; Zuo, M.J.; Qin, Y. Three new models for evaluation of standard involute spur gear mesh stiffness. Mech. Syst. Signal Process. 2018, 101, 424–434. [Google Scholar] [CrossRef]
  3. Liang, X.; Zuo, M.J.; Feng, Z. Dynamic modeling of gearbox faults: A review. Mech. Syst. Signal Process. 2018, 98, 852–876. [Google Scholar] [CrossRef]
  4. Leonard, J.W.; Nath, J.H. Comparison of finite element and lumped parameter methods for oceanic cables. Eng. Struct. 1981, 3, 153–167. [Google Scholar] [CrossRef]
  5. Liang, X.; Zuo, M.J.; Patel, T.H. Evaluating the time-varying mesh stiffness of a planetary gear set using the potential energy method. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2014, 228, 535–547. [Google Scholar] [CrossRef]
  6. Yang, D.C.H.; Lin, J.Y. Hertzian damping, tooth friction and bending elasticity in gear impact dynamics. J. Mech. Des. 1987, 109, 189–196. [Google Scholar] [CrossRef]
  7. Liang, X.; Zuo, M.J.; Pandey, M. Analytically evaluating the influence of crack on the mesh stiffness of a planetary gear set. Mech. Mach. Theory 2014, 76, 20–38. [Google Scholar] [CrossRef]
  8. Chen, Z.; Shao, Y. Mesh stiffness calculation of a spur gear pair with tooth profile modification and tooth root crack. Mech. Mach. Theory 2013, 62, 63–74. [Google Scholar] [CrossRef]
  9. Ma, H.; Song, R.; Pang, X.; Wen, B. Time-varying mesh stiffness calculation of cracked spur gears. Eng. Fail. Anal. 2014, 44, 179–194. [Google Scholar] [CrossRef]
  10. Pandya, Y.; Parey, A. Failure path based modified gear mesh stiffness for spur gear pair with tooth root crack. Eng. Fail. Anal. 2013, 27, 286–296. [Google Scholar] [CrossRef]
  11. Liang, X.; Zhang, H.; Liu, L.; Zuo, M.J. The influence of tooth pitting on the mesh stiffness of a pair of external spur gears. Mech. Mach. Theory 2016, 106, 1–15. [Google Scholar] [CrossRef]
  12. Weber, C. The Deformation of Loaded Gears and the Effect on Their Load Carrying Capacity; Technical Report No.3, Sponsored Research; Department of Scientific and Industrial Research: Delhi, India, 1951. [Google Scholar]
  13. Attia, A.Y. Deflection of spur gear teeth cut in thin rims. J. Eng. Ind. 1964, 86, 333–341. [Google Scholar] [CrossRef]
  14. Cornell, R.W. Compliance and stress sensitivity of spur gear teeth. J. Mech. Des. 1981, 103, 447–459. [Google Scholar] [CrossRef] [Green Version]
  15. Muskhelishvili, N. Some Basic Problems of the Mathematical Theory of Elasticity; P. Noordhoff Limited: Groningen, The Netherlands, 1953; pp. 230–235. [Google Scholar]
  16. Shao, Y.; Chen, Z. Dynamic features of planetary gear set with tooth plastic inclination deformation due to tooth root crack. Nonlinear Dyn. 2013, 74, 1253–1266. [Google Scholar] [CrossRef]
  17. Wang, Q.; Zhang, Y. A model for analyzing stiffness and stress in a helical gear pair with tooth profile errors. J. Vib. Control 2015, 23, 272–289. [Google Scholar] [CrossRef]
  18. Wan, Z.; Cao, H.; Zi, Y.; He, W.; Chen, Y. Mesh stiffness calculation using an accumulated integral potential energy method and dynamic analysis of helical gears. Mech. Mach. Theory 2015, 92, 447–463. [Google Scholar] [CrossRef]
  19. Ma, H.; Zeng, J.; Feng, R.; Pang, X.; Wang, Q.; Wen, B. Review on dynamics of cracked gear systems. Eng. Fail. Anal. 2015, 55, 224–245. [Google Scholar] [CrossRef]
  20. Cirelli, M.; Valentini, P.P.; Pennestrì, E. A study of the non-linear dynamic response of spur gear using a multibody contact based model with flexible teeth. J. Sound Vib. 2019, 445, 148–167. [Google Scholar] [CrossRef]
  21. Vivet, M.; Mundo, D.; Tamarozzi, T.; Desmet, W. An analytical model for accurate and numerically efficient tooth contact analysis under load, applied to face-milled spiral bevel gears. Mech. Mach. Theory 2018, 130, 137–156. [Google Scholar] [CrossRef]
  22. Cappellini, N.; Tamarozzi, T.; Blockmans, B.; Fiszer, J.; Cosco, F.; Desmet, W. Semi-analytic contact technique in a non-linear parametric model order reduction method for gear simulations. Meccanica 2018, 53, 49–75. [Google Scholar] [CrossRef]
  23. Xie, C.; Hua, L.; Lan, J.; Han, X.; Wan, X.; Xiong, X. Improved analytical models for mesh stiffness and load sharing ratio of spur gears considering structure coupling effect. Mech. Syst. Signal Process. 2018, 111, 331–347. [Google Scholar] [CrossRef]
  24. Xie, C.; Hua, L.; Han, X.; Lan, J.; Wan, X.; Xiong, X. Analytical formulas for gear body-induced tooth deflections of spur gears considering structure coupling effect. Int. J. Mech. Sci. 2018, 148, 174–190. [Google Scholar] [CrossRef]
  25. Chen, Z.; Zhang, J.; Zhai, W.; Wang, Y.; Liu, J. Improved analytical methods for calculation of gear tooth fillet-foundation stiffness with tooth root crack. Eng. Fail. Anal. 2017, 82, 72–81. [Google Scholar] [CrossRef]
  26. Ma, H.; Pang, X.; Feng, R.; Wen, B. Evaluation of optimum profile modification curves of profile shifted spur gears based on vibration responses. Mech. Syst. Signal Process. 2016, 70–71, 1131–1149. [Google Scholar] [CrossRef]
Figure 1. An elastic tooth subjected to a force F.
Figure 1. An elastic tooth subjected to a force F.
Applsci 10 02292 g001
Figure 2. Gear meshing scenarios.
Figure 2. Gear meshing scenarios.
Applsci 10 02292 g002
Figure 3. Illustration of acting forces on the tooth root.
Figure 3. Illustration of acting forces on the tooth root.
Applsci 10 02292 g003
Figure 4. Comparison of three methods for calculating gear body-induced tooth deflection.
Figure 4. Comparison of three methods for calculating gear body-induced tooth deflection.
Applsci 10 02292 g004
Figure 5. Deformation and stress along tooth root.
Figure 5. Deformation and stress along tooth root.
Applsci 10 02292 g005
Figure 6. Forces applied on different locations of a tooth.
Figure 6. Forces applied on different locations of a tooth.
Applsci 10 02292 g006
Figure 7. Effect of loading locations on tooth root normal stress.
Figure 7. Effect of loading locations on tooth root normal stress.
Applsci 10 02292 g007
Figure 8. Effect of loading locations on tooth root shearing stress.
Figure 8. Effect of loading locations on tooth root shearing stress.
Applsci 10 02292 g008
Figure 9. Effect of number of teeth on tooth root stress distribution.
Figure 9. Effect of number of teeth on tooth root stress distribution.
Applsci 10 02292 g009
Figure 10. Effect of gear bore size on tooth root stress distribution.
Figure 10. Effect of gear bore size on tooth root stress distribution.
Applsci 10 02292 g010
Figure 11. Comparison between 2-D and 3-D models.
Figure 11. Comparison between 2-D and 3-D models.
Applsci 10 02292 g011
Figure 12. Performance comparison in the single-tooth-pair meshing duration.
Figure 12. Performance comparison in the single-tooth-pair meshing duration.
Applsci 10 02292 g012
Figure 13. Illustration of meshing forces and the resulting deflections.
Figure 13. Illustration of meshing forces and the resulting deflections.
Applsci 10 02292 g013
Figure 14. Affiliated body stiffness of two gear pairs.
Figure 14. Affiliated body stiffness of two gear pairs.
Applsci 10 02292 g014
Figure 15. Comparison between the proposed and the finite element method.
Figure 15. Comparison between the proposed and the finite element method.
Applsci 10 02292 g015
Table 1. L, M, P and Q values used in four different papers.
Table 1. L, M, P and Q values used in four different papers.
LMPQ
Weber [12]5.211.40.294
Attia [13] 0.32
Cornell [14]5.3061.4 (plane stress)1.5340.32
1.14 (plane strain)
Sainsot and Velex [1]L(h, θf)M(h, θf)P(h, θf)Q(h, θf)
Table 2. Coefficients of the polynomial functions [1].
Table 2. Coefficients of the polynomial functions [1].
AiBiCiDiEiFi
L(h, θf)−5.574 × 10−5−1.9986 × 10−3−2.3015 × 10−44.7702 × 10−30.02716.8045
M(h, θf)60.111 × 10−528.100 × 10−3−83.431 × 10−4−9.9256 × 10−30.16240.9086
P(h, θf)−50.952 × 10−5185.50 × 10−30.0538 × 10−453.300 × 10−30.28950.9236
Q(h, θf)−6.2042 × 10−59.0889 × 10−3−4.0964 × 10−47.8297 × 10−3−0.14720.6904
Table 3. A pair of gears with the transmission ratio 1:1.
Table 3. A pair of gears with the transmission ratio 1:1.
Number of TeethGear ModuleBore RadiusPressure Angle
193.17517.520°
Table 4. Four loading positions.
Table 4. Four loading positions.
p1p2p3p4
Circumferential angle θ(o)0.00013.16314.62527.787
Height u (mm) (see Figure 1: An elastic tooth subjected to a force F)6.49773.40183.18602.1688
Table 5. Gear parameters used in the finite element analysis.
Table 5. Gear parameters used in the finite element analysis.
Number of Teeth2030405060708090100
h2.12.83.54.24.95.66.37.0
Table 6. Updated coefficients for L, M, P and Q.
Table 6. Updated coefficients for L, M, P and Q.
AiBiCiDiEiFi
L2.78263 × 10−3−1.96615 × 10−3−1.85863 × 10−44.77020 × 10−3−3.05544 × 10−36.80474
M−2.13880 × 10−32.90648 × 10−2−8.34310 × 10−3−9.91337 × 10−30.1625560.910564
P−4.98776 × 10−40.1882522.24233 × 10−45.41061 × 10−20.2895000.955057
Q−1.72430 × 10−40.0180115−9.69840 × 10−52.74944 × 10−2−0.1431200.690752
Table 7. Parameters of two standard involute spur gear pairs.
Table 7. Parameters of two standard involute spur gear pairs.
Number of TeethModulusPressure AngleBore RadiusTransmission Ratiohθf (°)
First pair703.17520°20.411:12.12.4660
Second pair1003.17520°29.481:12.11.7262
Table 8. Affiliated body stiffness represented by the average value (unit: N/m/m).
Table 8. Affiliated body stiffness represented by the average value (unit: N/m/m).
H = 2.1H = 2.8H = 3.5H = 4.2H = 4.9H = 5.6H = 6.3H = 7.0
Z = 2012.72007.85255.45314.01503.07492.42631.96051.6154
Z = 3011.62297.60885.46534.11283.19732.54992.07661.7212
Z = 4010.81287.33815.38314.10783.22472.58992.12031.7644
Z = 5010.21827.10625.28974.07663.22282.60172.13811.7845
Z = 609.76266.91195.20174.03913.21042.60202.14481.7942
Z = 709.40086.74815.12234.00123.19402.59712.14601.7986
Z = 809.10356.60755.05083.96493.17652.58982.14431.8001
Z = 908.85386.48534.98653.93083.15882.58132.14111.7998
Z = 1008.64026.37804.92863.89913.14162.57242.13701.7985
Table 9. Physical parameters of two spur gear pairs.
Table 9. Physical parameters of two spur gear pairs.
ParametersFirst Gear Pair [7]Second Gear Pair [26]
Driving GearDriven GearDriving GearDriven Gear
Number of teeth19315272
Module (mm)3.21.75
Pressure angle20°
Young’s modulus (GPa)206.8210
Poisson’s ratio0.3
h1.49682.58542.54783.0406

Share and Cite

MDPI and ACS Style

Zhang, H.; Pan, J.; Rao, M.; Li, Z.; Liang, X. A Detailed Investigation of Gear Body-Induced Tooth Deflections and Development of an Improved Analytical Solution. Appl. Sci. 2020, 10, 2292. https://doi.org/10.3390/app10072292

AMA Style

Zhang H, Pan J, Rao M, Li Z, Liang X. A Detailed Investigation of Gear Body-Induced Tooth Deflections and Development of an Improved Analytical Solution. Applied Sciences. 2020; 10(7):2292. https://doi.org/10.3390/app10072292

Chicago/Turabian Style

Zhang, Hongsheng, Jun Pan, Meng Rao, Zhengwei Li, and Xihui Liang. 2020. "A Detailed Investigation of Gear Body-Induced Tooth Deflections and Development of an Improved Analytical Solution" Applied Sciences 10, no. 7: 2292. https://doi.org/10.3390/app10072292

APA Style

Zhang, H., Pan, J., Rao, M., Li, Z., & Liang, X. (2020). A Detailed Investigation of Gear Body-Induced Tooth Deflections and Development of an Improved Analytical Solution. Applied Sciences, 10(7), 2292. https://doi.org/10.3390/app10072292

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop