Next Article in Journal
Kinetics of Dry-Batch Grinding in a Laboratory-Scale Ball Mill of Sn–Ta–Nb Minerals from the Penouta Mine (Spain)
Next Article in Special Issue
Process Condition Diagram Predicting Onset of Microdefects and Fracture in Cold Bar Drawing
Previous Article in Journal
Grouping Methods of Cluster Dynamics Model for Precipitation Kinetics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Research of Fracture Toughness of Aged Ferritic-Martensitic Steel

by
Remigijus Janulionis
1,
Gintautas Dundulis
1,2,* and
Albertas Grybėnas
3
1
Laboratory of Nuclear Installation Safety, Lithuanian Energy Institute, Breslaujos str. 3, LT-44403 Kaunas, Lithuania
2
Department of Mechanical Engineering, Faculty of Mechanical Engineering and Design, Kaunas University of Technology, Studentu str. 56, LT-51424 Kaunas, Lithuania
3
Laboratory of Material Research and Testing, Lithuanian Energy Institute, Breslaujos str. 3, LT-44403 Kaunas, Lithuania
*
Author to whom correspondence should be addressed.
Metals 2020, 10(12), 1686; https://doi.org/10.3390/met10121686
Submission received: 30 October 2020 / Revised: 8 December 2020 / Accepted: 12 December 2020 / Published: 17 December 2020
(This article belongs to the Special Issue Fracture and Damage Mechanics of Metals, Steels and Alloys)

Abstract

:
Generally, material properties such as the modulus of elasticity, yield strength or fracture toughness are determined by conducting an experiment. Sometimes experimental determination cannot be done due to specific experimental conditions, lack of testing material and so on. Also, experiments are time consuming and costly. Therefore, there arises the need for alternative determination methods. A numerical method for the fracture toughness determination of steel P91 is suggested in this paper. For this purpose, the universal finite element software ABAQUS was used. The numerical simulation of the C(T) specimen tension test was carried out using non-linear simulation for a conditional load PQ determination, and linear simulation for fracture toughness value KQ determination. The suggested method is validated by comparing numerical and experimental tests results. The secondary aim of the paper is the evaluation of the ageing effect on the fracture toughness of steel P91. Thermal ageing of the steel was carried out in an electric furnace at 650 °C up to 11,000 h. As the numerical results had a good coincidence with experimental data at room temperature, the prediction of fracture toughness at elevated temperature, i.e., 550 °C, using numerical method was carried out.

1. Introduction

