JAM Journal of Applied Mathematics 1687-0042 1110-757X Hindawi Publishing Corporation 709430 10.1155/2013/709430 709430 Research Article Comparison between Duncan and Chang’s EB Model and the Generalized Plasticity Model in the Analysis of a High Earth-Rockfill Dam Dong Weixin Hu Liming Yu Yu Zhen Lv He Liang Fayun State Key Laboratory of Hydro-Science and Engineering Department of Hydraulic Engineering Tsinghua University Beijing 100084 China tsinghua.edu.cn 2013 5 10 2013 2013 04 06 2013 19 08 2013 20 08 2013 2013 Copyright © 2013 Weixin Dong 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.

Nonlinear elastic model and elastoplastic model are two main kinds of constitutive models of soil, which are widely used in the numerical analyses of soil structure. In this study, Duncan and Chang's EB model and the generalized plasticity model proposed by Pastor, Zienkiewicz, and Chan was discussed and applied to describe the stress-strain relationship of rockfill materials. The two models were validated using the results of triaxial shear tests under different confining pressures. The comparisons between the fittings of models and test data showed that the modified generalized plasticity model is capable of simulating the mechanical behaviours of rockfill materials. The modified generalized plasticity model was implemented into a finite element code to carry out static analyses of a high earth-rockfill dam in China. Nonlinear elastic analyses were also performed with Duncan and Chang's EB model in the same program framework. The comparisons of FEM results and in situ monitoring data showed that the modified PZ-III model can give a better description of deformation of the earth-rockfill dam than Duncan and Chang’s EB model.

1. Introduction

The constitutive model of soil is the keystone in the finite element analyses of geotechnical structures. A suitable constitutive model can simulate the stress-strain relationships of soils under static or dynamic conditions. Numerical analysis, especially for finite element method incorporated with soil constitutive models, has played a very important role in geotechnical analyses which always include complex boundary conditions, nonlinearity of material, and geometry .

Biot presented the famous three-dimensional consolidation theory based on the effective stress theory, equilibrium equation, and continuity condition . However, it is quite difficult to give the theoretical solution of Biot’s consolidation theory except for few simple problems. Up to the 1960s, with the rapid development of electronic computer and constitutive models of soils, Biot’s consolidation theory was successfully implemented in finite element codes to study the behavior of geotechnical structures [3, 4]. So far, thousands of constitutive models have been proposed, which can be mainly grouped in two categories: nonlinear elastic models and elastoplastic models.

For nonlinear elastic model, the nonlinear characteristic of soil stress-strain relationship is considered by sectionalized linearization. A typical nonlinear elastic model is Duncan and Chang’s Model [5, 6], which has been widely used in the numerical analyses of earth-rockfill dams, as the model parameters are quite easy to be determined from conventional triaxial tests. And, a lot of experience of application has been accumulated for this model. However, nonlinear elastic models also have some inherent limitations to represent the stress-strain characteristics of soils, such as shear-induced dilatancy and stress path dependency.

Elastoplastic models would be very adequate in describing many key features of soils. Classical elastoplastic models are based on the plastic incremental theory composed of yield condition, flow rule, and hardening law. In the 1950s, Drucker et al. (1957)  suggested a cap yield surface controlled by volumetric strain. Roscoe et al. [8, 9] proposed the concepts of critical state line and state boundary surface, and then they built the Original Cam Clay Model based on triaxial tests. Burland  suggested a different energy equation and then established the Modified Cam Clay Model. Since the establishment of Cam Clay Model, some other types of elastoplastic constitutive models have also achieved great development . Among these models, the generalized plasticity model [16, 19, 20] can simulate the static and dynamic mechanical behaviors of clays and sands. This model is very flexible and convenient to extend, as the complicated yield or plastic potential surfaces need not to be specified explicitly. And the model has been used successfully in the static or dynamic analyses of some geotechnical structures . Furthermore, based on the framework of generalized plasticity theory , some limitations of the original model have been solved , such as pressure dependency, densification under cyclic loading. The details of the generalized plasticity theory and the original and proposed modified Pastor-Zienkiewicz-Chan’s models will be introduced in the sections below.

