基于遺傳算法的高速軋輥磨床磨頭液體動靜壓軸承的優(yōu)化設(shè)計

基于遺傳算法的高速軋輥磨床磨頭液體動靜壓軸承的優(yōu)化設(shè)計

吳懷超1,2令狐克均1孫官朝1李哲1張順風1

1.貴州大學,貴陽,5500252.清華大學摩擦學國家重點實驗室，北京,100084

摘要：液體動靜壓軸承性能的好壞直接影響著整臺磨床的磨削性能。針對一種可調(diào)節(jié)式的高速軋輥磨床磨頭液體動靜壓軸承，分析了其結(jié)構(gòu)特點和工作原理，并對其結(jié)構(gòu)尺寸和工作參數(shù)進行了設(shè)計計算。在此基礎(chǔ)上，以單位承載量下的總功率損失最小作為優(yōu)化設(shè)計的目標函數(shù)，運用遺傳算法對該軸承進行了優(yōu)化設(shè)計。優(yōu)化結(jié)果表明：優(yōu)化后目標函數(shù)值從0.0808下降到0.0678，軸承的總功率損失由初始設(shè)計的766.9W降低為642.6W，功耗降低約16%，而且軸承的結(jié)構(gòu)尺寸得到減小，其剛度亦得到明顯的提高。由此可見，通過優(yōu)化設(shè)計，不僅降低了軸承的制造成本，而且明顯改善了軸承的綜合性能。

關(guān)鍵詞：高速軋輥磨床；液體動靜壓軸承；優(yōu)化設(shè)計；遺傳算法

中圖分類號：TH133.3

收稿日期：2015-03-17

基金項目：國家自然科學基金資助項目(51165002,51465008);清華大學摩擦學國家重點實驗室開放基金資助項目(SKLTKF12A11);貴州省優(yōu)秀青年科技人才培養(yǎng)對象專項資金項目(黔科合人字(2013)17號)

作者簡介：吳懷超，男,1975年生。貴州大學機械工程學院教授,清華大學摩擦學國家重點實驗室博士后研究人員。主要研究方向為高效磨削及裝備。發(fā)表論文30余篇。令狐克均,男,1990年生。貴州大學機械工程學院碩士研究生。孫官朝,男,1989年生。貴州大學機械工程學院碩士研究生。李哲,男,1992年生。貴州大學機械工程學院碩士研究生。張順風,男,1988年生。貴州大學機械工程學院碩士研究生。

OptimizationDesignofLiquidHybridBearinginGrinding

HeadofHigh-speedRollGrinderBasedonGeneticAlgorithm

WuHuaichao1,2LinghuKejun1Sun Guanchao1Li Zhe1Zhang Shunfeng1

1.GuizhouUniversity,Guiyang,550025

2.StateKeyLaboratoryofTribology,TsinghuaUniversity,Beijing,100084

Abstract:The performances of liquid hybrid bearing affect directly grinding performances of the whole roll grinder.Aiming at a kind of adjustable liquid hybrid bearing in grinding head of high-speed roll grinder,its structural characteristics and working principles were analyzed, and its structural dimensions and working parameters were calculated.Then,the lowest total power loss per unit loading capacity was regarded as an objective function of optimization design, and optimization design of the bearing was performed based on genetic algorithm.Optimization results are as follows:objective function value decreases from 0.0808 to 0.0678,total power loss of the bearing decreases from 766.9W to 642.6W,and power consumption is reduced about 16%;moreover,structural dimensions of the bearing are reduced,and its stiffness is obviously improved.Therefore,through the optimization design,manufacturing cost of the bearing is reduced,and its comprehensive performances are improved obviously.

Keywords:high-speedrollgrinder;liquidhybridbearing;optimizationdesign;geneticalgorithm

0引言

高速、超高速磨削技術(shù)是面向21世紀的先進

制造技術(shù),它的廣泛應用使機械加工方法產(chǎn)生變革性的影響[1]。其中,磨頭的軸承技術(shù)是高速、超高速磨削技術(shù)中最為關(guān)鍵的技術(shù)，它越來越受到廣大學者的高度關(guān)注[2]。

