孫 科，江全元

(浙江大學(xué)電氣工程學(xué)院，浙江杭州310027)

0 引言

近年來(lái)，針對(duì)電氣化鐵路牽引供電系統(tǒng)的仿真分析研究也成為重要研究?jī)?nèi)容［1］。傳統(tǒng)的仿真分析研究主要集中在牽引網(wǎng)的潮流計(jì)算分析［2-5］，而目前針對(duì)牽引網(wǎng)諧波、短路以及操作過(guò)電壓等暫態(tài)過(guò)程需要的電磁暫態(tài)仿真研究主要都是通過(guò)現(xiàn)有電磁仿真工具實(shí)現(xiàn)的［6-7］。這種研究的局限性在于:①軟件仿真模塊基本針對(duì)電力系統(tǒng)設(shè)計(jì)，對(duì)于牽引網(wǎng)特性支持不完全;②仿真軟件針對(duì)大規(guī)模長(zhǎng)線路的牽引網(wǎng)詳細(xì)仿真模型建立步驟繁瑣，不能方便調(diào)節(jié)仿真參數(shù);③仿真軟件算法封裝，程序二次開(kāi)發(fā)能力差，無(wú)法針對(duì)牽引網(wǎng)進(jìn)行優(yōu)化。

本研究從牽引網(wǎng)各電氣設(shè)備基本模型出發(fā)，建立一套針對(duì)牽引網(wǎng)的電磁暫態(tài)仿真模型，并針對(duì)牽引網(wǎng)結(jié)構(gòu)進(jìn)行算法設(shè)計(jì)，簡(jiǎn)化建模步驟，為算法優(yōu)化提供接口。

1 牽引網(wǎng)電磁暫態(tài)模型概述

1.1 牽引網(wǎng)供電方式

電氣化鐵路牽引網(wǎng)通常的供電方式包括有:直接供電方式(TR)、帶回流線的直接供電方式(TRNF)、吸流變壓器供電方式(BT)、自耦變壓器供電方式(AT)、同軸電力電纜供電方式(CC)、混合供電方式(直接供電方式+AT 供電方式)等［8-10］，其中在我國(guó)高速鐵路中主要采用了TRNF 供電方式和AT 供電方式。在進(jìn)行牽引網(wǎng)的電磁暫態(tài)建模時(shí)，傳統(tǒng)研究方式需要按照牽引網(wǎng)所采用的供電接線方式來(lái)建立相關(guān)模型。

1.2 牽引網(wǎng)統(tǒng)一電磁暫態(tài)模型

上述所提及的牽引網(wǎng)供電方式，從整體上看，都可以認(rèn)為牽引網(wǎng)的主線路是平行多導(dǎo)體傳輸線，而其他電氣設(shè)備則可以看作并聯(lián)元件和串聯(lián)元件分布在線路中。其中，主要電氣設(shè)備劃分歸類如表1所示。

表1 牽引網(wǎng)主要電氣設(shè)備類型劃分

在進(jìn)行牽引網(wǎng)電磁暫態(tài)建模時(shí)，整個(gè)網(wǎng)絡(luò)模型都等效為電阻網(wǎng)絡(luò)和電流源組合的形式。本研究針對(duì)串聯(lián)元件和并聯(lián)元件，先忽略外部影響，得到獨(dú)立的元件模型。在進(jìn)行仿真計(jì)算時(shí)，筆者在系統(tǒng)相應(yīng)節(jié)點(diǎn)的導(dǎo)納矩陣和節(jié)點(diǎn)注入電流源處添加各元件對(duì)應(yīng)元素，再聯(lián)立求解，即可得全系統(tǒng)的解。假設(shè)整個(gè)牽引網(wǎng)有n條導(dǎo)線，按照牽引網(wǎng)串聯(lián)元件與并聯(lián)元件空間距離分布將牽引網(wǎng)導(dǎo)線以元件為切面分割成K 段，其中切面k 和k+1 間的線路模型示意圖如圖1所示。