For the prediction of the lifetime of construction or its components the ageing of the materials has to be taken into account. Various ageing mechanisms such as corrosion, fatigue and others can significantly reduce material capability to withstand an applied load by reducing its mechanical properties, i.e., tensile strength, yield strength, ductility, etc. Also, small defects, which appear in manufacturing process, can grow into cracks that could cause a component failure. The understanding of how ageing affects material resistance to cracking is important. Fracture toughness (KIC) [1] is the property which describes the ability of material to resist fracture. It has an important role in the fracture mechanics since it characterizes the stress field in the crack tip region, and helps to predict the fracture condition and remaining lifetime of a mechanical component with cracks [2]. It is a material property where different materials have a different value which depends on the environment, i.e., temperature [3] etc. Various standards such as ASTM and ISO describe the procedures for experimental determination of KIC. Sometimes experimental KIC determination cannot be done due to specific experimental conditions, such as lack of testing material and so on. Also, experiments are time consuming and costly. Therefore, there is a need for alternative KIC determination methods. This paper presents the methodology for numerical determination of fracture toughness, which is validated by experimental tests.
The finite element method (FEM) is the most widely used numerical method for fracture toughness and stress intensity factor (SIF) evaluation. Zare et al. [2] used the finite element method to calculate the stress intensity factor in the finite plate with an edge crack under uniaxial stress loading. It was showed that the neural network method should be used to predict the correlation of the SIF and the position of the edge crack along the length of a finite plate. The specimen size effect related to fracture toughness was investigated by Sakib and Adnan [4]. However, Sakib and Adnan presented the methodology for nanoscale specimens. The authors of the papers [5,6] deal with determination of the SIF using the finite element method and compare simulated results with the theoretical solution. The analysis shows that the numerical solution using the finite element method is more suitable than the analytical solution, because the analytical equations are limited in a range of parameters values; some parameters, such as plate length, are not considered and also it is difficult to determine the accurate SIF for cracks in complex structures. According to the good potential of FEM in fracture mechanics applications it was chosen as the numerical method for fracture toughness determination in this paper.
For the research object as-received and aged ferritic-martensitic steel P91 were chosen. Steel P91 was developed based on 1960s steels with 12% Cr content. 12% Cr content steel has a creep problem when it is exposed at higher temperatures and pressures, the conditions typically experienced at power plants, where temperatures can reach 565 °C. Therefore, steel with increased creep resistance, i.e., steel P91, was developed by the Oak Ridge National Laboratory in Tennessee, USA; it consists of 9% Cr and 1% Mo [7]. P91 has a tempered martensite structure after austenitization, air cooling and tempering [8]. At present, steel P91 is widely used in boiler components of ultra-supercritical power plants; also it is a common choice of pipe material for heat exchangers in nuclear power plants [9].
The main aim of this research is to present and validate the numerical method for the fracture toughness determination. The validation of the numerical determination methodology was done by comparing the obtained fracture toughness values with the values determined experimentally. As the secondary aim of the research is to investigate the ageing effect on fracture toughness of steel P91. For this purpose, ageing of the steel was performed in laboratory conditions at 650 °C for different durations and up to 11,000 h. Experimental evaluation of fracture toughness for the validation of the numerical method was done at room temperature. As the numerical results had a good coincidence with experimental results, the prediction of fracture toughness at elevated temperature, i.e., 550 °C, using numerical method was carried out. Steel P91 is being considered as a structural material in the fourth-generation reactors and fusion power plants due to its good thermo-physical properties, excellent stability under irradiation and good creep strength up to 550 °C. An elevated temperature of 550 °C was selected because it is designed to be an operational temperature of coolant in IV generation sodium cooled nuclear reactors [10].

2. Experimental Research of Steel P91

2.1. Material and Long-Term Ageing

Steel plate of 30 mm thickness and made of steel P91 was received for experimental testing. A metal analyser Q4 Tasman was used for determination of chemical composition of the as-received sample. These results are presented in Table 1.
The steel plate was cut into smaller samples. The samples were placed in steel boxes filled with calcium carbonate powder to minimize oxidation and isothermally annealed in an electric furnace at 650 °C. The ageing durations were 790, 4150 and 11,000 h. Later, aged steel samples were used for manufacturing of specimens for tensile and fracture toughness testing. The microstructures of P91 steel and its changes after ageing are presented in Figure 1. The figure shows that thermal ageing does have an influence on the microstructure of steel P91. The precipitates do increase in size and continuously degrade the material. Precipitates were observed at the grain boundaries as well as inside the ferrite grains. Coarsening was also noticed at longer ageing times. During ageing Cr and Mo partly replace Fe atoms in the M23C6 carbide lattice. Cr atoms dominate in the M23C6 carbide lattice, because Cr has a higher diffusion coefficient than Mo [11,12]; this was confirmed in previous research [13] where EDX analysis of extracted residue from steel P91 was done.

2.2. Experimental Research Methodology

An Instron Model 8801 test machine (Instron limited, Buckinghamshire, UK) was used for testing the mechanical properties and fracture toughness of steel P91. Room temperature tensile tests were carried out according to EN ISO 6892-1 [14] requirements and elevated temperature tests were carried out according to EN ISO 6892-2 [15] requirements. The cylindrical specimens with a threaded shoulder were selected for this research. The diameter of specimens for the testing was 8 mm and gauge length was 40 mm. The fracture toughness experiments were carried out according to the ASTM E399 [16] standard. The compact tension C(T) specimens were selected for this purpose. The drawing of manufactured specimens for fracture toughness testing is shown in Figure 2.
All experimental tests results are presented in Section 4.1.

3. Numerical Simulation of Fracture Toughness

3.1. Methodology of the Numerical Simulation of Fracture Toughness

The numerical research of fracture toughness was done by numerically simulating experimental tests according to the ASTM E 399 standard.
To determine fracture toughness KIC, at first, the so-called conditional stress intensity factor value KQ has to be found. Then, if all requirements and conditions described in the ASTM E399 standard are met, it can be stated that KQ = KIC. First of all, load vs. crack mouth opening displacement (CMOD) has to be created for KQ determination. Experimentally, CMOD is measured using a shoulder extensometer (Instron limited, Buckinghamshire, UK) placed in machined 90° V-shaped notches at the end of the specimen (see Figure 2). In numerical modelling, the CMOD can be determined by tracking the displacement of the notch edge. A constructed curve is used for determination of loads PQ and Pmax. Pmax is the load value at the highest load point of the curve. PQ is the load at which the SIF reaches the KQ value. For this purpose, additional lines must be drawn. The OA line is a tangent line to the linear part of the curve (see Figure 3), where point O is at the origin of the curve. Line OP5 is a secant line with a 5% slope to the OA line and found by (P/v)5 = 0.95(P/v)O. In the case of a Type I curve, PQ is found at the intersection point of the curve and the secant line. For Type II and Type III curves, the PQ load is found at the first vertex.
Knowing the PQ load, KQ can be analytically calculated by using the following equation for a standard C(T) specimen [16]:
K Q = ( P Q / B W 1 / 2 ) · f ( a / W ) ,
where:
f ( a / W ) = ( 2 + a W ) ( 1 a W ) 3 / 2 · [ 0.886 + 4.64 a W 13.32 ( a W ) 2 + 14.72 ( a W ) 3 5.6 ( a W ) 4 ] ,
where PQ is the load at which KQ is found, N; B is the thickness of the specimen (see Figure 4), m; W is the width of the specimen (see Figure 4), m; a is the length of the crack (see Figure 4), m.
Kmax is calculated using the same Equations (1) and (2) as for the KQ calculation, only replacing load PQ with load Pmax.
Numerical simulation of fracture toughness was thus done by applying the instructions described above to the finite element model. ABAQUS v6.11 [17] was chosen as the finite element code. The influence of the ageing effect of steel P91 was evaluated by applying material properties of aged steel (see Figure 5 and Table 2) in the FE model.

3.2. Numerical Model

The FE model of a C(T) specimen was created for numerical simulation of fracture toughness using finite element code ABAQUS v6.11 (Dassault Systèmes Simulia Corp, Johnston, RI, USA). The mesh of the model is shown in Figure 4. As the specimen is symmetric, in order to reduce the modelling and computational efforts, only one half of it was created. All dimensions of the numerical model correspond to the dimensions of the specimen used in the experiment (see Figure 2). C3D20R elements [17] were chosen for the FE mesh. These elements are three-dimensional, brick-shaped and they have 20 nodes. Fatigue crack length a1 (see Figure 4) was used for all numerical simulation cases and was equal to 1.8 mm. To describe the material of the model, the mechanical properties, presented in Section 4.1, were applied. The elastic-plastic material behaviour was described by using the stress–strain curves determined experimentally (see Figure 5). The actual values were entered up to the limit of Agt (see Table 2). The Poisson’s coefficient of ν = 0.3 was used in all cases.
Reference points (RPs), which were placed in the centre of the pin hole of the specimen model, was used for the boundary conditions. The displacement along X and Z axis and rotation around the Y axis were restricted for the RP point. However, the translation of the RP point along the Y axis was used as a load. In order to evaluate large strains at the front of the crack the option for nonlinear geometric effects Nlgeom [17] was selected. The numerical model has yielded the following results: reaction force at the RP point, crack mouth opening displacement (CMOD), stress intensity factor and J-integral at the crack front. CMOD was measured by tracking the position of point P1 (see Figure 4) in the Y direction (point P1 is situated on the specimen edge where an extensometer is attached in the experiment).
According to the user’s manual [17] the mesh around the crack tip should be created in a particular shape. The elements should be brick shaped but they have to be laid out in a circular pattern around the crack. It is possible to make the mesh like this by collapsing the element’s face nodes, which are the closest to the crack front, into the single line. That action makes prism shape elements which still have 20 nodes (see Figure 4). These types of elements are called singular elements.

