Journal of Applied Mathematics

Volume 2013 (2013), Article ID 385278, 7 pages

http://dx.doi.org/10.1155/2013/385278

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

^{1}Universiti Kebangsaan Malaysia (UKM), 43600 Bandar Baru Bangi, Selangor, Malaysia^{2}Department of Geotechnical Engineering, Faculty of Civil Engineering, K. N. Toosi University of Technology, 1346 Vali Asr Street, Vanak, Tehran 19697, Iran

Received 7 June 2013; Revised 31 July 2013; Accepted 15 August 2013

Academic Editor: Pengcheng Fu

Copyright © 2013 Homayoun Shaverdi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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

#### 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 to the anisotropic variable . 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: where is the second invariant of the stress tensor, is first invariant of the stress tensor, is Lode’s angle, and is a constant for isotropic materials and defined by the following equation for anisotropic materials: where is the constant material parameter, describes the bias in material microstructure spatial distribution, and and 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, is used to model the effect of inherent anisotropy. 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. 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: where shows the evolution of fabric anisotropy, is a constant, and is the ratio of the shear stress to the confining pressure, or . 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, , can be presented as follows: where is the magnitude of anisotropy and is the major principal direction of the fabric tensor. The variations of the parameters and 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 [1, 17, 25]. The following equation is used to consider the effect of the induced anisotropy:

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

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]: where 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 and the other was , in which and 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 and the other with particles which had been lubricated with an average friction angle of . The magnitude of the degree of anisotropy and the major direction of the fabric are calculated by the following equations:

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.

#### 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: where is the deviatoric stress and 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 is a function of the plastic distortional strain , and, hence, where is a limit value of stress ratio which is equal to at the critical state, ; 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 , in which is the void ratio and is the magnitude of the void ratio on the critical-state line, as follows: where 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: where is a material constant. Equation (7) can be modified as follows:

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:

The magnitudes of and approach a constant value in large shear strain [26, 28, 29]. The parameter is easily obtained by back calculation but as a rough estimation, its value is close to the magnitude of the bedding angle (for bedding angle between 15^{∘} and 45^{∘}).

Equation (10) can be shown in the following form for multiaxial direction (or in the general form):

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 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 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, may be obtained from the following equation [30]:

The parameters and may be obtained from the following equations in the presence of noncoaxiality [29]: where the dot over shows the variation. The most important parameter in the above equation is the interparticle mobilized friction angle, . This parameter is obtained from the following equation: where 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/cm^{2} and 2 kg/cm^{2} to study the effect of bedding angle with tilting angles = , and . For better modeling, the constant for the plastic shear strain less than 1% must be modified as follows:
where is a constant which depends on confining pressure and stands for the positive values only. In this simulation, is taken into account for the confining pressure 0.5 kg/cm^{2}, and for the confining pressure 2.0 kg/cm^{2}.

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.

#### 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 . Induced anisotropy was also included by a simple term in which and 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.

#### References

