Brittle or QuasiBrittle Fracture of Engineering Materials: Recent Developments and New Challenges
View this Special IssueResearch Article  Open Access
Application of Stress to Predict the Lower Bound Fracture Toughness for Increasing the Test Specimen Thickness in the Transition Temperature Region
Abstract
This work was motivated by the fact that although fracture toughness of a material in the ductiletobrittle transition temperature region exhibits the test specimen thickness (TST) effect on , frequently described as , experiences a contradiction that is deduced from this empirical formulation; that is, = 0 for large TST. On the other hand, our previous works have showed that the TST effect on could be explained as a difference in the outofplane constraint and correlated with the outofplane stress. Thus, in this work, the TST effect on for the decommissioned Shoreham reactor vessel steel A533B was demonstrated from the standpoint of outofplane constraint. The results validated that was effective for describing the decreasing tendency. Because the Shoreham data included a lower bound for increasing TST, a new finding was made that successfully predicted the lower bound of with increasing TST. This lower bound prediction with conquered the contradiction that the empirical predicts = 0 for large TST.
1. Introduction
The cleavage fracture toughness of a material in the ductiletobrittle transition (DBT) temperature region, which is important in the assessment of aging steel structures and reactor pressure vessels, has been known to exhibit test specimen size effects, even when tested using a standardized specimen [1–9]. For example, obtained using a shallow cracked specimen exhibits a higher value than that obtained using a deep cracked specimen. Another known size effect is the test specimen thickness (TST) effect on , hereafter abbreviated as the TST effect on , which is described as () [2, 10]. The two most physically logical explanations in general are the statistical weakest link (SWL) size effect and the loss of the cracktip constraint [2]. Both explanations lead to an increasing toughness with decreasing TST. The difference in obtained with a different planar specimen configuration, including the crack depth [4], has been explained as the differences in the cracktip constraint or the hydrostatic stress triaxiality, which fails to describe [3, 5–9]. However, the TST effect has been explained in terms of the SWL size effect being dominant [6, 11–13], even though does not decrease indefinitely with thickness [6], which contradicts the prediction from the SWL size effect [2].
Based on the above, the authors believed that the contribution of the cracktip constraint to the TST effect on could be demonstrated if the TST effect (especially the bounded nature of with increasing TST) was demonstrated using a series of nonstandard test specimens whose planar configurations are identical but whose thicknesstowidth ratios, , are changed to realize different thickness specimens and if the test results were reproduced using finite element analysis (FEA). This use of nonstandard test specimens was prompted by the inability to predict the bounded nature of using the SWL formulation. This prediction was thought to be enabled by these specimens because the outofplane cracktip constraint will increase and saturate with increasing , but the inplane cracktip constraint will not change. The fracture toughness tests for a series of nonstandard compacttension (CT) and threepointbend (3PB, also named as SE(B) specimen) specimens for 0.55% carbon steel S55C [14–16] and 0.40% carbon chromium molybdenum steel SCM440 [17] validated the noticeable contribution of the outofplane cracktip constraint to the TST effects on , and the constraint parameter stress was demonstrated to be effective for correlating this outofplane cracktip constraint with the TST effects on [14–17]. These results indicated a possibility of correlating the fracture toughness of a test specimen and the cracklike flaws in the structure more accurately by considering .
This work is an extension of our previous studies regarding the point that the contribution of the outofplane cracktip constraint to the TST effect on was demonstrated for the decommissioned Shoreham reactor vessel steel, ASTM A533 Grade B Class 1 (A533B) [1], which is experimentally formulated as [N/mm] = () to describe the decreasing tendency for increasing TST. Because the Shoreham data included a lower bound of for increasing TST, a new finding was made that successfully predicted the lower bound with increasing TST. This lower bound prediction with resolves the contradiction that the empirical predicts for large TST.
2. TST Effect on Described by the Stress
2.1. Stress
In an isotropic linear elastic body containing a crack subjected to symmetric (mode I) loading, the leading two terms in a series expansion of the stress field very near to the crack front are [18] where and are the inplane polar coordinates of the plane normal to the crack front, as shown in Figure 1, and is the local mode I stress intensity factor (SIF) at location A. Here is the direction formed by the intersection of the plane normal to the crack front and the crack plane. The terms and are the amplitudes of the secondorder terms in the threedimensional series expansion of the crack front stress field in the and directions, respectively.
2.2. TST Effect on Described by Stress
In our previous works [14, 15, 17], the following relationships were obtained for 0.55% carbon steel S55C [14, 15] and 0.40% carbon chromium molybdenum steel SCM440 [17] with both CT and 3PB specimens:
The object of these works was to demonstrate that the outofplane cracktip constraint has a noticeable contribution to the TST effect on and that the TST effect can be correlated with a mechanical parameter (expressing the outofplane cracktip constraint).
Because the bounded nature of with increasing TST could not be realized with the tested specimens of thicknesstowidth ratios , 0.4, and 0.5, the tested results with large were searched in the published documents, and the decommissioned Shoreham reactor vessel steel data [1] were found to fulfill our requirement. In the following, Shoreham’s data were compiled to validate the relationship (: material constant) and, in particular, to correlate the bounded nature of for increasing TST with .
3. Compilation of the Decommissioned Shoreham Reactor Vessel Steel Fracture Toughness Test Data from the Standpoint of OutofPlane Constraint
3.1. Prediction of a Lower Bound of for Increasing TST with
From our recent elastic FEA results for the nonstandard 3PB specimen with various values, as shown in Figure 2(a), the midplane normalized in the form of exhibited a strong dependence on [17]. was negative for , whereas it was positive and approached () for increasing TST. The negative recompiled in the loglog form, as shown in Figure 2(b), exhibited for and a bounded nature for in an engineering sense. This engineering onset of the bounded nature of was defined as the bounded value . Because the SIF corresponding to the fracture load exhibited a small change with TST [14–17], it was thought that the experimental formulation (, : material constants) together with could predict the lower bound value of with increasing TST as follows:
(a)
(b)
3.2. Compilation of the Decommissioned Shoreham Reactor Vessel Steel Fracture Toughness Test Data from the Standpoint of the OutofPlane Constraint
To determine whether the relationship is valid for other materials and especially whether the lower bound can be predicted by , the decommissioned Shoreham reactor vessel steel [1] A533B was selected in this work because a large amount of fracture toughness test data for A533B with various thickness 3PB specimens at a common temperature −91°C (located in the DBT temperature region) was published. A more detailed description for the fracture toughness tests can be found in [1].
Here, the fracture toughness test data for 3PB specimens with width and 50.8 mm whose thicknesses = 8, 15.9, 31.8, and 63.5 mm (thicknesstowidth ratio = 0.157~2.5) were recompiled from the published results [1] on the standpoint of the outofplane cracktip constraint. Although the eight replicate fracture toughness test results reported in [1] for these 3PB specimens were considered to be valid overall, some of the individual datum still appeared to deviate greatly from the remainder in each set. Considering the fact that the scatter from eight replicate tests always exceeded the guideline value as given in ASTM E1921 [19], we thought it was necessary to recompile these test results because the impact of the apparent deviated datum for each set was considered nonnegligible in studying the TST effect on the cleavage fracture toughness. Therefore, the cases with maximum and minimum values were excluded, with the test results of the remaining cases summarized in Tables 1 and 2.
 
Published results in [1]. 
 
Published results in [1]. 
The in the tables was obtained as the SIF corresponding to the fracture load from the following equation in ASTM E1921 [19]:
Here, is the support span, and is a function of , which is given in the ASTM E1921 [19].
in the table is the fracture toughness in terms of the SIF. was calculated from as = , where the value of Young’s modulus of GPa and the value of Poisson’s ratio of were used, as specified in [20]. , which reflects the fracture load and the actual crack length, was calculated from the solutions of elastic FEA, as summarized in the Appendix. and are the average and standard deviation of each value, respectively. is a reference value that was used to represent the magnitude of the data scatter.
It is seen from Tables 1 and 2 that, except for the case of mm with a very thin thickness mm (%), the reference value of was in the range from 33.1% to 45.6% for the selected specimens, which satisfied the guideline for given in ASTM E1921 [19] for . Here the guideline for 2 Σ/μ is 56(1–20/μ)% with the range from 40.7% to 47.4% for the data in Tables 1 and 2. As a result, it could be concluded that the scatter in the data of the selected specimens summarized in the tables was acceptable in an engineering sense.
One interesting fact was that the change in , that is, the SIF for the fracture load , exhibited a relatively small dependence on , although a significant change in the fracture toughness was observed. The average for each was in the range from 67.2 to 80.7 MPa m^{1/2} for mm and 75.8 to 95.7 MPa m^{1/2} for mm. This result was similar to the experience with S55C [14–16] and SCM440 [17], which validated one of the assumptions used to predict the lower bound of for large TST proposed in Section 3.1.
The relationship between and for A533B is shown in Figure 3; note that reflects the fracture load and the actual crack length for each , as summarized in Table 1 and 2. The solid marks represent the average for each . The difference in was distinguished by the color of the marks. As shown in Figure 3, all the data in Tables 1 and 2 are fitted to the power law expression for A533B tested using 3PB specimens at −91°C. seemed to be bounded for MPa. The bounded value of in Figure 3 for the case of mm was obtained from Table 1 as an average for the specimens of and 2.5. For the case of mm, the bounded was obtained from Table 2 as an average for .
On the other hand, if the method to predict the lower bound for increasing TST proposed in Section 3.1 is applied, for the case of mm as an example, first is calculated with for the case of and MPa m^{1/2} (the averaged SIF for B/W = 0.315~1.25 was used from Table 1, considering the fact that exhibited a very small dependence on TST) as . Then, the lower bound is predicted from (3) as , and it was close to experimental average 36.0 N/mm. In case of mm, by the same method, was obtained and was also very close to the experimental average 27.9 N/mm.
In summary, the TST effect on of A533B could be described by , as for . In addition, the lower bound value of = 31.8 N/mm was obtained for mm and = 27.2 N/mm for mm; both of them were close to the experimental average value, which indicated that can successfully predict the bounded nature of .
4. Discussion
In this work, the TST effect and the bounded nature of observed for the decommissioned Shoreham reactor vessel steel, A533B, at −91°C, which is in the DBT range [1], were compiled by stress in the general form of (5). In (5), the similar power law relationship between and was also valid for the combination of S55C [14, 15] and SCM440 [17] tested using both CT and 3PB specimens. In addition, , which seemed to be useful for predicting the bounded nature of for S55C [16], has also been proven to be valid for A533B. In these empirical equations, the TST effect and the bounded nature of were described with a single outofplane elastic parameter taken at the specimen midplane. Although the depicted relationship between the fracture toughness of a material and must be validated for other materials and other types of test specimen configurations, using as a relevant outofplane constraint parameter is definitely worthy of further investigation.
It could be argued that the relationship (Figure 3) is similar to the formulation deduced from the SWL model, but no more than what is predicted by the SWL model [2], because first approaches to 0 for large TST (Note: with increase in TST for 3PB specimen, negative first increases, crosses 0 and converges to ). As Anderson et al. indicated, as a contradiction of the SWL model, the “fracture toughness does not decrease indefinitely with thickness [6].” On the point that exhibits a saturating tendency for large TST, has also been proven to be able to predict the bounded behavior of (Figure 3). The advantage of using is that has the characteristic to not only describe the TST effect on but to also predict the bounded nature of . This advantage of successfully avoids the contradiction deduced from the SWL model; that is, for .
ASTM E1921 [19] presents a method to adjust for CT’s TST change by considering the empirical relationship , under the assumption that 1inch (1T) thickness CT toughness data exists. The presented method in this paper for a 3PB specimen can be generally applied to any type of test specimens, if a curve similar to Figure 2 is obtained. The fact that 1T CT test data are not necessary for our method can help practitioners in their works.
When the proposed general formulation of (3) is practically used for determining the lower bound of for a specific material tested with a fracture toughness test specimen, the material constants and should be first determined by conducting measurements on at least two differentsized specimens. Nevertheless, if measurements on only one size of specimen are conducted, (3) can also be simply but not accurately applied for predicting the lower bound fracture toughness just by assuming in the relationship for that one size of specimen considered, because the material constant has been verified for the materials S55C and SCM440 tested with both CT and 3PB specimens [14, 15, 17]; in addition, this work validated that the approximated which is close to was applicable for the material A533B tested using 3PB specimens.
The normalized stress solutions used in this work were taken at the specimen midplane. It is true that these values are distributed in the specimen thickness direction [21]. There are many possibilities to treat this 3D effect, but, considering the fact that the fracture tends to initiate at the specimen midplane, the values at the specimen midplane were chosen to represent the characteristic intensity of these values.
5. Conclusions
This paper demonstrated for the decommissioned Shoreham reactor vessel steel A533B [1] that the outofplane cracktip constraint has a noticeable contribution to the TST effect on and that the magnitude of this outofplane cracktip constraint can be described by the elastic stress. The experimental expression of the TST effect on using stress, which was proposed for 0.55% carbon steel S55C and 0.40% carbon chromium molybdenum steel SCM440 with both CT and 3PB specimens in our previous work [14, 15, 17], was shown to be a correct description for A533B. In concrete, the experimental relationship for A533B was compiled as () to describe the decreasing tendency for increasing TST. Because the Shoreham data included a lower bound for increasing TST, a new discovery was that successfully predicted the lower bound of with increasing TST. This lower bound of prediction with resolved the contradiction that the empirical predicts for large TST.
Appendix
The normalized solutions used to calculate in Tables 1 and 2 were obtained from the elastic FEA. In the present FEA, all the 3PB specimen dimensions were specified in accordance with those recorded in [1], and the material properties were set to be consistent with those specified in [20] for A533B.
The typical FEA model of the 3PB specimen used in the present elastic analysis is shown in Figure 4, with the details for the generated mesh being summarized in Table 3. The details of the elastic FEA procedure can be found in our recent work [17]. The normalized stress, , at the specimen midplane is summarized in Table 4, which is in a good agreement with the interpolated solutions from our previous results [22].


