La liquéfaction des sols qui est déclenchée par des mouvements sismiques forts peut modifier la réponse d’un site. Ceci occasionne des dégâts importants dans les structures comme a été mis en évidence lors des tremblements de terre récents tels que celui de Christchurch, Nouvelle-Zélande et du Tohoku, Japon. L’évaluation du risque sismique des structures nécessite une modélisation robuste du comportement non linéaire de sols et de la prise en compte de l’interaction sol-structure (ISS). En général, le risque sismique est décrit comme la convolution entre l’aléa et la vulnérabilité du système. Cette thèse se pose comme une contribution à l’étude, via une modélisation numérique, de l’apparition de la liquéfaction et à l’utilisation des méthodes pour réduire les dommages induits.A cet effet, la méthode des éléments finis(FEM) dans le domaine temporel est utilisée comme outil numérique. Le modèle principal est composé d’un bâtiment fondé sur un sable liquéfiable. Comme la première étape de l’analyse du risque sismique, la première partie de cette thèse est consacrée à la caractérisation du comportement du sol et à sa modélisation.Une attention particulière est donnée à la sensibilité du modèle à des paramètres numériques. En suite, le modèle est validé pour le cas d’une propagation des ondes 1D avec les mesures issus du benchmark international PRENOLIN sur un site japonais. D’après la comparaison, le modèle arrive à prédire les enregistrements dans un test en aveugle.La deuxième partie, concerne la prise en compte dans la modélisation numérique du couplage de la surpression interstitielle (Δpw)et de la déformation du sol. Les effets favorables ou défavorables de ce type de modélisation ont été évalués sur le mouvement en surface du sol lors de la propagation des ondes et aussi sur le tassement et la performance sismique de deux structures.Cette partie contient des éléments d’un article publié dans Acta Geotechnica (Montoya-Noguera and Lopez-Caballero, 2016). Il a été trouvé que l’applicabilité du modèle dépend à la fois du niveau de liquéfaction et des effets d’ISS.Dans la dernière partie, une méthode est proposée pour modéliser la variabilité spatiale ajoutée au dépôt de sol dû à l’utilisation des techniques pour diminuer le degré de liquéfaction. Cette variabilité ajoutée peut différer considérablement de la variabilité inhérente ou naturelle. Dans cette thèse, elle sera modélisée par un champ aléatoire binaire.Pour évaluer l’efficience du mélange, la performance du système a été étudiée pour différents niveaux d’efficacité, c’est-à-dire,différentes fractions spatiales en allant de non traitées jusqu’à entièrement traitées. Tout d’abord le modèle binaire a été testé sur un cas simple, tel que la capacité portante d’une fondation superficielle sur un sol cohérent.Après, il a été utilisé dans le modèle de la structure sur le sol liquéfiable. Ce dernier cas,en partie, a été publié dans la revue GeoRisk (Montoya-Noguera and Lopez-Caballero,2015). En raison de l’interaction entre les deux types de sols du mélange, une importante variabilité est mise en évidence dans la réponse de la structure. En outre, des théories classiques et avancées d’homogénéisation ont été utilisées pour prédire la relation entre l’efficience moyenne et l’efficacité. En raison du comportement non linéaire du sol, les théories traditionnelles ne parviennent pas à prédire la réponse alors que certaines théories avancées qui comprennent la théorie de la percolation peuvent fournir une bonne estimation. En ce qui concerne l’effet de la variabilité spatiale ajoutée sur la diminution du tassement de la structure, différents séismes ont été testés et la réponse globale semble dépendre de leur rapport de PHV et PHA. / Strong ground motions can trigger soil liquefaction that will alter the propagating signal and induce ground failure. Important damage in structures and lifelines has been evidenced after recent earthquakes such as Christchurch, New Zealand and Tohoku, Japanin 2011. Accurate prediction of the structures’ seismic risk requires a careful modeling of the nonlinear behavior of soil-structure interaction (SSI) systems. In general, seismic risk analysisis described as the convolution between the natural hazard and the vulnerability of the system. This thesis arises as a contribution to the numerical modeling of liquefaction evaluation and mitigation.For this purpose, the finite element method (FEM) in time domain is used as numerical tool. The main numerical model consists of are inforced concrete building with a shallow rigid foundation standing on saturated cohesionless soil. As the initial step on the seismic risk analysis, the first part of the thesis is consecrated to the characterization of the soil behavior and its constitutive modeling. Later on, some results of the model’s validation witha real site for the 1D wave propagation in dry conditions are presented. These are issued from the participation in the international benchmark PRENOLIN and concern the PARI site Sendaiin Japan. Even though very few laboratory and in-situ data were available, the model responses well with the recordings for the blind prediction. The second part, concerns the numerical modeling of coupling excess pore pressure (Δpw) and soil deformation. The effects were evaluated on the ground motion and on the structure’s settlement and performance. This part contains material from an article published in Acta Geotechnica (Montoya-Noguera andLopez-Caballero, 2015). The applicability of the models was found to depend on both the liquefaction level and the SSI effects.In the last part, an innovative method is proposed to model spatial variability added to the deposit due to soil improvement techniques used to strengthen soft soils and mitigate liquefaction. Innovative treatment processes such as bentonite permeations and biogrouting,among others have recently emerged.However, there remains some uncertainties concerning the degree of spatial variability introduced in the design and its effect of the system’s performance.This added variability can differ significantly from the inherent or natural variability thus, in this thesis, it is modeled by coupling FEM with a binary random field. The efficiency in improving the soil behavior related to the effectiveness of the method measured by the amount of soil changed was analyzed. Two cases were studied: the bearing capacity of a shallow foundation under cohesive soil and the liquefaction-induced settlement of a structure under cohesionless loose soil. The latter, in part, contains material published in GeoRisk journal (Montoya-Noguera and Lopez-Caballero, 2015). Due to the interaction between the two soils, an important variability is evidenced in the response. Additionally, traditional and advanced homogenization theories were used to predict the relation between the average efficiency and effectiveness. Because of the nonlinear soil behavior, the traditional theories fail to predict the response while some advanced theories which include the percolation theory may provide a good estimate. Concerning the effect of added spatial variability on soil liquefaction, different input motions were tested and the response of the whole was found to depend on the ratio of PHV and PHA of the input motion.

