Characterization of clastic rock using a biconcave bond model of DEM

This paper presents a biconcave bond model to investigate the effect of the cementation between grains on the mechanical behavior of rock. The proposed model considers the shape of the bonds among particles that have a biconcave cement form, based on observations of microscopic rock images. The general equations of the proposed model are based on Dvorkin theory. The accuracy and efficiency of the bond model is improved in three ways. After the biconcave bond model is implemented in the discrete element method software Particle Flow Code in 2 Dimensions, a series of numerical uniaxial compression tests were performed to investigate the relationships between the micro‐ to macro‐parameters. The simulations revealed that the biconcave bond model reflects the effect of micro‐parameters, such as the elastic modulus and Poisson's ratio of the cement, on the macroscopic deformation of cemented granular material. Variations in the bond geometry caused extremely diverse macro‐mechanical behaviors. Experimental results concerning rock demonstrate that the biconcave bond model accurately captures the mechanical behavior of intact rock and supports an innovative method for investigating the relationships between the micro‐ and macro‐parameters of cemented granular material. Copyright © 2016 John Wiley & Sons, Ltd.     

On incremental non‐linearity in granular media: phenomenological and multi‐scale views (Part I)

On the basis of fundamental constitutive laws such as elasticity, perfect plasticity, and pure viscosity, many elasto‐viscoplastic constitutive relations have been developed since the 1970s through phenomenological approaches. In addition, a few more recent micro‐mechanical models based on multi‐scale approaches are now able to describe the main macroscopic features of the mechanical behaviour of granular media. The purpose of this paper is to compare a phenomenological constitutive relation and a micro‐mechanical model with respect to a basic issue regularly raised about granular assemblies: the incrementally non‐linear character of their behaviour. It is shown that both phenomenological and micro‐mechanical models exhibit an incremental non‐linearity. In addition, the multi‐scale approach reveals that the macroscopic incremental non‐linearity could stem from the change in the regime of local contacts between particles (from plastic regime to elastic regime) in terms of the incremental macroscopic loading direction. Copyright © 2005 John Wiley & Sons, Ltd.     

基于细观试验和离散元法的盐岩力学特性

当前盐岩的宏观力学模型通常是唯象模型,不能很好地解释盐岩受力变形破坏的真正物理基础。盐岩是由于化学沉积而形成的矿物集合体,是一种主要由NaCl和少量杂质组成的多晶体,其变形机制主要由晶粒与晶界的力学特性控制。通过扫描电镜（SEM）,获得盐岩晶粒的微细观结构特征,采用分子动力学方法和纳米压痕技术,确定盐岩晶粒和晶界的微细观力学参数;将盐岩晶粒作为块体,基于Voronoi多边形技术,建立盐岩的微细观数值模型;利用离散元方法,对盐岩试件在单轴压缩和直剪条件下的宏观力学行为进行了数值模拟。数值模拟结果与宏观力学试验结果吻合度高,表明基于盐岩微细观晶粒结构特征并结合离散元数值模拟的方法能够较好地研究盐岩的宏观力学性能及其材料物理基础。     

Micromechanical Modelling of Granular Materials: Effect of Particle Size and Gradation

In this paper the effects of maximum particle size, particle gradation/sorting and fabric on bulk mechanical behaviour of granular materials such as coarse grained soils and rockfills are investigated" from micromechanical considerations starting from the grain scale level, using numerical" simulations based on Discrete Element Modelling (DEM). Hydrostatic compaction and biaxial tests on 2-dimensional assemblies of discs with varying particle sizes and gradations were modelled using DEM. An examination of the constitutive behaviour of granular media considering" the particulate nature of the medium has been attempted to explain the effect of particle size and gradation. Simulation results on perfectly parallel graded assemblies indicate that with increase in the size of the particles, a marginal increase (or no increase) in the angle of internal friction is observed during biaxial loading conditions. A change to a wider gradation (keeping the minimum grain size the same) results in a decrease in the angle of internal friction and an increase in volumetric strain to a considerable extent. Based on micromechanical force and fabric parameters, the basis for the physical behaviour was established. This helps in understanding the physics of parallel gradation techniques.     

A three‐dimensional micromechanically based model