圖1 牽引網(wǎng)統(tǒng)一模型示意圖

并聯(lián)元件和串聯(lián)元件模型的節(jié)點(diǎn)電壓方程通用表達(dá)式為:

式中:G—并聯(lián)元件的導(dǎo)納矩陣，n ×n 的方陣;J—等效節(jié)點(diǎn)注入電流源，n 維向量;串聯(lián)元件連接了2n 個(gè)節(jié)點(diǎn)，形成的導(dǎo)納矩陣G 為2n ×2n 的方陣，等效節(jié)點(diǎn)注入電流源J 為2n 維向量。

本研究采用上述建模方式對(duì)牽引網(wǎng)進(jìn)行建模，相比于利用現(xiàn)有的電磁暫態(tài)仿真工具，其優(yōu)點(diǎn)包括:①采用統(tǒng)一的建模方式，適用于不同牽引供電方式的牽引網(wǎng);②在處理AT 供電方式、CC 供電方式等需要對(duì)牽引網(wǎng)單一供電臂導(dǎo)線進(jìn)行多段分割的情形時(shí)，能夠便捷實(shí)現(xiàn);③在仿真過(guò)程中，能夠在機(jī)車實(shí)際運(yùn)行位置處動(dòng)態(tài)加入等效機(jī)車注入電流源，而不用修改仿真工具中的模型。

2 牽引網(wǎng)電氣元件電磁暫態(tài)模型

在牽引網(wǎng)的電氣元件電磁暫態(tài)模型中，特別考慮了牽引網(wǎng)導(dǎo)線、牽引變壓器、自耦變壓器的模型，其他電氣設(shè)備的模型可以參考穩(wěn)態(tài)模型類似推導(dǎo)。

2.1 牽引網(wǎng)導(dǎo)線電磁暫態(tài)模型

牽引網(wǎng)導(dǎo)線模型是牽引網(wǎng)模型研究的基礎(chǔ)，其采用的電磁暫態(tài)仿真模型參考了電力系統(tǒng)輸電線模型［11-15］。本研究根據(jù)多導(dǎo)體線路波過(guò)程建立基本波動(dòng)方程，將所得的多相耦合線路模型通過(guò)相模轉(zhuǎn)換解耦，等效為在模坐標(biāo)內(nèi)相互獨(dú)立的傳輸線，并進(jìn)行電磁暫態(tài)計(jì)算公式推導(dǎo)，最后在相坐標(biāo)下建立節(jié)點(diǎn)電壓方程。

2.1.1 單根輸電線電磁暫態(tài)模型

以送端s 端為例，建立的電磁暫態(tài)模型為:

式中:L，C—線路單位長(zhǎng)度的電感和電容;x—與送端的距離;Z—波阻抗，τ—線路上波傳播時(shí)間，l—線路長(zhǎng)度。

2.1.2 牽引網(wǎng)導(dǎo)線相模變換

利用相模變換矩陣將方程由相坐標(biāo)轉(zhuǎn)換到模坐標(biāo)進(jìn)行求解，再進(jìn)行后續(xù)計(jì)算。

設(shè)相模變換矩陣為Q，獲得相模變換為:

式中:Um，Im——模坐標(biāo)下電壓、電流;

可以參考單根輸電線電磁暫態(tài)模型獲得平行多導(dǎo)體導(dǎo)線在模坐標(biāo)下的計(jì)算模型，對(duì)于s 端有:

式中:Ym——模坐標(biāo)下導(dǎo)納矩陣，為各模量波阻抗倒數(shù)組成的對(duì)角陣，Ym=(Zm)－1;

2.1.3 牽引網(wǎng)導(dǎo)線節(jié)點(diǎn)電壓方程