However, little experience has as yet been accumulated in applying the generalized plasticity model to the simulation of rockfill materials. And we know that rockfill material is quite different from sands in mechanical properties . The rockfill material has large particle size and sharp edges and corners, which can result in remarkable particle breakage and change the shear-induced dilation [32, 33]. On the other hand, though the generalized plasticity model has gained great success in the modeling of soils, the application of this model in the large-scale finite element analyses of earth dams was less reported.

In this study, the original generalized plasticity model was modified to consider the stress-strain relationships of rockfill materials, as most of previous studies focused on sands and clays. Then, based on conventional triaxial test data, the model parameters for dam materials of the Nuozhadu high earth-rockfill dam in Southwest China are determined. Finally, the static simulation of this dam is carried out by using a finite element code incorporating with Duncan and Chang’s EB model and the modified generalize plasticity model. The comparison of numerical results and in situ monitoring data illustrates the advantages of modified generalized plasticity model in the simulation of earth-rockfill dams.

2. Constitutive Model Descriptions

Two constitutive models of soils were used in the finite element analyses. One is the Duncan and Chang’s EB model belonging to nonlinear elastic model, the other one is the generalized plasticity model.

2.1. Duncan and Chang’s Model

Duncan and Chang’s model  is a nonlinear elastic model, which has been widely used in the geotechnical engineering, especially in the numerical analyses of earth dams. It is attributed to Kondner  who proposed the hyperbolic stress-strain function below to describe the deviatoric stress-axial strain curve obtained from triaxial tests.

Consider (1)σ1-σ3=ε1a+bε1, in which a and b are model constants.

In this constitutive model, the tangential Young’s modulus Et and tangential bulk modulus Bt are used to simulate the nonlinear elastic response of soils, which are assumed to be (2)Et=KPa(σ3Pa)n(1-RfSl)2,Bt=KbPa(σ3Pa)m, where Pa is the atmospheric pressure, K and Kb are modulus numbers, n and m are exponents determining the rate of variation of moduli with confining pressure, and Rf is the failure ratio with a invariable value less than 1.

The Mohr-Coulomb failure criterion is adopted in the model, and Sl is a factor defined as shear stress level given by (3)Sl=(1-sinϕ)(σ1-σ3)2c·cosϕ+2σ3·sinϕ.

So far, the model has 8 parameters, c, φ, K, Kur, n, Rf, Kb, m. These parameters can be determined with a set of conventional triaxial tests.

In general, a curved Mohr-Coulomb failure envelop is adopted by setting c=0 and letting φ vary with confining pressure according to (5)φ=φ0-Δφlog(σ3Pa). Then parameters c and φ are replaced by φ0 and Δφ.

Although Duncan and Chang’s EB constitutive model is quite simple, it has gained significant success in geotechnical engineering. On one hand, it is easy to obtain the model parameters; on the other hand, much experience has been accumulated. Nevertheless, it cannot incorporate dilatancy which has an important influence in the mechanical behavior of soils. And furthermore, it can only consider unloading process in a crude way.

2.2. Generalized Plasticity Theory and Its Original Constitutive Model 2.2.1. Basic Theory

The generalized plasticity theory was proposed by Zienkiewicz and Mroz (1984)  to model the behaviors of sand under monotonic and cyclic loading. The key futures of this theory are that neither yield surface nor plastic potential surface needs to be defined explicitly, and consistency law is not required to determine plastic modulus. In the theory, the total strain increment is divided into elastic and plastic components.

Consider (6)dε=dεe+dεp, where dεe and dεp = elastic and plastic strain increments, respectively.

The relationship between strain and stress increments is expressed as (7)dσ=Dep:dε, where Dep is the elastoplastic stiffness tensor given as (8)Dep=De-De:ngL/U:nT:DeHL/U+nT:De:ngL/U, where De, ngL/U, n, and HL/U are elastic stiffness tensor, plastic flow direction vector, loading direction vector, and plastic modulus under loading or unloading conditions, respectively.

