董依帆， 柯世堂， 楊 青

(南京航空航天大學(xué) 土木工程系，南京 210016)

近年來隨著國(guó)內(nèi)能源結(jié)構(gòu)轉(zhuǎn)型步伐的加快，風(fēng)力發(fā)電設(shè)施陸續(xù)興建。其中我國(guó)西部風(fēng)能資源充沛[1]，為更好的捕捉風(fēng)能，風(fēng)力機(jī)必然朝著高大化方向發(fā)展。然而此地區(qū)戈壁沙漠面積廣大，沙粒受到近地邊界層風(fēng)場(chǎng)的驅(qū)動(dòng)，在一定高度范圍內(nèi)隨機(jī)運(yùn)動(dòng)，兩者之間相互耦合并發(fā)生能量遷移，使得沙粒以較大速度沖擊至風(fēng)力機(jī)表面，形成附加荷載及風(fēng)蝕效應(yīng)，對(duì)風(fēng)力機(jī)的安全性能和使用壽命產(chǎn)生顯著影響[2-3]。因此，開展氣-固兩相流理論研究，探討不同風(fēng)速和沙粒粒徑組合對(duì)風(fēng)力機(jī)氣動(dòng)載荷的作用機(jī)理，具有重要的工程意義和理論價(jià)值。

當(dāng)前國(guó)內(nèi)外對(duì)風(fēng)-沙共同作用的研究多集中在風(fēng)沙現(xiàn)象和風(fēng)沙運(yùn)動(dòng)本身，包括沙粒運(yùn)動(dòng)力學(xué)過程[4]、輸沙量的垂直分布[5]、輸沙率模型[6]、臨界起沙風(fēng)速[7]和風(fēng)速廓線規(guī)律[8]。相較于風(fēng)荷載作用，風(fēng)-沙共同作用下建筑物或構(gòu)筑物荷載效應(yīng)更加顯著[9]。文獻(xiàn)[10]通過風(fēng)洞測(cè)力試驗(yàn)，研究了不同沙濃度和風(fēng)速條件下風(fēng)沙對(duì)低矮建筑整體受力影響，主要考察風(fēng)沙環(huán)境中建筑物表面的風(fēng)壓特性；文獻(xiàn)[11]模擬了大氣邊界層中單體及兩相鄰對(duì)稱高層建筑的風(fēng)沙流場(chǎng)，并指出隨著沙粒相體積分?jǐn)?shù)的增大，風(fēng)壓系數(shù)不斷增大；文獻(xiàn)[12]對(duì)風(fēng)沙環(huán)境下水平軸風(fēng)力機(jī)進(jìn)行數(shù)值模擬，研究了在某一來流風(fēng)速和不同沙粒直徑組合工況下的風(fēng)力機(jī)風(fēng)輪轉(zhuǎn)矩。已有研究均未基于風(fēng)力機(jī)整機(jī)和沙粒真實(shí)環(huán)境分布特性展開多參數(shù)條件下大型風(fēng)力機(jī)風(fēng)沙特性的系統(tǒng)研究。

鑒此，本文以南京航空航天大學(xué)自主研發(fā)的5 MW水平軸三葉片風(fēng)力機(jī)[13]為研究對(duì)象，采用改進(jìn)的k-ε湍流模型和粉塵釋放模型(Dust Production Model，DPM)，設(shè)置不同風(fēng)速和不同沙粒粒徑形成9種組合計(jì)算工況，對(duì)處于典型停機(jī)狀態(tài)的塔架-葉片耦合體系展開數(shù)值模擬，以獲取沙粒在風(fēng)力機(jī)不同部位的分布密度及撞擊速度，并對(duì)比分析風(fēng)沙共同作用下風(fēng)力機(jī)體系繞流特性、表面壓力分布隨風(fēng)沙流速度和沙粒濃度的變化規(guī)律。

1 風(fēng)-沙雙相耦合算法

1.1 沙塵暴等級(jí)劃分

