基于CFD-DEM耦合的磨粒流微小孔加工數(shù)值分析與試驗(yàn)

李俊燁，蘇寧寧，胡敬磊，楊兆軍，盛 亮，張心明※

（1. 長(zhǎng)春理工大學(xué)機(jī)電工程學(xué)院，長(zhǎng)春 130022；2. 吉林大學(xué)機(jī)械科學(xué)與工程學(xué)院，長(zhǎng)春 130012）

0 引 言

磨粒流加工技術(shù)是20世紀(jì)60年代發(fā)展起來(lái)的一種精密加工方法，它利用磨粒流自身流動(dòng)性的特點(diǎn)可加工任何復(fù)雜零部件，幾乎沒(méi)有尺寸限制，可顯著降低被加工工件的表面粗糙度，得到理想的表面形貌[1-3]。磨粒流加工的基本原理是以液體為載體，以無(wú)數(shù)的顆粒為切削刀具，在外界作用下對(duì)所加工表面進(jìn)行微切削和微摩擦，從而達(dá)到精密加工的目的[4-7]。

目前國(guó)內(nèi)外學(xué)者對(duì)磨粒流加工技術(shù)進(jìn)行了大量研究，探討影響因素，從而獲取最佳加工參數(shù)。李俊燁等[8]以伺服閥噴嘴為研究對(duì)象，從沖蝕磨損的角度分析了不同磨粒硬度下磨粒流加工效果，結(jié)果顯示碳化硅顆粒加工效果優(yōu)于白剛玉，工件表面質(zhì)量得到明顯提高。計(jì)時(shí)鳴等[9]以表面覆有氧化層的硅片為研究對(duì)象，在磨粒流加工過(guò)程中利用超聲空化作用加快了材料的去除，并搭建試驗(yàn)平臺(tái)，驗(yàn)證了超聲空化可以改變磨粒的動(dòng)能，提高切削效率。Uhlmann等[10]采用麥克斯韋模型，對(duì)磨粒流加工過(guò)程進(jìn)行CFD模擬，對(duì)進(jìn)出口的壓力降和速度分布等參數(shù)進(jìn)行分析，發(fā)現(xiàn)材料去除率與壓力差以及速度大小成正比。Butola等[11]總結(jié)了循環(huán)次數(shù)、擠壓壓力和磨料濃度 3個(gè)參數(shù)對(duì)磨粒流加工的影響，并利用田口試驗(yàn)方法獲得最佳參數(shù)值。通過(guò)試驗(yàn)發(fā)現(xiàn)：當(dāng)循環(huán)次數(shù)為6、擠壓壓力為1.5 MPa和磨料濃度為100 g時(shí)獲得最佳加工表面質(zhì)量。周迪鋒等[12]采用雙入口磨粒流加工裝置，結(jié)合DEM和CFD方法，應(yīng)用Abaqus中Johnson-Cook模型分析了撞擊速度和撞擊次數(shù)對(duì)靶材質(zhì)量損失的影響并進(jìn)行試驗(yàn)分析，結(jié)果表明表面粗糙度明顯下降，加工效率提高。計(jì)時(shí)鳴等[13]運(yùn)用 CFD-DEM 耦合的方法探究了面約束軟性磨粒流加工的均勻性，發(fā)現(xiàn)低黏度流體下材料去除均勻性有明顯提升。喻黎明等[14-15]采用歐拉-拉格朗日湍流模型對(duì)水力旋流器和迷宮流道內(nèi)的水沙流動(dòng)進(jìn)行 CFD-DEM 模擬，分析單個(gè)沙粒的軌跡線、速度和沙粒群的運(yùn)動(dòng)規(guī)律等，模擬結(jié)果表明沙粒在迷宮流道中會(huì)增加沙粒的碰撞次數(shù)，降低運(yùn)動(dòng)速度，增加堵塞機(jī)率；在水力旋流器中發(fā)現(xiàn)沙粒越小越難以分離，通過(guò)單個(gè)沙粒和沙粒群的運(yùn)動(dòng)可知沙粒既做圓周運(yùn)動(dòng)又做直線運(yùn)動(dòng)。該研究方法成為迷宮流道設(shè)計(jì)和水力旋流器性能研究一種有效手段。

磨粒流加工的實(shí)質(zhì)是磨粒流中的顆粒對(duì)所加工零件表面的凸起部分的碰撞與剪切，從而得到平滑的表面形貌[16-19]。當(dāng)前關(guān)于磨粒流加工大部分的研究?jī)?nèi)容主要體現(xiàn)在機(jī)床參數(shù)和磨料配置[20-22]，而忽略了顆粒在整個(gè)加工過(guò)程中的決定性作用，因此很少有人將顆粒列為重點(diǎn)研究對(duì)象。本文基于前人研究成果，采用離散元方法，基于 CFD-DEM 耦合方法，以點(diǎn)膠頭為研究對(duì)象，考慮顆粒對(duì)壁面的碰撞作用，探討流體和顆粒在加工過(guò)程中的分布狀態(tài)，對(duì)不同入口速度條件下的流體和顆粒下的分布狀態(tài)進(jìn)行對(duì)比分析，揭示磨粒流微切削作用行為；建立顆粒對(duì)壁面的碰撞模型，闡明磨粒流對(duì)壁面的作用規(guī)律和材料去除機(jī)理，深入分析了顆粒對(duì)壁面的切削作用。

