351
351
351
AN INTRODUCTION TO
STRUCTURAL RELIABILITY
THEORY
.,
..1
The SHIP STRUCTURE COMMllTEE is constitutedto prosecutea r~earch programto improvethe hull
structuresof shipsand other marinestructuresby an extensionof knowledgepertainingto design,
materials,and methodsof construction.
The SHIP S~UCTURE SUBCOMMllTEE acts for the Ship StructureCommitteeon technicalmattersby
providingtechnicalcoordinaticmfor determinatingthe goatsand objdhws of the programand by
evaluatingand interpretingthe resultsin terms of structuraldesign,construction,and operation.
MHITARY
SEALIFT
COMMANR II s COAST GUARQ
Mr. AlbertJ. Attermeyer CART T. E. Thompson
Mr. Michael W. Touma CAPT DonaldS. Jensen
Mr. Jeffery E. Beach CDR Mark E. Nolt
MARITIME ADMINISTRATION
Mr. FrederickSeibold
Mr. Norman0. Hammer
Mr. Chao H, Lin
Dr. Walter M. Maclean
LT Bruce Mustain
Mr. AlexanderB. Stavovy
. . MERCHANT MARINE ACADFMY
NATIONAL CADFMY OF S E CES
Dr. C. B. Kim ~
L!. s, MVAL ACAD~ Mr. Stanley G. Stianaen
c2.-
Member Agencies:
Address Correspondent to:
SSC-351
December 3, 1990 SR-1310
J%. SIPES
Rear Admiral, U.S. Coast Guard
Chairman, Ship Structure Committee
Tmhnical Report Documentation Paaa
zr-
SSC-351
4. Title ond Subtitle 5. Rwport Dote
JANUARY
1989
ANINTRODUCTION
TO STRUCTURAL
RELIABILITYTHEORY 6. Performing Orgdniznti~n Code
SHIP STRUCTURE
COMMITTEE FINAL
U.S. COASTGUARD
2100 SECONDSTREET,SW, 14. Sponsoring A~mcy Code
This work was carried out under the sponsorship of the Ship Structure Committee
and its member agencies.
16. Abstruct
3 ,<,
L
METfttCCONVHWONFACTORS
*
wtwnYaJKrwm M
Swmlbd whmYou Kmw Muh*by To Fltd symbol
-Ilmws
cm
m MmMun
&m
d
d
d
19
TEMPE
2
t
cm
ABSTRACT i
ACKNOWLEDGEMENT ii
TABLE OF CONTENTS ...
Ill
NOTATION vi
iii
s (j
{
4. BASIC RELIABILITY CONCEPTS BASED ON FULLY
PROBAB~ISTIC METHODS - LEVEL 3: 1(M
4.1 Introduction - Reliability Levels 108
4.2 The Basic Problem - Level 3 la
4.3 The Normal Tail and Margin of Safety 112
4.4 Probability of Failure of a Ship Hull Girder 113
8. SYSTEM RELIABILITY 1%
8.1 Introduction 1%
8.2 General Formulation 1%
8.3 Bounds on the Probability of Failure of a Series
System 198
8.4 Bounds on the Probability of Failure of a Parallel
System
8.5 General Systems m
.,-
7 ~.,-
NOT ATIO~
The following list defines the main symbols appeming in this report.
L= ship length
x= random variable
Px = mean of X
ax = standard deviation of X
fx(.) = probability density function of X
Fx(.) = distribution function of X
EX = mean square value of X
A = parameter of exponential distribution
l,k = parameters of Weibull distribution
co = wave frequency
H(u) = frequency response function
&(w) = wave spectrum
S@) = response spectrum
N,n = number of records or encounters
Zn = random variable representing extreme amplitude of
total bending moment in n-records
(pzJ.) = probability density function of ~
Pf = probability of failure
p~n=l = probability of failure for n = 1
[A1UA2] = union of two events Al and Az
ii = complementary event of Ai
IAlnAz] = intersection of the two events Al and A2
Cg = generalized cost
Ci = initial cost of construction plus maintenance
Cf = cost due to failure
M= safety margin
s= random variable representing strength
z ,= random variable representing total bending moment
g(.) = limit state function or performance function
Xi* = coor ates of most likely failure point in original space
%
vi
Yi* = coordinates of most likely failure point in reduced space
Vx = coefficient of variation of a random variable X
0 = central safety factor
$( ),@{) = standard normal probability density and distribution
fuctions
correlation coefficient
direction cosine
safety index
partial safety factors associated with a random variable
lft = 0.3048 In
1 in. = 25.4 mm
1 psi = 6.894 kPa
1 ksi = 6.894MPa
1 lb-in. = 0.113 N-m
1ton-ft = 0.309 t-m
vii
L INTRODUCTION AND SUMMARY
... . .
11 L,v
vessels, for example, this is done either through the detemnination of a long-term
distribution of the wave loads [5,6,7] or through the evaluation of an extreme load
distribution [8,9,10,11] that may occur in a specific storm condition.
Referring back to Figure 1.1, other loads beside wave loads occur on a
marine structure. These loads may be important in magnitude, though usually
less random in nature (except possibly for wind loads on offshore stmctures). For
example, in the case of ocean-going vessels, these loads consist mainly of
stillwater loads and thermal loads.
Following Figure 1.1, the response of the marine structure to the total
combined loads is determined and compared with the resistance or capability of
the structure. This comparison may be conducted through one of several
relia bilitv methods. Based on these methods, safety indices or probabilities of
failure are estimated and compared with acceptable ones. A new cycle may be
necessary if the estimated indices are below the acceptable ones.
2
1
m
E
rationally including them in the design process. The load and the strength are
thus modelled as random variables.
Figure 1.2 shows the probability density functions of the load and the
strength of the beam in terms of applied bending moment and ultimate moment
capacity of the beam, respectively. Both, the load Z and the strength S are
assumed in this example to follow the normal (Gaussian) probability distribution
with mean values 1.LZ=20,000ft-tin and ps=30,000 ft-ton, respectively, and standard
deviations of 0Z=2,500 ft-ton and a~=3,000 f&tan, respectively.
We may now construct a simple function g(s,z), called the limit state
fuction, which desctibes the safety margin M between the strength of the beam
and the load acting on it, i.e.,
Both S and Z are random variables and may assume several values.
Therefore, the following events or conditions describe the possible states of the
beam,
4
,..
.,.,. :,..
where fs,z (s,z) is the joint probability density function of S and Z and the
domain of integration is over all values ofs and z where the margin M is not
positive, i.e., not in the safe state. If the applied load on the beam is statistically
independent from the beam strength the above equation can be simplified and
interpreted easily as:
where Fs (.) and fz (.) are the cumulative distribution fhnction of S and the
probability density function of Z, respectively, both in this example, are Gaussian.
5
Pf = @(-p) (1.4)
(1.5)
=*
Substituting in equation (1.5) the numerical values for ys, pz, ~s and Cz
given in our simple beam example results in a safety index ~ = 2.56. Equation
(1.4) can be then used in conjunction with tables of standard normal cumulative
distribution function to yield a probability of failure= 5.23x 10-3.
~formation
1. N for Reliability Ana Ivsis of Ma rine % r ucture~
The preceding example and Figure 1.1 indicate that certain specific load
and strength information are necessary for performing reliability analysis of
marine structures. It is mostly in this area that reliability analysis of marine
structures differs from typical civil engineering structures. In this report
emphasis is placed on developing the required load and strength information for
marine structures.
6
function is determined from strip theory using the equations of motion of the
vessel or from a towing tank experiment. In offshore structures, Monsons
equation is usually used to determine the wave load trwsfer fiction.
Short term prediction of the loads is not sufficient for the reliability
analysis. Extreme values and long-term prediction of the maximum loads and
their statistics are more valuable. For this purpose order statistics and statistics
of extremes play a very important role. Gumbels theory of asymptotic
distributions is often used in this regard. In the long-temn prediction, the fatigue
loads, i.e., the cyclic repetitive loads which cause cumulative damage to the
structure must also be considered.
The second major component in the reliability analysis is the strength (or
resistance) of the marine stmcture and the evaluation of its modes of failure. In
this regard several limit ~tates may be defined such as the ultimate limit state,
fatigue limit state and seficeability limit state. The first is related to the
maximum load carrying capacity of the structure, the second to the damaging
effect of repeated loading and the third to criteria governing normal use and .
durability. Each of these limit states include several modes; for example, the
ultimate limit state includes excessive yielding (plastic mechanisms) and
instability (buckling failure).
7
Chapters 4,5 and 6 present three different levels of reliability analysis based
on the load and strength information discussed in Chapters 2 and 3. Examples of
application to ships are also provided in these chapters. Chapter 7 introduces
simulation and Monte Carlo techniques as a tool for use in the reliability analysis.
System reliability, which deals with redundancy of structures and multiplicity of
failure modes is discussed in Chapter 8. Chapter 9 describes a procedure for
fatigue reliability which requires separate treatment from reliability under
extreme load. Several application examples to ships and offshore structures are
given in Chapter 10. The last chapter of the report discusses some shortcomings
and offers concluding remarks. Appendices are given at the end of the report one
of which includes some helpfd information and another describes a computer
program for performing reliability analysis.
The reader of this report can get full benefit of the material presented if
he/she has a background (one course) in basic probability theory and statistics
including probability distributions, random variables, expectation of a random
variable, sampling theory, and estimation methods.
8
1.1 Madsen, H. O., Krenk, S. and Lind, N. C., Methods of Structural Safety,
Prentice-Hall, Inc., Englewood Cliffs, 1986.
1.3 Theft-Christensen, P. and Baker, M., Structural Reliability Theory and Its
Application, Springer-Verlag, Berlin, 1982.
1.4 Konishi, I., Ang, A. H.-S. and Shinozuka, M. (editors), Structural Safety
and Reliability, Volumes I, 11 and 111, Proceedings of ICOSSAR, Kobe,
Japan, 1985.
1.7 Mansour, A., Probabilistic Design Concepts in Ship Stmctural Safety and
Reliability, Transactions, SNAME, 1972.
1.8 Ochi, M. and Motter, E., Prediction of Extreme Values of Impact Pressure
Associated with Ship Slamming, Journal of Ship Research, Vol. 13, No. 2,
June 1969.
9
1.11 Mansour, A. and Lozow, J., Stochastic Theory of Marine Vehicles
Slamming Response in Random Seas, Journal of Ship Research, Vol. 26,
No. 4, Dec. 1982.
1.15 Guesdes Soares, C. and Moan, T., Statistical Analysis of Stillwater Load
Effects in Ship Stmctures, Trans. SNAME, Vol, 96, 1988.
10
2. LOAD INFORMATION REQUIRED IN RELIABILrrY ANALYSIS OF
MARINE STRUCTURES
2.1,1. Definitions:
11
2/
predominantly sinusoidal time history of a measured quantity are
records of a deterministic process [2.2].
Figure 2.1
12
22
due to factors beyond ones control and ability to measure.
The
elevation of the water surface at any time is due to the entire history
of the meteorological conditions in that area and surrounding areas.
Therefore, under given macroscopic parameters such as wind
direction, speed, duration, etc., one cannot predict exactly the wave
surface at the given point. The wave elevation records can be thus
treated as records of a random process.
.+ 1
[
o 100 2CC 300 400 500 600 700 800 930 1000
TIME[SeC]
Figure 2.2
13
23
Each of the above records is called a sample. Some of these
samples are more probable than others, and in order to describe the
random process further it is necessary to give probability
information.
-2
-4
0
x(z) (t)
4 .
2
x(n) (t) O
-2
-1
o
TIME[SHC]
Figure 2.3.
A designer given the record ensemble shown in Figure 2.3 may:
a. Select the largest value in the ensemble and use it for his
new design with a factor of safety. In this case he will make no use
of all information he is given except for one value, i.e., the maximum
value, Orl
At a fixed value, t = tl, (see Figure 2.3) the values. of X(i)(t 1),
which represent wave elevation, can be described by a graph such as
in Figure 2.4. This graph shows the probability density functionl
(p.d.f.) f [ x (tl ) ] or simply f(x) which has the following properties:
F(%)
Figure 2.4.
l-s
1) The fraction of the ensemble members for which the wave
elevation X, treated as a random variable, lies between x and
X + dx is f(x) dx I i-e. ,
P[x<X~x+&]=f(x)dx
p[a<X~b] = f(x)dx
f
a
If, at two time instants t 1 and t2, the wave elevations treated
as random variables are denoted X(t 1) and X(t2 ), or simply XI and
X 2, the probability distribution of these (called joint probability
density function) is given by a surface f(xl ,x2) and has the
properties (Figure 2.5):
1) P[xl < Xl s Xl+ dxl and x2 < X2 <x2 + dx2] = f(xl,x2) dxldx2
~~
az al
f(xl ,x2) dxldx2
16
Ensemble Averages:
.
,.
x,
* .
OENSl T~ FuNCTION
JOINT PQOOAW--l T+
Figure 2.5.
17
For a given p.d.f. f(x), the following definitions will be used for
the wave elevation treated as a random variable X. When g(x) is
simply x, then
+m
E[x] =
--
J xfxdx
() = mean or ensemble average (2.1)
= expected value of X (analogous to
moment about origin or distance from origin to center of mass of a
rod of unit mass)
=
E[x1x2] - E[x11E[x21 (2.5)
(2.6)
19
27
and, thus, both the correlation coefficient and the covariance
A, $~
~i, ~z are zero. This means that independent random variables must
necessarily be also uncorrelated (h mnmrse is no~ necessarily true) .
This implies that the first order p.d.f. f(x) becomes a universal
distribution independent of time and E(x) and r k are also constants
independent of time. .
20
that the autocorrelation function is also a function of ~= tz - tl only
(see definition of the autocorrelation function given by equation 2.4).
z-
\
4
Figure 2.6. Autocorrelation function
21
various statistics of interest by averaging over time instead of
averaging over the ensemble.
f(t)
An ergodic process implies that E[x] =< f>and R(F) =#( &), An
ergodic process is necessarily stationary since < f } is a constant while
E[x] is generally a function of the time t = t ~at which the ensemble
average is performed except in the case of stationary process.
However, a random process can be stationary without being ergodic.
22
32
For ocean waves, it may be necessary to assume the ergodic
hypothesis if there is only one sample function of the process.
.
s(u) = *, f+m R(~)e-iwTd~ - JJ< m
.- (2.14)
* .
23
3?
The mean square (area under the spectral density curve) is the
average of the square of the wave elevation and the root mean
square (rms) is the square root of that value. Physically, the larger
the mean square (or the r.m.s. value), the larger is the wave
elevation and the higher is the sea state.
R (0) = E[x2] =
J
,
S(w ) dti; [ftk] ; therefore
w +(,d
24
4) In practice, a one-sided spectrum can be defined by simply
folding the S( @ ) cume about the w = O axis, that is,
m
J S(L )dw = ~2S(~)dU
a square of the p;ocess, where
= ~~ti
o
)dti = mean
= o (d<o
z, I
5) It can be shown that the area under the ocean waves spectral
density, that is, the mean square is proportional to the total energy
per unit area of the waves surface which is given by:
25
The energy in an increment Sw of wave frequencies at a cenbal
frequency u. is
Pg
(2.15)
Energy per unit area = ~~$ $:, where $a= wave amplitude.
therefore,
In tiis re~rt the energyrather than the anplitude spectral density will
k used.
-<
26
?6
It should be noted that, both the spectral density and the
autocorrelation function are measurable quantities that can be
indeed determined from time history records of ocean waves.
Qb w
27
Figure 2.11. Spectrum of a Wide Band Process,
(A broad range of frequencies are shown in the sample.)
(2.17)
28
x(t)
where & = d%
~d f%,~(., ., . ) is the j.p.d.f. of z and ~
Jk J
Similarly, the mean rate of crossing, that is, the average
number of crossing per unit time with a negative slope (from above)
is .
0
v;(t) = /11 ~ fX, ~ (a, i,t)di
. -m (2.18)
. .
If the threshold level a is zero, the corresponding mean rate of
crossing (from above and below) is
(2.19)
29
b. Peak distribution of a stationary narrow band random process:
dFP(a) du+
fp(a) = da = 1 (2.21)
~ d:
o
As wilI be discussed later, ocean waves can be considered as a
stationary narrow band Gaussian process with zero mean. It can be
shown that for such a process % and & are statistically independent,
a
i.e., the slope ~ is independent of the magnitude K . Therefore, the
j.p.d.f. is given by:
and
m
2
0.
x
=
fo (J2S ((II) da
(2.24)
30
+0
Thus, from equation (2.17) and (2.22) the mean rate of positive
crossing of an amplitude of level a is
1
a.
+ 1 x--r{ +) 2
v = e
a 21 (2.25)
x
Due to symmetry of the Gaussian process about zero v;= >:, therefore
O* a
u = 1
x J+ ( ) 2
a r x
x
(2.26)
and +
Vo= v;= ++--
x
-+(~) 2
fp(a) = ~ e x
(2.28)
x
Both equation (2.25) and (2.28) are important results for ocean
waves. Equation (2.25) gives the average number per unit time
(mean rate) of crossing a wave amplitude of level a and equation
(2.28) gives the p.d.f. of the peaks. In general, the following
important result was obtained: The peaks of a stationary,
31
+/
narrowband Gaussian process (e.g. ocean waves) follow a Rayleigh
distribution with parameter EX given by
m
2
= E [x2] = S(a) do
x = ax I o
= area under the energy (mean square) spectral density (2.29)
1) Over a short period of time (one hour) the wave records can be
assumed to be stationary, relatively narrow-band random process.
(2.30)
a%
2 = /+= X2 f(x)dx = E[x2]
-m
Notice that the variance is equal to the mean square since the mean
of the wave elevation E[x] is taken equal to zero.
2
f~a) = Z . e-~~ a>o (2.31)
E&
where EXis a parameter (see figure 2.14).
32
It has been shown in the previous section that for a narrow-band
normal process, with zero mean, the distribution of the peaks follows
a Rayleigh distribution with parameter
x
Figure 2.13. Probability Density Function of Wave Elevation
33
$(d)
1. Pierson-Moskowitz (1964)
2
s (w). = +.
#(g/v@4
w (2.32)
where:
34
-3
oi. 8.1 x 10
= 0.74
!
s (w)
p
i
I
/ \
I w
.
2. Bretschneider Spectrum
-5 =-AU-4 (2.33)
s(u) =au
35
where x and are given by
P
5
a=
16 4 ( iL/3 )2
P
P-=+
P
-B*+ (2.34)
S(u) = AB u-s e
where
A= 0.25 ( ~,,> )2
in
B= ( 0.817 ~ )+
T= mean wave period
36
46
1) Wind speed (most important parameter)
2) Wind duration
1,
0.
W
0.
trum
37
The angular integration of this spectrum will yield the same
one-directional spectrum as would be obtained from a record taken
at a fixed point.
s (Ld,
p) = s (h)) f(p) (2.35)
Thus,
38
4P
first obtained by Rice [2.1] and then used by Cartwright and
Longuet-Higgins [2.7] and their proof is given in reference [2.1].
fp(a,4) =
4-=
e
2camo
+ w-+ 0
e : .
.
+ kc
- *
[ 1 (2.37)
where
1
t
2= bandwidth parameter = 1 - ~
ml
o m,
9
*(U) =
.0
I
&
-+dz
m
m
=
1 o
n S(w) d- : nth moment of the spectrum n=0,2,4
39
positive fp(a,i)
maxima negative
maxima
i=l
(Gaussian
p.d.f.)
40
5(9
(2.39)
41
5/
1. If the excitation (wave input) is a stationary random
process, the response (output) is also a stationary random
process.
42
5.2
Equation (2.40) gives the input-output-relation in a frequency
domain, i.e., between the spectra of the input (waves) and the output
(vessel response). Similar relation can be obtained in the time
domain between the response time history y(t) and the wave time
history x(t). This relation as well as other relations in the time and
frequency domains are developed in [2.2]. The important results are
given here as follows.
y(t) = (2.41)
J x(t - 0) h(9) d9
D
and the mean of the response E[y(t)] in terms of the mean of the
stationary excitation E[x(t)] (if different from zero) is
43
amplitude. For example, if the response under consideration is the
bending moment, then the bending moment can be calculated for the
vessel in regular waves of different frequencies and for different
vessel speeds. The RAO curves would appear as shown in Figure
2.19.
Af 01 FFE12ENT
5H IP sPIEEIZ:)
where
f= water mass density
g = acceleration of gravity
B= vessel beam
L= vessel length
IH(w)IZ = (~}2
44
where V is the shear force.
and
iii
where
= velocity vector
= force vector
m= mass
n= moment acting on the body
G= angular velocity vector
I = moment of inertia about the coordinates axes
45
5s
2.15]. Once the vessel motions and accelerations are determined, the
shear force and bending moment (or any other loads) can be
computed. The values of these responses (including Ioads) due to
waves of unit amplitude and different frequencies give the required
RAOS. There are several computer programs to perform these rather
lengthy computations, for example, [2.16, 2.17, 2.18, 2.19].
sJd
46
stationary normal process, then the response spectrum characterizes
the process completely. If the resulting wave bending moment
spectrum is narrow-band as is the case usually, then the amplitudes
of the wave bending moment follow the Rayleigh distribution
(equation 2.28) with a parameter Ey related to the area under the
bending moment spectrum.
= J ;:)%,;.
o
Some statistics of the bending moment can be thus obtained:
average amplitude of bending moment = 1.25 ~y
47
Figure 2.21. Heading Angle in Regulti Waves.
Therefore, on integration
co m
Ey =
J
Q
Sy (UJ)&=
JSy
c
(@e)d@e
48
But from equation (2.46) $@e= [1 -@V cos e] $U
substituting in equation (2.47) a
i
s@) = Sj(ti) (2.48)
[i. T
~ v Lzuej
2.2.3 ~
The response of the vessel for all wave components can be then
obtained by integration over p
49
VESSEL VELOLITY VECTOR
/
WAVE COMPONENT
//
PREYAILIN6 WNO
4
DIRECTION
. (+f-J e
(x)= (++ x~o (2.50)
50
and - (%$
F(x) =l-e X>o (2.51)
51
random process which leads to an extreme distribution derived by
Vanmarcke [2.30] is the basis for the fourth method.
E = 1 - m22
o 4
+W
rnn =
J con S(u) do) , n = 0,2,4 (2.55)
%(z)
= [FZ(Z4N-* zzE) (2.56)
where fz (z,&) is the initial p.d.f. of the load peaks. For any band-
width load process, Rices distribution should be used as the initial
52
distribution, thus equations (2.37) and (2.38) hold for fz (z,&) and FZ
(z,E), respectively.
z=
a [(2moln N + h
[
l/ln(~
) 1}1+. (2.58)
which is independent of 5 (for mall a ).
53
n
F (z, G) = exp
%
that is, the asymptotic form is double exponential and the cumulative
distribution itself depends on N. ms is the mean value of the load
if different from zero.
F
%
(Z) =exp-e [
- aN{z I.@
1 (2.60)
- aN (z-UN) -%(2-%)
f (zJ=~Ne .exp -e
+ [ 1 (2.61)
54
The mean and standard deviation of the extreme value Z ~ are given,
respectively, by:
(2.62)
K
Zn = J-
6 aN (2.63)
N
=m+
s
(2.64)
(2.65)
where
N - S
a=.
Em
r o
and
(2.66)
55
The plus sign in equation (2.64) should be used if the mean value ~
is positive in order to obtain the larger characteristic value. It
should be noted that both A and ~ contain u ~ as defined in (2.66);
therefore, an iterative procedure must be used for determining u ~ .
To start the iterative procedure an initial value for u ~ is necessary
and may be taken as u ~= rn~t~~. The corresponding values of @ , ~ ,
+(-u) and @ ( ~ ) can then be determined. Equation (2.64) is
then checked to see if the right side is equal to the left side,
otherwise a new value of u ~ equals the right side of equation (2.64)
should be used in the second step of the iterative procedure. Three
or four steps are usually sufficient for convergence.
P max(X(t);
[
OSt S~)sx
1s (2.67)
where
56
2
--( )
lx
- m~
2%
r (2.68)
<=voe
and
v = ~~ ltsec
O ~ J m. (2.69)
\~wlj
FX(X) = exp
I -vOTe
(2.70)
57
Fx (x) = exp -v. T -e
(2.71)
L
(2.72)
()~q~l
Distributions given in (A), (B), (C) and (I)) above, are valid for a
load process represented by a stationary Gaussian process of any
band width. Reference [2.33] shows the corresponding equations for
the special cases of a narrow-band process ( E = o or q = o) and a
wide-band process ( & = 1 or q = 1). It should be noted that the
narrow-band case gives a conservative estimate of the extreme wave
load distribution and the resulting equations may be used for values
of ~ Up tO 0.60 shce they are insensitiveto E in the range O to 0.60.
58
b!?
conditions for a period of one hour. The following parameters were
computed for an earlier application given in [2.33]:
ms = 669,037 ft-tons
m. = 286,300 ft-tons
r
N=, 60 X 60
= 276.9
13a
59
Based on these results, one surprising conclusion can be drawn.
All extreme value distributions of the wave loads considered produce
similar results even though their basic assumptions and derivations
differ drastically. In fact, if one inspects the equations representing
the cumulative distribution functions of these distributions one sees
that these equations are not similar in form and may conclude
erroneously that they would produce very different results.
60
.001.01 .10.20 .30 ,50 ,70.m .90 .95.97.98 S .595 .997s 999
Ptiabilily CDF
co$
10 lhuONM, A C, ti O
00
09
47
0.6
0.$
04
02
02
0.1
0.0
--,- .-*,- -ti---
--m
61
pd. f.
xlo~
I DMnbuhonsA andC
5,0
............... DMribuuonB
4.0
1
3.0
20
10
0.0 I
700 900 llIJO 1300 1~ 1700 1goo 2100 2300 2500 27OO I1OOO
2.5 ~ ilities:
62
72
P[w=w] = FX w-l(x) [1- FX (x) ]
If R = L, then NE(x) ~ e- 1
In the latter case where failure during life of, say ,20 years, is
considered, and the return period is 100 years, then the non-
encounter probability, from (2.75) is 81.8%. If the return period is
increased to 1000 years; the non-encounter probability becomes
98.0%.
63
2.6 Stochastic co mbination of Loads on a Marine Structure
A closer look at the rapidly varying part shows that it also can
be decomposed into components. Figures 2.27 and 2.28 illustrate
64
74
records taken over shorter periods of time (larger scale). TWO main
central frequencies appear in these records. The smaller central
frequency is associated with the loads resulting from the motion of
the ship as a rigid body (primarily heave and pitch motions). This
lower central frequency is, therefore, close in magnitude to the wave
encounter frequencies for wave length nearly equal to ship length.
.* *
.>
65
Figure 2.28. Decomposition of a stress time history of an ocean going
bulk carrier.
The rigid body and the high frequency responses do not always
occur simultaneously in the same record. Quite often only the rigid
66
body response appears in a record; particularly, in that of a smaller
ship which has a high two-node mode frequency. Occasionally, only
the high frequency springing response appears in a record when a
ship is moving or resting in relatively calm water. In particular, long
flexible ships with low natural frequencies as those operating in the
Great Lakes do occasionally exhibit such records when operating in
calm water or in a low sea state composed mainly of short waves.
Under these conditions, a long ship will not respond as a rigid body
to the short waves, but the two-node mode frequency of the hull can
be sufficiently low to be. excited by the energy content of these short
waves. Figure 2.29 (from reference [2.36]) shows a measured
response spectrum of a large Great Lakes vessel where the response
is purely in the two-node mode S and higher frequencies with no
rigid body response appearing in the spectrum. The figure shows
that response at higher modes than the two-node mode can be
measured, although small and relatively unimportant in most cases.
13 ?
:
?.
z
\
10
~mMs. .*0*sl
5$
I* Mom
2-WOC ,
G~
o 0.3120 0,7215 09165
FREOUENCY!HERTZI
67
Two main steps should be used in the procedure for combining
primary responses of a vessel.
SwLL
r
--+_ ____ ____ ____ _l
68
-78
The components of the system represent the components of the
load response of the ship to waves, e.g., the low-frequency wave-
induced responses which consist of vertical, horizontal and torsional
moments, the high frequency responses such as the springing loads .
It is required now to determine the sum of these component
responses, i.e., to determine the output taking into consideration the
proper relations or the appropriate correlation of the response
components,
where
hc(~) = ~ aihi (T)
i=l
(2.77)
hi (z) is the impulse response function of each linear system, i.e., the
response of each linear system to unit excitation multiplied by time
(response to the Dirac delta function). hc (%) is a composite impulse
response function which sums the responses of the individual
components.
69
It should be noted that the impulse response functions of the
individual components hi (~) may or may not be easy to obtain
depending on the complexity of the system. With suitable
instrumentation, it is sometimes possible to obtain a good
approximation to hi (t) experimentally. For the ship system, hi (~)
can be determined for most load components.
Hi(u) = fhi(t)e-j~tdt
o (2.78)
H=(u) = ]rnhc(t)e-j(t
cit
o
=i~laiHi (u)
(2.79)
[Syyfw)li = sxx(~)H*i(~)Hi(~l
= SXX(W) lHi(w) 12
(2.80)
70
[R. A.O.li = [P(u)] = IHi(u)j
i
n
= sxx(~)~ ~ aiaj Hi (U)H*j (w).
(2.81)
i=l j~l
The double sum in equation (2,81) can be expanded such that a final
expression for the total response spectrum Sy y (OI), which combines
the individual response spectra, may be written in the form:
Hi(m)H*j (10) = o
(2.83)
71
81
then the second term in equation (2.82) becomes zero and the load
components are uncorrelated. In this case, the total spectrum is
simply the algebraic sum of the individual spectra of the load
components, modified by the factors ai . Furthermore, if the wave
input is considered to be a normal random process with zero mean,
as usually is the case, then the respective output load responses are
jointly normal and are independent. Thus, the total response, in this
special case, is a zero mean normal process with a mean square value
given by:
= (iyy (u) dw
Y 2
i~lai
;lHi
(w)I sXX(w)dw
o (2.84)
= (Syy (u)dw
aY2
(2.85)
72
In connection with equations (2.84) and (2.85) it should be
ei~ that the usual R@ig~ multiplier used to estimate certain
average quantities, such as average of the highest one third,one tenth
response, etc., are not generally applicable in the case of the
combined response, since these multipliers are associated with a
nmrow-band spectrum for which the amplitudes can be represented
by a Rayleigh distribution. It is only in the case where each of the
load component responses is narrow-band and each happen to be
closely concentrated around a common central frequency w;, that it
would be reasonable to conclude the combined response SYY (w ) is,
itself, a narrow-band process.
(2.86)
where
,,
s~,= correlation coefficients of individual load
d
components defined by,
73
U
pij =+rnHi(u)H*j (u)SxX(ti)dW
ijo
(2.88)
74
given by equation (2.77) form the basis for determining the
combined load response. The question of whether the time or the
frequency domain analysis should be used depends primarily on
what form the required input data is available. In general, most
wave data and practical analysis are done in a frequency domain,
although in some cases where slamming loads are a dominant factor,
it may be advisable to perform the analysis in time domain.
SEu
In this step, the stillwater and the thermal responses should be
combined to form the mean value for the rigid body motion and
higher frequency responses. The stillwater and thermal responses
are weakly time-dependent variables so that in a given design
extreme load condition they can be considered constants, say over
the duration of a design storm. Therefore, theses two responses can
be treated as static cases and can be combined for one or several
postulated design conditions without difficulty. Alternately, if
statistical data are available for each of these responses, the mean
and variances of the combined response can be easily determined.
75
gr
m. OF
Oxuulums
*UWO
nvc.
I
t
Lf e 20
!**
Figure 2.31. Histogram of stillwater bending moment of a
containership.
*5I Iol
76
These mean values together with a set of standard deviations,
which can also be estimated from the histogram, can be utilized to
determine the extreme total moment using a statistical approach.
77
Generally, large tankers traveling in oblique seas may
encounter horizontal bending moments of the same order of
magnitude of the vertical bending moments. Therefore, the
combined effect of the vertical and horizontal moments can be
critical under certain conditions. The distribution of the primary
stress in the deck as a result of the combined effect becomes non-
uniform and assumes a maximum value at one edge. The combined
stress at the deck edge, 8C is given- by
(2.89)
where,
.
GC = combined edge stress
.
where,
s
Iz=$
h
78
Applying equations (2.86), (2.87) and (2.88) and extending the
results for the case of a two-dimensional sea spectrum, the mean
square value of the combined response ~zM~ can be written as (see
equation (2.86)):
where G*M~ and ah h are the mean square values of the vertical and
horizontal bending moments, respectively, given by equation (2.87)
as6:
(2.92)
Mh2
= ;$
~ o
/mSxx(u, U)IHh(w) 1 dwdu
-7 (2.93)
= 1
vh
Mvw
77
~ &x(wrU)H, =(w)H*h(m)dwdu
--
2
(2.94)
79
The Response Amplitude Operators IHV(U )1 and IHh (co)l and
more generally, the system functions Hv (OI) and Hh (OI) can be
determined from any typical rigid-body ship motion computer
program.
s
=K=$
az h
.,
(rlns) = ~MC
c
v
(2.95)
where,
81
.
ql
m
= Sxx(m) I H~ (o)) 12 (2.98)
f
0
=
Y Vs correlation coefficient
and Hv(m) and H~(u) are the complex system functions of the vetical
wave moment and springing moment respectively.
The system
response functions can be computed using computer programs which
take into consideration the effects of ship flexibility in the response,
such as the springing -seakeeping program SPRINGSEA, [2.3 9].
Applications of equations (2.97), (2.98) and (2.99) to several Great
Lakes Vessels where springing is important is given in references
[2.38] and [2.39]. These equations together with equations (2.76) and
(2.77) for the time-domain analysis have a wide range of
applicability to any two or more dynamic random responses
including, combining primary and secondary~ responses [2.38],
vertical and torsional moments for ships where torsion is important,
high frequency loads with vertical and horizontal moments, etc. In
all of the cases, the coefficient ai must be appropriately determined.
When determining the various statistical averages from the
combined r.m.s. values, the more general distribution given by Rice
for the peaks should be used instead of the usual Rayleigh
distribution.
82
qz
20
r
15
10
5-
Head
0 45 90 135 180
Seas
HEADING (DEGREESI
w o
ul -
Z:
k
:u
15- HIGHSCAS
10 -
Q 5 -
. M&
Head
Seas
o~
o 30 60 90 120 150 180
HEAOING(DEGREES)
83
93
:2
v
I I 1
0[
() 10 20 30 40
SIGWICMJ1 WAVC MEICI-IT [FT)
84
2.1 Rice, S.0., Mathematical Analysis of Random Noise, Bell
System Technical Journal, Vol. 23, 1944; Vol. 24, 1945.
Reprinted in N. Wax, Selected Papers on Noise and Stochastic
Processes, Dover, New York, 1954.
2.4 Pierson, W.J. and Moskowitz, L., A Proposed Spectral Form for
Fully Developed Wind Seas Based on the Similarity Theory of
S.A. Kitaigorodiskii, Journal of Geophysical Research, Vol. 69,
December 1964.
2.9 Salvesen, N., Tuck, E. and Faltinsen, O., Ship Motions and Sea
Loads, Trans. SNAME, 1970, P.250.
85
2.10 Abkowitz, M., Stability and Motion Control of Ocean Vehicles,
M.I.T. Press, Cambridge Mass., 1969.
2.13 Tasai, F., Ship Motions in Beams Seas, Research Institute for
Applied Mechanics, Vol. XIII, No. 45, 1965.
2.14 Tasai, F., On the Swaying, Yawing and Rolling Motions of Ships
in Oblique Waves, International Shipbuilding Progress, Vol. 14,
No. 153, 1967.
2.19 Meyers, W. G., Sheridan, D.J. and Salvesen, N., Manual -NSRDC
Ship Motion and Seaload Computer Program. NSRDC Report No.
3376, February 1975.
86
2.20 Wahab, R., Amidships Forces and Moments an a CB = 0.8 Series
60 in Waves from Various Directions, Netherlands Ship
Research Center TND Report No. 1065, 1967.
87
2.30 Vanmarcke, E.H., On the Distribution of the First-Passage Time
for Normal Stationary Random Process, Journal of Applied
Mechanics, Trans. ASME, March 1975, pp. 215-220.
2.33 Mansour, A.E., A Note on the Extreme Wave Load and the
Associated Probability of Failure, Journal of Ship Research, Vol.
30, no. 2, June 1986, pp. 123-126.
2.35 Little, R. S., Lewis, E.V. and Bailey, F. C., A Statistical Studv of
Wave Induced Bending Moments on Large Oceangoing Tanker
and Bulk Carriers, Trans. SNAME, 1971.
.
2.37 Stiansen, S.G., Mansour, A. E., Jan, H.Y. and Thayamballi, A.,
Reliability Methods in Ship Structures, The Naval Architect,
Journal of RINA, July 1980.
88
2.39 Stiansen, S.G., Mansour, A.E. and Chen, Y.N., Dynamic Response
of Large Great Lakes Carriers to Wave-Excited Loads, Trans.
SNAME, 1978.
89
3. STRENGTH INFORMATION REQUIRED FOR RELIABILITY
ANALYSIS OF MARINE STRUCTURES:
90
po
1. Objective uncertainties. These are uncertainties associated
with random variables for which statistical data can be collected and
examined. They can be quantified by a coefficient of variation
derived from available statistical information. The variability in the
yield strength of steel is an example.
91
1 Uncertainties associated with the degree of effectiveness of plating
due to shear lag effects [3.3, 3.4].
92
@
quantities such as resistance of materials, weight of materials, and
geometric parameters of a section.
93
b. Failure due to instability or buckling of a panel longitudinal
stiffeners (flexural or tripping) or overall buckling of transverse and
longitudinal stiffeners of a grillage.
= x
6
x x
The estimated mean value may not be totally accurate in
comparison to the true mean value because of the sample size n. The
sampling or model emor in estimating ~ is
Ax=
;
and x -x
6
=*(X:+X:)
f(x)
m.
u-xl
xl x~
2= +(X1+X)
f(x
I
2
Xu-xl
u
Figure 3.2. Symmetric Triangular p.d.f.
96
[oL
If there is a bias towards higher values then the upper triangle
distribution shown in Fig. 3,3 is more appropriate, and in this case
z= +(X1+ 2x
u)
x x
6Y=A(
F 2XU
+X: )
f(xJ
Lower
Triangle
J
1 \
:
I
x
xl %
Figure 3.3. Upper and Lower Triangle p.d.f.
and
x -x
8
+:++ )
If a normal distribution is assumed (Fig. 3.4) with specified
limits of ~ k standard deviation, then
1
%= ( l.+X )
2
and
x-x
u
6X
-& ( ~ U+x
1
1
)
97
[01
f(x)
x
xl x~
tn) (3.1)
. . F*)
(3.2)
98
= k @s/tk,)t
i-l
g,,z
(3.3)
where ~: and ~fiAare the mean and the variance of the random variable
E: The p~tial derivatives in equation (3.3) are to be evaluated
at the mean value Et . Equation (3.3) can be normalized and written
in terms of the coefficient of variations (COV):
(3.4)
where 6%= ~ is the strength COV, and s~i=~: /~L are the COV of &i and
the partial derivatives are to be evaluated at the mean values.
6== z +
P
-
))* .-
. P.
as
lJ ( aci )(*)UU
J
E.
1 J (3.5) .
i J
99
Youngs modulus, ship steel plate thickness, ship steel corrosion
rates, residual stresses, and fabrication tolerance.
100
Fj is the specified yield stress for the steel grade used. The
weighted average of the COVS for the data presented in Table 3.1 is
0.089.
101
H1
method is prohibitively expensive, which precludes gathering
statistically significant amounts of data. Alpsten [3.14], for example,
reports that surveying a single plate for residual stress took 140
hours, not man-hours, but the time it took a team to collect the data.
Two representative papers which give results and empirical methods
for predicting residual stress are included in the list of references
[3.21, 3.22].
~ While no extensive
search was conducted on this subject, it seems safe to say that the
literature on the subject is limited. Basar [3.17] in his survey of
structural tolerances in the U.S. shipbuilding industry, states that
102
u-
Tablcg.d Ylald alrw@h data
No. of Me-n,
Test Samplcy psi Cov [)istribu;ion RHIIirrkh
~icld stress Mw 0.091 resumedhKnortd mill test 1 con~inmcnl vwi SA537 CrR
} iekl slrw :1964)0 0.10 eumcd Iutinormil mill tCSt 1 cmrtainmenl vessel SA36 :
Yield strem :15W1 0.042 extreme Type I cold straightenedsha~ HE~B
Yield stress
J 36583 0.0922 extremeTyIM I coldstrai~htenedshansHE$KW)B
.-
Yield she~ 19857 40073 0.103 extremeType I mill lmta
~itld ~trens 19217 43475 0.099 xtreme Type I mill laatn
\ irld strew 11170 57fr16 0.057 rxtreme Type I mill teats
Yield stress 2447 63.336 0.054 extreme Type I mill tsstn
upper ... mill tuts
yield poim :{974 40m 0.087 hrgnormal impi~ log nnrmal. refers tn Freudenthal,
AS(X, vol. 121, 1956
Yield strm 400 44mm 0.11 ...
Yield stress 33 36091 0.059 .. . 1948 tests, AM Class A platm 71$and VI in,
Yield pnint 79 34782 0.116 .. . 1948 tasta. ABS Class B pks /16.5/u,1l/l& 3/4,
% k %6,1h.
Yield puint 33831 0.081 1948 tests, ABS Class C pkk 1$E,13/ISi1/4,
1%. 11/1in.
Yield strength 39 34850 0.044 normal Yj-in. ABS B Steel
Yield strength 36 35(WI o.(mg .,, I /,.in. Af3S C Steel
Yield point :1124 39:160 o.07n . ASTM mill tests
Yield point 35 42900 ().10:{ mill hmt i - 1(XKIpin./inJscc
Yield .,.
S1r?mIrvd 35 41X4-I 0.1m snmt specim?ns w nlriIvr. wmIIliIled mill lPSI
Ywld strew+Ievd 35 :14m 0.121 wme SfMClmCllSas abe, stub dUInn lest.
stitic
strainraw
Yield stress 0.2% 50 61780 0.034 not normal cold-rolled
rmle
not lognormal
Yield stress0.2% 38 40530 0.045 not normal hot-roildrds
not ~gnormsl
Lnwer
yield 22 30923 0.110 .. . arbon structural ?4- nd I in. plata and /s-in.
Lower yield 20 42675 0.079 ... loflfi$y +4- and l-in. plata
bwer yield 10 M 290 0.068 ... low-alloy Yt.in. late
Yield point 120 35 OEO 0.038 ... ASTM A7-55T $JF barns, flanges
Yield point 39079 o.od4 .. . ASTM
A7-55T, WF hama, web
Yield point 38(WI 0.0261 .. . ASTMA7-55T,WF bea~, cover platea
Yield stre~ 41800 0.10 lo~normal Mill tek 4 different contimmcnt VISAS,
assumed SA516 GR70
Yield stress 0.066 Iognormal Mill tut. 1 con~inment VHS4 SA516GR70
assumed
P
*
. .
Tabla 3.lUtllmalo alrangthdda
No. of
T-1 Samples Mean, kai COV Distribution Remarks
Tcnaion 58.291 0.043 ,.. cold straighten shapa
Tens;on 3; 57.909 0.089 ,.. cold straighten shape
Tension 84.039 0.1124 . annealed, alloy Wed
Tension ; 124.9 0.1796 ... quenched, alloy steel
Ten~ion 60.405 0.0719 ,.. nominal maximum straas, varioua platex. st~ctural skel
Tcnaion ;: 73.525 0.074 . nominal maximum atreca, various platac, low-alloy s*1
Tension 80.39 0.109 ... nominal maximum stress, various plates, low-alloy 9tccl
Tension l;: 62M 0.0226 ... ASTMA7-MT,WF Iwama, flatrgm
Tension 64.33 0.0341 ,., ASTM A745T, WF barns, web
Temion u 60.84 0.0241 ... ASTMA7-55T,cover plates
Tension 39:; 64.27 0.0703 ... mill Wta
Tensile strength 59.27 0.044 ... 1W8testsABSClass A plataq h % in.
Tensile s~ength 79 0.091 ... 1948 tasta ABS Class B plataa, /1~ % 11/1~,~4,1%6,7A, 1%6, 1 in.
Tensile strength 13 R% 0.051 .,. 1948tsaraABSClass~ plata, 1Ylt. l?I~, 1~, 1%,11Ain,
Tension 6257 0.044 norm41
Tensile
strength % o.&7 normal ~;#A&B&&,
103
Ta&+a 3.3Ywtgs
modulus
dints
1N04
i}f
T=t Samples Mean, ksi Cov Distribution Remarks
Tcnaion 104 30.0x lW O.owln ... structural a~l from brid es
Tension 19 26.98x 102 0.0M9 ... variousM.A lloys,annm! ad and quenched,nd
drawn samplm
Tension 22 29.50xl@ o.(w2 ... structural steel
Compression 29.49x 1(P 0.0146 ... structural akl, same ssmpla as above
Tension :: 29.59x l@ 0.CN156 ... low alloy
ComprSion 20 29.Mx ltP 0.IJ370 ... low alloy, same sampl~ as abve
Tension 10 29.56x l@ O.(HH ... low Slloy
Compression 29.61X 10s 0.01106 ... low 4UOY,same aamplm u ake
Tension $ 29.42X I@ 0.01% ... variom-sixe a~ens from Y4,~ and 1in. plate,
sune WI& atn.rctural 3W1
Ten~ and stindard 94 31.20 x 10s 0.060 ... sttucturd steal
..-. .. -
N-1
REmmcEs
105
3.10 Pardine, C. and Privett, B,, Statistical Methods for Technologists,
the English University Press, London, 1964.
3.18 Purlee, E. L., Leyland, W. A., Jr., and McPherson, W.E., Economic
Analysis of Tank Coatings for Tankers in Clean Service, Marine
Technology, Vol. 2, No. 4, Oct. 1965.
106
[(L
3.19 Josefsson, A. and Lounamaa, K., comparative Comosion Tests
on Steel Plates Rolled from Continuously Cast Slabs and Rolled
from Mold Cast Ingots, Corrosion, Vol. 26, No. 5, May 1981.
3.21 Bjorhoude, R., Brozzetti, J., Alpsten, G.A., and Tall, L., Residual
Stresses in Thick Welded Plates, Welding Journal, Vol. 51, No.
8, Aug. 1972.
3.22 Nagaraja, N.R. and Tall, L., Residual Stresses in Welded Plates,
Welding Journal, Vol. 40, No. 10, Oct. 1961.
107
4. BASIC RELIABILITY CONCEPTS BASED ON FULLY PROBABILISTIC
METHODS - LEVEL 3:
108
(Rosenblatt transformation (1 969)) can be made to obtain equivalent
uncorrelated normal variables, thus, approximate probability of
failure may be determined. Similarly, certain approximation can be
made for nonlinear performance functions.
109
The probability of failure or the probability of reaching a
specific limit state is determined from:
fdxl, x 2 . . ..xn) ~~
(4.2)
110
pf z
I I s,2 (S,2) d= dz
[(s,2):s-250] (4.6)
s)
111
P~=l-pf
(4.8)
M= g (X1,X2) =s -z
(4.9)
o
Pf = PIM<O] = fM(m) dm FM(0)
f
.m (4.10)
fM(m)
1
area= pf
(4.11)
and
112
2 2
M=%+=: (4.12)
-
deviation, i.e., N(O,l). Equation (4.10) then yields
PM
Pf =FM(0)=#
[1 M
1 - 4 (p)
(4.13)
where PM
$z
M
is called the reliability or safety index, Note that the probability of
failure decreases as the safety index increases.
P
and
P
:=% (4.14)
P= q-~) =1-r(d)
f (4.15)
For example, if ~ = O, Pj= 0.50; if /?= 1,?} =().16 ~dif ~ = 3.1, PF= l(j3
The relations k- 6 and p ~ for other distributionof S and Z with
a mrgin gi~ by equation (4.9) are plottedin section11,1.2of @p31diX 1.
113
623
this case, the total bending moment) acting on the ship.
The total
bending moment consists of stillwater bending moment and the
extreme wave bending moment as developed by one of the methods
described earlier.
114
@zm(,, = [1- ,-(s%=)]9,> ~
. 0 otherwise (4.17)
The above results are extended to the case when the variability of
the stillwater bending moment is not neglected. The stillwater
bending moment is considered to follow a normal probability law
given by
(4.18)
Zn=T+Yn (4.19)
. .
(4.20)
115
Differentiating the last equation of (4.20) with respect to z to get the
pdf of Zn, we get
and
d--
(4.26)
116
FR(r) =J+. o
42.(4-F~(rz)& (4.27)
pf can be written either in the form (see also equations (4.4) and
(4.5))
(4.29)
or
+.
p/=l -
J0 @zm(z)Js(z)dz (4.30)
(4.31)
F~(s) -
d
I j~(8)d!l- *=
-- ()s u (4.32)
Two cases will be considered next. First, the case when the
stillwater bending moment - is regarded to be deterministic of value
m o. Using equation (4.17) for Ozn(z) and equation (4.31) for fs(z) in
equation (4.30), and noticing that z > m., the probability of failure
can be written in the form
117
121
A.s
0
[1 - #r+/w]~ ~oandpt~ 1
that is to say, as the ship encounters wave loads for an infinitely long
time, failure will eventually occur with probability equal to 1.
FREQUENCY
OF OCCURRENCE
f
f
SAGGING HOGGING
P*
LOAD BENDINGMOMENT OR ULTIMATE
) BENDINGSTRENGTH
p,,*d - [1-++)]
u-mo d
+e -T%.*
( P - ~
..
u
@
k ) (4.35)
= Pir+ Pf = Pf
= ;+= , S+-
f-+= Jo @/k) -l
*[(Z - $J/ukiz
(4.36)
Prl.-l - (1/%)v2/@
s J ~
.0
119
[27
w(z:-:)v(
=)& 38)
where equation (4.38) can be further reduced by letting c~ =0 and
equating it to equation (4.35) with m. = m. If the fiist term in (4.35)
representing the probability of failure under stillwater bending
moment is neglected, then
v =mf =1.0
( U* )
P m? u; + 4.O(U,+ u) (4.40)
For the more general- case when n >1, Bernoulli trials may be
assumed, and the total probability of failure, pf, can be written in the
form
n
Pf=l-[l-pfln=ll (4.42)
120
Equation (4.42) is applicable for both short and long terms, and
also for deterministic and random stillwater bending moment.
121
individual failure probabilities provided that the distributions are
determined separately. In practical application, however, the total
probability of failure is controlled by the direction of the stillwater
bending moment (just as in the deterministic approach). Thus if the
stillwater bending moment is a hogging moment, the total probability
of failure is simply equal to the probability of failure in that mode
since the probability of failure in the sagging mode will generally be
very small. Sum naml vessels ~r, may Ha- fmm this rule.
122
\32
4.1 Cornell, C.A., A Probability-based Structural Code, ACI Journal,
Dec. 1969.
4.4 Hasofer, A.M. and Lind, N. C., Exact and Invariant Second
Moment code Format, Journal of the Engineering Mechanics
Division, ASCE, Vol. 100, No. EMI, Feb. 1974, p. 111.
123
5. LEVEL 2 RELIABILITY ANALYSIS
M= S-Z (5.1)
124
M-pm
G=
u (5.4)
m
Pf = FM (0) = F~
[ ~: 1= G -d) 5-5)
12
Se section Al. 2of Appmdix 1 for the relationship h~ @and the
probahili~ of failurefOr severaldistributions.
*
where z ~ is the linearization point, and the partial derivatives are
evaluated at that point. In the MVFOSM method the linearization
point is set at the mean values ( El, -Zz -- --- ~n ).
(5.10)
(5.11)
~ = g(x],xd= g(s,z)= s - z
then applying equation (5.9) and (5. 11) for determining the mean
and variance, one immediately obtains identical results as given by
equations (5.3). Equations -(5.2) and (5.5) follow accordingly. This
method is called the ~VFOSM method because the linearization of
the limit state function takes place at the mean value (MV); only the
first-order (FO) terms are retained in Taylor series expansion, and up
to the second moment (SM) of the random variables (means and
126
131)
variances) are used in the reliability measure rather than their full
probability distributions.
M= UgS-UzZ+P~-lJz =0
, .
Failure
~gion
s
o
Figure 5.1. Limit State Function in the Space of Reduced Variates
The region on one side of the straight line which contains the
origin O represents the safe state (M > O) and the other region
represents the failure state (M ( O). Thus the distance from the
origin to the line M = O can be used as a measure of reliability. In
127
131
fact, from geometry, the minimum distance D shown on Figure 5.1
is given by
Notice that D is equal to the safety index ~ for the case of the
normal variates and linear limit state function discussed earlier, i.e.,
for this case
Pm
$ =D= a
1++
and the probability of failure is thus
Pf = 4 (-D)
128
then PF = FM (O)* still given as before by equation (5.5); however,
when r h and r. are computed from (5.9) and (5. 11) and
substituted in P
m
$ =
0
(5.12)
m
the following is obtained
!
P; - P:
$=
4p: 0; + 4p2Zz 02 05 (5.1s)
[ 1
129
13q
x. - ii.
1 1
Yi =
u
x.
1
(5.14)
(5.15)
(5.16)
(5.17)
G(.) is the
*
failure surface in the reduced space, and , are coordinates of the
Y
point closest to the origin in the reduced spkce (the checking point).
All partial derivatives are evaluated at the checking point. In effect,
this procedure is equivalent to linearizing the llmit state function in
-. the reduced variables space at the checking point and computing
P associated with that point.
130
Suppose that a simple beam is subjected to random loading that produces a
maximum stress with a mean value ~Z = 20,000 p.s.i. and a standard deviation Oz
= 2,500 p.s.i. The beam is made of material of mean yield strength P6 = 30,000
p.s.i. and a standard deviation cr~= 3,000 p.s.i.
The following three limit state functions are considered. They all represent
failure of the beam and, therefore, should yield the same value of the safety index
if the method used to determine the safety index is consistent (i.e., invariant to the
formulation of the limit state function). The three limit state functions are:
MI = S-Z (5.19)
w = ~+yz (5.21)
Where the margin means ~m and its standard deviation ~m are computed
for each limit state function using equations (5.9) and (5. 11), respectively. The
resulting ~s for the three limit state equations (5. 19) to (5.21) are:
131
k-vz
h = (5.23)
(0s2 + aza) 1/2
~2.p;
132 = (5.24)
(4ps26s2 + 4p~20.2)1/2
b p@ - hp~
133 = (5.25)
(%% 2 + GZVpz 2 )1/2
Substituting in the above equations the values of ps, ~z, as andaz, one
obtains the following results:
The values of Ps are not the same indicating that the MVFOSM method is
not invariant to mechanically equivalent formulation of the same problem. Notice
that we have not made use of the distribution information (S: Normal and Z:
Weibull) in the calculation. of the ~ values.
If the probability of failure for each limit state is to computed from pfi = @(-~j,
an emor will result since this equation is valid only if all the random variables are
normally distributed and the limit state fmction is linear. Let us, however, use
this equation in order to compare the results with those obtained by a more
accurate method described in the next section entitled Inclusion of Distribution
Information. Using,
(5.26)
132
b, The Hasofar/Lind Method
Equations (5.15) to (5.17) are applied to determine the ~ value for each limit
state function according to this method. For the non-linear limit state functions,
these equations must be solved iteratively since the evaluation of the derivatives
required for calculating ~ will depend on the coordinates of the most likely failure
point which are unknown. An iterative procedure would be simply to assume
values for the most likely failure point (e.g. mean values of the variables) and to
evaluate the derivatives of the limit state function at that point as required by
equation (5.17). Equation (5.17) is then substituted in (5.16) to obtain a set of
coordinates y? which will be a function of the unknown ~. These coordinates are
substituted in (5.15) and the resulting equation is solved for ~. The obkined ~ is
then used in (5.16) to obtain a new set of coordinates of the most likely failure
point. The procedure is repeated until convergence is obtained. The procedure
will be described in more detail in Chapter 6.
The procedure was applied to the limit state equations given by (5.19) to
(5.21). The linear limit state function (s.19) did not require any iteration; the
second limit state function given by (5.20) required five iterations and the last one
(5.21) required seven. The following results were obtained.
The probability of failure calculated from equation (5.26) for all ~ values
according to Hasofar&ind method is 0.005223. Here again no use is made of the
distribution information given in the problem. Therefore, unless all variables are
normally distributed (which is not the case in this problem), the probability of
failure computed using equation (5.26) will be in error. More accurate values of
the probability of failure will be given for this example fir the next section.
133
w?
5.2.2 Inclusion of Distribution Information:
(5.27 )
(5.28)
and
134
M%)
=* ++7-71
(5.29 )
~, = P{l-l[FX(x*)])
s (530 )
M*)
(531 )
f*(x)
A
\
ACTUAL DISTRIBUTION
w-
\
NORMAL DISTRIBUTION
\
\
X* x
LINEARIZATIONPOINT
135
In our simple beam example discussed earlier, the load follows a Weibull
distribution while the strength was assumed to be normal. An equivalent normal
distribution can be determined for the load using the procedure described above.
The three limit state functions describing the failure of the beam are given by
equations (5.19) to (5.21). II one uses Hasofar/Lind method, then one must
determine the parameters of the equivalent normal distribution in each step of the
iteration procedure for the nonlinear limit state functions (equations 5.20 and 5.21).
The results for the ~ values including the distribution information
are [5.11]:
The results indicate that the ~ values are invariant to the problem
formulation since ~1 = 13z= ~3. They also show that, in this case, inclusion of the
distribution information inmeased the ~ value from 2.5607 to 2.6690 and decreased
the corresponding probabli@ of ftilure from 0.005223 to 0.003802. The values of ~
and ~ deterrnined including the distribution information are more accurate. It
should be emphasized that inclusion of the distribution information does not
always yield larger safety index.
136
5.3 Correlated Random Variables
(5. 32)
where T
is an orthogonal transformation matrix and the superscript
U indicates the transpose. The transformation matrix ~ is such
that
Et c? J=A
[1 [1 (5. 33)
147
11
P~~
12 P13
1
-----
-----
Pln
~2n
1
23
c+ ------- p3n
[11 31 32 1
I -------------- ------- 1
- - --- - 1
1 nl n2 n3
J (5. 34)
where
Cov IXiixjl
P. =
lj a u
x.
1 J (5. 35)
*+ ag ag ag
~z .* .*.
[ ax; ax: ax; 1
* (5. 37)
Et =~~
(538 )
and
(5. 39)
where
[1rx = Q
x
2
0 %=,
u
x
n
139
Px
1
%= P
,X2
I
x
n
i 1
It can be shown (see reference [5.5] ) that the eigen values [~]
are also the variances of the variables YL .
g(x)=ao+ ai Xi
2 (5.40 )
i
*
where ~,, is given by (5..35).
d
The above procedure is valid for transforming a set of
correlated normal variables to a set of uncorrelated normal variables.
If the variables are non-normal, then as an approximation,
equivalent nokmal Yariables can be determined as described
previously under Inclusion of the distribution information. A more
exact procedure would be to use Rosenblatt transformation [5.6]
which requires information on the joint probability density function
of ~. The degree of approximation is illustrated through a numerical
example given in [5.5] and is usually small.
139
M= 9-2
140
w-, .
70
V ---
$0
+--
40[- 63
30+ 30
I 17
i
20- 20
-.
10L 10
+ 91 f--:A+ ~<%--- . 0
!0
)X9
. 20
40
30
44 43
,--- 40
Figure 5.3. Assumed ships route (Figure without the indicated route
was obtained from reference [5. 10] )
141
Table 5.1. General Characteristics of Eighteen Ships.
dwt
SHIP LBP (ft]B (ft) d (ft) Cb (approx. )
TankerNo.1 1076.IM 174.87 81.40 0.86
TankerNo.2 1069.25 163.25 Sa.o.i 0.R3
TankerNo.3 l$KJ.CUJ154.76 60.45 0.83
TankerNo.4 115.99 42.01O.fi
TankerNo.5 754:70 104.46U.40 0.82
TankerNo.6 754.70 105.654474 0.s09
TankerNo.7 7s.69 105.6544.74 O.MM
TankerNo.9 719.10 82.50 39.15 0.786 40; 970
Tuker No.10 620.81 35.72 0.784 31,500
TankerNo.11 .594.WI % E $.$ :.g 26, !ss0
TinkerNo.12* 775.W 105.50 75,500
ol&ep
700.65 9s.4340.70Oani 45,100
Cartiship No. 14 528.50 73.9929.SS 0.61.513,400
Cargo chip
No.15 520.00 75.
W 31,420.573 12,750
CargoshipNo.16 528.
W 76.0029.800.6W5 13,400
B::-cay
Uoo.
oo 106.00 44.550.840 74,200
Bu;ll
ce~
6ii6.2093.8042.630.793 ::,~:
T;nk;r
No.8 693.75 97.0039.170.77.5 , -
i Tanker and mariner ships from referenoeaf~,;~nd[s.~ $mWctivel y
1,000,000
900,000
800,000
i 700,000
4
600,000 TM4KER No. 1
12 KNOTS
500,000
400,000
300,000
6 KNOTS
200,000 i~
100,000
1 I I 1 I I
5 10 15 20 25 30 35 40 45 50
SIGNIFICANT WAVE HEIGHT [FT]
Figure 5.4. Rms of Wave Bending Moment for Three Large Tankers
142
~ 300,00( 6 KNOTS
100,000
So,ooo
.[
> 10 15 20 25 30 35 40 45 50
significant WAVE HEIGHT [FTI
I
. 300,000
OIL-OFX CARRIER MO. 13
c
CARGO SHIp No. 14
2- 15 KNOTS v CARI=O SHIP No. 15
250,000 4 BULK CARRIER No. 17
g :
+BULK CARRIER No. lE
x
+
u
z 200,000 b1
E 13.S KNOTS
z
:
150,000 20 KNOTS
:
+
m 100,000 4
/
i 1S KNOTS
i ; ~~
50,000
o 5 10 15 20 - 25 30 35 40 45 50
SIGNIFICANT wAw HEIGHT [ml
Figure 5.6. Rms of Wave Bending Moment for Cargo and Bulk
Carriers
143
The strength coefficient of variation V$ was investigated next.
A detailed procedure for determining this coefficient is described
and applied to three ships (tanker, cargo, and frigate) in references
[5.7, 5.8]. The strength coefficients of variation of these three ships
were in the range of 7 to 11.3 percent. These figures include the
estimated subjective and the computed objective uncertainties. The
eighteen ships under consideration were assumed to have a strength
coefficient of variation = 13 percent. This is rather pessimistic and is
on the conservative side. The mean values of the strength, m~ were
considered to be equal to the section modulus of the ship multiplied
by the yield strength, taken = 30 k/sq in.
With these parameters determined for all the ships, the safety
index is computed from the equation:
=+
where 0 = ~5/mZ and V5 and V2 are C.O.V. of the strength and load,
respectively.
144
MS = pgL2Bd, which is proportional
to the static bending moment.
Only 15 ships are shown in this figure as the three large tankers 1, 2,
and 3 have very large MS values (470 x 106, 310 X 106, and 267 X
106, respectively) and would change the abscissa scale considerably.
They do however follow the same general trend of the data, which is
higher safety ~ for higher Ms. The physical significance of this figure
is that designing ships on the basis of the static bending moment Ms
would lead, in general, to higher safety for larger ships; that is, the
seems to overestimate the load on large ships.
t.
t
I I
10
I
20
1
10
I
40
I
50
I I I I
60 70 [0 so (x 10$ %
145
Figure 5.9 shows the computed safety index ~ plotted versus
the parameter
(ml) Uy
F~ =
pg Lz Bd
: ~ %
6 1
f
)
t *
s _. :
s
: 1
.4
m
m :
n
I/
n
b
I I I I I I I 1 I I I
1.s 1.1 2.0 2.2 2.4 2.C 2.0 2.0 3.2 ).4 3.6 lx 10-% r,
146
656
bending moment and the strength in the traditional design procedures is a
possible reason. Another reason is the possibility ofinadequacy inthepredition
of the wave bending moment when using the traditional static balance
procedures. A third reason may be due to the slightly different niles of the
different classification societies used for dete rmining the section moduli of the
ships. In addition, the variation of the actual values of the section moduli of the
ships with respect to those specif5ed by the rules is also a contributing factor.
Finally, a more exact procedure for determining the parameters used in this
application may lead to a slightly different range of level of safety.
where
(SM) = minimum section modulus
fy = yield strength of the material
M Sw = stillwater bending moment
MW = wave bending moment
Notice that the product in the first term of (5.41) -makes the limit state
equation nonlinear. The limit state equation in the reduced space is
where the prime superscript indicates a reduced variable and p and c indicate the
mean and standard deviation, respectively.
147
Irl
%rn (lJfY+ fY$) = %m fy (5.43)
Wy(m (5.44)
Ggw (5.45)
-% (5.46)
The direction cosines ~i can then be calculated from (5.17). Using these,
the coordinates of the most likely failure point are calculated from (5.16) and
substitut&in the limit state equation in the reduced space (equation 5.15) to yield
the following resulti
(vsw-aswa*swpb(~w-a~a*wp) = o (5.47)
Equation (5.47) is to be solved for ~. Notice that the U*i are to be evaluated at
the most likely failure point, i.e., fYand (SM) values at that point must be inserted
in equations (5.43) and (5.44), respective y, when evaluating ~*i given by equation
(5. 17). Since these coordinates are unknown apriori, the mean values of these
variables can be used as initial values in an iterative procedure. fir solving for
~ from (5.47), a new set of coordinates can be determined from (5.16) and used ta
determine a new set of ~*i. The procedure is repeated until conve~.nce is
achieved.
The above procedure is accurate if all the random variables in the limit
state function (equation 5.41) are normally distributed. However, as discussed
earlier, the wave bending moment MW follows a Weibull or an exponential
distribution. In this case, the mean and standard deviation of an equivalent
normal distribution must be determined in each iteration step according to
equations (5.30) and (5.31) and used in equation (5.47)prior to solving for ~. This
149
transformation will produce additional non-lineari@ in the limit state equation in
the reduced space.
Numerics 1 Examp~
The values the means (@ and coefficients of variation (S) of the variable in
the Iimit state equation were determined as:
For the first iteration in the procedure described above the mean values
were taken as the most likely failure point. The parameters for the equivalent
normal distribution for the wave bending moment can be calculated directly from
(5.48) and (5.49). The derivatives and timction cosines are then calculated from
equation (5.43 to 5.46) and (5.17), respectively. The resultig direclion cosines are
substituted in (5.47) yielding an equation in ~. The first iteration solution for ~ is
5.95. Using the determined value of ~ a new set of coordinates of the most likely
failure point is determined according to (5.16). The procedure is repeated and the
second iteration result~~a ~ value of 5.04.More iterationsshould be performed until
convergence is achieved. For comparison, if the information on the wave bending
moment distribution was not included, the resulting ~ values are 5.83 and 5.73 for
the first and second iterations, respectively.
150
Ibo
lummcEs
5.2 Hasofer, A.M. and Lind, N. C., Exact and Invariant Second
Moment Format, Journal of the Engineering Mechanics
Division, ASCE, Vol. 100, No. EMI, Feb. 1974, pp. 111-121.
152
6. LEVEL I RELIABILITY 31WLYS1S :
where A; are psf and ~X: are the mean values of the variables. Since
equation (6. 1) represents the failure stmface, 4: ~z: must fall on the
,*
surface, preferably, the most probable failure point, i.e.,
$
x.
Ai=l
Px (6.2)
i
In the normalized variate space (see equation (5. 19 ) we may write
153
[0
Thus the original variate is
* *
x. = +axx?i=~ -af#ux
1 Pxi x.
i 1 i
= ~x(l-a~dvx) (6.3)
i i
Comparing equations (6.2) and (6.3) we conclude that
Ai=l-afpvx
i (6.4)
* *
3. Determine i lx+axxbi~x-m~ fl~.
i i i i
4. Use, the new values of X; from step 3 again in 1 until
convergence is achieved.
154
5* Using obtained in 4 above, reevaluate x ~ - mx
ip. *
P
or isl
n
a. + ai Ai Px =0
1 i
i=l
Because of linearity the partial derivatives ~~~ are independent
of xi, that is i
ag
axti za.o
1
i
M= 8.(X1,X2) = s - z
155
lb~
Z-pz
; z =
Z
and i3M
; az~ = - z
u a
therefore as = s z
i
and az = - 1
2 2 2
a + 0: T + a T
[ s 1 [ a z )
he=l-a~$~~=l-$v~ 1
2
a 2+a -r
[ s z 1
and u
z
AZ=l-UZ$VZ=l+$VZ 1
0 2+#;
s 1
T
[
156
/66
It is the objective of this section to review and assess the
implementation methods and use of reliability analysis in certain
design codes. Specifically, work attributed to or sponsored by the
American Petroleum Institute, the National Bureau of Standards, and
Comite Euro-International du Beton are examined in some detail.
Codes of several other organizations are also discussed and
compared.
R
SF z D + L + W + other load effect= (6.6)
where
157
Presently, nominal loads are all combined with factors of 1.0,
and constant safety factors of 1.67 and 1.25 are used for operating
and extreme loading, respectively. Note that there is a probabilistic
statement implicit in the given safety factors, in that since extreme
events are by nature, rare, the associated safety factor can be
reduced.
where
nominal strength or resistance of component i
partial resistance factor for component i
nominal gravity or dead load effect
load factor- for dead load
L= nominal live load effect
load factor for live load
w= nominal environmental force effects with
prescribed return period (usually 100 years)
3 w= load factor for environmental load
158
Each resistance factor 4 k; is calculated as the product of two
terms representing component strength uncertainty ( +: ) and system
consequence (+5 ), that is
*Ri = i = (6*8)
1.!59
for Prestressing (FIP). The design code is entitled CEB-FIF Model
Code for Concrete Structures [6.9]. Information explaining the use of
probability in the CEB Model Code can be found in references [6.9,
6.10]. It should be mentioned that CE13 ~has been studying structural
reliability for many years and, for this reason, is considered a leader
in the field of such code development.
(6.10)
The format used for defining the design load effects, following
CEB notation, is
(6.11)
where
160
xptfl~y= partial safety coefficients
poL = load combination factor
The CEB Model Code uses, in general, two types of partial safety
factors, x ~ and II , related to the strength and the loads [as in
equation (6. 11)], respectively.
161
design strength = f k /x (6.13)
n
162
The format recommended for use in the A58 Standard is a
combination of the CEB-FIP format described previously, and the
load and resistance factor design format proposed by Ravindra and
Galambos [6.12]. Since the CEB format has already been described ,
The LRFD format of Ravindra and Galambos is explained together
with how the two formats are combined in an optimum way.
Information concerning the LRFD method of Ravindra and Galambos
can be found accompanying reference [6.12] in the September 1978
Journal of the Structural Division of the American Society of Civil
Engineers.
where
resistance factor
u.= nominal resistance
partial safety coefficients
deterministic influence coefficients
mean dead, live, wind loads, etc.
The terms representing the load effects (that is, the right-hand
side of equation (6.14) are defined as follows:
163
The major difference between the CEB and the Ravindra and
Galambos load representations is that the live load is a separate case
with its own load factor in equation (6.14), but is a multiple of the
maximum load ~0~ ~i~ in equation (6.11). NBS believes that the
computational simplification realized by expressing the arbitrary-
point-in-time load as per equation (6.11) will outweigh certain
advantages due to the increased accuracy of having a separate
loading case in equation (6.14).
164
17+
In summary, the format recommended by NBS is
(6.15)
165
llj
where
s{... } = load effects due to all loads in the brackets (that is,
it is not a numerical operator)
TDfm? = load factors
D,L 9**. = are the loads (dead, live, . . .)
v= load combination probability factor equal to 1.0,
0.7, or 0.6 depending upon whether one, two, or
three loads are included within the brackets
m= 1.25 if D acts in the same way as the loads in the
brackets and is 0.85 if D acts in the opposite way.
The load factors account for variations in the load effects due to
model errors and uncertainties in the structural analysis. The v term
reflects the reduced probability of maximum dead, live, wind, etc.,
loads acting simultaneously. Note that if both live and wind loads
were present, equation (6. 16) would design using the entire wind
effect (depending, of course, upon the chosen load factor). The LRFD
procedure and the CEB method are bcth considered more flexible
than the format of equation (6.16).
(6.17)
where
166
Ym = material factor
Rk/ym = design resistance
K= factor depending on type of resistance
S[z F:Y~i} = design loading effect
Fi = characteristic load
Yfi = load factor
167
177
information on this procedure as well as some typical values of the
1oad and resistance factors corresponding to code formats developed
by the various organizations discussed previously. Examples and
additional information can be found in [6,8, 6.6, 6.3, 6.9, 6.11, 6.15].
168
11?
6.1 Ang, A. H.-S. and Tang, W. H., Probability Concepts in
Engineering Planning and Design, Vol. II, John Wiley & Sons,
New York, 1984.
6.3 Moses, F., Guidelines for Calibrating API RP2A for Reliability-
Based Design, API Prac. Project 80-22, American Petroleum
Institute, Oct. 1981.
169
6.9 Common Unified Rules for Different Types of Construction and
Material, Bulletin DInformation No. 124E, Comite Euro-
International du Beton, Paris, April 1978.
6.15 Yura, J., Galambos, T. and Ravindra, M., The Bending Resistance
of Steel Beams, Journal of the Structural Division, ASCE, Vol.
104, No. ST9, Sept. 1978, pp. 1355-1370.
169a
(THIS PAGE INTENTIONALLY LEFT BLANK)
.x
7, SIMULATION AND THE MONTE CARLO METHOD:
The name Monte Carlo method was introduced in 1944 by von Newmann
and Ulam as a code name for their secret work on neutron diffusion problems at
the Los Alamos Laboratory [7.11. The name was chosen apparently because of the
association of the town Monte Carlo with roulette which is one of the simplest
tools that can be used for generating random numbers.
The Monte Carlo method thus consists of the three basic steps:
170
Since the results from the Monte Carlo technique depend on the number of
samples used, they are not exact and are subject to sampling errors. Generally
the accuracy increases as the sample size increases.
The Monte Carlo method is considered now as one of the most powerful
techniques for analyzing complex problems. Since its chief constraint is
computer capability, it is expected to become even more commonly used in the
future as computer capabilities increase and become less expensive to use.
~ .2 fM r 1 li ili nli
171
(7.1)
where N is the tatal number of trials or simulations and n is the number of trials
in which g(~ ~ O.
The ratio nLN is usually very small and the estimated probability of failure
is subjected to considerable uncertainty. In particular the variance of n/N
. depends on the total number of trials N, decreasing as N increases. That is, the
uncertainty in estimating pf decreases as N increases. Statistical rules can be
used to establish the necessary number of trials which depends on the magnitude
of p~ Many variance reduction techniques have been developed to decrease the
variance of n/N with smaller number of trials than would have been necessary
otherwise.
b. In the second method, the probability of failure is estimated by, first fitting
an appropriate probability distribution for g(~ using the trial values described
earlier [7.3]. The moment or any other established statistical method may be used
in the fitting process. Elderton and Johnson [7.4] suggested some distributions
that are suitable for fitting the g(~ data. The probability of failure is then
determined from
where M = g(~) is a random variable representing the margin and fM(m) is its
probability density function as estimated from the fitting process.
7.3 Ge neration o f Rando m Numbers F or a IUmdom Van able With a Presc rib~
(-lontinuou~ Prob abilitv Distribution:
172
fw
specified probability distribution. Therefore, in this section we will first discuss
how to generate uniformly distributed random numbers then how to obtain the
corresponding random numbers with a specified probability distribution.
where a, b and m are nonnegative integers and the quotients xi/m constitute the
sequence of pseudo-random numbers [7.1, 7.21. Such a sequence
repeats itself
after almost m steps, i.e., cyclic. For this reason m must be set very-large e.g; 108
or larger.
173
Suppose that a set of random numbers are to be generated fora random
variable Y which follows a prescribed distribution with C.D.F. I?y(y). The value of
Y at F~y) = x is
Y= FY-l(x) (7.4)
Fx(x) = X O<x<l
This means that if (xl, X2 . . . xn) is a set of numbers of the r.v. X, the
corresponding number obtained from equation (7.4), i.e.,
will have the C.D.F. Fy(y) as required. As an example consider the r.v. Y to have
a Weibull distribution with C.D.F. given by
Yi = k[-lnxi]~? (7.8)
174
The Weibull distribution can be reduced to the Rayleigh and the exponential
distributions as discussed earlier. If the Weibull parameters k and ~ are equal to
~E and 2, respectively, the Weibull distribution reduces to the Rayleign
distribution. If k and 1 are equal to L and 1, it reduces to the exponential
distribution. Thus, substitution for k and! in equation (7.8) w-ill lead to a set of
random numbers for these two special distributions as well.
l-pfl12
error = (7.9)
z-
()
where N is the total number of simulations (sample size) and pf is the probability
of failure. There is a 95% chance that the actual error in the estimated probability
is less than that given by equation (7.9). It is seen that the error is dependent on
the number of simulations N and the probability of failure p~ it decreases by
increasing N or pf. Therefore, if the estimated probability pf is small which is
usually the case, N should be large enough to decrease the error.
There are techniques, however, which may reduce the error (or variance)
without increasing the sample size. These techniques are known as variance
reduction techniques, and the one that is used often in stmctural failure problems
is called Antithetic Variates.
175
estimations Ya = l/2 (Yl + Y2) is also an unbiased estimator since its expected
value E [Ya] is equal to Y. The vafiance ~~a2 of the new estimator Ya is
determined from the individual variances OYIZand OY2Zas:
If Y1 and Y2 are negatively correlated, i.e., the COV. (Yl, Y2) <0, it is seen
from equation (7.10) that the third term becomes negative and
That is, the accuracy of the estimator Ya can be improved (or its variance
can be reduced) if Y1 and Y2 are negatively correlated estimators. The antithetic
variates method is thus a procedure that ensures a negative correlation between
Y 1 and Y2. This can be accomplished in structural reliability problems as
follows.
7.5
~J
71mri
~ of Anal i 1 n im 1 Ev 1 ion f a
Function
176
188
The purpose of this example is to compare numerically simulated results
obtained using a standard random number generator with exact analytical
values. For this purpose, the random process x(i)(t) was examined in ref. [7.6]
where A and o are fixed (deternxinistic) constants and e(i) is a random phase
angle with uniform distribution shown in Figure 7.1.
Analytical Results
The analytical result for the first order and joint probability density
functions (jpdfl of x given that e is uniformly distributed can be derived by
standard statistical methods. These are given by
f(x) = l/ CHJQ2 - x~l -A<x<fq
= 0 otherwise
and
where ~ = ~ - tl. -0
,
The first order pdf of x is shown is figure 7.2. In general, the jpdf is ~ -
difficult to represent graphically. It is three-dimensional with spikes when the
argument of either of the two delta functions is equal to zero. The occurrence of
these spikes will depend on the values of A,o and z as well as the c,urrent values of
x 1 and X2. The factor in front of the delti functions modifies their sum so that the
total area underneath the jpdf will equal one.
Numerical Results
177
/$?
BASIC language was used to create the data file, because it contains a random
number generator. A separate routine was written in FORTRAN to compute the
pdfofeither the random phase angle ortherandomxs. This routine divided the
range from O to 2n for 6, or from -A to A for x, into n intervals of size A6 or Ax.
Then the number of 9s or xs in each intenal was counted. The value of the pdf
at the center of each intenal is given by,
The probability density function of theta compares well with the expected
uniform distribution, see figure 7.3. The pdf of x was computed for the case of
A= 1 magnitude unit
a = 1 rad / (time unit)
t = 1 time unit
N= 100000 samples
The numerically computed points compare extremely well with the analytical
curve, see figure 7.4. There is, however, disagreement between the computed
talues and the analytical curve at the singularities, x = Q.
Second Order PDF: The numerical procedure used to calculate the jpdf was
very similar to that for the first order pdf.
Some easily visualized cases of the jpdf were investigated. Consider the
case when tl = t2 so that w = O. Since sin(0) = O and COS(0)= 1, the argument of
both delta functions in the analytically derived result (equation 7.13) simplify to
(Xl - x2). This implies that the jpdf will consist of spikes along the line, xl = X2.
The numerical results for this case with
A= 1 magnitude unit
N= 100000 samples
Axl, Ax2 = 0.05 magnitude units
178
(fo
is given in figure 7.5.
Another interesting case is when a = 1 rad / (time unit), tl = 1 time unit and
~ = 2.57 time units, then a = tiz. Since sin(n/2) = 1 and cos(ti2) = O, the argument
of both delta functions will equal zero when X12 + X12 = A2. This means that the
jpdf will be a series of spikes along a circle centered at (xl, x2) = (O, O)and of radius
A. The numerically computed jpdf evaluated under this condition is shown in
figure 7.6 and it indeed appears as expected.
The trend shown in figure 7.7, as well as, in figure 7.8, shows that as N increases
for a constant Ax, the numerical values converge to the analytical result for the
pdf. In addition, comparing figure 7.7 to figure 7.8 for a constant N, shows that
decreasing the size of the interval Ax appears to increase the spread of the
fiumerical data. Therefore,. in order to get an accurate numerical representation
of the probability density function, not only the total number of samples and the
size of the increments are important, but mainly their relationship in
determining the number of samples per internal is important.
179
Probability Density Function of Theti
l/2rr I
--- . _..
j
1
(
+ --
0 Theta
I I
i
i
!
!
I i
(
I ) I
1
i I
I 1
I
; I
I I
I I
I I
I I., I
I
+-- L------- - --
1
-[
-~
-A A
~fi)
= .A(sinrt +- Gti))
Figure 7.2. Analytically Derived PDF of X
180
Probability Density Function of Theti
100000 Samples. A9 = 0.01
0.2
1/2n
0.1
3.0
0 in
Tineta
0.2
1. ... ., . .
l/21r
I * -- .-. . .
- -*--
0.1 +
I
I
!
0.0 t . 4
0 Theta 2%
181
Probability Density Function of x
(A, a, t) = (1,1, 1)
100000 %Ples Ax = 0.01
.
1.5
1.0 1
0.5
)
0.0 1 I [
-1.5 1.0
I 1 1
-0.5 0.0 1.0 1.s
Xfi) = Asjn (~ + ~;i))
1.5
1.0
0.s
I
I
0.0 1 I 1~1
. I -1.0 -0.5 0.0 . 1,5
xfi)= ,%in(vt + ~$
182
w
2.5
r
2.0
t
L
1.5
Ill
1.0
.1*OL
183
1
. Probability
1000
4%%- %?-
Ocn~iiy Funct~on
0.0!
of X
20
Probability
10000 $:;/.:
Density
W)-
Function
0.0,
of X
IS
Lo -
I .
.
.
1.3
1.0
os- . 03
.*. . -..-**.
H t --- J -*****-*~
.. ---- ***-**.
04 r T 4
0.0d r I 1
-15 -1,0 -0.s 0.0 03 1.0 13 -{3 -1.0 -0.s O.kl 1.IY 1.s
f
N = 1000 N= 10000
1.5-
. . .
.
.
1.0-
0.0
-13
1
-1.0 -0.5
I .
1
1.0 1.
0.5-
O.b1
-ls
L..J
-1.0
-6.s
*.
I
0.0
.
O.k
1.0
1
1.5
Xo) m A&t + ;0? X~ = hin(~t + Qq
N= 50000 N= 100000
[q6
Probability Density Function of X Probability Dermity Function of X
2.0
1,00 !!%;)*- vi)- ~.,
1
10000 !%;/.; ii)- 0.,
I
1.5i I
. ., .
I
t.o 7
0.5-
.
. j ,y_J
1
0.0 1 t 1 f
-1.s -1.0 -0.s 0.0 lh 1.3 -1.5 -1.0 -0.5 1.0 1.5
X(O _ ~i~(~ +0~0~ Xfi)s Aj~(a~ +~~i~
N = 1000 N= 10000
1.5-
.
1.0-
0.s -
0.0
1 I ( I I
1 0.0
-1.5 1.: r
-1.0 -0.5 0.0 1.0
-13 -1.0 -d.s 0.0 ok 1
1.0 1,1
X@)- fijn(~ +o~(i)) Xm . Asin@ + G(IO
N = 50000 N= 100000
185
177
?,5,2 A~nlication of Monte Ca rlo Method to Reliability of Shi~ Structures
(excerpted from White and Ayyub [7.7]):
In order to compare results of the Monte Carlo simulation method with the
other reliability methods discussed earlier in Chapters 4 and 5 an example
problem is solved using each method. The problem chosen is to determine
probability of ductile yielding of a vessels deck under extreme bending moments.
Any of the other possible modes of failure could have been chosen, for example,
plastic collapse or buckling, but the availability of data on this problem facilitated
comparison of methods. The vessel chosen for the analysis is a naval frigate, the
same one used by Mansour and Faulkner [7,8]. The principal dimensions are
given in Table 7.1 and the midship section is shown in Figure 7.9.
The problem is essentially a simple beam in bending and can be written as:
Mu = CY (7.14)
where Mu is the ultimate bending moment; C is the section modulus of the vessel;
Y is the tensile yields stress of the vessel material.
Z= R-Q (7.15)
where R is the resistance, given in tons/in2 and is equal to Yin Equation (7.14); Q
is the total load in tondinz .
Next a more complicated non-linear form is used. This form separates the
wave and still water bending moments, Mw and M. respectively; and Z is .
expressed in units of bending moment:
z =YC-MO-MW (7.16)
186
Ftl
The basic variables for each form are shown in Table 7.2 along with their
respective statistical properties.
First Order Reliability Method (Level 2): The method described in Chapter 5
was applied in [7.7] to the linear and non-linear limit state functions given by
equations (7. 15) and (7. 16), respectively. In both cases the distribution information
was included in the analysis. The results for the safety index ~ and the
corresponding probability of failure p~ are given in column 2 of Table 7.3.
step 1. Identify the basic variable with the most variability in the limit
state equation.
step 2. Condition the variable in Step 1 with respect to all the remaining
variables in the limit state equation.
step 3. Generate a uniformly distributed random deviate for each of the
conditional variables.
Step 4. Generate a second uniformly distributed random deviate which is
negatively comelated to the one from Step 3.
Step 5. Using the inverse transform method produce a random variable
for each deviate from Step 4.
step 6. Calculate the probability of failure using the probabilistic
characteristics of the variable identified in Step 1 for each set of
random variables.
187
step 7. Find the average probability of failure for the two pfs in Step 6.
Step 8. Repeat Steps 3 ta 7 N times.
step 9. Calculate the statistics of the N number of probabilities of failure
thus generated.
The results are shown in Table 7.3 for 2000 simulations cycles. Figures 7.10
and 7.11 from reference [7.7] show the simulation scheme converges on a solution
with increasing N.
188
200
..
.I
J*
.-
189
!.&,&,.,,,=
0.
0.s
0 12m 1-
190
202
0.1
0.
0.
0
0
~ w Sxunv%w Cvcum
191
Table 7.1. Vessel Characteristics
192
20$
Table 7.3. Example Problem Resul&
I
Linear Limit p = 4.75
Stab Equation (7.15) pf = 0.97x 10+ pf = 1.33 x 10+ M = 0.98X 10~
(C*O.V. = 0.0192)
[7.1] Elishakoff, L, J?robabilistic Methods in the Theom of St~c ture~, John Wiley
and Sons, 1983.
[7.2] Ang, A. H.-S. and Tang, W. H., Probability Concepts in Engineering
Planning and Design, Risk and Reliability, John Wiley and Sons, New
York, 1984.
[7.3] Theft-Christensen, P. and Baker, M., Structural Reliability Theory and Its
Application, Springer-Verlag, Berlin, 1982.
[7.4] Elderton, W. P. and Johnson, N. L., Systems of Frequency Curves,
Cambridge University Press, 1969.
[7.5] Shooman, M. L., Probabilistic Reliability: An Engineering Approach,
McGraw-Hill Book Co., New York, 1968.
[7.6] Ochs, J., Comparison of Analytical and Numerical Evaluation of the
Random Function x(i) = A sin(uz + e), Report for NA 240B, University of
California, Berkeley, May 1988.
[7.7] White, G. J. and Ayyub, B. M., Reliability Methods for Ship Structures,*
Naval Engineers Journal, Vol. 91, No. 4, May 1985.
[7.8] Mansour, A. and Faulkner, D., on Applying the Statistical Approach b
Extreme Sea Loads and Ship Hull Strength, Trans., Royal Institute of
Naval Architects, Vol. 115, 1973.
194
8. SYSTEM RELLM31LI~
8.1. Introduction
Two main sources of system effects are identified. The first is due to
possible multiplicity of failure modes of a component or a structural member. For
example, a beam under bending and axial loads may fail in buckling, flexure or
shear. Each one of theses modes can be defined by one limit state equation. Even
though in this case, we are dealing with a single member (beam), system
reliability methods must be used in order to combine the possible failure modes
and to obtain an assessment of the total risk of failure of the beam. The probability
of failure of one mode may be larger than the others, but the fact that there is a
possibility that the others may occur indicates that they must be included and
combined ta obtain-the tot@ probability of failure of the beam.
The exact system reliability problem taken into consideration possible time-
dependent random variables is an outcrossing problem. If the time-dependant
loads or response of the structure exceeds (outcrosses) one or more of several
possible failure modes (surfaces), failure of the structure occurs. The problem
formulated in terms of stochastic processes however is difficult to solve. Only a
few cases of very simply structures with certain load history models can be
196
Zog
evaluated in this manner and the reliability of the stn.mture at any time during its
life can be calculated. For a single time varying load it is possible to treat the
peaks as a random vaxiable and its extreme-value distribution may be formulated
to perform the reliability calculation.
gi@<O
i=l,2, . ..k
where fX (xl, X2 . . . xn) is the joint probability density function of the n random
variables and gi(~ are the k limit state fuctions. The domain of integration in
, equation (8.1) is over the entire space where each of the k limit state function is
negative or
. zero,
197
It should be noted that, in principle, simulation methods and the Monte
Carlo technique can be used to solve equation (8.1) in basically the same manner
discussed in Chapter 7. In this case numerical simulation of the random
variables is performed according to their prescribed joint distribution and all limit
state equations are checked to see if failure occurs. The ratio of failure
realizations to tital number of simulations gives an estimate of the probability of
system failure. Reduced variate techniques and other methods for improving
convergence may be used here. UsUally, for realistic stmctures the number of
simulations required for a reliable estimate of the system probability of failure is
still high, but these methods have potential for application in system reliability.
.
8,3. Bounds on the Pro babilitv of Failure of a Series SVS~
A series system is one which fails if any one or more of its components
fails. Such a system has no redundancy and is also known as weakest link
system. Schematically a series system is represented as in figure 8.1
----m--m+ 3t---
Figure 8.1. Schematic Representation of a Series System
If Fi denotes the i~ event of failure, i.e., the event that [gi W q 01, ~d Si
represents the corresponding safe event, i.e., [gi ~) a 01, then the combined
19Q
system failure event FS is determined as the union U of dl individual failure
events Fi as
+(y Fi) =
P(Fs) 1- P
()
nsi
i
(8.2)
199
211
On the other hand, if the events of failure are assumed to be statistically
independent, an upper bound (consemative) can be detemnined. In this case, for
independent, failure events of a series system, the right h~d side of equation (8.2)
reduces to
k k
1-P nsi =1- ~ P(Si)=l - II [1- P(Fi)l (8.4)
() i i=l i=l
.
k
where 11P( Si ) represents the product of the probabilities of suwival. The result
i= 1
given by equation (8.4) represents an upper bound on the true probabili~ of system
failure, i.e.,
k
P(Fs)~ 1- II [1- P(Fi)] (8.5)
i=l
Combining equations (8.3) and (8.5), one obtains an upper and lower
bounds, i.e.,
k
max. P(Fi)~P(Fs)~ 1- II [1- P(l?i)l (8.6)
i i=l
Although the upper bound in equation (8.6) is not dMcult to evaluate, it can
be further simplified by noticing that
k k
1- 11 [1- P(Fi)l s X P(Fi) (8.7)
i=l i=l
k
max. P(Fi) ~ P(Fs) ~ z P(Fi) (8.8)
i i=l
200
Equation (8.8) gives the final result for the bounds of a series system and
states the obvious conclusion that the actual probability of system ftiure lies
between the maximum of the individual probabilities and the sum of all
individual probabilities. These bounds are narrow if one mode of failure is
dominant, i.e., if one of the individual probabilities of failure is much larger than
the others. If not, these bounds may be tio wide to be useful. In such cases a
more narrow set of bounds should be considered (second order bounds).
These bounds were developed in references [8.2, 8.3, 8.4 and 8.5] and are
given in terms of pair-wise dependence between failure events, therefore, are
called second order bounds. The original bounds for k potential modes of ftilure
are given as [8.2, 8.31:
k i-1
P(F1) + z max. { [ P(Fi) - E P(Fi Fj) 1 ; O} ~ P(Fs) s
i=2 j=l
k, k-
Z P(Fi) - Z max. P(Fi Fj) (8.9)
i=l i=2 j<i
where P(F1) is the maximum of the individual probabilities of failure ~d P(Fi Fj)
is the probability of inte~section (mutual occurrence) of two events of failure, Fi
~d Fj. P
. .
The bounds given by equation (8.9) depend on the ordering of the failure
modes and different ordering may correspond to wider or narrower bounds.
Therefore, bounds corresponding to different ordering may have to be evaluated to
detmn.ine the narrowest bounds.
The evaluation of the joint probability P(Fi Fj) required in equation (8.9)
remains difficult. A weakened version of the of these bounds (more relaxed
bounds) was proposed by Ditlevsen in [8.4] as follows.
201
P(Fi Fj) s P(A) + P(B) (8.10)
.
whereas, in the upper bound, the same term is replaced by
(8.12)
(8.13)
and W) is the standard normal cumulative distribution fmction and pi are the
individual safety indices (Hasofar-Lind) as discussed in Chapter 5. p is the
correlation coefficient between two failure events (or modes). Such a correlation
coefficient between the failure events (Fi) = (gi ~s O) ~d (Fj) = (gj N ~ O) - be
(8.14)
where
Sg.
[1[1
agj
COV (gi, gj) = mgl#m. 8xrn* (8.15)
and
gi
= [511
[21 (8.16)
(8.17)
202
In equations .(8.14) to (8.17), ~, ~ ., .% are the reduced random variables
and the derivatives are evaluated at the most likely failure points as discussed in
Chapter 5. The proposed bound by Ditlevsen [8.4] apply only for normaliY
distributed random variables.
Narrower bounds than the second order bounds can be constructed, but,
they involve intersection of more than two failure events and are much more
complicated.
A parallel system is one which fails only if all its components fail, i.e.,
failure of one component only will not necessarily constitute ftilure of the system.
Schematically, such a system can be represented as shown in Figure 8.2
IZIEIEI
+
Figure 8.2. Schematic Representation of a Parallel System.
203
Generally, brittle failure implies that the member looses completely its load-
carrying capacity while in ductile failure, the member maintains a certain level
of load-carrying capacity after failure.
If Fi denotes again the i~ event of failure and & the corresponding safe
event, then the system failure event of a parallel system Fp of k components (i.e.,
failure events) is the intersection or mutual occurrence of all failure events, i.e.,
Simple first order lower and upper bounds can be constructed using
similar arguments as for the series system. Now however, perfect correlation
between all failure events (p = 1.0) corresponds to the upper bound and no
correlation between any pair corresponds to the lower bound. Thus, for positively
correlated failure events, these bounds are:
k
II P(Fi) s P(FP) & @n. P(FJ (8.20)
i=l 1
204
some special cases, however, the exact system failure can be evaluated. For
example, Theft-Christensen and Baker (see reference [8.7]) evaluated the
probability of parallel system failure under deterministic loading and other
restrictive conditions.
Fm- ~F31
%
--m-
Figure 8.3. Schematic Representation of Parallel Subsystems
Connected in a Series (Minimal Cut Set)
205
The failure event Fg of a genenl system consisting of parallel subsystems
connected together in a series (minimal cut set) is given by the union (seties) of
intersection (parallel) of individual failure events, i.e.,
Fg=j un
i (Fij) (8.21)
where (Fij) is the iti component failure in the ~ failure path. The probability of
failure of such a system is thus determined from
P(Fg) = P[y n
~ (Fij)
1 (8.22)
The PNET method is based on the fact that perfectly correlated (or, as
approximation, highly correlated) events of failure in a series system have a
system probability given by the lower bound of equation (8.6) or (8.8), whereas,
independent failure events have a system probability given by the upper bound.
Therefore, one may select a threshold value for the correlation coefficient and
206
assumes that failure events with correlation coefficient ~lj above or equal ta the
threshold value to be perfectly correlated, thus can be represented by a
representative event which is the event among them that gives the maximum
probability of failure, say P(Fr) (see lower bound of equation 8.8). M a set of n
failure events (modes) Fi where i = L . . . n is amanged in a decreasing order then
the failure probability of the representative event is P(F1). The remaining, events
with p Ij c p ~ are again rearranged in a decreasing order of their failure
probabilities. Let these be F2, F3, . . . Fk and the pair-wise correlation are p&, p&,
. . . P&- Those ev@nts with p~j >0 are represented by F;. The rem~~ng ones are
rearranged in a decreasing order and the procedure is repeated to search for
other representative events (modes) of failure. The mutual correlation between
the representative events will be low therefore, they may be assumed independent.
Thus, the probability of system failure maybe approximated by (see upper bound
of equation 8.6)
Fault tree analysis finds many applications in the design and operation of
nuclear power plants. It can also be applied to complex stmctural system such as
offshore platforms. In such analysis the fault tree diagram will help in
identifying in a systematic manner the various component failures that form a
failure path and the different possible failure paths that will lead to the top
event, the failure of the entire structure. In addition to identifying all potential
207
failure paths, the fault tree analysis may single out the critical events that
contribute significantly to the likelihood of failure of the structure.
Event tree analysis on the other hand starts with the top event (main
failure) and examines in a logical manner all possible consequences resulting
from the occurrence of such an event (for example loss of life, pollution, explosion,
fire, etc). The consequences of the top event (now called the initiating event) may
or may not be series depending on the possible occurrence of other adverse events
following the initiating event. The identdlcation of all possible subsequences and
scenarios is best accomplished through an event tree diagram. Each path in
the event tree represents a sequence of subsequent events leading ta a particular
consequence. The probability of occurrence of a spetic consequence depends on
the probabilities of the subsequent events and is simply the product of conditional
probabilities of all events along that path.
208
220
are calculated using the equations of motion of the ship if dynamic effects are to be
included otherwise by balancing the vessel on a wave configuration. The loads on
the vessel alternate from hogging which produce compression in the bottom
plating to sagging which induce compression in the deck. This hog/sag variation
must be considered in the hull reliability analysis.
Since hogging and sagging are mutually exclusive events, i.e., the vessel
can be either in hogging or in sagging condition (but not both at the same time),
then the union of the two events given in (8.24) is simply their sum. The
probability of combined event of failure is thus
As mentioned earlier, each of the hogging and sagging conditions will have
several possible modes of failure (or limit states). In each case these modes can be
modelled as a series system since any of them constitute a failure of the hull (or a
limit state to be prevented). Thus bounds on the probability of failure in hogging
condition P(Fh) and in sagging condition P(Fs) can be constructed using equation
(8.8) or (8.9) for first or second order bounds, respectively. The bounds on the
combined probability P(Fc) are simply the sum of the bounds on each condition as
implied by equation (8.25).
209
221
8.8.1 Numen cal Exa mde~
In reference [8.6], the probabilities of fkilure (or reaching a limit stite) were
calculated for a Mariner class vessel in hogging and sagging conditions using
level 3 reliability method. The stillwater moment for the Mariner is a hogging
moment. Several modes were considered in each condition and the results are as
follows:
Notice that these bounds are tight since one mode of failure is dominant
(post-buckling yield of bottom plating). There is no need to consider second order
bounds.
b. Saztinr Condition
210
:
.
.-
211
REFERENCES
8.2 Kounias, E. G., Bounds for the Probability of a Union with Applications,
Annals of Math. Stat., Vol. 39, No. 6,1968.
8.3 Hunter, D., An Upper Bound for the Probability of a Union, Journal of
Applied Probabili@, Vol. 3, No. 3,1976, pp. 597-603.
8.4 Ditlevsen, 0., Narrow Reliability Bounds for Structural Systems, Journal
of Structural Mechanics, Vol. 7, No. 4, 1979, pp. 453-472.
8.6 Mansour, A., Probabilistic Design Concepts in Ship Stmctural Safety and
Reliability, Trans. SNAME, Vol. 80, 1972, pp. 64-97.
~8.7 Theft-Christensen, P. and Baker, M., Structural Reliability Theory and Its
Application, Springer-Verlag, Berlin, 1982.
212
9. FATIGUE RELIABILITY
9.1 INTRODUCTION
Systems represent assemblages of elements, and can range from cargo ships
cargo) causes. The relentless cyclic forces are perhaps one of the most
design stage.
213
X2 !7
specifications and practices. This requires a high level ofQA
system.
The purpose of this discussion has been to outline the major factors and
fail safe design -- fatigue may occur, but the structure can continue
214
22 b
to function until the fatigue symptoms are detected and repairs are made.
occur and no repairs will be necessary. Safe life designs are difficult
structures.
215
9.2 FATIGUE
ANALYSIS
conditions.
structure (system).
(detail) of concern.
-
In development of the following simplified fatigue design procedures
[?.1], itwill be assumed that waves are the source of cyclic forces. It
will be assumed that the long-terni(e.g. T = 100 years) wave height
(9.1)
FH(h) bexp[.(;)k]
216
22t
Hunkane Cornponem
60 -
(g, )=l)
50
.
40 <0
\ \
\ \
30- b,
I \ \
20
O~rational
\
\
lo-
\
\
0 t t 1 i i
1 to 100 103 104 105 106 j(p @
w
The structural detail fatigue stress range (peak to peak) (Sf) will be
s, = CHa (9.3)
N = KS;m (9.4)
Accumulation
... of fatigue damage (D) will beassumed to be described by a
(9.5)
R(Sfi)
D= I1 N(Sfi)
*Footnote:
and, N = K Sf-m
218
TOTAL CYCLIC STRESS RANGE (ksi)
- 0.s 1 2 5 10 20 so Im
11 1[ I I
I I I I
I 111 I 111
o* 1,,,,1 1 1 I 1,,,,1 1 I 1
50 100 200 500 1000 2003 5000
HOT SPOT STRAIN RANGE X II)
CURVES X AND X
219
Where
Fatigue damage (DL) accumulated over the life (T) of the detail can be
( 9.6)
DL = *( YO+YJ
Where
-s (9.7)
. m
No am
Y. = ~~~m(lnNO)fO~
()
1+Eo
- am ( 9.8)
am
N]
Y, = ~H~(lnlJl)~
( )
l+
El
life damage:
220
232
DL (9.9)
D~=
F.,
T. = (F,,)OT (9.10)
where
The fatigue design stress (SfD) will be related to the fatigue design
(9.11)
s ID = CH;D
Thus,
(9.12)
s ID =
T,(YO+YJ
Note:
Use an HfD close to extreme condition design wave height, e.g. 100-year
height.
221
222
For example, for a detail life of 25years and a factor-of-safety of 2
(TS = 50 years), and a wave height exponent (a) of 1.0, the design stress
of 22 ksi (based on API X S-BN curves, Figure 9.2). For an a= 1.5, the
design stress range would increase to 150 ksi. For an SCF = 3, the
nominal stress range would be 50 ksi.
stress ratio, R, that is the ratio of the minimum stress, Smin, to the
maximum stress, Smax for the design wave height, HfD. Thus:
SjD (9 .13)
s PD =
(l-R)
? .
223
9.3 RELIABILITYllODELS
5* S-SCF relationship
Wirsching [9.3, 9,4] has made extensive fatigue reliability studies for
fixed offshore structures. Munse [9.5] has made similar studies for ship
224
TABLE q.1
ATIGUEBEHAVIOR
OF WELDED
JOINTS
Oeffnftfon of f~tlgue faflunfnS-Ndata
Size effect in S-N datS
Effect of ~ld pruffle
. Effect ofcormslon ard cathodic pmtectl~
Assumtlon of llnearmdel md lognomsl distribution for N
Classiffcatlon of jofnt on the bas~s of gemmtm rlther than load Ptttem
. Relationship ktwen stress at jointandstress used to obtain S-N curve
Ignoring possiblestnss endurance in 5-NCUIWC
Cmqatibility of determination of hot SpOt stress with S-N CU~C
WNUFACTURIN$
CONSIDERATIONS
9 F~bri cation uncsrtaintf es
Requi mmnts on wld contmm notmt
DEFINITION OF THE ENVIROWENT
s Use of full scitter diagram of ti~ - To
Variations in TB
Z occummceestimates
Wavedirectionality
Interaction of Waves and Currents
Theo~eti cal mdel used for ocean waves
HYDRODYNAMIC
LOAbSON STRUCTURE
Inertia and Drag Coefficient
Oimctional wave spectra uhich accounts for wave Spreadf ng
Marine growth
Sheltering effects
STRUCTURAL
RESPONSE
TO tiYMODYWIIC LOADS
. Ass@tions mde In spectml analysis
a. 1i near mponse dud ng transfer functf m devehmmnt.
b. 1t nearf zatlon of drag tem
c. at jofnts, 1. No flexibility
f1.ffectofan
iif.centertocenter
cmrdf nates
d. Sof 1 stf ffness in Dynamic Mo41
. D@ ng *f fects f n stntctiral response
f. Dynamic res~se mt accounted for fnanalysis
FATIGuE
STRESSES
ATJOINT
Hsthod of slmis to vtlwte stress concmtrstion
factors
{SCF)
. Paranmtrfc quatfm used for SCF
. Pofnt at lnte~eei~whitefaf
llt~0CCW5
FATIGuEDAMAGE
EWTIONS
. Assmption of Nfners Rule
Am,smtlon of nerrw bsnd @sage equation
fn spectral approach
Ammqtian of Uefbull di stributf on for stress ranges in stress distribution
pproach
WHERCONSIDERATIONS
Definition ~
w
M
M Have Melght Histo;y Spectral Analysis
u-l
1
I Stress Range History
I
* v T
N
Q-d
%7
TABLE %2 ~
FATIGUEDAMAGEAT TIME T,
D = TBmQ/K
n = stress parameter ~
n= fozcisim
i
f. = average frequency of stresses
si = stress range
Gi = fractionof total stress ranges that Si is acting
n= k(m)(247)mr(f+ l)~Yifiuim
accumulated damage (D = A), the S-N curve parameters (K, m, Eq. 9 .4), a
(9.14)
,.fi
Bf2
Pj = P[TST,] (9.15)
P, = @(-p) (9.16)
Index:
ln~/T. (9.17)
P ~,n7
where
p. = median value of T, and
228
T==- ( 9.18)
Bmfl
U,nT = [ln(l+cj)(l+c~)
(l+c~)m] (9.19)
variable.
range model error parameter (B). The errors are attributed to:
Thus,
B = BMB.B~BN.B~ (9.20)
and
for i =M, S, F, N, H.
229
Tables 9.3 and 9.4 summarize the statistical estimates on the B
Munses fatigue reliability analysis [9.5, 9.7] has been based on a two
9*5).
1!
(9.22)
F~(s) = P(S<S) = l-exp
[()]
- ~
where
(9.23)
~= Sm[llliVT]:
(9.24)
F,(s) = l.exp[-(:)flm,]
230
TABLE %3
OF VARIATION OF COMPONENTS OF 6
231
TABLE %4
I Data$et
Daslgnfactom AIB CID E F
[1) (6) (7)
-
S-N ame, inkips m 4*W 4.38 4.42 3*ocr 3.ti 3.00
per squm in~ K 4.6E12 4.6E12 L55F12 5.25E1O 1.29E11 1.46E1O
units . c~ 0.73 O*73 1.35 0.73 1.25 0.67
Rainflowcorrection k 0.79 0.79 0.79 0.86 0.85 0.86
Damage mtio A 1.(M) 1.00 1.00 1*(K) i.oo 1.00
c 0.30 0.30 ,0.30 0.30 0.30 0.30
Stress modeling BA O.w o.7d o*7cf 0.7(Y 0.7CY o.7cf
error c~ 0.17 0.50 0.50 0.50 0.50 0.50
Average frequen~, f., I
in hertz 0.25 0.25 0.25 0.25 0.25 0.25
Safety index implied bys
RP2A design wave peak
stress de (60 .ksi rl.lk)i B 5.34h 12.78 2.09
1
1262
2.57 1.83
Data from Conunentasy of API RF2A, ~. 81, Fig. (2.5.3-2
AWSX data, elastic range only.
T and K joint data providedbymemberof TechnicalAdvisoryCommittee
T and K joint data: an improved vemionof data set D.
VaIues provided by memk of Technical Advisory Committee. Value of C~
now thought to be 10W;
Values provided by memberof TechnicalAclvisory Committee. Numbers are
now considered reasonable for %orst ~ analysis in which wave spreading
and wave directionalityare not mnsidered.
As computed by solving for ~ in Eq. 7.16
Relatively high value due to srmdlvakeof C.. See supersaipt e above.
~~e T~e from UK K)C)ERULES
Note: For a 2@year life, S~ = 53.2 ks~ ~ = 0.69 for 20 year-wave. (Same stmc-
tum would have S~ = lUi ksi and ~ = 0.57 for I(Byear wave climate.)
232
Figure 9.5 [9.3] shows this distributionfor NT = 108 cycles. Figure 9.6
[9.7] shows measured 1ong term, low frequency, wave-induced ship hul1
girder stresses. Shape parameters (~) in the range ofO.7 to 1.3 model
these data well.
The Weibull shape parameter (E) will depend on a large number of factors
Guidelines for~ for platforms in the Gulf of Mexico are g ven in Figure
9.7 [9.8]. Waterline braces and floating marine structures may have ~ in
excess of 1.0. Munse reports [9.5] ES in the range ofO.7 to 1.3 for
hull girder stresses in tankers and cargo ships.
CN = [C;+my:+cyz (9.25)
where
CF = [S:N+C:R] (9.26)
233
&
I I
r
I L I
1//
I
l/vfi
r I I I I I I I I I
Lq
e-
d
E
234
.
4
k
g
a
I&l
> Kl=--ilii / fY .
. .
/4
.
.
//..
/
//.
/4
.
f-
m, o
o
235
.2y7
DEPENDS ON
SHOAUNG DEPENDS ON
1.0 (SHELF PRoFILE ) NAT. PERIOD
VERY NARROW
/
O*5
BROAD
a L
:
F f WATER 1
DEPTH ~ 100 H.
I I 1
w
o! I I 1 1 I I ! ! 1
& +~ MSL -20 -m -1oo -200
u
rm (B] POSITION IN STRUCTURE
FEET BELOW MSL
t DOMINATED
WATER DE~H ~ 100 FT.
o ! t
o I 2 3 4 5
(C) NATIJRAL PERIOO -
SECONDS.
AFTER REF. 9.8
236 Zjq
and
coefficients of variation:
, CSN = 0.62
0 CMR = 0.15
0 Cs = 0.10
0 cc = 0.40
0 CN = 0.96
(9.27)
E(sm) = ;
where R is the mean life, and E(Sm) is the mean or expected value of Sm,
gives the probability of failure (Pf) at the service life (N~) as:
237
~;l.ca
(9.28)
T, E(sm)r(l+c;lOs)
PI =
[
K
1
where
(9*29)
E(Sm) = A(mJSmF ~+ 1
()
Note that the Weibull shape parameter ~ has been approximated as:
E = c;108 (9.30)
N (9.31)
b
= r(:+l)
where
of fatigue data.
23s
9.4 DESIGNAPPLICA710NS
- DETAILS
s fD = S,PR~ (9.32)
where
@ = +(:)1 ( 9.34)
(9.35)
1
z1
(9.36)
1
z
R Fw =
E [ exp(~DulnT )
239
2s--1
where Pfd and BD are the design probability of failure and safety index,
respectively.
be optimized.
[9.7].
Step 3 - Obtain the mean fatigue stress ranges and slope of the S-N
curve based on the type of details and the design number of cycles
Step 4 - Compute the random load factorv from Eq. 9.34 based on NT,
~ and m.
Step 6 - Determine the design fatigue stress range from Eq. 9.32.
24o
-&?
Wirschings [9.3, 9.4] procedure for a fatigue design check proceeds
and then use Figure 9.8 to establish the design stress range,
SpDN = SpD/SCF. :
, *
the fatigue design safety index (P,) [9.3, 9.41. Table 9.5 summarizes
the data and method used to construct the design stress range curves of
Figure 9.8.
241
-2s2
TABLE 9.5
m 4.38 3.00
k 4.6E12 5.25E1O
cK .73 0.73
h .79 .86
2 1.00 1.0
c, .3 .3
B 7 0.7
CB .5 0.5
P, 2.78 2.62
Method for constructing curves: Consider Data Set B for which the target
safety index is B. = 2.78. SR forTs = 100 years is computed directly using
Eq. (a) with~ = 2.78. The 20-year curve we know must pass through the
reference. From Eq. (a), B = 1.63 for the 100-year wave condition. Using
this 8, the 20-year curve is establish from Eq. (a). For#. = 2.78 for a
Ts = 40 years, Eq. (a) fixed SR=52.4 ksi. This value Is scaled to SR = 78.7
ksi for the 100-year wave. Then.corresponding to~ = 0.57, Eq. (a) gives @ =
2.11 for the 100-year wave. This value of~ in Eq. (a) is used to construct
the 40-year curve; a B = 2.11 for 100 years ensures B = 2.78 for 40 years.
Eq. (a) is:
SR(mS,D) = [ln(fo~J]: AK
[ A~OT~Bmexp(pOd,n,)~(~+ 1)
1
242
0
U+a
-
DUO
00
z-
<
a
ma
I
Based on the fatigue analysis process outlined in Eqs. 9.1 through 9.13,
Geyer and Stahl [9.1] have developed a very useful simplified fatigue
feet, and a 20-year service life (T) with a factor of safety of2
(TS = 40-year design fatigue life), they developed the design fatigue
stress range curves shown in Figure g.9 as a function of the stress-wave
formulated the fatigue design as a design damage ratio (Eq. 9.5), d.:
yz (9.37)
AD =
13mexp(@Dalhr)
(9.38)
where
T,f2 (9.39)
Do=
K
244
320
280
200
and Hurricant
160
120
80
40
00 1 1 I I 1
2.4
.4 .8 1.2 1.6 2
Stress-Wave Ht Ex~ent (g)
245
Using reasonable fatigue reference data for calculation ofdo as a
function ofp,, Wirsching and Chen developed Figure 9.10. Then, 8,s
been targeted to a service life, Ts. For the lognormal formulation (Eq.
9.17), the safety index, P,, for any exposure period, (t) can be
expressed in terms of the design safety index, PD for a service life as:
In(t/T.) (9.40)
B, = flD- ~,nT
For t c Ts, the safety index is much larger than B, (Figure 9.11). This
This equation also points out how inspections and repairs might be
utilized to maintain the safety index above some value (Figure 9.12).
Inspections can be used to reduce the uncertainties that contribute to
a,,,(Eq. 9.17), and thus increase ~,. Repairs (if effective and well
done) can increase B, by erasing all or a large portion of the cyclic
damage.
246
Desfgn Lffe, TD, fora Service Life, TS, of 20 Yrs.
400 200 50
*
5 - \ u
v
3 c
u
.
,-
2
.
247
[,
.10 g
I I I I I I I I I I
_ 10-8
,,.
\ 10-7
- *
\ 8 (pl.o _ 10+
[
\ _ 10-6
- Y,*
_l o 4
+.+g
.,~_l 0-3
\ . . .
I I I 1 I I I I I
o 0.2 0.4 0.6 0.8 1.0
ma
EXPOSURE PERIOD/SERVICE LIFE
FIG. 9.11
248
6
I I I
EPAIR +1\ REPAIR +2-
4
,0-4
3
/?-*--:, ----
10-3
[ REPAIR
THRESHOLD
L
MINIMUM
$
2 ,0-2
1 I I I 1 1 I I I I
0246 a 10 12 14 16 18 so
TIME-YEARS
F~G 9.12
24!3
26/
9.5 FATIGUE RELIABILITY - SYST~
system reliability.
P fs [l-(l-Pf,)~] (9.41)
=
or approximately,
(9.42)
P fs = Pf, (9.43)
...
250
242
Correlation expresses the degree of dependence (or independence) between
dependence.
Cov(u,v) (9.44)
puv =
Uu,uv
where U and V are two random variables and the os are the standard
deviations of these variables. The covariance ofU and V, COV(U, V), is:
Cov(u,v)=q(u-pu)(v-pv)]
where the ps are the expected values, E[F],of the variables. The two
components as:
v: (9.45)
P =
v: + v:
where VS2 and VR2 are the squared coefficients of variation of the
(dependence).
251
Figure 9.13 [9.12] shows the Safety Indices for series systems, P,, as a
(B. = 3), the system safety index is approximately equal to the value
based on zero correlation (Eqs. 9.41 and 9.42) forP< 0.8. For high
PfsD (9.46)
P,,~ = ~
Such an approach has been used in developing fatigue design criteria for
(TLP) [9.13].
complex. Martingale and Wirsching [9.14] and Stahl and Geyer [9.15] have
252
Ig
~N+
.
-------
O.g
z
Ulo
=U
ul\
. -------
0.6
/ N=2
/#~ N=5
/
/ N=1O
./ +
0.4
LLiJJ
O* ,.
Ow &= 3.0
i= 0.2 * .
.<
a ~e= 2.0 ---
0.0 I 1 I 1 i 1 1 I I
253
4
3*: 1
3 -A - - .. .A.
2.5
1.5
1 2 3 4 5 7 8
COMPARISON TO PROGRESSIVE
COLLAPSE MODEL &s3.0, p =32%
254
a. First member failure
The effect of adding joints to the parallel elements for the condition of
behaves much more 1ike a series system (Figure 9.15). For example, for
an element correlation of 50 percent, Figure 9.15 indicates that for a
system of 10 elements (joints) in series, the B,tp. ratio is 0.5 and 0.75
number of joints per member, BS/F, = 0.5 to 0.6 and 0.75 to 0.80 for
B, = 2 and 3, respectively.
Given a target reliability for design of the system (B,.), Figures S.13
being developed [9.11, 9.16], and these methods are being extended to
255
P=30 -60%
N =1-4
PARALLEL
100% FAlLED
1.0
I
0.2
Q*O I I I I I I
2 4 6 8 10 12 14 16
FIG. 9.1!5
256
REFERENCES
257
~,11 Nordal, H., Cornell, C. A., and Karamchandani, A., A Structural
System Reliability Care Study of an Eight-leg Steel Jacket
Offshore Production Platform, Marine Structural Reliability
Symposium, Society of Naval Architects and Marine Engineers,
Arlington, Virginia, October 5-6, 1987.
9.18 Madsen, H. O., $kjong, R. K., Tallin, A. G., and Kirkemo, F.,
Probabilistic Fatigue Crack Growth Analysis of Offshore
Structures, with Reliability Updating Through Inspections,
Proceedings of the Marine Structural Reliability Symposium,
Society of Naval Architects and Marine Engineers, Arlington,
Virginia, October 5-6, 1987.
258
.qo
19. APPLICATIONS TO SHIPS AND MARINE STRUCTURES
259
27J
To further amplify on this point, the long-term distribution of
the wave loads acting on a ship may be determined by tracing the
expected route of the ship during its lifetime. Based on ocean wave
statistics along that route, the long-term wave load probability
distribution (usually taken as Weibull or exponential distribution) for
the entire history may be determined. Any lack or deficiency of data
on wave statistics over a period of time covering the ship life should
be corrected for. In the short-term analysis, a distribution of the
extreme load is predicted on the basis of criteria such as one
occumence in a lifetime, one extreme sea storm of a spectilc duration
or a short-term operation in a specific location under severe sea
conditions. For that purpose, one of the extreme wave load
distributions discussed under Prediction of Extreme Wave Loads in
Chapter II is used.
260
272
c* Number of days per year the ship is expected to be
at port and underway
d. Nominal cruising speed and maximum speed and
the corresponding fraction of time during operation.
11 In almost all the main areas where ships operate, statistical data
concerning wave heights and periods have been observed and
tabulated. The surface of the earth is divided into a grid of ten-
degree squares known as Marsden squares. These squares are
arranged into geographic areas over which wave conditions are fairly
uniform. The areas are given a code numbe~ see, for example
reference ~0,2].
261
6. Estimate the strength parameters (mean ~ and standard
deviation u or the coefficient of variation) including the
objective and subjective uncertainties. Each failure mode will
have its corresponding y and 6.
262
4. Calculate the probabilities of failure in each sea condition,
i.e.,
P [ R ~ lfith sea condition ]
263
evaluated. The duration of the presence of the ship in the storm is
limited by the stationarity condition to a short period of time.
Consider now a tanker of length = 763 ft, breadth =125 ft, and
depth = 54.5 ft. We will evaluate the probabilities of exceeding the
wa4Am
..t1b(ww
d4y*
Wd9m.y
, irramant 904 -~,im,
U aOErm,WAq
mAMl
-lm .X)
om, lrwnw d UWr9a I.
- mrawu *ZM
/ )
J
wti~e
--
Aiimmws-=w- 1
1
1A J
still
tmndlw
W.t.r
-*
L [ ,.
ltr*wth
?*b--f,
Figure 10.1.
Iz!@izz
Long-Term Procedure
.E!e
Tha R.M. S. V41U90 Of
wave bending -nt in
diffmrorrt soa conditions
fromstrip theory or X.1 test
I still Watwr
Bondinq mnne
conditional
f ilura in ach
,..
probability of
e= condition
-l
1
---4 Total probability
of failure
26A
initial yield limit state and the ultimate strength limit state in a
storm condition specified by a significant wave height of 29.0 ft and
an average wave period of 10.1 sec. The storm is assumed to be
stationary under these conditions for a period of one hour. The
following parameters were computed for the tanker using a typical
ship motion program:
The mean and standard deviation of the initial yield limit state
were computed to be P.r = 2,420,488 ft-ton ~d ~r = 314,663 ft-ton,
respectively. The corresponding values for the ultimate limit state
were computed to be ~r = 2,804,760 ft-ton and ~r = 392,666 ft-ton.
The mean of the initial yield limit state was simply computed as the
minimum section modulus of the hull amidship multiplied by the
average yield stress. The mean of the ultimate strength was
computed using USAS, an elaborate nonlinear finite-element
program (see 605]).
where FZN (z) is the extreme total bending moment (stillwater plus
extreme wave bending moment). Several distributions can be used
for FZN as discussed earlier in the section entitled Prediction of
Extreme Wave Load in Chapter II. We will use FZ ~ (z) as predicted
from order statistics given by equation (2.54) with Rice distribution
given by equation (2.38) as the initial distribution. In the latter
equation, ms , the stillwater moment is added . Substitution of these
equations in (10.1) and carrying out the integration numerically we
obtain the following values for the initial yield and ultimate strength
limit states:
266
next. The ship motion program computed values for the rms wave
moment ~ , the bandwidth parameter e, and the average period of
5
the wave moment are, respectively, 2.863 X 10 ft-tons, 0.364, and
13.0 sec. It should be noted that the rms value of the moment is on
the conservative side because of the linearity assumption. The
number of wave moment peaks in cme hour is thus 276.9. Based on
these values and the resistance parameters given previously, the
following probabilities of exceeding the limit states where
determined:
pf = 5P~Wfi;~fi=l
i=l i=l
26!3
280
Equation (IO.1) which gives the probability of failure has been
plotted as a function of non-dimensional variables for the case of
FZ ~ (z) estimated from order statistics with Raleigh distribution as the
initial distribution (i.e., Rice distribution with & = O). This
approximation leads to conservative estimates of the probability of
failure as discussed earlier. Figure 103,10.4,10.5,10.6 and 10.7 show the
value of pf as a function of G* = Q , N and v*= ~r - ms .
Fo no
r
As an example of the approximation involved, the probability of
failure computed from these figures for the initial yield limit state
and the second storm condition (significant wave height 38.75 ft.) of
the above example (with e = O) is 1.40x 10-2 compared to pf = 1.23 x
1()-2 obtained earlier. These figures thus will give slightly more
conservative values for pf but eliminate the necessity of numerical
integration.
.-
, .
26$1
10 :. T,,,,.,.,
Tr
:;:j
...,:..
;
i:~:~..;...
I ,.
.:
,.:.::
:.::
.,,,
.,..
. .,
J
1 .,,,..
.i
. . . . . .. . . ,, I.-i...
I.J: ,.,4 . 1, . . . . . .
- . -,.. .
:.. . ,.. .
II
;:::: :::
:.:::::: ::: . . . . .
. . . . . . . . .: .. -. . .
.;. . ,..
.... .
. . .. . . . . . . .. . . . ..-.
. . . . . . . . ,. . . . . .
,.-. ,.
II
.,1. . . . ..4 .
!- .-+.- ------ .. ,. . . . -.
.,...,=
+
4
. .
. .
.
.
. ----- . . ..-
----- ., ---
,
,..
...-
----
.,--- ----
. . ...
,.
4--.-
,-.-
. . .-
::. : -. =-=
,.-. .
!.. _
,..
. .
.
. -
I&
o .-- --
,,-
. . .
.+. - .
,+_ ~ .. -
,+- --dT- -
\ 1
A!&
,
,:
,1;
:= - ==7
--
=
-- :-. _ =
.
-- c
..
-
--
-.
:! -.--i- .T ;, -
S
,,, , :1
--if+ ~
:.) ,:,: 1,: !::!: .!
-
~
t-+ ~ ,, ill
10-7.
:1 :I ~! ~:~ :: I % 1,
.--+ .- n ==7
-
-
--
= ..
- $
: - --
---
..
~= . ...
.
-.
---
. ---- . .. . ~ - .+.
J
.- -
---
=F+l=l= ---
- . .. _..
: ---
,,
4+--
.--.]LJ. -..
2242 ,
m
1o-8- --.. . - ..-
~ ... :, ::
_. --.,. -,.
~ -- .+-- .- . 1 -+-- .,.
+ +-, .. .+ .- 1- L_
7-A
*.. . ., .
--- -4 .- .-*. u___ .4+ ,.- ,. A . ..- ,,
-.4- J .~. *
1= .-4 4+ .. :-.
:, ~:~
-- .: .:+ L+ :: . ... ~; :,4. :;
J
; p:
!;
1.-9
.!,
.,
-
:,;1
7 1
270
- - -
...
. .
::!:
.
.
.
T,,
:,,
. .
:.. .
.
..4. ;., .
. . . . ,.. .
.,..::: :.:: -
::.
.
:
. . . ,. ,., :.. :
.
... .. .,
. . .
... ,. .,: .
,..
.-
,.. . ..
-..
.
.
. .
.
,.,
.-
----
..
:.. ,
-.-.
....-
.
..-
:,
..
---
. .,. , ::
.. . .
,:
::
..
. ..
---
... ,
d
-..
-.. N. 539:.
:;:,:
2
.+.
. .. .40~ .._ ---
0--- . 300
~ 250
-___ =
----.
. .
-
. .._ . -. . . .
--
. -:. ----
----
. ~ ; >0---: ---;-
!.
~
~
=
7
--i-
II 1
! 4
;!:
. =
. ..-+...
.-........ .-----
.+-------._ +... =
. .. . -+. .. ... .. .. . .-
.+- ..... ... . ... . . .. -----
. . . . . ..+_+ .
-
. ----
.. . . ... . . . . .._. . .._ . . . . .-. ,
,.-. . ..+. ... . ,,. ..-< _._j
,.:: -.+,
:.: : .. . .
---
-4. . . .
T
.
; ,.
:,
.+ +:-ii.::::
::;+:
:;::
u *R
..+
-.+-.... .. ;~~jf:::.::
.-. :--
. .... ....
,,
.,-..,.. ,, ...._,., ,,., ,,,
,,
........ ...+.,, ., ....
10
5 9 )
271
-w
10- .,.,::,
::: ::; T
-
..::
T ,,. ....
.!, ,
:.. l! .,,
,.,4
.. . .
,1 . ...! ... ,..
!.. +
.
. ..-
., ., . . . . . .,.4 ... .
,.1 ,.. : 1 ..-
1 :14,..
I
.,4. l.. . . . . .
,.
,!
,, ,,
-
::: :.:::
.,,
- - -
. . . . .. ..
,,. . ,.. .
. .
..4
,.. . .,. .
::: ,.. - ,.-+
;+ -.. .
.,.
. ,.+4 ,.,
., !.,
-.. ,
.- .,. ,..
. . . ,.. -
.
:x ;:
. . . .
10-L ,. .,.
I I
.
I
. . .
1> d ,.. .
s+ .-.- WJ . . . ,..
-..
:::.
--.
- ..
x .- .
. . .. . . .. . . . ... . . . + !. :g
,.. ..- -.
+ Ilv ..- . ,
..m . .-, - .- .
!!i----
a
,... ...
......_
.
i --..-- .------ ---
- --+ ---
,-
---1. - -.
----
-. --
=:
IL
0 .+- .
--- - -
-
1
=
Y
1-4 =&
?- A
~_. = =
.
L ---
!,
--
.
... .
,-.:
, , .
,,
+ --
I ,, y ~ ,1, : l., !,
/1 ~
)
?1 : ,,
,,
10-8 . . . . +. .. -.. . . . . . . . . . .. .-_. - --..
.+.
-.
. .
I .-.
,..
.
-. ---
---- ---- -AL ----- . . . . . . .- .+.
,.
- ...: .,*=
.++.
!_. .
-+,--- .,,
. ..-.
-:--
-.. .
...
. ..
-
. . .
. . .
>
. . J.:.
,I
6
.+
::
.,,, .!--
7
. . . .
8
. ,!
1
9
.,. .
10
L 11
272
11,- -
.
::,
,,, ,,
!,.4
. .
.
.
.
.
.
;.
,,
........
..........
.11 ..4. . . . . .
. ... ~
,,
,.. . . . . . .,. . . . . . . ,.
I
,
,,
- -
::. :
.... ::::. .:;
.1-. ;:1; . .
,,.
::.1..: . . . . . ---
.. .
..-. . . .
. ..-
. ..-
10-4
::::
-. .. .. . ........
----- I
::;:
. . . . . .
. ... . ::::.
..
. ..... .. .. . . .: . --
:;l-
.,. .
.4. .
-----
.--.. . ... . .
----- . . . . . ..-
. . . .+.+-.. .
---- .._....-+ -.
. ..-. .. ----
tt
.-.
-
-.-4 ----
---1
:--1
,.
-.
. -
- 4--
-------
.
---
. ..-.
.
.- ..
- . .
x:
--
..
*
-:
?J$ ... ,
4 - .+ 7-
l-+ -. .
.,,
J-
=
-
,. 2%:
4= ---- e
%<
-
.-,,
-. L-
~-q7 -
y-
-,
-
.. .,
--.. + v=
-- - . .+
.
.
;___
..+.-
--- ~
-..++ Zt=
. ..
..
.
-..
+ -4..
.;
-
:::: . ,..-.
.. .. :: :: ...
----
---
. . .. .. . :::=. .v
.. . . . . . . . .---,
. . . . .- , ..---- ---, .
.. J- .-. .. . .... .... ...... --l
.
_,:i
1*2-.... .. ....
. ;.. ,
.... ...- .. .. .... .... .-.J
*
6 7 8 9 10 11
+ %-m,
P= Jq
273
,, ,. ,:
:
,
.,,::,
,,,
,.. ,,
,..
..,:
1:1/ T,,.+.
,.,
.,
,.,
;.,
.,,
,.
T....,
.,.
,,,
j:
:..1
:::...
. .
;;1,
.
:..
,.,
,,,
.
-
,.:.
.
.,. . :.
,,
,.
,.
. . .
It
;,!!
.. . .,
--_ -
::
-..
.-. .
..
----- . ..-. -.
w ,-. . ---
m . .-..
..-
,-- :
-
--
-.. .
.
---
.-.
d ,,
-
,,
----
u.
--- ---
,.,
-+. .
n
,
m
-+...
-.
o
Z ---
q!
-!19 .
.-. ~..- .
1, ----
Ia!iikE
,< ---- ,-
Eo
-.
_ .
.
-
:
-+
..! .- -
.-
- -..
-
&<
:Io-d : : ,
.+
-----
) : ;
la
T _
m
.-
,,
Zi.
a,, ,, 1<1
---
..
-
>--
.-.
,
.
..
:...
+
------
---
....--, .-
--
I , ; :!, l!;: T ;,,
,1 1, : 11! - -1
,,,1 t
-M!#&!
,!:
--
+.++-.
(,
---- --------
--.+.
%- .-.
.. . .
1==1=+t=l= I
.L
m
. . . . . ..
. . .. .
. ..
&
!:
10-0 -.. ., ---- . . ... . . . .. . . . -,. .
..- . . ..- -- _,-. __ --, -::
-. ..-. . ----- . . . . . . -.
....-;:.
----- . . . .. . ..
_+_r,----
. .. .... ...44:-....+ -- +
....._ .,. ...
.
;..
;,,
..+.
...4....
. .
..-...+
.. . . . .
!-
. .
-:--
.-
J:.-.
.
.~-
....
.... ....,.:.. .....
..+4
.: -4.-
,+
L1.,A
J..
. -----
. .
-
.
10-9
5
.- a;: !? :
4...+
b
.0:
i
7
I I
6
l. +--.+,.
,.
9
!.: .-:,
10
.4..
11
274
10.2.2. -- Level 3 Reliability Based on Four Different Extreme
Distributions.
Short-term Analysis:
N= 276.9 . = 0.35
* q
275
4.00 .
. _________ .______ -- __ . .
1
m
o
-3.00 .. _
.
x
C -t
&
+-
0
r /
/ .
n
n
0 -1
+rnmm
& o.o~O.do mlrlnn nT-n-l-?-r ~
50. 1m 0 K) 150 )0 200:00 25C )0 300. 00
Long-Term Analysis:
278
q o
The root mean square of the WBM ( ~ ) in each sea state are
calculated from standard sea-keeping program and are shown in
Table10 .3. In the same table the scale and location parameters ~Ni
and UNi of the asymptotic distribution as calculated from equations
(2.65) and (2.64) are shown (G = O).
g(~ = R- Ms-Mw
+tn
-UN-(m-M*-uN ) 1
P. l-
L
,lW
I
_@ Xp(-e
L i.-
T ( :
r
)) dm
279
In the second part of this example the SWBM is assumed to be
non-deterministic and correlated to the WBM. Only Level 2 method
w&s used to find ~ (or, pf), where,
1 0
II
1
R=o 1
23
0 P 32 1
00.3)
280
where Yr and Ymw are standard uncomelated variables of R and Mw
respectively and are expressed as,
r-p
r
Y=
r u
r
@4)
- UM(mW - u,)
Y= ? [ exp [ -e )1
mw
@Q5)
281
probability increases by about 25%. Such an increase in the
probability of failure is not considered to be significant and, in fact,
in terms of ~, it would be very small. The results indicate, therefore,
that the correlation between the stillwater and wave bending
moments is not important and may be neglected in future analysis.
282
.
Sea R.M.S. of Scale LOcatlon Exact Failure -
Approx.
State HBl?( m ) Parameter Parameter Probability Failure
(i) (ft-ton) (qj J ( UN;) Probability
pi&xact
p-~lpprox
au 21650 2.6055e-04 1.2212e05 3.2187e-06 1.1515e-g7
01 26700 2.1321e-04 1.5199e05 3.2187e-06 1.9185e-07
02 40415 1.4323e-04 2.3395e05 3.2187e-06 7.4390e-07
03 55570 1.0351e-04 3.1964e05 3.5167e-06 2.8659e-06
04 76500 7.3909e-05 4.3253e05 1.6212e-05 1.5214e-05
05 100315 5.5743e-05 5.6095e05 9.2506e-05 8.7141e-05
06 129905 4.2367e-05 7.1495e05 6.1077e-04 5.7095e-04
07 155885 3.4877e-05 8.4751e05 2.5651e-03 2.3866e-03
08 180420 2.9618e-05 9.6411e05 7.8703e-03 7.3035e-03
09 19918S 2.6752e-05 1.0614e06 1.7956e-02 1.6695e-02
10 221560 2.2548e-05 1.1068e06 2.6627e-02 2.4546e-02
11 241045 2.0779e-05 1.2073e06 5.3888e-02 4.9997e-02
12 262695 1.9320e-05 1.3332et16 1.1368e-01 1.0663e-01
13 278570 1.8063e-05 1.4017eE!6 1.6216e-01 1.5281e-01
294450 1.6522e-05 1.4325e06 1.9008e-01 1.7899e-01
:: 308880 1.5813e-05 1.5086eB6 2.9189e-01 2.4856e-01
16 324040 1.5C80e-05 1.5834e06 3.4412e-01 3.2866e-01
17 338470 1.3811e-05 1.5891e06 3.5483e-01 3.3817e-01
18 350740 1.3379e-05 1.6459e06 4.2291e-01 4.0512e-01
19 363730 1.3268e-OS 1.7553e06 5.5559e-01 4.6154e-01
91 398370 1.1054e-OS 1.7543e06 5.6263e-01 4.5673e-01
92 418580 9.5417e-06 1.6718e06 4.7341e-01 4.5091e-01
94 453940 7.8608e-06 1.6198e06 4.27,87e-01 4.0343e-~1
95 469100 1.5610e-05 3.4350e06 8.9814e-01 8.6712e-01
2.2.
2.L
3
p *lo
f 2.0
1.9.
1.8.
1.7.
1.6.
1.5,
0. .1 3 .5 .7 .9
Jij
285
$
,. 5
..
Linearized -- Exact Failure
Failure Sur~a.c/ . Surface.
Y
u
%N t
I
1 -1
%!! I -1
286
10.2.4. -- Application to Eighteen Ships Using Level 3. M.V.F.O.S.~
and the Improved First Order Methods:
Pf=[l-l[+]]+l[+-%l
+]+
e[-[ $x
@.6)
287
The mean value of the wave bending moment k was
determined by a long term procedure described in the application
example of Chapter 5. Similarly, an initial yield limit state was used
to determine the resistance parameters as described in the same
application example.
288
The lack of consistency in the degree of approximation when
using the M. V. F. O.S.M. or the improved method and the fact that they
lead to optimistic values of ship safety are matters of concern.
Vr Vc
ship t
P P
1 .13 O.aa 5.12 $.11
2 .13 .144 6.W 5.97
1 1 520[-7 1.611
(-7 1.541
C-7
8 $. Sac-? 7.178
[.7 1.784
[-6
9 9.965{-8 1.303
[-7 MU-l
10 5. M5E-7 7.178[:7 2. IW2-4
11 ?.60M-O 2.a52f-7 I.MM-b
12 }.21[-0 9.45[-8 1.wr-7
13 I .795E-a Z.sz[-a 1.1*7[-7
14 3. MC-6 6.212
[-4 1.%0[-$
15 3.732E-6 1,?32[-6 7.01
1[-6
16 2.ME4 5.770[-4 I.
OW[-J
17 I.U-lo 4.OC-10 7.430[-9
)a 1.075
[.0 2.01[-0 I.1mr-1
Table10.6
2 go
no
10
7
9
h 1
2
/
6
Is ,4
9t %
&
4
lb
3 -W
-,4,
I I I 1
Sog 6~ 7m am 1
9m Im Ilm
LEIP(flJ
l-d - 14.P.
Figure lQ 11
P,
10-.
IOti
-
u- -
10- -
m--
10- -
-*
10
m+ -
291
3!3
i)
Table 10.7. Partial Safety Factors.
AP= strength
reduction
fatir.
At= load
~nification
factor.
10.2.5.
Implied in ABS Rules for Ship Longitudinal Strength
292
and Classing Steel Vessels were used to design the ships and the
following three assumptions were made:
Ih = SM~ (107)
293
it is easy to find from standard tables that t = 1.64 satisfies @8) and
using the condition for the COV provides a second equation. Taking
the 300 ft ship as an example, the 9.1 percent COV gives ~s w =
2.316 X 104 ft-ton and Usw = 2.108 X 103 ft-ton.
P = (da-pH-?,
+ & + Ulu)/2 (lo.9)
294
{
p
Equation (LQ9) is the equation for ~ using both the mean-value
first-order second-moment (MVFOSM) as well as the Hasofer/Lind
Method (without fitting a normal distribution), since the limit state
function is linear. Using the values of the parameters for the 300 ft
ship and noticing the pw = Cw = ~, then
P = 4.70 &QIO)
295 -
Before calculating all the partial safety factors, one more
refinement was possible. This refinement is fitting a normal
distribution to the exponential wave load variate at the most likely
failure point. This is relevant since the calculation of the safety
index implicitly assumes that each variable involved has a normal
distribution and thus the safety index can be easily related to a
probability of failure.
(1012)
296
?08
partial safety factor; that is, both lead to the same safety index
values.
297
Table 10,8. Safety Indices of the Series 60 Ships
45
40
ADvANCED METHOD
MVFOSM METHOO
STRENGTH C OV. IO%
STRENGTH C.O.V * 12 %
b -
3.5
3.0
t
I I 1 I I 1 1 I I
2.51
300 400 500 600 700 800 900 I000 1100 I200
L. B.F! (FT)
298
w
TabIe10 .9. Partial Safety Factors and Safety Index of Series 60 Ships
299
10.1 Lartkford, W., The Structural Design of the ASR Catamaran
Cross Structure, Naval Engineers Journal, Aug. 1967.
10.5 Chen, Y, K., Kutt, L. M., Piaszczyk, C.M., and Bieniek, M-P.,
Ultimate Strength of Ship Structures, Trans. SNAME, Vol. 91,
1983.
300
10.10 Soares, C.G. and Moan, T., Statistical Analysis of Stillwater
Bending Moments and Shear Forces in Tankers, Ore and Bulk
Carriers, Proceedings, Third Iberoamerican Congress of Naval
Engineering, Madrid, 1982.
301
11. CONCUOXNG RIMRKS AND RKCM!ENDATIONS
.
302
and statistics. Practitioners and designers are gaining such
familiarity through seminars, symposia and special courses.
Educational institutions are also requiring more probability and
statistics courses to be taken by students at the graduate and
undergraduate levels. This, however, is a slow process that will take
at least one generation in order to produce the necessary infra-
structure for a routine use of reliability and probabilistic methods in
design.
303
In spite of the shortcomings stated above, use of reliability
analysis in design provides advantages and unique features. These
include:
304
2. In the calculations of the reliability indices and probabilities of
failure, the resulting values depend considerably on the
methods used in determining the loads acting on the strut ture
(e.g. extreme versus long-term waves loads) and on the method
of combining these loads. A need exists for standardizing
such procedures for use in design.
305
APPENDIX 1
HELPFUL INFORMATION
. ..
Weibull Tlstributi
AiA~Aq 1
x q-l
p.d.f. = &(x) = (V@
()
~ e-(X/k)~ (1)
where k and ~ are parameters to be determined from data, e.g., data of wave
amplitude or wave-bending moment amplitude. .
The first two moments (mean and variance) of the Weibull distribution are
given by:
306
r(t) = w rm
r (t+l) = t r (t) = t (t-1) . . . (t-r) 17(br)
i.e., generalization of the factorial function.
r(l) = 1
r (n+l) = n! for any n = integer
1
= Tx
z ()
rn+~
() .
1X3X5 . . . (2n-1)
2n T z; n = integer
rn+~
() s
lx3X5...
211/2
(1)l)
r
~ .
$
n=evenhtig~r
9. Exponential Distribution
(5)
From equations (3) and (4), the mean and variance of the exponential
distribution are given by
b
~~~
307
f
37!
here is the mean square value of the process, i.e., the area under the spectral
density of the process. From equations (1) and (2), the Rayleigh distribution is
given by
and, horn equations (3) and (4), its mean and variance are given by
Ei
~~ nfhi ri
a Method of Mometi
The mean and standard deviation of a data sample can be determined from
the usual equations:
308
i~l
ax = Sample slandard deviation =
( L
~
=
(q-Rp
)
*
The resulting values for ~ and ax can be used in conjunction with equations
(3) and (4) to determine the values of k and !, or more conveniently from the ratio:
(13)
which is a function off only. Thus ~ can be estimated from (13) and then inserted
in equation (3) or (4) to determine k.
Fx (X) = 1. ~ -(xik):
therefore,
or,
(1-M=~ogx ogk
1oglOg
Insert
w = log log
(~) ; v ogx
(14)
309
So if (l-FX(X))-l or (l-FX(x)) is plotted against x on log log versus log paper a
straight line is obtained (if the data fits the Weibull distribution closely). The slope
of the straight line is ~ or -1, respectively, and the intercept with the axis is -1log K.
Thus k and ~ can be determined.
kB 4.=1
1= 1,W
310
jw-
More detailed information on these tests and examples of their use in conjunction
with wave bending moment data fkom SL-7 containerships are given in reference
[Al]. A sample of Weibull probability paper which shows the fit of 12319 data
points from SL-7 containerships is shown in Figure Al obtained from reference
[Al].
A1.2. The Safetv of Index V ersus Probab ilitv of Failure for Normal and Othe~
Distributions;
In Chapter 4 and 5 it was shown that, for the simple margin M (or limit
state function) given by
.
w = ~FS(z), fz(z)dz (16)
o
= FG (-~) (17)
Figure A.2 obtained from reference [A.2] shows the relation between pf and
~ for some other distributions of S and Z and spetied values of their coefficients of
variation ( VS= 0.13 and VZ = 0.10). The plot shows that pfis sensitive to the type of
distributions of S and Z in the higher values of ~ (range of low probabilities of
failure).
311
32j
1
S - NORNAL //
v~ = 0.13 z - Exponential
+ NORMAL
v= = 0.10
10-7. s- NORMAL
z- ASYMPTOTIC
EXTREKE II
10-6:_
s.-NOWAL
2- EXTRH4E(n=3) + NO
10-5L
s- k3YMPTOTIC
EXTREME I
z- ASYMPTOTIC
EXTRHfE
10-4~
wvev~~
Nott4k L
10-3=
10-2
B 1 I I
L 2 3 4 5 6
SAFETY INDEX $
312
32$1
REFERENCES
313
Appendix 2: Computer Program CALREL for Performing Reliability
Analysis
Input Description: .
The input to CALREL consists of two parts: i) input data and ii) user
provided subroutines. The data input defines basic physical random
variables, i.e., their mean standard deviation, correlation, etc. and, or
parameters of the optional distribution functions. For level-II
methods only second moment characterization of random variables
are necessary. Both mean value first order second m o m e n t
(MVFOSM) and Hasofer-Lind first order second moment (FOSM)
reliability index can be calculated. For level 111 characterization of
the random variables two options are available: i) first order
marginal distribution method [FOMD, ref .A.1] or ii) first order full
distribution method [FOFD] using Rosenblatt transformation.
314
M
If the basic random variables are independent they can be specified
completely by their respective marginal distributions. For
dependent variables if only marginal distributions and pair wise
correlation among them is known then, a (non-unique) joint
distribution model is implicitly assumed that is consistent with the
specified marginal distributions and correlation structure [FOMD,
ref. All. The marginal distribution function can either be chosen from
the program library or can be specified in an analytical form through
a user specified subroutine called df, The correlation structure for
dependent variables is specified in the input data section.
Hi(xi I Xl, .... xi.1) = p(xi~ xi I Xl = xl, .... xi-l = xi.1) (Al)
If f~~ and F~(~, respectively, represent the joint density and joint
distribution function of & we have:
Xi
ai-1
Fxl, ....Xi(xl ....xi)
= axl .....~Xi-l
(A.2)
Xl, .....li.l (Xl,..., Xi-1)
315
U1 = @-l[Hl(Xl)] = @-l~xl (Xl)]
U3 = @-1[H3(X3&X2)] (A.3)
Output Description
(A.4)
vp(&*)
=~r D-1 (A.5)
317
where ~ = diagonal matrix of standard deviation
(A.6)
References
318
,
.,
.,
. .
User% Guide to
C6LREL - .
-h+
.
A First and Secfind-ctrder Structural Reliability &nalysis Pragram
:.
-. G
.-
INPUT DATA
321
9. Contrcll Flag fer Prclgram Executictn --- Fc4i-mat (15)
1- 5 NEXT Cantrc,l flag fclr program executian
NEXT=O stctp execution .
NEXT=l Restart a brand new analysis
Gll th~ afclrementiclned data
E shc~uld be input After this line.
~Lg
y NEXT=2 Re-analyze the oIU problem
.*. with a different set of parameters
in the performance function.
Only the values clf the parameters
in the perfcirmance function shcluld
be input after this line.
(1) If the nearest puin~ is not fctund in NIT1 steps- the final
status of the analysis will be stored in an unfc~rmatted
rile calrel.sav. This file mu~t remain unaltered if the
analysis is to be continued in an ensuing run.
(2) In a restart prablem? the pragran reads cmly the title and
the contrcll data frclm the input. file. The initial status
c,f the prctbl~m are read from ealrel.Sav*.
In ctrder te be ccmsistent, IFO, NX, NP, NU) N~? and NCORR
must be the same as the previ~us run.
.. - . -..
PF
p\s A 70
x*e
, ..
t
X=*
I ,
P2 = da >0
t3 X78
,..pz
=t=70
pl=u, -
-::
X7V ,
pz=v72
t:
I
. .
.
323
functlcin g(xsdp)
c
c..~..functicm subroutine te c~mpute the limit-state functian
C...*.X = vectar of basic variabl~s
C.i...dp = vectc~r of deterministic parameter%
c
implicit r~al@ (a-hla-~)
dimensictn x(l)~dp(l) .
. .
.
:.
9 E- ;
re**rn +
en~
.
=ubtioutine df(par?x~nd,cdf?pdf)
c
c . . . ..Subrcmt~ne tn cclrnpute pdf and cdf of user-defined distributions
c.....par = vectar ctf parameter distributions
C.....x = value of variable
c.....nd = distribution number (>20)
c.....cdf = computed cdf value
c.....pdf = camputed pdf value
c
implicit realw8 (a-h$a-z)
dimensien par(41
g= to (10?20,.. .) nd-20
10 calf= ...
pdf= ...
return
20 cdf = .. .
pdf= ...
return
.
.
return
end
subrclutine ,hfun(x~ih?ds~hi)
c
c.....Subreutine to compute conditional CDFS far Rasenblatt transformatic!n
c.....>: = v~ctcjr clf basic variables
c.....ih = rc,w number in Ra5Enblatt ~ransfarmaticsn
c.-...ds = v~ctclr cif deterministic distribution parameters
c.....hi = value of ith conditional cdf
c
implicit real+8(a-hsa-z)
dimension x(l),ds(l)
go to (10~20~...) ih
10 hi = ...
return
20 hi% ...
re*rn
*
~. =
.
return
end
324
k41e Usef*+Qr
.
furtctian g(x!dp)
implicit real*8 (a-h to-~)
dimensicm x(l), alp(l)
9 = dp(2)*x (2)** 2-alp(l)*>:(l)
return
end
subrc~utine dftpar?):?nd?cdf!pdf) ,.
implicit real*S (a-hSo-z)
di~n~icm par(4)
-. re.mrn
: end:
.
=ubroutin~ hfun{x?ih?ds~h)
implicit real+=(a-hsa-z)
dimen~ion x(l),ds(l)
return
end
r-
example 1
331 3 c1
2 2 00001
1 20 8 -C).c)ool
1 -11 90.9992 .CM41275 0. 0.
2 2 at>. !5. (). c1
1.0 0.!5
c1
0.
325
exaample 1
type clf first-clrder technique used .. ... ..ifa= 3
ifca=i ...mean-value! 1st order, 2nd moment (mvfclsm) method
ifc,=~ ...first-c,rder? seccind mc~ment (fctsm) methc~d
ifc,=3 . ..first-arder. marginal distribution (fclmd) methad
ifct=4 ...first-c.rderl full dist.(flasenblatt trans) method
type af secclnd-nrder technique used ... ..is~= 3
isct=~ ........ ... ..nci secc~nd-arder approximation
lsu=l . . . . . . . . . . . . . . . . . . . . . . peint fitting method
1sc~=2 . . . . . . . . . . . . . . . . ..curvature fitting method
i ~nz~ . . . . . ..point and curvature fitting methods
type c~f integration s.cfiemesused in seccmd-arder
analysls . ......... ... ..... ..... . ..... . . ltg= 1
itg.eq.C~ ....Hreitung? Tvedts 2-term and single integral
itg.ne.tj -.-all c.f th= abclve plus Tvedtss dauble integral
ty~e of sensitivity analysis required ...isv= 3
isv=Cl . ............... ..... net ~ensitivity analysis required
lsv=l .. ....... ... ..senzit ivity clf distribution parameters
i5v=2 . .. ....sensitivity cff perfclmance function parameters
isv=s . . ..dlStrlhltiCttT and perfclrmance functian parameters
number ctf randctm variables .... .. ..... ....nx= 2
number Clf dc~terministicparameters ... . ...np= 2
num5er clf us~r provided distributicm . ....nu= 0
mumtler Of pdramet=r= 1~ u=~r-defined full
di~tributians ....... .. .......... ..... ....ns= c}
ccirrelaticcn structure . . . . . . . . . . . ... . . . .ncc~rr= (1
nCctrr.eq.O . . . . . . . . . . . . . ..uncc!rrelated variables
ncctrr.n_e.0 . . . . . . . . . . . . . . . . . carrelat~d variablras
. .
initlall~~~tlctn flag ............. .. ... ..init= 6
init.eq10 ... .. ... =.initialization at mean paint
init.me.~ ...... .... ......init ialization by user
.
Clutput flag . .* . . . . . . . . . . . . . . . . . . . . .m. . . ipr= 1
lpr.eq.O . . . . . . . . . . . . . .cutput anly final results
ipr.ne.~ . ....= . .... ...autput at every ipr ~teps
.-.
326
Icop T=l . . . . . . . . . . . . . . . . . . . . . . . . . . . ..nr-flkli,~ *,iQu
,iclpt=2 ....................Mc.dified
RF-HL method
iopt=3 ...............Gradient Prejecticln method
maximum number ef iteratian cycles .....nitl= 20
ma%imum steps In lzne search .... ...... .ni t2= 4
ccfnverqence. tolerance . . . . . . . . . . . . ..tal= 1.CIOOE-C14
optimization parameter 1 ...... .. ..ctptl= 1.000E+OO
optimi=atlon parameter 2 ........ ..ctpt2= O*OCIOE-01
ctptimi:atlon parameter 3 . . . . . . . . ..apt3= C}.CIOOE-01
327
iteratic,n nc~. 4
var. lineari2atic1n paint unit narmal
x Y alpha
L 1.038E+02 4.636E-01 c1.3777
2 1.4E&E+01 -1.17CJE+C10 -CI.9239
reliability index beta = 1. 259!5
------------------------- -------------------- --
iteratictn no. 5
var. linearization unit nurmal paint ..
Y alpha
.
~~.c)A;E+02 4.7SSE-01 0.3781
: &4S6E+01 -1.166E+O0 -CI.9258 +
reliabil~~ index beta = 1. 2S89
----------------------- ------------------------- -------------------------- ---
iteratic~n no. &
var design puint unit sensitivity vectsrs nearest pt unit ficmmal
X* gamma delta eta y* alpha
1 lo040E+02 0.3702 -i~.271s -().C)841 4.761E-01 0.3782
2 1.4S&E+01 -cl.9257 c).9624 -0. ?965 -1.165E+CJCl -cl.9257
----------------------- ---------------------- -------------------------- ------
reliability index beta = 1. 23g9
failu~e probability = 1.04CIE-01
***************+*****************+********************
*
Sensitivity Analysis
* *
***************+*****************+********************
d(beta)/d(parameter) :
var mean ~td dev par 1 par 2 par 3 par 4
1 -1.773E-C12 -5.349E-03 -1.773E-02 4.157E+CM3
2 2.514E-01 -2.536E-01 3.7SJE+CKI -4.3S2E+CKI
d(pfl)fd (parameter) :
var mean !atd dev par 1 par ~ par 3 par 4
1 3.2C12E-03 9.661E-04 3.202E-CJ3 -7.5C17E-01
2 -4.54CIE-02 4.58CIE-02 -6.791E-01 7.914E-CI1
. . . .. .
*---
328
********+*********************************************
* *
* Second Order Structural Reliability Analysis *
* *
Point Fitting lleth~d
* *
******************************************************
. .
:.+
++ coordinates
-1
? and ave. main curvature= cif fitting pctints im ratated space +*+
+
Y 1 = ;.259944 yz= 1.216377
Y*1= -1.250944 y2= 1.e1750b
a 1 =-c}.2&47&lSE-01
------------------------- -.-_---&------------- -----_--_---------
failure generalized
probability reliability index
Breitung a=ymptcltic fc~rmula 1.077E-01 1.2390
Tvedt tbre~ term fclrmula 1.091E-01 1.2315
Iwedt single integral formula l.CW.IE-01 1.2316
Tvedt dcmble int~gral formula 1.090E-01 1.2316
**********+***************+ ***************************
*
Secand Order Structural Reliability Gnalysis *
* #
* Curvature Fitting Method *
*
*****+*******+**+*************************************
secclnd-arder apprnximatian **
failure generalized
probability reliability index
Breitung a5ymptGtic fctrmula . 1.078E-01 1.2303
Tvedt three term fc,rmula 1.092E-01 1.2306
Tve-dt single integral formula 1.C)92E-01 1.2306
Tvedt dcluble
.7= integral formula 1.092E-01 1.2306
~ :b
s
-- ~
a.
. . - - -.. -
329
4
examplE ~
type af first-order technique used .. .....ifo= 3
iffi=l ...rnean-valueq 1st carder, 2nd mcunent (mvfosm) methctd
ifc,=z .~.firgt-c{rder, secc,nd mament (fosm) method
ifo=3 ...first-cjrder~ marginal distribution (fomd) method
ifc,=4 . ..fir~t-order$ full dist. (Rc,senblatt trans) methed
type clf ~ecc~nd-order technique used .....iso= 3
i%a=(} ....... ......no second-arder appreximatien
isci=l . .... =.. ... ... ........paint fitting methed
iso=2 ....- .. ...........curvature fitting methed
ise=3 . . . . . ..paint and curvature fitting m=thed~
type c,f integrdtien schemes u~~d in seccmd-firder
analysis ..... .. . ..... .. . ......*1tg=
1
330
aptimizatic~n ~cheme used . . . . . . . . . . . . . ..lCipt= 1 4
icapt=l ............................ .RF-HL method
ictpt=~ .......... ..... .....Mcldified RF-HL mettmd
, ic,pt=3 . . . . . . . . . . . . . ..Gredient %ojecti~n method
maximum number Of iteratiun cycl~s . . . ..nitl= 20
ma>.lmurn =teps In line search ... .. ... .. .ni t2= 4
cctnvergsnc~ tal=rance . . . . . . . . . . . . ..tc.l= 1.000E-C)4
c~ptimizatlcm parameter 1 ..........optl= 1.000E+OO
Ciptimization parameter 2 ....... ...c*pt2= O.OCKIE-01
c~ptimizaticm parameter 3 . . . . ... . ..c.pt3= 0.000E-01
.-_:
availabl~ prabatility distributions:
,=
normal,
. .... ..............ndis=l
lngnarmal . . . . . . . . . . . . . . ..ndis=2
gamma . . . . . . . . . . . . . . . . . . ..ndi~=q
e::panential . . . . . . . . . . . . ..ndis=b
rayleigh . . . . . . . . . . . . . . . . ..n d is=S
unlfctrm . . . . . . . . . . . . . . . . . . .ndis=6
beta .......... ...... ......ndis=7
type i largest value .....ndis=ll
tfp~ i =mallest value ....ndiS=12
type ii lar~est value .. ..ndis=13
welbull . ..................ndis=14
. .- ~.
**** suluticm phase ****
*********~+**+**+*******++*****+4****+****************
~~1=
1.552342 Y2= 1.4&3146
y?l= -1 .!552B42 Y2= 1.43.5239
a i =--lJ.3d4sga3E-cll
------------------ ________ ______ ___ ______ _______ ______ ___
.-
332
function g(x~dp)
implicit real*S (a-hqa-~)
1 dimensicln x(l) !dp(l)
9 = alp(l) - dp(2)*x(l)
return
- dp(3)*x(2) ~ = 5-x, -xz
end
.
subrciutine df(par~x!nd~cdf?pdf)
implicit real*8 (a-hza-z)
difi~nsiun par(4)
rekrn
end,
subrcmtin~ hfun(x,ih?d~th)
implicit real+S(a-h9e-z)
dimensien x(l),ds(l)
ga to (10!2C1) ih
10h= l-dexp(-x(l))
r~turn
2~h= l-( l+ds(l)+x (s))+de:~p(-x (2)-d~( l)*~(l)*::(2))
return
end
233
+****+*******+***+*****+**s.**********+*******+**+***********
+ Un i v~r s i t y c, f Ca 1 i for n i a *
* Department cl
f Eivll Engineering *
* Division clf
+ Structural Engineering and Structural Mechanics *
*
* CALREL
* First and Second ~rder Reliability Analysis Prc~gram *
*
- *
*,.> Develc.p~d By . *
* Pei-Ling Liu and Armen Der Kiureghian *
*. *
Last Revision: September 1986
* *
* Extended fcrr F!3FD and SOSRAP by *
* HONG-ZDNG LIh!r Ilarch 19S6 *
**+**+*******+* ***+***************+*+***+*****++***+**++****
example 3
type cIf first-order tschnique used .. .....ifo= 4
i%c=l ...cIeari-value~ Ist cir-derz 2nd mclment (mvfclsm) m~tticid
ifc=2 . . .fir~t-c!rd~r$ s=cc~nd Mcim=nt (fosal) ,methctd
if~=3 ...fir~t-c~:-der~ marginal distributic!n (fc~ma) methctd
i~c~=b . . .flrst-cli-der> fuil dist. (Rcts~rhlatt tran=) methctd
type Cif SFcct.-lti -c,i-d~r tgchnique used ... ..isc= 3
isc~=O ............ .nc! secclnd-or~sr aFprciximatlan
ism=l .. . . . . . . . . . . . . . . . . . . ..paint fitting methctd
iscl=~ . . . . . . . . . . . . . . . . ..curvatur= fitting method
i=cl=~ . . . . . ..pc~lnt and curvature fitting metheds
type of ir,tegratic:n schemes used in secc!nd-~rder
dnalysls . ....*.. ... .. .. ....... . ..... .-. ltg= 1
]tg.r.q.(} . . . .H.-eitL~rlgsTvedtss 3-term and slnqle integral
i,tg.rl-.0 . . ll clf the abcrve plus Tv~dt~ C!cluble integral
type Clf %Snsitivit y ana.ly~i~. req~.;ired ...isv= 3
i~.V,=<) .. ........ . .. .....nci serlsitivity analysis r~quirec!
334
..
.
ic,pt=l ......... ............... . .. ..RF-HL methcld
ic,pt=2 . . . . . . . . . . . . . . . . . . . .flc,dified RF-HL methc,d
sc!pt=3 . . . . . . . . . . . . . ..(%adient Prajecticin method
q!imum number c[f iters,ticln cycles . . . ..nitl= 20
rna::imum steps in lirle search --.*... . ...ni t2= 4
ccnvergemce tcrler=nce . . . . . . . . . . . . . tc,l= 1.00CIE-134
.
c+ptlml:atlc~n parameter 1 . . . . . . . . . . c~ptl= 1. [:@t:~E+C)~)
ctptimizaticln p~ram=t~r 2 . . . . . . . . . . clpt~= fJ.~c~~)E-01
c,ptirr,izatic,n param~ter 3 . . . . . . . . ..r~pt3= O.OW)E-01
------------------------- ----------------------
iteraticln no. 1
var . linei~rizatictn paint unit nc~rmal
>: alpha
1 1.82SE-C)$ -9.&E-CJl -0.42!54
2 4.R19E+CK~ z.t)b4E+~cI c1.9135CI
r~liahjlty index keta = 2. 27s8
it~raticln rlcl. 2
var. linearizatic,n pc,int unit nc,rmal
.~ Y alpha
1 1.S17E-01 -9.65E-Cll -(2.4246
2 4.S18E+OC1 2.{}A2E+O(} Cl.9254
r~liabi.lty ind~x beta = 2. 2?E@
--------------------- ------------------------ --
iteraticin nci. 3
var. linearizaticin pctint unit nclr.mal
.A alpha
Y
1 1.823E-01 -9.475E-<11 -0. 42!52
2 4.G15E+O0 Z.(:163E+C)0 CI.9CE1
reliability index bcita = 2. Z78E!
iteraticlr m. 4
Var . linearizaticlri pclint unit normal
X Y alpha
i 1.SISE-01 -9.690E-c}I -CJ .4247
E 4.3:SE+O0 ~.(:)L.zE+{:)<i 0.?0!53
reliabzlty index beta = 2 27a8
------------------------- ----------------------
iteraticln ncl. 5
var . linearization pclint unit nurmal
:.: Y alpha
1 1.S%2E-01 -9.679E-111 -CI.42Eil
2 4.91BE+CKI 2.cI&3E+C)CI 0.9<)51
reliability inde:: b@*a = 2.2700
---------------------- -------------------------
iteratictn nc,. &
var. lirlearlzaticln pclint unit nc!r.mal
-- alpha
Y .
1 ti;91;E-1:)1 -9.68?E-01 -0.424S *.
2 4.81BE+OCI 2.C)63E+CM3 ().9()53
r~liabilty inde:: beta = 2. 27S8
-------=- --------------------- ----------------- . -. .
itfars.tic~n m,. 7
var . llrl~arizatinn
pciint unit nclrmal
X Y alpha
1 1.E)21E-01 -9.A91E-(:J1 -{1.42Z0
2 4.PISE+CWJ 2.C16SE+C}0 0. 90S2
reliability index beta = 2.2708
---- - . ----------- -------
--?
var. L1l-1e-1 liaLAL*Il p~.~i,+ -. .-.,.+ ,. -----
4
alpha
x
******************+*************************+*********
* *
* Sensitivity Analysis *
**+**************+************************+*******+***
.. . . .
336
~, j = 2.273824 ~.. ~ = 2.545[:)3<)
~s~= -Z.27W324 ~F~= 2.5%9305
s 1 = ~).=f:JL@2&pE-~~l
------ ------ ----- ----- ------ ----- ----------- _____ ----- _____ _____
337
Ta run ChLREL, da the following:
(5) Edit the three subrcwtines in u==r.far SUCh that they work
in the same manner as indicated in the cla=s handaut.
Tcc edit the subroutines, key in:
CC:\CE2493> FE USER.FOR
TO display the help file in PEJ press <Fl> key.
(9) Read c,utput file (if the output iS raut~d &a a file)
KC:\CE2491} TYPE OUT (read the file OUT on
screen)
tC:\CE2491> NETFRINT OUT (get a hardeapy of the output file)
Befclre the print queue is mbmitted~ make sure that the line
printer is enline.
..-. . , .. ....... . .
(10) Stare yaur files in a flctppy disk and delete ll files in
-.
...
,..
. the working direet~ry.
*.=
. [C:\CE249J> COPY (yaur file) A:
C
G
~- EC:\CE2491> ERASE *.* ~
&
A
(1$) Delete the temporary directery. *
CC:\CE2493> CD \
tC:\J> RD CE249
[12) la shut dawn the machine, take yaur disk~tt~ Gut and turn
...
C3
SHIP STRUCTURE COMMITTEE PUBLICATIONS