A theoretical–experimental approach to elastic and strength properties of artificially cemented sand

Evaluating the behavior parameters of soils and soil-binder mixes by means of theoretical models that are supported by laboratory tests still remains a key challenge in foundation design. In this context, the paper investigates some aspects of the mechanical behavior of artificially cemented sands (ACS) by means of experimental characterization and micromechanics-based modeling. Particular emphasis is given to the increase in elastic stiffness and strength brought by cementation. Based on the concept of a fictitious continuum medium and the homogenization theory, the effective elastic properties of ACS are evaluated using the Mori–Tanaka and self-consistent schemes. The elastic micromechanical approach is supported by bender element tests. Finally, the effective strength properties of ACS are assessed by means of micromechanics-based failure criterion formulated within the context of non-associated plasticity. Validation and calibration of the theoretical model are achieved by comparison with data from unconfined compression tests.     

Characterization of the mechanical properties of a claystone by nano-indentation and homogenization

Based on the analyses of mineralogical compositions by X-ray diffraction and microstructure by optical microscopy, the Young’ modulus and hardness of a claystone were characterized by the nano-indentation technique and homogenization method. Three distinct microstructural zones are identified in the claystone: clay matrix, a composite matrix of clay and small mineral grains and imbedded quartz grains. The elastic modulus and hardness of different zones were determined by nano-indentation testing. Based on the statistical analysis of nano-indentation results, the spatial mappings and frequency distributions of elastic modulus and hardness of the different zones were obtained. The elastic moduli of main constituent phases of the claystone are then estimated from the nano-indentation tests. These values were further used for the determination of the macroscopic elastic modulus of the claystone using two different homogenization schemes: the dilute scheme and Mori–Tanaka scheme. The predicted values by the homogenization schemes are compared with experimental data obtained from conventional uniaxial compression tests.     

A comparative micromechanical analysis of the effective properties of a geomaterial: Effect of mineralogical compositions

In order to provide a physical interpretation of the variation of the mechanical properties of Callovo-Oxfordian argillite with mineral composition, we implement three linear homogenization schemes. The argillite is modeled as a three phase material composed of a clay matrix and inclusions of quartz and calcite. It is shown that, unlike the dilute scheme and the self-consistent scheme, the Mori-Tanaka model describes the in situ experimental data well. The determined properties are finally used in a finite element computation. The aim is to evaluate the effect of mineral composition on the elastic response of the excavation of a vertical shaft in the context of the underground laboratory of Meuse/Haute Marne.     

裂隙岩体三维柔度张量的数值试验方法

在裂隙岩体二维柔度张量数值试验的基础上,建立了裂隙岩体三维柔度张量及其表征单元体积（REV）尺度研究的简化数值试验方法。首先根据岩体裂隙的统计参数及分布规律,运用蒙特卡罗方法在研究域内获得与实际岩体裂隙同分布的三维随机裂隙网络,提取位于不同方位的岩体试件,运用二维柔度张量的数值试验方法求出各个平面方向上的二维柔度张量,然后根据二维与三维柔度张量的拓扑关系,用数学方法求解裂隙岩体的三维柔度张量。对于含3组正交全贯通裂隙的岩体,通过数值试验获得的柔度张量与理论解相比,其主对角线上各参数的误差在5%以内,表明该数值试验方法较为可靠。最后将此方法运用到小湾水电站工程中,确定左岸坝区裂隙岩体的应力REV为11 m×11 m×11 m,并获得该区域裂隙岩体的三维柔度张量。     

Large‐scale poroelastic fractured reservoirs modeling using the fast multipole displacement discontinuity method

An effective approach to modeling the geomechanical behavior of the network and its permeability variation is to use a poroelastic displacement discontinuity method (DDM). However, the approach becomes rather computationally intensive for an extensive system of cracks, particularly when considering coupled diffusion/deformation processes. This is because of additional unknowns and the need for time‐marching schemes for the numerical integration. The Fast Multipole Method (FMM) is a technique that can accelerate the solution of large fracture problems with linear complexity with the number of unknowns both in memory and CPU time. Previous works combining DDM and FMM for large‐scale problems have accounted only for elastic rocks, neglecting the fluid leak‐off from the fractures into the matrix and its influence on pore pressure and stress field. In this work we develop an efficient geomechanical model for large‐scale natural fracture networks in poroelastic reservoirs with fracture flow in response to injection and production operations. Accuracy and computational performance of the proposed method with those of conventional poroelastic DDM are compared through several case studies involving up to several tens of thousands of boundary elements. The results show the effectiveness of the FMM approach to successfully evaluate field‐scale problems for the design of exploitation strategies in unconventional geothermal and petroleum reservoirs. An example considering faults reveals the impact of reservoir compartmentalization because of sealing faults for both geomechanical and flow variables under elastic and poroelastic rocks. Copyright © 2015 John Wiley & Sons, Ltd.     

