Effects of periodic fluctuations of micro‐polar boundary conditions on shear localizations in granular soil–structure interaction

In this work, the interface behavior between an infinite extended narrow granular layer and a rough surface of rigid body is investigated numerically, using finite element method in the updated Lagrangian (UL) frame. In this regard, the elasto‐plastic micro‐polar (Cosserat) continuum approach is employed to remove the limitations caused by strain‐softening of materials in the classical continuum. The mechanical properties of cohesionless granular soil are described with Lade's model enhanced by polar terms, including Cosserat rotations, curvatures, and couple stresses. Furthermore, the mean grain diameter as the internal length is incorporated into the constitutive relations accordingly. Here, the evolution and location of shear band, within the granular layer in contact with the rigid body, are mainly focused. In this regard, particular attention is paid to the effects of homogeneous distribution and periodic fluctuation of micro‐polar boundary conditions, prescribed along the interface. Correspondingly, the effects of pressure level, mean grain diameter, and stratified soil are also considered. The finite element results demonstrate that the location and evolution of shear band in the granular soil layer are strongly affected by the non‐uniform micro‐polar boundary conditions, prescribed along the interface. It is found that the shear band is located closer to the boundary with less restriction of grain rotations. Furthermore, the predicted thickness of shear band is larger for higher rotation resistance of soil grains along the interface, larger mean grain diameter, and higher vertical pressure. Regarding the stratified soil, comprising a thin layer with slightly different initial void ratio, the shear band moves towards the layer with initially higher void ratio. Copyright © 2011 John Wiley & Sons, Ltd.     

A general approach to the modelling of contaminant transport through composite landfill liners with intact or leaking geomembranes

Finite‐element models of contaminant transport through composite landfill liners require highly refined meshes around the interface between the geomembrane and the clay layer, especially if leakage through holes in the geomembrane is considered. In addition, no general formulation for transport through leaking geomembranes can be found in the literature. The paper develops a general approach to time‐dependent contaminant migration through composite liners with intact or leaking geomembranes. Equations are derived for various combinations of system conditions including Dirichlet and Neumann boundary conditions in the waste, constant mass of contaminants in the waste, steady state or transient transport in the geomembrane, and steady state or transient seepage velocities in the mineral liner. The effect of the geomembrane on transport in the soil is converted into an equivalent boundary condition applicable at the top of the clay layer. Hence, only the media underlying the top geomembrane are explicitly represented in the numerical model, yielding a computationally efficient algorithm. The new formulation is validated in conjunction with finite‐layer, finite‐element and boundary‐element methods, by comparing its predictions to those of more conventional approaches which represent the geomembrane explicitly. The scope of the method is illustrated by modelling a landfill liner with a geomembrane leaking in five locations. Copyright © 2007 John Wiley & Sons, Ltd.     

Exact solutions for one‐dimensional consolidation of single‐layer unsaturated soil

Based on the Fredlund consolidation theory of unsaturated soil, exact solutions of the governing equations for one‐dimensional consolidation of single‐layer unsaturated soil are presented, in which the water permeability and air transmission are assumed to be constants. The general solution of two coupled homogeneous governing equations is first obtained. This general solution is expressed in terms of two functions psi₁ and ψ₂, where ψ₁ and ψ₂, respectively, satisfy two second‐order partial differential equations, which are in the same form. Using the method of separation of variables, the two partial differential equations are solved and exact solutions for three typical homogeneous boundary conditions are obtained. To obtain exact solutions of nonhomogeneous governing equations with three typical nonhomogeneous boundary conditions, the nonhomogeneous boundary conditions are first transformed into homogeneous boundary conditions. Then according to the method of undetermined coefficients and exact solutions of homogenous governing equations, the series form exact solutions are put forward. The validity of the proposed exact solutions is verified against other analytical solutions in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

Semi‐analytical far field model for three‐dimensional finite‐element analysis

A challenging computational problem arises when a discrete structure (e.g. foundation) interacts with an unbounded medium (e.g. deep soil deposit), particularly if general loading conditions and non‐linear material behaviour is assumed. In this paper, a novel method for dealing with such a problem is formulated by combining conventional three‐dimensional finite‐elements with the recently developed scaled boundary finite‐element method. The scaled boundary finite‐element method is a semi‐analytical technique based on finite‐elements that obtains a symmetric stiffness matrix with respect to degrees of freedom on a discretized boundary. The method is particularly well suited to modelling unbounded domains as analytical solutions are found in a radial co‐ordinate direction, but, unlike the boundary‐element method, no complex fundamental solution is required. A technique for coupling the stiffness matrix of bounded three‐dimensional finite‐element domain with the stiffness matrix of the unbounded scaled boundary finite‐element domain, which uses a Fourier series to model the variation of displacement in the circumferential direction of the cylindrical co‐ordinate system, is described. The accuracy and computational efficiency of the new formulation is demonstrated through the linear elastic analysis of rigid circular and square footings. Copyright © 2004 John Wiley & Sons, Ltd.     

