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 [13]. 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 [68], 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.

Fig. 1
Simulated flaw shapes: (a) circumferential flaw and (b)axial flaw
Fig. 1
Simulated flaw shapes: (a) circumferential flaw and (b)axial flaw
Close modal

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 [68], 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).

Fig. 2
Probabilistic fracture mechanics analysis flow by PREFACE
Fig. 2
Probabilistic fracture mechanics analysis flow by PREFACE
Close modal
Fig. 3
Relationship between predicted flow stress, σf,p, and prediction error, σf-σf,p for centrifugal CASS
Fig. 3
Relationship between predicted flow stress, σf,p, and prediction error, σf-σf,p for centrifugal CASS
Close modal

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.

σy and σf were set as random variables in accordance with normal distribution. Their mean values were calculated using the TSS model for centrifugal CASS and the modified TSS model [11] for static CASS. The standard deviation was provided for each centrifugal CASS and static CASS based on the residual analysis of thermal aging CASS database. Figure 3 shows the relationship between the predicted σf and the prediction error, σf-σf,p, by the original TSS model for the centrifugal CASS. In this figure, the predicted σf are segmented for every 20 MPa and the standard deviation of each region is shown by the error bar. The variations in prediction errors depend on the σf levels. It seems that the larger the σf becomes, the larger the standard deviation becomes. From this relationship, the standard deviation was recalculated using parameters, which are the prediction errors normalized by the predicted σy,σy,p or the predicted σf, σf,p. Figure 4 shows the relationship between σf,p and the normalized errors, (σf-σf,p)/σf,p, for centrifugal CASS. Figure 5 shows the σf dependence on the standard deviations and the normalized standard deviations (Coefficient of variations). Coefficients of variation have less σf dependence and a lower correlation coefficient. Then, the probability distributions of σy and σf are defined by the following equations:
(1)
(2)
Fig. 4
Relationship between predicted flow stress, σf,p, and normalized prediction error, (σf-σf,p)/σf,p for centrifugal CASS
Fig. 4
Relationship between predicted flow stress, σf,p, and normalized prediction error, (σf-σf,p)/σf,p for centrifugal CASS
Close modal
Fig. 5
Relationship between predicted flow stress, σf,p, and standard deviations of prediction error for centrifugal CASS: (a)prediction error and (b) normalized prediction error
Fig. 5
Relationship between predicted flow stress, σf,p, and standard deviations of prediction error for centrifugal CASS: (a)prediction error and (b) normalized prediction error
Close modal
Fig. 6
Correlation of the prediction error of JIc and J6
Fig. 6
Correlation of the prediction error of JIc and J6
Close modal

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.

Table 1

Probabilistic distributions and standard deviations for random variables

Standard deviation, Sσy, Sσf, SσSS, SJIc, and SJ6
Random variableProbabilistic distributionCentrifugal CASSStatic CASS
σy/σy,pNormal0.06650.1304
σf/σf,pNormal0.03510.0815
σSSNormalRefer to Table 2 
JIcLognormal0.2396
J6Lognormal0.1781
Standard deviation, Sσy, Sσf, SσSS, SJIc, and SJ6
Random variableProbabilistic distributionCentrifugal CASSStatic CASS
σy/σy,pNormal0.06650.1304
σf/σf,pNormal0.03510.0815
σSSNormalRefer to Table 2 
JIcLognormal0.2396
J6Lognormal0.1781
The Stress–Strain curves (S–S curves) were determined in accordance with the Japan Society of Mechanical Engineers (JSME) rules on fitness-for-service (FFS). Their mean stress, σSS, at each strain, εSS, and its standard deviation, SσSS, are shown in Table 2 corresponding to sampled flow stress, σf,s, by Monte Carlo simulation following Eq. (2). σSS were determined as random variables in accordance with normal distribution. Their standard deviations were defined based on the standard deviations of the prediction errors in the tensile tests of CASS [8,13]. The probability distributions of σSS are defined by the following equations:
(3)
Table 2

Date table for sampling of S-S curve

