Ce travail de recherche se situe dans le cadre du développement de la méthode des éléments de frontière (BEM) pour les milieux poreux multiphasiques. À l'heure actuelle, l'application de la BEM aux pr oblèmes des milieux poreux non-saturés est encore limitée, car l'expression analytique exacte de la solution fondamentale n'a pas été obtenue, ni dans le domaine transformé ni dans le domaine réel. Ceci provient de la complexité du système d'équations régissant le comportement des milieux poreux non-saturés. Les développements de la BEM pour les sols non-saturés effectués au cours de cette thèse sont basés sur les modèles thermo-hydro-mécanique (THHM) et hydro-mécanique (HHM) présentés dans la première partie de ce mémoire. Ces modèles phénoménologiques basés sur la théorie de la poromécanique et les acquis expérimentaux sont obtenus dans le cadre du modèle mathématique présenté par Gatmiri (1997) et Gatmiri et al. (1998). Après avoir présenté les modèles THHM et HHM, on établit pour la première fois les équations intégrales de frontière et les solutions fondamentales associées pour un milieu poreux non-saturé sous chargement quasi-statique pour les deux cas isotherme (2D dans le domaine de Laplace) et non-isotherme (2D et 3D dans les domaines de Laplace et temporel). Aussi, les équations intégrales de frontière ainsi que les solutions fondamentales 2D et 3D (dans le domaine de Laplace) pour le modèle dynamique couplé des sols non-saturés sont obtenues. Ensuite, les formulations d'éléments de frontière (BEM) basées sur la méthode quadrature de convolution (MQC) concernant les milieux poreux saturé et non-saturé sous chargement quasi-statique isotherme et dynamique sont implémentées dans le code de calcul « HYBRID ». Ayant intégrées les formulations de BEM pour les problèmes de propagation d'ondes ainsi que pour les problèmes de consolidation dans les milieux poreux saturés et non-saturés, il semble que nous ayons fourni à l'heure actuelle le premier code de calcul aux éléments de frontière (BEM) qui modélise les différents problèmes dans les sols secs, saturés et non-saturés. Une fois le code vérifié et validé, des études paramétriques portant sur des effets de site sismiques sont effectuées. Le but recherché est d'aboutir à un critère simple, directement exploitable par les ingénieurs, combinant les caractéristiques géométriques et les caractéristiques du sol, permettant de prédire l'amplification du spectre de réponse en accélération dans des vallées sédimentaires aussi bien que vides / The purpose of this dissertation is to develop a boundary element method (BEM) for multiphase porous media. Nowadays, the application of the BEM for solving problems of unsaturated porous media is still limited, because no fundamental solution exists in the published literature, neither in the frequency nor time domain. This fact rises from the complexity of the coupled partial differential equations governing the behaviour of such media. The developments of the BEM for the unsaturated soils carried out during this thesis are based on the thermo-hydro-mechanical (THHM) and hydro-mechanical (HHM) models presented in the first part of this dissertation. These phenomenological models are presented based on the experimental observations and with respect to the poromechanics theory within the framework of the suction-based mathematical model presented by Gatmiri (1997) and Gatmiri et al. (1998). After having presented the THHM and HHM models, for the first time, one establishes the boundary integral equations (BIE) and the associated fundamental solutions for the unsaturated porous media subjected to quasi-static loading for both isothermal (2D in the Laplace transform domain) and non-isothermal (2D and 3D in Laplace transform and time domains) cases. Also, the boundary integral equations as well as the fundamental solutions (2D and 3D in the Laplace transform domain) are obtained for the fully coupled dynamic model of unsaturated soils.In the next step, the boundary element formulations (BEM) based on the convolution quadrature method (CQM) regarding the saturated and unsaturated porous media subjected to isothermal quasi-static and dynamic loadings are implemented via the computer code HYBRID. Having integrated the BEM formulations for the wave propagation, as well as the consolidation problems in the saturated and unsaturated porous media, it seems that now the first boundary element code is obtained that can model the various problems in dry, saturated and unsaturated soils. Once the code is verified and validated, parametric studies on seismic site effects are carried out. The aim is to achieve a simple criterion directly usable by engineers, combining the topographical and geological characteristics of the soil, to predict the amplification of acceleration response spectra in sedimentary as well as hollow valleys