趙曉龍 陳卓 李強勝 何俊賢

摘? 要： 針對實際電力系統(tǒng)運行和控制的復雜性與Matlab仿真非實時性的矛盾，提出跨平臺聯(lián)合仿真。采用基于Matlab/Simulink代碼生成技術與RTDS實時仿真系統(tǒng)相結合的方法，以高壓直流輸電系統(tǒng)模型為研究對象，利用Embedded Coder開發(fā)工具對系統(tǒng)的控制模塊進行代碼生成，并在RTDS中完成代碼移植實驗。通過對比仿真波形，驗證得到結果行為一致。由此表明，該方法能夠有效利用Simulink良好的開發(fā)環(huán)境和RTDS的實時仿真特性，在離線仿真成功的基礎上，完成實時仿真測試，使工作便捷高效，同時保證結果精準可信。

關鍵詞： 代碼生成技術; 跨平臺聯(lián)合仿真; Matlab/Simulink; 實時仿真測試; 高壓直流輸電系統(tǒng)建模; 模型參數計算; 控制模塊代碼移植; 仿真波形對比

中圖分類號： TN710?34; TM721.1? ? ? ? ? ? ? ? ? ? 文獻標識碼： A? ? ? ? ? ? ? ? ? 文章編號： 1004?373X（2020）09?0161?05

Application of Matlab/Simulink?based code generation technology in RTDS

ZHAO Xiaolong， CHEN Zhuo， LI Qiangsheng， HE Junxian

（College of Electrical Engineering， Guizhou University， Guiyang 550025， China）

Abstract： In view of the contradiction between the complexity of actual power system operation and control and the non?real?time characteristic of Matlab simulation， the cross?platform co?simulation is proposed. A method on the basis of the combination of Matlab/Simulink?based code generation technology and RTDS?based （real time digital simulator） real?time simulation system is adopted. By taking the high voltage DC transmission system model as the research object， the development tool Embedded Coder is used to generate code for the control module of the system and complete the code transplantation experiment in RTDS. It is verified that the result and the behavior are consistent by contrasting the simulation waveforms. Therefore， it shows that the proposed method can effectively utilize the excellent development environment of Simulink and the real?time simulation characteristics of RTDS to complete real?time simulation test on the basis of successful off?line simulation， which makes the work convenient and efficient and ensures the accuracy and credibility of the results.

Keywords： code generation technology; cross?platform co?simulation; Matlab/Simulink; real?time simulation test; high voltage DC transmission system modeling; model parameter calculation; control module code transplantation; simulation waveform contrast

0? 引? 言

Matlab/Simulink是一款功能強大的高效計算機仿真軟件，但其仿真大多為非實時仿真，在對于仿真實時性要求較高的任務中，單純的Matlab/Simulink仿真難以勝任。如文獻[1?2]中利用軟硬件協(xié)同控制，文獻[3?4]中采取重新建模的方式，都是可行的方法，但是缺點也很明顯，其中有的工程化不完全，有的步驟繁瑣，內容重復，并不是最佳方案。

在對仿真實時性要求高或者需要跨平臺硬件測試的情況下，使用Simulink Coder或者專門針對嵌入式系統(tǒng)的Embedded Coder工具箱，將Simulink模型轉換為C/C++語言模型，即通常所說的C代碼，根據目標配置能夠生成不同環(huán)境下的程序文件，載入到RTDS等平臺進行硬件擴展，可以與測試對象進行實時互動，由此實現動態(tài)系統(tǒng)的硬件調試，得到高效實時的仿真結果。

本文以高壓直流輸電模型為研究對象，進行Simulink建模與代碼轉換，并在RTDS上進行仿真測試，驗證C代碼跨平臺實驗的可行性。

1? 模型分析

1.1? 高壓直流輸電系統(tǒng)建模

本文以國際大電網會議高壓直流輸電標準測試系統(tǒng)（CIGRE HVDC標準測試系統(tǒng)）為研究對象[5]，如圖1所示。

此高壓直流輸電系統(tǒng)采用單極輸電，系統(tǒng)的構成包括交流系統(tǒng)、直流系統(tǒng)、等效阻抗、濾波器、電容器、變壓器和換流器等元件。整流側和逆變側都使用12脈波換流器，并且采用一組高通濾波器和一組雙調諧濾波器。為了使輸電模型更加簡潔直觀，直流輸電線路則由等效的T型電路代替。

1.2? 模型的控制策略

