基于解析试函数的广义协调边形膜元.pdfVIP

第 19卷第 4期 工 程 力 学 VoI．19No．4 2002年 8 月 ENGINEERING M匝CHANICS Aug． 2002 文章编号： 1000-4750(2002)04—012—05 基于解析试函数的广义协调四边形膜元 傅 向荣，龙驭球 (清华大学土木系，北京 100084) 摘 要 ：本文构造三个广义协调四边形膜元 。根据弹性力学平面 问题 的控制方程和艾雷应力函数 ，求出问题基 本解析解，然后用其作为试函数来构造单元：ATF—Q4a、 一Q4b、 一Q40。数值算例表 明，其中两个单 元ATF—Q4a和 ATF—Q4o对网格畸变不敏感 ，显示出良好 的性能 。 关键词：有 限元；试函数 ；艾雷应力函数；网格畸变敏感性 中图分类号 ：TBl15 文献标识码：A l 引言 一 组广义协调 四边形膜元 。探讨 了解析试函数法在 平面膜元 中的应用 。 在有 限元离散法 中吸取解析法 的优点而加 以 利用，在文献 中早有论述 。例如： 1982年文献 [1】 2 弹性平面 问题 的基本解析解 采用含应力奇点的解析解来构造奇异单元，为裂纹 在弹性理论 中，位移、应变、应力须满足几何 问题提出了有效解法，文献 [2】又作 了进一步探讨 。 方程、本构关系和平衡方程 。在平面弹性 问题 中， 1996年，文献 3『】提 出了理性有限元，并利用弹性 由双调和 的艾雷应力函数推导出的应力场、应变 平面 问题基本解析解构造 了一组优质平面单元 。文 场、位移场可 以满足以上基本物理控制方程。表 1 4【】利用厚板理论的基本解析解构造 了优质厚薄通 给出了弹性平面 问题的 l2个低阶解析解 。 用单元 ATF—PQ4 。本文利用弹性平面 问题基本解 析解作为试函数 ，应用广义协调元理论L5J，构造 了 表 1 弹性平面 问题的基本解析解 Table1 Fundamentalanalyticalsolutionsofplnaeelasticityproblems 项数 l 2 3 4 5 6 7 8 9 10 ll 12 ￡x O O O O —LL l 一21．tx 2y 一2py 2x 一61．txy 6xy ￡v O O O O l —LL 2x 一21．ty 2y 一2gx 6xy 一61．txy v O O O