4. Modélisation B#
4.1. Caractéristiques de la modélisation#
La méthode upwind est utilisée par PROPA_FISS pour résoudre les équations de propagation de la fissure. Une grille auxiliaire est utilisée.
Le même modèle que celui décrit pour la modélisation A est utilisé. Le domaine de calcul est localisé autour du fond de la fissure.
4.2. Caractéristiques du maillage#
On utilise le même maillage que celui de la modélisation A.
La grille auxiliaire utilisée est constituée de 1296 éléments HEXA8 réguliers de dimension \(0.25\times 0.25\times 0.25m\) (voir ).
Figure 4.2-a : maillage utilisé pour définir la grille auxiliaire
La grille est étendue à la seule zone de la structure intéressée par la propagation de la fissure.
4.3. Grandeurs testées et résultats#
On calcule les points d’intersection entre l’arête qui donne la position théorique du fond et les faces des éléments du maillage en utilisant l’opérateur MACR_LIGN_COUPE (voir ). Pour chacun de ces points, on calcule la valeur de la level set normale (\(\mathrm{LSN}\) ) et tangente (\(\mathrm{LST}\) ) en utilisant l’opérateur POST_RELEVE_T et on vérifie que les valeurs maximale et minimale sont presque nulles:
Propag. \(i\) |
Max \({\mathit{LSN}}_{i}\) |
Min \({\mathit{LSN}}_{i}\) |
Max \({\mathit{LST}}_{i}\) |
Min \({\mathit{LST}}_{i}\) |
1 |
9.02E-17 |
-2.77E-16 |
-3.88E-16 |
-5.7E-16 |
2 |
-1.94E-03 |
-2.08E-03 |
-2.09E-06 |
-2.11E-06 |
3 |
9.19E-04 |
-4.52E-05 |
-1.23E-03 |
-1.17E-03 |
Les valeurs obtenues sont calculées à partir des valeurs aux nœuds du maillage en utilisant les fonctions de forme des éléments. On s’attend donc à ce que ces valeurs soient affectées par une erreur qui dépend de la taille des éléments du maillage. En effet la précision de représentation des level set est elle même liée à la taille des éléments. En conséquence on utilise une tolérance pour vérifier si les level sets calculées sont presque nulles. En considérant que le maillage est grossier, on affecte une tolérance égale à 15% de la longueur de la plus grande arête du maillage dans la zone de propagation:
Tolérance utilisée = \(0.15\times 0.33=0.05m\)
4.4. Remarques#
Toutes les valeurs testées respectent la tolérance utilisée. Cela signifie que la position du fond de fissure calculée par la méthode upwind+grille auxiliaire est correcte et que la localisation du domaine marche bien pour la méthode upwind+grille auxiliaire.
Il faut enfin remarquer qu’après les trois propagations simulées, la fissure a dévié de \(90\) . Par contre l’avance totale est petite. Donc une propagation en mode mixte très sévère a été simulée. Les conditions utilisées sont plus sévères que les conditions qu’on trouve normalement pour des structures réelles. Cependant la méthode upwind+grille auxiliaire a bien calculé la position de la fissure et on a donc vérifié sa robustesse.