第二讲有限元法的理论基础演示

（优选）第二讲有限元法的理论基础现在是1页\一共有142页\编辑于星期一预备知识第一章有限单元法的理论基础

1.1微分方程的等效积分形式

1.2加权余量法

1.3变分原理主要内容现在是2页\一共有142页\编辑于星期一弹性力学的基本假设预备知识现在是3页\一共有142页\编辑于星期一一、连续性假设弹性理论同其他宏观物理学一样，不考虑实际工程材料细观粒子结构。1.物体抽象成连续密实的空间几何体，位移、应变、应力、能量等物理量作为空间点位置的函数定义在这个几何体上。2.物体在整个变形过程中始终保持连续，即：定义在该连续介质上的物理性质和物理量除了在某些孤立的点、线、面上可能奇异或间断外，在变形过程中始终保持为空间点位的连续函数。预备知识现在是4页\一共有142页\编辑于星期一二、弹性假设弹性体的变形与载荷在整个加载和卸载过程中存在一一对应的单值函数关系，且载荷卸去后变形完全消失。应力小于弹性极限时应力应变关系是线性的。服从虎克定律。小变形情况下，应变和位移导数间的关系是线性的。预备知识现在是5页\一共有142页\编辑于星期一三、均匀性假设物体在各点处的弹性性质都相同。四、自然状态假设假设物体不受外力作用和温度的影响，物体便没有应力和变形，即不考虑由于制造工艺引起的残余应力和装配应力。预备知识现在是6页\一共有142页\编辑于星期一弹性力学问题的矩阵表示预备知识现在是7页\一共有142页\编辑于星期一一、基本物理量位移：应变：应力：

预备知识现在是8页\一共有142页\编辑于星期一一、场方程几何方程：预备知识现在是9页\一共有142页\编辑于星期一预备知识现在是10页\一共有142页\编辑于星期一物理方程：这里假设材料是各向同性的。预备知识现在是11页\一共有142页\编辑于星期一注：

表示工程切应变，它们与张量切应变的关系为：预备知识现在是12页\一共有142页\编辑于星期一在平面问题中的弹性矩阵：平面应力问题：平面应变问题：预备知识现在是13页\一共有142页\编辑于星期一平衡方程：预备知识现在是14页\一共有142页\编辑于星期一边界条件：力边界：位移边界：预备知识现在是15页\一共有142页\编辑于星期一本章重点和应掌握的内容本章重点和应掌握的内容微分方程的等效积分形式及其“弱”形式的实质和构造方法，任意函数和场函数应满足的条件。不同形式加权余量法中权函数的形式和近似解的求解步骤，以及Galerkin法的特点。线性自伴随微分方程的变分原理的构造方法和泛函的性质，以及自然边界条件和强制边界条件的区别。第1章有限元法的理论基础

现在是16页\一共有142页\编辑于星期一经典Ritz方法的求解步骤、收敛条件及其局限性两种形式虚功原理（虚位移原理和虚应力原理）的实质和构造方法。从虚功原理导出最小位能原理和最小余能原理的途径，各自的性质以及场函数事先应满足的条件第1章有限元法的理论基础

现在是17页\一共有142页\编辑于星期一本章含盖三节内容：1.1微分方程的等效积分形式1.2加权余量法1.3变分原理第1章有限元法的理论基础

现在是18页\一共有142页\编辑于星期一1.1微分方程的等效积分形式第1章有限元法的理论基础

现在是19页\一共有142页\编辑于星期一1.1微分方程的等效积分形式微分方程：微分方程是联系自变量x，未知函数u(x)和它的某些阶导数的关系式：现在是20页\一共有142页\编辑于星期一1.1微分方程的等效积分形式求解微分方程的方法有：解析法；半解析法；数值法；现在是21页\一共有142页\编辑于星期一1.1微分方程的等效积分形式数值法主要包括：有限差分法——将微分方程化为差分形式，求近似解；加权余量法——将~转化为加权积分形式，求近似解；有限元法——将~转化为能量取驻值问题，并采用分片插值；边界元法——在边界上进行离散；无网格法——近似函数建立在离散点上，不需网格。现在是22页\一共有142页\编辑于星期一微分方程的等效积分形式一、连续介质问题微分方程的一般表达式且满足边界条件：

