边坡稳定分析刚体有限元上限法的锥规划模型

极限分析方法是土边坡稳定性分析的重要方法之一。刚体有限元上限法是其中的一类,此类方法仍旧存在一些关键问题需要完善。由于单元的刚性假设,系统的塑性变形内能耗散仅发生在单元间的界面上,故此类方法的性能主要取决于界面的布局,即采用非结构化三角形单元计算往往精度较差。为此,提出了基于滑动面摄动的刚体有限元上限法及临界滑动面的搜索方法。首先,在考虑刚体转动的基础上构造刚体有限元上限法的二阶锥规划模型,用于确定在给定试滑动面条件下的运动许可速度场。其次,将试滑动面的控制参数视为决策变量,建立搜索临界滑动面的非线性非凸优化问题模型,并采用非线性单纯形方法和粒子群方法求解此优化问题找出临界滑动面。通过经典边坡稳定问题的分析求解,验证了所提出的新方法,进一步证实了网格类型（即界面的布局）是影响刚体有限元上限法计算精度的主要因素。经过计算结果的对比发现,在刚体有限元上限法中考虑刚体转动是非常必要的,不仅可以提高刚体有限元上限法的计算精度,还可以克服此方法对界面布局的依赖性。     

填筑和开挖边坡的稳定性分析

在有限元应力应变分析的基础上,采用以剪应力定义的安全系数和模式搜索的方法确定边坡的安全系数及最危险滑动面。在平面应变条件下,对土层均质的填筑和开挖边坡的稳定性进行了较为全面的研究,并同天然边坡的安全系数和破坏面进行比较,计算结果表明,这3类边坡滑动面的形状基本一致,天然边坡安全系数要高于填筑和开挖边坡;同时比较了采用有限元稳定评价方法和采用极限平衡方法计算结果的差别;最后分析研究了计算参数变化对填筑和开挖边坡稳定性的影响。     

深层高压注浆加固古滑坡滑动带试验及效果分析

通过深层高压注浆法加固混碎石黏土古滑坡滑动带的成功工程实例,结合各种检测手段及有限元弹塑性分析,阐述对混碎石黏土古滑坡滑带土实施深层高压注浆加固的机制;测试和分析了滑动带压力注浆前后物理力学参数的变化特征;计算了注浆前后滑坡的剩余下滑力及安全系数。电镜观测表明,水泥浆液对于滑动带岩土体的充填置换率可达到5%～10%;土工试验结果显示,滑动带土体黏聚力由20.20 kPa提高到39.54 kPa,内摩擦角由16°升高到21.30°;声波CT扫描试验表明,滑动带土体的声波波速增长约为31%。有限元弹塑性分析计算表明,滑坡由注浆前的不稳定状态变为注浆后的基本稳定状态,边坡剩余下滑力大幅减小,为支挡结构设计提供了优化依据。试验及计算分析结果表明,高压注浆法用于加固滑动带岩土体是可行的,经注浆加固后的滑动带岩土体的力学性能得到了明显改善,边坡安全性得到了提高。     

蒙特卡洛法与有限元结合搜索边坡临界滑动面

把蒙特卡洛随机搜索法与有限元法相结合,搜索边坡的临界滑动面及其对应的最小安全系数.通过随机跳跃法生成一定数量的初始试算滑动面,根据有限元分析的应力和孔隙水压力结果计算给定滑动面的安全系数,利用随机走步法不断更新试算滑动面,使试算滑动面安全系数不断减小,最终确定边坡临界滑动面及其对应的最小安全系数.本方法通过随机搜索,克服了大多数常规优化方法易陷入安全系数局部极小的问题.数值算例分析说明了本文提出的确定边坡临界滑动面的方法的有效性和优越性.     

立节北山滑坡地震稳定性及破坏过程三维有限元分析

针对甘肃省舟曲县立节镇北山滑坡在地震情况下的稳定性及失稳后可能威胁的范围问题，笔者使用显式有限元及动力学大变形方法分析了该滑坡在地震情况下的稳定性及失稳后的破坏过程。首先基于滑坡纵剖面建立二维平面应变模型分析滑坡稳定性及破坏过程，之后使用航测点云数据建立三维滑坡模型进行计算分析。地震条件下三维计算确定的滑坡安全系数结果低于1.0，较二维计算结果低7.7%；三维情况计算出的滑动距离达1.73 km，高于二维情况17%，说明滑坡灾害问题使用二维小变形方法得出结果偏于危险。笔者采用显式有限元、三维及大变形计算方法分析更符合实际，可为滑坡灾害的工程防治提供科学理论依据。     

用有限元极限平衡法分析边坡的稳定性

