潘啟軍 孟慶云 張向明 薛高飛 張 磊

基于互感參數提取的多導體與多導線模型縮減方法研究

潘啟軍 孟慶云 張向明 薛高飛 張 磊

（海軍工程大學艦船綜合電力技術國防科技重點實驗室 武漢 430033）

金屬導體地電流通過地回路公共阻抗在不同電氣電子設備之間形成電磁干擾，嚴重時甚至使敏感設備無法正常工作。然而，地回路公共阻抗在不同回路中表現出不同的阻抗值，因此必須對鄰近效應情況下的金屬導體進行深入研究。本文利用 PEEC方法對鄰近效應時的金屬導體進行了建模，得到了系統矩陣方程。然后，將簡化等效電路和原始電路中導體與導線電流分別進行復頻率 s冪級數展開匹配，得到了電阻、電感及互感參數，給出了同時考慮趨膚效應與鄰近效應時獨立于頻率之外的簡化等效電路模型，并通過算例驗證了參數提取方法的正確性。

電磁兼容 矩量匹配方法 鄰近效應 互感參數提取 簡化等效電路

1 引言

隨著電力電子設備在艦船獨立系統中的廣泛應用，強電系統通過殼體鋼板共地耦合和磁場耦合對弱電系統形成電磁干擾（EMI），使得艦船獨立系統中電磁兼容（EMC）問題日趨嚴重[1-3]。因此，需要對金屬殼體地電流測試和金屬殼體部分單元等效電路（PEEC）建模兩大難點進行全面、深入的研究，形成分析大面積連續(xù)導體必需的實驗測試手段和較完整的理論體系，為全系統仿真、EMI分析等奠定較好的基礎。文獻[4-6]利用鐵心線圈對鋼板地電流測量進行了研究，然而這種測量裝置用于低頻測量時體積較大、鐵心線圈可能會改變鋼板局部地電流的分布以及在任何頻域范圍內均不能進行時域分析。為適應艦船鋼板表面地電流干擾研究的需要，文獻[7-9]研制了可以測量鋼板表面磁場分布的新型磁阻磁場探頭，利用鋼板表面切向磁場的分布來間接測量鋼板中的電流分布，克服了鐵心線圈測量裝置的缺陷。

對金屬導體建模是電磁兼容領域的一大難點，文獻[10-12]利用解析方法對導體阻抗進行了計算和實驗，但解析方法不能給出獨立于頻率之外的系統模型參數。Shizhong Mei在文獻[13-14]中通過漸近波形估算矩量顯式計算方法（Asymptotic Waveform Evaluation, AWE）得到了獨立于頻率之外金屬導體的縮減模型電阻與電感參數，但 AWE等直接矩量顯式計算會出現病態(tài)，當模型較大時會導致模型縮減的支路數有限[15]。為了避免矩量顯示計算所出現的病態(tài)問題，文獻[16-17]提出了基于Krylov子空間arnoldi和 Lanczos算法的模型縮減方法，將直接矩量顯式計算轉化為間接矩量隱式計算，將主導極點的抽取轉化為相對應的主導特征值的抽取，避免了AWE等直接矩量顯式計算方法只能提取有限幾條簡化支路的缺陷，極大提高了模型簡化技術的有效性和普適性。然而，上述文獻僅對孤立導體進行了模型縮減方法的嘗試，國內外對多導體與多導線情況下模型縮減方法的深入研究未見報道。此時，對金屬導體進行建模時不僅需要考慮趨膚效應，而且需要考慮鄰近效應，需給出獨立于頻率之外的等效電路模型，便于和電路仿真軟件兼容。在忽略鄰近效應的情況下，可以將金屬導體視為孤立元件進行處理，但當外部載流導線或導體靠近金屬導體時，將會改變金屬導體內部的電磁場分布，因此必須對鄰近效應時導體的模型進行分析。本文采用 PEEC方法對鄰近效應時的金屬導體進行了建模，將復雜的電磁場計算轉化為電路計算，在對孤立導體進行模型縮減的基礎上，基于互感參數提取得到了多導線與多導體等效縮減電路通用模型。

2 單導體與多導線時的PEEC建模及模型縮減

在矩形截面導體長度方向放置m根導線，其長度和矩形截面導體的長度一樣，并且導線兩端與矩形截面導體平齊如圖1所示（為簡明起見，圖中沒有標出多根導線的激勵電壓），圖 2所示為其對應的等效耦合電路圖。令矩形截面導體的驅動電壓為va，總電流為 ia；鄰近導線中電流為

