›› 2011, Vol. 32 ›› Issue (11): 3515-3520.

• 数值分析 • 上一篇    

基于Zienkiewicz-Pande屈服准则的弹塑性本构模型在FLAC3D中的二次开发及应用

左双英1, 2，肖 明1，陈俊涛1   

  1. 1. 武汉大学 水资源与水电工程科学国家重点实验室，武汉 430072；2. 贵州大学 资源与环境工程学院，贵阳 550003
  • 收稿日期:2010-03-02 出版日期:2011-11-01 发布日期:2011-11-09
  • 作者简介:左双英，女，1977年生，博士，副教授，主要从事地下工程方面的教学和研究工作
  • 基金资助:

    国家自然科学基金重大研究计划(No. 90715042)；国家杰出青年科学基金(No. 50725931)；国家科技支撑计划(No. 2008BAB29B01)

Secondary development and application of an elastoplastic constitutive model based on Zienkiewicz-Pande yield criterion in FLAC3D

ZUO Shuang-ying1, 2，XIAO Ming1，CHEN Jun-tao1   

  1. 1. State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China; 2. College of Resources and Environment Engineering, Guizhou University, Guiyang 550003, China
  • Received:2010-03-02 Online:2011-11-01 Published:2011-11-09

PDF (Chinese)

2015

摘要： Zienkiewicz-Pande屈服准则是Mohr-Coulomb准则的改进，在p-q子午面和π平面上都是光滑曲线，不存在尖点，不仅在数值迭代计算过程中易于处理，而且在一定程度上考虑了屈服曲线与静水压力的关系以及中主应力? 2对屈服的影响。该准则已应用于多个水电站地下洞室开挖支护计算中，但并未被大型数值计算软件所采用。根据弹塑性力学原理，详细推导了基于Zienkiewicz-Pande屈服准则的弹塑性本构模型在FLAC3D中的增量迭代计算格式；在FLAC3D提供的二次开发程序接口基础上，利用VC++程序语言实现了计算过程，并编译成动态链接库文件进行加载和调用。将自编的本构模型对一理想圆形隧洞进行模拟，计算结果与Mohr-Coulomb准则一致，验证了计算格式的正确性。将此本构模型应用于一大型地下洞室群的开挖围岩稳定分析中，得到的塑性区、位移场、应力场分布与Mohr-Coulomb准则计算结果分布规律相符，进一步验证了模型的适用性和可靠性。

关键词: Zienkiewicz-Pande准则, 本构模型, FLAC3D, 二次开发

Abstract: An elastoplastic constitutive model based on Zienkiewicz-Pande yield criterion is developed secondarily in FLAC3D and applied to stability analysis of surrounding rock during the process of excavation in a certain underground hydropower station. Zienkiewicz-Pande yield criterion is the one of forms of improved Mohr-Coulomb criterion，which is a smooth curve and has neither the singular vertex in the p-q plane nor the ridge corners in the π-plane. It is easy to handle not only in the numerical iterative calculation process, but also taking into account the relationship between yield curve and the hydrostatic pressure to some extent, as well as the impact role of the intermediate principal stress ? 2 on the yield surface. The criterion has been applied to calculation for excavation and supporting in several large-scale underground caverns of hydropower stations, but has not been adopted by any large numerical software. According to elastoplastic mechanics, incremental iterative calculation format of elastoplastic constitutive model on Zienkiewicz-Pande yield criterion is derived in detail. Based on the program interface of secondary development provided by FLAC3D, VC + + programming language is used to realize the calculation process and compile into dynamic-link library files loaded and invoked. In simulating to an ideal circular tunnel, the self-defined model is consistent with Mohr-Coulomb criterion to verify the correctness of the calculation format. Finally, this constitutive model is applied to stability analysis for excavated surrounding rock in large-scale underground caverns; and the distribution rules of calculated plastic zone, displacement field, stress field are accord with Mohr-Coulomb criterion, which validates the applicability and reliability of the model further.

Key words: Zienkiewicz-Pande yield criterion, constitutive model, FLAC3D, secondary development

中图分类号: 

  • TU 45