用有限元分析为基础的极限平衡法计算了边坡安全系数,并将结果与摩尔-库仑屈服准则、外接圆D-P屈服准则强度折减法得到的结果、以及Spencer（极限平衡）法的结果进行了比较,结果表明,有限元极限平衡法计算的边坡安全系数与摩尔-库仑屈服准则强度折减法和Spencer法的结果很一致,而D-P准则强度折减法得到的结果偏差较大。叙述了存在已知滑动面和不存在滑动面时有限元极限平衡法的计算方法。用大型通用有限元软件,如ANSYS等,可方便地进行这种方法的边坡安全系数的计算,并有足够的精度。     

修正有限元强度折减法与失稳判据在边坡稳定分析中的应用

现有的有限元强度折减法用于边坡稳定性分析,在土体变形特性的认识和失稳判据上有缺陷。基于土体的应力-应变曲线和邓肯-张非线性弹性模型进行有限元强度折减法计算分析,在现有方法仅对强度参数--黏聚力c和内摩擦角φ进行折减的基础上,提出了切线弹性模量Et的折减算法,更全面地体现折减的内涵;提出了以潜在滑动区域的能量积分变化作为判据,根据能量积分与折减系数关系曲线的特征,将其拟合成两段直线,以其交点所对应的折减系数作为边坡的安全系数。比较分析了修正Duncan有限元强度折减法、失稳判据与文献的计算结果,同时分析了边坡滑动区能量积分计算范围对结果的影响,并用于南水北调东线实际工程计算分析。结果表明：修正Duncan有限元强度折减法能较全面考虑土体的变形和强度特性,新的边坡滑动失稳判据可行且具有较好的计算稳定性,能保证工程建设的安全。     

基于数值应力场的边坡临界滑动场

数值分析方法（有限元、有限差分等）可计算边坡应力、变形和破坏,但直接计算边坡安全系数比较困难;而极限平衡法需对边坡体内力作出假设,计算精度受到限制。利用数值分析方法计算得到的应力场,重新定义推力最大原理中推力的含义,将常规边坡临界滑动场数值模拟方法作进一步的改进,提出基于数值应力场的边坡临界滑动场方法。该方法可充分发挥两大方法的优势,能够准确、方便地搜索出边坡临界滑动面位置,并给出更合理的边坡安全系数。通过两个算例的稳定性分析,并与其他方法进行比较,验证了该方法的准确性和合理性。     

乌江沙陀水电站厂房坝段深层抗滑稳定分析

通过有限元强度折减法计算了乌江沙陀水电站厂房坝段的深层抗滑稳定安全系数,与基于遗传算法的最危险滑动面搜索方法进行了比较,结果表明,两种方法得到的安全系数是接近的,滑动面的位置也基本一致。     

圆弧滑动有限元土坡稳定分析

提出一种基于有限元应力变形计算的边坡稳定分析方法,仍假定滑动面形状为圆弧形,有限元网格由一组同心圆作为纬线,一组竖向线为经线构成。对两相邻圆弧线所夹的弧形带分析滑动力和抗滑力,建立平衡方程,确定安全系数,其中小值安全系数对应的弧形带为可能的滑动带。变化圆心位置用优选方法寻找最小安全系数对应的圆心,从而得出真正的滑动面。算例分析表明,该方法计算所得安全系数与Bishop法接近,是合理的,其突出优点是由有限元计算直接得出滑面上的应力,而不须作近似的插值处理,因而应力更准确。该方法可以考虑土的非线性变形特性,也更符合土的实际情况。此外,用有限元计算得出位移,亦可将稳定分析和变形联系起来,为现场通过位移监测来估计边坡的稳定性提供了可能性,同时也为膨胀土边坡的稳定分析中考虑膨胀性的影响提供了可能性。     

An axisymmetric BNEF method of circular excavations taking into account soil-structure interactions

An axisymmetric BNEF method was proposed to analyze circular excavations. The governing equations of embedded cylindrical diaphragm walls were derived from thin-shell theory. A hyperbolic model was proposed to account for soil-structure interactions. Then, a FEM program was coded and verified using a case study. Overall, the method caught the main trend in the measured results and provided a reasonable estimation for engineering practice. Furthermore, the development of earth pressure and the hoop stress in the walls (ring beams) can be obtained using conventional parameters, extending the ability of the classic BEF method to aid the design of circular excavations.     

Elasto-plastic analysis of an underground opening by FEM and coupled FEBEM

The elasto-plastic analysis of a circular underground opening using the finite element method (FEM), the coupled finite element and the boundary element method (FEBEM) is presented. The coupling procedure and elasto-plastic formulation are discussed. The effect of in situ stress ratio on yielded zone, displaced shape and principal stresses is presented. The results of FEBEM analysis are compared with those obtained from FEM. The computation time and the number of iterations as required by FEBEM and FEM are compared.     

可调圆弧木模板的三维有限元计算分析

