肖行志，鄭 侃，廖文和

（南京理工大學(xué)機(jī)械工程學(xué)院，南京 210094）

超聲振動(dòng)輔助磨削牙科氧化鋯陶瓷切削力預(yù)測(cè)模型研究

肖行志，鄭 侃，廖文和

（南京理工大學(xué)機(jī)械工程學(xué)院，南京 210094）

通過(guò)單因素實(shí)驗(yàn)設(shè)計(jì)，開(kāi)展了牙科氧化鋯陶瓷的超聲振動(dòng)輔助磨削實(shí)驗(yàn)，分別建立了超聲振動(dòng)輔助磨削加工切削力指數(shù)預(yù)測(cè)模型、BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型以及理論預(yù)測(cè)模型。通過(guò)驗(yàn)證實(shí)驗(yàn)對(duì)比分析了三種模型的預(yù)測(cè)精度，并闡明了誤差產(chǎn)生的原因。結(jié)果表明，基于BP神經(jīng)網(wǎng)絡(luò)的切削力預(yù)測(cè)模型相對(duì)于指數(shù)和理論模型具有較高的預(yù)測(cè)精度，其平均相對(duì)誤差僅為9．60%，理論模型因未能考慮材料的塑性流動(dòng)去除，導(dǎo)致預(yù)測(cè)精度較低。

超聲振動(dòng)輔助磨削；牙科氧化鋯陶瓷；切削力預(yù)測(cè)；指數(shù)模型；BP神經(jīng)網(wǎng)絡(luò)；理論模型

氧化鋯陶瓷因其良好的生物相容性、耐磨性和美學(xué)特性而備受口腔修復(fù)領(lǐng)域的青睞［1－4］，然而由于其高硬度、低斷裂韌性，導(dǎo)致其難以直接加工。目前，全鋯牙冠是通過(guò)先對(duì)預(yù)燒結(jié)氧化鋯陶瓷進(jìn)行切削加工，再進(jìn)行二次燒結(jié)的方法來(lái)獲得。但此制作工藝復(fù)雜，燒結(jié)過(guò)程中牙冠各部位的收縮變形難以控制［5］。因此，將超聲振動(dòng)輔助磨削技術(shù)（原理如圖1）引入口腔修復(fù)領(lǐng)域，實(shí)現(xiàn)完全燒結(jié)氧化鋯陶瓷的直接加工，從而可提高全鋯牙冠的加工精度，同時(shí)縮短病人的等待時(shí)間。

圖1 超聲振動(dòng)輔助磨削原理圖Fig．1 Schematic illustration of UVAG

