孟令銳，何海建，楊 揚(yáng)，祝蘊(yùn)龍，晁貫良，董姝言

(許昌許繼風(fēng)電科技有限公司，河南 許昌 461000)

0 引 言

風(fēng)力發(fā)電機(jī)(以下簡(jiǎn)稱風(fēng)機(jī))主要由葉片、齒輪箱、發(fā)電機(jī)、塔筒、基礎(chǔ)等部分組成。塔筒作為風(fēng)機(jī)的重要組成部分，承擔(dān)著風(fēng)機(jī)的大部分重量，并為風(fēng)機(jī)運(yùn)行提供足夠強(qiáng)度、剛度等安全保障[1-2]。塔筒門(mén)框(以下簡(jiǎn)稱門(mén)框)作為塔筒的一部分，位于塔筒底部靠近基礎(chǔ)部分，主要供人員和設(shè)備進(jìn)出。由于門(mén)框的存在改變了塔筒原有結(jié)構(gòu)，使門(mén)框附近形狀發(fā)生突變，存在應(yīng)力集中現(xiàn)象，塔筒門(mén)框成為塔筒最易受到破壞的部位之一，因此很有必要對(duì)塔筒門(mén)框進(jìn)行優(yōu)化研究。

目前對(duì)于塔筒門(mén)框的分析，主流方法是有限元法。任瑞杰等[3]主要分析了塔筒門(mén)框的筋板厚度對(duì)門(mén)框的極限應(yīng)力及其出現(xiàn)位置的影響，但僅選取了幾個(gè)代表尺寸，忽略了筋板向內(nèi)外拉伸厚度的影響。晏紅文等[4]主要分析了門(mén)框形狀參數(shù)、門(mén)框離地高度等對(duì)門(mén)框應(yīng)力的影響，并對(duì)其進(jìn)行了優(yōu)化，但未考慮門(mén)框尺寸對(duì)其應(yīng)力的影響。

本文主要基于GL2010規(guī)范[5]分析門(mén)框尺寸變化對(duì)其應(yīng)力的影響，并對(duì)其進(jìn)行優(yōu)化分析。

1 有限元分析

1.1 有限元模型

塔筒模型坐標(biāo)系根據(jù)GL2010規(guī)范進(jìn)行確定。坐標(biāo)原點(diǎn)在塔筒底部法蘭下底面中心處，X軸由原點(diǎn)指向順風(fēng)向，Y軸由原點(diǎn)指向順風(fēng)向的左方，z軸由原點(diǎn)垂直指向塔筒上方，門(mén)框的開(kāi)口位置為塔筒迎風(fēng)面的正前方。

塔筒內(nèi)部平臺(tái)、爬梯、電纜支架等附屬件對(duì)門(mén)框結(jié)構(gòu)分析影響不大，在建模時(shí)對(duì)其進(jìn)行簡(jiǎn)化處理。為了提高計(jì)算效率，根據(jù)圣維南原理[6]，本次計(jì)算僅建立包含門(mén)框的塔筒底段模型。

根據(jù)某3MW風(fēng)力發(fā)電機(jī)組塔筒圖紙，本研究在Workbench軟件中建立塔筒底段的有限元模型。整個(gè)模型采用高階六面體單元進(jìn)行網(wǎng)格劃分，在門(mén)框及其附近區(qū)域?qū)W(wǎng)格進(jìn)行細(xì)化處理，門(mén)框以外的區(qū)域適當(dāng)?shù)販p小網(wǎng)格的密度。

門(mén)框及其附近區(qū)域網(wǎng)格劃分如圖1所示。

圖1 門(mén)框及其附近區(qū)域網(wǎng)格劃分示意圖

1.2 材料性能

塔筒與門(mén)框均采用Q345E結(jié)構(gòu)鋼焊接成型，其彈性模量為2.1×1011Pa，泊松比為0.3，密度為7 850 kg/m3。根據(jù)GB/T 1591-2008《低合金高強(qiáng)度結(jié)構(gòu)鋼》[7]，材料屬性如表1所示(根據(jù)GL2010規(guī)范，在進(jìn)行零部件靜強(qiáng)度計(jì)算時(shí)材料的安全系數(shù)應(yīng)取為1.1)。

表1 材料Q345屈服強(qiáng)度表

1.3 邊界條件

工程實(shí)踐表明門(mén)框強(qiáng)度對(duì)彎矩載荷較為敏感，因此為了計(jì)算方便僅考慮載荷分量對(duì)門(mén)框及其附近區(qū)域應(yīng)力分布的影響。施加完邊界條件的塔筒模型如圖2所示。