1 磨粒流數(shù)學(xué)模型

1.1 流體方程

磨粒流在加工過(guò)程中的流動(dòng)為不可壓縮流動(dòng)，其流體的連續(xù)性方程和動(dòng)量方程可如下描述[23]

式中q代表為流體相，qr為流體相密度，kg/m3；為流體相速度，m/s；為固體相到流體相的傳質(zhì)，p為混合相的壓力，Pa；為兩相之間的相互作用，為外部體積力，為壓力應(yīng)變張量。

1.2 顆粒方程

在磨粒流運(yùn)動(dòng)過(guò)程中，顆粒受到多個(gè)力的作用，比如浮力、慣性力、拖曳力和虛擬質(zhì)量力等[24]，所以在計(jì)算中會(huì)求解很多方程，導(dǎo)致了求解過(guò)程的復(fù)雜。由于在高速流動(dòng)過(guò)程中拖曳力和重力起主要作用，所以除了重力和拖曳力之外，其他的可以忽略[25]，簡(jiǎn)化后的方程為

式中m是流體速度，m/s；pm是顆粒速度，m/s；xg是補(bǔ)充加速度，m/s2（作用于顆粒每單位質(zhì)量力）；pt是顆粒弛豫時(shí)間，s；r是流體密度，kg/m3；pr是顆粒密度，kg/m3。

2 點(diǎn)膠頭物理模型和網(wǎng)格劃分

本文選擇全自動(dòng)噴膠機(jī)的配件點(diǎn)膠頭為研究對(duì)象，通過(guò)點(diǎn)膠頭將膠水噴在產(chǎn)品上，可以保證噴涂厚度和形狀，但由于點(diǎn)膠頭直徑最大為2 mm，傳統(tǒng)的加工方法如珩磨、精研和磨削等[26-28]很難進(jìn)行加工，且價(jià)格昂貴、效率低，耗時(shí)長(zhǎng)[29-30]，因此本文采用磨粒流加工技術(shù)進(jìn)行精密拋光，該拋光技術(shù)不會(huì)破壞工件，且能夠完全接觸被加工表面，降低表面粗糙度和殘余應(yīng)力，得到良好的表面性能。點(diǎn)膠頭材質(zhì)為不銹鋼，其實(shí)物圖如圖1a所示，幾何尺寸如圖1b所示。將其流道抽取出來(lái)得到三維計(jì)算流域，利用 ICEM 軟件進(jìn)行網(wǎng)格劃分，劃分后的網(wǎng)格如圖1c所示。液相流體選擇航空煤油，固相選擇碳化硅，體積分?jǐn)?shù)設(shè)為 10%[31]，基于離散元方法，通過(guò)CFD-DEM耦合對(duì)流場(chǎng)進(jìn)行求解分析。

由于點(diǎn)膠頭內(nèi)徑不是一致的，因此為了方便對(duì)后面的參數(shù)進(jìn)行分析和敘述，呈現(xiàn)點(diǎn)膠頭內(nèi)部流場(chǎng)的變化趨勢(shì)，需要對(duì)點(diǎn)膠頭進(jìn)行區(qū)域劃分，將 3個(gè)不同內(nèi)徑的位置分別劃分為1區(qū)、2區(qū)和4區(qū)，并將2區(qū)和4區(qū)的連接處劃分為3區(qū)，因?yàn)閺娜肟诘匠隹冢?區(qū)位置的截面變化最大，點(diǎn)膠頭內(nèi)部流場(chǎng)變化也最明顯，劃分結(jié)果如圖2所示。

圖1 點(diǎn)膠頭結(jié)構(gòu)尺寸和網(wǎng)格劃分Fig.1 Structure size and gridding division of dispensing head

圖2 點(diǎn)膠頭結(jié)構(gòu)區(qū)域劃分Fig.2 Structure division of dispensing head

3 磨粒流加工點(diǎn)膠頭的數(shù)值分析

3.1 參數(shù)設(shè)置

CFD設(shè)置：流體相為航空煤油，選用RNGk–ε湍流模型，采用速度進(jìn)口，自由出口，壁面條件選擇增強(qiáng)壁面函數(shù)和無(wú)滑移邊界條件，動(dòng)量和湍流動(dòng)能采用二階迎風(fēng)格式，采用SIMPLE算法進(jìn)行求解。

