An Anisotropic Model for Granular Material Based on Experiments

Soil is a heterogeneous material and most natural soil deposits show a definite stratification. The mechanical behaviour of such material is generally different in different directions, especially in the direction parallel and perpendicular to the stratification. A series of isotropic compression tests were carried out to study the behavior of granular material produced under controlled stratification in the laboratory. These tests were conducted both on cylindrical and square prismatic tri-axial specimens. It was observed that for hydrostatic loading, the strain response was different in different directions, especially in directions parallel and perpendicular to the direction of soil deposition. A definite trend of anisotropy was observed in the deformation pattern. The observed anisotropy is modeled in this paper by treating soil-dilatancy as a variable quantity. The equation of the plastic potential surface of the model which obeys a non-associated flow rule, is assumed to be dependent on three main variables confining pressure (\(\sigma_{3}\)), void ratio (e) and the angle of bedding plane orientation (δ) during deposition. The angle of bedding plane orientation (δ) was measured with respect to the direction of the major principal stress. The model has a cap yield surface in the isotropic stress direction, which is supplemented by a shear hardening Mohr–Coulomb surface in the deviator direction. This paper focuses on predicting the anisotropic strain response of stratified soil deposits subjected to isotropic compression. The proposed anisotropic model incorporates within an existing strain-hardening sand model, a modified cap yield surface and a modified plastic potential function related to the cap surface, to account for the anistropic response observed in isotropic compression tests. The two dimensional stress–strain model was extended to three dimensional Cartesian space. The strain anisotropy observed in the isotropic compression tests was predicted by the three dimensional anisotropic model proposed for granular materials.     

A note on thermomechanical anisotropy of clays

Deep clays exhibit a pronounced strain anisotropy both during mechanical loading as well as during heating and cooling at constant stress in drained isotropic conditions. During mechanical loading vertical strain is larger than the horizontal one. During heating the vertical strain is larger than the horizontal one within the elastic range; the opposite is observed in the elasto-plastic range. The above described response can be interpreted adopting a consistent rotational, kinematic hardening thermo-elasto-plastic constitutive law.     

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.     

基于旋转运动硬化理论的修正剑桥模型的改进

修正剑桥模型是最早建立和得到广泛承认的经典土体弹塑性模型之一,但不能模拟应力路径转折时土体的应力-应变特性以及应力引起的各向异性。将旋转运动硬化理论引入到剑桥模型中,给出了椭圆屈服面的旋转运动硬化机制,在不增加任何模型参数的情况下,把等向硬化的修正剑桥模型扩展为旋转运动硬化模型。扩展的新模型既保留了单调加载时的等向硬化的特性,又能反映应力路径转折时土体的本构特性与应力诱发的各向异性,初步验证了模型的有效性。     

Numerical modeling of changes in anisotropy during liquefaction using a generalized constitutive model

A constitutive model with rotation hardening was generalized from the triaxial compression state to a general stress state. With the generalized model, numerical simulations of cyclic and monotonic undrained triaxial tests were conducted to reproduce the phenomenon of continuous, orderly and rapid changes in anisotropy during liquefaction. The simulated results demonstrated that when sand enters the liquefaction stage, the yield surface in the stress space rotates quickly, causing continuous and rapid changes in anisotropy. Through comparison of the simulated and experimental results, the generalized constitutive model was shown to be able to capture the fundamental behaviors of sand demonstrated by the experimental data, and the rotational hardening rule adopted in the generalized model was proven suitable for describing the continuous, orderly and rapid changes in anisotropy that occur during liquefaction.     

Predictions of large stress reversals in true triaxial tests on cross‐anisotropic sand

A rotational kinematic hardening constitutive model with the capability of predicting the behavior of soil during three‐dimensional stress reversals has been developed. An existing elasto‐plastic constitutive model, the Single Hardening Model, utilizing isotropic hardening serves as the basic framework in these formulations. The framework of the kinematic hardening model was discussed in a companion paper. The previously proposed cross‐anisotropic Single Hardening Model is added to the present kinematic hardening mechanism to capture inherent anisotropy of sands in addition to the stress reversals. This model involves 13 parameters, which can be determined from simple laboratory experiments, such as isotropic compression, drained triaxial compression and triaxial extension tests. The results from a series of true triaxial tests with large three‐dimensional stress reversals performed on medium dense cross‐anisotropic Santa Monica Beach sand are employed for comparison with predictions. Copyright © 2008 John Wiley & Sons, Ltd.     

Localized failure of fibre‐reinforced elastic–plastic materials subjected to plane strain loading

