Research and Numerical Simulation of Lost Foam Casting for a Ductile Iron Long Fork Component

In our investigation, we focused on optimizing the lost foam casting process for a complex ductile iron lift arm used in agricultural machinery. This component, with overall dimensions of 375 mm × 142 mm × 552 mm and a weight of approximately 30 kg, presents significant challenges due to its intricate geometry featuring multiple bosses, recesses, transition fillets, and internal stepped through-holes. The material specified was QT600-3, a pearlitic ductile iron, requiring high internal soundness and surface quality for demanding load-bearing applications. The inherent difficulties of lost foam casting, such as foam decomposition dynamics and heat transfer complexities, necessitated a detailed numerical analysis to guide process design and eliminate shrinkage defects.

We employed ProCAST simulation software to analyze the filling and solidification stages under different gating system designs. The core of our methodology was the establishment of a detailed thermophysical model representing the coupled interactions between the molten metal, the evaporative foam pattern, the coating, and the unbonded sand mold. The finite element mesh was constructed with a refined element size of 4 mm for the complex casting surfaces, while coarser elements (15 mm and 100 mm) were used for the gating system and sand mold respectively, to balance accuracy and computational efficiency.

The thermophysical parameters for the lost foam casting simulation were carefully defined. The interfacial heat transfer conditions are critical in lost foam casting. We defined the interface between the hollow sprue and the foam as “EQUIV” (equivalent), while all other metal-foam interfaces were set as “CONIC”. The interfacial heat transfer coefficient (h) was assigned a value of 150 W/(m²·K), a representative value accounting for the insulating effect of the degrading foam and the coating layer. The initial temperature of the sand mold and foam pattern was set to room temperature (25°C), and the molten iron’s initial temperature was the pouring temperature of 1480°C. Boundary conditions included symmetric faces defined with “Symmetry” conditions, an ambient air cooling condition (“Heat”), and a pressure boundary condition to simulate the vacuum applied in lost foam casting. The pressure at the sand box outer surface was set to 0.05 MPa (vacuum), with atmospheric pressure (0.1 MPa) at the pouring cup surface, resulting in a pressure differential that aids filling. For the solidification and microstructural simulation, we coupled a thermal model with a microscopic model for ductile iron, accounting for graphite expansion during eutectic solidification. Key parameters for this micro-model were set as POROS=1, GRAPHITE=0.8, and FADING=0.8 within the ProCAST solver settings to model the self-feeding characteristics of ductile iron.

The governing equations for heat transfer during the lost foam casting process are fundamental. The energy conservation equation is expressed as:

$$ \rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T) + Q_{latent} + Q_{foam} $$

Where \( \rho \) is density, \( c_p \) is specific heat, \( T \) is temperature, \( t \) is time, \( k \) is thermal conductivity, \( Q_{latent} \) is the latent heat release rate during phase change, and \( Q_{foam} \) represents the heat sink/source term associated with the endothermic decomposition of the foam pattern. The interfacial heat flux between the metal and the foam is modeled as:

$$ q = h (T_{metal} – T_{foam}) $$

The filling process was simulated using a gravity-driven filling model coupled with fluid flow and temperature fields. The chemical composition of the QT600-3 ductile iron used in our study is summarized in the table below.

Table 1: Chemical Composition of QT600-3 Ductile Iron (wt.%)
Element C Si Mn P S Cr Ni Cu Mg Fe
Content 3.7 2.2 0.45 0.055 0.02 0.25 0.15 0.45 0.055 Bal.

Our initial process trials involved two classic gating schemes for lost foam casting: a bottom-gating system and a top-gating system. The bottom-gating design was relatively simple, consisting of a hollow sprue (Ø45 mm) and a horizontal runner transitioning to ingates at the base of the casting. A small slag trap (50 mm × 50 mm × 70 mm) was placed at the top of the upper arm, but no dedicated feeding riser was used, relying on the self-feeding capability of ductile iron. The top-gating system featured a central downsprue connected to a large riser at the top of the cylindrical section, with ingates branching to both arms. Additional side risers were placed on the upper and lower sections of the cylinder and the front of the upper arm.