切削力直接影響了切削加工過(guò)程中的穩(wěn)定性以及工件的表面質(zhì)量，特別是對(duì)于表面形貌復(fù)雜、壁厚較薄的全鋯牙冠而言，影響更加顯著。因此，開(kāi)展氧化鋯陶瓷超聲振動(dòng)輔助磨削過(guò)程中的切削力研究顯得尤為必要。目前，國(guó)內(nèi)外學(xué)者分別從實(shí)驗(yàn)分析和理論建模角度對(duì)超聲振動(dòng)輔助加工陶瓷等硬脆性材料的切削力開(kāi)展了大量的研究。在實(shí)驗(yàn)分析方面，Lauwers等［6］通過(guò)對(duì)比分析不同振幅下的超聲振動(dòng)輔助磨削氧化鋯陶瓷的切削力大小，指出超聲振動(dòng)的附加有助于減小切削力，且振幅越大效果越明顯。呂東喜等［7］通過(guò)光滑粒子流體動(dòng)力學(xué)方法對(duì)超聲振動(dòng)輔助磨削加工過(guò)程中，沖擊裂紋的形成機(jī)理及對(duì)加工過(guò)程的影響進(jìn)行數(shù)值模擬，研究表明沖擊裂紋及亞表層裂紋的產(chǎn)生有利于減小磨粒的切削力。文獻(xiàn)［8］研究了主軸轉(zhuǎn)速對(duì)超聲振動(dòng)輔助磨削牙科氧化鋯切削力的影響，指出隨著主軸轉(zhuǎn)速增大切削力逐漸減小。在理論建模方面，Pei等［9］通過(guò)對(duì)脆性材料超聲振動(dòng)輔助加工過(guò)程中裂紋的產(chǎn)生與擴(kuò)展進(jìn)行分析，提出了基于壓痕斷裂力學(xué)理論的超聲振動(dòng)輔助鉆削切削力理論模型，并通過(guò)實(shí)驗(yàn)進(jìn)行了驗(yàn)證。Zhang等［10］通過(guò)對(duì)單顆磨粒的切削過(guò)程進(jìn)行分析，提出了超聲振動(dòng)輔助磨削脆性材料的切削力理論模型。但到目前為止，將數(shù)值模型（指數(shù)模型、神經(jīng)網(wǎng)絡(luò)模型等）用于預(yù)測(cè)超聲振動(dòng)輔助磨削過(guò)程中切削力的研究卻鮮有報(bào)道。作為常用的切削力預(yù)測(cè)模型之一，數(shù)值模型在傳統(tǒng)切削加工研究中發(fā)揮著重要作用［11－12］。因此，有必要將數(shù)值模型引入到超聲振動(dòng)輔助加工領(lǐng)域，以便采用多種方法對(duì)加工過(guò)程中的切削力進(jìn)行預(yù)測(cè)研究。

本文將通過(guò)統(tǒng)計(jì)學(xué)習(xí)法，分別建立超聲振動(dòng)輔助磨削牙科氧化鋯陶瓷的切削力指數(shù)預(yù)測(cè)模型和BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型。并將此兩種模型與理論預(yù)測(cè)模型進(jìn)行對(duì)比分析，研究各模型的預(yù)測(cè)精度及預(yù)測(cè)誤差產(chǎn)生的原因。

1 超聲振動(dòng)輔助磨削實(shí)驗(yàn)

本次實(shí)驗(yàn)采用的機(jī)床為德國(guó)DMG Ultrasonic linear 20切削加工中心，其中刀具外徑為8 mm，內(nèi)徑為6．4 mm，切削液為乳化液。超聲振動(dòng)頻率為23．9 kHz，振幅為5μm。工件材料為完全燒結(jié)的牙科氧化鋯陶瓷，由秦皇島愛(ài)迪特高技術(shù)陶瓷有限公司提供，其主要力學(xué)性能見(jiàn)表1。工件尺寸為30 mm×15 mm×5 mm。旋轉(zhuǎn)超聲加工實(shí)驗(yàn)參數(shù)見(jiàn)表2。

表1 完全燒結(jié)氧化鋯陶瓷的力學(xué)性能Tab.1 Mechanicalproperty of full sintered zirconia ceram ics

表2 實(shí)驗(yàn)參數(shù)設(shè)置Tab.2 Design of experiment

以表2為依據(jù)設(shè)置切削參數(shù)，共進(jìn)行12組實(shí)驗(yàn)。在實(shí)驗(yàn)過(guò)程中，軸向切削力通過(guò)Kistler9257測(cè)力儀測(cè)得，具體實(shí)驗(yàn)數(shù)據(jù)見(jiàn)表3。

表3 軸向切削力實(shí)驗(yàn)結(jié)果Tab.3 Experimental results of axial cutting force

2 預(yù)測(cè)模型的建立

根據(jù)超聲振動(dòng)輔助磨削氧化鋯陶瓷實(shí)驗(yàn)中獲得的切削力數(shù)據(jù)，分別建立切削力指數(shù)預(yù)測(cè)模型和BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型。并通過(guò)對(duì)超聲振動(dòng)輔助磨削脆性材料的去除機(jī)理進(jìn)行分析，得出切削力理論預(yù)測(cè)模型。

