有限单元法及数学基础.docVIP

有限单元法的数学基础 1、引言 有限元方法归根结底是一种数值计算方法，它有严格的数学证明作为其近似的客观性和合理性的保证。力学问题最终归结为一组微分方程的边值问题或者初值问题抑或是混合问题。比如弹性静力学最终归结为L-N方程的微分提法。在很难或者根本不可能得到所得方程的理论解的情况下，究竟用什么样的方法才能得到方程的近似解（这种近似解已经能够满足实际工程的需要），在这种情况下，二十世纪五六十年代由结构力学家进而由数学家提出和证明了这种思想方法的合理性。有限元方法产生于力学计算，但是，它本质上并不是力学的专利。世间万物的变化过程很多都可以通过微分方程特别是偏微分方程来描述，也就是说，微分方程是很多现象和过程的数学结构，而大多数的微分方程是不能得到理论解的，这时候就可以使用有限元方法来求其近似解，因为有限元方法是求解微分方程（组）的数值计算方法。它适用于力学的微分方程，也同样适用于其它领域的相应的微分方程的数值求解。 2、有限元方法数学根源 对于一个给定的微分方程定解问题，为了求其近似解，我们可以使用Ritz方法和Galerkin方法。下面分别阐述这两种方法，然后讨论有限元方法和他们的关系。 Ritz法 Ritz法源于最小势能原理，设H是可分的Hilbert空间，在H中取有限维空间Sn，它是由N个线性无关向量张成，即： 用代替H，在上求泛函J(w)的极值,即求，使得 = 实际上寻求只需通过解一个线性方程组 D--------双线性形式 F--------线性泛函 因此，是一个以为未知数（自变量）的二次多项式，如果二次项的系数矩阵 是正定的，那么在N+1维空间是一个开口向上的椭球抛物面，它有且只有一个极（最）小值点，所谓在上求的极值，就是确定，使得： ＝ 极值条件： |＝0 （） 得: () 即： 适合方程组： KC=F 。。。。。。。 总结：我们可以发现，Ritz法主要体现在两点：第一，用有限维试解函数空间的基函数生成真解的一个近似函数，这个试验函数的系数是基本未知量；第二，在求解过程中进行了积分弱化，也就是说，由原来方程的高阶可微要求，通过一次积分和分步积分，得到了弱解形式；第三，使用势能驻值原理。 Galerkin法 Galerkin法源于虚位移原理，虚位移原理的本质只是一个数学恒等式。Galerkin方法也是用SN代替H,求UN∈SN，使得对任意WN∈SN 有 且满足 仍然只需要解一个代数方程组，事实上，设： 则： 0 = ＝ 由于是任意的，必须有： 因此得到与Ritz方法相同的代数方程组，此时U就是在取的Galerkin意义下的广义解。 总结：我们可以发现，Galerkin法主要体现在两点：第一，用有限维试解函数空间的基函数生成真解的一个近似函数，这个试验函数的系数是基本未知量；第二，由高阶可微变换成弱解积分形式；第三，使用虚位移原理。 有限元法 有限元方法和Ritz法、Galerkin法有什么样的关系呢？ 通过上面的分析可以看出，除去由于方法不同可能带来的数值运算差别之外，仅从原理上来说，Ritz法和Galerkin法位移的区别就在于前者使用了势能驻值原理，而后者使用了虚位移原理。而有限元法属于Ritz—Galerkin法的范畴，但又是改进，关键在于有限元法选取分块多项式函数类作为子空间SN。也就是说有限元方法可以从Ritz法导出，也可以从Galerkin法导出（即既可以基于势能驻值原理导出也可以基于虚位移原理导出），但是唯一的区别是它的试解函数空间选取不同。它采用分片插值函数空间，这正好体现了有限元剖分，而上面两种方法虽然也是近似，但是没有分片的概念。而这正是有限元的核心。 具体说，就是先把Ω区间剖分为一系列单元，然后在每个单元内做多项式插值，而使他们在单元的公共点、边、面上满足一定的连续性条件，以保证这种函数组成的有限维空间是H的子空间。 所以，有限元的基础是变分原理和分片多项式插值。 下面，以二维Poisson方程的第一边值问题为例来说明有限元方法的这种思想。 Poisson方程的边值问题： 满足边值条件： 第一边值条件 即 求广义解，就是求 ，使得： 其中： 对Ω进行三角剖分，记三角形单元为，Ω的内部节点为 ，边界节点为，所有单元的最大边长是h，对每一个单元e=△做线性插值，使它在三个顶点各取已知值,当节点满足上时，给定值为0，这样在Ω上就构成了一个多面体函数，当遍历空间所有向量时，函数就构成了一个N维函数空间。显然，是一个线性空间，且是（在意义下完备化）。 作为的NP维线性子空间，它在NP个基函数，使得可以表示成： 有限元方法就是求，使得：