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.

Dans ce travail on s'intéresse à la résolution des équations de Saint -Venant par la méthode des éléments finis. Deux tests numériques sont présentés en utilisant deux paires d 'éléments finis : P₁NC - P₁ et RT0 - P₀. Lorsque le systèrne d'équations de Saint-Venant admet une solution analytique, cette solution est comparée avec celle obtenue par la méthode numérique pour montrer la robustesse de ces éléments finis.

Malgré de nombreux travaux sur le sujet, la résolution des problèmes de contact constitue encore un défi pour le numéricien. Il existe deux types de problèmes de contact : soit le contact unilatéral et le contact frottant. Dans cette thèse, nous allons uniquement considérer le contact unilatéral. De plus, nous allons nous limiter au cas de l'élasticité linéaire. Malgré ces hypothèses, les principales difficultés du problème en grandes déformations y sont présentes. La difficulté majeure provient de la non différentiabilité engendrée par la contrainte d'inéquation du contact. Plusieurs auteurs utilisent des méthodes de régularisation afin d'obtenir un problème différentiable soluble par l'algorithme de Newton. Nous croyons que cette stratégie pose problème et qu'il est préférable de traiter le problème directement à partir de l'inéquation variationnelle. Dans la thèse, nous allons proposer plusieurs stratégies de résolution par éléments finis du problème de contact unilatéral. Tous les algorithmes sont basés sur les méthodes puissantes et efficaces de l'optimisation convexe. En premier lieu, on propose un algorithme de type gradient conjugué avec projection sur le cône positif. Cet algorithme baptisé GCP, s'est révélé fort efficace dans le cas du contact entre un solide déformable et une fondation rigide. On propose aussi un algorithme ALG3 basé sur une formulation à trois champs du problème de contact résolu par la méthode du Lagrangien Augmenté. Finalement on propose un algorithme ALGCP qui combine à la fois l'approche du Lagrangien Augmenté ainsi que la méthode du Gradient Conjugué Projeté (GCP). Les résultats numériques montreront la supériorité de ALG3 pour les problèmes de contact unilatéral à plusieurs corps.

La modélisation des écoulements en milieux naturels tels que les rivières, les lacs et les océans, fait souvent intervenir un système d'équations aux dérivées partielles dit de Saint-Venant. La plupart des méthodes numériques utilisées pour résoudre les équations de Saint-Venant générent des modes purement numériques en approximant les ondes de type inertie-gravité. Les modes parasites les plus dangereux sont les modes stationnaires. Ils conduisent généralement à des solutions erronées. Ce travail propose une étude des modes stationnaires résultants de la discrétisation des équations de Saint-Venant par des méthodes aux différences finies et d'éléments finis. Nous privilégions une approche de type algèbre linéaire au lieu de celle de type Fourier largement utilisée. Nous introduisons une nouvelle nomenclature des modes parasites qui dépasse largement celle utilisée jusqu'à présent et généralement restreinte aux seuls modes parasites pression. Enfin l'approche de type algèbre linéaire utilisée ici nous permet de tirer quelques conclusions préliminaires quant à la manifestation des modes parasites sur des maillages non structurés.

Une des applications des équations de réaction-diffusion est leur utilisation pour la génération de motifs, analogues en particulier à ceux sur le pelage de certains animaux. L'utilisation à cette fin de ce type d'équations remontent aux travaux d'Alan Turing dans les années 50. Dans ce mémoire nous considérons un système simple d'équations de réaction-diffusion pour étudier la formation de motifs. Nous approchons ces solutions en utilisant une discrétisation spatiale par éléments finis. Nous montrons comment choisir les paramètres pour atteindre un motif choisi sur un domaine rectangulaire fixe. Après avoir validé les méthodes numériques utilisées dans nos simulations, nous nous intéressons à l'évolution de motifs sur des domaines qui croissent dans le temps. En particulier, nous examinons l'influence de la condition initiale et de la vitesse de croissance du domaine sur les motifs apparaissant au terme de la croissance.

