韓 軍 曹龍凱 徐 睿 姚 晟 熊鳳生

（①內(nèi)蒙古科技大學(xué)機械工程學(xué)院，內(nèi)蒙古 包頭 014010；②山東港口集團青島港國際股份有限公司前港分公司，山東 青島 266000）

隨著航天科工領(lǐng)域的快速發(fā)展，薄壁零件由于質(zhì)量輕和比強度高等優(yōu)點被頻繁地使用[1]，薄壁件在生產(chǎn)中的銑削參數(shù)選擇常常對銑削力和工件變形等有著重要影響，而實際銑削參數(shù)常常是技術(shù)人員根據(jù)經(jīng)驗選擇的，所用參數(shù)不是最優(yōu)銑削組合，導(dǎo)致加工過程中因銑削力過大等原因造成工件變形，不易保證被加工零件的公差尺寸，生產(chǎn)效率也不高[2]。因而研究薄壁件銑削參數(shù)優(yōu)化具有重要意義。

隨著智能優(yōu)化算法不斷的發(fā)展與完善，越來越多的學(xué)者將其用于銑削參數(shù)的優(yōu)化中。曾莎莎[3]等采用了倒傳神經(jīng)網(wǎng)絡(luò)的方法，優(yōu)化工藝參數(shù)組合。Liang M等[4]采用模糊控制方法，對銑削時的參數(shù)進(jìn)行了優(yōu)化。鄧偉[5]以最小能耗為目標(biāo)，基于粒子群算法開展銑削參數(shù)的優(yōu)化。雖然上述研究取得了預(yù)期效果，但所用算法較復(fù)雜，優(yōu)化過程較長，而MIGA算法不僅優(yōu)化速度快，且全局收斂快，計算效率高。姜佳明[6]采用了響應(yīng)曲面法對蜂窩薄壁件開展高速切削參數(shù)優(yōu)化，提高了切削加工效率，但響應(yīng)面模型的建立需要大量的樣本點，較為繁瑣，而Kriging模型在搭建過程時，所需樣本點數(shù)較少，且全局精度高。

針對上述問題，本文以某薄壁殼零件為研究對象，通過Isight平臺，采用了最優(yōu)拉丁超立方試驗方法建立樣本點，研究主軸轉(zhuǎn)速、徑向切深和軸向切深對銑削力的影響規(guī)律，擬合銑削參數(shù)與銑削力間的Kriging近似模型，以銑削力最小為目標(biāo)，采用MIGA算法優(yōu)化得出最合適的銑削參數(shù)組合。

1 試驗方案

零件的尺寸要求如圖1所示。

圖1 薄壁框二維圖

該薄壁零件毛胚材料為6061鋁合金，材料剛性較差。同時，因為毛坯自身殘余應(yīng)力的影響，加工后，殘余應(yīng)力得到了相應(yīng)的釋放，會引起零件特別是薄壁件的變形，故為了避免有過多的外力因素作用，更好的研究銑削力與銑削參數(shù)對薄壁件變形的影響，需要先對裝夾方式進(jìn)行確定。由于零件較薄，壁厚僅為1.5 mm，故不能采用虎鉗、真空吸盤等傳統(tǒng)裝夾方式。結(jié)合6061鋁合金的材料特點和薄壁零件的結(jié)構(gòu)特征，提出了采用在毛坯背面覆膠的裝夾方式，使其固定在工作臺上，這種方式避免了夾具與零件的直接接觸，只需要在加工中預(yù)留0.1 mm的余量，也就是總厚度為8 mm的薄壁零件，在外輪廓加工時，下切深度到7.9 mm，待加工完畢后，將整個毛胚取下，剩余的0.1 mm余量做去毛刺處理，裝夾如圖2所示。

圖2 薄壁件裝夾圖

確定裝夾方案后，對薄壁零件進(jìn)行銑削，加工完成后，使用惟德Tomoscope S型復(fù)合式坐標(biāo)測量儀進(jìn)行測量，發(fā)現(xiàn)其輪廓尺寸分別為109.13 mm和58.16 mm，不符合圖紙公差要求，通過多次銑削實驗，發(fā)現(xiàn)在半精加工階段時（主軸轉(zhuǎn)速3 000 r/min，徑向切深1 mm，軸向切深1 mm，銑刀直徑為10 mm、齒數(shù)為 4、前角 5°及后角 15°）因銑削參數(shù)選擇不合理造成銑削力過大，是導(dǎo)致薄壁零件變形量過大的主要原因之一。

