ACMEG‐TS: A constitutive model for unsaturated soils under non‐isothermal conditions

This paper introduces an unconventional constitutive model for soils, which deals with a unified thermo‐mechanical modelling for unsaturated soils. The relevant temperature and suction effects are studied in light of elasto‐plasticity. A generalized effective stress framework is adopted, which includes a number of intrinsic thermo‐hydro‐mechanical connections, to represent the stress state in the soil. Two coupled constitutive aspects are used to fully describe the non‐isothermal behaviour. The mechanical constitutive part is built on the concepts of bounding surface theory and multi‐mechanism plasticity, whereas water retention characteristics are described using elasto‐plasticity to reproduce the hysteretic response and the effect of temperature and dry density on retention properties. The theoretical formulation is supported by comparisons with experimental results on two compacted clays. Copyright © 2008 John Wiley & Sons, Ltd.     

Review and enhancement of 3D concrete models for large‐scale numerical simulations of concrete structures

The present paper focuses on selected plasticity and damage‐plasticity models for describing the 3D material behavior of concrete. In particular, a plasticity model and a damage‐plasticity model are reviewed and evaluated. Based on the results of the evaluation, enhancements are proposed, aiming at improving the correspondence between predicted and observed material behavior and aiming at implementing a robust and efficient stress update algorithm in a finite element program for performing large‐scale 3D numerical simulations of concrete structures. The capabilities of the concrete models are demonstrated by 3D numerical simulations of benchmark tests with combined bending and torsional loading and combined compression and shear loading and by a large‐scale 3D finite element analysis of a model test of a concrete arch dam. Copyright © 2011 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.     

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.     

Directional response of a reconstituted fine‐grained soil—Part II: performance of different constitutive models

In this paper, the performance of different advanced constitutive models for soils is evaluated with respect to the experimentally observed behaviour of a soft reconstituted clay subject to a wide range of loading directions, see (presented in the companion paper). The models considered include a three‐surface kinematic hardening elastoplastic model; the CLoE hypoplastic model; a recently proposed K‐hypoplastic model for clays, and an enhanced version of the same model incorporating the concept of intergranular strain. A clear qualitative picture of the relative performance of the different models as a function of the loading direction is obtained by means of the incremental strain response envelopes. The definition of suitable error measures allows to obtain further quantitative information in this respect. For the particular initial conditions and loading programme considered in this study, the kinematic hardening and the enhanced K‐hypoplastic models appear to provide the best performance overall. Copyright © 2006 John Wiley & Sons, Ltd.     

Spectral analysis of nonlocal regularization in two‐dimensional finite element models

Strain‐softening in geomaterials often leads to ill‐posed boundary‐valued problems (BVP), which cannot be solved with finite element methods without introducing some kind of regularization such as nonlocal plasticity. Hereafter we propose to apply spectral analysis for testing the performance of nonlocal plasticity in regularizing ill‐posed BVP and producing mesh‐independent solutions when local plasticity usually fails. The spectral analysis consists of examining the eigenvalues and eigenvectors of the global tangential stiffness matrix of the incremental equilibrium equations. Based on spectral analysis, we propose a criterion for passing or failing the test of constitutive regularization in the context of BVP. If the eigenvalues of the tangential operator are all positive then the regularization succeeds, otherwise it fails and may not prevent artificial mesh‐dependent solutions from appearing. The approach is illustrated in the particular case of a biaxial compression with strain‐softening plasticity. In this particular case, local softening plasticity is found to produce negative eigenvalues in the tangential stiffness matrix, which indicates ill‐posed BVP. In contrast, nonlocal softening plasticity always produces positive eigenvalues, which regularizes ill‐posed BVP. The dominant eigenvectors, which generate localized deformation patterns, have a bandwidth independent of mesh size, provided that the mesh is fine enough to capture localization. These mesh‐independent eigenmodes explain why nonlocal plasticity produces numerical solutions that are mesh‐independent. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical simulations of simple shear with non‐coaxial soil models

This paper investigates the effects of a non‐coaxial model on simulated stress–strain behaviour of granular materials subject to simple shearing under various initial conditions. In most cases, a significant difference of predictions between coaxial and non‐coaxial modelling is found during the early stage in shearing. With the increase in shearing, non‐coaxial simulations approach and tend to coincide with coaxial simulations. It is also found that the roles of non‐coaxial modelling in simulating simple shear behaviour are considerably influenced by hardening rules, flow rules, initial static lateral pressure coefficients. In some cases, the non‐coaxial modelling gives a similar simulation as the coaxial modelling. In other cases, the non‐coaxial modelling decreases the hardening response or softening response of materials, compared with the coaxial modelling. Under certain conditions, the predicted peak strength of materials with non‐coaxial modelling is larger than that for coaxial modelling. Some of these observations can be attributed to the amount of principal stress rotation in various cases analysed. Others can be attributed to the difference between the directions of the non‐coaxial plastic flow and those for coaxial plastic flow. Copyright © 2005 John Wiley & Sons, Ltd.     

