航空長梁整體結(jié)構(gòu)件加工變形預(yù)測及校正分析

黃曉明,孫 杰,李劍峰(1.濱州學(xué)院 機電工程系,山東 濱州 56600; .山東大學(xué) 機械工程學(xué)院,濟南 50061)

航空長梁整體結(jié)構(gòu)件加工變形預(yù)測及校正分析

黃曉明1,2,孫 杰2,李劍峰2
(1.濱州學(xué)院 機電工程系,山東 濱州 256600; 2.山東大學(xué) 機械工程學(xué)院,濟南 250061)

航空整體結(jié)構(gòu)件的加工變形和變形校正是制約結(jié)構(gòu)件加工效率的瓶頸技術(shù),針對長梁零件加工后出現(xiàn)變形問題,構(gòu)建薄壁長梁整體結(jié)構(gòu)件加工變形有限元模型進行數(shù)值計算.研究表明,在毛坯初始應(yīng)力和加工應(yīng)力的作用下,長梁件出現(xiàn)彎曲扭轉(zhuǎn)復(fù)合變形,工件結(jié)構(gòu)剛度不對稱是導(dǎo)致工件彎曲變形的重要原因.以變形結(jié)果作為校正目標(biāo)模型,進行加載-卸載反彎和扭轉(zhuǎn)校正分析,對校正參數(shù)、校正載荷施加形式及校正效果等關(guān)鍵技術(shù)進行研究.根據(jù)長梁整體結(jié)構(gòu)件變形校正圖譜,當(dāng)在-Y方向施加9 050 N 4點反彎校正載荷,在垂直X順時針方向施加145×103N·m矯正扭矩時,可以得到理想的校正效果.

加工變形; 變形校正; 有限元數(shù)值模擬; 長梁整體結(jié)構(gòu)件

隨著航空制造工業(yè)的快速發(fā)展,為減輕飛機重量,提高飛機的各項機械性能,大型航空整體結(jié)構(gòu)件逐漸替代了傳統(tǒng)的螺栓連接和鉚接的飛機組合件.盡管整體薄壁件在航空工業(yè)中有著諸多的優(yōu)點,然而多年來,整體薄壁結(jié)構(gòu)件的加工變形問題一直困擾著航空工業(yè).即使采用高速加工方法加工大型梁框類整體結(jié)構(gòu)件,其變形問題也未能得到徹底解決.加工變形問題已經(jīng)成為幾乎所有大型飛機結(jié)構(gòu)件研制與生產(chǎn)的障礙.

國內(nèi)外各主要飛機制造公司都在探索有關(guān)整體結(jié)構(gòu)件制造的關(guān)鍵技術(shù),相關(guān)研究主要包括兩個方面:加工變形的機理和已產(chǎn)生加工變形零件的校正.如:美國的Nervi[1]建立了毛坯初始殘余應(yīng)力引起加工變形的數(shù)學(xué)預(yù)測模型,研究結(jié)構(gòu)件變形與初始毛坯應(yīng)力的分布狀態(tài)之間的關(guān)系;Scippa等[2]建立了薄壁件銑削加工變形的柔性預(yù)測模型,結(jié)合有限元分析方法,分析了不同切削速度對切削穩(wěn)定性的影響;浙江大學(xué)的董輝躍等[3]采用數(shù)值模擬的方法,研究了不同走刀路徑對工件殘余應(yīng)力分布變化以及加工變形的影響;南京航空航天大學(xué)武凱等[4]提出了參數(shù)化建立多框體模型的“砌墻法”,對某飛機192框結(jié)構(gòu)件進行了銑削加工模擬研究.

航空整體結(jié)構(gòu)件校正方法的研究一般是基于小變形校正原則的反變形校正[5],校正過程涉及復(fù)雜的力學(xué)過程,是多種非線性因素的綜合影響過程,包括幾何非線性、材料非線性等.山東大學(xué)的王中秋等[6]提出基于滾壓應(yīng)力調(diào)整來校正變形的理論和方法,該方法通過定量計算滾壓量,確定極限校正的載荷來保證變形校正的生產(chǎn)效率和安全.

