Next Article in Journal
Fresh, Mechanical, and Thermal Properties of Cement Composites Containing Recycled Foam Concrete as Partial Replacement of Cement and Fine Aggregate
Next Article in Special Issue
An Investigation of the Anisotropic Mechanical Properties of Additive-Manufactured 316L SS with SLM
Previous Article in Journal
Bearing Capacity Performance and Optimal Design of a Novel Aluminum Alloy Arch Gusset Joint
Previous Article in Special Issue
Achieving an Excellent Strength and Ductility Balance in Additive Manufactured Ti-6Al-4V Alloy through Multi-Step High-to-Low-Temperature Heat Treatment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterizing and Modeling Transformation-Induced Plasticity in 13Cr-4Ni Welds upon Cooling

by
Jean-Benoit Lévesque
1,2,*,
Carlo Baillargeon
1,
Daniel Paquet
1,
Jacques Lanteigne
1 and
Henri Champliaud
2
1
Institut de Recherche d’Hydro-Québec, 1800 Boul. Lionel-Boulet, Varennes, QC J3X 1S1, Canada
2
Department of Mechanical Engineering, École de Technologie Supérieure, 1100 Rue Notre-Dame O, Montréal, QC H3C 1K3, Canada
*
Author to whom correspondence should be addressed.
Materials 2023, 16(22), 7166; https://doi.org/10.3390/ma16227166
Submission received: 11 October 2023 / Revised: 6 November 2023 / Accepted: 7 November 2023 / Published: 15 November 2023

Abstract

:
Dilatometric experiments were conducted with the main purpose of measuring the transformation-induced coefficients of 13% chromium and 4% nickel, which are martensitic stainless steel base and filler materials used for hydraulic turbine manufacturing. To this end, a set of experiments was conducted in a quenching dilatometer equipped with loading capabilities. The measurement system was further improved by means of modified pushrods to allow for the use of specimens with geometries that are compliant with tensile test standards. This improvement allowed for the measurement of the materials’ phases and respective yield strengths. The dataset was further used to determine the relationship between the applied external stress and the martensitic start temperature ( M s ) upon cooling. The TRIP coefficient’s K values for both the S41500 steel and E410NiMo filler material were measured at 8.12 × 10 5 and 7.11 × 10 5 , respectively. Additionally, the solid phase transformation model parameters for both the austenitic and martensitic transformation of the filler material were measured. These parameters were then used to model austenitic-phase-transformation kinetics and martensite transformation, including transformation-induced plasticity effects. Good agreement was achieved between the calculation and the experiments.

1. Introduction

Hydroelectric turbines commonly fail due to fatigue damage, resulting in unexpected and often long-term shutdowns, productivity losses, and significant expenses for repairs. Francis-type turbine runners are manufactured from large cast parts, with the blades welded to the crown and ring belt. A common base and filler material used for their manufacture is martensitic stainless steel with a nominal content of 13% chromium and 4% nickel. Over the years, multiple incidents have revealed that fatigue cracks typically initiate and grow within or near the fillet weld in the heat-affected zone (HAZ), which is the region of the structure known to experience the highest operation stresses and high residual stress due to welding processes and potentially containing casting and welding defects. Once discovered, these cracks are usually repaired on-site using gouging and welding operations.
It is well known that residual stresses left by manufacturing and repairing operations, such as welding or heat treatment, have an effect on the initiation and propagation of fatigue cracks [1]. For example, Savaria et al. [2] demonstrated that the residual stress field left by induction gear tooth heat treatment had a significant effect on the fatigue life of MS6414 steel. Paquet et al. [3] showed that the evolution of residual stresses after loading could lead, in some cases, to unexpected fatigue life, stating that great care should be taken to consider the phenomenon in fatigue life calculations. Itoh et al. [4] conducted experiments to study the effect of welding residual stresses on fatigue crack propagation rate and found that the effective stress intensity range, Δ K e f f , must consider residual stress distribution and related crack opening effects. This has been confirmed by Deschênes et al. [5] on 13Cr-4Ni steel for the use of a specimen specially designed to include welding residual stress while avoiding microstructure alteration. Liljedahl et al. [6] performed in situ neutron diffraction on welding residual stress during fatigue tests to analyze residual stress redistribution and successfully predicted fatigue crack growth. Jones et al. [7], used plastically bent specimens to evaluate superposition concepts for fatigue crack growth rate estimation, concluding that the method is accurate but relies on precise knowledge of the residual stress field. This level of knowledge of residual stress distribution arising from turbine runner assembly operations cannot yet be known experimentally. One way to circumvent this lack of experimental capacity is to model the fabrication steps.
In the last few decades, finite element method frameworks for the simulation of residual stresses left by welding operations have proven effective [8,9]. Among others, the works of Deng et al. [10] and Rikken et al. [11] have shown the importance of considering solid-state phase transformation (SSPT) in the simulation of the welding processes of materials experiencing martensitic transformation upon cooling. Other studies [12,13,14,15] have also highlighted the need to consider transformation-induced plasticity (TRIP) during martensitic transformation in order to compute a reliable stress distribution. These plastic strains are induced as the SSPT progresses while stress, σ , is applied, even if the stress is below the yield strength of the weakest phase. This phenomenon is commonly explained by two mechanisms: (i) a mismatch of the volume expansion coefficient results of microplasticity in the weaker phase [16], and (ii) the application of stress, which promotes the formation of preferred martensite variants [17], resulting in anisotropic permanent strains.
For the purpose of evaluating TRIP strains, Leblond et al. [18,19] proposed an undisputed model. In this model, the total strain rate, ϵ ˙ , is decomposed as follows:
ϵ ˙ = ϵ ˙ e + ϵ ˙ t h m + ϵ ˙ p
where ϵ ˙ e is the elastic strain rate, ϵ ˙ t h m is the thermo-metallurgical strain rate, which includes thermal strains as well as the strains resulting from the evolution of phases in proportion, and ϵ ˙ p is the plastic strain rate. This last term is a summation of the contributions of the classical plasticity strain rate, which is time-independent in this case, ϵ ˙ c p , and relies on external stress and temperature rates, and the transformation-induced plasticity strain rate, ϵ ˙ t p , which depends on the phase transformation rate. This summation is decomposed as the sum of
ϵ ˙ p = ϵ ˙ c p + ϵ ˙ t p
The term ϵ ˙ t p can be obtained via a relation to the following form:
ϵ ˙ t p = 3 2 K S ϕ ( f ) f ˙
where K is the TRIP coefficient, S is the stress deviatoric tensor, ϕ ( f ) is a normalized function of the volume increase in the product phase formed, for instance, expressed as ϕ ( f ) = f ( 1 ln f ) , and f ˙ is the transformation rate. In their article series, Leblond et al. thoroughly derive each of the terms of Equations (1)–(3) based on material properties and constitutive relations. They even supply a theoretical background for the determination of the TRIP coefficient based on the difference between the expansion curves of the phases and the yield strength of the parent phase. However, the determination of the TRIP coefficient can be conducted experimentally [20,21], which provides an even better evaluation of the final TRIP strains.
This paper proposes the measurement of the TRIP coefficients of both the base and filler 13Cr-4Ni martensitic stainless steels used for turbine runner fabrication and assembly. Attention is also given to the effect of the application of external stress on the martensitic start temperature. Along with the help of the model proposed by Leblond et al., the values measured are then used to evaluate the strains that occur as martensitic transformation takes place and transformation-induced plasticity is triggered by the application of an external load. The results were then compared to the experimental results.