[1] 金青, 王艺霖, 崔新壮, 王成军 , 张珂, 刘正银, . 拉拔作用下土工合成材料在风化料-废弃轮胎 橡胶颗粒轻质土中的变形行为研究[J]. 岩土力学, 2020, 41(2): 408-418.
[2] 邓子千, 陈嘉帅, 王建伟, 刘小文, . 基于SFG模型的统一屈服面本构模型与试验研究[J]. 岩土力学, 2020, 41(2): 527-534.
[3] 李潇旋, 李涛, 彭丽云, . 控制吸力循环荷载下非饱和黏性土 的弹塑性双面模型[J]. 岩土力学, 2020, 41(2): 552-560.
[4] 程昊, 唐辉明, 吴琼, 雷国平. 一种考虑水力滞回效应的非饱和土弹塑性扩展 剑桥本构模型显式算法有限元实现[J]. 岩土力学, 2020, 41(2): 676-686.
[5] 何鹏飞, 马巍, 穆彦虎, 黄永庭, 董建华, . 黄土−砂浆块界面剪切特性试验及本构模型研究[J]. 岩土力学, 2019, 40(S1): 82-90.
[6] 杨文波, 邹涛, 涂玖林, 谷笑旭, 刘雨辰, 晏启祥, 何川. 高速列车振动荷载作用下马蹄形断面隧 道动力响应特性分析[J]. 岩土力学, 2019, 40(9): 3635-3644.
[7] 张超, 杨期君, 曹文贵. 考虑峰值后区应力跌落速率的 脆岩损伤本构模型研究[J]. 岩土力学, 2019, 40(8): 3099-3106.
[8] 刘斯宏, 沈超敏, 毛航宇, 孙 屹. 堆石料状态相关弹塑性本构模型[J]. 岩土力学, 2019, 40(8): 2891-2898.
[9] 张凌凯, 王睿, 张建民, 唐新军, . 考虑颗粒破碎效应的堆石料静动力本构模型[J]. 岩土力学, 2019, 40(7): 2547-2554.
[10] 穆锐, 浦少云, 黄质宏, 李永辉, 郑培鑫, 刘 旸, 刘 泽, 郑红超, . 土岩组合岩体中抗拔桩极限承载力的确定[J]. 岩土力学, 2019, 40(7): 2825-2837.
[11] 汪俊敏, 熊勇林, 杨骐莱, 桑琴扬, 黄强. 不饱和土动弹塑性本构模型研究[J]. 岩土力学, 2019, 40(6): 2323-2331.
[12] 邹佑学, 王睿, 张建民, . 可液化场地碎石桩复合地基地震动力响应分析[J]. 岩土力学, 2019, 40(6): 2443-2455.
[13] 汪 杰, 宋卫东, 谭玉叶, 付建新, 曹 帅, . 水平分层胶结充填体损伤本构模型及强度准则[J]. 岩土力学, 2019, 40(5): 1731-1739.
[14] 孙逸飞, 陈 成, . 无状态变量的状态依赖剪胀方程及其本构模型[J]. 岩土力学, 2019, 40(5): 1813-1822.
[15] 杨骐莱, 熊勇林, 张 升, 刘干斌, 郑荣跃, 张 锋, . 考虑温度影响的软岩弹塑性本构模型[J]. 岩土力学, 2019, 40(5): 1898-1906.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
[1] 张文杰，陈云敏. 垃圾填埋场抽水试验及降水方案设计[J]. , 2010, 31(1): 211 -215 .
[2] 宫伟力，安里千，赵海燕，毛灵涛. 基于图像描述的煤岩裂隙CT图像多尺度特征[J]. , 2010, 31(2): 371 -376 .
[3] 万 智，董 辉，刘宝琛. 基于正交设计下SVM滑坡变形时序回归预测的超参数选择[J]. , 2010, 31(2): 503 -508 .
[4] 孙曦源，栾茂田，唐小微. 饱和软黏土地基中桶形基础水平承载力研究[J]. , 2010, 31(2): 667 -672 .
[5] 王明年，郭 军，罗禄森，喻 渝，杨建民，谭忠盛. 高速铁路大断面黄土隧道深浅埋分界深度研究[J]. , 2010, 31(4): 1157 -1162 .
[6] 胡勇刚，罗 强，张 良，黄 晶，陈亚美. 基于离心模型试验的水泥土搅拌法加固斜坡软弱土地基变形特性分析[J]. , 2010, 31(7): 2207 -2213 .
[7] 谭峰屹，姜志全，李仲秋，颜惠和. 附加质量法在昆明新机场填料压实密度检测中的应用研究[J]. , 2010, 31(7): 2214 -2218 .
[8] 柴 波，殷坤龙，肖拥军. 巴东新城区库岸斜坡软弱带特征[J]. , 2010, 31(8): 2501 -2506 .
[9] 王光进，杨春和，张 超，马洪岭，孔祥云，侯克鹏. 超高排土场的粒径分级及其边坡稳定性分析研究[J]. , 2011, 32(3): 905 -913 .
[10] 胡海军，蒋明镜，赵 涛，彭建兵，李 红. 制样方法对重塑黄土单轴抗拉强度影响的初探[J]. , 2009, 30(S2): 196 -199 .
版权所有 © 《岩土力学》编辑部
地址：湖北武汉武昌小洪山中国科学院武汉岩土力学研究所(430071)　邮编：200042　电话：027-87198484, 87199252, 87199253, 87198752　E-mail: ytlx@whrsm.ac.cn
本系统由北京玛格泰克科技发展有限公司设计开发