A simple hypoplastic model for normally consolidated clay

The paper presents a simple constitutive model for normally consolidated clay. A mathematical formulation, using a single tensor-valued function to define the incrementally nonlinear stress–strain relation, is proposed based on the basic concept of hypoplasticity. The structure of the tensor-valued function is determined in the light of the response envelope. Particular attention is paid towards incorporating the critical state and to the capability for capturing undrained behaviour of clayey soils. With five material parameters that can be determined easily from isotropic consolidation and triaxial compression tests, the model is shown to provide good predictions for the response of normally consolidated clay along various stress paths, including drained true triaxial tests and undrained shear tests.     

In-plane and anti-plane strong shaking of soil systems and structures

The concept of in-plane and anti-plane shaking is introduced with a rigid block on a plane surface with Coulomb friction. Using a hypoplastic constitutive relation to model the mechanical behaviour of the soil, numerical solutions for a rigid block on a thin dry or saturated soil layer are obtained. The coupled nature of dynamic problems involving granular materials is shown, i.e. the motion of the block changes the soil state—skeleton stresses and density—which in turn affects the block motion. Motions of the block as well as soil response can be more realistically calculated by the new model. The same constitutive equation is applied to the numerical simulation of the propagation of plane waves in homogeneous and layered level soil deposits induced by a wave coming from below. Experiments with a novel laminar shake box as well as real seismic records from well-documented sites during strong earthquakes are used to verify the adequacy of the hypoplasticity-based numerical model for the prediction of soil response during strong earthquakes. The response of a homogeneous earth dam subjected to in-plane and anti-plane shaking is investigated numerically. In-plane and anti-plane shaking is shown to cause nearly the same spreading of a sand dam under drained conditions, whereas under undrained conditions anti-plane shaking causes stronger spreading of the dam. The dynamic behaviour of a breakwater founded on rockfill and soft clay during the 1995 Kobe earthquake is back-calculated to show the good performance of the proposed numerical model also with a structure. Section 9 deals with buildings on mattresses of densified cohesionless soils or fine-grained soils with granular columns, slopes with ‘hidden’ dams and structures on piles traversing clayey slopes to show the suitability of hypoplasticity-based models for the earthquake-resistant design and safety assessment of geotechnical systems.     

Effects of installation on the cyclic axial behaviour of suction buckets in sandy soils

Suction buckets differ with their easy and cost-efficient installation technique from other foundation types for offshore wind turbines. For successful completion of their installation process, suction is essential, but the imposed seepage leads to the changes in states of the soil in and around the bucket. Especially, a loosening of soil inside the bucket affects the load carrying behaviour of bucket subjected to repetitive loading resulting from environmental conditions. In this study, the behaviour of buckets under cyclic axial compressive loads with considering a possible loosening and related changes in permeability of soil inside the bucket is investigated numerically. In the framework of finite element analysis, a fully coupled two-phase model and a hypoplastic constitutive model are used to describe the saturated sandy soil behaviour under repetitive loading. The porosity-permeability variation is taken into account by Kozeny–Carman relationship. Special attention is dedicated to load carrying behaviour of bucket top plate, inner and outer skirt as well as base and their changes resulting from a loosening of soil inside the bucket with variable aspect ratio. For this purpose, cyclic axial compressive loads which cause an attenuation and progressive failure of soil-bucket system response are considered. The main findings on the changes in load carrying behaviour of bucket are presented and discussed.     

Influence of initial density of cohesionless soil on evolution of passive earth pressure

This paper focuses on the influence of the initial void ratio on the evolution of the passive earth pressure and the formation of shear zones in a dry sand body behind a retaining wall. For the numerical simulation a rigid and very rough retaining wall undergoing a horizontal translation against the backfill is considered. The essential mechanical properties of cohesionless granular soil are described with a micro-polar hypoplastic model which takes into account stresses and couple stresses, pressure dependent limit void ratios and the mean grain size as a characteristic length. Numerical investigations are carried out with an initially medium dense and initially loose sand using a homogeneous and random distribution of the initial void ratio. The geometry of calculated shear zones is discussed and compared with a corresponding laboratory model test.     

Numerical study of the 3D failure envelope of a single pile in sand

