[U4.34.03] Opérateur DEFI_LIST_INST
Contents
1. Si ACTION DECOUPE REAC PRECOND AUTRE PILTOAGE Q SUBD METHODE MANUEL DEFAUT AUTO Si SUBD METHODE MANUEL D SUBD PAS dus DEFAUT pas R D GUDD NIVEAU 3 DEFAUT niv T SUBD PAS MINI Or DEFAUT pasmin R Si SUBD METHODE AUTO D SUBD PAS MINI EN DEFAUT pasmin R Si ACTION ITER_SUPPL Q PCENT_ITER_PLUS 50 DEFAUT pcent I Q SUBD METHODE MANUEL DEFAUT AUTO Si SUBD METHODE MANUEL D GUDD PAS 4 DEFAUT pas R D GUDD NIVEAU 3 DEFAUT niv I D SUBD PAS MINI oye DEFAUT pasmin R Si SUBD METHODE AUTO D SUBD PAS MINI SS L DEFAUT pasmin R Si EVENEMENT COLLISION ACTION DECOUPE DEFAUT ARRET Si ACTION DECOUPE Q SUBD METHODE FAUTO DEFAUT MANUEL Si SUBD METHODE MANUEL D GUDD PAS 4 DEFAUT pas R D SUBD NIVEAU SO DEFAUT niv I D GUDD PAS MINI D DEFAUT pasmin R Si SUBD METHODE AUTO SUBD INST subd inst R SUBD DUREE SUBD DUREE R Si EVENEMENT INTERPENETRATION PENE MAXI penmax R Q ACTION ADAPT COEF PENA DEFAUT ARRET Si ACTION ADAPT_COEF_PENA Q COEF MAXI 1E12 DEFAUT coema2. Remarque Il n y a pas de valeur par d faut pour VALE I dans le catalogue de l op rateur DEFI LIST INST Cependant si sa valeur n est pas renseign e on prendra pour VALE I la moiti du nombre maximal d it rations de Newton d clar dans STAT NON LINE CONVERGENCE ITER GLOB MAXI 3 3 3 Op rande MODE CALCUL TPLUS MODE CALCUL TPLUS FIXE DEFAUT DELTA GRANDEUR ITER NEWTON IMPLEX Cet op rande sert sp cifier le mode de calcul du nouveau pas de temps en fait du coefficient c Si MODE CALCUL TPLUS FIXE alors le coefficient c est fixe Sa valeur est d termin e gr ce l op rande PCENT_AUGM Le nouveau pas de temps vaudra alors A f eil DE 4 f Si MODE CALCUL TPLUS DELTA GRANDEUR alors le coefficient c est li la variation au cours du pas de temps d une grandeur choisie On utilise la formule suivante E e c min CT avec N l ensemble des n uds ou des points de Gauss pour lesquels V n est pas ieN i nul La valeur de r f rence H est choisie par l utilisateur Si le champ AH est identiquement nul alors on consid re que le crit re l v nement n est pas v rifi Si MODE CALCUL TPLUS ITER NEWTION alors le coefficient c est li la variation au cours du pas de temps du nombre d it rations de Newton On utilise la formule suivante
3. Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Has Titre Op rateur DEFI LIS T_INST Date 26 11 2012 Page 9 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 max Alvi MIEN pg a valref L incr ment A est prendre au sens incr ment entre deux instants D Exemple n 1 avec une loi VMIS ISOT_LINI ECHEC F EVENEMENT DELTA GRANDEUR VALE REF 0 le 2 NOM CHAM VARI_ELGA NOM CMP SS PARK On aura chec si l incr ment de la plasticit cumul e entre deux pas de temps d passe 0 1 Exemple n 2 ECHEC _ F EVENEMENT DELTA GRANDEUR VALE REF 5 e 2 NOM CHAM DEPL NOM CMP DX _ F EVENEMENT DELTA GRANDEUR VALE REF 5 e 2 NOM CHAM DEPL NOM CMP DY On aura chec si la norme norme sup du d placement en 2D varie de plus de 5e 2 entre deux pas de temps Remarque dans les deux exemples pr c dents bien que EVENEMENT ERREUR ne soit pas pr sent il sera quand m me activ en dur 3 2 1 4 v nement COLLISION Cet v nement permet de d couper le pas de temps si le contact est d tect impact C est une fonctionnalit particuli r
4. ADAPTATION Les instants de calculs d finis par l utilisateur voir 3 1 2 sont cependant respect s 3 1 2 Op randes LIST _INST et VALE LIST INST list 1 r8 VALE vale R Je RESULTAT resu resultat SUBD PAS npas I Les instants de calcul peuvent tre d finis de 3 mani res diff rentes e soit en renseignant le mot cl LIST INST avec un concept d fini au pr alable par l op rateur DEFI LIST REEL U4 34 01 e soit en renseignant le mot cl VALE avec une liste de r els e soit en renseignant le mot cl RESULTAT avec un concept produit par STAT NON LINE ou DYNA NON LINE Ce mot cl permet de raffiner une liste d instants d j existante En effet on d finit alors une liste d instants npas fois plus fine que la liste des instants r ellement calcul s dans le STAT NON LINE ou DYNA NON LINE Il est noter que cette liste d instants r ellement calcul s peut tre diff rente de la liste des instants donn e en entr e du STAT_ NON LINE ou DYNA NON LINE par exemple cause du sous d coupage et de la liste des instants archiv s si on n archive pas tous les instants Cette fonctionnalit est utile pour v rifier l ind pendance des r sultats la discr tisation temporelle Par exemple pour obtenir une liste d instants 2 fois plus fine on choisit SUBD_PAS 2 Attention Si la SD r sultat resu provient d un calcul avec reprise il n est pas possible de la r utiliser dans la commande car on ne peut extr
5. d it rations est sup rieur pcent ITER_GLOB_ alors on consid re que l on a chou De m me si l extrapolation ne r ussit pas Dans ces deux cas il est possible d activer la d coupe du pas de temps mot clef SUBD METHODE Action ADAPT_COEF _PENA COEF MAXI 1E12 DEFAUT coemax R Si l v nement INTERPENETRATION se d clenche cette action consiste augmenter automatiquement la valeur du coefficient de p nalisation voir op rateur DEFI CONTACT U4 44 11 de mani re diminuer cette interp n tration Le param tre COEF MAXI permet de sp cifier la valeur maximale prise par le coefficient de p nalisation lorsque le m canisme d adaptation automatique est utilis ceci afin d viter des probl mes de conditionnement de la matrice voir les pr cautions d usage des m thodes p nalis es dans U2 04 04 Actuellement le m canisme ne fait qu augmenter ce coefficient et ne le fait jamais diminuer m me si l interp n tration n appara t plus apr s adaptation il a donc tendance transformer un probl me de contact avec interp n tration en probl me sans contact II convient donc de choisir un coefficient de p nalisation initial plus faible que n cessaire voir op rateur DEFI CONTACT U4 44 11 Action CONTINUE Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www
6. 3 2 1 5 v nement INTERPENETRATION nn ssnsrrersneesnese eus susonessennsee 9 E EEN E KT EE 9 3 2 2 D finition des actions effectuer en cas d chec 9 3 221 ee ARRE EE 10 e EE TEE 10 32 2 3 Action AUTRE PILOTAGE anus sonne nets EE 13 32 2A Acon IMPER EE ee 13 3 2 2 5 ACtIOn ADAPT COEF PENA oaoiiaaiiiai 13 Ee De CH CU 14 3 3 5 Op randes VALE_REF NOM _CHAM NOM CMP remmmnnnnnennnanneenne 16 33 60 Op rande NB ITER NEWTON Eege 16 34 Operande NEE eege 16 Manuel d utilisation Fascicule u4 34 Tables et fonctions 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 DEFI _LIST_INST Date 26 11 2012 Page 3 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 2 Syntaxe deflist list inst DEFI LIST INST 1 mots cl s pour la d finition a priori de la liste d instant DEFI LIST Ei METHODE MANUEL DEFAUT AUTO LIST_INST list 1_r8 VALE vale R e RESULTAT resu
7. N ig DE c D SS o N est le nombre d it rations de Newton effectu es et N pef une valeur de r f rence it choisie par l utilisateur On rappelle que la pr diction est consid r e comme l it ration de Newton n 0 ce qui signifie que si on converge d s la pr diction alors CN Si MODE CALCUL TPLUS IMPLEX alors le coefficient c est li aux incr ments de variables internes d formation plastique cumul e ou endommagement du pas pr c dent voir R5 03 81 Dans ce cas le mot cl ADAPTATION ne doit tre utilis qu une seule fois Le coefficient c sera forc ment compris entre 0 5 et 1 2 De plus par d finition de la m thode TMPLEX il n y aura aucune cause d chec 3 3 4 Op rande PCENT _AUGM 0 PCENT AUGM 100 DEFAUT pcent R Cet op rande sert sp cifier le pourcentage d augmentation du pas de temps si MODE CALCUL TPLUS FIXE Par d faut on augmente de 100 c est dire que l on double la valeur du pas de temps La valeur peut tre n gative et doit tre strictement sup rieure 100 3 3 5 Op randes VALE REF NOM CHAM NOM CMP Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Weder Titre Op rateur DEFI LIS T_INST Date 26 11 2012 Page 17 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 VALE REF valref R NOM CHAM DEP
8. SEUIL DEFAUT TOUT_INST AUCUN Cet op rande permet de sp cifier le crit re de d clenchement pour l adaptation du pas de temps Si EVENEMENT SEUIL alors l v nement n est v rifi que si le seuil est franchi Si EVENEMENT TOUT _INST alors l v nement est v rifi tous les instants ce cas est obligatoire avec la m thode IMPLEX Si EVENEMENT AUCUN alors l v nement n est jamais v rifi 3 3 2 Op randes NB_INCR SEUIL NOM PARA CRIT COMP VALE I o NB INCR SEUIL 2 DEFAUT nbincseuil I 0 NOM PARA NP ITER NEWTON DEFAUT o CRIT COMP LE DEFAUT g LT YGE nan GT 7 0 VALE I nbvalseuil I Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster K ier Titre Op rateur DEFI LS T_INST Date 26 11 2012 Page 16 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 Ces op randes ne sont renseigner qu en cas d v nement de type seuil EVENEMENT SEUIL L v nement est v rifi si on a nbincseuil fois de suite un feu vert un feu vert tant caract ris par la condition NOM PARA CRIT COMP VALE I Par exemple si NB INCR SEUIL 2 NOM PARA NB ITER NEWTON CRIT COMP LE et VALE Ip alors l v nement est v rifi si on a fait 2 fois de suite moins de 5 it rations de Newton
9. gnu org copyleft fal html Code Aster Weder Titre Op rateur DEFI LIS T_INST Date 26 11 2012 Page 15 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 Si l v nement INSTABILITE se d clenche cette action consiste tenter de poursuivre le calcul non lin aire C est la valeur par d faut pour cet v nement qui est le seul pouvoir se combiner avec cette action 3 3 Mot cl ADAPTATION L adaptation du pas de temps consiste calculer un nouveau pas de temps Art partir des informations du pas de temps actuel A f et des pas de temps pr c dents A7 1 A7 A Le but est donc de calculer Af c Atf o c est un coefficient r el Le proc d est le suivant Pour l occurrence k du mot cl ADAPTATION on d finit un v nement et un mode de calcul de AT Si l v nement est v rifi alors on calcule c avec le mode de calcul choisi On obtient un tableau r capitulatif de ce type Occurrence n 1 c 4 2 Occurrence n 2 v nement non v rifi Occurrence n 3 17 Occurrence n 4 c 3 9 Ensuite on choisit c min c parmi les occurrences dont l v nement est v rifi Dans cet exemple le nouveau pas de temps vaudra AT Le Pour utiliser MODE CALCUL TPLUS IMPLEX il faut que ce mot cl ne soit utilis qu une fois 3 3 1 Op rande EVENEMENT EVENEMENT
10. op rande PAS MINI n agit que pour le calcul d un nouveau pas de temps L op rande PAS MAXI permet de sp cifier le pas de temps maximal mais son fonctionnement est diff rent de celui de PAS MINI En effet si au cours de l adaptation du pas de temps on est amen traiter un pas de temps sup rieur pasmax alors le calcul ne s arr te pas contrairement au fonctionnement de PAS MINI mais le pas de temps est alors limit pasmax Si cet op rande n est pas renseign on ne fera pas cette tape s il n est pas renseign et que MODE CALCUL TPLUS IMPLEX alors le PAS MAXI vaut 10 fois le premier pas de calcul forc ment entr par l utilisateur via LIST NSL Si At gt PAS MAXI At PAS MAXI L op rande NB PAS MAXI permet de sp cifier le nombre de pas de temps au del duquel le calcul s arr te La valeur par d faut est un million c est aussi la valeur maximale autoris e 3 2 Mot cl ECHEC Ce mot clef facteur r p table permet de g rer la liste d instants en cas d chec de convergence chaque occurrence du mot cl ECHEC on d finit une cause d chec voir 3 2 1 et une action effectuer si cette cause est satisfaite voir 3 2 2 Les actions peuvent tre diff rentes entre chaque occurrence du mot cl ECHEC Si plusieurs causes d chec sont
11. 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster ch default Titre Op rateur DEFI_LIST_INST Date 26 11 2012 Page 12 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 e WE 1 an Erreur N N lt At gt lt At gt Bue Ton Zoe S Erreur No N N lt At At gt lt At t gt Lu eis 4 d a AG 45 B NA 1 ls r gt a At At gt ME e 3 4 D coupe r cursive en mode MANUEL avec SUBD_PAS 2 Apr s la premi re d coupe l chec se produit dans le nouveau deuxi me pas de temps On De At f f R red coupe celui ci avec A i A n SC Le pas de temps suivant reste inchang A Ces A CG le premier pas de temps est converg et donc reste galement inchang Ate AE Ce qui donne Pas 1 2 Initial 1 At AC D coupe 2 A ti A D A d D coupe 3 A rh A De A d Il ne peut y avoir un cart de plus d un niveau entre deux pas de temps D coupe en mode AUTO Le mode de d coupe automatique utilise deux possibilit s e Si l v nement d clencheur est une collision EVENEMENT COLLISION alors le pas de temps est d coup en ayant pour objectif d avoir un pas de temps de longueur subd inst maintenu pendant la dur e subd duree Si l instant de d coupe est trop pr s de la collision contact trop rasant la d coupe ne sera pas ef
12. Code Aster SE Titre Op rateur DEFI US T_INST Date 26 11 2012 Page 1 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 Op rateur DEFI_LIST_INST 1 But D finir la liste d instants de calcul ainsi que sa gestion pour les algorithmes de r solution it ratifs commandes STAT_NON_LINE et DYNA NON LINE e possibilit s de re d coupage du pas de temps en cas d chec e possibilit s d adaptation du pas de temps Produit une structure de donn es list_inst Manuel d utilisation Fascicule u4 34 Tables et fonctions 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 DEET LST INS Date 26 11 2012 Page 2 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 Table des Mati res 3M1 COp rand TER ei RTE 6 31 2 Op rand s LIST INST 6t EE 6 3 1 3 Op randes PAS_MINI PAS MAXI NB_PAS_MAXIl nnnsnnnnnnnnnnnsnnsnnnnnnnnnnnsnnnnrnnnnnnnrennnnne 6 912 MOI CI ECHEC EE 7 9 2 1 D finition des causes NEE 7 3 2 1 1 v nement ERREUR rrnernrnnenensanenerenanesnensansnesnonerenenrnenennue T 3 2 1 2 v nement DIVE_RESI nssssssssessssssnsnresnsrnssrerenennstreserierentrnnanenientrtranenererernanh 8 3 2 1 3 v nement DELTA GRANDEUR onns 8 3 2 14 v nement e RT e EE 9
13. L SIEF ELGA VARI ELGA NOM CMP cmp K8 Ces op randes servent sp cifier le calcul du coefficient c d augmentation du pas de temps si MODE CALCUL TPLUS DELTA GRANDEUR cf 3 3 3 avec VALE REF la valeur de r f rence V et NOM_CHAM le nom du champ de grandeur consid r 3 3 6 Op rande NB ITER NEWTON REF NB ITER NEWTON REF nitref I Cet op rande sert renseigner le nombre d it rations de Newton de r f rence not N nef cf 3 3 3 dans le cas o MODE CALCUL TPLUS ITER_NEWTON 3 4 Op rande INFO o INFO Pour le moment cet op rande ne sert pas Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html
14. T OUT OUT OUT O NON DIVE_RESI NON OUI NON NON ON O NON DELTA GRANDEUR OU OUI NON NON ON O NON COLLISION OU OUI NON NON ON O NON INTERPENETRATI OU NON NON NON ON OUI NON ON INSTABILITE OU NON NON NON ON O OUI Par exemple dans le cas le plus g n ral EVENEMENT ERREUR il y a cinq possibilit s e D couper le pas de temps ACTION DECOUPE e R actualiser le pr conditionneur du solveur it ratif ACTION REAC PRECOND e Autoriser des it rations de Newton suppl mentaires ACTION ITER_SUPPL e Changer la solution issue du pilotage ACTION AUTRE_PILOTAGE e S arr ter ACTION ARRET 3 2 2 1 Action ARRET Si l on choisit ACTION ARRET le calcul s arr te proprement la base est sauvegard e et sera ensuite r utilisable en poursuite 3 2 2 2 Action DECOUPE Q SUBD METHODE MANUEL AUTO Si ACTION DECOUPE ET SI SUBD METHODE AUTO SUBD INST subd_inst R SUBD DUREE SUBD DUREE R Si SUBD METHODE MANUEL D SUBD PAS 4 DEFAUT pas R D SUBD NIVEAU 3 DEFAUT niv T SUBD PAS MINT D DEFAUT pasmin R Si SUBD METHODE AUTO SUBD PAS MINT Ohi DEFAUT Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu or
15. aire de liste d instants uniques certains instants seront en double 3 1 3 Op randes PAS MINI PAS MAXI NB PAS MAXI 0 PAS MINI pasmin R 0 PAS MAXI pasmax R NB PAS MAXI nbpasmax I Ces op randes ne sont renseigner qu en cas de gestion automatique de la liste des instants METHODE AUTO Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Weder Titre Op rateur DEFI LIS T_INST Date 26 11 2012 Page 7 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 E 1 E 2 E At At a K p Figure 3 1 Adaptation automatique du pas de temps L op rande PAS MINI permet de sp cifier le pas de temps minimal Si au cours de l adaptation du pas de temps on est amen traiter un pas de temps inf rieur pasmin alors le calcul s arr te Si cet op rande n est pas renseign dans le cas g n ral PAS MINI vaut alors 107 17 c est d ailleurs la valeur minimale autoris e s il n est pas renseign et que MODE CALCUL TPLUS IMPLEX alors le PAS MINI vaut 1 1000 me du premier pas de calcul forc ment entr par l utilisateur via LIST INST Si At lt PAS_MINI STOP Attention la v rification de PAS MINI n est pas effectu pendant le sous d coupage ventuel du pas de temps en cas d chec c est SUBD PAS MINI qui agit pendant le sous d coupage L
16. ait appel au pilotage du chargement Lors du pilotage on peut avoir choisir entre deux solutions voir documentation R5 03 80 que l on discrimine souvent l aide d un crit re RESIDU Si l utilisateur s lectionne l action AUTRE PILOTAGE et que la convergence choue avec le choix de solution qui minimise le r sidu on recommence le pas de temps courant apr s l avoir r initialis correctement mais sans subdiviser la premi re it ration o le choix se pr sente on choisira la solution de r sidu maximal lors de cette seconde tentative On subdivise uniquement si les deux tentatives ont chou Cette action sera donc entreprise uniquement si le nombre maximal d it rations de Newton est atteint et s il s agit de la premi re tentative Action ITER SUPPL o PCENT ITER PLUS 50 PCENT MANUEL AUTO J EFAUT H SUBD METHODE J EFAUT L action ITER_SUPPL permet de faire des it rations de Newton au del du maximum donn par ITER GLOB MAXI ou ITER GLOB ELAS Elle vite que le calcul choue pour quelques it rations de Newton Pour cela gr ce l extrapolation lin aire des r sidus voir 3 2 2 2 il est possible d estimer le nombre d it rations de Newton n cessaire la convergence Le nombre d it rations suppl mentaires est limit par le mot cl PCENT_ITER_PLUS Si le nombre estim
17. alseuil I Q MODE CALCUL TPLUS z FIXE DEFAUT DE TA GRANDEUR ITER_NEWTON IMPLEX si MODE CALCUL TPLUS FIXE D DPCENT AUGM 100 DEFAUT pcent y R si MODE CALCUL TPLUS DELTA GRANDEUR VALE REF valref R NOM CHAM J DEPL SIEF ELGA VARI ELGA NOM CMP cmp TXT si MODE CALCUL TPLUS ITER NEWTON NB ITER NEWTON REF nitref I INFO T DEFAUT Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster K ier Titre Op rateur DEFI LIS T_INST Date 26 11 2012 Page 6 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 2 I Op randes 3 1 Mot cl DEFI LIST 3 1 1 Op rande METHODE METHODE MANUEL DEFAUT AUTO Cet op rande permet de choisir le mode de gestion de la liste d instants Par d faut la gestion est manuelle METHODE MANUEL ce qui signifie que l on va parcourir exactement la liste d instants d finie par l utilisateur voir 3 1 2 En cas de sous d coupage de nouveaux instants peuvent tre ins r s En gestion automatique METHODE AUTO le code calcule automatiquement les nouveaux instants de calcul Les modes de calcul des nouveaux pas de temps sont d finis sous le mot cl
18. d coupe Si la convergence est difficile et que l v nement se d clenche plusieurs fois de suite entra nant le d coupage r cursif du pas de temps on peut arr ter le processus de deux mani res en mode MANUEL e Avec le mot clef SUBD NIVEAU le d coupage s arr tera d s que le d coupage sera fait plus de niv fois e Avec le mot clef SUBD PAS MINI le d coupage s arr tera d s que le pas de temps sera inf rieur pasmin En mode AUTO seul le mot clef SUBD PAS MINI est disponible Remarque concernant le calcul de flambage Lors de calcul de flambage lastoplastique il peut arriver que la matrice tangente du syst me soit singuli re au cours des it rations de Newton En re d coupant le pas de temps on peut passer ces points durs Particularit s de l op rande SUBD NIVEAU e L op rande SUBD NIVEAU est une notion qui devrait tre globale toutes les occurrences du mot cl ECHEC Ainsi cet op rande ne devrait tre d fini qu une seule fois Dans la pratique ce n est pas le cas et rien n emp che de d finir une valeur de SUBD NIVEAU diff rente chaque occurrence du mot cl ECHEC Pour se pr munir de ce cas l l algorithme retiendra au final comme valeur globale et unique de SUBD NIVEAU le maximum des SUBD NIVEAU d finis par l utilisateur 3 2 2 3 Action REAC PRECOND L action REAC PRECOND n a de sens que si l on utilise un solveur lin air
19. e de contact est singuli re e La matrice du syst me est singuli re La d tection de singularit s appuie sur le param tre NPREC du mot clef facteur SOLVEUR voir U4 50 01 e Le solveur lin aire it ratif a atteint le nombre maximum d it rations de Krylov fix par NMAX ITER dans le mot clef facteur SOLVEUR voir U4 50 011 e A convergence des crit res physiques ne sont pas satisfaits lors de l int gration de la loi de comportement notion d event driven M me si on ne trouve aucune occurrence de ECHEC EVENEMENT ERREUR le code en rajoute automatiquement une avec les param tres actuels par d faut c est dire l action DECOUPE en mode MANUEL et avec SUBD NIVEAU A4 v nement DIVE_RESI Le choix de l v nement DIVE RESI permet de sp cifier une autre cause possible d chec event driven ai II y aura chec si le r sidu ne diminue pas suffisamment au cours des it rations de Newton Plus pr cis ment il y aura chec si min R i R i 1 gt R i 2 pour i gt 3 o R j est le r sidu REST GLOB MAXT l it ration de Newton j S il y a chec alors le pas de temps est sous d coup L int r t ce cette fonctionnalit est qu il n est pas n cessaire d attendre que l on arrive la fin des it rations de Newton autoris es pour re d couper le pas de temps en cas de non convergence notable On se rend compte de
20. e it ratif GCPC ou PETSC associ au pr conditionneur simple pr cision LDLT SP Dans ce cas elle permet la r actualisation du pr conditionneur apr s un chec de la r solution lin aire au premier chec on r actualise le pr conditionneur puis l on refait le m me pas de temps pas de d coupe si le pas de temps choue nouveau alors on r actualise le pr conditionneur et l on d coupe Remarque Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Kier Titre Op rateur DEFI LIS T_INST Date 26 11 2012 Page 14 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 Cette action permet de passer des calculs qui pr sentent de fortes non lin arit s ponctuelles d un pas de temps un autre mais qui le reste du temps peuvent se contenter d un pr conditionneur constant Elle n est en revanche pas adapt e des calculs pr sentant de fortes non lin arit s au sein d un m me pas de temps comme les calculs de contact Pour ce type de calcul il est en g n ral n cessaire de forcer la r actualisation du pr conditionneur en utilisant le mot cl REAC PRECOND de SOLVEUR cf U4 50 01 3 2 2 4 Action AUTRE PILOTAGE 3 2 2 5 3 2 2 6 3 2 2 7 L action AUTRE PILOTAGE permet de d finir une action dans un calcul qui f
21. ement utile en dynamique En effet dans ce cas une discr tisation temporelle inadapt e va faire appara tre des oscillations intempestives du d placement et de la vitesse au moment du choc Pour minimiser ces artefacts num riques on recommande en g n ral de d couper finement le pas de temps au moment de la collision ce qui n est possible que pour les cas particuliers o cet instant est d tectable a priori par l utilisateur Dans le cas g n ral on peut utiliser cette fonctionnalit pour que la d coupe se fasse automatiquement 3 2 1 5 v nement INTERPENETRATION PENE MAXI penmax R Cet v nement d tecte l interp n tration d un n ud dans une surface lorsque le contact unilat ral de type r gularis p nalisation est activ L utilisateur donne une valeur de la p n tration maximale qu il veut obtenir param tre PENE MAXI et le code adapte alors le coefficient de p nalisation ACTION ADAPT_COEF_PENA ou s arr te ACTION ARRET pour respecter cette contrainte Cette fonctionnalit n est disponible que pour le contact en FORMULATION DISCRETE voir DEFI CONTACT U4 44 11 3 2 1 6 v nement INSTABILITE Cet v nement permet de g rer le d roulement d un calcul non lin aire lorsqu une instabilit est d tect e lors de l appel CRIT_STAB cf U4 51 03 et U4 53 01 Si ce mot cl CRIT_STAB n est pas utilis alors cet v nement ne peut se p
22. ficace il faut alors changer la discr tisation initiale ou le param tre subd_ inst si le contact devient rasant apr s d coupe e Sinon on tente une extrapolation lin aire de la valeur des r sidus Si l on note iter l it ration de Newton et resi la valeur du r sidu extrapol RESI GLOB RELA et ou RESI GLOB MAXI alors on cherche a a et b tels que a iter a ext resi 1 b Voici l algorithme utilis Pour iter deb fin XX log reSi ter ect Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster EE Titre Op rateur DEFI_LIST_INST Date 26 11 2012 Page 13 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 Si iter gt fin 3 poids 2 Sinon poids 1 On calcule xn xn poids SX SX poids xx sy sy poids iter sxx sxx poids xx Syx syx poids xx iter Fin Pour A SXX SY SX SYX a sx sy Syx xn 2 bh ex sxx xn Pour faire cette extrapolation on demande au minimum trois it rations de Newton et l on consid re que les deux derni res it rations ont deux fois plus de poids que les premi res si ce n est pas possible ou si l extrapolation choue pour une autre raison on passe automatiquement en mode de d coupe MANUEL SUBD PAS 4 Contr le de la
23. g copyleft fal html Code Aster Weder Titre Op rateur DEFI LIS T_INST Date 26 11 2012 Page 11 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 pasmin R Cette action consiste refaire le pas de temps en le d coupant Il existe deux m thodes de d coupe du pas de temps le mode AUTO et le mode MANUEL D coupe en mode MANUEL En mode MANUEL le pas de temps est d coup de mani re uniforme Il faut donc pr ciser le nombre de d coupes l aide du mot clef SUBD PAS 1 2 E E m Erreur bu nm 2m a N At La A gt lt AtP AtY gt Figure 3 2 D coupe en mode MANUEL avec SUBD_PAS 2 On voit un exemple sur la figure 3 2 L v nement ERREUR se d clenche durant le pas de temps 1 initialement de longueur At Ce pas est alors d coup en deux pas SUBD PAS 2 de N 1 At longueurs gales telles que AREA Le pas suivant num ro 2 dans la configuration initiale devient alors le pas 3 sa longueur reste inchang e A t A E Pas 1 2 3 iti _ N N Initial 1 At At A N N N D coupe 2 At At At La notion de niveau de d coupe La d coupe du pas de temps est r cursive si v nement ERREUR se d clenche durant un pas d j d coup on red coupe On voit un exemple sur la figure 3 3 Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright
24. la non convergence plus vite un gain de temps est donc esp r Remarque L activation de cette fonctionnalit avec le contact discret n est pas conseill e car dans ce type de calcul le r sidu ne d cro t en g n ral pas de fa on monotone v nement DELTA GRANDEUR Le choix de l v nement DELTA GRANDEUR permet de sp cifier une autre cause possible d chec event driven ai Il y aura chec si l incr ment d une composante d un champ d passe un seuil fix convergence Il faut alors d finir le nom du champ par l op rande NOM CHAM champ de d placement de contraintes ou de variables internes le nom de la composante par l op rande NOM CMP et la valeur du seuil de r f rence par l op rande VALE REF Si EVENEMENT DELTA GRANDEUR VALE REF valref R NOM CHAM DEPL SIEF_ELGA VARI_ELGA cmp TXT NOM CMP Plus pr cis ment soient v la composante NOM CMP du champ NOM CHAM et valref la valeur de r f rence VALE REF Si le champ NOM CHAM est un champ aux n uds et V est l ensemble de tous les n uds du maillage alors il y aura chec si convergence max Alv VieN gt valref Si le champ NOM CHAM est un champ aux points de Gauss Re est l ensemble de tous les sous points des tous les points de Gauss de tous les l ments du mod le alors il y aura chec si convergence
25. resultat SUBD PAS npas I si METHODE AUTO 0 PAS MINI pasmin R 0 PAS MAXI pasmax R NB PAS MAXI 1000000 DEFAUT nbpasmax I 2 mots cl s pour le comportement en cas d chec ECHEC Fi EVENEMENT ERREUR DEFAUT DIVE_RESI DELTA GRANDEUR Z COLLISION INTERPENETRATION INSTABILITE Si EVENEMENT DELTA GRANDEUR VALE REF valref R NOM CHAM DEPL SIEF_ELGA VARI_ELGA NOM CMP cmp TXT o ACTION DECOUPE DEFAUT CARRET Si ACTION DECOUPE SUBD METHODE MANUEL DEFAUT AUTO Si SUBD METHODE MANUEL 0 SUBD PAS SA DEFAUT pas R SUBD NIVEAU 30 DEFAUT niv I o SUBD PAS MINI 0 DEFAUT pasmin R Si SUBD METHODE AUTO SUBD PAS MINI 0 DEFAUT pasmin R Si EVENEMENT ERREUR ACTION DECOUPE DEFAUT R D ARRET R EAC PRECOND Manuel d utilisation Fascicule u4 34 Tables et fonctions 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 DEFI L IST_INST Responsable Samuel GENIAUT ITER_SUPPL f AUTRE_PILOTAGE Date 26 11 2012 Page 4 17 Cl U4 34 03 R vision 10133
26. roduire Si une instabilit est d tect e une action sera alors d clench e suivant deux possibilit s Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster default Titre Op rateur DEFI_LIST_INST Date 26 11 2012 Page 10 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 ACTION CONTINU D c est la valeur par d faut et on va continuer le calcul m me en cas d instabilit eACTION ARRET si une instabilit est capt e alors on va arr ter la r solution non lin aire proprement la base sera sauvegard e et r utilisable en poursuite Plus de d tails sur le crit re d instabilit et son utilisation sont donn s dans les documentations R7 05 01 et U2 08 04 3 2 2 D finition des actions effectuer en cas d chec Oo ACTION DECOUPE DEFAUT RRET EAC PRECOND ER SUPPL RE PILOTAGE PT_COEF_PENA ONTINUE T U D DH Selon l v nement seules certaines actions sont possibles voir tableau ci dessous ACTION gt ARRET DECOUPE REAC_PRECOND ITER_SUPPL AUTRE_PILOTAGE ADAPT_COEF_PENA CONTINUE EVENEMENT ERREUR OU OU
27. satisfaites simultan ment l ordre de d clenchement de ces checs est pr d termin ERREUR est prioritaire sur DELTA GRANDEUR 3 2 1 D finition des causes d chec EVENEMENT ERREUR DEFAUT DIVE RESI DELTA GRANDEUR COLLISION INTERPENETRATION INSTABILITE L op rande EVENEMENT permet de d finir les causes d chec 3 2 1 1 v nement ERREUR Par d faut la cause d chec est ERREUR Cette cause regroupe les raisons suivantes Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster K ier Titre Op rateur DEFI LIS T_INST Date 26 11 2012 Page 8 17 Responsable Samuel GENIAUT Cl U4 34 03 R vision 10133 3 2 1 2 3 2 1 3 e D passement du nombre maximal autoris d it rations de Newton e chec lors de l int gration locale de la loi de comportement Les param tres locaux de l algorithme utilis pour l int gration de la loi de comportement se trouvent dans U4 51 11 e chec de l algorithme de Deborst pour les contraintes planes ou les mod les 1D On trouvera plus de d tails dans U4 51 11 e chec lors de la r solution de l quation de pilotage voir U4 51 03 e chec lors de la r solution du probl me de contact discret voir R5 03 50 II y a chec lorsque le nombre maximum d it rations de contact est d pass ou lorsque la matric
28. x R Manuel d utilisation Fascicule u4 34 Tables et fonctions Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Version default Titre Op rateur DEFI US T_INST Responsable Samuel GENIAUT Date 26 11 2012 Page 5 17 Cl U4 34 03 R vision 10133 Si EVENEMENT DIVE RESI R D ACTION DECOUPE DEFAUT D SUBD METHODE MANUEL DEFAUT AUTO Si SUBD METHODE MANUEL Q SUBD PAS 4 DEFAUT pas R Q SUBD NIVEAU S DEFAUT niv T Q SUBD PAS MINI AEN DEFAUT pasmin R Si SUBD METHODE AUTO Q SUBD PAS MINI E ay DEFAUT pasmin R Si EVENEMENT INSTABILITE D ACTION CONTINUE DEFAUT ARRET 3 mots cl s pour l adaptation si METHODE AUTO Q ADAPTATION EN D EVENEMENT SEUIL DEFAUT TOUT_INST AUCUN si EVENEMENT SEUIL 0 NB INCR SEUIL GE DEFAUT nbincseuil I D NOM PARA NB _ITER NEWTON DEFAUT CRIT COM E LE DEFAUT LT GE ZG VALE I nbv
Download Pdf Manuals
Related Search