构造几何不敏感四边形膜元的广义协调方法1.PDF

第　29 　卷　第　6 　期 力 　学 　学 　报 　 Vol. 29 , No. 6 1997 　年　11 　月 ACTA MECHANICA SINICA Sept . , 1997 构造几何不敏感四边形膜元的广义协调方法1) 龙驭球 　　须 　寅 (清华大学土木工程系 , 北京 100084) 摘要 　许多有限元模型在规则网格下具有很高的精度 , 但当网格畸变程度增大时 , 其计算精 度也随之迅速下降. 如何构造出对网格畸变不敏感的单元 , 长期以来是人们十分关注的课 题. 本文应用作者早期提出的广义协调方法 , 构造出具有平面内旋转自由度的任意四边形膜 元. 该单元不仅列式简洁 , 而且具有对网格畸变极不敏感的优异性能 , 为构造对网格畸变不 敏感的优质单元提供了一个通用方法. 关键词 　广义协调元 , 四边形膜元 , 几何不敏感 , 有限元 引 　言 协调位移型有限元是最简单、可靠的有限元形式 , 因而得到了广泛的应用. 然而 , 对网格 形状畸变敏感是这种单元模式比较突出的问题. 它导致单元计算精度随网格畸变程度的增加而 迅速下降 , 严重影响了这类单元的适应性. 这个问题一直受到国内外学者的重视. 国际著名有 限元学者 Pian [1 ] , Wilson[2 ]近年分别提出的杂交应力元和非协调位移元方案 , 虽然取得了较好 的效果 , 但其单元对网格畸变仍有不同程度的敏感. 从单元实用角度上讲 , 单变量位移元占有 很大优势. 因此 , 在协调位移模式中引入不协调的高次位移项 , 既能改善单元性能又可保持列 式的简洁. 但这种不协调的高次位移项的确定具有一定程度的任意性 , 而且不易保证单元的收 敛. Wilson[3 ]最早提出的非协调元就是一个典型的代表. 为了克服非协调位移元不能保证在任何情况下都能通过分片检验这一致命缺点 , 本文应用 广义协调方法[4 ] , 引入含 8 个内部自由度的低阶附加位移场 , 构造了广义协调四边形膜元 GQ12M8. 该单元在收敛性、计算精度均有出色的表现 , 尤其具有对网格畸变极不敏感的突出 特点. 1 　单元位移函数 本文所讨论的任意四边形膜元如图 1 所示, 其每个角点 i 含有两个线自由度 u , v 和一个 i i [ 5 , 6] 平面内附加刚体旋转自由度 θ . 单元结点自由度可用下列向量形式表示 i { q}e = [{ q } T 　{ q } T 　{ q } T 　{ q } T ]