La ALE est définie directement dans le jeu de données en python (et donc user-dependant). Le cas test ALE_frictionless permet de visualiser son principe de fonctionnement.

La ALE en X-FEM utilise une partition de l'opérateur éulérien-lagrangien. La phase de relocalisation des noeuds consiste à repositionner le maillage dans sa position initiale. Vu que l'on utilise des XFEM, le domaine maillé peut être un “carré” maillé de manière régulière. La level-set informe de la position courante de la frontière de la structure.

La phase de transfert des données est par contre plus complexe : une fois que l'on a relocalisé les noeuds et reconstruit la nouvelle level set sur le maillage (traduisant la position de la frontière dans la configuration courante), il faut redéfinir les cellules d'intégration sur les éléments découpés par la frontière. Le transfert des données aux points de Gauss est réalisé en utilisant les éléments non-découpés et les cellules d'intégration. Cette phase de transfert s'appuie sur la méthode volumes finis (voir les méthodes de transfert entre maillages)

La phase ALE n'est activée que lorsque le maillage devient trop mauvais.

Pour faire bref, le fonctionnement est le suivant :

• On fait un pas de calcul
• On vérifie la qualité du maillage
• Si le maillage a une qualité insuffisante, on enregistre les frontières actuelles de la structure (frontières matière et frontières de chargement)
• On reconstruit un nouvel objet metafor et un nouveau domain en utilisant ces frontières
• On transfert les données entre les deux maillages
• On reprend le calcul sur ce nouveau maillage

Etant donné qu'on utilise les opérateurs de transfert utilisés dans Metafor, il est nécessaire d'importer mtXFEM et les toolbax de dataTransfer en en-tête du fichier.

 from toolbox.importMetafor_XFEM import *
 from wrap.mtDataTransfer import *
 import toolbox.remeshedUtilities as remeshedUtilities

Comparativement à un cas test X-FEM classique, il faut définir ici les paramètres de l'opérateur de transfert (principalement le nombre de points d'intégration par élément de transfert, étant donné que seule une méthode de transfert FV est pour le moment disponible).

Par ailleurs, deux tolérances sont introduites, correspondant à la valeur maximum du squish index et de l'aspect ratio, pour déterminer si une opération de remaillage est nécessaire

Certaines variables, utilisées par une grande partie des fonctions, sont définies de manière globale :

 oldMetafor = None           # Metafor avant remaillage
 newMetafor = None           # Metafor après remaillage
 transferOperator = None     # Opérateur de transfert
 win = None                  # Fenêtre de visualisation
 no = 0                      # Paramètre de comptage des itérations