岳彩旭,黃 翠,劉二亮
(哈爾濱理工大學(xué)機(jī)械動(dòng)力工程學(xué)院,哈爾濱 150080)
永磁電機(jī)轉(zhuǎn)子護(hù)套用高溫合金Inconel625切削性能仿真與試驗(yàn)研究
岳彩旭,黃 翠,劉二亮
(哈爾濱理工大學(xué)機(jī)械動(dòng)力工程學(xué)院,哈爾濱 150080)
本文針對(duì)高速永磁電機(jī)轉(zhuǎn)子保護(hù)材料鎳基高溫合金Inconel625的切削性能展開研究。采用仿真軟件DEFORM-3D建立了該材料的三維有限元仿真模型,仿真模型采用了Johnson-Cook模型、Usui刀具磨損模型等關(guān)鍵技術(shù)。通過(guò)仿真得到了切削過(guò)程中切削力的變化趨勢(shì)、溫度場(chǎng)的分布特性、刀具磨損、切屑形狀等。在結(jié)合試驗(yàn)分析的基礎(chǔ)上,得到了切削參數(shù)對(duì)切削力、切削溫度、刀具磨損的影響規(guī)律。研究結(jié)果對(duì)Inconel625合金車削過(guò)程切削條件的合理選擇和刀具使用壽命的提高提供了理論依據(jù),為深入研究Inconel625合金切削機(jī)理提供了理論依據(jù)。
Inconel625;切削性能;有限元仿真;刀具磨損
高速永磁電機(jī)的轉(zhuǎn)子轉(zhuǎn)速高達(dá)十幾萬(wàn)轉(zhuǎn),轉(zhuǎn)子受離心力作用非常大,且工作溫度較高,因此,通常采用轉(zhuǎn)子護(hù)套對(duì)轉(zhuǎn)子進(jìn)行保護(hù)。轉(zhuǎn)子表面護(hù)套常采用機(jī)械強(qiáng)度高、熱穩(wěn)定性優(yōu)良、熱導(dǎo)率較高的材料,常用的護(hù)套材料如碳纖維、Inconel系列高溫合金以及鈦合金等[1]。與復(fù)合材料相比,高溫合金在熱穩(wěn)定性、熱強(qiáng)性和散熱能力等方面上具有較為明顯的優(yōu)勢(shì)[2]。Inconel625(相近牌號(hào)GH3625、GH625)鎳基高溫合金是以Mo、Nb為主要強(qiáng)化元素的固溶強(qiáng)化型鎳基變形高溫合金,在高溫下仍能保持其優(yōu)良的綜合機(jī)械性能,較為適合用作高速永磁電機(jī)的轉(zhuǎn)子保護(hù)材料[3]。鎳基高溫合金在切削加工過(guò)程中具有切削力大、切削溫度高、刀具磨損嚴(yán)重、加工硬化傾向大、不易斷屑等加工特點(diǎn),從而使得其切削加工的難度大大增加,屬于典型的難加工材料之一[3]。因此對(duì)Inconel625合金的切削過(guò)程進(jìn)行深入研究具有重要的理論和應(yīng)用價(jià)值。
國(guó)內(nèi)外學(xué)者對(duì)該類合金進(jìn)行了很多的研究,其中包括熱變形行為特征、材料損傷機(jī)理、加工工藝和LAMM加工等方面,為難加工材料鎳基高溫合金的加工問(wèn)題提供了有效的解決措施,并為其切削研究提供了大量的理論依據(jù)及技術(shù)支持。國(guó)外學(xué)者R.Ramanujam等在Inconel625干式車削中結(jié)合模糊函數(shù)和正交實(shí)驗(yàn)設(shè)計(jì)的混合方法,以最小表面粗糙度、最小能耗以及最大的材料去除率為指標(biāo)對(duì)加工參數(shù)進(jìn)行了優(yōu)化[4]。M.Hokka等利用試驗(yàn)獲得的J-C材料本構(gòu)模型建立了625鎳基高溫合金正交切削仿真模型,并進(jìn)一步修正了625合金的本構(gòu)模型。雖然仿真很好的模擬了切屑分離和鋸齒狀切屑的生成,但仿真結(jié)果中的切削應(yīng)力和切屑形狀與實(shí)驗(yàn)結(jié)果還沒(méi)有很好的吻合,為改進(jìn)仿真模型的精度提供良好的理論及實(shí)驗(yàn)基礎(chǔ)[5]。Bogdan S?odki利用WNMG 080404-23和VCMT 160404-SM兩種刀片對(duì)難加工材料Inconel625進(jìn)行了車削,通過(guò)觀察切屑的形狀,提出了一種正確選擇切削參數(shù)的算法,并將該算法應(yīng)用于仿真中,離線驗(yàn)證了槽型的填充過(guò)程[6]。Avik Samanta等通過(guò)實(shí)驗(yàn)及仿真對(duì)Inconel625進(jìn)行了激光輔助機(jī)械微加工(LAMM)以獲得實(shí)驗(yàn)過(guò)程表征并對(duì)切削力及殘余應(yīng)力進(jìn)行了預(yù)測(cè),建立了激光輔助機(jī)械微加工熱-機(jī)械耦合場(chǎng)的有限元模型,并研究了激光功率、切削速度、切削刃半徑、前角、激光位置和激光光斑直徑對(duì)切削力和已加工表面殘余應(yīng)力的影響[7]。
國(guó)內(nèi)學(xué)者陶琳等通過(guò)等溫?zé)釅嚎s實(shí)驗(yàn)獲得了Inconel625合金的真應(yīng)力-應(yīng)變曲線,通過(guò)非線性回歸建立了合金在高溫高速條件下的本構(gòu)模型[8]。馬妙,陸永浩研究了載荷、振幅和頻率對(duì)Inconel625鎳基合金在常溫常壓下摩擦磨損性能的影響,并對(duì)磨痕表面形貌進(jìn)行了觀察,探討了損傷機(jī)理[9]。何艷華分析了Inconel625的切削加工特點(diǎn),提出了選用合適的刀具材料、刀具參數(shù)以及優(yōu)化后的切削用量來(lái)解決難加工材料的加工問(wèn)題[10]。劉逢博通過(guò)對(duì)氣相質(zhì)譜分析儀中Inconel625零件的切削加工工藝的分析研究,總結(jié)并驗(yàn)證出了合理有效的工藝方案[11]。
近年來(lái),有限元仿真技術(shù)越來(lái)越廣泛地應(yīng)用在金屬切削加工的模擬中。精確的有限元仿真模型不僅可以大量減少試驗(yàn)時(shí)間及費(fèi)用,而且能為真實(shí)的加工提供較為可信的參考。綜合分析國(guó)內(nèi)外研究現(xiàn)狀,對(duì)Inconel625合金切削機(jī)理方面的相關(guān)研究較少,尤其是將有限元仿真分析方法以及試驗(yàn)分析方法相結(jié)合的研究方法更為少見(jiàn),因此,本文基于Deform-3D有限元仿真軟件建立了Inconel625合金切削過(guò)程的有限元仿真模型,得出不同切削參數(shù)下主切削力、切削溫度以及后刀具磨損量等切削過(guò)程變量的仿真結(jié)果,并結(jié)合相應(yīng)試驗(yàn)研究,分析了切削參數(shù)對(duì)切削過(guò)程的影響規(guī)律。
1.1 仿真模型的建立
利用Deform-3D軟件的金屬切削模塊進(jìn)行仿真分析,利用軟件的前處理器進(jìn)行有限元仿真模型的建立,運(yùn)行所建立的模型,通過(guò)后處理器觀察相關(guān)的仿真結(jié)果,所采用的有限元仿真技術(shù)路線圖如圖1所示。

