Approximate Methods For Fracture Analysis of Tubular Members Subjected To Combined Tensile and Bending Loads
Approximate Methods For Fracture Analysis of Tubular Members Subjected To Combined Tensile and Bending Loads
Introduction
This paper presents a very simple method for predicting the that methods for the purely elastic problem have been around
crack growth and fracture behavior of through-wall cracked for some time now, as summarized in Forman et al. (1985).
ductile tubular members when subjected to combined tensile Discussions of the conditions for achieving /-dominance and
and bending load. The method is general in that it may be the suitability of / a s a fracture parameter for combined bend-
applied in the full range between elastic and fully plastic con- ing and tensile loadings have been presented by Shih (1985)
ditions. The method is essentially a /-estimation scheme tech- and Shih and Hutchinson (1986) by studying the single-edge
nique. Because it is based upon /-tearing theory, it is subject notch specimen.
to the usual limitations imposed upon this theory; i.e., pro- Additional studies based on finite element analyses of the
portional loading is assumed. This of course means that, in single-edge notch specimen subjected to combined tension and
the strict sense, only very small amounts of crack growth are bending have recently appeared (Kaiser, 1985; Sonnerlind and
permitted. However, from a practical point of view, this method Kaiser, 1986). An important result obtained by Kaiser and
gives good predictions for even large amounts of crack growth Sonnerlind (1986) indicates that the value of / is essentially
and errors tend to be in the conservative direction (Brust, 1987). independent of whether tension is applied, then bending; bend-
The new estimation scheme represents an extension of the ing then tension; or both tension and bending is applied pro-
methods presented earlier by the authors for pure bending portionally. This is not intuitively obvious since such loading
(Brust, 1987; Gilles and Brust, 1991; and Brust and Gilles, clearly violates the hypothesis (necessary for valid /-tearing
1993). theory) of proportional loading. This result is important for
For circumferentially through-wall cracked pipes, elastic- what follows.
plastic analysis techniques which do not require full three-
dimensional finite element analysis for combined tensile and
bending loads are scarce. Paris and Tada (1983) have presented
a method which interpolates between the known elastic and The New Method
rigid plastic solutions by using a pseudo-plastic zone correction We begin by making some general comments with regard to
to the elastic solution. /-estimation techniques for combined tension and bending of
Klecker et al. (1986) introduced a method which is very through-wall cracked pipes. Assume that the total load point
similar to Paris and Tada (1983) approach, except it accounts rotation, <j>, and axial displacement (or stretch), 5, may be
for material strain hardening in a rather empirical way. Both written as
of these techniques require numerical integration. Recently, ,NC
Kumar and German (1988) presented a method which is based > = <b- + ( (lfl)
upon interpolating between a compilation of finite element c
S=o + 5 A
(lb)
solutions. The British R-6 method is a method to predict failure
loads for pipes subjected to combined tension and bending. In Eq. (1), the c superscript refers to displacement "due to
However, displacements are not provided. It should be noted crack," while NC refers to "no crack." The contributions to
displacement due to no crack may easily be found from ele-
mentary theory. Further
Contributed by the OMAE Division for publication in the JOURNAL OF OFF-
SHORE MECHANICS AND ARCTIC ENGINEERING. Manuscript received by the OMAE 4>c=ti + 4>cP (2a)
Division 1993; revised manuscript received June 29, 1994. Associate Technical
c c c
Editor: H. Chong Rhee. 8 + 6 e+ 5p (2b)
Journal of Offshore Mechanics and Arctic Engineering NOVEMBER 1994, Vol. 116 / 221
8p = L, a (8)
Tension. Consider replacing the actual circumferentially
through-wall cracked pipe with a pipe with a reduced thickness
since 8ceis linear in stress. bce, the elastic displacement, is assumed
(where te refers to "effective reduced thickness") extending
to be a known linear function of stress. Thus, the only unknown
for a distance d at the center of the pipe (see Fig. 2). The
in Eq. (8) is L\, which is only a function of geometry and
reduced thickness section, which actually results in a material
material properties. It is L't which is estimated via our analogy
discontinuity, is an attempt to simulate the reduced system
compliance due to the presence of the crack. It is assumed that suggested in Fig. 2. The following approach was originally
deformation theory plasticity prevails; i.e., given (Brust, 1987; Gilles and Brust, 1991) for pure bending,
and resulted in a very simple and accurate estimation scheme.
e= ee+ep (5a) The stress intensity factor may be written
e o o KlI=a,Fl(d)^f:^Re (9)
— = — + «( — (5b)
e0 o0 \o, where the t sub and superscript refers to tension, and
P
(5c) a,- 2-KRt (10)
where e0 = a0/E, a and n are Ramberg-Osgood coefficients, R in Eqs. (9) and (10), and throughout this paper, is the mean
and a0 is an arbitrary reference stress. Although this form is radius. A relation between the elastic stretch due to the crack
inconvenient, we use it because of its popularity. Note that if and P may be written, after utilizing Castigliano's theorem
the reference stress is changed from a0 to au a is accordingly (see, Brust, 1987; and Gilles and Brust, 1991, for example)
changed as
%=- 1,(0) (11)
AirtE
«i = a 0 l —I (6)
F,(6) and 1,(6) are provided in Appendix A for completeness.
Note that S£ given by Eq. (11) gives the displacement at the
where ai corresponds to ou a0 to a0. pipe centerline (Fig. 1(a), point O) rather than the rotation
Ilyushins' theorem assures us that, for a cracked body (Fig. point for perfect plasticity (Fig. 1(a), point O'); hence we
1), a unique relationship exists between a generalized displace- obtain a conservative prediction.
ment and a generalized load. For pure stretch, with defor- Now we estimate L', of Eq. (8) using the reduced thickness
mation theory plasticity prevailing (Eq. 5(c)), Ilyushins' section analogy, Fig. 2. If a is sufficiently large and strains
theorem provides small
Fig. 2 Pipe geometry used to estimate the load versus displacement relationship
in the cracked pipe to equal the limit load for the reduced cross
8Z=e„l — (12) section pipe, where ay-is the flow stress; whereupon
te = iP(6) (18)
&c ae Here,
^ °\i (13)
(14) which results if one calculates the rigid plastic limit load with
~a\on} U P applied through O ' in Fig. 1(a). Alternatively, a more con-
The c superscript in Eqs. (12) through (14) refer to the dis- servative estimate for t. is
placement which occurs due to the presence of the reduced e 2
thickness section pipe via the analogy of a cracked section of P(0)=1- -sin I (20)
pipe. 2
By comparing Eq. (14) to Eq. (8), it is obvious that which is the limit load for P applied through O. For the pre-
n-l dictions which follow, Eq. (20) is used.
L,= (15) Bending. Again, assuming a deformation theory plasticity
material law (Eq. 5(c)), that beam theory applies, and en-
for pure stretch. Thus, L', is estimated from the simple analogy forcing compatibility, it is shown (Brust, 1987; Gilles and Brust,
and used in Eq. (8), while the 8ce term in Eq. (8) is from the 1991) that the moment curvature relationship may be written
actual cracked pipe. For bending alone, this type of analogy as
gives an excellent estimate of the moment-rotation due to crack
relationship (Brust, 1987; Gilles and Brust, 1991). The value IB(0)
M" (21)
used for te will be discussed shortly. *P~LB{EO"0- (irR2t)"
To evaluate the tensile component of Jp, the second term with
of Eq. (4) is used. Placing (10) in (8), inverting, and integrating
Eq. (14), one obtains 7T
(22)
U =
1 irR 4k
4- Ea"-1 n+l 2
°HrL, .1,(6)
2wRt
(16) In Eq. (21), I (6) is a function of half-crack angle, 6, and is
B
provided in Appendix B, and R always refers to mean radius.
where For combined bending and pressure, the moment here equals
48Ff(6) J _ bLt_ M+Pe. Also,
H,(n,d) = + (17)
1,(0) Li
All other terms in Eq. (16) have been previously defined. Hence, K=- (23)
we have a simple closed-form relationship of the component
of Jp caused by pure stretch.
To evaluate L,, we force the rigid plastic tensile limit load where T represents the Gamma function. For our purposes
Journal of Offshore Mechanics and Arctic Engineering NOVEMBER 1994, Vol. 116 / 223
Brust (1987) and Gilles and Brust (1991) provide many more
details of these developments. It is sufficient to say that the
only assumption made with this method is the reduced thick- 3 4 5
ness section analogy represented by Fig. 2. This entire devel- Displacement, inch
opment is completely consistent with Illushins' theorem, except
Fig. 4 6-in. pipe, four-point bending load versus displacement; 304
we have reduced the complex three-dimensional problem to a stainles steel at 550°F <
one-dimensional problem, which results in a very simple closed-
form solution. This method was compared to finite element
results and experimental data (Brust, 1987; Gilles and Brust,
1991; Brust and Gilles, 1995) for a large number of cracked
nuclear pipes (including welded cracked pipes), and was found
to give very good predictions of crack initiation and maximum
moment, and very good moment-rotation predictions.
Combined Tension Plus Bending. J is evaluated as the sum
of Je plus Jp. Je is given by Eq. (3) where accurate values of
K'i, Kf are provided from the Sanders solution, and are given
in Appendix A. Jp is determined from the sum of Eqs. (16)
and (26), where the moment includes the Pe term. All dis-
Current Method (LBB.ENG2)
placements and rotations may also be calculated from the Paris/Tada Method
equations provided. Hence, we see that we have a very simple Klecker Method
./-estimation scheme for combined bending and pressure which Experimental Data
growth predictions are made using this method. Load-line Displacement, inch
Fig. 5 6-in. pipe, pressure and four-point bending load versus displace-
ment; 304 stainless steel at 550° F
Results
Figures 3,4, and 5 shows predicted load versus displacements in all respects, except for the method of loading, one can
compared to experimental results for a 6-in. pipe for pressure, compare the load interaction.
bending, and combined pressure plus bending, respectively. In Fig. 3, the displacement is axial and measured or predicted
All pipes are 304 stainless steel and were tested at 550°F. at the pipe centerline. Figure 4 represents loading applied via
Material properties along with experimental data are from four-point bending and displacement represents that at the load
Wilkowski et al. (1985-1988) and are listed in Appendix C. point (total load). The results presented in Fig. 5 are for a
All pipes had initial crack sizes of about 37 to 38 percent of pipe, originally loaded via internal pressure of 2.5 ksi, and
the pipe circumference. Because these pipes are nearly identical then loaded in four-point bending using displacement control.
Journal of Offshore Mechanics and Arctic Engineering NOVEMBER 1994, Vol. 116 / 225
Fl(8) = l+A,l- + B
+ C, (34)
'(*
A P P E N D I X C
Fb(d) = l+Ab + Bb -) +C„ - (35)
Pipe Dimensions and Material Properties for Example
with At through Cb given in Appendix A in functional form. Cases
Then
(1) Pressure Loading. This case represents a through-wall
6Ff (6) = ir\-+I- (2A,)+(2B,)[-J+(2C,)i- cracked pipe subjected to monotonic pressure to failure. A
bladder was placed over the crack to prevent leakage. This was
a nuclear Grade 304 stainless steel tested at 550°F.
+ | - 1 [A2 + (2A,Bt) I - ) + (2A,C, + B2]
The Pipe Dimensions:
Outer diameter 6.62 in.
+ (25,0 - +Ct (36)
\ir/ \7r Thickness 0.507 in.
1
1,(6) = 4 0F,2 (6)d6 Crack size ( - ] 0.386
+
>>c€ H c 'l;
1 a = 4.87
+ -(2A,C, + B1,) (37)
« = 3.88
Let The /-Resistance Curve:
2A, (2Bt)/e\ (2C,)/e J=Jlc+C(Aa)m
7 + + (38)
"-^r ~9~Uj ^rU Jlc = 8.1l k/in.
i- - (>*?) , ( 2 ^ A ) / e \ , (2Alc,+B2)(e C= 14.17 w = 0.6176
(39)
(2) Bending Load. This case represents a through-wall
(2B,Ct) 6\ C,2 6 cracked pipe subjected to four-point bending to failure. With
'o = - (40) four-point bending, the center section with the crack is sub-
4.5 \TT
jected to a constant bending moment. The load in Fig. 4 rep-
Then resents total applied load. This was a nuclear Grade 304 stainless
steel tested at 550°F.
1,(6) =262 1 + 4 | - '«,+ (Ih+Ih (41)
' See Fig. 1 for definitions.
Journal of Offshore Mechanics and Arctic Engineering NOVEMBER 1994, Vol. 116 / 227