Home

Code_Aster

1. Le r sultat de l archivage pour le premier U2 est le suivant num ro d archivage WW 2 3 4 5 instants correspondants 2 4 6 8 10 Quand on reprend le calcul deuxi me STAT NON LINE l instant de reprise 4 0 est inf rieur au dernier instant archiv 10 0 Le code met donc une erreur fatale pour avertir l utilisateur qu il risque de perdre des donn es archiv es pr c demment correspondant aux instants 6 0 8 0 et 10 0 Pour lever cette erreur fatale il suffit de pr ciser DE TR NUME SUIV OUI Le r sultat final de l archivage pour U2 est alors le suivant les num ros d ordre 3 4 et 5 sont recalcul s et archiv s num ro d archivage 1 2 3 4 5 instants correspondants 2 4 6 8 10 3 17 4 Op rande CHAM EXCLU Q CHAM EXCLU Permet de pr ciser les champs qui ne seront pas archiv s except au dernier pas de temps Le nom des champs exclus d pend des op rateurs 3 18 Mot cl AFFICHAGE Q AFFICHAGE F Ce mot clef facteur permet de personnaliser l affichage du tableau de convergence dans STAT NON LINE et DYNA NON LINE Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster defaut Titre Op rateur STAT NON LINE Date 18 04 2011 Page 33 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 Si ce mot clef n est pas rense
2. LISTI DEFI LIST REEL DEBUT 0 INTERVALLE F JUSQU A 10 NOMBRE 10 U STAT NON LINE INCREMENT F LIST INST LIST1 LIST2 DEFI LIST _ REEL DEBUT 20 INTERVALLE F JUSQU A 30 NOMBRE 10 U STAT NON LINE reuse U INCREMENT F LIST INST LIST2 ETAT INIT F EVOL NOLI U INST ETAT INIT 20 INE effectue le calcul des instants 1 10s INE effectue le calcul des instants 21 30s l tat initial correspondant E par d faut INST 10 Cet tat initial correspond pour ETAT INIT 20 C Exemple pour montrer l int r t de INST_ETAT_INIT pratique quand on fait du cyclique LISTI DEFI LIST REEL DEBUT 0 INTERVALLE F JUSQU A 10 NOMBRE 10 U1 STAT NON LINE INCREMENT F LIST INST LIST1 U2 STAT NON LINE INCREMENT F LIST INST LIST1 ETAT_INIT F EVOL NOLI Ul INST_ETAT_INIT 0 INE INE Premier STAT NON Second STAT NON l instant 10s du premier STAT NON LI ce second STAT NON _ LINE l instant 0s INST 1 effectue le calcul des instants 1 10s effectue le calcul des instants 1 10s l tat initial correspondant E par d faut INST 10 Cet tat initial correspond pour ETAT INIT 0 3 6 5 Op rande PRECISION CRITERE O PRECISION prec Cf U4
3. MIN MAX POINT pi I 0 SOUS POINT spi I Une fois valu les composantes ou la formule sur les composantes on peut e Extraire ces Valeurs sur les points et sous points d int gration avec EVAL ELGA VALE Dans ce cas il faut pr ciser explicitement le point et le sous point d int gration par POINT et SOUS POINT Les sous points d int g r ation apparaissent pour les l ments de structures poutres plaques coques tuyaux etc Demander d extraire le maximum EVAL ELGA MAX ou le minimum EVAL ELGA MIN sur tous les points et sous points d une maille Si on demande explicitement un point et un sous point on aura autant de r alisations que de points demand s multipli par le nombre de composantes demand es Par contre si on demande le maximum ou le minimum il y a aura une seule observation par composante demand e EVAL CHAM VALE DEFAUT MIN MAX MOY Une fois valu les composantes ou la formule sur les composantes ai nsi que le point sous point d extraction on peut e Extraire ces valeurs sur toutes les mailles avec EVAL CHAM VALI D Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 35 37 Responsable
4. Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur S TAT_NON LINE Responsable Mickael ABBAS ABSOLU Date 18 04 2011 Page 9 37 Cl U4 51 03 R vision 6167 Si CRITERE RELATIF PRECISION 1 0E 6 DEFAUT prec R Si CRITERE ABSOLU PRECISION prec R SUIVI DDL Fr NOM CHAM DEPL VITE ACCE DEPL ABSOLU VITE ABSOLU ACCE ABSOLU SIEF ELGA VART ELGA FORC NODA EVAL CHAM Zi VALE t DEFAUT MTN MAX MOY NOM CMP lnocmp l Kn EVAL CMP Z VALE t DEFAUT FORMULE S EVAL CMP FORMULE 9 FORMULE form formule aster Si CHAM est de type FLGA SIEF ELGA VARI _ ELGA TOUT OUI Y3 DEFAUT GROUP MA lgrma 1 gr maille Q MAILLE Ima 1 maille EVAL ELGA Zi VALE gt DEFAUT MIN Zil MAX r Si EVAL ELGA VALE POINT o Pi y I SOUS POINT spi I Si CHAM est de type NOEU TOUT OUI DEFAUT GROUP MA lgrma 1 gr maille Q MAILLE Ima 1 maille Q NOFUD no noeud Q GROUP NO PTG gr noeud Q TITRE LL ets s last k INFO DEFAUT rg TITRE tx Kn
5. Manuel d utilisation Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Fascicule u4 51 Analyse statique Code Aster Moar default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 10 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 3 Op randes 3 1 Op randes MODELE CHAM MATER CARA ELEM MODELE mo Nom du mod le dont les l ments font l objet du calcul m canique CHA ATER chmat Nom du champ de mat riau affect sur le mod le mo Attention toutes les mailles principales du mod le doivent tre associ es un mat riau sinon erreur fatale avec message peu explicite Q CARA ELEM carac Nom des caract ristiques carac des l ments de coque poutre tuyau barre c ble et l ments discrets affect s sur le mod le mo videmment ce mot cl est optionnel si le mod le ne contient pas de tels l ments il n est pas utile en revanche si le mod le contient de tels l ments il est obligatoire 3 2 Mot cl EXCIT EXCIT F Ce mot cl facteur permet de d crire chaque occurrence une charge sollicitations et conditions aux limites et ventuellement un coefficient multiplicateur et ou un type de charge 3 2 1 Op randes CHARGE CHARGE ch H ch est le chargement m canique comportant ventuellement l volution d un champ de temp rature pr cis la me occurrence de
6. NORM INCR DEPL DEFAUT ANGL INCR DEPL RESIDU Cet op rande permet de s lectionner la m thode permettant de choisir la valeur de pilotage dans le cas o plusieurs solutions sont fournies par la r solution de pilotage NORM_INCR_DEPL permet de s lectionner la valeur de pilotage par la plus petite norme de l incr ment de d placement sur le pas de temps consid r ANGL INCR DEPL permet de s lectionner la valeur de pilotage par le plus petit angle entre le d placement obtenu pour le pas de temps courant et le d placement obtenu pour le pas de temps pr c dent RESIDU permet de s lectionner la valeur de pilotage conduisant au plus petit r sidu Remarque Si on fait une reprise de calcul reuse avec le mot clef SELECTION ANGL INCR DEPL il est important de garder l esprit que ce crit re n cessite les deux pas de temps pr c dents Il faudra donc bien prendre soin d archiver correctement les r sultats du pr c dent calcul au risque d obtenir des r sultats faux Une alarme avertit l utilisateur 3 13 Mot cl SOLVEUR Q SOLVEUR EU La syntaxe de ce mot cl commun plusieurs commandes est d crite dans le document U4 50 01 3 14 Mot cl CONVERGENCE CONVERGENCE EI Si aucun des deux op randes suivants n est pr sent alors tout se passe comme si RESI GLOB RELA 1 E 6 3 14 1 Op rande RESI GLOB RELA RESI G
7. PRECISION 1 E 6 prec o CRITERE RELATIF ABSOLU DETR NUME SUIV OUT CHAM EXCLU list Ext AFFICHAGE ek UNITE unite INFO RESIDU NON OUI Q OBSERVATION F NOM CHAM EVAL CHAM NOM CMP EVAL CMP Date 18 04 2011 Page 8 37 Cl U4 51 03 DEPL R vision 6167 FAUT EFAUT DEFAUT VEE ACC DEPL VIT ACC SIE SE FOR VA J Y M M M Inocmp Pa Sal VARI EL FORMULI C NODA LE ALE IN AX OY DEFAUT l Kn VALE DEFAUT PJ L Si EVAL CMP FORMULE FORMULE form formule aster Si CHAM est de type ELGA SIEF ELGA VARI ELGA Q TOUT OD DEFAUT Q GROUP MA lgrma l gr maille OD MAILLE lma 1 maille EVAL ELGA ZY VALE DEFAUT MIN u MAX Si EVAL ELGA VALE POINT Dire 3 o SOUS POINT spi I Si CHAM est de type NOEU Q TOUT OD Ey DEFAUT Q GROUP MA lgrma l gr maille Q MAILLE lma 1 maille Q NOEUD no noeud Q GROUP NO rno gr noeud LIST INST linst listr8 O INST linst 1 R 0 PAS OBSE pas 1 D CRITERE RELATIF DEFAUT Manuel d utilisation Fascicule u4 51 Analyse statique
8. par un concept evol noli ETAT INIT avec EVOL NOLI il s agit de l instant dans le pr c dent calcul istetaini instini A Exemple simple comportement par d faut LISTI DEFI LIST REEL DEBUT 0 INTERVALLE F JUSQU A 4 NOMBRE 4 U STAT NON LINE INCREMENT F LIST_INST LIST1 LIST2 DEFI LIST REEL DEBUT 4 INTERVALLE F JUSQU A 10 NOMBRE 6 U STAT NON LINE reuse U INCREMENT F LIST INST LIST2 Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur S TAT_NON LINE Responsable Mickael ABBAS Premier STAT_NON_ Second STAT_NON_ correspondant au temps 4s ETAT INIT Date 18 04 2011 Page 14 37 Cl U4 51 03 R vision 6167 F EVOL NOLI U JINE effectue le calcul pour les instants 1 2 3 et 4s INE effectue le calcul pour les instants 5 6 7 8 9 et 10s l tat initial B Exemple pour montrer l int r t de INST_ETAT_INIT deux listes d instants diff rentes Premier STAT_NON_ Second STAT_NON_ l instant 10s du premier STAT NON LIN ce second STAT NON LINE l instant 20s INST
9. pour pouvoir effectuer une reprise 2 si on emploie un acc s par liste d instants alors les instants de calculs nouvellement cr s par red coupage automatique du pas de temps ne sont pas archiv s Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Ast a ogae_Aster default Titre Op rateur STAT NON UNE Date 18 04 2011 Page 32 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 3 L tat initial est syst matiquement archive sous le num ro d ordre 0 d s lors que l on n est pas en reprise de calcul pas de reuse 3 17 2 Op rande PRECISION CRITERE PRECISION 1 e 6 DEFAUT prec CRITERE RELATIF DEFAUT ABSOLU Cf U4 71 00 3 17 3 Op rande DETR_NUME SUIV DETR NUME SUIV NON DEFAUT Z OD 1 Cette op ration peut conduire craser des num ros d ordre pr existants le mot cl DETR NUME SUIV confirme cette destruction tandis que son absence met fin au calcul Illustration LIST DEFI LIST REEL DEBUT 0 INTERVALLE F JUSQU A 10 NOMBRE 5 U2 STAT NON LINE INCREMENT F LIST INST LIST LA U2 STAT NON LINE reuse U2 ETAT INIT F EVOL NOLI U2 INST 4 INCREMENT F LIST INST LIST ARCHIVAGE F LIST INST LIST DETR NUME SUIV OUI
10. contact R chi OUT lma fmult _F voir le document evol noli modele cham mater cara elem char mecal fonction formule DEFAUT char contact char mecal DEFAUT 1 maille cham elem SIEF R cham elem VARI R cham no DEPL R evol nolil EFAUT ow m H EFAUT H Ke Manuel d utilisation Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Fascicule u4 51 Analyse statique Code Aster Version default Titre Op rateur S TAT_NON LINE Responsable Mickael ABBAS ETHODE EWTON PREDICTION EVOL NOLI MATRICE REAC INCR REAC ITER REAC ITER ELAS PAS MINT ELAS TANGENTE ELASTIQUE EXTRAPOL DEPL CALCULE evol noli TANGENTE ELASTIQUE Jiz mf 0 fit 0 ZE 0 pasmini Ka bal Date 1 Cl U4 51 03 8 04 2011 Page 6 37 R vision 6167 EFAUT DEFAUT evol noli DEFAUT Ss als es Wels D D D C gd O IMPL EX Pi PREDICTION TANGENTE DEFAUT REAC INCR Ty DEFAU
11. es EVOL NOLI La structure de donn es EVOL NOLI contient la liste des champs archiv s au cours du calcul selon les diff rentes options de mot clef ARCHIVAGE Par d faut elle contient pour chaque instant la liste des champs suivants DEPL champ aux n uds des d placements SIEF ELGA champ aux points de Gauss des contraintes VARI ELGA champ aux points de Gauss des variables internes e COMPOR carte du comportement Selon certaines options de calcul d autres champs seront pr sents VALE CONT champ aux n uds des informations sur le contact frottement voir U4 44 11 pour plus de d tails sur le contenu de ce champ INDC ELGA champ aux points de Gauss des statuts de contact pour le cas XFEM avec contact COHE ELGA champ aux points de Gauss du param tre de coh sion pour le cas XFEM avec RELATION CZM SECO ELGA champ aux points de Gauss des statuts de frottement pour le cas XFEM avec contact et frottement En plus de ces champs la structure de donn es contient galement des param tres chaque instant on stocke au minimum Nom Mot clef origine Description Type INST Valeur de l instant de calcul R EXCIT EXCIT Informations sur les chargements K24 MODELE ODELE Mod le K8 CARAELEM CARA ELE Cara
12. D PLAN THM et la force centrifuge en grands d placements mot cl ROTATION dans AFFE CHAR MECA Si tchi vaut DIDI alors les conditions de Dirichlet d placements impos s conditions lin aires s appliqueront sur l incr ment de d placement partir de l instant donn sous ETAT _INIT NUME DIDI par d faut l instant de reprise du calcul et non sur le d placement total Par exemple pour un d placement impos mot cl DDL IMPO de AFFE CHAR MECA la condition sera de la forme u u d o u est le d placement d fini par NUME DIDI et non u d 3 3 Mot cl CONTACT CONTACT contact Ce mot cl simple permet d activer la r solution du contact frottement ou la prise en compte d une liaison unilat rale contact est un concept issu de l op rateur DEFI CONTACT U4 44 11 Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 12 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 Attention Ce mot cl simple n accepte qu un seul concept On ne peut donc pas m langer dans un m me calcul non lin aire la r solution du contact et la prise en compte d une liaison unilat rale On ne peut pas non plus m langer les diff rentes formulations DISCRETE CONTINUE et XFEM 3 4 Mot cl SOUS STRUC Q SOUS STRU
13. EXCIT Remarques 1 Dans un calcul thermo m canique si la temp rature initiale est diff rente de la temp rature de r f rence donn e dans l op rateur AFFE MATERIAU le champ de d formation associ l instant initial peut tre incompatible et donc conduire un tat de contraintes et de variables internes associ non nul Si l on utilise une relation de comportement incr mentale mot cl facteur COMP_INCR et si on ne d finit pas explicitement un tat de contraintes et de variables internes initial associ un champ de temp rature initiale diff rente de la temp rature de r f rence le champ de contraintes et de variables internes calcul au premier incr ment ne tiendra compte que de la seule variation de temp rature entre l instant initial et le premier instant et non des ventuelles contraintes de compatibilit associ es la temp rature initiale Pour prendre cet tat initial en compte il faut le donner explicitement par exemple gr ce aux mots cl s SIGM DEPL VARI et VARI NON LOCAL dans ETAT INIT Pour viter de telles situations qui peuvent conduire des erreurs de calculs il vaut mieux commencer un calcul en consid rant qu il faut partir d un tat vierge 2 Si on r alise un calcul en axisym trique et que l on impose des forces nodales ces efforts doivent tre divis s par 2Tr on travaille sur un secteur d un radian par rapport aux chargements r els De m me si l on souhaite calculer la r
14. MAX le minimum EVAL CHAM EVAL CHAM MOY MIN ou la moyenne Exemple Extrai r e e maximum de la composante DX du d placement sur GROUP NO TOT O OBSERVATION F NOM CHAM DEPL GROUP_NO TOTO EVAL CHAM MAX NOM CMP DX 3 20 Mot cl SUIVI _DDL SUIVI DDL _F Ce mot cl permet de post traiter certains champs aux n uds ou aux l ments sur des parties de mod le toutes les it rations de Newton et les afficher dans le tableau de convergence Le nombre simultan de SUIVI_DDL d pend des colonnes affich es et donc des fonctionnalti s activ es Le SUIVI DDL a la m me syntaxe que OBSERVATION pour l extraction des champs sauf que l on ne donne pas d informations sur les in s tants extraire pui s qu on le r alise chaque it ration de Newton il n y a pas les mots clefs LIST_INST INST PAS_OBSE CRITERE PRECISION Q TITRI D ltitre list k Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Maer default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 36 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 Ce mot clef attend une li s te de trois chaines au maximum et permet de nommer la colonne du tableau d affichage Les chaines sont tronqu s 16 caract res 3 21 Contenu de la structure de donn
15. NOM CHAM ainsi que ses composantes donn es par leur nom NOM CMP On ne peut d finir que 20 composantes maximum par occurrence du mot clef facteur OBSERVATION 3 19 4 Op randes TOUT NOEUD GROUP NOEUD MAILLE GROUP MA TOUT NON DEFAUT OUT NOEUD no no GROUP NO grno grno MAILLE lma ma GROUP MA lgrma grma Ces op randes permettent de d finir le support g om trique de post traitement pour des champs aux n uds DEPL VITE ACCE DEPL ABSOLU VITI ACCE ABSOLU VALE CONT FORC NODA on extrait la liste des noeuds e pour des champs aux points de Gauss SIEF ELGA VARI ELGA on extrait la liste des mailles Attention ne pas utiliser TOUT OUI sur de s gros maillages Kal ABSOLU 3 19 5 Observation d un champ ELGA O EVAL CMP VALE DEFAUT FORMULE Q FORMULE form formule aster On commence par choisir les composantes ou la formule entre les composantes Si EVAL CMP VALE on extrait simplement la liste des composantes donn e par NOM CMP eSI EVAL CMP FORMULE on value la formule donn e par le mot clef simple FORMULE Si on applique une formule sur les composantes on aura donc une valeur et donc une observation sinon on aura autant d observations que de composantes dans la liste NOM CMP EVAL ELGA VALE DEFAUT
16. NUME_INST_INIT le d faut est calcul de la mani re suivante 1 Si un tat initial est pr cis op rande ETAT_INIT et s il d finit un instant correspondant par EVOL_NOLI ou INST_ETAT_INIT alors l instant initial est celui d fini par cet tat initial 2 S il n y a pas d tat initial op rande ETAT_INIT absent ou qu il ne d finit pas d instant correspondant les champs sont donn s dans ETAT_INIT sans pr ciser INST ETAT INIT alors on prend le premier instant de la liste d instants HOME _INST_INIT 0 3 En cas d archivage voir mot clef ARCHIVAGE l instant initial en poursuite est le dernier pas archiv et non celui d fini dans INST_INIT NUME INST FIN nufin INST FIN instfin Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster defauit Titre Op rateur STAT_NON_LINE Date 18 04 2011 Page 16 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 L instant final dernier pas calcul est d sign de la m me mani re que l instant initial soit NUMI E INST_ FIN soit INST FIN sauf qu il n est pas possible de faire r f rence l instant de l tat initial Attention Si le re d coupage automatique du pas de temps est activ NUME INST FIN n en tient pas compte et travaille toujours sur la liste d instants initial NUME INST INIT et NUME INST FINne sont actifs qu l init
17. Page 20 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 3 9 3 3 10 3 10 1 3 11 3 11 1 On peut d finir une fr quence de r actualisation de la matrice de d charge avec le mot cl REAC ITER ELAS analogue de REAC ITER Si la matrice de d charge ne d pend pas de l tat de d formation ce qui est le cas pour les mat riaux plastiques mais pas pour les mod les d endommagement prendre REAC ITER ELAS 0 puisqu elle sera la m me au cours des it rations Utilit Cette option peut tre utile lorsque le red coupage automatique du pas de temps ne suffit pas faire converger un calcul Par exemple dans le cas de lois adoucissantes la matrice tangente peut devenir singuli re et il vaut donc mieux utiliser la matrice lastique pour converger O MATRICE ELASTIQUE La matrice utilis e correspond au calcul lastique elle n est valu e qu une fois l instant initial en d but d algorithme Cette matrice lastique est calcul e en utilisant le module d Young donn sous le mot cl ELAS de l op rateur DEFI MATERIAU et non pas la pente l origine de la courbe de traction donn e sous le mot cl TRACTION et qui sert elle dans l expression des relations de comportement VMIS ISOT TRAC ECMI ISOT TRAC VISC ISOT TRAC U4 51 11 Op rande EVOL NOLI O EVOL NOLI evol noli Nom du concept de type evol noli qui ser
18. de Newton de converger plus rapidement Lorsque le mot cl est absent c est la matrice tangente en vitesse option RIGI _ MECA TANG qui est utilis e si l on a choisi pour la m thode de Newton une MATRICE TANGENTE et c est la matrice lastique option RIGI MECA qui est utilis e si on a choisi MATRICE ELASTIQUE TANGENTE On utilise la matrice tangente du probl me en vitesse option RIGI MECH TANG ELASTIQUE On utilise la matrice lastique option RIGI MECA EXTRAPOL On calcule l estimation de l incr ment de d placement partir de l incr ment total obtenu comme solution au pas de temps pr c dent pond r par le rapport des pas de temps On projette cette estimation sur l ensemble des champs cin matiquement admissibles i e satisfaisant les conditions aux limites de Dirichlet selon la norme donn e par la matrice lastique qui doit donc tre calcul e Cette fonctionnalit est int ressante dans le cas de l utilisation de sch mas d int gration locale explicite de type Runge Kutta qui ne fournissent pas de matrice tangente dans ce cas la m thode de Newton utilise une matrice lastique mais le nombre d it rations n cessaires peut tre lev L utilisation de l extrapolation peut am liorer les performances DEPL_CALCULE Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster a Titre Op rateur STAT N
19. ind pendant de la g om trie mais sera pilot gr ce au mot cl PILOTAGE 21 Les charges pilotables doivent tre issues d AFFE CHAR MECA ou d AFFE CHAR MECA F si ce n est pas une fonction d pendant du temps et ne pas tre affect es du mot cl FONC MULT On ne peut pas piloter les chargements de pesanteur la force centrifuge les forces de Laplace les chargements thermiques ou de d formations initiales ou an lastiques et les conditions de liaison Si tch vaut SUIV le chargement est dit lt suiveur gt c est dire qu il d pend de la valeur des inconnues par exemple la pression tant un chargement s appliquant dans la direction normale une structure d pend de la g om trie actualis e de celle ci et donc des d placements Un chargement suiveur est r valu chaque it ration de l algorithme de r solution Un chargement fixe n est r valu qu chaque nouvel instant et seulement si ch d pend du temps d fini dans AFFE CHAR MECA F et param tr par l instant Actuellement les chargements qui peuvent tre qualifi s de SUIV sont le chargement de pesanteur pour l l ment de CABLE POULIE la pression pour les mod lisations 3D 3D SI D_PLAN D PLAN SI AXIS AXIS SI C PLAN C PLAN SI et pour toutes les mod lisations THM 3D HHM 3D HM 3D JOINT CT 3D THH 3D THHM 3D THM AXIS HHM AXIS HM AXIS THH AXIS THHM AXIS THM D PLAN HHM D PLAN HM D PLAN THH D PLAN THHM
20. qui sont piloter en d placement Ces param tres de pilotage permettent de stopper les calculs au m me chargement lorsque la structure est suffisamment adoucie Cette strat gie rend la comparaison plus ais e gr ce au contr le du dernier point de pilotage Avec PROJ BORNE NON on arrive dans certains cas d bloquer les calculs qui autrement ne convergent pas avec les conditions trop restrictives impos es via etarmin etarmax Soit on Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster W Titre Op rateur STAT NON LINE Date 18 04 2011 Page 27 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 pilote une pression impos e la structure et on s attend garder cette pression positive En fixant etarmin 0 le calcul s arr te en chec de pilotage En revanche en imposant etarmin l g rement n gatif on autorise de facto le passage par un tat non physique pendant les it rations de Newton ce qui facilite la convergence L tat converg dans ce cas peut aussi bien tre physique pression positive ou non physique C est la valeur de etarmin 0 qui gouvernera le comportement en cas de convergence hors borne Cette strat gie permet de ne conserver que les valeurs de pilotage positives si on trouve au moins une valeur de pilotage positive 3 12 110p rande SELECTION SELECTION
21. structure de donn es EVOL NOLLI ua ssssasssssasssssswasxawas 34 S 22 Opra nde IN EO EE 35 3 29 Opefande AE Ce 35 Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur STAT NON LINE Responsable Mickael ABBAS 2 Syntaxe Date 18 04 2011 Page 5 37 Cl U4 51 03 R vision 6167 statnllevol noli reuse MODELE CHA ATER Q CARA ELE EXCIT CHARGI FONC MULT al TYPE CHARG 9 CONTACT 9 SOUS STRUC CAS CHARGI TOUT Ia e SUPER MAILLE FONC MULT COMP_INCR COMP_ELAS _F voir le document o ETAT INIT FA ISIGM sig VARI vain DEPL depl EVOL_NOLI evol NUME ORDRE nuini INST instini PRECISTON 1 0E 3 prec CRITERE RELATIF ABSOLU NUME DIDI nudidi INST_ETAT_INIT istetaini INCREMENT Pi LIST INST litpsr8 litps NUME_INST_INIT nuini INST_INIT instini NUME_INST_ FIN uti INST FIN instfin PRECISION 1 0E 3 prec ERRE TEMPS NON OT STAT NON LINE statnl mo chmat carac ch fi FIXE CSTE FIXE PILO SUIV DIDI rJ E
22. sultante des efforts le r sultat est multiplier par jr pour avoir la r sultante totale sur la structure compl te De m me en contraintes planes ou en d formation plane on travaille sur une paisseur unit les efforts sur l paisseur appliqu s doivent tre divis s par l paisseur les efforts r els sont obtenus en multipliant par l paisseur les efforts du calcul Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster erau Titre Op rateur S TAT_NON LINE Date 18 04 2011 Page 11 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 3 2 2 Op rande FONC MULT 0 FONC MULT f Ke est la fonction du temps multiplicatrice du chargement pr cis la i e occurrence de EXCIT Le chargement et les conditions aux limites pour n occurrences du mot cl facteur EXCIT sont Jes fch i 1 Pour les conditions de Dirichlet bien entendu seule la valeur impos e est multipli e par f Par d faut f l i 3 2 3 Op rande TYPE CHARGE Q TYPE CHARGE FIXE CSTE DEFAUT SUIV EE FIXE PILO Par d faut tchi vaut FIXE CSTE cela correspond un chargement appliqu sur la g om trie initiale et non pilot II peut cependant tre une fonction et en particulier d pendre du temps Si tch vaut FIXE PILO le chargement est toujours fixe
23. tail DEFORMATION DEFORMATION garantit qu au moins un point de Gauss de la structure voit sa d formation voluer de fa on monotone On v rifie la relation cmult max Asz At point de Gauss Ce mode de pilotage est valable pour toutes les lois de comportement y compris en grandes d formations SIMO MIEHE PRED ELAS PRED ELAS assure qu au moins un point de Gauss de la structure sorte du seuil d lasticit lin aris n At DES l ni ein d une quantit Geh On v rifie la relation cmult k max el A t point de Gauss Ce mode de pilotage est valable uniquement pour les lois ENDO FRAGILE avec la version locale et la version non locale GRAD EPSI ENDO SCALAIRE avec la version non locale ENDO ISOT BETON et ENDO ORTH BETON avec la version locale et la version non locale BARENBLATT BETON DOUBLE DP CZM EXP avec les l ments discontinuit interne ELDI CZM OUV MIX et CZM TAC MIX l ments d interface INTERFACE CZM EXP REG l ments de joint JOINT ou mod lisation X FEM et CZM LIN REG l ments de joint La fixation du param tre cmult est difficile d finir du premier coup parce que la notion de sortie de crit re n est pas intuitive et varie selon les lois de comportement Pour les lois cmult ENDO FRAGILE ENDO SCALAIRE et ENDO ISOT BETON un
24. 67 tx est le titre du calcul Il sera imprim en t te des r sultats Voir U4 03 011 Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html
25. 71 00 pour la syntaxe d taill e Ce param tre sert rep rer le bon num ro d ordre NUME ORDRE quand l utilisateur renseigne l instant INST En effet les instants dans STAT NON LINE sont rep r s par un num ro d ordre un entier Si l utilisateur veut utiliser un instant un r el et non un num ro d ordre pour INST l op rande pr cision permet de s lectionner ce num ro d ordre Exemple NUME 1 2 3 4 5 6 7 INST 0 0010 0 0020 0 0030 0 0040 0 0050 0 0060 0 0070 Si l utilisateur veut s lectionner l instant correspondant NUME 4 il lui suffit de dire INST 0 004 Par contre pour le deuxi me exemple Manuel d utilisation Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Fascicule u4 51 Analyse statique Code Aster default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 15 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 NUME 1 2 3 4 5 6 7 INST 0 10000001 10 10000002 0 10000003 0 10000004 10 10000005 10 10000006 0 10000007 Si l utilisateur veut s lectionner l instant correspondant NUME 4 il ne lui suffit pas de dire 0 10000005 0 10000004 1E 7 0 10000004 qui est sup rieur la valeur de pr cision par d faut 1E 6 On ne pourra donc pas distinguer NUME 3 4 et 5 le code s arr te alors en erreur fatale Il suffit alors de changer le param tre PRECISION pour p
26. A iii 28 3 147 Op rande E 28 3 15 Motcl CRIT TEE 28 3 193 Op rand CHAR E WEE 29 3 14 Op rand NB EE 29 315 5 Op rande Oe RE EE 29 3 156 Op rande DDL E eege 30 3 16 Motile SENSIBILE EE 30 S L Motole ARCHIVAGE siteadi antaa aeaiiai aaa a ANa 30 3 17 1 Op rande LIST_INST INSTI PAS ARCH geesde 30 3 17 2 Op rand Tele E WR GE 30 EREECHEN 30 317A Opetrande ATC RE EN Kee LEE EN 3 18 2 Op rande INFO RESIDU uu uuu uuu uuu uuu alssasssesaqasiukusasckasssakaqushassQakanashianuqasaqasssu 31 E Motele OB SERVAWON E 32 3 19 1 Op erandes LIST_INST INST PAS OBSE Augen 32 3 19 2 Op randes PRECISION CRITERE u u uu u uuu uuu aaia aaa aaia aaie 32 3 19 3 Op randes NOM CHAM NOM GNP u uuu aa Dan ee ee Nee dee 32 3 19 4 Op randes TOUT NOEUD GROUP_NOEUD MAILLE GROUP_MA 32 3 19 5 Observation d un champ ELGA uuu uu uuu uuu a aaa aqa 33 3 19 6 ODsemation dun chiamp NOE uuu uuu nm aus eaan ariaa dus bias shila kdo 33 Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster default Titre Op rateur STAT_NON_LINE Date 18 04 2011 Page 4 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 S 20 Motele SUMO ge 34 3 21 Contenu de la
27. C Ce mot cl facteur permet de pr ciser quels sont les chargements utiliser pour les sous structures statiques qui font alors obligatoirement partie du mod le En son absence les chargements sur les sous structures sont nuls Ces chargements s ajoutent aux chargements l ments finis qui peuvent tre appliqu s sur le reste du mod le Pour plus de pr cision concernant l utilisation de sous structures lastiques lin aires dans une structure non lin aire on se reportera la documentation U2 07 021 et le cas test ssnv193a 3 4 1 Op rande CAS CHARGE CAS CHARGE nocas nocas est le nom du cas de charge utiliser Voir op rateur MACR ELEM STAT U4 62 01 3 4 2 Op randes TOUT SUPER MAILLE TOUT OUI Ce mot cl permet d affecter le chargement nocas toutes les sous structures du mod le SUPER MAILLE l mail Ce mot cl facteur permet de n affecter le chargement nocas qu certaines sous structures 3 4 3 Op rande FONC MULT Q FONC MULT f est la fonction du temps multiplicatrice du chargement pr cis la i me occurrence de SOUS STRUCT Le comportement de ce mot cl est le m me que pour son occurrence dans EXCIT 3 5 Mots cl s COMP_INCR et COMP_ELAS La syntaxe de ces mots cl s communs plusieurs commandes est d crite dans le document U4 51 11 3 6 Mot cl ETAT _INIT 9 ETAT_INIT Ce mot clef permet de d finir un tat initial de r f r
28. CONVERGENCE F RESI GLOB RELA 1 E 6 RESI GLOB MAXI resmax RESI GLOB RELA resrel RESI COMP RELA rescmp RESI_REFE_RELA resref SIGM _ REFE sigref EPST REFE sigref DEPL REFE depref FORC BEER forref VARI_REFE varref FLUX THER REFE sigref FLUX HYD1 REFE sigref FLUX HYD2 REFE sigref ITER GLOB ELAS 25 maxelas ITER GLOB MAXI 10 maglob TYPE PIC PLATEAU PLATEAU ITER 3 plaite PLATEAU RELA 1 E 3 plarel ARRET es OUL NON o CRIT FLAMB Pi NB _FREQ 135 nbfreq CHAR CRIT 10 10 Ke et eu LIST_INST list r8 INST see PAS CALC npas Q PRECISION 1 E 6 prec Q CRITERE RELATIF ABSOLU Q SENSIBILITI Kal F voir le document Date 18 04 2011 Page 7 37 Cl U4 51 03 R vision 6167 sd fiss xfem e ti rz D S l U4 50 01 EFAUT CHOHO J Z Z TT TT TT TU GZ alle DEFAUT I DEFAUT listr8 R I DEFAUT R DEFAUT U4 50 021 Manuel d utilisation Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Fascicule u4 51 Analyse statique Code Aster Version default Titre Op rateur STAT NON LINE Responsable Mickael ABBAS ARCHIVAGE E LIST_INST list r8 INST 1 r8 PAS_ARCH npas
29. Code Ast We oae_Aster default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 1 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 Op rateur STAT_ NON LINE 1 But Calculer l volution m canique ou thermo hydro m canique coupl e en quasi statique d une structure en non lin aire La non lin arit est li e soit au comportement du mat riau par exemple plastique soit la g om trie par exemple en grands d placements soit au contact frottement Pour avoir des d tails sur la m thode de r solution employ e on se reportera la documentation de r f rence R5 03 01 L volution peut tre tudi e en plusieurs travaux successifs concept r entrant soit en poursuite le dernier instant calcul est l instant initial du calcul suivant soit en reprise en partant d un instant ant rieur Si le temps n cessaire pour effectuer le calcul n est pas suffisant le programme s interrompt mais les r sultats d j calcul s sont sauvegard s si une base de donn es a t d finie dans le profil d tude de l utilisateur Produit une structure de donn es de type evol noli Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Katie default Titre Op rateur STAT_NON_LINE Date 18 04 2011 Page 2 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 Table des mati
30. ITER SIMULT U4 52 03 Ceci peut tre assez co teux pour les probl mes de grande taille Pour les grands d placements et les grandes d formations on r sout der KI A I 0 car KI contient alors K Le crit re est alors un crit re d instabilit quand change de signe donc passe par 0 le chargement est instable On stocke le mode propre correspondant la plus petite charge critique en valeur absolue dans l objet r sultat sous le nom MODE FLAMB Ce mode propre peut tre extrait et visualis comme un champ de d placements ou un mode propre classique ll est normalis 1 sur la plus grande composante de d placement L analyse de stabilit lin aire ne permettant pas de tenir compte de l aspect suiveur de certaines forces il faut alors utiliser CRIT FLAMB La documentation U2 08 04 pr sentent les diff rentes approches pour les analyses de stabilit dans Code_Aster donc CRIT_ FLAMB Op rande LIST_INST I INST PAS CALC LIST INST list r8 INST 1 r8 PAS CALC npas Les instants pour lesquels on veut faire un calcul de stabilit sont donn s par une liste d instants 1ist_r8 ou 1_r8 oupar une fr quence PAS CALC tous les npas de temps En l absence de ces mots cl s le crit re est calcul tous les pas de temps 3 15 2 Op rande PRECISION CRITERE PRECISION 1 e 6 DEFAUT prec CRITERE RELATIE DEFAUT ABSOLU Permet
31. LOB MAXI 9 IRESI GLOB RELA resrel R L algorithme continue les it rations globales tant que max F gt resrel max L id nbddl o F est le r sidu de l it ration n et L le vecteur du chargement impos et des r actions d appuis Cf R5 03 01 pour plus de d tails Lorsque le chargement et les r actions d appui deviennent nuls c est dire lorsque L est nul par exemple dans le cas d une d charge totale on passe du crit re de convergence relatif au crit re de convergence absolu RESI GLOB MAXI Cette op ration est transparente pour l utilisateur message d alarme mis dans le fichier messi Lorsque le vecteur L redevient diff rent de z ro on repasse automatiquement au crit re de convergence relatif REST GLOB RELA Si cet op rande est absent le test est effectu avec la valeur par d faut sauf si RESI GLOB MAXI est pr sent Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 28 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 IREST GLOB MAXI resmax R L algorithme continue les it rations globales tant que max F gt resmax i 1 nbddl o F est le r sidu de l it ration n Cf R5 03 01 pour plus de d tails Si cet op rande est absent le test n est pas effectu Si RESI GLOB RELA et
32. Mickael ABBAS Cl U4 51 03 R vision 6167 le minimum es Demander d extraire le maximum EVAL CHAM MOY EVAL CHAM MAX EVAL CHAM MIN ou la moyenne Exemple Extrai r e e maximum de la trace du tenseur des contraintes sur le GROUP MA TOT O trace FORMULE VALE 0 333 X SIXX SIYY SIZZ NOM PARA SIXX SIYY SIZZ OBSERVATION F NOM CHAM SIE F ELGA GROUP MA TOTO EVAL CHAM MAX NOM CMP SIXX SIYY SIZZ EVAL CMP FORMULE FORMULE trace EVAL ELGA MAX 3 19 6 Observation d un champ NOEU Q EVAL CMP VALE DEFAUT FORMULE FORMULE form formule_aster On commence par choisir les composantes ou la formule entre les composantes eSI EVAL_CMP VALE on extrait simplement la liste des composantes donn e par NO eSI EVAL CMP FORMULE on value la formule donn e par le mot clef simple FORMULI Si on applique une formule sur les composantes on aura donc une valeur et donc une observation sinon on aura autant d observations que de composantes dans la liste NOM_CMP E D D EVAL CHAM VALLE DEFAUT Zi MIN w Zi MAX O Jr MOY y e Une fois valu les composantes ou la formule sur les composantes on peut e Extraire ces valeurs sur toutes les mailles avec EVAL CHAM VALE Demander d extraire le maximum EVAL CHAM
33. OMBRE INCREMENT F LIST INST LIST INST_FIN 4 reuse U INCREMENT F LIST INST LIST INST_INIT 8 ETAT_INIT F R vision 6167 F JUSQU A 10 NOMBRE 10 U EVOL NOLI U 3 7 3 Op rande PRECISION Q PRECISION prec Cf U4 71 00 pour la syntaxe d taill e Ce param tre sert rep rer le bon num ro d ordre NUME INST FIN NUME INST_INIT quand l utilisateur renseigne l instant INST_FIN INST_INIT En effet les instants dans STAT NON LINE sont rep r s par un num ro d ordre un entier Si l utilisateur veut utiliser un instant un r el et non un num ro d ordre pour NUME INST_ l op rande pr cision permet de s lectionner ce num ro d ordre Exemple NUME 1 2 3 4 5 6 7 INST 0 0010 0 0020 0 0030 0 0040 0 0050 0 0060 0 0070 Si l utilisateur veut s lectionner l instant correspondant NUME 4 il lui suffit de dire INST 0 10000004 Par contre pour le deuxi me exemple NUME 1 2 3 4 5 6 7 INST 0 10000001 10 10000002 0 10000003 0 10000004 10 10000005 10 10000006 0 10000007 Si l utilisateur veut s lectionner l instant correspondant NUMI INST 0 10000004 car l cart relatif entre les instants vaut 0 10000004 qui est sup rieur la valeur de pr cision par d faut 1E 6 On ne pourra donc pas distinguer E 4 il ne lui suffit pas de dire 0 10000005 0 10000004 1E 7 Manuel d ut
34. ON LINE Date 18 04 2011 Page 19 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 Permet de proposer comme d placement pour la pr diction chaque pas de temps le d placement donn par une histoire m canique pr cis e sous le mot cl EVOL NOLI 13 Le d placement est projet sur l ensemble des champs cin matiquement admissible comme pour la m thode EXTRAPOLE Remarque Les m thodes EXTRAPOLE et DEPL CALCULE proc dent une projection de la solution sur l ensemble des champs cin matiquement admissibles On se sert pour cela des conditions aux limites de Dirichlet donn dans le mot clef EXCIT Dans ce cas il n est pas possible d utiliser des chargements de Dirichlet de type cin matique op rande AFFE CHAR CINE mais uniquement des chargements de Dirichlet par dualisation op rande AFFE CHAR MECA Une alarme pr vient l utilisateur dans le cas o Code Aster n aurait pas trouv de chargements de Dirichlet dualis s Le risque dans ce cas tant que le champ de d placement ne soit pas cin matiquement admissible Utilit 1 supposons qu on r alise un premier calcul avec un maillage grossier On souhaite r aliser le m me calcul mais sur un maillage plus fin On peut supposer que la solution en d placement pour ce second calcul n est pas loign e de celle du premier calcul et donc qu une bonne pr diction du d placement pour ce second calcul est la p
35. RESI GLOB MAXI sont pr sents tous les deux les deux tests sont effectu s Remarque Si les conditions limites de Dirichlet sont impos s par AFFE CHAR CINE limination et non par AFFE CHAR MECA dualisation les degr s de libert s portant ces conditions sont ignor es lors de l valuation du r sidu d quilibre Ce qui ne provoque pas de r sultats faux mais lorsque le chargement devient nul c est dire lorsque L est nul par exemple dans le cas d une d charge totale on passe du crit re de convergence relatif au crit re de convergence absolu RESI GLOB MAXI Cette op ration est transparente pour l utilisateur message d alarme mis dans le fichier mess Lorsque le vecteur L redevient diff rent de z ro on repasse automatiquement au crit re de convergence relatif RESI GLOB _ RELA 3 14 2 Op rande RESI_COMP_RELA IREST COMP RELA rescmp R Cet op rande conduit estimer la convergence de l algorithme de Newton en raisonnant composante par composante Pour cela on distingue dans le vecteur r sidu les sous vecteurs correspondant chaque composante cmp par exemple en THM emp DX DY DZ PREI PRE2 TEMP On norme ensuite ces sous vecteurs par la force interne correspondante Ainsi l algorithme continue les it rations globales tant que max JE max rbdd c 1 nbemp max IL RH H gt rescmp o Ei est la partie du r sidu Ei correspondant la composante c et ZL l
36. T 1 RECH LINEAIRE Ei METHODE CORDE DEFAUT MIXTE PILOTAGE RESI LINE RELA 1 E 1 DEFAUT reslin R ITER LINE MAXI 3 DEFAUT itelin I PAS MINI CRIT SS 0 DEFAUT pnicri R ITER LINE CRIT 20 DEFAUT itelic I RHO MIN 1 E 2 DEFAUT rmin R RHO MAX 1 E 1 DEFAUT rmax R RHO EXCL 9 E 3 DEFAUT rexc R Q PILOTAGE F TYPE DDL IMPO SAUT IMPO LONG ARC SAUT LONG ARC ANA LIM DEFORMATION PRED ELAS TOUT OUT DEFAUT GROUP MA lgrma 1 gr maille MAILLE Ima 1 maille NOEUD no noeud GROUP NO rno gr_noeud Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur S TAT_NON LINE Responsable Mickael ABBAS FISSUR NOM CM DIRE PILO COEF MULT FCO P ETA PILO R MAX ETA PILO R MIN ETA PILO MAX ETA PILO MIN PROJ_BORNES SELECTION Q SOLVEUR fiss nomcmp direpilo 1l cmult etarmax etarmin etamax etamin OUI U ONT ORM INCR DEP JUN JUN ANGL INCR DEP RESIDU _F voir le document 0
37. cements aux n uds pris l tat initial Si l un de ces champs n est pas pr cis il est pris nul par d faut Ils peuvent par exemple tre issus de la commande CREA CHAMP o bien avoir t lus dans un fichier au format I DEAS ou MED par la commande LIRE RESU 3 6 2 Op randes EVOL NOLI EVOL NOLI evol Nom du concept de type evol poli d o sera extrait l tat initial 3 6 3 Op rande NUME ORDRE INST NUME DIDI NUME ORDRE nuini INST instini Extraction de l tat m canique initial dans evol partir du num ro d archivage NUME ORDRE ou de l instant d archivage INST pour effectuer la poursuite du calcul Si NUME ORDRE ou INST ne sont pas remplis on prend le dernier num ro archiv existant dans evol NUME DIDI nudidi Dans le cas de chargements de type Dirichlet diff rentiel DIDI on donne sous NUME_DIDI le num ro d archivage de l tat m canique d placements qui sert de r f rence pour l application de ces conditions aux limites Par d faut on prend l tat m canique d fini sous NUME ORDRE ou INST 3 6 4 Op rande INST_ETAT_INIT INST ETAT INIT istetaini On peut associer une valeur d instant istetaini cet tat initial Par d faut 1 Lorsque l tat initial est d fini par la donn e des champs ETAT INIT avec DEPL SIGM VARI il n y a pas d instant associ 2 Lorsque l tat est donn
38. ct ristiques l mentaires K8 CHAMPMAT CHAM MATER Champ de mat riau K8 PARM THETA COMP INCR PARM THETA Param tre d int gration de la loi de comportement R ITER GLOB Nombre total d it rations de Newton l CHAR MINI Chargement minimum atteint au cours du pas de temps R ETA PILOTAGE Param tre de pilotage R Quand on recherche des modes de flambement ou des modes vibratoires on stocke le champ de d placement correspondant et la valeur du chargement critique ou la fr quence Nom Mot clef origine Description Type CHAR CRIT CRIT FLAMB Chargement critique du mode de flambement R FREQ MODE VIBR Fr quence du mode vibratoire E 3 22 Op rande INFO 9 INFO inf Permet d effectuer dans le fichier message diverses impressions interm diaires D autres impressions sont faites syst matiquement lors du calcul non lin aire ind pendamment de la valeur affect e au mot cl INFO ce sont les impressions des r sidus et des incr ments relatifs de d placement au cours des it rations de Newton Attention les fichiers mess peuvent devenir tr s importants avec INFO 2 3 23 Op rande TITRE O TITRE Ce Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Water default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 37 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 61
39. de s lectionner les isntants c f U4 71 00 3 15 3 Op rande CHAR CRIT CHAR CRIT 10 10 DEFAUT intcc Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 31 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 Le mot cl CHAR CRIT permet de gagner du temps en ne faisant qu un test de Sturm dans la bande de fr quence fournie Si on trouve au moins une fr quence alors on calcule r ellement les valeurs des charges critiques dans cet intervalle 3 15 4 Op rande NB FREQ NB FREQ 3 DEFAUT nbfreq Le mot cl NB_FREQ 3 par d faut d signe le nombre de charges critiques calculer En fait seule la premi re suffit mais il peut y avoir des modes multiples 3 15 5 Op rande RIGI_GEOM RIGI GEOM OUI DEFAUT NON Le mot cl RIGI GEOM OUI par d faut donne le choix l utilisateur entre effectuer une recherche de valeurs propres g n ralis es avec la matrice g om trique au second membre ou non cas des grandes d formations Choisir NON signifie que la matrice de raideur g om trique est remplac e par l identit 3 15 6 Op rande DDL EXCLUS DDL EXCLUS DX DY Le mot cl DDL EXCLUS liste vide par d faut d signe l ensemble des degr s de libert qu
40. diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 30 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 3 15 3 15 1 Si un des crit res de convergence globale choisis n est pas v rifi apr s maglob it rations alors le programme s arr te les r sultats pr c dents sont sauvegard s Q ARRET NON DEFAUT Si maglob est insuffisant pour v rifier les crit res de convergence donn s par l utilisateur on passe quand m me l instant suivant Cette option n est utilsiable qu en mode DEPL_CALCULE Cette option est utiliser avec pr caution car elle donne des r sultats faux Mot cl CRI T_FLAMB CRIT FLAMB EU Ce mot cl permet de d clencher le calcul la fin de chaque incr ment de temps d un crit re de stabilit Ce crit re est utile pour d celer au cours du chargement le point partir duquel on perd la stabilit par flambage par exemple Ce crit re est calcul de la fa on suivante la fin d un pas de temps en petites perturbations on r sout det KI A K 0 KT est la matrice tangente coh rente cet instant JK est la matrice de rigidit g om trique calcul e partir du champ de contraintes cet instant En pratique le chargement est instable si A lt en fait 1 lt A lt 0 On calcule les valeurs propres par la m thode de Sorensen cf MODE
41. e d un n ud ou groupe de n ud sur lequel on peut imposer diff rents modes de pilotage mot cl TYPE Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 22 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 Attention Avec FIXE PILO on ne peut pas utiliser pour le chargement de r f rence le mot cl FONC MULT Lorsque le chargement de r f rence est d fini par AFFE CHAR MECA F ce chargement peut tre fonction des variables d espace mais pas du temps De m me les changements issus de variables de commande comme la temp rature qui d pendent du temps ne sont pas utilisables avec le pilotage eLe mot cl PILOTAGE est interdit avec le contact en formulation hybride d conseill avec le contact en formulation p nalis e 3 12 1 Op rande TYPE T YPE DDL IMPO LONG ARC ANA LIM DEFORMATION PRED ELAS S S AUT IMPO AUT LONG ARC C est le type de pilotage effectu Sept modes de pilotage sont disponibles Cf R5 03 80 pour plus de d tails DDL IMPO Permet d imposer une valeur donn e d incr ment de d placement une seule composante possible en un unique n ud no ou d un groupe de n uds ne comportant qu un seul n ud chaque incr ment de temps on cherche l am
42. e l on souhaite mettre 0 dans le second membre de la recherche de valeurs propres g n ralis es Il ne peut tre utilis que sous la condition RIGI GEOM NON Cela permet d imposer des conditions suppl mentaires de compatibilit sur les modes propres et ainsi d effectuer une recherche s lective Cela est particuli rement adapt aux formulations mixtes Dans ce cas l limination des multiplicateurs de Lagrange au second membre permet d exclure les modes parasites dominantes Lagrangiennes et de valeurs propres n gatives 3 16 Mot cl SENSIBILITE La syntaxe de ce mot cl commun plusieurs commandes est d crite dans le document U4 50 02 3 17 Mot cl ARCHIVAGE Q ARCHIVAGE EI Permet d archiver des ou certains r sultats tous ou certains instants du calcul En l absence de ce mot cl tous les pas de temps sont archiv s y compris les instants de calculs nouvellement cr s par red coupage automatique du pas de temps L archivage permet de r duire sensiblement la taille des bases en s lectionnant les instants sauvegard s 3 17 1 Op rande LIST _INST I INST PAS ARCH Q LIST INST list r8 INST 1 r8 PAS ARCH npas La d signation des instants stocker est effectu e soit par une liste d instants List r8 ou 1_r8 ou alors par une fr quence d archivage tous les npas de temps En l absence de ces mots cl s tous les pas de temps sont archiv s Deux remarques 1 le dernier pas de calcul est toujours stock
43. e vecteur des forces internes au temps n correspondant cette m me composante c Cf R5 03 01 pour plus de d tails Les forces internes au temps n sont calcul es en d but de pas de temps partir du r sultat issu du pas de temps pr c dent Pour le premier pas de temps on passe automatiquement un crit re relatif de type RESI GLOB RELA Voir un crit re absolu pour les cas o le chargement est nul Ce choix n a d int r t que pour des probl mes de type volutif THM o r sident de forts contrastes entre les diff rentes inconnues 3 14 3 Op rande RESI REFE RELA RESI REFE RELA resref R SIGM _ REFE sigref R FORC _ REFE forref R VARI REF varref R EPSI REFE epsref R FLUX THER REF fthref R FLUX HYD1 REFE fhiref R FLUX HYD2 REFE fh2ref R DEPL REFE depref R LAGR REFE lagref R Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Ast oh ogae_Aster default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 29 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 Cet op rande conduit estimer la convergence de l algorithme de Newton de la mani re suivante A partir d une r f rence qui peut tre Une contrainte sigref Une d formation epsref e Une variab
44. e version diff rente de la d finition de I est utilis e o ce param tre est li l incr ment d endommagement voir R7 01 04 cmu Utilisation Attention Lorsqu on veut utiliser ces deux derniers modes de pilotage il est indispensable de faire un premier STAT NON LINE sans le mot cl PILOTAGE pour amorcer le probl me et obtenir un tat initial lt diff rent de z ro sinon division par z ro pour le pilotage par incr ment de d formation On effectue apr s une reprise partir de cet tat initial non nul et on utilise le pilotage De plus la r solution des deux quations pr c dentes permet d obtenir l intensit du chargement inconnue Dans certains cas la r solution de ces quations peut conduire plusieurs solutions pour l intensit On choisit alors toujours la solution qui est la plus proche de C est pourquoi lorsqu on veut imposer un chargement altern on est oblig chaque changement de signe du chargement de r aliser un premier STAT NON LINE sans le mot cl PILOTAGE afin d obtenir un Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Katie Titre Op rateur S TAT_NON LINE Date 18 04 2011 Page 24 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 El tat initial de traction ou de compression On effectue ens
45. ence Par d faut tous les champs sont identiquement nuls L tat initial peut tre d fini soit en pr cisant chaque champ de l tat initial soit en extraction depuis un concept de type evol noli pr existant La donn e d un tat initial n a de sens et n est donc prise en compte que pour la partie du domaine trait e en comportement incr mental COMP_INCR si le comportement est lastique COMP ELAS cela n a aucune incidence Si l on veut prendre en compte un tat initial en lasticit c est le mot cl ELAS situ sous COMP _INCR qu il faut utiliser Remarques Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster au default Titre Op rateur STAT_NON_LINE Date 18 04 2011 Page 13 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 e Dans le cas o l utilisateur a sp cifi que le concept r sultat est r entrant par le mot r serv reuse le mot cl ETAT_INIT est obligatoire Dans le cas o l on utilise la m thode continue du contact la reprise de calcul peut donner lieu des difficult s de convergence du fait de l oubli de l tat de contact pr c dent 3 6 1 Op randes SIGM VARI DEPL SIGM sig VARI vain DEPL depl sig est le champ de contraintes aux points de Gauss vain est le champ des variables internes aux points de Gauss et dep1 est le champ des d pla
46. et la matrice tangente coh rente option FULL MECA est r valu e toutes les it it rations de Newton pour un incr ment de temps donn pr cis ment aux it rations de num ro it 2it 3it Donc la premi re it ration de Newton on ne r assemble la matrice tangente que si it vaut 1 sinon on garde la matrice utilis e dans la phase de pr diction Par convention si it vaut O la matrice n est pas r valu e durant tout le pas de temps Q PAS MINI ELAS 0 DEFAUT pasmini R Q REAC ITER ELAS 0 DEFAUT it I Ces options permettent de passer de la matrice tangente la matrice de d charge i e en consid rant que les non lin arit s n voluent pas lorsque le pas de temps est inf rieur pasmini Cette matrice de d charge est la matrice lastique pour les mod les de comportement de type plastique pour les mod les d endommagement elle s identifie la matrice s cante Comme la convergence avec la matrice lastique est plus lente que celle avec la matrice tangente le mot cl ITER GLOB ELAS sous le mot cl facteur CONVERGENCE permet de d finir un nombre d it rations maximal sp cifique l utilisation de la matrice lastique et diff rent de celui associ l utilisation de la matrice tangente Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster default Titre Op rateur STAT NON LINE Date 18 04 2011
47. eur STAT NON LINE Date 18 04 2011 Page 21 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 La m thode MIXTE est plus labor e et utilise une m thode s cante avec des bornes variables Elle est plus efficace lorsque la fonctionnelle n est pas strictement concave probl mes avec endommagement ou THM par exemple La m thode PILOTAGE est r serv e au pilotage de type DEFORMATION PRED ELAS et LONG ARC voir 821 C est la seule m thode utilisable avec ce type de pilotage Pour le pilotage de type DDL IMPO on peut utiliser CORDE ou MIXTE 3 11 2 Op rande RESI LINE RELA ITER LINE MAXI Q RESI LINE RELA 1 E 1 DEFAUT reslin ITER LINE MAXI 3 DEFAUT itelin Ce sont les param tres de la recherche lin aire On donne le nombre d it rations maximum itelin effectuer et la pr cision reslin atteindre pour r aliser la convergence de la recherche lin aire Il est conseill de ne pas utiliser la recherche lin aire avec du contact Pour la m thode CORDE Il n est pas n cessaire de sp cifier une pr cision ni un nombre d it rations tr s lev s la pratique montrant que deux ou trois it rations de recherche lin aire sont suffisantes On peut donc se contenter de demander trois it rations avec la pr cision par d faut Vous ne pouvez mettre plu
48. ialisation Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur STAT NON LINE Responsable Mickael ABBAS A Exemple simple comportement par d faut LIST DEFI LIST RE Premier STAT NON LINE Second STAT NON LINE correspondant au temps 4s par d faut INST INIT INST STAT NON STAT NON INE INE EL DEBUT 0 INTERVALLE Date 18 04 2011 Page 17 37 Cl U4 51 03 II INCREMENT F LIST INST LIST INST FIN 4 reuse INCREMENT F LIST INST LIST ETAT INIT B Exemple pour montrer l int r t de INST_INIT LIST DEFI LIST RE Premier STAT NON LINE Second STAT NON LINE STAT NON L STAT NON INE effectue le calcul des instants 1 Ae effectue le calcul pour les instants 9 et 10s ne fait rien pour t 5 6 7 et 8s l tat initial correspondant au temps 4s par d faut INST 4 FF effectue le calcul pour les instants 1 2 3 et 4s effectue le calcul pour les instants 5 6 7 8 9 et 10s l tat initial ETAT INIT INST 4 EVOL NOLI EL DEBUT 0 INTERVALLE F JUSQU A 10 N
49. ign le tableau est construit suivant les diff rentes options de calcul recherche lin aire pilotage contact etc et avec INFO RESIDU NON 3 18 1 Op rande UNITE UNITE unit Le tableau de convergence sera dupliqu dans le fichier d unit unit au format csv le s parateur tant la virgule 3 18 2 Op rande INFO RESIDU INFO RESIDU NON DEFAUT OUI Cet op rande permet d ajouter une colonne pour or r sidu valu RESI GLOB RELA RESI GLOB MAXI RESI COMP RELA et RESI REFE RELA Cette colonne indiquera le n ud o le r sidu est maximum ce qui peut aider l utilisateur lorsqu il y a des difficult s de convergence Par exemple pour voir si le mat riau a t mal d fini avec une valeur incorrecte sur un l ment 3 19 Mot cl OBSERVATION Q OBSERVATION EI Ce mot cl permet de post traiter certains champs aux n uds ou aux l ments sur des parties de mod le des instants d une liste dite d observation g n ralement plus raffin e que la liste des instants archiv s d finie dans le mot cl ARCHIVAGE 31 o on stocke tous les champs sur tout le mod le II sert essentiellement des conomies de stockage mais aussi valuer des champs sur des parties r duites du maillage sans avoir besoin de post traiter apr s le calcul Il est possible par exemple de calculer la norme des contraintes a
50. ilisation Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Fascicule u4 51 Analyse statique Code Aster default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 18 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 NUME 3 4 et 5 le code s arr te alors en erreur fatale Il suffit alors de changer le param tre PRECISION pour pouvoir s lectionner l instant dans l exemple PRECISION 1E 8 conviendra 3 7 4 Op rande ERRE TEMPS Q ERRE TEMPS NON DEFAUT OUI Cet op rande permet d activer le calcul des indicateurs d erreur temporel pour les mod lisations HM instationnaires Voir R4 10 05 3 8 Op rande METHODE Q METHODE NEWTON IMPL EX Permet de choisir la m thode de r solution du probl me incr mental non lin aire NEWTON On utilise l algorithme de Newton Raphson pour r soudre le probl me voir R5 03 01 IMPL EX On utilise l algorithme IMPL EX pour r soudre le probl me voir R5 03 81 3 9 Mot cl NEWTON Pr cise les caract ristiques de la m thode de r solution du probl me incr mental non lin aire m thode de Newton Raphson 3 9 1 Op rande PREDICTION Q PREDICTION TANGENTE ELASTIQUE EXTRAPOL DEPL_CALCULE La phase de pr diction Cf R5 03 01 a pour but de calculer une estimation du champ de d placements afin de permettre la m thode
51. le interne varref si l on utilise des lois non locales gradient de d formation Un flux thermique fthref dans un cas THM e Deux flux hydriques fh1iref et fh2ref dans un cas HHM Un d placement depref si on utilise des l ments de joint avec un comportement de type CZM e Une force ou un moment g n ralis forref si on utilise des l ments de structure discrets barres poutres ou c bles e Un coefficient de Lagrange lagref pour les formulations mixtes en endommagement On calcule une r f rence de r sidu EI un vecteur de m me longueur que le vecteur r sidu La convergence sera r alis e si et seulement si Vie 1 nbadl F lt resref ET 3 14 4 Op rande ITER GLOB MAXI ITER GLOB MAXI 10 DEFAUT maglob Nombre d it rations maximum effectu pour r soudre le probl me global chaque instant 10 par d faut Ce test est toujours effectu sauf dans le cas du red coupage du pas de temps par la m thode EXTRAPOLE L augmentation excessive de ce param tre est g n ralement le signe d un probl me dans la mod lisation ou d une discr tisation temporelle inad quate 3 14 5 Op rande ITER_GLOB_ELAS ITER GLOB ELAS 25 DEFAUT maxelas Nombre d it rations maximum effectu avec la matrice lastique lorsqu on utilise le mot cl PAS MINI ELAS du mot cl facteur NEWTON voir 18 pour r soudre le probl me global chaque instant 25 par d faut On rappelle
52. ouvoir s lectionner l instant dans l exemple PRECISION 1E 8 conviendra INST 0 10000004 car l cart relatif entre les instants vaut 3 7 Mot cl INCREMENT INCREMENT EI D finit les intervalles de temps pris dans la m thode incr mentale Les instants ainsi d finis n ont de sens physique que pour des relations de comportement o le temps intervient explicitement visco lastiques ou visco plastiques par exemple Dans les autres cas ils permettent seulement d indicer les incr ments de charge et de param trer l volution d un ventuel champ de temp rature 3 7 1 Op rande LIST_INST LIST_INST litpsr8 listr8 litps list inst Si LIST INST litpsr8 listr8 Les instants de calcul sont ceux d finis dans le concept 1itpsr8 par l op rateur DEFI LIST REEL U4 34 01 Si LIST INST litps list inst Les instants de calcul sont ceux d finis dans le concept litps par l op rateur DEFI LIST INST U4 34 03 3 7 2 Op randes NUME INST_INIT INST_INIT NUME INST EIN INST_FIN NUME_INST_INIT nuini INST_INIT instini L instant initial du calcul qui donc n est pas re calcul est d sign soit par sa valeur INST_INIT soit par son num ro d ordre dans la liste d instants 1itps NUME INST INIT Pour pouvoir acc der par valeur il est n cessaire que la liste soit ordonn e En l absence des mots cl s INST_INIT ou
53. plitude n du chargement de r f rence qui permettra de satisfaire la relation incr mentale suivante cmult Au no At SAUT_IMPO Reprend le principe de DDL IMPO mais pour contr ler l incr ment du saut de d placement entre les l vres d une fissure X FEM Une seule direction est possible mais elle peut tre d finie dans une base locale normale ou tangente la fissure On contr le la moyenne de cet incr ment de saut sur un ensemble de points d intersection P de l interface avec les ar tes a du maillage Cet ensemble d crit toute la fissure si GROUP NO n est pas renseign comportement par d faut et seulement une partie s il l est Attention ce type de pilotage ne peut tre utilis qu en mod lisation X FEM 1 N ema 2 Au P a 1 At LONG ARC Permet de piloter l intensit n du chargement de r f rence par la longueur abscisse curviligne de la r ponse en d placement d un groupe de n uds utiliser par exemple lorsqu on veut contr ler le flambement d une prouvette On v rifie la relation suivante XX Au emt La ul A avec A u o n sont les n uds du pilotage et c les composantes du d placement des n uds consid r s M me si le groupe de n ud du pilotage est r duit un seul n ud il faut quand m me utiliser GROUP NO SAUT_ LONG ARC Reprend le principe de LONG ARC mais pour contr ler la norme de l incr ment du saut de d placement entre les l vres d une fi
54. que PAS_MINI_ELAS permet de passer de la matrice tangente la matrice lastique lorsque le pas de temps est ou devient par le red coupage inf rieur une certaine valeur pr cis e sous PAS MINI ELAS Contrairement ITER GLOB MAXI ce param tre peut facilement prendre des grandes valeurs plusieurs centaines car la convergence sur un probl me non lin aire avec la matrice lastique tr s raide est lente bien qu assur e du point de vue th orique pour toutes les lois d crivant les mat riaux standards g n ralis s 3 14 6 Op randes TYPE PLATEAU ITER PLATEAU RELA TYPE PIC DEFAUT PLATEAU PLATEAU ITER 3 DEFAUT plaite 1 PLATEAU RELA 1E 3 DEFAUT plarel R Cet op rande permet de contr ler le type d op rateur de convergence appliquer pour les mots clefs RESI Par d faut le mode PIC assure qu il y a convergence d s lors que la valeur seuil donn e par RESI est atteinte Le mode PLATEAU est plus s v re et exige que le crit re de convergence soit stable pendant plaite it rations autour dans un tunnel de largeur plarel autour de la valeur de r f rence donn e par RESI Ce mode est utile dans les mod lisations THM pour lesquelles le crit re de Newton est parfois pas assez s v re 3 14 7 Op rande ARRET Q ARRET OUI DEFAUT Manuel d utilisation Fascicule u4 51 Analyse statique Document
55. res D IT 1 2 OVO O ac ES DE ERP nn EDP aus 23 5 J OporandES enine iia Eege 10 3 1 Op randes MODELE CHAM_MATER CARA _ELEM snesen 10 S 2 MOLCI E XO HEEN 10 3 2 9 Op rande TYPE E 11 3 9 Oe CR EE 11 SE MOECI SOUS E 12 KE Maes ee e EE 12 3 42 Op randes TOUT SUPER MAILLE 28 uu ententes ane aiia 12 9 49 Operande FONC MULT uuu sement mesure 12 3 5 Mots cl s COMP INCR et COMP ELAS gesend 12 LO ee ETAT N S nee Rene as eee ua ti te Dem 12 3 6 1 Op randes SIGM VARMDEPE u uquna eege idee 13 36 2 Op randes EVOL NO 0 s ahua 13 3 6 5 Op rand PRECISION CRITERE uuu uuu luas asus aasssqbaqaaaqaawaqaqhanapasmapuqakaq ataqa hhiaaq 14 S ege e TC ME 15 S F 1 OP rANTE HST INST nana dede ten nana ue nee uen 15 3 7 2 Op randes NUME_INST_INIT INST_INIT NUME_INST FIN INST FIN 15 RETTEN E e Eed EE 16 S a Operando E EEN TEE 17 S S Op rand METIPEIODE E 17 S S Ee CE NEWTON BEE 17 39 T Operangde BREDIGTION uu EE 17 3 9 2 Op rand e MATRIC E uu uuu uuu auquaaskesadaqanapaptspahaqeaaqsqaqquqaqananaabal aqihakhua aqa aubskabuataunaua 18 S 9 3 Operaride Ee EE 19 O Motele IMPL EX eege EE 19 3 10 1 0p rande MATRICE EE 19 See PINEAIRENS u S uuu Ee 19 2 14 11 Op rande METHODE Q
56. rojection des d placements du premier calcul sur les n uds du nouveau maillage la projection des d placements sur le nouveau maillage doit tre r alis e pr alablement avec l op rateur PROJ CHAMP U4 72 05 Ce mot cl permet de r aliser ce mode de pr diction 2 cela permet de r duire la place m moire et de conserver ces r sultats en vue d une poursuite ult rieure Pour un gros calcul on peut stocker uniquement les d placements tous les instants aux formats IDEAS ou MED dans IMPR RESU Si on veut recalculer les contraintes et variables internes on fait un LIRE RESU au format ad quat puis on utilise DEPL CALCULE avec ITER GLOB MAXI 0 on effectue une seule it ration et ARRET NON il n y a pas convergence on ne v rifie pas l quilibre Il est toutefois n cessaire pour des raisons de syntaxe de donner un chargement viter les chargements Dirichlet qui imposent une r solution lin aire ainsi qu un crit re de convergence m me si ces informations ne sont pas prises en compte 3 9 2 Op rande MATRICE MATRICE TANGENTE REAC INCR 1 DEFAUT 8 mf 9 REAC ITER 0 DEFAUT it La matrice utilis e pour les it rations globales de la m thode est la matrice tangente R5 03 01 La matrice tangente de pr diction est r valu e tous les mf incr ments de temps mf positif ou nul
57. s de 999 it rations de recherche lin aire pour la m thode CORDE Par contre pour la m thode MIXTE sur des probl mes avec endommagement plusieurs dizaines d it rations sont souvent efficaces 3 11 3 Op rande PAS_MINI_CRIT I ITER_LINE CRIT Q PAS MINI CRIT 0 DEFAUT pmicri R Q ITER LINE CRIT 20 DEFAUT itelic I Lors de pas de temps o la convergence est delicate on peut vouloir augmenter le nombre maximum d it rations de recherche lin aire C est ce que permettent les mots cl s PAS MINI CRIT et ITER LINE CRIT Quand le pas de temps directement fix par l utilisateur ou cons quence de d coupages de pas de temps devient inf rieur la valeur pmicri le nombre d it rations de recherche de recherche lin aire passe de itelin renseign par ITER LINE MAXI itelic renseign par ITER LINE CRIT 3 11 4 Op randes RHO MIN RHO MAX RHO EXCL Q RHO MIN 1 E 2 DEFAUT rmin R Q RHO MAX 1 E 1 DEFAUT rmax R Q RHO EXCL 9 E 3 DEFAUT rexc R Ces mots cl s fixent l intervalle 1 dans lequel on calcule le coefficient RHO de la recherche lin aire sous la forme 7 rmin rmax rexc rexc R5 03 01 3 12 Mot cl PILOTAGE Q PILOTAGE F Lorsque l intensit n d une partie du chargement n est pas connue a priori chargement dit de r f rence d fini dans AFFE CHAR MECA ou AFFE CHAR MECA F avec charge de type FIXE PILO le mot cl PILOTAGE permet de piloter ce chargement par l interm diair
58. s tous les cas Le principe de fonctionnement est le suivant chaque it ration de Newton si l on trouve des valeurs de pilotage dans l intervalle etamin etamax toutes les valeurs de pilotage en dehors de cet intervalle sont ignor es En revanche si aucune valeur de pilotage n est trouv e dans cette intervalle toutes les valeurs de pilotage sont conserv es et c est l aide de mot cl PROJ BORNE qu on d cide comment proc der dans la suite Attention Avec la loi ENDO_ISOT_BETON ces deux mots cl s sont obligatoires 3 12 100p rande PROJ _BORNES PROJ BORNES OUI DEFAUT 8 NON En cas de d passement de l intervalle etamin etamax l utilisateur peut indiquer s il veut projeter la valeur de pilotage sur etamin etamax Avec PROJ BORNE OUI la projection sera effectu e si eta gt etamax alors eta etamax si eta lt etamin alors eta etamin ce qui permet en cas de convergence d arr ter le calcul pr cis ment sur etamin ou etamax Avec PROJ BORNE NON on ne modifie pas les valeurs de eta m me si pendant les it rations de Newton cette derni re a une valeur sup rieure etamax ou inf rieure etamin Par contre le calcul est arr t si la convergence eta d passent les bornes Une utilisation possible de l intervalle etamin etamax avec l option PROJ BORNE OUI est le suivant On d sire par exemple comparer plusieurs calcul pour un mod le adoucissant
59. ssure X FEM On contr le cette norme en moyenne sur un ensemble de points d intersection P de l interface avec les ar tes a du maillage Cet ensemble d crit toute la fissure si GROUP NO n est pas renseign comportement par d faut et seulement une partie s il l est Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Ast ia ogae_Aster default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 23 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 cmult All A avec Au A Gd Au lp o c sont les composantes du d placement M me si le groupe de nod du pilotage est r duit un seul n ud il faut quand m me utiliser GROUP NO ANA LIM Ce mode de pilotage est sp cifique au calcul de charge limite loi NORTON HOFF par approche cin matique cf R7 07 01 pour plus de d tail Si F d signe le chargement assembl pilot TYPE CHARGE FIXE PILO alors la fonction de pilotage s crit simplement P u F u 1 Except pour le calcul de charge limite cette fonctionnalit ne pr sente pas d int r t a priori Pour ce mode de pilotage aucun autre mot cl n est pr ciser L utilisation de lois de comportement adoucissantes peut conduire des snap backs brutaux qui rendent d licat le d roulement du calcul Les deux modes de pilotage suivants y rem dient Cf R5 03 80 pour plus de d
60. suivantes si deux noeuds sont les extr mit s respectives de deux ar tes non ind pendantes une seule sera retenue fig 3 12 2 2 si un noeud est extr mit de plusieurs ar tes on retient arbitrairement la premi re rencontr e par l algorithme si deux noeuds sont extr mit s d une m me fissure fig 3 12 2 3 une erreur sera renvoy e D une fa on g n rale il est conseill que tous les noeuds entr s soient du m me c t de la fissure si un noeud ne correspond aucune ar te fig 3 12 2 4 une erreur est renvoy e s 3 Noeuds utilisateur S lection automatique Figure 3 12 2 1 GROUP_NO non renseign S lection utilisateur Figure 3 12 2 2 Noeuds extr mit s d ar tes non ind pendantes Figure 3 12 2 3 erreur pour noeuds connect s la m me ar te Figure 3 12 2 4 erreur pour noeud non connect une ar te intersect e Tableau 3 1 Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster goe Titre Op rateur STA T NON UNE Date 18 04 2011 Page 25 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 3 12 3 Op randes TOUT MAI LLE GROUP MA TOUT OUT DEFAUT GROUP MA lgrma MAILLE lma 3 12 4 3 12 5 On donne les mailles ou groupes de mailles servant piloter le calcul A n utiliser qu a
61. u sens de Von Mises et de la stocker dans la table d observation Ce mot cl est r p table et permet la cr ation d une table d observation de m me nom que le concept r sultat de STAT NON LINE que l on pourra extraire l aide de la commande RECU TABLE On ne peut r aliser que 99 observations au maximum c est dire qu il y a au maximum 99 lignes dans la table chaque pas de temps observ 3 19 1 Op randes LIST _INST I INST PAS OBSE Q LIST INST list r8 INST 1 r8 PAS OBSE npas Ces op randes permettent de d finir aux choix une liste d instants d observation LIST INST INST et PAS OBSE ont la m me signification que les op randes de m me nom servant d finir une liste d archivage PAS OBSE jouant le m me r le que PAS ARCH dans ARCHIVAGE 31 3 19 2 Op randes PRECISION CRITERE Q PRECISION prec 9 CRITERE ABSOLU RELATIF Cf U4 71 00 pour la syntaxe d taill e Ces param tres permettent de g rer la pr cision de la s lection des instants pour l observation 3 19 3 Op randes NOM_CHAM NOM CMP NOM CHAM nomcham NOM CMP nomcmp Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Eeer default Titre Op rateur STAT NON LINE Date 18 04 2011 Page 34 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 Ces op randes permettent de d finir le champ post traiter
62. u uuu uuu uksa aa inaa anaa Aaaa 19 3 11 2 Op rande RESI LINE _ RELA7ITER LINE MAXI u u Uu assasi aq uawasbapaq ua ka 20 3 11 3 Op rande PAS _MINI_CRIT ITER_LINE_CRIT uu uuu ul uuu ennnnneeneennnnens 20 3 11 4 Op randes RHO_MIN RHO_MAX RHO_ EXCL rrasssssssssssssssssssssa 20 3 12 MOD Cl PILOTAGE dianian ironia aaie aaa EEA eet 20 312 Op rande EE 21 Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster W Titre Op rateur STAT_NON _LINE Date 18 04 2011 Page 3 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 3 12 2 Op randes NOEUD GROUP EE 22 3 12 3 Op rand s TOUT MAILLE GROUR MA L u susss2tedtcacanenoieiemennermuimtagnenarmmengiueres 23 3 124 Op rande Le Me EE 24 120 p rande DIRE E EE 24 3 12 6 Op rande FISSURE EE 24 2 7 Oporande COEF MILITE RSS ete aea nr 24 3 12 8 Op rande ETA PILO R_ MAX ETA PILO BR MIN een 24 Se deel 26 3 14 MOCCI CONVERGENCE EE 26 3 14 1 Op rande RESI GLOB RELA RESI GLOB MAXI 2 5 6 nnansrnns semis 26 3 142 Op rnde RESI GOMP RELA uuu de thin bahasa kaba ege 27 3 143 Op rande RT BEE 27 ST Open MER GLOB MAA nine 28 3 14 5 Op r nde ITER GLOB EE EE 28 3 14 6 Op randes TYPE PLATEAU_ITER PLATEAU_REL
63. uite un second STAT NON LINE en poursuite partir de l tat initial pr c dent avec le mot cl PILOTAGE z Remarque DEFORMATION et PRED ELAS ne sont pas disponibles pour les l ments de structures 3 12 2 Op randes NOEUD GROUP_NO NOEUD no GROUP NO grno A n utiliser qu avec DDL IMPO LONG ARC SAUT IMPO ou SAUT LONG ARC Pour DDL 1IMPO si on utilise l op rande GROUP NO le groupe de n uds en question ne doit contenir qu un seul n ud Dans les autres cas on utilise uniquement GROUP NO qui peut ventuellement ne contenir qu un seul n ud Pour SAUT _IMPO et SAUT LONG ARC l op rande est facultative Pour DDL IMPO et LONG ARC on donne le nom du n ud ou le nom du groupe de n uds sur lequel on va imposer le pilotage Pour SAUT IMPO et SAUT LONG ARC la d finition est plus subtile puisqu en mod lisation X FEM on ne pilote pas les valeurs sur des noeuds mais sur des points d intersection entre les ar tes du maillage et la fissure Dans la suite on d signe simplement par ar tes les ar tes intersect es L algorithme commence par construire un ensemble d ar tes ind pendantes qui couvre toute la fissure voir fig 3 12 2 1 Par d faut il pilote sur toutes ces ar tes Le mot cl GROUP NO permet l utilisateur de restreindre cet ensemble chaque noeud renseign correspondant alors l extr mit d une ar te que l on souhaite piloter Signalons alors les r gles
64. ul ne convergera pas si on est proche de la charge limite Attention la signification de COEF_MULT pour le pilotage de type PRED_ELAS Op rande ETA PILO R MAX ETA PILO R MIN ETA PILO R MAX etarmax R o ETA PILO R MIN etarmin R Ces deux mots cl s permettent de d finir l intervalle de recherche des valeurs de pilotage A chaque it ration de Newton toutes les valeurs de pilotage en dehors de etarmin etarmax sont ignor es Ceci peut emmener chec de pilotage si cet intervalle est trop restrictif Si on ne pr cise pas de valeurs c est co pour etarmin et pour etarmax Une utilisation possible de cet intervalle est le suivant On d sire par exemple piloter une pression impos e la structure et on s attend garder cette pression positive En fixant etarmin 0 cela permet d imposer les valeurs de pilotage positives Op rande ETA PILO MAX ETA PILO MIN ETA PILO MAX etamax R ETA PILO MIN etamin R Ces deux mots cl s permettent de pr ciser l intervalle de valeurs de pilotage souhait sans toutefois l imposer On l utilise soit pour guider le calcul lorsque plusieurs valeurs de param tre ETA PILOTAGE satisfont les quations de pilotage soit pour arr ter le calcul lors qu une d elles atteint les bornes de l intervalle L intervalle souhait doit tre plus restrictif que l intervalle de recherche d fini pr c demment comme ce derni re est appliqu dan
65. utilis pour le pilotage Par d faut cette valeur vaut 1 A ne pas utiliser avec ANA LIM Exemple avec DDL_ IMPO Supposons que l on veut conna tre la charge limite d une structure Le chargement impos sur la structure est la pression d intensit inconnue P n valeur de r f rence P sur le groupe de maille A Pour trouver la charge limite n on va piloter le d placement du n ud NO On veut que le d placement final suivant x de ce n ud soit gal 2 soit d apr s la liste d instants des pas de 0 2 soit un coefficient cmult 1 0 2 S5 PRESSION AFFE CHAR MECA PRES GROUP MA A PX 1 0 LIST DEFI LIST REEL DEBUT 0 INTERVALLE F JUSQU A 10 NOMBRE 10 RESU STAT NON LINE EXCIT F CHARGE PRESSION TYPE CHARGE FIXE_PILO PILOTAGE F TYPE DDL_IMPO NOEUD NO1 NOM CMP Lt COEF MULT 5 Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Es Titre Op rateur STAT NON LINE Date 18 04 2011 Page 26 37 Responsable Mickael ABBAS Cl U4 51 03 R vision 6167 3 12 8 3 12 9 Dans le fichier resu la valeur de n sera affich e chaque instant du calcul Pour conna tre la charge limite il suffit de faire P n xP lci P vaut 1 donc on a directement la charge limite Si on impose sur la structure une pression P proche de la charge limite sans utiliser le pilotage le calc
66. vec DEFORMATION ou PRED ELAS Int ressant pour all ger la r solution des quations de ces trois modes de pilotages Op rande NOM CMP 9 NOM CMP nomcmp C est le nom de la composante correspondant au degr de libert utilis e pour le pilotage DX par exemple A n utiliser qu avec DDL IMPO ou LONG ARC Op rande DIRE PILO 9 DIRE PILO direpilo C est le nom de la direction selon laquelle on contr le le saut de d placement Les valeurs possibles sont DX DY DZ DNOR pour la normale la fissure DTAN pour la premi re tangente produit vectoriel de la normale avec X DTAN2 pour la deuxi me tangente A n utiliser qu avec une mod lisation X FEM Utilisation pour les types SAUT IMPO SAUT LONG ARC ou avec PRED ELAS si la s lection sur le choix de la solution pilot e est ANGL INCR DEPL ou NORM INCR DEPL 3 12 6 Op rande FISSURE FISSURE fiss Nom de la sd fiss xfem A n utiliser qu avec une mod lisation X FEM Utilisation pour les types SAUT IMPO SAUT LONG ARC ou avec PRED ELAS si la s lection sur le choix de la solution pilot e est ANGL INCR DEPL ou NORM INCR DEPL 3 12 7 Op rande COEF MULT COEF MULT cmult C est la valeur not e Cut dans la formule de d finition par laquelle on multiplie le degr de libert
67. vira dans la pr diction par DEPL CALCULE Mot cl IMPL EX 9 IMPL EX Pr cise les caract ristiques de la m thode de r solution du probl me incr mental non lin aire m thode IMPL EX Op rande MATRICE Q PREDICTION TANGENTE Q REAC INCR 1 DEFAUT La matrice utilis e pour les it rations globales de la m thode est la matrice tangente R5 03 01 aucune autre n est permise Conform ment la m thode la matrice tangente de pr diction est r valu e tous les incr ments de temps REAC_INCR 1 Mot cl RECH_LINEAIRE Q BECH LINEAIRE F La recherche lin aire peut permettre d am liorer la convergence de la m thode de Newton Cf R5 03 01 pour plus de d tails Attention Il est d conseill d utiliser la recherche lin aire avec les d formations GROT GDEP pour les mod lisations COQUE 3D et en pr sence de contact Op rande METHODE Q METHODE CORDE DEFAUT MIXTE PILOTAGE Permet de choisir la m thode de recherche lin aire c est dire l algorithme de recherche du z ro de la fonctionnelle voir doc R5 03 01 La m thode CORDE par d faut est la m thode la plus simple c est une m thode s cante unidimensionnelle Manuel d utilisation Fascicule u4 51 Analyse statique Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster W Titre Op rat

Code_Aster

Contents

Download Pdf Manuals

Related Search

Related Contents