DEM設(shè)置：顆粒相為碳化硅，進(jìn)出口條件與流體相一致，模擬過(guò)程為瞬態(tài)模擬，開啟Track Collision，時(shí)間步長(zhǎng)為2E-7 s，模擬時(shí)間為1 s。

3.2 不同入口速度條件下的流體動(dòng)壓與顆?？偰芰狂詈蠄?chǎng)分析

根據(jù)實(shí)際加工條件，設(shè)置磨粒粒徑為200目（75μm），體積分?jǐn)?shù)設(shè)置為10%，選取入口速度為30、35、40和45 m/s進(jìn)行 CFD-DEM 耦合數(shù)值模擬分析，通過(guò)流體力學(xué)軟件FLUENT得到不同入口速度條件下CFD-DEM耦合場(chǎng)下的流體動(dòng)壓和顆?？偰芰咳鐖D3所示。

圖3 不同入口速度的流體動(dòng)態(tài)壓強(qiáng)與顆?？偰芰縁ig.3 Fluid dynamic pressure and total particle energy at different inlet velocities

從圖 3流體動(dòng)態(tài)壓強(qiáng)與顆粒總能量可以看出，磨粒流進(jìn)入 1區(qū)時(shí)，流體動(dòng)態(tài)壓強(qiáng)和顆粒總能量基本保持不變，流體動(dòng)態(tài)壓強(qiáng)和顆?？偰芰慷甲钚?，在磨粒流進(jìn)入點(diǎn)膠頭 2區(qū)時(shí)，流體動(dòng)態(tài)壓強(qiáng)和顆粒總能量開始增大，在3區(qū)，流體動(dòng)態(tài)壓強(qiáng)和顆?？偰芰窟M(jìn)一步增大，到達(dá)4區(qū)時(shí)，流體動(dòng)態(tài)壓強(qiáng)和顆?？偰芰炕颈３植蛔?，這是由于磨粒流速度因?yàn)辄c(diǎn)膠頭內(nèi)徑的變小而增大，流體動(dòng)態(tài)壓強(qiáng)逐漸增大，顆粒受到流體曳力作用，顆粒動(dòng)能增大，顆?？偰芰侩S之增大，顆?？偰芰吭酱髮?duì)點(diǎn)膠頭壁面碰撞越激烈，材料去除量越大，有利于對(duì) 3區(qū)進(jìn)行去毛刺，從而實(shí)現(xiàn)對(duì)點(diǎn)膠頭進(jìn)行光整加工。

圖 4為流體動(dòng)態(tài)壓強(qiáng)隨流道長(zhǎng)度的變化趨勢(shì)，可以看出動(dòng)態(tài)壓強(qiáng)隨著入口速度的增大而增大，尤其是在 3區(qū)更為明顯，然后趨于穩(wěn)定，表明磨粒流對(duì) 4區(qū)處的碰撞與摩擦更加劇烈且可取得表面均勻一致的形貌，可獲得較好的拋光效果。為了直觀準(zhǔn)確分析壁面處動(dòng)態(tài)壓強(qiáng)分布情況，對(duì) 4個(gè)分區(qū)中間位置的上下兩側(cè)壁面取平均值，得到不同入口速度條件下的點(diǎn)膠頭壁面處磨粒流動(dòng)態(tài)壓強(qiáng)分布情況，

所得數(shù)據(jù)匯合如表1所示。

圖4 不同入口速度的磨粒流動(dòng)態(tài)壓強(qiáng)Fig.4 Abrasive flow dynamic pressure at different inlet velocities

表1 不同入口速度下點(diǎn)膠頭壁面處磨粒流動(dòng)態(tài)壓強(qiáng)表Table 1 Abrasive flow pressure on dispensing head wall at different inlet velocities

由表1中的數(shù)據(jù)分析可以看出：1）在同一入口速度條件下，磨粒流動(dòng)態(tài)壓強(qiáng)在各個(gè)區(qū)域的變化趨勢(shì)為：3區(qū)>4區(qū)>2區(qū)>1區(qū)，3區(qū)處的流體動(dòng)態(tài)壓強(qiáng)最大，對(duì)3區(qū)處拋光效果最好，而4區(qū)的動(dòng)態(tài)壓強(qiáng)有所降低，這是由于4區(qū)的內(nèi)徑相同，從磨粒流進(jìn)入4區(qū)到流出4區(qū)，磨粒流與小孔壁面進(jìn)行相互作用，消耗大量能量，磨粒流的能量轉(zhuǎn)化為拋光壁面的切削能，從而對(duì)4區(qū)進(jìn)行光整加工。2）在不同入口速度條件下，隨著入口速度增大，磨粒流動(dòng)態(tài)壓強(qiáng)各區(qū)域都呈增大趨勢(shì)，在 3區(qū)增大較為明顯。以上說(shuō)明流道小的比流道大的拋光效果要好，因此在某些加工情況下，可借助于約束控制模塊減小流道體積，從而增大壓力，提升拋光效率；合理增大入口速度可以提高磨粒流的動(dòng)態(tài)壓強(qiáng)，從而提升點(diǎn)膠頭內(nèi)表面的光整加工效果。