表示对独立变量(时间，空间)的微分算子。1.1微分方程的等效积分形式现在是23页\一共有142页\编辑于星期一微分方程的等效积分形式1.1微分方程的等效积分形式现在是24页\一共有142页\编辑于星期一微分方程的等效积分形式1.1微分方程的等效积分形式现在是25页\一共有142页\编辑于星期一微分方程的等效积分形式1.1微分方程的等效积分形式现在是26页\一共有142页\编辑于星期一微分方程的等效积分形式1.1微分方程的等效积分形式现在是27页\一共有142页\编辑于星期一微分方程的等效积分形式1.1微分方程的等效积分形式现在是28页\一共有142页\编辑于星期一微分方程的等效积分形式1.1微分方程的等效积分形式现在是29页\一共有142页\编辑于星期一微分方程的等效积分形式1.1微分方程的等效积分形式现在是30页\一共有142页\编辑于星期一微分方程的等效积分形式1.1微分方程的等效积分形式现在是31页\一共有142页\编辑于星期一微分方程的等效积分形式1.1微分方程的等效积分形式现在是32页\一共有142页\编辑于星期一微分方程的等效积分形式1.1微分方程的等效积分形式现在是33页\一共有142页\编辑于星期一微分方程的等效积分形式1.1微分方程的等效积分形式现在是34页\一共有142页\编辑于星期一微分方程的等效积分形式例：图1：u为一个连续函数，满足C0连续图2：有一个一阶不连续点，但一阶导可积。图3：二阶导数在Δ区域内趋于无穷，使积分不能进行。现在是35页\一共有142页\编辑于星期一微分方程的等效积分的弱形式1.1微分方程的等效积分形式一、构造“～弱形式”目的降低对未知函数的连续性的要求假设：微分方程中，微分算子的最高阶导数为2m；现在是36页\一共有142页\编辑于星期一微分方程的等效积分的弱形式1.1微分方程的等效积分形式现在是37页\一共有142页\编辑于星期一微分方程的等效积分的弱形式1.1微分方程的等效积分形式现在是38页\一共有142页\编辑于星期一微分方程的等效积分的弱形式1.1微分方程的等效积分形式3)代价是提高对任意函数和的连续性要求。4)在物理上更符合实际问题对连续性的要求。5)若和取特定函数，则为加权余量法

的不同格式。现在是39页\一共有142页\编辑于星期一微分方程的等效积分的弱形式1.1微分方程的等效积分形式例：简支梁的弯曲问题微分方程和边界条件现在是40页\一共有142页\编辑于星期一微分方程的等效积分的弱形式1.1微分方程的等效积分形式微分方程的等效积分形式如下：对该等效积分形式要求在域内，w为三阶导数连续，很难实现。现在是41页\一共有142页\编辑于星期一微分方程的等效积分的弱形式1.1微分方程的等效积分形式等效积分弱形式：对等效积分弱形式要求在域内，w一阶导数连续即可。现在是42页\一共有142页\编辑于星期一微分方程的等效积分的弱形式1.1微分方程的等效积分形式方程的分类：1）稳态问题（平衡边值问题）场函数解只与位置坐标有关现在是43页\一共有142页\编辑于星期一微分方程的等效积分的弱形式1.1微分方程的等效积分形式方程的分类：1）瞬态问题（传播问题，初边值问题）场函数为空间与时间的函数Г、Ω可以理解为时-空域，t为开域（0，∞）t=0可以认为是初值条件现在是44页\一共有142页\编辑于星期一微分方程的等效积分的弱形式1.1微分方程的等效积分形式方程的分类：1）特征值问题若要有非零解某些参数取特定值取决于问题的物理、几何特性现在是45页\一共有142页\编辑于星期一1.2加权余量法第1章有限元法的理论基础

现在是46页\一共有142页\编辑于星期一1.2加权余量法加权余量法的基本思想加权余量法是：