Dans ces travaux de thèse on s’intéresse au développement d’un outil numérique pour résoudre le problème de conduction instationnaire avec changement de phase dans un milieu composite constitué d’une mousse de graphite infiltrée par un matériau à changement de phase tel que le sel, dans le contexte du stockage de l’énergie thermique solaire.Au chapitre 1, on commence par présenter le modèle sur lequel on va travailler. Il estséparé en trois sous-parties : un problème de conduction de chaleur dans la mousse, un problème de changement de phase dans les pores remplis de sel et une condition de résistance thermique de contact entre les deux matériaux qui est traduite par une discontinuité du champ de température.Au chapitre 2, on étudie le problème stationnaire de conduction thermique dans un milieu composite avec résistance de contact. Ceci permet de se focaliser sur la plus grande difficulté présente dans le problème qui est le traitement de la condition de saut à l’interface.Deux méthodes d’éléments finis sont proposées pour résoudre ce problème : une méthode basée sur les éléments finis Lagrange P1 et une méthode hybride-duale utilisant les éléments finis Raviart-Thomas d’ordre 0 et P0. L’analyse numérique des deux méthodes est effectuée et les résultats de tests numériques attestent des efficacités des deux méthodes [10]. Les matériaux à changement de phase qu’on étudie dans le cadre de cette thèse sont des matériaux pures, par conséquent le changement de phase s’effectue en une valeur de température fixe qui est la température de fusion. Ceci est modélisé par un saut dans la fonction fraction liquide et par conséquent dans la fonction enthalpie du matériau. Cette discontinuité représente une difficulté numérique supplémentaire qu’on propose de surmonter en introduisant un intervalle de régularisation autour de la température de fusion.Cette procédure est présentée dans le chapitre 3 où une étude analytique et numérique montre que l’erreur sur la température se comporte comme " en dehors de la zone de mélange, où " est la largeur de l’intervalle de régularisation. Cependant, à l’intérieur l’erreur se comporte comme p " et on montre que cette estimation est optimale. Cette diminution de vitesse de convergence est due à l’énergie qui reste bloquée dans la zone de mélange [58].Dans le chapitre 4 on présente quatre des schémas les plus utilisés pour le traitement de la non-linearité due au changement de phase: mise à jour du terme source, linéarisation de l’enthalpie, la capacité thermique apparente et le schéma de Chernoff. Différents tests numériques sont réalisés afin de tester et comparer ces quatre méthodes pour différents types de problèmes. Les résultats montrent que le schéma de linéarisation de l’enthalpie est le plus précis à chaque pas de temps tans dis que le schéma de la capacité thermique apparente donne de meilleurs résultats au bout d’un certain temps de calcul. Cela indique que si l’on s’intéresse aux états transitoires du matériaux le premier schéma est lemeilleur choix. Cependant, si l’on s’intéresse au comportement thermique asymptotique du matériau le second schéma est plus adapté. Les résultats montrent également que le schéma de Chernoff est le plus rapide parmi les quatre schémas en terme de temps de calcul et donne des résultats comparables à ceux des deux plus précis.Enfin, dans le chapitre 5 on utilise le schéma de Chernoff avec la méthode d’éléments finis hybride-duale Raviart-Thomas d’ordre 0 et P0 pour résoudre le problème non-linéaire de conduction thermique dans un milieu composite réel avec matériau à changement de phase. Le but étant de déterminer si un matériau composite avec une distribution uniforme de pores est assimilable à un matériau à changement de phase homogènes avec des propriétés thermo-physiques équivalentes. Pour toutes les expériences numériques exposées dans ce manuscrit on a utilisé le logiciel libre d’éléments finis FreeFem++ [41]. / In this thesis we aim to develop a numerical tool that allow to solve the unsteady heatconduction problem in a composite media with a graphite foam matrix infiltrated witha phase change material such as salt, in the framework of latent heat thermal energystorage.In chapter 1, we start by explaining the model that we are studying which is separated in three sub-parts : a heat conduction problem in the foam, a phase change problem in the pores of the foam which are filled with salt and a contact resistance condition at the interface between both materials which results in a jump in the temperature field.In chapter 2, we study the steady heat conduction problem in a composite media withcontact resistance. This allow to focus on the main difficulty here which is the treatment of the thermal contact resistance at the interface between the carbon foam and the salt. Two Finite element methods are proposed in order to solve this problem : a finite element method based on Lagrange P1 and a hybrid dual finite element method using the lowest order Raviart-Thomas elements for the heat flux and P0 for the temperature. The numerical analysis of both methods is conducted and numerical examples are given to assert the analytic results. The work presented in this chapter has been published in the Journal of Scientific Computing [10].The phase change materials that we study here are mainly pure materials and as a consequence the change in phase occurs at a single point, the melting temperature. This introduces a jump in the liquid fraction and consequently in the enthalpy. This discontinuity represents an additional numerical difficulty that we propose to overcome by introducing a smoothing interval around the melting temperature. This is explained in chapter 3 where an analytical and numerical study shows that the error on the temperature behaves like " outside of the mushy zone, where _ is the width of the smoothing interval. However, inside the error behaves like p " and we prove that this estimation is optimal due to the energy trapped in the mushy zone. This chapter has been published in Communications in Mathematical Sciences [58].The next step is to determine a suitable time discretization scheme that allow to handle the non-linearity introduced by the phase change. For this purpose we present in chapter 4 four of the most used numerical schemes to solve the non-linear phase change problem : the update source method, the enthalpy linearization method, the apparent heat capacity method and the Chernoff method. Various numerical tests are conducted in order to test and compare these methods for various types of problems. Results show that the enthalpy linearization is the most accurate at each time step while the apparent heat capacity gives better results after a given time. This indicates that if we are interestedin the transitory states the first scheme is the best choice. However, if we are interested in the asymptotic thermal behavior of the material the second scheme is better. Results also show that the Chernoff scheme is the fastest in term of calculation time and gives comparable results to the one given by the first two methods.Finally, in chapter 5 we use the Chernoff method combined with the hybrid-dual finiteelement method with P0 and the lowest order Raviart-Thomas elements to solve thenon-linear heat conduction problem in a realistic composite media with a phase change material. Numerical simulations are realised using 2D-cuts of X-ray images of two real graphite matrix foams infiltrated with a salt. The aim of these simulations is to determine if the studied composite materials could be assimilated to an equivalent homogeneous phase change material with equivalent thermo-physical properties. For all simulationsconducted in this work we used the free finite element software FreeFem++ [41].

