In modern engineering, particularly within the demanding realm of internal combustion engines, the pursuit of higher efficiency and lower emissions has placed unprecedented thermal and mechanical loads on critical components. Components such as exhaust manifolds and turbocharger housings operate under severe thermomechanical fatigue conditions, experiencing cyclic stresses during repeated engine start-ups and shut-downs. This environment makes the fatigue durability of these parts a paramount concern for design optimization. Consequently, a deep and thorough understanding of the stress-strain response of materials under both static and cyclic loading across their entire service temperature spectrum is fundamental for accurate fatigue life assessment and reliable component design.
Among the materials chosen for these challenging applications, high-silicon molybdenum ductile iron castings, such as the grade QTRSi4Mo1 (with approximately 4% Si and 1% Mo), are prevalent. These alloys are favored for their good castability, acceptable high-temperature strength, and growth resistance. They are routinely employed in parts where service temperatures can soar to 760°C and beyond. While there exists a body of international research on the mechanical performance of various ductile iron castings at room and elevated temperatures, detailed investigations covering the full range from ambient to peak operating temperatures—especially regarding cyclic deformation behavior—are less common. Such comprehensive data is crucial for developing robust constitutive and fatigue life prediction models. This article presents a systematic study on the effects of temperature on the tensile properties and cyclic stress-strain response of a cast QTRSi4Mo1 ductile iron casting, providing essential insights for the fatigue durability evaluation of high-temperature engine components.
Experimental Material and Methodology
The material investigated is a ferritic ductile iron casting designated QTRSi4Mo1. Its chemical composition, typical for such heat-resistant grades, is detailed in Table 1. Test specimens for both static tensile and strain-controlled low-cycle fatigue (LCF) tests were extracted from a cast block measuring 50 mm x 50 mm x 300 mm. The geometries of the tensile and fatigue specimens are illustrated in the referenced figures, featuring standardized gauge sections to ensure consistent stress states and failure locations.

