拉格朗日元法及其应用软件FLAC.pdf

维普资讯 。 一 磴拉格格朗朗日日元元法法及及其其应应用用软软件件FLAC 黑龙江矿院 幽 连金 李风仪 ／U J 黑龙江煤校 曼 丁／／4 J 摘要 拉格朗日元涤是一种新的敷值分析方法，特刺适合于求解非践性的太变彤问题， 在岩土力学中有重要的应用。本文舟绍了拉格朗日元岳的原理，应用软件的主要功能和特点。 关糊堂盟垂璧j ，皱形芝兰生 膨闻锨 本文介绍一种新的数值分析方法——拉格朗 日元法 (Lagrangianelementmethod)，它很适 合于求解非线性的大变形问题，在岩土力学中将有广泛的应用前景。结合一些典型的工程问 题．用此方法进行了分析和计算，结果与理论解和实测均吻合较好。 1 拉格朗日元法 拉格朗日元法源于流体力学。从所周知，在流体力学中，研究流体质点运动的方法有两 种，一种是定点观察法，亦称欧拉法 另一种是随机观察法，称为拉格朗日法。后者是研究 每个流体质点随时间而变化的状态，即研究某一流体质点在任一段时间内的运动轨迹、速度、 压力等特征。把拉格朗日法移值到固体力学中，把所研究的区域划分成网格，其结点就相当 于流体质点，然后按时步用拉格朗日法来研究网格结点的运动，这种方法就是拉格朗日元法。 它采用按时步的动力松驰进行求解，这与离散元法相同，求解时基于显示差分法，不需形成 刚度矩阵，不用求解大型方程组，因此，占用内存少，求解速度快，便于用微机求解较大规 模的工程问题。 拉格朗日元法采用差分方法求解，因此首先要将求 lr l 解的区域划分成四边形的网格，在边界和巷道周围等不 规则处也可用三角形网格拟合。图1是在大倾角煤层中 开掘的巷道状况，注意煤层、底板和顶板是不同的。 拉格朗 目元法的计算循环如图2所示，假定某一时 刻各个节点的速度为已知，则根据高斯定理可求得单元 的应变率，进而根据材料的本构关系就要求各单元的新 应力，进入下一个计算循环。 应变张量由增量形式表示为 ： ．S■ = 专[差一尝]△ () ’一’t ‘z‘ 式中 一 应为增量的张量，i， 1，2； ．——节点 图1 阿格的划分 的速度分量； 一 节点的坐标；＆一 时步 为提高求解的精度，一个四边形以左右两条对角线将之分为四个三角形 (图3a和3b中的 a、b、c和d)，每个三角形假定为常应变，于是四边形的应变为此四个三角形应变的平均值 ·224·1997 №3一 ‘矿 山压力与顶板菅理 维普资讯 根据高斯定理 ，对于函数 f有 ： 对每十单元 』 一』A ㈤ 式中 A～ 单元的面积{s——周边； ——外法线 的方向余弦。 则 一 (3) 由此可求得其它分量的值，将这些值代八 (1)式即 可求得应变增量，于是可以根据材料的本构关系求 图2 拉格朗日元法的计算循环 得应力增量为： N 由u一 ／( ，J|…) (4) 式中f表示本构关系的函 数，它与应变增量、原有的垒应 e w 力以及材料常数等有关。 作用在各节点的不平衡力R 4 3 4 3 可取