Search

Abstract In traditional finite element failure analyses of geotechnical structures, the micro grain rotations cannot be modelled and numerical solutions are mesh dependent. In this study, a user element including rotational degree of freedom has been developed based on micropolar theory (Cosserat theory), then an enhanced non-associated sand model is calibrated with laboratory data and used to model the plane strain tests. The simulated results demonstrate the polarized model is able to model reasonably the sand behavior as well as the grain rotations in the localized region. What’s more, with this enhanced model, the mesh independent numerical solutions in terms of mechanical responses, shear bands thickness and orientations have been obtained. Article highlights In failure analysis of geostructures, significant rotations of soil grains have been observed to occur in the strain localized regions, but the current commercial Finite Element tools cannot model the micro rotations. Therefore, a user defined element must be developed to include the rotational degree of freedom. The micropolar approach is proven to be effective to model the grain rations in present paper. More suitable than other classical soil or sand constitutive models, the selected non-associated sand model in present paper is capable of describing well the contraction and shear dilatancy behaviors of sand. Then the model has been enhanced by means of micropolar technique, in this way, the reasonable strain localization phenomena in laboratory tests could be predicted well. There are always the mesh dependent problems for traditional simulations of the strain localization phenomena in finite element analysis. It can be found in present paper that the mesh independent numerical solutions are obtained by means of micropolar technique. Furthermore, the micropolar approach can obviously improve convergence difficulties in finite element analyses.

International audience; We examine stationary regimes in granular materials from a dynamical systems theory (DST) perspective. The aim is to enrich the classical view of the critical state regime in granular materials, and more broadly, improve fundamental understanding of the underlying mesoscale mechanisms responsible for macroscopic stationary states in complex systems. This study is based on a series of discrete element method (DEM) simulations, in which 2D assemblies of non-uniformly sized circular particles are subjected to biaxial compression under constant lateral confining pressure. The lifespan and life expectancy of specific cluster conformations, comprising particles in force chains and grain loops, are tracked and quantified. Results suggest that these conformational clusters reorganize at similar rates in the critical state regime, depending on strain magnitudes and confining pressure levels. We quantified this rate of reorganization and found that the material memory rapidly fades, with an entirely new generation of force chains and grain loops replacing the old within a few percent strain.

International audience; We revisit the Critical State (CS) in a dynamic perspective at a mesoscopic scale to enrich the classical view of the critical state as a steady state. The dynamics of meso-structures has been characterized by recording generating and vanishing processes in discrete element simulations. Our results suggest that mesostructures have a short but regular lifespan at critical state along a biaxial compression loading path.

International audience; The aim of this paper is to offer a fresh perspective on the classic concept of critical state (CS) in granular materials by suggesting that CS can be defined through the use of a single proportional strain test. In classic conventional testing, CS manifests itself under constant lateral stress and controlled strain in one given direction whenever continuous shearing is applied without change being induced to material volume. However, a comparison between proportional strain tests and biaxial tests simulated with DEM has clearly shown that the CS line (CSL) characterized by stresses, void ratio and fabric indexes can act as an attractor. The mechanical responses and fabric metrics evolve along dilatant proportional strain loading paths according to similar values after the strain level has become large enough to wipe out the material memory in the homogeneous domains considered in this analysis, i.e., the shear band area in dense samples and the whole area in loose samples. This suggests that the micro-structure of a granular material subjected to any dilatant proportional strain loading paths evolves while preserving its ability to withstand shearing without volume change at any time. Therefore, the CS concept can be generalized to a wide class of loading paths which shows that CS acts as a general attractor irrespective of the loading path considered.

Due to the complexity of the hydro-mechanical behaviour of soils subjected to internal erosion, a high number of parameters are usually required for the erosion models and the constitutive models. This aspect makes it difficult to determine by trial-error their relevant values from laboratory tests. To address this issue, an efficient optimisation-based procedure for identifying the parameters of a recently developed hydro-mechanical model for internal erosion using an enhanced backtracking search algorithm (so-called MBSA-LS) has been proposed. The MBSA-LS incorporates two points: (1) modifying the mutation of the original Backtracking Search Algorithm (BSA) and (2) incorporating an efficient differential evolution (DE) as a local search to improve the optimisation performance. A mono-objective framework with six different criteria has been proposed to identify the parameters related to the interlocking effect and the erosion process. The proposed procedure was successfully applied to identify the parameters from the erosion tests of Hong Kong-Completely Decomposed Granite mixture (HK-CDG). All results demonstrated that coupling the MBSA-LS and the hydro-mechanical erosion model could efficiently solve the issue of parameter identification accounting for both the mechanical behaviour and internal erosion.

Currently, numerical studies at the real scale of an entire engineering structure considering internal erosion are still rare. This paper presents a three-dimensional (3D) numerical simulation of the effects of internal erosion within a linear dike located on a foundation. A two-dimensional (2D) finite element code has been extended to 3D in order to analyze the impact of internal erosion under more realistic hydromechanical conditions. The saturated soil has been considered as a mixture of four interacting constituents: soil skeleton, erodible fines, fluidized fine particles, and fluid. The detachment and transport of the fine particles have been modeled with a mass exchange model between the solid and the fluid phases. An elastoplastic constitutive model for sand-silt mixtures has been developed to monitor the effect of the evolution of both the porosity and the fines content induced by internal erosion upon the behavior of the soil skeleton. An unsaturated flow condition has been implemented into this coupled hydromechanical model to describe more accurately the seepage within the dike and the foundation. A stabilized finite element method was used to eliminate spurious numerical oscillations in solving the convection-dominated transport of fluidized particles. This numerical tool was then applied to a specific dike-on-foundation case subjected to internal erosion induced by a leakage located at the bottom of the foundation. Different failure modes were observed and analyzed for different boundary conditions, including the significant influence of the leakage cavity size and the elevation of the water level at the upstream and downstream sides of the dike.

