基于ABAQUS的航空葉輪銑削變形機(jī)理研究*

李　林，蔡安江

(西安建筑科技大學(xué) 機(jī)電工程學(xué)院，西安　710055)

?

基于ABAQUS的航空葉輪銑削變形機(jī)理研究*

李林，蔡安江

(西安建筑科技大學(xué) 機(jī)電工程學(xué)院，西安710055)

摘要：高速銑削制造航空整體葉輪時(shí)，葉輪葉片的變形問題是影響葉輪加工精度的主要原因。而銑削過程中的銑削力、切削熱、殘余應(yīng)力則是導(dǎo)致葉片變形的直接因素。為了探究葉片的變形機(jī)理，基于葉片的幾何特征和加工工況，提出了將葉輪葉片簡(jiǎn)化為懸臂梁結(jié)構(gòu)的分析方法，同時(shí)采用有限元分析軟件ABAQUS建立反映葉片銑削過程高溫、高應(yīng)變率狀態(tài)的銑削模型，模擬分析不同銑削參數(shù)下整體立銑刀銑削航空鋁合金7075-T7451葉片的過程。葉片銑削模擬過程揭示了銑削速度、每齒進(jìn)給量、徑向銑削深度三個(gè)主要銑削參數(shù)對(duì)銑削力、切削熱、葉片表面殘余應(yīng)力的影響，為制造整體葉輪時(shí)加工參數(shù)的選取、加工變形及震顫控制提供依據(jù)。

關(guān)鍵詞：航空葉輪；銑削變形；殘余應(yīng)力；有限元分析

0引言

復(fù)雜航空整體結(jié)構(gòu)件中的典型工件為整體葉輪，在葉輪銑削過程中葉片的變形控制及獲得較高的加工精度是加工過程中考慮最多且最重要的兩大因素[1]。傳統(tǒng)的基于幾何形狀規(guī)劃的刀具軌跡、選取工藝參數(shù)的方法，無法前瞻性的預(yù)測(cè)葉片的變形必然導(dǎo)致實(shí)際加工表面輪廓與理論輪廓之間存在較大偏差，難以保證葉片加工的精度和效率[2]。為了能從根本上解決這一問題，必須探究葉片的變形機(jī)理，從而選取合理的加工參數(shù)及加工補(bǔ)償。分析葉輪的幾何結(jié)構(gòu)可知，葉片根部與輪轂固定、頂端處無約束，與懸臂梁結(jié)構(gòu)類似，故將葉片簡(jiǎn)化為懸臂梁結(jié)構(gòu)。同時(shí)將簡(jiǎn)化葉片導(dǎo)入ABAQUS有限元分析軟件建立了整體葉輪葉片的銑削模型，通過不同銑削參數(shù)組合，揭示了葉片銑削過程中銑削力、切削熱及殘余應(yīng)力的變化規(guī)律，進(jìn)而為后序優(yōu)化關(guān)鍵的工藝參數(shù)及葉片加工變形補(bǔ)償提供了依據(jù)，最終達(dá)到改善葉片加工精度和提高整體葉輪加工效率的目的。

1有限元模擬關(guān)鍵技術(shù)處理與模型定義

有限元模擬的關(guān)鍵技術(shù)包括Johnson-Cook塑性本構(gòu)模型、材料模型的設(shè)置、切屑與工件的分離準(zhǔn)則、工件與刀具的接觸設(shè)置等。上述技術(shù)的解決是建立有限元銑削仿真模擬的關(guān)鍵所在，同時(shí)對(duì)于模擬結(jié)果是否準(zhǔn)確有著直接的影響[3]。

1.1塑性材料本構(gòu)模型

材料本構(gòu)模型一般表示為流動(dòng)應(yīng)力與應(yīng)變、應(yīng)變率、溫度等變形參數(shù)之間的數(shù)學(xué)函數(shù)關(guān)系。當(dāng)前常用的塑性材料本構(gòu)模型主要有：Bodner-Paton、Follansbee-Kocks、Johnson-Cook、 Zerrilli-Armstrong等模型[4]，Johnson—Cook模型認(rèn)為材料在高應(yīng)變速率下表現(xiàn)為應(yīng)變硬化、應(yīng)變速率硬化和熱軟化效應(yīng)，符合葉片銑削過程中的工件狀態(tài)[5]，故選擇Johnson-Cook模型作為材料本構(gòu)模型, 同時(shí)采用Johnson-Cook剪切失效準(zhǔn)則作為刀屑分離的準(zhǔn)則。

(1)

1.2刀具和切屑接觸摩擦模型

