基于風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量補(bǔ)償策略的濾波參數(shù)選取方法
【專(zhuān)利摘要】本發(fā)明公開(kāi)了一種考慮通信時(shí)滯的風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量補(bǔ)償策略的濾波參數(shù)選取方法，首先確定風(fēng)力機(jī)模擬器系統(tǒng)的基本參數(shù)，之后確定一階數(shù)字濾波器的濾波參數(shù)；并判斷此時(shí)模擬器系統(tǒng)是否滿(mǎn)足穩(wěn)定性條件；然后根據(jù)判斷的結(jié)果對(duì)一階數(shù)字濾波器的濾波參數(shù)進(jìn)行自增；濾波參數(shù)選取結(jié)束后，轉(zhuǎn)動(dòng)慣量補(bǔ)償有效，風(fēng)力機(jī)模擬器系統(tǒng)能夠穩(wěn)定運(yùn)行。本發(fā)明解決了含通信時(shí)滯的風(fēng)力機(jī)模擬器系統(tǒng)中轉(zhuǎn)動(dòng)慣量補(bǔ)償策略的一階低通濾波器參數(shù)選取問(wèn)題，通過(guò)判斷含通信時(shí)滯的風(fēng)力機(jī)模擬器的系統(tǒng)穩(wěn)定性來(lái)確定風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量補(bǔ)償策略中一階低通濾波器的濾波參數(shù)選取，使WTS能安全穩(wěn)定地模擬大轉(zhuǎn)動(dòng)慣量風(fēng)力機(jī)的機(jī)械動(dòng)態(tài)。
【專(zhuān)利說(shuō)明】
基于風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量補(bǔ)償策略的濾波參數(shù)選取方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明屬于風(fēng)力機(jī)模擬器領(lǐng)域，特別是一種基于風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量補(bǔ)償策略 中低通濾波器參數(shù)的選取方法。
【背景技術(shù)】
[0002] 開(kāi)展風(fēng)力機(jī)試驗(yàn)的目的將助力風(fēng)輪葉片氣動(dòng)設(shè)計(jì)、風(fēng)力機(jī)伺服控制、電氣變流技 術(shù)等的深入研究，對(duì)風(fēng)力發(fā)電系統(tǒng)的開(kāi)發(fā)、優(yōu)化和現(xiàn)場(chǎng)應(yīng)用具有重要的意義。
[0003] 在實(shí)地開(kāi)展風(fēng)力機(jī)試驗(yàn)的代價(jià)十分高昂且耗費(fèi)時(shí)間，故風(fēng)力機(jī)模擬器（Wind Turbine Simulator，WTS)已經(jīng)成為風(fēng)力機(jī)試驗(yàn)的重要測(cè)試設(shè)備。
[0004] 由于風(fēng)力機(jī)模擬器的物理轉(zhuǎn)動(dòng)慣量遠(yuǎn)小于實(shí)際風(fēng)力機(jī)的轉(zhuǎn)動(dòng)慣量，模擬器須采取 轉(zhuǎn)動(dòng)慣量補(bǔ)償策略以還原實(shí)際風(fēng)力機(jī)的慢機(jī)械動(dòng)態(tài)特性。傳統(tǒng)轉(zhuǎn)動(dòng)慣量補(bǔ)償策略中一階低 通濾波器的參數(shù)選擇對(duì)于模擬性能的優(yōu)劣具有很大影響。
[0005] 當(dāng)模擬大轉(zhuǎn)動(dòng)慣量風(fēng)力機(jī)時(shí)，應(yīng)用一階濾波器的轉(zhuǎn)動(dòng)慣量補(bǔ)償策略由于通信時(shí)滯 而存在失效現(xiàn)象，濾波器參數(shù)選取不合適一方面會(huì)影響模擬效果(參數(shù)選取過(guò)大），另一方 面可能易造成系統(tǒng)失穩(wěn)崩潰(參數(shù)選取過(guò)小）。
[0006] 基于上述情況，目前迫切需要一種風(fēng)輪模擬器轉(zhuǎn)動(dòng)慣量補(bǔ)償策略中濾波器的參數(shù) 選取方法，使得風(fēng)力機(jī)模擬器能夠穩(wěn)定可靠地模擬大轉(zhuǎn)動(dòng)慣量風(fēng)機(jī)的機(jī)械動(dòng)態(tài)。

【發(fā)明內(nèi)容】

[0007] 本發(fā)明的目的在于提供一種基于風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量補(bǔ)償策略的濾波參數(shù)選 取方法。
[0008] 實(shí)現(xiàn)本發(fā)明目的的技術(shù)解決方案為:一種基于風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量補(bǔ)償策略的 濾波參數(shù)選取方法，包含以下步驟：
[0009] 步驟1、確定風(fēng)力機(jī)模擬器系統(tǒng)的基本參數(shù)，包括通信時(shí)滯步長(zhǎng)kQ，轉(zhuǎn)動(dòng)慣量補(bǔ)償 倍數(shù)Jt/J s，其中Jt為風(fēng)力機(jī)轉(zhuǎn)動(dòng)慣量，Js為風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量;時(shí)滯步長(zhǎng)ko的確定公式 為：
[0010] k〇=[i/T]
[0011] 式中，τ是通信時(shí)滯時(shí)長(zhǎng)，T是系統(tǒng)采樣時(shí)長(zhǎng)。
[0012] 步驟2、確定一階數(shù)字濾波器的濾波參數(shù)α;所述一階數(shù)字濾波器的濾波參數(shù)α的初 值為0.5。
[0013] 步驟3、根據(jù)步驟1和步驟2確定的參數(shù)，判斷此時(shí)模擬器系統(tǒng)是否滿(mǎn)足穩(wěn)定性條 件;如果滿(mǎn)足則執(zhí)行步驟5，否則執(zhí)行步驟4;判斷模擬器系統(tǒng)是否滿(mǎn)足穩(wěn)定性條件具體包括 以下步驟：
[0014] 步驟3-1、構(gòu)建方程組
[0015]
[0016] 式中，r是極徑，Θ是極角；
[0017] 步驟3-2、對(duì)方程組進(jìn)行初始化，具體是:將極角Θ設(shè)置為〇;判斷此時(shí)時(shí)滯步長(zhǎng)ko是 否為奇數(shù)，若為奇數(shù)，則執(zhí)行步驟3-3;否則執(zhí)行步驟3-4;
[0018] 步驟3-3、判斷Jt/JS〈2/(l-a)是否滿(mǎn)足，若滿(mǎn)足則執(zhí)行步驟3-4,若不滿(mǎn)足則執(zhí)行 步驟3-9;
[0019] 步驟3-4、極角Θ進(jìn)行自增，具體為:θ = θ+ Δ θ，Δ Θ取值為〇 · 〇1;
[0020] 步驟3-5、根據(jù)下述公式計(jì)算極徑r，
[0021]
[0022] 步驟3-6、對(duì)r值進(jìn)行校驗(yàn)，校驗(yàn)公式為：
[0023]
[0024] 式中ε是校驗(yàn)精度，取值為IX ΚΓ3;若校驗(yàn)通過(guò)，則執(zhí)行下一步，若不通過(guò)返回執(zhí)行 步驟3-4;
[0025] 步驟3-7、判斷極徑r是否小于1，若是則執(zhí)行步驟3-8,否則直接執(zhí)行步驟3-9;
[0026]步驟3-8、判斷極角Θ是否大于231，若大于231，則風(fēng)力機(jī)模擬器系統(tǒng)的穩(wěn)定性條件成 立，否則返回執(zhí)行步驟3-4;
[0027]步驟3-9、風(fēng)力機(jī)模擬器系統(tǒng)的穩(wěn)定性條件不成立。
[0028] 步驟4、對(duì)一階數(shù)字濾波器的濾波參數(shù)α進(jìn)行自增，具體為:α = α+ Δ α，Δ α取〇 . 〇1， 之后返回步驟3;
[0029] 步驟5、濾波參數(shù)選取結(jié)束，轉(zhuǎn)動(dòng)慣量補(bǔ)償有效，風(fēng)力機(jī)模擬器系統(tǒng)能夠穩(wěn)定運(yùn)行。
[0030] 本發(fā)明與現(xiàn)有技術(shù)相比，其顯著優(yōu)點(diǎn)為：1)本發(fā)明首次提出了一個(gè)考慮通信時(shí)滯 的一階濾波器參數(shù)選取方法，解決了因一階濾波器參數(shù)選取不當(dāng)而造成的模擬器失穩(wěn)問(wèn) 題;2)本發(fā)明公開(kāi)了一階濾波器參數(shù)選取的詳細(xì)步驟，簡(jiǎn)單易行且能夠保證風(fēng)力機(jī)模擬器 能夠穩(wěn)定運(yùn)行，使WTS能模擬更大轉(zhuǎn)動(dòng)慣量的風(fēng)力機(jī)模型，并且能更加有效和方便地復(fù)現(xiàn)實(shí) 際風(fēng)力機(jī)的慢動(dòng)態(tài)。
【附圖說(shuō)明】
[0031] 圖1為本發(fā)明的風(fēng)力機(jī)模擬器中一階低通濾波器參數(shù)選取方法的流程圖。
[0032] 圖2為本發(fā)明的具體實(shí)施例的模擬器結(jié)構(gòu)圖及濾波器設(shè)置示意圖。
[0033] 圖3為對(duì)比例I的軌跡圖，其中圖（a)是應(yīng)用傳統(tǒng)策略的模擬器轉(zhuǎn)速軌跡圖，圖（b) 是應(yīng)用傳統(tǒng)策略的模擬器轉(zhuǎn)矩軌跡圖。
[0034] 圖4為對(duì)比例II的軌跡圖，圖（a)是應(yīng)用本發(fā)明后的模擬器轉(zhuǎn)速軌跡圖，圖（b)是應(yīng) 用本發(fā)明后的模擬器轉(zhuǎn)矩軌跡圖。
[0035] 圖5為實(shí)施例1II中應(yīng)用本發(fā)明所選取的濾波器參數(shù)后在湍流風(fēng)速下WTS穩(wěn)定運(yùn)行 示意圖。
【具體實(shí)施方式】
[0036] 本發(fā)明公開(kāi)了一種基于風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量補(bǔ)償策略的濾波參數(shù)選取方法，包 含以下步驟：
[0037] 步驟1、確定風(fēng)力機(jī)模擬器系統(tǒng)的基本參數(shù)，包括通信時(shí)滯步長(zhǎng)ko,轉(zhuǎn)動(dòng)慣量補(bǔ)償 倍數(shù)Jt/Js，其中Jt為風(fēng)力機(jī)轉(zhuǎn)動(dòng)慣量，J s為風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量;時(shí)滯步長(zhǎng)ko的確定公式 為：
[0038] ko=[i/T]
[0039] 式中，τ是通信時(shí)滯時(shí)長(zhǎng)，T是系統(tǒng)采樣時(shí)長(zhǎng)。
[0040] 步驟2、確定一階數(shù)字濾波器的濾波參數(shù)α;所述一階數(shù)字濾波器的濾波參數(shù)α的初 值為0.5。
[0041] 步驟3、根據(jù)步驟1和步驟2確定的參數(shù)，判斷此時(shí)模擬器系統(tǒng)是否滿(mǎn)足穩(wěn)定性條 件;如果滿(mǎn)足則執(zhí)行步驟5，否則執(zhí)行步驟4;
[0042] 判斷模擬器系統(tǒng)是否滿(mǎn)足穩(wěn)定性條件具體包括以下步驟：
[0043] 步驟3-1、構(gòu)建方程組
[0044]
[0045] 式中，r是極徑，Θ是極角；
[0046] 步驟3-2、對(duì)方程組進(jìn)行初始化，具體是:將極角Θ設(shè)置為〇;判斷此時(shí)時(shí)滯步長(zhǎng)ko是 否為奇數(shù)，若為奇數(shù)，則執(zhí)行步驟3-3;否則執(zhí)行步驟3-4;
[0047]步驟3-3、判斷Jt/Js〈2/(l-a)是否滿(mǎn)足，若滿(mǎn)足則執(zhí)行步驟3-4,若不滿(mǎn)足則執(zhí)行 步驟3-9;
[0048] 步驟3-4、極角Θ進(jìn)行自增，具體為:θ = θ+ Δ θ，Δ Θ取值為〇. 〇1;
[0049] 步驟3-5、根據(jù)下述公式計(jì)算極徑r，
[0050]
[00511步驟3-6、對(duì)r值進(jìn)行校驗(yàn)，校驗(yàn)公式為：
[0052]
[0053] 式中ε是校驗(yàn)精度，取值為IX ΚΓ3;若校驗(yàn)通過(guò)，則執(zhí)行下一步，若不通過(guò)返回執(zhí)行 步驟3-4;
[0054]步驟3-7、判斷極徑r是否小于1，若是則執(zhí)行步驟3-8,否則直接執(zhí)行步驟3-9;
[0055]步驟3-8、判斷極角Θ是否大于231，若大于231，則風(fēng)力機(jī)模擬器系統(tǒng)的穩(wěn)定性條件成 立，否則返回執(zhí)行步驟3-4;
[0056] 步驟3-9、風(fēng)力機(jī)模擬器系統(tǒng)的穩(wěn)定性條件不成立。
[0057] 步驟4、對(duì)一階數(shù)字濾波器的濾波參數(shù)α進(jìn)行自增，具體為:α = α+ Δ α，Δ α取〇 . 〇1， 之后返回步驟3;
[0058] 步驟5、濾波參數(shù)選取結(jié)束，轉(zhuǎn)動(dòng)慣量補(bǔ)償有效，風(fēng)力機(jī)模擬器系統(tǒng)能夠穩(wěn)定運(yùn)行。
[0059] 下面結(jié)合實(shí)施例對(duì)本發(fā)明作進(jìn)一步詳細(xì)的描述：
[0060] WTS實(shí)驗(yàn)平臺(tái)建立在實(shí)驗(yàn)室內(nèi)，其主要的構(gòu)成部分如圖2有：1)一臺(tái)三相異步電動(dòng) 機(jī)與一臺(tái)永磁同步發(fā)電機(jī)構(gòu)建的對(duì)拖系統(tǒng)，2)1024脈沖/轉(zhuǎn)的旋轉(zhuǎn)編碼器用于轉(zhuǎn)速測(cè)量，3) 基于VAC0N變頻器的電機(jī)驅(qū)動(dòng)系統(tǒng)，4)基于Beckhoff PLC的實(shí)時(shí)數(shù)字控制系統(tǒng)(real-time digital control system,RTDCS),一階低通數(shù)字濾波器處于該系統(tǒng)中，5)Profibus_DP現(xiàn) 場(chǎng)總線的通訊網(wǎng)絡(luò)。
[0061] RTDCS根據(jù)包括模擬的湍流風(fēng)速，氣動(dòng)轉(zhuǎn)矩和補(bǔ)償轉(zhuǎn)矩計(jì)算電動(dòng)機(jī)的轉(zhuǎn)矩參考值， 一階低通濾波器在補(bǔ)償轉(zhuǎn)矩的計(jì)算回路中。然后將濾波后的轉(zhuǎn)矩參考值通過(guò)Prof ibus-DP 總線傳輸?shù)诫妱?dòng)機(jī)驅(qū)動(dòng)系統(tǒng)。通過(guò)工業(yè)驅(qū)動(dòng)技術(shù)，電機(jī)的運(yùn)行轉(zhuǎn)矩得以精確控制去跟隨變 化的轉(zhuǎn)矩參考值。
[0062]風(fēng)力機(jī)模擬器平臺(tái)參數(shù)如表1所示。電磁轉(zhuǎn)矩計(jì)算公式為：
[0063] Tg = k〇pt · ω g2
[0064] 這是風(fēng)力發(fā)電機(jī)一種廣泛應(yīng)用的最大功率點(diǎn)跟蹤（maximum power point tracking，MPPT)控制方法，被稱(chēng)為最優(yōu)轉(zhuǎn)矩控制，其中k〇Pt是最優(yōu)轉(zhuǎn)矩增益（（Nm/(rad/s )2)。
[0065] 表1 WTS實(shí)驗(yàn)平臺(tái)參數(shù)
[0066]
[0067]對(duì)比例I:首先，令WTS的轉(zhuǎn)動(dòng)慣量補(bǔ)償倍數(shù)為Jt/Js = 3，選取一階濾波器參數(shù)α = 0.5。此時(shí)風(fēng)力機(jī)模擬器在階躍風(fēng)速下運(yùn)行的轉(zhuǎn)速軌跡分別如圖3(a)所示，因?yàn)橥ㄐ艜r(shí)滯， 一階濾波器存在失效現(xiàn)象(參數(shù)選取不當(dāng)），使得風(fēng)力機(jī)模擬器出現(xiàn)補(bǔ)償轉(zhuǎn)矩振蕩的失穩(wěn)現(xiàn) 象，如圖3(b)所示。
[0068]對(duì)比例II:通信時(shí)滯的存在導(dǎo)致一階濾波器失效并引發(fā)模擬器系統(tǒng)失穩(wěn)，故需要 對(duì)其重新進(jìn)行參數(shù)選取，保證風(fēng)力機(jī)模擬器能夠安全穩(wěn)定運(yùn)行。
[0069]確定通信時(shí)滯τ約為60ms，即取得時(shí)滯步長(zhǎng)ko = 6,轉(zhuǎn)動(dòng)慣量補(bǔ)償倍數(shù)Jt/Js = 3, 一 階數(shù)字濾波器的濾波參數(shù)選取過(guò)程如下，并可參照?qǐng)D1的具體流程圖：
[0070] 1)選擇一階數(shù)字濾波器的濾波參數(shù)α初始值0.5;
[0071] 2)針對(duì)給定的參數(shù)ko，Jt/Js&及濾波器參數(shù)α，判斷此時(shí)風(fēng)力機(jī)模擬器系統(tǒng)的是否 穩(wěn)定。具體判斷過(guò)程如下，并可參照?qǐng)D1的判斷流程圖。
[0072] 2-1)初始化極角θ = 〇,判斷此時(shí)時(shí)滯步長(zhǎng)ko是否為奇數(shù)。若為奇數(shù)，則直接執(zhí)行步 驟2-2;若ko為偶數(shù)，則判斷Jt/Js〈2/(l-a)是否滿(mǎn)足。若滿(mǎn)足則執(zhí)行步驟2-2,否則直接進(jìn)入 步驟2-6;
[0073] 2_2)θ自增：Θ = Θ+Δ θ，Δ Θ可取0.01;
[0074] 2-3)根據(jù)式(ii)計(jì)算此時(shí)Θ對(duì)應(yīng)的極徑r，對(duì)極徑值進(jìn)行校驗(yàn)，校驗(yàn)公式為：
[0075]
[0076] 式中ε是校驗(yàn)精度，可取IX 10-3。若校驗(yàn)通過(guò)，則進(jìn)入下一步，若不通過(guò)返回執(zhí)行步 驟 2-2;
[0077] 2-4)判斷極徑r是否小于1，若是則執(zhí)行步驟2-5,否則直接執(zhí)行步驟2-6;
[0078] 2-5)判斷Θ是否大于231，若大于231，則風(fēng)力機(jī)模擬器系統(tǒng)的穩(wěn)定性條件成立，否則 返回執(zhí)行步驟2-2;
[0079] 2-6)風(fēng)力機(jī)模擬器系統(tǒng)的穩(wěn)定性條件不成立；
[0080] 3)若步驟2)中穩(wěn)定性條件不成立，一階數(shù)字濾波器的濾波參數(shù)α自增，α = α+Δα， Δα 取 0.01;
[0081] 4)若步驟2)中穩(wěn)定性條件成立，則循環(huán)結(jié)束。經(jīng)上述循環(huán)步驟后取得最佳的濾波 器系數(shù)α = 0.87;
[0082] 5)應(yīng)用有上述步驟獲得的一階濾波器系數(shù)后，風(fēng)力機(jī)模擬器在階躍風(fēng)速下轉(zhuǎn)速運(yùn) 行軌跡如圖4(a)所示。從圖中可見(jiàn)，風(fēng)力機(jī)模擬器能夠穩(wěn)定運(yùn)行，補(bǔ)償轉(zhuǎn)矩收斂穩(wěn)定，如圖4 (b)所示，不存在原先的失穩(wěn)振蕩現(xiàn)象。這說(shuō)明了本發(fā)明所提出的濾波器參數(shù)選取能夠解決 因一階濾波器參數(shù)選取不當(dāng)而造成的WTS失穩(wěn)問(wèn)題。
[0083]實(shí)施例1II:最后，對(duì)湍流風(fēng)速輸入下驗(yàn)證WTS在不同轉(zhuǎn)動(dòng)慣量補(bǔ)償倍數(shù)下的穩(wěn)定 性。實(shí)驗(yàn)結(jié)果如圖5所示，利用本發(fā)明中設(shè)計(jì)的一階濾波器參數(shù)選取規(guī)則，WTS能夠穩(wěn)定且準(zhǔn) 確的模擬實(shí)際風(fēng)力機(jī)的機(jī)械動(dòng)態(tài)。
[0084]由上可知，本發(fā)明的方法簡(jiǎn)單易行且能夠保證風(fēng)力機(jī)模擬器能夠穩(wěn)定運(yùn)行，使WTS 能模擬更大轉(zhuǎn)動(dòng)慣量的風(fēng)力機(jī)模型，并且能更加有效和方便地復(fù)現(xiàn)實(shí)際風(fēng)力機(jī)的慢動(dòng)態(tài)。
【主權(quán)項(xiàng)】
1. 一種基于風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量補(bǔ)償策略的濾波參數(shù)選取方法，其特征在于，包含 以下步驟： 步驟1、確定風(fēng)力機(jī)模擬器系統(tǒng)的基本參數(shù)，包括通信時(shí)滯步長(zhǎng)ko,轉(zhuǎn)動(dòng)慣量補(bǔ)償倍數(shù) Jt/Js，其中Jt為風(fēng)力機(jī)轉(zhuǎn)動(dòng)慣量，Js為風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量； 步驟2、確定一階數(shù)字濾波器的濾波參數(shù)α; 步驟3、根據(jù)步驟1和步驟2確定的參數(shù)，判斷此時(shí)模擬器系統(tǒng)是否滿(mǎn)足穩(wěn)定性條件；如 果滿(mǎn)足則執(zhí)行步驟5，否則執(zhí)行步驟4; 步驟4、對(duì)一階數(shù)字濾波器的濾波參數(shù)α進(jìn)行自增，具體為:α = α+Δ α，Δ α取〇.〇1，之后 返回步驟3; 步驟5、濾波參數(shù)選取結(jié)束，轉(zhuǎn)動(dòng)慣量補(bǔ)償有效，風(fēng)力機(jī)模擬器系統(tǒng)能夠穩(wěn)定運(yùn)行。2. 根據(jù)權(quán)利要求1所述的基于風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量補(bǔ)償策略的濾波參數(shù)選取方法， 其特征在于，步驟1中時(shí)滯步長(zhǎng)ko的確定公式為： ko= [τ/Τ] 式中，τ是通信時(shí)滯時(shí)長(zhǎng)，Τ是系統(tǒng)采樣時(shí)長(zhǎng)。3. 根據(jù)權(quán)利要求1所述的基于風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量補(bǔ)償策略的濾波參數(shù)選取方法， 其特征在于，步驟2中一階數(shù)字濾波器的濾波參數(shù)α的初值為0.5。4. 根據(jù)權(quán)利要求1所述的基于風(fēng)力機(jī)模擬器轉(zhuǎn)動(dòng)慣量補(bǔ)償策略的濾波參數(shù)選取方法， 其特征在于，步驟3中判斷模擬器系統(tǒng)是否滿(mǎn)足穩(wěn)定性條件具體包括以下步驟： 步驟3-1、構(gòu)建方程組式中，r是極徑，Θ是極角； 步驟3-2、對(duì)方程組進(jìn)行初始化，具體是:將極角Θ設(shè)置為〇;判斷此時(shí)時(shí)滯步長(zhǎng)ko是否為 奇數(shù)，若為奇數(shù)，則執(zhí)行步驟3-3;否則執(zhí)行步驟3-4; 步驟3-3、判斷Jt/Js〈2/(l-a)是否滿(mǎn)足，若滿(mǎn)足則執(zhí)行步驟3-4,若不滿(mǎn)足則執(zhí)行步驟 3-9； 步驟3-4、極角Θ進(jìn)行自增，具體為：Θ = Θ+Δ θ，Δ Θ取值為〇.〇1; 步驟3-5、根據(jù)下述公式計(jì)算極徑r，步驟3-6、對(duì)r值進(jìn)行校驗(yàn)，校驗(yàn)公式為：式中ε是校驗(yàn)精度，取值為IX 10-3;若校驗(yàn)通過(guò)，則執(zhí)行下一步，若不通過(guò)返回執(zhí)行步驟 3-4; 步驟3-7、判斷極徑r是否小于1，若是則執(zhí)行步驟3-8,否則直接執(zhí)行步驟3-9; 步驟3-8、判斷極角Θ是否大于231，若大于231，則風(fēng)力機(jī)模擬器系統(tǒng)的穩(wěn)定性條件成立， 否則返回執(zhí)行步驟3-4; 步驟3-9、風(fēng)力機(jī)模擬器系統(tǒng)的穩(wěn)定性條件不成立。
【文檔編號(hào)】F03D17/00GK105952587SQ201610272758
【公開(kāi)日】2016年9月21日
【申請(qǐng)日】2016年4月28日
【發(fā)明人】殷明慧, 徐洋超, 李偉杰, 何立君, 李群, 劉建坤, 周前, 陳兵, 汪成根, 卜京, 謝云云, 姚娟, 鄒云, 陳哲, 張寧宇, 衛(wèi)鵬
【申請(qǐng)人】南京理工大學(xué), 國(guó)網(wǎng)江蘇省電力公司電力科學(xué)研究院