2 試驗設(shè)計與結(jié)果

為獲取銑削參數(shù)與銑削力之間的近似模型，采用Isight試驗設(shè)計算法庫中的最優(yōu)拉丁超立方設(shè)計方法先行設(shè)計樣本點，該方法設(shè)計的樣本點具有均勻性好、可靠性高等優(yōu)點。根據(jù)現(xiàn)場加工經(jīng)驗選取銑削參數(shù)范圍：主軸轉(zhuǎn)速2 500～6 000 r/min，徑向切深0.5～2 mm，軸向切深0.5～3 mm，設(shè)計的30組樣本點和對應(yīng)的銑削力，如表1所示。

表1 試驗設(shè)計樣本點

3 基于Kriging-MIGA算法的優(yōu)化方案

基于Kriging-MIGA算法的優(yōu)化步驟如下：

（1）采用最優(yōu)拉丁超立方試驗設(shè)計方法設(shè)計主軸轉(zhuǎn)速、徑向切深和軸向切深3個銑削參數(shù)的樣本點，通過實驗得出樣本點對應(yīng)的銑削力。

（2）根據(jù)樣本點，搭建銑削參數(shù)與銑削力之間的Kriging近似模型。

（3）以銑削力最小為優(yōu)化目標(biāo)，采用MIGA算法對銑削參數(shù)進(jìn)行優(yōu)化。

（4）獲得銑削力的最小值，輸出對應(yīng)最優(yōu)銑削參數(shù)組合。

Kriging-MIGA算法流程如圖3所示：

圖3 Kriging-MIGA算法流程圖

4 銑削參數(shù)的優(yōu)化

4.1 近似模型的建立

Isight平臺中的Kriging近似模型屬于半?yún)?shù)化插值技術(shù)，在樣本點中對未知點插值具有很高的精度，對非線性關(guān)系擬合效果較好，常由回歸部分與隨機函數(shù)組合而成[7]。同時，由于Kriging的半?yún)?shù)化性,因而不需要搭建特定的函數(shù)模型[8]，相比其他近似模型技術(shù)，會更加方便與靈活。

給出樣本點x集合，對應(yīng)的響應(yīng)為y，樣本與響應(yīng)間的Kriging近似模型數(shù)學(xué)表達(dá)式如式（1）所示。

式中：β為回歸系數(shù)；fT(x)為默認(rèn)為高斯函數(shù)；z(x)為 隨機誤差；z(x)的統(tǒng)計特性[9]如式（2）所示。

式中：xi、xj為 樣本中任意兩點；R(θ，xi，xj)為表示樣本的空間相關(guān)性。

根據(jù)最優(yōu)拉丁超立方試驗設(shè)計的樣本點，通過Kriging近似模型方法，構(gòu)建銑削參數(shù)與銑削力之間的近似模型。

構(gòu)建Kriging近似模型時，銑削參數(shù)與銑削力之間至少需要2n+1個設(shè)計樣本點[10]，其中變量n為銑削參數(shù)的個數(shù)，由表1可知n=3，而本文設(shè)計了30個樣本點，故滿足分析要求，近似模型如圖4所示。

圖4 近似模型圖

為保證近似模型的精度，根據(jù)試驗方法選取樣本點開展誤差分析，采用可決系數(shù)評估擬合精度[11]，其表達(dá)式為

式中：N為試驗次數(shù)；為 預(yù)測值，yi為模擬值。

根據(jù)計算可知R2=0.935，當(dāng)可決系數(shù)大于0.9的時候[12]，Kriging近似模型達(dá)到精度要求，即模型是可靠的。

4.2 MIGA算法的優(yōu)化

MIGA算法屬于改進(jìn)后的并行遺傳算法，相比傳統(tǒng)的遺傳算法，不同之處在于MIGA算法把種群比作“島”，多個種群便有多個島，每個島上獨立地進(jìn)行遺傳算法尋優(yōu)，大大提升了全局優(yōu)化能力，其特點有：