我國(guó)西北地區(qū)沙塵暴天氣頻發(fā)，主要集中在春季4～5月，其發(fā)生時(shí)間集中、頻率高、強(qiáng)度大。根據(jù)探空資料，沙塵暴發(fā)生時(shí)形成的沙塵壁高達(dá)300 m，沙塵暴影響高度在2 100 m以內(nèi)[14-15]。沙塵暴天氣根據(jù)最大風(fēng)速與最小能見度劃分為4個(gè)等級(jí)，如表1所示。本文選取中、強(qiáng)、特強(qiáng)沙塵暴天氣典型風(fēng)速18 m/s，22 m/s，25 m/s作為風(fēng)場(chǎng)計(jì)算風(fēng)速。

表1 沙塵暴天氣強(qiáng)度劃分標(biāo)準(zhǔn)

1.2 風(fēng)沙流密度

風(fēng)沙流密度是指空氣運(yùn)動(dòng)速度在大于起沙風(fēng)速情況下，單位體積空氣氣流中所含沙的質(zhì)量，單位為kg/m3。在近地層中，風(fēng)受地面摩擦阻力的影響而速度降低，一般而言，摩擦力隨著高度的增加而減小，故風(fēng)速隨著高度的增加而增大，風(fēng)沙流密度與高度、風(fēng)力等均有密切關(guān)系。

據(jù)文獻(xiàn)[16]統(tǒng)計(jì)，當(dāng)風(fēng)速在40 m/s以下時(shí)，風(fēng)沙流密度的數(shù)量級(jí)一般約為10-5，且隨高度的增加而下降，基本以3 m高度作為分界點(diǎn)。在3 m高度以下，風(fēng)沙流密度隨高度的增加而大幅下降，在3 m高度以上風(fēng)沙流密度隨高度增加逐步下降，且沙粒運(yùn)動(dòng)以飛揚(yáng)為主?；诠こ贪踩钥紤]，采取極端工況進(jìn)行數(shù)值模擬，故風(fēng)場(chǎng)入口處3 m以下高度風(fēng)沙流密度取為10-4kg/m3，速度與風(fēng)速相同；3 m以上高度風(fēng)沙流密度取為10-5kg/m3，速度為0。

1.3 計(jì)算模型與顆粒運(yùn)動(dòng)方程

標(biāo)準(zhǔn)k-ε模型在科學(xué)研究及工程實(shí)際中得到了廣泛的應(yīng)用，但用于強(qiáng)旋流、彎曲壁面流動(dòng)或彎曲流線流動(dòng)時(shí)，會(huì)產(chǎn)生一定的失真。為彌補(bǔ)其缺陷，本文采用Realizablek-ε模型，Realizablek-ε模型中k與ε的模數(shù)化輸運(yùn)方程組為

(1)

(2)

與標(biāo)準(zhǔn)k-ε模型比較發(fā)現(xiàn)，Realizablek-ε模型主要變化是：湍流黏度計(jì)算公式發(fā)生了變化，引入了與旋轉(zhuǎn)和曲率有關(guān)的內(nèi)容；ε方程發(fā)生了很大變化，方程中的產(chǎn)生項(xiàng)不再包含有k方程中的產(chǎn)生項(xiàng)Gk；ε的倒數(shù)第二項(xiàng)不具有任何奇異性，即使值很小或?yàn)榱悖膊粫?huì)為零。

相比于視為連續(xù)相的空氣來流，沙粒在空氣中所占比例較小，即使在特強(qiáng)沙塵暴(風(fēng)速≥25 m/s)時(shí)，沙粒占空氣體積分?jǐn)?shù)仍遠(yuǎn)小于10%?；诖祟悓?shí)際情況，DPM模型被廣泛應(yīng)用于多相流聯(lián)合研究。本文在風(fēng)場(chǎng)計(jì)算穩(wěn)定后，將離散相模型作為第二相，介入至連續(xù)相中展開風(fēng)-沙雙向耦合運(yùn)算。沙粒在風(fēng)場(chǎng)中運(yùn)動(dòng)平衡方程為

(3)

(4)

