Cette étude est consacrée à la modélisation du comportement thermo-mécanique des argiles raides saturées et au développement d'un algorithme d'intégration efficace de contrainte correspondant. Le comportement mécanique de l'argile de Boom naturelle dans des conditions isothermes a été caractérisé. Le modèle Cam Clay modifié (MCC) a été ensuite appliquée pour simuler le comportement de l'argile de Boom naturel. Il a été constaté que le MCC donne des prédictions de mauvaise qualité pour le comportement de l'argile de Boom naturel. Ainsi, un modèle Cam Clay (ACC-2) adapté a été développé en introduisant une nouvelle surface de charge et un nouveau potentiel plastique ainsi que d'un mécanisme plastique de Deux surfaces. Ce modèle permet la description satisfaisante des caractéristiques principales du comportement mécanique de l'argile de Boom naturelle. De plus, les équations de ce modèle peuvent être formulées mathématiquement comme dans un modèle élasto-plastique classique. L'algorithme d'intégration de contrainte classique peut donc être appliqué. Les effets thermiques ont été examinés par l'évaluation de la pertinence de trois lois thermomécaniques avancées (Cui et al, 2000; Abuel-Naga et al, 2007; Laloui et François, 2008; 2009). Il apparaît que tous les trois modèles peuvent décrire les caractéristiques principales du comportement thermo-mécanique des argiles saturées. Cependant, chaque modèle a ses limites ou des points peu clairs du point de vue théorique. L'algorithme d'intégration de contrainte du modèle thermo-mécanique de Cui et al. (2000) au point de contrainte a également été développé spécifiquement en utilisant une méthode adaptive du pas de temps. Le temps de calcul nécessaire pour obtenir une précision donnée est ainsi largement réduit pour des chemins de chargements thermiques et mécaniques. Un modèle thermo-mécanique à Deux surfaces (modèle TEAM) a été développé en se basant sur le mécanisme plastique de Deux surfaces. Le modèle proposé a étendu le modèle de Cui et al. (2000) à une formulation de Deux surfaces considérant le couplage entre les déformations plastiques des chemins de chargements thermiques et mécaniques. La simulation des essais drainés montre que ce modèle peut décrire les caractéristiques principales thermo-mécaniques de l'argile de Boom naturelle le long de différents chemins de chargements. Le modèle TEAM a finalement été étendu à des conditions non drainées. Après la clarification du concept des contraintes effectives et la définition d'une condition de déformation volumique, le processus d'échauffement non drainé est analysé. La validité des équations thermo-hydro-mécaniques de ce modèle a été examinée en se basant sur des résultats d'essais typiques / This study is devoted to the thermo-mechanical constitutive modeling for saturated stiff clays and the development of a corresponding efficient stress integration algorithm. The mechanical behavior of natural Boom Clay in isothermal conditions was first characterized. The Modified Cam Clay model (MCC) was then applied to simulate the natural Boom Clay behavior. It has been found that the MCC gives poor-quality predictions of the natural Boom Clay behavior. Thereby, an adapted Cam Clay model (ACC-2) was developed by introducing a new yield surface and a new plastic potential as well as a Two-surface plastic mechanism. This model allows satisfactory prediction of the main features of the mechanical behavior of natural Boom Clay. Moreover, the constitutive equations of this model can be formulated mathematically as in a classic elasto-plastic model. Thus, the classic stress integration algorithm can be applied. The thermal effects were considered by assessing the performance of some advanced thermo-mechanical models (Cui et al., 2000; Abuel-Naga et al., 2007; Laloui and François, 2008; 2009). It appears that all the three models can capture the main features of the thermo-mechanical behavior of saturated clays. However, each constitutive model has its own limitations or unclear points from the theoretical point of view. The stress integration algorithm of the thermo-mechanical model proposed by Cui et al. (2000) at the stress point level was also developed using a specifically designed adaptive time-stepping scheme. The computation time required to achieve a given accuracy is largely reduced with the adaptive sub-stepping considered for both mechanical and thermal loadings. A Two-surface thermo-mechanical model (TEAM model) was developed based on the Two-surface plastic mechanism. The proposed model extends the model of Cui et al. (2000) to a Two-surface formulation, considering the plastic strain coupling between the thermal and the mechanical loading paths. The simulation of drained tests shows that this model can capture the main thermo-mechanical features of natural Boom Clay along different loading paths. The TEAM model was finally extended to undrained conditions. After setting up an appropriate effective stress principle and defining a volumetric strain condition, the undrained heating process was analyzed. The validity of the thermo-hydro-mechanical constitutive equations was examined based on the data from typical tests