田榮鑫，姚倡鋒，武導俠

（西北工業(yè)大學現(xiàn)代設計與集成制造技術教育部重點實驗室，西安 710072）

高強度鋁合金由于其密度低、比強度高、耐磨性好、加工性能好以及焊接性能良好等特點，被廣泛應用于航空航天等領域[1]。但是鋁合金材質(zhì)較軟，加工易產(chǎn)生熱變形，容易粘刀，形成積屑瘤；線膨脹系數(shù)大，裝夾和加工時容易引起變形，影響尺寸精度。使用鋁合金制造飛機蒙皮、翼梁、隔框、長桁和起落架時，加工量非常大，而且許多精密部件需減小結構上的加工變形，鋁合金大型壁板的加工變形問題已成為制約航空業(yè)高速發(fā)展的一大問題。而加工變形主要是由切削過程中的切削力和切削溫度引起的，所以研究切削過程中的切削力和切削溫度很有必要。

國內(nèi)外學者對鋁合金高速銑削過程中的銑削力和銑削溫度進行了大量的試驗研究[2-4]。然而，隨著數(shù)值分析方法的發(fā)展及其在制造業(yè)中應用的推廣，Ng等[5]運用有限元分析方法研究了刀具磨損對工件的溫度和剪應力分布的影響。由Strenkowski等[6]提出的正交切削的歐拉有限元模型可以對工件、切屑和刀具的溫度進行預測。Xie等[7]通過ABAQUS有限元軟件研究了連續(xù)切屑的形成過程，進行熱量傳導分析后，獲得了穩(wěn)定切削狀態(tài)下的溫度分布。成群林等[8]對7050-T745鋁合金的高速銑削過程進行了有限元模擬，發(fā)現(xiàn)切削力和切削溫度的仿真數(shù)據(jù)和試驗數(shù)據(jù)具有很好的一致性。此外，一些學者通過切削過程的三維仿真研究[9-10]來解決二維仿真精度較低的問題。趙云峰等[11]采用有限元分析軟件AdvantEdge對鋁合金A12024銑削加工過程中銑削力變化和切削溫度的分布情況進行了研究，通過模擬獲得了三向銑削力的變化曲線及工件和刀具上的溫度分布規(guī)律，并進行了試驗驗證，發(fā)現(xiàn)仿真分析獲得結果與試驗結果吻合良好。成群林等[12]結合立銑加工的實際情況，建立了面向立銑加工的可用于預測不同切削參數(shù)組合下切削力的螺旋齒單刃斜角切削有限元模型。畢運波等[13]針對螺旋齒立銑刀加工過程，模擬分析了工件溫度場的分布和變化， 研究發(fā)現(xiàn)銑削熱只影響到加工表層附近很薄的一部分材料的溫升。

鋁合金易發(fā)生加工變形影響加工精度，而在高速切削范圍內(nèi)，切削力與切削速度成反比，因此高速銑削是減小鋁合金加工變形的有效手段。本文以高速銑削為背景，在使用新型仿真軟件AdvantEdge的基礎上，重點對高速銑削鋁合金7055的銑削力和銑削溫度進行了有限元仿真研究。

1 有限元模型的建立

運用美國Third Wave Systems 公司開發(fā)的AdvantEdge FEM有限元仿真軟件，仿真分析了7055鋁合金高速銑削過程中的力和溫度，銑削力和銑削溫度仿真模型如圖1所示。

圖1 有限元模型Fig.1 Finite element model

7055鋁合金假定為可延展的，根據(jù)冪強化法則建立的材料本構方程如下[9,14]：如果

如果

式中，為等效Mises應力，g為流動應力，ep為累積塑性應變，為累積塑性應變率；為基準塑性應變率，為低應變率敏感指數(shù)，m2為高應變率敏感指數(shù)，為塑性應變率分界點。

溫度軟化方程為：

式中，n為強化系數(shù)，T為材料表面溫度，T0為基準溫度，α為軟化系數(shù)，s0為基準溫度T0時的屈服應力，ep為累積塑性應變，為基準塑性應變。

