Optimization of Casting Process for Thick-Section Steel Castings Based on Numerical Simulation

In the manufacturing of heavy machinery, the production of defect-free, high-integrity steel castings with substantial cross-sectional thickness remains a formidable challenge for foundry engineers. These thick-section steel castings, often serving as critical load-bearing components, are subjected to stringent non-destructive testing standards. Internal defects such as shrinkage porosity and cavities can lead to catastrophic failures under operational stress. The traditional approach to process design, heavily reliant on empirical knowledge and iterative trial-and-error methods, proves to be inefficient, costly, and unreliable for such complex geometries. The advent of computational numerical simulation has fundamentally transformed this landscape. By virtually modeling the filling, solidification, and defect formation processes, simulation software enables a scientific, predictive approach to process optimization before any metal is poured. This article delves into a comprehensive case study involving a main bearing pedestal, a quintessential thick-section casting, illustrating the systematic application of simulation technology to diagnose inherent flaws in an initial process and to develop and validate a superior, defect-minimizing design.

The component under investigation is a main bearing pedestal for a heavy-duty drilling pump. Its function necessitates exceptional pressure tightness and structural integrity, mandating a very high degree of denseness with zero tolerance for shrinkage defects in its critical zones. With a maximum sectional dimension of 210mm x 200mm, this part classifies as a thick-walled steel casting. The material specified is ZG230-450 (a cast carbon steel), with a finished casting weight of approximately 666 kg. The initial, conventionally designed casting process is schematically represented in the provided figures. This scheme employed a vertical pouring orientation with a single large-diameter (φ300mm x 450mm) exothermic sleeved riser placed atop the casting. To promote directional solidification, external chills were positioned at the bottom arc and side regions. A bottom gating system was utilized to ensure a tranquil fill. Despite these seemingly sound traditional measures, ultrasonic testing (UT) per stringent Grade II requirements revealed unacceptable levels of shrinkage discontinuities concentrated in the core regions at both ends of the casting, farthest from the single riser.

The imperative for numerical simulation in modern foundry engineering cannot be overstated. For critical steel castings, the financial and temporal costs of physical trials are prohibitive. Simulation provides a virtual prototyping environment. We employed a standard simulation workflow: the 3D geometry of the casting, rigging system (risers, gates, chills), and mold was created using CAD software. This model was then imported into the ProCAST finite-element-based simulation suite. The simulation solves the fundamental governing equations of heat, mass, and momentum transfer during the casting process. The key outputs analyzed are the temporal evolution of the temperature field and the prediction of shrinkage defects using recognized criteria. For steel castings, the Niyama criterion is a widely accepted and reliable indicator for predicting the onset of microporosity. It is a function of the local thermal gradients (G) and cooling rates (R) during the final stages of solidification, expressed as:

$$
Niyama = \frac{G}{\sqrt{\dot{T}}} \approx \frac{G}{\sqrt{R}}
$$

where $\dot{T}$ is the cooling rate. Regions with a Niyama value below a critical threshold (specific to the alloy) are prone to shrinkage porosity. A lower G/√R ratio indicates a region that is cooling slowly with a shallow temperature gradient, creating conditions ideal for pore nucleation and growth as interdendritic feeding ceases.

The simulation of the original process vividly revealed its shortcomings. The temperature field analysis during solidification showed that while the chilled regions solidified rapidly, the thermal centers at both longitudinal ends of the casting, distant from the single top riser, remained as isolated liquid pools for an extended period. These hot spots solidified last, disconnected from the feeding source (the riser), leading to internal shrinkage. The Niyama criterion plot conclusively predicted severe shrinkage porosity precisely in these end-core regions, which matched perfectly with the actual UT findings and physical slice inspection of the defective castings. This correlation validated the accuracy of the simulation model and its parameters for this specific steel casting.

