Opérateur POST_DYNA_ALEA
Contents
1. Dans cette exemple on effectue un r chantillonnage N nbtr 50 pour estimer les fractiles de la courbe 5 50 m diane et 95 et on d termine les enveloppes 100 et 0 5 2 Mot cl INTE SPEC Premier exemple POSTALEA POST_DYNA_ALEA INTE SPEC INTERS TOUT ORDRE OUI OPTION DIAG Deuxi me exemple POSTALEA POST DYNA ALEA INTE SPEC INTERS NOEUD I N1 NOM CMP _I DX NOEUD J N1 NOM CMP _J DX Manuel d utilisation Fascicule u4 84 Dynamique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html2. R sultats fournis Manuel d utilisation Fascicule u4 84 Dynamique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster so Titre Op rateur POS T_DYNA ALEA Date 31 05 2011 Page 6 8 Responsable Irmela ZENTNER Cl U4 84 04 R vision 6350 4 1 Mot cl FRAGILITE Les param tres de la table en sortie sont PARAMETRES TYPE DESCRIPTION xM Titre de la table Param tre estim par maximum de vraisemblance partir de l chantillon original Param tre B estim par maximum de vraisemblance partir de l chantillon original PARA NOCI Valeurs du param tre de nocivit pour o lesquelles on value les courbes PFA Valeurs de la courbe de fragilit param tres AM et BETA FRACTILES R Valeurs des courbes pour le fractile f Manuel d utilisation Fascicule u4 84 Dynamique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster ba Titre Op rateur POS T_DYNA_ALEA Date 31 05 2011 Page 7 8 Responsable Irmela ZENTNER Cl U4 84 04 R vision 6350 4 2 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 e les moments spectraux e des param tres statistiques utiliser s il s agit d un autospectre e cart type facteur d irr gularit nombre moyen d extrema par seconde nombre
3. consulter la documentation de la commande DYNA AL pour davantage d informations sur le sens des param tres du mot cl FA MODAL U4 53 22 NUME VITE FLUI nume nume est le num ro d ordre lorsque la table de linterspectre contient plusieurs matrices interspectrales cas des calculs de couplage fluide structure param tr es par une vitesse du fluide incident TOUT ORDRE OUl Toutes les tables d interspectre sont prises en compte On d finit ensuite les termes de la des matrice s dont les fonctions vont subir le traitement L ORDRE J lnumi lnumj NUME ORDRE NUME Manuel d utilisation Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Fascicule u4 84 Dynamique Code Aster oo Titre Op rateur POS T_DYNA_ ALEA Date 31 05 2011 Page 5 8 Responsable Irmela ZENTNER Cl U4 84 04 R vision 6350 3 2 1 3 3 3 4 4 Lorsque les autospectres ou les interspectres sont calcul s sur les modes e Inumi est la liste des num ros d ordre des modes i Exemple 2 3 1 e Inumj 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 e 2 2 correspond l autospectre sur le mode 2 e 3 1 correspond l interspectre entre le mode 3 et le mode 1 lnumi et 1numj doivent contenir le m me nombre de termes NOEUD I
4. listr8 3 1 3 Op rande AM INI et BETA INI ALPHA INI BETA INI On donne des valeurs initiales des param tres et B estimer point de d marrage pour l algorithme d optimisation 3 1 4 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 par la m thode de r chantillonnage m thode dite de bootstrap L op rande FRACTILES permet de donner les fractiles qu on souhaite calculer 0 FRACTILES fract listr8 Par d faut on tire autant d chantillons bootstrap qu on dispose de donn es c est le nombre N de simulation de Monte Carlo effectu au pr alable et dont les r sultats sont stock es dans la table Manuel d utilisation Fascicule u4 84 Dynamique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur POS T_DYNA ALEA Responsable lrmela ZENTNER Date 31 05 2011 Page 4 8 Cl U4 84 04 R vision 6350 TABL RESU La commande NB TIRAGE permet n anmoins de diminuer le nombre de tirage effectuer NB TIRAGE nbt 1 Il faut que nbt soit inf rieur ou gale au nombre de valeurs dans TAI BL RI ESU nbt lt N Cette fonctionnalit permet de r duire le temps de calcul mais est d conseill dans le cas g n ral car les r s
5. r alis s Cet op rateur produit une table de type table sdaster imprimable par IMPR TABLE U4 91 03 Manuel d utilisation Fascicule u4 84 Dynamique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur POST_DYNA ALEA Date 31 05 2011 Page 2 8 Responsable lrmela ZENTNER Cl U4 84 04 R vision 6350 2 Syntaxe table sdaster POST DYNA ALEA FRAGILITE F TABL RESU tabres VALE liste LIST PARA laster ALPHA INI amO0 0 4 O BETA INI beta0 f 053 FRACTILES fract NB TIRAGE nbt INTE SPEC inter INTE SPEC est renseign NUME VITE FLUI nume TOUT ORDRE OUI NUME ORDRE I Inumi NUME ORDRE J lnumj NOEUD I Inoeudi NOEUD J Inoeud NOM CMP I lcmpi NOM CMP J lcmpj OPTION DIAG MOMENT Imom INFO T La 2 TITRE titre table sdaster 1 R listr8 R DEFAUT R DEFAUT listr8 T1 table fonction DEFAUT l Kn Manuel d utilisation Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Fascicule u4 84 Dynamique Code Aster ba Ti
6. Code Aster AO Titre Op rateur POS T_DYNA _ALEA Date 31 05 2011 Page 1 8 Responsable Irmela ZENTNER Cl U4 84 04 R vision 6350 Op rateur POST DYNA ALEA 1 But POST DYNA ALEA permet d effectuer deux types de post traitement l issue d un calcul de dynamique stochastique Calcul de courbes de fragilit partir d une table contenant les r sultats d une simulation de Monte Carlo A partir d une table table sdaster contenant l information sur les niveaux d excitation en analyse sismique on choisit en g n ral le PGA ainsi que sur la d faillance ou non de la structure pour cette excitation POST DYNA ALEA permet de d terminer les param tres d une courbe de fragilit selon le mod le lognormale et de calculer des valeurs de cette courbe Le lecteur peut consulter U2 08 05 pour une description plus d taill e Post traitement statistiquement des r sultats de type interspectre POST DYNA ALEA permet sur des fonctions s lectionn es dans un concept de type table fonction de calculer des param tres statistiques moments spectraux cart type distribution des pics fr quence centrale Les tables d interspectres sont obtenues par diff rents op rateurs LIRE INTE SPEC U4 36 01 CALC INTE SPEC U4 36 03 DEFI INTE SPEC U4 36 021 DYNA ALEA MODAL U4 53 221 DYNA SPEC MODAL U4 53 23 ou REST SPEC PHYS U4 63 22 On se reportera R7 10 01 pour la description des traitements
7. de passages par z ro par seconde fr quence centrale Les param tres de cette table sont PARAMETRES TYP DESCRIPTION E 5 x I I N N CR EE ECART Si INFO 1 on imprime dans le fichier MESSAGE e le nom utilisateur de la table e les deux indices les 2 n uds ou les 2 modes de la fonction s lectionn e e le type de r sultat calcul e les options de calculs choisies ou prises par d faut e les valeurs des fonctions s lectionn es Manuel d utilisation Fascicule u4 84 Dynamique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster so Titre Op rateur POS T_DYNA ALEA Date 31 05 2011 Page 8 8 Responsable lrmela ZENTNER Cl U4 84 04 R vision 6350 5 Exemple 5 1 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 50000E 01 0 00000E 01 0 00000E 01 1 50000E 01 0 50000E 01 0 00000E 01 1 00000E 01 OO NN H U e Exemple du calcul d une courbe de fragilit TAB POST POST_DYNA_ALEA FRAGILITE F TABL RESU TAB1 LIST PARA Ir E AM INI 0 3 BETA INI 0 1 FRACTILE 0 0 0 05 0 5 0 95 1 0 NB TIRAGE 50 s TITRE courbe 1 INFO 2
8. lnoeudi NOEUD J Inoeudj 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 e Inoeudi est la liste des n uds suivant i NO92 NO95 NO98 e _lnoeudj est la liste des n uds suivant j NO92 NO92 NO92 e lcmpi est la liste des composantes suivant i DX DX DY e lcmpj est la liste des composantes suivant j DX DX DX Les n uds et composantes sont appair s suivant le m me rang e NO92 DX NO92 DX correspond l autospectre au n ud NO92 dans la direction DX e NO98 DY NO92 DX correspond l interspectre entre le n ud NO92 dans la direction DX et le noeud NO95 dans la direction DY Inoeudi Inoeudj lcmpi et lcpmj doivent contenir le m me nombre de termes OPTION DIAG Les calculs sont effectu s sur l ensemble des autospectres de la matrice et uniquement pour ceux l Op rande MOMENT MOMENT Imom Imom 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
9. tre Op rateur POST _ DYNA_ALEA Date 31 05 2011 Page 3 8 Responsable Irmela ZENTNER Cl U4 84 04 R vision 6350 3 Op randes 3 1 Mot cl FRAGILITE FRAGILITE Le mot_cl FRAGILITE permet de d terminer les param tres 4 et B m diane et cart type logarithmique d une courbe de fragilit selon le mod le log normale U2 08 05 _ In a 4 P ATF Ro On peut galement calculer les valeurs de la courbe pour les valeurs de param tres Am et B 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 3 1 1 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 NOCT c est l indicateur 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 3 1 2 Op randes LIST VAL 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 dis a 0 VALE liste 1 R ou en donnant le nom du concept de type listr8 contenant la liste des valeurs 0 LIST VAL liste
10. ultats sont peu fiables 3 2 Mot cl INTE SPEC INTE SP EC inter inter est le nom utilisateur de la table d interspectres La table d interspectres 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 Les param tres n cessaires de la table sont FONCTION K24 si la table contient au moins une matrice interspectrale NOM CHAM K16 VITE FLUIDE IR si les autospectres ou les interspectres sont calcul s sur les modes YNUME ORDRE I I NUME ORDRE J I si les autospectres ou les interspectres sont calcul s sur les n uds YNOEUD I K8 YNOEUD J K8 NOM CMP l K8 NOM CMP J K8 Remarque Les moments spectraux sont d finis comme des int grales de la densit spectrale de puissance DSP lwl S ylw d w 2 f w S ylw do 0 0 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 w gt 0 Par ailleurs les DSP sont d finis en fonction de la fr quence naturelle f 2rrw Hz dans POST DYNA ALEA On utilise les formules suivantes cf R7 010 01 re Rele 4e 1 S al S xle Le lecteur est invit
Download Pdf Manuals
Related Search