On hypo-elastic analogues of the dilatant double-sliding model

The relation between the hypo-elastic constitutive law and Mehrabadi and Cowin's dilatant double-sliding model for cohesionless granular materials is studied. Conditions that must be satisfied by hypo-elastic analogues of the double-sliding model are derived constructively, and a simple example is given.     

FE-modeling of shear resistance degradation in granular materials during cyclic shearing under CNS condition

The behaviour of dry and cohesionless granular material during quasi-static cyclic shearing under a constant normal stiffness (CNS) condition is theoretically studied. A particular attention is laid to the volumetric strain change and the degradation of the shear resistance in the course of shearing. Numerical calculations are carried out for several shear cycles under boundary conditions which are relevant to investigate the shear interface behaviour. The global and local evolution of deformation, stress and density within the granular material is investigated with a finite element method on the basis of a hypoplastic constitutive model extended by micro-polar quantities: rotations, curvatures and couple stresses. A mean grain diameter is used as a characteristic length of micro-structure. The constitutive equations for stresses and couple stresses take also into account the effect of the evolution of the void ratio, pressure dependent relative density, direction of rate of deformation and rate of curvature. The numerical results are qualitatively compared with corresponding laboratory tests on direct wall shearing performed by DeJong, Randloph and White. In addition, the results for cyclic shearing of an infinite granular layer between two very rough boundaries under CNS conditions are also enclosed and discussed.     

A simple hypoplastic constitutive model for sand

The paper presents a hypoplastic constitutive model for the three-dimensional non-linear stress–strain and dilatant volume change behaviour of sand. The model is developed without recourse to the concept in elastoplasticity theory such as yield surface, plastic potential and decomposition into elastic and plastic parts. Benefited from the non-linear tensorial functions available from the representation theorem the model possesses simple mathematical formulation and contains only four material parameters, which can be easily identified with triaxial compression tests. Comparison of the predictions with the experimental results shows that the model is capable of capturing the salient behaviour of sand under monotonic loading and is applicable to both drained and undrained conditions.     

Modelling the time‐dependent behaviour of granular material with hypoplasticity

This paper presents a constitutive model for time‐dependent behaviour of granular material. The model consists of 2 parts representing the inviscid and viscous behaviour of granular materials. The inviscid part is a rate‐independent hypoplastic constitutive model. The viscous part is represented by a rheological model, which contains a high‐order term denoting the strain acceleration. The proposed model is validated by simulating some element tests on granular soils. Our model is able to model not only the non‐isotach behaviour but also the 3 creep stages, namely, primary, secondary, and tertiary creep, in a unified way.     

BEYOND FAILURE IN GRANULAR MATERIALS

Recent investigations on the hypoplastic constitutive model for granular materials show that the failure surface can be surpassed by some stress paths. This is contradictory to the conventional definition of failure surface in plasticity, according to which the stress is allowed to move on the failure surface but never across it. In the present paper, the interrelations among the different constitutive models are discussed with special reference to failure and stability. For the hypoplastic constitutive equation, the accessible stress states and the stable stress states are found to be enclosed by a bound surface and a stability surface in the stress space, respectively. Theoretical findings about the bound surface and the stability surface are verified qualitatively by presenting results of triaxial tests on dry sand. © 1997 by John Wiley & Sons, Ltd.     

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.     

堆石料Gudehus-Bauer亚塑性本构模型改进及参数确定方法

研究了Gudehus-Bauer亚塑性本构模型和模型参数的求取方法。采用侧限压缩试验曲线求取模型参数颗粒硬度hs和指数n。根据模型方程的推导,建立了拟合指数 和 与围压之间的关系,并提出了新的拟合参数。考虑到堆石料具有明显的剪胀、剪缩性,在Gudehus-Bauer模型线性项中增加了主要控制体积应变项 ,以改善模型对堆石料体积应变曲线的描述。采用堆石料大型侧限压缩试验、常规三轴试验分别验证了新的拟合参数和改进后的Gudehus-Bauer亚塑性本构模型。与堆石料试验成果比较,提出的新拟合参数与改进后的Gudehus-Bauer亚塑性本构模型可以较好地模拟堆石料的应力-应变特性,并较好地改善了堆石料体变曲线的模拟结果。对改进后的模型作了常规三轴加、卸载模拟,模拟结果反映了改进的Gudehus-Bauer亚塑性本构模型具有一定的卸载适应性。     