No.σf,s (MPa)Strain point1234567
1281εSS (%)0.20.51361220
σSS (MPa)134166185227274366473
SσSS (MPa)10.812.512.115.919.926.634.0
2318εSS (%)0.20.5124815
σSS (MPa)149192222254297371484
SσSS (MPa)13.613.113.315.416.618.922.9
3377εSS (%)0.20.30.512410
σSS (MPa)178204235269308360487
SσSS (MPa)13.914.613.817.121.125.736.8
4423εSS (%)0.20.30.40.5125
σSS (MPa)182211235255314365462
SσSS (MPa)14.315.717.818.124.330.539.0
No.σf,s (MPa)Strain point1234567
1281εSS (%)0.20.51361220
σSS (MPa)134166185227274366473
SσSS (MPa)10.812.512.115.919.926.634.0
2318εSS (%)0.20.5124815
σSS (MPa)149192222254297371484
SσSS (MPa)13.613.113.315.416.618.922.9
3377εSS (%)0.20.30.512410
σSS (MPa)178204235269308360487
SσSS (MPa)13.914.613.817.121.125.736.8
4423εSS (%)0.20.30.40.5125
σSS (MPa)182211235255314365462
SσSS (MPa)14.315.717.818.124.330.539.0

Fracture Toughness.

The J-R curves were defined by the following equations:
(4)
(5)
JIc and J6 are the fracture toughness, which are the JR value when the flaw extension, Δa, becomes ΔaIc and 6 mm on J-R curves. C1 and C2 are constant and calculated so as to satisfy (JR, Δa) = (JIc,s, ΔaIc) and (J6,s, 6 mm). JIc and J6 are defined as random variables in accordance with a lognormal distribution. Their mean values, JIc,p and J6,p, are predicted using H3T model. Their standard deviations, SJIc and SJ6, are defined as the standard deviation of the logarithmic prediction errors for the test database of CASS, which were used in the development of the H3T model [10]. The probability density functions of JIc and J6 are expressed by the following equations:
(6)
(7)

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.

A set of the mechanical properties (σy,s, σf,s, σSS,s, JIc,s, J6,s) was sampled N times by Monte Carlo simulation. For each parameter set, condition of load and flaw size, fracture evaluations were performed. When the number of failures, Nf, is calculated, the failure probability pf is defined by the following equation:
(8)

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.

Fig. 7
Correlation of the prediction error of σf and JIc (centrifugal CASS)
Fig. 7
Correlation of the prediction error of σf and JIc (centrifugal CASS)
Close modal
Fig. 8
Comparison of example of sampling results and prediction error data of fracture toughness
Fig. 8
Comparison of example of sampling results and prediction error data of fracture toughness
Close modal
Table 3

Correlation coefficient of the prediction error in the material properties

Combination of material propertiesCasting processCorrelation coefficient
(a) σy-σf, σf-σSS and JIc-J6
Tensile propertiesσy/σy,p and σf/σf,pCentrifugal0.61
Static0.78
σf/σy,p and σSS/σSS,avgCommon0.70
Fracture toughnesslog(JIc) and log(J6)Common0.78
Combination of material propertiesCasting processCorrelation coefficient
(a) σy-σf, σf-σSS and JIc-J6
Tensile propertiesσy/σy,p and σf/σf,pCentrifugal0.61
Static0.78
σf/σy,p and σSS/σSS,avgCommon0.70
Fracture toughnesslog(JIc) and log(J6)Common0.78
(b) Tensile property-fracture toughness
σy and JIcCentrifugal−0.16
σy and J6−0.14
σf and JIc−0.45
σf and J6−0.33
σy and JIcStatic−0.11
σy and J60.05
σf and JIc−0.09
σf and J60.08
(b) Tensile property-fracture toughness
σy and JIcCentrifugal−0.16
σy and J6−0.14
σf and JIc−0.45
σf and J6−0.33
σy and JIcStatic−0.11
σy and J60.05
σf and JIc−0.09
σf and J60.08
Based on these results, the correlations between the mechanical properties were considered in PFM analysis when the absolute value of R was larger than 0.40. When two random variables, X and Y, have a correlation coefficient R and a value of X is sampled in PFM, the mean value of Y, E(Y|X), and its standard deviation, E(Var(Y|X)), are expressed by the following equations:
(9)
(10)

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.

Fig. 9
Cumulative probability distributions of yield strength before/after thermal aging under the representative material of static CASS and thermal aging condition
Fig. 9
Cumulative probability distributions of yield strength before/after thermal aging under the representative material of static CASS and thermal aging condition
Close modal

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.

