Next Article in Journal
Research of Resistance of Selected Materials to Abrasive Wear to Increase the Ploughshare Lifetime
Next Article in Special Issue
Failure Rate Model of Materials under Uncertain Constant Amplitude Cyclic Load
Previous Article in Journal
Use of Low Melting Point Metals and Alloys (Tm < 420 °C) as Phase Change Materials: A Review
Previous Article in Special Issue
In-Situ Thermography Investigation of Crack Growth in Armco Iron under Gigacycle Fatigue Loading
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simplified Elastoplastic Fatigue Correction Factor Analysis Approach Based on Minimum Conservative Margin

Science and Technology on Reactor System Design Technology Laboratory, Nuclear Power Institute of China, Chengdu 610213, China
*
Author to whom correspondence should be addressed.
Metals 2022, 12(6), 943; https://doi.org/10.3390/met12060943
Submission received: 28 April 2022 / Revised: 17 May 2022 / Accepted: 23 May 2022 / Published: 30 May 2022
(This article belongs to the Special Issue Fatigue Design of Steel and Composite Structures)

Abstract

:
ASME and RCC-M codes specify an elastoplastic fatigue analysis technique: a simplified elastoplastic fatigue analysis method based on linear elastic analysis. In this method, the elastic strain range is multiplied by the elastoplastic correction factor (Ke) to envelope the actual plastic strain range for fatigue evaluation. The ASME or RCC-M provide the Ke parameters of typical materials, such as austenitic stainless steel and low alloy steel. However, how can the parameters of the material not included in the codes be determined? Based on the existing material Z2CND18.12 (nitrogen control) in the codes and taking into account various sensitive factors, the minimum conservative margin of Ke for this material is calculated, and then the parameters of nonstandard materials are determined iteratively based on the conservative margin. The sensitive factors include the different structure model, load types, the loading control mode, temperature value and the material constitutive model. Based this approach, the Ke parameters of TA16 are determined and verified by the transient with drastic change in temperature and pressure. The results of the case show that the simplified elastoplastic fatigue analysis can envelope the results of cyclic plastic fatigue analysis. The minimum margin approach established in this paper can reasonably determine the Ke value of materials beyond the ASME and RCC-M codes.

1. Introduction

The simplified elastoplastic fatigue analysis method is mainly used in the fatigue calculation of nuclear codes ASME [1] and RCC-M [2]. The main purpose of this method is to simplify the calculation process, multiplying the elastic strain range by the elastoplastic correction factor Ke [3] to envelope the actual strain range. It is adopted to prevent fatigue and fracture failure of the components in actual operation [4,5,6,7]. Because the simplified elastoplastic fatigue analysis has the characteristic of high efficiency, it is widely used in engineering [8,9,10]. A number of Ke factors from nuclear and non-nuclear codes is also evaluated and briefly discussed in reference [11,12,13], and an alternative improved plasticity correction method is proposed. Different from the ASME code, for austenitic stainless steel, RCC-M and ASME Code Case N-779 [14] explicitly distinguish between stresses arising due to mechanical and thermal loads.
The ASME or RCC-M provide Ke expressions and parameters for typical materials, such as austenitic stainless steel and low alloy steel. For determining the Ke of new material, the experimental method is too expensive. Shao et al. [15,16,17] proposed a numerical determination approach of Ke, but it did not give the determination process of the minimum margin of materials in the codes in detail, including the influence of sensitive factors on Ke.
In this paper, a numerical calculation method for determining the Ke parameters of new materials is established, based on the Ke parameters of austenitic stainless steel Z2CND18.12 (nitrogen control) given in the code, and taking into account various sensitive factors. Titanium alloy TA16, as a new material, needs to be used in equipment, but lacks Ke parameters. The uniaxial tension, cyclic strain and stress experiments of TA16 are conducted to establish the constitutive models. The Ke expression and parameters of TA16 are determined by the Ke numerical calculation method. Then, the new Ke parameters of TA16 are applied to the simplified fatigue assessment of the feedwater nozzle. By comparing the simplified elastoplastic fatigue analysis and cyclic-plastic fatigue analysis of the structure, it is found that the elastic strain is larger than the elastoplastic strain after Ke modification. The correctness and conservativeness of the Ke of TA16 are verified. Accordingly, a simplified elastoplastic fatigue analysis technology of titanium alloy material is developed.

2. Analysis Method of Ke

2.1. Ke in ASME and RCC-M

For the case of mechanical load, ASME [1] and RCC-M [2] calculate Ke according to the range of primary plus secondary stress Se and allowable stress intensity Sm by the formulas as follows:
K e ( m e c h ) = 1.0     f o r   S n   3 S m 1.0 + ( 1 n ) n ( m 1 ) S n 3 S m 1 1 / n     f o r   S n   3 m S m f o r   3 S m < S n < 3 m S m
where n is the material hardening index, and m is the ratio when the value of Sn/3Sm reaches the maximum. The values for m and n are given in ASME and RCC-M for low alloy steel, martensitic stainless steel, carbon steel, austenitic stainless steel, and nickel-chromium-iron alloy.
Different from the ASME, for austenitic stainless steel, RCC-M [2] adopts the following formula to calculate the elastoplastic correction factor under thermal load:
K e ( t h e r m ) = max 1.86 1 1 1.66   +   ( S n / S m ) 1
S P = K e ( m e c h ) S P ( m e c h ) + K e ( t h e r ) S P ( t h e r m )
Ke(RCCM) can be obtained according to Equation (3), based the Sp(mech) and Sp(therm), where Sp is the range of total stresses. The RCC-M code provides the parameters for some typical materials, such as austenitic stainless steel, low alloy steel, and so on. It is known from Equations (1) and (2) that the value of Ke(mech) is determined by parameters m and n, and the value of Ke(therm) is determined by parameters 1.86 and 1.66, defined as A and B, respectively.

2.2. Numerical Method for Determining Ke Parameters

Because the same expression of Ke(mech) is used in RCC-M and ASME codes, the expression of elastoplastic correction factor in RCC-M and ASME will also be assumed for the new material. The Ke(thermal) will also adopt the same expression as that of RCC-M austenitic stainless steel. Once the values of A, B, m, and n of the new materials are determined, the Ke expression can be established.
The numerical method needs to consider the following factors for the Ke:
(1)
Influence of material strength and constitutive model
The experimental error will affect the yield stress strength and the allowable stress. Equation (1) shows that Ke is related to Sm, so the fluctuation in the strength value will affect the determination of Ke.
Different constitutive models can be used in finite element elastoplastic calculation. The multilinear elastoplastic constitutive model and nonlinear follow-up hardening Chaboche model are adopted. Table 1 shows the stress and strain relationship of Z2CND18.12 (nitrogen control) at different temperatures, which can be used in the multilinear elastoplastic constitutive model. These properties come from Section 3—Tom 1, subsection z-appendix A3.3s of the RCC-MRx code [18].
(2)
Load effect
(a)
Load type: thermal load; mechanical load; mixed load (thermal load plus mechanical load)
(b)
Loading control mode: stress control and displacement control
(c)
Temperature: two typical temperatures (T = 30 °C and T = 350 °C)
In order to consider the influence of different load types, loading control mode, and temperature, the working condition parameters shown in Table 2 are set for the finite element calculation.
(3)
Geometry
According to Equation (1), it is necessary to make the calculated primary plus secondary stress range distributed in the interval of 3Sm~6Sm, and some adjustment methods such as the stress distribution and load value are adopted.
As shown in Figure 1, by adjusting the relevant size parameters, the stress concentration of the model can be changed. For the notched sheet, R is notch radius, and the stress distribution at the notched sheet will change if the R value is adjusted. The Model-II is a typical nozzle, and L3 is the length of the thickness region. The L3 can control the curvature of geometrical discontinuities, so as to change the stress distribution.
Figure 2 shows the x-direction strain nephogram of the finite element calculation of the notched sheet. It can be seen from Figure 2 that the maximum x-direction strain is located at node 125 or 127 at the center of the notch, whether it is elastic analysis or elastoplastic analysis. The stress and strain values of these two nodes will be taken as the discussion object in the subsequent finite element calculation.
Figure 3 and Figure 4 show the finite element calculation results of the nozzle under elastic analysis and elastoplastic analysis under the displacement control and load control, respectively. Among them, the load and displacement are applied at the right end. A pressure load and thermal load are applied to the inner surface of the nozzle together. It can be seen from the figures that the maximum deformation occurs at the transition section of node 409 and 298, marked in Figure 1, and the stress-strain of these nodes will be taken as the discussion object in the subsequent mechanical finite element calculation.
(4)
Loading waveform (Ratio)
A different ratio to control the loading waveform is selected (see Figure 5), such as in the four cases where the ratio is 0, −0.5, −0.8 and −1. The ratio is equal to the ratio of the minimum load to the maximum load during the cycle. Even under the same load, different results can be obtained by adjusting the ratio.
Figure 6 shows the stress-strain curves of different stress peak and valley ratios of 409 nodes simulated by the Chaboche model. It can be seen from the figure that the cyclic stress-strain curve is closely related to the ratio, that is, with the increase in ratio, the cyclic strain (strain rate) increases. However, no matter what the ratio is, the cyclic stress-strain curve is stable after a certain number of cycles (for example, when the number of cycles is 10, as shown in Figure 6c). Therefore, the number of cycles is controlled between 10 and 15 when the Chaboche model is used for finite element calculation later.
In this part, the notched sheet is taken as the representative to give detailed analysis results. Figure 7 shows the influence of the loading control mode on the elastoplastic correction factor at T = 20 °C and T = 350 °C, respectively. As shown in the figure:
(1) When T = 20 °C, there is little difference between the Ke calculated by the two loading control methods;
(2) When T = 350 °C, the calculated Ke of the load control loading mode is higher than that of the displacement control loading mode;
(3) The elastoplastic correction factors calculated under the two loading control modes are below the Ke (RCC-M).
Figure 7. Ke with Se under different loading control modes with notched sheet model.
Figure 7. Ke with Se under different loading control modes with notched sheet model.
Metals 12 00943 g007
Figure 8 shows the Ke with Se calculated by the multilinear elastoplastic constitutive model and the nonlinear follow-up hardening model at T = 20 °C and T = 350 °C, respectively. It can be seen from the figure that under the two temperatures, the Ke calculated by the Chaboche model is significantly lower than that calculated by the EP model.
Figure 9 shows the Ke with Se calculated at two temperatures (T = 20 °C and T = 350 °C). It can be seen from the figure that temperature is the key factor affecting the elastoplastic correction factor. Under the same Se, the Ke value calculated at high temperature is much higher than that calculated at room temperature.
Based on the two typical structures, the effects of loading control mode, constitutive model, temperature, and load type on the elastoplastic correction factor Ke are calculated. The results of the nozzle model under thermal load are listed in Table 3. The results show that the above factors affect the value of Ke, and Ke(EP) is completely enveloped by Ke (RCC-M).
Finally, the conservative margin of Z2CND18.12 (nitrogen control) under various sensitive factors can be determined as follows:
η = min K e ( R C C M ) K e ( E P ) K e ( R C C M ) 1 ,   when   S n > 3 m S m
K e ( R C C M ) is obtained by Equations (1)–(3), and K e ( E P ) is obtained by dividing the elastoplastic strain range Δεep by the elastic strain range Δεee. Based on a large number of calculations under various factors, Table 3 and Table 4 give detailed calculation results of the condition where the minimum margin is located. It can be seen from Table 3 that the minimum conservative margin of Ke under thermal load is 11.3%. The minimum conservative margin under mechanical load is 14.4% in Figure 9, and the minimum conservative margin under thermal mechanical coupling is 7.7% in Table 4. According to the previous description, the flowchart of the numerical method for determining the Ke parameters is shown in Figure 10.
Table 4. Results of Ke under thermal (ΔT = 280 °C) and pressure load of Z2CND18.12 (nitrogen control) for nozzle model.
Table 4. Results of Ke under thermal (ΔT = 280 °C) and pressure load of Z2CND18.12 (nitrogen control) for nozzle model.
NodeP
(MPa)
Sn
(MPa)
Sp
(Therm)
Sp
(Mech)
Δee
(%)
Δep
(%)
Ke
(EP)
Ke
(Therm)
Ke
(Mech)
Ke
(Mixed)
η
(%)
1576−30/0609.51073.9120.00.550.891.621.593.331.7719.58
−38/0616.11073.9152.00.560.931.661.593.331.8117.96
−60/0646.01073.9200.00.580.981.691.603.331.8821.33
−70/0720.11073.9280.00.611.081.761.633.331.9822.20
−80/0775.71073.9340.00.641.251.941.643.332.0510.02
298−30/0655.01088.8120.60.560.911.621.613.331.7820.87
−38/0661.31088.8152.70.560.941.671.613.331.8218.96
−60/0675.41088.8200.90.580.991.701.613.331.8820.02
−70/0749.91088.8281.30.621.101.781.633.331.9820.75
−85/0805.91088.8341.60.651.271.971.653.332.057.66
Figure 10. The flowchart of numerical method for determining Ke parameters.
Figure 10. The flowchart of numerical method for determining Ke parameters.
Metals 12 00943 g010

3. Experiments and Constitutive Model

3.1. TA16 Axial and Cyclic Experiments

Firstly, the monotonic tension test, cyclic strain test, and stress test of TA16 [19] thin-walled circular tube at room and high temperatures are carried out. The chemical composition of TA16 is the same as in reference [20].
The monotonic tensile test adopts displacement control and the loading rate is 0.02 mm/s. The cyclic test is carried out under two control modes: strain control and stress control. The loading strain rate of the strain control is 0.2%/s, and the loading rate of the stress control is 100 MPa/s. The test equipment is MTS809. Figure 11 shows the uniaxial tensile stress-strain curves of TA16 at 30 °C and 350 °C. Figure 12a,c show the cyclic stress-strain curves of TA16 at different strain amplitudes at 30 °C and 350 °C, showing obvious cyclic stability characteristics. Figure 12b,d show the stable cyclic stress-strain curve of TA16. By comparison with the uniaxial tensile stress-strain curve, it is found that the uniaxial tensile curve of TA16 at 350 °C is slightly lower than that of the stable cyclic stress-strain curve.
Figure 13 shows the cyclic stress-strain curves of TA16 under stress control. It can be seen from these figures that TA16 has obvious ratcheting behavior. When TA16 is applied to the steam generator, it will bear a cyclic thermal load, which will inevitably lead to the production of ratcheting behavior. The cyclic accumulation of plastic strain occurs in the high stress region of the material, which will lead to structural deformation and damage, and then affect the service life of structures. Therefore, the ratcheting behavior must be reasonably considered in the structural design of titanium alloys.

3.2. Cyclic Constitutive Model

The Chaboche constitutive model can describe the dynamic hardening and deformation behavior of materials in the stress-strain cycle; the back-stress evolution equation is as follows:
X = i = 1 M X i
X ˙ i = 2 3 C i ε ˙ p γ i X i p ˙   ( i = 1 , 2 )
The stress-strain curves under the strain control and stable cycling at different temperatures are fitted by the least square method. The parameters used in the Chaboche model are shown in Table 5, and the determination method of parameters can be found in [21].

4. The Ke for New Material TA16

4.1. Determining Ke Parameters

Assuming the initial values of m and n for titanium alloy TA16, there are two suggestions for reference:
(1) In the RCC-M, the n value of low alloy steel, carbon steel, austenitic steel, and martensitic steel is 0.2~0.3 and the m value is 1.7~3.0.
(2) The value of m can be calculated according to when the elastic strain reaches 0.3%; n can be obtained by fitting the tensile stress-strain curve of the TA16 with the constitutive relationship σ = A ε n .
According to the above principles, for titanium alloy TA16, the given initial value is m = 2.0 and n = 0.2. With reference to austenitic stainless steel, the initial values of A and B are A = 1.86 and B = 1.66, respectively.
The elastic and elastoplastic calculations under different loads are carried out by using the constitutive parameters of TA16 and considering the geometric model and load type in Section 2. The minimum conservative margin of TA16 under thermal load is 76.1%; under mechanical load, it is 64.7%; and under thermal mechanical coupling, it is 52.1%. These are much larger than the minimum conservative margin of austenitic stainless steel. It shows that the initial parameters are too conservative for calculating the elastoplastic correction factor of TA16, so the correlation coefficients A, B and m, n must be optimized.
The parameter values are adjusted repeatedly. Finally, the elastoplastic correction factor is recalculated to obtain the minimum conservative margin shown in Figure 14, which is close to and slightly greater than that of Z2CND18.12 (nitrogen control), indicating that the parameters of TA16 meet the requirements. Considering the fluctuation increase in TA16 material parameters, a conservative amount of 4.4% is added to the minimum conservative margin of Z2CND18.12 (nitrogen control). The final values are A = 1.37, B = 1.26, m = 2.0 and n = 0.25.

4.2. Verification Calculation

Because of the drastic and frequent temperature transient of the steam generator, the stress exceeds the yield limit increases in some parts, such as the feed water nozzle. In this section, the feed water nozzle model is used to verify the conservativeness and correctness of the elastoplastic correction factor of TA16, and to realize the engineering application of the alloy elastoplastic correction factor. The transient temperature and pressure of the structure are calculated step by step. In the cyclic elastoplastic fatigue analysis, the transient combination with the greatest fatigue usage factor is selected.
The maximum principal strain range and the maximum alternating stress of the structure are calculated by the constitutive model in Section 1, and then the fatigue usage factor of the transient combination is obtained by using the fatigue curve. The simplified elastoplastic calculation results of nodes N1 and N2 under the combination of transients are shown in Table 6. The results of cyclic plastic fatigue analysis are shown in Table 7. From the comparison of the fatigue usage factor U in Table 6 and Table 7, it can be seen that the fatigue usage factor of elastoplastic fatigue analysis is smaller than that of simplified elastoplastic analysis, and the elastoplastic correction factor calculated and optimized by the numerical analysis method can be reasonably applied to simplified elastoplastic fatigue analysis.
Based on the determination of Ke for TA16 and the elastoplastic constitutive model, the simplified elastoplastic strain analysis technology is improved. The engineering application of Ke for TA16 is realized, which makes the analysis method provided by the code become implemented in structure design.

5. Conclusions

(1)
Considering many kinds of sensitive factors, including different structure models, load types, loading control modes, temperature values and material constitutive models, the minimum conservative margin of Ke for Z2CND18.12 (nitrogen control) is 7.7%;
(2)
According to the minimum conservative margin of Z2CND18.12 (nitrogen control), the elastic and elastoplastic calculations for TA16 are carried out, and the Ke of TA16 is determined as A = 1.37, B = 1.26, m = 2.0 and n = 0.25;
(3)
For the feed water nozzle with dramatic temperature transient, the fatigue usage factor of the elastoplastic fatigue analysis is smaller than that of simplified elastoplastic analysis. It verifies the correctness and conservativeness of Ke for titanium alloy TA16.

Author Contributions

Methodology, J.T. and H.L.; supervision, X.F.; writing—original draft, X.S.; writing—review and editing, J.D. and F.X.; data curation, X.L.; validation, H.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (12102416).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

AThe parameter of Ke(therm)
BThe parameter of Ke(therm)
CiChaboche constitutive model parameter
KeElastoplastic correction factor
mThe parameter of Ke
nThe material hardening index, Ke parameter
SnThe range of primary plus secondary stress
SpThe range of total stresses
SmThe allowable stress intensity
γiChaboche constitutive model parameter
εStrain
ε ˙ p Plastic strain rate
σStress

References

  1. ASME. Asme Boiler and Pressure Vessel Code: Section III: Division 1 E Subsection Nb: Class 1 Components Rules for Construction of Nuclear Power Plant Components; ASME: New York, NY, USA, 2010. [Google Scholar]
  2. RCC-M. Design and Construction Rules for Mechanical Components of PWR Nuclear Islands, Section I: Subsection B, Class 1 Components; AFCEN: Paris, France, 2007. [Google Scholar]
  3. Hubel, H. Simplified Elasto-Plastic Fatigue Analysis of Smooth Structures; European Commission: Madrid, Spain, 1995. [Google Scholar]
  4. Branco, R.; Costa, J.; Borrego, L.; Berto, F.; Razavi, S.; Macek, W. Comparison of different one-parameter damage laws and local stress-strain approaches in multiaxial fatigue life assessment of notched components. Int. J. Fatigue 2021, 151, 106405. [Google Scholar] [CrossRef]
  5. Cao, Y.; Nie, W.; Yu, J.; Tanaka, K. A novel method for failure analysis based on three-dimensional analysis of fracture surfaces. Eng. Fail. Anal. 2014, 44, 74–84. [Google Scholar] [CrossRef]
  6. Kamal, M.; Rahman, M.M. Advances in fatigue life modeling: A review. Renew. Sustain. Energy Rev. 2018, 82, 940–949. [Google Scholar] [CrossRef]
  7. Macek, W.; Branco, R.; Korpyś, M.; Łagoda, T. Fractal dimension for bending–torsion fatigue fracture characterisation. Measurement 2021, 184, 109910. [Google Scholar] [CrossRef]
  8. Fu, X.; Wang, D. Plastic correction in simplified elasto-plastic fatigue analysis of nuclear components. Mech. Eng. 2017, 1, 67–70. [Google Scholar]
  9. Zhang, J.; Li, W.; Dai, H.; Liu, N.; Lin, J. Study on the elasticoplastic correlation of low-cycle fatigue for variable asymmetric loadings. Materials 2020, 13, 2451. [Google Scholar] [CrossRef] [PubMed]
  10. Fu, X.L.; Zhang, L.P.; Tian, J.; Shao, X.J.; Du, J.; Zhang, Y.; Liu, Z.Y. Plastic Strain Correction in Fatigue Analysis of Nuclear Power Plant Equipment. Strength Mater. 2021, 53, 183–188. [Google Scholar] [CrossRef]
  11. Clarkson, D.M.; Bell, C.D.; Mackenzie, D. Critical review of ASME III plasticity correction factors for fatigue design-by-analysis of nuclear power plant components. J. Press. Vessel. Technol. 2021, 143, 061201. [Google Scholar] [CrossRef]
  12. Lang, H.; Rudolph, J.; Ziegler, R. Performance study of Ke factors in simplified elastic plastic fatigue analyses with emphasis on thermal cyclic loading. Int. J. Press. Vessel. Pip. 2011, 88, 330–347. [Google Scholar] [CrossRef]
  13. Slagis, G.C. Meaning of Ke in design-by-analysis fatigue evaluation. J. Press. Vessel. Technol. 2006, 128, 8–16. [Google Scholar] [CrossRef]
  14. ASME. ASME Boiler and Pressure Vessel Code, Case N-779: Alternative Rules for Simplified Elastic Plastic Analysis; ASME: New York, NY, USA, 2010. [Google Scholar]
  15. Shao, X.; Du, J.; Zhang, L.; Xie, H.; Tian, J.; Shi, K.; Zhang, Y.; Liu, Z. Determination method of elasto-plastic correction factor of TA16 titanium alloy based on RCC-M code methodology. J. Press. Vessel. Technol. 2021, 143, 88–93. [Google Scholar] [CrossRef]
  16. Du, J.; Shao, X.; Zhang, L.; Kan, Q. Research on Elastoplastic Strain Correction Factor of Titanium Alloy. Nucl. Power Eng. 2014, 35, 101–105. [Google Scholar]
  17. Kan, Q.H.; Guo, S.J.; Li, J.; Kang, G.Z.; Yan, W.Y. Numerical investigation on plastic strain correction factor in simplified elasto-plastic fatigue analysis. Appl. Mech. Mater. 2017, 853, 226–230. [Google Scholar] [CrossRef]
  18. RCC-MRx. Design and Construction Rules for Mechanical Components of Nuclear Installations: High Temperature, Research and Fusion Reactors; AFCEN: Paris, France, 2018. [Google Scholar]
  19. He, K.; Zhao, Y.; Dang, Y.; Wang, L.; Jiang, F. Fatigue Behavior of TA16 Tube with Surface Defect under High Cycle Loading. Nucl. Power Eng. 2016, 37, 48–51. [Google Scholar]
  20. Luo, Q.; Wang, L.; Chen, X. Studies on the corrosion behavior of TA16 and TA17 titanium alloys in high temperature and high pressure water. Qing Jinshu 2012, 2, 56–59. [Google Scholar]
  21. Shao, X.; Xie, H.; Xiong, F.; Shi, K.; Yu, X. Comparative study on the nonlinear calculation of ratcheting deformation using different constitutive model. In Proceedings of the ASME Pressure Vessels and Piping Conference, Virtual, Online, 3 August 2020. PVP2020-21182. [Google Scholar]
Figure 1. Notched sheet and nozzle finite element model.
Figure 1. Notched sheet and nozzle finite element model.
Metals 12 00943 g001
Figure 2. Finite element simulated von Mises strain contour of notched sheet (X-direction, T = 20 °C, b = 20 mm, r = 25 mm, ratio = −1). (a) Elastic analysis; (b) elastoplastic analysis.
Figure 2. Finite element simulated von Mises strain contour of notched sheet (X-direction, T = 20 °C, b = 20 mm, r = 25 mm, ratio = −1). (a) Elastic analysis; (b) elastoplastic analysis.
Metals 12 00943 g002
Figure 3. Finite element simulated von Mises strain contour of nozzle under axial load. (a) Elastic analysis; (b) elastoplastic analysis.
Figure 3. Finite element simulated von Mises strain contour of nozzle under axial load. (a) Elastic analysis; (b) elastoplastic analysis.
Metals 12 00943 g003
Figure 4. Finite element simulated von Mises strain contour of nozzle under mechanical and thermal loads. (a) Elastic analysis; (b) elastoplastic analysis.
Figure 4. Finite element simulated von Mises strain contour of nozzle under mechanical and thermal loads. (a) Elastic analysis; (b) elastoplastic analysis.
Metals 12 00943 g004
Figure 5. Loading waveforms under different ratios. (a) Ratio = 0; (b) Ratio = −0.5; (c) Ratio = −0.8; (d) Ratio = −1.
Figure 5. Loading waveforms under different ratios. (a) Ratio = 0; (b) Ratio = −0.5; (c) Ratio = −0.8; (d) Ratio = −1.
Metals 12 00943 g005
Figure 6. Stress-strain curves of 409 nodes in notched sheet simulated by Chaboche model under different ratios. (a) Ratio = −1; (b) ratio = −0.8; (c) ratio = −0.5.
Figure 6. Stress-strain curves of 409 nodes in notched sheet simulated by Chaboche model under different ratios. (a) Ratio = −1; (b) ratio = −0.8; (c) ratio = −0.5.
Metals 12 00943 g006
Figure 8. Ke with Se under different constitutive models with notched sheet model.
Figure 8. Ke with Se under different constitutive models with notched sheet model.
Metals 12 00943 g008
Figure 9. Ke with Se under different temperatures with notched sheet model.
Figure 9. Ke with Se under different temperatures with notched sheet model.
Metals 12 00943 g009
Figure 11. Uniaxial tensile stress-strain curves at 30 °C and 350 °C of TA16 titanium alloy.
Figure 11. Uniaxial tensile stress-strain curves at 30 °C and 350 °C of TA16 titanium alloy.
Metals 12 00943 g011
Figure 12. Experimental results under strain control: (a) cyclic stress-strain curve at room temperature; (b) stable cyclic curve and uniaxial tensile curve at room temperature; (c) cyclic stress-strain curves at 350 °C; (d) stable cyclic curve and uniaxial tensile curve at 350 °C.
Figure 12. Experimental results under strain control: (a) cyclic stress-strain curve at room temperature; (b) stable cyclic curve and uniaxial tensile curve at room temperature; (c) cyclic stress-strain curves at 350 °C; (d) stable cyclic curve and uniaxial tensile curve at 350 °C.
Metals 12 00943 g012
Figure 13. Cyclic stress-strain curve under stress control at room temperature (a) and 350 °C (b).
Figure 13. Cyclic stress-strain curve under stress control at room temperature (a) and 350 °C (b).
Metals 12 00943 g013
Figure 14. The minimum conservative margin of TA16 and Z2CND18.12 (nitrogen control).
Figure 14. The minimum conservative margin of TA16 and Z2CND18.12 (nitrogen control).
Metals 12 00943 g014
Table 1. Stress and strain relationship of Z2CND18.12 (nitrogen control) at different temperatures.
Table 1. Stress and strain relationship of Z2CND18.12 (nitrogen control) at different temperatures.
T = 20 °CT = 100 °CT = 200 °CT = 250 °CT = 300 °CT = 350 °C
εσ
(MPa)
εσ
(MPa)
εσ
(MPa)
εσ
(MPa)
εσ
(MPa)
εσ
(MPa)
0.00112200.00091750.00081500.00081380.00081300.0008127
0.00232500.00161900.00191700.00191600.00131400.0012133
0.00272600.00192000.00261800.00271700.00181500.0016140
0.00332700.00242100.00371900.00391800.00261600.0024150
0.00412800.00312200.00522000.00571900.00391700.0036160
0.00512900.00412300.00752100.00852000.00581800.0055170
0.00643000.00742500.01072200.01032050.00881900.0085180
0.00923150.01002600.01292250.01252100.01332000.0131190
Table 2. Load type and peak load range of each model.
Table 2. Load type and peak load range of each model.
ModelMechanical LoadMixed Load
Temperature (°C)Axial Pressure (MPa)Axial Displacement (mm)Temperature (°C)Internal Pressure (MPa)
Model-I3040~1800.2~10--
35040~1200.1~7
Model-II30215~4609.5~2830~450−20~31
350140~2806.5~17
Table 3. Results of Ke under thermal load of Z2CND18.12 (nitrogen control) for nozzle model.
Table 3. Results of Ke under thermal load of Z2CND18.12 (nitrogen control) for nozzle model.
NodeTemp_Load
(°C)
Sn_th
(MPa)
Δεee_th
(%)
Δεep_th
(%)
Ke_th
(EP)
Ke_th
(RCCM)
η (%)
1576350/20693.236.276 × 10−39.762 × 10−31.5551.62811.580
300/20584.815.219 × 10−37.835 × 10−31.5011.58414.211
250/20476.684.190 × 10−36.030 × 10−31.4391.52816.819
200/20398.683.194 × 10−34.373 × 10−31.3691.46620.901
125/20212.001.800 × 10−32.074 × 10−31.1521.26041.625
298350/20748.986.368 × 10−39.993 × 10−31.5691.64211.331
300/20631.305.292 × 10−38.071 × 10−31.5251.61013.981
250/20514.534.251 × 10−36.235 × 10−31.4671.56817.905
200/20369.463.239 × 10−34.486 × 10−31.3851.49221.716
125/20228.731.825 × 10−32.146 × 10−31.1761.36551.851
Table 5. Material parameters of Chaboche model at different temperatures of TA16.
Table 5. Material parameters of Chaboche model at different temperatures of TA16.
Temperature c 1   ( MPa ) γ 1 c 2   ( MPa ) γ 2
30 °C13,0001020,00050
350 °C20,0005100050
Table 6. Simulated results by simplified elastoplastic analysis.
Table 6. Simulated results by simplified elastoplastic analysis.
NodeΔP (MPa)ΔT (°C)Sn (MPa)Sp (MPa)Sm (MPa)KeUmax
N10/16.684/279329.1428.2234.92.200.66
N20/16.684/279350.3349.8234.92.470.57
Table 7. Simulated results of cyclic elastoplastic fatigue analysis.
Table 7. Simulated results of cyclic elastoplastic fatigue analysis.
NodePrincipal Strain RangeSa (MPa)NnU
N10.00516280.34001000.25
N20.00434235.9555.61000.18
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shao, X.; Du, J.; Fu, X.; Xiong, F.; Li, H.; Tian, J.; Lu, X.; Xie, H. Simplified Elastoplastic Fatigue Correction Factor Analysis Approach Based on Minimum Conservative Margin. Metals 2022, 12, 943. https://doi.org/10.3390/met12060943

AMA Style

Shao X, Du J, Fu X, Xiong F, Li H, Tian J, Lu X, Xie H. Simplified Elastoplastic Fatigue Correction Factor Analysis Approach Based on Minimum Conservative Margin. Metals. 2022; 12(6):943. https://doi.org/10.3390/met12060943

Chicago/Turabian Style

Shao, Xuejiao, Juan Du, Xiaolong Fu, Furui Xiong, Hui Li, Jun Tian, Xifeng Lu, and Hai Xie. 2022. "Simplified Elastoplastic Fatigue Correction Factor Analysis Approach Based on Minimum Conservative Margin" Metals 12, no. 6: 943. https://doi.org/10.3390/met12060943

APA Style

Shao, X., Du, J., Fu, X., Xiong, F., Li, H., Tian, J., Lu, X., & Xie, H. (2022). Simplified Elastoplastic Fatigue Correction Factor Analysis Approach Based on Minimum Conservative Margin. Metals, 12(6), 943. https://doi.org/10.3390/met12060943

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