遲玉倫 李郝林

(上海理工大學(xué)機(jī)械工程學(xué)院，上海200093)

高精密外圓磨削加工過(guò)程中，磨削工件的表面質(zhì)量往往受多方面因素影響。特別是難加工材料奧氏體不銹鋼Cr19Ni10 的磨削，在不同加工工藝參數(shù)下磨削的工件表面質(zhì)量會(huì)有很大不同。如何有效分析優(yōu)化各工藝參數(shù)大小，即調(diào)整控制砂輪線速度、磨削深度、進(jìn)給速度和工件線速度等參數(shù)數(shù)值提高工件質(zhì)量和磨削效率，是實(shí)現(xiàn)難加工材料高精度高效率磨削急需解決的重要問(wèn)題［1-3］。

隨著高精密磨削技術(shù)廣泛應(yīng)用于機(jī)械加工領(lǐng)域，許多國(guó)內(nèi)外學(xué)者對(duì)磨削工藝參數(shù)優(yōu)化進(jìn)行大量研究工作。Kuo-Ming Lee 等提出了一種新進(jìn)化算法對(duì)平面磨削工藝參數(shù)進(jìn)行優(yōu)化，取得了很好計(jì)算結(jié)果［4］;Kishalay Mitra 等提出了基于遺傳算法的各磨削工藝參數(shù)不確定性的多目標(biāo)優(yōu)化模型對(duì)磨削過(guò)程進(jìn)行仿真優(yōu)化［5］;合肥工業(yè)大學(xué)孫林等在基于最小二乘支持矢量機(jī)的成形磨削表面粗糙度預(yù)測(cè)模型基礎(chǔ)上，提出了磨削工藝參數(shù)優(yōu)化設(shè)計(jì)的可行性方案［6］。由于磨削過(guò)程影響因素較為復(fù)雜，其預(yù)測(cè)模型與磨削工藝參數(shù)優(yōu)化設(shè)計(jì)仍有許多問(wèn)題需要解決。

隨著灰色系統(tǒng)理論技術(shù)在機(jī)械加工領(lǐng)域的創(chuàng)新應(yīng)用，為解決傳統(tǒng)加工問(wèn)題提供了有效解決辦法?；疑到y(tǒng)是一種少數(shù)據(jù)不確定性的系統(tǒng)，針對(duì)很多信息不完備、不確定，數(shù)據(jù)較少的問(wèn)題，可以應(yīng)用灰色系統(tǒng)理論進(jìn)行研究。灰關(guān)聯(lián)分析模型可有效解決因子間關(guān)聯(lián)的相對(duì)度量，提供相對(duì)客觀的評(píng)價(jià)指標(biāo)權(quán)重量度的方法［7-8］。所以，本文基于灰色系統(tǒng)理論方法，選用外圓磨削中的砂輪線速度、進(jìn)給速度、磨削深度和工件線速度4 個(gè)工藝參數(shù)作為實(shí)驗(yàn)研究因素，將磨削工件的表面粗糙度、磨削工件圓度作為試驗(yàn)評(píng)價(jià)指標(biāo)。通過(guò)正交試驗(yàn)方法來(lái)合理安排各因素與評(píng)價(jià)指標(biāo)的實(shí)驗(yàn)設(shè)計(jì)，應(yīng)用灰色理論方法對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行處理計(jì)算。分析工件表面粗糙度與磨削工件圓度之間的對(duì)應(yīng)關(guān)系、不同磨削工藝參數(shù)對(duì)磨削工件質(zhì)量的影響顯著性大小，及滿足上述實(shí)驗(yàn)中工件表面粗糙度值最小時(shí)的磨削工藝參數(shù)組合。

1 灰色理論

應(yīng)用灰關(guān)聯(lián)度分析計(jì)算步驟如下:

(1)由原始序列集確定參考序列(k)和比較序列(k)。

(2)初值化生成。由于不同的評(píng)價(jià)指標(biāo)使用的量綱和單位不一致，在比較的時(shí)候很難得到正確的結(jié)果。為了便于分析，以保證各因素具有等效性和同序性，應(yīng)將評(píng)價(jià)指標(biāo)無(wú)量綱化，即做規(guī)范化處理。方法如下。

可將評(píng)價(jià)指標(biāo)集劃分為極大型和極小型兩大類指標(biāo)，對(duì)于級(jí)大型指標(biāo)。其值越大越好，如下式［3］:

對(duì)于極小型指標(biāo)，其值越小越好，如下式:

此外，一種簡(jiǎn)單常用初值化生成方法，是用序列x的初始值x(1)去除序列x中每一個(gè)數(shù)，以獲得初值化生成序列，如下式:

式中:i=1，2，…，m;k=1，2，…，m為實(shí)驗(yàn)的次數(shù);n為比較序列的個(gè)數(shù);(k)為原始比較序列，(k)為處理后的比較序列;max(k)、min(k)分別為(k)中的最大值和最小值。

(3)灰色關(guān)聯(lián)系數(shù)計(jì)算［4］。比較序列(k)對(duì)于參考序列(k)在第k個(gè)點(diǎn)的灰色關(guān)聯(lián)系數(shù)計(jì)算，如下式:

其中:

(4)計(jì)算灰關(guān)聯(lián)度。令γi表示灰關(guān)聯(lián)度，如下式計(jì)算

在實(shí)際應(yīng)用中，由于各指標(biāo)之間重要程度的差異，必須考慮相應(yīng)的權(quán)重?？捎弥饔^或客觀賦權(quán)法給定評(píng)價(jià)指標(biāo)間的權(quán)值ωk，則灰關(guān)聯(lián)度計(jì)算如下式:

式中:ωk表示給定第k個(gè)指標(biāo)的權(quán)值;假設(shè)每個(gè)指標(biāo)的權(quán)值相等，則式(6)等同于式(5)。

灰關(guān)聯(lián)度是比較序列與參考序列之間的比較測(cè)度，如果比較序列與參考序列相同，則灰關(guān)聯(lián)度的值為1?；谊P(guān)聯(lián)度也反映了各比較序列之間對(duì)參考序列的重要程度，如果一個(gè)比較序列比其他序列對(duì)參考序列的影響程度大，則該比較序列計(jì)算的灰關(guān)聯(lián)度也越大［8，11］。

2 實(shí)驗(yàn)研究

如圖1 所示，實(shí)驗(yàn)機(jī)床選用為Schleifring K-C33 精密數(shù)控外圓磨床，磨削砂輪為單晶剛玉(SA)砂輪，磨削液選用水溶性磨削油TY101S，其濃度配比為4%，實(shí)驗(yàn)所使用的工件材料為奧氏體不銹鋼Cr19Ni10。該實(shí)驗(yàn)?zāi)康氖欠治鲅芯烤ルA段不同磨削深度、進(jìn)給速度、砂輪線速度和工件線速度等工藝參數(shù)對(duì)磨削工件表面粗糙度和工件圓度的影響，進(jìn)而對(duì)各工藝參數(shù)進(jìn)行優(yōu)化。實(shí)驗(yàn)中所選定的工藝參數(shù)如表1 所示。

圖1 磨削試驗(yàn)裝置

表1 工藝參數(shù)水平表

本文設(shè)置不同工藝參數(shù)磨削后的工件評(píng)價(jià)指標(biāo)為:(1)表面粗糙度;(2)工件圓度。磨削后的工件粗糙度通過(guò)使用Talysurf PGI1240 粗糙度儀測(cè)量，如圖2所示，選擇磨削工件不同位置測(cè)量3 次后求其粗糙度平均值Ra。磨削工件圓度測(cè)量，如圖3 所示，所使用測(cè)量?jī)x器為Tayloround 595，設(shè)置倍率為20 000，波數(shù)為1 ～15，選擇磨削工件不同位置測(cè)量3 次后求其圓度平均值Rness。

圖2 測(cè)量工件表面粗糙度

圖3 測(cè)量磨削工件圓度

如果將上述實(shí)驗(yàn)中每一種參數(shù)組合都進(jìn)行磨削測(cè)試，共需要81 次實(shí)驗(yàn)，如此多實(shí)驗(yàn)就需要很大經(jīng)濟(jì)成本和很長(zhǎng)時(shí)間來(lái)完成，顯然不符合實(shí)際情況。本文使用一種科學(xué)安排與分析多因素的正交試驗(yàn)設(shè)計(jì)方法，合理地構(gòu)造參數(shù)的不同組合進(jìn)行實(shí)驗(yàn)測(cè)試。根據(jù)因素與各因素水平分析，本試驗(yàn)中選擇L9(34)正交試驗(yàn)方案進(jìn)行安排各實(shí)驗(yàn)測(cè)試，如表2 所示。

表2 正交試驗(yàn)設(shè)計(jì)

表3 處理后序列

3 實(shí)驗(yàn)結(jié)果計(jì)算與分析

3.1 實(shí)驗(yàn)結(jié)果計(jì)算

