周建平，冷洪坤，溫廣瑞，冉祥鋒，姜 宏

(1.新疆大學(xué) 機(jī)械工程學(xué)院，烏魯木齊 830047；2.西安交通大學(xué) 機(jī)械制造系統(tǒng)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，西安 710049)

多目標(biāo)差分進(jìn)化算法在軸系全息動(dòng)平衡中的應(yīng)用*

周建平1，冷洪坤1，溫廣瑞2，冉祥鋒1，姜 宏1

(1.新疆大學(xué) 機(jī)械工程學(xué)院，烏魯木齊 830047；2.西安交通大學(xué) 機(jī)械制造系統(tǒng)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，西安 710049)

傳統(tǒng)全息動(dòng)平衡方法對(duì)轉(zhuǎn)子進(jìn)行動(dòng)平衡時(shí)，僅以殘余振動(dòng)平方和作為標(biāo)準(zhǔn)對(duì)平衡效果進(jìn)行評(píng)價(jià)，平衡效果受到平衡面和平衡精度及效率的約束。文章利用不同維度的標(biāo)準(zhǔn)benchmark測(cè)試函數(shù)驗(yàn)證了差分進(jìn)化算法的優(yōu)越性，研究了差分進(jìn)化算法配重優(yōu)化的原理，融合模擬退火因子提高算法的尋優(yōu)能力，利用該算法對(duì)軸系進(jìn)行全息動(dòng)平衡優(yōu)化，綜合殘余振動(dòng)平方和、殘余振動(dòng)極值及殘余振動(dòng)均值等多個(gè)目標(biāo)建立評(píng)價(jià)函數(shù)對(duì)軸系動(dòng)平衡狀況進(jìn)行評(píng)價(jià)，通過(guò)實(shí)驗(yàn)驗(yàn)證了該方法的軸系動(dòng)平衡效果良好。

評(píng)價(jià)函數(shù)；差分進(jìn)化；平衡優(yōu)化

0 引言

透平機(jī)械中軸系不平衡使附加載荷增大，是設(shè)備和零部件損壞的最常見的四大故障之一。振動(dòng)作為不平衡的主要表現(xiàn)形式，不僅產(chǎn)生噪音加速軸承磨損、縮短機(jī)械壽命，嚴(yán)重時(shí)會(huì)造成破壞性事故[1]?，F(xiàn)場(chǎng)動(dòng)平衡的常用方法有影響系數(shù)法、振型平衡法和全息動(dòng)平衡法[2-3]。影響系數(shù)法啟停車次數(shù)較多，平衡面選擇不當(dāng)易出現(xiàn)病態(tài)矩陣，導(dǎo)致平衡量不合理，且對(duì)高階振型靈敏度很低。振型平衡法需要預(yù)知轉(zhuǎn)子系統(tǒng)的振型，必須在臨界轉(zhuǎn)速附近進(jìn)行測(cè)量、振型做正交分解或仿真計(jì)算，對(duì)操作者要求較高，不宜實(shí)現(xiàn)計(jì)算機(jī)輔助。全息動(dòng)平衡法利用全息譜技術(shù)集成測(cè)點(diǎn)、截面和整個(gè)軸系的信息，可以在任意轉(zhuǎn)速下將失衡響應(yīng)分解為力和力偶兩部分并分別加以平衡，有效避免了以上兩種傳統(tǒng)平衡方法的不足，但受制于平衡面和平衡精度，使得傳統(tǒng)求解所得配重的平衡效果不能保證最優(yōu)。把轉(zhuǎn)子動(dòng)平衡視為一線性問(wèn)題，將智能優(yōu)化算法與全息動(dòng)平衡相融合，可以有效提高平衡精度。

目前平衡效果評(píng)估多以單目標(biāo)(殘余振動(dòng)平方和最小)為主，難以保證殘余振動(dòng)均勻性[4-5]。本文以殘余振動(dòng)平方和最小、殘余振動(dòng)均方差最小以及殘余振動(dòng)極差值最小三個(gè)方面為約束目標(biāo)嚴(yán)格控制殘余振動(dòng)的大小，基于融合模擬退火因子的差分進(jìn)化算法研究解決全息軸系現(xiàn)場(chǎng)動(dòng)平衡多目標(biāo)優(yōu)化問(wèn)題。