In the framework of generalized plasticity theory, all the components of the elastoplastic constitutive matrix are determined by the current state of stress and loading/unloading condition.

2.2.2. Pastor-Zienkiewicz-Chan Model

This model was presented by Pastor et al. . The relationships between elastic volumetric and shear strain increments and stress increments are defined as (10)dp=Kevdεve,dq=3Gesdεse, where Kev, Ges are tangential bulk and shear moduli, respectively, and they are assumed to be (11)Kev=Kesoppo,Ges=Gesoppo, where Keso, Geso, and po are model parameters.

In order to determine the plastic stiffness tensor, variables ngL/U, n, and HL/U need to be defined. ngL/U and n are expressed as follows: (12)ngL=(dg1+dg2,11+dg2)T,n=(df1+df2,11+df2)T.

The dilatancy dg and stress ratio η=q/p are related as follows: (13)dg=dεvpdεsp=(1+αg)(Mg-η). And df has a similar expression as (14)df=(1+αf)(Mf-η), where αf, αg are model parameters and Mg/Mf is equal to relative density. If df = dg, associated flow rule is used, otherwise nonassociated flow rule is used.

The loading plastic modulus HL is proposed as (16)HL=H0pHf(Hv+Hs)HDM, where Hf=(1-η/ηf)4 limits the possible state and ηf=(1+1/αf)Mf, Hv=1-η/Mg accounts for phase transformation; Hs=β0β1exp(-β0ξ) considers soil degradation and ξ is the accumulated plastic shear strain; HDM=(ςMAX/ς)γ accounts for past history and ς=p[1-αfη/(1+αf)/Mf]f(-1/α) which is the mobilized stress function; and H0, β0, β1, γ are model parameters.

Under unloading condition, the plastic modulus is defined as (17)HU=Hu0(Mgηu)γu,  Mgη>1,HU=Hu0,Mgη1, respectively, where Hu0, γu are model parameters and ηu is the stress ratio from which unloading takes place.

2.2.3. Modified Model

The Pastor-Zienkiewicz-Chan model (PZ-III for short) has gained considerable success in describing the behavior of sands and clays under monotonic and cyclic loadings. But it still has some shortcomings to predict the static or dynamic responds of sands, especially for rockfill materials which are widely used in earth-rockfill dams. The Original PZ-III model has serious limitation in reflecting pressure dependency of soils.

Figure 1 shows the stress-strain relationships of a rockfill material under drained conventional triaxial tests using a set of parameters under different confining pressures, but PZ-III model gives the same ε1-εv curve, where ε1, εv are axial strain and volumetric strain, respectively. As confining pressure ranges from 0 kPa to several MPa for a rockfill dam with height of 200–300 m, the original PZ-III model cannot be used to describe the mechanical behavior of rockfill dams.

Simulation of stress-strain relationships for Original PZ-III model.

Some relations of the original model are modified to take into account the influence of confining pressure as (18)Kev=Ke0pa(ppa)m,Ges=Ge0pa(ppa)n,HL=H0pa(ppa)mHf(Hv+Hs)HDM, where Ke0 and Ge0 are elastic constants, m and n are model parameters to consider the effect of pressure dependency.

As sand behavior is dependent on densities or void ratio, a state pressure index, Ip, proposed by Wang et al.  was introduced in the PZ-III model and (13) was modified as (19)dg=dεvpdεsp=(1+αg)(MgIpmp-η), where mp is a model parameter and Ip=p/pc in which pc is the mean pressure at critical state. The critical state line is given by (20)ec=Γ-λlog(pc).

Nuozhadu hydropower project is located in the Lancang River which is also named Mekong River in the downstream in Yunnan Province, Southwest China, as shown in Figure 2(a). The installed capacity of the powerstation is 5850 MW. The most important part of Nuozhadu hydropower project is the high earth-rockfill dam with a maximum height of 261.5 m, which is the highest one with the same type in China and the fourth highest in the world. The reservoir has a storage capacity of 237.0 × 108 m3, with the normal storage water level of 812.5 m and dead water level of 765 m.

