1、毕业设计文献综述 信息与计算科学 线性有限元法的稳定性和误差分析 有限元 方法的基本思想 是用较简单的问题代替复杂问题后再求解 .它将求解域看成是由许多称为有限元的小的互连子域组成 , 对每一单元假定一个合适的 、 较简单的近似解 , 然后推导求解这个域总的满足条件 (如结构的平衡条件） , 从而得到问题的解 .这个解不是准确解 , 而是近似解 , 因为实际问题被较简单的问题所代替 .由于大多数实际问题难以得到准确解 , 而有限元不仅计算精度高 , 而且能适应各种复杂形状 , 因而成为行之有效的工程分析手段 .和 每一项新技术的推出 的背景一样 , 有限元方法的产生也 是由于时代的迫切需要 ,
2、 而新技术的出现后也需要经历历史的重重考验 .在上个世纪 40 年代 , 由于航空事业的快速发展 , 对飞机内部结构设计提出了越来越高的要求 , 即重量轻、强度高、刚度好 , 人们不得不进行精确的设计和计算 .正是在这一背景下 , 有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术领域 , 成为一种丰富多彩、应用广泛并且实用高效的 数值分析 方法 1,2. 关于有限元方法 早期的一些成功的实验求解方法与专题论文 , 完全或部分的内容对有限元技术的产生做出的贡献 , 首先在应用数学界第一篇有限元论文是 1943年 Courant R发表的 , 文中描述了他使用三角形区域的多项式函数来求
3、解扭转问题的近似解 , 由于当时计算机尚未出现 , 这篇论文并没有引起应有的注意 . 1956 年 , M.J.Turner (波音公司工程师 ), R.Clough, H.C.Martin 以 及 L.J.Topp 等四位共同在航空科技期刊上发表一篇采用有限元技术计算飞机机翼的强 度 的论文 , 文中 把这种解法称为刚性法 (Stiffness), 一般认为这是工程学界上有限元法的开端 . 1960 年 , RayClough 教授在美国土木工程学会 (ASCE)会议上 , 发表一篇名为 The Finite Element in Plane Stress Analysis 的论文 , 将应
4、用范围扩展到飞机以外之土木工程上 , 同时有限元法 (Finite Element Method)的名称也第一次被正式提出 .由此之后 , 有限元法的理论迅速地发展起来 , 并广泛地应用于各种力学问题和非线性问题 , 成为分析大型、复杂工程结构的 强有力手段 .并且随着计算机的迅速发展 , 有限元法中人工是难以完成的大量计算工作能够由计算机来实现并快速地完成 .因此 , 可以说计算机的发展很大程度上促进了有限元法的建立和发展 . 有限元方法在国内的产生和发展情况大致如下 , 我国的力学工作者为有限元方法的初期发展做出了许多贡献 , 其中比较著名的有 : 陈伯屏（结构矩阵方法） , 钱令希（余能
5、原理） , 钱伟长（广义变分原理） , 胡海昌（广义变分原理） , 冯康（有限单元法理论） 3. 遗憾的是由于当时环境所致 , 我国有限元方法的研究工作受到阻碍 , 有限元理论的发展也逐 渐与国外拉开了距离 . 20 世纪 60 年代初期 , 我国的老一辈计算科学家较早地将计算机应用于土木、建筑和机械工程领域 .当时黄玉珊教授就提出了 “ 小展弦比机翼薄壁结构的直接设计法 ” 和“ 力法应力设计法 ” ; 而在 70 年代初期 , 钱令希教授提出了 “ 结构力学中的最优化设计理论与方法的近代发展 ” . 这些理论和方法都为国内的有限元技术指明了方向 . 1964 年初崔俊芝院士研制出国内第一个
6、平面问题通用有限元程序 , 解决了刘家峡大坝的复杂应力分析问题 . 20 世纪 60 年代到 70 年代 , 国内的有限元方法及有限元软件诞生之后 , 曾计算过 数十个大型工程 , 应用于水利、电力、机械、航空、建筑等多个领域 4,5,6. 关于有限元法的一些基本知识 , 文献 2,7,8,9对本文的研究起到至关重要的作用 ,本文首先 引入两点边值问题 2 dd( ) ( ) , ,dd( ) 0 , ( ) 0 ,uL u p qu f x a x bxxu a u b 其中 1m in0( ) ( ) ,( ) , ( ) 0 ,( ) , ( ) 0 .f x C Ip C I p x
7、pq C I q x 使用有限元方法求解的稳定性 . 这一部分 , 可以参考文献 3, 将上述问题变成等价的变分问题 相应的矩阵表达形式 bAU , 我们可以讨论 对有限元方程 bAU 的 系数矩阵、右端向量发生扰动对解向量的影响 , 特别地 , 当 A 对称时 , 取m inm ax)( Acond , 我们有结论 7:存在正常数 C , 使得 2m in)( ChAcond , h 为剖分单元的步长 . 我们还可以对上述两点边值问题 , 讨论其有限元解法的收敛性 , 参考文献 7, 可得该问题的变分形式即为 求线性有限元解 )()( 1 IHVxu EhEh ,使 hEhhhh Vvvfv
8、ua ),(),( , 其中 dd( , ) ( ) d ,dd( , ) d ,b hhh h h habhhauva u v p q u v xxxf v f v x 我们可以得到 正交投影性质 hEhhh Vvvuua ,0),( , 以及最佳逼近性质 hh Vv hh vuuu 0in f1 , 其中 nk ehh kxuxuxuxu 1 2 ,121 )()()()( , 2 221,( ) ( ) ( ( ) ( ) ) ( ( ) ( ) ) dk kh h he eu x u x u x u x u x u x x , 特别地 , 参考文献 10, 我们还可以得到 11 Ih
9、uuCuu , hEI Vxu )( , 为 )(xu 的分段线性插值多项式 , 即满足插值条件 nkxuxu kkI )1(0),()( , 于是 , 我们可以得到线性有限元解函数 )(xuh 与真解函数 )(xu 之间的误差在空间 1H 的范数意义下有如下误差估计式 011 uChuuCuu Ih , 其中 2 20 ( ( ) dbau u x x , 利 用 Nitsche 技巧 , 可得误差估计式 020 uChuu h . 参考文献 1 R. A. Adams. Sobolev spaces. Academic Press, New York, 1975. 2 李荣华 . 偏微分方
10、程数值解法 . 北京 : 高等教育出版社 .2005. 3 冯康 . 基于变分原理的差分格式 . 应用数学与计算数学 , 1965, 2(4):237261. 4 夏道行 , 吴卓人 , 严绍宗等 . 实变函数论与泛函分析 M. 北京 : 高等教育出版社 , 1985. 5 Wu Haijun and Li Ronghua. Error estimate for finite volume element methods for general second elliptic problems. NM for PDE, 2002, 693708. 6 舒适 . 偏微分方程典型离散化方法的基本理
11、论与算法分析 . 内部讲义 , 2007, 568. 7 陈传淼 , 黄云清 . 有限元高精度理论 . 湖南科技出版社 , 1995. 8 C. Bernardi. and R. Verfurth. Adaptive finite element for elliptic equations with nonsmooth coefficients. Numerische Mathematik, 2000.85: 579608. 9 A. Bowyer. Computing Dirichlet Tessellations. Computer Journal, 1981, 24(2):162166. 10 李开泰 , 黄庆怀 . 有限元方法及其应用 . 科学出版社 , 2006.