The root cause analysis was clear: the original process failed to establish a controlled directional solidification pattern towards an effective feeding source across the entire casting, particularly across its long span. The single riser was insufficient to feed the massive thermal mass at its far end. The end chills, while accelerating surface cooling, were inadequate to overcome the large thermal inertia at the core of the 200mm thick sections. The solidification sequence was not unidirectional but resulted in multiple isolated hot spots. The fundamental requirement for sound steel castings is summarized by the volume balance during phase change:

$$
V_{liquid} \cdot \rho_{liquid} = V_{solid} \cdot \rho_{solid} + V_{defect}
$$

To achieve $V_{defect} \to 0$, the volumetric shrinkage ($\beta = (\rho_s – \rho_l)/\rho_s$) must be compensated by continuous liquid metal feed from the riser(s) until the casting section is completely solid. This requires a positive pressure gradient from the riser to the solidifying region, maintained by a properly oriented temperature gradient (G). The original design failed to maintain this gradient in the end zones.

Comparison of Key Parameters in Original vs. Optimized Processes
Parameter Original Process Optimized Process 1 Optimized Process 2
Pouring Orientation Vertical Horizontal
Number of Riser(s) 1 (Top, φ300mm) 1 (Top, φ300mm) 2 (End, φ325mm)
Riser Type Exothermic Sleeve Exothermic Sleeve
Chill Placement Bottom & Sides Both Ends (100mm thick) Center Arc (120mm thick)
Feeding Distance Coverage Insufficient for ends Improved but incomplete Fully Covered
Simulated Shrinkage in Critical Zone Severe Porosity & Cavities Reduced Porosity No Porosity

Guided by the simulation diagnosis, we formulated and tested two sequential optimization strategies. The overarching goal was to enforce a definitive directional solidification sequence, ensuring every part of the casting solidified progressively toward a dedicated, potent feeding source.

Optimization Strategy 1: The first logical step was to enhance cooling at the problematic ends to try and pull the solidification front from the ends towards the central riser. The casting was re-oriented to a horizontal position for improved filling stability and thermal symmetry. Substantial external chills (100mm thick) were applied to both end faces. The single top riser was retained. Simulation of this modified design showed a marked improvement. The severe shrinkage cavity was eliminated, but the Niyama criterion still indicated a risk of dispersed microporosity in the core of the end sections. The chills were effective at the surface but could not sufficiently depress the solidification time of the deep thermal center. The effective feeding distance from the single central riser was still being exceeded for these thick sections, a classic limitation in steel castings. The feeding distance (L) for a plain section can be empirically related to its thickness (T) and the efficiency of the end effect (chill), often expressed as $L = k \cdot \sqrt{T}$, where k is a multiplier affected by chill efficiency. Here, the effective L was still less than the half-length of the casting.

Optimization Strategy 2: The lesson from Strategy 1 was pivotal: for very thick sections, relying solely on chills to extend feeding distance from a remote riser is unreliable. The solution was to bring the feeding source directly to the problem area. We adopted the principle of “placing the riser on the hottest spot.” Using the modulus method—a foundational technique for steel castings design—we calculated the required riser size to feed the heavy end sections. The modulus (M) of a section is its volume (V) to surface area (A) ratio, $M = V/A$, which governs its solidification time. For a riser to feed a casting section effectively, its modulus must be greater than that of the section, typically $M_{riser} > 1.2 \cdot M_{casting}$. Based on this, two large exothermic-sleeved risers (φ325mm x 325mm) with high moduli were designed and placed directly on the two end flanges, the locations of the original hot spots. A chill was placed on the center arc section to ensure it solidified before the end sections, thus establishing a clear solidification sequence: Center Chill Zone → End Sections → End Riser Feed Path.

The simulation results for this final configuration were exemplary. The solidification fronts progressed cleanly from the chilled center towards the end risers. The Niyama plot confirmed the complete elimination of shrinkage porosity in all critical, mechanically stressed zones of the casting. Any minor predicted porosity was relegated to non-critical, heavily fed regions adjacent to the riser necks. The optimization was deemed a simulation success.

