Elasto‐plastic seepage‐induced stresses due to tunneling

When tunneling is carried out beneath the groundwater table, hydraulic boundary is altered, resulting in seepage entering into the tunnel. The development of flow into the tunnel induces seepage stresses in the ground and the lining is subjected to additional loads. This can often cause fine particles to move, which clog the filter resulting in the long‐term hydraulic deterioration of the drainage system. However, the effect of seepage force is generally not considered in the analysis of tunnel. While several elastic solutions have been proposed by assuming seepage in an elastic medium, stress solutions have not been considered for the seepage force in a porous elasto‐plastic medium. This paper documents a study that investigates the stress behavior, caused by seepage, of a tunnel in an elasto‐plastic ground and its effects on the tunnel and ground. New elasto‐plastic solutions that adopt the Mohr–Coulomb failure criterion are proposed for a circular tunnel under radial flow conditions. A simple solution based on the hydraulic gradient obtained from a numerical parametric study is also proposed for practical use. It should be noted that the simple equation is useful for acquiring additional insight into a problem on a tunnel under drainage, because only a minimal computational effort is needed and considerable economic benefits can be gained by using it in the preliminary stage of tunnel design. The proposed equations were partly validated by numerical analysis, and their applicability is illustrated and discussed using an example problem. Comments on the tunnel analysis are also provided. Copyright © 2010 John Wiley & Sons, Ltd.     

Elasto‐plastic model for cement‐treated sand

This paper presents an elasto‐plastic model for non‐linear analyses of cement‐treated sand. Various laboratory tests were systematically carried out to investigate the pre‐peak and post‐peak behaviours of a cement‐treated sand. On the basis of these experimental results, the new model was built within the framework of a relatively simple elasto‐plastic theory. Two failure criteria are employed to express tensile and shear failure characteristics observed in the experimental results of the cement‐treated sand. The proposed model can describe strain‐hardening and strain‐softening responses under both failure modes. In the strain‐softening rules, the smeared crack concept is used, and a characteristic length is considered to avoid the issue of mesh‐size dependency. Since the failure criterion and strain‐hardening/softening rules are based on the experimental evidences, the model is relatively easy to understand and the parameters used in the model have clear physical meaning. The proposed model was applied to simulate the behaviour of cement‐treated sand in various laboratory tests, allowing for a reasonable comprehensive evaluation. It was demonstrated that the proposed model is suitable for describing both the tensile and shear failure behaviours of cement‐treated sand. Copyright © 2006 John Wiley & Sons, Ltd.     

Bifurcation and instability modelling by a multimechanism elasto‐plastic model

The bifurcation and instability conditions in geomechanics are closely related to the elasto‐plastic behaviour. In this paper the potential of a multimechanism elasto‐plastic model to predict various modes of failure is examined. First, a brief overview for the essential aspects of the constitutive model and the development of the elasto‐plastic constitutive matrix for this model are presented. Then, numerical simulations of different drained and undrained paths in the axisymmetric and plane‐strain conditions for the Hostun sand are illustrated. These examples confirm the capacity of the model to reproduce instability and strain localization phenomena. The obtained response is in agreement with experimental observations, theoretical developments and numerical analyses existing in the literature. Copyright © 2007 John Wiley & Sons, Ltd.     

Delayed plastic model for time‐dependent behaviour of materials

A delayed plastic model, based on the theory of plasticity, is proposed to represent the time‐dependent behaviour of materials. It is assumed in this model that the stress can lie outside the yield surface and the conjugate stress called static stress is defined on the yield surface. The stress–strain relation is calculated based on the plastic theory embedding the static stress. Thus, the stress–strain relation of the model practically corresponds to that of the inviscid elastoplastic model under fairly low rate deformation. The delayed plastic model is coupled with the Cam‐clay model for normally consolidated clays. The performance of the model is then examined by comparing the model predictions with reported time‐dependent behaviour of clays under undrained triaxial conditions. It is shown that the model is capable of predicting the effect of strain rate during undrained shear and the undrained creep behaviour including creep rupture. Copyright © 2001 John Wiley & Sons, Ltd.     