圖2 施加完邊界條件的塔筒模型

本研究根據(jù)GL2010規(guī)范，計(jì)算得到極限工況下塔筒頂部為1.934E8 Nm，載荷施加在塔筒頂部中心點(diǎn)，如圖2中B點(diǎn)所示。對(duì)塔筒底部進(jìn)行全約束，如圖2中A點(diǎn)所示。

1.4 門(mén)框強(qiáng)度計(jì)算結(jié)果

經(jīng)計(jì)算門(mén)框及其附近區(qū)域的最大等效應(yīng)力為280.83 MPa，位于門(mén)框下底部門(mén)框筋板與筒壁結(jié)合處，應(yīng)力云圖如圖3所示。

圖3 門(mén)框等效應(yīng)力云圖

查表1可知當(dāng)部件厚度為88 mm時(shí)，材料的屈服強(qiáng)度為305 MPa。材料極限強(qiáng)度安全系數(shù)計(jì)算公式如下：

(1)

式中:γf—載荷安全系數(shù)，取1.0；γm—材料安全系數(shù)，取1.1；σeqvmax—等效應(yīng)力的最大值，MPa；Re—材料屈服強(qiáng)度，MPa。

經(jīng)計(jì)算得到門(mén)框?yàn)?.987，無(wú)法滿足設(shè)計(jì)要求，因此需要對(duì)門(mén)框結(jié)構(gòu)進(jìn)行優(yōu)化來(lái)確保門(mén)框及其附近區(qū)域應(yīng)力滿足設(shè)計(jì)要求。

2 計(jì)算數(shù)學(xué)模型優(yōu)化

要建立優(yōu)化計(jì)算的數(shù)學(xué)模型需要先確定模型的設(shè)計(jì)變量、約束條件、優(yōu)化目標(biāo)。

2.1 設(shè)計(jì)變量

影響門(mén)框及其附近區(qū)域應(yīng)力變化的因素眾多，例如門(mén)框形狀、門(mén)框離地高度、筒壁厚度等，本文僅假設(shè)在其他條件均不變的情況下，研究門(mén)框筋板厚度t1、筋板向外拉伸厚度h1、筋板向內(nèi)拉伸厚度h2三者對(duì)門(mén)框及其附近區(qū)域應(yīng)力的影響。因此本研究將t1、h1、h2設(shè)為優(yōu)化設(shè)計(jì)變量，記為X=(t1,h1,h2)T。

各變量示意圖如圖4所示。

圖4 門(mén)框設(shè)計(jì)變量示意圖

各設(shè)計(jì)變量的優(yōu)化取值范圍如表2所示。

表2 設(shè)計(jì)變量取值范圍

2.2 目標(biāo)函數(shù)與約束條件

目標(biāo)函數(shù)是判斷優(yōu)化方案優(yōu)劣的數(shù)學(xué)表達(dá)式，通過(guò)找出目標(biāo)函數(shù)的極值，從而找出最優(yōu)的方案[8]。因門(mén)框的靜強(qiáng)度無(wú)法滿足設(shè)計(jì)要求，需要對(duì)塔筒門(mén)框及其附近區(qū)域的最大等效應(yīng)力進(jìn)行優(yōu)化。同時(shí)在確保塔筒門(mén)框強(qiáng)度滿足設(shè)計(jì)要求的前提下，為了降低門(mén)框的成本，也要對(duì)門(mén)框的質(zhì)量進(jìn)行優(yōu)化。

由于門(mén)框筋板厚度在優(yōu)化時(shí)是變化的，本研究為了保守起見(jiàn)，取其最大厚度120 mm計(jì)算得到其許用應(yīng)力為259 MPa。優(yōu)化目標(biāo)是使門(mén)框及其附近區(qū)域的最大等效應(yīng)力在不大于259 MPa情況下，盡量減小門(mén)框的等效應(yīng)力和門(mén)框的質(zhì)量。因此在優(yōu)化時(shí)應(yīng)將等效應(yīng)力優(yōu)化優(yōu)先級(jí)設(shè)為最高，門(mén)框質(zhì)量的優(yōu)化優(yōu)先級(jí)次之。

設(shè)等效應(yīng)力的優(yōu)化目標(biāo)函數(shù)為F(X)，質(zhì)量的優(yōu)化目標(biāo)函數(shù)為G(X)。

2.3 優(yōu)化模型的數(shù)學(xué)表達(dá)式

構(gòu)造好的優(yōu)化模型數(shù)學(xué)表達(dá)式如下：

(2)

