Dans le présent travail, nous mettons en œuvre un développement numérique d’une méthode de discrétisation mixte MEF/Meshless pour la résolution de problème mécanique fortement non-linéaire. Une attention particulière est attribuée à la construction des fonctions de forme par approximation diffuse. Dans le but de traiter des problèmes de la mécanique des solides en transformations finies, nous développons une méthode numérique dite « mixte » unissant à la fois la méthode numérique Meshless afin de discrétiser les zones à fort gradient de déformation, et la méthode des Eléments Finis (MEF) pour les zones les moins sollicitées. Nous veillons donc à assurer le couplage entre ces deux méthodes de discrétisation à travers la mise en œuvre de techniques spécifiques, notamment pour la continuité du domaine et sa consistance en premier lieu, mais aussi pour le choix de la méthode d’intégration numérique appropriée. Nous testons par la suite la fiabilité et la performance du modèle par la mise en place d’une étude comparative avec un modèle MEF standard conçu pour les mêmes conditions numériques, afin de tester la performance des techniques numériques attribuées au modèle « mixte ». Nous développons, par la suite, un modèle de comportement à travers une description des transformations finies. Nous adoptons ainsi une discrétisation spatiale en éléments «assumed strain», et une discrétisation temporelle adaptées. Pour valider notre modèle « mixte » retenu, nous réalisons une étude comparative avec des exemples simples de calcul non linéaire / In this work, we implement a development of mixed discretization MEF-Meshless for solving strongly nonlinear mechanical problem. Particular attention is given to the construction of the shape functions by diffuse approximation. In order to deal with problems of solid mechanics with large strain, we develop a so-called “mixed” numerical method combining both Meshless to discretize areas of high deformation gradient, and Finite Element Method (MEF) for non-concerned areas. We ensure coupling between both discretization methods through the implementation of specific techniques, including the continuity and consistency of the field and the choice of the appropriate method of numerical integration. We test later the reliability and performance of the model by the introduction of a comparative study with a standard FEM model designed for the same numerical conditions to evaluate numerical techniques attributed to our “mixed” model. We develop, thereafter, a model of behavior through a large strain description. We adopt spatial discretization elements “assumed strain” and a suitable time discretization. To validate our “mixed” model, we perform a comparative study of nonlinear simple calculation examples

