u4.84.04 Opérateur POST_DYNA_ALEA#

Syntaxe#

Opérandes#

Mot clé FRAGILITE#

FRAGILITE =

Le mot_clé FRAGILITEpermet de déterminer les paramètres \({A}_{m}\) et \(\beta\) (médiane et écart-type logarithmique) d’une courbe de fragilité selon le modèle log-normale[U2.08.05] :

\({P}_{f\mid a}=\Phi (\frac{\ln(a/{A}_{m})}{\beta})\)

On peut également calculer les valeurs de la courbe pour les valeurs de paramètres Am et \(\beta\) obtenues. L’option FRACTILES (facultatif) permet en outre de déterminer des fractiles pour la courbe par une méthode de rééchantillonnage de l’échantillon original qu’on a renseigné dans TABL_RESU.

Opérande TABL_RESU#

♦TAB_RESU = tabres            [table_sdaster]

On donne le nom de la table [table_sdaster] qu’on doit avoir crée auparavant à l’aide de CREA_TABLE [U4.33.02] . Cette table doit avoir au moins deux colonnes avec clés d’accès (nom de label de colonne): PARA_NOCI (c’est l’indicate fract = 0.5 ur caractérisant le niveau de l’excitation) et DEFA (les valeurs de cette colonne sont 0 si on n’a pas observé de défaillance ou 1 s’il y a eu défaillance) ou DEMANDE(les valeurs de la variable d’intérêt réelle caractérisant la défaillance ou l’endommagement, appelé aussi demande sismique dans ma littérature).

Opérande METHODE#

♦ METHODE=/ “EMV”

/ 'REGRESSION'

On choisit entre les deux méthodes pour le calcul de la courbe de fragilité lognormale: EMVpour l’estimation par maximum de vraisemblance ou REGRESSIONpour la régression linéaire. On trouve plus de détails sur ces deux méthodes dans la documentation [U2,08.05].Si l’on choisit REGRESSION, alors la tabletabresdoit contenir une colonne DEMANDErenseignant la demande sismique (variable d’intérêt comme le drift, une contrainte maximale, …) et il faut renseigner le seuil de défaillance via le mot-clé SEUIL.

Opérande SEUIL#

◊    SEUIL=    seuil                    [R]

Si la table TAB_RESU contient une colonneDEMANDE, alors il faut renseigner le seuil de cette variable à partir duquel la structure est considérée défaillante.

Opérandes LIST_PARA et VALE#

On peut donner une liste de réels, valeurs pour lesquelles on évalue la courbe de fragilité.

Ceci peut se faire sous forme d’une liste contenant les valeurs de calcul \(({a}_{1,}{a}_{2,}\mathrm{...},{a}_{n})\) :

◊VALE =liste [l_R]

ou en donnant le nom du concept de type listr8 contenant la liste des valeurs:

◊LIST_PARA =liste listr8

Opérande AM_INI et BETA_INI#

♦  AM_INI

◊ BETA_INI

Si on a choisi METHODE=”EMV”, alors il est impératif de donner une valeur initiale pour l’estimation du paramètre \({A}_{m}\) et il est conseillé de donner une estimation initiale pour \(\beta\) (points de démarrage pour l’algorithme d’optimisation).

Opérandes FRACTILES et NB_TIRAGE#

Ces opérandes doivent être renseignées si on souhaite déterminer des intervalles de confiances ou plus précisément des fractiles pour la courbe de fragilité estimé par la méthode du maximum de vraisemblance ('EMV'). La méthode de rééchantillonnage (dite aussi «bootstrap» dans la littérature anglo-saxonne) est utilisée pour cela. L’opérande FRACTILES permet de donner les fractiles qu‘on souhaite calculer.

◊FRACTILES = fract                    listr8

Par défaut, on tire autant d’échantillons «bootstrap» qu’on dispose de données (c’est le nombre :math:`N`de simulation de Monte Carlo effectué au préalable et dont les résultats sont stockées dans la table TABL_RESU). La commande NB_TIRAGE permet néanmoins de diminuer le nombre de tirage à effectuer:

◊NB_TIRAGE = nbt                [I]

Il faut que nbt soit inférieur ou égale au nombre de valeurs dans TABL_RESU (\(\mathrm{nbt}\le N\) ). Cette fonctionnalité permet de réduire le temps de calcul mais est déconseillé dans le cas général car les résultats sont peu fiables.

Mot clé INTERSPECTRE#

Opérande INTE_SPEC#

♦    INTE_SPEC = inter

inter est le nom utilisateur de la matrice interspectrale.

La matrice interspectrale peut être obtenue par différents opérateurs: lire_inte_spec [U4.36.01], calc_inte_spec [U4.36.03], defi_inte_spec [U4.36.02], dyna_alea_modal [U4.53.22], DYNA_SPEC_MODAL [U4.53.23] ou REST_SPEC_PHYS [U4.63.22].

Remarque :

Les moments spectraux sont définis comme des intégrales de la densité spectrale de puissance (DSP):

