岩石水力劈裂有限元模拟.pdfVIP

4J0 以科学发展观促进科技剖新(上) 5徐兴永．石学法，于洪军，李萍崂山顶、涧、沟、坡、麓、滩、岬一带巨砾成因研究[J]．海洋科学，2004， 28(6)：10～13 6赵松龄，张宏才北京西山灵岳寺古冰川沉积物的初步研究中国更新世冰川与更新世地质论文集 (一)．北京：地质出版社，1979．69～72 7施雅风，崔之久，牵吉均等．中国东部更新世冰Jll与环境问题、北京：科学出版社，1989 8朱诚等．华东山地更新世沉积环境研究南京：南京大学出版社，2000，1～354 9李乃胜，石学法，赵松龄，于洪军崂山地质与古冰川研究．2003，北京：海洋出版社 10徐馨，沈志迭．对长江中下游更新世冰川发育的新认识．贵州师范大学学报(自然科学版)，1999 (1)：1～6 xuxingyong@ms．qdio．ac．cn，xuxingyong@mail．rio．gov．cn。 岩石水力劈裂有限元模拟① 谢兴华张瑞凯 江苏省南京市广州路223号，南京水利科学研究院，210029 摘要本文将细观的观点引入岩石工程问题的力学模拟。在应力集中的局部单元引八 断裂力学的应力强度因子判据判断是否达到开裂条件。若是已经开裂，根据能量释放率最大 的方向就是裂纹开裂方向的结论．确定裂纹的扩展方向。已经达到开裂条件的单元．做开裂标 记，认为单元位置成为裂隙，已经失去抗拉、剪的能力。随着新裂隙的产生，材料的渗透性套发 生变化。渗透系数的变化与村料变形有关，损伤变量描述材料变形对力学性质的损伤。本文 以损伤变量作为桥梁，在材料变形与渗透性变化之间建立联系，作为渗流应力耦合的纽带。最 后文章通过算倒模拟在水压力驻动下的裂纹扩展过程，验证了本文模型的正确性。 关键词 岩土工程 水力劈裂损伤 裂纹扩展 水力劈裂问题的数值模拟，需要解决这样几个问题：水的作用方式、开裂的起始点、裂缝的扩展以及 水的渗透与裂缝开裂的耦合关系。 关于水的作用方式，在渗流应力耦合的文献中，多为基于拟连续介质假设下，作为体力在应力场中 叠加，相当于孔隙压力的作用。实际上，水对岩体的作用不仅仅是孔隙压力的作用，还有面力作用和渗 透力(水头梯度)的作用。面力可以作为力学计算的载荷，渗透力(水头梯度)作用在数值分析时的处理 还未见到相关的报道。 开裂的起始点问题。岩体是非均质体，在力的作用下，裂纹首先从弱点开始起裂。如何在拟连续假 设之下模拟这种非均质性是一个困难的问题。文引入概率分布的方法来描述岩石的非均匀性。认为岩 石内部各点的力学性质(比如弹模)是随机变量，在统计概念之下服从一定的分布(比如Weib。11分布)， 并且在此基础上开发了计算程序RFPA”，计算结果接近试样破坏状态。 ①本文受到南京水利科学研究院青年科学基金资助(Y10402)。 以科学发展观促进科技创新(二)411 许多研究者曾经开发过一些二维、假三维或全三维的数值工具用来模拟水力劈裂。康奈尔大学综 Element 合FRANSYS(Fracture Syslem)两个程序开发了用于水力劈 AnalysisSystem)和BES(Boundary 裂分析的三维程序HYFSYS(HydraulicFracturing 单驱动的、图形化分析软件。它采用拓扑数据结构，初始条件、边界条件、材料属性等几何描述和模拟属 domain)、三维直接边 性与拓扑属性相关联，裂缝的开裂要重新划分网格。BES是一个单分析域(Single 界元方法。采用裂隙流与固体骨架耦合的分析方法，裂隙网络的变化要靠网格的重新划分，这样计算工 作量很大，限制了计算规模，不适合计算大型区域的水力劈裂问题。这也是有限元方法模拟裂纹开裂问 题面临的一大挑战。非连续变形(DDA)数值方法模拟裂纹开裂有一定优势。近十多年发展起来的数 Free 值流彤方法、无单元方法(ElementMethod)分析水力劈裂问题也具有相当的优势。研究开