整流側采用定電流控制，取測量的實際電流值與整定電流值的偏差進行PI控制。其中，整定電流值是由低壓限流裝置VDCOL和主控層的電流調制控制器共同決定的。PI控制環(huán)節(jié)進行計算后，輸出控制閥所需的延遲觸發(fā)角指令，閥控裝置根據這個指令生成相應的控制脈沖，控制框圖如圖2所示。

為使直流電流能夠流通，整流側的直流電壓需要高于逆變側的直流電壓，同時，整流器的定電流控制也具備抑制直流電流增加的作用。

逆變側采用定關斷角控制為主，定電流控制為輔的控制方式。從直流系統(tǒng)的逆變側換流器測得上一周期關斷角的最小值[γ]，然后將逆變側測得的直流電流實際值與整定值之差通過電流偏差控制，得到電流偏差角度[Δγ]。電流偏差控制能夠使逆變側的控制方式在定關斷角控制和定電流控制之間進行平滑過渡。CIGRE標準一般規(guī)定直流系統(tǒng)的整定[γ]角為[15°]，最大偏差限制為[-31°]。通過PI控制后，輸出逆變側延遲觸發(fā)角指令，控制框圖如圖3所示。

通過這種控制方式，其目的一方面是為了防止逆變器換相失敗，同時又保證無功需求最小。通過對實際關斷角和設定值差的補正量，確保所需最小的關斷角。補正量會讓關斷角向增大的方向移動。

逆變側定電流控制與整流側原理類似，區(qū)別在于逆變側整定電流值比整流側整定電流值要小一個電流裕度，通常取為0.1 p.u.，所以整定值比實際電流小，一般作為輔助控制手段。

1.3? 模型的參數計算

在CIGRE標準中，濾波器、電容器和等效阻抗等數據已經給出，只需計算出換流器的參數。由于每個12脈波換流器由兩個6脈波換流器組成，所以采用雙繞組變壓器，結合CIGRE標準給出的系統(tǒng)參數，可以得到整流側電路數據，并計算出每個變壓器的容量[Tmva]和正序電抗[XL]（單位：p.u.）：

[Vo=V?3?2π] （1）

[Vdc=Vo?cos α-3ωLπ?Idc] （2）

[Tmva=Pdc?π3cos α-3ωLπ?IdcVo] （3）

[Zbase=V2Tmva] （4）

[XL（p.u.）=ωLZbase] （5）

式中：[Vdc]，[Idc]和[Pdc]分別為整流側的直流電壓、直流電流和直流功率（電壓和功率只有實際值的一半）;[V]為變壓器二次側電壓;[α]為觸發(fā)角。

晶閘管的緩沖電路電阻和電容取決于額定直流電壓，RC時間常數至少是仿真步長的兩倍，當仿真步長為50 μs時，需要[Rs?Cs>100 μs]。最后，利用CIGRE標準提供參考的PI參數，結合文獻[6?8]完成PI參數的優(yōu)化。

2? 代碼轉換

2.1? 代碼生成流程

使用Build Model將已經通過離線仿真檢驗的CIGRE HVDC標準測試模型轉換為rtw文件，然后通過Simulink Coder中的目標語言編譯器將rtw文件轉換為一系列的源文件，這個過程中TLC所使用的文件包括3類：系統(tǒng)目標文件、模塊的目標文件和支持代碼生成的TLC函數庫等文件。模型的源代碼全部生成之后，就可以導入目標硬件中進行實機運行。代碼生成流程如圖4所示。

2.2? 直流輸電模型代碼生成的配置方法

在進行代碼生成時，需要對Simulink模型進行參數的設置，在Configuration Parameter中集中管理著Simulink模型的代碼生成方法、格式等約束條件。為了生成能夠在RSCAD中正確運行的代碼，需要修改以下參數。

1） 在solver頁面的solver options一欄中，type參數設置為Fixed?step，選擇ode3（Bogacki?Shampine）解法器，同時，設置Fixed?step size為50 μs，這是因為代碼生成支持的是固定步長的模型，而且控制模塊計算只需用到常微分方程。

2） 在Code Generation頁面的Target selection一欄，System target file參數配置為選擇生成ert.tlc文件，這是Embedded Coder提供的系統(tǒng)目標文件，針對嵌入式進行性能與內存的優(yōu)化，同時保證結構簡潔，方便調用。

3） 由于模型存在連續(xù)時間函數模塊，但是嵌入式器件的運行是由固定時鐘進行驅動，必須為離散模塊，因此，需要在Code Generation的Interface頁面中的Software environment一欄勾選continuous time選項，確保代碼成功生成。

