A new contact algorithm in the material point method for geotechnical simulations

Contact between stiff structural elements and soil is encountered in many applications in geotechnical engineering. Modelling of such contact is challenging as it often involves impact that would lead to large deformation and failure of the soil. The Material Point Method (MPM) is a mesh‐free method that has been applied to simulate such phenomena. However, the frictional contact algorithm commonly used in MPM only supports Coulomb friction and cannot model fully or partially rough contact conditions in terms of geotechnical engineering. Moreover, because of very different stiffness of contacting materials, the contact force predicted by the previous frictional contact algorithms usually suffers from severe oscillation when applied in structure–soil interaction. This paper presents a new contact algorithm, termed Geo‐contact, designed for geotechnical engineering. In Geo‐contact, a penalty function is incorporated to reduce the oscillation in contact computation, and a limited shear stress is specified along the contact interface. The proposed Geo‐contact algorithm has been implemented to simulate smooth, partially rough and rough contact in typical large deformation penetration problems. The resistance–displacement curves obtained using the Geo‐contact are compared with analytical solutions of limit analysis and large deformation finite element results to verify the accuracy and robustness of the proposed contact algorithm. Copyright © 2014 John Wiley & Sons, Ltd.     

Biaxial test simulations using a packing of polygonal particles

The mechanical response of cohesionless granular materials under monotonic loading is studied by performing molecular dynamic simulations. The diversity of shapes of soil grains is modelled by using randomly generated convex polygons as granular particles. Results of the biaxial test obtained for dense and loose media show that samples achieve the same void ratio at large strains independent of their initial density state. This limit state resembles the so‐called critical state of soil mechanics, except for some stress fluctuations, which remain for large deformations. These fluctuations are studied at the micro‐mechanical level, by following the evolution of the co‐ordination number, force chains and the fraction of the sliding contacts of the sample. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

Soil‐pile interaction in the pile vertical vibration considering true three‐dimensional wave effect of soil

The dynamic response of an end bearing pile embedded in a linear visco‐elastic soil layer with hysteretic type damping is theoretically investigated when the pile is subjected to a time‐harmonic vertical loading at the pile top. The soil is modeled as a three‐dimensional axisymmetric continuum in which both its radial and vertical displacements are taken into account. The pile is assumed to be vertical, elastic and of uniform circular cross section. By using two potential functions to decompose the displacements of the soil layer and utilizing the separation of variables technique, the dynamic equilibrium equation is uncoupled and solved. At the interface of soil‐pile system, the boundary conditions of displacement continuity and force equilibrium are invoked to derive a closed‐form solution of the vertical dynamic response of the pile in frequency domain. The corresponding inverted solutions in time domain for the velocity response of a pile subjected to a semi‐sine excitation force applied at the pile top are obtained by means of inverse Fourier transform and the convolution theorem. A comparison with two other simplified solutions has been performed to verify the more rigorous solutions presented in this paper. Using the developed solutions, a parametric study has also been conducted to investigate the influence of the major parameters of the soil‐pile system on the vertical vibration characteristics of the pile. Copyright © 2013 John Wiley & Sons, Ltd.     

Dynamic characteristics of a large‐diameter pile in saturated soil and its application

This study theoretically investigates the dynamic response of an end‐bearing pile embedded in saturated soil considering the transverse inertial effect of the pile. The saturated soil surrounding the pile is described by Biot poroelastic theory, and the pile is represented by a Rayleigh‐Love rod because both the vertical and radial displacements at the soil‐pile interface are considered. The potential function decomposition method and variable separation method are introduced to solve the governing equations of the soil, in which the vertical and radial displacement components are coupled. The governing equation of the pile is solved using the continuity conditions at the pile‐soil interface. Next, the velocity admittance in the frequency domain and the velocity response in the time domain at the pile top are presented based on the Laplace transform and inverse Fourier transform, respectively. Subsequently, the reduced solution is compared with a 1‐dimensional model solution to verify the validity, and the influences of the slenderness ratio of the pile on the transverse inertial effect of the pile are analyzed. Moreover, Poisson ratio, the slenderness ratio of the pile, and the pile‐soil modulus ratio are studied. Finally, the theoretical and measured curves in the engineering project are compared, and the results demonstrate the good application prospects of the solution presented in this article.     

Study on low‐strain integrity testing of pipe‐pile using the elastodynamic finite integration technique