式中：μ為流體黏性系數(shù)；dp為顆粒直徑；Re為相對(duì)雷諾數(shù)； 可表示為

(5)

考慮沙粒離散相影響后，風(fēng)連續(xù)相基本控制方程可表示為

(6)

(7)

(8)

式中：I為單位張量； 等式右邊第二項(xiàng)為體積膨脹作用。

1.4 顆粒碰撞模型

沙粒沖擊到風(fēng)力機(jī)表面過程服從動(dòng)量守恒定律，求解沖擊力的關(guān)鍵在于碰撞時(shí)間。計(jì)算中忽略沙粒在沖擊過程中可能發(fā)生的破裂現(xiàn)象，認(rèn)為沙粒與結(jié)構(gòu)間相互作用遵循牛頓第二定律并假設(shè)反彈后的速度與撞擊前的速度一致。由動(dòng)量定理

(9)

式中：f(t)為單個(gè)沙粒沖擊力矢量， N；vs為沙粒速度矢量。

沙粒在單位時(shí)間內(nèi)對(duì)結(jié)構(gòu)的沖擊力F(τ)為

(10)

將撞擊風(fēng)力機(jī)的沙粒近似看作球體，則

(11)

由于沙粒直徑在1 mm及以下，撞擊前水平末速度相對(duì)較大，可將風(fēng)力機(jī)表面簡(jiǎn)化為無限大平面[17]，將碰撞時(shí)間τ取為

(12)

(13)

式中：m為沙粒的質(zhì)量；φ為速度恢復(fù)系數(shù)取1；k1，k2取值與沙粒、風(fēng)力機(jī)材料楊氏模量、泊松比有關(guān)。

則沙粒對(duì)結(jié)構(gòu)的沖擊力可簡(jiǎn)化為

(14)

2 數(shù)值模擬

2.1 工程簡(jiǎn)介

表2給出了該5 MW風(fēng)力機(jī)的主要結(jié)構(gòu)尺寸及示意圖，葉片長(zhǎng)度為60 m，機(jī)艙尺寸18 m×6 m×6 m，各葉片沿周向成120°夾角分布，塔架通長(zhǎng)變厚，高為124 m，塔底半徑3.5 m，塔頂半徑3.0 m，塔頂厚度90 mm，塔底厚度200 mm。

表2 5 MW風(fēng)力發(fā)電結(jié)構(gòu)體系主要參數(shù)列表

2.2 工況設(shè)置

該風(fēng)力機(jī)位于B類地貌，綜合考慮我國(guó)沙塵暴天氣等級(jí)劃分，對(duì)比研究3種風(fēng)速和3種沙粒粒徑組合下極端工況對(duì)風(fēng)力機(jī)整機(jī)的作用規(guī)律。其中，小風(fēng)、中風(fēng)和大風(fēng)分別以沙塵暴等級(jí)中、強(qiáng)和特強(qiáng)的基本風(fēng)速進(jìn)行劃分；沙粒粒徑分別取為最大統(tǒng)計(jì)概率區(qū)間內(nèi)的0.10 mm，0.25 mm以及極端工況下的1 mm，3種風(fēng)速和沙粒粒徑組合共計(jì)9種工況，如表3所示。

表3 風(fēng)力機(jī)風(fēng)-沙耦合計(jì)算工況劃分表

2.3 計(jì)算域設(shè)置

如圖1所示，本文研究的風(fēng)力機(jī)均處于0°來流風(fēng)向角。為保證風(fēng)力機(jī)尾流能夠充分發(fā)展，計(jì)算域尺寸設(shè)置為12D×5D×5D(流向X×展向Y×豎向Z，D為風(fēng)輪直徑)，風(fēng)力機(jī)置于距離計(jì)算域入口3D處。為兼顧計(jì)算效率與精度，同時(shí)考慮到葉片表面扭曲復(fù)雜，網(wǎng)格劃分采用混合網(wǎng)格離散形式，將整個(gè)計(jì)算域劃分為局部加密區(qū)域和外圍區(qū)域。局部加密區(qū)域內(nèi)含風(fēng)力機(jī)模型，采用非結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分，外圍區(qū)域形狀規(guī)整，采用高質(zhì)量的結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分。

