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

加权残量法在建立多变量有限元模型中的应用Ξ Vol. 21 No. 1 第 2 1卷 第 1 期 固 体 力 学 学 报 2000年 3月 March 2000 ACTA MECHANICA SOLIDA SINICA 加权残量法在建立多变量有限元 X 模型中的应用 史宝军 鹿晓阳 山东建筑工程学院计算力学研究所 ,济南 ,250014 摘 要 应用加权残量法基本原理 ,给出了建立多变量有限元模型的简明格式和新的列 式 ,讨论了该列式与几种典型多变量有限元模型的对应、 退化和等价关系.从而阐明了各类多变 量有限元模型之间的内在联系 ,为各种多变量有限元模型的相互促进与发展提供了有效途径. 并探讨了多场变量的选择原则与单元刚度阵秩的关系.文后给出了若干典型算例. 关键词 加权残量法 ,多变量有限元 ,物理2几何关系 ,对应性 1 引言 应用加权残量法求解问题近似解的方法 ,具有方法原理的统一性 ,应用的广泛性 ,简便 [1~4 ] 准确等优点 .加权残量法在有限元法中也获得了应用 ,如邱吉宝将有限元法和最小二乘 [3 ] 法相结合 ,提出了最小二乘有限元法 ;赵俊卿、 龙驭球用变分法和加权残量法相结合 ,构造 [4 ] 非协调元 等. 为了改善单变量有限元的计算精度和几何适应性 ,许多学者提出了各种类型的多变量 [5 ] [6 ,7 ] [8 ] [9 ,10 ] 有限元模型.如杂交应力元 ,拟协调元 ,广义杂交元 ,杂交/混合元 等.这些单元 是多变量有限元模型中的典型代表.多变量有限元在计算精度和适应性上表现出明显的优 越性.那么 ,各种多变量有限元模型之间是否存在着某种对应关系呢是否存在着统一的数 学基础呢这些问题的研究对各种多变量有限元的相互促进与发展具有重要意义. [11~13 ] 目前 ,已有一些关于单元之间对应性和等价性的研究 ,但多限于 Wilson 元及其修 [14 ] 正模型 为代表的非协调位移元与杂交元/广义杂交元之间的讨论.而关于前述几种典型 多变量有限元模型之间的对应性还没有人进行过系统的研究. 另外 ,在建立多变量有限元模型时 ,还必须考虑多场变量之间的匹配问题.场变量选择 的是否合理 ,直接关系到单元性能的好坏.其中 ,应力 应变 场插值矩阵的选择是一个很重 要的问题 ,它直接关系到单元刚度阵的秩及单元的力学性能. 本文首先应用加权残量法基本原理 ,导出了建立多变量有限元模型的简明格式和新的 列式.然后从有限元列式的角度 ,对各种多变量有限元模型之间的对应关系及退化关系和等 价关系进行了探讨.讨论了多场变量的合理选择与单元刚度阵秩的关系.为各种多变量有限 元模型的相互借鉴和相互促进与发展提供了有效途径. X 1998205205收到第 1稿 ,1999205214收到修改稿 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd All rights reserved.第 1期 史宝军等 : 加权残量法在建立多变量有限元模型中的应用 ? 4 1 ? 2 建立多变量有限元模型的简明格式 在有限元法应用于结构和固体力学的问题中 ,定义节点力和节点位移关系的单元刚度 阵表征单元的特性.因此 ,最关键的步骤是单元刚度阵的形成.下面应用加权残量法基本原 理 ,从物理2几何关系的弱形式 1 式出发 ,直接离散单元应力场 ,导出建立多变量有限元模 型的简明格式. T Φ σ ν S - Du d 0 1 ∫ ν σ 这里 S - Du 0 2 σ σ 为变形体的物理2几何关系.其中 S 和 Du 可分别看作由假定应力场 和位移场 u 而产生的 σ 应变场.如果此假定的应力场 和位移场 u 恰好是变形体的真实应力场和位移场 ,则物理2 几何关系 2 式应精确满足 ,这对以分片试函数为基础的有限元法来说是很困难的.故宜从 其弱形式出发离散单元应力场. σ σ β Φ 设单元应力场 和位移场 u 分别为 P , u Nq ,取权函数 P ,则由式 1 可得 Hβ Gq 3 T H P S Pdν 4 ∫ ν T ν G P DN d 5a ∫ ν 0 对易于构造协调位移场的 C 类问题 ,可直接应用 3 ~ 5a 式离散单元应力场.而对具有 1 C 类连续性要求的问题 ,如板壳元模型的建立等 ,构造协调位移场是很困难的 ,为此 ,借鉴 [6 ,7 ] 拟协调元 的做法 ,对式 3右端分部积分 T T T ν Gq - D P ud + M P ~ u d s 6 ∫ ∫ ν ν 设域内位