Strength anisotropy of granular material consisting of perfectly round particles

The strength anisotropy of granular materials deposited under gravity has mostly been attributed to elongated particles' tendency to align long axes along the bedding plane direction. However, recent experiments on near‐spherical glass beads, for which preferred particle alignment is inapplicable, have exhibited surprisingly strong strength anisotropy. This study tests the hypothesis that certain amount of fabric anisotropy caused by the anisotropic stress during deposition under gravity can be locked in a circular‐particle deposit. Such locked‐in fabric anisotropy can withstand isotropic consolidation and leads to significant strength anisotropy. 2D discrete element method simulations of direct shear tests on circular‐particle deposits are conducted in this study, allowing for the monitoring of both stress and fabric. Simulations on both monodispersed and polydispersed circular‐particle samples generated under downward gravitational acceleration exhibit clear anisotropy in shear strength, thereby proving the hypothesis. When using contact normal‐based and void‐based fabric tensors to quantify fabric anisotropy in the material, we find that the intensity of anisotropy is discernible but low prior to shearing and is dependent on the consolidation process and the dispersity of the sample. The fact that samples with very low anisotropy intensity measurements still exhibit fairly strong strength anisotropy suggests that current typical contact normal‐based and void‐based second‐order fabric tensor formulations may not be very effective in reflecting the anisotropic peak shear strength of granular materials. Copyright © 2017 John Wiley & Sons, Ltd.     

控制应力路径散体材料真三轴试验强度特征的离散元模拟与细观机制分析

通过一系列真三轴离散元数值试验,模拟了不同应力路径下的等b试验中散体材料的强度特征。根据模拟结果详细地分析了三维应力条件下中主应力和应力路径对散体材料峰值强度的影响,研究了峰值摩擦角、峰值应力比的变化规律,并根据真应力的概念和组构张量的演化结果分析了散体材料的强度成因。研究表明,在不同类型的数值试验中峰值偏应力随b参数的变化规律不同,但采用初始围压归一化后的应力-应变曲线规律一致。峰值强度线的斜率只与b值有关而与应力路径无关,且随着b值的增加,峰值应力比qf /pf逐渐减小,数值模拟结果与室内试验结果吻合较好;随着应变的发展,数值试样的组构也随之发生变化,产生了明显的应力诱发各向异性;散体的强度为颗粒摩擦及材料各向异性共同作用的结果;理论上,组构比-应力比坐标系中破坏点位置仅取决于颗粒摩擦角 ,而数值模拟结果与理论值的差异源于颗粒间咬合和滚动摩擦的影响,其影响与颗粒表面摩擦系数有关,也受空间应力状态的影响。     

二维粒状材料双轴压缩试验研究

为了研究粒状材料的各向异性力学行为与细观组构演化之间的关系,采用自主研发的双轴压缩试验系统,以圆形和椭圆形截面的金属棒状材料组成的二维堆积体为试验对象,对不同大主应力方向角?（沉积面与大主应力作用面的夹角）的试样进行了各向等压、常侧向压力、等p剪切3种应力路径试验,并通过分析试样在不同变形阶段的数字照片得到了其细观组构演化规律。发现对于椭圆形截面的试样存在一个卓越剪切方向,随剪应变增大,颗粒长轴呈现出向该方向偏转的趋势,并且在大变形条件下沿该方向形成剪切带;卓越剪切方向与沉积面方向关系不大,而与大主应力作用面方向夹角约为45°+ /2, 为残余内摩擦角;随卓越剪切方向与沉降方向夹角的不同,颗粒偏转程度的不同是导致剪胀特性和峰值强度各向异性的主要原因     

初始组构各向异性对砂土力学特性及临界状态的影响

采用离散元方法,利用半径扩展法和重力沉积法分别生成具有初始各向同性和各向异性内结构的试样,并开展三轴不排水压缩和拉伸试验,研究不同制样方法产生的初始各向异性对砂土宏微观力学特性及其临界状态的影响。运用组构张量对砂土的各向异性进行量化,分析不同初始组构各向异性对组构张量演化的影响并确定了组构张量的临界值。试验结果表明：初始组构各向异性对试样的剪胀性有重要影响,由于受重力影响形成初始各向异性,其各向异性程度越大、组构方向与加载方向越一致,剪胀性越显著;初始组构各向异性对试样的临界状态没有影响,砂土的组构张量具有唯一的临界状态值。     

人工制备初始应力各向异性结构性土方法探讨