In order to study the impact of internal erosion at the scale of an engineering structure, a hydro-mechanical continuous modelling approach considering suffusion is needed. It requires a relevant mechanical model for granular soils considering the-dependency (fines content) and a hydraulic model for suffusion to control the changes in the fines content. To this purpose, the mechanical models for granular soil, the unified modelling approaches for-dependency of granular soils and the hydraulic modelling of suffusion are first reviewed. Then, a hydro-mechanical model considering both suffusion and mechanical loading is developed by combining the three components. For each component of the model, alternative choices are provided. Simulations of laboratory tests as well as an example of a dike-on-foundation problem demonstrated the reliability and applicability of this coupled numerical approach for internal erosion problems.

International audience; This experimental work aimed to identify local mechanisms in remoulded clays under different triaxial loading conditions. The behaviour at a given stress level reached by two different stress paths, i.e. the conventional constant σ′3 stress path and the purely deviatoric stress path, and the behaviour at two different stress levels along a given stress path were examined. The investigation concerned both normally consolidated and overconsolidated saturated specimens of remoulded kaolin clay. The microstructural states of the clay specimens induced by different triaxial loading paths were examined at two different scales using the scanning electron microscopy and the X-ray microtomography. The results highlighted five particle orientation modes that could be activated. These modes seem to be highly dependent on the stress level, the overconsolidation ratio and the stress path. Mesocracks were found to develop within strongly overconsolidated specimens accompanied by the occurrence of dilatancy. The proposed conceptual modes provide an interesting approach to understand the mechanical behaviour of remoulded clays, especially the macroscale dilatancy phenomenon, through the possible activation of different mechanisms at micro- and mesoscales.

This paper attempts to formulate a coupled practical model in the framework of continuum mechanics to evaluate the phenomenon of internal erosion and its consequences on the mechanical behavior of soils. For this purpose, a four-constituent numerical approach has been developed to describe the internal erosion process. The detachment and transport of the fine particles have been described by a mass exchange formulation between the solid and fluid phases. The stress–strain relationship of the soil is represented by a nonlinear incremental model. Based on experimental data, this constitutive model has been enhanced by the introduction of a fines content–dependent critical state, which allows accounting for the influence of fines on soil deformation and strength. The applicability of the practical approach to capture the main features of the internal erosion process and its impact on the mechanical behavior of the eroded soil have been validated by comparing numerical and experimental results of internal erosion tests on Hong Kong completely decomposed granite (HK-CDG) mixtures, which demonstrated that the practical model was able to reproduce, with reasonable success, the experimental data. Furthermore, the influence of the stress state, the initial soil density, and the initial fraction of fines have been analyzed through numerical simulations using the proposed model.

Suffusion occurs when fines are plucked off by seepage forces and transported throughout the pores of the matrix constituted by coarser soil particles. Natural or human made soils are seldom homogeneous, which makes suffusion more complex. Suffusion is usually combined with the self-filtration of the fine particles and the transport of these fines may cause the soil structure to become looser. At the same time a clogging may occur which could reduce the permeability leading to an increase of excess pore pressure. The combination of these two phenomena will result in strength degradation. Currently, most suffusion analyses are performed without taking into account the soil's spatial variability. In this paper, a four-constituent continuum finite difference model for suffusion has been extended through the self-filtration process. A random field theory was at this point introduced into the finite difference code to investigate soil suffusion with a randomly distributed initial porosity and fines content. A probabilistic study using the Monte Carlo method was conducted to analyze the effect of the variance, the spatial correlation length, and the cross correlation of the randomly distributed initial porosity and fines content on the eroded mass and on the evolution of the hydraulic conductivity under 1D and 2D conditions. Based on all the simulations, it was possible to quantify the probability of particle blockage during erosion.

One of the major causes of instability in geotechnical structures such as dikes or earth dams is the phenomenon of suffusion including detachment, transport and filtration of fine particles by water flow. Current methods fail to capture all these aspects. This paper suggests a new modeling approach under the framework of the porous continuous medium theory. The detachment and transport of the fine particles are described by a mass exchange model between the solid and the fluid phases. The filtration is incorporated to simulate the filling of the inter-grain voids created by the migration of the fluidized fine particles with the seepage flow, and thus, the self-filtration is coupled with the erosion process. The model is solved numerically using a finite difference method restricted to one-dimensional (1-D) flows normal to the free surface. The applicability of the model to capture the main features of both erosion and filtration during the suffusion process has been validated by simulating 1-D internal erosion tests and by comparing the numerical with the experimental results. Furthermore, the influence of the coupling between erosion and filtration has been highlighted, including the development of material heterogeneity induced by the combination of erosion and filtration.

One of the major causes of instability in geotechnical structures such as dikes or earth dams is internal erosion, an insidious process that occurs over a long period of time. Research on this topic is still fairly new and much more needs to be understood in order to solve the problems posed by this phenomenon. This paper proposes a hydromechanical model based on porous continuous medium theory to assess how internal erosion impacts the safety of earthen structures. The saturated soil is considered as a mixture of four interacting constituents: soil skeleton, erodible fines, fluidized fine particles, and fluid. The detachment and transport of the fine particles are described by a mass exchange model between the solid and the fluid phases. An elastoplastic constitutive model for sand-silt mixtures has been developed to monitor the effect of the evolution of both porosity and fines content induced by internal erosion upon the behavior of the soil skeleton. The model has been numerically solved with the finite element method. It has then been applied to the specific case study of a dike foundation subjected to internal erosion induced by the presence of a karstic cavity beneath the alluvium layer. The numerical results show the onset of erosion, the time-space evolution of the eroded zone, and the hydromechanical response of the soil constituting the dike, all of which highlights the effects of the cavity location, the erosion rate, and the fines content.