長梁承重結(jié)構(gòu)件是航空整體結(jié)構(gòu)件中最常見的結(jié)構(gòu)之一,其變形情況往往是伸縮、彎曲、扭轉(zhuǎn)等多種形式的組合,并且與具體的尺寸和結(jié)構(gòu)特點相關(guān).由于零件變形過大,不得不對變形過大的零件進行校形.加工變形與校正是航空整體結(jié)構(gòu)件生產(chǎn)的一體化工作,目前尚缺乏基于加工變形結(jié)果的變形校正的整體研究.本文通過建立薄壁整體結(jié)構(gòu)件加工有限元仿真模型，并進行了計算,獲得變形部位的變形量.以此為基礎(chǔ),采用有限元仿真技術(shù),對相應(yīng)的校正參數(shù)進行模擬計算,計算結(jié)果將用于指導(dǎo)校形生產(chǎn).

1 基于應(yīng)力施加的加工變形有限元分析

1.1長梁整體結(jié)構(gòu)件

圖1為某飛機長梁件結(jié)構(gòu)示意圖,為非對稱雙面框薄壁梁,其結(jié)構(gòu)包含了槽、筋條、緣條、孔等整體結(jié)構(gòu)件的典型特征,研究該零件對長梁結(jié)構(gòu)件具有一定的代表意義.該長梁沿X方向長為1 785 mm,下緣條高60 mm,上緣條高40 mm,上下緣條高度不相等,整體呈梯形.工件整體Y方向?qū)挾葹?0 mm,壁板厚度為2 mm,正面有9個隔框,反面有7個半隔框.

圖1 飛機長梁整體件結(jié)構(gòu)圖Fig.1 Beam aerospace monolithic componentsschematic drawing

1.2長梁結(jié)構(gòu)件應(yīng)力分布規(guī)律

應(yīng)力的存在是工件加工變形的根本原因,其中應(yīng)力又分為毛坯初始應(yīng)力和工件薄層加工應(yīng)力兩部分.

本文所研究的長梁整體結(jié)構(gòu)件是由7050-T7451鋁合金預(yù)拉伸板材通過高速銑削制成.該材料在生產(chǎn)過程中要經(jīng)歷淬火、預(yù)拉伸、人工時效等溫度、組織轉(zhuǎn)變等復(fù)雜的非線性彈塑性應(yīng)變過程,內(nèi)部形成不均勻的應(yīng)力場.圖2為基于裂紋柔度法測試、計算獲得的7050-T7451鋁合金預(yù)拉伸板材殘余應(yīng)力分布圖[7].從圖中可以看出,鋁合金預(yù)拉伸板材軋制方向和橫向方向殘余應(yīng)力沿板厚的中性面對稱分布,呈拉壓周期震蕩“M”型分布.

圖2 7050-T7451鋁合金板材初始殘余應(yīng)力分布Fig.2 7050-T7451 aluminium alloy plate initialresidual stress distribution

軋制方向殘余應(yīng)力呈現(xiàn)拉應(yīng)力—壓應(yīng)力—拉應(yīng)力變化趨勢.殘余拉應(yīng)力的最大值約為26.43 MPa,在距板材的表面處;最大的殘余壓應(yīng)力約為-29.08 MPa,在距板中面約22 mm板厚處.橫向方向在板中心區(qū)域為拉應(yīng)力區(qū)域,最大值約為7.48 MPa,殘余壓應(yīng)力區(qū)域為表層至厚度為10 mm區(qū)域,最大值約為-14.57 MPa.長梁整體結(jié)構(gòu)件的X,Y和Z方向分別對應(yīng)毛坯材料的預(yù)拉伸方向、寬度方向和厚度方向.

工件加工表面殘余應(yīng)力的形成過程是一個復(fù)雜的高應(yīng)變率、塑性大變形過程.研究表明,加工應(yīng)力的大小、特性與銑削加工參數(shù)存在以下非線性關(guān)系[8]:

(1)

1.3加工變形有限元建模過程

7050-T7451航空鋁合金預(yù)拉伸毛坯板材的基本物理屬性參數(shù)如表1所示.

