裂隙水流与混凝土开裂相互作用的耦合模型.pdfVIP

水 利 学 报 2007年 12月 SHUILI XUEBAO 第38卷 第12期 文章编号：0559—9350(2007)12—1466—09 裂隙水流与混凝土开裂相互作用的耦合模型 方修君，金 峰 (清华大学水沙科学与水利水电工程国家重点实验室，北京 100084) 摘要：利用扩展有限元法对裂隙水压作用下混凝土结构的开裂过程进行了模拟，裂纹中的流体运动在立方定理的 基础上求解，混凝土开裂力学行为采用了黏聚裂纹模型来描述，采用松散型的耦合方式考虑裂隙水流和开裂结构 间的相互作用。推导了在基于虚节点法的扩展有限元法框架下的耦合模型控制方程，通过双层迭代结构进行总体 方程求解，对混凝土试件在裂纹水压作用下的开裂过程进行了数值模拟。不同加载速率条件下的计算结果反映的 裂纹水压变化规律均与试验相近，结果也说明了扩展有限元法是处理连续一非连续问题的一种有效方法。 关键词：扩展有限元；黏聚裂纹模型；水力劈裂；松散耦合；数值模拟 中图分类号： rv313 文献标识码：A 石油开采过程中，人们常利用水力致裂的办法来提高产量…。水利工程中，随着筑坝技术的提高， 人们已经有能力兴建一批200m甚至300m级的超高坝，对如此高的混凝土坝进行坝体安全评估时，考 虑高压水劈裂的影响是非常重要的 。奥地利的Kolnbrein拱坝在发生坝体开裂时检测到坝基扬压力 水头的上升。还有学者 指出，在动力作用下(如地震荷载)，坝体裂隙中的水压会发生显著变化，产生 附加的上抬力从而对坝体的动力稳定产生影响。在这些有压力水作用的环境中，混凝土裂缝或者岩体 裂隙的失稳破坏与水在其中的流动有密切的关系，开裂过程往往伴随着与裂隙水流的相互耦合作用。 鉴于水力劈裂问题的重要性，国内外许多学者对此进行了一些试验 “ 和数值模型¨ 研究。在一 般的有限元和差分等数值方法的框架下，对非连续界面的描述依赖于网格边界，因而模拟这类问题时需 要预先设置裂纹的路径，或者采用重划分网格的技术，计算量大。在模拟裂隙中水压力分布时，裂隙水 压力常考虑为裂纹张开度的函数 ，忽略裂隙中水流动的时间效应。因此，这些模型一般只适用于静态 或准静态问题。扩展有限元法(Extended Finite Element Method，XFEM) 通过引入非连续位移模式，使 得不连续位移场的描述独立于网格划分，可以在不重新划分网格的前提下，显式地描述裂纹扩展路径， 为描述裂纹扩展与裂纹中水压耦合作用提供了方便。本文在扩展有限元方法的框架下建立裂隙水流与 混凝土结构相互作用的耦合模型，考虑裂纹的扩展以及裂隙中流体流动的时间效应，并通过前人的混凝 土劈裂试验进行模拟验证。 1 扩展有限元法 XFEM是 1999年以来由Belytschko等… 。 提出和发展的一种新兴的数值方法，该方法是基于单位 分解 理论对传统有限元法进行的扩展，其核心思想是通过富集具有特殊性质的函数来增强对待求场 函数的描述能力。一个常见的例子是在XFEM中引入Heaviside函数来构造非连续位移场，可表示为 u ( )=∑Ⅳ ( )u +∑ ( )(ttj+dyH( )) (1) iEl／J l∈3 收稿日期：2007．03．22 基金项目：国家自然科学基金资助项目50539020) 作者简介：方修君(1982一)，男，安徽广德人，博士生，从事水工结构研究。E-mail：farIsxj98@mails．tsinghua．edu．ca — l466 — 式中： ( )表示位置 处的近似场函数值；“ 、 为常规自由度；d 为与富集函数相关的附加自由度； ，为离散域中所有节点的集合；J为影响域被非连续界面穿过的节点的集合；N ( )、Ⅳ，( )分别表示节 点 i、J影响域上的形函数；H( )为Heaviside函数，它的富集使得式(1)具有描述非连续位移场的能力。 在XFEM中，节点影响域由所有共享该节点的单元组成，也即单位分解法中的子域 (subdomains)或 分片(patches)。XFEM