Dans le cadre de cette thèse, une modélisation par éléments finis continue renforcée discrète est développée. Le modèle est basé sur le concept de la cellule élémentaire, avec utilisation de connecteurs pour prendre en compte la rigidité en tension de la structure, et modéliser l'anisotropie du tissu et un élément coque pour décrire la rigidité en cisaillement et gérer le frottement et le contact avec les outils lors de la simulation du procédé de préformage. Les différents éléments de ce maillage spécifique sont générés automatiquement à l'aide des scripts pythons. Associés à ce choix de modélisation, les caractéristiques du comportement sont déterminées à l'aide d'essais expérimentaux conduits sur le renfort et d'une démarche de méthode inverse. Le comportement de ces différents éléments est choisi d'être linéaire élastique en premier lieu, et par la suite à évoluer pour devenir non linéaire. Ces développements numériques sont réalisés dans le logiciel Abaqus/ Explicit. Dans un second temps le modèle développé est utilisé dans le cadre de simulation de l'étape de préformage des renforts secs. A ce titre un démonstrateur expérimental développé dans le cadre de la thèse sert à obtenir des résultats expérimentaux pour corréler ces travaux de simulations. Des comparaisons avec des simulations mono-couche, hémisphérique, permettre de mettre en évidence l'influence des paramètres de simulations mises en jeu. Dans le cadre de préformages de plusieurs plis, l'influence des caractéristiques du contact / frottement est mise en évidence. Finalement une procédure d'automatisation afin de rendre la simulation plus ergonomique dans un contexte industriel est présentée. Cette procédure consiste à générer le maillage spécifique et mettre la simulation en donné dans Catia. Les simulations peuvent ainsi se réaliser sans le passage par Abaqus/CAE. / In this work, a discrete approach for the simulation of the preforming of dry woven reinforcementis proposed. A “unit cell” is built using elastic isotropic shells and axial connectors instead of bars and beams used in previous studies. Shell elements are used to take into account the in-plane shear stiffness and to manage contact phenomenon with the punch and die. Connectors reinforce thestructure in the yarn directions and naturally capture the specific behavior of the fabric. To identify the material parameters, uniaxial tensile tests and bias tests have been employed. A numerical algorithm,coupling Matlab and Abaqus/Explicit, is used to determine the shear parameters by an Inverse method. The model has been implemented in Abaqus to simulate hemispherical stamping.Experimental results are compared to numerical simulations, good agreement between both resultsis shown for the case of single layer forming.This model has been then used to simulate a multilayer performing and a numerical study of the effect of fiction coefficient has been performedFinally a tool which interfaces a CAD Software to the Abaqus/code has been developed.

La modélisation multi-échelle des matériaux hétérogènes est un challenge en mécanique numérique. Dans le contexte non linéaire, les propriétés effectives des matériaux hétérogènes ne peuvent pas être obtenues par les techniques utilisées pour les milieux linéaires car le principe de superposition n'est plus valable. Ainsi, dans le contexte des éléments finis, une alternative au maillage de l'ensemble de la structure avec la prise en compte de toutes les hétérogénéités, est l'utilisation de la méthode d'éléments finis multi-échelles (EF2). Les techniques de ce type offrent de nombreux avantages, tels que la prise en compte : des grandes déformations au niveau micro et macro sont souvent résolus par les procédures classiques de Newton-Raphson, qui sont généralement adaptées à la résolution des problèmes non linéaires mais qui présentent des difficultés en présence d'instabilités. Dans cette thèse, la combinaison de la méthode des éléments finis multi-échelles (EF2) et la méthode asymptotique numérique (MAN), surnommée MAN multi-échelle, permet de mettre en œuvre une technique numérique efficace pour traiter les problèmes d'instabilités dans le cadre des matériaux hétérogènes. Ces instabilités peuvent survenir à la fois au niveau micro et au niveau macro. Différentes classes de comportement des matériaux ont été implantées dans notre procédure. Pour améliorer le conditionnement du problème multi-échelle à résoudre, une technique d'homogénéisation du second ordre a été également adaptée dans le cadre de la technique MAN multi-échelle. Par ailleurs, afin de réduire le temps de calcul, quelques techniques ont été proposées dans ce travail / The multiscale modelling of the heterogeneous materials is a challenge in computational mechanics. In the nonlinear case, the effective properties of heterogeneous materials cannot be obtained by the techniques used for linear media because the superposition principle is no longer valid. Hence, in the context of the finite element method, an alternative to mesh the whole structure, including all heterogeneities, is the use of the multiscale finite element method (FE2). These techniques have many advantages, such as taking into account : large deformations at the micro and macro scales, the nonlinear constitutive behaviors of the material, and microstructure evolution. The nonlinear problems in micro and macro scales are often solved by the classical Newton-Raphson procedures, which are generally suitable for solving nonlinear problems but have difficulties in the presence of instabilities. In this thesis, the combination of the multiscale finite element method (FE2) and the asymptotic numerical method (ANM), called Multiscale-ANM, allows one to obtain a numerical effective technique for dealing with the instability problems in the context of heterogeneous materials. These instabilities can occur at both micro and macro levels. Different classes of material constitutive relation have been implemented within our procedure. To improve the multiscale problem conditioning, a second order homogenization technique was also adapted in the framework of Multiscale-ANM technique. Furthermore, to reduce the computational time, some techniques been proposed in this work