2. Materials and Methods

2.1. Materials

The two materials under study are (1) a base metal, 13Cr-4Ni soft martensitic hot-rolled stainless steel, denominated as UNS S41500, as per ASTM-A240/A240M-18 [22], and (2) the filler material, a homogeneous electrode, designated E410NiMo according to ASME BPVC.II.C [23]. The base material used in this study was in the shape of a 60 mm thick plate (ArcelorMittal, Saint-Denis, France). The filler metal (Hobart Welding Products, Appleton, WI, USA) was deposited in several layers on a substrate plate in order to form an 80 mm thick weldment and allow further specimen sampling. Both materials were fully austenitized at 1050 ° for 1 h for microstructure homogenization and then furnace-cooled to room temperature in order to obtain fully martensitic microstructures. The micrographs of both materials, obtained at a 200× magnification, are presented in Figure 1. UNS S41500 shows a fully martensitic microstructure with few solid inclusions located at the prior austenitic grains boundaries. The E410NiMo filler material also presents a fully martensitic microstructure, with smaller prior austenitic grains and smaller martensite blocs. Many micrometric inclusions of flux and the microporosities left by the welding operations can be observed.
An analysis of the chemical content was conducted on both materials to ensure compliance with the standards. The results, provided in weight percentage (wt%), are shown in Table 1. The chemical contents were measured by glow discharge atomic emission spectrometry, except for carbon and phosphorus, which were measured by using the mean of the combustion infra-red detection technique.
Figure 2 shows the results of the dilatometric experiments carried out in the preamble of this study. It shows that both materials follow nearly the same thermal expansion in both the fully austenitic and martensitic states. However, significant differences arise when focusing on the regions where solid phase transformations occurred. During heating at the same rate, austenitic transformation starts at the same temperature but behaves differently as austenite grows more quickly in the E410NiMo filler material. Upon cooling, the temperature at which martensitic transformation starts also differs. A martensitic transformation start temperature of E410NiMo, M s 410 , was measured at 307.9 ± 8.4 ° C, which is noticeably higher than that of S41500, M s 415 , at 261.7 ± 7.0 ° C. This trend was also observed for the end of the martensitic transformation temperatures, with M f 410 measured at 161.1 ± 19.6 ° C and M f 415 measured at 127.2 ± 8.1 ° C.

2.2. Experimental Setup