對(duì)整個(gè)牽引供電系統(tǒng)進(jìn)行電磁暫態(tài)仿真時(shí)，需要將得到的模坐標(biāo)下的多導(dǎo)線模型轉(zhuǎn)換到相坐標(biāo)下，進(jìn)而建立節(jié)點(diǎn)電壓方程求解。

仍以s 端為例，根據(jù)相模轉(zhuǎn)換公式，模坐標(biāo)下的公式轉(zhuǎn)換到相坐標(biāo)后為:

式中:Y—相坐標(biāo)下導(dǎo)納矩陣，Y=QYmQT。

對(duì)于電氣化鐵路牽引網(wǎng)平行多導(dǎo)體導(dǎo)線，由于其在每隔一段距離就會(huì)存在線路間的短接或者接地，需要按照實(shí)際線路連接步長(zhǎng)將牽引網(wǎng)導(dǎo)線劃分為多段進(jìn)行處理。針對(duì)傳輸距離并不長(zhǎng)的導(dǎo)線，研究者可以用集中在線路兩端的電阻網(wǎng)等效模擬沿線路均勻分布的電阻，得到集中等效到線路兩端的電阻網(wǎng)絡(luò)矩陣Rl:

仍以s 端為例，假設(shè)有n 條牽引網(wǎng)導(dǎo)線，則可以得到牽引網(wǎng)導(dǎo)線的節(jié)點(diǎn)電壓方程為:

2.2 牽引變壓器電磁暫態(tài)模型

通常情況下，電氣化鐵路采用的是非對(duì)稱式接線方法，將上、下行的牽引供電網(wǎng)接入電力系統(tǒng)三相供電網(wǎng)。牽引變壓器最常用的接線方式是三相V/v 接線，廣泛應(yīng)用于京津、京滬等高鐵線路中。

2.2.1 單相變壓器電磁暫態(tài)模型

針對(duì)一阻抗支路建立電磁暫態(tài)模型，這是組成變壓器勵(lì)磁支路和等效漏電感支路的模型基礎(chǔ)。利用隱式梯形積分法得到離散化計(jì)算方程為:

式中:in+1，un+1—n +1 時(shí)刻的支路電流及支路電壓;G—支路等效導(dǎo)納，G=h/(2L +hR);Jn—n 時(shí)刻歷史注入電流，Jn=Gun+(2L－h(huán)R)in/(2L +hR);h—計(jì)算步長(zhǎng)。

對(duì)于單項(xiàng)變壓器，通常使用等效電路進(jìn)行建模，將變壓器一次側(cè)和二次側(cè)的漏電阻和漏電抗都等效集中到二次側(cè)，將變壓器的勵(lì)磁電感等效并聯(lián)在一次側(cè)。

本研究推導(dǎo)其電磁暫態(tài)模型主要采用受控源模型來(lái)等效理想變壓器［16-18］，可以獲得其等效電路圖如圖2所示。

圖2 基于受控源的單相雙繞組變壓器模型

可以得到變壓器的矩陣模型方程為:

2.2.2 三相V/v 接線變壓器電磁暫態(tài)模型

三相V/v 接線中一次側(cè)中有一相進(jìn)線與另一相出線首尾相連，采用和單相變壓器相同處理方法，可以獲得三相V/v 接線變壓器基于受控源的等效電路圖如圖3所示。

圖3 基于受控源的三相V/v 接線變壓器等效電路

由此可以得到變壓器的矩陣模型方程。

2.3 自耦變壓器電磁暫態(tài)模型

本研究采用基于受控源的自耦變壓器模型，同時(shí)考慮了自耦變壓器磁耦合和直接電耦合的特點(diǎn)［19-20］。

本研究將單相雙繞組自耦變壓器的公共繞組以及串聯(lián)繞組的漏阻抗都等效到一次側(cè)，建立的模型等效電路圖如圖4所示。

圖4 基于受控源原理的單相雙繞組自耦變壓器模型

得到變壓器的矩陣模型方程為:

2.4 其他并聯(lián)元件