圖1 有限元仿真技術(shù)路線圖
利用三維制圖軟件UG對(duì)試驗(yàn)過(guò)程中使用的刀具進(jìn)行建模,為了保證仿真效率,刀片選截刀尖部分并導(dǎo)出.stl格式,再將其導(dǎo)入Deform-3D中作為仿真對(duì)象,如圖2所示。為了權(quán)衡模擬速度及計(jì)算精度,將刀具與工件的網(wǎng)格數(shù)目均劃分為10萬(wàn)左右,且對(duì)刀尖以及工件與刀尖接觸部分進(jìn)行網(wǎng)格的局部劃分,并設(shè)置加密窗口跟蹤刀具的運(yùn)動(dòng)軌跡,以軟件的自適應(yīng)網(wǎng)格劃分方法來(lái)進(jìn)行網(wǎng)格的重劃分。刀具及工件的網(wǎng)格劃分情況如圖3所示。

圖2 刀具模型的建立

圖3 刀具及工件的網(wǎng)格劃分
刀具基體材料為WC,涂層為coating-TiAlN,其材料屬性設(shè)置為剛性。工件材料為Inconel625,材料的主要物理性能見(jiàn)表1,其材料屬性設(shè)置為塑性。由于Deform-3D軟件的材料庫(kù)中并未封裝Inconel625此種材料,因此采用文獻(xiàn)[12]通過(guò)霍普金森壓桿實(shí)驗(yàn)獲得的可以反映在高溫、大應(yīng)變、大應(yīng)變率條件下的Inconel625的Johnson-Cook流動(dòng)應(yīng)力模型(1),各項(xiàng)參數(shù)見(jiàn)表2,再將其結(jié)合物理性能生成該材料的.key文件并導(dǎo)入仿真模型中應(yīng)用。

