多顆磨粒微磨削硬脆材料BK7磨削力仿真研究

周 慧，梁志勇

(東華大學(xué)，上海 201620)

多顆磨粒微磨削硬脆材料BK7磨削力仿真研究

周 慧，梁志勇

(東華大學(xué)，上海 201620)

將金剛石磨粒簡(jiǎn)化為球形磨粒，從單顆微磨削力模型入手，結(jié)合動(dòng)態(tài)磨刃分布模型建立了多顆磨粒微磨削力學(xué)模型；結(jié)合FME和SPH的方法建立磨粒隨機(jī)分布的多顆磨粒微磨削力模型；利用所建立的力學(xué)預(yù)測(cè)模型和仿真模型，對(duì)不同進(jìn)給速度和磨削深度時(shí)所對(duì)應(yīng)的微磨削力進(jìn)行研究，并將模擬結(jié)果與試驗(yàn)進(jìn)行對(duì)比。結(jié)果表明：所建仿真模型可實(shí)現(xiàn)多顆磨粒立式微磨削的過(guò)程，且結(jié)果與試驗(yàn)結(jié)果吻合較好，為后續(xù)研究多顆磨粒微磨削仿真奠定了基礎(chǔ)。

多顆磨粒；微磨削；FME；SPH

隨著現(xiàn)代尖端科技朝著微型化、精密化發(fā)展，高質(zhì)量的微零件和微結(jié)構(gòu)元件的制造對(duì)我國(guó)具有非常重要的戰(zhàn)略意義。硬脆材料因其高硬度、高強(qiáng)度、高耐磨性、高熱強(qiáng)性以及耐腐蝕性和能承受復(fù)雜應(yīng)力等一系列優(yōu)良的機(jī)械物理性質(zhì)，廣泛應(yīng)用于航空航天、光學(xué)等領(lǐng)域。細(xì)微磨削指采用超微顆粒的砂輪在機(jī)械表面上進(jìn)行微納米機(jī)械去除的加工方法，是目前實(shí)現(xiàn)微小精密硬脆材料零件的重要方法。因此研究硬脆材料的微磨削力的產(chǎn)生及影響具有深刻的意義。

許多學(xué)者從單顆磨粒入手建立了磨削力模型。黃輝和林思煌[1]等對(duì)不同形狀的金剛石磨粒進(jìn)行單顆磨粒磨削實(shí)驗(yàn)，得到了磨削力的分布曲線：法向力呈駝峰形，切向力有所波動(dòng)，并且得到最大磨削力的位置；順磨與逆磨對(duì)磨削力的分布曲線沒(méi)有明顯的影響。余劍武和劉智康[2]等建立了單顆CBN磨粒磨削合金鋼20CrMo的力學(xué)和仿真模型，研究了微觀角度的切削成屑機(jī)理，并得到磨削速度、前角變化及切削深度對(duì)磨削力的影響。

上述研究主要集中在單顆磨粒微磨削力建模研究，針對(duì)多顆磨粒微磨削建模研究較少。本文以硬脆材料BK7為研究對(duì)象，在考慮到磨削時(shí)產(chǎn)生裂紋和切屑的過(guò)程是一個(gè)大變形過(guò)程，存在易產(chǎn)生網(wǎng)格畸變導(dǎo)致計(jì)算中止的問(wèn)題，采用了SPH與FEM結(jié)合的算法計(jì)算。本文基于立式微磨削進(jìn)行建模分析。

1 多顆磨粒微磨削力學(xué)模型的建立

本文從單顆磨粒微磨削力學(xué)的模型入手，結(jié)合動(dòng)態(tài)磨刃分布模型來(lái)推導(dǎo)出立式微磨削過(guò)程中多顆磨粒磨削力力學(xué)模型。

1.1 動(dòng)態(tài)磨刃密度模型

