有限元法的理论基础 .PPT

求解线性方程组 有限元法的理论基础-变分原理 解的收敛性 1) 连续性要求 满足Cm-1 阶连续性 2) 完备性要求 取自完备的函数序列 有限元法的理论基础-变分原理 特点 1) 近似解对全域而言 2) 试探函数要求满足一定的边界条件，近似解 的精度与试探函数的选择有密切关系。 3) 待定系数不表示特定的物理意义。 4) 如果我们对问题了解比较清楚，能找到合适 的试函数，可以说事半功倍，但缺乏一般性。 有限元法的理论基础-变分原理 提示 经典意义上的泛函变分理论只适应于线性自伴随 微分方程。 2) 收敛性有严格的理论基础（泛函分析）。 3) 事先满足强制边界条件，则解有明确的上下界性质。如不事先满足，需要进行处理(约束变分原理)。 有限元法的理论基础-变分原理 有限元法的理论基础-变分原理 关于强制边界条件与自然边界条件 若微分算子是线性自伴随的， Galerkin法的等效积分形式 问题泛函 近似场函数 应满足强制边界条件 假如微分算子是2m阶 0至m-1阶导的边界条件称为强制边界条件 m至2m-1阶导的边界条件称为自然边界条件 未知场函数无需事先满足自然边界条件 有限元法的理论基础-变分原理 关于泛函取极值 根据Galerkin格式或变分原理， 微分算子线性自伴随： 假设微分算子L的最高阶导数是2m偶数阶， 则： 有限元法的理论基础-变分原理 有限元法的理论基础-变分原理 1.2 加权余量法 2.最小二乘法 最小二乘法是加权余量法的一种。 标准最小二乘法是： 要使域 Ω内每一点的残数（或误差）的平方和最小，或平方的积分最小。 1.2 加权余量法 2.最小二乘法 1.2 加权余量法 2.最小二乘法 1.2 加权余量法 2.最小二乘法 1.2 加权余量法 2.最小二乘法 1.2 加权余量法 2.最小二乘法 可见：矩阵对称，但需要数值积分 1.2 加权余量法 3.伽辽金（Galerkin ）法 非对称，系数矩阵含积分运算。 若自伴随问题 利用格林公式 可以构造有限元格式 1.2 加权余量法 3.伽辽金（Galerkin ）法 说明： 1 如果要形成有限元格式，则希望得到对称 系数矩阵，同时希望积分中的微分阶数降低。 2 Galerkin加权余量法（见后） 1.2 加权余量法 3.伽辽金（Galerkin ）法 1.2 加权余量法 3.伽辽金（Galerkin ）法 如果L为二阶微分算子，则C、D均为一阶。 如果L为四阶微分算子，则C、D均为二阶。 如果L为自伴随算子，第一项将得到对称系数矩阵。 1.2 加权余量法 3.伽辽金（Galerkin ）法 例：二维稳态热传导方程（ Galerkin 格式） 1.2 加权余量法 3.伽辽金（Galerkin ）法 1.2 加权余量法 3.伽辽金（Galerkin ）法 1.2 加权余量法 3.伽辽金（Galerkin ）法 1.2 加权余量法 利用格林公式分部积分 1.2 加权余量法 不考虑温度边界条件，上式整理得： 其中： 1.2 加权余量法 说明： （1）由 Galerkin 法得到与变分法相一致的方程形式，与有限元格式类似。 （2）如离散后采用上法，即可得到有限元格式。 （3）如果一个问题存在变分泛函，则采用加权 余量法 Galerlin 格式与变分方法可得相同结果的方程。 变分原理 自然变分原理 修正泛函的变分原理 有限元法的理论基础 线性、自伴随微分算子 如果微分方程具有线性、自伴随的性质，则： 不仅可以建立它的等效积分形式，并可利用加权余量法求其近似解； 还可建立与之相等效的变分原理，基于它的另一种近似求解方法——Ritz法 自然变分原理 有限元法的理论基础-变分原理 线性、自伴随微分方程的定义： 微分方程： 为微分算子 若 具有性质： 则称 为线性微分算子。 有限元法的理论基础-变分原理 自然变分原理 有限元法的理论基础-变分原理 自然变分原理 泛函的构造 设有微分方程： 有限元法的理论基础-变分原理 自然变分原理 有限元法的理论基础-变分原理 自然变分原理 有限元法的理论基础-变分原理 自然变分原理 自然变分原理 有限元法的理论基础-变分原理 变分原理是针对以下积分形式定义的标量泛函而言， 有限元法的理论基础-变分原理 自然变分原理 原问题微分方程和边界条件的等效积分Galerkin 提法等效于泛函取驻值。反之泛函取驻值则等效于微分方程和边界条件。 这里泛函可以通过等效积分的Galerkin提法得到。 这种变分原理称为自然变分原理。 例如，弹性力学中的最小位能原理、粘性流体中最小能力耗散原理，称为自