We simulated both initial lost foam casting designs. The bottom-gating system showed a calm, progressive filling sequence from the bottom upward, with a total fill time of 24.08 seconds. The simulation indicated good surface quality potential with minimal slag entrapment. However, the solidification simulation revealed a critical flaw: the runner and sprue solidified prematurely, isolating the casting from external feed metal early in the solidification process. The last regions to solidify were the thick sections near the lower arm socket and the inner walls of the cylinder. The coupled micro-model predicted severe shrinkage porosity in these areas due to inadequate feeding, which was confirmed by actual trial castings that exhibited significant internal shrinkage defects.

The top-gating lost foam casting scheme filled the mold from the top downward in a controlled manner, with a final fill point at the front of the lower arm. The total fill time was slightly shorter. The solidification analysis showed that the large central riser/sprue remained liquid longer, providing some external feeding. However, the smaller side risers solidified before the adjacent casting hot spots, particularly at the upper cylinder junction. This resulted in isolated liquid pools and predicted shrinkage porosity in the upper thermal centers, again validated by physical trials.

The analysis of the temperature gradients and solidification sequences can be quantified. The local solidification time \( t_f \) at a point is crucial for predicting shrinkage. The Niyama criterion, often adapted for cast iron, relates thermal parameters to porosity risk. A simplified form considers the temperature gradient \( G \) and cooling rate \( \dot{T} \):

$$ \text{Niyama-like parameter} = \frac{G}{\sqrt{\dot{T}}} $$

Lower values of this parameter indicate a higher propensity for shrinkage porosity. Our simulations calculated this field for both initial designs, confirming the high-risk zones. To formalize the comparison, key simulation results from the initial lost foam casting trials are tabulated below.

Table 2: Comparison of Initial Lost Foam Casting Schemes
Parameter Bottom-Gating Scheme Top-Gating Scheme
Total Fill Time (s) 24.08 ~22.0
Filling Character Calm, bottom-up Controlled, top-down
Last Fill Point Top of arms/cylinder Front of lower arm
Primary Shrinkage Location (Simulation) Lower arm socket, inner cylinder wall Upper cylinder junction, arm junctions
Main Feeding Source Self-feeding only (inadequate) Central riser (partial), self-feeding
Major Defect in Trial Severe shrinkage porosity in lower thick sections Severe shrinkage porosity in upper thick sections

Driven by these results, we developed a modified step-gating system for the lost foam casting process. This design aimed to create a more favorable temperature gradient and ensure prolonged liquid connectivity to the critical hot spots. The system employed a two-tier ingate arrangement connected to a substantial vertical runner acting as a feeder channel. The hollow sprue (Ø45 mm x 250 mm high) fed into a solid foam vertical runner (60 mm x 60 mm x 240 mm). From this runner, two sets of ingates (60 mm x 40 mm x 40 mm, tapered to 40 mm x 30 mm cross-section at the casting) introduced metal at two levels: one near the lower cylinder hot spot and one near the upper cylinder hot spot. Small slag traps (30 mm x 50 mm x 20 mm) were placed at the top of the cylinder and the front of the upper arm solely for slag collection, designed to solidify quickly and avoid draining metal from the casting.

The simulation of this improved lost foam casting process showed marked improvement. The fill time was reduced to 16.64 seconds due to the shorter flow paths from the two ingate levels, ensuring higher metal temperature during mold filling and better foam degradation. More importantly, the solidification simulation demonstrated a significant change in the thermal field. The large vertical runner, acting as a feeder, remained liquid for an extended period, maintaining a feeding path to both the upper and lower thermal centers of the casting almost until the end of solidification. The solidification sequence became more directional toward this feeder channel. The predicted shrinkage porosity was drastically reduced. Only minor isolated liquid regions, potentially leading to minimal micro-shrinkage, were indicated in the extreme ends of the cylinder thick sections, with the only significant shrink predicted in the feeder runner itself, which is inconsequential as it is part of the gating system.

The effectiveness of the step-gating system can be analyzed through the thermal modulus. The thermal modulus \( M \), proportional to volume-to-surface area ratio \( V/A \), dictates the local cooling rate:

$$ M \propto \frac{V}{A} $$

Regions with a higher modulus solidify later. Our step-gating design effectively made the feeder runner the region with the highest modulus, ensuring it solidified last and fed the casting’s hot spots. The temperature difference \( \Delta T_{feed} \) between the feeder and the casting hot spot at the critical solidification fraction \( f_s \) is key for feeding:

$$ \Delta T_{feed}(f_s) = T_{feeder}(f_s) – T_{hotspot}(f_s) $$

In the improved lost foam casting design, \( \Delta T_{feed} \) remained positive for a much larger fraction of solidification compared to the initial schemes, enabling effective mass feeding to compensate for shrinkage.

To quantify the improvement in feeding efficiency, we can define a Feeding Efficiency Parameter \( \eta \) for a critical hot spot as the ratio of the time it is connected to a liquid feeder to its total local solidification time:

$$ \eta = \frac{t_{connected}}{t_{local}} $$

Where \( t_{local} \) is the local solidification time for the hot spot (time between liquidus and solidus), and \( t_{connected} \) is the duration within \( t_{local} \) during which a continuous liquid path exists to a feeder with sufficient pressure head. For the lower arm socket hot spot, \( \eta \) was nearly 0 for the bottom-gating scheme (no external feeder), very low for the top-gating scheme (distant, indirect feeder), and approached 0.9 for the step-gating scheme due to the direct, nearby lower-tier ingate connection to the large feeder runner.

The final validation involved producing actual castings using the optimized lost foam casting step-gating design. The results confirmed the simulation predictions. The castings were sound, with no macroscopic shrinkage cavities visible upon sectioning. Only a small area of concentrated micro-porosity was found in the lower hot spot region, which was acceptable for the application and far less severe than the defects in the initial trials. The successful outcome underscored the reliability of the numerical simulation in optimizing the lost foam casting process when a comprehensive thermophysical model is used.

Table 3: Summary of Key Parameters and Results for the Optimized Step-Gating Lost Foam Casting Process
Category Parameter Value or Description Role/Impact
Process Design Gating Type Two-tier Step Gating Reduces fill time, improves thermal gradient.
Vertical Feeder Runner Size 60 x 60 x 240 mm High thermal modulus ensures last solidification.
Ingate Cross-section (at casting) 40 x 30 mm (tapered) Prevents reverse feeding, promotes rapid filling.
Slag Trap Size 30 x 50 x 20 mm Small size ensures early solidification, collects slag only.
Simulation Input Interfacial Heat Transfer Coefficient (h) 150 W/(m²·K) Models foam/metal/coating interfacial resistance.
Pouring Temperature 1480 °C Initial condition for molten metal.
Vacuum Pressure -0.05 MPa (gauge) Boundary condition for lost foam casting filling.
Simulation Output Total Fill Time 16.64 s Faster than initial designs, better thermal state.
Predicted Major Shrinkage Confined to feeder runner Indicates effective feeding of casting.
Critical Solidification Sequence Directional toward feeder runner Achieved by design.
Niyama Parameter in Critical Zones Increased significantly Indicates lower porosity risk.
Experimental Result Internal Shrinkage No major cavities, minor micro-porosity in one hot spot Confirms simulation accuracy and process success.
Casting Soundness Accepted for application Objective of lost foam casting optimization achieved.

In conclusion, our study demonstrates the powerful synergy between advanced numerical simulation and practical foundry engineering for lost foam casting. By iteratively designing and simulating different gating systems for a complex ductile iron component, we identified the shortcomings of conventional bottom and top-gating in lost foam casting when applied to parts with dispersed, heavy sections. The proposed step-gating system, evolved from simulation insights, successfully addressed the feeding problem by creating a controlled solidification pattern directed toward a large, persistent feeder. The close correlation between the simulated predictions and the actual casting results validates the fidelity of our thermophysical model for the lost foam casting process. This approach, combining ProCAST simulation with an understanding of ductile iron solidification and lost foam casting characteristics, provides a robust methodology for designing defect-free castings, reducing trial-and-error costs, and improving yield in lost foam casting production.

The broader implications for lost foam casting technology are significant. The ability to accurately model the complex interplay of fluid flow, foam decomposition heat sink, heat transfer through coating layers, and solidification with phase transformations allows for the optimization of even highly intricate castings. Future work could involve refining the foam decomposition model parameters further or exploring the sensitivity of the process to variations in coating thickness and permeability. Nevertheless, this case study firmly establishes that numerical simulation is an indispensable tool for advancing the reliability and application range of lost foam casting, particularly for demanding materials like high-grade ductile iron.

Scroll to Top