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].     

Simulated configurational temperature of particles and a model of constitutive relations of rapid-intermediate-dense granular flow based on generalized granular temperature

In gas–solid flat-base spout bed with a jet, the flow of particles must go through an intermediate regime where both kinetic/collisional and frictional contributions play a role. In this paper, the statistical framework is proposed to define the generalized granular temperature which sums up the configurational temperature and translational granular temperature. The configurational temperature, translational and rotational granular temperatures of particles are simulated by means of CFD-DEM (discrete element method) in a 3D flat-base spout bed with a jet. The configurational temperatures of particles are calculated from instantaneous overlaps of particles. The translational and rotational granular temperatures of particles are calculated from instantaneous translational and angular velocities of particles. Roughly, the simulated translational and rotational granular temperatures increase, reach maximum, and then decrease with the increase of solids volume fractions. However, the configurational temperature increases with the increase of solids volume fractions. At high solid volume fraction, the predicted configurational temperatures are larger than the translational and rotational granular temperatures, indicating that the rate of energy dissipation do contributes by contact deformation of elastic particles. The generalized granular temperature is proposed to show the relation between the variance of the fluctuation velocity of deformation and the variance of the translational fluctuation velocity of particles. The constitutive relations of particle pressure, viscosity, granular conductivity of fluctuating energy and energy dissipation in rapid-intermediate-dense granular flows are correlated to the generalized granular temperature. The variations of particle pressure, shear viscosity, energy dissipation and granular conductivity are analyzed on the basis of generalized granular temperature in a flat-base spout bed with a jet. The axial velocities of particles predicted by a gas–solid two-fluid model of rapid-intermediate-dense granular flows agree with experimental results in a spout bed.     

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.     

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.     

Three-dimensional constitutive relations of ideal plasticity and the flow on the coulomb-tresca prism edge

In the present paper, we consider basic relations of the mathematical theory of plasticity for the spatial state corresponding to the edge of the Coulomb-Tresca prism, which follow from the generalized associated flow law restricting the plastic flow freedom for the above states to the minimal possible extent. We found that the spatial relations of the theory of plasticity, formulated by A. Yu. Ishlinsky in 1946, can be derived from the above version of the theory of flow. We show that the A. Yu. Ishlinsky constitutive relations for states on the Coulomb-Tresca prism edge express the commutativity of the stress tensor and the tensor of plastic strain increments. We obtained one explicit form of the constitutive relation relating the stress tensor to the plastic strain increments for the stressed states corresponding to the Coulomb-Tresca prism edge.     

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.     

Structure of the constitutive relations for hereditarily elastic materials reinforced by hard fibers

The problem of simplifying the nonlinear hereditary elasticity relations is considered for strongly anisotropic materials such as fiber-reinforced composites. This is done using their property that the material stiffness is high along the reinforcement and is low in the cross-sectional direction. The material is considered transversally isotropic. The simplification is performed by analyzing asymptotic representations of creep relations. Relations of various degrees of accuracy for various types of composites and stress states are obtained. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 3, pp. 120–127, May–June, 2005.     

Three-dimensional constitutive model for shape memory alloys based on microplane model

A new model for the behavior of polycrystalline shape memory alloys (SMA), based on a statically constrained microplane theory, is proposed. The new model can predict three-dimensional response by superposing the effects of inelastic deformations computed on several planes of different orientation, thus reproducing closely the actual physical behavior of the material. Due to the structure of the microplane algorithm, only a one-dimensional constitutive law is necessary on each plane. In this paper, a simple constitutive law and a robust kinetic expression are used as the local constitutive law on the microplane level. The results for SMA response on the macroscale are promising: simple one-dimensional response is easily reproduced, as are more complex features such as stress-strain subloops and tension-compression asymmetry. A key feature of the new model is its ability to accurately represent the deviation from normality exhibited by SMAs under nonproportional loading paths.     

基于扰动状态概念与E-B模型的粗粒料力学行为模拟

考虑到基于复杂弹塑性模型与扰动状态概念的本构模型在数值积分时的困难,以Duncan-Chang(E-B)模型描述材料在相对完整状态(RI)下的响应,用临界状态模型描述材料在完全调整状态(FA)下的力学行为。依据常规三轴试验曲线判断材料状态的转变并建议了相应的扰动因子计算方法。利用Abaqus软件实现了相应的程序代码,数值算例模拟了堆石料的三轴试验,结果表明邓肯-张模型与扰动状态概念的结合在模拟应变软化、剪缩与剪胀方面具有良好能力。在此基础上进一步模拟了某混凝土面板堆石坝的施工、蓄水、退水的过程,数值结果表明变形分布符合一般规律,且在定量上与原型观测吻合较好。     

On a thermodynamical consistent constitutive law based on the concept of internal variables

A multiaxial constitutive law describing the behaviour of polycrystalline ice in a phenomenological way will be introduced. Internal variables represent the hardening and softening of the material and also a modified inelastic work. The equations are implemented in a thermodynamical frame.To derive the unknown functions and parameters a strategy is developed, facilitating the fitting in a simple and fast manner.Finally, the comparison between computations and experiments shows a good agreement.     

A model-reduction approach in micromechanics of materials preserving the variational structure of constitutive relations