為確保數(shù)值模擬結(jié)果的可信度，本文增加了網(wǎng)格無關(guān)性驗(yàn)證，表4給出了不同網(wǎng)格方案下網(wǎng)格質(zhì)量和迎風(fēng)面壓力系數(shù)。由表4可知，隨著網(wǎng)格總數(shù)的增加，網(wǎng)格質(zhì)量逐漸增加，網(wǎng)格歪斜度和迎風(fēng)面風(fēng)壓系數(shù)呈現(xiàn)逐漸減小的趨勢(shì)，而840萬網(wǎng)格數(shù)和1 100萬網(wǎng)格數(shù)的網(wǎng)格質(zhì)量和計(jì)算結(jié)果無明顯差異，綜合計(jì)算精度和效率，本文選取840萬網(wǎng)格總數(shù)的方案。計(jì)算域及具體的網(wǎng)格劃分如圖1所示。

表4 不同網(wǎng)格方案下網(wǎng)格質(zhì)量和迎風(fēng)面壓力系數(shù)

圖1 計(jì)算域及網(wǎng)格劃分示意圖Fig.1 Schematic diagram of computational domain and encrypted mesh

計(jì)算域左側(cè)和頂部邊界條件為速度入口，右側(cè)為壓力出口，數(shù)值計(jì)算采用3D單精度、分離式求解器，流場(chǎng)流速為絕對(duì)速度，空氣模型等效為理想不可壓縮流體，計(jì)算域入口采用冪指數(shù)為0.15的風(fēng)廓線模型，離地10 m高度處的風(fēng)速分別設(shè)置為“2.2”節(jié)中3種基準(zhǔn)風(fēng)速。流場(chǎng)求解采用SIMPLEC算法實(shí)現(xiàn)速度與壓力之間的耦合，對(duì)流項(xiàng)求解格式為二階，計(jì)算過程中設(shè)置了網(wǎng)格傾斜校正以提高混合網(wǎng)格計(jì)算效果，控制方程的計(jì)算殘差設(shè)置為1×10-6，最后初始化風(fēng)場(chǎng)進(jìn)行迭代計(jì)算。圖2給出了平均風(fēng)速、湍流度剖面模擬結(jié)果與理論值的對(duì)比曲線，結(jié)果表明平均風(fēng)速和湍流度剖面均與理論值吻合良好，風(fēng)場(chǎng)模擬標(biāo)準(zhǔn)滿足工程要求。

圖2 不同風(fēng)速下風(fēng)場(chǎng)速度及湍流度剖面示意圖Fig.2 Velocity and turbulence profiles under different wind speeds

2.4 有效性驗(yàn)證

圖3給出了3種風(fēng)速下塔架30 m高度環(huán)向壓力系數(shù)分布圖。對(duì)比可知各工況平均風(fēng)壓系數(shù)分布曲線的負(fù)壓極值點(diǎn)和分離點(diǎn)對(duì)應(yīng)角度與規(guī)范曲線[18]一致，迎風(fēng)和側(cè)風(fēng)區(qū)域風(fēng)壓系數(shù)數(shù)值吻合較好，僅在背風(fēng)區(qū)負(fù)壓系數(shù)略大于規(guī)范值，對(duì)比結(jié)果驗(yàn)證了風(fēng)場(chǎng)數(shù)值模擬的有效性。

圖3 3種風(fēng)速下風(fēng)力機(jī)塔架30 m高度壓力系數(shù)與規(guī)范對(duì)比曲線Fig.3 Pressure coefficients on the 30 m height tower under three wind speeds and domestic codes

3 結(jié)果分析

為研究風(fēng)-沙荷載特征值隨環(huán)境參數(shù)變化的作用規(guī)律，首先分析風(fēng)-沙雙向耦合作用下風(fēng)場(chǎng)分布特性、沙粒著點(diǎn)及其分布規(guī)律；在此基礎(chǔ)上，通過計(jì)算沙粒沖擊荷載、等效外壓系數(shù)定量分析不同工況下沙致效應(yīng)差異。

3.1 風(fēng)場(chǎng)分布特性

圖4給出了不同風(fēng)速下大型風(fēng)力機(jī)體系風(fēng)速流線圖。由圖4可知，由于葉片的遮擋效應(yīng)，塔架迎風(fēng)面壓力明顯減小，葉片與塔架重合部位內(nèi)側(cè)出現(xiàn)附著于塔架迎風(fēng)面的自由剪切層，并在塔架背風(fēng)面形成成熟的旋渦脫落。

圖4 風(fēng)力機(jī)體系風(fēng)速流線圖Fig.4 Wind speed streamline of wind turbine

圖5和圖6分別給出了未干擾及干擾區(qū)段塔架渦量云圖。對(duì)比分析可知：未干擾區(qū)段氣流在風(fēng)力機(jī)塔架迎風(fēng)面前緣發(fā)生流動(dòng)分離，并在側(cè)邊出現(xiàn)加速效應(yīng)；而在干擾區(qū)段由于葉片對(duì)塔架的干擾作用，塔架和葉片之間的區(qū)域隨著風(fēng)速的增加出現(xiàn)明顯的能量積聚，持續(xù)發(fā)展后在背風(fēng)面形成尾流渦旋以及回流，且隨風(fēng)速的增加而加強(qiáng)。

圖5 風(fēng)力機(jī)塔架未干擾段湍動(dòng)能圖Fig.5 Wind turbine tower undisturbed turbulent kinetic energy

圖6 風(fēng)力機(jī)塔架干擾段湍動(dòng)能圖Fig.6 Turbulence kinetic energy of wind turbine tower interference section

3.2 風(fēng)沙運(yùn)動(dòng)特征

風(fēng)沙流場(chǎng)中沙粒沖擊風(fēng)力機(jī)表面，易產(chǎn)生極端荷載效應(yīng)。為清晰展示沙粒在風(fēng)力機(jī)不同部位的分布密度，基于顆粒水平向速度對(duì)沙粒軌跡進(jìn)行追蹤，同時(shí)定義單位時(shí)間內(nèi)沖擊塔架沙粒的位置分布為沙粒著點(diǎn)，圖7給出了9種工況風(fēng)沙場(chǎng)中沙粒著點(diǎn)圖，并對(duì)沙粒的密集程度進(jìn)行了等比例粗化處理。

由圖7可知，由于葉片的對(duì)塔架的遮擋效應(yīng)，各工況沙粒撞擊位置多集中分布風(fēng)力機(jī)塔架迎風(fēng)區(qū)域0～0.6H高度范圍內(nèi)，H為風(fēng)力機(jī)塔架高度。受氣流旋渦驅(qū)動(dòng)作用，塔架遮擋區(qū)和背風(fēng)區(qū)壁面僅有少量沙粒附著；塔架、葉片收集到的沙粒均以工況7最多，且隨風(fēng)速的增加而增加，隨沙粒粒徑的增大逐漸變小。

圖7 風(fēng)沙場(chǎng)中沙粒著點(diǎn)示意圖Fig.7 Sand-reflect motion trajectory in sand and wind field

3.3 風(fēng)沙荷載分布

3.3.1 沙粒空間分布規(guī)律

圖8和圖9分別給出了9種工況下風(fēng)力機(jī)塔架和葉片收集到的沙粒數(shù)量、撞擊速度及速度占有率(各速度分布區(qū)間內(nèi)沙粒數(shù)量與總數(shù)量的比值)對(duì)比曲線。由圖8和圖9對(duì)比可知：

