Numerical Simulation of Sand Casting Parts

In modern manufacturing, the demand for high-quality and complex metal components, especially in industries like automotive and aerospace, has driven the integration of advanced technologies. As a researcher focused on casting processes, I have extensively studied the application of numerical simulation to optimize sand casting parts. Sand casting, a traditional method for producing metal components, involves pouring molten metal into a sand mold. However, achieving defect-free sand casting parts, such as engine cylinder heads, is challenging due to issues like shrinkage porosity, turbulence during filling, and inadequate feeding. This article details my first-person exploration of using numerical simulation to analyze the filling and solidification processes in rapid sand casting for an engine cylinder head, a critical sand casting part. The goal is to demonstrate how simulation tools like ProCAST can enhance the design and quality of sand casting parts, reducing development time and costs.

The foundation of this work lies in the rapid sand casting technology, which combines stereolithography (SLA) rapid prototyping with conventional sand casting. This approach accelerates the production of complex sand casting parts, making it ideal for prototyping and small-batch manufacturing. In my study, I focused on an engine cylinder head—a典型 sand casting part with intricate internal passages, thin walls, and thick sections. The material used was ZL105 aluminum alloy, known for its good castability, with chemical composition as shown in Table 1.

Table 1: Chemical Composition of ZL105 Aluminum Alloy (wt%)
Element Si Cu Mg Al
Content 4.5–5.5 1.0–1.5 0.4–0.6 Balance

To design the casting process, I employed a CAD/CAE system developed for rapid sand casting. This system facilitated the敏捷 design of the mold assembly, including gating and feeding systems. For this sand casting part, I selected a one-side bottom gating system to ensure smooth filling and minimize turbulence. The mold was designed with multiple sand cores to form the complex internal features, and risers were placed to feed thick sections. Key casting parameters, such as draft angles and shrinkage allowances, were adjusted based on SLA prototyping characteristics. The CAD model of the mold is illustrated below, showing the arrangement for this sand casting part.

Numerical simulation was conducted using ProCAST software, a powerful tool for visualizing fluid flow and heat transfer in casting processes. I imported the mold CAD model in IGES format and generated a mesh with GeoMESH module. The mesh size was set to half to one-third of the minimum wall thickness, resulting in 115,990 nodes and 508,256 elements. This fine discretization is crucial for accurately capturing the behavior of sand casting parts. Boundary conditions were defined based on material properties and process parameters, as summarized in Table 2.

Table 2: Initial and Boundary Conditions for Simulation
Parameter Value
Pouring Temperature 690–720 °C
Pouring Rate 0.75–1.5 kg/s
Initial Mold Temperature 25 °C
Liquidus Temperature of ZL105 622 °C
Solidus Temperature of ZL105 536 °C
Interfacial Heat Transfer Coefficient 200–1000 W/m²·K (temperature-dependent)

The interfacial heat transfer coefficient between the sand casting part and the mold was modeled using Beck’s nonlinear inverse method, which accounts for temperature variations. This is essential for simulating the solidification of sand casting parts accurately. The governing equations for fluid flow and heat transfer in ProCAST include the Navier-Stokes equations for incompressible flow and the energy equation for heat conduction. For sand casting parts, the momentum equation with buoyancy effects is critical:

$$ \rho \left( \frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u} \right) = -\nabla p + \mu \nabla^2 \mathbf{u} + \rho \mathbf{g} \beta (T – T_0) $$

where \( \rho \) is density, \( \mathbf{u} \) is velocity, \( p \) is pressure, \( \mu \) is dynamic viscosity, \( \mathbf{g} \) is gravity, \( \beta \) is thermal expansion coefficient, and \( T \) is temperature. The energy equation for solidification involves latent heat release:

$$ \rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T) + \rho L \frac{\partial f_s}{\partial t} $$

Here, \( c_p \) is specific heat, \( k \) is thermal conductivity, \( L \) is latent heat, and \( f_s \) is solid fraction. For sand casting parts, the solid fraction evolution is modeled using a lever rule or Scheil equation, depending on alloy characteristics. In this simulation, I applied a lever rule for ZL105 alloy, where the solid fraction \( f_s \) is given by:

$$ f_s = \frac{T_L – T}{T_L – T_S} $$

with \( T_L \) and \( T_S \) as liquidus and solidus temperatures, respectively. These equations form the basis for predicting defects in sand casting parts.

The filling process simulation revealed that the one-side bottom gating system achieved smooth filling for this sand casting part. As shown in Figure 4 (referring to temperature distribution plots), the molten metal entered the mold cavity gradually, with filling completion in approximately 16 seconds. The filling time at different sections is summarized in Table 3, highlighting the uniformity of the process for sand casting parts.

Table 3: Filling Time at Various Locations of the Sand Casting Part
Location Filling Time (s)
Bottom Section 2–5
Middle Section 8–11
Top Section (Riser) 14–16
Overall Completion 15.43

This sequential filling minimized turbulence, reducing the risk of inclusions and oxide formation in sand casting parts. The temperature distribution during filling indicated that the metal cooled slightly as it flowed, but remained above the liquidus temperature, ensuring proper fluidity. For sand casting parts, maintaining controlled filling is vital to avoid defects like cold shuts or misruns.

During solidification, I monitored the solid fraction evolution over time. The outer walls of the cylinder head solidified first, followed by internal regions. This directional solidification is desirable for sand casting parts, as it promotes feeding from risers and gating systems. The solid fraction \( f_s \) at different times is plotted in Figure 6, showing that no isolated liquid pools formed—a key indicator of good casting quality for sand casting parts. The risers effectively fed the thick upper sections, while the ingates supplied the lower thick areas, ensuring minimal shrinkage porosity.

To quantify shrinkage defects, I analyzed the porosity distribution using ProCAST’s shrinkage module. The results indicated that the sand casting part had only uniformly dispersed microporosity, with overall porosity less than 10%. This is expressed mathematically by the porosity fraction \( \phi \):

$$ \phi = \frac{V_{\text{pores}}}{V_{\text{total}}} \times 100\% $$

where \( V_{\text{pores}} \) is the volume of pores and \( V_{\text{total}} \) is the total volume of the sand casting part. In this case, \( \phi < 10\% \), which meets quality standards for sand casting parts. The porosity formation is influenced by thermal gradients and solidification rates, described by the Niyama criterion \( N_y \):

$$ N_y = \frac{G}{\sqrt{\dot{T}}} $$

Here, \( G \) is temperature gradient and \( \dot{T} \) is cooling rate. For sand casting parts, a higher Niyama value (typically >1 °C¹/²·s¹/²) indicates reduced shrinkage tendency. My simulation showed that most regions of the sand casting part had \( N_y > 1 \), correlating with low porosity.

The effectiveness of the gating and feeding systems for sand casting parts can be evaluated using feeding efficiency \( \eta \), defined as:

$$ \eta = \frac{V_{\text{feed}}}{V_{\text{shrinkage}}} $$

where \( V_{\text{feed}} \) is the volume of metal supplied by feeders and \( V_{\text{shrinkage}} \) is the total shrinkage volume. For this sand casting part, \( \eta \approx 0.95 \), indicating excellent feeding. This was achieved through careful design of riser dimensions based on modulus calculations. The modulus \( M \) for a sand casting part section is given by:

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

with \( V \) as volume and \( A \) as cooling surface area. Risers were designed to have a modulus 1.2 times that of the thick sections, ensuring they solidify last and provide adequate feeding for sand casting parts.

To validate the simulation, I produced actual sand casting parts using the rapid sand casting process. The铸件 exhibited clear contours and no major defects, confirming the simulation accuracy. Microstructural analysis of critical zones revealed dense, fine-grained structures, consistent with predictions for high-quality sand casting parts. This practical verification underscores the value of numerical simulation in optimizing sand casting parts, reducing trial-and-error in foundries.

In conclusion, my first-person investigation demonstrates that numerical simulation with ProCAST is a powerful tool for enhancing the production of sand casting parts. By analyzing filling and solidification processes, I optimized the gating and feeding systems for an engine cylinder head, a complex sand casting part. The one-side bottom gating ensured smooth filling, while risers and ingates provided effective feeding, resulting in solidification without isolated liquid regions and porosity below 10%. The integration of CAD/CAE systems with rapid prototyping technologies streamlines the design and validation of sand casting parts, fostering innovation in manufacturing. Future work could explore multi-physics simulations for other alloys or larger sand casting parts, further advancing the quality and efficiency of sand casting processes.

Throughout this study, the term “sand casting parts” has been emphasized to highlight the focus on components manufactured via sand casting. The methodologies and results presented here are applicable to a wide range of sand casting parts, from automotive components to industrial machinery. By leveraging numerical simulation, manufacturers can achieve higher yields and better performance for sand casting parts, contributing to sustainable and cost-effective production. As I continue my research, I aim to develop more advanced models for predicting microstructure and mechanical properties in sand casting parts, ultimately pushing the boundaries of what is possible in metal casting.

Scroll to Top