沈 鈺， 白海清

(陜西理工大學 機械工程學院， 陜西 漢中 723000)

鈦合金因具有熱強度高、抗腐蝕性好、耐熱性強等特點，被廣泛應用于航空航天制造及醫(yī)療器械等領(lǐng)域。但是，鈦合金作為一種典型的難加工材料，其導熱性差、彈性模量小、化學活性強致使其易與刀具材料產(chǎn)生化學反應，切削溫度高導致刀具磨損嚴重、加工效率低[1]。麻花鉆鉆削鈦合金材料時，鉆頭的幾何參數(shù)、刀具材料和切削用量是影響鈦合金鉆孔質(zhì)量的主要因素，鉆削力、鉆削溫度及刀具磨損則是反映鉆削性能及加工過程的重要依據(jù)。

對于如何提高麻花鉆鉆削鈦合金的鉆孔質(zhì)量，一直以來都是諸多學者研究的重點。韓榮第等[2]分析了鉆頭直徑、鉆削參數(shù)和刀具材料對鈦合金鉆削力的影響，結(jié)果表明鉆頭直徑對鉆削力影響最大，進給量次之，切削速度最??；南成根等[3]分析了碳纖維復合材料/鈦合金疊層鉆孔質(zhì)量損傷的原因，結(jié)果表明鉆削鈦合金時產(chǎn)生的切屑是導致CFRP孔徑超差的主要原因；胡立湘等[4]采用TiAlN涂層硬質(zhì)合金麻花鉆對鈦合金(TC4)進行鉆削試驗，獲得了切削用量對軸向力、加工孔壁表面粗糙度和形態(tài)的影響規(guī)律。

然而，以上研究采用的試驗方式均為傳統(tǒng)方式，存在周期長、成本高、場地受限等局限。故本文以鈦合金(Ti6Al4V)為研究對象，采用高效便捷的有限元仿真技術(shù)來彌補傳統(tǒng)試驗方式的不足。

本文基于DEFORM-3D軟件，并結(jié)合正交試驗方法，以鉆頭直徑d、切削速度vc、進給量f這3個變量為影響因素，研究麻花鉆鉆削鈦合金過程中，麻花鉆軸向力的變化。通過極差和方差分析，確定各因素的主次順序及最優(yōu)參數(shù)組合，并檢驗各因素的顯著性。根據(jù)正交試驗數(shù)據(jù)，結(jié)合MATLAB軟件進行多元線性回歸分析，建立鉆削力的數(shù)學模型，并對經(jīng)驗公式計算與仿真數(shù)據(jù)對比分析。

1 DEFORM-3D鉆削有限元仿真

1.1 麻花鉆模型的建立

對于麻花鉆模型的建立，本文采用DEFORM-3D軟件自身的三維建模功能，通過文獻[5]確定麻花鉆的幾何參數(shù)和刃磨參數(shù)后默認建模，其參數(shù)設(shè)置界面如圖1所示。對于DEFROM-3D軟件中標準麻花鉆默認建模的合理性，可通過麻花鉆的橫刃方程公式(1)和后角方程公式(3)計算驗證[6]，驗證其建模后的標準麻花鉆后角和橫刃斜角是否在合理值范圍內(nèi)，如果角度值合理，則進一步證明了標準麻花鉆模型的準確性。

橫刃方程：

(1)

(2)

后角方程：

(3)

(4)

1.2 預處理參數(shù)設(shè)置

1.2.1 刀具與工件參數(shù)設(shè)置

本文機械加工類型選擇鉆孔采用國際單位制SI標準，刀具設(shè)為剛性模型，材料選用硬質(zhì)合金；工件設(shè)為塑性模型，其形狀選擇直徑為16 mm，高度為6 mm的圓柱體，材料選用鈦合金(Ti6Al4V)，主要性能參數(shù)如表1所示，金屬材料的本構(gòu)模型采用Johnson-Cook(J-C)模型。DEFORM-3D軟件采用自適網(wǎng)格劃分(AMG)技術(shù)，其目的是為了保證在有限元分析過程中出現(xiàn)的不合格單元形狀的網(wǎng)格能夠即時重新劃分，確保網(wǎng)格的收斂性，從而確保仿真結(jié)果的準確性[7]。刀具采用相對網(wǎng)格劃分，設(shè)置網(wǎng)格數(shù)為20 000，size ratio為4；工件采用絕對網(wǎng)格劃分，設(shè)置最小單位邊長小于進給量的1/2，即最小網(wǎng)格尺寸小于單刃進給量，size ratio為7。在邊界條件中，設(shè)置工件側(cè)面與底面在X、Y、Z方向上的速度均為0，其目的是確保工件靜止；刀具沿-Z方向做進給運動并繞自身軸線旋轉(zhuǎn)。

表1 鈦合金的主要性能參數(shù)

圖2 鉆削有限元仿真模型

1.2.2 工作狀態(tài)參數(shù)設(shè)置

在DEFORM-3D軟件中，環(huán)境溫度設(shè)為20 ℃；對流系數(shù)為0.02 N/s/mm/C；定義摩擦類型為剪切摩擦，摩擦系數(shù)設(shè)為0.6；熱傳導系數(shù)為45 N/s/mm/C。