A hypoplastic constitutive model for debris materials

Debris flow is a very common and destructive natural hazard in mountainous regions. Pore water pressure is the major triggering factor in the initiation of debris flow. Excessive pore water pressure is also observed during the runout and deposition of debris flow. Debris materials are normally treated as solid particle–viscous fluid mixture in the constitutive modeling. A suitable constitutive model which can capture the solid-like and fluid-like behavior of solid–fluid mixture should have the capability to describe the developing of pore water pressure (or effective stresses) in the initiation stage and determine the residual effective stresses exactly. In this paper, a constitutive model of debris materials is developed based on a framework where a static portion for the frictional behavior and a dynamic portion for the viscous behavior are combined. The frictional behavior is described by a hypoplastic model with critical state for granular materials. The model performance is demonstrated by simulating undrained simple shear tests of saturated sand, which are particularly relevant for the initiation of debris flows. The partial and full liquefaction of saturated granular material under undrained condition is reproduced by the hypoplastic model. The viscous behavior is described by the tensor form of a modified Bagnold’s theory for solid–fluid suspension, in which the drag force of the interstitial fluid and the particle collisions are considered. The complete model by combining the static and dynamic parts is used to simulate two annular shear tests. The predicted residual strength in the quasi-static stage combined with the stresses in the flowing stage agrees well with the experimental data. The non-quadratic dependence between the stresses and the shear rate in the slow shear stage for the relatively dense specimens is captured.     

Finite deformation analysis of geomaterials

The mathematical structure and numerical analysis of classical small deformation elasto–plasticity is generally well established. However, development of large deformation elastic–plastic numerical formulation for dilatant, pressure sensitive material models is still a research area. In this paper we present development of the finite element formulation and implementation for large deformation, elastic–plastic analysis of geomaterials. Our developments are based on the multiplicative decomposition of the deformation gradient into elastic and plastic parts. A consistent linearization of the right deformation tensor together with the Newton method at the constitutive and global levels leads toward an efficient and robust numerical algorithm. The presented numerical formulation is capable of accurately modelling dilatant, pressure sensitive isotropic and anisotropic geomaterials subjected to large deformations. In particular, the formulation is capable of simulating the behaviour of geomaterials in which eigentriads of stress and strain do not coincide during the loading process. The algorithm is tested in conjunction with the novel hyperelasto–plastic model termed the B material model, which is a single surface (single yield surface, affine single ultimate surface and affine single potential surface) model for dilatant, pressure sensitive, hardening and softening geomaterials. It is specifically developed to model large deformation hyperelasto–plastic problems in geomechanics. We present an application of this formulation to numerical analysis of low confinement tests on cohesionless granular soil specimens recently performed in a SPACEHAB module aboard the Space Shuttle during the STS‐89 mission. We compare numerical modelling with test results and show the significance of added confinement by the thin hyperelastic latex membrane undergoing large stretching. Copyright © 2001 John Wiley & Sons, Ltd.     

Effect of grain crushing on shear localization in granular bodies during plane strain compression

This article deals with the effect of grain crushing on shear localization in granular materials during plane strain monotonic compression tests under constant lateral pressure. The grain diameter and the initial void ratio were stochastically distributed using a spatial correlation. To describe the mechanical behavior of cohesionless granular materials during a monotonic deformation path in plane strain compression, we used a micropolar hypoplastic constitutive model that is able to describe the salient properties of granular bodies including shear localization. The model was extended by introducing changes to the grain diameter with varying pressure using formulae from breakage mechanics proposed for crushable granulates. The initial void ratios and grain diameters took the form of correlated random spatial fields described by both symmetric and nonsymmetric random distributions using a homogeneous correlation function. The field realizations were generated with the help of an original conditional rejection method. A few representative samples of the random fields selected from the generated set were taken into account in numerical calculations. Copyright © 2011 John Wiley & Sons, Ltd.     

Non‐coaxiality and stress–dilatancy rule in granular materials: FE investigation within micro‐polar hypoplasticity

This paper deals with FE investigations of shear localization in dilatant granular bodies. The calculations were carried out with a hypoplastic constitutive law enhanced by micro‐polar terms to properly model the shear zone evolution. The behaviour of an initially medium dense sand specimen with very smooth and very rough horizontal boundaries was analyzed during a plane strain compression test. A stochastic distribution of the initial void ratio was assumed to be spatially correlated. Attention was focused on the non‐coaxiality of the directions of the principal strain increments and principal stresses in the shear zone and on the stress–dilatancy rule. Copyright © 2008 John Wiley & Sons, Ltd.     

Kinematic models for non‐coaxial granular materials. Part I: theory

The purpose of this paper is to present a physically based plasticity model for non‐coaxial granular materials. The model, which we shall call the double slip and rotation rate model (DSR² model), is a pair of kinematic equations governing the velocity field. The model is based on a discrete micro‐analysis of the kinematics of particles in contact, and is formulated by introducing a quantity called the averaged micro‐pure rotation rate (APR) into the unified plasticity model which was proposed by one of the authors. Our macro–micro mechanical analysis shows that the APR is a non‐linear function of, among other quantities, the macro‐rotation rate of the major principal axis of stress taken in the opposite sense. The requirement of energy dissipation used in the double‐sliding free‐rotating model appears to be unduly restrictive as a constitutive assumption in continuum models. In the DSR² model the APR tensor and the spin tensor are directly linked with non‐coaxiality of the stress and deformation rate tensors. We also propose a simplified plasticity model based on the DSR² model for a class of dilatant materials, and analyse its material stability. Copyright © 2005 John Wiley & Sons, Ltd.     

亚塑性理论简介

简要介绍了亚塑性理论和亚塑性本构模型以及在颗粒材料力学分析上的应用。详细地介绍了亚塑性模型参数的确定方法。亚塑性本构模型不是从弹塑性理论发展而来,它是以非线性张量为基础,没有屈服面、塑性势、流动法则、硬化定律以及把变形分解为弹性和塑性部分等概念。亚塑性在很多方面具有一定的优越性,特别是涉及到应力水平非常低或者非常高的情况。     

Unified modelling of granular media with Smoothed Particle Hydrodynamics

In this paper, we present a unified numerical framework for granular modelling. A constitutive model capable of describing both quasi-static and dynamic behaviours of granular material is developed. Two types of particle interactions controlling the mechanical responses, frictional contact and collision, are considered by a hypoplastic model and a Bagnold-type rheology relation, respectively. The model makes no use of concepts like yield stress or flow initiation criterion. A smooth transition between the solid-like and fluid-like behaviour is achieved. The Smoothed Particle Hydrodynamics method is employed as the unified numerical tool for both solid and fluid regimes. The numerical model is validated by simulating element tests under both quasi-static and flowing conditions. We further proceed to study three boundary value problems, i.e. collapse of a granular pile on a flat plane, and granular flows on an inclined plane and in a rotating drum.     

基于ABAQUS的颗粒材料亚塑性模型的数值实现及探讨

颗粒材料的亚塑性模型以Jaumann应力率张量及变形率张量描述本构关系,基于Cauchy应力的Jaumann速率及变形率给出了有限元的切线刚度矩阵,由此指出在严格意义上基于ABAQUS二次开发的亚塑性模型数值实现必需借助UEL接口。为了简化程序开发,文中建议了与有限元相结合的两种近似切线模量矩阵,即基于矩阵广义逆一致性切向模量矩阵与近似切向模量矩阵,从而形成了通过UMAT接口实现亚塑性模型的数值方案。由此降低了程序开发的难度,同时也可借助ABAQUS的后处理功能,提高了工作效率。数值算例表明了所开发程序的正确性以及所建议方案的可行性。     

