The manufacturing of large-scale, structurally intricate components for heavy machinery presents a significant challenge within the domain of steel casting. Components such as the rocker arm housing for shearers, which feature variable cross-sections and multi-stage wall thicknesses, are particularly prone to defects like shrinkage porosity, hot tearing, and residual stress concentration if the casting process is not meticulously designed and controlled. Traditional trial-and-error methods for process development are costly and time-intensive, especially for such massive castings. This work details a comprehensive numerical investigation and subsequent optimization of the casting process for a ZG20SiMn steel rocker arm shell, leveraging simulation technology to enhance quality and reliability.
The core objective was to diagnose and mitigate casting defects through a virtual prototyping approach. Two primary gating system designs—a top-gating scheme and a bottom-gating scheme—were conceptualized for the housing. The entire filling, solidification, and cooling process was simulated using a dedicated casting simulation software. By analyzing the resulting temperature fields, solidification sequences, and predicted defect locations, the more promising process was identified and subsequently refined. The optimization was guided by quantitative criteria, including the Niyama criterion for shrinkage and the distribution of thermal stress, culminating in a validated process that significantly reduces defect propensity and improves the structural integrity of the final steel casting.
Foundational Methodology and Numerical Modeling
The fidelity of a casting simulation hinges on the accurate mathematical representation of the underlying physical phenomena. For this steel casting analysis, the governing equations for fluid flow, heat transfer, and stress development were solved numerically.
Governing Equations for Filling and Solidification
The flow of molten ZG20SiMn steel during mold filling is treated as a transient, incompressible flow. The conservation of mass (continuity equation) is expressed as:
$$
\frac{\partial (\rho u_x)}{\partial x} + \frac{\partial (\rho u_y)}{\partial y} + \frac{\partial (\rho u_z)}{\partial z} + \frac{\partial \rho}{\partial t} = 0
$$
For incompressible flow, $\frac{\partial \rho}{\partial t} = 0$, simplifying to:
$$
\frac{\partial (u_x)}{\partial x} + \frac{\partial (u_y)}{\partial y} + \frac{\partial (u_z)}{\partial z} = 0
$$
where $u_x$, $u_y$, and $u_z$ are the velocity components, and $\rho$ is the fluid density. The energy equation, accounting for heat transfer during filling and solidification, is given by:
$$
\frac{dT}{dt} + u\frac{\partial T}{\partial x} + v\frac{\partial T}{\partial y} + w\frac{\partial T}{\partial z} = \frac{1}{\rho} \cdot \frac{\partial}{\partial x}\left[ \left( \frac{\lambda}{c} + \frac{\mu_t}{\sigma_t} \right) \frac{\partial T}{\partial x} \right] + \frac{1}{\rho} \cdot \frac{\partial}{\partial y}\left[ \left( \frac{\lambda}{c} + \frac{\mu_t}{\sigma_t} \right) \frac{\partial T}{\partial y} \right] + \frac{1}{\rho} \cdot \frac{\partial}{\partial z}\left[ \left( \frac{\lambda}{c} + \frac{\mu_t}{\sigma_t} \right) \frac{\partial T}{\partial z} \right]
$$
Here, $T$ is temperature, $\lambda$ is thermal conductivity, $c$ is specific heat, $\mu_t$ is turbulent dynamic viscosity, and $\sigma_t$ is the turbulent Prandtl number.
Upon filling, the solidification process is dominated by heat conduction. The governing transient heat conduction equation is:
$$
\rho c_p \frac{\partial T}{\partial x} = \frac{\partial}{\partial x} \left( \lambda \frac{\partial T}{\partial x} \right) + \frac{\partial}{\partial y} \left( \lambda \frac{\partial T}{\partial y} \right) + \frac{\partial}{\partial z} \left( \lambda \frac{\partial T}{\partial z} \right) + \dot{E}
$$
where $c_p$ is the specific heat at constant pressure and $\dot{E}$ is the volumetric heat source rate.
Defect and Stress Prediction Models
Shrinkage porosity, a critical defect in steel casting, forms when liquid metal cannot adequately feed volumetric shrinkage in isolated liquid pools during solidification. The Niyama criterion is a widely accepted metric for predicting this defect, defined as:
$$
\frac{G}{\sqrt{\dot{T}}} \leq C_{Niyama}
$$
where $G$ is the local temperature gradient, $\dot{T}$ is the local cooling rate, and $C_{Niyama}$ is a material-dependent critical value. Locations where this criterion is violated are prone to microporosity. The temperature gradient $G$ is calculated from the nodal temperatures in the discretized model.
Thermal stresses arise due to differential cooling and constraints within the casting and mold. The evolution of stress was modeled considering the elastic-plastic behavior of the material at high temperatures. The equivalent (von Mises) stress, used to assess stress concentration and potential failure, is calculated as:
$$
\sigma_e = \sqrt{ \frac{1}{2} \left[ (\sigma_1 – \sigma_2)^2 + (\sigma_2 – \sigma_3)^2 + (\sigma_3 – \sigma_1)^2 \right] } = \sqrt{ \frac{1}{2} \left[ (\sigma_x – \sigma_y)^2 + (\sigma_y – \sigma_z)^2 + (\sigma_z – \sigma_x)^2 + 6(\tau_{xy}^2 + \tau_{yz}^2 + \tau_{zx}^2) \right] }
$$
where $\sigma_1, \sigma_2, \sigma_3$ are the principal stresses and $\sigma_e$ is the equivalent stress.
Simulation Setup: Model, Materials, and Initial Processes
The subject of this steel casting study is a shearer rocker arm housing with a complex geometry, featuring thin-walled sections (10-20 mm) and thick bosses (70-80 mm). Its overall envelope dimensions are 1742 mm x 1040 mm x 524 mm, with a net casting weight of approximately 1462 kg. The material specified is ZG20SiMn cast steel, chosen for its favorable combination of strength and toughness in heavy-load applications.
| C | Si | Mn | Mo | Cr | Ni | S | P |
|---|---|---|---|---|---|---|---|
| 0.18 | 0.71 | 1.13 | 0.12 | 0.09 | 0.05 | 0.021 | 0.017 |
| Temperature (°C) | Density (g/cm³) | Enthalpy (kJ/kg) | Thermal Conductivity (W/m·K) |
|---|---|---|---|
| 25 | 7.81 | 111.31 | 34.80 |
| 653 | 7.65 | 358.67 | 30.89 |
| 1053 | 7.48 | 650.01 | 38.51 |
| 1453 | 7.20 | 984.52 | 38.59 |
| 1853 | 6.83 | 1508.21 | 39.05 |
The initial process design involved two distinct gating systems for this vertical steel casting:
- Top-Gating System: Molten metal is introduced from the top of the mold cavity. This design is simple but can lead to high impact turbulence and oxide formation.
- Bottom-Gating System: Metal enters from the base of the mold through multiple ingates, promoting a more tranquil, upward fill that minimizes turbulence and aids in slag trapping.
In both designs, feeder heads (risers) were placed on the thicker upper sections to compensate for solidification shrinkage. The mold material was defined as resin-bonded sand. The key process parameters for the initial simulation are summarized below.
| Parameter | Value |
|---|---|
| Pouring Temperature | 1580 °C |
| Pouring Rate | 35 kg/s |
| Mold Preheat Temperature | 250 °C |
| Heat Transfer Coefficient (Cast/Mold) | 500 W/(m²·K) |
Comparative Analysis of Initial Casting Processes
Filling Behavior and Temperature Evolution
The simulation of the filling stage revealed markedly different behaviors between the two gating systems. The top-gating system exhibited a chaotic fill pattern with multiple metal streams colliding, leading to potential oxide entrapment and an unsteady temperature distribution. The bottom-gating system, conversely, showed a calm, predictable rise of the metal front, resulting in a more favorable thermal gradient from the bottom (first to fill) to the top (last to fill). This stable fill is a significant advantage in complex steel casting to ensure soundness.
Solidification Sequence and Defect Prediction
The solidification analysis was critical. In both schemes, the solidification front progressed from the thinner sections at the bottom towards the thicker top sections and feeder heads. However, at a solid fraction of 90%, both processes showed a dangerous isolated liquid region in the thick-walled section at the base of the output shaft housing. This “hot spot” solidified last, after the feeder heads had already frozen, severing the feeding path. The Niyama criterion clearly flagged this area as highly susceptible to macro- and micro-shrinkage.
The quantitative defect prediction confirmed this visual analysis:
| Gating Scheme | Shrinkage Volume (cm³) | Percentage of Casting Volume |
|---|---|---|
| Top-Gating | 143.41 | 0.0775% |
| Bottom-Gating | 133.57 | 0.0721% |
While the bottom-gating process showed a marginal improvement (~6.9% less shrinkage volume), the defect level was still unacceptable for a critical load-bearing component. The stress analysis after complete cooling to 50,000 seconds further highlighted critical areas. High residual stresses, in some cases exceeding the material’s yield strength at room temperature, were concentrated at geometrical discontinuities such as thin-to-thick section transitions and around bore openings (e.g., motor shaft hole, output shaft hole). These stresses are the root cause of subsequent distortion or cracking during machining or service.
Process Optimization Strategy and Results
Based on the analysis, the bottom-gating scheme was selected for optimization due to its superior filling characteristics. The optimization strategy was multi-faceted, targeting both defect reduction and stress mitigation in the final steel casting:
- Enhanced Feeding: The size of the main feeder heads was increased to extend their solidification time. Furthermore, a strategically placed side feeder was added near the problematic output shaft housing thick section to provide a direct, localized source of feed metal.
- Controlled Cooling with Chills: External chills (metal inserts with high thermal conductivity) were added to specific thick sections adjacent to the isolated liquid zone. These chills extract heat rapidly, accelerating solidification in those areas and eliminating the last-to-freeze hot spot, thereby promoting directional solidification towards the feeders.
- Reduced Thermal Gradients: The mold preheat temperature was increased from 250°C to 300°C. This simple change slows down the initial cooling rate, reduces the severity of thermal gradients during early solidification, and consequently lowers the developed thermal stresses.
Quantitative Outcomes of Optimization
The impact of the optimized steel casting process was profound, as evidenced by the simulation results.
| Metric | Initial Bottom-Gating | Optimized Process | Improvement |
|---|---|---|---|
| Shrinkage Volume | 133.57 cm³ | 9.01 cm³ | 93.3% reduction |
| % of Casting Volume | 0.0721% | 0.0049% | — |
The solidification pattern was transformed. The previously isolated liquid pool was completely eliminated, and a clear, progressive solidification front from the casting extremities towards the enhanced feeder heads was established. The most significant achievement was the drastic reduction in residual stress concentrations.
| Monitoring Point (Location) | Initial Stress (MPa) | Optimized Stress (MPa) | Stress Optimization |
|---|---|---|---|
| 1 (Main Bearing Boss) | 696 | 351 | 49.57% |
| 2 (Motor Hole Thin Wall) | 850 | 523 | 38.47% |
| 3 (Idler Gear Bore) | 625 | 449 | 28.16% |
| 4 (Output Shaft Hole) | 632 | 56.4 | 91.08% |
The stress at the critical output shaft hole region was reduced by over 91%, moving it from a high-risk zone to a benign one. This not only minimizes the risk of casting cracks but also results in a more dimensionally stable component that is less prone to distortion during subsequent machining operations.
Conclusion
This detailed simulation study successfully demonstrated the power of virtual prototyping in optimizing a complex steel casting process. By comparing top and bottom gating, the fundamental superiority of a bottom-filled system for this complex geometry was established, leading to calmer filling and a marginally better defect profile. However, the initial bottom-gating design was still insufficient.
The systematic optimization—combining enlarged and strategically placed feeders, the application of chills to control local solidification rates, and an increased mold preheat temperature—proved highly effective. The optimized process achieved a near-elimination of shrinkage porosity (93.3% volume reduction) and a dramatic decrease in harmful residual stresses, particularly at critical stress concentration points where reductions of 38% to over 91% were realized. This holistic approach ensures the production of a sound, reliable ZG20SiMn steel casting with enhanced structural integrity, directly contributing to the performance and longevity of the heavy mining equipment it serves. The methodology underscores that advanced simulation is an indispensable tool for modern, quality-driven steel casting production, enabling first-time-right manufacturing even for the most challenging components.