[U4.81.01] Opérateur CALC_ELEM
Contents
1. TYPE OPTION SIGM_MASSIF OPTION SIEF ELNO ELGA SIGM ELNO DEPL SIEF ELGA DEPL ARCO ELNO SIGM options de calcul des contraintes l ments de structure poutres tuyaux coques ct S3 6 11 TYPE OPTION SIGM STRUCT OPTION SIEF ELNO ELGA SIGM ELNO DEPL SIEF ELGA DEPL SIGM ELNO TUYO SIPO ELNO DEPL EFGE ELNO DEPL EFGE ELNO CART SIGM_ELNO CART SIGM_ELNO SIEF SIPO ELNO SIEF options de calcul des d formations cf 3 6 2 TYPE OPTION EPSI OPTION EPSI ELNO DEPL EPSI ELGA DEPL EPME ELNO DEPL EPME ELGA DEPL DEGE_ELNO_DEPL EPSI ELNO TUYO EPVC_ELNO EPVC_ELGA options de calcul d nergies cf S3 6 4 TYPE OPTION MENER OPTION EPOT_ELEM_DEPL ECIN_ELEM_DEPL ENE E GA ENEL ELNO ELGA ETOT ELGA ETOT ELNO ELGA DISS ELGA DISS ELNO ELGA ETOT _ELEM 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 5 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision2. 2 3 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 OUl Toutes les mailles porteuses d l ments finis seront trait es C est la valeur par d faut GROUP MA 1 grma MAILLE 1l maille Seules les mailles incluses dans 1_grma et ou 1 maille seront trait es 2 4 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 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 TOQUI Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http
3. ia M evol_noli VMIS DEFAUT TOTAL VMIS_CINE TOTAL CIN Ra autres options cf S3 6 7 TYP E OPTION AUTRES OPTION a VALE NCOU MAXI 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 8 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 NOM CHAM ch cham elem NOM CMP cmp TXM options thermiques OPTION FLUX ELNO TEMP FLUX ELGA TEMP DEUL ELGA TEMP DETE_ELNO DLTE ERRE ELEM TEMP ERRE ELNO ELEM SOUR ELGA ELEC 4 DURT_ELGA_META d DURT_ELNO META A i HYDR ELNO ELGA r RESULTAT resu evol_ther 0 SENSIBILITE 1l parasensi theta theta_geom listpara para_sensi options acoustiques OPTION PRES ELNO DBEL PRES_ELNO_REEL PRES_ELNO_IMAG INTE_ELNO ACTI INTE_ELNO REAC RESULTAT resu acou harmo mode_acou TITRE titre 1 Kn INFO 1 DEFAUT 2 SENSIBILITE l parasensi
4. 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 Moyenne d une composante du d placement Moyenne d une composante du tenseur des contraintes 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 24 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 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
5. 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 l DEUL ELGA TEMP l DETE ELNO DLTE D riv e Eul rienne du champ de temp rature aux points de Gauss ou aux n uds R4 03 01 N cessite la connaissance de la d riv e Lagrangienne des temp ratures donc d avoir activ l option SENSIBILITE dans THER LINEAIRE et d utiliser le mot cl SENSIBILITE dans CALC ELEM DURT_ELGA META DURT ELNO META Calcul de duret aux points de Gauss ou aux n uds partir des phases m tallurgiques cf R4 04 01 HYDR ELNO ELGA Calcul de l hydratation 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 8 Op randes pour les options acoustiques 2 8 1 Op rande OPTION PRES ELNO DBEL Calcul de la pression aux n uds en d cibels PRES ELNO REEL Calcul des parties r elles du champ de pression aux n uds PRES ELNO IMAG Calcul des parties imaginaires r elles du champ de pression aux n uds INTE ELNO ACTI Calcul de l intensit acoustique active aux n uds INTE ELNO RE
6. theta theta_geom listpara para_sensi 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 9 Titre Op rateur CALC_ELEM Date 05 01 2010 Page 9 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 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 504 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 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 qu
7. 2102 options de calcul de crit res cf S3 6 5 TYP E OPTION GRITS gt OPTION E UI ELNO SIGM UI ELGA SIGM UI ELNO EPSI UI ELGA UI ELNO EPME UI ELGA EPME ENDO ELNO SIGA ENDO ELNO SINO ENDO ELGA ENDO ELN ELGA TUYO TUYO RUPT O SIEQ ELNO O O EPEQ ELN CRIT_ELN options de calcul d indicateurs d erreur cf S3 6 6 TYP OPTION INDI ERRE OPTION SIGM_NOZ1_ELGA ERRE ELEM NOZ1 SIGM NOZ2 ELGA ERRE ELEM NOZ2 SIRE ELNO DEPL ERRE ELGA NORI ERRE ELNO ELG zal D autres options cf S3 6 7 TYPE OPTION YAUTRES OPTION VALE NCOU MAXI NOM CHAM ch cham elem NOM CMP cmp TXM PRES DBEL DEPL VNOR_ELEM DEPL options de calcul de d riv es li es la sensibilit cf S3 2 TYPE OPTION DERIVEES OPTION DEU ELGA _ DEPL i DEDE ELNO DLDE DEST E NO DLSI 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 6 29 Responsable Aimery ASSIRE
8. GAMP PCR SEU IND HYD PCOHI COMP ROC SEUI ANG DEV nri x22 n33 x12 X13 x23 DIST_DEV DEV SUR CRIT DIST_ISO NB_ITER ARRET NB REDE SIGNE Remarque Date 05 01 2010 Page 17 29 Cl U4 81 01 R vision 2102 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 Rapport entre le seuil d viatoire et le seuil d viatorique critique Distance normalis e au seuil isotrope Nombre d it rations internes Valeur du test local d arr t du processus it ratif Nombre de red coupage local du pas de temps Signe du produit contract de la contrainte d viatorique par la d
9. L emploi des mots cl s TOUT ORDRE Date 05 01 2010 Page 12 29 Cl U4 81 01 R vision 2102 est d crit dans le document NUM ORDRE INST FREQ U4 71 00 2 6 2 6 1 Op randes pour les options m caniques Option de calcul des contraintes EPL SIEF ELGA D Calcul de l tat de contrainte par l ment aux points d int gration de l l ment points de GAUSS ou points 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 SIEF ELNO ELGA Calcul de l tat de contrainte aux n uds par l ment partir de l tat de contrainte aux points de Gauss EPL SIGM_ELNO D Calcul des contraintes par l ment aux n uds partir des d placements lasticit lin aire SIGM ELNO COQU NUME COUCHE i DEFAUT NIV COUCHE 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 GREEN GR cette opti
10. formation plastique d viatorique 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 R8VIDE 2 6 4 Options de calcul d nergie ECIN ELEM DEPL nergie cin tique d un l ment ENEL ELNO ELGA 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 DEPL qui calcule l nergie d 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 nergi lastique qui s crit Er 50 o 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 9 Titre Op rateur CALC_ELEM Date 05 01 2010 Page 18 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 Ce calcul s appuie sur le champ de contraintes aux points de Gauss obtenu par SIEF ELGA OUSIEF ELGA DEPL l DISS_ELNO ELGA 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 do
11. 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 4 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 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 S3 6 1
12. Date 05 01 2010 Page 23 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 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 ECA 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 II 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 S
13. 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 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 EFGE ELNO DEPL De plus cette option n est utilisable que pour les DKT DST 04G GRILLE 0 ANGLE delta I Ao Oes DEFAUT edelta angle en degr s valeur enti re compt partir de la position de la g n ratrice de l l ment tuyau ECTEUR x y z 1 R pard faut VECTEUR 1 0 0 VV ANGL REP a B 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 ELEM U4 42 01 2 5 S lection des num ros d ordre 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 9 Titre Op rateur CALC_ELEM Responsable Aimery ASSIRE
14. 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 ELEM Estimateur d erreur en quantit s d int r t bas s sur les r sidus calcul aux n uds RADI ELGA SIGM Indicateur de perte de radialit aux points de Gauss R4 20 01 RADI_ ELNO SIGM Indicateur de perte de radialit aux n uds R4 20 01 SIGM_NOZ1 ELGA 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 SIGM NOZ2 ELGA 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 01 Remarque SING ELEM Pour les options RADI ELNO SIGM et RADI ELGA SIGM 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 L ERA ve lool 0 L indicateur de perte de radialit est calcul par IP 1 lo 01 0 Par d faut le calcul se fait pour les num ros d ordre 1 n 1 Si on pr cise la liste d instant avec des trous ventuellement le c
15. fdl html Cod e AS ter Version 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 10 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 Ce mot cl est suivi d une liste de param tres sensibles Il pr cise que l on ne s int resse pas au r sultat en lui m me mais la d riv e du r sultat par rapport un param tre Ainsi une s quence du type RESULTAT resu ENSIBILITE ps OPTION SIEF ELGA _ DEPL uv Signifie que l on veut calculer aux points de Gauss la d riv e des contraintes par rapport au param tre ps Voir U4 50 02 pour les d tails sur les param tres associ s aux mots cl DEUL ELGA DEPL DEDE ELNO DLDE D riv e Eul rienne du champ de d placements aux points de Gauss ou aux n uds R4 03 01 disponible en lin aire seulement N cessite la connaissance de la d riv e Lagrangienne des d placements donc d avoir activ l option SENSIBILITE dans MECA STATIQUE et d utiliser le mot cl SENSIBILITE dans CALC ELEM DESI ELNO DLSI D riv e Eul rienne du champ de contraintes aux n uds R4 03 01 disponible en lin aire seulement N cessite la connaissance de la d riv e Lagrangienne des contraintes en lasticit lin aire donc d avoir activ l option SENSIBILIT dans MECA STATIQUE et d utiliser le mot cl SENSIBILITE dans CALC ELEM
16. 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 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 e 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 e la nouvelle taille des l ments finis est donn e un remailleur en l occurrence HOMARD pour le 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 calcu
17. www gnu org copyleft fdl html Cod e AS ter Version 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 11 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 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 T 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 aa Pour la couche nume d finie par NUME COUCH effectuer le calcul l mentaire permet de pr ciser l ordonn e o l on d sire 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 MAIL DEFAUT MOY INE SUP e MAIL plan du maillage MOY plan moyen e
18. AC Calcul de l intensit acoustique r active aux n uds Les d finitions se trouvent dans R4 02 01 2 9 Op rande TITRE 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 29 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 0 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 OUI OPTION FLUX ELNO TEMP 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 SIGM 3 3 Contraintes aux points de GAUSS pour un calcul thermo m canique evolas CALC ELEM reuse evolas RESULTAT evolas TOUT ORDRE OUL OPTION SIEF ELGA DEPL 3 4 Calcul des nergies potentielles pour un mode propre mode CALC ELEM reuse mode RESULTAT mode NUME MODE 3 OPTION EPOT ELEM DEPL 3 5 Calcul de la d riv e des contraintes evolas CALC ELEM reuse evolas RESULTAT evolas
19. Cl U4 81 01 R vision 2102 options pour les r sultats non lin aires produits par STAT NON LINE ou DYNA NON LINE RESULTAT resu evol noli Q TYPE OPTION TOUTES DEFAUT OPTION toutes les options ci dessous options de calcul des contraintes l ments d milieux continus 2D et 3D cf S3 6 1 TYPE OPTION SIGM_MASSIF OPTION SIEF ELNO ELGA ARCO ELNO SIGM options de calcul des contraintes l ments d structure poutres tuyaux coques cf S3 6 1 TYPE OPTION SIGM STRUCT OPTION SIEF ELNO ELGA EFGE ELNO CART SIGM ELNO TUYO SIGM ELNO COQU SIGM ELNO SIEF SIPO ELNO SIEF options de calcul des d formations cf S3 6 2 TYPE OPTION EPSI OPTION EPST ELNO DEPL EPSI ELGA DEPL EPSG E iNO DEP jl EPSG ELGA DEPL EPME ELNO DEPL EPME ELGA DEPL EPMG_ELNO_ DEPL EPMG ELGA DEPL EPSP ELNO EPSP_ELGA EPFD_ELNO EPFD_ELGA EPFP ELNO EPFP ELGA EPVC_ELNO EPVC_ELGA EPSI ELNO TUYO DEGE ELNO DEPL options d interpolation et d extraction des variables internes cf S3 6 3 TYP E OPTION VARI OPTION VARI ELN
20. Cod e AS ter Version 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 1 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 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 DEPL FLUX ELGA TEMP VARI ELNO ELGA 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 2 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 Table des mati res 2 1 Op randes RESULTAT MODELE CHAM MATER CARA _ELEM EXCIT SOLVEUR 9 211LOpP rANdes RESULTAT Es end tdi ditnde dede donr diam 9 2 1 2 Op randes MODELE CHAM _MATER CARA ELEM ni anneeeennnenmenns 9 2 1 3 Mote EXC IT 5888 drone gg ere vnton ragnpo
21. EG2 3 TRIA3 6 QUADA4 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 enqu rit tort d une ventuelle CL ERRE ELNO ELEM Estimateur d erreur en r sidu calcul aux n uds R4 10 02 ERRE ELGA ELEM Estimateur d erreur en r sidu calcul aux points de Gauss R4 10 021 QIZ1 ELEM SIGM respectivement OIZ2 ELEM SIGM Estimateur d erreur en quantit s d int r t bas sur la m thode de Zhu Zienkiewicz lasticit lin aire 2D QIRE ELEM SIGM Estimateur d erreur en quantit s d int r t bas sur les r sidus en m canique calcul par l ment Conseil d utilisation des options QIZ1 ELEM SIGM QIZ2 ELEM SIGM QIRE ELEM _ SIGM Le domaine d utilisation des options QIZ1_ ELEM SIGM et OIZ2 ELEM SIGM est le m me que pour les options ERZ1 ELEM SIGM et ERZ2 ELEM SIGM et celui de l option OIRE ELEM SIGM est le m me que celui de l option ERRE ELEM SIGM en m canique
22. O ELGA VARI ELNO TUYO VARI ELNO COQU T EXTR ELGA VARI EXTR ELNO VARI 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 7 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 options de calcul d nergies cf S3 6 4 TY P E OPTION MENER OPTION options de calcul de crit res cf S3 6 5 TYP E OPTION CRET OPTION T Z 2O00O0O0OTO T T Q T an z Q GOG G GGG PRE EE pe T P PME IGA INO UPT Q F T J O Z T el O Z sS S R Q Del H J 7 z T 2 J O a ELGA EF EF GA TUYO TUYO l lt J U w wW m HE Q Z l I I l J O EAE 5 EE 0O00O OOx O OOpPODP S S H O H W HE D HZ H HO mt H O 224 Q 3 options de calcul d indicateurs d erreur cf 53 6 6 TYPI E OPTION INDI ERRE OPTION ERRE_ ELEM SIGM ERZ1 ELEM SIGM ERZ2 ELEM SIGM QIZ1 ELEM SIGM
23. SENSIBILITE psl ps2 OPTION SIEF ELGA DEPL 3 6 Calcul de l endommagement de Lema tre ou de Lema tre Sermage evolas CALC ELEM reuse evolas OPTION ENDO ELGA Fu ENDO ELNO ELGA E RESULTAT evolas Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html
24. Valeurs de PM et PMPB aux points de Gauss calcul es partir de SIEF ELGA PMPB ELNO SIEF valeurs de PM et PMPB aux n uds calcul es partir de SIEF ELNO ELGA En toute rigueur ces crit res sont appliquer aux contraintes primaires Cette distinction est faire par l utilisateur 2 6 6 Options de calcul d indicateurs d erreur DCHA ELGA SIGM Indicateur local de d charge aux points de Gauss R4 20 01 DCHA ELNO SIGM Indicateur local de d charge aux n uds R4 20 01 Remarque Pour les options DCHA ELGA SIGM et DCHA ELNO SIGM 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 lol llo lo 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 ERZ1 ELEM SIGM respectivement ERZ2 ELEM SIGM Calcul de l estimateur d erreur de ZHU ZIENKIEWICZ lasticit lin aire 2D partir de l option SIGM NOZ1 ELGA respectivement SIGM NOZ2 ELGA Si ce dernier champ n existe pas dans resu il est automatiquement const
25. alcul 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 Minna 28 3 9 Calcul dela d riv e des COntraINtes s seu miunsestansnnermmentounesapsutannenienthennmenmonnadesaes 28 3 6 Calcul de l endommagement de Lema tre ou de Lema tre Sermage ss 28 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 9 Titre Op rateur CALC_ELEM Date 05 01 2010 Page 3 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 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 MAILLE
26. alcul 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 PREC ERR err R O TYPE ESTI ERRE_ ELEM SIGM ERZ1 ELEM SIGM ERZ2 ELEM SIGM QIRE_ ELEM SIGM QIZ1 ELEM SIGM QIZ2 ELEM SIGM 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 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 9 Titre Op rateur CALC_ELEM Date 05 01 2010 Page 25 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 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
27. amps ont t cr s Le principe de ces champs est d extraire du champ VARI ELGA ou VARI ELNO ELGA pour le cham_elem aux n uds la variable interne qui nous int resse via un mot cl plus parlant que V1 V2 NOM VARI nom vari TXM Le nom des nouveaux champs est EXTR ELGA VARI et EXTR ELNO VARI pour les cham_elem et EXTR NOEU VARI pour le cham_no En tant que post traitement ces champs sont calcul s par CALC EI syntaxe utiliser est la suivante EM et CALC NO La 1 pour un cham _elem GAMP CALC ELEM RESULTAT UI OPTION EXTR ELNO VARI NOM VARI GAMP gt nouveau mot cl pour indiquer quelle variable on souhaite extraire via un nom cod 1 pour un cham_ no GAMP CALC NO reuse GAMP RESULTAT GAME OPTION EXTR NOEU VARI 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 variation de la porosit du mat riau 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 9 Titre Op rateur CALC_ELEM Responsable Aimery ASSIRE DRHOLQ DPVP SATLIQ MEVP IND ETA LL D IND END TEMP MAX
28. c R de la loi de comportement EPOT_ELEM_DEPL 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 1 th 1 th th EPOT gt f E U AE U dy f Ee U AE U dv 5 f E U A U dv element element element e pour les l ments de poutres EPOT U K U U B Ae 7 s 4e 2 e et pour les l ments de plaques et coques EPOT ZU K U U B As 2 6 5 Options de calcul de crit res CRIT_ELNO_RUPT 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_DEPL et mots cl s NUME COUCHE et NIVE COUCHE et des contraintes limites fournies sous ELAS ORTH dans DEFI _ MATERIAU le champ CRIT ELNO RUPT contient 6 composantes 0 mes CRIL crit re de rupture en traction dans le sens L si o gt 0 Xr Or ax CRILP crit re de rupture en compression dans le sens L si o lt 0 Xe Or ix CRIT crit re de rupture en traction dans le sens T si o gt 0 Y Or 7 i CRITP crit re de rupture en compression dans le sens T si o lt 0 Yo 2 Alal Eae CRILT crit re de rupture en cisaillement S LT et CRITH crit re de Tsai Hill voir exemple dans test SSLS121 V3 03 121 Toutes ces
29. e 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 encore traiter STOP_SINGULIER NON 2 2 Op rande SENSIBILITE SENSIBILITE 7 theta theta geom listpara para_sensi Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft
30. e 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 6 8 Op rande NORME Pour les options DCHA ou RADI_ on doit choisir une norme pour le tenseur des contraintes R4 20 01 le choix est fait gr ce au mot cl NORME VMIS DEFAUT On prend le second invariant du d viateur du tenseur des contraintes o VMIS_CINE Avant de prendre le second invariant du d viateur on retire au tenseur des contraintes le tenseur X caract risant l tat d crouissage cin matique TOTAL On prend le second invariant du tenseur en tenant compte de la trace TOTAL CINI za Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Version 9 Code Aster Titre Op rateur CALC_ELEM Responsable Aimery ASSIRE Date 05 01 2010 Page 27 29 Cl U4 81 01 R vision 2102 On retire d abord X avant de prendre le second invariant du tenseur total des contraintes tenant compte de la trace Remarque Les normes VMIS CINE et TOTAL CINE ne sont utilisables que si le calcul lastoplastique a t fait avec le c
31. entes aux n uds VMIS contrainte de von Mises VMIS SG contrainte de von Mises sign e par la trace de g PRIN_1 PRIN 2 PRIN 3 contraintes principales TRESCA contrainte de Tresca Pour les l ments de milieux continus 2D et 3D elles sont partir de la version 7 2 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 DEPL 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 EQUI NOEU SIGM 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 DEPL en lin aire et SIGM ELNO COQU en non lin aire INDI LOCA 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 INDI LOCA ELGA contient les composant
32. es suivantes INDICE Indicateur de localisation valant 0 si det N H N gt O et valant 1 sinon ce qui correspond a l initiation de la 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 EPEQ ELNO TUYO SIEQ ELNO TUYO Calcul des d formations g n ralis es et des contraintes pour un l ments tuyau Ce sont des valeurs quivalentes de type EQUI ELGA SIGM et EQUI FELGA EPSI en un point de la section C est une extraction effectu e suivant le m me principe que l option d j existante SIGM ELNO TUYO Le calcul des d formations s effectue dans une couche et un secteur angulaire d l ments tuyau PMPB_ELGA_SIEF l PMPB_ELNO SIEF 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 S N M R pups el 7 avec M M M 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 22 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 Ceci correspond la valeur maximum de SIXX dans une section circulaire R3 08 01 PMPB ELGA SIEF
33. l 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 e _ ERRE ELEM SIGM 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_ SIGM 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 ERRE_ELEM SIGM est choisi par d faut message d alarme mis Si les deux indicateurs de Zhu Zienkiewicz sont pr sents on choisit ERZ1 ELEM SIGM 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 DEPL qui est l nergie potentielle de d formation lastique
34. n de fluage propre de Granger 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 16 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 EPSP ELNO D formations an lastiques obtenues par extrapolation aux n uds des valeurs aux points de Gauss cf EPSP ELGA 2 6 3 Options d interpolation et d extraction des variables internes VARI ELNO ELGA 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 VARI ELNO TUYO 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 VART ELNO COQU ES Calcul des variables internes dans une couche d l ments coque d finie par NUME COUCH et NIVE COUCHE E EXTR_ELNO VARI EXTR ELGA VARI 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 ch
35. nes er daurs ee trs ate totem A 9 2 114MOt cl SOLVE UR sen hrs datest ermraimmncheme annee rentree entame 9 2 2 Op rand SENSIBILIT ccris AAAA EEEE 9 2 3 S lection des mailles concern es par le calcul 10 2 4 D finition d un rep re local de d pouillement op rande REPE_COQUE iii 10 2 9 Selection des NUM rOS T Ordr e rsicissiii inatia ironiaa near de a aa None aa a iaaa aaia 11 2 6 Op randes pour les options M CANIQUES iiiirannnaanaannnnnnnnmnnennennnnss 11 26 1 Option d ecalc ldes Contrai S misrin aiaiai 11 2 6 2 Options de calcul des d fOrMAtIONS 242hnh00 6 uonanirsss sens sarsnameete ne drenaarate mare 14 2 6 3 Options d interpolation et d extraction des variables internes 15 2 0 4 Optone de calkuldietetg e HR sd ie teen 17 2 6 5 Options de calc lde CHHST S 2080 2422 nitro sndane semer mme ada ame iaia 18 20 L AUO SODAONS a aaa te states 2e ann DE eai aa tin 25 2 6 8 Op rand NORME 28500200008rm084uamesanebacsbnrshunnendurs sde ladensenumelateies saaunespiedentatennii 25 2 7 Op randes pour les options thermiques 27 2 1 HOp rand e OPH NESESS SN en enr ne ni toi tre aniame screen 27 2 8 Op randes pour les Options ACOUSTIQUES 222 2 22808 222 sua rai ina datane inerte a samenta ta ad aa 27 2 8 1 Op rande OPTION 8eme inner nnuade tetes lient shine te tentent data 27 29 Opetande TRE asiingie Aaa aaa 27 3 1 Calcul du fux pour un evol ANET sikioni aa iaaea ai aa a aaa daa 28 3 2 C
36. omportement VMIS CINE LINE 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 9 Titre Op rateur CALC_ELEM Date 05 01 2010 Page 28 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 2 7 Op randes pour les options thermiques 2 7 1 Op rande OPTION FLUX ELGA TEMP Calcul des flux de chaleur aux points d int gration de GAUSS partir de la temp rature FLUX ELNO TEMP Calcul des flux de chaleur aux n uds partir de la temp rature ERRE ELEM TEMP TERRE BE NO E EM Estimateurs d erreur en r sidu en thermique R4 10 03 Il faut pr alablement effectuer dans CALC ELEM le calcul des flux aux n uds via FLUX ELNO TEMP Le mot cle INFO procure tous les affichages interm diaires du calculs connectivit s normales diam tres valeurs des champs jacobien L option ERRE ELNO ELEM permet de ramener le champ par l ment ERRE ELEM TEMP a un champ aux n uds par r l ment ce qui permet de faire des relev s de valeurs ou des impressions Visualisations SOUR ELGA ELEC 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 021
37. on 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 SIGM ELNO TUYO 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 SIGM EFGE ELNO CART ELNO CART 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 EPL 0 REP EFGE ELNO D E COQUE 0 PLAN MAIL DEFAUT MOY YINF 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 Manuel d utilisation Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copylef
38. ons 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 1 th EPMG ELNO DEPL EPMG ELGA DEPL 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 E u 5 ui U it TETE EPSG ELGA DEPL D formations de Green Lagrange aux points de Gauss 1 Ep u 5l i u tu itk j EPSG ELNO DEPL D formations de Green Lagrange aux n uds EPSI ELNO DEPL EPSI ELGA DEPL Calcul des d formations par l ment aux n uds ou aux points de Gauss partir des d placements 1 Ey u 5 U 0 EPST ELNO TUYO 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 He et de variables internes on calcule chaque instant e elu A o e T e e o est la d formatio
39. ontrainte de von Mises VMIS 4 Zs S avec s 0y atr 0 6 IS SG contrainte de von Mises sign e par la trace de g INVA_2 second invariant de INVA 256 second invariant de sign par la trace de PRIN_1 PRIN 2 PRIN 3 contraintes principales ESCA contrainte de Tresca ECT_ 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 D 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 EQUI ELNO EPSI EQUI ELNO EPME D formations quivalentes aux n uds calcul es partir des champs EPSI ELNO DEPL OUEPME ELNO DEPL INVA 2 second invariant de INVA 2SG second invariant de sign par la trace de lt 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 21 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 PRIN_1 PRIN 2 PRIN 3 d formations principales EQUI ELNO SIGM Contraintes quival
40. quantit s sont calcul es dans le rep re d orthotropie de la coche consid r e ENDO ELNO SIGA ENDO ELNO SINO 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 19 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 Calcul du taux de triaxialit et de la contrainte quivalente d endommagement aux n uds _SINO ou aux points de Gauss _ SIGA partir des contraintes Soit s le d viateur du tenseur des contraintes s o itr o 3 Oy 55s o tala 3 Le taux de triaxialit est d fini par La contrainte quivalente d endommagement est o 0 2 1 v 3 1 2v 06 ENDO ELGA Calcul du dommage qd 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 YP a z p sl P Di 2 go Y z 2E 1 D o S et s sont des coefficients caract ristiques du mat riau et p 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 Calcul syst matique du taux de triaxialit et de la contrainte quivalente d endommagement o SI_ENDO CRETE Tr TRIAX O Oog 1 s o 3
41. quation valable uniquement en lasticit L utilisation de cette option en lastoplasticit est donc manipuler avec prudence SING ELNO ELEM 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 ELEM n est pas calcul e 2 6 7 Autres options VNOR ELEM DEPL 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 l VALE_ NCOU MAXI NOM CHAM ch eham 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 DEPL les champs de d formations EPSI ELGA DEPL les champs de valeurs quivalentes EQUI_ELGA SIGM EQUI ELGA EPSI les champs de variables internes VARI_ELGA Le champ cr e de nom VALE NCOU MAXI contient pour chaque instant les composantes E MIN valeur minimum MAX valeur maximum NCOUMIN num ro de la couche pour la valeur min NCOUMAX num ro d
42. ruit au pr alable voir R4 10 01 ERRE t HA EM SIGM 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 ERRE ELEM SIGM 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 OUl 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 DEPL OU SIRE ELNO DEPL en lin aire par SIEF ELNO ELGA 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 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 9 Titre Op rateur CALC_ELEM
43. s vo tes Le but de cette option est d estimer le champ de contraintes sur les parements amont et aval de l ouvrage surface alors que la structure est mod lis e en volumique Ce champ est valu partir d un champ de contraintes aux n uds par l ment option SIGM_ELNO_DEPL dans le cas lin aire ou option SIEF ELNO ELGA dans le cas non lin aire calcul sur les mailles volumiques MODELISATION 3D ou 3D SI Les composantes du nouveau tenseur des contraintes sont 1 SI_NORM La composante suivant la normale la face 2 SI ARC La composante de la contrainte d arc elle est perpendiculaire la normale de l l ment et parall le l horizontale direction fix e en imposant conventionnellement la troisi me ascendante 3 SI CONS La composante de la contrainte de console c est la composante suivant la troisi me direction dans le plan tangent et donc ascendante Pour calculer ces contraintes on effectue un changement de rep re du rep re global XYZ au rep re local X1X2X3 lors du traitement de chaque n ud de la maille de peau Le rep re local en un n ud d une maille est d fini par 1 X3 vecteur normal 2 X1 un vecteur horizontal calcul partir d une direction verticale e d finie soit par ses composantes cart siennes l aide du mot cl VECTEUR soit par ses angles nautiques l aide du mot cl ANGL REP dans l option CALC ELEM cf Figure sui
44. sur un l ment fini et int gr e partir des d placements et de la temp rature valable uniquement pour un comportement lastique T 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 Fascicule u4 81 Outils g n raux Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version 9 Titre Op rateur CALC_ELEM Date 05 01 2010 Page 26 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 En toute rigueur le calcul de l ordre de la singularit est obtenu partir de l nergie th orique en pointe de fissure
45. t fdl html Code Aster Version 9 Titre Op rateur CALC_ELEM Responsable Aimery ASSIRE Date 05 01 2010 Page 13 29 Cl U4 81 01 R vision 2102 Dans le cas des mod lisations de plaques avec excentrement DKT DST Q4G GRILLI PLAN permet de d finir le plan de calcul ma x 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 SIPO ELNO DEPL Contraintes dans la section de poutre d compos e en contributions de chaque effort g n ralis SN uns due l effort normal MY z SMFY O x7 7 due au moment de flexion MY Rd MZ y SMFZ O x7 due au moment MZ L Vya BE n SVY O y F due l effort tranchant Vy a coefficient de cisaillement dans la direction y Feu SVZ Gi 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 o dues aux deux moments de flexion sont les valeurs maximum de celles calcul es en Ymin Ymax d une part et en Zmin Zmax d autre part pour une section g n rale cf AFFE CARA ELEM U4 42 01 T Pour une section rectangulaire e on calcule la valeur de SMFY en z HZ 2 e on calcule la vale
46. tr o I avec Oog 2 1 o 3 o Calcul du dommage total par cumul lin aire 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 20 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 V V T VI D CUMULI CH DST D 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 1 valeur du dommage de Lema tre Sermage D_CUMULE valeur du dommage de Lema tre Sermage cumul el E e Es ENDO ELNO ELGA Endommagement de Lemaitre Sermage obtenues par extrapolation aux n uds des valeurs aux points de Gauss cf ENDO ELGA EQUI ELGA EPSI EQUI ELGA EPME D formations quivalentes aux points de Gauss calcul es partir des champs EPSI ELGA DEPL OUEPME ELGA DEPL INVA_2 second invariantde s nva 2 2 EE INVA 256 second invariant de sign par la trace de PRIN_1 PRIN_ 2 PRIN 3 d formations principales Pour les l ments TUYAU les composantes calcul s sont INVA_2 et INVA 25G EQUI ELGA SIGM Contraintes quivalentes aux points de Gauss IS c
47. ur de SMFZ en y HY 2 Pour une section circulaire on calcule les valeurs de SMFY et SMFZ pour y et z valant R SIGM E SIPO E T NO SIEF NO SIEF T Calcul des contraintes lin aris es par l ment aux n uds partir des efforts g n ralis s contenus dans le champ SIEF ELNO ELGA Les expressions sont les m mes que pour les options SIGM ELNO DEPL ou SIPO ELNO DEPL utilisables seulement en lasticit lin aire Ici les options SIGM ELNO SIEF et SIPO ELNO SIEF effectuent les m mes calculs que N M yR M yR S L L partir du champ SIEF ELNO ELGA qui peut tre calcul pour des comportements non lin aires Ces contraintes locales ne sont pas les contraintes r elles mais une estimation des contraintes dues aux efforts g n ralis s sous l hypoth se d une r partition lin aire dans la section de la poutre U N SIGM_ELNO_DE U OU SIPO T NO DEPL par exemple o XX 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 14 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 ARCO ELNO SIGM Calcul du champ de contraintes d arc et de console pour l analyse m canique des ouvrage
48. vante 3 X2 vecteur compl tant le tri dre v n Th e 12 calcul du tri dre local au n ud 1 Le vecteur horizontal h est d abord valu par produit vectoriel de n normale port e par le n ud consid r de l l ment par e v4 vecteur vertical compl tant le tri dre local se d duit alors automatiquement 2 6 2 Options de calcul des d formations DEGE ELNO DEPL Calcul des d formations g n ralis es 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 EPFP 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 9 Titre Op rateur CALC ELEM Date 05 01 2010 Page 15 29 Responsable Aimery ASSIRE Cl U4 81 01 R vision 2102 EPFP ELGA Calcul aux n uds ou aux points de Gauss des d formations de fluage propre associ es au mod le GRANGER FP o 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 EPME ELNO DEPL EPME ELGA DEPL Calcul aux n uds ou aux points de Gauss des d formati
Download Pdf Manuals
Related Search