We consider discontinuous bifurcations as the indicator of a localized failure for a class of composites that are characterized by elastic fibres reinforcing an elastic–plastic matrix. A macroscopic tangent stiffness tensor for the fibre‐reinforced composite is developed by consistently homogenizing the contribution of fibres in a spherical representative volume element. Analytical solutions are derived for the critical hardening modulus and corresponding bifurcation directions for the case of plane strain loading. Properties of the solutions are further illustrated on the example of the non‐associated Drucker–Prager model at onset of yielding. Results show that presence of fibres decreases the critical hardening modulus, thus inhibiting the onset of strain localization. The rate of decrease in the critical hardening modulus is the highest for pure shear, followed by uniaxial tension, uniaxial compression, biaxial tension and biaxial compression. The main fibre parameters that control the onset of strain localization are their volumetric content and their stiffness modulus whereby very stiff fibres can produce the most significant decrease in the critical hardening modulus, especially for the state of biaxial tension. The critical hardening modulus for the non‐associated Drucker–Prager model exhibits a full range of localization modes including compaction bands, dilation bands, and transition in the form of shear bands regardless of the presence of fibres. Presence of fibres affects bifurcation directions, except in the case when Poisson's ratio of the matrix is equal to 0.25. The results demonstrate stabilizing effects of fibres by which they provide the control against the onset of strain localization. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulation of yielding and stress–stain behavior of shanghai soft clay

In this paper, a simple bounding surface plasticity model is used to reproduce the yielding and stress–strain behavior of the structured soft clay found at Shanghai of China. A series of undrained triaxial tests and drained stress probe tests under isotropic and anisotropic consolidation modes were performed on undisturbed samples of Shanghai soft clay to study the yielding characteristics. The degradation of the clay structure is modeled with an internal variable that allows the size of the bounding surface to decay with accumulated plastic strain. An anisotropic tensor and rotational hardening law are introduced to reflect the initial anisotropy and the evolution of anisotropy. Combined with the isotropic hardening rule, the rotational hardening rule and the degradation law are incorporated into the bounding surface formulation with an associated flow rule. Validity of the model is verified by the undrained isotropic and anisotropic triaxial test and drained stress probe test results for Shanghai soft clay. The effects of stress anisotropy and loss of structure are well captured by the model.     

Use of acoustic emissions to investigate localization in high-porosity sandstone subjected to true triaxial stresses

Acoustic emissions were used to investigate the evolution of damage and strain localization in Castlegate sandstone specimens subjected to true triaxial states of stress, where the intermediate principal stress ranged from equal to minimum compression to equal to maximum compression. Three failure modes were observed: shear band formation at low mean stresses, compaction localization at intermediate stresses, and no localization at high mean stresses. For shear bands, the onset of localization consistently occurred prior to peak stress, while compaction localization initiated at the beginning of a stress plateau. The band angle (defined as the angle between the band normal and the maximum compression direction) determined by fitting a plane through the localized acoustic emission events corresponded well with the shear band angle expressed on the specimen jacket. As expected, the band angle decreased with increasing mean stress. The theoretically predicted dependence of band angle on deviatoric stress state was not confirmed; however, data scatter due to natural variations in material could obscure such a trend. Each failure mode displayed a unique acoustic emission rate response; therefore, this rate response alone can be used to determine the failure mode and the onset of localization in Castlegate sandstone.     

不同应力路径下剪切带的数值模拟

采用回映应力更新算法,编写了基于伏斯列夫面的超固结黏土本构关系模型子程序,嵌入非线性有限元软件ABAQUS。通过对单元试验进行三轴压缩、三轴伸长及平面应变等问题的模型预测,再现了超固结黏土在不同初始超固结比和应力路径时的变形和强度特性,从而验证了子程序的正确性。借助该本构模型,对三轴压缩、三轴伸长及平面应变应力路径下超固结黏土体变形局部化问题,进行了三维数值模拟。分析结果表明：超固结黏土在三轴压缩及伸长状态时,土体变形局部化在应力-应变关系软化时出现,而平面应变状态时,在应力-应变关系硬化阶段出现,其超固结黏土的剪胀特性在剪切带的形成过程中起重要作用。     

Kinematic hardening predictions of large stress-reversals in 3-D test on loose sand

A rotational kinematic hardening constitutive model with the capability of predicting the behavior of soil during 3-D stress-reversals has been developed. An existing elasto-plastic constitutive model, the Single Hardening Model, utilizing isotropic hardening serves as the basic framework in these formulations. To this framework is added the capability of handling cross-anisotropic behavior as well as the kinematic hardening mechanism to capture inherent anisotropy of the sand in addition to the large stress-reversals. The model involves thirteen parameters, which can be determined from simple laboratory experiments, such as isotropic compression, drained triaxial compression and triaxial extension tests. The results from a series of true triaxial tests with large stress-reversals performed on loose cross-anisotropic Santa Monica Beach sand are employed for comparison with predictions.     