2 銑削過程仿真

2.1 模擬參數(shù)設定

仿真過程中具體參數(shù)設定選擇：工件材料是7055鋁合金，切削類型為順銑，刀具選擇φ12mm 的三齒硬質(zhì)合金刀具，前角14°，后角10°，螺旋角40°，刀尖圓角半徑為0.04mm，冷卻條件為干切削，初始溫度20℃。銑削參數(shù)范圍：銑削速度vc400～1300m/min，每齒進給量fz0.02～0.2mm/z，銑削深度ap0.1～1.9mm，銑削寬度ae2～11mm。

2.2 仿真過程分析

銑削時影響切削力的主要因素是切削層的厚度，銑削過程中的切削層厚度直接影響金屬彈、塑性變形產(chǎn)生的銑削變形抗力，以及切屑與刀具的摩擦阻力。圖2為當vc=800m/min、fz=0.04mm/z、ap=0.5mm、ae=6mm 時，一個刀齒切入材料到切出材料獲得的切削力（Fx、Fy）和切削溫度（T）隨時間（t）變化的曲線。在切削開始階段，切削力隨著刀齒切入材料厚度的增加而迅速從0開始增大；當?shù)洱X完全切入材料時，切削力達到最大值，隨后又逐漸減小，這是由于切削厚度越來越??；當?shù)洱X完全切出材料時，切削力降為0，切削溫度和切削力的變化趨勢一致。

圖2 切削力和切削溫度隨時間變化Fig.2 Distribution of cutting force and temperature with time change

由于摩擦熱沿前刀面逐漸增加的緣故，銑削過程中刀具的最高溫度點處于與刀尖一定距離的前刀面處。vc=800m/min、fz=0.04mm/z、ap=0.5mm、ae=6mm條件下刀具前刀面和后刀面的溫度分布見圖3，可以看出，刀具的最高溫度點出現(xiàn)在距離刀尖點0.05～0.1mm的前刀面上。

vc=800m/min、fz=0.04mm/z、ap=0.5mm、ae=6mm 條件下，不同切削階段工件、刀具和切屑上的溫度場分布和變化見圖4，可以看出：（1）切削的起始階段，工件和刀具上溫度場分布范圍較窄，隨著切削過程的進行，切削熱逐漸由刀尖向刀具內(nèi)部擴散，刀具的溫度場范圍逐漸增大，并達到一個穩(wěn)定的狀態(tài)。（2）高速加工過程中，切屑溫度明顯高于工件和刀具的溫度，這是由于切屑帶走的熱量遠多于傳入刀具和工件中的熱量。（3）高速加工過程中，工件基體溫度變化不大，只有已加工表面的一個薄層溫度變化比較明顯。這是因為高速切削過程中，刀具和工件間的相對速度很大，刀具會很快離開已加工表面，使得傳入工件的熱量快速對流散失，工件已加工表面溫度場分布快速趨于穩(wěn)定。因此，針對鋁合金切削加工，高速切削是一種通過減小工件溫升來抑制工件熱變形進而獲取良好加工表面質(zhì)量的方法。

圖3 刀尖點附近溫度分布Fig.3 Temperature distribution near nose of tool

3 銑削參數(shù)對銑削力和銑削溫度的影響

3.1 銑削參數(shù)對銑削力的影響

