有限元第四章一些数学和结论材料.pptVIP

§4-6 协调位移单元的收敛性和误差估计 u Sh uh H1 πu 0 图4-15 1. 有限元空间Sh 对于任何实际需要用有限元分析的问题首先要做到: (i) 将区域Ω划分为m个单元（在平面情况下可以是三角元或矩形元）。 (iii) 选择每个单元内的插值多项式。多项式满足协调性和可微性要求，且多项式的系数要能由单元所包括的结点参数唯一确定。 (ii) 配置好每个单元的结点，选择好每个结点参数。对于二阶问题总以结点函数值为结点参数，以下我们仅限于讨论二阶问题。 得到一组基函数φ1、φ2、…φn。它们满足协调性和可微性要求；每个基函数只在一个结点上为１，在其余结点上为零。 我们把形如u =∑uiφi　（其中φi为定义的基函数，ui 任意常数）的函数全体称为由基函数φ1、φ2、…φn 张成的有限元空间，记作 Sh。 Πu —以精确的结点参数值进行插值所得到的分片插值函数。 Sh 是一个有限维空间。 uh —有限元解，结点参数由δπP= 0 定出，一般仅是近似值。 ２. 有限元解的投影关系 (4-6-1) (4-6-2) (4-6-3) 的势能驻值条件可表示为：对任意 有 二阶问题的真实解 u 可以表述为：寻找一个 u ∈H1(Ω), 使得对任何　δu ∈H1(Ω)（记作 ?δu∈H1(Ω)） 都有 有限元解 uh 则可表述为：寻找一个uh∈Sh，使得对? δu∈Sh 都有 将（4-6-1）与（4-6-2）相减 （4-6-3）表明：真实解 u 与有限元解 uh 之差与Sh 的任一元素 在能量内积的意义下正交。（图4-15）。换句话说：uh 是 u 在 Sh 上的投影。 u Sh uh H1 πu 0 图4-15 即在能量模意义下，uh 是 Sh 的所有元素中与 u “最接近”的一个元素。 对于πu （由结点的精确值构造的试探函数）有 (4-6-4) 证明如下： (4-6-5) 所以有 由于 上式说明了这样的事实：πu（精确的结点参数值构成的插值函数）在能量模的意义下并不比uh（由近似的结点参数值构成的插值函数）在整个求解域内更好（图4-16）。 图4-16 x uh x1 u x2 x3 0 ３ . 收敛性证明 设真实解 u ∈HK(Ω)，各单元的分片插值多项式完全到 k－１次 （k－1次称为有限元空间Sh 的次数），有限元解为uh。 对于二阶问题，能量模最多同时包括函数及函数的一阶导数。对每个单元 ei 都有 (4-6-6) (4-6-7) 误差的能量模估计: (4-6-8) (4-6-9) 当h→0 时‖u-uh‖→0 函数的精度高于导数的精度 §4-7 有限元（协调位移单元）解的特点 设 u 为真实解 ( u ∈H1(Ω)) ，则总势能为 (4-7-1) (4-7-2) 由对 ? δu ∈H1(Ω)有 当取δu=u 时也应该成立，故有 类似对于有限元解 uh∈Sh 也有 (4-7-5) (4-7-4) (4-7-3) (4-7-5) (4-7-4) (4-7-3) 当外载荷仅是一个集中力（假如这种载荷是允许的话），那么力作用点的位移的有限元解一般比真实解小，这就是通常所说的，“有限元解的位移偏小、刚度偏大）的确切涵义。 有限元解的总势能高于真实解的总势能 有限元解的变形能小于真实解的变形能 外力在有限元解uh上做的功小于外力在真实解u上做的功 §4-8 结束语 ! 本章的结论对任何一种协调位移单元都适用，从理论角度讲，我们已经“相当彻底”地解决了(插值逼近，收敛性，有限元解的特点)问题。 ! 借助于一些数学概念对协调位移单元进行了进一步讨论(插值逼近，收敛性，有限元解的特点)。(仅限于讨论二阶问题)。 ! 有限元解的精度不但依赖于插值多项式的次数，而且依赖于广义解的可微性。 ! 就平面问题而言，我们仅构造出三角元和矩形元两种单元。换句话说，我们构造单元的实际能力相对说来还显得薄弱。 第四章　一些数学概念和结论 关于有限元方法的一些数学概念和结论； 有限元解的收敛性以及单元精度问题； 本章的主要对象是函数： 真实解是一个函数； 基函数是一组函数； 试探函数是某一类函数； 有限元解是某类函数中使πP 取最小值的那一个函数。 本章中： “元素”指的就是函数； “空间” 就是具有某种性质的函数的集合，即函数空间。 §4-1 线性空间（向量空间） 1. 线性空间的定义 满足下列条件的空间E为线性空间