SURFACE DISPLACEMENTS OF A NON-HOMOGENEOUS ELASTIC HALF-SPACE SUBJECTED TO UNIFORM SURFACE TRACTIONS. PART I: LOADING ON ARBITRARILY SHAPED AREAS

A numerical technique is presented for the analysis of surface displacements of a non-homogeneous elastic half-space subjected to vertical and/or horizontal surface loads uniformly distributed over an arbitrarily shaped area. The non-homogeneity considered is a particular form of power variation of Young's modulus with depth. Since the exponent which determines the degree of non-homogeneity may vary from zero to unity, both the homogeneous half-space and the Gibson soil may be included as limiting cases in a single numerical scheme. In order to account for the arbitrary shape of the loading, the boundary of the loaded area is linearized piecemeal. This enables the modeling of any load pattern according to the desired degree of accuracy. Special attention is focused on the integration scheme, since the singularity associated with the Green's function becomes progressively more pronounced the greater the non-homogeneity parameter gets. The performance of the numerical procedure is studied using analytical solutions for rectangular shaped areas. Further comparisons with well-known solutions based on integral transform techniques for a uniformly distributed load acting on a circular area of the non-homogeneous soil mass show excellent agreement as well. © 1997 John Wiley & Sons, Ltd.     

The P-y Response of Laterally Loaded Flexible Piles in Residual Soil

This paper describes the main features related to lateral displacements with depth after successive lateral loading–unloading cycles applied to the top of reinforced-concrete flexible bored piles embedded in naturally bonded residual soil. The bored piles under study have a cylindrical shape, with 0.40-m in diameter and 8.0-m in length. Both bored piles types (P1 and P2) include an embedded steel pipe section in their center as longitudinal steel reinforcements: pile type P1 has another 16 steel rods as steel reinforcement to concrete while pile type P2 has no further steel reinforcement. Pile type P1 has three times as much stiffness (EI) and four and a half times the plastic moment (M_y) than pile type P2. A similar load–displacement performance was observed at initial loads as for small displacements of both piles. At this initial loading stage, the response of the reinforced concrete piles is a function of the soil characteristics and of a linear elastic pile deformation. During this stage, piles can even be understood as probes for evaluating soil reactions. For larger horizontal displacements, after the concrete section starts undergoing large deformations, approaching the ultimate bending moment, pile behavior and consequently the load–displacement relation starts to diverge for both piles. For pile P1 the values of relevant lateral displacements are extended to about 2.5-m in depth, while for pile P2 lateral displacements are mostly constrained to about 2.0-m in depth. Measurements of horizontal displacements of pile P1 against depth recorded with a slope indicator show that, after unloading, lateral loads at distinct stages (small and near failure loads), exhibits a much higher elastic phase of the system response. An analytical fitting model of soil reaction is proposed based on the measured displacements from slope indicator. The integration of a continuous model proposed for the soil reaction agrees fairly well with the measured displacements up to moments close to plastic limit. Results of load–displacement show that the stiffer pile (P1) was able to mobilize twice as much lateral load compared to pile P2 for a service limit displacement of about 20 mm. The paper shows results that enable the isolation of the structural variable through real scale pile load tests, thus granting understanding of its importance and enabling its quantitative visualization in examples of piles embedded in residual soil sites.

     

TBM护盾−围岩相互作用荷载识别方法

随着隧道掘进机（tunnel boring machine, TBM）开挖隧道埋深逐渐加大,穿越的地层结构更加复杂,如何降低卡机风险成为亟待解决的问题。提出了一种基于牛顿迭代和有限元理论的TBM护盾?围岩相互作用荷载反演方法。该反演方法以护盾内表面测点应变测量值与计算值的误差作为目标函数,并引入Moore-Penrose广义逆对反演方程进行求解。采用了护盾外表面有限分区的策略,定义各分区节点荷载值为反演参数,可以方便地控制反演参数的数量。各分区内部荷载通过节点值插值得到,实现了护盾表面任意分布荷载的离散拟合。算例结果表明：该反演方法计算效率高,对观测值误差具有很好的免疫效果;能够有效地获取护盾外表面的荷载分布规律和摩阻力的大小,在一定程度上能够为现场安全掘进、卡机预警以及事故处理提供指导。同时,该反演算法也可以扩展到低维情况,对其他结构受力或者物理参数的反演具有一定的借鉴意义。     

Stresses and displacements in an elastic mass due to distributed loading by discretization technique

Numerical solutions have been obtained for stresses and displacements in a linear elastic half space due to distributed loads of circular, rectangular and elliptical shapes. The technique primarily involves use of a multi-dimensional numerical integration technique to integrate point load solutions over the distributed loading after discretizing the area into a finite number of elements. Both uniform vertical and shear loads have been considered as well as vertical conical loads and inward shear loads. The technique evolved facilitates the determination of stresses and displacements by the use of mini-computers and is neither as tedious and cumbersome as the use of tables and charts nor as costly as FEM solutions. A detailed comparison has been presented between the results obtained by the numerical solutions and those of the existing analytical solutions wherever they are available. It is found that the agreement between the two is within one per cent for displacements at all depths for the different cases studied. The matching is also good in the case of stresses, except at shallow depths.     

Dynamic characteristics of multiscale longitudinal stress and particle rotation in ballast track under vertical cyclic loads

Loads transfer in ballast track through contacts among randomly distributed ballast particles and have strong heterogeneity. Since the size ratio between ballast track and ballast particles is generally small, using averaged stress to describe the internal mechanical state in ballast track faces practical difficulties. For example, particle movements and high local concentration stress tend to be ignored. The inter-particle contact stress is crucial to evaluate the particle behaviors, such as abrasion, movements, and furtherly the performance of ballast track. However, the contact stress on ballast particles is hard to predict or measure. We conduct a full size model test to investigate the dynamic characteristics of longitudinal stress on ballast particles as well as different lateral regions under vertical cyclic loads with various loading magnitudes and frequencies. An obvious seesaw effect of longitudinal contact stress is observed: the stresses at some contact areas have the same phase with applied cyclic load while at other contact areas have an opposite phase. The seesaw effect of contact stress is then used to evaluate the rotational movements of ballast particles. The variation of contact area and stress of the ballast particles with loading magnitudes demonstrates that the rigid contact assumption is appropriate when analyzing the contact behavior of ballast particles. The cumulative probability distribution of contact stress with stress level can be described by an inversely proportional function, based on which the maximum contact stress can be estimated according to the longitudinal average stress. Besides, the lateral dispersion angle of the vertical loads in the ballast track is about 35°, which is independent of the given loading magnitudes and frequencies.

     

Cyclic shakedown of piles subjected to two‐dimensional lateral loading

Pile foundations are frequently subjected to cyclic lateral loads. Wave and wind loads on offshore structures will be applied in different directions and times during the design life of a structure. Therefore, the magnitude and direction of these loads in conjunction with the dead loads should be considered. This paper investigates a loading scenario where a monotonic lateral load is applied to a pile, followed by two‐way cycling in a direction perpendicular to the initial loading. This configuration is indicative of the complexity of loading that may be considered and is referred to in the paper as ‘T‐shaped’ loading. The energy‐based numerical model employed considers two‐dimensional lateral loading in an elasto‐plastic soil, with coupled behaviour between the two perpendicular directions by local yield surfaces along the length of the pile. The behaviour of the soil–pile system subjected to different loading combinations has been divided into four categories of shakedown previously proposed for cyclic loading of structures and soils. A design chart has been created to illustrate the type of pile behaviour for a given two‐dimensional loading scenario. Copyright © 2009 John Wiley & Sons, Ltd.     

横观各向同性层状地基动应力响应的求解与分析

动力响应问题的求解对于地基在外荷载作用下引起的弹性波动问题研究有重要的意义。本文提出了一种求解横观各向同性层状地基在施加时间简谐荷载作用下任意点的应力响应的算法。此算法利用傅里叶变换将广义平面应变问题频率-空间域的动力方程转化到频率-波数域内,结合对偶变量的引入,利用高精度的精细积分算法对状态方程进行求解,在得到频率-波数域内的位移响应的基础上,利用傅里叶逆变换得到任意点的动应力响应。简谐荷载不仅可以施加在地基表面,而且可以施加在地基内部。对比算例验证了本文算法的准确性,同时对地基各向异性特性、激励频率和阻尼比对动应力响应的影响进行了参数分析,为工程实际提供可靠的数值依据。     

基于两阶段法的堆载对公路桥梁桩基础影响分析

首先基于布西奈斯克解,将表面作用有集中荷载时半无限弹性体的应力变形解在荷载作用区域利用复合辛普森公式进行数值积分,得到了表面处于不同形式的分布荷载作用下土体自由场水平及竖直方向应力和变形的计算公式。结合有限差分法及Mindlin位移解分别推导出弹性地基中桩顶作用有集中荷载时单桩桩、土单元的位移,两者结合便得到主动桩的分析方法。通过将堆载作用下的土体自由场位移施加于桩上把上述自由场分析及主动桩分析结合起来,推导出堆载作用下被动单桩的竖向和水平向承载特性的分析方法。选取地表附近作用有矩形均布荷载的计算工况,通过与有限元计算结果的对比验证了所提出的分析方法的正确性;将所提出的方法计算结果与某处工程案例中实测数据结果进行对比,分析表明利用该方法分析堆载对邻近桩基的水平及竖向影响是合理可靠的。     

Flexural Analysis in Dynamic Pinned Head Pile Testing

A Laplace transform is used to solve the problem of the steady state and transient response of a pinned head pile embedded into a viscoelastic Winkler soil medium. The pile is modeled as an Euler–Bernoulli beam while the soil medium is modeled using a Winkler subgrade approach. Two analytical solutions are developed to specifically address both steady state and transient loads encountered during dynamic pile testing. After choosing a proper contour integration in the complex plane, inverse integration is evaluated. The steady state solutions are associated to the residues of the integration around the poles while the transient solutions are associated to the integration paths along the contour integration. The derived solutions are applied to a case history for which results of dynamic pile tests are available. Dynamic pile flexion is generated by delivering eccentric impact using a dynamic loading test module. Validity of the proposed solution is discussed basing on geotechnical campaign and recorded pile head bending moment and rotation rate.     

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.     

堆载下单桩负摩阻力工作性状非线性数值分析

利用三维非线性数值方法对一算例在堆载作用下桩侧负摩阻力进行了计算分析。计算结果揭示了摩擦型、端承型桩负摩阻力工作性状的异同;分析了堆载速度、桩顶荷载对负摩阻力的影响:堆载越快,负摩阻力越小;在无桩顶荷载作用下,由负摩阻力引起的下拉荷载最大,可视为常规方法设计的上限值;桩顶荷载与堆载施工顺序对负摩阻力的影响也很大,先施加桩顶荷载,后进行堆载所产生的负摩阻力最大,反之最小。     

构造应力场边界载荷反演的有限元逆逼近法

构造应力场问题的边界条件、加载方式是地应力场分析的难题,由于构造运动未知和地质构造复杂,无法直接求解。边界力反演是进行构造应力场研究的重要手段。从弹性力学基本方程出发,针对油气储层地应力场反问题,结合对区域构造应力场认识,利用阻尼最小二乘法建立了应力场反演的约束模型,提出了一种用于反演构造应力场边界力的反分析方法,该方法根据部分实测点的应力值进行边界力反演,并与有限元方法和有限差分法相结合,以边界力为参量逼近目标函数。根据研究区少量地应力实测资料采用该方法进行最优化计算,使得计算应力场与实测应力达到最优拟合,以弥补目前的回归反演方法和边界载荷调整法的不足来提高计算精度。算例表明了该方法的有效性,可用于解决构造应力场边界力的反演问题。     

多向荷载作用下层状地基刚性桩筏基础分析

提出一种多向荷载作用下层状地基中刚性桩筏基础的计算方法。基于剪切位移法,采用传递矩阵形式分析了竖向荷载下桩顶面-桩顶面相互作用;引入修正桩侧地基模量,采用有限差分法分析了水平荷载下桩顶面-桩顶面相互作用;基于层状弹性半空间理论,分析了多向荷载下桩顶面-土表面、土表面-桩顶面、土表面-土表面的相互作用关系。建立了桩土体系柔度矩阵,得到了多向荷载下层状地基中刚性桩筏基础的受力和变形的关系以及桩的内力和变形沿桩身分布规律。通过与有限元对比,验证了该方法的合理性和修正地基模量的优越性,并对多向荷载作用下的桩筏基础进行了计算分析,计算结果表明,水平力将会引起桩筏基础的倾斜。     

Lateral force and centroid location caused by horizontal and vertical surcharge strip loads on a cross‐anisotropic backfill

This work presents analytical solutions for determining lateral force (force per unit length) and centroid location caused by horizontal and vertical surcharge surface loads acting on a cross‐anisotropic backfill. The surcharge loading types are point load, line load, uniform strip load, upward linear‐varying strip load, upward nonlinear‐varying strip load, downward linear‐varying strip load, and downward nonlinear‐varying strip load. The planes of cross‐anisotropy are assumed parallel to the backfill ground surface. The proposed solutions, derived by integrating the lateral stress solutions (Int. J. Numer. Anal. Meth. Geomech. 2005; 29 :1341–1361), do not exist in literature. Clearly, the type and degree of material anisotropy, loading distance from the retaining wall, and loading types markedly impact the proposed solutions. Two examples are utilized to illustrate the type and degree of soil anisotropy, and the loading types on the lateral force and centroid location in the isotropic/cross‐anisotropic backfills generated by the horizontal and vertical uniform, upward linear‐varying and upward nonlinear‐varying strip loads. The parametric study results demonstrate that the lateral force and centroid location accounting for soil anisotropy, loading distance from the retaining wall, dimension of the loading strip, and loading directions and types differ significantly from those estimated using existing isotropic solutions. The derived solutions can be added to other lateral pressures, such as earth pressure or water pressure, required for stability and structural analysis of a retaining wall. Additionally, they can simulate realistically actual surcharge loading problems in geotechnical engineering when backfill materials are cross‐anisotropic. Copyright © 2007 John Wiley & Sons, Ltd.     

RESPONSE OF A LAYERED ELASTIC MEDIUM TO A MOVING STRIP LOAD

Analytical determination of stresses and deformations caused by moving loads is vital to foundation and pavement designs. In current applications, moving loads are often approximated to be vertical impact loads. In this work, however, a live load is modelled as a uniform distribution of normal or shear stresses in actual motion. Then, a layer stiffness approach utilizing linear elasticity is followed in determining the surface and interior deformations due to the live load. By superimposing the two solutions for normal and shear surface stresses, the new approach can be made to provide an approximate solution to the problem of evaluating stresses and deformations caused by a wide wheel load rolling on a layered elastic system. Although elastic solutions in general are inadequate to explain the more significant consequences of pore pressure generation and dissipation in the soil subgrade, these results can certainly be useful to examine the shearing effects of wide rolling wheels on the asphalt layer and immediate settlement of the subgrade. It is found that the dynamic effects of a smoothly rolling wide load are significant at relatively low wheel velocities compared to those of shear waves in the subgrade and base.     

Lateral stress caused by horizontal and vertical surcharge strip loads on a cross‐anisotropic backfill

This study derives analytical solutions for estimating the lateral stress caused by horizontal and vertical surcharge strip loads resting on a cross‐anisotropic backfill. The following loading types are employed in this work: point load, line load, uniform strip load, upward linear‐varying strip load, upward nonlinear‐varying strip load, downward linear‐varying strip load and downward nonlinear‐varying strip load. The cross‐anisotropic planes are assumed to be parallel to the horizontal surface of the backfill. The solutions proposed herein have never been mentioned in previous literature, but can be derived by integrating the point load solution in a Cartesian co‐ordinate system for a cross‐anisotropic medium. The calculations by the presented solutions are quick and accurate since they are concise and systematized. Additionally, the proposed calculations demonstrate that the type and degree of material anisotropy and the horizontal/vertical loading types decisively influence the lateral stress. This investigation presents examples of the proposed horizontal and vertical strip loads acting on the surface of the isotropic and cross‐anisotropic backfills to elucidate their effects on the stress. The analytical results reveal that the stress distributions accounting for soil anisotropy and loading types are quite different from those computed from the available isotropic solutions. Restated, the derived solutions, as well as realistically simulating the actual surcharge loading circumstances, provide a good reference for the design of retaining structures for the backfill materials are cross‐anisotropic. Copyright © 2005 John Wiley & Sons, Ltd.     

动载确定方法对岩石动态断裂韧度测试的影响

为了考察不同方法确定的动态载荷对测试岩石动态断裂韧度的影响,在SHPB压杆系统上动态冲击直径80 mm的大理岩圆孔裂缝平台巴西圆盘,获得了弹性压杆上的应力波形,间接计算得到3种不同的作用在圆盘端部的动态载荷。将载荷输入ANSYS动态有限元模型中,求得了相应的动态应力强度因子,并根据试验-数值分析方法确定了岩石的动态断裂韧度测试值。结果表明,在加载速率约为4.0×104 MPa•m1/2/s的条件下,采用三波法确定的大理岩的平均动态断裂韧度为 3.92 MPa•m1/2,采用一波法比三波法计算的结果偏低11.22%,采用二波法比三波法计算的结果偏高20.15%,3种方法得到的结果差异较大。应力波在传播过程中,通过圆盘表面和预制裂缝面发生散射,部分能量不断发生释放是造成圆盘试件两端加载载荷不相等的主要原因。三波法是3种方法中比较理想的动态载荷确定方法,但需要考察试件的动态应力平衡性。     

The influence of load inclination on the undrained bearing capacity of strip footings on slopes

The effect of inclined loading on the bearing capacity of foundations on horizontal ground surface is well established and both the exact solution and simpler empirical equations are available for the calculation of the failure loads. However, for footings on or near slopes complete solutions are available only for vertical loading. This paper investigates the influence of inclined loading on the horizontal and vertical failure loads. The finite element, upper bound plasticity and stress field methods are used to examine a wide range of geometries and soil properties. The methods are first validated against known solutions for two special cases and are subsequently employed to investigate the effect of the geometrical and material properties on the failure loads and the bearing capacity load interaction diagram. Based on this investigation an empirical equation is proposed for the load interaction diagram for undrained inclined loading of footings on or near slopes.     

饱和半空间瞬态加载衬砌隧道的动力响应研究

工程中,地下衬砌隧道会遇到水压破裂压力、爆炸及突然开挖等瞬态荷载作用,将这些荷载理想化为作用在衬砌内边界上的均布瞬态荷载,研究圆柱形衬砌隧道在突加荷载、阶跃荷载和三角形脉冲荷载作用下的动力响应规律。根据Biot波动理论推导出半空间饱和介质的控制方程;视衬砌结构为弹性材料导出衬砌结构的控制方程。用极大半径凸圆弧近似半空间直边界,采用Graff加法公式进行坐标变换,将直角坐标表示的通解转化为极坐标表示的通解。根据边界条件导出衬砌和土体的位移、应力和孔隙压力的Laplace变换域的解答。利用反Laplace变换数值计算方法,将解答转换为时域解,得出3种瞬态荷载作用下衬砌隧道地面位移峰值、衬砌应力和孔隙压力的分布规律。     

Elastic solutions for stresses in a transversely isotropic half‐space subjected to three‐dimensional buried parabolic rectangular loads

In many areas of engineering practice, applied loads are not uniformly distributed but often concentrated towards the centre of a foundation. Thus, loads are more realistically depicted as distributed as linearly varying or as parabola of revolution. Solutions for stresses in a transversely isotropic half‐space caused by concave and convex parabolic loads that act on a rectangle have not been derived. This work proposes analytical solutions for stresses in a transversely isotropic half‐space, induced by three‐dimensional, buried, linearly varying/uniform/parabolic rectangular loads. Load types include an upwardly and a downwardly linearly varying load, a uniform load, a concave and a convex parabolic load, all distributed over a rectangular area. These solutions are obtained by integrating the point load solutions in a Cartesian co‐ordinate system for a transversely isotropic half‐space. The buried depth, the dimensions of the loaded area, the type and degree of material anisotropy and the loading type for transversely isotropic half‐spaces influence the proposed solutions. An illustrative example is presented to elucidate the effect of the dimensions of the loaded area, the type and degree of rock anisotropy, and the type of loading on the vertical stress in the isotropic/transversely isotropic rocks subjected to a linearly varying/uniform/parabolic rectangular load. Copyright © 2002 John Wiley & Sons, Ltd.