固液相耦合的有限元法及应用.pdfVIP

第19卷第4期 地 F 审 涮 V01l 9 No4 1999年】o月 UNDER(；R()RNDsPACE ()cf999 固液相耦合的有限元法及其应用 盛建龙 刘新波 欧阳j台华 30081) (武汉科技大学，4 State PennsylvaniaUniversity，USA)) [摘要]本文探讨1岩体内固液相耦合的有f暑元分析方法．考虑了地下水孔隙压力与砉体变形之 问的相互耦合作用。建立孔隙水压的计算模型．并计算分析了在重力场和孔隙水压力共同作用下某露无边坡 的应力分布特点。结果表明，固液相耦合的有限无法能反映地下水对岩体工程的应力分市和位移变彤的影 响。 [关键词] 目液杷耦台 有限无法 边垃 l 引 言 大量的岩体工程实例已经证明，地下水魁影响岩体工程稳定性的主要因素之一。地下水侉 流的特性及其对岩体的应力分布、变形特征的影响越来越受到国内外学者的高度重视。由于岩 体内地下水的分布和渗流很下均匀．并且等体布身是·种非均质、不连续的各向异性体，造成 研究地下水渗流问题极其复杂。因此，在芎虑岩体骨架(固体)与渗流(液体)的相互耦合作崩 下．分析岩体工程的应力和变形分布规律具商P分重要的意义。 众所周知，岩体是由岩石和不连续面构成的，根据岩体内裂隙、孔隙的发育程度，一般采用 两类解析方法来模拟岩体的不连续性。一类是采用节理单元模拟规模较大的节理裂隙等不连 续面，但计算机能够模拟的不连续面个数有限；另一‘类是从宏观上采用连续介质力学理论来分 析节理裂隙发育且规模较小、或岩块同体介质匿水性较好的岩体。对于含有大量孔隙的岩体。 因此隙在岩体中所占的体积较小，宏观r电可看成连续介质。 本文将岩体固体骨架变形与地F水7L隙压力相互耦合，建立相应的有限元计算公式井再 来分析岩体的渗流问题。在考虑自重应力场和渗流水压的共同作用下，研究一矿山岩质边坡的 应力分布和变形特征。计算时考虑了边坡体为完全豌于和完全饱和两种情形。 2 固液相耦合的有限元法 对于裂隙发育或透水性较好的岩体，在考虑地下水的影响时，假发孔隙中充满r各向同性 的液体。通过分别建立岩体田体骨架的位移方程和孔液体的渗流方程．井将两者耦合分衬二，得 出多孔隙介质固坡两4H辋台的有限元方程。 2．1岩体固体骨架的位移方程 假设单f_；)=体积弹性多孔介质内充满各向同性的液体，则总应力的分布满足平衡微分方程： 磐一o (。乩2) 324 地下空间 第19卷 式中％为总应力；F，为单位固体的体积力。 根据弹性力学的符号约定，总应力Zi，可表示为： L，一％一pd。 (z，j一1，2) (2) 式中：巩，作为在固体骨架上的有效应力；户为孔隙压力；d，，一1(当i=，时)，屯=0，(当j≠， 时)。 经过时间t后，有效应力为： oi，=日十云， (3) 式中：砚为初始有效应力；毛为有效应力增量。 如果假定孔隙介质各向同性，由线弹陛应力一应变关系得出： 己，一2G∈v-4-一ckja， (4) 其中：i，为固体骨架的应变增量；G和^为拉梅弹性常数，G又称剪切弹性模量， G—E／[z(1+u)]，^一Ev／[(1+口)(1—2v)]，E为弹性模量，w为泊松比。 一般， i。一e。一ej (5) 其中：s：为初始应变，由温度变化等因素产生．在此