In the automotive industry, the clutch pressure plate is a critical component responsible for transmitting engine power to the transmission system. Its structural integrity and dimensional accuracy are paramount for vehicle safety and performance. Typically manufactured through casting processes, pressure plates, especially those made from nodular cast iron, are susceptible to residual stresses due to non-uniform temperature fields during solidification. These stresses arise from variations in the yield limit and elastic modulus of the metal material, as well as uneven thermal plasticity at high temperatures, ultimately affecting dimensional precision, service life, and forming quality. Therefore, investigating the forming quality and residual stress distribution in clutch pressure plate castings is essential. This study employs the finite element software ProCAST to simulate the filling and solidification processes of an automotive clutch pressure plate casting. By analyzing temperature and stress evolution, we aim to elucidate the residual stress distribution and magnitude, thereby providing insights for optimizing casting工艺 parameters and enhancing product quality and production efficiency.
The pressure plate under consideration is a small-scale casting with a disk-like structure, featuring six lugs on the outer circumference—three larger and three smaller ones arranged symmetrically. The overall dimensions include an outer diameter of 215.9 mm, an inner hole diameter of 116 mm, and a total length of 254.7 mm. The average wall thickness is 10 mm, with a maximum of 11.9 mm at the lugs and a minimum of 5.5 mm at the plate surface. The material is nodular cast iron, specifically grade GGV30, which is akin to ductile iron but with vermicular graphite morphology; however, for consistency in this analysis, we refer to it as nodular cast iron due to its similar mechanical properties and widespread application. The chemical composition of this nodular cast iron is summarized in Table 1.
| C | Si | Mn | Cr | Cu | Mg |
|---|---|---|---|---|---|
| 3.65 | 2.65 | 0.27 | 0.019 | 0.097 | 0.018 |
The three-dimensional model of the pressure plate, along with the gating system, was assembled in UG software and imported into ProCAST’s Visual-Mesh module for preprocessing. Mesh generation utilized a tetrahedral element approach, with a global element size of 3 mm to balance accuracy and computational efficiency. After repairing surface connectivity issues and optimizing the mesh, the final model comprised 157,390 surface elements, 3,225,544 volume elements, and 569,502 nodes, ensuring robust numerical simulation. The meshed model is depicted below, highlighting the intricate geometry of the pressure plate and gating system.
Simulation parameters were defined to replicate industrial casting conditions. The mold material was resin sand, and the boundary conditions included an interface heat transfer coefficient of 500 W/m²K. The pouring temperature was set at 1419°C, with a pouring speed of 0.55 m/s and a pouring time of 7.5 s under gravity pouring. The mold initial temperature was 25°C, and cooling occurred via natural air convection. The gravitational acceleration was 9.8 m/s² in the negative y-direction. The simulation ran for 50,000 steps until the temperature dropped below 500°C. These parameters are consolidated in Table 2.
| Parameter | Value |
|---|---|
| Material (Casting) | Nodular Cast Iron GGV30 |
| Material (Mold) | Resin Sand |
| Interface Heat Transfer Coefficient | 500 W/m²K |
| Pouring Temperature | 1419°C |
| Pouring Time | 7.5 s |
| Pouring Method | Gravity Pouring |
| Cooling Method | Natural Air Cooling |
| Mold Temperature | 25°C |
| Gravitational Direction | Negative Y |
| Gravitational Acceleration | 9.8 m/s² |
| Simulation Steps | 50,000 |
| Stopping Criterion | Temperature < 500°C |
The filling process was simulated first, revealing the temperature distribution over time. Metal fluid entered through the pouring cup, filled the sprue, runner, and ingates, and then flowed into the mold cavity via the lugs. The entire filling was completed in 7.35 s, with no defects such as cold shuts or misruns. However, significant temperature gradients were observed; at 5.09 s, the lower part near the inner circle and small lugs had cooled to approximately 1300°C, while at the end of filling, the maximum temperature difference across the casting reached about 124°C, as shown in Figure 3 (not labeled per instructions). This non-uniform cooling is a primary driver of residual stresses in nodular cast iron components.
To quantify the thermal behavior, the heat conduction equation governs the temperature evolution during solidification:
$$ \frac{\partial T}{\partial t} = \alpha \nabla^2 T $$
where \( T \) is temperature, \( t \) is time, and \( \alpha \) is thermal diffusivity, which for nodular cast iron can be expressed as \( \alpha = \frac{k}{\rho c_p} \), with \( k \) being thermal conductivity, \( \rho \) density, and \( c_p \) specific heat. The latent heat release during phase change further complicates this, modeled using the enthalpy method:
$$ H = \int_{T_{ref}}^T \rho c_p \, dT + f_s L $$
Here, \( H \) is enthalpy, \( f_s \) is solid fraction, and \( L \) is latent heat. For nodular cast iron, the solidification range and graphite formation influence \( f_s \), leading to variable shrinkage and stress development.
The solidification process was analyzed through solid fraction contours. The sequence initiated at the inner and outer circles due to chilling effects from the mold, progressed toward the central plate region, and concluded at the large lugs connected to risers. This pattern resembles a convergent solidification from both ends toward the middle, which can induce thermal strains. To detail this, characteristic points were selected on the pressure plate, as illustrated in Figure 5, and their solid fraction over time is plotted in Figure 6. Points 1 (inner circle) and 6 (small lug) solidified first, while points 3 (large lug) and 5 (central plate) solidified later, confirming the non-uniform solidification. The solid fraction \( f_s \) as a function of time \( t \) for these points can be approximated by a sigmoidal curve:
$$ f_s(t) = \frac{1}{1 + e^{-k(t – t_0)}} $$
where \( k \) is a rate constant and \( t_0 \) is the time at 50% solidification. For nodular cast iron, the eutectic reaction during solidification alters this curve, but the general trend holds.
Residual stress analysis followed the thermal simulation. The equivalent stress (von Mises stress) distribution on the pressure plate surface, shown in Figure 7, indicates that tensile stresses dominate, with higher values near the inner circle and large lugs—up to 360 MPa—posing a cracking risk. The stress generally decreases from the inner to outer regions. This stress pattern stems from differential contraction during cooling; areas that solidify earlier constrain later-solidifying regions, generating internal stresses. The constitutive relation for stress in nodular cast iron during cooling can be expressed using thermo-elasto-plastic theory:
$$ \sigma = E \epsilon + \sigma_{th} $$
where \( \sigma \) is stress, \( E \) is Young’s modulus (temperature-dependent for nodular cast iron), \( \epsilon \) is elastic strain, and \( \sigma_{th} \) is thermal stress contribution from temperature gradients. The equivalent stress \( \sigma_{eq} \) is calculated as:
$$ \sigma_{eq} = \sqrt{\frac{1}{2}\left[(\sigma_x – \sigma_y)^2 + (\sigma_y – \sigma_z)^2 + (\sigma_z – \sigma_x)^2 + 6(\tau_{xy}^2 + \tau_{yz}^2 + \tau_{zx}^2)\right]} $$
For nodular cast iron, the graphite nodules act as stress concentrators, but overall, the matrix behavior dictates bulk stress.
To delve deeper, stress evolution at specific nodes (points 1-6 in Figure 8) was tracked over time, as graphed in Figure 9. The curves exhibit three phases: rapid increase during initial solidification due to constrained shrinkage, a decrease during eutectic expansion, and a gradual increase post-solidification from differential cooling. This behavior aligns with the properties of nodular cast iron, where graphite expansion during eutectic reaction can temporarily relieve stresses. The stress rate can be modeled as:
$$ \frac{d\sigma}{dt} = E \frac{d\epsilon}{dt} – \beta \frac{dT}{dt} $$
with \( \beta \) as the thermal expansion coefficient. For nodular cast iron, \( \beta \) varies with temperature, contributing to residual stress buildup.
Deformation analysis focused on the z-direction (normal to the pressure plate working surface). The displacement contour, Figure 10, reveals an “outer convex, inner concave” pattern, reducing flatness accuracy. Points along a horizontal line (points 1-14 in Figure 11) were analyzed for z-displacement, with data in Table 3 and a plot in Figure 12 showing symmetric displacement variation: negative (outward) at outer edges and positive (inward) near the inner hole. This distortion arises from uneven cooling; the outer region cools faster, contracting more and causing convex deformation, while the inner region is constrained by the core, leading to concave deformation. The displacement \( u_z \) can be related to strain via:
$$ \epsilon_z = \frac{\partial u_z}{\partial z} $$
and for small deformations, the total strain comprises elastic and thermal components:
$$ \epsilon_{total} = \epsilon_{elastic} + \epsilon_{thermal} = \frac{\sigma}{E} + \beta \Delta T $$
In nodular cast iron, the thermal strain \( \epsilon_{thermal} \) is significant due to high \( \beta \) values, exacerbating distortions.
| Point Number | Displacement (10⁻³ cm) |
|---|---|
| 1 | -1.4 |
| 2 | -1.0 |
| 3 | -0.4 |
| 4 | 0.3 |
| 5 | 1.08 |
| 6 | 1.78 |
| 7 | 2.5 |
| 8 | 2.18 |
| 9 | 1.58 |
| 10 | 0.74 |
| 11 | 0.07 |
| 12 | -0.41 |
| 13 | -0.78 |
| 14 | -0.83 |
The simulation results underscore the critical role of process parameters in residual stress management for nodular cast iron castings. To mitigate stresses, several strategies can be explored: optimizing riser design to ensure directional solidification, modifying pouring temperature to reduce thermal gradients, or implementing post-casting heat treatments like stress relief annealing. For instance, the stress relief process for nodular cast iron often involves heating to 500-600°C, holding, and slow cooling, which can reduce residual stresses by up to 80% without affecting mechanical properties. The effectiveness of such treatments can be predicted using kinetic models:
$$ \sigma_{relief} = \sigma_0 e^{-kt} $$
where \( \sigma_0 \) is initial stress, \( k \) is a temperature-dependent rate constant, and \( t \) is time.
Furthermore, the material properties of nodular cast iron play a pivotal role. Nodular cast iron, also known as ductile iron, exhibits high strength and ductility due to its spherical graphite nodules, which blunt crack propagation. However, its thermal properties, such as conductivity and expansion, vary with composition and microstructure. For example, the thermal conductivity \( k \) of nodular cast iron can be estimated as:
$$ k = k_{matrix} (1 – V_g) + k_{graphite} V_g $$
with \( V_g \) being graphite volume fraction. Typical values range from 30 to 50 W/mK, influencing cooling rates and stress generation. Additionally, the elastic modulus \( E \) of nodular cast iron decreases with temperature, approximated by:
$$ E(T) = E_0 \left(1 – \frac{T}{T_m}\right) $$
where \( E_0 \) is modulus at room temperature and \( T_m \) is melting point. This temperature dependence must be accounted for in accurate stress simulation.
In practice, the casting of nodular cast iron components like pressure plates requires careful control of inoculation and magnesium treatment to ensure proper nodularization, as deviations can lead to flake graphite formations and increased stress concentrations. The quality of nodular cast iron is often assessed by nodule count and shape, which affect mechanical performance. For instance, a higher nodule count generally improves toughness and reduces stress risers. The relationship between nodularity and tensile strength \( \sigma_t \) can be empirical:
$$ \sigma_t = A + B \cdot N $$
where \( N \) is nodule count per unit area, and \( A, B \) are material constants.
To summarize the findings, Table 4 consolidates key simulation outcomes for the nodular cast iron pressure plate, highlighting areas of concern and suggesting optimizations.
| Aspect | Observation | Implication | Potential Improvement |
|---|---|---|---|
| Temperature Field | Max ΔT ~124°C at filling end | High thermal gradients | Lower pouring temperature or preheated mold |
| Solidification Sequence | Inner/outer circles first, large lugs last | Non-uniform shrinkage | Riser redesign for directional solidification |
| Residual Stress | Tensile stresses up to 360 MPa at inner circle/large lugs | Cracking risk, reduced fatigue life | Stress relief heat treatment, optimized cooling |
| Deformation (Z-direction) | Outer convex (-0.014 cm), inner concave (+0.025 cm) | Reduced flatness accuracy | Stiffening ribs, uniform wall thickness |
| Material Behavior | Nodular cast iron shows significant thermal expansion | Exacerbates stresses and distortions | Alloy adjustments (e.g., lower Si content) |
From a broader perspective, numerical simulation tools like ProCAST enable virtual prototyping, reducing trial-and-error in foundries. For nodular cast iron castings, integrating microstructural models with thermo-mechanical simulations could further enhance accuracy. For example, the Johnson-Mehl-Avrami-Kolmogorov (JMAK) equation can describe phase transformations in nodular cast iron:
$$ f = 1 – e^{-(kt)^n} $$
where \( f \) is transformed fraction, \( k \) is rate constant, and \( n \) is Avrami exponent. Coupling this with stress analysis would allow prediction of transformation-induced stresses, common in nodular cast iron due to austenite-to-ferrite changes.
In conclusion, this study demonstrates the efficacy of ProCAST in analyzing residual stresses and deformations in automotive clutch pressure plate castings made from nodular cast iron. The simulations reveal that surface tensile stresses concentrate near the inner circle and large lugs, with a cracking tendency, while the working surface exhibits non-uniform residual stress distribution and z-direction deformation compromising flatness. These insights underscore the need for process optimizations, such as adjusted pouring parameters, modified gating systems, or post-casting treatments, to enhance the quality and reliability of nodular cast iron components. Future work could explore real-time monitoring or machine learning integration for predictive control in casting nodular cast iron parts. Ultimately, mastering the behavior of nodular cast iron through simulation paves the way for safer, more efficient automotive systems.