表1 Inconel625的物理性能

表2 通過(guò)霍普金森壓桿實(shí)驗(yàn)獲得的Johnson-Cook本構(gòu)模型參數(shù)
本文采用的Johnson-Cook本構(gòu)模型為:

摩擦模型定義為軟件自帶的shea(r剪切)摩擦模型,摩擦系數(shù)為0.6,模型為:

式中:fs為摩擦力,m為摩擦系數(shù),k為剪切屈服應(yīng)力。
刀具磨損模型采用軟件自帶的并廣泛應(yīng)用于金屬切削仿真中的Usui磨損模型:

式中:W為刀具磨損量,p為正壓力,V為工件材料相對(duì)于刀具的滑動(dòng)速度,T為刀面絕對(duì)溫度,a與b為特征常數(shù)(主要由切削參數(shù)及材料決定),根據(jù)文獻(xiàn)[13]查得:a=3.6e-9,b=1200。
1.2 仿真過(guò)程切削參數(shù)的確定
Inconel625鎳基高溫合金(下文簡(jiǎn)稱625合金)具有較好的高溫強(qiáng)度、塑性,良好的抗高溫氧化性、抗高溫腐蝕性等綜合性能,與此同時(shí)伴隨著較差的切削加工性,為了優(yōu)選出適合其半精加工的切削參數(shù),進(jìn)行了單因素切削參數(shù)設(shè)計(jì)進(jìn)而研究切削參數(shù)對(duì)切削過(guò)程變量的影響。在高速、低進(jìn)給以及小切深的條件下對(duì)625合金進(jìn)行仿真以及后續(xù)的試驗(yàn)研究,詳細(xì)切削參數(shù)見(jiàn)表3。

表3 切削參數(shù)的取值用量
1.3 仿真結(jié)果提取及分析
刀具壽命是切削高溫合金的關(guān)鍵問(wèn)題,而選用合理的切削用量對(duì)提高刀具耐用度又是十分重要的[14],這里我們主要提取625合金切削過(guò)程中不同切削參數(shù)下的切削力、刀具及工件的溫度、刀具磨損等仿真結(jié)果。圖4所示為在切削條件vc=80m/min、f=0.2mm/r、ap=0.5mm下,刀具、工件表面溫度、切削力變化趨勢(shì)、后刀面磨損形貌的仿真結(jié)果。圖5、6、7分別為切削參數(shù)對(duì)刀具及工件的溫度、刀具磨損的影響,其中由于刀具磨損是隨著切削時(shí)間的累積而不斷變化的,因此,此處取相同的切削時(shí)間來(lái)提取各個(gè)切削條件下的后刀面磨損量的結(jié)果。


