張西寧 趙 明 溫廣瑞

西安交通大學(xué)機(jī)械制造系統(tǒng)工程國家重點(diǎn)實(shí)驗(yàn)室，西安，710049

0 引言

回轉(zhuǎn)機(jī)械軸系不平衡會(huì)引起轉(zhuǎn)子撓曲和內(nèi)應(yīng)力，加速軸承和軸封等零件的磨損，降低機(jī)器的工作效率，嚴(yán)重時(shí)會(huì)引起各種事故［1］。對(duì)于軸系失衡常采用影響系數(shù)法［2］、模態(tài)平衡法［3］、全息動(dòng)平衡法［1］等方法進(jìn)行現(xiàn)場校正。由于軸系動(dòng)平衡本身與單個(gè)轉(zhuǎn)子動(dòng)平衡存在差異，同時(shí)平衡還受到測量面選擇、軸系各支承狀況等因素的影響，故直接求解的平衡配重往往不能保證達(dá)到最優(yōu)。本文基于粒子群優(yōu)化算法研究解決全息軸系現(xiàn)場動(dòng)平衡的多目標(biāo)優(yōu)化問題。

1 全息軸系動(dòng)平衡

為了消除傳統(tǒng)動(dòng)平衡方法由于各向支承剛度相等假設(shè)所帶來的計(jì)算誤差［4］，本文以全息軸系動(dòng)平衡方法為基礎(chǔ)展開平衡優(yōu)化。全息軸系動(dòng)平衡計(jì)算需要用到轉(zhuǎn)頻橢圓、全息初相點(diǎn)和遷移矩陣等概念。轉(zhuǎn)頻橢圓是轉(zhuǎn)子各測量面上兩個(gè)相互垂直振動(dòng)信號(hào)一倍頻分量合成的橢圓［5］。當(dāng)轉(zhuǎn)子上的鍵相槽對(duì)準(zhǔn)鍵相傳感器時(shí)轉(zhuǎn)子中心在轉(zhuǎn)頻橢圓上的位置為全息初相點(diǎn)。轉(zhuǎn)頻橢圓中心到全息初相點(diǎn)的向量為全息初相矢。遷移矩陣是在某加重面上添加1000g∠0°試重在各個(gè)測量面上引起的振動(dòng)響應(yīng)矩陣。獲得機(jī)組遷移矩陣后，通過調(diào)整各加重面上試重的大小和角度，可減小或消除機(jī)組各測量面上的原始振動(dòng)。

設(shè)軸系具有M個(gè)平衡面、N個(gè)測振截面，平衡方案中只有一個(gè)平衡轉(zhuǎn)速，那么我們期望尋找一組合適的配重向量（pij），使得配重后轉(zhuǎn)子系統(tǒng)滿足以下的振動(dòng)方程：

式中，Ai0為i測振截面的原始振動(dòng)；cij為遷移矩陣中j平衡面對(duì)i測振截面的影響系數(shù)；pj為在平衡面j上施加的配重。

傳統(tǒng)的影響系數(shù)平衡法根據(jù)M與N的不同，通過求解方程得到配重的解，或按照Goodman提出的最小二乘方法［3］求解殘余振動(dòng)平方和最小意義下的配重解。全息軸系動(dòng)平衡按照軸系失衡響應(yīng)的對(duì)偶性原則來初步確定平衡配重。

一般平衡的加重面數(shù)M都小于測振截面數(shù)N，因此，現(xiàn)場平衡不可能將各測量面的振動(dòng)都完全抵消。不同的加重方案會(huì)得到不同的殘余振動(dòng)，有必要對(duì)確定的平衡配重作進(jìn)一步的優(yōu)化。在初步確定的平衡配重的基礎(chǔ)上進(jìn)行優(yōu)化，進(jìn)而可得到最優(yōu)平衡配重。

2 軸系平衡配重優(yōu)化

2.1 平衡優(yōu)化目標(biāo)選擇