(1) 當(dāng)沙粒粒徑較小時(shí)，風(fēng)力機(jī)塔架收集到的沙粒數(shù)量隨風(fēng)速的增大而增加，當(dāng)粒徑達(dá)到中等以上時(shí)，風(fēng)力機(jī)塔架收集到的沙粒數(shù)量隨風(fēng)速的增大而減少；然而，粒徑和風(fēng)速對(duì)風(fēng)力機(jī)葉片的沙粒著點(diǎn)影響并不顯著，這是由于葉片位置較高，處于風(fēng)沙流密度影響較小區(qū)域。

圖8 風(fēng)力機(jī)塔架沙粒數(shù)量與水平末速度分布曲線Fig.8 Wind turbine tower sand number and horizontal velocity distribution curve

圖9 風(fēng)力機(jī)葉片沙粒數(shù)量與水平末速度分布曲線Fig.9 Wind turbine blade grit and horizontal velocity distribution curve

(2) 撞擊風(fēng)力機(jī)塔架的沙粒，隨著風(fēng)速增大，小粒徑和大粒徑沙粒的平均速度均增大，而中等粒徑的沙粒平均速度隨風(fēng)速的增大而減??；風(fēng)力機(jī)葉片上追蹤的沙粒平均速度特性也呈現(xiàn)相似規(guī)律，不同在于風(fēng)力機(jī)葉片上的大粒徑沙粒平均速度隨風(fēng)速的增大而減小。

(3) 風(fēng)力機(jī)塔架、葉片收集到的沙粒水平速度占有率分布規(guī)律基本一致，速度占有率隨著水平末速度(順風(fēng)向平均速度)的增大先增加后減小。此外，收集到的沙粒平均水平速度均遠(yuǎn)小于基準(zhǔn)風(fēng)速，且沙粒水平撞擊速度均隨著沙粒直徑的增加而增大。

3.3.2 葉片、塔架沙粒沖擊荷載

定義與塔架重合的葉片為A、迎風(fēng)面左側(cè)葉片為B、右側(cè)葉片為C。采用式(12)分別進(jìn)行9種工況下風(fēng)力機(jī)塔架、葉片沙荷載計(jì)算，同時(shí)給出了塔架、葉片不同高度(葉展長(zhǎng)度)范圍沙荷載以及風(fēng)力機(jī)不同高度范圍沙荷載與該區(qū)域風(fēng)荷載的比值，如圖10、圖11和表5所示。經(jīng)對(duì)比發(fā)現(xiàn)：

(1) 大粒徑沙粒與不同風(fēng)速的組合工況(工況3、工況6、工況9)在風(fēng)力機(jī)表面產(chǎn)生顯著荷載，由于位置和形狀的不同，風(fēng)力機(jī)塔架、葉片沙荷載分布呈現(xiàn)不同的規(guī)律；

(2) 葉片A、葉片B、葉片C沙荷載均沿葉展方向呈減小趨勢(shì)，由于所處高度不同，風(fēng)沙流密度影響程度不同，葉片A承受最大沙荷載約為葉片B、葉片C的兩倍；

(3) 不同工況對(duì)風(fēng)力機(jī)塔架表面沙荷載分布影響較大，但均在0～0.05H高度范圍內(nèi)產(chǎn)生該工況下極大沙荷載，工況3下該高度范圍內(nèi)沙荷載與風(fēng)荷載比值達(dá)22.57%；

(4) 顯著工況(工況3、工況6、工況9)下風(fēng)力機(jī)塔架極大沙荷載隨風(fēng)速的變大而增大，沙荷載沿塔架高度先減小后增大再減小，0.2H～0.6H高度范圍內(nèi)沙荷載極大值也隨風(fēng)速的增大而增大。由于葉片A的遮擋效應(yīng)，塔架0.5H高度以上沙荷載減小。

圖10 各工況風(fēng)力機(jī)葉片沙荷載示意圖Fig.10 Sand load distribution on different height of wind turbine blades

圖11 風(fēng)力機(jī)塔筒沿子午線高度荷載示意圖Fig.11 Comparison curves of equivalent pressure in typical meridian of wind turbine tower

表5 工況1～工況9塔架不同高度范圍風(fēng)沙荷載比值列表

3.3.3 沙壓系數(shù)

