The range of problems that geotechnical engineers must face is increasing in complexity and scope. Often, complexity arises from the interaction between the soil and the environment – the topic of this lecture. To deal with this type of problem, the classical soil mechanics formulation is progressively generalised in order to incorporate the effects of new phenomena and new variables on soil behaviour. Recent advances in unsaturated soil mechanics are presented first: it is shown that they provide a consistent framework for understanding the engineering behaviour of unsaturated soils, and the effects of suction and moisture changes. Building on those developments, soil behaviour is further explored by considering thermal effects for two opposite cases: high temperatures, associated with the problem of storage and disposal of high-level radioactive waste; and low temperatures in problems of freezing ground. Finally, the lecture examines some issues related to chemical effects on soils and rocks, focusing in part on the subject of tunnelling in sulphate-bearing rocks. In each case new environmental variables are identified, enhanced theoretical formulations are established, and new or extended constitutive laws are presented. Particular emphasis is placed on mechanical constitutive equations, as they are especially important in geotechnical engineering. The lecture includes summary accounts of a number of case histories that illustrate the relevance and implications of the developments described for geotechnical engineering practice.

A bonded elasto-plastic soil model is adopted to simulate self-boring pressuremeter tests. Calibration of the model with results from several natural clays identifies a meaningful parameter range for the simulations. Then a systematic sensitivity analysis of the loading SBPM test curve is carried out with an infinite cylindrical cavity analogue, focusing on the variables and parameters that are specific to the bonded material. This analysis reveals that the effects of mechanical overconsolidation and bonding on the pressuremeter loading curves are very similar. Unloading and strain-holding phases of the test show a greater capability to separate mechanical overconsolidation from bonding. The limited effects of finite pressuremeter length and partial drainage on the simulated curves do not change those conclusions.

The Yucca Mountain Drift Scale Test (DST) is a multiyear, large-scale underground heating test designed to study coupled thermal–hydrological–mechanical–chemical behavior in unsaturated fractured and welded tuff. As part of the international cooperative code-comparison project DEvelopment of COupled models and their VALidation against EXperiments, four research teams used four different numerical models to simulate and predict coupled thermal–hydrological–mechanical (THM) processes at the DST. The simulated processes included heat transfer, liquid and vapor water movements, rock-mass stress and displacement, and stress-induced changes in fracture permeability. Model predictions were evaluated by comparison to measurements of temperature, water saturation, displacement, and air permeability. The generally good agreement between simulated and measured THM data shows that adopted continuum model approaches are adequate for simulating relevant coupled THM processes at the DST. Moreover, thermal-mechanically induced rock-mass deformations were reasonably well predicted using elastic models, although some individual displacements appeared to be better captured using an elasto-plastic model. It is concluded that fracture closure/opening caused by change in normal stress across fractures is the dominant mechanism for thermal-stress-induced changes in intrinsic fracture permeability at the DST, whereas fracture shear dilation appears to be less significant. This indicates that such changes in intrinsic permeability at the DST, which are within one order of magnitude, are likely to be mostly reversible.

This paper presents thermo-hydro-mechanical (THM) analyses simulating the Drift Scale Test (DST) performed at Yucca Mountain. A double structure approach based on two superimposed domains is adopted. Intrinsic permeability changes with deformations imply full THM coupling. Temperatures and gas permeabilities were measured during 4 years and are used to validate the model. Measured gas permeability variations show patterns that are successfully explained by the model calculations. These gas permeability variations may be attributed to thermo-hydraulic effects, and also to mechanical effects. Different cases of intrinsic permeability variations have been considered in the model and their influence on the calculated temperatures, degree of saturations and gas permeabilities are presented. Volumetric deformation, in contraction or dilatancy, implies changes in the aperture of rock fractures that in turn lead to changes in intrinsic permeability. Dilatancy, caused by shear stresses, increases intrinsic permeability. Consideration of this factor contributes significantly to improve the agreement of calculated gas permeability with the measured values obtained during the DST experiment.

