有限元分析中的最优离散化有限元分析中的最优离散化.pdf

维普资讯 第23卷 第 3期 力 学 进 展 VO1．23 № 3 1993年 8月25日 ADVANCES IN M ECHANICS Aug． 25， 1993 垄 塑 重虎大学(邮政编码6300,1：4) 0 2、z／f 提要 系统的最优离散化具有很大的实用价值和理论意义．采用这种方法可 以用最 少的 自由度得到最精确的数值解 ．本文给出了两种有代表性的方法，并收集了一些 应用实例，其中包括笔者和刘永成，朱建中在 Zienkiewicz指 导下所作的一些工 作．这种方法作为一种重分析方法不仅可用于线性 问题，而且也可 以用于具有表面 非线性的弹性接触问题． 关键词 立坠垂 ±垂优离墼 叠层单元 1 引 言 在有限元分析中，已经广泛采用各种网格 自动生成方法．然而这类方法一般只能给出均 匀分布或逐步加密的有限元网格，常常不能反映出不同求解区域能量密度的明显差异．从工 程分析的需要出发，总希望能在能量密度大的区域局部加密网格或局部增加高阶 自由度，但 不必重新划分全部网格，也不需要逐步过渡．这样的问题不仅具有很大的实用价值，而且具 有重要 的理论意义．因为从数学上说，有限元法是根据变分原理求解数学物理 问题的一种数 值方法，归根到底是一个泛函极值 问题 ．如果人为地任意划分网格，设置节点 自由度，则只 能在既定的离散系统下变动场函数在各点上的值，使其能量泛函达到极值．而最优离散化则 不仅使场函数 (节点变量)的值可变，而且使节点坐标、节点交量的类型及节点 自由度数在 一 定范围内可变，从而更精确地求得极值．这是一个高度非线性的极值 问题． McNeice和Marcal“提出了最优离散化命题．Babuska 首先在这方面提出了实用方 法．Turcke和 McNeice ，Pedersonn。。，Carroll ¨等也研究了一些经典问题 阻引． 最优离散化一开始即沿着两个方 向发展 ：局部加密网格即 收敛方法}局部增加高阶 自 由度即 收敛方法．前者的代表是Peano提出的分级节点来实现单元增殖的方法，后者的代 表是Irons提出的叠层单元 (hierarchicalelement)方法．但70年代 中期提出的这些方法当 时并未获得广泛应用．主要原因是难于判定在什么区域需要加密网格或增加高阶 自由度 ． O1iveira“ 首先提出在应变能密度变化最大的区域增设单元 (或增设高阶 自由度)的 方法，即首先按照初始 网格，求出场函数节点值，根据其分布情况计算其应变能密度，并绘 制出等能量线，从而判定哪些区域需要加密网格．有时甚至可 以利用这种等能量线及其正交 ·336 · 维普资讯 线族 自动形成新的单元．Molinari和 Vivieni… 认为这种单元就是近似的最优离散化的有 限元．Melosh和 Marcal 认为比较简单的方法是直接在能量密度大的地方加密网格．也 有人主张用主应力线、最大剪应力线等来判断．然而所有这些方法都必须经过多 次反 复 试 算，因而存在计算工作量太大的缺点．而且这种方法只适合于仅有一种外载的情 况． ，对 于组合加载的情况就很难确定网格加密区域． 同时，采用这类方法，在 同格加密 以后，一般 要全部或局部重算刚度矩阵，重新求解新的方程组．由于这些原因，网格局部加密和其他最 优离散化方法没有得到广泛应用． 这些方法在8O年代重新引起人们重视并成为热门话题，其原因是出现了有限元误差估计 方法．根据这种方法，在经过l一2次试算后即能找到局部加密网格或者局部增加高阶 自由度 的求解区域，从而有条件地实现最优离散化．在这方面，zienkiewicZ，Babuska，Keiiy和 Gago等人作 了许多工作 ．，‘”． 笔者及朱建中、刘永成在ZienklewicZ的直接帮助下分别就基于 p收敛和 收敛方法的 最优离散化问题作了一些工作，并且首次应用于接触 问题，取得了较为满意的结果 C7F82