一阶偏微分方程求解方法PPT课件.ppt

1 偏微分方程求解 有限元法的原理 加权余量法和变分法 解析法应用范围有限 适用于理论求解 但有强烈的物理含义 常系数微分方程 某些复杂问题 很考虑根本找不到解析解2 数值法工程实际中应用广泛 复杂场域问题 但物理含义不很清楚 任何问题总可以找到数值解 数学方法 1 2 数值求解方法 2 4 1 基本思想 以偏微分方程的近似解来代替其真解 只要近似解与真解足够接近 就可以近似解作为问题的解 并满足足够的精度 2 基本方法 假设一个近似解 该解一组 形式上 简单函数的线性组合来表示 线性组合的系数就是一组待定系数然后建立一种考虑了微分方程和边界条件的关于真解和近似解间误差的目标函数F用适当的算法使得该目标函数最小化 最小化的过程就确定了待定系数 从而也就得到了问题的近似解 尝试函数 基函数 形函数 2 2 数值求解方法 2 4 目标函数最小化的目的 一方面 使得近似解最大程度接近真解 另一方面 求得构成近似解的待定系数 数学上 构成目标函数的方法很多 不同的构成方法就形成了不同的数值解法 电磁场中就常见的是 加权余量法和变分法 3 3 电磁场位函数偏微分方程的数值求解方法 加权余量法 电磁场问题总可以用位函数的偏微分方程和相应的边界条件表述 两个偏微分方程形式相同 故以电位方程的求解过程为例 磁位矢量的方程可以分解到个分量上变为标量方程 4 在求解场域内 偏微分方程的真解为 近似解为它由一组简单函数的线性组合表达 表达中有待定系数即 3 电磁场位函数偏微分方程的数值求解方法 加权余量法 加权余量法 简单函数 一般选用简单形式的函数 一旦选定就是已知的了 待定系数是真正的求解目标 问题的自由度 近似解 5 3 电磁场位函数偏微分方程的数值求解方法 加权余量法 加权余量法就是一种定义近似解与真解之间误差 即余数 并设法使其最小的方法 加权余量法误差 即余数 的定义 注意 一般余数并不表示近似解与真解间的差 场域内 加权余量法的采用拉普拉斯算子作用后的差别 即余数 来代表近似解接近偏微分方程真解的程度 问题的自由度 6 3 电磁场位函数偏微分方程的数值求解方法 加权余量法 当余数小于要求的精度时 就可以认为近似解就是偏微分方程的解 要减少余数 我们可以通过寻求适当的待定系数来实现 为有效表达减小余数的效果 还选取适当的加权函数 以使余数和该加权函数的积分为0 加权余量法 的来由 7 3 电磁场位函数偏微分方程的数值求解方法 加权余量法 加权余数的定义 加权函数的选取方法很多 如点重合 子域重合 最小二乘法 迦辽金法 效果较好的 运用较多的是迦辽金法 即 迦辽金法选取尝试函数本身为加权函数 8 3 电磁场位函数偏微分方程的数值求解方法 加权余量法 由此构建加权量法的目标函数 上述过程中 已经将偏微分方程转化为j个代数方程组 便于计算机求解 关于函数是函数 称为 泛函数 或泛函 9 3 加权余量法 例 例1 两极电容板内部电场分布问题 根据问题特点将3维问题简化为2维 进一步简化为1维 该问题是静态电场问题 偏微分方程和边界条件 10 加权余量法求解 1 选取尝试函数 构造近似解 2 结合问题 写出余数表达式 3 加权余量法 例 理论上任意选取 操作中越简单越好 11 2 结合问题 写出余数表达式 3 加权余量法 例 12 3 加权余数表达式 3 加权余量法 例 13 3 加权余数表达式 3 加权余量法 例 14 4 求解上述两个代数方程组 得到待定系数 从而确定近似解 3 加权余量法 例 加权余量法求解流程 1 选取尝试函数 构造近似解2 结合问题 写出余数表达式3 写出加权余数表达式4 令各加权余数表达式为0 得到代数方程组 解之得到待定系数 从而确定近似解 15 该静态电场问题的真解 解析解 3 加权余量法 例 真解与近似解相同是由于尝试函数选择的刚好 通常是有差别的 如选用三角函数 但求解过程会复杂 可见尝试函数的选取是有技巧的 16 17 4 加权余量法求解一般化偏微分方程的归纳 一般化偏微分方程 线性微分算子 则其余数为 令加权余数为0 构建代数方程 18 2020 1 8 19 4 加权余量法求解一般化偏微分方程的归纳 由于是线性微分算子 故微分 求和 积分次序可调换 代数方程变形 有j个代数方程 通常等于待定系数个数 20 4 加权余量法求解一般化偏微分方程的归纳 代数方程写成矩阵形式 系数 激励 边界条件 系数矩阵n n 待定系数矩阵 源矩阵 边界矩阵n 1 矩阵元素值 虽然元素值还需要积分 微分的求得 还难以借助计算机求解 但至少化为了代数方程组 通过选择合适的加权函数和尝试函数可以大大简化矩阵元素的矩阵方程 有限元方法就是如此 21 5 加权余量法的进一步优化 边界条件的处理 适当的选取加权函数 并对加权余数积分进行处理 可使某些边界条件从加权余数的表达式中消失 从而简化矩阵方程及其系数的求解 以有源静电场问题为例 帕松方程 22 由近似解表述的加权余数为 5 加权余量法求解一般化方法的进一步优化 注意余数的实质 23 通过尝试函数 简化加权余数后 5 加权余量法求解一般化方法的进一步优化 上式第一项 由格林第一定律得 降了微分阶数 等于降了近似解 尝试函数 的连续性要求 从而扩展了其选择范围 24 代入后 5 加权余量法求解一般化方法的进一步优化 由于近似解在1类边界上常数 所以此项为0 选取特殊加权函数后 两项和为0 第二类边界条件也消失了 说明已经自动满足了 25 令加权余数为0即可得到求解原微分方程的一组代数方程 5 加权余量法求解一般化方法的进一步优化 这里加权函数只有一个了 进一步 用迦辽金法 选加权函数为尝试函数本身 26 5 加权余量法求解一般化方法的进一步优化 由于是线性微分算子 故微分 求和 积分次序可调换 代数方程变形 对比简化前的代数方程 已经大大简化 关键是边界条件项全部消失 微积分计算也降阶 简化 27 5 加权余量法求解一般化方法的进一步优化 代数方程写成矩阵形式 对称矩阵 简化计算 还有积分 求和 梯度 差分 有限元将作处理 小结 简化后1 2类边界条件自动满足 尝试函数 加权函数选取 微分降阶 简化计算对称矩阵 简化计算根据情况源矩阵 边界矩阵可能为0 28 6 简化后加权余量法例2 例1中的静电场问题 变为两电极板接地 中间充满电荷 帕松方程 29 加权余量法求解 1 初选尝试函数 构造近似解 利用问题 对近似解进行简化 对尝试函数进行优化 6 简化后加权余量法例2 通过尝试函数的选取 近似解满足1类边界条件 使得1类边界条件在方程中消失 由此 尝试函数和近似解优化为 30 2 修正尝试函数 以满足1类边界条件 6 简化后加权余量法例2 31 3