(a) Global mesh
(b) Details of cracktip for elastic analysis
Nomenclature
Specimen thickness  
:  Material constant (see (3)) 
:  Young’s modulus 
:  integral 
and :  Fracture toughness and its average 
:  Lower bound fracture toughness 
:  Local mode I stress intensity factor (SIF) 
:  Fracture toughness 
:  SIF corresponding to fracture load 
:  Fracture load 
:  Crack tube radius 
:  Support span for 3PB specimen 
and :  stresses 
:  stress corresponding to fracture load 
:  Bounded value of stress 
:  Specimen width 
:  Crack length 
:  Inplane polar coordinates 
:  Cracktip local coordinates () 
:  Singular element size 
:  Standard deviation 
, :  Normalized forms of the stresses 
:  Bounded value of 
:  Material constant (see (3)) 
:  Average value 
:  Poisson’s ratio 
:  Stress components (). 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported in part by JSPS KAKENHI Grant no. 24561038. Their support is greatly appreciated.
References
 H. J. Rathbun, G. R. Odette, T. Yamamoto, and G. E. Lucas, “Influence of statistical and constraint loss size effects on cleavage fracture toughness in the transition—a single variable experiment and database,” Engineering Fracture Mechanics, vol. 73, no. 1, pp. 134–158, 2006. View at: Publisher Site  Google Scholar
 K. Wallin, “The size effect in K_{IC} results,” Engineering Fracture Mechanics, vol. 22, no. 1, pp. 149–163, 1985. View at: Google Scholar
 M. Nevalainen and R. H. Dodds Jr., “Numerical investigation of 3D constraint effects on britttle fracture in SE(B) and C(T) specimens,” International Journal of Fracture, vol. 74, no. 2, pp. 131–161, 1995. View at: Google Scholar
 R. H. Dodds Jr., T. L. Anderson, and M. T. Kirk, “A framework to correlate a/W ratio effects on elasticplastic fracture toughness (J_{c}),” International Journal of Fracture, vol. 48, no. 1, pp. 1–22, 1991. View at: Publisher Site  Google Scholar
 M. T. Kirk, R. H. Dodds, and T. L. Anderson, “An approximate technique for predicting size effects on cleavage fracture toughness (J_{c}) using the elastic T stress,” in Fracture Mechanics, J. D. Landes, D. E. Mccabe, and J. A. M. Boulet, Eds., vol. 24, pp. 62–86, American Society for Testing and Materials, Philadelphia, Pa, USA, 1994. View at: Google Scholar
 T. L. Anderson, D. Stienstra, and R. H. Dodds, “A theoretical framework for addressing fracture in the ductilebrittle transition region,” in Fracture Mechanics, J. D. Landes, D. E. Mccabe, and J. A. M. Boulet, Eds., vol. 24, pp. 186–214, American Society for Testing and Materials, Philadelphia, Pa, USA, 1994. View at: Google Scholar
 T. J. Theiss and J. W. Bryson, “Influence of crack depth on the fracture toughness of reactor pressure vessel steel,” in Constraint Effects in Fracture, E. M. Hackett, Ed., pp. 104–119, American Society for Testing and Materials, Philadelphia, Pa, USA, 1993. View at: Google Scholar
 J. D. G. Sumpter, “An experimental investigation of the T stress approach,” in Constraint Effects in Fracture, E. M. Hackett, Ed., pp. 495–502, American Society For Testing and Materials, Philadelphia, Pa, USA, 1993. View at: Google Scholar
 M. T. Kirk, K. C. Koppenhoefer, and C. F. Shih, “Effect of constraint on specimen dimensions needed to obtain structurally relevant toughness measures,” in Constraint Effects in Fracture, E. M. Hackett, Ed., pp. 79–103, American Society for Testing and Materials, Philadelphia, Pa, USA, 1993. View at: Google Scholar
 J. P. Petti and R. H. Dodds Jr., “Coupling of the Weibull stress model and macroscale models to predict cleavage fracture,” Engineering Fracture Mechanics, vol. 71, no. 1314, pp. 2079–2103, 2004. View at: Publisher Site  Google Scholar
 K. Wallin, T. Saario, and K. Torronen, “Statistical model for carbide induced brittle fracture in steel,” Metal Science, vol. 18, no. 1, pp. 13–16, 1984. View at: Google Scholar
 T. Lin, A. G. Evans, and R. O. Ritchie, “A statistical model of brittle fracture by transgranular cleavage,” Journal of the Mechanics and Physics of Solids, vol. 34, no. 5, pp. 477–497, 1986. View at: Google Scholar
 C. Ruggieri, “Probabilistic treatment of fracture using two failure models,” Probabilistic Engineering Mechanics, vol. 13, no. 4, pp. 309–319, 1998. View at: Google Scholar
 T. Meshii and T. Tanaka, “Experimental T_{33}stress formulation of test specimen thickness effect on fracture toughness in the transition temperature region,” Engineering Fracture Mechanics, vol. 77, no. 5, pp. 867–877, 2010. View at: Publisher Site  Google Scholar
 T. Tanaka and T. Meshii, “Formulating test specimen thickness effect on fracture toughness with T_{33}stress: case of 3PB test specimen,” in Proceedings of the ASME Pressure Vessels and Piping Division/KPVP Conference (PVP '10), pp. 1213–1219, Seattle,Wash, USA, July 2010. View at: Google Scholar
 T. Meshii, K. Lu, and R. Takamura, “A failure criterion to explain the test specimen thickness effect on fracture toughness in the transition temperature region,” Engineering Fracture Mechanics, vol. 104, pp. 184–197, 2013. View at: Google Scholar
 T. Meshii and T. Tanaka, “Framework to correlate test specimen thickness effect on fracture toughness with T_{33}stress,” in Proceedings of the ASME Pressure Vessels and Piping Conference, pp. 1109–1117, Baltimore, Md, USA,Baltimore, Md, USA, 2011. View at: Google Scholar
 T. Nakamura and D. M. Parks, “Determination of elastic Tstress along threedimensional crack fronts using an interaction integral,” International Journal of Solids and Structures, vol. 29, no. 13, pp. 1597–1611, 1992. View at: Google Scholar
 ASTM, E192110 Standard test method for determination of reference temperature, T_{0}, for ferritic steels in the transition range. Annual Book of ASTM Standards. American Society for Testing and Materials, Philadelphia, Pa, USA, 2010.
 H. J. Rathbun, G. R. Odette, M. Y. He, and T. Yamamoto, “Influence of statistical and constraint loss size effects on cleavage fracture toughness in the transition—a model based analysis,” Engineering Fracture Mechanics, vol. 73, no. 18, pp. 2723–2747, 2006. View at: Publisher Site  Google Scholar
 J. Qu and X. Wang, “Solutions of Tstresses for quarterelliptical corner cracks in finite thickness plates subject to tension and bending,” International Journal of Pressure Vessels and Piping, vol. 83, no. 8, pp. 593–606, 2006. View at: Publisher Site  Google Scholar
 K. Lu and T. Meshii, “Threedimensional Tstresses for threepointbend specimens with large thickness variation,” Engineering Fracture Mechanics, pp. 197–203, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Kai Lu and Toshiyuki Meshii. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.