$$$$ FIABILI NOTICE CHAT 11/09/12 21:16:11 7124
DATE 11/09/12
Procedure FIABILI Voir aussi : EXCENATAF
----------------- FDENSREPARTFINVREPA
FIAB TAB1 ;
TAB1 . param_optimisation . methode
. param_optimisation . t0
. param_optimisation . s0
. param_optimisation . vxmin
. param_optimisation . vxmax
. param_optimisation . vcmax
. param_optimisation . maxiteration
. noms_des_variables
. max_iteration
. fct_limite
. grad_fct_limite
. param_va . k . typva
. param_va . k . A
. param_va . k . B
. param_va . k . LAMBDA
. param_va . k . MU
. param_va . k . MOYENNE
. param_va . k . ECART_TYPE
. param_va . k . TAU
. param_va . k . K
. param_va . k . W
. param_va . k . MIN
. param_va . k . MAX
. param_va . k . U
. matcov
. points_initiaux
. critere
. resu . i . indfiab
. resu . i . <>
. resu . i . proba_defaillance
. resu . i . facteurs_de_sensibilite
. resu . i . vecteurs_des_sensibilites
Description :
_____________
The procedure FIABILI search the failure probability of a structure.
This probability is estmated by the FORM method. The procedure give
the Hasofer-Lind reliability index. The procedure give the sensibility
of each random variable.
Contents :
__________
TAB1 . 'PARAM_OPTIMISATION' : is a table that hold the optimisation's
method parameter. We use the operator EXCE of castem 2000 so you have to
see the notice of this operator for more details. Attention the value
are given here in listreels.
TAB1 . 'PARAM_OPTIMISATION' . 'METHODE' : (optional) is a MOT. The choice
is between 'STA', 'MOV', 'LIN'.
TAB1 . 'PARAM_OPTIMISATION' . 'T0' : (optional) is a real between 0 and 1.
TAB1 . 'PARAM_OPTIMISATION' . 'S0' : (optional) is a real between 0 and 1.
TAB1 . 'PARAM_OPTIMISATION' . 'VXMIN' : listreels of minimal values that
can take the random variables.
TAB1 . 'PARAM_OPTIMISATION' . 'VXMAX' : listreel of maximal values that
can take the random variables.
TAB1 . 'PARAM_OPTIMISATION' . 'VCMAX' : listreel of the Cjmax.
TAB1 . 'PARAM_OPTIMISATION' . 'MAXITERATION' :(optional) integer.
Change the maximal value of the iteration number in EXCE.
TAB1 . 'NOMS_DES_VARIABLES' : listmots containing the name of each
random variables.
TAB1 . 'MAX_ITERATION' : (optional) changes the maximum value
permitted for the number of times we use EXCE.
TAB1 . 'FCT_LIMITE' : table ranging from one to number of limit function,
which contains the character chains whose are the name of the procedure
computing the limit functions.
These procedures have in input a listreel of values taken by the random
value. they give as output a real which is the value of the limit
function in that point.
TAB1 . 'GRAD_FCT_LIMITE' :(optional) table ranging from one to number
of limit
function, which contains the character chains whose are the name of the
procedure computing the limit functions gradient.
These procedures have in input a listreel of values taken by the random
value. they give as output a real which is the value of the limit
function 's gradient in that point. This is optional. Don't give it,
you'll do less call to the limit function.
TAB1 . 'PARAM_VA' : table ranging from one to random variables number.
TAB1 . 'PARAM_VA' . k : table containing the parameters needed for the
kth random variables.
TAB1 . 'PARAM_VA' . k . 'TYPVA' : chain of character containing the type
of the kth variable.
The types are :
'LOI_UNIFORME'
'LOI_DE_LAPLACE'
'LOI_NORMALE_STANDARD' (i.e. centree,reduite)
'LOI_EXPONENTIELLE'
'LOI_LOGNORMALE'
'LOI_NORMALE'
'LOI_WEIBULL_MIN'
'LOI_NORMALE_TRONQUEE'
'LOI_EXPONENTIELLE_TRONQUEE'
'LOI_GUMBEL_MAX'
'LOI_NORMALE_TRONQUEE_INF'
'LOI_DE_FRECHET'
In the uniform distribution case :
TAB1 . 'PARAM_VA' . k . 'A'
TAB1 . 'PARAM_VA' . k . 'B' : are the real between wich range the kth
random variable. (A= mu
0 else
In the case of the lognormal distribution :
TAB1 . 'PARAM_VA' . k . 'MOYENNE'
TAB1 . 'PARAM_VA' . k . 'ECART_TYPE'
are the mean value and the standard deviation of the kth random variable.
In the case of the normal distribution :
TAB1 . 'PARAM_VA' . k . 'MOYENNE'
TAB1 . 'PARAM_VA' . k . 'ECART_TYPE'
are the mean value and the standard deviation of the kth random variable.
In the case of the truncated normal distribution :
TAB1 . 'PARAM_VA' . k . 'MOYENNE'
TAB1 . 'PARAM_VA' . k . 'ECART_TYPE'
TAB1 . 'PARAM_VA' . k . 'MIN'
TAB1 . 'PARAM_VA' . k . 'MAX'
the first and the second parameter are the mean value and the
standard deviation of the kth random variable.
MIN et MAX are two reals between which, range the random variable.
In the case of the truncated shifted exponential :
TAB1 . 'PARAM_VA' . k . 'LAMBDA'
TAB1 . 'PARAM_VA' . k . 'MU'
TAB1 . 'PARAM_VA' . k . 'MIN'
TAB1 . 'PARAM_VA' . k . 'MAX'
MIN et MAX are two reals between which, range the random variable.
In the case of the Gumbel max distribution :
TAB1 . 'PARAM_VA' . k . 'LAMBDA'
TAB1 . 'PARAM_VA' . k . 'MU'
The density is :
lambda*exp(-lambda*(x-mu)-exp(-lambda*(x-mu)))
In the case of the Frechet distribution :
TAB1 . 'PARAM_VA' . k . 'U'
TAB1 . 'PARAM_VA' . k . 'K'
TAB1 . 'PARAM_VA' . k . 'B'
The density is :
((u - b)/(x - b))**k * exp(- ((u - b)/(x - b))**k) * k / (x - b)
TAB1 . 'MATCOV' : listreel containing the correlation matrix when the
random variables are dependant. It's the Nataf transformation which
is used.
For a matrix | a b c |
| b d e |
| c e f |
you have to give (prog a b d c e f).
The distributions permitted are :
UNIFORM, Standard normal, Normal, Lognormal, Shifted exponential.
TAB1 . 'POINTS_INITIAUX' : table ranging from 1 to initial points number
containing listreels containing the coordinates in the physical space
of the initial points.
TAB1 . 'CRITERE' : (optional) flottant(real). convergence criterion
on the reliability index(by default 1e-4)
TAB1 . 'RESU' : table ranging from 1 to initial points number
containing the results of the calculus.
TAB1 . 'RESU' . i . 'INDFIAB' : listreel containing the Hasofer-Lind index
at each iteration for the ith initial point.
TAB1 . 'RESU' . i . '<>' : where
<> is a chain of 4 character corresponding
with the name of a random variable in the
listmot TAB1 . 'NOMS_DES_VARIABLES'. Listreels containing the value
of the random variable at each iteration.
TAB1 . 'RESU' . i . 'FACTEURS_DE_SENSIBILITES' : listreel containing
the sensibilities of each variable at the last iteration.
TAB1 . 'RESU' . i . 'VECTEURS_DES_SENSIBILITES' : In the case where the
variables are dependant, listreel containing
the sensibilities of each variable at the last iteration.
You could do a restart. You have to give TAB1.
TAB1 . RESU . REPRISE . DERNIER : integer. This integer is the index at which
iteration will restart. You can modify it to your needs.
TAB1 . EXPLORATION : logical. If vrai, in the case of a restart
the program will search which value to give to vxmin and vxmax from the value
in TAB1.
TAB1 . RESU . i . INDICE_FIABILITE : Float. the fiability index must be
lower than it.
TAB1 . RESU . i . ATTEINT_EN : Integer.the index of the preceding result
in TAB1 . RESU . i . indfiab.