Numerical Simulation of Filling and Solidification in Rapid Sand Casting: Enhancing Quality for Sand Casting Products

As a researcher in manufacturing engineering, I have extensively explored the integration of rapid prototyping with traditional sand casting to produce high-quality sand casting products. This approach, particularly for complex components like engine cylinder heads, leverages numerical simulation to optimize processes and reduce development time. In this article, I delve into the methodology, simulation details, and outcomes of using ProCAST software for analyzing the filling and solidification phases in rapid sand casting. The focus is on improving the integrity and performance of sand casting products, which are critical in industries such as automotive and aerospace. Throughout this discussion, I will emphasize the advancements in sand casting products through computational tools.

The rapid sand casting process combines stereolithography (SLA) rapid prototyping with conventional sand casting techniques. This synergy enables the fast production of intricate sand casting products, especially for low-volume or prototype applications. The core of this method lies in the CAD/CAE system, which facilitates agile design and validation of casting molds. By simulating metal flow and solidification, we can predict defects like shrinkage porosity and optimize gating systems, ultimately enhancing the reliability of sand casting products. The following sections outline the entire workflow, from design to validation, with a strong emphasis on numerical aspects.

To begin, the three-dimensional model of the engine cylinder head—a representative sand casting product—is created using CAD software like UG NX. The material chosen is ZL105 aluminum alloy, known for its excellent casting properties, which is common in sand casting products for automotive engines. Its chemical composition is summarized in Table 1. This alloy has a液相线 temperature of 622°C and a固相线 temperature of 536°C, key parameters for simulation accuracy.

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

The casting process design involves several critical steps. First, the parting line is selected to ensure easy core removal and dimensional accuracy. For the cylinder head, multiple sand cores are used, including oil passage cores, water jacket cores, and intake/exhaust port cores. The gating system employs a one-side bottom-gating approach, where molten metal enters from the bottom to minimize turbulence and oxidation. This design is crucial for defect-free sand casting products. Riser design follows modulus calculations to provide adequate feeding for thick sections. The CAD model of the mold assembly, as shown in the image above, incorporates all these elements, highlighting the complexity of sand casting products.

In numerical simulation, the mesh generation is a pivotal step. Using ProCAST’s Geomesh module, the CAD model is discretized into finite elements. The mesh size is typically set to one-half to one-third of the minimum wall thickness, ensuring precision. For this cylinder head, the grid comprises approximately 115,990 nodes and 508,256 elements. The governing equations for fluid flow and heat transfer during filling and solidification are solved numerically. The momentum equation for incompressible flow is given by:

$$ \rho \left( \frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot \nabla \mathbf{v} \right) = -\nabla p + \mu \nabla^2 \mathbf{v} + \mathbf{f}_b $$

where \( \rho \) is density, \( \mathbf{v} \) is velocity, \( p \) is pressure, \( \mu \) is dynamic viscosity, and \( \mathbf{f}_b \) represents body forces like gravity. The energy equation accounts for phase change:

$$ \rho c_p \frac{\partial T}{\partial t} + \rho c_p \mathbf{v} \cdot \nabla T = \nabla \cdot (k \nabla T) + Q $$

Here, \( c_p \) is specific heat, \( T \) is temperature, \( k \) is thermal conductivity, and \( Q \) is latent heat release during solidification. The solidification kinetics for aluminum alloys can be modeled using the Scheil equation for microsegregation:

$$ C_s = k C_0 (1 – f_s)^{k-1} $$

where \( C_s \) is solid composition, \( C_0 \) is initial composition, \( k \) is partition coefficient, and \( f_s \) is solid fraction. These equations underpin the simulation of sand casting products, allowing us to predict temperature distributions and defect formation.

Boundary conditions are essential for realistic simulation. The mold material is phenolic urethane resin-bonded silica sand, with interface heat transfer coefficients varying nonlinearly with temperature. Initial conditions include pouring temperature, mold temperature, and pouring speed. Table 2 summarizes these parameters, which are critical for reproducing actual casting conditions for sand casting products.

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

The filling process simulation reveals that the one-side bottom-gating system ensures smooth filling. Figure 4 in the original work (not shown here due to restrictions) illustrates the temperature field at different times. At t = 2 s, filling begins; by t = 16 s, it completes. The sequential filling minimizes turbulence, which is vital for high-integrity sand casting products. The filling time distribution across the mold shows uniform progression, with no stagnant zones. This can be quantified using the filling efficiency \( \eta_f \), defined as:

$$ \eta_f = \frac{V_{\text{filled}}}{V_{\text{total}}} \times 100\% $$

where \( V_{\text{filled}} \) is filled volume and \( V_{\text{total}} \) is total cavity volume. For this case, \( \eta_f \) approaches 100% without issues, demonstrating the robustness of the gating design for sand casting products.

