Final report for WP4.3: Enhancement of design methods ... - Upwind
Final report for WP4.3: Enhancement of design methods ... - Upwind
Final report for WP4.3: Enhancement of design methods ... - Upwind
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
A Fracture Mechanical (FM) modelling <strong>of</strong> the crack growth is applied assuming that the crack can be<br />
modelled by a 2-dimensional semi-elliptical crack. It is assumed that the fatigue life may be represented<br />
by a fatigue initiation life and a fatigue propagation life. It is there<strong>for</strong>e:<br />
N N N + =<br />
I<br />
P<br />
Where N is the number <strong>of</strong> stress cycles to failure, NI is the number <strong>of</strong> stress cycles to crack propagation<br />
and NP is the number <strong>of</strong> stress cycles from initiation to crack through.<br />
The number <strong>of</strong> stress cycles from initiation to crack through is determined on the basis <strong>of</strong> a twodimensional<br />
crack growth model. The crack is assumed to be semi-elliptical with length 2c and depth a.<br />
The crack growth can be described by the following two coupled differential equations.<br />
da<br />
= C<br />
dN<br />
dc<br />
= C<br />
dN<br />
A<br />
C<br />
m ( ∆K<br />
) a(<br />
N )<br />
A<br />
m ( ∆KC<br />
) c(<br />
N0<br />
) = c0<br />
0<br />
= a<br />
0<br />
Where CA, CC and m are material parameters, a0 and c0 describe the crack depth a and crack length c,<br />
respectively, after NI cycles and where the stress intensity ranges are ∆KA(∆σ) and ∆KC(∆σ). ∆KA and ∆KC<br />
are obtained based on the models in [49] and [50].<br />
The stress range ∆σ is obtained from<br />
∆σ = XWave XSCF ·Y·∆σ e<br />
where XWave, XSCF are model uncertainties, Y is the model uncertainty related to geometry function and<br />
∆σ e is the equivalent stress range:<br />
e 1<br />
∆σ<br />
=<br />
⎣n<br />
⎡ nσ<br />
⎢ ∑ ni∆<br />
i=<br />
1<br />
1/<br />
m<br />
m ⎤<br />
i<br />
σ ⎥<br />
⎦<br />
48<br />
= σ n<br />
The total number <strong>of</strong> stress ranges per year is n ∑ ni<br />
.<br />
i=<br />
1<br />
In the assessment <strong>of</strong> the equivalent constant stress range the effect <strong>of</strong> a possible lower threshold value<br />
∆KTH on the crack growth inducing stress intensity factor ∆K has not been taken into account explicitly.<br />
This effect is assumed implicitly accounted <strong>for</strong> by evaluation <strong>of</strong> the equivalent stress range using the appropriate<br />
SN-curve exponent m.<br />
The crack initiation time NI is modelled as Weibull distributed with expected value µ0 and coefficient <strong>of</strong><br />
variation equal to 0.35, see e.g. [51].<br />
The limit state function is written<br />
g<br />
( X)<br />
= N − n t<br />
where t is time in the interval from 0 to the service life TL.<br />
To model the effect <strong>of</strong> different weld qualities, different values <strong>of</strong> the crack depth at initiation a0 can be<br />
used. The corresponding assumed length is 5 times the crack depth. The critical crack depth ac is taken<br />
as the thickness <strong>of</strong> the tubular member.