Les modèles de comportement mécanique, dits macroscopiques, sont développés à la fois pour leur légèreté, permettant le calcul d’éléments structuraux pouvant atteindre d’importantes dimensions, et pour leur finesse de représentation des phénomènes mécaniques observés par le matériau à des échelles plus fines. Le développement de tels modèles est ici effectué dans le cadre de la sollicitation sismique, donc des chargements cycliques alternés, appliquée à des ouvrages en matériaux quasi-fragiles, et plus précisément en béton. À ce jour, les modèles macroscopiques, effectivement applicables au calcul de structures, et représentatifs du comportement cyclique du béton sont encore rares. En conséquence de la complexité du problème de fissuration à homogénéiser, les modèles macroscopiques existants affichent une robustesse limitée ou ne permettent pas de reproduire l’ensemble des phénomènes mécaniques observés par le matériau. Une des barrières à la résolution de ces deux problématiques est le manque de données expérimentales relatives aux phénomènes à modéliser. En effet, en cause de la difficulté technique de les réaliser, peu de résultats d’essais cycliques alternés sur du béton sont disponibles dans la littérature.
 Une démarche d’expérimentation numérique a donc été élaborée sur la base d’un modèle fin du matériau, dit microscopique, capable de fournir les résultats nécessaires à la formulation et à l’identification d’un modèle macroscopique. Dans le modèle microscopique le matériau est considéré comme une structure à part entière, il a été développé afin de ne nécessiter qu’une quantité réduite de résultats d’essais, maîtrisés, pour être mis en oeuvre. Le modèle microscopique, un modèle particulaire lattice, a été développé sur la base d’un modèle lattice existant, enrichi pour être en mesure de simuler le comportement des matériaux quasi-fragiles sous chargements multi-axiaux et cycliques. Le modèle microscopique a alors été validé en tant qu’outil d’expérimentation numérique, et exploité afin d’établir les équations constitutives du modèle macroscopique fondées sur les théories de l’endommagement et de la plasticité. La régularité de la relation de comportement proposée, intégrant un effet unilatéral progressif, a notamment été garantie par l’utilisation d’un modèle d’élasticité non-linéaire. Le modèle macroscopique a finalement été calibré, entièrement, à l’aide du modèle microscopique, et mis à l’oeuvre dans la simulation de la réponse d’un voile en béton armé soumis à un chargement de cisaillement cyclique alterné. Cette simulation a permis de mettre en avant la robustesse numérique du modèle développé, ainsi que la contribution significative du comportement uni-axial cyclique alterné du béton à l’amortissement de telles structures. / Macroscopic mechanical behavior models are developed for their light computational costs, allowing the simulation of large structural elements, and the precise description of mechanical phenomena observed by the material at lower scales. Such constitutive models are here developed in the seismic solicitation framework, therefore implying cyclic alternate loadings at the material scale, and applied to civil engineering buildings, often made of concrete, or more generally of quasi-brittle materials. To date, macroscopic models applicable to structural computations, while representing the cyclic mechanical behavior are rare. In consequence of the intricacy of the fracture processes to homogenize, macroscopic constitutive models either do not present sufficient robustness or miss on important phenomena. One of the limitations to the resolution of this issue is the lack of experimental data. Indeed, because of the complexity of the experiments to set up, few results on alternate cyclic tests on concrete are available in the literature.A virtual testing approach has therefore been established on a microscopic model of the material, able to provide results needed to the formulation and the calibration of a macroscopic model. In the microscopic model, the material is considered as structure itself, it is developed so as to only necessitate a reduced amount of results from controlled experimental tests, in order to be used. The microscopic model, a lattice discrete element model, has been developed on the basis of an existing lattice model and extended to the simulation of multi-axial and cyclic loadings. The microscopic model has then been validated as a virtual testing tool and used to establish equations of the macroscopic model, on the basis of damage and plasticity theories. The consistency of the proposed constitutive relation, embedding progressive unilateral effect, has been achieved using non-linear elasticity. The macroscopic model has finally been calibrated, entirely with the microscopic model, and employed to simulate the response of a reinforced concrete wall under alternate shear loading. This simulation has served to showcase the numerical robustness of the proposed model, as well as the significant contribution of the uni-axial alternate behavior of concrete to the structural damping of such structures.

La hauteur des arbres est un trait écologique majeur représentant l'intensité de la compétition pour la lumière. De plus, la croissance des arbres est le résultat de multiples compromis afin de maintenir leur orientation verticale et leur stabilité mécanique tout en assurant les autres fonctions écophysiologiques. Le contrôle de l'orientation de la croissance est réalisé par deux mécanismes: la croissance différentielle au niveau du méristème apical et la formation de bois de réaction au cours de la croissance secondaire. Cependant, la modélisation simultanée de la croissance et des rétroactions biomécaniques dépasse le cadre classique de la mécanique des structures. En effet, le concept de configuration de référence devient imprécis dû à l'apparition de nouveaux points matériels libres de contraintes sur une surface déformée au cours de la croissance. Dans cette thèse, un nouveau formalisme mathématique est proposé à partir de la théorie des poutres et modélise simultanément les effets de la croissance et de la biomécanique de l'arbre. Afin de résoudre le système d'équations aux dérivées partielles, de nouvelles méthodes numériques sont développées et tiennent compte de la dépendance entre l'espace et le temps. Enfin, deux problèmes de contrôle optimal sont analysés, modélisant les stratégies d'allocation dynamique de la biomasse pour la croissance primaire et secondaire, en fonction de différents contextes écologiques. Ce travail offre de nouvelles perspectives sur les mathématiques de la mécanique de la croissance et ses applications en biologie. / Height is a major ecological trait for growing trees, representing the intensity for light competition. Moreover, tree height results from a trade-off between different functions, including tree mechanical stability. Trees develop growth strategies to maintain their vertical orientation and mechanical stability, in addition to other ecophysiological functions, through differential primary growth near the shoot apical meristem and formation of reaction wood during secondary growth. However, this coupling is a problematic issue since the progressive addition of new material on an existing deformed body makes the definition of a reference configuration unclear. This thesis presents a new mathematical framework for rod theory modelling simultaneously the interactions between the growth process and tree biomechanics. In order to solve the obtained system of partial differential equations, new numerical methods are developed and take into account the dependence between space and time, which is a specific feature of surface growth problems. Finally, the present work addresses the mathematical formulation of two optimal control problems characterising tree's growth strategies. Growth strategies are analysed with respect to the ecological context, through two variables, which are the ratio of biomass allocated to primary growth and the distribution of biomass allocated to secondary growth along the growing stem. This work gives new insights to the mathematical framework of surface growth mechanics and its applications in biology.

