Gravity-driven dry granular slow flows down an inclined moving plane: a comparative study between two concepts of the evolution of porosity

Two concepts in modeling the effects of the evolution of porosity in dry granular flows are investigated to illuminate their performance and limitations. To this end, the thermodynamic analysis, based on the Müller-Liu entropy principle, and the quasi-linear theory, are employed to deduce the ultimate constitutive models and the restrictions on their thermodynamic consistencies. The models are employed to study an isothermal dry granular slow flow down an inclined moving plane, of which the results are compared with the experimental outcomes. Results show that, while the two models deliver appropriate equilibrium expressions of the Cauchy stress tensor for compressible grains, the model in which the evolution of porosity is treated kinematically yields a spherical stress tensor for incompressible grains. Only the model with a dynamic evolution of porosity can give rise to a non-spherical stress tensor at equilibrium. Moreover, whilst the former model can better capture the characteristics of flows with slow to moderate speeds, the latter model is more able to describe the features of very rapid flows like avalanches. The present study illustrates the essential difference between the two concepts in modeling the effects of the evolution of porosity, and can be extended for further studies on other microstructural effects in granular flows.     

Constitutive models for granular materials including quasi-static frictional behaviour: Toward a thermodynamic theory of plasticity

This work deals with the thermodynamic formulation of constitutive models for materials whose quasi-static behaviour is governed by internal friction, e.g., dry granular materials. The process of internal friction is represented here phenomenologically with the help of a second-order, symmetric-tensor-valued internal variable. A general class of models for the evolution of this variable is considered, including as special cases a hypoelastic-like form for this relation as well as the hypoplastic form of Kolymbas (1991). The thermodynamic formulation is carried out in the context of the Müller-Liu entropy principle. Among other things, it is shown that for the hypoelastic-type models, a true equilibrium inelastic Cauchy stress exists. On the other hand, such a stress does not exist for the hypoplastic model due to its rate-independence and incremental non-linearity. With the help of a slight generalization of the notion of thermodynamic equilibrium, i.e., to thermodynamic “quasi-equilibrium,” however, such a Cauchy stress can be formulated for the hypoplastic model. As it turns out, this quasi-equilibrium for the Cauchy stress represents a thermodynamic generalization of the so-called quasi-static stress postulated for example by Goddard (1986) in the context of his viscoplastic model for a frictional-dissipative, and in particular for granular, materials. Received January 22, 1999     

A generalized Mohr–Coulomb criterion implied by the revised Goodman–Cowin theory

In the present study, the constitutive relations derived in the revised Goodman–Cowin theory for granular matter are shown to imply a generalized Mohr–Coulomb criterion for impending flows. Due to the concept of microcontinuum and the incorporation of the internal friction into the expression of the Cauchy stress tensor, a constrained equilibrium stress state characterized by the Mohr–Coulomb criterion is yet obtained under uniform distributions of the grains.     

A generalized continuum theory for granular materials

A generalized continuum theory for granular media is formulated by allowing for the possibility of rotation of granules. The basic balance laws are presented and based on thermodynamical consideration a set of constitutive equations are derived. The theory naturally gives rise to the generation of antisymmetric stress tensor and existence of couple stresses. The basic equations of motion are derived and it is shown that the theory contains Mohr-Coulomb criterion of limiting equilibrium as a special case. The problem of coupled porosity and microrotational wave propagation is investigated and the rectilinear shear flow of granular materials is discussed.     

基于深度强化学习算法的颗粒材料应力?应变关系数据驱动模拟研究

颗粒材料的宏观力学行为受颗粒组分等材料参数, 孔隙率、配位数等状态参数的影响, 同时又具备复杂的加载路径和加载历史相关性, 建立包含多个内变量以及各变量间相互关联的颗粒材料本构模型是一个重要的科学难题. 不同于传统的基于屈服面、流动法则和硬化函数框架下的唯象本构模型, 本文基于颗粒物质力学的研究基础, 以颗粒材料平均孔隙率、细观组构参数和弹性刚度参数作为内变量, 结合深度学习方法建立以有向图表征的数据本构模型. 有向图中以不同的链接网络表示不同的内变量信息流动方向, 各个内变量间的映射关系采用循环神经网络来建立, 将各个神经网络相互组合, 形成包含不同内变量且具有不同预测能力的本构模型. 该本构模型的建立过程等价于在众多可能的内变量链接关系空间中寻找最能描述实际材料宏观应力应变行为的优化问题. 因此, 可将有向图本构模型的建立过程看作“马尔可夫决策过程”, 采用深度强化学习算法构建有向图的内变量链接组合优化过程, 具体采用AlphaGo Zero算法自动寻找最优的颗粒材料数据驱动本构模型建模路径. 研究结果表明, 采用有向图和深度强化学习算法可建立起完全依靠“数据驱动”的颗粒材料应力?应变关系. 此外, 本方法提供了一种将不同理论模型从数据角度统一起来, 且基于人工智能算法发展更优模型的研究思路, 可为相似问题的研究提供借鉴.      

