The goal of this work is to establish a life prediction methodology for thermomechanical loading of the Ni-base superalloy Mar-M247 over a larger temperature range than previous work. The work presented in this thesis extends the predictive capability of the Sehitoglu-Boismier unified thermo-viscoplasticity constitutive model and thermomechanical life model from a maximum temperature of 871C to a maximum temperature of 1038C. The constitutive model, which is suitable for predicting stress-strain history under thermomechanical loading, is adapted and calibrated using the response from isothermal cyclic experiments conducted at temperatures from 500C to 1038C at different strain rates with and without dwells. In the constitutive model, the flow rule function and parameters as well as the temperature dependence of the evolution equation for kinematic hardening are established. In the elevated temperature regime, creep and stress relaxation are critical behaviors captured by the constitutive model. The life model accounts for fatigue, creep, and environmental-fatigue damage under both isothermal and thermomechanical fatigue. At elevated temperatures, the damage terms must be calibrated to account for thermally activated damage mechanisms which change with increasing temperature. At lower temperatures and higher strain rates, fatigue damage dominates life prediction, while at higher temperatures and slower strain rates, environmental-fatigue and creep damage dominate life prediction. Under thermomechanical loading, both environmental-fatigue and creep damage depend strongly on the relative phasing of the thermal and mechanical strain rates, with environmental-fatigue damage dominating during out-of-phase thermomechanical loading and creep damage dominating in-phase thermomechanical loading. The coarse-grained polycrystalline microstructure of the alloy studied causes a significant variation in the elastic response, which can be linked to the crystallographic orientation of the large grains. This variation in the elastic response presents difficulties for both the constitutive and life models, which depend upon the assumption of an isotropic material. The extreme effects of a large grained microstructure on the life predictions is demonstrated, and a suitable modeling framework is proposed to account for these effects in future work.

O presente trabalho aborda o estudo do comportamento de túneis em maciços de solo. É dada ênfase na aplicação de análises numéricas com emprego de diferentes modelos constitutivos elasto-plásticos para solos. São apresentadas análises numéricas tridimensionais de um túnel amplamente instrumentado pertencente à Companhia do Metropolitano de São Paulo. As análises são realizadas com auxílio de um programa computacional de elementos finitos. O comportamento do maciço em pontos situados ao redor da escavação é minuciosamente estudado e a capacidade dos modelos em representar adequadamente o comportamento verificado na obra pelas instrumentações é avaliada. Além das análises numéricas o trabalho aborda os principais conceitos relacionados com escavações de túneis em maciços de solo. Conceitos relacionados com a engenharia prática de túneis são apresentados de maneira qualitativa, sem formulações teóricas e matemáticas. Uma revisão bibliográfica com publicações recentes das principais revistas e periódicos que tratam do tema de análise numérica aplicada a túneis é apresentada. O trabalho também apresenta uma revisão dos principais conceitos relacionados com os modelos constitutivos comumente utilizados para análise de problemas de geotecnia. Além do modelo elástico são apresentados os modelos elasto-plásticos de Tresca, von Misses, Drucker-Prager e Mohr-Coulomb. Uma breve introdução aos conceitos básicos de estado crítico, juntamente com as equações do modelo Cam-Clay original e Cam-Clay modificado são apresentadas. Antes da apresentação das equações desses modelos constitutivos, são introduzidos os conceitos básicos relacionados com o comportamento dos materiais elasto-plásticos. Os conceitos de material elástico perfeitamente plástico e de material com endurecimento (ou hardening) e amolecimento (ou softening) são apresentados. Os conceitos de superfície de plastificação e de superfície de potencial plástico também são apresentados. Por fim, são sintetizados os pontos mais relevantes da pesquisa realizada, apontando as limitações do trabalho com sugestões de novos estudos a serem realizados nessa mesma linha de pesquisa. / The present research approaches the study of the behaviour of tunnels in soil. It is given emphasis in the application of numerical analyses using different elasto-plastic constitutive models for soils. Three-dimensional numerical analyses of a widely instrumented tunnel belonging to the Company of the Metropolitan of São Paulo are presented. The analyses are carried through with aid of a computational program of finite elements. The behaviour of the soil mass in points located around the excavation is thoroughly studied and the capacity of the models in adequately representing the field behavior verified by the instrumentations is evaluated. Furthermore, the work approaches the main concepts related to tunneling in soils. Concepts related to practical engineering of tunnels are presented in a qualitative way, without theoretical and mathematical formulations. A literature review of recent publications of the most important periodic magazines and that deal with the subject of numerical analysis applied to tunnels is presented. The work also presents a revision of the main concepts related to the constitutive models normally used for analysis of geotechnical problems. Beyond the elastic model the elasto-plastics models of Tresca, von Misses, Drucker-Prager and Mohr-Coulomb are presented. Brief introductions to the basic concepts of critical state, together with the equations of the (original) Cam-Clay original and (modified) Cam-Clay modified models are presented. Before the presentation of the equations of these constitutive models, the basic concepts of the behaviour of the elasto-plastics materials are introduced. The concepts of perfectly plastic elastic material and material with hardening and softening are presented. The concepts of plastic surface and plastic potencial surface are also presented. Finally, the most relevant points of the research are synthesized, pointing the limitations of the developed work along with suggestions for new studies to be carried through in this line of research.

