Le renforcement des sols par géosynthétique est appliqué dans de nombreux types d’ouvrage : remblais sur sol compressible, talus sur fondations stables, remblais sur des cavités et ouvrages de soutènement. La stabilité de ces ouvrages dépend entre autres de l’efficacité des ancrages des nappes géosynthétiques. Les ancrages droit et avec retour sont les plus couramment utilisés. Afin d'améliorer les connaissances actuelles sur le comportement des systèmes d'ancrage, des études expérimentales et numériques ont été développées conjointement. Ce travail de thèse concerne dans une première partie, la modélisation physique tridimensionnelle du comportement des géosynthétiques pour deux types ancrages (ancrage droit et ancrage avec retour). Ces essais ont été réalisés dans une chambre d’étalonnage en conditions contrôlées et instrumentées en laboratoire. Dans une deuxième partie de cette thèse, les paramètres d’interaction sol/géosynthétique déduits à partir de l’étude expérimentale ont été implémentés dans le code de calcul numérique bidimensionnel en milieu continu FLAC2D pour une meilleure compréhension du comportement des géosynthétique en ancrage. L’influence de plusieurs paramètres sur le comportement du géosynthétique en ancrage avec et sans retour a été traitée dans cette étude numérique. Parallèlement à cette modélisation, une autre modélisation numérique (discontinue) par la méthode des éléments discrets (PFC2D) a été réalisée. Ces modélisations ont donné des résultats proches de ceux obtenus expérimentalement et confirme l'analyse faite au sujet des mécanismes d'ancrage. / The soil reinforcement by geosynthetic is used in many types of structures: embankments on compressible soil, slope on a stable foundation, embankments on cavities and retaining structures. The stability of these structures specially depends on the efficiency of the anchors holding the geosynthetic sheets. The simple run-out and anchorage with wrap around are most commonly used. In order to improve the available knowledge of the anchorage systems behaviour, experimental and numerical studies were developed jointly. This thesis work concerns in the first part a three-dimensional physical modelling of the behaviour of geosynthetics in two anchors (simple run-out and anchorge with wrap around). The pull-out tests were performed in a test tank under controlled conditions in the laboratory. In the second part, the parameters of the interaction soil/geosynthetic found from the experimental study were used into the numerical code “FLAC2D” (continuous) for a better understanding of the behaviour of geosynthetics in anchorage. The influence of several parameters on the behaviour of geosynthetic was treated in this numerical study. In parallel with this model, another numerical modelling (discontinuous) by the discrete element method (PFC2D) was carried out. The results of these models are closely to experimental results which confirm the analysis about the anchoring mechanisms.

Ce travail est une contribution à l'étude numérique de la propagation des ondes sismiques dans les formations géologiques superficielles. Dans la première partie, un programme de calcul combinant les méthodes d'Aki-Larner et de Thomson-Haskell a été développé. Il permet de calculer la réponse de configurations géologiques géométriquement complexes à des ondes SH planes. Une étude systématique de quelques configurations géologiques-type montre que les effets dûs aux hétérogénéités verticales (contraste d'impédance: amplification) et latérales (irrégularités géométriques: diffraction) sont fortement couplés : il est le plus souvent impossible de décomposer le calcul d'une structure complexe en 2 sous-problèmes, plus simples. En génie parasismique, il est donc nécessaire de traiter chaque site cas par cas pour les études détaillées. Dans la deuxième partie, on étudie l'influence de la saturation en eau dans les formations poreuses très perméables. La représentation du milieu saturé est celle de Biot, avec quelques compléments issus de la théorie de l'homogénéisation. Ce modèle est appliqué au calcul de la réponse sismique de structures simples tricouches (sol poreux sec sur sol poreux saturé sur bedrock élastique imperméable) en 1 D ou 2D, dans la gamme 0 - 25 HZ. On montre des cas théoriques pour lesquels l'onde P2 (générée aux interfaces du milieu poreux, puis très atténuée et non observable en surface) intervient dans le mécanisme d'atténuation des ondes de compression. On insiste sur l'importance de l'interaction fluide-solide aux frontières du milieu poreux: dès que la perméabilité est supérieure à 10-10 m2 environ, il s'avère nécessaire de prendre en compte explicitement la présence de fluide. En particulier, il est alors impossible (sous peine d'erreurs supérieures à 100 %) d'approximer le matériau poreux par un modèle monophasique viscoélastique équivalent.