Three dimensional elasto‐plastic interface for embedded beam elements with interaction surface for the analysis of lateral loading of piles

This paper presents a numerical formulation for a three dimensional elasto‐plastic interface, which can be coupled with an embedded beam element in order to model its non‐linear interaction with the surrounding solid medium. The formulation is herein implemented for lateral loading of piles but is able to represent soil‐pile interaction phenomena in a general manner for different types of loading conditions or ground movements. The interface is formulated in order to capture localized material plasticity in the soil surrounding the pile within the range of small to moderate lateral displacements. The interface is formulated following two different approaches: (i) in terms of beam degrees of freedoms; and (ii) considering the displacement field of the solid domain. Each of these alternatives has its own advantages and shortcomings, which are discussed in this paper. The paper presents a comparison of the results obtained by means of the present formulation and by other well‐established analysis methods and test results published in the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

Discrete modelling of time‐dependent rockfill behaviour

This modelling study deals with the time‐dependent behaviour of rockfill media, which is of particular interest during the life of rockfill dams. Breakage of rock blocks and crack propagation are the main processes responsible for rockfill creep and collapse. The modelling procedure presented here is performed on two scales: on the rock block scale, where the grain is taken to be an assembly of rigid particles initially endowed with cohesive bonds, and on the rockfill scale, which is taken to involve a set of breakable grains interacting via contact and friction processes. The grain breakage process is described in term of a thermodynamically consistent damage interface model, where the damage is a gradual delayed process. This model was implemented in a non‐smooth contact dynamics code. The effects of the main parameters involved were analysed by performing numerical studies. The ability of the model to predict the creep behaviour of rockfill media is confirmed by presenting several simulations. Copyright © 2008 John Wiley & Sons, Ltd.     

Poro‐elasto‐plastic response of an unconsolidated formation confined with stiff seal rocks under radial injection

In this article, we evaluate geomechanics of fluid injection from a fully penetrating vertical well into an unconsolidated formation confined with stiff seal rocks. The coupled behavior of an isotropic, homogeneous sand layer is studied under injection pressures that are high enough to induce plasticity yet not fracturing. Propagation of the significant influence zone surrounding the injection borehole, quantified by the extent of the plastic domain in the elasto‐plastic model, is examined for the first time. First, a new fully coupled axisymmetric numerical model is developed. A comprehensive assessment is performed on pore pressures, stresses/strains, and failure planes during the entire transient period of an injection cycle. Results anticipate existence of five distinctive zones in terms of plasticity state: liquefaction at wellbore; two inner plastic domains surrounding the wellbore, where failure occurs along two planes and major principal stress is in vertical direction; remaining of the plastic domain, where formation fails along one plane and major principal stress is in radial direction; and a non‐plastic region. Failure mechanism at the wellbore is found to be shear followed by liquefaction. Next, a novel methodology is proposed based on which new weakly coupled poro‐elasto‐plastic analytical solutions are derived for all three stress/strain components. Unlike previous studies, extension of the plastic zone is obtained as a function of injection pressure, incorporating plasticity effects on the subsequent elastic domain. Solutions, proven to be a good approximation of numerical simulations, offer a huge advantage as the run time of coupled numerical simulations is considerably long. Copyright © 2016 John Wiley & Sons, Ltd.     

Stress‐dependent hardening and failure surfaces of dry sand

