颗粒材料的微观临界状态理论模型

存在临界状态是颗粒材料的一个重要特性。基于孔隙胞元的颗粒离散元方法对二维颗粒体进行双轴加载数值试验,在详细分析数值模拟结果的基础上,从微观几何组构的角度揭示了临界状态的存在机制。基于剪胀性原理,提出了以接触价键表征的微观临界状态理论模型,得到了接触价键与塑性剪切应变的关系表达式,理论模型的结果和二维离散元数值模拟得到的结果吻合较好。通过比较不同情况下数值结果和理论模型中的参数,得到以下结论：表征微观临界状态的参数（临界接触价键和达到临界状态所需要的塑性剪切应变）依赖于颗粒体的微观特性,如颗粒形状、表面摩擦性质、颗粒体的围压和初始孔隙比。     

颗粒材料剪胀性的微观力学分析

剪胀性是颗粒材料在加载过程中表现出来的重要变形特性。以孔隙胞元描述颗粒材料内部结构的最小单元,通过对单个孔隙胞元进行剪切受力分析,探讨了剪切过程中颗粒材料体积的改变对应力比和单个孔隙胞元形状的依赖关系,解释了排列密实的颗粒材料在剪切过程中先压缩后剪胀的微观机制。用离散元数值模拟得到了在双轴剪切过程中单个孔隙胞元形状以及孔隙胞元体积变形的演化过程。离散元数值结果表明,加载过程中孔隙胞元形状由初始各向同性到沿大主应力方向变大变长、体积变形先压缩后膨胀,并且体积变形在加载过程中存在局部化现象,体积变化大的孔隙胞元在较大变形时,排列成倾斜的窄带。综合孔隙胞元的受力分析和离散元数值结果表明,致密排列颗粒材料的剪胀性与微观尺度上孔隙胞元的几何结构及其内部的力链传递方式密切相关。     

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.     

散粒介质临界状态细观力学结构特征的数值模拟与分析

基于离散单元法对不同密实度理想散粒体进行了双轴剪切试验的宏微观数值模拟,通过网格剖分将Voronoi多边形表征的loop单元作为散粒介质细观力学结构的基本单元,模拟了剪切过程中不同类型loop单元数量、几何形态和力学特征的演化过程,并重点分析了临界状态时散粒介质的细观力学结构特征。模拟结果显示,初始密实度不同的试样在向临界状态发展的过程中,高阶单元与低阶单元的发展规律完全不同,不同初始密实度试样中同阶loop单元的发展规律也不相同,但同阶loop单元的数量比例、几何形态、颗粒接触力及单元内滑动率最终均达到了各自的临界状态。从细观角度分析,散体介质的临界状态是高阶和低阶loop单元在荷载作用下相互转化的结果,是所有loop单元物理力学状态的综合平均与外在表现,临界状态时不同阶数的loop单元处于一个动态平衡状态,宏观上表现为常剪应力和常体积下剪切变形的不断发展。数值模拟结果也表明,loop细观结构单元包含了丰富的信息,其数量、几何形态、受力特征及接触稳定性的发展与散粒体的强度、剪胀以及临界状态的发展密切相关,可以将其作为散粒介质细观尺度的分析单元。     

Bounding surface plasticity model for stress-strain and grain-crushing behaviors of rockfill materials

Crushing of grains can greatly influence the strength,dilatancy,and stress-strain relationship of rockfill materials.The critical state line(CSL)in the void ratio versus mean effective stress plane was extended to the breakage critical state plane(BCSP).A state void-ratio-pressure index that incorporated the effect of grain crushing was proposed according to the BCSP.Rowe’s stress-dilatancy equation was modified by adding the breakage voidratio-pressure index,which was also incorporated into the formulations of the bounding stress ratio and plastic modulus.A BCSP-based bounding surface plasticity model was proposed to describe the state-dependent stressstrain behaviors and the evolution of grain crushing during shearing process of rockfill materials,and was shown to sufficiently capture the breakage phenomenon.     

Analysis of soil‐structural interface behavior using three‐dimensional DEM simulations

