4. Modélisation B#

4.1. 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.

4.2. 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 \(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.

4.3. Grandeurs testées et résultats#

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

Maillage

Localisation

Identification

Référence (Analytique)

Tolérance

1

Nœud \(B\)

Contrainte \(\mathrm{SIYY}\)

303.0

5.0%

1

Nœud \(A\)

Contrainte \(\mathit{SIXX}\)

-100.0

15.0%

2

Nœud \(B\)

Contrainte \(\mathrm{SIYY}\)

303.0

2.0%

2

Nœud \(A\)

Contrainte \(\mathit{SIXX}\)

-100.0

3.0%

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