Abstract
The material of primary coolant pipe of nuclear power plants is cast austenitic stainless steel (CASS), which has a coarse crystal structure and its anisotropy causes low permeability of ultrasonic waves and low detectability. Also, the material strength and fracture toughness of CASS are affected by thermal aging depending on ferrite content, chemical compositions and casting process. The ASME Section XI Code Case N-838 provides the target flaw sizes for performance demonstration (PD) examination system of ultrasonic testing (UT) based on probabilistic fracture mechanics (PFM) analysis incorporated a thermal aging model. In this study, the target flaw sizes for Japanese pressurized water reactor (PWR) plants were investigated by using the PFM analysis code “PREFACE.” Thermal aging prediction models for tensile strength and fracture toughness of CASS were incorporated into the PREFACE Code depending on not only ferrite content and chemical compositions but also the casting process, centrifugal casting, or static casting. Those parameters were input as probabilistic variables, which were based on the material database from the Japanese PWR plants. The PFM analysis revealed that static CASS pipes have smaller target flaw sizes than centrifugal CASS pipes and the obtained target flaw sizes were very close to those in ASME Code Case N-838.
Introduction
It is well known that cast austenitic stainless steel (CASS) used for primary coolant pipes of pressurized water reactor (PWR) plants is affected by thermal aging [1–3]. That is, thermal aging causes increases in tensile strength and decreases in fracture toughness during plant operation, which affects failure modes and allowable flaw depths of the CASS pipes. The ASME Section XI Code Case N-838 (CC N-838) [4] provides the target flaw sizes for Nondestructive Inspection in performance demonstration (PD) for ultrasonic examination system, which is based on fracture evaluations of flaws taking into consideration the effect of thermal aging.
A flaw is detected by volumetric examination during plant operation such as ultrasonic testing (UT). However, the coarse crystalline structure of CASS and its anisotropy cause low permeability of ultrasonic waves and low detectability. In recent years, the introduction of PD for ultrasonic examination systems in Japan was being examined by the Japanese Society for Nondestructive Inspection following the United States. For the introduction of the PD examination system, it is necessary to determine the target flaw sizes for the nondestructive testing. As a method for determining the target flaw sizes corresponding to the maximum flaw size without fracture, deterministic fracture mechanics is one candidate for application. However, there are several overconservatisms in the procedure, for example, use of fatigue flaw growth rate with 95% upper bound or a J-R curve from a high constraint laboratory specimen. As a result, the obtained target flaw sizes are small and difficult to detect by UT. Therefore, CC N-838 applied the probability fracture mechanics (PFM) [5] procedure to determine the target flaw sizes instead of deterministic fracture mechanics.
In the previous studies [6–8], a PFM analysis code, “PREFACE” was applied, and determination of target flaw sizes was attempted. Hayashi et al. [7] performed sensitivity analyses on the ferrite contents only in the fully thermal aged condition using PREFACE to obtain the target flaw sizes in the fully thermal aged condition. As a result, the target flaw sizes became the smallest at the minimum ferrite content, which is a condition of relatively larger fracture toughness and lower tensile strength. This was because the lower tensile strength caused plastic collapse. This report suggested that sensitivity analyses in arbitrary thermal aging conditions were needed to determine the actual minimum target flaw sizes.
Nishi et al. [8] performed sensitivity analyses using PREFACE with parameters of the ferrite contents and the thermal aging time. As a result, it was confirmed that the target flaw sizes became a minimum at lower ferrite content and relatively short thermal aging time, which is the condition with higher toughness and lower tensile strength than the initial values.
In the previous study, there were several points to be improved in the analysis conditions of PREFACE. First, the probabilistic distributions of the mechanical properties used in the analysis did not correspond to those of the CASS pipe of operating plants. The probabilistic distributions of each coefficient in prediction models for tensile properties and fracture toughness were applied. As a result, the mechanical properties had larger distribution than the actual distribution. Second, the correlation between tensile strength and fracture toughness was ignored and the target flaw sizes were calculated under the condition of low tensile strength and also low fracture toughness. In those cases, PFM analysis conditions were considered to be more conservative than the actual plant conditions. In order to set more realistic target flaw sizes, it is necessary to use probabilistic distributions of the mechanical properties considering correlations between tensile strength and fracture toughness of aged CASS.
On the other hand, there may be a less conservative condition in the previous studies using the PREFACE code. CASS straight parts of primary coolant pipe for PWR were manufactured by centrifugal casting and the other parts, such as elbows, were manufactured by static casting. Chopra [9] indicated that the thermal aging degradation behavior of CASS depends on the casting process. However, the existing thermal aging prediction models, the H3T, and true stress–true strain curve prediction method (TSS) models proposed by Kawaguchi et al. [10] are based on the database of the material properties of CASS manufactured by both casting processes. Similarly, Griesbach et al. [5] did not consider the effect of the differences of casting processes when they set the probabilistic distributions of the mechanical properties as random variables.
Momotani et al. [11] applied the H3T model and the TSS model to aged CASS manufactured by centrifugal casting and static casting, separately, and pointed out that the TSS model sometimes predicted much higher values for tensile properties than the actual aged material of the static CASS. In this case, the fracture evaluations of static CASS using the tensile properties predicted by the TSS model will be unconservative. Then, they proposed a modified TSS model, which is applicable to static CASS. By applying the different models to predict the tensile properties for centrifugal CASS and static CASS, it is expected that reasonable target flaw sizes will be obtained for each casting process.
In this study, PFM analyses for centrifugal and static CASS were carried out in order to obtain more reasonable target flaw sizes for the introduction of PD system in Japan. In the PFM analysis, more reasonable probabilistic distributions of mechanical properties were set by introducing the correlation between them and the lower limits of the tensile properties.
Approach for Determination of the Target Flaw Sizes
In order to determine the target flaw sizes for the PD examination system, PFM analysis was performed. The target flaw sizes were defined as the maximum flaw sizes, which had sufficiently low probability of fracture. In the PFM analysis, material properties were set as random variables. For more reasonable target flaw sizes, the following analysis conditions were adopted. First, to take into consideration the effect of the casting process on the material properties, thermal aging prediction methods were used for each casting process, centrifugal and static. Second, the probabilistic distributions of material properties were set reasonably based on residual analysis of the thermal aging CASS database (Refer to “Prediction Methods for the Mechanical Properties”). In addition, the correlations between the mechanical properties and the lower limits of the tensile properties were considered (Refer to “Advanced Sampling Methods for Mechanical Properties”). Finally, analysis conditions for Japanese CASS pipe and their chemical conditions in the PFM analysis were selected conservatively from among actual Japanese PWR plants (Refer to “Selection of Analysis Conditions for Japanese CASS Pipe”).
Methodology of Probabilistic Fracture Mechanics Analysis
The fracture analysis was carried out by assuming a semi-elliptical surface flaw in the circumferential or axial direction on the inner surface of the pipe as shown in Fig. 1. In the analysis, the PFM analysis code, PREFACE [6–8], was used. The analysis flow in PREFACE is shown in Fig. 2. In PREFACE, the mechanical properties, i.e., tensile properties and fracture toughness, are set as random variables, and the failure probability is calculated by the Monte Carlo method. Flaw sizes, loads, chemical compositions, and thermal aging condition are set as fixed parameters. Each analysis case calculates failure probabilities by changing the flaw depth. The target size is defined as the flaw depth at which the failure probability is equal to a specified value (10−6 or 10−4).

Relationship between predicted flow stress, σf,p, and prediction error, σf-σf,p for centrifugal CASS
Prediction Methods for the Mechanical Properties.
In the fracture evaluation, the tensile properties (σy, σf, and Stress–Strain curve) and the fracture toughness (JIc, J6, and J-R curve) were used as random variables. The probability distribution of each random variable was determined based on the prediction error for test data. The details are as follows.
Ferrite Content and Chemical Compositions.
The H3T model and TSS model were functions of the ferrite content and the chemical compositions. The ferrite content corresponding to the chemical composition was calculated based on ASTM A800 [12]. The ferrite content and the chemical composition were set as fixed parameters in PREFACE.
Tensile Properties.

Relationship between predicted flow stress, σf,p, and normalized prediction error, (σf-σf,p)/σf,p for centrifugal CASS