1 優(yōu)化算法選擇

Differential Evolution(DE)算法用于解決連續(xù)變量的全局優(yōu)化問(wèn)題[6]。工作步驟與遺傳算法等其他進(jìn)化算法基本相同，主要包括變異(Mutation)、交叉(Crossover)、選擇(Selection)三種操作。基于種群的全局搜索策略，DE采用實(shí)數(shù)編碼、基于差分的簡(jiǎn)單變異操作和一對(duì)一的競(jìng)爭(zhēng)生存策略，降低了遺傳操作的復(fù)雜性。同時(shí)，DE可以動(dòng)態(tài)跟蹤當(dāng)前的搜索情況，以調(diào)整其搜索策略，具有較強(qiáng)的全局收斂能力和魯棒性。圖1所示為標(biāo)準(zhǔn)差分進(jìn)化算法的流程。

圖1 標(biāo)準(zhǔn)差分進(jìn)化算法流程

一個(gè)良好的算法應(yīng)滿足：初始階段應(yīng)有較強(qiáng)的全局搜索能力，以盡可能發(fā)現(xiàn)較多的全局最優(yōu)解；后期階段應(yīng)具有較強(qiáng)的局部搜索能力，這樣可以提高算法的求解精度和收斂速度。模擬退火算法通過(guò)賦予搜索過(guò)程一種時(shí)變(退火因子)且最終趨于零的概率突跳性，從而可有效避免陷入局部極小并最終趨于全局最優(yōu)的串行結(jié)構(gòu)的優(yōu)化算法。差分進(jìn)化算法基于全局搜索策略，可以融合模擬退火因子來(lái)提高算法的計(jì)算性能[7]。

為驗(yàn)證差分進(jìn)化算法的計(jì)算收斂性能，此處選取幾個(gè)典型的benchmark測(cè)試函數(shù)分別利用不同優(yōu)化算法進(jìn)行求解，根據(jù)求解情況，對(duì)比各函數(shù)的優(yōu)化效果。這幾種算法分別為融合模擬退火因子的差分進(jìn)化算法、帶壓縮因子粒子群算法、綜合學(xué)習(xí)策略粒子群算法、序列前向選擇算法。

(1)單峰函數(shù)

(1)

表1 單峰函數(shù)優(yōu)化對(duì)比

(2)多局部極小值多峰函數(shù)

(2)

表2 多局部極小值多峰函數(shù)優(yōu)化對(duì)比

(3)多局部極小值低維多峰函數(shù)

(3)

表3 多局部極小值低微多峰函數(shù)優(yōu)化對(duì)比

由以上三個(gè)典型的函數(shù)模型求解情況可知，差分進(jìn)化算法適用于不同測(cè)試函數(shù)且求解過(guò)程耗時(shí)最少，相比于其它幾種智能算法具有較好的計(jì)算性能。

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

傳統(tǒng)柔性轉(zhuǎn)子動(dòng)平衡的振動(dòng)信息都是單個(gè)傳感器從轉(zhuǎn)子單向采集，忽略了轉(zhuǎn)子系統(tǒng)各向剛度的差異性，必然給平衡帶來(lái)誤差，降低平衡精度。本文基于將幅、頻、相信息高度綜合的全息軸系動(dòng)平衡方法研究多目標(biāo)平衡優(yōu)化。全息軸系動(dòng)平衡采用三維全息譜表示機(jī)組的振動(dòng)，使軸系振型和加重響應(yīng)直觀化[8]，通過(guò)試重獲得的遷移矩陣具有良好的穩(wěn)定性和可靠性，有利于提高平衡精度。在軸系某個(gè)加重面上添加1000g∠0°配重時(shí)在其它測(cè)重面上引起的振動(dòng)響應(yīng)矩陣即為遷移矩陣。根據(jù)遷移矩陣調(diào)整軸系的配重大小和相位，可顯著減輕軸系不平衡。

