Initial planar deformation of dilatant granular materials

A theory for the initial planar deformation of dilatant granular materials based on a kinematic proposal of R. Butterfield and R.M. Harkness (1972) is presented. The theory introduces an additional parameter called the angle of dilatancy into the traditional structure of plasticity theories for granular materials and soils. When the angle of dilatancy is zero, the present theory reduces to the theory introduced by A.J.M. Spencer in 1964. When the angle of dilatancy is equal to the angle of internal friction, the present theory reduces to the planar form of the theory introduced by D. C. Drucker and W. Prager in 1952. The properties of the theory presented here include coincidence of the stress and velocity characteristics, realistic energy dissipation predictions, and, in general, non-coincidence of the principal axes of stress and strain-rate. However, the angle of dilatancy is assumed to be a constant in this analysis and it does not decrease to zero with increased monotonic shearing deformation as experiment requires that it should, the theory therefore being limited to the initial deformation of dilatant granular materials.     

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

Continuum modeling of a porous solid with pressure-sensitive dilatant matrix

The pressure-sensitive plastic response of a material has been studied in terms of the intrinsic sensitivity of its yield stress to pressure and the presence and growth of cavities. This work focuses on the interplay between these two distinctly different mechanisms and the attendant material behavior. To this end, a constitutive model is proposed taking both mechanisms into account. Using Gurson's homogenization, an upper bound model is developed for a voided solid with a plastically dilatant matrix material. This model is built around a three-parameter axisymmetric velocity field for a unit sphere containing a spherical void. The void is also subjected to internal pressure; this can be relevant for polymeric adhesives permeated by moisture that vaporizes at elevated temperatures. The plastic response of the matrix material is described by Drucker–Prager's yield criterion and an associated flow rule. The resulting yield surface and porosity evolution law of the homogenized constitutive model are presented in parametric form. Using the solutions to special cases as building blocks, approximate models with explicit forms are proposed. The parametric form and an approximate explicit form are compared against full-field solutions obtained from finite element analysis. They are also studied for loading under generalized tension conditions. These computational simulations shed light on the interplay between the two mechanisms and its enhanced effect on yield strength and plastic flow. Among other things, the tensile yield strength of the porous solid is greatly reduced by the internal void pressure, particularly when a liquid/vapor phase is the source of the internal pressure.     

压力敏感材料中的空洞长大研究及其在页岩及高分子材料中的应用

主要研究压力敏感材料中含内压的空洞长大,如页岩或者高分子材料。采用数值方法研究含内压空洞的对称和非对称球形和柱形胞元的宏观力学行为。结果表明,压力敏感性及其空洞内压将极大影响空洞的形核与长大。在球形胞元情形中未出现柱形胞元的单轴拉伸现象。将胞元有限变形的数值计算结果与基于近期提出的考虑压力敏感材料中空洞长大的塑形力学模型的分析结果进行了对比。     

Dynamics and wave propagation in dilatant granular materials

The equations of motion for dilatant granular material are obtained from a Hamiltonian variational principle of local type in the conservative case. The propagation of nonlinear waves in a region with uniform state is studied by means of an asymptotic approach that has already appeared useful in an investigation on wave propagation in bubbly liquids and in fluid mixtures. When the grains are assumed to be incompressible, it is shown that the material behaves as a continuum with latent microstructure.

---

Sommario Si ricavano le equazioni di moto per i materiali granulari dilatanti da un principio variazionale Hamiltoniano di tipo locale nel caso conservativo. Si studia la propagazione delle onde non lineari in una regione di stato costante per mezzo di un approccio asintotico già rivelatosi utile nello studio della propagazione di onde nei liquidi con bolle e nelle miscele di fluidi. Quando si supponga che i granuli siano incomprimibili, si dimostra che il materiale si comporta come un continuo con microstruttura latente.

     

On the spatial localization of the laminar boundary layer in a non-Newtonian dilatant fluid