Les sols non saturés jouent un rôle important dans de nombreux domaines. En Génie Civil (stabilité des talus, fondations superficielles), en Ingénierie de l'environnement (transport de polluants par les eaux souterraines, enfouissement de déchets radioactifs), en Ingénierie des Réservoirs (subsidence suite aux pompages dans les champs pétroliers ou dans les nappes aquifères), etc. L'objet de la thèse est l'étude du comportement hydromécanique des sols fin non saturés, tels les limons et les argiles. Les phénomènes typiques du comportement de ces sols sont d'abord présentés avant d'aborder l'étude de leurs lois constitutives. Divers modèles sont passés en revue et un modèle unifié original est développé. Ensuite, une vaste gamme d'essais de laboratoire est analysée. Ils portent sur des sols saturés, quasi saturés et non saturés et concernent quatre matériaux: l'argile de Kruibeke, le limon de Gembloux, l'argile de Boom et la bentonite d'Almeria. Cette partie de la recherche conduit à une critique détaillée des procédures classiques d'essais de laboratoires sur les sols supposés saturés. Elle fournit également une base de données importante pour la calibration de la loi constitutive unifiée développée plus haut. Enfin, des simulations numériques, réalisées en aveugle, d'expériences n'ayant pas servi à la calibration permettent de valider ce modèle.

This thesis is dedicated to the study of structural and petroleum geology problems. To this purpose, a frictional elastoplastic law based on the Van Eekelen criterion is formulated, which avoids the classical drawbacks of the Drücker Prager criterion. Also, a 2D automatic adaptive re-meshing algorithm is developed for complex multidomains configurations, in order to overcome the limitation of the Lagrangian mesh. Details of the hydromechanical formulation implemented in the LAGAMINE FE code in a large strain context are presented. Application cases (reproduction of sandbox simulation, study of a hydrocarbon trap) concentrate on the study of the strain localisation and potential fracturation using different criteria. In the first case, re-meshing technique allowed to reproduce successfully analogue experiment of thrusting propagation. In the second case, a detailed study based on different initial conditions has brought new insight to the reactivated origin of some faults and has allowed to obtain information on the potential fracturing of the hydrocarbon reservoir unit. As an academic case, the study of anorthosite diapirism is carried out, which confirms the validity of the petrological model of diapirism. Eventually, the hydromechanical coupling effects between a layered porous medium and a fault are illustrated on a simple case.