A model of strain localization in porous sandstone as a function of tectonic setting,burial and material properties; new insight from Provence (southern France)

The analysis of three cataclastic band sets from Provence (France) reveals that the band density, their conjugate angles, their ratio of shear displacement to compaction, and the amount of cataclasis within the bands differ and can be expressed as functions of tectonic setting and petrophysical properties. We identify (1) a dense and closely spaced network of shear-enhanced (reverse) compaction bands; (2) a regularly spaced less dense network of reverse compactional shear bands; and (3) a localized network of normal shear bands. The field data show that strain localization is favored in an extensional regime and is characterized by shear bands with a large shear to compaction ratio and a small conjugate band angle. In contrast, distributed strain is favored in a contractional regime and is characterized by compactional bands with a low ratio of shear to compaction and a large conjugate band angle. To explain the mechanical origin of this strain localization, we quantify the yield strength and the stress evolution in extensional and contractional regimes in a frictional porous granular material. We propose a model of strain localization in porous sands as a function of tectonic stresses, burial depth, material properties, strain hardening and fluid pressure. Our model suggests that stress reduction, inherent to extensional regime, favors strain localization as shear bands, whereas stress increase during contraction favors development of compactional bands.     

干密度对粗粒料力学特性的影响

粗粒料的力学特性不仅取决于应力状态,还与粗粒料本身的松密程度密切相关,即与粗粒料的材料状态相关。通过不同初始干密度的粗粒料常规大型三轴各向等压固结排水剪切试验研究了干密度对粗粒料力学特性的影响。试验结果表明：粗粒料应力-应变曲线的形态取决于密度和围压的共同作用,破坏状态之前密度和围压共同决定了粗粒料的强度和变形,而当应变足够大处于渐进状态或临界状态时,粗粒料的应力和体积应变受初始干密度的影响逐渐减小直至消失。对于软化型曲线,相变状态的应力小于渐进状态的应力,而渐进状态的应力小于破坏状态的应力。无论围压大小如何,随着初始干密度的增大,应力-应变曲线的硬化性逐渐减弱,而软化性则逐渐增强。     

On the undrained strain-induced anisotropy of loose sand

An experimental study was carried out to investigate the effects of previous deviatoric strain histories on the undrained behaviour of loose and saturated Toyoura sand and compared with known results of Hostun RF sand. From an initial isotropic stress state, recent deviatoric strain histories in the compression side of the triaxial plane were generated by a standard drained presheared cycle up to a specified mobilized stress ratio. Mainly, the fully liquefied, contractive, unstable and softening behaviour of loose sand was progressively transformed into the non-liquefied, dilative, fully stable and hardening behaviour of dense-like sand, while remaining within a narrow range of loose density. The paper validates and extends the current understanding of strain-induced anisotropy of loose sand. New experimental data support the directional dependency of the instability cone on the stress increment direction, suggest the bifurcation characteristics of loose sand and evidence the important role of past deviatoric strain histories.     

Thermoplasticity and strain localization in transversely isotropic materials based on anisotropic critical state plasticity

Geomaterials such as soils and rocks are inherently anisotropic and sensitive to temperature changes caused by various internal and external processes. They are also susceptible to strain localization in the form of shear bands when subjected to critical loads. We present a thermoplastic framework for modeling coupled thermomechanical response and for predicting the inception of a shear band in a transversely isotropic material using the general framework of critical state plasticity and the specific framework of an anisotropic modified Cam–Clay model. The formulation incorporates anisotropy in both elastic and plastic responses under the assumption of infinitesimal deformation. The model is first calibrated using experimental data from triaxial tests to demonstrate its capability in capturing anisotropy in the mechanical response. Subsequently, stress‐point simulations of strain localization are carried out under two different conditions, namely, isothermal localization and adiabatic localization. The adiabatic formulation investigates the effect of temperature on localization via thermomechanical coupling. Numerical simulations are presented to demonstrate the important role of anisotropy, hardening, and thermal softening on strain localization inception and orientation. Copyright © 2016 John Wiley & Sons, Ltd.     

煤岩体孔隙结构应力特征的数值模拟研究

为了探究不同加载方式下煤岩体孔隙周围的应力分布规律,通过CT三维重建技术构建含有不同孔隙形状的煤岩体骨架模型,并利用ABAQUS软件进行单轴、三轴压缩实验模拟。结果表明,球状孔隙结构在单轴压缩条件下,上下区域表现为拉应力集中,左右区域表现为压应力集中。不同倾角的椭球状孔隙结构其长短轴区域的应力集中类型不同。单轴压缩过程中,轴向加载速度影响球状孔隙周围的Mises应力峰值和σ₁应力的变化;三轴压缩过程中,孔隙结构依然经历了压密、弹性、塑性和破坏4个阶段,较低的围压条件使得弹性阶段“应力-应变”曲线与“应力-时间”曲线高度重合。从微观角度为煤岩体力学研究提供了一种新的方法和思路。     

