A nonlinear elastic approach to modelling the hydro-mechanical behaviour of the SEALEX experiments on compacted MX-80 bentonite

Hydraulic seals using compacted sand–bentonite blocks are an important part of the closure phase of deep geological disposal facilities for the isolation of many categories of radioactive wastes. An understanding of the hydro-mechanical behaviour of these seals and the ability to model their behaviour is a key contribution to safety cases and licence applications. This work reports the development of a hydro-mechanically coupled model and its application to the simulation of a range of test conditions investigated in the SEALEX experiments conducted by IRSN at Tournemire URL. The work has been conducted as part of the recently completed DECOVALEX-2015 project. Richards’ equation for unsaturated fluid flow is coupled to a nonlinear elastic strain-dependent mechanical model that incorporates a moving finite element mesh, and calibrated against laboratory experiments. Stress and volumetric dependencies of the water retention behaviour are incorporated through the Dueck suction concept extended to take into account permanent changes in water retention behaviour during consolidation. Plastic collapse in laboratory results is modelled with the application of a source term activated by a threshold defined in terms of the net axial stress and net suction. The model is used to simulate both a 1/10 scale mock-up laboratory test and full-scale in situ performance test and is capable of reproducing the major trends in the data with just nine mechanical parameters and an experimentally defined stress threshold.