如表2 所示，選用工件表面粗糙度和工件圓度作為各磨削工藝參數(shù)每次試驗(yàn)的評(píng)價(jià)指標(biāo)。因?yàn)閷?shí)驗(yàn)中工件圓度粗糙度和工件圓度值越小越好，屬于極小型指標(biāo)，所以采用式(2)進(jìn)行初值化處理計(jì)算。設(shè)定最小的粗糙度和工件圓度為參考序列(k)，k=1，2;9 次試驗(yàn)的結(jié)果為比較序列(k)，i=1，…，9，k=1，…，9。計(jì)算結(jié)果及確定的參考序列如表3 所示。

第一次實(shí)驗(yàn)的比較序列與參數(shù)序列在第k點(diǎn)的絕對(duì)差計(jì)算如下:

用上述方法計(jì)算每次實(shí)驗(yàn)的比較序列與參數(shù)序列的絕對(duì)差，如表4 所示。其中兩級(jí)絕大差和兩級(jí)絕小差的計(jì)算值如下:

表4 絕對(duì)差值表

設(shè)該實(shí)驗(yàn)中各評(píng)價(jià)指標(biāo)的權(quán)重相同，即分辨系數(shù)ζ=0.5，將該分辨參數(shù)帶入式(4)和式(5)中計(jì)算每次實(shí)驗(yàn)的比較序列與參考序列在k點(diǎn)的灰關(guān)聯(lián)系數(shù)及灰關(guān)聯(lián)度，如表5 所示。

3.2 分析討論

灰關(guān)聯(lián)度是反應(yīng)比較序列與參考序列之間的比較測(cè)度，如果灰關(guān)聯(lián)度越大表示該序列與參考序列之間的相互關(guān)聯(lián)度越大。在本實(shí)驗(yàn)中，灰關(guān)聯(lián)度越大則表示工件表面粗糙度及工件圓度值越小，則工件磨削質(zhì)量越好。

根據(jù)上述計(jì)算結(jié)果如圖4 所示，不同工藝參數(shù)下9 次正交試驗(yàn)所對(duì)應(yīng)的粗糙度與工件圓度的灰關(guān)聯(lián)系數(shù)的曲線圖，表面粗糙度與工件圓度數(shù)值變化總趨勢(shì)具有一定的相似性。在圖6 中9 次實(shí)驗(yàn)中，第1 次正交試驗(yàn)的工藝參數(shù)下對(duì)應(yīng)的灰關(guān)聯(lián)度最大，即工件質(zhì)量最好。

圖4 9 次正交試驗(yàn)的粗糙度與工件圓度的灰關(guān)聯(lián)度曲線

根據(jù)上述正交試驗(yàn)的灰關(guān)聯(lián)度計(jì)算結(jié)果表5，可計(jì)算出每個(gè)磨削工藝參數(shù)不同水平所對(duì)應(yīng)的灰關(guān)聯(lián)度。該計(jì)算方法是計(jì)算每個(gè)工藝參數(shù)水平對(duì)應(yīng)所有灰關(guān)聯(lián)度的平均值，如計(jì)算切削深度水平1 的灰關(guān)聯(lián)度:

使用上述方法可將每個(gè)磨削工藝參數(shù)不同水平的灰關(guān)聯(lián)度計(jì)算出來(lái)，如表6 所示。得到每個(gè)工藝參數(shù)各水平最大灰關(guān)聯(lián)度與最小灰關(guān)聯(lián)度之差，該差值越大表示該工藝參數(shù)對(duì)磨削質(zhì)量的影響越顯著，如磨削深度的灰關(guān)聯(lián)度差值為0.259 1、工件進(jìn)給速度為0.353 6、砂輪線速度為0.033 5 及工件線速度為0.029 8。通過(guò)不同工藝參數(shù)的灰關(guān)聯(lián)度差值比較得到磨削工件進(jìn)給速度對(duì)工件質(zhì)量的影響最為顯著和工件線速度的影響顯著最小。

由表6 可知，磨削深度A和工件進(jìn)給速度B越小，灰關(guān)聯(lián)度就越大，則磨削質(zhì)量越好;而砂輪線速度C和工件線速度D分別是在2 水平上對(duì)應(yīng)的灰關(guān)聯(lián)度最大，說(shuō)明砂輪線速度和工件線速度過(guò)大或過(guò)小都會(huì)降低工件的表面質(zhì)量。將上述工藝參數(shù)優(yōu)化后得到最優(yōu)參數(shù)組合為磨削深度(水平1)、工件進(jìn)給速度(水平1)、砂輪線速度(水平2)和工件線速度(水平2)，即A1、B1、C2 和D2。選取上述優(yōu)化后的磨削工藝參數(shù)進(jìn)行磨削實(shí)驗(yàn)，如表7 所示。該實(shí)驗(yàn)結(jié)果得到較好的工件粗糙度和圓度值，驗(yàn)證了上述方法的可行性和有效性。

