A hypoplastic constitutive model for debris materials

Debris flow is a very common and destructive natural hazard in mountainous regions. Pore water pressure is the major triggering factor in the initiation of debris flow. Excessive pore water pressure is also observed during the runout and deposition of debris flow. Debris materials are normally treated as solid particle–viscous fluid mixture in the constitutive modeling. A suitable constitutive model which can capture the solid-like and fluid-like behavior of solid–fluid mixture should have the capability to describe the developing of pore water pressure (or effective stresses) in the initiation stage and determine the residual effective stresses exactly. In this paper, a constitutive model of debris materials is developed based on a framework where a static portion for the frictional behavior and a dynamic portion for the viscous behavior are combined. The frictional behavior is described by a hypoplastic model with critical state for granular materials. The model performance is demonstrated by simulating undrained simple shear tests of saturated sand, which are particularly relevant for the initiation of debris flows. The partial and full liquefaction of saturated granular material under undrained condition is reproduced by the hypoplastic model. The viscous behavior is described by the tensor form of a modified Bagnold’s theory for solid–fluid suspension, in which the drag force of the interstitial fluid and the particle collisions are considered. The complete model by combining the static and dynamic parts is used to simulate two annular shear tests. The predicted residual strength in the quasi-static stage combined with the stresses in the flowing stage agrees well with the experimental data. The non-quadratic dependence between the stresses and the shear rate in the slow shear stage for the relatively dense specimens is captured.


Introduction
Debris flow is a very common natural hazard in the mountainous areas of many countries. It represents the gravity-driven flow of a mixture of various sizes of sediment, water and air, down a steep slope, often initiated by heavy rainfall and landslides [17]. The highest velocity of debris flows can be more than 30 m/s; however, typical velocities are less than 10 m/s [24]. The fast debris flows may cause significant erosion, while increasing the sediment charge and destructive potential. Such mass flows cause serious casualties and property losses in many countries around the world. The initiation mechanisms of debris flow and the predicted possible velocity are essential information for the design of protective measures. Numerical analysis plays an important role to obtain this information, where a competent constitutive model for debris materials is required. The main factors influencing the initiation of debris flow are, among others, the topography, material parameters, water and the initial stress state in the affected slope [22]. Earth slopes with inclinations ranging from 26 to 45 have been generally identified as most prone to debris flow initiation [40]. The  volume fraction of debris materials, defined as the ratio between the solid volume and the total volume of a representative volume element, varies between about 30 and 65 %. The water from heavy rainfall or snow melting makes the unconsolidated superficial deposit on a steep hillside saturated, thereby leading to a reduced shear strength due to the decreasing of matric suction, and further triggering a landslide. Such an upland landslide may develop into a hillside debris flow when the water in the sliding mass cannot be discharged quickly and therefore gives rise to excessive pore water pressure. In this case, based on the principles of soil mechanics, the effective stresses between solid particles will decrease to cause the reduction or complete loss of shear strength. Upon initiation of debris flow, debris material shows fluid-like behavior. As concluded by Iverson [19], debris flow can be mobilized by three processes: (i) widespread Coulomb failure along a rupture surface within a saturated soil or sediment mass, (ii) partial or complete liquefaction of a sliding mass due to high pore-fluid pressure and (iii) conversion of landslide translational energy to internal vibrational energy. In these processes, the development of high pore water pressure is likely the most significant triggering factor. In addition, experimental observation [18] shows that an almost constant excess pore water pressure persists during the runout and depositing of debris flows. Thus, a suitable constitutive model which can capture the solid-like behavior before failure and the fluid-like behavior after failure should has the capability to describe the developing of pore water pressure (or effective stresses) in the initiation stage and determine the residual effective stresses exactly. Some important material parameters such as solid volume fraction (or void ratio in soil mechanics) and the internal friction coefficient need to be taken into account. Actually, debris materials are normally simplified as solid spherical particle-viscous fluid mixture and treated as a fluid continuum with microstructural effect in the constitutive modeling [10,11]. In most conventional models, constitutive equations for the static and dynamic regimes are formulated and applied separately, such as the models for the solid-like behaviors of granular materials [8,27,41,43] and that for the fluid-like behaviors [1,6,21]. Although some models for granular-fluid flows have taken the stress state of the quasi-static stage into account, the employed theories for the static regime, such as Mohr-Coulomb criterion [34] and extended von Mises yield criterion [32], still fail to determine the changing of pore water pressure from the deformation directly. Hypoplasticity was proposed as an alternative to plasticity for the description of solid-like behavior of granular materials [41,43]. The distinctive features of hypoplasticity are its simple formulation and capacity to capture some salient features of granular materials, such as non-linearity, dilatancy and yielding [42].
It may be the suitable choice for the description of solid-like behavior of debris materials. In this paper, a framework which consists of a static portion for the frictional behavior and a dynamic portion for the viscous behavior is introduced at first. Bagnold's constitutive model for a gravity-free suspension [1] is chosen as the dynamic portion in the framework. Then, the applicability of a specific hypoplastic model in the description of granular-fluid flows is studied by using this model to simulate the undrained simple shear test of saturated granular materials as shown in Fig. 1, which is in analogy to the initiation of a debris flow. The dynamic model, which is modified by fitting Bagnold's experimental data and taking a parameter termed critical solid volume fraction into account [15], is combined with the hypoplastic portion to obtain a new complete constitutive model for debris flows. The performance of the proposed model is demonstrated by some element tests in which the new model is used to simulate two annular shear tests with different materials and apparatus.