In 2003 the authors proposed a model-reduction technique, called the Nonuniform Transformation Field Analysis (NTFA), based on a decomposition of the local fields of internal variables on a reduced basis of modes, to analyze the effective response of composite materials. The present study extends and improves on this approach in different directions. It is first shown that when the constitutive relations of the constituents derive from two potentials, this structure is passed to the NTFA model. Another structure-preserving model, the hybrid NTFA model of Fritzen and Leuschner, is analyzed and found to differ (slightly) from the primal NTFA model (it does not exhibit the same variational upper bound character). To avoid the “on-line” computation of local fields required by the hybrid model, new reduced evolution equations for the reduced variables are proposed, based on an expansion to second order (TSO) of the potential of the hybrid model. The coarse dynamics can then be entirely expressed in terms of quantities which can be pre-computed once for all. Roughly speaking, these pre-computed quantities depend only on the average and fluctuations per phase of the modes and of the associated stress fields. The accuracy of the new NTFA-TSO model is assessed by comparison with full-field simulations. The acceleration provided by the new coarse dynamics over the full-field computations (and over the hybrid model) is then spectacular, larger by three orders of magnitude than the acceleration due to the sole reduction of unknowns.     

Hyperelastic constitutive model for rubber-like materials based on the first Seth strain measures invariant

The mechanical behaviour of isotropic and incompressible vulcanized natural rubbers (NR's) and that of quasi-incompressible carbon black filled vulcanized natural rubbers (NR 70) are considered both theoretically and experimentally. We start by generalising the neo-Hookean model to derive an original form of the strain energy density function (W). W satisfies the hypothesis of the Valanis–Landel function, which allows reducing the number of needed experimental tests to identify the parameters of the model. In the present study, in the order to identity the analytical form of W, we undertake only simple tension tests. The two-dimensional field of in-plane homogeneous displacements is determined here using a home-developed image analysis cross-correlation technique. Our model is also identified using results taken from the literature in the case of (NR's) for different types of solicitations, including simple tension, equibiaxial tension and pure shear deformation. Comparison of numerical results with the experimental data indicates that the present model can characterise the hyperelastic behaviour of NR's and that of NR 70 for all the tested modes of deformation. Moreover, it seems to be valid over a wide range of deformation intervals.     

基于近场动力学理论的准脆性材料的本构力函数的构建

近场动力学理论(PD)是基于非局部思想的连续介质力学新理论,用于研究材料破坏问题。根据准脆性材料破坏的线性和非线性的力学行为,在初始微观弹脆性材料(PMB)的本构力函数中引入了键的损伤模型,将键的断裂过程分成了线性的弹性变形阶段和非线性的损伤变形阶段,以此构建了准脆性材料的本构力函数的基本形式。以典型的准脆性材料为例构建了其本构力函数,通过在压缩载荷下对含预制不同角度单裂纹缺陷的类岩材料的裂纹扩展进行PD数值模拟仿真,裂纹起裂位置和扩展方向与试样试验结果在一定程度上保持了一致,证明了该基于近场动力学理论的典型准脆性材料的本构力函数可用于该类材料的破坏分析。     

A constitutive model in linear thermoviscoelasticity of polymers based on the concept of cooperative relaxation

New constitutive relations are derived for amorphous glassy polymers based on the concept of cooperative relaxation. A polymer is treated as a system of rearranging regions (flow units) embedded into a homogeneous elastic matrix. The viscoelastic (time-dependent) response of a medium is explained by rearrangements of segments of long chains in relaxing regions which occur at random instants. The kinetics of rearrangement is described in the framework of the Eyring concept of thermally activated processes, whereas the energy of any flow unit is assumed to randomly change at the instant of its reformation. Based on experimental data, phenomenological formulas are proposed for material functions. Adjustable parameters are found by fitting observations for mixtures of nylon with lithium halides in isothermal tensile relaxation tests. The thermoviscoelastic response in other tests is studied numerically. It is demonstrated that the material behavior predicted by the constitutive model in non-isothermal tests substantially differs from that predicted by conventional models whose adjustable parameters are determined by using the same experimental data. Received September 30, 1998     

An integral constitutive law for viscoelastic fluids based on the concept of evolving natural configurations: stability analysis

A general conceptual framework has recently been developed within the context of large deformations for the study of materials (solids and fluids) that either undergo microstructural changes or for which stress states are related to the presence of several relaxation mechanisms. The corner stone of this new approach is that all materials exhibit an infinity of stress-free natural configurations that are evolving in a thermodynamically admissible process. In this paper, that is exploratory in nature, we study the behavior of a non-linear integral viscoelastic constitutive equation (CE) for a fluid assumed to posses an infinity of evolving natural configurations. The CE stability pattern with respect to small perturbations about the rest state is also addressed.     

Scale-dependent constitutive relations and the role of scale on nominal properties

Size effects in strength and fracture energy of heterogeneous materials is considered within a context of scale-dependent constitutive relations. Using tools of wavelet analysis, and considering the failure state of a one-dimensional solid, constitutive relations which include scale as a parameter are derived from a ‘background’ gradient formulation. In the resulting theory, scale is not a fixed quantity independent of deformation, but rather directly dependent on the global deformation field. It is shown that strength or peak nominal stress (maximum point at the engineering stress–strain diagram) decreases with specimen size while toughness or total work to fracture per nominal area (area under the curve in the engineering stress–strain diagram integrated along the length of the considered one-dimensional specimen) increases. This behavior is in agreement with relevant experimental findings on heterogeneous materials where the overall mechanical response is determined by variations in local material properties. The scale-dependent constitutive relations are calibrated from experimental data on concrete specimens.