吉林大学学报(地球科学版) ›› 2021, Vol. 51 ›› Issue (5): 1347-1355.doi: 10.13278/j.cnki.jjuese.20200292

• 岩土防灾与减灾 • 上一篇    下一篇

基于FLAC3D的三维边坡稳定性强度折减法计算效率改进算法及其应用

徐文刚1, 余旭荣1, 年廷凯2, 曹琦1, 曹爱武3, 裴振伟2   

  1. 1. 中国电建集团华东勘测设计研究院有限公司交通市政工程院, 杭州 311122;
    2. 大连理工大学土木工程学院, 辽宁 大连 116024;
    3. 浙江中科依泰斯卡岩石工程研发有限公司, 杭州 311122
  • 收稿日期:2020-12-07 出版日期:2021-09-26 发布日期:2021-09-29
  • 通讯作者: 年廷凯(1971-),男,教授,博士生导师,主要从事岩土工程方面的研究,E-mail:tknian@dlut.edu.cn E-mail:tknian@dlut.edu.cn
  • 作者简介:徐文刚(1985-),男,高级工程师,博士,主要从事公路边坡设计与滑坡治理方面的研究,E-mail:xu_wg@ecidi.com
  • 基金资助:
    国家自然科学基金项目（U1765107）

Optimization and Application of FLAC3D Strength-Reduction Computation in Three-Dimension Slope Stability Analysis

Xu Wengang1, Yu Xurong1, Nian Tingkai2, Cao Qi1, Cao Aiwu3, Pei Zhenwei2   

  1. 1. Transportation and Municipal Engineering Institute, Power China Huadong Engineering Corporation Limited, Hangzhou 311122, China;
    2. School of Civil Engineering, Dalian University of Technology, Dalian 116024, Liaoning, China;
    3. HydroChina ITASCA R&D Center, Hangzhou 311122, China
  • Received:2020-12-07 Online:2021-09-26 Published:2021-09-29
  • Supported by:
    Supported by the National Natural Science Foundation of China (U1765107)

PDF (PC)

119

摘要： 强度折减数值分析方法在边坡工程中应用广泛，其中FLAC3D已成为三维边坡稳定性分析中的主流计算工具。FLAC3D自带强度折减计算模块，因其计算过程复杂且耗时，一些学者提出了基于"二分法"的强度折减法，但依然不能令人满意。本文利用Python编程语言对FLAC3D进行二次开发，实现了"二分法"强度折减法的全面优化，包括初始应力场载入优化和计算时步上限优化，建议计算时步上限取初始应力场计算时步的1.50倍，此时既可以保证计算结果的准确性又可以节省计算时间。利用优化后方法对某公路复杂三维填方路基边坡进行了三维强度折减法计算，结果表明，在保证计算准确性的同时，本文方法计算效率较其他改进的强度折减法至少提升50%，较程序自带模块提升了300%以上。

关键词: 三维边坡, 强度折减法, FLAC3D, 公路边坡, Python

Abstract: The strength-reduction method (SRM) is widely applied in slope engineering projects, and FLAC3D has become the main tool for 3D slope stability analysis. Although FLAC3D has SRM intrinsic modular to complete this task, the computing process is complicated and time-consuming, hence the bisectional method is then involved to simplify the process by some researchers, which, however,is defactive in some aspects. An optimized bisectional method is supplied in this paper based on Python, including in-situ stress field reloading and selection of the upper limit of time-steps, which obviously improves computing efficiency. By using this optimized method, a fill slope in road engineering project was analyzed, showing that the strength reduced factor is consistent with those of the other methods and with much higher computing efficiency (ca. 50% and 300% higher than the other improved methods and the intrinsic modular, respectively).

Key words: three-dimension slope, strength-reduction method, FLAC3D, road slope, Python

中图分类号: 

  • U417.2
