3DEC中节理法向循环加卸载模型开发.docVIP

3DEC中节理法向循环加卸载模型开发 3DEC中节理法向循环加卸载模型开发 关 莉 莉，张 雪 锋 (黄河水利职业技术学院，河南 开封 475004) 摘要：节理在循环加卸载作用下，其法向变形表现出了很强的非线性，Bandis根据室内试验成果，提出了用双曲线方程来表征节理法向应力与闭合量的非线性关系。鉴于目前3DEC程序在模拟节理非线性变形特征中的不足，结合3DEC提供的接口，对Bandis节理模型进行了二次开发，并给出了运用C++开发3DEC节理模型的步骤。同时，通过内部FISH语言开展了节理在循环加卸载条件下的单轴压缩数值试验，并与室内试验及理论结果进行对比。结果表明：开发的Bandis模型能较好地表征出循环加卸载条件下节理法向应力与法向闭合量之间的非线性关系；数值试验计算结果与理论值完全吻合，验证了所开发模型的正确性。研究成果为循环荷载下岩质边坡的稳定分析提供了基础。 关 键 词：节理; 法向循环加卸载; 非线性变形； Bandis节理模型； 3DEC程序 1 研究背景 节理岩体是不连续介质，较完整岩石其力学特性相对复杂，力学性质具有非线性。在水利枢纽工程、边坡稳定及采矿工程中，大多涉及到了节理岩体的变形稳定性分析。在高应力及循环荷载等复杂环境下，节理的力学性更为复杂。节理及其力学性质影响着岩体中的应力分布，并显著影响整个岩体工程的稳定性。节理在循环荷载下的非线性变形特征成为了节理岩体变形稳定分析的重点。 美国ITASCA公司开发了三维离散单元软件3DEC(3 Dimensional Distinct Element Code)。该软件是基于Cundall于1971年提出来的离散元方法开发的，允许模型中不同的块体之间在计算时产生大位移、旋转，甚至块体之间完全脱离，并能自动判断块体之间的接触状态变化[1]。因此，在岩体边坡滑坡、煤矿开采、水利枢纽等工程的岩体稳定性分析中，3DEC有着独特的应用优势。目前，在3DEC中节理模型主要有库仑滑移节理模型和由Cundall、Hart提出的连续屈服节理模型[2]。但是，这两种模型均存在着一些不足：库仑滑移节理模型无法合理地模拟出循环加卸载下的节理法向变形特征；库仑滑移节理模型中法向变形为线弹性，无法合理地模拟节理的非线性力学特性。连续屈服节理模型采用指数方程来表征节理模型法向变形特征，其无法合理地模拟出循环加卸载下的节理法向变形特征；且连续屈服节理模型中的一些参数难以确定，需要依靠经验来确定。 从唯象的角度出发，国内外学者通过对室内试验结果进行数据拟合，提出了多种节理变形曲线，主要有：用双曲线方程表征的Goodman模型[3]、Bandis模型[4]，用对数方程表征的Shehate模型[5]，用指数方程表征的经典指数模型、Malama-Kulatilake模型[6]、周建民等[7]、俞缙等[8]修正的指数模型等节理法向模型。目前，Bandis提出的双曲线模型应用较为广泛。 3DEC本构模型的主要作用是根据给出的应变增量按照相应本构方程计算出新的应力增量，进而得到新的应力。目前3DEC为用户提供了两种方法进行本构模型开发：第一种方法是在原有本构模型基础上，通过3DEC内部FISH语言命令在3DEC计算过程中对相应的参数进行动态的修改；第二种方法是通过3DEC为用户提供的接口，利用C++语言编写3DEC程序可以直接加载调用的DLL文件(动态链接库文件)，3DEC程序自带的本构模型均是DLL文件的形式。通过FISH在3DEC原有的本构模型进行编写比较方便，但其计算速度比较慢，直接影响计算效率。吴顺川等[9]根据Bandis等[4]的室内法向循环加卸载试验数据，通过3DEC中的FISH语言动态修正3DEC内置的连续屈服节理模型的参数，对3DEC提供的模型进行开发。把节理法向循环加卸载下非线性变形特征通过C++编写相应的节理模型，并将其引入到3DEC中，能更好地模拟出循环加卸载条件下节理非线性变形特征。 为了弥补3DEC中节理本构模型的不足，引入了Bandis关于节理的法向循环加载试验研究成果，基于3DEC提供的User-defined joint constitutive models 接口，直接采用C++语言编写了相应的节理本构模型的 DLL文件(动态链接库文件)。基于编写的jmodelbandis节理本构模型进行了法向循环加载试验，并与室内试验对比，验证了该模型的正确性。同时，实现了将建立的节理法向循环加卸载非线性变形本构模型内嵌于三维离散元软件3DEC中，为后续研究循环荷载作用下岩质边坡工程变形稳定提供了基础。 2 节理本构模型确定 国内外许多学者对节理法向变形规律进行过研究，主要是根据室内试验得出岩体和相应完整岩石的法向应力-法向变形曲线(见图1)。其中，岩体是由节理和岩石组