This paper describes an on-going experimental and numerical modelling research project on salt rock specimens. The experimental part of the study consists of a number of Mode I fracture tests with the WST (Wedge-Splitting Test) configuration, which are performed at different loading rates and complemented by a series of standard uniaxial creep tests. The preliminary WST results show a greater mechanical fracture work accompanied with lower force peaks, for the slower tests. As a first attempt to represent the experimental results, an in-house Finite Element model has been used, which combines an inviscid discrete fracture approach with a Maxwell chain model for the continuum material. The simulations show a decrease of the mechanical work needed for opening the fracture and higher peak force, as foreseen by the ongoing experimental results, but not with the same intensity, which seems to indicate that work dissipation may not be caused exclusively by the bulk viscosity.

In recent years, the authors and co-workers have developed a 3D finite element model for coupled thermo-hydro-mechanical (THM) problems in fractured rock masses. Zero-thickness interface elements are used for taking into account explicitly the effect of fractures and discontinuities in the fluid flow as well as the effect of fluid pressure in the crack propagation. Furthermore, the use of zero-thickness elements as a discrete modelling approach for fractures and discontinuities makes it possible to account for the heat transport taking place within these elements, even when advection dominates over diffusion (high Peclet number) [1]. The model has been implemented in the finite element code DRAC5, which is equipped with fracture-based interface elements and MPI parallel capabilities [2]. The code was originally developed considering water-saturated porous medium and fractures. The new developments described in the present paper, include the extension of the original formulation to the case of two-phase (liquid and gas) flow within the porous medium and discontinuities. The liquid includes only liquid water species, while the gas phase includes water vapour and gas species. The formulation includes the equilibrium equation, the mass balance of water and gas species and the energy balance equation. The parameters of the retention and relative permeability curves for the interface elements, such as the gas entry value and the residual water saturation, are updated with the variation of the normal aperture. The new capabilities of the model are illustrated with some academic verification examples.

In the framework of the Finite Element Method, zero-thickness interface elements have been widely used to model fracturing processes in quasi-brittle materials in a broad variety of problems. In particular, interface elements equipped with elastoplastic constitutive laws that account for the softening of the material strength parameters due to the fracturing mechanical work has been proved to accurately reproduce observed fracture propagation behaviour in concrete. Along this line, this paper presents the extension of an existing constitutive law of this kind to include the effect of chemical degradation of the material in the formation of fractures. The law is defined in terms of the normal and shear stresses on the average plane of the crack and the corresponding normal and shear relative displacements. A hyperbolic cracking (plastification) surface in the stress state determines the crack initiation. The softening of the cracking surface is governed by two history variables: an internal variable that accounts for the dissipated fracturing (plastic) work, and an external variable to be provided by a chemical degradation model that accounts for the effect of chemical degradation on the strength parameters. After a detailed discussion of the formulation, the main characteristics of the proposed law are illustrated with a number of academic examples for different combinations of mechanical loading and chemical degradation sequences. The model is finally validated against experimental results from the literature consisting of three-point bending tests performed on mortar samples previously exposed to an aggressive solution for different time periods.

In Carbon Capture Storage (CCS) sites, an important element of risk to be considered is the integrity of the cement seals of the abandoned wells in the reservoir [1]. The main goal of abandonment procedure once the life of a well is completed is to provide an effective isolation of the reservoir fluids in order to reduce environmental risk of contamination. In the case that the site has been reconverted to CCS, this is even more essential to prevent CO2 leaks from the storage site. It is important to note that the cracking conditions of the well cement seal can be affected by the long term changes in pore pressures that take place after the oil exploitation activities have stopped [2]. For example, slow pressure return around extraction wells (where the pore pressure had been subject to a sustained reduction during long extraction periods) may cause a progressive reduction of the effective stresses acting on the cement casing and plug, while the opposite can happen at injection wells. And these effects may be partially modified by the overall structural response due to the volume changes implied by the effective stress changes [3]. In this paper, a preliminary study of the effects of such stress changes on the potential integrity of a 2D cross-section of the sealed oil well system (caprock-external cement sheath-steel casingcement plug) during its service-life (injection/production activities and abandonment) has been performed by means of FE method including zero-thickness interface elements to represent potential cracks. In particular, these elements are pre-inserted in the analysis in between the contacts of caprock-external cement sheath, external cement sheath-casing and casing-cement plug. The results presented show that, depending on the initial state and range of pressure evolution, the different interfaces considered may open or close in a non-trivial manner during the pressure return process. This seems to indicate the importance of considering carefully the pressure return process and subsequent effective stresses evolution in abandoned reservoirs recycled to CCS, in order to avoid that new cracks in well cement seals may lead to potential CO2 leakage in the storage site.