Dans cette thèse, on utilise des outils mathématiques et numériques pour modéliser les écoulements tridimensionnels incompressibles à surface libre instationnaires. L'application industrielle visée est l'étude de la phase de remplissage des moules dans une fonderie. On développe un algorithme pour le transport de l'interface par la vitesse du fluide pour un fluide diphasique incompressible visqueux non-miscible de rapport de densité important en utilisant la méthode de caractéristiques pour traiter le problème de convection. Il y a des défis majeurs dans le contexte de la modélisation des fluides diphasiques. Tout d'abord, on doit prendre en considération l'évolution de l'interface et de ses changements topologiques. Deuxièmement, on doit traiter la non-linéarité convective de l'interface et de l'écoulement. Troisièmement, les équations de Navier-Stokes et du transport doivent être munies des conditions aux bords appropriées. En outre, il faut traiter soigneusement les singularités géométriques et topologiques à travers l'interface en particulier dans le cas de rapport de densité et viscosité important. On doit également maintenir la résolution d'une interface d'épaisseur nulle durant les cas du pliage, la rupture et la fusion de l'interface. Quatrièmement, on doit respecter les propriétés physiques telles que la conservation de la masse pour tout écoulement d'un fluide incompressible. Cinquièmement, il faut toujours penser aux limitations du temps de calcul et de mémoire pour résoudre ce genre de problème dans les cas pratiques. Notre but est de trouver un schéma fiable capable de modéliser le remplissage des moules tridimensionnelles industrielles. La première partie de cette thèse est dédiée à la description mathématique du schéma de transport de l'interface par la vitesse du fluide. Le mouvement des fluides est décrit par les équations de Navier-Stokes. L'interface est capturée par la fonction Level-Set. Le problème est discrétisée en espace par la méthode des éléments finis et en temps par la méthode de caractéristiques.Des conditions aux bords appropriées pour le problème du remplissage d'un moule sont introduites et un algorithme de calcul de la solution est présentée. Finalement,des résultats numériques montrent et valident l'efficacité duschéma proposé. Dans la deuxième partie de cette thèse, on introduit une méthode de décomposition de domaine qui correspond à la discrétisation par la méthode des caractéristiques dans le but d'améliorer la performance de l'algorithme proposé lors de la modélisation du remplissage des moules industrielles à moyennes séries. Des résultats numériques de comparaison valident la précision du code parallèle. / In this thesis, we use mathematical and numerical tools to model three dimensional incompressible laminar flows with free surface. The described industrial application is the study of the mould filling phase in a foundry. We develop an algorithm for the transport of the interface by the fluid velocity for a viscous incompressible immiscible fluids of high density ratio in two-phase flow using the method of characteristics for the convection problem.There are, however, major challenges in the context of two-phase flow modeling.First, we have to take into account the evolution of the interface and its topological changes. Second, we have to deal with the non-linearity for the convection of the flow and the interface. Third, we must assign appropriate boundary conditions to the flow and transport equations.In addition, care must be taken in treating the geometrical and topological singularities across the interface.We also have to maintain a sharp interface resolution, including the cases of interface folding, breaking and merging.Furthermore, we should respect the physical properties such as the mass conservation for any incompressible fluid flows.Finally, we have to keep in mind the limitations in the time of computation and memory to solve this kind of problem in practical cases. Our purpose is to find a reliable scheme able to model the filling of three dimensional industrial moulds.The first part of the thesis is devoted to the mathematical description of the interface transport scheme by the fluid velocity. The fluids motion is described by the Navier-Stokes equations. The interface is captured by the Level-Set function. The problem isdiscretized by the characteristics method in time and finiteelements method in space. The interface is captured by the Level-Setfunction. Appropriate boundary conditions for the problem ofmould filling are investigated, a new natural boundary conditionunder pressure effect for the transport equation is proposed andan algorithm for computing the solution is presented. Finally,numerical experiments show and validate the effectiveness of theproposed scheme.In the second part of the thesis, we introduce a domain decomposition method that suits the discretization by the method of characteristics in order to improve the performance of the proposed algorithm to model the filling phase for moulds of medium series. Numerical results of comparison validate the precision of the parallel code.