通过对结构性黏土的研究可以掌握天然土受荷过程中的变形破坏过程,从而为考虑土结构性的结构物的设计、地基加固等提供依据。近年来随着高、深和大型建筑的兴建,结构性黏土的研究变得尤为重要。天然土都具有结构性和各向异性。发展了一种能够考虑初始应力各向异性影响的结构性土的人工制样方法。通过对原料粉质黏土中添加水泥形成颗粒之间的胶结作用,添加盐粒并溶解后形成大孔隙组构分布和在水化过程中的侧限应力状态的结构性土样的端部施加竖向荷载,从而人工制备了具有初始应力各向异性的结构性土样。然后对初始均质结构性土样、初始应力各向异性结构性土样和重塑土样进行了三轴固结排水剪切试验,初步分析了初始应力对结构性土样的应力-应变特性的影响和初始应力各向异性结构性土的破损机制。     

A micromechanics-based model for estimating localized failure with effects of fabric anisotropy

This paper presents a micromechanics-based approach to investigate the effects of fabric anisotropy on the behavior of localized failure in granular materials. Based on a micromechanical analysis, the origin of deviatoric stress is decomposed into two components: contact force anisotropy and fabric anisotropy. Using a micro–macro approach, the back stress is interpreted as an contribution to the change of the fabric’s principal direction. The evolution of the back stress is deduced from the stress–fabric relationship and determined with reference to the deviation of the principal directions between the rate of the reduced stress tensor and the actual reduced stress tensor. With this micro–macro framework, a mixed (isotropic–kinematic) hardening model is developed based on the classical isotropic hardening theory. A laboratory simple shear test is first analyzed to validate the proposed model and illustrate the kinematic-hardening mechanism which is usually displayed under non-proportional loading. The analysis further focuses on the anisotropic aspect of localized failure. It has been discovered that the fabric anisotropy can play an important role in the occurrence of shear banding. An increasing degree of fabric anisotropy tends to delay the initiation of the strain localization and result in higher failure strength. The effects of fabric anisotropy have also been illustrated by comparing the theoretical predictions and measured results on the shear band inclination angle, shear strain level and dilatancy at bifurcation.     

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.     

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.

     

Discrete element modelling of material non‐coaxiality in simple shear flows

We investigate the quasi‐static simple shear flow of a two‐dimensional assembly of cohesionless particles using discrete element method (DEM) simulations. We focus on the unsteady flow regime where the solid would experience significant evolution of stresses, mobilised shear strength and dilation. We construct the DEM model using a discretised‐wall confined granular cell where the apparent boundary is allowed to dilate or contract synchronously with the confined solid. A rather uniform simple shear field is achieved across the whole assembly, which benefits rheological studies in generalising constitutive laws for continuum methods. We examine two aspects of the simple shear behaviour: macroscopic stress and strain rate evolution, particularly the non‐coaxiality between the principal directions of the two; and micromechanics such as evolution of fabric. For an initially anisotropic specimen sheared under constant normal pressure, the direction of principal stress rotates towards that of the principal strain rate, gradually reducing the degree of non‐coaxiality from about 45° to fluctuating around 0°. The rate in approaching coaxiality is slower in samples with larger initial porosity, stress ratio and mean stress. Generally, a faster rate in approaching coaxiality in simple shear is observed in a more dilatant sample, which often shows a larger degree of mobilised fabric anisotropy, suggesting the possible important role of instantaneous internal friction angle. The evolution of principal fabric direction resembles that of the principal stress direction. © 2013 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons, Ltd.     

复杂应力条件下土体各向异性及其建模思路

由于加荷方式不同,土体在复杂应力状态下在各主应力方向上应力-应变关系表现出显著应力各向异性,在常规三轴试验基础上,采用经典弹塑性理论各向同性土体模型对此不能合理描述。通过真三轴试验,总结应力各向异性柔度矩阵规律,结合试验规律进行相应理论研究,用非线性各向异性弹性矩阵代替弹塑性模型的弹性矩阵,用具有各向异性屈服准则的弹塑性模型描述塑性部分,建立非线性各向异性弹性-塑性模型,可以改善柔度矩阵矩阵形态,反映复杂应力状态下土体应力各向异性特征。     

Undrained expansion of a cylindrical cavity in clays with fabric anisotropy: theoretical solution

This paper presents a novel, exact, semi-analytical solution for the quasi-static undrained expansion of a cylindrical cavity in soft soils with fabric anisotropy. This is the first theoretical solution of the undrained expansion of a cylindrical cavity under plane strain conditions for soft soils with anisotropic behaviour of plastic nature. The solution is rigorously developed in detail, introducing a new stress invariant to deal with the soil fabric. The semi-analytical solution requires numerical evaluation of a system of six first-order ordinary differential equations. The results agree with finite element analyses and show the influence of anisotropic plastic behaviour. The effective stresses at critical state are constant, and they may be analytically related to the undrained shear strength. The initial vertical cross-anisotropy caused by soil deposition changes towards a radial cross-anisotropy after cavity expansion. The analysis of the stress paths shows that proper modelling of anisotropic plastic behaviour involves modelling not only the initial fabric anisotropy but also its evolution with plastic straining.     

