Dans le contexte de la biomécanique, les outils numériques constituent des moyens puissants et indispensables dans la compréhension des mécanismes de blessures. Ils permettent de pallier les freins que sont les expérimentations sur l'humain, liés à des raisons d'éthique qui limitent la possibilité d'essais sur des SHPM (Sujets Humain Post Mortem). Le développement de ces outils numériques a conduit à celui de plusieurs mannequins numériques permettant de stimuler diverses sollicitations (civiles ou militaires), nous donnant ainsi accès à des limites de tolérances.En vue d'explorer la réponse dynamique du corps humain soumis à des sollicitations diverses, un modèle de mannequin numérique a été développé au sein du laboratoire. Ce travail de thèse tente donc d'apporter une contribution dans la recherche sur la définition d'un critère de blessure et l'établissement de limites de tolérance du corps humain soumis aux chargements violents de la partie thoracique dans des contextes militaires. / The development of computer science has allowed an increase in the use of numerical approaches such as finite elements method in order to understand physical mechanisms. These numerical tools are often used to extend and complete experimental investigations wich are limited because of high financial cost and ethical issues. Thus, the use of simulation to avoid thes limitations becomes essential in biomechanics investigations. Many numericalmodels of the thorax/abdomen system have been developped over the last two decades. In that framework, a finite element model of the human body, dedicated to high speed loadings, has been developed in the laboratory. In this context, the objective of this Ph.D Thesis is to investigate the consequences of such loadings on the human body and to contribute to the research of injuries criteria and tolerance limit definition.

La magnéto-tellurique (MT) (Cagniard 1953, Tikhonov 1950) est une technique d'exploration de la Terre basée sur des mesures de champs électromagnétiques (EM). Une source naturelle (non artificielle) harmonique en temps et située dans l'ionosphère (Weaver 1994) produit un champ EM régi par les équations de Maxwell. Les champs électromagnétiques sont enregistrés par plusieurs récepteurs placés sur la surface de la Terre. Ces mesures sont utilisées pour produire une image du sous-sol à partir d'un procédé d'inversion utilisant des méthodes numériques. Nous utilisons la méthode hp-FEM résultant d'une extension du travail de Demkowicz 2005. Nous avons développé un logiciel qui résout, pour la première fois, le problème MT avec des éléments finis auto-adaptatifs. La méthode hp-FEM permet des raffinements locaux, à la fois en taille h et en ordre p sur les éléments, ce qui est un avantage notoire puisque la combinaison de ces deux types de critères permet de mieux capter la présence de singularités, fournissant ainsi des erreurs de discrétisation faible. C'est donc une méthode très précise dont la convergence est exponentielle (Gui and Babuska 1986, Babuska and Guo 1996). En raison des défis d'implémentation encore non résolus (Demkowicz et al. 2002) et de la complexité technique des calculs hp-FEM en 3D, nous nous limitons, dans ce travail, à des calculs en 1D et 2D.Le domaine de calcul est tronqué par un matériau absorbant (Perfectly Matched Layer PML, Berenger 1994), qui est conçu pour s'adapter automatiquement aux propriétés physiques des matériaux. En particulier, il s'ajuste efficacement à l'interface air-sol, où le contraste entre la conductivité des matériaux atteint jusqu'à seize ordres de grandeur. Dans cette thèse, nous présentons également des résultats préliminaires pour la mise en place d'une technique dimensionnelle adaptative plus connue sous le nom de DAM (Dimensionally Adaptive Method (DAM)). Lorsque la distribution de la résistivité du sous-sol dépend de multiples variables spatiales, une analyse correcte de la dimensionnalité (Ledo 2005, Martí et al. 2009, Weaver and Agarwal 2000) rend parfois possible de considérer les différentes régions avec des dimensions spatiales différentes. Par exemple, il est parfois possible d’interpréter la distribution comme une formation unidimensionnelle plus quelques hétérogénéités en 2D (ou 3D). Basée sur cette interprétation, la DAM tire profit d’une telle situation. Ainsi, l'idée principale de cette méthode est d'effectuer l'adaptativité sur la dimension spatiale en commençant par un problème de faible dimension et en utilisant les résultats obtenus pour minimiser le coût des problèmes de dimension supérieure. Nous commençons l'inversion avec un modèle 1D. Les résultats de ce problème d'inversion 1D sont utilisés comme information a priori sur les modèles de dimension supérieure. Un avantage fondamental de cette approche est que nous pouvons utiliser les solutions des problèmes de dimension inférieure précédemment calculées comme composantes du terme de régularisation associé à un problème de dimension supérieure afin d'augmenter la robustesse de l'inversion. Cette thèse propose également une analyse numérique rigoureuse de divers aspects des problèmes MT. En particulier, nous avons: (a) étudié l'effet de la source, (b) effectué une analyse fréquentielle de sensibilité, (c) illustré l'augmentation du taux de convergence lorsque l'adaptativité hp est employée, (d) séparé les effets 1D et 2D dans la solution numérique et (e) exploré l'intérêt de considérer différentes variables pour effectuer l'inversion. / The magnetotelluric (MT) method is a passive exploration technique that aims at estimating the resistivity distribution of the Earth's subsurface, and therefore at providing an image of it. This process is divided into two different steps. The first one consists in recording the data. In a second step, recorded measurements are analyzed by employing numerical methods. This dissertation focuses in this second task. We provide a rigorous mathematical setting in the context of the Finite Element Method (FEM) that helps to understand the MT problem and its inversion process. In order to recover a map of the subsurface based on 2D MT measurements, we employ for the first time in Mts a multi-goal oriented self adaptive hp-Finite Element Method (FEM). We accurately solve both the full formulation as well as a secondary field formulation where the primary field is given by the solution of a 1D layered media. To truncate the computational domain, we design a Perfectly Matched Layer (PML) that automatically adapts to high-contrast material properties that appear within the subsurface and on the air-ground interface. For the inversion process, we develop a first step of a Dimensionally Adaptive Method (DAM) by considering the dimension of the problem as a variable in the inversion. Additionally, this dissertation supplies a rigorous numerical analysis for the forward and inverse problems. Regarding the forward modelization, we perform a frequency sensitivity analysis, we study the effect of the source, the convergence of the hp-adaptivity, or the effect of the PML in the computation of the electromagnetic fields and impedance. As far as the inversion is concerned, we study the impact of the selected variable for the inversion process, the different information that each mode provides,and the gains of the DAM approach.