Les céramiques piézoélectriques, comme le Titatano-Zirconate de Plomb (PZT), peuvent produire une tension électrique quand elles sont soumises à une contrainte mécanique et, inversement, se déforment sous l'effet d'un champ électrique. Ce couplage électromécanique peut être décrit par des équations de comportement linéaires pour des chargements modérés. Cependant, au-delà de certaines valeurs de champ électrique ou de contrainte mécanique, ce couplage devient fortement non linéaire à cause des phénomènes de réorientation de la polarisation électrique. Dans ce travail de thèse, un modèle phénoménologique, qui tient compte des réorientations ferroélectrique (par un champ électrique) et ferroélastique (par une contrainte mécanique) de la polarisation électrique, est proposé. Deux variables internes sont considérées pour décrire l'histoire du chargement et deux surfaces de charges électrique et mécanique sont définies pour déterminer les débuts des écoulements ferroélectrique et ferroélastique. Une version bi-dimensionnelle de ce modèle est développée également pour faire l'étude des structures piézoélectriques minces. Les deux versions 2D et 3D du modèle phénoménologique sont intégrées implicitement en adoptant la méthode de retour radial (prédiction/correction). Deux éléments finis coque et hexaédrique de premier ordre, qui intègrent ce modèle phénoménologique non linéaire, sont ensuite développés et implémentés dans le code de calcul par éléments finis Abaqus via la routine utilisateur UEL (User ELement) / Piezoceramics like lead zirconate titanate or PZT can produce an electric potential when they are subjected to a mechanical stress and deform in the presence of an electric field. This electromechanical coupling can be described by linear constitutive equations for moderate loadings. Nevertheless, this coupling becomes highly non linear when piezoceramics are subjected to high electromechanical loadings due to the electric polarization switching. In this thesis work, a phenomenological material constitutive model that describe the electric polarization ferroelectric switching (by an electric field) and ferroelastic switching (by a mechanical stress) is proposed. To describe the loading history, two internal variables are considered and two electric and mechanical loading surfaces are defined to indicate the onset of domain switchings. A bi-dimensional version of this model is developed to study thin piezoelectric structures. The phenomenological model 2D and 3D versions are implicitly integrated by adopting the return-mapping algorithm. Two shell and hexahedral first-order finite elements are then formulated and implemented into the commercial finite element code Abaqus via the user subroutine UEL (User ELement)