Fig. 10
Predicted average values of the material properties of CASS pipe of Japanese PWR plants: (a) centrifugal and (b) static
Fig. 10
Predicted average values of the material properties of CASS pipe of Japanese PWR plants: (a) centrifugal and (b) static
Close modal
Table 4

Standard deviations of normalized prediction error of tensile properties before and after thermal aging

Casting processParameterThermal agingRandom variableStandard deviation
CentrifugalσyBefore agingσy0/σy0,p0.1636
After agingσy/σy,p0.0665
σfBefore agingσf0/σf0,p0.0654
After agingσf/σf,p0.0351
StaticσyBefore agingσy0/σy0,p0.0726
After agingσy/σy,p0.1304
σfBefore agingσf0/σf0,p0.0398
After agingσf/σf,p0.0815
Casting processParameterThermal agingRandom variableStandard deviation
CentrifugalσyBefore agingσy0/σy0,p0.1636
After agingσy/σy,p0.0665
σfBefore agingσf0/σf0,p0.0654
After agingσf/σf,p0.0351
StaticσyBefore agingσy0/σy0,p0.0726
After agingσy/σy,p0.1304
σfBefore agingσf0/σf0,p0.0398
After agingσf/σf,p0.0815

Recalculation of J6,s When JIc,s Was Higher Than J6,s.

When JIc,s became higher than J6,s in a sampling result in PFM, the J-R curve could not be defined by the Eq. (4). In the previous study [7,8], in such cases, both JIc,s and J6,s were resampled. As a result, the probability of sampling lower JIc,s became higher than the originally assumed probability, and failure probabilities could be estimated higher due to the lower JIc,s. Therefore, in order not to change the probability distributions of JIc,s from actual distributions, only J6,s is recalculated by the following equation in the case of JIc,s > J6,s to set J6,s slightly larger than JIc,s
(11)

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.

Table 5

Selection of analysis conditions for Japanese CASS pipe in PFM analysisa,b,c,d

Thermal aging times at 325 °C (h)
Casting processCASS pipe materialFerrite contenta (vol. %)20,000b25,000c50,000c525,600dNotes
CentrifugalC-I9.08C-I-ShortC-I-MediumC-I-LongC-I-Short has the lowest tensile properties.
C-II9.90C-II-ShortC-II- MediumC-II-Long
C-III16.38C-III-ShortC-III- MediumC-III-Long
C-IV22.24C-IV-ShortC-IV- MediumC-IV-LongC-IV has highest ferrite content in Japanese PWR plants.
C-IV-Long has the lowest fracture toughness.
StaticS-I9.44S-I-ShortS-I- MediumS-I-LongS-I-Short has the lowest tensile properties.
S-II10.18S-II-ShortS-II-MediumS-II-Long
S-III15.85S-III-ShortS-III-MediumS-III-Long
S-IV20.54S-IV-ShortS-IV-MediumS-IV-LongS-IV-Long has the lowest fracture toughness.
Thermal aging times at 325 °C (h)
Casting processCASS pipe materialFerrite contenta (vol. %)20,000b25,000c50,000c525,600dNotes
CentrifugalC-I9.08C-I-ShortC-I-MediumC-I-LongC-I-Short has the lowest tensile properties.
C-II9.90C-II-ShortC-II- MediumC-II-Long
C-III16.38C-III-ShortC-III- MediumC-III-Long
C-IV22.24C-IV-ShortC-IV- MediumC-IV-LongC-IV has highest ferrite content in Japanese PWR plants.
C-IV-Long has the lowest fracture toughness.
StaticS-I9.44S-I-ShortS-I- MediumS-I-LongS-I-Short has the lowest tensile properties.
S-II10.18S-II-ShortS-II-MediumS-II-Long
S-III15.85S-III-ShortS-III-MediumS-III-Long
S-IV20.54S-IV-ShortS-IV-MediumS-IV-LongS-IV-Long has the lowest fracture toughness.
a

Predicted by ASTM A800 [12].

b

Corresponding to the shortest operating time in Japanese PWR plants.

c

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.

d

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.

Table 6

Probabilistic fracture mechanics analysis conditions of loads and flaw lengths