表1 7050-T7451基本物理屬性Tab.1 Mechanical properties of 7050-T7451

本文采用有限元軟件ABAQUS,對隔框梁因材料去除產(chǎn)生的變形進行仿真分析,將鋁合金預(yù)拉伸板材的初始應(yīng)力施加到有限元模型中.目前,初始應(yīng)力場施加一般采用如下處理方法:先將初始應(yīng)力場離散化,然后將離散后的應(yīng)力數(shù)值逐層施加到網(wǎng)格單元中,這樣處理會產(chǎn)生應(yīng)力場離散誤差.為了減小應(yīng)力場離散化導(dǎo)致的計算誤差,根據(jù)應(yīng)力分布特點,將毛坯初始應(yīng)力進行多峰高斯曲線擬合,擬合公式如式(2)所示.將初始應(yīng)力場與毛坯厚度之間的函數(shù)關(guān)系,通過子程序SIGINI施加到鋁合金預(yù)拉伸毛坯板材中,這樣保證了初始應(yīng)力的連續(xù)性,提高計算精度.

σx= 70.74-71.51e-0.003 93Z2-

(2)

采用“生死單元技術(shù)”來實現(xiàn)材料的去除過程.有限元軟件通過程序?qū)卧匠痰膭偠染仃?或其他量)乘以一個約為0的減縮因數(shù),來達到殺死單元的目的,表現(xiàn)為與失效單元相關(guān)的單元屬性(質(zhì)量、阻尼、比熱以及其他的量)全部設(shè)置為0.利用單元生死技術(shù)可以快速、有效模擬銑削加工過程材料的切除過程,避免了從材料變形的微觀角度來討論被切除單元的分離標(biāo)準(zhǔn).在ABAQUS中,將要去除的隔框內(nèi)的材料設(shè)置為一個單元集合,每去掉一個單元集合的材料,有限元將其作為一個載荷步進行一次計算,隨著模擬過程中工件材料的去除,殘余應(yīng)力逐層釋放.在去除材料后的隔框表面施加加工應(yīng)力,由于加工應(yīng)力分布在工件表面淺表層,為了“捕捉”該加工應(yīng)力,保證仿真精度,將加工表面網(wǎng)格細化,如圖3所示.

圖3 長梁結(jié)構(gòu)件網(wǎng)格剖分Fig.3 Beam aerospace monolithic componentsmesh generation

加工參數(shù)選擇為:v=628 m/min,fz=0.06 mm,ap=5 mm,ae=12 mm.

采用6點定位原理對工件進行定位約束,當(dāng)零件加工完成后,釋放約束,工件處于自由狀態(tài).有限元模型求解應(yīng)力平衡和變形應(yīng)滿足以下2個條件:① 零件的整體和任一元素在節(jié)點上都必須保持靜力平衡;② 變形協(xié)調(diào)條件,工件變形后必須保持連續(xù)性.

1.4變形預(yù)測結(jié)果與分析

圖4為采用“生死單元技術(shù)”去除材料后長梁結(jié)構(gòu)件的變形云圖.為了便于觀察,將變形量放大200倍.從模擬的結(jié)果來看,長梁主要表現(xiàn)為在Y方向和Z方向的彎曲變形.工件沿Y正方向兩端向上翹曲,Y方向變形量為0.664 mm.沿Y方向彎曲的同時,長梁結(jié)構(gòu)件兩端沿Z負方向下垂,中間部位沿Z正方向凸起變形.在彎曲的變形的同時,長梁薄壁工件出現(xiàn)垂直X方向截面的扭轉(zhuǎn)變形,工件截面最大扭轉(zhuǎn)角為2.47°.

長梁兩端壁板橫截面為梯形,梯形下底長度為60 mm,上端為40 mm,因此,下翼板區(qū)域剛度大于上翼板區(qū)域剛度,剛度非對稱分布特征導(dǎo)致了長梁沿Y方向兩端上翹,中間部位下凹.該復(fù)雜長梁為非對稱雙面框結(jié)構(gòu),其中正面框加強筋鏈接上下翼板,而反面框為非封閉半框結(jié)構(gòu),因此,正面框區(qū)域剛度大于反面框區(qū)域剛度,這就導(dǎo)致了長梁沿Z方向出現(xiàn)彎曲.