自 20 世纪 90 年代我国开始使用可调圆弧木模板以来,定型圆弧钢模板在建筑中的地位被逐步取代。由于传统的以梁代替面进行的计算,存在着过度简化计算模型的缺陷,造成计算结果与实际相差较大。所以,采用三维有限元对模板 进行计算分析,得到的计算结果经检验与实际能吻合的较好。     

顶管工程工作井位移与土抗力分析

为保证顶管工程的施工安全和经济设计,必须深入分析顶力作用下工作井及周围土体的应力、位移特性。以两个实际顶管工程为工程背景,针对工作井的浅埋、深埋圆形沉井,采用三维有限元分析,给出了浅埋沉井土抗力沿圆周分布的拟合方程。分别采用《规程》[1]和《手册》[2]推荐的计算方法和三维有限元分析,对顶力作用下深埋、浅埋工作井的位移和新增土抗力进行对比分析,结果表明：①由于只考虑了顶力后背一侧半圆范围内土体抗力的作用,规范法和手册法将导致土抗力计算结果偏大;②顶力作用位置对深埋沉井的井壁变位、土抗力大小和分布情况影响显著;③规范法和手册法仅适用于顶力作用于沉井底部的浅埋沉井。     

基于连续介质模型的隧道衬砌结构可靠性研究

基于连续介质力学模型,提出了一种简单易行的可靠性分析方法,该方法能利用现行通用的有限元程序,能得到显式的功能函数,但又不同于响应面法。以公路隧道圆形隧洞的开挖支护为例,采用此方法将围岩参数c (黏聚力)、E,? 和二次衬砌厚度D作为随机变量,采用有限元程序得到单一变量变化下的衬砌轴应力值,进而拟合出衬砌轴应力关于单一随机变量的关系式,并通过多元线性或非线性回归得到衬砌轴应力关于所有随机变量的显示表达式,再对圆形隧洞的衬砌结构进行了敏感性和可靠性分析,得到了一些有意义的结论。     

轴对称饱和土体结构在多向荷载下的液固耦合问题

以轴对称系统液固二相耦合饱和介质为研究对象,在假定所有物理量沿径向和轴向可以变化但沿环向不变化、变形满足小变形、孔压流满足达西定律及有效应力原理成立的前提下,建立相应的增量形式的有限元分析方法。以一饱和轴对称试样为例,应用修正剑桥弹塑性模型,进行了应变控制式常规三轴固结不排水剪切试验过程的有限元模拟,以表明所建立方法的可行性。     

深层搅拌桩复合地基沉降计算理论研究

为了寻找深层搅拌桩复合地基沉降计算的简便方法，根据已有的研究成果，进行了桩与桩周土间摩擦力分布规律的假设和简化，运用单位元法推导了深层搅拌桩复合地基沉降计算公式和设计用曲线，同时采用轴对称有限元法，对一实际工程的实测成果进行了分析。计算结果表明：文中理论计算的中级荷载下的沉降量值介于分层总和法与有限元法计算值之间，但文中方法比有限元计算法简单，比分层总和法更合理。     

混合法计算三维非均匀介质中的电磁场

介绍了联合运用积分方程法与有限元法(简称混合法)来计算三维非均匀介质中电磁场分布的理论方法,并进行了数值模拟实验。混合法的原理是引入一个包围非均匀目标体的虚构边界,在边界内部的场用有限元法模拟,在边界外部(包括边界)的场用积分方程表达,二者在边界上通过场的连续性耦合起来。数值实验结果表明,混合法既能显著地减小网格规模,又能灵活地模拟复杂的介质情况,且计算精度较高。     

GIS在边坡有限元数值分析前处理中的应用

有限元数值分析是边坡稳定性研究中较为常用的定量分析方法,但有限元分析的前处理工作通常比较繁琐,从而影响了该方法的整体工作效率。利用GIS软件实现了从计算剖面确定、计算剖面自动绘制到网格自动剖分的全部前处理过程。有效提高了边坡有限元计算的前处理工作效率,同时也提高了计算剖面图的科学性和边坡评价的准确性。     

边坡变形破坏过程的大变形有限元分析

在用有限元法分析边坡稳定性时,引入计算大变形问题的更新的拉格朗日方法,推导了边坡大变形弹塑性有限元分析的方程式。采用边坡某一幅值的等效塑性剪应变区,从坡脚到坡顶贯通前的折减系数作为边坡安全系数。在此基础上,采用弹塑性大变形有限元分析软件计算了均质土坡不同坡角的安全系数,将其与小变形分析的结果进行了对比分析,结果表明：用弹塑性大变形有限元分析边坡失稳破坏的过程中,既考虑了岩土材料的非线性,又考虑了边坡的几何非线性,使计算结果更趋合理。并结合东深供水改造工程BIII2边坡进行了大变形有限元分析,计算结果与勘查到的实际边坡的滑动面分布位置比较接近。研究表明：该方法尤其适宜于软土类边坡或基坑的稳定性分析。