Fatigue Life Simulation of Rack and Pinion Gear in Welding Robots

In modern industrial automation, the rack and pinion gear mechanism plays a critical role in converting rotational motion into linear displacement, particularly in welding robots where precision and reliability are paramount. As a key component of the gantry assembly, the rack and pinion system ensures smooth and accurate movement along the Y-axis. However, prolonged operation under cyclic loads often leads to surface damage, wear, or even tooth fracture in the rack and pinion gear, resulting in transmission instability or failure. This study focuses on simulating the fatigue life of a rack and pinion gear system under dynamic conditions using finite element analysis (FEA) and nCode DesignLife software. By investigating the effects of various parameters such as friction coefficient, rotational speed, and load, I aim to provide insights for optimizing the rack and pinion design to enhance durability and performance in welding applications.

The rack and pinion gear mechanism in this study is designed as a spur gear system to facilitate easy manufacturing and maintenance. The material selected for both the rack and pinion is 42CrMo, known for its high strength and fatigue resistance. Key material properties include a density of 7,850 kg/m³, elastic modulus of 2.12 × 10¹¹ Pa, Poisson’s ratio of 0.28, yield strength of 930 MPa, and ultimate tensile strength of 1,080 MPa. The design parameters, summarized in Table 1, were calculated based on the maximum operational requirements of the welding robot, including a module of 2 mm, 23 teeth on the pinion, a pressure angle of 20°, and a pinion shaft diameter of 22 mm. Using SolidWorks, a parametric 3D model of the rack and pinion gear was developed, ensuring accurate representation of the meshing behavior.

To streamline the finite element analysis, the rack and pinion model was simplified while preserving simulation accuracy. Given that the spur gear system has a contact ratio of approximately 1.981, two pairs of teeth were retained in the model to maintain the required overlap during meshing. The mesh was generated using tetrahedral elements with a size of 0.5 mm, ensuring a mesh quality above 0.6 for reliable nonlinear contact analysis. The finite element model was then imported into ANSYS Workbench for transient dynamics simulation, where boundary conditions were applied to replicate real-world operational scenarios. The rack was fixed to the frame, while the pinion was allowed to rotate on the rack with a defined friction coefficient. A rotational velocity of 2.65 rad/s (equivalent to 200 RPM after gear reduction) was applied to the pinion, and a load of 2,190 N was imposed in the opposite direction to simulate the resistance from the gantry assembly’s weight and friction in linear guides. This load was derived using the formula:

$$ f = K F_N f_s $$

where $K = 0.15$ is the friction coefficient of the linear guides, $F_N = 3,650$ N is the gravitational force from the 365 kg gantry, and $f_s = 8$ is the static safety factor accounting for impact loads during startup.

Table 1: Design Parameters of the Rack and Pinion Gear
Parameter Value
Module 2 mm
Number of Pinion Teeth 23
Pressure Angle 20°
Pinion Width 25 mm
Material 42CrMo

Transient dynamics analysis was conducted over a time period of 0.206 seconds to capture a complete meshing cycle of the rack and pinion gear. The analysis revealed a maximum penetration of 3.2984 × 10⁻⁴ mm at the pitch circle of the pinion, indicating negligible risk of tooth gouging or adhesion. The equivalent stress distribution showed that the highest stresses occurred at the tooth root and meshing regions, with peak values observed during the startup phase due to impact loads. For instance, at 0.01648 seconds, the contact stress reached 234.42 MPa, while during stable operation, periodic peaks such as 322.31 MPa at 0.07004 seconds and 279.96 MPa at 0.16230 seconds were recorded. All stress values remained well below the material’s yield strength of 930 MPa, confirming that the rack and pinion gear meets strength requirements and is susceptible to high-cycle fatigue failure under repeated loading.

To evaluate the fatigue life of the rack and pinion gear, I employed the nCode DesignLife software, integrating the transient dynamics results from ANSYS. The fatigue analysis framework was built on the SN TimeSeries module, which converts stress tensors into scalar values for cycle counting and applies the S-N curve of the material to predict life. Since 42CrMo was not available in the nCode material library, I created a custom S-N curve, adjusted for ultimate tensile strength (UTS) corrections, as shown in Figure 1. The Goodman method was used for mean stress correction, and the Palmgren-Miner linear cumulative damage theory was applied to calculate fatigue damage. According to this theory, the total damage $D$ is given by:

$$ D = \sum_{i=1}^{K} \frac{n_i}{N_i} $$

