# Micromechanical Formulation of the Yield Surface in the Plasticity of Granular Materials

Authors: Homayoun Shaverdi; Mohd. Raihan Taha; Farzin Kalantary
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

An equation is proposed to unify the yield surface of granular materials by incorporating the fabric and its evolution. In microlevel analysis by employing a Fourier series that was developed to model fabric, it is directly included in the strength of granular materials. Inherent anisotropy is defined as a noncoaxiality between deposition angle and principal compressive stress. Stress-induced anisotropy is defined by the degree of anisotropyα and the major direction of the contact normals. The difference between samples which have the same density (or void ratio) but different bedding angles is attributed to this equation. The validity of the formulation is verified by comparison with experimental data.


---

## Body

## 1. Introduction

There are numerous experimental observations showing that the shape of the failure surface for soils is influenced by the microstructural arrangement (or fabric) (e.g., [1–3]). It has long been known that the failure condition is influenced by the microstructural arrangement of the constituent particles. Several expressions for failure criteria have been proposed to include the effect of fabric and its evolution. Baker and Desai [4] proposed the so-called joint isotropic invariants of stress and appropriate anisotropic tensorial entities. Pastor [5], by using this method, proposed a constitutive model to account for fabric anisotropy.Pietruszczak and Mroz [6] related inherent anisotropy to the microstructural arrangement within the representative volume of material. They used a second-order tensor whose eigenvectors specify the orientation of the axes of the material symmetry. The failure criteria proposed by Pietruszczak and Mroz [6] were formulated in terms of the stress state and a microstructure tensor. Lade [3], by using the method proposed by Pietruszczak and Mroz [6], related the loading directions to the principal directions of the cross-anisotropic microstructure arrangement of the particles.In order to connect the microscopic character of the granular materials with overall macroscopic anisotropy, various quantities have been proposed; for example, Oda [1], Oda et al. [2], and Oda [7] defined the fabric of anisotropy by using the distribution of the unit contact normals. Mehrabadi et al. [8] defined another microstructural arrangement and connected these parameters to the overall stress and other mechanical characteristics of granular materials. Gao et al. [9] and Gao and Zhao [10] proposed a generalized anisotropic failure criterion through developing an isotropic failure criterion by introducing two variables to account for fabric anisotropy. The first one is the fabric anisotropy that was proposed by Oda and Nakayama [11] and the second one is the joint invariants of the deviatoric stress tensor and the deviatoric fabric tensor to characterize the relative orientation between stress direction and fabric anisotropy. They related the frictional coefficient ηp to the anisotropic variable A. Fu and Dafalias [12] showed that there is a difference between friction angle in the isotropic and anisotropic cases. In the isotropic case, friction angle would be a direction-independent constant, while in the anisotropic case, it is a function of the bedding angle with respect to the shear plane (in the Mohr-Coulomb failure criterion). Fu and Dafalias [13] by using discrete element method (DEM) investigated the effect of fabric on the shear strength of granular materials. They proposed an anisotropic shear failure criterion on the basis of noncoaxiality between the bedding plane orientation and the shear plane. The inherent fabric anisotropy was taken into account by considering the orientation of the bedding plane with respect to the principal stress axes.The specification of the condition at failure for anisotropic granular soils constitutes an important problem and numerous criteria have been proposed in the past. In this paper, we endeavor to incorporate the effect of inherent and induced anisotropy in the yield surface. The inherent and induced anisotropies are expressed as explicit functions of the bedding angleβ and the magnitude of anisotropy α (in the distribution of contact normals). These two elements (inherent and induced anisotropy) are combined, and the Mohr-Coulomb yield surface which is modified to account for the kinematic yield surface [14–16] is developed by including the fabric and its evolution. The equation of the yield surface that is proposed for granular soils is compared with the experimental results from Oda et al. [17]. It shows that the equation is able to capture the shearing behavior of soils with different bedding angles.


