In this contribution, a hypoplastic viscous soil model incorporating viscous and grain-inertia rate effects observed in granular flows is enhanced by means of the coordination number, which allows to characterize the evolution of the microstructure of the granular material. Near a critical solid volume fraction and critical coordination number there exists a transition from a solid-like (rate-independent) to a flow-like (rate-dependent) behavior of granular materials. This information is used to calibrate a so-linear concentration factor, which controls viscous and grain-inertial rate effects in terms of the rate dependent evolution of the coordination number. The influence of this proposed hypoplastic model considering the coordination number is investigated in this paper by means of a simple undrained shear test . The viscous hypoplastic model proposed by Guo et al. [1] is used as the basis for the hypoplastic model proposed in this work, which is characterized by an additive structure,where σ H and σ D are the static and dynamic contributions to the total Cauchy stress tensor σ, respectively. The dynamic partdepends on the second invariant of the rate of deformation gradient II d , on the deviatoric part of the rate of the deformation gradient d dev and on the viscous and grain-inertia scaling factorswhere µ is the dynamic viscosity of the interstitial fluid surrounding the grains, C is the mean solid volume fraction, ρ s the density of the solid grains, and d is the mean diameter. λ N is the proposed linear concentration factor given bywhereN max is the maximum coordination number corresponding to the maximum possible volume fraction C max , N * is the rate dependent coordination number N , which is computed using a Prandtl-type equation [2], and α, β are parameters that are calibrated according to experiments. It should be noted that in Bagnold's theory [3], the linear concentration factor is written as λ = Cmax C 1 3 − 1 −1 . Figure 1 a) displays λ N − C plots for different values of α (β = −1).
Coordination number and evolution of soil microstructureDiscrete Element Method (DEM) simulations of steady and unsteady shear flows performed by Vescovi et al. [4] show, that there exists a critical volume fraction C c and critical coordination number N c below which the granular material exhibits a rate dependent (flow-like) behavior , whereas a rate independent (solid-like) behavior is found above these critical values. A curve fitting (performed by the authors) of the N -C plots obtained from DEM simulations for steady shear flows [4] at different shearing rates is shown in Figure 1 b). Steady state values provided by DEM results are subsequently combined with the model proposed by Rothenburg et al. [5], describing the strain-driven evolution of the coordination number. This allows its determination for a given strain level. The proposed non-linear equation relating the shear strain γ and N is given bywhere c is a constant independent of the coordination number, N 0 is the initial value of the coordination num...