圖5為銑削參數(shù)對銑削力和銑削溫度的影響曲線。vc=900m/min、ap=0.5mm、ae=6mm條件下，每齒進給量對銑削力和銑削溫度的影響曲線見圖5（a）。其中，x為刀具進給方向，y為切削寬度方向，可以看出，當每齒進給量小于0.07mm/z時，x方向切削力小于y方向；隨著每齒進給量的增大，x、y方向切削力呈增大趨勢，且x方向切削力大于y方向。這是由于隨著每齒進給量的增加，切削厚度增大，銑削變形抗力、切屑與刀具的摩擦阻力增大，故銑削力增大。fz=0.04mm/z、ap=0.5mm、ae=6mm條件下，銑削速度對銑削力和銑削溫度的影響曲線見圖5（b），其中，x和y向的切削力隨著切削速度的增大略微減小。這主要是因為：一方面，隨著切削速度的提高，剪切區(qū)域溫度升高，導致材料抗剪強度、摩擦系數(shù)和材料變形系數(shù)減小，進而使x和y向的切削力呈減小趨勢；另一方面，切削速度的提高使剪切角φ增大，導致剪切面積減小，變形減小，x和y向的切削力呈減小趨勢。同時，切削速度的提高導致高頻沖擊力逐漸增大，但在本文涉及的仿真切削速度范圍內(nèi)，單位切削面積切削力的減小程度大于高頻沖擊力的增大程度，所以x和y向的切削力呈減小趨勢。vc=900m/min、fz=0.04mm/z、ap=0.5mm條件下，銑削寬度對銑削力和溫度的影響曲線見圖5（c）。可以看出，x、y方向切削力隨著銑削寬度的增大而緩慢增大，這是因為當銑削寬度小于刀具半徑時，切削長度隨著銑削寬度的增加而增大，材料去除總體積增大，從而使工件變形產(chǎn)生的銑削變形抗力及切屑與刀具的摩擦阻力增大，故銑削力增大。vc=900m/min、fz=0.04mm/z、ae=6mm條件下，銑削深度對銑削力和銑削溫度的影響曲線見圖5（d）?？梢钥闯觯S著銑削深度的增大，x、y方向切削力線性增大，這是因為隨著銑削深度的增加，切削刃接觸長度增大，剪切面積增加，因而增大了切削力。

圖4 不同時刻工件和刀具上的溫度場分布Fig.4 Temperature field distribution on the tool and workpiece at different time

圖5 銑削參數(shù)對銑削力和銑削溫度的影響Fig.5 Effect of cutting parameters on cutting force and temperature

3.2 銑削參數(shù)對銑削溫度的影響

從圖5（a）中可以看出，隨著每齒進給量的增大，切削溫度變化范圍為210～422℃，銑削進給量的增加會使得溫度緩慢上升。這是由于進給量的增加會使得單位時間金屬切除量增加，切削熱增加，切削溫度升高；但與此同時切屑的變形系數(shù)隨著進給量增加會減小，單位體積材料去除量的切削功降低；此外隨著進給量增大，銑刀與切屑接觸長度增大，增大了切削熱量的傳出面積，切屑帶走更多的熱量，幾個因素相互影響，使得溫度增加。從圖5（b）中可以看出，切削溫度變化范圍為230～309℃，切削溫度與銑削速度成正相關性，這是由于隨著銑削速度提高，單位時間切除的金屬量越多，克服金屬變形與摩擦所需要的功也越大，產(chǎn)生切削熱也越多，故溫度上升。從圖5（c）中可以看出，切削溫度變化范圍為233～319℃，切削溫度隨著銑削寬度的增加而升高。原因在于銑削寬度增加時，去除的材料體積更大，導致克服切屑變形和摩擦消耗功增加，切削熱增加，進而切削溫度增大。從圖5（d）中可以看出，隨著銑削深度的增大切削溫度基本不變。雖然隨著銑削深度的增大，切削熱增大，但是切削刃切觸長度增加，提高了散熱條件，切削的接觸面積增大，散熱程度增大，因此銑削溫度變化不大。

4 結論

根據(jù)模擬鋁合金7055高速銑削過程中銑削力和銑削溫度的結果，得到以下結論：刀具上的最高溫度點位于距離刀尖大約0.05～0.1mm的前刀面上；每齒進給量增大，刀具進給方向上的銑削力顯著增大；銑削深度增大，銑削力顯著增大；銑削速度和銑削寬度對銑削力的影響不明顯；銑削溫度隨著每齒進給量、銑削速度、銑削寬度的增大而升高，銑削深度對切削溫度影響不大。

[1]FRIDLYANDER I N. Aluminum alloys in aircraft in the periods of 1970-2000 and 2001-2015[J]. Metal Science and Heat Treatment，2001,43(1)：6-10.