La mayoría de las estructuras de hormigón armado se someten a solicitaciones combinadas de esfuerzos axiles, flexión, cortante y torsión. La fisuración del hormigón, plastificación de las armaduras y otros efectos no-lineales hacen que las secciones transversales de estos elementos presenten un comportamiento anisótropo que deriva en el acoplamiento de los esfuerzos normales y tangenciales. Es decir, esfuerzos normales o momentos flectores pueden producir deformaciones de corte y vice versa. Aunque en algunas ocaciones, esta interacción es considerada de forma simplificada en el dimensionamiento de estructuras, hasta el momento no se ha realizado un análisis profundo de los efectos acoplados en secciones de forma arbitraria bajo cargas 3D generales utilizando modelos de fibras.El objetivo principal de esta tesis es generalizar el análisis de secciones de hormigón armado mediante fibras, de forma que se pueda reproducir la res-puesta no-lineal acoplada frente a esfuerzos normales y tangenciales bajo solicitaciones tridimensionales generales. De igual forma, se pretende obtener, para los esfuerzos cortantes y torsión, la misma capacidad de representación de geometrías y combinación de materiales que ofrecen los modelos de fibras para esfuerzos de flexo-compresión.La primera problemática estriba en representar adecuadamente la cinemática de la sección transversal. Con la excepción de las deformaciones normales contenidas en el plano de la sección, no existe una teoría cinemática que a priori pueda dar la distribución del resto de deformaciones o tensiones en la sección, sin dejar de satisfacer las condiciones de equilibrio interno o continuidad entre las fibras que componen la misma.Por otra parte, para materiales anisótropos, como el hormigón fisurado, en general todos los esfuerzos internos pueden estar acoplados. Además, es preciso considerar la distorsión de la sección transversal para satisfacer el equilibrio entre fibras.El problema se aborda de forma general, considerando una sección de forma y materiales cualesquiera. Se parte del problema diferencial de equilibrio de un sólido con el que se ha podido deducir un sistema de equilibrio entre fibras (equilibrio a nivel sección). Se puede demostrar que éste es complementario al problema estándar de vigas. El sistema complementario permite recuperar información tridimensional que normalmente se pierde al resolver un problema de vigas.Posteriormente, se propone una solución interna del problema complementario, en la que el alabeo y la distorsión de la sección quedan expresados como una función de las deformaciones generalizadas de una viga: deformaciones axil y cortantes, curvaturas de flexión y torsión. No son necesarios grados de libertad adicionales a nivel estructura ni hipótesis a-priori sobre la forma de los campos de deformación o tensión interna.A partir de la formulación teórica, se desarrolla un modelo de elementos finitos plano de la sección transversal. El modelo está preparado para servir como respuesta constitutiva de cualquier tipo de elemento viga en sus puntos de integración. %Se evita así la necesidad de realizar un modelo de elementos sólidos de toda la barra para estudiar la respuesta frente a una combinación general de esfuerzos normales y tangenciales.Se implementan una serie de modelos constitutivos para distintos materiales. En particular, se implementa un modelo constitutivo triaxial para hormigón fisurado, considerando la anisotropía inducida por la fisuración e incluyendo la superficie de rotura según un criterio multiaxial.La formulación seccional es validada mediante varios casos de estudio teóricos y experimentales. La respuesta no-lineal acoplada bajo diversas combinaciones de esfuerzos normales y tangenciales es reproducida con precisión, lo cual queda patente tanto en las curvas esfuerzo-deformación obtenidas como en las matrices de rigidez seccionales.Finalmente, se recopilan las conclusiones derivadas de la presente investigación y seofren recomendaciones para futuros trabajos. / Most RC structures are subjected to combined normal and tangential forces, such as bending, axial load, shear and torsion. Concrete cracking, steel yielding and other material nonlinearities produce an anisotropic sectional response that results in a coupling between the effects of normal and shear forces, i.e. normal force or bending moments may produce shear strains and vice versa. Although this interaction is sometimes taken into account, in a simplified manner, in the design of RC structures, a deep analysis of the coupling effects of RC sections using fiber models has not yet been made for arbitrary shape sections under general 3D loading.The main objective of this thesis is to generalize the fiber-like sectional analysis of reinforced concrete elements, to make it capable of considering the coupled non-linear response under tangential and normal internal forces, from a general 3D loading.Similarly, it is desired to obtain, for torque and shear forces, the same capacity and versatility in reproducing complex geometry and materials combination that fiber-like sectional representations offers for bending and stretching.The first problematic lies in finding a proper representation of the section's kinematics under such general loading. Except for in-plane normal strains, there is no single kinematical theory capable of a-priori representing the correct distribution of the others strains or stresses satisfying, at the same time, inter-fiber equilibrium and continuity. On the other hand, for rather anisotropic materials, such as cracked concrete, all internal forces are, in general, coupled. It is also required that distortion is allowed for the section's kinematics in order to guarantee satisfaction of internal equilibrium.The problem is dealt in a general form considering arbitrary shaped sections and any material behaviour. Starting from the differential equilibrium of a solid, an inter-fiber equilibrium system (equilibrium at the sectional level) was deduced. This system, which is complementary to the standard equilibrium problem of a beam-column, allows to recuperate information of the three-dimensional problem that is generally lost when solving a beam problem.Further, a solution of the equilibrium at the sectional level is proposed in which the section's warping and distortion are posed as a function of the generalized beam-column strains (axial and shear strains, bending and torsion curvatures). No additional degrees of freedom are required at the structural level nor a-priori hypotheses on the distribution of the internal strains or stresses.After the theoretical formulation, a planar finite element model for cross-sectional analysis is developed. The model can be used as a constitutive law for general beam column elements at their integration points.A series of constitutive models have been implemented for several materials. In particular, a triaxial constitutive model for cracked concrete is implemented considering crackinduced anisotropy and a multiaxial failure criterion.The sectional formulation is validated by means of various theoretical and experimental case studies. Non-linear coupled response under normal and tangential internal forces is reproduced with accuracy, as can be seen both in the predicted internal force-strain curves and in the sectional stiffness matrixes.Finally, the conclusions drawn from the current research are summarized andrecomendations for future works are given.