設(shè)軸系具有M個(gè)平衡面、N個(gè)測(cè)振面，添加配重后軸系的遷移矩陣為：

(4)

不同配重下的轉(zhuǎn)化算子為：

(5)

添加mg∠α°后機(jī)組的振動(dòng)響應(yīng)為：

(6)

配重后的殘余振動(dòng)方程為：

V0+AW1×C1+AW2×C2+…+AWi×

Ci+…+AWM×CM=V1

(7)

為消除軸系振動(dòng)，為軸系添加一組合適的配重向量pij，平衡后的軸系振動(dòng)方程應(yīng)滿足：

(8)

其中，Ai0為測(cè)振面i的原始振動(dòng)；qi j為配重面j對(duì)測(cè)振面i的影響系數(shù)；pi為平面j上所加配重。一般情況下，軸系中配重面與測(cè)振面的個(gè)數(shù)并不相同，配重面數(shù)小于測(cè)振面數(shù)，現(xiàn)場(chǎng)動(dòng)平衡完全消除各測(cè)振面的振動(dòng)是不現(xiàn)實(shí)的。不同的配重對(duì)應(yīng)不同的遷移矩陣，得到的殘余振動(dòng)也不盡相同，利用智能優(yōu)化算法對(duì)確定的平衡配重進(jìn)行優(yōu)化可以得到最優(yōu)的軸系配重方案。

3 差分進(jìn)化配重優(yōu)化

3.1 平衡優(yōu)化目標(biāo)

傳統(tǒng)軸系現(xiàn)場(chǎng)動(dòng)平衡以殘余振動(dòng)平方和最小為目標(biāo)求得軸系配重方案，不能保證該配重方案下各測(cè)振面的振動(dòng)均得到很好的消除，即有的測(cè)振面可能存在較大的振動(dòng)。針對(duì)這種情況，可以引入轉(zhuǎn)子振動(dòng)的最大值和極差作為評(píng)價(jià)目標(biāo)來(lái)評(píng)判平衡效果，該多目標(biāo)優(yōu)化配重方案的平衡效果可以達(dá)到最優(yōu)[5]。

(9)

該下降半梯形模糊分布函數(shù)中：Mj為Fj(X)的最大值；mj是Fj(X)最小值?？芍?，當(dāng)解愈接近函數(shù)最小值(最優(yōu)解)時(shí)，其隸屬度愈高。

利用線性加權(quán)法并給各目標(biāo)函數(shù)分配不同的權(quán)重系數(shù)，得到各目標(biāo)函數(shù)的平衡優(yōu)化函數(shù)：

H(X)=λ1μ(F1(X))+λ2μ(F2(X))+λ3μ(F3(X))

(10)

其中，λ1、λ2、λ3為線性加權(quán)中各隸屬度函數(shù)不同的權(quán)重系數(shù)，此處仍然選取殘余振動(dòng)平方和最小為主要目標(biāo)評(píng)判平衡標(biāo)準(zhǔn)，殘余振動(dòng)的極值和極差為輔助目標(biāo)，根據(jù)文獻(xiàn)可知在選取λ1=0.7，λ2=0.2，λ3=0.1時(shí)所得配重方案可以取得較優(yōu)的平衡效果。

3.2 差分進(jìn)化配重優(yōu)化

差分進(jìn)化算法的目標(biāo)向量設(shè)為xi=(x0,i,…,xN-1,i)T，試驗(yàn)向量設(shè)為vi=(v0,i,…,vN-1,i)，其中i=0,1,…,NP-1 ，N是目標(biāo)函數(shù)的維數(shù)，一般常用的變異算子有五種，此處列舉兩種：

(11)

其中，G為迭代次數(shù)，r1,r2,r3∈(0,NP-1)為取值范圍內(nèi)各不相等的整數(shù)，F(xiàn)為變異因子。

