有限元法基础讲义
- 格式:ppt
- 大小:2.90 MB
- 文档页数:351
有限元方法讲义第1讲抛物问题有限元方法1、椭圆问题有限元方法考虑椭圆问题边值问题:(1)问题(1)的变分形式:求使满足(2)的性质,广义解的正则性结果。
区域的剖分,矩形剖分,三角剖分,剖分规则,正则剖分条件,拟一致剖分条件。
剖分区域上分片次多项式构成的有限元空间。
的逼近性质,逆性质:这里,为的插值逼近。
问题(2)的有限元近似:求使满足(3)(3)的解唯一存在,且满足。
(3)的解所满足的矩阵方程(离散方程组)形式:(4)刚度矩阵的由单元刚度矩阵组装而成。
模误差分析:由(2)-(3)可得(5)由(5)可首先得到则得到(6)-模误差分析设满足用与此方程做内积,由(5)式和插值逼近性质得到再利用模误差估计结果,得到(7)最优阶误差估计和超收敛估计概念。
当与时间相关时(为抛物问题准备),由(5)式得(8)利用(7),类似分析可得(9)2、抛物问题半离散有限元方法考虑抛物型方程初边值问题:(10)(10)的变分形式:求使满足(11)(11)的半离散有限元近似:求使满足(12)令,代入(12),依次取可导出常微分方程组:(13)其中为质量矩阵,K为刚度矩阵。
求解常微分方程组(13),得到代回的表达式,即得半离散有限元解。
定理1.问题(12)的解唯一存在且满足稳定性估计:(14)证明:在(12)中取得到整理为(注意是正定的)对此式积分,证毕。
误差分析。
引进解的椭圆投影逼近:满足(15)根据椭圆问题的有限元结果可知(16)分解误差:的估计由(16)式给出,只须估计。
由(11),(12)和(15)知,满足取,类似稳定性论证可得可取为的投影,插值逼近等。
由(17)式,三角不等式和(16),得到(18)3、抛物问题全离散有限元近似剖分时间区间:。
引进差分算子:规定,当为连续函数时,,则有由此得到(19)(20)定义问题(11)的全离散向后Euler有限元近似:求,使满足(21)将代入(21)可导出全离散方程组(22)其中。