Search

The structural assessment of segmental grey cast iron (GCI) tunnel linings to nearby construction is challenging due to the presence of the joints affecting the stiffness of the tunnel lining. This paper presents an extensive investigation, using 3D finite element (FE) analyses, into the bending moment-rotation (M-θ) behaviour of two GCI tunnel joint geometries. These two geometries correspond to standard running and station tunnels of the London Underground network. The contribution of this study is two-fold. i) The novel characterisation of the M-θ response enables the development of new models for simulating the mechanical response of GCI tunnel joints with structural elements which can be used in simplified, 2D geotechnical analysis for tunnel safety assessments. ii) The analyses provide insight into the behaviour of GCI tunnel linings that would be difficult to achieve through experimental and field observations alone. More specifically, the analyses show that when the bolts are removed from the joints the possibility of tensile failure can be disregarded; that the initial bolt preload influences the rotational stiffness only after some rotation has taken place and does not alter the bending moment of opening; and that the out-of-plane displacement restraint has little influence on the joint response.

Soil–vegetation–atmosphere interaction has long been known to induce significant pore pressure variations at shallow depths and associated superficial slope movements. Recent findings suggest that the effect of this interaction may also extend to large depths in natural clay slopes. Multiple examples of weather-induced deep landslide mechanisms can be found in the Southern Apennines (Italy), where slopes are often formed of fissured clays. The relationship between the activity of these landslides and the hydro-mechanical processes due to soil–vegetation–atmosphere interaction was investigated herein by means of a two-dimensional coupled hydro-mechanical finite element analysis. A constitutive model capable of simulating the behaviour of highly overconsolidated clays, in both saturated and unsaturated states, was adopted in the analysis, in conjunction with a boundary condition capable of reproducing the combined effects of rainfall infiltration, evapo-transpiration, and run-off. The results of the analysis corroborate the connection between weather conditions, pore pressure variations, and slope movements in natural clay slopes. The importance of adequately reproducing the geological history of a natural slope to define its current state is also demonstrated.

Seawalls are constructed by covering and protecting the sloping seashore ground with concrete plates or blocks. Their purpose is to sustain high waves induced by strong winds and prevent ground erosion, but they often collapse, mobilising different modes of failure, including that of the ground. Nevertheless, limited research has been conducted on ground failure caused by high waves. In this study, a series of novel centrifuge model tests was first conducted to investigate the failure mechanisms of seawalls due to wave propagation, focusing on the failure of the ground. Finite-element analyses were subsequently conducted to explore the failures observed in the model tests and to provide further insight as to the state of the ground leading to failure. Two failure modes were demonstrated to prevail: floating of the covering panel and sliding failure of the ground. In addition, of the possible causes of failure, the following three were identified in the current study: increased unit weight and reduced suction from wetting; enhanced seepage force under the panel and around the toe block during backwash; water pressure on the back of the panel and the landward side of the toe block during backwash.

In this paper, 2-dimensional, hydro-mechanically coupled finite element analyses are conducted to assess the performance of an engineered barrier, constructed from natural geomaterials, aimed at reducing flood risk in urban environments. The barrier consists of an unsaturated compacted soil layer with water holding properties and a drainage layer of a coarse granular material, that acts as a capillary break, and is constructed on top of the natural soil, in this case London clay. The barrier is vegetated so that its water storage capacity is renewed after each rainfall event. Sophisticated boundary conditions are used to simulate the effect of precipitation and evapotranspiration. The evolution of the rainfall infiltration and runoff rate is simulated both for a treated soil column with an engineered barrier and an untreated one consisting solely of in-situ London Clay. The percolation rate of rainfall water from the bottom of the barrier is also estimated. This comparison highlights the effectiveness of the engineered barrier in reducing the risk of fast flooding, in preventing excessive deformations and in protecting underground infrastructure during wetting and drying cycles. The effect of the hydraulic properties and geometry of the barrier is investigated by means of an extensive parametric analysis. Finally, recommendations for the design of barrier systems are made.

Soil-atmosphere interaction has been attracting increasing interest as the seasonal variation of pore water pressures (pwp) has been linked to a variety of geotechnical problems (e.g. slope stability and serviceability, foundation subsidence or swelling, desiccation cracking etc.) or has been identified as part of the solution of geotechnical problems (e.g. in sustainable urban drainage systems). Prediction of how the pwp will change within soils of low permeability under the combined effect of evapotranspiration and precipitation requires adequate knowledge of the soil permeability and how it varies spatially (e.g with depth) and temporally (e.g. with suction or degree of saturation, void ratio or due to the opening and closing of desiccation cracks). Nonetheless, in-situ measurements of permeability that satisfy both the spatial and temporal variation are difficult. In order to clarify the importance of variable permeability in predicting pwp variations under atmospheric loads, a series of one- and two-dimensional finite element analyses was performed, where the permeability model and the variation of permeability were parametrically studied. The results demonstrated that the variation of permeability, as well as the model employed in the analysis, e.g. allowing or not for desiccation cracking, influenced the values of suction calculated as well as the pwp profile with depth, highlighting the importance of estimating the spatial and temporal variation of permeability with some level of confidence.