A set of experiments was designed for both materials. These experiments needed to allow for a measurement of material strain while changing temperature. The apparatus also needed to permit a controlled tensile or compressive stress to be applied to the specimen in order to trigger transformation-induced plasticity upon cooling the specimen. The device used for this purpose was a DIL805 A/D/T™ quenching dilatometer (TA Instruments, New Castle, DE, USA). The quenching dilatometer was set to the tension mode, which enables the use of a hydraulic actuator to apply both tensile and compressive stresses. The actuator was linked to the sample through a load cell, allowing for the measurement of the applied force, F. Heating of the samples was carried out by using an induction coil, well-known for its fast heating capability. Cooling was controlled by the admission of an argon flow around the sample. Strain measurements were accomplished by using a linear variable differential transformer (LVDT). In this study, an improvement to the measuring system is proposed in order to allow for the use of a specimen geometry that complies with the ASTM E8/E8M [24] standard for tensile tests. A modification of the extensions of the LVDT, namely push-rods, was made to introduce wedges pointing toward the specimen’s reduced section. A small force, radial to the specimen, was applied on the wedge during the installation and was maintained during the test. The movement associated with the specimen strain was then followed as the wedge’s tip stayed on the sample due to friction. Prior to mounting a specimen, the LVDT must be tared to the initial gauge length, with a calibration specimen consisting of a tensile sample with 2 radial grooves precisely distanced by 10.0 mm. The use of this configuration removes the stress concentration in the reduced section of the usual geometry and thereby prevents the formation of plastic strain in the early loading stages, which, in turn, enables the measurement of stress-strain curves with small-size specimens at different temperatures. Figure 3 shows a schematic of the system used for the experiments.

2.3. Specimen Design

The specimen is designed to meet the ASTM E8/E8M [24] geometric requirements for tensile tests on a cylindrical specimen (testing device loading attachment geometry and capacity). The maximum load that can be applied by the actuator is 10 kN, and the room temperature ultimate tensile stress of the steels under study falls just below 1200 MPa. Therefore, a reduced section with a 3.25 mm nominal diameter was selected, allowing the materials to reach their ultimate strength during a tensile test. The length of the reduced section is 21 mm, with a fillet radius of 10.0 mm. A finite element simulation was performed to validate that no buckling would occur at the compression stress level in this study. Specimen fabrication was accomplished by using conventional machining, with the specimen axis oriented along the rolling direction of the S41500 steel and along the welding direction for the E410NiMo deposit.

2.4. Thermal Cycle

Since one of the goals of this work was to determine the materials’ transformation-induced plasticity parameters during martensitic transformation, the temperature of the specimen needed to be modulated for SSPT to occur. The thermal history is shown in Figure 4 and described here. The specimen was heated to a sufficient temperature to allow the martensite-to-austenite transformation to be completed. From room temperature, R T , the specimen was heated at a continuous heating rate of 10 ° C/s until it reached an austenitization temperature, T γ , of 1000 ° C. Once this temperature was attained, it was maintained for a 30-s soaking time before starting the cooling process. This short soaking time ensured the complete austenitization of the specimen’s reduced section and conformed to the thermal histories encountered during the welding process. From the thermal finite element method simulation of welding performed in previous work [25], it is known that the cooling rate that occurs in the weld, following the flux-cored arc welding (FCAW) process, is around 10 ° C/s as martensitic transformation starts. Therefore, this cooling rate was maintained from T γ until the specimen temperature dropped to 350 ° C, a temperature slightly above martensitic transformation, which is the M s of the studied steels. Both materials are known to experience martensitic transformation independently of the cooling rate, allowing the temperature to be held just before the martensitic transformation starts. A 55-s dwell time was maintained while a load was gradually applied to the specimen. The cooling was then resumed at the same previous rate and continued until R T was reached, allowing martensitic transformation to be completed.

2.5. Loading for TRIP Experiments

The amount of TRIP strain encountered during transformation is related to the deviatoric stress tensor and the yield stress of the parent phase. Multiple stress levels must be performed to determine the materials’ transformation-induced plasticity parameters. These stress levels must be lower than the parent austenitic phase yield stress, σ γ y , to prevent the plasticity of the parent phase prior to phase transformation. The yield stress of the austenitic phase was measured using the previously described apparatus and specimens for both materials. The temperature at which to perform the tensile test has to be higher than the martensitic start temperature, M s , which was measured at 307.9 ° C for E410NiMo, as stated in Section 2.1. Hence, a temperature of 350 ° C was chosen to perform the tensile tests after a full austenitization at 1000 ° C. The tensile tests were performed at a constant loading rate of 10 N/s. According to these tests, yield strengths of 90 MPa and 110 MPa were measured for S41500 and E410NiMo, respectively.
The tested stress levels were then determined from these values. Stress levels of the same magnitude in tension and compression were tested. This strategy allows for the transformation-induced strain to be isolated. As suggested by Neubert et al. [20], by averaging the absolute value of the final total strains of the paired stress levels, the strains due to phase thermal coefficient mismatches and prestresses vanished. Table 2 summarizes the stress levels tested for both materials. Two specimens were tested for each stress level. The bottom of Figure 4 shows the loading history over time. The magnitude of the force was determined prior to the experiment in order to obtain the intended stress value. This force was gradually applied during the 350 ° C isotherm in a 5 s time lapse before sample cooling was resumed.

3. Results and Discussion

3.1. Effect of the Application of Stress on the Total Strain upon Martensitic Transformation

The results of the experiments described in Section 2 are shown in Figure 5 for S41500, and Figure 6 shows the results of E410NiMo. These figures show the evolution of strain as the specimens are cooled from the 350 ° C isotherm after the imposition of the load. The corresponding load-to-yield stress ratios are given in the inset. The strain signals are filtered by means of a moving average to remove noise. Since the obtained data are similar for a given stress level, the test results of both specimens are averaged for each load level for the sake of clarity. At first, it can be observed that the final total strains are distributed roughly linearly with applied stress. The effect of applied stress on M s can also be seen in these figures. With increasing stress magnitude for both tension and compression, the martensitic transformation starts at higher temperatures.

3.2. Experimental Determination of the TRIP Parameter

From these experimental results, and with the assumption that the TRIP strain rate is a function of the transformation rate f ˙ , the TRIP parameter K can be determined, as indicated in (3). Figure 7 shows the final transformation-induced strains as a function of the magnitude of the applied stress for both materials. The absolute final strain, which is the strain left at the end of the experiment when reaching room temperature, is used for this purpose. During the tests, other sources, such as prestress from specimen preparation or internal thermal forces, contribute to plastic strain. However, it is assumed that these remain constant between the experiments. Therefore, the TRIP strain for a given stress magnitude can be obtained by calculating the difference between the final strain of the tests with and without applied stress. Then, the results of the experiments with the same equivalent stress, both in tension and compression, are averaged. This procedure reveals the amount of strain attributable to the transformation-induced plasticity as the other plastic strain sources vanish. The stress-strain couples are used to perform a linear regression, assuming that TRIP strain is zero when no stress is applied. The slope of this regression is the transformation-induced plasticity parameter.
The values obtained are K 415 = 8.12 × 10 5 and K 410 = 7.11 × 10 5 for both materials. The TRIP parameter for S41500 steel is slightly higher than that for E410NiMo, which complies with their respective measured yield stress. As stated in Section 2, the S41500 austenitic phase yield stress is lower than that of E410NiMo, which means that for an equivalent external stress ratio, the austenitic parent phase of S41500 stainless steel is more likely to develop microplastic strains as the harder martensitic phase grows, resulting in a volumetric expansion. The TRIP parameter value obtained for S41500 steel in the current study is in accordance with the 7.7 × 10 5 value obtained in the preamble of this work, where the tests were performed on a tensile test frame equipped with an electric furnace and a high-temperature extensometer. The value obtained for E410NiMo also compares well with the work of Neubert et al. [20], who obtained 5.3 × 10 5 using another type of equipment and a slightly different material.

3.3. Relationship between M s and Applied Stress

In addition to the TRIP parameters, the experimental results, along with the results of some higher stress ratios not suitable for TRIP parameters measurements, were used further to determine a relationship between the applied external stress and the martensitic transformation start temperature. The temperature and strain measurements are processed in the same manner as in previous work [25]. Polynomial regressions were performed on the strain recording: one on the part before and another on the part after the martensitic transformation had started. The point where the extrapolated regression curves cross each other is then referenced as the theoretical martensite start temperature, T K M . This process is repeated on each set of measurements to obtain scattered data on the effect of applied stress on the martensitic transformation of both materials at distinct plots. The regressions are constrained to reduce to the measured T K M with respect to the material when no stress was applied, which was measured as being 254.6 ± 4.2 ° C and 288.3 ± 3.6 ° C for S41500 and E410NiMo, respectively. The results are shown in Figure 8, along with a linear regression of those data points. The coefficients of determination, R 2 , are 0.875 and 0.836 for S41500 and E410NiMo, respectively, with S41500 steel showing better agreement. This might be because E410NiMo is a weld metal and, therefore, shows a columnar microstructure more prone to heterogeneity. Equations (4) and (5) give the above-mentioned relations for each material.
T K M 415 = 254.6 + 0.147 | σ |
T K M 410 = 288.3 + 0.197 | σ |

3.4. Validation of the Transformation-Induced Plasticity Model