3.3 不同入口速度下的流體與顆粒動(dòng)能耦合場(chǎng)分析

在同樣的初始條件下，對(duì)不同入口速度下的流體動(dòng)能與顆粒動(dòng)能進(jìn)行 CFD-DEM 耦合數(shù)值模擬分析，其動(dòng)能分布及變化趨勢(shì)如圖5和圖6所示。

圖5 不同入口速度的湍流動(dòng)能分布Fig.5 Turbulent kinetic energy distribution at different inlet velocities

從圖 5可以看出，流體湍流動(dòng)能和顆粒動(dòng)能在點(diǎn)膠頭 1區(qū)內(nèi)基本保持不變，此時(shí)的流體湍流動(dòng)能和顆粒動(dòng)能最小，當(dāng)磨粒流到達(dá)點(diǎn)膠頭 3區(qū)時(shí)，流體湍流動(dòng)能和顆粒動(dòng)能開始增大，隨著磨粒流進(jìn)入點(diǎn)膠頭 4區(qū)，流體湍流動(dòng)能和顆粒動(dòng)能繼續(xù)增大，顆粒在3區(qū)和4區(qū)處的動(dòng)能越大，對(duì)點(diǎn)膠頭碰撞越激烈，拋光效果越好。

從圖6可以看出流體動(dòng)能在點(diǎn)膠頭1區(qū)內(nèi)有所降低，是因?yàn)槟チＥc點(diǎn)膠頭進(jìn)行隨機(jī)碰撞，消耗部分能量，磨粒流動(dòng)能轉(zhuǎn)化為切削能，從而對(duì)點(diǎn)膠頭進(jìn)行微磨削；磨粒流到達(dá) 2區(qū)時(shí)，由于點(diǎn)膠頭型腔孔徑開始縮小，磨料速度增大，湍流動(dòng)能開始增大，有利于對(duì) 2區(qū)進(jìn)行去毛刺；磨粒流到達(dá) 4區(qū)時(shí)，湍流動(dòng)能開始減弱，原因與 1區(qū)減弱相同。對(duì) 4個(gè)分區(qū)中間位置的上下兩側(cè)壁面取平均值，得到不同入口速度條件下的點(diǎn)膠頭壁面處磨粒流湍流動(dòng)能分布情況，結(jié)果匯合如表2所示。

圖6 不同入口速度的湍流動(dòng)能Fig.6 Turbulent kinetic energy at different inlet velocities

從表2可以看出：1）在同一入口速度條件下，磨粒流湍流動(dòng)能在各個(gè)區(qū)域的變化趨勢(shì)為：4區(qū)>3區(qū)>2區(qū)>1區(qū)，說(shuō)明點(diǎn)膠頭 1區(qū)壁面磨粒流的湍流動(dòng)能最小，點(diǎn)膠頭4區(qū)壁面處湍流動(dòng)能最大，有利于對(duì)3區(qū)處進(jìn)行去毛刺、倒圓角，對(duì)4區(qū)光整加工效果最好。2）在不同入口速度條件下，隨著入口速度增大，點(diǎn)膠頭壁面處湍流動(dòng)能逐漸增大，磨粒流與壁面發(fā)生的能量交換越劇烈，拋光效果越好，且點(diǎn)膠頭壁面紋理更為均勻。

表2 不同入口速度條件下壁面處磨粒流湍流動(dòng)能分布Table 2 Turbulent kinetic energy distribution of abrasive flow at different inlet velocities

4 不同時(shí)刻的顆粒運(yùn)動(dòng)狀態(tài)及材料去除機(jī)理數(shù)值模擬分析

為了對(duì)顆粒動(dòng)力學(xué)進(jìn)行相關(guān)分析，將顆粒單獨(dú)作為研究對(duì)象，觀測(cè)顆粒在不同時(shí)刻的運(yùn)動(dòng)狀態(tài)，分析顆粒對(duì)壁面的作用效果。由于不同速度下的顆粒對(duì)壁面的作用效果是相同的，這里僅僅是為了體現(xiàn)顆粒對(duì)壁面的材料去除情況，所以選擇其中一種速度分析即可，本文選擇速度45 m/s下的顆粒運(yùn)動(dòng)狀態(tài)進(jìn)行分析，不同時(shí)刻下的顆粒運(yùn)狀態(tài)如圖7所示。