It has long been established that pore water pressure (pwp) variations affect the stability and serviceability of slopes. Pwp variations may be due to consolidation/swelling processes within soils of low permeability or may be due to seasonal evapotranspiration and precipitation processes. The simultaneous study of such phenomena and of hydro-mechanical coupling in sloping ground requires use of advancednumerical methods. Often, infrastructure slopes are considered in plane strain (two-dimensional) conditions for simplicity and to date this is the case for most numerical analyses considering soil-atmosphere interaction. However, this approach makes it impossible to study the longitudinal extent of a possible slip surface forming following vegetation removal. Three-dimensional, fully-coupled numerical analyses of acut slope were performed herein to study the effect of vegetation clearance on stability and explore numerically ways of implementing effective vegetation management.

The structural response of a segmental grey cast iron (GCI) tunnel lining ring under distortion was investigated by means of finite element (FE) analysis. Building on previous experimental investigations, a 3D numerical model, capable of reproducing accurately the behaviour observed in the laboratory, was developed with the aim of providing guidelines for the structural assessment of GCI linings in engineering practice. A comprehensive validation of the segmental ring model with the laboratory data was first completed. Subsequently, a parametric study was conducted using a set-up that replicated the widely adopted elastic continuum method, so that differences between the numerical and the analytical solution could be attributed to the presence of the longitudinal joints. In this manner, the influence of the joints on the ring response was quantitatively established and recommendations for routine engineering calculations developed. A set of bending stiffness reduction factors are proposed as a function of the tunnel ovalisation, providing upper and lower limits of the bending stiffness, as well as a global reduction factor which is an average measure of the bending stiffness reduction. These factors can be integrated into the calculation procedure of closed-form solutions in order to account for the segmental nature of GCI linings.

Atmospheric phenomena such as rainfall and evapotranspiration contribute to slope movements in unsaturated soils, the study of which requires fully coupled numerical methods, combined with realistic boundary conditions and appropriate mechanical and hydraulic soil properties. This paper focuses on the effect of the hydraulic behaviour, and in particular of the modelling of the soil-water retention curve and the permeability on slope movements, with the aim of identifying which model parameters are critical and, therefore, require careful experimental identification.

The stability of cut slopes is greatly influenced by seasonal pore-water pressure variations under the combined effect of rainfall and vegetation. However, predicting soil–atmosphere interaction is not straightforward, due to the complexity of both the boundary conditions involved and the hydromechanical behaviour of soils, which is coupled and highly nonlinear, rendering the use of numerical tools, such as finite element analysis, necessary. This paper discusses the numerical modelling of soil–atmosphere interaction and presents the analysis of a slope cut in London clay in a highly vegetated area. The whole life cycle of the slope is considered with phases of low and high water demand vegetation and vegetation clearance. The analysis results indicate that dense vegetation is associated with high factors of safety, but may induce large differential displacements, which are likely to affect the serviceability of the slope. Vegetation clearance, however, may initiate instability, highlighting the need for effective vegetation management to achieve a balance between serviceability and ultimate limit states. Although the case considered is representative of southeast England, it introduces the necessary tools for realistic numerical analysis of soil–atmosphere interaction.

A temperature rise in soils is usually accompanied by an increase in excess pore fluid pressure due to the differential thermal expansion coefficients of the pore fluid and the soil particles. To model the transient behaviour of this thermally induced excess pore fluid pressure in geotechnical problems, a coupled THM formulation was employed in this study, which accounts for the non-linear temperature-dependent behaviour of both the soil permeability and the thermal expansion coefficient of the pore fluid. Numerical analyses of validation exercises (where there is an analytical solution), as well as of existing triaxial and centrifuge heating tests on Kaolin clay, were carried out in the current paper. The obtained numerical results exhibited good agreement with the analytical solution and experimental measurements respectively, demonstrating good capabilities of the applied numerical facilities and providing insight into the mechanism behind the observed evolution of the thermally induced pore fluid pressure. The numerical results further highlighted the importance of accounting for the temperature-dependent nature of the soil permeability and the thermal expansion coefficient of the pore fluid, commonly ignored in geotechnical numerical analysis.