Nos travaux se sont intéressés à la modélisation de l'apparition des contraintes résiduelles dans le procédé d'empilement par soudage et consolidation en continu développé dans l'industrie aéronautique. Dans les conditions standard d'élaboration, les pièces réalisées par ce type de procédé sont le siège d'importantes contraintes résiduelles. L'ignorance de leur origine et développement est un frein important à la validation industriel de ce procédé prometteur. Dans cette perspective, l'originalité de l'étude a été de développer une méthode de modélisation numérique multi-échelle et multi-physique permettant de réaliser une modélisation couplée aux échelles macroscopique et microscopique du phénomène de l'apparition des contraintes résiduelles, au cours du procédé. L'échelle microscopique, décrite à l'aide de la méthode des éléments naturels contraints (CNEM), apporte à l'échelle macroscopique les propriétés thermomécaniques homogénéisées du matériau à chaque pas de temps. L'échelle macroscopique apporte à l'échelle microscopique les conditions aux limites (températures, déplacements), qui permettent de déterminer les champs de température, de déformation et de contrainte microscopiques dans le matériau au cours du temps. Les résultats obtenus en terme de validation et d'application de la méthode au procédé d'empilement par chauffage et consolidation en continu sont satisfaisants et prometteurs. La méthode développée peut de plus être aisément appliquée à d'autres types de procédé de mise en forme des composites thermoplastiques.

L'utilisation de stents actifs (DES) a révolutionné le traitement de l'athérosclérose. Le relargage contrôlé de médicaments anti-prolifératifs dans la paroi artérielle (PA) a permis de réduire fortement le taux de resténose intra-stent. Mais le risque de thromboses intra-stents tardives demeure un enjeu majeur des DES en partie lié au retard de cicatrisation de la PA endommagée lors de l'implantation. Cette thèse présente une méthode d'optimisation du design des DES afin d'inhiber la resténose sans affecter la cicatrisation. Pour quantifier la performance des différents designs, un modèle numérique décrivant l'écoulement sanguin et le transport de médicaments dans les artères stentées a été développé. Il prend en compte la structure multi-couches de la PA et les interactions du médicament avec les cellules. Un algorithme d'optimisation est couplé au modèle afin d'identifier les DES optimaux. L'optimisation du temps de relargage ainsi que de la concentration initiale du médicament dans le revêtement du DES ont un effet significatif sur la performance. Lorsque le médicament utilisé est le paclitaxel, les solutions optimales consistent à relarguer le produit à des concentrations nettement inférieures à celles des DES actuels soit pendant quelques heures, soit pendant une durée d'un an. Pour le sirolimus, un relargage lent est nécessaire. Les formes optimales des spires du DES sont toujours allongées mais profilées seulement lorsque le relargage est rapide. Ces résultats permettent d'expliquer en partie les performances des différents DES récents et révèlent un fort potentiel d'amélioration dans la conception des DES par rapport aux dispositifs commerciaux actuels.

Le travail présenté dans ce manuscrit s'articule en trois axes distincts. 1) Dérivation de modèles mathématiques de phénomènes électrophysiologiques en cardiologie. Nous utilisons des méthodes d'analyse asymptotique pour dériver un modèle simplifié à partir d'un modèle de tissu auriculaire tridimensionnel, tout en contrôlant l'erreur d'approximation. Ces méthodes ont permis de dériver un modèle bisurfacique qui permet de simuler des comportements tridimensionnels dans les oreillettes pour un coût numérique bidimensionnel afin d'étudier des phénomènes entrant en jeu lors d'arythmies auriculaires, tels que la dissociation électrique ou des hétérogénéités transmurales. La preuve de la convergence du modèle bisurfacique est apportée, et une stratégie d'optimisation du modèle en dehors du régime symptotique est formalisée. Une méthode d'homogénéisation est également utilisée pour construire un modèle continu homogénéisé de l'activité des myocytes incluant le comportement non linéaire des gap junctions. 2)Processus déclencheurs d'arythmie. Des preuves de concepts de mécanismes arythmogènes sont apportées à l'aide de modèles numériques des veines pulmonaires. Le premier mécanisme repose sur un bloc de conduction unidirectionnel engendré par une discontinuité dans la structure fibreuse. Le second est basé sur une dynamique différente lors de la dépolarisation et de la repolarisation lorsque deux couches de fibres de directions différentes sont superposées. 3)Perpétuation des arythmies auriculaires. A partir d'un modèle bicouche des oreillettes, nous étudions l'influence d'hétérogénéités transmurales de fibrose sur la perpétuation des arythmies. Plusieurs protocoles d'ablation sont ensuite testés. Enfin, une méthode de personnalisation du modèle auriculaire est formalisée.

