静态场边值问题的解法.ppt

同理可得 在节点 0 有 节点 0 处电位的有限差分方程 这样在0点就用一个有限差分方程近似地代替了该点的泊松方程。对于ρv=0的区域, 电位函数满足拉普拉斯方程, 节点0 上的有限差分方程为 (5 - 67) 5.5.2 有限差分方程的常用数值解法 1. 迭代法 (1) 先对每一网格点任意设初始值。 ; (2) 根据有限差分方程(5 - 67), 规定域中各点每次用围绕其周围的邻接点的电位平均值作为新算出的电位值。如果已知第n次电位值, 则计算第n+1 次电位值的表示式是 ? (3) 连续逐点计算, 直至每一点前一次的电位值与后一次新算出的电位值相差满足要求的精度为止。 这时各电位值就是所求的解。 2. 松弛法(高斯—塞德尔Seidel迭代法) 上述迭代法的主要缺点是收敛速度慢。 每点算得新值后不能立即冲掉旧值, 后者在计算后面邻接点的电位值时还要用到。 如果将前述迭代法稍加改进, 在计算每个网点电位时充分利用已经得到的新值, 以使得到每点新值时立即把旧值冲掉, 这种方法称为松弛法。 松弛法计算网点电位公式是 3. 超松弛法 式中α称为迭代参数或松弛因子。取0＜α＜2 均能收敛, 但当1＜α＜2称为超松弛, 且在该范围内有一个最优松弛因子, 选择最优因子能使收敛速度大大加快。 5.5.3 边界条件的有限差分关系式 图 5 - 35 场域D的边界 1. 第一类边界条件(狄里赫利条件) 所以能直接代入有限差分方程式(5 - 69)进行计算。如果边界线和网格线重合, 则点O的差分关系式是 式中UB是在S面上点B的电位。 图 5 - 36 边界线与网格线不重合 当(M)0=0 时, 则 0 点的差分关系式是 泊松方程的有限差分形式是 若p=q=1 时, 上式简化为 2. 第二类边界条件(诺诶曼条件) 如果边界线和网格线重合, 可在边界外虚设一节点N, 则有差分关系式 或 式中h是网格的步长, G0是边界面上O点已知的电位梯度。 图 5 - 37 不同媒质分界面 的差分处理 首先将媒质εb全部换成εa, 这时0点位函数的差分表示式是 ? 然后将εa全部换成εb, 则 0 点的位又可表示为 (5 - 80) 差分式是 ? 上两式相加, 并将式(5 - 80)代入之, 得到 例 5.19 设如图 5 - 38 所示的长方形截面长导体槽, 宽边为 4h,窄边为3h, 槽的导体盖板与两侧面绝缘, 电位为 100V。试计算槽内电位分布。 图 5 - 38 导体槽截面的分格 表 5 - 2 迭代法计算 6 个内节的电位值(V) 内节点序号 迭代次数 1 2 3 4 5 6 赋初始值 0.0 0.0 0.0 0.0 0.0 0.0 第 1次 25.0 0.0 25.0 0.0 25.0 0.0 第 2 次 31.25 6.25 37.5 6.25 31.25 6.25 表 5 - 3 超松弛法计算 6 个内节点的电位值(V) 也可写为极坐标的形式: 由z=x+jy, W=u+jv确定的复平面分别为Z平面和W平面, 该两平面上的点有一一对应关系, Z平面上的图形在W平面上也有对应的图形。 这种几何上的对应关系称为映射或变换。 例如下列变换: 即 图5-26 W=z2的变换 5.4.2 解析函数, 保角变换 柯西 - 黎曼方程的直角坐标和极坐标表示式分别是 1. 由柯西 - 黎曼方程可以得到解析函数W=u+jv中的两个实变函数u和v都满足拉普拉斯方程 同理可得 上两式相加, 可得 (5-50a和b) 用类似方法, 由式(5-50a和b)可得到极坐标中u , v满足的拉普拉斯方程: ? 所以u和v是两个调和函数, 在复变函数理论中也称为共轭函数。 2. 解析函数具有保角变换的重要性质 图 5 - 27 保角变换 在z′→z0和z″→z0的极限时, 如果只考虑复数的幅角关系, 则 ? 上两式相减, 得 在Z平面上有两条曲线u(x, y)=常数和v(x, y)=常数相交于z0点, 根据梯度性质, 在z0点▽u与u(x, y)=常数的曲线垂直, 而▽ v与v(x, y)=常数的曲线垂直, 由梯度公式 5.4.3 解析函数的应用, 计算位函数和场分布 解析函数W(z)=u(x, y)+jv(x, y)中的两个实变函数u , v都满足拉普拉斯方程。而且在Z平面上u(x, y)=常数和v(x, y)=常数是两组相互正交的曲线族, 故