A combined geological and rock mechanics approach to tunnel face behaviour prediction, based on improved understanding of brittle fracture processes during TBM excavation, was developed to complement empirical design and performance prediction for TBM tunnelling in hard rock geological conditions. A major challenge of this research was combining geological and engineering terminology, methods, and objectives to construct a unified Geomechanical Characterisation Scheme. The goal of this system is to describe the spalling sensitivity of hard, massive, highly stressed crystalline rock, often deformed by tectonic processes. Geological, lab strength testing and TBM machine data were used to quantify the impact of interrelated geological factors, such as mineralogy, grain size, fabric and the heterogeneity of all these factors at micro and macro scale, on spalling sensitivity and to combine these factors within a TBM advance framework. This was achieved by incorporating aspects of geology, tectonics, mineralogy, material strength theory, fracture process theory and induced stresses. Three main approaches were used to verify and calibrate the Geomechanical Characterisation Scheme: geological and TBM data collection from tunnels in massive, highly-stressed rock, interpretation of published mineral-specific investigations of rock yielding processes, and numerical modelling the rock yielding processes in simulated strength tests and the TBM cutting process. The TBM performance investigation was used to identify the mechanism behind the chipping processes and quantify adverse conditions for chipping, including tough rock conditions and stress induced face instability. The literature review was used to identify the critical geological parameters for rock yielding processes and obtain strength and stiffness values for mineral-specific constitutive models. A texture-generating algorithm was developed to create realistic rock analogues and to provide user control over geological characteristics such as mineralogy, grain size and fabric. This methodology was applied to investigate the TBM chipping process to calibrate the Geomechanical Characterisation Scheme. A Chipping Resistance Factor was developed to combine the quantified geological characteristic factors and laboratory strength values to predict conditions with high risk of poor chipping performance arising from tough rock. A Stress-Related Chip Potential Factor was developed to estimate conditions with high risk of advance rate reduction arising from stress-induced face instability. / Thesis (Ph.D, Geological Sciences & Geological Engineering) -- Queen's University, 2008-09-25 23:58:58.071

Design of landfills must consider both stability and integrity of the lining system. Therefore, stresses and strains in both mineral and geosynthetic lining materials must be controlled. Interaction between waste and barrier system is of particular importance for assessing the stability and structural integrity of steep non-self supporting barrier systems. The most appropriate approach to assess the interaction is the use of numerical modelling techniques, and therefore an appropriate constitutive model for waste material is required to represent its mechanical behaviour. In a literature review the key aspects of mechanical behaviour of municipal solid waste (MSW) were investigated, including the influence of compressible and reinforcing particles on compression and shear behaviour of MSW were identified. Constitutive modelling of both MSW and soil material were reviewed, based on which the methodology for this study have been developed. In addition, requirements of an appropriate constitutive model for MSW have been suggested from the numerical modelling experience, and a framework to develop a constitutive model for MSW was produced. A one-dimensional compression model was developed by including the influence of compressible particles on MSW compression behaviour. One-dimensional compression tests on both real and synthetic waste samples were modelled and the results have shown that the compression model can reproduce the measured behaviour. A fibre reinforcing model was developed by including the influence of reinforcing particles on MSW shear behaviour. A triaxial compression test on fibre reinforced sand was modelled and the results have shown that the reinforcing model can predict its shear strength. A constitutive model for MSW has been developed by combining the Modified Cam-Clay with the one-dimensional compression and the fibre reinforcing models. Typical MSW triaxial compression tests have been modelled and the results have shown that the MSW model can reproduce the stress-strain behaviour in specific strain ranges. The constitutive model for MSW has been coded into a non-linear elasto-plastic finite element method program. Comparisons between the finite element analysis results and the analytical solutions have been performed and good agreements have been obtained.

O presente trabalho aborda o estudo do comportamento de túneis em maciços de solo. É dada ênfase na aplicação de análises numéricas com emprego de diferentes modelos constitutivos elasto-plásticos para solos. São apresentadas análises numéricas tridimensionais de um túnel amplamente instrumentado pertencente à Companhia do Metropolitano de São Paulo. As análises são realizadas com auxílio de um programa computacional de elementos finitos. O comportamento do maciço em pontos situados ao redor da escavação é minuciosamente estudado e a capacidade dos modelos em representar adequadamente o comportamento verificado na obra pelas instrumentações é avaliada. Além das análises numéricas o trabalho aborda os principais conceitos relacionados com escavações de túneis em maciços de solo. Conceitos relacionados com a engenharia prática de túneis são apresentados de maneira qualitativa, sem formulações teóricas e matemáticas. Uma revisão bibliográfica com publicações recentes das principais revistas e periódicos que tratam do tema de análise numérica aplicada a túneis é apresentada. O trabalho também apresenta uma revisão dos principais conceitos relacionados com os modelos constitutivos comumente utilizados para análise de problemas de geotecnia. Além do modelo elástico são apresentados os modelos elasto-plásticos de Tresca, von Misses, Drucker-Prager e Mohr-Coulomb. Uma breve introdução aos conceitos básicos de estado crítico, juntamente com as equações do modelo Cam-Clay original e Cam-Clay modificado são apresentadas. Antes da apresentação das equações desses modelos constitutivos, são introduzidos os conceitos básicos relacionados com o comportamento dos materiais elasto-plásticos. Os conceitos de material elástico perfeitamente plástico e de material com endurecimento (ou hardening) e amolecimento (ou softening) são apresentados. Os conceitos de superfície de plastificação e de superfície de potencial plástico também são apresentados. Por fim, são sintetizados os pontos mais relevantes da pesquisa realizada, apontando as limitações do trabalho com sugestões de novos estudos a serem realizados nessa mesma linha de pesquisa. / The present research approaches the study of the behaviour of tunnels in soil. It is given emphasis in the application of numerical analyses using different elasto-plastic constitutive models for soils. Three-dimensional numerical analyses of a widely instrumented tunnel belonging to the Company of the Metropolitan of São Paulo are presented. The analyses are carried through with aid of a computational program of finite elements. The behaviour of the soil mass in points located around the excavation is thoroughly studied and the capacity of the models in adequately representing the field behavior verified by the instrumentations is evaluated. Furthermore, the work approaches the main concepts related to tunneling in soils. Concepts related to practical engineering of tunnels are presented in a qualitative way, without theoretical and mathematical formulations. A literature review of recent publications of the most important periodic magazines and that deal with the subject of numerical analysis applied to tunnels is presented. The work also presents a revision of the main concepts related to the constitutive models normally used for analysis of geotechnical problems. Beyond the elastic model the elasto-plastics models of Tresca, von Misses, Drucker-Prager and Mohr-Coulomb are presented. Brief introductions to the basic concepts of critical state, together with the equations of the (original) Cam-Clay original and (modified) Cam-Clay modified models are presented. Before the presentation of the equations of these constitutive models, the basic concepts of the behaviour of the elasto-plastics materials are introduced. The concepts of perfectly plastic elastic material and material with hardening and softening are presented. The concepts of plastic surface and plastic potencial surface are also presented. Finally, the most relevant points of the research are synthesized, pointing the limitations of the developed work along with suggestions for new studies to be carried through in this line of research.

The assessment of ageing tunnels requires a deeper understanding of the long-term behaviour of twin tunnels, whilst lack of permeability data limits the accuracy of long-term predictions. This thesis therefore investigates long-term twin-tunnel behaviour through finite-element parametric analyses, and provides additional pereability data through laboratory studies. Permeability tests are performed on fissured London Clay, exploring the effect of isotropic stress cycles on the permeability of fissures. A model explaining the permeability-stress relationship is proposed to explain irrecoverable changes observed in fissure permeability, and is formulated mathematically for numerical implementation. Laboratory investigations are performed on grout from the London Underground tunnels, investigating permeability, porosity, microstructure and composition. A deterioration process is proposed to explain observations, consisting of acid attack and leaching. The deterioration had appeared to transform the grout from impermeable to permeable relative to the soil. The change in grout permeability with time would strongly influence long-term movements. The long-term behaviour of single tunnels is investigated in a finite-element parametric study. A new method is formulated to predict long-term horizontal and vertical surface displacements after excavation of a single tunnel, and incorporates an improved measure of relative soil-lining permeability. The study also predicts significant surface movements during the consolidation period, contradicting the lack of further building damage observed in the field. A further parametric study also investigates the long-term behaviour of twin tunnels. Key interaction mechanisms are identified, leading to the postulation of the long-term interaction behaviour under different tunnelling conditions. Long-term interaction is found to be complex and significant, and should be accounted for in numerical simulations.