基于等效积分形式或等效积分弱形式的近似方法。现在是47页\一共有142页\编辑于星期一1.2加权余量法设：定解问题现在是48页\一共有142页\编辑于星期一1.2加权余量法1.构造近似解现在是49页\一共有142页\编辑于星期一1.2加权余量法那么，当n有限时，方程存在偏差（余量）即：在域Ω内在边界Г上现在是50页\一共有142页\编辑于星期一1.2加权余量法等效积分形式：现在是51页\一共有142页\编辑于星期一1.2加权余量法2.以加权意义上为零，形成求解方程组（等效积分的解析式）即：或：为权函数，（预先设定）线性无关。作用：强迫余量在某种平均意义上等于零现在是52页\一共有142页\编辑于星期一1.2加权余量法现在是53页\一共有142页\编辑于星期一1.2加权余量法3.加权余量法的关键(两种函数的选择)1）与等效积分形式不同：一个是精确解，而加权余量法得到的为是近似解。a.近似表达式为有限项。b.对某些特定的权函数（非任意）2）试函数：如能满足一定的域内条件或边界条件，使问题简化，且有一定的精确度。3）权函数：不同的权函数，涉及不同的计算格式。例如：现在是54页\一共有142页\编辑于星期一1.2加权余量法采用使余量的加权积分为零的等效积分的“弱”形式，来求得微分方程近似解的方法称为加权余量法。它是求微分方程近似解的一种有效方法。现在是55页\一共有142页\编辑于星期一1.2加权余量法加权余量法常用的几种常用方案为了讨论方便，不失一般性，认为已满足边界条件，因此仅剩域内积分项；为线性微分算子，可用表示。现在是56页\一共有142页\编辑于星期一1.2加权余量法1.配点法取:则有:注：现在是57页\一共有142页\编辑于星期一1.2加权余量法1.配点法现在是58页\一共有142页\编辑于星期一1.2加权余量法1.配点法现在是59页\一共有142页\编辑于星期一1.2加权余量法1.配点法现在是60页\一共有142页\编辑于星期一1.2加权余量法1.配点法现在是61页\一共有142页\编辑于星期一1.2加权余量法1.配点法这种方法相当于简单地强迫若干个在域内的点上余量等于零。说明：Kij

非对称，不用求积分。现在是62页\一共有142页\编辑于星期一1.2加权余量法2.最小二乘法最小二乘法是加权余量法的一种。标准最小二乘法是：要使域Ω内每一点的残数（或误差）的平方和最小，或平方的积分最小。现在是63页\一共有142页\编辑于星期一1.2加权余量法2.最小二乘法现在是64页\一共有142页\编辑于星期一1.2加权余量法2.最小二乘法现在是65页\一共有142页\编辑于星期一1.2加权余量法2.最小二乘法现在是66页\一共有142页\编辑于星期一1.2加权余量法2.最小二乘法现在是67页\一共有142页\编辑于星期一1.2加权余量法2.最小二乘法可见：矩阵对称，但需要数值积分现在是68页\一共有142页\编辑于星期一1.2加权余量法3.伽辽金（Galerkin

）法非对称，系数矩阵含积分运算。若自伴随问题利用格林公式可以构造有限元格式现在是69页\一共有142页\编辑于星期一1.2加权余量法3.伽辽金（Galerkin

）法说明：<1>如果要形成有限元格式，则希望得到对称系数矩阵，同时希望积分中的微分阶数降低。<2>Galerkin加权余量法（见后）现在是70页\一共有142页\编辑于星期一1.2加权余量法3.伽辽金（Galerkin

）法现在是71页\一共有142页\编辑于星期一1.2加权余量法3.伽辽金（Galerkin

）法如果L为二阶微分算子，则C、D均为一阶。如果L为四阶微分算子，则C、D均为二阶。如果L为自伴随算子，第一项将得到对称系数矩阵。现在是72页\一共有142页\编辑于星期一1.2加权余量法3.伽辽金（Galerkin

）法例：二维稳态热传导方程（Galerkin格式）现在是73页\一共有142页\编辑于星期一1.2加权余量法3.伽辽金（Galerkin

）法现在是74页\一共有142页\编辑于星期一1.2加权余量法3.伽辽金（Galerkin

）法现在是75页\一共有142页\编辑于星期一1.2加权余量法3.伽辽金（Galerkin