This study focuses on the three‐dimensional (3‐D) characteristics of wave propagation in pipe‐pile using elastodynamic finite integration technique. First, a real 3‐D pile‐soil model in cylindrical coordinate system is presented. Then, the governing equations are established. With the boundary and initial conditions, the numerical solution is obtained. The accuracy and feasibility of the self‐written code are further verified via comparing with the measured data. Velocity histories at different angles of pile top and pile tip are illustrated, and the snapshots reflecting the 3‐D characteristics of wave propagation are given. It shows that the interferences of Rayleigh waves can confuse the result interpretation for pile integrity testing. The increase of hammer contact time can effectively mitigate the interferences, and the interferences of Rayleigh waves are weakest at an angle of 90° from where hammer hits. Besides, surrounding soil can partly mitigate the wave interferences. Copyright © 2012 John Wiley & Sons, Ltd.     

Inverse Analysis of Deep Excavation Using Differential Evolution Algorithm

This paper presents the applications of the differential evolution (DE) algorithm in back analysis of soil parameters for deep excavation problems. A computer code, named Python‐based DE, is developed and incorporated into the commercial finite element software ABAQUS, with a parallel computing technique to run an FE analysis for all trail vectors of one generation in DE in multiple cores of a cluster, which dramatically reduces the computational time. A synthetic case and a well‐instrumented real case, that is, the Taipei National Enterprise Center (TNEC) project, are used to demonstrate the capability of the proposed back‐analysis procedure. Results show that multiple soil parameters are well identified by back analysis using a DE optimization algorithm for highly nonlinear problems. For the synthetic excavation case, the back‐analyzed parameters are basically identical to the input parameters that are used to generate synthetic response of wall deflection. For the TNEC case with a total of nine parameters to be back analyzed, the relative errors of wall deflection for the last three stages are 2.2, 1.1, and 1.0%, respectively. Robustness of the back‐estimated parameters is further illustrated by a forward prediction. The wall deflection in the subsequent stages can be satisfactorily predicted using the back‐analyzed soil parameters at early stages. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical simulation of the response of a reinforced wall to a high speed train passage

The dynamic response of a mechanically stabilized earth wall to the passing of a high‐speed train is modelled using the finite element method. A three‐dimensional analysis is carried out, using a specific framework that allows performing the analysis with a moderate computational effort. In the first place, a so‐called multiphase approach is used to take into account the reinforcing strips. The moving load is taken into account by performing the calculation in a mobile referential using the properties of symmetry of the train cars and a simplifying assumption of periodicity for the whole train. We also assume a steady state. A partial validation of the approach is obtained by means of a comparison with an analytical solution. The quick increase in displacements induced by the train passing when the speed comes close to the celerity of Rayleigh waves clearly appears in the results. The vertical displacements, vertical stresses in the backfill, tensile forces in the strips and the influence of the stiffness of the soil are discussed. Copyright © 2017 John Wiley & Sons, Ltd.     

Using TDR and modeling tools to investigate effects of interactive factors on preferential flow and transport in field sandy clay soil

Investigation of the extent of preferential flow and transport affected by several factors and their interactions in the vadose zone using the advanced measurement and modeling techniques is crucial for protection of groundwater from agricultural chemicals like fertilizers and pesticides. The objectives of this study were to investigate the interactive effects of factors like soil structure, initial soil water content (SWC), and application rate on preferential flow and transport using the time domain reflectometry (TDR) measurements of SWC and electrical conductivity (EC) in the plots of 12 treatments in a sandy clay field soil, applying the models (HYDRUS-1D and MACRO) to the measured data, performing Tukey test statistical analysis, and relating model parameters to basic soil properties, consequently, flow and transport characteristics. Analysis of response times and the changes in SWC and EC with time during the experiments in the profiles of the treatments confirmed the existence of preferential flow and transport at the site. As long as the other factors or conditions were constant; undisturbed versus disturbed soil, wet versus dry initial SWC, and high versus low application rate caused preferential flow and transport in the soil based on the TDR measurements. Overall HYDRUS-1D had better performance than MACRO in the simulations of the measured data. These results suggest that different scales of these factors in different field soils need to be further studied for better understanding the flow and transport processes in the vadose zone.     

Equivalent-Linear Seismic Ground Response Analysis of Some Typical Sites in Mumbai

Any earthquake event is associated with a rupture mechanism at the source, propagation of seismic waves through underlying rock and finally these waves travel through the soil layers to the particular site of interest. The bedrock motion is significantly modified at the ground surface due to the presence of local soil layers above the bedrock beneath the site of interest. The estimation of the amplifications in ground response due to the local soil sites is a complex problem to the designers and the problem is more important for mega cities like Mumbai in India, where huge population may get affected due to devastations of earthquake. In the present study, the effect of local soil sites in modifying ground response is studied by performing one dimensional equivalent-linear ground response analysis for some of the typical Mumbai soil sites. Field borelog data of some typical sites in Mumbai city viz. Mangalwadi site, Walkeswar site, BJ Marg near Pandhari Chawl site are considered in this study. The ground responses are observed for range of input motions and the results are presented in terms of surface acceleration time history, ratio of shear stress to vertical effective stress versus time, acceleration response spectrum, Fourier amplitude ratio versus frequency etc. The typical amplifications of ground accelerations considering four strong ground motions with wide variation of low to high MHA, frequency contents and durations are obtained. Results show that MHA, bracketed duration, frequency content have significant effects on the amplification of seismic accelerations for typical 2001 Bhuj motion. The peak ground acceleration amplification factors are found to be about 2.50 for Mangalwadi site, 2.60 for Walkeswar site and 3.45 for BJ Marg site using 2001 Bhuj input motion. The response spectrum along various soil layers are obtained which will be useful for designers for earthquake resistant design of geotechnical structures in Mumbai for similar sites in the absence of site specific data.     

A multidirectional semi‐analytical method for analysis of laterally loaded pile groups in multi‐layered elastic strata

A semi‐analytical method for calculating the response of single piles and pile groups subjected to lateral loading is developed in this paper. Displacements anywhere in the soil domain are tied to the displacements of the piles through decay functions. The principle of virtual work and the calculus of variations are used to derive the governing differential equations that describe the response of the piles and soil. The eigenvalue method and the finite difference technique are used to solve the system of coupled differential equations for the piles and soil, respectively. The proposed method takes into account the soil surface displacement along and perpendicular to the loading direction and produces displacement fields that are very close to those produced by the finite element method but at lower computational effort. Compared with the previous method that considered only the soil displacement along the loading direction, accounting for the multi‐directional soil displacement field produces responses for the piles and soil that are closer to those approximated by the finite element method. Copyright © 2016 John Wiley & Sons, Ltd.     

岩体抗剪强度参数现场测试新方法及工程应用

岩体抗剪强度参数的确定是进行岩体工程稳定性分析中最基本和最重要的,也是最困难的一项工作。由于室内和现场试验方法的局限性以及经验准则的人为性,快速、简易和较可靠确定岩体抗剪强度参数的原位测试手段一直是岩石工程界关注的热点。采用引进美国生产的岩石钻孔剪切仪,在实验室混凝土空心圆柱体和向家坝水电站坝基岩体中进行了尝试,取得了有益的试验成果。将主要从钻孔剪切仪的工作原理、数值模拟分析和工程应用等几方面对其进行全面的阐述和总结。试验成果和数值分析结果显示,岩石钻孔剪切仪是一种适用于软岩~中硬岩的原位试验设备,其试验成果总体钻孔剪切试验得到的c、 小于传统三轴和直剪试验成果,但钻孔剪切试验的 值在 <40?的范围内略低于传统三轴和直剪试验成,而钻孔剪切试验的c值较传统三轴和直剪试验的c 值降低较多,一般为25%～50%。     

A coupling method for soil structure interaction problems

This paper describes a soil‐structure coupling method to simulate blast loading in soil and structure response. For the last decade, simulation of soil behavior under blast loading and its interaction with semi buried structure in soil becomes the focus of computational engineering in civil and mechanical engineering communities. In current design practice, soil‐structure interaction analysis often assumes linear elastic properties of the soil and uses small displacement theory. However, there are numerous problems, which require a more advanced approach that account for soil‐structure interaction and appropriate constitutive models for soil. In simplified approaches, the effect of soil on structure is considered using spring‐dashpot‐mass system, and the blast loading is modeled using linearly decaying pressure–time history based on equivalent trinitrotoluene and standoff distance, using ConWep, a computer program based on semi‐empirical equations. This strategy is very efficient from a CPU time computing point of view but may not provide accurate results for the dynamic response of the structure, because of its significant limitations, mainly when soil behavior is strongly nonlinear and when the buried charge is close to the structure. In this paper, both soil and explosive are modeled using solid elements with a constitutive material law for soil, and a Jones–Wilkins–Lee equation of state for explosive. One of the problems we have encountered when solving fluid structure interaction problems is the high mesh distortion at the contact interface because of high fluid nodal displacements and velocities. Similar problems have been encountered in soil structure interaction problems. To prevent high mesh distortion for soil, a new coupling algorithm is performed at the soil structure interface for structure loading. The coupling method is commonly used for fluid structure interaction problems in automotive and aerospace industry for fuel sloshing tank, and bird impact problems, but rarely used for soil structure interaction problems, where Lagrangian contact type algorithms are still dominant. Copyright © 2012 John Wiley & Sons, Ltd.     

Investigation of wetting front behavior due to rainfall and ponding head effects in arid region wadi alluvium

Intensive experimental as well as numerical investigations have been carried out in order to develop empirical relationships to predict wetting front travel time and infiltration rate in deep arid region alluvium due to the effect of variable rainfall intensity and ponding head, respectively. A rainfall simulator that is capable of producing variable rainfall intensity up to 160 mm/h and a large diameter infiltrometer were constructed in the study area and then used to test the effect of rainfall and ponding head up to a depth of 3 m. Results of these experiments were utilized to calibrate an unsaturated zone model which was then used to test different scenarios for wetting front movement at depths beyond the scope of these experiments. Empirical infiltration relationships for rainfall and ponding head, respectively, were derived for shallow and deep soil columns from these experiments.     

降雨入渗下非饱和坡残积土湿润锋运移试验研究

我国东南沿海地区丘陵山地发育，降雨入渗到土体中的水分是导致滑坡灾害频发的关键因素。以福建省泉州市德化县地质灾害点为主要研究对象，考察典型地质灾害点具有代表性的坡积土与花岗岩残积土的渗透特性。利用自行研制的土体入渗装置，分别在降雨强度为15、30、60 mm/h条件下，考虑降雨历时一致（180 min）与过程降雨量一致（90 mm）两种工况开展一维土柱入渗试验，得到相应的各个土柱含水率、湿润锋、入渗率随时间变化的响应规律。试验结果表明：（1）土体渗透系数越大、雨强越大，土体湿润蔓延距离越深、速度越快。（2）降雨入渗过程中，土体含水率由浅及深逐次对降雨进行响应。不同雨强对含水率的影响主要体现在第一次响应时间以及饱和速度上，雨强越大，响应时间越快，饱和速度也越快。（3）提出可表征在不同雨强作用下，德化县马坪滑坡与崩土岭滑坡的湿润锋入渗公式。该研究成果对台风暴雨型滑坡的孕灾机制分析以及精细化监测预警具有重要的理论及实际意义。     

Vertical dynamic response of a rigid foundation embedded in a poroelastic soil layer

A simplified analytical method is presented for the vertical dynamic analysis of a rigid, massive, cylindrical foundation embedded in a poroelastic soil layer. The foundation is subjected to a time‐harmonic vertical loading and is perfectly bonded to the surrounding soil in the vertical direction. The soil underlying the foundation base is represented by a single‐layered poroelastic soil based on rigid bedrock while the soil at the side of the foundation is modeled as an independent poroelastic layer composed of a series of infinitesimally thin layers. The behavior of the soil is governed by Biot's poroelastodynamic theory and its governing equations are solved by the use of Hankel integral transform. The contact surface between the foundation base and the soil is smooth and fully permeable. The dynamic interaction problem is solved following standard numerical procedures. The accuracy of the present solution is verified by comparisons with the well‐known solutions obtained from other approaches for both the elastodynamic interaction problem and poroelastodynamic interaction problem. Numerical results for the vertical dynamic impedance and response factor of the foundation are presented to demonstrate the influence of nondimensional frequency of excitation, soil layer thickness, poroelastic material parameters, depth ratio and mass ratio on the dynamic response of a rigid foundation embedded in a poroelastic soil layer. Copyright © 2008 John Wiley & Sons, Ltd.     

A semi-analytical solution to consolidation of unsaturated soils with the free drainage well

Based on Fredlund’s one-dimensional consolidation equation for unsaturated soil, Darcy’s law and Fick’s law, a semi-analytical solution was presented to the free drainage well with a finite thickness under application of uniform vertical loading and the boundary of the top and bottom surfaces impermeable to water and air. According to the polar governing equations of water and air phases and the boundary and initial conditions, the excess pore-air and pore-water pressures and the soil layer settlement in the Laplace transformed domain are obtained by performing the Laplace transform and utilizing the Bessel functions. Crump’s method is used to perform the inversion of Laplace transform in order to obtain numerical solutions in the real time domain. Finally, a typical example is given to illustrate the changes in the excess pore-air and pore-water pressures and soil layer settlement with time factor at different ratios of air–water permeability coefficient and/or different distances from the well.     

Direct back analysis by the meshless local Petrov‐Galerkin method and Bayesian statistics

This study presents a preliminary development of a direct back analysis procedure by the meshless local Petrov–Galerkin (MLPG) method and Bayesian statistics and the application of resulting procedure to characterize soil properties using laboratory results. As compared to direct back analyses that are based on the finite element method (FEM), it is intended to show that MLPG‐based direct back analyses may be more suitable for some kinds of characterization problems; for example, involving a complex subsurface stratification or the characterization of soil properties of just an inclusion of a soil profile. The existing MLPG method is first slightly modified to analyse time‐dependent problems. Using the resulting method, quantities to be characterized are evaluated so that they give numerical results as close to measured data as possible. The Akaike information criterion is introduced for simplifying the evaluation. A one‐dimensional finite strain consolidation problem is introduced to do an error analysis for prediction by the proposed MLPG method. Another example illustrates experiences of performing an MLPG‐based direct back analysis. Comparison of MLPG‐based and FEM‐based direct back analyses is taken. Copyright © 2006 John Wiley & Sons, Ltd.     

Effects of boundary conditions and partial drainage on cyclic simple shear test results—a numerical study

Owing to imperfect boundary conditions in laboratory soil tests and the possibility of water diffusion inside the soil specimen in undrained tests, the assumption of uniform stress/strain over the sample is not valid. This study presents a qualitative assessment of the effects of non‐uniformities in stresses and strains, as well as effects of water diffusion within the soil sample on the global results of undrained cyclic simple shear tests. The possible implications of those phenomena on the results of liquefaction strength assessment are also discussed. A state‐of‐the‐art finite element code for transient analysis of multi‐phase systems is used to compare results of the so‐called ‘element tests’ (numerical constitutive experiments assuming uniform stress/strain/pore pressure distribution throughout the sample) with results of actual simulations of undrained cyclic simple shear tests using a finite element mesh and realistic boundary conditions. The finite element simulations are performed under various conditions, covering the entire range of practical situations: (1) perfectly drained soil specimen with constant volume, (2) perfectly undrained specimen, and (3) undrained test with possibility of water diffusion within the sample. The results presented here are restricted to strain‐driven tests performed for a loose uniform fine sand with relative density D_r=40%. Effects of system compliance in undrained laboratory simple shear tests are not investigated here. Copyright © 2004 John Wiley & Sons, Ltd.     

Radial Permeability Measurements for Shale Using Variable Pressure Gradients

Shale gas is becoming an important component of the global energy supply, with permeability a critical controlling factor for long-term gas production. Obvious deviation may exist between helium permeability determined using small pressure gradient(SPG) methods and methane permeability obtained under actual field production with variable pressure gradients(VPG). In order to more accurately evaluate the matrix permeability of shale, a VPG method using real gas(rather than He) is established to render permeability measurements that are more representative of reservoir conditions and hence response. Dynamic methane production experiments were performed to measure permeability using the annular space in the shale cores. For each production stage, boundary pressure is maintained at a constant and the gas production with time is measured on the basis of volume change history in the measuring pump. A mathematical model explicitly accommodating gas desorption uses pseudo-pressure and pseudo-time to accommodate the effects of variations in pressuredependent PVT parameters. Analytical and semi-analytical solutions to the model are obtained and discussed. These provide a convenient approach to estimate radial permeability in the core by nonlinear fitting to match the semi-analytical solution with the recorded gas production data. Results indicate that the radial permeability of the shale determined using methane is in the range of 1×10-6– 1×10-5 mD and decreases with a decrease in average pore pressure. This is contrary to the observed change in permeability estimated using helium. Bedding geometry has a significant influence on shale permeability with permeability in parallel bedding orientation larger than that in perpendicular bedding orientation. The superiority of the VPG method is confirmed by comparing permeability test results obtained from both VPG and SPG methods. Although several assumptions are used, the results obtained from the VPG method with reservoir gas are much closer to reality and may be directly used for actual gas production evaluation and prediction, through accommodating realistic pressure dependent impacts.