Opérateur CALC_ELEM
Contents
1. 2 1 2 Op randes MODELE CHAM MATER CARA ELEM MODELE mo Nom du mod le sur lequel sont calcul s les efforts les contraintes les d formations Il est optionnel car peut tre extrait du r sultat CHAM MATER chmater Champ de mat riau associ au mod le mo Ce mot cl est optionnel et ne doit tre fourni que dans des cas exceptionnels modification volontaire du mat riau par exemple 0 CARA ELEM carac Caract ristiques l mentaires associ es au mod le mo s il contient des l ments de structure ou si les l ments iso param triques sont affect s par un rep re local d anisotropie Ce mot cl est optionnel 2 1 3 Mot cl EXCIT Ce mot cl facteur optionnel permet de sp cifier les chargements thermiques ou m caniques utiliser pour le calcul des options en lieu et place de ceux qui ont servi dans le calcul de la structure de donn es sp cifi e sous le mot cl RESULTAT La d finition de ce mot cl est identique celle des commandes qui ont construit la structure de donn es resu voir les commandes MECA STATIQUE U4 51 01 STAT NON LINE U4 51 03 DYNA LINE HARM U4 53 11 et DYNA LINE TRAN U4 53 02 2 1 4 Mot cl SOLVEUR La syntaxe de ce mot cl commun plusieurs commandes est d crite dans le document U4 50 01 Remarque dans la commande le solveur n est utilis que pour l estimateur d e2. D_END TEMP_MAX GAMP PCR SEU IND PCO eo SE ANG_ Dar X22 x33 X12 xig est DIST DEV U variation de la porosit du mat riau variation de la masse volumique du mat riau variation de la pression de vapeur saturation du liquide d formation plastique volumique cumul e Indicateur d tat m canique Valeur de l endommagement Indicateur d endommagement Temp rature maximale D formation d viatoire plastique cumul e Pression critique Seuil hydrique Indicateur d irr versibilit hydrique Pression de coh sion Comportement de la roche Seuil isotrope Angle du seuil d viatoire Composantes du tenseur d crouissage cin matique Composantes du tenseur d crouissage cin matique Composantes du tenseur d crouissage cin matique Composantes du tenseur d crouissage cin matique Composantes du tenseur d crouissage cin matique Composantes du tenseur d crouissage cin matique Distance normalis e au seuil d viatoire Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster an Titre Op rateur CALC_ELEM Date 20 06 2011 Page 17 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 DEV SUR CRIT Rapport entre le seuil d viatoire et le seuil d viatorique critique DIST 1SO Distance normalis e au
3. ERRE ELGA ELEM Estimateur d erreur en r sidu calcul aux points de Gauss R4 10 02 QIZ1 ELEM respectivement QIZ2 ELEM Estimateur d erreur en quantit s d int r t bas sur la m thode de Zhu Zienkiewicz lasticit lin aire 2D QIRE ELEM Estimateur d erreur en quantit s d int r t bas sur les r sidus en m canique calcul par l ment Le domaine d utilisation des options QIZ1 ELEM et QIZ2_ELEM est le m me que pour les options ERZ1 ELEM et ERZ2 ELEM et celui de l option OIRE ELEM est le m me que celui de l option ERME ELEM en m canique Il est n cessaire de d finir en plus du probl me initial probl me primal un second probl me probl me dual Ce probl me d finit de mani re sous jacente la quantit d int r t sur laquelle on veut obtenir une erreur A ce jour seulement deux quantit s d int r t sont disponibles e Moyenne d une composante du d placement e Moyenne d une composante du tenseur des contraintes Le probl me dual diff re du probl me primal uniquement par son chargement celui ci tant la quantit d int r t les conditions de bords restant les m mes Ainsi le chargement imposer sur le sous domaine voulu par le biais de la commande AFFE CHAR MECA est FORCE INTERNE effort unitaire
4. Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Pat Titre Op rateur CALC_ELEM Date 20 06 2011 Page 25 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 comprise strictement entre O et 1 un message fatal est mis si cette condition n est pas v rifi e Le p rim tre d utilisation est le m me mais plus r duit que celui de l indicateur d erreur choisi savoir e Pour l indicateur en r sidu l ments finis des milieux continus en 2D triangles et quadrangles ou 3D uniquement les t tra dres pour un comportement lastoplastique e Pour l indicateur de Zhu Zienkiewicz l ments finis des milieux continus en 2D triangles et quadrangles pour un comportement lastique En toute rigueur le calcul de l ordre de la singularit est obtenu partir de l nergie th orique en pointe de fissure quation valable uniquement en lasticit L utilisation de cette option en lastoplasticit est donc manipuler avec prudence SING ELNO D tection des singularit s et carte de modification de tailles aux n uds par l ments Le calcul pr alable de SING ELEM est donc n cessaire Si SING ELEM est absent un message d alarme est mis et l option SING ELNO n est pas calcul e 2 5 7 Autres options VNOR ELEM DE
5. EXTR Date 20 06 2011 Page 16 28 Cl U4 81 01 R vision 6579 ELGA Extraction des variables internes en THM uniquement respectivement aux n uds par l ments et aux points de Gauss Pour pouvoir post traiter les variables internes en THM de fa on plus conviviale des champs ont t cr s Le principe VARI cl plus parlant que V1 V2 NOM _VARI En tant que post traitement ces champs sont calcul s par CALC E nom vari Le nom des nouveaux champs est EXTR NOEU pour le cham no de ces champs est d extraire du champ VARI ELGA ou ELNO pour le cham_elem aux n uds la variable interne qui nous int resse via un mot qu TXM ELGA et EXTR EXTR_ELNO pour les cham_elem et EM et CALC_NO La syntaxe utiliser est la suivante pour un cham_elem GAMP CALC_EL EM R ESULTAT U1 OPTION EXTR E NOM VARI GAMP pour un cham_no NO nouveau mot cl pour indiquer quelle variable on souhaite extraire via un nom cod GAMP CALC_NO reuse GAMP RESULTAT GAMP OPTION EXTR_NOEU Puisqu il s agit juste d extraire une et une seule variable interne les cham elem correspondants doivent avoir t calcul s au pr alable La liste des diff rents noms symboliques des variables internes est DPORO DRHOLQ DPVP SATLIQ EVP IND D IN ETA
6. DEFAUT MOY YINE SUP Calcul des efforts g n ralis s par l ment aux n uds partir des d placements lasticit lin aire cette option n a de sens que pour un mod le avec des l ments de structure poutre coque Dans le cas des mod lisations de plaques avec excentrement DKT DST Q4G GRILLE PLAN permet de d finir le plan de calcul lt MAIL plan du maillage MOY plan moyen e CINE plan sup rieur dans le sens de la normale e SUP plan inf rieur dans le sens de la normale SIPO ELNO Contraintes dans la section de poutre d compos e en contributions de chaque effort g n ralis SN D due l effort normal MY z SMFY O 1 due au moment de flexion MY 5 MZ y SMFZ Tamy due au moment MZ Wya sai SVY T F due l effort tranchant Vy a coefficient de cisaillement dans la direction y Vza gt ie SVZ Dg F due l effort tranchant Vz a coefficient de cisaillement dans la direction z MXR SMT o due au moment de torsion MX YZ x Tout ceci en rep re local rep re principal d inertie de la section droite R3 08 01 Les valeurs de dues aux deux moments de flexion sont les valeurs maximum de celles calcul es en Ymin Ymax d une part eten Zmin Zmax d autre part pour une section g n rale cf AFFE CARA ELEM U4 42 01 Manuel d utilisation Fascicule u4 81
7. ETOT _E EM ENEL ELGA ENEL ELNO options de calcul de crit res cf 2 5 5 TYPI E OPTION CRIT OPTION E SIE Q ELNO Q ELGA Q ELNO EPEQ ELGA Q Q T ELNO ELGA ELNO ENDO ELGA ENDO ELNO B ELNO B ELGA INDL ELGA S go w 7 39 w pa CI ITQ ELNO PTQ ELNO options de calcul d indicateurs d erreur cf 2 5 6 TRI E OPTION INDI ERRE OPTION ERME ELEM ERZ1 ELEM ERZ2 ELEM QIZ1 ELEM OIZ2 ELEM SING ELEM PREC ERR err R SING ELNO ERME ELNO RESU DUAL rd evol_noli QIRE_ELNO DERA ELGA DERA ELNO autres options cf 2 5 7 TYPE OPTION AUTRES OPTION E l SPMX_ELGA NOM CHAM ch cham elem NOM CMP cmp TXM options de calcul des flux hydrauliques l ments THM cf 82 5 8 TYP E OPTION FLUX OPTION FLEN ELGA Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur CA
8. TYP E OPTION CRITE OPTION SI SI EP T ELNO ELGA ELNO LGA ELNO ELGA ELGA DO _ELNO SITQ ELNO EPTQ ELNO CRIT_ELNO l T l T l l l l J El EI EL E O O O l E P P P EN N I options de calcul d indicateurs d erreur cf 2 5 6 TYP E OPTION INDI ERRE OPTION SIZ1 NOEU ERZ1 ELEM SIZ2 NOEU ERZ2 ELEM ERME E EM ERME ELNO autres options cf 2 5 7 TYPE OPTION YAUTRES OPTION SPMX ELGA NOM CHAM ch cham elem NOM CMP cmp TXM PRES DBEL DEPL VNOR ELEM DEPL VARC ELGA Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Fo Titre Op rateur CALC_ELEM Date 20 06 2011 Page 6 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 options pour les r sultats non lin aires produits par STAT_NON LINE ou DYNA NON LINE RESULTAT resu evol noli TYPE OPTION TOUTES DEFAUT OPTION toutes les options ci dessous options de calcul des contraintes l ments de milieux continus 2D et 3D cf 2 5 1 TYPE
9. MOY ordonn e moyenne de la couche feuillet moyen 0 PLAN ISMAIL DEFAUT MOY ESINE SUP e MAIL plan du maillage MOY plan moyen e INF plan sup rieur dans le sens de la normale e SUP plan inf rieur dans le sens de la normale Cet op rande permet de sp cifier le plan de calcul des champs l mentaires pour un mod le avec des l ments de plaques en tenant compte de l excentrement ventuel Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur CALC_ELEM Responsable Aimery ASSIRE Date 20 06 2011 Page 11 28 Cl U4 81 01 R vision 6579 Limitations cette option n est disponible que pour le calcul des efforts g n ralis s par l ments aux n uds partir des d placements lasticit lin aire option De plus EFGE ELNO cette option n est utilisable que pour les DKT DST Q4G GRILLE ANGLE delta I Ho 0 DEFAUT e delta angle en degr s valeur enti re compt partir de la position de la g n ratrice de l ment tuyau 0 VECTEUR x y z 1 R par d faut VECTEUR 1 0 0 ANGL REP a amp 1 _R Mots cl s permettant la construction d un rep re local aux l ments de coques ou de plaques afin de calculer les champs produits par
10. Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster default Titre Op rateur CALC_ELEM Date 20 06 2011 Page 24 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 QIZ2 ELEM Cette option R4 10 04 vise am liorer le traitement des singularit s dans les strat gies d adaptation de maillage en l occurrence avec HOMARD En pratique les indicateurs d erreur sont lev s dans les zones singuli res si bien que rapidement seules les zones singuli res sont raffin es et masquent donc les autres zones sensibles zones fort gradient que l on souhaiterait raffiner Cette option est un champ constant par l ment et comporte trois composantes 1 DEGRE qui correspond la d tection des l ments finis singuliers En pratique cette composante vaut le degr d interpolation des l ments finis choisis si l l ment fini n est connect aucune singularit et vaut l ordre de la singularit si l l ment fini est connect un n ud consid r par la m thode comme singulier par exemple pour un l ment voisin de la pointe d une fissure cette valeur vaut 0 5 2 RAPPORT qui correspond la carte de modification de taille des l ments finis en cas de remaillage pour une erreur globale donn e Cette composante est gale au rapport entre la nouvelle taille de l l ment fini et la taille actuelle 3 TAILLE qui corr
11. EPOT _ELEM 3 5 Calcul de l endommagement de Lema tre ou de Lema tre Sermage evolas CALC ELEM reuse evolas O pol ESULTAT evolas PTION ENDO ELGA ENDO ELNO Manuel d utilisation Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Fascicule u4 81 Outils g n raux
12. OPTION SIGM MASSIF OPTION SIEF E NO PROJ ELEM options de calcul des contraintes l ments de structure poutres tuyaux coques cf 2 5 1 TYPE OPTION SIGM STRUCT OPTION SIEF_ELNO EFCA ELNO _ SITU ELNO SICO ELNO options de calcul des d formations cf 2 5 2 TY P E OPTION EPS OPTION EPSI ELNO EPSI ELGA EPSG ELNO EPSG ELGA EPME ELNO EPME ELGA EPMG ELNO EPMG ELGA EPSP ELNO EPSP ELGA EPFD ELNO ELGA ELNO EPFP ELGA EPVC ELNO EPVC ELGA EPTU ELNO DEGE ELNO EPE EPEF Q options d interpolation et d extraction des variables internes cf 2 5 3 TYPE OPTION VAR I OPTION VART ELNO VATU ELNO VACO ELNO EXTR ELGA EXTR ELNO Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster E Titre Op rateur CALC ELEM Date 20 06 2011 Page 7 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 options de calcul d nergies cf 2 5 4 TYPE OPTION MENER OPTION ETOT ELGA ETOT_ELNO
13. cl s GROUP MA MAILLE TOUT facultatif avec TOUT OUI en d faut permet de d finir le rep re local de post traitement et la localisation groupe de mailles par groupe de mailles exemple faire le calcul des contraintes sur le tuyau coude en une seule fois MAILLE MAILLE Imail Ce mot cl permet d appliquer les hypoth ses sur les mailles dont la liste est indiqu e en argument TOUT TOUT OUI Ce mot cl permet d appliquer les hypoth ses sur toutes les mailles du maillage GROUP MA GROUP MA gma Ce mot cl permet d appliquer les hypoth ses sur le groupe de mailles gma indiqu en argument 0 NUME COUCHE nume Dans le cas d un mat riau multicouche coque multicouche d finie par DEFI COQU MULT ou d un l ment de structure avec comportement non lin aire local int gr par couches NUME_COUCHE est la valeur enti re comprise entre 1 et le nombre de couches n cessaire pour pr ciser la couche o l on d sire effectuer le calcul l mentaire Par convention la couche 1 est la couche inf rieure dans le sens de la normale dans le cas des l ments de coque m canique ou de coque thermique et correspond la couche interne dans le cas d un l ment TUYAU NIVE COUCHE Pour la couche nume d finie par NUME COUCHE permet de pr ciser l ordonn e o l on d sire effectuer le calcul l mentaire INF ordonn e inf rieure de la couche peau interne SUP ordonn e sup rieure de la couche peau externe
14. du mat riau et P e le seuil d endommagement li l nergie stock e dans le mat riau si s 1 on obtient la loi de Lema tre classique avec Y Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster M Titre Op rateur CALC_ELEM Date 20 06 2011 Page 19 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 Calcul syst matique du taux de triaxialit et de la contrainte quivalente d endommagement g SI ENDO o o 2 1 v 3 1 2v e On TRIAX Q Oog ll s o tr o I 3 avec g 3 S S 2 l o 3 tr o Calcul du dommage total par cumul lin aire D_CUMULE D D 3 TRIAX valeur du taux de triaxialit SI_ENDO valeur de la contrainte d endommagement de Lema tre Sermage COENDO valeur de la contrainte d endommagement de Lema tre Sermage normalis e DOM LEM valeur du dommage de Lema tre Sermage valeur du dommage de Lema tre Sermage cumul w a G a Ra ENDO ELNO Endommagement de Lemaitre Sermage obtenues par extrapolation aux n uds des valeurs aux points de Gauss cf ENDO ELGA l EPEQ ELGA EPMQ ELGA D formations quivalentes aux points de Gauss calcul es partir des champs EPSI ELGA 0U EPME ELGA INVA 2 d formation quivalen
15. ments de poutres EPOT UT K U UT B Ae 7e Ae e et pour les l ments de plaques et coques EPOT U K U U B 4e 2 5 5 Options de calcul de crit res CRIT ELNO Calcul des crit res de rupture pour les coques en mat riaux composites R4 01 01 A partir des contraintes calcul es pour une couche donn e option SIGM ELNO et mots cl s NUME COUCHE et NIVE COUCHE et des contraintes limites fournies sous ELAS ORTH dans DEFI MATERIAU le champ CRIT _ELNO contient 6 composantes ex Fa CRIL crit re de rupture en traction dans le sens L si o gt 0 Xr T tou CRILP crit re de rupture en compression dans le sens L si o lt 0 C Opo a CRIT crit re de rupture en traction dans le sens T si gt 0 Y O A 3 CRITP crit re de rupture en compression dans le sens T si op lt 0 Yo o lorr us oni CRILT crit re de rupture en cisaillement SET et CRITH crit re de Tsai Hill voir exemple dans test SSLS121 V3 03 121 Toutes ces quantit s sont calcul es dans le rep re d orthotropie de la coche consid r e l ENDO ELGA Calcul du dommage d aux points de Gauss partir du tenseur des contraintes et de la d formation plastique cumul e p La cin tique d endommagement est donn e par la loi de Lema tre Sermage ne a P SI P Psenil Ge 2E 1 D o et s sont des coefficients caract ristiques
16. C est une extraction effectu e suivant le m me principe que l option d j existante SITU ELNO Le calcul des d formations s effectue dans une couche et un secteur angulaire d l ments tuyau PMPB_ELGA PMPB_ELNO Calcul de crit res du RCC M G3000 pour les l ments de poutres POU D E et POU D T Deux quantit s sont calcul es PM et PMPB N PM 5l purs nl MER avec M M M Ceci correspond la valeur maximum de SIXX dans une section circulaire R3 08 01 PMPB_ELGA valeurs de PM et PMPB aux points de Gauss calcul es partir de SIEF ELGA PMPB_ELNO valeurs de PM et PMPB aux n uds calcul es partir de SIEF ELNO En toute rigueur ces crit res sont appliquer aux contraintes primaires Cette distinction est faire par l utilisateur 2 5 6 Options de calcul d indicateurs d erreur DERA ELGA Indicateur local de d charge et indicateur de perte de radialit aux points de Gauss R4 20 01 DERA ELNO Indicateur local de d charge et indicateur de perte de radialit aux n uds R4 20 01 Remarque Pour les options DERA ELGA et DERA ELNO il faut savoir que le calcul n cessite de comparer les champs de contraintes aux instants t et t Le r sultat est rang au num ro d ordre associ l instant t lell llo loll Par d faut le calcul se fait pour les num ros d ordr
17. Code Aster PA Titre Op rateur CALC_ELEM Date 20 06 2011 Page 1 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 Op rateur CALC ELEM 1 But Cr er ou compl ter un r sultat en calculant des champs par l ments contraintes d formations Chaque champ l mentaire d sir est caract ris par le mot cl OPTION SIGM ELNO FLUX ELGA VARI ELNO Le concept r sultat produit est soit cr soit modifi c est dire que l appel CALC ELEM se fait de la fa on suivante resu CALC ELEM RESULTAT resu reuse resu ou bien resul CALC ELEM RESULTAT resu Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster ui Titre Op rateur CALC_ELEM Date 20 06 2011 Page 2 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 Table des Mati res a E T nn eue italiennes intimiste E AAE 1 PRE Ale E A A E A E ET EL A A A A SEA 3 2 1 Op randes RESULTAT MODELE CHAM_MATER CARA_ELEM EXCIT SOLVEUR 9 21 1 Op randes RESULTAT reniris aaa aaa aA onnaa 9 2 1 2 Op randes MODELE CHAM_MATER CARA ELEM isssseeeennnnneennnnns 9 2 1 3 Motell EXCH occidi a a a a aaa 9 2 1 4Motcl SOLVEUR seb uen fetr detesmenu pen ne genres thin tt
18. LC_ELEM Date 20 06 2011 Page 8 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 options thermiques OPTION FLUX ELNO FLUX ELGA ERTH ELEM ERTH ELNO SOUR ELGA DURT ELGA META DURT ELNO HYDR ELNO RESULTAT resu evol_ther options acoustiques OPTION PRAC ELNO INTE ELNO RESULTAT resu acou harmo __ mode acou TITRE titre 1 Kn INFO 1 DEFAUT J 2 Manuel d utilisation Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Fascicule u4 81 Outils g n raux Code Aster Ho Titre Op rateur CALC_ELEM Date 20 06 2011 Page 9 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 2 1 Op randes RESULTAT MODELE CHAM MATER CARA ELEM EXCIT SOLVEUR 2 1 1 Op randes RESULTAT RESULTAT resu Nom de la structure de donn es r sultat enrichir Cet argument peut tre le m me que celui utilis pour le concept enrichi par l op rateur ou un nom diff rent ce qui cr era une nouvelle structure de donn es r sultat voir par exemple le test SSLS504 V3 03 5041 Remarque dans la majorit des situations la structure de donn es resu contient toutes les informations n cessaires au calcul des options le mod le le champ de mat riau les caract ristiques l mentaires les chargements Les mots cl s MODELE CHAM MATER CARA ELEM et EXCIT sont donc inutiles
19. N naine tent i a tient 27 SHEXEMPISS nn haleine Na aa aaa aeaaee ian iaaa 27 31 Calcul d fl x po r un evol TAB sur ninoi naen ai i Aa 27 3 2 Calcul de l estimateur d erreur ZZ2 pour quelques instants d un concept de type evol_elas 28 3 3 Contraintes aux points de GAUSS pour un calcul thermo m canique 28 3 4 Calcul des nergies potentielles pour un mode propre Mine 28 3 5 Calcul de l endommagement de Lema tre ou de Lema tre Sermage 28 Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster A Titre Op rateur CALC _ELEM Date 20 06 2011 Page 3 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 2 Syntaxe resu CALC_ELEM reuse resu MODELE mo modele CHAM MATER chmater cham mater 0 CARA ELEM carac cara elem 0 SOLVEUR F voir le document U4 50 01 0 EXCIT F CHARGE l charge l char meca COEF MULT cm R COEF MULT C cmc C FONC_MULT fm fonction FONC MULT C fmc fonction C PHAS DEG pd R Q PUIS_PULS n E TYPE CHARGE FIXE O S lection des mailles concern es par le calcul TOUT vourr DEFAUT GROUP MA l grma 1 gr mail
20. ON SIEF ELNO SIGM ELNO SIEF ELGA PROJ ELEM options de calcul des contraintes l ments de structure poutres tuyaux coques CE CS2 5 11 TYPE OPTION SIGM STRUCT OPTION SIEF ELNO SIGM ELNO SIEF ELGA SITU ELNO SIPO ELNO EFGE ELNO EFCA ELNO SICA ELNO options de calcul des flux hydrauliques l ments THM ct 11 S2 55811 TYP E OPTION FLUX OPTION FLEN ELGA options de calcul des d formations cf 2 5 2 TYP E OPTION EPSI OPTION EPSI ELNO EPSI ELGA EPME ELNO EPME ELGA DEGE ELNO EPTU_ELNO EPVC_ELNO EPVC_ELGA options de calcul d nergies cf 2 5 4 TYP E OPTION MENER OPTION EPOT ELEM ECIN ELEM ENE E GA ENE E NO ETOT ELGA ETOT ELNO DISS ELGA DISS ELNO Manuel d utilisation Fascicule u4 81 Outils g n raux CI Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster M Titre Op rateur CALC_ELEM Date 20 06 2011 Page 5 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 ETOT ELEM options de calcul de crit res cf 2 5 5
21. Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur CALC_ELEM Responsable Aimery ASSIRE Cl U4 81 01 R vision Pour une section rectangulaire e on calcule la valeur de SMFY en z HZ 2 e on calcule la valeur de SMFZ en y HY 12 Date 20 06 2011 Page 13 28 6579 Pour une section circulaire on calcule les valeurs de SMFY et SMFZ pour y et z valant R PROJ ELEM Calcul des champs de contraintes par l ments sur les parements surfaces amont et aval d un ouvrage hydraulique alors que la structure est mod lis e en volumique Les composantes de ces champs sont les suivantes Contraintes normales aux faces des l ments calcul es aux centres des faces SIG NX composante suivant X dans le rep re global de la contrainte no rmale SIG NY composante suivant Y dans le rep re global de la contrainte normale SIG NZ composante suivant Z dans le rep re global de la contrainte normale SIG N contrainte normale Contraintes tangentes aux faces des l ments calcul es aux centres des faces SIG TX composante suivant X dans le rep re global de la contrainte tangentielle SIG TY composante suivant Y dans le rep re global de la contrainte tangentielle SIG TZ composante suivant Z dans le rep re global de la contrainte tangentielle Rosettes des contrainte
22. P pour les b tons EPFD ELNO EPFD ELGA Calcul aux n uds ou aux points de Gauss des d formations de dessiccation des b tons pour le mod le BETON UMLV FP EPME EPME ELNO ELGA Calcul aux n uds ou aux points de Gauss des d formations m caniques partir des d placements Ce calcul est fait en th orie des petits d placements Les d formations calcul es sont gales aux d formations totales moins les d formations thermiques m th E u z t yt uja E E PMG EPMG ELNO ELGA Calcul aux n uds ou aux points de Gauss des d formations m caniques partir des d placements Ce calcul est fait en th orie des grands d placements Les d formations calcul es sont gales aux d formations totales moins les d formations thermiques m th u u ptu ttg tig j mu Le iT Uk Ug jT E ELGA D formations de Green Lagrange aux points de Gauss EPSG E u z Mi j TUy itU iUe j ELNO D formations de Green Lagrange aux n uds EPSI EPSI EPSG ELNO ELGA Calcul des d formations par l ment aux n uds ou aux points de Gauss partir des d placements 1 ey u 5 ytu Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl ht
23. PL Projection d un champ de vitesse sur la normale des l ments de type coque ou plaque Cette option sert notamment au cha nage avec le code VARIA VARC ELGA Calcul des variables de commandes ayant servi un calcul m canique 7 variables sont syst matiquement calcul es TEMP HYDR SECH CORR IRRA NEUTI NEUT2 Remarque Les variables qui n ont pas t d finies sont initialis es la valeur R8VIDI nombre r el tr s grand de l ordre de 1 D308 CE SPMX ELGA NOM CHAM ch cham elem NOM CMP cmp TXM Extraction des valeurs extr males en chaque point de Gauss lin ique d un l ment de tuyau de la composante cmp du champ ch sur tous les points d int gration de la section Les champs possibles sont les champs de contraintes SIEF ELGA SIEF ELGA les champs de d formations EPSI ELGA les champs de valeurs quivalentes SIEQ ELGA EPEQ ELGA les champs de variables internes VARI ELGA Le champ cr e de nom SPMX ELGA contient pour chaque instant les composantes IN valeur minimum AX valeur maximum NCOUMIN num ro de la couche pour la valeur min NCOUMAX num ro de la couche pour la valeur max NSEGMIN num ro du secteur angulaire pour la valeur min NSEGMAX num ro du secteur angulaire pour la valeur max NPCOUMIN num ro du point de la couche NCOUMIN NPCOUMAX num ro du poin
24. acoustiques 2 7 1 Op rande OPTION PRAC_ELNO Calcul de la pression aux n uds en partie r elle partie imaginaire et d cibels INTE ELNO Calcul de l intensit acoustique active et r active aux n uds Les d finitions se trouvent dans R4 02 01 2 8 Op rande TITRE Q TITRE titre Titre que l on veut donner au r sultat de la commande U4 02 01 3 Exemples 3 1 Calcul du flux pour un evol_ther evoth CALC ELEM reuse evoth RESULTAT evoth TOUT_ORDRE COUT Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur CALC ELEM Responsable Aimery ASSIRE OPTION FLUX ELNO Date 20 06 2011 Page 28 28 Cl U4 81 01 R vision 6579 3 2 Calcul de l estimateur d erreur ZZ2 pour quelques instants d un concept de type evol_elas evolas CALC ELEM reuse evolas RESULTAT evolas INST 1 10 20 OPTION ERZ2 ELEM 3 3 Contraintes aux points de GAUSS pour un calcul thermo m canique evolas CALC ELE reuse evolas RESULTAT evolas TOUT ORDRE OUI OPTION SIEF ELGA 3 4 Calcul des nergies potentielles pour un mode propre mode CALC ELE reuse mode RESULTAT mode NUME MODE 3 OPTION
25. de chaleur aux n uds partir de la temp rature ERTH_ELEM ERTH_ELNO Estimateurs d erreur en r sidu en thermique R4 10 03 II faut pr alablement effectuer dans CALC ELEM le calcul des flux aux n uds via FLUX ELNO Le mot cle INFO procure tous les affichages interm diaires du calculs connectivit s normales diam tres valeurs des champs jacobien L option ERTH_ELNO permet de ramener le champ par l ment ERTH_ELEM a un champ aux n uds par l ment ce qui permet de faire des relev s de valeurs ou des impressions visualisations SOUR_ELGA Calcul d une source de chaleur pouvant tre introduite dans un calcul thermique via le mot cl SOURCE SOUR_CALCULEE de la commande AFFE CHAR THER U4 44 02 Cette source est calcul e partir d un potentiel lectrique via la loi d Ohm Ce potentiel lectrique doit avoir t calcul par l op rateur THER LINEAIRE U4 54 01 en faisant les analogies n cessaires DURT ELGA META DURT ELNO Calcul de duret aux points de Gauss ou aux n uds partir des phases m tallurgiques cf R4 04 011 HYDR ELNO Calcul de lhydratation aux n uds partir de l hydratation aux points de Gauss calcul e par THER NON LINE pour la mod lisation du b ton R7 01 121 2 7 Op randes pour les options
26. e 1 n 1 L indicateur de d charge est calcul par ID Si on pr cise la liste d instant avec des trous ventuellement le calcul ne concernera que les instants demand s mais il comparera toujours l instant t avec l instant f dans la liste des instants ayant servi faire le calcul non lin aire ERZ1 ELEM respectivement ERZ2 ELEM Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster default Titre Op rateur CALC_ELEM Date 20 06 2011 Page 22 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 Calcul de l estimateur d erreur de ZHU ZIENKIEWICZ lasticit lin aire 2D partir de l option SIZ1_ NOEU respectivement SIZ2 NOEU Si ce dernier champ n existe pas dans resu il est automatiquement construit au pr alable voir R4 10 01 ERME ELEM Estimateur d erreur en r sidu en m canique R4 10 02 et en hydro m canique stationnaire R4 10 04 calcul par l ment Conseils d utilisation de l option ERME ELEM Pour bien effectuer l estimation d erreur du calcul m canique dans les limites th oriques de la formule mise au point dans le cadre elliptique avec fronti re r guli re il faut l effectuer sur tout le mod le TOUT OUI Valeur par d faut A noter que le mod le n est pas forc ment d fin
27. e d un l ment entre l instant courant et l instant pr c dent Lorsque cette option est appliqu e tout le maillage et pour des comportements diff rents de VMIS ISOT_ LINE et VMIS ISOT TRAC la somme de chacun des incr ments d nergie donne le m me le r sultat que celui obtenu par la commande ENER TOTALE de l op rateur POST _ELEM Pour les deux comportements VMIS ISOT_ LINE et VMIS ISOT TRAC le r sultat est diff rent puisque l op rateur POST ELEM utilise un traitement sp cifique pour le calcul de l nergie de d formation totale DISS ELNO DISS ELGA Calcul de l nergie de dissipation aux points de Gauss ou aux n uds de chaque l ment Valable uniquement pour les l ments DKTG et la loi GLRC DM Leurs expressions sont donn es dans la doc R de la loi de comportement EPOT ELEM Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster nn Titre Op rateur CALC_ELEM Date 20 06 2011 Page 18 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 Calcul de l nergie potentielle de d formation int gr e sur un l ment partir des d placements U et des temp ratures T e pour les l ments de milieux continus 2D et 3D EPOT gt f E U AE U dy f e U Ae U dv 3 f U A U dv element element element e pour les l
28. e ee atse tn taimetue ann e uut 9 2 2 S lection des mailles concern es par le calcul nn 9 2 3 D finition d un rep re local de d pouillement op rande REPE_COQUE 10 24 S lection des Num ros d OFUFS 278 hunbistmennnbenetrermnnedeni anaana a ae a t nee het onted es 11 2 5 1 Option d calcul des ContraiNteS isisi iiianuasii tiaina aana iaaa iranienne 11 2 5 2 Options de calcul des d formations iiisisiciniuinnininonini ninina aaae aidia 14 2 5 3 Options d interpolation et d extraction des variables internes 15 2 54 Options de caleul d nergie nta nas Penebnnaats trame nann ia tnt ientmnpia errant nes 16 255 Options d Calsulide CnT fes sn inner naet eitotnennnteeDabhend din nieduner saura de 17 2 5 6 Options de calcul d indicateurs Qem eUT smriti aiaiai ataa 21 2 9 l AUTTES OPDNS ariiraa niai aieiaa aiaia aada aiaiai 24 2 5 8 Option de calcul des flux hydrauliques l ments THM 25 2 5 JOPD ARE NORMES Reste ana ensneaean nn eme ete idae re 25 2 6 Op randes pour les options 1hermIqUes 112126480800120mavmnhednitpsesiahansrrenteelanantanueanc das 27 2 6 1 Op rand OPTION si snntandirtinennhuetaeasithnnnesnsintingencantatanaiarennteneri alien aia 27 2 7 Op ra des pour les options ACOUSUIQU SE eisini none trade iiiaae aaa daa 27 2 1 Op rand OPTION cucina i 27 2 8 Operando TRE aisis ten nee nt nets A
29. es T tant la temp rature et le coefficient de dilatation e le retrait de s chage EPSECE utilis pour les lois d crivant le comportement du b ton EU K yselS S S tant la variable de commande s chage et K le coefficient de retrait de dessiccation e le retrait d hydratation EPEYDR utilis pour les lois d crivant le comportement du b ton E B ondog gt h tant la variable de commande hydratation et B tant le coefficient de retrait endog ne dessic hydr _ E endog 2 5 3 Options d interpolation et d extraction des variables internes VARI ELNO Calcul des variables internes aux n uds des l ments partir des points de Gauss Le nombre et le type de ces variables internes sont sp cifiques chaque mod le de comportement cf doc U4 de STAT NON LINE par exemple VATU ELNO Calcul des variables internes dans une couche et pour un secteur angulaire d l ments tuyau mots cl s REPE COQUE NUME COUCHE NIVE COUCHE et ANGLE VACO ELNO Calcul des variables internes dans une couche d l ments coque d finie par NUME COUCHE et NIVE COUCHE Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur CALC ELEM Responsable Aimery ASSIRE EXTR ELNO
30. espond la carte des nouvelles tailles des l ments finis en cas de remaillage pour une erreur globale donn e Cette donn e est directement utilisable par certains mailleur GMSH par exemple Cette option peut s utiliser selon deux sch mas Les l ments finis consid r s comme singuliers par la m thode peuvent tre exclus du processus de d coupage en leur affectant par exemple une erreur nulle ela nouvelle taille des l ments finis est donn e un remailleur en occurrence HOMARD pour Code _Aster pour que celui ci construise le nouveau maillage en respectant au mieux cette nouvelle carte de taille Actuellement le logiciel HOMARD d coupe une fois l l ment par exemple en 2D un triangle est divis en 4 mais pas plus Pour continuer le d coupage il faut faire appel de nouveau HOMARD Une volution est donc pr voir pour qu on puisse diviser plusieurs fois un l ment et donc respecter au mieux la carte de taille du nouveau maillage Le calcul de cette option n cessite au pr alable le calcul d un indicateur d erreur c est la composante absolue qui est utilis e et c est cod en dur dans Aster et de l nergie de d formation totale Dans le cas o l une de ces options n est pas calcul e un message d alarme est mis et l option SING ELEM n est pas calcul e Pour l indicateur d erreur quatre choix sont possibles e ERME_ELEM pour l indicateur en r sid
31. i sur toute la g om trie Il faut aussi effectuer pr alablement dans CALC ELEM le calcul des contraintes aux n uds cf R3 06 03 par SIGM ELNO en lin aire par SIEF ELNO en non lin aire Sinon une alarme est mise et le calcul d erreur n est pas effectu sans provoquer l arr t de l ex cution Si le champ de contraintes aux n uds existe d j dans la structure de donn es resultat il n est pas recalcul En ce qui concerne les chargements Il faut fournir CALC ELEM les chargements utilis s pour le calcul m canique EXCIT F CHARGE en prenant bien garde aux r gles de surcharges diff rentes pour le solveur m canique et pour cette option de CALC ELEM Ainsi le calcul m canique MECA STATIQUE STAT NON LINE agr ge les conditions aux limites alors que le calcul de l erreur ne va retenir pour un type de conditions aux limites donn que la derni re list e dans le EXCIT de CALC ELEM L ordre a donc une importance cruciale Il ne faut donc pour un type de conditions aux limites qu une seule occurrence dans les AFFE CHAR On ne tient compte que des chargements de type PESANTEUR ROTATION FORCE INTERNE PRES REP FORCE FACE FORCE ARETE Seules les trois derniers peuvent tre Variables Il est conseill d utiliser des l ments finis d ordre deux dans le cas de f
32. ini par l utilisateur dans la commande AFFE CARA ELEM Dans le cas des coques en grands d placements et grandes rotations COQUE 3D avec DEFORMATION GROT GDEP cette option int gre galement le calcul des contraintes de Cauchy partir des contraintes de Piola Kirchhoff Les contraintes issues de cette option sont donc des contraintes de Cauchy dans une couche SITU ELNO Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster on Titre Op rateur CALC_ELEM Date 20 06 2011 Page 12 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 Calcul des contraintes dans une couche et pour un secteur angulaire d l ments tuyau mots cl s REPE COQUE NUME COUCHE NIVE COUCHE et ANGLE SICA ELNO EFCA ELNO Changement de rep re des contraintes ou des efforts g n ralis s par l ment aux n uds du rep re local au rep re global de description du maillage cette option consiste convertir un champ de contrainte ou d efforts g n ralis s pour un mod le avec des l ments de structure attach s au rep re de r f rence d un ensemble de plaques ou de coques ou du rep re d inertie principal d un l ment de poutre pour les exprimer dans le rep re global EFGE ELNO 0 REPE COQUE 0 PLAN MAIL
33. le MAILLE 1 mail 1 maille S lection des num ro d ordre TOUT ORDRE OUT y NUME ORDRE 1 nuor LI LIST ORDRE 1 nuor listis NUME MODE 1 numo L T NOEUD CMP l_nomo 1 K16 NOM CAS nocas K16 INST l inst L R FREQ l freq LR LIST_INST l inst listr8 LIST_FREQ 1 freq listr8 O P RECISION prec 1 0E 3 DEFAUT CRITERE RELATIF DEFAUT ABSOLU 0 REPE COQUE TOUT OUT z DEFAUT MAILLE mail l maille GROUP MA gma F group ma ANGLE delta 1 0 DEFAUT Q PLAN MAIL DEFAUT MOY YINF 7 SUR NUME COUCHE nume I T DEFAUT NIVE COUCHE A INE SUP MOY DEFAUT 0 ANGLE REP p 1 R VECTEUR x y 2 1 R Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster ss Titre Op rateur CALC_ELEM Date 20 06 2011 Page 4 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 options pour des r sultats m caniques lin aires RESULTAT resu TYPE OPTION TOUTES DEFAUT OPTION toutes les options ci dessous options de calcul des contraintes l ments de milieu continu 2D et 3D cf S255 17 TYPE OPTION SIGM MASSIF OPTI
34. les contraintes locales en un point de l paisseur aux n uds SIGM ELNO en lin aire et SICO ELNO en non lin aire INDL ELGA Indicateur de localisation bas sur le tenseur acoustique crit re de RICE d fini par det N H N lt 0 o H d signe l op rateur tangent et N la normale aux directions de localisation Cet indicateur d fini un tat partir duquel le probl me local d int gration du comportement perd son caract re d unicit La m thode n est d velopp e que dans le cas 2D et pour la loi de comportement de type DRUCKER PRAGER L option INDL ELGA contient les composantes suivantes INDICE Indicateur de localisation valant O si det N HN gt 0 et valant 1 sinon ce qui correspond a l initiation de la Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Mat Titre Op rateur CALC_ELEM Date 20 06 2011 Page 21 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 localisation DIRI correspond la premi re normale la zone de localisation DIR2 la deuxi me normale DIR3 la troisi me normale DIRA4 la quatri me normale EPTQ ELNO l SITQ ELNO Calcul des d formations g n ralis es et des contraintes pour un l ments tuyau Ce sont des valeurs quivalentes de type SIEQ ELGA et EPEQ ELGA en un point de la section
35. les options demand es contraintes efforts dans ce rep re local La d finition de ce rep re local est identique celle de l op rateur AFFE CARA E EM U4 42 01 2 4 S lection des num ros d ordre L emploi des mots cl s TOUT ORDRE NUM ORDRE INST FREQ est d crit dans le document U4 71 00 2 5 Op randes pour les options m caniques 2 5 1 Option de calcul des contraintes l SIEF_ELGA Calcul de l tat de contrainte par l ment aux points d int gration de l l ment points de GAUSS ou points d d int grations pour chaque couche des l ments de coque et chaque secteur des l ments tuyaux partir des d placements lasticit lin aire voir U2 01 05 ST EE _ELNO Calcul de l tat de contrainte aux n uds par l ment partir de l tat de contrainte aux points de Gauss n SIGM_ELNO Calcul des contraintes par l ment aux n uds partir des d placements lasticit lin aire SICO ELNO NUME COUCHE nume 1 DEFAUT NIV_COUCHE INF SUP MOY DEFAUT Calcul des contraintes dans une couche d l ments de coque mots cl s NUME COUCH NIVE COUCHE partir des contraintes aux points d int gration de chaque couche E et SIEF ELGA calcul es lors d un calcul non lin aire Ces contraintes sont calcul es dans le rep re local de la coque d f
36. ml Code Aster Mat Titre Op rateur CALC_ELEM Date 20 06 2011 Page 15 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 Pour les l ments de structure ces d formations sont obtenues dans le rep re local de l l ment Pour les plaques et coques elles sont calcul es dans la couche et l altitude demand e dans le mot cl REPE COQUE cf 2 3 EPTU ELNO Calcul des d formations dans une couche et pour un secteur angulaire d l ments tuyau voir le mot cl REPE COQUE EPSP_ELGA D formations an lastiques aux points de Gauss A partir du champ de d placements u de contraintes o de temp ratures T de d formations an lastiques ventuelles et de variables internes on calcule chaque instant e u 4 o e T e e est la d formation de fluage propre de Granger o EPSP ELNO D formations an lastiques obtenues par extrapolation aux n uds des valeurs aux points de Gauss cf EPSP ELGA EPVC ELNO EPVC ELGA Calcul des d formations aux noeuds ou aux points de Gauss li es aux variables de commande Pour l instant ne sont d finies que les composantes suivantes e d formations thermiques EPTHER L EPTHER T EPTHER N telle que E a T T ie L T N sile mat riau est isotrope les 3 composantes sont gal
37. non e calcul sur les mailles volumiques MODELISATION 3D ou 3D_ SI de suivante Calcul des champs de contraintes aux n uds des faces des l ments 3D Ces champs sont valu s partir d un champ de contraintes aux n uds par l ment lin aire la fa on Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Version default Date 20 06 2011 Page 14 28 Cl U4 81 01 R vision 6579 Code Aster Titre Op rateur CALC ELEM Responsable Aimery ASSIRE Moyennation de chacune des composantes du tenseur des contraintes au centre des faces d l ments Projection du tenseur des contraintes suivant un vecteur normal 51G N la face et suivant un vecteur le plan tangent la face SIG T Diagonalisation du tenseur des contraintes dans le plan tangent SIG T1 et SIG T2 2 5 2 Options de calcul des d formations D ELNO EG E Calcul des d formations g n ralis es par l ment aux n uds partir des d placements cette option n a de sens que pour les l ments de structure poutre plaque tuyau Les d formations g n ralis es sont obtenues dans le rep re local de l l ment EPFP ELNO EPFP ELGA Calcul aux n uds ou aux points de Gauss des d formations de fluage propre associ es au mod le GRANGER FP ou au mod le BETON_UMLV_F
38. orces volumiques sinon ce terme est tr s mal calcul puisque DIV SIGMA est quasi nul Pour prendre en compte l erreur relative une CL nulle il faut l imposer en tant que fonction via UN AFFE CHAR MECA F Via une constante elle ne sera pas prise en compte Maillage Le maillage doit tre triangulaire quadrangle t tra drique ou hexa drique avec aucun GROUP NO si on veut remailler ensuite via HOMARD En 2D il ne prend en compte que les erreurs sur et entre les l ments isoparam triques SEG2 3 TRIA3 6 QUAD4 8 9 En 3D idem avec FACE3 4 6 8 9 TETRA4 10 PENTA6 13 15 et HEXA8 20 27 donc pas les PYRAM ni les l ments de structure coque plaque poutre D autre part il faut veiller ne pas intercaler de segments entre deux quadrangles ou deux triangles resp quad ou triangle entre deux hexa sinon on ne peut pas calculer le terme de saut relatif ce voisinage A la place on s enquiert tort d une ventuelle condition aux limites ERM E ELNO Estimateur d erreur en r sidu calcul aux n uds R4 10 02 Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Mat Titre Op rateur CALC_ELEM Date 20 06 2011 Page 23 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 Conseil d utilisation des options 01Z1 ELEM QIZ2 ELEM QIRE ELEM
39. pour la composante voulue du d placement EPSI _INIT d formation unitaire pour la composante voulue du tenseur des contraintes Une fois les deux probl mes r solus on calcule pour chacun des deux l estimateur d erreur classique d sir le m me pour les deux et enfin il faut d finir un nouveau CALC ELEM avec une des options de calcul d estimateur d erreur en quantit d int r t Un exemple d utilisation du calcul de l estimateur en quantit s d int r t bas sur les r sidus peut tre trouv dans le test ssiv113c et d QIRE ELNO Estimateur d erreur en quantit s d int r t bas s sur les r sidus calcul aux n uds SIZ1 NOEU Calcul des contraintes aux n uds lasticit lin aire 2D les contraintes sont obtenues par un lissage global au sens des moindres carr s des contraintes aux points de GAUSS Voir Estimation d erreur par lissage des contraintes R4 10 01 S1Z2 NOEU Calcul des contraintes aux n uds lasticit lin aire 2D les contraintes sont obtenues par un lissage local un patch d l ments au sens des moindres carr s des contraintes aux points de GAUSS voir R4 10 011 SING ELEM PREC ERR err R TYPE ESTI ERME ELEM ERZ1 ELEM ERZ2 ELEM OIRE ELEM OIZ1 ELEM Manuel d utilisation Fascicule u4 81
40. rreur zz1 Les 3 solveurs autoris s sont LDLT MULT FRONT et MUMPS d faut MULT FRONT En toute rigueur le solveur MUMPS n est pas recommand car il ne sait pas encore traiter STOP SINGULIER NON 2 2 S lection des mailles concern es par le calcul Les mots cl s TOUT GROUP MA et MAILLE permettent l utilisateur de choisir les mailles sur lesquelles il souhaite faire ses calculs l mentaires de post traitement TOUT OUI Toutes les mailles porteuses d l ments finis seront trait es C est la valeur par d faut Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Fo Titre Op rateur CALC_ELEM Date 20 06 2011 Page 10 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 GROUP MA 1 grma MAILLE 1l maille Seules les mailles incluses dans 1_grma et ou 1 maille seront trait es 2 3 D finition d un rep re local de d pouillement op rande REPE COQUE Ce mot cl facteur est r p table Il regroupe les mots cl s simples utilis s pour le post traitement des coques et tuyaux NUME COUCHE NIVE COUCHE ANGLE et PLAN et les mot cl s d finissant le rep re local ANGL REP et VECTEUR du mot cl facteur COQUE de la commande AFFE CARA ELEM De plus la pr sence des mot
41. s partir des champs EPSI ELNO ou EPME ELNO INVA 2 d formation quivalente de Von Mises INVA 256 d formation quivalente de Von Mises sign par la trace de PRIN_1 PRIN_ 2 PRIN 3 d formations principales On note que les d formations quivalentes obtenues partir de EPSI ELNO et EPME ELNO sont identiques En effet la diff rence entre les deux tenseurs est un tenseur sph rique d formation thermique Comme la d formation quivalente est obtenue partir du second invariant du d viateur le tenseur sph rique disparait lorsque l on prend le d viateur SIEQ ELNO Contraintes quivalentes aux n uds VMIS contrainte de von Mises VMIS SG contrainte de von Mises sign e par la trace de g TRESCA contrainte de Tresca PRIN_1 PRIN 2 PRIN 3 contraintes principales Pour les l ments de milieux continus 2D et 3D elles sont extrapol es aux n uds partir des contraintes quivalentes calcul es aux points de Gauss elles m mes calcul es partir des champs de contraintes aux points de Gauss SIEF ELGA en lin aire et SIEF ELGA en non lin aire Dans le cas o on calcule ensuite les valeurs moyenn es aux n uds par l option SIEQ NOEU de CALC NO du fait des interpolations on n a pas forc ment VMIS ABS VMIS SG Pour les l ments de coques elles sont calcul es directement sur
42. s valeurs principales des contraintes projet es sur le plan tangent la face d l ment calcul es aux centres des faces premi re valeur des contraintes projet es SIG T1X composante suivant X dans le rep re global de la premi re valeur principale des contraintes projet es sur le plan tangent SIG T1Y composante suivant Y dans le rep re global de la premi re valeur principale des contraintes projet es sur le plan tangent SIG T1Z composante suivant Z dans le rep re global de la premi re valeur principale des contraintes projet es sur le plan tangent SIG T1 premi re valeur principale des contraintes tangentielle Rosettes des contraintes valeurs principales des contraintes projet es sur le plan tangent la face d l ment calcul es aux centres des faces deuxi me valeur des contraintes projet es SIG T2X composante suivant X dans le rep re global de la deuxi me valeur principale des contraintes projet es sur le plan tangent SIG T2Y composante suivant Y dans le rep re global de la deuxi me valeur rincipale des contraintes projet es sur le plan tangent DT SIG T2Z composante suivant Z dans le rep re global de la deuxi me valeur principale des contraintes projet es sur le plan tangent SIG T2 deuxi me valeur principale des contraintes tangentielle option SIGM ELNO dans le cas lin aire ou option SIEF ELNO dans le cas
43. seuil isotrope NB ITER Nombre d it rations internes ARRET Valeur du test local d arr t du processus it ratif NB REDE Nombre de red coupage local du pas de temps SIGNE Signe du produit contract de la contrainte d viatorique par la d formation plastique d viatorique Remarque Lorsque la variable extraire ne fait pas partie des variables internes des lois concern es une alarme est mise mais le champ est tout de m me affect R amp VIDE 2 5 4 Options de calcul d nergie ECIN ELEM nergie cin tique d un l ment ENEL E NO ENEL E GA Calcul de la densit d nergie lastique aux points de Gauss ou aux n uds de chaque l ment Cette option diff re de l option EPOT ELEM qui calcule l nergie de d formation lastique int gr e dans chaque l ment cette nergie tant un scalaire pour un l ment donn Ici on calcule la densit d nergie lastique qui s crit E 50 40 Ce calcul s appuie sur le champ de contraintes aux points de Gauss obtenu par SIEF ELGA OU SIEF ELGA ETOT_ELGA ETOT ELNO Calcul de l incr ment de densit d nergie de d formation totale aux points de Gauss ou aux n uds de chaque l ment entre l instant courant et l instant pr c dent ETOT ELEM Calcul de l incr ment d nergie de d formation total
44. t de la couche NCOUMAX NPSECMIN num ro du point sur le secteur NSECMIN NPSECMAX num ro du point sur le secteur NSECMAX Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Version Code Aster default Titre Op rateur CALC_ELEM Date 20 06 2011 Page 26 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 2 5 8 Option de calcul des flux hydrauliques l ments THM TYPE OPTION FLUX OPTION E FLHN ELGA Calcul des flux hydrauliques p 5M y aux points de Gauss sur les l ments de bord 2D ou 3D partir du vecteur flux aux noeuds l option SIEF_NOEU doit avoir t calcul e au pr alable O M est le vecteur flux hydraulique du composant ij U2 04 05 ij L int grale des flux sur une surface est effectu e dans POST_ELEM par int gration de ce champ ri Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code _ Aster etait default Titre Op rateur CALC_ELEM Date 20 06 2011 Page 27 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 2 6 Op randes pour les options thermiques 2 6 1 Op rande OPTION FLUX ELGA Calcul des flux de chaleur aux points d int gration de GAUSS partir de la temp rature FLUX ELNO Calcul des flux
45. te de Von Mises NVA _2 2 dev e dev ji avec dev 5 tr e INVA 256 d formation quivalente de Von Mises sign e par la trace de amp PRIN_1 PRIN_ 2 PRIN 3 d formations principales Pour les l ments TUYAU les composantes calcul s sont INVA_2 et INVA 25G SIEQ ELGA Contraintes quivalentes aux points de Gauss 1 VMIS contrainte de von Mises VMIS s avec s5073 lo ji vVMIS_SG contrainte de von Mises sign e par la trace de g TRESCA contrainte de Tresca PRIN 1 PRIN 2 PRIN 3 contraintes principales Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster nn Titre Op rateur CALC_ELEM Date 20 06 2011 Page 20 28 Responsable Aimery ASSIRE Cl U4 81 01 R vision 6579 VECT_1 X VECT 1 Y VECT 3 Z contraintes d formations et directions principales uniquement pour les mod lisations ci dessous 3D 3D SI 3D GRAD VARI SHB8 seulement pour les contraintes AXIS AXIS SI AXIS GRAD VARI eD PLAN D PLAN SI D PLAN GRAD EPSI D PLAN GRAD VARI C PLAN C PLAN SI C PLAN GRAD EPSI C PLAN GRAD VARI Pour les l ments TUYAU les composantes calcul s sont VMIS et VMIS SG et leur version sign e SG EPEQ ELNO EPMQ ELNO D formations quivalentes aux n uds calcul e
46. us e ERZ 1 ou 2 ELEM SIGM pour l indicateur de Zhu Zienkiewicz versions 1 ou 2 e QIRE_ELEM pour l indicateur en quantit d int r t bas sur les r sidus e QIZ 1 ou 2 ELEM SIGM pour l indicateur en quantit d int r t bas sur Zhu Zienkiewicz versions 1 ou 2 e Si les six indicateurs sont pr sents et que rien n est pr cis avec TYPE ESTI l indicateur en r sidu ERME ELEM est choisi par d faut message d alarme mis Si les deux indicateurs de Zhu Zienkiewicz sont pr sents on choisit ERZ1 ELEM Pour l nergie de d formation totale on utilise e Avec STAT NON LINE ETOT ELEM qui est l nergie de d formation totale sur un l ment fini valable pour un comportement lastique et pour un comportement lastoplastique VMIS ISOT XXX e Avec MECA STATIQUE EPOT ELEM qui est l nergie potentielle de d formation lastique sur un l ment fini et int gr e partir des d placements et de la temp rature valable uniquement pour un comportement lastique L utilisateur doit galement renseigner le mot cl PREC ERR un message fatal est mis en cas d absence qui permet de calculer la pr cision souhait e sur l erreur globale pour d terminer la carte de modification de taille cf R4 10 04 La valeur de PREC ERR est
Download Pdf Manuals
Related Search