Interface elements are frequently employed in finite element (FE) analyses to represent soil-structure interfaces or rock joints. The modelling of coupled thermo-hydro-mechanical (THM) problems in geotechnical engineering requires equally a coupled and robust THM formulation for interface elements. This paper presents such a formulation which is capable of reproducing the coupled THM behaviour of discontinuities and soil-structure boundaries, and is compatible with other types of finite elements used to discretise the soil and structural domains (e.g. solid and shell elements). The coupled THM three-dimensional (3D) zero-thickness interface element is implemented into the bespoke FE code employed in this research and its features are verified using a number of numerical exercises. To demonstrate their performance, the proposed interface elements are employed in the simulation of the coupled THM behaviour of a fissured triaxial sample subjected to a thermal load and the influence of the presence of fissures on soil behaviour is presented.

In this paper, 2-dimensional, hydro-mechanically coupled finite element analyses are performed to assess the performance of an engineered barrier aimed at reducing flood risk in urban environments. The barrier consists of an unsaturated compacted soil layer and a drainage layer of a coarse granular material, constructed on top of the natural soil, in this case London clay. The barrier is vegetated so that its water storage capacity is renewed after each rainfall event. Sophisticated boundary conditions are used to simulate the effect of precipitation and evapotranspiration. The rainfall water infiltration and the initiation of water run-off during intense precipitation events are simulated. The effect of the hydraulic properties of the unsaturated soil layer on the performance of the system is investigated by means of parametric analyses. The effect of precipitation rate and geometry of the barrier is also discussed. Design recommendations for the properties of the compacted layer and the dimensions of the system are given at the end of the paper.

Buffer materials for nuclear waste disposal applications generally consist of blocks made of highly expansive compacted clay. However, high-density pellets of bentonite are being evaluated as an alternative buffer material for waste isolation. The material response of pellet-based buffers may be quite different from that of compacted buffers, because of the peculiar discontinuous porosity presented. An overview of the literature available on pellet-based buffers is presented and, in particular, two main topics are discussed: firstly, the characteristics of the fabric of the pellets that can be observed through techniques of micro-structural investigation, secondly, the most important behavioural features that can be seen during material testing. Additionally, the constitutive frameworks that have already been developed specifically for pellets are also reviewed. The overall objective of the paper is to highlight the differences between compacted and pellet-based bentonite buffers, in order to propose suitable assumptions to start developing a constitutive model for the latter.

The structural testing and finite element (FE) analysis described in this paper were part of a major research project undertaken at Imperial College London to investigate the deformation of bolted segmental grey cast iron (GCI) tunnel linings. A key aim was to quantify how joints influence the behaviour of the lining, through a three-path approach comprising physical experiments, finite element modelling, and field instrumentation. The laboratory results have been used to assess the validity of the tunnel assessment methods used by industry. This study examined joint articulation under the serviceability limit state in the absence of hoop force focussing on factors such as applied bolt preload, the loading direction and the freedom of the circumferential flange to deflect. Two half-scale GCI lining segments were bolted together at the longitudinal flanges to form a bolted arch in a similar fashion to the tests performed by Thomas (1977). Modern instrumentation was implemented to gain detailed measurements quantifying changes in global displacements of the two segments, bolt forces and joint opening under applied loading. For the first time, the physical experiments were conducted contemporaneously with the development of a three-dimensional FE model of the joint. The experimental data and the results from the FE analysis indicate a reduction in joint stiffness as the joint articulates under applied load. It is shown that the presence of a joint has far greater influence on the behaviour of the ‘arch’ than the level of preload applied to the bolts in the joint. The FE analysis allowed the deformation behaviour of the joint under positive and negative bending to be investigated: its response under the two modes differs significantly.

This paper describes a calibration procedure for the double structure constitutive model ICDSM (Imperial College Double Structure Model), developed for highly expansive clays, when the model is applied to MX-80 bentonite. Firstly, the model calibration process is discussed and organised in a number of hierarchical steps. These steps involve the estimation of the macrostructural parameters that can be derived from oedometer, isotropic and triaxial laboratory data. Estimation of the microstructural parameters is more challenging due to the limited knowledge of an expansive clay’s fabric and of the physico-chemical phenomena that control its evolution upon wetting. Nevertheless, this paper discusses the available sources of data and identifies the appropriate information that is needed to characterise the micro-structural behaviour of the bentonite. Finally, through the simulation of a swelling pressure test on a bentonite plug, the hydration of the material is studied as a hydro-mechanical coupled process. Particular attention is devoted to the evolution of the stress state of the sample, which is compared to the experimental measurements in order to demonstrate that the constitutive model accurately reproduces the expansive behaviour of MX-80 bentonite.

