大位移分析中DDA%2f流形元计算模型的研究.pdfVIP

大位移分析qaDOA／流形元计算模型研究 朱以丈周颖徐金虎 (武汉水利电力大学) 【提要】奉文分析了非连续变形分析(DDA)／流形』亡法的计算援型在大变形问题中两类误差，并 给出了改进办法，用新的计算模型的算例与实验结果较接近。 实际岩体中含有软弱结构面，它们往往是岩十结构稳定性的控制因素。在开挖、蓄 水、地震和爆破等工作条件下，岩士结构经常会进入局部或瞬态的人变形状态，如滑 移、开裂和脱离，甚至导致失去稳定，进入运动状态，如整体滑动、倾倒和下坍。结构 面具有强烈的不连续性，建立这类介质的计算模型，尤其是在大变形问题中的计算模型 是岩体结构工程中迫切需要解决的问题。这已经引起许多学者的关注，提出了多种计算 模型和数值分析方法，如节理单元法、等效连续体模型、损伤模型、离散单元法、块体 理论、刚性元法和不连续分析(DDA)等。近年来，又在块体理论}flDDA法的基础上发展 出流形元方法“J。目前，DDA／流形元法取得较大进展。在j：程中开始得到应『{；j。本文讨 论大变形分析中的DDA／流形元计算模型的若干问题。 一、刚体转动项的线性化误差及其解决办法 在DDA／流形元方法中所建立的块体单元位移函数如下[1m]： 在每一时间步或荷载增量步内，第，价块体内的位移为(为简明起见，以平面问题为 例) 0(少 ，O 眇 ” j卜 “ 喘砧．协 ) 扛 蛳¨荫鼻毋～ ‰，v。，酿为块体单元形心处之刚体位移。 毒 占，，占。y。，为块体单元之常应变分姑。 苎 为了提高求解精度，不少学者对上述线性位移模式进行了改进，典型的工作如二 阶、三阶多项式的位移模式懵』，及Fourier级数作为位移模式…。值得注意的是，无论怎 样提高位移模式之阶数，‰。vo，oo，‘，占。，y。这六个量总是不可缺少的，而且 由于DDA／流形元法主要用下岩体1：程分析，所谓人变形主要是指火位移小鹿变类型．上 述这六个位移分量起着重要作用，在麻崩丁1人变形问题时，必须解决其刚体转动项玩线 ·298· 性化引起的误差。 刚体转动皖造成的水平位移和垂直位移的线性化的结果。 一(XI，，】)．令以p尼则： cz， {：：}=RJlc。oins。(aa十+岛00，)～-。cions口a} 肖吼值很小时，N令-cosOo 力 ㈦=尺雠枷熙(Y-¨Yo)严1． 线性化后，，谁不是到达尸1，而是剑达岛(见图2)， 而 PB=OP‘喀岛=RtgO()z R‰。月岛 所以 AP=ROo COSa=“自． 4咒=ROo sina=vn 从而我仃J可以看到，每转过…个00所导致的线膨胀误著