As a researcher in the field of advanced manufacturing, I have focused on the critical issue of residual stress in machine tool castings, which are essential components for high-precision CNC machine tools. The demand for ultra-precision machining requires machine tool castings to exhibit minimal casting stress and exceptional dimensional stability. However, the complex geometry of bed castings, characterized by long guideways and significant wall thickness variations, often leads to non-uniform temperature fields during solidification and cooling. This inhomogeneity induces thermal stress, phase transformation stress, and structural stress, resulting in deformation that severely compromises dimensional accuracy and retention. Despite its importance, recent studies on residual stress and dimensional stability in machine tool castings are scarce, making deformation and even cracking due to residual stress a common technical challenge in the foundry industry. Therefore, in this work, we employed numerical simulation and experimental measurements to investigate the residual stress and deformation in machine tool castings, aiming to provide insights for optimizing casting processes and enhancing performance.
The core of our approach involves using JSCAST casting simulation software to model the entire process for gray iron bed castings from two precision machining centers: the HTC2050 precision turning center and the HMC50e precision horizontal machining center. We complemented this with experimental residual stress measurements using the blind-hole method on both gray iron bed castings and ductile iron ram box castings. This integrated methodology allows for a comprehensive analysis of stress generation and mitigation strategies in machine tool castings.
In the numerical simulation phase, we utilized SolidWorks to create detailed three-dimensional models of the machine tool casting structures. These models were exported in STL format and imported into JSCAST for analysis. The software enables prediction of flow patterns, defect formation tendencies, solidification sequences, temperature distributions, and stress development. By setting appropriate computational conditions and material property parameters, we performed filling, solidification, and stress analysis. The governing equations for fluid flow, heat transfer, and stress evolution are integral to this process. For instance, the heat transfer during solidification can be described by the transient heat conduction equation:
$$ \rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T) + Q_{latent} $$
where \( \rho \) is density, \( c_p \) is specific heat, \( T \) is temperature, \( t \) is time, \( k \) is thermal conductivity, and \( Q_{latent} \) represents the latent heat release during phase change. The stress development, crucial for understanding deformation in machine tool castings, is governed by the thermo-elastic-plastic constitutive model, which accounts for thermal strain, elastic strain, and plastic strain. The total strain \( \epsilon_{total} \) can be expressed as:
$$ \epsilon_{total} = \epsilon_{thermal} + \epsilon_{elastic} + \epsilon_{plastic} $$
with thermal strain given by \( \epsilon_{thermal} = \alpha (T – T_{ref}) \), where \( \alpha \) is the coefficient of thermal expansion and \( T_{ref} \) is a reference temperature. These equations underpin the simulation predictions for machine tool casting behavior.
The filling process simulation revealed significant insights into the design of gating systems for machine tool castings. For the HTC2050 bed casting, a gating system without a choke design resulted in a smooth and stable filling process. In contrast, the HMC50e bed casting employed a horizontal choke design, which led to higher fluid velocity in the runner, increased turbulence tendency, and poor slag removal. This turbulent flow can entrain gas, causing porosity and slag inclusions in the machine tool casting, which may lead to stress imbalances and subsequent deformation during machining. Optimizing the gating system is thus critical for reducing defects and ensuring uniform stress distribution in machine tool castings.
Temperature field and solidification process analysis are paramount for assessing stress origins in machine tool castings. The simulation results indicated that the guideway regions of the bed casting exhibited higher temperature fields compared to other areas, with locations farther from the sprue and with greater wall thickness retaining heat longer. This non-uniform temperature distribution, especially between the upper and lower parts of the casting, is a primary source of thermal stress. Specifically, the upper part of the machine tool casting solidified faster due to higher heat transfer coefficients at core contacts, while the lower bed body, particularly the guideways, solidified more slowly. The substantial variation in solidification times across different sections, driven by wall thickness differences, promotes the development of phase transformation stress and structural stress. The temperature gradient \( \nabla T \) is a key factor, as stress \( \sigma_{thermal} \) can be approximated by:
$$ \sigma_{thermal} \approx E \alpha \Delta T $$
where \( E \) is Young’s modulus and \( \Delta T \) is the temperature difference. This highlights how localized hot spots in a machine tool casting contribute to residual stress buildup.
Defect prediction through simulation identified shrinkage porosity in critical areas of the machine tool casting, such as the cross rails, vertical rails, and junctions. These regions, experiencing significant temperature variations, are prone to contraction stress. The presence of shrinkage defects can act as stress concentrators, exacerbating deformation risks. The simulation output categorized defects, allowing for targeted process improvements like using chills to promote directional solidification or optimizing feeder placement to enhance feeding. The shrinkage volume \( V_{shrinkage} \) can be related to the solidification contraction rate \( \beta \) and the volume of the casting section \( V_{section} \):
$$ V_{shrinkage} \propto \beta \cdot V_{section} \cdot (1 – f_s) $$
where \( f_s \) is the solid fraction. This relationship underscores the need for careful control of solidification parameters in machine tool casting production.
Stress and deformation simulation provided quantitative data on residual stress magnitude and distortion patterns in the machine tool casting. For the HTC2050 bed, the cross rails exhibited pronounced bending deformation, with tensile stress at the rail roots reaching approximately 160–180 MPa and compressive stress on the rail surfaces ranging from -70 to -180 MPa. The vertical rails, being shorter, experienced lower stresses (tensile: 14–26 MPa, compressive: -27 to -52 MPa) and less deformation. The stress distribution \( \sigma(x,y,z) \) across the machine tool casting geometry was visualized, revealing high-stress concentrations at geometry transitions. The deformation displacement field for the HMC50e bed casting was analyzed along three axes: X (width), Y (length), and Z (height). The results showed that deformation in the Y-direction (length) was significantly larger than in the X and Z directions, with greater contraction observed toward the ends. This anisotropic deformation behavior is critical for designing machine tool castings to maintain dimensional accuracy over time. The displacement \( \delta \) can be expressed as a function of the strain field \( \epsilon \):
$$ \delta = \int \epsilon \, dl $$
where the integration is along the path of interest.
To validate and complement the simulation findings, we conducted experimental residual stress measurements using the blind-hole method. This technique involves bonding a strain rosette at specific points on the machine tool casting surface, drilling a small hole (1.5 mm diameter, 2 mm depth) at the center, and measuring the released strain to calculate residual stress. The equipment included a ZDL-II drilling device, YC-III stress instrument, and specialized strain rosettes. Measurements were taken on gray iron bed castings shaken out at different temperatures (200°C for HMC50e and 500°C for HTC2050) and on ductile iron ram box castings before and after rough machining. The residual stress \( \sigma_{res} \) is derived from the measured strains \( \epsilon_1 \), \( \epsilon_2 \), and \( \epsilon_3 \) at orientations of 0°, 45°, and 90° using the following equations for principal stresses:
$$ \sigma_{1,2} = \frac{E}{2} \left( \frac{\epsilon_1 + \epsilon_3}{1 – \nu} \pm \frac{\sqrt{2}}{1 + \nu} \sqrt{(\epsilon_1 – \epsilon_2)^2 + (\epsilon_2 – \epsilon_3)^2} \right) $$
where \( E \) is Young’s modulus and \( \nu \) is Poisson’s ratio. The results are summarized in the table below, which provides a comprehensive comparison of residual stresses at various measurement points on the machine tool castings.
| Point No. | Gray Iron HMC50e (200°C Shakeout) | Gray Iron HTC2050 (500°C Shakeout) | Ductile Iron Ram Casting (Before Rough Machining) | Ductile Iron Ram Casting (After Rough Machining) | ||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| σ₁ (MPa) | σ₂ (MPa) | σₓₓ (MPa) | σᵧᵧ (MPa) | σ₁ (MPa) | σ₂ (MPa) | σₓₓ (MPa) | σᵧᵧ (MPa) | σ₁ (MPa) | σ₂ (MPa) | σₓₓ (MPa) | σᵧᵧ (MPa) | σ₁ (MPa) | σ₂ (MPa) | σₓₓ (MPa) | σᵧᵧ (MPa) | |
| 1 | -60.9 | -75.9 | -62.6 | -74.2 | -151.0 | -172.9 | -160.0 | -164.2 | 46.8 | 6.1 | 6.4 | 46.5 | -131.2 | -130.7 | -30.3 | -110.4 |
| 2 | -1.4 | -95.7 | -51.7 | -45.3 | -93.1 | -139.6 | -134.0 | -98.6 | 37.8 | 24.0 | 34.6 | 27.2 | -156.3 | -153.2 | -110.4 | -105.3 |
| 3 | -7.9 | -52.7 | -49.6 | -11.0 | -55.3 | -91.4 | -65.1 | -81.5 | 54.9 | 37.7 | 40.5 | 52.1 | -156.5 | -155.0 | -105.3 | – |
| 4 | -19.7 | -63.1 | -25.2 | -57.5 | -26.4 | -91.6 | -83.6 | -34.4 | 39.1 | -65.5 | 23.8 | -50.2 | – | – | – | – |
| 5 | -49.0 | -66.8 | -49.2 | -66.6 | -78.4 | -71.4 | -73.0 | -25.8 | -59.1 | -25.8 | -59.1 | -29.8 | – | – | – | – |
| 6 | -56.1 | -65.2 | -62.2 | -59.1 | -109.2 | -91.3 | -42.1 | -30.5 | -65.5 | -51.4 | -44.5 | -107.3 | – | – | – | – |
The data from the table elucidates several key trends regarding residual stress in machine tool castings. First, the measured residual stresses are predominantly compressive, aligning with the simulation predictions. The maximum residual stress in the gray iron bed casting shaken out at 500°C was -172.9 MPa, significantly higher than the -95.7 MPa observed at 200°C shakeout temperature. This confirms that a lower shakeout temperature is beneficial for reducing residual stress in machine tool castings, as it minimizes thermal gradients during cooling. Second, the ductile iron ram casting exhibited lower residual stress before rough machining (maximum of -65.5 MPa) compared to gray iron castings, highlighting the material’s influence on stress development. However, after rough machining, the residual stress in the ductile iron machine tool casting increased markedly, reaching -156.5 MPa. This underscores that mechanical processing introduces additional stresses, potentially exacerbating deformation in finished machine tool castings. The relationship between machining-induced stress \( \sigma_{mach} \) and initial residual stress \( \sigma_{res} \) can be modeled as:
$$ \sigma_{total} = \sigma_{res} + \sigma_{mach} + \sigma_{thermal\_mach} $$
where \( \sigma_{thermal\_mach} \) accounts for heat generated during cutting. This cumulative effect necessitates stress-relief treatments post-machining for critical machine tool castings.
Further analysis of the simulation and experimental data allows us to formulate recommendations for optimizing machine tool casting processes. The filling system design should avoid horizontal chokes to reduce turbulence and gas entrapment; instead, a choke-free or vertical choke design is preferable for stable filling. Temperature field management during solidification is crucial—using chills or cooling channels in sand molds can promote uniform cooling, thereby mitigating thermal stress in machine tool castings. The placement of feeders and risers should be optimized based on solidification simulation to eliminate shrinkage defects that act as stress raisers. Additionally, the anisotropic deformation behavior, particularly the larger contraction in the length direction, must be accounted for in the pattern design stage by incorporating appropriate allowances or compensation curves. The overall distortion \( D \) of a machine tool casting can be approximated as a function of its geometry and stress state:
$$ D = f(L, W, H, \sigma_{avg}, E, \rho) $$
where \( L, W, H \) are length, width, and height dimensions, and \( \sigma_{avg} \) is the average residual stress.
In conclusion, this integrated study of numerical simulation and experimental measurement provides a robust framework for understanding and controlling residual stress in machine tool castings. The JSCAST software effectively predicted filling patterns, temperature distributions, defect formations, and stress states, while the blind-hole method offered quantitative validation of residual stresses. Key findings indicate that gray iron machine tool castings exhibit higher residual stresses than ductile iron counterparts, lower shakeout temperatures reduce residual stress, and rough machining significantly increases stress levels. To enhance the dimensional stability and performance of machine tool castings, future work should focus on advanced heat treatment techniques like vibration stress relief or aging treatments, real-time monitoring of casting parameters, and multi-scale modeling that couples macro-scale stress with microstructural evolution. The continuous improvement of machine tool casting processes is essential for meeting the escalating demands of high-precision manufacturing industries worldwide.
Moreover, the implications of residual stress extend beyond initial deformation; they affect the long-term fatigue life and dynamic performance of machine tool castings under operational loads. Therefore, a holistic approach encompassing design, simulation, process control, and post-processing is vital. As we advance, the integration of artificial intelligence for predictive analysis and optimization of machine tool casting parameters could revolutionize the field, enabling the production of near-net-shape castings with minimal residual stress and superior dimensional integrity. The journey toward stress-free machine tool castings is challenging but imperative for the next generation of ultra-precision machine tools.