Home

[U4.51.03] Opérateur STAT_NON_LINE

Contents

1. OUT lorma lma VALE MIN MAX D ELGA SIEF ELGA VARI DEFAUT DEFAUT DEFAUT DI p FAUT DI FAUT I DEFAUT DEFAUT list k DEFAUT 1 En DEFAUT formule _ aster ELGA DEFAUT 1 gr maille l maille DEFAUT Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version 12 Titre Op rateur S TAT_NON LINE Responsable Mickael ABBAS SUIVI DDL Si EVAL ELGA POINT o SOUS POINT VALI D Si CHAM est de type NOEU Q TOUT GROUP_MA MAILLE NOEUD GROUP_NO 9 9 9 9 S OBSE LIST_INST INST PAS OBS CRITERE SSS S posi CRITERE Q PRECISION CRITERE PRECISION Q TITRI NOM CHAM D Q EVAL CHAM ETAT INIT DEE lgrm lma no rno OUI NON linst LINSTE pas RELATIF ABSOLU spi a RELATIF 1 0E 6 prec ABSOLU prec Po titre DEPL VITI ACC SIEF E VARI D DI LE Ga GA FORC NODA
2. 3 6 5 Op rande PRECISION CRITERE O PRECISION prec Confer U4 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 A 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 NUME 1 2 3 4 5 6 7 INST 0 10000001 10 10000002 10 10000003 0 10000004 0 10000005 0 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 pouvoir s lectionner l instant dans l exemple PRECISION 1E 8 conviendra INST 0 10000004 car l cart relatif entre les instants vaut Manuel d utilisation Fascicule
3. Table des mati res D BU dream eege eege EES deed dE tard 1 2 OVO O a 5 S LOS LE LES LE LOS enine ieia EE Eege 11 3 1 Op randes MODELE CHAM_MATER CARA _ELEM ii nnnnnnnmmnnnnmns 1i S 2 MOLCI E XG HE 11 Spee TYPE CHARGE eebe 12 3 9 Oe CR EE 13 SE MOECI SOUS E 13 KE ae ee e E 13 3 42 Op randes TOUT SUPER MAILLE LLL Sua un asas iaaii 13 34 3 Op rahde FONC AT EE 13 5 Mot cl COMPORTEMEN T EE 13 LO ee TT EE 14 3 6 1 Op randes SIGM VARI 7 DEPL STRX isiin diiin iaaii nine Ee tulle 14 26 2 Op randes EVOL el EE 15 3 6 5 Op rand PRECISION CRITERE uuu 2h anus unusasussaasssbaqaaaqkaswaqaqhanapasmapuqakaq atuqhina 16 3 1 ege e TC ME 17 S F T OP rANTE HST INST nana ne den sua eue nee ue tan 17 3 7 2 Op randes NUME_INST_INIT INST_INIT NUME_INST FIN INST FIN 17 5 1 0 Op rande PRECISION E 18 38 Operande GIRI Ee TI EE 19 S S Operalide METHODE uy uuu nes en nn iaiia shua a aiaia R 19 SNO Motiele E RE 19 3 0 T Operanmde PREDICTION DEE 19 10 2 Op rande MATRICE E 20 3 10 3 Operande EVOL NOLDI eet ee 21 3M Motele RECH IN eege eege Ee EE hayaya ushaq kus 21 911 1Op rand e ME THODE ege 21 3 11 2 Op rande RESI UNE _RELA ITER_LINE MAXI rssss
4. ITER GLOB ELAS 125 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 19 pour r soudre le probl me global chaque instant 25 par d faut On rappelle 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 rande ARRET Q ARRET OUI DEFAUT Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e A ster Version 12 Titre Op rateur STAT_NON_LINE Date 12 10 2015 Page 31 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 3 15 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 es
5. SELECTION Q SOLVEUR Q CONVERGENCE RESI GLOB RELA LA F voir le document D I RESI GLOB MAXI RESI GLOB REL RESI COMP REL RESI REFE REL SIGM REFE EPSI REFE DEPL REFE FORC_REFE VARI REFF FLUX THER REFE FLUX HYD1 REFE FLUX HYD2 REFE ITER_GLOB_ELAS ITER_GLOB_MAXI ARRET O CRIT _STAB TYP es ND _FREQ COEF DIM ESPACE L RIGI GEOM MODI RIGI CHAR CRIT DDL EXCLUS DDL STAB LIST _INST INST PAS CALC PRECISION DS E cmult etarmax etarmin etamax etamin OUI NON SANS DECROISSANT CROISSANT N A R ORM INCR DEPL NGL INCR DEPL ESIDU LA E 6 resmax resrel rescmp resref sigref sigre depre forre varre sigre sigre sigre 25 maxelas 10 maglob OUI NON syyn a x Eh Hy Eh Eh Eh Eh Eh T Ei FLAMBEMENT STABILITE 3 nbfreq listl ddl 1 E 6 prec U4 50 011 Date 12 10 2015 Page 7 40 Cl U4 51 03 R vision 13949 EFAUT O 0 0 J tJ FAUT J EFAUT DEFAUT DEFAUT R EFAUT EFAUT 4 U J Z J J gt 9 EFAUT
6. Le mot cl COEF_DIM_ESPACE 5 est sa valeur par d faut permet l utilisateur de contr ler la taille du sous espace dans la m thode de Sorensen la taille du sous espace est gale la multiplication de ce coefficient par la valeur nbfrea renseign e pr c demment L int r t tant de pouvoir r duire cet espace dans le cas o l on utilise en plus l op rande DDL STAB 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 Op rande MODI_RIGI MODI RIGI NON DEFAUT Z OUT Le mot cl MODI RIGI NON par d faut permet de pr ciser si la matrice de rigidit globale et la matrice de rigidit g om trique si elle est utilis e doit tre modifi e au niveau des degr s de libert que l on liste avec DDL EXCLUS Cela permet par exemple pour des mod les mixtes de ne mener l analyse de stabilit en excluant certains type de degr de libert et en corrigeant les matrices de rigidit globales pour que les termes li s ces degr s de libert ne viennent pas perturber la recherche d instabilit On donne d autres d tails dans le paragraphe consacr DDL EX
7. L autre utilisation possible est celle de efamax en tant que charge limite maximale Attention Avec la loi ENDO ISOT BETON ces deux mots cl s sont obligatoires car ils sont utilis s pour fixer les bornes de pilotage au niveau l mentaire 3 12 100p rande PROJ BORNES PROJ BORNES OUl DEFAUT 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 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
8. 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 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 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 Sion r alise un calcul en axisym trique et que l on impose des forces nodales ces efforts doivent tre divis s par 27 on travaille sur un secteur d un radian par rapport aux chargements r els De m me si l on souhaite calculer la r sultante des efforts le r sultat est multiplier par Jm pour avoir la r sultante totale sur la structure compl te De m me en contraintes planes ou en d formation plan
9. ud soit gal 2 soit d apr s la liste d instants des pas de 0 2 soit un coefficient cmult 1 0 2 5 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 EI TYPE DDL IMPO NOEUD NO1 NOM CMP DX COEF MULT 5 PILOTAGI D 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 pne 5 nX P Ici PP 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 calcul 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 Q ETA PILO R MAX etarmax R Q ETA PILO R MIN etarmin R Ces deux mots cl s permettent de d finir l intervalle de recherche des valeurs de pilotage 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 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp
10. 2015 Page 16 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 LIST2 DEFI LIST REEL DEBUT 20 INTERVALLE F JUSQU A 30 NOMBRE 10 euse U STAT NON LINE r INCREMENT F LIST INST LIST2 I ETAT INIT F EVOL NOLI U NST ETAT INIT 20 Premier STAT NON LINE effectue le calcul des instants 1 10s Second STAT NON LINE effectue le calcul des instants 21 30s l tat initial correspondant l instant 10s du premier STAT NON LINE par d faut INST 10 Cet tat initial correspond pour ce second STAT NON LINE l instant f 20s INST 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 LISTI U2 Il Wei TAT NON LINE INCREMENT F LIST INST LISTI ETAT INIT F EVOL NOLI Ul INST ETAT INIT 0 Premier STAT NON LINE effectue le calcul des instants 1 10s Second STAT NON LINE effectue le calcul des instants 1 10s l tat initial correspondant l instant 10s du premier STAT NON LINE par d faut INST 10 Cet tat initial correspond pour ce second STAT NON LINE l instant f 0s INST ETAT INIT 0
11. ETAT INIT avec DEPL SIGM VARI il n y a pas d instant associ 2 Lorsque l tat est donn 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 ETAT INIT F EVOL NOLI U Premier STAT NON LINE effectue le calcul pour les instants 1 2 3 et 4s Second STAT NON LINE effectue le calcul pour les instants 5 6 7 8 9 et 10s l tat initial correspondant au temps 4s DI Di B Exemple pour montrer l int r t de INST_ETAT_INIT deux listes d instants diff rentes LISTI DREI LIST REEL DEBUT 0 INTERVALLE F JUSQU A 10 NOMBRE 10 U STAT NON LINE INCREMENT F LIST_INST LISTI Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Co d e A ster Version 12 Titre Op rateur STAT_NON_LINE Date 12 10
12. O HO HO o Ia rj D C 1 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version 12 Titre Op rateur S TAT_NON LINE Responsable Mickael ABBAS 9 CRITERE SIGNI ESI PREC INSTAB O ENERGIE 0 ARCHIVAG LIST_INST INST PAS ARCH Q PRECISION L 9 CRITERE Q CHAM EXCLU Q AFFICHAGE Q UNITI PAS Q INFO RESIDU Ka Q INFO TEMPS Q OBSERVATION Q TITRE NOM CHAM Q EVAL CHAM NOM CMP Q EVAL CMP Si EVAL CMP FORMULE FORMULE y SE SC S T 7 F 1 1 n Z Re 2 of d D I d Si CHAM est de type 9 TOUT Q GROUP Q MAILLE Q EVAL ELGA A Date 12 10 2015 Page 8 40 Cl U4 51 03 R vision 13949 RELATIF ABSOLU POSITIF NEGATIF POSITIF NEGATIF 1 E 6 prec instab E ist_r8 C8 pas 1 E 6 prec RELATIF ABSOLU E unite unite NON OUT NON OUT itre DEPL PANE ACCE DEPL ABSOLU VITE ABSOLU ACCE ABSOLU SIEF ELGA VARI ELGA FORC NODA VALE CONT VALE es L
13. SIE F ELGA GROUP MA TOTO EVAL CHAM MAX NOM CMP SIXX SIYY SIZZ EVAL CMP FORMULE FORMULE trace EVAL ELGA MAX Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Version 12 Titre Op rateur STAT NON LINE Responsable Mickael ABBAS Date 12 10 2015 Page 38 40 Cl U4 51 03 R vision 13949 3 19 6 Observation d un champ NOEU 9 EVAL CMP 9 FORMULE VALE DEFAUT f FORMULE form formule aster On commence par choisir les composantes ou la formule entre les composantes SiEVAL CMP VAL E on extrait simplement la liste des composantes donn e par NOM CMP SiEVAL 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 9 EVAL CHAM VALE DEFAUT MIN MAX Jr MOY t 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 MAX le minimum EVAL CHAM MIN ou la moyenne EVAL CHAM MOY Exemple Extraire le
14. TYP g m GROUP MA S z D H pa E t FISSUR HOM CMP DIRE PILO al NON OUI EWTON PLEX EWTON KRYLOV TANGENTE ELASTIQUE EXTRAPOLE DEPL CALCULE evol noli TANGENTE ELASTIQUE 1 mf 0 TE 0 VUE 0 pasmini DC CORDI MIXTI PILOTAG 1 E 1 reslin 3 itelin 1 E 2 rmin 1 E 1 rmax 9 E 3 rexc E ea Ia DDL IMPO SAUT IMPO LONG ARC SAUT_ LONG ARC ANA LIM DEFORMATION PRED ELAS OUI lgrma lma no rno fiss nomcmp direpilo Date 12 10 2015 Page 6 40 Cl U4 51 03 R vision 13949 DEFAUT DEFAUT DEFAUT evol _ noli DEFAUT G lle HO 4H O0 ele m Hj gt C J HOW WOHOXTE DEFAUT 1 gr maille l maille noeud gr noeud sd fiss xfem Kn Kn Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version 12 Titre Op rateur S TAT_NON LINE Responsable Mickael ABBAS COEF MULT ETA PILO TA PILO TA PILO MA ETA PILO MIN PROJ BORNES R MAX R MIN x Eg EJ EVOL PARA
15. cas o l on peut d douaner a priori une partie du domaine Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Co d e A ster Version 12 Titre Op rateur STAT NON LINE Date 12 10 2015 Page 34 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 d instabilit alors il est pertinent de modifier le crit re avec le mot cl SIGNE On rappelle que la charge critique calcul e par CRIT STAB dans le cas o la matrice de rigidit g om trique est prise en compte est l inverse du coefficient multiplicateur du chargement impos qui rend le probl me instable Donc si la valeur calcul e par CRIT STAB vaut 1 cela signifie qu on est instable pour la charge impos e Si on obtient la valeur 1 alors l instabilit se produira pour une charge impos e de m me valeur mais de signe oppos Donc pour des chargements impos s connus et voluant de fa on monotone il est ais de restreindre le domaine d instabilit car on sait que le chargement ne peut changer de signe En revanche pour des chargements cycliques ou quelconques il est plus s r de ne pas restreindre le domaine d instabilit 3 15 110p rande PREC_INSTAB 3 16 3 17 3 17 1 3 17 2 Q PREC INSTAB 1 E 6 DEFAUT prec_instab R Le mot cl PREC_IN STAB permet de d finir la tol rance relative ave
16. d une quantit On v rifie la relation t mult cmult max S pred e as 74t 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 FISS EXP non locale uniquement 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 Ia UJ mult ENDO FRAGILE crit re n est pas intuitive et varie selon les lois de comportement Pour les lois une version ENDO SCALAIRE ENDO FISS_EXP et ENDO ISOT BETON diff rente de la d finition de est utilis e o ce param tre est li l incr ment cmult d endommagement voir R7 01 04 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 diff rent de
17. de Newton Raphson 3 10 1 Op rande PREDICTION Q PREDICTION TANGENTE ELASTIQUE EXTRAPOLE 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 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 MECA TANG ELASTIQUE On utilise la matrice lastique option RIGI MECA EXTRAPOLE 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 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Version 12 T
18. et 4s effectue le calcul pour les instants 5 6 7 8 9 et 10s l tat initial correspondant au temps 4s par d faut INST INIT INST_ ETAT INIT INST 4 USQU A 10 NOMBR IST INST SST NOT EIN 4 IST INST LIST NST_INIT 8 ETAT INIT F EVOL NOLI R vision 13949 10 U EVOL NOLI U E_INST INIT quand 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 NUMI E INST zl 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 0 10000005 0 10000006 0 10000007 Si l utilisateur veut s lectionner l instant correspondant NUMI INST 0 10000004 car l cart relatif entre les instants vaut E 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 NU PRI FCISION 1 E 3 4 et 5 le code s arr te alors en erreur fatale Il suffit alors de changer
19. le param tre ECISION pour pouvoir s lectionner l instant dans l exemple PRI E 8 conviendra Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e A ster Version 12 Titre Op rateur STAT_NON_LINE Date 12 10 2015 Page 19 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 3 8 Op rande CRIT_QUALITE 9 ERRE TEMPS THM NON DEFAUT OUI Cet op rande permet d activer le calcul des indicateurs d erreur ERRE_TEMPS_ OM est l indicateur d erreur temporel pour les mod lisations HM instationnaires Voir R4 10 05 3 9 Op rande METHODE Q METHODE NEWTON IMPLEX NEWTON KRYLOV 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 IMPLEX On utilise l algorithme IMPLEX pour r soudre le probl me voir R5 03 81 NEWTON KRYLOV On utilise une version inexacte de l algorithme de Newton Raphson la pr cision des r solutions de syst mes lin aires par une m thode it rative est adapt e au cours de chaque pas de chargement voir R5 03 01 3 10 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
20. le volume des impressions particuli rement en dynamique explicite Il n a pas d influence sur la r actualisation du fichier au format csv mot clef UNITE 3 18 3 Op rande INFO RESIDU INFO RESIDU NON DEFAUT OUI Cet op rande permet d ajouter une colonne pour chaque 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 18 4 Op rande INFO TEMPS INFO TEMPS NON DEFAUT OUI Cet op rande permet d ajouter une colonne qui donner le temps pass dans l it ration de Newton 3 19 Mot cl OBSERVATION Q OBSERVATION F TITRE titre 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 34 o on stocke tous les champs sur tout le mod le ll sert essentiellement des conomies de stockage mais aussi valuer des champs sur Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http
21. liste d instants litps 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 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 L instant final dernier pas calcul est d sign de la m me mani re que l instant initial soit NUME_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 red 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 initialisation A Exemple simple comportement par d faut Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R a
22. mais sera pilot gr ce au mot cl PILOTAGE 221 Les charges pilotables doivent tre issues de l op rateur AFFE CHAR MECA ou 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 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 ou bien affect par FONC MULT Si TYPE CHARGE 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 Actuellement les chargements qui peuvent tre qualifi s de suiveur sont le chargement de pesanteur pour l l ment de CABLE POULIE la pression pour les mod lisations 3D 3D ST 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
23. maximum de la composante DX du d placement sur GROUP NO TOTO 3 19 7 Contenu de la table OBSERVATION F NOM CHAM DEPL GROUP NO TOTO EVAL CHAM MAX NOM CMP DX La table contiendra au maximum 16 colonnes NOM OBSERVATION K80 Nom donn automatiquement ou par le mot clef TITRE TYPE OBJET K16 La table ne contient que des valeurs r elles donc R NOM SD K24 3 NUME REUSE I Indice de r utilisation de la table en cas de REUSE NUME OBSE I Num ro d ordre de l observation INST R Instant de l observation NOM CHAM K16 Nom du champ observ EVAL CHAM K8 Type d valuation du champ HOM CMP K8 Nom de la composante observ e EVAL CMP K8 Type d valuation de la composante NOEUD K8 N ud o se r alise l observation champ aux n uds MAILLE K8 Maille o se r alise l observation champ aux mailles EVAL ELGA K8 Type d valuation des champs aux points de Gauss POINT I Point d int gration o se r alise l observation champs aux mailles Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e A ster Version 12 Titre Op rateur STAT_NON_LINE Date 12 10 2015 Page 39 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 SOUS_POI
24. 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 22 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 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Co d e A ster Version 12 Titre Op rateur STAT_NON_LINE Date 12 10 2015 Page 22 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 Q RESI LINE RELA 1 E 1 DEFAUT reslin 9 ITER LINE MAXI L3 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 L utilisateur ne peut pas mettre plus de 999 it rations de recherche lin aire pour la m thode CORDE Par con
25. tests sont effectu s Remarque Si les conditions limites de Dirichlet sont impos s par AFFE CHAR CINE limination et non bar 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 cmp DX DY DZ DREI 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 nbddl c 1 nbemp max E j i 1 nbddl gt rescmp o Ei est la partie du r sidu J correspondant la composante c et L le vecteur des forces internes au temps n cor
26. u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e A ster Version 12 Titre Op rateur STAT_NON_LINE Date 12 10 2015 Page 17 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 3 7 Mot cl INCREMENT INCREMENT F 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 011 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 INUME_INST_FIN 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
27. 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 tat initial de traction ou de compression On effectue ensuite un second STAT NON LINE en poursuite partir de l tat initial pr c dent avec le mot cl PILOTAGE D al Remarque DEFORMATION et PRED ELAS ne sont pas disponibles pour les l ments de structures 3 12 2 Op randes NOEUD GROUP NO NOEUD GROUP NO no grno 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 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version 12 Titre Op rateur S TAT_NON LINE Date
28. 12 10 2015 Page 25 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 3 12 3 3 12 4 3 12 5 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 n uds 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 n ud renseign correspondant alors l extr mit d une ar te que l on souhaite piloter Signalons alors les r gles suivantes si deux n uds sont les extr mit s respectives de deux ar tes non ind pendantes une seule sera retenue fig 3 12 2 2 si un n ud est extr mit de plusieurs ar tes on retient arbitrairement la premi re rencontr e par l algorithme sl deux n uds sont extr mit s d une m me fissure fig 3 12 2 3 une erreur sera renvoy
29. C MULT fmult COMPORTEMENT F voir le document U4 51 11 o ETAT INIT EI SIG sig cham elem VARI vain cham elem DEPL depl cham no STRX strx cham elem COHE cohe cham elem EVOL NOLI evol evol noli NUME ORDRE nuini I INST instini R PRECISION 1 0E 3 DEFAUT prec R CRITERE RELATIF DEFAUT ABSOLU NUME DIDI nudidi I INST_ ETAT INIT istetaini R INCREMENT Pi LIST INST litpsr8 listr8 litps list inst NUME_INST INIT nuini I INST INIT instini R NUME_INST FIN nufin I INST FIN instfin R PRECISION 1 0E 3 DEFAUT prec R Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version 12 Titre Op rateur S TAT_NON LINE Responsable Mickael ABBAS Q CRIT QUALITE Q ERRE TEMPS THM Q METHODE Q NEWTON PREDICTION EVOL NOLI MATRICE REAC INCR REAC ITER REAC ITER ELAS PAS MINI ELAS Q PECH LINEAIRE METHODE D ELA RESI LINE ITER LINE MAXI RHO MIN RHO MAX RHO EXCL Q PILOTAGI
30. CLUS Si la liste de degr s de libert exclus est vide alors MODI_RIGI ne sert donc rien Op rande CHAR CRIT CHAR CRIT 10 10 DEFAUT intcc 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 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e A ster Version 12 Titre Op rateur STAT_NON_LINE Date 12 10 2015 Page 33 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 des charges critiques dans cet intervalle Ce mot cl ne peut tre utilis que sous la condition RIGI _ GEOM OUI valeur par d faut de l option RIGI GEOM 3 15 8 Op rande DDL EXCLUS 3 15 9 9 DDL EXCLUS DX DY Le mot cl DDL EXCLUS liste vide par d faut d signe l ensemble des degr s de libert que 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 ousiMODI RIGI OUI Dans le cas RIGI GEOM NON et MODI RIGI 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 f
31. Co d e A ster Version 12 Titre Op rateur STAT NON LINE Date 12 10 2015 Page 1 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 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 Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Co d e As ter Version 12 Titre Op rateur STAT NON LINE Date 12 10 2015 Page 2 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949
32. D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version 12 Titre Op rateur S TAT_NON LINE Date 12 10 2015 Page 27 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 3 12 9 Si on ne pr cise pas de valeurs c est o pour etarmin et o 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 On l utilise pour arr ter proprement le calcul lorsque ETA PILOTAGE atteint les bornes de cet intervalle Cet intervalle doit tre plus restrictif que l intervalle de recherche d fini pr c demment car ce dernier est appliqu dans tous les cas Le principe de fonctionnement est le suivant convergence des it rations de Newton si l on a atteint l une des bornes le calcul s arr te Une utilisation possible de cet intervalle est la suivante Dans le cas de pr sence de snap back en fixant efamin une faible valeur cela permet d arr ter le calcul avant une d chirure endommagement complet de l chantillon et viter ainsi la divergence au dernier pas de temps
33. D PLAN HHM D PLAN HM D PLAN THH D PLAN THHM D PLAN THM et la force centrifuge en grands d placements mot cl ROTATION dans AFFE CHAR MECA Important La pression peut tre d fini e par une fonction d pendant de la g om trie AFFE CHAR MECA F Dans le cas suiveur on peut choisir si cette d pendance se fait par rapport la g om trie initiale avec les param tres X Y et Z de DEFI FONCTION ou par rapport la g om trie r actualis e avec les param tres XF YF et ZF Si TYPE CHARGE 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 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e A ster Version 12 Titre Op rateur STAT_NON LINE Date 12 10 2015 Page 13 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 3 3 3 4 3 4 1 3 4 2 3 4 3 3 5 total Par exemple pour un d placement impos mot cl DDL IMPO de l op rateur AFFE CHAR MECA la condition sera de la forme u u d o u est le d placement d fini par NUME_DIDI etnon u d Mot cl CONTACT CONTACT contact Ce m
34. DEPL CALCULE soient issus de calcul utilisant le m me maillage car les conditions limites doivent tre coh rentes Du fait de cette projection sur les conditions limites ces deux options sont incompatibles avec les fonctionnalit s du PILOTAGE Attention du fait de impossibilit de projeter correctement les conditions limites d un maillage l autre il est d sormais fortement d conseill d utiliser EXTRAPOLE et DEPL CALCULE partir d un maillage diff rent du calcul courant Utilit eSi on a obtenu une premi re solution sur le m me maillage avec d autres param tres mat riaux ou un autre comportement les champs de d placements peuvent tre r utilis s dans le calcul 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 ny 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 inform
35. DEPL VIBR MODE _VIBR Mode vibratoire juste disponible dans Champ de DYNA NON LINE type DEPL 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 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 Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html
36. GLOB RELA Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e A ster Version 12 Titre Op rateur STAT_NON LINE Date 12 10 2015 Page 29 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 Toutefois ce m canisme de basculement ne peut pas fonctionner au premier pas de temps En effet pour trouver une valeur de RESI GLOB MAXI raisonnable de mani re automatique puisque l utilisateur ne l a pas renseign on a besoin d avoir eu au moins un pas converg sur un mode RESI GLOB RELA D s lors si le chargement est nul d s le premier instant le calcul s arr te L utilisateur doit d j alors v rifier que le chargement nul est normal du point de vue de la mod lisation qu il r alise et si tel est le cas trouver un autre crit re de convergence RESI GLOB MAXI par exemple Si cet op rande est absent le test est effectu avec la valeur par d faut sauf si RESI GLOB MAXI est pr sent RESI GLOB MAXI resmax R L algorithme continue les it rations globales tant que max F gt resmax Ce D 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 RESI GLOB MAXI sont pr sents tous les deux les deux
37. 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 cmult dE u P 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 S Au emult Au At avec Au 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 fissure 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 emt Il All Ar avec Au BI Gd Au c a 1 o c sont les composantes du d placement M me si le groupe de n ud du pilotage est r duit un seul n ud il faut quand m me utiliser GROUP MO ANA LIM Ce mode de pilotage est sp cifique au calcul de charge limite loi NORTON HOFF par
38. NT I Sous point d int gration o se r alise l observation champs aux mailles VALE R Valeur Le param tre NUME REUSE sert en cas d enrichissement de la structure de donn es r sultat En effet si la reprise crase d anciens num ros d ordre dans la structure de donn es r sultat voir mot clef ETAT INIT ce n est pas le cas des valeurs dans la table d observation qui n est jamais modifi e r troactivement On peut donc avoir deux valeurs diff rentes pour le m me instant dans la table la distinction se fera alors sur NUME REUSE 3 20 Mot cl SUIVI _DDL o SUIVI DDL EI 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 fonctionnalit s activ es Le mot cl facteur 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 instants extraire puisqu 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 TITRE ltitre list k Ce mot cl attend une liste de tr
39. VALE MIN MAX MOY MINI ABS MAXI ABS NOM CMP lnocmp EVAL CMP VALE FORMULE Si EVAL CMP FORMULE 9 FORMULE SS form Si CHAM est de type ELGA SIEF ELGA o TOUT OUl Q GROUP MA lgrma Q MAILLE lma 9 EVAL ELGA VALE MIN MAX U e Si EVAL ELGA VALE POINT pi Q SOUS POINT spi Si CHAM est de type NOEU 9 TOUT OUI Date 12 10 2015 Page 9 40 Cl U4 51 03 VARI ELGA R vision 13949 1 EFAUT list k DEFAUT Il Kn DEFAUT formule _ aster DEFAUT 1 gr maille 1 maille DEFAUT DEFAUT Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version 12 Titre Op rateur STAT NON LINE Responsable Mickael ABBAS 9 9 9 9 9 INFO Q TITRE GROUP_MA MAILLE NOEUD GROUP_NO 1 2 tx lgrma lma rno Date 12 10 2015 Page 10 40 Cl U4 51 03 R vision 13949 1 gr maille 1 maille noeud gr_ noeud DEFAUT Kn Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e A ster Ve
40. 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 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 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Version 12 Titre Op rateur STAT NON LINE Date 12 10 2015 Page 24 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 AE Io Ce mode de pilotage est valable pour toutes les lois de comportement y compris en grandes d formations SIMO_MIEHE cmult max At point de Gauss PRED ELAS PRED ELAS assure qu au moins un point de Gauss de la structure sorte du seuil d lasticit lin aris f pred e A
41. 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 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 Q SELECTION NORM INCR DEPL DEFAUT ANGL_INCR_DEPL RESIDU MIXTE Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Co d e A ster Version 12 Titre Op rateur STAT NON LINE Date 12 10 2015 Page 28 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 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 d
42. ations ne sont pas prises en compte 3 10 2 Op rande MATRICE MATRICE TANGENTE REAC INCR vail DEFAUT mf REAC ITER 0 DEFAUT Jit 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 m positif ou nul 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 1 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e A ster Version 12 Titre Op rateur STAT_NON LINE Date 12 10 2015 Page 21 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 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 vo
43. c laquelle on souhaite v rifier le crit re d instabilit qui est param tr par le mot cl pr c dent SIGNE Mot cl ENERGIE 9 ENERGIE EI Ce mot cl permet d activer le calcul du bilan d nergie son affichage en cours de calcul et son stockage dans la table de nom PARA CALC Le bilan d nergie peut tre extrait de cette table l aide de la commande RECU TABLE U4 71 021 Mot cl ARCHIVAGE Q ARCHIVAGE F 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 Op rande LIST_INST I INST PAS ARCH Q LIST INST list r8 INST GES PAS ARCH npas La d signation des instants stocker est effectu e soit par une liste d instants 1ist r8 ou l 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 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 3 L tat initial est
44. chier par la commande LIRE RESU Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version 12 Titre Op rateur S TAT_NON LINE Date 12 10 2015 Page 15 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 3 6 2 Op randes EVOL NOLI 3 6 3 3 6 4 EVOL NOLI evol Nom du concept de type evol poli d o sera extrait l tat initial Op rande NUME ORDRE INST NUME DIDI NUME ORDRE nuini INST instini Extraction de l tat m canique initial dans evo1 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 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
45. cule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version 12 Titre Op rateur S TAT_NON LINE Date 12 10 2015 Page 37 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 pour des champs aux points de Gauss SIEF ELGA VARI ELGA on extrait la liste des mailles Attention ne pas utiliser TOUT OUT sur des gros maillages 3 19 5 Observation d un champ ELGA Q EVAL CMP VALE DEFAUT FORMULE FORMULE form formule_aster On commencoe par choisir les composantes ou la formule entre les composantes SiEVAL CMP VALE on extrait simplement la liste des composantes donn e par NOM CMP SiEVAL 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 Zt MIN t Zt MAX t POINT pi r Q SOUS_ POINT spi I Une fois valu les composantes ou la formule sur les composantes on peut 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
46. de ITER GLOB MAXI uuu uuu aiaiai kadadin dons frame aana 30 214 Opanda MER Ee 30 AE OPT CEARR EE 31 215 Morde ENEE 31 3 15 1 Op rande LIST INST INST PAS CAICG ksanssssazansanakakanzananaassssakanss 31 3 15 2 Op rande PRECISION CRITERE EE 32 3153 Op rand NB PREO ett 32 3 19 6 Op rand MODI RIG ege 32 3 15 7 Op rand CHAR CRIT EE 32 3 15 8 Op rande DDL EXCLUS uuu uuu aana aa a EAA AA 33 3139 Oprande DDL STAB uu uuu utu asua tn anti ta tata en ant 33 ST Rat E ee EE 33 3 15 11 Op ra de PREC INS RE 34 Sei ee ET EE 34 Bretelle 34 3 17 1 Op rande LIST INST INST PAS _ARCH8 uu uu asss aasasa n wassaasinaausanasas 34 317 2 Op rande Ee Ca EE 34 317 3 Operande CHAM EE 35 S IS Mot CIE Ee LEE 35 318 1 Op rarde RTR NEE 35 3 18 2 Op rande EE 35 3 18 JOpetande INEO RESIDU rr iiianoe eege 35 3 184 Op rande INEO TEMPS riroriro a anaa a iaaa 35 S S Mot Cle OBSERVATION EE 35 3 19 1 Op randes LIST_INST INST PAS_OBSE OBSE_ETAT INIT 36 3 19 2 Op randes PRECISION CRITERE EE 36 3 19 3 Op randes NOM CHAM NOM CNP uuuuuu enee ege iaiaaeaia aa aaie 36 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL
47. e D une fa on g n rale il est conseill que tous les n uds entr s soient du m me c t de la fissure siun n ud ne correspond aucune ar te fig 3 12 2 4 une erreur est renvoy e S lection automatique Figure 3 12 2 1 GROUP MO 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 Op randes TOUT MAILLE GROUP MA TOUT OUI DEFAUT GROUP MA lgrma MAILLE lma On donne les mailles o groupes de mailles servant piloter le calcul n utiliser qu avec 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 n utiliser qu avec DDL IMPO o LONG ARC Op rande DIRE PILO Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version 12 Titre Op rateur STAT NON LINE Date 12 10 2015 Page 26 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 Q DIRE PILO direp
48. e 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 3 On d conseille l utilisation du chargement de type LIAISON SOLIDE en grandes transformations Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Version 12 Titre Op rateur S TAT_NON LINE Date 12 10 2015 Page 12 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 3 2 2 Op rande FONC MULT 0 FONC MULT f f est la fonction du temps multiplicatrice du chargement pr cis la i e occurrence de EXC IT Le chargement et les conditions aux limites pour n occurrences du mot cl facteur EXCIT sont ey IT ch i 1 Pour les conditions de Dirichlet bien entendu seule la valeur impos e est multipli e par f Par d faut f l 3 2 3 Op rande TYPE CHARGE Q TYPE CHARGE FIXE CSTE DEFAUT SUIV DIDI FIXE PILO Par d faut TYPE CHARGE 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 TYPE CHARGE vaut FIXE PILO le chargement est toujours fixe ind pendant de la g om trie
49. e 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 MIXTE permet de s lectionner la valeur de pilotage en s appuyant sur plusieurs strat gies On commence d abord avec la strat gie NORM INCR DEPL ci dessus Si les r sultats de la fonction objectif la norme de l incr ment de d placement sont trop proches on bascule pour cette it ration sur la strat gie RESIDU L encore si les r sidus sont trop proches on revient la strat gie NORM INCR DEPL et on examine si la liste des r sidus RESI GLOB MAXI du pas de temps courant pr sente des cycles Si c est le cas c est la moins bonne solution de NORM INCR DEPL qui est choisie pour cette it ration Sinon on choisit simplement la meilleure des deux m me si elles ne sont pas suffisamment contrast es 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 correct
50. ement les r sultats du pr c dent calcul au risque d obtenir des r sultats faux Une alarme avertit l utilisateur 3 12 120p rande EVOL PARA EVOL PARA SANS DEFAUT DECROISSANT CROISSANT Cet op rande permet d imposer le la croissance ou la d croissance du param tre de pilotage 3 13 Mot cl SOLVEUR Q SOLVEUR F La syntaxe de ce mot cl commun plusieurs commandes est d crite dans le document U4 50 01 3 14 Mot cl CONVERGENCE 9 CONVERGENCE F 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 GLOB MAXI IRESI GLOB RELA resrel R L algorithme continue les it rations globales tant que max F gt resrel max L bddl 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 essaie de passer du crit re de convergence relatif RESI GLOB RELA 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
51. ersions voir U4 51 11 La version lastique lt pure gt est choisie automatiquement par le code lorsqu il n y a pas d tat initial sinon on s lectionne la version incr mentale Or seule cette derni re est efficace si l on a de fortes variations des caract ristiques mat riaux lastiques dans ce cas par rapport aux variables de commande Remarque Si le calcul a du mal converger dans ce cas le plus simple est de cr er artificiellement un tat de contraintes initiales nul pour forcer la s lection du mode incr mental ou d utiliser une matrice purement lastique op rande NEWTON qui permet une convergence certaine au prix d un grand nombre d it rations 3 6 1 Op randes SIGM VARI DEPL STRX SIGM sig VARI vain STRX strx DEPL depl COHE cohe sig estle champ de contraintes aux points de Gauss vain estle champ des variables internes aux points de Gauss et dep1 est le champ des d placements aux n uds pris l tat initial et strx est le champ d efforts et de d placements correspondant aux l ments de structures Le champ cohe est le champ de variables internes coh sives de la loi CZM LIN MIX dans le cas d un calcul XFEM sa structure particuli re justifie l emploi d un mot cl distinct de VART 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 ou bien avoir t lus dans un fi
52. es K8 CHAMPMAT CHAM MATER Champ de mat riau K8 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Version 12 Titre Op rateur S TAT_NON_LINE Date 12 10 2015 Page 40 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 PARM THETA 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 d instabilit avec STAT NON LINE ou DYNA NON LINE ou des modes vibratoires avec DYNA NON LINE seulement 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 STAB avec Chargement critique du mode de flambement R TYPE FLAMBEMENT MODE _FLAMB CRIT_STAB avec Mode de flambement Champ de TYPE FLAMBEMENT type DEPL CHAR STAB CRIT_STAB avec Valeur d instabilit R TYPE STABILITE MODE STAB CRIT STAB avec Mode d instabilit Champ de TYPE STABILITE type DEPL FREQ MODE VIBR Fr quence du mode vibratoire juste disponible R dans DYNA NON LINE
53. http www gnu org copyleft fdl html Co d e As ter Version 12 Titre Op rateur STAT NON LINE Date 12 10 2015 Page 4 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 3 19 4 Op randes TOUT NOEUD GROUP_NOEUD MAILLE GROUP_ MA 36 3 195 Observatiom d un champ ELGA EE 37 3 19 6 Observation d un champ NOEEU E 38 3 I9 TER TEE 38 3 20 Motole SUIVI RTR EE 39 3 21 Contenu de la structure de donn es EVOL_NOLT dd aaanaenennmnnanne 39 S 22 ele EC DEE 40 S 23 Operalide OR re ne Pen Ne 40 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Version 12 Code Aster Titre Op rateur STAT NON LINE Date 12 10 2015 Page 5 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 2 Syntaxe T_NON LINE statnllevol noli STA reuse statnl evol noli MODELE mo modele CHAM MATER chmat cham mater Q CARA ELEM caras cara elem FXCIT E CHARGE Chi char meca FONC MULT fi fonc formule TYPE CHARGE FIXE CSTE DEFAUT FIXE_PILO SUIV DIDI Q CONTACT contact char contact SOUS STRUC PI CAS CHARGE chi char meca TOUT OUI DEFAUT SUPER MAILLE Ima 1 maille FON
54. ilo 3 12 6 3 12 7 3 12 8 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 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 u E Op rande FISSURE 9 FISSURE fiss Nom de la sd fiss xfem 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 Op rande COEF MOL 9 COEF MULT cmult C est la valeur not e 7 dans la formule de d finition par laquelle on multiplie le degr de libert utilis pour le pilotage Par d faut cette valeur vaut 1 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 X valeur de r f rence P sur le groupe de maille A Pour trouver la charge limite P on va piloter le d placement du n ud NOI On veut que le d placement final suivant x de ce n
55. ines sous structures Op rande FONC_MULT Q FONC MULT f Jf 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 Mot cl COMPORTEMENT La syntaxe de ce mot cl commun plusieurs commandes est d crite dans le document U4 51 11 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Co d e A ster Version 12 Titre Op rateur STAT NON LINE Date 12 10 2015 Page 14 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 3 6 Mot cl ETAT _INIT 9 ETAT_INIT Ce mot clef permet de d finir un tat initial de r f rence 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 Remarques 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 e Dans le cas o l on utilise la m thode continue du contact la reprise de calcul prend automatiquement en compte les tats de contact et les seuils de frottement issus du calcul pr c dent Si on utilise un tat initial dont le MODELE est diff ren
56. it dans le cas o le minimum est n gatif et l estimation du crit re de stabilit associ sont stock s dans l objet r sultat sous le nom MODE STAB CHAR STAB C x 3 15 1 Op rande LIST_INST INST PAS CALC 9 LIST INST list cp INST 1 r8 PAS CALC npas Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e A ster Version 12 Titre Op rateur STAT NON LINE Date 12 10 2015 Page 32 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 3 15 2 3 15 3 3 15 4 3 15 5 3 15 6 3 15 7 Les instants pour lesquels on veut faire un calcul de stabilit sont donn s par une liste d instants 1ist_r8 ou 1_r8 ou par 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 Op rande PRECISION CRITERE PRECISION 1 e 6 DEFAUT prec Q CRITERE RELATIF DEFAUT ABSOLU Permet de s lectionner les instants confer U4 71 00 Op rande NB_FREQ Q NB FREQ 7293 xg DEFAUT nbfreq Le mot cl NB FRE0 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 Op rande COEF DIM ESPACE COEF DIM ESPACE 5 DEFAUT coef
57. itre Op rateur S TAT NON LINE Date 12 10 2015 Page 20 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 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 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 15 Le d placement est projet sur l ensemble des champs cin matiquement admissible comme pour la m thode EXTRAPOLE Remarques e 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 l est n cessaire que les d placements utilis s dans EXTRAPOLE et
58. l 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 VMIS ECMI TRAC VISC ISOT TRAC U4 51 11 3 10 3 Op rande EVOL NOLI Q EVOL NOLI evol noli Nom du concept de type evol_noli qui servira dans la pr diction par DEPL CALCULE 3 11 Mot cl RECH_LINEAIRE Q BECH LINEAIRE Ei 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 20 et en pr sence de contact 3 11 1 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 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
59. luent 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 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 Q MATRICE ELASTIQUE La matrice utilis e correspond au calcul lastique elle n est valu e qu une fois l instant initia
60. lus 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 Il 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_STAB La documentation U2 08 04 pr sentent les diff rentes approches pour les analyses de flambement dans Code Aster Par une tude de signe sur la d riv e seconde de l nergie en respectant l accroissement des degr s de libert irr versibles dans le cas d une m canique dissipative TYPE STABILITE Pour traiter ce cas particulier on impose de prendre comme matrice de rigidit g om trique la matrice identit K l On recherche ensuite le minimum de la fonctionnelle quadratique suivante Cie K x o i est la matrice tangente coh rente l instant tudi et le vecteur X X transpos du champ des inconnus nodaux x sous des contraintes de positivit sur les degr s de libert de de nature irr versible Le signe du minimum permet ensuite de conclure sur la stabilit du chargement Dans le cas o celui ci est n gatif la solution est instable Dans le cas contraire la solution obtenue num riquement est stable Le mode obtenu qui est le vecteur minimisant C x dit mode d instabil
61. mp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Version 12 Titre Op rateur STAT NON LINE Responsable Mickael ABBAS Date 12 10 2015 Page 18 40 Cl U4 51 03 F JUSQU A 10 NOMBRE LIST INST L INST FIN 4 LIST DEFI LIST REEL DEBUT 0 INTERVALLE U STAT NON LINE INCREMENT F U STAT NON LINE reuse U INCREMENT F Premier STAT NON LINE Second STAT NON LINE B Exemple pour montrer l int r t de INST_INIT LIST DEFI LIST REEL DEBUT 0 INTERVALLE F U STAT NON LINE INCREMENT F L I U STAT NON LINE reuse U INCREMENT F L I Premier STAT NON LINE Second STAT NON LINE LIST INST L ETAT INIT 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 3 7 3 Op rande PRECISION Q PR ECISION prec Cf U4 71 00 pour la syntaxe d taill e Ce param tre sert rep rer le bon num ro d ordre NUMI l utilisateur renseigne l instant INST_FIN INST_INIT En effet les instants dans STAT NON LINE Pi E_INST FIN NUMI IST IST effectue le calcul pour les instants 1 2 3
62. nts incompressibles des l ments de grille et de membrane Une variable interne varref si l on utilise des lois non locales gradient de d formation Un flux thermique fthref dans un cas THM Deux flux hydriques fh1ref et fh2ref dans un cas HHM Un d placement depref si on utilise des l ments de joint avec un comportement de type CZM Une force forref et un moment momref si on utilise des l ments de structure discrets barres poutres ou c bles Un coefficient de Lagrange lagref pour les formulations mixtes en endommagement On calcule une r f rence de r sidu 77 un vecteur de m me longueur que le vecteur r sidu La convergence sera r alis e si et seulement si Viec 1 nbadl F lt resref FZ 3 14 4 Op rande ITER GLOB MAXI Q 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 Dans le cas de la r solution d un probl me de contact frottement par la formulation CONTINU Newton g n ralis il est souvent n cessaire d augmenter le nombre d it rations de Newton FI en 3 14 5 Op rande ITER_GLOB_ELAS
63. ois chaines au maximum et permet de nommer la colonne du tableau d affichage Les chaines sont tronqu es 16 caract res 3 21 Contenu de la structure de donn 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 e 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 Caract ristiques l mentair
64. ormulations 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 Dans le cas MODI RIGI OUI cela permet de modifier la matrice de rigidit et si besoin est la matrice de rigidit g om trique de mani re mener l analyse de stabilit en ne tenant pas compte des degr s de libert exclus Par exemple on doit utiliser cette option pour les mod les fluide structure coupl s formulation u D confer documentation R4 02 02 qui est utilisable avec DYNA NON LINE mais pas STAT NON LINE pour exclure les degr s de libert fluide car la matrice de rigidit assembl e globale est singuli re pour ces degr s de libert Pour plus de d tails l utilisateur pourra utilement se reporter aux documentations U2 06 11 et U2 08 041 Op rande DDL STAB 9 DDL STAB DAMG Le mot cl DDL _STAB d signe l ensemble des degr s de libert irr versibles dans l tude de stabilit que l on souhaite r aliser avec CRIT STAB Il ne peut tre utilis que sous les conditions TYPE STABILITE et RIGI GEOM NON Cela permet d effectuer une tude de signe sur la d riv e seconde de l nergie au chargement consid r en ne regardant que les perturbations susceptibles d augmenter les degr s de libert d clar s dans DDL STAB Ceci afin de respecter les conditions m caniques d irr versibili
65. ot 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 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 Mot cl SOUS_STRUC Q SOUS STRUC 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 02 et le cas test ssnv193a Op rande CAS CHARGE CAS CHARGE nocas nocas est le nom du cas de charge utiliser Voir op rateur MACH ELEM STAT U4 62 01 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 certa
66. r les champs l instant initial car l instant initial n est pas g rable par la liste d instants mot clefs LIST INST INST 3 19 2 Op randes PRECISION CRITERE S PRECISION prec 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 Ces op randes permettent de d finir le champ post traiter NOM CHAM ainsi que ses composantes donn es par leur nom NOM CMP On ne peut d finir que 20 composantes au maximum par occurrence du mot clef facteur OBSERVATION On peut post traiter les champs sortis par d faut par l op rateur voir 39 et la documentation U4 53 01 pour DYNA NON LINE ou demander la sortie d autres champs pour l instant seulement EPSI ELGA Pour d autres champs il faut faire une demande d volution 3 19 4 Op randes TOUT NOEUD GROUP NOEUD MAILLE GROUP MA TOUT NON DEFAUT OUI 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 VITE ABSOLU ACCE ABSOLU VALE CONT FORC NODA on extrait la liste des n uds Manuel d utilisation Fasci
67. respondant 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 BEER RELA RESI REFE RELA resref R Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Version 12 Titre Op rateur STAT NON LINE Date 12 10 2015 Page 30 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 SIGM REFE sigref R FORC BEER forref momref ER VARI _ REFE varref R EPSI REFE epsref R FLUX THER REFE fthref R FLUX HYD1 REFE fhiref R FLUX HYD2 REFE fh2ref R DEPL REFE depref R LAGR REFE lagref R 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 pour des l me
68. rsion 12 Titre Op rateur STAT_NON LINE Date 12 10 2015 Page 11 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 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 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 EI 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 ch CHARGE ch H est le chargement m canique comportant ventuellement l volution d un champ de temp rature pr cis la me occurrence de EXCIT Remarques 1 Dans un calcul thermo m canique si la temp rature initiale est diff rente de la temp rature de r f
69. sous point d int gration par POINT et SOUS POINT Les sous points d int gration apparaissent pour les l ments de structures poutres plaques coques tuyaux etc Demander d extraire le maximum EVAL ELGA MAX oule 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 MINI_ABS MAXI_ABS MINI_ABS est la Valeur minimale en absolu MINI ABS 1 3 4 12 0 1 0 1 MAXI ABS est la valeur maximale en absolu MAXI ABS 1 3 4 12 0 1 12 Une fois valu les composantes ou la formule sur les composantes ai nsi que le point sous point d extraction on peut Extraire ces valeurs sur toutes les mailles avec EVAL CHAM VALE Demander d extraire le maximum EVAL CHAM MAX le minimum EVAL CHAM MIN ou la moyenne EVAL CHAM MOY Exemple Extraire e maximum de la trace du tenseur des contraintes sur le GROUP MA TOTO trace FORMULE VALE 0 333 SIXX SIYY SIZZ NOM PARA SIXX SIYY SIZZ OBSERVATION F NOM CHAM
70. ssssssssssssssssssa 22 3 11 3 Op randes RHO MIN RHO MAX RHO EXCL uuu uuu aienti 22 S 2 Motele PILOTAGE EE 22 2 12 T Opernaride dE 22 3 12 2 Operandes NOEUD ZGROUP NO uuu r uuu urus aisukumuahasaakaauhbaqas hu qaa qaa maqaaasnba 24 3 12 3 Op randes TOUT MAILLE GROUP MA ere mme taime anal 25 23 12 2 Op rande NOM CMP iiiaio eege 25 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e As ter Version 12 Titre Op rateur STAT NON LINE Date 12 10 2015 Page 3 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 125 Op rande DIRE E Eeer 26 3 126 0p rand e EEN 26 3 12 7 Op rande COEF MULT estate eege 26 3 12 8 Op rande ETA BIO R MAX ETA PILO R MIN nn nnnenennnnnnnnnnnnnnnnn 26 3 12 9 Op rande ETA PILO MAXETA PILO MIN u doe sos tisn tee genie 27 312 10 Op rande PROJ e 27 3 12 11 Op rande SELECTION EE 27 3 12 12 Op rande EVOL Eege 28 Ee EE 28 3 14 Motele CONVERGENCE EE 28 3 14 1 Op rande RESI GLOB_RELA RESI _GLOB_MAXI ni snnnnenennnnnnnnnnnnnnnns 28 3 14 2 Op rande RESI COMP RBELA uuuuu u uuaassssassssskakasstasqsykanaakansskasakanakakakasassquss 29 3 14 3 Op rande RESI REFE RELA onsu aaaea 30 3144 Op ran
71. 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 Op rande PRECISION CRITERE PRECISION 1 e 6 DEFAUT prec CRITERE ES RELATIF DEFAUT ABSOLU Cf U4 71 00 Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e A ster Version 12 Titre Op rateur STAT_NON_LINE Date 12 10 2015 Page 35 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 3 17 3 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 Si ce mot clef n est pas renseign 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 9 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 PAS Q PAS pas Fr quence de r actualisation de l affichage dans le fichier message Cet op rande permet de r duire
72. t 3 15 100p rande SIGNE SIGNE POSITIF NEGATIF DEFAUT POSITIF NEGATIF Le mot cl SIGNE permet de sp cifier quel type de crit re d instabilit sera utilis Ce crit re permettra de d clencher un un arr t propre base sauvegard e du calcul non lin aire en cas d instabilit si l utilisateur le pr cise sous DEFI LIST INST confer documentation U4 34 03 avec la syntaxe suivante ECHEC F EVENEMENT INSTABILITE ACTION ARRET Sans cette d claration sous DEFI LIST INST m me en cas d instabilit d tect e le calcul non lin aire tentera de se poursuivre c est le mode par d faut Pour les analyses de stabilit sans matrice de rigidit g om trique le crit re d instabilit c est qu une charge critique tende vers 0 ou change de signe Dans ce cas le mot cl STGNE ne sert pas En revanche pour les cas o la matrice de rigidit g om trique est utilis e ce mot cl SIGNE est utile Avec la Valeur par d faut SIGNE POSITIF NEGATIF la soluti on sera d clar e instable dans les cas o une charge critique devient comprise entre 1 et 1 Si l utilisateur choisit l option NEGATIF alors le domaine d instabilit sera born par les valeurs 1 et 0 Inversement l option POSITIF d finira les valeurs 0 et 1 comme limites du domaine d instabilit Le choix par d faut est le plus conservatif mais dans certains
73. t du MODELE de calcul renseign dans l op rateur Code Aster proc de automatiquement au passage entre les deux mod les Si le mod le de calcul est inclus dans le mod le donn dans l tat initial les donn es de l tat initial sont simplement recopi es maille maille pour tous les champs e Si le mod le donn dans l tat initial est inclus dans le mod le de calcul Code Aster commence par copier les valeurs de l tat initial pour les mailles communes puis compl te avec la valeur nulle Dans le dernier cas il faut prendre garde au sens que peut avoir une variable interne nulle initiale dans le calcul En grandes d formations l utilisateur souhaitant utiliser le formalisme GDEF LOG avec un champ de contrainte initial ETAT _INIT se reportera au cas test SSNP159B En effet pour imposer un champ de contrainte initial l utilisateur doit donner en entr e le tenseur de contrainte d fini dans l espace logarithmique T et non celui de Cauchy Les composantes de ce dernier tant stock es en tant que variables internes il faut utiliser les op randes VARI et DEPL d crites ci dessous ces champs peuvent par exemple tre obtenus par la commande CREA CHAMP U4 72 04 Si les caract ristiques mat riaux d pendent fortement des variables de commande il convient d tre prudent sur la s lection de l tat initial au risque de ne pas converger En effet les comportements lastiques FLAS existent en deux v
74. t 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 utilisable qu en mode DEPL CALCULE Cette option est utiliser avec pr caution car elle donne des r sultats faux Mot cl CRIT_STAB 9 CRIT STAB EI 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 dans le cas de ph nom ne m canique r versibles TYPE FLAMBEMENT Ce crit re est alors calcul de la fa on suivante la fin d un pas de temps en petites perturbations on r sout det K A K 0 KI est la matrice tangente coh rente cet instant KE est la matrice de rigidit g om trique calcul e partir du champ de contraintes cet instant En pratique le chargement est instable si JA lt 1 en fait _ lt lt 0 On calcule les valeurs propres par la m thode de Sorensen cf CALC MODES U4 52 021 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 det KI A 1 0 car gT contient alors KE 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 p
75. tre pour la m thode MIXTE sur des probl mes avec endommagement plusieurs dizaines d it rations sont souvent efficaces 3 11 3 Op randes RHO MIN RHO MAX RHO EXCL RHO MIN J 1 E 2 DEFAUT rmin R RHO MAX 1 E 1 DEFAUT rmax R RHO EXCL 9 E 3 DEFAUT rexc R Ces mots cl s fixent l intervalle dans lequel on calcule le coefficient RHO de la recherche lin aire sous la forme T rmin rmax rexc rexc R5 03 01 3 12 Mot cl PILOTAGE 9 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 diaire d un n ud ou groupe de n ud sur lequel on peut imposer diff rents modes de pilotage mot cl TYPE Attention s 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 e Le mot cl PILOTAGE est interdit avec le contact sauf dans le cas du contact XFEM sl n est pas possible de faire du PILOTAGE a
76. vec PREDICTION DEPL CALCULE ou PREDICTION EXTRAPOLE voir 19 3 12 1 Op rande TYPE 9 TYPE DL IMPO ONG ARC NA LIM EFORMATION ED ELAS UT IMPO UT LONG ARC C est le type de pilotage effectu Sept modes de pilotage sont disponibles Confer R5 03 80 pour plus de d tails Manuel d utilisation Fascicule u4 51 Analyse statique Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e A ster Version 12 Titre Op rateur STAT_NON_LINE Date 12 10 2015 Page 23 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 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 amplitude 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
77. www gnu org copylett fdl html Co d e A ster Version 12 Titre Op rateur STAT_NON_LINE Date 12 10 2015 Page 36 40 Responsable Mickael ABBAS Cl U4 51 03 R vision 13949 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 au 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 utiliser que gt 99 occurrences du mot clef OBSERVATION au maximum II est possible de nommer une occurrence de l observation colonne NOM OBSERVATION en utilisant le mot clef TITRE S il n est pas utilis la colonne NOM OBSERVATION contient OBSERVATION xx avec xx Variant de 1 99 3 19 1 Op randes LIST _INST INST PAS OBSE OBSE ETAT _INIT LIST INST list r8 INST 1 r8 PAS OBSE npas OBSE ETAT INIT OUI DEFAUT NON 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 34 L op rande OBSE ETAT INIT pr cise si l on doit observe

Download Pdf Manuals

Related Search

Related Contents

Samsung Galaxy K zoom Käyttöopas  81-121802-00 Jupiter 3W Manual(EN)俄罗斯_R(0.1)  Samsung Gear VR Uživatelská přiručka  Formation E-Learning Ennéagramme Manuel du participant  Dear Exhibitor, On behalf of the entire Tissue World Management  AS-i 3.0 EtherNet/IP+Modbus TCP    MultiLink™ - Emerson Network Power  ランチングホイール取扱説明書 - JOYCRAFT ジョイクラフト  

Copyright © All rights reserved. Failed to retrieve file