[U4.53.21] Opérateur DYNA_TRAN_MODAL
Contents
1. Date 10 07 2015 Page 6 36 pour les rotors Cl U4 53 21 GROUP NO 1 grnol 0 NOEUD 2 no2 GROUP NO 2 grno2 MAILLE ma GROUP MA grma OBSTACLE obs NORM OBST nor 0 ORIG OBST ori 0 JEU aL eu ANGL VRIL gamma 0 DIST 1 distl 0 DIST 2 dist2 0 SOUS STRUC 1 ssl O SOUS STRUC 2 ss2 REPERE GLOBAL nom sst RIGI NOR kn 0 AMOR NOR 10 y 0 RIGI TAN 0 kt 0 AMOR TAN ct FROTTEMENT NON COULOMB COULOMB mu COULOMB STAT DYNA COULOMB STAT mus COULOMB DYNA mud R R vision 13175 up_no ud up no D K8 K8 EFAUT K8 mmama R EFAUT EFAUT R R mamanman EFAUT R R Op randes sp cifiques la prise en compte d un transitoire de vitesse vitesse de rotation variable 0 VITESSE VARIABLE NON DEFAUT f OUT si VITESSE VARIABLE OUI VITE ROTA vrota fonction ATR GYRO gyro matr asse gene R 0 ACCE ROTA arota fonction B MATR RIGY gyro matr_asse_gene_R si VITESSE VARIABLE NON VITE ROTA 7 O07 DEFAUT vrota R Mot clef sp cifique a la prise en compte d une fissure dans un rotor ROTOR2. COUPLAGE EDYOS F VITE ROTA vrota PAS TPS EDYOS dtedyos PALIER EDYOS F 7 UNITE uled GROUP NO grnoed NOEUD noed 0 TYPE EDYOS PAPANL PAFINL PACONL PAHYNL 7 1 group_no noeud Fin des mots cl s facteurs associ uniquement au couplage avec le code EDYOS O ARCHIVAGE F E 0 CRITE 7 0 LIST _INST INST PAS ARCH RE SION O PRECI CRITERE PRECISION SOLVEUR F voir U4 50 01 0 INFO 1 J 2 0 IMPRESSION F 0 TOUT QUI NIVEAU INST_INIT ti INST FIN tf 0 UNITE DIS VISC 0 TITRE titre DEPL LOC VITE LOC FORC LOC list in ipa RELATIF ABSOLU prec TAUX CHOC listr8 R 1 DEFAUT DEFAUT R R DEFAUT DEFAUT R R I 1 Kn Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Eos Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 9 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 3 Op randes 3 1 Matrices g n ralis es Dans le cas d un calcul par recombinaison modale les matrices g n ralis es doivent tre tablies par
3. Titre Op rateur DYNA TRAN MODAL Date 10 07 2015 Page 28 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 3 13 3 Op randes li s au comportement 3 13 4 3 14 Le comportement DIS VISC est un comportement rh ologique visco lastique non lin aire de type Zener tendu permettant de sch matiser le comportement d un amortisseur uni axial entre deux n uds Pour la direction locale x et seulement celle l du dispositif on fournit cinq coefficients Leurs unit s doivent tre en accord avec l unit des efforts l unit des longueurs et l unit de temps du probl me K1 raideur lastique de l l ment 1 du mod le rh ologique K2 raideur lastique de l l ment 2 du mod le rh ologique e K3 raideur lastique de l l ment 3 du mod le rh ologique e UNSUR K1 souplesse lastique de l l ment 1 du mod le rh ologique e UNSUR K2 souplesse lastique de l l ment 2 du mod le rh ologique e UNSUR K3 souplesse lastique de l l ment 3 du mod le rh ologique e PUIS ALPHA puissance du comportement visqueux de l l ment q e C coefficient du comportement visqueux de l l ment Il existe des conditions respecter sur les valeurs des coefficients pour que la tangente soit toujours d finie ki gt 1078 uskI0 k3210 usk3 gt 0 usk2 gt 10 k2 gt 0 C gt 10 10 lt 0 lt 1 On ne peut donc pas avoir la fois uskIi 0 usk3 0 et k2 0 c est dire le ca
4. 3 12 Mot cl ANTI_SISM Le mot cl ANTI SISM est incompatible avec un calcul par sous structuration dynamique ll permet de calculer la force non lin aire qui existe si un dispositif anti sismique est plac entre les deux n uds antagonistes dont les noms sont pr cis s par les mots cl s NOEUD 1 ou GROUP NO 1 et NOEUD 2 OU GROUP NO 2 K K x i A a it xX E max 1 K 1 TAn P 0 RIGI K1 RIGI K2 SEUIL FX C PUIS ALPHA et DX MAX Param tres de la force due la pr sence d un dispositif anti sismique A titre d exemple les valeurs des param tres pour un dispositif anti sismique de type JARRET sont K1 6 E 06N m K2 053E 06N m Py 1200 C 0 07E 05 Nmls alpha 0 2 et xmax 0 03m si le probl me est pos en USI 3 13 Mot cl DIS_VISC C est un comportement visco lastique non lin aire entre deux n uds cf R5 03 17 Ce comportement n affecte que le degr de libert DX local de l l ment La direction x locale de l l ment va du n ud 1 au n ud 2 Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster ul Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 27 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 Remarque Les r sultats concernant l effort les d placements visqueux et relatifs entre les deux n uds ainsi que la dissipation du dispositif non lin ai
5. Pour les syst mes tr s faiblement amortis cette option permet d viter un calcul co teux de la phase lin aire pr c dant le premier choc Cette phase appel e par la suite phase transitoire pr c de l tablissement d un r gime constitu d une succession de phases non lin aires de chocs et ou de phases lin aires dites de vol selon les fonctions d excitation du syst me m canique appliqu es Le temps de transitoire correspond un d placement gal au jeu d une but e Il peut tre relativement important 50 100 secondes ETAT STAT OUI permet le passage en un seul pas de temps de calcul de la phase transitoire Le passage de la phase transitoire est r alis en supposant le syst me m canique en vol Le temps n cessaire au passage du transitoire est estim par l algorithme en fonction des caract ristiques m caniques du syst me en vol Cette estimation est bas e sur un crit re o interviennent le param tre PREC_ DUREE et les dur es d excitations dues aux efforts turbulents Remarque Si Pon demande une simulation avec calcul en un pas de temps de la phase transitoire il Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Version E default Titre Op rateur DYNA_TRAN_MODAL Date 10 07 2015 Page 14 36 Responsable Albe
6. Ce sont les op randes suivants du mot cl facteur INCREMENT NB POIN PERIODE N Nombre de points par p riode apparente C est ce param tre qui fixe la pr cision du calcul Il doit tre au moins gal 20 sa valeur par d faut 50 garantit une pr cision satisfaisante de l ordre de 1 dans la plupart des cas VITE MIN M thode de calcul de la vitesse de r f rence utilis e pour valuer la fr quence apparente Quand le d nominateur de la fr quence apparente X x _ devient faible celle ci peut devenir tr s lev e ce qui conduit un raffinement injustifi du pas de temps Pour y rem dier l algorithme utilise le crit re suivant Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster peras Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 11 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 rl l Xa Kni lt V e a At min far 271 VA V in peut tre calcul de deux fa ons diff rentes selon la valeur de VITE MIN Vit NORM V t EE pour tous les degr s de libert Peut tre utilis e si le syst me poss de plusieurs degr s de libert e si l ordre de grandeur du d placement n est pas trop diff rent selon les degr s de libert Max 1V t T MAXI
7. Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster ou Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 16 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 Si l utilisateur ne donne pas de valeur ce param tre facultatif les sch mas ADAPT ORDREl ADAPT ORDRE2 estimeront une valeur not e dts partir de la fr quence de coupure de la base ventuellement corrig e par les raideurs de chocs En revanche les sch mas de Runge Kutta n auront aucune limitation sup rieur en termes de pas de temps Pour retrouver le fonctionnement des versions ant rieures du code il suffit d imposer dtmax dt donc la m me valeur au param tre PAS qu PAS MAXI Si l utilisateur donne une valeur sup rieure dts une alarme sera mise pr venant d un risque de perte de pr cision PAS MINI dtmin Valeur minimale du pas de temps Si les conditions de diminution du pas de temps sont remplies le pas de temps courant pourra alors diminuer jusqu cette valeur limite Si l utilisateur ne donne pas de valeur ce param tre facultatif alors le code calculera le pas de temps minimal partir de PAS LIMI RELA Pour retrouver le fonctionnement des versions ant rieures du code il suffit donc de ne pas d finir PAS MINI 3 4 2 Op randes INST_INIT INST FIN NUME FIN 0 INST_INIT to e M thodes EULER D
8. NEWMARK sticitasiii A A At 34 4 3 Phase d ex cution pour les m thodes ADAPT_ORDRE et ADAPT_ORDRE2 34 4 4 Phase d ex cution pour la m thode TM 35 Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Fo Titre Op rateur DYNA_TRAN MODAL Date 10 07 2015 Page 4 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 2 Syntaxe tranmo tran gene DYNA TRAN MODAL 0 reuse tranmo MATR MASS ma matr asse gene R MATR RIGI ri matr asse gene R 0 MATR AMOR am matr asse gene R AMOR MODAL FE AMOR REDUIT la 1 R LIST AMOR l amor listr8 E 0 SCHEMA TEMPS F SCHEMA NEWMARK DEFAUT EULER 7 DEVOGE LV 21 Er El ADAPT ORDR RUNGE_KUTTA 54 RUNGE KUTTA 32 ADAPT_ORDR ITMI Mots cl s associ s uniquement au sch ma NEWMARK O BETA 0 25 DEFAUT beta R 0 GAMMA 0 5 DEFAUT gamma R Mots cl s associ s uniquement aux sch mas RUNGE KUTTA 0 TOLERANCE 1 E 3 DEFAUT tol R O ALPHA 1 E 3 DEFAUT alpha R Mots cl s associ s uniquement au sch ma I
9. intervalle de recherche 1 prec instant 1 prec instant ABSOLU intervalle de recherche instant prec instant prec Le crit re est RELATIF par d faut Op rande PRECISION PRECISION Pa prec E 06 DEFAUT R8 Indique avec quelle pr cision la recherche de l instant doit se faire Description du chargement sous variables s par es mot cl EXCIT 0 EXCIT Mot cl d finissant le chargement Ce mot cl doit tre r p t autant de fois qu il y a de vecteurs chargement g n ralis f Le chargement total est la somme de ces vecteurs chargement Ce mot cl permet de d finir le chargement sous a forme de vecteurs g n ralis s multipli s par des fonctions multiplicatrices Op randes VECT_ASSE GENE NUME ORDRE Le chargement est pris en compte sous forme de vecteur projet sur la base modale EXCIT F VECT ASSE GENE ou sous forme de composante modale EXCIT F NUME MODE ou les deux la fois VECT_ASSE GENE v Vecteur g n ralis permettant de d crire la r partition spatiale du chargement Concept de type vect_asse_ gene Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster cu Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 18 36 Responsable Albert ALARCON Cl
10. pas adaptatif Ce sch ma fait partie de la famille des sch mas d int gration de type Runge Kutta En particulier il s agit du sch ma d int gration explicite de Dormand Prince 54 R5 06 04 pas de temps adaptatif Le sch ma RUNGE_KUTTA_54 supporte la prise en compte de toutes les non lin arit s disponibles dans l op rateur Le calcul du pas de temps optimal se fait par contr le de l erreur entre les approximations d ordre 5 et 4 de la pr diction du vecteur d tat concat nation des vecteurs de d placement et vitesse Ce sch ma s appuie sur la condition de contr le de l erreur relative suivante err lt tol avec erro y 22n et sc MAX y y 11 Q ni SC 1 o e y est la valeur de la pr diction d ordre 5 de la composante du vecteur d tat y Yu est la valeur de la pr diction d ordre 4 de la composante du vecteur d tat y e n est la taille du vecteur d tat y e Jo est la valeur de la composante du vecteur d tat y l tat actuel 0 TOLERANCE tol Valeur de contr le d erreur relative donn e par l utilisateur Par d faut elle vaut 1 E 3 0 ALPHA alpha Valeur de regularisation donn e par l utilisateur intervenant dans l expression de sc Par d faut elle vaut 1 E 3 3 3 1 6 SCHEMA RUNGE KUTTA 32 sch ma explicite pas adaptatif Comme le sch ma RUNGE_KUTTA_54 le sch ma RUNGE KUTTA 32 fait partie de la famille des sch mas d int grat
11. thodes EULER DEVOGE NEWMARK ITMI Entier d finissant la p riodicit d archivage de la solution du calcul transitoire dans le concept r sultat tran gene Siipa 5 on archive tous les 5 pas de calcul Quelle que soit l option d archivage choisie on archive le dernier pas de temps et tous les champs associ s pour permettre une ventuelle reprise Par d faut on archive tous les pas de calcul e M thodes ADAPT ORDREl et ADAPT ORDRE Entier qui permet de calculer l intervalle entre deux instants d archivage dans le concept r sultat gal PAS ARCH PAS Avec cette convention le pas d archivage est toujours sup rieur ou gal au pas maximal utilis par le calcul Avec un pas variable les instants d archivage ne correspondent pas exactement des pas de calcul L algorithme archive donc les grandeurs aux pas de calcul les plus proches des instants d archivage indiqu s par l utilisateur en Tn sur ce sch ma Pas de calcul Pas d archivage Instants d archivage 3 19 3 Op rande CRITERE 0 CRITERE Indique avec quelle pr cision la recherche de l instant archiver doit se faire RELATIF intervalle de recherche 1 prec instant 1 prec instant ABSOLU intervalle de recherche instant prec instant prec La valeur par d faut du crit re de recherche est RELATIF 3 19 4 Op rande PRECISION 0 PRECISION 1 E 06 DEFAUT prec R Indiq
12. BASE U4 63 13 partir de la base modale utilis e pour les matrices de rigidit g n ralis es ou par l op rateur RECU GENE U4 71 03 appliqu un calcul pr c dent Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Version default Titre Op rateur DYNA TRAN_MODAL Responsable Albert ALARCON 3 5 1 3 5 2 3 5 3 3 5 4 3 6 3 6 1 Date 10 07 2015 Page 17 36 Cl U4 53 21 R vision 13175 Op randes RESU DEPL VITE RESULTAT tran Concept de type tran gene issu d un calcul pr c dent avec DYNA TRAN MODAL l DEPL do Concept de type vect_asse_gene d placements g n ralis s initiaux l VITE vo Concept de type vect_asse_gene vitesses g n ralis es initiales Op randes INST_INIT NUME ORDR INST_INIT to Instant du calcul pr c dent extraire et prendre comme tat initial dans le cas d une reprise En l absence de cet op rande l instant de reprise est pris gal au dernier instant de calcul pr c dent archiv NUME_ORDRE nuord D signe le num ro d archivage du calcul pr c dent extraire et prendre comme tat initial dans le cas d une reprise Op rande CRITERE Q CRITERE Indique avec quelle pr cision la recherche de l instant doit se faire RELATIF
13. FISS F NOEUD G nog noeud GROUP NO G grnog group_no NOEUD D nod noeud GROUP NO D grnod group_no ANGL _INIT O Ol DEFAUT O ANGL ROTA 0 0 fonction K PHI kphi fonction DK DPHI dkdphi fonction O VERI CHOC F Q STOP CRITERE OUL DEFAUT E NON Q SEUIL 0 51 DEFAUT S R Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Version default Titre Op rateur DYNA_TRAN_MODAL Responsable Albert ALARCON 0 0 ANTI_SISM E 7 NOEUD 1 nol GROUP NO 1 grnol NOEUD2 no2 GROUP NO 2 grno2 0 RIGI K1 Ebo E kn 0 RIGI K2 Des kn SEUIL FX 0 Py 0 C Os Cs 0 PUIS ALPHA 0 alpha 0 DX MAX Re dx DIS VISC F NOEUD 1 nol GROUP NO 1 grnol NOEUD 2 no2 GROUP NO 2 grno2 K1 kl UNSUR K1 uskl K2 k2 UNSUR K2 usk2 K3 k3 UNSUR_K3 usk3 E C PUIS ALPHA 05 alpha 0 ITER INTE MAXI 20 iter 0 RESI_INTE_RELA 1 0E 06 resi FLAMBAGE E NOEUD 1 nol GROUP NO 1 grnol 0 NOEUD 2 no GROUP NO 2 grno2 4 OBSTACLE obs Q ORIG OBST ori NORM OBST nor A
14. JEU O JEU jeu Dans le cas d un choc entre une structure mobile et un obstacle ind formable l op rande 91 repr sente e la demi distance inter plans pour des obstacles de type PLAN Y et PLAN Z e le rayon de l obstacle circulaire pour un obstacle de type CERCLE Ce mot cl est inutilis dans le cas d obstacles discr tis s par segments de type DISCRET El G Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster ul Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 22 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 Remarque L obstacle de type PLAN Y ou PLAN Z comporte en fait deux obstacles plans Ainsi dans le cas o l utilisateur souhaite mod liser le choc sur un plan unique pour ne pas tre g n par le rebond de la structure tudi e sur le plan sym trique on conseille l utilisateur de le repousser tr s loin cf Figure 3 6 1 6 a j repr sente le jeu r el entre la structure tudi e et l obstacle Zloc orig_obs Figure 3 6 1 6 a Syst me masse ressort impactant une paroi fixe Remarque Le mot cl JEU n est pas utilis dans le cas de choc entre structures mobiles Les diff rents cas de jeux sont repr sent s dans la documentation de DEFI OBSTACLE U4 44 21 3 10 1 7 Op rande ANGL_VRIL 0 ANGL VRIL g
15. compte de non lin arit s localis es aux n uds de type chocs et frottement ainsi que des non lin arit s de type dispositif anti sismique Pour pouvoir calculer la r ponse d une machine tournante sur paliers non lin aires un coupleur avec le code EDYOS a t d velopp et n cessite l utilisation conjointe des mot cl s facteurs COUPLAGE EDYOS et PALIER EDYOS La structure de donn es resultat contient pour diff rents instants de calcul les r sultats g n ralis s et les forces de choc calcul es La conversion des r sultats g n ralis s dans l espace physique est possible par les op rateurs REST GENE PHYS U4 63 31 ou pour une composante par RECU FONCTION U4 32 03 Produit un concept de type tran gene Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Version default Titre Op rateur DYNA_TRAN_MODAL Date 10 07 2015 Page 2 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 Table des mati res Bl bie pea 1 A E O El tester a rss te tte 4 3 Op randes tt ina mel Rene A AS 9 31 Matrices generalista nee mit ans tele pe nl A nr Pare eaan 9 3 2 Mot cl AMOR MODAL scenie e EAEE Ea EE E A EEA EN 9 3 2 1 Op randes AMOR_REDUIT LIST_AMOR coccconcccnccnncn
16. e Ence qui concerne les autres types d obstacles la valeur par d faut de gamma est z ro 3 10 1 8 Op randes DIST_1 DIST_2 DIST_1 disti Distance caract ristique de mati re entourant NOEUD 1 nol ou GROUP NO 1 Op rande sp cifique au contact entre deux structures mobiles DIST 2 dist2 Distance caract ristique de mati re entourant NOEUD 2 no2 ou GROUP NO 2 Op rande sp cifique au contact entre deux structures mobiles Remarques Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster dl Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 24 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 e DIST 1 et DIST 2 sont d finies au sens des normales sortantes des deux solides en vis vis DIST_1 et DIST 2 sont positives car elles repr sentent l paisseur des structures tudi es e Du fait du calcul de la distance normale de choc la somme de DIST 1 et de DIST 2 doit tre suffisamment grande par rapport l amplitude suppos e du d placement relatif des n uds de chocs cf R5 06 03 3 10 1 9 Op randes SOUS_STRUC_1 SOUS_STRUC_2 0 SOUS STRUC_1 ssl Nom de la sous structure qui contient le n ud de choc renseignant le mot cl NOEUD 1 ou GROUP NO 1 SOUS STRUC 2 ss2 Nom de la sous structure qui contient le n ud de choc
17. pend pas de la g om trie du rotor mais uniquement de la forme de la fissure et d un coefficient de dimension 0 ROTOR FISS F NOEUD G l vre gauche de la fissure GROUP NO G groupe de n ud nommant la l vre gauche de la fissure NOEUD D l vre droite de la fissure GROUP NO D groupe de n ud nommant la l vre gauche de la fissure ANGL INIT angle initial du fond de fissure par rapport sa d finition dans la loi de comportement de fissure par d faut 0 0 0 ANGL ROTA fonction donnant la loi impos e de position angulaire du fond de fissure par rapport sa d finition dans la loi de comportement de fissure transitoires de vitesse K PHI loi de comportement en raideur de la fissure DK DPHI d riv e de la loi de comportement en raideur Orientation de l axe du rotor Pour respecter le sens de rotation trigonom trique il est important de bien orienter le rotor l axe du rotor est automatiquement orient par la fissure en allant du bord gauche vers le bord droite de la fissure 3 11 Mot cl VERI CHOC Mot cl qui permet d valuer a posteriori l aptitude de la base modale repr senter correctement les impacts Si VERI CHOC est pr sent on calcule en chaque n ud de choc et pour chaque mode le taux de n Cd y reconstitution de la solution statique t K 77 et pour information le taux de s statique k i 1 i Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyri
18. revanche la m thode de calcul transitoire sur les bases des sous structures ne supporte que le sch ma d Euler et les sch mas pas de temps adaptatif 3 3 1 1 SCHEMA NEWMARK sch ma implicite Ce sch ma ne permet que l int gration de probl mes lin aires C est le sch ma par d faut pour la r solution On peut pr ciser les param tres d int gration B et y 0 BETA beta Valeur du param tre f pour la m thode de NEWMARK Par d faut B 0 0 25 0 GAMMA gamm Valeur du param tre y pour la m thode de NEWMARK Par d faut y 0 5 3 3 1 2 SCHEMA EULER sch ma explicite d ordre 1 Ce sch ma supporte le calcul avec prise en compte de l ensemble des non lin arit s localis es disponibles 3 3 1 3 SCHEMA DEVOGE sch ma explicite d ordre 4 Le sch ma de DEVOGELAERE supporte le calcul avec prise en compte de l ensemble des non lin arit s localis es disponibles 3 3 1 4 SCHEMA ADAPT_ORDRE2 sch ma explicite d ordre 2 Ce sch ma appel ADAPT dans les versions ant rieures du code supporte le calcul avec prise en compte de l ensemble des non lin arit s localis es disponibles Cette m thode utilise le sch ma des diff rences centr es l algorithme d adaptation du pas de temps s appuie sur le calcul d une fr quence apparente XT Xii aL Fa MT XX On pr cise ci apr s les op randes sp cifiques la m thode d int gration par pas de temps adaptatifs
19. 1 soit 99 53 de souplesse locale TAUX DE RESTIT EFFORT TRANCHANT 1 8979E 02 soit 1 89 de l effort tranchant On imprime ces quantit globalement pour l ensemble des modes et pour chaque mode Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster dl Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 34 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 On imprime en outre e pour chaque n ud de choc les rapports souplesse locale sur souplesse de choc et souplesse statique moins souplesse locale sur souplesse de choc e pour chaque mode sa participation sur les d form es statique aux n uds de choc Elle vaut le rapport du nombre de conditionnement de la matrice ferm e par le vecteur modal et les d form es statiques sur le nombre de conditionnement de la matrice des d form es statiques 3 21 Op rande IMPRESSION 0 IMPRESSION Mot cl facteur qui permet d imprimer dans le fichier RESULTAT des grandeurs non imprimables par un op rateur d impression telles que le d placement local la vitesse locale les forces de contact aux n uds de choc et la valeur cumul e sur tous les modes de la base modale de projection du taux de reconstitution de la solution statique 3 21 1 Op randes TOUT NIVEAU Le mot cl NIVEAU perm
20. Code Aster un Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 1 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 Op rateur DYNA TRAN MODAL 1 But Calculer la r ponse dynamique transitoire d un syst me amorti ou non en coordonn es g n ralis es Le calcul est effectu par superposition modale ou par sous structuration Des conditions initiales non nulles peuvent tre introduites permettant entre autre d utiliser les r sultats d un calcul ant rieur Le chargement est donn sous la forme d une combinaison lin aire de vecteurs g n ralis s et de fonctions du temps d crivant l volution temporelle de ces vecteurs Six m thodes d int gration explicites EULER DEVOGE ADAPT ORDRE1 et ADAPT ORDRE2 m thodes d int gration pas de temps adaptatif d ordre 1 ou d ordre2 RUNGE KUTTA 32 et RUNGE_KUTTA_54 m thodes d int gration pas de temps adaptatif de la famille de Runge Kutta ordres 54 et 32 respectivement une m thode int grale ITMI et une m thode d int gration implicite NEWMARK sont disponibles Les algorithmes explicites et ITMI supportent le calcul avec prise en compte de non lin arit s localis es aux n uds de type chocs et frottement Les m thodes EULER ADAPT ORDREl et ADAPT ORDRE2 supportent la prise en compte de non lin arit s de type dispositif anti sismique Les m thodes RUNGE KUTTA 32 et RUNGE KUTTA 54 supportent la prise en
21. ESSE VARIABLE OU I pour les transitoires de vitesse ou constante VITESSE VARIABLE NON 3 9 2 Op randes VITE ROTA MATR GYRO ACCE ROTA et MATR RIGY Ces op randes sont les param tres d finissant le transitoire de vitesse de rotation Si VITESSE VARIABLE OUI alors il faut renseigner les param tres suivants VITE ROTA fonction donnant la loi de vitesse angulaire impos e au rotor MATR GYRO matrice d amortissement gyroscopique 0 ACCE ROTA fonction donnant la loi d acc l ration angulaire impos e au rotor MATR RIGY matrice de raideur gyroscopique Remarque Si ACCE _ROTA n est pas fournie elle est n glig e ie pas de d rivation num rique partir de VITE ROTA Si VITESSE VARIABLE NON il faut renseigner la valeur de la vitesse de rotation constante VITE ROTA vit esse de rotation de l arbre par d faut 0 0 3 10 Prise en compte de non lin arit s localis es de type choc frottement et rotor fissur 3 10 1 Non lin arit s localis es de type choc et frottement mot cl CHOC 0 CHOC Ce mot cl facteur est utilis pour l tude de la r ponse de structures g n ralement lanc es dont les d placements sont limit s en un ou plusieurs point s pr cis s a priori par Putilisateur par la pr sence d un obstacle les diff rents types d obstacles disponibles sont d crit dans la documen
22. EVOGE NEWMARK YADAPT ORDRE1 et YADAPT ORDRE2 Instant de d but du calcul transitoire En cas de reprise on utilise le mot cl ETAT_INIT sous ce mot cl l instant initial est r cup r avec l op rande INST INIT ou pris gal au dernier instant de calcul pr c dent archiv L op rande INST_INIT doit donc tre utilis uniquement s il n y a pas de reprise d un calcul pr c dent e Sch ma ITMI D signe linstant de d but de simulation Lorsque le calcul en un pas de temps de la phase transitoire est demand le simulation d bute INS_INIT temps de calcul du transitoire INST_FIN tf Instant de fin de la simulation NUME_FIN tf Num ro de l instant de fin de calcul dans LIST_INST 3 5 Mot cl ETAT_INIT Mot cl facteur qui permet une poursuite d un calcul transitoire en prenant comme tat initial e soit un r sultat issu d un calcul par synth se modale pr c dent EXCIT RESULTAT e soit des d placements et vitesses exprim s sous forme de vecteurs assembl s g n ralis s EXCIT DEPL et VITE Remarques e Cette fonctionnalit n est pas disponible pour un calcul par sous structuration transitoire sans double projection ni pour la m thode ITMI e Lors d une poursuite l tat d adh rence ou de choc n est pas sauvegard e Les d placements et vitesses g n ralis s doivent tre tablis par l op rateur PROJ VECT
23. I 1 E PAFINL PACONL PAHYNL Nom du mod le de palier associ selon la terminologie du code EDYOS Seuls les types suivants de paliers non lin aires sont autoris s PAPANL paliers patins mobiles PAFINL palier patins fixes PACONL palier avec prise en compte de possibilit s de contact PAHYNL palier hybrides hydrodynamiques hydrostatiques avec poches de soul vement 3 19 Mot cl ARCHIVAGE 0 ARCHIVAGE Mot cl facteur d finissant l archivage Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster ul Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 32 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 Remarque Pour les sch mas d int gration de type Runge Kutta l archivage se fait de mani re syst matique sur tous les instants de calcul Ainsi les valeurs donn es dans les op randes LIST ARCH INST et PAS ARCH ne sont pas prises en compte 3 19 1 Op rande LIST _ARCH INST 0 LIST INST 1 arch Liste d entiers d finissant les instants de calcul pour lesquels la solution doit tre archiv e dans le concept r sultat tran_ gene 0 INST Instants de calcul pour lesquels la solution doit tre archiv e dans le concept r sultat tran gene 3 19 2 Op rande PAS ARCH 0 PAS ARCH ipa M
24. LIMI_RELA H 1 E 6 DEFAUT per R o NB_POIN PERIODE 50 DEFAUT ON I NMAX ITER PAS 7 16 DEFAUT N I Op randes sp cifiques une int gration par pas de temps adaptatifs ordre 1 et 2 ainsi que Runge Kutta PAS MAXI dtmax R 0 PAS MINI dtmin R EF RESULTAT res tran gene Si RESULTAT INST_INIT to R NUME_ORDR no I O CRITERE RELATIF DEFAUT Y PRECISION 1 E 06 DEFAUT prec R CRITERE ABSOLU PRECISION prec R DEPL do vect asse gene VITE vo vect asse gene F _F VECT_ASSE GENE v vect_asse gene NUME ORDRE nmordr I 0 FONC MULT fonction COEF MULT a R ACCE ac fonction O VITE vi fonction 0 DEPL dp fonction 0 MULT APPUI NON DEFAUT f OUL DIRECTION dx dy dz drx dry drz 1 R NOEUD lno 1 noeud GROUP NO lgrno 1 groupe no 0 CORR STAT OUI D FONC DT dfdt fonction D FONC DT2 dfdt2 fonction STAT psi mode meca CORR modcor mult elas mode meca RESULTAT resuforc tran gene COEF_MULT ai R COEF MULT C aci C 0 INTITULE int 1 Kn NOEUD 1 nol noeud Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur DYNA_TRAN_MODAL Responsable Albert ALARCON
25. NGL VRIL 0 gamma 0 JEU 1 jeu DIST 1 distl 0 DIST 2 dist2 0 REPERE GLOBAL nom sst o RIGI NOR kn 0 FNOR CRIT flim Q FNOR POST FL fseuil Q RIGI NOR POST FL k2 RELA EFFO DEPL F e B NOEUD noe Date 10 07 2015 Page 7 36 Cl U4 53 21 R vision 13175 roup no EFAUT EFAUT EFAUT EFAUT EFAUT EFAUT noeud group_no noeud group_no R R R R R u R R t R defa defaut E defaut R group_no noeud group_no obstacle listr8 listr8 DEFAUT EFAUT K8 R R R R R DEFAUT R noeud Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur DYNA_TRA N_MODAL Responsable Albert ALARCON 0 0 SOUS_STRUC NOM CMP RELATION ELA EFFO VITI e oo Ba NOEUD SOUS_STRUC NOM_CMP RELATION y ss nomcmp noe ss nomcmp Date 10 07 2015 Page 8 36 Cl U4 53 21 R vision 13175 K8 K8 fonction noeud K8 K8 fonction Mots cl s facteurs associ s uniquement au couplage avec le code EDYOS
26. N_MODAL Date 10 07 2015 Page 30 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 La relation non lin aire doit tre d finie sur sur 00 0o La phase non lin aire dans les post traitements correspond la plage d instants quand la relation non lin aire tait non nulle L quation d quilibre pour une structure soumise une acc l ration de sol horizontale a dans la direction x et ayant des termes de correction provenant de non lin arit s s crit Mxi Cx Kx Ma F o F est la force corrective due la non lin arit du sol Elle peut tre par exemple d finie par la relation suivante cf cas test SDND103 F kx f x avec Si f x k Lxl Xo Dans l exemple ci dessus on impose donc sous l op rande RELATION la fonction Y a 2 Il k Aloe x0 pour x gt x F x 0 pour x lt x 3 16 Mot cl RELA EFFO VITE 0 RELA EFFO VITE Mot cl facteur permettant de d finir une relation force vitesse sur un degr de libert d un n ud donn sous la forme d une fonction non lin aire Les op randes NOEUD SOUS STRUC NOM CMP et RELATION ont le m me sens pour les mots cl s RELA EFFO DEPL et RELA EFFO VITE lls ne sont donc pas d taill s dans ce paragraphe 3 17 Mot cl COUPLAGE EDYOS 0 COUPLAGE EDYOS Mot cl facteur obligatoire lors du couplage avec le code EDYOS pour calculer la r ponse de
27. PAS pour permettre de franchir la difficult avec un pas de temps plus petit COEF MULT PAS cmp Coefficient d augmentation du pas lorsque l erreur est suffisamment faible 0 75 Nf apn Sa valeur par d faut cmp 1 1 garantit stabilit et pr cision mais il peut en g n ral tre augment au plus jusqu 1 3 pour acc l rer l int gration A fa lt gt At cmp t COEF DIVI PAS cdp Coefficient de raffinement du pas de temps gt 1 lorsque l erreur est sup rieure 1 que le nombre d it rations maximales N_MAX ITER PAS n est pas atteint et que le pas de temps minimal n est pas atteint Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl htmil Code Aster ul Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 12 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 l A tas N ipn N iter Ne max et A ta plr A E enitial At cdp La valeur par d faut est de 1 33333334 soit une r duction d un facteur 0 75 PAS LIMI RELA plr Coefficient appliqu au pas de temps initial pour d finir la limite de raffinement et donc le pas de temps minimal La valeur par d faut est de 1 33333334 soit une r duction d un facteur 0 75 A T min plr A f initial 3 3 1 5 SCHEMA RUNGE_KUTTA_54 sch ma explicite
28. STAT 0UI 3 8 2 Prise en compte du multi appuis mots cl s MODE STAT MULTI APPUI et ACCE VITE DEPL Dans le cas d une structure multi support e afin de restituer les grandeurs calcul es dans le rep re absolu ou prendre en compte des non lin arit s localis es il faut calculer la r ponse g n ralis e en prenant en compte la composante d entra nement Pour plus de d tails on se reportera la r f rence R4 05 01 Les mots cl s MODE STAT MULT APPUI ACCE VITE DEPL DIRECTION et NOEUD ou GROUP NO sp cifiques la prise en compte du caract re multi support doivent tre simultan ment pr sents Un cham no issu de la projection d un CALC CHAR SEISME repr sente le vecteur d excitation sur l appui Il ne doit pas tre oubli m me si l information peut para tre redondante avec la donn e de l appui et de la direction de s isme 0 MODE STAT psi Concept de type mode meca produit par la commande MODE STATIQUE U4 52 14 qui correspond aux 3 ou 6 nb supports modes statiques o nb supports est le nombre de supports qui subissent une acc l ration diff rente EXCIT F MULT APPUI Si on calcule la r ponse sismique d une structure multi support e MULT APPUI OUI on compare chaque instant le vecteur des d placements absolus de chacun des points de choc consid r s afin de d terminer si il y a choc et de calculer les f
29. SULTAT resu gene Structure de donn es resu gene d finissant le chargement g n ralis 3 7 2 Op rande COEF MULT COEF MULT f Coefficient multiplicateur vaut 1 0 par d faut 3 8 Cas particulier de l analyse sismique 3 8 1 Prise en compte des modes n glig s par correction statique mots cl s CORR_STAT MODE CORR et D _FONC Lors de l analyse sismique d une structure mono excit e il est possible de prendre en compte a posteriori l effet statique des modes n glig s Dans ce cas lors du retour sur la base physique les d placements relatifs calcul s respectivement les vitesses et acc l rations relatives sont corrig s par un pseudo mode On trouvera les d tails de ce type de correction dans R4 05 01 l int rieur du mot cl facteur EXCIT CORR STAT OUI permet la prise en compte des modes n glig s par correction statique il est alors obligatoire de renseigner les mots cl s MODE CORR D FONC DTetD FONC DT2 9 MODE CORR modcor Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster ul Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 19 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 Concept de type mult_elas produit par la macro commande MACRO ELAS MULT U4 51 02 ou mode _meca qui correspond la r ponse
30. TMI O BASE ELAS FLUI meles melasflu NUME VITE FLUI Nvitf I O ETAT STAT NON DEFAUT OUI PREC DUREE 1 E 2 DEFAUT prec R 0 CHOC FLUI NON DEFAUT Z TOUT 0 NB MODE Nmode 1 NB MODE FLUI Nmodef I TS REG ETAB tsimu R INCREMENT F LIST INST litps listr8 PAS dt R INST_INIT ti R 0 INST FIN tf R NUME_FIN nufin I O VERI PAS P OUT DEFAUT NON Op randes sp cifiques une int gration par pas de temps adaptatifs d ordre 1 et 2 VITE MIN NORM DEFAUT MAXI Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Version default Titre Op rateur DYNA TRAN_MODAL Responsable Albert ALARCON ETAT_INIT Op randes et mots cl s sp cifiques l analyse sismique EXCIT RE _F MOD MOD E uN G Date 10 07 2015 Page 5 36 Fin des op randes et mots cl s sp cifiques l analyse sismique 0 CHOC Fi Cl U4 53 21 R vision 13175 0 COEF MULT PAS 1 1 DEFAUT cmp R 0 COEF DIVI_PAS 1 33333334 DEFAUT fred R 0 PAS
31. U4 53 21 R vision 13175 Les vecteurs g n ralis s doivent tre tablis par l op rateur PROJ_VECT BASE U4 63 13 partir de la base modale utilis e pour les matrices g n ralis es Dans le cas d un calcul par sous structuration dynamique les vecteurs g n ralis s doivent tre tablis par l op rateur ASSE VECT GENE U4 65 05 partir de la num rotation g n ralis e utilis e pour les matrices g n ralis es NUME ORDRE nmordr Num ro d ordre du mode d excitation de la structure Attention Il ne faut pas confondre le num ro d ordre du mode donn par le calcul modal dans l ordre o ils ont t calcul s et le num ro du mode intitul dans Code Aster NUME_ MODE 3 6 2 Op rande FONC MULT COEF MULT FONC MULT f Fonction du temps fonction permettant de d crire l volution temporelle du vecteur chargement COEF MULT a Coefficient multiplicateur du vecteur g n ralis valeur r elle constante par rapport au temps 3 7 Mot cl EXCIT_RESU Mot cl permettant de d finir un chargement sous la forme d une volution temporelle g n ralis e sans s paration de variables cas le plus g n ral Cette volution temporelle peut avoir t calcul e partir de l op rateur PROJ BASE option RESU GENE qui r alise la projection d un r sultat dynamique transitoire dyna trans 3 7 1 Mot cl RESULTAT RE
32. a que doter Li aldea ia ds 28 3 15 4 Op rande RELATION ranoni eera a a a a aoa eia 28 3 16 Mot cl RELA EFFO VITE 0000 ea dd Re a a e doo 29 3 17 Mot cl COUPLAGE EDY O S r a i aaa a aa a a aa a a a aada aa aiai aa aaaea 29 3 17 1 Op rande VITE ROTA 00d ida 29 3 17 2 0p rand PAS TPS EDO Sii ias 29 3 18 Mot cl PALIER EDYOS omita andan Rana deter 29 318 1 Op rande UNITE 300 30 3 18 2 Op rande GROUP NO Hama at asar 30 3 18 30p rande NOEUD ane cr meme A tt arr a 30 3 18 4 Op rande TYPE_EDYOS iii 30 3 19 Mot cl ARCHIVAGE states aladdin dto Le At asa 30 3 19 1 Op rande LIST_ARCH INST iii 31 3 19 2 Op rand PAS ARCH armies eatae A nent tennis 31 3 19 3 0p rande CRITERE 00 a o en tt ee NRA 31 3 19 4 Op rande PRECIS O N ie a e aaa aa Eaa a A aa OAAS 31 3 20 Op rand INFO oi A AA A id 32 3 21 Op rande IMPRESS O N nie lee tee dE a et feat ti nain tn dre tte 33 3 21 1 Op randes TOUT NIVEAU oococccccccccccccccinncnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnrnnnninnnnnnnnnnnnnnnnnnne 33 3 21 2 Op randes INST_INIT INST_FIN 00 ie 33 3 21 3 Op rande UNITE_DIS_VISC iii 33 3 22 Op rande TIRE nee nee nt Re RE ns de enr rs E 33 4 Phase d x cution anus entente re one rer tete li A A A dent One Rent 34 4 1 V rification sur les matrices 34 4 2 V rification et conseil sur le choix du pas de temps pour les sch mas EULER DEVOGE et
33. amma y angle en degr s d finissant la position angulaire du rep re local de l obstacle dans son plan Par convention la normale n au plan de coupe de l obstacle NORM OBST d finit l axe Xe du rep re local On passe du rep re global X Y Z au rep re du plan de l obstacle n y Z par un produit de deux rotations d angles autour de Z puis B autour du transform y de Y La position de l obstacle dans ce plan est obtenue par une rotation d angle B autour de la direction normale X e cf Figure 3 6 1 7 a Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Version Code Aster default Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 23 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 Y XI 22 vi x2 Zloc wa B_ Al Y dE Xi Z2 Yloc al v TE ls Bonzi Y Y2 Z Z1 Y1 Y2 X2 Xloc Figure 3 6 1 7 a Rotations permettant de passer du rep re global au rep re local de l obstacle Les angles a et B sont d termin s automatiquement partir de la normale l obstacle n Le rep re local X pe Yc Z se d duit ensuite du rep re M y Z par rotation d un angle de vrille ANGL_VRIL autour de n Remarque e Si l utilisateur ne pr cise rien l angle de vrille est calcul par le code dans le cas de chocs entre structures mobiles avec des obstacles de type BT PLAN
34. aut ona TS REG ETAB INST FIN INST INIT temps estim transitoire Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster peores Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 15 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 Remarques Ce sch ma d int gration n est pas utilisable en poursuite et ne permet pas le calcul par sous structuration dynamique La pr sence du mot cl CHOC est imp rative m me pour des simulations de phases sans chocs dites phases de vol 3 4 Mot cl INCREMENT 3 4 1 Op randes LIST_INST PAS VERI_PAS PAS MINI PAS MAXI LIST INST 1 temp Concept liste de r els de type listr8 Liste de r els d finissant les instants f de calcul de la solution e Sch mas RUNGE_KUTTA_54 et RUNGE_KUTTA32 Pour les sch mas de type Runge Kutta le mot cl LIST _INST n est pas pris en compte 7 PAS dt e Sch mas EULER DEVOGE NEWMARK Pas de temps du calcul transitoire e Sch mas ADAPT_ORDRE1 et ADAPT_ORDRE2 D signe le pas de temps initial utilis par l algorithme Ce param tre doit tre suffisamment faible e pour permettre le calcul des phases statiques qui utilise toujours le pas de temps maximal e pour d marrer correctement l algorithme Il doit
35. cependant tre suffisamment lev pour ne pas p naliser l ensemble du calcul e Sch ma ITMI D signe le pas de temps retenu pour le premier pas de calcul apr s passage ventuel du transitoire Par la suite l algorithme g re automatiquement le pas de calcul en fonction de la rigidit de la structure et des zones de transition vol choc Sch mas RUNGE_KUTTA_54 et RUNGE_KUTTA32 D signe le pas de temps initial propos par l utilisateur Si l erreur de pr diction entre les ordres de contr le v rifie err lt 1 alors il s agit du premier pas de calcul Sinon l algorithme choisit automatiquement le pas de temps n cessaire afin de v rifier cette condition Par la suite le choix du pas de temps dans les algorithmes de Runge Kutta est g r automatiquement 0 VERI PAS rep V rification du pas de temps de calcul relativement au pas de temps limite d termin en fonction de la fr quence la plus lev e des modes de la base modale consid r e ou des bases des sous structures Op randes sp cifiques une int gration par pas de temps adaptatifs avec les sch mas ADAPT ORDRE1 ADAPT ORDRE2 ainsique RUNGE KUTTA 54 et RUNGE KUTTA 32 0 PAS MAXI dtmax Valeur maximale du pas de temps Si les conditions d augmentation du pas de temps sont remplies le pas de temps courant pourra alors augmenter jusqu cette valeur limite Manuel d utilisation Fascicule u4 53 Analyse dynamique
36. enant la deuxi me sous structure Remarque On v rifie que les groupes de n uds contiennent bien un et un seul n ud GROUP MA OUMAILLE On peut aussi entrer les n uds de chocs en vis vis sous la forme de mailles SEG2 dessin es dans le maillage Ainsi on conserve une m me fa on de d crire les chocs que pour DYNA NON _ LINE avec les l ments discrets de choc DIS CHOC On peut entrer une liste de MAILLE ou de GROUP MA 3 10 1 3 Op rande OBSTACLE OBSTACLE obs Nom du concept de type obstacle d finissant la g om trie d un obstacle ind formable ou la forme enveloppe du jeu entre deux structures antagonistes ll est produit par l op rateur DEFI OBSTACLE U4 44 21 3 10 1 4 Op rande NORM OBST NORM OBST nor Liste de 3 r els d finissant la normale au plan de coupe de l obstacle c est dire le vecteur X we On conseille que Xe soit la direction de la fibre neutre ou d une g n ratrice de la structure tudi e 3 10 1 5 Op rande ORIG OBST ORIG OBST ori Liste de 3 r els d finissant la position de l origine de l obstacle dans le rep re global mot cl obligatoire dans le cas de chocs entre une structure mobile et une paroi fixe Dans le cas de chocs entre deux structures mobiles le code consid re par d faut que l origine est situ e au milieu des deux n uds de choc NOEUD 1 ou n ud du GROUP NO 1 et NOEUD 2 ou n ud du GROUP NO 2 3 10 1 6 Op rande
37. et d imprimer un ou plusieurs tableau x parmi DEPL LOC VITE LOC FORC LOC et TAUX CHOC Avec TOUT OUI valeur par d faut on imprime les quatre tableaux 3 21 2 Op randes INST_INIT INST FIN Ces deux mots cl s permettent l utilisateur de filtrer les impressions dans chaque boucle sur les pas de temps 3 21 3 Op rande UNITE DIS VISC 0 UNITE DIS VISC unit Les r sultats concernant l effort les d placements visqueux et relatifs entre les deux n uds ainsi que la dissipation du dispositif non lin aire peuvent tre sauvegard dans un fichier directement exploitable par les commandes de Code Aster 3 22 Op rande TITRE O TITRE titre Titre de la structure de donn es r sultat U4 03 01 Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster At Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 35 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 4 Phase d ex cution 41 V rification sur les matrices Dans le cas d un calcul par recombinaison modale on v rifie que les matrices g n ralis es sont bien issues d une projection sur une base commune et avec le m me nombre de vecteurs de base Dans le cas d un calcul par sous structuration dynamique on v rifie que les
38. ght 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Al Titre Op rateur DYNA TRAN MODAL Date 10 07 2015 Page 26 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 n T P i E o T reconstitution de l effort tranchant t oi F impo K On calcule ensuite les i 1 i 1 valeurs cumul es sur l ensemble des modes qui constituent la base modale utilis e On v rifie que le rapport de la souplesse n glig e souplesse statique moins souplesse statique reconstitu e sur la souplesse de choc reste inf rieur la valeur donn e par l op rande SEUIL SEUIL vaut 0 5 par d faut sinon e siSTOP CRITERE OUI on arr te l ex cution du programme c est le cas par d faut e si STOP CRITERE NON on continue l ex cution du programme avec mission d une alarme Remarques e Cette fonctionnalit n est disponible que pour des obstacles de type plan ou bi plan e Sile taux de reconstitution de la solution statique est inf rieur la valeur du seuil on conseille l utilisateur de compl ter la base modale par les modes locaux aux points de choc qui ont une souplesse locale importante e La formule n est pas applicable en cas de modes statiques matrice de rigidit non inversible Le calcul se poursuit alors sans v rification des crit res de choc et l utilisateur en est averti
39. ion de type Runge Kutta Dans ce cas il s agit du sch ma d int gration explicite de Bogacki Shampine 32 R5 06 04 pas de temps adaptatif Comme le sch ma pr c dent il supporte la prise en compte de l ensemble des non lin arit s disponibles dans l op rateur Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster ln Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 13 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 Pour ce sch ma le calcul du pas de temps optimal se fait par contr le de l erreur entre les approximations d ordre 3 et 2 de la pr diction du vecteur d tat Le calcul du pas de temps optimal quant lui se fait de mani re analogue au sch ma pr c dent 3 3 1 7 SCHEMA ADAPT_ORDRE1 sch ma explicite d ordre 1 Ce sch ma est une variante du sch ma pr c dent ADAPT ORDRE2 C est en fait une version du sch ma d Euler avec pas de temps adaptatif En dehors de cette diff rence ce sch ma s utilise de la m me mani re que le sch ma adaptatif d ordre 2 la syntaxe des mot cl s est la m me et les m thodes de pilotage du pas de temps aussi 3 3 1 8 SCHEMA ITMI sch ma int grale pour le calcul de la r ponse de syst mes m caniques tr s faiblement amortis avec couplages fluid lastiques Ce sch ma d int gration par m thode i
40. isateur La rigidit normale de choc apr s flambage k2 est ensuite diff rente de la rigidit avant flambage kn Seuls les op randes sp cifiques au mot cl FLAMBAGE sont d taill s Les autres mots cl s permettent de d finir les lieux de choc et sont identiques aux op randes du mot cl CHOC FNOR CRIT flim Force normale limite qui entra ne le flambage de la structure 0 FNOR POST FL fseuil Force normale limite apr s flambage qui provoque une d formation r siduelle de la structure RIGI NOR POST FL k2 Valeur de la rigidit normale apr s flambage 3 15 Mot cl RELA EFFO DEPL RELA EFFO DEPL Mot cl facteur permettant de d finir une relation force d placement ou moment rotation sur un degr de libert donn sous la forme d une courbe non lin aire 3 15 1 Op rande NOEUD NOEUD no Nom du n ud de la structure sur lequel porte la relation 3 15 2 Op rande SOUS_STRUC o SOUS STRUC ss Nom de la sous structure contenant le n ud renseignant l op rande NOEUD 3 15 3 Op rande NOM CMP O NOM CMP nomcmp Nom de la composante du n ud de la structure sur laquelle porte la relation 3 15 4 Op rande RELATION RELATION f Nom de la fonction non lin aire Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster ou Titre Op rateur DYNA TRA
41. ise en compte d un transitoire de vitesse de rotation 19 3 9 1 Op rande VITESSE _ VARIABLE iiiiiiiiiieeeeeeeeeeeneeeeeeeeeeeeeneeee 19 3 9 2 Op randes VITE ROTA MATR GYRO ACCE_ROTA et MATR_RIGY 20 3 10 Prise en compte de non lin arit s localis es de type choc frottement et rotor fissur 20 3 10 1 Non lin arit s localis es de type choc et frottement mot cl CHOC ocicnciciciccccncnno 20 3 10 2 Non lin arit s localis es de rotor fissur ROTOR FISS 24 3 11 Motek VER CHO Ces d 25 3 12 Mot Cl ANTES IS Mii A ia 25 3 13 Mot cl DIS VISC aeaa td tre an em eee den da da dd nn 25 II RS LE PC E A E E 26 3 13 2 Op randes li s la position du dispositif 26 3 13 3 Op randes li s au comportement 26 3 13 4 Op randes li s la convergence du comportement du dispositif 27 Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster E Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 3 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 3 14 Mot cl FLAMBA GE siitin onii a AA e A des due 27 3 19 M t cl RELA EFFO DEP E 200 a ido dada ed 28 31510p rande NOEUD 22 OR aa ae IRA aaa 28 3 15 2 0Op rande SOUS STRUC a a ter ne aleteo ici rt pierdes 28 3 15 3 Op rande NOM CMP a
42. ji ee pour le degr de libert i min 100 Peut tre utilis e sile syst me poss de un petit nombre de degr s de libert de 1 3 e pour un syst me plusieurs degr s de libert dans le cas o l ordre de grandeur du d placement est tr s diff rent selon les degr s de libert par exemple en pr sence de degr s de libert de Lagrange en sous structuration e __ si l ordre de grandeur de la vitesse ne varie pas trop au cours du temps 0 NMAX ITER PAS N Nombre maximal de r ductions du pas de temps par pas de calcul Il est par d faut gal 16 ce qui limite le coefficient de r duction du pas 0 75 10 par it ration lorsque le pas de temps est trop lev on reprend le calcul avec un pas plus faible Af 0 75 Af NMAX_ITER PAS peut tre e augment pour permettre au pas de temps de chuter de fa on plus brutale e diminu si le pas de temps semble excessivement raffin par exemple en pr sence de discontinuit s frottement sec excitation discontinue Si un instant donn on atteint ce nombre maximal de r ductions successives du pas de temps alors le code va quand m me consid rer que le pas final est correct et passer au pas suivant Un message d alarme est alors mis qui signale un ventuel risque de perte de pr cision et qui conseille l utilisateur de relancer le calcul avec des param tres modifi s en jouant sur PAS NMAX_ITER PAS et ou COEF DIVI
43. l op rateur PROJ_MATR_BASE U4 63 12 ou par la macro commande PROJ BASE U4 63 11 partir de la m me base modale Dans le cas d un calcul par sous structuration dynamique les matrices g n ralis es doivent tre tablies par l op rateur ASSE MATR GENE U4 65 04 partir de la m me num rotation g n ralis e MATR MASS ma Matrice de masse du syst me g n ralis Concept de type matr asse gene R MATR RIGI ri Matrice de rigidit du syst me g n ralis Concept de type matr asse gene R MATR AMOR am Matrice d amortissement du syst me g n ralis Concept de type matr asse gene R Cette option n est pas disponible avec la m thode DEVOGE 3 2 Mot cl AMOR MODAL Ce mot cl permet de prendre en compte un amortissement quivalent de l amortissement modal d compos sur une base de modes pr calcul e sous forme de concept de type mode meca Cet amortissement est globalement pris en compte dans l quation d quilibre dynamique comme une force correctrice au second membre C X 3 2 1 Op randes AMOR REDUIT LIST_AMOR 0 AMOR REDUIT ln Liste des amortissements r duits 7 n2 N pourcentages de l amortissement critique correspondants chaque mode du syst me sous forme de liste de r els Cette option n est pas disponible en sous structuration dynamique car les amortissements r duits doivent tre d finis pour chaque sous s
44. machines tournantes sur paliers non lin aires Le mot cl facteur PALIER EDYOS doit alors aussi tre renseign cf le paragraphe suivant 3 17 1 Op rande VITE_ROTA VITE ROTA vrota D finition de la vitesse de rotation du rotor en tours par minute 3 17 2 Op rande PAS TPS EDYOS PAS TPS EDYOS dtedyos Pas de temps initial pour le code EDYOS 3 18 Mot cl PALIER EDYOS 0 PALIER EDYOS Deuxi me mot cl facteur obligatoire lors du couplage avec le code EDYOS Il permet de d finir la position et le type des paliers non lin aires dont les r ponses seront calcul es par EDYOS Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster en Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 31 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 Les donn es pr ciser pour chaque palier non lin aire sont e position du palier soit par le nom du n ud du maillage soit par le nom du groupe contenant ce n ud seul e type de palier choisi correspondant donc un mod le EDYOS particulier Ces donn es peuvent tre entr es suivant deux modes de fonctionnement Soit toutes les donn es sont dans un fichier texte externe dont l unit logique sera pr cis e par le mot cl UNITE Soit on rentre les donn es dans le fichie
45. matrices g n ralis es sont bien issues de la m me num rotation g n ralis e 4 2 V rification et conseil sur le choix du pas de temps pour les sch mas EULER DEVOGE et NEWMARK On s assure que le pas de temps choisi v rifie les conditions de stabilit du sch ma num rique crit re de CFL e dans le cas de NEWMARK la stabilit est toujours assur e mais le d passement du crit re peut induire un manque de pr cision sur le r sultat et est signal par un message le calcul se poursuit au risque de produire un r sultat peu pr cis ou faux e dans le cas des sch mas d EULER et DEVOGE si l op rande VERI PAS vaut OUI valeur par d faut l ex cution est arr t e un pas de temps minimum est propos Si l op rande VERI_PAS vaut NON ou s il s agit d un sch ma adaptatif un message d alarme est mis et le calcul se poursuit au risque de produire un r sultat peu pr cis ou faux Dans une analyse transitoire sans non lin arit il faut veiller ce que le pas de temps soit tel que dt lt 0 1 f pour NEWMARK et DEVOGE dt lt 0 05 f pour EULER f tant la fr quence la plus lev e des modes de la base modale consid r e Remarque On mentionne qu avec des non lin arit s localis es le pas de temps choisi doit tre parfois tr s inf rieur cette valeur conseill e 4 3 Phase d ex cution pour les m thodes ADAPT ORDREl et ADAPT_ORDRE2 AE L ex cution est interrompue lorsque le pas de temps atteint un
46. nt grale permet pour les syst mes faiblement amortis d obtenir une r ponse exacte en tenant compte des variations de forces fluid lastiques obtenues en pr sence de chocs On d crit ci dessous les mots cl s sp cifiques au calcul de la r ponse de syst mes m caniques lin aires tr s faiblement amortis avec couplages fluid lastiques associ s ventuellement des non lin arit s localis es aux n uds de type chocs et frottements 0 BASE ELAS FLUI meles Base modale utilis e pour le calcul Concept de type melasflu produit par l op rateur CALC FLUI STRU U4 66 02 qui contient l ensemble des bases modales calcul es pour les diff rentes vitesse d coulement d finies Ce mot cl est obligatoire pour la m thode ITMI Le calcul transitoire sur base modale modifi e par le couplage fluid lastique s effectue en prenant en compte les valeurs des amortissements ajout s dus l coulement du fluide qui sont pr sents dans le concept mlasflu d entr e Les amortissements modaux r cup r s de la base fluid lastique remplacent ceux renseign s sous le mot cl global AMOR REDUIT de l op rateur DYNA TRAN MODAL NUME VITE FLUI Nvitf Vitesse d coulement retenue pour le calcul num ro d ordre Permet d extraire dans le concept melasflu la base modale correspondant la vitesse d coulement retenue cf U4 66 02 Ce mot cl est obligatoire pour la m thode ITMI ETAT STAT
47. ononnnnnnncnnnnnnnnnnnnnnnnnnononnnnnnnnnnnnnaannnns 9 3 3 Sch mas d int gration Mot cl SCHEMA_ TEMPS iii 9 3 3 1 Op rande SCHEMA cococooccccnccnonooonnnnnoncnnonnnnnnononnnnnnnnnnnnnnnnn nn nan nnnnn nn nn nn nr nnnn nn nn 9 3 4 Mot cle INCREMEN Poca a A be daa 14 3 4 1 Op randes LIST_INST PAS VERI_PAS PAS_MINI PAS_ MAXI 14 3 4 2 Op randes INST_INIT INST_ FIN NUME_ FIN iii 16 3 5 Mot cl ETAT IN iso a oda 16 3 5 1 Op randes RESUT DEPL VITE erara a e dd ai 16 3 5 2 Op randes INST_INIT NUME_ORDR occcccccooooncnncononononcnnccoononnnnonccononononnncconannnnnnncnnnnnnos 16 3 5 3 0p rande CRITERE setii a e tos ruedo O o do 17 3 5 4 0perande PRECISO No ui Aia 17 3 6 Description du chargement sous variables s par es mot cl EXCIT 17 3 6 1 Op randes VECT_ASSE_GENE NUME_ORDRE iii 17 3 6 2 Op rande FONC_MULT COEF_MULT oooooooccncccnonononcnnccnononnnnnnccnnnnnnnnnnncnnnnnnnnnnncccnnnnninnns 17 3 1 Mot cle EXCIT RESU ua a o barides 18 3 1 Motcl RESULTA Tic ts 18 3 1 2 Op erande COEFR MULTA indie at ein 18 3 8 Cas particulier de l analyse sismique 18 3 8 1 Prise en compte des modes n glig s par correction statique mots cl s CORR STAT MODE CORR et D FONO rascar ere rte re dames nee ter est 18 3 8 2 Prise en compte du multi appuis mots cl s MODE _STAT MULTI APPUI et ACCE VITE DEPL ronnan aneta os te o o dedo o o ts 19 3 9 Pr
48. orces de choc correspondantes Sinon MULT APPUI NON on compare chaque instant le vecteur des d placements relatifs de chacun des n uds susceptibles de choquer 0 EXCIT F ACCE ac 0 VITE vi 0 DEPL dp Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster dl Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 20 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 Noms des fonctions acc l ration ACCE vitesse VITE et d placement DEPL impos es lors du calcul de la r ponse sismique de structures multi support es Remarque Si la structure est mono excit e l acc l rogramme est d fini par le mot cl FONC_MULT EXCIT F DIRECTION dx dy dz drx dry drz Composantes du vecteur donnant la direction du s isme dans le rep re global EXCIT Fi NOEUD lno GROUP_NO lgrno Liste des noms de n uds ou de groupe de n uds correspondants aux appuis concern s o le s isme est impos EXCIT F VECT ASSE GENE v Vecteur projet de l excitation sismique issu de CALC CHAR SEISME U4 63 01 3 9 Prise en compte d un transitoire de vitesse de rotation 3 9 1 Op rande VITESSE VARIABLE Pr cise si la vitesse de rotation du rotor est variable en fonction du temps VIT
49. pas minimal gal PAS X PAS LIMI RELA Remarques Le sch ma des diff rences centr es ne restitue pas de fa on exacte les pulsations propres d un Syst me ce qui conduit d importantes erreurs de calcul dans les deux cas suivants e Calcul d un tr s grand nombre de p riodes d oscillations libres gt Calcul des oscillations d un syst me tr s faiblement amorti lt 10 excit sur une fr quence de r sonance Dans ces deux cas il est souvent n cessaire d augmenter le param tre NB POIN PERIODE Les m thodes ADAPT ORDREl et ADAPT ORDRE2 peuvent tre utilis es en sous structuration Le pas de temps peut tre r cup r par l op rateur RECU FONCTION avec la syntaxe suivante pas RECU FONCTION RESU GENE dynamoda NOM CHAM PTEM Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster dl Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 36 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 4 4 Phase d ex cution pour la m thode ITMI L ex cution est interrompue e lorsque la dur e d excitation choisie par l utilisateur est incompatible avec le temps de simulation souhait r gime tabli simulation apr s obtention du r gime tabli Dans ce cas l utilisateur en est inform avec pr cision du temps minimum d excitation requis pour
50. ptatifs nombre de points pas p riode nbp coefficient de remont e du pas de temps cmp coefficient de division du pas de temps cdp pas de temps minimal dtmini pas de temps maximal dtmaxi nombre maximal de r ductions du pas nbred vitesse minimale variable gt m thode de calcul de la vitesse de r f rence lt fin des informations sp cifiques aux sch mas adaptatifs nombre de pas de calcul est nbe Le pas de temps du calcul est valeur du pas de temps de calcul le nombre de RELA EFFO DEPL est nbrelaed uniquement si le nombre de relations est non nul le nombre de RELA EFFO VITE est nbrelaev uniquement si le nombre de relations est non nul Si INFO 2 on imprime en plus des informations crites dans le cas o INFO vaut 1 les informations suivantes dans le fichier MESSAGE Pour chaque obstacle eLe num ro et type de l obstacle Le nom et les coordonn es dans le rep re global du n ud de choc des n uds de choc dans le cas d un choc entre structures mobiles L orientation dans le rep re global de la normale l obstacle La valeur de l angle de vrille La valeur du jeu initial Et pour chaque n ud de choc et pour chaque mode le num ro du mode les valeurs des raideurs locales de choc et du taux de flexibilit locale et de la souplesse locale On imprime galement la fin pour chaque n ud de choc TAUX DE RESTIT FLEXIBILITE 9 9539E 0
51. r de commande directement via GROUP NO ou NOEUD et TYPE EDYOS 3 18 1 Op rande UNITE UNITE uled D finition de l unit logique du fichier texte externe qui contient la d finition des positions et des types de tous les paliers non lin aires Voici un exemple de ce fichier 2 1 NOEO PACONL 2 NOE10 PACONL La premi re ligne indique le nombre total de paliers non lin aires ici deux Ensuite chaque ligne correspond un palier avec trois param tres num ro arbitraire de palier position du palier nom du n ud ou du groupe de n ud correspondant et mod le de palier EDYOS ici PACONL Pour le deuxi me argument le code va chercher dans tous les groupes de n uds du maillage puis si le nom indiqu n est pas trouv la recherche se poursuit dans les noms de n uds et en cas d chec un message d erreur est mis et pr cise que le nom indiqu par l utilisateur n a pas pu tre trouv Le code s arr te alors en erreur fatale 3 18 2 Op rande GROUP _NO GROUP NO grnoed Groupe de n ud ne contenant qu un seul n ud et qui correspond la position du palier non lin aire consid r 3 18 3 Op rande NOEUD NOEUD noed Nom du n ud qui correspond la position du palier non lin aire consid r On conseille classiquement de privil gier l utilisation de GROUP NO NOEUD pour des raisons de nommage plus pratique 3 18 4 Op rande TYPE EDYOS TYPE EDYOS PAPAN
52. re peuvent tre sauvegard s dans un fichier directement exploitable par les commandes de Code Aster La d finition du fichier se fait par le mot clef simple UNITE DIS VISC quiest sous le mot clef facteur IMPRESSION de la commande 2 3 13 1 Syntaxe 0 DIS VISC F NOEUD 1 nol noeud GROUP NO 1 grnol group no NOEUD 2 no2 noeud GROUP NO 2 grno2 group no K1 Ki R UNSUR K1 uskl R K2 k2 R UNSUR K2 usk2 R K3 k3 R UNSUR K3 usk3 R Cc c R PUIS ALPHA 7 045 defaut alpha R 9 ITER INTE MAXI 20 defaut iter 1 RESI INTE RELA 1 0E 06 defaut resi R 3 13 2 Op randes li s la position du dispositif NOEUD 10u GROUP NO 1 NOEUD 2 OU GROUP NO 2 N uds ou nom des groupes de n ud de la structure entre lesquels le dispositif non lin aire est plac Si on utilise un GROUP NO pour d finir l un des points d accrochage du dispositif le groupe de n ud ne doit contenir qu un seul n ud Lors du calcul il est n cessaire de conna tre la direction du dispositif non lin aire car il ne fonctionne que dans son axe faut donc que la distance entre les deux n uds soit non nulle Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster ui
53. renseignant le mot cl NOEUD 2 ou GROUP_NO_2 3 10 1 10 Op rande REPERE 0 REPERE rep Pr cise le rep re dans lequel la position de l obstacle est d finie GLOBAL La position absolue de l obstacle est d finie ind pendamment des rotations et translations auxquelles sont soumises les diff rentes sous structures nom sst Nom d une sous structure La position et la normale de l obstacle sont d termin es dans le rep re utilis pour d finir les coordonn es des n uds de la sous structure nom sst la position et la normale finales de l obstacle tant le r sultat de la rotation et de la translation auxquelles est soumise la sous structure 3 10 1 11 Op rande RIGI_NOR RIGI NOR kn Valeur de la rigidit normale de choc unit N m en USI 3 10 1 12 Op rande AMOR NOR AMOR NOR cn Valeur de l amortissement normal de choc unit N m s en USI 3 10 1 13 Op rande RIGI_ TAN RIGI TAN kt Valeur de la rigidit tangentielle de choc unit N m en USI 3 10 1 14 Op rande AMOR TAN AMOR TAN ct Valeur de l amortissement tangentiel de choc unit N m s en USI Remarque Si une raideur K est sp cifi e et que le mot cl AMOR TAN est absent le code calcule un amortissement optimis de fa on minimiser les oscillations r siduelles en adh rence selon la formule Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Documen
54. rt ALARCON Cl U4 53 21 R vision 13175 faudra veiller introduire une dur e d excitation suffisamment longue Cette dur e doit correspondre la dur e n cessaire au passage du transitoire augment e de la dur e de simulation en r gime tabli souhait e Cette dur e totale de simulation sera renseign e via les deux op randes INST_INIT et INST FIN sous le mot cl facteur INCREMENT ETAT STAT NON La simulation ne distingue pas l tat transitoire du r gime tabli PREC DUREE prec Permet de d finir la pr cision choisie pour d terminer la dur e de la phase transitoire selon la formule n prec n e PC Ge o amp et w d signent respectivement l amortissement r duit et la pulsation o Wo de chaque mode consid r La valeur par d faut de ce param tre est 1 CHOC_FLUI D termine le traitement r alis par l algorithme pendant les phases de choc vis vis des forces fluid lastiques Par d faut la variation des forces fluid lastiques en phase de choc li e la modification de la rigidit et de l amortissement du syst me m canique impact sur la but e n est pas prise en compte NB MODE Nmode Nombre de modes de la base modale retenus pour le calcul dynamique Les modes conserv s correspondent des fr quences croissantes premiers modes Si NB_MODE n est pas pr cis on prend tous les modes de la base modale du concept de type melasflu
55. s de l amortisseur seul Op randes li s a la convergence du comportement du dispositif O ITER INTE MAXI 20 defaut ater I 0 RESI INTE RELA 1 0E 06 defaut resi R Ces op randes ont la m me signification que lorsqu ils sont utilis s avec la commande STAT_NON_LINE COMPORTEMENT U4 51 11 La relation de comportement DIS_VISC n cessite de r soudre un syst me non lin aire par une m thode de Runge Kutta d ordre 5 pas adaptatif Le contr le de l algorithme nombre d it ration et r sidu sont utilis s pour tester la convergence et adapter le pas si besoin Mot cl FLAMBAGE Ce mot cl est utilis pour la d tection de flambage ventuel et pour l valuation de la d formation r siduelle d un l ment lors d un choc entre deux structures mobiles ou entre une structure mobile et une paroi fixe La force de r action lors d un choc avec prise en compte du flambage peut tre r sum e par le sch ma suivant Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster MA Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 29 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 F Flim Fseuil compression On consid re qu il y a flambage si la force de r action F atteint la valeur limite F d finie par l util
56. sal NB MODE FLUI Nmodef Nombre de modes de la base modale perturb s par les ph nom nes de couplage fluid lastiques en phase de choc inf rieur au nombre de modes retenus pour le calcul dynamique Les modes conserv s correspondent aux Nmodef premi res fr quences croissantes premiers modes Si NB MODE FLUI n est pas pr cis on prend le nombre de modes retenus pour le calcul dynamique TS REG ETAB tsimu Dur e de la simulation souhait e Dans le cas d une simulation sans calcul pr alable et en un pas de temps de la phase transitoire ETAT STAT NON cette dur e correspond la dur e de simulation quel que soit l tat du syst me entre les instants de d but et de fin de simulation Par cons quent on devra s assurer que TS REG ETAB lt INST FIN INST INIT Par d faut on aura TS REG ETAB INST FIN INST INIT Dans le cas d une simulation avec calcul de la phase transitoire ETAT STAT OUI cette dur e correspond la dur e de simulation r ellement souhait e lorsque la phase de chocs est tablie du point de vue num rique Par cons quent on devra s assurer que TS REG ETAB lt INST FIN INST INIT temps estim transitoire Dans le cas o cette derni re condition n est pas respect e l utilisateur en est inform avec pr cision du temps minimum d excitation requis pour son calcul INST_ FIN INST INIT Par d f
57. son calcul e lorsque l algorithme ne r ussit pas trouver une solution converg e lors de la diagonalisation de la matrice de raideur e lorsque les phases de transition vol choc ne peuvent tre d termin es avec une pr cision suffisante Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html
58. statique lin aire de la structure un chargement unitaire de type force impos e acc l ration uniforme dans la direction du s isme consid r e On note qu il y a autant de cas de charge que de direction de s isme EXCIT F CORR STAT Si MODE CORR est pr sent CORR_STAT 0UI permet de prendre en compte la contribution de la correction modale a posteriori pour chaque occurrence du mot cl EXCIT EXCIT F D FONC DT et D FONC DT2 D_FONC DT et D _FONC DT2 sont respectivement les d riv es premi res et d riv es secondes du temps de l acc l rogramme d fini dans chaque direction sismique consid r e par l op rande FONC MULT Elles pond rent la contribution de la correction modale a posteriori pour chaque occurrence du mot cl EXCIT afin d obtenir respectivement les corrections de vitesse et d acc l ration sur la base physique Remarques e La prise en compte de la correction statique exclue celle du multi appuis e Le concept mult elas doit s appuyer sur une num rotation des quations coh rente m me profil et m me option de renum rotation avec celle du syst me r solu dans l op rateur DYNA TRAN MODAL e A la i me occurrence du mot cl EXCIT correspond la i me solution lastique de MODCOR e Pour que la correction statique soit effectivement prise en compte lors du retour vers les coordonn es physiques dans l op rateur REST_GENE_ PHYS ou RECU FONCTION il faut pr ciser CORR
59. t diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster peris Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 25 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 c 24 k k m 2E k m o i est l indice du mode pr pond rant dans la r ponse de la structure 3 10 1 15 Op rande FROTTEMENT FROTTEMENT NON La condition de contact est sans frottement COULOMB COULOMB mu Valeur du coefficient de frottement sans dimension COULOMB STAT DYNA COULOMB STAT mus Valeur du coefficient d adh rence sans dimension COULOMB DYNA mud Valeur du coefficient de glissement sans dimension 3 10 2 Non lin arit s localis es de rotor fissur ROTOR FISS Les op randes suivants sont sp cifiques au calcul transitoire avec non lin arit localis e de type rotor fissur pour les calculs de ligne d arbre mod lis e en 1D poutre La fissure est consid r e totalement incluse dans une section du rotor Elle est d limit e par deux n uds distincts mais de coordonn es confondues l un reli e la partie gauche de la ligne d arbre l autre la partie droite Ils repr sentent respectivement la l vre gauche et la l vre droite de la fissure Le comportement de la fissure est donn par une loi de raideur de fissure et sa d riv e Cette loi est d termin e par ailleurs par des calculs 3D en quasi statique Elle ne d
60. tation U4 44 21 de l op rateur DEFI OBSTACLE d une autre structure antagoniste 3 10 1 1 Op rande INTITULE 0 INTITULE int Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster E Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 21 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 Intitul huit caract res au maximum permettant de nommer la non lin arit Si rien n est pr cis par l utilisateur l intitul est le nom du NOEUD _1 3 10 1 2 Op randes NOEUD_1 NOEUD_2 GROUP_NO_1 GROUP_NO_2 GROUP_MA NOEUD 1 ou GROUP NO 1 Noeud ou nom du groupe de n ud de la structure sur lequel porte la condition de non lin arit Dans le cas d un calcul non lin aire par sous structuration dynamique on indique sous ce mot cl le n ud de choc appartenant la premi re sous structure les diff rentes sous structures n appartiennent pas au m me maillage NOEUD 2 ou GROUP NO 2 N ud ou nom du groupe de n ud de la seconde structure sur lequel porte la condition de non lin arit Cette op rande est sp cifique la d finition d un contact entre deux structures mobiles Dans le cas d un calcul non lin aire par sous structuration dynamique on pr cise le n ud de choc coincidant avec le n ud indiqu dans NOEUD 1 ou GROUP NO 1 mais appart
61. tructure s par ment op rateur MACR ELEM DYNA U4 65 01 Remarque Si le nombre d amortissements r duits donn s est inf rieur au nombre de vecteurs de base utilis s dans la base modale les amortissements des vecteurs suppl mentaires sont pris gaux au dernier amortissement de la liste LIST AMOR c formule n Nom du concept de type 1istrg8 contenant la liste des amortissements r duits 3 3 Sch mas d int gration Mot cl SCHEMA TEMPS Sous ce mot cl on peut renseigner un sch ma d int gration avec ventuellement ses param tres Les sch mas disponible sont d clarer sous l op rande SCHEMA Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster ou Titre Op rateur DYNA TRAN_MODAL Date 10 07 2015 Page 10 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 3 3 1 Op rande SCHEMA 0 SCHEMA Choix de la m thode num rique de r solution Dans le cas d un calcul classique par recombinaison modale Putilisateur dispose de six m thodes de type explicite d une m thode int grale et d une m thode de type implicite Dans le cas d un calcul par sous structuration dynamique R4 06 04 la m thode de calcul transitoire sur base modale calcul e par sous structuration supporte tous les sch mas d int gration voqu s sauf la m thode int grale En
62. ue avec quelle pr cision la recherche de l instant archiver doit se faire Manuel d utilisation Fascicule u4 53 Analyse dynamique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code_Aster default Titre Op rateur DYNA_TRAN_MODAL Date 10 07 2015 Page 33 36 Responsable Albert ALARCON Cl U4 53 21 R vision 13175 3 20 Op rande INFO INFO imp Entier permettant de pr ciser le niveau d impression dans le fichier MESSAGE Si INFO 1 on imprime les informations suivantes dans le fichier MESSAGE lt I gt lt nom de la routine o sont crites les informations suivantes gt Si lt I gt lt MDTR74 gt on rappelle que c est un calcul transitoire sur base modale classique sinon lt I gt lt SSDT74 gt c est un calcul transitoire sur base modale par sous structuration dynamique lt gt Calcul par superposition modale la base de projection est un SNPE NE T PASE OS PRIE le nb d quations est nb la m thode utilis e est gt nom de la m thode d int gration lt la base utilis e est gt nom de la base modale lt le nb de vecteurs de bas st nbb le pas de temps initial est valeur du pas de temps initial le nb de pas d archive est nba le nombre de lieu x de choc est nbchoc d but des informations sp cifiques aux sch mas ada
Download Pdf Manuals
Related Search