通過正確配置，完成高壓直流輸電控制模塊嵌入式C代碼的生成。同時可以選擇自動生成報告（Create Code Generation Report），Simulink自帶注釋（Include comments），參數內聯(lián)（Inline parameters）等，使得代碼易于修改，并增強可讀性。

最后，為了提高生成的C代碼的可靠性，需要進行軟件在環(huán)（SIL）測試，即在Code Generation的Verification頁面，將Create block選為SIL，用于生成控制模塊的[S]函數。由于SIL測試時要求模型中不能有連續(xù)時間函數，所以需要將采樣時間設置成離散，同時由于模型中有非內聯(lián)函數，需要在Code Generation?>Interface頁面勾選non?Inlined S?functions項，保證成功生成[S]函數。最后給予[S]函數和原控制模型相同的輸入，并將兩組輸出結果做減法運算，驗證代碼與模型的行為一致性。

3? 硬件驗證

3.1? 實時仿真機RTDS

實時仿真機（Real Time Digital Simulator，RTDS）是一種專門用于研究電力系統(tǒng)電磁暫態(tài)現象的裝置。其算法繼承了H.W.Dommel建立的經典電磁暫態(tài)計算理論，將那些復雜的電力網絡進行等值變化，能夠獲得等值計算網絡，然后對其進行方程求解，獲得所需要的結果。由于它的電力系統(tǒng)運算法則能被計算得足夠快，能產生連續(xù)的輸出結果，這些輸出代表在實際網絡中的情形[9]。

RTDS基本組成單元稱為Rack，裝置由IRC，GPC，WIF等多種板卡和處理器組成，一個CPU模擬一個電力系統(tǒng)元件，CPU之間采用并行?串行?并行的通信方式。本次RTDS高壓直流輸電仿真模型使用了1個RACK，共計6塊PB5處理器，將4塊PB5分給控制模塊用于快速計算，其余2塊用于模擬交流電源、變壓器和換流器等器件模型。

3.2? RSCAD仿真

RSCAD軟件是RTDS自帶的專業(yè)仿真軟件，作為電磁暫態(tài)仿真的一種，通常會將仿真步長設置為50 μs。在進行高壓直流輸電代碼移植的仿真實驗時，需要在RSCAD的Draft模塊搭建高壓直流輸電的基本電路模型，通過CBuilder中的Tools/Import Simulink Component工具，將Simulink生成的嵌入式C代碼移植到RSCAD中，形成用戶自定義模塊[10?12]，如圖5所示。

在C File Associations頁面，進入Edit中對代碼的接口名稱進行修改，方便用戶識別，此時還能夠添加用戶自行編程的代碼，增強自定義模塊的功能。另外，在Parameters和Graphics頁面，以表格化的固定形式編輯模型，定義模型數據格式與可變參數，方便在模型使用時進行調節(jié)與管理。

最后在進行Draft建模時，通過器件庫的Add Component功能添加用戶自定義模塊到電路中去，這樣就完成控制模塊從Simulink模型到RTDS模型的移植。值得注意的是，RTDS中電壓與電流的單位都是kV和kA。

完成Draft設計之后，對模型進行編譯，通過compile生成相關文件。其中.sib文件集顯示仿真結果與調節(jié)相關參數為一體，打開.sib文件，然后關聯(lián)所需信號與對應的顯示面板。同時，還要設置關聯(lián)控制按鈕和滑輪，這樣可以在實時仿真中隨時調節(jié)所關聯(lián)的電壓參數，以及控制模型整流和逆變的啟動和停止。

3.3? 仿真結果分析

CIGRE高壓直流輸電標準測試系統(tǒng)的額定直流電流為2 kA，額定電壓為500 kV，整流側交流電壓為345 kV，逆變側交流電壓為230 kV。圖6和圖7為模型在Simulink中仿真的直流電流和電壓波形，在主控制器的調控下，系統(tǒng)根據設定在[0.02 s]啟動，電流與電壓能夠按照額定的速率平緩地上升到額定值，并保持穩(wěn)定運行。

為了突出控制模塊的作用，對比Simulink仿真與RTDS仿真啟動時無主控的情況，現除去主控制器，由選定的控制策略單獨控制，圖8和圖9為Simulink仿真的直流電流和電壓波形圖，可以看出，在啟動的極短時間內，電流和電壓都會有較大的波動，隨后在控制模塊的控制下，模型會進入并維持在穩(wěn)定狀態(tài)。