[2]DE LACALLE L N L，LAMIKIZ A，SáNCHEZ J A，et al.Cutting conditions and tool optimization in the high-speed milling of aluminium alloys[J]. Proceedings of the Institution of Mechanical Engineers， Part B ： Journal of Engineering Manufacture，2001，215(9)：1257-1269.

[3]JEONG D K，YOUN H K. High-speed machining of aluminium using diamond endmills[J]. International Journal of Machine Tools &Manufacture，1997,37(8)：1155-1165.

[4]楊巧鳳，張明賢，史興寬，等. 鋁合金高速銑削溫度變化規(guī)律試驗研究[J].西北工業(yè)大學學報，1998(3)：439-443.

YANG Qiaofeng，ZHANG Mingxian，SHI Xingkuan，et al. The critical speed for high-speed milling of aluminum alloy[J]. Journal of Northwestern Polytechnical University，1998(3)：439-443.

[5]NG E G，SZABLEWSKI D，DUMITRESCU M，et al. High speed face milling of a aluminium silicon alloy casting[J]. CIRP Annals -Manufacturing Technology，2004,53(1)：69-72.

[6]STRENKOWSKI J S，MOON K J. Finite element prediction of chip geometry and tool/workpiece temperature distributions in orthogonal metal cutting[J]. Journal of Engineering for Industry，1990,112(4)：313-318.

[7]XIE L J，SCHMIDTA J，SCHMIDTA C，et al. 2D FEM estimate of tool wear in turning operation[J]. Wear，2005,258(10)：1479-1490.

[8]成群林，柯映林，董輝躍. 航空鋁合金高速銑削加工的有限元模擬 [J].浙江大學學報（工學版），2006，40(1)：113-117

CHENG Qunlin，KE Yinglin，DONG Huiyue. Simulation of highspeed milling process of aerospace aluminum alloy[J]. Journal of Zhejiang University(Engineering Science)，2006，40(1)：113-117.

[9]路冬，李劍峰，融亦鳴，等. 航空鋁合金7075 - T7451三維銑削過程有限元仿真[J]. 中國機械工程，2008，19(22)： 2708-2710.

LU Dong，LI Jianfeng，RONG Yiming，et al. 3D milling process simulation of aero-aluminum alloy 7075-T7451[J]. China Mechanical Engineering，2008，19(22)： 2708-2710.

[10]SAFFAR R J，RAZFAR M R，ZAREI O，et al. Simulation of three-dimension cutting force and tool deflection in the end milling operation based on finite element method[J]. Simulation Modelling Practice and Theory，2008,16(10)： 1677-1688.

[11]趙云峰，張建華，邢棟梁，等. 鋁合金A12024銑削仿真分析與實驗研究[J].制造技術與機床，2011(1):101-103.

ZHAO Yunfeng，ZHANG Jianhua，XING Dongliang，et al. Simulation analysis and experimental study on milling of aluminum alloy A12024[J].Manufacturing Technology & Machine Tool，2011(1):101-103.

[12]成群林，柯映林，董輝躍. 航空鋁合金銑削加工中切削力的數(shù)值模擬研究[J].航空學報，2006，27(4):724-727.

CHENG Qunlin，KE Yinglin，DONG Huiyue. Numerical simulation study on milling force for aerospace aluminum alloy[J]. Acta Aeronautica et Astronautica Sinica，2006，27(4):724-727.

[13]畢運波，方強，董輝躍，等. 航空鋁合金高速銑削溫度場的三維有限元模擬及試驗研究[J].機械工程學報，2010，46(7):160-165.

BI Yunbo，F(xiàn)ANG Qiang，DONG Huiyue，et al. Research on 3D numerical simulation and experiment of cutting temperature for high speed milling of aerospace aluminum alloy[J]. Journal of Mechanical Engineering，2010，46(7):160-165.

[14]MARUSICH T D, ORTIZ M. Modelling and simulation of high-speed machining[J]. International Journal for Numerical Methods in Engineering, 1995, 38(21): 3675-3694.