粗粒料接触面模型及其在土石坝工程中的应用

在土石坝工程中越来越重视岸坡与坝料之间的接触特性,将土体本构模型中使用较多的双曲线型硬化规律运用至接触面模型,提出了一个简洁的接触面本构模型,推导了接触面模型刚度矩阵表达式。通过对4组试验的预测结果与试验结果的对比,表明该模型可以较好地预测粗粒料的界面剪切试验。将提出的接触面模型嵌入有限元程序,应用于如美心墙坝河谷与坝料的接触分析。三维有限元计算结果表明：坝体两侧岸坡的剪切位移变化规律都是中间部分剪切滑移量较大,岸坡边缘部分相对较小,陡坡一侧最大滑移量大于较缓一侧。这些均符合粗粒料滑移的基本规律,可为岸坡?坝料接触特性研究以及土石坝工程计算提供参考。     

一类新型的散粒型土体本构理论

详细地介绍了国际上近30才发展起来的具有独特建模思想的一类新型的散粒型土体本构建模理论--亚塑性本构理论。将该理论与经典的弹塑性理论进行比较后发现,亚塑性理论的本构表达式自动隐含了弹塑性理论中相对应的一些基本概念和假定,无需额外引入,因此增加了模型的客观性。以Gudehus-Bauer亚塑性模型为例,分析了模型在临界状态时的本构特征。此外,通过三轴试验和循环剪切试验的亚塑性数值模拟表明,亚塑性本构模型能很好地反映无黏性散粒型土体的非弹性、非线性及剪胀（剪缩）性等主要应力-应变特性。     

粗粒土的破碎耗能计算及影响因素

粗粒土的颗粒破碎直接改变了土体本身结构,对粗粒土的剪胀和内摩擦角都会产生影响。在土体剪切过程中,体积应力和剪切应力在体积应变和剪切应变上做功,这部分能量在剪切过程中转化为颗粒的弹性储能、颗粒间的摩擦耗能、颗粒剪胀时对外做功和颗粒破碎耗能4部分。准确计算剪切过程中粗粒土破碎耗能的目的是：从能量角度分析颗粒破碎对土体本构关系的影响,为建立考虑颗粒破碎的粗粒土本构关系创造条件。通过分析粗粒土的常规三轴试验数据,计算得到了剪切过程中的粗粒土破碎耗能。计算结果表明,常规三轴试验条件下粗粒土破碎耗能主要受固结应力、土体摩擦系数M等因素的影响。     

Deformation analysis of sands in cubical and hollow cylinder devices

For discrete materials like sands, the continuum field variables, stress and strain, are defined in terms of micro-level quantities by considering the deformation mechanism of granular soils from a microscopic point of view. Under the application of load, soil is considered to deform due to the movement relative to each other of clusters of particles. Based on this deformation mechanism, the kinematics of soils are developed and a strain tensor for granular soils, in terms of local displacements and geometric measures, is introduced. A local constitutive law relating local displacements and local tractions is defined. Using the local constitutive law, the relationships between stress and strain for the media are developed. The developed model incorporates the influence of strain hardening and material anisotropy on the deformation behaviour of the media. Comparisons of the model predictions and experimental results from tests conducted in cubical and hollow cylinder devices are presented.     

Numerical investigations of shear localization in a micro‐polar hypoplastic material

In this paper a micro‐polar continuum approach is proposed to model the essential properties of cohesionless granular materials like sand. The model takes into account the influence of particle rotations, the mean grain size, the void ratio, the stresses and couple stresses. The constitutive equations for the stresses and couple stresses are incrementally non‐linear and based on the concept of hypoplasticity. For plane strain problems the implementation of the model in a finite element program is described. Numerical studies of the evolution of micro‐polar effects within a granular strip under plane shearing are presented. It is shown that the location and evolution of shear localization is strongly influenced by the initial state and the micro‐polar boundary conditions. For large shearing the state quantities tend towards a stationary state for which a certain coupling between the norm of the stress deviator and the norm of the couple stress tensor can be derived. Copyright © 2003 John Wiley & Sons, Ltd.