土的可恢复剪胀的一种解释

从土的各向异怀角度对土的可恢复剪胀现象进行了解释。基于各向异性情况下的土体弹性本构关系理论分析，认为土的可恢复剪胀现象可部分归因于土的各向异性引起的弹性剪胀。借助有关土体弹性参数实验结果，研究了应力诱导各迥异性对土体弹性剪胀的影响，结果表明：随土体应力诱导各向异性的增大，土体的弹性剪胀也增大。从土体弹性剪胀角度研究了土的卸荷体缩条件，认为土体卸荷体缩取决于加载应力路径的应力增量比，给出了土体出现卸荷体缩的区域。     

宏-细观结合的砂土单剪试验非共轴特性分析

针对传统本构理论无法描述土体单剪试验非共轴变形的不足,采用非共轴修正模型进行改进。模型基于材料状态相关临界状态理论,采用宏-细观结合的方法,将1个新的各向异性状态变量引入本构模型来描述砂土的各向异性。考虑细观组构张量和应力张量的几何关系的变化,模型可以描述砂土在主应力轴旋转条件下材料状态的变化,材料状态变化直接导致模型的硬化规律和剪胀性发生变化,因此,模型可以描述该条件下原生向异性对砂土变形的影响。引入非共轴理论对本构模型进行修正,建立了三维非共轴各向异性模型。单剪试验的加载条件会造成主应力轴相对土体沉积面发生旋转,修正模型不但能够描述砂土在主应力轴旋转条件下其原生各向异性对变形的影响,而且可以描述主应力轴旋转造成的应力诱发各向异性对土体变形的影响,因此,该模型能够对整个单剪试验的变形规律进行描述,而且物理意义清晰。通过铝棒堆积体和Toyoura砂单剪试验验证表明,非共轴修正各向异性模型能对单剪试验的整个变形过程进行较好的模拟。     

A hypoplastic model for granular soils incorporating anisotropic critical state theory

It is well known that soil is inherently anisotropic and its mechanical behavior is significantly influenced by its fabric anisotropy. Hypoplasticity is increasingly being accepted in the constitutive modeling for soils, in which many salient features, such as nonlinear stress-strain relations, dilatancy, and critical state failure, can be described by a single tensorial equation. However, within the framework of hypoplasticity, modeling fabric anisotropy remains challenging, as the fabric and its evolution are often vaguely assumed without a sound basis. This paper presents a hypoplastic constitutive model for granular soils based on the newly developed anisotropic critical state theory, in which the conditions of fabric anisotropy are concurrently satisfied along with the traditional conditions at the critical state. A deviatoric fabric tensor is introduced into the Gudehus-Bauer hypoplastic model, and a scalar-valued anisotropic state variable signifying the interplay between the fabric and the stress state is used to characterize its impact on the dilatancy and strength of the soils. In addition, fabric evolution during shearing can explicitly be addressed. Modifications have also been undertaken to improve the performance of the undrained response of the model. The anisotropic hypoplastic model can simulate experimental tests for sand under various combinations of principle stress direction, intermediate principal stress (or mode of shearing), soil densities, and confining pressures, and the associated drastic effect of different principal stress orientations in reference to the material axes of anisotropy can be well captured.     

The signature of shear-induced anisotropy in granular media

This paper presents a micro-mechanical study on the characteristics of shear-induced anisotropy in granular media. Based on three-dimensional Discrete Element Method (DEM) simulations, the distinct features associated with the evolution of internal granular structure and different anisotropy sources during drained/undrained shearing of granular samples are carefully examined. The study finds that static liquefaction occurs when the geometrical anisotropy in a sample dominates the mechanical anisotropy in the overall shear strength, and the weak force network features an exceptionally high proportion of sliding contacts and develops certain degree of anisotropy. Phase transformation corresponds to a transitional, unstable state associated with a dramatic change in both coordination number and the proportion of sliding contacts in all contacts. The critical state in a granular material is always associated with a highly anisotropic fabric structure wherein both the critical void ratio and critical fabric anisotropy are uniquely related to the mean effective stress. The relations provide a more comprehensive definition for the critical state in granular media with proper reference to the critical fabric anisotropy.     

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.     

Three-dimensional random network model for thermal conductivity in particulate materials

