Sand casting is a widely used manufacturing process for producing complex metal components, particularly in the automotive industry. The integration of rapid prototyping technologies with traditional sand casting methods has led to the development of rapid sand casting, which enables the quick production of prototypes and small batches. This approach combines the advantages of additive manufacturing, such as flexibility and speed, with the robustness of sand casting. In this study, we focus on the numerical simulation of the filling and solidification processes in rapid sand casting for engine cylinder heads, utilizing advanced software tools to optimize the process and ensure high-quality outcomes.
The engine cylinder head is a critical component in internal combustion engines, characterized by intricate geometries, numerous curved surfaces, and internal isolated hot spots. These features pose significant challenges in achieving precise casting with minimal defects. Traditional trial-and-error methods for process design are time-consuming and costly. Therefore, numerical simulation has become an essential tool for predicting and analyzing the behavior of molten metal during casting, thereby reducing development cycles and improving product quality. This paper explores the application of numerical simulation in rapid sand casting, emphasizing the use of ProCAST software to visualize and optimize the filling and solidification stages.
Rapid sand casting involves several key steps: creating a three-dimensional model of the cast part, designing the mold and gating system using CAD tools, performing numerical simulations to validate the design, fabricating patterns and core boxes via stereolithography (SLA), assembling the SLA prototypes, producing sand molds, and finally pouring the metal. The CAD/CAE system developed for this purpose allows for agile and accurate mold design, coupled with simulation-based verification. This integrated approach ensures that potential issues, such as shrinkage porosity and misruns, are identified and addressed early in the design phase.
The material used in this study is ZL105 aluminum alloy, known for its excellent castability and mechanical properties. Its chemical composition is detailed in Table 1. The alloy’s properties, such as thermal conductivity and solidification characteristics, play a crucial role in the simulation. The sand mold material is phenolic urethane resin-bonded silica sand, which provides the necessary thermal and mechanical stability during the casting process.
| Element | Composition (wt%) |
|---|---|
| Si | 4.5–5.5 |
| Cu | 1.0–1.5 |
| Mg | 0.4–0.6 |
| Al | Balance |
The governing equations for the numerical simulation of sand casting processes include the conservation of mass, momentum, and energy. The fluid flow during filling is described by the Navier-Stokes equations, while heat transfer during solidification is modeled using the Fourier heat conduction equation. The energy equation accounts for phase change effects, as the metal transitions from liquid to solid. The general form of the heat conduction equation is:
$$ \frac{\partial (\rho c_p T)}{\partial t} = \nabla \cdot (k \nabla T) + S $$
where \( \rho \) is the density, \( c_p \) is the specific heat capacity, \( T \) is the temperature, \( t \) is time, \( k \) is the thermal conductivity, and \( S \) represents source terms, such as latent heat release during solidification. For the solidification process, the enthalpy method is often employed to handle the phase change. The latent heat \( L \) is incorporated as:
$$ H = c_p T + f L $$
where \( H \) is the enthalpy and \( f \) is the liquid fraction, which varies from 0 (solid) to 1 (liquid). The liquid fraction is a function of temperature, typically defined based on the alloy’s solidus and liquidus temperatures. For ZL105 alloy, the solidus temperature is 536°C and the liquidus temperature is 622°C.
In sand casting, the interface heat transfer coefficient between the cast metal and the sand mold is critical for accurate simulation. This coefficient varies with temperature and is determined using nonlinear inverse methods, such as the Beck technique. The boundary conditions for the simulation are summarized in Table 2. These parameters ensure that the model replicates real-world conditions, enabling reliable predictions of defect formation.
| Parameter | Value |
|---|---|
| Pouring Temperature | 690–720°C |
| Pouring Rate | 0.75–1.5 kg/s |
| Initial Mold Temperature | 25°C |
| Interface Heat Transfer Coefficient | 200–1000 W/m²·K |
The mold design for the engine cylinder head involves selecting an appropriate parting plane to facilitate core making and removal, minimizing the number of parting planes to enhance dimensional accuracy, and preferring flat parting surfaces over curved ones to reduce deformation errors. Based on the cylinder head structure, the cores are divided into several sections, including lower left and right oil passage cores, water passage cores, upper oil passage cores, intake port cores, and exhaust port cores. The gating system is designed as a one-sided bottom gating system, which allows molten metal to enter from the bottom, minimizing turbulence and oxidation. The cross-sectional areas of the gating channels are calculated to ensure smooth filling. Riser design involves automatic computation of modulus, heat dissipation volume, and surface area to achieve effective feeding. The process parameters, such as draft angles and shrinkage rates, are adjusted according to the characteristics of SLA prototyping.
For the numerical simulation, the CAD model of the mold is exported in IGES format and imported into ProCAST. The GeoMesh module is used for grid generation, with element sizes typically set to half or one-third of the minimum wall thickness. The mesh consists of approximately 115,990 nodes and 508,256 elements, ensuring sufficient resolution to capture detailed phenomena. The simulation models the entire process from filling to solidification, outputting results such as temperature distributions, fluid flow patterns, and defect predictions.
The filling process simulation reveals that the one-sided bottom gating system enables smooth and sequential filling. Figure 4 illustrates the temperature distribution at different stages of filling. At t = 2 s, filling begins; by t = 5 s, 15% of the mold is filled; at t = 8 s, 40% is filled; at t = 11 s, 70% is filled; and at t = 16 s, filling is complete. The initial slow filling prevents sand erosion and gas entrapment, while the final stage, though slower due to gravity and metal pressure, successfully fills the riser without requiring an increase in hydrostatic head height. The filling time distribution across the mold, as shown in Figure 5, indicates that after 15% filling, the filling times at同一水平面 are nearly identical, demonstrating stable and uniform filling that avoids defects like inclusions and turbulence.
The solidification process is analyzed through the evolution of solid fraction over time. Figure 6 displays the solid fraction at various stages: t = 20 s, 45 s, 130 s, and 240 s. The outer walls of the cylinder head solidify first, followed by the internal regions. The risers effectively feed the thick sections at the upper part of the cylinder head, while the ingates supply the lower thick sections. No significant isolated liquid regions form, indicating that the design promotes directional solidification. The absence of such regions reduces the risk of macro-shrinkage.
Shrinkage porosity and voids are common defects in sand casting, often occurring in areas where solidification isolates liquid pools. The simulation predicts the location and extent of these defects, as shown in Figure 7. The results indicate that only uniformly dispersed micro-porosity forms within the casting, with overall porosity levels below 10%. This low porosity level signifies good casting quality, as confirmed by experimental validation. The actual cast cylinder head, produced based on the optimized process, exhibits clear contours and no visible defects. Microstructural examination of critical sections reveals dense and fine-grained structures, aligning with the simulation predictions.
The effectiveness of numerical simulation in sand casting is further demonstrated through process optimization. By adjusting parameters such as pouring temperature and gating design, potential defects can be mitigated. For instance, increasing the pouring temperature within the specified range can enhance fluidity but may also increase shrinkage risks. Thus, a balance must be struck using simulation-guided insights. The integration of CAD/CAE systems in rapid sand casting not only accelerates design but also reduces material waste and costs associated with physical trials.
In conclusion, numerical simulation using ProCAST software provides a powerful means to analyze and optimize the filling and solidification processes in rapid sand casting for engine cylinder heads. The one-sided bottom gating system ensures smooth filling, while the riser and ingate design facilitates effective feeding during solidification. The simulation results show no isolated liquid regions and minimal porosity, leading to high-quality castings. This approach underscores the importance of sand casting in modern manufacturing, especially when combined with rapid prototyping technologies. Future work could explore the application of machine learning algorithms to further refine simulation accuracy and automate process optimization in sand casting environments.
The mathematical models used in this simulation can be extended to other sand casting applications. For example, the heat transfer equation can be modified to account for different mold materials or alloy systems. Additionally, the fluid flow equations can incorporate turbulence models for more complex gating designs. The general form of the continuity and momentum equations for incompressible flow is:
$$ \nabla \cdot \mathbf{v} = 0 $$
$$ \frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla) \mathbf{v} = -\frac{1}{\rho} \nabla p + \nu \nabla^2 \mathbf{v} + \mathbf{g} $$
where \( \mathbf{v} \) is the velocity vector, \( p \) is pressure, \( \nu \) is the kinematic viscosity, and \( \mathbf{g} \) is gravitational acceleration. In sand casting, these equations are solved coupled with the energy equation to capture the interplay between flow and thermal fields.
Overall, the success of this numerical simulation highlights the transformative potential of digital tools in advancing sand casting processes. By leveraging computational power, manufacturers can achieve higher efficiency, better quality control, and faster time-to-market for complex components like engine cylinder heads. As sand casting continues to evolve, the integration of simulation technologies will play an increasingly vital role in meeting the demands of industries such as automotive and aerospace.