International audience; The degradation of the effective normal stress in soil surrounding the pile caused by undrained cyclic loading affects the shaft capacity of the pile and can lead to structural instability. In order to investigate this phenomenon, a series of constant volume monotonic and cyclic simple shear tests on Fontainebleau sand has been conducted. Based on the experimental results, the shear stress at phase transformation state is first determined for different initial void ratios and initial normal effective stresses. Then, the number of cycles to liquefaction is estimated as a function of both the cyclic and the average shear stresses normalised by the shear stress at phase transformation. An empirical equation to evaluate the normal stress degradation is formulated and the procedure of parameter identification is presented. The performance of the suggested formulation has been analysed, based on simple shear test results on Fraser River sand and Quiou sand, and triaxial tests on Karlsruhe sand, and it has been validated by a series of additional tests on Fontainebleau sand. All the results indicate that the proposed formulation is able to estimate with good accuracy the effective normal stress degradation in sand subjected to undrained cyclic shearing.

Internal erosion occurs when fines are pulled off by seepage forces and transported throughout the matrix of soil particles, which results in the degradation of soil strength. This phenomenon has been widely studied experimentally and numerically over the last few decades. In fact, soil nature is highly variable in space and then heterogenous in magnitude [10]. At the present time, few studies have been conducted to analyze the influence of variations of soil initial state (i.e. the initial porosity and the initial fraction of fines content) on the internal erosion behaviors. In this paper, a thermodynamically consistent continuum model of internal erosion is combined with the random field theory to investigate the erosion behavior of soil with spatially variable initial state of soil in one-dimensional condition.

The tunnel leakage in sand layer plays an important role in the evolution of the lining stresses and the ground movement. The phenomenon of internal erosion occurs when underground water leaks into the tunnel from the damaged joints or the cracks of the linings. Fines in sand are pulled off by seepage forces and transported throughout the soil matrix into the tunnel. The loss of fines due to seepage flow affects the mechanical behavior of the soil around the tunnel. Conversely, the change of porosity influences the permeability of the soil and, therefore, its hydraulic conductivity. In this study, the effect of tunnel leakage in sand layer is investigated numerically with a coupled hydro-mechanical approach, considering the internal erosion phenomenon induced by the local flows of underground water, using the continuous porous medium theory and a critical-state based constitutive model. Simulation results show that the loss of fines as well as the leakage locations have significant impacts on the tunnel lining and the ground movement.

International audience; In geotechnical and particularly in dam engineering, observations in the field and, in some cases, dramatic failures of structures have often incited research and increased knowledge. For instance, the first rockfill dams were elevated without any prior knowledge of the mechanical capacities of the rockfill mass. As natural and local resources are usually preferred for economic and environmental issues, their non-optimal properties can lead to potential severe damage to the hydraulic structure. Similarly, a basic principle of construction of rockfill dams or embankments consists in sluicing rockfill after deposition and before compaction: this recommendation has been adopted for a long time. However, the beneficial effects of such a practice have been understood only in light of recent studies on unsaturated soil and multi-scale analyses.

Ferroic domain walls could play an important role in microelectronics, given their nanometric size and often distinct functional properties. Until now, devices and device concepts were mostly based on mobile domain walls in ferromagnetic and ferroelectric materials. A less explored path is to make use of polar domain walls in nonpolar ferroelastic materials. Indeed, while the polar character of ferroelastic domain walls has been demonstrated, polarization control has been elusive. Here, we report evidence for the electrostatic signature of the domain-wall polarization in nonpolar calcium titanate (CaTiO3). Macroscopic mechanical resonances excited by an ac electric field are observed as a signature of a piezoelectric response caused by polar walls. On the microscopic scale, the polarization in domain walls modifies the local surface potential of the sample. Through imaging of surface potential variations, we show that the potential at the domain wall can be controlled by electron injection. This could enable devices based on nondestructive information readout of surface potential., Comment: 30 pages, 12 figures

[Departement_IRSTEA]Eaux [TR1_IRSTEA]RIVAGE; International audience; This paper presents an experimental study of mechanical response of natural Camargue silty sand. The analysis of test results used the equivalent intergranular void ratio instead of the global void ratio. The calculation of equivalent intergranular void ratio requires the determination of parameter b which represents, physically, the fraction of active fines participating on the chain forces network, hence the strength of the soil. A new formula for determining the parameter b by using an approach based on the coordination number distribution and probability calculation is proposed. The validation of the developed relationship was done through back-analysis of published datasets in literature on the effect of fines content on silty sand behaviour. It is shown that the equivalent intergranular void ratio calculated with the b value obtained by the new formula is able to provide strong correlation to not only the critical state of but also the onset of instability of various silty sands, in different terms as peak deviator stress, peak stress ratio or cyclic resistance. Therefore, it is suggested that the use of the equivalent void ratio concept and the new b calculating formula is highly desirable in predicting of the silty sand behaviour.