基于介观结构的饱和与非饱和多孔介质有效应力

基于描述含液颗粒材料介观结构的Voronoi 胞元模型和离散颗粒集合体与多孔连续体间的介-宏观均匀化过程, 定义饱和与非饱和多孔介质有效应力. 导出了计及孔隙液压引起之颗粒体积变形的饱和多孔介质广义有效应力. 用以定义广义有效应力的Biot 系数不仅依赖于颗粒材料的多孔连续体固体骨架及单个固体颗粒的体积模量(材料参数),同时与固体骨架当前平均广义有效应力及单个固体颗粒的体积应变(状态量) 有关. 提出了描述非饱和多孔介质中非混和固体颗粒、孔隙液体和气体等三相相互作用的具介观结构的Voronoi 胞元模型.具体考虑在低饱和度下双联(binary bond) 模式的摆动(pendular) 液桥系统介观结构. 导出了基于介观水力-力学模型的非饱和多孔介质的各向异性有效应力张量与有效压力张量. 考虑非饱和多孔介质Voronoi 胞元模型介观结构的各向同性情况,得到了与非饱和多孔连续体理论中唯象地假定的标量有效压力相同的有效压力形式.但本文定义的与确定非饱和多孔介质有效应力和有效压力相关联的Bishop 参数由基于三相介观水力-力学模型, 作为饱和度、孔隙度和介观结构参数的函数导出,而非唯象假定.      

Towards a constitutive law for the unsteady contact stress in granular media

In order to build a unified modelling for granular media by means of Eulerian averaged equations, it is necessary to study two contributions in the effective Cauchy stress tensor: the first one concerns solid and fluid matter, including contact and collisions between grains; the second one focuses on the random movements of grains and fluid, similar to Reynolds stress for turbulent flows. It is shown that the point of view of piecewise continuous media already used for two phase flows allows one to derive a constitutive equation for the first contribution, under the form of a partial differential equation. Similarly as for the Reynolds stress in turbulent flows, this equation cannot be written only in terms of averaged quantities without adequate approximations. The structure of the closed equation is discussed with respect to irreversible thermodynamics, and in connection with some already existing models. It is emphasised that numerical simulations by the discrete elements method can be used in order to validate these approximations, through numerical experiments statistically considered. Finally an extension of this approach to the second contribution of the effective Cauchy stress tensor is briefly discussed, showing how the modelling of both contributions have to be coupled.      

A thermo-mechanical continuum theory with internal length for cohesionless granular materials

This article continues Part I. Here the non-equilibrium responses of the constitutive variables t (Cauchy stress tensor), q (heat flux vector), h (equilibrated stress vector), Γ (flux term associated with the internal length ℓ), Π (production term associated with ℓ) and f (equilibrated intrinsic body force) as well as the Helmholtz free energy Ψ are postulated by use of a quasi-linear theory for three of four models deduced in Part I. In so doing, together with the equilibrium responses gained in Part I, a complete set of constitutive equations for the constitutive quantities for each model is obtained. The implemented models are applied to investigate typical isothermal steady granular shearing flows with incompressible grains, namely, simple plane shear flow, inclined gravity-driven flow and vertical channel-flow. The emphasis is on the models in which ℓ is considered a material constant (Model I) and an independent dynamic field quantity (Model III). Numerical results show that Model III is more appropriate than Model I since in the former model the effect of the motion of an individual grain can better be taken into account. Such a result is in particular significant for avalanches, since it verifies the existence of a thin layer immediately above the base of an avalanche, in which the grains are colliding strongly with one another, and provides a quantitative means to measure such a thin layer.     

Some Properties of a New Model for Slow Flow of Granular Materials

