长江科学院院报 ›› 2016, Vol. 33 ›› Issue (12): 61-67.DOI: 10.11988/ckyyb.20150680

• 岩土工程 • 上一篇    下一篇

基于二阶锥规划的边坡稳定上限有限元分析

杨昕光1,2,周密1,张伟1,潘家军1   

  1. 1.长江科学院 水利部岩土力学与工程重点实验室,武汉 430010;
    2.水利部土石坝破坏机理与防控技术重点实验室, 南京 210029
  • 收稿日期:2015-08-12 出版日期:2016-12-20 发布日期:2016-12-13
  • 作者简介:杨昕光(1983-),男,内蒙古赤峰人,工程师,博士,研究方向为土石坝地震工程与土工数值分析,(电话)18694049883(电子信箱)yyfreshman@163.com。
  • 基金资助:
    国家自然科学基金项目(51509019);长江科学院创新团队项目(CKSF2015051/YT);水利部土石坝破坏机理与防控技术重点实验室开放研究基金项目(YK914017)

Upper Bound Finite Element Limit Analysis of Slope Stability Using Second-order Cone Programming

YANG Xin-guang1,2 ,ZHOU Mi1, ZHANG Wei1, PAN Jia-jun1   

  1. 1.Key Laboratory of Geotechnical Mechanics and Engineering of the Ministry of Water Resources, Yangtze River Scientific Research Institute, Wuhan 430010, China;
    2.Key Laboratory of Failure Mechanism and Safety Control Techniques of Earth-rock Dam of the Ministry of Water Resources, Nanjing 210029, China
  • Received:2015-08-12 Published:2016-12-20 Online:2016-12-13

摘要: 边坡稳定上限极限分析通常采用三节点三角形单元,并对屈服准则线性化,因而其计算精度较低。针对此问题,基于六节点三角形单元和二阶锥规划,假设材料严格满足Mohr-Coulomb屈服准则并考虑孔隙水压力与地震荷载的作用,发展了一种应用于边坡稳定性分析的上限有限元法,扩大了上限极限分析的应用范围。根据上限定理,在满足屈服条件、流动法则、功能平衡以及相应的边界条件的基础上,将边坡稳定的上限分析形成二阶锥规划数学模型,并用先进的内点法进行求解。通过2个算例的计算分析,并与已有计算方法进行对比,验证了本文方法的正确性与可行性,并表明该方法的计算结果并不明显依赖于有限元网格的密度,且材料内摩擦角较大时,仍具有较高的计算精度。

关键词: 边坡稳定性, 极限分析, 上限定理, 二阶锥规划, 有限元

Abstract: Three-node triangular elements together with the linearized yield criterion are commonly used in the finite element upper bound limit analysis. Therefore, the method is of low calculation precision. Aiming at this problem, a method of upper bound limit analysis using six-node triangular elements and second-order cone programming is developed to investigate the slope stability subjected to pore water pressure and earthquake loads. The proposed method formulates the slope stability problem as a second-order cone programming with constraints based on the yield criterion, flow rule, boundary conditions, and the energy-work balance equation. The optimization problem is solved by a state-of-the-art interior-point method, and the strict upper bound solutions can be obtained. Finally, the results of two numerical examples are compared with published solutions, which demonstrate the validity of the proposed method. The results also indicate that the mesh-dependence phenomenon is overcome and the calculation precision is improved even for large internal friction angel of materials.

Key words: slope stability, limit analysis, upper bound theorem, second-order cone programming, finite element

中图分类号: