摘 要為實(shí)現(xiàn)壓氣機(jī)葉片的排序優(yōu)化設(shè)計(jì)，根據(jù)葉片重量和頻率的工藝要求，提出了逐步調(diào)整優(yōu)化算法，實(shí)質(zhì)是先考慮一個(gè)限制條件，再利用另一個(gè)約束條件進(jìn)行逐步調(diào)優(yōu)。建立了兩個(gè)優(yōu)化模型并利用MATLAB來(lái)實(shí)現(xiàn)，得到了一組最優(yōu)結(jié)果。以一實(shí)例驗(yàn)證了算法的時(shí)間復(fù)雜度。結(jié)果表明，逐步調(diào)整優(yōu)化算法運(yùn)算簡(jiǎn)單、容易理解，為壓氣機(jī)葉片排序提供了一種新的方法，在壓氣機(jī)葉片排序上有很好的應(yīng)用前景。

關(guān)鍵詞 逐步調(diào)整優(yōu)化法;最優(yōu)解;順、逆序插值法;算法的時(shí)間復(fù)雜度

中圖分類(lèi)號(hào) TK473文獻(xiàn)標(biāo)識(shí)碼A文章編號(hào)1673-9671-(2009)111-0074-02

葉片是發(fā)動(dòng)機(jī)及壓縮機(jī)上的重要部件，由于葉片都是經(jīng)手工打磨加工的。因而各葉片之間難免有質(zhì)量和振蕩頻率的差別，安裝時(shí)需對(duì)各葉片進(jìn)行排序，其排序過(guò)程是很難通過(guò)手工完成的。而排序不當(dāng)，將會(huì)對(duì)發(fā)動(dòng)機(jī)及壓縮機(jī)的安全造成很大的隱患。國(guó)內(nèi)有關(guān)壓氣機(jī)葉片排序的文獻(xiàn)很少，目前較好的方法是:先測(cè)出葉輪上的每張葉片繞葉輪中心轉(zhuǎn)動(dòng)時(shí)的力矩以及葉輪不平衡力矩，然后根據(jù)這些量用計(jì)算機(jī)算出葉片在葉輪上的排序。而算葉片排序的方法，一般考慮最優(yōu)化算法，但因?yàn)槿~片數(shù)量多，故最優(yōu)化算法中自變量太多，以致最優(yōu)化算法(如分枝定界法[1])難以應(yīng)用。

本文針對(duì)現(xiàn)有排序方法過(guò)于復(fù)雜、難以應(yīng)用的缺點(diǎn)提出了逐步調(diào)整優(yōu)化算法。并以具體的工程中的壓氣機(jī)葉片排序問(wèn)題為例詳細(xì)介紹了該算法。

1 逐步調(diào)整優(yōu)化法

壓氣機(jī)葉片因加工出的重量和頻率不同，安裝時(shí)需要按工藝要求重新排序。某壓氣機(jī)葉片數(shù)值如表1所示，重量單位g ，頻率單位Hz。壓氣機(jī)24片葉片均勻分布在一圓盤(pán)邊上，分成六個(gè)象限，每象限4片葉片的總重量與相鄰象限4片葉片的總重量之差不允許超過(guò)一定值，設(shè)定為8g;兩相鄰葉片頻率差盡量大，使相鄰葉片頻率差不小于一定值，設(shè)定為6Hz。當(dāng)葉片確實(shí)不滿(mǎn)足上述要求時(shí)，允許更換少量葉片。按照逐步調(diào)整優(yōu)化算法，能快速的找出符合上述條件的最優(yōu)解，即葉片的最佳排列順序，并按照排列順序輸出葉片號(hào)。

1.1 理論基礎(chǔ)

逐步調(diào)整優(yōu)化法即采用先考慮一個(gè)限制條件，再利用另一個(gè)約束條件進(jìn)行逐步調(diào)優(yōu)的方法。在該算法中認(rèn)為葉片的重量和頻率間存在相關(guān)關(guān)系，認(rèn)為重量較大的頻率也較大，它們的關(guān)系如圖1所示。

重量和頻率之間的相關(guān)關(guān)系和圖1是可以互相證明的，所以為了同時(shí)滿(mǎn)足相鄰葉片頻率之差大于一定值，各象限間重量之差小于一定值，只要使得頻率或者重量是按照大小大小順序排列的，另一種也必定是滿(mǎn)足這種關(guān)系。所以無(wú)論先考慮重量還是頻率都應(yīng)該盡量遵循“大小大小”排列順序。

1.2 頻率約束模型:

先只考慮頻率的限制條件，即相鄰葉片間的頻率差要大于一定值，并使頻率差盡可能大。在此約束條件下，從最優(yōu)解(即相鄰葉片間差值盡可能地大)開(kāi)始，求出單約束條件下的最優(yōu)解，對(duì)求出的最優(yōu)解檢驗(yàn)重量的分布是否符合條件。這樣可以大大減少尋找符合條件及最優(yōu)解排序的求解循環(huán)次數(shù)，并且可以同時(shí)尋找可能解和最優(yōu)解。

1.2.1頻率約束條件下最佳排列順序的求解

s.t

1.2.2逆序插值法

在頻率約束條件下葉片的排列順序:

首先與順序插值法相同將24個(gè)數(shù)據(jù)分成兩組，不同的是，把第二組數(shù)據(jù)反向插入到第一組數(shù)據(jù)中，這時(shí)排列在前面的數(shù)據(jù)之間的差值最大，隨著序列數(shù)的減少，差值隨之減少。把不符合的數(shù)據(jù)插入到差值最大的一組數(shù)據(jù)中。即“最大最小”組合，如此循環(huán)，直到每相鄰的兩組數(shù)據(jù)之間的差值均大于一定值(本題中設(shè)定定值為6Hz)。由于頻率相同的葉片重量有可能不同，故在只考慮頻率約束條件下先將頻率相同的葉片按其重量由大到小進(jìn)行排列，示意圖如圖2。

當(dāng)該組排序不滿(mǎn)足重量要求時(shí)，繼續(xù)采用上面所述的“最大最小”差值辦法尋找合適的解，如圖3所示。

對(duì)逆序插值法得出的頻率約束下得到排序進(jìn)行檢驗(yàn)，由MATLAB程序得到的一組最優(yōu)解:

1.3重量約束模型

先只需要考慮各相鄰象限內(nèi)葉片的重量之差不超過(guò)定值(設(shè)定為8g)，而不考慮其他約束條件，及各象限內(nèi)葉片的排列順序。按照重量約束條件把圓盤(pán)分為6個(gè)象限后，再對(duì)所有符合條件的排列順序按照頻率約束進(jìn)行檢驗(yàn)，這樣可以大大減少循環(huán)次數(shù)和搜尋時(shí)間。

s.t (2)，式中，為各象限葉片的

總質(zhì)量

1.3.1重量約束排列順序

(Ⅰ)交替排序法:

把重量按從大到小的順序排列，每一葉片的重量離平均值 m的距離為m- m?？芍?，越在序列兩端，距離越大。左邊的是正向最大，右邊是負(fù)向最大。正向偏移量與對(duì)應(yīng)位置負(fù)向偏移量組合可以使其綜合距離較近。并避免大大小小相遇。依此取出六組數(shù)據(jù)的重量之和相差較小。由MATLAB程序得到的一組最優(yōu)解，如表3所示，表中各象限內(nèi)葉片沒(méi)有排序。

表3 交替排序法給出的葉片排列順序

第一象限第二象限第三象限第四象限第五象限第六象限

2 10 4 813 11 16 91 20 18717 12 5 2214 23 3 246 21 15 19

2 結(jié)果與討論

2.1頻率約束模型結(jié)果討論

對(duì)得出的各象限重量和頻率的數(shù)據(jù)進(jìn)行方差和均值分析:

均值:(3)，方差: (4)

在對(duì)重量進(jìn)行分析時(shí)，xi表示各象限重量之和，i=1，……，6，方差越小，表示重量之差越小，得出的數(shù)據(jù)越符合期望。在對(duì)頻率進(jìn)行分析時(shí)，xi表示每葉片的頻率，i=1，……，24，方差越大，表示頻率之差越大，得出的數(shù)據(jù)越符合期望。兩種方法所得的重量和頻率期望和方差結(jié)果如圖4所示。

從圖4可以看出，頻率約束模型下的逆序插值法可以達(dá)到葉片的工藝性要求，各象限間重量差較小、頻率差較大，且經(jīng)比較得到逆序插值法得出的解更優(yōu)

2.2重量約束模型結(jié)果討論

表5表示了兩種排序法得到的結(jié)果。

從圖5中我們可以看到重量約束模型下的交替排序法可以達(dá)到葉片的工藝性要求，各象限間重量差較小、頻率差較大，且經(jīng)比較得到交替排序法得出的解更優(yōu)。

3 結(jié)語(yǔ)

本文引用了兩組工程中的實(shí)際壓氣機(jī)葉片數(shù)據(jù)來(lái)解釋逐步調(diào)整優(yōu)化算法，并利用其中一組數(shù)據(jù)對(duì)該算法進(jìn)行了實(shí)例驗(yàn)證，結(jié)果表明逐步調(diào)整優(yōu)化算法求解簡(jiǎn)便，結(jié)果優(yōu)化，具有極大的推廣價(jià)值，值得在工程中廣泛實(shí)用。但本文所探討的數(shù)據(jù)量有限，該算法在一些葉片數(shù)量較多(如在100片以上)的機(jī)器上是否可以應(yīng)用的問(wèn)題有待進(jìn)一步探討。

參考文獻(xiàn)

[1] 傅國(guó)耀. 透平機(jī)械葉片排序的一個(gè)實(shí)用算法[J]. 應(yīng)用科學(xué)學(xué)報(bào)， 1988 6(4):335-339.

[2] 西北工業(yè)大學(xué)2003年數(shù)學(xué)建模B題，accessed 28th July，2008.

[3] 岳健民，等. 葉片排序的逐步調(diào)優(yōu)模擬搜索算法[J]. 四川理工學(xué)院學(xué)報(bào)(自然科學(xué)版). 第18卷第3期，2005.3 .

[4] 吳建國(guó)，等編. 數(shù)學(xué)建模案例精編[M]. 北京:中國(guó)水利水電出版社，2005.

[5] 王沫然. Matlab6.0與科學(xué)計(jì)算[M]. 北京:電子工業(yè)出版社，2001.

[6] 唐紹軍，等. 遺傳算法對(duì)壓氣機(jī)葉片排序的應(yīng)用[J]. 航空動(dòng)力學(xué)報(bào)， 2005 20(3):518-522.

[7] 李平等. 汽輪機(jī)動(dòng)葉片安裝中的力矩平衡及計(jì)算機(jī)優(yōu)化排序[J]. 西北電力技術(shù)，1999(2).