目前，高速、超高速機床較多采用的是陶瓷球軸承、磁懸浮軸承以及液體滑動軸承。陶瓷球軸承滾珠制造難度大、成本高,對拉伸應力和缺口應力敏感，在高速旋轉(zhuǎn)時滾珠會產(chǎn)生很大的離心力和陀螺力矩，限制了它在高速精密機床上的應用[3-5];磁懸浮軸承因為無機械接觸,所以無磨損、能耗小、噪聲小、壽命長、無需潤滑,但造價昂貴,承載能力相對較低，常用于高速輕切削機床主軸系統(tǒng)[6-8];液體滑動軸承分為液體靜壓軸承、液體動壓軸承以及液體動靜壓軸承，其中，液體動靜壓軸承因兼有液體動壓軸承和液體靜壓軸承的優(yōu)點，在全速度范圍內(nèi)能保持很高的承載能力、油膜剛度、主軸回轉(zhuǎn)精度以及很好的抗振性，而且，在主軸系統(tǒng)中采用冷卻措施后，軸承溫升和熱變形都能得到有效控制，因而，在高速軋輥磨床上得到了較多的應用[9-11]。

傳統(tǒng)液體滑動軸承的設(shè)計較多依賴經(jīng)驗和近似計算，并通過試驗對比來尋求較優(yōu)方案。對于高速軋輥磨床磨頭液體動靜壓軸承的設(shè)計，本文首先設(shè)計該軸承各部分的基本結(jié)構(gòu)，并確定其主要結(jié)構(gòu)尺寸和工作參數(shù)，在此基礎(chǔ)上，建立該軸承的優(yōu)化設(shè)計數(shù)學模型，并運用遺傳優(yōu)化算法對其進行優(yōu)化設(shè)計。

1高速液體動靜壓軸承的結(jié)構(gòu)

1.1軸承的基本結(jié)構(gòu)及工作原理

根據(jù)高速軋輥磨床磨頭的功能及結(jié)構(gòu)需求,其液體動靜壓軸承采用軸承和軸承套配合的結(jié)構(gòu)，并為了提高其應用范圍和運行效率，軸承上設(shè)計有調(diào)節(jié)機構(gòu)。該軸承的結(jié)構(gòu)如圖1所示。

1.前調(diào)節(jié)蓋　2.前調(diào)節(jié)推板　3.軸承過油孔 4.軸承本體　5.軸承套外螺紋　6.軸承套 7.密封圈　8.后調(diào)節(jié)蓋內(nèi)螺紋　9.后調(diào)節(jié)蓋 10.后調(diào)節(jié)推板　11.軸承套回油孔　12.軸承套進油孔 13.軸承進油孔　14.軸承靜壓腔 15.前調(diào)節(jié)推板回油孔 圖1　高速液體動靜壓軸承的結(jié)構(gòu)

如圖1所示,該軸承主要由軸承本體4、軸承套6、前調(diào)節(jié)推板2、后調(diào)節(jié)推板10、前調(diào)節(jié)蓋1和后調(diào)節(jié)蓋9組成。軸承4上均布有4個錐形肋板,與具有同樣錐度內(nèi)孔的軸承套6配合,從而強制軸承本體在每個錐形肋板處發(fā)生凹陷變形,于是在軸承本體內(nèi)壁和主軸之間形成8個楔形動壓腔,此凹陷變形量和動壓腔的形狀和大小可以通過軸承本體和軸承套的配合程度來進行調(diào)節(jié)。軸承套兩端均加有外螺紋5,分別與前調(diào)節(jié)蓋1和后調(diào)節(jié)蓋9上的內(nèi)螺紋8配合,在此螺紋副的作用下,調(diào)節(jié)主要是通過擰動前調(diào)節(jié)蓋和后調(diào)節(jié)蓋，并通過前調(diào)節(jié)推板7或后調(diào)節(jié)推板10的推動來實現(xiàn)的。

圖1所示的高速液體動靜壓軸承中,最關(guān)鍵的部件是軸承本體，圖2是軸承本體三維結(jié)構(gòu)圖。

1.靜壓腔　2.導流槽　3.錐形肋板 4.矩形肋板　5.過油孔　6.進油孔 圖2　高速液體動靜壓軸承本體的結(jié)構(gòu)

如圖2所示,在軸承本體內(nèi)部開設(shè)有4個靜壓腔和4個導流槽。其中,靜壓腔主要起兩方面的作用：一是在啟動時,向靜壓腔強力注入高壓油,使主軸懸浮于軸承中間,避免發(fā)生干摩擦;二是在停車時,因4個靜壓腔存在有高壓油，能有效減緩動壓突然消失而在重力作用下出現(xiàn)的干摩擦。導流槽主要起加強油液進入軸承內(nèi)壁后的流動和擴散作用。

1.2軸承的結(jié)構(gòu)尺寸和工作參數(shù)

圖1所示的高速液體動靜壓軸承的基本結(jié)構(gòu)尺寸及主要工作參數(shù)如表1所示。

