离散流变模型.doc

离散元构造面流变模型2004年5月31日 摘 要 建立适用于离散元方法的岩体构造面流变模型,在传统的离散元基本框架的基础上,引入积分型流变本构方程和粘性等效接触力增量并嵌入到离散元的计算流程中,成功模拟了离散介质结构面的流变性质。通过对一个简单算例计算值和理论解的对比,验证了本模型有良好的精度,具有应用于岩质高边坡和离散块体系统流变分析的重要前景。关键词 离散单元法,流变,接触离散单元法(DEM)是一种适用于模拟离散介质的数值方法。自Cundall[1]于70年代提出以来,离散元方法在岩石介质数值分析中得到了迅速的发展[2,3],它与有限元、DDA等方法一起,成为岩体变形与稳定分析的有力工具。由于离散元方法既能模拟接触面的大变形,又能模拟块体内部的连续变形,它特别适合于模拟被节理、断层等构造面切割而成的不连续岩体介质。大量工程实践表明,岩体系统的变形和稳定问题与时间因素密切相关,在长期荷载持续作用下,其流变性态即使在不太高的应力水平下也十分明显,甚至成为最终变形破坏的主导因素。因此,在进行岩质块体结构的变形和稳定分析时,考虑介质的流变属性是十分必要的。当岩块本身坚硬完整时,其时效变形主要体现在节理、裂隙等构造面的流变行为上。离散单元法由于计算时步很小,计算中所施加的阻尼主要是加快收敛速度,不能直接用于考虑较长时间效应的流变计算,作为采用离散元模型分析岩体流变效应的初步成果,本文建立了分析构造面流变的可变形体离散元数值模型。,又能同时模拟接触面的张开和滑移等力学行为,由于采用显式时步迭代,对应于应力增量,本构方程采用增量形式,离散元计算方法对材料非线形本构关系和大变形问题具有内在的适应性。如图1(a)所示,岩体被裂隙、断层等构造面切割成块体系统。构造面采用法向和切向弹簧来模拟,构造面的接触力Fc可以根据构造面的法向和切向变形,通过接触弹簧的力-变形关系来确定。由于在不平衡力作用下各个块体不断运动,这些构造面上的接触力的大小、位置、方向均在不断变化,需要在计算过程中不断搜索、判断,这是离散元法的重要技术,正在不断地改进中[1],具体内容可以参见文献[1,3]。可变形体离散元法中,每个块体本身为可变形体,它的变形通过将块体内部划分成为三角形常应变差分单元模拟,如图1(b)所示,网格节点N为差分单元B1,B2…Bn的公共节点。三角形差分单元把惯性质量平均分配给其三个节点,即把虚线所围起的多边形的质量分配给N节点。因为每个三角形差分单元的应力、应变均为常量,可以由三角形各节点的位移和块体本身的本构关系来确定,再通过积分可以得到节点弹性力。离散元计算采用显式时步迭代的方法,静力问题亦采用类似动力分析的松弛迭代求解,阻尼力Fd作为一种吸收能量加速系统收敛的手段被引入方程。离散元系统每个网格节点都应满足Newton第二定律,即 m=Fe+Fc+Fd+Fb+Fg (1)其中m为网格节点所分配的质量, 是该节点的加速度。Fe,Fc,Fd,Fb,Fg分别为作用在该节点上的弹性力,接触力,阻尼力,给定的外力和重力。离散元法用显式时步迭代进行平衡计算,对于每一计算时步,将对所有节点进行循环,计算步骤(如图2虚框所