On mechanics of saturated granular materials

A continuum theory of saturated granular materials is formulated. The basic balance laws for the solid phase as well as for the fluid phase are presented. The constitutive equations are derived and the basic equations of motion of the solid and fluid continua are obtained. Several cases of interest, such as incompressible granules saturated with liquids are discussed. It is shown that the theory contains, as its special cases, the Mohr-Coulomb criterion for a granular material as well as Darcy's law of flow through porous media.     

Numerical solution to the shearing flow of granular materials between two plates

We study the shearing flow of granular materials between two horizontal flat plates where the top plate is moving with a constant speed. The constitutive relation used for the stress is based on the continuum model proposed by Rajagopal and Massoudi (DOE Report, DOE/PETC/TR-90/3, 1990). The material coefficients such as viscosity and normal stress coefficients are based on the model of Boyle and Massoudi (Int. J. Eng. Sci 28 (1990) 1261). The governing equations are non-dimensionalized and the resulting system of non-linear differential equations is solved numerically using finite difference technique.     

New constitutive equations for the rapid flow of granular materials—II

This paper compliments a previous paper, (J. Non-Newtonian Fluid Mech.,9(1981) 147), which discussed new constitutive equations for rapid collisionally maintained flow of granular materials as a non-Newtonian microfluid in which the gradient of microrotation played an important role as a kinematic variable. In this article we discuss another variation of such constitutive equations; one for which we do not consider the effects of gradients of the mean microrotation of grains. Derived are expressions for the dispersive normal and shear stresses for a plane shear rapid flow of a granular material in the presence of intergranular slip and friction. These expressions reduce to classical results if friction is set equal to zero. A graph of the variation of stresses versus the friction factor is also presented which reveals a kind of choking phenomenon at larger values of μ_k.     

Plane mechanics and kinematics of compressible ideal granular materials

Summary In a previous paper (1) proposals were made for the equations which govern the mechanical behaviour in plane strain of an ideal incompressible granular material. In this paper tentative suggestions are made for the extension of this theory to compressible granular materials. The material is envisaged as a mixture of an ideal gas and solid particles. The state of an element of the mixture is determined by its overall density and the mass and volume concentrations of the constituents. It is proposed that the material satisfies a flow condition which relates the two principal invariants of the stress (in two dimensions) and the density. The condition for the stress equations to be hyperbolic is obtained, and for the hyperbolic case a natural interpretation is obtained for the angle of internal friction.For the kinematic behaviour it is proposed that the deformation consists of superposed shearing deformations on the stress characteristic surfaces, as described in (1), together with a superposed dilatation. The equations describing this behaviour are expressed in characteristic form. They reduce to the equations given in (1) when the density is constant.     

散料在料仓内流动特性的实验研究

用振动筛测量多种散料的筛分尺寸，测量散料的容积密度、安息角，用应变直剪仪进行散料内摩擦角、壁面摩擦角的测量，分析了散料颗料的形状、尺寸对散料的安息角、内摩擦角、壁面摩擦角的影响。用示踪法进行了圆筒内散料重力场中颗粒流线和等时线的试验，以及出口直径与流量的关系、筒仓内散料高度对流量的影响的试验，流量与筒仓内散料高度关系的试验。分析了散料各物性参数对散料流动性的影响，为散料存贮和运输设备的结构设计提供实验依据。     

煤仓内煤散料流动状态与力学行为影响因素

针对煤仓内煤散料流动问题及其力学行为,采用三维颗粒流模拟程序PFC3D建立了某型号煤仓与某种煤散料的离散元模型,简述了其力学模型与求解步骤,模拟分析了煤仓内煤散料卸料流动状态。通过分析水平向侧压力、颗粒速度场和接触力场,重点讨论了煤仓下部锥体内壁面摩擦系数、锥仓倾角和卸料口径等对煤散料颗粒流动状态和力学行为的影响。结果显示,深仓卸料流动为整体流动与中心流动混合状态,煤仓内壁摩擦系数、锥体倾角和卸料孔开口半径均对煤散料流动和水平侧压力有较大影响。     