[1] Arunav C, Diganta G. State of the Art:Three Dimensional (3D) Slope-Stability Analysis[J]. International Journal of Geotechnical Engineering, 2016, 10(5):493-498.
[2] Lin H D, Wang W C, Li A J. Investigation of Dilatancy Angle Effects on Slope Stability Using the 3D Finite Element Method Strength Reduction Technique[J]. Computers and Geotechnics, 2020, 118:103295.
[3] Yu G, Xie M W, Liang J, et al. A GIS-Based 3D Slope Stability Analysis Method Based on the Assumed Normal Stress on the Slip Surface[J]. Scientific Reports, 2020, 10:4442.
[4] Xu J S, Yang X L. Three-Dimensional Stability Analysis of Slope in Unsaturated Soils Considering Strength Nonlinearity Under Water Drawdown[J]. Engineering Geology, 2018, 237:102-115.
[5] 张国祥, 刘宝琛. 边坡滑动面三维空间分析[J]. 中国公路学报, 2003, 16(4):25-29. Zhang Guoxiang, Liu Baochen. Analysis of Slope Slip Surface by Three-Dimension Finite Element Method[J]. China Journal of Highway and Transport, 2003, 16(4):25-29.
[6] 王林峰, 唐红梅, 唐芬, 等. 复杂缓倾角岩体边坡的三维稳定性分析[J]. 中国公路学报, 2018, 31(2):57-66. Wang Linfeng, Tang Hongmei, Tang Fen, et al. Three-Dimension Stability Analysis of Complex Gentle Incline Rock Slope[J]. China Journal of Highway and Transport, 2018, 31(2):57-66.
[7] 张标, 王晅, 张家生, 等. 地震效应下二级边坡三维动态稳定性上限分析[J]. 中国公路学报, 2018, 31(2):86-96. Zhang Biao, Wang Xuan, Zhang Jiasheng, et al. Three-Dimension Seismic Stability Limit Analysis of Two-Stage Slope by Kinematical Approach[J]. China Journal of Highway and Transport, 2018, 31(2):86-96.
[8] 侯超群, 邓欣, 孙志彬, 等. 非线性准则下三维加筋边坡稳定性的上限分析[J]. 中国公路学报, 2018, 31(2):124-132. Hou Chaoqun, Deng Xin, Sun Zhibin, et al. Upper Bound Analysis of Stability of Three-Dimension Reinforced Slopes Based on Nonlinear Failure Criterion[J]. China Journal of Highway and Transport, 2018, 31(2):124-132.
[9] Nian T K, Huang R Q, Wan S S, et al. Three-Dimensional Strength Reduction Finite Element Analysis of Slopes:Geometric Effects[J]. Canadian Geotechnical Journal, 2012, 49(5):574-588.
[10] Arellano D, Stark T D. Importance of Three-Dimensional Slope Stability Analysis in Practice[M]. New York:American Society of Civil Engineers, 2000:18-32.
[11] 韩同春, 林博文, 何露, 等. 基于GIS与数值模拟软件耦合的三维边坡建模方法及其稳定性研究[J]. 岩土力学, 2019, 40(7):2855-2865. Han Tongchun, Lin Bowen, He Lu, et al. Three-Dimensional Slope Modelling Method and Its Stability Based on Coupled GIS and Numerical Simulation Software[J]. Rock and Soil Mechanics, 2019, 40(7):2855-2865.
[12] 邓东平, 李亮. 基于滑动面应力假设下的三维边坡稳定性极限平衡法研究[J]. 岩土力学, 2017, 38(1):189-196. Deng Dongping, Li Liang. Three-Dimensional Limit Equilibrium Method for Slope Stability Based on Assumption of Stress on Slip Surface[J]. Rock and Soil Mechanics, 2017, 38(1):189-196.
[13] Wang L, Sun D A, Li L. Three-Dimensional Stability of Compound Slope Using Limit Analysis Method[J]. Canadian Geotechnical Journal, 2019, 56(1):116-125.
[14] Griffiths D V, Lane P A. Slope Stability Analysis by Finite Elements[J]. Geotechnique, 1999, 49(3):387-403.
[15] Usluogullari O F, Temugan A, Duman E S. Comparison of Slope Stabilization Methods by Three-Dimensional Finite Element Analysis[J]. Natural Hazards, 2016, 81:1027-1050.
[16] 郑颖人, 赵尚毅. 有限元强度折减法在土坡与岩坡中的应用[J]. 岩石力学与工程学报, 2004, 23(19):3381-3388. Zheng Yingren, Zhao Shangyi. Application of Strength Reduction FEM to Soil and Rock Slope[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(19):3381-3388.
[17] 杨有成, 李群, 陈新泽, 等. 对强度折减法若干问题的讨论[J]. 岩土力学, 2008, 29(4):1103-1106. Yang Youcheng, Li Qun, Chen Xinze, et al. Discussion on Strength Reduction Using FLAC[J]. Rock and Soil Mechnics, 2008, 29(4):1103-1106.
[18] 年廷凯, 张克利, 刘红帅, 等. 基于强度折减法的三维边坡稳定性与破坏机制[J]. 吉林大学学报(地球科学版), 2013, 43(1):178-185. Nian Tingkai, Zhang Keli, Liu Hongshuai, et al. Stability and Failure Mechanism of Three-Dimensional Slope Using Strength Reduction Method[J]. Journal of Jilin University(Earth Science Edition), 2013, 43(1):178-185.
[19] 郑颖人, 陈祖煜, 王恭先, 等. 边坡与滑坡工程治理[M]. 第二版, 北京:人民交通出版社, 2010:194-196. Zheng Yingren, Chen Zuyu, Wang Gongxian, et al. Engineering Treatment of Slope and Landslide[M]. 2nd ed. Beijing:China Communications Press, 2010:194-196.
[20] 栾茂田, 武亚军, 年廷凯. 强度折减有限元法边坡失稳的塑性区判据及其应用[J]. 防灾减灾工程学报, 2003, 23(3):1-8. Luan Maotian, Wu Yajun, Nian Tingkai. A Criterion for Evaluating Slope Stability Based on Development of Plastic Zone by Shear Strength Reduction FEM[J]. Journal of Disaster Prevention and Mitigation Engineering, 2003, 23(3):1-8.
[21] Itasca Consulting Group. Fast Lagrange Analysis of Continua in 3 Dimensions, User's Manual[M]. USA, Minneapolis:Itasca Consulting Group, 2005.
[22] Wu J, Li Y, Cheng Q, et al. A Simplified Method for the Determination of Vertically Loaded Pile-Soil Interface Parameters in Layered Soil Based on FLAC3D[J]. Frontiers of Structural and Civil Engineering, 2016, 10:103-111.
[23] Hao F S, Wang L. Application Study of FLAC in Analysis of Slope Stability[J]. Physical and Numerical Simulation of Geotechnical Engineering, 2016, 23(1):3-8.
[24] 宋民崇, 余云燕. 三维数值模拟在某高速公路边坡稳定性分析中的应用[J]. 黑龙江工程学院学报, 2016, 30(1):4-23. Song Minchong, Yu Yunyan. Application of Three-Dimensional Numerical Simulation to the Stability Analysis of a Highway Slope[J]. Journal of Heilongjiang Institute of Technology, 2016, 30(1):4-23.
[25] 陈育民, 徐鼎平.FLAC/FLAC3D基础与工程实例[M]. 北京:中国水利水电出版社, 2008:267-269. Chen Yumin, Xu Dingping. FLAC/FLAC3D Foundation and Engineering Application[M]. Beijing:China Water & Power Press, 2008:267-269.
[26] 苏飞, 李博. 基于FLAC3D的强度折减法程序改进研究[J]. 工程与材料科学, 2017, 35:72-78. Su Fei, Li Bo. Improvement of Strength Reduction Program Based on FLAC3D[J]. Engineering and Material Sciences, 2017, 35:72-78.
[27] Cai F, Ugai K. Numerical Analysis of the Stability of a Slope Reinforced with Piles[J]. Soils and Foundations, 2000, 40(1):73-84.
[28] Wei W B, Cheng Y M. Strength Reduction Analysis for Slope Reinforced with One Row of Piles[J]. Computers and Geotechnics, 2009, 36:1176-1185.
[29] Cheng Y M, Yip C J. Three-Dimensional Asymmetrical Slope Stability Analysis Extension of Bishop's and Janbu's Techniques[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2007, 133(12):1544-1555.
[30] Chugh K A. On the Boundary Conditions in Slope Stability Analysis[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2003, 27(11):905-926.
[31] Arai K, Tagyo K. Determination of Noncircular Slip Surface Giving the Minimum Factor of Safety in Slope Stability Analysis[J]. Soils and Foundations, 1985, 25(1):43-51.
[32] Greco V R. Efficient Monte Carlo Technique for Locating Critical Slip Surface[J]. Journal of Geotechnical Engineering, 1996, 122(7):517-525.
[33] Kim J, Salgado R, Lee J. Stability Analysis of Complex Soil Slopes Using Limit Analysis[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2002, 128(7):546-557.
[34] Tandjiria V, Low B K, The C I. Effect of Reinforcement Force Distribution on Stability of Embankments[J]. Geotextiles and Geomembranes, 2002, 20:423-443.
[1] 魏家斌, 王卫东, 吴江斌. 免共振沉桩过程对地表振动影响的FLAC3D数值模拟[J]. 吉林大学学报(地球科学版), 2021, 51(5): 1514-1522.
[2] 张海清, 贾会会, 聂庆科. 土体抗拉强度对均质边坡稳定性的影响[J]. 吉林大学学报(地球科学版), 2021, 51(5): 1324-1337.
[3] 熊良宵，李天斌，杨昌斌，虞利军. 层状岩体单轴和双轴压缩蠕变特性的数值试验[J]. 吉林大学学报(地球科学版), 2013, 43(5): 1539-1548.
[4] 王科，王常明，王彬, 姚康，王天佐. 基于MorgensternPrice法和强度折减法的边坡稳定性对比分析[J]. 吉林大学学报(地球科学版), 2013, 43(3): 902-907.
[5] 于子望，张延军，张庆，许天福. TOUGHREACT搭接FLAC3D算法[J]. 吉林大学学报(地球科学版), 2013, 43(1): 199-206.
[6] 刘红帅, 年廷凯, 万少石. 三维边坡稳定性分析中的边界约束效应[J]. J4, 2010, 40(3): 638-644.
[7] 张晨，陈剑平,肖云华. 基于神经网络对有限元强度折减法分析[J]. J4, 2009, 39(1): 114-0118.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
[1] 李春柏，张新涛，刘 立，任延广,孟 鹏. 布达特群热流体活动及其对火山碎屑岩的改造作用--以海拉尔盆地贝尔凹陷为例[J]. J4, 2006, 36(02): 221 -0226 .
[2] 邹新宁，孙 卫，张盟勃,万玉君. 地震属性分析在岩性气藏描述中的应用[J]. J4, 2006, 36(02): 289 -0294 .
[3] 郭洪金,李勇,钟建华,王海侨. 山东东辛油田古近系沙河街组一段碳酸盐岩储集特征[J]. J4, 2006, 36(03): 351 -357 .
[4] 杜业波，季汉成，朱筱敏. 川西前陆盆地上三叠统须家河组成岩相研究[J]. J4, 2006, 36(03): 358 -364 .
[5] 刘家军，李志明，刘建明，王建平，冯彩霞，卢文全. 自然界中的辉锑矿-硒锑矿矿物系列[J]. J4, 2005, 35(05): 545 -553 .
[6] 苏继军，殷 琨，郭同彤. 金刚石绳索取心钻杆接头螺纹的优化研究[J]. J4, 2005, 35(05): 677 -680 .
[7] 唐健生，夏日元，邹胜章,梁 彬. 新疆南天山岩溶系统介质结构特征及其水文地质效应[J]. J4, 2005, 35(04): 481 -0486 .
[8] 熊 彬. 大回线瞬变电磁法全区视电阻率的逆样条插值计算[J]. J4, 2005, 35(04): 515 -0519 .
[9] 杜春国,邹华耀,邵振军，张俊. 砂岩透镜体油气藏成因机理与模式[J]. J4, 2006, 36(03): 370 -376 .
[10] 许盛伟,王明常,白亚辉,张学明. 基于J2EE的分布式海量影像分发服务研究和实现[J]. J4, 2006, 36(03): 491 -496 .
版权所有 © 《吉林大学学报(地球科学版)》编辑部
地址：长春市西民主大街938号(130026)　电话：+86-431-88502374；13756165364　E-mail: xuebao1956@jlu.edu.cn
本系统由北京玛格泰克科技发展有限公司设计开发