In underground salt and potash mines, such as the ones in the region of Bages in central Catalonia, long-term deformations and stability of the mine tunnels are strongly influenced by the viscous (creep) behavior of saline rock. Tunnel excavation causes deviatoric stresses which in turn trigger creep strain in the salt rock. It has been often considered that the main consequences of creep are time-dependent convergences with potential consequences for the serviceability of the mine. However, creep deformations may be accompanied by stress redistribution, and potentially also by cracking and fracture, especially if rock exhibits preexisting discontinuities and/or layers separated by weaker contacts. In this way, tunnel crosssections that are perfectly stable after excavation, with time may not only deform, but also approach failure collapse. Traditionally, salt rock has been mainly characterized with regard to the creep behavior, while strength issues have been addressed independently in ways similar to other rock materials, via strength criteria such as Mohr-Coulomb, etc. However, there seem to be very few studies considering the fracture mechanics of salt rock. In this paper, the on-going experimental and numerical research along this line at ETSECCPBUPC is described. This includes standard creep tests as well as mode I fracture Wedge Splitting Tests (WST). WSTs have been performed at different loading rates with the purpose of assessing the time dependency of parameters such as the tensile strength and the specific apparent fracture energy. Preliminary test results show that as the loading rate is increased, the tensile strength seems higher but the apparent fracture decreases. Numerical calculations include finite element simulations of the WST, as it has been done before for other rock types by some of the authors1. Continuum elements with visco-elastic behavior were used to represent the salt rock material, while the fracture path was represented via preinserted zero-thickness interface elements. As a first attempt, the constitutive model used for the interface elements was a (time-independent) elasto-plasticity framework with fracture energybased evolution laws 2 . The preliminary numerical simulations qualitatively reproduce the experimental results, however, it seems that in order to quantitatively fit the experimental results, a new time-dependent constitutive model for the interface elements would be required in which fracture properties and evolution should interact with creep and time-dependent behavior.

Alkali-Silica Reaction (ASR) is a particular type of chemical reaction in concrete, which produces cracking and overall expansion of the affected structural element due to the formation of expansive reaction products within the cracks. This paper develops the formulation of a coupled Chemo-Mechanical (C-M) Finite Element (FE) model for simulating ASR expansions in soda-lime glass concrete at the meso-scale. The model considers several C-M coupling mechanisms, including a reaction-expansion mechanism qualitatively proposed by the authors elsewhere on the basis of experimental results, which is introduced in order to reproduce the effect of compressive stresses on the development of ASR expansions. The model has the characteristic ingredient of using zero-thickness interface elements for modelling the C-M mechanisms leading to the propagation of cracks due to the formation of ASR products within them. This fact has required the development of: (i) a new FE formulation for diffusion-reaction processes occurring within discontinuities represented by interface elements, and (ii) a new mechanical constitutive law for interface elements, which is able to reproduce the propagation of a crack induced by the development of an internal pressure exerted by solid reaction products formed within it. In addition, the numerical implementation of the diffusion-reaction formulation has been advantageously performed with clear separation between the boundary-value or ‘structural’ governing equations (i.e. continuity and concentration gradient relations), and the ‘constitutive’ (i.e. chemical) equations. The model is illustrated with some application examples.