阻滑桩加固土坡稳定性分析的上限解法

从极限分析机动学方法出发,利用土的抗剪强度折减系数概念,建立了土坡的极限平衡状态方程,由此确定土坡的临界稳定安全系数及其相应的潜在破坏模式.对于典型问题,通过与现有极限平衡解和有限元数值解的对比分析,验证了这种上限解法的合理性.进而对于在给定的荷载条件下不能满足抗滑稳定性要求的土坡,考虑采用阻滑桩加固方式,根据桩侧有效土压力的合理分布模式确定桩体与滑动面相交的截面上等效抗滑力和抗滑力矩,利用极限分析上限定理建立了阻滑桩加固土坡的极限平衡状态方程,将桩侧土压力作为目标函数,运用数学规划方法确定了极限平衡状态时的临界桩侧土压力,以此为土坡加固中阻滑桩设计提供依据.通过数值计算与分析探讨了阻滑桩加固位置的优化布置等问题.     

地震效应下三维非均质土坡稳定性极限分析

基于极限分析上限定理,采用拟动力法研究了地震效应下非均质土坡的三维稳定性,并通过重度加大法（GIM）推导出边坡安全系数的显式表达式。此外,采用遗传算法进行优化,将所得结果进行对比验证,详细分析了地震条件下相关参数对边坡稳定性的影响。分析结果表明：边坡高度一定时,宽高比、边坡倾角的增大以及土体内摩擦角、非均质系数的减小均会造成安全系数的降低;拟静力法获得的结果相比于拟动力法较大,两种方法结果的差值会随着水平地震系数、有效内摩擦角的增大而增大,随着边坡倾角的增大而减小;土体放大系数的增大会导致安全系数减小,而地震波周期及波速的变化几乎对安全系数没有影响;滑动面轨迹受水平地震系数的影响较大,而受非均质系数的影响相对较小。     

非均质与各向异性黏土基坑抗隆起稳定分析

经典的Prandtl 机构不能很好地反映支护结构入土深度的影响,也不能反映土体强度非均质和各向异性的影响。根据塑性极限分析上限定理,基于修正的Prandtl机构,推导了既可以合理反映支护结构入土深度影响,又可以反映土体强度非均质和各向异性特性的基坑抗隆起稳定的上限解公式,同时基于商业化程序的二次开发,建立了可以考虑土体强度非均质和各向异性影响的强度折减弹塑性有限元方法。通过提出的修正Prandtl机构上限解与强度折减有限元方法以及多块体上限解的对比,验证了修正机构上限解以及强度折减有限元法的合理性和可行性,最后讨论了非均质系数、各向异性系数、支护机构入土深度等因素对基坑抗隆起稳定性的影响,并通过了实例的验证,相关结论可供工程参考。     

基于变分法的均质土坡稳定性分析

在平面应变条件下,从考虑强度折减的滑动土体极限平衡方程出发,推导出一般情况下均质土坡稳定分析的泛函极值等周模型。依据泛函取极值时必须满足的欧拉方程,通过坐标变换对变分问题求解,得到滑裂面函数及其上分布的法向应力函数。对于给定的边界条件及横截条件,土坡稳定问题转化为求解以两个拉格朗日乘子、安全系数和滑裂面一个端点的x坐标为4个未知量的函数极值问题,可以通过数值计算得到最小安全系数及所对应的最危险滑裂面。考虑坡顶拉裂缝,通过算例分析了不同的滑裂面位置及对应的安全系数。通过对ACADS边坡考题算例分析,在不考虑坡顶开裂的条件下,采用该方法的分析结果与依据极限平衡法得到的各裁判答案非常接近,在考虑坡顶开裂时,安全系数略有减小。     

非均质边坡稳定性上限分析评价研究

边坡稳定性评价研究是土力学经典问题。为研究非均质边坡稳定性问题,以上限分析法为理论基础,建立非均质边坡分段对数螺旋线破坏机构。应用拟合多项式近似描述对数螺旋线,计算破坏机构外力功率和内能耗散率,应用设计的计算流程进行强度折减循环计算,使边坡逐步过渡至极限平衡状态,确定边坡安全系数。以某露天矿南帮边坡为工程背景,经过反复强度折减,最终确定安全系数,通过变换初值的方式验证计算流程的稳定性,应用严格Morgenstern-Price法验证计算结果的准确性。研究结果表明：该流程性能稳定,不受初始值影响。相比于严格Morgenstern-Price法,二者安全系数计算结果偏差小于5%,可充分满足工程实践要求。同时,上限分析法求得的最危险滑面满足速度分离关系,具有更强的工程参考价值。     

考虑土体非均质和各向异性的锚索极限抗拔力研究