The shear behavior at the interface between the soil and a structure is investigated at the macroscale and particle‐scale levels using a 3‐dimensional discrete element method (DEM). The macroscopic mechanical properties and microscopic quantities affected by the normalized interface roughness and the loading parameters are analyzed. The macro‐response shows that the shear strength of the interface increases as the normalized roughness of the interface increases, and stress softening and dilatancy of the soil material are observed in the tests that feature rough interfaces. The particle‐scale analysis illustrates that a localized band characterized by intense shear deformation emerges from the contact plane and gradually expands as shearing progresses before stabilizing at the residual stress state. The thickness of the localized band is affected by the normalized roughness of the interface and the normal stress, which ranges between 4 and 5 times that of the median grain diameter. A thicker localized band is formed when the soil has a rough shearing interface. After the localized band appears, the granular material structuralizes into 2 regions: the interface zone and the upper zone. The mechanical behavior in the interface zone is representative of the interface according to the local average stress analysis. Certain microscopic quantities in the interface zone are analyzed, including the coordination number and the material fabric. Shear at the interface creates an anisotropic material fabric and leads to the rotation of the major principal stress.     

Critical state of granular materials at mesoscale: Transition from local to global

Critical state theory is widely accepted for understanding and interpreting the behavior of granular materials such as sandy soils. Recently, motivated by the development of a multiscale analysis method, this paper presents a discrete element investigation on the critical state behavior of granular materials on two scales level, that is, the macroscale and mesoscale. Loops are taken as the basic mesostructures, and the concept of mesoscale critical state is proposed on terms of the local stress and void ratio. Based on series of numerical simulations with varying densities and confining pressures, the characters of critical state lines (CSLs) for each loop type are analyzed. An analytical relationship of the transition from local to global is then obtained for critical state. Besides, the influence of interparticle friction coefficient on CSL is also discussed on both scales. The prominent observation is that each type of loops achieves its own critical state, from which a unique global CSL can be obtained. The CSLs of different loops exhibit same shapes but move upwards with the increasing order.     

胶结结构性土统一硬化模型

结构性土颗粒间的胶结使试样剪切破坏最终应力比高于相应重塑土,也限制了试样剪切时体积应变的自由发挥。在考虑结构垮塌为主的结构性土统一硬化(UH)模型基础上,将应力空间中静止的临界状态线扩展为动态的移动临界状态线。据此,通过建立新的屈服面方程并修正剪胀方程,将结构性土统一硬化(UH)模型扩展为胶结结构性土统一硬化(UH)模型。相对于原模型,新模型增加了1个模型参数,即初始胶结应力,反映土颗粒之间的初始胶结作用。通过4种结构性土试验数据与模型预测对照表明：所提模型能够较合理地描述结构性土等向压缩、常规三轴排水与不排水剪切等特性。     

A predictive deep learning framework for path-dependent mechanical behavior of granular materials

As we transition into an era of data generation and collection, empirical summaries in the classical continuum modeling of granular materials cannot take full advantage of the increasingly larger data sets. This work presents a data-driven model for modeling granular materials, with the material data being extracted from discrete element method (DEM) simulations. A long short-term memory (LSTM) network is then employed to learn the mechanical behaviors of granular materials from the material dataset. Particular emphasis is placed on three elements: modification of LSTM unit cell, phase space sampling, and material history parameterization. The LSTM unit cell is modified so that the initial hidden state can be specified as the initial states of granular materials. Massive DEM simulations are performed to consider the effects of particle size distribution, initial density, confining pressure, and loading path on the mechanical behaviors of granular materials. The history-dependency of the granular materials is well represented by the architecture of the LSTM network and internal variable-based history parameterization. We compare the model predictions against DEM simulations to assess the performance of the proposed data-driven model. The results demonstrate that the model can predict the material behaviors of granular materials with different microstructures and initial states and reproduce the material responses under complex nonmonotonic loading paths. This data-driven model exhibits good generalization ability and high prediction accuracy in various situations.

     

Cavity expansion analyses of crushable granular materials with state‐dependent dilatancy

Crushability is one of the important behaviors of granular materials particularly under high stress states, and affects both the deformability and strength of the materials that are in essence associated with state‐dependent dilatancy. In this presentation, first, a new critical state model is proposed to take into account the three different modes of compressive deformation of crushable granular materials, i.e. particle rearrangement, particle crushing and pseudo‐elastic deformation. Second, the governing equations for cavity expansion in crushable granulates are introduced, in which the state‐dependent dilatancy as well as the bounding surface plasticity model are used. Then, the procedure to obtain semi‐analytical solutions to cavity expansion in the material is described in detail, in which a commercial differential equation solver is employed. Finally, cavity expansion analyses are carried out on Toyoura sand, a well‐documented granular material, to demonstrate the effects of crushability and state‐dependent dilatancy. The study shows that particle crushing does occur at both high stress and critical states and affects the stress fields and the deformation behavior of the material surrounding the cavity in association with state‐dependent dilatancy. This leads to conclusion that particle crushing and state‐dependent dilatancy have to be taken into account when cavity expansion theory is used to interpret cone penetration tests and pressuremeter tests. Copyright © 2011 John Wiley & Sons, Ltd.     

颗粒大小对颗粒材料力学行为影响初探

利用一种特殊颗粒材料-玻璃珠进行了一系列室内直剪试验,研究颗粒大小对颗粒材料力学行为的影响。试验一共考虑了3条近乎平行的级配曲线和4种颗粒摩擦情况：干燥状态、水浸润状态、水淹没状态和油浸润状态。试验结果表明,颗粒大小对颗粒材料的力学行为有显著影响,剪胀性随着粒径的增大而增强。为考虑颗粒大小对剪胀性的影响,提出了一种新的剪胀关系式。在该剪胀关系式中,剪胀系数为依赖于颗粒大小和颗粒摩擦等颗粒基本性质的变量。试验研究同时表明临界状态摩擦角随着颗粒大小的增加而增加。此外,从颗粒细观运动的角度提出了颗粒滑动的功能模型,推导出了功能方程,并以此揭示了颗粒大小对临界状态摩擦角影响的细观机制。     

Stress–dilatancy based modelling of granular materials and extensions to soils with crushable grains

Stress–dilatancy relations have played a crucial role in the understanding of the mechanical behaviour of soils and in the development of realistic constitutive models for their response. Recent investigations on the mechanical behaviour of materials with crushable grains have called into question the validity of classical relations such as those used in critical state soil mechanics. In this paper, a method to construct thermodynamically consistent (isotropic, three‐invariant) elasto‐plastic models based on a given stress–dilatancy relation is discussed. Extensions to cover the case of granular materials with crushable grains are also presented, based on the interpretation of some classical model parameters (e.g. the stress ratio at critical state) as internal variables that evolve according to suitable hardening laws. Copyright © 2004 John Wiley & Sons, Ltd.     

Microscopic contact model of lunar regolith for high efficiency discrete element analyses

The grains of lunar regolith are characterized with rough surfaces, angular shapes and mutual adhesions due to short-range interactions. These features control the macroscopic mechanical behavior of lunar regolith but have not been completely captured by contact models in previous Discrete Element Method (DEM) analyses. In this paper, a simplified two-dimensional microscopic contact model is proposed for high efficiency DEM analyses of lunar regolith. The model consists of three components in the normal, tangential and rolling directions respectively, plus two new parameters. A shape parameter is used to control the rolling resistance ability at the contact area between two particles to capture the features of grain shape and interlocking. The second parameter, micro-separation, which denotes the nominal minimum distance between the molecules of the two contacting particles, is introduced to account for van der Waals force as the major component of the short-range interactions that contribute to the adhesion of regolith grains in lunar environment conditions. The novel model has been implemented in a two-dimensional DEM code for numerical simulations of biaxial compression tests on lunar regolith. The effects of interparticle friction, grain shape, lunar environment conditions and void ratio on the strength of lunar regolith were numerically investigated. The results show that soils in the simulated lunar environment exhibit greater strength and more apparent strain-softening and shear dilatancy than on the Earth. The proposed model can capture the main features of the mechanical behavior of lunar regolith (apparent cohesion and high peak friction angle) and a wide range of strength indices can be obtained by the contact model.     

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.     

Induced fabric under cyclic and rotational loads in a strain space multiple mechanism model for granular materials