The ultimate validation came from the foundry floor. The optimized process (Strategy 2) was used to produce multiple units of different pump models. Non-destructive ultrasonic inspection of all castings confirmed 100% compliance with the demanding Grade II quality standard. This transitioned a problematic casting with a high scrap risk into a reliable, repeatable production item. The table below summarizes the production outcome:

Production Results Using the Final Optimized Casting Process
Pump Model Quantity Produced UT Pass Rate (Grade II) Key Improvement
F-800 / F-1300 / F-2200 156 units (total) 100% (156/156) Elimination of shrinkage in end-core zones.

The science behind a sound casting process, especially for heavy steel castings, hinges on mastering thermal management. The following principles, reinforced by this case study, are universal:

  1. Directional Solidification: This is the cardinal rule. Solidification must progress from the extremities of the casting towards the riser(s). The temperature gradient (G) must be positive towards the feeder. This can be expressed as a requirement for the temperature field T(x,t):

$$
\frac{\partial T}{\partial x} \text{ (in direction towards riser) } > 0 \quad \text{during solidification.}
$$

  1. Adequate Feeding Volume and Feed Path: The riser must contain sufficient liquid metal to compensate for the total volumetric shrinkage of the casting section it feeds: $V_{riser, min} \ge \frac{\beta \cdot V_{casting}}{η}$, where η is the feeding efficiency of the riser (enhanced by exothermic sleeves). Crucially, the feed path—the mushy zone between the riser and the solidifying region—must remain open (i.e., permeable) until solidification is complete.
  2. Modulus Guidance: The modulus method provides an excellent first approximation for riser sizing in steel castings, ensuring the riser solidifies later than the casting.
  3. Strategic Use of Chills: Chills are powerful tools to locally increase the solidification rate and manipulate the temperature gradient. They are most effective in creating directional solidification when used in conjunction with risers, not as a substitute for them in feeding thick sections.

The numerical simulation’s role is to quantify and visualize these principles for a specific geometry. It calculates the evolving G and R values everywhere, allowing engineers to assess the directional solidification criterion and the Niyama criterion in detail. It transforms qualitative rules into quantitative, color-coded maps that pinpoint risk areas.

Common Defect Prediction Criteria in Casting Simulation
Criterion Mathematical Form Physical Interpretation Application in Steel Castings
Temperature Gradient (G) $G = |\nabla T|$ Rate of temperature change in space. High G promotes directional solidification. Used to assess feed path viability and misrun risks.
Cooling Rate (R or $\dot{T}$) $R = \left| \frac{dT}{dt} \right|$ Rate of temperature change in time. Affects microstructure fineness. Linked to grain size and mechanical properties.
Niyama Criterion (Ny) $Ny = G / \sqrt{R}$ Ratio of gradient to sqrt(cooling rate). Low values indicate shrinkage porosity risk. Primary tool for predicting microshrinkage and macroshrinkage locations.
Solid Fraction (f_s) $f_s = \frac{V_{solid}}{V_{total}}$ Fraction of material solidified. Tracks from 0 to 1. Identifies isolated liquid pools and last-to-freeze zones.

In conclusion, the journey from a defective to a robust casting process for the thick-section main bearing pedestal underscores a modern paradigm in foundry engineering. The integration of numerical simulation is no longer a luxury but a necessity for the economically and technically sound production of high-value steel castings. It moves the process development cycle from the physical shop floor to the digital desktop, enabling rapid, cost-effective iteration and optimization. By accurately predicting thermal behavior and defect formation, tools like ProCAST allow engineers to design rigging systems that enforce the laws of solidification science—directional solidification, adequate feeding, and thermal gradient control. This case study exemplifies that for challenging thick-section geometries, a scientifically designed process using multiple, strategically placed risers based on modulus calculations, aided by chills for sequence control, and thoroughly validated through simulation, is the definitive path to achieving consistently sound, high-quality steel castings that meet the most demanding performance criteria.

Scroll to Top