其兩端電壓為當導體與導線相距很近時，鄰近效應的作用顯著，每個導體用自身的解耦簡化等效電路來描述趨膚效應，同時使用導體各支路與導線間的互感來描述鄰近效應。

圖1 鄰近效應時矩形截面導體Fig.1 A rectangular conductor with proximity effect

圖2 鄰近效應時等效耦合電路Fig.2 Equivalent coupling circuit with proximity effect

根據歐姆定律可得：式中，Mabji= Mabij（i=1, 2, …, m; j=1, 2, …, N）為矩形截面導體與鄰近導線的互感矩陣，ia1、ia2、…、iaN為矩形導體支路的電流。

將式（1）經過變換后用矩陣表示為

式中，EN和Em分別為N維和m維單位陣，且有

導體和導線中的總電流可分別表示為

式中，1N為N維全1列向量，0m為m維全0列向量，1m為m維全1列向量，0N為N維全0列向量。

將圖 2所示的等效耦合電路進行簡化，如圖 3所示。比較圖 3和圖 2可知，圖 3左半部分變?yōu)闆]有鄰近效應時的簡化等效支路，解耦電阻與電感可由 Arnoldi算法或 Lanczos算法計算得到[16,17]，簡化等效支路數為n條，n＜＜N。相應的，圖3中左半部分與右半部分之間的互感矩陣也和圖2中不同。

圖3 鄰近效應時的簡化電路Fig.3 The reduced circuit of the system with proximity

根據圖3所示可描述趨膚與鄰近效應的簡化r-l等效電路，和式（3）、式（4）類似，簡化電路的總電流為

式中，1n為n維全1列向量，式中En為n維單位陣，且有

要準確捕獲趨膚與鄰近效應，式（5）應與式（3）相等，式（6）應與式（4）相等，但這一要求太過苛刻，實際上只需令各對應式的 s冪級數的有限項相等即可。由于電壓的取值任意，將上述關系進一步簡化，因此只需令

和

展開成s的冪級數后的對應項相匹配，即可解出式（7）中互感矩陣M所有互感的值。

經過矩陣變換，式（7）可展開成s的麥克勞林級數為

同理，式（8）可展開成s的麥克勞林級數為

由于式（7）中矩陣M的維數為n×m維，則未知互感參數有nm個，因此只需令式（9）和式（10）中的si的nm個系數相等得到nm個方程，就可以求得nm個互感參數，其中可取i=1,2,3,…,nm。

令式（10）中si的系數為yi，則有

由式（12）可得

將式（13）代入式（11）可得

設

則式（14）可寫為

則由式（17）可得

將式（18）中矩陣M的各互感參數提取出來可得

取i=1,2,3,…,nm，可得nm個方程。

式（13）、式（14）、式（16）及式（20）中使用了廣義逆矩陣，式（20）所求解為極小最小二乘解。

3 多導體情況

對于多導體情況，可采用與2類似的方法求取互感參數，但互感參數求取過程更為復雜。設導體總數為m個，每個導體的剖分數分別為N1,N2,…,Nm，每個導體的解耦簡化等效支路數分別為n1,n2,…,nm。

設

其中 R1,R2,…,Rm分別為各導體解耦簡化電路的對角電阻矩陣，其維數分別為L1,L2,…,Lm分別為各導體解耦簡化電路的對角電感矩陣，其維數分別為；Mkj為第k個導體簡化等效支路與第 j個導體簡化等效支路間的互感矩陣，其維數為。

經過矩陣變換可以得到

其中

將式（23）寫成矩陣形式為

其中 Mequ為所有互感參數組成的維列向量，即所需求取的列向量；H為維矩陣；Z為由原始矩陣計算得到的維列向量，求解方法和第2部分類似。

對于m個導體和n根導線情況，將n根導線的集合看成不解耦的第（m+1）個導體，則此時互感的求取方法和多個導體情況基本相同，只需將 A0對角線上的第（m+1）個分塊矩陣略做修改即可。對于磁性導體，首先根據磁性部分元等效電路（MagPEEC）建立系統方程，然后按照第2和第3部分介紹的方法來研究。

4 鄰近效應時互感算例

4.1 單導體與雙導線

金屬導體長20μm，寬10μm，厚1μm，兩鄰近導線位于寬度正中且平行于長度方向，線徑為0.1μm，兩導線與導體的距離為分別為 0.1μm 和1μm，導線圓截面的幾何平均距離為半徑的 0.778 8倍，頻率范圍為 100MHz～30GHz。表 1為不同支路數的互感值及對應的展開系數，表中互感矩陣的行對應于導體簡化支路，列對應于第1和第2根導線。由模型縮減方法得到的導體簡化等效電路各支路的電阻與電感見表 2。

