Validation of a three‐dimensional constitutive model for nonlinear site response and soil‐structure interaction analyses using centrifuge test data

The capability of a bounding surface plasticity model with a vanishing elastic region to capture the multiaxial dynamic hysteretic responses of soil deposits under broadband (eg, earthquake) excitations is explored by using data from centrifuge tests. The said model was proposed by Borja and Amies in 1994 (J. Geotech. Eng., 120, 6, 1051‐1070), which is theoretically capable of representing nonlinear soil behavior in a multiaxial setting. This is an important capability that is required for exploring and quantifying site topography, soil stratigraphy, and kinematic effects in ground motion and soil‐structure interaction analyses. Results obtained herein indicate that the model can accurately predict key response data recorded during centrifuge tests on embedded specimens—including soil pressures and bending strains for structural walls, structures' racking displacements, and surface settlements—under both low‐ and high‐amplitude seismic input motions, which was achieved after performing only a basic material parameter calibration procedure. Comparisons are also made with results obtained using equivalent linear models and a well‐known pressure‐dependent multisurface plasticity model, which suggested that the present model is generally more accurate. The numerical convergence behavior of the model in nonlinear equilibrium iterations is also explored for a variety of numerical implementation and model parameter options. To facilitate broader use by researchers and practicing engineers alike, the model is implemented as a “user material” in ABAQUS Standard for implicit time stepping.     

A stepwise damping‐solvent extraction method for large‐scale dynamic soil–structure interaction analysis in time domain

Time‐domain analysis of dynamic soil–structure interaction based on the substructure method plays an increasing role in practical applications as compared with the frequency‐domain analysis. Efficient and accurate modelling of the unbounded soil or rock medium has been a key issue in such an analysis. This paper presents a subregional stepwise damping‐solvent extraction formulation for solving large‐scale dynamic soil–structure problems in the time domain. Accuracy and efficiency of the formulation are evaluated in detail for a classical problem involving a rigid strip foundation embedded in a half‐space. A practical large‐scale soil–structure interaction problem, which represents a high concrete gravity dam subjected to seismic load, is then analysed using the proposed method. Various responses of the dam, including time histories of the crest displacement and acceleration and contours of the peak principal stresses within the dam body, are presented. Comparisons are also made between these results with those obtained using other models for the unbounded medium. Copyright © 2007 John Wiley & Sons, Ltd.     

An efficient finite–discrete element method for quasi‐static nonlinear soil–structure interaction problems

An efficient finite–discrete element method applicable for the analysis of quasi‐static nonlinear soil–structure interaction problems involving large deformations in three‐dimensional space was presented in this paper. The present method differs from previous approaches in that the use of very fine mesh and small time steps was not needed to stabilize the calculation. The domain involving the large displacement was modeled using discrete elements, whereas the rest of the domain was modeled using finite elements. Forces acting on the discrete and finite elements were related by introducing interface elements at the boundary of the two domains. To improve the stability of the developed method, we used explicit time integration with different damping schemes applied to each domain to relax the system and to reach stability condition. With appropriate damping schemes, a relatively coarse finite element mesh can be used, resulting in significant savings in the computation time. The proposed algorithm was validated using three different benchmark problems, and the numerical results were compared with existing analytical and numerical solutions. The algorithm performance in solving practical soil–structure interaction problems was also investigated by simulating a large‐scale soft ground tunneling problem involving soil loss near an existing lining. Copyright © 2011 John Wiley & Sons, Ltd.     

Analysis of soil‐structural interface behavior using three‐dimensional DEM simulations

The shear behavior at the interface between the soil and a structure is investigated at the macroscale and particle‐scale levels using a 3‐dimensional discrete element method (DEM). The macroscopic mechanical properties and microscopic quantities affected by the normalized interface roughness and the loading parameters are analyzed. The macro‐response shows that the shear strength of the interface increases as the normalized roughness of the interface increases, and stress softening and dilatancy of the soil material are observed in the tests that feature rough interfaces. The particle‐scale analysis illustrates that a localized band characterized by intense shear deformation emerges from the contact plane and gradually expands as shearing progresses before stabilizing at the residual stress state. The thickness of the localized band is affected by the normalized roughness of the interface and the normal stress, which ranges between 4 and 5 times that of the median grain diameter. A thicker localized band is formed when the soil has a rough shearing interface. After the localized band appears, the granular material structuralizes into 2 regions: the interface zone and the upper zone. The mechanical behavior in the interface zone is representative of the interface according to the local average stress analysis. Certain microscopic quantities in the interface zone are analyzed, including the coordination number and the material fabric. Shear at the interface creates an anisotropic material fabric and leads to the rotation of the major principal stress.     

Dynamic soil–structure interaction analysis using Lanczos vectors with adaptive time integration technique