在粘結(jié)的條件下，刀屑接觸面的摩擦狀態(tài)表征為切屑與刀具的粘結(jié)層與其上層金屬的內(nèi)摩擦發(fā)生的剪切滑移，與材料的粘結(jié)面積大小以及材料的流動(dòng)應(yīng)力有關(guān)[6]。滑移區(qū)的摩擦力與接觸壓力成正比，切屑在滑移區(qū)滑動(dòng)最終離開前刀面[7]。如圖1所示,二維模擬切削過程中的摩擦類型及分布情況。

圖1　二維模擬切削摩擦類型及分布

利用庫倫摩擦定律定義刀具與切屑間的摩擦，通過下式來判斷摩擦的類型[8]

(2)

1.3三維銑削幾何模型建立與網(wǎng)格劃分

ABAQUS軟件的三維建模能力有限，本文在UG中建立銑刀的三維模型，轉(zhuǎn)換為中性IGS格式后導(dǎo)入ABAQUS的part模塊中進(jìn)行模型修正。銑刀直徑為12 mm，螺旋角為30o，銑刀刃數(shù)為4，工件的尺寸為300mm×50mm×5mm。本文采用自適應(yīng)網(wǎng)格劃分技術(shù)選定工件單元類型為C3D8R(六面體8節(jié)點(diǎn)三維實(shí)體縮減單元)，網(wǎng)格尺寸為0.1 mm，刀具單元設(shè)定為三角形單元Tri，單元類型為R3D3,網(wǎng)格尺寸為0.5mm。

1.4材料參數(shù)與運(yùn)動(dòng)參數(shù)設(shè)置

工件材料為航空鋁合金7050-T7451，密度ρ=2820kg/m3,彈性模量E=71.7GPa，泊松比μ=0.33, 轉(zhuǎn)換到Johnson-Cook 模型材料參數(shù)分別為A=490MPa、B=206.9MPa、n=0.344、c=0.005、m=1.80[9]。刀具的材料為硬質(zhì)合金，密度ρ=15000kg/m3，彈性模量E=210GPa，泊松比μ=0.22，線膨脹系數(shù)al=0.22mm/℃,導(dǎo)熱率K=4.6E+01W/(m·K)。

運(yùn)動(dòng)參數(shù)的設(shè)置是在載荷模塊中固定工件下底面的六個(gè)自由度(U1=U2=U3=UR1=UR2=UR3=0)，在局部坐標(biāo)系下，選中銑刀中心參考點(diǎn)，進(jìn)行進(jìn)給速度(V1—X軸負(fù)方向)數(shù)值設(shè)定和旋轉(zhuǎn)速度(VR3—繞Z軸旋轉(zhuǎn))數(shù)值設(shè)定，完成刀具沿著X軸負(fù)方向進(jìn)給線速度以及繞Z軸旋轉(zhuǎn)的角速度的運(yùn)動(dòng)軌跡定義。

各項(xiàng)參數(shù)設(shè)定完成后，提交Job命令進(jìn)行葉片銑削過程模擬，如圖2所示。圖中所示的是銑刀沿x軸方向行至100mm處時(shí)的銑削狀態(tài)。

圖2　葉片銑削過程模擬

2銑削加工過程分析

在實(shí)際加工過程中高速銑削相較于常規(guī)加工具有提高加工效率、改善工件的加工精度、優(yōu)化整體結(jié)構(gòu)零件加工表面質(zhì)量的優(yōu)點(diǎn)[10]，故將銑削速度定義在高速銑削狀態(tài)，揭示葉片變形規(guī)律，在DMU210FD銑車復(fù)合加工中心進(jìn)行葉輪的加工試驗(yàn)。選定的銑削速度為v=400、450、500m/min，fz=0.05、0.08、0.10、0.12mm/z，ae=0.5、1.0、2.0、3.0mm銑削深度ap=10mm，順銑加工。

圖3　銑削加工過程

3計(jì)算結(jié)果分析

圖4為v=400m/min，fz=0.1mm/z，ae=3 mm時(shí)得到的銑削力，可以看出銑削力呈周期性分布，這是由于銑刀有四個(gè)刀齒，所以當(dāng)?shù)毒哌\(yùn)動(dòng)一定時(shí)間后,切削力出現(xiàn)有規(guī)律性的波峰和低谷現(xiàn)象。將銑削力的系列波峰和波谷的最大值求和取平均值，通過對(duì)比發(fā)現(xiàn)當(dāng)徑向深度為1mm時(shí)隨每齒進(jìn)給量的增加銑削力沒有太大變化，當(dāng)徑向深度增加為2mm和3mm時(shí)，銑削力隨每齒進(jìn)給量的增加明顯加大，這是由于徑向銑削深度的增加，使瞬時(shí)參與切削的刀刃數(shù)開始多于1。此外在其他銑削參數(shù)不變的情況下，隨著銑削速度的增大，銑削力呈減小的趨勢(shì)。所以在相同的參數(shù)條件下，應(yīng)盡可能選擇高的銑削速度，以提高銑削效率，減小銑削力。

