›› 2014, Vol. 35 ›› Issue (S2): 626-637.

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

岩石弹塑性应力-渗流-损伤耦合模型研究(Ⅰ):模型建立及其数值求解程序

王军祥1,姜谙男1,宋战平2   

  1. 1. 大连海事大学 道路与桥梁工程研究所,辽宁 大连 116026;2. 西安建筑科技大学 土木工程学院,西安 710055
  • 收稿日期:2014-02-14 出版日期:2014-10-31 发布日期:2014-11-12
  • 通讯作者: 姜谙男,男,1971年生,博士,教授,从事岩土工程稳定性分析和多场耦合方面的研究。E-mail: jiangannan@163.com E-mail:w.j.xgood@163.com
  • 作者简介:王军祥,男,1985年生,博士研究生,从事岩石多场耦合机理及数值模拟方面的研究
  • 基金资助:

    国家自然科学基金(No. 51079010);中央高校基本科研业务费专项资金(No. 2013YB03,No. 3132014326,);大连市交通科技项目 (No. 2011-10);吉林省交通厅交通运输科技项目(No. 2012-1-6)。

Study of the coupling model of rock elastoplastic stress-seepage-damage (Ⅰ) : Modelling and its numerical solution procedure

WANG Jun-xiang1,JIANG An-nan1,SONG Zhan-ping2   

  1. 1. Institute of Road and Bridge Engineering, Dalian Maritime University, Dalian, Liaoning 116026, China; 2. College of Civil Engineering, Xi’an University of Architecture and Technology, Xi’an 710055, China
  • Received:2014-02-14 Online:2014-10-31 Published:2014-11-12

摘要: 在地下水渗流场、应力场、损伤场的耦合作用下更易造成隧道围岩坍塌或涌水等灾害。首先,将围岩材料视作各向同性连续介质,基于Drucker-Prager准则建立岩石弹塑性损伤本构模型,采用完全隐式返回映射算法实现弹塑性损伤本构方程的数值求解。其次,以上述研究为基础根据岩石处于弹塑性状态时渗透系数动态演化公式,建立岩石弹塑性应力-渗流-损伤耦合模型,并给出三场耦合情况下的数值求解迭代方法。针对耦合模型中涉及参数较多且不易测定的问题,基于差异进化算法原理建立智能反分析方法,对耦合模型中的损伤参数进行反演。最后,利用C++语言编制相应的岩石弹塑性应力-渗流-损伤耦合程序和参数反演程序,利用所编程序进行以下计算:(1)对智能反分析程序的性能、正确性进行分析,对比不同差异策略、交叉因子、变异因子的反演精度和收敛速度。(2)分别采用弹性模型和弹塑性损伤模型进行隧道围岩位移场、应力场的计算。(3)不考虑力学作用的情况下进行孔隙水压力、渗流量的计算。(4)采用所建耦合模型计算得到隧道围岩应力场、渗流场以及损伤场的相互影响规律。研究结果表明,基于差异进化算法的智能反分析程序能够较好地解决耦合模型中损伤参数不易确定的难题,为实际工程中获得不易测定的计算参数提供了有效的方法,同时所建立的耦合模型通过应力、渗流和损伤的相互作用更能够真实地反映出岩石材料的宏观破坏现象,所编计算程序能够模拟地下水渗流场、应力场、损伤场之间的耦合特性,为受地下水影响严重的工程建设提供了方法,研究结论为后期对实际隧道工程进行耦合计算奠定基础。

关键词: 弹塑性损伤本构模型, 渗透系数, 弹塑性应力-渗流-损伤耦合模型, 损伤参数反演

Abstract: It is more likely to cause the collapse of tunnel surrounding rock and water disasters under the action of the coupling effect of the groundwater seepage field, the stress field and the damage field. First of all, surrounding rock is regarded as the isotropic continuum material, and the elastoplastic damage constitutive model of rock is established based on the Drucker Prager criterion. The fully implicit return mapping algorithm is adopted to achieve the numerical solution of elastoplastic damage constitutive equations. Secondly, according to the dynamic evolution formula of permeability coefficient in elastoplastic state based on the above research, the elasoplastic stress-seepage-damage model of rock is established; and the coupling solving iterative method of many fields is presented. Back analysis program is compiled on the basis of the principle of differential evolution algorithm in order to solve the problem that elastoplastic damage model involves many parameters and not easy to determine, and the damage parameters are inversion; Finally, the corresponding coupling program of the elastoplastic stress-seepage-damage is compiled and parameters inversion program are compiled using C++ language respectively. The following calculations are conducted by using the programs: (1) The performance and the correctness of the intelligent back analysis program are analyzed, comparing the inversion accuracy and the convergence rate of the differences between the different strategies, the cross factors and the variation factors. (2) The calculation results of the displacement field and the stress field by the elastic model and the elastoplastic model are compared. (3) The pore water pressure and the seepage flow are calculated without considering the effect of mechanics. (4) The interaction rules of the stress field, the seepage field and the damage field of tunnel surrounding rock are calculated using the coupling model. The results show that the problem of damage parameters which is difficult to determine can be better to solve in damage model by the back analysis program based on differential evolution algorithm. It provides effective method for the calculation parameters are acquired not easily determined in practical engineering. At the same time, the macro failure of rock material can be more realistically reflect through the interaction of the stress, the seepage and the damage of the built coupling model. The calculation program which provides a method for the engineering construction affected by groundwater severely can simulate the coupling characteristics of the groundwater seepage field, the stress field and the damage field. The study of this paper lays a solid foundation for the coupling calculation of the later period for actual tunnel engineering.