Nuozhadu dam. (a) Nuozhadu dam location, (b) project blueprint, (c) Nuozhadu dam under construction, and (d) dam site geomorphology.

Figure 3 shows the material zoning and construction stages of the maximum cross-section. The elevation of the earth core bottom and the crest of the dam are 562.6 m and 824.1 m, respectively. The dam crest has a longitudinal length of 630 m with a width of 18 m. The upstream and downstream slopes are at 1.9 : 1 and 1.8 : 1, respectively. The dam body is composed of several different types of materials. The shells of upstream and downstream are composed of decomposed rock materials. Anti-seepage material in the earth core is clay mixed with gravel. Adding gravel to the clay can improve the strength of clay and reduce the arching effect between shells and earth core. The gravel material consists of fresh crushed stone of breccia and granite with a maximum diameter of 150 mm. In addition to these, the fine rockfill and filter materials are filled against the earth core to prevent the fine particle from being washed away.

The maximum cross-section. (a) Material zoning and (b) construction stage.

The dam construction was started in 2008 and was completed at the end of 2012. Figure 2(c) shows the dam under construction. Figure 3(b) demonstrates the practical construction process.

4. Experimental Validation of Model Parameters

The modified PZ-III model was implemented in a finite element code which has been successfully used to analyze earth dams with Duncan and Chang’s EB model and some other constitutive models. A set of triaxial test data was used to make sure that the model has been incorporated into the FEM code accurately.

The proposed generalized plasticity model totally needs 17 parameters. The model parameters used in the computation of the earth-rockfill dam were obtained by fitting the triaxial test results. Drained triaxial tests under different confining pressures were conducted to test the rockfill materials and mixed gravel clay, which are the main parts of the dam body.

Duncan and Chang’s EB model parameters are shown in Table 1 and the modified PZ-III model parameters in Table 2. As shown in Figures 4, 5, 6, 7, 8, and 9, the modified PZ-III model presents a better ability to simulate the mechanics behavior of rockfill materials and mixed gravel clay, especially for dilatancy. With the reduction of confining pressure, the rockfill materials tend to dilate as the experimental volumetric strain curve shows. Especially for the rockfill materials under low confining pressure, negative volumetric strain rapidly develops after a short stage of volumetric contraction. Due to the intrinsic limitation, Duncan and Chang’s EB model cannot simulate the dilatancy which is a crucial feature of rockfill materials.

Material parameters of Duncan and Chang’s EB model.

Material Rockfill I Rockfill II Mixed gravel clay
φ 55.82 54.33 39.30
Δφ 12.29 12.07 9.80
R f 0.73 0.74 0.77
K 1450 1360 520
K b 550 600 250
K u r 2800 2500 900
n 0.30 0.43 0.42
m 0.13 0.08 0.25

Material parameters of the modified PZ-III model.

Material Rockfill I Rockfill II Mixed gravel clay
K 0 500 1000 300
G 0 1500 3000 900
m 0.50 0.50 0.50
n 0.50 0.50 0.50
α f 0.45 0.45 0.45
α g 0.45 0.45 0.45
M f c 1.05 0.90 0.60
M g c 1.60 1.35 1.10
β 0 0.00 0.00 0.00
β 1 0.00 0.00 0.00
Γ 0.34 0.31 0.34
λ 0.10 0.09 0.03
m p 0.35 0.40 0.0
H 0 800 1200 900
γ 5 5 5
γ u 5 5 5
H u 0 /MPa 9 9 10

Comparison between fittings of Duncan and Chang’s EB model and experimental triaxial tests results for rockfill material I.

Comparison between fittings of the modified PZ-III model and experimental triaxial tests results for rockfill material I.

Comparison between fittings of Duncan and Chang’s EB model and experimental triaxial tests results for rockfill material II.

Comparison between fittings of the modified PZ-III model and experimental triaxial tests results for rockfill material II.

Comparison between fittings of Duncan and Chang’s EB model and experimental triaxial tests results for clay.