2.4.1 并聯(lián)阻抗元件

在牽引網(wǎng)計(jì)算中，可以用并聯(lián)阻抗元件描述網(wǎng)絡(luò)中的并聯(lián)補(bǔ)償裝置、無(wú)源濾波裝置等?；镜慕Y(jié)構(gòu)為兩導(dǎo)線間并聯(lián)阻抗或者一導(dǎo)線與地并聯(lián)阻抗。其他并聯(lián)形式都可以分解為基本形式后疊加得到。

當(dāng)兩導(dǎo)線間并聯(lián)阻抗時(shí)，可以得到元件節(jié)點(diǎn)電壓方程為:

當(dāng)一導(dǎo)線與地并聯(lián)阻抗時(shí)，可以得到元件節(jié)點(diǎn)電壓方程為:

2.4.2 橫向連接和短路故障

橫向連接包括了牽引供電方式中線路的金屬性連接，AT 供電方式中的CPW 線、BT 供電方式中的吸上線等都可用該模型描述。短路故障在仿真計(jì)算中也采用和橫向連接類似的方法。

實(shí)質(zhì)上，橫向連接和短路故障的金屬性連接都是用很小的電阻等效，屬于并聯(lián)阻抗元件的特殊情況，其元件節(jié)點(diǎn)電壓方程為:

當(dāng)兩導(dǎo)線間發(fā)生連接時(shí):

當(dāng)一導(dǎo)線與地發(fā)生連接時(shí):

式中:Δ—1 很大的實(shí)數(shù)，可取1012。

2.4.3 電力機(jī)車

在牽引網(wǎng)仿真分析過(guò)程中，機(jī)車一般等效為在接觸線和鋼軌之間并聯(lián)的元件，隨運(yùn)行時(shí)間在系統(tǒng)切面上移動(dòng)。通常工程中采用等效電流源來(lái)模擬機(jī)車，而更精確的計(jì)算需要按照詳細(xì)的建模仿真得出。本研究主要討論牽引網(wǎng)的模型，故機(jī)車選取了電流源模型來(lái)模擬。

2.5 其他串聯(lián)元件

2.5.1 串聯(lián)阻抗元件

串聯(lián)阻抗元件可以用來(lái)模擬線路中串聯(lián)的無(wú)源RLC 元件，形成的節(jié)點(diǎn)導(dǎo)納矩陣為2n ×2n 的方陣，節(jié)點(diǎn)電壓方程為:

2.5.2 斷線故障

斷線故障可以看做阻抗很小的串聯(lián)阻抗元件，等效模型中，等效注入電流源為0，等效導(dǎo)納為很小的實(shí)數(shù)，可以取10－6。

3 牽引網(wǎng)電磁暫態(tài)仿真計(jì)算

3.1 牽引網(wǎng)電磁暫態(tài)節(jié)點(diǎn)電壓方程建立

筆者已經(jīng)介紹了牽引網(wǎng)整體模型的建立方法，現(xiàn)本研究以AT 供電方式為例介紹牽引網(wǎng)電磁暫態(tài)節(jié)點(diǎn)電壓方程建立過(guò)程。

在建模時(shí)，因牽引網(wǎng)中某些導(dǎo)線持續(xù)性并聯(lián)，可簡(jiǎn)化計(jì)算將這些導(dǎo)線合并。AT 供電方式中通常是將左、右平行的鋼軌合并、接觸線和承力索合并。

方程建立首先根據(jù)牽引網(wǎng)電氣元件SS、AT 等的空間位置分布以及CPW、接地裝置間距和機(jī)車位置，分割牽引網(wǎng)導(dǎo)線，使每段導(dǎo)線兩端都連接不同空間位置的牽引網(wǎng)電氣元件。整個(gè)系統(tǒng)計(jì)算節(jié)點(diǎn)由此確定劃分。

