夹硬层土中单桩贯入挤土效应的数值模拟

桩贯入土体产生的挤土效应问题较为复杂。利用ABAQUS软件建立了单桩贯入夹硬层土和均质土的二维轴对称有限元模型,经过分析比较,得出了单桩贯入夹硬层土体所特有的位移场及应力场的变化规律。分析表明：桩贯入夹硬层土过程中,软硬土层交界处土体水平位移变化剧烈;硬土层的存在,会使土体水平及竖向位移受到约束;夹硬层土的水平挤压应力要远大于均质土情况;与水平应力相比,竖向挤压应力在硬土层处明显偏小。     

埋置简谐扭转荷载作用下广义Gibson饱和地基动力响应

考虑地基为饱和半空间,研究了广义Gibson饱和地基内作用简谐扭转动荷载时地基的动力响应问题。从Biot饱和地基固结理论出发,结合扭转振动的特点,建立了剪切模量随深度线性变化的饱和地基扭转振动的动力微分方程,通过Hankel变换求解此微分方程,给出了Hankel变换域内的切向位移和剪应力关于待定系数的表达式。根据饱和地基表面为自由表面,荷载作用面位移连续、剪应力差等于动荷载大小,波的辐射条件等边界条件求解出待定系数,借助Hankel逆变换给出地基内的位移和应力的表达式。通过数值算例研究发现:在同一水平面内,地基内的切向位移和剪应力曲线的实部和虚部都呈现出非常明显的波动变化规律;在竖向平面内,动荷载作用面上部区域内随深度逐渐增大时,地基内切向位移和剪应力曲线的实部逐渐增大,而在动荷载作用面下部区域则正好相反;扭转动荷载的影响范围主要是荷载作用面上下2倍半径区域。     

水平简谐荷载作用下层状饱和土体动力响应

根据Biot波动理论,采用传递、反射矩阵（TRM）方法研究了水平简谐荷载作用下层状饱和土动力响应问题。由Helmholtz矢量分解求出基本解,再利用TRM法推导了层状饱和土动力响应,并由数值Hankel逆变换得到层状土地基位移、应力及孔压在空间域内的解。利用计算结果与已有结果相比较,二者相吻合,验证了算法的正确性。算例分析表明,水平简谐荷载作用在有软弱夹层的层状土体中比均质土中具有更显著的动力响应,尤其是软夹层上下有硬土层时,会引起软弱夹层土体孔隙水压升高、位移幅值增大、土体波动性增强;而荷载作用硬夹层及夹层上下有软土层时,情况则相反。     

Research on Construction Technology and Formation Deformation of Shield Tunneling Under River

In view of the construction technology and formation deformation of the existing rivers under the subway shield interval tunnel, relying on Zhengzhou rail transit line 17 shield tunnel through the south-north water transfer channel project, using the trial tunneling method of shield tunnel test section to optimize all kinds of construction parameters, and the formation analysis and calculation of the determined digging pressure, noting the final parameters that are determined by comparison between the extractive parameters such as pulp pressure, and they are used in MIDAS GTS/NX finite unit analysis software for numerical simulation analysis, combined with on-site testing, to study the variation of different depth formation subsidence values and horizontal displacement values when passing through the south-north water transfer canal under water conditions, and put forward the reinforcement scheme of shield tunneling through the general river bed and hole. The results show that the maximum subsidence of the trunk canal shield through the bottom is 11.3 mm, and the verification parameters are reasonable and feasible. The horizontal subsidence trough of the soil above the over-soil layer at the top of the main canal embankment is distributed by Gauss. The sinking distribution can be estimated using the Peck formula, the middle of the layer subsidence trough is distributed horizontally, and the soil subsidence trough below it is inverted hump-like distribution, with the peak point directly above the center line of the two tunnels, and the vertical zero horizontal displacement surface is located in the middle line position between the two tunnels, and the horizontal zero horizontal displacement surface is located. The horizontal displacement curve of the soil layer above the soil layer is inverted S-type, the maximum displacement appears at the anti-bending point of the surface sedimentation trough curve, the horizontal displacement curve of the soil layer on the horizontal zero horizontal displacement surface is zero, and the lower soil layer curve is double-inverted S-shaped. In the construction, we should pay attention to the monitoring of the vertical subsidence and horizontal displacement of the deep soil, adjust the shield parameters and reduce the horizontal and vertical shearing effect of the deep soil on the surrounding pipeline, inner pile base and other structures during the construction of the shield.

     

