The article presents a multi-scale modelling approach of cohesive granular materials, its numerical implementation and its results. At microscopic level, Discrete Element Method (DEM) is used to model dense grains packing. At the macroscopic level, the numerical solution is obtained by a Finite Element Method (FEM). In order to bridge the micro- and macro-scales, the concept of Representative Elementary Volume (REV) is applied, in which the average REV stress and the consistent tangent operators are obtained in each macroscopic integration point as the results of DEM’s simulation. In this way, the numerical constitutive law is determined through the detailed modelling of the microstructure, taking into account the nature of granular materials. We first elaborate the principle of the computation homogenisation (FEM × DEM), then demonstrate the features of our multiscale computation in terms of a biaxial compression test. Macroscopic strain location is observed and discussed. 相似文献
The paper describes the basic ideas and the main features of a new class of constitutive laws, in the framework of incrementally non-linear constitutive equations. CLoE is a generic name for that new class of laws, with reference to consistency at the limit surface, and explicit localization analysis. A top–down analysis of the model is presented, and illustrated by examples. 相似文献
Double-scale numerical methods constitute an effective tool for simultaneously representing the complex nature of geomaterials and treating real-scale engineering problems such as a tunnel excavation or a pressuremetre at a reasonable numerical cost. This paper presents an approach coupling discrete elements (DEM) at the microscale with finite elements (FEM) at the macroscale. In this approach, a DEM-based numerical constitutive law is embedded into a standard FEM formulation. In this regard, an exhaustive discussion is presented on how a 2D/3D granular assembly can be used to generate, step by step along the overall computation process, a consistent Numerically Homogenised Law. The paper also focuses on some recent developments including a comprehensive discussion of the efficiency of Newton-like operators, the introduction of a regularisation technique at the macroscale by means of a second gradient framework, and the development of parallelisation techniques to alleviate the computational cost of the proposed approach. Some real-scale problems taking into account the material spatial variability are illustrated, proving the numerical efficiency of the proposed approach and the benefit of a particle-based strategy. 相似文献
This paper reports some results of a large experimental program on Boom Clay conducted in Grenoble in the framework of the European project SELFRAC. The program included isotropic compression up to relatively high stress, drained triaxial compression tests at different cell pressures, as well as permeability measurements under isotropic and deviatoric stress. Local measurement of axial and radial displacements allowed the detection of strain localization during deviatoric loading. The permeability of Boom Clay is found to depend on the mean effective stress. The response of Boom Clay during deviatoric loading appears to be strongly affected by the swelling experienced during the isotropic stage preceding triaxial compression. The rate of swelling decreases with isotropic stress. The longer the swelling before shear, more the response under shear becomes ductile and the lower the initial stiffness. Permeability depends on the mean effective stress and it is found to decrease of about two orders of magnitude when the mean stress increases from 1 to 32 MPa. Permeability during shear loading is essentially constant and does not seem to be affected by strain localization. These results are complemented by a few observations obtained using X-ray microtomography in the framework of the more recent European project TIMODAZ. These findings illustrate the impact of pre-existing inclusions and fissures on specimen deformation upon deviatoric loading in the laboratory. 相似文献
Shear bands with characteristic spatial patterns observed in an experiment for a cubic or parallelepiped specimen of dry dense sand were simulated by numerical bifurcation analysis using the Cam‐clay plasticity model. By incorporating the subloading surface concept into the plasticity model, the model became capable of reproducing hardening/softening and contractive/dilative behavior observed in the experiment. The model was reformulated to be compatible with the multiplicative hyperelasto‐plasticity for finite strains. This enhanced constitutive model was implemented into a finite‐element code reinforced by a stress updating algorithm based on the return‐mapping scheme, and by an efficient numerical procedure to compute critical eigenvectors of elastoplastic tangent stiffness matrix at bifurcation points. The emergence of diamond‐ and column‐like diffuse bifurcation modes breaking uniformity of the materials, followed by the evolution of shear bands through strain localization, was observed in the analysis. In the bifurcation analysis of plane strain compression test, unexpected bifurcation modes, which broke out‐of‐plane uniformity and led to 3‐dimensional diamond‐like patterns, were detected. Diffuse bifurcations, which were difficult to observe by experiments, have thus been found as a catalyst creating diverse shear band patterns. 相似文献
This paper presents an analysis of two large rock toppling/sliding events which occurred in January 2014 and February 2019 at the Cliets unstable slope (Savoie, French Alps). To understand the mechanism involved and its control by external forcings, a multi-technique analysis approach is used combining geological observations, meteorological data analysis, topographic measurements and simple physical modeling. The pre-failure stage of the events is more particularly analyzed. No direct relationships are found between triggering factors and surface motion though a kinematics analysis highlights the transition toppling-sliding. It showed that, at first order, this transition occurred 4 years before the first failure of 2014, while it happened 2 months before the second failure of 2019. From this date, the environment is considered like a block sliding on an inclined plane. By applying a frictional model (Helmstetter et al. in Journal of Geophysical Research: Solid Earth 109(B2), 2004), we illustrated that the two events belong to an unstable velocity-weakening sliding regime. The time to failure (Voight in Science 243(4888):200–203, 1989) is forecasted with the model, and the results are consistent with the observations. They confirm that the gravitational factor is predominant over the triggering factors for the two events.