The framework of constitutive modeling for debris materials
As stated in the preceding section, debris materials show solid-like behavior before failure and fluid-like behavior after failure. This particular phenomenon cannot be modeled only within the framework of statics or dynamics. An applicable model may need to combine a static and a dynamic portion and make the transition from solid-like to fluid-like behavior turns out as an outcome [42]. In our former work [15], based on the velocity analysis of dry sand flow [4,26] and the force balance of an inclined plane supporting a uniform layer of sand-water mixture beneath a uniform layer of pure water [34], a framework for the constitutive model of debris materials was developed as the following form, Fig. 1 Schematic of undrained simple shear tests where P and T are the normal and shear stresses for the solid phase; P 0 and T 0 are the normal and shear stress caused by prolonged contact between particles; T v , T i and P i are slightly modified Bagnold's constitutive relations for a gravity-free dispersion of solid spheres sheared in Newtonian liquids. The stresses P 0 and T 0 are the static portion of the framework which satisfy a generalized Mohr-Coulomb type yield criterion [9,31]. Thus, where / denotes the residual friction angle after failure. They correspond to the residual stresses of debris materials in the quasi-static stage. For a simple shearing, the shear stress for the so-called macro-viscous regime, T v , has the following expression where the coefficient K 1 is related to the material property and expressed as U is the shear velocity as shown in Fig. 1 and dU=dy denotes the shear rate changing along the depth direction; C is the mean solid volume fraction and C c is the maximum solid volume fraction to assure a full shearing to occur; n is a fitting parameter and l is the dynamic viscosity of the interstitial fluid; k is a dimensionless parameter termed linear concentration. For perfectly spherical particles, k is defined as where s is the mean free distance between two particles; C 1 is the asymptotic limit of the maximum measured solid volume fraction as the container dimensions approach infinity, which is also related to the size of the particles [16]. The shear stress for the 'grain-inertia' regime, T i , is formulated as in which is also a coefficient related to the material property and is a correction factor based on the experimental results in [1]; q s and d denote the material density and mean diameter of the grains, respectively; the tangent of the angle a i corresponds to the ratio between the shear and normal stress in the 'grain-inertia' regime. Therefore, the expression of the normal stress in the 'grain-inertia' regime is T v , T i and P i are termed the dynamic portion of the framework (1). This framework implies that the contributions of contact friction, fluid viscosity and particle collisions coexist in the entire flow process. Bagnold's tests [1] for two different interstitial fluids with different viscosities but the same density show significant differences in the slow shear stage and tend to the same stress-strain relation when the shear velocity is large enough. In the rapid shear stage, the particle collisions become very fierce; the bulk behavior and dissipation of the flow kinetic energy are dominated by the inelastic and frictional particle collisions. An impact between two particles in a viscous liquid approximates a dry impact since the fluid effect is insignificant in comparison with the collision force in this stage [45]. Therefore, the linear term T v , rather than the quadratic terms T i and P i , makes the models based on the framework (1) capable of distinguishing granular-fluid mixtures with different interstitial fluid. The dry granular flow can be treated as a particular case where air is the interstitial fluid. For a free surface dry granular flow shown in Fig. 2, the viscous terms T v are normally much less than the residual strength T 0 in the beginning of the flow since the shear rate is very small in this stage. It is also negligible in the fast shearing stage since the viscous effect of air is insignificant compared to the frictional and collisional effect of particles. Thus, the framework (1) will be reduced to the following form in the case of free surface dry granular flow.
As stated in the literature [26], the relation (6) predicts a steady uniform flow only when the slope h is equal to the angle a i . However, experimental results [3] show that such steady flow can be obtained not only at a single slope but over a slope range. This experimental observation can be predicted by the reduced framework (10) [25,31]. According to force balance when a steady uniform flow is obtained in the free surface dry granular flow, we have where g is the gravity acceleration; h is the depth along the y axis which is normal to the flow bed. Then we get the stress ratio Let us assume that a i is greater than the residual friction angle /, which is consistent with the experimental observations of dry granular flows [30]. The normal stress P i is zero in the critical state of triggering the flow since the flow velocity is almost null at that time point. Thus, from (11) and (12), we obtain and where h 1 is the critical inclination for the granular material start flowing. With the increasing of inclination, another critical state will be reached. In this state, the component of gravity perpendicular to the flowing bed is totally supported by P i since the flow velocity is large enough at this inclination. From (12), we get where h 2 is the maximum inclination for equation (12) holds. It indicates that the framework (10), in which the stresses are divided into a static portion generated by prolonged contact of particles and a dynamic portion produced in particle collisions, can predict steady uniform flows over a slope range h 2 ½/; a i . By taking the effect of the interstitial fluid into account, a constitutive model developed within the complete framework (1) can describe not only dry granular flows but also granular-fluid flows.
In the above analysis, the simple formula for the initial value of P 0 , (13), is only applicable for free surface dry granular flows. As pointed out in the preceding section, debris materials are saturated solid-fluid mixtures which will be partially or fully liquefied in the initiation of debris flows. The normal stress P 0 is the effective stress and obtained by subtracting the excess pore water pressure from the total normal stress in this case. A proper theory is required to capture the partial or complete liquefaction, and further determine the residual strength P 0 and T 0 . As introduced before, hypoplasticity may be the suitable choice for the description of solid-like behavior of debris materials. In the following section, we study the capability of a specific hypoplastic model for capturing the main properties of debris materials in the quasi-static stage.
3 The applicability of hypoplastic models for debris materials Hypoplastic constitutive equations are based on nonlinear tensorial functions with the major advantages of simple formulation and few parameters. Two hypoplastic models, the one developed by Wu et al. [41] and the one by Gudehus [13], are compared in the selection of the static portion for the framework (1). In the more recent models by Gudehus [13], mainly the stiffness is modified by the two factors, f b and f e , which take into account the influence of stress state and density, respectively. In modeling debris flow, however, the strength is very important and the stiffness is not important. Moreover, his model makes use of the exponential functions for the dependence of critical void ratio and minimum void ratio on pressure. For each function the parameters reduce from 3 to 2. However, there are only few data in the literature for the exponential functions. Therefore, in this paper, we will embark on the model proposed by Wu et al. [41] which is the first hypoplastic model with critical state to verify that, by employing an appropriate hypoplastic model as the static portion, the combined model based on the framework (1) can fulfill an entire and quantitative description of stress state for debris materials from quasi-static stage to fast flow stage. It is worth mentioning that the hypoplastic model with critical state is just one of the choices for describing the initiation of debris flows. Recently some improved models have been available, e.g. [12,23,35], which are developed from some widely used versions of hypoplastic model [28,37] and aim to improve the dependence of stiffness on pressure and density. However, the capability of these models for capturing the phenomenon of liquefaction and the stability in the cases of large deformation or low confining pressure still need to be verified. A more concise hypoplastic model with the former mentioned capability and stability can be employed to determine the stress state in the quasi-static stage of debris materials.
The hypoplastic model with critical state is an improvement of a basic hypoplastic model for sand developed by Wu and Bauer [43] as where c i ði ¼ 1; . . .; 4Þ are dimensionless material parameters; T h and D denote the stress tensor and the strain rate tensor, respectively; T h Ã is the deviatoric stress tensor expressed by k D k¼ ffiffiffiffiffiffiffiffiffiffiffiffi ffi trðD 2 Þ p stands for the Euclidean norm and 1 is unit tensor. The Jaumann stress rate tensor T h in (16) is defined by where _ T h is the stress rate tensor (material time derivative of T h ); W denotes the spin tensor. The hypoplastic model (16) possesses simple mathematical formulation and contains only four material parameters, c 1 $ c 4 . The specific determination process of c 1 $ c 4 can be obtained in the literatures [5,41,43]. Two stress states, the initial hydrostatic and the state at failure, are chosen for the identification of c 1 $ c 4 based on a triaxial test with constant confining pressure, i. e. _ T h ð2; 2Þ ¼ _ T h ð3; 3Þ ¼ 0. And then, the following parameters are introduced: the stress ratio, R ¼ T h ð1; 1Þ=T h ð3; 3Þ; the initial tangent modulus, the initial Poisson ratio, t i ¼ ½Dð3; 3Þ=Dð1; 1Þ R¼1 ; the failure stress ratio, the failure Poisson ratio, t f ¼ ½Dð3; 3Þ=Dð1; 1Þ R¼R f : The failure stress ratio R f and the failure Poisson ratio t f are related to the friction angle / 0 and the dilatancy angle w, respectively, through the following relations [43]: and Taking the four material constants c 1 $ c 4 as unknowns, a system of four linear equations can be obtained by substituting the corresponded stress and strain rate of the two stress states into the model (16). Therefore, the material constants are determined as functions of the wellestablished parameters in soil mechanics, the initial tangent modulus E i , the initial Poisson ratio t i , the friction angle / 0 and the dilatancy angle w. It should be pointed out that these parameters are related to a specific confining pressure, all the sets of material constants used in this paper are obtained with a confining pressure T h ð3; 3Þ ¼ 100 kPa. In addition, the deviatoric loading in the initial hydrostatic state is considered to be zero, i. e. the initial Poisson ratio t i ¼ 0.
By taking the effect of void ratio and stress level into account, the model (16) was slightly modified to the following form [41]. where is a factor called density function. a is a material parameter related to the stress level and is the modified relative density; e is the void ratio; e min and e crt are the minimum and the critical void ratio, respectively. The effect of void ratio and stress level on the behavior of granular materials is taken into account in the model (21) by using the following expressions, and a ¼ q 1 þ q 2 expðq 3 j trT h jÞ ð25Þ where p i ði ¼ 1; . . .; 3Þ and q i ði ¼ 1; . . .; 3Þ are material parameters and can be determined by fitting the experimental data of drained triaxial tests under different confining pressure; j Á j denotes absolute value. It is shown that the model (21) is applicable to both initially and fully developed plastic deformation of granular materials with drained or undrained conditions [41,43]. It will reduce to the original one (16) when the void ratio e is equal to the critical value e crt from (22) and (23). It means, for same material, same constants c 1 $ c 4 will be obtained for the original and extended models in the case of e ¼ e ecrt . Thus, the material constants emerging in the model (21) can be determined by the same way as done for (16). The dilatancy angle w is equal to zero since there is no volume deformation in this case [44]. About the material parameters p i ði ¼ 1; . . .; 3Þ and q i ði ¼ 1; . . .; 3Þ, some theoretical and experimental analyses are presented in [41]. p 1 is the critical void ratio when the confining pressure approaches infinity, since p 3 is negative. The value of p 1 should be close to the minimum void ratio under a high confining pressure. For the case of zero confining pressure, the critical void ratio is equal to p 1 þ p 2 which may close to the maximum void ration measured with very low confining pressure. q 1 is assumed to be always equal to 1 and q 3 is a negative value. For the case of trT h ! 1, the difference between dense and loose packing tends to disappear since the parameter a ! 1. Based on the numerical parametric study [41], q 2 is suggested to lie in the range (-0.3, 0.0). p 3 and q 3 for quartz sand are assumed to be -0.0001 kPa. In the case of very low confining pressure, such as the state of liquefaction, relatively higher values of q 2 , p 3 and q 3 may be needed to keep the sensitivity of I e to the stress level. The hypoplastic model (21) may be a proper choice for describing the shear softening (liquefaction) and the residual strength in the beginning of a debris flow. During debris flow, the material is subjected to large shear deformation. For developing and evaluating constitutive models, the planar simple shear motion is particularly relevant [14]. Therefore, we try to verify the applicability of the hypoplastic model (21) in the simulation of debris flows initiation by using this model to reproduce the typical experimental results of granular materials in undrained simple shear tests. As presented in the literatures [7,46], saturated sand specimens with different initial void ratios demonstrate three types of stress-strain behavior in undrained simple shear tests as indicated in Fig. 3: (i) the dense specimens have tendency of dilation and show shear hardening to reach a ultimate steady state (USS) finally; (ii) the very loose specimens demonstrate shear softening to obtain constant residual strength or complete liquefaction in the critical steady state (CSS); (iii) the specimens with medium void ratio first soften, then harden and reach also a ultimate steady state [47]. The shear softening is considered to be the main mechanism in the mobilization of debris flows. Now we intend to reproduce these three types of stressstrain behavior in the element tests. The experimental results will be reproduced qualitatively rather than precisely, due to some important material parameters are not presented in the literature [46]. In order to obtain the material constants c 1 $ c 4 for sand in the critical state with I e ¼ 1, the initial tangent modulus E i is determined approximately by the following relation [20,36] in which P a is the atmospheric pressure (101.3 kPa) and r 33 is the effective confining stress, given as 100 kPa in the experiments. Thus, the initial tangent modulus is obtained approximately 15 MPa. The friction angle / 0 is assumed equal to a relatively small value, 25 , for saturated loose sand in the critical state with e ¼ e crt . Both the initial Poisson ratio t i and the dilatancy angle w are assumed to be 0 as stated before. The determined material constants for the model (21) are presented in Table 1.
The three types of stress-strain behavior are reproduced as shown in Fig. 4, when the values in Table 2 are employed for p i and q i in the relations (24) and (25).
It is indicated that, when the hardening arises, the increasing of the stress level reduces the critical void ratio e crt and increases the parameter a. Both changes increase the density function I e and then limit the developing of hardening. Conversely, when the softening occurs, I e will decrease to restrict softening and liquefaction. Due to the regulatory function of I e , the model (21) can describe the shear softening and the residual strength of very loose  Fig. 5, the normal stresses r ii ði ¼ 1; 2; 3Þ of the very loose specimen with e ¼ 0:876 tend to be isotropic when the shear strain is large enough, no matter what is the initial stress state. The isotropic normal stress in the large deformation stage corresponds to the former mentioned thermodynamic pressure P 0 .