Les volcans boucliers basaltiques se déforment sous leur propre poids, tout particulièrement en période d'éruption. A Hawaii, cette déformation co-éruptive est en majeure partie expliquée par un plan de décollement à la base des édifices, dont le glissement serait associé à la force de poussée des intrusions dans les " rift zones " (i.e. concentration de dykes subverticaux dans les zones de faiblesse de l'édifice). Cependant, cette explication n'est pas valable pour beaucoup d'autres volcans basaltiques, où l'existence d'un décollement est peu probable. Nous avons profité de l'incision intense du Piton des Neiges (La Réunion) par l'érosion pour observer la structure interne d'un volcan basaltique, et comprendre, par le terrain et la modélisation numérique, comment la déformation d'un tel édifice s'articule avec son activité magmatique. Notre étude structurale révèle que trois populations d'intrusions perpendiculaires (deux subverticales et une subhorizontale) coexistent, et que les principales zones de faiblesses sont des " sill zones " (i.e. concentrations d'intrusions subhorizontales), Parallèlement, notre étude des déformations, essentiellement cassantes, révèle que deux extensions perpendiculaires dominent dans l'édifice, avec l'apparition de régimes localement décrochants et même compressifs à proximité des sill zones. Les directions de déformations et les directions d'intrusions sont compatibles et coïncident aussi avec la direction d'écoulement des grandes avalanches de débris qui ont affectée le Piton des Neiges au cours de son histoire. A partir de ces données de terrain, nous proposons un modèle de déformation du Piton des Neiges où les intrusions dans l'édifice génèrent des cycles de permutations de contraintes. Le stade ultime de chaque cycle serait la mise en place de sills en régime compressif. En testant numériquement l'effet de telles injections sur la déformation d'un édifice, nous montrons que les sills peuvent être un facteur majeur d'instabilité, capable de conduire à des déstabilisations de flanc telles que les avalanches observées. Notre modèle conceptuel, déduit du terrain et quantifié par la modélisation, constitue ainsi une alternative au modèle de déformation hawaiien, applicable sur des édifices tels que le Piton de la Fournaise (La Réunion), Tenerife (Canaries) et Fernandina (Galapagos). Notre étude démontre enfin l'intérêt essentiel d'étudier les systèmes anciens pour mieux comprendre les volcans actifs.

Les mécanismes d'érosion sont la principale cause de rupture des ouvrages hydrauliques en remblai, c'est pourquoi il est capital de pouvoir quantifier la résistance des sols à l'érosion. Divers appareillages permettent d'en obtenir une estimation, notamment le Jet Erosion Test (JET) dont le modèle d'interprétation est empirique. L'objectif de ce travail est de statuer sur la pertinence du modèle d'interprétation du JET. Pour cela, nous avons développé un modèle numérique 2D de type Navier-Stokes turbulent avec déplacement d'interface et remaillage, permettant de modéliser l'érosion d'un sol cohésif par un écoulement turbulent. En injectant dans notre modèle les paramètres d'érosion trouvés à la suite d'essais de JET, on retrouve numériquement l'évolution de la profondeur d'affouillement obtenue expérimentalement. Nous avons appliqué cette méthode à trois différents essais de JET, les résultats obtenus sont en bon accord avec les résultats expérimentaux. En plus d'une importante étude paramétrique, ces résultats ont permis d'apporter d'importants éléments de validation au modèle d'interprétation du JET. Pour étendre le champ d'applicabilité du modèle, nous avons également appliqué cette méthode de modélisation à la configuration de l'érosion de conduit. Trois essais de Hole Erosion Test (HET) ont été modélisés et nous avons également obtenus des résultats en bon accord avec les résultats expérimentaux. Une analyse de la loi d'érosion et des paramètres d'érosion obtenus à la suite d'essais de JET et de HET a ensuite été initiée. Les bases d'une étude portant sur l'influence de l'angle d'incidence de l'écoulement, sur l'efficacité de l'érosion, ont été posées.

Une onde sismique se propageant dans un milieu poreux contenant un uide crée un mouvement relatif entre ce uide et la matrice solide. En transportant les ions présents en excès dans la double couche électrique, ce mou- vement relatif crée un courant d'électro ltration. Ce phénomène électrocinétique, survenant à l'échelle microscopique, est à l'origine d'au moins deux types de conversions sismo-électromagnétiques observables à l'échelle macroscopique. D'une part, il donne naissance à des champs électromagnétiques cosismiques accompagnant l'onde sismique perturba- trice. D'autre part, lorsqu'une onde sismique traverse une interface séparant deux milieux de propriétés mécaniques, électriques ou hydrologiques di érentes, elle génère une discontinuité du potentiel de part et d'autre de cette interface, qui agit comme une source électromagnétique dipolaire. Les réponses interfaciales ainsi créées voyagent à des vitesses électromagnétiques et sont observées presque simultanément par tous les récepteurs. Cette thèse s'attache à la com- préhension de ces phénomènes et au développement des méthodes d'imagerie géophysique basées sur les conversions sismoélectromagnétiques. Deux axes de travail sont investigués. Le premier concerne le traitement des données sismoélectromagnétiques et en particulier l'extraction des réponses interfaciales, qui peuvent renseigner sur le sous-sol en profondeur. Nous comparons l'impact de plusieurs méthodes de ltrage sur les amplitudes et les formes d'onde des réponses interfaciales. Nous introduisons également une nouvelle technique de séparation d'ondes en domaine des curvelets, qui s'appuie sur les similarités entre les données sismiques et cosismiques. D'autre part, la théorie de la propagation couplée des ondes sismiques et électromagnétiques est limitée aux milieux poreux saturés. Nous l'étendons ici à des milieux partiellement saturés, pour le cas d'un mélange air/eau, en calculant les propriétés mécaniques e ectives du mélange et en écrivant la conductivité électrique, la permittivité diélectrique et le couplage sismoélectrique en fonction de la saturation en eau. Nous adaptons aux conditions non-saturées un programme de simulation numérique des ondes sismo-électromagnétiques basé sur la méthode de la ré ectivité générale. Mots clés : Géophysique de proche surface - Signaux sismo-électromagnétiques - Milieux poreux - Ondelettes - Séparation d'ondes - Modélisation numérique - Zone non-saturée

Les travaux de recherche que je mène depuis une dizaine d'années, au Laboratoire d'Ingénierie des Matériaux de Bretagne, portent principalement sur deux domaines : la modélisation des procédés de mise en forme des matériaux métalliques et la simulation du comportement de matériaux polymères. Sachant que les conditions d'emplois et de mise en forme des matériaux sont de plus en plus sévères, l'objectif de ces thématiques est d'améliorer la caractérisation et la modélisation du comportement mécanique des matériaux. Tout d'abord, en cherchant à intégrer des phénomènes mécaniques mal maîtrisés (effets de la vitesse de chargement élevée, sensibilité à la température et dépendance du chemin de chargement) ensuite, en développant des méthodes expérimentales ad hoc permettant de mesurer ces phénomènes, et enfin, en mettant en oeuvre des stratégies numériques originales, en tentant de simuler au mieux ces phénomènes. Ces travaux, s'effectuant dans des domaines de transformations finies, un formalisme en grandes transformations s'appuyant sur la description du mouvement du solide par les coordonnées matérielles entraînées a d'abord été précisé. L'étude de l'intégration temporelle du taux de déformation dans différents repères objectifs a montré l'importance du repère en rotation logarithmique et que la dérivée dans le repère corotationnel n'en était qu'une approximation. Ces restrictions à l'esprit, la modélisation de la prédiction du retour élastique lors de la mise en forme des alliages d'aluminium et de la prise en compte de la température au cours de ces étapes a pu être évaluée à travers le test d'ouverture d'un anneau découpé dans un godet cylindrique. L'influence combinée des modèles de comportement, des paramètres numériques, des modes de découpe numérique et du type de modélisation thermomécanique a été analysée, montrant qu'il reste encore des progrès à réaliser pour la modélisation du retour élastique. Enfin, un modèle de comportement d'Hyperélasto-Visco-Hystérétique, discernant les phénomènes mécaniques observés grâce à plusieurs contributions en contrainte, a été appliqué sur deux matériaux polymères. Une méthode d'identification simple, rapide et nécessitant un seul essai de chargement interrompu par des phases de relaxation a permis de montrer l'existence dans ces matériaux des effets de la viscosité, de l'hyperélasticité et d'une partie hystérétique non-visqueuse. La prise en compte des effets de la température s'est également avérée concluante.

Mes activités de recherche se sont centrées sur les procédés de transformation des polymères, essentiellement le procédé d'injection-moulage mais aussi le soufflage et maintenant la peinture. J'ai participé à plusieurs projets de recherche appliquée pluridisciplinaire donc à finalité industrielle. Dans ces projets ma contribution porte sur la modélisation du procédé, que ce soit par l'utilisation de code de calculs ou le développement de modèles numériques à l'aide de logiciels libres. L'originalité de mon travail réside dans la pluridisciplinarité des projets qui mêlent des aspects théoriques et numériques aux aspects expérimentaux; la science des polymères à la thermique, à la mécanique et à la rhéologie. Les modélisations s'appuient sur des mesures expérimentales qui sont, dans les procédés de transformation, particulièrement délicats à isoler et à interpréter. Ces dispositifs expérimentaux représentent une part considérable de l'effort de recherche dans ces projets. Durant ces travaux, j'ai tenté de répondre à des problématiques distinctes : Dans les procédés de plasturgie, la viscoélasticité est finalement assez absente. Lorsqu'elle se manifeste, c'est sous la forme de défauts (instabilités, recirculations) qu'il faut éviter. C'est pourquoi, la plupart des écoulements impliqués dans les procédés sont de type lubrification~: une dimension au moins est petite devant les autres; ce qui atténue de façon effective les effets viscoélastiques. Une illustration frappante en est donnée par les recirculations viscoélastiques créées par la courbure du chenal d'une vis : lorsque le rapport d'aspect de la section augmente, elles s'atténuent et se décalent vers les parois latérales, précisément là où dans la section, aucune dimension n'est petite devant l'autre. Seul le procédé d'extrusion soufflage s'appuie totalement sur les propriétés viscoélastiques du polymère en écoulement, au travers des phénomènes de gonflement et de fluage. Cet écoulement est particulièrement discriminant et doit permettre à terme de déterminer les paramètres matériels d'une loi viscoélastique intégrale en comparant les prédictions des simulation numériques (par le code FEM Polyflow ou la méthode des tubes de courant) aux mesures expérimentales du gonflement et de fluage, devenues réalisables désormais avec notre technique de visualisation par illumination laser. En ce qui concerne le phénomène de plastification dans les monovis d'injection, ce thème de recherche est celui qui a nécessité le plus d'efforts et c'est celui qui est le plus ardu. La difficulté première, celle des observations, est en passe d'être surmontée avec le fourreau à fenêtres qui permet une visualisation inégalée des phénomènes de plastification. Nous allons pouvoir reprendre et vérifier toutes les hypothèses usuelles de la plastification : accélération du lit solide, fonte contiguë ou dispersée, position et épaisseur des films fluides. En ce qui concerne la modélisation numérique du convoyage solide du polymère dans le chenal de la vis. Une autre problématique, très importante dans le domaine industriel de l'automobile ou de la cosmétique, concerne les défauts de surfaces des pièces plastiques injectées. De nombreux facteurs peuvent contribuer à l'apparition de tels défauts. Nous avons étudié les effets thermiques en montrant que même des revêtements métalliques extrêmement fins pouvaient avoir une influence thermique. Qu'en est-il de l'aspect dynamique de l'écoulement et du contact (tension de surface) polymère-métal ? Ce sujet est d'importance pour l'injection moulage classique mais aussi la micro-plasturgie. La simulation numérique de la projection de peinture représente une direction nouvelle de recherche, active en Europe, mais inédite en France, et pour laquelle nous sommes bien placés à l'insa-lyon avec la création récente d'une cabine de peinture pilote. Cela reste toujours de la mécanique des fluides et des polymères mais avec des outils différents (code volumes finis, OpenFOAM, StarCCM). Encore une fois, une analyse des phénomènes prépondérants, notamment dans l'aérodynamique et le champ électrostatique, associée à des mesures expérimentales, devraient aboutir à un modèle simplifié des trajectoires des gouttes et donc de la projection de peinture.