考虑土体的非均质性与各向异性,结合预应力锚索锚固段土体的破裂形式,利用非线性强度准则,根据极限状态下力的平衡原理,推导出一个可考虑土体非均质性、各向异性与非线性破坏性质等多种因素影响的锚索极限抗拔力计算公式,并通过理论计算,分析土体非均质性、各向异性以及非线性强度系数对锚索极限抗拔力的影响。计算研究表明,土体的非均质性与各向异性对锚索抗拔力的影响显著,随着土体非均质常数的增加,锚索极限抗拔力不断提高,而随着各向异性系数与非线性强度系数的增加,抗拔力则有所下降;当土体非均质常数与各向异性系数提高时,土体破裂范围不断增大,研究结果可为锚索锚固体的设计提供一定参考。     

考虑滑体法向应力突变的阻滑桩成拱效应研究

依据阻滑桩的受力特点,研究了桩间土拱效应的传力机理,分析指出桩间土拱效应是法向应力和拱脚剪切力共同作用而产生的土颗粒挤压变形现象,提出采用突变效应替代常用的桩-土承载比衡量桩间土拱效应的强度及作用范围更为直接可靠.模型桩推力试验结果和三维有限元模型计算表明,桩间水平土拱具一定的纵向效应,据此对比分析了阻滑桩三维与二维计算结果之间的异同,指出桩间土拱效应不仅受桩-土相对位移的影响,同时也受滑体自重应力和滑面摩擦效应的制约.运用二维模型重点研究了桩-土相对位移、桩间净距、推力荷载、土体抗剪强度对桩间土拱效应的影响,将所得规律与已有文献作了进一步地对比,继而总结了上述因素对拱效应的敏感度.     

地震效应和坡顶超载对均质土坡稳定性影响的拟静力分析

基于强度折减技术和极限分析上限定理,假定机动容许的速度场破坏面,考虑坡顶超载、水平和竖向地震效应影响推导了边坡稳定性安全系数的计算表达式。采用序列二次规划迭代方法（和内点迭代方法）对边坡安全系数目标函数进行能量耗散最小化意义上的优化计算,与多个算例的对比验证了其方法和程序计算的正确性;对影响土质边坡动态稳定性的一些因素进行了参数分析,分析表明：随着边坡倾角?、坡顶超载q、水平和竖向地震效应影响系数的增大,边坡稳定性安全系数显著下降;随着坡顶超载q、水平地震效应影响系数kh的增大、竖向地震效应影响系数kv的减小,边坡的潜在滑动面越来越深,潜在破坏范围越来越大。竖向地震效应对边坡稳定性也有一定影响,强震条件下的设计计算必须考虑竖向地震效应的影响。     

非均质各向异性软基上管桩基础承载力分析

对于非均质各向异性软土地基上管桩基础的承载力,基于塑性极限分析理论考虑软土强度的非均质与各向异性效应,提出了一种改进的三维极限分析上限解法。利用有限元分析软件ABAQUS,采用Hill理想弹塑性模型,对软基上管桩水平承载力进行了数值模拟,由此所得到的破坏模式与极限分析方法中的假设失稳机制基本一致。通过变动参数对比分析,探讨了软黏土不排水抗剪强度的各向异性和非均质性对于管桩基础承载力的影响。     

考虑强度各向异性的边坡稳定有限元分析

天然沉积的土层总是表现出一定程度的强度各向异性,但现有的边坡稳定有限元分析极少考虑各向异性的影响。对大型有限元软件ABAQUS进行二次开发,使其能够考虑土体黏聚力随大主应力方向的变化,动态更新增量迭代求解过程中边坡不同位置处的抗剪强度,进而提出具备安全系数自动搜索功能的各向异性边坡稳定分析方法。计算结果表明,均匀边坡的有限元解与极限分析上限解相差很小。如果采用土体固结方向的黏聚力并按各向同性评价缓坡的稳定性,可能严重高估安全系数,尤其是在黏聚力较高的情况下。与极限分析不同,建立的强度更新有限元模型能够分析成层边坡的稳定性。     

土坡稳定极限分析上限法中的多解性分析

土坡稳定极限分析上限法中采用了完全塑性区假定,即假定土坡内存在一个处处位于屈服面上的滑动体。将滑动体内部条块界面上的强度发挥系数视为优化变量,利用功能平衡方程得到了对于给定滑动面的多个安全系数,基于和声算法和潘家铮极值原理选择其中最大的安全系数作为该滑动面的抗滑稳定安全系数,利用该方法对于两个典型土坡算例进行了分析,将其结果与常规极限分析上限法结果进行了比较,给出了完全塑性区假定对计算结果的影响。     