Some properties of a new continuum model for the bulk flow of a dense granular material in which neighbouring grains are in contact for a finite duration of time and in which the contact force is non-impulsive – the so called slow flow regime – are presented. The model generalises both the plastic potential and double-shearing models and contains an additional kinematic quantity – the intrinsic spin. The stress tensor is, in general, non-symmetric and separate yield conditions govern translational and rotational yield. We consider homogeneous, quasi-static loadings for the symmetric part of the stress and dynamic loading for the anti-symmetric part of the stress. A solution for the stress state in terms of a single parameter, namely the major principal direction of the symmetric part of the stress, is presented. This direction itself is determined by a consideration of the flow equations in the context both dilatant and isochoric simple shear flows. These simple flows are used to complete the characterisation of the relationship between the anti-symmetric part of the stress and the intrinsic spin.     

Multi-slip gradient formulation for modeling microstructure effects on shear bands in granular materials

This paper presents a higher order gradient multi-slip formulation to model the effect of inhomogeneous deformation in granular materials. The effects of heterogeneity and porosity anisotropy within the multi-slip formulation are taken into consideration through the modification of the mobilized friction. The mobilized friction is assumed to be a direct function of either the gradient of the porosity distribution or the fabric tensor. The formulation with two active slip planes was implemented into a finite element code and used to simulate biaxial shear tests on dry sand. The analysis quantifies most of the shear band characteristics observed by past experimentation. It is shown that the localization and shear band characteristics in granular materials are very much dependent on the initial fabric and slip system arrangement.     

A unified evolution equation for the Cauchy stress tensor of an isotropic elasto-visco-plastic material

In the present study an evolution equation for the Cauchy stress tensor is proposed for an isotropic elasto-visco-plastic continuum. The proposed stress model takes effects of elasticity, viscosity and plasticity of the material simultaneously into account. It is ascribed with some scalar coefficient functions and, in particular, with an unspecified tensor-valued function N, which is handled as an independent constitutive quantity. It is demonstrated that by varying the values and the specific functional forms of these coefficients and N, different known models in non-Newtonian rheology can be reproduced. A thermodynamic analysis, based on the Müller–Liu entropy principle, is performed. The results show that these coefficients and N are not allowed to vary arbitrarily, but should satisfy certain restrictions. Simple postulates are made to further simplify the deduced general results of the thermodynamic analysis. They yield justification and thermodynamic consistency of the existing models for a class of materials embracing thermoelasticity, hypoelasticity and in particular hypoplasticity, of which the thermodynamic foundation is established successively for the first time in literature. The study points at the wide applicability and practical usefulness of the present model in different fields from non-Newtonian fluid to solid mechanics. In this paper the thermodynamic analysis of the proposed evolution-type stress model is discussed, its applications are reported later.      

A model for elastic-plastic pressure sensitive materials subjected to large deformations

Development of a finite deformation elasto-plastic model for pressure sensitive materials is presented. The chosen model, which has its roots in the MRS-Lade material model is influenced by recent developments. The thermodynamic consequences of introducing non-associative yielding (both deviatoric and volumetric) and hardening⧸softening characteristics are assessed. The consistently linearized Algorithmic Tangent Stiffness (ATS) tensor is presented. This tensor is used in the constitutive driver as a key feature of the efficient iterative procedure for satisfying equilibrium in the case of stress (or mixed) control.The chosen model is calibrated using data from experiments conducted in a Directional Shear Cell (DSC) , which has been used extensively at the University of Colorado at Boulderto investigate the behavior of pressure sensitive materials under deformations of large magnitude.     

Statics and kinematics of discrete Cosserat-type granular materials

A theoretical framework is presented for the statics and kinematics of discrete Cosserat-type granular materials. In analogy to the force and moment equilibrium equations for particles, compatibility equations for closed loops are formulated in the two-dimensional case for relative displacements and relative rotations at contacts. By taking moments of the equilibrium equations, micromechanical expressions are obtained for the static quantities average Cauchy stress tensor and average couple stress tensor. In analogy, by taking moments of the compatibility equations, micromechanical expressions are obtained for the (infinitesimal) kinematic quantities average rotation gradient tensor and average Cosserat strain tensor in the two-dimensional case. Alternatively, these expressions for the average Cauchy stress tensor and the average couple stress tensor are obtained from considerations of the equivalence of the continuum force and couple traction vectors acting on a plane and the resultant of the discrete forces and couples acting on this plane. In analogy, the expressions for the average rotation gradient tensor and the average Cosserat strain tensor are obtained from considerations of the change of length and change of rotation of a line element in the two-dimensional case. It is shown that the average particle stress tensor is always symmetrical, contrary to the average stress tensor of an equivalent homogenized continuum. Finally, discrete analogues of the virtual work and complementary virtual work principles from continuum mechanics are derived.     

Assessment of plastic heterogeneity in grain interaction models using crystal plasticity finite element method