2.1 指數(shù)預(yù)測(cè)模型

超聲振動(dòng)輔助磨削實(shí)質(zhì)上是普通磨削加工和超聲振動(dòng)的復(fù)合。但在本實(shí)驗(yàn)中振動(dòng)參數(shù)為定值，因此，可將指數(shù)預(yù)測(cè)模型表示為：

式中：F為切削力；C為材料性能、刀具以及機(jī)床本身等對(duì)切削力的影響；n為主軸轉(zhuǎn)速；ap為切削深度；vf為進(jìn)給速度；a1、a2、a3為待定系數(shù)，其大小代表各因素對(duì)切削力的影響程度。

根據(jù)式（2），反求出C＝103．8442＝6 986．23。將各參數(shù)代入式（1）中，得到最終的超聲振動(dòng)輔助磨削切削力指數(shù)預(yù)測(cè)模型為：

2.2 BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型

人工神經(jīng)網(wǎng)絡(luò)（Artificial Neural Networks，ANN）是由大量類似于生物神經(jīng)元的處理單元互相連接而構(gòu)成的非線性復(fù)雜網(wǎng)絡(luò)系統(tǒng)。它用某種相對(duì)簡(jiǎn)單的數(shù)學(xué)模型來(lái)描述生物神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)，并讓其在一定算法指導(dǎo)下在某種程度上能夠模擬生物神經(jīng)網(wǎng)絡(luò)的智能行為，來(lái)解決傳統(tǒng)算法所不能完成的智能信息處理問(wèn)題［13］。誤差逆向傳播算法的多層前饋神經(jīng)網(wǎng)絡(luò)，即BP神經(jīng)網(wǎng)絡(luò)，是目前應(yīng)用最多的一種人工神經(jīng)網(wǎng)絡(luò)模型［14］。

典型的BP神經(jīng)網(wǎng)絡(luò)通常由輸入層、隱含層和輸出層構(gòu)成。每層的神經(jīng)元個(gè)數(shù)不確定，根據(jù)訓(xùn)練需要自行設(shè)置。同層的神經(jīng)網(wǎng)絡(luò)之間互不干擾，相鄰層間通過(guò)權(quán)值連接。輸入的訓(xùn)練樣本由輸入層到達(dá)隱含層，通過(guò)一定的算法對(duì)樣本進(jìn)行訓(xùn)練，最終通過(guò)輸出層輸出訓(xùn)練結(jié)果。根據(jù)已開(kāi)展的超聲振動(dòng)輔助磨削氧化鋯陶瓷實(shí)驗(yàn)數(shù)據(jù)，選定輸入層神經(jīng)元數(shù)目為3，隱含層層數(shù)為1，神經(jīng)元數(shù)目為5，輸出神經(jīng)元數(shù)目為1，最終的BP神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)見(jiàn)圖2。

BP神經(jīng)網(wǎng)絡(luò)的基本訓(xùn)練流程見(jiàn)圖3。在將訓(xùn)練樣本即實(shí)驗(yàn)數(shù)據(jù)輸入到BP神經(jīng)網(wǎng)絡(luò)之前，需要對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行歸一化處理，將各個(gè)參數(shù)的數(shù)據(jù)均歸一化到區(qū)間（0，1）。本文通過(guò)matlab編寫B(tài)P神經(jīng)網(wǎng)絡(luò)訓(xùn)練程序，采用動(dòng)量梯度下降法對(duì)所建立的神經(jīng)網(wǎng)絡(luò)模型進(jìn)行訓(xùn)練。其中，訓(xùn)練目標(biāo)最小誤差為0．001 5，允許的最大訓(xùn)練步長(zhǎng)為50 000步。通過(guò)該訓(xùn)練模型對(duì)12組實(shí)驗(yàn)數(shù)據(jù)進(jìn)行訓(xùn)練，當(dāng)訓(xùn)練誤差小于訓(xùn)練目標(biāo)最小誤差時(shí)，訓(xùn)練結(jié)束，形成超聲振動(dòng)輔助磨削氧化鋯陶瓷切削力的BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型。