Le principal objectif de nos travaux était l’établissement de modèles suffisamment précis pour étudier le comportement des machines électriques en présence d’un défaut de court-circuit entre spires et d’en déduire les signatures pertinentes pour la détection de ce type de défaut. L’autre objectif était de dimensionner des machines électriques à courant de court-circuit d’amplitude limitée pour réduire le risque de propagation du défaut. La première approche de modélisation consiste à effectuer une étude complète en utilisant la méthode d’éléments finis pas à pas dans le temps. Les résultats obtenus par ce modèle "éléments finis" concernant une MSAP et une MAS, saines et aussi en présence de plusieurs défauts "entre spires" de niveaux de sévérité différents, concordent avec ceux obtenus expérimentalement sur deux bancs d’essai. La seconde approche a consisté à mettre au point un modèle "circuits électriques" dont la complexité dépend du type de structure magnétique et du type de bobinage de la machine étudiée. Nous avons proposé deux méthodes de détermination des paramètres : 1- des méthodes numériques (éléments finis) ; 2- l’établissement des nouvelles expressions analytiques. Dans le dernier chapitre, une méthode basée sur la segmentation des aimants sous un pôle qui n’est en fait qu’une démultiplication du nombre de pôles au rotor sans modification du bobinage statorique est proposée est utilisé pour réduire le courant de défaut / The main objective of this research is to establish the sufficiently precise models to study the behavior of electrical machines in the presence of inter-turn short circuit fault and then find the relevant signatures to detect this type of fault. The other objective is to design a limited short-circuit current electrical machines to reduce the risk of fault development. The first modeling approach is a comprehensive study using the time stepping finite element method. The results obtained by this model "finite element" on a MSAP and MAS, healthy and faulty, for different levels of fault severity, are close with those obtained experimentally by two test benches. The second approach is to develop a model circuit electric, whose complexity depends on the type of magnetic structure and the type of machine winding. We have proposed two methods for determining the model parameters: 1 - numerical methods (FEM) which require long time bur very precise; 2 – establish new analytical expressions which is fast but less precise. In the last part, a method based on segmentation of the magnet is presented in order to reduce the short circuit current. The segmented PM motor contains the reduced fault current and can be used in the application which requires high degree of reliability