## 2. Definition of Inherent Anisotropy

Inherent anisotropy is attributed to the deposition and orientation of the long axes of particles [1, 2, 7]. Oda et al. [17] and Yoshimine et al. [18] showed that the drained and undrained response of sand and approaching the critical state failure are actually affected by the direction of the principal stress relative to the orientation of the soil sample. Pietruszczak and Mroz [6] included the effect of fabric by the following equation:
(1)F=τ-ηg(θ)po,
where τ=J21/2 is the second invariant of the stress tensor, po=trσ/3 is first invariant of the stress tensor, g(θ) is Lode’s angle, and η is a constant for isotropic materials and defined by the following equation for anisotropic materials:
(2)η=ηo(1+Ωijlilj),
where ηo is the constant material parameter, Ωij describes the bias in material microstructure spatial distribution, and li and lj are the loading directions. Lade [3] by using these formulations proposed a failure criterion for anisotropic materials. Wan and Guo [19] accounted for the effect of inherent anisotropy in microlevel analysis by the ratio of projection of major-to-minor principal values of the fabric tensor along the direction of the principal stresses. Li and Dafalias [20, 21] incorporated this effect by the fabric tensor which was proposed by Oda and Nakayama [11]. These two methods used the same basic approach; they used the principal values of the fabric tensor in their formulations. However, micromechanical studies [2, 11] have shown that in the shearing process, the preferred orientation of the particles in a granular mass may undergo only small changes. Its value may well endure after the onset of the critical state; hence, the fabric anisotropy renders the locus of the critical state line. In this paper, cos2(βi-β∘) is used to model the effect of inherent anisotropy. βi indicates the variation of the long axes of particles with respect to the major principal stress; β∘ is the angle of deposition with respect to the major principal stress. Hence,
(3)(σ1σ2)f∝cos2(βi-β∘).


## 3. Definition of Stress-Induced Anisotropy

With increasing shear loads, the contact normals tend to concentrate in the direction of the major compressive stress. Contacts are generated in the compressive direction and disrupted in the tensile direction. These disruption and generation of the contact normals are the main causes of the induced anisotropy in the granular materials [2]. In order to include the fabric evolution (or induced anisotropy), a function in which changes of the contact normals are included must be defined. Wan and Guo [19] used the following equation:
(4)F˙ij=xη˙ij,
where F˙ij shows the evolution of fabric anisotropy, x is a constant, and η˙ij is the ratio of the shear stress to the confining pressure, or η=(q/p). Dafalias and Manzari [22] related the evolution of fabric to the volumetric strain in the dilatancy equation. The evolution of fabric comes to play only after dilation. Based on DEM simulation presented by Fu and Dafalias [12], Li and Dafalias [23] developed an earlier model (yield surface) to account for fabric and its evolution in a new manner by considering the evolution of fabric tensor towards its critical value.By using Fourier series, Rothenburg and Bathurst [24] showed that the contact normals distribution, E(n), can be presented as follows:
(5)E(n)=(12π)(1+αcos2(θ-θf)),
where α is the magnitude of anisotropy and θf is the major principal direction of the fabric tensor. The variations of the parameters α and θf represent the evolution of anisotropy in the granular mass. Experimental data shows that the shear strength of the granular material is a function of the magnitude of α and θf [1, 17, 25]. The following equation is used to consider the effect of the induced anisotropy:
(6)(σ1σ2)f∝(1+(12)αcos2(θσ-θf)).As previously mentioned, the shear strength in the granular medium is a function of inherent and induced anisotropy. The equation can predict the difference between samples due to the fabric which is a combination of the inherent and induced anisotropy as follows [26]:
(7)(σ1σ2)f∝[(1+(12)αcos2(θ-θf))cos2(βi-β∘)].Another parameter that must be added to the above relation is the rolling strength of the granular material. Oda et al. [25] and Bardet [27] showed the importance of the rolling strength of the particles, especially in a 2D case. This effect is incorporated in the following form [26]:
(8)(σ1σ2)f∝[(1+(12)αcos2(θ-θf))×cos2(βi-β∘)mexp(cos2(βi-β∘))(1+(12)αcos2(θ-θf))],
where m is a constant that depends on the interparticle friction angle ϕμ and the shape of the particles. When the samples with equal densities are subjected to the shear loads, the difference in the shear strength due to the fabric can be attributed to (8).