- M. Oda, “Initial fabrics and their relations to mechanical properties of granular materials,”
*Soils and Foundations*, vol. 12, no. 1, pp. 17–36, 1972. View at Google Scholar - M. Oda, S. Nemat-Nasser, and J. Konishi, “Stress-induced anisotropy in granular mass,”
*Soils and Foundations*, vol. 25, no. 3, pp. 85–97, 1985. View at Google Scholar · View at Scopus - P. V. Lade, “Failure criterion for cross-anisotropic soils,”
*Journal of Geotechnical and Geoenvironmental Engineering*, vol. 134, no. 1, pp. 117–124, 2008. View at Publisher · View at Google Scholar · View at Scopus - R. Baker and C. S. Desai, “Induced anisotropy during plastic straining,”
*International Journal for Numerical and Analytical Methods in Geomechanics*, vol. 8, no. 2, pp. 167–185, 1984. View at Google Scholar · View at Scopus - M. Pastor, “Modelling of anisotropic sand behaviour,”
*Computers and Geotechnics*, vol. 11, no. 3, pp. 173–208, 1991. View at Google Scholar · View at Scopus - S. Pietruszczak and Z. Mroz, “Formulation of anisotropic failure criteria incorporating a microstructure tensor,”
*Computers and Geotechnics*, vol. 26, no. 2, pp. 105–112, 2000. View at Publisher · View at Google Scholar · View at Scopus - M. Oda, “Fabric tensor and its geometrical meaning,” in
*Mechanics of Granular Materials: An Introduction*, M. Oda and K. Iwashita, Eds., pp. 19–33, Balkema, Rotterdam, The Netherlands, 1999. View at Google Scholar - M. Mehrabadi, S. Nemat-Nasser, and M. Oda, “On statistical description of stress and fabric in granular materials,”
*International Journal For Numerical and Analytical Methods in Geomechanics*, vol. 6, no. 1, pp. 95–108, 1982. View at Google Scholar - Z. Gao, J. Zhao, and Y. Yao, “A generalized anisotropic failure criterion for geomaterials,”
*International Journal of Solids and Structures*, vol. 47, no. 22-23, pp. 3166–3185, 2010. View at Publisher · View at Google Scholar · View at Scopus - Z. W. Gao and J. D. Zhao, “Efficient approach to characterize strength anisotropy in soils,”
*Journal of Engineering Mechanics*, vol. 138, no. 12, pp. 1447–1456, 2012. View at Google Scholar - M. Oda and H. Nakayama, “Yield function for soil with anisotropic fabric,”
*Journal of Engineering Mechanics*, vol. 115, no. 1, pp. 89–104, 1989. View at Google Scholar · View at Scopus - P. Fu and Y. F. Dafalias, “Fabric evolution within shear bands of granular materials and its relation to critical state theory,”
*International Journal for Numerical and Analytical Methods in Geomechanics*, vol. 35, no. 18, pp. 1918–1948, 2011. View at Publisher · View at Google Scholar · View at Scopus - P. Fu and Y. F. Dafalias, “Study of anisotropic shear strength of granular materials using DEM simulation,”
*International Journal for Numerical and Analytical Methods in Geomechanics*, vol. 35, no. 10, pp. 1098–1126, 2011. View at Publisher · View at Google Scholar · View at Scopus - D. M. Wood, K. Belkheir, and D. F. Liu, “Strain softening and state parameter for sand modelling,”
*Geotechnique*, vol. 44, no. 2, pp. 335–339, 1994. View at Google Scholar · View at Scopus - A. Gajo and D. Muir Wood, “Severn-Trent sand: a kinematic-hardening constitutive model: the $q-p$ formulation,”
*Geotechnique*, vol. 49, no. 5, pp. 595–614, 1999. View at Google Scholar · View at Scopus - D. Muir Wood,
*Geotechnical Modelling*, Taylor and Francis, 2004. - M. Oda, I. Koishikawa, and T. Higuchi, “Experimental study of anisotropic shear strength of sand by plane strain test,”
*Soils and Foundations*, vol. 18, no. 1, pp. 25–38, 1978. View at Google Scholar - M. Yoshimine, K. Ishihara, and W. Vargas, “Effects of principal stress direction and intermediate principal stress on undrained shear behavior of sand,”
*Soils and Foundations*, vol. 38, no. 3, pp. 179–188, 1998. View at Google Scholar · View at Scopus - R. G. Wan and P. J. Guo, “Drained cyclic behavior of sand with fabric dependence,”
*Journal of Engineering Mechanics*, vol. 127, no. 11, pp. 1106–1116, 2001. View at Publisher · View at Google Scholar · View at Scopus - X. S. Li and Y. F. Dafalias, “Constitutive modeling of inherently anisotropic sand behavior,”
*Journal of Geotechnical and Geoenvironmental Engineering*, vol. 128, no. 10, pp. 868–880, 2002. View at Publisher · View at Google Scholar · View at Scopus - X. S. Li and Y. F. Dafalias, “A constitutive, framework for anisotropic sand including non-proportional loading,”
*Geotechnique*, vol. 54, no. 1, pp. 41–55, 2004. View at Publisher · View at Google Scholar · View at Scopus - Y. F. Dafalias and M. T. Manzari, “Simple plasticity sand model accounting for fabric change effects,”
*Journal of Engineering Mechanics*, vol. 130, no. 6, pp. 622–634, 2004. View at Publisher · View at Google Scholar · View at Scopus - X. S. Li and Y. F. Dafalias, “Anisotropic critical state theory: role of fabric,”
*Journal of Engineering Mechanics*, vol. 138, no. 3, pp. 263–275, 2012. View at Publisher · View at Google Scholar · View at Scopus - L. Rothenburg and R. J. Bathurst, “Analytical study of induced anisotropy in idealized granular materials,”
*Geotechnique*, vol. 39, no. 4, pp. 601–614, 1989. View at Google Scholar · View at Scopus - J. Konishi, M. Oda, and S. Nemat-Nasser, “Inherent anisotropy and shear strength of assembly of oval cross-section rods,” in
*Proceedings of the IUTAM Symposium on deformation and failure of granular materials*, Delft, The Netherland, 1982. - H. Shaverdi, M. R. Taha, and F. Kalantary, “Effect of fabric on the shear strength of granular materials in biaxial compression,” in
*Proceedings of the 7th International Conference on Micromechanics of Granular Media*, pp. 615–618, Sydney, Australia, July 2013. View at Publisher · View at Google Scholar - J. P. Bardet, “Observations on the effects of particle rotations on the failure of idealized granular materials,”
*Mechanics of Materials*, vol. 18, no. 2, pp. 159–182, 1994. View at Google Scholar · View at Scopus - M. R. Kuhn, “Micro-mechanics of fabric and failure in granular materials,”
*Mechanics of Materials*, vol. 42, no. 9, pp. 827–840, 2010. View at Publisher · View at Google Scholar · View at Scopus - H. Shaverdi, M. R. Taha, and F. Kalantary, “Evolution of fabric in the shearing process based on micromechanics,” in
*Constitutive Modeling of Geomaterials*, Q. Yang, J. Zhang, H. Zheng, and Y. Yao, Eds., pp. 469–475, Springer, Berlin, Germany, 2013. View at Google Scholar - H. Matsuoka and H. Geka, “A stress-strain model for granular materials considering mechanism of fabric change,”
*Soils and Foundations*, vol. 23, no. 2, pp. 83–97, 1983. View at Google Scholar · View at Scopus