沙致壓力系數(shù)(簡(jiǎn)稱沙壓系數(shù))計(jì)算方法為：①將表面各監(jiān)控點(diǎn)沙粒撞擊荷載轉(zhuǎn)化成沙致壓強(qiáng)；②計(jì)算監(jiān)控點(diǎn)沙致壓強(qiáng)與對(duì)應(yīng)參考高度處風(fēng)壓比值，即沙壓系數(shù)。圖12給出了各工況下風(fēng)力機(jī)塔架沙壓系數(shù)三維分布示意圖。由圖12可知：

(1) 各工況風(fēng)力機(jī)塔架沙壓系數(shù)均主要集中于迎風(fēng)面兩側(cè)各60°范圍內(nèi)，其余范圍數(shù)值基本為0，各工況沙壓系數(shù)最大值均發(fā)生在0～0.05H高度范圍內(nèi)，各工況沙壓系數(shù)最大值為0.517；

(2) 沙壓系數(shù)分布呈現(xiàn)明顯三維效應(yīng)，其隨塔架高度增加而迅速減小，塔架與葉片重合部分沙壓系數(shù)趨于0；

(3) 不同組合工況下沙壓系數(shù)大小差異明顯，且沙粒粒徑的增大而顯著增大，此外，其最大值隨風(fēng)速的增大呈減小趨勢(shì)，而其總體分布更趨于平均化。

圖12 各工況風(fēng)力機(jī)塔架沙致壓力系數(shù)三維分布示意圖Fig.12 Three-dimensional distribution of sand-reflect and sand pressure coefficient on surface of wind turbine tower

3.3.4 風(fēng)-沙共同作用等效壓力系數(shù)

為定量比較不同組合工況下的塔架風(fēng)-沙致壓力分布，將沙壓系數(shù)與風(fēng)壓系數(shù)矢量加和，定義等效壓力系數(shù)作為沙致效應(yīng)的衡量參數(shù)。

考慮沙粒著點(diǎn)分布特性及葉片與塔筒之間的氣動(dòng)干擾效應(yīng)，選取塔架未干擾區(qū)段沙粒分布密度差異較大，干擾區(qū)段塔架與葉片重合面積不同的兩個(gè)斷面為典型斷面。圖13給出各工況下風(fēng)力機(jī)塔架4個(gè)典型斷面等效壓力系數(shù)對(duì)比曲線，分析可知：

圖13 風(fēng)力機(jī)塔架典型斷面環(huán)向等效壓力系數(shù)對(duì)比曲線Fig.13 Comparison curves of equivalent pressure coefficient in typical section of wind turbine tower

(1) 不同工況下相同高度截面處等效壓力系數(shù)分布規(guī)律及數(shù)值基本一致，位于未干擾區(qū)(0.025H，0.30H)塔筒斷面環(huán)向風(fēng)壓系數(shù)呈現(xiàn)良好的對(duì)稱性；而位于顯著干擾區(qū)(0.70H，0.90H)塔筒斷面，其環(huán)向表面壓力系數(shù)分布曲線不再保持對(duì)稱，塔筒迎風(fēng)面中心點(diǎn)處呈現(xiàn)負(fù)壓；

(2) 同一工況下4個(gè)典型斷面等效壓力系數(shù)數(shù)值均略有差異，其中最大負(fù)壓隨著高度的增加先增大后減小，在葉片遮擋區(qū)域附近約為-1.5；

(3) 不同工況下0.025H高度塔筒斷面迎風(fēng)面兩側(cè)0～30°內(nèi)壓力系數(shù)差異明顯，工況3、工況6、工況9大于其余工況下等效壓力系數(shù)，且最大壓力系數(shù)隨風(fēng)速的增大而減小，背風(fēng)面呈負(fù)壓；

(4) 位于干擾區(qū)段塔架斷面，隨高度增加迎風(fēng)面、側(cè)風(fēng)面和背風(fēng)面負(fù)壓分別呈增大、減小和增大的趨勢(shì)。