表2 導體簡化電路各支路的電阻與電感Tab.2 Resistance and inductance values in the reduced RL circuits of a conductor

根據表1可知，簡化電路與原電路的s展開系數誤差隨冪次數的增加而增加，當冪次數小于 10時，兩者的相對誤差小于10%，當冪次數為12時，最大相對誤差約為22%。

4.2 雙導體

兩金屬導體均長 20μm，厚 1μm，第 1個金屬導體寬5μm，第2個金屬導體寬7μm，兩金屬導體在寬度方向距離0.5μm，在厚度與長度方向上對齊，頻率范圍為 100MHz～30GHz。表 3為不同支路數的互感值及對應的展開系數，表4為由模型縮減方法得到的導體簡化等效電路各支路的電阻與電感，表中互感矩陣的行對應于第1個導體，列對應于第2個導體。

表3 不同支路數時的互感值及對應的展開系數Tab.3 Mutual inductance values between the reduced RL circuits and the wire and the corresponding expanded coefficients

表4 導體簡化電路各支路的電阻與電感Tab.4 Resistance and inductance values in the reduced RL circuits of a conductor

根據表3可知，簡化電路與原電路的s展開系數隨冪次數的增加而增加，當冪次數小于20時，兩者的相對誤差小于 10%，當冪次數為 32時，最大相對誤差約為 15.5%。當冪次數繼續(xù)增加時，由于計算得到的系數太小，因此全為 0。由此可見，使用本文介紹的方法，s冪次數不能取得過高。為降低s冪次數的值，可只考慮導體主導簡化支路之間、導體主導簡化支路和導線之間和互感，降低所求解方程的維數，也降低了所需s冪次數的值。另外，可基于矩陣變換的概念來求解互感，從而避免求解s冪級數系數的值。

5 結論

在孤立導體PEEC建?；A上，本文建立了鄰近效應時金屬導體的PEEC模型。首先利用模型縮減技術得到金屬導體簡化等效電路的電阻及電感參數后，然后提取出簡化等效電路各支路和鄰近導線或導體之間的互感參數，并通過算例驗證了該縮減方法的正確性。本文建立的簡化等效模型和計算得到的電阻、電感及互感參數獨立于頻率之外，為全系統仿真、EMI分析等奠定了較好的基礎。

[1] 趙治華, 馬偉明. 鋼板地阻抗的頻率特性[J]. 中國電機工程學報, 2003, 23(8): 119-124.

Zhao Zhihua, Ma Weiming. Frequency characteristic of ground impedance for a flat steel plane[J]. Proceeding of CSEE, 2003, 23(8): 119-124.

[2] Zhao Zhihua, Yuan Jiansheng, Ma Weiming. Coupling model and impedance calculation of the steel ground loops with proximity effect. IEEE Transactions on Electromagnetic Compatibility, 2006, 48(3): 522-529.

[3] 咸哲龍, 鐘玉林, 孫旭東, 等. 用于傳導電磁干擾分析的接地回路模型與參數[J]. 中國電機工程學報, 2005, 25(7): 156-161.

Xian Zhelong, Zhong Yulin, Sun Xudong et al. Grounding circuit model and parameters for conducted EMI analysis[J]. Proceedings of CSEE, 2005, 25(7): 156-161.

[4] Tang Shiping, Jiang Quanxin. Calibration method and device of surface current probe[C]. International Sysposium on Electromagnetic Compatibility Preceedings, 1997: 130-133.

[5] 蔣全興, 湯仕平. 表面電流探頭的校準和應用[J].電子技術應用, 1997(11): 35-37.

Jiang Quanxin, Tang Shiping. Calibration and application of surface current probe[J]. Application of Electronic Technique, 1997(11): 35-37.

[6] 宋文武, 徐曉明, 張坤元, 等. 地電流測試裝置研究[J]. 電磁兼容技術, 2005, 2: 46-49.

Song Wenwu, Xu Xiaoming, Zhang Kunyuan et al. A study of ground current testing apparatus[J]. Electromagnetic Compatibility Technology, 2005, 2: 46-49.

[7] 潘啟軍, 馬偉明, 趙治華, 等. 適用于艦船鋼板地電流檢測的新型磁場測量儀器設計[J]. 電工技術學報, 2004, 19(11): 82-86.