Granular materials react with complicated mechanical responses when subjected to external loading paths. This leads to sophisticated constitutive formulations requiring large numbers of parameters. A powerful and straightforward way consists in developing micro‐mechanical models embedding both micro‐scale and meso‐scale. This paper proposes a 3D micro‐mechanical model taking into account an intermediate scale (meso‐scale) that makes it possible to describe a variety of constitutive features in a natural way. The comparison between experimental tests and numerical simulations reveals the predictive capability of this model. Particularly, several simulations are carried out with different confining pressures and initial void ratios, based on the fact that the critical state is quantitatively described without requiring any critical state formulations and parameter. The model mechanism is also analyzed from a microscopic view, wherein the evolution of some key microscopic parameters is investigated. Copyright © 2017 John Wiley & Sons, Ltd.     

A simple three‐dimensional distinct element modeling of the mechanical behavior of bonded sands

This paper presents a simple three‐dimensional (3D) Distinct Element Method (DEM) for numerical simulation of the mechanical behavior of bonded sands. First, a series of micro‐mechanical tests on a pair of aluminum rods glued together by cement with different bond sizes were performed to obtain the contact mechanical responses of ideally bonded granular material. Second, a 3D bond contact model, which takes into account the influences of bond sizes, was established by extending the obtained 2D experimental results to 3D case. Then, a DEM incorporating the new contact model was employed to perform a set of drained triaxial compression tests on the DEM bonded specimens with different cement contents under different confining pressures. Finally, the mechanical behavior of the bonded specimens was compared with the available experimental results. The results show that the DEM incorporating the simple 3D bond contact model is able to capture the main mechanical behavior of bonded sands. The bonded specimen with higher cement content under lower confining pressure exhibits more pronounced strain softening and shear dilatancy. The peak and residual strengths, the apparent cohesion and peak/residual friction angles, and the position and slope of the critical state line increase with increase in cement content. Microscopically, bond breakage starts when the system starts to dilate and the maximum rate of bond breakage coincides with the maximum rate of dilation. Bond breakage is primarily due to tension‐shear failure and the percentage of such failures is independent of both confining pressure and cement content. Copyright © 2015 John Wiley & Sons, Ltd.     

Modeling the particle breakage of rockfill materials with the cohesive crack model

A practical combined finite–discrete element method was developed to simulate the breakage of irregularly shaped particles in granular geomaterials, e.g., rockfill. Using this method, each particle is discretized into a finite element mesh. The potential fracture paths are represented by pre-inserted cohesive interface elements (CIEs) with a progressive damage model. The Mohr–Coulomb model with a tension cut-off is employed as the damage initiation criterion to rupture the predominant failure mode occurs at the particle scale. Two series of biaxial tests were simulated for both the breakable and unbreakable particle assemblies. The two assemblies have identical configurations, with the exception that the former is inserted with CIEs and is breakable. The simulated stress–strain–dilation responses obtained for both assemblies are in agreement with experimental observations. We present a comprehensive study of the role of particle breakage on the mechanical behavior of rockfill materials at both the macroscopic and microscopic scales. The underlying mechanism of particle breakage can be explained by the force chain in the assemblies.     

Polygonal grain‐based distinct element modeling for mechanical behavior of brittle rock

The microstructure of rock was numerically reproduced by a polygonal grain‐based model, and its mechanical behavior was examined by performing the uniaxial compression test and Brazilian tests via the Universal Distinct Element Code. The numerical results of the model demonstrated good agreement with the experimental results obtained with rock specimens in terms of the stress–strain behavior, strength characteristics, and brittle fracture phenomenon. An encouraging result is that the grain‐based model‐Universal Distinct Element Code model can reproduce a low ratio of tensile to compressive strength of 1/20 to 1/10 without the need for an additional process. This finding is ascribed to the fact that the geometrical features of polygons can effectively capture the effects of angularity, finite rotation, and interlocking of grains that exist in reality. A numerical methodology to monitor the evolution of micro‐cracks was developed, which enabled us to examine the progressive process of the failure and distinguish the contribution of tensile cracking to the process from that of shear cracking. From the observations of the micro‐cracking process in reference to the stress–strain relation, crack initiation stress, and crack damage stress, it can be concluded that the failure process of the model closely resembles the microscopic observations of rock. We also carried out a parametric study to examine the relationships between the microscopic properties and the macroscopic behavior of the model. Depending on the micro‐properties, the model exhibited a variety of responses to the external load in terms of the strength and deformation characteristics, the evolution of micro‐cracks, and the post‐peak behavior. Copyright © 2016 John Wiley & Sons, Ltd.     