圖4 長梁結(jié)構(gòu)件加工變形云圖Fig.4 Beam aerospace monolithic componentsmachining deformation

2 變形校正有限元分析

根據(jù)上節(jié)所得到的長梁結(jié)構(gòu)件的變形預(yù)測結(jié)果,該零件變形的主要特點是彎曲和扭轉(zhuǎn)的組合變形,因此,需要針對這兩種變形方式進行反彎校正和扭轉(zhuǎn)校正.變形校正是彈塑性力學(xué)行為,7050-T7451鋁合金的塑性力學(xué)特性如圖5所示.

2.1彎曲矯正

變形校正中,需要計算外載荷作用下,工件產(chǎn)生彈塑性變形后總的變形位移量,以及卸載后的殘余位移量.其中載荷作用下工件總的變形位移量為校正壓力機壓頭的總行程,卸載后殘留變形量代表該載荷下可校正的工件變形量.為了達到校正效果,應(yīng)使卸載后的殘留變形量與工件的初始變形量一致.對于結(jié)構(gòu)復(fù)雜,無法獲得初始變形參數(shù)和校正參數(shù)之間直接對應(yīng)關(guān)系,采用二分法迭代進行校正參數(shù)計算.為了比較不同校正載荷量及校正方式對變形校正效果的影響,分別計算-Y方向不同載荷量作用下3點彎曲校正和4點彎曲校正的校形效果.

圖5 鋁合金7050-T7451應(yīng)力-應(yīng)變曲線Fig.5 7050-T7451 aluminium alloy stress-strain curve

壓彎校正過程中,整體結(jié)構(gòu)件采用簡支梁約束,3點壓彎采用最大變形處作為校正時的加力點.表2為3點壓彎校正作用下,不同載荷量與反彎壓下最大位移量及校正后殘留變形的數(shù)據(jù).

4點彎曲校正時,以長梁件三等分處作為校正加力點,反彎載荷平分施加在兩個反彎點位置.長梁兩端采用簡支梁約束方式.表3為4點壓彎校正作用下,不同載荷量與反彎壓下最大位移量及校正后殘留變形的數(shù)據(jù).

表2 3點壓彎載荷量和下壓量與校正變形量數(shù)值關(guān)系Tab.2 The relationship of three-point bending load and the correction value

表3 4點壓彎載荷量和下壓量與校正變形量數(shù)值關(guān)系Tab.3 The relationship of four-point bending load and the correction value

彎曲變形校正的載荷位移關(guān)系就是反彎下壓載荷量與下壓位移量之間的關(guān)系,在考察彎曲變形校正特點時,建立載荷量與位移量之間的對應(yīng)關(guān)系,繪制反彎矯正關(guān)系圖.圖6為根據(jù)表2和表3數(shù)值繪制的載荷與下壓量關(guān)系圖.

圖6 反彎校正載荷與下壓量關(guān)系圖Fig.6 The relationship of anti-bendingload and down shifting value

從圖6可以看出,在隔框梁件外加載荷和總下壓量之間的關(guān)系曲線中,兩者基本呈線性關(guān)系,即使工件材料進入了塑性階段,也無法明顯地反映出工件的屈服特點.在相同載荷條件下,4點壓彎校正下壓量小于3點壓彎校正下壓量,換言之,如果要達到相同的反彎效果,4點壓彎校正需要施加更大的載荷.

在彎曲變形校正中,更加關(guān)注的是工件的殘留變形,拋開工件材料的彈性效應(yīng),考察載荷量與殘留變形量之間的對應(yīng)關(guān)系,得到反彎矯正關(guān)系圖.圖7為校正載荷與殘留變形量之間的關(guān)系曲線.從圖中可以看出,校正載荷與殘留變形之間的非線性特征非常明顯.在3點壓彎校正過程中,當(dāng)外加載荷為6 750 N時,校正工件開始進入塑形變形區(qū),根據(jù)校正關(guān)系圖譜,對于本例彎曲變形,校正載荷約為7 300 N;在4點壓彎校正過程中,當(dāng)外加載荷為8 350 N時,校正工件開始進入塑形變形區(qū),根據(jù)校正關(guān)系圖譜,4點校正載荷約為9 050 N.

