王 巖，董 雯，李成振

(1.長春工程學院，長春 130012； 2.吉林省水文水資源局，長春 130000； 3.中水東北公司科研院，長春 130000)

0 引言

風能是資源潛力大、技術(shù)基本成熟的可再生能源，在減排溫室氣體、應(yīng)對氣候變化的新形勢下，越來越受到世界各國的重視，并已在全球大規(guī)模開發(fā)利用。到2015年底，全球累計風電裝機容量約為2億kW，最近5年年平均增長率約為30%左右，隨著風電的大規(guī)模開發(fā)，風電設(shè)備正朝著特性化和大型化方向發(fā)展，2～3 MW風電機組已是市場主流產(chǎn)品，5 MW風機也已經(jīng)開始應(yīng)用。

隨著風電的發(fā)展，葉片展長在增加，而葉片作為主要的能量轉(zhuǎn)換載體，其使用壽命就顯得極為重要，因此，需要對不同工作狀態(tài)下的葉片進行性能分析。除了借助風洞試驗的方法外，隨著高性能計算機的普及，采用現(xiàn)有軟件進行數(shù)值模擬，可以克服試驗條件不足和實體葉片截面復(fù)雜的問題，具有一定的現(xiàn)實意義。

但虎在ANSYS中建立復(fù)合材料模型，并對葉片模態(tài)進行分析計算；張鎖懷等借助Solidwork平臺，建立MW級風力發(fā)電機的三維幾何模型；衣傳寶等在研究風電機組變槳距控制時在Simulink里建立數(shù)學模型。

綜合查閱文獻后得出目前風機建模存在有以下兩種問題：1)在有限元軟件中建模機組運轉(zhuǎn)困難或收斂較慢；2)以變槳距為主要研究對象的建模，忽略了葉片和塔架的柔性結(jié)構(gòu)，只進行簡化建模。

針對上述情況，筆者通過有限元軟件abaqus建立葉片及塔架的實體模型，在adams里建立齒輪箱、傳動軸、變槳機構(gòu)等，并在adams中進行整體組裝，建立完整的風機模型，實現(xiàn)柔性機構(gòu)設(shè)計與機電控制的耦合。

1 風機建?；緟?shù)

以目前應(yīng)用廣泛的Vestas V80型2 MW風機為例進行模擬，風機基本參數(shù)為：風機功率調(diào)節(jié)類型為變速變槳距；額定功率2 MW；切入風速4 m/s；額定風速16 m/s；切出風速25 m/s；齒輪箱為兩級行星齒輪和一級斜齒輪；機艙質(zhì)量6.8×104kg；葉輪質(zhì)量3.8×104kg；葉片展長39 m；單葉片質(zhì)量5.05 t；額定轉(zhuǎn)速16.7 r/min，葉片最大弦長2.889 m；塔架高度78 m。

2 葉片參數(shù)及實體建模

風力機葉片作為主要的能量捕捉裝置，其性能大大影響風能的利用效率。風力機葉片翼型是由飛機機翼發(fā)展而來的，因此，有標準翼型供參考。早期風力機葉片設(shè)計首選比較成熟、升阻特性較好的傳統(tǒng)航空翼型，隨著葉片向大而長的方向發(fā)展，傳統(tǒng)翼型的劣勢逐漸突出。目前，國外各風電生產(chǎn)企業(yè)已經(jīng)加大對新翼型的研發(fā)力度，而國內(nèi)翼型研究才剛剛起步，而且具體數(shù)據(jù)極難獲得，筆者主要針對整機結(jié)構(gòu)性能進行分析，因此，翼型研究適當簡化。

2.1 實體建模

目前，葉片建模常采用NACA4412、NACA23012等兩種常見的翼型。本次建模采用NACA4412翼型，使得葉片具有良好的氣動性能，葉片截面翼型數(shù)據(jù)通過profili軟件得到，具體參數(shù)見表1。設(shè)計過程中，根據(jù)葉片沿翼展方向的弦長數(shù)據(jù)，對NACA4412單位參數(shù)進行縮放。

表1 翼型上下弦單位參數(shù)

葉片展長39 m，最大弦長2.889 m，葉片尖部弦長0.05 m，葉片沿翼展方向分成24個葉素，不同位置剖面需依據(jù)實際弦長對基本參數(shù)進行縮放，葉片根部為同槳距軸承配合須做部分修改，使之漸變?yōu)閳A形。建立相應(yīng)的空間三維坐標，排序并生成能以keypoint型式導入ansys的.txt文件。在ansys軟件中導入坐標點，依實際線形將相鄰序號連成線，形成翼型上弧線和下弧線，將相鄰弧線生成曲面，依照點成線、線成面的方式，建立三維模型，具體模型如圖1所示。

圖1 葉片三維模型

2.2 材料參數(shù)設(shè)置

整個葉片在徑向和軸向的鋪層結(jié)構(gòu)和厚度存在不同。大型風力機葉片多采用玻璃纖維復(fù)合材料或碳/?；祀s復(fù)合材料，同時加設(shè)防雷電網(wǎng)，聚氨酯膠衣等，蒙皮由雙向玻纖織物增強，空腔內(nèi)加設(shè)承剪腹板。為簡化建模過程，材料設(shè)置為展向彈性模量24.2 MPa，徑向彈性模量9.87 MPa，泊松比0.39，單葉片質(zhì)量5.05 t。

2.3 網(wǎng)格劃分