3.3. Numerical Determination of SIF

As described in the methodology (see Section 3.1) tension load vs. CMOD curves were determined after simulation. Non-linear tension simulation with elastic-plastic material properties was used for this purpose. The comparison of experimental and numerical simulation results for not aged material at 20 °C is presented in Figure 6. The figure shows three experimentally determined curves and a curve determined by numerical simulation. As can be seen, experimental results are very close. According to Figure 6 the numerically determined curve quite well matches the curves determined experimentally, where the highest deviation of results is less than 5%.
When the load vs. CMOD curves are determined PQ and Pmax loads can be found. These loads were determined using ASTM E399 procedures [16] as described in Section 3.1. As an example, the determination of the PQ load for not aged steel at room temperature is shown in Figure 7.
According to the ASTM E399 procedures [16] load PQ can be used for KQ determination. For this purpose, the load vs. stress intensity factor curve has to be constructed. The stress intensity factor can be calculated analytically or numerically. To calculate the SIF analytically the Equations (1) and (2) presented in Section 3.1 are used. Numerically, the stress intensity factor can be calculated by running the linear simulation of C(T) specimen tension. The other quite common way to calculate the SIF is by determination of the J-integral and recalculation of the SIF from the relation between J and K under plane strain condition. The stress intensity factor KJ obtained from the J-integral can be found by equation:
K J = E J 1 υ 2 ,
where E is the modulus of elasticity, Pa and ν is Poisson’s ratio.
The comparison of the analytically calculated stress intensity factor with the numerically simulated KJ (see Figure 8) shows that at lower levels of loading KJ is close to the analytically calculated stress intensity factor. However, in our case, it starts to deviate when the load reaches 30 kN. Also, the patterns of the curves are different: the numerically simulated curve is exponential while the analytically calculated curve is linear. It means that this stress intensity factor calculation approach can be used only at low loading levels. And as PQ load is around 40 kN for the as-received and aged steel P91, this approach for numerical simulation of fracture toughness is not suitable.
The other approach is direct calculation of the stress intensity factor by the procedures implemented in ABAQUS [17]. However, as the stress intensity factor is the parameter of linear fracture mechanics [1] it can be determined only by linear FE analysis [17]. Therefore, the second linear elastic specimen tension simulation has to be carried out, where only elastic material properties are used. Figure 9 shows that the direct stress intensity factor calculation does yield close results to the analytical approach. Numerically determined SIF values at load of 50 kN deviate from the analytical results only up to 2%. This small deviation was received due to the fact that the 3D FE model does evaluate strains in three directions and strains along the thickness of the specimen have a considerable effect on the specimen of 25 mm thickness. However, the ASTM E399 standard, where the equation for the analytical SIF calculation is described, presents a test procedure for plane–strain fracture toughness where strains along the thickness should be neglected.

4. Results and Discussion

4.1. Results of Experimental Research

The averaged mechanical properties of the as-received and aged steel P91 determined experimentally are shown in Table 2.
Experimentally determined stress–strain curves of the as-received and aged steel P91 at room and elevated temperatures are presented in Figure 5. These curves were used for the description of material behaviour in the numerical model for fracture toughness simulation.
Experimentally determined fracture toughness results at room temperature are shown in Table 3.

4.2. Results of Numerical Fracture Toughness Simulation

