Next Article in Journal
Gemological Characteristic Difference between Colorless CVD Synthetic Diamonds and Natural Diamonds
Next Article in Special Issue
Study on Microstructure and Fatigue Properties of FGH96 Nickel-Based Superalloy
Previous Article in Journal
Dialkyl Succinates and Adipates as Alternative Plasticizers—Even More Efficient Synthesis
Previous Article in Special Issue
Temperature Dependence of Fracture Characteristics of Variously Heat-Treated Grades of Ultra-High-Strength Steel: Experimental and Modelling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Incorporation of Temperature and Plastic Strain Effects into Local Approach to Fracture

1
G. V. Kurdyumov Institute for Metal Physics, National Academy of Sciences of Ukraine 36, Academician Vernadsky Blvd., UA-0380 Kyiv, Ukraine
2
E. O. Paton Institute for Materials Science and Welding, National Technical University of Ukraine “Igor Sikorsky Kyiv Polytechnic Institute” 37, Peremohy Ave., UA-03056 Kyiv, Ukraine
3
Institute of Physics of Materials, Academy of Sciences of the Czech Republic, Zizkova 22, 61669 Brno, Czech Republic
4
Institute of Materials Science and Engineering, Brno University of Technology, Technicka 2, 61669 Brno, Czech Republic
*
Author to whom correspondence should be addressed.
Submission received: 23 August 2021 / Revised: 2 October 2021 / Accepted: 13 October 2021 / Published: 19 October 2021
(This article belongs to the Special Issue Recent Advances in Mechanisms of Fracture and Fatigue)

Abstract

:
An unjustified simplification of the local quantitative criterion regarding cleavage nucleation is a key problem in the utilisation of the Local Approach to Fracture (LA), particularly to predict the fracture toughness within the ductile-to-brittle transition (DBT) region. The theoretical concept of the effect of both temperature and the plastic strain value on the crack nuclei (CN) generation rate in iron and ferritic steels is presented. It is shown how the plastic strain and temperature affect CN formation rate and, as a consequence, govern the shape of the temperature dependence of fracture toughness KJc and its scatter limits. Within the framework of the microscopic model proposed, dependences of the CN bulk density on the plastic deformation value and temperature are predicted. Convenient approximation dependences for incorporating this effect into the LA are suggested. The experimental data of reactor pressure vessel steel and cast manganese steel demonstrate that the use of these dependences enables one to predict, with sufficient accuracy, the effect of temperature on the value of fracture toughness and its scatter limits over the DBT region. It is shown that accounting for both the temperature and strain dependence of CN bulk density gives rise to the invariance of parameters of the Weibull distribution to temperature.

1. Introduction

The Local Approach (LA) to fracture was introduced in the 1980s as a probabilistic method accounting for microstructural effects in the evaluation of global parameters. It aimed to solve key problems of fracture mechanics such as the prediction of the specimen geometry and statistical size effect on fracture toughness (transferability problem), as well as the effects of test temperature, operation temperature, and degradation of the material due to neutron irradiation, etc. The possibility of solving these complex problems was seen in the employment of statistical local criterion for initiating fracture in the vicinity of a crack, because this approach accounts for the key features of the mechanism initiating fracture at the micro-scale.
The Local Approach is based on the postulation that the probability of cleavage initiation within the vicinity of a crack/notch is described by the Weibull distribution. Accordingly, the parameters of this distribution are material constants that do not depend on its stress–strain state and temperature [1]. This approach describes the temperature dependence and the scatter limits of the fracture toughness K I C . In this case, the temperature dependence of the value K I C is governed by the temperature dependence of the magnitude of yield strength. However, in the region of the ductile-to-brittle-transition, where there is a rapid increase in the fracture toughness against the background of a slight decrease in the yield strength, especially for high-strength steels, this approach underestimates the slope of the temperature dependence of the fracture toughness curve. In order to compensate for this shortcoming, one has to agree that shape parameter m and scale factor σ u are not constants; their values depend on the temperature [2,3,4,5,6,7,8]. Moreover, the values of these parameters depend on specimen geometry [2,9,10] and local plastic strain value [6]. All this contradicts the initial statements of the Local Approach. A number of studies have attempted to solve this problem [11,12,13,14,15,16]. These works showed the need to address two key issues, namely: (i) the need to incorporate the effect of temperature and the magnitude of local plastic strain on the CN bulk density ρ; (ii) the need to account for the value of the threshold stress σ th . As for the latter, in fact, this is a methodical problem. It consists of the development of a technique for the experimental determination of σ th . Therefore, a simplified method for determining σ th for structural steels was proposed [17]. According to this technique, the value of σ th is determined based on the results of uniaxial tensile tests of 5–6 notched cylindrical specimens with special geometries (maximum diameter is 5 mm; minimum notched diameter is 5.2 mm; notch radius is 2 mm) at the boiling point of liquid nitrogen (−196 °C).
Ruggieri and Jivkov et al. [15,16] tried to incorporate the effect of plastic strain and test temperature on the number of CN forming in the local plastic zone ahead of a macrocrack tip. This was achieved by revising the expression for the Weibull stress. The nature of this revision was to consider the effect of both temperature and equivalent plastic strain on the CN (“cleavage initiators”) number. This has been accomplished by introducing an empirical “thinning function” that makes it possible to predict “how many cleavage initiators must be generated relative to a reference case (for example all particles)” depending on temperature and the plastic strain value. The values of constants in this function were determined by a calibration procedure using experimental data. This enables one to predict, with high accuracy, the critical values of the J I C integral, both for low and high temperatures within the DBT region.
These studies have convincingly exhibited that the employment of LA within the DBT region requires the consideration of the temperature and deformation dependence of the CN bulk density. Here, in this contribution, the authors emphasised the need to explain the reasons for the significant effect of temperature on the number of formed “cleavage initiators”. As attempts to answer this question, one considers the research results presented in [17,18,19,20], where a physical version of the LA to cleavage fracture was proposed and evidenced. This version of LA does not employ any empirical dependences for the probability of both CN formation and instability. These functions are derived from a detailed analysis of the processes of CN formation and unstable equilibrium in a polycrystalline aggregate. This approach enabled us to ascertain the influence of both the metal structure and the conditions of its loading on the probability of fracture. However, it turned out to be quite sophisticated for engineering calculations. At the same time, this approach can be used as a tool to analyse the key effects that control the initiation of cleavage in the vicinity of the macrocrack, and in particular to analyse the effect of temperature and the magnitude of plastic deformation on the CN formation rate inside of the “process zone”.
The paper presents ideas about the physical nature of the effect of temperature and plastic strain value on the intensity of CN generation during the plastic deformation preceding cleavage. Physically substantiated dependences of CN density on temperature and the plastic strain value are obtained. To account for this effect within the LA, approximate dependences of CN density on temperature and the value of equivalent plastic strain are proposed. By the example of two structural steels, the suitability of using these dependences in predicting the value of fracture toughness within the DBT region is validated.

2. Theory

The inhomogeneity of microplastic deformation, which gives rise to plastic deformation incompatibility on grain and interphase boundaries, is a general reason for the CN formation in polycrystalline solids. A generalised model of CN formation in a polycrystalline aggregate was proposed [19,20].
The specific feature of this model is that it describes the formation and instability of the CN, accounting for both the main structural features of the polycrystal (inhomogeneous distribution of grain sizes and statistical distribution of their orientations) and for the action of microstresses generated by the incompatibility of microplastic deformations in polycrystal.
This model can be adopted for the prediction of the cleavage fracture in structural steels with ferritic basic microstructures arising from the expression for the bulk density ρ of CN formed at a given value of local plastic strain:
ρ = 2 ρ c b   t C t r g ( t )   d t
where
g ( t ) = 1 2 π exp t 2 2
and ρ c b is the carbide particle density.
In this equation, t is the magnitude of normalised shear microstresses ξ n s acting in slip systems:
t = ξ n s D ξ n s
Here, D ξ n s is the variance of shear microscopic stresses ξ n s . Additionally t C is the critical value of normalised microstresses at which the CN are formed:
t C = 1 k σ M + 1 σ _ ξ C C d max β e _ d
where σ _ and e ¯ are the equivalent macroscopic stresses and strains; d and d max are the average and maximum (with a given probability) ferritic grain sizes; ξ C is the critical stress of CN formation as a result of the carbide particle cleavage; M is the orientation factor (for α F e  M = 0.36); k σ , C, β are the coefficients (for ferritic steels k σ = 0.225; C = 0.0336 N/m; β ≈ 2.57 MPa m0.5 [19]. Furthermore, t r is the critical value of normalised microstresses in which there is relaxation of incompatibilities in intergranular boundaries:
t r = 1 k σ M + 1 σ _ ξ Y m b r d max β e _ d
where ξ Y is the critical stress of the relaxation onset; m b is the orientation factor for relaxation slip systems in intergranular boundaries; r is the distance from the grain boundary to the origin of microstress relaxation.
In Equation (5), the expression β e ¯ / d specifies the value of shear microstresses caused by the interaction of a grain of average orientation M with the surrounding matrix plastically deformed to strain e ¯ . In polycrystals, these stresses arise as a result of the accumulation of microplastic strain incompatibilities at grain boundaries. The magnitude of the microplastic strain incompatibility changes nonmonotonically with an increase in macroplastic strain e ¯ . At small macroplastic strains, it grows and then begins to decrease as a result of the predominance of relaxation processes in the near-boundary regions. A transition from growth to a decrease in these incompatibilities occurs at the critical strain e C . At the macroscale, the value of e C can be estimated based on the dependence of the cleavage fracture stress of smooth (unnotched) specimens, σ f , on the value of strain preceding fracture. When the strain e C is achieved, this stress σ f reaches its minimum value. For common structural steels with basic ferritic microstructures, e C ≈ 0.02 [10]. With an increase in ferrite grain size by annealing, the critical strain e C can increase to ≈ 0.05 [10]. Thus, dependences (4) and (5) are valid for strains not exceeding the critical one e C .
At strains e ¯ , greater than e C , the expressions for t C and t r may be presented as follows:
t C = 1 k σ M + 1 σ _ ξ C C d max β e C d + k e e _ e C 1
where k e is the coefficient, characterising the intensity of relaxation process in the near-boundary regions. An estimate of this coefficient based on the [10] data (Figure 1) gives the value of k e   1.52 MPa m0.5.
Or, as the case may be, t r :
t r = 1 k σ M + 1 σ _ ξ Y m b r d max β e C d + k e e _ e C 1
The theoretical dependences of CN bulk density ρ on the magnitude of plastic strain and temperature for reactor pressure vessel (RPV) RPV steel are shown in Figure 1.
The following values of microscopic parameters were used in their building: r = 1 µm (typical distance from grain boundary to stress relaxation source; this value is commensurate with the substructure parameters); m b = 0.1 (the average value of orientation factor for the relaxation source); d = 10 µm and d max = 30 µm (mean and maximum ferrite grain values); ξ C = 7 GPa ( ξ C 0.1 G where G 70   GPa is the iron carbide shear modulus) [21].
In the first approximation, the mean value of ξ Y may be estimated by the value of the thermally activated component of the yield strength ξ Y τ Y . The expression for τ Y is the following [22]:
τ Y = 0.5 C 1 exp C 2 C 3 ln e ˙ T
where e ˙ is the plastic strain rate (for quasi-static tension e ˙ = 10−4·s−1); C1, C2 and C3 are the constants, which values for typical ferritic steels are: C1 = 1033 MPa, C2 = 0.0068 K−1, C3 = 0.000415 K−1 [22].
The equivalent stress value σ ¯ was calculated according to equation
σ ¯ = σ 0.2 e ¯ 0.002 n
where n is the work hardening exponent (value n = 0.05 was used) and σ 0.2 is the yield strength:
σ 0.2 = σ a + C 1 exp C 2 + C 3 ln e ˙ T
where σ a is the athermal component of σ 0.2 .
For the investigated reactor pressure vessel (RPV) steel at room temperature (T = 293 K), the yield strength was measured to be σ 0.2 = 610   MPa . Accounting for Equation (10) and given the above C1, C2 and C3 values, the athermal component of the yield strength is σ a = 564 MPa. The temperature dependences of both yield strength σ 0.2 and the work hardening exponent n for RPV steel utilised in calculations are presented in Figure 2a.
According to the calculation results (Figure 1), the non-monotonic dependence of ρ on the plastic strain value is a characteristic feature of the CN formation in polycrystalline metals. As was shown above, this is due to non-monotonic change in shear microstresses, induced by microplastic strain incompatibility at the grain boundaries, with macroplastic strain e ¯ growth.
As the calculations show, the temperature increase should cause a monotonic decrease in the value of the CN density ρ. This is due to the thermally activated relaxation processes at the grain boundaries. In terms of the current model, this means a decrease in the value t r (dependences (7) and (5)) against the background of increasing the critical value of normalised stress of CN formation t C (dependences (6) and (4)).
Figure 1 shows that the effect of plastic strain magnitude and temperature on the CN bulk density ρ is a non-additive one. The maximum influence of plastic deformation on ρ is observed at strains close to the critical one e C . With this deformation, the maximum sensitivity of ρ to temperature changes is also observed. At large strains, the sensitivity of ρ to changes in temperature decreases. In general, a comparison of the contribution of plastic strain and temperature to change in the intensity of the crack nuclei formation shows that the prevailing contribution is made by the change in temperature. An increase in temperature from +100 °C to −196 °C can lead to an almost fourfold increase in ρ. A decrease in temperature also increases the sensitivity of ρ to the plastic strain value. The dependences of CN density described above were obtained based on the analysis of the process of CN forming at microscale. Therefore, they enabled us to ascertain the fundamental regularities of the influence of temperature and plastic strain on the value of CN bulk density ρ, as well as determine the physical nature of this effect. However, these dependences cannot be directly used in LA since they contain microscopic parameters, the values of which can be estimated only with an accuracy of an order of magnitude. To incorporate these regularities into the LA, it is important to find the appropriate functions to approximate these dependences ρ on the temperature and the magnitude of plastic strain. The constants in these functions characterise the influence of microstructure parameters of specific steel on the sensitivity of ρ to both temperature and the magnitude of plastic strain. The values of these constants can be estimated from experimental evidence by means of calibration procedures.
In this case, as a first approximation, the following expressions can be used:
ρ = ρ C a 1 e _ e C   for   e ¯ e C
and
ρ = ρ C b e _ e C 1   for   e max e ¯ e C
where
ρ C = ρ 0 1 exp α τ Y
In the above equation, α and ρ 0 are coefficients, the values of which are calculated using experimental evidence by the calibration procedure, e C = 0.02 and e max 0.3 ÷ 0.5 .
The dependences in Figure 1 give the following values: a = 1.50 × 1013 m−3; b 0.12 × 1013 m−3.

3. Experiment

Experimental investigations were performed on a reactor pressure vessel steel and on a low–alloyed manganese steel. To determine the value of fracture toughness, standard (ASTM E1921-05 [23]) specimens 1T CT (for the reactor pressure vessel steel) and 1T SENB (single edge notched bending) (for the low–alloyed manganese steel) were used. Chemical composition of the steels is presented in Table 1. More details on materials used can be found in our previous works, e.g., [17,24,25]. Both steels were investigated as received and had a ferritic structure. Mechanical tests of the pre-cracked specimens were executed in the DBT temperature range (Figure 3). The specimens were tested under crosshead speed of 1 mm/min in either environmental and/or cryogenic chamber. Experiments were carried out on a servo-hydraulic testing machine, Instron 8802. The value of K J c was estimated according to the ASTM E1921-05 standard. According to this standard, the maximum limiting the valid values of elastic plastic fracture toughness, K J c limit , was calculated according to the equation:
K J c limit = E b 0 σ 0.2 30 1 ν 2
where E is the Young modulus, (E = 210 GPa); ν is Poisson’s ratio (ν = 0.27); b 0 is the unbroken ligament ( b 0 = 22.4 mm); σ 0.2 is the yield strength (the temperature dependences for σ 0.2 of the investigated steels are presented in Figure 2). The temperature values of K J c limit are presented in Figure 4. To determine the mechanical properties of the investigated steels and obtain true stress vs. true deformation curves, low-temperature tests for uniaxial tension of smooth (unnotched) specimens were performed. According to the results of these tests, the values of the yield stress and the strain hardening exponent were determined (Figure 2).
According to the results of these tests, the magnitudes of Weibull parameters were determined by calibration (shape parameter m and scale parameter σ u ) as well as the values of ρ 0 and α. To determine local parameters, a program was created that allows the calculation of local parameters using the maximum likelihood method for a three-parameter modification of the Beremin approach. This is a typical procedure for the Local Approach. The difference is only that V 0 is not constant but depends on ρ ( V 0 = 1 / ρ ),the value of which is determined by dependences (11) and (12).
The magnitude of threshold stress σ t h was estimated by the technique described in [17] (Table 2).

4. Results and Discussion

Relatively small m values (m ≤ 10) are the consequence of accounting for the values of threshold in the calculation [26]. During the calibration procedure, it was possible to find such a set of Weibull parameters, at which the values of ρ 0 and α for both steels are almost the same (Table 3). This validates the ascertained values ρ 0 and α as those inherent to ferritic steels. In this case, the features of the microstructure of steels are accounted for by the values of the Weibull parameters. It should be emphasised that accounting for the ρ temperature dependence ρ (in terms of the Weibull distribution—for the parameter V 0 ( ρ = V 0 1 ) gives rise to the independence (invariance) of the Weibull parameters σ u and m on temperature. This allows us to consider σ u and m as material parameters. Contrarily, in the conventional versions of Local Approach, it is assumed that the quantity V 0 is constant. This leads to temperature dependences of parameters σ u and m, which contradicts the physical sense and complicates calculations.
Based on these data, temperature dependences of fracture toughness K J c for the studied steels were plotted (Figure 3a,b). The algorithm of calculations of fracture toughness and fracture probabilities was as follows:
1. For each jth stage of loading K I j , the values of equivalent plastic strain e ¯ i and normal tensile stresses σ Y Y i in each i t h finite element were calculated using the finite elements method (software ABAQUS, version 2018).
2. Further, the probability of the cleavage initiation in the i t h element was calculated as
P i = 1 exp V i V 0 i × σ Y Y i σ t h σ u m
where V i is the finite element volume and V 0 i is the volume per one crack nucleus ( ρ = V 0 1 ). The value of ρ i was calculated according to Equations (11)–(13).
3. Then, the total probability of cleavage initiation in specimen is determined as
P Σ = 1 n = 1 n = N 1 P i n
where N is the number of finite elements belonging to the yield region.
The step size and the number of loading stages were chosen in such a way as to obtain K I values when the probabilities of cleavage initiation P Σ are 5%, 50% and 95%.
In accordance with the calculation results, the theoretical curve describes the temperature dependence of fracture toughness K J c with a 90% probability scatter band (Figure 3). It should be noted that such a coincidence is due to the fact that not only a decrease in yield strength, but also a reduction in the CN generation rate is the reason for K J c growth with increasing temperature. If we do not take into account this fact, i.e., if we assume that ρ = const, as the case may be, V 0 = const ( V 0 = 1/ ρ), then the slope of K J c temperature dependence becomes unbelievably low (Figure 3, dashed lines). A particularly significant error occurs at low values of the failure probability, i.e., when determining K J c threshold values). This is especially important for applications. It should be mentioned that accounting for ρ temperature dependence provides invariance to temperature of the Weibull distribution parameters, which agrees well with the LA concept that considers m and σ u as material parameters.
The proposed approach can also be used to analyse the Master Curve method. For this purpose, the values of reference temperatures T 0 were determined for the studied steels by the multitemperature method, and, according to the ASTM E1921-05 standard, the temperature dependences of fracture toughness were plotted (Figure 4).
According to these data, the Master Curve method overestimates the value K J c at temperatures above the reference one T 0 . This means the need to “correct” the width of the temperature range near T 0 , within which the approximation of experimental evidence by the Master Curve method, is valid. The approach accounting for the temperature dependence of CN density formation may be used as a tool to fix this problem.
Fracture toughness is a stochastic characteristic. This is most pronounced within the DBT region. Therefore, it is necessary to pay attention to errors not only of K J c absolute values, but also to errors of probabilities of their realisation. Figure 4 exhibits that the error in the K J c value due to neglecting of ρ temperature dependence is in the tens of percent. At the same time, the error in probabilities of realisation of the given K J c values can be quite high.

5. Conclusions

  • In a general case, the incompatibility of microplastic deformations at grain or interphase boundaries causes the crack nuclei (CN) formation. The value of this incompatibility depends on the magnitude of plastic strain and temperature. This is the reason for the influence of these factors on the intensity of the crack nuclei generation. The microscopic model presented in the work enables one to describe, quantitatively, the regularities of this effect.
  • The strain and temperature dependence of the CN bulk density results in the invariance of the Weibull distribution parameters both to the temperature and the plastic strain value, i.e., σ u and m, may be considered as material parameters. This is in line with the original concepts underlying the Local Approach.
  • Decreasing CN generation intensity with temperature growth is one of the key factors that governs the slope of K J c temperature dependence and its scatter limits. Ignoring this effect underestimates the slope of K J c temperature dependence and their absolute values. Especially large, up to 20–30%, is the error at low probabilities P Σ = 5% fracture (lower fracture toughness threshold).
  • The Master Curve method, vice versa, overestimates K J c sensitivity to temperature. This may result in the overestimation of the fracture toughness value at temperatures higher than the reference one T 0 .

Author Contributions

Conceptualization, S.K. and I.D.; methodology, S.K.; software, V.K.; validation, O.Z. and G.Z.; formal analysis, N.S.; investigation, O.Z.; data curation, G.Z.; writing—original draft preparation, S.K.; writing—review and editing, I.D.; visualization, N.S.; project administration, S.K.; funding acquisition, I.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Academy of Sciences of Ukraine, grant number 0121U107569; the APC is funded by Institute of Materials Science and Engineering, Brno University of Technology, Brno, Czech Republic.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data supporting the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Beremin, F.M.; Pineau, A.; Mudry, F.; Devaux, J.-C.; D’Escatha, Y.; Ledermann, P. A local criterion for cleavage fracture of a nuclear pressure vessel steel. Met. Trans. A 1983, 14, 2277–2287. [Google Scholar] [CrossRef]
  2. Wiesner, C.S.; Goldthorpe, M.R. The effect of temperature and specimen geometry on the parameters of the ‘Local Approach’ to cleavage fracture. J. Phys. IV 1996, 6, C6–295. [Google Scholar] [CrossRef]
  3. Wasiliuk, B.; Petti, J.R.; Dodds, R.H. Temperature dependence of Weibull stress parameters: Studies using euro-material. Eng. Fract. Mech. 2006, 73, 1046–1059. [Google Scholar] [CrossRef]
  4. Pineau, A. Development of the local approach to fracture over the past 25 years: Theory and applications. Int. J. Fract. 2006, 138, 139–166. [Google Scholar] [CrossRef]
  5. Ruggieri, C.; Dodds, R.H. A local approach to cleavage fracture modelling: An overview of progress and challenge for engineering applications. Eng. Fract. Mech. 2018, 187, 381–403. [Google Scholar] [CrossRef]
  6. Gao, X.; Dodds, R.H., Jr. Constraint effects on the ductile-to-brittle transition temperature of ferritic steels: A Weibull stress model. Int. J. Fract. 2000, 102, 43–69. [Google Scholar] [CrossRef]
  7. Qian, G.; González-Albuixech, V.F.; Niffenegger, M. Calibration of Beremin model with the Master Curve. Eng. Fract. Mech. 2015, 136, 15–25. [Google Scholar] [CrossRef]
  8. Gao, X.; Zhang GSrivatsan, T.S. A probabilistic model for prediction of cleavage fracture in the ductile-to-brittle transition region and the effect of temperature on model parameters. Mat. Sci. Eng. 2006, 415, 264–272. [Google Scholar] [CrossRef]
  9. Cao, Y.; Hui, H.; Wang, G.; Xuan, F.-Z. Inferring the temperature dependence of Beremin cleavage model parameters from the Master Curve. Nucl. Eng. Des. 2011, 241, 39–45. [Google Scholar] [CrossRef]
  10. Kotrechko, S.; Strnadel, B.; Dlouhý, I. Fracture toughness of cast ferritic steel applying local approach. Appl. Fract. Mech. 2007, 47, 171–181. [Google Scholar] [CrossRef]
  11. Bordet, S.R.; Karstensen, A.D.; Knowles, D.M.; Wiesner, C.S. A new statistical local criterion for cleavage fracture in steel. Part I: Model presentation. Eng. Fract. Mech. 2005, 72, 435–452. [Google Scholar] [CrossRef]
  12. Bordet, S.R.; Karstensen, A.D.; Knowles, D.M.; Wiesner, C.S. A new statistical local criterion for cleavage fracture in steel. Part II: Application to an offshore structural steel. Eng. Fract. Mech. 2005, 72, 453–474. [Google Scholar] [CrossRef]
  13. Gao, X.; Zhang, G.; Srivatsan, T.S. Prediction of cleavage fracture in ferritic steel: A modified Weibull stress model. Mat. Sci. Eng. A 2005, 394, 210–219. [Google Scholar] [CrossRef]
  14. Ruggieri, C.; Savioli, R.G.; Dodds, R.H., Jr. An engineering methodology for constraint corrections of elastic–plastic fracture toughness—Part II: Effects of specimen geometry and plastic strain on cleavage fracture predictions. Eng. Fract. Mech. 2015, 146, 185–209. [Google Scholar] [CrossRef]
  15. Ruggieri, C.; Jivkov, A.P. A local approach incorporating the measured statistics of microcracks to assess the temperature dependence of cleavage fracture for a reactor pressure vessel steel. Proc. Struct. Int. 2019, 18, 28–35. [Google Scholar] [CrossRef]
  16. Jivkov, A.P.; Burgos, D.S.; Ruggieri, C.; Beswick, J.; Savioli, R.G.; James, P.; Sherry, A. Use of local approaches to calculate changes in cleavage fracture toughness due to pre-straining and constraint effects. Appl. Fract. Mech. 2019, 104, 102380. [Google Scholar] [CrossRef]
  17. Kotrechko, S.; Zatsarna, O.; Kozák, V.; Dlouhý, I. Threshold fracture stress: Theory and application. Proc. Struct. Int. 2019, 23, 413–418. [Google Scholar] [CrossRef]
  18. Kotrechko, S. Physical Fundamentals of Local Approach to Analysis of Cleavage Fracture. In Transferability of Fracture Mechanical Characteristics; Dlouhý, I., Ed.; NATO Science Series (Series II: Mathematics, Physics and Chemistry); Springer: Dordrecht, The Netherlands, 2002; Volume 78, pp. 135–150. [Google Scholar] [CrossRef]
  19. Kotrechko, S. The key problems of local approach to cleavage fracture. J. Theor. Appl. Mech. 2013, 51, 75–89. [Google Scholar]
  20. Kotrechko, S.; Mamedov, S. Multi-scale local approach to cleavage fracture and its applications. In Proceedings of the 19th European Conference on Fracture (ECF19), Kazan, Russia, 26–31 August 2012; Curran Associates, Inc.: New York, NY, USA, 2016; pp. 971–982. [Google Scholar]
  21. Mizubayashi, H.; Li, S.J.; Yumoto, H.; Shimotomai, M. Young’s modulus of single phase cementite. Scr. Mat. 1999, 40, 773–777. [Google Scholar] [CrossRef]
  22. Zerilli, F.J.; Armstrong, R.W. Dislocation-mechanics-based constitutive relations for material dynamics calculations. J. Appl. Phys. 1987, 61, 1816–1825. [Google Scholar] [CrossRef] [Green Version]
  23. Standard Test Method for Determination of Reference Temperature, Т0, for Ferritic Steels in the Transition Range. Annual book of ASTM Standards 2005; Section 3: Metals Test Methods and Analytical Procedures; v.03.01: Metals–Mechanical Testing, Elevated and Low-Temperature Tests, Metallography. Available online: https://www.astm.org/Standards/E1921 (accessed on 18 October 2021).
  24. Holzmann, M.; Jurášek, L.; Dlouhý, I. Master Curve Methodology and Data Transfer from Small on Standard Specimens. In Transferability of Fracture Mechanical Characteristics; Dlouhý, I., Ed.; NATO Science Series (Series II: Mathematics, Physics and Chemistry); Springer: Dordrecht, The Netherlands, 2002; Volume 78, pp. 225–242. [Google Scholar] [CrossRef]
  25. Dlouhy, I.; Chlup, Z.; Kozak, V. Constraint effects at brittle fracture initiation in a cast ferritic steel. Eng. Fract. Mech. 2004, 71, 873–883. [Google Scholar] [CrossRef]
  26. Ruggieri, C. Influence of threshold parameters on cleavage fracture predictions using the Weibull stress model. Int. J. Fract. 2001, 110, 281–304. [Google Scholar] [CrossRef]
Figure 1. (Colour online). The effect of temperature T and plastic strain value e ¯ on the CN density ρ for the RPV steel: e ¯ is the equivalent plastic strain; e C is the critical value of strain ( e C = 0.02).
Figure 1. (Colour online). The effect of temperature T and plastic strain value e ¯ on the CN density ρ for the RPV steel: e ¯ is the equivalent plastic strain; e C is the critical value of strain ( e C = 0.02).
Materials 14 06224 g001
Figure 2. The temperature dependences of yield strength σ 0.2 and work hardening exponent n for RPV steel (a) and cast manganese steel (b).
Figure 2. The temperature dependences of yield strength σ 0.2 and work hardening exponent n for RPV steel (a) and cast manganese steel (b).
Materials 14 06224 g002
Figure 3. (Colour online). The temperature dependences of fracture toughness K J c for RPV steel (a) and cast low-alloyed manganese steel (b): circles represent the experimental values; lines are the calculation results (solid black lines are the results of K J c calculation accounting for temperature dependence for the CN density ρ, dashed blue lines are the same at the constant values of ρ).
Figure 3. (Colour online). The temperature dependences of fracture toughness K J c for RPV steel (a) and cast low-alloyed manganese steel (b): circles represent the experimental values; lines are the calculation results (solid black lines are the results of K J c calculation accounting for temperature dependence for the CN density ρ, dashed blue lines are the same at the constant values of ρ).
Materials 14 06224 g003
Figure 4. (Colour online). The temperature dependences of fracture toughness K J c for RPV steel (a) and cast low-alloyed manganese steel (b): circles show the experimental evidence, lines are the calculation results (solid black lines are the results of K J c calculation accounting for the CN density ρ and dashed blue lines are the results of K J c calculation by the Master Curve method); T 0 is the reference temperature (short-dashed red lines are the temperature dependences of elastic plastic fracture toughness validity limit K J c limit ), calculated according to Equation (14).
Figure 4. (Colour online). The temperature dependences of fracture toughness K J c for RPV steel (a) and cast low-alloyed manganese steel (b): circles show the experimental evidence, lines are the calculation results (solid black lines are the results of K J c calculation accounting for the CN density ρ and dashed blue lines are the results of K J c calculation by the Master Curve method); T 0 is the reference temperature (short-dashed red lines are the temperature dependences of elastic plastic fracture toughness validity limit K J c limit ), calculated according to Equation (14).
Materials 14 06224 g004
Table 1. Chemical composition of the investigated steels (wt %, remainder is Fe).
Table 1. Chemical composition of the investigated steels (wt %, remainder is Fe).
CMnSiCrNiMoCuSPVAsCoSnSb
RPV Steel
0.170.460.272.181.190.540.13<0.010.010.10<0.003<0.030.006<0.003
Cast Ferritic Steel
0.091.180.370.120.290.030.290.0250.01N/AN/AN/AN/AN/A
Table 2. Mechanical properties of the tested steels ( σ 0.2 is the yield strength determined as 0.2 pct proof stress; σ u l is the ultimate tensile strength; e u is the uniform elongation; ψ is the reduction in area—all at room temperature; n is the work hardening exponent).
Table 2. Mechanical properties of the tested steels ( σ 0.2 is the yield strength determined as 0.2 pct proof stress; σ u l is the ultimate tensile strength; e u is the uniform elongation; ψ is the reduction in area—all at room temperature; n is the work hardening exponent).
Material σ 0.2 ( MPa ) σ u l ( MPa ) e u (---) ψ (---)n (---)
RPV steel6107140.070.750.06
Cast steel3194810.260.740.15
Table 3. Calculation findings ( σ t h is the threshold stress; m and σ u are the shape and scale parameters; ρ 0 and α are the coefficients; T 0 is the reference temperature).
Table 3. Calculation findings ( σ t h is the threshold stress; m and σ u are the shape and scale parameters; ρ 0 and α are the coefficients; T 0 is the reference temperature).
Material σ t h ( MPa ) σ u ( MPa ) m ρ 0 ×   10 13   ( m 3 ) α (MPa−1) T 0 ( ° C )
RPV steel110068355.440.024−70
Cast steel72037008.040.020−145
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kotrechko, S.; Kozák, V.; Zatsarna, O.; Zimina, G.; Stetsenko, N.; Dlouhý, I. Incorporation of Temperature and Plastic Strain Effects into Local Approach to Fracture. Materials 2021, 14, 6224. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14206224

AMA Style

Kotrechko S, Kozák V, Zatsarna O, Zimina G, Stetsenko N, Dlouhý I. Incorporation of Temperature and Plastic Strain Effects into Local Approach to Fracture. Materials. 2021; 14(20):6224. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14206224

Chicago/Turabian Style

Kotrechko, Sergiy, Vladislav Kozák, Oleksandra Zatsarna, Galyna Zimina, Nataliya Stetsenko, and Ivo Dlouhý. 2021. "Incorporation of Temperature and Plastic Strain Effects into Local Approach to Fracture" Materials 14, no. 20: 6224. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14206224

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