表1　軸承的基本結(jié)構(gòu)尺寸及主要工作參數(shù)

基于表1所示數(shù)據(jù),為了進一步給后續(xù)軸承的優(yōu)化設(shè)計提供所需參數(shù),下面給出軸承的一些主要結(jié)構(gòu)尺寸和工作參數(shù)的計算公式[11]。

(1)每個靜壓油腔的有效承載面積Ae:

Ae=2R(I+Ia)sin((θ1+θ2)/2)

(1)

其中,R為軸承的半徑,I為油腔的軸向?qū)挾?θ1為軸承油腔包角的一半,θ2=45°。根據(jù)各參數(shù)值,可計算出Ae=5478mm2。

(2)軸承的剛度J:

J=6Aeps(β-1)cosθ1/(h0β2)

(2)

其中,β為節(jié)流比,毛細管節(jié)流β取為2。根據(jù)各參數(shù)值,可計算出J=9.5×108N/m。

(3)軸承的流量Q:

(3)

其中,pR為供油壓力ps的一半;η為潤滑油動力黏度,對于10號主軸油,在50℃時其動力黏度為η=0.0086Pa·s。根據(jù)各參數(shù)值，可計算出Q=3.32L/min。

(4)軸承的溫升Δt:

Δt=Ht/(QρCv)

(4)

其中,Ht為總功率;ρ為潤滑油的密度;Cv為潤滑油的質(zhì)量定容熱容。Ht的計算式為

Ht=Hf+Hp=(1+K)Hp

(5)

式中,Hf為摩擦功率;Hp為泵功率,且Hp=psQ;K為功率比,K=Hf/Hp。

將式(5)代入式(4),可得

Δt=(1+K)ps/(ρCv)

(6)

根據(jù)各參數(shù)值,可計算出Δt=7.6℃。

2高速液體動靜壓軸承的優(yōu)化設(shè)計

上述對高速軋輥磨床磨頭液體動靜壓軸承的設(shè)計僅是從功能的角度以及結(jié)構(gòu)方面的需求來進行的。為了提高該軸承的使用性能，需要對其結(jié)構(gòu)進行優(yōu)化設(shè)計，為此，首先要構(gòu)造該軸承優(yōu)化設(shè)計的數(shù)學模型。

2.1優(yōu)化設(shè)計的數(shù)學模型

2.1.1設(shè)計變量的選擇

圖1所示高速液體動靜壓軸承中,相對寬度L/D直接影響其靜態(tài)和動態(tài)性能;相對封油面尺寸Ia/L和It/D直接影響其有效承載面積、封油面上的動壓效應以及摩擦功耗，并且影響通過封油面的流量大小和油腔之間的周向內(nèi)部流動，從而對軸承的性能影響較大；軸承半徑間隙h0對于泵功率Hp和摩擦功率Hf影響較大,當h0增大時,泵功率Hp增大，Hf減小。綜上所述,選擇L/D、Ia/L、It/D和h0作為設(shè)計變量。因而,設(shè)計變量可表述為

X=[L/DIa/LIt/Dh0]T

(7)

2.1.2目標函數(shù)的確定

液體動靜壓軸承的發(fā)熱和溫升是影響磨頭正常運行的關(guān)鍵因素,協(xié)調(diào)好軸承的承載性能和總功率損失之間的關(guān)系是其優(yōu)化設(shè)計的核心,因此,以單位承載量下的總功率損失最小作為優(yōu)化設(shè)計的目標函數(shù)，即

F(X)=min(Ht/W)

(8)

式中,W為承載量。

Hf的計算公式如下：

Hf=FfU

(9)

其中,Ff為黏滯阻力,其計算式為

Ff=ηAfU/h0

(10)

U為軸頸處線速度,其計算式為

U=2πnR/60

(11)

式(10)中,Af為軸承的有效摩擦面積,其計算式為[2]

Af=(AL+AR/4)N

(12)

式中,AL為一個油腔四周封油面積;AR為一個油腔面積;N為油腔數(shù)目。

將式(10)～式(12)代入式(9),可得

Hf=π2ηD2n2Af/(3600h0)

(13)

泵功率Hp的計算公式如下：

(14)

承載量W即為軸承在穩(wěn)態(tài)運行時候承受的外載荷,圖1所示的高速軋輥磨床磨頭液體動靜壓軸承在穩(wěn)態(tài)運行下的承載量W為9488N。

2.1.3約束條件

根據(jù)液體動靜壓軸承一般的設(shè)計規(guī)則,前述選取的設(shè)計變量的取值范圍如下。

(1)軸承寬徑比L/D：