A micro-mechanics-based elastic–plastic model for porous rocks: applications to sandstone and chalk

A micro-mechanics-based elastic–plastic model is proposed to describe mechanical behaviors of porous rock-like materials. The porous rock is considered as a composite material composed of a solid matrix and spherical pores. The effective elastic properties are determined from the classical Mori–Tanaka linear homogenization scheme. The solid matrix verifies a pressure-dependent Mises–Schleicher-type yield criterion. Based on the analytical macroscopic yield criterion previously determined with a nonlinear homogenization procedure (Shen et al. in Eur J Mech A/Solids 49:531–538, 2015), a complete elastic–plastic model is formulated with the determination of a specific plastic hardening law and plastic potential. The micro-mechanics-based elastic–plastic model is then implemented for a material point in view of simulations of homogeneous laboratory tests. The proposed model is applied to describe mechanical behaviors of two representative porous rocks, sandstone and chalk. Comparisons between numerical results and experimental data are presented for triaxial compression tests with different confining pressures, and they show that the micro-mechanical model is able to capture main features of mechanical behaviors of porous rock-like rocks.     

From micron-sized needle-shaped hydrates to meter-sized shotcrete tunnel shells: micromechanical upscaling of stiffness and strength of hydrating shotcrete

Knowledge on the stresses in shotcrete tunnel shells is of great importance, as to assess their safety against severe cracking or failure. Estimation of these stresses from 3D optical displacement measurements requires shotcrete material models, which may preferentially consider variations in the water–cement and aggregate–cement ratios. Therefore, we employ two representative volume elements within a continuum micromechanics framework: the first one relates to cement paste (with a spherical material phase representing cement clinker grains, needle-shaped hydrate phases with isotropically distributed spatial orientations, a spherical water phase, and a spherical air phase; all being in mutual contact), and the second one relates to shotcrete (with phases representing cement paste and aggregates, whereby aggregate inclusions are embedded into a matrix made up by cement paste). Elasticity homogenization follows self-consistent schemes (at the cement paste level) and Mori–Tanaka estimates (at the shotcrete level), and stress peaks in the hydrates related to quasi-brittle material failure are estimated by second-order phase averages derived from the RVE-related elastic energy. The latter permits upscaling from the hydrate strength to the shotcrete strength. Experimental data from resonant frequency tests, ultrasonics tests, adiabatic tests, uniaxial compression tests, and nanoindentation tests suggest that shotcrete elasticity and strength can be reasonably predicted from mixture- and hydration-independent elastic properties of aggregates, clinker, hydrates, water, and air, and from strength properties of hydrates. At the structural level, the micromechanics model, when combined with 3D displacement measurements, predicts that a decrease of the water–cement ratio increases the safety of the shotcrete tunnel shell.     

Effective properties of a cemented or an injected granular material

In this paper the macroscopic elastic properties of injected or cemented sands are derived from the characteristics of the constituents and the analysis of the microstructure using a multi‐scale modelling approach. Particular interest is given to the choice of the representative elementary volume, by relying on existing microstructural data. The periodic homogenization is adopted and required numerical solutions are performed by the finite element method. An assessment of the validity of the multi‐scale approach is achieved through comparison with theoretical and experimental results on cemented and injected granular media reported in the literature. The capabilities of the model are also used to investigate the influence of geometrical and mechanical microscale parameters on the macroscopic behaviour of the treated materials. Copyright © 2004 John Wiley & Sons, Ltd.     

A multiscale homogenization model for strength predictions of fully and partially frozen soils

Soil freezing is often used to provide temporary support of soft soils in geotechnical interventions. During the freezing process, the strength properties of the soil–water–ice mixture change from the original properties of the water-saturated soil to the properties of fully frozen soils. In the paper, a multiscale homogenization model for the upscaling of the macroscopic strength of freezing soil based upon information on three individual material phases—the solid particle phase (S), the crystal ice phase (C) and the liquid water phase (L)—is proposed. The homogenization procedure for the partially frozen soil–water–ice composite is based upon an extension of the linear comparison composite (LCC) method for a two-phase matrix–inclusion composite, using a two-step homogenization procedure. In each step, the LCC methodology is implemented by estimating the strength criterion of a two-phase nonlinear matrix–inclusion composite in terms of an optimally chosen linear elastic comparison composite with a similar underlying microstructure. The solid particle phase (S) and the crystal ice phase (C) are assumed to be characterized by two different Drucker–Prager strength criteria, and the liquid water phase (L) is assumed to have zero strength capacity under drained conditions. For the validation of the proposed upscaling strategy, the predicted strength properties for fully and partially frozen fine sands are compared with experimental results, focussing on the investigation of the influence of the porosity and the degree of ice saturation on the predicted failure envelope.     