）法现在是76页\一共有142页\编辑于星期一1.2加权余量法利用格林公式分部积分现在是77页\一共有142页\编辑于星期一1.2加权余量法不考虑温度边界条件，上式整理得：其中：现在是78页\一共有142页\编辑于星期一1.2加权余量法说明：（1）由Galerkin法得到与变分法相一致的方程形式，与有限元格式类似。（2）如离散后采用上法，即可得到有限元格式。（3）如果一个问题存在变分泛函，则采用加权余量法Galerlin格式与变分方法可得相同结果的方程。现在是79页\一共有142页\编辑于星期一变分原理自然变分原理修正泛函的变分原理有限元法的理论基础

现在是80页\一共有142页\编辑于星期一线性、自伴随微分算子如果微分方程具有线性、自伴随的性质，则：不仅可以建立它的等效积分形式，并可利用加权余量法求其近似解；

还可建立与之相等效的变分原理，基于它的另一种近似求解方法——Ritz法自然变分原理有限元法的理论基础-变分原理

现在是81页\一共有142页\编辑于星期一线性、自伴随微分方程的定义：微分方程：为微分算子若具有性质：则称为线性微分算子。有限元法的理论基础-变分原理

自然变分原理现在是82页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

自然变分原理现在是83页\一共有142页\编辑于星期一泛函的构造

设有微分方程：

有限元法的理论基础-变分原理

自然变分原理现在是84页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

自然变分原理现在是85页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

自然变分原理现在是86页\一共有142页\编辑于星期一自然变分原理有限元法的理论基础-变分原理

现在是87页\一共有142页\编辑于星期一变分原理是针对以下积分形式定义的标量泛函而言，有限元法的理论基础-变分原理

自然变分原理现在是88页\一共有142页\编辑于星期一原问题微分方程和边界条件的等效积分Galerkin提法等效于泛函取驻值。反之泛函取驻值则等效于微分方程和边界条件。这里泛函可以通过等效积分的Galerkin提法得到。这种变分原理称为自然变分原理。例如，弹性力学中的最小位能原理、粘性流体中最小能力耗散原理，称为自然变分原理。有限元法的理论基础-变分原理

自然变分原理现在是89页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

自然变分原理现在是90页\一共有142页\编辑于星期一最小位能原理：真实位移使体系总位能取极小值，即：有限元法的理论基础-变分原理

自然变分原理现在是91页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

自然变分原理现在是92页\一共有142页\编辑于星期一自然变分原理有限元法的理论基础-变分原理

现在是93页\一共有142页\编辑于星期一Ritz（里兹）法——基于变分原理的近似解法求解步骤假设近似解：为待定参数，满足强制边界条件。将代入泛函的极值问题（求函数u），转化为求多元（）函数的极值问题。有限元法的理论基础-变分原理

现在是94页\一共有142页\编辑于星期一求解线性方程组有限元法的理论基础-变分原理

现在是95页\一共有142页\编辑于星期一解的收敛性1)

连续性要求满足Cm-1

阶连续性2)完备性要求取自完备的函数序列有限元法的理论基础-变分原理

现在是96页\一共有142页\编辑于星期一

特点1)

近似解对全域而言2)试探函数要求满足一定的边界条件，近似解

的精度与试探函数的选择有密切关系。3)待定系数不表示特定的物理意义。4)如果我们对问题了解比较清楚，能找到合适的试函数，可以说事半功倍，但缺乏一般性。有限元法的理论基础-变分原理

现在是97页\一共有142页\编辑于星期一提示经典意义上的泛函变分理论只适应于线性自伴随微分方程。2)收敛性有严格的理论基础（泛函分析）。3)事先满足强制边界条件，则解有明确的上下界性质。如不事先满足，需要进行处理(约束变分原理)。有限元法的理论基础-变分原理

现在是98页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

关于强制边界条件与自然边界条件若微分算子是线性自伴随的，Galerkin法的等效积分形式问题泛函近似场函数应满足强制边界条件现在是99页\一共有142页\编辑于星期一假如微分算子是2m阶0至m-1阶导的边界条件称为强制边界条件m至2m-1阶导的边界条件称为自然边界条件未知场函数无需事先满足自然边界条件有限元法的理论基础-变分原理

现在是100页\一共有142页\编辑于星期一关于泛函取极值根据Galerkin格式或变分原理，微分算子线性自伴随：假设微分算子L的最高阶导数是2m偶数阶，则：有限元法的理论基础-变分原理

现在是101页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