ui=(u0,i,…,uN-1,i)T

(12)

利用上述公式可以增加參數(shù)向量的多樣性，其中CR為交叉概率。根據(jù)相關(guān)文獻(xiàn)取CR=0.1,F=0.5。

為保證算法前期能具有較高的全局搜索能力以便得到較多的合適解，在后期能提高局部搜索能力加快收斂速度，文章引入模擬退火策略：

(13)

設(shè)軸系有K個(gè)配重面，差分進(jìn)化算法的算法流程如下：

(2)根據(jù)轉(zhuǎn)子振動(dòng)信號(hào)求得轉(zhuǎn)子的原始振動(dòng)及遷移矩陣，計(jì)算轉(zhuǎn)子平衡配重質(zhì)量。根據(jù)所得轉(zhuǎn)子平衡配重及軸系允許的配重閾值確定個(gè)體搜索的區(qū)間。

(3)初始化差分進(jìn)化種群，確定多目標(biāo)下各子目標(biāo)函數(shù)的最大值和最小值。

(4) 根據(jù)原始振動(dòng)、遷移矩陣和個(gè)體所對(duì)應(yīng)的配重計(jì)算出轉(zhuǎn)子殘余振動(dòng)。

(5)根據(jù)多目標(biāo)適應(yīng)度函數(shù)H(X)=λ1μ(F1(X))+λ2μ(F2(X))+λ3μ(F3(X))，計(jì)算出種群個(gè)體的適應(yīng)度。

(6)按照差分策略調(diào)整個(gè)體的交叉和變異，得到新的個(gè)體。

(7)返回第(4)步進(jìn)行迭代，求得最佳配重質(zhì)量與角度。

4 軸系動(dòng)平衡差分進(jìn)化優(yōu)化算法的驗(yàn)證

某電廠2號(hào)汽輪機(jī)發(fā)電機(jī)組進(jìn)行修整，對(duì)轉(zhuǎn)子進(jìn)行動(dòng)平衡，汽輪機(jī)組的結(jié)構(gòu)和傳感器安裝如圖所示。電廠試圖通過(guò)動(dòng)平衡降低高中壓轉(zhuǎn)子的振動(dòng)。

圖2 汽輪機(jī)發(fā)電機(jī)組結(jié)構(gòu)示意圖

平衡過(guò)程中，選擇1#軸承端面法蘭、2#/3#軸承和4#/5#軸承間聯(lián)軸器凸緣作為平衡面，即圖中所示平衡面A、B、C。通過(guò)人工對(duì)機(jī)組進(jìn)行現(xiàn)場(chǎng)動(dòng)平衡，根據(jù)優(yōu)化配重方案，可明顯改善高中壓轉(zhuǎn)子振動(dòng)。

平衡機(jī)組中各配重面的遷移矩陣可以根據(jù)機(jī)組的歷史平衡數(shù)據(jù)計(jì)算得出，也可以利用現(xiàn)場(chǎng)動(dòng)平衡過(guò)程中的配重和軸振信號(hào)計(jì)算得出，利用融入模擬退火因子的差分進(jìn)化算法計(jì)算機(jī)組軸系的動(dòng)平衡配重并對(duì)其進(jìn)行優(yōu)化。表4所示為多目標(biāo)差分進(jìn)化算法優(yōu)化后得到的機(jī)組轉(zhuǎn)子現(xiàn)場(chǎng)動(dòng)平衡配重方案。經(jīng)配重平衡后的機(jī)組轉(zhuǎn)子振動(dòng)及平衡前的振動(dòng)對(duì)比如表5所示，可以看出，經(jīng)多目標(biāo)平衡優(yōu)化后的機(jī)組高中壓缸的振動(dòng)明顯改善，各軸瓦的振動(dòng)也相對(duì)一致，滿足電廠對(duì)機(jī)組軸振控制的要求。

表4 差分進(jìn)化計(jì)算出的平衡配重方案

表5 機(jī)組現(xiàn)場(chǎng)動(dòng)平衡前后及優(yōu)化后軸振對(duì)比