1.2.3 模擬參數(shù)設(shè)置

模擬運算總步數(shù)10 000步，每25步存儲1次；求解器采用共軛梯度求解，迭代方法為直接迭代法；刀具磨損模型采用經(jīng)驗公式Usui模型，設(shè)置參數(shù)a為1×10-6，b為855。

在Database generation對話框中，單擊Check按鈕，數(shù)據(jù)庫檢查確認無誤后，單擊Generate按鈕生成DB文件，生成的鉆削有限元仿真模型見圖2所示。

2 正交試驗設(shè)計

2.1 試驗內(nèi)容及方法

本文的目的是研究鉆頭直徑d(mm)、鉆削速度vc(m/min)、進給量f(mm/r)3個因素對麻花鉆軸向力的影響，故以搖臂鉆床Z3040×16為試驗目標，采用3水平3因素的試驗因素水平表，如表2所示。

表2 試驗因素水平表

采用L9(34)的正交試驗表，利用DEFORM-3D軟件對9組數(shù)據(jù)進行有限元仿真分析，并以麻花鉆的軸向力Fz(kN)作為鉆削性能指標，其鉆削正交試驗安排及結(jié)果如表3所示。

2.2 鉆削軸向力分析

對于麻花鉆軸向力在鉆削過程中的變化趨勢如圖3所示，其中X軸代表時間(Time，單位為s)，Y軸代表軸向力(Z Load，單位為kN)。由圖分析可知，軸向力變化趨勢先是不斷上升，后逐漸趨于平穩(wěn)，并在平穩(wěn)階段內(nèi)上下波動。分析其原因可能是由于網(wǎng)格即時重新劃分時造成仿真過程的不連續(xù)和切屑的變形或斷裂而引起的[6]。

表3 鉆削正交試驗安排及結(jié)果

(a) d=6 mm，vc=12 m/min， (b) d=6 mm，vc=18 m/min， (c) d=10 mm，vc=12 m/min， (d) d=10 mm，vc=24 m/min， f=0.20 mm/r f=0.24 mm/r f=0.28 mm/r f=0.24 mm/r圖3 鉆削軸向力變化曲線圖

2.3 正交試驗結(jié)果的極差與方差分析

表4 極差分析表

注：kij為Kij的平均值；極差R為kij中最大值與最小值之差。

根據(jù)正交試驗的仿真結(jié)果，本文對于軸向力仿真數(shù)據(jù)進一步進行了極差分析，以圖表的形式來確定各因素的主次順序及最優(yōu)鉆削參數(shù)組合。然而，各因素對試驗結(jié)果的影響程度，即各因素作用的顯著性，還無法準確地體現(xiàn)出來[8]。因此，還需對試驗結(jié)果進行方差分析，從而彌補上述分析的不足。分析結(jié)果如表4、表5、圖4所示，表5中各參數(shù)計算公式如下[9]：

fj=n-1,σj=SSj/fj,Fj=σj/σ, (6)

式中：SSj為平方和；Kij為第j因素第i水平所對應的試驗指標和；T為試驗結(jié)果之和；fj為自由度；n為各因素水平個數(shù)；σj為均方；σ為總誤差均方。

表5 方差分析表

(a) 鉆頭直徑的影響 (b) 鉆削速度的影響 (c) 進給量的影響圖4 鉆頭直徑、鉆削速度和進給量對鉆頭軸向力的影響趨勢圖

由圖4分析可知，在麻花鉆鉆削鈦合金的過程中，鉆頭的軸向力隨鉆頭直徑d的增大，先增大后減??；隨鉆削速度vc的增大，先減小后增大；隨進給量f的增大而增大。根據(jù)極差分析表(表4)及方差分析表(表5)可知，極差R的大小可判斷選定的各因素的主次順序，R值越大，表示該因素的水平變化對試驗結(jié)果的影響越大。因此，進給量f(因素C)對鉆頭軸向力的影響最大且最為顯著，鉆頭直徑d(因素A)次之，鉆削速度vc(因素B)影響最小。根據(jù)平均值大小選定各因素的最優(yōu)水平，平均值越大，影響越大。在鉆削過程中，麻花鉆的軸向力主要是由橫刃產(chǎn)生，而軸向力越小，越有利于提高橫刃的定心精度，同時也能提高刀具的耐用度，延長使用壽命。因此，確定最優(yōu)鉆削參數(shù)組合為A1B2C1，即d=6 mm，vc=18 m/min，f=0.20 mm/r。

3 鉆削力模型的建立

根據(jù)本文正交試驗數(shù)據(jù)，利用MATLAB軟件對其進行多元線性回歸分析，建立鉆削力的數(shù)學模型。根據(jù)金屬切削原理，切削力與切削參數(shù)之間存在復雜的指數(shù)函數(shù)關(guān)系，故軸向力的通用形式為[10]：

(7)

式中：Fz為軸向力；CF為系數(shù)；XF、YF、ZF為指數(shù)。

