Combined Finite-Discrete Element Modelling of the Influence of Bedding Planes on the Mechanical Behaviour of Opalinus Clay

The mechanical behaviour of Opalinus Clay, an indurated over-consolidated argillaceous rock, is heavily influenced by the presence of preferably oriented bedding planes. The objective of this study is to illustrate a new modeling approach aimed at capturing the bedding-induced anisotropy of Opalinus Clay as observed in the laboratory and in situ. The combined finite-discrete element method (FEM/DEM) is employed as the main simulation tool. FEM/DEM is a numerical method that models material damage processes by explicitly considering fracture nucleation and propagation. A transversely isotropic elastic constitutive law is used to account for different elastic property values in the direction parallel and perpendicular to the layering, while strength anisotropy is captured using a directional cohesive fracture model. Quantitative calibration and validation of the model is based on laboratory results of indirect tensile tests and uniaxial compression tests. Emergent strength and deformation properties, together with the simulated damage mechanisms, are shown to be in good agreement with experimental observations. Preliminary results relative to the analysis of a circular tunnel show that the approach is able to reproduce the main Excavation Damaged Zone (EDZ) formation mechanisms as observed in the field.