Flaw direction
ItemParameterCircumferentialAxial
LoadDefinition of load parameter(Pm + Pb + Pe)/SmPh/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 valueSm = 116 MPaSm = 116 MPa
Pm = 58 MPa
Pe = 0 MPa
Flaw lengthDefinition of flaw length parameterFlaw angle, 2θ (deg.)Normalized flaw length, l/√ (Rmt)
Analysis condition (range)2θ = 10–60l/√ (Rmt)=0.1–3.0
Flaw direction
ItemParameterCircumferentialAxial
LoadDefinition of load parameter(Pm + Pb + Pe)/SmPh/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 valueSm = 116 MPaSm = 116 MPa
Pm = 58 MPa
Pe = 0 MPa
Flaw lengthDefinition of flaw length parameterFlaw angle, 2θ (deg.)Normalized flaw length, l/√ (Rmt)
Analysis condition (range)2θ = 10–60l/√ (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.

Table 7

Probabilistic fracture mechanics analysis conditions of allowable failure probability and sampling number

Service condition
ParameterA & BC & D
Allowable failure probability<1 × 106<1 × 104
Sampling number5 × 1075 × 105
Service condition
ParameterA & BC & D
Allowable failure probability<1 × 106<1 × 104
Sampling number5 × 1075 × 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:

Fig. 11
The target flaw sizes for circumferential flaws in each stress level (service condition: C & D, flaw angle: 2θ = 30 deg., material: S-IV)
Fig. 11
The target flaw sizes for circumferential flaws in each stress level (service condition: C & D, flaw angle: 2θ = 30 deg., material: S-IV)
Close modal
  • 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.

Fig. 12
Comparison of the target flaw sizes for circumferential flaws in each static CASS pipe material (service condition: C&D, flaw angle: 2θ = 30 deg)
Fig. 12
Comparison of the target flaw sizes for circumferential flaws in each static CASS pipe material (service condition: C&D, flaw angle: 2θ = 30 deg)
Close modal

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:

Fig. 13
The target flaw sizes for axial flaws in each stress ratio (service condition: C&D, nondimensional flaw length: ℓ/√(Rmt) = 1.0, material: S-IV)
Fig. 13
The target flaw sizes for axial flaws in each stress ratio (service condition: C&D, nondimensional flaw length: ℓ/√(Rmt) = 1.0, material: S-IV)
Close modal
  • 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.

Fig. 14
Comparison of the target flaw sizes for axial flaws in each static CASS pipe material (service condition: C&D, flaw length: ℓ/√(Rmt) = 1.0)
Fig. 14
Comparison of the target flaw sizes for axial flaws in each static CASS pipe material (service condition: C&D, flaw length: ℓ/√(Rmt) = 1.0)
Close modal

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.

Fig. 15
Example of the determination of the target flaw size based on the PFM results
Fig. 15
Example of the determination of the target flaw size based on the PFM results
Close modal
Table 8

Target flaw size, a/t, of circumferential flaw in centrifugal CASS pipe

Flaw angle, 2θ (deg.)
Stress ratio, (Pm + Pb + Pe)/Sm102030405060
(a) Allowable condition: A&Ba
1.60.740.730.680.510.430.38
1.40.740.740.720.710.620.54
1.20.740.740.740.720.710.70
1.00.750.740.740.740.730.72
0.80.750.750.740.740.740.74
0.60.750.750.750.750.740.74
(b) Allowable condition: C&Db
2.60.080.060.050.050.040.04
2.40.350.120.100.090.090.09
2.20.540.260.190.160.150.14
2.00.710.610.350.280.240.22
1.80.740.720.620.450.380.34
1.60.740.740.720.680.550.48
1.40.750.740.740.730.710.65
1.20.750.750.740.740.730.73
1.00.750.750.750.740.740.74
Flaw angle, 2θ (deg.)
Stress ratio, (Pm + Pb + Pe)/Sm102030405060
(a) Allowable condition: A&Ba
1.60.740.730.680.510.430.38
1.40.740.740.720.710.620.54
1.20.740.740.740.720.710.70
1.00.750.740.740.740.730.72
0.80.750.750.740.740.740.74
0.60.750.750.750.750.740.74
(b) Allowable condition: C&Db
2.60.080.060.050.050.040.04
2.40.350.120.100.090.090.09
2.20.540.260.190.160.150.14
2.00.710.610.350.280.240.22
1.80.740.720.620.450.380.34
1.60.740.740.720.680.550.48
1.40.750.740.740.730.710.65
1.20.750.750.740.740.730.73
1.00.750.750.750.740.740.74
a

Bold: C-IV-Long, Normal: C-III-Medium have the minimum target flaws sizes.

b

Bold: C-IV-Long, Under line: C-III-Long, Italic: C-III-Short, have the minimum target flaws sizes.

Table 9

Target flaw size, a/t, of circumferential flaw in static CASS pipe

Flaw angle, 2θ (deg.)
Stress ratio, (Pm + Pb + Pe)/Sm102030405060
(a) Service condition: A&Ba
1.60.580.330.200.170.150.14
1.40.710.680.440.330.270.24
1.20.740.740.700.570.460.40
1.00.750.750.740.730.680.58
0.80.750.750.750.740.740.74
0.60.750.750.750.750.750.74
(b) Service condition: C&Db
2.60.010.010.010.010.010.01
2.40.030.030.020.020.020.02
2.20.200.070.060.060.060.05
2.00.450.180.130.120.100.10
1.80.640.440.250.200.170.16
1.60.720.690.500.350.290.26
1.40.740.740.700.570.460.40
1.20.750.750.740.730.670.57
1.00.750.750.750.740.740.73
Flaw angle, 2θ (deg.)
Stress ratio, (Pm + Pb + Pe)/Sm102030405060
(a) Service condition: A&Ba
1.60.580.330.200.170.150.14
1.40.710.680.440.330.270.24
1.20.740.740.700.570.460.40
1.00.750.750.740.730.680.58
0.80.750.750.750.740.740.74
0.60.750.750.750.750.750.74
(b) Service condition: C&Db
2.60.010.010.010.010.010.01
2.40.030.030.020.020.020.02
2.20.200.070.060.060.060.05
2.00.450.180.130.120.100.10
1.80.640.440.250.200.170.16
1.60.720.690.500.350.290.26
1.40.740.740.700.570.460.40
1.20.750.750.740.730.670.57
1.00.750.750.750.740.740.73
a

Bold: S-IV-Long, Under line: S-II-Short, Normal: S-I-Short have the minimum target flaw sizes.

b

Bold: S-IV-Long, Under line: S-II-Short, Normal: S-IV-Short have the minimum target flaw sizes.

Table 10

Target flaw size, a/t, of axial flaw in centrifugal CASS pipe

Normalized flaw length ℓ/√(Rmt)
Stress ratio, Ph/Sm0.10.20.40.60.81.01.52.02.53.0
(a) Service condition: A&Ba
1.60.550.550.550.260.190.170.140.130.120.12
1.40.700.700.700.690.440.330.240.210.200.19
1.20.740.740.740.740.700.590.390.330.300.28
1.00.750.750.750.750.740.720.570.470.420.40
0.80.750.750.750.750.750.740.720.640.570.53
0.60.750.750.750.750.750.750.740.740.730.69
0.40.750.750.750.750.750.750.750.750.750.74
0.20.750.750.750.750.750.750.750.750.750.75
Normalized flaw length ℓ/√(Rmt)
Stress ratio, Ph/Sm0.10.20.40.60.81.01.52.02.53.0
(a) Service condition: A&Ba
1.60.550.550.550.260.190.170.140.130.120.12
1.40.700.700.700.690.440.330.240.210.200.19
1.20.740.740.740.740.700.590.390.330.300.28
1.00.750.750.750.750.740.720.570.470.420.40
0.80.750.750.750.750.750.740.720.640.570.53
0.60.750.750.750.750.750.750.740.740.730.69
0.40.750.750.750.750.750.750.750.750.750.74
0.20.750.750.750.750.750.750.750.750.750.75
(b) Service condition: C&Db
2.40000000000
2.20.160.030.030.030.020.020.020.020.020.02
2.00.300.300.080.070.060.060.050.040.040.04
1.80.480.480.480.180.140.130.110.100.090.09
1.60.640.640.640.510.320.260.200.180.170.16
1.40.730.730.730.730.600.460.330.280.260.24
1.20.740.740.740.740.730.670.500.410.370.35
1.00.750.750.750.750.750.740.650.570.500.47
0.80.750.750.750.750.750.750.740.700.640.60
0.60.750.750.750.750.750.750.750.750.740.73
(b) Service condition: C&Db
2.40000000000
2.20.160.030.030.030.020.020.020.020.020.02
2.00.300.300.080.070.060.060.050.040.040.04
1.80.480.480.480.180.140.130.110.100.090.09
1.60.640.640.640.510.320.260.200.180.170.16
1.40.730.730.730.730.600.460.330.280.260.24
1.20.740.740.740.740.730.670.500.410.370.35
1.00.750.750.750.750.750.740.650.570.500.47
0.80.750.750.750.750.750.750.740.700.640.60
0.60.750.750.750.750.750.750.750.750.740.73
a

Bold: C-IV-Long, Under line: C-III-Long, Italic: C-III-Short, Normal: C-III-Medium have the minimum target flaw sizes.

b

Bold: C-IV-Long, Italic: C-III-Short, Normal: C-I-Short have the minimum target flaw sizes.

Table 11

Target flaw size, a/t, of axial flaw in static CASS pipe

Normalized flaw length ℓ/√(Rmt)
Stress ratio, Ph/Sm0.10.20.40.60.81.01.52.02.53.0
(a) Service condition: A&Ba
1.60.210.040.030.030.030.030.030.020.020.02
1.40.440.440.200.120.100.090.080.070.070.07
1.20.650.650.650.540.330.260.200.170.160.15
1.00.740.740.740.740.670.550.380.320.290.27
0.80.750.750.750.750.740.730.590.500.440.41
0.60.750.750.750.750.750.750.740.700.640.59
0.40.750.750.750.750.750.750.750.750.740.74
0.20.750.750.750.750.750.750.750.750.750.75
(b) Service condition: C&Db
2.40000000000
2.20000000000
2.00000000000
1.80.240.070.040.040.040.030.030.030.030.03
1.60.450.450.250.140.120.100.090.080.070.07
1.40.640.640.640.480.300.240.190.160.150.14
1.20.730.730.730.730.610.480.340.290.260.24
1.00.750.750.750.750.740.700.530.440.400.37
0.80.750.750.750.750.750.740.700.610.550.51
0.60.750.750.750.750.750.750.750.740.710.67
Normalized flaw length ℓ/√(Rmt)
Stress ratio, Ph/Sm0.10.20.40.60.81.01.52.02.53.0
(a) Service condition: A&Ba
1.60.210.040.030.030.030.030.030.020.020.02
1.40.440.440.200.120.100.090.080.070.070.07
1.20.650.650.650.540.330.260.200.170.160.15
1.00.740.740.740.740.670.550.380.320.290.27
0.80.750.750.750.750.740.730.590.500.440.41
0.60.750.750.750.750.750.750.740.700.640.59
0.40.750.750.750.750.750.750.750.750.740.74
0.20.750.750.750.750.750.750.750.750.750.75
(b) Service condition: C&Db
2.40000000000
2.20000000000
2.00000000000
1.80.240.070.040.040.040.030.030.030.030.03
1.60.450.450.250.140.120.100.090.080.070.07
1.40.640.640.640.480.300.240.190.160.150.14
1.20.730.730.730.730.610.480.340.290.260.24
1.00.750.750.750.750.740.700.530.440.400.37
0.80.750.750.750.750.750.740.700.610.550.51
0.60.750.750.750.750.750.750.750.740.710.67
a

In all condition, S-IV-Long have the minimum target flaw sizes.

b

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.

Table 12

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 directionStress level (value of the stress ratio)A & BC & D
CircumferentialHigh (over 2.2)S-IV-Long
Intermediate (around 1.2–2.0)S-II-ShortS-II-Short
Low (under 1.0)S-IV-LongS-IV-Long
AxialHigh (over 2.0)S-I-Short
Intermediate (around 1.4–1.8)S-IV-LongS-IV-Long
Low (under 1.2)
Service condition
Flaw directionStress level (value of the stress ratio)A & BC & D
CircumferentialHigh (over 2.2)S-IV-Long
Intermediate (around 1.2–2.0)S-II-ShortS-II-Short
Low (under 1.0)S-IV-LongS-IV-Long
AxialHigh (over 2.0)S-I-Short
Intermediate (around 1.4–1.8)S-IV-LongS-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.

Table 13

Probabilistic fracture mechanics analysis conditions for Identification of failure materials (Sampling number: 50,000 times)

(a) Circumferential flaw
Case(Pm + Pb + Pe)/SmFlaw angle 2θ (deg.)Normalized flaw depth a/tAnalysis object of CASS pipe material
High stress/small flaw2.4300.026S-IV-Long
Intermediate stress/medium flaw1.60.501S-II-Short
(a) Circumferential flaw
Case(Pm + Pb + Pe)/SmFlaw angle 2θ (deg.)Normalized flaw depth a/tAnalysis object of CASS pipe material
High stress/small flaw2.4300.026S-IV-Long
Intermediate stress/medium flaw1.60.501S-II-Short
(b) Axial flaw
CasePh/SmNormalized flaw length ℓ/√(Rmt)Normalized flaw depth a/tAnalysis object of CASS pipe material
High stress/small flaw2.01.00.010S-I-Short
Intermediate stress/medium flaw1.80.038S-IV-Long
(b) Axial flaw
CasePh/SmNormalized flaw length ℓ/√(Rmt)Normalized flaw depth a/tAnalysis object of CASS pipe material
High stress/small flaw2.01.00.010S-I-Short
Intermediate stress/medium flaw1.80.038S-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.

Fig. 16
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)
Fig. 16
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)
Close modal

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.

