平面弹性力学问题的离散元法.pdf

第27卷专辑 固体力学学报 V01．27S．Issue MECHANICA 2006年12月 ACTA SOLIDASINlCA December2006 平面弹性力学问题的离散元法 凌道盛 吕铁墩 何淳健 (浙扛大学建筑工程学院．杭州，$10027) 摘要 根据离散元的基本原理，基于变形体的理论提出了适用于平面弹性力学问题的界面位移、应变和应 力模式．建立了求解平面弹性力学问题的离散元方程和相应的迭代求解方法．通过界面位移可以简洁地将位移和 力的边界条件引入离散系统的控制方程，也可以方便地求解节点位移．数值算例表明，与具有相同网格的有限元结 果相比，离散元能同时给出精度相对较高的应力解和精度相当的位移解． 关键词 离散元法．平面弹性力学同题，界面应力 0 引育 质、不连续介质和两类介质混合结构的数值分析方 大多数实际工程问题的破坏都是一个渐进的过 法．鲍鹏印]等提出了可变形体离散元模型，并对计 程，如边坡的稳定破坏，混凝土构件在外荷载作用下的 算参数的选取进行了分析研究． 材料破坏等．这类问题都有一个共同的特点，即分析对 建立具有一般离散元法的优点且能考虑连续介 象的力学行为从连续介质向非连续介质转化．基于传 质变形特点的离散元分析方法是将离散元法应用于 统的连续介质力学的有限单元法在求解这类问题时存 边坡等破坏过程分析的前提和关键之一．本文从离散 在一些困难，而离散元法则在计算和模拟材科具体的 元的基本理论出发，以单元形心点位移为基本未知 破坏形式和破坏过程方面显示出巨大的生命力． 量，提出了考虑单元体变形的交界面位移和应变模 离散元法是20世纪70年代初由CundaUcl]首式，建立了适用于平面弹性力学问题的离散元方法． 先提出的．最初，它的研究对象主要是岩石、粉体(颗 1 单元位移模式 粒散体)等非连续介质．近年来，离散元法的应用领 如图1示，建立直角坐标系：my，并规定单元节 域扩展到求解连续介质及连续介质向非连续介质转 点按逆时针编号．图1中四边形单元定义为i。i2f；厶， 化的力学问题．刘凯欣、高凌天[2]给予离散元法一 称线段fl如、f2i3、f3厶、il{l为单元界面． 个较为宽松的定义，并进行了比较详细的评述． 离散单元本身一般假定为刚体，单元形心的刚 体位移是离散元法求解的基本未知量，通常也称单 元形心为物理节点．另外，离散单元法本身对单元形 状没有特殊的要求．单元间的相对位移等变形行为 由单元交界面的连接来实现，连接的形式在力学机 理上可分成接触型和连结型两大类．接触型是散体 特有的连接形式，适于分析散体的力学行为．连结型 连接考虑单元间没有缝隙且符合变形协调条件，主 图1 离散元坐标系 要用来处理连续介质力学问题．相比较而言，无论在 理论上还是应用上。连结型连接形式方面的研究还 在平面内取一点G作为单元位移描述的参照 不充分，研究成果相对较少．川井[31最先提出了块 点，G点的三个刚体位移分量为。，％和∞。，三个应 体离散元模型，并命名为刚体一弹簧模型(RBSM)， 变分量为e可、eI和y．．‰，％分别为沿z，y方向的位 钱令希、张雄[．]从虚功原理出发，赋予刚体一弹簧模 移分量，以坐标轴正向为正；埘l为绕G点的转动位 型以理论基础．由于与有限单元法极为相似，也称之