Relationship between predicted flow stress, σf,p, and standard deviations of prediction error for centrifugal CASS: (a)prediction error and (b) normalized prediction error
Note that N (μ, S2) is a normal distribution with the mean of μ and the standard deviation of S. Table 1 shows the standard deviations of the normalized prediction errors of yield stress, Sσy, and flow stress, Sσf.
Probabilistic distributions and standard deviations for random variables
Standard deviation, Sσy, Sσf, SσSS, SJIc, and SJ6 | |||
---|---|---|---|
Random variable | Probabilistic distribution | Centrifugal CASS | Static CASS |
σy/σy,p | Normal | 0.0665 | 0.1304 |
σf/σf,p | Normal | 0.0351 | 0.0815 |
σSS | Normal | Refer to Table 2 | |
JIc | Lognormal | 0.2396 | |
J6 | Lognormal | 0.1781 |
Standard deviation, Sσy, Sσf, SσSS, SJIc, and SJ6 | |||
---|---|---|---|
Random variable | Probabilistic distribution | Centrifugal CASS | Static CASS |
σy/σy,p | Normal | 0.0665 | 0.1304 |
σf/σf,p | Normal | 0.0351 | 0.0815 |
σSS | Normal | Refer to Table 2 | |
JIc | Lognormal | 0.2396 | |
J6 | Lognormal | 0.1781 |
Date table for sampling of S-S curve
No. | σf,s (MPa) | Strain point | 1 | 2 | 3 | 4 | 5 | 6 | 7 |
---|---|---|---|---|---|---|---|---|---|
1 | 281 | εSS (%) | 0.2 | 0.5 | 1 | 3 | 6 | 12 | 20 |
σSS (MPa) | 134 | 166 | 185 | 227 | 274 | 366 | 473 | ||
SσSS (MPa) | 10.8 | 12.5 | 12.1 | 15.9 | 19.9 | 26.6 | 34.0 | ||
2 | 318 | εSS (%) | 0.2 | 0.5 | 1 | 2 | 4 | 8 | 15 |
σSS (MPa) | 149 | 192 | 222 | 254 | 297 | 371 | 484 | ||
SσSS (MPa) | 13.6 | 13.1 | 13.3 | 15.4 | 16.6 | 18.9 | 22.9 | ||
3 | 377 | εSS (%) | 0.2 | 0.3 | 0.5 | 1 | 2 | 4 | 10 |
σSS (MPa) | 178 | 204 | 235 | 269 | 308 | 360 | 487 | ||
SσSS (MPa) | 13.9 | 14.6 | 13.8 | 17.1 | 21.1 | 25.7 | 36.8 | ||
4 | 423 | εSS (%) | 0.2 | 0.3 | 0.4 | 0.5 | 1 | 2 | 5 |
σSS (MPa) | 182 | 211 | 235 | 255 | 314 | 365 | 462 | ||
SσSS (MPa) | 14.3 | 15.7 | 17.8 | 18.1 | 24.3 | 30.5 | 39.0 |
No. | σf,s (MPa) | Strain point | 1 | 2 | 3 | 4 | 5 | 6 | 7 |
---|---|---|---|---|---|---|---|---|---|
1 | 281 | εSS (%) | 0.2 | 0.5 | 1 | 3 | 6 | 12 | 20 |
σSS (MPa) | 134 | 166 | 185 | 227 | 274 | 366 | 473 | ||
SσSS (MPa) | 10.8 | 12.5 | 12.1 | 15.9 | 19.9 | 26.6 | 34.0 | ||
2 | 318 | εSS (%) | 0.2 | 0.5 | 1 | 2 | 4 | 8 | 15 |
σSS (MPa) | 149 | 192 | 222 | 254 | 297 | 371 | 484 | ||
SσSS (MPa) | 13.6 | 13.1 | 13.3 | 15.4 | 16.6 | 18.9 | 22.9 | ||
3 | 377 | εSS (%) | 0.2 | 0.3 | 0.5 | 1 | 2 | 4 | 10 |
σSS (MPa) | 178 | 204 | 235 | 269 | 308 | 360 | 487 | ||
SσSS (MPa) | 13.9 | 14.6 | 13.8 | 17.1 | 21.1 | 25.7 | 36.8 | ||
4 | 423 | εSS (%) | 0.2 | 0.3 | 0.4 | 0.5 | 1 | 2 | 5 |
σSS (MPa) | 182 | 211 | 235 | 255 | 314 | 365 | 462 | ||
SσSS (MPa) | 14.3 | 15.7 | 17.8 | 18.1 | 24.3 | 30.5 | 39.0 |
Fracture Toughness.
Methodology of Fracture Analysis.
The methods of fracture analysis for assumed flaws on the inner surface of pipe were based on the JSME rules on FFS. Two methods were applied in the fracture analysis: Two-parameters method and Limit-load method. Ductile fracture and plastic collapse due to assumed flaws were evaluated by using each method. The details of these two methods are described in Appendix. The safety factors, SF, were not considered in the PFM analysis (SF = 1). The J-integral for the two-parameters method was calculated by the reference stress method [14]. It's important to note that flow stress affects both methods of fracture analysis. In the Two-parameter method, fracture become less likely with higher flow stress, because it causes higher stress of S-S curves, σSS, and lower J-integral. Also, in the Limit-load method, plastic collapse stress become higher with higher flow stress. On the other hand, yield stress affects only Two-parameter method. What these tensile properties have in common is that lower value causes higher fracture probability.
The Failure Probability.
Advanced Sampling Methods for the Mechanical Properties.
In order to determine more reasonable probability distributions of the mechanical properties, new sampling methods were introduced as described below.
Correlation of the Mechanical Properties.
When a material has a high σy, it is likely to have a higher σf. High strength materials generally have lower toughness, which is a negative correlation between the variation in tensile properties and fracture toughness. In previous analyses using the PREFACE code, those properties were independently generated as probabilistic variables despite those kinds of correlations.
To confirm these correlations, relations of each prediction error used for the calculation of the standard deviation of the mechanical properties were arranged. As a typical example, Figs. 6 and 7 show the relationship between the prediction errors in JIc-J6 and σf-JIc of centrifugal CASS. The correlation coefficients in each combination are arranged in Table 3. The combinations of JIc-J6, σy-σf, and σf-σSS gave positive correlation with the correlation coefficients, R, ranging from 0.61 to 0.78. On the other hand, the correlation between the tensile properties and the fracture toughness were relatively weaker, and the values of R ranged from −0.45 to 0.08.
Correlation coefficient of the prediction error in the material properties
Combination of material properties | Casting process | Correlation coefficient | |
---|---|---|---|
(a) σy-σf, σf-σSS and JIc-J6 | |||
Tensile properties | σy/σy,p and σf/σf,p | Centrifugal | 0.61 |
Static | 0.78 | ||
σf/σy,p and σSS/σSS,avg | Common | 0.70 | |
Fracture toughness | log(JIc) and log(J6) | Common | 0.78 |
Combination of material properties | Casting process | Correlation coefficient | |
---|---|---|---|
(a) σy-σf, σf-σSS and JIc-J6 | |||
Tensile properties | σy/σy,p and σf/σf,p | Centrifugal | 0.61 |
Static | 0.78 | ||
σf/σy,p and σSS/σSS,avg | Common | 0.70 | |
Fracture toughness | log(JIc) and log(J6) | Common | 0.78 |
(b) Tensile property-fracture toughness | ||
---|---|---|
σy and JIc | Centrifugal | −0.16 |
σy and J6 | −0.14 | |
σf and JIc | −0.45 | |
σf and J6 | −0.33 | |
σy and JIc | Static | −0.11 |
σy and J6 | 0.05 | |
σf and JIc | −0.09 | |
σf and J6 | 0.08 |
(b) Tensile property-fracture toughness | ||
---|---|---|
σy and JIc | Centrifugal | −0.16 |
σy and J6 | −0.14 | |
σf and JIc | −0.45 | |
σf and J6 | −0.33 | |
σy and JIc | Static | −0.11 |
σy and J6 | 0.05 | |
σf and JIc | −0.09 | |
σf and J6 | 0.08 |
Figure 8 shows the comparison of prediction error and sampling results considering the correlation by using Eqs. (9) and (10). The variation and the correlation in prediction errors data could be sampled equivalently by considering the correlation, and more reasonable samplings for mechanical properties were possible.

