劉光祖,李俊峰,楊延濤,曹學(xué)文,孫于晴
(西安航天動力機械有限公司,陜西 西安 710025)
旋壓是制造薄壁筒形件的有效工藝方法之一。旋壓所得的各類產(chǎn)品精度高,優(yōu)于沖壓加工精度;材料利用率高(最高可達80%),損耗低,經(jīng)濟效益顯著;力學(xué)性能優(yōu)良,經(jīng)過旋壓,材料的強度可提高35%~45%。旋壓加工廣泛應(yīng)用于航空、航天、軍工等領(lǐng)域[1-3]。
在旋壓加工中,通常將多個(兩個以上)旋輪沿工件環(huán)向均勻分布,相鄰旋輪之間沿軸向相互錯開,在垂直工件方向上又有不同的壓下量進行減薄,這種加工方式就是所謂的多旋輪錯距旋壓法。錯距旋壓生產(chǎn)效率高,加工過程中能有效避免工件震動,產(chǎn)品精度更高[4-6]。
7055是目前7×××系鋁合金中強度最高的鋁合金,具有比強度高,熱加工性能好,耐腐蝕和抗疲勞等優(yōu)點。用7055鋁合金制備的薄壁筒體,滿足輕質(zhì)高強的需求,具有廣闊的應(yīng)用前景[7-8]。本研究利用有限元模擬軟件,構(gòu)建仿真模型,模擬7055鋁合金筒形件錯距旋壓的加工方式,探究溫度、減薄率、進給比等工藝參數(shù)對成形規(guī)律的影響。
采用三旋輪錯距反旋壓模型,原理如圖1、圖2所示(圖中的t1、t2、t3為減薄量,C1、C2為錯距量)。將旋輪和芯模約束視為剛體;運動方式采用相對運動,即旋輪繞芯模旋轉(zhuǎn)并進給,同時依照線速度匹配的原則旋輪進行自轉(zhuǎn);為簡化模型,不考慮工件與芯模和空氣之間的熱量流動。在軟件中導(dǎo)入對應(yīng)溫度條件下的應(yīng)力-應(yīng)變數(shù)據(jù),進行動力學(xué)顯式模擬[9]。
圖1 反旋加工示意圖Fig.1 Schematic diagram of reverse spinning
圖2 三旋輪錯距旋壓示意圖Fig.2 Schematic diagram of staggered spinning of three spinning wheel
模型中工件的尺寸:內(nèi)直徑200 mm,壁厚10 mm,長300 mm;旋輪直徑200 mm,成形角25°,圓角半徑10 mm;芯模外直徑200 mm,長400 mm。裝配后如圖3所示。
圖3 仿真模型裝配圖Fig.3 Assembly drawing of simulation model
在仿真的準(zhǔn)備階段,選用鍛造退火態(tài)7055鋁合金進行熱壓縮實驗,獲得不同應(yīng)變速率和不同溫度條件下的真實應(yīng)力-應(yīng)變曲線,如圖4所示。
圖4 7055鋁合金不同應(yīng)變速率條件下真實應(yīng)力應(yīng)變曲線Fig.4 True stress-strain curve of 7055 aluminum alloy under different strain rate
7055鋁合金的密度為2 850 kg/m3,不同溫度下的楊氏模量和泊松比如表1所示。
表1 不同溫度條件下7055鋁合金楊氏模量和泊松比Table 1 Young′s modulus and Poisson′s ratio of 7055 aluminum alloy at different temperatures
本次模擬只進行一個道次,針對錯距旋壓的關(guān)鍵參數(shù)進行模擬計算分析,其中關(guān)鍵參數(shù)如表2所示。不同減薄率條件下各旋輪壓下量安排:減薄率20%(t1=1.0 mm,t2=0.5 mm,t3=0.5 mm);減薄率30%(t1=1.0 mm,t2=1.0 mm,t3=1.0 mm);減薄率40%(t1=2.0 mm,t2=1.0 mm,t3=1.0 mm)。其他參數(shù)固定不變,芯模轉(zhuǎn)速為6.283 r/s,旋輪間錯距量為C1=C2=4 mm,旋輪成形角為25°,旋輪和工件之間的摩擦因數(shù)為0.08,工件和芯模之間的摩擦因數(shù)為0.12。
表2 參數(shù)模擬方案Table 2 Parameter simulation scheme
鑒于鋁合金旋壓容易產(chǎn)生材料堆積、不貼膜等缺陷,將材料的堆積情況及貼模狀況作為評價旋壓效果的標(biāo)準(zhǔn),溫度、減薄率和進給比等都是控制材料堆積的關(guān)鍵參數(shù)。常見的材料堆積及不貼膜缺陷原因及預(yù)防措施如表3所示[3]。
表3 旋壓缺陷分析Table 3 Spinning defect analysis
工件和旋輪以及工件和芯模均設(shè)置為罰函數(shù)運動接觸,對模型進行網(wǎng)格劃分,采用C3D8R六面體結(jié)構(gòu)化網(wǎng)格類型,使用中性軸算法進行掃掠網(wǎng)格劃分,得到的網(wǎng)格更加規(guī)整,便于應(yīng)用ALE(任意拉格朗日-歐拉分析)網(wǎng)格自適應(yīng)技術(shù),在變形過程中網(wǎng)格自動調(diào)整形狀和尺寸,避免過度變形使模擬計算結(jié)果產(chǎn)生偏差[1,4]。工件的網(wǎng)格尺寸為3 mm,單元數(shù)量為87 600;芯模的網(wǎng)格尺寸為10 mm,單元數(shù)目為4 480;旋輪的網(wǎng)格尺寸為8 mm,單元數(shù)目為3 410。模型的網(wǎng)格劃分如圖5所示。由于使用了減縮積分的單元類型,為避免沙漏現(xiàn)象的產(chǎn)生,對工件網(wǎng)格進行沙漏控制[10-11]。
圖5 模型網(wǎng)格劃分Fig.5 Model mesh generation
本文對錯距反旋壓過程采用動態(tài)顯式算法模擬,模型網(wǎng)格數(shù)量較多,仿真過程引入質(zhì)量縮放因子為100 000,為了保證模型的可靠性,從能量角度進行驗證。常用的方法是構(gòu)建工件的動能/內(nèi)能(Kinetic energy/Internal energy,即ALLKE/ALLIE)比值變化曲線,比值低于5%時,模擬過程符合準(zhǔn)靜態(tài)條件,模型有效[6]。
如圖6所示,在模擬的開始階段,旋輪剛接觸工件,動能/內(nèi)能比值波動較大;隨著模擬的進行,動能/內(nèi)能比值逐漸穩(wěn)定在遠低于5%的水平,符合準(zhǔn)靜態(tài)模擬,模型是有效的。
圖6 工件的動能/內(nèi)能比值隨時間變化曲線Fig.6 Curve of kinetic energy/internal energy ratio of workpiece with time
鑒于模型的網(wǎng)格采用的是三維八節(jié)點減縮積分單元(C3D8R),模擬過程中可能出現(xiàn)沙漏現(xiàn)象,沙漏現(xiàn)象是指線性單元本身的積分點數(shù)比單元的節(jié)點數(shù)少,減縮積分單元在單元每個方向上的積分點數(shù)又會減少一個,所以極有可能出現(xiàn)沒有剛度的零能模式,也就是所謂的沙漏現(xiàn)象。評價沙漏現(xiàn)象對模擬準(zhǔn)確性的影響程度,采用的是偽應(yīng)變能/內(nèi)能(Artificial strain energy/Internal energy,即ALLAE/ALLIE)比值變化曲線,比值低于1%時,說明沙漏現(xiàn)象對模擬計算的影響可以忽略,計算結(jié)果沒有問題;比值在1%~5%之間,證明沙漏現(xiàn)象對模擬計算有影響,但計算結(jié)果仍可以接受;比值在5%~10%之間,說明沙漏現(xiàn)象對模擬計算影響較大,需要對模型做出調(diào)整,比如網(wǎng)格細(xì)化、沙漏控制等;比值大于10%時,說明模擬過程是不合理的,需要重新設(shè)置參數(shù)[12-14]。
輸出工件的偽應(yīng)變能/內(nèi)能比值如圖7所示??梢钥闯?,模擬過程中,偽應(yīng)變能/內(nèi)能比值迅速降低至5%以下,并處于穩(wěn)定狀態(tài),說明本文計算模型是可靠的。
圖7 工件的偽應(yīng)變能/內(nèi)能比值隨時間變化曲線Fig.7 Curve of artificial strain energy/internal energy ratio of workpiece with time
考慮進一步的網(wǎng)格細(xì)化和沙漏控制,比值還能降到更低水平,但是細(xì)化網(wǎng)格使網(wǎng)格數(shù)量成倍增加,穩(wěn)定時間增量急劇減小,計算效率會有所下降。不同尺寸網(wǎng)格的劃分效果如圖8所示。不同尺寸網(wǎng)格的穩(wěn)定時間增量如表4所示。
圖8 不同尺寸網(wǎng)格劃分效果Fig.8 Mesh division effect of different sizes
表4 不同尺寸網(wǎng)格的時間增量Table 4 Time increment of meshes with different sizes
模擬溫度對隆起缺陷的影響時,固定減薄率為30%,進給比為1 mm/r,模擬結(jié)果如圖9所示。由圖9可以看出,溫度越高,工件材料堆積的現(xiàn)象越明顯。分析其原因,可能是溫度對材料的軟化作用,材料流動性提升,位于旋輪前端的材料容易堆積而產(chǎn)生隆起;結(jié)合工件隆起高度變化曲線,實際旋壓時對該型鋁合金進行200 ℃左右的加熱,能有效降低旋壓時材料過度堆積而產(chǎn)生破壞的可能性。
圖9 工件隆起高度隨溫度變化曲線Fig.9 Variation curve of workpiece bulging height with temperature
模擬減薄率對旋壓加工的影響時,固定溫度為300 ℃,進給比為1 mm/r。工件隆起高度隨減薄率變化曲線如圖10所示。可以看出,隨著減薄率增加,工件材料的堆積隆起現(xiàn)象逐漸明顯。原因是,較大的減薄率意味著較大的變形量,材料的加工硬化導(dǎo)致反旋時材料向旋輪運動相反方向流動的阻力變大,材料大量堆積而產(chǎn)生隆起。鋁合金的極限減薄率能達到70%左右[1],但為了保證工件的尺寸精度和加工效率,實際生產(chǎn)中單道次旋壓的減薄率控制在30%~40%。
圖10 工件隆起高度隨減薄率變化曲線Fig.10 Variation curve of workpiece bulging height with thinning rate
研究進給比對旋壓的影響時,固定溫度條件為300 ℃,減薄率為30%。由隆起高度隨進給比變化曲線(圖11)可以發(fā)現(xiàn),大進給比導(dǎo)致材料隆起高度變大。其原因是大的進給比會使材料的環(huán)向流動分量減小,材料更多地沿軸向流動,工件易于堆積產(chǎn)生隆起。實際加工過程中,選擇合適的進給比,既能保證工件的加工精度,又能保證生產(chǎn)效率。
圖11 工件隆起高度隨進給比變化曲線Fig.11 Variation curve of workpiece bulging height with feed ratio
根據(jù)有限元模擬所得結(jié)論,選用江蘇豪然生產(chǎn)的噴射成形7055鋁合金筒形毛坯進行旋壓工藝實驗。旋壓采用三旋輪錯距熱反旋,加熱方式采用乙炔-氧火焰槍局部烘烤,全程使用激光測溫儀進行溫度監(jiān)測,減薄流程為10 mm→6 mm→4 mm→2 mm,預(yù)期壁厚2 mm。質(zhì)量指標(biāo):壁厚差為±0.20 mm,直線度為1.5 mm,圓度為2 mm。具體旋壓所用工藝參數(shù)如表5所示。各道次三個旋輪的壓下量如表6所示。
表5 旋壓工藝參數(shù)Table 5 Spinning process parameters
表6 各道次三個旋輪的壓下量Table 6 Reduction of three spinning wheels for each pass
首先對筒形毛坯進行精車,確保內(nèi)外表面光滑,起旋端加工出5 mm×30°倒角,方便旋壓加工時旋輪咬入,最終得到內(nèi)徑200 mm,長度300 mm,厚10 mm的旋壓毛坯。接下來安裝模具,模具外徑199 mm,并對模具進行80 ℃左右預(yù)熱;套裝毛坯并進行60 ℃左右預(yù)熱,達到溫度后開始旋壓;旋壓過程中對旋壓前區(qū)進行補熱,控制溫度在60 ℃~70 ℃,對成形區(qū)進行溫度測量,溫度為150 ℃左右。工藝流程如圖12所示。旋壓加工過程分三個道次進行,下一道次開始旋壓前都進行補充加熱,消除部分應(yīng)力,弱化加工硬化效應(yīng),最終得到旋壓成品如圖13所示。
圖12 加工工藝流程Fig.12 Processing flow
圖13 成品件Fig.13 Finished products
對所的成品件的加工精度進行檢測,外徑平均值為203.69 mm,壁厚平均值為1.94 mm,直線度為0.54 mm,圓度為0.36 mm。對比預(yù)期目標(biāo),壁厚差均為負(fù)值且低于0.20 mm的目標(biāo)值,平均壁厚差0.06 mm,直線度和圓度遠優(yōu)于指標(biāo)要求,整體成形精度較好。
針對7055鋁合金的多旋輪錯距旋壓進行有限元模擬,著重研究溫度、減薄率和進給比對材料堆積隆起的影響。熱反旋加工實驗所得成品件精度整體較好,研究結(jié)果對7055鋁合金的強力旋壓工藝有重要指導(dǎo)意義。
1)模型構(gòu)建思路能夠保證仿真計算的可靠性。動能/內(nèi)能比值處在極低的水平,偽應(yīng)變能/內(nèi)能比值能夠穩(wěn)定控制在5%以下,模型有效。
2)溫度、減薄率和進給比對旋壓加工質(zhì)量有重要影響。溫度升高,材料流動性好,容易出現(xiàn)堆積;隨著減薄率和進給比的增加,材料的隆起高度也增加。
3)7055鋁合金的熱反旋加工精度較好。有限元仿真及旋壓加工實驗結(jié)果對今后的工藝探索有重要參考意義。