以往動(dòng)平衡以殘余振動(dòng)平方和最小來求解平衡配重。平衡結(jié)果雖然殘余振動(dòng)平方和最小，但并不能保證每個(gè)測點(diǎn)的殘余振動(dòng)都比較小，可能存在振動(dòng)比較大的測點(diǎn)。為此我們選用最小化的殘余振動(dòng)平方和殘余振動(dòng)最大值minF2（X）＝min（maxxi）及各測點(diǎn)殘余振動(dòng)差異性minF3（X）＝min（maxxi－minxi）來衡量機(jī)組殘余振動(dòng)。這樣最優(yōu)配重方案的求解就成為一個(gè)多目標(biāo)優(yōu)化問題。

對(duì)于多目標(biāo)規(guī)劃問題可根據(jù)各個(gè)目標(biāo)的重要程度設(shè)置不同的權(quán)重，通過線性加權(quán)求和法將多目標(biāo)問題轉(zhuǎn)化為單目標(biāo)問題求解。然而由于各目標(biāo)函數(shù)項(xiàng)大小及量綱的差異，故容易造成一些目標(biāo)函數(shù)總是處于支配地位，而其他的目標(biāo)函數(shù)得不到有效的優(yōu)化［6］。為此，基于模糊數(shù)學(xué)原理建立了各目標(biāo)函數(shù)的隸屬函數(shù)。各目標(biāo)函數(shù)的隸屬函數(shù)構(gòu)造為下降半梯形模糊分布函數(shù)：

其中，Mj為目標(biāo)函數(shù)Fj（X）在可行域內(nèi)的最大值；mj為可行域內(nèi)的最小值。由式（1）可知：某個(gè)解越接近目標(biāo)函數(shù)的最優(yōu)解（最小值），其隸屬度越高。

采用線性加權(quán)和法對(duì)各目標(biāo)的隸屬函數(shù)進(jìn)行加權(quán)求和，得到總的平衡優(yōu)化目標(biāo)函數(shù)：

其中，λ1、λ2、λ3分別為各隸屬度函數(shù)的權(quán)重系數(shù)，其取值取決于各目標(biāo)函數(shù)的重要程度，在這里考慮到優(yōu)化目標(biāo)F1（X）的主導(dǎo)作用，將權(quán)重系數(shù)定為λ1＝0.7，λ2＝0.2，λ3＝0.1，這能使H（X）取得較大值的配重方案是較優(yōu)的。

2.2 粒子群配重優(yōu)化算法

粒子群優(yōu)化算法［7-8］是基于對(duì)鳥群捕食行為模擬的演化算法。粒子群優(yōu)化求解時(shí)，問題解對(duì)應(yīng)于搜索空間中一只鳥的位置，稱這些鳥為“粒子”。每個(gè)粒子都有自己的位置、速度和一個(gè)由被優(yōu)化函數(shù)決定的適應(yīng)值。各個(gè)粒子記憶、追隨當(dāng)前的最優(yōu)粒子在解空間中搜索。每次迭代過程不是完全隨機(jī)的，尋優(yōu)中將會(huì)以當(dāng)前找到的較好解為依據(jù)來尋找最優(yōu)解。每次迭代中粒子按照公式，通過跟蹤兩個(gè)“極值”來更新自己的速度和位置：一個(gè)是粒子自身找到的最好解，叫做個(gè)體極值點(diǎn)；另一個(gè)是整個(gè)種群目前找到的最好解，稱為全局極值點(diǎn)。跟蹤公式如下：