Cumulative probability distributions of yield strength before/after thermal aging under the representative material of static CASS and thermal aging condition
Lower Limits of the Tensile Properties.
Thermal aging increases tensile strength. Caused by larger variations in tensile properties after thermal aging, in the sampling, tensile properties after thermal aging sometimes become smaller than the value assumed by variation before thermal aging, which differs from the actual situation because yield stress cannot decrease by thermal aging. In order to avoid these excessive conservatisms, it was considered that the probabilistic distributions of σy and σf after thermal aging have lower limits by the initial tensile properties before thermal aging, σy0 and σf0.
The standard deviations of the tensile properties before thermal aging are shown in Table 4. These were calculated by using the data of the mill sheets for Japanese PWR plants. Table 4 also shows the standard deviations after aging from Table 1. In static CASS, the standard deviations of σy and σf before thermal aging were smaller than after aging. The cumulative probability of the tensile properties before and after thermal aging under the representative material of static CASS and thermal aging conditions is as shown in Fig. 9. It indicates that if the probability distributions of the tensile properties only after thermal aging are used for sampling in PFM analysis, excessively low tensile properties may be sampled, which is lower than before thermal aging. In the PFM analysis, the cumulative probability density of the tensile properties of static CASS was defined by combinations of cumulative probability densities before and after thermal aging such as the red line in Fig. 9. This simplified combination model of cumulative probability densities does not correspond to actual variation in tensile properties, but it can avoid excessively lower sampling tensile properties and too conservative failure probabilities.

Predicted average values of the material properties of CASS pipe of Japanese PWR plants: (a) centrifugal and (b) static
Standard deviations of normalized prediction error of tensile properties before and after thermal aging
Casting process | Parameter | Thermal aging | Random variable | Standard deviation |
---|---|---|---|---|
Centrifugal | σy | Before aging | σy0/σy0,p | 0.1636 |
After aging | σy/σy,p | 0.0665 | ||
σf | Before aging | σf0/σf0,p | 0.0654 | |
After aging | σf/σf,p | 0.0351 | ||
Static | σy | Before aging | σy0/σy0,p | 0.0726 |
After aging | σy/σy,p | 0.1304 | ||
σf | Before aging | σf0/σf0,p | 0.0398 | |
After aging | σf/σf,p | 0.0815 |
Casting process | Parameter | Thermal aging | Random variable | Standard deviation |
---|---|---|---|---|
Centrifugal | σy | Before aging | σy0/σy0,p | 0.1636 |
After aging | σy/σy,p | 0.0665 | ||
σf | Before aging | σf0/σf0,p | 0.0654 | |
After aging | σf/σf,p | 0.0351 | ||
Static | σy | Before aging | σy0/σy0,p | 0.0726 |
After aging | σy/σy,p | 0.1304 | ||
σf | Before aging | σf0/σf0,p | 0.0398 | |
After aging | σf/σf,p | 0.0815 |
Recalculation of J6,s When JIc,s Was Higher Than J6,s.
In such cases, C1 value is almost same with JIc,s and C2 value is almost zero in Eq. (4).
Analysis Conditions
Selection of Analysis Conditions for Japanese Cast Austenitic Stainless Steel Pipe.
The previous studies [8] reported that failure probability is affected by the ferrite content calculated by chemical compositions and the thermal aging condition. This means that chemical compositions and the thermal aging time should be selected appropriately in order to conservatively determine the target flaw sizes.
Figure 10 shows the average values of σf and JIc predicted by the TSS model or the modified TSS model and the H3T model for static or centrifugal CASS pipe of Japanese PWR plants. In Fig. 10, the thermal aging times were changed in the range of 20,000–525,600 h (corresponding to 2–60 years) at the thermal aging temperature 325 °C. Because the CASS pipe of each plant has different chemical compositions, σf and JIc in Fig. 10 have large variations.
In this study, in order to obtain conservative target flaw sizes for the Japanese PWR plants except plants planed for decommissioning, several CASS pipes, which are assumed to have the highest failure probability, were selected for PFM analysis as described below.
Four CASS pipes from each casting process were selected for the PFM analysis, that is, material C-I to C-IV for centrifugal, and S-I to S-IV for static. Three conditions of the thermal aging time were set, short time (20,000 h), medium time (25,000 or 50,000 h), and long time (525,600 h) at the temperature of 325 °C. In the fracture analysis, ductile fracture occurs at lower fracture toughness (JIc, J6) and lower tensile properties (σy, σf, σSS), and plastic collapse occurs at lower flow stress (σf). Therefore, the failure probability was assumed to become higher with a condition of lower tensile properties and lower fracture toughness. As shown in Fig. 10, σf and JIc values of selected CASS pipes designated by the colored symbols envelope the lower limits of those of for all pipes. It was expected that one of the selected pipes provided the highest failure probability and the smallest target flaw sizes among all CASS pipes. Table 5 summarizes the selected CASS pipes for PFM analysis.
Thermal aging times at 325 °C (h) | |||||||
---|---|---|---|---|---|---|---|
Casting process | CASS pipe material | Ferrite contenta (vol. %) | 20,000b | 25,000c | 50,000c | 525,600d | Notes |
Centrifugal | C-I | 9.08 | C-I-Short | — | C-I-Medium | C-I-Long | C-I-Short has the lowest tensile properties. |
C-II | 9.90 | C-II-Short | — | C-II- Medium | C-II-Long | ||
C-III | 16.38 | C-III-Short | — | C-III- Medium | C-III-Long | ||
C-IV | 22.24 | C-IV-Short | — | C-IV- Medium | C-IV-Long | C-IV has highest ferrite content in Japanese PWR plants. | |
C-IV-Long has the lowest fracture toughness. | |||||||
Static | S-I | 9.44 | S-I-Short | S-I- Medium | — | S-I-Long | S-I-Short has the lowest tensile properties. |
S-II | 10.18 | S-II-Short | S-II-Medium | — | S-II-Long | ||
S-III | 15.85 | S-III-Short | S-III-Medium | — | S-III-Long | ||
S-IV | 20.54 | S-IV-Short | S-IV-Medium | — | S-IV-Long | S-IV-Long has the lowest fracture toughness. |
Thermal aging times at 325 °C (h) | |||||||
---|---|---|---|---|---|---|---|
Casting process | CASS pipe material | Ferrite contenta (vol. %) | 20,000b | 25,000c | 50,000c | 525,600d | Notes |
Centrifugal | C-I | 9.08 | C-I-Short | — | C-I-Medium | C-I-Long | C-I-Short has the lowest tensile properties. |
C-II | 9.90 | C-II-Short | — | C-II- Medium | C-II-Long | ||
C-III | 16.38 | C-III-Short | — | C-III- Medium | C-III-Long | ||
C-IV | 22.24 | C-IV-Short | — | C-IV- Medium | C-IV-Long | C-IV has highest ferrite content in Japanese PWR plants. | |
C-IV-Long has the lowest fracture toughness. | |||||||
Static | S-I | 9.44 | S-I-Short | S-I- Medium | — | S-I-Long | S-I-Short has the lowest tensile properties. |
S-II | 10.18 | S-II-Short | S-II-Medium | — | S-II-Long | ||
S-III | 15.85 | S-III-Short | S-III-Medium | — | S-III-Long | ||
S-IV | 20.54 | S-IV-Short | S-IV-Medium | — | S-IV-Long | S-IV-Long has the lowest fracture toughness. |
Predicted by ASTM A800 [12].
Corresponding to the shortest operating time in Japanese PWR plants.
Medium thermal aging times at which σf and JIc become medium values between at 20,000 h and 525,600 h in thermal aging times.
Corresponding to 60 years operation at 325 °C.
Conditions of Loads, Flaw Lengths, and Allowable Failure Probability.
CC N-838 provides the target flaw sizes for different flaw lengths and load levels for each flaw direction (circumferential or axial) and service condition (A & B or C & D), in which load levels and safety factors are different. In this study, the analysis conditions were set so as to make tables similar to those of the JSME rules on FFS. The flaw length and the load were set as parameters shown in Table 6. The load parameters were defined as (Pm + Pb + Pe)/Sm with changing Pb for circumferential flaws and Ph/Sm with changing Ph for axial flaws. Here, Sm is the design stress, which is defined as one-third of ultimate stress or two-thirds of yield stress, for example,. The PFM analysis were carried out for each combination of load parameters and flaw length parameters.
Probabilistic fracture mechanics analysis conditions of loads and flaw lengths
Flaw direction | |||
---|---|---|---|
Item | Parameter | Circumferential | Axial |
Load | Definition of load parameter | (Pm + Pb + Pe)/Sm | Ph/Sm |
Analysis condition (range) | (Pm + Pb + Pe)/Sm = Service condition A&B: 0.6–1.6Service condition C&D: 1.0–2.6(by changing Pb) | Ph/Sm = Service condition A&B: 0.2–1.6 Service condition C&D: 0.6–2.4(by changing Ph) | |
Constant value | Sm = 116 MPa | Sm = 116 MPa | |
Pm = 58 MPa | |||
Pe = 0 MPa | |||
Flaw length | Definition of flaw length parameter | Flaw angle, 2θ (deg.) | Normalized flaw length, l/√ (Rmt) |
Analysis condition (range) | 2θ = 10–60 | l/√ (Rmt) = 0.1–3.0 |
Flaw direction | |||
---|---|---|---|
Item | Parameter | Circumferential | Axial |
Load | Definition of load parameter | (Pm + Pb + Pe)/Sm | Ph/Sm |
Analysis condition (range) | (Pm + Pb + Pe)/Sm = Service condition A&B: 0.6–1.6Service condition C&D: 1.0–2.6(by changing Pb) | Ph/Sm = Service condition A&B: 0.2–1.6 Service condition C&D: 0.6–2.4(by changing Ph) | |
Constant value | Sm = 116 MPa | Sm = 116 MPa | |
Pm = 58 MPa | |||
Pe = 0 MPa | |||
Flaw length | Definition of flaw length parameter | Flaw angle, 2θ (deg.) | Normalized flaw length, l/√ (Rmt) |
Analysis condition (range) | 2θ = 10–60 | l/√ (Rmt) = 0.1–3.0 |
As shown in Table 7, the allowable failure probabilities in each service condition of A & B or C & D were assumed with reference CC N-838 [4,5]. The sampling numbers in the PFM analysis were set to be N = 5 × 107 or 5 × 105 according to the allowable failure probabilities in the service condition of A & B or C & D.
Probabilistic fracture mechanics analysis conditions of allowable failure probability and sampling number
Service condition | ||
---|---|---|
Parameter | A & B | C & D |
Allowable failure probability | <1 × 106 | <1 × 104 |
Sampling number | 5 × 107 | 5 × 105 |
Service condition | ||
---|---|---|
Parameter | A & B | C & D |
Allowable failure probability | <1 × 106 | <1 × 104 |
Sampling number | 5 × 107 | 5 × 105 |
Common Analysis Condition.
Common analysis conditions such as the pipe geometry, the thermal aging temperature, and the material constants were set as follows:
Outer pipe diameter, Do: 882 mm, Pipe thickness, t: 72.7 mm
Thermal aging temperature, Tage: 325 °C
Young's modules, E: 174,000 MPa, Poisson ratio, ν: 0.3
Analysis Results
The PFM analysis were carried out by the procedure shown in Fig. 2, and the thickness ratio of the maximum flaw depth satisfying the allowable failure probability, which is called the target flaw size, a/t, were calculated for each load and flaw length. Since the target flaw sizes of centrifugal CASS exceeded those of the static CASS under almost all conditions, the details of the latter results are described below.
Circumferential Flaw.
As a typical example of the results of PFM for circumferential flaws, Fig. 11 shows the relationship between target flaw sizes and stress ratios (Pm + Pb + Pe)/Sm of S-IV for 30 deg. flaw angle in the service condition C & D. The thermal aging time at which the target flaw sizes are the minimum values was 525,600 h for S-IV. The same analyses were performed for the other materials, and the thermal aging times which had the minimum target flaw sizes were determined as below:

