Due to machining errors, location inaccuracies, human error, and various other factors, it is challenging to avoid assembly errors during the production of spiral bevel gears (SBGs). When SBG assembly errors occur, it can cause the appearance of edge contact and may even lead to severe tooth tip chipping. In this study, we propose an improved method based on loaded tooth contact analysis (LTCA) to examine mesh characteristics, including time-varying mesh stiffness (TVMS), unloaded transmission error, and contact stress. Furthermore, we explore the effects of assembly errors and tooth tip chipping. Moreover, it is observed that assembly errors can alter the contact area of SBGs and potentially reduce the peak-to-peak value of TVMS. Additionally, the occurrence of tooth tip chipping decreases TVMS within the chipping region, lowers transmission error, and increases maximum contact stress. Notably, when assembly errors are present, the reduction in TVMS due to tooth tip chipping exceeds that of a properly assembled SBG pair.
Keywords: spiral bevel gear; mesh characteristic; assembly error; tooth tip chipping
Spiral bevel gears (SBGs) are widely used in aerospace, automobile, and maritime industries. With the progress of science and technology, and the increasing demands in industrial applications, the study of tooth contact performance and fault characteristics of SBGs is becoming increasingly important. There are a large number of scholars who research the mesh characteristics of SBGs. Ding H et al. [[
Extensive research has been conducted both domestically and internationally to investigate the impact of assembly errors on SBGs. Wang B et al. [[
The majority of current studies on gear faults focus on spur or helical gears, with only a limited number of researchers exploring gear faults in SBGs. Furthermore, tooth tip chipping (see Figure 1), a significant aspect of SBG faults, has received very little attention in the existing literature. At present, there are four main methods for studying gear faults: (a) Finite element method (FEM). Jia S et al. [[
When unloaded, SBGs have a local point contact. However, when loaded, they transition to elliptical contact. This leads to numerous challenges in modeling and analysis, particularly when considering assembly errors and tooth tip chipping. Undesirable assembly errors can result in the movement of the contact area from the center of the tooth surface to the tooth tip which will cause edge contact. An increase in edge contact may worsen the mesh characteristics of SBGs, and even lead to tooth tip chipping. Hence, conducting research on SBG considering assembly errors and tooth tip chipping is deemed highly essential. Despite the existence of research focusing on assembly errors, there remains a scarcity of scholars who have explored the combined effects of assembly errors and tooth tip chipping, notably in SBG. Hence, this paper primarily focuses on two aspects: (
The structure of this paper is outlined as follows: Section 2 develops an enhanced method which is based on LTCA and incorporates considerations for assembly errors and tooth tip chipping. Section 3 validates the proposed method with publications, and analyzes the effects of assembly errors and tooth tip chipping. Finally, Section 4 provides a summary of the findings presented in this paper.
In this section, we will introduce an enhanced method which is a LTCA of SBGs considering assembly error and tooth tip chipping.
The method proposed in this paper is based on loaded tooth contact analysis of spiral bevel gear pairs. The detailed processes of TCA and LTCA are stated in the paper published by Li et al. [[
The TCA is divided into two parts: (
The LTCA can be generalized into four parts: (
Assembly errors do not directly impact the tooth, therefore, here we will only introduce the method for bringing four assembly errors into TCA.
As shown in Figure 2, S
(
(
where ΔAP is the axial assembly error of pinion; ΔAG is the axial assembly error of gear; ΔE is the offset error; ΔΣ is the shaft angle error; and Φ
Then, the tooth surface equation and the normal vector of pinion and gear in the meshing coordinate system considering assembly errors can be calculated as:
(
where,
Given the rotor angle of the pinion (Φ
(
where v
The unloaded transmission error can be calculated as follows:
(
where z
We did not consider lubrication in our analysis because the presented method was a quasi-static approach. According to LTCA, the TVMS can be calculated by iterating the equation of compatibility which takes normal forces of contact points and the angle of pinion as a solution. The equation governing compatibility at an arbitrarily chosen meshing moment is as follows:
(
where λ
Following this, we can obtain the contact area (see Figure 3), which consists of potential contact points for which the distributed normal force exceeds zero for various meshing moments of SBG. Figure 3a is the rotating projection of a healthy pinion without assembly error and tooth tip chipping, and Figure 3b depicts the contact area on the rotating projection of the pinion considering the gear axial error ΔAG = 0.5 mm. Observations reveal that when the gear axial error ΔAG is 0.5 mm, the contact area shifts towards the tooth tip and heel.
As shown in Figure 4, a three-tooth model is established in MATLAB (R2021b) software. The proposed method needs only to get the global compliances obtained by the traversal loading method. Therefore, we constrain all degrees of nodes at the bottom surface of the pinion and gear when simulating the working conditions of SBG pairs. The elements used are 8-node hexahedral elements. Moreover, it is worth noting that the grids in the contact area are not fine, as satisfactory global compliances can be achieved through spatial interpolation and contact compliances can be estimated using empirical formulas instead of direct contact calculations.
TVMS of SBGs can be expressed as:
(
where, R is the mean cone distance; the δ
Hertz contact theory plays a significant role in making out the local stress and strain between two subjects contacting due to pressure. This paper utilizes Hertz contact theory to calculate the contact stress of SBG pairs. To apply Hertz contact theory, it is essential to obtain the curvature of each potential contact point. There is a three-point method to obtain the curvature. According to three spatial points P
(
where the (x, y, z) is the coordinate of the center of the circle across the three points, and κ
When having the coordinate of the center of a circle, the curve normal vector n can be computed by the positional relationship between three points and the center point. The surface normal curvature κ
(
According to Hertz contact theory, the contact stress at potential contact point can be given as:
(
where F
In the absence of load, there is local point contact. The contact points can form a line known as the initial contact trajectory during a mesh cycle. To simplify calculations, this trajectory is represented by dispersed points which are called initial contact points. When under load, SBGs transition to elliptical contact. To simplify calculations, we utilize a line which is long enough but cannot exceed the tooth profile along the long axis of the elliptical contact region instead of elliptical contact surfaces. Additionally, the line is represented by dispersed points, referred to as potential contact points, for further simplification.
The proposed method can accurately incorporate and take into consideration the actual shapes of the tooth tip chipping. Initially, it needs to extract the shape (see Figure 5a) of tooth tip chipping from the gear affected by the chipping fault. When tooth tip chipping occurs, there is a reduction in the clearances between the pinion and gear within the chipping area (see Figure 5b). For simulating the occurrence of tooth tip chipping, this paper adjusts the clearances (see Figure 5c) between the potential contact points within the chipping area to represent the chip. The increased values in clearance reflect the depth of tooth tip chipping, whereas the extent of the chipped area is presented by the potential contact points for which clearances are altered. Subsequently, the equation governing compatibility is solved anew using the updated clearances, enabling the analysis of the mesh characteristics of the gear under tooth tip chipping fault.
The method presented in this paper is outlined in a simplified flowchart depicted in Figure 6. Firstly, based on the gear and machine tool parameters, the assembly errors are incorporated into the TCA to acquire the data for LTCA and to calculate the unloaded transmission error. Subsequently, LTCA is executed. Prior to assessing the compatibility, tooth tip chipping is indicated by noting the changes in clearance within the tooth tip chipping area. Following this, the TVMS can be computed using the compatibility equation. Finally, contact stress is determined using the Hertz contact formula, leveraging the curvature and normal force at potential contact points obtained from LTCA.
This section will validate the LTCA considering assembly error with literature and analyze the effect of assembly error and tooth tip chipping for TVMS, transmission error, and contact stress of SBGs.
For validating the correctness of the proposed method in this paper, geometry parameters and machining parameters of the example SBGs [[
On the rotating projection of the gear tooth surface, the contact trajectories of the example SBGs obtained by the presented method under different assembly errors are compared in Figure 7. In addition, the comparisons of contact trajectories due to assembly errors in Han's paper [[
The effects of assembly error are investigated in this section. There are three aspects to illustrate illustrating the effects of assembly error. They are the TVMS, unloaded transmission error, and contact stress. There are four values for ΔAP, ΔAG and ΔE, which are −0.5 mm, −0.2 mm, 0.2 mm, 0.5 mm, and four values for ΔΣ, which are −0.4°, −0.2°, 0.2°, 0.4°. Following this, the two chip areas are compared and the effects of tooth tip chipping discussed further.
The TVMS of SBGs considering four types of assembly errors are shown in Figure 9. From this, the following four main points can be concluded: (
From the variations and trends above, it can be seen that all exhibit reductions in TVMS and peak-to-peak values of TVMS. This discrepancy arises from the contact of the tooth root, which can lead to extreme vibrations, loud noises, significant contact forces, and even crack formation, severely impacting the lifespan of SBGs. In practical applications, many SBGs undergo modifications to the tooth root to prevent contact. Consequently, the presented method solely considers the contact of the tooth surface. If there is contact in the tooth root area, the method will not distribute force to potential points in the tooth root area, resulting in a reduction in TVMS due to the decrease in bearing area.
The unloaded transmission errors considering different assembly errors are displayed in Figure 10. The range from the x-coordinate of the lowest point on the left side to zero is called the left meshing location and the range on the right side is called the right meshing location. There is a common phenomenon that the unloaded transmission errors will gradually equalize at the left meshing location and gradually reach their lowest point at the right meshing location. Due to the complexity of the SBG's tooth surface, assembly errors will lead to changes in the contact area, resulting in unpredictable variations in unloaded transmission error. As of now, no specific rules can be summarized regarding these variations. Assembly errors can induce changes in unloaded transmission errors. Favorable changes, such as ΔAP = −0.5 mm and so forth, can increase the lowest point of unloaded transmission error.
The contact stresses of SBGs considering different assembly errors are shown in Figure 11. It is evident that the contact area, which consists of potential contact points wherein the contact stress exceeds zero, and the initial contact trajectory of ΔAP and ΔΣ move from heel and tip to toe and root. On the contrary, those of ΔAG and ΔE move from the toe and root to the heel and tip. In addition, with the movement of the contact area, the edge contact and the increase of tip stress appear gradually. This phenomenon alters the load distribution across the tooth surface and elevates the risk of tooth tip chipping. Moreover, it diminishes the service life of SBGs.
As Section 3.2 found, some assembly errors can result in edge contact and an increase in tip stress, ultimately leading to tooth tip chipping. Consequently, this section will discuss the effects of tooth tip chipping from three key aspects: TVMS, unloaded transmission error, and contact stress. Each assembly error selects a value characterized by edge contact and high tooth tip stress to introduce an identical area of chipping. They conclude ΔAP = −0.5 mm, ΔAG = 0.5 mm, ΔE = 0.5 mm, and ΔΣ = −0.4°.
In general, assembly errors will decrease TVMS in the influence range of tooth tip chipping (see Figure 12). A, B, D, E are the moments which have the maximal reduction of TVMS at ΔAP, ΔAG, ΔE and ΔΣ, respectively. The reduction of A (see Figure 12a), B (see Figure 12b), D (see Figure 12c), and E (see Figure 12d), is 3.28%, 3.92%, 4.91, 3.88%, respectively. It is evident that the effect of ΔE is the most significant when compared to ΔAP and ΔAG at the same error magnitude. ΔAP exhibits the smallest effect, although its influence range on tooth tip chipping is the widest. An intriguing observation is that the TVMS of ΔAP and ΔE may exhibit increases at certain points in time. For example, the TVMS of moment C (see Figure 12c) increased 2.31%. That is because tooth tip chipping causes a change in the smallest clearance (see Figure 13a). Due to tooth tip chipping, the clearances at the chipping region increase, and the smallest clearance changes from M (0 mm) in tooth 1 to P (0.0113 mm) in tooth 1. It is worth noting that the clearance at P is nearly equal to N located in tooth 2. Tooth 1 initially bears a load of 12.45 kN (see Figure 13b), whereas tooth 2 carries 14.76 kN before tooth tip chipping occurs. However, following tooth tip chipping, there is a shift in load distribution, with tooth 2 now bearing a heavier load of 18.48 kN, surpassing by far tooth 1's load of 8.64 kN.
As shown in Figure 14, The unloaded transmission error experiences varying degrees of reduction as a result of tooth tip chipping. The reduction of ΔE = 0.5 mm (see in Figure 14c) is 131.3402 μrad which is the most significant and the reduction of ΔΣ = −0.4° (see in Figure 14d) is 62.1242 μrad which is the smallest relative to the other three assembly errors. The effect of ΔAP = −0.5 mm (see in Figure 14a) is similar to ΔAG = 0.5 mm (see in Figure 14b). The comparison of contact stress considering the four assembly errors under tooth tip chipping is displayed in Figure 15. When tooth tip chipping occurs, the area of the chip does not bear the load force, thereby increasing the contact stress near the edge of the chipping area. For instance, with ΔAG = 0.5 mm, the maximum contact stress reaches 1356.5 MPa. In contrast, the maximum contact stress for a healthy assembly is only 1186.9 MPa.
The effects of different sizes of chips may lead to different changes in mesh characteristics. Figure 16 showcases the TVMS for different chip areas under healthy installation and E = 0.5 mm. Figure 17 illustrates the unloaded transmission error and contact stress for different chip areas. There are two sizes of tooth tip chipping to compare and analyze, Chip_1 and Chip_2 (see Figure 17). Firstly, when the installation is healthy, the maximum reduction of Chip_1 is 0.91% and this appeared at moment K (see in Figure 16a), and the maximum reduction of Chip_2 is 3.56% and this appeared at moment J (see in Figure 16a) under healthy installation. Following this, when the installation is ΔE= 0.5 mm, the reductions at moments F, G, H, I are listed in Table 3 under ΔE = 0.5 mm. The max reduction of Chip_1 is 0.8% and this appeared at moment H (see in Figure 16b) and the max reduction of Chip_2 is 4.92% and this appeared at moment G (see in Figure 16b). It is simple to deduce that assembly errors increase the reduction of TVMS under tooth tip chipping. Why Chip_1 decreases the healthy installation is due to single tooth contact when Chip_1 happens.
As shown in Figure 17a, the larger the chip area, the greater the reduction in unloaded transmission errors. The OQ and RT are the moments of the affected regions by Chip_1 and Chip_2, respectively. For instance, the maximum reduction in unloaded transmission errors in Chip_1 is 48.3884 μrad, which is smaller than the 131.3402 μrad observed in Chip_2. The maximum contact stress of Chip_1 (see in Figure 17b) is 1204.7 MPa and Chip_2 (see in Figure 17b) is 1324.7 MPa. Additionally, it is apparent that the contact moment affected by Chip_1 (OQ) is less than that of Chip_2 (RT).
In this study, an enhanced method is proposed for analyzing the mesh characteristics of spiral bevel gears (SBGs), taking into account assembly errors and tooth tip chipping. This method is primarily based on tooth contact analysis (TCA) and loaded tooth contact analysis (LTCA) and has been validated using existing literature. The study explores the effects of assembly errors and tooth tip chipping and the conclusions can be summarized as follows:
- (
1 ) The enhanced method has been validated through the consistency of variation trends observed between this paper and previously published literature, considering different assembly errors. - (
2 ) When the contact area moves to the tip, the peak-to-peak value of time-varying meshing stiffness (TVMS) will reduce. In addition, the TVMS has different kinds of increases or decreases with the increase of assembly errors. - (
3 ) Assembly errors lead to an escalation in contact stress at the tooth tip, amplifying the risk of tooth tip chipping. Tooth tip chipping decreases the TVMS and unloaded transmission error within the affected region, whilst concurrently increasing the maximum contact stress. Furthermore, for identical chips, the TVMS considering assembly errors exhibits a more pronounced decline compared to a healthy installation.
Graph: Figure 1 Tooth tip chipping of spiral bevel gear.
Graph: Figure 2 Coordinate systems of SBG with assembly errors.
Graph: Figure 3 Contact area of (a) healthy installation and (b) ΔAG = 0.5 mm.
Graph: Figure 4 A three-tooth FEM model of SBG.
Graph: Figure 5 Contact area and chipping area on the rotating projection of the pinion tooth surface, the actual shape of tooth tip chipping and the clearances change: (a) actual chip; (b) contain chipping and ΔAG = 0.5 mm (c) clearances change due to the chipping.
Graph: Figure 6 Flow chart of LTCA considering assembly error and tooth tip chipping.
Graph: Figure 7 Contact trajectories of different assembly errors on the rotating projection of the gear tooth surface: (a) ΔAP; (b) ΔAG; (c) ΔE; (d) ΔΣ.
Graph: Figure 8 Comparison of contact trajectories due to the assembly errors in Han's paper [[
Graph: Figure 9 TVMS of SBG considering the different assembly errors: (a) ΔAP; (b) ΔAG; (c) ΔE; (d) ΔΣ.
Graph: Figure 10 Unloaded transmission error of SBGs considering different assembly errors: (a) ΔAP; (b) ΔAG; (c) ΔE; (d) ΔΣ.
Graph: Figure 11 Contact stress of SBGs considering different assembly errors.
Graph: Figure 12 TVMS of different assembly errors considering tooth tip chipping: (a) ΔAP = −0.5 mm; (b) ΔAG = 0.5 mm; (c) ΔE = 0.5 mm; (d) ΔΣ = −0.4°.
Graph: Figure 13 The clearances and distributed force between potential contact points at moment C: (a) the clearances; (b) the distributed force.
Graph: Figure 14 Unloaded transmission error of different assembly errors considering tooth tip chipping: (a) ΔAP = −0.5 mm; (b) ΔAG = 0.5 mm; (c) ΔE = 0.5 mm; (d) ΔΣ = −0.4°.
Graph: Figure 15 Contact stress of different assembly errors considering tooth tip chipping: (a) ΔAP = −0.5 mm; (b) ΔAG = 0.5 mm; (c) ΔE = 0.5 mm; (d) ΔΣ = −0.4°.
Graph: Figure 16 TVMS of different chip areas under healthy installation and E = 0.5 mm: (a) healthy installation; (b) E = 0.5 mm.
Graph: Figure 17 Unloaded transmission errors and contact stress of different chip areas under E = 0.5 mm: (a) unloaded transmission errors; (b) contact stress.
Table 1 Basic geometry parameters of the example gear pair.
Parameters Pinion Gear Tooth number 26 31 Module (mm) 8.654 Shaft angle (°) 90 Mean spiral angle (°) 35 Face width (mm) 57 Mean cone distance (mm) 146.57 Hand of spiral Left Right Pitch angles (°) 39.98 50.01 Root angles (°) 37.13 46.17 Face angles (°) 43.83 52.87 Addendum (mm) 8.27 6.44 Dedendum (mm) 8.47 10.31
Table 2 Machining parameters of the example gear pair.
Parameter Pinion Gear Concave Convex Concave Convex Cutter point radius (mm) 108.58 112.05 116.45 112.15 Pressure angle (°) −17.50 22.50 −18.50 21.50 Root fillet radius (mm) 1.10 1.10 2.50 2.50 Machine center to back (mm) −4.34 6.90 0 0 Sliding base (mm) 2.14 −4.68 1.46 1.46 Blank offset (mm) −2.50 2.5 0 0 Radial distance (mm) 116.73 127.5 123.81 123.81 Machine root angle (°) 37.38 37.38 46.17 46.17 Cradle angle (°) −51.07 −53.00 49.13 49.13 Velocity ratio 1.50 1.62 1.30 1.30 Modified Roll Coefficient C 0.0197 −0.0231 0 0 Modified Roll Coefficient D −0.0172 0.0477 0 0
Table 3 Reduction of TVMS of two chip areas at different moments under E = 0.5 mm.
Moment F Moment G Moment H Moment I Chip_1 0 −0.43% −0.8% −0.24% Chip_2 2.31% −4.92% −4.39% −1.94%
Conceptualization, Z.L. and R.L.; Software, J.Z.; Formal analysis, Y.W.; Investigation, J.Z.; Resources, Z.L.; Project administration, Y.W.; Funding acquisition, Z.L., R.Z. and H.M. All authors have read and agreed to the published version of the manuscript.
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Author Ying Wang was employed by the company AECC Harbin Dongan Engine Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Modulus of elasticity of pinion and gear, respectively Distributed normal force at potential contact point Distributed normal force of potential contact points Distance between the two adjacent potential contact points First three-order transformation matrix of matrix Coordinate transformation matrix from Normal vector of one point at pinion and gear coordinate system, respectively and Normal vector of the unloaded contact point on the pinion and gear tooth face, respectively Three-point curve normal vector Tooth surface normal vector of potential contact point at meshing coordinate system Coordinates of the unloaded contact point at pinion and gear coordinate system, respectively and Coordinates of the unloaded contact point on the pinion and gear tooth face, respectively R Mean cone distance ste Static transmission error Auxiliary coordinate system of pinion and gear, respectively Pinion and gear coordinate system, respectively T Input torque Relative velocity of the virtual generating gear face and machined gear face represented in S2 Tooth numbers of the pinion and gear, respectively Greek symbols ΔAP Axial assembly error of pinion ΔAG Axial assembly error of gear ΔE Offset error ΔΣ Shaft angle error Σ Shaft angle Rotor angle of the pinion and gear, respectively Pressure angle Mean spiral angle Pitch angle of the spiral bevel gear Clearance of the Clearance between corresponding contact points at pinion and gear at meshing moments Rotation angles of driving pinion at different meshing moments Three-point curvature and surface normal curvature, respectively and Surface normal curvature at potential contact point of pinion and gear, respectively Global compliance matrix Contact compliance matrix Poisson ratio of pinion and gear, respectively Abbreviations ETCA Error Tooth Contact Analysis FEA Finite Element Analysis FEM Finite Element Method LTCA Loaded Tooth Contact Analysis SBG Spiral bevel gear TCA Tooth contact analysis TVMS Time-varying meshing stiffness
By Ying Wang; Juntao Zhang; Zhanwei Li; Ruijun Liang; Rupeng Zhu and Hui Ma
Reported by Author; Author; Author; Author; Author; Author