Cette thèse s’intéresse à la modélisation et à la commande de robots déformables, c’est à dire de robots dont le mouvement se fait par déformation. Nous nous intéressons à la conception de lois de contrôle en boucle fermée répondant aux besoins spécifiques du contrôle dynamique de robots déformables, sans restrictions fortes sur leur géométrie. La résolution de ce défi soulève des questions théoriques qui nous amènent au deuxième objectif de cette thèse: développer de nouvelles stratégies pour étudier les systèmes de grandes dimensions. Ce manuscrit couvre l’ensemble du développement des lois de commandes, de l’étape de modélisation à la validation expérimentale. Outre les études théoriques, différentes plateformes expérimentales sont utilisées pour valider les résultats. Des robots déformables actionnés par câble et par pression sont utilisés pour tester les algorithmes de contrôle. A travers ces différentes plateformes, nous montrons que la méthode peut gérer différents types d’actionnement, différentes géométries et propriétés mécaniques. Cela souligne l’un des intérêts de la méthode, sa généricité. D’un point de vue théorique, les systèmes dynamiques à grande dimensions ainsi que les algorithmes de réduction de modèle sont étudiés. En effet, modéliser des structures déformables implique de résoudre des équations issues de la mécanique des milieux continus, qui sont résolues à l’aide de la méthode des éléments finis (FEM). Ceci fournit un modèle précis des robots mais nécessite de discrétiser la structure en un maillage composé de milliers d’éléments, donnant lieu à des systèmes dynamiques de grandes dimensions. Cela conduit à travailler avec des modèles de grandes dimensions, qui ne conviennent pas à la conception d’algorithmes de contrôle. Une première partie est consacrée à l’étude du modèle dynamique à grande dimension et de son contrôle, sans recourir à la réduction de modèle. Nous présentons un moyen de contrôler le système à grande dimension en utilisant la connaissance d’une fonction de Lyapunov en boucle ouverte. Ensuite, nous présentons des algorithmes de réduction de modèle afin de concevoir des contrôleurs de dimension réduite et des observateurs capables de piloter ces robots déformables. Les lois de contrôle validées sont basées sur des modèles linéaires, il s’agit d’une limitation connue de ce travail car elle contraint l’espace de travail du robot. Ce manuscrit se termine par une discussion qui offre un moyen d’étendre les résultats aux modèles non linéaires. L’idée est de linéariser le modèle non linéaire à grande échelle autour de plusieurs points de fonctionnement et d’interpoler ces points pour couvrir un espace de travail plus large. / This thesis focuses on the design of closed-loop control laws for the specific needs of dynamic control of soft robots, without being too restrictive regarding the robots geometry. It covers the entire development of the controller, from the modeling step to the practical experimental validation. In addition to the theoretical studies, different experimental setups are used to illustrate the results. A cable-driven soft robot and a pressurized soft arm are used to test the control algorithms. Through these different setups, we show that the method can handle different types of actuation, different geometries and mechanical properties. This emphasizes one of the interests of the method, its genericity. From a theoretical point a view, large-scale dynamical systems along with model reduction algorithms are studied. Indeed, modeling soft structures implies solving equations coming from continuum mechanics using the Finite Element Method (FEM). This provides an accurate model of the robots but it requires to discretize the structure into a mesh composed of thousands of elements, yielding to large-scale dynamical systems. This leads to work with models of large dimensions, that are not suitable to design control algorithms. A first part is dedicated to the study of the large-scale dynamic model and its control, without using model reduction. We present a way to control the large-scale system using the knowledge of an open-loop Lyapunov function. Then, this work investigates model reduction algorithms to design low order controllers and observers to drive soft robots. The validated control laws are based on linear models. This is a known limitation of this work as it constrains the guaranteed domain of the controller. This manuscript ends with a discussion that offers a way to extend the results towards nonlinear models. The idea is to linearize the large-scale nonlinear model around several operating points and interpolate between these points to cover a wider workspace.

