二维离散单元法中流量守恒方程的求解方法.pdfVIP

二维离散单元法中流量 守恒方程的求解方法函 肖裕行1 王泳喜2 王思敬。 杨天鸿3张寿全1 (2东北大学f：七昭110005, 3。01 中科院地赢所北京i0002．q) 、：抚’砭搿利分靛抚嗽：I 摘要本文给}【：了二维离散单元{主中块及块体矧空隙的流馈守恒方程的矗式差分格式及其稳 定条件，并扔步建立了T系统水力耦台行为分析的计算循环。这些结茉扦。键丁二维离散单元法樽 拟问题的能力。 关键调水力行为分析 离散单元藏 显式差分 1 引言 离散单元法因能模拟的问题与工程岩体结梅的几何构形和力学行为的相似，自创立以 来得到了进一步发展和深入广泛应用“““，并有丁商业鞠件，如二维中功能强大的UDEC (Universal1)istincl ElementCode)日]，但是它对水力辐台叫题的考虑迈不完善．目前假 UDEC能够处理水沿块体间空隙流动。实际上离散单元法模拟的系统巾块体具有渗透怙的 情况较常见．它可能是具有渗透性的土体孔隙舟贡，或出于可操作性与精度两万酉的考虑。 还可能是裂隙岩体中数量众多的低级别结构面1：力贡献的总体反映。圈此发展二维离散单 元法使之能考虑块体的渗透性是必要的。文献[4]绐出了水力鞘台情况下采统的梳量{j：恒万 程，本文初步给出它们的显式差分格式和求解方‘式。但左交啦F沦垮蛱体视为士{车孔隙介 盾。 2 二维离散单元法基本原理和系统的流量守恒方程 为讨论方便，首先简要介绍二维离散_嘞元法的原理和系统的流量守恒宵程。 离散单元法中系统的几何构形是由块体及块问空隙捣成的系统(图1)，单个块体的变 形特征和运动方式没有人为限制。块俸闻的相互作用力通过它们闯的接触传递，接触由块体 间的相对位置确定(图1中块体元间的短线表示接触)．已有的接触可以脱离，新接触也可以 产生；块闻空隙的几何构形通过接触确定(图1中的充填部分)，它也随块体同相对位黄的变 化而变化。对于刚性块体因自由度有限可直接求解运动方程，对于可变形体是把无限自由度 的求解转化为有限自由度的求解，即：首先把块体剖分成一系列的常应变三角形单元，然后 对每一个节点在它的专有子域中建立它的运动方程(图2)。节点n的专有于城是指：对于以 ① 中国科学院资源与生态环境研究重点项目(KZ952一Jj一416)史持 图1 离散元法中系统的Jl,#1图形 尉2块体的三角形剖分和节点n的专有子城 它为顶点的所有围绕三角形单元，依次把笆们的中心和边上的中点连接起来形成的围绕节 点n的多边形(图2中的充填部分)。系统采用显式差分动态方法求解“。1。 文献口]给出了考虑块体渗透性时系统的流量守恒方程，结果如下： 块体的流量守恒方程是 JKJ7=e，等…e：莩 (1) 速度；口，，以是与介质变形有关的常数， 块闻空隙的流量守恒方程是 W4-Q：怕==扣㈨一半i8／4十e^警(2) 其中，k是流体的体积模量；Ⅳ为源汇项；H，V分别是块间空隙的面积和几何中心的水头} e为块间空隙的充填度，无充填时e一1．充填时e为充填物的孔隙率；Q，为它与其它块闻空 隙问的流量交换；Q。为块闻空隙与相邻块体间的流量交换(图1)，即 Q：；=∑q’Ii 在上述方程中，块问空隙与块体边界面上边界条件的处理方法是(图1)：水头连续；对于流 量交换，块体必第二类边界条件考虑，块间空隙以源汇项考虑．．即 q·nllr-=el，巩·n2I一=E2，……，12。·．．!，=‘ (3) 3块体流量守恒方程的显式差分格式 进行水力耦合行为的数值模拟要求力学方程与流量守恒方程定义在相同的予空间，以 确保耦合行为在空间的一致馅：。为此，两种场需采用相同的三角形制分网络，同时节点的流 量守恒方程也壶在它的专有子域申建立。当不考虑方程(1)右端第二项，即不考虑耦合行为 程(1)右边第二项可看着源汇项，它的出现并没有改变不考虑该项时对应方程的数学结构， 仍为抛物线