A micro‐mechanical study of peak strength and critical state

We present a micro‐mechanical analysis of macroscopic peak strength, critical state, and residual strength in two‐dimensional non‐cohesive granular media. Typical continuum constitutive quantities such as frictional strength and dilation angle are explicitly related to their corresponding grain‐scale counterparts (e.g., inter‐particle contact forces, fabric, particle displacements, and velocities), providing an across‐the‐scale basis for a better understanding and modeling of granular materials. These multi‐scale relations are derived in three steps. First, explicit relations between macroscopic stress and strain rate with the corresponding grain‐scale mechanics are established. Second, these relations are used in conjunction with the non‐associative Mohr–Coulomb criterion to explicitly connect internal friction and dilation angles to the micro‐mechanics. Third, the mentioned explicit connections are applied to investigate, understand, and derive micro‐mechanical conditions for peak strength, critical state, and residual strength. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical simulation of liquefaction behaviour of granular materials using Discrete Element Method

In this paper, numerical simulation of 3-dimensional assemblies of 1000 polydisperse sphere particles using Discrete Element Method (DEM) is used to study the liquefaction behaviour of granular materials. Numerical simulations of cyclic triaxial shear tests under undrained conditions are performed at different confining pressures under constant strain amplitude. Results obtained in these numerical simulations indicate that with increase in confining pressure there is an increase in liquefaction resistance.     

考虑颗粒棱角影响的直剪试验的离散元模拟

为了研究颗粒棱角对颗粒材料力学行为的影响,建立了具有不同棱角度的对称多面体颗粒,采用了一种简单并适合任意颗粒形状的接触本构模型,对三维离散元开源程序YADE进行了修改,研究了颗粒棱角度在模拟直剪试验中的影响以及接触力各向异性在剪切过程中的演化规律。研究结果表明,颗粒棱角度越小,颗粒间相互咬合自锁的作用越小,颗粒受剪更易转动,致使颗粒体系的剪切强度和剪胀性下降;竖向加载力越大,颗粒棱角度的影响越明显;法向接触力的各向异性在剪切过程中表现为先增后减最后趋向稳定的趋势;法向接触力的各向异性变化程度随颗粒棱角度的增大而增大。     

Study on the mechanical behavior of sands using 3D discrete element method with realistic particle models

Presented is a study on the geometrical characteristics of sand particles and the mechanical behavior of sand material under external loading. Based on computed tomography technique, a reconstruction method of granular particles was developed and used to build a database of 3D geometrical models for sand particles. The studied sand particles showed good regularities in morphological characteristics and thus were suitable to be used for the random generation of numerical samples. DEM tests using realistically shaped particles were proven to better simulate the mechanical behavior of the sample during elastoplastic loading stage, which was an issue for the simplified spherical particles. The generation, extension, and breakage of the force chains controlled the strain softening behavior of sands. Anisotropy analysis using the spherical harmonic series showed that the evolution of anisotropy directions and parameters corresponded well with the macroscopic mechanical behavior of the material. Pore volume computation based on Voronoi diagram was performed to illustrate the formation and evolution of localized shear zone. The mesoscopic analysis showed that particle shape significantly influences the mechanical behavior of sands and thus should be properly modeled in numerical simulations.

     

Adaptation of an existing constitutive model for describing mechanical behaviour of cemented granular soils

In this paper, an existing elastoplastic constitutive model, originally developed for granular soils, is adapted to describe the stress–strain behaviour of cemented granular soils. The existing model (CJS), due to its modular formulation, can be easily developed to take into account different supplementary behavioural aspects in soil mechanics. In the present study, the failure mechanism of the CJS model is modified by introducing the essential aspects in the behaviour of cemented granular soils in its formulation. All of the model parameters have clear physical meaning and can be identified using classical laboratory tests. A set of direct relations between model parameters and famous mechanical parameters of soils such as internal friction angle and cohesion at peak and residual states is presented. In order to validate the model, the results of triaxial and uniaxial tests in the compression and extension performed on cemented granular materials are used. The validation results indicate the good capability of the proposed model.     

Numerical simulation of drained triaxial test using 3D discrete element modeling