Analysis of stone-column reinforced foundations

A numerical model is proposed to analyse elastic as well as elastoplastic behaviour of stone-column reinforced foundations. The stone-columns are assumed to be dispersed within the in situ soil and a homogenization technique is invoked to establish equivalent material properties for in situ soil and stone-column composite. The difficulties encountered in carrying out elastoplastic analyses of composite materials are overcome by adopting a separate yield function for each of the constituent materials and a sub-iteration procedure within an implicit backward Euler stress integration scheme. In the proposed procedure, equilibrium as well as kinematic conditions implied in the homogenization procedure are satisfied for both elastic as well as elastoplastic stress states. The proposed model is implemented in an axi-symmetric finite element code and numerical prediction is made for the behaviour of model circular footings resting on stone-column reinforced foundations. This prediction indicates good agreement with experimental observation. Finally, a new scheme in which the length of stone-column is variable is proposed and its behaviour is examined through a numerical example. © 1998 John Wiley & Sons, Ltd.     

Elastic wave velocities and permeability of cracked rocks

Cracks play a major role in most rocks submitted to crustal conditions. Mechanically, cracks make the rock much more compliant. They also make it much easier for fluid to flow through any rock body. Relying on Fracture Mechanics and Statistical Physics, we introduce a few key concepts, which allow to understand and quantify how cracks do modify both the elastic and transport properties of rocks. The main different schemes, which can be used to derive the elastic effective moduli of a rock, are presented. It is shown from experimental results that an excellent approximation is the so-called non-interactive scheme. The main consequences of the existence of cracks on the elastic waves is the development of elastic anisotropy (due to the anisotropic distribution of crack orientations) and the dispersion effect (due to microscopic local fluid flow). At a larger scale, macroscopic fluid flow takes place through the crack network above the percolation threshold. Two macroscopic fluid flow regimes can be distinguished: the percolative regime close to the percolation threshold and the connected regime well above it. Experimental data on very different rock types show both of these behaviors.     

Micromechanical approach to effective viscoelastic properties of micro‐fractured geomaterials

The aim of this paper is to formulate a micromechanics‐based approach to non‐aging viscoelastic behavior of materials with randomly distributed micro‐fractures. Unlike cracks, fractures are discontinuities that are able to transfer stresses and can therefore be regarded from a mechanical viewpoint as interfaces endowed with a specific behavior under normal and shear loading. Making use of the elastic‐viscoelastic correspondence principle together with a Mori‐Tanka homogenization scheme, the effective viscoelastic behavior is assessed from properties of the material constituents and damage parameters related to density and size of fractures. It is notably shown that the homogenized behavior thus formulated can be described in most cases by means of a generalized Maxwell rheological model. For practical implementation in structural analyses, an approximate model for the isotropic homogenized fractured medium is formulated within the class of Burger models. Although the approximation is basically developed for short‐term and long‐term behaviors, numerical applications indicate that the approximate Burger model accurately reproduce the homogenized viscoelastic behavior also in the transient conditions.     

Compressibility and permeability of sand–kaolin mixtures. Experiments versus non‐linear homogenization schemes

This study deals with the behaviour of mixtures of sand and saturated kaolin paste considered as composite materials made of permeable and deformable (with non‐linear behaviour) matrix (the kaolin paste) with rigid and impervious inclusions (the sand grains). Oedometric and permeability tests highlight the key role of the state of the clay paste, and show the existence of a threshold of sand grain concentration above which a structuring effect influences both compressibility and permeability. At the light of these experiments two homogenization schemes (with simplifying assumptions to make the problem manageable) are considered to model these two parameters. Qualitative and quantitative comparisons with experimental data point out their respective domain of interest and limitations: a tangent homogenization scheme is shown to be sufficient to describe the macroscopic properties for dilute sand concentration; above the concentration threshold, the structuring effect is captured by the new homogenization scheme developed in this paper. Copyright © 2010 John Wiley & Sons, Ltd.     

劈裂注浆复合土体三维等效弹性模型理论研究

