有限元法的误区和拟协调元-NSFC.PDF

中国科学: 物理学 力学 天文学 2014 年 第44 卷 第7 期: 1– 13 《中国科学》杂志社 SCIENTIA SINICA Physica, Mechanica Astronomica SCIENCE CHINA PRESS 论 文 有限元法的误区和拟协调元 ① ②* ② 唐立民 , 胡平 , 夏阳 ① 大连理工大学工程力学系, 大连 116024 ② 大连理工大学工业装备结构分析国家重点实验室, 运载工程与力学学部汽车工程学院, 大连 116024 *联系人, E-mail: pinghu@ 收稿日期: 2013-05-13; 接受日期: 2014-04-01 国家自然科学基金重点项目(批准号:和国家自然科学基金(批准号:资助项目 摘要 现行有限元法的理论和作法存在误区, 对有限元函数空间和弱形式导数认识不足, 因而阻碍了 有限元法的进一步发展. 本文给出了拟协调元的理论和作法作为对照. 利用函数序列表示有限元函数空 间, 提出了单元独立性原理, 证明协调条件不是单元函数构造时必须考虑的因素. 强调多项式基函数在 构造单元函数中的作用, 讨论了单元函数的构造, 指出单元函数应随着单元细化收敛于相应的泰勒级数. 证明了使用弱形式平衡方程时, 必须同时使用弱形式的协调方程. 提出网线弱导数, 拓广了有限元中弱 导数的内涵. 关键词 网线弱导数, 有限元函数空间, 单元独立性原理, 拟协调元, 有限元 PACS: 46.15.-x, 45.10.-b, 02.70.-c, 46.70.-p doi: 数值计算包括两个部分: 第一, 用什么算子, 例 方法的研究滞留了一百多年. 直到有限元法在二十 如微分方程、变分和差分等; 第二, 用什么插值函数, 世纪五、六十年代得到大量应用后, 对一般弹性力学 例如解析的幂级数、三角函数或是离散的差分格式等. 问题的求解才得到解决. 算子和插值函数的选择决定了数值计算方法的性能. 有限元法巧妙的采用了所谓“分片多项式”, 将 例如, 离散的差分法计算格式简单, 但是计算量较大, 连续的整体分割成许多个单元, 把解析分析和离散 而且难以处理不规则形状的边界条件. 如果在整个 的数值分析结合起来. 由于单元变小, 所以单元域的 区域用一个解析的幂级数多项式作为基函数, 那就 近似函数只需几项多项式, 简单方便, 可以适用于复 要用很多项进行逼近, 增加了计算量, 而且也存在难 杂的边界条件. 有限元的“分片多项式”方法是了不 以处理复杂边界条件的问题. 在工程实际中, 需要计 起的一大创举. 自二十世纪五十年代初, 伴随着电子 算量小、精度高, 可以处理复杂边界条件的数值方法. 计算机运用的出现, 有限元方法取得了辉煌成功, 对 尽管弹性力学的微分方程在十九世纪二十年代就建 弹性力学和机械、船舶、汽车、航空航天等学科及行 立了, 但是只能停留在理想的理论阶段, 有效的求解 业产生了巨大的影响, 至今仍然作为基础的现代计 引用格式: 唐立民, 胡平, 夏阳. 有限元的误区和拟协调元. 中国科学: 物理学 力学 天文学, 2014, 44: 1– 13 Tang L M, Hu P, Xia Y. The misconception in finite element metho