## 4. Verification of (8) with the Experimental Data

In order to show the ability of (8) to represent the effect of the fabric on the shear strength, the predictions are compared with the experimental tests presented by Konishi et al. [25]. They conducted an experimental study on biaxial deformation of two-dimensional assemblies of rod-shaped photoelastic particle with oval cross section. The samples were confined laterally by a constant force of 0.45 kgf and then compressed vertically by incremental displacement. Two types of particle shapes were used; one was r1/r2=1.1 and the other was r1/r2=1.4, in which r1 and r2 are the major and minor axes of cross section respectively. To consider the influence of friction, two sets of experiments were performed on these two particle shapes, one with nonlubricated particles of average friction angle of 52∘ and the other with particles which had been lubricated with an average friction angle of 26∘. The magnitude of the degree of anisotropy α and the major direction of the fabric θf are calculated by the following equations:
(9)A=∫02πE(θ)sin2θdθ,(10)B=∫02πE(θ)cos2θdθ,(11)θf=(12)arctan(AB).To show the ability of (8), the proportion of fabric with the shear strength variations is shown in Figure 1. The differences in the shear strength ratio at failure for different bedding angles are attributed to the differences in the developed anisotropic parameters. In other words, the combination of anisotropic parameters (for inherent and induced anisotropy) is proportional to the shear strength. The variation of right-hand side of (8) is proportional to the variation of shear strength ratio for different bedding angles. The right-hand side of (8) is shown by fabric anisotropy in Figure 1. The effect of bedding angle on stress ratio at failure for the different interparticle friction angle ϕμ is also shown in Figure 1.Samples withϕμ=52∘ and r1/r2=1.1 (a), ϕμ=26∘ and r1/r2=1.1 (b), ϕμ=52∘ and r1/r2=1.4 (c), and ϕμ=26∘ and r1/r2=1.4 (d).
(a)
(b)
(c)
(d)


## 5. Incorporation of the Fabric and Its Evolution in the Yield Surface

Muir Wood et al. [14] proposed the kinematic version of the Mohr-Coulomb yield surface as follows:
(12)f=q-ηyfpo,
where q is the deviatoric stress and ηyf is the size of the yield surface. Muir Wood et al. [14] and Muir Wood [16] assumed that the soil is a distortional hardening material; hence, the current yield surface ηyf is a function of the plastic distortional strain εqp, and, hence,
(13)ηyf=εqpc+εqpηp,
where ηp is a limit value of stress ratio which is equal to M at the critical state, ηp=M=q/p; c is a soil constant.Wood et al. [14] and Gajo and Muir Wood [15] developed the above equation to include the effect of state parameter ψ=e-ecr, in which e is the void ratio and ecr is the magnitude of the void ratio on the critical-state line, as follows:
(14)ηyf=εqpc+εqp(M-kψ),
where k is a constant.Li and Dafalias [20] modified the effect of state parameter ψ to account for a wide range of stress and void ratio as follows:
(15)ηp=Mexp(-nbψ),
where nb is a material constant. Equation (7) can be modified as follows:
(16)ηyf=εqpc+εqpMexp(-nbψ).In the previous section, the shear strength was shown to be a function of inherent and induced anisotropy (see (8)). Thus, the effect of inherent and induced fabric anisotropy for triaxial case can be expressed as follows:
(17)ηyf=εqpc+εqp(1+(12)αcos2(θf-θσ))×cos2(βi-β∘)Mexp(-nbψ).The magnitudes ofα and θf approach a constant value in large shear strain [26, 28, 29]. The parameter cos2(βi-β∘) is easily obtained by back calculation but as a rough estimation, its value is close to the magnitude of the bedding angle cosδ (for bedding angle δ between 15∘ and 45∘).Equation (10) can be shown in the following form for multiaxial direction (or in the general form):
(18)f=τ-ηyfg(θ)po.It is similar to the equation proposed by Pietruszczak and Mroz [6] and Lade [3] but in this formulation, another function is used for fabric and its evolution.