Constitutive models for unsaturated soils are most commonly based on conventional critical state type models for saturated soils. The ellipse of the modified Cam-clay model is usually adopted for the formulation of the yield and the plastic potential surfaces in the mean stress – deviatoric stress plane. Despite the wide use of this ellipse it has long been shown that the peak deviatoric stress observed for highly overconsolidated states may be largely overestimated. For fully saturated conditions, a planar surface, termed the Hvorslev surface, has been employed in the literature in order to overcome this shortcoming on the dry side of the critical state. In addition, by employing a no-tension cut-off, stress states are restricted within the compressive region. In unsaturated conditions and due to the effect of suction the soil may be subjected to tensile total stresses which should not, however, exceed a model dependent limit. Nonetheless, adopting a planar surface and in the absence of any provision, illegal tensile stresses may be obtainable. For this reason, a new expression for the Hvorslev surface, which restricts the stress state within acceptable boundaries, is proposed, based on a small number of model parameters. The new surface is curved rather than planar and can also be used in combination with critical state type constitutive models developed for saturated conditions. The new surface is calibrated based on fully saturated undrained triaxial compression and extension tests and on unsaturated drained triaxial tests. The superiority of the new surface in comparison with commonly adopted shapes for the yield and the plastic potential functions – such as the ellipse and the bullet shape of the modified and the original Cam clay models – is demonstrated through numerical analyses of the unsaturated drained triaxial laboratory tests.

While the evolution with time of the factor of safety against slope failure is well understood for excavations performed under fully saturated conditions, this evolution has not been extensively studied for excavations carried out in unsaturated soils. The objective of this paper is to investigate the relationship between the factor of safety and time for an excavation performed in an unsaturated silty soil, employing the finite element method. A hypothetical boundary value problem is considered and two types of analyses performed; in the first one, unsaturated soil behaviour is modelled through appropriate constitutive and soil–water retention curve models, whereas in the second one full saturation is assumed. The effect that the saturated soil permeability, its variation with suction, the increase of apparent cohesion due to suction, the depth of the groundwater table and the hydraulic hysteresis have on the results of the unsaturated analysis is examined in a parametric study. The analyses results demonstrate that for unsaturated soils the factor of safety may increase with time, in contrast to what is commonly accepted to be the case in fully saturated soils. Furthermore, it is not possible to know in advance which one of the two types of analysis performed (i.e. the saturated or the unsaturated) will produce conservative results. It is, therefore, advisable to perform unsaturated analysis in geotechnical practice, when dealing with excavations in such soils.

The paper presents the Finite Element formulation of the equations proposed by Tsiampousi et al. (2016) for coupled consolidation in unsaturated soils. Their coupling is discussed in relation to a conceptual model which divides soil behaviour into zones ranging from fully saturated to dry states. The numerical simulation of a laboratory experiment involving drainage of water from a vertical column of sand is used to validate the equations. Finally, the example of rainfall infiltration into a cut slope highlights how aspects of the conceptual model are reflected in the numerical analysis of boundary value problems involving unsaturated soils.

The equations governing coupled consolidation in unsaturated soils are known to contain additional parameters when compared to the equations for saturated soils. Nonetheless, the variation of these parameters with suction or degree of saturation is not generally agreed upon. The paper introduces a novel approach to deriving general equations for each of these parameters and their variation, and explains that, for consistency with the constitutive and soil-water retention curve models adopted, these general equations need to be transformed into case-specific expressions. Finally, a conceptual model is presented highlighting how the behaviour of unsaturated soil reflects aspects of its water content.

A technique for controlling relative humidity (RH) is presented, which involves supplying a sealed chamber with a continuous flow of air at a computer-regulated RH. The desired value of RH is achieved by mixing dry and wet air at appropriate volumes and is measured for servo-control at three locations in the chamber with capacitive RH sensors and checked with a sensitive VAISALA sensor. The setup is capable of controlling RH steadily and continuously with a deviation of less than 0.2% RH. The technique was adopted to determine wetting soil-water retention curves (SWRC) of statically compacted London Clay, under both free-swelling and constant volume conditions. The RH within the chamber was increased in a step-wise fashion, with each step maintained until vapour equilibrium between the chamber atmosphere and the soil samples was established. Independent filter paper measurements further validate the method, while the obtained retention curves complement those available in the literature for lower ranges of suction.