然后對(duì)系統(tǒng)節(jié)點(diǎn)電壓方程的導(dǎo)納矩陣和等效注入電流源，按照各電氣元件的獨(dú)立模型，添加對(duì)應(yīng)項(xiàng)，最后建立節(jié)點(diǎn)電壓方程:

式中:Un+1，Un+1—n +1 時(shí)刻的節(jié)點(diǎn)電壓和電源節(jié)點(diǎn)注入電流;Y—系統(tǒng)導(dǎo)納矩陣;Jn—n 時(shí)刻的等效節(jié)點(diǎn)注入電流，由各電氣設(shè)備和牽引網(wǎng)導(dǎo)線模型中等效電流源計(jì)算得出。

3.2 牽引網(wǎng)電磁暫態(tài)計(jì)算方法

牽引網(wǎng)電磁暫態(tài)仿真計(jì)算的步驟如下:

步驟1。根據(jù)牽引網(wǎng)導(dǎo)線空間分布和參數(shù)，計(jì)算其分布參數(shù)，同時(shí)求解相模轉(zhuǎn)換矩陣，獲得轉(zhuǎn)換后的參數(shù)矩陣;

步驟2。根據(jù)牽引網(wǎng)上電氣設(shè)備及接地裝置分布，將牽引網(wǎng)導(dǎo)線分割，根據(jù)牽引網(wǎng)元件模型，建立牽引網(wǎng)全線的節(jié)點(diǎn)電壓方程;

步驟3。利用牽引網(wǎng)潮流計(jì)算程序獲得牽引網(wǎng)各計(jì)算節(jié)點(diǎn)的初值;

步驟4。判斷是否有短路、斷路、開(kāi)關(guān)閉合等改變系統(tǒng)拓?fù)浣Y(jié)構(gòu)的操作。若有，則按此修改系統(tǒng)的導(dǎo)納矩陣和等效注入電流源項(xiàng)。

步驟5。根據(jù)歷史電壓電流數(shù)據(jù)，計(jì)算各節(jié)點(diǎn)處的等效電流源注入電流，包括電力系統(tǒng)電網(wǎng)、牽引網(wǎng)導(dǎo)線、牽引變壓器、自耦變壓器等，同時(shí)根據(jù)機(jī)車當(dāng)前運(yùn)行位置，在對(duì)應(yīng)計(jì)算節(jié)點(diǎn)處添加機(jī)車注入電流。

步驟6。求解當(dāng)前時(shí)刻點(diǎn)節(jié)點(diǎn)電壓方程，記錄計(jì)算結(jié)果，同時(shí)針對(duì)牽引網(wǎng)導(dǎo)線計(jì)算模坐標(biāo)下的電壓電流。

步驟7。重復(fù)步驟4、5、6，其中步驟5 中的歷史電壓電流數(shù)據(jù)采用上一時(shí)步中步驟6 計(jì)算所得結(jié)果。循環(huán)計(jì)算，直至仿真時(shí)刻到達(dá)設(shè)定的仿真結(jié)束時(shí)間為止。

4 算例分析

根據(jù)本研究提出的牽引網(wǎng)電磁暫態(tài)仿真計(jì)算方法，編寫(xiě)了牽引網(wǎng)電磁暫態(tài)計(jì)算程序。采用一AT 供電方式的牽引網(wǎng)進(jìn)行仿真算例計(jì)算，該牽引網(wǎng)導(dǎo)線型號(hào)和空間分布如圖5所示，電氣設(shè)備位置分布如表2所示。牽引變壓器一次側(cè)接入220 kV 電網(wǎng)，二次側(cè)電壓為2 ×27.5 kV，牽引變壓器容量為63 MVA，牽引變壓器短路阻抗8.6%，鋼軌泄漏電阻為10 Ω/km，每隔5 km 鋼軌-保護(hù)線經(jīng)CPW 線連接，同時(shí)保護(hù)線經(jīng)小電阻接地，小電阻阻值1 Ω［21-22］。