Stability determination for layered soil slopes using the upper bound limit analysis

By using the upper bound limit analysis, a method has been proposed in this article to compute stability numbers for layered soil slopes. The rupture surface was assumed to be a combination of different logarithmic spiral arcs with a common focus. The analysis ensures the kinematic admissibility of the collapse mechanism with respect to the rigid rotation of the bounded soil mass about the focus of the logarithmic spiral. For the sake of illustration, the computations, however, were exclusively performed only for a two layered soil slope. The effect of the pore water pressure and horizontal earthquake body forces was also incorporated in the computations. The computational procedure was validated by making a comparison of the obtained stability numbers with those reported in literature.     

顺层岩质边坡稳定性极限分析上限法

基于极限分析上限法理论,运用体积力增量法,考虑单层滑动面极限分析模型的缺陷,建立考虑含结构面的多岩层错动的任意块体模型,利用岩块体在外荷载作用下达到极限破坏时,外荷载做的功与岩体消耗的功相等为基础建立等式方程,由此推导得到此类边坡稳定系数的计算公式。在分析典型边坡算例的基础上,对比离散元数值模拟结果,验证了该方法的可行性和适用性。对云南省普宣高速公路某顺层岩质边坡进行计算,其稳定系数接近离散元数值模拟结果,数值稍有偏大,边坡稳定。其研究结果为此类边坡稳定性评价提供了一种新的计算方法。     

基于Excel数据表和极限分析法的滑坡抗剪强度参数反演分析

针对采用传统的确定性边坡稳定性分析方法进行滑裂面抗剪强度参数反演分析所需滑坡断面不少于两个的不足,基于概率可靠度理论,提出了一种结合Excel数据表和极限分析上限法的滑坡强度参数反演新方法。依托湘西朱雀洞特大滑坡进行算例分析,研究结果表明：结合Excel数据表法和极限分析上限法反演得到的滑裂面抗剪强度参数与工程现场实测数据吻合良好,验证了新分析方法的正确性和有效性;运用该法反演分析只需一个滑坡断面,可有效解决单一滑裂面抗剪强度参数反演问题,减少了计算工作量,具有广泛的适用性和良好的工程应用前景。     

抗滑桩预加固边坡的能量分析方法

基于强度折减技术和岩土塑性极限分析理论,考虑抗滑桩预加固岩土边坡的加固效应,从功能平衡原理探讨了岩土边坡抗滑桩加固位置、桩长和多抗滑桩共同加固边坡的分析方法。假定机动容许的对数螺旋线速度机构,推导了抗滑桩预加固边坡稳定性安全系数的计算表达式。采用内点迭代方法对安全系数目标函数进行能量耗散最小化意义上的优化计算,算例对比验证了方法和程序的正确性。分析了单抗滑桩加固措施条件下的边坡局部稳定性特性,并由此推导了多抗滑桩加固边坡的稳定性分析算式。研究表明,抗滑桩设置位置改变对边坡安全系数、临界滑裂面位置和桩长均有显著影响,抗滑桩加固边坡的最优设置位置应位于边坡中下部区域     

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

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

基于MATLAB的边坡稳定性极限上限分析程序开发

基于系统能耗最小原理的极限分析上限定理,利用MATLAB强大的数据处理功能,对边坡稳定系数Ncr采用序列二次规划法（SQP）进行优化计算,并与枚举法相结合,从而解决SQP算法求解优化解容易陷入局部最优解的问题。同时开发了一款简单实用的边坡稳定性分析软件,实现了计算结果与临界滑裂破坏面图形的可视化输出。与经典算例对比分析,计算结果具有较好的一致性,表明了程序的可靠性。分析了内摩擦角?、边坡倾斜角?、? 3个参数对边坡稳定系数Ncr的影响。计算结果表明：边坡稳定性系数Ncr对?、? 较为敏感,而对? 的敏感次之。其次,当?、? 比较接近或在小范围内波动时,对稳定性系数Ncr的影响特别敏感;当?、?值较小时,滑动面通过坡趾下方,而当?值比较大时,滑动面通过坡趾。     

基于下限原理有限元的强度折减法

对于岩土工程中常用的强度折减系数,其规划问题是非线性的,不能直接利用线性规划进行求解。基于四边形单元的下限原理有限元法,根据强度折减系数与超载系数近似符合双曲函数的特点,通过调整强度参数使得超载系数逼近于1[1],可将边坡稳定性分析中常用的强度折减系数的非线性规划求解转化为线性规划求解问题。分析表明,采用拟合双曲线插值法求解强度折减系数的计算效率高于常规的二分法及割线法,且具有较好的收敛性;该方法能够充分利用当前高效的线性规划算法,便于工程应用。