Pan Qijun, Ma Weiming, Zhao Zhihua, et al. A novel magnetic-field measurement apparatus design applied in ground current testing of the steel plate in vessels [J]. Transactions of China Electrotechnical Society, 2004, 19(11): 82-86.

[8] Pan Qijun, Ma Weiming, Zhao Zhihua et al. A modified design of two-stage amplification ground current testing apparatus in large steel plane[C]. ICEMI, 2005, 2: 416-420.

[9] 潘啟軍, 馬偉明, 趙治華, 等. 鋼板地電流新型測量儀器校準系統設計及分析[J]. 電工技術學報, 2006, 21(2): 114-118.

Pan Qijun, Ma Weiming, Zhao Zhihua, et al. An approach to calibrate the novel apparatus applied in ground current testing of steel plane[J]. Transactions of China Electrotechnical Society, 2006, 21(2): 114-118.

[10] Coperich K M, Ruehli A E．Enhanced skin effect for partial-element equivalent-circuit(PEEC) models[J]. IEEE Transactions on Microwave Theory and Techniques, 2000, 48(9): 1435-1442.

[11] Tsuk M J, Kong A J．A hybrid method for the calculation of the resistance and inductance of transmission lines with arbitrary cross sections[J]. IEEE Transactions on Microwave Theory and Techniques, 1991, 39(8): 1338-1347.

[12] Zhao Zhihua, Ma Weiming．AC impedance of an isolated flat conductor[J]．IEEE Transactions on Electromagnetic Compatibility, 2002, 44(3): 482-486.

[13] Mei Shizhong, Ismail Y I. Modeling skin and proximity effects with reduced realizable RL circuits [J]. IEEE Transactions on Very Large Scale Integration(VLSI) System, 2004, 12(4): 437-447.

[14] Mei Shizhong, Ismail Y I．Modeling skin effect with reduced decoupled R-L circuits[C]．The International Symposium on Circuits and Systems, Bangkok, Thailand, 2003.

[15] Pillage L T, Rohrer R A. Asymptotic waveform evaluation for timing analysis[J]. IEEE Transactions Computer-Aided Design, 1990, 9(4): 352-366.

[16] 潘啟軍, 馬偉明, 趙治華, 等. 基于Arnoldi算法的金屬導體 PEEC 模型縮減方法[J]. 電工技術學報, 2008, 23(3): 1-7.

Pan Qijun, Ma Weiming, Zhao Zhihua, et al. PEEC model reduction of a conductor based on Arnoldi algorithm[J]. Transactions of China Electrotechnical Society, 2008, 23(3): 1-7.

[17] 潘啟軍, 馬偉明, 趙治華, 等. 基于 Krylov子空間中 Lanczos算法的金屬導體PEEC模型縮減[J]. 中國電機工程學報, 2008, 28(24): 130-135.

Pan Qijun, Ma Weiming, Zhao Zhihua, et al. PEEC model reduction of a conductor based on Lanczos algorithm in Krylov subspace[J]. Proceedings of the CSEE, 2008, 28(24): 130-135.

A Model Reduction Study of Multiple Conductors and Multiple Wires Based on Extraction of Mutual Inductance Parameters

Pan Qijun Meng Qingyun Zhang Xiangming Xue Gaofei Zhang Lei

（National Key Laboratory for Vessel Integrated Power System Technology Naval University of Engineering Wuhan 430033 China）

The Electromagnetic Interference(EMI) among various electrical and electronic equipment can be formed by ground current in metal plate via common ground loop impedance. It will even result in the malfunction of the sensitive devices under some serious conditions. However, the common impedance has different value in different loops. Therefore, it is necessary to analyze and study a conductor with proximity. This paper adopts partial element equivalent circuit(PEEC) to model a conductor with proximity, and to obtain the system matrix equation. After matching the corresponding current coefficients expanded in powers of s between the reduced circuit and the original circuit, the resistance, inductance, and mutual inductances of the reduced circuit are provided. And then, the reduced equivalent circuit model regarding of skin and proximity effects is presented independent of frequency. Furthermore, the mutual inductance between the reduced branches and the proximity wires are obtained, which is validated by examples.

Electromagnetic compatibility(EMC), moment matching method, proximity effect, extraction of mutual inductance parameters, reduced equivalent circuit

TM53

潘啟軍 男，1972年生，博士，副教授，研究方向為電力電子應用與電磁兼容性研究。

國家重點基礎研究發(fā)展計劃（2013CB035601）和國家自然科學基金（51077129）資助項目。

2013-06-25 改稿日期 2013-08-26

孟慶云 男，1979年生，講師，研究方向為電力電子大容量電能變換技術。