Le comportement thermomécanique des géomatériaux est un sujet de recherche qui touche un nombre important de domaines d'application : le stockage des déchets radioactifs en couche géologique profonde, le comportement des structures géothermiques, le stockage de chaleur, le comportement des matériaux à proximité des câbles à haute tension, le comportement saisonnier des matériaux asphaltiques, l'origine des tremblements de terre, etc. Dans le contexte du stockage des déchets radioactifs en couches géologiques profondes, un intérêt particulier a récemment été accordé à l'étude du comportement thermomécanique de l'argile de Boom, l'argile à Opalines, l'argilite du Callovo-Oxfordien et la diorite d'Äspö. L'ensemble de ces études montre qu'une exposition à des températures élevées peut altérer les propriétés physiques et mécaniques des géomatériaux : l'angle de frottement, la cohésion, la perméabilité-porosité, les modules élastiques, la résistance, le comportement volumique, le comportement post-pic, le comportement à long terme, etc. Depuis une trentaine d'années, un nombre important de modèles de comportement ont été développés dans le but d'intégrer les effets de la température à des modèles élasto-plastiques dérivés de la théorie de l'état critique. Cependant, peu de ces modèles intègrent les effets du temps sur le comportement thermomécanique des géomatériaux ; lequel apparaît essentiel au dimensionnement de structures à long terme. L'objectif de ces travaux de thèse est de répondre à cette problématique en intégrant les effets de la température au modèle de comportement mécanique L&K. Le modèle de comportement L&K est un modèle à deux mécanismes, l'un dit "plastique" permettant de décrire le comportement mécanique instantané des matériaux, l'autre dit "viscoplastique" permettant de décrire le comportement différé. Après un état des lieux détaillé concernant les comportements thermomécaniques instantané et différé des géomatériaux, les effets de la température ont été pris en compte à travers l'évolution des paramètres d'écrouissage des deux mécanismes. Ce modèle thermomécanique a été validé sur des applications expérimentales en support de ce travail de thèse. Les résultats montrent que la formulation proposée permet de reproduire fidèlement le comportement thermomécanique des géomatériaux. Les applications à venir concernent des calculs sur ouvrage dans le contexte du stockage des déchets radioactifs / The effect of temperature on the behaviour of geomaterials is a crucial issue in geotechnical and underground engineering. The thermo-mechanical behaviour of rocks and soils contains many applications in the fields of high-level nuclear waste disposal, heat storage, geothermal structures, petroleum drilling, zones around buried high-voltage cables, bituminous materials, and geological research. In the context of nuclear waste disposal at great depths, the thermo-mechanical behaviour of Boom clay, Opalinus clay, Callovo-Oxfordian argillite and Äspö diorite has recently received special attention in Europe. Research in these areas has demonstrated that rocks and soils may suffer from changes in their mechanical properties during short-to long-term exposure to an elevated temperature. These changes include effects on the friction angle, permeability/porosity, elastic moduli, shear strength, dilatancy, softening, brittle-to-ductile transition, creep, etc... Since Prager's first works on the modelling of non-isothermal plastic deformation, many constitutive models have been developed to include these phenomena in computational inelasticity. Most models generalize the critical-state model to include the effects of temperature on the short-term behaviour of clays and rocks. However, the effect of time on the thermo-mechanical behaviour is often not coupled to the rate-independent plasticity even if the long-term behaviour appears essential for ensuring the safety and stability during the design and construction analysis in many fields, such as the storage of nuclear waste and more generally underground excavations. Thus, it appears important to combine both instantaneous and delayed thermo-mechanical effects to obtain appropriate constitutive equations to model such problems. In this Ph.D thesis, a unified thermo-plastic/viscoplastic constitutive model has been developed for this purpose. This model is a straightforward extension of the unified elasto-plastic/viscoplastic L&K constitutive model which was developed in previous Ph.D works. The updated thermo-mechanical model includes the evolutions of the two yield limits and the fluidity coefficient with temperature. The model was validated under several thermo-mechanical conditions on clayey rocks. The typical features of the thermo-mechanical behaviour of geomaterials were well reproduced. The numerical predictions of the triaxial compression tests and creep tests clearly indicate that the model can predict the overall behaviour of geomaterials under deviatoric and non-isothermal stress paths