圖5 牽引網(wǎng)導(dǎo)線空間分布圖

表2 測(cè)試牽引網(wǎng)變電所位置分布

4.1 沿線電壓分布分析

測(cè)試算例利用等效注入電流源模擬機(jī)車電流注入，機(jī)車運(yùn)行至距起始點(diǎn)27.84 km 處，注入電流為950 A。

牽引網(wǎng)接觸線沿線的電壓有效值圖如圖6所示，包括了采用本研究方法的電磁暫態(tài)計(jì)算程序和PSCAD 仿真所得結(jié)果。筆者分析了采用本研究提供的電磁暫態(tài)計(jì)算方法和采用PSCAD 仿真獲得的沿線接觸線電壓分布的誤差，分析結(jié)果如表3所示。

圖6 牽引網(wǎng)接觸線沿線分布電壓

表3 牽引網(wǎng)接觸線沿線分布電壓誤差比較

結(jié)果表明，在機(jī)車運(yùn)行處牽引網(wǎng)接觸線有明顯的電壓降。采用本研究方法建立的牽引網(wǎng)電磁暫態(tài)仿真計(jì)算與PSCAD 計(jì)算在穩(wěn)態(tài)時(shí)電壓空間分布基本一致，相對(duì)誤差均在0.2%以內(nèi)，說(shuō)明本研究方法能夠準(zhǔn)確反映牽引網(wǎng)運(yùn)行狀況。

4.2 短路故障分析

測(cè)試算例分析了牽引網(wǎng)接觸線發(fā)生接地短路和發(fā)生牽引線負(fù)饋線線間短路時(shí)的短路電流。機(jī)車運(yùn)行到距起始點(diǎn)85.3 km 處，在距離起始點(diǎn)95.3 km 處發(fā)生短路故障。

4.2.1 接觸線接地短路

當(dāng)發(fā)生接觸線接地短路時(shí)，牽引變TPS2 接觸線右臂出口處的短路電流計(jì)算結(jié)果與PSCAD 仿真結(jié)果比較如圖7所示。

圖7 發(fā)生接觸線接地短路時(shí)電流

結(jié)果表明，在發(fā)生牽引網(wǎng)接地短路瞬間電流增大而后趨于平穩(wěn)。本研究通過(guò)對(duì)圖7 的仿真結(jié)果誤差進(jìn)行分析，得到兩種方法的均方根差值為34.17 A，以1 kA為基準(zhǔn)時(shí)相對(duì)誤差3.42%，兩種方法仿真結(jié)果基本一致。基于本研究算法的程序能夠反映牽引網(wǎng)發(fā)生接觸線接地短路時(shí)系統(tǒng)運(yùn)行狀況。

4.2.2 接觸線和負(fù)饋線線間短路

當(dāng)發(fā)生接觸線負(fù)饋線接地短路時(shí)，牽引變TPS2接觸線右臂出口處的短路電流計(jì)算結(jié)果與PSCAD 仿真結(jié)果比較如圖8所示。

圖8 發(fā)生接觸線負(fù)饋線線間短路時(shí)電流

結(jié)果表明，在發(fā)生牽引網(wǎng)接觸線線間短路瞬間電流增大而后趨于平穩(wěn)。通過(guò)對(duì)圖8 的仿真結(jié)果誤差分析，得到兩種方法的均方根差值為33.23 A，以1 kA 為基準(zhǔn)時(shí)相對(duì)誤差3.23%，兩種方法仿真結(jié)果基本一致?；诒狙芯克惴ǖ某绦蚰軌蚍从碃恳W(wǎng)發(fā)生接觸線負(fù)饋線線間短路時(shí)系統(tǒng)運(yùn)行狀況。

4.3 斷路故障分析