A new constitutive model for debris materials
Based on the above analysis, the hypoplastic model (21) and the tensor form of the modified Bagnold's model are employed as the static and dynamic portions of the new constitutive model, respectively. The structure of the new model is proposed as In our former work about the constitutive model of granular-fluid flows [15], the three-dimensional form of the dynamic portion was obtained based on a simplest model structure for describing non-Newtonian fluid (see, for example, [38]). The general three-dimensional form of the dynamic portion is where Table 1 Material constants for the model (21) in the simulation of the experiments in [46] [-] [-] [-] -50.0 -629.6 -629.6 1220.8 Fig. 4 Simulation results of (21) for saturated sand with different initial void ratio in undrained simple shear tests: a shear strain versus shear stress, b mean principal stress versus shear stress Table 2 Parameters for e crt and a in the simulation of the experiments in [46]  is the strain rate deviator tensor and is the second invariant of D Ã . It can be easily shown that for a simple shear flow where the strain rate tensor takes the form D ¼ the dynamic stress (28) is reducible to the dynamic portion of the framework (1). From (27), the concrete model for debris materials is determined as where _ T h can be determined by (21).
The structure of the new model is demonstrated by simulating an undrained simple shearing flow. As shown in Fig. 6, a static portion obtained by the hypoplastic model (21) is combined with a dynamic portion to get the total effective stress.
What need to be mentioned is that the static portion, T h , is rate independent. It is varying due to the accumulation of the shear strain rather than the changing of the shear rate. By merging with the dynamic portion, the total effective stress (32) becomes rate dependent. As shown in Fig. 5, the normal stresses reach a residual constant when the shear strain is approximately 0.4. This process is finished with very small shear velocity in the so-called quasi-static stage. Thus, in the simulation, the shear rate must be kept at a small value before the failure of the granular-fluid mixture to make sure that the static portion is the dominant part in the total effective stress T. The static portion should be much greater than the dynamic portion at the point A in Fig. 6. One approach to meet this requirement in numerical calculations is using small shear strain acceleration and increasing the time steps for the stage before failure.
It is worth mentioning that Wu [42] developed a rate form framework by combining a hypoplastic model and a rate-dependent dynamic model as where T is the total Jaumann stress rate tensor and T d is the dynamic part of the Jaumann stress rate tensor. The models developed within this framework may have the capability to account for the different behaviors for loading and unloading. However, the Jaumann strain acceleration tensor, D, makes the implementation of these models in some numerical methods more difficult. It will be an interesting exploration to solve this problem in our future work.

