Formulation faible, formulation forte et méthodes de calcul numérique du champ magnétique en vue de la modélisation d'une prothèse innovante pour l'assistance circulatoire

Gérard Vinsard

Description

Niveau: Supérieur, Doctorat, Bac+8
THÈSE En vue de l'obtention du DOCTORAT DE L'UNIVERSITÉ DE TOULOUSE Délivré par L'Institut National Polytechnique de Toulouse Discipline ou spécialité : Génie Electrique JURY Pr. CAPOLINO Gérard-André (Rapporteur et Président du jury) Dr. VINSARD Gérard (Raapporteur) Pr. LEPRINCE Pascal Pr. NOGAREDE Bertrand Dr. LEFEVRE Yvan Ecole doctorale : GEET Unité de recherche : Laboratoire LAPLACE Directeur(s) de Thèse : Dr. LEFEVRE Yvan Présentée et soutenue par Nicolas MARTINEZ Le 6 Novembre 2008 Titre : Formulation faible, formulation forte et méthodes de calcul numérique du champ magnétique en vue de la modélisation d'une prothèse innovante pour l'assistance circulatoire

• code de calcul


• thoracique de l'ap-hp pitié-salpêtrière de paris


• simulation numérique


• champ magnétique en potentiel scalaire


• répartition du champ magnétique sur le corps déformable

Sujets

• Institut national polytechnique de Toulouse
• Martínez
• Lefèvre
• Université de Toulouse
• Simulation informatique

Informations

• Publié par: profil-feym-2012
• Publié le: 01 novembre 2008
• Nombre de lectures: 114
• Langue: Français
• Poids de l'ouvrage: 3 Mo

Extrait

THÈSE


En vue de l'obtention du

DOCTORAT DE L’UNIVERSITÉ DE TOULOUSE

Délivré par L'Institut National Polytechnique de Toulouse
Discipline ou spécialité : Génie Electrique


Présentée et soutenue par Nicolas MARTINEZ
Le 6 Novembre 2008

Titre : Formulation faible, formulation forte et méthodes de calcul numérique du champ
magnétique en vue de la modélisation d'une prothèse innovante pour l'assistance
circulatoire

JURY
Pr. CAPOLINO Gérard-André (Rapporteur et Président du jury)
Dr. VINSARD Gérard (Raapporteur)
Pr. LEPRINCE Pascal
Pr. NOGAREDE Bertrand
Dr. LEFEVRE Yvan