其中，Xi＝ （xi1，xi2，…，xiD）和Vi＝ （vi1，vi2，…，viD）分別是粒子i的位置和速度是粒子i在第k次迭代中第d維的速度；c1、c2是加速系數(shù)，分別調(diào)節(jié)向全局最好粒子和個(gè)體最好粒子方向飛行的最大步長，若太小則粒子可能遠(yuǎn)離目標(biāo)區(qū)域，若太大則會(huì)導(dǎo)致突然向目標(biāo)區(qū)域飛去，或飛過目標(biāo)區(qū)域［7］，本文參考相關(guān)文獻(xiàn)取c1＝c2＝2.05；rand1、rand2是［0，1］之間的隨機(jī)數(shù)是粒子i在第k次迭代中第d維的當(dāng)前位置；pbestid是粒子i在第d維的個(gè)體極值點(diǎn)的位置；gbestd是整個(gè)群在第d維的全局極值點(diǎn)的位置。為防止粒子遠(yuǎn)離搜索空間，粒子的每一維速度vd都會(huì)被鉗位在［－vdmax，vdmax］之間，vdmax太大，粒子將飛離最好解，太小將會(huì)陷入局部最優(yōu)［7］，本文將vdmax設(shè)計(jì)為0.1xdmax。采用將慣性因子從1到0.2隨迭代過程線性減小的方法，使其在前期具有較高的搜索能力以得到合適的種子，在后期具有較高的收斂能力以加快收斂速度。

當(dāng)加重面數(shù)為K時(shí)，粒子群優(yōu)化算法流程如下：

（2）根據(jù)遷移矩陣和原始振動(dòng)計(jì)算配重。再根據(jù)計(jì)算配重和允許配重范圍確定粒子搜索空間。

（3）對(duì)粒子群體進(jìn)行初始化，并確定各個(gè)子目標(biāo)函數(shù)的最大值Mj和最小值mj。

（4）由原始振動(dòng)、遷移矩陣和各粒子所對(duì)應(yīng)的配重計(jì)算出殘余振動(dòng)。

（5）根據(jù)適應(yīng)度函數(shù)H（X）＝λ1μ（F1（X））＋λ2μ（F2（X））＋λ3μ（F3（X）），計(jì)算種群中各粒子的適應(yīng)度。

（6）按照規(guī)則調(diào)整各粒子的位置和速度，產(chǎn)生新一代群體。

（7）返回第（4）步進(jìn)行迭代，最終求出模糊最優(yōu)解。

3 軸系動(dòng)平衡粒子群優(yōu)化方法的驗(yàn)證

2001年10 月某電廠2號(hào)汽輪發(fā)電機(jī)組振動(dòng)實(shí)施動(dòng)平衡，機(jī)組結(jié)構(gòu)和傳感器安裝示意見圖1。平衡目標(biāo)是降低高中壓轉(zhuǎn)子的振動(dòng)。使用的平衡面為1號(hào)軸承端面法蘭、2號(hào)3號(hào)軸承和4號(hào)5號(hào)軸承間聯(lián)軸器凸緣，依次記為A、B、C面。對(duì)機(jī)組進(jìn)行人工現(xiàn)場動(dòng)平衡，經(jīng)過多次加重，使得高中壓轉(zhuǎn)子的振動(dòng)明顯減小，遠(yuǎn)離報(bào)警水平。

圖1 300MW機(jī)組結(jié)構(gòu)示意圖

這里根據(jù)該機(jī)組當(dāng)時(shí)動(dòng)平衡過程中的加重?cái)?shù)據(jù)和振動(dòng)信號(hào)，計(jì)算出了機(jī)組各加重面的遷移矩陣，并采用粒子群方法進(jìn)行平衡配重計(jì)算和優(yōu)化，表1是經(jīng)多目標(biāo)優(yōu)化得到的平衡配重方案。機(jī)組現(xiàn)場平衡前后的振動(dòng)及經(jīng)多目標(biāo)優(yōu)化后的振動(dòng)如表2和圖2所示。經(jīng)多目標(biāo)優(yōu)化后機(jī)組的高中壓缸的振動(dòng)顯著降低，同時(shí)各個(gè)軸瓦的振動(dòng)也比較均勻一致，符合人們對(duì)平衡結(jié)果的期望。因此，使用多目標(biāo)優(yōu)化，可以得到平衡效果更好、更滿足人們期望的配重方案。

