干输电线路设计这些年,我遇到最频繁的问题就是建设单位在可研阶段眼巴巴地问:线路要过这片规划区,地面电场到底多大?雨天会不会嗡嗡响?验收会不会照样超标?这些问题在杆塔立起来之前,唯一的回答手段就是电磁场耦合仿真。把工频电场、工频磁场、无线电干扰、可听噪声一起算出来,拿着一沓曲线去跟环评、规划乃至居民解释,比任何口头承诺都管用。
这篇内容就是围绕输电线路电磁环境的电磁场耦合仿真展开的。适用于刚接触线路电磁环境评估的工程师、做输变电工程环评的技术人员,以及想在毕业论文里用仿真的研究生。我会把从算法选型、建模过程、限值判断到实测校验的完整思路都过一遍,按我实际干活的经验来讲,不绕弯子。
1. 为什么要同时算电场、磁场、无线电干扰和可听噪声四个量
1.1 四个物理量,三种物理来源
很多新手上来就问"输电线路电磁环境到底仿真什么"。答案是四个量:工频电场、工频磁场、无线电干扰、可听噪声。这四个量看起来一大堆,归类之后其实只有三条物理链。
工频电场主要来自导线对地的电压。导线带电后,在周围空间建立电场,电压越高、导线离地越近,地面附近的电场强度就越大。决定它的不是你传送了多少功率,而是系统电压等级和几何布置。
工频磁场则完全换了一个来源,电流是主角。导线流过负荷电流时,周围产生磁场,电流越大,磁场越强。这解释了为什么同样的500kV线路,空载和满载时磁场可以差出好几倍,而电场几乎没有变化。
无线电干扰和可听噪声则要复杂一些,它们都不是导线"正常"工作时的直接产物,而是高电压导线表面发生局部放电——即电晕放电——时的副产品。导线表面电场强度超过起晕场强,空气被局部击穿,产生脉冲电流,辐射出频谱很宽的干扰信号;同时空气分子在放电过程中被激发、振动,产生人耳能听到的嘶嘶声、噼啪声,特别是在雨水、雾天这样的高湿天气下更明显。
1.2 耦合仿真在工程预测里的不可替代性
这四个量实际上是相互牵连的,这就是标题里"耦合"两个字的含义。导线表面电场强度由电压和几何参数决定;表面电场强度超过临界值,电晕放电开始出现;电晕放电强度直接决定无线电干扰水平和可听噪声水平。也就是说,你改变导线型号、分裂数、相间距离、对地高度等任何一个几何参数,电场会变,磁场相应也会变,电晕放电水平跟着变。这四个量不是四个孤立计算,而是一条物理链上的多个输出端。
这种耦合关系在规划阶段的预测价值非常大。运行中的线路你可以拿仪器去测,但规划中的线路,设计阶段就必须给出一份"未来电磁环境预测报告"。测量的局限在于只能回答"现在、这里、这种天气下是多少",而仿真可以回答"如果线路加高3米会怎样""如果采用4分裂而不是6分裂会怎样""如果跨越段改到那边的山脊上会怎样"。建设单位要决策,环评要依据,方案要优化,这些都必须靠仿真来给答案。
所以做输电线路电磁环境评估,前提就是把四个量放进一个模型里统筹计算,而不是今天算电场、明天翻规范套一个磁场经验值、后天再翻手册查个干扰公式。割裂计算的问题是参数不一致,几何条件各有各的说法,最后结果根本对不上。耦合仿真至少能保证:所有物理量的计算基于同一套导线参数、同一个塔型、同一段弧垂。
| 物理量 | 主导来源 | 与运行参数的关系 | 主要影响因素 |
|---|---|---|---|
| 工频电场 | 导线对地电压 | 几乎不随负荷变化 | 电压等级、对地高度、相间距离、相序 |
| 工频磁场 | 导线负荷电流 | 随负荷电流线性变化 | 电流大小、导线空间位置、相位差 |
| 无线电干扰 | 电晕放电 | 天气敏感,雨天显著增加 | 导线表面电场强度、导线直径、天气 |
| 可听噪声 | 电晕放电 | 天气敏感,雨天显著增加 | 导线表面电场强度、分裂数、子导线直径 |
2. 三种主流仿真路线的选型逻辑与我的建议
2.1 路线一:模拟电荷法加毕奥-萨伐尔自编脚本
第一种路线是完全自编程序,电场上使用模拟电荷法,磁场上使用毕奥-萨伐尔定律离散叠加,再用激发函数法处理电晕产生的干扰和噪声。语言一般选MATLAB或者Python,网上有大量公开代码可以改。
这条路线的最大优势是完全可控。每一行代码都知道在算什么,参数怎么改都不心疼,适合做方案扫参——比如要比较导线高度从12米变到20米场强的变化曲线,脚本循环跑一遍,一张图全出来。商业软件当然也能做参数扫描,但设置一个参数扫描工况的时间,我的经验是比写脚本一次循环慢得多。
缺点也明显。第一,二维平面下模拟电荷法算平直段线路没有问题,但遇到大跨越、转角塔、线路从山谷穿出来这种空间走向复杂的情况,二维模型就撑不住了,需要写成三维。第二,大地不是理想导体时,土壤分层、接地方式的影响很难在简单脚本里考虑周全。第三,干扰和噪声的经验公式系数各家有各家的版本,脚本里一旦选错了系数,结果可能与实测差出10dB以上。
2.2 路线二:通用电磁场商业软件
第二类是用通用电磁场仿真软件,比如COMSOL、ANSYS Maxwell、CDEGS。这类软件的特点是建模能力强、材料属性丰富、后处理图形漂亮,适合处理复杂地形下的三维问题,也能把人体、管道、构架这类实际物体放进模型里,计算感应电压和人体接触电流。
但这类软件用在架空线路上有个典型问题:模型边界难取。架空线路是半无限域问题,导线延伸到视线尽头,档距几百米,场点只关心线路附近几十米。你用有限元法去剖分,就把一块巨大的空气域全部剖了,虽然对称性可以简化,但网格量大,调一次参数重剖一次网格,效率非常低。CDEGS这类基于矩量法的软件会好一点,但学习曲线也不平缓。
另外还有一个"软件陷阱":商业软件模型一旦建得复杂,内部参数几十个,哪个物理场、哪个边界条件的设置一旦不对,结果就是错的,而新手几乎很难察觉。我见过不止一个同学用COMSOL算出来的电场云图整片都是均匀的,原因是激励相的电压没填对,相位设成了零。
2.3 路线三:线路设计专用软件的内置模块
第三类是在线路设计软件——比如PLS-CADD,或者国内各大电力设计院自研的线路设计软件中,使用内置的电磁环境计算模块。这些模块大多内置了成熟算法,输入杆塔、导线、弧垂数据之后可以一键出结果。
这类工具的优势在于数据链完整。它直接从线路平断面和杆塔明细表里读参数,弧垂、档距、导线型号都是现成的,不需要你手工输入和建模,对设计院做标准化交付来说最方便。缺点是可定制性差,你想要调整某个经验系数、改一种导线排列方式、算一个非标工况,模块不见得支持。同时算法是黑盒,写报告时"计算过程"这部分往往说不清楚用了什么方法。
2.4 选型对比表
| 路线 | 精度 | 灵活性 | 上手成本 | 典型场景 | 成本 |
|---|---|---|---|---|---|
| 自编脚本 | 中高 | 极高 | 中 | 方案扫参、论文分析、快速评估 | 仅时间成本 |
| 通用商业软件 | 高 | 中高 | 高 | 复杂地形、人体耦合、精确建模 | 高 |
| 线路专用软件 | 高 | 低 | 低 | 标准线路设计、交付文件 | 中高 |
2.5 我的组合打法
我自己实际做项目时不太押注单一工具。常规做法是:初步方案阶段,用自编脚本快速扫参,确定导线高度范围、相序排列和走廊宽度;到了可研或者环评阶段的关键工况,再用商业软件做精确建模复核一遍;交付时如果需要线路设计标准的成果格式,就把参数导入线路设计软件出正式报告。
这个组合打法看上去多了一道工序,实际上效率最高。脚本负责快速找趋势,商业软件负责精确校核,两者结果互相印证,任何一个环节出bug都能及时发现。纯靠脚本,审核工程师会质疑算法可靠性;纯靠商业软件,方案比选阶段会累死。两条腿走路,才是这行的常态。
3. 建模实操:从原始参数表到场强云图
3.1 第一步:整理一份合格的输入参数表
建模的第一步永远不是打开软件,而是把参数表整理清楚。我吃过亏,拿着一个参数说导线4分裂,结果算出来场强死活不对,最后发现分裂间距450毫米记成了450厘米。这类低级错误看着可笑,实际工作中极其常见。
一份合格的输入参数表至少要包括以下内容:
- 系统标称电压和最高运行电压。很多新手用标称电压算电场,结果偏小。工程上一般取最高运行电压,比如500kV系统取550kV来计算。
- 导线型号和几何参数。以JL/G1A-400/35为例,需要记下的是计算直径26.8毫米(不是包括铝股间隙的名义直径,是计算直径),直流电阻0.07389欧姆每千米,这两个参数一个影响电场、一个影响磁场和损耗。
- 分裂结构和分裂间距。4分裂、分裂间距450毫米,这是常见的500kV配置。注意分裂间距也直接影响等效半径。
- 杆塔上各相导线的挂点高度、水平偏移、相序排列。三相水平排列、垂直排列、三角排列,结果差异巨大。
- 代表档距下的弧垂值。弧垂最低点的导线高度,才是某个断面实际的导线离地高度。
- 地线参数。地线的直径、高度也要,虽然地线对地面场强贡献不大,但会影响无线电干扰的计算裕度。
我用一个典型工况举例:500kV单回线路,三相水平排列,相间距12米,导线4分裂,子导线直径26.8毫米,分裂间距450毫米,最低弧垂点导线对地高度15米,计算地面1.5米高度处的电场分布。这个例子在后面会反复用到。
3.2 第二步:建立二维计算模型,用模拟电荷法求电场
架空线路平直段内,电场问题可以简化成二维问题来求解。二维简化基于一个基本事实:平直段线路沿档距方向的电场变化远小于垂直档距方向的电场变化,取一个横断面计算就足够代表整段线路。
模拟电荷法的思路是这样的:导线表面是等位面,电位就是导线电压。我们在导线内部放置若干虚拟线电荷,虚电荷的数值未知,但其共同作用在导线表面产生的电位必须等于给定的导线电压。用这个条件建立起代数方程组,解出电荷量,之后空间中任意一点的电位和电场就是这些电荷以及它们在地中的镜像电荷的叠加。
步骤说通俗一点:
- 对每一根子导线,在其轴心位置放置一根无限长线电荷。导线有4根分裂子导线,就放4根线电荷。
- 大地用镜像法处理,在每根线电荷关于地面的镜像位置放一个等量异号的镜像电荷,这样大地电位为零的边界条件自动满足。
- 写出电位系数矩阵。矩阵的第i行第j列元素表示第j根电荷对第i个匹配点(取各子导线表面)的电位贡献系数。自项即电荷本身作用于自身表面的电位系数,要用分裂导线的等效半径来计算。
- 求解线性方程组,得到每根线电荷的电荷密度。
- 空间任意一点的电位和场强,由所有实电荷和镜像电荷叠加得到。
有一个细节非常关键:分裂导线的等效半径不能直接用子导线半径。四分裂导线,每根子导线半径只有13.4毫米,但作为一个整体对外呈现的电场特性,应该按几何平均半径来处理,等效半径约为分裂间距半径乘以一个与子导线数、子导线半径有关的系数。忽略这一点,导线表面场强会算得离谱地偏大,电晕计算也随之出错。
得到电荷分布后,地面附近的电场强度分布就水到渠成了。关心的是离地1.5米高度处的场强绘制曲线,这个高度对应人体躯干所处的空间位置,是工程和环保关注的核心高度。
3.3 第三步:用毕奥-萨伐尔离散叠加算磁场
磁场计算相对电场要直观一些,核心是毕奥-萨伐尔定律。无限长直导线在空间中一点产生的磁感应强度大小等于真空磁导率乘以电流除以两倍圆周率乘以该点到导线的距离,方向由右手定则确定。
实际线路没那么理想。三相导线中流过的电流幅值相等、相位各差120度,每一相产生的磁场在空间某一点的方向都不同,必须矢量叠加。叠加时还要把三相电流的复数相位考虑进去,最终算出来的磁场是一个旋转矢量,工程上一般关心其轨迹长轴,也就是有效值最大的那个方向。
对于弧垂和档距的影响,二维近似下可以把导线视为无限长直线,这在距杆塔较远的档距中央区域是够用的。但如果杆塔附近、跳线区域或者大跨越段,导线走向明显偏离水平直线,就要把导线沿弧垂曲线离散成若干小段,对每一小段做三维毕奥-萨伐尔积分再叠加。
磁场和电流的强耦合关系值得单独强调。电流是随负荷变化的运行参数,所以磁场的预测结果必须在报告中注明对应的负荷工况。环评阶段一般取线路设计输送容量的典型电流来算,再给出重载工况的上限包络线。
3.4 第四步:把电晕放电的链条接上来
电场的计算没有结束在电场本身,它还要为无线电干扰和可听噪声提供输入。这一步是整个耦合仿真中概念上最绕的一环。
电晕放电的强弱取决于导线表面电场强度。当表面场强超过起晕场强,导线附近空气发生局部电离。工程上常用皮克公式估算导线的起晕场强,对于圆柱形导线,起晕场强约等于30倍的空气相对密度乘以一个与导线半径有关的系数。半径越小、表面越脏(比如附着雨滴或者灰尘),起晕越容易发生。
无线电干扰的计算思路是把电晕放电产生的脉冲电流当作源,用激发函数法求解干扰场强在空间中的分布。激发函数本质上是描述电晕源强度的经验函数,不同研究者(比如国际大电网会议CIGRE、IEEE工作组)给出的系数差异较大。计算结果通常表示为某个参考点(我国标准取边相导线投影外20米、频率0.5MHz)的干扰场强,单位是dB(μV/m)。
可听噪声的工程计算也不复杂,各大设计院一般用经验公式估算,输入是导线表面最大电位梯度和子导线直径、分裂数,输出是雨天或者晴天的A计权声压级。BPA公式、Cheung公式等方法都是可以用的,工程上差别不算大。计算的关键在于确认表面电位梯度是不是按分裂导线等效半径算的,这个参数一错,噪声误差会到5dB以上。
3.5 第五步:结果输出与判读技巧
仿真结果最常用的输出形式有两种:一种是地面不同高度处的场强分布曲线,横坐标是距离线路中心的水平距离,纵坐标是场强;另一种是走廊断面云图,把电场或磁场的空间分布变成彩色图,供报告和汇报使用。
判读结果时有一个经验:工频电场的最大值通常出现在边相导线投影位置附近或相间区域,而不是线路正下方正中位置。如果你算出来的电场曲线最高点不在边相投影附近,那就要回头检查相序、相间距离或者导线坐标是否输入错误。
另外,对于500kV典型线路,导线最低点离地15米时,地面1.5米高度处的最大场强通常落在6到10kV/m之间;一旦超过这个范围,要么导线不够高,要么相间距离和排列方式有特殊情况,需要进入优化环节。这些数值区间可以作为快速自查的参考标准。
4. 四个约束条件,直接决定线路方案能否落地
4.1 工频电场:公众曝露4kV/m
工频电场的控制限值是我国电磁环境管理中最被关注的指标。按GB 8702的要求,公众曝露控制限值为工频电场强度4kV/m。输变电工程环境影响评价中,线路经过居民区、学校、医院等敏感目标时,线下及建筑物所在位置的地面1.5米高度处电场强度要满足这个限值。
工程上还有一个额外概念是"线路走廊"的划分。电场强度较高的区域被划入线路保护区,这就是为什么线路不能直接跨过居民屋顶,必须预留足够距离。实际操作中,电场超标一般通过加高杆塔来缓解——高度每增加3米,线下最大场强能下降25%到40%,效果很显著,但代价是铁塔钢材和基础造价同步上升。
4.2 工频磁场:公众曝露100μT
磁场限值同样来自GB 8702,工频磁感应强度公众曝露控制限值为100μT。坦白说,实际运行线路要达到这个限值并不容易——标准工况下即使满载,线路附近的磁感应强度一般也就几十微特斯拉。所以磁场在环评中通常不是制约因素。
但对仿真来说,磁场不能因为"容易达标"就不算。在线路与地下管线、光缆、通信线交叉跨越时,磁场和静电感应引起的干扰问题依然是要拿出计算依据的。我做过一个项目,线路跨越一段天然气管道,对方要求提供管道沿线任一点的磁感应强度分布,施工阶段还要求数值不要超过某个阈值,这时候磁场计算就是唯一的依据。
4.3 无线电干扰:0.5MHz处53dB
无线电干扰的参考频率是0.5MHz,参考位置是边相导线投影外20米处,限值是53dB(μV/m),这是针对110kV到500kV高压交流架空送电线路的标准要求。考量的本质是保护中波广播频段的接收质量。
这个指标很有意思,因为它与电场计算结果直接挂钩——导线的表面电场强度越高,电晕放电越强,无线电干扰越严重。控制它最有效的办法不是加高杆塔,而是改变导线表面电场状态:增大子导线直径、增加分裂数、加大分裂间距,都可以降低表面场强。这个规律在方案优化时特别有用:电场不达标要加高,干扰不达标要换导线,两个措施的指向完全不同。
4.4 可听噪声:雨天50dB(A)
可听噪声的控制指标一般取边相导线投影外20米处,雨天测量,A计权声级不大于50dB(A)(对于500kV线路)。取雨天作为控制工况是因为导线表面附着水滴后电晕强烈,噪声会显著高于晴天。
判断可听噪声达不达标,同样看导线表面场强。实际项目里我发现一个规律:导线表面最大场强控制在18kV/cm左右,500kV四分裂线路的雨天噪声通常能压在50dB(A)以内;一旦超过20kV/cm,噪声超标概率就大幅上升。所以这个指标本质上也是在约束你的导线选型和分裂方案。
4.5 超标时的调整手段
这四个限值并非同时起作用,具体项目里哪个卡脖子要算过才知道。但调整手段可以分类整理:
- 加高杆塔:对降低地面电场最有效,对磁场也有帮助,对电晕噪声和干扰帮助不大。
- 更换导线截面或分裂形式:对降低表面场强、减少电晕最有效,间接改善噪声和干扰,对地面电场影响有限。
- 优化相序排列:可以改变地面电场分布形态,把最大值移出敏感区域,但改变不了场强量级。
- 增加走廊宽度或调整线路路径:属于规划层面手段,把敏感建筑物移出高场强区域,最彻底但代价最大。
实际方案调整永远是组合拳,先看是哪个指标不达标,再对症下药。四个指标连算出来的优势就在这里:你一次就能看清短板在哪个参数上,不用做完电场发现超标加高,加完再做噪声,又发现噪声因为加高没改善,再回头换导线,白走一大圈。
5. 实测与仿真对不上的六个原因
仿真做得再漂亮,到了现场验收阶段都得接受实测的检验。我踩过的坑和见过的对不上案例,归纳起来跑不出六个原因。
5.1 天气状态与仿真工况不一致
这是最容易被忽略也最伤脑筋的一条。无线电干扰的可听噪声实测,你是在晴天去测的还是雨天测的,结果可以差出10到20dB。仿真模型里如果算的是"好天气",你拿着雨天数据来对比,必然对不上。
处理方式是提前约定对比工况。在设计阶段就把好天气和雨天两种工况都算出来,现场实测时阴雨天测噪声和干扰,晴天测电场和磁场,各取所需。没有人会在暴雨天测电场,但偏偏有人会在雨刚停地还湿漉漉的时候去测电场,拿回来对比晴天仿真,那当然对不上。
5.2 弧垂状态没有取对
导线对地距离不是固定的,它随温度、负荷和覆冰状态变化。同一个档距,夏天高温大负荷时弧垂最大,导线最低点离地最近,地面电场最大;冬天低温轻载时导线绷得更紧,电场变小。
现场实测时如果恰逢夏季用电高峰,那么与"年均温度"工况下的仿真值相比,最大场强差出20%以上是正常的。正确做法是实测当天记录导线弧垂状态,或者实测导线最低点对地距离,把仿真模型里对应改掉再对比,而不是拿一个固定设计工况去硬碰。
5.3 线下物体造成的局部畸变
仿真模型很多时候是"干净"的——平坦地面,没有树,没有房子,没有围墙。但现场实测区域往往存在杂物。一辆停在边相下面的金属厢式货车,就能让附近的电场测量值产生明显变化。树木、低矮建筑物、甚至潮湿的灌木丛都会改变局部电场分布。
现场实测遇到这种情况,要尽量选择周围净空的测点,避开明显物体。如果无法避开,那么仿真时就应该把该物体简化建模放进去。这不是降低仿真精度的问题,而是根本不在同一个物理场景里。
5.4 地表条件和土壤参数差异
仿真默认大地是无限大导体平面,土壤电阻率无穷小。实测场地可能是干燥砂石地、水泥路面、或者湿草地,地表条件不同,地面附近电场的垂直分量测量结果会有差异。水泥地面和湿草地相比,在相同线路条件下的场强读数可能差10%上下。
土壤参数对磁场几乎没影响,对电场的影响在中性点有效接地系统下也不大,但在地下有金属管线、接地网密集的区域,地面电场的畸变不能忽视。实测选点时要尽量远离这些埋地金属体。
5.5 运行电压和负荷电流不是设计值
设计阶段算电场用的是最高运行电压,现场实测时线路可能运行在90%甚至更低的电压水平。电场与电压成正比,这一条比较好修正,记录实测时刻的母线电压即可。但磁场是完全跟着负荷电流走的,现场轻负荷时测出来的磁场只有设计满载工况的零头,拿这个数据说仿真不准,是很冤枉的。
解决方法是建立对比前先把仿真模型的电压、电流输入改成实测时刻的运行值,而不是拿着设计工况的报告去对照。这个道理很多实测人员不太了解,需要仿真工程师主动提醒。
5.6 测量仪器的频率响应与标定问题
工频电场测量仪的表头朝向、探头离地高度、仪器的频率响应范围都要检查。有的仪器从0Hz到100kHz都能响应,有的只对标50Hz做了校准,读数会有细微差异。无线电干扰测量对天线位置和极化方向敏感,天线放偏一点,读数就可能差几dB。
所以实测对比仿真,仪器分辨率、校准证书和测量方法都要作为记录项写进报告。排除掉仪器问题之后,再谈模型误差才有意义。
6. 自编脚本时我常用的一组模型参数和边界处理细节
6.1 典型输入参数表
分享一组我常用作验证算例的输入参数,来自一个典型500kV单回水平排列线路,计算电场和磁场都非常顺手。
| 参数 | 取值 |
|---|---|
| 系统电压 | 500kV(计算取550kV最高运行电压) |
| 导线型号 | JL/G1A-400/35 |
| 导线计算直径 | 26.8mm |
| 每相分裂数 | 4 |
| 分裂间距 | 450mm |
| 相间距离 | 12m |
| 排列方式 | 水平排列 |
| 边相离地高度(弧垂最低点) | 15m |
| 计算场点高度 | 1.5m |
| 典型电流 | 1000A(单相) |
基于这组参数,二维模拟电荷法算出来的地面1.5米高度最大场强大约在7到8kV/m之间,边相投影位置附近的场强明显高于线路中心正下方。磁场在电流1000A时,地面附近的磁感应强度大约在15到25μT量级。你可以拿这个量级范围来复核自己的脚本。
6.2 镜像与截断边界处理
大地的镜像处理是模拟电荷法的核心边界条件,这里有一个容易忽略的细节:镜像电荷的电位系数公式里,距离项用的是镜像电荷到场点的距离。导线离地高度相对于导线半径来说非常大,所以自项距离约等于2倍导线高度,这个近似在常规高度下足够精确。但如果遇到特高压线路那种导线离地三四十米的情况,仍然可以这样处理。
计算域的水平范围,我的习惯是从线路中心向两侧各取80米,电场在这个距离上通常已经衰减到背景噪声量级,如果算到80米场强还超过1kV/m,那说明线路参数或者计算设定有问题,需要回头检查。步长方面,导线投影附近加密到0.5米,再往两边按1米步长扩展,既能保证最大场强的捕捉精度,又不至于让计算量失控。
6.3 网格和收敛判据
二维模拟电荷法用得比较多的是直接求解,不涉及迭代,所以收敛判据主要体现在电荷数量上。每根子导线一般放一个线电荷就够,但如果要做更精细的表面场强分布(比如判断子导线表面哪一侧场强最大),就需要在子导线表面放多个匹配点和对应的电荷。
匹配点数量增加后要验证一个东西:电位还原误差。解法是把解出来的电荷代回各匹配点,看计算电位与给定电压的偏差,一般要求相对误差小于1%。如果达不到,不是增加电荷数量就是检查有没有把分裂导线的等效半径用错。这个验证步骤虽然简单,却是脚本正确性的第一道保险。
6.4 一个验证脚本正确性的土办法
最后分享一个我每次换环境、换语言重写脚本时都要做一遍的土办法:用单根无限长直导线做解析解对比。
单根导线在无限大平地上方,地面某一点的电场强度是可以用解析公式直接算出来的。这时候把脚本里的分裂数改成1,导线条数改成1,算出来的场强跟解析公式对比,能对上,说明核心求解器没问题;对不上,那问题一定出在电荷设置或者镜像系数上。通过这个测试之后,再改回4分裂,结果就是可信的。
这个方法我到现在都在用,它比任何复杂的验证算例都管用。做电磁场仿真,先确保最简情况正确,再逐步加复杂条件,是电气工程师最朴素的调试哲学。