圖7 不同校正方式載荷量與殘留彎曲量關(guān)系曲線Fig.7 The relationship of anti-bending loadand residual bending value

為了比較不同的校正形式對工件校正精度的影響,采用上述兩種校正載荷,將工件的最大變形處校正到理想位置,比較沿長梁長度方向不同位置校正前(即初始變形量)、校正后偏離理想位置的情況,校正效果數(shù)值如表4所示.

表4 不同校形載荷作用的校形結(jié)果Tab.4 The correction resulets of three-point and four-point anti-bending

從表4中看出,盡管對于工件最大變形處(工件的中間部位)采用3點彎曲校正和4點彎曲校正均可使該區(qū)域的變形取得很好的校正效果,但是在離開中間部位,如X=535.5 mm和X=1 249.5 mm,4點彎曲的校正效果顯然優(yōu)于3點彎曲的校正效果.這主要是因為4點彎曲校正時,在兩集中載荷作用的中間區(qū)域形成純彎矩的作用效果,使工件獲得了較好的校正效果.而3點彎曲中,只是在集中載荷作用的中間區(qū)域產(chǎn)生顯著的塑性變形,該區(qū)域相當(dāng)于塑性鉸.在偏離該集中載荷作用位置處,產(chǎn)生繞這一塑性鉸的剛性轉(zhuǎn)動,形成剛性位移,因此,工件其他區(qū)域的校正變形沒有適應(yīng)工件初始變形的曲率特點,無法獲得理想校形效果.

2.2扭轉(zhuǎn)校正

扭轉(zhuǎn)變形校正采用一端約束6個自由度,另一端施加校正扭矩.當(dāng)校正扭矩增大到一定程度,工件產(chǎn)生塑性變形,外載扭矩卸掉以后,工件將產(chǎn)生不可恢復(fù)的塑性變形.表5為對于扭轉(zhuǎn)作用下,不同外載荷扭矩與最大扭轉(zhuǎn)角量及校正后殘留扭轉(zhuǎn)角的數(shù)據(jù).根據(jù)表5繪制總校正扭矩與總扭轉(zhuǎn)角關(guān)系曲線,如圖8所示.從圖8中可以看出,即使工件材料進入塑性階段,也無法反映出工件的屈服特點.

同理,拋開工件材料的彈性效應(yīng),考察外加扭矩載荷量與殘留扭轉(zhuǎn)變形量之間的關(guān)系,得到扭轉(zhuǎn)矯正關(guān)系圖譜.圖9為校正扭轉(zhuǎn)載荷與殘留扭轉(zhuǎn)變形量之間的關(guān)系曲線,其反映的校正載荷與殘留變形之間的非線性特征是非常明顯的.從圖中可以看出,對于1 785 mm的懸臂長度在扭轉(zhuǎn)校正過程中,當(dāng)外加載荷為100×103N·m時,校正工件開始進入塑形變形區(qū),根據(jù)校正關(guān)系圖譜,對于本例彎曲變形,校正載荷約為145×103N·m.

表5 不同扭轉(zhuǎn)載荷量與總扭轉(zhuǎn)角及殘留扭轉(zhuǎn)角數(shù)值關(guān)系Tab.5 The relationship of torsion load and the correction value

圖8 長梁結(jié)構(gòu)件扭矩與總扭轉(zhuǎn)角關(guān)系曲線Fig.8 The relationship of anti-torsionload and torsion angle

圖9 長梁結(jié)構(gòu)件扭矩與殘留扭轉(zhuǎn)角關(guān)系曲線Fig.9 The relationship of anti-torsion loadand residual torsion value

3 結(jié)語