During several triaxial compression experiments on plastic hardening, softening, and failure properties of dense sand specimens, it was found on various stress paths that the size of the failure surface was not constant. Instead, it changed depending on the current state of hydrostatic pressure. This finding is in contrast to the standard opinion consisting of the fact that the failure surface remains constant, once it has been reached during an experiment or in situ. In general, the behaviour of cohesionless granular‐material‐like sand is somehow characterised in between fluid and solid, where the solid behaviour results from the angle of internal friction and the confining pressure. Although the friction angle is an intrinsic material property, the confining pressure varies with the boundary conditions, thus defining different solid properties like plastic hardening, softening, and also failure. Based on our findings, it was the goal of the present contribution to introduce an improved setting for the plastic strain hardening and softening behaviour including the newly found yield properties at the limit state. For the identification of the material parameters, a complete triaxial experimental analysis of the tested sand is given. The overall elasto‐plasticity concept is validated by numerical computations of several laboratory foundation‐ and slope‐failure experiments. The performance of the proposed approach is compared with the standard concept of a constant failure surface, where the corresponding yield surfaces are understood as contours of equivalent plastic work or plastic strain. Copyright © 2012 John Wiley & Sons, Ltd.     

Non‐isothermal plasticity model for cyclic behaviour of soils

On the one hand, it has been observed that liquefaction‐induced shear deformation of soils accumulates in a cycle‐by‐cycle pattern. On the other hand, it is known that heating could induce plastic hardening. This study deals with the constitutive modelling of the effect that heat may have on the cyclic mechanical properties of cohesive soils, a relatively new area of interest in soil mechanics. In this paper, after a presentation of the thermo‐mechanical framework, a non‐isothermal plasticity cyclic model formulation is presented and discussed. The model calibration is described based on data from laboratory sample tests. It includes numerical simulations of triaxial shear tests at various constant temperatures. Then, the model predictions are compared with experimental results and discussed in the final section. Both drained and undrained loading conditions are considered. The proposed constitutive model shows good ability to capture the characteristic features of behaviour. Copyright © 2007 John Wiley & Sons, Ltd.     

General coupling extended multiscale FEM for elasto‐plastic consolidation analysis of heterogeneous saturated porous media

This paper presents a general coupling extended multiscale FEM (GCEMs) for solving the coupling problem of elasto‐plastic consolidation of heterogeneous saturated porous media. In the GCEMs, the numerical multiscale base functions for the solid skeleton and fluid phase of the coupling system are all constructed on the basis of the equivalent stiffness matrix of the unit cell, which not only contain the interaction between the solid and fluid phases but also consider the time effect. Furthermore, in order to improve the computational accuracy for two‐dimensional problems, a multi‐node coarse element strategy for the GCEMs is proposed, and a two‐scale iteration algorithm for the elasto‐plastic consolidation analysis is developed. Some one‐dimensional and two‐dimensional homogeneous and heterogeneous numerical examples are carried out to validate the proposed method through the comparison with the coupling multiscale FEM and standard FEM. Numerical results show that the newly developed GCEMs can almost preserve the same convergent property as the standard FEM and also possesses the advantages of high computational efficiency. In addition, the GCEMs can be easily applied to other coupling multifield and multiphase transient problems. Copyright © 2014 John Wiley & Sons, Ltd.     

Dominant grains network and behavior of sand–silt mixtures: stress–strain modeling

Stress–strain modeling of sand–silt mixtures is important in the analysis and design of earth structures. In this paper, we develop a stress–strain model that can predict the behavior of sand–silt mixtures with any amount of fines content. This model is based on a micromechanics approach, which involves mean‐field assumptions. For the mixtures with low amount of fines, the mechanical behavior is dominated by sand grains network. On the other hand, for the mixtures with high amount of fines, the mechanical behavior is dominated by silt grains network. Using this concept of dominant grains network, the behavior of mixtures with any amount of fines can be predicted from knowing the behavior of sand and silt, alone. We also modeled the critical state friction angle, critical state void ratio, and elastic stiffness for the mixtures as a function of fines content. The applicability of this developed stress–strain model is shown by comparing the simulated and measured results for two different types of sand–silt mixtures with full range of fines content. Copyright © 2012 John Wiley & Sons, Ltd.     

Localized failure of fibre‐reinforced elastic–plastic materials subjected to plane strain loading