切屑形成過(guò)程中，磨具的表面形貌參數(shù)在磨具與工件表面接觸中起到重要作用。在磨削過(guò)程中，控制切屑形成的最相關(guān)的刀具表面參數(shù)是參與磨削的動(dòng)態(tài)磨刃數(shù)。在已知單顆磨粒所產(chǎn)生效應(yīng)的情況下，可以通過(guò)磨刃數(shù)目估算整體磨具的情況。

磨粒隨機(jī)分布在磨具表面上。磨粒的密度、尺寸和形狀均具有隨機(jī)性。由于磨具中參與磨削的磨粒數(shù)目對(duì)磨具所受的整體磨削力產(chǎn)生重要的影響，因此，磨具的表面形貌參數(shù)對(duì)建立一個(gè)可靠的磨削力模型具有重要的意義。參與磨削的磨粒數(shù)目不僅由靜態(tài)磨刃密度所決定，還取決于磨削的運(yùn)動(dòng)參數(shù)。因此，需要建立磨刃密度模型。

動(dòng)態(tài)磨刃數(shù)量是由靜態(tài)磨刃分布情況和磨削的運(yùn)動(dòng)條件及動(dòng)態(tài)磨削效果共同決定的。靜態(tài)磨刃的數(shù)目是由磨粒的尺寸，砂輪的孔隙和修整條件所決定。結(jié)合磨輪法向力的作用，引入修正系數(shù)δ，可得修正靜態(tài)磨刃密度公式為：

Cs(z′)=A(Z′)k

(1)

Z′=Z+δ

(2)

其中，Cs(z)表示修正靜態(tài)磨刃密度(每平方毫米)；A表示與靜態(tài)磨刃數(shù)相關(guān)的經(jīng)驗(yàn)常數(shù)；磨輪修正徑向距離Z′(μm)[3]。

在微磨削過(guò)程中，不是所有的磨粒都參與磨削。因此，Hecker等人認(rèn)為實(shí)際參與磨削過(guò)程的磨粒密度公式為[4]：

Cd(z′)Cs(z′)(1-Vsh/Vt)

(3)

其中，Vt表示砂輪與工件接觸的總體體積；Vsh表示由實(shí)際工作磨粒產(chǎn)生的動(dòng)態(tài)陰影體積。

磨輪工作表面的磨粒是否參與切削主要取決于工件與磨輪表面的接觸面粗糙度。工作磨粒數(shù)隨著接觸面積的增加而增加。此外，微磨削因其在磨削過(guò)程中單位面積的顆粒密度保持不變，從而達(dá)到更均勻去除工件材料的效果，所以磨削效果比常規(guī)磨削更好。

假設(shè)磨具中每個(gè)晶粒都作為一個(gè)磨削刃，晶粒上的磨削邊緣都可以等同于球形邊緣。等效球形磨粒的直徑都可由工件和磨具之間的幾何關(guān)系來(lái)確定。磨粒的接觸圓直徑由下式計(jì)算得到[5]：

(4)

其中，ac表示切削方向上磨粒的最大接觸長(zhǎng)度；bc表示磨刃的切削寬度；dc表示磨粒的接觸元直徑。

當(dāng)磨粒切削深度t已知，磨粒等效直徑dgeq公式為：

(5)

等效磨粒直徑的大小取決于磨粒的形狀和磨粒切削深度[4]。

切屑生成模型如圖1所示。圖中弧長(zhǎng)l′可由下式計(jì)算得到：

(6)

其中t是磨削深度。

圖1 切屑生成模型示意圖Fig.1 Chip formation

微磨削過(guò)程中，當(dāng)工件的寬度等于磨具的直徑D時(shí)，切屑體積可由下式計(jì)算：

(7)

其中Am為切屑的橫截面積，l為切屑的長(zhǎng)度。

因此，在時(shí)間t內(nèi)，切屑的體積Volchip為：

(8)

在同樣的磨削時(shí)間τ內(nèi)，工件的材料去除體積Vwp可由下式計(jì)算：

Vwp=vtBτ

(9)

由于切削點(diǎn)產(chǎn)生的切屑與整體材料去除率相等，得到未變形切屑厚度[6]：

(10)