L’objectif de cette thèse est d’appréhender la complexité des mécanismes biologiques de l’acupuncture afin de construire un modèle mathématique multi-échelle. Ce modèle est étudié théoriquement et numériquement. L’acupuncture est une des plus vielles pratiques de l’histoire de la médecine et une partie intégrante de la médecine traditionnelle chinoise. Dans sa pratique la plus classique, une ou plusieurs aiguilles sont placées à des endroits spécifiques, nommés points d’acupuncture. L’aiguille est ensuite manipulée en utilisant des mouvements de rotation et de translation de façon à stimuler le point d’acupuncture. Les effets cliniques de l’acupuncture pourraient être le résultat d’effet de cascades de réactions produites par les interactions entre l’hypoderme et les systèmes nerveux, endocrinien et immunitaire. Le travail présenté s’articule sur la modélisation de l’insertion d’une aiguille dans le tissu conjonctif de l’hypoderme. Un modèle d’écoulement en milieu poreux du liquide interstitiel de l’hypoderme a permis d’étudier numériquement les composantes de contrainte qui agissent sur les récepteurs à la surface des cellules du tissu et notamment des mastocytes.Un modèle mathématique de la réponse chimiotactique des mastocytes à une contrainte physique créée par le traitement d’acupuncture est développé. Ce modèle prend en compte les mécanismes de signalisation cellulaire. La contrainte physique induit la libération rapide et continue, grâce au recrutement chimotactique de mastocytes, d’attractants et de médiateurs chimiques. Le modèle est basé sur le modèle de chimiotaxie de type Keller-Segel. / The objective of this thesis is to comprehend the complexity of the underlying basis of acupuncture. Acupuncture needling is investigated in order to establish a multiscale model that takes into account the complexity of biology but is mathematically simple enough to run simulations.Acupuncture is one of the oldest practices in the history of medicine and is the core of Traditional Chinese Medicine. Once needles are inserted in the right locations, called acupoints, they are manipulated via manual needling to stimulate the acupoint. The physiological reactions of acupuncture needling lead to therapeutic effects which can be explained by a series of interactions between the skin and the nervous, the endocrine, and the immune systems.In the present work, the thrusting and lifting of an acupuncture needle inserted in subcutaneous connective tissue is modeled. A porous media model is used to run simulations and compute the pressure and shear stress affecting the organization of fibers and of isolated cells in their matrix. A mathematical model was conceived to take into account cell signaling. There is ample evidence that needle manipulation in acupuncture can cause degranulation of mastocytes directly through a physical stress to occur. Activated mastocytes rapidly release granules containing chemical mediators. These chemical mediators play a key role recruiting mastocytes in their environment and are known to affect the excitability of nerve endings as well as local microcirculation permeability and size for the appropriate transfer of long-term acting endocrine signals. The process is sustained by the recruitment of mastocytes through chemotaxis.

L’identification de certaines lois de comportement demeure délicate surtout pour des comportements non linéaires. Plusieurs stratégies d’identification ont été développées en s’appuyant principalement sur le dialogue essai-calcul. La partie expérimentale consiste à mesurer une ou plusieurs quantités physiques qui caractérisent le comportement du matériau testé et la partie théorique permet de calculer ces mêmes quantités. Pour remonter aux propriétés matérielles recherchées, un accord entre les quantités mesurées et calculées doit être visé. Dans cette thèse, on s’intéresse au zinc qui un matériau polycristallin. Celui-ci est modélisé à travers des lois de plasticité cristalline qui permettent de prédire son comportement mécanique sous chargements complexes ainsi que l’évolution de sa microstructure. Les lois de plasticité cristalline implémentées dans le logiciel commercial Abaqus par Huang et Marin que nous avons adapté au cas du zinc ont été utilisées. L’étape d’identification de ces lois non linéaires est basée sur une approche inverse qui tient compte des hétérogénéités des déformations à l’échelle des grains. Pour cela, nous avons utilisé l’essai de nanoindentation qui est un essai hétérogène permettant d’extraire plus d’informations que les essais homogènes. A cet effet, et pour caractériser la plasticité cristalline du zinc, nous avons utilisé les courbes charge - profondeur de pénétration et les profils de déplacements mesurés sur l’empreinte résiduelle sur des grains de différentes orientations cristallographiques mesurées par EBSD. La confrontation des résultats expérimentaux et montrent la bonne adéquation entre les données expérimentales et les modèles identifiés. Les résultats obtenus ont ainsi permis de caractériser les mécanismes de déformation des monocristaux de zinc. / The identification of certain behavior laws of remains delicate especially for nonlinear behaviors. Several identification strategies have been developed based mainly on the experiment-calculation dialogue. The experimental part consists in measuring one or more physical quantities which characterize the behavior of the tested material and the theoretical part makes it possible to calculate these same quantities. To go back to the material properties sought, an agreement between the quantities measured and calculated must be aimed at. In this thesis, we are interested in zinc, which is a polycrystalline material. It is modeled through crystal plasticity laws that predict its mechanical behavior under complex loadings and the evolution of its microstructure. The crystal plasticity laws implemented in the commercial software Abaqus by Huang and Marin that we adapted to the case of zinc were used. The step of identifying these non-linear laws is based on an inverse approach that takes into account the heterogeneities of the deformations at the grain scale. For this, we used the nanoindentation test which is a heterogeneous test to extract more information than the homogeneous tests. For this purpose, and to characterize the crystal plasticity of zinc, we used the load - penetration depth curves and the displacement profiles measured on the residual imprint on grains of different crystallographic orientations measured by EBSD. The confrontation of the experimental and numerical results show the good agreement between the experimental data and the identified models. The results obtained made it possible to characterize the deformation mechanisms of zinc single crystals.