We consider discontinuous bifurcations as the indicator of a localized failure for a class of composites that are characterized by elastic fibres reinforcing an elastic–plastic matrix. A macroscopic tangent stiffness tensor for the fibre‐reinforced composite is developed by consistently homogenizing the contribution of fibres in a spherical representative volume element. Analytical solutions are derived for the critical hardening modulus and corresponding bifurcation directions for the case of plane strain loading. Properties of the solutions are further illustrated on the example of the non‐associated Drucker–Prager model at onset of yielding. Results show that presence of fibres decreases the critical hardening modulus, thus inhibiting the onset of strain localization. The rate of decrease in the critical hardening modulus is the highest for pure shear, followed by uniaxial tension, uniaxial compression, biaxial tension and biaxial compression. The main fibre parameters that control the onset of strain localization are their volumetric content and their stiffness modulus whereby very stiff fibres can produce the most significant decrease in the critical hardening modulus, especially for the state of biaxial tension. The critical hardening modulus for the non‐associated Drucker–Prager model exhibits a full range of localization modes including compaction bands, dilation bands, and transition in the form of shear bands regardless of the presence of fibres. Presence of fibres affects bifurcation directions, except in the case when Poisson's ratio of the matrix is equal to 0.25. The results demonstrate stabilizing effects of fibres by which they provide the control against the onset of strain localization. Copyright © 2006 John Wiley & Sons, Ltd.     

High‐P tectono‐metamorphic evolution of mylonites from the Variscan basement of the Northern Apennines,Italy

Strain localization within shear zones may partially erase the rock fabric and the metamorphic assemblage(s) that had developed before the mylonitic event. In poly‐deformed basements, the loss of information on pre‐kinematic phases of mylonites hinders large‐scale correlations based on tectono‐metamorphic data. In this study, devoted to a relict unit of Variscan basement reworked within the nappe stack of the Northern Apennines (Italy), we investigate the possibility to reconstruct a complete pressure (P)–temperature (T)–deformation (D) path of mylonitic micaschist and amphibolite by integrating microstructural analysis, mineral chemistry and thermodynamic modelling. The micaschist is characterized by a mylonitic fabric with fine‐grained K‐white mica and chlorite enveloping mica‐fishes, quartz, and garnet pseudomorphs. Potassic white mica shows Mg‐rich cores and Mg‐poor rims. The amphibolite contains green amphibole+plagioclase+garnet+quartz+ilmenite defining S₁ with a superposed mylonitic fabric localized in decimetre‐ to centimetre‐scale shear zones. Garnet is surrounded by an amphibole+plagioclase corona. Phase diagram calculations provide P–T constraints that are linked to the reconstructed metamorphic‐deformational stages. For the first time an early high‐P stage at >11 kbar and 510°C was constrained, followed by a temperature peak at 550–590°C and 9–10 kbar and a retrograde stage (<475°C, <7 kbar), during which ductile shear zones developed. The inferred clockwise P–T–D path was most likely related to crustal thickening by continent‐continent collision during the Variscan orogeny. A comparison of this P–T–D path with those of other Variscan basement occurrences in the Northern Apennines revealed significant differences. Conversely, a correlation between the tectono‐metamorphic evolution of the Variscan basement at Cerreto pass, NE Sardinia and Ligurian Alps was established.     

Experimental study and constitutive modelling of elasto‐plastic damage in heat‐treated mortar

This study investigates the effect of a heat‐treatment upon the thermo‐mechanical behaviour of a model cement‐based material, i.e. a normalized mortar, with a (w/c) ratio of 0.5. First, a whole set of varied experimental results is provided, in order to either identify or validate a thermo‐mechanical constitutive model, presented in the second paper part. Experimental responses of both hydraulic and mechanical behaviour are given after different heating/cooling cycling levels (105, 200, 300, 400^°C). The reference state, used for comparison purposes, is taken after mass stabilization at 60^°C. Typical uniaxial compression tests are provided, and original triaxial deviatoric compressive test responses are also given. Hydraulic behaviour is identified simultaneously to triaxial deviatoric compressive loading through gas permeability K_gas assessment. K_gas is well correlated with volumetric strain evolution: gas permeability increases hugely when ε_v testifies of a dilatant material behaviour, instead of contractile from the test start. Finally, the thermo‐mechanical model, based on a thermodynamics approach, is identified using the experimental results on uniaxial and triaxial deviatoric compression. It is also positively validated at residual state for triaxial deviatoric compression, but also by using a different stress path in lateral extension, which is at the origin of noticeable plasticity. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