Comparison between fittings of the modified PZ-III model and experimental triaxial tests results for clay.

5. Three-Dimensional Finite Element Analyses 5.1. Computation Model

The numerical analyses were performed to simulate the performance of the dam during construction and impounding periods with effective stress finite element analysis.

First, the 2D finite element mesh of the maximum cross-section of the dam was discretized according to the material zoning and construction design (see Figure 3). Then, the 2D mesh was extended to 3D mesh in accordance with contour line of the river valley. Figure 10 shows the 3D mesh of the Nuozhadu dam with 8095 brick and degenerated brick elements and 8340 nodes.

3D FEM mesh of Nuozhadu dam.

The numerical simulations contain two stages, filling and impounding. During the filling stage, the dam body mainly subjects to body weight. Then, at the end of construction, upstream water level goes up to the normal storage water level. The interaction between pore water and soil skeleton was considered through the whole numerical computation.

5.2. Results and Analyses 5.2.1. Numerical Results Analyses

Figures 11 and 12 show the numerical results of finite element analyses with Duncan and Chang’s EB model and the modified PZ-III model, respectively.

Displacement and stress contour of the maximum section for Duncan and Chang’s EB model: (a) displacement along river (m), (b) vertical displacement (m), (c) major principle stress (MPa), and (d) minor principle stress (MPa).

Displacement and stress contour of the maximum section for the modified PZ-III model: (a) displacement along river (m), (b) vertical displacement (m), (c) major principle stress (MPa), and (d) minor principle stress (MPa).

Through the comparison and analysis of the numerical results (Figures 11 and 12), we can find some similarities and differences for these two models.

On one hand, we can see many similar places in the distributions of displacements and stresses.

After the reservoir impounding, due to the huge water pressure on upstream dam, horizontal displacement develops toward the downstream, and the largest displacement is about 1.05 m for EB model and 0.74 m for modified PZ-III model.

The maximum settlement occurs in the middle of core wall due to lower modulus of clayey soil.

Because of the tremendous differences of modulus between rockfill material and clayey soil, there exists obvious arching effect in the core wall.

Effective stress in upstream shell is less than the downstream shell due to the water pressure in the upstream shell.

On the other hand, some differences also exist, which illustrate the advantages of modified PZ-III model.

After the reservoir is impounded, upward displacement as large as 0.7 m (see Figure 11(b)) develops on the upstream shell near dam crest for EB model and nearly 0 m for modified PZ-III model (see Figure 12(b)). In fact, monitoring data of practical engineering projects shows that no large upward displacement happened after impounding. This is due to its weakness of EB model to distinguish the loading and unloading condition during the water impounding.

In the distribution of minor principle stress (Figures 11(d) and 12(d)), negative stress (i.e., tensile stress) occurs in the upstream shell for EB model, whereas very little tensile stress exists for modified PZ-III model. As we know, rockfill material is a typical kind of cohesionless coarse-grained soil, which means that it has no tensile strength. Therefore, the existence of large area of tensile stress in the upstream shell is unreasonable.

5.2.2. Comparison between Numerical and <italic>In Situ</italic> Monitoring Data

Settlement is a key indicator to assess the safety of an earth dam. Figures 13 and 14 show the in situ monitoring data and FEM results of settlement in the maximum cross-section. The in situ data were obtained from electromagnetism type settlement gauges which were embedded during construction in the dam (as shown in Figure 3(a)). Through the comparisons of in situ monitoring and numerical results, we can see that the modified PZ-III model gave a better prediction than the EB model. However, as deformation induced by wetting of rockfill materials was not considered, the FEM result of settlement was below than the in situ monitoring data.

Comparison between in situ monitoring settlement and FEM results.

Comparison between in situ monitoring settlement and FEM results.

As an elastoplastic model, the PZ-III model is capable of representing the mechanical behavior of soils better than nonlinear elastic model such as Duncan and Chang’s EB model. And the above finite element analyses also proved it.

6. Conclusions