Performance of the proposed model
In this section, the new model, (32), will be used to predict the stress-strain relations of granular-fluid flows with different materials and experimental apparatus in some element tests. The experimental data of two annular shear tests as undrained simple shear tests are employed to verify the applicability of the new model. In our former work [15], these two experiments are also simulated by a constitutive model which cannot capture the shear softening of granular-fluid materials in the quasi-static stage. The former simulation results can be used as a control group to highlight the function of the hypoplastic portion in the new model.

Dry granular materials
The experimental data of dry granular materials sheared in a annular shear cell were reported by Savage and Sayed [33]. The data for 1:0 mm spherical polystyrene beads are selected for the element tests. The loads applied by the upper disk range from 100 to 1500 N/m 2 which is normal to the flow surface. By checking the measured normal stress for 1:0 mm beads, we assume that the initial confining pressure of an element at the upper surface of the specimen has a value around 500 N/m 2 . The exact value of C 1 was not reported in the literature [33] and here is assumed equal to 0.64 which is a typical value for monosized spheres [2,16]. Thus, the corresponding minimum void ratio is determined equal to 0.563. The critical volume fraction C c is approximately 0.62 [34]. The internal friction angle / 0 of 1:0 mm spherical beads is 23 and the initial tangent modulus E i is assumed equal to 15 MPa as a typical value of loose granular materials with the confining pressure is 100 kPa. Based on the identification of parameters introduced in Sect. 3, the material constants c 1 $ c 4 and the parameters for the density function I e are determined and listed in Table 3.
The parameters for the dynamic portion are listed in Table 4.
As shown in Fig. 7, the predicted results are in good agreement with the experimental data of different solid volume fractions when some typical values are employed for the unstated parameters. The non-quadratic dependence between the stresses and strain rate in the slow shear stage for the samples with C ¼ 0:524 is captured by the new model. In the experiments [33], the shear velocity was adjusted to keep the height constant, thereby keeping the volume of the samples unchanged. It is equivalent to the undrained condition in the tests of saturated granular materials. The mean effective stress would decrease from the initial confining pressure to a residual value or zero to offset the tendency of volume compression in the quasistatic stage of very loose granular-fluid mixture. The residual normal and shear stresses, corresponding to the stresses P 0 and T 0 in the framework (1), are determined by the hypoplastic portion and presented in Table 5. Only the test with C ¼ 0:524 demonstrates residual strength. It is consistent with the experimental observation [46] that granular materials will be fully liquefied when the initial void ratio exceeds a threshold value. For the looser specimens where C ¼ 0:504, 0.483 and 0.461, the stress-strain rate curves in the rapid shear stage show a slope of about 2 in the logarithmic coordinates. It means the linear term T v which characterizes the effect of the interstitial fluid is insignificant while the quadratic law T i is dominant in this case as analyzed in the Sect. 2. It proves that the proposed model (32) can describe the shear softening of dry granular materials in the quasi-static stage and the stress-shear rate relation throughout the shear process from the quasi-static stage to the fast shearing stage.

Granular-water mixture
For the case of a granular-fluid mixture, we take Hanes and Inman's experiments [16] about spherical particles sheared in water as an example. The data for particles with diameter 1:85 mm which were stated as good quality ones are chosen to verify the new model. The maximum measured volume fraction for 1:85 mm particles was reported to be 0.55. Thus, the asymptotic limit C 1 is presumed to be approximately 0.61 and the minimum void ratio is 0.64. The critical volume fraction is assumed to be 0.52 since partial shearing was observed in the test of the specimen with C ¼ 0:53. The load from the upper disk is almost 500 N/m 2 . The internal friction angle / 0 is stated to be 28 and the initial tangent modulus E i is 15 MPa. The determined material constants c 1 $ c 4 and the parameters for the density function I e are presented in Table 6.
The parameters of the dynamic portion are listed in Table 7.
The simulation results are shown in Fig. 8. The stress states of the two specimens are reproduced based on the prediction of the residual stresses. The specific values are presented in Table 8. The sample with C ¼ 0:51 results in residual strength after failure in the undrained simple shearing. The shear stress-shear rate curves of C ¼ 0:49 has a slope less than 2 in the stage with shear rate between 1 and 10 in the logarithmic coordinates. Comparing to the dry granular flows, the effect of the interstitial fluid in a granular-fluid flow is non-negligible. The slight difference between the slopes of predicted curves and the experimental data for C ¼ 0:49 in the rapid shear stage implies a nonzero residual strength of this specimen. It may be attributed to that the employed parameters for the hypoplastic portion are not reasonable for this case.

