圆柱坐标系中恒定电磁场边值问题的数值模拟.pdfVIP

圆柱坐标系中恒定电磁场边值问题的数值模拟 卢学山，张敏，谭俊杰，王学德，刘少培，许彬 南京理工大学动力工程学院，南京 (210094) E-mail ：learnhill@ 摘 要：用商用软件 ANSYS 对恒定电磁场中的磁位势分布进行数值计算，并同理论分析的精确解 相比较。在数值模拟中，同时采用二维模型和三维模型，并比较所得的结果。通过完整的模拟过程， 展示了ANSYS 在电磁场分析中的正确性和实用性。 关键词：有限容积法，ANSYS ，电磁场 恒定电场和磁场都可引入（标量或矢量）位函数来表示。在均匀媒介中，这些位函数都满足拉 普拉斯方程或泊松方程，并且在场域的边界面上，位函数还满足相应的边界条件。所以边值问题的 求解，可归结为在给定边值条件下，对拉普拉斯方程或泊松方程的求解[1-4] 。 本文利用商用软件 ANSYS 对恒定电磁场进行数值模拟，讨论和分析计算结果。通过算例，将 数值解与精确解进行比较，并得到相似的结果，以此证实本文 ANSYS 运用的正确性与实用性。 1. 基本控制方程 对于电磁场的计算，为了使问题得到简化，通过定义两个变量把电场和磁场分离开来，分别形 成独立的电场和磁场偏微分方程，这样有利于数值求解[3] 。这两个量一个是标量电位势φ ，另一个 r 是矢量磁势A 。在柱坐标系中，标量电位势φ控制方程的表达式为， 2 2 1 ∂ ∂φ 1 ∂φ ∂φ ρ (r ) + + + 0 (1-1) 2 2 2 r ∂r ∂r r ∂ϕ ∂z ε0 其中，ρ 为电荷体密度，ε0 为真空介电常数。电场强度为， ⎛ ∂ rφ ⎞ 1 ( )r 1 ∂φ r ∂φ r E =−∇φ =−⎜ ar + 2 aϕ + az ⎟ (1-2) ⎝r ∂r r ∂ϕ ∂z ⎠ r 矢量磁势A 在柱坐标中的表达式为， r r r r A A a =+A a +A a (1-3) r r ϕ ϕ z z r 其中A 、A 、A 表示矢量磁势A 在r、ϕ、z 方向上的分量。 r ϕ z 2 2 1 ∂ ∂A 1 ∂ A ∂ A ⎛ r ⎞ r r