從政,曹巖,宋泠霞,汪文婷
西安工業(yè)大學(xué)機(jī)電工程學(xué)院
TC11鈦合金材料是Ti-Al-Mo-Zr-Si系馬氏體型α+β兩相鈦合金,具有抗拉強(qiáng)度和屈服強(qiáng)度高、耐腐蝕和抗腐蝕性好的特點(diǎn),屬于難加工材料,在加工過(guò)程中存在切削力大、切削溫度高以及刀具磨損嚴(yán)重等問(wèn)題。
結(jié)合某型軍用盤(pán)類(lèi)零件的車(chē)削過(guò)程,考慮具體加工條件和經(jīng)濟(jì)成本等問(wèn)題,決定采用有限元法模擬加工環(huán)境。陶亮等[1]采用AdvantEdge軟件,研究了基于二維模型的Inconel 718高溫合金切削溫度分布和工藝參數(shù)優(yōu)化。馬浩騫等[2]采用ABAQUS軟件,基于二維模型研究了切削TC4鈦合金過(guò)程中的力-熱分布規(guī)律,并進(jìn)行了參數(shù)優(yōu)化。韓甲棟[3]研究了不同切削參數(shù)切削TC11鈦合金材料時(shí)的表面粗糙度和刀具磨損變化規(guī)律,并得出最優(yōu)切削參數(shù)組。朱曉偉[4]采用硬質(zhì)合金刀具切削TC11鈦合金,得到了切削力和刀具使用壽命的經(jīng)驗(yàn)公式。Weiwei Ming等[5]評(píng)估了TC11鈦合金精硬車(chē)削的可加工性,對(duì)比了精車(chē)TC4和TC11合金過(guò)程中切削力、表面完整性和刀具磨損的變化規(guī)律。Sulin Chen等[6]使用不同涂層刀具、潤(rùn)滑方式和切削參數(shù),以切削力、切削溫度和表面粗糙度為對(duì)象,得出了最佳切削參數(shù)。許松等[7]研究了摩擦系數(shù)對(duì)微銑削TC4鈦合金的影響規(guī)律,研究發(fā)現(xiàn),當(dāng)摩擦系數(shù)為0.1~0.9時(shí),切削力隨著摩擦系數(shù)增大而增大,其中,當(dāng)摩擦系數(shù)為0.1~0.4時(shí),切削力明顯增大;當(dāng)摩擦系數(shù)為0.4~0.9時(shí),切削力變化平緩。湯祁[8]研究了材料本構(gòu)模型和摩擦系數(shù)對(duì)殘余應(yīng)力的影響關(guān)系,結(jié)合有限元法和車(chē)削試驗(yàn)求出摩擦系數(shù),確定了45鋼的殘余應(yīng)力仿真模型,并進(jìn)行了實(shí)驗(yàn)驗(yàn)證。
目前,切削仿真模型多使用二維簡(jiǎn)化模型,三維車(chē)削模型較少。在切削TC11鈦合金時(shí),對(duì)大切削深度、小進(jìn)給量和低切削速度工藝參數(shù)下的切削力研究較少。本文深度結(jié)合實(shí)際加工條件,通過(guò)建立同等比例的三維車(chē)削仿真模型,研究了切削速度、切削深度和摩擦系數(shù)對(duì)切削力的影響規(guī)律。
ABAQUS有限元仿真軟件具有豐富的材料庫(kù)、單元庫(kù)和強(qiáng)大的非線性處理能力,能夠有效模擬和處理金屬切削過(guò)程,是目前研究金屬切削主要軟件之一。試驗(yàn)基于ABAQUS軟件建立如圖1和圖2所示的工件和刀具簡(jiǎn)化模型,工件外徑為80mm,內(nèi)徑為66mm,軸向厚度為5mm。
圖1 工件模型
圖2 刀具模型
對(duì)工件進(jìn)行分割處理,其中,切削區(qū)域軸向厚度為1mm,網(wǎng)格大小為0.18mm,工件基體軸向厚度為4mm,網(wǎng)格大小為1mm,網(wǎng)格類(lèi)型均為C3D8R,網(wǎng)格總數(shù)分別為37萬(wàn)和7450個(gè)。對(duì)切削區(qū)域和工件基體施加綁定約束,對(duì)工件整體施加旋轉(zhuǎn)載荷,真實(shí)模擬實(shí)際工件轉(zhuǎn)動(dòng)。
刀片型號(hào)為CNMG120404,網(wǎng)格類(lèi)型為C3D8R,總數(shù)為4540個(gè)。程序分析類(lèi)型為動(dòng)態(tài)-顯式,分析時(shí)間0.1s,放大倍數(shù)1000倍。刀具定義為剛體,工件為彈塑性變形體。接觸定義為主面刀具,從面切削區(qū)域。接觸屬性采用傳統(tǒng)的罰函數(shù)法。圖3為刀具和工件的裝配實(shí)例。
圖3 刀具與工件裝配模型
圖4為檢驗(yàn)?zāi)P偷那邢髁ψ兓瑱z驗(yàn)參數(shù)來(lái)自于文獻(xiàn)[6],建立了直徑120mm、長(zhǎng)度90mm的工件和CNMG120408刀具模型。采用文獻(xiàn)[6]中的C1切削條件,以第五組正交實(shí)驗(yàn)為檢驗(yàn)條件,求得平均切削力的仿真值為166.24N,文獻(xiàn)中的切削力為149.8N,誤差為10.97%。產(chǎn)生誤差的主要原因是仿真環(huán)境較為理想,模型、材料屬性、接觸和網(wǎng)格密度均會(huì)產(chǎn)生一定誤差,但誤差在工程仿真可接受范圍內(nèi)。
圖4 切削力隨切削時(shí)間的變化曲線
金屬切削是復(fù)雜的彈塑性變形過(guò)程,其中包含了大應(yīng)變、大變形、力-熱耦合等非線性因素。在切削過(guò)程中,材料的變形、擠壓以及摩擦?xí)a(chǎn)生較大的切削力和大量熱量。結(jié)合ABAQUS軟件設(shè)置現(xiàn)有材料屬性,選用Johnson-Cook本構(gòu)模型(以下簡(jiǎn)稱(chēng)J-C本構(gòu)模型)進(jìn)行研究。TC11鈦合金工件的J-C本構(gòu)模型參數(shù)見(jiàn)表1[9]。 J-C本構(gòu)模型是將影響流動(dòng)應(yīng)力的應(yīng)變硬化效應(yīng)、應(yīng)變率效應(yīng)和溫度效應(yīng)相互聯(lián)系起來(lái),可以有效模擬切削過(guò)程[9]。
J-C本構(gòu)模型表達(dá)式為
(1)
表1 TC11鈦合金的J-C參數(shù)
仿真分析中,基于大切削深度、低進(jìn)給量和低切削速度的加工工況選取工藝參數(shù),在考察實(shí)際加工條件并與工藝人員進(jìn)行溝通后確定。選取切削速度、切削深度和摩擦系數(shù)作為研究變量,切削過(guò)程中的切削力作為研究對(duì)象。
采用單因素實(shí)驗(yàn)法研究三因素對(duì)切削力的影響規(guī)律。表2為單因素實(shí)驗(yàn)表。第一組實(shí)驗(yàn)是在切削深度為2mm、摩擦系數(shù)為0.1時(shí)進(jìn)行的切削速度單變量分析;第二組實(shí)驗(yàn)是在切削速度70m/min、摩擦系數(shù)0.1時(shí)進(jìn)行的切削深度單變量分析;第三組實(shí)驗(yàn)是在切削速度為70m/min、切削深度為2mm時(shí)進(jìn)行的摩擦系數(shù)單變量分析。三組單因素實(shí)驗(yàn)進(jìn)給量均為0.07mm/r。在單因素實(shí)驗(yàn)中,分別選取變化單因素組的四個(gè)水平值進(jìn)行實(shí)驗(yàn)。
表2 單因素實(shí)驗(yàn)安排
結(jié)合具體實(shí)際工藝參數(shù),使用正交實(shí)驗(yàn)法探尋在三因素交互影響下的切削力變化規(guī)律,以彌補(bǔ)單因素實(shí)驗(yàn)法的不足。正交實(shí)驗(yàn)法具有均勻分散、齊整可比的特點(diǎn),可以有效減少實(shí)驗(yàn)次數(shù)。表3為正交實(shí)驗(yàn)因素水平表,正交實(shí)驗(yàn)進(jìn)給量選為0.1mm/r。
表3 正交實(shí)驗(yàn)因素水平
采用單因素實(shí)驗(yàn),選取刀具穩(wěn)態(tài)切削階段的平均切削力(由于刀具定義為剛體,因此刀具上每點(diǎn)切削力相同),依據(jù)表2中數(shù)據(jù)繪制的切削力變化曲線如圖5所示。
圖5 單因素實(shí)驗(yàn)中的切削力變化曲線
由表2和圖5可以看出,隨著切削深度增大,切削力顯著增大,增長(zhǎng)達(dá)510N;隨著切削深度增大,單位切削體積和切削力增大。由切削速度變化對(duì)切削力的影響可以看到,切削速度為70m/min時(shí)切削力最大;切削速度為50~70m/min時(shí),切削力緩慢增加,增大3N;切削速度為70~110m/min時(shí),切削力緩慢降低,減少18N,拐點(diǎn)為70m/min時(shí),因此在實(shí)際加工中應(yīng)盡量避免選取這個(gè)切削速度。
隨著摩擦系數(shù)增大,切削力緩慢增大,共增大36N,增長(zhǎng)趨勢(shì)與文獻(xiàn)[7]相符。當(dāng)摩擦系數(shù)增大時(shí),同一切削時(shí)刻與刀具接觸的切屑網(wǎng)格材料變多,切屑網(wǎng)格更難從工件剝離,網(wǎng)格變得更粘,但由于切-屑分離準(zhǔn)則一定,因此摩擦導(dǎo)致的切削力增長(zhǎng)幅度緩慢。降低摩擦系數(shù)可以減少切削力,研究摩擦系數(shù)對(duì)切削力的影響規(guī)律可為選用切削液和改進(jìn)切削方式提供另一種思路。
對(duì)于上述單因素影響規(guī)律,在實(shí)際加工過(guò)程中,不同因素往往交互影響切削力。根據(jù)表3進(jìn)行正交實(shí)驗(yàn)的結(jié)果見(jiàn)表4,正交實(shí)驗(yàn)極差分析見(jiàn)表5。圖6為正交實(shí)驗(yàn)因素對(duì)切削力的影響規(guī)律。
表4 正交實(shí)驗(yàn)結(jié)果
表5 切削力極差分析
圖6 正交實(shí)驗(yàn)因素對(duì)切削力的影響
由表5的切削力極差分析可知,正交實(shí)驗(yàn)因素對(duì)切削力的影響排序?yàn)榍邢魃疃?摩擦系數(shù)>切削速度。切削深度和摩擦系數(shù)的增大均使切削力增大,變化趨勢(shì)與單因素影響規(guī)律相符,分別增大534.4N和76.6N。隨著切削速度增大,切削力的變化幅度在20N以?xún)?nèi),切削力峰值落在60~70m/min區(qū)間??紤]到圖5中的切削力最大值出現(xiàn)在切削速度為70m/min時(shí),因此在實(shí)際加工中,可以盡量避免切削速度落在70m/min附近。
利用SPSS軟件對(duì)表4的正交實(shí)驗(yàn)數(shù)據(jù)進(jìn)行處理,建立切削力的回歸模型,并對(duì)切削力的回歸模型進(jìn)行方差檢驗(yàn)。切削力公式為
F=-58.6-0.385x1+118.78x2+248.0x3
(2)
式中,F(xiàn)為切削力;x1為切削速度;x2為切削深度;x3為摩擦因素。
切削力回歸模型調(diào)整后的R2為99.83%,模型公式高度擬合數(shù)據(jù),切削力公式可對(duì)三因素變化區(qū)間切削力進(jìn)行準(zhǔn)確預(yù)測(cè)。表6為切削力回歸模型檢驗(yàn)[10]。
表6 切削力回歸模型方差檢驗(yàn)
在實(shí)際加工中,切削力過(guò)大會(huì)產(chǎn)生刀具快速磨損、工件表面質(zhì)量降低和破壞機(jī)床自身剛性等問(wèn)題,因此選取最小切削力為目標(biāo),對(duì)式(2)進(jìn)行工藝參數(shù)優(yōu)化。
根據(jù)正交實(shí)驗(yàn)的極差分析,只能求得水平組上的最優(yōu)參數(shù)組合,而采用算法工具可以在因素的整個(gè)變化區(qū)間內(nèi)進(jìn)行尋優(yōu)。因此使用MATLAB軟件,借助遺傳算法進(jìn)行優(yōu)化。算法設(shè)置為迭代次數(shù)100,交叉因子0.8,變異概率0.2,交叉分布指數(shù)20,變異分布指數(shù)20。
圖7為個(gè)體間的平均距離圖??芍诘螖?shù)超過(guò)30次后,迭代結(jié)果趨于穩(wěn)定。經(jīng)過(guò)多次重復(fù)求解,發(fā)現(xiàn)最終求得的最優(yōu)參數(shù)組為切削速度80m/min,切削深度2mm,摩擦系數(shù)0.1,對(duì)應(yīng)的切削力為172.96N。可見(jiàn),算法工具可以有效優(yōu)化工藝參數(shù),對(duì)實(shí)際生產(chǎn)加工提供切削用量選取指導(dǎo)。
圖7 個(gè)體間的平均距離
針對(duì)TC11鈦合金在車(chē)削過(guò)程中存在的難加工問(wèn)題進(jìn)行了仿真模擬分析,并對(duì)工藝參數(shù)進(jìn)行了算法優(yōu)化,主要得出以下結(jié)論。
(1)采用實(shí)際工件尺寸和刀具型號(hào)建立了同等比例的三維車(chē)削模型。在不影響仿真結(jié)果的前提下對(duì)三維模型進(jìn)行了優(yōu)化簡(jiǎn)化,提高了仿真計(jì)算效率。提出了基于ABAQUS有限元軟件的三維車(chē)削方法,并對(duì)三維車(chē)削模型進(jìn)行檢驗(yàn),檢驗(yàn)結(jié)果中,切削力誤差為10.97%,可見(jiàn)該模型可以準(zhǔn)確預(yù)估車(chē)削試驗(yàn)。
(2)采用單因素實(shí)驗(yàn)法和正交實(shí)驗(yàn)法對(duì)三維車(chē)削切削力進(jìn)行了仿真分析,得出了單因素條件和正交實(shí)驗(yàn)條件下的切削深度、切削速度和摩擦因素對(duì)切削力的影響規(guī)律??芍瑔我蛩胤ê驼粚?shí)驗(yàn)法切削力的變化規(guī)律相符,切削速度對(duì)切削力的影響均為先增大后減少,極值出現(xiàn)在70m/min附近,摩擦系數(shù)和切削深度的增大均使得切削力增大。正交試驗(yàn)因素對(duì)切削力的影響排序?yàn)榍邢魃疃?摩擦系數(shù)>切削速度。不同參數(shù)水平下的切削力可以有效地為實(shí)際加工提供指導(dǎo)。
(3)以最小切削力為前提,使用MATLAB軟件遺傳算法對(duì)正交實(shí)驗(yàn)因素水平進(jìn)行尋優(yōu),得到的最優(yōu)參數(shù)組為切削速度80m/min,切削深度2mm,摩擦系數(shù)0.1,切削力172.96N。優(yōu)化后的切削參數(shù)可以有效減少切削力、降低刀具磨損和提高加工表面質(zhì)量。