A discrete element modeling of granular material was carried out using a 3D spherical discrete model with a rolling resistance, in order to take into account the roughness of grains. The numerical model of Labenne sand was generated, and the desired porosity was obtained by a radius expansion method. Using numerical triaxial tests the micro-mechanical properties of the numerical material were calibrated in order to match the macroscopic response of the real material. Numerical simulations were carried out under the same conditions as the physical experiments (porosity, boundary conditions and loading). The pre-peak, peak and post-peak behavior of the numerical material was studied. The calibration procedure revealed that the peak stress of the sand sample does not only depend on local friction parameters but also on the rolling resistance. The larger the value of the applied rolling resistance, the higher the resulting stress peak. Furthermore, the deformational response depends strongly on local friction. The numerical results are quantitatively in agreement with the laboratory test results.     

不同含水率砂岩单轴压缩力学特性及波速法损伤

为了探究不同含水率砂岩单轴压缩下力学特性及损伤变化情况,首先对5种不同含水率砂岩进行静态单轴压缩实验,获得其物理力学参数;而后利用超声波检测仪测量该5种不同含水率砂岩的波速,将该批次砂岩加载到80%峰值强度后卸载到0,测量此时砂岩声波波速,利用声波法定义损伤,研究含水率对砂岩损伤的影响。结果表明:单轴压缩下,砂岩的峰值强度及弹性模量随含水率升高逐渐降低,而峰值应变变化呈相反趋势;总功及弹性能随含水率增加而下降,耗散能随含水率增加而增大。实验前声波波速随含水率升高逐渐减小,实验后声波波速随含水率升高而下降,且下降趋势较实验前更迅速;砂岩的残余塑性变形和内部损伤均随含水率升高逐渐增大。研究结果可为现场岩石损伤测试提供依据。      

The mechanical behaviour of synthetic, poorly consolidated granular rock under uniaxial compression

In order to isolate the effect of grain size and cementation on the mechanical behaviour of poorly consolidated granular rock, we prepared synthetic rock samples in which these two parameters were varied independently. Various proportions of sand, Portland cement and water were mixed and cast in a mold. The mixture was left pressure-free during curing, thus ensuring that the final material was poorly consolidated. We used two natural well-sorted sands with grain sizes of 0.22 and 0.8 mm. The samples were mechanically tested in a uniaxial press. Static Young's modulus was measured during the tests by performing small stress excursions at discrete intervals along the stress–strain curves. All the samples exhibited nonlinear elasticity, i.e., Young's modulus increased with stress. As expected, we found that the uniaxial compressive strength increased with increasing cement content. Furthermore, we observed a transition from grain size sensitivity of strength at cement content less than 20–30% to grain size independence above this value. The measured values of Young's modulus are well explained by models based on rigid inclusions embedded in a soft matrix, at high cement content, and on cemented grain-to-grain contacts, at low cement content. Both models predict grain size independence in well-sorted cemented sands. The observed grain size sensitivity at low cement content is probably due to microstructural differences between fine- and coarse-grained materials caused by small differences in grain sorting quality.     

A micromorphic elastoplastic model and finite element simulation on failure behaviors of granular materials

One of the purposes in this study is to develop a modified micromorphic continuum model for granular materials on the basis of a micromechanics approach. A symmetric curvature tensor is proposed in this model, and a symmetric couple stress tensor is derived conjugating the symmetric curvature tensor. In addition, a correct derivation is presented to obtain the symmetric stress tensor conjugated with the symmetric strain tensor. The modified model provides a complete deformation mode for granular materials by considering the decomposition for motions (displacement and rotation) of particles. Consequently, the macroscopic constitutive relationships and constitutive moduli are derived in expressions of the microstructural information. Furthermore, the balance equations and boundary conditions are obtained for the modified micromorphic model. By considering the extended Drucker-Prager yield criterion, the micromorphic elastoplastic model is developed. Another purpose of this study is to derive the finite element formulation for the developed micromorphic elastoplastic model. Based on the ABAQUS user element (UEL) interface, numerical simulations investigated the load-displacement relationship and the strain localization behavior of granular materials and investigated the influence of microscopic parameters in the micromorphic model on these macroscopic mechanical responses. Numerical results illustrate the presented model's capability of simulating the strain-softening and strain localization behaviors, and the capability of considering the influence of microstructural information on the macroscopic mechanical behaviors of granular materials.