輸電線路風(fēng)偏瞬態(tài)響應(yīng)計(jì)算方法
【專利摘要】本發(fā)明提出了一種考慮氣動(dòng)力系數(shù)隨風(fēng)攻角變化的輸電線路風(fēng)偏瞬態(tài)響應(yīng)計(jì)算方法；通過風(fēng)洞試驗(yàn)測(cè)量裸導(dǎo)線及覆冰導(dǎo)線氣動(dòng)力系數(shù)；考慮輸電線路各點(diǎn)的空間相關(guān)性以及同一點(diǎn)縱向和豎向風(fēng)速的相關(guān)性，利用諧波疊加法模擬三維脈動(dòng)風(fēng)速；建立連續(xù)多跨多分裂輸電線路有限元模型，利用每一時(shí)刻的合風(fēng)速及隨風(fēng)攻角變化的氣動(dòng)力系數(shù)精確計(jì)算風(fēng)荷載，對(duì)導(dǎo)線風(fēng)偏進(jìn)行時(shí)程分析。
【專利說明】
輸電線路風(fēng)偏瞬態(tài)響應(yīng)計(jì)算方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明屬于輸電線路技術(shù)領(lǐng)域，具體涉及一種輸電線路風(fēng)偏瞬態(tài)響應(yīng)計(jì)算方法。
【背景技術(shù)】
[0002] 風(fēng)荷載是輸電線路的主要荷載之一，由于導(dǎo)線屬于柔性結(jié)構(gòu)，在強(qiáng)風(fēng)作用下產(chǎn)生 較大變形，當(dāng)導(dǎo)線與桿塔或不同相導(dǎo)線之間的空氣間隙不滿足絕緣條件時(shí)便會(huì)被擊穿，發(fā) 生風(fēng)偏閃絡(luò)事故，嚴(yán)重影響線路正常運(yùn)行，造成較大的經(jīng)濟(jì)損失。為有效預(yù)防風(fēng)偏閃絡(luò)事故 的發(fā)生，國內(nèi)外相關(guān)領(lǐng)域的學(xué)者對(duì)導(dǎo)線風(fēng)偏進(jìn)行了大量研究，通常建立連續(xù)多跨有限元模 型，利用計(jì)算機(jī)模擬脈動(dòng)風(fēng)速，與平均風(fēng)速疊加后得到總風(fēng)速，根據(jù)準(zhǔn)定常假設(shè)轉(zhuǎn)換為風(fēng)荷 載并施加在模型上，利用有限單元法求解風(fēng)偏瞬態(tài)響應(yīng)。
[0003] 對(duì)于風(fēng)速，以往研究只考慮縱向（順風(fēng)向）風(fēng)速，不考慮豎向脈動(dòng)風(fēng)速的影響。風(fēng)流 經(jīng)導(dǎo)線分別產(chǎn)生與來流方向相同的阻力、與來流方向垂直的升力和與導(dǎo)線截面垂直的扭 矩，它們對(duì)應(yīng)的無量綱系數(shù)分別稱為阻力系數(shù)、升力系數(shù)和扭轉(zhuǎn)系數(shù)。以往研究和規(guī)范在計(jì) 算導(dǎo)線風(fēng)荷載時(shí)只考慮阻力，且阻力系數(shù)取固定值，我國《110kV~750kV架空輸電線路設(shè)計(jì) 規(guī)范》(GB 50545-2010)和《架空輸電線路桿塔結(jié)構(gòu)設(shè)計(jì)技術(shù)規(guī)定》(DL/T 5154-2012)均規(guī) 定，線徑小于17mm時(shí)阻力系數(shù)取1.2,大于或等于17mm時(shí)取1.1;國外規(guī)范阻力系數(shù)一般取 1. 〇或根據(jù)可靠的實(shí)測(cè)數(shù)據(jù)和風(fēng)洞試驗(yàn)結(jié)果取值;對(duì)于分裂導(dǎo)線，在計(jì)算風(fēng)荷載時(shí)，各子導(dǎo) 線均按單導(dǎo)線計(jì)算，不考慮子導(dǎo)線間相互遮擋對(duì)阻力系數(shù)的影響。以往研究表明，導(dǎo)線氣動(dòng) 力系數(shù)隨風(fēng)攻角不斷變化，尤其是對(duì)于覆冰導(dǎo)線，這種變化更加明顯；由于上游子導(dǎo)線對(duì)下 游子導(dǎo)線的尾流干擾效應(yīng)，多分裂各子導(dǎo)線氣動(dòng)力系數(shù)隨風(fēng)攻角變化規(guī)律存在一定差異。

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

