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.     

On the role of pre-existing discontinuities on the micromechanical behavior of confined rock samples: a numerical study

The deformation process and failure mechanism of rock mass with increased density of initial joints subjected to confined stress state are investigated in this study using discrete element method (DEM). A numerical model of standard size granite samples is developed and validated using experimental data for both intact and jointed rocks. The micro-parameters of the rock material are first determined, and the effects of the rock discontinuity on strength, deformability, stress–strain relationship, and failure modes are then investigated at the macro-scale level. Analyses are also performed to examine the tensile and shear crack distributions, fragmentation characteristics, particle kinematics, and energy dissipation to advance the current understanding of the deformation processes and failure mechanisms of jointed rock masses. The microscopic evolutions in the fabric and force anisotropy during loading and distributions of contact forces provide insights into the influence of increasing initial jointing on the macroscopic deformational behavior of the rock. The results show how the deceleration in the growth of fabric and contact force anisotropies develops and confirms that the increase in initial jointing and the associated changes in microstructure can restrain the development of anisotropy, thereby reducing significantly the strength of the rock samples.

     

On the role of sphericity of falling rock clusters—insights from experimental and numerical investigations

In this study, the dynamic behavior of rock clusters falling on rough slopes and impacting a vertical barrier is investigated experimentally and numerically using discrete element analysis. A specially designed laboratory setup that involves a flume of adjustable slope lined with a bumpy surface and equipped with an instrumented wall at the toe is used in the experimental investigation. The velocity profiles and impact forces were measured for three inclination angles using two different rock clusters. Three-dimensional discrete element analysis is then conducted to investigate the mechanical behavior of the rockfall and examine the role of sphericity of the rock cluster on the overall behavior of the system. This was achieved by explicitly simulating the complex shapes of the used rocks and the rough surface of the slope. The material coefficient of friction was measured using heap tests, and the results are compared with those obtained numerically using four different particle sphericities. Conclusions are made regarding the effect of slope inclination angle and the volume of the cluster on the impact forces exerted on rigid barriers. This study suggests that rock sphericity plays important roles on the dynamic behavior of the system and should be taken into consideration in simulating rockfall problems.     

On the role of joint roughness on the micromechanics of rock fracturing process: a numerical study

The deformation process and failure mechanism of rock mass with increased joint roughness subjected to unconfined compression are investigated in this study using discrete element method. A numerical model is developed using soft-bonded particle and validated to realistically replicate the mechanical response of the rock mass. The micro-parameters of the rock material are first determined, and the effects of the joint roughness on the macromechanical response and fracture growth mechanism are then investigated. Analyses are also performed to examine the tensile and shear crack distributions, acoustic emission (AE) characteristics, coordination number, and crack anisotropy to advance the current understanding of the role of joint roughness on the mechanical behavior and deformability of rock mass. The results show that strength and deformability of the jointed rocks are highly dependent on the joint orientation and roughness. Joint roughness is found to restrain the propagation and coalescence of microcracks and AE events from the interlocking of asperities. In addition, the spatial distribution of the contact forces allows for better understanding of the effect of joint inclination angle on the response of the investigated rock samples.

     

Investigating the Effects of Local Contact Loss on the Earth Pressure Distribution on Rigid Pipes

In this paper, the earth pressure distribution acting on a buried pipe with localized support loss is investigated experimentally and numerically in this study. A laboratory setup has been designed to facilitate the simulation of the local wall separation and to track the changes in earth pressure at selected locations along the pipe circumference. Validated by the experimental results, two-dimensional finite element analysis has been conducted to examine the role of soil-pipe interaction on the pressure distribution around the pipe before and after the contact loss is introduced. Experimental and numerical results revealed that the presence of a gap between the pipe wall and the surrounding backfill can lead to significant changes in contact pressure and bending moment in the pipe wall in the immediate vicinity of the gap. This study suggests that efforts to detect and repair areas experiencing support loss should be made before significant changes in pressure develop causing stress concentration in the pipe wall as it may lead to pipe damage.     

Experimental evaluation of the performance of earth levees deteriorated by wildlife activities

In this study, an earth levee model is constructed to investigate the impact of animal burrows on the integrity and performance of earthen structures. A series of centrifuge experiments are conducted on homogenous scaled-down 1H:1V levee models built from the natural Kasama soil. Both intact and deteriorated models were subject to a 35g acceleration level. Invasive animal intrusions were introduced in the form of horizontal array of idealized cylindrical burrows at the mid-height of the levee. The water level was gradually increased during the centrifuge flight, and the response of the levee was monitored throughout the test. Pore pressures were recorded using pressure transducers placed at preselected locations within the model. Surface displacements were measured using laser LVDTs and supplemented with three digital cameras for tracking the overall deformation pattern of the levee model. A summary of the test procedure and selected results is presented herewith. The observed deformation mechanism due to the presence of animal burrows is also described. As compared with the intact levee, the presence of burrows is found to alter the pattern of the water flow through the deteriorated levee structure—leading to a notable increase in the exit hydraulic gradient, internal erosion, and subsequently slope failure.     

Evaluating the performance of an explicit dynamic relaxation technique in analyzing non-linear geotechnical engineering problems

Explicit dynamic relaxation is an efficient tool that has been used to solve problems involving highly non-linear differential equations. The key feature of this method is the ability to use explicit dynamic algorithms in solving static problems. Few attempts have been made to date to apply this technique in conventional geotechnical engineering. In this study, an algorithm that incorporates the application of a stiffness dependent time step scheme is proposed. The algorithm has been successfully used to solve 2D and 3D non-linear geotechnical engineering problems. To calibrate the developed algorithm, numerical simulations have been conducted for a strip and square footings supported by Mohr–Coulomb material. Performance of four different types of brick elements used in collapse load calculation is examined in terms of convergence speed and accuracy. In addition, the role of employing adaptive time steps in reducing the number of iterations needed for convergence is also evaluated.