本文以復(fù)雜航空整體長梁結(jié)構(gòu)件為例,綜合研究加工變形—校正過程.基于毛坯初始應(yīng)力和加工應(yīng)力的施加獲得工件加工變形規(guī)律,以加工變形的有限元數(shù)值作為校正目標(biāo),進行有限元加載-卸載校正分析.研究結(jié)果表明,零件在Y方向彎曲變形、沿Z方向的彎曲和垂直X方向的扭轉(zhuǎn)變形是該零件的主要變形形式,基于有限元數(shù)值仿真可獲得實際校正操作需要的校正參數(shù),如位移下壓量、反彎校正載荷等,具有一定的工程意義.

[1] NERVI S.A mathematical model for the estimation of effects of residual stresses in aluminum parts[D].Washington:Washington University,2005.

[2] SCIPPA A,GROSSI N,CAMPATELLI G.FEM based cutting velocity selection for thin walled part machining[J].Procedia CIRP,2014,14:287-292.

[3] DONG H Y,KE Y L.Study on machining deformation of aircraft monolithic component by FEM and experiment[J].Chinese Journal of Aeronautics,2006,19(3):247-254.

[4] 武凱,何寧,廖文和,等.薄壁腹板加工變形規(guī)律及其變形控制方案的研究[J].中國機械工程,2004,15(8):670-674.

WU K,HE N,LIAO W H,et al.Study on machining deformations and their control approaches of the thin-web in end milling[J].China Mechanical Engineering,2004,15(8):670-674.

[5] 孫杰,柯映林.隔框類航空整體結(jié)構(gòu)件變形校正關(guān)鍵技術(shù)研究[J].浙江大學(xué)學(xué)報(工學(xué)版),2004,38(3):351-356. SUN J,KE Y L.Research on key technologies for straightening pocketed aircraft frame parts[J].Journal of Zhejiang University(Engineering Science),2004,38(3):351-356.

[6] 王中秋,李劍峰,孫杰,等.航空整體結(jié)構(gòu)件側(cè)壁滾壓校正的有限元分析[J].中國機械工程,2009,20(5):612-616.

WANG Z Q,LI J F,SUN J,et al.FEM analysis of deformation correction by side-wall rolling of aircraft monolithic components[J].China Mechanical Engineering,2009,20(5):612-616.

[7] 周長安,孫杰,李衛(wèi)東,等.航空鋁合金7050-T7451毛坯初始殘余應(yīng)力測試[J].機械科學(xué)與技術(shù),2014,33(8):1251-1254.

ZHOU C A,SUN J,LI W D,et al.The measurement of original residual stress in aircraft aluminum alloy[J].Mechanical Science and Technology for Aerospace Engineering,2014,33(8):1251-1254.

[8] HUANG X M,SUN J,LI J F,er al.An experimental investigation of residual stresses in high-speed end milling 7050-T7451 aluminum alloy[J].Advances in Mechanical Engineering,2013(2):435-447.

Analysisofbeamaerospacemonolithiccomponentsmachiningdeformationpredictionandcorrection

HUANGXiaoming1,2,SUNJie2,LIJianfeng2
(1.Mechatronics Engineering Department,Binzhou University,Binzhou 256600, Shandong, China; 2.School of Mechanical Engineering,Shandong University,Ji’nan 250061, China)

The machining deformation of aerospace monolithic components is bottleneck problem for aviation manufacturing industry.We construct the beam machining deformation and distortion correction finite element model to study the problem.According to the finite element method(FEM)result,deformation characteristics were obtained.Studies show that the bending and torsion deformation is a universal phenomenon for the beam aerospace monolithic components under the coupling effect of blank original residual stresses and machining residual stress,and the asymmetry structure stiffness of the beam is an important reason leading to deformation.Deformation characteristics were used to evaluate the correcting effect.Correcting key technology such as determination of straightening parameters,anti-bend and anti-twist load were studied.For typical deformational beam,correcting characteristics and correcting rule were given.The ideal correction results were 9 050 N anti-bend load in -Ydirection and 145×103N·m anti-twist load were applying to the example.

machining deformation; straightening; finite element method(FEM)numerical simulation; beam aerospace monolithic components

TP 391.73

： A

： 1672-5581(2017)03-0210-06

國家自然科學(xué)基金資助項目(51275277)

黃曉明(1981—),男,副教授.博士.E-mail:hxm2552@163.com