FE-Ch01.3变分原理与里兹法.pptVIP

泛函： 变分： 相互独立， 所以， 或 由： 得到矩阵形式 其中 共有 3?n 个方程， 若 为完备的函数系列 则， 时， 收敛于精确解， 若 n 为有限项，则 为近似解。 上述方法为Ritz法 2）将 代入 Ritz（里兹）法——基于变分原理的近似解法 1.求解步骤： 1）假设近似解： 为待定参数， 满足强制边界条件。 泛函 的极值问题（求函数u），转化为求 多元（ ）函数的极值问题。 3）求解线性代数方程组 u的近似解 2.解的收敛性 1）连续性要求 满足 阶连续性 2）完备性要求 取自完备的函数序列 3.特点 1) 近似解对全域而言 2) 试探函数要求满足一定的边界条件，近似解的 精度与试探函数的选择有密切关系。 3) 待定系数 任意，不表示特定的物理意义。 4) 如果我们对问题了解比较清楚，能找到合适的 试函数，可以说事半功倍，但缺乏一般性。 4.讨论： 1）经典意义上的泛函变分理论只适应于线性 自伴随微分方程。 2）收敛性有严格的理论基础（泛函分析）。 3）事先满足强制边界条件，则解有明确的上 下界性质。如不事先满足，需要进行处理 （约束变分原理）。 但是未知函数 往往还需要服从一些附加条件， 1.修正（约束）变分原理 建立了自然变分原理后，问题的解为泛函?取驻值。 约束条件 我们把这些变分原理称之为“具有附加条件的变分原理”。 可以将附加条件引入泛函，重新构造一个“修正泛函”，把问题转化为求修正泛函的驻值问题。常用方法： Lagrange 乘子法，罚函数法。 原泛函?的约束变分问题，转化为修正泛函?*的 无约束变分，代价是修正泛函增加了附加未知函 数 。 2. Lagrange 乘子法(?乘子法) 修正泛函 ?* ： 单变量泛函 双变量修正泛函. 近似解： ， 线性 修正泛函的变分： 其中： 可得： 即 的系数阵为0。 所以 方程中不含 项， 对线性问题，得线性方程组； 因为： 整理得到 讨论（放松约束条件的代价）： 1）很明显方程的阶数增加了。 2）方程的系数矩阵主元（对角元素）出现零元素，对求解方程增加了困难。（不能用一般的消元法） 3）一般的物理问题中得到的自然变分问题是一极值问题。而对修正的泛函，由于附加项的积分性质不清，一般为驻值问题。（不再有极值性质） 4）利用?乘子法，弹性力学各种变分原理的转换。 3.罚函数法 修正泛函 其中?称为罚数 正定的， 对?为极小值问题，? 取正数； ? 值越大，约束条件满足的越好。(近似性越好) 这种方法好处很明显，不增加任何未知函数。 ( ? 是事先给定的 ) 例：极值问题（函数极值问题） 约束条件 ，所以 ： 解方程，得： 上述方程可写为矩阵形式 分析方程： 来自原来的泛函， 来自约束条件。 * §1.3 变分原理与里兹法 “最速落径问题”---质量为m的小环从A处自由滑下， 试选择一条曲线使所需时间最短。（不计摩擦） A B X Y 设路径为 y=y（x） 所需时间 a y 称T为y（x）的泛函， y（x）为自变函数。 即以函数作自变量以积 分形式定义的函数为泛函。 一. 变分的一些基本概念 X A Y 变分运算在形式上与微分运算相同。 y=y（x） x+dx dy x 称 为y（x）的变分，它是一个无穷小的任意函数。 微分与变分运算次序可以交换。 积分与变分运算次序也可以交换。 二. 函数的定义和泛函的定义 1. 函数的定义： 若对于自变量x域中的每一个值，y有一值与之对应，或数y对应于数x的关系成立。则称变量y是变量x的函数，即： y= y( x)。 2. 泛函的定义： 若对于某一类函数{y(x)}中的每一函数 y(x)，? 有一值与之对应，或数 ? 对应于函数 y(x) 的关系成立。则称变量 ? 是函数 y(x) 的泛函，即： ?= ?(y(x))。 3. 微分和变分 微分 : x 的增量 △x 是指某两值之差 △x=x-x1 .如果 x 的微分用 dx 表示，则 dx 也是增量的一种，即当这种增量很小很小时， dx= △x。 变分: y(x)的增量在它很小时称为变分，用?y(x)或?y表示， ?y(x)是指y(x)和与它相接近的y1(x)之差，即 ?y(x)=y(x)-y1(x);这里： ?y(x)也是x的函数，只是 ?y(x)在指定的x域中都是微量。（假定y(x)在接近y1(x)的一类函数中是任意改变的 ）。 4. 函数的微分和泛函的变分 函数的微分:函数的增量 △y=y(x+ △x)-