圖2 (a)立式微磨削示意圖；(b)實(shí)際切屑形狀；(c)理想切屑形狀Fig.2 (a)Schematic diagram of the vertical micro-grinding；(b)Actual chip shape;(c)Idealized chip shape

(11)

(12)

線間距L也可由下式計(jì)算：

(13)

1.2 單顆磨粒切屑生成力模型

Basuray等提到耕犁力與切屑生成力的轉(zhuǎn)變與臨界前角有關(guān)。根據(jù)幾何關(guān)系，臨界前角可由最小未變形切屑厚度和切削刃半徑得到。如下式所示[7]：

(14)

在微磨削中，瞬時(shí)前角通過(guò)幾何關(guān)系計(jì)算得到。如下式所示：

(15)

當(dāng)大于最小切削厚度時(shí)，微切削機(jī)理可以用Merchant模型來(lái)描述。本文考慮未變形切削厚度的有限元模型dh。在磨削時(shí)，刀具前角位置附近，切屑形成力與摩擦角和刀具前角有關(guān)，如圖3所示。因此，在二維簡(jiǎn)化模型中，微切削力增量dfcgx可由下式求得[8，9]：

(16)

其中，dfcgx為材料剪切強(qiáng)度；σs為材料的剪切強(qiáng)度；β為摩擦角；φ為剪切角。

圖3 微磨削過(guò)程中單顆磨粒作用機(jī)理Fig.3 Mechanism of single abrasive particles in micro-grinding process

為了建立單顆磨粒接觸模型，將單顆磨粒形狀假設(shè)為球形。理想的球形磨粒的橫截面與切削深度增長(zhǎng)率及切削前角的關(guān)系如下[10]：

(17)

因此，每單顆磨粒的切屑形成力如下：

(18)

單顆磨粒的切屑形成力法向分力如下：

(19)

1.3 單顆磨粒的耕犁力模型

在微磨削過(guò)程中，當(dāng)切削深度小于產(chǎn)生切屑的最小切削深度時(shí)，材料只發(fā)生塑性變形而不產(chǎn)生切屑，這個(gè)力稱為耕犁力。耕犁力在單顆磨粒相互作用的情況下的變化可通過(guò)布氏硬度測(cè)試。HB定義為負(fù)載fbrnl與壓痕彎曲面積的比值[11]：

(20)

其中，d為直徑；b為壓痕直徑。

由于耕犁力的影響使得壓痕力朝著αcr方向發(fā)展。單顆磨粒的總切向力和法向耕犁力結(jié)合壓痕力和摩擦力的影響來(lái)計(jì)算。如下式所示：

fpgy=fbml(sinα+μpcosα)

(21)

fpgx=fbml(cosα-μpsinα)

(22)

式中μp為耕犁摩擦系數(shù)。

耕犁摩擦系數(shù)μp可由Goddard，Wilman和Sin等的預(yù)測(cè)模型求出。如下式所示[12，13]：

(23)

其中，w為壓痕寬度；dgep為球形磨粒半徑。

1.4 磨削力理論模型

單顆磨粒的磨削力由切屑生成力和耕犁力組成。計(jì)算公式如下：

fn=fcgy+fpgy

ft=fcgx+fpgx

(24)

本模型中，剪切角的大小由Emst和Merchant模型來(lái)預(yù)測(cè)。計(jì)算公式如下[14]：

(25)

最終所要求的磨削力Fres可以由三個(gè)部分組成：法向力Fn、徑向力Fl和切向力Ft，如圖4所示。這些力可認(rèn)為是磨削過(guò)程中，單顆磨粒所受的力的合力。計(jì)算公式如下：

(26)

(27)

(28)

其中，Ng為動(dòng)態(tài)切削刃數(shù)；Fn為單顆磨粒所受法向力；Ft為單顆磨粒所受切向力。

圖4 微磨削力示意圖Fig.4 Micro-grinding force

2 多顆磨粒的微磨削力仿真模型

2.1 多顆磨粒幾何及分布模型的建立

