[U4.53.02] Opérateur DYNA_LINE_TRAN
Contents
1. PRECISION 1 E 06 DEFAUT prec R Indique avec quelle pr cision la recherche de l instant doit se faire 4 7 Mot cl EXCIT 0 EXCIT Op rande permettant de d finir plusieurs excitations spatio temporelles Soit en indiquant un vecteur assembl correspondant un chargement soit des charges qui conduiront au calcul et l assemblage d un second membre Le vecteur assembl peut tre associ une fonction d volution temporelle ou un coefficient multiplicateur constant Le chargement total est la somme des chargements d finis par toutes les occurrences du mot cl EXCIT cf 84 7 2 4 7 1 Op randes VECT_ASSE CHARGE VECT ASSE vecti Vecteur assembl correspondant un chargement concept de type cham no DEPL R COEF MULT ci Coefficient multiplicatif du vecteur assembl vecti FONC _ MULT Q Voir 4 7 2 CHARGE chi chi est le chargement comportant ventuellement l volution d un champ de temp rature pr cis par la me occurrence de EXCIT Voir 4 7 2 4 7 2 Op rande FONC_MULT 0 FONC_MULT Q a est la fonction du temps multiplicative du vecteur assembl ou du chargement pr cis la L me occurrence de EXCIT Le chargement ch et les conditions aux limites pour y occurrences du mot cl facteur EXCIT sont n ch gt a t ch i 1 Le ou les champs de temp rature ne sont pas multip2. Op rateur DYNA LINE TRAN Date 05 04 2013 Page 13 16 Responsable Albert ALARCON Cl U4 53 02 R vision 10811 Pour le sch ma adaptatif il d signe la fois le pas de temps initial et le pas de temps maximal utilis s par l algorithme Ce param tre doit tre suffisamment faible e pour permettre le calcul des phases statiques qui utilisent toujours le pas maximal e pour d marrer correctement l algorithme Il doit cependant tre suffisamment lev pour ne pas p naliser l ensemble du calcul 4 12 2 Op randes INST_INIT INST_FIN NUME_FIN Pour les sch mas des diff rences centr es et pas de temps adaptatif 9 INST INIT ti En cas de reprise on utilise le mot cl ETAT_INIT 84 6 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 sous INCREMENT doit donc tre utilis e uniquement s il n y a pas reprise d un calcul pr c dent 9 INST FIN tf Instant de fin du calcul transitoire Obligatoire pour les sch mas des diff rences centr es et pas de temps adaptatif NUME FIN nufin Num ro de l instant de fin de calcul dans LIST_INST uniquement pour les sch mas de Newmark et Wilson Si INST_INIT n est pas pr sent l instant initial est z ro 4 12 3 Op randes VITE MIN COEF MULT_PAS COEF_DIV_PAS PAS LIMI RELA NB_POIN_PERIODE NMAX ITE
3. gt 1 lorsque l erreur est suffisamment faible 0 75 At lt depuis plus de 5 pas cons cutifs gt A t min cmp Af A t pax Nf AP avec At_ Af max initial n 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 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 NMAX ITER PAS n est pas atteint et que le pas de temps minimal n est jes atteint f At At ne Niter lt Niter max et AT pir xAt At 7 AP cdp Sa valeur par d faut est de 1 3334 soit une r duction d un facteur 0 75 9 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 A min plr A initial NB POIN PERIODE N Nombre de points par p riode apparente C est ce param tre qui fixe la pr cision du calcul doit tre au moins gal 20 sa valeur par d faut 50 garantit une pr cision satisfaisante de l ordre de 1 2 dans la plupart des cas NMAX ITER PAS Nombre maximal de r ductions du pas de temps par pas de calcul si err gt l et Nre lt N remax At cdpxAt Il est par d faut gal 16 ce qui limite le coefficient de r duction du pas 171 33 10 par
4. 00 A a 10 4 9 10p randes MODE_MECA AMOR_REDUIT NB_MODE 10 4 10Motcl ENERGIE hs a o e e o n ba ak ai 10 4 11Moteel SOLVEUR 55 asa Laa Ma aa a a lat tn aaa a La a ns 11 4 12Motcl INCREMENT vaciones dana a an a a ga a ange sa Ka a alan di 11 4 12 10p randes LIST INST PAS aaa eaaa eaaa eaaa aaa aana iiiiiiiiiiesseseeeeeeeeseeeeeeeee 11 4 12 20p randes INST_INIT INST_FIN NUME FIN 12 4 12 30p randes VITE MIN COEF MULT PAS COEF DIV PAS PAS LIMI RELA NB_POIN_ PERIODE NMAX_ITER_PAS PAS MINI 12 4 13Mot cl ARCHIVAGE tab ee te und te ss fete ln ht ak En net LS 13 4 13 10p randes LIST_INST INST a a a e ag dd aa na aaa a aana A na St Nga 14 4 13 20p rande PAS ARGE san apa awa Bana nana aa baa A Tan a a Mg it 14 4 13 30p rande CRITERE 2 28 000 a anah aa a E a ag Sa a GA da 14 4 13 40p rande PRECISION ooooooocccccccccncccncnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnncnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnanicines 14 4 13 50p rande CHAM EXCLU asasaran saa a aa a a aaa ga ae a ga an aas eka aga a saa nan naa 14 4 140p rand TITRE kk aa b naak nada NGA Nak a ge tt et a Da saba A aa b aa a abab pa aa Da Aak ng a 14 EoI pood produits sense aaa 14 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
5. MODELE coi nr nur Minnie hein ab bawa atan a te 6 4 20p rande CHAM MATER sa aan Ba a aan Mann ne Me ne AA ee Me a E a 6 4 3Op rande CARA ELEM 5 35 agan ad tt Laa sadana a Eu mare meme Aa a eee as A Na a Ta 6 4 4Matrices du Probl me siba aa aa A aaa aa a ga Gad a a aa aa a aa NUN iaa 6 4 5Sch mas d int gration Mot cl SCHEMA_ TEMPS 7 4 5 10p rande SCHEMA coil en dre na element na MEA gg E Te teen nn ad 7 4 GMot Cl ETAT AN Tsani aab aa a trie nl nai im ne me nn ann Reda en 8 4 6 16p randes RESULTA Titus A aire nid a lan da TA de 8 4 6 20p randes DEPL VITE ACCEoooocccccccccnnncnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnininnns 8 4 6 30p randes NUME_ORDRE INST_INlT ooooocccccccccnncnncnnnnnnnnnnnnnnnnnnnnconcnnnnnnnnnnnnnnnnnnnnnnncinncnns 8 4 6 40p rande CRITERE m e 5546 edak eaaa ea bana ag aa aga a a a gra aa ke gawakna ea baa aga a aaa E a ataga 8 4 6 50p rande PRECISION ii aa aana anana aana nana aana aana near 9 4 7Motcl EXQI Taraa ah A AA a a aa ag a aaa dane baa a saanak aaa a 9 4 7 10p randes VECT_ASSE CHARGE anana aana anana nana nana a aana anana ena aaa 9 4 7 20p rande FONC MU TT 5 sisng aa ap panak ak ah a PE naa a gag paha ada sa ka E ka ban aaa na NG nahan ag a a a 9 4 7 30p randes MULTI_APPUI ACCE VITE DEPL DIRECTION NOEUD GROUP _ NO MODESSTA Tica a Ii A AA a dada daa NA daa 10 4 8Mot cl EXCIT RESU LAA a a a a gae 10 4 9Motecl AMOR MODAL
6. PARA CALC Le bilan d nergie peut tre extrait de cette table l aide de la commande RECU TABLE U4 71 02 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 LINE TRAN Date 05 04 2013 Page 12 16 Responsable Albert ALARCON Cl U4 53 02 R vision 10811 4 11 Mot cl SOLVEUR La syntaxe de ce mot cl commun plusieurs commandes est d crite dans le document U4 50 01 Cependant il convient de faire attention quelques singularit s de l utilisation des solveurs dans DYNA LINE TRAN La renum rotation RENUM RCMK utilisable pour LDLT et GCPC modifie l ordre des inconnues NUME_EQUA alors que la renum rotation pour MULT_FRONT METIS MDA MD est interne au solveur les inconnues des champs solutions sont rang es dans l ordre naturel des n uds du maillage Or dans DYNA LINE TRAN les matrices tant fournies en entr e elles ont d j une num rotation SANS ou RCMK et on ne pas toujours changer cette num rotation dans DYNA LINE TRAN On r capitule ci dessous les diff rents cas de figure 1 si matrices K et M sont num rot es avec METHODE MULT FRONT RENUM METIS MDA MD ou METHODE LDLT RENUM SANS alors on pe
7. 4 13 2 Op rande PAS ARCH PAS ARCH ipa Entier d finissant la p riodicit d archivage de la solution du calcul transitoire dans le concept r sultat dyna_ trans Si ipa 5 on archive tous les 5 pas de calcul 4 13 3 Op rande CRITERE 9 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 4 13 4 Op rande PRECISION PRECISION 1 E 06 DEFAUT prec R Indique avec quelle pr cision la recherche de l instant archiver doit se faire 4 13 5 Op rande CHAM_EXCLU 9 CHAM EXCLU l DEPL l VITE ACCE LA Permet d exclure l archivage d un ou plusieurs champs parmi DEPL VITE et ACCE Cette exclusion est ignor e pour le dernier instant de calcul les trois champs sont n cessaires pour une POURSUITE 4 14 Op rande TITRE 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 copylett fdl html Code Aster dl Titre Op rateur DYNA LINE TRAN Date 05 04 2013 Page 16 16 Responsable Albert ALARCON Cl U4 53 02 R v
8. A LINE TRAN Elle permet aussi de d finir des conditions initiales de type champs aux noeuds Nota bene Pour les sch mas d ordre sup rieur NEWMARK ou WILSON l acc l ration initiale acce init joue un r le important dans l initialisation du sch ma 4 6 1 Op randes RESULTAT RESULTAT dy Concept de type dyna trans issu d un calcul pr c dent avec DYNA LINE TRAN et d finissant les conditions initiales pour le nouveau calcul 4 6 2 Op randes DEPL VITE ACCE DEPL do Concept correspondant aux d placements initiaux champ aux noeuds de grandeur DEPL R VITE vo Concept correspondant aux vitesses initiales champ aux n uds de grandeur DEPL R ACCE ao Concept correspondant aux acc l rations initiales champ aux n uds de grandeur DEPL R Si le mot clef est pr sent on utilise le champ d acc l ration entr pour initialiser les diff rents sch mas d int gration en temps selon les algorithmes d crits dans le document R5 05 02 S il est absent on calcule une acc l ration initiale par la formule suivante M ao Fext t to C vo K xo Remarque importante Lorsque l tat initial du syst me dynamique est d fini par des champs de DEPL VITE et ou ACCE les composants de ces champs qui n ont pas t explicitement renseign s lors de la cr ation des champs sont consid r s nuls lors du calcul dynamique transitoire 4 6 3 Op randes NUME
9. Code Aster Co Titre Op rateur DYNA LINE TRAN Date 05 04 2013 Page 1 16 Responsable Albert ALARCON Cl U4 53 02 R vision 10811 Op rateur DYNA LINE TRAN 1 But Calculer la r ponse dynamique transitoire une excitation temporelle quelconque Le chargement temporel doit tre donn sous la forme d une combinaison lin aire de vecteurs forces assembl s constants dans le temps Ils peuvent tre fournis directement sous forme de vecteurs assembl s ou sous forme de charges qui seront assembl es dans l algorithme Seuls les coefficients de la combinaison lin aire sont fonction du temps Les m thodes d int gration implicites disponibles sont WILSON theta et NEWMARK et les m thodes d int gration explicites disponibles sont le sch ma aux diff rences centr s et une version pas adaptatif de ce m me sch ma Des instants d archivage peuvent tre sp cifi s Produit un concept r sultat de type dyna trans 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 LINE TRAN Date 05 04 2013 Page 2 16 Responsable Albert ALARCON Cl U4 53 02 R vision 10811 Table des mati res B tan aaea o aio a a ale 1 SN PE aa AT PG aaa a E a Ta Kampak ETN 4 3Equations du comportement sous excitation transitoire 6 LO A Aa 6 4 10p rande
10. DYNA LINE TRAN Responsable Albert ALARCON 6Phase d ex cution 7Bibliographie Date 05 04 2013 Page 3 16 Cl U4 53 02 R vision 10811 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_LINE_TRAN Date 05 04 2013 Page 4 16 Responsable Albert ALARCON Cl U4 53 02 R vision 10811 2 Syntaxe dyn dyna trans DYNA LINE TRAN reuse dyn ODELE mo modele CHAM MATER chmat cham mater 0 CARA ELEM carac cara elem ATR MASS m matr asse DEPL R ATR RIGI k matr asse DEPL R ATR AMOR c matr asse DEPL R 9 ODE_STAT modestat mode meca SCHEMA TEMPS F Y SCHEMA NEWMARK DEFAUT WILSON DIFF CENTRE ADAPT ORDRE2 0 BETA 0 25 DEFAUT beta R 0 GAMMA 0 5 DEFAUT gamma R 0 THETA 1 4 DEFAUT th R 9 ETAT_INIT LE RESULTAT dy dyna trans NUME ORDRE nuord I INST_INIT to R 0 CRITERE RELATIF DEFAUT PRECISION 1 E 06 DEFAUT prec R CRITERE ABSOLU PRECISION prec R DEPL depl ch
11. ECA U4 44 01 ils peuvent tre fournis directement sous forme de vecteurs assembl s ou sous forme de charges qui seront assembl es dans l algorithme La solution X X X est calcul e sur une discr tisation temporelle 1 de l intervalle d tude pr cis par l utilisateur 4 Op randes 4 1 Op rande MODELE 9 MODELE mo Nom du mod le dont les l ments font l objet du calcul dynamique Cet op rande est obligatoire lorsque l on applique une excitation de type charge avec le mot cl EXCIT cf 84 7 4 2 Op rande CHAM MATER 0 CHAM MATER chmat Nom du champ de mat riau affect sur le mod le mo n cessaire lorsque l on applique une excitation de type charge avec le mot cl EXCIT 4 3 Op rande CARA ELEM 9 CARA ELEM carac Nom des caract ristiques des l ments de poutre coque etc n cessaire une excitation de type charge avec le mot cl EXCIT 4 4 Matrices du probl me 4 MATR MASS m lorsque l on applique Concept matrice assembl e de type matr _ asse DEPL R correspondant la matrice de masse du syst me MATR RIGI k Concept matrice assembl e de type matr asse DEPL R correspondant la matrice de rigidit du syst me 0 MATR AMOR c Concept matrice assembl e de type matr asse DEPL R correspondant la matrice d amortissement du syst me Manuel d utilisation Fascicule u4 53 Analyse dynamique Co
12. R PAS PAS MINI Ces op randes ne concernent que le sch ma pas de temps adaptatif 0 VITE MIN NORM DEFAUT MAXI 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 la fr quence apparente 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 pour chaque degr de libert i Ia al Li ley z fap E j At 2m f vnt Vn peut tre calcul e de deux fa ons diff rentes selon la valeur de VITE_MIN 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 cc a default Titre Op rateur DYNA_LINE_TRAN Date 05 04 2013 Page 14 16 Responsable Albert ALARCON Cl U4 53 02 R vision 10811 el Mas kh X y 1 2 Kris 100 libert de m me nature que le degr de libert les plus proches de i dans la num rotation DX ou DY ou DZ NORM y ft Max 10 ms ot k et sont les degr s de min bte 100 Peut tre utilis si l ordre de grandeur de la vitesse ne varie pas trop au cours du temps MAXI vi t Max min 0 lt t lt t 10 ms pour le degr de libert COEF MULT PAS cmp Coefficient de d raffinement du pas de temps
13. XI NORM DEFAUT O COEF MULT PAS cmp R Lt DEFAUT 0 COEF DIVI PAS cdp R 1 3334 DEFAUT 0 PAS LIMI RELA plr R 1 D 06 DEFAUT 0 NB POIN PERIODE npp I 50 DEFAUT 0 NMAX ITER PAS nip 16 DEFAUT 9 PAS MINI dtmin R 4 ENERGIE F 0 ARCHIVAGE F LIST INST list listr8 INST in R PAS ARCH ipa 1 0 CRITERE RELATIF DEFAUT Y PRECISION 1 E 06 DEFAUT prec R CRITERE ABSOLU PRECISION prec R 0 CHAM EXCLU DEPL VITE ACCE LA 4 0 TITRE titre 1 Kn 0 INFO 1 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 fdl html Code Aster Version default Titre Op rateur DYNA LINE TRAN Date 05 04 2013 Page 6 16 Responsable Albert ALARCON Cl U4 53 02 R vision 10811 3 Equations du comportement sous excitation transitoire L op rateur r alise l int gration temporelle directe d un probl me m canique lin aire transitoire de la forme Mx Cx Kx a 1 F x o les matrices M C K sont les matrices r elles assembl es du probl me l ments finis respectivement de masse d amortissement et de rigidit du syst me Les a sont des fonctions du temps cf DEFI FONCTION U4 31 02 et les F sont des vecteurs assembl s issus de chargements en force impos e cf AFFE CHAR M
14. _ORDRE INST_INIT 0 NUME_ORDRE nuord 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 INST_INIT to Instant du calcul pr c dent extraire et prendre comme tat initial dans le cas d une reprise En Pabsence de NUME_ORDRE et INST_INIT l instant de reprise est pris gal au dernier instant de calcul pr c dent archiv 4 6 4 Op rande CRITERE 0 CRITERE 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 Po Titre Op rateur DYNA LINE TRAN Date 05 04 2013 Page 9 16 Responsable Albert ALARCON Cl U4 53 02 R vision 10811 Indique avec quelle pr cision la recherche de l instant 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 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 LINE TRAN Date 05 04 2013 Page 10 16 Responsable Albert ALARCON Cl U4 53 02 R vision 10811 4 6 5 Op rande PRECISION 0
15. am no sdaster VITE vite cham no sdaster ACCE acce cham no sdaster EXCIT F 4 CHARGE chi char meca O FONC MULT fi fonction VECT ASSE va cham no sdaster 0 COEF MULT ai R FONC MULT fi fonction DEPL depl fonction VITE vite fonction ACCE acce fonction MULT APPUI OUT NON DEFAUT DIRECTION dl d2 d3 1 El 9 NOEUD lno 1 noeud 9 GROUP NO lgrno 1 gr noeud 9 EXCIT_RESU TEN Y RESULTAT resuforc dyna_trans COEF MULT ai R Manuel d utilisation Fa Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Scicule u4 53 Analyse dynamique Code Aster Version default Titre Op rateur DYNA_LINE_TRAN Date 05 04 2013 Page 5 16 Responsable Albert ALARCON Cl U4 53 02 R vision 10811 0 AMOR MODAL F AMOR _REDUIT 1 amor 1_R MODE MECA mode mode meca 0 NB MODE nbmode I 9999 DEFAUT 4 0 SOLVEUR F voir le document U4 50 01 et le paragraphe correspondant dans la pr sente doc 4 INCREMENT F LIST_INST litps listr8 PAS dt R 0 INST_INIT ti R 0 INST_FIN tf R NUME FIN nufin I 9 PAS CALCUL Ts DEFAUT ipas I VITE MIN MA
16. d autres valeurs WILSON Sch ma d int gration implicite de type WILSON Avec ce sch ma on peut renseigner THETA th Valeur du param tre Y pour la m thode de WILSON Par d faut 9 1 4 Ce sch ma ne doit pas tre utilis lorsque l on impose des d placements non nuls par l interm diaire d un vecteur assembl Voir R5 05 02 DIFF CENTRE Sch ma d int gration explicite par diff rences centr es L utilisation de ce sch ma impose certaines restrictions d utilisation num r es au 86 3 La description th orique du sch ma est faite dans bib 2 ADAPT ORDRE2 Sch ma d int gration explicite pas de temps adaptatif variante du sch ma des diff rences centr es L utilisation de ce sch ma impose certaines restrictions d utilisation num r es au 6 3 voir bib 2 Nota bene On ne peut pas utiliser les sch mas explicites DIFF CENTRE ADAPT ORDRE avec les l ments de plaque et coque sauf SHB 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 LINE TRAN Responsable Albert ALARCON 4 6 Mot cl ETAT _INIT Date 05 04 2013 Page 8 16 Cl U4 53 02 R vision 10811 Cette fonctionnalit permet une poursuite d un calcul transitoire en prenant comme tat initial un r sultat obtenu par un calcul pr c dent avec DYN
17. ision 10811 5 Concept produit dyn est un concept produit de type dyna_trans qui contiendra de 1 3 champs aux n uds pour chaque instant archiv Ces cham_no ont pour nom symbolique DEP d placement VITI vitesse ACCI my F an my F an acc l ration 6 Phase d ex cution L utilisation des sch mas des diff rences centr es et adaptatifs impose certaines restrictions d utilisation e ces deux sch mas n cessitent l utilisation d une matrice de masse diagonale Un test v rifie que la matrice de masse a t cr e avec l option MASS MECA DIAG de CALC MATR ELEM D autre part la matrice de masse doit tre stock e en ligne de ciel il ne doit pas y avoir d autres conditions aux limites que des degr s de libert bloqu s Un test v rifie qu il n y a pas de conditions aux limites de type liaisons entre degr s de libert Il n est pas non plus possible d imposer des d placements non nuls par l interm diaire d un vecteur assembl pour le sch ma des diff rences centr es on s assure que le pas de temps choisi v rifie les conditions de stabilit dt lt 0 05 f ma avec fyx MAX 1 lt i lt nddl matrices de raideur et de masse Tk E et k et m termes diagonaux des Mi 7 Bibliographie 1 2 BATHE K J Finite Element Procedures in engineering Analysis Prentice Hall 1982 LEGER A C Introduction des sch mas explicites l
18. it ration 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 0 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 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 un Titre Op rateur DYNA LINE TRAN Date 05 04 2013 Page 15 16 Responsable Albert ALARCON Cl U4 53 02 R vision 10811 4 13 Mot cl ARCHIVAGE O ARCHIVAGE Mot cl facteur d finissant l archivage En l absence de ce mot cl facteur tous les pas de temps sont archiv s 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 poursuite 4 13 1 Op randes LIST_INST INST 0 LIST INST list Liste de r els d finissant les instants de calcul pour lesquels la solution doit tre archiv e dans le concept r sultat dyna tran 9 INST Instants de calcul pour lesquels la solution doit tre archiv e dans le concept r sultat dyna tran
19. li s par en analyse thermom canique Remarque importante Les conditions aux limites de type d placement impos non nul peuvent tre impos es avec un vecteur assembl ou une charge il faut alors utiliser imp rativement le sch ma de Newmark 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 LINE TRAN Date 05 04 2013 Page 11 16 Responsable Albert ALARCON Cl U4 53 02 R vision 10811 4 7 3 Op randes MULTI_APPUI ACCE VITE DEPL DIRECTION NOEUD 4 8 4 9 4 9 1 4 10 GROUP NO MODE STAT Dans le cas d une excitation multi appuis MULT APPUI OUI les autres op randes ont exactement la m me signification que dans le mot cl facteur EXCIT de l op rateur DYNA TRAN MODAL U4 53 21 Mot cl EXCIT_RESU Mot cl permettant de d finir plusieurs compl ments de chargement sous forme d une volution transitoire de vecteurs assembl s seconds membres 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 CX N B Cette fa on d introduire l amor
20. pyright 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 LINE TRAN Date 05 04 2013 Page 7 16 Responsable Albert ALARCON Cl U4 53 02 R vision 10811 N B Les trois matrices doivent s appuyer sur la m me num rotation et tre construites avec le m me mode de stockage C est vrai aussi d une matrice d amortissement construite comme combinaison lin aire des matrices de rigidit et de masse par la m thode de Rayleigh utiliser la matrice de la matrice de masse compl te pour construire la matrice d amortissement et la matrice de masse diagonale sch mas explicites tels que DIFF CENTRE OU ADAPT pour l int gration en temps peut mener une instabilit num rique 4 5 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 4 5 1 Op rande SCHEMA NEWMARK Sch ma d int gration implicite de type NEWMARK C est le sch ma par d faut pour l analyse transitoire sur base physique On peut pr ciser les param tres d int gration B et y 9 BETA beta Valeur du param tre B pour la m thode de NEWMARK Par d faut B 0 235 GAMMA gamm Valeur du param tre y pour la m thode de NEWMARK Par d faut y 0 5 Voir R5 05 02 pour le choix
21. t diff rences centr es et pas de temps adaptatif dans l op rateur DYNA LINE TRAN du Code Aster Note EDF HP51 97 067 A 1997 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
22. tissement modal dans un probl me calcul sur base physique peut r duire les propri t s de stabilit des sch mas en temps En particulier pour le sch ma d int gration NEWMARK elle peut conduire r duire le pas de temps par rapport au pas de temps sans amortissement pour viter des divergences num riques Op randes MODE MECA I AMOR _ REDUIT I NB_MODE MODE MECA mode AMOR REDUIT 1 amor 0 NB_MODE nbmode Le concept mode de type mode meca entr par lop rande MODE MECA repr sente la base de modes pr calcul e sur laquelle on d compose l amortissement modal Cette base doit imp rativement avoir le m me profil de num rotation que celui du syst me dynamique d fini par les param tres du mot cl SOLVEUR 84 11 Il es possible de tronquer la base modale un nombre de modes d fini par NB MODE d faut on prend tous les modes de la base modale Les amortissements modaux sous forme r duite sont donn s sous forme d une liste de r els dont le nombre de termes est inf rieur ou gal au nombre de modes pris en compte Si le nombre de termes de la liste est strictement inf rieur on tend cette liste avec la valeur de son dernier terme jusqu ce que sa taille atteigne le nombre de modes calcul s Mot cl ENERGIE 9 ENERGIE gt Ed Ce mot cl permet d activer le calcul du bilan d nergie son affichage en cours de calcul et son stockage dans la table de nom
23. ut utiliser dans DYNA LINE TRAN METHODE MULT FRONT RENUM METIS MDA MD METHODE LDLT RENUM SANS on NE PEUT PAS utiliser dans DYNA LINE TRAN METHODE LDLT RENUM RCMK en r alit on peut le faire mais RCMK est ignor 2 si matrices K et M sont num rot es avec METHODE LDLT RENUM RCMK alors on peut utiliser dans DYNA LINE TRAN METHODE MULT FRONT RENUM METIS MDA MD METHODE LDLT RENUM RCMK on NE PEUT PAS utiliser dans DYNA LINE TRAN METHODE LDLT RENUM SANS en r alit on peut le faire mais SANS est ignor 4 12 Mot cl INCREMENT Mot cl facteur d finissant les instants de calcul 4 12 1 Op randes LIST _INST PAS e Pour les sch mas de Newmark et Wilson LIST INST 1 temp Concept liste de r els de type listr8 Liste de r els d finissant les instants de calcul de la solution e Pour les sch mas des diff rences centr es et pas de temps adaptatif PAS dt D signe le pas de temps utilis par l algorithme Ce mot cl est obligatoire pour le sch ma des diff rences centr es et pour le sch ma adaptatif et non disponible pour les sch mas de Newmark et Wilson 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
Download Pdf Manuals
Related Search