偏微分方程数值解讲稿.ppt

* （1）弦平衡的平衡原理与极小位能原理 2 两点边值问题的变分原理 考察一根长为 l 的弦，两端固定在点 A(0,0)和B(l,0)。当没有外力作用时，它的位置沿水平方向与X轴重合。设有强度为f(x)的外荷载垂直向下作用在弦上，于是弦发生形变。假定荷载很小，因而发生的形变也很小。用 u(x) 表示在荷载f(x)的作用下弦的平衡位置。 ? ? * 求弦的平衡位置归结为求解两点边值问题： 设弦处于某一位置u=u(x)，可得到其总位能为 极小位能原理： 其中T是弦的张力。 平衡原理 弦的平衡位置 (记为 )将在满足边值条件 u(0)=0，u(l)=0 的一切可能位置中，使位能取极小值。 弦的平衡位置 是下列变分问题的解 * 在数学上，要将某个微分方程的定解问题转化为一个变分问题求解，必须针对已给的定解问题构造一个相应的泛函，并证明定解问题的解与泛函极值问题的解等价。 有限元方法正是利用这种等价性（边值问题与变分问题的等价性），先将微分方程定解问题转化为变分问题（或变分方程）的求解问题，然后再设法近似求解变分问题（或变分方程）。 * （2）两点边值问题的变分原理 ① 构造泛函 考察二阶常微分方程边值问题： 引入泛函算子 则 * ② 变分问题 与前述二阶常微分方程边值问题相应的变分问题是 其中 求 ，使 * ③ 变分原理（变分问题与边值问题的等价性） 设 ， 是边值问题 的解，则 使 J(u) 达到极小值； 反之，若 使 J(u) 达到极小值，则 是边值问题的解。 其中 是强制边界条件， 是自然边界条件，区别这两类边界条件在用有限元方法求解边值问题时很重要。 * （3）虚功原理 对两点边值问题： 其中 虚功原理 ，且满足变分方程： 设 ，以v乘方程两端，沿[a,b]积分，并利用 ，得变分方程 对任意 在力学里， 表示虚功 设 ，则 是边值问题解的充要条件是： * 对于复杂的边界条件，边值问题的求解一般是困难的。若将微分方程化为相应的变分问题或变分方程，则只需处理强加边界条件，无需处理自然边界条件（自然边界条件已包含于变分问题中泛函的构造或已包含于给出的变分方程之中）。这一特点对研究微分方程离散化方法及其数值解带来了极大的方便。 * 3 二阶椭圆边值问题的变分原理 （1）极小位能原理 模型方程 其中G是平面有界区域。 ① 构造泛函 引入泛函算子 则 * ② 变分问题 与前述二阶椭圆边值问题相应的变分问题是 求 ，使 其中 * ③ 变分原理（变分问题与边值问题的等价性） 对第一边值问题，无论齐次或非齐次边界条件，泛函是一样的，只是边界条件要作为强加边值条件加在所取的函数类上。 设 ， 是二阶椭圆边值问题的解，则 使 J(u) 达到极小值； 反之，若 使 J(u) 达到极小值，则 是二阶椭圆边值问题的解。 其中 对第二、三类边值问题，无论齐次或非齐次边界条件，二次泛函形式相对于第一边值问题有所改变，但函数类的选取与边界条件无关。 * （2）虚功原理 问题 其中 设 ，以v乘方程两端后在G上积分，并利用Green公式，得变分方程 * 虚功原理 在力学里， 表示虚功 设 是边值问题的解，则对任意 ， 满足变分方程。 反之，若 ，且对任意 满足变分方程，则 为边值问题的解。 与极小位能原理类似，第一类边界条件为强加边界条件，第二、三类边界条件为自然边界条件。 虚功原理比极小位能原理应用更广。 * 目的：求解相应的变分问题或相应的变分方程。 Ritz方法是近似求解变分问题(即二次泛函极小值)的算法。Galerkin方法是近似求解变分方程的算法，这两种算法统称为R