This paper presents a modified PZ-III model based on the generalized theory and original Pastor-Zienkiewicz-Chan model to simulate the stress-strain relationship of rockfill materials.

Triaxial test results of the filling materials of Nuozhadu dam were used to validate the proposed model and determine the model parameters of Duncan and Chang’s EB model and the modified PZ-III model, respectively. The simulations of triaxial stress-strain response show that the modified PZ-III model is capable of representing the key features of cohesionless soil, such as nonlinearity, dilatancy, and pressure dependency.

The proposed model has been incorporated into a finite element code to simulate the static response of a high earth-rockfill dam in China. The results were compared with those of Duncan and Chang’s EB model. The two set of results have both similarities and differences and the differences illustrate the advantages of the modified PZ-III model. The comparisons of FEM results, and in situ monitoring data showed that the modified PZ-III model can give a better description of deformation of the earth-rockfill dam than Duncan and Chang’s EB model.

Acknowledgments

This work was supported by the National Nature Science Foundation of China (51179092) and the State Key Laboratory of Hydroscience and Engineering Project (2012-KY-02 and 2013-KY-4).

Duncan J. M. State of the art: limit equilibrium and finite-element analysis of slopes Journal of Geotechnical and Geoenvironmental Engineering 1996 122 7 577 596 2-s2.0-0001690750 Biot M. A. General theory of three-dimensional consolidation Journal of Applied Physics 1941 12 2 155 164 2-s2.0-0345643551 10.1063/1.1712886 Sandhu R. S. Wilson E. L. Finite element analysis of seepage in elastic media Journal of the Engineering Mechanics Division 1969 95 3 641 652 Christian J. T. Boehmer J. W. Plane strain consolidation by finite elements Journal of Soil Mechanics & Foundations Division 1970 96 4 1435 1457 Duncan J. M. Chang C.-Y. Nonlinear analysis of stress and strain in soils Journal of the Soil Mechanics and Foundations Division 1970 96 5 1629 1653 Duncan J. M. Byrne P. M. Wong K. S. Mabry P. Strength, stress-strain and bulk modulus parameters for finite element analyses of stresses and movements in soil masses 1980 UCB/GT/80-01 Berkeley, Calif, USA University of California Drucker D. C. Gibson R. E. Henkel D. J. Soil mechanics and work-hardening theories of plasticity Transactions of the American Society of Civil Engineers 1957 122 338 346 Roscoe K. Schofield A. Wroth C. On the yielding of soils Geotechnique 1958 8 1 22 53 Roscoe K. Schofield A. Thurairajah A. Yielding of clays in states wetter than critical Geotechnique 1963 13 3 211 240 Burland J. Correspondence on ‘The yielding and dilation of clay’ Geotechnique 1965 15 211 214 Lade P. V. Duncan J. M. Elastoplastic stress-strain theory for cohesionless soil Journal of the Geotechnical Engineering Division 1975 101 10 1037 1053 2-s2.0-0016568353 Sandler I. S. DiMaggio F. L. Baladi G. Y. Generalized cap model for geological materials Journal of the Geotechnical Engineering Division 1976 102 7 683 699 2-s2.0-0016974674 Li X.-S. Dafalias Y. F. Wang Z.-L. State-dependent dilatancy in critical-state constitutive modelling of sand Canadian Geotechnical Journal 1999 36 4 599 611 2-s2.0-0033237201 Yao Y.-P. Sun D. Application of Lade's criterion to Cam-clay model Journal of Engineering Mechanics 2000 126 1 112 119 2-s2.0-0342472149 Baladi G. Y. Rohani B. Elastic-plastic model for saturated sand Journal of the Geotechnical Engineering Division 1979 105 4 465 480 2-s2.0-0018455636 Zienkiewicz O. Mroz Z. Desai C. S. Gallagher R. H. Generalized plasticity formulation and applications to geomechanics Mechanics of Engineering Materials 1984 New York, NY, USA John Wiley & Sons 655 679 Desai C. S. Faruque M. O. Constitutive model for geological materials Journal of Engineering Mechanics 1984 110 9 1391 1408 2-s2.0-0021643829 Murthy S. B. R. Vatsala A. Nagaraj T. S. Revised Cam-clay model Journal of Geotechnical Engineering 1991 117 6 851 871 2-s2.0-0026173736 Pastor M. Zienkiewicz O. C. Chan A. H. C. Generalized plasticity and the modelling of soil behaviour International Journal for Numerical & Analytical Methods in Geomechanics 1990 14 3 151 190 2-s2.0-0025573626 Mroz Z. Zienkiewicz O. Desai C. S. Gallangher R. H. Uniform formulation of constitutive equations for clays and sands Mechanics of Engineering Materials 1984 New York, NY, USA John Wiley & Sons 415 449 Wang G. Zhang J.-M. Dynamic consolidation finite element analysis of a sediment-protecting dyke under ocean wave loading Rock and Soil Mechanics 2006 27 4 555 560 2-s2.0-33646384961 Alyami M. Rouainia M. Wilkinson S. M. Numerical analysis of deformation behaviour of quay walls under earthquake loading Soil Dynamics and Earthquake Engineering 2009 29 3 525 536 2-s2.0-58649083204 10.1016/j.soildyn.2008.06.004 Li H. Manuel P. Li T. Application of an generalized plasticity model to ultra-high rockfill dam Proceedings of the 12th International Conference on Engineering, Science, Construction, and Operations in Challenging Environments—Earth and Space March 2010 Honolulu, Hawaii, USA 385 398 2-s2.0-78649247993 10.1061/41096(366)40 Li T. Zhang H. Dynamic parameter verification of P-Z model and its application of dynamic analysis on rockfill dam Proceedings of the 12th International Conference on Engineering, Science, Construction, and Operations in Challenging Environments—Earth and Space March 2010 Honolulu, Hawaii, USA 2706 2713 2-s2.0-78649241014 10.1061/41096(366)255 Pastor M. A generalized plasticity model for anisotropic behaviour of sand Computer Methods and Advances in Geomechanics 1991 1 661 668 2-s2.0-0026305315 Bolzon G. Schrefler B. A. Zienkiewicz O. C. Elastoplastic soil constitutive laws generalized to partially saturated states Geotechnique 1996 46 2 279 289 2-s2.0-0030165092 Ling H. I. Liu H. Pressure-level dependency and densification behavior of sand through generalized plasticity model Journal of Engineering Mechanics 2003 129 8 851 860 2-s2.0-0042823994 10.1061/(ASCE)0733-9399(2003)129:8(851) Ling H. I. Yang S. Unified sand model based on the critical state and generalized plasticity Journal of Engineering Mechanics 2006 132 12 1380 1391 2-s2.0-33751106660 10.1061/(ASCE)0733-9399(2006)132:12(1380) Marschi N. D. Chan C. K. Seed H. B. Evaluation of properties of rockfill materials Journal of the Soil Mechanics and Foundations Division 1972 98 1 95 114 Marsal R. J. Large scale testing of rockfill materials Journal of the Soil Mechanics and Foundations Division 1967 93 2 27 43 Marsal R. J. Mechanical properties of rockfill Embankment Dam Engineering 1973 New York, NY, USA John Wiley & Sons 109 200 Lade P. V. Yamamuro J. A. Bopp P. A. Significance of particle crushing in granular materials Journal of Geotechnical and Geoenvironmental Engineering 1996 122 4 309 316 2-s2.0-0030390147 Hardin B. O. Crushing of soil particles Journal of Geotechnical Engineering 1985 111 10 1177 1192 2-s2.0-0022267649 Kondner R. L. Hyperbolic stress-strain response: cohesive soils Journal of the Soil Mechanics and Foundations Division 1963 89 1 115 143 Wang Z.-L. Dafalias Y. F. Li X.-S. Makdisi F. I. State pressure index for modeling sand behavior Journal of Geotechnical and Geoenvironmental Engineering 2002 128 6 511 519 2-s2.0-0036606489 10.1061/(ASCE)1090-0241(2002)128:6(511)