The paper presents the results of an experimental study on the effects of suction and temperature on the volumetric behaviour of an unsaturated soil. Statically compacted Boom clay samples have been used in the testing programme, which has been performed using an oedometer apparatus with simultaneous control of suction and temperature. Suctions up to 0·45 MPa and temperatures up to 80°C have been applied. Special measures have been taken to control and minimise the problems associated with phase changes and vapour migration at high temperatures. A comprehensive experimental programme has been carried out, including isothermal tests with control of vertical stress, suction and temperature, non-isothermal tests under constant suction, and isothermal swelling pressure tests. Testing samples at two compacted soil densities provides observations concerning a wide scope of behaviour ranging from collapse to swelling behaviour on wetting. Test results obtained at 80°C have been compared with those of tests performed at room temperature (22°C). Based on these comparisons, thermal effects on a variety of features of behaviour such as swelling potential, swelling pressure, collapsibility, shrinkage strains, compressibility, stress path dependence and degree of strain reversibility are presented and discussed. Although the basic behaviour is similar at high and low temperatures, various significant differences are identified concerning mainly swelling strains developed during wetting and compressibility with respect to loading. Swelling pressure test results indicate that the high-temperature samples yield at a lower stress than that observed in room-temperature samples, suggesting that the stress domain bounded by the main yield surface reduces as temperature increases—a behaviour feature shared with saturated soils.

The paper presents an elasto-plastic model for unsaturated soils that takes explicitly into account the mechanisms with which suction affects mechanical behaviour as well as their dependence on degree of saturation. The proposed model is formulated in terms of two constitutive variables directly related to these suction mechanisms: the average skeleton stress, which includes the average fluid pressure acting on the soil pores, and an additional scalar constitutive variable, ξ, related to the magnitude of the bonding effect exerted by meniscus water at the inter-particle contacts. The formulation of the model in terms of variables closely related to specific behaviour mechanisms leads to a remarkable unification of experimental results of tests carried out with different suctions. The analysis of experimental isotropic compression data strongly suggests that the quotient between the void ratio, e, of an unsaturated soil and the void ratio e_s, corresponding to the saturated state at the same average soil skeleton stress, is a unique function of the bonding effect due to water menisci at the inter-particle contacts. The same result is obtained when examining critical states at different suctions. Based on these observations, an elasto-plastic constitutive model is developed using a single yield surface the size of which is controlled by volumetric hardening. In spite of this simplicity, it is shown that the model reproduces correctly many important features of unsaturated soil behaviour. It is especially remarkable that, although only one yield surface is used in the formulation of the model, the irreversible behaviour in wetting–drying cycles is well captured. Because of the behaviour normalisation achieved by the model, the resulting constitutive law is economical in terms of the number of tests required for parameter determination.

The paper examines the interaction between host rock and a clay-based engineered barrier in the context of deep underground disposal of high-level radioactive waste. A large scale “in situ” heating test currently under way in the underground Grimsel Rock Laboratory (GRL; Switzerland) is adopted as a representative case for study. The main features of the behaviour of the rock/barrier system have been examined using a comprehensive programme of coupled thermo–hydro–mechanical analyses. A parametric study has been performed to assess the effects of individual phenomena and parameters. The characterisation activities carried out in relation with the in situ test provide the information on material properties and test conditions required. Results highlighting the influence of vapour diffusion, permeability, and retention curve of rock and bentonite are presented. The interaction between clay barrier and host medium exhibits a high degree of complexity.

The estimation of soil and rock parameters based on field instrumentation data is a common procedure in geomechanics. The use of system identification and optimization techniques allows the performance of this type of analyses in a more rational and objective manner. In this paper a probabilistic formulation for the backanalysis problem is presented. The procedure described involves the evaluation of the measurement covariance matrices, which are derived for some geotechnical instruments used in field instrumentation. The algorithm used to solve the mathematical problem of optimization is also presented, as well as its coupling to a finite element code. The algorithm requires the computation of the sensitivity matrix, which can be evaluated “exactly” in terms of the finite element method. Finally, a synthetic example, based on the excavation of a tunnel, is presented in which the elastic modulus E and the K0 parameter of the material are identified from measured displacements. The effect of the number of measurements and their error structure is also discussed.

A general statistical framework to perform backanalysis in geotechnical problems from field instrumentation has been presented in a companion paper. Here, an application to a real case involving the excavation of a tunnel in stiff overconsolidated clay is described. Both, extensometer and inclinometer measurements are used as input data and elastic moduli of the ground and the value of the K0 coefficient are estimated. The finite element method is used as the computational procedure to solve the direct problem, and has been coupled to the identification algorithm as described in the companion paper. In addition, a discussion on the reliability of the parameters identified is presented.