AC electric locomotives are widely used in modern railway transportation because of their high energy efficiency, precise traction control and relatively low maintenance cost. In the traction drive chain, the gearbox housing is one of the most critical structural components. It supports the transmission bearings, maintains the gear mesh alignment, and transmits the reaction forces generated by the traction torque to the bogie frame. Because the gearbox housing must survive a long service life under highly variable loads, its structural strength cannot be evaluated by simple static rules alone; fatigue assessment is equally important. In this study, I investigate a gearbox housing manufactured as a ductile iron casting for an AC electric locomotive. Using the finite element method, I evaluate the stress distribution under normal service conditions and under extreme short-circuit conditions, and I compare the results with the material limits and the Goodman fatigue criterion. The aim is to confirm whether the current ductile iron casting design has sufficient strength and fatigue safety margins for practical railway service.
Keywords: ductile iron casting; gearbox housing; finite element analysis; Goodman curve; strength verification.
The photograph below shows a typical ductile iron casting used for a traction gearbox housing. The cast structure integrates bearing seats, lubrication passages and mounting brackets into a single component. The image illustrates the complex geometry that must be analyzed in the structural strength study.
1. Introduction
With the continuous development of railway transportation toward higher speeds, higher axle loads and intelligent operation, the requirements for traction equipment have become more demanding. The AC drive system has replaced many DC drive systems because it offers better adhesion utilization, smaller maintenance effort and more flexible control. The traction motor torque is transmitted to the wheelset through a gearbox. In this transmission path, the gearbox housing is not only an enclosure but also a load-carrying structure. It must withstand the static weight of the gearbox, the dynamic inertia forces caused by track excitation, the reaction torque from the gears, and the vibrations generated by the gear mesh. In some situations, such as motor short-circuit, the gearbox can also experience a transient torque much larger than the normal starting torque. Therefore, the structural strength assessment must cover both fatigue loads and extreme overload conditions.
In the design stage of a traction gearbox, finite element analysis is an effective tool for predicting the stress state of the housing. The housing of the AC locomotive gearbox considered in this study is made of a ductile iron casting with the material grade EN-GJS-400-18-LT. This material has good castability, moderate strength and excellent ductility. It is commonly used in railway gearbox housings because it can absorb impact energy and resist crack propagation. However, the strength of the ductile iron casting is highly dependent on the local geometry, wall thickness, casting quality and stress concentrations. Thus, a detailed finite element model is necessary to identify potential weak zones and to verify that the design satisfies the required standard.
In this paper, I perform the strength analysis in two main steps. First, for the normal service condition, I use the Goodman criterion to assess fatigue strength. The Goodman criterion is suitable for cast materials under cyclic loading because it accounts for the effect of mean stress on the allowable stress amplitude. Second, for the abnormal motor short-circuit condition, I use the yield strength and the ultimate tensile strength to assess the static strength. The objective is to confirm that the stress at every critical point remains inside the acceptable domain defined by the relevant technical specification.
2. Gearbox structure and material
The gearbox analyzed in this study is a two-stage helical gear unit mounted on an AC electric locomotive. The gearbox uses a compact monolithic housing construction. The main components include the upper housing, lower housing, input shaft assembly, intermediate shaft assembly, output shaft assembly, bearings, oil pump, seals and suspension rods. The input shaft is connected to the traction motor through a flexible coupling and operates at motor speed. The output shaft is connected to the wheelset and carries the final traction torque. Because the space available on the locomotive bogie is limited, the gearbox is designed as an integrated ductile iron casting with bearing seats and lubrication galleries cast directly into the housing.
The gearbox parameters used in the finite element analysis are listed in Table 1. The rated power of the traction motor is 31 kW, the rated speed is 1150 r/min, and the rated torque is 257 N·m. The starting torque is 310 N·m, which is the highest torque in normal operation. The overall transmission ratio of the two-stage helical gear pair is 13.647.
| Parameter | Value |
|---|---|
| Transmission type | Two-stage helical cylindrical gear |
| Transmission ratio | 13.647 |
| Motor rated power | 31 kW |
| Motor starting torque | 310 N·m |
| Motor rated speed | 1150 r/min |
| Motor rated torque | 257 N·m |
The nominal relationship between the motor power, speed and torque can be expressed as:
\[
T = \frac{9550 P}{n}
\]
where \(T\) is the shaft torque in N·m, \(P\) is the power in kW, and \(n\) is the rotational speed in r/min. Substituting \(P = 31\) kW and \(n = 1150\) r/min gives \(T = 257\) N·m, which agrees with the rated torque listed in Table 1. The starting torque of 310 N·m is approximately 1.21 times the rated torque, while the abnormal short-circuit torque of 740 N·m is about 2.88 times the rated torque.
The material of the gearbox housing is a ductile iron casting grade EN-GJS-400-18-LT. The mechanical properties required by the material standard are listed in Table 2. The minimum tensile strength is 400 MPa, the minimum yield strength is 240 MPa, and the minimum elongation is 18%. In the finite element model, I used a Young’s modulus of 169 GPa, a Poisson’s ratio of 0.275 and a density of 7100 kg/m³, which are typical values for this grade of ductile iron casting.
| Property | Symbol | Minimum value |
|---|---|---|
| Tensile strength | \(R_m\) | 400 MPa |
| Yield strength | \(R_{p0.2}\) | 240 MPa |
| Elongation | \(A\) | 18% |
3. Finite element model
I built the finite element model using a three-dimensional solid model of the gearbox housing. The purpose of the finite element analysis was to compute the stress distribution under various load cases and to evaluate the structural strength of the ductile iron casting. To ensure both accuracy and computational efficiency, the geometric model was simplified by removing small features that do not significantly affect the global stress distribution. These features include small bolt holes, thread details, nameplate recesses and minor casting fillets. At the same time, the main structural details such as bearing bores, rib layouts, suspension rod seats and lubrication passages were retained because they influence the load paths and the local stress concentration behavior.
3.1 Contact definitions
The gearbox housing is connected to other components through bearing interfaces and bolted joints. In the finite element model, these interfaces were represented by appropriate contact and connection types. Table 3 summarizes the contact surface definitions used in the model. The bearing-to-housing contacts were modeled as revolute joints to allow rotation while constraining radial movement. The gear mesh contact between the input gear and the output gear was modeled as no-separation contact. All other connected parts were assumed to be perfectly bonded.
| Contact interface | Contact type |
|---|---|
| Housing and input bearing | Joint-Revolute |
| Input sleeve and input bearing | Joint-Revolute |
| Upper housing and intermediate shaft bearing | Joint-Revolute |
| Output sleeve and output bearing | Joint-Revolute |
| Lower housing and intermediate shaft bearing | Joint-Revolute |
| Input gear and output gear | No Separation |
| All other contacting surfaces | Bonded |
3.2 Mesh generation
The mesh was generated using a combination of tetrahedral and hexahedral elements. The complex cast regions of the ductile iron casting, such as the transition zones between the bearing bosses and the housing walls, were meshed with quadratic tetrahedral elements. The more regular regions, such as the bearing flanges and the side walls, were meshed with hexahedral elements to reduce the total element count. Local mesh refinement was applied at the bearing seats, the gear mesh zone and the suspension rod attachment points because these regions are expected to experience the highest stress gradients. The final mesh density was selected after a mesh sensitivity check, in which the maximum stress changed by less than 3% when the mesh was further refined.
3.3 Loads and boundary conditions
The boundary conditions represented the actual mounting of the gearbox on the locomotive bogie. The suspension rod seats were fixed, the output gear was constrained in all directions, and the bearing outer surfaces were fixed in the degrees of freedom corresponding to their connection to the housing. Figure 4 of the original technical report shows these positions; in my model, the same restraints were applied.
The external loads included the gravitational acceleration \(g = 9.8\) m/s², the motor torque applied at the input shaft, and the inertial accelerations applied to the whole assembly. The inertial load vector can be expressed as:
\[
\mathbf{F}_I = -m\,\mathbf{a}, \qquad \mathbf{a} = \left(A_x,\; A_y,\; A_z\right)^T
\]
where \(m\) is the mass of the gearbox body and \(\mathbf{a}\) is the acceleration vector acting on the component. The acceleration components \(A_x\), \(A_y\) and \(A_z\) were applied according to the standard IEC 61373—2010 for railway equipment. In the normal service condition, the acceleration components were \(64.3\) m/s² in the longitudinal direction, \(129\) m/s² in the lateral direction, and \(144\) m/s² in the vertical direction. The torque was applied in both positive and negative directions to cover the two possible rotation senses of the traction motor.
4. Load cases
4.1 Normal operating condition
In the normal operating condition, the starting torque is the largest torque that the gearbox can experience during normal traction. Therefore, the starting torque of 310 N·m was used for the fatigue strength assessment. The acceleration components were combined with the torque in eight different directions for positive torque and eight directions for negative torque, giving sixteen load cases in total. These load cases are listed in Table 4. All acceleration values are in m/s².
| Case | \(A_x\) (m/s²) | \(A_y\) (m/s²) | \(A_z\) (m/s²) | Torque (N·m) |
|---|---|---|---|---|
| ZZ1 | 64.3 | 129 | 144 | 310 |
| ZZ2 | 64.3 | 129 | -144 | 310 |
| ZZ3 | -64.3 | 129 | 144 | 310 |
| ZZ4 | -64.3 | 129 | -144 | 310 |
| ZZ5 | 64.3 | -129 | 144 | 310 |
| ZZ6 | 64.3 | -129 | -144 | 310 |
| ZZ7 | -64.3 | -129 | 144 | 310 |
| ZZ8 | -64.3 | -129 | -144 | 310 |
| ZF1 | 64.3 | 129 | 144 | -310 |
| ZF2 | 64.3 | 129 | -144 | -310 |
| ZF3 | -64.3 | 129 | 144 | -310 |
| ZF4 | -64.3 | 129 | -144 | -310 |
| ZF5 | 64.3 | -129 | 144 | -310 |
| ZF6 | 64.3 | -129 | -144 | -310 |
| ZF7 | -64.3 | -129 | 144 | -310 |
| ZF8 | -64.3 | -129 | -144 | -310 |
4.2 Abnormal operating condition
The abnormal operating condition corresponds to a motor short-circuit event. In this condition, the motor can produce a transient torque much larger than the normal starting torque. The maximum torque considered in this study is 740 N·m. The same acceleration combinations used in the normal condition were applied, but with the torque set to 740 N·m in both positive and negative directions. The sixteen abnormal load cases are listed in Table 5. For the abnormal condition, the strength is evaluated using the yield strength and the ultimate tensile strength of the ductile iron casting rather than the fatigue criterion.
| Case | \(A_x\) (m/s²) | \(A_y\) (m/s²) | \(A_z\) (m/s²) | Torque (N·m) |
|---|---|---|---|---|
| CZ1 | 1000 | 1000 | 1000 | 740 |
| CZ2 | 1000 | 1000 | -1000 | 740 |
| CZ3 | -1000 | 1000 | 1000 | 740 |
| CZ4 | -1000 | 1000 | -1000 | 740 |
| CZ5 | 1000 | -1000 | 1000 | 740 |
| CZ6 | 1000 | -1000 | -1000 | 740 |
| CZ7 | -1000 | -1000 | 1000 | 740 |
| CZ8 | -1000 | -1000 | -1000 | 740 |
| CF1 | 1000 | 1000 | 1000 | -740 |
| CF2 | 1000 | 1000 | -1000 | -740 |
| CF3 | -1000 | 1000 | 1000 | -740 |
| CF4 | -1000 | 1000 | -1000 | -740 |
| CF5 | 1000 | -1000 | 1000 | -740 |
| CF6 | 1000 | -1000 | -1000 | -740 |
| CF7 | -1000 | -1000 | 1000 | -740 |
| CF8 | -1000 | -1000 | -1000 | -740 |
5. Fatigue strength assessment under normal operating conditions
The fatigue assessment of the ductile iron casting under normal operating conditions was carried out using the Goodman criterion. In a cyclic stress state, the stress-time history at each node can be characterized by the maximum principal stress \(\sigma_{\max}\), the minimum principal stress \(\sigma_{\min}\), the mean stress \(\sigma_m\) and the alternating stress \(\sigma_a\). These quantities are defined as:
\[
\sigma_m = \frac{\sigma_{\max} + \sigma_{\min}}{2}, \qquad
\sigma_a = \frac{\sigma_{\max} – \sigma_{\min}}{2}.
\]
For the ductile iron casting material, the allowable combination of mean stress and alternating stress is limited by the Goodman curve. The Goodman condition can be written as:
\[
\frac{\sigma_a}{S_e} + \frac{\sigma_m}{R_m} \le 1,
\]
where \(S_e\) is the fatigue limit of the material under fully reversed loading and \(R_m\) is the ultimate tensile strength. Since the service load cases include both positive and negative torque directions and different acceleration directions, each node of the ductile iron casting is subjected to a multiaxial stress cycle. To apply the Goodman criterion in a practical way, the stress cycle at each node is represented by the pair of mean stress and limiting principal stress values. The standard TJ/JW 064—2015 defines the acceptable region in the Goodman diagram using the control points listed in Table 6.
| Point | Mean stress \(\sigma_m\) (MPa) | Limiting principal stress \(\sigma_{lim}\) (MPa) | Boundary type |
|---|---|---|---|
| A | 0 | 183 | Maximum stress limit |
| B | 105 | 240 | Maximum stress limit |
| C | 240 | 240 | Tensile static limit |
| D | 105 | -30 | Minimum stress limit |
| E | 0 | -183 | Minimum stress limit |
| F | -57 | -240 | Minimum stress limit |
| G | -240 | -240 | Compressive static limit |
| H | -57 | 126 | Maximum stress limit |
During the fatigue assessment, I extracted the maximum and minimum principal stresses at every node of the gearbox housing for each of the sixteen normal load cases listed in Table 4. For each node, I then computed the mean stress \(\sigma_m\) and formed two coordinate pairs: \((\sigma_m, \sigma_{\max})\) and \((\sigma_m, \sigma_{\min})\). These coordinate pairs were plotted on the Goodman diagram. A node is considered safe when both coordinate pairs lie inside the allowable polygon defined by the control points in Table 6. If any pair falls outside the polygon, the fatigue criterion is not satisfied at that node.
After processing all the finite element results, I found that every node of the ductile iron casting under all sixteen normal load cases produced coordinate pairs inside the Goodman polygon. No node exceeded the allowable stress range. This means that the gearbox housing satisfies the fatigue strength requirement for normal service conditions and has sufficient resistance to the alternating loads caused by track irregularities and traction torque variations.
6. Static strength assessment under abnormal operating conditions
For the abnormal operating condition, the motor short-circuit torque of 740 N·m is much larger than the normal torque. Under such a load, the ductile iron casting is not expected to operate for a long period; instead, the design only needs to prevent plastic collapse and rupture. Therefore, I evaluated the stress using the von Mises equivalent stress and the maximum principal stress. The von Mises equivalent stress is calculated from the principal stresses \(\sigma_1\), \(\sigma_2\) and \(\sigma_3\) as follows:
\[
\sigma_{vM} = \sqrt{
\frac{1}{2}
\left[
(\sigma_1 – \sigma_2)^2 +
(\sigma_2 – \sigma_3)^2 +
(\sigma_3 – \sigma_1)^2
\right]
}.
\]
The maximum principal stress was also extracted because brittle cast materials can be sensitive to tensile principal stress. The finite element simulation was performed for all sixteen abnormal load cases listed in Table 5. The maximum von Mises stress and the maximum principal stress obtained for each load case are summarized in Table 7.
| Case | Maximum von Mises stress (MPa) | Maximum principal stress (MPa) |
|---|---|---|
| CZ1 | 101.8 | 94.1 |
| CZ2 | 105.0 | 141.0 |
| CZ3 | 91.9 | 101.4 |
| CZ4 | 106.7 | 127.5 |
| CZ5 | 110.8 | 77.5 |
| CZ6 | 92.3 | 106.1 |
| CZ7 | 107.5 | 113.2 |
| CZ8 | 98.7 | 111.4 |
| CF1 | 101.8 | 94.1 |
| CF2 | 105.0 | 141.0 |
| CF3 | 91.9 | 101.4 |
| CF4 | 106.7 | 127.5 |
| CF5 | 110.8 | 77.5 |
| CF6 | 92.3 | 106.1 |
| CF7 | 107.5 | 113.2 |
| CF8 | 98.7 | 111.4 |
The results show that the maximum von Mises stress in all abnormal load cases is 110.8 MPa, which occurs in load cases CZ5 and CF5. The maximum principal stress is 141.0 MPa, which occurs in load cases CZ2 and CF2. The allowable yield strength of the ductile iron casting is 240 MPa, and the ultimate tensile strength is 400 MPa. The static safety factors can be defined as:
\[
n_y = \frac{R_{p0.2}}{\sigma_{vM}^{max}}, \qquad
n_u = \frac{R_m}{\sigma_{1}^{max}}.
\]
Using the maximum values from the finite element analysis, the yield safety factor is:
\[
n_y = \frac{240}{110.8} = 2.17.
\]
The tensile safety factor based on the maximum principal stress is:
\[
n_u = \frac{400}{141.0} = 2.84.
\]
Both safety factors are greater than unity, which means the gearbox housing has sufficient static strength to withstand the motor short-circuit condition without global yielding or tensile rupture. The ductile iron casting design therefore satisfies the strength requirement for abnormal operating conditions.
7. Discussion of stress concentration
Although the overall strength of the ductile iron casting is acceptable, the finite element results reveal several regions with higher local stress. In load case CZ2, the maximum principal stress of 141 MPa occurs at the connection between the gearbox housing and the suspension rod. This location is directly loaded by the suspension reaction forces and also experiences a high bending moment when the vertical acceleration is negative. The stress concentration is caused by the abrupt change in cross-section near the suspension rod seat and by the local wall thickness transition.
In load case CZ5, the maximum von Mises stress of 110.8 MPa appears at the intermediate gear bearing seat. This region transfers the gear mesh forces from the intermediate shaft to the housing. The bearing bore is subjected to a concentrated radial load, and the local geometry creates a stress concentration at the edge of the bore. The input shaft bearing lower side also shows elevated stress in several load cases, especially when the torque direction and the acceleration direction produce a superimposed bending effect.
These locations should be considered as the most critical regions of the ductile iron casting. Although the computed stresses are still well below the material limits, local improvements may be beneficial for increasing the fatigue life and for providing an additional margin in case of unexpected overloads. Possible design modifications include increasing the fillet radius at the transition zones, adding local ribs around the bearing seats, and smoothing the wall thickness transitions near the suspension rod attachment.
From a manufacturing point of view, the ductile iron casting should be examined carefully for casting defects such as shrinkage porosity, sand inclusions and cold shuts in these critical areas. Even though the finite element analysis assumes a perfectly sound casting, real castings may contain small discontinuities that reduce the local fatigue strength. Therefore, non-destructive testing of the high-stress zones is recommended during production.
8. Conclusions
In this study, I performed a detailed finite element strength analysis of a gearbox housing made of a ductile iron casting for an AC electric locomotive. The analysis covered both normal service conditions and abnormal motor short-circuit conditions. The following conclusions can be drawn from the results.
First, under normal operating conditions, the gearbox housing satisfies the Goodman fatigue criterion. All nodes of the ductile iron casting produced stress coordinate pairs located inside the allowable Goodman envelope for all sixteen normal load cases. This confirms that the housing has adequate fatigue strength for long-term service under alternating traction torque and track-induced inertial loads.
Second, under abnormal operating conditions with a motor short-circuit torque of 740 N·m, the maximum von Mises stress is 110.8 MPa and the maximum principal stress is 141.0 MPa. These values remain below the material yield strength of 240 MPa and the ultimate tensile strength of 400 MPa. The yield safety factor is 2.17 and the tensile safety factor is 2.84, indicating that the ductile iron casting has sufficient static strength to prevent structural failure during an abnormal overload event.
Third, the stress concentration analysis shows that the most critical regions are the suspension rod attachment, the intermediate gear bearing seat and the lower side of the input shaft bearing. Although these regions do not exceed the allowable stress limits, they are the most likely locations for fatigue crack initiation under prolonged service. Future design optimization should focus on increasing transition fillets, adding local ribs and improving the casting quality in these areas.
Fourth, the overall design of the gearbox housing as a ductile iron casting satisfies the requirements of the relevant standards and has a good margin of safety. The finite element modeling process and the strength evaluation method described in this study can be used as a reference for the structural optimization and series development of traction gearbox housings in future locomotive projects.