測(cè)試算例時(shí)，分區(qū)所開(kāi)關(guān)閉合，牽引變SP1 和SP2同時(shí)為兩牽引變間公共線路供電。機(jī)車運(yùn)行至距起始點(diǎn)65.07 km 處，在距離起始點(diǎn)77.03 km 處發(fā)生接觸線短路故障，原來(lái)由SP1 和SP2 雙牽引變供電的機(jī)車變?yōu)橹挥蒚PS1 單牽引變供電，機(jī)車處的接觸線電壓的計(jì)算結(jié)果和PSCAD 仿真結(jié)果比較如圖9所示。

圖9 發(fā)生斷線故障時(shí)接觸線電壓

結(jié)果表明，當(dāng)機(jī)車由雙牽引變供電因斷線故障變?yōu)橛奢^遠(yuǎn)的牽引變供電時(shí)，機(jī)車所處位置接觸線電壓幅值降低。通過(guò)對(duì)圖9 的仿真結(jié)果誤差分析，得到兩種方法的均方根差值為0.603 kV，以27.5 kV 為基準(zhǔn)時(shí)相對(duì)誤差2.19%，兩種方法仿真結(jié)果基本一致?；诒狙芯克惴ǖ某绦蚰軌蚍从碃恳W(wǎng)斷路故障時(shí)系統(tǒng)運(yùn)行狀況。

4.4 算法優(yōu)勢(shì)

通過(guò)上述算例分析表明，采用本研究算法進(jìn)行牽引網(wǎng)電磁暫態(tài)仿真，能夠準(zhǔn)確有效地反映牽引網(wǎng)電磁暫態(tài)過(guò)程的運(yùn)行狀態(tài)。而相比于傳統(tǒng)地利用PSCAD進(jìn)行仿真計(jì)算，本研究算法的優(yōu)勢(shì)包括:

(1)針對(duì)牽引網(wǎng)開(kāi)發(fā)，建立模型以及修改模型參數(shù)操作方便，適用于多種牽引網(wǎng)絡(luò)結(jié)構(gòu)，特別是針對(duì)分段較多的長(zhǎng)牽引網(wǎng)線路，相比于PSCAD 拖動(dòng)模型模塊搭建過(guò)程，極大地減少了操作量;

(2)能夠支持牽引線導(dǎo)線采用不同參數(shù)的情形，克服了PSCAD 模型上的不足;

(3)算法過(guò)程可控，建立的節(jié)點(diǎn)電壓方程為對(duì)角分塊矩陣，各分塊矩陣間相互解耦，可作為后續(xù)進(jìn)行計(jì)算并行化的基礎(chǔ)進(jìn)行二次開(kāi)發(fā)。

5 結(jié)束語(yǔ)

本研究針對(duì)牽引網(wǎng)提出了統(tǒng)一電磁暫態(tài)模型體系，建立了包括牽引網(wǎng)導(dǎo)線、牽引變壓器、自耦變壓器在內(nèi)的電氣設(shè)備的電磁暫態(tài)模型。并提出了基于牽引網(wǎng)電氣設(shè)備分布，分割牽引網(wǎng)導(dǎo)線，建立牽引網(wǎng)全線的電磁暫態(tài)節(jié)點(diǎn)電壓方程的計(jì)算方法。該方法能夠根據(jù)針對(duì)不同供電接線方式的牽引網(wǎng)系統(tǒng)建立全線的電磁暫態(tài)模型。但本研究的牽引網(wǎng)電磁暫態(tài)模型主要集中在理論性的探討，所采用的牽引網(wǎng)電器元件模型相對(duì)比較簡(jiǎn)單，與實(shí)際牽引供電系統(tǒng)還有一定的差距。

［1］李 芾，安 琪，王 華.高速動(dòng)車組概論［M］.成都:西南交通大學(xué)出版社，2008.

［2］何正友，方 雷，郭 東，等.基于AT 等值電路的牽引網(wǎng)潮流計(jì)算方法［J］.西南交通大學(xué)學(xué)報(bào)，2008(2):1-6.