Les géo-structures telles que les barrages et les digues sont soumises à des écoulements hydrauliques variant dans le temps et dans l'espace. L'eau qui traverse ces milieux poreux peut entraîner le détachement et le transport de certaines particules des sols constituant les structures et leurs fondations. Ce problème est généralement appelé "érosion interne". Le terme suffusion, un type d'érosion interne, se réfère au détachement et au transport de particules les plus fines à traversune matrice de sol poreuse plus grossière en raison d'un écoulement hydraulique. L'évolution temporelle de la suffusion peut modifier les propriétés hydrauliques et mécaniques des sols et peut entraîner des changements importants dans le comportement de telles structures pouvant aller jusqu’à leur effondrement. Ce travail de thèse tente de contribuer à la conception et à la durabilité des ouvrages en ingénierie géotechnique et hydraulique en mettant un accent particulier sur les barrages, les levées et les digues. Il a été consacré à développer un modèle numérique de suffusion en introduisant d’une part le couplage des phénomènes hydrauliques et mécaniques et d’autre part le couplage des phénomènes d'érosion et de filtration. / Geo-structures such as dams and levees or dikes are subjected to seepage varying in time and space. The water flowing through these porous media can lead to the detachment and transport of part of the soil particles within the structures or their foundations. This problem is usually called internal erosion. The term suffusion, one type of internal erosion, refers to the detachment and transport of finer particles through a coarser porous soil matrix due to seepage flow. Suffusion can modify with time the hydraulic and mechanical properties of the soils and may trigger significant damage on such structures and lead eventually to their collapse. This research attempts to contribute to the design and sustainability of geotechnical and hydraulic engineering structures, with a particular focus on embankment dams, levees, and dikes. It aims to develop a numerical model of suffusion by introducing, on the one hand, the coupling of the hydraulic and mechanical phenomena and, on the other hand, the coupling of erosion and filtration.