圖2 BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)圖Fig．2 Structure diagram of BP ANN

圖3 BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練基本流程圖Fig．3 Flow diagram of BP ANN

2.3 理論預(yù)測(cè)模型

目前，主要是基于壓痕斷裂力學(xué)理論對(duì)超聲振動(dòng)輔助磨削脆性材料的去除過(guò)程進(jìn)行理論分析，材料的去除被看作為大量金剛石壓頭對(duì)材料進(jìn)行沖壓作用的結(jié)果［15］。在這個(gè)過(guò)程中，材料通過(guò)橫向裂紋的產(chǎn)生及最終擴(kuò)展到工件表面而形成脆性斷裂來(lái)實(shí)現(xiàn)去除，其中裂紋的形成見(jiàn)圖4。

圖4 單顆模型形成的裂紋示意圖Fig．4 Schematic illustration of cracks of the single abrasive particle

金剛石磨粒隨著刀具做軸向的超聲頻振動(dòng)，從而對(duì)材料產(chǎn)生沖壓作用，在磨粒的正下方，最開(kāi)始產(chǎn)生塑性變形區(qū)域，當(dāng)壓力逐步增大時(shí)，從塑性變形區(qū)開(kāi)始形成中央裂紋、徑向裂紋和橫向裂紋。中央裂紋自壓痕底部向下延伸，導(dǎo)致材料的強(qiáng)度下降。而橫向裂紋自壓痕底部塑性變形區(qū)交界處向上及兩側(cè)擴(kuò)展，當(dāng)裂紋擴(kuò)展到工件表面時(shí)，最終形成材料的脆性斷裂去除［16］。其中，橫向裂紋的寬度CL和深度Ch可分別通過(guò)式（5）和式（6）求得。

式中：C1為無(wú)量綱常量；a為金剛石磨粒的頂角；E為工件材料的彈性模量；HV為材料的維氏硬度；v為材料的泊松比；Fn為金剛石磨粒受到的正壓力。

由式（5）和式（6）可知，橫向裂紋的寬度和深度與金剛石磨粒的正壓力有著直接關(guān)系。當(dāng)綜合考慮參與磨削加工的全部磨粒時(shí)，可以求出切削力與材料性能、刀具參數(shù)、切削參數(shù)以及振動(dòng)參數(shù)的關(guān)系?；趬汉蹟嗔蚜W(xué)理論建立了超聲振動(dòng)輔助磨削脆性材料的切削力理論模型，具體形式見(jiàn)式（7）。

式中：F為旋轉(zhuǎn)超聲加工脆性材料時(shí)的軸向切削力；K為比例系數(shù)，可通過(guò)實(shí)驗(yàn)獲得；Ca為金剛石刀具的磨粒濃度；wc為切削寬度；Do和Di分別為刀具的外徑和內(nèi)徑；x與金剛石磨粒的幾何形狀有關(guān)，此處取x＝1．85；Sa為磨粒的尺寸。

根據(jù)式（7）可知，比例系數(shù)K是未知量，其他參數(shù)取決于加工過(guò)程中的工件材料、刀具以及所選擇的切削參數(shù)和振動(dòng)參數(shù)。因此，當(dāng)求得比例系數(shù)K后，即可得到符合本實(shí)驗(yàn)條件的切削力理論預(yù)測(cè)模型。

比例系數(shù)K由式（7）反求可得：

在本實(shí)驗(yàn)中，金剛石磨粒的頂角a＝π／2，磨粒的濃度Ca＝100，磨粒的尺寸Sa＝126 mm。將此刀具參數(shù)、表3中的加工數(shù)據(jù)、表1中工件材料的力學(xué)性能指標(biāo)以及振動(dòng)參數(shù)代入式（8）中，求出每組實(shí)驗(yàn)所對(duì)應(yīng)的K值（見(jiàn)表4）。

