岩质边坡关键块体的搜索方法及工程应用

工程边坡开挖所揭露的某些结构面组合块体在破坏原有静力平衡状态后,进行应力重新分配,进而块体会发生失稳、滑动等变形,甚至影响到整个边坡发生破坏,因此,结构面组合块体的稳定性是决定岩质边坡稳定的重要因素。本文通过对某工程开挖边坡地形和结构面三维实体模型的建立,实现了边坡结构面与开挖地形实体模型在AutoCAD平台上的可视化操作,剖切不同位置的二维图形预测分析了边坡开挖后可能存在的结构面组合块体,采用块体计算程序(Swegde)验证关键块体并进行稳定性评价。结果表明：三维实体模型及剖切技术与块体计算方法的结合,易于搜索关键块体。事实证明,该方法操作简便,计算结果准确,搜索方法具备可行性。     

Jennings边坡稳定性分析的概率方法(英文)

Jennings边坡 (简称J边坡 )是一种常见于格状与层状岩体中的典型岩质边坡模型 ,论文在岩质边坡中对J边坡运用蒙特卡罗方法进行概率分析 ,考虑陡、缓倾角两组节理随机组合 ,对计算范围内节理分布进行位置搜索 ,并考虑水压力及地震灾害条件下 ,计算相应坡坏概率及其可靠度指标。使用该方法对江西某铜矿高边坡工程进行了计算 ,取得令人满意的结果     

引入退火机制的智能单粒子算法在复杂边坡最危险滑动面搜索中的应用

复杂边坡的安全系数可能存在多个局部极小值点,如何确定边坡的最小安全系数是复杂边坡稳定性分析中的一个关键问题。本文结合简化Bishop法,采用一种新的启发式全局优化算法——智能单粒子算法(ISPO)来搜索复杂边坡的最危险滑动面。为帮助该算法快速跳出局部极值点,本文将模拟退火(SA)机制引入到智能单粒子算法中,结合了两种算法各自的优点,提出了引入退火机制的智能单粒子算法SA-ISPO。将本文提出的SA-ISPO算法用于搜索两个复杂边坡的最危险滑动面,并与其它方法相比较,验证了SA-ISPO算法的优越性,该算法搜索效率高,计算结果不受搜索范围的影响,是一种较好的全局优化算法。     

混沌优化算法在土质边坡稳定性分析中的应用

确定边坡最危险滑动面并计算与之相对应的安全系数是边坡支护的重要任务。本文结合简化Bishop法,用一种新的智能优化算法——混沌优化算法来搜索全局最优解。该方法利用混沌运动本身具有遍历性、随机性、“规律性”等内在特点,能在一定范围内按其自身“规律”不重复地遍历所有状态,易于跳出局部最优解,具有很强的全局搜索能力。通过坡高为12.3m的某电厂三层土质边坡的典型算例分析,并和遗传算法、枚举法计算结果对比可知,计算结果超于一致,其差值接近于0,因此混沌优化算法能在很高精度下搜索到全局最优解,能很好地解决边坡稳定性分析中的优化问题。     

基于代理模型的边坡稳定可靠度算法

针对基于极限平衡模式而建立的边坡稳定可靠度功能函数的隐式特征及其导致的可靠度求解等问题,首先利用Kriging插值预测方法,以极限平衡理论的Janbu算法为例,建立了边坡稳定可靠度功能函数的新型显式代理模型;然后将拉丁超立方试验抽样设计方法、Janbu算法、几何可靠度指标计算三者结合起来,研制了该代理模型的求解程序;构建了基于Kriging插值预测的边坡稳定可靠度近似分析新方法。选用某著名边坡算例,将本文方法分析结果与基于蒙特卡洛法的精确解进行了对比,两者的相对误差为5.21%,但前者计算工作量不到后者的万分之一,表明了代理模型方法的实用性和有效性。最后用代理模型方法分析了衡桂高速某段路堑边坡的稳定可靠性,表明本文方法计算结果能够满足工程精度要求,具有一定的应用价值。     

基于无网格—图论方法的岩体边坡稳定性分析研究

将无网格伽辽金法应用于岩体边坡稳定性分析,发展了基于无网格模型和有向加权图Bellman-Ford最短路径搜索算法相结合的无网格-图论边坡滑移面搜索方法,以搜寻节理岩体边坡失稳时的临界滑移面并得出其相应的安全系数.区别于传统的滑移面搜索算法,本文方法无需假定滑移面形状,更适用于具有复杂滑移线形状的节理岩体边坡的稳定性分析与计算,具有稳定和高效的算法特点.文中详细论述了无网格-图论最短路径算法的理论、方法和程序实现,并通过算例说明该方法在岩体边坡稳定性分析中的适用性.     

基于无网格一图论方法的岩体边坡稳定性分析研究

将无网格伽辽金法应用于岩体边坡稳定性分析,发展了基于无网格模型和有向加权图Bellman—Ford最短路径搜索算法相结合的无网格-图论边坡滑移面搜索方法,以搜寻节理岩体边坡失稳时的临界滑移面并得出其相应的安全系数。区别于传统的滑移面搜索算法,本文方法无需假定滑移面形状,更适用于具有复杂滑移线形状的节理岩体边坡的稳定性分析与计算,具有稳定和高效的算法特点。文中详细论述了无网格-图论最短路径算法的理论、方法和程序实现,并通过算例说明该方法在岩体边坡稳定性分析中的适用性。     

基于最大最小原理的三维边坡稳定性探讨

最大最小原理在边坡稳定性分析中的应用,不仅可以检验所建模型的准确性,还可以进一步获得改进计算模型的相关条件。本文利用三维简化Sarma法,分析了滑面剪切力方向采用不同的分布形式引起的稳定系数的变化,计算结果表明,不管采用何种分布形式,当它发生变化时,对应的稳定系数总是存在一个极大值,这就证明了边坡稳定性最大原理在三维边坡稳定分析中的适用性; 同时,本文还在一种严格的二维临界滑面搜索算法的基础上,建立了三维临界滑面搜索的算法,利用ZhangXing算例详细分析了目标函数中各变量对稳定系数的影响,结果表明,二维上寻找临界滑面是切实可行的,这也证明最小原理在二维边坡稳定分析中是满足的,三维上则由于三维极限平衡法未考虑三维边坡的边界条件,导致基于三维极限平衡法的三维临界滑面理论上不存在,三维边坡稳定分析结果与最小原理的要求有一定的差异。     

框架预应力锚杆边坡支护结构的稳定性分析方法及其 应用

以边坡极限平衡理论及框架预应力锚杆边坡支护结构的圆弧滑动破坏模式为基础,在考虑框架、锚杆对土体边坡稳定性影响的情况下,基于圆弧滑动条分法的思想利用积分法建立了内部稳定性安全系数计算模型和最危险滑移面搜索模型。确定滑移面圆心坐标为几何控制参数,推导了安全系数计算模型中各变量与滑移面圆心坐标之间的函数关系式,从而获得了滑移面几何控制参数与内部稳定性安全系数之间的函数关系。利用网格法对框架预应力锚杆支护结构最危险滑移面的圆心坐标进行动态搜索和求解。最后采用Matlab语言编制了框架预应力锚杆边坡支护结构的稳定性分析程序,并结合工程实例进行了验算,讨论了此种方法的适用条件,结果表明该方法简明适用,较传统的经验分析方法更加合理。     

昆明绕城高速公路杨林隧道出口的边坡稳定性研究

针对严格Janbu法不易收敛的缺点,采用Matlab曲线拟合工具箱里的三次样条函数csapi拟合条间力矩,简化了计算程序;采用改进的严格Janbu法等极限平衡方法和考虑地震动荷载时程分析的拟静力法对昆明绕城高速公路杨林隧道出口段某滑坡体的抗滑稳定性进行了系统地分析和研究,计算了该边坡在自然工况和地震工况条件下的抗滑稳定安全系数值。结果表明:在自然工况下边坡抗滑稳定安全系数在1.447～1.452之间,边坡是稳定的;在地震工况下,通过综合计算分析所得到的边坡抗滑稳定安全系数在0.55～1.017之间,边坡存在失稳风险。     

边坡稳定分析的虚功率法

基于有限元计算所得应力场的改进极限平衡法, 对渗流与有效应力耦合作用强烈的边坡或地基的稳定性分析具有优势. 本文提出了边坡稳定分析的虚功率法, 即基于极限分析的上限定理, 利用机动许可的组合刚体滑动机构和有限元应力场, 用滑动机构的速度间断面上的抗滑力功率与滑动力功率的比值计算安全系数. 通过分步优化方法, 获得边坡给定滑动机制的稳定安全系数. 对2个典型折线滑动机构的边坡案例进行了分析, 比较了采用静力平衡应力场和静力许可应力场对安全系数的影响. 指出基于土体线弹性本构模型所得的有效应力场计算的稳定安全系数, 也是边坡稳定安全性的一个不错的度量. 算例中本文计算所得的边坡稳定安全系数, 与文献中推荐的答案很接近, 其滑动机制与有限元强度折减法分析所得的滑动机制基本一致, 安全系数也接近, 表明本文提出的方法是合理的边坡稳定分析新方法, 为边坡和地基的稳定性分析提供了新的选择.      

A stability analysis of an oscillating body located in fluid flow using automatic differentiation

This paper presents a stability analysis of an oscillating body subjected to fluid forces located in a transient incompressible viscous flow. If the body is supported by elastic springs, oscillation will begin. If the characteristic period of the body and the excited oscillating period due to fluid forces match each other, resonance can occur. Stability analysis is therefore needed to determine the nonlinear behavior of the body. This paper presents an analysis of the changing stability of bodies by the numerical computation. To implement the computation, the motion of fluid around a body is expressed by the Navier–Stokes equation described in the arbitrary Lagrangian–Eulerian form. The fluid influence on the body is discretized by the finite element method based on a mixed interpolation by the bubble function in space. The motion of the body is assumed to be expressed by the equations of motion. To evaluate stability, stability function is defined by the total energy of the oscillating body. The stability is judged according to a stability index, obtained by the use of the automatic differentiation (AD) of the stability function. AD is a derivative computation method that gives high accuracy. By the use of AD, the second‐order derivative matrix, which is needed to compute the stability index, can be obtained exactly. For the numerical studies, analyses of one degree of freedom and two degrees of freedom (2DOF) for a circular cylinder and 2DOF for a rectangular cylinder are carried out. A combination of a cylinder and supporting elastic spring can produce stable, neutral and unstable states. It is shown that the stability of the cylinder can be determined by the stability index. This paper shows new possibilities for stability analysis of bodies located in a fluid flow. Copyright © 2008 John Wiley & Sons, Ltd.     