In dilatant fluids the shear perturbation propagation rate is finite, in contrast to Newtonian and pseudoplastic fluids in which it is infinite [1]. Therefore, in certain dilatant fluid flows, frontal surfaces separating regions with zero and nonzero shear perturbations may be formed. Since, in a sense, the boundary layer is a “time scan” of the nonstationary shear perturbation propagation process, in dilatant fluids the boundary layer should definitely be spatially localized. This was first mentioned in [2] where, however, it was mistakenly asserted that boundary layer spatial localization does not take place in all dilatant fluids and is absent in so-called “hardening” dilatant fluids. In [3], the solutions of the laminar boundary-layer equations for speudoplastic and “hardening” dilatant fluids were investigated qualitatively. The formation of frontal surfaces in dilatant fluid flows is usually mathematically related with the existence of singular solutions of the corresponding differential equations [4]. However, since the analysis performed in [3] was inaccurate, in that study singular solutions were not found and it was incorrectly concluded that in “hardening” dilatant fluids there is no spatial boundary layer localization. The investigation performed in [5] showed that in fact in “hardening” dilatant fluids boundary layers are spatially localized, since there exist singular solutions of the corresponding differential equations. Subsequently, this result was reproduced in [6], where an attempt was also made to carry out a qualitative investigation of the solutions of the laminar boundary-layer equations for other types of dilatant fluids. The author did not find singular solutions in this case and mistakenly concluded that in these fluids there is no spatial boundary layer localization. This misunderstanding was due to the fact that in [6] it was not understood that in dilatant fluid flows the formation of frontal surfaces can be mathematically described not only in relation to the existence of singular solutions.     

Research note on a simple model for pressure-sensitive strain-hardening materials

The true stress-plastic strain curves for pressure-sensitive high strength steels are successfully predicted using a simple constitutive model. Both the strength-differential and plastic compressibility effects are taken into account in the modeling. All material constants appearing in the constitutive relations have been obtained. The procedures outlined can be readily incorporated in standard elastoplastic finite element codes.     

A simplified second gradient model for dilatant materials: Theory and numerical implementation

This paper deals with the development of a new second gradient model, its numerical implementation and its validation. In order to remedy to the spurious mesh dependency of the post localized computation enhanced models incorporating some internal length are necessary. These models are very time consuming. In this paper we present a simplified theory within the framework of constrained micromorphic models involving only the micro volumetric strain. Provided the use of an additional penalty term in the numerical treatment, this model is quite efficient to regularize problems modelling behaviors exhibiting plastic volumetric strain such as the ones of geomaterials. More over this model is notably less time consuming than the more general ones.     

A gradient approach to localization of deformation. I. Hyperelastic materials

By utilizing methods recently developed in the theory of fluid interfaces, we provide a new framework for considering the localization of deformation and illustrate it for the case of hyperelastic materials. The approach overcomes one of the major shortcomings in constitutive equations for solids admitting localization of deformation at finite strains, i.e. their inability to provide physically acceptable solutions to boundary value problems in the post-localization range due to loss of ellipticity of the governing equations. Specifically, strain-induced localized deformation patterns are accounted for by adding a second deformation gradient-dependent term to the expression for the strain energy density. The modified strain energy function leads to equilibrium equations which remain always elliptic. Explicit solutions of these equations can be found for certain classes of deformations. They suggest not only the direction but also the width of the deformation bands providing for the first time a predictive unifying method for the study of pre- and post-localization behavior. The results derived here are a three-dimensional extension of certain one-dimensional findings reported earlier by the second author for the problem of simple shear.     

Three-dimensional constitutive relations for granular materials based on the dilatant double shearing mechanism and the concept of fabric

