A 2.5-D coupled FE–BE model for the dynamic interaction between saturated soil and longitudinally invariant structures

A two-and-a-half-dimensional (2.5-D) coupled finite element–boundary element (FE–BE) model is presented to simulate the three-dimensional dynamic interaction between saturated soils and structures with longitudinally invariant geometries. A regularized 2.5-D boundary integral equation for saturated porous media is derived that avoids the evaluation of singular traction integrals. The 2.5-D coupled FE–BE model is established by using the continuity conditions on the soil–structure interface. The developed model is verified through comparison with an existing semi-analytical method. Two case studies of a tunnel embedded in a poroelastic half-space and the efficiency of a vibration isolating screen are presented.     

Method of deformation analysis for composite structures of soils and masonry stones

The coupled numerical manifold method (NMM) and discontinuous deformation analysis (DDA) are enhanced to simulate deformations of continuous soil and discontinuous masonry structures. An elasto-plastic NMM-DDA is formulated that incorporates elasto-plastic constitutive laws into incremental forms of the equation of motion. A node-based uniform strain element is applied to avoid volumetric locking, which often occurs in conventional NMM-DDA. The proposed method is applied to three fundamental boundary value problems: a beam bending problem, a bearing capacity problem of a footing, and a bearing capacity problem of a masonry structure. The method is verified through comparisons with conventional solutions.     

Application of the third medium method for frictionless contact problems in geomechanics

Traditional approaches in contact mechanics demand complicated search algorithms at the interface between the contacting bodies. Recently, a new contact method based on the concept of a third medium has been developed, which overcomes the drawbacks of conventional contact mechanics techniques. This new scheme is based on a space filling mesh, in which the contacting bodies can move and interact. The ability and accuracy of this method in predicting displacements, as well as the contact forces, is validated by solving selected numerical examples. The potential merits of this method for analysing geotechnical problems by the finite element method are addressed.     

Empirical assessment of the critical time increment in explicit particulate discrete element method simulations

This contribution considers the critical time increment (Δt_crit) to achieve stable simulations using particulate discrete element method (DEM) codes that adopt a Verlet-type time integration scheme. The Δt_crit is determined by considering the maximum vibration frequency of the system. Based on a series of parametric studies, Δt_crit is shown to depend on the particle mass (m), the maximum contact stiffness (K_max), and the maximum particle coordination number (C_N,max). Empirical expressions relating Δt_crit to m, K_max, and C_N,max are presented; while strictly only valid within the range of simulation scenarios considered here, these can inform DEM analysts selecting appropriate Δt_crit values.     

A comparison of discrete element and micromechanical methods for determining the effective elastic properties of geomaterials

Data generated from scanning electron microscopy images of oil-shale geomaterials are used to generate configurations and to acquire parameters required for use in homogenization schemes for the determination of the effective elastic properties of the samples. Two alternative homogenization methods are employed: numerical simulation using the Discrete Element Method and the Polycrystalline Self-Consistent Method from micromechanics. The schemes give rise to predictions of the effective elastic properties that are in very good agreement.     

Assessing the impacts of a large slope failure using 3DEC: The Chiu-fen-erh-shan residual slope

This paper assesses the impact of Chiu-fen-erh-shan residual slope failure using a three-dimensional distinct element program. The simulation results indicate that rocks will severely damage the Lung-Nan path. The Taanshan syncline generates a depression zone adjoining to the slope toe and acts as a buffer zone to trap failing rocks. Some rocks will travel through the syncline, which poses a risk to visitors at the monument and the Shihmen observation deck. Few rocks will pass though the memorial park near the northern boundary of the slope. Visitors should be evacuated out of the impact area during a heavy rainfall event.     

Dimensionless input parameters in discrete element modeling and assessment of scaling techniques

A set of dimensionless input parameters were defined for DEM using a characteristic time which is a function of density and elastic modulus of particles and an arbitrary characteristic length. Dimensionless strain rate and mass damping ratio are inversely proportional to the characteristic time, and stress is normalized by elastic modulus to give dimensionless stress. It was demonstrated that the response of a model in the dimensionless scale is invariant with the choice of density, elastic modulus and the characteristic length if dimensionless strain rate and mass damping ratio are kept constant. Small time step is a prohibitive aspect of DEM. Scaling techniques are widely employed to enlarge the time step. Using the dimensionless scheme, it was learned that density scaling is equivalent to the use of a higher strain rate, and stiffness scaling results in a higher strain rate and an elevated stress state in the dimensionless scale.     

Predicting bending failure of CDM columns under embankment loading

The mechanism of bending failure and the magnitude of the bending moment in a column formed by cement deep mixing (CDM) under an embankment load were investigated by a series of three-dimensional (3D) finite element analyses. Based on the numerical results, a design method to consider the bending failure of the CDM column has been established. The usefulness of the proposed methods was verified using the results of two centrifuge model tests of the embankments on the CDM column improved soft model ground model reported in the literature, and the columns had tensile cracks or showed complete failure.     

Reliability bearing capacity analysis of footings on cohesive soil slopes using RFEM

Reliability analysis of bearing capacity of a strip footing at the crest of a simple slope with cohesive soil was carried out using the random finite element method (RFEM). Analyses showed that the coefficient of variation and the spatial correlation length of soil cohesion can have a large influence on footing bearing capacity, particularly for slopes with large height to footing width ratios. The paper demonstrates cases where a footing satisfies a deterministic design factor of safety of 3 but the probability of design failure is unacceptably high. Isotropic and anisotropic spatial variability of the soil strength was also considered.     

Seepage simulation using pipe network flow model in a discrete element system

The pipe network flow model can simulate the seepage process with DEM conveniently because of its simple algorithm. However, whether it can recover the correct seepage process has not been verified. In this paper, the equation to update the fluid pressure is rebuilt according to the flow conservation. Through the steady seepage simulation, this algorithm is verified to be able to recover Darcy’s law, and the equation to calibrate the aperture according to macro permeability is derived. Furthermore, the modified algorithm is used to simulate the unsteady seepage process, and the results show good agreement with the analytical solutions.