Joints in unsaturated rocks: Thermo-hydro-mechanical formulation and constitutive behaviour

A formulation for the coupled analysis of thermo-hydro-mechanical (THM) problems in joints is first presented. The work involves the establishment of equilibrium and mass and energy balance equations. Balance equations were formulated taking into account two phases: water and air. The joint element developed was implemented in a general purpose finite element computer code for THM analysis of porous media (Code_Bright). The program was then used to study a number of cases ranging from laboratory tests to large scale in situ tests. A numerical simulation of coupled hydraulic shear tests of rough granite joints is first presented. The tests as well as the model show the coupling between permeability and the deformation of the joints. The experimental investigation was focused on the effects of suction on the mechanical behaviour of rock joints. Laboratory tests were performed in a direct shear cell equipped with suction control. Suction was imposed using a vapour forced convection circuit connected to the cell and controlled by an air pump. Artificial joints of Lilla claystone were prepared. Joint roughness of varying intensity was created by carving the surfaces in contact in such a manner that rock ridges of different tip angles were formed. These angles ranged from 0° (smooth joint) to 45° (very rough joint profile). The geometric profiles of the two surfaces in contact were initially positioned in a “matching” situation. Several tests were performed for different values of suctions (200, 100, and 20 MPa) and for different values of vertical stresses (30, 60, and 150 kPa). A constitutive model including the effects of suction and joint roughness is proposed to simulate the unsaturated behaviour of rock joints. The new constitutive law was incorporated in the code and experimental results were numerically simulated.