The fracture toughness was investigated after a long-term exposure of P91 steel at 650 °C temperature. To evaluate the ageing impact on the fracture toughness of the P91 steel it was calculated for several specimens with different ageing time. The numerical simulation was carried out at room temperature and the determined results were compared with the experiment data. Numerically and experimentally determined PQ and Pmax loads are presented in Figure 10. According to the data ageing of the steel does have an influence and the biggest effect is for the Pmax load. The highest deviation of numerically simulated results comparing to the experiment data was up to 9% and it was received for the specimen aged for 11,000 h. In general, simulation results are in good agreement with the experiment.
Experimentally and numerically determined KQ and Kmax values at room temperature for as-received and aged materials are presented in Figure 11. According to the figure, with an increase of the ageing time both fracture toughness values are decreasing. Comparing numerical and experimental results it can be seen, that at all ageing time points numerically determined KQ values are lower and Kmax values are higher than experimentally determined values. However, the difference between numerical and experimental results is not large, where the highest deviation is less than 9%. The decrease of fracture toughness after an ageing time of 11,000 h is not significant, where experimentally determined KQ and Kmax values decreased respectively by 4% and 7% and numerically determined values decreased by 6% and 5%.
As the numerical simulation results of fracture toughness at room temperature have a good coincidence with experiment data, the numerical prediction of the fracture toughness at elevated temperature, i.e., 550 °C, was carried out. Accordingly [18] the ageing time of steel can be recalculated into the IV generation nuclear reactor operation time regarding reactor design temperature. Hence, steel ageing time of 11,000 h at the 650 °C temperature does correspond approximately to 60 years of operation in an IV generation nuclear reactor at 550 °C.
Numerical simulation results of fracture toughness dependencies on ageing time at elevated temperature are presented in Figure 12. The influence of temperature in the fracture toughness determination was estimated using material properties determined at the given temperature, i.e., the C(T) specimen behaviour is described by the material properties determined at the elevated temperature. Also, it should be mentioned that the PQ load is determined to have close to a linear behaviour of C(T) specimen (the secant line OP5 only has a 5% slope) and for determination of the PQ load no crack growth or initiation modeling is needed. In order to be sure that numerical prognoses are accurate the experimental test for the as-received steel at elevated temperature have been carried out. The comparison shows that experimental and numerical results are close and do not deviate more than 6%. This small deviation confirms that numerical prediction of fracture toughness at elevated temperature is accurate. According to the obtained results, the decrease of KQ and Kmax values are respectively 18% and 6%. However, as only one experimental test was done at higher temperature, the extended investigation on the difference of fracture mechanism, crack mouth opening and closure or crack mouth deformation could not be done reliably. Therefore, it cannot be stated with high confidence that the presented numerical methodology is suitable for fracture toughness determination at higher temperatures.
The study revealed that ageing process has an influence on mechanical properties as well as on fracture toughness of the steel P91. However, the decrease of KQ values for aged material at room temperature was not big. That shows a good resistance to thermal ageing of steel P91 which was confirmed by Baek et al. [19] and other researchers. Baek et al. estimated that the microhardness and strength of modified 9Cr-1Mo steel after long-term ageing at 600 °C for 50,000 h changed insignificantly. Researches determined that yield strength decreased by 28 MPa (equiv. to 5% reduction) and the ultimate tensile strength decreased by 7 MPa (equivalent to a 1% reduction) by comparing aged steel with that of the as-received one. Zieliński and Dobrzański [20] also confirms that ageing at 600 °C and 650 °C do not result in significant changes in mechanical properties. However, Cumino et al. [21] conducted an experiment for steel P91 where it was aged at 625 and 650 °C temperatures and they estimated a bigger influence to material properties. After experiments the authors concluded that tensile properties of steel P91 reduced by 10% after 1000 h of ageing. The authors claimed that after 3000 h of ageing the impact properties were reduced and brittle fracture was increased. Junek et al. [22] evaluated the mechanical properties and microstructures of similar and dissimilar welds of martensitic steels P91 and P92 after long-term laboratory ageing at elevated temperatures. The research of mechanical properties of base metals, i.e., steel P91 and P92, had similar results, i.e., insignificantly changes in mechanical properties. The ultimate strength of aged metals even up to 30,000 h had no significant decrease (approx. 20 MPa) when the ageing temperature was 650 °C. However, the ageing effect was strong at an ageing temperature of 700 °C. Here ultimate strength of steel P91 decreased by 95 MPa. In our study the ageing effect is evaluated by application of material properties of the aged steel in the FE model. Consistent with the findings that aged material properties (see Table 2) have changed insignificantly the decrease of fracture toughness for aged material was also small. Current research revealed that KQ values at room temperature decreased only up to 4% determined experimentally and 6% determined numerically.

5. Conclusions

This research covers the presentation and validation of a numerical method for the fracture toughness determination and evaluation of ageing influence on the fracture toughness of steel P91.
Numerical simulation was done by running non-linear and linear tension simulations of C(T) specimen. Non-linear simulation was used for PQ and Pmax load determination while linear simulation was used for determination of KQ and Kmax. To validate the numerical 3D FE model for fracture toughness analysis the numerical simulation results were compared to the experimental data.
Experimental tests showed that after long time exposure at 650 °C the mechanical properties of steel P91 decreased. However, the change of tensile and yield strength after ageing for 11,000 h was relatively small and did not exceed 5%. Fracture toughness test results showed that KQ and Kmax values also had a small decrease. The highest decrease was for the Kmax value and it reached 7%.
Numerical simulation of ageing influence on fracture toughness of steel P91 was carried out at room and elevated temperature (550 °C). The comparison of the numerically determined fracture toughness results with the experimental results show a good coincidence which demonstrates the suitability of the numerical simulation method for fracture toughness determination at room temperature.
Experimental testing and numerical simulation results obtained in this study show that steel P91 has good resistance to thermal ageing and still has high ductility even after long thermal exposure (11,000 h at 650 °C).

Author Contributions

Conceptualization, G.D. and R.J.; methodology, G.D., R.J. and A.G.; experiment, A.G.; numerical simulation, R.J.; data curation, A.G.; result analysis, G.D., R.J. and A.G., writing—original draft, G.D., R.J. and A.G.; writing—review & editing, G.D., R.J. and A.G.; visualization, R.J.; supervision, G.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Anderson, T.L. Fracture Mechanics, Fundamentals and Applications; CRC Press Inc.: Boca Raton, FL, USA; Ann Arbor, MI, USA; Boston, MA, USA, 1991; 793p, ISBN 0-8493-4277-5. [Google Scholar]
  2. Zare, A.; Kosari, E.S.M.; Asadi, I.; Bigham, A.; Bigham, Y. Finite element method analysis of stress intensity factor in different edge crack positions, and predicting their correlation using neural network method. Res. J. Recent Sci. 2014, 3, 69–73. [Google Scholar]
  3. Peiyuan, L.; Li, C.; Xiaojun, Y.; Dawei, H.; Xiaoyu, Q.; Xiaoyong, Z. A temperature-dependent model for predicting the fracture toughness of superalloys at elevated temperature. Theor. Appl. Fract. Mec. 2018, 93, 311–318. [Google Scholar] [CrossRef]
  4. Sakib, A.R.N.; Adnan, A. On the size-dependent critical stress intensity factor of confined brittle nanofilms. Eng. Fract. Mech. 2012, 86, 13–22. [Google Scholar] [CrossRef]
  5. Mohsin, N.R. Comparison between theoretical and numerical solutions for center, single edge and double edge cracked finite element plate subjected to tension stress. Int. J. Mech. Prod. Eng. Res. Dev. 2015, 5, 11–20. [Google Scholar]
  6. Kacianauskas, R.; Zenon, M.; Zarnovskij, V.; Stupak, E. Three-dimensional correction of the stress intensity factor for plate with a notch. Int. J. Fract. 2005, 136, 75–98. [Google Scholar] [CrossRef]
  7. Silva, F.J.G.; Pinho, A.P.; Pereira, A.B.; Paiva, O.C. Evaluation of welded joints in P91 steel under different heat-treatment conditions. Metals 2020, 10, 99. [Google Scholar] [CrossRef] [Green Version]
  8. Blinn, B.; Görzen, D.; Fischer, T.; Kuhn, B.; Beck, T. Analysis of the thermomechanical fatigue behavior of fully ferritic high chromium steel Crofer®22 H with cyclic indentation testing. Appl. Sci. 2020, 10, 6461. [Google Scholar] [CrossRef]
  9. Yang, K.; Zhang, Y.; Zhao, J. Elastoplastic fracture analysis of the P91 steel welded joint under repair welding thermal shock based on XFEM. Metals 2020, 10, 1285. [Google Scholar] [CrossRef]
  10. Bertrand, F.; Marie, N.; Bachrata, A.; Brun-Magaud, V.; Droin, J.B.; Manchon, X.; Herbreteau, K.; Farges, B.; Carluec, B.; Poumerouly, S.; et al. Status of severe accident studies at the end of the conceptual design of ASTRID: Feedback on mitigation features. Nucl. Eng. Des. 2018, 326, 55–64. [Google Scholar] [CrossRef] [Green Version]
  11. Sasikala, G.; Ray, S.K.; Mannan, S.L. Kinetics of transformation of delta ferrite during creep in a type 316(N) stainless steel weld metal. Mater. Sci. Eng. A 2003, 359, 86–90. [Google Scholar] [CrossRef]
  12. Thomas Paul, V.; Saroja, S.; Vijayalakshmi, M. Microstructural stability of modified 9Cr–1Mo steel during long term exposures at elevated temperatures. J. Nucl. Mater. 2008, 378, 273–281. [Google Scholar] [CrossRef]
  13. Baltušnikas, A.; Lukošiūtė, I.; Makarevičius, V.; Kriūkienė, R.; Grybėnas, A. Influence of thermal exposure on structural changes of M23C6 carbide in P91 steel. J. Mater. Eng. Perform. 2016, 25, 1945–1951. [Google Scholar] [CrossRef]
  14. ISO 6892-1 Metallic Materials–Tensile Testing–Part 1: Method of Test at Room Temperature; European Committee for Standardization: Brussels, Belgium, 2009.
  15. ISO 6892-2 Metallic Materials–Tensile Testing–Part 2: Method of Test at Elevated Temperature; European Committee for Standardization: Brussels, Belgium, 2011.
  16. ASTM E399-12e3 Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness KIc of Metallic Materials; ASTM International: West Conshohocken, PA, USA, 2012.
  17. ABAQUS/Standard User’s Manual, Version 6.11. 2011. Available online: http://130.149.89.49:2080/v6.11/index.html (accessed on 10 December 2020).
  18. Grybėnas, A.; Makarevičius, V.; Baltušnikas, A.; Lukošiūtė, I.; Kriūkienė, R. Correlation between structural changes of M23C6 carbide and mechanical behaviour of P91 steel after thermal aging. Mat. Sci. Eng. A Struct. 2017, 696, 453–460. [Google Scholar] [CrossRef]
  19. Baek, J.H.; Kim, S.H.; Lee, C.B.; Hahn, D.H. Mechanical properties and microstructural evolution of modified 9Cr-1Mo steel after long-term ageing for 50,000 h. Met. Mater. Int. 2009, 15, 565–573. [Google Scholar] [CrossRef]
  20. Zieliński, A.; Dobrzański, J. Characteristics of changes in properties and structure of X10CrMoVNb9-1 steel due to long-term impact of temperature and stress. Arch. Mater. Sci. 2013, 60, 72–78. [Google Scholar]
  21. Cumino, G.; Di Cuonzo, S.; Di Gianfrancesco, A.; Tassa, O. Advanced high chromium ferritic steels for boiler components operating at high temperature. Lat. Am. Appl. Res. 2002, 32, 229–235. [Google Scholar]
  22. Junek, M.; Svobodova, M.; Janovec, J.; Horvath, J.; Duchacek, P. Long-term thermal degradation of narrow gap orbital welded P91 and P92 steels. Int. J. Press. Vessels Pip. 2020, 185, 104133. [Google Scholar] [CrossRef]
Figure 1. SEM (CARL ZEISS MICROSCOPY LIMITED, Cambridge, UK) microstructure of the P91 steel in the as-received condition (a) and after 11,000 h of ageing at 650 °C (b).
Figure 1. SEM (CARL ZEISS MICROSCOPY LIMITED, Cambridge, UK) microstructure of the P91 steel in the as-received condition (a) and after 11,000 h of ageing at 650 °C (b).
Metals 10 01686 g001
Figure 2. C(T) specimen drawing (dimensions are in mm), according to the ASTM E 399 standard [16].
Figure 2. C(T) specimen drawing (dimensions are in mm), according to the ASTM E 399 standard [16].
Metals 10 01686 g002
Figure 3. Principal types of force-displacement (CMOD) records (ASTM E399-12e3 2012).
Figure 3. Principal types of force-displacement (CMOD) records (ASTM E399-12e3 2012).
Metals 10 01686 g003
Figure 4. Finite element mesh of the C(T) specimen.
Figure 4. Finite element mesh of the C(T) specimen.
Metals 10 01686 g004
Figure 5. Stress–strain curves of the as-received and aged steel P91 at room and elevated temperatures.
Figure 5. Stress–strain curves of the as-received and aged steel P91 at room and elevated temperatures.
Metals 10 01686 g005
Figure 6. Experimentally and numerically determined tension load vs. CMOD curves for as-received steel P91 at room temperature.
Figure 6. Experimentally and numerically determined tension load vs. CMOD curves for as-received steel P91 at room temperature.
Metals 10 01686 g006
Figure 7. Determination of PQ load for the as-received specimen at room temperature.
Figure 7. Determination of PQ load for the as-received specimen at room temperature.
Metals 10 01686 g007
Figure 8. Comparison of KJ with the analytically calculated stress intensity factor.
Figure 8. Comparison of KJ with the analytically calculated stress intensity factor.
Metals 10 01686 g008
Figure 9. Stress intensity factor vs. load.
Figure 9. Stress intensity factor vs. load.
Metals 10 01686 g009
Figure 10. Numerically and experimentally determined PQ and Pmax loads of steel P91 at different ageing times.
Figure 10. Numerically and experimentally determined PQ and Pmax loads of steel P91 at different ageing times.
Metals 10 01686 g010
Figure 11. Numerically and experimentally determined fracture toughness values vs. ageing time of steel P91 at room temperature.
Figure 11. Numerically and experimentally determined fracture toughness values vs. ageing time of steel P91 at room temperature.
Metals 10 01686 g011
Figure 12. Numerically predicted fracture toughness values vs. ageing time of steel P91 at elevated temperature.
Figure 12. Numerically predicted fracture toughness values vs. ageing time of steel P91 at elevated temperature.
Metals 10 01686 g012
Table 1. Chemical composition of the as-received P91 steel plate (wt.%).
Table 1. Chemical composition of the as-received P91 steel plate (wt.%).
CSiMnPSNAlCrMoCuNbNiV
0.0940.380.350.006<0.00050.0350.0128.780.960.070.080.160.18
Table 2. Mechanical properties of the P91 steel plate.
Table 2. Mechanical properties of the P91 steel plate.
Test Temperature (°C)Ageing Time at 650 °C (h)Rp (MPa)Rm (MPa)Agt (%)At (%)E (GPa)
Room temperature05136669.823214
7905076589.923212
41505016519.522212
11,0004876379.221209
55003503852.027166
7903503882.527165
41503493812.028165
11,0003373672.028162
Table 3. Averaged fracture toughness test results at room temperature for the as-received and aged P91 steel.
Table 3. Averaged fracture toughness test results at room temperature for the as-received and aged P91 steel.
Ageing at 650 °C (h)PQ (kN)Pmax (kN)KQ (MPa√m)Kmax (MPa√m)
042.482.775.1146.4
79041.481.373.4143.5
415041.978.073.9137.5
11,00041.177.070.6135.8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Janulionis, R.; Dundulis, G.; Grybėnas, A. Numerical Research of Fracture Toughness of Aged Ferritic-Martensitic Steel. Metals 2020, 10, 1686. https://doi.org/10.3390/met10121686

AMA Style

Janulionis R, Dundulis G, Grybėnas A. Numerical Research of Fracture Toughness of Aged Ferritic-Martensitic Steel. Metals. 2020; 10(12):1686. https://doi.org/10.3390/met10121686

Chicago/Turabian Style

Janulionis, Remigijus, Gintautas Dundulis, and Albertas Grybėnas. 2020. "Numerical Research of Fracture Toughness of Aged Ferritic-Martensitic Steel" Metals 10, no. 12: 1686. https://doi.org/10.3390/met10121686

APA Style

Janulionis, R., Dundulis, G., & Grybėnas, A. (2020). Numerical Research of Fracture Toughness of Aged Ferritic-Martensitic Steel. Metals, 10(12), 1686. https://doi.org/10.3390/met10121686

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