The target flaw sizes for circumferential flaws in each stress level (service condition: C & D, flaw angle: 2θ = 30 deg., material: S-IV)
S-I and S-II: at 20,000 h (S-I-Short and S-II-Short)
S-III: 20,000 h for intermediate stress ratio, or 525,600 h for high and low stress ratios (S-III-Short or S-III-Long)
S-IV: at 525,600 h (S-IV-Long)
The comparison of the target flaw sizes of each material in these thermal aging times (S- III-Short are shown as a representative) is shown in Fig. 12. The materials which had the minimum target flaw sizes were S-II-Short for medium stress ratios ((Pm + Pb + Pe)/Sm = 1.4–2.0) or S-IV-Long for lower and higher stress ratios ((Pm + Pb + Pe)/Sm = 1.0–1.2 and 2.2–2.6). In almost other flaw angles, the same materials had the minimum target flaw sizes. S-II-Short was the material which had the second lowest tensile property, and S-IV-Long was the material, which had the lowest fracture toughness. The dominant failure modes were ductile fracture in all conditions, and failure probabilities caused by plastic collapse were relatively lower. In the service condition A & B, the dominant thermal aging time and failure modes for minimizing the target flaw sizes were similar to those in the service condition C & D.

Comparison of the target flaw sizes for circumferential flaws in each static CASS pipe material (service condition: C&D, flaw angle: 2θ = 30 deg)
Axial Flaw.
As a typical example of the results of PFM for axial flaws, Fig. 13 shows the result of the target flaw sizes in S-IV with normalized flaw length ℓ/√ (Rmt) = 1.0 in service condition C & D. The trend of the target flaw sizes is similar to Fig. 11, and S-IV-Long had the minimum target flaw sizes. The thermal aging time at which the target flaw sizes were the minimum are shown as below:

The target flaw sizes for axial flaws in each stress ratio (service condition: C&D, nondimensional flaw length: ℓ/√(Rmt) = 1.0, material: S-IV)
S-I and S-II: at 20,000 h (S-I-Short and S-II-Short)
S-III and S-IV: at 525,600 h (S-III-Long and S-IV-Long)
A comparison of the target flaw sizes of each material at these thermal aging times is shown in Fig. 14. The materials having the minimum target flaw sizes were S-IV-Long for lower stress ratios (Ph/Sm = 0.6–1.8) or S-I-Short for higher stress ratios (Ph/Sm = 2.0–2.4). S-I-Short had the lowest tensile properties, and S-IV-Long had the lowest fracture toughness. These tendencies were the same in the other flaw lengths. Plastic collapse was the dominant failure mode only for the higher stress ratio, and ductile fracture was dominant for the other stress ratios.