微磨削過(guò)程從本質(zhì)上來(lái)說(shuō)是一種瞬時(shí)動(dòng)態(tài)的過(guò)程，伴隨著大變形、瞬態(tài)沖擊等，ABAQUS作為一種常用的有限元分析軟件，能很好的求解瞬態(tài)動(dòng)力學(xué)問(wèn)題，適用于微磨削。

目前的研究中，主要將磨粒簡(jiǎn)化為規(guī)則的形狀，例如：球形、多棱錐形、球形等，本文采用球形模型。

磨粒一般隨機(jī)分布在磨具表面上，本文在簡(jiǎn)化磨輪模型時(shí)采用隨機(jī)函數(shù)生成隨機(jī)位置，再將磨粒分布在隨機(jī)位置上。最終建立立式磨削模型，如圖5所示。

圖5 立式微磨削建模示意圖Fig.5 Finite element modeling for vertical micro- grinding

2.2 多顆磨粒微磨削模型的建立

金剛石磨粒因其較工件材料BK7，具有硬度高、強(qiáng)度大的特點(diǎn)，在磨削過(guò)程中不會(huì)發(fā)生大的形變，遂假設(shè)其為剛體。磨粒與工件的材料力學(xué)性能如表1所示。

表1 磨粒與工件的力學(xué)性能Tab.1 Mechanical properties of abrasive grains and workpieces

考慮到微磨削形成切屑時(shí)變形較大，采用有限元法容易產(chǎn)生網(wǎng)格畸變，導(dǎo)致計(jì)算終止。SPH采用純拉格朗日方法，避免了網(wǎng)格畸變問(wèn)題，適用于求解動(dòng)態(tài)大變形問(wèn)題，但由于其計(jì)算效率較低，且邊界處理不如有限元法成熟，所以采用工件材料失效時(shí)轉(zhuǎn)化為SPH的方法來(lái)建模。劃分網(wǎng)格時(shí)，將工件與磨具作用區(qū)域的網(wǎng)格進(jìn)行細(xì)分，保證仿真精度，遠(yuǎn)離磨削作用的區(qū)域網(wǎng)格進(jìn)行粗化。在仿真模型中設(shè)置相應(yīng)的參數(shù)為：磨具直徑為670 μm、進(jìn)給速度v(1～25 μm/mint)、磨削深度t(0.1～0.5 μm)，對(duì)硬脆材料BK7微磨削過(guò)程進(jìn)行仿真。仿真模型如圖6所示。

圖6 多顆磨粒微磨削的Abaqus仿真模型Fig.6 Abaqus simulation model of multi-abrasive micro-grinding

3 多顆磨粒微磨削的磨削力分析驗(yàn)證

將ABAQUS模型中計(jì)算所得的x、y、z三個(gè)方向的磨削力輸出導(dǎo)出，利用MATLAB進(jìn)行數(shù)據(jù)擬合，得到各工藝下的磨削基本穩(wěn)定的磨削力的大小，與力學(xué)模型MATLAB值相對(duì)比。并參照文獻(xiàn)[15]中實(shí)驗(yàn)所得數(shù)據(jù)，得到圖7和圖8。

圖7 切削深度對(duì)微磨削力的影響Fig.7 Effect of depth of cut on the micro-grinding force

圖8 進(jìn)給速度對(duì)微磨削力的影響Fig.8 Effect of feed rate on the micro-grinding force

由圖7和圖8可看出，力學(xué)預(yù)測(cè)模型和仿真模型在不同切削深度和不同的進(jìn)給速度下和實(shí)驗(yàn)數(shù)據(jù)具有相同的規(guī)律且微磨削力值相仿?？勺C明仿真模型的正確性。

4 結(jié)論

第一，從單顆磨粒微磨削力學(xué)的力學(xué)模型入手，結(jié)合動(dòng)態(tài)磨刃分布模型來(lái)推導(dǎo)出立式微磨削過(guò)程中多顆磨粒磨削力力學(xué)模型。