圖7為顆粒從進(jìn)口到出口的5個(gè)瞬時(shí)狀態(tài)。從圖7中可以看出，顆粒在2×10–3s充滿整個(gè)流道，并均勻分散。在顆粒未到達(dá)橫截面變化之前，顆粒速度穩(wěn)定在45 m/s左右；顆粒與壁面發(fā)生碰撞然后反彈，無(wú)數(shù)顆粒循環(huán)往復(fù)的碰撞之后完成對(duì)內(nèi)壁面的加工，由于速度較穩(wěn)定，因而可得到均勻一致的表面形貌；當(dāng)顆粒運(yùn)動(dòng)到橫截面變化處，顆粒顏色由藍(lán)到黃再到紅，顆粒速度明顯增大，顆粒與壁面間的剪切與摩擦隨之增大，最高以110 m/s以上的速度流出，實(shí)現(xiàn)對(duì)零件內(nèi)表面的精加工。

圖7 不同時(shí)刻的顆粒運(yùn)動(dòng)Fig.7 Particles movement at different time

為了闡明單個(gè)顆粒對(duì)微小孔內(nèi)表面的材料去除機(jī)理，本文隨機(jī)選取10個(gè)相同的磨粒，設(shè)置磨粒入射角度為45°，采用Abaqus軟件，觀察磨粒連續(xù)碰撞工件表面的材料去除情況，探討材料去除機(jī)理及規(guī)律，如圖8所示。

從圖8可以看出，由于磨粒入射角度為45°，所以左側(cè)塑性變形要大于右側(cè)的塑性變形。初期磨粒碰撞工件表面，有明顯的變形，說(shuō)明磨粒對(duì)工件表面的作用力超過(guò)了材料的屈服應(yīng)力，工件表面已發(fā)生較為明顯的不可逆的塑性變形，工件表面形成了一個(gè)較為明顯的壓坑。隨著碰撞次數(shù)的增加，工件表面的壓坑不斷增大，當(dāng)磨粒第 6次撞擊工件表面時(shí)，工件表面開始出現(xiàn)變形唇，這說(shuō)明工件表面材料開始被去除，隨著碰撞次數(shù)的累積，材料去除量也越來(lái)越大，工件表面的材料去除是一個(gè)長(zhǎng)時(shí)間的碰撞過(guò)程，即單個(gè)磨粒碰撞工件表面只能對(duì)工件表面產(chǎn)生塑性變形而未出現(xiàn)材料去除，磨粒持續(xù)對(duì)工件表面進(jìn)行撞擊，才能逐步對(duì)工件表面進(jìn)行材料去除。

圖8 碰撞次數(shù)對(duì)材料去除的數(shù)值模擬Fig.8 Numerical simulation of collision times on material removal

為了更直觀的研究磨粒對(duì)工件材料的去除量，利用Abaqus軟件輸出材料磨損量隨碰撞次數(shù)的變化，如圖 9所示。

由圖9可以看出，磨粒碰撞工件前4次時(shí)，工件不發(fā)生材料去除，隨著磨粒碰撞次數(shù)不斷增加，工件材料的去除量也不斷的增加。

為了觀察多磨粒碰撞壁面后的運(yùn)動(dòng)情況，設(shè)置每個(gè)磨粒速度都為45 m/s，選取前4個(gè)磨粒進(jìn)行分析，得到連續(xù)碰撞工件表面的磨粒運(yùn)動(dòng)情況，如圖10所示。

圖9 材料磨去除隨碰撞次數(shù)的變化Fig.9 Change of material wear amount with collisions times

由圖10可以看出，從第1次碰撞到第4次碰撞，可以清楚地看到磨粒碰撞壁面前后的運(yùn)動(dòng)情況。其中第1次碰撞與其他3次碰撞后的磨粒運(yùn)動(dòng)方向有所偏差，這是因?yàn)榈?次碰撞，磨粒到達(dá)工件最深處時(shí)，磨粒速度變?yōu)?，隨后磨粒就會(huì)反彈，但是由于反彈過(guò)程中磨粒還會(huì)撞擊工件附近的待加工表面，迫使磨粒運(yùn)動(dòng)方向有所偏移，而其他3次磨粒會(huì)按照第一次磨粒撞擊后的槽進(jìn)行運(yùn)動(dòng)，經(jīng)過(guò)多次撞擊，工件材料產(chǎn)生斷裂，從而實(shí)現(xiàn)對(duì)工件切削作用。

圖10 顆粒連續(xù)碰撞工件表面的仿真模擬Fig.10 Simulation of particles continuously collision on workpiece surface

5 驗(yàn)證試驗(yàn)

通過(guò)以上數(shù)值分析發(fā)現(xiàn)，磨粒流速度越大，拋光效果越好。為了證實(shí)仿真分析可靠性和磨粒流加工對(duì)點(diǎn)膠頭加工質(zhì)量提升的有效性，有必要進(jìn)行試驗(yàn)驗(yàn)證。根據(jù)所選取的數(shù)值分析條件，利用自行設(shè)計(jì)的磨粒流加工設(shè)備對(duì)點(diǎn)膠頭進(jìn)行磨粒流加工，磨粒流的流體材料為航空煤油，顆粒材料為碳化硅，碳化硅顆粒粒徑為75μm，體積分?jǐn)?shù)為10%，以入口速度30、35、40和45 m/s進(jìn)行磨粒流加工試驗(yàn)。