Conclusions
In the initiation of debris flows, the development of excess pore water pressure is considered as the most significant triggering factor. Debris materials are normally simplified as granular-fluid mixture for constitutive modeling. Therefore, a theory which can be used to describe the solid-like behavior of debris materials should have the ability to capture the changing of pore water pressure. Moreover, a constitutive relation for the debris materials in the flowing stage should be rate dependent, in which some important material parameters in a granular-fluid flow, such as solid volume fraction, fluid viscosity and particle density, are taken into account. In a former developed framework [15], a static portion for the friction component and a dynamic portion for the viscous component are combined. The dynamic portion is composed of a linear term for drag force of the fluid and a quadratic term for the collisional force. For a dry granular flow on an inclined plane, the linear term is negligible since the viscous effect of air is insignificant compare to the frictional and collisional effect of particles. In this case, the model predicts a steady uniform flow over a slope range, which is consistent with the experimental observation [3]. The constitutive relations based on this framework can describe the stress state throughout the shear process from yielding to high-speed shearing. Moreover, a smooth transition is obtained between the so-called macro-viscous and grain-inertia regimes. The applicability of hypoplastic model in describing debris flows before failure is studied by simulating the undrained simple shear test of saturated granular material. Such test condition is particularly relevant to the initiation mechanism of debris flow. Three types of stressstrain behavior in which the 'liquefaction' is regarded as the main factor of debris flow mobilization are reproduced by the hypoplastic model. It is shown that the hypoplastic model has the capability to describe the changes of pore water pressure and further capture the shear softening and hardening behavior of granular-fluid mixtures. Therefore, it Then, this static part is combined with the tensor form of the modified Bagnold's dynamic model to obtain a new complete model for the modeling of debris materials from static to dynamic state. The new model is employed to simulate two annular shear tests with dry and water-saturated granular materials. In the case of dry granular flow with constant volume, the hypoplastic portion predicts that only the densest one of the four specimens has residual strength. It implies a non-quadratic dependence between the stresses and the shear rate in the slow flowing stage which was observed in the experiments. Similar conclusion is also obtained in the case of water-saturated granular flow.
Comparing to the dry granular flow, the linear term T v , which characterizes the effect of the interstitial fluid, is nonnegligible in the granular-fluid flow. The element test results show that the new model is applicable to the modeling of granular materials with different interstitial fluid. The predicted stress-strain curves agree well with the experimental data. Further verification is still needed for the new model. It is our intention to implement this model in some numerical codes for large deformation, such as SPH and computational fluid dynamics (CFD) codes, to simulate granularfluid flows in an inclined channel or a rotating drum. As mentioned before, a hypoplastic model developed by Wang and Wu [39] has been implemented in SPH for large deformation analysis [29]. Therefore, SPH will be the preferred choice for further verification of the new model. As mentioned before, the models in the rate form may have the capability to account for the different behaviors for loading and unloading. It will be an interesting exploration to develop the rate form expression for the dynamic portion in which the loading and unloading process can be distinguished.