劈裂注浆可以有效改善土体的变形参数,大大降低土体在受力状态改变时的变形量,对劈裂注浆后复合土体的等效变形参数进行研究十分重要。在综合分析劈裂注浆扩散机制和工程应用实际的基础上,基于均质化理论提出了劈裂注浆后复合土体的三维单元体几何模型,按等效原则给出了浆-土体积及受力分配关系模型图;接着基于横向各向同性本构关系推导了模型的等效弹性模量和等效泊松比的解析解。然后采用有限元方法取得了模型特定条件下的等效弹性模量和等效泊松比,并与解析结果进行对比分析。最后把模型和相应的解析结果引入Flac3D岩土工程专业分析软件,结合一个热力隧道工程实例对隧道劈裂注浆后关键位置的沉降进行预测分析,并与实测值进行了对比。研究表明：对所提出的计算模型,解析计算与有限元方法计算结果吻合度较高,说明了解析结果的正确性;基于该模型及其解析结果得到的隧道开挖后的沉降预测值与实测值具有良好的一致性,说明所提出的模型和相应的解析计算方法具有一定的可靠性和实用性。     

Equivalent properties of rock strata: static and dynamic analysis

The main purpose of this investigation is to study the state of stress of layered rocks forming slopes of deep river valleys. For this purpose averaging technique and a variant of the Variation-Difference Method are used. Because of the averaging method, equivalent homogeneous properties of layered elastic medium are determined. The paper has two parts. The first one is devoted to the analysis of a static stress–strain state of the slopes under gravity. The rock mass in the second part is subjected to dynamic loading caused by an earthquake. As a result of the numerical solution of the raised problems, the stress distribution in slopes and at the base of a deep canyon-like river valley was obtained. © 1997 John Wiley & Sons, Ltd.     

Homogenization framework for three‐dimensional elastoplastic finite element analysis of a grouted pipe‐roofing reinforcement method for tunnelling

This paper deals with the grouted pipe‐roofing reinforcement method that is used in the construction of tunnels through weak grounds. This system consists on installing, prior to the excavation of a length of tunnel, an array of pipes forming a kind of ‘umbrella’ above the area to be excavated. In some cases, these pipes are later used to inject grout to strengthen the ground and ‘connect’ the pipes. This system has proven to be very efficient in reducing tunnel convergence and water inflow when tunnelling through weak grounds. However, due to the geometrical and mechanical complexity of the problem, existing finite element frameworks are inappropriate to simulate tunnelling using this method. In this paper, a mathematical framework based on a homogenization technique to simulate ‘grouted pipe‐roofing reinforced ground’ and its implementation into a 3‐D finite element programme that can consider stage construction situations are presented. The constitutive model developed allows considering the main design parameters of the problem and only requires geometrical and mechanical properties of the constituents. Additionally, the use of a homogenization approach implies that the generation of the finite element mesh can be easily produced and that re‐meshing is not required as basic geometrical parameters such as the orientation of the pipes are changed. The model developed is used to simulate tunnelling with the grouted pipe‐roofing reinforcement method. From the analyses, the effects of the main design parameters on the elastic and the elastoplastic analyses are considered. Copyright © 2004 John Wiley & Sons, Ltd.     

A 4th order fabric tensor approach applied to granular media

In this paper, we propose a micromechanical approach of the behavior of granular media, which takes into account the anisotropy by means of a fourth order fabric tensor. The proposed approach is implemented in an homogenization scheme based on Voigt and Reuss localization assumption. The fabric tensor-based approach is then combined with a new kinematic localization rule and yields a general homogenization scheme for anisotropic granular media.     

Elastoplastic deformation around underground openings under biaxial initial stress field

Estimation of elastoplastic deformation around an underground opening induced by the excavation of it, especially displacement and strain field in plastic region, is presented in this paper, as well as the formulation for calculating the displacement and strain in the plastic region around the underground opening by the coupled Boundary Element Method - Characteristics Method (BEM-CM). In this method, the non-associated flow rule is adopted to calculate the displacement and strain field in the plastic region, which is determined by the integration of the displacement along characteristics lines under the boundary condition of the elastic displacement on an elastoplastic interface analysed. It is shown that this method is one of the accurate and effective methods for estimating not only the shape and extent of the plastic region but also the state of the displacement and strain in the plastic region around the underground opening, comparing the theoretical solution with numerical results by this method for a circular opening under hydrostatic initial stress condition. Furthermore, this method is applied to rectangular and horse-shoe shaped openings and the characteristics of the strain field in the plastic region are discussed. © 1997 by John Wiley & Sons, Ltd.