Introduction
Deep geological disposal facilities (GDF) comprising engineered and geological barriers are internationally considered as the most feasible option for long-term isolation and disposal of radioactive wastes (Kim et al. 2011). Regional groundwater flow could provide a mechanism for radionuclides to be transported from disposed waste to the biosphere, so it is fundamental to many disposal concepts that a GDF is effectively sealed. Hydraulic seals, designed to perform this function in current research and development programmes, are primarily composed of unsaturated compacted blocks of bentonite (Andra 2005). This is an expansive clay that swells to fill engineered voids on exposure to water, and has low permeability to ensure diffusion-dominated flow and advantageous sorptive properties to inhibit radionuclide migration (Hansen et al. 2013). The specific role of bentonite in a GDF may vary from concept to concept, e.g. in addition to hydraulic seals, bentonite is also proposed as a buffer material around high-level waste canisters (e.g. Oy 2010) and as a backfill material for excavation drifts (e.g. NAGRA 2009;Oy 2010).
In order to support safety case and licence applications, a thorough understanding of bentonite behaviour and predictive capability is required. The engineering lifetime of a GDF is expected to be in excess of 100,000 years, and, as such, numerical modelling is required in order to provide confidence in the safety performance of a hydraulic seal. Therefore, development of numerical models for the simulation and prediction of bentonite behaviour under the range of conditions relevant to its application within a GDF is an active area of research (Navarro et al. 2014(Navarro et al. , 2015Saba et al. 2014a, b;Ferrari and Seiphoori 2015).
Numerical models of unsaturated bentonite behaviour under a range of conditions typically couple Richards' equation or multiphase flow to a chosen mechanical model, such as nonlinear elasticity (Fredlund and Rahardjo 1993;Cui et al. 2002) or elasto-plasticity (Gens and Alonso 1992;Wheeler and Sivakumar 1995). This is achieved through the dependence of hydraulic parameters on the mechanical stresses and has led to the extension of mechanical elastoplastic models such as the Cam Clay and Mohr-Coulomb models to include the dependence of soil suction or saturation in the mechanical formulation (Gens and Alonso 1992;Börgesson et al. 1995;Wheeler and Sivakumar 1995;Alonso and Vaunat 1999;Sánchez et al. 2005;Navarro et al. 2014;Zhou and Sheng 2015).
In recent years, elasto-plastic models, such as the Barcelona Basic Model (BBM) or the Barcelona Expansive Model (BExM), have become predominant in the simulation of unsaturated soils and expansive materials (Alonso et al. 1990;Gens and Alonso 1992;Alonso and Vaunat 1999;Navarro et al. 2014). Elasto-plastic approaches based on different combinations of stress variables (e.g. Wheeler and Sivakumar 1995) and generalised plasticity models addressing the double structure of bentonite have also been proposed (Sánchez et al. 2005). The BBM and BExM are extensions of the modified Cam Clay model that incorporate soil suction in the mechanical solution through an additionally defined failure surface, known as the loadingcollapse curve (Alonso et al. 1990;Gens and Alonso 1992;Alonso and Vaunat 1999). The BExM is an extension of the BBM in which the effect of the micro-and macrostructure of bentonite on the mechanical behaviour has been addressed by additional fitting functions (Gens and Alonso 1992;Alonso and Vaunat 1999). The effect of the double structure on hydraulic behaviour has also recently been included in these formulations through parameterisation of two water retention curves for each structural scale (Alonso et al. 2011;Navarro et al. 2015).
However, the increasing complexity of the models has led to increasing numbers of parameters used for calibration. For example, parameterisation of the BBM for a constant volume test in Navarro et al. (2015) requires 14 mechanical parameters, and Rutqvist et al. (2011) report 18 mechanical parameters in the implementation of the BBM in TOUGH-FLAC. While this provides flexibility in reproducing experimental data, the final parameterisation may not be unique and may require an extensive experimental programme to calibrate the model to. Consequently, the predictive potential of models with large parameter sets depends on the completeness of the calibration data set.
In this paper, we present a comparatively simple model based on nonlinear elasticity coupled to Richards' equation for unsaturated fluid flow. The motivation of this study is to investigate an alternative method to that currently used, in an effort to minimise the parameter set required to simulate the general hydro-mechanical behaviour of bentonite. Nonlinear mechanical behaviour is simulated with an incremental function in which the mechanical properties of the bentonite are dependent on the strain history of the sample. An Updated Lagrangian finite element mesh is developed within the open source code OpenGeoSys  for use within the mechanical solution. Wetting-induced collapse is modelled through an applied load activated when a threshold defined in terms of axial stress and net suction is exceeded. Water retention behaviour in unsaturated expansive soils has been shown to be both stress and volume dependent (Gallipoli et al. 2003;Dueck 2004;Dueck and Börgesson 2007). A simple extension to the method proposed by Dueck and Börgesson (2007) has been implemented in this model to complete the hydro-mechanical coupling.
The model described in this paper has been developed as part of the DECOVALEX-2015 project (Bond et al. 2014a(Bond et al. , 2015a and is used to simulate laboratory experiments undertaken by IRSN before application to more complex hydro-mechanical conditions posed by the presence of a technological void (a designed annular gap around the sample) and loss of confinement. Finally, a fullscale in situ performance test at Tournemire underground rock laboratory with a non-uniform technological void is simulated. We show that the general features of the experimental results, such as final swelling stress and water uptake, can be matched with this simple model with nine mechanical parameters but, in order to capture the full range of behaviour, more complex modelling methods incorporating elasto-plasticity may be required.

Mathematical model
In this work, we couple Richards' equation to a nonlinear elastic mechanical model incorporating a moving finite element mesh, in order to model the complex hydro-mechanical behaviour of bentonite with a limited parameter set.

The hydraulic process
The SEALEX experiments used in Task A of DECOVA-LEX-2015 are primarily concerned with the saturation of the bentonite seal from a compacted, partially saturated state. As a result, the description of fluid flow within the medium requires the consideration of more than one fluid phase. Multiphase flow in porous media is commonly described by the mass balance equations for each phase and Darcy's law of pressure driven fluid flux. A common simplification used in soil mechanics and hydro-geology for systems where the two fluids are a liquid and a gas is that of Richards' equation (Celia et al. 1990). The assumption made in Richards' equation is that the pressure change of the gas phase is negligible and therefore can be considered constant, i.e. the gas phase is infinitely mobile. This results in a direct relationship between the capillary pressure or total suction, and the degree of saturation. A pressure-based form of Richards' equation is given for the single fluid phase as derived by Thorenz (2001): where S w is the degree of saturation of the medium (-), S t is the storativity of the soil (-), n is the porosity of the medium (-), s is the suction (Pa) defined as the difference between the air pressure and the liquid pressure and therefore negative in unsaturated conditions, t is time (s), q w is the fluid density (kg/m 3 ), k rel is the relative permeability (-), k is the intrinsic permeability (m 2 ), l w is the dynamic viscosity of the wetting fluid (Pa s), g is the acceleration due to gravity in vector form of (0, 0, 1) to ensure it acts only in the z axis (m/s 2 ), and Q w represents a source term of fluid into the system (m 3 /s). This form of Richards' equation in OpenGeoSys is solved for the primary variable of fluid pressure, i.e. suction, with the degree of saturation and fluid velocity calculated as secondary variables. Two relationships are required to satisfy this equation, the constitutive relationship between the suction and degree of saturation, known as the water retention curve, and the dependence of permeability on the saturation of the sample, known as the relative permeability.
The accuracy of the coupled solution has been shown to rely heavily on the water retention curve, so it is important to choose an appropriate constitutive relationship for the material and anticipated stress conditions (Fredlund 2002;Fredlund et al. 2012). Commonly, a direct relationship between the degree of saturation and suction derived by constant volume water retention tests is fitted with water retention curve functions such as the van Genuchten relationship (van Genuchten 1980). However, these were developed for non-expansive soils and their applicability to expansive soils has been questioned due to the volumetric behaviour of expansive clays in free swelling stress conditions (Gallipoli et al. 2003). Specifically, it has been observed that variations in sample volume can lead to differences in void ratio that affect the pore dimensions and connectivity. Consequently, application of the van Genuchten formula to expansive soils requires the water retention properties to be cast in terms of degree of saturation, suction, and specific volume (Gallipoli et al. 2003).
Alternatively, the water retention curve can also be formulated in terms of degree of saturation, suction (s), and stress by directly incorporating the stress into the degree of saturation-suction relationship. This forms a state surface between suction, degree of saturation, and stress (Gallipoli et al. 2003;Dueck 2004;Dueck and Börgesson 2007). The free swell water retention experimental data provide the maximum water content at zero confining stress over the full range of applied suctions and associated changes in void ratio. Under constant volume conditions, a finite volume of water can be taken up by the sample over the full range of applied suctions, but a swelling pressure develops acting on the confining apparatus. This swelling pressure, also known as the confining pressure, is the modified effective stress (''r''): where r is the mechanical stress, u a is the air pressure (assumed to be zero in Richards' equation), b is the internal strain factor (discussed in the description of the hydromechanical coupling), Ds rel is the change in relative suction (as discussed below), and S w is the degree of saturation. The concept of Dueck suction states that the confining stress (p in Fig. 1) is equal to the suction difference at a given water content (w) between the free swell water retention curve (retention curve in Fig. 1) and constant volume water retention curve (confined retention curve in Fig. 1). This is shown in Fig. 1 where the water content at points F and G is equal on the two water retention curves, but corresponds to different values of suction (s F and s G ). This implies that the free swell water retention behaviour represents a limit to the maximum potential water uptake of the sample, and the constant volume water retention is a finite water uptake relative to the maximum potential. The suction that drives flow is therefore the relative suction, which is a stress-dependent parameter defined as: where s rel is the relative suction (Pa) that is available to drive flow, s fs is the free swell suction (Pa). This allows the determination of the hydration state of the sample under a range of confining stress conditions from constant volume to free swell conditions (Dueck and Börgesson 2007). Furthermore, it provides a direct coupling between the hydraulic and mechanical processes and provides a simple method of incorporating stress-dependent water retention behaviour into the model. Here we adopt the Dueck suction method and include a minor adaptation for consolidation stress conditions to complete the volumetric dependency of the water retention curve. Experimental data in the form of water content versus suction are converted to a degree of saturation for use in the pressure-based Richards' equation (Eq. 1). Water content by mass (w) is converted to a volumetric water content (h) defined as: by multiplication with the specific gravity of the sample: where SG is the specific gravity of the medium defined as: where q ref is the reference density of water = 1000 kg/m 3 and q d is the dry density of the material (kg/m 3 ). Dry density is thought to be the main control over the magnitude of swelling pressure and has therefore become a design criterion for bentonite applications as compacted blocks in a buffer or as pellets in backfill (Gens et al. 2011). Degree of saturation is calculated from: where n is the porosity. It is possible to define a theoretical maximum water content a specific element can contain at saturation by substituting Eq. (6) into Eq. (7) and rearranging to give: Assuming that saturation is achieved when all voids are taken up by fluid, degree of saturation is given by: The evolution of porosity during volumetric changes in the bentonite sample leads to a dependency of the degree of saturation on volume. The state surface defined by Dueck suction is applicable for stress conditions permitting constant volume and free swell behaviour (Dueck and Börgesson 2007). However, it has been noted that this approach may lead to unrealistic representations of saturation during loading-unloading experiments on clays, because irreversible void ratio changes are not taken into account by a unique water retention curve, i.e. the water content is equal for a given stress condition irrespective of whether this is achieved before or after loading (Gallipoli et al. 2003).
In this current work, the use of the free swell curve as the unique water retention curve coupled to the confining stress conditions is chosen to provide the appropriate volume-dependent water retention properties for constant volume and expansive conditions (Fig. 2,left). The retention curves shown in Fig. 1 are plotted on the state surface-the retention curve at zero confining stress (labelled, black) and the constant volume retention curve (red). In this figure, the saturation line in the stress-water content plane represents the maximum water content in the many possible degrees of confinement between constant volume conditions and free swell conditions. However, under stress conditions that lead to consolidation behaviour, i.e. external loading, the saturation line cannot remain static. In order to incorporate void ratio change in consolidation and its influence on the water retention properties, the degree of saturation limit is used to calculate the new location of the water retention curve for a given void ratio based on the theoretical local maximum water content. Consequently, the water content at a given suction is dependent on the stress state and the void ratio (or dry density), such that with decreasing void ratio the maximum water content reduces (Fig. 2). The location of the saturation line in the left-hand image of Fig. 2 is moved to a lower water content (i.e. moved towards the confining stress axis) if stress conditions leading to volume reduction prevail. This takes into account changing hydraulic properties with changing void ratio in compression and avoids a reversible water content behaviour for a loading-unloading stress path. Hysteresis of water retention properties, as experimentally shown by Ferrari and Seiphoori (2015), is not considered in this simple formulation of the stress state and volume-dependent water retention properties.
Data points are often limited in the water retention curve results at low suctions due to the large change in water content over a small range of suction. Therefore, to provide a smoothed relationship between suction pressure and water content, an empirical function proposed by (Bond and Benbow 2009) is used to match the free swell water retention data where M b is the mass fraction of bentonite (-), w is the water content by mass (-), and a, b, c, and d are all fitting constants. The free swell suction given by this equation is in MPa. The second relationship required to fulfil Richards' equation is the relative permeability (k rel ). A simple power law is used in this work: where S w is the degree of saturation and c is the power to which the degree of saturation is raised. The power is used to calibrate the relative permeability model to the infiltration test data.
In addition to the relative permeability, in stress conditions that give rise to volumetric changes the intrinsic permeability is likely to change, especially in active clays, such as bentonite (Karnland et al. 2006). Common relationships used to determine the saturated permeability for porous media relate the permeability to porosity, e.g. the Kozeny-Carmen relationship. However, it has been shown that the K-C relationship is most applicable to non-expansive materials and is less applicable to bentonites (Liu 2010;Liu et al. 2011). There are a variety of permeability models proposed for bentonites including the K-C relationship, cluster model (Achari et al. 1999), gel model (Pusch and Yong 2003), macroscopic Brookes and Corey function, statistical approaches (Agus 2005), and empirical formulae (Karnland et al. 2006). Liu et al. (2011) propose an edited K-C equation which aims to take into account the shape of the montmorillonite sheets and the montmorillonite content with respect to the solid grains and has been successfully used to calculate the saturated permeability of pure MX-80 bentonite: where d p is the particle thickness, C k is considered as the scaled pore shape factor, # p is the fraction of montmorillonite within the solid grains, u is the volume fraction of the solid, and m and l are calibrating model parameters, as in the Kozeny-Carmen relationship. The volume fraction of the solids is related to the porosity (n) by: This model can be used to simulate bentonites of varying montmorillonite content, including bentonite-sand mixtures, through the variation of the parameter # p in line with the montmorillonite fraction of the material solids. This has the effect of increasing or decreasing the permeability to reflect the incorporation of accessory minerals with larger grain sizes. d p 2 and C k are estimated as the average values from Liu et al. (2011), and m and l are calibrated to fit the infiltration test. This formulation introduces the dependency of the permeability on the porosity of the sample and therefore calculates permeability for a range of dry densities. In addition, it provides a second mechanical-hydraulic coupling mechanism through the volumetric change associated with a change in intrinsic permeability.

The mechanical process
Deformations of a medium can be described in terms of stress by the momentum balance equation: where r is the Cauchy stress tensor, q is the solid density (kg/m 3 ), g is a tensor with the acceleration due to gravity in the leading diagonal associated with the vertical direction (m/s 2 ). In this model, the primary variable of displacement is solved for by substituting the constitutive stress-strain relationship for a linear elastic material, i.e. Hooke's Law: where r 00 is here the modified effective stress tensor because suction is included in the mechanical solution as a body force, e is the strain, and D is the elastic constitutive matrix that defines the elastic material properties of the sample (Lewis and Schrefler 1998): E is the nonlinear Young's modulus (Pa), and m is the Poisson's ratio (-). Strain is calculated from the displacement by multiplication with the differential operator (L): e ¼ Lu: A global system of equation is built comprising the material properties in the stiffness matrix ([K]), the unknown displacements ({u}), and the known body forces or surface tractions in the force vector {f}. To solve for displacements, the stiffness matrix is inverted and multiplied with the force vector: A nonlinear elastic mechanical model based on the strain history is developed to describe the mechanical behaviour of bentonite under varying stress conditions, e.g. constant volume, free swell, and consolidation. An explicit assumption for the calculation of the material properties is used, such that Hooke's Law is calculated for each timestep, but non-linearity is achieved by updating the Young's modulus as a function of the strain on the sample. The Young's modulus is updated at each time-step and is given by the following strain-dependent empirical formula: where E t=0 is the reference initial Young's modulus (Pa), a is the hardening or softening factor depending on the strain direction, e v is the total volumetric strain, and d is a calibration parameter. This function leads to an increase in elastic modulus with a decrease in sample volume (i.e. strain hardening) and also a decrease in elastic modulus with an increase in sample volume (i.e. strain softening). The initial void ratio of the sample is used as a control to switch between the softening and hardening factors in the calculation of the Young's modulus. This allows softening and hardening parameters that reflect the different stress states and strain directions to be defined separately. As a result, four different conditions are considered: • volume reduction at void ratios greater than the initial void ratio, • volume expansion at void ratios greater than the initial void ratio, • volume reduction for void ratios less than the initial value, and • volume expansion for void ratios less than the initial value.
For the first two cases, the hardening and softening parameters are equal to the user-defined softening parameter. For the third case, the hardening parameter takes the user-defined value, and for the final case linear elastic rebound is assumed by calculating a constant Young's modulus.
Analysis of this empirical formula shows that calculations of large expansive strains could lead to an elastic modulus B 0.0. This is because the formula was initially developed in line with the experimental consolidation tests. Therefore, it requires extension to account for stress conditions that permit large expansive strains in order to ensure that the elastic modulus is always [ 0.0. This extension only applies in expansion and when the calculated Young's modulus is below a defined threshold. The Young's modulus used in the simulation step is then: where E min is a user-defined Young's modulus threshold (Pa), below which the effect of volumetric strain on reducing Young's modulus is reduced, and k is a negative exponent that controls the rate of reduction in Young's modulus with continued volumetric strain. As water content is closely linked to expansive volumetric strains in free swell conditions, this function implies that at higher suctions a small change in water content has a larger effect on the mechanical material properties than at low suctions. This aims to reflect the transition from dry, stiff conditions towards a gel so that at low suctions and large volumetric strains further expansion occurs with a less pronounced change in mechanical material parameters. The mechanical formulation is solved to give nodal displacements, stresses, and strains using the Galerkin finite element method in the open source code Open-GeoSys . In this current work, the nodal displacements are then used to update the mesh coordinates to produce a moving mesh known as an Updated Lagrangian mesh. The new element volume is used to calculate the volumetric strain, which is then used to calculate the new element void ratio, porosity, permeability, and dry density, as well as the new Young's modulus.
Void ratio (g) is determined from the incremental volumetric strain as: where g 0 is the void ratio of the element before the mesh coordinates were updated and De v is the increment of volumetric strain of the element. The void ratio is used to determine the element porosity (n): A new element dry density is given by: where q d is the dry density in kg/m 3 , M s e is the mass of solids of the element (kg), and V e is the volume of the element (m 3 ).
This empirical formula contains six parameters: the initial Young's modulus, the hardening and softening parameters, the calibration parameter (d), the minimum Young's modulus threshold, and the exponent k. In order to fully describe the elastic behaviour, the Poisson's ratio is also required. Other necessary inputs, such as porosity and initial void ratio, are determined from the experimental data.

Coupling
The hydraulic and mechanical processes governing the saturation of bentonite are extremely strongly coupled; a change in hydration state contributes to a change in the effective stress calculated by the mechanical process, and the mechanical stresses exert an influence over the amount of change of hydration that can occur (Laloui and Nuth 2005;Gens et al. 2006;Nuth and Laloui 2008). As a result, this coupling must be taken into account in both hydromechanic (HM) and mechano-hydraulic (MH) directions.
The coupling of the mechanical process to the hydraulic process (MH) is achieved in two parts: an evaluation of the suction-water content relationship as a function of the confining stress at a given void ratio and an evaluation of suction at the start of hydraulic iteration based on the confining stress conditions. Both parts of the coupling are required for the implementation of the Dueck suction concept (Eq. 3) to allow for a change in mechanical stress conditions to directly affect the hydraulic potential of the sample with an explicit calculation method.
The coupling of the hydraulic process to the mechanical process (HM) is incorporated through the calculated displacements being a function of both the hydraulic stress caused by the change in hydration state of the sample (change in suction) and the externally applied mechanical stress. Consequently, the effective stress can be expressed as: where the subscripts M and H refer to the mechanical and hydration stresses, respectively, and r 0 is the effective stress. The hydration stress is a function of the suction change and is determined from the change in water content calculated by the hydraulic solution. In the mechanical solution, the suction is treated as a body force as it simultaneously acts to hold the sample together and produce a hydraulic gradient that drives fluid flow. Therefore, a change in suction leads to a change in the stress state of the sample, e.g. an increase in hydration reduces the suction pressure holding the sample together and leads to a volumetric expansion. However, not all of suction change contributes to effective stress and subsequent deformation (Likos and Lu 2006). At low water contents, a reduction in suction causes swelling of clay clusters that is initially accommodated by crystalline swelling of aggregates into the macro-void spaces (Likos and Lu 2006; Wayllace 2008). As a result, only a small portion of the suction change contributes to effective stress. However, at higher water contents the bentonite aggregates have expanded to fill the macro-void spaces (under constant volume conditions) and the majority of further suction change contributes to effective stress, i.e. there is no space for internal deformation, and therefore, suction change translates directly to stress build-up. A similar process occurs if the sample is unconstrained, but the suction change contributes to swelling (Likos and Lu 2006). At low water contents, bentonite aggregates swell into macro-voids and the bulk volume change for a given suction change is less than at high water contents where large volume change occurs for a small suction change (Likos and Lu 2006). This is in part because the water adsorption mechanism transitions from crystalline swelling, which is dominant from high suction to around 20 MPa, to osmotic swelling which is dominant from around 0.8 MPa suction (Wayllace 2008). Typically for expansive clays, it is assumed that the intra-aggregate voids are saturated, i.e. crystalline swelling is complete, before osmotic swelling predominates at water contents of [40% (Madsen and Müller-Vonmoos 1989;Wayllace 2008).
In order to represent both the accommodation of aggregate swelling by the macro-porosity and the transition from crystalline to osmotic swelling processes, a factor, termed here the internal strain factor (b), is applied to determine the contribution of hydraulic stress to effective stress. The contribution is also a function of the degree of saturation, such that hydraulic stress is given by: where Ds rel is the change in suction within the sample (a negative value in hydration), b is the internal strain factor, and S w is the degree of saturation. Substituting Eq. (24) into Eq. (23) and using the notation of Bishop's effective stress, this formulation can be shown to be similar to the use of Biot's constant and gives the modified effective stress (as in Eq. 2): where r 00 is the modified effective stress, (ru a ) is the net stress, and v is the Bishop's parameter that is here assumed to be effectively equivalent to the degree of saturation (S w ) (Lewis and Schrefler 1998).
In unsaturated non-expansive materials, Biot's constant is a proportionality factor between the suction of the fluid in the pores caused by the capillarity effect and the effective stress that is responsible for volumetric deformation (Lewis and Schrefler 1998). However, Agus et al. (2010) indicate that the dry density of bentonite, and therefore pore space size and shape, has a minimal contribution to the total suction. They conclude that suction in an expansive soil is dominated by the physico-chemical hydration forces within a clay aggregate. Following the interpretation that hydrating fluid is taken up into the clay aggregates in preference to filling the macro-porosity (Kröhn 2003), the influence of pore fluid pressure is minimal. Instead, the influence of the intra-aggregate suction is dominant and the swelling processes of aggregates define the contribution of hydraulic stress to the total effective stress.
Instead of representing an effect of pore fluid pressure, the internal strain factor therefore represents the coupled hydro-mechanical effect of the double structure of porosity observed in bentonites and the swelling mechanism operating at different water contents (Madsen and Müller-Vonmoos 1989;Lloret et al. 2003;Agus 2005;Pusch and Yong 2006;Wayllace 2008;Navarro et al. 2015). Furthermore, this process is dry density dependent because the volume of macro-void space that can accommodate expanding clay clusters is reduced for a highly compacted sample in comparison with a low dry density sample (Fig. 3) Agus 2005;Agus et al. 2010;Wang et al. 2012Wang et al. , 2014Wang 2012a;Saba et al. 2014a).
Due to the non-linearity in the relationship between suction and mechanical stress, the coupling procedure is solved iteratively via a Picard iteration between the processes. The staggered, explicit iterative procedure employed is shown in Fig. 4, where s rel is the net suction (Pa), s fs is the free swell suction (Pa), w is the gravimetric water content (%), r 00 is the modified effective stress (Pa), (ru a ) is the external confining mechanical stress (Pa), E is the Young's modulus (Pa), n is the porosity (%), g is the void ratio, q d is the dry density of the material (kg/m 3 ), and superscript i refers to the iteration number.
The suction relative to the maximum potential suction at a given water content is used from the previous timestep (or initial conditions) as input to the hydraulic model as calculated from Dueck suction (Eq. 2). The hydraulic process calculates a new fluid pressure and water content via the saturation-dependent relative permeability function. The resulting change in hydration state is used to calculate any displacements and a new modified effective stress in the mechanical solution, taking into account any external forces (ru a ). If the system is perturbed by a mechanical stress, this leads to change in the hydraulic properties in the pre-iteration calculation of the relative suction via the Dueck suction concept, i.e. a change in mechanical conditions causes a change to the initial hydraulic conditions for the next iteration. As a result of the dependence of the hydraulic properties on the modified effective stress (Eq. 2), an iterative procedure is required between these processes such that the modified effective stress between two iterations is within a predefined tolerance: Convergence Calculated displacements lead to a change in the volume of a sample and a change in porosity, dry density, and permeability. Therefore, the hydraulic solution needs updating with each iteration to ensure the consistency of the hydromechanical coupling. To achieve this, we calculate the hydraulic process over an intra-time-step Updated Lagrangian mesh defined with each iteration (Fig. 4). Once the iterations have converged, the calculated displacements are used to update the nodal coordinates and the volumedependent parameters for both processes for the start of the next time-step. In the following sections, this model is applied to a series of experiments, known as the SEALEX experiments, conducted by the Institut de Radioprotection et de Sûreté Nucléaire (IRSN) at the Tournemire underground research laboratory.

SEALEX experiments and model results
The SEALEX experimental programme comprises three main stages: a series of laboratory tests (that form modelling Steps 0 and 1), a hydro-mechanical performance test of the host rock (that forms modelling Step 2), and a fullscale in situ performance test of a bentonite core in the host rock (that forms modelling Step 3). This work is primarily concerned with the development of a hydro-mechanical bentonite model and therefore presents results for Step 0, Step 1, and Step 3.
The SEALEX experiments, as described in Barnichon and Deleruyelle (2010), Barnichon et al. (2012) and Wang (2012), investigate the hydro-mechanical behaviour of a potential hydraulic seal. A 70:30 MX-80 bentonite-sand mixture compacted to dry densities of 1.67 and 1.97 Mg/ m 3 with an approximate water content of 11% by mass were used in both the laboratory tests and field performance test. MX-80 bentonite is a sodium-dominated bentonite that has a montmorillonite component of 92% in this experimental programme (Wang 2012).
The laboratory tests consist of two main stages: standard laboratory tests that make up the results for the Step 0 modelling task and a 1:10 scale mock-up test of the in situ water injection test that makes up Step 1. The Step 0 tests were designed to be used to parameterise the chosen or developed model and comprised confined and unconfined water retention tests, constant volume infiltration, and constant suction oedometer consolidation tests (Barnichon et al. 2012;Wang 2012).
For the consolidation tests, three samples were compacted to a dry density of 1.67 Mg/m 3 , emplaced in the oedometer with no annular void, and hydrated to an imposed suction of 38, 12.6, and 4.2 MPa, respectively, before consolidation was performed. The fourth consolidation test was performed on a bentonite made up of the same composition compacted to 1.97 Mg/m 3 dry density and hydrated to a suction of 0 MPa (Wang 2012). For Step 1, the experimental set-up comprised a 0.12-mlong bentonite sample compacted to 1.97 Mg/m 3 dry density with a diameter of 0.0555 m emplaced in a hydration cell with a diameter of 0.06 m. This produced a symmetric technological void which was initially flooded with water. The 1:10 scale mock-up test was then divided into three main phases, with Phase 2 further subdivided into 2a and 2b: • Phase 1: Initial saturation phase: flooding of the technological void and subsequent hydration from both the void and the base of the sample. During this phase, axial deformation was restrained and the build-up of axial stress was measured. The injected water volume was also measured. • Phase 2a: The void recovery phase was initiated after axial stress equilibrium and involved the release of the vertical constraint with hydration from the base only. • Phase 2b: Additional hydration from the top surface after 2.8% axial strain had accumulated. • Phase 3: The confinement phase began after 20% axial strain when constant height conditions were imposed again and the stress build-up observed.
The third modelling step was the full-scale in situ performance test PT-A1. This consisted of a 1.2-m-long compacted bentonite-sand core with a diameter of 0.56 m emplaced in a horizontal drift hole with a diameter of 0.6 m leading to a non-symmetric technological void. Relative humidity was recorded at specific intervals along the core, axial stress at each end, and radial stress half way along the core at evenly spaced intervals around the circumference of the core (Barnichon et al. 2012;Wang 2012). In the following sections, we present the results of the model from each of the three steps.
Step 0 The Step 0 modelling step consists of the water retention test, the constant volume infiltration test, and consolidation tests at imposed suctions.

Water retention curve
The water retention curve is calibrated to the free swell water retention test results using the empirical formula defined in Eq. (9). The HM coupling scheme based on the Dueck suction method leads to a confined water retention curve calculated by the model as a function of the confining stress (and void ratio during consolidation). Two experiments under free swell and constant volume conditions were conducted on pre-compacted bentonite at a dry density of 1.67 Mg/m 3 with an initial water content of 11% and initial suction of 65 MPa calculated from the relative humidity. Relative humidity is related to the suction (s) through Kelvin's law: where R is the gas constant (8.314 J/mol/K), T is the absolute temperature (K), M w is the molecular weight of water (18.016 kg/kmol), q w is the unit weight of water in kg/m 3 , and RH is the relative humidity (Blatz et al. 2008). Initial relative humidity is 62 ± 1% which corresponds to 64.5 ± 1.5% MPa suction. Figure 5 shows the calibration This calibrated curve was then used throughout the simulations with the Dueck suction concept to determine the water content at any given stress state and suction condition. Figure 5 also shows the confined water retention data calculated by the model for a constant volume simulation (blue) as a demonstration of the applicability of the approach.

Infiltration test
The infiltration test in Step 0 was performed in a metallic cylinder 250 mm high with a diameter of 50 mm and capped to maintain constant volume conditions. A precompacted bentonite sample with dry density of 1.67 Mg/ m 3 and a water content of 11% was placed in the cylinder and hydrated from the base. Relative humidity was measured at four sensors through the sample: 50, 100, 150, and 200 mm.
A 2D axisymmetric fully coupled HM model was calibrated to the experimental data. The results are shown in Fig. 6. The (Liu et al. 2011) permeability model and relative permeability are calibrated as shown in Table 1. The coupled hydro-mechanical model is able to reproduce the general trends of the simple hydration test at the sensors from 100 to 200 mm from the water inlet source in the SEALEX experiments (left of Fig. 6). Final hydration states at all sensors are calculated to within 2% relative humidity. However, the rapid hydration observed at the lowest sensor (red) is under predicted. This is a common problem throughout the modelling teams in the DECO-VALEX-2015 task (see summary paper this issue, Millard et al. 2016) and could be due to the permeability model not capturing the full range of permeability change associated with the hydration or the use of a single unique water retention curve to describe the water content-suctionstress relationship rather than specifically addressing the different scales of porosity.
The applicability of this model is shown in Fig. 6b by HM modelling of the infiltration test performed on the Gao-Miao-Zi (GMZ) bentonite by Ye et al. (2009). The GMZ bentonite is compacted to a different dry density, is a pure bentonite (i.e. M b = 1.0), and has a significantly different chemical make-up [less Na ? and more Ca 2? and Mg 2? (Ye et al. 2009)]. Nevertheless, the permeability model is able to give a good fit to the experimental data with only a very minor change to the parameterisation to account for all these differences (Table 1). The results indicate that this model can be applied to model bentonite and provides a good calibration to the SEALEX Step 0 data as a good starting point for the later modelling steps.

Consolidation tests
The consolidation tests in Step 0 were performed at different imposed suctions of 0, 4.2, 12.6, and 38 MPa. The 0 MPa suction case was performed on a bentonite compacted to a dry density of 1.97 Mg/m 3 and trimmed to diameter of 35.13 mm to leave a technological void around the edge of the sample. The higher suction tests were conducted on samples prepared in the same manner as the Fig.  infiltration test and had a diameter of 38 mm equal to that of the oedometer ring. There was no technological void for these samples. All the samples were hydrated from the base (and the technological void in the 0 MPa suction case) and allowed to swell until the imposed suction. On stabilisation of the axial strain, the samples were subjected to an axial loading-unloading cycle with a maximum net axial stress of 49 MPa. A 2D axisymmetric fully coupled HM model using the hydraulic parameters above was calibrated to the three consolidation tests conducted on the 1.67 Mg/m 3 dry density samples. The models were run from a consistent starting material and follow the experimental procedure of swelling to an imposed suction, followed by consolidation and then unloading. The calibration presented is consistent with the mechanical model used for the later modelling steps (Step 1 and Step 3) and as such does not necessarily represent the best possible fit to the individual consolidation test results. The internal strain factor for the 1.67 Mg/m 3 dry density samples is given in (Fig. 7), and the full mechanical calibration is given in Table 2.
The internal strain factor is low for low water contents where crystalline swelling is dominant and macro-pores accommodate volume expansion of aggregates (Likos and Lu 2006;Wayllace 2008). Osmotic swelling is thought to predominate at water contents [40% once crystalline swelling is complete (Madsen and Müller-Vonmoos 1989;Wayllace 2008) and as such macro-pores are filled and   Nguyen et al. (2005). The nonlinear elastic approach is capable of reproducing the irrecoverable strains observed in the oedometer tests (Fig. 8). The model is able to adequately reproduce the behaviour of the lower suction samples, but the model calibration is unable to recreate the increasing stiffness of the higher suction sample (38 MPa). The experimental data indicate that the pre-consolidation pressure increases with increasing suction in a linear fashion (Wang 2012) and is thought to reflect the movement of the failure curve with respect to suction in elasto-plastic models, i.e. suction hardening (Gens and Alonso 1992). The nonlinear elastic model is not able to recreate this with a calibration consistent with the later modelling steps. However, this model requires just 8 parameters to be able to account for irreversible strains and reproduce the major trends of the data. As the aim of this study is to investigate an alternative model of bentonite with an emphasis on keeping the number of parameters to a minimum, the results obtained by this nonlinear elastic model are considered sufficient to move forward to the application of the model to more complicated hydro-mechanical conditions.
Step 1 The laboratory experiment used for the Step 1 modelling task was a 1:10 scale mock-up of the SEALEX in situ test. The experimental set-up is shown in Fig. 9 and comprises a 55 mm diameter by 120 mm height compacted bentonitesand mixture sample (with a ratio of 70:30) at a dry density of 1.97 Mg/m 3 . This was placed in a rigid hydration cell with a diameter of 60 mm to simulate a technological void around the bentonite sample. The technological void is symmetric in this mock-up test which differs slightly to the full-scale in situ test. The top of the cell is a piston 60 mm in diameter and 150 mm long that enables control over the axial swelling.

Model set-up
A challenge of this modelling stage was how to represent the presence and effect of the technological void in the simulation. Here we use the moving finite element mesh to keep track of the bentonite swelling into the void space and when the void has been filled by the expanding bentonite. Displacement-dependent boundary conditions for both the hydraulic and mechanical processes are implemented. Once the maximum displacement of a boundary node has been reached, the hydraulic boundary condition of the technological void representing liquid hydration is removed and a zero displacement mechanical boundary introduced. These boundary conditions can be expressed as: where BC H and BC M are the hydraulic and mechanical boundary conditions, respectively, s is the suction (Pa), Q is the flux (m 3 /s), r c is the confining stress (Pa), u is the calculated displacement (m), and u max is the displacement limit (m) that represents the technological void in Phase 1 and the axial displacement in Phase 2b. This representation of the technological void defines the void to be filled with liquid water until the bentonite expansion fills the void space. It therefore does not take into account the formation of colloidal suspensions or gel that has been observed in experiment (Saba et al. 2014b). A 2D axisymmetric mesh, shown in Fig. 10, is used with the fully coupled HM model and the Step 0 calibration presented earlier. However, the Step 1 experiments are conducted on samples with a higher dry density than the calibrated model. The difference in dry density is reflected in the model through the internal strain factor. The higher dry density material has less macro-voids, and as a result clay aggregates have less space to swell into as hydration proceeds (Likos and Lu 2006). This results in a greater transfer of suction to effective stress and requires an increase in the internal strain factor (Fig. 11). This function coupled with the Dueck suction concept has the effect of ensuring that a higher dry density sample takes up less water than a lower dry density sample under constant volume conditions for a given change in suction and a higher swelling pressure for a given water content. The Step 1 axial stress results for Phase 1 exhibit a wetting-induced collapse phenomenon that causes a reduction in axial stress as hydration proceeds. This has been observed in other experimental programmes Romero et al. 2003;Gens et al. 2011) during hydration under constant volume conditions, although  Fig. 11 Comparison of internal strain factor used for the different dry density materials more commonly in Ca 2? bentonites such as Calcigel . The process has been hypothesised in a number of ways regarding the collapse of the macrostructure, e.g. as a result of cementing bonds between the grains in initially dry samples being broken by the hydrating water (Agus 2005); by collapse of the macrostructure in response to the developed swelling pressure/confining load ; or by structural rearrangement of clay clusters within the sample as a result of the frictional coefficient of the soil being overcome by wetting-induced shear forces at the micro-scale (Agus 2005). These mechanisms all consider the wetting-induced plastic collapse phenomenon to be as a result of the interaction between micro-structural processes and the collapse of the macro-structure.
Under experimental constant volume conditions, the collapse or rearrangement of the macro-structure leads to a reduction in axial stress, and in order to maintain the experimental constant volume conditions the applied load is reduced . Collapse of the macrostructure is associated with macro-structure compressive strains which are partially offset by expansion of the microstructure to fill the voids Gatabin et al. 2016). As a result, the stress reduces until the microstructure strain leads to a further build-up in stress. This is interpreted as a plastic phenomenon, and therefore, due to the elastic assumption in this work, the simple formulation presented is not able to predict this behaviour. In order to represent these internal compressive strains and reproduce the expansive collapse phenomenon, a source term model to account for the reduction in axial stress is proposed.
The source term is a mechanical stress applied to the model and is activated by a failure curve that can be derived from experimental results. Wang (2012) identified a linear relationship between the yield stress and net suction in the consolidation tests (Fig. 12). Considering the yield stress to be the axial stress at which plastic behaviour begins for the imposed net suction, the linear trend is used to define an axial stress-suction relationship that determines the activation of the source term. This is similar to the loading-collapse curve in elasto-plastic models that represents the change in yield curve location with respect to net mean stress and suction (Alonso et al. 1990;Gens and Alonso 1992;Alonso and Vaunat 1999). The calculated source term acting in the opposing direction to the build-up of axial stress is then a function of the proportion of the sample with stress-suction conditions in excess of the yield curve (Fig. 13). The dependency of the onset of the source term on the axial stress and suction provides a method to calculate the gradual development of the expansive collapse phenomenon seen in experimental results.

Results of Step 1
Experimental axial stress measurements were taken from the 60-mm-diameter piston at the top of the sample (Wang 2012) and as such represent the mean axial stress across the top surface. The model results are therefore calculated from the top surface of the 2D axisymmetric model and presented along with water uptake in Fig. 14. Model results are presented for both the nonlinear elastic model (NLE) without accounting for wetting-induced collapse and the nonlinear elastic source term model (NLE-ST).
The axial stress results for the two models (NLE and NLE-ST) show that the nonlinear elastic model over-predicts the axial stress build-up during Phase 1 as the plastic collapse behaviour is not captured. With the source term addition to the model, the final stress can be well matched and the magnitude of energy lost in the process is quantified by the gap between the NLE and NLE-ST model results. The NLE-ST model provides a good fit to the experimental data, capturing a gradual collapse process and the subsequent build-up in stress. The water uptake behaviour is well matched during the early stages of the model, but the maximum water is overpredicted by the end of Phase 1. Due to the application of the boundary conditions as a mechanism to represent the technological void, the water uptake and swelling of the bentonite are calculated to occur over a short time period. The technological void is predicted to be closed after 32 h. This rate of expansion is slower than that observed during experiments (Saba et al. 2014b) on a 35-mm diameter 1.97 Mg/m 3 dry density sample with a 14% technological void, where radial stress builds up after less than 1 h. It does, however, indicate that the technological void is filled quickly before the process of homogenisation proceeds due to continued swelling as the dry centre is gradually hydrated (Saba et al. 2014b).
During Phase 2, displacement develops at two distinct rates relating to the hydration conditions: slower deformation occurs initially as hydration is from the base only (Phase 2a), but more rapid deformation proceeds when hydration is from both the base and the top (Phase 2b). Displacement was measured throughout, and the model calculation from the NLE-ST model is presented in Fig. 15. The end of Phase 2a was temporally defined in the model, but the end of Phase 2b/start of Phase 3 is a displacement-dependent boundary condition that becomes active after 20% axial strain and is therefore a predicted result.
The elastic rebound calculated by this model is overpredicted due to the excess of stress in the system and the strain-softened material properties associated with c. 14% swelling into the technological void. The source term model does not remove any of the stress from the system that can be released from the unloading of the piston, and as a result the initial elastic rebound is too large. The activation of the source term as a function of axial stress and suction leads to the removal of the source term during Phase 2a. The release of this energy leads to a large elastic rebound which also causes the 20% axial strain to be reached twice as quickly as the experimental results.
The stress build-up during Phase 3 is shown in Fig. 16. The model predicts an excellent fit to the final stress build-up during Phase 3. Furthermore, the evolution of the stress build-up is very similar to that of the experimental results. However, due to the over-predicted rebound in Phase 2 the onset of stress build-up is too early. Comparing the stress build-up without considering the temporal prediction of Phase 2 shows an excellent match to the experimental data (Fig. 17). The main control over the development of this axial stress in Phase 3 is the suction remaining within the sample that can be alleviated under constant volume conditions. Although the NLE-ST model is not able to reproduce all the experimental data, it is able to represent all the major steps of the complex hydro-mechanical behaviour with just six mechanical parameters and an experimentally defined stress threshold. Final axial stress calculations and water uptake calculations match well with experimental data, but the temporal evolution of the deformation is found harder to predict without the consideration of plasticity.
Step 3 The third modelling step of the SEALEX experiments is of the full-scale in situ test PT-A1 (Barnichon et al. 2012). The PT-A1 performance test consisted of a bentonite-sand core 1.2 m in length and 0.56 m in diameter emplaced in a horizontal drift hole with a diameter of 0.60 m. Constant volume conditions were maintained by a confining plug in the gallery end and water inlets at both ends connected to the same water tank. Due to the horizontal emplacement, the technological void was non-uniform around the sample as shown in Fig. 18. Water was injected into the system under 0.2 MPa pressure to ensure the technological void was flooded. Hydration proceeded from the technological void and through the sample ends. The bentonite-sand core contained relative humidity sensors at 0.22 (RH22-1, RH22-2) and 0.52 m (RH52-1, RH52-2), axial stress sensors at both ends, and radial stress sensors evenly spaced around the edge of the sample at 0.6 m (Fig. 18).

Model set-up
The non-uniform nature of the technological void presents a challenge for a 2D axisymmetric approach so this is supplemented with a 2D slice model at 0.60 m. The Step 3 2D axisymmetric model is an upscaled version of the Step 1 mock-up test. The respective finite element meshes are presented in Fig. 19.
The displacement-dependent boundary conditions used for Step 1 are again employed in Step 3 and extended to the 2D slice model by defining the displacement limit from the equation for a circle with the same radius as the excavated drift hole. The boundary conditions for the 2D slice model can be expressed as: where x max is the calculated coordinate on the circumference of a circle with radius r for the nodal coordinate y n . This value is compared to the mesh x coordinate (x n ). The mechanical model employed is the nonlinear elastic model without the source term approximation of plasticity. The hydraulic model is consistent with the previous steps except for an increase in permeability as shown in Table 3.

Results of modelling the PT-A1 performance test
The axisymmetric model is used to calculate the axial stress at each end of the sample: 1.2 and 0.0 m and the relative humidity at sensors RH22-2 and for both sensors at RH52. An axisymmetric model is not capable of predicting the relative humidity of RH22-1 due to the non-uniform technological void in the experiment so the 2D slice is used to investigate the potential difference in relative humidity due to the non-symmetric technological void.
The axial stress results in Fig. 20 indicate that the model provides a good fit to the experimental data for the first 180 days. After this, the experimental results plateau, but the model calculates a continued increase. The experimental results do not show the development of stress until 20 days due to the lag in hydration during the experimental procedure (Wang 2012). The gradient of the stress development matches sensor AS_000, but does not predict the non-symmetry stress development between the two sample ends. However, the physical reason for the dissymmetry in axial stress is not known (Millard et al. 2016). For example, it could be an experimental difference in stress conditions at each end of the sample, such as an unaccounted for void into which the bentonite was able to expand, rather than a physical process occurring within the sample during hydration.
The model results for relative humidity at sensors RH52-2, RH22-1, and RH22-2 are presented in Fig. 21. Relative humidity is calculated to within 1% RH at RH52-2 with the elevated permeability model. The final magnitude of relative humidity is also closely matched at sensor RH52-2 with the imposed displacement-dependent boundary conditions suggesting a lack of fluid input once the technological void is filled. The 2D slice model is able to reproduce the difference in relative humidity between Fig. 18 Schematic of the experimental set-up for the PT-A1 performance test. The core contained multiple sensors whose location is shown in the schematic diagrams on the right. The non-uniform technological void is shown as a result of the horizontal experimental set-up  Fig. 21. This is because the sensors RH22-1 and RH22-2 are located closer to the source of the forced hydration than RH52. The axisymmetric model suggests the hydrating front from the base reaches 0.22 m along the sample within a time frame of 80 days, after which Fig. 21 shows the model results diverging from the data. The greater divergence of the results for RH22-1 compared to RH22-2 suggests that the non-symmetric technological void results in a more significant contribution of fluid from the end hydration at RH22-1. Model results in which hydration continued through the radial boundary were not able to recreate the final relative humidity or the dissymmetry between sensors RH22-1 and RH22-2. Therefore, the results indicate that fluid from the host rock does not contribute significantly to the hydration of the sample. These results highlight the relative contributions of hydration from the technological void and from the forced end hydration. Radial stress at sensor S60-2 is calculated with the 2D slice model and provides a good match to the rate of stress increase in comparison with experimental data (Fig. 22). However, the model calculates the technological void to be fully closed after 126 days and radial stress increases from this point. Similarly to the Step 1 model results, the model  calculates a closure of the technological void later than observed in the experiment (40 days). Although the model does not capture the gradual increase in radial stress, which is likely due to the formation of mechanically weak gel in the technological void as observed in small-scale experiments (Saba et al. 2014b), the later development of stress is well matched.

Conclusions
The aim of this study was to investigate the complex hydro-mechanical behaviour of bentonite with as simple a constitutive model and as few parameters as possible. In order to do so, a nonlinear elastic model was coupled to an unsaturated fluid flow simulation and used to model the hydro-mechanical behaviour of compacted MX-80 bentonite in the SEALEX experiments. An Updated Lagrangian mesh application was implemented in OpenGeoSys and used to inform the material non-linearity in the mechanical solution.
A series of parameterisation laboratory tests were reproduced by the simple HM-coupled model with an adequate fit achieved for consolidation at different suctions and constant volume infiltration. The coupling mechanism of Dueck suction using the free swell water retention curve and mean net stress to determine constant volume water retention properties was shown to be suitable for modelling the water retention behaviour at different void ratios and stresses. This model was then extended to simulate the complex behaviour of bentonite swelling into a technological void in both small-scale laboratory experiments and a full-scale in situ performance test.
The technological void was modelled using displacement-dependent boundary conditions for the hydro-mechanical conditions of the gap without considering the formation of a colloid or gel. The implemented Updated Lagrangian mesh was used to determine the closure of the technological void. Calculated water uptake from the technological void indicates that the void is filled within 2 days.
The simple model requires few parameters to describe the hydro-mechanical behaviour of the simple laboratory tests and the full-scale in situ performance test PT-A1. However, the wetting-induced collapse hydro-mechanical phenomenon seen in the 1:10 scale mock-up test cannot be modelled using only a nonlinear elastic formulation. Instead, an extension to the model incorporating a source term to account for plastic strains is presented. This enables the main features of the experimental results, such as axial stress development, water uptake, and relative humidity in the different experimental conditions to be reproduced with just nine mechanical parameters.
This empirical model is able to simulate the reduction in axial stress during wetting-induced collapse, but the elastic assumption in the constitutive framework results in an over prediction of the rebound displacement upon loss of confinement that could potentially be better predicted with more complex elasto-plastic models.
Model simulations of the PT-A1 performance test reproduced the axial stress evolution at one end of the sample and non-symmetric water up take due to the geometry of the technological void. Furthermore, relative humidity results indicate that hydration of the hydraulic seal from the host rock is negligible.