The strain space multiple mechanism model idealizes the behavior of granular materials on the basis of a multitude of virtual simple shear mechanisms oriented in arbitrary directions. Within this modeling framework, the virtual simple shear stress is defined as a quantity dependent on the contact distribution function as well as the normal and tangential components of interparticle contact forces, which evolve independently during the loading process. In other terms, the virtual simple shear stress is an intermediate quantity in the upscaling process from the microscopic level (characterized by contact distribution and interparticle contact forces) to the macroscopic stress. The stress space fabric produces macroscopic stress through the tensorial average. Thus, the stress space fabric characterizes the fundamental and higher modes of anisotropy induced in granular materials. Herein, the induced fabric is associated with monotonic and cyclic loadings, loading with the rotation of the principal stress, and general loading. Upon loading with the rotation of the principal stress axis, some of the virtual simple shear mechanisms undergo loading whereas others undergo unloading. This process of fabric evolution is the primary cause of noncoaxiality between the axes of principal stresses and strains. Although cyclic behavior and behavior under the rotation of the principal stress axis seem to originate from two distinct mechanisms, the strain space multiple mechanism model demonstrates that these behaviors are closely related through the hysteretic damping factor. Copyright © 2011 John Wiley & Sons, Ltd.     

A constitutive model for granular materials with grain crushing and its application to a pyroclastic soil

A constitutive model for granular materials is developed within the framework of strain–hardening elastoplasticity, aiming at describing some of the macroscopic effects of the degradation processes associated with grain crushing. The central assumption of the paper is that, upon loading, the frictional properties of the material are modified as a consequence of the changes in grain size distribution. The effects of these irreversible microscopic processes are described macroscopically as accumulated plastic strain. Plastic strain drives the evolution of internal variables which model phenomenologically the changes of mechanical properties induced by grain crushing by controlling the geometry of the yield locus and the direction of plastic flow. An application of the model to Pozzolana Nera is presented. The stress–dilatancy relationship observed for this material is used as a guidance for the formulation of hardening laws. One of the salient features of the proposed model is its capability of reproducing the stress–dilatancy behaviour observed in Pozzolana Nera, for which the minimum value of dilatancy always follows the maximum stress ratio experienced by the material. Copyright © 2002 John Wiley & Sons, Ltd.     

Constitutive modelling of granular materials using a contact normal-based fabric tensor

This paper presents a fabric tensor-based bounding surface model accounting for anisotropic behaviour (e.g. the dependency of peak strength on loading direction and non-coaxial deformation) of granular materials. This model is developed based on a well-calibrated isotropic bounding surface model. The yield surface is modified by incorporating the back stress which is proportional to a contact normal-based fabric tensor for characterising fabric anisotropy. The evolution law of the fabric tensor, which is dependent on both rates of the stress ratio and the plastic strain, rules that the material fabric tends to align with the loading direction and evolves towards a unique critical state fabric tensor under monotonic shearing. The incorporation of the evolution law leads to a rotational hardening of the yield surface. The anisotropic critical state is assumed to be independent of the initial values of void ratio and fabric tensor. The critical state fabric tensor has the same intermediate stress ratio (i.e. b value) and principal directions as the critical state stress tensor. A non-associated flow rule in the deviatoric plane is adopted, which is able to predict the non-coaxial flow naturally. The stress–strain relation and fabric evolution of model predictions show a satisfactory agreement with DEM simulation results under monotonic shearing with different loading directions. The model is also validated by comparing with laboratory test results of Leighton Buzzard sand and Toyoura sand under various loading paths. The comparison results demonstrate encouraging applicability of the model for predicting the anisotropic behaviour of granular materials.

     

Evolution of induced fabric in a strain space multiple mechanism model for granular materials

The strain space multiple mechanism model idealizes the behavior of granular materials based on a multitude of virtual simple shear mechanisms oriented in arbitrary directions. Within this modeling framework, the virtual simple shear stress is defined as a quantity that depends on the contact distribution function as well as the normal and tangential components of inter‐particle contact forces, which evolve independently during the loading process. In other terms, the virtual simple shear stress is an intermediate quantity in the upscaling process from the microscopic level (characterized by the contact distribution and inter‐particle contact forces). The stress space fabric (i.e. the orientation distribution of the virtual simple shear stress) produces macroscopic stress through the tensorial average. Thus, the stress space fabric characterizes the fundamental and higher modes of anisotropy induced in granular materials. Comparing an induced fabric associated with the biaxial shear of plane granular assemblies obtained via a simulation using Discrete Element Method to the strain space multiple mechanism model suggests that the strain space multiple mechanism model has the capability to capture the essential features in the evolution of an induced fabric in granular materials. Copyright © 2012 John Wiley & Sons, Ltd.