The transformation-induced plasticity parameters calculated in Section 3.2 were used, along with the martensite start temperature equations defined in Section 3.3, to model the materials’ behavior under the experimental conditions described in Section 2. For this purpose, the evolution of the phase proportion must be modeled as temperature and time vary. The S41500 stainless steel material parameters for modeling phase transformation during heating and cooling were experimentally determined [25].
The austenitic transformation was modeled using the so-called Johnson–Mehl–Avrami–Kolmogorov (JMAK) model [26,27,28,29] by considering the isochronal heating of dilatometric specimens. The same methodology was used here to characterize the parameters for the E410NiMo filler material, which, in turn, allowed us to determine the activation energy ( E a ) by using the Kissinger method [30], the chemical reaction pre-exponential factor of the Arrhenius function, k 0 , and the JMAK exponent, n. These parameters allow the next equations to be used to model the proportion of the austenite phase during austenitic transformation:
f γ = 1 e ( k 0 s ) n
with
s = i = 1 m e E a R T i Δ t i
where Δ t i and T i are the i t h intervals of time and temperature, respectively, m is the total number of increments, and R is the universal gas constant.
The martensitic transformation was modeled by using a Koistinen–Marburger (K-M)-type model [31,32]. Therefore, the same set of experiments enabled the determination of the E410NiMo theoretical martensite start temperature without applied stress, T K M , and the transformation rate, a m . These parameters are used in Equation (8) to determine the proportion of the martensite phase upon cooling.
f α = 1 e a m T K M T
The parameters for both materials’ phase transformation models are listed in Table 3. However, in this study, the values of T K M refer to the previously defined Equations (4) and (5). The martensitic and austenitic phase coefficients of thermal expansion (CTE) of S41500 and E410NiMo are also listed in Table 3. These material parameters enable the modeling of the thermal and transformation strains associated with the proportion of the austenitic and martensitic phases and the transformation from one phase to another.
Figure 9 compares the modeled total strains to the strains measured during the transformation-induced plasticity experiments for S41500 stainless steel. Each of the strain rate components of Leblond’s model was computed using the previously defined parameters and by using an explicit scheme. A good agreement between the modeled behavior and experimental data was achieved. For stress ratios between 0.36 and 0.36 , the maximum difference in the strain is 0.025 when room temperature is reached. For the higher stress ratio of 0.73 , a maximum difference of 0.08 % is recorded for the tensile test as the cooling ends. An additional experiment was conducted with a stress ratio of 1.09 , exceeding the austenite yield strength. This last test makes it possible to visualize the behavior of the material when it undergoes a classical plastic strain component as the transformation occurs. For this test, the model overestimates the strains in the early transformation stage. A discrepancy is formed as the strains accumulate at a faster rate in the model. However, as the transformation continues, the gap gradually narrows and ends with a difference of 0.075 % between the two strains.
Figure 10 shows the results obtained using the model, along with the experimental values for the E410NiMo filler metal. Again, the model is in good agreement with the experimental data, especially for the lower stress ratios between 0.60 and 0.60 , where the largest strain gap is 0.05 % for the higher applied tensile stress ratio ( 0.60 ) test at room temperature. Two more tests with high tensile stress ratio levels were performed at 0.90 and 1.20 for the cases where the applied stress was close to or above the yield strength. A discrepancy is noticeable at the end of the 0.90 stress ratio test. For the 1.20 stress ratio, the model overestimates the strains at the start of martensitic transformation. However, as the cooling proceeds, the gap between the measured and modeled strains vanishes.
Feeding the model with the parameters presented here enables a consistent simulation of the material behavior during martensitic transformation under the effect of stresses for both materials. However, a discrepancy is observed for stress ratios close to or above the material yield strength. This disparity is due to different sources. One of these is the difficulty associated with measuring material properties for phases separately. In the present study, the austenite and martensite stress-strain curves were measured at 350 ° C to obtain a realistic elastic limit. However, these properties are used throughout the entire cooling process as this is the only temperature at which the material properties were evaluated. Additionally, although this temperature has been chosen to be close to M s , it may not be high enough to completely avoid triggering the martensitic transformation at the beginning of plasticity. Moreover, the plasticity model considers a perfectly plastic material for both phases. This choice has been made since the TRIP coefficients are determined experimentally and, therefore, account mainly for the parent phase yield. Finally, in this study, the transformation parameter α m of the K-M model was kept constant throughout the modeling of martensitic transformation. Recently, Liu et al. [33] reported better agreement with experiments when α m is defined based on phase proportion and applied stress, especially for the early transformation steps.

4. Conclusions

In this paper, the material parameters required to model the transformation-induced plasticity of the 13Cr-4Ni base metal, as well as those of the filler material, during martensitic transformation were determined experimentally. For this purpose, the experiments were performed in a quenching dilatometer. The following conclusions can be drawn:
  • An improvement to the measurement system was made, allowing for the use of specimen geometry that met the requirements of the tensile test standards.
  • The TRIP coefficient’s K values for both S41500 steel and the E410NiMo filler material were determined as 8.12 × 10 5 and 7.11 × 10 5 , respectively.
  • The set of experiments was used to determine a linear relationship between applied stress and the martensitic start temperature ( M s ) for both materials. This feature allows for martensitic transformation in the model to be triggered consistently with that in the experiments.
  • The solid state phase transformation model parameters for the E410NiMo filler material were also determined following the same methodology used in the author’s previous work [25] for S41500.
  • When fed with the material parameters found in this study, the Leblond model has been successfully compared with the experiments. The model showed very good agreement with the experiments and has proven effective in reproducing the transformation-induced plastic strain behavior of the stress levels used for TRIP coefficient determination. Although more discrepancy is observed at higher magnitudes of stress close to or beyond the austenitic phase yield strength, the model still provides results that are consistent with the experiments.
The presented work, thus, allows for the strains that occur during martensitic transformation in 13Cr-4Ni steels upon cooling with applied external loads to be predicted accurately. It also allows for the simulation of the manufacturing and repairing processes of hydraulic turbines, which will enhance residual stress calculations and fatigue life assessments.