\({\lambda}_{i}=\underset{-\infty }{\overset{+\infty }{\int}}{\mid \omega \mid }^{i}{S}_{\text{XX}}(\omega )d\omega =2\underset{0}{\overset{+\infty }{\int}}{\omega}^{i}{S}_{\text{XX}}(\omega )d\omega\)

Ainsi, si la DSP est donnée pour les fréquences positives uniquement, POST_DYNA_ALEA multiplie par 2 les intégrales des DSP calculés pour \(\omega >0\) *. Par ailleurs, les DSP sont définis en fonction de la fréquence naturelle* \(f=2\pi \omega\) ( \(\mathit{Hz}\) ) dans POST_DYNA_ALEA . On utilise les formules suivantes [cf. R7.010.01]:

\(\begin{array}{}{S}_{\text{XX}}(f)=\underset{-\infty }{\overset{+\infty }{\int}}{R}_{\text{XX}}(\tau ){e}^{-\mathrm{2i}\pi f\tau }d\tau ;\\ S{'}_{\text{XX}}(\omega )=\frac{1}{2\pi }{S}_{\text{XX}}(f)\end{array}\)

Le lecteur est invité à consulter la documentation de la commande DYNA_ALEA_MODAL [U4.53.22] pour davantage d’informations sur le sens des paramètres du mot-clé.

Opérandes NUME_ORDRE_I, NUME_ORDRE_J#

♦ / ♦ NUME_ORDRE_I = lnumi ♦ NUME_ORDRE_J = lnumj

Ces mots-clés permettent de définir les termes de la matrice dont les fonctions vont subir le traitement.

Lorsque les autospectres ou les interspectres sont calculés sur les modes :

lnumi est la liste des numéros d’ordre des modes “i”. Exemple : (2,3,1).
lnumj est la liste des numéros d’odre des modes “j”. Exemple : (2,1,4)

Les indices sont appairés suivant le même rang.

(2,2) correspond à l’autospectre sur le mode 2,
(3,1) correspond à l’interspectre entre le mode 3 et le mode 1.

lnumi et lnumj doivent contenir le même nombre de termes.

Opérandes NOEUD_I, NOEUD_J, NOM_CMP_I, NOM_CMP_J#

♦    /    ♦    NOEUD_I = lnoeudi

♦ NOEUD_J = lnoeudj ♦ NOM_CMP_I = lcmpi ♦ NOM_CMP_J = lcmpj

Lorsque les autospectres ou les interspectres sont calculés sur les nœuds dans une direction donnée :

lnoeudi est la liste des nœuds suivant « i » : \((\mathit{NO92},\mathit{NO95},\mathit{NO98})\)
lnoeudj est la liste des nœuds suivant « j » : \((\mathit{NO92},\mathit{NO92},\mathit{NO92})\)
lcmpi est la liste des composantes suivant « i » : \((\mathit{DX},\mathit{DX},\mathit{DY})\)
lcmpj est la liste des composantes suivant « j » : \((\mathit{DX},\mathit{DX},\mathit{DX})\)

Les nœuds et composantes sont appairés suivant le même rang :

\((\mathit{NO92}\mathit{DX},\mathit{NO92}\mathit{DX})\) correspond à l’autospectre au nœud \(\mathit{NO92}\) dans la direction \(\mathit{DX}\) ,
\((\mathit{NO98}\mathit{DY},\mathit{NO92}\mathit{DX})\) correspond à l’interspectre entre le nœud \(\mathit{NO92}\) dans la direction \(\mathit{DX}\) et le noeud \(\mathit{NO95}\) dans la direction \(\mathit{DY}\) .

lnoeudi, lnoeudj, lcmpi et lcpmj doivent contenir le même nombre de termes.

Opérande OPTION#

♦ / OPTION = “TOUT”

Les calculs sont effectués sur l’ensemble des interspectres de la matrice.

/ OPTION = “DIAG”

Les calculs sont effectués sur l’ensemble des autospectres de la matrice et uniquement pour ceux-là.

fract = 0.5

Mot clé DUREE et FRACTILE#

◊    DUREE = duree

Si le mot-clé duree est renseigné, alors on détermine le maximum médian ou tout autre fractile ainsi que le facteur de pic du processus stochastique stationnaire Gaussien selon les formules de Vanmarcke. duree désigne alors l'intervalle du temps considéré pour estimer ces quantités. A titre d'exemple, dans le cadre d'une analyse sismique, duree peut être pris égal à la durée de la phase forte du signal sismique.

 ◊    FRACTILE    =     /frac            [R]

/0.5 [DEFAUT]

Permet de définir le fractile considéré pour les calculs du maximum.

Par défaut, on prend fract = 0.5, ce qui correspond à la médiane des maxima.

Opérande MOMENT#

◊    MOMENT = lmom

lmom est la liste des ordres des moments spectraux qui seront calculés. Par défaut, les moments spectraux d’ordres 0, 1, 2, 3 et 4 sont toujours calculés. Il convient donc de mentionner dans cette liste les moments d’ordre supérieur à 4. Exemple: (5,7,8).

Opérande INFO#

◊    INFO=

1	impression des résultats demandés.
2	comme 1 mais avec plus de détails.

Opérande TITRE#

◊    TITRE= titre

titre est le titre du calcul. Il sera imprimé en tête des résultats. Voir [U4.03.01].

Résultats fournis#

Mot-clé FRAGILITE#

Les paramètres de la table en sortie sont:

PARAMETRES	TYPE	DESCRIPTION
TITRE	TXM	Titre de la table
AM	R	Paramètre \({A}_{m}\) estimé par maximum de vraisemblance à partir de l’échantillon original
BETA	R	Paramètre \(\beta\) estimé par maximum de vraisemblance à partir de l’échantillon original
PARA_NOCI	R	Valeurs du paramètre de nocivité pour lesquelles on évalue les courbes
PFA	R	Valeurs de la courbe de fragilité (paramètres AMet BETA)
FRACTILES	R	Valeurs des courbes pour le fractile \(f\)

Mot-cle INTE_SPEC#

Pour chaque fonction choisie dans l’interspectre, POST_DYNA_ALEA stocke dans une table accessible par IMPR_TABLE [U4.91.03]

les moments spectraux
des paramètres statistiques (à utiliser s’il s’agit d’un autospectre):
écart-type,
facteur d’irrégularité,
nombre moyen d’extrema par seconde,
nombre de passages par zéro par seconde,
fréquence centrale
le facteur de pic selon la formule de Vanmarcke
le maximum moyen selon la formule de Vanmarcke

Les paramètres de cette table sont:

PARAMETRES	TYPE	DESCRIPTION
NUME_ORDRE_I	I	numéro d’ordre des modes \(i\)
NUME_ORDRE_J	I	numéro d’ordre des modes \(j\)
NOEUD_I	NO	Noeud \(i\)
NOEUD_J	NO	Noeud \(j\)
NOM_CMP_I	TXM	Nom de la composante au noeud \(i\) \((\mathit{DX},\mathit{DY},\mathit{DY})\)
NOM_CMP_J	TXM	Nom de la composante au noeud \(j\) \((\mathit{DX},\mathit{DY},\mathit{DY})\)
LAMBDA_00	R	moment spectral d’ordre 0
LAMBDA_01	R	moment spectral d’ordre 1
LAMBDA_02	R	moment spectral d’ordre 2
LAMBDA_03	R	moment spectral d’ordre 3
LAMBDA_04	R	moment spectral d’ordre 4
ECART	R	écart-type
NB_EXTREMA_P_S	R	nombre moyen d’extrema par seconde
NB_PASS_ZERO_P_S	R	nombre de passages par zéro par seconde
FREQ_APPAR	R	fréquence centrale
FACT_IRRE	R	facteur d’irrégularité
MAX_MOY	R	Maximum moyen
FACT_PIC	R	Facteur de pic

Si INFO = 1 on imprime dans le fichier MESSAGE

le nom utilisateur de la table,
les deux indices (les 2 nœuds ou les 2 modes) de la fonction sélectionnée,
le type de résultat calculé,
les options de calculs choisies ou prises par défaut,
les valeurs des fonctions sélectionnées.

Exemple#

Mot-clé FRAGILITE#

Exemple d’une table générées au préalable, en faisant appel à CALC_TABLE, lors de la simulation de Monte Carlo (voir aussi [U2.08.05]):

#TABLE_SDASTER

PARA_NOCI DEFA

5.00000E-01 1

4.50000E-01 0

3.00000E-01 0

3.00000E-01 1

1.50000E-01 0

2.50000E-01 0

9.00000E-01 1

4.00000E-01 1

\(\mathrm{⋮}\) \(⋮\)

Exemple du calcul d’une courbe de fragilité:

TAB_POST=POST_DYNA_ALEA( FRAGILITE=(_F(TABL_RESU=TAB1,

LIST_PARA=lr,

AM_INI =0.3 ,

BETA_INI=0.1 ,

FRACTILE = (0.0,0.05,0.5,0.95,1.0),

NB_TIRAGE =50,

),),

TITRE = “courbe 1”,

INFO=2,);

Dans cette exemple, on effectue un rééchantillonnage (\(N=\mathit{nbtr}=50\) ) pour estimer les fractiles de la courbe 5%, 50% (médiane) et 95% et on détermine les enveloppes (100% et 0%).

Mot-clé INTERSPECTRE#

Premier exemple :

POSTALEA =POST_DYNA_ALEA( INTERSPECTRE = _F( INTE_SPEC= INTERS,

OPTION=”DIAG” )

)

Deuxième exemple :

POSTALEA=POST_DYNA_ALEA(INTERSPECTRE = _F( INTE_SPEC=INTERS,

NOEUD_I=”N1”, NOM_CMP_I=”DX”, NOEUD_J=”N1”, NOM_CMP_J=”DX”, )

)