圖4　銑削力隨時(shí)間變化曲線

為了分析切削熱與銑削參數(shù)之間的關(guān)系，在此分別改變銑削速度、每齒進(jìn)給量、徑向深的大小進(jìn)行仿真分析結(jié)果如圖5所示。從銑削模擬過程來看，銑刀刀刃作用位置附近的工件表面溫度呈非均勻分布狀態(tài)。刀尖處最高溫度均隨銑削速度增加而升高，但溫度升高的速率是逐漸減小的。發(fā)生這一現(xiàn)象的原因是銑削速度的增大使得在單位時(shí)間內(nèi)產(chǎn)生的熱量增多，致使切削溫度升高，同時(shí)切屑產(chǎn)生的速度加快，切屑帶走的熱量增多，導(dǎo)致傳遞給銑刀切削熱的比例減小，從而使得刀尖溫度升高的速率減小。另外刀尖處最高溫度隨每齒進(jìn)給量、徑向切深的增加也有升高的趨勢(shì)，但是升高的幅度較小，這是由于每齒進(jìn)給量或徑向切深的增加均使得單位時(shí)間內(nèi)材料的切除率增大，導(dǎo)致切屑帶走的熱量也增加，故對(duì)刀具溫度升高幅度影響較小。

圖5　不同銑削參數(shù)對(duì)銑刀溫度的影響

在銑削加工中，殘余應(yīng)力的產(chǎn)生與切削力熱載荷和材料內(nèi)部微觀結(jié)構(gòu)等因素有很大關(guān)系[11]。圖6是在不同的銑削參數(shù)下，工件表面殘余應(yīng)力分布規(guī)律。

圖6　表面殘余應(yīng)力變化曲線

圖6a中是在銑削速度為450m/min、切削深度為1mm時(shí)，不同的每齒進(jìn)給量下殘余應(yīng)力的曲線圖。可以看出: 隨著每齒進(jìn)給量的增大，出現(xiàn)最大殘余壓應(yīng)力的深度由0.2mm增大到0.7mm，由壓應(yīng)力轉(zhuǎn)換為拉應(yīng)力的深度即塑性變形深度由0.56 mm增大到1.90mm。當(dāng)每齒進(jìn)給量為0.08mm/z時(shí)整體表面殘余應(yīng)力最小，每齒進(jìn)給量偏大和偏小都會(huì)使兩者增大。進(jìn)給量增大，工件塑性變形區(qū)域增大從而導(dǎo)致產(chǎn)生的表面殘余應(yīng)力減小，同時(shí)銑削溫度隨之增大，熱應(yīng)力引起的表面殘余應(yīng)力增大。進(jìn)給量減小，塑形變形區(qū)域會(huì)減小，產(chǎn)生的殘余應(yīng)力變大，銑削溫度降低，故綜合導(dǎo)致了圖中所示的情況。

圖6b中是在銑削速度為450m/min、每齒進(jìn)給量為0.1mm時(shí)，不同的徑向銑削深度下殘余應(yīng)力的曲線圖。可以看出：隨著徑向銑削深度的增加，工件殘余壓應(yīng)力的深度也隨之增大，工件發(fā)生塑性變形的深度加深。產(chǎn)生殘余壓應(yīng)力增大這一現(xiàn)象的原因是: 由于徑向銑削深度增大從而增加了銑削區(qū)域，增大了切削過程中的切削力，刀具對(duì)工件的作用力增大，導(dǎo)致工件殘余壓應(yīng)力的范圍變大。

圖6c是在每齒進(jìn)給量為0.1mm、切削深度為1mm時(shí)，不同的銑削速度下殘余應(yīng)力的曲線圖。可以看出工件表面的殘余應(yīng)力隨著銑削速度的增加有減小的趨勢(shì)，其原因在于隨著銑削速度增加，銑削力增大導(dǎo)致銑刀刀刃后方工件上的塑性變形增大，壓應(yīng)力增大，拉應(yīng)力減小; 而同時(shí)銑削溫度升高，熱應(yīng)力使表層產(chǎn)

生殘余拉應(yīng)力，拉應(yīng)力增大，壓應(yīng)力減小。上述兩方面原因疊加導(dǎo)致表面殘余應(yīng)力的變化。