表5 灰關(guān)聯(lián)度

表6 各參數(shù)不同水平的平均灰關(guān)聯(lián)度

表7 優(yōu)化磨削參數(shù)實(shí)驗(yàn)結(jié)果

4 結(jié)語(yǔ)

本文選用正交試驗(yàn)方法對(duì)選定的工藝參數(shù)和評(píng)價(jià)指標(biāo)進(jìn)行實(shí)驗(yàn)設(shè)計(jì)，使用灰理論方法研究了各工藝參數(shù)與工件磨削質(zhì)量之間的關(guān)系，主要研究結(jié)論如下:

(1)分析了各工藝參數(shù)與表面粗糙度及工件圓度的對(duì)應(yīng)關(guān)系，研究得出上述9 次正交試驗(yàn)設(shè)計(jì)中的最優(yōu)參數(shù)組合為磨削深度(水平1)、工件進(jìn)給速度(水平1)、砂輪線速度(水平2)和工件線速度(水平2)。

(2)通過(guò)灰理論研究該實(shí)驗(yàn)中各工藝參數(shù)對(duì)工件表面粗糙度和工件圓度的影響顯著程度，得出對(duì)磨削質(zhì)量影響的顯著程度依次是:工件進(jìn)給速度、磨削深度、砂輪線速度及工件線速度。

(3)通過(guò)上述灰色關(guān)聯(lián)度優(yōu)化的磨削工藝參數(shù)進(jìn)行磨削實(shí)驗(yàn)，如表7 所示，該實(shí)驗(yàn)結(jié)果得到較好的工件粗糙度和圓度值，驗(yàn)證了上述方法的可行性和有效性。

［1］Xu W，Wu Y，Sato T，et al Effects of process parameters on workpiece roundness in tangential-feed centerless grinding using a surface grinder［J］. Journal of Materials Processing Technology.2010，210 (5):759 -766.

［2］Antti Mantynen，Alexey Zakharov，Sirkka-Liisa Jamsa-Jounela.Optimization of grinding parameters in the production of colorant paste［J］. Power Technology，2012，217:216 -222.

［3］Sanjay Agarwal，Venkateswara Rao P. Grinding characteristic，material removal and damage formation mechanisms in high removal rate grinding of silicon carbide［J］. International Journal of Machine Tools and Manufacture，2010，50:1137 -1146.

［4］Lee Kuo-Ming，Hsu Ming-Ren，Chou Jyh-Hong，et al. Improved differential evolution approach for optimization of surface grinding process［J］.Expert Systems with Applications，2011，38:5680 -5686.

［5］Kishalay Mitra. Multiobjective optimization of an industrial grinding operation under uncertainty［J］. Chemical Engineering Science，2009(64):5043 -5056.

［6］孫林，楊世元. 基于最小二乘支持矢量機(jī)的成形磨削表面粗糙度預(yù)測(cè)及磨削用量?jī)?yōu)化設(shè)計(jì)［J］. 機(jī)械工程學(xué)報(bào)，2009，45(10):254-258.

［7］Lin Yan-Cherng，Lee Ho-Shiun. Optimizationof machining parameters using magnetic-force-assisted EDM based on gray relational analysis［J］.Int J Adv Manuf Technol，2009(42):1052 -1064.

［8］李郝林，王健. 基于灰色關(guān)聯(lián)分析的平面磨削工藝參數(shù)優(yōu)化［J］.中國(guó)機(jī)械工程，2011，22(6):631 -635.

［9］鄧聚龍. 灰理論基礎(chǔ)［M］. 武漢:華中科技大學(xué)出版社，2002.

［10］李雁，李玉忱，任磊. 利用灰關(guān)聯(lián)理論定量評(píng)價(jià)選擇供應(yīng)商的算法研究［J］. 計(jì)算機(jī)應(yīng)用與軟件，2004，21(8):90 -98.

［11］Ghani J A，Choudhury I A，Hassan H H. Applicationof taguchi method in the optimization of end milling parameters［J］. Journalof Materials Processing Technology，2004(145):84 -92.