表6 優(yōu)化后的軸振情況

5 結(jié)論

差分進(jìn)化算法與其它智能算法相比，具有計(jì)算量小、收斂速度快的特點(diǎn)。文章將融合模擬退火因子的差分進(jìn)化算法應(yīng)用在石化機(jī)組現(xiàn)場(chǎng)動(dòng)平衡中，以多個(gè)目標(biāo)為約束對(duì)轉(zhuǎn)子動(dòng)平衡效果進(jìn)行評(píng)價(jià)，滿足機(jī)組現(xiàn)場(chǎng)動(dòng)平衡的實(shí)際要求，可以解決軸系動(dòng)平衡的多目標(biāo)優(yōu)化。

[1] 劉春英. 理論淺析回轉(zhuǎn)件轉(zhuǎn)子動(dòng)平衡[J].科技與企業(yè)，2015(22):198-200.

[2] 康成亮. 用于電動(dòng)工具轉(zhuǎn)子的全自動(dòng)平衡機(jī)的影響系數(shù)標(biāo)定方法研究[J].組合機(jī)床與自動(dòng)化技術(shù)，2007(6):23-25.

[3] 劉石,屈梁生. 全息譜技術(shù)在軸系現(xiàn)場(chǎng)動(dòng)平衡方法中的應(yīng)用[J]．熱能動(dòng)力工程，2009,24(1)：24-30.

[4] 瞿紅春,崔秀峰. 基于遺傳算法的動(dòng)平衡優(yōu)化研究[J]．中國(guó)民航大學(xué)學(xué)報(bào)，2013,31(5):24-26.

[5] 張西寧,趙明,溫廣瑞. 一種軸系全息動(dòng)平衡配重的多目標(biāo)優(yōu)化方法[J]．中國(guó)機(jī)械工程，2010，21(16)：1974-1977.

[6] 張春美. 差分進(jìn)化算法理論與應(yīng)用[M].北京：北京理工大學(xué)出版社，2014.

[7] 楊艷霞. 一種基于模擬退火操作的混合差分進(jìn)化算法[J].智能系統(tǒng)學(xué)報(bào)，2014，9(1)：109-114.

[8] 屈梁生. 機(jī)械故障的全息譜診斷原理[M].北京：科學(xué)出版社，2007.

(編輯 李秀敏)

Multi-objective Differential Evolution Algorithm in Holo-balancing

ZHOU Jian-ping1，LENG Hong-kun1，WEN Guang-rui2，RAN Xiang-feng1，JIANG Hong1

(1.College of Mechanical Engineering,Xinjiang University ,Urumqi 830047, China; 2.State Key Laboratory of Mechanical Manufacture System Engineering, Xi′an Jiao Tong University，Xi′an 710049, China)

when the conventional holo-balancing method balances the rotors, only using the sum of residual vibration square as the standard to evaluate the effect of balancing which the effect is limited to the balancing area and accuracy. To solve the problem above, this paper uses different dimensions of standard benchmark test functions to verify the superiority of the differential evolution algorithm，and research on the principle of differential evolution algorithm to optimize counterweight, fusion simulated annealing algorithm factors to increase optimization capacity. Using this algorithm to optimize holographic fields balancing and integrating the maximum ,uniformity and sum of square of residual vibration residual vibration squares and residual vibration to establish evaluation function to evaluate the status of the fields balancing with multiple targets. The effectiveness of the proposed method is verified by a practical balancing application.

evaluation function;differential evolution;balancing optimization

1001-2265(2017)05-0017-04

10.13462/j.cnki.mmtamt.2017.05.005

2016-07-28；

2016-09-06

新疆維吾爾自治區(qū)科技支疆項(xiàng)目(項(xiàng)目201404071245)

周建平(1974—)，男，江蘇贛榆人，新疆大學(xué)教授，博士，研究方向?yàn)樘胤N加工與數(shù)控技術(shù)、信號(hào)采集與分析，(E-mail)rxf_9106@163.com。

TH124；TG659

A