在abaqus里進行網(wǎng)格劃分，對于殼單元可劃分為四邊形和三角形。四邊形網(wǎng)格劃分對于局部扭曲面網(wǎng)格劃分過大，不滿足計算要求，因此劃分為三角形，采用三結(jié)點三角形通用殼S3單元。由于整機模型，自有度較多，網(wǎng)格劃分過細會產(chǎn)生龐大的計算工作量，因此，在滿足精度的基礎(chǔ)上，可以加大網(wǎng)格劃分，邊緣及全局網(wǎng)格尺寸控制在0.5 m。網(wǎng)格劃分后的葉片如圖2所示。

圖2 葉片網(wǎng)格劃分

3 塔架建模

塔架是重要的承重結(jié)構(gòu)，隨著風力機單機容量的增加，與之配套的圓筒形塔架也向高聳化發(fā)展。當風機運行時，隨著風荷載、機組自重、機組振動激勵等的影響，塔架頂會產(chǎn)生較大的位移，其剛度直接影響著風力機組的穩(wěn)定性。

塔架除應(yīng)滿足對機組的支承作用，還設(shè)置一些輔助設(shè)備，入門洞、爬梯、通風口等，過于復(fù)雜的結(jié)構(gòu)會降低仿真步長，延長仿真時間，因此需對結(jié)構(gòu)進行適當簡化。

1)建模及材料設(shè)置：本文研究對象為變截面錐桶型塔架，78 m高，材料為Q345鋼，底部直徑4.23 m，頂端直徑2.98 m，頂部壁厚18 mm，底部壁厚28 mm，彈性模量2.06×105MPa，泊松比0.3，密度7 800 kg/m3。鑒于塔架壁厚逐段變化，將塔架從上至下分為4部分獨立建模：第一段長度33 m，鋼管厚度18 mm；第二段長度25 m，鋼管厚度23 mm；第三段長度20 m，底部基礎(chǔ)20 m，厚度均為28 mm。在有限元軟件中建模，具體模型如圖3。

2)網(wǎng)格劃分：塔架形狀較規(guī)則，適合劃分為四邊形網(wǎng)格。網(wǎng)格尺寸控制在1.0 m以內(nèi)，單元類型為線性、有限薄膜應(yīng)變、四邊形殼單元S4R。

圖3 塔架網(wǎng)格剖分

4 傳動鏈及變槳機構(gòu)建模

為模擬風機在不同工作狀態(tài)下的工作特性，需建立傳動鏈及變槳距模型。風力發(fā)電機的運行可分為3個階段，當風速介于切入風速和額定風速之間時，為保證最優(yōu)葉尖速比以獲得最大風能，發(fā)電機通過齒輪箱實現(xiàn)對葉輪旋轉(zhuǎn)速度的控制，為最佳葉尖速比運行區(qū)；當風速高于額定風速時，首先進入恒轉(zhuǎn)速運行區(qū)，葉輪轉(zhuǎn)速保持恒定；當風速繼續(xù)加大，需減小葉片槳距角，使風能利用率降低，從而實現(xiàn)輸出功率的恒定。

建模所用齒輪箱傳動比為100.5∶1。在adams所建模型分為兩級傳動，直徑比分別為10∶1和10.05∶1，利用圓柱體代替實際齒輪，通過建立齒輪副使齒輪咬合點產(chǎn)生相同速度，實現(xiàn)兩級齒輪的連動，此時右側(cè)齒輪旋轉(zhuǎn)時，左側(cè)發(fā)電機以傳動比放大后的轉(zhuǎn)速進行旋轉(zhuǎn)。

在adams軟件中，柔性體葉片與剛性輪轂之間設(shè)置傘狀剛性過渡圓盤及連接軸，連接軸與輪轂之間位置固定，并能以自身軸線為中心帶動葉片旋轉(zhuǎn)，從而實現(xiàn)葉片槳距角的變化。

圖4 齒輪箱模型

5 整機模型建立

大功率風力發(fā)電機的葉片及塔架呈細長狀，運行過程中變形明顯，屬柔性體。adams軟件側(cè)重于剛性構(gòu)件間的相對運動，因此，建立組合模型前，葉片及塔架應(yīng)在有限元軟件中進行模態(tài)分析并生成可以讀入adams軟件中的MNF文件，以便進行構(gòu)建間的組合安裝。進行模態(tài)分析時，葉片及塔架取前六階模態(tài)，塔架基礎(chǔ)在生成MNF文件前，應(yīng)將基礎(chǔ)做固定約束。組合模型如圖5所示。

圖5 風力發(fā)電機組合模型

6 結(jié)語

1)通過對原型機進行簡化，建立傳動軸模型及變槳距模型；建立葉片及塔架的有限元模型，并生成可導入adams的MNF文件。

2)集有限元軟件及多體動力學軟件的優(yōu)點于一體的組合模型，避免有限元軟件中模擬時動態(tài)速度過慢的缺點，以及采取數(shù)學建模時，葉片過于復(fù)雜，簡化困難的缺點。

3)便于開展氣動分析，軟件adams與風荷載軟件aerodyn有良好接口，便于開展接續(xù)工作。

4)便于實現(xiàn)變速及變槳距的機電控制，adams軟件與simulink具有兼容的接口，可以實現(xiàn)基于simulink的轉(zhuǎn)輪轉(zhuǎn)速及槳距角的實時控制。

此次剛?cè)峄旌辖?，也便于對運行中的風力機進行進一步動態(tài)力學性能研究。