Research Article  Open Access
Weixin Dong, Liming Hu, Yu Zhen Yu, He Lv, "Comparison between Duncan and Chang’s EB Model and the Generalized Plasticity Model in the Analysis of a High EarthRockfill Dam", Journal of Applied Mathematics, vol. 2013, Article ID 709430, 12 pages, 2013. https://doi.org/10.1155/2013/709430
Comparison between Duncan and Chang’s EB Model and the Generalized Plasticity Model in the Analysis of a High EarthRockfill Dam
Abstract
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 stressstrain 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 earthrockfill 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 PZIII model can give a better description of deformation of the earthrockfill 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 stressstrain 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 [1].
Biot presented the famous threedimensional consolidation theory based on the effective stress theory, equilibrium equation, and continuity condition [2]. 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 stressstrain 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 earthrockfill 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 stressstrain characteristics of soils, such as shearinduced 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) [7] 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 [10] 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 [11–18]. 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 [21–24]. Furthermore, based on the framework of generalized plasticity theory [16], some limitations of the original model have been solved [25–28], such as pressure dependency, densification under cyclic loading. The details of the generalized plasticity theory and the original and proposed modified PastorZienkiewiczChan’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 [29–31]. The rockfill material has large particle size and sharp edges and corners, which can result in remarkable particle breakage and change the shearinduced 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 largescale finite element analyses of earth dams was less reported.
In this study, the original generalized plasticity model was modified to consider the stressstrain 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 earthrockfill 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 earthrockfill 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 [5] 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 [34] who proposed the hyperbolic stressstrain function below to describe the deviatoric stressaxial strain curve obtained from triaxial tests.
Consider in which and are model constants.
In this constitutive model, the tangential Young’s modulus and tangential bulk modulus are used to simulate the nonlinear elastic response of soils, which are assumed to be where is the atmospheric pressure, and are modulus numbers, and are exponents determining the rate of variation of moduli with confining pressure, and is the failure ratio with a invariable value less than 1.
The MohrCoulomb failure criterion is adopted in the model, and is a factor defined as shear stress level given by
In the unloading and reloading stage, the tangential Young’s modulus is defined as
So far, the model has 8 parameters, , , , , , , , . These parameters can be determined with a set of conventional triaxial tests.
In general, a curved MohrCoulomb failure envelop is adopted by setting and letting vary with confining pressure according to Then parameters and are replaced by 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) [16] 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 where and = elastic and plastic strain increments, respectively.
The relationship between strain and stress increments is expressed as where is the elastoplastic stiffness tensor given as where , , , and are elastic stiffness tensor, plastic flow direction vector, loading direction vector, and plastic modulus under loading or unloading conditions, respectively.
The loading direction vector is used to judge the loading and unloading conditions:
Then, the elastoplastic stiffness tensor can be obtained corresponding to the loading and unloading conditions.
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. PastorZienkiewiczChan Model
This model was presented by Pastor et al. [19]. The relationships between elastic volumetric and shear strain increments and stress increments are defined as where , are tangential bulk and shear moduli, respectively, and they are assumed to be where , , and are model parameters.
In order to determine the plastic stiffness tensor, variables , , and need to be defined. and are expressed as follows:
The dilatancy and stress ratio are related as follows: And has a similar expression as where , are model parameters and / is equal to relative density. If = , associated flow rule is used, otherwise nonassociated flow rule is used.
In the case of unloading, the unloading plastic flow direction vector is defined as
The loading plastic modulus is proposed as where limits the possible state and , accounts for phase transformation; considers soil degradation and is the accumulated plastic shear strain; accounts for past history and which is the mobilized stress function; and , , , are model parameters.
Under unloading condition, the plastic modulus is defined as respectively, where , are model parameters and is the stress ratio from which unloading takes place.
2.2.3. Modified Model
The PastorZienkiewiczChan model (PZIII 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 earthrockfill dams. The Original PZIII model has serious limitation in reflecting pressure dependency of soils.
Figure 1 shows the stressstrain relationships of a rockfill material under drained conventional triaxial tests using a set of parameters under different confining pressures, but PZIII model gives the same  curve, where , 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 PZIII model cannot be used to describe the mechanical behavior of rockfill dams.
(a)
(b)
Some relations of the original model are modified to take into account the influence of confining pressure as where and are elastic constants, and are model parameters to consider the effect of pressure dependency.
As sand behavior is dependent on densities or void ratio, a state pressure index, , proposed by Wang et al. [35] was introduced in the PZIII model and (13) was modified as where is a model parameter and in which is the mean pressure at critical state. The critical state line is given by
3. Nuozhadu Hydropower Project
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 earthrockfill 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 × 10^{8} m^{3}, with the normal storage water level of 812.5 m and dead water level of 765 m.
(a)
(b)
(c)
(d)
Figure 3 shows the material zoning and construction stages of the maximum crosssection. 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. Antiseepage 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.
(a)
(b)
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 PZIII 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 earthrockfill 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 PZIII model parameters in Table 2. As shown in Figures 4, 5, 6, 7, 8, and 9, the modified PZIII 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.


(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
5. ThreeDimensional 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 crosssection 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.
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 PZIII model, respectively.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
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. (1) 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 PZIII model.(2) The maximum settlement occurs in the middle of core wall due to lower modulus of clayey soil.(3) Because of the tremendous differences of modulus between rockfill material and clayey soil, there exists obvious arching effect in the core wall. (4) 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 PZIII model.(1) 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 PZIII 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.(2) 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 PZIII model. As we know, rockfill material is a typical kind of cohesionless coarsegrained 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 In Situ 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 crosssection. 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 PZIII 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.
(a)
(b)
(c)
As an elastoplastic model, the PZIII 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 PZIII model based on the generalized theory and original PastorZienkiewiczChan model to simulate the stressstrain 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 PZIII model, respectively. The simulations of triaxial stressstrain response show that the modified PZIII 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 earthrockfill 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 PZIII model. The comparisons of FEM results, and in situ monitoring data showed that the modified PZIII model can give a better description of deformation of the earthrockfill 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 (2012KY02 and 2013KY4).
References
 J. M. Duncan, “State of the art: limit equilibrium and finiteelement analysis of slopes,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 122, no. 7, pp. 577–596, 1996. View at: Google Scholar
 M. A. Biot, “General theory of threedimensional consolidation,” Journal of Applied Physics, vol. 12, no. 2, pp. 155–164, 1941. View at: Publisher Site  Google Scholar
 R. S. Sandhu and E. L. Wilson, “Finite element analysis of seepage in elastic media,” Journal of the Engineering Mechanics Division, vol. 95, no. 3, pp. 641–652, 1969. View at: Google Scholar
 J. T. Christian and J. W. Boehmer, “Plane strain consolidation by finite elements,” Journal of Soil Mechanics & Foundations Division, vol. 96, no. 4, pp. 1435–1457, 1970. View at: Google Scholar
 J. M. Duncan and C.Y. Chang, “Nonlinear analysis of stress and strain in soils,” Journal of the Soil Mechanics and Foundations Division, vol. 96, no. 5, pp. 1629–1653, 1970. View at: Google Scholar
 J. M. Duncan, P. M. Byrne, K. S. Wong, and P. Mabry, “Strength, stressstrain and bulk modulus parameters for finite element analyses of stresses and movements in soil masses,” Tech. Rep. UCB/GT/8001, University of California, Berkeley, Calif, USA, 1980. View at: Google Scholar
 D. C. Drucker, R. E. Gibson, and D. J. Henkel, “Soil mechanics and workhardening theories of plasticity,” Transactions of the American Society of Civil Engineers, vol. 122, pp. 338–346, 1957. View at: Google Scholar
 K. Roscoe, A. Schofield, and C. Wroth, “On the yielding of soils,” Geotechnique, vol. 8, no. 1, pp. 22–53, 1958. View at: Google Scholar
 K. Roscoe, A. Schofield, and A. Thurairajah, “Yielding of clays in states wetter than critical,” Geotechnique, vol. 13, no. 3, pp. 211–240, 1963. View at: Google Scholar
 J. Burland, “Correspondence on ‘The yielding and dilation of clay’,” Geotechnique, vol. 15, pp. 211–214, 1965. View at: Google Scholar
 P. V. Lade and J. M. Duncan, “Elastoplastic stressstrain theory for cohesionless soil,” Journal of the Geotechnical Engineering Division, vol. 101, no. 10, pp. 1037–1053, 1975. View at: Google Scholar
 I. S. Sandler, F. L. DiMaggio, and G. Y. Baladi, “Generalized cap model for geological materials,” Journal of the Geotechnical Engineering Division, vol. 102, no. 7, pp. 683–699, 1976. View at: Google Scholar
 X.S. Li, Y. F. Dafalias, and Z.L. Wang, “Statedependent dilatancy in criticalstate constitutive modelling of sand,” Canadian Geotechnical Journal, vol. 36, no. 4, pp. 599–611, 1999. View at: Google Scholar
 Y.P. Yao and D. Sun, “Application of Lade's criterion to Camclay model,” Journal of Engineering Mechanics, vol. 126, no. 1, pp. 112–119, 2000. View at: Google Scholar
 G. Y. Baladi and B. Rohani, “Elasticplastic model for saturated sand,” Journal of the Geotechnical Engineering Division, vol. 105, no. 4, pp. 465–480, 1979. View at: Google Scholar
 O. Zienkiewicz and Z. Mroz, “Generalized plasticity formulation and applications to geomechanics,” in Mechanics of Engineering Materials, C. S. Desai and R. H. Gallagher, Eds., pp. 655–679, John Wiley & Sons, New York, NY, USA, 1984. View at: Google Scholar
 C. S. Desai and M. O. Faruque, “Constitutive model for geological materials,” Journal of Engineering Mechanics, vol. 110, no. 9, pp. 1391–1408, 1984. View at: Google Scholar
 S. B. R. Murthy, A. Vatsala, and T. S. Nagaraj, “Revised Camclay model,” Journal of Geotechnical Engineering, vol. 117, no. 6, pp. 851–871, 1991. View at: Google Scholar
 M. Pastor, O. C. Zienkiewicz, and A. H. C. Chan, “Generalized plasticity and the modelling of soil behaviour,” International Journal for Numerical & Analytical Methods in Geomechanics, vol. 14, no. 3, pp. 151–190, 1990. View at: Google Scholar
 Z. Mroz and O. Zienkiewicz, “Uniform formulation of constitutive equations for clays and sands,” in Mechanics of Engineering Materials, C. S. Desai and R. H. Gallangher, Eds., pp. 415–449, John Wiley & Sons, New York, NY, USA, 1984. View at: Google Scholar
 G. Wang and J.M. Zhang, “Dynamic consolidation finite element analysis of a sedimentprotecting dyke under ocean wave loading,” Rock and Soil Mechanics, vol. 27, no. 4, pp. 555–560, 2006. View at: Google Scholar
 M. Alyami, M. Rouainia, and S. M. Wilkinson, “Numerical analysis of deformation behaviour of quay walls under earthquake loading,” Soil Dynamics and Earthquake Engineering, vol. 29, no. 3, pp. 525–536, 2009. View at: Publisher Site  Google Scholar
 H. Li, P. Manuel, and T. Li, “Application of an generalized plasticity model to ultrahigh rockfill dam,” in Proceedings of the 12th International Conference on Engineering, Science, Construction, and Operations in Challenging Environments—Earth and Space, pp. 385–398, Honolulu, Hawaii, USA, March 2010. View at: Publisher Site  Google Scholar
 T. Li and H. Zhang, “Dynamic parameter verification of PZ model and its application of dynamic analysis on rockfill dam,” in Proceedings of the 12th International Conference on Engineering, Science, Construction, and Operations in Challenging Environments—Earth and Space, pp. 2706–2713, Honolulu, Hawaii, USA, March 2010. View at: Publisher Site  Google Scholar
 M. Pastor, “A generalized plasticity model for anisotropic behaviour of sand,” Computer Methods and Advances in Geomechanics, vol. 1, pp. 661–668, 1991. View at: Google Scholar
 G. Bolzon, B. A. Schrefler, and O. C. Zienkiewicz, “Elastoplastic soil constitutive laws generalized to partially saturated states,” Geotechnique, vol. 46, no. 2, pp. 279–289, 1996. View at: Google Scholar
 H. I. Ling and H. Liu, “Pressurelevel dependency and densification behavior of sand through generalized plasticity model,” Journal of Engineering Mechanics, vol. 129, no. 8, pp. 851–860, 2003. View at: Publisher Site  Google Scholar
 H. I. Ling and S. Yang, “Unified sand model based on the critical state and generalized plasticity,” Journal of Engineering Mechanics, vol. 132, no. 12, pp. 1380–1391, 2006. View at: Publisher Site  Google Scholar
 N. D. Marschi, C. K. Chan, and H. B. Seed, “Evaluation of properties of rockfill materials,” Journal of the Soil Mechanics and Foundations Division, vol. 98, no. 1, pp. 95–114, 1972. View at: Google Scholar
 R. J. Marsal, “Large scale testing of rockfill materials,” Journal of the Soil Mechanics and Foundations Division, vol. 93, no. 2, pp. 27–43, 1967. View at: Google Scholar
 R. J. Marsal, “Mechanical properties of rockfill,” in Embankment Dam Engineering, pp. 109–200, John Wiley & Sons, New York, NY, USA, 1973. View at: Google Scholar
 P. V. Lade, J. A. Yamamuro, and P. A. Bopp, “Significance of particle crushing in granular materials,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 122, no. 4, pp. 309–316, 1996. View at: Google Scholar
 B. O. Hardin, “Crushing of soil particles,” Journal of Geotechnical Engineering, vol. 111, no. 10, pp. 1177–1192, 1985. View at: Google Scholar
 R. L. Kondner, “Hyperbolic stressstrain response: cohesive soils,” Journal of the Soil Mechanics and Foundations Division, vol. 89, no. 1, pp. 115–143, 1963. View at: Google Scholar
 Z.L. Wang, Y. F. Dafalias, X.S. Li, and F. I. Makdisi, “State pressure index for modeling sand behavior,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 128, no. 6, pp. 511–519, 2002. View at: Publisher Site  Google Scholar
Copyright
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.