Reliability Based Design Optimization
Reliability Based Design Optimization
Reliability Based Design Optimization
Optimization
Outline
RBDO formulation
min F (d,p) objective
fi (d,p) > 0 soft constraints
P (gj (d,p ) > 0) > Pt hard constraints
Comparison b/w RBDO and Deterministic
Optimization
Feasible Region
Deterministic Optimum
Basic reliability problem
Reliabilty Calculation
Probability of failure
Reliability Index
Reliability index
Formulation of structural reliability problem
Vector of basic random variables
Safe domain
f X ( x ) = const. f
Failure domain
g(x) = 0
s
Limit state surface
0 x1
Geometrical interpretation
uS
safe domain
S
0 uR
Distance from the origin [uR, uS] to the Cornell reliability index
linear limit state surface
Hasofer-Lind reliability index
• Lack of invariance, characteristic for the Cornell reliability
index, can be resolved by expanding the Taylor series
around a point on the limit state surface. Since alternative
formulation of the limit state function correspond to the same
surface, the linearization remains invariant of the formulation.
• The point chosen for the linearization is one which has the
minimum distance from the origin in the space of transformed
standard random variables . The point is known as the
design point or most probable point since it has the highest
likelihood among all points in the failure domain.
Geometrical interpretation
For the linear limit state function, the absolute value of the
reliability index, defined as , is equal to the distance
from the origin of the space (standard normal space) to
the limit state surface.
Hasofer-Lind reliability index
u2
failure domain
f
safe domain u*
S
u1
0 G3 ( u ) = 0
G2 ( u ) = 0
G1( u ) = 0
RBDO formulations
RBDO Methods
Objective
function
Reliability Reliability
Evaluation Evaluation
For 1st For mth
constraint constraint
Decoupled method(SORA)
Deterministic optimization loop
Objective function : min F(d,µx)
Subject to : f(d,µx) < 0
g(d,p,µx-si,) < 0
k = k+1
si = µxk – xkmpp dk,µxk
xkmpp ,pmpp
inactive
constraint
PMA MPP
PMA MPP
RIA MPP
U Space
Methods for reliability computation
Numerical computation of the integral in definition for
large number of random variables (n > 5) is extremely difficult
or even impossible. In practice, for the probability of failure
assessment the following methods are employed:
f
n ( u,0,I ) = const
u* region of most
contribution to
probability integral
*
G (u) = 0
0 u1
l(u) = 0
s
SORM – Second Order Reliability Method
u2 v n v n
f
~
v
v*
Gv(v) = f v ( ~
v ) – vn = 0
v n = sv ( ~
v)
0 u1
s
Gradient Based Method for finding MPP
find α = -grad(uk)/||grad(uk)||
uk+1=βt* α
If |uk+1-uk|<ε, stop
uk+1 is the mpp point
else goto start
If g(uk+1)>g(uk), then perform an arc search
which is a uni-directional optimization
Abdo-Rackwitz-Fiessler algorithm
find
subject to
x1 β1
Reliability based
β2
optimum
x2
Scope for Future Research