Les constructeurs automobiles utilisent des aciers auto-trempant à haute résistance au bore manganèse (22MnB5) pour des questions de sécurité et environnementale. La mise en forme à chaud de tôles constituées de cet acier impose aux outillages des contraintes mécaniques sévères et un échauffement intense. L'intégration d'une fonction de cisaillage sous presse à suivre permet de répondre à des problématiques rentabilité de production, mais les lames sont confrontées à des problèmes similaires de durabilité. L'objectif de cette thèse est de comprendre la dégradation des outils de cisaillage à chaud afin de proposer des guides de choix de matériaux d'outil (X38CrMoV5-3 ou X70CrMoV5-2) et de paramètres du procédé. Cela requiert l'estimation des sollicitations mécaniques et thermiques dans la partie active de l'outil. La méthodologie développée passe par quatre étapes : une recherche bibliographique, des essais de cisaillage à chaud sur un module spécifique du pilote MEFISTO de l'ICA, le développement d'un modèle thermomécanique éléments finis de l'opération et des analyses des microstructures. Un état de l'art portant sur la définition de l'opération de découpage de tôles, sur les propriétés physiques et mécaniques du matériau de tôle et des lames, le comportement à l'interface tôle/outil et sur les techniques permettant de simuler le découpage a été proposé. Des essais sur le module de découpage ont permis d'accéder à des résultats d'effort et des observations des endommagements des lames de cisaillage. Ces informations sont utilisées pour valider le modèle de calcul éléments finis. Ce modèle montre que les contraintes mécaniques sont locales, intenses et se déplacent dans l'arête de coupe. De plus, l'échauffement de l'outil engendre une température de surface proche des températures de revenu des matériaux d'outil. La comparaison des résultats de simulation numérique et de l'étude expérimentale permettent de comprendre le lien entre les dégradations observées et les sollicitations thermomécaniques. / Car manufacturers are building using a self-hardening high strength boron steel (22MnB5) for safety and environmental concerns. But the shaping of sheets made of these kind of steel imposes on tools severe mechanical stresses and strong thermal transfer. The integration of a blanking function in transfer presses is a response to reach high profitability, but blades are subjected to the same durability problems. The aim of this thesis is to understand how hot blanking tools are wearing in order to give to manufacturers tool material guidelines (X38CrMoV5-3 or X70CrMoV5-2) and process parameters guidelines. To do so, an estimation of mechanical stresses and thermal solicitation is required. The methodology is based on four steps: literature search, hot blanking trials on a specific module of ICA laboratory industrial pilot MEFISTO, finite element simulations of the operation, and microstructural analysis. The state of art is focused on the description of the operation, the physical and mechanical properties of the blank material and tools materials, the blank/tool interface behavior, and the techniques used simulate the hot blanking process. Trials on the hot blanking module provided blanking force and observations of worn blades. These informations are needed to validate the finite element model. This model shows that mechanical stresses are local, intense and slides in the blade cutting edge. Moreover, the blade heating generates a surface temperature close to tool steel tempering temperature. Experimental results compared to numerical ones allow to understand the link between thermo-mechanical stresses and how the damage occurs.

