有限元法课程学习

（5）伽辽金法简单地说，将近似解的试探函数作为权函数。更简洁的形式：伽辽金法的一般表达式引入变分等效积分形式第1页/共50页第一页，共51页。8/9/20132静态线弹性有限元定解问题真实位移的变分，连续可导。在给定位移的边界上，虚应变高斯定律张量形式矩阵形式第2页/共50页第二页，共51页。等效积分形式：与原有微分方程和定解条件完全等价。加权余量法：对场函数进行近似，令加权余量等于零。伽辽金法：加权函数与场函数的试探函数（基函数、形函数）相同。小结：伽辽金法是有限元法中使用最为普遍的。第3页/共50页第三页，共51页。

基本概念偏微分方程和偏微分方程组：一个未知函数及其偏导数组成的方程叫偏微分方程，两个以上未知函数及其偏导数组成的方程组叫偏微分方程组。方程组中未知函数和方程个数相等，叫封闭的偏微分方程组（或完全的）。偏微分方程的阶和偏微分方程组的阶：方程中偏导数的最高阶次叫偏微分方程的阶;偏微分方程组的阶是方程组中各偏微分方程的阶数之和。第4页/共50页第四页，共51页。线性、非线性和拟线性偏微分方程:a)方程中所有出现未知函数或其偏导数的项都是未知函数的一次式的方程叫线性方程b)未知函数项或未知函数偏导数项不是一次式的方程叫非线性方程；c)非线性方程中所有未知函数的最高阶偏导数是一次式的方程叫拟线性方程。第5页/共50页第五页，共51页。齐次和非齐次偏微分方程偏微分方程分类：Auxx+2Buxy+Cuyy+Dux+Euy+f=0Δ（=B2–AC）>0：双曲型，波动方程.Δ（=B2–AC）=0：抛物型，热传导方程.Δ（=B2–AC）<0：椭圆型，

位势方程.第6页/共50页第六页，共51页。定解问题：偏微分方程+定解条件（边界条件+初始条件）a)初值问题n阶方程有n个初始条件，初始条件偏导数的最高阶次是n–1.b)边值问题i:第一类边界条件（Dirichlet）,ii:第二类边界条件（Neuman）,iii:第三类边界条件（Robin）.第7页/共50页第七页，共51页。

偏微分方程的解a)解：使偏微分方程两端恒等的有定义的函数叫偏微分方程的解.b)通解:对于n阶方程，未知函数有m个自变量，其通解由n个独立的满足一定可微要求的函数组成，且每个函数有m–1个自变量.

uxyz=0,u=f1(x,y)+f2(x,z)+f3(y,z)

utt–c2uxx=0,u=f1(x+ct)+f2(x–ct).

c)定解问题的解：满足边界条件和初始条件的通解.第8页/共50页第八页，共51页。

线性叠加原理解的存在、唯一性和稳定性解的性质a）椭圆形方程的极值只能在边界达到。解在内部没有弱间断，解在边界上间断，在内部也是充分光滑的，边界条件是封闭的.b）双曲型方程没有像椭圆形方程那样的极值原理，解在内部可以有弱间断，边界条件不是封闭的。第9页/共50页第九页，共51页。

泛函与变分（Functionalandvariation）1．泛函函数的函数a)两端固定的曲线长度：b)弹性杆的总势能：c)温度场泛函：式中f,u,T叫做泛函的容许函数：满足一定边界条件和连续性的所有函数有限元法的基本原理曲线长度总势能温度场泛函第10页/共50页第十页，共51页。

变分定义

a)容许函数的变分有限元法的基本原理第11页/共50页第十一页，共51页。第12页/共50页第十二页，共51页。i）

泛函的值由1个自变量的函数确定ii）泛函的值由有3个自变量的函数确定iii）泛函的值由有3个自变量的2个函数确定第13页/共50页第十三页，共51页。d）变分运算第14页/共50页第十四页，共51页。3．变分问题a)函数的极值问题（无约束和约束）b)变分问题：求泛函的极值函数c)泛函极值函数的必要条件d)欧拉方程第15页/共50页第十五页，共51页。有限元法的基本原理第16页/共50页第十六页，共51页。e)定解问题与变分问题此式即杆的平衡方程，它就是变分的欧拉方程。有限元法的基本原理i）