將磨粒流加工前后的點(diǎn)膠頭內(nèi)表面用光柵表面粗糙度測(cè)量?jī)x進(jìn)行測(cè)量，測(cè)量后的結(jié)果如圖11所示。

從圖11可以看出，隨著磨粒流加工速度增加，磨粒流加工前后點(diǎn)膠頭的表面粗糙度有明顯變化，加工前的表面粗糙度較大，為2.03μm，且表面高低不平，當(dāng)入口速度為30、35、40和45 m/s時(shí)，表面粗糙度值分別降為1.67、1.15、0.96和0.65μm，可以看出，隨著磨粒流入口速度的增大，點(diǎn)膠頭表面粗糙度顯著降低，表面更加平滑，加工后的表面趨于均勻一致，與數(shù)值分析結(jié)果一致。為了進(jìn)一步分析磨粒流速度對(duì)點(diǎn)膠頭表面加工質(zhì)量的影響，利用掃描電鏡對(duì)點(diǎn)膠頭表面的微觀形貌進(jìn)行檢測(cè)，磨粒流加工前后的表面微觀形貌如圖12所示。

由圖12看到，隨著磨粒流入口速度的增大，表面材料去除量明顯增加，零件表面更加光滑和平整，表面質(zhì)量得到顯著提高。表面粗糙度降低的原因是在磨粒流加工過(guò)程中，流體與壁面的毛刺能夠形成無(wú)縫接觸，而顆粒又是以流體為載體，因而顆粒也可以與壁面上的毛刺緊密接觸，隨著磨粒流的連續(xù)流動(dòng)，來(lái)自各個(gè)方向的顆粒與毛刺連續(xù)碰撞摩擦，從而將毛刺剪切掉，達(dá)到光整加工的目的。最終結(jié)果與仿真分析保持一致。

圖11 不同入口速度下加工前后的表面粗糙度Fig.11 Surface roughness before and after processing at different inlet velocities

圖12 不同入口速度下加工前后的微觀形貌對(duì)比Fig.12 Comparison of micro-morphology before and after processing at different inlet velocities

6 結(jié) 論

在磨粒流光整加工過(guò)程中，主要運(yùn)動(dòng)為顆粒與毛刺的碰撞和剪切，流體只是作為一個(gè)載體，引導(dǎo)顆粒的運(yùn)動(dòng)，進(jìn)而對(duì)壁面進(jìn)行拋光，因此研究顆粒在流場(chǎng)的表現(xiàn)行為具有重要的意義。本文運(yùn)用離散元方法，基于CFD-DEM耦合方法，考慮入口速度對(duì)磨粒流加工性能的影響，觀察顆粒在運(yùn)動(dòng)過(guò)程中的力學(xué)行為，對(duì)固液兩相磨粒流進(jìn)行數(shù)值分析，得到以下結(jié)論：

1）數(shù)值分析結(jié)果表明，隨著入口速度的增大，流體與顆粒的湍流作用更加激烈，流體與顆粒的湍流動(dòng)能與入口速度成正比，顆粒與壁面的碰撞更加劇烈，加快了表面材料的去除，提高了磨粒流加工效率。

2）通過(guò)試驗(yàn)研究發(fā)現(xiàn)，點(diǎn)膠頭表面粗糙度隨著入口速度的增大而顯著降低，表面粗糙度Ra值由磨粒流加工前的2.03μm降低到磨粒流加工后的0.65μm，經(jīng)磨粒流加工表面變得光滑平整。

[1]丁金福，劉潤(rùn)之，張克華，等. 磨粒流精密光整加工的微切削機(jī)理[J]. 光學(xué)精密工程，2014，22(12)：3324－3331.Ding Jinfu, Liu Runzhi, Zhang Kehua, et al. Micro cutting mechanism of abrasive flow precision machining[J]. Optics and Precision Engineering, 2014, 22(12): 3324－3331. (in Chinese with English abstract)

[2]李琛，計(jì)時(shí)鳴，譚大鵬，等. 軟性磨粒流加工特性及近壁區(qū)域微切削機(jī)理[J]. 機(jī)械工程學(xué)報(bào)，2014，50(9)：161－168.Li Chen, Ji Shiming, Tan Dapeng, et al. Study of near wall area micro-cutting mechanism and finishing characteristics for softness abrasive flow finishing[J]. Journal of Mechanical Engineering, 2014, 50(9): 161－168. (in Chinese with English abstract)

[3]李俊燁，喬澤民，楊兆軍，等. 介觀尺度下磨料濃度對(duì)磨粒流加工質(zhì)量的影響[J]. 吉林大學(xué)學(xué)報(bào)(工學(xué)版)，2017，47(3)：837－843.Li Junye, Qiao Zemin, Yang Zhaojun, et al. Influence of abrasive concentration on processing quality of abrasive flow in mesoscopic scale[J]. Journal of Jilin University (Engineering and Technology Edition), 2017, 47(3): 837－843. (in Chinese with English abstract)

[4]Gorana V K, Jain V K, Lal G K. Forces prediction during material deformation in abrasive flow machining[J]. Wear,2006, 260(1): 128－139.

[5]Jain V K, Kumar R, Dixit P M, et al. Investigations into abrasive flow finishing of complex workpieces using FEM[J].Wear, 2009, 267(1): 71－80.

[6]李俊燁，衛(wèi)麗麗，張心明，等. 介觀尺度下磨粒流拋光溫度對(duì)發(fā)動(dòng)機(jī)噴油嘴質(zhì)量影響研究[J]. 兵工學(xué)報(bào)，2017,38(10)：2010－2018.Li Junye, Wei Lili, Zhang Xinming, et al. Impact of abrasive flow polishing temperature on nozzle quality under mesoscopic scale[J]. Acta Armamentarii, 2017, 38(10):2010－2018. (in Chinese with English abstract)

[7]李俊燁，董坤，王興華，等. 顆粒微切削表面創(chuàng)成的分子動(dòng)力學(xué)仿真研究[J]. 機(jī)械工程學(xué)報(bào)，2016，52(17)：94－104.Li Junye, Dong Kun, Wang Xinghua, et al. Molecular dynamics simulation research into generative mechanism of particles micro-cutting surface[J]. Journal of Mechanical Engineering, 2016, 52(17): 94－104. (in Chinese with English abstract)

[8]李俊燁，衛(wèi)麗麗，尹延路，等. 磨粒流研拋伺服閥閥芯噴嘴的沖蝕磨損分析[J]. 光學(xué)精密工程，2017，25(07)：1857－1865.Li Junye, Wei Lili, Yin Yanlu, et al. Analysis to erosive wear of abrasive flow polishing servo valve core nozzle[J]. Optics and Precision Engineering, 2017, 25(7): 1857－1865. (in Chinese with English abstract)

