A Generalized Backward Euler algorithm for the numerical integration of an isotropic hardening elastoplastic model for mechanical and chemical degradation of bonded geomaterials

An extended version of the classical Generalized Backward Euler (GBE) algorithm is proposed for the numerical integration of a three‐invariant isotropic‐hardening elastoplastic model for cemented soils or weak rocks undergoing mechanical and non‐mechanical degradation processes. The restriction to isotropy allows to formulate the return mapping algorithm in the space of principal elastic strains. In this way, an efficient and robust integration scheme is developed which can be applied to relatively complex yield surface and plastic potential functions. Moreover, the proposed algorithm can be linearized in closed form, thus allowing for quadratic convergence in the global Newton iteration. A series of numerical experiments are performed to illustrate the accuracy and convergence properties of the algorithm. Selected results from a finite element analysis of a circular footing on a soft rock layer undergoing chemical weathering are then presented to illustrate the algorithm performance at the boundary value problem level. Copyright © 2002 John Wiley & Sons, Ltd.     

Stress–dilatancy based modelling of granular materials and extensions to soils with crushable grains

Stress–dilatancy relations have played a crucial role in the understanding of the mechanical behaviour of soils and in the development of realistic constitutive models for their response. Recent investigations on the mechanical behaviour of materials with crushable grains have called into question the validity of classical relations such as those used in critical state soil mechanics. In this paper, a method to construct thermodynamically consistent (isotropic, three‐invariant) elasto‐plastic models based on a given stress–dilatancy relation is discussed. Extensions to cover the case of granular materials with crushable grains are also presented, based on the interpretation of some classical model parameters (e.g. the stress ratio at critical state) as internal variables that evolve according to suitable hardening laws. Copyright © 2004 John Wiley & Sons, Ltd.     

A hypoplastic macroelement formulation for single batter piles in sand

Batter piles are widely used in geotechnical engineering when substantial lateral resistance is needed or to avoid the interference with existing underground constructions. Nevertheless, there is a lack of fast numerical tools for nonlinear soil‐structure interactions problems for this type of foundation. A novel hypoplastic macroelement is proposed, able to reproduce the nonlinear response of a single batter pile in sand under monotonic and cyclic static loadings. The behavior of batter piles (15°, 30°, and 45°) is first numerically investigated using 3D finite element modeling and compared with the behavior of vertical piles. It is shown that their response mainly depends on the pile inclination and the loading direction. Then, starting from the macroelement for single vertical piles in sand by Li et al (Acta Geotechnica, 11(2):373‐390, 2016), an extension is proposed to take into account the pile inclination introducing simple analytical equations in the expression describing the failure surface. 3D finite element numerical models are adopted to validate the macroelement that is proven able to reproduce the nonlinear behavior in terms of global quantities (forces‐displacements) and to significantly reduce the necessary computational time.     

Time-dependent ground movements induced by shield tunneling in soft clay: a parametric study

In this work, the effects of coupled hydromechanical (consolidation) processes associated with shield tunneling excavation in soft clays are investigated with particular attention to the prediction of ground movements at the ground surface. A series of 2d FE analyses have been carried out in parametric form in order to investigate the effects of tunnel excavation velocity relative to the soil consolidation rate and the hydraulic boundary conditions at the tunnel boundary. The shield advancement process has been simulated with a simplified procedure incorporating both volume loss and ovalization of the tunnel section. In order to investigate the relative importance of soil consolidation during the excavation process, different characteristic times for the tunnel face advancement and for the consolidation process around the tunnel have been considered, for the two limiting conditions of fully permeable liner and impervious liner. The potential damage induced by the tunnel excavation on existing structures, based on computed ground surface distortions and horizontal deformations, has been found to vary significantly with time during the consolidation process. The results of the simulations allowed to obtain useful information on the minimum tunnel face advancement speed for which the assumption of fully undrained conditions for the soil during the excavations is acceptable, as well as on the speed range for which solving the fully coupled hydromechanical problem is necessary.     

A constitutive model for bonded geomaterials subject to mechanical and/or chemical degradation