一种改进的颗粒移动床μ(I)拟流体模型及应用

颗粒移动床在工业领域应用广泛, 发展实用可靠的颗粒移动床模型具有理论和应用价值. 本文基于颗粒流μ(I)模型, 补充局部颗粒体积分数与颗粒局部压力和局部颗粒流密度的关系式, 将移动床内密集颗粒处理成可压缩拟流体, 建立了颗粒流单相可压缩流μ(I)模型, 并建立了颗粒流?壁面摩擦条件, 在计算中对颗粒流拟黏度和拟压力项进行正则化处理. 采用上述模型与方法对3种典型散料在移动床缩口料仓内的流动进行模拟, 与实验对比, 得到了玻璃珠、刚玉球和粗沙的μ(I)模型参数, 分析了3种不同散料在料仓内的颗粒速度、体积分数等分布特性, 模拟结果较好地揭示了料仓内不同物料的整体流和漏斗流特性; 进而以玻璃珠为例, 对移动床颗粒单管绕流流动进行了模拟, 所得结果合理揭示了管流附近的流动特性. 计算结果表明, 对于本文的计算工况, 颗粒体积分数变化最大范围为0.510 ~ 0.461, 绝大部分区域流动惯性数小于0.1, 改进的单相μ(I)模型能合理预测出密集颗粒流移动床内的流动特性, 方法可行且较多相流算法能明显减小计算量.      

Incorporating the effects of fabric in the dilatant double shearing model for planar deformation of granular materials

The objective of this paper is to incorporate the effects of fabric and its evolution into the Dilatant Double Shearing Model [Mehrabadi, M.M., Cowin, S.C., 1978. Initial planar deformation of dilatant granular materials. J. Mech. Phys. Solids 26, 269–284] for granular materials in order to capture the anisotropic behavior and the complex response of granular materials in cyclic shear loading. An important consequence of considering the fabric is that one can have unequal shearing rates along the two slip directions. This property leads to the non-coaxiality of the principal axes of stress and strain rate, which is more appropriate for a material that exhibits initial and induced anisotropy. In addition, we employ a fabric-dependent elasticity tensor with orthotropic symmetry. The model developed in this paper also predicts one of the experimentally observed characteristics of granular materials: the gradual concentration of the contact normals towards the maximum principal stress direction.We implement the constitutive equations into ABAQUS/Explicit by writing a user material subroutine in order to predict the strength anisotropy of granular materials in a plane strain biaxial compression test and investigate the mechanical behavior of granular materials under the cyclic shear loading conditions. The predictions from this model show good quantitative agreement with the experiments of [Park, C.S., 1990. Anisotropy in deformation and strength properties of sands in plane strain compression, Masters Thesis, University of Tokyo; Park, C.S., Tatsuoka, F., 1994. Anisotropic strength and deformation of sands in plane strain compression. In: XIII ICSMFE, New Delhi, India; Okada, N., 1992. Energy dissipation in inelastic flow of cohesionless granular media. Ph.D. Thesis, University of California, San Diego].     

Theory of Short-Term Microdamageability of Granular Composite Materials

The theory of microdamageability of granular composites is outlined through the simulation of microdamages in the components by pores filled with compression-resisting particles of a destroyed material. The damage criterion for a microvolume of a component is taken in the Schleicher–Nadai form, which allows for the difference between the ultimate tensile and compressive loads. The ultimate strength is a random function of Weibull-distributed coordinates. The stress–strain state and the efficient properties of the material are determined from the stochastic equations of elastic theory for a granular composite with porous components. The equations of deformation and microdamage are closed by the equations of porosity balance in the components. Nonlinear diagrams of the concurrent processes of deformation in the granular material and microdamage in the matrix are plotted. The effect of the physical and geometrical parameters on them is studied     

Solution existence conditions for elastoplastic constitutive models of granular materials