Deep and slow landslide processes are frequently observed in clay slopes located along the Southern Apennines (Italy). A case study representative of these processes, named Pisciolo case study, is discussed in the paper. The geo-hydro-mechanical characteristics of the materials involved in the instability phenomena are initially discussed. Pluviometric, piezometric, inclinometric and GPS monitoring data are subsequently presented, suggesting that rainfall infiltration constitutes the main factor inducing slope movements. The connection between formation of landslide bodies and slope-atmosphere interaction has been demonstrated through a hydro-mechanical finite element analysis, whose results are finally reported in the work. This analysis has been conducted employing a constitutive model that is capable of simulating both saturated and unsaturated soil behaviour, as well as a boundary condition able to simulate the effects of the soil-vegetation-atmosphere interaction.

© The authors and ICE Publishing: All rights reserved, 2015.Experimental and numerical studies of the behaviour of expansive clays have been attracting increasing interest, due to their good sealing properties, which render them ideal to be used as engineered barriers (buffers) in both active (e.g. nuclear) and non-active waste disposal facilities. Both large scale and laboratory scaled experiments indicate that the sealing capabilities of the buffer are fundamentally governed by its volumetric behaviour when wetted. In this paper, a constant volume column infiltration test, performed under isothermal conditions on compacted MX80 bentonite, is modelled numerically using the Imperial College Finite Element Program (ICFEP). A modified version of the Barcelona Basic Model is used to simulate the behaviour of the buffer, which is inherently partly saturated. The numerical results agree well with the observed experimental data, especially with regard to the advancement of the wetting front. A detailed interpretation of the computed evolutions with time of stress state, suction and void ratio at different elevations along the samples axis is carried out, providing insight into the complex hydro-mechanical response of the buffer during the experiment. Indeed, even though the overall volume of the sample was kept constant, a region of localised dilation, which induced the contraction of other zones of the material, was observed to advance simultaneously with the wetting front along the height of the soil column.

Interaction between atmosphere and soil has only recently attracted significant interest. Soil-atmosphere interaction takes place under dynamic climatic conditions, which vary throughout the year and are expected to suffer considerable alterations due to climate change. However, Geotechnical Analysis has traditionally been limited to simplistic approaches, where winter and summer pore water pressure profiles are prescribed. Geotechnical Structures, such as cut slopes, are known to be prone to large irreversible displacements under the combined effect of water uptake by a complex vegetation root system and precipitation. If such processes take place in an unsaturated material the complexity of the problem renders the use of numerical analysis essential. In this paper soil-atmosphere interaction in cut slopes is studied using advanced, fully coupled partially saturated finite element analyses. The effect of rainfall and evapotranspiration is modelled through sophisticated boundary conditions, applying actual meteorological data on a monthly basis. Stages of low and high water demand vegetation are considered for a period of several years, before simulating the effect of vegetation removal. The analysis results are presented with regard to the serviceability and stability of the cut slope.

One of the most important features in unsaturated soil mechanics is the soil-water retention curve and its coupling to the mechanical component of soil behaviour. It has long been recognised that the retention curve exhibits significant hysteresis, and that it is affected by the specific volume. Several attempts have been made in the past to model this behaviour. A novel approach is proposed herein, which accounts for both the hydraulic hysteresis and the specific volume dependence of the retention relationship in a three-dimensional formulation. The primary and the scanning paths are simple geometric curves, which have a common tangent at the point of intersection, ensuring a smooth transition from scanning to primary paths. A small number of parameters are required to define the primary paths, and no fitting parameters are necessary for generation of the scanning paths. As knowledge of the specific volume and its variation is required, the retention model needs to be employed in conjunction with a constitutive model capable of reproducing the complex behaviour of unsaturated soils. To guarantee consistency with the retention model, the degree of saturation needs to be incorporated in the specific volume–suction relationship adopted within the constitutive model. To accommodate such a feature when absent, a new expression for the soil compressibility with suction as a function of the degree of saturation is proposed. Simulations of laboratory experiments on unsaturated soils, involving cyclic changes of applied suction, demonstrate the effectiveness of the proposed modelling approach.

Despite the advances in both understanding and modelling of the behaviour of unsaturated soils, numerical analyses considering partial saturation are rare in engineering practice. The purpose of the current paper is to compare the factors of safety against slope failure obtained from different types of analyses (unsaturated and effective stress analyses), using the finite element method, and to investigate the degree of inaccuracy introduced when partial saturation is neglected. The results of the numerical analyses demonstrate that, although similar factors of safety are calculated for steady-state conditions, significant disparity is exhibited in the short-term.