[9]計(jì)時(shí)鳴，陳凱，譚大鵬，等. 超聲空化對(duì)軟性磨粒流切削效率和質(zhì)量的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(12)：82－90.Ji Shiming, Chen Kai, Tan Dapeng, et al. Effect of ultrasonic cavitation on maching efficiency and quality of soft abrasive flow[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(12): 82－90. (in Chinese with English abstract)

[10]Uhlmann E, Schmiedel C, Wendler J. CFD Simulation of the Abrasive Flow Machining Process [J]. Procedia Cirp, 2015,31: 209－214.

[11]Butola R, Jain R, Bhangadia P, et al. Optimization to the parameters of abrasive flow machining by Taguchi method[J].Materials Today Proceedings, 2018, 5(2): 4720－4729.

[12]周迪鋒，劉冬玉. 耦合DEM-CFD法雙入口磨粒流動(dòng)力學(xué)模擬及加工試驗(yàn)[J]. 機(jī)電工程，2013，30(12)：1467－1471.Zhou Difeng, Liu Dongyu. Coupling DEM-CFD double inlet abrasive flow dynamics of numerical simulation and test device[J]. Journal of Mechanical & Electrical Engineering,2013, 30(12): 1467－1471. (in Chinese with English abstract)

[13]計(jì)時(shí)鳴，葛江勤，高濤，等. 基于CFD-DEM耦合的面約束軟性磨粒流加工特性研究[J]. 機(jī)械工程學(xué)報(bào)，2018，54(5)：129－141.Ji Shiming, Ge Jiangqin, Gao Tao, et al. Study on machinability of surface-constrained softness abrasive flow based on CFD-DEM coupled method[J]. Journal of Mechanical Engineering, 2018, 54(5): 129－141. (in Chinese with English abstract)

[14]喻黎明，譚弘，鄒小艷，等. 基于CFD-DEM耦合的迷宮流道水沙運(yùn)動(dòng)數(shù)值模擬[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2016，47(8)：65－71.Yu Liming, Tan Hong, Zou Xiaoyan, et al. Numerical simulation of water and sediment flow in labyrinth channel based on coupled CFD-DEM[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(8): 65－71. (in Chinese with English abstract)

[15]喻黎明，鄒小艷，譚弘，等. 基于CFD-DEM耦合的水力旋流器水沙運(yùn)動(dòng)三維數(shù)值模擬[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2016，47(1)：126－132.Yu Liming, Zou Xiaoyan, Tan Hong, et al. 3D Numerical simulation of water and sediment flow in hydrocyclone Based on coupled CFD-DEM[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(1): 126－132.(in Chinese with English abstract)

[16]Kenda J, Duhovnik J, Tav?ar J, et al. Abrasive flow machining applied to plastic gear matrix polishing[J].International Journal of Advanced Manufacturing Technology,2014, 71(1/2/3/4): 141－151.

[17]Sushil M, Vinod K, Harmesh K. Experimental investigation and optimization of process parameters of Al/SiC MMCs finished by abrasive flow machining[J]. Advanced Manufacturing Processes, 2015, 30(7): 902－911.

[18]Sankar M R, Jain V K, Ramkumar J. Nano-finishing of cylindrical hard steel tubes using rotational abrasive flow finishing (R-AFF) process[J]. International Journal of Advanced Manufacturing Technology, 2016, 85(9): 1－9.

[19]陳鐵牛，郭鐘寧，曾柏文，等. 激光誘導(dǎo)空泡微孔拋光機(jī)理及實(shí)驗(yàn)研究[J]. 中國(guó)機(jī)械工程，2018，29(3)：273－278+285.Chen Tieniu, Guo Zhongning, Zeng Baiwen, et al. Study on micro-hole polishing mechanism and experiments based on laser-induced cavitation bubbles[J]. China Mechanical Engineering, 2018, 29(3): 273－278+ 285. (in Chinese with English abstract)

[20]劉薇娜，王尹琛，張英芝，等. 磨料體積分?jǐn)?shù)對(duì)噴油嘴噴孔加工影響分析[J]. 重慶理工大學(xué)學(xué)報(bào)(自然科學(xué))，2018，32(1)：81－85.Liu Weina, Wang Yinchen, Zhang Yingzhi, et al. Abrasive flow machining analysis and research of nozzle[J]. Journal of Chongqing University of Technology(Natural Science), 2018,32(1): 81－85. (in Chinese with English abstract)

[21]李俊燁，胡敬磊，楊兆軍，等. 離散相磨粒粒徑對(duì)磨粒流研拋共軌管質(zhì)量的影響[J]. 吉林大學(xué)學(xué)報(bào)(工學(xué)版)，2018，48(2)：492－499.Li Junye, Hu Jinglei, Yang Zhaojun, et al. Effect of the size of discrete phase abrasive particles on the abrasive flow polishing quality of common rail pipe[J]. Journal of Jilin University (Engineering and Technology Edition), 2018,48(2): 492－499. (in Chinese with English abstract)

[22]高航，李世寵，付有志，等. 金屬增材制造格柵零件磨粒流拋光[J]. 航空學(xué)報(bào)，2017，38(10)：231－239.Gao Hang, Li Shichong, Fu youzhi, et al. Abrasive flow machining of additively manufactured metal grilling parts[J].Acta Aeronautica Et Astronautica Sinica, 2017, 38(10):231－239. (in Chinese with English abstract)

[23]Ji S M, Xiao F Q, Tan D P. A new ultraprecision machining method with softness abrasive flow based on discrete phase model[J]. Advanced Materials Research, 2010, 97－101:3055－3059.

[24]李俊燁，胡敬磊，董坤，等. 固液兩相磨粒流研拋工藝優(yōu)化及質(zhì)量影響[J]. 光學(xué)精密工程，2017，25(06)：1534－1546 Li Junye, Hu Jinglei, Dong Kun, et al. Technological parameter optimization and quality effects on solid-liquid phase abrasive flow polishing[J]. Optics and Precision Engineering, 2017, 25(6): 1534－1546. (in Chinese with English abstract)

[25]唐學(xué)林，余欣，任松長(zhǎng)，等. 固-液兩相流體動(dòng)力學(xué)及其在水力機(jī)械中的應(yīng)用[M]. 鄭州：黃河水利出版社，2006.

[26]Jain R K, Jain V K. Specific energy and temperature determination in abrasive flow machining process[J].International Journal of Machine Tools & Manufacture, 2001,41(12): 1689－1704.

[27]Jain V K, Adsul S G. Experimental investigations into abrasive flow machining (AFM)[J]. International Journal of Machine Tools & Manufacture, 2000, 40(7): 1003－1021.

[28]李晨光，蘇鐵熊，張艷崗. 珩磨工藝研究發(fā)展現(xiàn)狀綜述[J].煤礦機(jī)械，2018，39(3)：1－3.

[29]Kenda J, Pusavec F, Kermouche G, et al. Surface integrity in abrasive flow machining of hardened tool steel AISI D2[J].Procedia Engineering, 2011, 19(1): 172－177.

[30]Gov K, Eyercioglu O, Cakir M V. Hardness effects on abrasive flow machining[J]. Journal of Mechanical Engineering,2013, 59(10): 626－631.

[31]計(jì)時(shí)鳴，李琛，譚大鵬，等. 基于Preston方程的軟性磨粒流加工特性[J]. 機(jī)械工程學(xué)報(bào)， 2011，47(17)：156－163.Ji Shiming, Li Chen, Tan Dapeng, et al. Study on machinability of softness abrasive flow based on preston equation[J].Journal of Mechanical Engineering, 2011, 47(17): 156－163.(in Chinese with English abstract)