In this paper, the conditions of solution existence for stress rates under given strain rates are investigated. The focus of the solution existence investigation is on the non-associated flow rule and elastic stress–strain relationship. Granular materials characterized with strong non-associated plastic flows are used as a particular example for analysis. Various flow rules for granular materials are analyzed, including Rowe’s, Roscoe’s flow rules and their modified versions. In the elastic stress–strain relationships of materials, the effects of Poisson’s ratio on solution existence are investigated. Both isotropic and anisotropic elasticity are considered. Given a granular material and its states, it is found that there exists a critical Poisson’s ratio for a particular non-associated flow rule. When the Poisson’s ratio of a material is above this critical Poisson’s ratio, its constitutive model is susceptible to solution non-existence. It is suggested that special attentions should be paid to the selection of material Poisson’s ratio and non-associated flow rule to ensure the existence of elastoplastic solutions.     

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.     

Volume-weighted mixture theory for granular materials

In the present work we treat granular materials as mixtures composed of a solid and a surrounding void continuum, proposing then a continuum thermodynamic theory for it. In contrast to the common mass-weighted balance equations of mass, momentum, energy and entropy for mixtures, the volume-weighted balance equations and the associated jump conditions of the corresponding physical quantities are derived in terms of volume-weighted field quantities here. The evolution equations of volume fractions, volume-weighted velocity, energy, and entropy are presented and explained in detail. By virtue of the second law of thermodynamics, three dissipative mechanisms are considered which are specialized for a simple set of linear constitutive equations. The derived theory is applied to the analysis of reversible and irreversible compaction of cohesionless granular particles when a vertical oscillation is exerted on the system. In this analysis, a hypothesis for the existence of a characteristic depth within the granular material in its closely compacted state is proposed to model the reversible compaction.     

Influence of porosity on the non-linear deformation of granular materials

An algorithm is suggested to determine non-linear deformative properties of granular composite materials with non-linearly deformed porous components. Stochastic differential equations of physically non-linear elasticity theory with the following use of conditional averaging are assumed as a basis. Solution of non-linear equations with respect to conditional moment is constructed by the method of successive approximation. The volume content of the components and pores in the components have been studied for their effect on the character of the strain curves of porous granular composite.     

A kinetic model for rapid flows of granular materials

A kinetic model for rapid flows of granular material is considered. An evolution equation for the first order distribution function is developed which includes the effects of intergranular friction force in an average sense. The collision operator is approximated by the BGK relaxation model. The fundamental equations of motion including an equation for dynamics of the fluctuation energy are derived and discussed. The gravity and the plane Couette flows of granular materials are treated as examples of the applications of the present theory.     

A three-invariant cap model with isotropic–kinematic hardening rule and associated plasticity for granular materials

In this paper, a three-invariant cap model is developed for the isotropic–kinematic hardening and associated plasticity of granular materials. The model is based on the concepts of elasticity and plasticity theories together with an associated flow rule and a work hardening law for plastic deformations of granulars. The hardening rule is defined by its decomposition into the isotropic and kinematic material functions. The constitutive elasto-plastic matrix and its components are derived by using the definition of yield surface, material functions and non-linear elastic behavior, as function of hardening parameters. The model assessment and procedure for determination of material parameters are described. Finally, the applicability of proposed plasticity model is demonstrated in numerical simulation of several triaxial and confining pressure tests on different granular materials, including: wheat, rape, synthetic granulate and sand.     

激波与堆积粉尘相互作用的实验和理论研究

通过对激波与堆积粉尘相互作用的实验和理论分析，可得到堆积粉尘本构方程。基于该方程，本文对激波与堆积粉尘相互作用现象进行了数值模拟。计算所反映的流场结构与实验图像一致。     

Ram extrusion force for a frictional plastic material: model prediction and application to cement paste