Multiscale modeling and analysis of compaction bands in high-porosity sandstones

We present a multiscale investigation on the initiation and development of compaction bands in high-porosity sandstones based on an innovative hierarchical multiscale approach. This approach couples the finite element method and the discrete element method (DEM) to offer direct, rigorous linking of the microscopic origins and mechanisms with complex macroscopic phenomena observed in granular rocks such as strain localization and failure. To simulate compaction band in granular cementitious sandstone, we adopt a bonded contact model with normal and tangential interparticle cohesions in the DEM and propose a dual-porosity structure consisting of macro-pores and interstitial voids for the representative volume element to mimic the typical meso-structure of high-porosity sandstones. In the absence of particle crushing, our multiscale analyses identify debonding and pore collapses as two major contributors to the formation of compaction bands. The critical pressures predicted by our simulations, corresponding to surges of debonding and pore collapse events, agree well with the estimations from field data. The occurrence patterns of compaction band are found closely related to specimen heterogeneity, porosity and confining pressure. Other deformation band patterns, including shear-enhanced compaction bands and compactive shear bands, were also observed under relatively low confining pressure conditions with a rough threshold at \(0.55P^{*}\) (\(P^{*}\) is the critical pressure) on the failure envelop. Key microscopic characteristics attributable to the occurrence of these various deformation patterns, including fabric anisotropy, particle rotation, debonding and pore collapse, are examined. Shear-enhanced compaction bands and pure compaction bands bear many similarities in terms of these microscopic characteristics, whereas both differ substantially from compactive shear bands.     

各向异性弹黏塑性模型在ABAQUS中的研发

通过将Perzyna过应力理论与临界状态理论相结合,并引入Wheeler旋转硬化法则,提出一个能描述土体初始各向异性及应力诱发各向异性的三维弹黏塑性本构模型。模型考虑流变发生的下限,在三维应力空间,模型存在形状相似的静屈服面及动态加载面。采用缩放形式的幂函数。本构模型数值算法采用回映算法,借助ABAQUS软件UMAT子程序接口实现。并通过对三轴不排水蠕变试验的模拟,确定合适的积分步长。此后,分别对三轴不排水蠕变试验及常应变率三轴不排水剪切试验进行了模拟。模拟中通过设置不同参数值,可将模型退化为各向同性模型,并对这两种模拟结果进行了比较。模拟结果表明：(1) 对于三轴不排水蠕变,在低剪应力水平下,各向同性模型和各向异性模型模拟的结果相差不大,而在高剪应力水平下,各向异性模型模拟结果更接近试验结果;(2) 对于常应变率加载试验的模拟,模型合理反映了土体不排水强度随着加载速率的增大而增大现象。     

Induced anisotropy in cohesionless soil: Experiments and modelling

Induced anisotropy in cohesionless material is an up to date topic: since 1980, some papers present experimental results which clearly demonstrate that the effects of induced anisotropy are quite pronounced, but, at this time few authors have checked the predictions of their constitutive relations in such experiments.

First, this paper describes experimental results performed on a true triaxial apparatus which allows 90°-rotations of principal stress directions. Therefore induced anisotropy created by an initial loading may be studied. Then, an elastoplastic model with isotropic and kimematic hardening is developed. Its predictions for anisotropic effects are in good agreement with experiments.     

颗粒物质在等比例应变加载下的分散性失稳模式

由地下水引起的静力液化可能是边坡失稳的隐含机制之一,松砂在不排水剪切条件下可能发生静力液化,密实的颗粒集合体在特定的应变路径下也会出现相似的现象,即试样整体发生急剧的失稳,应力状态尚处于峰值强度线以内。该种失稳模式称为分散性失稳,是为了强调失稳模式中没有出现应变局部化或者剪切带。采用连续-离散耦合分析方法,研究由不规则形状颗粒组成的密实集合体在等比例应变加载路径下的力学特性。根据Hill的材料失稳理论,当试样的应力增量 和应变增量 对应的2阶功 为负时,试样即发生不可逆的整体失稳破坏。以根据不同等比例应变路径得到 曲线为界,在 平面内将试样的应力状态分为剪缩区、剪胀-稳定区和剪胀-非稳定区,连接不同围压下试样发生分散性失稳时的应力状态形成失稳线发现,峰值强度线高于临界状态线,临界状态线高于失稳线。