固定边界变分问题与基本边界条件两端约束的弹性杆问题：第17页/共50页第十七页，共51页。ii）非固定边界变分问题与自然边界条件边界条件：泛函：此式即杆的平衡方程第18页/共50页第十八页，共51页。一端约束（指定位移）的弹性杆解法1：Lagrange乘子法构造新泛函iii）含有约束条件的变分问题Lagrange乘子法与原定解问题完全等效，代价就是引入了新的求解变量第19页/共50页第十九页，共51页。解法2：罚函数法构造新泛函有限元法的基本原理当k无穷大时，则满足第一类边界条件。不引入新的求解变量。第20页/共50页第二十页，共51页。小结：1.定解问题（微分方程加定解条件）等价于相应泛函取极值。2.泛函取极值就是有限元方法的理论基础，将微分形式变成了积分形式。3.不是所有的定解问题都存在相应的泛函。4.不存在泛函的定解问题，可以直接用更广义的加权余量法。第21页/共50页第二十一页，共51页。泛函的变分弹性长杆的定解问题微分方程定解条件对应泛函第22页/共50页第二十二页，共51页。近似函数表示的微分方程的残差边界条件的残差其思想是使由近似函数表示的微分方程的残差和边界条件的残差的加权积分为零

2．加权余量法直接从微分方程出发的一种积分方法。假设未知函数采用近似函数表达：有限元法的基本原理第23页/共50页第二十三页，共51页。伽辽金加权余量法弱形式与变分结果一致。分部积分伽辽金加权余量法

3．加权余量法弱形式：泛函取极值变分等于零可以看出，第24页/共50页第二十四页，共51页。取加权函数的试探函数与近似解的试探函数相同

加辽金有限元法解题过程1）构造加权残差积分方程2）离散化3）单元分析4）总体分析5）建立和求解有限元方程组4．加辽金加权余量法第25页/共50页第二十五页，共51页。弱形式1）构造加权余量积分方程第26页/共50页第二十六页，共51页。b）单元的总体节点编号和局部节点编号，单元e=I，总体节点编号：1，2，局部节点编号：1，2；单元e=III，总体节点编号：3，4，局部节点编号：1，2。2）离散化第27页/共50页第二十七页，共51页。a)在自然坐标系中构造单元近似解：b）构造加权函数c）单元坐标：3）单元分析加权函数、近似解试探函数、坐标插值函数的类型一致第28页/共50页第二十八页，共51页。d）单元平衡方程第29页/共50页第二十九页，共51页。a）建立选择矩阵：4）总体分析第30页/共50页第三十页，共51页。b)组集单元刚度矩阵c)组集等效节点载荷d）解以节点为未知量的方程组第31页/共50页第三十一页，共51页。热传导问题的有限元方法第32页/共50页第三十二页，共51页。1）傅里叶定律热传导方程1.一维问题热流密度与温度梯度成正比。q:单位时间、单位面积流过的热量单位：W/(m·K)第33页/共50页第三十三页，共51页。2）平衡方程Q=cpmΔT比热容：cp1千克的物质的温度上升（或下降）1摄氏度所需的能量。单位：W·s/(kg·K)扩散率的单位：m^2/s第34页/共50页第三十四页，共51页。第35页/共50页第三十五页，共51页。本构关系第36页/共50页第三十六页，共51页。第37页/共50页第三十七页，共51页。3）三维各向同性导热材料第38页/共50页第三十八页，共51页。第39页/共50页第三十九页，共51页。第40页/共50页第四十页，共51页。第41页/共50页第四十一页，共51页。第42页/共50页第四十二页，共51页。第43页/共50页第四十三页，共51页。第44页/共50页第四十四页，共51页。第45页/共50页第四十五页，共51页。第46页/共50页第四十六页，共51页。第47页/共50页第四十七页，共51页。第48页/共50页第四十八页，共51页。第49页/共50页第四十九页，共51页。感谢您的观看！第50页/共50页第五十页，共51页。内容总结