式中：F0(X)—門(mén)框許用等效應(yīng)力；G0(X)—門(mén)框原始質(zhì)量；Fi(X)—應(yīng)力函數(shù)響應(yīng)值，i=1,2,……n；Gj(X)—質(zhì)量函數(shù)響應(yīng)值，j=1,2,……n；n—樣本個(gè)數(shù)。

3 響應(yīng)面分析

響應(yīng)面法(response surface method， RSM)是通過(guò)對(duì)設(shè)計(jì)點(diǎn)的試驗(yàn)分析得到目標(biāo)函數(shù)在約束條件下的響應(yīng)面模型，并以此來(lái)預(yù)測(cè)非設(shè)計(jì)點(diǎn)響應(yīng)值的方法[9]。在工程應(yīng)用當(dāng)中，響應(yīng)面模型一般采用二階多項(xiàng)式來(lái)表達(dá)[10]，其數(shù)學(xué)表達(dá)式如下：

(3)

式中；y(x)—響應(yīng)值；xi—設(shè)計(jì)變量，i=1,2,……,k；β0，βi，βii，βij—待確定的回歸系數(shù)，可以通過(guò)最小二乘法來(lái)估計(jì)多項(xiàng)式的回歸系數(shù)β。

3.1 實(shí)驗(yàn)設(shè)計(jì)(DOE)

樣本點(diǎn)選取對(duì)構(gòu)造的響應(yīng)面模型精確度影響很大，因此需要一個(gè)很好的實(shí)驗(yàn)設(shè)計(jì)來(lái)確保樣本點(diǎn)選取的合理性。本文采用Workbench實(shí)驗(yàn)設(shè)計(jì)中的中心復(fù)合實(shí)驗(yàn)設(shè)計(jì)(CCD)來(lái)選取樣本點(diǎn)，該算法能夠在二階實(shí)驗(yàn)設(shè)計(jì)階段重復(fù)利用一階實(shí)驗(yàn)設(shè)計(jì)，它以最小的實(shí)驗(yàn)循環(huán)提供了諸多設(shè)計(jì)變量與實(shí)驗(yàn)誤差之間的相關(guān)信息，從而以最少的樣本數(shù)量來(lái)建立響應(yīng)面模型[11]，節(jié)約了計(jì)算時(shí)間。CCD實(shí)驗(yàn)設(shè)計(jì)結(jié)果如表3所示。

表3 CCD實(shí)驗(yàn)設(shè)計(jì)結(jié)果

根據(jù)實(shí)驗(yàn)設(shè)計(jì)結(jié)果，通過(guò)最小二乘法計(jì)算出式(3)的回歸系數(shù)矩陣，構(gòu)建各響應(yīng)值的二次多項(xiàng)式響應(yīng)面函數(shù)如下：

G(X)=1.49+0.0214t1+0.000 9h1+0.000 8h2+

(4)

F(X)=0.225 9-0.117 1t1+0.057 6h1-0.085 1h2+

0.811 8t1h1+3.134 4t1h2+0.526 7h1h2-

(5)

3.2 響應(yīng)面模型檢核

(6)

(7)

(8)

(9)

利用式(6,9)對(duì)響應(yīng)面模型進(jìn)行方差分析，計(jì)算結(jié)果如表4所示。

表4 響應(yīng)面模型顯著性校核結(jié)果

3.3 響應(yīng)面分析

經(jīng)響應(yīng)面分析，得出各設(shè)計(jì)變量在變化區(qū)間內(nèi)與目標(biāo)函數(shù)之間的相應(yīng)關(guān)系。由于主要考察的是設(shè)計(jì)變量與門(mén)框的最大等效應(yīng)力之間的關(guān)系，而質(zhì)量只是一個(gè)次要的優(yōu)化目標(biāo)，本文只研究各設(shè)計(jì)變量之間交互作用對(duì)門(mén)框最大等效應(yīng)力的影響。t1、h1交互作用對(duì)等效應(yīng)力的影響如圖5所示。

圖5 t1、h1交互作用對(duì)等效應(yīng)力影響

t1、h2交互作用對(duì)等效應(yīng)力的影響如圖6所示。

圖6 t1、h2交互作用對(duì)等效應(yīng)力影響

h1、h2交互作用對(duì)等效應(yīng)力的影響如圖7所示。

圖7 h1、h2交互作用對(duì)等效應(yīng)力影響

可以看出，3個(gè)設(shè)計(jì)變量之間的交互作用當(dāng)中，筋板厚度t1與筋板向內(nèi)拉伸厚度h2的交互作用對(duì)等效應(yīng)力的影響最大，筋板向內(nèi)拉伸厚度h2與筋板向外拉伸厚度h1的交互作用對(duì)等效應(yīng)力的影響次之，筋板厚度t1與筋板向外拉伸厚度h1的交互作用對(duì)等效應(yīng)力的影響最小。后續(xù)在對(duì)門(mén)框進(jìn)行優(yōu)化時(shí)，可以優(yōu)先針對(duì)筋板厚度和筋板向內(nèi)拉厚度進(jìn)行調(diào)整。分析表明，隨著設(shè)計(jì)變量值的增加，筋板厚度t1、筋板向外拉伸厚度h1與門(mén)框所受最大等效應(yīng)力之間的函數(shù)關(guān)系為單調(diào)遞增的凸拋物線形式，而筋板向內(nèi)拉伸厚度h2與門(mén)框所受最大等效應(yīng)力之間的函數(shù)關(guān)系為凹拋物線形式。t1所對(duì)應(yīng)的等效應(yīng)力變化幅值是最大的，單個(gè)變量與門(mén)框所受最大等效應(yīng)力之間的函數(shù)關(guān)系式會(huì)隨著其他兩個(gè)變量的改變而變化。

4 靈敏度分析

靈敏度分析是指關(guān)注的變量對(duì)整體結(jié)構(gòu)響應(yīng)的影響程度，數(shù)值越大影響越大。這樣可以為后續(xù)提供一個(gè)有效的優(yōu)化思路，當(dāng)要控制結(jié)構(gòu)的響應(yīng)時(shí)，可以改變對(duì)此結(jié)構(gòu)響應(yīng)靈敏度較大的變量來(lái)快速實(shí)現(xiàn)。門(mén)框及其附近區(qū)域最大等效應(yīng)力對(duì)各設(shè)計(jì)變量變化的靈敏度如圖8所示。

圖8 門(mén)框靈敏度分析結(jié)果

由圖8可知，對(duì)門(mén)框及其附近區(qū)域等效應(yīng)力最大值影響最大的是塔架門(mén)框厚度t1，筋板向內(nèi)拉伸厚度h2次之，筋板向外拉伸厚度影響最小。同時(shí)可以發(fā)現(xiàn)門(mén)框的等效應(yīng)力敏感度分析結(jié)果與響應(yīng)面分析的結(jié)果基本一致，這也說(shuō)明了響應(yīng)面模型可靠有效。

5 基于遺傳算法的多目標(biāo)優(yōu)化

對(duì)于多目標(biāo)優(yōu)化問(wèn)題，各優(yōu)化目標(biāo)之間一般不能夠同時(shí)達(dá)到最優(yōu)，因此得出來(lái)的是一系列的有效解，此類解稱為pareto解[14-15]。多目標(biāo)優(yōu)化問(wèn)題就是在一系列的pareto解當(dāng)中找到最有效的幾個(gè)解，之后根據(jù)實(shí)際工程需要，從中再選出最滿意的解。本文基于多目標(biāo)遺傳算法(MOGA)，利用Workbench軟件中目標(biāo)驅(qū)動(dòng)優(yōu)化(GDO)模塊，根據(jù)2.2節(jié)制定的優(yōu)化目標(biāo)函數(shù)及優(yōu)先級(jí)，對(duì)目標(biāo)函數(shù)進(jìn)行優(yōu)化求解。

經(jīng)過(guò)優(yōu)化，得出的3組最優(yōu)pareto解如表5所示。

表5 3組最佳候選設(shè)計(jì)點(diǎn)

由表5可知，雖然候選點(diǎn)1的門(mén)框質(zhì)量與等效應(yīng)力都是最小的，但考慮到制造因素與實(shí)用性等原因，應(yīng)選取候選點(diǎn)2作為最優(yōu)設(shè)計(jì)點(diǎn)。

將候選點(diǎn)2數(shù)據(jù)進(jìn)行圓整并作為最優(yōu)設(shè)計(jì)點(diǎn)，對(duì)該設(shè)計(jì)參數(shù)下的門(mén)框進(jìn)行靜強(qiáng)度計(jì)算，優(yōu)化前后結(jié)果如表6所示。

表6 優(yōu)化前后結(jié)果對(duì)比

由表6可知，優(yōu)化后的門(mén)框質(zhì)量減少73.82%，門(mén)框最大等效應(yīng)力減小26.69%，門(mén)框的極限強(qiáng)度SRFext為1.48，優(yōu)化后門(mén)框性能指標(biāo)有所加強(qiáng)。

6 結(jié)束語(yǔ)

