电磁场与电磁波(西电)第4章 静态场的解

电磁场与电磁波第4章静态场的解——边值问题分类、唯一性定理与三大经典求解方法Contents本章目录电磁场理论核心方法与求解体系概览01边值问题分类与唯一性定理02镜像法原理与典型应用03分离变量法与多坐标系求解04有限差分法与数值计算CHAPTER01边值问题分类与唯一性定理建立静态场求解的数学框架：从泊松方程到三类边界条件，再到解的唯一性保证ELECTROMAGNETICFIELDTHEORY静态场的基本方程与边值问题静态场的求解本质是在给定边界条件下求解拉普拉斯方程或泊松方程。电位函数φ满足∇²φ=-ρ/ε（泊松方程），无源区退化为∇²φ=0（拉普拉斯方程），这一数学框架统一适用于静电场、恒定电场和恒定磁场。01泊松方程与静电场静电场电位φ满足∇²φ=-ρ/ε，其中ρ为自由电荷密度，ε为介电常数。该方程揭示了场源与电位分布的定量关系，是分析有源区域静电问题的核心数学工具。02拉普拉斯方程无源区域（ρ=0）泊松方程退化为∇²φ=0，这是静态场分析中最核心的偏微分方程。该方程的解具有唯一性，为边值问题的求解奠定了数学基础。03电磁对称性恒定磁场的矢量磁位A同样满足∇²A=-μJ，体现了电场与磁场在数学形式上的深刻对称性。这种对称性使得电磁场理论具有统一而优美的结构。04边值问题定解边值问题由控制方程（偏微分方程）和边界条件两部分组成，二者缺一不可才能构成完整的定解问题。常见的边界条件包括狄利克雷、诺依曼和混合边界条件。BOUNDARYVALUEPROBLEMS三类边值问题的分类静态场边值问题按边界条件的给定方式分为三类，分类决定了求解策略的选择。狄利克莱问题在整个边界Γ上给定每一点的电位值φ|Γ=f(r)，即已知边界上的函数值。典型场景：已知导体电极的形状和所加电压，求解空间电场分布。DIRICHLET·第一类诺伊曼问题在整个边界Γ上给定每一点的电位法向导数∂φ/∂n|Γ=g(r)。物理含义：已知导体表面电荷密度σ，因为∂φ/∂n=−σ/ε₀；给定导体总电量也属此类。NEUMANN·第二类混合边值问题一部分边界给定电位值，另一部分给定法向导数，两者混合出现。工程实际中最常见，如部分边界为已知电压导体、部分为已知电荷分布的介质面。MIXED·第三类UniquenessTheorem唯一性定理唯一性定理保证：在给定边界条件下，满足泊松方程或拉普拉斯方程的解是唯一的。这意味着无论采用何种方法，只要所得解同时满足控制方程和全部边界条件，即为正确且唯一的解。01定理表述在区域V内给定电荷分布ρ，在边界S上给定φ或∂φ/∂n，则V内满足泊松方程的解φ唯一确定。02证明思路假设存在两个解φ₁和φ₂，令差值u=φ₁−φ₂，利用格林第一恒等式可证u在整个区域为常数或零。03对镜像法的支撑只要构造的等效电荷系统能使边界条件满足，由唯一性定理保证该解就是真实解，无需关心求解过程。04对数值法的支撑有限差分法或有限元法的数值解只要收敛且满足边界条件，就逼近唯一真解，保证了数值方法的可靠性。UniquenessTheorem唯一性定理的证明与解的性质唯一性定理的证明依赖格林第一恒等式。第一类和第三类边值问题的解完全唯一；第二类（纯诺伊曼）问题的解可差一个常数，但电场（电位梯度）仍唯一确定，不影响物理结果。格林第一恒等式∫ᵥ(φ∇²ψ+∇φ·∇ψ)dV=∮ₛφ(∂ψ/∂n)dS是连接体积分与面积分的关键数学工具第一类边值问题：边界上φ已知，差值函数u在边界为零，由格林恒等式推出∫|∇u|²dV=0，故u为常数且等于零，解完全唯一第二类边值问题：边界上∂φ/∂n已知，差值函数的法向导数为零，推出u为常数但不一定为零，即解可差一常数，电场仍唯一三类边值问题解的性质对比边值问题类型边界条件解的唯一性物理说明第一类（狄利克莱）给定φ|S完全唯一已知电压分布第二类（诺伊曼）给定∂φ/∂n|S差一常数电场E=−∇φ仍唯一第三类（混合）部分给φ，部分给∂φ/∂n完全唯一最常见工程场景第一类和第三类问题解完全唯一，第二类问题解差一常数但电场唯一StaticFieldMethods静态场求解方法总览静态场求解方法分为解析法和数值法两大类。解析法包括镜像法和分离变量法，适用于规则边界和简单几何；数值法以有限差分法为代表，适用于复杂边界和工程实际。镜像法利用等效镜像电荷替代导体表面感应电荷，将有限区域问题转化为无限大均匀空间的叠加问题，物理概念清晰直观。适用于边界为平面、球面、圆柱面等规则几何的场景，数学推导简洁优美，是理解边值问题的经典入门方法。物理直观性强解析解形式简洁解析法分离变量法在直角、圆柱、球坐标系中将偏微分方程分离为多个常微分方程，逐一求解本征值问题后叠加得到通解。适用范围最广的解析方法，可处理矩形槽、圆柱腔、球形区域等典型几何结构，是电磁理论课程的核心内容。适用范围广泛数学体系完备解析法有限差分法将连续求解区域离散为网格节点，用差分近似替代微分算子，转化为大规模线性方程组迭代求解。不受边界几何限制，适合计算机编程实现，是现代电磁仿真软件的理论基础，工程应用价值极高。边界适应性强易于编程实现数值法CHAPTER02镜像法原理与典型应用用虚设的等效电荷替代边界感应电荷，将复杂边界问题转化为自由空间的简单叠加ElectromagneticMethods镜像法的基本原理镜像法用位于待求场域外的虚设等效电荷替代导体表面感应电荷，将有限区域问题转化为无限大均匀空间中真实电荷与镜像电荷的叠加。理论依据是唯一性定理：只要边界条件满足，解即为正确。核心思想在待求场域边界外放置镜像电荷，替代导体表面感应电荷的效应，使边界条件自然满足MirrorCharge理论依据唯一性定理保证——只要构造的解满足控制方程和边界条件，就是唯一正确解，无需关心求解途径Uniqueness三个关键约束①镜像电荷必须在待求区域之外②将空间处理为无限大均匀媒质③实际电荷与镜像电荷共同满足原边界条件3Constraints方法优势与局限物理直观、计算简洁，但仅适用于平面、球面、圆柱面等少数规则边界，复杂几何需借助其他方法ScopeELECTROMAGNETICS·METHODOFIMAGES点电荷对无限大接地导体平面的镜像点电荷q位于接地导体平面上方d处时，其镜像电荷为-q，位于平面下方对称位置d处。上半空间任意点电位为真实电荷与镜像电荷产生电位的叠加，导体平面电位恒为零，边界条件自然满足。⚡镜像法物理模型真实电荷+q位于z=+d导体平面接地，位于z=0镜像电荷−q位于z=−d点电荷q与镜像电荷−q关于导体平面对称分布01物理模型：点电荷q位于z=d处，无限大接地导体平面位于z=0平面，待求场域为z>0的上半空间02镜像电荷设置：在z=-d处放置电荷量为-q的镜像电荷，与真实电荷关于导体平面对称分布03电位表达式：φ=q/(4πε₀)·(1/r₁−1/r₂)，其中r₁、r₂分别为场点到真实电荷和镜像电荷的距离04边界条件验证：在z=0平面上r₁=r₂，故φ=0恒成立，与接地导体平面电位为零的边界条件完全吻合InducedCharge&Field导体平面的感应电荷与电场分布导体表面的感应电荷密度在正对点电荷处最大，向外递减。总感应电荷积分结果恰好等于镜像电荷-q，验证了镜像法的物理自洽性。01感应面电荷密度σ=−ε₀(∂φ/∂z)|z=0=−qd/[2π(r²+d²)3/2]，负号表示与q异号02面电荷密度在r=0（正对点电荷处）取最大值−qd/(2πd³)，随r增大按(r²+d²)−3/2衰减03总感应电荷Q=∫₀∞σ·2πrdr=−q，恰好等于镜像电荷，证明镜像电荷完整等效了感应电荷04上半空间电场线从点电荷q出发，终止于导体平面的感应负电荷，电场线与平面垂直（导体表面边界条件）Electrostatics·MethodofImages点电荷对接地导体球面的镜像接地导体球面（半径a）外距球心d处的点电荷q，其镜像电荷q'=-qa/d位于球心到q连线上距球心a²/d处。镜像位置与真实电荷关于球面满足反演关系d·d'=a²，球面电位恒为零。01镜像电荷大小q'=−qa/d，电荷量与点电荷q成正比，与距离d成反比，负号表示异号q'=−qa/d02镜像电荷位置d'=a²/d，位于球心到点电荷的连线上且在球内，与真实电荷关于球面成反演对称d'=a²/d03反演几何关系d·d'=a²，真实电荷和镜像电荷关于球面互为反演点，这是球面镜像的核心几何特征d·d'=a²04不接地孤立导体球需额外在球心放置第三个镜像电荷q''=qa/d以保证球面总电荷为零，电位为常数q''=qa/dElectrostatics·MethodofImages导体圆柱面与介质平面的镜像镜像法可推广至圆柱面和介质分界面。导体圆柱面的镜像与球面类似，满足反演关系；介质平面镜像则需分别处理两侧区域，因极化效应导致镜像电荷量不同于导体情况，体现了镜像法的灵活性。01·CYLINDRICALCONDUCTOR导体圆柱面的镜像线电荷λ距接地导体圆柱轴线d处，镜像线电荷λ'=−λ位于轴线另一侧d'=a²/d处与球面镜像完全类似的反演关系，但用线电荷替代点电荷，适用于二维场问题d'=a²/d02·DIELECTRICINTERFACE介质平面的镜像介质分界面存在极化面电荷，需分别对待求区域两侧设置不同的镜像电荷来处理折射和反射效应区域1等效电荷q'=q(ε₁−ε₂)/(ε₁+ε₂)，区域2为q''=2ε₂q/(ε₁+ε₂)，体现介质极化程度的影响ε₁·ε₂ELECTROMAGNETICS·METHODREVIEW镜像法总结与适用性评估镜像法是静态场求解中物理直觉最强的方法，通过等效替代将边界问题化为自由空间叠加。适用范围限于平面、球面、圆柱面等少数规则边界，复杂几何需转向分离变量法或数值方法。01核心解题步骤：①确定待求区域和边界②找镜像电荷的位置和大小③验证边界条件④叠加求电位和电场4步02常见错误：镜像电荷放入待求区域内、忘记验证边界条件、混淆接地与不接地导体的处理方式3类03与分离变量法的互补：镜像法适合简单边界和少量点电荷，分离变量法适合连续边界条件和复杂几何互补镜像法在不同场景的适用性评估导体平面镜像在直观性和简便性上更优，但两者在几何适用范围上都有限CHAPTER03分离变量法与多坐标系求解将偏微分方程分解为常微分方程组，在直角、圆柱、球坐标系中系统求解拉普拉斯方程MathematicalPhysics分离变量法的基本思想分离变量法假设解可表示为各坐标变量独立函数的乘积，代入拉普拉斯方程后分离为多个常微分方程。分离常数的正负号由边界条件决定，这是方法的核心技巧。01基本假设设φ(x,y,z)=X(x)·Y(y)·Z(z)，将多变量函数分解为各坐标变量独立函数的乘积形式02分离过程代入∇²φ=0后除以XYZ，各项仅含单一变量，令其分别等于分离常数且所有常数之和为零03分离常数选择根据边界条件决定正负号——齐次边界方向取负常数得三角函数解，非齐次方向取正常数04一般解与特解常微分方程的通解乘以叠加系数构成一般解，利用边界条件确定系数得到满足全部条件的特解BOUNDARYVALUEPROBLEM直角坐标系中的分离变量法直角坐标系中分离变量法将二维拉普拉斯方程分解为X和Y方向的常微分方程。通过齐次边界条件确定本征值和本征函数，再利用非齐次边界条件的傅里叶展开确定叠加系数，是解析求解的标准化流程。01典型问题：无限长矩形槽（0≤x≤a,0≤y≤b），三面接地φ=0，顶面y=b处φ=V₀，求槽内电位分布02分离过程：设φ=X(x)Y(y)，得X''+k²X=0和Y''-k²Y=0，x方向取正弦解、y方向取双曲正弦解03本征值确定：由x=0和x=a处φ=0得kₙ=nπ/a（n=1,2,3...），本征函数为sin(nπx/a)04系数求解：利用y=b处φ=V₀做傅里叶正弦展开，得系数Cₙ=4V₀/(nπ·sinh(nπb/a))，n为奇数MATHEMATICALPHYSICS圆柱坐标系中的分离变量法圆柱坐标分离变量产生三个方向的常微分方程：z方向为指数/三角函数，φ方向为周期三角函数，r方向为贝塞尔方程。贝塞尔函数是处理圆柱形边界问题的核心工具。01方程分离拉普拉斯方程在圆柱坐标(r,φ,z)中分离为三个ODE：R(r)、Φ(φ)、Z(z)各自满足独立的常微分方程R·Φ·Z02φ方向解Φ(φ)=A·cos(mφ)+B·sin(mφ)，m为整数以保证方位角的2π周期性条件2πPeriodic03r方向·贝塞尔含z变化时R(r)满足贝塞尔方程，解为m阶贝塞尔函数Jₘ(kr)和诺伊曼函数Nₘ(kr)的线性组合Jₘ+Nₘ04r方向·幂函数无z变化时退化为R(r)=Crᵐ+Dr⁻ᵐ的幂函数形式，适用于无限长圆柱或二维圆形区域问题Crᵐ+Dr⁻ᵐELECTROMAGNETICS·EXAMPLE例题：均匀外场中的导体圆柱无限长导体圆柱置于均匀外电场E₀中，圆柱外电位等于原均匀场电位与等效偶极子场电位的叠加。感应电荷使圆柱附近电场增强、远区恢复均匀，体现了边界对场的局部扰动效应。01问题设定半径为a的接地导体圆柱置于均匀外场E₀（沿x方向），求r≥a区域的电位分布02边界条件①r→∞时φ→-E₀r·cosφ（恢复均匀场）②r=a时φ=0（接地导体面等电位）03求解结果φ=-E₀(r-a²/r)cosφ，可分解为均匀场项-E₀r·cosφ和偶极子项E₀a²cosφ/r04物理意义导体圆柱的感应电荷等效为一个线偶极子，表面电场增强为2E₀（场强集中效应）SPHERICALCOORDINATES球坐标系中的分离变量法球坐标分离变量产生径向幂函数、角向勒让德函数和方位角三角函数三类解。轴对称情况下退化为勒让德多项式Pₙ(cosθ)，是处理球形边界问题的核心工具，在物理和工程中有广泛应用。分离变量拉普拉斯方程在球坐标(r,θ,φ)中分离为径向方程、极角方程和方位角方程三个常微分方程。三个ODE方位角解Φ(φ)Φ(φ)=A·cos(mφ)+B·sin(mφ)，m为整数保证周期性；轴对称时m=0，解与φ无关。m=0轴对称极角解Pₙᵐ(cosθ)θ方向满足勒让德方程，解为缔合勒让德函数Pₙᵐ(cosθ)；轴对称时退化为勒让德多项式Pₙ(cosθ)。勒让德多项式径向解R(r)R(r)=Arⁿ+Br⁻⁽ⁿ⁺¹⁾，正幂次项在球外发散，负幂次项在球心发散，需根据区域选取。rⁿ/r⁻⁽ⁿ⁺¹⁾Electrostatics·BoundaryValue例题：均匀外场中的导体球接地导体球在均匀外场E₀中的解表明：球外电位为均匀场与等效点偶极子场的叠加，偶极矩p=4πε₀a³E₀。球面极处电场增强至3E₀，揭示了球形导体表面的场强集中效应。01问题设定半径a的接地导体球置于均匀外场E₀（沿z轴），轴对称条件下φ与方位角无关，仅含Pₙ(cosθ)项02边界条件①r→∞时φ→−E₀r·cosθ②r=a时φ=0，利用P₁(cosθ)=cosθ确定只需n=1项03解的结果φ=−E₀(r−a³/r²)cosθ，分解为均匀场项−E₀r·cosθ和偶极子项E₀a³cosθ/r²04球面电场E_r(a,θ)=3E₀cosθ，极处(θ=0)场强最大为3E₀，是高压绝缘设计中球面放电的理论基础3E₀COMPARISON三种坐标系分离变量法对比分离变量法的核心在于选择与边界几何匹配的正交坐标系。直角、圆柱、球坐标系分别对应矩形、圆柱、球形边界，各自的本征函数体系决定了级数展开的形式。坐标系选择原则：边界形状应与坐标面重合——矩形边界选直角坐标，圆柱边界选柱坐标，球面边界选球坐标边界↔坐标面重合本征函数的完备性：三种坐标系的基函数各自构成完备正交系，保证任意边界条件都可用级数展开来匹配完备正交系→级数收敛三种坐标系分离变量法核心要素对比坐标系分离变量本征函数适用边界直角坐标x,y,z三角函数sin/cos、指数函数e±kx矩形槽、矩形波导、平板结构圆柱坐标r,φ,z贝塞尔函数Jm(kr)、三角函数圆柱波导、同轴电缆、圆柱腔球坐标r,θ,φ勒让德多项式Pn(cosθ)、球谐函数球形导体、球腔天线、原子轨道三种坐标系的基函数体系各不同，但都构成完备正交系，保证级数解的收敛性CHAPTER04有限差分法与数值计算将连续区域离散为网格节点，用差分近似替代微分，实现静态场的计算机数值求解NUMERICALMETHODS有限差分法的基本原理有限差分法用泰勒展开将微分算子离散化为差分格式。二维拉普拉斯方程在正方形网格上离散后，每个节点的电位等于四个相邻节点电位的算术平均，体现了调和函数的中值性质。差分近似利用泰勒展开将导数离散化，如∂²φ/∂x²≈[φ(i+1,j)−2φ(i,j)+φ(i−1,j)]/h²，h为网格步长。泰勒展开拉普拉斯离散φ(i,j)=[φ(i+1,j)+φ(i−1,j)+φ(i,j+1)+φ(i,j−1)]/4，即中心节点等于四邻均值。四邻均值物理诠释离散形式体现了调和函数的中值定理——无源区域内任一点电位等于其周围电位的平均值。中值定理截断误差中心差分格式的截断误差为O(h²)，精度随步长减小呈二次方提高，加密网格可提升计算精度。O(h²)FINITEDIFFERENCEMETHOD差分方程建立与边界条件处理有限差分法将偏微分方程转化为大型线性代数方程组Ax=b。内部节点用标准五点差分格式，边界节点根据边界条件类型特殊处理。网格密度决定了精度与计算量的平衡。01网格划分N×Mnodes将求解区域均匀离散为N×M个网格节点，步长h越小精度越高但未知数按N×M增长。合理的网格划分是数值求解的第一步，直接影响后续方程组的规模与求解效率。均匀离散·步长控制02内部节点方程Ax=b每个内部节点产生一个线性方程，N个内部节点构成N元线性方程组Ax=b。五点差分格式将拉普拉斯方程离散为相邻节点电位的加权平均关系。五点格式·稀疏矩阵03第一类边界处理Dirichlet边界节点电位已知，直接代入相邻内部节点方程中作为已知量，减少未知数个数。Dirichlet条件是最常见的边界类型，物理上对应导体表面或接地电极。已知电位·降维处理04不规则边界处理不等距差分靠近曲边界的节点采用不等距差分格式或插值方法，保证边界条件的准确施加。实际工程中的复杂几何形状需要特殊的数值处理技术来维持计算精度。曲线拟合·局部加密NumericalMethods差分方程组的迭代求解方法大规模差分方程组通常采用迭代法求解。从雅可比迭代到高斯-赛德尔再到逐次超松弛（SOR），每步改进都显著提升收敛速度。基本迭代法雅可比迭代每次迭代用上一轮全部旧值计算新值，φ_new=四邻旧值之和/4，实现简单但收敛慢φ/4基本迭代法高斯-赛德尔按顺序更新节点，已更新节点立即使用新值，收敛速度约为雅可比法的两倍2×加速收敛方法逐次超松弛法引入松弛因子ω∈(1,2)，新值=旧值+ω×修正量，最优ω可大幅减少迭代次数ω∈(1,2)加速收敛方法收敛判据相邻两次迭代最大节点电位差小于预设精度ε时停止迭代计算εNUMERICALCASESTUDY有限差分法算例：方形区域电位分布以正方形区域（三边接地、顶边V₀=100V）为例，10×10网格离散后用SOR迭代法求解。迭代约50次收敛，最优松弛因子可将迭代次数降至20次左右，展示了有限差分法的完整计算流程和实用性。问题设定：1m×1m正方形区域，三边φ=0，顶边φ=100V，网格步长h=0.1m，内部81个未知节点初始条件：所有内部节点电位初始值设为0V，边界节点按已知条件固