International audience; Grain crushing is in part responsible for soil plasticity and dissipates most of the input energy at high stresses. This energy dissipation, as well as the energy dissipated by friction, occurs through heat exchanges. Real-time observations of temperature fields on the surface of a breaking specimen are performed by the infrared thermography technique to determine the heat generated at the crack tip when a crack is propagating. This research was motivated by recent trends in soil mechanics where constitutive models have been developed based on thermodynamic potentials including energy dissipation. A high-speed infrared camera was used to monitor the cracking of multiple specimens of rock and mortar in a diametral compression configuration. In the case of rock specimens, cracking was accompanied by a sudden rise in temperature. No temperature changes were observed in the case of mortar specimens. The maximum temperature reached and the heat dissipation profile attained depend on the nature of the rock specimen.

International audience; A novel approach has been developed to predict the mechanical behavior of very coarse granular materials with a constitutive model, which considers both grain breakage and size effect. The behavior of granular assemblies is significantly influenced by particle breakage. A critical-state double-yield-surface model incorporating the change in the critical-state line and in elastic stiffness caused by grain breakage during loading has been adopted. The amount of grain breakage was estimated by extending the size effect theory on individual grains to granular assemblies. The results from earlier studies on granular materials with parallel gradations have been usefully exploited to calibrate and to validate the model. Comparisons between experiments and simulations suggest that this approach can predict the mechanical behavior of very coarse granular materials from test results performed on a finer fraction with a homothetic gradation.

The main objective of this study is to understand the effects of the coordination number on the behavior of a single grain, before and during fracture. A new apparatus is designed, with the purpose of conducting multipoint crushing tests, and to investigate the effects of the contacts number, type and position on the failure of a particle. Some of the tests were monitored using digital image correlation (DIC), in order to get an in-depth view into the mechanics of fracture. The fragmentation of the grains was studied, and the cracks were shown to follow specific paths. The failure forces of multiple contact configurations were compared to show the effect of the number and position of contact forces on the primary and secondary cracks. It was shown that the position and magnitude of the contact forces plays a significant role in the fragmentation of the grain.