Ce travail de thèse a concerné la modélisation mathématique et l’approximation numérique de l’influence d‘une inclusion très fine sur un substrat élastique de différente rigidité. L’étude est motivée par les applications dans les pneumatiques et ne se base pas sur des techniques d’homogénéisation classiques. En effet, l’objectif a été de traiter l’interaction entre une seule inclusion et son milieu élastique et non une densité d’inclusions. L’étude a comporté trois volets, le premier concernant une modélisation mathématique pour des lois de comportement linéaires aboutissant à une expression de la contribution de l’inclusion sous la forme du champ sans inclusion corrigé par des correcteurs à différents ordres. Ces correcteurs sont indépendants de la taille caractéristique de l’inclusion, Le second a concerné l’approximation numérique de cette influence moyennant la méthode des éléments finis et celle des éléments finis inversés. Une stratégie numérique de prise en compte de l’influence de plusieurs inclusions y est aussi présentée. Le dernier volet est prospectif et discute de la possibilité de l’extension de l’approche pour des lois de comportement non linéaires. / This work focused on mathematical modeling and numerical approximation of the influence of a very thin inclusion on an elastic substrate of different stiffness. The study is motivated by applications in tires and is not based on conventional homogenization techniques. Indeed, the objective was to treat the interaction between a single inclusion and its elastic medium and not a density of inclusions. The study consisted of three parts, the first concerning mathematical modeling for linear behavior laws leading to an expression of the contribution of the inclusion in the form of the inclusion-free field corrected by correctors at different orders. These correctors are independent of the characteristic size of the inclusion. The second relates to the numerical approximation of this influence by means of the finite element method and that of the inverted finite elements. A numerical strategy for taking into account the influence of several inclusions is also presented. The last part is prospective and discusses the possibility of extending the approach for nonlinear behavioral laws.