Comparison of the target flaw sizes for axial flaws in each static CASS pipe material (service condition: C&D, flaw length: ℓ/√(Rmt) = 1.0)
In service condition A & B, S-IV-Long had the minimum target flaw sizes in any flaw lengths and loads, and the dominant fracture mode was ductile unstable fracture.
Target Flaw Size Tables.
Based on the results of the PFM analysis in all conditions, target flaw size tables were arranged. The target flaw sizes were determined as the values corresponding to the combination of the minimum target flaw sizes for each material, for example, S-II-Short and S-IV-Long shown in Fig. 15. Tables 8 and 9 show target flaw size tables of circumferential flaws for centrifugal and static CASS. Tables 10 and 11 show those of axial flaws for centrifugal and static CASS. The centrifugal CASS had larger target flaw sizes than the static CASS in almost conditions. This is because of the higher averages and smaller standard deviations of the tensile properties in centrifugal CASS.
Target flaw size, a/t, of circumferential flaw in centrifugal CASS pipe
Flaw angle, 2θ (deg.) | ||||||
---|---|---|---|---|---|---|
Stress ratio, (Pm + Pb + Pe)/Sm | 10 | 20 | 30 | 40 | 50 | 60 |
(a) Allowable condition: A&Ba | ||||||
1.6 | 0.74 | 0.73 | 0.68 | 0.51 | 0.43 | 0.38 |
1.4 | 0.74 | 0.74 | 0.72 | 0.71 | 0.62 | 0.54 |
1.2 | 0.74 | 0.74 | 0.74 | 0.72 | 0.71 | 0.70 |
1.0 | 0.75 | 0.74 | 0.74 | 0.74 | 0.73 | 0.72 |
0.8 | 0.75 | 0.75 | 0.74 | 0.74 | 0.74 | 0.74 |
0.6 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.74 |
(b) Allowable condition: C&Db | ||||||
2.6 | 0.08 | 0.06 | 0.05 | 0.05 | 0.04 | 0.04 |
2.4 | 0.35 | 0.12 | 0.10 | 0.09 | 0.09 | 0.09 |
2.2 | 0.54 | 0.26 | 0.19 | 0.16 | 0.15 | 0.14 |
2.0 | 0.71 | 0.61 | 0.35 | 0.28 | 0.24 | 0.22 |
1.8 | 0.74 | 0.72 | 0.62 | 0.45 | 0.38 | 0.34 |
1.6 | 0.74 | 0.74 | 0.72 | 0.68 | 0.55 | 0.48 |
1.4 | 0.75 | 0.74 | 0.74 | 0.73 | 0.71 | 0.65 |
1.2 | 0.75 | 0.75 | 0.74 | 0.74 | 0.73 | 0.73 |
1.0 | 0.75 | 0.75 | 0.75 | 0.74 | 0.74 | 0.74 |
Flaw angle, 2θ (deg.) | ||||||
---|---|---|---|---|---|---|
Stress ratio, (Pm + Pb + Pe)/Sm | 10 | 20 | 30 | 40 | 50 | 60 |
(a) Allowable condition: A&Ba | ||||||
1.6 | 0.74 | 0.73 | 0.68 | 0.51 | 0.43 | 0.38 |
1.4 | 0.74 | 0.74 | 0.72 | 0.71 | 0.62 | 0.54 |
1.2 | 0.74 | 0.74 | 0.74 | 0.72 | 0.71 | 0.70 |
1.0 | 0.75 | 0.74 | 0.74 | 0.74 | 0.73 | 0.72 |
0.8 | 0.75 | 0.75 | 0.74 | 0.74 | 0.74 | 0.74 |
0.6 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.74 |
(b) Allowable condition: C&Db | ||||||
2.6 | 0.08 | 0.06 | 0.05 | 0.05 | 0.04 | 0.04 |
2.4 | 0.35 | 0.12 | 0.10 | 0.09 | 0.09 | 0.09 |
2.2 | 0.54 | 0.26 | 0.19 | 0.16 | 0.15 | 0.14 |
2.0 | 0.71 | 0.61 | 0.35 | 0.28 | 0.24 | 0.22 |
1.8 | 0.74 | 0.72 | 0.62 | 0.45 | 0.38 | 0.34 |
1.6 | 0.74 | 0.74 | 0.72 | 0.68 | 0.55 | 0.48 |
1.4 | 0.75 | 0.74 | 0.74 | 0.73 | 0.71 | 0.65 |
1.2 | 0.75 | 0.75 | 0.74 | 0.74 | 0.73 | 0.73 |
1.0 | 0.75 | 0.75 | 0.75 | 0.74 | 0.74 | 0.74 |
Bold: C-IV-Long, Normal: C-III-Medium have the minimum target flaws sizes.
Bold: C-IV-Long, Under line: C-III-Long, Italic: C-III-Short, have the minimum target flaws sizes.
Target flaw size, a/t, of circumferential flaw in static CASS pipe
Flaw angle, 2θ (deg.) | ||||||
---|---|---|---|---|---|---|
Stress ratio, (Pm + Pb + Pe)/Sm | 10 | 20 | 30 | 40 | 50 | 60 |
(a) Service condition: A&Ba | ||||||
1.6 | 0.58 | 0.33 | 0.20 | 0.17 | 0.15 | 0.14 |
1.4 | 0.71 | 0.68 | 0.44 | 0.33 | 0.27 | 0.24 |
1.2 | 0.74 | 0.74 | 0.70 | 0.57 | 0.46 | 0.40 |
1.0 | 0.75 | 0.75 | 0.74 | 0.73 | 0.68 | 0.58 |
0.8 | 0.75 | 0.75 | 0.75 | 0.74 | 0.74 | 0.74 |
0.6 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 |
(b) Service condition: C&Db | ||||||
2.6 | 0.01 | 0.01 | 0.01 | 0.01 | 0.01 | 0.01 |
2.4 | 0.03 | 0.03 | 0.02 | 0.02 | 0.02 | 0.02 |
2.2 | 0.20 | 0.07 | 0.06 | 0.06 | 0.06 | 0.05 |
2.0 | 0.45 | 0.18 | 0.13 | 0.12 | 0.10 | 0.10 |
1.8 | 0.64 | 0.44 | 0.25 | 0.20 | 0.17 | 0.16 |
1.6 | 0.72 | 0.69 | 0.50 | 0.35 | 0.29 | 0.26 |
1.4 | 0.74 | 0.74 | 0.70 | 0.57 | 0.46 | 0.40 |
1.2 | 0.75 | 0.75 | 0.74 | 0.73 | 0.67 | 0.57 |
1.0 | 0.75 | 0.75 | 0.75 | 0.74 | 0.74 | 0.73 |
Flaw angle, 2θ (deg.) | ||||||
---|---|---|---|---|---|---|
Stress ratio, (Pm + Pb + Pe)/Sm | 10 | 20 | 30 | 40 | 50 | 60 |
(a) Service condition: A&Ba | ||||||
1.6 | 0.58 | 0.33 | 0.20 | 0.17 | 0.15 | 0.14 |
1.4 | 0.71 | 0.68 | 0.44 | 0.33 | 0.27 | 0.24 |
1.2 | 0.74 | 0.74 | 0.70 | 0.57 | 0.46 | 0.40 |
1.0 | 0.75 | 0.75 | 0.74 | 0.73 | 0.68 | 0.58 |
0.8 | 0.75 | 0.75 | 0.75 | 0.74 | 0.74 | 0.74 |
0.6 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 |
(b) Service condition: C&Db | ||||||
2.6 | 0.01 | 0.01 | 0.01 | 0.01 | 0.01 | 0.01 |
2.4 | 0.03 | 0.03 | 0.02 | 0.02 | 0.02 | 0.02 |
2.2 | 0.20 | 0.07 | 0.06 | 0.06 | 0.06 | 0.05 |
2.0 | 0.45 | 0.18 | 0.13 | 0.12 | 0.10 | 0.10 |
1.8 | 0.64 | 0.44 | 0.25 | 0.20 | 0.17 | 0.16 |
1.6 | 0.72 | 0.69 | 0.50 | 0.35 | 0.29 | 0.26 |
1.4 | 0.74 | 0.74 | 0.70 | 0.57 | 0.46 | 0.40 |
1.2 | 0.75 | 0.75 | 0.74 | 0.73 | 0.67 | 0.57 |
1.0 | 0.75 | 0.75 | 0.75 | 0.74 | 0.74 | 0.73 |
Bold: S-IV-Long, Under line: S-II-Short, Normal: S-I-Short have the minimum target flaw sizes.
Bold: S-IV-Long, Under line: S-II-Short, Normal: S-IV-Short have the minimum target flaw sizes.
Target flaw size, a/t, of axial flaw in centrifugal CASS pipe
Normalized flaw length ℓ/√(Rmt) | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
Stress ratio, Ph/Sm | 0.1 | 0.2 | 0.4 | 0.6 | 0.8 | 1.0 | 1.5 | 2.0 | 2.5 | 3.0 |
(a) Service condition: A&Ba | ||||||||||
1.6 | 0.55 | 0.55 | 0.55 | 0.26 | 0.19 | 0.17 | 0.14 | 0.13 | 0.12 | 0.12 |
1.4 | 0.70 | 0.70 | 0.70 | 0.69 | 0.44 | 0.33 | 0.24 | 0.21 | 0.20 | 0.19 |
1.2 | 0.74 | 0.74 | 0.74 | 0.74 | 0.70 | 0.59 | 0.39 | 0.33 | 0.30 | 0.28 |
1.0 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.72 | 0.57 | 0.47 | 0.42 | 0.40 |
0.8 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.72 | 0.64 | 0.57 | 0.53 |
0.6 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.74 | 0.73 | 0.69 |
0.4 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 |
0.2 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 |
Normalized flaw length ℓ/√(Rmt) | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
Stress ratio, Ph/Sm | 0.1 | 0.2 | 0.4 | 0.6 | 0.8 | 1.0 | 1.5 | 2.0 | 2.5 | 3.0 |
(a) Service condition: A&Ba | ||||||||||
1.6 | 0.55 | 0.55 | 0.55 | 0.26 | 0.19 | 0.17 | 0.14 | 0.13 | 0.12 | 0.12 |
1.4 | 0.70 | 0.70 | 0.70 | 0.69 | 0.44 | 0.33 | 0.24 | 0.21 | 0.20 | 0.19 |
1.2 | 0.74 | 0.74 | 0.74 | 0.74 | 0.70 | 0.59 | 0.39 | 0.33 | 0.30 | 0.28 |
1.0 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.72 | 0.57 | 0.47 | 0.42 | 0.40 |
0.8 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.72 | 0.64 | 0.57 | 0.53 |
0.6 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.74 | 0.73 | 0.69 |
0.4 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 |
0.2 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 |
(b) Service condition: C&Db | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
2.4 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
2.2 | 0.16 | 0.03 | 0.03 | 0.03 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 |
2.0 | 0.30 | 0.30 | 0.08 | 0.07 | 0.06 | 0.06 | 0.05 | 0.04 | 0.04 | 0.04 |
1.8 | 0.48 | 0.48 | 0.48 | 0.18 | 0.14 | 0.13 | 0.11 | 0.10 | 0.09 | 0.09 |
1.6 | 0.64 | 0.64 | 0.64 | 0.51 | 0.32 | 0.26 | 0.20 | 0.18 | 0.17 | 0.16 |
1.4 | 0.73 | 0.73 | 0.73 | 0.73 | 0.60 | 0.46 | 0.33 | 0.28 | 0.26 | 0.24 |
1.2 | 0.74 | 0.74 | 0.74 | 0.74 | 0.73 | 0.67 | 0.50 | 0.41 | 0.37 | 0.35 |
1.0 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.65 | 0.57 | 0.50 | 0.47 |
0.8 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.70 | 0.64 | 0.60 |
0.6 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.73 |
(b) Service condition: C&Db | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
2.4 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
2.2 | 0.16 | 0.03 | 0.03 | 0.03 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 |
2.0 | 0.30 | 0.30 | 0.08 | 0.07 | 0.06 | 0.06 | 0.05 | 0.04 | 0.04 | 0.04 |
1.8 | 0.48 | 0.48 | 0.48 | 0.18 | 0.14 | 0.13 | 0.11 | 0.10 | 0.09 | 0.09 |
1.6 | 0.64 | 0.64 | 0.64 | 0.51 | 0.32 | 0.26 | 0.20 | 0.18 | 0.17 | 0.16 |
1.4 | 0.73 | 0.73 | 0.73 | 0.73 | 0.60 | 0.46 | 0.33 | 0.28 | 0.26 | 0.24 |
1.2 | 0.74 | 0.74 | 0.74 | 0.74 | 0.73 | 0.67 | 0.50 | 0.41 | 0.37 | 0.35 |
1.0 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.65 | 0.57 | 0.50 | 0.47 |
0.8 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.70 | 0.64 | 0.60 |
0.6 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.73 |
Bold: C-IV-Long, Under line: C-III-Long, Italic: C-III-Short, Normal: C-III-Medium have the minimum target flaw sizes.
Bold: C-IV-Long, Italic: C-III-Short, Normal: C-I-Short have the minimum target flaw sizes.
Target flaw size, a/t, of axial flaw in static CASS pipe
Normalized flaw length ℓ/√(Rmt) | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
Stress ratio, Ph/Sm | 0.1 | 0.2 | 0.4 | 0.6 | 0.8 | 1.0 | 1.5 | 2.0 | 2.5 | 3.0 |
(a) Service condition: A&Ba | ||||||||||
1.6 | 0.21 | 0.04 | 0.03 | 0.03 | 0.03 | 0.03 | 0.03 | 0.02 | 0.02 | 0.02 |
1.4 | 0.44 | 0.44 | 0.20 | 0.12 | 0.10 | 0.09 | 0.08 | 0.07 | 0.07 | 0.07 |
1.2 | 0.65 | 0.65 | 0.65 | 0.54 | 0.33 | 0.26 | 0.20 | 0.17 | 0.16 | 0.15 |
1.0 | 0.74 | 0.74 | 0.74 | 0.74 | 0.67 | 0.55 | 0.38 | 0.32 | 0.29 | 0.27 |
0.8 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.73 | 0.59 | 0.50 | 0.44 | 0.41 |
0.6 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.70 | 0.64 | 0.59 |
0.4 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.74 |
0.2 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 |
(b) Service condition: C&Db | ||||||||||
2.4 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
2.2 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
2.0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
1.8 | 0.24 | 0.07 | 0.04 | 0.04 | 0.04 | 0.03 | 0.03 | 0.03 | 0.03 | 0.03 |
1.6 | 0.45 | 0.45 | 0.25 | 0.14 | 0.12 | 0.10 | 0.09 | 0.08 | 0.07 | 0.07 |
1.4 | 0.64 | 0.64 | 0.64 | 0.48 | 0.30 | 0.24 | 0.19 | 0.16 | 0.15 | 0.14 |
1.2 | 0.73 | 0.73 | 0.73 | 0.73 | 0.61 | 0.48 | 0.34 | 0.29 | 0.26 | 0.24 |
1.0 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.70 | 0.53 | 0.44 | 0.40 | 0.37 |
0.8 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.70 | 0.61 | 0.55 | 0.51 |
0.6 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.71 | 0.67 |
Normalized flaw length ℓ/√(Rmt) | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
Stress ratio, Ph/Sm | 0.1 | 0.2 | 0.4 | 0.6 | 0.8 | 1.0 | 1.5 | 2.0 | 2.5 | 3.0 |
(a) Service condition: A&Ba | ||||||||||
1.6 | 0.21 | 0.04 | 0.03 | 0.03 | 0.03 | 0.03 | 0.03 | 0.02 | 0.02 | 0.02 |
1.4 | 0.44 | 0.44 | 0.20 | 0.12 | 0.10 | 0.09 | 0.08 | 0.07 | 0.07 | 0.07 |
1.2 | 0.65 | 0.65 | 0.65 | 0.54 | 0.33 | 0.26 | 0.20 | 0.17 | 0.16 | 0.15 |
1.0 | 0.74 | 0.74 | 0.74 | 0.74 | 0.67 | 0.55 | 0.38 | 0.32 | 0.29 | 0.27 |
0.8 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.73 | 0.59 | 0.50 | 0.44 | 0.41 |
0.6 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.70 | 0.64 | 0.59 |
0.4 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.74 |
0.2 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 |
(b) Service condition: C&Db | ||||||||||
2.4 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
2.2 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
2.0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
1.8 | 0.24 | 0.07 | 0.04 | 0.04 | 0.04 | 0.03 | 0.03 | 0.03 | 0.03 | 0.03 |
1.6 | 0.45 | 0.45 | 0.25 | 0.14 | 0.12 | 0.10 | 0.09 | 0.08 | 0.07 | 0.07 |
1.4 | 0.64 | 0.64 | 0.64 | 0.48 | 0.30 | 0.24 | 0.19 | 0.16 | 0.15 | 0.14 |
1.2 | 0.73 | 0.73 | 0.73 | 0.73 | 0.61 | 0.48 | 0.34 | 0.29 | 0.26 | 0.24 |
1.0 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.70 | 0.53 | 0.44 | 0.40 | 0.37 |
0.8 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.70 | 0.61 | 0.55 | 0.51 |
0.6 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.75 | 0.74 | 0.71 | 0.67 |
In all condition, S-IV-Long have the minimum target flaw sizes.
Bold: S-IV-Long, Normal: S-I-Short have the minimum target flaw sizes.
Discussion
Identification of the Failure Materials.
As a result of the PFM analysis, the static CASS pipes, which had the minimum target flaw sizes, are summarized in Table 12. Among them, S-I-Short or S-II-Short were materials having the lower tensile properties, and S-IV-Short or S-IV-Long were materials having the lower fracture toughness. These results indicate that the dominant mechanical properties for failure were different in each condition.
Static CASS pipe materials mainly having the minimum target flaw sizes in each flaw direction and service condition in the result of PFM analysis
Service condition | |||
---|---|---|---|
Flaw direction | Stress level (value of the stress ratio) | A & B | C & D |
Circumferential | High (over 2.2) | — | S-IV-Long |
Intermediate (around 1.2–2.0) | S-II-Short | S-II-Short | |
Low (under 1.0) | S-IV-Long | S-IV-Long | |
Axial | High (over 2.0) | — | S-I-Short |
Intermediate (around 1.4–1.8) | S-IV-Long | S-IV-Long | |
Low (under 1.2) |
Service condition | |||
---|---|---|---|
Flaw direction | Stress level (value of the stress ratio) | A & B | C & D |
Circumferential | High (over 2.2) | — | S-IV-Long |
Intermediate (around 1.2–2.0) | S-II-Short | S-II-Short | |
Low (under 1.0) | S-IV-Long | S-IV-Long | |
Axial | High (over 2.0) | — | S-I-Short |
Intermediate (around 1.4–1.8) | S-IV-Long | S-IV-Long | |
Low (under 1.2) |
Then, simple PFM analyses were performed and the dominant mechanical properties for failure were confirmed. In order to identify the combination of mechanical properties failed in PFM samplings, failure or not were analyzed in the typical conditions of stress levels and flaw depths for materials which had the minimum target flaw sizes. Conditions of stress levels and flaw depths in these analyses are shown in Table 13. In these analyses, the sampling number is set to 50,000 times for the convenience of the analysis.
Probabilistic fracture mechanics analysis conditions for Identification of failure materials (Sampling number: 50,000 times)
(a) Circumferential flaw | ||||
---|---|---|---|---|
Case | (Pm + Pb + Pe)/Sm | Flaw angle 2θ (deg.) | Normalized flaw depth a/t | Analysis object of CASS pipe material |
High stress/small flaw | 2.4 | 30 | 0.026 | S-IV-Long |
Intermediate stress/medium flaw | 1.6 | 0.501 | S-II-Short |
(a) Circumferential flaw | ||||
---|---|---|---|---|
Case | (Pm + Pb + Pe)/Sm | Flaw angle 2θ (deg.) | Normalized flaw depth a/t | Analysis object of CASS pipe material |
High stress/small flaw | 2.4 | 30 | 0.026 | S-IV-Long |
Intermediate stress/medium flaw | 1.6 | 0.501 | S-II-Short |
(b) Axial flaw | ||||
---|---|---|---|---|
Case | Ph/Sm | Normalized flaw length ℓ/√(Rmt) | Normalized flaw depth a/t | Analysis object of CASS pipe material |
High stress/small flaw | 2.0 | 1.0 | 0.010 | S-I-Short |
Intermediate stress/medium flaw | 1.8 | 0.038 | S-IV-Long |
(b) Axial flaw | ||||
---|---|---|---|---|
Case | Ph/Sm | Normalized flaw length ℓ/√(Rmt) | Normalized flaw depth a/t | Analysis object of CASS pipe material |
High stress/small flaw | 2.0 | 1.0 | 0.010 | S-I-Short |
Intermediate stress/medium flaw | 1.8 | 0.038 | S-IV-Long |
Circumferential Flaw.
Figure 16 shows the combination of the sampled mechanical properties, σf and JIc, and those of failure. In the case of High stress/Small flaw as shown in Fig. 16(a), ductile fracture occurred when both σf and JIc were located near lower limits. On the other hand, in the case of Intermediate stress/Medium flaw as shown in Fig. 16(b), ductile fractures occurred when σf were located at the lower limit. These results indicated that dominant mechanical properties are different in each stress or flaw size. That is, it is supposed that only tensile properties are dominant at intermediate stress and medium flaw size, but both tensile properties and fracture toughness affect the ductile fracture at high stress and small flaw size.