Age and P–T conditions of the Gridino‐type eclogite in the Belomorian Province,Russia

Although eclogites in the Belomorian Province have been regarded as Archean in age and among the oldest in the world, there are also multiple studies that have proposed a Paleoproterozoic age. Here, we present new data for the Gridino‐type eclogites, which occur as boudins and metamorphosed dykes within tonalite–trondhjemite–granodiorite gneisses. Zircon from these eclogites has core and rim structures. The cores display high Th/U ratios (0.18–0.45), negative Eu anomalies and strong enrichment in HREE, and have Neoarchean U–Pb ages of c. 2.70 Ga; they are interpreted to be magmatic in origin. Zircon cores have δ¹⁸O of 5.64–6.07‰ suggesting the possibility of crystallization from evolved mantle‐derived magmas. In contrast, the rims, which include the eclogite facies minerals omphacite and garnet, are characterized by low Th/U ratios (<0.035) and flat HREE patterns, and yield U–Pb ages of c. 1.90 Ga; they are interpreted to be metamorphic in origin. Zircon rims have elevated δ¹⁸O of 6.23–6.80‰, which was acquired during eclogite facies metamorphism. Based on petrography and phase equilibria modelling, we recognize a prograde epidote amphibolite facies mineral assemblage, the peak eclogite facies mineral assemblage and a retrograde high‐P amphibolite facies mineral assemblage. The peak metamorphic conditions of 695–755°C at >18 kbar for the Gridino‐type eclogites suggest an apparent thermal gradient of <39–42°C/kbar for the Lapland–Kola collisional orogeny.     

Formation of bed‐scale spatial patterns in dolomite abundance during early dolomitization: Part I. Mechanisms and feedbacks revealed by reaction–transport modelling

Dolomites of varied ages exhibit metre‐scale nested patterns of lateral periodic variation in permeability and porosity and, by inference, dolomite abundance as most examples are 100% dolomite. Two‐dimensional reaction–transport modelling simulations of bed‐scale dolomitization were used to assess whether those patterns in dolomite abundance could form during near‐surface replacement dolomitization. Simulations used a 2 m high and 18 m long model domain, a low‐Mg calcite grainstone precursor and an evaporated Mississippian seawater brine (430 parts per thousand salinity) as the dolomitizing fluid. The domain was initially populated with random variations in porosity and/or grain size. Results reveal that spatial patterns in dolomite abundance emerge when there is as little as 1% dolomite formed, with similarities between the modelled patterns and outcrop‐documented patterns. The nested patterns include a near‐random component that constitutes ≤40% of the total variance, short‐range correlation ranging from 1·5 to 3·3 m and a longer‐range periodic trend with a wavelength up to 6·5 m. The emergence of pattern in dolomite abundance is the result of an autogenic self‐organizing phenomenon. It is triggered by variation in initial calcite reactive surface area that occurs due to the random heterogeneities in initial porosity and/or grain sizes. The pattern develops due to a combination of kinetic disequilibrium reactions (dolomite precipitation and calcite dissolution) and positive feedbacks between dolomite growth, calcite dissolution and fluid flow. Flow is around loci of higher dolomite, lower porosity and higher reactive surface areas, but through loci of lower dolomite, higher porosity and lower reactive surface areas. The resulting less porous/more dolomite and more porous/less dolomite structures at the metre‐scale arise from those localized interactions. This self‐organizing mechanism for pattern formation constitutes a new model for geochemical self‐organization during dolomitization and is the only self‐organization model that is proven applicable to the formation of metre‐scale patterns during early, near‐surface dolomitization.