In the realm of advanced manufacturing, high-precision computer numerical control (CNC) machine tools are pivotal, with their accuracy heavily reliant on the dimensional stability and low residual stress of key components like bed castings. As a researcher focused on casting technology, I have observed that machine tool casting structures, such as those for beds, are often complex—featuring long guideways and significant wall thickness variations. During solidification and cooling, inhomogeneous temperature fields induce thermal stress, phase transformation stress, and structural stress, leading to deformation that compromises precision retention. Despite this, recent studies on residual stress and dimensional stability in machine tool casting are scarce, making stress-induced deformation and cracking a pervasive technical challenge in the foundry industry. Therefore, in this work, we employed numerical simulation and experimental methods to investigate residual stress and deformation in machine tool casting components, aiming to provide insights for optimizing casting processes and enhancing performance.
The core of our approach involved using JSCCAST simulation software to model the casting process for gray iron bed castings from precision turning centers (e.g., HTC2050) and horizontal machining centers (e.g., HMC50e). We complemented this with experimental residual stress measurements via 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 in machine tool casting applications.

From a theoretical perspective, residual stresses in machine tool casting arise from thermal gradients and material transformations during casting. The fundamental equations governing these phenomena include the heat conduction equation for temperature distribution and the stress-strain relationships for elastic and plastic deformation. For instance, thermal stress can be approximated using the formula: $$\sigma_{thermal} = E \cdot \alpha \cdot \Delta T$$ where \(E\) is the Young’s modulus, \(\alpha\) is the coefficient of thermal expansion, and \(\Delta T\) is the temperature difference across the casting section. In machine tool casting, these stresses are superposed with phase transformation stresses, given by: $$\sigma_{phase} = \beta \cdot \Delta V$$ where \(\beta\) is a material constant and \(\Delta V\) is the volume change during phase transformation. The total residual stress \(\sigma_{residual}\) in a machine tool casting can be expressed as: $$\sigma_{residual} = \sigma_{thermal} + \sigma_{phase} + \sigma_{structural}$$ where \(\sigma_{structural}\) accounts for stresses due to geometric constraints and cooling rates.
To model these effects, we first created 3D models of the bed castings in SolidWorks and exported them as STL files for import into JSCAST. The simulation parameters included material properties, boundary conditions, and process variables specific to machine tool casting. Key material properties for gray iron and ductile iron used in the simulations are summarized in Table 1, which highlights differences affecting stress generation. These properties are critical for accurate prediction of behavior in machine tool casting applications.
| Material | Density (kg/m³) | Thermal Conductivity (W/m·K) | Young’s Modulus (GPa) | Coefficient of Thermal Expansion (10⁻⁶/K) | Specific Heat (J/kg·K) |
|---|---|---|---|---|---|
| Gray Iron | 7100 | 46 | 110 | 12 | 540 |
| Ductile Iron | 7100 | 36 | 169 | 11 | 460 |
In the simulation of filling processes, we analyzed fluid flow patterns to assess turbulence and defect formation. For the HTC2050 machine tool casting, a gating system without choke design resulted in stable filling, whereas the HMC50e casting used a horizontal choke, leading to increased turbulence and potential gas entrapment. The governing fluid dynamics equation for incompressible flow during filling is the Navier-Stokes equation: $$\rho \left( \frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot \nabla \mathbf{v} \right) = -\nabla p + \mu \nabla^2 \mathbf{v} + \rho \mathbf{g}$$ where \(\rho\) is density, \(\mathbf{v}\) is velocity, \(p\) is pressure, \(\mu\) is dynamic viscosity, and \(\mathbf{g}\) is gravity. This simulation helps optimize gating systems to minimize defects in machine tool casting.
The temperature field and solidification process simulations revealed significant insights. As shown in Table 2, the temperature distribution varied across casting sections, with slower solidification in lower box regions leading to higher thermal stresses. The heat transfer during solidification can be described by the Fourier equation: $$\frac{\partial T}{\partial t} = \alpha_t \nabla^2 T$$ where \(\alpha_t\) is thermal diffusivity. For the machine tool casting, areas like guideways exhibited higher temperatures and longer solidification times, increasing susceptibility to stress.
| Casting Section | Maximum Temperature (°C) | Minimum Temperature (°C) | Solidification Time (s) | Stress Concentration Factor |
|---|---|---|---|---|
| Upper Box (HTC2050) | 1149 | 20 | 120 | 1.2 |
| Lower Box Bed (HTC2050) | 960 | 208 | 300 | 1.8 |
| Cross Guideway (HMC50e) | 1186 | 20 | 280 | 2.0 |
| Vertical Guideway (HMC50e) | 991 | 214 | 200 | 1.5 |
Defect prediction, particularly for shrinkage porosity, was another critical aspect. Using JSCCAST, we identified shrinkage defects in guideways and junctions, which contribute to contraction stress. The propensity for shrinkage can be quantified by the Niyama criterion: $$N_y = \frac{G}{\sqrt{T}}$$ where \(G\) is temperature gradient and \(T\) is local solidification time. Lower Niyama values indicate higher shrinkage risk in machine tool casting. Our simulations showed values below critical thresholds in guideway regions, necessitating design modifications like chill placement to promote directional solidification.
Stress and deformation simulations provided detailed maps of residual stress distributions. For the HTC2050 machine tool casting, the cross guideways experienced significant bending deformation, with tensile stresses up to 180 MPa and compressive stresses around -180 MPa. The displacement field analysis for HMC50e castings, summarized in Table 3, indicated larger deformations in the length direction (Y-axis) compared to width (X-axis) and height (Z-axis), aligning with the anisotropic shrinkage behavior of machine tool casting. The deformation \(\delta\) can be related to stress via Hooke’s law for plane stress: $$\delta = \int \frac{\sigma}{E} \, dx$$ where integration is over the casting length.
| Direction | Maximum Displacement (mm) | Minimum Displacement (mm) | Average Displacement (mm) | Standard Deviation (mm) |
|---|---|---|---|---|
| X (Width) | 15.324 | -12.180 | 1.572 | 8.5 |
| Y (Length) | 12.337 | -20.605 | -4.134 | 10.2 |
| Z (Height) | 8.673 | -7.506 | 0.583 | 5.3 |
Experimentally, we measured residual stresses using the blind-hole method, which involves drilling a small hole and measuring strain relief. The residual stress components \(\sigma_{xx}\) and \(\sigma_{yy}\) are calculated from measured strains \(\epsilon_1\), \(\epsilon_2\), and \(\epsilon_3\) using the equations: $$\sigma_{xx} = \frac{E}{2(1-\nu^2)} \left( \epsilon_1 + \epsilon_3 \right)$$ $$\sigma_{yy} = \frac{E}{2(1-\nu^2)} \left( \epsilon_1 – \epsilon_3 \right)$$ where \(\nu\) is Poisson’s ratio. For our machine tool casting samples, measurements were taken at different shakeout temperatures and before/after rough machining, as summarized in Table 4. This data underscores the impact of process parameters on residual stress in machine tool casting.
| Sample Condition | Measurement Point | Principal Stress σ₁ (MPa) | Principal Stress σ₂ (MPa) | Stress Component σₓₓ (MPa) | Stress Component σᵧᵧ (MPa) |
|---|---|---|---|---|---|
| Gray Iron HMC50e (200°C Shakeout) | 1 | -60.9 | -75.9 | -62.6 | -74.2 |
| 2 | -1.40 | -95.7 | -51.7 | -45.3 | |
| 3 | -7.90 | -52.7 | -49.6 | -11.0 | |
| Gray Iron HTC2050 (500°C Shakeout) | 1 | -19.7 | -63.1 | -25.2 | -57.5 |
| 2 | -49.0 | -66.8 | -49.2 | -66.6 | |
| 3 | -56.1 | -65.2 | -62.2 | -59.1 | |
| Ductile Iron Ram Box (Before Rough Machining) | 1 | -66.1 | -151 | -26.4 | -109.2 |
| 2 | -24.3 | -93.1 | -78.40 | -172.9 | |
| Ductile Iron Ram Box (After Rough Machining) | 1 | -55.3 | -139.6 | -91.40 | -160.0 |
| 2 | -91.60 | -71.40 | -134.0 | -65.10 |
Analyzing these results, we found that residual stresses in machine tool casting are predominantly compressive, with values aligning closely with simulation predictions. The gray iron castings exhibited higher residual stresses (up to -172.9 MPa at 500°C shakeout) compared to ductile iron castings (up to -65.5 MPa before machining), highlighting the material-dependent stress behavior in machine tool casting. Lower shakeout temperatures (e.g., 200°C) reduced residual stresses, whereas rough machining significantly increased them due to induced mechanical stresses. This underscores the importance of post-casting treatments, such as aging, to relieve stresses in machine tool casting.
Further discussion involves the implications for machine tool casting design and process optimization. The stress concentration factors derived from simulations can guide geometric modifications to mitigate stress risers. For instance, the stress intensity factor \(K\) for a crack-like defect in a machine tool casting can be estimated as: $$K = \sigma \sqrt{\pi a}$$ where \(\sigma\) is applied stress and \(a\) is defect size. Reducing residual stress through controlled cooling or alloy adjustments can enhance the fatigue life of machine tool casting components. Additionally, the use of simulation software like JSCAST allows for iterative design improvements, reducing trial-and-error in foundry practices for machine tool casting.
In conclusion, this study demonstrates the efficacy of combining numerical simulation and experimental measurement to address residual stress challenges in machine tool casting. Key findings include the identification of critical stress zones in guideways, the benefits of lower shakeout temperatures, and the adverse effects of rough machining on stress levels. Future work should explore advanced materials and hybrid manufacturing techniques to further improve the dimensional stability of machine tool casting. By integrating these insights, manufacturers can enhance the accuracy and longevity of high-precision machine tools, ultimately advancing industrial capabilities.