第二，利用ABAQUS軟件，將FME和SPH結(jié)合起來(lái)建立了磨粒隨機(jī)分布的多顆磨粒微磨削模型。

第三，力學(xué)模型預(yù)測(cè)結(jié)果及ABAQUS仿真模擬結(jié)果與實(shí)驗(yàn)結(jié)果吻合良好，表明提出的建模方法是合理可行的，為研究多顆磨粒微磨削仿真奠定了基礎(chǔ)。

[1] 黃輝，林思煌，徐西鵬. 單顆金剛石磨粒磨削玻璃的磨削力研究[J]. 中國(guó)機(jī)械工程，2010，21(11)：1278-1282.

[2] 余劍武，劉智康，吳耀，等. 合金鋼20CrMo的單顆磨粒高速磨削仿真研究[J]. 制造技術(shù)與機(jī)床，2015，(12)：97-102.

[3] Shaw MC.Principles of abrasive processing[M]. Clarendon, Oxford，1996.

[4] Hecker RL，Ramoneda IM，Liang SY.Analysis of wheel topography and grit force for grinding process modeling[J].J Manuf Proc，2003，(05)：13-23.

[5] Verkerk J, Peters J .Final report concerning CIRP cooperative work on the characterization of grinding wheel topography[J].Ann CIRP，1977，(26)：385-395.

[6] Srihari G，Lal GK.Mechanics of vertical surface grinding[J].J Mat Proces Tech，1994，(44)：14-28.

[7] Basuray PK，Misra BK，Lal GK.Transition from ploughing to cutting during machining with blunt tools[J].Wear，1977，(43)：341-349.

[8] Liu K，Li XP.Modeling of ductile cutting of tungsten carbide[J]. Trans NAMRI/SME，2001，(29)：251-258.

[9] Son SM，LimHS，Ahn JH.Effects of the friction coefficient on the minimum cutting thickness in micro cutting[J]. Int J Mach Tools Manuf，2005， 45(4-5)：529-535.

[10] Park HW, Liang SY.Force modeling of microscale grinding process incorporating thermal effects[J]. Int J Adv Manuf Techn，2009，44(5-6)：476-486.

[11] Shaw, MC. Fundamentals of grinding. In: Proc of the int’1grinding conf. on new developments in grinding[J].Pittsburgh，1972，PA,April,：221-258.

[12] M Goddard J,WilmanH.A theory of friction and wear during theabrasion of metals[J]. Wear，1962，(05)：114-135.

[13] Sin H, Saka N, Suh NP .Abrasive wear mechanisms and the grit size effect[J]. Wear，1979，(55)：163-190.

[14] MerchantME.Mechanics of metal cutting process, and type 2 chip[J]. J Appl Phys，1945，(16)：318-324.

[15] Asma P，M.rahman, Y. S. Wong.Modeling and simulation of cutting forces generated during vertical micro-grinding[J]. Int J Adv Manuf Techn ，2014，(71)：1539-1548.

Simulationresearchonmicro-grindingforceofBK7basedonmultiabrasivegrainscutting

ZHOU Hui, LIANG Zhi-yong

(Donghua University, Shanghai 201620, China)

The typical diamond abrasive grain is simplified into sphere. Micro-grinding force model of multi abrasive grains is set up based on combine the grinding force of single abrasive grain and the model of dynamic cutting edges. A multi abrasive micro grinding force simulation model based on FME and SPH was established. Based on the established mechanical prediction model and simulation model, the micro grinding force corresponding to different feed rate and grinding depth is studied, and the simulation results are compared with the experimental results. The results show that the simulation model by built can realize the vertical micro grinding with multiple abrasive grains process; and the whole process is in good agreement with the experimental results. It can lay the foundation for the simulation of multi abrasive grains micro grinding.

Multi abrasive grains; Micro-grinding; FEM; SPH

TG580.614

A

1674-8646(2017)21-0001-05

2017-10-20

周慧(1993-)，女，研究生。

梁志勇(1966-)，男，工學(xué)博士，副教授，e-mail：zhyliang@dhu.edu.cn。