[工学]有限元原理加权余量法和变分法.ppt

第三讲 1.偏微分方程求解－－有限元法的原理（加权余量法和变分法） 解析法 应用范围有限，适用于理论求解，但有强烈的物理含义（常系数微分方程） 某些复杂问题，很考虑根本找不到解析解 2. 数值法 工程实际中应用广泛，复杂场域问题，但物理含义不很清楚。任何问题总可以找到数值解（数学方法） 2.数值求解方法 2/4 1. 基本思想： 以偏微分方程的近似解来代替其真解，只要近似解与真解足够接近，就可以近似解作为问题的解，并满足足够的精度。 2. 基本方法： 假设一个近似解，该解为一组（形式上）简单函数 的线性组合来表示，线性组合的系数就是一组待定系数 然后建立一种考虑了微分方程和边界条件的关于真解 和近似解间误差的目标函数 F 用适当的算法使得该目标函数最小化――最小化的过程就确定了待定系数，从而也就得到了问题的近似解。 尝试函数，基函数，形函数 2.数值求解方法 2/4 目标函数最小化的目的：一方面，使得近似解最大程度接近真解； 另一方面，求得构成近似解的待定系数。 数学上，构成目标函数的方法很多，不同的构成方法就形成了不同的数值解法，电磁场中就常见的是：加权余量法和变分法。 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 电磁场问题总可以用位函数的偏微分方程和相应的边界条件表述 两个偏微分方程形式相同，故以电位方程的求解过程为例。磁位矢量的方程可以分解到各个分量上变为标量方程。 在求解场域内，偏微分方程的真解为 ，近似解为 它由一组简单函数 的线性组合表达，表达中有待定系数 即： 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 加权余量法 简单函数，一般选用简单形式的函数，一旦选定就是已知的了 待定系数是真正的求解目标 问题的自由度 近似解 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 加权余量法就是一种定义近似解与真解之间误差（即余数），并设法使其最小的方法。 加权余量法误差（即余数）的定义： 注意：一般余数并不表示近似解与真解间的代数差（场域内），加权余量法的采用拉普拉斯算子作用后的差别（即余数），来代表近似解整体接近偏微分方程真解的程度。 问题的自由度 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 当余数小于要求的精度时，就可以认为近似解就是偏微分方程的解。 要减少余数，我们可以通过寻求适当的待定系数来实现。 为有效表达减小余数的效果，还选取适当的加权函数，以使余数和该加权函数的积分为0。－－“加权余量法”的来由。 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 加权余数的定义： 加权函数的选取方法很多：如点重合、子域重合、最小二乘法、迦辽金法。 效果较好的、运用较多的是迦辽金法： 即：迦辽金法选取尝试函数本身为加权函数 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 由此构建加权量法的目标函数： 上述过程中，已经将偏微分方程转化为j个代数方程组，便于计算机求解。 关于函数的函数，称为：泛函数，或泛函 3. 加权余量法－－例1 例1.两极电容板内部电场分布问题： 根据问题特点将3维问题简化为2维， 进一步简化为1维。 该问题是静态电场问题， 偏微分方程和边界条件： 加权余量法求解： 1.选取尝试函数、构造近似解： 2.结合问题，写出余数表达式： 3. 加权余量法－－例1 理论上任意选取，操作中越简单越好 2.结合问题，写出余数表达式： 3. 加权余量法－－例1 3. 加权余数表达式： 3. 加权余量法－－例1 3. 加权余数表达式： 3. 加权余量法－－例1 4. 求解上述两个代数方程组，得到待定系数，从而确定近似解 3. 加权余量法－－例1 加权余量法求解流程： 1.选取尝试函数、构造近似解 2.结合问题，写出余数表达式 3. 写出加权余数表达式 4. 令各加权余数表达式为0，得到代数方程组，解之得到待定系数，从而确定近似解 该静态电场问题的真解（解析解：） 3. 加权余量法－－例1 真解与近似解相同是由于尝试函数选择的刚好，通常是有差别的，如选用三角函数，但求解过程会复杂，可见尝试函数的选取是有技巧的。 4. 加权余量法求解一般化偏微分方程的归纳 一般化偏微分方程： 线性微分算子 则其余数为： 令加权余数为0，构建代数方程： 4. 加权余量法求解一般化偏微分方程的归纳 由于是线性微分算子，故微分、求和、积分次序可调换,代数方程变形： 有j个代数方程，通常等于待定系数个数 4. 加权余量法求解一般化偏微分方程的归纳 代数方程写成矩阵形式： 系数 激励 边界条件 系数矩阵n×n 待定系数矩阵、源矩阵、边界矩阵n