An analytical procedure to obtain the response of soil–structure interaction problems, time domain is described. The procedure makes use of large domain for descritization along with co-ordinate transformation using Lanczos vectors. The responses are obtained in time domain using an adaptive direct integration method. The scheme has the ability to estimate errors due to temporal discretization as well as co-ordinate transformation. The procedure has been applied to half-space problems and non-convex domains for validation of the scheme, and the scheme obeys causality condition in both the situations. The present method has all the advantages of time domain scheme which is local both in space and time with small computational effort. Copyright © 1999 John Wiley & Sons, Ltd.     

Cyclic macro‐element for soil–structure interaction: material and geometrical non‐linearities

This paper presents a non‐linear soil–structure interaction (SSI) macro‐element for shallow foundation on cohesive soil. The element describes the behaviour in the near field of the foundation under cyclic loading, reproducing the material non‐linearities of the soil under the foundation (yielding) as well as the geometrical non‐linearities (uplift) at the soil–structure interface. The overall behaviour in the soil and at the interface is reduced to its action on the foundation. The macro‐element consists of a non‐linear joint element, expressed in generalised variables, i.e. in forces applied to the foundation and in the corresponding displacements. Failure is described by the interaction diagram of the ultimate bearing capacity of the foundation under combined loads. Mechanisms of yielding and uplift are modelled through a global, coupled plasticity–uplift model. The cyclic model is dedicated to modelling the dynamic response of structures subjected to seismic action. Thus, it is especially suited to combined loading developed during this kind of motion. Comparisons of cyclic results obtained from the macro‐element and from a FE modelization are shown in order to demonstrate the relevance of the proposed model and its predictive ability. Copyright © 2001 John Wiley & Sons, Ltd.     

Discrete element method analysis of non‐coaxial flow under rotational shear

This study focuses on non‐coaxial flow behavior of cohesionless soil undergoing cyclic rotational shear, with a special interest in the effects of particle‐scale characteristics. To this end, we perform a series of 2D discrete element simulations with various particle shapes, inter‐particle coefficient of friction, initial density, and stress ratios. The validity and efficacy of the numerical model is established by systematically comparing numerical simulation results with existing laboratory testing results. Such comparison shows that the numerical simulations are capable of capturing mechanical behavior observed in laboratory testing under rotational shear. We further demonstrate and quantify a strong yet simple relationship between the deviatoric part of the normalized strain increment and the non‐coaxial angle, denoted by and ψ, respectively. This quantitative correlation between ψ and is independent of applied stress ratio, initial and current void ratio, and the number of cycles applied, but dependent on the principal stress orientation and particle‐scale characteristics. At the same , specimens with higher inter‐particle friction angle or smaller particle aspect ratio show greater non‐coaxial angles. A simple model is able to fit this ψ‐ relationship well, which provides a useful relationship that can be exploited in developing constitutive models for rotational shearing. Copyright © 2014 John Wiley & Sons, Ltd.     

A macro‐element for a circular foundation to simulate 3D soil–structure interaction

This paper presents a non‐linear interface element to compute soil–structure interaction (SSI) based on the macro‐element concept. The particularity of this approach lies in the fact that the foundation is supposed to be infinitely rigid and its movement is entirely described by a system of global variables (forces and displacements) defined in the foundation's centre. The non‐linear behaviour of the soil is reproduced using the classical theory of plasticity. Failure is described by the interaction diagram of the ultimate bearing capacity of the foundation under combined loads. The macro‐element is appropriate for modelling the cyclic or dynamic response of structures subjected to seismic action. More specifically, the element is able to simulate the behaviour of a circular rigid shallow foundation considering the plasticity of the soil under monotonic static or cyclic loading applied in three directions. It is implemented into FedeasLab, a finite element Matlab toolbox. Comparisons with experimental monotonic static and cyclic results show the good performance of the approach. Copyright © 2007 John Wiley & Sons, Ltd.     

地震分析中人工边界处理与地震动输入方法研究

基于柱面波波动方程,推导建立了适用于土-结构地震动力相互作用分析的地震动输入和人工边界的处理方法.其中,地震动的输入是通过在人工边界上施加等效节点力来实现的,等效节点力的大小与入射地震波波速成正比;而人工边界的处理方法使得人工边界条件不仅在时间上是局部的,而且在空间上也是局部的.这种处理方法简单、有效,物理意义清晰,且很容易在有限元法中实现,结合Newmark时间积分是无条件稳定的.为了验证方法的有效性和精度,给出了两个算例,分别用于检验人工边界条件的性能以及地震动输入方法的正确性.算例分析结果表明,所提出的方法是十分有效的.     

Granular element method for three‐dimensional discrete element calculations