對公式兩邊取對數(shù)換算，使之變換成線性函數(shù)：lgFz=lgCF+XFlgd+YFlgvc+ZFlgf，進一步化簡，令y=lgFz，b0=lgCF，b1=XF，b2=YF，b3=ZF，x1=lgd，x2=lgvc，x3=lgf，則軸向力通用公式對應的回歸方程為：

y=b0+b1x1+b2x2+b3x3。

(8)

將正交試驗數(shù)據(jù)全部取對數(shù)后，遵循相應編程格式輸入MATLAB軟件中建立M文件，并采用多元線性回歸函數(shù)regress函數(shù)對數(shù)據(jù)進行處理，具體操作如下[11]：

x1=[0.7782,0.7782,0.7782,0.9031,0.9031,0.9031,1,1,1];

x2=[1.0791,1.2553,1.3802,1.0791,1.2553,1.3802,1.0791,1.2553,1.3802];

x3=[-0.6990,-0.6198,-0.5528,-0.6198,-0.5528,-0.6990,-0.5528,-0.6990,-0.6198];

y=[3.0273,3.0045,3.1359,3.1781,3.2394,3.1014,3.2243,2.9613,3.0572];

X=[ones(length(y),1),(x1)’,(x2)’,(x3)’];

Y=y’;

[b,bint,r,rint,stats]=regress(Y,X);

b,bint,stats

程序運行后可得軸向力的回歸方程為：y=3.890 3+0.151 7x1-0.167 8x2+1.145 9x3，

故軸向力的數(shù)學模型為：

Fz=7 767.835d0.151 7vc-0.167 8f1.145 9。

(9)

將選定的最優(yōu)鉆削參數(shù)組合代入DEFORM軟件中進行仿真，獲得仿真數(shù)據(jù)0.83，利用經(jīng)驗公式計算得0.99，比較兩者的誤差大小為16.16%，證明了經(jīng)驗公式的可行性。

4 結(jié) 語

(1)利用DEFORM-3D軟件對鈦合金材料鉆削加工進行的有限元仿真試驗，相比傳統(tǒng)試驗方式更加便捷、高效,且隨著計算機及有限元分析軟件的不斷發(fā)展與完善，其仿真數(shù)據(jù)的精度將進一步提高；

(2)對于鉆削參數(shù)對軸向力的影響，采用正交試驗法，選取代表性試驗點，并結(jié)合極差分析和方差分析，直觀確定了各因素的主次順序，得到進給量對麻花鉆的軸向力影響最為顯著，獲得了最佳鉆削參數(shù)組合，即d=6 mm，vc=18 m/min，f=0.20 mm/r；

(3)根據(jù)正交試驗結(jié)果，利用MATLAB軟件對其進行了多元線性回歸分析，建立了鉆削力的數(shù)學模型，通過利用經(jīng)驗公式與仿真數(shù)據(jù)的對比誤差分析，進一步證明了經(jīng)驗公式的可行性；

(4)通過鈦合金鉆削過程的有限元分析，為麻花鉆鉆削鈦合金材料的實際鉆削加工過程中提高鉆孔的質(zhì)量、合理選擇鉆削參數(shù)提供了參考依據(jù)，為鉆削參數(shù)的優(yōu)化打下了基礎(chǔ)。

[參考文獻]

[1] 李克強.鈦合金車削用PCD刀具優(yōu)選及切削參數(shù)優(yōu)化[D].哈爾濱:哈爾濱理工大學,2013.

[2] 韓榮第,吳健.鈦合金TC4的鉆削力試驗研究[J].工具技術(shù),2009,43(1):25-27.

[3] 南成根,吳丹,馬信國,等.碳纖維復合材料/鈦合金疊層鉆孔質(zhì)量研究[J].機械工程學報,2016,52(11):177-185.

[4] 胡立湘,李鵬南,陳明,等.TiAlN涂層鉆頭鉆削鈦合金的實驗研究[J].宇航材料工藝,2016,46(2):48-52.

[5] 戴俊平.麻花鉆內(nèi)錐面刃磨試驗[J].陜西理工學院學報(自然科學版),2003,19(1):7-9.

[6] 朱超,白海清.采用DEFORM-3D對鉆削過程的有限元仿真[J].煤礦機械,2015,36(6):286-289.

[7] 高興軍,李萍,閆鵬飛,等.基于Deform 3D不銹鋼鉆削機理的仿真研究[J].工具技術(shù),2011,45(4):17-20.

[8] 金浩,程寓,高超,等.玻璃鋼復合材料夾層結(jié)構(gòu)鉆削試驗研究[J].機械設(shè)計與制造,2011(4):133-135.

[9] 張利新,沈興全,張曉,等.基于正交試驗的麻花鉆鉆削鈦合金的刃磨參數(shù)選擇[J].現(xiàn)代制造工程,2013(6):8-11.

[10] 陸建中,孫家寧.金屬切削原理與刀具[M].北京:機械工程出版社,2011.

[11] 朱超,白海清.基于DEFORM-3D和遺傳算法的鉆削用量優(yōu)化研究[J].工具技術(shù),2016,50(1):48-51.