圖4 仿真過(guò)程中的物理量結(jié)果提取
圖5所示為固定f=0.2mm/r,ap=0.5mm,vc依次取40、60、80、100m/min下的刀具、工件溫度以及刀具磨損的仿真結(jié)果,分別提取的是刀尖處最高溫度以工件切削區(qū)域的平均溫度,以及同一切削時(shí)間下的后刀面磨損值。隨著切削速度的增大,刀具及工件的溫度呈上升趨勢(shì),原因在于當(dāng)切屑沿前刀面流出時(shí),切屑底層與前刀面發(fā)生強(qiáng)烈的摩擦,這些摩擦產(chǎn)生的熱量在很短的時(shí)間內(nèi)生成,來(lái)不及向切屑內(nèi)部傳導(dǎo),因此切削溫度會(huì)明顯的提高。
圖6所示為固定vc=40m/min,ap=0.5mm,f依次取0.10、0.15、0.20、0.25mm/r下的刀具、工件溫度以及刀具磨損的仿真結(jié)果,隨著進(jìn)給量的增大,刀具及工件的溫度呈上升趨勢(shì),原因在于進(jìn)給量的增大使得單位時(shí)間內(nèi)金屬切削量變多,所以產(chǎn)生的熱量多,但同時(shí)切屑變厚,傳熱面積增加,由此導(dǎo)致切削溫度稍緩上升的趨勢(shì)。

圖5 切削速度對(duì)切削溫度及刀具磨損的影響

圖6 進(jìn)給量對(duì)切削溫度及刀具磨損的影響
圖7所示為固定vc=80m/min,f=0.2mm/r,ap依次取0.5、0.75、1、1.25mm下的刀具、工件溫度以及刀具磨損的仿真結(jié)果。隨著切削深度的增大,刀具及工件的溫度呈緩慢上升的趨勢(shì),原因在于切削深度的增大雖然會(huì)增加大量的切削熱,但同時(shí)也會(huì)使切削刃參與工作的長(zhǎng)度增長(zhǎng),改善散熱條件,故溫升不明顯。

圖7 切削深度對(duì)切削溫度及刀具磨損的影響
由仿真所得到的結(jié)果可以看出,隨著切削速度、進(jìn)給量、切削深度的增大,刀具及工件的溫度以及后刀面磨損值均增大,并且三個(gè)參數(shù)對(duì)切削溫度及刀具磨損的影響程度是非常相似的。其影響程度為:切削速度的影響最大、進(jìn)給量次之、切削深度的影響最小。在切削過(guò)程中,刀-工接觸區(qū)里發(fā)生著強(qiáng)烈的摩擦,這樣會(huì)產(chǎn)生很高的溫度和壓力。因此,前、后刀面隨著切削的進(jìn)行會(huì)逐漸產(chǎn)生磨損。隨著前、后刀面的磨損,推擠力和摩擦力逐漸增大,功耗增加,產(chǎn)生的熱量增多,以使得切削溫度上升,由此可看出刀具磨損與切削溫度二者互相影響并相互作用,因此,二者受切削參數(shù)變化的影響程度也較為相似。
對(duì)前面所建立的仿真模型進(jìn)行了試驗(yàn)驗(yàn)證,提取了切削過(guò)程中的切削力、切屑、后刀面磨損形貌并與仿真結(jié)果進(jìn)行了對(duì)比分析。
2.1 試驗(yàn)系統(tǒng)
切削所用刀具為三菱公司生產(chǎn)的CNMG120408-MJ刀片,刀桿采用山特維克公司生產(chǎn)的DCLNR型刀桿,刀片刀尖部分前角為13°、主切削刃部分的前角為9°,后角為0°,刀尖圓弧半徑為0.8mm。切削參數(shù)同仿真參數(shù),切削試驗(yàn)在數(shù)控車床CA6140上進(jìn)行,采用瑞士KISTLER測(cè)力儀進(jìn)行切削力的測(cè)量,利用超景深顯微鏡對(duì)刀具磨損形貌進(jìn)行觀察,試驗(yàn)裝置如圖8所示。

圖8 切削性能試驗(yàn)裝置
2.2 試驗(yàn)與仿真結(jié)果對(duì)比分析
圖9和圖10所示為在切削參數(shù):vc=60m/min,f=0.2mm/r,ap=0.5mm條件下仿真與試驗(yàn)的后刀面磨損形貌以及切屑形狀的對(duì)比。

