Hadrien Laubie

Fracture toughness of calcium-silicate-hydrate from molecular dynamics simulations

Abstract Concrete is the most widely manufactured material in the world. Its binding phase, calcium–silicate–hydrate (C–S–H), is responsible for its mechanical properties and has an atomic structure fairly similar to that of usual calcium silicate glasses, which makes it appealing to study this material with tools and theories traditionally used for non-crystalline solids. Here, following this idea, we use molecular dynamics simulations to evaluate the fracture toughness of C–S–H, inaccessible experimentally. This allows us to discuss the brittleness of the material at the atomic scale. We show that, at this scale, C–S–H breaks in a ductile way, which prevents one from using methods based on linear elastic fracture mechanics. Knowledge of the fracture properties of C–S–H at the atomic scale opens the way for an upscaling approach to the design of tougher cement paste, which would allow for the design of slender environment-friendly infrastructures, requiring less material.

Atomic-scale modelling of elastic and failure properties of clays

The elastic and failure properties of a typical clay, illite, are investigated using molecular simulation. We employ a reactive (ReaxFF) and a non-reactive (ClayFF) force field to assess the elastic properties of the clay. As far as failure is concerned, ReaxFF was used throughout the study; however, some calculations were also performed with ClayFF. A crack parallel to the clay layers is found to have low fracture resistance when submitted to a tensile loading perpendicular to the crack. The mechanism of both yield and fracture failures is decohesion in the interlayer space. In contrast, under shear loading, the nanoscale failure mechanism is a stick-slip between clay layers. No fracture propagation is observed as the clay layers slide on top of each other. The low fracture resistance in mode I and the stick-slip failure in mode II are both the consequence of the lack of chemical bonds between clay layers where the cohesion is provided by non-covalent interactions. This work, which provides a description of the failure of clays at the microscopic scale, is the first step towards describing the failure of clays at a larger scale where the polycrystalline distribution of clay grains must be taken into account.

Fracture Mechanisms in Organic-Rich Shales: Role of Kerogen

In this work we study role of kerogen in the fracture properties of organic-rich shales and, in particular, in the ductility of shales. The presence of kerogen and clays in shale is known to increase the ductility. We propose here a multiscale approach to develop a fine understanding of shale ductility from the molecular scale. We develop and validate a methodology at the molecular scale that can capture the toughness and ductility of a material. We apply this methodology successfully to a silica polymorph and to a kerogen analog, and we confirm the significant ductility of kerogen. Interestingly the silica-kerogen interface exhibits a similar ductility, which is central for the properties of the heterogeneous shale. Finally, we consider a tentative upscaling considering the pull out phenomenon as a likely mechanism of fracture of the shale.

Plane-Strain Crack Problem in Transversely Isotropic Solids for Hydraulic Fracturing Applications

AbstractThis paper aims at understanding and predicting how pressurized cracks propagate in anisotropic brittle solids, a situation frequently encountered in hydraulic fracturing. Special attention is paid to transverse isotropy, often used to model shale. Although the theory of linear elastic fracture mechanics of anisotropic solids is well established at present, this paper shows that the application of Muskhelishvili’s formalism to Lekhnitskii’s anisotropic complex potentials provides a powerful tool to extend the validity of the classical tools of isotropic fluid-driven crack models to the anisotropic case, provided that the appropriate elastic constants are used. These elastic constants are identified and derived in closed form for transversely isotropic solids. The constants are shown to be directly related to quantities easily measured in a laboratory at macroscopic scale through indentation tests and acoustic measurements. Moreover, several crack-kinking criteria are compared. Contrary to the isot...

Stress Transmission and Failure in Disordered Porous Media

By means of extensive lattice-element simulations, we investigate stress transmission and its relation with failure properties in increasingly disordered porous systems. We observe a non-Gaussian broadening of stress probability density functions under tensile loading with increasing porosity and disorder, revealing a gradual transition from a state governed by single-pore stress concentration to a state controlled by multipore interactions and metric disorder. This effect is captured by the excess kurtosis of stress distributions and shown to be nicely correlated with the second moment of local porosity fluctuations, which appears thus as a (dis)order parameter for the system. By generating statistical ensembles of porous textures with varying porosity and disorder, we derive a general expression for the fracture stress as a decreasing function of porosity and disorder. Focusing on critical sites where the local stress is above the global fracture threshold, we also analyze the transition to failure in terms of a coarse-graining length. These findings provide a general framework which can also be more generally applied to multiphase and structural heterogeneous materials.

Effective Potentials and Elastic Properties in the Lattice-Element Method: Isotropy and Transverse Isotropy

AbstractLattice approaches have emerged as a powerful tool to capture the effective mechanical behavior of heterogeneous materials using harmonic interactions inspired from beam-type stretch and ro...

Mesoscale Poroelasticity of Heterogeneous Media

AbstractThe poromechanics of heterogeneous media is reformulated in a discrete framework using the lattice element method (LEM) that accounts for the presence of interfaces as well as local microte...