L’apparition des fissures macroscopiques peut manifestement se déclencher à plusieurs étapes du processus de la vie d’une cathode dans une cuve d’électrolyse. Diverses causes peuvent engendrer un tel surgissement. Au préchauffage électrique, par exemple, la répartition spatiale irrégulière des gradients thermiques peut conduire à la naissance d’une ou plusieurs fissures. L’objectif principal de ce travail se résume à la modélisation et la simulation de la propagation de fissures dans un contexte multiphysique, thermo-électro-mécanique. L’approche se base sur une philosophie de développement numérique monolithique, du couplage fortement non linéaire avec la prise en compte des échanges thermoélectriques à travers les lèvres de fissures, par la méthode des éléments finis étendus (XFEM). Ce projet a été dressé sous forme de trois sous-projets. En premier lieu, un nouveau modèle de la propagation de fissures en contexte thermomécanique dans un domaine fissuré a été développé. Une stratégie d’enrichissement géométrique par XFEM a été élaborée avec succès pour modéliser la propagation de fissures ainsi que pour atteindre le taux optimal de convergence espéré sur les réponses physiques et sur le calcul des facteurs d’intensité de contraintes. Trois benchmarks ont été examinés et validés à partir des résultats existants dans la littérature. Ensuite, un nouvel exemple de propagation de fissures et défauts multiples a été conçu. La comparaison des résultats obtenus montre la robustesse ainsi que la précision du développement numérique. En second lieu, une approche originale de la modélisation du couplage fort thermoélectrique, en présence d’une discontinuité matérielle forte statique et en tenant compte des échanges thermique et électrique à travers ses lèvres, a été développée. Tout d’abord, un modèle a été élaboré et validé dans le cas d’une fissure capacitive, où les propriétés du diélectrique, entre les lèvres de la fissure, jouent un rôle déterminant pour la quantification de son intensité d’échange. De même, un modèle a été confirmé pour le choix d’un milieu gazeux remplissant la région entre les deux lèvres, assurant l’aspect résistif d’échange thermique. Le modèle thermoélectrique a été mis en œuvre par XFEM, en implémentant la version améliorée, et en tenant compte du traitement des parasites numériques dus aux éléments de transition. Trois cas d’études ont été proposés, le premier a été appliqué pour valider la mise en œuvre numérique du modèle proposé via un patch test développé. Une bonne convergence et précision des résultats numériques ont été observées. Le second a porté sur le cas d’une fissure courbée, moins fréquente dans la littérature, qui peut être rencontré dans de nombreuses applications, et représente un défi numérique notable. Le troisième est une étude élaborée sur la cathode. L’analyse de l’impact de la présence des fissures sur l’efficacité énergétique a été soulevée aussi pour le cas de fissures multiples. En troisième lieu, une originalité numérique de la propagation de fissures en couplage thermo-électro-mécaniques (TEM) avec les échanges thermique et électrique entre les lèvres de la fissure a été exposée. La fissure n’est pas statique, cette fois-ci. L’aire générée par le déplacement des lèvres est prise en considération dans les expressions physiques des coefficients d’échanges thermique et électrique. Cela est transcrit par la quantification du saut de déplacement à travers les lèvres. Cet aspect rend le problème mécanique non linéaire, aussi, par le biais des échanges, et par conséquent le système global TEM est fortement non linéaire. Une approche originale du point de vue de la technicité d’intégration a été implémentée. Elle est fondée sur une technique de sous-triangulation barycentrique par une ’toile d’araignée’ pour les éléments de surface coupés par la fissure. Une autre technique basée sur l’intégration par ’bras’ de fissure réservée pour les éléments d’échanges thermique et électrique. Deux cas d’études essentiels ont été envisagés. Le premier est un exemple de propagation en mode mixte d’une fissure inclinée en TEM avec les échanges thermique et électrique en fonction du saut de déplacement. Le deuxième a été réservé pour un cas d’étude industrielle d’une cuve en opération, en contexte TEM. L’impact de la présence de la fissure sur les différentes réponses physiques est analysé. Comme constatation, en expansion thermoélectrique du bloc cathodique, la fissure a plus tendance à rejoindre la surface supérieure pour mettre, probablement, fin à la vie de la cathode. / The outbreak of visible cracks can be triggered at several stages in the life of a cathode block in an electrolysis cell. Various matters can prompt such an upsurge. Under electrical pre-heating, for example, the random spatial distribution of thermal gradients can lead to the rise of one or more cracks. The main objective of this work is to model and simulate the crack propagation in a multi-physical, thermo-electro-mechanical, context. The approach is based on a monolithic numerical development philosophy of a strongly nonlinear coupling, with the consideration of thermoelectric exchanges through the crack lips using the extended finite element method (XFEM). This project was sketched essentially on three sub-projects. Firstly, a new model of crack propagation in a thermomechanical fashion over a cracked domain has been developed. A geometrical enrichment strategy by XFEM has been successfully performed to model crack growth as well as to achieve the expected optimal convergence rate in physical responses and the computation of stress intensity factors. Three benchmarks were examined and validated from existing results in the literature. Then, anew example of the propagation of multiple cracks and multiple defects was designed. The comparison of the results obtained shows a good agreement with the reported works as well as the robustness and the accuracy of the numerical development. Secondly, an original approach to the modeling of full thermoelectric coupling, in the presence of a strong static material discontinuity and taking into account the thermal and electrical exchanges through its lips, has been elaborated. First, a model has been established and validated in the case of a capacitive crack, where the properties of the dielectric, between the inter-crack-lips, play a determining role for the ascertainment of its exchange intensity. Similarly, a model has been settled for the choice of a gaseous medium filling the gap between the two rims, ensuring the resistive heat exchange aspect. The thermoelectric model has been implemented via XFEM by performing the enhanced version and taking into account the treatment of numerical noise due to the blending elements. Three case studies were intended, the first one was employed to validate the numerical implementation of the stated model via a developed patch test. Good convergence and accuracy of numerical outcomes have been noted. The second focused on the case of a curved crack, less considered in the literature, which can be encountered in many applications, and represents a significant numerical challenge. The third is an elaborate study on the cathode. The analysis of the impact of cracks on energy efficiency was also raised for the case of multiple cracks. Thirdly, numerical originality of crack propagation in the context of the full thermo-electro-mechanical (TEM) coupling combined with thermal and electrical exchanges between the crack lips has been displayed. The area formed by the displacement of the two crack bank sis carried out in the physical expressions of the heat and electrical coefficients. The quantification of this gap is rendered by the resulting displacement jump across the lips. This aspect makes the mechanical problem nonlinear as well through exchanges, and therefore the global TEM system is strongly nonlinear. An original integration strategy, from a technical point of view, has been realized. It is based on a technique of barycentric sub-triangulation through a ’spider-web’ for the surface elements cut by the crack. Another procedure based on integration by crack ’arms’ reserved for heat and electrical exchanges elements. Two case studies have been examined. The first one is an example of the mixed-mode growth of a sloped crack in TEM context with thermal and electrical exchanges as a function of the displacement jump. The second was held for an industrial case of an electrolysis cell in operation and TEM circumstances. The impact of the presence of crack on several physical responses is analyzed. As an outcome, due to the thermoelectric expansion of the cathode block, the crack is expected to join the upper surface and lead to the end of the cathode’s life.