La thèse porte sur la simulation, en régime fréquentiel, du rayonnement acoustique en écoulement subsonique quelconque et dans un domaine infini. L'approche choisie s'appuie sur la résolution d'un système équivalent aux équations d'Euler linéarisées : le modèle de Galbrun. Ce modèle repose sur une représentation mixte Lagrange-Euler et aboutit à une équation dont l'unique inconnue est la perturbation du déplacement Lagrangien. Une des difficultés de l'approche de Galbrun est qu'une discrétisation directe de cette équation par une méthode d'éléments finis standard n'est pas stable. Un moyen de contourner cet obstacle est d'écrire une équation augmentée en ajoutant une nouvelle inconnue, le rotationnel du déplacement, appelée par abus vorticité. Cette approche conduit à un système qui couple une équation de type équation des ondes avec une équation de transport en régime fréquentiel. Et elle permet l'utilisation de couches parfaitement adaptées (PML) pour borner le domaine de calcul. La première partie du manuscrit est dédiée à l’étude de l’équation de transport harmonique et de sa résolution numérique, en particulier par un schéma de type Galerkin discontinu. Un des points délicats est lié au caractère oscillant des solutions de l'équation. Une fois cette étape franchie, la résolution du problème de propagation acoustique a été abordée. Une approximation basée sur l'utilisation d'éléments finis mixtes continus-discontinus avec couches parfaitement adaptées (PML) a été étudiée. En particulier, les caractères bien posés des problèmes continu et discret ainsi que la convergence du schéma numérique ont été démontrés sous certaines conditions sur l'écoulement porteur. Enfin, une mise en œuvre a été effectuée. Les résultats montrent la validité de cette approche mais aussi sa pertinence dans le cas d'écoulements complexes, voire d'écoulements dits instables / This thesis deals with the numerical simulation of time harmonic acoustic propagation in an arbitrary mean flow in an unbounded domain. Our approach is based on an equation equivalent to the linearized Euler equations called the Galbrun equation. It is derived from a mixed Eulerian-Lagrangian formulation and results in a single equation whose only unknown is the perturbation of the Lagrangian displacement. A direct solution using finite elements is unstable but this difficulty can be overcome by using an augmented equation which is constructed by adding a new unknown, the vorticity, defined as the curl of the displacement. This leads to a set of equations coupling a wave like equation with a time harmonic transport equation which allows the use of perfectly matched layers (PML) at artificial boundaries to bound the computational domain. The first part of the thesis is a study of the time harmonic transport equation and its approximation by means of a discontinuous Galerkin scheme, the difficulties coming from the oscillating behaviour of its solutions. Once these difficulties have been overcome, it is possible to deal with the resolution of the acoustic propagation problem. The approximation method is based on a mixed continuous-Galerkin and discontinuous-Galerkin finite element scheme. The well-posedness of both the continuous and discrete problems is established and the convergence of the approximation under some mean flow conditions is proved. Finally a numerical implementation is achieved and numerical results are given which confirm the validity of the method and also show that it is relevant in complex cases, even for unstable flows