薄板分析的一种离散条法

建立薄板分析的一种离散条法，该方法通过引入满足边界条件的正交形函数，将目前用于梁，杆等一维结构的系统的离散元法扩展应用于薄板的力学分析。     

Comparative study of rock slope stability analysis methods for hydropower projects

This paper presents a comparative study of methods used for the analysis of the stability of rock slopes for hydropower projects. The comparison concerns the application of the limit equilibrium method, the block element method, and the finite element method on the stability analysis of two hydropower engineering projects in China. The study results show that the limit equilibrium method and the block element method give very close result, while the finite element method gives more conservative safety factor.     

索穹顶结构的找力分析方法--不平衡力迭代法

索穹顶的刚度和承载能力都是通过施加预应力来获取的,预应力设计是对索穹顶进行任何其它分析的基础,所以它的找力分析(Force finding)是一项很重要的工作。本文提出一种找力分析方法———不平衡力迭代法,它可以克服整体可行预应力法的一些不足:即使在索穹顶构件分组发生错误时仍可找到满足已知形状的预应力分布,并且该方法具有良好的稳定性和运算效率,算例表明不平衡力迭代法非常适合大型复杂索穹顶结构的找力分析。     

无单元法在薄板稳定问题中的应用

用无单元法研究了薄板的弹性稳定问题，从滑动最小二乘法和变分原理出发导出了薄板的无单元法几何刚度矩阵，编制了相应的计算程序，并给出了算例，结果表明，方法合理可行，且精度高于有限元。     

三维刚体极限平衡法程序的编制及其应用

西安市黑河水库右坝肩变形体是在旋转变形条件下形成的, 这样, 对它的稳定性评价必须按空间问题来考虑, 要用三维分析方法来计算它的稳定性。鉴于目前的三维计算方法比较复杂, 得出的结果实用性不大, 本文首次编制了三维刚体极限平衡法程序, 该程序既能计算稳定系数, 又能计算出相应稳定系数下的滑坡推力。应用该程序计算黑河水库右坝肩变形体的稳定性, 得出的结果与实际相吻合, 同时又为工程治理提供了依据。     

多自由度参激系统稳定性分析的数值解法

现有参激系统的动力稳定性问题研究主要集中在主不稳定区域上。为获得组合不稳定区域,基于Floquet方法,采用Bolotin方法在不同周期数下设解形式,结合特征值分析法得到确定多自由度参激系统动力不稳定区域的数值解法。对一个两自由度受周期轴向力的旋转轴系算例的稳定性分析,发现通过增加设解近似项数可获得高阶不稳定区域,且各阶不稳定区域边界随近似次数的增加逐渐趋于稳定,此外,增大阻尼可使各不稳定区域边界变得更加平滑。本文方法可用于一般多自由度周期参激阻尼系统,是一种简明易操作的直接数值解法。     

STABILITY ANALYSIS FOR TANK-ROOF BASED ON RELIABILITY

根据加肋拱顶稳定性校核公式,应用结构可靠度理论,建立了油罐拱顶可靠性分析模型,对拱顶稳定性进行分析计算,以确定油罐拱顶的使用安全性. 运用最小距离法和有限差分法对该模型的可靠度进行求解.结合工程实例,分析了顶板自支承作用、肋板的加强作用在外载荷作用下对拱顶可靠指标的影响规律,拟合出可靠指标与拱顶各参数间数值关系. 通过对现役油罐拱顶稳定性进行可靠性评定,验证了该模型在现役条件下分析拱顶稳定性的可行性.     

基于模态缺陷的受压球壳屈曲特性研究

针对受压球壳非线性屈曲过程，对含初始缺陷受压球壳的稳定性进行研究。根据EN1993-1-6（2007）规范，给出不同制造等级壳体的等效缺陷值计算方法；基于线性特征值分析的模态构型给出初始缺陷的分布。利用非线性有限元弧长法对球壳受压失稳过程进行数值模拟，得到屈曲前后球壳变形情况及全过程载荷-位移曲线。计算一致缺陷模态法和N阶缺陷模态法对应的球壳屈曲临界载荷，结果表明，受压球壳对缺陷较敏感，承载能力随缺陷值增大而降低；一致缺陷模态法计算便捷，在工程应用上具有合理性，N阶缺陷模态法考虑高阶模态缺陷构型，结果更加全面，可以为工程中缺陷结构稳定性设计提供参考。