表4 每組實(shí)驗(yàn)所對(duì)應(yīng)的K值Tab.4 The value ofK of each experim ent

從表4可知，個(gè)別比例系數(shù)K的值過(guò)大或過(guò)小。需要通過(guò)去掉最大值和最小值再求平均值的方法來(lái)減少誤差，而后得到切削力理論模型中的比例系數(shù)K＝52 474．76。因此，超聲振動(dòng)輔助磨削牙科氧化鋯陶瓷切削力理論預(yù)測(cè)模型為：

3 預(yù)測(cè)結(jié)果與討論

基于上述實(shí)驗(yàn)數(shù)據(jù)，建立了超聲振動(dòng)輔助磨削牙科氧化鋯陶瓷的軸向切削力指數(shù)預(yù)測(cè)模型、BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型和切削力理論預(yù)測(cè)模型。為了分析各模型的預(yù)測(cè)精度，在同等實(shí)驗(yàn)條件下設(shè)計(jì)了3組驗(yàn)證實(shí)驗(yàn)，實(shí)驗(yàn)過(guò)程中切削參數(shù)及軸向切削力實(shí)測(cè)結(jié)果見(jiàn)表5。

表5 驗(yàn)證實(shí)驗(yàn)切削參數(shù)及結(jié)果Tab.5 Cutting parameters and experimentalresults of the verification experim ent

將驗(yàn)證實(shí)驗(yàn)中的切削參數(shù)分別代入切削力指數(shù)預(yù)測(cè)模型、BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型和理論預(yù)測(cè)模型中，得出不同切削參數(shù)下的軸向切削力預(yù)測(cè)值見(jiàn)表6。

表6 三種切削力預(yù)測(cè)模型的預(yù)測(cè)值及相對(duì)誤差Tab.6 Predictive value and relative error of three prediction models

從表6可知，BP神經(jīng)網(wǎng)絡(luò)模型的預(yù)測(cè)精度最高，其平均相對(duì)誤差僅為9．60%，其次為指數(shù)預(yù)測(cè)模型，理論模型預(yù)測(cè)精度最低，其平均相對(duì)誤差達(dá)到48．12%。這是因?yàn)?，理論模型的建立是以材料的脆性斷裂去除為前提，認(rèn)為在旋轉(zhuǎn)超聲加工脆性材料時(shí)，材料的去除完全是通過(guò)裂紋的產(chǎn)生與擴(kuò)展實(shí)現(xiàn)的。但目前有些文獻(xiàn)［17－18］指出，在脆性材料的超聲振動(dòng)輔助磨削過(guò)程中，同樣存在塑性流動(dòng)去除。材料的塑脆性轉(zhuǎn)變主要取決于切削深度，當(dāng)切削深度大于臨界切深時(shí)，材料去除方式由塑性流動(dòng)轉(zhuǎn)變?yōu)榇嘈詳嗔选Ｔ趯?shí)際加工過(guò)程中，由于軸向超聲振動(dòng)的作用，磨粒切入切出工件，導(dǎo)致實(shí)際的切削深度處于周期性的變化中，從而使得在材料去除過(guò)程中發(fā)生脆－塑性去除轉(zhuǎn)變。其次，理論模型在建模過(guò)程中，沒(méi)有考慮到磨粒無(wú)規(guī)則排布帶來(lái)的相鄰磨粒的切削軌跡的干涉問(wèn)題，致使在建模過(guò)程中對(duì)理論材料去除體積的計(jì)算不準(zhǔn)確，從而影響模型預(yù)測(cè)精度。此外，理論模型中忽略的刀具磨損、切削溫度等也對(duì)模型的精度有一定的影響。因此，完全基于壓痕斷裂力學(xué)建立的切削力理論模型必然會(huì)與實(shí)際切削過(guò)程存在偏差。圖5為三種切削力模型預(yù)測(cè)值與驗(yàn)證實(shí)驗(yàn)值的對(duì)比分析曲線，從圖5可更直觀的看出BP神經(jīng)網(wǎng)絡(luò)模型的預(yù)測(cè)值與實(shí)驗(yàn)值的擬合效果最好。

