A new nonlinear unified strength (NUS) criterion is obtained based on the spatially mobilized plane (SMP) criterion and Mises criterion. New criterion is a series of smooth curves between SMP curved triangle and Mises circle in the π plane and thereby unifies the strength criteria. The new criterion can reflect the effect of the intermediate principal stress and consider the strength nonlinearity of a material. Based on the fabric tensor, the anisotropic parameter A is defined, and the anisotropic equation is proposed and introduced into the NUS criterion to form a nonlinear unified anisotropic strength criterion. The new criterion can be used to predict the strength variation of granular materials and cohesive materials under three-dimensional stress and can present the strength anisotropy of the geomaterials. The validity of the new criterion was checked using rock and soil materials. It is shown that the prediction results for the criterion agree well with the test data.
Natural Science Foundation of Ningbo2019A610394Ningbo University2140011540012Systematic Project of Guangxi Key Laboratory of Disaster Prevention and Structural Safety2019ZDK005Ningbo Public Welfare Science and Technology Planning Project2019C500121. Introduction
The research of strength criterion theory [1–11] is an important topic in geotechnical engineering. There are usually three types of methods for establishing strength criteria, namely, theoretical methods, empirical methods, and combined theoretical-empirical methods. Theoretical methods are usually based on a hypothesis, and there is a model to interpret the failure mechanism of geomaterials. All the model parameters in the strength criteria established by this method have clear physical meanings. This type of strength criterion uses the law of friction to explain the failure of geomaterials, including the Mohr–Coulomb strength theory, the D-P criterion [11], the spatially mobilized plane (SMP) criterion [2], and the twin shear unified strength theory [9, 12, 13]. Strength criteria established by empirical methods are usually based on the fitting of test data and include the famous Willam–Warnke criterion [14] and the Hoek–Brown criterion [6, 15]. Compared to the criteria that are based on theoretical methods, the physical meanings of some of the parameters of the strength criteria established by empirical methods are often not very clear. Finally, combined theoretical-empirical methods are usually used to establish special strength criteria for specific material properties. This type of criterion applies the Mohr–Coulomb strength theory to jointed rock materials [16] by introducing additional parameters [17–19].
However, the above strength criteria are associated with a single shape in the π plane and are thus unable to reflect the factors of material strength that may vary with the change of internal factors. To address this problem, researchers [10, 11, 20–26] have proposed to unify the strength criteria. The common method is to introduce some parameters to change the size of the shape function. The above criteria are all isotropic strength criteria and cannot depict the strength anisotropy of geomaterials. However, numerous geomaterials exhibit transverse isotropy due to deposition; that is, the strength of the material is the same in the deposition surface but different in the depositional surfaces of different directions. The directional shear test of most soils shows that the change in the direction of the principal stress also leads to strength anisotropy, reflecting the initial anisotropy of the material. Therefore, it is necessary to establish a uniform strength criterion that considers the depositional characteristics or initial anisotropy of geomaterials.
The above previous studies show that the strength criterion is researched through a simple single form to a unified strength theory system with a relatively wide scope of application. Consequently, it is important to explore a widely applicable strength theory system. In the present paper, a new anisotropic nonlinear unified strength criterion is established based on Mises criterion and SMP criterion. The new criterion is approximated as a series of smooth curves between the SMP curved triangle and the von Mises circle in the π plane to unify the strength criteria. Based on the fabric tensor, a nonlinear anisotropic strength criterion is proposed. The correctness of the criterion is verified using sand, clay, and rock materials.
2. Form of the NUS Criterion2.1. Linear Unified Strength (LUS) Criterion
SMP criterion (Matsuoka and Nakai [2]) has been widely used to check the strength of geomaterials because of its simple form and clear physical meaning. In the present study, the SMP criterion is modified to make it more widely applicable. The SMP criterion is shown in the following equation:(1)I1I2I3=C,where I1, I2, and I3 are stress invariants:(2)I1=σ1+σ2+σ3I2=σ1σ2+σ2σ3+σ1σ3I3=σ1σ2σ3.
C is a model parameter. Under conventional triaxial compression (σ2 = σ3), the sine of the friction angle φ0 of a granular material is(3)sinφ0=σ1−σ3σ1+σ3.
Substituting equations (2) and (3) into equation (1) gives(4)I13I3=3−sinφ031−sinφ0−sin2φ0+sin3φ0.
Using the trigonometric identity, the expression of q is obtained as(5)3p3−L2q2p+L3cos3θq3=0,where(6)L2=3+sin2φ012sin2φ0,(7)L3=9−sin2φ0108sin2φ0,(8)p=σ1+σ2+σ33,(9)q=12σ1−σ22+σ2−σ32+σ1−σ32,where p is the mean stress and q is generalized shear stress. Thus, the shape function of the SMP criterion is(10)g0θ=qq0=cos1/3arccos−33L3/2L23/2cos1/3arccos−33L3cos3θ/2L23/2,where q0 is the data for q when θ = 0. Therefore, the forms of p and q of the SMP criterion are(11)q=Mg0θp,where g0θ is obtained by equation (10) and M is the coefficient of friction:(12)M=6sinφ03−sinφ0.
The shape function of the SMP criterion cannot reflect the influence of internal factors on material; therefore, the application of the criterion is limited. In this report, the shape function is rewritten as the following expression:(13)q=Mgθp+σ0,where M is obtained from equation (12), σ0 which is the triaxial tensile strength of the material is given by(14)σ0=ccotφ0,and p and q are the average stress and the generalized shear stress, as defined by equations (9) and (10), respectively.
Therefore,(15)gθ=rθr0,where(16)rθ=31−α2L2cosδθ+6αsinφ03−sinφ0r0=31−α2L2cosδ0+6αsinφ03−sinφ0,and α is the model parameter determined by the strength ratio under the triaxial compression and extension conditions. The value of α ranges from 0 to 1. One has(17)δθ=13arccos−33L3cos3θ2L23/2δ0=13arccos−33L32L23/2,where L2 and L3 are obtained from equations (6) and (7), respectively.
Figure 1 illustrates the influences of parameters φ0 and α on linear criterion in the deviatoric plane. As Figure 1(b) shows, the shape of the criterion is shifted from the SMP curved triangle to the von Mises circle when α gradually increases from 0 to 1.
Variations of criterion with φ0 and α on the octahedral plane.
2.2. Nonlinear Isotropic Strength (NIS) Criterion
In the triaxial test of the rockfill material tests [27], the large confining pressure causes the granules to break [28–30], and the material exhibits remarkable strength nonlinearities. Therefore, to obtain the nonlinear of the rockfill material, the above new criterion equation (13) need be further modified as expressed in the following equation:(18a)q=Mfgθprp+σ0prn.
That is,(18b)lnqpr=nlnp+σ0pr+lnMfgθ,where pr is reference stress usually taken as pr = 101 kPa; Mf is usually obtained from the intercept of equation (18b) and the shape function; and n is a curvature parameter.
In the space of lnq/pr−lnp+σ0/pr, lnMfgθ is the intercept and n is the slope.
When n = 1,(19)Mf=M=6sinφ03−sinφ0,where n is the curvature parameter, and the shape function in the deviatoric plane is(20)gθ=31−α/2L2cosδθ+6αsinφ0/3−sinφ031−α/2L2cosδ0+6αsinφ0/3−sinφ0.
2.3. Nonlinear Unified Strength (NUS) Criterion
Numerous geomaterials exhibit transverse isotropy due to deposition, as shown in Figure 2; that is, the strength of the materials is the same within the depositional surface but different in the depositional surfaces of different directions, reflecting the initial anisotropy of the material. Therefore, it is necessary to establish a strength criterion that considers the depositional characteristics or initial anisotropy of the geomaterial. The premise of establishing the anisotropic strength criterion is to quantify the degree and direction of the initial anisotropy of the geomaterial. The fabric tensor proposed by Brewer is a good choice. The anisotropic function in the present study refers to the form in [31–33]:(21)fA=expη1A−A0+η2A−A02,where A is the anisotropic parameter and A0 is the value that A takes from b = 0 and α = 0. η1 and η2 in the anisotropic function are selected to depict the anisotropy of the material. The anisotropic function and the nonlinear isotropic strength criterion are combined to obtain an anisotropic strength criterion, as shown in the following equation:(22)q=fAMfgθprp+σ0prn.
Application of true triaxial tests on cross-anisotropic soils in the octahedral plane.
2.4. Effects of η1 and η2 on ANUS
The effect of isotropic parameters on the shape of the strength criterion in the deviatoric plane has been discussed above. This section focuses on the effect of anisotropic parameters on the shape of the strength criterion. According to the characteristics of the test data, the parameters η1 and η2 in the anisotropic function are selected to depict the anisotropy of the material. Different formulas can be selected for the anisotropic parameter A in equation (22) according to the true triaxial test or the hollow torsional shear test. In the present study, the anisotropic strength in the true triaxial test is mainly discussed.
In the true triaxial test, when the stress tensor and the fabric tensor are coaxial (β = 0), as shown in Figure 2, the shape of the strength criterion in the π plane varies with the parameters η1 and η2, as shown in Figures 3(a), 3(c), and 3(e). In this section, with η2 fixed at −0.25, −0.333, and 0, the effect of varying η1 on the strength curve in the π plane is discussed. As shown in Figures 3(a) and 3(b), with η2 = −0.25, the strength values of the anisotropic criterion and the isotropic criterion at θ = 180° are equal with different η1 values, and the criterion is symmetric about the σz axis. The criterion expands in the π plane when η1 > 0 and shrinks in the π plane when η1 < 0, and the degrees of expansion or shrinkage are related to the value of η1. The anisotropic function variation in the π plane is consistent with that of the strength criterion, as shown in Figure 3(b). With η2 = −0.333, as shown in Figure 3(c), the strength values of the anisotropic criterion and the isotropic criterion at θ = 120° and 240°, respectively, are equal with different η1 values, and the criterion is symmetric about the σz axis. Relative to the isotropic criterion, the anisotropic criterion expands in regions I and II and shrinks in region III of the π plane when η1 > 0, and it shrinks in regions I and II and expands in region III of the π plane when η1 < 0, with the degrees of expansion or shrinkage being related to the value of η1. The anisotropic function variation in the π plane is consistent with that of the strength criterion, as shown in Figure 3(d). When η2 = 0, as shown in Figure 3(e), the criterion at different η1 values is symmetric about the σz axis; when η1 > 0, the anisotropic criterion expands in the π plane relative to the isotropic criterion; and when η1 < 0, the anisotropic criterion shrinks in the π plane relative to the isotropic criterion. The degrees of expansion or shrinkage are related to the value of η1. The anisotropic function variation in the π plane is consistent with that of the strength criterion, as shown in Figure 3(f).
Effect of η1 on the failure loci and f(A) in the deviatoric plane with fixed η2 in the octahedral plane: (a, c, e) the failure curves and (b, d, f) the f(A).
When the stress tensor and the fabric tensor are not coaxial (β ≠ 0), the changes in the shape of the strength criterion in the π plane at different included angles are shown in Figures 4(a) and 4(c). With fixed η1 and η2 values (η1 = −0.08, η2 = 0), the effect of the change in the included angle between the deposition surface and the vertical stress on the strength curve in the π plane is discussed. When β = 0, as shown in Figure 4, the strength curve is symmetric about the σz axis in the π plane. When the deposition surface is rotated in the x−z plane, as shown in Figure 4(a), the strength curve is symmetric with respect to the σy axis in the π plane when β = 90°. The changes in the strength curve when β = 30° and 60° are shown in Figure 4(a). As the figure shows, when the deposition surface is rotated in the x−z plane, the strength values of the anisotropic strength criterion are equal at θ = 60° and 240°. The included angle has a large effect on the position and size of the strength curve. The method proposed in the present study can predict the trend and position of the strength curve at different included angles. When the deposition surface is rotated in the y−z plane, as shown in Figure 4(c), the strength curve is symmetric about the σx axis in the π plane when β = 90°. The changes in the strength curve when β = 30° and 60° are shown in Figure 4(c). As the figure indicates, when the deposition surface is rotated in the y−z plane, the strength values of the anisotropic strength criterion are equal at θ = 120° and 300°. The variations of the anisotropic function in the π plane in both cases are shown in Figure 4(b) and 4(d), where a similar situation with the change in the strength criterion is observed.
The variations of the failure surface and f(A) in the octahedral plane: (a, b) the y−z plane and (c, d) the x−z plane.
3. Test Verification3.1. Parameter Determination
To facilitate the application of the criteria, the model parameters should be determined as much as possible using conventional triaxial compression or extension tests. In this section, the proposed ANUS criterion is applied to various geomaterials. The geomaterial model parameters need to be determined only by conventional triaxial compression or extension tests. In this section, the strength parameters and anisotropy parameters in the anisotropic criterion are mainly determined. Callisto et al. [34] performed a large number of conventional triaxial tests and true triaxial tests on Pietrafitta clay to determine all the parameters in the anisotropic criterion. In the present study, the parameters are determined by conventional triaxial test data, and other test data are used for model verification. The specific determination steps are as follows:
Determination of Mf, n, and σ0. The p−q curve is drawn using the failure data of the conventional triaxial test, as shown in Figure 5(a). As shown in the figure, p−q exhibits a linear relationship. Therefore, n = 1, the slope of the straight line is M(M=6sinφ0/3−sinφ0), and the intercept is C (C=6ccosφ0/3−sinφ0), where c is the cohesion, and φ0 is the friction angle. From the fitted linear relation in the figure, it can be seen that, for the Pietrafitta clay, M = 0.5472 and C = 105.59 kPa. Therefore, it can be determined that c = 49.98 kPa and φ0 = 14.478°, and σ0 = ccotφ0 = 193.58 kPa can be obtained from the values of c and φ0, taking ξ = 0.
Determination of η1 and η2, η1, and η2 can be determined by the q value of other stress paths with b ≠ 0. Usually, the strength of triaxial extension is taken as the reference value for calibration. In most triaxial tests, only one of the terms in equation (21) is needed to satisfy the prediction accuracy requirement. For this reason, only the η1 term is taken in this test. In the triaxial extension test, when b = 1 and θ = 180°A–A0 = 4. Using qe=fBMfgθp+σ0, qe = 254.41 kPa (θ = 180° p = 250 kPa) can be obtained using the same method. Therefore, g(θ) = 0.8462, f(A) = exp(4η1), and η1 = 0.0417.
Using the above steps, all parameters of the anisotropic criterion for Pietrafitta clay are determined. The test parameter determination method is shown in Figure 5(a). The comparison of the prediction results by the isotropic criterion and the anisotropic criterion in Figure 5(b) shows that the anisotropic criterion proposed in the present study better predicts the test data.
Illustration on calibrating the parameters Mf, c, σ0 (b) Comparison of the two criteria for natural Pietrafitta clay (test data from [34]).
3.2. Strength Nonlinearity
A series of large triaxial tests for rockfill materials were conducted to research strength nonlinearity of this material [27]. In the meridional plane (p−q), the test results show the strength nonlinearities. Therefore, it is necessary to use nonlinear strength criteria to depict the strength development pattern. The two strength criteria are compared by using them to predict the test data, as shown in Figure 6. As the figure shows, the nonlinear strength criterion proposed this paper can depict the strength nonlinearity of this material.
Verifications of the two criteria. (a) Parameter determination. (b) Test data in the p−q diagram.
3.3. Anisotropy3.3.1. True Triaxial Test of Soil
In addition to the above true triaxial test of Pietrafitta clay, fine glass-bead sand is then used to verify the anisotropic criterion. Haruyama [35] conducted a large number of conventional triaxial tests and true triaxial tests to study the initial anisotropy of fine glass-bead sand. The effective stress of the test is p = 294 kPa. The conventional triaxial compression test and the triaxial tensile test are used to determine the parameters. The parameters determined by the test data are as follows: M = 1.142, σ0 = 0, n = 1, η1 = −0.038, and η2 = 0. The prediction results and test data in the deviatoric plane and b-φ plane are shown in Figure 7. It can be seen from the figure that the anisotropic criterion of the present study can reasonably predict the strength of the fine glass sand in the π plane. However, the isotropic criterion fails to predict the strength of regions II and III in the deviatoric plane and overestimates the material strength. In addition, in the b-φ plane, the isotropic criterion can only predict one variation curve of the friction angle, whereas the anisotropic criterion results in a continuous prediction curve in each region, reflecting the anisotropy of the material. Therefore, the anisotropic criterion can reasonably predict the peak strength of the fine glass-bead sand compared to the isotropic criterion.
Verifications of the two criteria for glass beads in (a) the octahedral plane; (b) the b–φ diagram (test data from [35]; this figure is reproduced from [36]).
3.3.2. Torsional Shear Tests for Soils
The failure friction angle measured from torsional shear tests on Leighton Buzzard sand of different densities and on spherical glass beads conducted by Yang et al. [37] was analyzed using the proposed ANUS criterion. The effective stress of the test is p = 200 kPa. The parameters determined by the test data are as follows: for glass ballotini, M = 1.43, σ0 = 0, n = 1, η1 = −0.368, and η2 = 0.1187; for medium dense LB sand, M = 1.805, σ0 = 0, n = 1, η1 = −0.2463, and η2 = 0.0593; for dense LB sand, M = 1.91, σ0 = 0, n = 1, η1 = −0.2668, and η2 = 0.0655. The prediction results (ANUS criterion and ALD criterion [38]) and test data in the α-φ plane are shown in Figure 8. As shown in Figure 8, under torsional shear tests, the sand shows obvious anisotropy, the new anisotropic criterion can predict the failure friction angle well, and the prediction curve is smooth.
Verifications of ANUS criterion for different granular materials (data from [37]).
3.3.3. Rock Triaxial Test considering Deposition Surface
This test is mainly used to verify the peak strength of the rock considering the different angles between the rock deposition surface and the principal stress. The isotropic strength criterion and the anisotropic strength criterion were used for prediction, and a comparison of the prediction results with the test data shows the superiority of the anisotropic strength criterion. Here, based on the schist triaxial tests with different deposition surface angles in the literature, the parameters of the anisotropic strength criterion are M = 2.36, σ0 = 8 Mpa, n = 0.74, η1 = −2.47, and η2 = −0.33. It can be seen from the test results that the two types of rock exhibit strong strength anisotropy in different deposition directions. In comparison, the use of the isotropic strength criterion to predict strength variations can lead to large errors, making it difficult to promote the isotropic criterion in engineering applications.
As shown in Figure 9, the anisotropic criterion can reasonably predict the strength failure line with different deposition surface angles in the p−q plane. In contrast, the isotropic strength criterion can only predict the failure line when β = 0°, cannot reflect the strength anisotropy, and overestimates the strength values with different deposition surface angles, especially for β = 45°. Meanwhile, the anisotropic criterion can reasonably predict the trend of strength variation under different confining pressures at different deposition surfaces, as shown in Figure 9(a). It can be seen from Figure 9(b) that the strengths of rocks with different deposition surface angles differ considerably. The peak strength of the rock is the largest when β = 0° or 90° and the smallest when the deposition surface angle is approximately 45°. In contrast, the isotropic criterion is a straight line in the β−q plane and, hence, cannot reflect the strength anisotropy caused by the direction of the deposition surface.
Comparison of the two criteria for Tournemire shale in (a) the p−q diagram and (b) the β−q diagram (test data from [39]; this figure is reproduced from [40]).
4. Conclusion
In this paper, a new strength criterion is proposed on the basis of the SMP criterion. In the π plane, the new criterion is a series of smooth curves and can unify the strength criteria.
The new criterion can depict the 3D strength variation of the material and reflect the effect of the intermediate principal stress and the nonlinear characteristics of the material strength.
The unified anisotropic strength criterion based on the pattern of the fabric evolution of granular materials can reflect the strength anisotropy caused by the deposition characteristics of the material. This criterion can be applied to the true triaxial test to consider the direction angle of the deposition surface. The correctness of the criterion was verified using sand, clay, and rock materials.
Data Availability
The underlying data used in the presented study were obtained from the literature.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was sponsored by Ningbo Natural Science Foundation Project (2019A610394) and the Initial Scientific Research Fund of Young Teachers in Ningbo University of Technology (2140011540012) and supported by the Systematic Project of Guangxi Key Laboratory of Disaster Prevention and Structural Safety (no. 2019ZDK005) and the Ningbo Public Welfare Science and Technology Planning Project (no. 2019C50012). These financial supports are gratefully acknowledged.
DruckerD. C.PragerW.GreenbergH. J.Extended limit design theorems for continuous media19749438138910.1090/qam/45573MatsuokaH.NakaiT.Stress-deformation and strength characteristics of soil under three different principal stresses19741974232597010.2208/jscej1969.1974.232_592-s2.0-0016352892MatsuokaH.On the significance of the “spatial mobilized plane”19761619110010.3208/sandf1972.16.912-s2.0-0016929665LadeP. V.DuncanJ. M.Elastoplastic stress-strain theory for cohesionless soil1975101GT1010371053LadeP. V.Elasto-plastic stress-strain theory for cohesionless soil with curved yield surfaces197713111019103510.1016/0020-7683(77)90073-72-s2.0-0017722509HoekE.BrownE.Empirical strength criterion for rock masses1980106910131035ChenW. F.Strength theory: applications. developments and prospects for 21st century1998Beijing, New YorkScience Press748LiuM. D.CarterJ. P.General strength criterion for geomaterials20033225325910.1061/(asce)1532-3641(2003)3:2(253)2-s2.0-24044506388YuM.-H.ZanY.-W.ZhaoJ.A unified strength criterion for rock material200239897598910.1016/s1365-1609(02)00097-72-s2.0-0036891643YoshimineY.LuD.ZhouA.ZouB.Generalized non-linear strength theory and transformed stress space200447669170910.1360/04ye01992-s2.0-26644445142XiaoY.LiuH.LiangR. Y.Modified Cam-Clay model incorporating unified nonlinear strength criterion201154480581010.1007/s11431-011-4313-42-s2.0-79954616079YuM.HeL.SongL.Twin shear stress theory and its generalizeation1985281111741183YuM.-h.Advances in strength theories for materials under complex stress state in the 20th Century200255316921810.1115/1.14724552-s2.0-0011700810WillamK. J.WarnkeE. P.Constitutive model for the triaxialbehavior of concrete19Proceedings of the International Association for Bridge and Structural Engineering, Seminar on Concrete Structure Subjected to Triaxial Stresses, ISMESMay 1975Bergamo, Italy130HoekE.BrownE. T.Practical estimates of rock mass strength19973481165118610.1016/s1365-1609(97)80069-x2-s2.0-0031432421MahendraS.BhawaniS.Modified Mohr-Coulomb criterion for non-linear triaxial and polyaxial strength of jointed rocks201251435210.1016/j.ijrmms.2011.12.0072-s2.0-84859507626GaoZ. W.ZhaoJ. D.YaoY. P.A generalized anisotropic failure criterion for geomaterials20124722–233166318510.1016/j.ijsolstr.2010.07.0162-s2.0-77956227334TongZ.FuP.ZhouS.DafaliasY. F.Experimental investigation of shear strength of sands with inherent fabric anisotropy20149225727510.1007/s11440-014-0303-62-s2.0-84896403971WangJ.SongZ.ZhaoB.LiuX.LiuJ.LaiJ.A study on the mechanical behavior and statistical damage constitutive model of sandstone201843105179519210.1007/s13369-017-3016-y2-s2.0-85050292141LiL.AubertinM.SimonR.BussièreB.Formulation and application of a general inelastic locus for geomaterials with variable porosity42260162310.1139/t04-1132-s2.0-24044520042BussièreSimonG.A new yield and failure criterion for geomaterials200858212513210.1680/geot.2008.58.2.1252-s2.0-41349094817MortaraG.A hierarchical single yield surface for frictional materials200936696096710.1016/j.compgeo.2009.03SuD.WangZ.-L.XingF.A two-parameter expression for failure surfaces200936351752410.1016/j.compgeo.2008.09.0012-s2.0-60849087331LiuM.GaoY.LiuH.A nonlinear Drucker-Prager and Matsuoka-Nakai unified failure criterion for geomaterials with separated stress invariants20125011010.1016/j.ijrmms.2012.01.0022-s2.0-84862786370LuD. C.MaC.DuX. L.JinL.GongQ. M.Development of a new nonlinear unified strength theory for geomaterials based on the characteristic stress concept20161720401605810.1061/(asce)gm.1943-5622.00007292-s2.0-85009812500OsswaldP. V.OsswaldT. A.A strength tensor based failure criterion with stress interactions20173982826283410.1002/pc.242752-s2.0-85008220545LiuM. C.GaoY. F.HuangX. M.Study on elasto-plastic constitutive model of rockfills with nonlinear strength characteristics2005273294298in ChineseVaradarajanA.SharmaK. G.AbbasS. M.DhawanA. K.Constitutive model for rockfill materials and determination of material constants20066422623710.1061/(asce)1532-3641(2006)6:4(226)2-s2.0-33746049163GuptaA. K.2000Delhi, IndiaIndian Institute of TechnologyLiuX. R.TuY. L.. ,WangP.ZhongZ. L.TangW. L.DuL. B.Particle breakage of soil-rock aggregate based on large-scale direct shear tests201739814251434in ChineseDafaliasY. F.PapadimitriouA. G.LiX. S.Sand plasticity model accounting for inherent fabric anisotropy2004130111319133310.1061/(asce)0733-9399(2004)130:11(1319)2-s2.0-10244229732PietruszczakS.MrozZ.Formulation of anisotropic failure criteria incorporating a microstructure tensor200026210511210.1016/s0266-352x(99)00034-82-s2.0-0033884159PietruszczakS.MrozZ.On failure criteria for anisotropic cohesive-frictional materials200125550952410.1002/nag.1412-s2.0-0035946388CallistoL.RampelloS.Shear strength and small-strain stiffness of a natural clay under general stress conditions200252854756010.1680/geot.2002.52.8.5472-s2.0-0036802509HaruyamaM.Anisotropic deformation-strength characteristics of an assembly of spherical particles under three dimensional stresses1981214415510.3208/sandf1972.21.4_412-s2.0-0019713876WangS.ZhongZ.LiuX.Development of an anisotropic nonlinear strength criterion for geomaterials based on SMP criterion20202030401918310.1061/(asce)gm.1943-5622.0001588YangL.-T.WanatowskiD.LiX.YuH.-S.CaiY.Laboratory and micromechanical investigation of soil anisotropy20162941141610.3208/jgssp.oth-122-s2.0-84971449469CaoW.WangR.ZhangJ. M.Formulation of anisotropic strength criteria for cohesionless granular materials20171771100401615110.1061/(asce)gm.1943-5622.00008612-s2.0-85018828966NiandouH.ShaoJ. F.HenryJ. P.FourmaintrauxD.Laboratory investigation of the mechanical behaviour of Tournemire shale199734131610.1016/s1365-1609(97)80029-92-s2.0-0030917493WangR.FuP.TongZ. X.Discussion of “description of inherent and induced anisotropy in granular media with particles of high sphericity”20181810701701410.1061/(asce)gm.1943-5622.0001026