5 采用矩量法计算输电线路杆塔再辐射场强(讲稿).docVIP

第5篇 采用矩量法计算输电线路杆塔再辐射场强 清华大学电机系 邹 军，袁建生 求解细长圆柱导体端口电压的积分方程 在已知外加场的作用下，用迟后位的积分来表示和产生的散射场，则在导体表面有如下边界条件。当待求的导体为有源的发射天线时，，则边界条件即为导体表面处：。 另外，对于任意的圆导线、杆件等细长导体，如直径－长度比0.01，可作如下的近似[]： 电流只沿导体的轴向流动； 导体的电流和电荷密度近似为，位于轴线上的线电流I和线电荷线密度σ。 只对导线表面上的轴向切向分量使用边界条件，如合成电场切向分量为0 (忽略表面阻抗时)。 求解空间中的电场强度的公式为： 式~中： 为导体表面轴向上的导体产生的电场强度 为导体表面轴向上的矢量磁位，为空间中的矢量磁位 为导体表面处的电位标量 为空间磁导率 为轴线上的源点到导体表面的场点的距离 为导体上的线电流 为导体上的线电荷密度 称为相位常数，单位1/m，为电源频率 为介电常数 为角频率 2. 矩量法原理简介 矩量法（Method Of Moments，简写为MOM）是求解算符方程的一种普遍方法。20世纪70年代，R.F.Harrington将矩量法全面引入对电磁场问题的求解中。假设有非齐次方程 式中为线性算符，为待求的场函数，为已知源函数。在函数空间中，将展开为一组完备基函数之和，即令 式中，是系数，被称为展开函数或者基函数。当式取无穷项之和时，可获得精确解。当取有限项时，则为近似解，代入方程，得 规定一个适当的内积，在的值域上定义n个权函数或称检验函数，，，…，。根据算符方程等号两边同时做内积相等的原则，对每个取式两边的内积，形成矩阵方程 写成矩阵形式为 式中 ，， 求解矩阵方程，得到待定系数（），由此即可得到近似解。 求解一个特定问题时，采用矩量法的关键在于选择适当的基函数和权函数。通常，在实际问题中，求矩阵元素中的积分是很困难的。因此，选择基函数时，除可采用全域基函数外，一种更实用的近似方法是分域基函数。在分域基法中，每个基函数只在定义域的分域上存在。这种方法可以简化矩阵元素的计算。 3. 矩量法求解A-φ积分方程 式~式的积分方程组描述了线天线的电流与表面电场的关系，在已知边界条件的情况下，求解式~式即可得到天线的电流和电荷密度分布。采用矩量法求解式~式的步骤如下： 3.1选取基函数和匹配函数 如图1所示，对线天线进行离散化，将导体沿轴向剖分单元，设置节点。在单元内分别取电流和电荷的一组基函数，即用分段插值的方法近似表示电流和电荷的分布。本文对电流分布采用线性插值，对电荷分布采用恒值插值，所以电流的基函数沿轴向为三角波，而线电荷密度沿轴向为矩形波。电流即表示为： 其中，基函数为、定义在区间的线性函数（也称为形状函数），此时，待求量即为节点电流。通过式，电荷密度可由电流表示 因此方程的待求量仅为节点电流。 图1 线天线的剖分和插值示意图 将、分别带入、得到 完成单元剖分和基函数选取之后，选取匹配线段和检验函数。如图2所示，设置匹配线段于导线表面，对匹配线段进行剖分，将匹配线段的起始点选在相邻单元的中点，匹配函数选为。 图2 匹配线段的示意图 3.2形成矩阵方程 式两端同时与匹配函数做内积即得到 可进一步整理得到 式中，为第k段匹配线段，为第k段匹配线段两端的电压差，为第k段匹配线段两端的标量电位差。将各节点电流和匹配线段两端的电压差整理成矩阵形式，定义矩阵 将式、带入式，经过整理之后，式可写成矩阵形式 至此，式~的积分方程已变换为一组联立的代数方程。其中，为待求的未知量。由式给出的边界条件可知： 即 求出矩阵的元素即可解出，矩阵的元素中包含了积分： 、、 以上积分的求解须借助数值积分的方法，本文采用高斯积分的方法来求解。所谓高斯积分法，就是在单元内选择某些积分点（称为高斯采样点），求出积分点处的被积函数值，然后用一些权因子乘这些函数值，再求总和，便可得到积分值[]。一般每个单元设置3个高斯点，将积分区间归一化为[0，1]后，3个高斯采样点的高斯积分公式为： 求解矩阵的部分元素时，源点单元与匹配线段产生重合，于是在式、中，产生奇异积分。遇到奇异积分时，应增加高斯点数量，例如将高斯点采样数量增至30个，高斯点和全因子的求解可参考文献[]。 3.3矩阵方程的修正 在解出矩阵的全部元素后，还需根据电流注入情况和导体的联结情况等，对线性方程进行必要的修正。 假设注入电流为，注入点位于导体第k个节点处。此时，第k节点处的电流即为注入电流，即 为保证求解方程时，式成立，令矩阵中的元素 同时