International audience; It has been widely shown that particle crushing increases the compressibility of granular materials. For a particular crushable material and given test conditions, an empirical relation can be established between the breakage ratio and the plastic work. Along these lines, constitutive models have been developed based on the effect of grading evolution during crushing. In parallel, due to corrosive attacks of the humid environment at the tip of microcracks within solid grains, the mechanical behavior of crushable granular materials depends also on the water content: the higher the material humidity, the higher the particle crushing. However, the experimental data on the relation between loading–wetting conditions and the breakage ratio are still quite scarce. In this paper, we present experimental results on crushable sand to study the effect of flooding under isotropic, oedometric, and triaxial stress paths. The main objective of this study is to obtain a consistent framework for the effect of water based on the breakage ratio. Our results have shown that, for a given initial density and stress path, the dry material after flooding reaches the equivalent behavior of the initially wetted material in terms of compression curve, particle crushing, and creep compressibility index, regardless of the point of flooding. Moreover, the relation between the breakage ratio and the final void ratio is unique and depends neither on the stress path, the water content, the point of flooding, nor the loading condition (time of creep or relaxation), but exclusively on the initial density and on intrinsic parameters. These findings could improve the prediction of the effect of water and time on the mechanical response of crushable granular materials through constitutive models based on grading evolution.; Il est reconnu que la rupture des grains augmente la compressibilité des matériaux granulaires. Pour des conditions d'essai données, une relation empirique peut être établie entre le taux de rupture des grains et le travail plastique. Ainsi, des modèles constitutifs ont été développés sur la base de l'évolution de la granulométrie liée a ` cette rupture des grains. Parallèle-ment, du fait des attaques corrosives de l'environnement humide en pointe des micro-fissures dans les grains solides, le comportement mécanique des matériaux granulaires broyables dépend également de la teneur en eau: plus l'humidité du matériau est élevée, plus la rupture semble favorisée. Toutefois, les données expérimentales sur la relation entre les conditions de chargement, les conditions hydriques et le taux de rupture sont encore très rares. Dans cet article, on présente des résultats expérimentaux sur un sable affecté par la rupture des grains pour étudier l'effet de l'imbibition sur des chemins de contraintes isotropes, eudiométriques, et triaxiales. L'objectif principal de cette étude est d'obtenir un schéma de comportement cohérent décrivant l'effet de l'eau sur le taux de rupture des grains. Les résultats montrent que, pour une densité initiale et un chemin de contraintes donné, après l'imbibition d'un échantillon sec, il suit un comportement identique à celui d'un échantillon saturé dès le début du chargement en termes de compressibilité, ruptures des grains et d'indice de compressibilité au fluage, quel que soit le point d'imbibition. Par ailleurs, la relation entre le taux de rupture et l'indice des vides final est unique et ne dépend pas du chemin de contraintes, ni de la teneur en eau, ni du point initial de l'imbibition, ni des conditions de chargement (temps de fluage ou de relaxation), mais uniquement de la densité initiale et de paramètres intrinsèques. Ces résultats pourraient améliorer la prédiction de l'effet de l'eau et du temps sur la réponse mécanique des matériaux granulaires broyables à l'aide des modèles constitutifs basés sur la granulométrie.

International audience; This paper discusses the notion of failure in a granular assembly by examining the key microstructural mechanisms which are most likely to trigger the nucleation and propagation of instabilities within a granular material. For this purpose, the key variable to predict the occurrence of failure, known as second-order work, is expressed from variables on the grain scale. The local behaviour incidents (where contacts may open or slide), compared to the global response of the assembly, are analysed by two approaches. First, numerical computations made by a discrete element model confirm the microscopic definition of the second-order work. Secondly, a micromechanical model, based on a homogenization procedure, relating the macroscopic behaviour to microscopic ingredients, namely contact planes, points to a close link between the occurrence of failure on the macroscopic scale as well as on the contact planes.

[Departement_IRSTEA]Eaux [TR1_IRSTEA]RIVAGE; International audience; This paper presents the results of an experimental investigation of the mechanical behaviour of silty sand regarding instability phenomena using conventional triaxial tests under drained and undrained conditions. First, a series of monotonic triaxial tests was carried out on reconstituted clean sand samples with emphasis placed on the effect of basic parameters such as relative density and confining pressure. To highlight the effect of fines content on the initiation of instability and on the critical state, a second series of monotonic triaxial tests was conducted on sand/silt mixture samples reconstituted with different percentages of fines. In this series, one group of samples was prepared with the same overall void ratio, while in the second one, all the specimens had the same intergranular void ratio. The effect of fines on the evolution of the maximum and the minimum void ratios was also investigated. The results obtained show that fines content has a significant influence both on resistance at the onset of instability and on the position of the critical state line of silty sands. In addition, they also revealed that the critical fine content for this soil is about 20%. We guess that, beyond this value, the effect of fines content on the mechanical response of the soil is reversed as reported by several authors in the literature.

The aim of this study is to characterise along various stress paths the microstructural evolution of a deep-water marine sediment from the Gulf of Guinea (GoG). Results obtained at the macroscopic level have shown that the phenomenological behaviour of the GoG sediment is explained by the properties of remoulded samples given by the isotropic and critical state lines, and by the properties of the intact samples given by the damage evolution of the inter-aggregate bonds. Microstructural investigations conducted by mercury intrusion porosimetry and scanning electron microscopy demonstrate that the volumetric strains measured at the macroscopic level come from inter-aggregate pore variations associated with the degrading of the inter-aggregate bonds. The mechanism responsible for reorienting the particles appears to apply in the case of the remoulded sediment. As for the natural sediment, it is first necessary to increase the isotropic consolidation stress towards higher values in order to break the bonds between the aggregates before it becomes possible to activate the particle reorientation mechanism and thus approach the remoulded behaviour.

International audience; The aim of this paper is to investigate the effect of soil grading on the stress-strain and the critical-state behavior of granular materials from idealized spheres to natural soils as well as from discrete to continuous modeling. The three-dimensional discrete-element method has been applied to study the mechanical behavior of idealized granular materials. The results confirm that the critical-state line (CSL) shifts down-ward as grading broadens with an increase of the coefficient of uniformity Cu. An exponential relationship between the critical-state parameters and the coefficient of uniformity Cu can then be established. Furthermore, experimental investigations on an artificial material (glass beads) and a natural material (Hostun sand) were carried out. The experimental results confirm the grading dependency of the critical-state behavior, and a similar relationship between CSL and Cu, for both glass beads and natural sand, extending previous literature results for high Cu values up to 20. Moreover, the experimental results show that the CSL tends to stabilize when the coefficient of uniformity Cu becomes higher than 10. Based on all test results, a simple elastoplastic model accounting for the influence of the grain size distribution (GSD) on the mechanical behavior of granular materials has been developed within the framework of critical-state soil mechanics. Comparisons between experimental results and simulations show that the model can accurately reproduce the mechanical behavior of granular materials with different GSDs.

International audience; Experimental observations have shown the significant impact of fines or coarse grains on the behavior of sand-silt mixtures. To describe the behavior of sand-silt mixtures under both drained and undrained conditions, this paper presents a mathematical model based on a micromechanical approach. The novelty of this model is the introduction of the equivalent mean size and the evolution of the position of the critical state line with fines content for various sand-silt mixtures. The predictive capability of the model was evaluated by comparing the model simulations with experimental results on undrained triaxial tests of Foundry sand-silt mixtures with fines content, f(c) = 0-100% and Ottawa sand-silt mixtures with fines content f(c) = 0-50%, and on drained triaxial tests of Hong Kong Completely Decomposed Granite (HK-CDG) mixtures before and after erosion. The predicted local behavior in the contact planes has also been examined. It shows that all local contact planes are mobilized to different degrees in terms of local stress and strain and that a few active contact planes contribute dominantly to the deformation of the assembly, leading to an anisotropic global behavior when the soil is subjected to external loading.

International audience; Testing the mechanical response of coarse granular materials requires very large and expensive laboratory equipments. During the 1960s, pioneering experimental programs were carried out on several rockfill dam materials, and those results are still a reference for engineers and researchers. However, only few experimental works have been reported to this day, and due to the scarcity of empirical data, the role of the size effect caused by grain crushing is not well known. To improve understanding of this rarely studied issue and the influence of individual particle strength, this paper analyzes the size effect on rock aggregate crushing strength and its connection with the shear envelope of rockfills. The suitability of the 4-parameter Weibull equation to describe size effects on the crushing strength reported in the literature is discussed. Furthermore, a Weibull statistical analysis was carried out for a wide number of experimental results on rock aggregates, where it has been observed that strength decreases with particle size. In parallel, the results of large triaxial tests on homothetic scaled rockfill samples of 250 and 1,000 mm in diameter reveal that the coarser the material, the higher the amount of grain breakage and the lower the shear strength. The impact of size effects obtained from the experiments is analyzed and discussed in terms of the factor of safety of rockfill slope stability. Furthermore, the results are compared with the only existing theoretical method that links the rock aggregate with the strength of the granular assembly. Good agreement between the empirical results and this theoretical method has been confirmed.

International audience; The paper summarizes the results from undrained triaxial compression tests on Merville clay, a highly overconsolidated clay at Merville in the north of France. The specimens were prepared in four different overconsolidation ratios (OCR=7, 14, 28, 56), and for each OCR, undrained shear tests were performed using three or four axial strain rates varying from 0.26 %/h to 25.89 %/h. The strain-rate dependency of the undrained stress level q/q f , the pore pressure, the effective stress path and the stress ratio q/p' were investigated. The undrained shear strength was found to increase with the strain rate. The shear band formation during loading was observed for all strain rates. For each OCR test series, a linear relationship between the normalized peak deviatoric stress and the strain rate was obtained in the semi-logarithm plane. As OCR increasing, the slope of best-fitting line for the linear relationship increases. Therefore, the strain rate influence on the normalized undrained shear strength is greater for higher OCRs.

International audience; This paper aims at investigating the influence of grading on the undrained behavior of granular materials. Series of undrained triaxial tests were carried out with two different materials, glass balls and Hostun sand. For each material, samples with different gradings and similar relative densities were prepared. Experimental results show that the undrained shear strength decreases when the coefficient of uniformity Cu=d60/d10Cu=d60/d10 increases from 1.1 to 20. The conditions of instability for the selected granular materials were also analyzed, based on the sign of the second-order work during undrained triaxial loading. The results demonstrate a significant influence of the grading: increasing the coefficient of uniformity heightens the potential of static liquefaction and the materials become more unstable. Furthermore, numerical tests using the three-dimensional discrete element method (DEM) were conducted on assemblies of spheres. The DEM inter-particle parameters were derived from the experimental test results on glass balls. The DEM simulations showed similar behaviors compared to experimental results and confirmed the influence of the grain size distribution on the stress–strain relationship and instability phenomena.; Cet article vise à étudier l'influence de l'étalement granulométrique sur le comportement non drainé des matériaux granulaires. Des séries d'essais triaxiaux non drainés ont été réalisés sur deux matériaux différents, des billes de verre et du sable Hostun. Pour chaque matériau, des échantillons de différentes granulométries et de densités relatives similaires ont été préparés. Les résultats expérimentaux montrent que la résistance au cisaillement diminue lorsque le coefficient d'uniformité Cu=d60/d10Cu=d60/d10 augmente de 1,1 à 20. Les conditions d'instabilité pour les matériaux granulaires sélectionnés ont également été analysées, sur la base du signe du travail du second ordre au cours d'un chargement triaxial non drainé. Les résultats montrent une influence significative de l'étalement granulométrique : l'augmentation du coefficient d'uniformité accroît le potentiel de liquéfaction statique et les matériaux deviennent plus instables. En outre, des essais numériques en utilisant la méthode des éléments discrets en trois dimensions (DEM) ont été effectués sur des assemblages de sphères. Les paramètres inter-particules DEM ont été calculés à partir des résultats d'essais expérimentaux sur des billes de verre. Les simulations DEM ont montré des comportements similaires comparés aux résultats expérimentaux et confirment l'influence de l'étalement granulométrique sur la relation contrainte–déformation et les phénomènes d'instabilité.

International audience; Depending on stress state and loading conditions, instability in loose sand can develop at a shear stress level much lower than the Mohr-Coulomb failure criterion. This type of failure mode, very different from strain-softening in dense sand, can lead to catastrophic collapse of earth structures. An appropriate constitutive model is therefore needed that can capture the correct failure modes and accurately predict the stress-strain response for both dense and loose sands. In this paper, a micromechanical approach developed earlier for the modelling of granular material behaviour is adopted. The constitutive model is utilised to simulate undrained triaxial, constant-q and proportional strain tests on loose sand, and numerical results are compared with experimental results. The model is capable of accurately reproducing the macroscopic experimental tests, and is also able to provide microstructural information. This enables the observation of instability at the particle scale, and provides a better understanding of the linkage between particle-scale instability and instability at the assembly level.

International audience; Since the pioneering work of Griffith, Linear Elastic Fracture Mechanics has been widely experimentally validated and successfully developed in solid mechanics modeling. However, recent theoretical models applying the energy balance found in Griffith theory specifically for quasi-static confined comminution have until now not been systematically confronted to experiments. In this study, we analyze data of compression tests on crushable sand, where grain breakage has been triggered by flooding the initially dry material at constant stresses. We consider a partition of the dissipation between surface fracture energy and the rearrangement of fragments and grains surrounding crushed particles. Our results show that the role of the surface fracture energy is stressdependent and that its influence becomes less significant at high stresses.

The parameters obtained from oedometric consolidation tests are commonly used in the development of constitutive modeling and for engineering practice. This paper focuses on the influence of the natural deposition plane orientation on oedometric consolidation behavior of three natural clays from the southeast coast of China. Oedometer tests were conducted on intact specimens prepared by sampling at a series of angles relative to the natural deposition plane. For each specimen, yield stress, compressibility indexes, secondary compression, and permeability coefficients were determined. The influence of the sampling angle on these properties was investigated, revealing that yield stress, compression index, swelling index, creep index, ratio of secondary compression coefficient to compression index (C αe/C c) and permeability coefficient were all dependent to some extent on the sampling angle. These findings indicate the role of the anisotropy due to the natural deposition on the oedometric consolidation behavior.

International audience; The critical state concept has been widely used in soil mechanics. The purpose of this study is to apply this concept in the framework of multi-mechanism elastoplasticity. The developed model has two yield surfaces: one for shear sliding and one for compression. In this model, the location of the critical state line is explicitly considered and related to the actual material density to control the peak strength and the phase transformation characteristics. The stress reversal technique is incorporated into the model for describing clay behavior under complex loading including changes of stress direction. The determination of the model parameters is discussed; it requires only one drained or undrained triaxial test up to failure with an initial isotropic consolidation stage. The model is used to simulate drained and undrained tests under monotonic loading with different over-consolidation ratios on various remolded and natural clays, including true triaxial tests with different Lode's angles. Drained and undrained tests under cyclic loadings are also simulated by using the set of parameters determined from monotonic tests. The comparison between experimental results and numerical simulations demonstrate a good predictive ability of this new simple model.

International audience; Experimental results show the significant influence of cement content on the mechanical properties of cement-treated clays. Cementation is produced by mixing a certain amount of cement with the saturated clay. The purpose of this paper is to model the cementation effect on the mechanical behavior of cement-treated clay. A micromechanical stress-strain model is developed considering explicitly the cementation at inter-cluster contacts. The inter-cluster bonding and debonding during mechanical loading are introduced in two ways: an additional cohesion in the shear sliding and a higher yield stress in normal compression. The model is used to simulate isotropic compression and undrained triaxial tests under various confining stresses on cement-treated Ariake clay and Singapore clay with various cement contents. The applicability of the present model is evaluated through comparisons between numerical and experimental results. The evolution of local stresses and local strains in inter-cluster planes are discussed in order to explain the induced anisotropy due to debonding at contact level under the applied loads. The numerical simulations demonstrate that the proposed micromechanical approach is well adapted for taking into account the main physical properties of cement-treated clay, including damage and induced anisotropy under mechanical loading.

This paper aims to investigate the influence of grain size distribution (GSD) on the critical state of Hostun sand. Drained triaxial tests were conducted on saturated specimens with different GSD (Cu=1.4, 5, 10, 20) under different confining pressures (100, 200, 400 kPa). The experimental results highlight that the critical state of the tested sand predominantly depends on its GSD.

International audience; The shear strength of coarse granular materials, widely used in civil works such as rockfill in dam construction, is seldom measured because of severe practical experimental limitations. This paper presents an original method for evaluating the shear strength of such materials, based on size effects in granular materials affected by grain breakage according to fracture mechanics. Through a general size effect relation operating on shear strength envelopes, this method makes it possible to determine the shear strength of a coarse-grained granular material from the measured properties of a finer-grained granular material made of the same mineral. In the paper, the method is explicitly proved, taking into account the statistical distribution of breakage resistances of particles with different sizes, within materials that can be considered as physically similar, particularly having parallel grain size distributions and the same compactness. A wide set of independent experimental results is shown to validate the method consistently.

International audience; A particle breakage has a significant impact on the mechanical behavior of granular materials. In this paper, we present an elasto-plastic model with two yield surfaces to which the influence of particle breakage has been introduced. The main feature of this model is to incorporate the change in the critical state line (CSL) consequent to the grain breakage induced by isotropic and deviatoric stresses during loading. For this purpose we propose a breakage function which connects the evolution of the CSL to the energy consumed. Results from earlier studies on drained and undrained compression and extension triaxial tests were used to calibrate and validate the model. Comparison between earlier results and our simulations indicates that the model can reproduce with good accuracy the mechanical behavior of crushable granular materials and predict the evolution of the grain size distribution during loading.

International audience; Natural soils or coarse granular materials containing large particles such as rockfill or ballast are often an issue for geotechnical projects. However, since it is usually complicated or impossible to set up a big enough experimental apparatus to test these materials, their mechanical behavior remains poorly known.To circumvent this difficulty, we have developed a set of triaxial cells of various sizes to study the experimental behavior of granular soils containing particle sizes of up to 160 mm. This paper presents the first results of a two-part experimental study. The first part consists of triaxial tests on sand samples of different sizes in order to examine specimen scale effects. Pre-peak behavior is not affected by the specimen size, whereas post-peak behavior depends on the test conditions that control the development of strain localizations. In the second part, we present drained triaxial compression tests on large samples of calcareous rockfill materials that have parallel grain size distributions.The results show a nonstandard evolution of the friction angle in relationship to grain size if compared to similar studies on coarse granular materials. Although grain crushing occurs during compression, the strain-stress curves of the different materials with different grain sizes do remain close. We explain this phenomenon by the strength evolution of the individual grains with respect to their sizes

The compressibility of engineered barrier clays is, to a large extent, controlled by microstructure change due to the presence of chemical ions in clay-water system. This paper aims to investigate the change of microstructure of clays due to the presence of heavy metal ions in pore water. We use two pure clays (kaolinite and bentonite) in the study. One-dimensional consolidation tests were performed on reconstituted samples, which are prepared with distilled water and three types of heavy metal solutions (Pb(NO3)2, Cu(NO3)2, Zn(NO3)2,). In order to better understand the impact of chemical pore fluid on microstructure of the two clays, following the consolidation test, scanning electron microscope (SEM) observations and mercury intrusion pore size distribution measurements (MIP) were conducted. Due to the measurement range of MIP, which is only allowed to measure the minimal pore size 20 Å, BET method by gas sorption, whose measurement pore size range is from 3.5 Å to 500 Å, is used to measure the micropore size distribution. By this method, specific surface area of the soils can be also determined. It can be employed to demonstrate the difference of creep performance between the soils. Furthermore, a series of batch equilibrium tests were conducted to better understand the physical-chemical interactions between the particles of soils and the heavy metal ions. With the further consideration of the interparticle electrical attractive and repulsive force, an attempt has been made to predict the creep behaviour by using the modified Gouy-Chapman double layer theory. The results of calculation were compared with that of tests. The comparison shows that the prediction of compressibility of the clays according to the modified double diffuse layer theory can be reasonably agreement with the experimental data.

International audience; An elasto-plastic model with two yield surfaces has been developed for the simulation of granular materials with consideration of grain ruptures. Grain breakage is induced by deviatoric as well as isotropic stresses. One yield function is based on a Mohr-Coulomb criterion with a hyperbolic hardening function of the plastic deviatoric strain. A second yield function is introduced to describe the plastic behaviour under compression, in which a hardening function of the plastic volumetric strain is introduced. The main assumption of this model is that, upon loading, the position of the critical state changes as a consequence of grain breakage. The effect of the grain size distribution is introduced in the relationship between the void ratio at critical state and the mean effective stress. Triaxial and oedometer tests have been performed on crushable granular materials. Comparison of experimental results and numerical simulations shows that the new model can reproduce with good accuracy the behaviour of granular materials subjected to grain ruptures during mechanical loading.

Because of the lack of experimental devices for testing granular materials with large particles, the mechanical behavior of rockfill is not very well known. This lack of knowledge can lead to improper design of rockfill dams, and as a consequence, to severe accidents including dam failures. A method has been developed by Frossard (2005) for evaluating the shear strength of rockfill based on the estimation of size effect in granular materials, affected by grain breakage. Comparison between experimental results and numerical predictions shows that this method can estimate with good accuracy the strength of granular materials made of large particles.

Granular materials with large particles have been widely used to construct embankments or rockfills dams. However, due to limitations in the size of experimental set-ups, their behavior remains relatively less investigated than the behavior of sands for instance. A new large triaxial apparatus, with a maximum specimen diameter of 1000mm, has been designed to study the behavior of such coarse granular materials. This device is described in this paper. The results of a triaxial test on calcareous particles with a maximum particle size of 160 mm are presented. Special attention has been paid on grain crushing under shear loading.

International audience; Under certain loading conditions, instability in loose sand can develop at a shear stress level much lower than the critical state failure line. To analyse these types of problems, we adopt the micromechanical approach developed earlier for the modelling of granular material behaviour. The stress-strain relationship for a granular assembly is determined by integrating the behaviour of the inter-particle contacts in all orientations. The constitutive model is applied to simulate undrained triaxial and constant-q tests on loose sand. Numerical simulations are compared to experimental results in order to evaluate the model's ability to predict the different modes of instability of the granular assembly. With the model we have also analysed instability at inter-particle level to observe how it is connected to instability at the assembly level.; Dans certaines conditions de chargement, des ruptures par instabilité peuvent prendre naissance dans les sables lâches pour des niveaux de contraintes très nettement inférieurs à ceux correspondant à la plasticité parfaite. Pour analyser ce type de problèmes, nous adoptons une approche d'homogénéisation développée précédemment pour la modélisation du comportement des matériaux granulaires. La relation contrainte-déformation macroscopique est obtenue en intégrant le comportement défini au niveau du contact inter-granulaire. Le modèle est utilisé pour simuler des essais non drainés et des essais à déviateur constant sur un sable lâche. Les résultats de la simulation à l'échelle macroscopique sont comparés à des résultats expérimentaux de façon à apprécier la capacité du modèle à prédire la naissance de l'instabilité sur ces chemins particuliers. Une analyse à l'échelle des différents plans de contact inter-granulaire a également été entreprise afin de relier les conditions d'instabilité aux différentes échelles.

International audience; In light of environmental differences (such as gravitational fields, surface temperatures, atmospheric pressures, etc.), the mechanical behavior of the subsurface soil on the Moon is expected to be different from that on the Earth. Before any construction on the Moon can be envisaged, a proper understanding of soil properties and its mechanical behavior in these different environmental conditions is essential. This paper investigates the possible effect of surface-energy forces on the shear strength of lunar soil. All materials, with or without a net surface charge, exhibit surface-energy forces, which act at a very short range. Although, these forces are negligible for usual sand or silty sand on Earth, they may be important for surface activated particles under extremely low lunar atmospheric pressure. This paper describes a constitutive modeling method for granular material considering particle level interactions. Comparisons of numerical simulations and experimental results on Hostun sand show that the model can accurately reproduce the overall mechanical behavior of soils under terrestrial conditions. The model is then extended to include surface-energy forces between particles in order to describe the possible behavior of lunar soil under extremely low atmospheric pressure conditions. Under these conditions, the model shows that soil has an increase of shear strength due to the effect of surface-energy forces. The magnitude of increased shear strength is in reasonable agreement with the observations of lunar soil made on the Moon's surface.

Grouting is generally used for ground reinforcement or for reducing soil permeability. However, choosing grouting parameters or predicting improvements is at present rather empirical. A numerical approach may therefore enable grouting to be more efficient. The purpose of this paper is to model the evolution, in time and space, of two major variables: grout concentration and fluid pressure in grouted soil. After a brief description of a first tried diffusive model, this study aims to describe the performance of a software program, Athos. This program was initially developed for the oil industry and uses a piston model. Results are given in one- and three-dimensional approaches. They underline a rapid evolution of concentration from grout- saturated zones to non-saturated zones. The evolution of pressure depending on this degree of saturation is also highlighted. The second part of the paper compares the numerical model with a series of experiments carried out on grouted sand columns and in saturated sand tanks. In the latter tests, the grout is injected into the sand using a sleeved grout pipe, and the progress of grout propagation in the soil is monitored by acoustic emission (AE). This comparison shows the validity of the chosen model. The permeation distance detected by AE is compared with that predicted by simple grout propagation models.