[U4.81.01] Opérateur CALC_ELEM
Contents
1. e 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 cet 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 comprise strictement entre 0 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 Manuel d utilisation Fascicule2. R vision 9958 evol_ther Cl U4 81 01 FLUX_ELNO FLUX ELGA ERTH ELEM ERTH ELNO SOUR ELGA DURT ELGA META DURT ELNO HYDR ELNO PRAC ELNO INTE ELNO acou harmo __ mode acou 1_Kn D EFAUT 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 Version 10 Titre Op rateur CALC_ELEM Date 23 10 2012 Page 9 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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 2 1 2 Op randes MODELE CHAM MATER CARA ELEM MODELE mo Nom du mod le sur lequel sont calcul s les efforts les contrain
3. QIRE ELNO Estimateur d erreur en quantit s d int r t bas s sur les r sidus calcul aux n uds SIZI 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 SIZ2 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 QIRE ELEM QIZ1 ELEM OIZ2 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 me
4. CRIT OPTION E SIE Q ELNO Q ELGA Q ELNO EPEQ ELGA Q Q T ELNO ELGA ELNO ENDO ELGA ENDO ELNO PMPB ELNO F ELGA INDL ELGA E U U B9 w pa options de calcul d indicateurs d erreur cf 2 5 6 TRI E OPTION INDI_ERRE OPTION ERME_ELEM ERZ1_ELEM ERZ2_ELEM OIZ1 ELEM OIZ2 ELEM SING ELEM PREC ERR SING ELNO ERME ELNO err R 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 10 Titre Op rateur CALC ELEM Responsable Aimery ASSIRE options thermiques OPTION RESULTAT resu options acoustiques OPTION RESULTAT resu TITRE titre INFO T 2 Date 23 10 2012 Page 8 27
5. 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 Cod e AS ter Version 10 Titre Op rateur CALC ELEM Date 23 10 2012 Page 4 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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
6. OPTION 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 AS ter Version 10 Titre Op rateur CALC_ELEM Date 23 10 2012 Page 5 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 ETOT ELEM options de calcul de crit r
7. 23 10 2012 Page 13 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 Pour une section rectangulaire on calcule la valeur de SMFY en z HZ12 on calcule la valeur de SMFZ en y HY 2 Pour une section circulaire on calcule les valeurs de SMFY et SMFZ pour y et z valant R PROJ EI EM 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 normale S1G 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 S1G 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 S1G TZ composante suivant Z dans le rep re global de la contrainte 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 premi re valeur des contraintes projet es SIG T1X comp
8. 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 AS ter Version 10 Titre Op rateur CALC_ELEM Date 23 10 2012 Page 7 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 options de calcul d nergies cf 2 5 4 TYPE OPTION MENER OPTION ETOT ELGA ETOT_ELNO ETOT _E EM ENEL ELGA ENEL ELNO options de calcul de crit res cf 2 5 5 TYPI E OPTION
9. Les champs possibles sont les champs de contraintes SIEF ELGA SIFF champs de d formations EPSI ELGA les champs de valeurs quivalentes SI EPEQ ELGA les champs de variables internes VARI ELGA ELGA les FQ 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 point de la couche NCOUMAX NPSECMIN num ro du point sur le secteur NSECMIN NPSECMAX num ro du point sur le secteur NSECMAX 2 5 8 Option de calcul des flux hydrauliques l ments THM TYPE OPTION FLUX OPTION FLHN_ELGA Calcul des flux hydrauliques P M v aux points de Gauss sur les l ments de bord 2D ou 3D partir du vecteur flux aux noeuds l option STIEF NOEU doit avoir t calcul e au pr alable O M est le vecteur flux hydraulique du composant ij U2 04 05 L int grale des flux sur une surface est effectu e dans POST _ ELEM par int gration de ce champ Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org co
10. 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 COUCHE et NIVE COUCHE partir des contraintes aux points d int gration de chaque couche 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 fini 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 galeme
11. 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 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
12. 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 Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Cod e AS ter Version 10 Titre Op rateur CALC ELEM Date 23 10 2012 Page 22 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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 fini 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
13. 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 E F F DISS E NO GA 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 Cod e AS ter Version 10 Titre Op rateur CALC ELEM Date 23 10 2012 Page 18 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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 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 EL
14. 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 AS ter Version 10 Titre Op rateur CALC_ELEM Date 23 10 2012 Page 17 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 DEV SUR CRIT Rapport entre le seuil d viatoire et le seuil d viatorique critique DIST 1SO Distance normalis e au 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 cont
15. 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 AS ter Version 10 Titre Op rateur CALC _ELEM Date 23 10 2012 Page 3 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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 maille
16. 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 AS ter Version 10 Titre Op rateur CALC_ELEM Date 23 10 2012 Page 21 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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 l PMPB ELGA l 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 pur X S purs LR avec M M M Ceci correspond la valeur maximum de SIXX dans une section circulaire R3 08 01 HA PMPB ELGA Valeurs de PM et PMPB aux points de Gauss ca
17. 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 e Pour l indicateur d erreur quatre choix sont possibles ERME_ELEM pour l indicateur en r sidus 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
18. 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 AN naine tent
19. u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Version 10 Code Aster Titre Op rateur CALC_ELEM Responsable Aimery ASSIRE Date 23 10 2012 Page 25 27 Cl U4 81 01 R vision 9958 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 l 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 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 al 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
20. A6 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 ERME ELNO Estimateur d erreur en r sidu calcul aux n uds R4 10 021 ERRE ELGA EL EM 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 Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code AS ter Version 10 Titre Op rateur CALC_ELEM Date 23 10 2012 Page 23 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 Conseil d utilisation des options 01Z1 ELEM Q1Z2 ELEM 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 01Z2 ELEM est le m me que pour les options ERZ1 EL
21. Code Aster Version 10 Titre Op rateur CALC_ELEM Date 23 10 2012 Page 1 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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 Cod e AS ter Version 10 Titre Op rateur CALC ELEM Date 23 10 2012 Page 2 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 Table des Mati res TE leonina a intimiste rates ten ie 1 PARA Ale et 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 tte ee atse tn
22. EM 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 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 sslv113c et d
23. 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 EPMG 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 html Code AS ter Version 10 Titre Op rateur CALC_ELEM Date 23 10 2012 Page 15 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 Pour les l ments de structure ces d formations sont
24. NF 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 10 Titre Op rateur CALC_ELEM Responsable Aimery ASSIRE Date 23 10 2012 Page 11 27 Cl U4 81 01 R vision 9958 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 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
25. NO 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 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 utilisati
26. N_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 Cod e AS ter Version 10 Titre Op rateur CALC ELEM Date 23 10 2012 Page 20 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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 es partir des champs EPSI ELNO ou EPME ELNO INVA 2 d formation quivalente de Von Mises INVA 256 d formation quivalente de
27. e Calcul des champs de contraintes aux n uds des faces des l ments 3D Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Version 10 Code Aster Titre Op rateur CALC ELEM Responsable Aimery ASSIRE Date 23 10 2012 Page 14 27 Cl U4 81 01 R vision 9958 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_FP 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
28. e 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 _ a 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 10 Titre Op rateur CALC ELEM Responsable Aimery ASSIRE EXTR ELNO EXTR Date 23 10 2012 Page 16 27 Cl U4 81 01 R vision 9958 ELGA Extraction des variables internes en THM uniquement respectivement aux n uds par l ments et aux points de Gau
29. es cf 2 5 5 TYPE OPTION CRETE OPTION SIEQ ELNO SIEQ ELGA EPEQ ELNO EPEQ ELGA EPMO ELNO EPMO ELGA ENDO ELGA ENDO ELNO CRIT ELNO 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 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 Cod e AS ter Version 10 Titre Op rateur CALC ELEM Date 23 10 2012 Page 6 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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 OPTION SIGM MASSIF OPTION SIEF E NO PROJ ELEM options de calcul des contraintes l ments de structure poutres
30. lcul 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 f Le r sultat est rang au num ro d ordre associ l instant t lol llel loll Par d faut le calcul se fait pour les num ros d ordre 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 t dans la liste des instants ayant servi faire le calcul non lin aire ERZL ELEM respectivement ERZ2 ELEM 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
31. logies 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 01 l 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 12 2 7 Op randes pour les options 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 TITRE en titre Titre que l on veut donner au r sultat de la commande U4 02 01 Exemples 3 1 Calcul du flux pour un evol ther evoth CALC ELEM reuse evoth RESULTAT evoth TOUT ORDRE OUI 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 10 Titre Op rateur CALC_ELEM Responsable Aimery ASSIRE OPTION FLUX ELNO Date 23 10 2012 Page 27 27 Cl U4 81 01 R vision 9958 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 poin
32. nt fini n est connect aucune singularit et vaut l ordre de la singularit si l l ment fini est Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Co d e A ster Version 10 Titre Op rateur CALC_ELEM Date 23 10 2012 Page 24 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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 correspond 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 l occurrence HOMARD pour Code Aster pour que celui ci construise le nouveau maillage en
33. nt 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 Cod e AS ter Version 10 Titre Op rateur CALC ELEM Date 23 10 2012 Page 12 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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 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 op
34. 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 gales 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 command
35. on Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code AS ter Version 10 Titre Op rateur CALC_ELEM Date 23 10 2012 Page 19 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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 quivalente 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 PRI
36. osante 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 Ces champs sont valu s partir d un champ de contraintes aux n uds par l ment option SIGM ELNO dans le cas lin aire ou option SIEF ELNO dans le cas non lin aire e calcul sur les mailles volumiques MODELISATION 3D ou 3D SI de la fa on suivant
37. oupe 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 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 I
38. 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 forces 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 QUADA4 8 9 En 3D idem avec FACE3 4 6 8 9 TETRA4 10 PENT
39. pyleft fdl html Co d e A ster Version 10 Titre Op rateur CALC_ELEM Date 23 10 2012 Page 26 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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 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 E le calcul des flux aux n uds via FLUX ELNO Le mot de 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 ana
40. ract 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 totale 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
41. re 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 Cod e AS ter Version 10 Titre Op rateur CALC_ELEM Date 23 10 2012 Page 10 27 Responsable Aimery ASSIRE Cl U4 81 01 R vision 9958 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 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 gr
42. ss 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 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
43. tes 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 erreur 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 enco
44. tion 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 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code AS ter Version 10 Titre Op rateur CALC_ELEM Date
45. ts de GAUSS pour un calcul thermo m canique evolas CALC ELEM reuse evolas RESULTAT evolas TOUT ORDRE OUI OPTION SIEF ELGA 3 4 Calcul des nergies potentielles pour un mode propre mode CALC ELEM reuse mode RESULTAT mode NUME MODE 3 OPTION EPOT_ELEM 3 5 Calcul de l endommagement de Lema tre ou de Lema tre Sermage evolas CALC ELEM reuse evolas OPTION ENDO ELGA ENDO ELNO RESULTAT evolas Manuel d utilisation Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Fascicule u4 81 Outils g n raux
Download Pdf Manuals
Related Search