Author Contributions

Conceptualization, J.-B.L. and C.B.; methodology, J.-B.L. and C.B.; validation, J.-B.L.; formal analysis, J.-B.L.; investigation, J.-B.L.; resources, J.L.; data curation, J.-B.L.; writing—original draft preparation, J.-B.L.; writing—review and editing, D.P., J.L., H.C. and J.-B.L.; visualization, J.-B.L.; supervision, J.L. and H.C.; project administration, D.P., J.L. and H.C.; funding acquisition, J.L. and H.C. All authors have read and agreed to the published version of the manuscript.

Funding

The authors want to thank the Natural Sciences and Engineering Research Council of Canada (NSERC) for their financial support of this project.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors also wish to express their gratitude to Mathieu Paquin for his precious contribution to the welding carried out for sample preparation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hensel, J.; Nitschke-Pagel, T.; Tchoffo Ngoula, D.; Beier, H.-T.; Tchuindjang, D.; Zerbst, U. Welding Residual Stresses as Needed for the Prediction of Fatigue Crack Propagation and Fatigue Strength. Eng. Fract. Mech. 2018, 198, 123–141. [Google Scholar] [CrossRef]
  2. Savaria, V.; Bridier, F.; Bocher, P. Predicting the Effects of Material Properties Gradient and Residual Stresses on the Bending Fatigue Strength of Induction Hardened Aeronautical Gears. Int. J. Fatigue 2016, 85, 70–84. [Google Scholar] [CrossRef]
  3. Paquet, D.; Lanteigne, J.; Bernard, M.; Baillargeon, C. Characterizing the Effect of Residual Stresses on High Cycle Fatigue (HCF) with Induction Heating Treated Stainless Steel Specimens. Int. J. Fatigue 2014, 59, 90–101. [Google Scholar] [CrossRef]
  4. Itoh, Y.Z.; Suruga, S.; Kashiwaya, H. Prediction of Fatigue Crack Growth Rate in Welding Residual Stress Field. Eng. Fract. Mech. 1989, 33, 397–407. [Google Scholar] [CrossRef]
  5. Deschênes, P.-A.; Lanteigne, J.; Verreman, Y.; Paquet, D.; Lévesque, J.-B.; Brochu, M. A New Experimental Method to Study the Influence of Welding Residual Stresses on Fatigue Crack Propagation. Int. J. Fatigue 2017, 100, 444–452. [Google Scholar] [CrossRef]
  6. Liljedahl, C.D.M.; Zanellato, O.; Fitzpatrick, M.E.; Lin, J.; Edwards, L. The Effect of Weld Residual Stresses and Their re-Distribution with Crack Growth During Fatigue Under Constant Amplitude Loading. Int. J. Fatigue 2010, 32, 735–743. [Google Scholar] [CrossRef]
  7. Jones, K.W.; Dunn, M.L. Fatigue Crack Growth Through a Residual Stress Field Introduced by Plastic Beam Bending. Fatigue Fract. Eng. Mater. Struct. 2008, 31, 863–875. [Google Scholar] [CrossRef]
  8. Lindgren, L.-E. Numerical Modelling of Welding. Comput. Methods Appl. Mech. Engrg. 2006, 195, 6710–6736. [Google Scholar] [CrossRef]
  9. Yaghi, A.H.; Hyde, T.H.; Becker, A.A.; Sun, W.; Wen, W.; Hilson, G.; Simandjuntak, S.; Flewitt, P.E.J.; Pavier, M.J.; Smith, D.J. Comparison of Measured and Modelled Residual Stresses in a Welded P91 Steel Pipe Undergoing Post Weld Heat Treatment. Int. J. Press. Vessels Pip. 2020, 181, 104076. [Google Scholar] [CrossRef]
  10. Deng, D. FEM Prediction of Welding Residual Stress and Distortion in Carbon Steel Considering Phase Transformation Effects. Mater. Des. 2009, 30, 359–366. [Google Scholar] [CrossRef]
  11. Rikken, M.; Pijpers, R.; Slot, H.; Maljaars, J. A Combined Experimental and Numerical Examination of Welding Residual Stresses. J. Mater. Process. Technol. 2018, 261, 98–106. [Google Scholar] [CrossRef]
  12. Jiang, W.; Chen, W.; Woo, W.; Tu, S.-T.; Zhang, X.-C.; Em, V. Effects of Low-Temperature Transformation and Transformation-Induced Plasticity on Weld Residual Stresses: Numerical Study and Neutron Diffraction Measurement. Mater. Des. 2018, 147, 65–79. [Google Scholar] [CrossRef]
  13. Sun, J.; Hensel, J.; Klassen, J.; Nitschke-Pagel, T.; Dilger, K. Solid-State Phase Transformation and Strain Hardening on the Residual Stresses in S355 Steel Weldments. J. Mater. Process. Technol. 2019, 265, 173–184. [Google Scholar] [CrossRef]
  14. Zhang, K.; Dong, W.; Lu, S. Experimental and Numerical Investigation of Stress and Distortion in AF1410 Steel Under Varying Quenching Conditions. J. Mater. Eng. Perform. 2022, 31, 6858. [Google Scholar] [CrossRef]
  15. Coroas, C.; Viéitez, I.; Martín, E.; Román, M. Numerical Modeling for the Prediction of Microstructure and Mechanical Properties of Quenched Automotive Steel Pieces. Materials 2023, 16, 4111. [Google Scholar] [CrossRef]
  16. Greenwood, G.W.; Johnson, R.H. The Deformation of Metals Under Small Stresses During Phase Transformations. Proc. R. Soc. A 1965, 283, 403–422. [Google Scholar]
  17. Magee, C.L. Transformation Kinetics, Microplasticity and Aging of Martensite in Fe-31NI. Ph.D. Thesis, Carnegie Mellon University, Pittsburgh, PA, USA, 1966. [Google Scholar]
  18. Leblond, J.B.; Devaux, J.; Devaux, J.C. Mathematical Modelling of Transformation Plasticity in Steels I: Case of Ideal-Plastic Phases. Int. J. Plast. 1989, 5, 551–572. [Google Scholar] [CrossRef]
  19. Leblond, J.B. Mathematical Modelling of Transformation Plasticity in Steels I: Coupling with Strain Hardening Phenomena. Int. J. Plast. 1989, 5, 573–591. [Google Scholar] [CrossRef]
  20. Neubert, S.; Pittner, A.; Rethmeier, M. Experimental Determination of Trip-Parameter K for Mild- and High-Strength Low-Alloy Steels and a Super Martensitic Filler Material. SpringerPlus 2016, 5, 754. [Google Scholar] [CrossRef]
  21. Huo, Y.; Liu, Y.; Zuo, S.; Wang, B.; Li, Q. Prediction and experimental verification of phase transformation plasticity. Mater. Today Commun. 2023, 36, 106769. [Google Scholar] [CrossRef]
  22. ASTM A240/A240M-18; Standard Specification for Chromium and Chromium-Nickel Stainless Steel Plate, Sheet, and Strip for Pressure Vessels and for General Applications. ASTM International: West Conshohocken, PA, USA, 2018.
  23. ASME BPVC.II.C-2021; ASME Boiler and Pressure Vessel Code SECTION II Part C Specifications for Welding Rods, Electrodes, and Filler Metals. American Society of Mechanical Engineers: New York, NY, USA, 2021.
  24. ASTM E8/E8M-22; Standard Test Methods for Tension Testing of Metallic Materials. ASTM International: West Conshohocken, PA, USA, 2018.
  25. Lévesque, J.-B.; Lanteigne, J.; Champliaud, H.; Paquet, D. Modeling Solid-State Phase Transformations of 13Cr-4Ni Steels in Welding Heat-Affected Zone. Metall. Mater. Trans. A 2020, 51, 1208–1220. [Google Scholar] [CrossRef]
  26. Johnson, W.A.; Mehl, R.F. Reaction Kinetics in Processes of Nucleation and Growth. Trans. AIME 1939, 135, 416–442. [Google Scholar]
  27. Avrami, M. Kinetics of Phase Change. I General Theory. J. Chem. Phys. 1939, 7, 1103–1112. [Google Scholar] [CrossRef]
  28. Kolmogorov, A.N. On The Static Theory of Metal Crystallization. Izv. Akad. Nauk SSR Ser. Fiz. Mat. Nauk. 1937, 1, 355–360. [Google Scholar]
  29. Barmak, K. A Commentary on: “Reaction Kinetics in Processes of Nucleation and Growth”. Metall. Mater. Trans. A 2010, 41A, 2711–2775. [Google Scholar] [CrossRef]
  30. Kissinger, H.E. Reaction Kinetics in Differential Thermal Analysis. Anal. Chem. 1957, 29, 1702–1706. [Google Scholar] [CrossRef]
  31. Koistinen, D.P.; Marburger, R.E. A General Equation Prescribing the Extent of the Austenite-Martensite Transformation in Pure Iron-Carbon Alloys and Plain Carbon Steels. Acta Mater. 1959, 7, 59–60. [Google Scholar] [CrossRef]
  32. van Bohemen, S.M.C.; Sietsma, J. Effect of Composition on Kinetics of Athermal Martensite Formation in Plain Carbon Steels. Mater. Sci. Technol. 2009, 25, 1009–1012. [Google Scholar] [CrossRef]
  33. Liu, L.; Guo, B. Dilatometric Analysis and Kinetics Research of Martensitic Transformation under a Temperature Gradient and Stress. Materials 2021, 14, 7271. [Google Scholar] [CrossRef]
