裂隙岩体变渗透系数流固耦合分析.PDFVIP

裂隙岩体变渗透系数流固耦合分析 1 2 1 陈敬运 ，刘仲秋 ，刘明 1 河海大学工程力学系，南京（210098 ） 2 山东农业大学水利土木工程学院，山东泰安（271018 ） E-mail ：hhu-cjy@ 摘 要：本文将裂隙岩体等效为连续多孔介质，将岩体渗透系数考虑为应力和应变的函数， 利用ABAQUS 的二次开发接口编制了相应的子程序，实现了渗透系数的动态变化，相对合 理地模拟渗透系数的演化规律，较好的反映了渗流-应力耦合中的“ 间接耦合”效应，并结合 该子程序讨论了隧洞周边围岩在渗流-应力耦合机制作用下渗透系数不变与变化两种工况的 耦合响应。 关键词：渗透系数；渗流-应力耦合；弹脆塑性；ABAQUS ；二次开发 中图分类号：TU45 1 引 言 裂隙岩体渗流研究主要是考虑流体在岩石裂隙中的运动规律，尤其是流体渗流和岩体应 库 力之间的相互耦合作用机理。随着工程应用的发展，渗流场或应力场—渗流场的耦合研究， 已成为水利工程以及其它涉及地下水渗流的工程设计中热点问题。 [1] 例 渗流应力耦合效应分为直接耦合和间接耦合 ，直接耦合是指力学变形和流体渗流之间 的相互作用，其实质就是岩体的缝隙（孔隙与裂隙）体积的变化和缝隙水之间的相互影响； 案 间接耦合是指由于岩体变形和渗流的影响，岩体的水—力性质发生相互改变，即当岩体的缝 E 隙体积改变后，它的力学和水力性质也会发生变化，比如渗透系数会随着岩体的变形和应力 A 状态而变化。早期的研究多把岩石中水的渗透系数或渗透率定义为常数，在渗流场或应力 C [2] 。事实上在各类地下工程中，对 场—渗流场的耦合分析中，渗透系数也常被当作常数处理 于低渗透率岩体，工程可能是稳定的，直接耦合作用更明显；但对于裂隙岩体，在应力作用 模 下岩体渗透率可能成数十倍增加，工程的稳定性就可能出问题，间接耦合作用则显得更重要。 因此，岩体渗透性规律的研究是渗流应力耦合研究中的一个重要课题。 凯 ABAQUS 是由美国HKS 公司研制的目前国际上最为先进的大型通用有限元计算分析软 O 件之一，其非线性（材料、几何、接触）力学分析功能处于世界领先水平。本文将裂隙介质 岩体等效为连续多空介质，利用ABAQUS 提供的二次开发用户子程序接口，尝试开发了岩 E M 体变渗透系数渗流耦合模型，选择深埋引水隧洞周边围岩作为研究对象，开展了考虑渗透系 A 数动态演化的隧洞围岩水力耦合响应研究，以期指导工程实践。 C 2 渗流-应力耦合的等效连续介质模型 等效连续介质模型是渗流–应力耦合问题数值模拟中常采用的模型，该模型将裂隙岩体 等效为非均质各向异性的连续孔隙介质进行模拟，用连续介质理论描述渗流方程。耦合模型 的基本方程大多数采用经典多孔介质力学的Biot 非稳定饱和渗流耦合方程[3]，渗流-应力的 基本耦合方程可以表示为： 平衡方程： σ +ρX 0 (i, j 1,2,3) （1） ij , j j 几何方程： ε (1/ 2)(u =+u ) ε ε +ε +ε