The failed sampling results in circumferential flaw for static CASS: (a) high stress/small flaw (S-IV-long) and (b)intermediate stress/medium flaw (S-II-short)
Axial Flaw.
Figure 17 shows the combination of the sampled mechanical properties and those of failure. In the case of High stress/Small flaw as shown In Fig. 17(a), sampling results failed only when σf were located at the lower limit, and in this case, the dominant failure mode was plastic collapse. On the other hand, in cases of Intermediate stress/Medium flaw as shown in Fig. 17(b), the fracture toughnesses are located at the lower limit of failure sampling results. These results indicate that failure modes were different at higher stress or intermediate stress. That is, the fracture toughness was dominant in intermediate stress because dominant failure mode was ductile fracture, while only σf affects failure at high stress because plastic collapse could easily occur.

The failed sampling results in axial flaw for static CASS: (a) high stress/small flaw (S-I-short) and (b) intermediate stress/medium flaw (S-IV-long)
Comparison of the Target Flaw Sizes With ASME Code Case N-838.
In CC N-838 [4,5], the target flaw sizes were calculated by PFM analysis for primary coolant pipe made of CASS in the United States, and the target flaw size tables were arranged. These were compared with the tables analyzed in this study.
The comparison of the PFM analysis conditions between this study for static CASS and CC N -838 is shown in Table 14. There are many differences in the analysis conditions, for example, material conditions (chemical composition, casting process, thermal aging time, and temperature), fracture evaluation methods, and variations in mechanical properties and size of CASS pipes. The comparison of the PFM analysis results between this study for static CASS and CC N-838 is shown in Fig. 18 for circumferential flaws and Fig. 19 for axial flaws. Compared with CC N-838, the target flaw sizes in this study had the tendencies described below.