Figure 1. Micrographs of the materials under study: (a) UNS S41500 base metal; (b) E410NiMo base metal.
Figure 1. Micrographs of the materials under study: (a) UNS S41500 base metal; (b) E410NiMo base metal.
Materials 16 07166 g001
Figure 2. Dilatation as a function of temperature during the heating and cooling of S41500 and E410NiMo soft martensitic stainless steels.
Figure 2. Dilatation as a function of temperature during the heating and cooling of S41500 and E410NiMo soft martensitic stainless steels.
Materials 16 07166 g002
Figure 3. Schematic of the test chamber of the quenching dilatometer. The force, F, compressive or tensile, is applied to the sample through the Load cell. The other end of the sample is fixed to the wall.
Figure 3. Schematic of the test chamber of the quenching dilatometer. The force, F, compressive or tensile, is applied to the sample through the Load cell. The other end of the sample is fixed to the wall.
Materials 16 07166 g003
Figure 4. TRIP experiments: thermal and loading history.
Figure 4. TRIP experiments: thermal and loading history.
Materials 16 07166 g004
Figure 5. Total strain measured during the cooling of the S41500 specimens for various magnitudes of uniaxial stress (compressive and tensile).
Figure 5. Total strain measured during the cooling of the S41500 specimens for various magnitudes of uniaxial stress (compressive and tensile).
Materials 16 07166 g005
Figure 6. Total strain measured during the cooling of the E410NiMo specimens for various magnitudes of uniaxial stress (compressive and tensile).
Figure 6. Total strain measured during the cooling of the E410NiMo specimens for various magnitudes of uniaxial stress (compressive and tensile).
Materials 16 07166 g006
Figure 7. Final transformation-induced strains as a function of the applied stress magnitude for S41500 and E410NiMo.
Figure 7. Final transformation-induced strains as a function of the applied stress magnitude for S41500 and E410NiMo.
Materials 16 07166 g007
Figure 8. Martensitic start transformation as a function of applied stress.
Figure 8. Martensitic start transformation as a function of applied stress.
Materials 16 07166 g008
Figure 9. Comparison of modeled total strains (solid) and experimental total strains data (dashed) for S41500 stainless steel.
Figure 9. Comparison of modeled total strains (solid) and experimental total strains data (dashed) for S41500 stainless steel.
Materials 16 07166 g009
Figure 10. Comparison of modeled total strains (solid) and experimental total strains data (dashed) for E410NiMo stainless steel.
Figure 10. Comparison of modeled total strains (solid) and experimental total strains data (dashed) for E410NiMo stainless steel.
Materials 16 07166 g010
Table 1. Chemical compositions of UNS S41500 stainless steel and E410NiMo filler material.
Table 1. Chemical compositions of UNS S41500 stainless steel and E410NiMo filler material.
Materialwt%CMnPSSiCrNiMo
ASTMStd. [22]<0.050.50–1.00<0.030<0.030<0.6011.5–14.03.5–5.50.50–1.00
S41500Meas.0.0340.680.0180.0010.4412.74.00.57
ASMEStd. [23]<0.06<1.00<0.04<0.03<1.011.0–12.54.0–5.00.40–0.70
E410NiMoMeas.0.0190.370.0100.0080.4911.94.50.62
Table 2. Selected uniaxial loading ratio for the transformation-induced plasticity parameter determination experiments.
Table 2. Selected uniaxial loading ratio for the transformation-induced plasticity parameter determination experiments.
MaterialLoading Ratio σ / σ γ y
S41500−0.73−0.36−0.1800.180.360.73
E410NiMo−0.60−0.30−0.1500.150.300.60
Table 3. Material parameters for thermal strain calculations and phase transformation models.
Table 3. Material parameters for thermal strain calculations and phase transformation models.
MaterialCTE,  ° C 1 JMAK,  α γ K-M,  γ α
α γ E a , kJ mol n K 0 a m ° C 1 T KM , ° C
S415002.2  × 10 7 T  + 0.00117.0  × 10 7 T  +  0.0017 509.90.664.04 e 24 0.0267Equation (4)
E410NiMo3.5  × 10 7 T  + 0.00117.4  × 10 7 T  +  0.0016 468.90.587.31 e 22 0.0280Equation (5)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lévesque, J.-B.; Baillargeon, C.; Paquet, D.; Lanteigne, J.; Champliaud, H. Characterizing and Modeling Transformation-Induced Plasticity in 13Cr-4Ni Welds upon Cooling. Materials 2023, 16, 7166. https://doi.org/10.3390/ma16227166

AMA Style

Lévesque J-B, Baillargeon C, Paquet D, Lanteigne J, Champliaud H. Characterizing and Modeling Transformation-Induced Plasticity in 13Cr-4Ni Welds upon Cooling. Materials. 2023; 16(22):7166. https://doi.org/10.3390/ma16227166

Chicago/Turabian Style

Lévesque, Jean-Benoit, Carlo Baillargeon, Daniel Paquet, Jacques Lanteigne, and Henri Champliaud. 2023. "Characterizing and Modeling Transformation-Induced Plasticity in 13Cr-4Ni Welds upon Cooling" Materials 16, no. 22: 7166. https://doi.org/10.3390/ma16227166

APA Style

Lévesque, J. -B., Baillargeon, C., Paquet, D., Lanteigne, J., & Champliaud, H. (2023). Characterizing and Modeling Transformation-Induced Plasticity in 13Cr-4Ni Welds upon Cooling. Materials, 16(22), 7166. https://doi.org/10.3390/ma16227166

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop