Calculation and comparison on fracture toughness of specific reliability between ASTM and ISO standards

Fracture toughness, as a crucial mechanical property, is essential in the research of fatigue crack propagation and fatigue life. Considering the enormous data scatter of the fracture toughness in experiments, a data processing method is proposed to calculate the fracture toughness for the specified reliability. Based on the traditional data processing method shown in the relative standards, the two-dimensional, one-sided tolerance factor is introduced to calculate the influence of the reliability. In addition, by comparing the data processing methods of the fracture toughness presented in the ISO12135:2016 and ASTM E1820-2013 standards, the differences between them were mainly the range of the qualified data, the function of the construction line, and the equation of the regression line. Because the construction lines of the two standards are based on two different material constitutive models, if the constitutive model of material is more aligned with that of the ideal elastic plastic material, the ASTM E1820-2013 standard is recommended; otherwise, the ISO12135:2016 standard is recommended. For 30CrNi2MoVA, the ISO 12135:2016 standard was more suitable, and the error of the JIC obtained by the two standards was about 5%.


Introduction
The fracture toughness, as one of the most important mechanical properties, indicates a material's resistance to crack propagation. It is used widely in the structural design and fatigue life predictions of structures with defects [1,2]. An accurate and reliable determination of the fracture toughness is of great importance for the failure assessment of equipment.
At present, the fracture toughness tests are performed mainly following the ISO 12135:2016 and ASTM E1820-2013 standards, which are the most accurate methods for measuring the fracture toughness [3,4]. For different materials and working conditions, different kinds of specimens and experimental procedures have been proposed in the literature to measure the fracture toughness. In 1981, the small punch (SP) technique was proposed for radiation embrittlement studies in the United States. This technique was extensively developed by many researchers to estimate the potential for residual life assessment of components during service [5][6][7][8]. The Charpy impact test was proposed to avoid the disadvantages of traditional CT(Compact Tension) and SENB (Single Edge Notch Beam) samples. In contrast, this test requires smaller specimens and simpler experimental equipment. However, its result only shows the susceptibility of steel to brittle fractures instead of providing the value of fracture toughness [9]. Chen et al. proposed a method to accurately estimate the fracture toughness for 16MnDR steel in the transition region under different test conditions using the master curve method [1,10]. The cracked Brazilian disk (BD) method was proposed by Atkinson to calculate the fracture toughness data for brittle materials [11]. The spherical indentation test (SIT) was proposed to measure the fracture toughness for ductile metals [2]. Jelitto et al proposed an experimental method to measure the fracture toughness values of porous materials. However, the CT and three-point bending specimens are the two most common specimens [12]. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
In particular, CT samples require smaller dimensions and are more suitable for tests of high fracture toughness values.
However, due to defects in the material itself and randomness during the sample processing and experimental loading, even though the experiment process is strictly followed and controlled according to the experiment standard, the test results of the fracture toughness inevitably exhibit considerable statistical dispersion [13]. With the development of a reliable design, the influence of the uncertainties has attracted growing amounts of attention of researchers [14]. The scattered behavior of the fracture toughness of fiberreinforced hybrid bio-composites has been investigated in detail by Chandramohan et al. [15,16]. Tagawa et al. researched the in the ductile-brittle transition of a structural steel, and the Weibull distribution was applied to describe the scatter characteristics of the fracture toughness [17,18]. Lia and Zhou proposed a modified analytical model of the fracture toughness of composite materials, which allowed the possible range of fracture toughness values to be predicted as a function of the microstructure [19]. Wallin proposed a new statistical assessment method to analyze the Euro fracture toughness data set, and the results could be used to form a new basis for the micro-mechanistic modelling of cleavage fracture [20]. Xu et al proposed a method to calculate the scatter of the fracture toughness using the normal function. However, this method was not able to account for historical data, and a vast amount of data were required [21].
Although many studies have examined the reliability and scatter of fracture toughness, the studies on how to obtain the fracture toughness with a specified reliability are rare. However, in the design of the structural damage tolerance considering the reliability, this parameter is necessary. Generally, a high reliability of the fracture toughness cannot be obtained without a large number of samples and tests. Herein, a method to obtain the value of the fracture toughness with a specific reliability is proposed by introducing a two-dimensional one-sided tolerance factor.
Tests of the fracture toughness can be divided into two types: tests on linearly elastic materials and tests on elastoplastic materials. The former is usually used for brittle materials, and for highly ductile metal materials, such as those used in pressure vessels, the latter test is more suitable. The fracture toughness is usually characterized by the J-integral. Because the material tested in this study is used for high-pressure vessels, the elastoplastic fracture toughness J IC is more reasonable.
The outline of this article is as follows: In the 'Materials and Methods' section, the basic material properties and test procedure for the fracture toughness are presented, and the J-Δa data set of four specimens was obtained. Data processing was conducted following the ISO 12135:2016 standard, and the reliability was considered. A method to calculate the fracture toughness J IC at the specified reliability is proposed. Furthermore, a comparison between the J IC values obtained by the ASTM E1820-18ae1 and ISO 12135:2016 methods was performed and discussed. Finally, conclusions are presented.

Material properties
The material tested in this study was 30CrNi2MoVA, a high-strength steel alloy. The main components, heat treatment process, and basic mechanical property are presented in the table 1-3, respectively [22].

. Specimen preparation
For stable crack growth, the J-Δa curve is generally used to describe the ability of the material to resist the crack extension, and it is also called the J−R resistance curve or J-R curve. The critical J-integral, J IC , at the onset of cleavage fracture is considered to be the fracture toughness, and it is a point on the J−R curve. The J-Δa curve is the foundation of the test of the fracture toughness, J IC . In this study, J IC was obtained by a single sample method. Namely, one whole J-Δa curve can be obtained by one specimen. To improve the reliability of the test results, four sets of J-Δa curves were obtained using four specimens. As shown in figure 1, a CT specimen was used in the test, and the processed specimen is shown in figure 2. Their specific dimensions are shown in table 4.

Test procedure
The test was conducted strictly following the standard, and the test procedure can be divided into four stages: (1) fatigue pre-cracking procedure, (2) side grooves, (3) loading and testing, and (4) data recording and analysis.
The main machine and gauge involved in this test are shown in figure 3.   (1) The fatigue pre-cracking was conducted under a controlled force. The load type was a sinusoidal cyclic load with a stress amplitude of 10.00 kN and stress ratio of 0.1. As the crack grew, the stress amplitude gradually decreased. However, over the entire loading procedure, the loading amplitude was less than 1 kN.
(2) In the test of J IC , side grooves are highly recommended when the compliance method of the crack length prediction is used. The side grooves also ensure that a straight crack front can form. The specimens with side grooves are shown in figure 4. They are marked as J-01, J-02, J-03, and J-04.
In the test of the J-Δa curve, the specimen should be broken by tension using the fatigue machine after finishing the loading. A fracture surface examination was conducted to determine the initial crack length a 0 and the stable crack growth length Δa during the testing process. The macro-fracture of the broken specimen is shown in figure 5.
The macrofracture of the broken specimen can be divided into four zones: (a) Fatigue crack growing zone. In this zone, the fracture surface was smooth and flat. It was mainly produced in the fatigue pre-cracking stage. (b) Stretch zone. The front of the zone was arch shaped, and it was mainly formed during the crack initiation stage. (3) Stable crack zone. This zone was produced during the stable crack propagation stage under the effect of stable cyclic loading. (4) Final fracture zone. This zone was formed in the last tension broken stage of the experiment, and the surface of this zone was very coarse.

Test results
The Δa and J data of the four specimens were calculated and recorded by the gauges automatically. J-Δa curves of the specimens are shown in the figure 6 and are distinguished using different marker shapes.
As mentioned above, the ISO 12135:2016 and ASTM E1820-18ae1 standards are the two most common methods for testing the material fracture toughness. They involve similar experimental processes to obtain the J-Δa data but have a few differences in the data processing procedure. Thus, in the next section, the data processing procedure of the fracture toughness is performed using the two standards.   3. Data processing considering the data reliability (ISO 12135:2016) In the ISO 12135:2016 standard, a data processing method is provided. The data for a crack extension with Δa>0.1 is used to calculate the J IC and is fit by the following equation: where α, C 1 0, and 0C 2 1.
During the data fitting process, it was found that the value of α was close to zero, so its value was set to α=0 for fitting convenience. Finally, the fitting results of the four specimens are shown in table 5. Their fitting curves are shown in figure 6 and expressed by the different line types. Figure 6 shows that the J-Δa curves of the four specimens had the same variation trend. However, for the different specimens, there were subtle differences between the J-Δa curves. This mainly resulted from the differences in the material microstructure, processing   technology, and heat treatment process. Generally, the average fracture toughness values of the different specimens are determined to be the final measured values. However, with the wide application of the reliability design method, it would be better if the highly reliable fatigue fracture parameters could be used in the damage tolerance design and the fatigue life prediction. Thus, in this study, the two-dimensional, one-sided tolerance factor was used to obtain a J-Δa curve considering the reliability. It is marked as the p-γ-J-Δa curve. Statistical data has shown that the dispersion of the fracture toughness is greater than that of the other mechanical property parameters. The fracture toughness approximately follows the three-parameter Weibull distribution, lognormal distribution, and normal distribution [13]. The Weibull distribution is usually applied in studies of the relationship between the fracture toughness and temperature. With fewer specimens, it is difficult to accurately estimate the three parameters of the Weibull distribution. According to the research of Liu and Xiong, the lg(J)-(Δa) data follows a normal distribution [23,24], which can be expressed as follows: m s lg J N , As equation (3) shows, y is linearly related to x. As lg(J) follows a normal distribution, when the reliability is p, lg(J) can be expressed as follows: where y p is the value of lg(J) when the reliability is p, and u p is the p quantile of the standard normal distribution, i.e., The value of u p can be determined using the normal distribution table. However, the true value of [lg(J)] p cannot be obtained in practice. Hence, the upper and lower confidence limits of [lg(J)] p are usually applied in the reliability analysis. In this study, a conservative value of the fracture toughness was sought, so the lower confidence limit of [lg(J)] p was used, denoted as [lg(J)] pl .
The sample average J lg( ) and sample variance S 2 are the unbiased estimators of population mean μ 0 and population variance σ 2 , respectively. where k is the two-dimensional, one-sided tolerance factor. It is calculated as follows: where uγ is a quantile of the standard normal distribution, n is the degree of the sample mean J lg , ( ) and f is the degree of sample variance S 2 . J lg( ) and S 2 may come from different samples, so there is no certain relationship between f and n. When J lg( ) and S 2 come from the same sample, f=n − 1. In this case, k is transformed into a one-dimensional, one-sided tolerance factor. If k is determined, then [lg(J)] pl can be calculated using equation (6). The lower confidence limit of the J-Δa curve with the given reliability p and confidence coefficient γ can be obtained.
In this study, four sets of experimental data were obtained, namely n=4 and f=3. Substituting the given reliability p and confidence coefficient γ into equation (8), the value of k can be obtained. According to equation (6), the key for the calculation of [lg(J)] pl lies in the calculation of J lg( ) and S 2 . For the four specimens, the crack extension Δa recorded during the experiment process was not always the same. To obtain the value of lg(J) corresponding to the different crack extensions Δa of the different specimens, the change range of Δa was set to 0.1-1.6 mm. For a given crack extension Δa, the corresponding value lg(J) of the specific specimen is marked as lg(J) ji . It is calculated as follows: where, j is the number of the specimens (e.g., 1,2,3, or 4), and i is the data point corresponding to the crack extension Δa. For every data point, the relative sample average J lg i ( ) and sample variance S i 2 can be calculated as follows: Here, p=95%, γ=95% and p=99%, γ=95% were considered. In addition, the average J-Δa curve of the four specimens was also calculated for comparison. The curves are shown in figure 7, and their corresponding J-Δa equations are presented in table 6.
As shown in figure 7, the solid line with the rectangle markers represents the average J-Δa curve of the four specimens, which is marked as (J-Δa) ave . The solid line with the circle markers represents the J-Δa curve when the reliability p=95%, which is marked as (J-Δa) p=95% . The (J-Δa) p=99% corresponding to a reliability of p=99% is shown by the solid line with triangle markers. When the crack extension Δa was the same, (J) ave > (J) p=95% >(J) p=99% . This means the if the data reliability is taken into consideration, the fracture toughness J IC obtained from the J resistance curve is smaller, and a more conservative result is obtained.
Next, according to the ISO 12135:2016 standard, J versus Δa is plotted in figure 8 and expressed by the solid line with solid symbols. For this standard, the construction line was determined based on that the material's true stress-strain relationship satisfies a power law [25]. The construction can be expressed as follows: where R m is the ultimate tensile strength in MPa. As shown in figure 8, the black straight line is the construction line. The data located in the zone marked by red is the qualified data. The left and right boundaries of the qualified data region are shown by lines parallel to the construction line intersecting the abscissa at 0.15 and 1.5 mm, respectively. They are marked as l 0.15(ISO) and l 0.5(ISO) , respectively. A line parallel to the construction and exclusion lines is plotted at an offset value of 0.2 mm. It was marked as l 0.2(ISO) , and its equation can be expressed as follows: The intersection of the line l 0.2(ISO) and the J-Δa curve is defined as J 0.2BL , which can be calculated by solving equations (6) and (14). However, the validity must be judged in accordance with the ISO 12135:2016 7.6.2.3, 7.6.2.4, and 7.6.2.5 (ASTM E1820-2013 A9.8 Qualification of Data) standards. The results show that the values of J 0.2BL with different reliabilities all satisfied the requirements of this standard. As the J 0.2BL is independent of the dimensions of the specimens, in engineering applications, it is usually taken as J IC . The J IC value determined while considering the reliability is marked as J IC,p as shown in table 7.

.1. Comparison on the J resistance curves
In addition to the ISO 12135:2016 standard, the ASTM E1820-18ae1 standard is also a standard for the testing of the fracture toughness. They are similar, but there are differences in the data processing procedures. Guan and Li presented a comparison on these test methods [25]. However, the comparison was performed on different  samples. The data scatter of the fracture toughness likely increased the differences between the results obtained by the two standards. Thus, in this paper, the comparison was conducted using the same sample to avoid the influence of randomness.
In the ASTM E1820-18ae1 standard, the construction line is determined based on the assumption of the elastic-perfectly plastic material [25]. It is expressed as follows: where M2, and it can be determined from the (J, Δa) data as well. Here, M=2. σ Y is the effective yield strength, which is an assumed value of the uniaxial yield strength. It can be calculated as follows: where σ YS is the 0.2% offset yield strength, and σ TS is the ultimate tensile strength.
The qualified data region in the ASTM E1820-18ae1 standard is confined by J limit , Δa min , and Δa limit . J limit can be calculated as follows: Two parallel lines with an offset of 0.15 and 1.5 mm from the construction line were plotted and marked as l 0.15 and l 1.5 , respectively, as shown in figure 9. They are plotted as two red solid lines.
Then, similar to the ISO 12135:2016 standard, the data that fell in the region between the l 0.15 and l 1.5 were used to determine a linear regression line of the following form: If k=1.0 mm and α in equation (1) is 0, equation (19) can be transformed into equation (1). However, for the two standards, the data used to fit the regression line were slightly different, so the regression line was slightly different. For the ASTM standard, the regression lines of the specimens were expressed by different types of lines, and their corresponding equations are shown in table 5.
As shown in table 5, the J-Δa equations of the same sample obtained by the two standards were different. As the ASTM had higher refinements on the qualified data, the J-Δa curves of the specimens obtained by the ASTM were below those obtained by the ISO. However, this does not indicate that the material fracture toughness J IC obtained by the ASTM is more conservation, because the value of the fracture toughness J IC is decided by the regression and construction lines. The construction line in the ISO standard has a bigger slope and influences the value of J IC as well.

Comparison of J resistance curves considering the reliability
A similar method to that used in the section 3 is applied in this part. The J-Δa curves for different reliabilities were obtained as shown in figure 8. They are expressed by the dashed line with hollow symbols. In the ASTM standard, the qualified data are limited to the region confined by the l 0.15(ASTM) and l 1.5(ASTM) , which are expressed by the black dashed-dotted lines in this figure.
As shown in the figure, the average J-Δa curves obtained by the different standards almost coincided, while the J-Δa curves considering the reliability were slightly different. When the crack extension was very small, the (J-Δa) p(ISO) and (J-Δa) p(ASTM) were slightly different. As the crack extension Δa increased, when the crack extension Δa>0.6 mm, the differences between (J-Δa) p(ISO) and (J-Δa) p(ASTM) gradually increased.
This was mainly because the average J-Δa curve was the average value of four specimens without considering the influence of the sample variance S 2 . The average values of J lg( ) obtained by the different standards were almost the same.
Because the J resistance equations obtained using the two standards were different, the sample variances calculated by the J resistance equations were different as well. Therefore, when the reliability was taken into account, the (J-Δa) p curves obtained by the two standards were different. As the crack extension Δa increased, the differences between the value of J obtained by the two types of J resistance equations increased exponentially. The differences shown in figure 8 are increasingly evident.
Similar to the ISO 12135:2016 standard, based on the equation of (J-Δa) p(ASTM) , the material fracture toughness J IC can also be obtained by the ASTM. The result was also validated by the requirement described in the ASTM E1820-2013 A9.8 Qualification of Data standard. The values of J IC obtained by the different standards are shown in table 7. The value of Δa IC is the crack extension length corresponding to J IC in figure 8.
As shown in table 7, the value of the fracture toughness J IC obtained by the ASTM standard was always larger than the corresponding value obtained by the ISO 12135:2016 standard. Because the construction lines of the two standards are based on two different material constitutive models, if the material constitutive model is more aligned with the stress-strain curve of the ideal elastic plastic material, the standard ASTM is recommended; otherwise the ISO is recommended. For 30CrNi2MoVA, the true stress-strain relationship was more aligned with a power function, and the the ISO 12135:2016 standard calculation of the fracture toughness was more suitable. Compared with the standard ASTM, this is a more conservative result, and the error of the fracture toughness J IC obtained by the two standards was about 5%. However, for the different reliabilities, the change of the error is subtle.
In addition to the material 30CrNi2MoVA, the calculation method on the fracture toughness of specific reliability is applied into three more other materials, and the results are shown in table7. The J-Δa data of these materials are respectively derived from the refrence [26][27][28]. As shown in table 7, the fracture toughness J IC obtained by the ISO standard is usually smaller than that obtained by the ASTM standard. Similar to the material 30CrNi2MoVA, the constitutive model of Alloy6061-T6511 and Ductile Cast Iron are not the ideal elastic plastic material. For them, the ISO standard is recommended. But for the material APL 5 L X80, the ASTM standard is more suitable. As described in table 7, the error of the fracture toughness measured by the two standards can even reach up to 11.63%.That means for different materials, reasonable choices of the test standard are of great importance.
For 30CrNi2MoVA and Ductile Cast Iron, the change of the error between the different reliabilities is subtle. But for the other two materials, APL 5 L X80 and Alloy6061-T6511, the change of the error is great. It is mainly related with the degree of the dispersion of the J-Δa data. For different specimens used in the test of one material, the greater are the differences between the data of different samples, the larger is the value of the sample variance S2, then the smaller is the sample average J lg( ) calculated by equation (6), and more significant are the differences between the fracture toughness JIC of different specific reliability. That means, higher scatter of the experimental data is supposed to generate a more conservative fracture toughness of the specific reliability obtained by this method.

Conclusions
Fracture toughness, as a crucial parameter, is used widely in the research of fatigue life predictions and failure analysis. The test of the fracture toughness is one of the most common actions in such research. Because of the randomness in the test, the fracture toughness J IC considering the reliability was studied. In addition, the J IC values obtained by the standard ISO and ASTM were compared: (1) A new data processing method for the fracture toughness considering the reliability was proposed. The twodimensional, one-side tolerance factor was applied to calculate the J IC for a given reliability. This method can take advantage of the current test data and historical data to decide the J IC,p value and can reduce experimental costs. By taking the reliability of the fracture toughness J IC into the account, a more conservative result was obtained, and a more reliable guarantee of the performance of the equipment was provided.
(2) Based on the comparison of the data processing method for the fracture toughness, J IC , presented in the ISO and ASTM standards, their differences were mainly due to the range of qualified data, function of the construction line, and equation of the regression line. These differences resulted in different fracture toughness, J IC , values using the two standards, and the J IC tested by the ISO was smaller. Because the construction lines of the two standards are based on two different material constitutive models, if the material constitutive model is more aligned with the stress-strain curve of the ideal elastic plastic material, the ASTM standard is recommended; otherwise the ISO standard is recommended. For 30CrNi2MoVA, the true stress-strain relationship is more aligned with a power function, and the ISO 12135:2016 standard calculation of the fracture toughness is more suitable. For different materials, the change of the error between the different reliabilities is mainly related with the degree of the dispersion of the J-Δa data. The higher is the scatter degree of the experimental data, a more conservative value of the fracture toughness of specific reliability can be obtained by this method.