We developed a model to predict the ram extrusion force of frictional plastic materials such as cement-based pastes. The extrusion of cement-based materials has already been studied, but the interaction between shaping force and paste behaviour still have to be understood. Our model is based on the plastic frictional behaviour of cement-based materials and integrates the physical mechanisms that govern material extrusion flow and extrusion force increase. When the process starts, a pressure gradient is created in the extruder due to wall friction of the paste that is submitted to plug flow. It induces a consolidation of the material. As a result, a large increase of extrusion force appears. A Coulomb law is used to model cement-based materials, which is considered as consolidating granular media. Such modelling is compared with experimental results. Tests were carried out on extrudible cement pastes. Modelling and experimental results are in good agreement.Paper presented at the Annual European Rheology Conference, Grenoble, April 2005     

Inertial steady flow of a bed of granular material down a surface with microrelief with allowance for finite time of intergranular contacting

An inertial flow of a granular material can be described by the laws of conservation of mass, momentum, and energy of random motion of solid particles by invoking some closing relations. In this work, these closing relations are inferred from the dimensional theory. The system of equations obtained is used to determine characteristics of a steady flow of a bed of a granular material down an inclined surface with a microrelief for various Richardson numbers and finite contact times of the particles during their collisions. Novosibirsk Military Institute, Novosibirsk 630103. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 6, pp. 128–132, November–December, 1999.     

Numerical simulations of granular free-surface flows using smoothed particle hydrodynamics

This paper presents a two-dimensional SPH model designed to simulate free-surface flows of dense granular materials. Smoothed particle hydrodynamics (SPH) is a mesh-free numerical method based on a Lagrangian discretization of the continuum mass and momentum conservation equations. The rheology of dense granular materials is modelled using a new local constitutive law recently proposed by Jop et al. (Nature, 2006). Of the viscoplastic class, this law is characterized by an apparent viscosity depending both on the local strain rate and local pressure. Validation test cases of the model in steady and unsteady configurations are presented. For steady cases (vertical chute flow and uniform free-surface layers on inclines), excellent agreement with analytical predictions is obtained. In the unsteady case, the simulations satisfactorily capture the dynamics of gravity-driven surges observed in experiments, including behaviours that are very specific to granular materials. Among the various parameters involved in the computations, the influence of SPH particle configuration within the flow and of the threshold viscosity used in the regularization of the constitutive yield criterion are particularly discussed.     

Analytical solutions for tapering quadratic and cubic rat-holes in highly frictional granular solids

New exact analytical solutions are presented for both stress and velocity fields for a Coulomb–Mohr granular solid assuming non-dilatant double-shearing theory. The solutions determined apply to highly frictional materials for which the angle of internal friction φ is assumed equal to 90°. This major assumption is made primarily to facilitate exact analytical solutions, and it is discussed at length in the Introduction, both in the context of real materials which exhibit large angles of internal friction, and in the context of using the solutions derived here as the leading term in a regular perturbation solution involving powers of 1−sinφ. The analytical velocity fields so obtained are illustrated graphically by showing the direction of the principal stress as compared to the streamlines. The stress solutions are also exploited to determine the static stress distribution for a granular material contained within vertical boundaries and a horizontal base, which is assumed to have an infinitesimal central outlet through which material flows until a rat-hole of parabolic or cubic profile is obtained, and no further flow takes place. A rat-hole is a stable structure that may form in storage hoppers and stock-piles, preventing any further flow of material. Here we consider the important problems of two-dimensional parabolic rat-holes of profile y=ax², and three-dimensional cubic rat-holes of profile z=ar³, which are both physically realistic in practice. Analytical solutions are presented for both two and three-dimensional rat-holes for the case of a highly frictional granular solid, which is stored at rest between vertical walls and a horizontal rigid plane, and which has an infinitesimal central outlet. These solutions are bona fide exact solutions of the governing equations for a Coulomb–Mohr granular solid, and satisfy exactly the free surface condition along the rat-hole surface, but approximate frictional conditions along the containing boundaries. The analytical solutions presented here constitute the only known solutions for any realistic rat-hole geometry, other than the classical solution which applies to a perfectly vertical cylindrical cavity.