Micromechanical models aimed at simulating deformation textures and resulting plastic anisotropy need to incorporate local plastic strain heterogeneities arising from grain interactions for better predictions. The ALAMEL model [Van Houtte, P., Li, S., Seefeldt, M., Delannay, L. 2005. Deformation texture prediction: from the Taylor model to the advanced Lamel model. Int. J. Plasticity 21, 589–624], is one of the models in which the heterogeneous nature of plastic deformation in metals is introduced by accounting for the influence of a grain boundary on the cooperative deformation of adjacent grains. This is achieved by assuming that neighbouring grains undergo heterogeneous shear rates parallel to the grain boundary. The present article focuses on understanding the plastic deformation fields near the grain boundaries and the influence of grain interaction on intra-grain deformations. Crystal Plasticity Finite Element Method (CPFEM) is employed on a periodic unit cell consisting of four grains discretised into a large number of elements. A refined study of the local variation of strain rates, both along and perpendicular to the grain boundaries permits an assessment of the assumptions made in the ALAMEL model. It is shown that the ALAMEL model imbibes the nature of plastic deformation at the grain boundaries very well. However, near triple junctions, the influence of a third grain induces severe oscillations of the stress tensor, reflecting a singularity. According to CPFEM, such singularity can lead to grain subdivision by the formation of new boundaries originating at the triple junction.     

Chemical potential tensor for a two-phase continuous medium model

Relations for jumps of thermodynamic variables with allowance for inertial terms are derived under conditions of thermal equilibrium and in the absence of dissipation on the interphase surface. The notion of the chemical potential tensor is generalized for this case within the framework of the elastic continuous medium model. A thermodynamically well-posed definition of the chemical potential tensor is proposed for a class of two-phase models of deformable solids. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 3, pp. 12–22, May–June, 2005.     

颗粒材料三维应力路径下的接触组构特性

颗粒材料的宏观应力变形特征与其微观接触力、组构等紧密相关.一般而言,强接触系统属于颗粒内部体系的传力结构,其对应的组构张量是影响宏观应力性质的重要因素.细观数值方法(如离散单元法)能够反映物理试验的基本规律,并且可以方便地提取宏微观数据来研究颗粒体系的应力变形机制.采用离散单元法(discrete element method,DEM)进行一系列等$p$等$b$应力路径下颗粒材料的真三轴试验,在此基础上研究了三维应力路径下颗粒材料的宏微观力学参数的演化过程、三维组构张量与应力张量多重联系以及强接触体系反映的宏观应力特征.研究表明:颗粒体系偏应力峰值状态和临界状态均存在与加载路径无关的宏微观特征;三维应力路径下组构张量与应力张量存在非共轴性,但其联合不变量演化过程表现出加载路径无关的特征;与弱接触系统的组构张量相比,强接触系统的组构张量更能反映宏观应力张量的特征;强弱接触体系的组构张量对颗粒体系宏观响应的贡献不同,其分界点存在一定取值范围,但采用平均接触力较为简单合理.      

Mean stress in a granular medium in dynamics

In this paper we propose to construct a mean stress tensor for a granular medium, valid in static and in dynamics, which takes into account the contact reactions and the body forces acting at the grain level. A simple analytical example shows that taking account of inertia forces is essential to insure the symmetry of the mean stress tensor in dynamics. Finally, numerical simulations illustrating this definition of the mean stress tensor are presented for a granular medium ensiled.     

Investigation of proper modeling of very dense granular flows in the recirculation system of CFBs

The aim of this paper is the development of new models and/or the improvement of existing numerical models, used for simulating granular flow in CFB (circulating fluidized bed) recirculation systems. Most recent models follow the TFM (two-fluid model) methodology, but they cannot effectively simulate the inter-particle friction forces in the recirculation system, because the respective stress tensor does not incorporate compressibility of flow due to change of effective particle density. As a consequence, the induced normal and shear stresses are not modeled appropriately during the flow of the granular phase in the CFB recirculation system. The failure of conventional models, such as that of von Mises/Coulomb, is mainly caused by false approximation of the yield criterion which is not applicable to the CFB recirculation system. The present work adopts an alternative yield function, used for the first time in TFM Eulerian modeling. The proposed model is based on the Pitman–Schaeffer–Gray–Stiles yield criterion. Both the temporal deformation of the solid granular phase and the repose angle that the granular phase forms are more accurately simulated by this model. The numerical results of the proposed model agree well with experimental data, implying that frictional forces are efficiently simulated by the new model.