The mechanical behaviour of bonded geomaterials is described by means of an elastoplastic strain‐hardening model. The internal variables, taking into account the ‘history’ of the material, depend on the plastic strains experienced and on a conveniently defined scalar measure of damage induced by weathering and/or chemical degradation. For the sake of simplicity, it is assumed that only internal variables are affected by mechanical and chemical history of the material. Despite this simplifying assumption, it can be shown that many interesting phenomena exhibited by weathered bonded geomaterials can be successfully described. For instance, (i) the transition from brittle to ductile behaviour with increasing pressure of a calcarenite with collapsing internal structure, (ii) the complex behaviour of chalk and other calcareous materials in oedometric tests, (iii) the chemically induced variation of the stress and strain state of such kind of materials, are all phenomena that can be qualitatively reproduced. Several comparisons with experimental data show that the model can capture the observed behaviour also quantitatively. Copyright © 2003 John Wiley & Sons, Ltd.     

The effects of Soil–Structure Interaction on a reinforced concrete viaduct

This paper presents a numerical strategy to model a three‐pier viaduct made of prestressed concrete. The viaduct was tested pseudodynamically in ELSA laboratory (JRC Ispra, Italy). During the experimental campaign, only the three piers were tested, whereas the behaviour of the deck was simulated using the finite element method. The first part of the paper presents a numerical model of the viaduct based on the Timoshenko multifibre beam elements and non‐linear constitutive laws. Comparisons with the experimental results show the good performance of the approach. In the second part, a parametric study is carried out showing the influence of Soil–Structure Interaction (SSI). Various types of soils are considered using a recently developed macro‐element representing a rigid shallow foundation. The macro‐element is suitable for dynamic (seismic) loadings and it takes into account the plasticity of the soil, the uplift of the foundation, P?θ effects and the radiative damping. Finally, the numerical results are compared with the ones coming from a classical engineering approach using linear elastic springs at the base of the piers. This comparison shows that SSI is a complex phenomenon inducing displacements and internal forces in the structure that are difficult to predict with the linear approach. Based on the results obtained in this paper, it seems now possible to use this approach to investigate numerically the behaviour of a wider variety of configurations. Copyright © 2010 John Wiley & Sons, Ltd.     

Directional response of a reconstituted fine‐grained soil—Part I: experimental investigation

This paper discusses the results of a large experimental program designed to investigate in a systematic manner the main features of the incremental response of fine‐grained soils. The results are obtained from triaxial stress probing experiments carried out on a French silty clay (Beaucaire Marl). All the tests have been performed on reconstituted specimens, normally consolidated to an initial state which is either isotropic or anisotropic. In the interpretation of the experimental results, extensive use is made of the concept of strain response envelope. The response envelopes obtained for different stress increment magnitudes are remarkably consistent with each other and indicate an inelastic and irreversible material response, i.e. a strong dependence on the stress increment direction, also at relatively small strain levels. A companion paper (Int. J. Numer. Anal. Meth. Geomech., this issue, 2006) assesses the performance of some advanced constitutive models in reproducing the behaviour of reconstituted Beaucaire Marl as observed in this experimental program. Copyright © 2006 John Wiley & Sons, Ltd.     

On the effects of pore water pressure buildup and dissipation on the seismic performance of a propped r.c. diaphragm wall in sand

Acta Geotechnica - The study concerns the analysis of a retaining structure composed by a couple of r.c. diaphragm walls propped at the crest in loose and medium-dense, variably saturated sand...     

Spatially distributed rainfall thresholds for the initiation of shallow landslides

The evaluation of the combined influence of rainfall patterns (in terms of mean intensity and duration) and the geomorphological and mechanical characteristics of hillslopes on their stability conditions is a major goal in the assessment of the shallow landslide triggering processes. Geographic Information Systems (GIS) represent an important tool to develop models that combine hydrological and geomechanical analyses for the evaluation of slope stability, as they allow to combine information concerning rainfall characteristics with topographic and mechanical properties of the slopes over wide areas. In this paper, a GIS-based code is developed to determine physically based intensity/duration rainfall thresholds at the local scale. Given the rainfall duration and the local geometric, hydrological and mechanical characteristics of the slopes, the code evaluates the spatial distribution of the minimum rainfall intensity that triggers shallow landslides and debris flows over a given area. The key feature of the code is the capability of evaluating the time t _p required to reach the peak pore pressure head on the failure surface and computing the corresponding critical intensity/duration thresholds based on post-event peak pore pressures. The reliability of the model is tested using a set of one-dimensional analyses, comparing the physically based thresholds obtained for three different slopes with some empirical rainfall thresholds. In a log–log scale, the thresholds provided by the model decrease linearly with increased rainfall duration and they are bracketed by the empirical thresholds considered. Finally, an example of application to a study area of the Umbria region in central Italy is presented, describing the capability of the model of providing site-specific thresholds for different rainfall scenarios.