Ecole doctorale : GEET
Unité de recherche : Laboratoire LAPLACE
Directeur(s) de Thèse : Dr. LEFEVRE Yvan
Résumé
L’objectif de ce travail est d’améliorer les méthodes numériques de calcul de
champ magnétique en vue de l’estimation de la répartition des forces magnétiques
dans les actionneurs, notamment ceux conçus sur de nouveaux concepts à base de
matériaux magnétoactifs. Cette étude se place dans le cadre d’une collaboration
entre le groupe de recherche en Electrodynamique du laboratoire LAPLACE et le
service de chirurgie Cardiovasculaire et Thoracique de l’AP-HP Pitié-Salpêtrière
de Paris.
En effet, au cours de ces dernières décennies de nombreux progrès ont été réa-
lisés dans le domaine de l’assistance circulatoire. Toutefois, l’implantation per-
manente de ces nouveaux systèmes se heurte encore à de nombreuses barrières :
encombrement, difficulté et lourdeur de l’implantation, fiabilité, systèmes non phy-
siologiques n’assurant que des débits continus, efficacité insuffisante. Ainsi, de nou-
veaux concepts de pompes sont aujourd’hui à l’étude, en vue d’assurer un débit
pulsé physiologique à l’aide d’un système totalement implantable. L’approche pri-
vilégiée dans notre étude est basée sur le "morphing" électroactif d’un corps défor-
mableremplidefluidemagnétorhéologique,dontledéplacementetlaconformation
spatiale sont contrôlés, par l’intermédiaire du champ magnétique, par des bobines
externes. Des simulations par éléments finis ont été réalisées pour montrer la fai-
sabilité des principaux concepts et ont permis le dimensionnement d’un premier
démonstrateur.
Lors de ces simulations numériques réalisées avec ANSYS Multiphysics, nous
avons eu besoin de connaître très précisément la répartition du champ magnétique
surle corpsdéformable pourcalculer la distribution des forcess’y exerçant etainsi,
déterminer sa déformation. Une analyse approfondie des résultats de calcul a mis
en évidence une discontinuité importante de la composante tangentielle du champ
d’excitationmagnétiqueauniveaudel’interfacededeuxmatériauxdeperméabilité
différente.
Bien que constituant une limite bien connue de la formulation faible utilisée
courammentpourlecalculparélémentsfinisduchampmagnétique,ladiscontinuité
de cette composante n’a pas été étudiée de manière systématique jusqu’à présent.
Or, dans la plupart des méthodes de calcul de force magnétique, l’expression de
la densité de force s’exerçant à la surface du corps polarisé fait intervenir cette
composante et la discontinuité constatée peut entraîner des erreurs importantes.
iPourcomprendreceproblème,nousavonsdéveloppéuncodedecalculparéléments
finis permettant la détermination du champ, soit en termes de potentiel scalaire,
soit en termes de potentiel vecteur. Grâce à ce code, une étude a pu être menée
de façon systématique pour analyser et quantifier cette discontinuité ainsi que son
influence sur le calcul des densités de force magnétique selon le maillage et le
potentiel utilisés.
Cettediscontinuité étantdirectement liéeàl’utilisationdelaformulationfaible,
la solution proposée pour palier ce problème est l’élaboration et la mise en oeuvre
d’un code de calcul utilisant une formulation forte. A l’instar du code précédem-
mentdéveloppé enéléments finis,cetoutilpermetdecalculerlechampmagnétique
en potentiel scalaire ou en potentiel vecteur. Une étude comparative des deux for-
mulations,faibleetforte,apuêtremenée,montrantlasupérioritédelaformulation
forte en ce qui concerne le respect des conditions de passage aux interfaces de deux
milieux magnétiques deperméabilité très différente, etce,sans dégraderlasolution
dans le reste du domaine d’étude.
Dans une dernière partie, les résultats de ce code ont été confrontés à
l’expérience. Cette dernière étape conforte les conclusions précédentes quant à
l’amélioration des procédures de calcul numérique du champ. La mise en oeuvre
des solutions préconisées constitue un outil précieux pour engager les études
d’optimisation du nouveau concept d’assistance circulatoire présenté au début de
la thèse.
Mots clés : Formulation faible, Formulation forte, Méthode des éléments finis,
Méthodedesdifférences finies,Fluidemagnétorhéologique,Assistancecirculatoire.
iiAbstract
This thesis focuses on the improvement of numerical methods to compute
the magnetic field in order to estimate the distribution of magnetic forces in
actuators, including those designed on innovative concepts based on the use of
magnetoactive materials. Thus, in the field of mechanical circulatory support,
new implantable cardiac prosthesis are under consideration to ensure a pulsating
flow. Here, the approach is based on the electroactive "morphing" of a deformable
body, so the displacement and the deformation can be modified at distance. The
repartition of magnetic field around the deformable body must be precisely known
to compute the repartition of forces and to determine its deformation. An analysis
of the results obtained from finite element method, highlighted a discontinuity
of the tangential component of magnetic field excitation at the interface of two
materials of different magnetic permeability. To understand this problem due to
the use of a weak formulation, a precise study was led in order to analyze this
discontinuity and its influence on the calculation of the magnetic force density. To
overcome this problem, the proposed solution is to develop a code using a strong
formulation. A comparative study between the two formulations and experimental
value has been completed to show the superiority of the strong formulation. The
solution is a valuable tool to initiate the optimization studies of the new concept
of prosthesis presented at the beginning of the thesis.
Keywords : Weak formulation, strong formulation, Finite Element method,
Finite Difference method, Continuity conditions, Mechanical circulatory support,
Cardiac Prothesis, Magnetorheological Fluid.
iiiA ma femme Amélie,
à mes enfants, David, Abigaëlle et Alissa,
à mes parents,
vRemerciements
Le travailprésenté dans ce mémoire a été réalisé au sein du groupede recherche
en Electrodynamique - GREM3 du LAboratoire PLAsma et Conversion de l’Ener-
gie (LAPLACE) à l’ENSEEIHT de Toulouse. Il a été réalisé en collaboration
avec le service de Chirurgie Cardiovasculaire et Thoracique du centre Hospitalier
Pitié-Salpêtrière APHP, Paris.
Je tiens tout d’abord à remercier le Pr. Bertrand Nogarède, Professeur des
universités à l’INPT (Institut National Polytechnique de Toulouse) pour la
confiance qu’il m’a accordée en me proposant de réaliser cette thèse au sein de
son groupe de recherche en Electrodynamique. De plus, son enthousiasme, son
humanité et son dynamisme m’ont donné les motivations nécessaires pour réaliser
ces travaux.
Je tiens également à remercier M. Yvan Lefèvre, Chargé de Recherche au
CNRS, mon directeur de thèse qui m’a accompagné tout au long de ces 3 ans,
et qui a toujours été présent pour m’aider, me supporter et surtout m’ouvrir les
portes du merveilleux domaine de la modélisation numérique. Merci.
Je souhaite aussi remercier les membres du jury :
– M. Capolino, Professeur à l’université de Picardie - Jule Verne, pour avoir
accepté le rôle de rapporteur et contribué à améliorer la qualité de ce mé-
moire.
– M. Vinsard, Maitre de conférence à l’Institut National Polytechnique de Lor-
raine - INPL, pour avoir également accepté la charge de rapporteur et p