The paper presents a comprehensive study of the failure envelope (or capacity diagram) of a single elastic pile in sand. The behavior of a pile subjected to different load combinations is simulated using a large number of finite element numerical calculations. The sand is modeled using a constitutive law based on hypoplasticity. In order to find the failure envelope in the three-dimensional space (i.e. horizontal force H, bending moment M and vertical force V), the radial displacement method and swipe tests are numerically performed. It is found that with increasing vertical load the horizontal bearing capacity of the pile decreases. Furthermore, the presence of bending moment on the pile head significantly influences the horizontal bearing capacity and the capacity diagram in the H–M plane manifests an inclined elliptical shape. An analytical equation providing good agreement with the 3D numerical results is finally proposed. The formula is useful for design purposes and the development of simplified modeling numerical strategies such as macro-element.     

Coupled hydro-mechanical model for partially saturated soils predicting small strain stiffness

In the paper, we present newly developed hydro-mechanical hypoplastic model for partially saturated soils predicting small strain stiffness. Hysteretic void ratio dependent water retention model has been incorporated into the existing hypoplastic model. This required thorough revision of the model structure to allow for the hydro-mechanical coupling dependencies. The model is formulated in terms of degree of saturation, rather than of suction. Subsequently, the small strain stiffness effects were incorporated using the intergranular strain concept modified for unsaturated conditions. New features included degree of saturation-dependent size of the elastic range and an updated evolution equation for the intergranular strain. The model has been evaluated using two comprehensive data sets on completely decomposed tuff from Hong-Kong and Zenos Kaolin from Iran. It has been shown that the modified intergranular strain formulation coupled with the hysteretic water retention model correctly reproduces the effects of both the stress and suction histories on small strain stiffness evolution. The model can correctly predict also different other aspects of partially saturated soil behaviour, starting from the very small strain range up to the asymptotic large-strain response.     

Capability of constitutive models to simulate soils with different OCR using a single set of parameters

Incorporation of void ratio as a state variable into constitutive models allows, in principle, to use a single set of parameters for soils with different OCRs. Two sets of experimental data on reconstituted clays are used for evaluation of three constitutive models of different complexity (Modified Cam clay model, 3SKH model, hypoplastic model for clays). Although all the models predict the influence of OCR correctly from the qualitative point of view, quantitative comparison using a suitable scalar error measure reveals merits and shortcomings of different models.     

Three-dimensional topology optimization for geotechnical foundations in granular soil

This article presents three-dimensional structural optimization in geotechnical engineering for foundations in granular soil. The general design (topology) of a shallow foundation is optimized with respect to its deformational behaviour within the service limit state. The SIMP (solid isotropic material with penalization) method is applied to optimize the distribution of foundation material. The soil is modelled as a hypoplastic material with a constitutive model suitable for optimization using finite element analysis. Two load cases are examined. The optimized topology is validated against two-dimensional optimization and 1g-model test results. The present study proves the applicability and shows the potential of topology optimization in geotechnical engineering.     

Liapunov instability of the hypoplastic model for soils

The hypoplastic constitutive model is known for its numerous application to the problems of soil mechanics and also for its excessive ratcheting. The paper shows that this deficiency can be interpreted as Liapunov instability in dynamic case. This is demonstrated with a simple one-dimensional swinger. Infinitesimally small as well as finite oscillations are analytically examined and the variability of stiffness upon a single cycle is considered. Several methods to circumvent the problem of ratcheting are discussed: implementation of ‘elastic’ range by means of so-called intergranular strain, usage of overlay modelling with parallel coupling of hypoplastic materials as originally proposed by Valanis in his endochronic theory and a special coupling of hypoplastic models based on comparison of partial stiffnesses.     

Adaptive integration of constitutive rate equations

In finite element calculations the constitutive model plays a key role. The evaluation of the stress response of the constitutive relation for a given strain increment, which is a time integration in the case of models of the rate type, is a typical sub task in such calculations. Adaptive behaviour of the time integration is essential to assure numerical stability and to control the accuracy of the solution. An adaptive second order semi-implicit method is developed in this paper. Its numerical behaviour is compared with an adaptive second order explicit scheme. The two proposed methods control the local error and guarantee numerical stability of the time integration. We include several numerical geotechnical element tests using hypoplasticity with intergranular strain. The element tests simulate the behaviour of a finite element method based on the displacement formulation.