本研究通過(guò)Workbench軟件建立了風(fēng)機(jī)塔筒有限元模型，對(duì)其進(jìn)行靜力分析得到了門(mén)框原始力學(xué)性能；采用中心復(fù)合實(shí)驗(yàn)設(shè)計(jì)方法建立門(mén)框優(yōu)化的響應(yīng)面模型，回歸計(jì)算證明其具有很高的擬合度；通過(guò)響應(yīng)面分析方法研究了門(mén)框尺寸變化對(duì)其所受等效應(yīng)力的影響。

研究結(jié)果表明：隨著設(shè)計(jì)變量值的增加，筋板厚度t1、筋板向外拉伸厚度h1與門(mén)框所受最大等效應(yīng)力之間的函數(shù)關(guān)系為單調(diào)遞增的凸拋物線形式，而筋板向內(nèi)拉伸厚度h2與門(mén)框所受最大等效應(yīng)力之間的函數(shù)關(guān)系為凹拋物線形式，結(jié)果與工程實(shí)際相吻合。

基于響應(yīng)面分析結(jié)果，本文運(yùn)用遺傳算法對(duì)門(mén)框進(jìn)行了多目標(biāo)優(yōu)化。優(yōu)化結(jié)果表明：優(yōu)化后的門(mén)框質(zhì)量減少73.82%，門(mén)框最大等效應(yīng)力減小26.69%，門(mén)框的極限強(qiáng)度SRFext為1.48，表明優(yōu)化是有效的。

參考文獻(xiàn)(References):

[1] 丁小川，李 憶，呂渤林，等.大型風(fēng)電機(jī)組塔架門(mén)洞結(jié)構(gòu)強(qiáng)度分析與優(yōu)化[J].機(jī)械制造，2014，49(11):34-36.

[2] 朱鍇年.風(fēng)力發(fā)電機(jī)組塔架制造工藝與質(zhì)量控制[J].機(jī)電工程技術(shù)，2016，45(10)：136-140.

[3] 任瑞杰，王 林，麻 磊，等.影響風(fēng)機(jī)塔筒門(mén)設(shè)計(jì)的主要因素[J].現(xiàn)代機(jī)械，2011(6):45-47.

[4] 晏紅文，田紅旗，梁裕國(guó)，等.基于Isight的塔架門(mén)框結(jié)構(gòu)的優(yōu)化設(shè)計(jì)[J].機(jī)械強(qiáng)度，2015，37(4):646-650.

[5] Germanischer Lloyd. Rules and guidelines IV Industrial Services Part I: guideline for the certification of wind turbines[M]. Hamburg: Germanischer Lloyd,2010.

[6] 徐芝綸.彈性力學(xué)[M].北京：高等教育出版社，2006.

[7] GB/T 1591-2008.低合金高強(qiáng)度結(jié)構(gòu)鋼[S].北京：中國(guó)標(biāo)準(zhǔn)出版社，2008.

[8] 蔡 新，郭興文，張旭明，等.工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)[M].北京：中國(guó)水利水電出版社，2003.

[9] 胡俊峰，徐貴陽(yáng)，郝亞洲，等.基于響應(yīng)面法的微操作平臺(tái)多目標(biāo)優(yōu)化[J].光學(xué)精密工程，2015,23(4):1096-1103.

[10] 王永菲，王成國(guó).響應(yīng)面法的理論與應(yīng)用[J].中央民族大學(xué)學(xué)報(bào):自然科學(xué)版，2005,14(3):236-239.

[11] 李 莉，張 賽，何 強(qiáng)，等.響應(yīng)面法在試驗(yàn)設(shè)計(jì)與優(yōu)化中的應(yīng)用[J].實(shí)驗(yàn)室研究與探索，2015,34(8):41-45.

[12] 呂 輝，于德介，謝 展，等.基于響應(yīng)面法的汽車盤(pán)式制動(dòng)器穩(wěn)定性優(yōu)化設(shè)計(jì)[J].機(jī)械工程學(xué)報(bào)，2013,49(9):55-30.

[13] 代 翠，孔繁余，董 亮，等.基于響應(yīng)面法的離心泵作透平水力和聲學(xué)性能優(yōu)化[J].農(nóng)業(yè)工程學(xué)報(bào)，2015,13(15):40-46.

[14] 姜 衡，管貽生，邱志成，等.基于響應(yīng)面法的立式加工中心動(dòng)靜態(tài)多目標(biāo)優(yōu)化[J].機(jī)械工程學(xué)報(bào)，2011,47(11):125-132.

[15] 魏心友，曾國(guó)輝，王衛(wèi)軍.基于ABAOUS及modeFRONTIER的多目標(biāo)優(yōu)化方法[J].輕工機(jī)械，2016，34(1)：59-63.