In the mining industry, belt conveyors serve as critical transportation equipment, and the reducer is a pivotal component within the transmission system. The performance of the reducer ensures the reliable operation of the entire conveyor. The input shaft, particularly a bevel gear shaft, plays a vital role in connecting the motor and the load, transmitting torque while withstanding various forces. However, frequent failures and fractures at the shaft extension of the input bevel gear shaft in reducers have been observed in practice, leading to operational downtime and safety concerns. These failures often occur under dynamic conditions, such as during startup, shutdown, or unexpected loading scenarios where impact loads are prevalent. While static analyses are common, the influence of transient impact loads on the structural integrity of the bevel gear shaft remains underexplored. Therefore, this article employs finite element analysis (FEA) via ANSYS Workbench to conduct a comprehensive investigation. Initially, a static structural analysis is performed. Subsequently, the dynamic response under impact loading is examined to assess stress and strain characteristics. Finally, a modal analysis is conducted to evaluate vibration behavior and resonance risks. The findings aim to provide a theoretical foundation for optimizing the design of such bevel gear shafts, enhancing their durability and reliability in harsh mining environments.
The bevel gear shaft in question is part of a reducer used in a belt conveyor drive system. The motor transmits power through a splined coupling to the reducer’s input shaft, which then drives the conveyor via speed reduction and torque multiplication. The bevel gear shaft is subjected not only to transmitted torque but also to unbalanced forces, such as the gravitational load of the suspended coupling, which can induce bending moments. To facilitate an efficient and accurate finite element simulation, certain geometric features of the bevel gear shaft are simplified in the three-dimensional model. These simplifications include omitting small fillets, grooves, detailed tooth profiles of the bevel gear, and threads, as their influence on global stress distribution is considered secondary for this study. The primary focus remains on the shaft body, keyways, and critical shoulder regions. The simplified geometry is essential for meshing and computational economy while retaining the essential load-bearing characteristics of the bevel gear shaft.
The three-dimensional model of the bevel gear shaft is created using UG NX software. It is then imported into ANSYS Workbench in IGES format. The material assigned is 17CrNiMo6 alloy steel, a high-strength carburizing steel commonly used for gears in heavy-duty applications like mining and transportation. Its properties include high surface hardness, excellent wear resistance, and good core toughness. The material properties used in the simulation are summarized in Table 1.
| Property | Value | Unit |
|---|---|---|
| Density (ρ) | 7850 | kg/m³ |
| Young’s Modulus (E) | 2.1e11 | Pa |
| Poisson’s Ratio (ν) | 0.3 | – |
| Yield Strength (σ_y) | 585e6 | Pa |
| Tensile Strength (σ_u) | 980e6 | Pa |
For meshing, the Solid186 element is selected. This is a higher-order 3D 20-node solid element with three degrees of freedom per node (translations in x, y, z directions), well-suited for modeling complex geometries and stress gradients. A mesh independence study is conducted to ensure solution accuracy without excessive computational cost. The final mesh, shown in Figure 2 of the original context (though not referenced directly here), consists of approximately 377,383 elements and 527,505 nodes. This refined mesh is particularly dense in areas of expected stress concentration, such as shaft shoulders and keyway regions.
The boundary conditions and loads applied simulate the actual operating environment of the bevel gear shaft. Gravity is applied as a standard earth acceleration of 9.8066 m/s² in the Y-direction. A rotational velocity of 1480 rpm (approximately 155 rad/s) is applied about the Z-axis, representing the input speed from the motor. Constraints are applied to simulate bearing supports: cylindrical supports are applied on the cylindrical surfaces where bearings are mounted, restricting radial and axial displacements as appropriate. The shoulder face and the face where the locking nut sits are constrained against displacement in the Z-direction. The pitch circle surface of the bevel gear is constrained to have zero rotation about the Z-axis, simulating the reaction from the meshing gear. The input torque is converted into a pressure load and applied on the side faces of the key. The weight of the coupling is modeled as a remote force acting at the coupling’s geometric center, connected to the shaft end via a remote point. The applied operational loads are detailed in Table 2.
| Parameter | Value | Unit |
|---|---|---|
| Input Shaft Diameter | 70 | mm |
| Input Rotational Speed (ω) | 1480 | rpm |
| Input Power (P) | 238 | kW |
| Rated Input Torque (T_rated) | 1535.7 | Nm |
| Remote Force (Coupling Weight) | 1569 | N |
The static structural analysis provides the baseline stress and deformation under steady-state operating conditions. The governing equilibrium equation solved is:
$$ [K]\{u\} = \{F\} $$
where [K] is the global stiffness matrix, {u} is the nodal displacement vector, and {F} is the applied force vector. The results indicate that the maximum Von Mises stress is 302 MPa, located at a shaft shoulder, which is below the material’s yield strength of 585 MPa. The maximum shear stress is 157.9 MPa at the same location. The total deformation has a maximum value of 0.31 mm, which is within acceptable limits for such applications. The stress concentration factor at the shoulder can be theoretically approximated for a stepped shaft under bending and torsion. The nominal stress components need evaluation. For torsion, the shear stress (τ) is given by:
$$ \tau = \frac{T \cdot r}{J} $$
where T is the torque, r is the radius, and J is the polar moment of inertia. For the shaft diameter of 70 mm at the input, J is:
$$ J = \frac{\pi d^4}{32} = \frac{\pi (0.07)^4}{32} \approx 2.357 \times 10^{-6} \, m^4 $$
The shear stress due to torque is then:
$$ \tau_{torque} = \frac{1535.7 \times 0.035}{2.357 \times 10^{-6}} \approx 22.8 \times 10^6 \, Pa = 22.8 \, MPa $$
Bending stresses arise from the remote force (coupling weight) and possibly misalignment. The bending moment (M) at the shoulder can be significant. The bending stress (σ_b) is:
$$ \sigma_b = \frac{M \cdot c}{I} $$
where c is the distance from the neutral axis (d/2), and I is the area moment of inertia (I = πd⁴/64). The Von Mises stress (σ_vm) under combined stress states is calculated as:
$$ \sigma_{vm} = \sqrt{\sigma_x^2 + \sigma_y^2 + \sigma_z^2 – \sigma_x\sigma_y – \sigma_y\sigma_z – \sigma_z\sigma_x + 3(\tau_{xy}^2 + \tau_{yz}^2 + \tau_{zx}^2)} $$
For a shaft primarily under torsion and bending, this simplifies to:
$$ \sigma_{vm} = \sqrt{\sigma_b^2 + 3\tau_{torque}^2} $$
The static FEA results confirm that the bevel gear shaft is safe under normal loading, but the stress concentration at the shoulder is notable.
The dynamic analysis under impact loading is crucial, as real-world operations involve transients. The impact load is modeled as a torque pulse with a magnitude twice the rated torque (3071.4 Nm). The time history of the applied torque is defined over 0.5 seconds: a rapid increase to 3071.4 Nm from 0 to 0.1 s, a drop back to rated torque from 0.1 to 0.2 s, a constant rated torque from 0.2 to 0.3 s, and a repetition of the pulse from 0.3 to 0.5 s. This simulates severe start-up or sudden jamming conditions. The equation of motion for the transient structural analysis is:
$$ [M]\{\ddot{u}\} + [C]\{\dot{u}\} + [K]\{u\} = \{F(t)\} $$
where [M] is the mass matrix, [C] is the damping matrix (assumed proportional for simplicity), and {F(t)} is the time-varying load vector. The Newmark time integration method is typically used in Workbench for such dynamic problems.
The results reveal a dramatic increase in stress. The maximum Von Mises stress follows the load profile, peaking at 650.41 MPa at 0.1 s and 0.4 s. This exceeds the yield strength of 585 MPa, indicating that plastic deformation or fracture initiation is likely at the shaft shoulder under such impact conditions. The maximum shear stress reaches 363.57 MPa. The time-history data for stress at the critical point is summarized in Table 3.
| Time (s) | Applied Torque (Nm) | Max Von Mises Stress (MPa) | Max Shear Stress (MPa) |
|---|---|---|---|
| 0.0 – 0.1 | 0 → 3071.4 | Increasing to 650.41 | Increasing to 363.57 |
| 0.1 – 0.2 | 3071.4 → 1535.7 | Decreasing to ~350 | Decreasing to ~200 |
| 0.2 – 0.3 | 1535.7 (constant) | ~350.6 | ~195 |
| 0.3 – 0.5 | Repeat of 0-0.2 s | Peak at 650.41 again | Peak at 363.57 again |
The stress concentration factor (K_t) for the shoulder can be estimated from empirical charts based on geometry (fillet radius, diameter ratio). The high dynamic stress underscores the vulnerability of the bevel gear shaft to transient overloads. Repeated such events would lead to fatigue failure, even if single events cause yield.
Modal analysis is performed to determine the natural frequencies and mode shapes of the bevel gear shaft. This is essential to avoid resonance, where the operating frequency coincides with a natural frequency, leading to excessive vibrations and accelerated fatigue. The free vibration equation is:
$$ ([K] – \omega_i^2 [M]) \{\phi_i\} = 0 $$
where ω_i is the i-th natural frequency (rad/s), and {φ_i} is the corresponding mode shape vector. The Lanczos method is used to extract the first four modes. The results are listed in Table 4, and the mode shapes are described qualitatively.
| Mode Number | Natural Frequency (Hz) | Description of Mode Shape |
|---|---|---|
| 1 | 1132.6 | Lateral translation (bending) of the shaft extension in the Y-direction. |
| 2 | 1736.7 | Lateral translation (bending) of the shaft extension in the X-direction. |
| 3 | 2241.3 | Torsional vibration of the shaft extension about the Z-axis. |
| 4 | 2718.4 | Global torsional vibration of the entire shaft about the X-axis. |
The operating frequency (f_op) is determined from the input speed:
$$ f_{op} = \frac{1480 \, \text{rpm}}{60} \approx 24.67 \, \text{Hz} $$
This is far below the first natural frequency (1132.6 Hz), so resonance at the operating speed is not a concern. However, excitation sources might include gear meshing frequencies. The meshing frequency (f_mesh) for the bevel gear depends on the number of teeth (Z) and rotational speed. If we assume a typical tooth count, for example Z=20, then:
$$ f_{mesh} = Z \times f_{op} = 20 \times 24.67 \approx 493.4 \, \text{Hz} $$
This is still below the first natural frequency. However, higher harmonics or external impacts could excite these modes. The analysis shows that the shaft extension is the most flexible part, exhibiting the largest displacements in the first two bending modes, making it susceptible to fatigue if excited at or near its natural frequencies.
Further discussion on the design implications is warranted. The stress concentration at the shoulder is a critical design flaw. The theoretical stress concentration factor for bending and torsion in a shaft with a shoulder fillet can be expressed using empirical formulas. For example, for a shaft in torsion, the stress concentration factor K_ts is a function of the ratio of fillet radius (r) to smaller diameter (d), and the ratio of the two diameters (D/d). Increasing the fillet radius is a standard method to reduce K_t. The design of the bevel gear shaft must incorporate generous fillets at all shoulders. Additionally, surface treatments like shot peening can introduce compressive residual stresses, improving fatigue life. The material choice, 17CrNiMo6, is appropriate, but its heat treatment process (carburizing depth, tempering) must be optimized to ensure a tough core and a hard, wear-resistant surface, especially for the bevel gear teeth.
Regarding dynamic loads, the system should incorporate soft-start mechanisms or torque limiters to mitigate impact loads during startup or overloads. The analysis clearly shows that even short-duration overloads can induce stresses beyond yield, leading to plastic deformation and crack initiation. A factor of safety based on dynamic loading should be applied. The factor of safety (FS) under impact can be defined as the ratio of material yield strength to the maximum dynamic Von Mises stress:
$$ FS_{impact} = \frac{\sigma_y}{\sigma_{vm,max}} = \frac{585}{650.41} \approx 0.9 $$
This value less than 1 confirms the risk of failure under the simulated impact. For safe design, the factor of safety should typically be above 1.5 or higher for dynamic loading.
For vibration control, the modal analysis indicates that the bevel gear shaft’s natural frequencies are high, but the mode shapes highlight the shaft extension as a vulnerable area. Stiffening the shaft extension or modifying the support stiffness (bearing preload, housing design) can shift these frequencies. Damping elements, such as elastomeric couplings, can also help absorb vibrations and reduce transient load transmission to the bevel gear shaft.
In conclusion, this comprehensive finite element analysis of a coal mine reducer bevel gear shaft incorporating static, dynamic, and modal studies provides valuable insights. The static analysis confirms adequate strength under normal loads, but identifies stress concentrations at geometric discontinuities. The dynamic analysis under impact loading reveals that transient overloads can induce stresses exceeding the material yield limit, posing a high risk of fracture at the shoulder regions. This underscores the importance of considering impact scenarios in the design phase. The modal analysis shows that the shaft’s natural frequencies are sufficiently separated from the operating frequency, but the shaft extension is the most flexible part and prone to resonant vibrations if excited. Therefore, design optimizations should focus on reducing stress concentrations through improved geometry, selecting appropriate materials and treatments, implementing measures to mitigate impact loads, and considering vibration damping. These measures will enhance the reliability and service life of bevel gear shafts in demanding mining applications. Future work could involve fatigue life prediction based on the dynamic stress history and experimental validation of the FEA models.