（1）在使用MIGA算法時，需要將優(yōu)化問題的解進(jìn)行編碼，目的是為了形成類似染色體的鏈碼，每一個鏈碼即為1個個體，多個個體便組成種群。MIGA算法的編碼為格雷碼編碼[13]，如式（4）所示。

（2）MIGA算法將種群分成多個島，為了能選出島中最優(yōu)的個體，增加島內(nèi)個體的多樣性，島與島之間常會產(chǎn)生遷移、選擇等行為。

基于銑削參數(shù)與銑削力之間的近似模型，在Isight 平臺的Optimization模塊里采用MIGA算法對模型進(jìn)行尋優(yōu)。同時，算法優(yōu)化之前，以銑削力最小為優(yōu)化目標(biāo)，對相關(guān)銑削參數(shù)進(jìn)行約束，約束條件如式（5）所示。

式中：n為主軸轉(zhuǎn)速；ap為軸向切深；ae為徑向切深。

MIGA算法參數(shù)設(shè)置如表2所示。

表2 多島遺傳（MIGA）算法參數(shù)配置

相關(guān)參數(shù)與約束設(shè)置成功后，進(jìn)行優(yōu)化，迭代過程如圖5所示。

圖5 MIGA優(yōu)化迭代圖

MIGA算法在全局尋優(yōu)的過程中，相應(yīng)的目標(biāo)變量在不斷地收斂，雖然初始跳動較大，但達(dá)到80代后，目標(biāo)變量逐漸平穩(wěn)，偶爾因交叉變異出現(xiàn)較小的波動[14]，90代后多島遺傳算法出現(xiàn)最優(yōu)解。

5 優(yōu)化結(jié)果分析

銑削參數(shù)與銑削力在優(yōu)化前后的結(jié)果對比如表3所示。

表3 優(yōu)化結(jié)果對比

將優(yōu)化后的銑削參數(shù)組合（主軸轉(zhuǎn)速5 331 r/min，徑向切深0.91 mm，軸向切深0.75 mm）應(yīng)用到現(xiàn)場實驗加工中，銑刀與工件材料的選擇與優(yōu)化前一致，所得結(jié)果為269.79 N，而通過MIGA算法優(yōu)化所得的銑削力為244.15 N，兩者誤差為9.5%，誤差產(chǎn)生的原因是因為薄壁件實際銑削中會有熱變形以及機床振動等影響因素的存在，只要誤差在合理范圍內(nèi)，就可以認(rèn)為Kriging近似模型準(zhǔn)確擬合出了銑削參數(shù)與銑削力之間的關(guān)系。此時，銑削力從原先的444.95 N減小到269.79 N，減小了39.37%。

6 實驗及結(jié)果分析

將毛坯裝夾于德瑪吉DMU50數(shù)控加工中心，在半精加工階段采用優(yōu)化后的銑削參數(shù)組合進(jìn)行銑削，加工后工件如圖6所示。

圖6 加工工件

使用惟德Tomoscope S型復(fù)合式坐標(biāo)測量儀對工件進(jìn)行驗證，結(jié)果如表4所示。

表4 檢測結(jié)果 mm

由表4可知，鋁合金薄壁殼零件尺寸相比優(yōu)化前的尺寸均在公差要求范圍以內(nèi)。

7 結(jié)語

通過Isight軟件平臺，本文提出了一種基于Kriging-MIGA算法的薄壁零件銑削參數(shù)優(yōu)化方法。以6 061鋁合金薄壁殼零件為研究背景，搭建以銑削力最小為優(yōu)化目標(biāo)，主軸轉(zhuǎn)速、徑向切深和軸向切深為優(yōu)化變量的多目標(biāo)優(yōu)化模型，采用了MIGA算法對模型尋優(yōu)，分別將優(yōu)化前后的銑削力數(shù)值進(jìn)行對比分析，研究發(fā)現(xiàn)本文采用的優(yōu)化方法有效地降低了銑削力，減小了薄壁件的變形量，為現(xiàn)場技術(shù)人員選擇合理的銑削參數(shù)提供了理論依據(jù)。