圖14給出了風(fēng)力機(jī)塔架0°，90°，180°及270°四條典型子午線的等效壓力系數(shù)對(duì)比曲線，分析可得：

(1) 風(fēng)力機(jī)塔架90°與270°子午線等效壓力系數(shù)隨高度先增大后減小再增大，均在0.5H，0.8H高度處出現(xiàn)拐點(diǎn)；0°和180°子午線迎風(fēng)面等效壓力系數(shù)發(fā)生數(shù)值分離多集中在沙粒著點(diǎn)數(shù)目較多的0.5H高度以下的未干擾區(qū)段，0.1H～0.6H高度范圍內(nèi)背風(fēng)面等效壓力系數(shù)出現(xiàn)明顯差異；

(2) 90°和270°子午線等效壓力系數(shù)最大值均產(chǎn)生在0.5H高度處，最大值約為為-1.890；受風(fēng)沙流影響，0°子午線最大等效壓力系數(shù)發(fā)生0.005H高度處，約為1.092；270°子午線方向最大等效壓力系數(shù)發(fā)生在塔架頂部附近，最大值為-0.282；

(3) 0°子午線底部等效壓力系數(shù)受風(fēng)沙流影響顯著，大粒徑沙粒對(duì)塔架的正面沖擊作用加強(qiáng)，沙粒著點(diǎn)隨風(fēng)速的增大而減少，故等效壓力系數(shù)隨著風(fēng)速的增大而減小。

圖14 風(fēng)力機(jī)塔筒典型子午線等效壓力系數(shù)對(duì)比曲線Fig.14 Comparison curves of equivalent pressure coefficient in typical meridians of wind turbine tower

4 結(jié) 論

(1) 僅考慮單向風(fēng)作用時(shí)風(fēng)力機(jī)體系風(fēng)場(chǎng)和風(fēng)壓系數(shù)均呈現(xiàn)明顯的三維效應(yīng)，考慮風(fēng)-沙雙向耦合后沙粒主要沖擊塔架迎風(fēng)面中下部，僅有少量沙粒隨氣流進(jìn)入塔筒背風(fēng)面；沙粒水平方向作用力與分布范圍大小受風(fēng)速和粒徑的影響，在塔架和葉片等不同部位呈現(xiàn)不同的規(guī)律。

(2) 風(fēng)力機(jī)塔架、葉片沙粒沖擊數(shù)量隨風(fēng)速的增大而增加，隨粒徑的增大而減少，沙粒統(tǒng)計(jì)數(shù)量最多的為工況7(特強(qiáng)沙塵暴發(fā)生風(fēng)速+典型工況下較小沙粒粒徑)。塔架上的沙粒速度占有率最大值隨風(fēng)速、粒徑的增大而降低，速度分布更加均勻，而沖擊葉片沙粒速度占有率分布規(guī)律不受風(fēng)速、粒徑等因素干擾，主要表現(xiàn)為峰值右移。

(3) 不同工況下沙粒撞擊位置主要集中在塔架0.6H高度以下，受兩側(cè)流動(dòng)分離渦裹挾，集中于迎風(fēng)區(qū)域兩側(cè)各60°范圍，且大粒徑沙粒沖擊塔架表面更易形成顯著荷載，沙致壓力系數(shù)最大值為0.517，發(fā)生在工況3(中等沙塵暴發(fā)生風(fēng)速+極端工況下沙粒粒徑)的塔架底部迎風(fēng)面。

(4) 不同工況下風(fēng)力機(jī)塔架最大負(fù)壓值隨著高度的增加先增大后減小，背風(fēng)區(qū)域負(fù)壓值隨著高度的增加先減小后增大，在葉片遮擋區(qū)域附近約為-1.5。各子午線等效壓力系數(shù)隨著高度的增加變化趨勢(shì)不同，0°子午線底部等效壓力系數(shù)受風(fēng)沙流影響顯著，大粒徑沙粒對(duì)塔架的正面沖擊作用加強(qiáng)，最大等效壓力系數(shù)發(fā)生0.005H高度處，數(shù)值為1.092。