Fig. 17
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)
Fig. 17
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)
Close modal

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.

Fig. 18
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)
Fig. 18
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)
Close modal
Fig. 19
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
Fig. 19
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
Close modal
Table 14

Comparison of PFM analysis condition between this study and CC N-838

PFM analysis conditionThis studyCC N-838 [4,5]
Size of pipeOuter diameter (mm)882840
Thickness (mm)72.760
Condition of materialMaterial typeCF8M (JIS-SCS14A)CF3, CF3A, CF3M, CF8, CF8A, CF8M
Casting processSeparated static and centrifugal CASSMixed static and centrifugal CASS
Ferrite content (vol. %)Static CASS: 9.44–20.54 Centrifugal CASS: 9.08–22.2420–25
Thermal aging temperature (°C)325260–350
Thermal aging time (h)20,000–525,600Fully aged
Advanced Correlation of each material propertiesConsidered
sampling methodLower limit of tensile propertiesConsidered
Evaluation methodFracture evaluation methodTwo-parameters method and limit load methodJ-T method
Calculation method for J-integralReference stress methodSimplified method uniquely developed in the United States
Allowable failure probabilityService condition A10−610−6
Service condition B10610−5
Service condition C & D10−410−4
PFM analysis conditionThis studyCC N-838 [4,5]
Size of pipeOuter diameter (mm)882840
Thickness (mm)72.760
Condition of materialMaterial typeCF8M (JIS-SCS14A)CF3, CF3A, CF3M, CF8, CF8A, CF8M
Casting processSeparated static and centrifugal CASSMixed static and centrifugal CASS
Ferrite content (vol. %)Static CASS: 9.44–20.54 Centrifugal CASS: 9.08–22.2420–25
Thermal aging temperature (°C)325260–350
Thermal aging time (h)20,000–525,600Fully aged
Advanced Correlation of each material propertiesConsidered
sampling methodLower limit of tensile propertiesConsidered
Evaluation methodFracture evaluation methodTwo-parameters method and limit load methodJ-T method
Calculation method for J-integralReference stress methodSimplified method uniquely developed in the United States
Allowable failure probabilityService condition A10−610−6
Service condition B10610−5
Service condition C & D10−410−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.

  1. 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.

  2. 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.

  3. 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.

  4. 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.

