DEM of triaxial tests on crushable cemented sand

Using the discrete element method, triaxial simulations of cemented sand consisting of crushable particles are presented. The triaxial model used features a flexible membrane, allowing realistic deformation to occur, and cementation is modelled using inter-particle bonds. The effects of particle crushing are explored, as is the influence of cementation on the behaviour of the soil. An insight to the effects that cementation has on the degree of crushing is presented

Discrete Element Modeling of Compound Rockfall Fence Nets

Compound mesh panels are structures in which two different nets geometries are employed: a main mesh that provides the bearing capacity and a weaker mesh with a thin sieve size to catch smaller blocks that can pass through otherwise. Typically, only the effect of the main mesh is investigated, and the weaker mesh is considered to provide negligible structural resistance. In this paper, after a calibration procedure, numerical simulations of quasi-static punch tests and a dynamic block impact on a composite double-twist and strand rope mesh are performed. The results show that, under dynamic conditions, the presence of the finer mesh lowers the peak force acting on the main mesh. This effect is not found under quasi-static conditions and has important repercussions on the overall structural resistance as the energy dissipation mechanism reduces the stress on the mesh fence posts

3D Voronoi Tessellation for the Study of Mechanical Behavior of Rocks at Different Scales

Numerical investigation of crack damage development and microfracturing in brittle rocks is a widely studied topic, given the number of applications involved. In the framework of the Discrete Element Method (DEM) formulation, the grain-based distinct element model with random polygonal blocks can represent an alternative to the Bonded-Particle Model (BPM) based on particles. Recently, a new engine called Neper has been made available for generating 3D Voronoi grains. The aim of this study is to investigate the applicability of a Neper-based 3D Voronoi tessellation technique for the simulation of the mechanical macro response of rocks. Simulation of unconfined compression tests on synthetic specimens is conducted and a calibration procedure tested. The issue related to scale effects is also addressed, with an application to the case study of a deep geothermal reservoir

Evaluation of Mode I Fracture Toughness Assisted by the Numerical Determination of K-Resistance

The fracture toughness of a rock often varies depending on the specimen shape and the loading type used to measure it. To investigate the mode I fracture toughness using semi-circular bend (SCB) specimens, we experimentally studied the fracture toughness using SCB and chevron bend (CB) specimens, the latter being one of the specimens used extensively as an International Society for Rock Mechanics (ISRM) suggested method, for comparison. The mode I fracture toughness measured using SCB specimens is lower than both the level I and level II fracture toughness values measured using CB specimens. A numerical study based on discontinuum mechanics was conducted using a two-dimensional distinct element method (DEM) for evaluating crack propagation in the SCB specimen during loading. The numerical results indicate subcritical crack growth as well as sudden crack propagation when the load reaches the maximum. A K-resistance curve is drawn using the crack extension and the load at the point of evaluation. The fracture toughness evaluated by the K-resistance curve is in agreement with the level II fracture toughness measured using CB specimens. Therefore, the SCB specimen yields an improved value for fracture toughness when the increase of K-resistance with stable crack propagation is considered

The influence of fines content and size-ratio on the micro-scale properties of dense bimodal materials

This paper considers factors influencing the fabric of bimodal or gap-graded soils. Discrete element method simulations were carried out in which the volumetric fines content and the size ratio between coarse and fine particles were systematically varied. Frictionless particles were used during isotropic compression to create dense samples; the coefficient of friction was then set to match that of spherical glass beads. The particle-scale data generated in the simulations revealed key size ratios and fines contents at which transitions in soil fabric occur. These transitions are identified from changes in the contact distributions and stress-transfer characteristics of the soils and by changes in the size of the void space between the coarse particles. The results are broadly in agreement with available experimental data on minimum void ratio and contact distributions. The results have implications for engineering applications including assessment of the internal stability of gap-graded soils in embankment dams and flood embankments

Modelling of hydraulic fracturing process by coupled discrete element and fluid dynamic methods

A three-dimensional model is presented and used to reproduce the laboratory hydraulic fracturing test performed on a thick-walled hollow cylinder limestone sample. This work aims to investigate the implications of the fluid flow on the behaviour of the micro-structure of the rock sample, including the material strength, its elastic constants and the initialisation and propagation of fractures. The replication of the laboratory test conditions has been performed based on the coupled Discrete Element Method (DEM) and Computational Fluid Dynamics scheme. The numerical results are in good agreement with the experimental data, both qualitatively and quantitatively. The developed model closely validates the overall behaviour of the laboratory sample, providing a realistic overview of the cracking propagation towards total collapse as well as complying with Lame’s theory for thick-walled cylinders. This research aims to provide some insight into designing an accurate DEM model of a fracturing rock that can be used to predict its geo-mechanical behaviour during Enhanced Oil Recovery applications