This paper describes a three-dimensional random network model to evaluate the thermal conductivity of particulate materials. The model is applied to numerical assemblies of poly-dispersed spheres generated using the discrete element method (DEM). The grain size distribution of Ottawa 20–30 sand is modeled using a logistic function in the DEM assemblies to closely reproduce the gradation of physical specimens. The packing density and inter-particle contact areas controlled by confining stress are explored as variables to underscore the effects of micro- and macro-scales on the effective thermal conductivity in particulate materials. It is assumed that skeletal structure of 3D granular system consists of the web of particle bodies interconnected by thermal resistor at contacts. The inter-particle contact condition (e.g., the degree of particle separation or overlap) and the particle radii determine the thermal conductance between adjacent particles. The Gauss–Seidel method allows evaluation of the evolution of temperature variation in the linear system. Laboratory measurements of thermal conductivity of Ottawa 20–30 sand corroborate the calculated results using the proposed network model. The model is extended to explore the evolution of thermal conduction depending on the nucleation habits of secondary solid phase as an anomalous material in the pore space. The proposed network model highlights that the coordination number, packing density and the inter-particle contact condition are integrated together to dominate the heat transfer characteristics in particulate materials, and allows fundamental understanding of particle-scale mechanism in macro-scale manifestation.     

A thermomicromechanical approach to multiscale continuum modeling of dense granular materials

A new method is proposed for the development of a class of elastoplastic thermomicromechanical constitutive laws for granular materials. The method engenders physical transparency in the constitutive formulation of multiscale phenomena from the particle to bulk. We demonstrate this approach for dense, cohesionless granular media under quasi-static loading conditions. The resulting constitutive law—expressed solely in terms of particle scale properties—is the first of its kind. Micromechanical relations for the internal variables, tied to nonaffine deformation, and their evolution laws, are derived from a structural mechanical analysis of a particular mesoscopic event: confined, elastoplastic buckling of a force chain. It is shown that the constitutive law can reproduce the defining behavior of strain-softening under dilatation in both the mesoscopic and macroscopic scales, and reliably predict the formation and evolution of shear bands. The thickness and angle of the shear band, the distribution of particle rotation and the evolution of the normal contact force anisotropy inside the band, are consistent with those observed in discrete element simulations and physical experiments.     

基于微形态理论的各向异性散体波动分析

基于散粒体微观力学理论,忽略颗粒转动引起的相对位移,考虑颗粒接触的组构各向异性,根据宏微观能量守恒推导得到了散体材料各向异性微形态本构关系,进而通过单位接触方向积分的递推公式推导出了各向异性本构张量表达式;在此基础上,根据哈密顿原理得到了各向异性散体材料的运动平衡方程和边界条件,从而求得了平面波在各向异性散粒体中的传播规律和频散关系,最后对波的频散关系和频率带隙进行了详细的参数分析。研究表明：该模型预测了散体中包含3类12种位移波：3种纵波、6种横波和3种平面内横向剪切波;横观各向同性条件下,接触各向异性参数a20越大,纵波和横波的频率越大,而平面内横向剪切波的频率越小;正交各向异性条件下,随着接触各向异性参数a22的增大,与2方向运动相关的横波频率增大,而与3方向运动相关的横波频率则减小;但a22的变化对纵波频率影响很小。材料各向异性程度对横波带宽影响不大,但对纵波带宽影响较大：a20的增大使得声?光学波间的带宽减小,而光学波间的带宽增大,当a20>0.84时,声?光学波间的带隙消失;但是a22的增大则使得声?光学波间的带宽增大,而光学波间的带宽减小。退化为各向同性模型后,预测3类波的频散曲线与其他各向同性模型的结果基本一致。     

基于3D雕刻技术的岩体结构面剪切各向异性研究

相同形貌结构面重复性剪切试验和各向异性剪切试验一直是岩体结构面剪切力学特性试验研究中的难点,而该问题对于工程岩体的开挖力学响应和稳定性分析、评价与控制至关重要,问题的关键在于同一种岩石的结构面形貌的再现。为此,基于3D扫描和3D雕刻技术,重复制备了3种不同粗糙度的大理岩结构面试样,开展了不同法向应力下相同形貌结构面的各向异性剪切试验。试验结果发现：（1）同一法向应力、不同剪切方向下,相同形貌结构面的剪切强度、剪胀和剪切破坏特征均呈现明显各向异性;（2）结构面粗糙度各向异性很大程度上决定了剪切强度的各向异性,两者具有较好的正相关性;（3）随着法向应力的升高,结构面剪切特征的各向异性有逐渐弱化的趋势。同时,研究还充分表明,3D雕刻技术是系统开展结构面剪切力学特征各向异性研究的可靠手段,可在将来研究中发挥更为重要的作用。