An extension of a three-dimensional model proposed by Anand and Gu (2000) for amorphous granular materials to include the effects of initial and induced anisotropy is presented in this paper. The proposed model can also be considered as a three-dimensional generalization of a model recently developed by Zhu et al. (2005) for the planar deformation of granular materials. The main ingredients of the model include the dilatant double shearing mechanism (Spencer, 1964, Mehrabadi and Cowin, 1978), the concept of fabric (Oda, 1972), and an extension of the Mohr–Coulomb yield criterion (Shield, 1955, Spencer, 1982) to three dimensions.The constitutive equations are implemented in the finite element program ABAQUS/Explicit (ABAQUS, 2001) by developing a user-material subroutine to conduct numerical triaxial compression tests for samples of granular materials with different initial anisotropy. The numerical results agree with the observed behavior and show that the extended constitutive model is capable of capturing the strength anisotropy of granular materials. Employing the anisotropic model developed here, we have also repeated the numerical simulation of the stress state in a static conical sand pile conducted earlier by Anand and Gu (2000). We find that fabric has little or no influence on the vertical stress distribution except at the base of the sand pile where the peak value of this stress is slightly higher than that predicted by the model of Anand and Gu (2000) which does not include the effects of fabric. We also find that the direction of the principal compressive stress changes from vertical at points away from the center of the pile to almost horizontal at points close to the center of the pile. This result provides a possible explanation for the observed dip in the vertical stress distribution in sand piles.     

Predicting the deformation characteristics of polymer materials

Theoretical methods of predicting the deformation characteristics of polymers are considered. The principles of reduced variables are examined. The relation between characteristics functions of polymers is analyzed. A prediction method is considered for polymers; the method is based upon the use of macroscopical mechanical models having a basis in terms of physics.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 2, pp. 147–152, March–April, 1973.     

Plastic deformation localization in commercial Zr-base alloys

The issues concerning the localization of plastic deformation in commercial Zr alloys used in the nuclear power industry are addressed. The possible types of deformation localization pictures corresponding with the respective stages of plastic flow are described. These are shown to be various kinds of self-excited wave processes of plastic flow. The dislocation structure of the material occurring within and in between the nuclei of localized deformation is investigated. The use of the self-excited wave patterns of plastic flow localization as an additional source of information on the mechanical properties of metals and alloys is substantiated.     

关于材料力学中小变形问题的讨论

着重讨论了材料力学中小变形假设下采用原尺寸原理和以直代曲方法所引起的相对误差,得出小变形假设的适用范围. 主要通过对称结构、非对称结构和超静定对称结构,分析应用原尺寸原理所得的结果与根据实际变形计算的结果之区别,得出小变形的定义是个相对的概念,与结构的初始角度和允许的应变有关;小变形假设的适用范围是由工程中的允许误差来确定的.     

Conditions for the onset of elastic and material instabilities in hyperelastic materials

Several constitutive inequalities have been proposed in the literature to quantify the notion that ‘stress increases with strain’ in an elastic material. Due to some inherent shortcomings in them, which we discuss, we propose a new tensorial criterion for isotropic materials. We also present necessary conditions in terms of elasticity tensors for the onset of elastic instabilities.     

Analysis of shear localization in porous materials

A set of constitutive equations are derived based on the authors' lower bound yield loci for porous materials. By using these equations, the conditions for shear localization in porous materials are then investigated and the results are compared with those of Gurson's equations and the finite element analysis. The advantages of the present constitutive equations are fully illustrated.     

On the deformation and fracture of porous materials

The constitutive behavior of porous materials (including the yield loci, the void growth rate, the macro stress-strain relation and the strain to localization instability) is examined based on the lower bound approach proposed by the present authors. These results are then compared with the experimental and the finite element results as well as those predicted by Gurson's equations. Emphasis is placed on approaching the real behavior from the upper and the lower bound analysis. Calculation is also made on the influence of void nucleation on the critical strain to instability and a modified strain-controlled nucleation criterion is proposed. Finally the instability and fracture of AISI4340 steel in plane strain tension is examined and comparison is made between theoretical and experimental results.     

Propagation of shear perturbations in dilatant liquids

Processes of shear-perturbation propagation in dilatant liquids are investigated. The evolution over time of the carrier of the Cauchy-problem solution for the quasilinear parabolic equation describing such processes is analyzed. Nonlinear effects of finite velocity and spatial localization of the shear perturbations are discussed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 60–65, November–December, 1978.