## 6. Fabric Evolution

The parametersα and θf show the status of the fabric and its evolution. These parameters have a great influence on the behavior of the dilatancy equation. Shaverdi et al. [29] proposed an equation which can predict the magnitude of α and θf in the presence of the noncoaxiality between stress and fabric. This equation is obtained from the microlevel analysis. To calculate the α parameter, the magnitude of the shear to normal stress ratio on the spatially mobilized plane (SMP) must be determined. In the triaxial case, for example, τ/p may be obtained from the following equation [30]:
(19)τp=σ1σ3-σ3σ1.The parametersα and θf may be obtained from the following equations in the presence of noncoaxiality [29]:
(20)α=(τ/p)cosϕμmob-sinϕμmobsin(2θf+ϕμmob)-((τ/p)cos(2θf+ϕμmob)),(21)θ˙f=θ˙σ+(12)·dη·(θσ-θf),
where the dot over θ shows the variation. The most important parameter in the above equation is the interparticle mobilized friction angle,  ϕμmob. This parameter is obtained from the following equation:
(22)tan-1(τp)=θσ-θfz+λ(ε˙vε˙q)+ϕμmob,
where z and λ are material constants. Kuhn [28] and Shaverdi et al. [29] showed that the variation of α with the shear strain is similar to the variation of shear to normal stress ratio with shear strain.


## 7. Verification of (16) with Experimental Data

Oda et al. [17] conducted some experimental tests on Toyoura sand with an initial void ratio 0.67-0.68 and the confining pressures 0.5 kg/cm2 and 2 kg/cm2 to study the effect of bedding angle with tilting angles δ= 0∘,30∘,60∘, and 90∘. For better modeling, the constant c for the plastic shear strain less than 1% must be modified as follows:
(23)c=0.001+0.001〈1-m′εqp〉,
where m′ is a constant which depends on confining pressure and 〈⋯〉 stands for the positive values only. In this simulation, m′=1 is taken into account for the confining pressure 0.5 kg/cm2, and m′=6.5 for the confining pressure 2.0 kg/cm2.Equation (17) can model the different behavior of the granular soils with the same confining pressure and initial void ratio (density) in which the only difference is due to the fabric and its evolution, as shown in Figures 2 and 3.Comparison between experimental data and simulation by using (16) for the confining pressure 0.5 kg/cm2.
(a)
(b)
(c)
(d)Comparison between experimental data and simulation by using (16) for the confining pressure 2.0 kg/cm2.
(a)
(b)
(c)
(d)


## 8. Conclusion

An equation was proposed to include the effect of inherent and induced anisotropy. This relation was obtained by combining the effect of inherent and induced anisotropy. Rolling resistance is also included in this equation. The differences between the samples due to inherent and induced anisotropy were well captured by applying (8). Verifying the experimental data shows that this equation can predict the ratio of the shear strength at failure of granular materials in the presence of inherent anisotropy as good as possible. The effect of inherent anisotropy was incorporated by a single term cos2(βi-β∘). Induced anisotropy was also included by a simple term (1+(1/2)αcos2(θf-θσ)) in which α and θf can be easily calculated and obtained. The extended Mohr-Coulomb was developed to incorporate the effect of fabric and its evolution. Verification with the experimental tests demonstrated the validity of this formulation.



---

*Source: 385278-2013-09-22.xml*