[0004] 鑒于以往風(fēng)偏計(jì)算中存在的問題，本發(fā)明提出了一種更加精確的輸電線路風(fēng)偏瞬 態(tài)響應(yīng)計(jì)算方法，計(jì)算結(jié)果更加符合真實(shí)情況，可對(duì)容易發(fā)生較大風(fēng)偏的位置進(jìn)行預(yù)測(cè)，指 導(dǎo)輸電線路防風(fēng)偏設(shè)計(jì)，從而降低風(fēng)偏閃絡(luò)事故發(fā)生的概率。
[0005] 為了解決上述技術(shù)問題，本發(fā)明采用如下技術(shù)方案:輸電線路風(fēng)偏瞬態(tài)響應(yīng)計(jì)算 方法，所述方法包括以下步驟：
[0006] 步驟1:制作精確模擬導(dǎo)線表面粗糙度的裸導(dǎo)線和覆冰導(dǎo)線模型，通過高頻天平測(cè) 力風(fēng)洞試驗(yàn)測(cè)定不同風(fēng)攻角下裸導(dǎo)線和覆冰導(dǎo)線各子導(dǎo)線的氣動(dòng)力系數(shù)；
[0007] 步驟2:根據(jù)某真實(shí)輸電線路的動(dòng)力學(xué)特性和結(jié)構(gòu)形式，選擇合適的單元類型，建 立包含各子導(dǎo)線、絕緣子串、間隔棒、金具在內(nèi)的精細(xì)化有限元模型；
[0008] 步驟3:確定地貌類型，根據(jù)導(dǎo)線所在位置處的平均風(fēng)速以及湍流特性，考慮輸電 線路各點(diǎn)的空間相關(guān)性以及同一點(diǎn)縱向和豎向風(fēng)速的相關(guān)性，利用諧波疊加法模擬三維脈 動(dòng)風(fēng)速；
[0009] 步驟4:根據(jù)某時(shí)刻的合風(fēng)速以及該時(shí)刻風(fēng)攻角對(duì)應(yīng)的氣動(dòng)力系數(shù)，精確計(jì)算各子 導(dǎo)線的阻力和升力，并施加在有限元模型上；
[0010] 步驟5:利用有限單元法求解輸電線路各點(diǎn)風(fēng)偏瞬態(tài)響應(yīng)。
[0011] 所述步驟1在邊界層風(fēng)洞中完成，通過風(fēng)洞試驗(yàn)測(cè)定分裂導(dǎo)線各子導(dǎo)線阻力系數(shù) 〇〇」(9」)和升力系數(shù)〇^(0」），其中偽子導(dǎo)線編號(hào)，0」為子導(dǎo)線」的風(fēng)攻角。
[0012] 所述步驟2根據(jù)某真實(shí)輸電線路的導(dǎo)線、絕緣子串、間隔棒、金具的物理參數(shù)，建立 多分裂精細(xì)化有限元模型，所述物理參數(shù)包括幾何尺寸、密度、剛度。
[0013] 所述步驟4考慮豎向脈動(dòng)風(fēng)速后，合風(fēng)速方向時(shí)刻改變，合風(fēng)速與多分裂導(dǎo)線的夾 角即風(fēng)攻角時(shí)刻改變，相應(yīng)的氣動(dòng)力系數(shù)亦隨風(fēng)攻角變化，可精確計(jì)算得到導(dǎo)線風(fēng)荷載。
[0014] 所述步驟5重復(fù)步驟4中的過程，求解輸電線路模型各點(diǎn)縱向和豎向瞬態(tài)位移，進(jìn) 而計(jì)算懸垂絕緣子串的瞬態(tài)風(fēng)偏角。
[0015] 本發(fā)明具有以下有益效果：
[0016] 1.導(dǎo)線風(fēng)荷載計(jì)算準(zhǔn)確。氣動(dòng)力系數(shù)是導(dǎo)線風(fēng)荷載計(jì)算的重要參數(shù)，通過風(fēng)洞試 驗(yàn)準(zhǔn)確測(cè)定?？紤]客觀存在的豎向脈動(dòng)風(fēng)速以及隨風(fēng)攻角變化的氣動(dòng)力系數(shù)，精確計(jì)算風(fēng) 偏過程中各時(shí)刻風(fēng)荷載。
[0017] 2.導(dǎo)線風(fēng)偏響應(yīng)計(jì)算準(zhǔn)確。根據(jù)真實(shí)輸電線路物理參數(shù)建立包含各子導(dǎo)線、絕緣 子串、間隔棒、金具等在內(nèi)的精細(xì)化有限元模型，施加風(fēng)荷載后得到的風(fēng)偏響應(yīng)與真實(shí)情況 接近。
[0018] 3.具有普遍適用性。可根據(jù)不同地貌、不同輸電線路、導(dǎo)線的覆冰情況等實(shí)時(shí)差異 調(diào)整物理和數(shù)值模型的計(jì)算參數(shù)。
【附圖說明】
[0019] 圖1為本發(fā)明的步驟流程圖。
[0020] 圖2為八分裂導(dǎo)線子導(dǎo)線編號(hào)和風(fēng)攻角定義。
[0021 ]圖3為B類地貌下某1000kV八分裂輸電線路。
[0022]圖4為空間任意一點(diǎn)風(fēng)速場(chǎng)的方向定義方法。
[0023]圖5為三維脈動(dòng)風(fēng)速時(shí)程。
[0024]圖6為八分裂導(dǎo)線風(fēng)偏狀態(tài)示意圖。
[0025]圖7為懸垂絕緣子串風(fēng)偏角計(jì)算方法示意圖。
[0026]圖8為懸垂絕緣子串風(fēng)偏角時(shí)程。
【具體實(shí)施方式】
[0027]下面結(jié)合附圖對(duì)本發(fā)明作進(jìn)一步詳細(xì)說明。
[0028] 如圖1所示，輸電線路風(fēng)偏瞬態(tài)響應(yīng)計(jì)算方法：
[0029] 步驟1:通過風(fēng)洞試驗(yàn)測(cè)定分裂導(dǎo)線各子導(dǎo)線阻力系數(shù)和升力系數(shù)αχθ」）。
[0030] 其中j為子導(dǎo)線編號(hào)，0」為子導(dǎo)線j的風(fēng)攻角。實(shí)際測(cè)量時(shí)，可僅在某一根子導(dǎo)線j 下方安裝測(cè)力天平，直接測(cè)量其阻力FDj㈨)和升力凡辦丄并計(jì)算其氣動(dòng)力系數(shù)，根據(jù)對(duì)稱 性和周期性，得到其他子導(dǎo)線的氣動(dòng)力系數(shù)。以八分裂導(dǎo)線為例，子導(dǎo)線編號(hào)和風(fēng)攻角定義 如圖2所示。
[0031] 步驟2:建立連續(xù)多跨多分裂輸電線路有限元模型。
[0032] 以如圖3所示的B類地貌下某1000kV八分裂輸電線路為例，根據(jù)導(dǎo)線、絕緣子串、間 隔棒和金具的物理參數(shù)，利用有限元軟件進(jìn)行建模，導(dǎo)線初始形狀為懸鏈線，施加重力進(jìn)行 找形分析并更新各節(jié)點(diǎn)坐標(biāo)。
[0033]步驟3:通過諧波疊加法模擬三維脈動(dòng)風(fēng)速。
[0034]空間任意一點(diǎn)風(fēng)速都可以分解為相互垂直的三個(gè)分量，如圖4所示，其中U為縱向 (順風(fēng)向）風(fēng)速，W為豎向風(fēng)速，V為橫向風(fēng)速。對(duì)于平坦地貌，豎向和橫向平均風(fēng)速為0，因此：
[0036] 式中：(6/為縱向平均風(fēng)速;u、w、v分別為縱向、豎向和橫向脈動(dòng)風(fēng)速。
[0037]本發(fā)明采用諧波疊加法中基于互譜密度矩陣特征正交分解的譜表示法(P0D)來模 擬三維脈動(dòng)風(fēng)速。自譜密度函數(shù)采用ESDU推薦的Karman譜，縱向平均風(fēng)速按我國《建筑結(jié)構(gòu) 荷載規(guī)范》(GB 5009-2012)B類地貌風(fēng)剖面計(jì)算，縱向湍流度按日本《AIJ Recommendations for Loads on Buildings》Π 類地貌（對(duì)應(yīng)我國B類地貌）計(jì)算；空間相干函數(shù)采用 Davenport推薦的形式;對(duì)于空間任意一點(diǎn)，只考慮縱向和豎向脈動(dòng)風(fēng)速之間的相關(guān)性。計(jì) 算互譜密度矩陣并進(jìn)行P0D分解，利用得到的特征向量、特征值以及隨機(jī)生成的相位角可以 得到脈動(dòng)風(fēng)速的表達(dá)式。為提高計(jì)算效率，引入快速傅里葉變換。圖3中輸電線路第二跨跨 中處的三維脈動(dòng)風(fēng)速時(shí)程如圖5所示。
[0038] 步驟4:精確計(jì)算風(fēng)偏過程中各時(shí)刻輸電線路風(fēng)荷載。
[0039] 由于橫向脈動(dòng)風(fēng)速與輸電線路的走向相同，對(duì)風(fēng)偏影響較小，因此風(fēng)偏計(jì)算中可 以忽略。以八分裂導(dǎo)線為例，導(dǎo)線某斷面在風(fēng)偏過程中運(yùn)動(dòng)狀態(tài)的改變?nèi)鐖D6所示?？紤]導(dǎo) 線與來流之間的相對(duì)運(yùn)動(dòng)，t時(shí)刻子導(dǎo)線j處（圖6中以子導(dǎo)線8為例，其他子導(dǎo)線同理）的縱 向風(fēng)速山、豎向風(fēng)速^、合風(fēng)速V句按下式計(jì)算：
[0043] 式中：〇氣+和內(nèi),分別為子導(dǎo)線j縱向和豎向運(yùn)動(dòng)速度，利用NewmarkP法計(jì)算;uj和wj 分別為子導(dǎo)線j處縱向和豎向脈動(dòng)風(fēng)速;@7為子導(dǎo)線j處縱向平均風(fēng)速。
[0044] 風(fēng)攻角為合風(fēng)速與導(dǎo)線的夾角，t時(shí)刻子導(dǎo)線j處風(fēng)攻角~可表示為：
[0045] θ」=Ω-φ」(5)
[0046] 式中：Ω為導(dǎo)線斷面轉(zhuǎn)過的角度；為子導(dǎo)線j處縱向風(fēng)速與合風(fēng)速夾角。
[0047]子導(dǎo)線j處阻力FDj與合風(fēng)速方向同向，升力Fu與合風(fēng)速方向垂直，根據(jù)準(zhǔn)定常假 設(shè)可分別表示為：
[0048] FDj(0j，t)=O.5pV合j(t)2CDj(0j)Aj (6)
[0049] FLj(0j,t)=〇.5pV-j(t)2CLj(0j)Aj (7)式中：P 為空氣密度;Aj 為迎風(fēng)面積。
[0050] 將阻力和升力分別沿X軸和Z軸投影，得到縱向荷載和豎向荷載：
[0051] Fxj = FDj · cos Φ j+Fij · sin Φ j (8)
[0052] Fzj = FDj · sin Φ j-Fij · cos Φ j_Gj (9)
[0053] 式中:Gj為導(dǎo)線重力。
[0054] 為了消除在導(dǎo)線上施加突加荷載后的沖擊放大效應(yīng)，在模擬的風(fēng)速時(shí)程前增加一 段風(fēng)速由0線性增至第一個(gè)時(shí)間點(diǎn)風(fēng)速的過程。
[0055] 對(duì)于某一斷面，= 0時(shí)刻，八分裂導(dǎo)線各子導(dǎo)線初始狀態(tài)坐標(biāo)已知，t時(shí)刻子導(dǎo)線 j風(fēng)攻角~(ω=〇,對(duì)應(yīng)的阻力系數(shù)為％%(ω]，升力系數(shù)為ακθχω]，利用公式⑷~ (9)計(jì)算tl時(shí)刻荷載，施加 tl時(shí)刻荷載后，得到t2= A t時(shí)刻的各子導(dǎo)線位移，進(jìn)而得到t2時(shí) 刻各子導(dǎo)線坐標(biāo)，利用ti、t2時(shí)刻坐標(biāo)可以計(jì)算出各子導(dǎo)線兩兩之間連線轉(zhuǎn)過的角度，分別 記為 Ω I2(t2)、Ω 23(t2)、Ω 34(t2)、Ω 45(t2)、Ω 56(t2)、Ω 67(t2)、Ω 78(t2)、Ω 81(t2)，由于t2時(shí) 刻八分裂導(dǎo)線之間的相對(duì)位置與初始狀態(tài)不同，不再是正八邊形，t2時(shí)刻導(dǎo)線斷面轉(zhuǎn)過的 角度Ω(?2)可表示為：
[0056] Ω (t2) = [ Ω I2(t2) + Ω 23(t2) + Ω 34(t2) + Ω 45(t2) + Ω 56(t2) + Ω 67(t2) + Ω 78(t2) + Ω 81 (t2)]/8 (10)
[0057] 利用公式(5)計(jì)算t2時(shí)刻子導(dǎo)線j風(fēng)攻角~(t2)，對(duì)應(yīng)的阻力系數(shù)為Cw%(t 2)]，升 力系數(shù)為(U0j(t2)]，利用公式(6)~(9)計(jì)算t2時(shí)刻荷載，施加 t2時(shí)刻荷載后，得到t3 = 2A t時(shí)刻的各子導(dǎo)線位移，再利用ti、t3時(shí)刻坐標(biāo)計(jì)算Ω (t3)，對(duì)應(yīng)的氣動(dòng)力系數(shù)為(^[θ」 (t3) ]、CLJ [ Θj (t3)]，進(jìn)而得到t3時(shí)亥IJ荷載，以此類推。
[0058]步驟5:重復(fù)步驟4中的過程，求解輸電線路模型各點(diǎn)縱向和豎向瞬態(tài)位移，進(jìn)而計(jì) 算懸垂絕緣子串的瞬態(tài)風(fēng)偏角。
[0059] 利用有限元軟件求解時(shí)，考慮大變形和應(yīng)力剛化效應(yīng)，得到輸電線路模型各點(diǎn)縱 向位移Ux(t)和豎向位移Uz(t)。如圖7所示，懸垂絕緣子串長度為Lj，則風(fēng)偏角祆?)可表示 為：
[0060] (p{!) = arcs\n(Ux(t)/Lj) (11)
[0061] 圖3中輸電線路掛點(diǎn)2處懸垂絕緣子串風(fēng)偏角時(shí)程如圖8所示。
[0062] 在強(qiáng)風(fēng)作用下，輸電線路產(chǎn)生較大變形，當(dāng)導(dǎo)線與桿塔或不同相導(dǎo)線之間的空氣 間隙不滿足絕緣條件時(shí)便會(huì)被擊穿，發(fā)生風(fēng)偏閃絡(luò)事故，嚴(yán)重影響線路正常運(yùn)行，造成較大 的經(jīng)濟(jì)損失，準(zhǔn)確計(jì)算導(dǎo)線風(fēng)偏對(duì)有效預(yù)防風(fēng)偏閃絡(luò)事故的發(fā)生具有重要意義。以往風(fēng)偏 計(jì)算方法中，導(dǎo)線阻力系數(shù)取固定值，且不考慮豎向脈動(dòng)風(fēng)速和升力，但實(shí)際風(fēng)偏過程中， 豎向脈動(dòng)風(fēng)速是客觀存在的，合風(fēng)速的方向不斷變化，分裂導(dǎo)線也會(huì)發(fā)生轉(zhuǎn)動(dòng)，因此，合風(fēng) 速與導(dǎo)線夾角，即風(fēng)攻角處于動(dòng)態(tài)變化中。考慮到多分裂導(dǎo)線的阻力和升力系數(shù)會(huì)隨風(fēng)攻 角不斷變化，尤其是對(duì)于覆冰導(dǎo)線而言，該變化尤為明顯，可見，以往的風(fēng)偏計(jì)算方法不夠 準(zhǔn)確。為此，本發(fā)明提出了一種考慮氣動(dòng)力系數(shù)隨風(fēng)攻角變化的輸電線路風(fēng)偏瞬態(tài)響應(yīng)計(jì) 算方法;通過風(fēng)洞試驗(yàn)測(cè)量裸導(dǎo)線及覆冰導(dǎo)線隨風(fēng)攻角變化的氣動(dòng)力系數(shù);考慮輸電線路 各點(diǎn)的空間相關(guān)性以及同一點(diǎn)縱向和豎向風(fēng)速的相關(guān)性，利用諧波疊加法模擬三維脈動(dòng)風(fēng) 速;建立連續(xù)多跨多分裂輸電線路有限元模型，利用每一時(shí)刻的合風(fēng)速及隨風(fēng)攻角變化的 氣動(dòng)力系數(shù)精確計(jì)算風(fēng)荷載，對(duì)導(dǎo)線風(fēng)偏進(jìn)行時(shí)程分析。
【主權(quán)項(xiàng)】
1. 輸電線路風(fēng)偏瞬態(tài)響應(yīng)計(jì)算方法，其特征在于:所述方法包括以下步驟： 步驟1:制作精確模擬導(dǎo)線表面粗糙度的裸導(dǎo)線和覆冰導(dǎo)線模型，通過高頻天平測(cè)力風(fēng) 洞試驗(yàn)測(cè)定不同風(fēng)攻角下裸導(dǎo)線和覆冰導(dǎo)線各子導(dǎo)線的氣動(dòng)力系數(shù)； 步驟2:根據(jù)某真實(shí)輸電線路的動(dòng)力學(xué)特性和結(jié)構(gòu)形式，選擇合適的單元類型，建立包 含各子導(dǎo)線、絕緣子串、間隔棒、金具在內(nèi)的精細(xì)化有限元模型； 步驟3:確定地貌類型，根據(jù)導(dǎo)線所在位置處的平均風(fēng)速以及湍流特性，考慮輸電線路 各點(diǎn)的空間相關(guān)性以及同一點(diǎn)縱向和豎向風(fēng)速的相關(guān)性，利用諧波疊加法模擬三維脈動(dòng)風(fēng) 速； 步驟4:根據(jù)某時(shí)刻的合風(fēng)速以及該時(shí)刻風(fēng)攻角對(duì)應(yīng)的氣動(dòng)力系數(shù)，精確計(jì)算各子導(dǎo)線 的阻力和升力，并施加在有限元模型上； 步驟5:利用有限單元法求解輸電線路各點(diǎn)風(fēng)偏瞬態(tài)響應(yīng)。2. 根據(jù)權(quán)利要求1所述的輸電線路風(fēng)偏瞬態(tài)響應(yīng)計(jì)算方法，其特征在于:所述步驟1在 邊界層風(fēng)洞中完成，通過風(fēng)洞試驗(yàn)測(cè)定分裂導(dǎo)線各子導(dǎo)線阻力系數(shù)(^(θ〇和升力系數(shù) Cu (9j)，其中j為子導(dǎo)線編號(hào)，為子導(dǎo)線j的風(fēng)攻角。3. 根據(jù)權(quán)利要求1所述的輸電線路風(fēng)偏瞬態(tài)響應(yīng)計(jì)算方法，其特征在于:所述步驟2根 據(jù)某真實(shí)輸電線路的導(dǎo)線、絕緣子串、間隔棒、金具的物理參數(shù)，建立多分裂精細(xì)化有限元 模型，所述物理參數(shù)包括幾何尺寸、密度、剛度。4. 根據(jù)權(quán)利要求1所述的輸電線路風(fēng)偏瞬態(tài)響應(yīng)計(jì)算方法，其特征在于:所述步驟4考 慮豎向脈動(dòng)風(fēng)速后，合風(fēng)速方向時(shí)刻改變，合風(fēng)速與多分裂導(dǎo)線的夾角即風(fēng)攻角時(shí)刻改變， 相應(yīng)的氣動(dòng)力系數(shù)隨風(fēng)攻角變化，可精確計(jì)算得到導(dǎo)線風(fēng)荷載。5. 根據(jù)權(quán)利要求1所述的輸電線路風(fēng)偏瞬態(tài)響應(yīng)計(jì)算方法，其特征在于:所述步驟5重 復(fù)步驟4中的過程，求解輸電線路模型各點(diǎn)縱向和豎向瞬態(tài)位移，進(jìn)而計(jì)算懸垂絕緣子串的 瞬態(tài)風(fēng)偏角。
【文檔編號(hào)】G06F17/50GK106096105SQ201610389916
【公開日】2016年11月9日
【申請(qǐng)日】2016年6月2日 公開號(hào)201610389916.X, CN 106096105 A, CN 106096105A, CN 201610389916, CN-A-106096105, CN106096105 A, CN106096105A, CN201610389916, CN201610389916.X
【發(fā)明人】樓文娟, 徐海巍, 李天昊, 梁洪超
【申請(qǐng)人】浙江大學(xué)