La rentabilité des cellules à base de silicium est un point essentiel pour le marché photovoltaïque et cela passe notamment par l'amélioration du rendement électrique, la baisse des coûts de production ainsi que le renforcement de la fiabilité/durabilité des wafers. Des procédés innovants émergent, qui permettent d'obtenir des wafers ultra minces avec moins de perte de matière première. Cependant il est nécessaire de mettre en place des méthodes de caractérisation afin d’analyser la rigidité et la tenue mécanique de ces matériaux. Dans ce travail, des essais de flexion ont été effectués pour caractériser à la fois la rigidité et la rupture. Afin d’étudier la rupture fragile, une caméra rapide a été utilisée, des analyses fractographiques ont été menées. La diffraction d'électrons rétrodiffusés et la diffraction par rayon X de Laue ont été utilisées afin d'explorer le lien entre les orientations cristallographiques et les comportements observés. Conjointement, des simulations numériques EF ont été mise en place. Grâce à ce couplage expériences-simulations numériques, une caractérisation fiable de la rigidité des wafers a été effectuée. Une stratégie d'identification de l'origine de la rupture est également proposée. L'étude de la rupture du silicium monocristallin a mis en évidence la stabilité du clivage (110), la grande vitesse d'amorçage de la fissure, la dépendance de la forme du front de fissure à la vitesse de propagation ainsi que l'apparition de "Front Waves" pour les fissures à très grande vitesse. L'étude de la rupture des wafers multi-cristallins démontre une fissuration intra-granulaire. Des éprouvettes jumelles ont permis d’étudier la répétabilité du chemin de fissuration : une attention particulière a été portée à la nature des plans de clivage ainsi que l'effet des joints de grains. Enfin, une modélisation par la méthode des éléments finis étendus est proposée. Elle permet de reproduire le chemin de fissuration expérimentalement observé. / The profitability of silicon solar cells is a critical point for the PV market and it requires improved electrical performance, lower wafer production costs and enhancing reliability and durability of the cells. Innovative processes are emerging that provide thinner wafers with less raw material loss. But the induced crystallinity and distribution of defects compared to the classical wafers are unclear. It is therefore necessary to develop methods of microstructural and mechanical characterization to assess the rigidity and mechanical strength of these materials. In this work, 4-point bending tests were performed under quasi-static loading. This allowed to conduct both the stiffness estimation and the rupture study. A high speed camera was set up in order to track the fracture process thanks to a 45° tilted mirror. Fractographic analysis were performed using confocal optical microscope, scanning electron microscope and atomic force microscope. Electron Back-Scatter Diffraction and Laue X-Ray diffraction were used to explore the relationship between the microstructural grains orientations/textures of our material and the observed mechanical behavior. Jointly, finite element modeling and simulations were carried out to provide auxiliary characterization tools and help to understand the involved fracture mechanism. Thanks to the experiment-simulation coupled method, we have assessed accurately the rigidity of silicon wafers stemming from different manufacturing processes. A fracture origin identification strategy has been proposed combining high speed imaging and post-mortem fractography. Fracture investigations on silicon single crystals have highlighted the deflection free (110) cleavage path, the high initial crack velocity, the velocity dependent crack front shape and the onset of front waves in high velocity crack propagation. The investigations on the fracture of multi-crystalline wafers demonstrate a systematic transgranular cracking. Furthermore, thanks to twin multi-crystalline silicon plates, we have addressed the crack path reproducibility. A special attention has been paid to the nature of the cleavage planes and the grain boundaries barrier effect. Finally, based on these observations, an extended finite element model (XFEM) has been carried out which fairly reproduces the experimental crack path.

Dans ce projet, le matériau est l'alliage de magnésium AZ31-O en tôle. L'épaisseur de tôles est de 1,2 mm. Les essais de traction à chaud sont réalisés afin d'étudier la ductilité de l'alliage de magnésium AZ31-O, la température et l'influence la vitesse de déformation sont incluses dans tous les tests. Le résultat d'analyse montre que la ductilité est renforcée avec une température croissante et une vitesse de déformation décroissante, le phénomène d'adoucissement est évident à la température élevée. La propriété anisotrope n'est pas considérée dans ce projet. Les essais Nakazima à chaud avec le poinçon d'hémisphère sont réalisés pour étudier la formabilité de l'alliage de magnésium AZ31-O. Enfin, la FLD (Forming Limit Diagram) est identifiée et les comparaisons montrent que la formabilité est préférable à une température plus élevée. En outre, les prédictions des limites de formage sont effectuées dans le modèle M-K. La comparaison montre clairement avec la prédiction théorique ne convient pas avec l'expérience. Les simulations des éléments finis sont effectuées pour un emboutissage par poinçon hémisphérique et un emboutissage en croix. Tout d'abord, les simulations d'emboutissage de poinçon hémisphérique sont réalisés sur FORGE® et sur ABAQUS®. Les résultats des simulations de FORGE et de ABAQUS sont comparés afin d'étudier la différence de divers codes de simulation des éléments finis. Deuxièmement, le comportement de d'endommagent est étudié dans FORGE par modèle d'endommagement Lemaitre. Enfin, la simulation d'emboutissage en croix qui est un benchmark de la conférence 2011 NUMISHEET est réalisée avec FORGE. La charge de poinçon, l'épaisseur et la distribution de température sont obtenues et comparées pour chaque simulation. En outre, ces résultats de la simulation de benchmark (FORGE) sont également comparés à d'autres logiciels de simulation en conférence. Les résultats des analyses détaillées sont présentés dans cette thèse.

