Importation dans Comsol des données hydrodynamiques Matlab

La résolution des équations de poroélasticité sur la digue est réalisée à l'aide du solveur par éléments finis Comsol 3.5. La résolution de ces équations nécessite la connaissance des champs de pressions et de vitesses de la phase fluide en fonction du temps, qui ont déjà été calculés précédemment en utilisant Matlab.

Le problème ici est donc d'importer les résultats de nos simulations de Dupuit-Boussinesq de Matlab vers Comsol.

Pour cela, nous avons adopté le protocole suivant:

 

   1) Écriture des résultats du programme Matlab dans un fichier texte sous la forme:

Pas de temps Pas d’espace Charge hydraulique
t(i) x(1) H(i,1)
t(i) x(2) H(i,2)
t(i) x(Nx) H(i,Nx)
t(i+1) x(1) H(i+1,1)
t(i+1) x(2) H(i+1,2)

   2) Lecteur du fichier texte en entrée de Comsol:

Onglet Options/ Functions.../ New/  Interpolation - Use data from file

 

   3) Les maillages utilisés sous Matlab et Comsol sont différents: on utilise sous Matlab un maillage 1D  uniforme alors que Comsol utilise un maillage 2D triangulaire (éléments finis) dont la forme et la taille varient dans l'espace. Ainsi, il est nécessaire d'interpoler en espace, et pour chaque pas de temps, les hauteurs d'eau calculées avec Matlab pour les faire correspondre avec le maillage de Comsol.

Interpolation spatiale des hauteurs d'eau Matlab/Comsol

 


Précédent: Mise en équations

Suivant: Création du modèle de digue

Accueil - Etude poroélastique Comsol