偏心荷重下转动浅基承载力的严密解法

根据偏心荷载作用下地基失稳的特征,构建了偏心荷重下地基极限状态时土体中不连续应力场和速度场,分析了土体中不连续应力场和速度场的特点。基于上、下限原理的推论,建立计算模型,构建求解地基极限承载力的算法,通过解三类边值问题,求得偏心荷重下转动浅基极限承载力的严密解,该解满足刚性基础转动的速度边界和其他已知的边界条件。算例表明解垂直偏心荷载下地基极限平衡问题的解法是可行的,且结论可靠。     

条形荷载下不排水土坡破坏模式判定及极限承载力估算

为了确定边坡放置地基的极限承载力,采用弹-理想塑性有限元,分析了不排水土坡(u=0)的坡肩上放置条形基础的地基极限承载力。结果表明,不排水土坡地基极限承载力仍可沿用Terzaghi水平地基极限承载力计算公式qu=cuNc。不考虑重度条件下,土坡地基的破坏模式类似于水平地基破坏模式,承载力系数Nc与极限分析法计算的结果一致;考虑重度条件下,可用土坡实际状态的Taylor数与其极限状态的Taylor数的差值N作为判别指标,当该值较小时将发生边坡破坏模式,较大时发生地基破坏模式,找出了两种破坏模式的N分界点,确定了Nc与N的关系,建立起坡肩作用条形荷载下不排水土坡极限承载力计算公式。     

基于二阶锥规划有限元增量加载法的条形浅基础极限承载力分析

地基极限承载力分析是土力学研究中的一个经典课题。基于Hellinger-Reissner混合变分原理和有限元方法,将岩土体弹塑性问题构造成基于有限元框架的二阶锥规划（second-order cone programming,SOCP）问题,进而提出一种基于二阶锥规划理论的增量有限元法,即FEM-SOCP法。将岩土体弹塑性问题构造成二阶锥规划的数学优化问题,可以避免采用传统弹塑性计算中复杂的应力点积分等算法和屈服面棱角的平滑处理。此外,对于二阶锥规划问题,可以采用具有原始?对偶内点求解法的标准数学规划求解器MOSEK进行求解。将增量加载FEM-SOCP法应用于经典的基底粗糙的条形浅基础地基极限承载力分析中,分别考虑了关联和非关联塑性条件下的Mohr-Coulomb屈服准则。数值结果表明：所提出的增量加载FEM-SOCP法获得的地基承载力系数及地基承载力与传统FEM计算结果基本一致,而与常规有限元计算结果相比,基于增量加载的FEM-SOCP法所获得的屈服区更加平滑。     

含下卧软层地基极限分析数值试验及承载力安全性研究

含下卧软层的地基在工程中比较常见,但目前对其破坏机制和承载力计算还依赖于很多的假设条件和经验。基于极限分析和快速拉格朗日数值模拟方法,对含下卧层地基的加载破坏过程进行应变控制的数值试验。首先用Prandtl的极限承载力解析解验证了数值试验的有效性,然后对一个含下卧层的地基模型按照不同加载尺寸进行加载数值试验,发现模拟结果存在明显的荷载尺寸效应。模拟获得不同荷载条件下的极限承载力以及相应的应力-应变曲线、渐进破坏模式、下卧软层的沉降规律。进一步将下卧层的地基承载力加载安全系数与强度折减系数进行对比,发现折减系数与加载系数随荷载尺寸变化的趋势相同,但前者的数值偏小。     

关于偏心荷载下地基承载力的验算方法

针对偏心荷载作用下的地基承载力验算,以条形基础为例,介绍了Meyerhof提出的考虑荷载偏心的地基承载力计算方法,并通过与有限元上下限极限分析这一极限荷载精细数值计算方法的比较验证了其正确性和准确性。所考虑的地基既包括排水条件下的一般地基,也包括不排水强度随深度近似线性增大的饱和黏性土地基。依据有关机理深入分析了基础埋深对偏心荷载下地基承载力计算精度的影响,明确了可以考虑基础埋深有利作用的偏心距上限。进而采用本课题组所提出计算地基承载力的改进公式,通过计算并结合理论分析比较了国内规范验算方法与按Meyerhof方法的验算,表明采用后一方法更为合理。     

倾斜荷载作用下层状非均质地基的极限承载力

通过一定的室内模型试验对倾斜荷载作用下层状非均匀地基的极限承载力进行了试验研究.引入荷载倾斜影响系数,对于层状非均质地基的极限承载力计算公式进行了修正,给出了倾斜荷载作用下的地基极限承载力计算公式.根据模型试验结果,对所建立的计算公式进行验证.结果表明,对于倾斜荷载作用下的层状非均匀地基,根据修正后的Meyerhof承载力公式所得到的计算结果能够较好地与试验结果吻合.     

各向异性损伤对地基承载力影响分析

