Modelisation et identification de la loi de comportement des elastomeres magneto-rheologiques

In this thesis, we study a class of “active materials” called Magnetorheological elastomers (MRE) which are ferromagnetic impregnated rubbers whose mechanical properties are altered by the application of external magnetic fields. With the purpose of characterizing the behavior of MREs up to large strains and high magnetic fields, this work brings a completely novel experimental, theoretical and numerical approach.The first part of this study focuses on an experimental investigation of MRE where multiple microstructures (isotropic and transversely isotropic materials) and multiple particles’ volume fraction are tested. A special sample geometry is designed in order to increase the uniformity of internal magnetic and mechanical fields measured during coupled-field experiments. The interfacial adhesion between the iron fillers and the silicone matrix is investigated and we show that when specimens are subjected to external magnetic fields, a silane primer treatment of the particles is needed to prevent debonding at the interface particle/matrix. Then, we present the magneto-mechanical testing setup that allows simultaneous 3D mechanical and magnetic measurements before discussing the results. Even if is found that instabilities are ubiquitous in MREs, lots of useful data are collected and will be used to compute the parameters proposed in the material model.The second part of the thesis is dedicated to the modeling of isotropic MREs. The continuum description proposed by Kankanala, Triantafyllidis and Danas (2004, 2012, 2014) to derive constitutive laws that account for finite strains is used and, in particular, the energetic approach (that requires an energy density function) is chosen. Multiple equivalent variational formulation alternatives (based on different choices of the independent magnetic variable used in the energy function: B, H or M) are given and implemented into 3D finite element (FEM) codes. Based on the use of FEM simulation in combination with least square optimization methods, the previously collected experimental data are fitted and all three energy functions ψB , ψH and ψM are computed. The obtained material model proves to have excellent predictive capabilities when compared to other experiments not used in the fitting process. The use of numerical tools is necessary to make sure that the calculated material parameters are not influenced by the shape of experimental specimens.The last part of this work details the numerical implementation of the different variational formulations. For each one of them, it is found that isoparametric elements are well suited to simulate coupled magneto-mechanical boundary value problems. We show that special care is needed when implementing variational formulations using the displacement vector and the magnetic vector potential as independent variables. Indeed, ensuring the uniqueness of the vector potential requires to numerically enforce the Coulomb gauge, which leads to numerical complications that are addressed in this thesis. Before describing the different patch tests that have been considered to validate the numerical codes, we show which are the valid boundary conditions for the magnetic vector potential and how to use the symmetry properties of a given boundary value problem to reduce its complexity and the computational resources needed to solve it.Ce travail de thèse porte sur une catégorie de matériaux actifs dénommés Elastomères Magnéto-Rhéologiques (EMR). Ces derniers sont composés de particules micrométriques et magnétisables imprégnées dans une matrice élastomère isolante. Il est possible de modifier les propriétés mécaniques de tels matériaux en les soumettant à un champ magnétique externe. Avec pour objectif d’aboutir à une caractérisation couplée (magnéto-mécanique) du comportement des EMRs en grandes déformations et en présence de champs magnétiques élevés, ce travail propose une approche à la fois expérimentale, théorique et numérique.La première partie de ce travail s’intéresse à des aspects expérimentaux où l’influence de la microstructure (isotrope et transverse isotrope) et l’influence de la fraction volumique de particules sont étudiées. Un échantillon dédié est développé afin d’obtenir simultanément des champs mécaniques et magnétiques les plus homogènes possibles dans celui-ci lors d’une caractérisation couplée. La question de l’adhésion interfaciale entre les particules de fer doux et la matrice en silicone est également traitée et il est montré qu’un traitement chimique des particules est nécessaire afin d’éviter toute décohésion avec la matrice lorsque le matériau est soumis à un champ magnétique externe. Avant d’analyser les données obtenues, le dispositif expérimental permettant d’obtenir de manière simultanée une mesure du champ de déformation en trois dimensions et une mesure des champs magnétiques internes, est décrit. Malgré l’ensemble des difficultés expérimentales en grande partie dûes à des phénomènes d’instabilité qui sont omniprésents chez les EMRs, de nombreuses données sont collectées et serviront à la calibration des lois de comportement.La seconde partie de cette thèse couvre la modélisation couplée magnéto-mécanique des EMRs en s’appuyant sur le cadre théorique général des solides magnéto-élastiques proposé par Kankanala, Triantafyllidis et Danas (2004, 2012, 2014). En particulier, la méthode énergétique (qui s’appuie sur l’utilisation d’une fonction d’énergie libre) est préférée et des formulations variationnelles équivalentes (qui diffèrent entre elles simplement par le choix de la variable magnétique indépendante utilisée pour décrire le problème : B, H ou M) sont proposées et implémentées dans des codes numériques 3D s’appuyant sur la méthode des éléments finis. Ces outils numériques sont combinés à la méthode de minimisation des moindres carrés afin d’obtenir l’ensemble des paramètres matériaux du modèle de comportement des EMRs. L’utilisation de simulations numériques est nécessaire car une approche purement analytique ne permettrait pas de modéliser « l’effet de forme » observé expérimentalement. En effet, il est primordial de modéliser ce dernier car dans le cas contraire les paramètres identifiés dépendraient de la forme de l’échantillon expérimental et ne décriraient pas uniquement le matériau.La troisième partie de cette étude décrit en détail l’implémentation numérique des différentes formulations variationnelles proposées précédemment. Dans chacun des cas, il est prouvé que l’utilisation d’éléments isoparamétriques est bien adaptée. De nombreuses difficultés numériques ont été observées dans le cas des formulations variationnelles utilisant le champ de déplacement et le potentiel vecteur magnétique comme variables indépendantes. L’ensemble de ces difficultés (comme par exemple la minimisation de l’énergie potentielle sous la contrainte imposée par la jauge de Coulomb) est surmonté dans ce travail. Avant de décrire les différents problèmes tests utilisés pour s’assurer de la validité et de la précision des codes numériques, les différentes étapes nécessaires à la simulation d’un problème aux limites sont expliquées. Plus précisément, les questions liées aux spécificités des conditions aux limites à appliquer sur le potentiel vecteur magnétique ou encore aux conditions de symétries, sont traitées