La prise en compte des pertes fer dans les machines électriques est une étape cruciale lors du processus de conception. L'objectif de ce travail a donc été de valider l'implémentation d'un modèle de pertes fer tout en effectuant une modélisation précise des machines électriques étudiées. Le choix s'est porté sur l'utilisation de modèles de calcul des pertes en post-traitement d'un code par éléments finis. Ainsi, différents modèles analytiques ont été comparés en régime sinusoïdal et non sinusoïdal (sans cycles mineurs) en termes de précision, afin de choisir le plus adapté pour une implantation dans code_Carmel. Par ailleurs, nous avons montré que la prise en compte des cycles mineurs peut être améliorée en remplaçant le terme classique de Steinmetz par un modèle d'hystérésis ou par un modèle analytique adapté. L'implémentation et la précision du modèle analytique ont été évaluées d'abord sur des systèmes simples puis sur plusieurs machines tournantes. La comparaison des résultats simulés avec les mesures issues des essais à vide montre une bonne concordance et démontre la capacité du modèle à estimer les pertes fer avec une bonne précision pour une large gamme de machines électriques et pour les deux types d'alimentations utilisés (sinusoïdale et MLI).

Cette thèse porte sur la modélisation des ondes électromagnétiques en milieu complexe et à haute fréquence. Ceci reste un défi, de part les besoins en temps de calcul et en place mémoire que cela suscite. L'idée de la méthode Finite Element Tearing and Interconnect (FETI) est de diviser le domaine en sous-domaines distincts dans lesquels des problèmes électromagnétiques peuvent être résolus de manière indépendante. Une modification de la méthode classique, appelée FETI-DPEM2-full, est ici proposée afin d'améliorer le processus numérique. Les champs diffractés par différents objets ont ainsi été simulés et comparés avec succès à des champs expérimentaux. A partir de ces champs, les techniques d'imagerie micro-onde permettent de déterminer les paramètres électromagnétiques des objets illuminés. Ce problème inverse peut être réécrit sous la forme d'un problème d'optimisation faisant appel à la résolution de deux problèmes directs à chaque itération. Une méthode de type quasi-Newton a donc été couplée efficacement avec la méthode FETI-DPEM2-full. Des reconstructions obtenues à partir de champs mesurés sont présentées et valident la chaîne méthodologique proposée. / This PhD work is devoted to the electromagnetic modeling of threedimensional large-scale wave propagation problems, which is very challenging in terms of memory and computation time. The main idea of the Finite Element Tearing and Interconnect (FETI) method is to divide the domain into non-overlapping subdomains where each electromagnetic problem can be independently solved. A modification of the classical FETI method, called the FETI-DPEM2-full method, is proposed herein to improve the numerical process. The fields scattered by various structures have thus been computed and successfully compared with measured fields. From these experiments, quantitative microwave imaging algorithms attempt to estimate the physical parameters of the illuminated target. This inverse problem is recasted into a minimization problem where two forward problems are required at each iteration step. A quasi-Newton algorithm has thus been efficiently coupled with the FETI-DPEM2-full method. Reconstructions of various targets from measured scattered fields have been successfully performed, which validates the effectiveness of the proposed methodology.

Les Surfaces à Haute Impédance (SHI) ont été largement étudiées pour améliorer toutes sortes de performances des antennes, comme le gain, le facteur de qualité, les formes et dimensions. L'objectif de cette thèse est de modéliser les structures de SHI et de caractériser leurs performances en vue de futures applications aux antennes.Après une brève introduction aux structures SHI et une étude de quelques modèles analytiques fréquemment traités dans la littérature, deux nouvelles méthodes numériques sont proposées pour calculer l'impédance de surface de structures SHI. Ces deux méthodes (dites « méthode du flux de Poynting » et « méthode <E>/<H> ») sont validées sur des structures symétriques, puis mises en service sur des structures de SHI asymétriques. Elles sont également validées par comparaison de résultats analytiques, numériques et expérimentaux.Nous présentons ensuite un modèle équivalent basé sur l'idée de remplacer les structures hétérogènes de SHI par une surface homogène, caractérisée par son impédance surfacique. Ce modèle nous permet d'avoir une prédiction avec un temps de calcul et une occupation de mémoire PC largement réduits. / High impedance surfaces (HIS) have proved good candidates for antenna miniaturization or antennas performance improvement. Within a certain frequency band, they can enhance the gain of an antenna while simultaneously suppressing the unwanted surface waves.In this thesis, the focus of our work is on numerical modeling of these structures by using the finite element method (FEM) based on edge elements.One of our contributions is that we propose two new numerical methods (the Poynting flux method and <E>/<H> method) to calculate the surface impedance not only for HIS structures with symmetric geometries, but also for those with asymmetric geometries. These two numerical methods have been validated by comparing analytical, numerical and experimental results.Another significant contribution of the thesis is that we introduce an equivalent model, based on the idea of replacing the heterogeneous HIS structures by a homogeneous surface, characterized by its surface impedance. Compared with the normal model, this equivalent model can save computing time and memory space.