煤层地下开采地表沉陷预测的边值方法

采空区上覆岩、土体对地表沉陷的耦合行为可归结为弹性理论的边值问题,由不同性状松散表土层确定的应力边界和岩、土体之间形成的位移边界构成了相应边值问题的定解条件。在已有研究成果的基础上,提出了地表下沉预测研究的边值提法,以探讨地下开采引起的地表沉陷规律。给出了采空区上覆表土层内任一点的位移和地表下沉的计算公式,确定了地表沉陷的范围,并以神东矿区大柳塔1203工作面开采引发的地表沉陷为例给出了工程实例。研究结果表明,采空区上覆岩体对地表及表土层内任一点下沉的影响取决于上部岩层破断后的结构形态;采空区上覆土体对地表下沉量的作用与土体性质及表土层厚度有关,土体对表土层下沉的影响由下而上增大,呈非线性变化,至地表达最大值;地表沉陷区域随表土层厚度的增加而增大,其影响范围远大于采空区。     

Analytical layer‐element method for non‐axisymmetric consolidation of multilayered soils

A numerically efficient and stable method is developed to analyze Biot's consolidation of multilayered soils subjected to non‐axisymmetric loading in arbitrary depth. By the application of a Laplace–Hankel transform and a Fourier expansion, the governing equations are solved analytically. Then, the analytical layer‐element (i.e. a symmetric stiffness matrix) describing the relationship between generalized displacements and stresses of a layer is exactly derived in the transformed domain. Considering the continuity conditions between adjacent layers, the global stiffness matrix of multilayered soils is obtained by assembling the inter‐related layer‐elements. Once the solution in the Laplace–Hankel transformed domain that satisfies the boundary conditions has been obtained, the actual solution can be derived by the inversion of the Laplace–Hankel transform. Finally, numerical examples are presented to verify the theory and to study the influence of the layered soil properties and time history on the consolidation behavior. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical solution for consolidation and desiccation of soft soils

The consolidation and desiccation behaviour of soft soils can be described by two time‐dependent non‐linear partial differential equations using the finite strain theory. Analytical solutions do not exist for these governing equations. In this paper, we develop efficient numerical methods and software for finding the numerical solutions. We introduce a semi‐implicit time integration scheme, and show numerically that our method converges. In addition, the numerical solution matches well with the experimental result. A boundary refinement method is also developed to improve the convergence and stability for the case of Neumann type boundary conditions. Interface governing equations are derived to maintain the continuity of consolidation and desiccation processes. This is useful because the soil column can undergo desiccation on top and consolidation on the bottom simultaneously. The numerical algorithms has been implemented into a computer program and the results have been verified with centrifuge test results conducted in our laboratory. Copyright © 2001 John Wiley & Sons, Ltd.     

基于传递矩阵法的层状土中管桩水平动力阻抗分析

在考虑土体分层特性的基础上,分别建立了管桩桩周土体和桩芯土体的水平振动控制方程。通过引入势函数并考虑桩周土和桩芯土径向位移和环向位移的边界条件及其奇偶性,求得了管桩-土动力相互作用的刚度系数和阻尼系数。将土体模拟为连续分布的弹簧-阻尼器,并考虑桩芯土和桩周土的作用,建立了层状土中管桩的水平振动方程。借助初参数法和传递矩阵法求解了管桩的水平振动,得到了管桩桩顶的水平动力阻抗。通过数值分析,得到了土层剪切模量、管桩壁厚、桩周土和桩芯土剪切模量比、土层厚度等对管桩桩顶动力阻抗的影响规律。土层剪切模量、管桩壁厚、桩周土和桩芯土剪切模量比对层状土中管桩水平振动的影响主要在低频处,土层厚度在较宽的频率范围内对管桩水平振动有影响;管桩壁越厚,桩周土的剪切模量越大时,管桩水平动力阻抗的绝对值越大。     

Consolidation analysis of saturated multi‐layered soils with anisotropic permeability caused by a point sink

This paper presents the analytical layer‐element method to analyze the consolidation of saturated multi‐layered soils caused by a point sink by considering the anisotropy of permeability. Starting from the governing equations of the problem, the solutions of displacements and stresses for a single soil layer are obtained in the Laplace–Hankel transformed domain. Then, the analytical layer‐element method is utilized to further derive the solutions for the saturated multi‐layered soils in the transformed domain by combining with the boundary conditions of the soil system and continuity conditions between adjacent layers. The actual solutions in the physical domain can be acquired by the inversion of Laplace–Hankel transform. Numerical results are carried out to show the accuracy and stability of the proposed method and evaluate the influence of sink depth and anisotropic permeability on excess pore pressure and surface settlement. Copyright © 2011 John Wiley & Sons, Ltd.     

Three‐dimensional finite element analyses of passive pile behaviour

Piles may be subjected to lateral soil pressures as a result of lateral soil movements from nearby construction‐related activities such as embankment construction or excavation operations. Three‐dimensional finite element analyses have been carried out to investigate the response of a single pile when subjected to lateral soil movements. The pile and the soil were modelled using 20‐node quadrilateral brick elements with reduced integration. For compatibility between the soil–pile interface elements, 27‐node quadrilateral brick elements with reduced integration were used to model the soil around the pile adjacent to the soil–pile interface. A Mohr–Coulomb elastic–plastic constitutive model with large‐strain mode was assumed for the soil. The analyses indicate that the behaviour of the pile was significantly influenced by the pile flexibility, the magnitude of soil movement, the pile head boundary conditions, the shape of the soil movement profile and the thickness of the moving soil mass. Reasonable agreement is found between some existing published solutions and those developed herein. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

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.     

Paleo‐fires and Atmospheric Oxygen Levels in the Latest Permian: Evidence from Maceral Compositions of Coals in Eastern Yunnan,Southern China

Inertinite maceral compositions of the Late Permian coals from three sections in the terrestrial and paralic settings of eastern Yunnan are analyzed in order to reveal the paleo-fire events and the atmospheric oxygen levels in the latest Permian. Although the macerals in the studied sections are generally dominated by vitrinite, the inertinite group makes up a considerable proportion. Its content increases upward from the beginning of the Late Permian to the coal seam near the Permian-Triassic boundary. Based on the microscopic features and the prevailing theory that inertinite is largely a by-product of paleo-fires, we suggest that the increasing upward trend of the inertinite abundance in the latest Permian could imply that the Late Permian peatland had suffered from frequent wildfires. Since ignition and burning depend on sufficient oxygen, a model-based calculation suggests that the O2 levels near the Wuchiapingian/Changhsingian boundary and the Permian-Triassic boundary are 27% and 28% respectively. This output adds supports to other discoveries made in the temporal marine and terrestrial sediments, and challenges the theories advocating hypoxia as a mechanism for the Permian-Triassic boundary crisis.     

基于2.5维有限元饱和土地基中空沟隔振性能研究

基于Biot饱和多孔介质U-W格式动力控制方程,采用Galerkin法和Fourier变换,推导了饱和土体2.5维有限元方程及黏弹性人工边界,建立了饱和土地基中空沟分析模型,并在波数域中进行求解,通过快速Fourier变换（FFT）进行波数展开,获得三维空间域中结果。算例分析了移动荷载作用下均质饱和土地基、分层饱和土地基、上覆单相弹性层饱和土地基3种饱和土地基模型中空沟的隔振效果。结果表明：饱和土地基中空沟的隔振效果不仅与空沟自身深度有关,还与地基中成层土体的分界面以及土体参数有关,波在不同土体分界面上的透射和反射会影响空沟的隔振效果;饱和土地基中上覆单相弹性层厚度对空沟的隔振效果影响显著,随着上覆单相弹性层厚度的增加,饱和土地基中空沟的隔振效果变好。     

PML吸收边界条件影响因素分析

影响完全匹配层方法吸收效果的主要因素有吸收层厚度和衰减系数,笔者通过地震波场数值模拟,讨论了吸收层厚度和衰减系数对吸收效果的影响。研究表明,在衰减系数一定时,吸收层厚度越大,边界反射越弱。吸收层的厚度一般在12～20道较为适当,在吸收层的厚度比较小时,随着衰减系数的增大,边界反射逐渐减弱甚至完全没有边界反射,但是若再进一步增大衰减系数则又会逐步出现边界反射。因此,在实际应用中应该注意衰减系数的选择,衰减系数一般在500～2 000之间,衰减系数较大反而影响其效果。     

Modelling strain localization in granular materials using micropolar theory: numerical implementation and verification

Implementation and applications for a constitutive numerical model on F‐75 silica sand, course silica sand and two sizes of glass beads compressed under plane strain conditions are presented in this work. The numerical model is used to predict the stress versus axial strain and volumetric strain versus axial strain relationships of those materials; moreover, comparisons between measured and predicted shear band thickness and inclination angles are discussed and the numerical results compare well with the experimental measurements. The numerical model is found to respond to the changes in confining pressure and the initial relative density of a given granular material. The mean particle size is used as an internal length scale. Increasing the confining pressure and the initial density is found to decrease the shear band thickness and increase the inclination angle. The micropolar or Cosserat theory is found to be effective in capturing strain localization in granular materials. The finite element formulations and the solution method for the boundary value problem in the updated Lagrangian frame (UP) are discussed in the companion paper. Copyright © 2006 John Wiley & Sons, Ltd.     

Time‐dependent behaviour of a rigid foundation on a transversely isotropic soil layer

An analytical solution is presented in this paper to study the time‐dependent settlement behaviour of a rigid foundation resting on a transversely isotropic saturated soil layer. The governing equations for a transversely isotropic saturated soil, within Biot's poroelasticity framework, are solved by means of Laplace and Hankel transforms. The problem is subsequently formulated in the Laplace transform domain in terms of a set of dual integral equations that are further reduced to a Fredholm integral equation of the second kind and solved numerically. The developed analytical solution is validated via comparison with the existing analytical solution for an isotropic saturated soil case, and adopted as a benchmark to examine the sensitivities of the mesh refinement and the locations of truncation boundaries in the finite element simulations using ABAQUS. Particular attention is paid to the influences of the degree of soil anisotropy, boundary drainage condition, and the soil layer thickness on the consolidation settlement and contact stress of the rigid foundation. Copyright © 2009 John Wiley & Sons, Ltd.     

A quasistatic analysis of a plate on consolidating layered soils by analytical layer‐element/finite element method coupling

In this paper, a coupling method between finite element and analytical layer‐elements is utilized to analyze the time‐dependent behavior of a plate of any shape and finite rigidity resting on layered saturated soils. Based on the integral transform techniques together with the aid of an order reduction method, an analytical layer‐element solution is derived from the governing equations for three‐dimensional Biot consolidation with respect to a Cartesian coordinate system and then extended to be the fundamental solution for the layered saturated soil under a point load. The Mindlin plate is modeled by eight‐noded isoparametric elements. The governing equations of the interaction between soil and plate in the Laplace‐Fourier transformed domain are deduced by referring to the coupling theory of FEM/BEM, and the final solution is obtained by applying numerical inversion. Numerical examples concerned with the time‐dependent response of a plate are performed to demonstrate the influence of soil and plate properties on the interaction process. Copyright © 2014 John Wiley & Sons, Ltd.     

应用亮兰染色剂指示溶质迁移边界层的研究

应用边界层方法确定溶质迁移参数是一种简单可靠的方法,但如何确定边界层的位置是限制边界层方法应用的关键因素。而亮兰染色剂在土壤中迁移时,具有边界层清晰、容易确定的特点。基于此,以Cl-为例,研究了该染色剂指示非反应性溶质迁移边界层的可行性。研究结果表明:在同一土柱中,亮兰染色剂的边界层与Cl-的边界层的运动规律相似,亮兰染色剂和Cl-的边界层随时间的变化在弥散度很小的情况下都可以用线性关系式近似描述。根据这一规律,得出在风沙土和黄绵土中Cl-边界层与亮兰染色剂边界层之间的关系表达式。两者边界层之间转换关系的求得,为应用亮兰染色剂准确指示非反应性溶质迁移边界层成为可能提供了一条思路。     

Three‐dimensional nonlinear finite element analysis of pile groups in saturated porous media using a new transmitting boundary

The dynamic behaviour of pile groups subjected to an earthquake base shaking is analysed. An analysis is formulated in the time domain and the effects of material nonlinearity of soil, pile–soil–pile kinematic interaction and the superstructure–foundation inertial interaction on seismic response are investigated. Prediction of response of pile group–soil system during a large earthquake requires consideration of various aspects such as the nonlinear and elasto‐plastic behaviour of soil, pore water pressure generation in soil, radiation of energy away from the pile, etc. A fully explicit dynamic finite element scheme is developed for saturated porous media, based on the extension of the original formulation by Biot having solid displacement (u) and relative fluid displacement (w) as primary variables (u–w formulation). All linear relative fluid acceleration terms are included in this formulation. A new three‐dimensional transmitting boundary that was developed in cartesian co‐ordinate system for dynamic response analysis of fluid‐saturated porous media is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb surface waves as well as body waves. The pile–soil interaction problem is analysed and it is shown that the results from the fully coupled procedure, using the advanced transmitting boundary, compare reasonably well with centrifuge data. Copyright © 2007 John Wiley & Sons, Ltd.