Displacement and stress distributions under a uniform inclined rectangular load on a cross‐anisotropic geomaterial

In practical engineering, an applied rectangular area load is not often horizontally or vertically distributed but is frequently inclined at a certain angle with respect to the horizontal and vertical axes. Thus, the solutions of displacements and stresses due to such a load are essential to the design of foundations. This article yields the analytical solutions of displacements and stresses subjected to a uniform rectangular load that inclines with respect to the horizontal and vertical axes, resting on the surface of a cross‐anisotropic geomaterial. The planes of cross‐anisotropy are assumed to be parallel to the horizontal ground surface. The procedures to derive the solutions can be integrated the modified point load solutions, which are represented by several displacement and stresses elementary functions. Then, upon integrations, the displacement and stress integral functions resulting from a uniform inclined rectangular load for (1) the displacements at any depth, (2) the surface displacements, (3) the average displacements in a given layer, (4) the stresses at any depth, and (5) the average stresses in a given layer are yielded. The proposed solutions are clear and concise, and they can be employed to construct a series of calculation charts. In addition, the present solutions clarify the load inclinations, the dimensions of a loaded rectangle, and the analyzed depths, and the type and degree of geomaterial anisotropy profoundly affect the displacements and stresses in a cross‐anisotropic medium. Parametric results show that the load inclination factor should be considered when an inclined rectangular load uniformly distributed on the cross‐anisotropic material. Copyright © 2008 John Wiley & Sons, Ltd.     

土钉墙开挖性状的有限元分析

土钉墙是粉砂土、黏性土等地区基坑开挖围护的主要形式之一。采用ABAQUS有限元软件建立土钉-面层-土体的相互作用模型,土体应用Mohr-Coulomb本构模型,土体与土钉间的接触性质为摩擦小滑移,在此基础上,得到了土体深层水平位移、面层位移、面层土压力和土钉轴力的分布形式和基本规律。计算结果表明：深层土体水平位移的最大值发生在墙顶靠下或者墙顶的位置;坑底水平位移从坡脚向坑中逐渐变小,在基坑开挖深度范围内,减小的速率较大,以后趋于某一稳定值;土钉轴力两头较小,中间较大,作用在面层上的作用力值远小于土钉轴力,该力与作用在面层上的作用力之和等于土钉轴力;土钉墙的面层土压力随着深度的增加,先增加到最大值后,再逐渐减小,每开挖一步面层土压力就会增加,且最大值向下移动。     

基坑工程平面竖向弹性地基梁法中土的水平抗力比例系数反分析研究

提出了一种结合Ucode反分析软件和Abaqus有限元分析软件,并根据围护墙实测变形来反分析基坑工程平面竖向弹性地基梁法中土的水平抗力比例系数的方法。重点讨论了反分析的原理和实现流程,给出了反分析过程中迭代控制参数的确定方法及待估参数取值的确定方法。通过参数敏感分析可找到对围护结构变形影响较大的土层,针对起主导作用的土层m值进行反分析,既可提高迭代效率,又能保证反分析结果的准确性。实际工程反分析的结果表明,计算值与实测围护墙变形吻合得很好。在此基础上,根据多个测孔的反分析结果给出了各层土的水平抗力比例系数的建议取值。     

自钻式原位剪切旁压模型PFC3D细观模拟

自钻式原位剪切旁压试验(self-boring in-situ shear pressure-meter)以其独特的多级加载方式, 能够直接测出土体的强度和变形参数, 然而, 目前对探头周围土体变形机理研究较少.为揭示自钻式原位剪切旁压仪试验过程中测定器周围土颗粒变形机理, 应用PFC3D(particle flow code in three dimensions)颗粒流程序对自钻式原位剪切旁压试验进行了仿真数值模拟, 对多级加载过程中探头周围土体的位移场和应力场发展变化以及数值试样各阶段变形模量和土颗粒运动轨迹进行了分析研究.试验结果表明: 随着剪应力的逐级施加, 中间区域颗粒的位移量不断增大, 且位移矢量方向性更加显著.径向应力在探头附近两侧形成近似呈对称分布的数个"应力核"; 竖向应力在探头两侧形成扁平状应力带, 且在肩部形成应力集中区.中部区域球颗粒的运动轨迹成台阶状, 且随距探头距离的增大由近到远可分为3个特征区域; 球颗粒的Z向和XY向位移量亦随之呈负指数形式衰减, Z向位移量衰减速率更快.      

Three-dimensional finite element modelling for consolidation due to groundwater withdrawal in a desaturating anisotropic aquifer system

A poroelastic numerical model is presented to evaluate three-dimensional consolidation due to groundwater withdrawal from desaturating anisotropic porous media. This numerical model is developed based on the fully coupled governing equations for groundwater flow in deforming variably saturated porous media and the Galerkin finite element method. Two different cases of unsaturated aquifers are simulated for the purpose of comparison: a cross-anisotropic soil aquifer, and a corresponding isotropic soil aquifer composed of a geometrically averaged equivalent material. The numerical simulation results show that the anisotropy has a significant effect on the shapes of three-dimensional hydraulic head distribution and displacement vector fields. Such an effect of anisotropy is caused by the uneven partitioning of the hydraulic pumping stress between the vertical and horizontal directions in both groundwater flow field and solid skeleton deformation field. Copyright © 1999 John Wiley & Sons, Ltd.     

Vertical impedance of an end‐bearing pile in viscoelastic soil

This paper presents a new method to derive the analytical solution for the vertical impedance of an end‐bearing pile in viscoelastic soil. The soil is assumed as a homogeneous and isotropic layer, and the pile is considered as a one‐dimensional Euler rod. Considering both the vertical and radial displacements of soil and soil–pile coupled vibration, the governing equations of the soil and pile are established. The volumetric strain of soil is obtained by transformation on the equations of soil and variable separation method. Then the vertical and radial displacements of soil are obtained accordingly. The displacement response and impedance function of pile are derived based on the continuity assumption of the displacement and stress between the pile and soil. The solution is verified by being compared with an existing solution obtained by introducing potential functions. Furthermore, a comparison with two other simplified solutions is conducted. Numerical examples are presented to analyze the vibration characteristics of the pile. Copyright © 2014 John Wiley & Sons, Ltd.     

倾斜土层地震液化诱发大变形数值分析方法

基于有限元理论,进行了缓坡土层地震液化引起大变形的数值方法研究,即采用二维有效应力动力有限元方法进行分析,在液化分析过程的每一时段考虑地震液化和振动软化得到土单元的模量,通过非线性静力方法计算每时段地震液化引起的大变形,得到土层各深度处的水平和竖向位移。由算例分析了地震动和土层坡度等因素的影响,通过对比分析表明了该方法的有效性,可为工程场地地震地质灾害评价提供参考数据。     

考虑基坑卸载影响的改进弹性抗力法及其应用

弹性抗力法是目前实际工程中基坑支护结构变形内力计算的一种主要方法,该方法借鉴的是计算水平荷载作用下竖直放置弹性桩的地基反力系数法,但没有考虑开挖卸载对土抗力系数和坑内土体的变形影响,因而导致变形计算值通常小于实测值。首先在开挖面施加反向自重应力模拟基坑的开挖卸荷,利用Mindlin解得到坑内土体的卸荷应力分布;在前人的实验研究基础上,推导了不同卸荷比下考虑卸荷影响的土体水平抗力系数;将支护结构的变形分解为开挖瞬时的卸载作用和基坑内外水土压力差作用2个阶段,考虑卸载应力路径对坑内土体水平抗力系数以及开挖瞬时变形的影响,对传统的弹性抗力法进行了改进,最后将该方法应用于典型三维基坑算例。研究发现,改进法得到的支护结构水平变形大于传统方法,而对内力结果影响不大;基坑平面尺寸越大,两者的变形结果差别越明显,改进法更能反映平面尺寸对基坑变形的影响     

地震动力作用下土-地铁隧道模型分析

针对地铁建设中典型的马蹄形断面隧道建立比例尺为1∶20的分析模型,并采用数值分析方法研究马蹄形隧道处于单一土层及工程所处区域典型的成层土体中时的动力响应。分析结果表明,在单一土层中由于土体的约束作用,结构产生的位移以整体沉降为主,在成层土体中除产生一定的整体变形外还伴随一定的扭转变形。在两种地层情况下马蹄形地铁隧道在地震动力作用下的动力加速度响应、竖向位移均在拱顶处产生最大值,其中在单一土层中的加速度响应最大值为结构中部加速度的2.29倍。结构在顶部和侧板处所产生的动应力响应值也较大。研究表明,地震动力荷载作用下顶板、侧板均为受力较大部位,在设计和施工中应予以充分重视。     

直角域中凸起和孔洞对SH波的散射与地震动

按照波函数展开法和镜像方法,对直角域中半圆形凸起和圆形孔洞对SH波的散射进行了分析,得到其稳态解。对含孔洞和凸起的直角域做分区,等效为一个含孔洞与凹陷的直角域和一个圆域的契合,其在分界面上满足位移和应力的连续性条件,即契合条件,分别构造两个区域内的位移波函数,按照孔洞边界柱面上的应力自由和契合条件定解波函数展开式的系数。按Fourier级数展开法,得到定解条件的线性代数方程组,截断求解,进而得到问题的解析解。数值算例给出圆形孔洞边沿动应力和地表位移幅值的分布情况,得到直角域自由边界、凸起、孔洞对散射和地震动的影响。     

Elastic solutions for a transversely isotropic half-space subjected to buried asymmetric-loads

Elastic closed-form solutions for the displacements and stresses in a transversely isotropic half-space subjected to various buried loading types are presented. The loading types include finite line loads and asymmetric loads (such as uniform and linearly varying rectangular loads, or trapezoidal loads). The planes of transverse isotropy are assumed to be parallel to its horizontal surface. These solutions are directly obtained from integrating the point load solutions in a transversely isotropic half-space, which were derived using the principle of superposition, Fourier and Hankel transformation techniques. The solutions for the displacements and stresses in transversely isotropic half-spaces subjected to linearly variable loads on a rectangular region are never mentioned in literature. These exact solutions indicate that the displacements and stresses are influenced by several factors, such as the buried depth, the loading types, and the degree and type of rock anisotropy. Two illustrative examples, a vertical uniform and a vertical linearly varying rectangular load acting on the surface of transversely isotropic rock masses, are presented to show the effect of various parameters on the vertical surface displacement and vertical stress. The results indicate that the displacement and stress distributions accounted for rock anisotropy are quite different for those calculated from isotropic solutions. Copyright © 1999 John Wiley & Sons, Ltd.     

Elastic solutions for a transversely isotropic half-space subjected to a point load

We rederive and present the complete closed-form solutions of the displacements and stresses subjected to a point load in a transversely isotropic elastic half-space. The half-space is bounded by a horizontal surface, and the plane of transverse isotropy of the medium is parallel to the horizontal surface. The solutions are obtained by superposing the solutions of two infinite spaces, one acting a point load in its interior and the other being free loading. The Fourier and Hankel transforms in a cylindrical co-ordinate system are employed for deriving the analytical solutions. These solutions are identical with the Mindlin and Boussinesq solutions if the half-space is homogeneous, linear elastic, and isotropic. Also, the Lekhnitskii solution for a transversely isotropic half-space subjected to a vertical point load on its horizontal surface is one of these solutions. Furthermore, an illustrative example is given to show the effect of degree of rock anisotropy on the vertical surface displacement and vertical stress that are induced by a single vertical concentrated force acting on the surface. The results indicate that the displacement and stress accounted for rock anisotropy are quite different for the displacement and stress calculated from isotropic solutions. © 1998 John Wiley & Sons, Ltd.     

Neotectonic and modern stresses of South Sakhalin

Tectonophysical studies are conducted in South Sakhalin for identification of temporal–spatial changes in the geodynamic settings of the formation of the local structures. Analysis of the field data reveals 11 local stress state (LSS) in the large newest megastructures, which were formed on geological basement of different ages. The parameters of the tectonic stresses are significantly distinct in each LSS, especially the orientations (up to reindexation) of the compression and extension axes in different fault wings. Tectonic stresses of two ages and constant latitudinal and horizontal compression axis are reconstructed. The earlier cofold shear stress field with a horizontal and longitudinal extension axis is post-Miocene and the later stress field of the reversed fault with a vertical extension axis is orogenic. The LSSs reconstructed for the first time by the displacement vectors on slickensides, along with the data on the earthquake mechanisms, substantiate the reindexation of the horizontal extension axis with the vertical intermediate axis of major normal stresses at the postfold orogenic stage of evolution of the territory. These results are in agreement with previous data on the transformation of the dextral to reverse thrust displacements along the longitudinal fault systems. The young stress field is more confidently interpreted in the activation fault zones, which limit the orogenic blocks, whereas the traces of cofold deformations without younger orogenic stress fields better remain inside the blocks which are composed of older and strongly dislocated Mesozoic rocks.     

多次冻融后格栅加筋黏土路堤工作性能分析

利用多次冻融后加筋黏土路堤相关试验参数对其进行有限元计算, 确定了加筋黏土结构各初始因素对其工作性能的影响规律, 并利用ABAQUS有限元计算软件对各种计算条件下的位移和应力进行求解, 确定了冻融循环后位移和应力的最大值. 结果显示: 增加加筋黏土路堤中的格栅层数可以减小冻融后加筋路堤的竖向和水平位移及土体剪应力最大值, 且格栅层数的增加对土体变形的影响强于对土体剪应力值的影响. 土体压实度的增大可以减小加筋路堤的竖向和水平位移, 同时引起土中剪应力数值的增大; 土体初始含水率的减小可以减小加筋路堤的竖向和水平位移, 同时引起土中剪应力数值的减小.     

Subsidence due to pumping from layered soil—a perturbation theory

Persistent pumping of groundwater is known to cause subsidence of the ground surface. When the soil is composed of alternating layers of sand and clay, subsidence is, in general, a consequence of pressure change in and deformation of all layers. Even with various simplifications, the mathematical task for the mechanics of a layered soil is complex and, in the past, many idealizations have been introduced in the literature. Among the existing approximations, a well-known approximation by Hantush¹ for well-hydraulics in leaky aquifers is intuitively attractive. He used the fact known to be strictly true for simple flows in a soil system with horizontal layers, that the flow must be mainly horizontal in a relatively porous aquifer and mainly vertical in a highly impervious aquitard, and invoked the same approximation for transient flows induced by well-pumping. In this paper we shall apply a perturbation analysis to give a mathematical confirmation of the hydrological approximation and to show its degree of accuracy. In particular, we treat a three-layered system, where the middle layer is a relatively soft aquitard, and the well withdraws water from the artesian aquifer at the bottom. The upper aquifer has a phreatic surface. Physically, we examine the often ignored effects of self-weight which is important for thick aquitards. The perturbation theory is worked out for sufficiently weak pumping rates and small soil deformation. Subsidence in the vertical direction due mainly to the deformation of the soft aquitard is studied. Results show that the common assumption of constant total stress with respect to depth is valid only when the storage coefficient is very small.     

Deformation of a layered half-space due to a very long tensile fault

The problem of the coseismic deformation of an earth model consisting of an elastic layer of uniform thickness overlying an elastic half-space due to a very long tensile fault in the layer is solved analytically. Integral expressions for the surface displacements are obtained for a vertical tensile fault and a horizontal tensile fault. The integrals involved are evaluated approximately by using Sneddon’s method of replacing the integrand by a finite sum of exponential terms. Detailed numerical results showing the variation of the displacements with epicentral distance for various source locations in the layer are presented graphically. The displacement field in the layered half-space is compared with the corresponding field in a uniform half-space to demonstrate the effect of the internal boundary. Relaxed rigidity method is used for computing the postseismic deformation of an earth model consisting of an elastic layer of uniform thickness overlying a viscoelastic half-space.