现在是102页\一共有142页\编辑于星期一关于解的下限性：有限元法的理论基础-变分原理

现在是103页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

现在是104页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

现在是105页\一共有142页\编辑于星期一最小余能原理：真实解使得系统的总余能最小。考虑平衡方程：有限元法的理论基础-变分原理

现在是106页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

现在是107页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

现在是108页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

现在是109页\一共有142页\编辑于星期一最小势能原理解的下限性：由能量守恒定理知：变形过程中的功等于弹性体变形后的应变能。有限元法的理论基础-变分原理

现在是110页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

现在是111页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

现在是112页\一共有142页\编辑于星期一同样的分析得到：由最小余能原理得到的近似应力场，总体偏大。有限元法的理论基础-变分原理

现在是113页\一共有142页\编辑于星期一修正泛函变分原理

建立了自然变分原理后,问题的解为泛函Π取驻值。有限元法的理论基础-变分原理

修正泛函（约束）变分原理现在是114页\一共有142页\编辑于星期一但是未知函数往往还需要服从一些附加条件，

约束条件把这些变分原理称之为：

“具有附加条件的变分原理”。修正泛函（约束）变分原理有限元法的理论基础-变分原理

现在是115页\一共有142页\编辑于星期一解决的办法

可以将附加条件引入泛函，重新构造一个“修正泛函”，把原问题转化为求修正泛函的驻值问题。常用方法：Lagrange乘子法，罚函数法。修正泛函（约束）变分原理有限元法的理论基础-变分原理

现在是116页\一共有142页\编辑于星期一2.

Lagrange乘子法（λ乘子法）修正泛函（约束）变分原理有限元法的理论基础-变分原理

现在是117页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

修正泛函（约束）变分原理2.

Lagrange乘子法（λ乘子法）现在是118页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

修正泛函（约束）变分原理2.

Lagrange乘子法（λ乘子法）现在是119页\一共有142页\编辑于星期一修正泛函（约束）变分原理2.

Lagrange乘子法（λ乘子法）有限元法的理论基础-变分原理

现在是120页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

修正泛函（约束）变分原理2.

Lagrange乘子法（λ乘子法）现在是121页\一共有142页\编辑于星期一讨论（放松约束条件的代价）：1)很明显方程的阶数增加了。2)方程的系数矩阵主元（对角元素）出现零元素，对求

解方程增加了困难。（不能用一般的消元法）3)一般的物理问题中得到的自然变分问题是一极值问题。

而对修正的泛函，由于附加项的积分性质不清，一般

为驻值问题。（不再有极值性质）4)利用乘子法，可以得到弹性力学各种λ变分原理的转换。修正泛函（约束）变分原理有限元法的理论基础-变分原理

现在是122页\一共有142页\编辑于星期一3.罚函数对Π为极小值问题，α取正数；α值越大，约束条件满足的越好。(近似性越好)这种方法好处很明显，不增加任何未知函数。(α是事先给定的)有限元法的理论基础-变分原理

修正泛函（约束）变分原理现在是123页\一共有142页\编辑于星期一例：极值问题（函数极值问题）有限元法的理论基础-变分原理

现在是124页\一共有142页\编辑于星期一有限元法的理论基础-变分原理

现在是125页\一共有142页\编辑于星期一讨论：有限元法的理论基础-变分原理

修正泛函（约束）变分原理现在是126页\一共有142页\编辑于星期一修正泛函（约束）变分原理有限元法的理论基础-变分原理

讨论：现在是127页\一共有142页\编辑于星期一关键概念：等效积分形式等效积分“弱”形式加权余量法Galerkin方法线性自伴随算子泛函和变分原理强制边界条件自然边界条件泛函的驻值和极值Ritz方法

虚位移原理虚应力原理最小位能原理最小余能原理有限元法的理论基础

现在是128页\一共有142页\编辑于星期一有限元法的理论基础-课后作业

等效积分形式和等效积分“弱”形式的区别何在？为何后者在数值分析中得到更多的应用？不同形式的加权余量法之间的区别，你能提出其它形式的加权余量法吗？加权余量法的Galerkin方法特点，自然边界条件和强制边界条件的区别里兹法的特点及优缺点，与Galerkin方法的异同虚功原理有哪两种不同形式，和弹性力学