A unified approach to the implicit integration of standard,non‐standard and viscous plasticity models

It is shown how modern concepts to integrate the elasto‐plastic rate equations of standard plasticity via an implicit algorithm can be generalized to plasticity without an explicitly defined yield surface and to overstress‐type models of viscoplasticity, where the stress point can be located outside the loading surface. For completeness, a tangent operator is derived that is consistent with the update algorithm. Copyright © 2002 John Wiley & Sons, Ltd.     

A coupled chemo‐thermo‐hygro‐mechanical model of concrete at high temperature and failure analysis

A hierarchical mathematical model for analyses of coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes at high temperature is presented. The concretes are modelled as unsaturated deforming reactive porous media filled with two immiscible pore fluids, i.e. the gas mixture and the liquid mixture, in immiscible–miscible levels. The thermo‐induced desalination process is particularly integrated into the model. The chemical effects of both the desalination and the dehydration processes on the material damage and the degradation of the material strength are taken into account. The mathematical model consists of a set of coupled, partial differential equations governing the mass balance of the dry air, the mass balance of the water species, the mass balance of the matrix components dissolved in the liquid phases, the enthalpy (energy) balance and momentum balance of the whole medium mixture. The governing equations, the state equations for the model and the constitutive laws used in the model are given. A mixed weak form for the finite element solution procedure is formulated for the numerical simulation of chemo‐thermo‐hygro‐mechanical behaviours. Special considerations are given to spatial discretization of hyperbolic equation with non‐self‐adjoint operator nature. Numerical results demonstrate the performance and the effectiveness of the proposed model and its numerical procedure in reproducing coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes subjected to fire and thermal radiation. Copyright © 2005 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.     

An original semi‐discrete approach to assess gas conductivity of concrete structures

For civil engineering structures with a tightness role, structural permeability is a key issue. In this context, this paper presents a new proposition of a numerical modelling of leakage rate through a cracked concrete structure undergoing mode I cracking. The mechanical state of the material, considered in the framework of continuum mechanics based on finite element modelling, is described by means of the stress‐based nonlocal damage model which takes into account the stress state and provides realistic local mechanical fields. A semi‐discrete method based on the strong discontinuity approach to estimate crack opening is then considered in the post‐treatment phase. Using a Poiseuille's like relation, the coupling between the mechanical state of the material and its dry gas conductivity is performed. For validation purposes, an original experimental campaign is conducted on a dry concrete disc loaded in a splitting setup. During the loading, gas conductivity and digital image correlation analysis are performed. The comparison with the 3D experimental mechanical global response highlights the performance of the mechanical model. The comparison between crack openings measured by digital image correlation and estimated by the strong discontinuity method shows a good agreement. Finally, the results of the semi‐discrete approach coupled with the gas conductivity compared with experimental data show a good estimation of the structural conductivity. Consequently, if the mechanical problem is well modelled at the global scale, then the proposed approach provides good estimation of gas conductivity. Copyright © 2016 John Wiley & Sons, Ltd.     

A novel computational approach for large deformation and post‐failure analyses of segmental retaining wall systems

Segmental retaining wall (SRW) systems are commonly used in geotechnical practice to stabilize cut and fill slopes. Because of their flexibility, these systems can tolerate minor movements and settlements without incurring damage or crack. Despite these advantages, very few numerical studies of large deformations and post‐failure behavior of SRW systems are found in the current literature. Traditional numerical methods, such as the finite element method, suffer from mesh entanglement, thus are unable to simulate large deformations and flexible behavior of retaining wall blocks in SRW systems. To overcome the above limitations, a novel computational framework based on the smoothed particle hydrodynamics (SPH) method was developed to simulate large deformations and post‐failure behavior of soils and retaining wall blocks in SRW systems. The proposed numerical framework is a hybrid continuum/discontinuum approach that can model soil as an elasto‐plastic material and retaining wall blocks as independent rigid bodies associated with both translational and rotational degrees of freedom. A new contact model is proposed within the SPH framework to simulate the interaction between the soil and the blocks and between the blocks. As an application of the proposed numerical method, a two‐dimensional simulation of an SRW collapse was simulated and compared to experimental results conducted under the same conditions. The results showed that the proposed computational approach provided satisfactory agreement with the experiment. This suggests that the new framework is a promising numerical approach to model SRW systems. Copyright © 2014 John Wiley & Sons, Ltd.     

Interaction between different internal length scales in strain localization analysis of fully and partially saturated porous media—the 1‐D case

The paper analyses the interaction between two internal length scales during dynamic strain localization in multiphase porous materials. The first internal length is introduced in the mathematical model by the gradient‐dependent plasticity for the solid skeleton, while the second one is naturally contained in the multiphase model and is due to the seepage process of the water via Darcy's law, which induces a rate‐dependent behaviour of the solid skeleton. Numerical results of a one‐dimensional example of water saturated porous medium demonstrate the competing effect between these two length scales. The porous medium is here treated as a multiphase continuum, with the pores filled by water and air, the last one at constant atmospheric pressure. Copyright © 2005 John Wiley & Sons, Ltd.     