Key words: elastoplastic damage constitutive model, permeability coefficient, elastoplastic stress-seepage-damage coupling model, damage parameters inversion

中图分类号: 

  • TU 452
[1] 李红坡, 陈征, 冯健雪, 蒙宇涵, 梅国雄, . 双层地基水平排水砂垫层位置优化研究[J]. 岩土力学, 2020, 41(2): 437-444.
[2] 彭家奕, 张家发, 沈振中, 叶加兵, . 颗粒形状对粗粒土孔隙特征和渗透性的影响[J]. 岩土力学, 2020, 41(2): 592-600.
[3] 王刚, 韦林邑, 魏星, 张建民, . 压实黏土三轴压缩变形过程中的渗透性变化规律[J]. 岩土力学, 2020, 41(1): 32-38.
[4] 刘丽, 吴羊, 陈立宏, 刘建坤, . 基于数值模拟的湿润锋前进法测量精度分析[J]. 岩土力学, 2019, 40(S1): 341-349.
[5] 徐浩青, 周爱兆, 姜朋明, 刘顺青, 宋苗苗, 陈亮, . 不同砂−膨润土垂直防渗墙填筑土料的掺量研究[J]. 岩土力学, 2019, 40(S1): 424-430.
[6] 张玉国, 万东阳, 郑言林, 韩帅, 杨晗玥, 段萌萌. 考虑径向渗透系数变化的真空预压 竖井地基固结解析解[J]. 岩土力学, 2019, 40(9): 3533-3541.
[7] 胡明鉴, 崔 翔, 王新志, 刘海峰, 杜 韦, . 细颗粒对钙质砂渗透性的影响试验研究[J]. 岩土力学, 2019, 40(8): 2925-2930.
[8] 李 贤, 汪时机, 何丙辉, 沈泰宇, . 土体适用MICP技术的渗透特性条件研究[J]. 岩土力学, 2019, 40(8): 2956-2964.
[9] 范日东, 刘松玉, 杜延军, . 基于改进滤失试验的重金属污染 膨润土渗透特性试验研究[J]. 岩土力学, 2019, 40(8): 2989-2996.
[10] 余良贵, 周建, 温晓贵, 徐杰, 罗凌晖, . 利用HCA研究黏土渗透系数的标准探索[J]. 岩土力学, 2019, 40(6): 2293-2302.
[11] 陶高梁, 吴小康, 甘世朝, 肖衡林, 马 强, 罗晨晨, . 不同初始孔隙比下非饱和黏土渗透性 试验研究及模型预测[J]. 岩土力学, 2019, 40(5): 1761-1770.
[12] 张 昭, 程靖轩, 刘奉银, 齐吉琳, 柴军瑞, 李会勇, . 基于土颗粒级配预测非饱和 渗透系数函数的物理方法[J]. 岩土力学, 2019, 40(2): 549-560.
[13] 刘一飞, 郑东生, 杨 兵, 祝 兵, 孙明祥. 粒径及级配特性对土体渗透系数影响的细观模拟[J]. 岩土力学, 2019, 40(1): 403-412.
[14] 朱长歧, 崔 翔, 胡明鉴, 王新志, 唐健健, . 钙质土电导率和渗透性的相关研究[J]. 岩土力学, 2018, 39(S2): 142-148.
[15] 马瑞男, 郭红仙, 程晓辉, 刘景儒, . 微生物拌和加固钙质砂渗透特性试验研究[J]. 岩土力学, 2018, 39(S2): 217-223.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 孙树林,李 方,谌 军. 掺石灰黏土电阻率试验研究[J]. , 2010, 31(1): 51 -55 .
[2] 李英勇,张顶立,张宏博,宋修广. 边坡加固中预应力锚索失效机制与失效效应研究[J]. , 2010, 31(1): 144 -150 .
[3] 王丽艳,姜朋明,刘汉龙. 砂性地基中防波堤地震残余变形机制分析与液化度预测法[J]. , 2010, 31(11): 3556 -3562 .
[4] 李秀珍,王成华,邓宏艳. DDA法和Fisher判别法在潜在滑坡判识中的应用比较[J]. , 2011, 32(1): 186 -192 .
[5] 吉武军. 黄土隧道工程问题调查分析[J]. , 2009, 30(S2): 387 -390 .
[6] 陈力华 ,林 志 ,李星平. 公路隧道中系统锚杆的功效研究[J]. , 2011, 32(6): 1843 -1848 .
[7] 陈立文,孙德安. 不同应力路径下水土耦合超固结黏土分叉分析[J]. , 2011, 32(10): 2922 -2928 .
[8] 郑 刚 张立明 刁 钰. 开挖条件下坑底工程桩工作性状及沉降计算分析[J]. , 2011, 32(10): 3089 -3096 .
[9] 赵明华,雷 勇,张 锐. 岩溶区桩基冲切破坏模式及安全厚度研究[J]. , 2012, 33(2): 524 -530 .
[10] 马 刚 ,常晓林 ,周 伟 ,周创兵 . 基于Cosserat理论的重力坝深层抗滑稳定分析[J]. , 2012, 33(5): 1505 -1512 .