李波

摘 要：本文通過極大似然法、雙線性回歸法、相關(guān)系數(shù)法及概率權(quán)重矩法的威布爾估計算法對風(fēng)機載荷相關(guān)變量（葉根面內(nèi)彎矩、葉根面外彎矩及葉尖撓度）進行載荷外推，并通過極大似然值對該幾種方法進行比較，提出適合風(fēng)機載荷變量的估計算法，為使用三參數(shù)威布爾分布對風(fēng)機載荷外推提供了一定的參考。

關(guān)鍵詞：威布爾分布；風(fēng)機載荷；極大似然；外推

中圖分類號：O213 文獻標(biāo)志碼：A

風(fēng)機載荷計算結(jié)果通常有兩種表達形式：極限載荷和疲勞載荷。極限載荷可以通過風(fēng)機在其整個服役期間所經(jīng)歷的最大載荷再乘以相應(yīng)的安全系數(shù)得到。當(dāng)來流風(fēng)是確定性風(fēng)時，產(chǎn)生的載荷為確定性載荷，其極值就是提取所有工況中的載荷最大值；當(dāng)來流風(fēng)是隨機性風(fēng)時（湍流模型NTM），因其特有的隨機性和不穩(wěn)定性，有必要在統(tǒng)計學(xué)的基礎(chǔ)上分析載荷的極值，以便得到一個適當(dāng)?shù)奶匦暂d荷。按照IEC61400-1：2005第三版的要求，我們需要得到50年一遇的載荷極值，而計算機仿真的有限次數(shù)很難直接獲得該條件下的理想數(shù)值，為了得到合理的極限載荷值，我們需要借助概率的方法進行載荷外推從而得到該值。

隨機變量分布中的威布爾分布，在可靠性工程領(lǐng)域中被廣泛應(yīng)用，由于它可以利用概率值很容易地推斷出它的分布參數(shù)，因而被廣泛應(yīng)用于各種壽命試驗的數(shù)據(jù)處理。三參數(shù)威布爾分布應(yīng)用于材料的壽命試驗中，通常情況下，它具有比對數(shù)正態(tài)分布更大的適用性。目前，對三參數(shù)威布爾分布的形狀參數(shù)、尺度參數(shù)及位置參數(shù)的估計方法有多種。在風(fēng)機載荷外推中，根據(jù)IEC的標(biāo)準(zhǔn)，我們往往需要得到的參數(shù)是葉根面內(nèi)彎矩Mx，面外彎矩My及葉尖撓度。本文將三參數(shù)威布分布應(yīng)用于風(fēng)機載荷外推中，并采用極大似然法、概率權(quán)重矩法、雙線性回歸法及相關(guān)系數(shù)法分別進行參數(shù)估計，為了對各個估計算法進行比較，我們通過對比最大似然函數(shù)值從而選出適合每個變量的參數(shù)估計法，即找出適合于葉根面內(nèi)彎矩Mx，面外彎矩My及葉尖撓度的參數(shù)估計法，為風(fēng)機載荷外推提供一定的參考。

1 三參數(shù)威布分布模型

三參數(shù)威布爾分布函數(shù)是：

（1）

其中γ、η、m分別是位置參數(shù)、尺度參數(shù)和形狀參數(shù)。

2 計算分析流程

計算分析流程如圖1所示，首先根據(jù)工況要求并基于Bladed平臺，將提取所有仿真工況的極值即載荷最大值，將這些數(shù)值根據(jù)葉根面內(nèi)彎矩、葉根面外彎矩及葉尖撓度進行分類，再運用不同的估計算法如極大似然法等進行參數(shù)估計，需要說明的是，本文中所用到的數(shù)據(jù)與相關(guān)系數(shù)法是引用文獻，然后求出各個參數(shù)估計法對應(yīng)的參數(shù)即形狀參數(shù)、尺度參數(shù)和位置參數(shù)。根據(jù)IEC標(biāo)準(zhǔn)，50年一遇的超越概率是3.8×10-7，求此概率下的載荷極值，同時，為了進行結(jié)果比較，本文也給出了各個算法的極大似然值。

Excel中的規(guī)劃求解是一組命令的組成部分，這些命令有時也稱作假設(shè)分析工具。 使用規(guī)劃求解，可求得工作表上某個單元格中公式的最優(yōu)值。規(guī)劃求解將對直接或間接與目標(biāo)單元格中公式相關(guān)聯(lián)的一組單元格中的數(shù)值進行調(diào)整，最終在目標(biāo)單元格公式中求得期望的結(jié)果。在創(chuàng)建模型過程中，可以對規(guī)劃求解模型中的可變單元格數(shù)值應(yīng)用約束條件，而且約束條件可以引用其他影響目標(biāo)單元格公式的單元格。本文中，極大似然法和概率權(quán)重矩法對Excel中的規(guī)劃求解功能得到了很好的應(yīng)用，同時利用規(guī)劃求解，已知超越概率而反求極值變得非?？旖莘奖?。

3 結(jié)果匯總

各個參數(shù)估計法匯總見表1。

結(jié)語

在本算例風(fēng)機載荷外推中，對于葉根面內(nèi)彎矩、面外彎矩，推薦使用極大似然法，而葉尖撓度推薦使用雙線性回歸法。需要說明的是，隨著科技的進步，載荷外推的算法會越來越多，本文給出的算例所使用的方法畢竟是有限的，需要在實際應(yīng)用中不斷地驗證，最終的目的是為風(fēng)機載荷外推提供一定的參考。

參考文獻

[1]史景釗，任學(xué)軍，陳新昌，等.一種三參數(shù)Weibull分布極大似然估計的求解方法[J].河南科學(xué)，2009，27（7）：832-834.

[2]莊渭峰.用微機實現(xiàn)威布爾參數(shù)的雙線性回歸最小二乘法估計[J].電子產(chǎn)品可靠性與環(huán)境試驗，1999（5）：2-7.

[3]李成本，高德忠.外推法在風(fēng)力發(fā)電機載荷計算中的應(yīng)用[J].一重技術(shù)，2015（1）：9-13.