In this thesis, I focus on the tooth modification design theory and methodology for automotive straight bevel gears, aiming at improving the meshing performance, reducing vibration and noise, and optimizing the stress distribution on tooth flanks. The research combines precise three-dimensional modeling, dynamic finite element simulation, and gear geometry analysis. I first derived the spherical involute equation in a Cartesian coordinate system and built accurate solid models of straight bevel gears using a combined Matlab and SolidWorks approach. Then, I established a dynamic finite element model based on the ANSYS/LS-DYNA platform to simulate the meshing process. Furthermore, I investigated profile modification and tooth trace modification techniques. A novel concept called isometric tooth trace modification is proposed, and its effectiveness is verified by numerical simulation. The main results show that both profile modification and tooth trace modification can significantly reduce meshing impact and improve the uniformity of tooth load distribution. Throughout the study, the key phrase straight bevel gear is repeatedly examined from modeling, simulation, and design aspects.
The automotive differential is one of the most critical components of a vehicle. Its performance directly affects the overall driving quality, fuel efficiency, and safety. The gear transmission system inside the differential is mainly composed of straight bevel gear pairs. Unlike cylindrical gears, the tooth flanks of a straight bevel gear are generated by spherical involutes, which are spatial curves. Because of the complexity of spherical involute geometry, many existing models approximate the tooth profile by circular arcs or planar involutes. In my opinion, such approximations lead to significant modeling errors and may mislead the subsequent finite element analysis. To overcome this problem, I developed a parameterized modeling method based on the exact spherical involute equation. In this method, the gear model is not built by Boolean operations of simple primitives but by constructing the true tooth flank surfaces from the mathematical definition of the spherical involute. Therefore, the model precision is greatly improved, and any modification to the design parameters, such as pressure angle, number of teeth, module, and cone distance, can be implemented quickly by regenerating the coordinate data.
The dynamic meshing of a straight bevel gear pair is characterized by time-varying contact stiffness, nonlinear contact force, and possible edge contact due to misalignment or deformation. Traditional static analysis cannot capture these transient behaviors. In this work, I adopted the explicit dynamic finite element code LS-DYNA to simulate the meshing process. The simulation model includes the entire gear body, not just a single tooth, because the compliance of the gear body has a non-negligible influence on the local tooth deformation. I carefully defined the material properties, element types, contact algorithm, and boundary conditions for the straight bevel gear pair. The gear material is 20CrMnTi steel with density \(7.8\times 10^3 \,\mathrm{kg/m^3}\), elastic modulus \(2.07\times 10^5 \,\mathrm{MPa}\), and Poisson’s ratio \(0.3\). These parameters are typical for automotive bevel gears. I selected SOLID164 elements for the gear bodies and SHELL163 elements for the inner rigid rings. The inner rings are used to apply the rotational velocity and torque while sharing nodes with the gear body elements. The contact between the two gears is modeled by the surface-to-surface penalty algorithm with a static friction coefficient of 0.5 and a dynamic friction coefficient of 0.3. In order to reduce the computation time without losing accuracy, I refined the mesh in the tooth flank and tooth root regions, while the rest of the gear body uses a coarser mesh. The final model contains about 260,000 solid elements and about 3,500 shell elements. The meshed straight bevel gear pair is shown in the following figure.
In the simulation, the driving gear (planetary gear) is given a constant angular velocity of \(157 \,\mathrm{rad/s}\), and the driven gear (half shaft gear) is subjected to a constant resistance torque of \(10 \,\mathrm{N\cdot m}\). The simulation time is set to 0.004 seconds, which is sufficient to observe the meshing process of one or two teeth. The load is ramped up gradually within a short time to avoid initial shock waves. The output files include the binary result files d3plot and the time-history files, which are post-processed by LS-PREPOST. Through the simulation, I can extract the instantaneous contact stress, the von Mises stress distributions, the nodal displacement, velocity, and acceleration histories. These quantities are essential for evaluating the dynamic behavior of the straight bevel gear pair before and after tooth modification.
One of the most challenging tasks in the finite element analysis of bevel gears is the accurate generation of the tooth flank geometry. For a standard straight bevel gear, the tooth flank is a part of the spherical involute surface. The derivation of the surface equation starts from the rolling of a tangent plane over the base cone. Consider a base cone with base cone angle \(\theta_b\). When a plane rolls over the cone without slipping, a line in the plane generates the so-called spherical involute surface. To write the equation in a fixed Cartesian coordinate system, I define a moving frame \((x’, y’, z’)\) attached to the plane. The conversion from the moving frame to the fixed frame yields:
\[
\begin{cases}
x = x’\cos\phi – y’\sin\phi, \\[4pt]
y = (x’\sin\phi + y’\cos\phi)\cos\theta – z’\sin\theta, \\[4pt]
z = (x’\sin\phi + y’\cos\phi)\sin\theta + z’\cos\theta,
\end{cases}
\]
where \(\phi\) is the rotation angle of the plane around the cone axis, and \(\theta\) is the cone angle. For a line \(OA\) in the moving plane with an inclination angle \(\varphi\) relative to the instantaneous rolling axis, the coordinates of the line in the moving frame are \(x’ = z’\tan\varphi\), \(y’ = 0\). Substituting these values into the above transformation, I obtain the spherical involute surface equation. By imposing the sphere condition \(x^2 + y^2 + z^2 = l^2\), I get the explicit Cartesian equation of the spherical involute on a sphere of radius \(l\):
\[
\begin{cases}
x = l\left(\cos(\phi\sin\theta_b)\sin\theta_b\cos\phi + \sin(\phi\sin\theta_b)\sin\phi\right), \\[4pt]
y = l\left(\cos(\phi\sin\theta_b)\sin\theta_b\sin\phi – \sin(\phi\sin\theta_b)\cos\phi\right), \\[4pt]
z = l\cos(\phi\sin\theta_b)\cos\theta_b,
\end{cases}
\]
where \(\theta_b\) is the base cone angle, and \(\phi\) is the parameter varying from 0 to a maximum value corresponding to the tooth tip. The base cone angle is related to the pitch cone angle \(\delta\) and the pressure angle \(\alpha\) by the relationship \(\sin\theta_b = \sin\delta\cos\alpha\). For the differential gears used in this study, the planetary gear has 10 teeth and the half shaft gear has 15 teeth, with a module of 4.438 mm, pressure angle \(22.5^\circ\), and shaft angle \(90^\circ\). The pitch cone angles are \(33.69^\circ\) for the planetary gear and \(56.31^\circ\) for the half shaft gear. The cone distance is \(R = 40 \,\mathrm{mm}\). These parameters are listed in Table 1.
| Parameter | Planetary gear | Half shaft gear |
|---|---|---|
| Number of teeth | 10 | 15 |
| Module (mm) | 4.438 | |
| Pressure angle (deg) | 22.5 | |
| Shaft angle (deg) | 90 | |
| Pitch diameter (mm) | 44.38 | 66.57 |
| Cone distance (mm) | 40 ± 0.025 | |
| Pitch cone angle (deg) | 33.69 | 56.31 |
| Face cone angle (deg) | 46 | 66 |
| Root cone angle (deg) | 24 | 44 |
| Addendum (mm) | 4.98 | 3.28 |
| Whole depth (mm) | 9.17 | 9.169 |
With the exact spherical involute equation, I can generate the coordinates of points on both the large-end and small-end spheres. In Matlab, I wrote a script that takes the base cone angle, the sphere radius, and the angular range as inputs, and outputs the point coordinates. I then imported these coordinates into SolidWorks via the “Curve through XYZ points” function. The large-end and small-end spherical involutes were connected by the tooth addendum and dedendum lines to form a closed boundary. By using the loft or boundary surface operation, I constructed the tooth flank surface. The opposite flank is obtained by mirroring the first flank with respect to a symmetry plane that contains the tooth center line. The tooth root surface is generated by a ruled surface connecting the root lines and being tangent to the flank surface. The addendum surface is a conical surface. The combination of all these surfaces forms a closed volume, which represents a tooth space. The complete gear is obtained by multiplying the tooth space around the axis by the number of teeth and subtracting it from the main gear body. The final solid model of the half shaft gear and the planetary gear is shown in Figure 2. This modeling process is fully parametric, because the only free variables are the gear parameters. When the parameters are changed, the Matlab script regenerates the coordinates and the SolidWorks model can be rebuilt automatically through a macro, or by simply updating the curve data. This method avoids the common approximations that use planar involutes or circular arcs for the tooth flank, and ensures that the finite element model is geometrically correct.
The design of tooth modification for the straight bevel gear is a key issue in this work. Tooth modification can be divided into two categories: profile modification and tooth trace modification. Profile modification refers to removing a small amount of material from the tip or root of the tooth along the involute profile direction. This is necessary because the actual tooth base pitch cannot be perfectly equal due to manufacturing errors, assembly errors, and elastic deformations under load. When the base pitches of the two meshing gears are not equal, the meshing point will not move along the theoretical line of action, causing instantaneous changes in the transmission ratio and resulting in additional dynamic loads, vibration, and noise. In Figure 3, I illustrate the concept of profile modification. The line BC represents the starting line of profile modification, S1 and S2 are the modification amounts at the small end and large end, respectively, and L1 and L2 are the modification heights. The dashed line shows the modified tooth profile, which is a circular arc in this study.
To determine the profile modification amount, I used the formula \( \delta_a = \frac{F_t}{b} \cdot c_r\), where \( \delta_a\) is the elastic deformation of the tooth flank in micrometers, \(F_t\) is the tangential force, \(b\) is the tooth width, and \(c_r\) is the meshing stiffness per unit width, approximately \(20 \, \mathrm{N/(mm\cdot \mu m)}\). The deformation can also be computed from the finite element results by taking the difference of the circumferential displacements between a loaded tooth and an unloaded tooth at the same radial position. For the large end, the modification amount \(S_2\) is obtained directly. The small end modification amount \(S_1\) is then calculated from the geometric relation \(S_1 = S_2 \cdot (R-b)/R\), where \(R\) is the cone distance and \(b\) is the face width. For the planetary gear and half shaft gear, I obtained the profile modification amounts as shown in Table 2.
| Gear | S1 (small end) μm | S2 (large end) μm |
|---|---|---|
| Half shaft gear | 17.33 | 29.03 |
| Planetary gear | 18.06 | 27.69 |
The profile modification height is also important. There are two types: long modification and short modification. Long modification removes the entire double-pair contact zone, which is only suitable for gears with large helix angles and large overlap ratios. For a conventional straight bevel gear with a small overlap ratio, short modification is preferred. In short modification, only a part of the double-pair contact zone is removed, and the load transition becomes smoother. Based on the gear parameters, I selected the short modification height as listed in Table 3.
| Gear | L1 (small end) mm | L2 (large end) mm |
|---|---|---|
| Planetary gear | 1.79 | 3.0 |
| Half shaft gear | 2.4 | 3.0 |
In addition to the modification amount and height, the shape of the modification curve determines the smoothness of the transition. Many researchers have proposed different curves, such as the Walker curve \(x\Delta = \Delta_{\max}(x/l)^{1.5}\), the Hidaka-Terauchi curve \(x\Delta = \Delta_{\max}(x/l)^{1.22}\), and the simple linear curve. In general, they can be expressed as:
\[
\Delta(x) = \Delta_{\max} \left(\frac{x}{l}\right)^n
\]
where \(\Delta_{\max}\) is the maximum modification amount, \(l\) is the modification height, \(x\) is the distance from the tooth tip, and \(n\) is the curve exponent. For the straight bevel gear, the involute is a spherical curve, so planar involute modification curves cannot be directly used. Circular arcs are easy to model and manufacture, and they provide a smooth transition. I therefore used a circular arc as the modification curve. In the SolidWorks environment, I created the modified tooth flank by deleting the original surface region and replacing it with a surface generated by the arc. The procedure is as follows: first, determine the starting point B of the modification on the involute. Second, determine the tip point A based on the modification amount S1. Third, construct the perpendicular bisector of line AB and find its intersection with the normal line of the involute at point B. This intersection is the center of the circular arc. The arc from A to B is the modification curve. This curve is then used to generate the new tooth flank surface. The modified gear model was then imported into ANSYS/LS-DYNA for the dynamic simulation.
After applying profile modification to the straight bevel gear pair, I compared the stress distributions and dynamic responses before and after modification. The contact stress on the half shaft gear decreased from about \(3.5 \times 10^9 \, \mathrm{Pa}\) to \(2.9 \times 10^9 \, \mathrm{Pa}\), which corresponds to a reduction of 17.1%. For the planetary gear, the contact stress dropped from \(2.4 \times 10^9 \, \mathrm{Pa}\) to \(1.7 \times 10^9 \, \mathrm{Pa}\), a reduction of 29.1%. More importantly, the angular acceleration and axial acceleration of the gear teeth were significantly reduced, indicating that the meshing impact was mitigated. The time-history curves of the angular acceleration for a node on the tooth tip are shown in Figure 5. The peak-to-peak amplitude of the angular acceleration was reduced by about 30% after profile modification. The axial acceleration of the planetary gear also decreased, which is beneficial for reducing the axial vibration of the differential. However, the axial acceleration of the half shaft gear was relatively small to begin with, so the improvement was not as obvious. These results confirm that profile modification is an effective way to reduce vibration and noise in straight bevel gear transmissions.
In the second part of my work, I studied tooth trace modification, also known as longitudinal modification. Due to the elastic deflections of the gear shafts and the deformations of the gear bodies, the contact along the tooth width direction is not uniform in practice. The load tends to concentrate at one end of the tooth, causing edge loading and premature failure. The conventional solution is to crown the tooth, i.e., to make the tooth flank slightly barrel-shaped. A crowned tooth has a reduction in the tooth thickness from the center to the ends, so that the contact area is initially concentrated at the center and gradually spreads as the load increases. The crowning amount is usually determined by empirical formulas or by calculating the composite deformation of the gear system. For cylindrical gears, crowning is widely used. However, for a straight bevel gear, crowning is difficult to achieve by ordinary machine tools, and the position of the crown is hard to control. In the differential, the contact pattern is normally desired near the small end of the tooth, but traditional crowning often places the crown at the middle of the face width, which is not optimal.
To solve this problem, I proposed a new method called isometric tooth trace modification for the straight bevel gear. The principle is to remove a uniform layer of material along the normal direction of the original spherical involute surface in a selected local region of the tooth flank. The modified surface is parallel to the original surface but shifted inward. Since the shift is uniform in the normal direction, the new tooth profile remains a spherical involute; only the effective thickness is reduced. The contact area between the meshing teeth will then be limited to the unmodified region or to the region where the two surfaces first come into contact. By choosing the location and the size of the isometric modification, one can control the contact pattern on the tooth flank. In a gear pair, the pinion (planetary gear) is usually more flexible and more prone to bending deformation, so I applied the isometric modification only to the planetary gear, leaving the half shaft gear unchanged. This is a practical assumption because the pinion has a smaller number of teeth and a thinner tooth shape.
To define the isometric modification region, I used a rectangular patch on the tooth flank, described by the parameters \(a\), \(b\), \(c\), and \(d\) as shown in Figure 6. The parameter \(a\) is the distance from the small end to the edge of the modified region, \(b\) is the length of the modified region along the tooth width, \(c\) is the distance from the tooth root to the lower boundary, and \(d\) is the distance from the upper boundary to the tooth tip. According to the bevel gear accuracy standard and the practical experience, the ideal contact area for the automotive differential should occupy about 50% to 70% of the tooth length and 55% to 75% of the tooth height. I chose the modified region to be 60% of the tooth length and 50% of the tooth height, placing it close to the small end. Specifically, I set \(a = 2.3 \, \mathrm{mm}\), \(b = 3.2 \, \mathrm{mm}\), \(c = 2.3 \, \mathrm{mm}\), and \(d = 2.4 \, \mathrm{mm}\).
The amount of isometric modification, denoted by \(h\), must be chosen carefully. If \(h\) is too large, the tooth will not have enough material to sustain the load; if it is too small, the contact pattern will not be sufficiently controlled. I determined \(h\) using the finite element method by comparing the angular displacement of a node on the loaded tooth with that of a node on an unloaded tooth. The difference in the circumferential displacement gives the elastic deformation of the tooth in the normal direction. For the planetary gear, this value was found to be 30 μm. To verify that this value is compatible with the gear kinematic accuracy, I derived a relationship between the isometric modification amount and the resulting angular error. For a gear rotating about its axis, if a point on the tooth flank is moved inward by a distance \(h\) along the normal, the angular error of the gear can be approximated by:
\[
\Delta\Phi = \frac{h}{R_b}
\]
where \(R_b\) is the base circle radius of the equivalent cylindrical gear. Substituting \(h = 30\,\mu\mathrm{m}\) and \(R_b \approx 38\,\mathrm{mm}\), I got \(\Delta\Phi \approx 7.9 \times 10^{-4} \,\mathrm{rad}\), which is within the allowable one-tooth tangential composite error for the 7th accuracy grade according to GB/T11365-1989. Therefore, the selected modification amount is reasonable.
The modified tooth flank is generated by offsetting the original spherical involute surface along its normal by a constant distance \(h\) in the selected region. Practically, I used SolidWorks to create a 3D sketch of the boundary of the modified region on the tooth flank. Then I created a surface offset by 0.03 mm from the original flank and trimmed it using the boundary sketch. The modified tooth flank is the trimmed offset surface. The rest of the tooth flank remains unmodified. This process is easy to implement and does not require special machine tool adjustments, which is a major advantage over crowning. Moreover, the position of the contact area can be adjusted by simply moving the rectangular boundary patch. This flexibility is very useful for different load conditions and accuracy grades.
To evaluate the effect of the isometric tooth trace modification, I performed a dynamic finite element simulation of the modified straight bevel gear pair. The same boundary conditions and material properties were used as in the unmodified case. The contact stress distribution along the tooth width direction was extracted at the moment when the tooth is in the single-tooth contact zone. Figure 7 presents the von Mises stress contour on the planetary gear before and after modification. In the unmodified case, the stress concentrated near the large end of the tooth, with a maximum value of about 2400 MPa at the edge. After isometric modification, the stress concentration at the large end was significantly relieved. The contact area moved to the central part of the modified region, and the stress distribution became more symmetric with respect to the middle of the contact patch. It is interesting to note that the maximum contact stress on the planetary gear slightly increased after modification, because the contact area became smaller. The maximum value rose from about 2400 MPa to about 2500 MPa. However, the stress concentration factor, defined as the ratio of the maximum stress to the average stress along the tooth width, decreased from 1.5 to 1.15. This indicates that the load distribution is more uniform, which is beneficial for avoiding local damage such as pitting and spalling.
For the half shaft gear, which was not modified, the contact stress distribution along the tooth width also changed accordingly. Because the contact area on the pinion was moved toward the small end, the half shaft gear received a more balanced load. The maximum stress on the half shaft gear decreased slightly, and the edge loading at the large end disappeared. I also compared the stress histories at the tooth root and at the contact region. The results are summarized in Table 4.
| Gear | Location | Before modification | After modification | Change |
|---|---|---|---|---|
| Planetary gear | Tooth root | 960 | 650 | -32% |
| Contact region | 1000 | 2100 | +110% | |
| Half shaft gear | Tooth root | 850 | 780 | -8.2% |
| Contact region | 1200 | 1100 | -8.3% |
The increase of the contact stress on the planetary gear may be acceptable if the contact stress is still below the allowable limit. In practice, the hardened case on the tooth flank can withstand high contact stress. The significant reduction of the tooth root stress is very important because tooth root bending fatigue is a common failure mode for straight bevel gear teeth. The reduction of 32% in the tooth root stress indicates that the isometric modification not only improves the contact pattern but also reduces the bending stress. This is because the contact area is moved to a more favorable position, where the load arm is shorter and the tooth is thicker.
In addition to the stress analysis, I also studied the dynamic response of the modified gear pair. The angular acceleration history of a node on the planetary gear tooth tip is shown in Figure 8. The peak angular acceleration was reduced from \(1500 \,\mathrm{rad/s^2}\) to \(1000 \,\mathrm{rad/s^2}\) after isometric modification. The axial acceleration also decreased slightly. These results demonstrate that the isometric tooth trace modification is effective in reducing the vibration of the straight bevel gear transmission. The method is particularly suitable for precision forged bevel gears, because the modification can be directly applied to the electrode gear used in the forging die. By modifying the electrode gear, the resulting die cavity has the desired tooth flank shape, and the forged gear product will have the same modified surface. This is much more efficient and cost-effective than hiring special grinding or lapping machines.
Let me now summarize the theoretical derivation for the influence of the isometric modification on the gear kinematics. In a gear pair, the meshing point on the tooth flank is normally determined by the intersection of the line of action and the tooth surface. If the tooth surface is offset along its normal by a small constant \(h\), the new surface is parallel to the original. In differential geometry, the normal vector at any point on the offset surface remains the same as the original surface. The position vector of a point on the modified surface can be written as:
\[
\mathbf{r}’ = \mathbf{r} + h\mathbf{n}
\]
where \(\mathbf{r}\) is the position vector on the original surface and \(\mathbf{n}\) is the unit normal vector. Since \(h\) is constant, the derivative of \(\mathbf{r}’\) with respect to the surface parameters is \(\partial \mathbf{r}’ / \partial u = \partial \mathbf{r} / \partial u + h \partial \mathbf{n} / \partial u\). The second fundamental form of the offset surface changes by \(h\) times a term related to the curvature. For a small modification amount, the change in the local curvature is negligible. Therefore, the new surface is still a valid involute surface that can mesh with a conjugate gear, provided that the offset does not cause any interference. This is an important advantage of isometric modification over arbitrary profile modification, because it maintains the correctness of the tooth engagement.
I also derived the influence of the modification on the transmission error. The transmission error of a gear pair is defined as the difference between the actual angular position of the driven gear and the theoretical position based on the ideal gear ratio. When the tooth surface of the driving gear is offset by \(h\), the normal deviation along the line of action causes an angular error on the driven gear. To a first-order approximation, the angular error is given by:
\[
\theta_e = \frac{h}{R_b \cos\alpha}
\]
where \(\alpha\) is the pressure angle. This formula can be used to check whether the modification amount violates the accuracy requirement. For our parameters, \(\theta_e\) is less than 20 arc-seconds, which is acceptable for automotive differentials.
It is worth noting that the proposed isometric modification is different from the conventional crowning. In crowning, the tooth thickness varies continuously from the center to the ends, and the contact point is at the center of the crown. In isometric modification, the tooth flank is locally recessed in a patch, and the contact area is controlled by the boundary of the patch. The contact point is initially at the boundary of the modified region, and as the load increases, the contact area expands elastically. This allows the designer to place the contact area at the desired position, such as near the small end of the tooth. For a straight bevel gear in a differential, the small end is usually where the gear meshes with the pinion in the most favorable manner. The method is also easier to implement in the forging process. In the production of precision forged straight bevel gear, the die cavity is produced by electrical discharge machining using an electrode gear. If the electrode gear has the isometric modified tooth flank, the cavity will have the negative form of the modification, and the desired product geometry is obtained.
To fully understand the influence of the modification parameters, I performed a series of simulations with different modification amounts and positions. The results show that the optimal modification amount depends on the load level. For a lower torque, a smaller modification amount is sufficient to avoid edge contact, while for a higher torque, a larger modification amount is necessary to accommodate the larger elastic deformation. The position of the modified region should be chosen such that the remaining unmodified area is large enough to support the full load. In my simulations with a torque of 10 N·m, the optimal amount was found to be 30 μm. I also tested 20 μm and 40 μm. With 20 μm, the edge contact at the large end was not completely eliminated, and the stress distribution still showed a peak at the large end. With 40 μm, the contact area became too small, causing a very high local stress and an increase in the transmission error. Therefore, 30 μm is a good compromise for this gear pair under this load.
Another important factor is the transition between the modified and unmodified regions. In the ideal case, the offset surface should be tangent to the original surface at the boundary to avoid a sharp edge. If the boundary is sharp, it creates a stress concentration and may also cause noise during the meshing. In my model, I introduced a small fillet or a smooth transition region by making the offset amount gradually decrease to zero at the boundary. This can be achieved by defining the offset amount as a function of the position along the boundary. For example, the offset height \(h(u,v)\) could be constant in the central part and then decrease linearly or parabolically to zero within a small transition band. This smooth transition reduces the stress concentration at the edge of the modified region. In the finite element model, I realized the smooth transition by creating a set of surfaces with varying offset distances, or by using a tetrahedral mesh with a variable thickness. The simulation results confirmed that the smooth transition significantly reduces the local stress at the boundary. The maximum stress in the transition zone was about 15% lower than that in the sharp-edge model.
The dynamic simulation of the straight bevel gear pair also revealed the phenomenon of “end contact” in the unmodified gear. At the beginning of the meshing, the contact line is not parallel to the tooth trace but is inclined, because the gear axis is not perfectly aligned due to the shaft deflection. This inclination causes the load to concentrate at one end of the tooth. The isometric modification effectively moves the contact area away from the edge and toward the central part of the tooth. The contact pattern observed in the simulation matches well with the theoretical design: the contact area is centered at the middle of the modified region, and the length of the contact area is about 60% of the tooth width, which is within the recommended range for 7th-grade bevel gears. The height of the contact area is about 50% of the tooth height. This contact pattern is stable throughout the meshing cycle, and the fluctuation of the contact area is reduced.
In addition to the isometric modification, I also investigated the combination of profile modification and tooth trace modification. Since the gear tooth is subjected to both the bending deformation in the radial direction and the torsional deformation in the tooth width direction, a comprehensive modification strategy is often required. In my simulations, I applied both profile modification and isometric tooth trace modification to the planetary gear. The combined modification resulted in a further reduction of the tooth root stress and a smoother acceleration curve. The maximum contact stress did not increase significantly compared with the tooth-trace-only modification because the profile modification reduced the edge impact at the tooth tip. The combined modification is recommended for practical applications where both the meshing impact and the load distribution need to be optimized.
The finite element model established in this thesis can be further extended to study the influence of manufacturing and assembly errors on the performance of the straight bevel gear. In the current model, I assumed an ideal gear pair without any misalignment. However, in real differentials, the gear shafts may not be perfectly perpendicular, and there may be axial play and radial runout. These errors can be introduced in the model by adjusting the position and orientation of the two gears. For example, I can rotate the half shaft gear by a small angle around its axis to simulate the axial clearance, or translate it along the axis to simulate the assembly error. The dynamic simulation with these errors will provide more realistic insights into the tooth modification design. Some preliminary simulations with an angular misalignment of 0.1° showed that the unmodified gear suffers from severe edge contact, while the isometric modified gear can tolerate the misalignment because the contact area is already away from the edge. This is a practical advantage of the isometric modification.
Another extension is the multi-body dynamic analysis of the entire differential, including the gear shafts, bearings, and the differential housing. In that case, the finite element model of the gear pair can be coupled with a lumped-parameter model of the bearings and shafts. The reaction forces at the gear supports can be extracted from the gear contact simulation and then fed into the system model. This is useful for predicting the vibration level and noise radiation of the differential. However, this is beyond the scope of the current thesis and will be addressed in future work.
In conclusion, my research provides a comprehensive study of the tooth modification technology for the straight bevel gear. The main contributions are summarized as follows:
- I derived the exact spherical involute equation and developed a parameterized modeling method using Matlab and SolidWorks. This method eliminates the modeling errors caused by approximate tooth profiles and can quickly generate models for different gear parameters.
- I established a dynamic finite element model of the straight bevel gear pair based on ANSYS/LS-DYNA. The model uses a suitable combination of solid and shell elements, a refined mesh in the critical areas, and realistic boundary conditions. The simulation results provide the stress, strain, and acceleration histories of the gear meshing process.
- I studied profile modification for the straight bevel gear and determined the modification parameters based on elastic deformation and dynamic simulation. The results show that profile modification can reduce the contact stress and the angular acceleration oscillations, thereby improving the meshing smoothness and reducing noise.
- I proposed a new isometric tooth trace modification method. This method creates a parallel offset surface on the tooth flank and is easy to implement in the forging process. I derived the kinematic influence of the isometric modification and verified its effectiveness by finite element simulation. The isometric modification effectively improves the load distribution along the tooth width and reduces the tooth root stress.
- I also showed that the combination of profile modification and isometric tooth trace modification yields the best overall performance. The combined modification is recommended for high-performance automotive differential gears.
There are still some limitations in my work. The finite element model assumes a linear elastic material and does not consider the elastic-plastic behavior of the tooth surface under very high loads. In reality, the gear tooth may undergo plastic deformation at the contact area, especially when the load exceeds the endurance limit. In future work, I plan to include an elastic-plastic material model and a fatigue life prediction algorithm. The thermal deformation of the gear and the effect of lubrication are also neglected in the current model. These factors can affect the actual contact pattern and the temperature rise of the tooth surface. A full thermo-mechanical coupling analysis would be necessary for a more accurate prediction.
Another limitation is that the dynamic simulation was performed for a relatively short time duration of 0.004 seconds, which only covers a few meshing cycles. To study the steady-state response and the statistical properties of the vibration, a longer simulation duration is needed. However, the computational cost increases proportionally. In my future work, I will use a multi-scale approach or a reduced-order model to extend the simulation time while maintaining a reasonable computational cost.
Despite these limitations, the results of this thesis provide a solid foundation for the design and manufacturing of high-quality straight bevel gear with tooth modification. The proposed isometric tooth trace modification is particularly attractive because it is simple, flexible, and compatible with the precision forging process. I hope that the findings of this research will be useful for the automotive industry, especially for the development of lighter and quieter differentials.
In summary, this thesis systematically addressed the key aspects of straight bevel gear modification: precise modeling, dynamic simulation, profile modification, and tooth trace modification. The combination of theoretical derivation, geometric modeling, and finite element analysis provides a comprehensive design tool for engineers. The methods and results presented in this thesis can be directly applied to the design of other bevel gears, such as hypoid gears and spiral bevel gears, with appropriate modifications. The isometric modification concept can also be extended to other types of gear pairs, such as cylindrical gears, by replacing the spherical involute with a planar involute. The underlying principle is universal: by locally removing a uniform layer of material from the tooth flank, one can control the contact pattern and improve the load distribution. This is a valuable alternative to the conventional crowning method, especially when the desired contact area is not at the center of the tooth.
I believe that the technology of straight bevel gear modification will continue to evolve with the advancement of numerical simulation and additive manufacturing. The precise control of the tooth flank geometry opens up new possibilities for optimizing the gear performance under different operating conditions. The methods developed in this thesis are an important step toward the intelligent design of gear transmissions.