0.8≤L/D≤1.2

(15)

(2)軸向封油面尺寸Ia/L:

0.1≤Ia/L≤0.5

(16)

(3)周向封油面尺寸It/D:

0.1≤It/D≤0.25

(17)

(4)半徑間隙h0:

0.02mm≤h0≤0.04mm

(18)

2.2數(shù)學模型的求解

目前優(yōu)化算法主要有復合形法、可變?nèi)莶罘ê蛻土P函數(shù)法等傳統(tǒng)算法以及一些現(xiàn)代優(yōu)化計算方法，包括模擬退火優(yōu)化算法、遺傳優(yōu)化算法和神經(jīng)網(wǎng)絡(luò)優(yōu)化算法等[12]。由于動靜壓軸承的數(shù)學模型是屬于有約束的非線性數(shù)學模型，且目標函數(shù)問題比較復雜，因此本文采用遺傳優(yōu)化算法。

遺傳優(yōu)化算法是一種求解復雜問題并具有自適應、自組織和隨機優(yōu)化性質(zhì)的算法，它的基本數(shù)學模型可表示為[13]

F=f(C,E,P0,M,φ,Γ,Ψ,T)

(19)

式中,C為個體的編碼方法；E為個體的適應度評價函數(shù)；P0為初始種群；M為種群大小；φ為選擇算子；Γ為交叉算子；Ψ為變異算子；T遺傳運算終止條件。

根據(jù)遺傳算法的基本思想和上述基本數(shù)學模型，結(jié)合圖1所示的高速軋輥磨床磨頭液體動靜壓軸承的結(jié)構(gòu)特點，對該軸承采用遺傳算法進行優(yōu)化設(shè)計，流程如圖3所示。

圖3　基于遺傳算法的優(yōu)化設(shè)計流程圖

根據(jù)上述流程圖,在MATLAB中編寫程序，實現(xiàn)對該軸承的優(yōu)化,其核心程序代碼如下：

maxgen=2000;

sizepop=20;

pcross=[0.4];

pmutation=[0.1];

lenchrom=[1 1 1 1];

bound=[0.8 1.2;0.1 0.5;0.1 0.25;0.02 0.04];

individuals=struct('fitness',zeros(1,sizepop),'chrom',[]);

avgfitness=[];

bestfitness=[];

bestchrom=[];

fori=1:sizepop；

individuals.chrom(i,:)=Code(lenchrom,bound);

x=individuals.chrom(i,:);

individuals.fitness(i)=fun(x);

end

[bestfitnessbestindex]=min(individuals.fitness);

bestchrom=individuals.chrom(bestindex,:);

avgfitness=sum(individuals.fitness)/sizepop;

trace=[avgfitnessbestfitness];

fori=1:maxgen

individuals=select(individuals,sizepop);

avgfitness=sum(individuals.fitness)/sizepop;

individuals.chrom=Cross(pcross,lenchrom,individuals.chrom,sizepop,bound);

individuals.chrom=Mutation(pmutation,lenchrom,individuals.chrom,sizepop,[imaxgen],bound);

forj=1:sizepop

x=individuals.chrom(j,:);

individuals.fitness(j)=fun(x);

end

[newbestfitness,newbestindex]=min(individuals.fitness);

[worestfitness,worestindex]=max(individuals.fitness);

ifbestfitness>newbestfitness

bestfitness=newbestfitness;

bestchrom=individuals.chrom(newbestindex,:);

end

individuals.chrom(worestindex,:)=bestchrom;

individuals.fitness(worestindex)=bestfitness;

avgfitness=sum(individuals.fitness)/sizepop;

trace=[trace;avgfitnessbestfitness];

end

在上述主程序中,首先對遺傳算法的參數(shù)進行了初始化設(shè)置,設(shè)置進化代數(shù)為2000,種群規(guī)模為20,交叉概率選擇0.4,變異概率選擇0.1;然后設(shè)置優(yōu)化設(shè)計變量的取值范圍;最后初始化種群,調(diào)用目標函數(shù)子程序計算適應度,進行迭代尋優(yōu)。在調(diào)用目標函數(shù)子程序時，選擇、交叉和變異遺傳算子被迭代計算，從而計算出適應度，進而代替上一次進化中的最好的染色體。如此循環(huán)迭代計算，記錄了每一代進化中的最好的適應度，從而最終得到設(shè)計參數(shù)的優(yōu)化結(jié)果。

3高速液體動靜壓軸承優(yōu)化設(shè)計結(jié)果及其分析

