第二章--有限单元法的数学基础

第二章有限单元法的数学基础,有限单元法方程的基本思想,按其推导方法分:(1)直接法根据单元的物理意义，建立有关场变量表示的单元性质方程。(2)加权余量法直接从控制方程中得到有限单元方程，是一种近似解法。(3)变分法直接从求解泛函的极值问题入手，把泛函的极植问题规划成线性代数方程组，然后求其近似解的一种计算方法。,有限单元法的数学基础,有限元最初是由结构工程师们基于物理概念提出来的，但这种被称为直接刚度法的方法仅适用于力学领域，很难将其推广到其它领域。,我们今天要面临的问题，不仅限于力学问题，要涉及到其它许多相关领域的问题，如，当我们运用数值方法模拟降雨型滑坡的成因和机理时，我们不只是考虑应力场，还需考虑地下渗流场、地表径流等诸多因素的共同作用。,变分法是建立有限元法的通用方法，同时变分法还是其它现代数值方法，如边界元法、数值流形法、有限体积法等的基础。,有限单元法的数学基础,自然界中的许多规律都是以微分方程的形式提出的，微分方程是关于场变量及其导数的一个方程，微分方程和相应的初边值条件一起构成了所谓的微分方程的定解问题。,三大类微分方程：椭圆型方程用于描述与时间无关的稳态问题；抛物型方程的典型应用是热传导问题；双曲型方程的典型对象则是振动问题,由于问题的复杂性，能够获得解析解的问题是非常有限的。为了求解一般情况下的定解问题，就需要借助于某种数值方法。,有限单元法的数学基础,数值方法，粗略的理解就是一种近似方法，但这种近似方法不是任意的，而是能够通过数值方法本身所提供的某种手段，通过加大运算量而使得近似解越来越逼近解析解。,有限元法就是求解微分方程的一种数值方法，而且它特别适用求解椭圆型边值问题。,在有限元法问世之前所广泛采用的一种数值方法是有限差分法，有限差分法是一切微分方程数值方法的基础，至今仍在某些领域如流体力学问题中起着主导作用。有限差分法的基本思想就是用差商代替微商，而使得原来的定解问题的求解变为线性方程组的求解。所有的微分方程数值解法都是化微分方程定解问题为线性方程组的求解问题。,有限单元法的数学基础,有限元法的先行课程(按依赖程度)：线性代数、数理方程、数值分析+一门连续统力学，这门连续统力学对不同学科有不同的要求，对于土木工程学科，可以是弹性力学或结构力学。,Einstein求和约定,引入目的：书写和推导公式的简便；阅读其它采用张量记法的有关文献。,首先将三维空间的坐标x，y，z分别记为,Einstein求和约定：如果表达式中的乘积项中含有两个相同指标，则表示这个指标遍历1，2，n，n=2或3，是空间的维数，而后将这n个同类型的表达式相加而得到的总和。如果相同的指标有好几对，则表示每一对指标做上述求和。,Einstein求和约定,相同的指标是求和指标，称为哑指标，而其余的指标称为自由指标。,求和约定中求和指标遍历的是空间维数n，若在一个求和表达式中的求和指标不是遍历空间维数n，则我们不应省略求和号,求和约定的好处示例：三维变形体的平衡方程,变分法的基本引理,子空间,如果存在,在区间a,b上有连续n阶导数的全体函数空间为,假设,则在区间a,b上,变分法的基本引理,证明：采用反证法。,设存在一点,根据函数的连续性存在,构造连续函数,但,微分方程的等效积分形式,工程和物理学中的许多问题，都可以以未知场函数应满足的微分方程和边界条件的形式表达。未知场函数u,微分方程的等效积分形式,根据前面的引理微分方程组可以采用完全等效的积分形式进行表达,微分方程的等效积分形式,分部积分：假设,Gauss公式：区域积分向边界积分转换,分部积分公式,微分方程的等效积分形式,等效积分的“弱”形式,运用分部积分法：,为什么叫“弱”形式？,原来的方程中对于V并无连续性要求，现在的V连续性要求提高了，同时对U的连续性要求降低了。,微分方程的等效积分形式,示例：二维稳态热传导方程,热源密度,热传导系数,温度,边界给定温度,边界给定热流,微分方程的等效积分形式,其等效积分形式为：,通过选择函数自动满足,强制边界条件,微分方程的等效积分形式,其等效积分弱形式的建立过程：（1）分部积分,（2）按积分区域合并,微分方程的等效积分形式,（3）简化,边界上,函数的任意性,令,则等效积分弱形式为：,其边界条件不见了？,自然边界条件,如果前面已经执行强制边界条件此处可取零去掉积分项,加权余量法思想,将未知场函数u采用近似函数表达，一般取：,例如当u是位移时，可以,N：形函数,加权余量法思想,取,则,W：权函数,余量的加权积分为零即可得到一组求解方程，因此称为加权余量法。,加权余量法思想,权函数的选择有多种方法：,（1）配点法,（2）子域法,（3）最小二乘法,强迫余量在某些点处为0,强迫余量在某些子域积分为0,强迫,取最小值,加权余量法思想,权函数的选择有多种方法：,（4）力矩法,（5）伽辽金法（Galerkin法）,强迫余量的各次矩为0，又称为积分法。,采用这种方法得到的求解方程系数矩阵是对称的，所以应用最为广泛。,加权余量法示例,示例：二阶常微分方程,边界条件：,数值求解：首先取近似解,该近似解已经满足了边界条件。强制边界条件如何执行？,该近似解在域内产生余量R，余量的加权积分应为0：,加权余量法示例,近似解的项数取得越多，则计算精度越高，本例中我们只讨论一项和两项近似解。,一项近似解：,代入微分方程，余量R为：,二项近似解：,代入微分方程，余量R为：,加权余量法示例,A.配点法求解,配点法即强迫余量在某些点处为0.,二项近似解：配点取x1=1/3，x2=2/3,一项近似解：配点取x=1/2,解得：,解得：,加权余量法示例,B.子域法求解,子域法即强迫余量在某些子域内积分为0.,二项近似解：子域取0,1/2,1/2,1,一项近似解：子域取全域，,解得：,解得：,加权余量法示例,C.最小二乘法求解,最小二乘法实质即余量的二次方积分后取最小值.,对I的表达式求导：,余量的二次方积分：,要让积分值取最小值，则,权函数,加权余量法示例,二项近似解：,一项近似解：,解得：,解得：,加权余量法示例,D.力矩法求解,力矩法即强迫余量在区域内的各次矩为0.,二项近似解：W1=1，W2=x,一项近似解：W1=1,解得：,解得：,加权余量法示例,E.伽辽金法求解,伽辽金法实质即取近似函数作为权函数.,一项近似解：,则权函数为：,解得：,加权余量法示例,二项近似解：,则权函数为：,解得：,加权余量法示例,该问题的解析解为,什么是泛函,如果对某一类函数中的每一个函数v(x)，按照一定的变换规则J，都有一个确定的实数，记为J(v)，与之对应，那么就称该变换是一个泛函。泛函就是由函数空间到一维欧氏空间R的一个变换，泛函的“自变量”或者“原象”取自某一类函数，而变换的结果是一个实数。,例如,由于其满足了,所以可以称为线性泛函。,下面两个属于非线性泛函。,什么是泛函的变分,如果泛函J的增量J可表示为,则称泛函J可微，并称,假设J可微，则对于任一实数a,有,其中,关于v是线性的,为一阶变分，简称变分，记为,等式两端除a并令,泛函极值的充分条件和必要条件,如果函数u(x)M使得对应的J(u)取极值，则,仿照函数的Taylor展开,的极值点为,若且0(或0),则u(x)是J(v)的一个极小或极大函数,泛函取极值的充分条件,若但随u有正有负，则u(x)是J(v)