岩土力学 ›› 2019, Vol. 40 ›› Issue (10): 4084-4091.doi: 10.16285/j.rsm.2019.0226

• 数值分析 • 上一篇    下一篇

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

刘锋涛1,张绍发1,戴北冰2,张澄博1,林凯荣2   

  1. 1. 中山大学 地球科学与工程学院 广东省地球动力作用与地质灾害重点实验室,广东 广州 510275; 2. 中山大学 土木工程学院,广东 广州 510275
  • 收稿日期:2019-01-28 出版日期:2019-10-11 发布日期:2019-10-20
  • 通讯作者: 戴北冰,男,1981年生,博士,副教授,主要从事岩土工程数值计算方法方面的研究。E-mail: daibb@mail.sysu.edu.cn E-mail: liuft@mail.sysu.edu.cn
  • 作者简介:刘锋涛,男,1979年生,博士,讲师,主要从事边坡稳定分析及数值方法方面的研究。
  • 基金资助:
    国家自然科学基金(No. 51209237,No.51779279);广州市科技计划项目(No. 201707010082);中央高校基本科研业务费(No. 18lgp16)。

Upper bound limit analysis of soil slopes based on rigid finite element method and second-order cone programming

LIU Feng-tao1, ZHANG Shao-fa1, DAI Bei-bing2, ZHANG Cheng-bo1, LIN Kai-rong2   

  1. 1. Guangdong Provincial Key Lab of Geodynamics and Geohazards, School of Geoscience and Engineering, Sun Yat-Sen University, Guangzhou, Guangdong 510275, China; 2. School of Civil Engineering, Sun Yat-Sen University, Guangzhou, Guangdong 510275, China
  • Received:2019-01-28 Online:2019-10-11 Published:2019-10-20
  • Supported by:
    This work was supported by the National Natural Science Foundation of China (51209237, 51779279), the Science and Technology Program of Guangzhou, China (201707010082) and the Chinese Universities Scientific Fund (18lgp16).

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

关键词: 边坡稳定, 刚体有限元上限法, 二阶锥规划, 临界滑动面

Abstract: The upper bound approach based on rigid finite element method (RFEM) is one of the limit analysis methods, and it is an important method for stability analysis of soil slopes. However, there are still some key problems to be overcome in this method. Due to the assumption of rigidity, the plastic deformation energy of the system is stored only at the interfaces of all elements. Therefore, the accuracy of the approach is highly dependent on the alignment of the interfaces, which indicates that the accuracy is poor if an unstructured mesh is employed. To overcome this limitation, we propose a novel RFEM-based upper bound approach using perturbation method, according to the sequential limit analysis. Firstly, considering the rotation of elements, a novel second-order cone programming (SOCP) is put forward to construct the kinematically admissible velocity field based on RFEM. Secondly, a model for searching the critical slip surface is built by using sequential SOCP of RFEM-based upper bound approach and then solved by the nonlinear simplex and particle swarm optimization algorithm. From the analysis on stability problems of two benchmark soil slopes, the proposed method is verified. The main influence factors on computational efficiency and accuracy are also investigated in this paper. We have found that the type of mesh is a significant factor affecting the accuracy of the proposed method. It is necessary to consider the rotation of element in RFEM-based upper bound approach, which can not only improve the accuracy of the approach but also overcome the dependence of RFEM-based upper bound limit analysis on the alignment of interfaces.

Key words: slope stability, RFEM-based upper bound limit analysis, second-order cone programming, critical slip surface

中图分类号: 

  • TU 433
[1] 史振宁, 戚双星, 付宏渊, 曾铃, 何忠明, 方睿敏, . 降雨入渗条件下土质边坡含水率分 布与浅层稳定性研究[J]. 岩土力学, 2020, 41(3): 980-988.
[2] 苏永华, 李诚诚. 强降雨下基于Green-Ampt模型的边坡稳定性分析[J]. 岩土力学, 2020, 41(2): 389-398.
[3] 孙锐, 杨峰, 阳军生, 赵乙丁, 郑响凑, 罗静静, 姚捷, . 基于二阶锥规划与高阶单元的 自适应上限有限元研究[J]. 岩土力学, 2020, 41(2): 687-694.
[4] 刘红岩. 宏细观缺陷对岩体力学特性及边坡稳定影响研究[J]. 岩土力学, 2019, 40(S1): 431-439.
[5] 韩同春, 林博文, 何露, 苏钰钦, . 基于GIS与数值模拟软件耦合的 三维边坡建模方法及其稳定性研究[J]. 岩土力学, 2019, 40(7): 2855-2865.
[6] 余 国, 谢谟文, 郑正勤, 覃事河, 杜 岩, . 基于GIS的边坡稳定性计算方法研究[J]. 岩土力学, 2019, 40(4): 1397-1404.
[7] 王冬勇, 陈曦, 于玉贞, 吕彦楠, . 基于二阶锥规划有限元增量加载法的条形浅基础极限承载力分析[J]. 岩土力学, 2019, 40(12): 4890-4896.
[8] 夏侯云山, 张抒, 唐辉明, 刘晓, 吴琼, . 考虑参数空间变异结构的结构化交叉约束 随机场模拟方法研究[J]. 岩土力学, 2019, 40(12): 4935-4945.
[9] 唐洪祥, 韦文成. 耦合强度各向异性与应变软化的边坡稳定 有限元分析[J]. 岩土力学, 2019, 40(10): 4092-4100.
[10] 刘素锦, 郭明伟, 李春光, . 三维边坡整体主滑方向的确定[J]. 岩土力学, 2018, 39(S2): 37-44.
[11] 代仲海,胡再强,尹小涛,吴振君,. 工程荷载作用下缓倾角反倾似层状岩质边坡变形稳定性分析[J]. , 2018, 39(S1): 412-418.
[12] 秦雨樵,汤 华,冯振洋,尹小涛,王东英, . 基于聚类分析的边坡稳定性研究[J]. , 2018, 39(8): 2977-2983.
[13] 郭重阳,李典庆,曹子君,高国辉,唐小松. 考虑空间变异性条件下的边坡稳定可靠度高效敏感性分析[J]. , 2018, 39(6): 2203-2210.
[14] 吴震宇,陈建康. 土坡体系可靠度分析方法及在高土石坝工程中的应用[J]. , 2018, 39(2): 699-704.
[15] 郑 刚,于晓旋,杜 娟,尹 鑫,周海祚,杨新煜, . 临近边坡的条形基础地基极限承载力数值分析[J]. , 2018, 39(10): 3812-3820.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!