La compréhension du fonctionnement du cerveau est un des défis majeurs des neurosciences. Pour appréhender cet organe in vivo de nombreux dispositifs sont développés, parmi eux l'électroencéphalographie (EEG). L'EEG mesure directement et d’une manière non invasive les signaux du cerveau avec une haute résolution temporelle. L’activation électrique d’une région dans le cerveau est modélisée par un dipôle de courant. Le principal objectif de l’EEG est le suivi de l'activité cérébrale pendant une tâche spécifique et la localisation/reconstruction des dipôles sources dans le cerveau. Il est utilisé à la fois dans la pratique et dans les recherches en neuroscience. Le processus de localisation des sources implique la résolution d'un problème inverse qui nécessite des modèles de la solution du problème direct. Le problème direct concerne la prédiction du potentiel électrique, résultant de l'activité cérébrale, sur le scalp à partir d'une répartition des sources données dans le cerveau. La précision de la localisation des sources dépend largement des performances de la solution du problème direct, qui est liée à la précision du modèle de tête, du modèle des sources ainsi que de la méthode de résolution du problème direct. Actuellement en pratique, pour des raisons de simplicité de calcul et de rapidité, les algorithmes EEG utilisent principalement le modèle de tête sphérique multicouche ou la méthode des éléments de frontière pour résoudre le problème direct. L’amélioration de la résolution en EEG nécessite l’utilisation de modèles de têtes plus réalistes qui prennent en considération d’avantages de paramètres tels que l’inhomogénéité et l’anisotropie des tissus. C’est pour ces raisons que la méthode des éléments finis (FEM) est la mieux adaptée et a attiré l’attention de plusieurs chercheurs. Néanmoins, avec la discrétisation de la FEM, les sources peuvent présenter une singularité numérique, ceci impacte négativement la solution du problème direct. Afin de remédier à ce problème des techniques tel que la méthode directe, la méthode de soustraction et la méthode de Saint Venant sont utilisées pour traiter la singularité. Cependant, toutes ces méthodes montrent des instabilités numériques dans le cas de sources situées à proximité des interfaces des tissus du cerveau. Pour y remédier une amélioration a été apportée à la méthode de Saint Venant au cours de cette thèse. Le rapport donne un aperçu sur l’activité cérébrale, un rappel de la FEM, son application pour la résolution du problème direct en EEG, le traitement de la singularité numérique et des instabilités près des interfaces, puis une présentation de la version modifiée de Saint-Venant, proposée et mise en œuvre dans ce travail. Nous validerons les résultats de la méthode modifiée de Saint-Venant dans des modèles de têtes sphériques multicouches et des modèles à géométries réelles. Pour finir, on testera cette méthode dans les outils et logiciels actuellement utilisés dans la pratique pour la localisation des zones activé dans le cerveau et nous montrerons les améliorations qui peuvent être apportées en utilisant notre approche de la méthode de Saint Venant modifiée. / Understanding the brain functions is one of the major challenges of neuroscience. To apprehend this organ in vivo many devices are developed, among them the electroencephalography (EEG). EEG measures directly and noninvasively the brain signals with high time resolution. The electrical activation of a region in the brain is modeled by a current dipole. The main objective of the EEG is monitoring the brain activity during a specific task and the location/reconstruction of the dipole sources within the brain. It is used both in practice and in research in neuroscience. The source’s localization process consists of solving an inverse problem that requires models of the solution of the direct or forward problem. The forward problem is the prediction of the electric potential on the scalp from a distribution sources in the brain. The accuracy of the source’s localization depends largely on the performance of the forward problem’s solution, which is related to the accuracy of the head model, the source model and the related method used to solve the forward problem. Currently in practice for reasons of computational simplicity and speed, EEG codes use mainly the multilayered spherical head model or the boundary element method for solving the forward problem. Improving the EEG’s source localization requires the use of more realistic head models which take into account more parameters such as inhomogeneity and anisotropy of the tissues. For these reasons the finite element method (FEM) is best suited and attracted the attention of many researchers. Nevertheless, with the FEM discretization, the sources may bring a numerical singularity; this negatively impacts the solution of the forward problem. To treat this problem some techniques such as the direct method, the subtraction method and the method of Saint-Venant are developed. However, all these methods show numerical instabilities in the case from sources close to the interfaces of brain tissue. To remedy these instabilities a modification was made on the Saint-Venant’s method. The report provides an overview on brain activity, a reminder of the FEM, its application for solving the forward problem in EEG, processing to treat the numerical singularity of the source and the instability near interfaces using the modified method of the Saint-Venant. We validate the results of the modified Venant’s method in models of multilayered spherical head and models with real geometries. Finally, we will test this approach in tools and software currently used in practice for locating areas activated in the brain and shows the improvements that can be made using the method of modified Saint-Venant.

Le travail présenté dans ce mémoire s’intéresse à la prédiction des vibrations d'un noyau de transformateur multicouche, constitué d'un assemblage de tôles ferromagnétiques. Le problème couplé magnéto-mécanique est résolu par une approche séquentielle progressive : la résolution magnétique est suivie d'une résolution mécanique. Un modèle multi-échelle simplifié 3D décrivant les anisotropies magnétiques et magnétostrictives, et considérant les non-linéarités magnétiques et de magnétostriction, est utilisé comme loi de comportement du matériau. La structure du noyau du transformateur est modélisée en 2D. Une technique d'homogénéisation permet de tenir compte du comportement anisotrope de chaque couche afin de définir un comportement moyen pour chaque élément du maillage éléments finis.. Des mesures expérimentales sont ensuite effectuées, permettant d’une part la validation des lois de comportement matériau utilisées, et d’autres part des modèles de comportement structurel statique, du comportement structurel dynamique et de l'estimation du bruit. Différents matériaux et différentes géométries de prototypes de transformateurs sont considérés pour ce travail. Des optimisations structurelles sont finalement proposées grâce à des simulations numériques s’appuyant sur le modèle développé, afin de réduire les vibrations et les émissions de bruit du noyau du transformateur. / This thesis deals with the prediction of the vibration of a multi-layer transformer core made of an assembly of electrical sheets. This magneto-mechanical coupled problem is solved by a stepping finite element method sequential approach: magnetic resolution is followed by mechanical resolution. A 3D Simplified Multi-Scale Model (SMSM) describing both magnetic and magnetostrictive anisotropies is used as the constitutive law of the material. The transformer core structure is modeled in 2D and a homogenization technique is implemented to take the anisotropic behavior of each layer into consideration and define an average behavior at each element of the finite element mesh. Experimental measurements are then carried out, allowing the validation of the material constitutive law, static structural behavior, dynamic structural behavior, and the noise estimation. Different materials geometries are considered for this workStructural optimizations are finally achieved by numerical simulation for lower vibration and noise emission of the transformer cores.

Ce mémoire entame le développement d’une méthode d’évaluation d’assemblages en aluminium liée par adhésifs structuraux dans un contexte d’ingénierie du secteur industriel. Les travaux ont pour but de fournir à l’entreprise Services Précicad inc. des outils de conception applicables sur le véhicule électrique Kargo. Plus spécifiquement, des séries d’essais expérimentaux sont réalisées sur des échantillons ASTM au moyen de plans d’expériences employant la méthode Taguchi pour étudier l’incidence des conditions de préparation manufacturière sur la résistance mécanique et pour faire un dépistage de produits commerciaux. La famille des méthacrylates est alors privilégiée pour répondre aux besoins de flexibilité d’utilisation et de préparation de surface. Ensuite, différents modèles de calculs analytiques sont étudiés et une analyse paramétrique est réalisée sur le joint en cisaillement à simple recouvrement (SLS). Le modèle du joint sandwich simplifié de Bigwood et Crocombe est choisi comme outil de calculs de prédimensionnement de joints génériques, cependant une bonne connaissance des conditions limites est nécessaire. L’approche de simulation par la méthode des éléments finis en mécanique des milieux continus est alors recommandée. D’importantes difficultés de convergence des résultats des contraintes sont identifiées avec cette théorie, entre autres par la présence de singularités numériques aux interfaces du volume de l’adhésif. Des techniques de modélisation simplifiées, utilisant les outils du logiciel NX Nastran, sont ensuite testées et l’outil Surface-Gluing est alors proposé pour permettre la simulation de joints qui reflètent les caractéristiques de modélisation « industrielles » de Précicad. Finalement, une méthode de mesure des propriétés élastiques in situ pour adhésifs flexibles, compensée au moyen d’analyses par éléments finis, est développée afin de pouvoir réaliser une caractérisation avec un extensomètre conventionnel, tout en préservant la précision souhaitée. Ainsi, il a été possible d’isoler les modules de cisaillement et de tension de l’adhésif H8000. / This thesis initiates the development of a structural assessment method for aluminium adhesively bonded joints, intended for the industrial engineering businesses sector. The present work is devoted to the company Services Précicad inc, who needs design tools that can be used for its Kargo electric vehicle assemblies. In particular, the Taguchi design of experiment method is used to study the effects of joint conditioning on the static strength of different ASTM samples configurations. Commercially available products are screened out during the process. Thus, methacrylates is found to be the best adhesive chemical category for Précicad design and manufacturing needs. Subsequently, a review of different analytical models is carried and the Bigwood & Crocombe sandwich model is selected to perform a parametric analysis on the Single Lap Shear joint (SLS). The simplified version of the sandwich model was found to be the best suited for pre-sizing generic forms of adhesive joints, although a good knowledge of the boundary conditions is mandatory to make a realistic assessment. Consequently, the standard continuum mechanics finite element method is recommended for most joints analysis. However, major convergence difficulties can be found when using this theory, partly due to the presence of numerical singularities in the adhesive domain. Several simplefied modelling techniques using NX Nastran are tested on the SLS joint and the Surface-Gluing tool is proposed for evaluating complex “industrial” meshes such as the ones used by Precicad. Finally, a method for measuring in situ elastic properties of flexible adhesives with finite element analysis compensation is developed for achieving a complete characterization with a conventional clip gauge. Thus, it has been possible to measure the shear and tensile modulus of the H8000 adhesive.