在RTDS中仿真的波形如圖10和圖11所示，電流和電壓的波形與Simulink仿真中的基本一致，只是在啟動時的波動有偏差，這是由兩個平臺使用的脈波發(fā)生器的差異以及仿真算法的不同導致的。另外，RTDS仿真由于多CPU并行計算的速度足夠快，因此產生了連續(xù)的輸出結果，基本能夠代表實際運行網絡中的情形，這也是實時仿真的優(yōu)勢，并且RTDS仿真在運行過程中能夠方便地調節(jié)系統(tǒng)基本參數，以此展現出動態(tài)波形的變化，方便對系統(tǒng)整體運行的分析與研究，這也是單純的Simulink仿真無法比擬的。

通過對比發(fā)現，CIGRE高壓直流輸電系統(tǒng)在相同的控制模塊下，于Matlab平臺和RTDS平臺上都能夠持續(xù)穩(wěn)定地運行，證明了利用代碼生成技術將系統(tǒng)的控制模塊向RTDS上的移植是可行的。

4? 結? 語

通過Simulink的代碼移植，實現跨平臺仿真，由于Matlab應用的廣泛性，從而擁有豐富的Simulink模型庫和大量的實例。另外，Matlab的離線性可以方便工作人員進行仿真設計，在完成離線驗證之后，再移植到對應的平臺，省去重復查找錯誤的步驟，提高工作效率。而且C代碼具備通用性，不僅限于在RTDS上的移植，還可以應用于多個平臺[13?14]，針對不同平臺硬件的差異進行編譯，從而實現多平臺互聯(lián)，這使得代碼生成技術具有很強的實用性。

參考文獻

[1] 錢振天，徐曉軼，諶平平，等.基于Matlab/Simulink的軟件在環(huán)仿真技術研究[J].電力電子技術，2016，50（10）：5?7.

[2] LAI C K， TSAI C C. The modeling and simulation for current regulation of PMSM drive system using Verilog HDL and variable structure control [C]// 2016 International Automatic Control Conference. Taichung， China： IEEE， 2016： 218?219.

[3] 關紅兵，黃立濱，周鵬鵬，等.基于RTDS的高壓直流控制系統(tǒng)鎖相環(huán)仿真研究[J].電力系統(tǒng)保護與控制，2017，45（12）：128?132.

[4] 武霄，崔勇，黃志龍，等.基于RTDS的高壓直流輸電控制功能研究[J].華北電力大學學報（自然科學版），2014，41（6）：60?66.

[5] FARUQUE M O， ZHANG Y， DINAVAHI V. Detailed mo?deling of CIGRE HVDC benchmark system using PSCAD/EMTDC and PSB/SIMULINK [J]. IEEE transactions on power deli?very， 2005， 21（1）： 378?387.

[6] 李天云，劉智銘，黨震宇，等.基于Simulink的HVDC PI控制器參數優(yōu)化方法[J].電力自動化設備，2012，32（1）：45?48.

[7] 于飛，鐘秋海，張望，等.基于Matlab的超高壓直流輸電系統(tǒng)仿真研究[J].繼電器，2004，32（13）：31?34.

[8] KWON D H， MOON H J， KIM R G， et al. Modeling of CIGRE benchmark HVDC system using PSS/E compared with PSCAD [C]// International Youth Conference on Energy. Pisa， Italy： IEEE， 2015： 1?8.

[9] 翟凱.基于RTDS的雙單元背靠背直流輸電系統(tǒng)仿真與研究[D].北京：華北電力大學，2017.

[10] 徐志恒，錢珞江.RTDS用戶自定義元件模型（UDC）的應用方法[J].電力系統(tǒng)保護與控制，2009，37（22）：141?145.

[11] 楊洋，肖湘寧，陶順，等.多端口電磁?機電暫態(tài)混合實時仿真接口技術與應用[J].電工技術學報，2017，32（4）：231?240.

[12] 楊柳林，李德奎，陳延明，等.基于RTDS的電弧接地故障自定義建模及仿真分析[J].電力系統(tǒng)保護與控制，2016，44（16）：137?142.

[13] 劉濤，蔡燕，張海華，等.基于代碼生成技術的Matlab/Simulink實驗教學改革[J].實驗技術與管理，2018，35（8）：157?159.

[14] 王夢茹，周珊，王金波，等.基于Matlab的智能測試系統(tǒng)PMSM物理實時仿真平臺[J].系統(tǒng)仿真學報，2018，30（6）：2225?2235.