Search

Based on the meso-structured Voronoi cell model and the meso-macro homogenization procedure between the discrete particle assembly and the porous continuum for wet granular materials, meso-hydro-mechanically informed effective pressure and effective stress measures for saturated and unsaturated porous media are defined. The meso-hydro-mechanically informed generalized effective stress for saturated porous continua taking into account the volumetric deformation of solid grains due to pore liquid pressure is derived. The Biot coefficient associated to the meso-hydro- mechanically informed generalized effective stress for saturated porous media is formulated. The differences of the definitions for proposed generalized effective stress and Biot coefficient compared with those defined in the generalized Biot theory of saturated porous continua and in averaging theories are discussed. The wet meso-structured Voronoi cell model, consisting of three immiscible and interrelated (i.e. solid grains, interstitial liquid and gas) phases, at low bulk saturation (below about 30%) is proposed. A meso-structural pattern with the binary bond mode of pendular liquid bridges is assumed in particular to derive the meso-hydro- mechanically informed macroscopic anisotropic effective pressure and effective stress tensors for unsaturated porous media. As the isotropic case of the wet meso-structured Voronoi cell model is considered, the meso-hydro-mechanically informed effective pressure tensor degrades to the scalar variable in the same form as in the theory of macroscopic unsaturated porous continua. The proposed meso-hydro-mechanically informed Bishop's parameter is derived and obtained as a function of saturation, porosity, and meso-structural parameters, without need to introduce any macroscopic phenomenological assumptions for the description of hydro-mechanical constitutive behavior.

Based on the meso-structured Voronoi cell model for discrete particle assembly and the derived meso-mechanically informed constitutive relations of anisotropic Cosserat continuum, thermodynamic framework of isothermal meso-mechanically informed damage-healing and plastic process for granular materials is presented. The accumulated net (effective) damage factor tensor combining both material damage and healing effects is defined in terms of the initial (undamaged) and current (damaged) elastic moduli tensors of the meso-structured Voronoi cell attributed to the material point. According to the non-negativity of thermodynamic energy dissipations, the net damage variable is separated into the two component internal state variables; i.e. the damage and healing variables, which are accumulated in terms of incremental damage and healing variables, respectively. The meso-mechanically informed macroscopic damage-healing and plastic characterization are achieved without the need to specify macroscopic phenomenological damage, healing and plastic criteria, and their evolution laws. The merit of the proposed tensorial net damage and healing variables in modeling healing effects on initial weakened elastic stiffness (i.e. initial material defects) is demonstrated in terms of their isotropic scalar forms and integrated into the continuum damage-healing mechanics. The numerical results conceptually illustrate the performance of the proposed definitions of meso-mechanically informed net damage, damage, and healing variables. The coupled damage-healing and plastic process in anisotropic Cosserat continuum for granular materials is characterized in terms of densities of thermodynamic dissipations that make effects of the damage-healing and the plastic component processes on the material failure quantitatively comparable.

A mixed finite element (FE) procedure of the gradient Cosserat continuum for the second-order computational homogenisation of granular materials is presented. The proposed mixed FE is developed based on the Hu---Washizu variational principle. Translational displacements, microrotations, and displacement gradients with Lagrange multipliers are taken as the independent nodal variables. The tangent stiffness matrix of the mixed FE is formulated. The advantage of the gradient Cosserat continuum model in capturing the meso-structural size effect is numerically demonstrated. Patch tests are specially designed and performed to validate the mixed FE formulations. A numerical example is presented to demonstrate the performance of the mixed FE procedure in the simulation of strain softening and localisation phenomena, while without the need to specify the macroscopic phenomenological constitutive relationship and material failure model. The meso-structural mechanisms of the macroscopic failure of granular materials are detected, i.e. significant development of dissipative sliding and rolling frictions among particles in contacts, resulting in the loss of contacts.

The microstructures of granular materials are represented by Voronoi cells generated with a Voronoi tessellation of discrete particle assembly. A Voronoi cell model including not only the reference particle laid inside the Voronoi cell but also its intermediate neighboring particles is presented to formulate micromechanically based macroscopic constitutive relations and constitutive modular tensors of effective Cosserat continuum. The anisotropy of effective Cosserat continuum due to intrinsic characters and deformation-induced evolutions of microstructure of granular materials of the Voronoi cell is quantitatively demonstrated. The derived micromechanically informed macroscopic constitutive relation of effective Cosserat continuum reveals that the Cauchy stresses are not only constitutively related to the strains but also to the curvatures defined in Cosserat continuum, likewise, the couple stresses are not only constitutively related to the curvatures but also to the strains. The derived modular tensors are verified by comparisons of them with those given for classical isotropic Cosserat continuum and are used to identify the elastic constitutive parameters of isotropic Cosserat continuum. The micromechanically informed macroscopic damage factor tensor to characterize anisotropic material damage of effective Cosserat continuum is formulated with no need specifying macroscopic phenomenological damage criterion and damage evolution law. The principal directions and values of the derived damage factor tensor along with numerical results reveal the microscopic mechanisms of macroscopic damage phenomenon, i.e. loss of contacts, re-orientation of contacts of the reference particle with its intermediate neighboring particles and concomitant volumetirc dilatation of the Voronoi cell.

The paper reports recent advances in multiscale modeling of granular materials, particularly in the second-order computational homogenization method and corresponding global–local mixed FEM-DEM nested analysis scheme. The gradient Cosserat continuum and the classical Cosserat continuum are assumed for modelling granular media at the macro- and meso- scales, respectively. According to the generalized Hill’s lemma formulated for the adopted meso-macro continuum modeling the non-uniform macroscopic strain field with macroscopic strain gradients is downscaled to each representative volume element (RVE), while satisfaction of the generalized Hill–Mandel condition is ensured. The advantage of the gradient Cosserat continuum model in capturing the meso-structural size effect and the performance of the proposed second-order computational homogenization in the simulation of strain softening and localization phenomena are demonstrated, without need to specify macroscopic phenomenological constitutive relationship and material failure model.