La géométrie et la flexibilité d'un canyon sont les paramètres qui affectent grandement la valeur des périodes naturelles dans les barrages en terre. Le canyon entourant des barrages peut être considéré comme un domaine illimité. Pour prendre en compte ces deux effets, le canyon a été modélisé par SBFEM et le barrage en terre, à géométrie limitée, par FEM. La technique hybride SBFEM-FEM pour l'analyse tridimensionnelle dynamique de l'interaction sol-barrage a été validée avec les résultats disponibles dans la littérature. Comme la matrice de rigidité dynamique du domaine non borné est complexe et dépendante de la fréquence, la méthode classique de superposition de modes n'est pas simple pour le système d'interaction sol-structure. Ainsi, pour obtenir la fréquence propre fondamentale, le barrage a été excité en direction amont-aval. Les périodes naturelles du barrage de terre pour des canyons de formes géométriques et de coefficient de impédance différents ont été obtenues. Ils se sont avérés avoir des effets significatifs sur la période naturelle. Les résultats ont été comparés aux données enregistrées réelles. Il a été constaté que les graphiques proposés dans cette étude peuvent être utilisés par des concepteurs de barrages pour l'estimation des périodes naturelles des barrages en terre dans des canyons de formes et de propriétés matérielles différentes. Plusieurs fonctions d'amplification correspondant à différentes conditions de canyon ont été obtenues en appliquant un déplacement uniforme à la limite du canyon. Une étude approfondie a été réalisée pour examiner les effets de la géométrie et de la flexibilité du canyon sur la réponse en régime permanent du barrage. Ces deux effets ont influencé de manière importante les fonctions d'amplification. Alors que la flexibilité du canyon affecte de manière significative la valeur de la fonction d'amplification maximale, cette valeur ne change pas pour les barrages en terre dans lesquels les canyons ont des formes différentes et la même longueur. De plus, la réponse latérale du barrage de terre dans le domaine temporel a été calculée pour analyser les effets susmentionnés lors d'un tremblement de terre réel. Les fonctions d'amplification proposées ont été utilisées pour comparer les spectres de réponse enregistrés du barrage d'El Infiernillo lors des tremblements de terre de 1966 avec la fonction d'amplification calculée. Un accord raisonnable a été observé entre eux. La méthode linéaire équivalente (EQL) a été implémentée dans le FEM. La technique FEMSBFEM a été étendue pour prendre en compte l'effet du comportement non linéaire des barrages en terre. Il a été observé que le comportement non linéaire affecte grandement la fréquence naturelle, la fonction d'amplification et l'accélération de crête maximale du barrage de terre situé dans les canyons. Les effets de la géométrie et de la flexibilité du canyon sur le comportement non linéaire ont été examinés, et on a vu qu'en augmentant la flexibilité du canyon, l'effet de la non-linéarité était diminué. Le barrage d'El Infiernillo a été modélisé par FEM-SBFEM non linéaire 3D, et une comparaison de la fonction d'amplification de crête obtenue par la méthode proposée avec les données enregistrées montre la précision du FEM-SBFEM non linéaire. / The canyon surrounding a dam can be assumed as an unbounded domain, and the geometry and flexibility of a canyon are parameters that greatly affect the values of natural periods in earth dams. In this thesis, in order to take into account these two effects, canyons are modeled by SBFEM, and earth dams, which have limited geometries, are modeled by FEM. The hybrid FEM-SBFEM technique used for the dynamic three-dimensional analysis of soil-earth dam interactions is validated with results available in the literature. Because the dynamic-stiffness matrix of the unbounded domain is complex and frequency-dependent, the classical mode-superposition method is not straightforward for a soil-structure interaction system, and thus, to obtain their fundamental natural frequencies, the modeled dams were excited in the upstream-downstream direction. The natural periods of earth dams in canyons with different geometries shapes and impedance ratios are obtained, and are found to have significant effects on the dams' natural periods. The results are compared with actual recorded data, and it is found that the graphs put forward in this study may be used by practical engineers for the estimation of natural periods of earth dams in canyons with different shapes and material properties. Several amplification functions corresponding to different canyon conditions are obtained by applying a uniform displacement at the canyons' boundaries. A comprehensive study is performed to examine the effects of canyon geometry and flexibility on the steady-state responses of the dams, and it is found that these two effects significantly influence the amplification functions. While the flexibility of the canyon does affect the maximum amplification function value, this value does not change for earth dams in canyons that have different shapes but the same length. In addition, the lateral responses of earth dams in the time domain are computed in order to analyze the aforementioned effects under an actual earthquake. The proposed amplification functions are used to compare the recorded response spectra of the El Infiernillo dam under the two 1966 earthquakes with the calculated amplification function, and a reasonable agreement is observed between them. The equivalent linear method (EQL) is implemented into the FEM, and the FEM-SBFEM technique is extended in order to take into consideration the effect of earth dams' nonlinear behavior. It is observed that such nonlinear behavior greatly affects the natural frequency, the amplification function, and peak crest acceleration of earth dams located in canyons. The effects of canyon geometry and flexibility on the nonlinear behavior are examined, and it is found that by increasing canyon flexibility, the effect of nonlinearity is decreased. The El Infiernillo dam is modeled by the 3D nonlinear FEM-SBFEM, and comparison of the crest amplification function obtained by the proposed method with the recorded data shows the accuracy of the nonlinear FEM-SBFEM.