［3］郭 東，楊建維，何正友，等.一種基于牛頓法的交流高速鐵路牽引供電潮流計(jì)算方法的研究［J］.繼電器，2007，35(18):16-29.

［4］吳命利.電氣化鐵道牽引網(wǎng)的統(tǒng)一鏈?zhǔn)诫娐纺Ｐ停跩］.中國(guó)電機(jī)工程學(xué)報(bào)，2010，30(28):52-58.

［5］陳宏偉，耿光超，江全元.電氣化鐵路牽引供電系統(tǒng)車網(wǎng)耦合的潮流計(jì)算方法［J］.電力系統(tǒng)自動(dòng)化，2012，36(3):76-80.

［6］LEE H M，LEE C，JANG G.Harmonic analysis of the korean high-speed railway using the eight-port representation model［J］.IEEE Trans on Power Delivery，2006，21(2):979-986.

［7］張 揚(yáng)，劉志剛.基于電磁暫態(tài)分析的高速鐵路牽引網(wǎng)諧波模型及諧波特性分析［J］.電網(wǎng)技術(shù)，2011，35(5):70-75.

［8］李群湛，賀建閩.牽引供電系統(tǒng)分析［M］.成都:西南大學(xué)出版社，2007.

［9］李保霞.行李車DC 600 V 供電系統(tǒng)改造設(shè)計(jì)［J］.2014(S1):157-158.

［10］譚秀炳.交流電氣化鐵道牽引供電系統(tǒng)［M］.成都:西南交通大學(xué)出版社，2006.

［11］DOMMEL H W.Computation of electromagnetic transients［J］.Proc IEEE，1974，62(7):983-993.

［12］DOMMEL H W.Digital computer solution of electromagnetic transients in single and multiphase networks［J］.IEEE Trans on Power Apparatus and Systems，1969，88(4):388-399.

［13］吳維韓，張芳榴，刁頤民.多導(dǎo)線輸電線路上波過(guò)程的貝杰龍計(jì)算方法［J］.高電壓技術(shù)，1981(4):9-24.

［14］徐 政.耦合長(zhǎng)線電磁暫態(tài)分析的擴(kuò)張Bergeron 模型［J］.中國(guó)電機(jī)工程學(xué)報(bào)，1996，16(5):354-357.

［15］NGUYEN H V，DOMMEL H W，MARTI J R.Direct phase-domain modelling of frequency-dependent overhead transmission lines［J］.IEEE Trans actions on Power Delivery，1997，12(3):1335-1342.

［16］LEON F D，SEMLYEN A.Complete Transformer Model for Electromagnetic Transients［J］.IEEE Trans actions on Power Delivery，1994，9(1):231-239.

［17］彭龍生.暫態(tài)保護(hù)AD 模塊開(kāi)發(fā)［J］.機(jī)電工程技術(shù)，2014(12):111-113.

［18］趙亮亮.電力系統(tǒng)電磁暫態(tài)數(shù)字仿真中自耦變壓器模型的研究［D］.天津:天津大學(xué)電氣工程學(xué)院，2004.

［19］朱 翊.電力變壓器電磁暫態(tài)仿真模型與算法研究［D］.天津:天津大學(xué)電氣工程學(xué)院，2012.

［20］王慶平，董新洲，周雙喜.基于自適應(yīng)模型的變壓器暫態(tài)全過(guò)程數(shù)值計(jì)算［J］.電力系統(tǒng)自動(dòng)化，2004，28(18):54-58.

［21］余丹萍.電氣化鐵路牽引供電系統(tǒng)的仿真及影響研究［D］.杭州:浙江大學(xué)電氣工程學(xué)院，2011.

［22］余丹萍，周 盛，江全元.CRH3 型動(dòng)車組牽引傳動(dòng)系統(tǒng)的直接轉(zhuǎn)矩控制研究［J］，機(jī)電工程，2010，27(10):62-67.