In this work, I focus on a practical and difficult engineering problem: how to analyze, with repeatable accuracy, the bearing capacity of the drive gear set used in the main bearing of a large boring machine. The main bearing is the central power output component of the machine, and the drive gear set is one of its most critical transmission assemblies. During operation, the pinion gear and the large ring gear are subjected to strongly alternating loads. If the load-bearing capacity of the pinion gear or the ring gear is underestimated during design, the result can be tooth-root fracture, severe pitting, or complete loss of the main drive function. Because underground repair is extremely difficult, I treat the bearing capacity analysis of the pinion gear and the ring gear as a design-stage requirement rather than a later verification step.
I propose a normative method that combines engineering standard calculation with refined finite element simulation. The method is intended for large gears whose dimensions and tooth geometry make ordinary mesh generation unreliable. In my approach, the pinion gear, the ring gear, the contact region, the tooth-root region, and the supporting structure are all treated in a unified workflow. I use a standard gear load-capacity formulation as the first screening tool, and then I use a refined finite element model to capture contact stress, tooth-root bending stress, and strain distribution. The two results are compared, and the deviation is used as a consistency check. I apply the method to a 6000-type boring machine drive gear set and show that the pinion gear and the ring gear satisfy the required contact and bending strength conditions.
1. Engineering Context and Motivation
Large boring machines are used in metro construction, long-distance tunnels, and other projects where the excavation face must be supported while the cutterhead is driven. The main bearing is the interface between the stationary support structure and the rotating cutterhead. The drive gear set transmits torque from several motors and reducers to the main bearing inner ring. In many designs, the ring gear is integrated with the main bearing inner ring, while multiple pinion gears are arranged around the ring gear. This architecture distributes the drive torque, but it also creates a complex load-sharing problem. Each pinion gear must carry a portion of the total torque, and the ring gear must carry the sum of all pinion gear loads.
The bearing capacity of the drive gear set depends on several factors: gear geometry, material properties, heat treatment, contact ratio, misalignment, support stiffness, load distribution among the pinion gears, and the actual meshing position. In large gears, the tooth thickness difference between the addendum and the dedendum is significant. This makes automatic mesh generation difficult. A poor mesh can produce stress values that differ from the true values by several times. Therefore, I consider the mesh preprocessing method to be as important as the solver settings. My objective is to create a general mesh generation procedure for large gears and a stable constraint and load application procedure for the drive gear set.
2. Analytical Framework for the Drive Gear Set
I begin with the fundamental kinematic and dynamic relations of the drive gear set. The gear ratio between the ring gear and the pinion gear is determined by the tooth numbers. If the ring gear has \(z_g\) teeth and the pinion gear has \(z_p\) teeth, then the ratio is
$$ i = \frac{z_g}{z_p} = \frac{198}{16} = 12.375. $$
The rotational speed of the pinion gear is therefore
$$ n_p = i n_g = 12.375 \times 0.5 = 6.1875 \ \text{r/min}, $$
where \(n_g\) is the ring gear speed. Because the drive system uses multiple pinion gears arranged around the ring gear, the torque per pinion gear is obtained by dividing the total ring gear torque by the number of pinion gears and by the gear ratio. If \(N_p\) is the number of pinion gears and \(T_g\) is the total ring gear torque, then
$$ T_p = \frac{T_g}{N_p i} = \frac{36830}{16 \times 12.375} = 186.010101 \ \text{kN}\cdot\text{m}. $$
This value is the nominal torque carried by one pinion gear under the assumed equal load-sharing condition. The pitch diameter of the pinion gear is
$$ d_p = m z_p = 25 \times 16 = 400 \ \text{mm}, $$
and the pitch diameter of the ring gear is
$$ d_g = m z_g = 25 \times 198 = 4950 \ \text{mm}. $$
The tangential force at the pinion gear pitch circle is
$$ F_t = \frac{2000 T_p}{d_p}. $$
Using \(T_p = 186.010101 \ \text{kN}\cdot\text{m}\) and \(d_p = 400 \ \text{mm}\), I obtain
$$ F_t = \frac{2000 \times 186.010101}{400} = 930.0505 \ \text{kN}. $$
The total tangential force on the ring gear from all pinion gears is approximately
$$ F_{t,\text{total}} = N_p F_t = 16 \times 930.0505 = 14880.808 \ \text{kN}. $$
This total force is consistent with the total ring gear torque and the ring gear pitch diameter:
$$ F_{t,\text{total}} = \frac{2000 T_g}{d_g} = \frac{2000 \times 36830}{4950} = 14880.808 \ \text{kN}. $$
The radial and axial components are zero for a spur gear set, but in a general drive gear set they would be calculated from the pressure angle and helix angle. For completeness, I write the general expressions as
$$ F_r = F_t \tan \alpha_t, $$
$$ F_a = F_t \tan \beta, $$
where \(\alpha_t\) is the transverse pressure angle and \(\beta\) is the helix angle. For the present spur gear case, \(\beta = 0\), so \(F_a = 0\).
For the standard contact stress calculation, I use the form
$$ \sigma_H = Z_H Z_E Z_\epsilon Z_\beta \sqrt{\frac{F_t}{b d_p} \frac{u+1}{u} K_A K_V K_{H\beta} K_{H\alpha}}, $$
where \(\sigma_H\) is the contact stress, \(Z_H\) is the zone factor, \(Z_E\) is the elasticity factor, \(Z_\epsilon\) is the contact ratio factor, \(Z_\beta\) is the helix angle factor, \(b\) is the face width, \(u\) is the gear ratio, \(K_A\) is the application factor, \(K_V\) is the dynamic factor, \(K_{H\beta}\) is the face load distribution factor, and \(K_{H\alpha}\) is the transverse load distribution factor. For the tooth-root bending stress, I use
$$ \sigma_F = \frac{F_t}{b m_n} Y_{Fa} Y_{Sa} Y_\epsilon Y_\beta K_A K_V K_{F\beta} K_{F\alpha}, $$
where \(\sigma_F\) is the tooth-root bending stress, \(m_n\) is the normal module, \(Y_{Fa}\) is the form factor, \(Y_{Sa}\) is the stress correction factor, \(Y_\epsilon\) is the contact ratio factor, \(Y_\beta\) is the helix angle factor, \(K_{F\beta}\) is the face load distribution factor, and \(K_{F\alpha}\) is the transverse load distribution factor. The permissible contact stress is
$$ \sigma_{HP} = \frac{\sigma_{H,\lim}}{S_H} Z_N Z_L Z_R Z_V Z_W Z_X, $$
and the permissible bending stress is
$$ \sigma_{FP} = \frac{\sigma_{F,\lim}}{S_F} Y_{ST} Y_{NT} Y_{\delta relT} Y_{RrelT} Y_X. $$
Here \(S_H\) and \(S_F\) are the safety factors for contact and bending, respectively. The material used in the drive gear set is a carburized and quenched 17CrNiMo6-type steel. Its allowable contact stress is taken as 3000 MPa, and its yield limit is taken as 1350 MPa. These values are used only as screening limits in my standard calculation. The final acceptance criterion is that the calculated contact stress and tooth-root bending stress remain below the corresponding permissible values with adequate safety margin.
3. Engineering Standard Calculation Workflow
I organize the standard calculation as a repeatable workflow. First, I establish the gear geometry from the design parameters. Second, I define the load cases, including the extreme working condition. Third, I apply the constraints that represent the actual support conditions. Fourth, I run the load spectrum and extract the contact stress and tooth-root bending stress. Fifth, I compare the results with the material limits. The design parameters I used for the 6000-type drive gear set are summarized in Table 1.
| Parameter | Pinion Gear | Ring Gear |
|---|---|---|
| Module \(m\) | 25 mm | 25 mm |
| Number of teeth \(z\) | 16 | 198 |
| Pressure angle \(\alpha\) | 20° | 20° |
| Profile shift coefficient \(x\) | +0.5 | -1.0 |
| Face width \(b\) | 275 mm | 275 mm |
| Material | 17CrNiMo6, carburized and quenched | 17CrNiMo6, carburized and quenched |
| Number of pinion gears | 16, approximately equally spaced | — |
| Center distance | 2275 mm | 2275 mm |
The extreme working condition loads are listed in Table 2. I use these values to calculate the nominal pinion gear torque and the corresponding tangential force. The ring gear torque is 36830 kN·m, and the local torque for one meshing position is 2301.875 kN·m. The ring gear speed is 0.5 r/min, and the pinion gear speed is 6.1875 r/min. The pinion gear torque is 186.010101 kN·m.
| Quantity | Value |
|---|---|
| Ring gear torque | 36830 kN·m |
| Ring gear local torque | 2301.875 kN·m |
| Ring gear speed | 0.5 r/min |
| Pinion gear torque | 186.010101 kN·m |
| Pinion gear speed | 6.1875 r/min |
For the standard calculation, I build the complete drive gear set with all 16 pinion gears. I do not simplify the model to a single pinion gear at this stage, because the load-sharing among the pinion gears influences the ring gear stress distribution and the local contact conditions. The ring gear and the main bearing inner ring are treated as an integrated component. To represent the support, I introduce a drive shaft component with an inner bore, and I mount the ring gear into the inner bore. The drive shaft and the ring gear are set as one integrated body. I add five equally spaced conceptual rigid bearings on the outer side of the drive shaft to support rotation. These conceptual bearings can carry radial force, axial force, and tilting moment. I release only the rotational degree of freedom of the ring gear and constrain the other directions. I apply a distributed load over the face width of the ring gear.
For each pinion gear, I create a pinion shaft and mount the pinion gear on it. The pinion shaft and the pinion gear are set as one integrated body. I add one bearing on each side of the pinion shaft to support rotation. I release only the rotational degree of freedom of the pinion gear and constrain the other directions. I apply a distributed load over the face width of the pinion gear. The same procedure is used for all 16 pinion gears. The complete model is then placed in the gearbox coordinate system. Using the ring gear coordinates as the reference, I calculate the coordinates of each pinion gear and assemble the pinion gears around the ring gear. This complete assembly is used for the standard calculation.
The standard calculation results under the extreme working condition are summarized in Table 3. The maximum contact stress is 1171 MPa, which is well below the allowable contact stress of 3000 MPa. The maximum tooth-root bending stress of the pinion gear is 374 MPa, and the maximum tooth-root bending stress of the ring gear is 534 MPa. Both values are well below the yield limit of 1350 MPa. Therefore, the drive gear set satisfies the design bearing capacity requirements according to the standard calculation.
| Result | Value | Allowable | Status |
|---|---|---|---|
| Maximum contact stress | 1171 MPa | 3000 MPa | Satisfied |
| Pinion gear tooth-root bending stress | 374 MPa | 1350 MPa | Satisfied |
| Ring gear tooth-root bending stress | 534 MPa | 1350 MPa | Satisfied |
The standard calculation provides a fast engineering estimate. However, it does not fully resolve the local contact pressure, the stress concentration at the tooth root, or the strain distribution inside the carburized layer. For this reason, I use the standard calculation as the first stage and the refined finite element simulation as the second stage.
4. General Mesh Preprocessing for Large Gears
The main difficulty in the finite element analysis of large gears is that the pinion gear and the ring gear have a large tooth thickness variation from the addendum to the dedendum. If I use a conventional tetrahedral mesh, the number of elements becomes very large, and the mesh quality near the tooth root and the contact surface is often poor. A poor mesh can cause convergence failure or inaccurate stress values. I therefore propose a general mesh preprocessing method for large gears. The method is based on a two-dimensional mapped quadrilateral mesh and a subsequent three-dimensional hexahedral extrusion along the face width.
The mesh generation logic is as follows. First, I create a two-dimensional gear profile for the pinion gear or the ring gear. Second, I divide the tooth profile into a contact strain region, a transition region, and a hub region. Third, I generate a mapped quadrilateral mesh on each region. Fourth, I extrude the two-dimensional mesh along the face width to create a three-dimensional hexahedral structured mesh. Fifth, I check the mesh quality and adjust the local sizing if necessary.
The contact strain region includes the tooth flank and the region immediately beneath it. Because the contact stress and the strain gradient are highest in this region, I refine the mesh here. The transition region connects the refined contact region to the coarser hub region. The hub region is far from the contact and can use a coarser mesh. This graded mesh strategy preserves the tooth geometry while controlling the total number of elements. I use 10-node solid185-type second-order elements to ensure convergence and stability. The element size and layer details are given in Table 4.
| Region | Element size or layer | Purpose |
|---|---|---|
| Contact strain region | 1.5 mm, 5 refined layers | Resolve contact stress and subsurface strain |
| Transition region | Up to 10 mm | Gradual transition from fine to coarse mesh |
| Hub region | Up to 50 mm | Reduce total element count |
| Element type | 10-node solid185 second-order | Improve convergence and bending accuracy |
The mesh quality is evaluated using the Jacobian ratio, the aspect ratio, the warping factor, and the minimum angle. In my method, the contact strain region maintains a mesh quality above 0.9. This is important because the pinion gear and the ring gear contact problem is highly nonlinear. A high-quality mesh in the contact region reduces the risk of contact oscillation and improves the accuracy of the contact pressure distribution. The mesh is also checked for independence. I refine the mesh in the contact region until the maximum contact stress changes by less than a small tolerance. This ensures that the result is not controlled by the mesh size.
The general mesh preprocessing method can be applied to different types of large gears. It is not limited to the pinion gear or the ring gear of a specific machine. The key idea is to separate the contact strain region from the non-contact region and to use a mapped quadrilateral mesh that can be extruded into a hexahedral mesh. This avoids the poor-quality elements that often appear near the tooth root when a fully automatic tetrahedral mesh is used. It also avoids the excessive element count that would result from refining the entire gear. For the pinion gear, the refined contact region captures the local contact pressure. For the ring gear, the same method captures the contact pressure and the bending stress at the tooth root. The pinion gear tooth root is especially important because the pinion gear has fewer teeth and a higher number of load cycles than the ring gear.
5. Refined Finite Element Analysis of the Drive Gear Set
For the refined finite element analysis, I simplify the drive gear set to one pinion gear and a portion of the ring gear. This simplification is justified because the extreme working condition is a single-tooth-pair contact condition at the most heavily loaded position. The simplified model greatly reduces the computational cost while preserving the local contact and tooth-root behavior. I use a three-dimensional solid model. The pinion gear and the ring gear are modeled with the actual tooth profile, including the profile shift. The material is modeled as an elastic-plastic material with the yield limit of 1350 MPa and the allowable contact stress of 3000 MPa. The carburized layer is represented by the surface material properties and by the strain results in the subsurface region.
The contact between the pinion gear and the ring gear is defined as frictional contact. The pinion gear tooth surface is the contact surface, and the ring gear tooth surface is the target surface. I use the adjust-to-touch option to eliminate the initial gap. This option moves the contact nodes to the target surface if the initial gap is small, which improves convergence. The contact stiffness is updated at each iteration. The normal penalty stiffness is selected to balance penetration and convergence. The friction coefficient is set according to the lubrication condition. In the present analysis, I use a moderate friction coefficient to represent the mixed lubrication condition. The contact formulation is a surface-to-surface contact formulation, which is more accurate than a node-to-surface formulation for gear contact.
The boundary conditions are chosen to represent the actual support. I create a cylindrical coordinate system based on the inner bore surface of the pinion gear. I select the inner bore surface and apply a displacement constraint in this coordinate system. I release the rotational degree of freedom about the pinion gear axis and constrain the other two directions. This allows the pinion gear to rotate under the applied torque while preventing rigid-body motion. For the ring gear, I apply a displacement constraint on the outer surface and constrain all directions. This represents the support provided by the main bearing and the gearbox structure. I apply a moment on the inner bore surface of the pinion gear. The moment magnitude is taken from the extreme working condition. The torque direction is set so that the pinion gear drives the ring gear in the specified direction.
The load application is important for convergence. I apply the torque gradually using load steps. The first load step establishes the contact. The second load step applies the full torque. The third load step may be used for verification or for extracting the final results. I use the Newton-Raphson method with adaptive descent. The convergence criteria are based on force residuals and displacement increments. I also monitor the contact penetration and the contact pressure. If the contact pressure oscillates, I reduce the normal penalty stiffness or refine the contact mesh. If the force residual does not decrease, I check the contact status and the boundary conditions. The pinion gear and the ring gear contact problem is highly nonlinear, so a stable constraint and load setup is essential.
The finite element model uses the following material properties: elastic modulus \(E = 206000 \ \text{MPa}\), Poisson’s ratio \(\nu = 0.3\), yield strength \(\sigma_y = 1350 \ \text{MPa}\), and allowable contact stress \(\sigma_{HP} = 3000 \ \text{MPa}\). The material model is isotropic and homogeneous. The carburized layer is not modeled as a separate layer in the simplified model, but the strain results are interpreted with respect to the case depth. In a more detailed model, the case depth can be represented by a gradient material property. For the present bearing capacity analysis, the homogeneous model is sufficient because the maximum stress is below the yield limit and the contact stress is below the allowable contact stress.
6. Finite Element Results
The finite element simulation gives the contact stress, the tooth-root bending stress, and the strain distribution. The maximum contact stress is 1072.5 MPa. This is lower than the standard calculation value of 1171 MPa, but the difference is within ten percent. The contact region is approximately rectangular, and the maximum contact stress occurs in the central region of the contact area. This is consistent with the expected Hertzian contact behavior. The maximum contact stress is well below the allowable contact stress of 3000 MPa. Therefore, the pinion gear and the ring gear satisfy the contact strength requirement.
The maximum tooth-root bending stress of the pinion gear is 338.65 MPa. The maximum value is located at the center of the tooth-root transition fillet. The maximum tooth-root bending stress of the ring gear is 490.84 MPa. The maximum value is located at the center of the tooth-root chamfer region. Both values are well below the yield limit of 1350 MPa. Therefore, the pinion gear and the ring gear satisfy the bending strength requirement. The pinion gear has a lower tooth-root bending stress than the ring gear in this load case, but the pinion gear experiences more load cycles. This is an important consideration for fatigue life. The pinion gear tooth root should be monitored in the subsequent optimization stage.
| Result | Finite element value | Allowable | Status |
|---|---|---|---|
| Maximum contact stress | 1072.5 MPa | 3000 MPa | Satisfied |
| Pinion gear tooth-root bending stress | 338.65 MPa | 1350 MPa | Satisfied |
| Ring gear tooth-root bending stress | 490.84 MPa | 1350 MPa | Satisfied |
| Maximum strain | 0.0042 mm/mm | — | Within elastic range |
The strain simulation shows that the maximum strain is 0.0042 in the normalized sense. The maximum strain occurs inside the carburized layer of the ring gear tooth surface, at a certain distance from the contact surface. This is typical for gear contact: the maximum shear strain or equivalent strain can occur below the surface because of the combined normal and tangential contact loads. The tooth-root regions on both sides of the contacting teeth show stress concentration. The bending deformation region on the side away from the contact area is larger. This indicates that the tooth root is subjected to a bending moment that produces tension on one side and compression on the other. The pinion gear tooth root and the ring gear tooth root both require careful design of the fillet and the chamfer to reduce the stress concentration.
7. Comparison of Standard Calculation and Finite Element Simulation
I compare the standard calculation and the finite element simulation to validate the normative method. The maximum contact stress from the standard calculation is 1171 MPa, and the maximum contact stress from the finite element simulation is 1068 MPa or 1072.5 MPa depending on the extraction method. The relative difference is
$$ \Delta = \frac{|\sigma_{H,\text{std}} – \sigma_{H,\text{FE}}|}{\sigma_{H,\text{std}}} \times 100\%. $$
Using the representative finite element value of 1068 MPa, I obtain
$$ \Delta = \frac{|1171 – 1068|}{1171} \times 100\% = 8.80\%. $$
Using the value of 1072.5 MPa, I obtain
$$ \Delta = \frac{|1171 – 1072.5|}{1171} \times 100\% = 8.41\%. $$
Both deviations are within ten percent. This is a reasonable range for gear load-capacity analysis, especially for a large gear with complex support conditions. The difference is mainly due to the load distribution among the pinion gears, the contact stiffness, the mesh stiffness, and the local geometry. The standard calculation assumes a simplified load distribution. The finite element simulation resolves the local contact and the support stiffness. The agreement between the two methods gives confidence in the result.
| Quantity | Standard calculation | Finite element simulation | Deviation |
|---|---|---|---|
| Maximum contact stress | 1171 MPa | 1068 MPa | 8.80% |
| Pinion gear tooth-root bending stress | 374 MPa | 338.65 MPa | 9.46% |
| Ring gear tooth-root bending stress | 534 MPa | 490.84 MPa | 8.08% |
The comparison also shows that the finite element results are slightly lower than the standard calculation results. This is not a general rule, but it is reasonable in this case because the standard calculation uses conservative load distribution factors and does not account for the actual support flexibility. In other cases, the finite element result could be higher if the load concentration is more severe. Therefore, I do not use the finite element result to replace the standard calculation. I use both results together. The standard calculation provides the code-based safety margin, and the finite element simulation provides the local stress and strain distribution.
8. Parametric and Sensitivity Considerations
The bearing capacity of the pinion gear and the ring gear is sensitive to several parameters. I express the load-sharing factor for the \(j\)-th pinion gear as
$$ K_{\gamma,j} = \frac{T_{p,j}}{T_{p,\text{nom}}}, $$
where \(T_{p,j}\) is the actual torque carried by the \(j\)-th pinion gear and \(T_{p,\text{nom}}\) is the nominal torque. In an ideal system, \(K_{\gamma,j} = 1\) for all pinion gears. In a real system, the load-sharing factor depends on the manufacturing tolerance, the support stiffness, the ring gear flexibility, and the pinion gear position. The maximum load-sharing factor is
$$ K_{\gamma,\max} = \max_j K_{\gamma,j}. $$
The contact stress and the tooth-root bending stress scale with the load. For a linear elastic contact model, the contact stress scales approximately with the square root of the load, while the tooth-root bending stress scales approximately linearly with the load. Therefore,
$$ \sigma_H \propto \sqrt{F_t}, $$
$$ \sigma_F \propto F_t. $$
If the load-sharing factor increases by ten percent, the contact stress increases by approximately \(\sqrt{1.1} – 1 = 4.9\%\) and the tooth-root bending stress increases by approximately ten percent. This means that the tooth-root bending stress is more sensitive to load-sharing than the contact stress. For the pinion gear, which has a smaller tooth number and a higher number of load cycles, this sensitivity is important.
The mesh stiffness also influences the load distribution. The mesh stiffness per unit face width can be approximated as
$$ k_m = \frac{F_t}{b \delta}, $$
where \(\delta\) is the normal deflection. The total mesh stiffness is
$$ K_m = k_m b. $$
The load distribution along the face width depends on the mesh misalignment. If the misalignment is \(\Delta\), the load distribution factor can be approximated as
$$ K_{H\beta} = 1 + \frac{c \Delta}{F_t / b}, $$
where \(c\) is a stiffness coefficient. This equation shows that a larger face width or a smaller tangential force reduces the relative effect of misalignment. However, a larger face width also increases the sensitivity to misalignment because the absolute misalignment over the face width can be larger. This is one reason why large gears require careful alignment and support design.
The tooth-root stress concentration factor is defined as
$$ K_t = \frac{\sigma_{\max}}{\sigma_{\text{nom}}}, $$
where \(\sigma_{\max}\) is the maximum local stress and \(\sigma_{\text{nom}}\) is the nominal bending stress. The fillet radius, the profile shift, and the chamfer geometry influence \(K_t\). For the pinion gear, a larger fillet radius reduces \(K_t\) but may reduce the effective tooth thickness. For the ring gear, the chamfer and the root relief influence the stress distribution. The finite element method can resolve these local effects, which is why I use it in the normative method.
9. Fatigue and Life Implications
The bearing capacity analysis is not only a static strength check. The pinion gear and the ring gear are subjected to cyclic loads. The contact fatigue life and the bending fatigue life must be considered. The contact stress life equation can be written as
$$ N_H = N_{H0} \left( \frac{\sigma_{H0}}{\sigma_H} \right)^{m_H}, $$
where \(N_H\) is the number of cycles to contact fatigue, \(N_{H0}\) is the reference number of cycles, \(\sigma_{H0}\) is the reference contact stress, and \(m_H\) is the contact fatigue exponent. The bending fatigue life equation can be written as
$$ N_F = N_{F0} \left( \frac{\sigma_{F0}}{\sigma_F} \right)^{m_F}, $$
where \(N_F\) is the number of cycles to bending fatigue, \(N_{F0}\) is the reference number of cycles, \(\sigma_{F0}\) is the reference bending stress, and \(m_F\) is the bending fatigue exponent. For carburized gears, the contact fatigue exponent is often in the range of 6 to 10, and the bending fatigue exponent is often in the range of 6 to 12. These exponents depend on the material, heat treatment, and surface condition.
The cumulative damage can be estimated using Miner’s rule:
$$ D = \sum_{i=1}^{n} \frac{n_i}{N_i}, $$
where \(n_i\) is the number of cycles at stress level \(i\), \(N_i\) is the fatigue life at stress level \(i\), and \(D\) is the cumulative damage. The design is acceptable if \(D \leq 1\). For the pinion gear, the number of cycles is higher than that of the ring gear because the pinion gear rotates faster. If the pinion gear torque is \(T_p\) and the pinion gear speed is \(n_p\), the number of cycles over a time \(t\) is
$$ n = n_p t. $$
For the ring gear, the number of cycles is
$$ n = n_g t. $$
Because \(n_p = i n_g\), the pinion gear experiences \(i\) times as many cycles as the ring gear for the same operating time. This makes the pinion gear the critical component for fatigue life in many cases. The pinion gear tooth root should therefore be designed with a sufficient fillet radius and a suitable surface treatment. The pinion gear contact surface should also be designed to avoid pitting and micropitting.
10. Design Guidelines and Normative Procedure
Based on the results, I propose the following normative procedure for the bearing capacity analysis of a large boring machine drive gear set.
| Step | Action | Output |
|---|---|---|
| 1 | Define the gear geometry, material, and load cases | Design parameter table |
| 2 | Calculate the nominal pinion gear torque and tangential force | Load table |
| 3 | Build the complete drive gear set model | Assembly model |
| 4 | Apply constraints and loads according to the support conditions | Boundary condition set |
| 5 | Run the standard calculation | Contact stress and bending stress |
| 6 | Create the simplified finite element model | Refined solid model |
| 7 | Generate the general hexahedral mesh | High-quality mesh |
| 8 | Apply frictional contact and stable constraints | Converged contact model |
| 9 | Run the finite element simulation | Local stress and strain |
| 10 | Compare the standard and finite element results | Validation table |
| 11 | Evaluate fatigue life and optimization potential | Design recommendation |
I also recommend the following design checks. First, the maximum contact stress should be below the allowable contact stress with a safety factor greater than one. Second, the maximum tooth-root bending stress of the pinion gear and the ring gear should be below the yield limit with a safety factor greater than one. Third, the load-sharing factor among the pinion gears should be estimated and used to correct the nominal torque. Fourth, the mesh misalignment should be controlled by the support design and the assembly tolerance. Fifth, the fillet radius and the chamfer geometry should be optimized to reduce the tooth-root stress concentration. Sixth, the carburized case depth should be selected to ensure that the maximum subsurface strain occurs within the hardened layer. Seventh, the pinion gear should be checked for a higher number of load cycles than the ring gear.
11. Discussion of the Proposed Method
The proposed method has several advantages. It combines the speed of the standard calculation with the detail of the finite element simulation. The standard calculation provides a code-based reference and a conservative estimate. The finite element simulation provides the local contact pressure, the tooth-root stress, and the strain distribution. The general mesh preprocessing method solves the mesh quality problem for large gears. The stable constraint and load application method solves the convergence problem for the contact analysis. The method is repeatable and can be used for different large gear designs.
The method also has limitations. The simplified finite element model uses only one pinion gear and a portion of the ring gear. This is sufficient for the local contact and tooth-root stress, but it does not capture the global load-sharing among all pinion gears. The standard calculation uses the complete model, so the two methods complement each other. The material model is homogeneous and does not explicitly represent the carburized layer. This is acceptable for the bearing capacity check because the maximum stress is below the yield limit. For a more detailed fatigue analysis, a gradient material model could be used. The friction coefficient is assumed to be constant, but in reality it depends on the lubrication condition, the sliding speed, and the surface roughness. A sensitivity study could be performed to evaluate the effect of the friction coefficient.
Despite these limitations, the method provides a reliable engineering basis for the design of the pinion gear and the ring gear. It can be used to compare different design options, such as profile shift, face width, fillet radius, and surface treatment. It can also be used to evaluate the effect of load-sharing on the pinion gear stress. In my view, the most important contribution is the combination of a general mesh preprocessing method and a stable contact simulation procedure. Without these two elements, the finite element results for large gears can be unreliable. With them, the pinion gear and the ring gear can be analyzed with confidence.
12. Conclusion
I have presented a bearing capacity analysis method for the drive gear set of a large boring machine. The method combines engineering standard calculation with refined finite element simulation. I applied the method to a 6000-type drive gear set with 16 pinion gears and one large ring gear. The standard calculation gave a maximum contact stress of 1171 MPa, a pinion gear tooth-root bending stress of 374 MPa, and a ring gear tooth-root bending stress of 534 MPa. The finite element simulation gave a maximum contact stress of approximately 1068 MPa to 1072.5 MPa, a pinion gear tooth-root bending stress of 338.65 MPa, and a ring gear tooth-root bending stress of 490.84 MPa. The deviations between the two methods are within ten percent. The pinion gear and the ring gear satisfy the contact and bending strength requirements.
The general mesh preprocessing method produces a high-quality hexahedral mesh for large gears. It refines the contact strain region, transitions gradually to a coarser mesh, and maintains a mesh quality above 0.9 in the critical region. The stable constraint and load application method ensures convergence of the frictional contact analysis. The method can be applied to different large gear types and different load cases. It provides an important basis for the integrated design of the main bearing, the pinion gear, and the ring gear. It also provides a reference for tooth modification, root fillet optimization, and fatigue life improvement. I conclude that the proposed method is a practical and reliable normative approach for the bearing capacity analysis of large boring machine drive gear sets.