Ce projet est le résultat d'une collaboration entre le laboratoire de vision et systèmes numériques (LVSN) de l'Université Laval et le département de Radio-oncologie du CHUQ. Il consistait à développer une méthode de simulation de la résorption d'un oedème de la prostate. La méthode développée est basée sur l'algorithme des masses-tenseurs élastique linéaire et sur une discrétisation en éléments finis. Elle consiste à étirer artificiellement le contour de la prostate pour simuler l'oedème, puis à relâcher et observer le comportement. Cette méthode a l'avantage de donner un modèle tridimensionnel de la résorption de l'oedème contrairement aux modèles existants qui sont des modèles de modélisation de la décroissance du volume seulement, comme le modèle exponentiel par exemple. Les résultats obtenus avec cette méthode sont assez proches du comportement réel. Ces résultats seront éventuellement exploités dans la boucle de calcul de dosimétrie pour le traitement du cancer de la prostate par curiethérapie avec implants permanents. Le modèle visco-élastique linéaire servira également à la réalisation d'un simulateur de traitement en réalité virtuelle.

Le comportement de l'acier est significativement différent à des températures élevées par rapport à la température ambiante. Avec l'augmentation de la température, la relation contrainte-déformation devient non linéaire et la résistance et la rigidité de l'acier diminuent. Les normes de conception au feu de l'acier existantes ne parviennent pas à prédire avec précision la résistance à des températures élevées en suivant les mêmes approches de conception à température ambiante avec de petites modifications. Cette recherche vise à proposer une approche de conception innovante des profilés creux en acier en situation d'incendie au moyen du Overall Interaction Concept (O.I.C.). Le Overall Interaction Concept (O.I.C.) est basé sur l'interaction résistance-instabilité en utilisant un rapport d'élancement généralisé et s'est avéré plus précis et cohérent que les normes actuelles. Un modèle d'éléments finis a été développé pour étudier le comportement des sections carrées et rectangulaires en acier à haute température. La précision du modèle a été confirmée par la validation avec des résultats expérimentaux. Le modèle validé a été utilisé pour effectuer des études paramétriques incluant différentes dimensions géométriques, températures et cas de charge. Les résultats des études paramétriques ont ensuite été utilisés pour identifier le paramètre de conception principal et proposer O.I.C. équations en situation d'incendie. Pour présenter un meilleur niveau de précision et de cohérence de l'O.I.C. propositions, leur performance a été comparée aux approches actuelles de conception anti-incendie, y compris les normes européennes, américaines et canadiennes. / The behavior of steel is significantly different at elevated temperatures compared to the ambient temperature. With increasing temperature, the stress-strain relationship gets nonlinear and the strength and stiffness of steel reduce. Existing steel fire design standards fail to accurately predict the resistance at high temperatures by following the same design approaches at room temperature with small modifications. This research is aimed to propose an innovative design approach for steel hollow sections in fire situations by means of the Overall Interaction Concept (O.I.C.). The Overall Interaction Concept (O.I.C.) is based on the resistance-instability interaction using a generalized slenderness ratio and has been proved to be more accurate and consistent than the current standards. A finite element (F.E.) model was developed to study the behavior of steel square and rectangular sections at high temperatures. The accuracy of the model was confirmed through validation with experimental results. The validated model was used to perform parametric studies including different geometrical dimensions, temperatures, and load cases. The results of the parametric studies were then used to identify the leading design parameter and propose O.I.C. equations in fire situations. To exhibit a better level of accuracy and consistency of the O.I.C. proposals, their performance was compared to the current fire design approaches including European, American, and Canadian Standards.

L'objectif de ce mémoire est de présenter une méthode d'adaptation de maillage pour la méthode des éléments finis plus efficace que les méthodes dont on dispose actuellement. C'est d'autant plus important en dimension 3 puisque les problèmes rencontrés sont naturellement plus complexes et demandent de grandes ressources computationnelles. On présentera d'abord quelques notions essentielles et on étudiera plus particulièrement les méthodes d'adaptation existantes avant d'approfondir notre méthode. Celle-ci se base sur un estimateur d'erreur hiérarchique et montre des différences marquées avec les méthodes classiques. On fera donc la construction de l'estimateur en 2D avant de généraliser l'idée en 3D et d'expliquer le cadre dans lequel cet estimateur guide l'adaptation. Après avoir abordé certains aspects spécifiques de notre méthode, on fournira enfin des résultats numériques permettant entre autres de la comparer avec l'adaptation avec métrique.