(Circumferential flaws)
(A1)
(A2)
(A3)
(A4)
(Axial flaws)
(A5)
(A6)
(A7)
Kr and Kr′ are defined as following equation:
(A8)
(A9)
In this study, Japp is calculated by reference stress method as shown below:
(A10)
(A11)
(A12)

ε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, Pb0 or Ph0.

(Circumferential flaw)
(A13)
(A14)
(A15)
(Axial flaw)
(A16)
(A17)
(A18)

References

1.
Chopra
,
O. K.
, and
Chung
,
H. M.
,
1985
, “
Aging of Cast Duplex Stainless Steels in LWR Systems
,”
Nucl. Eng. Des.
,
89
(
2–3
), pp.
305
318
.10.1016/0029-5493(85)90069-X
2.
Chopra
,
O. K.
,
1993
, “
Long-Term Embrittlement of Cast Duplex Stainless Steels in LWR Systems
,” Argonne National Laboratory, Argonne, IL, Report No.
NUREG/CR-4744, 7(2) ANL-93/11
.https://inis.iaea.org/collection/NCLCollectionStore/_Public/25/010/25010513.pdf
3.
Chung
,
H. M.
,
1992
, “
Aging and Life Prediction of Cast Duplex Stainless Steel Components
,”
Int. J. Pressure Vessels Piping
,
50
(
1–3
), pp.
179
213
.10.1016/0308-0161(92)90037-G
4.
American Society of Mechanical Engineers
,
2015
, “
Flaw Tolerance Evaluation of Cast Austenitic Stainless Steel Pipes Components
,”
ASME
Paper No. PVT-20-1162.10.1115/PVT-20-1162
5.
Griesbach
,
T. J.
,
Harris
,
D. O.
,
Cofie
,
N. G.
,
Chockie
,
A.
, and
Dedhia
,
D.
,
2013
, “
Basis for ASME Section XI Code Case for Flaw Tolerance Evaluation of CASS Pipes
,”
ASME
Paper No. PVP2013-97712.10.1115/PVP2013-97712
6.
Hojo
,
K.
,
Hayashi
,
S.
,
Nishi
,
W.
,
Kamaya
,
M.
,
Katsuyama
,
J.
,
Masaki
,
K.
,
Nagai
,
M.
,
Okamoto
,
T.
,
Takada
,
Y.
, and
Yoshimura
,
S.
,
2016
, “
Benchmark Analysis of Probabilistic Fracture Mechanics for Cast Stainless Steel Pipe
,”
Mech. Eng. J.
,
3
(
4
), p. 16–00083.10.1299/mej.16-00083
7.
Hayashi
,
S.
,
Ochi
,
M.
,
Hojo
,
K.
,
Yamane
,
T.
, and
Nishi
,
W.
,
2017
, “
Allowable Flaw Size of Japanese Cast Stainless Steel Pipe Using Probabilistic Fracture Mechanics Method
,”
ASME
Paper No. PVP2017-65633.10.1115/PVP2017-65633
8.
Nishi
,
W.
,
Hirota
,
T.
,
Ochi
,
M.
,
Takagoshi
,
D.
, and
Hojo
,
K.
,
2018
, “
Improvement of Target Flaw Sizes of CASS Pipe for PD Approval Using PFM Code PREFACE
,”
ASME
Paper No. PVP2018-85015.10.1115/PVP2018-85015
9.
Chopra
,
O. K.
,
1991
, “
Estimation of Fracture Toughness of Cast Stainless Steels During Thermal Aging in LWR Systems
,” Argonne National Laboratory, Argonne, IL, Report No. NUREG/CR-4513, ANL-90/42.
10.
Kawaguchi
,
S.
,
Nagasaki
,
T.
, and
Koyama
,
K.
,
2005
, “
Prediction Method of Tensile Properties and Fracture Toughness of Thermally Aged Cast Duplex Stainless Steel Pipes
,”
ASME
Paper No. PVP2005-71528.10.1115/PVP2005-71528
11.
Momotani
,
Y.
,
Takagoshi
,
D.
,
Nomura
,
Y.
,
Nishi
,
W.
,
Hirota
,
T.
, and
Hojo
,
K.
, “
Modified Prediction Model of Tensile Properties for Thermally Aged Cast Austenitic Stainless Steel of Static Cast
,”
ASME J. Pressure Vessel Technol.
, epub.
12.
ASTM A 800/A 800M
,
1989
,
Standard Practice for Steel Casting, Austenitic Alloy, Estimating Ferrite Content Thereof
,
American Society of Testing and Materials
,
West Conshohocken, PA
.
13.
Hojo
,
K.
,
Ochi
,
M.
,
Ogawa
,
N.
,
Hayashi
,
S.
,
Nishi
,
W.
, and
Yamane
,
T.
,
2017
, “
Flaw Evaluation of Cast Stainless Steel Pipe for Revision of JSME Rules on Fitness for Service
,”
Mech. Eng. J.
,
4
(
4
), p. 16-00642.10.1299/mej.16-00642
14.
Ainsworth
,
R. A.
,
1984
, “
The Assessment of Defects in Structures of Strain Hardening Material
,”
Eng. Fract. Mech.
,
19
(
4
), pp.
633
642
.10.1016/0013-7944(84)90096-1