Static uniaxial tensile tests were conducted on a mechanical high-temperature testing machine. The tests spanned a wide temperature range: Room Temperature (RT ~25°C), 200°C, 300°C, 400°C, 500°C, 600°C, 700°C, and 760°C. A constant strain rate of $$2.5 \times 10^{-4} \, \text{s}^{-1}$$ was applied for all temperatures to assess the isothermal tensile properties including yield strength, ultimate tensile strength, elongation, and reduction of area.
Strain-controlled LCF tests were performed on a servo-hydraulic testing system. Prior to testing, the gauge length of each specimen was meticulously polished and cleaned to eliminate surface irregularities and contaminants that could serve as premature crack initiation sites. Testing was conducted at five distinct temperatures: RT, 200°C, 400°C, 500°C, and 760°C. A fully-reversed strain ratio (R = -1) with a triangular waveform was employed at a constant strain rate of $$5 \times 10^{-3} \, \text{s}^{-1}$$. Various total strain amplitudes were applied at each temperature, as summarized in Table 2. The fatigue life, $$N_f$$, was defined as the number of cycles to complete fracture or a 5% drop in the maximum tensile stress from its stabilized value.
| C | Si | Mo | Mg | Mn | S | P | Fe |
|---|---|---|---|---|---|---|---|
| 3.18 | 4.43 | 1.23 | 0.03 | 0.31 | 0.008 | 0.027 | Bal. |
| Temperature (°C) | Strain Amplitude (%) | Strain Ratio (R) | Strain Rate (s-1) |
|---|---|---|---|
| 25 (RT) | 0.25 | -1 | 5 × 10-3 |
| 0.3 | |||
| 0.4 | |||
| 0.5 | |||
| 200 | 0.25 | ||
| 0.3 | |||
| 0.4 | |||
| 400 | 0.25 | ||
| 0.3 | |||
| 0.4 | |||
| 500 | 0.15 | ||
| 0.25 | |||
| 0.4 | |||
| 760 | 0.15 | ||
| 0.3 | |||
| 0.4 |
Temperature Effects on Monotonic Tensile Properties
The results from the static tensile tests at various temperatures are consolidated in Table 3. To clearly visualize the influence of temperature, key properties are normalized with respect to their room temperature values and plotted in Figure 1. The data reveals distinct trends for strength and ductility parameters.
| Temperature (°C) | 0.2% Yield Strength (MPa) | Ultimate Tensile Strength (MPa) | Elongation (%) | Reduction of Area (%) |
|---|---|---|---|---|
| 25 | 500 | 627 | 13.5 | 7.7 |
| 200 | 451 | 588 | 8.7 | 7.7 |
| 300 | 406 | 558 | 8.6 | 4.2 |
| 400 | 386 | 525 | 11.4 | 11.6 |
| 500 | 306 | 325 | 18.9 | 21.1 |
| 600 | 144 | 157 | 31.3 | 34.8 |
| 700 | 63 | 70 | 36.8 | 44.7 |
| 760 | 39 | 42 | 38.4 | 45.8 |
Both the yield strength ($$ \sigma_{0.2} $$) and ultimate tensile strength (UTS) exhibit a continuous decline with increasing temperature, as expected due to thermal activation of dislocation motion and other softening mechanisms. A more detailed analysis of the yield strength data reveals a bilinear relationship with temperature, as illustrated in Figure 2. This relationship can be accurately described by a piecewise linear function:
$$
\sigma_{0.2}(T) = \begin{cases}
481.80 – 0.2624T & \text{for } T \leq 419^\circ\text{C} \\
806.15 – 1.0365T & \text{for } T > 419^\circ\text{C}
\end{cases}
$$
The transition point at approximately 419°C signifies a change in the dominant strengthening or softening mechanism. In contrast, the ductility metrics—percentage elongation and reduction of area—display a non-monotonic trend. They initially decrease, reaching a minimum in the 300°C to 400°C range, before increasing significantly at temperatures above 500°C. This ductility trough is indicative of a well-documented phenomenon in some ductile iron castings known as “400°C embrittlement” or intermediate temperature brittleness. This embrittlement can be attributed to dynamic strain aging, where solute atoms (like carbon or silicon) interact with dislocations, pinning them and leading to localized plastic instability and reduced ductility. The coincidence of this ductility minimum with the kink in the yield strength-temperature plot is noteworthy. Above 500°C, the material exhibits excellent ductility, which is beneficial for accommodating thermal strains but also indicates a significant loss of strength.
Cyclic Stress-Strain Response and Hysteresis Behavior
The cyclic deformation behavior of ductile iron castings is complex and highly temperature-dependent, as revealed by the strain-controlled LCF tests. The evolution of the stress amplitude with the number of cycles at different temperatures and strain amplitudes is shown in Figure 3.
At room temperature, the material exhibits a three-stage behavior: primary cyclic hardening, followed by a period of cyclic softening, and concluding with secondary hardening prior to failure. The softening stage becomes less pronounced at lower strain amplitudes. A similar but distinct three-stage evolution is observed at 200°C, where an initial rapid hardening leads to a peak stress amplitude, followed by softening and a final secondary hardening stage where the stress amplitude remains below the initial peak.
The behavior changes markedly at 400°C and 500°C. At these intermediate temperatures, the ductile iron castings display persistent cyclic hardening throughout most of the fatigue life, with no apparent softening regime. The hardening is monotonic at 400°C, whereas at 500°C, the stress amplitude saturates after the initial rapid increase before a final drop due to crack propagation. At the highest service temperature of 760°C, the behavior is dominated by cyclic softening from the first cycle onwards. This is a characteristic feature of high-temperature LCF, where time-dependent deformation mechanisms like creep and microstructural evolution (e.g., coarsening of precipitates or spheroidization of graphite) lead to a progressive loss of load-bearing capacity.
Cyclic Stress-Strain Curves and Constitutive Modeling
The cyclic stress-strain curve (CSSC) represents the stable or half-life stress-strain response of a material under cyclic loading. It is effectively the locus of the tips of stable hysteresis loops from tests at different strain amplitudes. The CSSC is commonly described by the Ramberg-Osgood relationship, which partitions the total strain amplitude ($$\varepsilon_a$$) into elastic and inelastic components:
$$
\varepsilon_a = \varepsilon_a^e + \varepsilon_a^{in} = \frac{\sigma_a}{E} + \left(\frac{\sigma_a}{K’}\right)^{1/n’}
$$
Here, $$\sigma_a$$ is the stable stress amplitude (taken at half-life, $$N_f/2$$), $$E$$ is the elastic modulus, $$K’$$ is the cyclic strength coefficient, and $$n’$$ is the cyclic strain hardening exponent. The hysteresis loops at half-life for various temperatures and strain amplitudes are plotted in Figure 4, along with the CSSC fitted using the above equation. The fitted parameters for each temperature are listed in Table 4.
| Temperature (°C) | Elastic Modulus, E (GPa) | Cyclic Strength Coefficient, K’ (MPa) | Cyclic Hardening Exponent, n’ |
|---|---|---|---|
| 25 | 166 | 1231.18 | 0.1417 |
| 200 | 127 | 637.67 | 0.0716 |
| 400 | 114 | 648.19 | 0.0517 |
| 500 | 111 | 618.17 | 0.0778 |
| 760 | 51 | 79.26 | 0.0516 |
Several key observations can be made. First, the hysteresis loops are noticeably asymmetric in tension and compression for a given strain amplitude, with the compressive stress magnitude often exceeding the tensile one. This asymmetry is a common feature in cast materials due to microstructural heterogeneities and the presence of graphite nodules. Second, the Ramberg-Osgood model provides an excellent fit to the experimental CSSC data across all temperatures, confirming its utility for describing the cyclic response of these ductile iron castings. The parameters $$K’$$ and $$n’$$ show a clear temperature dependence, with $$K’$$ decreasing significantly at high temperatures, reflecting the material’s softening.
Analysis of Masing and Non-Masing Behavior
A material is said to exhibit Masing behavior if the ascending branches of stable hysteresis loops from tests at different strain amplitudes coincide when translated to a common origin (typically the compressive tip). This implies that a single master curve can describe the material’s cyclic plastic flow. Non-Masing behavior, where the loops do not superimpose, indicates that the cyclic deformation resistance depends on the strain history or amplitude.
To investigate this, the half-life hysteresis loops at each temperature were translated to a common origin. The analysis reveals a significant temperature dependence in the cyclic deformation characteristics of these ductile iron castings. At lower temperatures (RT, 200°C, and 400°C), the translated loops do not superimpose, indicating clear non-Masing behavior. This is consistent with the complex cyclic hardening-softening-hardening sequences observed, suggesting evolving internal stresses and dislocation structures that depend on the applied strain amplitude.
In contrast, at 500°C and 760°C, the ascending branches of the translated loops align remarkably well, as shown in Figure 5, demonstrating classical Masing behavior. This transition likely relates to changes in the dominant deformation mechanisms. At high temperatures, enhanced dislocation climb and cross-slip, along with other recovery processes, may lead to a more stable and amplitude-independent dislocation structure, resulting in Masing behavior. This distinction has important implications for fatigue life prediction models, as non-Masing materials often require more sophisticated damage parameters that account for the mean stress evolution during cycling.
Microstructural Considerations and Implications for Component Design
The mechanical behavior described above is rooted in the microstructure of ductile iron castings. The matrix of QTRSi4Mo1 is primarily ferritic, which provides good ductility and thermal conductivity. The addition of silicon strengthens the ferrite and improves oxidation resistance, while molybdenum enhances high-temperature strength by forming stable carbides and solid solution strengthening. The spherical graphite nodules act as stress concentrators but also provide sites for energy dissipation during deformation.
The 400°C embrittlement is a critical microstructural phenomenon. It is often associated with the interaction of interstitial atoms with dislocations in the ferritic matrix. This dynamic strain aging leads to the formation of solute atmospheres (Cottrell atmospheres) around dislocations, pinning them and causing a yield point phenomenon and reduced uniform elongation. This embrittlement zone is a potential weak link for components subjected to thermomechanical fatigue (TMF), as the low ductility can accelerate crack initiation under the combined thermal and mechanical strains experienced during engine thermal cycles.
The transition in cyclic behavior from complex non-Masing to simpler Masing type with increasing temperature also reflects microstructural stability. The initial complex hardening at lower temperatures is associated with dislocation generation, interaction, and entanglement around the graphite nodules and second-phase particles. At high temperatures, recovery processes become active, dynamically balancing dislocation generation and annihilation, leading to a steady-state condition (saturation or softening) that is less dependent on the specific strain history.
For engineers designing exhaust manifolds or similar components from ductile iron castings, these findings underscore several key points:
- Material Selection and Processing: The chemistry and heat treatment must be optimized to mitigate the intermediate temperature brittleness as much as possible, perhaps by controlling the levels of certain interstitial elements or by modifying the matrix structure.
- Finite Element Analysis (FEA): Accurate constitutive models are needed for stress simulation. The piecewise yield strength function (Equation 1) and the temperature-dependent Ramberg-Osgood parameters (Table 4) should be implemented in FEA software to capture the correct stress response under monotonic and cyclic thermal loads.
- Fatigue Life Prediction: Life prediction models, especially for TMF, must account for the non-Masing behavior and mean stress evolution at lower temperatures (below 500°C). At higher temperatures, where Masing behavior and cyclic softening dominate, time-dependent effects like creep and oxidation become critical and must be integrated into the damage accumulation model.
- Operational Limits: Understanding that the material loses most of its strength above 500°C, while gaining ductility, is crucial. Designs must ensure that operational stresses remain within safe limits, often relying on the material’s ability to relax stresses through plasticity at these temperatures rather than its strength.
Conclusion
This comprehensive investigation into the temperature-dependent mechanical behavior of QTRSi4Mo1 ductile iron castings from room temperature to 760°C has yielded critical insights for high-temperature component design. The monotonic tensile properties show a continuous decrease in strength and a ductility minimum around 300-400°C due to dynamic strain aging embrittlement, followed by excellent ductility at higher temperatures. The cyclic deformation response is highly temperature-sensitive, evolving from complex hardening-softening sequences with non-Masing characteristics at lower temperatures to simpler saturation/softening behavior with Masing characteristics above 500°C. The constitutive behavior under cyclic loading is well-described by the Ramberg-Osgood equation, with parameters that are strongly temperature-dependent.
These results form a vital database for the accurate simulation and fatigue life assessment of critical engine components like exhaust manifolds. They highlight the necessity of using temperature-specific material data and models that reflect the underlying microstructural mechanisms, such as intermediate temperature embrittlement and the transition from non-Masing to Masing cyclic response. Future work could focus on developing unified viscoplastic constitutive models that seamlessly capture this full spectrum of behavior and on conducting out-of-phase thermomechanical fatigue tests to validate life prediction methodologies under service-relevant conditions for ductile iron castings.