where $n_i$ is the number of cycles at stress level $i$, and $N_i$ is the fatigue life at that stress level. For nonlinear damage accumulation, the Corten-Dolan theory was considered, expressed as:

$$ D = \sum_{i=1}^{p} n_i m_i^c r_i^d $$

where $m_i$ is the number of damage nuclei, $r_i$ is the damage propagation rate, and $c$ and $d$ are material constants. However, for this study, the linear approach was deemed sufficient due to the high-cycle nature of the fatigue.

The fatigue life results indicated a minimum life of 2.134 × 10⁷ cycles at node 40,082, located near the pitch circle where stress concentration is highest. The fatigue damage cloud plot showed maximum damage values of 4.687 × 10⁻⁸, primarily in the meshing regions due to prolonged sliding and stress accumulation. This underscores the importance of surface treatments, such as heat treatment, to reduce friction and stress concentrations in the rack and pinion gear.

Table 2: Effects of Input Parameters on Rack and Pinion Gear Performance
Parameter Range Effect on Equivalent Stress (MPa) Effect on Fatigue Life (Cycles) Effect on Fatigue Damage
Friction Coefficient 0.1–0.8 Increase from 150 to 350 MPa Decrease to 1.31 × 10⁴ Increase to 6.86 × 10⁻⁶
Rotational Speed (rad/s) 1.5–3.5 Moderate increase Significant decrease to 2.85 × 10⁵ Moderate increase
Load (N) 1,000–3,000 Linear increase (~62.11 MPa per 500 N) Decrease to 2.29 × 10⁶ Minor increase

To further investigate the dynamic behavior of the rack and pinion gear, I analyzed the influence of varying input parameters on fatigue performance. As summarized in Table 2, the friction coefficient between the rack and pinion teeth had the most pronounced effect on equivalent stress, with a sharp increase from 150 MPa to 350 MPa as the coefficient rose from 0.1 to 0.8. This is attributed to increased frictional forces and wear during meshing. Rotational speed exhibited the greatest impact on fatigue cycle counts, reducing life to 2.85 × 10⁵ cycles at 3.5 rad/s due to higher meshing frequencies and impact loads. Load variations showed a linear relationship with stress, with each 500 N increment raising stress by approximately 62.11 MPa, but the effect on fatigue damage was relatively minor compared to friction and speed. These findings highlight that optimizing the rack and pinion gear for fatigue resistance should prioritize reducing friction through surface finishing and controlling operational speeds.

The transient dynamics analysis also involved evaluating the contact stress over time using the time integration method. The stress history curve displayed periodic fluctuations, with peaks corresponding to tooth engagement events. The maximum von Mises stress of 322.31 MPa occurred at 0.07004 seconds, validating the robustness of the rack and pinion design under extreme conditions. The structural integrity was further confirmed by the low penetration values, ensuring no adhesive wear occurs in the rack and pinion mechanism.

In terms of fatigue life prediction, the S-N curve for 42CrMo was fitted using the following power-law relationship:

$$ S = A N^B $$

where $S$ is the stress amplitude, $N$ is the number of cycles to failure, and $A$ and $B$ are material constants derived from experimental data. For 42CrMo, the curve parameters were calibrated to reflect the high-cycle fatigue regime, with stress levels decreasing exponentially as cycle counts increase. This model was essential for accurate life estimation in nCode, especially when dealing with the variable amplitude loading experienced by the rack and pinion gear in welding applications.

Additionally, I explored the sensitivity of the rack and pinion system to different loading scenarios by simulating multiple load cases. For instance, under a constant load of 2,000 N and a friction coefficient of 0.15, the fatigue life was computed to be over 10⁷ cycles, indicating good durability. However, when dynamic loads with fluctuations were introduced, the life reduced significantly due to the accumulation of damage at stress concentrations. This emphasizes the need for comprehensive load spectrum analysis in real-world rack and pinion applications to prevent unforeseen failures.

In conclusion, this study demonstrates the effectiveness of combining FEA and fatigue analysis tools to assess the performance of rack and pinion gear systems in welding robots. The rack and pinion mechanism was found to satisfy strength criteria under transient loads, with fatigue life primarily limited by surface interactions and operational parameters. To enhance the longevity of rack and pinion gears, future work should focus on material enhancements, such as carburizing or nitriding to improve surface hardness, and operational adjustments to minimize friction and speed variations. The insights gained from this simulation provide a foundation for anti-fatigue optimization in rack and pinion designs, ensuring reliable and efficient operation in industrial automation.

Scroll to Top