圖5 模型預(yù)測(cè)值與實(shí)驗(yàn)值的對(duì)比曲線Fig．5 Comparison ofpredictive value and experimental value

圖6 模型相對(duì)誤差與主軸轉(zhuǎn)速關(guān)系的對(duì)比曲線Fig．6 Relationship between relative error and spindle speed

圖6為各模型預(yù)測(cè)值的相對(duì)誤差隨主軸轉(zhuǎn)速的變化曲線。從圖6可知，當(dāng)主軸轉(zhuǎn)速由2 000 r／min增加8 000 r／min時(shí)，指數(shù)預(yù)測(cè)模型和BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型的相對(duì)誤差整體上表現(xiàn)出相同的變化趨勢(shì)，先逐漸減小而后緩慢增大。相對(duì)誤差的最大值均出現(xiàn)主軸轉(zhuǎn)速為2 000 r／min時(shí)，其值分別為19．71%和16．97%，且在主軸轉(zhuǎn)速為5 000 r／min時(shí)，該兩種模型的相對(duì)誤差均達(dá)到最小，分別為4．06%和4．92%。但是理論預(yù)測(cè)模型的相對(duì)誤差隨著主軸轉(zhuǎn)速的增大而持續(xù)減小，其中，在主軸轉(zhuǎn)速處于2 000～5 000 r／min時(shí)，相對(duì)誤差由117．53%迅速減小為18．24%，而當(dāng)主軸轉(zhuǎn)速＞5 000 r／min后，減小較為緩慢，直到轉(zhuǎn)速為8 000 r／min時(shí)，模型的相對(duì)誤差達(dá)到最小，其值為8．59%。在整個(gè)變化過(guò)程中，指數(shù)預(yù)測(cè)模型、BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型和理論預(yù)測(cè)模型的相對(duì)誤差最大值與最小值之差分別為15．56%、12．05%以及108．94%。以上分析表明，BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型的穩(wěn)定性最高，指數(shù)預(yù)測(cè)模型的穩(wěn)定性略低于BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型，且該兩種模型均在主軸轉(zhuǎn)速為5 000 r／min時(shí)，模型的預(yù)測(cè)精度最高。理論模型預(yù)測(cè)值的相對(duì)誤差波動(dòng)較大，只有當(dāng)主軸轉(zhuǎn)速達(dá)到8 000 r／min時(shí)，模型的預(yù)測(cè)精度才屬于可接受范圍。

4 結(jié) 論

將超聲振動(dòng)輔助磨削技術(shù)應(yīng)用于口腔修復(fù)領(lǐng)域，實(shí)現(xiàn)全鋯牙冠的直接加工成型。在加工過(guò)程中，為了減小切削力對(duì)全鋯牙冠質(zhì)量的影響，開(kāi)展了超聲振動(dòng)輔助磨削牙科氧化鋯陶瓷切削力預(yù)測(cè)模型的研究，分別建立的切削力指數(shù)預(yù)測(cè)模型、BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型和理論預(yù)測(cè)模型，通過(guò)將各模型的預(yù)測(cè)值與驗(yàn)證實(shí)驗(yàn)值進(jìn)行對(duì)比分析，得出以下結(jié)論：

（1）基于BP神經(jīng)網(wǎng)絡(luò)的切削力模型整體預(yù)測(cè)精度最高，其預(yù)測(cè)值與實(shí)驗(yàn)值的平均相對(duì)誤差僅為9．60%。指數(shù)預(yù)測(cè)模型的精度略低于BP神經(jīng)網(wǎng)絡(luò)模型，其平均相對(duì)誤差為11．09%。理論模型因未能考慮材料的塑性流動(dòng)去除，導(dǎo)致其預(yù)測(cè)精度較低，其平均相對(duì)誤差達(dá)到48．12%。

