簡介面向軌道交通電氣系統(tǒng)設(shè)計(jì)、維護(hù)及相關(guān)研究人員的MATLAB/Simulink建模仿真資源圍繞地鐵牽引供電系統(tǒng)展開涵蓋外部電源、24脈波整流變電站、直流母線、接觸網(wǎng)與鋼軌大地模型以及含逆變電路和牽引負(fù)荷的牽引傳動系統(tǒng)?；诔啥嫉罔F實(shí)際參數(shù)進(jìn)行空載、負(fù)載與故障工況仿真并與現(xiàn)場故障錄波對比驗(yàn)證模型有效性可幫助工程師理解建模方法、優(yōu)化系統(tǒng)設(shè)計(jì)并開展故障分析。資源為單個(gè)PDF文件共671KB內(nèi)容包含論文復(fù)現(xiàn)所需的詳細(xì)可運(yùn)行代碼及逐段中文解釋涉及24脈波整流移相變壓器組、接觸網(wǎng)PI型分布參數(shù)模型、SVPWM逆變控制及系統(tǒng)級仿真框架等關(guān)鍵技術(shù)并體現(xiàn)了分層驗(yàn)證體系與多時(shí)間尺度仿真策略。已有244人學(xué)習(xí)下載適合需要對照論文落地仿真、掌握地鐵牽引供電建模與故障仿真方法的工程技術(shù)及研究人員。1. 地鐵牽引供電仿真為什么值得啃下這篇論文復(fù)現(xiàn)地鐵牽引供電系統(tǒng)是城市軌道交通的能量大動脈而基于MATLAB/Simulink的建模與仿真一直是碩博論文和工程驗(yàn)證里出現(xiàn)頻率最高的方向之一。這篇標(biāo)題對應(yīng)的復(fù)現(xiàn)工作核心不是把別人的圖重新畫一遍而是用Simulink把牽引變電所、整流機(jī)組、接觸網(wǎng)、鋼軌回流和列車負(fù)載串成一個(gè)能跑的動態(tài)系統(tǒng)再人為注入短路、整流器故障等工況從波形里讀出保護(hù)配合邏輯和故障特征。對做供電專業(yè)的學(xué)生來說它是打通電路理論、電力電子和繼電保護(hù)三塊知識的最佳實(shí)踐載體對現(xiàn)場工程師而言一套可調(diào)參的仿真模型能提前驗(yàn)證保護(hù)定值和運(yùn)行方式變更的影響遠(yuǎn)比直接動真設(shè)備安全。我最早做這個(gè)方向時(shí)被Simulink里整流橋的數(shù)值震蕩折磨過也因忽略了鋼軌散流模型導(dǎo)致故障電流和實(shí)測對不上。這篇文章會把建模流程、關(guān)鍵參數(shù)、故障注入方法和常見坑一次講透代碼段都能直接復(fù)制到你的Simulink環(huán)境里改參數(shù)跑通希望幫你在復(fù)現(xiàn)這條路上少走幾天彎路。2. 從牽引變電所到接觸網(wǎng)Simulink模型拓?fù)渑c最小可跑通方案2.1 牽引供電系統(tǒng)的基本組成與仿真模型選型地鐵牽引供電系統(tǒng)的主干結(jié)構(gòu)并不復(fù)雜高壓交流電源通常110kV或35kV進(jìn)入主變電所降壓后供給牽引變電所牽引變電所內(nèi)通過整流機(jī)組把交流變成直流1500V部分線路用750V輸出到接觸網(wǎng)或第三軌列車通過受電弓取流電流經(jīng)鋼軌和回流線流回牽引變電所。模型里必須把這幾個(gè)環(huán)節(jié)都放進(jìn)去否則故障分析時(shí)回路不完整短路電流路徑會少一條。在MATLAB/Simulink里做這個(gè)系統(tǒng)選型上有三條路一是用Simscape Electrical舊稱SimPowerSystems里的基礎(chǔ)電氣元件自己搭整流橋、變壓器和直流母線靈活度最高二是用Simscape Electrical的Specialized Power Systems庫里的“Universal Bridge”和“Three-Phase Transformer”模塊搭起來快內(nèi)部開關(guān)邏輯已經(jīng)封裝好三是用Simscape基礎(chǔ)元件庫做物理域建模但和傳統(tǒng)Simulink控制信號交互稍麻煩。我一般選第二種原因是論文復(fù)現(xiàn)時(shí)整流機(jī)組的脈波數(shù)、變壓器聯(lián)結(jié)組別都需要顯式設(shè)置Universal Bridge可以直接指定橋臂數(shù)和開關(guān)器件類型最貼近技術(shù)論文里的描述。模型拓?fù)湮医ㄗh按“交流電源—整流變壓器—整流橋—直流母線—接觸網(wǎng)—列車負(fù)載—鋼軌回流”的單線圖來搭。為了后續(xù)能注入故障直流母線正極和接觸網(wǎng)之間加一個(gè)分段開關(guān)接觸網(wǎng)末端預(yù)留故障點(diǎn)。仿真步長不用全局固定Simscape Electrical會自動變步長只需要把最大步長限制在10微秒以內(nèi)避免整流橋換相時(shí)數(shù)值突變引起步長收縮過慢。2.2 用Simscape Electrical搭建整流機(jī)組與直流母線第一步是在Simulink里新建模型拖入以下模塊Three-Phase Source三相交流源、Three-Phase Transformer整流變壓器、Universal Bridge整流橋、Series RLC Branch直流平波電抗器、DC Voltage Measurement和Scope。交流源設(shè)為35kV/50Hz變壓器二次側(cè)線電壓設(shè)為1180V左右——這個(gè)值是為了讓空載直流電壓達(dá)到1500V的1.35倍關(guān)系即約2025V空載電壓帶負(fù)載后自然跌落到1500V左右。整流橋選“Universal Bridge”Power Electronic Device選“Diode”橋臂數(shù)選“3”這樣就是典型的三相6脈波整流。變壓器聯(lián)結(jié)組別是關(guān)鍵工程上常用Dd0或Dy11。6脈波整流用Dd0和Dy11均可但12脈波必須用兩組變壓器分別接Dd0和Dy11利用相位差30度來抵消5、7次諧波。如果論文里討論24脈波就要在變壓器二次側(cè)加延邊三角形移相繞組。復(fù)現(xiàn)時(shí)先跑通6脈波再擴(kuò)展到12脈波驗(yàn)證步驟更清晰。下面是最小模型的參數(shù)腳本可以寫在模型初始化回調(diào)里% 牽引供電系統(tǒng)最小模型參數(shù)初始化 % 交流電源部分 V_line_rms 35e3; % 電源線電壓有效值單位V f_grid 50; % 電網(wǎng)頻率單位Hz R_source 0.5; % 電源等效內(nèi)阻單位ohm L_source 5e-3; % 電源等效電感單位H % 整流變壓器參數(shù) V_primary 35e3; % 一次側(cè)線電壓單位V V_secondary 1180; % 二次側(cè)線電壓單位V空載直流電壓約2025V R_leak 0.02; % 變壓器短路電阻標(biāo)幺值或折算到二次側(cè)的歐姆值 L_leak 0.15e-3; % 變壓器短路電感單位H % 整流橋與直流側(cè) cnv_pulses 6; % 整流脈波數(shù)6脈波基礎(chǔ)模型 L_smooth 5e-3; % 直流平波電抗器單位H R_smooth 0.01; % 平波電抗器內(nèi)阻單位ohm C_dc 10e-3; % 直流支撐電容單位F論文中有時(shí)不裝按需開啟這段腳本的關(guān)鍵是變壓器二次側(cè)電壓的取值邏輯三相橋式整流電路空載直流電壓是二次側(cè)線電壓有效值的1.35倍所以要讓空載電壓約2025V二次側(cè)線電壓就是2025/1.35≈1500V這里其實(shí)有個(gè)常見誤解。工程上牽引變壓器二次側(cè)額定電壓取1180V~1200V是考慮負(fù)載壓降和整流管壓降后的設(shè)計(jì)值空載直流電壓約1600V上下滿載時(shí)直流電壓約1500V。不同論文取值有差異復(fù)現(xiàn)時(shí)最好先按論文給定的空載電壓折算回去。如果論文沒給就用二次側(cè)1200V起步跑負(fù)載看直流電壓是否落在1500V附近再微調(diào)。2.3 列車負(fù)載與線路參數(shù)的模擬方法列車負(fù)載在牽引供電仿真里是最容易“湊合”的部分但它直接決定直流電流波形。最簡單的做法是恒定電阻負(fù)載用直流電壓除以期望電流算電阻值更接近實(shí)際的用受控電流源或Simscape的DC-DC變換器模型模擬列車牽引時(shí)的恒功率特性和制動時(shí)的回饋特性。論文復(fù)現(xiàn)建議分步走第一步用恒定電阻跑通系統(tǒng)第二步換成可變電阻模擬列車啟動—加速—惰行—制動四象限工況。線路參數(shù)不能只給一個(gè)純電阻接觸網(wǎng)和鋼軌都有單位長度電阻、電感和對地電容。1500V直流牽引網(wǎng)通常用走行軌回流鋼軌對地有泄漏電導(dǎo)這個(gè)參數(shù)在故障分析中影響站間距離與短路電流的對應(yīng)關(guān)系。在Simulink里可以用“PI Section Line”模塊或者多段串聯(lián)的RLC分支來模擬。每段線路長度按站間距取值例如1.5km一段每段電阻0.01Ω/km、電感0.8mH/km、對地電導(dǎo)0.5S/km潮濕隧道環(huán)境取值會更高。注意鋼軌回流路徑必須和接觸網(wǎng)正極一起構(gòu)成回路否則故障時(shí)仿真會提示“Ground connection”問題。列車負(fù)載模擬我常用一個(gè)受控電流源輸入信號是列車牽引電流曲線。這個(gè)電流曲線可以用Lookup Table導(dǎo)入實(shí)測或論文給定的數(shù)據(jù)點(diǎn)。比如某列車啟動電流從0平緩升到1500A維持到速度拐點(diǎn)后線性下降。生成這個(gè)曲線的腳本如下% 生成列車牽引電流曲線存入Lookup Table t [0 5 10 15 20 25 30]; % 時(shí)間軸單位s I_train [0 600 1500 1500 900 300 0]; % 牽引電流單位A % 在Simulink中創(chuàng)建Lookup Table時(shí)直接使用這兩組數(shù)據(jù) plot(t, I_train, o-); xlabel(時(shí)間 (s)); ylabel(牽引電流 (A)); title(列車牽引電流參考曲線); grid on;電流源的輸出接在接觸網(wǎng)正極和鋼軌回流之間方向是從接觸網(wǎng)流向鋼軌正號代表牽引取流負(fù)號代表制動回饋。Simulink里的受控電流源需要并聯(lián)一個(gè)大電阻比如1MΩ避免數(shù)值孤立這個(gè)并聯(lián)電阻不影響正常工況但在短路故障瞬間會產(chǎn)生一個(gè)微小分流分析故障電流時(shí)需要減掉這部分誤差。列車位置的移動性可以用多個(gè)電流源加開關(guān)在不同時(shí)間接入來實(shí)現(xiàn)但初期模型里固定一個(gè)位置就夠了論文復(fù)現(xiàn)的重點(diǎn)是波形的“形態(tài)”不是精確的時(shí)空分布。3. 整流機(jī)組與列車負(fù)載兩個(gè)決定仿真精度的關(guān)鍵技術(shù)3.1 整流機(jī)組等效模型與脈波數(shù)選擇整流機(jī)組是地鐵牽引供電系統(tǒng)里最核心的電力電子環(huán)節(jié)它的脈波數(shù)直接決定了直流側(cè)電壓脈動和注入交流系統(tǒng)的諧波含量。6脈波整流是最基礎(chǔ)的結(jié)構(gòu)直流電壓紋波較大5次、7次諧波突出12脈波由兩組6脈波整流器通過移相30度并聯(lián)或串聯(lián)組成5、7次諧波被大幅抵消直流紋波也明顯變小24脈波則是在12脈波基礎(chǔ)上再加移相繞組工程上用于大型主變電所。在Simulink里實(shí)現(xiàn)12脈波常見做法是使用兩臺三繞組變壓器一次側(cè)并聯(lián)二次側(cè)一組Y接、一組Δ接分別帶各自的整流橋然后直流側(cè)串聯(lián)輸出。注意兩套整流橋的直流輸出電壓必須極性一致且幅度接近否則會形成環(huán)流。參數(shù)上兩組變壓器的變比必須嚴(yán)格一致短路阻抗誤差盡量控制在±5%以內(nèi)。我用Universal Bridge分別搭兩套整流橋橋臂數(shù)和器件類型保持一致只改變變壓器聯(lián)結(jié)組別。為了在論文里定量說明脈波數(shù)的影響可以在模型里加一個(gè)FFT Analysis工具穩(wěn)態(tài)運(yùn)行后測量直流電壓紋波系數(shù)和諧波頻譜。紋波系數(shù)的定義是直流電壓交流分量有效值與直流平均值之比。6脈波理論上紋波系數(shù)約4.2%12脈波降到1.0%左右24脈波進(jìn)一步降到0.26%。如果你復(fù)現(xiàn)的論文給出了紋波系數(shù)對比表仿真值應(yīng)該落在理論值附近差太多就檢查整流橋觸發(fā)角是否為零或變壓器漏抗是否過大。整流機(jī)組的等效模型還有一層IGBT或晶閘管開關(guān)的動態(tài)特性。在故障暫態(tài)分析里開關(guān)器件的關(guān)斷延遲和反向恢復(fù)過程會影響短路電流的初始上升率。Universal Bridge默認(rèn)的Snubber電路參數(shù)也需要注意默認(rèn)Rs1e5Ω、Csinf這相當(dāng)于把開關(guān)理想化。若要更真實(shí)的故障電流需要把Snubber電路的電容設(shè)為有限值例如10nF電阻降到10Ω但這會顯著增加仿真時(shí)長建議只在最終驗(yàn)證階段開啟。3.2 列車牽引/制動特性曲線怎么進(jìn)模型列車負(fù)載的動態(tài)特性直接影響直流母線電壓的波動范圍。地鐵列車的運(yùn)行工況分為牽引、巡航、惰行、制動四象限其中牽引時(shí)列車從接觸網(wǎng)吸取功率制動時(shí)一部分能量通過再生制動回饋電網(wǎng)。Simulink里最直觀的建模方式是用受控電壓源和受控電流源組合牽引時(shí)電流方向?yàn)檎苿訒r(shí)為負(fù)電壓始終由直流網(wǎng)側(cè)決定。我把列車模型拆成兩部分功率計(jì)算和電流注入。功率計(jì)算根據(jù)列車速度、加速度和質(zhì)量按簡化力學(xué)公式算出瞬時(shí)牽引力再乘以車速得到機(jī)械功率除以效率后得到電功率電流注入則把這個(gè)電功率除以直流母線電壓得到參考電流送給受控電流源。這個(gè)模型能體現(xiàn)恒功率特性——母線電壓下降時(shí)電流會上升和實(shí)際列車特性一致。下面是一個(gè)最小實(shí)現(xiàn)function I_ref train_current(v, a, m, eta, Vdc) % 列車牽引電流計(jì)算函數(shù) % v 車速 m/sa 加速度 m/s^2m 質(zhì)量 kgeta 效率Vdc 直流母線電壓 V % 基本阻力簡化為 A B*v C*v^2 A 4500; % 滾動阻力系數(shù) N B 90; % 線性阻力系數(shù) N·s/m C 15; % 空氣阻力系數(shù) N·s^2/m^2 F m * a A B*v C*v^2; % 所需牽引力 N if F 0 P_mech F * v; % 機(jī)械功率 W P_elec P_mech / eta; % 牽引時(shí)電功率大于機(jī)械功率 else P_mech F * v; % 制動力F為負(fù)則P_mech為負(fù) P_elec P_mech * eta; % 回饋時(shí)電功率小于機(jī)械功率且為負(fù) end I_ref P_elec / Vdc; % 直流側(cè)電流參考值 A end在Simulink里用MATLAB Function模塊調(diào)用這個(gè)函數(shù)輸入車速、加速度和直流電壓測量值輸出I_ref給受控電流源。注意初始時(shí)刻車速為0時(shí)理論上牽引力公式里v0會導(dǎo)致功率為0但實(shí)際列車啟動時(shí)會有靜止?fàn)恳λ孕枰谀Ｐ屠锛右粋€(gè)啟動邏輯當(dāng)v0.1m/s且F0時(shí)按恒力矩處理即I_ref用F乘以一個(gè)等效半徑換算的常數(shù)代替。忽略這個(gè)邏輯會導(dǎo)致仿真啟動瞬間電流參考值跳變到0列車根本沒法“動起來”這個(gè)坑我踩過兩次。3.3 直流側(cè)短路故障的Simulink仿真實(shí)現(xiàn)直流側(cè)短路是地鐵供電系統(tǒng)最嚴(yán)重的故障類型接觸網(wǎng)正極對鋼軌直接短接短路電流可達(dá)數(shù)千安到數(shù)十千安。仿真中實(shí)現(xiàn)短路最簡單的方式是用一個(gè)斷路器Three-Phase Breaker或Ideal Switch加一個(gè)很小電阻在指定時(shí)間閉合。為了模擬不同短路位置的過渡電阻可以用Breaker串聯(lián)一個(gè)可變電阻。故障時(shí)序設(shè)定先讓系統(tǒng)空載或帶負(fù)載穩(wěn)定運(yùn)行到0.1s再在0.1s時(shí)閉合短路支路持續(xù)0.15s后由保護(hù)動作斷開。保護(hù)動作在仿真里用Breaker的打開時(shí)間模擬也可以用Simulink的Step模塊觸發(fā)開關(guān)。關(guān)鍵是要記錄短路發(fā)生前后直流電壓和電流的波形觀察沖擊電流峰值、電壓跌落幅度和恢復(fù)到正常的時(shí)間。注入故障時(shí)我建議檢查三個(gè)細(xì)節(jié)第一短路支路的閉合電阻不能設(shè)置為0否則Simulink會出現(xiàn)代數(shù)環(huán)或數(shù)值奇異一般取1e-3Ω第二故障點(diǎn)與整流機(jī)組之間必須有線路阻抗否則短路電流峰值會遠(yuǎn)超真實(shí)值因?yàn)樽儔浩髀┛购椭绷麟娍蛊魇窍拗贫搪冯娏鞯闹饕绻€路阻抗被忽略了峰值電流可能高估30%以上第三故障后的恢復(fù)過程要保留平波電抗器它決定短路電流上升率。下面是一個(gè)簡單的故障控制腳本可以用作Model Callback也可以放在Constant模塊的Signal Builder里% 直流側(cè)短路故障控制 % 使用Simulink Signal Builder或Step模塊本腳本生成故障開關(guān)時(shí)間序列 t_step [0 0.1 0.1 0.25 0.25]; % 時(shí)間節(jié)點(diǎn) Brk_state [0 0 1 1 0]; % 0表示打開1表示閉合 % 對應(yīng)的斷路器控制信號在Simulink中通過Repeating Sequence模塊導(dǎo)入 % 短路電阻取0.005ohm閉合后經(jīng)歷0.15s的保護(hù)動作時(shí)間實(shí)際運(yùn)行時(shí)會發(fā)現(xiàn)短路瞬間直流電壓從1500V驟降到接近0直流電流迅速攀升到峰值隨后由于交流側(cè)阻抗和直流電感的作用逐漸衰減到穩(wěn)態(tài)短路電流。這個(gè)波形的上升沿斜率主要受直流回路電感限制而穩(wěn)態(tài)值受回路電阻限制。論文里如果給出了短路電流的解析計(jì)算式仿真峰值和解析值對不上時(shí)優(yōu)先排查故障點(diǎn)位置和直流側(cè)總電感是否正確。4. 故障分析把短路電流波形和繼電保護(hù)動作講清楚4.1 整流器內(nèi)部短路故障分析整流器內(nèi)部故障常見有整流二極管擊穿短路、同一橋臂上下管直通、交流側(cè)相間短路等。在Simulink里模擬二極管擊穿最簡單的做法是在被擊穿的二極管兩端并聯(lián)一個(gè)受控電壓源或小電阻開關(guān)但這種建模太理想化不能反映故障對交流側(cè)的影響。更接近實(shí)際的建模方法是在整流橋外面并聯(lián)故障支路故障點(diǎn)設(shè)在整流橋直流輸出端模擬逆變失敗導(dǎo)致的直流母線短路。另一種常見故障是整流變壓器二次側(cè)匝間短路。由于變壓器模型在Simscape Electrical里是基于理想變壓器的匝間短路需要把變壓器拆成多繞組模型在二次繞組上并聯(lián)一個(gè)小電阻來模擬局部短路。復(fù)現(xiàn)論文里討論整流器故障時(shí)重點(diǎn)看的是故障后的特征交流側(cè)電流出現(xiàn)嚴(yán)重不平衡直流電壓跌落且伴有較大諧波直流電流出現(xiàn)周期性尖峰。我在仿真整流器內(nèi)部故障時(shí)發(fā)現(xiàn)一個(gè)需要注意的現(xiàn)象整流橋內(nèi)部故障如果不改變觸發(fā)條件仿真器會自動報(bào)“Singularity”錯(cuò)誤因?yàn)楣收现钒央娫粗苯佣搪妨?。解決辦法是在故障支路串入一個(gè)0.01Ω小電阻并把仿真器最大步長放寬到1e-5s。如果仍然不收斂把整流橋的Snubber電路電阻調(diào)大比如從1e5Ω調(diào)到1e6Ω。4.2 接觸網(wǎng)短路與保護(hù)動作配合接觸網(wǎng)短路是地鐵系統(tǒng)中最需要保護(hù)配合的故障因?yàn)榱熊嚥粩嘁苿庸收宵c(diǎn)位置不同導(dǎo)致短路電流大小不同。仿真里通常設(shè)置幾個(gè)典型位置牽引變電所出口、線路中點(diǎn)、末端。出口短路時(shí)短路電流最大末端短路時(shí)電流最小保護(hù)需要通過電流大小和距離信息來區(qū)分故障區(qū)段避免越級跳閘。Simulink里實(shí)現(xiàn)多位置短路很簡單把接觸網(wǎng)分成若干段每段末端接一個(gè)斷路器到鋼軌不同時(shí)刻閉合不同的斷路器就行。但要注意各短路點(diǎn)不應(yīng)該同時(shí)閉合否則會重合短路。合理的做法是用一個(gè)故障序列控制邏輯例如0.1s在第一個(gè)點(diǎn)短路0.3s斷開后又切換第二個(gè)點(diǎn)。保護(hù)動作配合仿真可以用簡單的繼電器邏輯實(shí)現(xiàn)測量直流牽引網(wǎng)電流超過設(shè)定值后延時(shí)動作斷開牽引變電所出口斷路器。為了模擬整定值的配合關(guān)系可以分別設(shè)置兩段保護(hù)——瞬時(shí)速斷和定時(shí)過流。在Simulink里用Relay模塊加延時(shí)環(huán)節(jié)搭輸出觸發(fā)斷路器的信號。關(guān)鍵是把保護(hù)動作時(shí)間設(shè)定得比故障持續(xù)時(shí)間短否則斷路器不會斷開仿真會一直維持短路狀態(tài)。比較有指導(dǎo)意義的仿真實(shí)驗(yàn)是調(diào)整保護(hù)定值和線路參數(shù)觀察保護(hù)動作的覆蓋范圍。比如末端短路電流為1500A如果速斷保護(hù)定值設(shè)在3000A那么末端故障只能靠過流保護(hù)定值1200A延時(shí)0.5s動作導(dǎo)致故障切除時(shí)間更晚。論文里討論保護(hù)配合不合理時(shí)通常就是拿這種仿真結(jié)果作對比說明不同定值下保護(hù)范圍的變化。4.3 故障波形怎么看電壓、電流、功率的異常特征故障分析最終要落到波形判讀上。直流電壓、電流波形是最直觀的但功率波形、諧波頻譜、交流側(cè)電流也要一起看。整流器出口短路時(shí)直流電壓瞬間塌陷電流沖高后回落功率波形會出現(xiàn)一個(gè)正向尖峰后持續(xù)低值交流側(cè)電流則出現(xiàn)明顯的不對稱某相電流幅值異常增大。我還有幾個(gè)讀波形的習(xí)慣第一不要只看瞬時(shí)值用RMS轉(zhuǎn)換模塊計(jì)算故障期間的電流有效值注意Simulink里到的RMS是滑動窗口值窗口長度要設(shè)置為工頻周期20ms否則波形平滑度不夠無法和工程上的保護(hù)動作值對應(yīng)第二用Power Spectral Density或FFT工具看故障前后諧波變化整流器內(nèi)部故障通常會產(chǎn)生明顯的偶次諧波而外部短路主要產(chǎn)生直流分量和低頻分量第三把仿真數(shù)據(jù)和短路電流解析公式的計(jì)算結(jié)果畫在同一張圖里驗(yàn)證模型是否在可接受誤差范圍內(nèi)。如果論文中提到了“直流牽引網(wǎng)短路電流上升率”這一指標(biāo)你需要在仿真里測量的就是di/dt。Simulink里可以對電流信號用Derivative模塊求導(dǎo)然后輸出去找最大值。這個(gè)值受直流回路電感影響仿真中如果發(fā)現(xiàn)上升率比論文給的大很多去看看平波電抗器參數(shù)是否設(shè)置正確、接觸網(wǎng)線路電感是否被忽略了。電感量每減少1mHdi/dt可能上升近一倍參數(shù)敏感性很高。5. 避坑指南MATLAB/Simulink牽引供電仿真的5個(gè)常見坑與排查5.1 整流橋仿真不收斂報(bào)“Singularity”錯(cuò)誤現(xiàn)象仿真跑到整流橋換相瞬間突然報(bào)錯(cuò)提示“Derivative of state is infinite”或“Singularity”步長反復(fù)收縮到極小值后崩潰。原因整流橋開關(guān)器件切換時(shí)理想開關(guān)兩端電感電流突變Simulink的變步長求解器無法處理不連續(xù)狀態(tài)或者故障支路電阻設(shè)為了0直接把電容或電壓源短路。解決給所有開關(guān)器件和故障支路串入一個(gè)0.001~0.01Ω的小電阻并檢查直流側(cè)是否串有平波電抗器同時(shí)在求解器設(shè)置里把“Max Step Size”設(shè)為1e-5s把“Relay”模塊的切換容差調(diào)大一些。我一般還會把Universal Bridge的Snubber電阻從1e5Ω改成100Ω這能大幅緩解突變。5.2 直流電壓偏低空載值只有1000V而不是1500V現(xiàn)象模型搭建完成后交流側(cè)參數(shù)都對但直流電壓測量值遠(yuǎn)低于設(shè)計(jì)值而且波形的脈波數(shù)量不對。原因變壓器二次側(cè)電壓設(shè)置錯(cuò)誤或者整流橋橋臂數(shù)選錯(cuò)。比如應(yīng)該選3橋臂三相全橋卻選了1橋臂更常見的是把變壓器二次側(cè)線電壓和相電壓搞混計(jì)算空載直流電壓時(shí)應(yīng)該用線電壓乘1.35而不是相電壓。解決先斷開所有負(fù)載和故障支路單獨(dú)測試整流機(jī)組空載輸出。確認(rèn)變壓器二次側(cè)額定電壓是否正確用萬用表模塊分別測量二次側(cè)相電壓和線電壓。如果線電壓是1180V空載直流電壓應(yīng)該約1590V1180×1.35在正常范圍。若實(shí)測值只有約1400V檢查整流橋是否選成了“半控”模式或觸發(fā)電平有問題。5.3 故障電流峰值比論文小很多甚至只有一半現(xiàn)象用斷路器模擬直流側(cè)短路理論估算短路電流峰值約9kA仿真測量只有4.5kA。原因常見原因有三個(gè)一是直流側(cè)平波電抗器電感值設(shè)置過大限制了電流上升率二是交流電源等效內(nèi)阻和變壓器漏抗被隱式放大了比如變壓器勵(lì)磁電感默認(rèn)值過小三是故障點(diǎn)位置離整流機(jī)組太遠(yuǎn)線路阻抗分掉了大量壓降。解決把平波電抗器電感從5mH降到2mH看峰值是否上升再檢查變壓器“Leakage inductance”參數(shù)是否與實(shí)際短路阻抗匹配。還有一個(gè)隱性參數(shù)是“Three-Phase Source”內(nèi)部的“Source internal resistance”默認(rèn)0.01Ω偏大時(shí)也會限制短路電流。把這個(gè)參數(shù)改成0.001Ω同時(shí)確認(rèn)故障斷路器的導(dǎo)通電阻不超過1mΩ。5.4 仿真速度極慢跑一秒要十分鐘現(xiàn)象模型加上整流橋、列車負(fù)載、線路分段后仿真速度急劇下降步長始終小于1e-7s。原因電氣系統(tǒng)中存在極小的時(shí)間常數(shù)比如大電容和極小電阻串聯(lián)或開關(guān)器件的緩沖電路阻容時(shí)間常數(shù)太小迫使求解器采用極小步長。另一個(gè)常見原因是過度使用了理想開關(guān)開關(guān)頻率太高。解決檢查所有串聯(lián)電阻和并聯(lián)電容組合算出時(shí)間常數(shù)盡量讓它們不小于1e-6s把Universal Bridge的Snubber電路電容設(shè)為Inf即無緩沖可以顯著提速前提是建模重點(diǎn)是穩(wěn)態(tài)和保護(hù)配合而不是開關(guān)過程精確波形。另外把線路模型從分布式參數(shù)改成集中參數(shù)用RLC串聯(lián)分支代替PI Section Line速度能提升數(shù)倍。5.5 三相交流側(cè)電流波形嚴(yán)重畸變諧波含量異?，F(xiàn)象交流側(cè)電流波形不是預(yù)期的近似正弦波而是尖峰密集的畸波FFT結(jié)果里諧波高達(dá)20%以上與論文的仿真圖完全不符。原因整流變壓器二次側(cè)繞組聯(lián)結(jié)錯(cuò)誤或者整流橋觸發(fā)邏輯用錯(cuò)了PWM信號。牽引供電整流是不控整流二極管應(yīng)該自然換相如果你錯(cuò)誤接入了觸發(fā)脈沖等于給二極管強(qiáng)行加上了開關(guān)信號導(dǎo)致電流斷續(xù)。解決檢查Universal Bridge的“Power Electronic Device”是否確實(shí)選為“Diode”而不是“IGBT/Diode”。IGBT混合模型默認(rèn)允許外部觸發(fā)信號未接線時(shí)可能處在錯(cuò)誤的門控狀態(tài)。另一個(gè)原因是交流側(cè)電源內(nèi)阻太小且無濾波電抗器整流器換相時(shí)產(chǎn)生尖峰。在交流側(cè)每個(gè)相線上串入一個(gè)小電感約0.1mH來模擬實(shí)際變壓器漏抗波形會立即平滑下來。6. 進(jìn)階把仿真結(jié)果變成論文復(fù)現(xiàn)證據(jù)的三種做法到這里模型已經(jīng)能跑但復(fù)現(xiàn)的關(guān)鍵是讓仿真曲線和論文結(jié)果“對上話”。我的習(xí)慣是先做“參數(shù)敏感性驗(yàn)證”再做“故障特征對比表”最后生成“三圖合一”的復(fù)合波形。先做一個(gè)簡單的腳本批量修改平波電抗器電感值記錄直流電壓紋波系數(shù)和短路電流峰值生成對比表。這個(gè)表直接放進(jìn)論文里能說明模型對參數(shù)的響應(yīng)符合物理規(guī)律比單張波形更有說服力。腳本如下% 批量掃描平波電抗器電感記錄關(guān)鍵指標(biāo) L_values [1e-3, 2e-3, 5e-3, 8e-3]; results zeros(length(L_values), 3); for i 1:length(L_values) set_param(traction_sim/L_smooth, L, num2str(L_values(i))); sim(traction_sim); results(i, 1) L_values(i); results(i, 2) Vdc_ripple_ratio; % 從工作區(qū)讀取紋波系數(shù) results(i, 3) Isc_peak; % 短路電流峰值 end disp(table(results(:,1), results(:,2), results(:,3), ... VariableNames, {L_H, ripple_ratio, Isc_peak}));第二種做法是驗(yàn)證相對誤差。找到論文里給出的穩(wěn)態(tài)直流電壓、短路電流峰值或諧波含量數(shù)據(jù)把仿真值算成相對誤差控制在±10%以內(nèi)即可認(rèn)為模型可信。如果超過這個(gè)范圍優(yōu)先懷疑模型參數(shù)與論文不一致特別是變壓器短路阻抗、線路單位長度參數(shù)和列車負(fù)載特性曲線。第三種做法是把仿真波形、解析計(jì)算波形和論文截圖如果你有版權(quán)許可放在同一張圖里對比。沒有論文原圖時(shí)至少要把仿真波形和解析公式曲線對比比如空載直流電壓、短路穩(wěn)態(tài)電流都可以手算。這樣審查老師會認(rèn)為你的模型經(jīng)過了雙向驗(yàn)證不是“黑匣子”跑出來的。我做這類復(fù)現(xiàn)的一個(gè)深刻教訓(xùn)是仿真波形“像不像”論文是次要的“為什么像”才是重點(diǎn)。某次我在整流器故障分析里得到一條和論文幾乎一致的電壓波形后來仔細(xì)檢查才發(fā)現(xiàn)是因?yàn)槲以谡鳂蜉敵龆瞬⒘艘粋€(gè)大電容恰好把這個(gè)故障特征“濾”得和論文一樣了。這個(gè)電容在真實(shí)系統(tǒng)里不存在純粹是我為了波形好看加上去的。后來我刪掉電容波形露出了原本的尖峰和振蕩反而更像真實(shí)故障?，F(xiàn)在我的習(xí)慣是每加一個(gè)元件前先問自己真實(shí)系統(tǒng)里有沒有這個(gè)東西如果沒有那就說明模型建錯(cuò)了而不是論文畫錯(cuò)了。希望這篇筆記能讓你在MATLAB/Simulink牽引供電建模與故障分析這條路上少一點(diǎn)玄學(xué)多一點(diǎn)可復(fù)現(xiàn)的確定性。如果你跑通了自己的第一個(gè)模型再去對比論文的波形、保護(hù)動作時(shí)序和故障特征表你會發(fā)現(xiàn)那些看似復(fù)雜的圖其實(shí)都是同一套原理在不同參數(shù)下的投影祝搭建順利。本文還有配套的精品資源點(diǎn)擊獲取