Comparison of the target flaw size for circumferential flaw between this study (static CASS) and ASME CC N-838: (a)failure probability: 10−6 (corresponding to service condition A& B in this study) and (b) failure probability: 10−4 (corresponding to service condition C & D in this study)

Comparison of the target flaw size for axial flaw in the case of 10−6 in failure probability between this study (static CASS) and ASME CC N-838
Comparison of PFM analysis condition between this study and CC N-838
PFM analysis condition | This study | CC N-838 [4,5] | |
---|---|---|---|
Size of pipe | Outer diameter (mm) | 882 | 840 |
Thickness (mm) | 72.7 | 60 | |
Condition of material | Material type | CF8M (JIS-SCS14A) | CF3, CF3A, CF3M, CF8, CF8A, CF8M |
Casting process | Separated static and centrifugal CASS | Mixed static and centrifugal CASS | |
Ferrite content (vol. %) | Static CASS: 9.44–20.54 Centrifugal CASS: 9.08–22.24 | 20–25 | |
Thermal aging temperature (°C) | 325 | 260–350 | |
Thermal aging time (h) | 20,000–525,600 | Fully aged | |
Advanced | Correlation of each material properties | Considered | — |
sampling method | Lower limit of tensile properties | Considered | — |
Evaluation method | Fracture evaluation method | Two-parameters method and limit load method | J-T method |
Calculation method for J-integral | Reference stress method | Simplified method uniquely developed in the United States | |
Allowable failure probability | Service condition A | 10−6 | 10−6 |
Service condition B | 10−6 | 10−5 | |
Service condition C & D | 10−4 | 10−4 |
PFM analysis condition | This study | CC N-838 [4,5] | |
---|---|---|---|
Size of pipe | Outer diameter (mm) | 882 | 840 |
Thickness (mm) | 72.7 | 60 | |
Condition of material | Material type | CF8M (JIS-SCS14A) | CF3, CF3A, CF3M, CF8, CF8A, CF8M |
Casting process | Separated static and centrifugal CASS | Mixed static and centrifugal CASS | |
Ferrite content (vol. %) | Static CASS: 9.44–20.54 Centrifugal CASS: 9.08–22.24 | 20–25 | |
Thermal aging temperature (°C) | 325 | 260–350 | |
Thermal aging time (h) | 20,000–525,600 | Fully aged | |
Advanced | Correlation of each material properties | Considered | — |
sampling method | Lower limit of tensile properties | Considered | — |
Evaluation method | Fracture evaluation method | Two-parameters method and limit load method | J-T method |
Calculation method for J-integral | Reference stress method | Simplified method uniquely developed in the United States | |
Allowable failure probability | Service condition A | 10−6 | 10−6 |
Service condition B | 10−6 | 10−5 | |
Service condition C & D | 10−4 | 10−4 |
Bold: Different condition with CC N-838.
For circumferential flaws: Fig. 18
In the case of 10−6 in failure probability, corresponding to service condition A & B in this study, this study had nearly same sizes in the case of 2θ = 30–40 deg. and smaller sizes in the case of 2θ = 60 deg.
In the case of 10−4 in failure probability, corresponding to service condition C & D in this study, this study had smaller sizes.
For axial flaws: Fig. 19
This study had larger sizes in the case of ℓ/√(Rmt) = 0.2–1.0 and nearly the same sizes in the case of ℓ/√(Rmt) = 2.0.
The target flaw sizes in this study tended to be smaller for circumferential flaws and larger for axial flaws. The difference in fracture evaluation methods is considered to be the main reason of this tendency. The methods of fracture evaluation and calculating J-integral applied in this study were relatively severe for circumferential flaw, compared with the J-T method in CC N-838. Therefore, especially in the case of 10−4 in failure probability, target flaw sizes for circumferential flaws were smaller in this study.
On the other hand, in the case of 10−6 in failure probability, the difference in target flaw sizes for circumferential flaws between this study and CC N-838 becomes smaller. This result is considered to be mainly caused by the variations in mechanical properties. Especially, the lower limit of the tensile properties was set in this study as described in the section “Advanced Sampling Methods for the Mechanical Properties,” so very small tensile property was not sampled in the PFM analysis, and then failure became less likely to happen and target flaw sizes became larger. The effect became more significant when increasing sampling number to calculate smaller failure probability.
Consequently, in the case of 10−6 in failure probability, the target flaw sizes for circumferential flaws in this study became relatively larger.
The size of the pipe in the PFM analysis also affects the target flaw size. The pipe thickness in this study was about 1.2 times thicker than CC N-838 and flaw depth, a, was 1.2 times larger for the same normalized target flaw sizes, a/t. Deep flaws have higher J-integral, which causes higher failure probability and smaller target flaw sizes.
Although the magnitude of target flaw sizes was different between this study and CC N-838, the trends in the relationship between the target flaw sizes and the stress ratios are similar despite many differences in the analysis conditions. It indicates that these PFM analysis methods, conditions and results are appropriate.
Conclusion
In order to obtain reasonable target flaw size tables for the PD examination system in Japan, fracture evaluation of flaws in centrifugal or static CASS pipes by PFM analysis code, PREFACE, was carried out. The main results are as shown below.
For PFM analysis, the variations in mechanical properties for the analysis (such as standard deviations, correlations, lower limits) were arranged based on prediction errors for actual test data. In addition, the TSS model and the modified TSS model were applied as the thermal aging prediction model in order to consider the differences in tensile properties for each casting process, centrifugal or static.
In order to estimate the conservative target flaw sizes for all PWR plants planned to be restarted in Japan, PFM analyses were performed for typical CASS pipes, and the target flaw sizes were calculated respectively for centrifugal CASS and static CASS. Based on the results of PFM analyses, the target flaw size tables for PD examination system were arranged by using the combination of the minimum target flaw sizes of each pipe. In almost all conditions, the target flaw sizes are smaller in static CASS than centrifugal CASS.
Depending on the flaw direction and stress ratio, the CASS pipe material, which is characterized by chemical composition, thermal aging time and casting process, having the minimum target flaw sizes were different. It was because the failure mode and dominant mechanical properties could be switched depending on the flaw direction and stress ratio. It indicates that PFM evaluations under the conditions of different thermal aging time are necessary to provide conservative failure probabilities and target flaw sizes.
This study and the ASME Code Case N-838 have similar trends in the relationship between the target flaw sizes and the stress ratios. The magnitude of the relation of the target flaw sizes between this study and CC N-838 are different in each flaw direction, flaw length, or allowable failure probability. It is affected by differences in the conditions of material (chemical composition, casting process, thermal aging time and temperature), fracture evaluation methods, the variations in mechanical properties, and size of CASS pipe.
Data Availability Statement
The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.
Nomenclature
- a =
flaw depth
- C1 =
constant of J-R curve after thermal aging
- C2 =
constant of J-R curve after thermal aging
- Do =
outer diameter
- E =
Young's modulus
- JIc =
fracture toughness which is JR value at ΔaIc on a J-R curve
- JIc,p =
JIc after thermal aging predicted by H3T model
- JIc,s =
sampled value of JIc by the Monte Carlo method
- JR =
fracture toughness on a J-R curve after thermal aging
- J6 =
JR value at Δa = 6 mm
- J6,p =
J6 after thermal aging predicted by H3T model
- J6,s =
sampled value of J6 by the Monte Carlo method
- l =
flaw length of an axial flaw
- N =
total sampling number
- Nf =
sampling number of failures
- Pb =
primary bending stress
- Pe =
secondary bending stress
- Ph =
hoop stress
- Pm =
membrane stress
- pf =
failure probability
- R =
correlation coefficient
- Ri =
inner radius (=Ro-t)
- Rm =
medium radius (=(Ro + Ri)/2)
- Ro =
outer radius (=Do/2)
- SJIc =
standard deviation of JIc
- SJ6 =
standard deviation of J6
- Sm =
design stress
- Sσf =
standard deviation of σf
- SσSS =
standard deviation of σSS
- Sσy =
standard deviation of σy
- SF =
safety factor
- SX, SY =
standard deviation of X or Y
- t =
thickness
- Tage =
thermal aging temperature
- γf =
increasing ratio of σf by thermal aging (=σf/σf0)
- γy =
increasing ratio of σy by thermal aging (=σy/σy0)
- Δa =
flaw extension on a J-R curve
- ΔaIc =
flaw extension at JR = JIc
- εSS =
strain on a stress-strain curve after thermal aging
- θ =
half flow angle of a circumferential flaw
- μX, μY =
average value of X or Y
- ν =
Poisson's ratio
- σf =
flow stress, which is average of tensile strength and yield stress
- σf,p =
σf after thermal aging predicted by TSS model or modified TSS model
- σf0 =
initial flow stress before thermal aging
- σf0,p =
σf0 predicted by TSS model or modified TSS model
- σf,s =
sampled value of σf by the Monte Carlo method
- σy =
yield stress
- σy,p =
σy after thermal aging predicted by TSS model or modified TSS model
- σy,s =
sampled value of σy by the Monte Carlo method
- σSS =
stress on a stress-strain curve after thermal aging
- σSS,avg =
average value of σSS predicted by Table 2
- σSS,s =
sampled value of σSS by the Monte Carlo method
- σy0 =
initial yield stress before thermal aging
- σy0,p =
σy0 predicted by TSS model or modified TSS model
- CASS =
cast austenitic stainless steel
- FFS =
fitness-for-service
- H3T model =
fracture toughness prediction method
- JSME =
The Japan Society of Mechanical Engineers
- PD =
performance demonstration
- PFM =
probability fracture mechanics
- PWR =
pressurized water reactor
- TSS model =
true stress-true strain curve prediction method
- UT =
ultrasonic testing
Appendix
The details of methods of fracture analysis based on the JSME rules on FFS are shown below.
Two-Parameter Method.
Ductile fracture is evaluated by Two-parameter method. This method uses two kinds of parameters Kr and Sr relating to fracture mechanics and mechanics of materials, and is fundamentally equivalent to the J-T criterion for evaluation of ductile instability. Fracture assessment curve and ductile crack growth line are calculated, and fracture will occur when these don't intersect. Fracture assessment curve is defined by Sr and Kr and is calculated by increasing stresses for the initial flaw size. Ductile crack growth line is defined by Sr′ and Kr′ and is calculated by increasing the flaw depth for constant stresses.
Sr and Sr′ are defined as following equations.
εSS,ref is a strain value on the S-S curve corresponding to σSS,ref. K is a stress intensity factor calculated by an equation for a semi-elliptical inner surface crack in a cylinder.
Limit-Load Method.
Plastic collapse is evaluated by Limit-load method. Fracture will occur when the applied stress, Pbapp for circumferential flaws or Phapp for axial flaws, exceeds plastic collapse stress, or .