The effect of pre‐tectonic reaction and annealing extent on behaviour during subsequent deformation: insights from paired shear zones in the lower crust of Fiordland,New Zealand

Granulite facies pargasite orthogneiss is partially to completely reacted to garnet granulite either side of narrow (<20 mm) felsic dykes, in Fiordland, New Zealand, forming ~10–80 mm wide garnet reaction zones. The metamorphic reaction changed the abundance of minerals, and their shape and grain size distribution. The extent of reaction and annealing (temperature‐related coarsening and nucleation) is greatest close to the dykes, whereas further away the reaction is incomplete. As a consequence, grain size and the abundance of garnet decreases away from the felsic dykes over a few centimetres. The aspect ratios of clusters of S1 pyroxene and pargasite in the orthogneiss, which are variably reacted to post‐S1 garnet, decrease from high in the host, to near equidimensional close to the dyke. Post‐reaction deformation localized in the fine‐grained partially reacted areas. This produced a pattern of ‘paired’ shear zones located at the outer parts of the garnet reaction zone. Our study shows that grain size sensitive deformation occurs where the grain size is sufficiently reduced by metamorphic reaction. The weakening of the rock due to the change in grain size distribution outweighs the addition of nominally stronger garnet to the assemblage.     

Numerical modeling of hydraulic fracture propagation,closure and reopening using XFEM with application to in‐situ stress estimation

In this paper, a fully coupled model is developed for numerical modeling of hydraulic fracturing in partially saturated weak porous formations using the extended finite element method, which provides an effective means to simulate the coupled hydro‐mechanical processes occurring during hydraulic fracturing. The developed model is for short fractures where plane strain assumptions are valid. The propagation of the hydraulic fracture is governed by the cohesive crack model, which accounts for crack closure and reopening. The developed model allows for fluid flow within the open part of the crack and crack face contact resulting from fracture closure. To prevent the unphysical crack face interpenetration during the closing mode, the crack face contact or self‐contact condition is enforced using the penalty method. Along the open part of the crack, the leakage flux through the crack faces is obtained directly as a part of the solution without introducing any simplifying assumption. If the crack undergoes the closing mode, zero leakage flux condition is imposed along the contact zone. An application of the developed model is shown in numerical modeling of pump‐in/shut‐in test. It is illustrated that the developed model is able to capture the salient features bottomhole pressure/time records exhibit and can extract the confining stress perpendicular to the direction of the hydraulic fracture propagation from the fracture closure pressure. Copyright © 2016 John Wiley & Sons, Ltd.     

A high strain‐rate constitutive model for sand and its application in finite‐element analysis of tunnels subjected to blast

The paper presents a constitutive model for simulating the high strain‐rate behavior of sands. Based on the concepts of critical‐state soil mechanics, the bounding surface plasticity theory and the overstress theory of viscoplasticity, the constitutive model simulates the high strain‐rate behavior of sands under uniaxial, triaxial and multi‐axial loading conditions. The model parameters are determined for Ottawa and Fontainebleau sands, and the performance of the model under extreme transient loading conditions is demonstrated through simulations of split Hopkinson pressure bar tests up to a strain rate of 2000/s. The constitutive model is implemented in a finite‐element analysis software Abaqus to analyze underground tunnels in sandy soil subjected to internal blast loads. Parametric studies are conducted to examine the effect of relative density and type of sand and of the depth of tunnel on the variation of stresses and deformations in the soil adjacent to the tunnels. Copyright © 2012 John Wiley & Sons, Ltd.     

Chemo‐mechanics of cemented granular solids subjected to precipitation and dissolution of mineral species

This paper studies the chemo‐mechanics of cemented granular solids in the context of continuum thermodynamics for fluid‐saturated porous media. For this purpose, an existing constitutive model formulated in the frame of the Breakage Mechanics theory is augmented to cope with reactive processes. Chemical state variables accounting for the reactions between the solid constituents and the solutes in the pore fluid are introduced to enrich the interactions among the microstructural units simulated by the model (i.e., grains and cement bonds). Two different reactive processes are studied (i.e., grain dissolution and cement precipitation), using the chemical variables to describe the progression of the reactions and track changes in the size of grains and bonds. Finally, a homogenization strategy is used to derive the energy potentials of the solid mixture, adopting probability density functions that depend on both mechanical and chemical indices. It is shown that the connection between the statistics of the micro‐scale attributes and the continuum properties of the solid enables the mathematical capture of numerous mechanical effects of lithification and chemical deterioration, such as changes in stiffness, expansion/contraction of the elastic domain, and development of inelastic strains during reaction. In particular, the model offers an interpretation of the plastic strains generated by aggressive environments, which are here interpreted as an outcome of chemically driven debonding and comminution. As a result, the model explains widely observed macroscopic signatures of geomaterial degradation by reconciling the energetics of the deformation/reaction processes with the evolving geometry of the microstructural attributes. Copyright © 2015 John Wiley & Sons, Ltd.