（2）BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型的穩(wěn)定性最高，指數(shù)預(yù)測(cè)模型略低于BP神經(jīng)網(wǎng)絡(luò)模型，理論模型預(yù)測(cè)值的相對(duì)誤差波動(dòng)較大。BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型和指數(shù)預(yù)測(cè)模型的相對(duì)誤差隨著轉(zhuǎn)速的增大先減小后增大，均在轉(zhuǎn)速為5 000 r／min時(shí)達(dá)到最小，分別為4．06%和4．92%。理論預(yù)測(cè)模型的相對(duì)誤差隨轉(zhuǎn)速的增加而迅速減小，在8 000 r／min時(shí)達(dá)到最小，其值為8．59%。

（3）在后續(xù)工作中，將會(huì)對(duì)超聲振動(dòng)輔助磨削脆性材料的切削力預(yù)測(cè)模型進(jìn)行進(jìn)一步的研究，以期能夠建立考慮振動(dòng)參數(shù)的數(shù)值預(yù)測(cè)模型以及綜合考慮塑性流動(dòng)去除和脆性斷裂去除的理論預(yù)測(cè)模型，為全鋯牙冠的加工提供指導(dǎo)。

［1］Piconi C，Maccauro G．Zirconia as a ceramic biomaterial［J］．Biomaterials，1999，20：1－25．

［2］Preis V，Behr M，Handel，etal．Wear performance of dental ceramics after grinding and polishing treatments［J］．Journal of the Mechanical Behavior of Biomedical Materials，2012，10：13－22．

［3］Manicone P F，Iommetti P R，Raffaelli L．An overview of zirconia ceramics：Basic properties and clinical applications［J］．Journal of Dentistry，2007，35：819－826．

［4］Kosovka O D，Vesna M，Slobodan D，et al．Dilemmas in zirconia bonding：A review［J］．Srpski Arhivza Celokupno Lekarstvo，2013，141：395－401．

［5］石連水，朱洪水，廖嵐．前牙烤瓷冠燒結(jié)前后收縮變化的測(cè)定［J］．口腔頜面修復(fù)學(xué)，2000，1（3）：147－148．

SHILian-shui，ZHUHong-shui，LIAO Lan．Determination of linear sintering shrinkage rate of PFMcrown［J］．Chinese Journal of Prosthodontics，2000，1（3）：147－148．

［6］Lauwers B，Bleicher F，Ten Haaf P，et al．Investigation of the process-material interaction in ultrasonic assisted grinding of ZrO2 based ceramicmaterials［J］．4th CIRP International Conference on High Performance Cutting，2010，2：59－64．

［7］呂東喜，黃艷華，唐永健，等．基于SPH算法的磨粒沖擊工件表面過(guò)程數(shù)值模擬［J］．振動(dòng)與沖擊，2013，32（7）：169－174．

LüDong-xi，HUANG Yan-hua，TANG Yong-jian et al．Simulating process of abrasive impacting a workpiece surface based on SPH method［J］．Journal of Vibration and Shock，2013，32（7）：169－174．

［8］鄭侃，肖行志，廖文和，等．完全燒結(jié)氧化鋯牙科陶瓷材料的旋轉(zhuǎn)超聲加工實(shí)驗(yàn)研究［J］．人工晶體學(xué)報(bào)，2013，42（9）：1864－1869．

ZHENG Kan，XIAO Xing-zhi，LIAO Wen-he．Study on rotary ultrasonic machining of sintered zirconia dental ceramics［J］．Journal of Synthetic Crystals，2013，42（9）：1864－1869．

［9］Liu D F，CongW L，Pei Z J，etal．A cutting forcemodel for rotary ultrasonic machining of brittle materials［J］．International Journal ofMachine Tools＆Manufacture，2012，52：77－84．