This paper endows the recently‐proposed granular element method (GEM) with the ability to perform 3D discrete element calculations. By using non‐uniform rational B‐Splines to accurately represent complex grain geometries, we proposed an alternative approach to clustering‐based and polyhedra‐based discrete element methods whereby the need for complicated and ad hoc approaches to construct 3D grain geometries is entirely bypassed. We demonstrate the ability of GEM in capturing arbitrary‐shaped 3D grains with great ease, flexibility, and without excessive geometric information. Furthermore, the applicability of GEM is enhanced by its tight integration with existing non‐uniform rational B‐Splines modeling tools and ability to provide a seamless transition from binary images of real grain shapes (e.g., from 3D X‐ray CT) to modeling and discrete mechanics computations.Copyright © 2013 John Wiley & Sons, Ltd.     

Discrete element modeling of tool‐rock interaction II: rock indentation

The failure mechanisms induced by a wedge‐shaped tool indenting normally against a rock surface are investigated using the discrete element method (DEM). The main focus of this study is to explore the conditions controlling the transition from a ductile to a brittle mode of failure. The development of a damage zone and the initiation and propagation of a brittle fracture is well captured by the DEM simulations. The numerical results support the conjecture that initiation of brittle fractures is governed by a scaled flaw length Λ, a ratio between the flaw size λ and the characteristic length (where K_Ic is the toughness and σ_c the uniaxial compressive strength). The size of the damage zone agrees well with analytical predictions based on the cavity expansion model. The effects of a far‐field confining stress and the existence of a relief surface near the indenter are also examined.Copyright © 2012 John Wiley & Sons, Ltd.     

Discrete element modeling of tool‐rock interaction I: rock cutting

Tool‐rock interaction processes can be classified as indentation or cutting depending on the direction of motion of the tool with respect to the rock surface. The modes of failure induced in the rock by an indenting or a cutting tool can be ductile and/or brittle. The ductile mode is associated with the development of a damage zone, whereas the brittle mode involves the growth of macrocracks. This is the first part of a series of two papers concerned with an analysis of the cutting and the indentation processes based on using the discrete element method. In this paper, numerical simulations of the cutting process are conducted to reproduce the transition from a ductile to a brittle failure mode with increasing depth of cut, which is observed in experiments. The numerical results provide evidence that the critical depth of cut d _* controlling the failure mode transition is related to the characteristic length ? = (K_Ic ∕ σ_c)² with K_Ic denoting the material toughness and σ_c its unconfined compressive strength. The nature of frictional contact between the cutter face and the rock in the ductile failure mode is also examined. It is shown that the inclination of the total cutting force is controlled by a multi‐directional flow mechanism ahead of the cutter that is related to the formation of a wedge of failed material, intermittently adhering to the cutter. As a result, the inclination of the total cutting force varies with the rake angle of the cutter and cannot be considered an intrinsic measure of the interfacial friction between the cutter and the rock. Copyright © 2012 John Wiley & Sons, Ltd.     

土-地下结构体系地震反应的简化分析方法

基于Penzien提出的土-结构动力相互作用分析的集中质量模型,考虑等价土体的层间剪切刚度与阻尼效应,提出了土-地下结构动力相互作用体系地震反应分析的简化分析方法,选取具有不同地震动特性的Taft波、汶川地震什邡八角波和松潘波作为基岩水平向输入地震动,采用该简化方法和二维有限元法对土-地铁地下车站结构体系的地震加速度反应特性进行了对比分析,结果表明：简化方法计算的地铁地下车站结构峰值加速度反应大于二维有限元法计算的地铁地下车站结构峰值加速度反应,两者的差异与输入地震动特性有关,但其随地铁地下车站结构高度变化的总体趋势较为一致;随着输入地震动强度的增大,其差异程度也有所加大。该简化方法可合理反映土-地下结构体系的动力相互作用效应,可作为地下结构抗震设计分析的一种辅助方法。     

Three‐dimensional finite element analysis of the interaction between tunneling and pile foundations

This paper concerns analysis of the impact of construction of urban tunnels on adjacent pile foundations. It is carried out using an elastoplastic three‐dimensional finite element modelling. Numerical simulations are performed in two stages, which concern, respectively, the application of the pile axial loading and the construction of the tunnel in presence of the pile foundations. Analysis is carried out for both single piles and groups of piles. Results of numerical simulations show that tunneling induces significant internal forces in adjacent piles. The distribution of internal forces depends mainly on the position of the pile tip regarding the tunnel horizontal axis and the distance of the pile axis from the centre of the tunnel. Analysis of the interaction between tunneling and a group of piles reveals a positive group effect with a high reduction of the internal forces in rear piles. Copyright © 2002 John Wiley & Sons, Ltd.     

层状地基-非线性结构系统地震响应特性分析

根据从系统能量导出的地震动力方程,采用渐近解析法对层状地基 非线性结构系统的地震响应进行了分析,给出了共振区的振幅-频率以及相位-频率的解析公式,并对地基参数的不同取值,给出了一系列相应的计算结果。进而根据计算结果分析了地基参数对结构地震响应特性的影响,讨论了共振区畸变、振幅的分叉和突变现象以及“路径”效应等非线性动力特性。分析结果表明,同一类别地基的参数变化对结构的动力特性有显著的影响,按现行抗震设计方法对层状地基 非结构系统进行抗震分析时,可能产生较大的偏差,需考虑相互作用和非线性特性,以使设计参数更合理和优化。     

Discrete element method modeling of inherently anisotropic rocks under uniaxial compression loading

A new numerical approach is proposed in this study to model the mechanical behaviors of inherently anisotropic rocks in which the rock matrix is represented as bonded particle model, and the intrinsic anisotropy is imposed by replacing any parallel bonds dipping within a certain angle range with smooth‐joint contacts. A series of numerical models with β = 0°, 15°, 30°, 45°, 60°, 75°, and 90° are constructed and tested (β is defined as the angle between the normal of weak layers and the maximum principal stress direction). The effect of smooth‐joint parameters on the uniaxial compression strength and Young's modulus is investigated systematically. The simulation results reveal that the normal strength of smooth‐joint mainly affects the behaviors at high anisotropy angles (β > 45°), while the shear strength plays an important role at medium anisotropy angles (30°–75°). The normal stiffness controls the mechanical behaviors at low anisotropy angles. The angle range of parallel bonds being replaced plays an important role on defining the degree of anisotropy. Step‐by‐step procedures for the calibration of micro parameters are recommended. The numerical model is calibrated to reproduce the behaviors of different anisotropic rocks. Detailed analyses are conducted to investigate the brittle failure process by looking at stress‐strain behaviors, increment of micro cracks, initiation and propagation of fractures. Most of these responses agree well with previous experimental findings and can provide new insights into the micro mechanisms related to the anisotropic deformation and failure behaviors. The numerical approach is then applied to simulate the stress‐induced borehole breakouts in anisotropic rock formations at reduced scale. The effect of rock anisotropy and stress anisotropy can be captured. Copyright © 2015 John Wiley & Sons, Ltd.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Parallel finite element analysis of seismic soil structure interaction using a PC cluster

This paper describes an implementation of a highly scalable parallel computational facility with high speedup efficiency using relatively low-cost hardware, which consists of a cluster of desktop personal computers (PCs) connected via a 10-Gigabit Ethernet. Two-levels of parallelization were implemented. Communication between different PCs was achieved using message passing interface (MPI) protocol. Domain decomposition was automated and based on element numbering. Domain continuity was assured largely by re-numbering the elements using a “front squasher” code prior to decomposition. Within each PC, the shared memory parallelization was implemented using either the open-multiprocessing (OpenMP) or the MPI protocol. Analysis of three different problems with number of degrees-of-freedom ranging from about 129,000 to about 2,260,000 shows a speedup efficiency generally above 70%. Super-linear speedup was achieved in several of the cases examined in this study, with the hybrid MPI-OpenMP approach generally performing better compared to the pure MPI method for parallelization. The results demonstrate the feasibility of acquiring a parallel computing facility with relatively modest outlay that is within the reach of consulting or engineering offices.     

Discrete element modelling of material non‐coaxiality in simple shear flows

We investigate the quasi‐static simple shear flow of a two‐dimensional assembly of cohesionless particles using discrete element method (DEM) simulations. We focus on the unsteady flow regime where the solid would experience significant evolution of stresses, mobilised shear strength and dilation. We construct the DEM model using a discretised‐wall confined granular cell where the apparent boundary is allowed to dilate or contract synchronously with the confined solid. A rather uniform simple shear field is achieved across the whole assembly, which benefits rheological studies in generalising constitutive laws for continuum methods. We examine two aspects of the simple shear behaviour: macroscopic stress and strain rate evolution, particularly the non‐coaxiality between the principal directions of the two; and micromechanics such as evolution of fabric. For an initially anisotropic specimen sheared under constant normal pressure, the direction of principal stress rotates towards that of the principal strain rate, gradually reducing the degree of non‐coaxiality from about 45° to fluctuating around 0°. The rate in approaching coaxiality is slower in samples with larger initial porosity, stress ratio and mean stress. Generally, a faster rate in approaching coaxiality in simple shear is observed in a more dilatant sample, which often shows a larger degree of mobilised fabric anisotropy, suggesting the possible important role of instantaneous internal friction angle. The evolution of principal fabric direction resembles that of the principal stress direction. © 2013 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons, Ltd.