表1 多目標(biāo)優(yōu)化后的平衡配重方案

表2 機(jī)組現(xiàn)場平衡前后及多目標(biāo)優(yōu)化后計(jì)算的振動(dòng)

為了進(jìn)一步驗(yàn)證粒子群優(yōu)化方法的性能，這里又分別采用遺傳算法和粒子群優(yōu)化算法對(duì)該機(jī)組的平衡方案進(jìn)行了優(yōu)化。遺傳算法優(yōu)化時(shí)的雜交率、變異率根據(jù)文獻(xiàn)［9］建議取為0.4和0.01。兩種方法的種群規(guī)模均設(shè)定為150。為了消除隨機(jī)性影響，結(jié)果是30次實(shí)驗(yàn)數(shù)據(jù)的平均值。圖3是兩種方法優(yōu)化過程的收斂曲線。從圖3中可以看出，在種群規(guī)模相同時(shí)粒子群優(yōu)化算法的收斂速度明顯高于遺傳算法的收斂速度。

圖2 平衡前后機(jī)組的振動(dòng)三維全息譜

圖3 粒子群優(yōu)化算法與遺傳算法的收斂過程比較

4 結(jié)論

（1）本文提出了評(píng)價(jià)軸系平衡效果的三個(gè)子目標(biāo)：機(jī)組殘余振動(dòng)平方和、殘余振動(dòng)最大值、殘余振動(dòng)差異性，滿足和符合現(xiàn)場機(jī)組平衡的實(shí)際要求。

（2）基于模糊數(shù)學(xué)方法構(gòu)造的平衡優(yōu)化總體目標(biāo)可實(shí)現(xiàn)機(jī)組子目標(biāo)間的優(yōu)化協(xié)調(diào)。

（3）實(shí)際平衡優(yōu)化案例驗(yàn)證結(jié)果表明本文研究的全息軸系動(dòng)平衡粒子群優(yōu)化方法可實(shí)現(xiàn)軸系動(dòng)平衡多目標(biāo)優(yōu)化；與遺傳算法比較，粒子群優(yōu)化算法具有計(jì)算量小、收斂速度快的特點(diǎn)。

［1］ 吳松濤.大型機(jī)組現(xiàn)場動(dòng)平衡原理與技術(shù)［D］.西安：西安交通大學(xué)，2003.

［2］ Grobel L P.Balancing Turbines－generator Rotors［J］.General Electric Review，1953，56（4）：22.

［3］ Goodman T P.A Least－squares Method for Computing Balance Corrections［J］.ASME Transactions，Journal of Engineering for Industry，1964，86（3）：273－279.

［4］ 屈梁生，邱海，徐光華.全息動(dòng)平衡技術(shù)：原理與實(shí)踐［J］.中國機(jī)械工程，1998，9（1）：60－63.

［5］ 屈梁生，王琇峰.全息動(dòng)平衡原理［J］.陜西電力，2007（1）：1－5.

［6］ 賈兆紅，陳華平，孫耀暉.多目標(biāo)粒子群優(yōu)化算法在柔性車間調(diào)度中的應(yīng)用［J］.小型微型計(jì)算機(jī)系統(tǒng)，2008，29（5）：885－889.

［7］ Eberhart R C，Shi Y.Particle Swarm Optimization：Developments，Applications and Resources［C］／／Proceedings of the IEEE Congress on Evolutionary Computation.Piscataway，NJ，2001：81－86.

［8］ 楊維，李歧強(qiáng).粒子群優(yōu)化算法綜述［J］.中國工程科學(xué)，2004，16（15）：87－94.

［9］ 黎鈞琪，石國楨.遺傳算法交叉率與變異率關(guān)系的研究［J］.武漢理工大學(xué)學(xué)報(bào)，2003，27（1）：97－99.