应用耦合损伤的滑移线法,求解平面应变条件下的极限荷载问题,推导出考虑各向异性损伤的平面应变问题应力场的滑移线解,针对平顶钝角楔体极限荷载的滑移线场进行了详细推导,绘出滑移线场图,推导出极限荷载的公式,分析了各向异性损伤对浅基础极限荷载的影响,即损伤对地基承载力系数Nq影响很大;损伤越大,地基承载力系数Nq值越小,且随着D1,D2的增加,Nq越接近最小值。     

压实填土地基承载力深度修正系数研究

软弱土层上的压实填土地基作为一种非均质地基广泛应用,但目前关于这种地基的破坏模式和承载力机制的研究还不够深入,工程中仍采用现行规范深度修正系数取值为2的方式进行设计。采用有限差分方法建立数值模型,研究了影响压实填土地基承载力深度修正系数的因素,并与现行规范推荐方法的取值进行对比。结果表明：填土越深、基础外超载越大,参与地基剪切破坏运动的土体体积越大,地基极限承载力越高;而随着填土宽度的增大,地基极限承载力先提高,后趋于稳定;砂土内摩擦角的增大也会提高地基极限承载力。总体来说,规范中的推荐值是偏于保守的,会造成承载力的浪费;而在下卧层强度很小、填土深度较大时该推荐值可能高估压实填土地基承载力,存在偏于不安全的可能,在工程设计时应予以考虑。     

极限分析有限元法讲座Ⅲ--增量加载有限元法求解地基极限承载力

利用有限元法,通过增量加载的方式来求解地基极限承载力。随着荷载的逐渐增加,地基由初始的线性弹性状态逐渐过渡到塑性流动的极限破坏状态,此时有限元的计算将不收敛。它不但可以获得地基的极限荷载的值及荷载-位移关系,而且还能得到经典极限分析法中所采用到的破坏机构。当采用关联流动法则或采用剪胀角为 φ/2 的非关联流动法则时,获得的破坏机构与 Prandtl 的破坏机构一样。对 Prandtl 解的经典算例进行了分析,结果表明:屈服准则的选用对计算结果的影响很大,选择与实际问题相匹配的屈服准则方能得到比较精确的结果。在求解平面应变问题时,在关联流动法则条件下,采用Mohr-Coulomb 内切圆屈服准则,或在非关联流动法则下采用 Mohr-Coulomb 匹配 DP 准则所得结果与 Prandtl 精确解极为接近,可望应用于实际工程分析中。     

空洞上方浅基础地基破坏模式与极限承载力分析

利用极限分析上限法求解地基极限承载力问题的关键在于构造合适的破坏模式。当地基下方存在空洞时,地基的破坏模式变得相当复杂。通过分析空洞存在时地基的受力特点及破坏形态,将地基破坏范围划分成为不同的刚性区和过渡区,构造了空洞上方条形基础地基的破坏模式。利用上限法,建立与破坏模式对应的速度场,推导了破坏模式不同区域内的耗散功率和外力功率,得到地基极限承载力的目标函数,并采用数学优化方法进行求解,获得了极限承载力的上限解。通过算例分析,讨论了空洞顶板厚度、空洞大小与地基极限承载力的关系,并与无空洞条件下地基极限承载力进行对比分析。结果表明,随着空洞顶板厚度增加,地基极限承载力增加,破坏模式也由地基与空洞之间扩散到地基两侧;空洞顶板厚度存在临界值,当超过此临界值时,空洞对地基极限承载力的影响可忽略。     

复合加载情况下双层地基极限承载力研究

在复合加载情况下精确求解层状非均质地基的极限承载力,具有很强的工程实用与理论参考价值。基于土体极限平衡理论与通用有限元软件ABAQUS,针对复合加载情况下上硬下软的双层不排水饱和软黏土地基的极限承载力,进行了大量的数值计算,得出了上层土临界深度Hcr的计算公式、竖向极限承载力Pv的计算公式以及复合加载情况下地基破坏时的破坏包络面方程。研究结果表明：上层土临界深度Hcr取决于土层间强度比Su1/Su2;竖向极限承载力Pv与破坏包络面取决于土层间强度比Su1/Su2、上层土深度H1与基础型式。     

水平受荷桩的弹性有限元虚拟加载上限分析

基于将塑性上限分析等效为弹性迭代计算的总量虚拟加载上限分析理论,在商业化有限元软件ABAQUS中实现了弹性有限元虚拟加载上限方法（弹性有限元T-EMSD）。应用弹性有限元T-EMSD法分析了不排水黏土中的二维水平受荷桩,其获得的荷载-位移曲线与弹塑性有限元分析结果一致,其极限承载力与塑性解相近。在极限位移加载量下弹性有限元T-EMSD法对应的上限机构从弹性始速度场开始随迭代逐渐演化,迭代收敛后的速度场和解析塑性破坏机构相似。与其他基于可变强度概念（MSD）的方法相比,弹性有限元T-EMSD法对水平受荷桩桩身的分析具有更高的精度。弹性有限元T-EMSD法最大的优势在于可在计算中自然地获得塑性机构,因而可被用于研究一些塑性机构难以构造的复杂问题,并对弹塑性数值方法进行验证。     

Upper-bound multi-rigid-block solutions for bearing capacity of two-layered soils

The limit analysis method has been widely used in the stability analysis of geotechnical problems including the bearing capacity of foundations. Two main approaches have been followed in the limit analysis to improve the calculation of the bearing capacity of foundations. One approach is to combine limit analysis with the finite element method and linear/nonlinear programming. The other is to use a multi-rigid-block mechanism to obtain an upper-bound solution. In this paper, the multi-rigid-block upper-bound method with a modified failure mechanism from that of Florkiewicz [Florkiewicz A. Upper bound to bearing capacity of layered soils. Can Geotech J 1989;26(4):730–6.] was employed to calculate the bearing capacity of foundations. Attention was paid particularly to the bearing capacity of strip footings over a two-layered soil. In order to verify the effectiveness of the modified mechanism, comparisons were made with other well-known solutions. The results showed improvements over the best available multi-rigid-block upper-bound solutions given by Michalowski and Shi [Michalowski RL, Shi L. Bearing capacity of footings over two-layer foundation soils. J Geotech Eng ASCE 1995;121(5):421–8.], and fair consistence with the results from the finite element limit analysis in Shiau et al. [Shiau JS, Lyamin AV, Sloan SW. Bearing capacity of a sand layer on clay by finite element limit analysis. Can Geotech J 2003;40(5):900–15.].     

Undrained lower bound solutions for end bearing capacity of shallow circular piles in non-homogeneous and anisotropic clays

The undrained bearing capacity of shallow circular piles in non-homogeneous and anisotropic clay is investigated by the lower bound (LB) finite element limit analysis (FELA) under two-dimensional (2D) axisymmetric condition using second-order cone programming, and the new solution of the problem is presented. Modified from the isotropic von Mises yield criterion, a cross-anisotropic undrained strength criterion of clays under the axisymmetric state of stress requiring three input shear strengths in triaxial compression, direct simple shear, and triaxial extension is employed in the 2D axisymmetric LB FELA. Parametric studies on the effects of pile embedment ratio, dimensionless strength gradient, anisotropic strength ratio, and pile roughness are investigated extensively, while the predicted failure mechanisms associated with these parameters are discussed and compared. Numerical results of undrained end bearing capacity of shallow circular piles are summarized in the form of design tables that are useful for design practice and represent a new contribution to the field of pile capacity considering the combined effects of undrained strength non-homogeneity and anisotropy.     

Stability of unsupported conical excavations in non-homogeneous clays

New plasticity solutions for the undrained stability of unsupported conical excavations in homogeneous and non-homogeneous clays were solved by axisymmetric finite element limit analysis. Three parametric studies were performed on excavated height ratios, slope inclinations and dimensionless strength gradients. In all cases, the exact stability factors were accurately bracketed by computed bound solutions within 0.6%. An accurate closed-form equation of the stability factor was proposed from nonlinear regression analysis of lower bound solutions. New conical stability factors for soil cohesion, strength gradient, and coupling effect of these components were deduced to conveniently and accurately predict a safe solution in practice.     

三维塑性极限分析下限法原理及应用

介绍了三维塑性极限分析下限法的基本原理。将计算区域离散为有限个三维空间单元,利用极限分析中的下限定理,借助有限元思想建立数学规划模型,并引入非线性规划的数值方法寻求问题的下限解。以均匀土质边坡的极限荷载作为算例,比较了所得数值解与经典塑性力学理论解,论证了该方法的正确性。     

Undrained stability of wide rectangular tunnels

Upper bound rigid block methods and finite element limit analysis are applied to investigate the undrained stability of rectangular tunnels. Solutions for a range of soil properties and tunnel geometries are obtained for situations in which the tunnel is wider than it is high. The upper and lower solutions for tunnel stability, which accurately bracket the true collapse load, are summarised in a series of stability charts for use in tunnel design.     

Two- and three-dimensional bearing-capacity solutions for footings on two-layered clays

In this paper, finite element analysis is used to predict the undrained bearing capacity of strip, square and circular footings resting on layered clays. The soil profile consists of two clay layers with different thicknesses and properties. The results are compared with previous solutions for strip footings on layered clays. The bearing-capacity behaviour is discussed and the bearing-capacity factors are given for various cases involving a range of layer thicknesses and properties of the two clay soil layers.     

各向异性和非均质黏土粗糙地基承载力上限计算

通过构建一种网格状多刚性块破坏机制,利用上限法研究了黏土地基各向异性和非均质对粗糙条形基础地基承载力的影响。该破坏机制允许滑动面和速度矢量沿破坏区域过渡区的径向和切向发生变化,形成了更为精确的破坏机制和塑性流动速度场。根据上限定理得到各向异性和非均质性黏土上条形基础极限承载力的目标函数,将其转化为非线性规划问题,并通过编程对该计算模型进行求解。与已有研究成果对比分析表明,两者具有较好一致性,且优于现有同类方法。文中亦探讨了黏土地基非均质、各向异性以及其他相关参数对地基承载力的影响。