Modélisation B
==============


Déroulement du TP
-----------------

Il s’agit de mettre en œuvre l'adaptation de maillage à la suite d'un calcul élastique.

On repartira donc de l'étude réalisée lors de la modélisation A.


Adaptation de maillage
-----------------------

Après avoir lancé le module AsterStudy, on modifie le fichier de commande pour ajouter le calcul de l'indicateur d'erreur:


*   Ajouter les résultats d’erreur avec la commande CALC_ERREUR dans la catégorie Post Processing: choisir l’option ERME_ELEM en utilisant reuse.


*   Relancer le calcul.

Pour visualiser les champs de contrainte et d’erreur (composante ERREST, erreur totale), on peut effectuer le post-traitement dans l’onglet Results (fortement recommandé) ou dans le module ParaVis. Les étapes dans le module ParaViS sont comme suivant:


*   Cocher les champs à visualiser pour les charger avant d’appliquer GenerateVectors: SIEF_ELGA, ERME_ELEM


*   Pour visualiser la composanteSIYYdu champ SIEF_ELGA notamment au bord du trou: Menu Filters → Mechanics → ELGA field To Surface et passer en représentation Surface With Edges. Pour relever la position de la valeur maximale (de l'ordre de :math:`270\mathit{MPa}` ): Menu Edit → Find Data et choisir l’option Cells dans le champ crée ELGAfieldToSurface, selectionner SIEF_ELGA_Vector(1) représentant SIYY et la formule is max.


*   Visualiser la cartographie de l'indicateur d'erreur (champ ERME_ELEM, composante ERREST, erreur totale).

On retourne au module AsterStudy, on modifie le fichier de commande pour effectuer une adaptation de maillage et calculer ensuite directement la nouvelle valeur d'indicateur d'erreur:


*   Catégorie Post Processing / commande MACR_ADAP_MAIL pour une adaptation de maillage de type RAFFINEMENT: on nomme le nouveau maillage par MAILLAGE_NP1=MAIL2, sélectionner le critère des éléments (option CRIT_RAFF_PE) à 10% pour le paramètre d’adaptation: RESULTAT_IN / NOM_CHAM=ERME_ELEM / sa composante NOM_CMP = ERREST.


*   Enchaîner ensuite les commandes AFFE_MODELE, AFFE_MATERIAU, AFFE_CHAR_MECA, MECA_STATIQUE,CALC_CHAMP et CALC_ERREUR sur ce nouveau maillage MAIL2 (Duplicate les commandes en modifiant le maillage et d’autres termes concernants). **Attendtion** : les noms des concepts de ces commandes ne doivent pas être identiques à ceux d’avant, et ils ne dépassent pas 8 caractères.


*   Imprimer les nouveaux résultats de calcul au format MED: Commande IMPR_RESU.


*   Lancer le calcul.

Pour visualiser les champs de contrainte et d’erreur, on peut effectuer le post-traitement dans l’onglet Results(fortement recommandé) ou dans le module ParaVis.


Grandeurs testées et résultats
------------------------------

Valeur des composantes de contraintes (SIGM_NOEU)avec le maillage initial et après une adaptation:


.. csv-table::

    "**Maillage**", "**Localisation**", "**Identification**", "**Référence (Analytique)**", "**Tolérance**"
    "1", "Nœud :math:`B`", "Contrainte :math:`\mathrm{SIYY}`", "303.0", "5.0%"
    "1", "Nœud :math:`A`", "Contrainte :math:`\mathit{SIXX}`", "-100.0", "15.0%"
    "2", "Nœud :math:`B`", "Contrainte :math:`\mathrm{SIYY}`", "303.0", "2.0%"
    "2", "Nœud :math:`A`", "Contrainte :math:`\mathit{SIXX}`", "-100.0", "3.0%"


Ces tests sur des valeurs analytiques sont doublés par des tests de non-régression.