簡介直驅(qū)風(fēng)機永磁直驅(qū)式風(fēng)力發(fā)電系統(tǒng)完整Simulink仿真模型適用于風(fēng)電控制、電力電子及新能源發(fā)電方向的研究者與工程師。模型涵蓋從風(fēng)輪捕獲風(fēng)能到電能輸出的整體仿真流程并提供380V與690V兩種電壓等級方案分別對應(yīng)不同規(guī)模應(yīng)用場景。仿真波形可直觀反映風(fēng)速變化下的功率、轉(zhuǎn)速等關(guān)鍵參數(shù)便于性能評估與參數(shù)優(yōu)化。資源共13個文件包含txt模型說明與分析文檔、jpg波形截圖、doc文獻(xiàn)文檔及html輔助頁面整體壓縮包約506KB。支持2018與2021兩個Matlab版本并可導(dǎo)出2015b至2022版本兼容性良好。目前已有136人學(xué)習(xí)下載通過配套說明和文獻(xiàn)可快速掌握建模思路、模塊功能與仿真操作方法適合用于課題研究、課程設(shè)計或系統(tǒng)方案驗證。 做永磁直驅(qū)風(fēng)電Simulink整體仿真這個需求我估計在風(fēng)電、電力電子方向的圈子里出現(xiàn)頻率很高。很多剛接觸風(fēng)電仿真的人第一步就是找一個能跑通全過程的模型——從風(fēng)速輸入、風(fēng)輪機捕獲、永磁同步發(fā)電機PMSG出力到機側(cè)變流器、直流母線、網(wǎng)側(cè)逆變器最后并網(wǎng)或者帶負(fù)載整個鏈路完整跑起來。我自己也幫人調(diào)試過好幾套類似的模型一個比較典型的要求就是同一套系統(tǒng)分別給出380V和690V兩個電壓等級下的仿真一套接實驗室樣機驗證另一套用來演示兆瓦級并網(wǎng)場景。這篇就把這類模型的整體結(jié)構(gòu)、兩套電壓等級的差異、仿真調(diào)試中容易翻車的點以及我自己實測下來的操作細(xì)節(jié)一次說清楚。1. 直驅(qū)系統(tǒng)整體在仿什么從風(fēng)輪到電網(wǎng)的完整鏈路先明確一個概念永磁直驅(qū)風(fēng)力發(fā)電系統(tǒng)和傳統(tǒng)的雙饋異步風(fēng)電系統(tǒng)最大的區(qū)別在傳動鏈。雙饋電機轉(zhuǎn)子側(cè)通過滑環(huán)接入變頻器定子直接并網(wǎng)齒輪箱是標(biāo)配直驅(qū)系統(tǒng)則沒有齒輪箱風(fēng)輪機直接帶動永磁同步電機的轉(zhuǎn)子。因為轉(zhuǎn)速低所以電機極對數(shù)多、體積大但好處是機械結(jié)構(gòu)簡單、維護(hù)成本低、發(fā)電效率高低風(fēng)速下表現(xiàn)尤其好。Simulink里的整體仿真說直白點就是把這條能量鏈路的每一環(huán)都用數(shù)學(xué)模型搭出來模擬風(fēng)速變化時整個系統(tǒng)的動態(tài)響應(yīng)。一個完整模型大致包含這些部分風(fēng)輪機模型輸入風(fēng)速根據(jù)風(fēng)輪半徑、葉尖速比、槳距角計算空氣動力學(xué)轉(zhuǎn)矩把風(fēng)的動能變成機械能。永磁同步發(fā)電機PMSG把機械能變成電能定子側(cè)輸出的變頻變幅交流電頻率隨轉(zhuǎn)速變化。機側(cè)變流器MSC一個不可控整流橋或者PWM整流器把PMSG輸出的交流轉(zhuǎn)成直流。在高級一點的模型里用PWM整流可以實現(xiàn)最大功率追蹤。直流母線電容支撐電壓隔離機側(cè)和網(wǎng)側(cè)。網(wǎng)側(cè)變流器GSC把直流逆變成工頻交流控制功率流向。濾波器和變壓器濾除諧波電壓匹配最后并網(wǎng)或接負(fù)載。這套模型的價值在于它把空氣動力學(xué)、電機學(xué)、電力電子和控制理論全部揉在一起能觀察的不只是某一個模塊而是整個系統(tǒng)的功率平衡和動態(tài)響應(yīng)。2. PMSG和變流器建模時我建議你這樣搭參數(shù)2.1 風(fēng)輪機這塊不要只填一個風(fēng)速常數(shù)很多簡化的模型喜歡把風(fēng)輪機固定成恒功率源這種做法在電網(wǎng)側(cè)仿真里夠用但做整機控制就露餡了。真正要體現(xiàn)直驅(qū)系統(tǒng)的響應(yīng)特性風(fēng)輪機模型必須有點風(fēng)的動態(tài)。我一般推薦在基速以下用最大功率追蹤MPPT的思路建模。核心關(guān)系是風(fēng)能利用系數(shù) Cp 和葉尖速比 λ 之間的關(guān)系葉尖速比公式λ ωr·R / Vw其中 ωr 是風(fēng)輪角速度R 是風(fēng)輪半徑Vw 是來流風(fēng)速。機械功率Pm 0.5·ρ·π·R2·Vw3·Cp(λ, β)ρ 是空氣密度β 是槳距角。在Simulink里搭這個部分時建議直接用MATLAB Function或者Fcn模塊寫 Cp 表達(dá)式不要用查表法糊弄。這樣后續(xù)改風(fēng)剪效應(yīng)、塔影效應(yīng)或者陣風(fēng)模型都很方便換一組風(fēng)速序列就能變化出不同的運行工況。2.2 永磁同步電機的參數(shù)別照抄別的電機直驅(qū)風(fēng)機的PMSG和普通伺服電機建模不太一樣核心參數(shù)是極對數(shù) p、永磁磁鏈 ψf、d軸和q軸電感 Ld、Lq以及定子電阻 Rs。在Simulink的Permanent Magnet Synchronous Machine模塊里你可以選擇SI單位或者標(biāo)幺值pu單位。這里有個實際的建議如果想仿真結(jié)果在分析時比較直觀用SI單位國際單位制填就行這樣你要觀察轉(zhuǎn)速rad/s、電磁轉(zhuǎn)矩N·m、電流A都是真實值不用來回?fù)Q算。如果做快速控制原型或者硬件在環(huán)再考慮用pu進(jìn)行歸一化。2.3 變流器的IGBT模型和開關(guān)頻率機側(cè)和網(wǎng)側(cè)變流器的實現(xiàn)方式直接決定了仿真精度和速度。如果你只關(guān)心控制邏輯用平均值模型Average-Model就夠——速度極快沒有開關(guān)紋波適合長時間仿真和策略驗證。如果你要觀察電流諧波、電壓紋波甚至做并網(wǎng)電能質(zhì)量分析就得用IGBT模型。需要注意一點IGBT模型的開關(guān)頻率會直接限制仿真步長。比如開關(guān)頻率10kHz仿真步長至少要在1e-5秒量級才能看到比較合理的PWM波形。若步長設(shè)到1e-4秒PWM狀態(tài)變化大概率被跳過波形會非常難看而且容易莫名其妙發(fā)散。3. 380V和690V兩套模型不只是改個電壓數(shù)字的事拿到380V和690V兩個仿真這個需求時最忌諱的做法就是生成兩個模型然后把一個的電壓常量從380改到690其他什么都不動。這兩套系統(tǒng)背后其實對應(yīng)著完全不同的功率等級和應(yīng)用場景。3.1 先搞清楚這兩套系統(tǒng)分別代表什么380V的系統(tǒng)在我的理解里更偏向?qū)嶒炇覙訖C或者小功率測試臺比如幾kW到幾十kW的直驅(qū)模擬平臺。690V的系統(tǒng)則是兆瓦級風(fēng)電機組常見的定子側(cè)交流電壓。實際工程中690V是低壓三電平或者多電平變流器的常見水平再往上就是中壓1140V或者10kV/35kV并網(wǎng)了。所以同一個課題380V用來驗證控制算法690V用來做整機并網(wǎng)性能評估場景完全不同。3.2 參數(shù)縮放的對應(yīng)關(guān)系項目380V系統(tǒng)690V系統(tǒng)典型功率等級10kW~50kW1MW~3MW額定線電壓380V690V直流母線電壓約620V~650V約1100V~1200V永磁磁鏈小按低電壓設(shè)計大按額定并網(wǎng)設(shè)計定子相電流幾十A幾百到上千A濾波電感較大標(biāo)幺值高較小標(biāo)幺值按工程標(biāo)準(zhǔn)電壓等級確定之后首先受影響的是變流器調(diào)制比。直流母線電壓至少要高于交流側(cè)線電壓峰值SVPWM才能輸出所期望的交流電壓。380V系統(tǒng)的交流線電壓峰值大約是537V考慮到開關(guān)壓降和阻抗壓降母線定在650V給一定裕量690V系統(tǒng)的線電壓峰值約976V母線一般會取1100V到1200V。如果母線電壓不夠高逆變器會進(jìn)入過調(diào)制區(qū)電流波形失真嚴(yán)重功率也送不出去。3.3 控制器的PI參數(shù)跟著一起換這個是最多人踩的坑。同一套PI參數(shù)在380V的10kW樣機上調(diào)得好好的搬到690V的兆瓦級模型上直接震蕩發(fā)散。原因是電流環(huán)和電壓環(huán)的增益都和電壓、功率等級相關(guān)。我在實際調(diào)試中一般先做電流環(huán)的內(nèi)環(huán)設(shè)計按照電機電感和電阻計算時間常數(shù)再用零極點對消法整定。額定電壓變了額定電流變了電感阻抗的標(biāo)幺值也變了PI參數(shù)必然要重新算。網(wǎng)側(cè)電壓環(huán)的時間常數(shù)還要比電流環(huán)慢5~10倍保證內(nèi)外環(huán)時間尺度拉開不然整個系統(tǒng)會耦合在一起動態(tài)響應(yīng)一團(tuán)亂。所以拿到模型后在改電壓等級的同時一定要同步檢查電流環(huán)的比例積分增益是否按照新系統(tǒng)的阻抗參數(shù)重算電壓外環(huán)帶寬是否比電流環(huán)低一個數(shù)量級功率外環(huán)如果有MPPT轉(zhuǎn)速環(huán)是否能匹配新系統(tǒng)的機械時間常數(shù)4. 仿真運行流程與波形驗證這樣跑才能確認(rèn)模型是對的4.1 從初始化到穩(wěn)態(tài)我的標(biāo)準(zhǔn)流程第一件事建議把求解器設(shè)置成定步長固定步長步長取1e-5或者更小。連續(xù)狀態(tài)離散PWM混合系統(tǒng)用變步長求解器在開關(guān)動作點上會遇到很多不連續(xù)問題經(jīng)常出現(xiàn)仿真速度時快時慢甚至卡死的情況。定步長雖然耗時一點但勝在穩(wěn)定結(jié)果可控。然后按照這個順序遞進(jìn)跑空載啟動強制直流母線預(yù)充電到額定電壓機側(cè)和網(wǎng)側(cè)變流器全都封鎖脈沖檢查母線電壓能不能穩(wěn)定。網(wǎng)側(cè)掛負(fù)載/并網(wǎng)放開網(wǎng)側(cè)逆變器控制讓直流母線電壓閉環(huán)穩(wěn)定觀察母線電壓是否紋波大、是否有低頻振蕩。機側(cè)啟動發(fā)電給風(fēng)輪機一個恒定風(fēng)速機側(cè)變流器開始整流觀察發(fā)電機轉(zhuǎn)速、電磁轉(zhuǎn)矩、功率輸出是否平滑上升。整機運行風(fēng)速加入階躍或者斜坡變化觀察各個變量能否平滑過渡是否出現(xiàn)超調(diào)或者震蕩。極限工況額定風(fēng)速以上突降風(fēng)速、負(fù)載切換看看系統(tǒng)的過流、過壓保護(hù)邏輯是否起作用。4.2 最該盯的幾個波形仿真跑起來之后別只看功率或者直流母線電壓這些量平順不代表控制系統(tǒng)沒問題。我更建議盯這幾個發(fā)電機電流波形是否正弦、是否有明顯畸變、電流環(huán)調(diào)節(jié)是否有靜差。直流母線電壓是否存在6倍頻脈動機側(cè)不控整流常見的諧波來源在負(fù)載突變時電壓下墜幅度是否在允許范圍內(nèi)。網(wǎng)側(cè)三相電流是否和電網(wǎng)電壓保持正確的相位關(guān)系有功無功解耦是否干凈。風(fēng)輪轉(zhuǎn)速與風(fēng)能利用系數(shù)這一步最能驗證MPPT有沒有起作用Cp能不能被穩(wěn)定控制在最大值附近。5. 仿真翻車現(xiàn)場發(fā)散、跑不動、波形畸變的排查鏈路5.1 模型莫名發(fā)散問題多半在求解器我在給模型調(diào)參時遇到的仿真發(fā)散十次里有六次不是控制參數(shù)的問題而是求解器配置不對。Simulink默認(rèn)的變步長ode45算法求解剛性系統(tǒng)時可能效率很低甚至反復(fù)步進(jìn)失敗。如果遇到仿真一步就開始發(fā)散的情況排查順序是把所有PWM脈沖封掉只留直流母線電容和電源看能否收斂。封掉交流源只留整流器和電機空載跑看能否建立電壓。逐步解鎖控制環(huán)節(jié)每解開一層就觀察一層波形。如果每一步都沒有發(fā)散但合在一起就發(fā)散多數(shù)是因為某個閉環(huán)的反饋極性反了。PID控制器輸出負(fù)反饋變成正反饋系統(tǒng)自然瞬間爆炸。這種問題在Simulink里可以用示波器跟變量曲線判斷比如觀察控制器的輸出指令是否往反方向越調(diào)越大。5.2 電流波形畸變PWM扇區(qū)和死區(qū)相關(guān)的坑如果你用SVPWM模塊波形畸變先檢查扇區(qū)判斷邏輯和調(diào)制比是否超限。SVPWM的線性調(diào)制范圍最大是母線電壓除以根號3超過這個值就進(jìn)入過調(diào)制區(qū)波形會嚴(yán)重失真。另外一個容易被忽視的點是死區(qū)時間。IGBT模塊里的死區(qū)時間設(shè)置為0時理論上波形最干凈但實際轉(zhuǎn)換器模型不會因為死區(qū)產(chǎn)生明顯的問題可一旦死區(qū)時間設(shè)置和開關(guān)頻率不匹配比如10kHz的載波死區(qū)設(shè)成10微秒波形就會因為有效占空比損失而畸變。仿真的目的是驗證控制邏輯死區(qū)時間設(shè)小一點比如1微秒是合理的操作既保證現(xiàn)象符合實際工程規(guī)律又不會讓波形難看。5.3 仿真跑得太慢或者直接卡死直驅(qū)風(fēng)機整體模型包含風(fēng)輪機、PMSG、兩個變流器和控制系統(tǒng)加上IGBT這種帶高頻開關(guān)的模塊仿真速度確實是個真問題。如果你的模型一個幾十秒的風(fēng)速變化工況要跑一個晚上肯定是配置有問題。我實測比較好用的優(yōu)化組合是用離散電力電子模型不要用連續(xù)狀態(tài)。在電力系統(tǒng)庫中選擇離散化版本比如用Discrete模式的變流器可以顯著提高運算速度。把開關(guān)頻率降低到實際應(yīng)用的下限。對于驗證MPPT和并網(wǎng)策略來說5kHz和10kHz在波形形狀上的差異遠(yuǎn)沒有你想象的大但計算量差距巨大。風(fēng)速和機械參數(shù)用連續(xù)域電氣部分用離散域混合求解。這是Simulink里性能最優(yōu)的一種做法需要你對模型的代穩(wěn)定性有一定控制能力。如果模型實在太大可以考慮在參數(shù)初始化時用matlab腳本批量賦值避免每次仿真都在命令行堆一大堆set_param操作。5.4 690V系統(tǒng)特有的收斂問題690V系統(tǒng)因為功率等級大變頻器模塊里的開關(guān)瞬態(tài)電流變化率很高仿真中更容易出現(xiàn)數(shù)值振蕩表現(xiàn)為波形末尾出現(xiàn)高頻毛刺。這種毛刺在低功率的380V模型上幾乎看不到但在690V模型上很常見。原因是電感電流變化率大數(shù)值積分誤差被放大。應(yīng)對辦法很簡單把仿真步長再降低一檔比如從1e-5降到5e-6或2e-6同時確認(rèn)開關(guān)器件模塊的默認(rèn)阻抗參數(shù)沒有被改動過。部分情況下在變流器的理想開關(guān)模型里增加一個很小的緩沖電阻Snubber Resistance也能顯著收斂數(shù)值振蕩。最后再分享一個我自己固定下來的習(xí)慣每次修改模型參數(shù)前先把仿真結(jié)果的穩(wěn)態(tài)值記錄在Excel或者.mat文件里包括母線電壓、輸出功率、電流幅值。改完參數(shù)后對照這些基準(zhǔn)值能快速發(fā)現(xiàn)參數(shù)改動導(dǎo)致的雪崩式漂移。尤其是380V和690V兩套模型幾套參數(shù)來回切換的時候這樣管理參數(shù)才不會把兩套系統(tǒng)的結(jié)果搞混。永磁直驅(qū)風(fēng)力發(fā)電這個方向模型跑通只是第一步能把兩套電壓等級下的行為差異講清楚才算真的把這個仿真吃透了。本文還有配套的精品資源點擊獲取