4結(jié)束語

利用ABAQUS有限元分析軟件，通過針對(duì)不同銑削速度、每齒進(jìn)給量和徑向切削深度組合的模擬實(shí)例，得到了簡(jiǎn)化葉片模型在高速銑削過程中表面的殘余應(yīng)力在進(jìn)給方向上沿深度的分布規(guī)律以及銑削力和切削熱的變化情況。為進(jìn)一步研究葉片加工變形補(bǔ)償和葉片加工震顫問題提供了依據(jù)。

[參考文獻(xiàn)]

[1] 湯愛君.薄壁零件銑削加工三維穩(wěn)定性及參數(shù)優(yōu)化[M].北京：國(guó)防工業(yè)出版社，2014.

[2] 程鳳軍，陳國(guó)定，王濤.Deform3D軟件支持的球頭銑刀銑削仿真分析[J].機(jī)械科學(xué)與技術(shù)，2012，31(3)：412-416.

[3] 姬偉. 薄壁件銑削加工誤差預(yù)測(cè)與快速銑削仿真平臺(tái)的開發(fā)[D].南昌：南昌航空大學(xué)，2012.

[4] tasi J S,Liao C L. Finite-element modeling of staic surface errors in the peripheral milling of thin walled workpieces[J]. Journal of Materials Processing Technology,1999,94:235-246.

[5] Ratchev S, Nikov S, Moualek I. Material removal simulation of peripheral milling of thin wall low-rigidity structures using FEA[J].Adv.Eng.Soft,2004,35(8-9):481-491.

[6] 金浩.高速切削銑削力的有限元分析與試驗(yàn)研究[D].西安：西安理工大學(xué)，2014.

[7] Zorev N N. Inter-relationship between shear processes occurring along tool face and shear plane in metal cutting[A].International Research in Production Engineering,1963.

[8] 江丙云.ABAQUS工程實(shí)例詳解[M].北京：人民郵電出版社，2014.

[9] 中國(guó)航空材料手冊(cè)編輯委員會(huì).中國(guó)航空材料手冊(cè)(第四卷)[M]. 2版.北京:中國(guó)標(biāo)準(zhǔn)出版社，2002.

[10] Tlusty J, Smith S,Winfough W R. Techniques for the use of long slender end mills in high speed milling[J].Annal of the CIRP,1996,45(1):393-396.

[11] 曾紅，韓笑.銑削加工曲面殘余應(yīng)力有限元分析[J].機(jī)床與液壓，2013，41(5)：49-52.

(編輯趙蓉)

Research on the Deformation Mechanism of the Aviation Impeller Milling Based on ABAQUS

LI Lin，CAI An-jiang

(School of Mechanical and Electrical Engineering，Xi’an University of Architecture and Technology，Xi ’an 710055，China)

Abstract：The deformation of the impeller blade is the main reason for the machining accuracy of the impeller when high speed milling the whole impeller. The milling process of milling force, cutting heat, residual stress is a direct factor leading to blade deformation. In order to explore the mechanism of deformation of blade, based on the geometrical characteristics of the blade and the processing condition, puts forward a design of simplifying the impeller blades to the cantilever beam structure。Using finite element analysis software ABAQUS to establish a model of the milling blade to reflect the state at high temperature and high strain rate. The simulation of Al 7075 -T7451 blade milling process is under different milling parameters and by overall vertical milling-tool. This milling machining simulation process reveals the influence of the milling speed, feed per tooth, radial milling depth, three main processing parameters ,on milling force, cutting heat and residual stress of the blade. The research can provide helps in the selection of milling parameters, machining deformation and controlling milling vibration.

Key words：aviation impeller; milling deformation; residual stress; finite element analysis

文章編號(hào)：1001-2265(2016)05-0036-03

DOI:10.13462/j.cnki.mmtamt.2016.05.010

收稿日期：2015-06-04；修回日期：2015-07-10

*基金項(xiàng)目：陜西省工業(yè)攻關(guān)項(xiàng)目(2012K09-24、2014K07-06)；西安市科技計(jì)劃項(xiàng)目(CXY1439(3))

作者簡(jiǎn)介：李林(1990—)，男，山東棗莊人，西安建筑科技大學(xué)碩士研究生，研究方向?yàn)槎噍S銑削加工工藝，(E-mail) willsendme@163.com；蔡安江(1965—)，男，上海人，西安建筑科技大學(xué)教授，博士生導(dǎo)師，研究方向?yàn)閿?shù)字化設(shè)計(jì)制造技術(shù)。

中圖分類號(hào)：TH166;TG659

文獻(xiàn)標(biāo)識(shí)碼：A