圖9 刀具磨損形貌對(duì)比

圖10 切屑形狀對(duì)比
從圖9、圖10中,我們可以看出仿真的后刀面磨損形貌與試驗(yàn)所得的后刀面磨損形貌基本一致,靠近刀尖處的后刀面磨損比較嚴(yán)重。仿真得出的切屑的卷屑形狀與試驗(yàn)所得的螺旋卷屑形狀基本相似。
圖11所示為仿真與試驗(yàn)后經(jīng)處理得到的主切削力平均值的數(shù)值對(duì)比以及所考察的三個(gè)切削參數(shù)對(duì)主切削力的影響。
如圖11所示,主切削力的仿真與試驗(yàn)結(jié)果在變化趨勢(shì)上是相同的。從圖11(a)可看出,隨著切削速度的提高,切削力呈減小趨勢(shì)。主要原因在于隨著切削速度的增加,切削溫度逐漸升高,前刀面摩擦系數(shù)減小,變形系數(shù)減小,使切削力逐漸減小,且漸趨穩(wěn)定。從圖11(b)、(c)可看出,隨著進(jìn)給量、切削深度的增大,主切削力增大,其中切削深度對(duì)切削力的影響程度比進(jìn)給量f要大。主要原因?yàn)殡m然隨著切削深度及進(jìn)給量的增大,切削面積都成正比增加,但是,切削深度增大時(shí),單位切削力不變,而進(jìn)給量增大時(shí),變形系數(shù)減小使得單位切削力減小,所以切削深度對(duì)切削力的影響更明顯??梢钥闯龇抡娴恼`差在10%~15%之間,究其主要原因?yàn)榉抡娴那疤幚韰?shù)如傳熱系數(shù)、摩擦系數(shù)等與真實(shí)的試驗(yàn)條件存在一定的差距,但誤差范圍在一定程度上是可以接受的,因此我們可以通過(guò)仿真來(lái)預(yù)測(cè)試驗(yàn)的可行性、相關(guān)切削條件的選擇,以及判斷整個(gè)切削過(guò)程物理量的變化趨勢(shì)等等。


圖11 主切削力的仿真與試驗(yàn)結(jié)果對(duì)比
本文針對(duì)永磁電機(jī)轉(zhuǎn)子護(hù)套用高溫合金Inconel625展開了切削性能研究。采用有限元軟件Deform-3D建立了高溫合金Inconel625三維切削過(guò)程有限元仿真模型,并通過(guò)試驗(yàn)驗(yàn)證了該仿真模型的精度。結(jié)合試驗(yàn)與仿真研究得出以下結(jié)論:
(1)車削Inconel625合金時(shí),隨著切削速度的增大,切削力逐漸減小。隨著切削深度、進(jìn)給量的增大,切削力均呈增大趨勢(shì),其中,切削深度對(duì)切削力的影響程度更大。仿真與試驗(yàn)得到的主切削力在變化趨勢(shì)上是一致的,數(shù)值上的誤差在15%以內(nèi),證明了所建立的仿真模型較為準(zhǔn)確。
(2)車削Inconel625合金時(shí),隨著切削速度、切削深度、進(jìn)給量的增大,切削溫度均呈增大趨勢(shì),隨著切削速度的提高,切削溫度上升幅度最大,進(jìn)給量次之,切削深度對(duì)切削溫度的影響最小。刀具磨損隨切削參數(shù)的變化趨勢(shì)與刀、工溫度隨切削參數(shù)的變化趨勢(shì)相似。
(3)車削Inconel625合金時(shí),中等的切削速度、較大的進(jìn)給量以及較小的切削深度為該材料較為適宜的半精加工切削條件。
[1]董劍寧,黃允凱,金龍,等.高速永磁電機(jī)設(shè)計(jì)與分析技術(shù)綜述[J].中國(guó)電機(jī)工程學(xué)報(bào),2014,(27):011.
[2]張鳳閣,杜光輝,王天煜,等.高速永磁電機(jī)轉(zhuǎn)子不同保護(hù)措施的強(qiáng)度分析[J].中國(guó)電機(jī)工程學(xué)報(bào), 2013,(33):195-202.
[3]趙秀芬,王玉華,劉陽(yáng),等.鎳基高溫合金的切削加工[J].航空制造技術(shù),2010,(11):46-50.
[4]Ramanujam R,Venkatesan K,Saxena V,et al. Optimization of Machining Parameters Using Fuzzy Based Principal Component Analysis during dry turning operation of Inconel 625-A hybrid approach [C].12th Global Congress on Manufacturing and Management,India:Elsevier Ltd,2014:668-676
[5]Hokka M,Gomon D,Shrot A,et al.Dynamic Behavior and High Speed Machining of Ti-6246 and Alloy 625 Superalloys:Experimental and Modeling Approaches[J].Experimental Mechanics,2014,54(2): 199-210.
[6]S?odki B.Chipformers reliability in Inconel 625 longitudinal turning[J].Management and Production Engineering Review,2013,4(2):93-101.
[7]Samanta A,TeliM,Singh R.Experimental characterization and finite element modeling of the residual stresses in laser-assisted mechanical micromachining of Inconel 625[J].Proceedings of the Institution of Mechanical Engineers Part B:Journal of EngineeringManufacture,2015,37(8):745-8.
[8]陶琳,程明,張偉紅,等.Inconel 625合金高溫高速熱變形行為[J].材料熱處理學(xué)報(bào),2012,33(9): 55-59.
[9]馬妙,陸永浩.Inconel625鎳基合金微動(dòng)磨損性能研究[J].摩擦學(xué)學(xué)報(bào),2012,(5):007.
[10]何艷華.鎳基合金的車削加工[J].金屬加工:冷加工,2010,(9):36-37.
[11]劉逢博.基于鎳基高溫合金加工工藝的研究與應(yīng)用[J].機(jī)電產(chǎn)品開發(fā)與創(chuàng)新,2010,(5):143-145.
[12]趙娜.高溫合金Inconel625本構(gòu)建模及切削過(guò)程仿真研究[D].哈爾濱理工大學(xué),2016
[13]Lotfi M,Jahanbakhsh M,FaridAA.Wear estimation of ceramic and coated carbide tools in turning of Inconel 625: 3D FE analysis[J]. Tribology International,2016,99:107-116.
[14]沈興東,黃景衡.高溫合金的可切削性分析與切削刀具設(shè)計(jì)[J].工具技術(shù),2006,40(9):50-54.
Simulation and Experimental Study on Cutting Performance of High Temperature Alloy Inconel625 for PM Motor Rotor Retaining Sleeve
YUE Caixu,HUANG Cui,LIU Erliang
(Mechanical Manufacture and Automation,Harbin University of Science and Technology, Harbin 150080,China)
In this paper,the cutting performance of high speed permanent magnet motor rotor protection material nickel base superalloy Inconel625 is studied.A 3D finite element simulation model is established by using the finite element simulation software DEFORM-3D.The Johnson-Cook model is used to establish the constitutive relation of workpiece material,and the tool wear model is established based on the Usui model.The change of cutting force,the distribution characteristics of temperature field,tool wear and chip formation are obtained by finite element simulation.Combined with the experimental results,the influence of cutting parameters on the cutting force,cutting temperature and tool wear are analyzed,and the simulation model is verified.The research results provide a theoretical basis for the reasonable selection of cutting conditions and the improvement of the tool life,they also offer some reference dates for further research on the mechanism of Inconel625 alloy cutting.
Inconel625;finite element simulation;cutting force;cutting temperature;tool wear
TM351
A
1000-3983(2017)01-0045-06
2016-09-15
1.國(guó)家科技重大專項(xiàng):航空發(fā)動(dòng)機(jī)盤環(huán)軸零件國(guó)產(chǎn)化成套刀具產(chǎn)品開發(fā)及應(yīng)用(2014ZX04001020-011);2.哈爾濱理工大學(xué)青年拔尖創(chuàng)新人才培養(yǎng)計(jì)劃資助(201507)