L’Interaction Fluide-Structure (IFS) est un sujet d’intérêt dans beaucoup de problèmes pratiques aussi bien pour les recherches académiques ainsi que pour les applications industrielles. Différents types d’approches de simulation numérique peuvent être utilisés pour étudier les problèmes d’IFS afin d’obtenir de meilleurs conceptions et d’éviter des incidents indésirables. Dans ce travail, le domaine du fluide est simulé par une méthode hybride sans maillage (SPH-ALE), et la structure est discrétisée par la méthode d’ ´ Eléments Finis (EF). Considérant le fluide comme un ensemble de particules, on peut suivre l’interface entre le fluide et la structure d’une manière naturelle. Une stratégie de couplage conservant l’énergie est proposée pour les problèmes d’IFS transitoires où différents intégrateurs temporels sont utilisés pour chaque sous-domaine: 2nd ordre schéma de Runge-Kutta pour le fluide et schéma de Newmark pour le solide. En imposant la continuité de la vitesse normale à l’interface, la méthode proposée peut assurer qu’il n’y a ni injection d’énergie ni dissipation d’énergie à l’interface. L’énergie de l’interface est donc nulle (aux erreurs de troncature près) durant toute la période de simulation numérique. Cette méthode de couplage assure donc que la simulation de couplage est numériquement stable en temps. Les expérimentations numériques montrent que le calcul converge en temps avec l’ordre de convergence minimal des schémas utilisés dans chaque sous-domaine. Cette méthode proposée est d’abord appliquée `a un problème de piston mono-dimensionnel. On vérifie sur ce cas qu’elle ne dégrade pas l’ordre de précision en temps des schémas utilisés. On effectue ensuite les études des phénomènes de propagation d’ondes de choc au travers de l’interface fluide-structure. Un excellent accord avec la solution analytique est observé dans les cas de teste de propagation d’onde en 1-D. Finalement, les exemples multi-dimensionnels sont présentés. Ses résultats sont comparés avec ceux obtenus par d’autres méthodes de couplage. / The Fluid-Structure Interaction (FSI) effects are of great importance for many multi-physical problems in academic researches as well as in engineering sciences. Various types of numerical simulation approaches may be used to investigate the FSI problems in order to get more reliable conception and to avoid unexpected disasters. In this work, the fluid sub-domain is simulated by a hybrid mesh-less method (SPH-ALE), and the structure is discretized by the Finite Element (FE) method. As the fluid is considered as a set of particles, one can easily track the fluid structure interface. An energy-conserving coupling strategy is proposed for transient fluid-structure interaction problems where different time integrators are used for each sub-domain: 2nd order Runge-Kutta scheme for the fluid and Newmark time integrator for the solid. By imposing a normal velocity constraint condition at the interface, this proposed coupling method ensures that neither energy injection nor energy dissipation will occur at the interface so that the interface energy is rigorously zero during the whole period of numerical simulation. This coupling method thus ensures that the coupling simulation shall be stable in time, and secondly, the numerical simulation will converge in time with the minimal convergence rate of all the time integrators chosen for each sub-domain. The proposed method is first applied to a mono-dimensional piston problem in which we verify that this method does not degrade the order of accuracy in time of the used time integrators. Then we use this coupling method to investigate the phenomena of propagation of shock waves across the fluidstructure interface. A good agreement is observed between the numerical results and the analytical solutions in the 1-D shock wave propagation test cases. Finally, some multi-dimensional examples are presented. The results are compared with the ones obtained by other coupling approaches.

La conception des bâtiments selon le code sismique européen ne prend pas en compte les effets de l'interaction sol-structure (ISS). L'objectif de cette recherche est de proposer une technique de modélisation pour prendre en compte l’ISS et l'interaction structure-sol-structure (ISSS). L'approche de propagation unidirectionnelle d’une onde à trois composantes (1D-3C) est adoptée pour résoudre la réponse dynamique du sol. La technique de modélisation de propagation unidirectionnelle d'une onde à trois composantes est étendue pour des analyses d'ISS et ISSS. Un sol tridimensionnel (3D) est modélisé jusqu'à une profondeur fixée, où la réponse du sol est influencée par l’ISS et l’ISSS, et un modèle de sol 1-D est adopté pour les couches de sol plus profondes, jusqu'à l'interface sol-substrat. Le profil de sol en T est assemblé avec une ou plusieurs structures 3-D de type poteaux-poutres, à l’aide d’un modèle par éléments finis, pour prendre en compte, respectivement, l’ISS et l’ISSS dans la conception de bâtiments. La technique de modélisation 1DT-3C proposée est utilisée pour étudier les effets d’ISS et analyser l'influence d'un bâtiment proche (l'analyse d’ISSS), dans la réponse sismique des structures poteaux-poutres. Une analyse paramétrique de la réponse sismique des bâtiments en béton armé est développée et discutée pour identifier les paramètres clé du phénomène d’ISS, influençant la réponse structurelle, à introduire dans la conception de bâtiments résistants aux séismes. La variation de l'accélération maximale en haut du bâtiment avec le rapport de fréquence bâtiment / sol est tracée pour plusieurs bâtiments, chargés par un mouvement à bande étroite, excitant leur fréquence fondamentale. Dans le cas de sols et de structures à comportement linéaire, une tendance similaire est obtenue pour différents bâtiments. Cela suggère l'introduction d'un coefficient correcteur du spectre de réponse de dimensionnement pour prendre en compte l’ISS. L'analyse paramétrique est répétée en introduisant l'effet de la non-linéarité du sol et du béton armé. La réponse sismique d'un bâtiment en béton armé est estimée en tenant compte de l'effet d'un bâtiment voisin, pour un sol et des structures à comportement linéaire, dans les deux cas de charge sismique à bande étroite excitant la fréquence fondamentale du bâtiment cible et du bâtiment voisin. Cette approche permet une analyse efficace de l'interaction structure-sol-structure pour la pratique de l'ingénierie afin d'inspirer la conception d'outils pour la réduction du risque sismique et l'organisation urbaine. / Building design according to European seismic code does not consider the effects of soil-structure interaction (SSI). The objective of this research is to propose a modeling technique for SSI and Structure-Soil-Structure Interaction (SSSI) analysis. The one-directional three-component (1D-3C) wave propagation approach is adopted to solve the dynamic soil response. The one-directional three-component wave propagation model is extended for SSI and SSSI analysis. A three-dimensional (3-D) soil is modeled until a fixed depth, where the soil response is influenced by SSI and SSSI, and a 1-D soil model is adopted for deeper soil layers until the soil-bedrock interface. The T-soil profile is assembled with one or more 3-D frame structures, in a finite element scheme, to consider, respectively, SSI and SSSI in building design. The proposed 1DT-3C modeling technique is used to investigate SSI effects and to analyze the influence of a nearby building (SSSI analysis), in the seismic response of frame structures. A parametric analysis of the seismic response of reinforced concrete (RC) buildings is developed and discussed to identify the key parameters of SSI phenomenon, influencing the structural response, to be introduced in earthquake resistant building design. The variation of peak acceleration at the building top with the building to soil frequency ratio is plotted for several buildings, loaded by a narrow-band motion exciting their fundamental frequency. In the case of linear behaving soil and structure, a similar trend is obtained for different buildings. This suggests the introduction of a corrective coefficient of the design response spectrum to take into account SSI. The parametric analysis is repeated introducing the effect of nonlinear behaving soil and RC. The seismic response of a RC building is estimated taking into account the effect of a nearby building, for linear behaving soil and structures, in both cases of narrow-band seismic loading exciting the fundamental frequency of the target and nearby building. This approach allows an easy analysis of structure-soil-structure interaction for engineering practice to inspire the design of seismic risk mitigation tools and urban organization.

Ce travail présente quelques résultats concernant des méthodes numériques déterministes et probabilistes pour des équations aux dérivées partielles à coefficients aléatoires, avec des applications à l'hydrogéologie. On s'intéresse tout d'abord à l'équation d'écoulement dans un milieu poreux en régime stationnaire avec un coefficient de perméabilité lognormal homogène, incluant le cas d'une fonction de covariance peu régulière. On établit des estimations aux sens fort et faible de l'erreur commise sur la solution en tronquant le développement de Karhunen-Loève du coefficient. Puis on établit des estimations d'erreurs éléments finis dont on déduit une extension de l'estimation d'erreur existante pour la méthode de collocation stochastique, ainsi qu'une estimation d'erreur pour une méthode de Monte-Carlo multi-niveaux. On s'intéresse enfin au couplage de l'équation d'écoulement considérée précédemment avec une équation d'advection-diffusion, dans le cas d'incertitudes importantes et d'une faible longueur de corrélation. On propose l'analyse numérique d'une méthode numérique pour calculer la vitesse moyenne à laquelle la zone contaminée par un polluant s'étend. Il s'agit d'une méthode de Monte-Carlo combinant une méthode d'élements finis pour l'équation d'écoulement et un schéma d'Euler pour l'équation différentielle stochastique associée à l'équation d'advection-diffusion, vue comme une équation de Fokker-Planck. / This work presents some results about probabilistic and deterministic numerical methods for partial differential equations with stochastic coefficients, with applications to hydrogeology. We first consider the steady flow equation in porous media with a homogeneous lognormal permeability coefficient, including the case of a low regularity covariance function. We establish error estimates, both in strong and weak senses, of the error in the solution resulting from the truncature of the Karhunen-Loève expansion of the coefficient. Then we establish finite element error estimates, from which we deduce an extension of the existing error estimate for the stochastic collocation method along with an error estimate for a multilevel Monte-Carlo method. We finally consider the coupling of the previous flow equation with an advection-diffusion equation, in the case when the uncertainty is important and the correlation length is small. We propose the numerical analysis of a numerical method, which aims at computing the mean velocity of the expansion of a pollutant. The method consists in a Monte-Carlo method, combining a finite element method for the flow equation and an Euler scheme for the stochastic differential equation associated to the advection-diffusion equation, seen as a Fokker-Planck equation.

Cette recherche vise à résoudre un problème particulier du contrôle de vibration des structures intelligentes. Notre objectif est de réduire les vibrations dans une zone spécifique de la structure intelligente avec une perturbation qui couvre une large gamme de fréquences. De plus, dans cette zone spécifique, ni l'actionnement ni la détection ne sont possibles.Ici, nous faisons face à plusieurs défis principaux. Premièrement, nous devons contrôler les vibrations d’une zone spécifique de la structure, alors que nous n’avons accès aux mesures que dans d’autres zones. Deuxièmement, la large bande passante de la perturbation implique que nombreux modes doivent être contrôlés au même temps, ce qui nécessite l'utilisation de plusieurs actionneurs et capteurs. Cela conduit à un contrôleur MIMO difficile à obtenir avec les méthodes classiques de conception de contrôleur. Troisièmement, il faut éviter le problème de propagation, qui consiste à garantir la stabilité en boucle fermée lorsque le contrôleur basé sur un modèle est appliqué à la configuration réelle. Pour relever ces défis, nous étudions deux stratégies de contrôle: le contrôle centralisé et le contrôle distribué.Pour le contrôle centralisé, nous proposons une méthodologie qui nous permet d’obtenir un contrôleur MIMO simple permettant de relever ces défis. Tout d'abord, plusieurs techniques de modélisation et d’identification sont appliquées pour obtenir un modèle précis d'ordre faible de la structure intelligente. Ensuite, une méthode de synthèse basée sur le contrôle H_∞ avec un critère H_∞ particulièrement proposé est appliquée. Ce critère H_∞ intègre plusieurs objectifs de contrôle, y compris les défis principaux. En particulier, le problème de débordement se transforme en un problème de stabilité robuste et sera garanti en utilisant ce critère. Le contrôleur H_∞ obtenu est une solution standard du problème H_∞. Le contrôleur final est obtenu en simplifiant ce contrôleur H_∞ sans perdre la stabilité en boucle fermée ni dégrader les performances. Cette méthodologie est validée sur une structure de poutre avec des transducteurs piézoélectriques et la zone centrale est celle où les vibrations devraient être réduites. L'efficacité du contrôleur obtenu est validée par des simulations et des expériences.Pour le contrôle distribué, on considère la même structure de poutre et les mêmes objectifs de contrôle. Il existe des méthodes visant à concevoir des contrôleurs distribués pour les systèmes spatialement interconnectés. Cette recherche propose une méthode basée sur la FEM, associée à plusieurs techniques de réduction de modèle, permettant de discrétiser spatialement la structure de poutre et d'en déduire les modèles d’espace d'état des sous-systèmes interconnectés. La conception des contrôleurs distribués ne sera pas abordée dans cette recherche. / This research aims at solving a particular vibration control problem of smart structures. We aim at reducing the vibration in a specific zone of the smart structure under the disturbance that covers a wide frequency band. Moreover, at this specific zone, neither actuation nor sensing is possible.Here we face several main challenges. First, we need to control the vibration of a specific zone of the structure while we only have access to measurements at other zones. Second, the wide bandwidth of the disturbance implies that numerous modes should be controlled at the same time which requires the use of multiple actuators and sensors. This leads to a MIMO controller which is difficult to obtain using classical controller design methods. Third, the so-called spillover problem must be avoided which is to guarantee the closed-loop stability when the model-based controller is applied on the actual setup. To tackle these challenges, we investigate two control strategies: the centralized control and the distributed control.For centralized control, we propose a methodology that allows us to obtain a simple MIMO controller that accomplishes these challenges. First, several modeling and identification techniques are applied to obtain an accurate low-order model of the smart structure. Then, an H_∞ control based synthesis method with a particularly proposed H_∞ criterion is applied. This H_∞ criterion integrates multiple control objectives, including the main challenges. In particular, the spillover problem is transformed into a robust stability problem and will be guaranteed using this criterion. The obtained H_∞ controller is a standard solution of the H_∞ problem. The final controller is obtained by further simplifying this H_∞ controller without losing the closed-loop stability and degrading the performance. This methodology is validated on a beam structure with piezoelectric transducers and the central zone is where the vibration should be reduced. The effectiveness of the obtained controller is validated by simulations and experiments.For distributed control, we consider the same beam structure and the same control objectives. There exist methods aiming at designing distributed controllers of spatially interconnected system. This research proposes a FEM based method, combined with several model reduction techniques, that allows to spatially discretize the beam structure and deduce the state-space models of interconnected subsystems. The design of distributed controllers will not be tackled in this research.