劉海亮 賈翠玲 王東輝
(內(nèi)蒙古工業(yè)大學工程訓練教學部,內(nèi)蒙古 呼和浩特010062)
數(shù)控加工技術是在傳統(tǒng)制造技術基礎上集成了自動控制技術、計算機技術和網(wǎng)絡通信技術而成的產(chǎn)物。與傳統(tǒng)機械制造比,具有自動化程度高、生產(chǎn)效率高、質(zhì)量高等諸多優(yōu)勢[1]。工藝參數(shù)作為數(shù)控加工的可控參數(shù),參數(shù)數(shù)值的設置極大地影響數(shù)控加工質(zhì)量,對參數(shù)進行優(yōu)化可以提高生產(chǎn)效率和生產(chǎn)質(zhì)量,同時降低生產(chǎn)成本,因此研究參數(shù)優(yōu)化問題具有重要的實際意義。
依據(jù)優(yōu)化目標的不同,數(shù)控加工技術的優(yōu)化可以分為能耗優(yōu)化、生產(chǎn)效率/成本優(yōu)化、刀具磨損/壽命優(yōu)化、生產(chǎn)質(zhì)量優(yōu)化等。文獻[2]面向數(shù)控機床能耗進行優(yōu)化,根據(jù)機床閑置時間制定機床待命、系統(tǒng)待命和關機3種控制方案,實現(xiàn)了數(shù)控機床的智能化低能耗控制。文獻[3]使用響應面法對金屬加工過程中的能耗和刀具使用壽命進行優(yōu)化,得到了刀尖半徑對目標參數(shù)的影響規(guī)律。文獻[4]為了降低機床熱誤差對數(shù)控機床加工精度的影響,使用布谷鳥算法優(yōu)化灰色GM(1,N)模型的系數(shù),提高了灰色模型預測精度,為熱誤差補償提供了依據(jù)。文獻[5]使用灰色關聯(lián)分析、主成分分析、響應面等方法對工件表面粗糙度和能耗進行優(yōu)化,得到了優(yōu)化的進給速度、切削深度、切削速度等參數(shù)。當前關于數(shù)控加工工藝優(yōu)化面臨2個難題:一是優(yōu)化目標經(jīng)驗公式精度略低,且不具有個體差異的適應性,難以滿足高要求的優(yōu)化問題;二是優(yōu)化模型的求解方法需進一步研究,力求得到最優(yōu)參數(shù),而不是次優(yōu)參數(shù)。
本文針對汽油機活塞的數(shù)控車削加工過程,以提高加工質(zhì)量和生產(chǎn)效率為目標建立了多目標優(yōu)化模型。使用非支配排序差分進行算法進行模型求解,并提出了基于加權相對距離的決策算法。經(jīng)過優(yōu)化和決策,實現(xiàn)了提高生產(chǎn)質(zhì)量和生產(chǎn)效率的目的。
本文以某汽油機活塞的數(shù)控車削加工過程為研究背景,對數(shù)控加工的產(chǎn)品質(zhì)量和加工速率進行優(yōu)化。汽油機活塞材料為鋁硅合金,活塞尺寸為φ81.5 mm×50 mm,經(jīng)車削和銑削后的最終產(chǎn)品如圖1所示,本文只研究車削過程的優(yōu)化問題。加工前毛坯為圓柱體,尺寸為φ85 mm×55 mm。
實驗用車床為FeelErftc-350型車床,如圖2所示。刀具為PCD活塞外圓專用刀具,刀具前角為5°,后角為11°,刀尖圓弧直徑為1.0 mm。
對于多目標優(yōu)化問題,其數(shù)學模型包括建立目標模型、設計優(yōu)化參數(shù)、設置約束條件等3個方面。
(1)建立目標模型。對1.1節(jié)設置的研究背景,本文以零件表面粗糙度Ra、材料去除速率Q為優(yōu)化目標。零件的表面粗糙度越小則零件質(zhì)量越好,因此零件表面粗糙度是一個最小化目標;材料去除速率越大則車削加工速率越大,因此材料去除速率是一個最大化目標。綜合以上分析,得到多目標優(yōu)化函數(shù)為:
obj.f=(minRa)∩(maxQ)
(1)
零件表面粗糙度有經(jīng)驗公式可以使用,為:
(2)
式中:af為每轉(zhuǎn)進給量,r為刀具半徑。
材料去除速率是衡量數(shù)控車床加工的重要指標,是指單位時間內(nèi)切除掉的材料體積,其完全由切削速度、背吃刀量和每轉(zhuǎn)進給量決定,即:
Q=ae·ap·af
(3)
式中:ae為切削速度,ap為背吃刀量。
在此需要強調(diào)的是,材料去除率計算公式為精確計算公式,可以直接應用。但是零件表面粗糙率為經(jīng)驗計算公式,計算精度有限。
(2)設計優(yōu)化參數(shù)。選擇工藝參數(shù)切削速度ae、背吃刀量ap、每轉(zhuǎn)進給量af為優(yōu)化參數(shù),則3個參數(shù)形成的優(yōu)化向量表示為x=(ae,ap,af)。
(3)設置約束條件。本次優(yōu)化問題的約束條件包括優(yōu)化參數(shù)取值范圍約束、零件質(zhì)量參數(shù)取值范圍約束、數(shù)控車床加工能力約束等,即:
(4)
式中:Fmax為最大切削力,Wmax為最大功率。
在1.2節(jié)中提到,材料去除速率為精確計算公式,可以直接使用;但是零件表面粗糙度是經(jīng)驗計算公式,其計算精度有限,且對每次加工、不同設備不具有自適應性和個體差異性。為了更加精確地對參數(shù)進行優(yōu)化,本文使用神經(jīng)網(wǎng)絡對工藝參數(shù)與目標參數(shù)間的關系進行映射。
簡單地講,BP神經(jīng)網(wǎng)絡地構(gòu)建需要3個步驟,分別為數(shù)據(jù)預處理、數(shù)據(jù)前向傳播、基于誤差反向傳播的參數(shù)修正。
(1)數(shù)據(jù)預處理。輸入輸出數(shù)據(jù)的量綱相差較大時會降低訓練速度,因此需要對輸入輸出數(shù)據(jù)進行歸一化處理,即
(5)
式中:xmin為參數(shù)最小值,xmax為參數(shù)最大值,x為參數(shù)原始數(shù)據(jù),x′為參數(shù)歸一化后數(shù)據(jù)。
(2)數(shù)據(jù)的正向傳播。數(shù)據(jù)正向傳播即為神經(jīng)網(wǎng)絡的輸出過程[7]。
(6)
式中:O為神經(jīng)網(wǎng)絡輸出;xi(i=1,2,3)為神經(jīng)網(wǎng)絡輸入,x1表示ae,x2表示ap,x3表示af;wij為輸入神經(jīng)元i與隱藏神經(jīng)元j的傳遞權值;θj為隱藏神經(jīng)元j的閾值;vj為隱藏神經(jīng)元j向輸出神經(jīng)元的權值,r為輸出神經(jīng)元閾值。f1為隱藏神經(jīng)元激活函數(shù),使用Tansig函數(shù);f1為輸出神經(jīng)元激活函數(shù),使用Tansig函數(shù)。
(3)參數(shù)訓練。首先計算擬合誤差,為:
(7)
式中:E為訓練誤差,U為數(shù)據(jù)長度,yi為樣本輸出,Oi為擬合輸出。
以傳遞權值wij為例,參數(shù)的訓練方法為:
(8)
式中:Δwij(t)為迭代t次時的參數(shù)修正值,η為學習效率。傳統(tǒng)神經(jīng)網(wǎng)絡一般使用固定的學習效率,但是學習效率過小容易陷入局部最優(yōu),過大則收斂速度過慢,因此本文設置自適應學習效率為:
(9)
使用Box-Behnken實驗設計法[8]在優(yōu)化空間內(nèi)進行抽樣,按照抽樣點對應的參數(shù)對汽油機活塞進行車削加工,而后測量加工件的表面粗糙度。得到32組實驗結(jié)果,如表1所示。
表1 實驗數(shù)據(jù)
使用式(8)和式(9)給出的訓練方法對神經(jīng)網(wǎng)絡參數(shù)進行訓練,訓練完畢后對網(wǎng)絡擬合準確性進行檢驗。將擬合準確性定義為:
(10)
式中:Acc為擬合準確性。
分別計算BP神經(jīng)網(wǎng)絡的擬合精度、經(jīng)驗公式(2)的計算精度,如表2所示。
表2 精度對比
由表2可以看出,BP神經(jīng)網(wǎng)絡的擬合準確度為99.11%,經(jīng)驗公式的計算精度為91.25%,遠小于BP神經(jīng)網(wǎng)絡的擬合精度。說明將BP神經(jīng)網(wǎng)絡應用于參數(shù)優(yōu)化??梢缘玫礁訙蚀_的優(yōu)化結(jié)果。
本文使用融合非支配排序的差分進化算法對多目標優(yōu)化模型進行求解。多目標優(yōu)化與單目標優(yōu)化的區(qū)別體現(xiàn)在基于非支配排序的選擇過程。
(1)染色體初始化。由于設置了3個優(yōu)化參數(shù),因此染色體具有3個基因位,即x=(ae,ap,af)。使用隨機方式進行初始化,為:
(11)
(2)變異操作。以差分的方式[9]對染色體執(zhí)行變異操作,即:
xij(n+1)=xOj(n)+F·[xPj(n)-xQj(n)]
(12)
式中:i≠O≠P≠Q(mào);F為縮放因子,n為迭代次數(shù)。當變異操作產(chǎn)生的個體超出取值范圍時,則重新選擇差分個體進行變異。
(3)交叉操作。在變異操作后執(zhí)行交叉操作,從而提高染色體多樣性。對于兩個相異的染色體xpj和xOj,交叉操作為:
(13)
式中:rand(0,1)為(0,1)間的隨機數(shù),CR為交叉概率。
(4)選擇操作。對于多目標優(yōu)化問題,基于非支配排序[10]和擁擠度因子進行選擇操作。其步驟為:首先,將父代染色體和子代染色體進行混合;而后,對混合染色體進行非支配排序;最后,基于非支配排序結(jié)果和擁擠度因子[11]對染色體進行選擇,保持染色體規(guī)模不變。
經(jīng)過3.1節(jié)的求解,得到多目標優(yōu)化的Pareto前沿解集,那么接下來存在解的選擇與決策問題。為了解決這一問題,本文提出了基于加權相對距離的決策方法,具體分為以下步驟:
Step1:目標參數(shù)歸一化,方法同式(5),以2維目標參數(shù)為例,歸一化結(jié)果寫為矩陣形式:
(14)
式中:N為Pareto解集中解的數(shù)量,(z11,z12)為解集中解1的2個目標取值,其余zij含義與此相同。
Step2:設置權重。為相對偏重的目標設置較大的權重ω,得到加權后的矩陣形式為:
本文將2個優(yōu)化目標視為同等重要,即取ω=1。
Step3:確定理想解與不理想解。所謂理想解是指由各目標最優(yōu)值組成的解,理論上此解無法實現(xiàn)。所謂非理想解是指由各目標最差值組成的解。以2維最小化目標為例,理想解與最不理想解如圖4所示。
Step4:計算相對距離。計算Pareto解集中某解i與理想解、非理想解的相對距離。
(15)
Step5:計算解i的優(yōu)劣指標參數(shù),為:
(16)
式中:Ci為解i的優(yōu)劣指標參數(shù),其值越大表示解的質(zhì)量越高,則選擇Ci值最大的解作為決策結(jié)果。
2.3節(jié)已經(jīng)對BP神經(jīng)網(wǎng)絡的擬合精度進行了驗證,本節(jié)只對多目標優(yōu)化結(jié)果和決策結(jié)果進行驗證。
根據(jù)1.1節(jié)的實驗方案和表1中的實驗數(shù)據(jù),使用非支配排序差分進化算法對多目標優(yōu)化模型進行求解,算法參數(shù)設置為:染色體規(guī)模為60,算法最大迭代次數(shù)為500。得到Pareto前沿解集如圖5所示。
按照3.2節(jié)的步驟計算Pareto解集中各解的優(yōu)劣指標參數(shù)Ci,所得結(jié)果如表3所示。
表3 Pareto中各解的指標值
根據(jù)表3的計算結(jié)果,選擇第17組解作為決策結(jié)果。另外,為了對本文的優(yōu)化方案進行說明,設置了2組比較方案,一是工廠現(xiàn)有加工方案,二是使用經(jīng)驗公式計算結(jié)合多目標求解與決策的方法?;?種生產(chǎn)方案均生產(chǎn)10個汽油機活塞,測量表面粗糙度并計算材料去除速率,對比結(jié)果如表4所示。
表4 不同生產(chǎn)方案比較
由表4可知,基于經(jīng)驗公式的優(yōu)化結(jié)果在表面粗糙度和材料去除速率2個方面均不如本文的優(yōu)化結(jié)果,這是因為經(jīng)驗公式的計算精度差于BP神經(jīng)的擬合精度,因此經(jīng)驗公式無法進行高精度的優(yōu)化,或者說基于經(jīng)驗公式的優(yōu)化結(jié)果可信度小于基于BP神經(jīng)網(wǎng)絡優(yōu)化的可信度。工廠現(xiàn)用生產(chǎn)方案中,生產(chǎn)工件的表面粗糙度均值為9.894 μm,材料去除速率為144.000 cm3/min,而本文優(yōu)后生產(chǎn)工件的表面粗糙度均值為4.267 μm,比工廠現(xiàn)用方案減小了1倍以上,材料去除速率為227.500 cm3/min,生產(chǎn)效率提高了57.98%。這是因為工廠的生產(chǎn)工藝確定與優(yōu)化是基于專家經(jīng)驗給出的,基于專家經(jīng)驗的優(yōu)化深度有限,且各專家的生產(chǎn)經(jīng)驗區(qū)別較大,因此工廠現(xiàn)用的生產(chǎn)工藝差于基于智能算法的優(yōu)化結(jié)果。
本文研究了汽油機活塞數(shù)控車削加工過程的優(yōu)化問題,建立了表面粗糙度、生產(chǎn)效率的多目標優(yōu)化模型,給出了多目標優(yōu)化與決策算法。經(jīng)驗證得出以下結(jié)論:(1)車削加工件表面粗糙度的經(jīng)驗公式計算精度差于BP神經(jīng)的擬合精度。(2)經(jīng)過優(yōu)化與決策,本文給出的生產(chǎn)方案不僅降低了工件表面粗糙度,而且提高了生產(chǎn)效率。