具有任意形状单元离散模型的界面元法.pdfVIP

第 15 卷第 2 期 工 程 力 学 Vol. 15 No.2 1998 年 5 月 ENGINEERING MECHANICS May. 1998 方义琳 卓家寿 章 青 ( 河海大学土木学院, 南京 210098) 文[ 1] 推出了求解 连续介质体的界面应力元法, 并在理论和实例验算两方面论 证了该法具有精度好、效率高、适用范围广等优点。通过进一步研究, 该法还可适应任 意网络离 散模型的计算。本文引用 自然坐标, 实施了任意形状单元 的界面元算法, 并给 出了算例, 成果表 明,本文的算法可行, 其精度也是令人满意的。 界面元, 自然坐标, 面积坐标,体积坐标, 坐标转换 、 众所周知, 在有限元位移法中离散模型单元形状是受限制的, 即使采用 同形状单元组 合模型时, 也有如下几点 便之处: ( 1) 对 同的单元必须采用 同的位移模式及坐标变换 式; ( 2) 考虑网格局部加密, 采用形状相同而大小 同的单元时也会遇到在单元连接方面较 大的困难; (3) 采用任意曲面的等参单元时, 要使坐标变换一一对应, 必要条件是单元 能 过于歪斜或出现凹角。这些问题的存在给应用带来了许多 便, 而用本文界面元法却可以 很方便地解决。 文[ 1] 提出的界面应力元法是在 Ka ai 求解均质、线性静力问题的刚体弹簧元的基础 上, 以反映变形特性的界面元件( 如弹性元件, 塑性元件, 粘性元件和接触元件) 取代弹簧而 建立的适用于求解非均匀、各向异性、非线性的 连续介质静动力问题的一种新模型。该模 式实质上为基于分片的刚体位移场、积累变形于单元界面的非协调位移元, 以分块单元交界 面上的界面应力表征结构的应力场, 由相邻单元的相对变形和 同材料特性得出界面应力 公式, 根据广义变分原理导出求解结构位移场的支配方程, 最终得到结构的变形( 含界面上 本文收稿 日期: 1997 年 3 月 28 工 程 力 学 相对变形) 和众多离散界面的应力以及单元内的表征应力。 界面元法具有很多优点, 如: ( 1) 支配方程带宽小, 耗时少, 效率高; ( 2) 界面应力依赖 于联接相邻单元微分条的相对变形量, 故其应力精度 低于位移精度, 从而提高了应力状态 判据的可靠性, 使其非线性解 致出现漂移; ( 3) 由于它是非协调元, 可在一定程度内消除 位移协调元模型太刚的弊病, 有望获得较高的位移精度; (4) 可以描述各界面的工作性态, 如是连续的或开裂的或滑移的; (5) 由于其劲度矩阵[ K ] 是通过界面积分形成的, 因此其主 要的计算依赖于每个界面, 与分块单元的形态无等等。本文着重讨论第 5 点的实施方法, 通 过对各类型块体单元采用相同的位移模式, 坐标变换采用自然坐标表示, 给出考虑任意形状 单元的界面元法的算法和考例。 、 从最小势能原理出发, 放松单元交界面的位移连续条件, 引用 Lagrange 乘子法得到分 区 完全的广义变分原理, 即真解 u 和 T 使泛函 为驻值, 即 i i mp = [ A ( e ) d- f u d- pu ds - T u ds] =