加权残量法在建立多变量有限元模型中的应用.PDF

Vol. 21 No. 1 　　 　　　　　　 　 　　　　　　　　　 　 第21卷 第1期 固 体 力 学 学 报 2000年 3月 　 ACTA MECHANICA SOLIDA SINICA March 2000 加权残量法在建立多变量有限元 模型中的应用 史宝军 　　鹿晓阳 ( 山东建筑工程学院计算力学研究所,济南 ,250014) 摘　要 　应用加权残量法基本原理 ,给出了建立多变量有限元模型的简明格式和新的列 式 ,讨论了该列式与几种典型多变量有限元模型的对应、退化和等价关系. 从而阐明了各类多变 量有限元模型之间的内在联系 ,为各种多变量有限元模型的相互促进与发展提供了有效途径. 并探讨了多场变量的选择原则与单元刚度阵秩的关系. 文后给出了若干典型算例. 关键词 　加权残量法 ,多变量有限元 ,物理几何关系 ,对应性 1 　引言 应用加权残量法求解问题近似解的方法 ,具有方法原理的统一性 ,应用的广泛性 ,简便 准确等优点[1～4 ] . 加权残量法在有限元法中也获得了应用 ,如邱吉宝将有限元法和最小二乘 法相结合 ,提出了最小二乘有限元法[3 ] ;赵俊卿、龙驭球用变分法和加权残量法相结合 ,构造 [4 ] 非协调元 等. 为了改善单变量有限元的计算精度和几何适应性 ,许多学者提出了各种类型的多变量 有限元模型. 如杂交应力元[5 ] ,拟协调元[6 ,7 ] ,广义杂交元[8 ] ,杂交/ 混合元[9 ,10 ] 等. 这些单元 是多变量有限元模型中的典型代表. 多变量有限元在计算精度和适应性上表现出明显的优 越性. 那么 ,各种多变量有限元模型之间是否存在着某种对应关系呢 ? 是否存在着统一的数 学基础呢 ? 这些问题的研究对各种多变量有限元的相互促进与发展具有重要意义. 目前, 已有一些关于单元之间对应性和等价性的研究[11～13 ] ,但多限于 Wilson 元及其修 正模型[14 ]为代表的非协调位移元与杂交元/ 广义杂交元之间的讨论. 而关于前述几种典型 多变量有限元模型之间的对应性还没有人进行过系统的研究. 另外 ,在建立多变量有限元模型时 ,还必须考虑多场变量之间的匹配问题. 场变量选择 ( ) 的是否合理 ,直接关系到单元性能的好坏. 其中 ,应力 应变 场插值矩阵的选择是一个很重 要的问题 ,它直接关系到单元刚度阵的秩及单元的力学性能. 本文首先应用加权残量法基本原理 ,导出了建立多变量有限元模型的简明格式和新的 列式. 然后从有限元列式的角度 ,对各种多变量有限元模型之间的对应关系及退化关系和等 价关系进行了探讨. 讨论了多场变量的合理选择与单元刚度阵秩的关系. 为各种多变量有限 元模型的相互借鉴和相互促进与发展提供了有效途径. 收到第 1 稿 收到修改稿. © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 第 1 期 　　　　　　史宝军等 : 　加权残量法在建立多变量有限元模型中的应用 　　　　　　　　　　　 ·41 · 2 　建立多变量有限元模型的简明格式 在有限元法应用于结构和固体力学的问题中 ,定义节点力和节点位移关系的单元刚度 阵表征单元的特性. 因此 ,最关键的步骤是单元刚度阵的形成. 下面应用加权残量法基本原