During solidification, the temperature gradient and cooling rate dictate microstructure and defects. The simulation outputs solid fraction contours over time. Initially, the outer walls solidify first due to higher heat extraction, followed by internal regions. The risers effectively feed thick sections at the top, while ingates supplement bottom regions. This directional solidification prevents isolated liquid pools, a common issue in sand casting products. The solidification time \( t_s \) can be estimated using Chvorinov’s rule:

$$ t_s = B \left( \frac{V}{A} \right)^n $$

where \( B \) is a mold constant, \( V \) is volume, \( A \) is surface area, and \( n \) is an exponent typically around 2. For the cylinder head, the modulus \( V/A \) varies across sections, influencing local solidification times. Table 3 compares moduli for different regions, highlighting feeding requirements.

Table 3: Modulus Calculations for Key Sections of the Cylinder Head
Section Volume (cm³) Surface Area (cm²) Modulus (cm)
Upper Thick Wall 150.2 85.3 1.76
Lower Thick Wall 120.8 72.1 1.68
Thin Wall Region 50.5 60.2 0.84

Shrinkage porosity prediction is a key outcome of simulation. The Niyama criterion is often used to identify regions prone to shrinkage defects in sand casting products. It is expressed as:

$$ G / \sqrt{\dot{T}} \leq C $$

where \( G \) is temperature gradient, \( \dot{T} \) is cooling rate, and \( C \) is a material-dependent constant. For ZL105 aluminum, values below 1 °C¹/²·s¹/² may indicate shrinkage. In our simulation, porosity is uniformly dispersed with overall porosity less than 10%, as shown in Figure 7 (referenced but not displayed). This low porosity level ensures that sand casting products meet stringent quality standards. The porosity fraction \( \phi \) can be calculated from simulation data:

$$ \phi = \frac{V_{\text{porosity}}}{V_{\text{casting}}} \times 100\% $$

where \( V_{\text{porosity}} \) is the volume of pores and \( V_{\text{casting}} \) is the casting volume. For this cylinder head, \( \phi < 10\% \), confirming good quality.

The production validation phase involves actual casting using the optimized design. The resultant cylinder head exhibits clear contours and no major defects, as seen in Figure 8 (not shown). Microstructural analysis reveals fine grains and dense morphology, aligning with simulation predictions. This success underscores the value of numerical simulation in developing reliable sand casting products. Furthermore, the rapid sand casting process reduces lead time from weeks to days, making it ideal for prototyping and small batches of sand casting products.

To generalize the findings, I have derived a set of best practices for simulating sand casting products. First, always use fine meshing in critical areas to capture thermal gradients. Second, calibrate boundary conditions with experimental data to improve accuracy. Third, incorporate alloy-specific solidification models, such as those for aluminum-silicon systems. The latent heat release during solidification can be modeled as:

$$ Q = L \frac{\partial f_s}{\partial t} $$

where \( L \) is latent heat per unit volume. This affects the cooling curve and defect formation. Additionally, for sand casting products with complex geometries, multiple iterations of simulation may be needed to optimize riser placement and gating design.

In terms of material properties, Table 4 lists thermal properties for ZL105 aluminum and silica sand, essential for simulation inputs. These properties influence heat transfer and solidification rates in sand casting products.

Table 4: Thermal Properties of Materials Used
Material Density (kg/m³) Thermal Conductivity (W/m·K) Specific Heat (J/kg·K)
ZL105 Aluminum 2700 150 900
Silica Sand Mold 1600 0.8 1100

The economic impact of this approach is significant. By reducing trial-and-error in foundries, numerical simulation cuts costs and material waste. For sand casting products like engine blocks or pump housings, even a 5% reduction in scrap rate can save thousands of dollars annually. Moreover, the ability to rapidly prototype sand casting products accelerates product development cycles, giving manufacturers a competitive edge.

Looking ahead, advancements in simulation software may include AI-driven optimization and real-time monitoring. However, the core principles remain: accurate modeling of fluid dynamics and heat transfer is paramount for high-quality sand casting products. I recommend further research on multi-scale modeling to link macro-scale simulation with microstructural evolution, enhancing the predictability of mechanical properties in sand casting products.

In conclusion, the integration of rapid prototyping and numerical simulation via ProCAST software effectively optimizes the filling and solidification processes for sand casting products. The case study of an engine cylinder head demonstrates that a one-side bottom-gating system promotes sequential filling and solidification, minimizing defects like shrinkage porosity. With porosity levels below 10% and no isolated liquid regions, the sand casting products achieve excellent quality. This methodology not only improves the reliability of sand casting products but also shortens development time, making it a valuable tool for modern manufacturing. As the demand for complex and durable sand casting products grows, such computational approaches will become indispensable in foundry operations worldwide.

Scroll to Top