[U4.81.06] Opérateur CALC_ERREUR
Contents
1. e Pour l nergie de d formation totale on utilise Manuel d utilisation Fascicule u4 81 Outils g n raux Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster pe es default Titre Op rateur CALC_ERREUR Date 21 02 2013 Page 9 9 Responsable Josselin DELMAS Cl U4 81 06 R vision 10471 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 O et 1 un message fatal est mis si cette condition n est pas v rifi e Le p rim tre d utilisation est le m me mais plus r duit que celui de l indicateur d erreur choisi savoir e Pour l indicateur en r sidu l ments finis des milieux2. sur les r sidus calcul s aux n uds SIZ1 NOEU Calcul des contraintes aux n uds lasticit lin aire 2D et 3D 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 Manuel d utilisation Fascicule u4 81 Outils g n raux Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Fi default Titre Op rateur CALC_ERREUR Date 21 02 2013 Page 8 9 Responsable Josselin DELMAS Cl U4 81 06 R vision 10471 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 SING ELEM PREC ERR err R TYPE ESTI ERME ELEM ERZ1 ELEM ERZ2 ELEM QIRE E EM QIZ1 ELEM QIZ2 ELEM Cette option R4 10 04 vise am liorer le traitement des singularit s dans les strat gies d adaptation de maillage en l occurrence avec HOMARD En pratique les indicateurs d erreur sont lev s dans les zones singuli res si bien que rapidement seules les zones singuli res sont raffin es et masquent donc les autres zones sensibles zones fort gradient que l on souhaiterait raffine
3. continus en 2D triangles et quadrangles ou 3D uniquement les t tra dres pour un comportement lastoplastique e Pour l indicateur de Zhu Zienkiewicz l ments finis des milieux continus en 2D triangles et quadrangles pour un comportement lastique En toute rigueur le calcul de l ordre de la singularit est obtenu partir de l nergie th orique en pointe de fissure quation valable uniquement en lasticit L utilisation de cette option en lastoplasticit est donc manipuler avec prudence l SING_ELNO D tection des singularit s et carte de modification de tailles aux n uds par l ment 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 4 Op randes pour les options thermiques 2 4 1 Op rande OPTION ERTH_ELEM ERTH_ELNO ERTH_NOEU Estimateurs d erreur en r sidu en thermique calcul s partir des flux aux n uds R4 10 03 S le calcul des flux aux n uds n est pas pr alablement effectu dans la commande CALC CHAMP via FLUX ELNO il est automatiquement r alis dans CALC ERREUR Le mot cl 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 fair
4. Code Aster Pat Titre Op rateur CALC_ ERREUR Date 21 02 2013 Page 1 9 Responsable Josselin DELMAS Cl U4 81 06 R vision 10471 Op rateur CALC ERREUR 1 But Cr er ou compl ter un r sultat en calculant des champs par l ment d erreur de discr tisation Manuel d utilisation Fascicule u4 81 Outils g n raux Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Version Aster default Titre Op rateur CALC ERREUR Date 21 02 2013 Page 2 9 Responsable Josselin DELMAS Cl U4 81 06 R vision 10471 Table des Mati res MBU bander non quote tnte initie tintin iles eine tete lineal ein inNtn tete nleniei site i ere i 1 PRE A le E A A ALE EEE E ee 3 2 1 Op randes RESULTAT MODELE CHAM_MATER CARA_ELEM EXCIT SOLVEUR 5 21 1 Op randes RESULTAT 2 252 u2tin aana oebpapene gagna oaa 5 2 1 2 Op randes MODELE CHAM_MATER CARA ELEM sise 5 2 1 3 Moteli EXCIT 5 su diarnre sponsors a a aaia 5 2 14 Motole SOLVEUR robe ruine setramenuponne teur se prise nee danse tatainetue ane eee 5 2 2 Selecion des NUMEROS MOLITE enr tete ndta near tite nee ta taie tenants tetes non 5 2 3 Op randes pour les options m Caniques MMM dnnnnnannnnnnnneneennnns 5 2 3 1 Options de calcul d indicateurs d erreur 6 2 4 1 Op r nd OPTION inerte a tet
5. DELE 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 0 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 0 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 chargeme
6. e des relev s de valeurs ou des impressions visualisations L option ERTH_NOEU permet de transformer le champ aux n uds par l ment en un champ aux n uds 2 5 Op rande TITRE 0 TITRE titre Titre que l on veut donner au r sultat de la commande U4 02 01 Manuel d utilisation Fascicule u4 81 Outils g n raux Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html
7. es donn que la derni re list e dans le EXCIT de CALC ERREUR 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 la divergence du tenseur des contrainte div o est quasi nulle Pour prendre en compte l erreur relative une condition aux limites nulle il faut l imposer en tant que fonction via un AFFE CHAR MECA F Via une constante elle ne sera pas prise en compte Maillage Le maillage doit tre triangulaire quadrangle t tra drique ou hexa drique avec aucun GROUP _ NO si on veut remailler ensuite via HOMARD En 2D il ne prend en compte que les erreurs sur et entre les l ments isoparam triques SEG2 3 TRIA3 6 QUAD4 8 9 En 3D idem avec FACE3 4 6 8 9 TE TRA4 10 PENTA6 13 15 HEXA8 20 27 et PYRAM5 13 mais pas les l ments de structure coque plaque poutre Manuel d utilisation Fascicule u4 81 Outils g n raux Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http
8. es ttes 9 25 Opcrande TITRE tee nee nn ons tnt aiaa Duo unaeente fenetre tte Aa 9 Manuel d utilisation Fascicule u4 81 Outils g n raux Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster Version default Titre Op rateur CALC_ERREUR Date 21 02 2013 Page 3 9 Responsable Josselin DELMAS Cl U4 81 06 R vision 10471 2 Syntaxe resu CALC ERREUR 0 reuse resu MODELE mo modele 0 CHAM MATER chmater cham mater 0 CARA ELEM carac cara_elem 0 SOLVEUR F voir le document U4 50 01 EXCIT F CHARGE l charge 1 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 PUIS PULS n 1 TYPE CHARGE FIXE O S lection des mailles concern es par le calcul TOUT VOUI DEFAUT S lection des num ro d ordre TOUT_ORDRE TOUT y NUME ORDRE l_ nuor 1 I LIST ORDRE l nuor listis NUME MODE 1 numo 1 I NOEUD CMP 1 nomo 1 K16 NOM CAS nocas K16 1 E INST l
9. in aire 2D pour la version 2 partir de l option SIZ1 NOEU respectivement S1Z2 NOEU Si ce dernier champ n existe pas dans resu il est automatiquement construit au pr alable voir R4 10 01 ERME ELEM Estimateur d erreur en r sidu en m canique R4 10 02 et en hydro m canique stationnaire R4 10 04 calcul par l ment Conseils d utilisation de l option ERME ELEM Pour bien effectuer l estimation d erreur du calcul m canique dans les limites th oriques de la formule mise au point dans le cadre elliptique avec fronti re r guli re il faut l effectuer sur tout le mod le TOUT OUI Valeur par d faut A noter que le mod le n est pas forc ment d fini sur toute la g om trie On a galement besoin du calcul des contraintes aux n uds confer R3 06 03 par SIGM ELNO Si ce champ de contraintes aux n uds n existe pas d j dans la structure de donn es resultat il est automatiquement calcul dans CALC ERREUR En ce qui concerne les chargements Il faut fournir CALC ERREUR les chargements utilis s pour le calcul m canique EXCIT F CHARGE en prenant bien garde aux r gles de surcharges diff rentes pour le solveur m canique et pour cette option de CALC ERREUR 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 limit
10. inst 1 R FREQ l freq 1 R LIST_INST l inst listr8 LIST_FREQ 1l freq listr8 P RECISION prec 1 0E 3 DEFAUT CRITERE RELATIF DEFAUT ABSOLU options m caniques RESULTAT resu OPTION SIZ1 NOEU ERZ1 ELEM SIZ2 NOEU ERZ2 ELEM ERME_ELE 1 ERME_ELNO ERME_NOEU QIZ1 ELEM QIZ2 ELEM QIRE_ELE 1 QIRE_ELNO QIRE_NOEU SING ELEM PREC ERR err R TYPE ESTI ERME ELEM ERZ1 ELEM ERZ2 ELEM OIRE ELEM QIZ1 ELEM QIZ2 ELEM Manuel d utilisation Fascicule u4 81 Outils g n raux Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Version default Date 21 02 2013 Page 4 9 Code Aster Cl U4 81 06 R vision 10471 Titre Op rateur CALC ERREUR Responsable Josselin DELMAS options thermiques RESULTAT resu evol ther OPTION ERTH_ELEM ERTH_E NO ERTH_NOEU TITRE titre 7 1_Kn INFO 1 DEFAUT f 2 Fascicule u4 81 Outils g n raux Manuel d utilisation Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster ou Titre Op rateur CALC_ ERREUR Date 21 02 2013 Page 5 9 Responsable Josselin DELMAS Cl U4 81 06 R vision 10471 2 1 Op randes RESULTAT MO
11. lle 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 dans le source d 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 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
12. nitial 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 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 PRE EPSI 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 ERREUR avec une des options de calcul d estimateur d erreur en quantit d int r t Un exemple d utilisation du calcul de l estimateur en quantit s d int r t bas sur les r sidus peut tre trouv dans le test ssiv113c et sslv113d QIRE ELNO Estimateur d erreur en quantit s d int r t bas sur les r sidus calcul s aux n uds par l ment QIRE NOEU Estimateur d erreur en quantit s d int r t bas
13. nts 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 solveurs autoris s sont GCPC LDLT MULT FRONT d faut MUMPS et PETSC 2 2 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 001 Manuel d utilisation Fascicule u4 81 Outils g n raux Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster or Titre Op rateur CALC_ERREUR Date 21 02 2013 Page 6 9 Responsable Josselin DELMAS Cl U4 81 06 R vision 10471 2 3 Op randes pour les options m caniques 2 3 1 Options de calcul d indicateurs d erreur ERZ1 ELEM respectivement ERZ2 ELEM Calcul de l estimateur d erreur de ZHU ZIENKIEWICZ lasticit lin aire 2D et 3D pour la version 1 lasticit l
14. r Cette option est un champ constant par l ment et comporte trois composantes 1 DEGRE qui correspond la d tection des l ments finis singuliers En pratique cette composante vaut le degr d interpolation des l ments finis choisis si l l ment fini n est connect aucune singularit et vaut l ordre de la singularit si l l ment fini est connect un n ud consid r par la m thode comme singulier par exemple pour un l ment voisin de la pointe d une fissure cette valeur vaut 0 5 2 RAPPORT qui correspond la carte de modification de taille des l ments finis en cas de remaillage pour une erreur globale donn e Cette composante est gale au rapport entre la nouvelle taille de l l ment fini et la taille actuelle 3 TAILLE qui 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 mailleurs 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 Code Aster pour que celui ci construise le nouveau maillage en respectant au mieux cette nouve
15. www gnu org copyleft fal html Code Aster ut Titre Op rateur CALC_ERREUR Date 21 02 2013 Page 7 9 Responsable Josselin DELMAS Cl U4 81 06 R vision 10471 e D autre part il faut veiller ne pas intercaler de segments entre deux quadrangles ou deux triangles respectivement QUAD ou triangle entre deux HEXA sinon on ne peut pas calculer le terme de saut relatif ce voisinage la place on s enquiert tort d une ventuelle condition aux limites ERME ELNO Estimateur d erreur en r sidu calcul aux n uds par l ment R4 10 02 ERME NOEU Estimateur d erreur en r sidu calcul aux n uds T QIZ1_ELEM respectivement QIZ2 LEM Estimateur d erreur en quantit s d int r t bas sur la m thode de Zhu Zienkiewicz lasticit lin aire 2D QIRE ELEM Estimateur d erreur en quantit s d int r t bas sur les r sidus en m canique calcul par l ment Conseil d utilisation des options QIZ1_ ELEM OIZ2 ELEM QIRE ELEM Le domaine d utilisation des options 01Z1 ELEM et 01Z2 ELEM est le m me que pour les options ERZ1 ELEM et ERZ2 ELEM et celui de l option OIRE ELEM est le m me que celui de l option ERME ELEM en m canique Il est n cessaire de d finir en plus du probl me i
Download Pdf Manuals
Related Search