第四章静态场的解详解.pptVIP

三个重要的数学定理 定理1：如果函数y1 x 和y2 x 是方程y”+py’+qy 0的两个特解，则y c1y1 x +c2y2 x 也是方程的特解。 定理2：如果函数y1 x 和y2 x 满足条件：y1 x /y2 x ≠常数，则函数y1 x 和y2 x 线性无关。 定理3：如果函数y1 x 和y2 x 是方程y”+py’+qy 0的两个线性无关的特解，则y c1y1 x +c2y2 x 是方程的通解。 直角坐标系中的分离变量法 应有条件：界面形状适合用直角坐标系表示。 分析方法：先用分离变量法求通解，再重点利用边界 条件求定解。 直角坐标系中的拉普拉斯方程： 变量分离： 设 ，将其代入上式，得 除以XYZ，得 上式成立的唯一条件是三项中每一项都是常数，故可分解为下列三个方程： 其中： α,β,γ为分离常数，但不能全为实数或全为虚数 常微分方程的解： 以X”/X α2式为例，说明X的形式与α的关系 当α2 0时，则 当α2 0时，令α jkx kx为正实数 ，则 当α2 0时，令α kx，则 a，b，c，d为积分常数，由边界条件决定Y y Z z 的解和X x 类似 或 或 积分常数的大致确定方法： 若在某一个方向（如x方向）的边界条件是周期的，则分离常数是虚数，其解选三角函数； 若在某一个方向的边界条件是非周期的，则分离常数是实数，其解选双曲函数或者指数函数。其中：有限区域选双曲函数，无限区域选指数衰减函数； 若位函数与某一坐标无关，则沿该方向的分离常数为零，其解为常数。 分离变量法的求解步骤： 建立正确的坐标系，确定变量的个数； 求方程的通解； 利用边界条件求方程的定解，即求出待定系数 例4-7 横截面为如图所示的导体长槽，上方有一块与槽相互绝缘的导体盖板，截面尺寸为a×b，槽体的电位为零，盖板的电位为U0， 求此区域内的电位。 解：导体槽内为无源区，故电位满足拉普拉斯方程和边界条件： 用分离变量法求解过程： 0 设 代入 通过引入分离常数k，将二维拉普拉斯方程分解为两个齐次常微分方程。分别解这两常微分方程可得原问题的通解 解常微分方程（k值取值不同解形式不同） 当k 0时： 当k≠0时： 由于三角函数具有周期性，因此解中的分离变量k可以取一系列特定的值kn n 1,2,3…… ，即 将所有特解的线性组合起来，得到电位函数的通解 解中所有未知系数和分离变量kn由边界条件确定 已知边界条件： 由条件 1 得： 由条件 2 得： 由条件 3 得： 由条件 4 得： 将此式按傅立叶级数展开，即等式两边同乘以 再对x从0到a积分，得 等式左边 利用三角函数的正交性质 等式右边 可得 即 所以，接地导体槽内部电位分布为 即 §4.7 有限差分法 引入 分离变量法和镜像法都是求边值问题得解析解的方法。但是在许多实际问题中由于边界条件过于复杂而无法求得解析解。这就需要借助于数值法来求电磁场的数值解。有限差分法便是一种比较容易的数值解法。 原理 把求解的区域划分成网格，把区域内连续的场分布用网络节点上的离散的数值解来代替。网格划分得充分细，才能够达到足够的精度。应用有限差分法计算静态场边值问题时，需把微分方程用差分方程替代。 差分表示式 当h很小时，忽略四阶以上高次项，得 同理，有 两式相加并考虑 得： 方法 简单迭代法 步骤：先对每一网格点 设初值。然后按固定顺 序（从左到右，从下到 上），利用二维拉普拉 斯方程的有限差分形式 用围绕它的四个点的电 位的平均值作其新值，当所有点计算完后，用它们的新值代替旧值，即完成了一次迭代。然后再进行下一次迭代，直到新值和旧值之差小于指定的范围为止。为提高精度，必须进一步划分出更多更小的网格数。 塞德尔（Seidel）迭代法 上式为异步迭代法。由于更新值提前使用，异步迭代法比简单迭代法收敛速度加快一倍左右，存储量也小 为节约计算时间，对简单迭代法进行改进，每当算出一个节点的高一次的近似值，就立即用它参与其它节点的差分方程迭代。 超松弛法： 2、采用松弛因子法：α称为松驰因子，其值介于1和2之间。当其值为1时，超松驰迭代法就蜕变为赛德尔迭代法。 简单迭代法在解决问题时收敛速度比较慢，实用价值不大。实际中常采用超松弛法，相比之下它有两点重大的改进 ： 1、计算每一网格点时，把刚才计算得到的临近点的新值代入，即在计算 j,k 点的电位时，把它左边的点 j-1,k 和下面的点 j,k-1 的电位用刚才算过的新值代入。 α的选取一般只能依经验进行。但是对矩形区域，当M、N都很大时，可以由如下公式计算最佳松驰因子 其中，M、N分别为沿x，y两个方向的内节点数。对其 它形状的实际区域，α0应用其近似值，处理方法如下： 》将区域等效为近似的矩形区域，再依上式计算α0