［10］Zhang C L，Zhang J F，F(xiàn)eng P F．Mathematical model for cutting force in rotary ultrasonic face milling of brittle materials［J］．International Journalof Advanced Manufacture Technology，2013，69（1／4）：161－170．

［11］王剛，萬(wàn)敏，劉虎，等．粒子群優(yōu)化模糊系統(tǒng)的銑削力建模方法［J］．機(jī)械工程學(xué)報(bào)，2011，47（13）：123－130．

WANG Gang，WAN Min，LIU Hu，etal．Modeling ofmilling force by using fuzzy system optimized by particle swarm algorithm［J］．Journal of Mechanical Engineering，2011，47（13）：123－130．

［12］Chandrasekaran M，Muralidhar M，Murali Krishna C，et al．Application of soft computing techniques in machining performance prediction and optimization：a literature review［J］．International Journal of Advanced Manufacturing Technology，2010，46（5－8）：445－464．

［13］Haykin S．神經(jīng)網(wǎng)絡(luò)原理［M］．北京：機(jī)械工業(yè)出版社，2004．

［14］張德豐．MATLAB神經(jīng)網(wǎng)絡(luò)應(yīng)用設(shè)計(jì)［M］．北京：機(jī)械工業(yè)出版社，2009．

［15］Malkin S，Hwang TW．Grinding mechanisms for ceramics［J］．Annals of CIRP，1996，45：569－580．

［16］馮冬菊，趙福令，徐占國(guó)，等．超聲波銑削加工材料去除率的理論模型［J］．中國(guó)機(jī)械工程，2006，17（13）：1399－1403．

FENG Dong-ju，ZHAO Fu-ling，XU Zhan-guo，et al．Mathematic model of material removal rate for ultrasonic milling［J］．Chinese Journal of Mechanical Engineering，2006，17（13）：1399－1403．

［17］Pei Z J，F(xiàn)erreira P M，Haselkorn M．Plastic flow in rotary ultrasonic machining of ceramics［J］．Journal of Materials Processing Technology，1995，48：771－777．

［18］王棟，馮平法，張承龍，等．KDP晶體各向異性對(duì)劃痕特性影響的實(shí)驗(yàn)研究［J］．人工晶體學(xué)報(bào)，2012，41（3）：568－572．

WANG Dong，F(xiàn)ENG Ping-fa，ZHANG Cheng-long，et al．Experimental research on the influence of the KDP crystal anisotropy on scratch characteristics［J］．Journal of Synthetic Crystals，2012，41（3）：568－572．

Prediction model of cutting force in ultrasonic vibration assisted grinding of zirconia ceram ics

XIAO Xing-zhi，ZHENG Kan，LIAOWen-he
（School of Mechanical Engineering，Nanjing University of Science and Technology，Nanjing 210094，China）

Ultrasonic vibration assisted grinding（UVAG）experiments on dental zirconia ceramicswere conducted by using single-factormethod．An index predictionmodel，a BP neural networks predictionmodel，aswell as a theoretical prediction model of cutting force in UVAG were proposed．A verification experimentwas carried out to study the prediction accuracy of these threemodels，and the causes of errorswere revealed．The results indicate that the BP neural networks prediction model has higher precision compared with the two others，and the relative error is only 9．60%．The prediction accuracy of theoreticalmodel is poor due to the neglect of plastic flow removal during UVAG．

ultrasonic vibration assisted grinding；dental zirconia ceramics；cutting force prediction；index model；BP neural networks；theoreticalmodel

TH－16

A

10．13465／j．cnki．jvs．2015．12．024

國(guó)家自然科學(xué)基金項(xiàng)目（51305206）；中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金（2012XQTR001）

2014－01－08 修改稿收到日期：2014－06－10

肖行志男，博士生，1990年生

鄭侃男，博士，講師，1983年生