在MATLAB軟件中運行優(yōu)化算法程序，得到了如表2所示的優(yōu)化計算結(jié)果。由表2的優(yōu)化結(jié)果可以看出：軸承的寬度優(yōu)化后,其值變小,結(jié)構(gòu)變小;半徑間隙減小,軸承的流量由原來的3.32L/min變化到2.4L/min,軸承的流量減小，從而降低了泵功率的損失；摩擦功率也由原來的545.1W降低到482.4W,從而導致軸承的總功率由初始設(shè)計的766.9W降為642.6W,功耗降低約16%;優(yōu)化后目標函數(shù)從0.0808下降到了0.0678。另外,軸承的剛度由原來的950MN/m變化為1030MN/m,軸承的剛度有所提升,使得軸承的穩(wěn)定性有所提高;軸承的溫升由原來的7.6℃變化為8.8℃,優(yōu)化前后溫升變化為1.2℃,對軸承的性能影響較小。

表2　高速液體動靜壓軸承優(yōu)化計算結(jié)果

4結(jié)論

針對一種高速軋輥磨床磨頭液體動靜壓軸承，為了改善其結(jié)構(gòu)和使用性能，采用遺傳優(yōu)化算法對其進行了優(yōu)化設(shè)計。優(yōu)化結(jié)果表明：軸承的結(jié)構(gòu)尺寸得到減小，從而節(jié)約了材料，降低了成本；半徑間隙得到減小，流量也得到減小，從而降低了泵功率損耗，達到了節(jié)能的目的；軸承的剛度得到提高，從而提升了軸承的穩(wěn)定性；軸承的溫升在優(yōu)化前后變化不大，對軸承的性能影響不大；目標函數(shù)值下降，從而達到了優(yōu)化設(shè)計的目的。由此可見，通過對軸承進行基于遺傳算法的優(yōu)化設(shè)計，其性能得到了明顯的提高。

參考文獻:

[1]周志雄,鄧朝暉,陳根余.磨削技術(shù)的發(fā)展關(guān)鍵技術(shù)[J].中國機械工程,2000,11(1/2):186-189.

ZhouZhixiong,DengZhaohui,ChenGenyu,etal.DevelopmentandKeyTechnologyofGrinding[J].ChinaMechanicalEngineering,2000,11(1/2):186-189.

[2]朱海港.超高速液體動靜壓混合軸承參數(shù)優(yōu)化設(shè)計研究[D]. 沈陽：東北大學,2009.

[3]喻煒.高速混合陶瓷球軸承分析與優(yōu)化設(shè)計[D].天津:天津大學,2011.

[4]WemhonerJ.Ceramics-AMilestoneontheWaytotheHigh-performanceRollingBearing[J].CFICeramicForumInternational,2002,79(4):21-24.

[5]HeQ,LiuHZ,ZhangYB,etal.TheCharacteristicsofHybridCeramicBallBearingforHigh-speedSpindle[C]//2009IEEEInternationalConferenceonMechatronicsandAutomation,Changchun,2009,2489-2494.

[6]劉迎澍,黃田.磁懸浮軸承研究綜述[J].機械工程學報,2000,36(11):5-9.

LiuYingshu,HuangTian.SurveyoftheResearchofMagneticBearings[J].ChineseJournalofMechanicalEngineering,2000,36(11)：5-9.

[7]SamantaP,HiraniH.MagneticBearingConfigurations:TheoreticalandExperimentalStudies[J].IEEETransactionsonMagnetics,2008,44(2):292-300.

[8]MusolinoA,RizzoR,TucciM,etal.ANewPassiveMaglevSystemBasedonEddyCurrentStabilization[J].IEEETransactionsonMagnetics,2009,45(3):984-987.

[9]MeruaneV,PascualR.IdentificationofNonlinearDynamicCoefficientsinPlainJournalBearings[J].TribologyInternational,2008,41(8):743-754.

[10]MatsudaK,KanemitsuY,KijimotoS.OptimalClearanceConfigurationofFluid-filmJournalBearingsforStabilityImprovement[J].JournalofTribology-transactionsoftheASME,2004,126(1):125-131.

[11]鐘洪,張冠坤.液體靜壓動靜壓軸承設(shè)計使用手冊[M].北京:電子工業(yè)出版社,2007.

[12]TavazoeiMS,HaeriM.AnOptimizationAlgorithmBasedonChaoticBehaviorandFractalNature[J].JournalofComputationalandAppliedMathematics,2007,206(2):1070-1081.

[13]雷英杰，張善文，李續(xù)武，等.MATLAB遺傳算法工具箱及應用[M].西安：西安電子科技大學出版社，2005.

(編輯袁興玲)