大型飛機(jī)機(jī)身雙曲度蒙皮縱向拉伸成形加載優(yōu)化設(shè)計(jì)與試驗(yàn)研究

文松濤，曾斌，汪洋華，彭靜文，李衛(wèi)東，羅華

（1.江西洪都航空工業(yè)集團(tuán)有限責(zé)任公司，南昌330024;2.北京航空航天大學(xué)機(jī)械工程及自動(dòng)化學(xué)院，北京 100191）

蒙皮零件是構(gòu)成和維持飛機(jī)氣動(dòng)外形的外表零件。蒙皮零件的制造水平和產(chǎn)品質(zhì)量直接影響著飛機(jī)的氣動(dòng)性能和使用壽命。我國(guó)已開始自主研制大型民用客機(jī)，因此深入研究和發(fā)展先進(jìn)蒙皮數(shù)字化制造技術(shù)對(duì)于我國(guó)航空制造業(yè)具有十分重要的意義［1］。

大型飛機(jī)的蒙皮零件有著尺寸大、曲率小、剛性弱等特點(diǎn)，主要是通過(guò)縱向拉伸成形方式進(jìn)行成形。縱拉一般針對(duì)雙曲率或者縱向曲率大而橫向曲率小的零件?？v向拉形機(jī)的鉗口由一組子夾鉗組成，夾鉗間可以轉(zhuǎn)動(dòng)形成一定曲率的形狀，以成形雙曲度中的一個(gè)曲度，另外一個(gè)曲度在拉伸方向上形成。縱向不同于橫拉，其夾鉗的運(yùn)動(dòng)方式非常復(fù)雜，采用簡(jiǎn)單的通過(guò)模具截面線來(lái)確定加載軌跡的設(shè)計(jì)方法，通常不能得到能夠良好加載的軌跡，導(dǎo)致了需要通過(guò)大量的嘗試設(shè)計(jì)，來(lái)獲得一組比較好的加載參數(shù)，進(jìn)行生產(chǎn)試驗(yàn)。

近年來(lái)，為了滿足大型飛機(jī)的研制，我國(guó)各大主機(jī)廠分別引進(jìn)了數(shù)控縱向蒙皮拉伸機(jī)，但缺乏配套的工藝參數(shù)設(shè)計(jì)和優(yōu)化比較解決方案。現(xiàn)有工藝人員在縱向拉形方面經(jīng)驗(yàn)較少，靠傳統(tǒng)的以經(jīng)驗(yàn)分析和試驗(yàn)為主的飛機(jī)蒙皮成形方法，很難成形復(fù)雜蒙皮，材料和人工的浪費(fèi)非常大，數(shù)控拉形設(shè)備對(duì)變形的精確控制功能未能發(fā)揮，難以滿足大型飛機(jī)大尺寸蒙皮零件的生產(chǎn)要求，同時(shí)也嚴(yán)重制約了鈑金廠蒙皮制造水平的提高［2—3］。

本文針對(duì)典型的大型飛機(jī)雙曲率蒙皮零件，應(yīng)用蒙皮縱拉工藝設(shè)計(jì)與制造系統(tǒng)［4］，對(duì)夾鉗的加載軌跡進(jìn)行了優(yōu)化設(shè)計(jì);采用有限元模擬的方法，判斷設(shè)計(jì)軌跡是否合理;根據(jù)優(yōu)化后的加載軌跡，生成驅(qū)動(dòng)數(shù)控縱向拉形機(jī)運(yùn)動(dòng)的數(shù)控代碼進(jìn)行零件的成形;通過(guò)對(duì)成形后的零件進(jìn)行測(cè)量，獲得滿足零件交付要求的蒙皮零件。

1 典型蒙皮零件工藝分析模型的建立

圖1為某大型飛機(jī)機(jī)身蒙皮零件數(shù)模。該零件為雙曲率零件，除了通常具有的在拉伸方向垂直方向上的較大曲率外，拉伸方向的曲率也較明顯。另外，根據(jù)該零件形狀設(shè)計(jì)的模具如圖2所示，從中可清楚地看到該零件在右前方和左后方相比另一方有較大的傾斜度，該零件明顯不對(duì)稱。

圖1 雙曲率蒙皮零件Fig.1 The double-curved skin part

圖2 雙曲率零件縱拉模具Fig.2 The longitude stretch forming tool of double-curved skin part

零件材料為鋁合金2524-T3，厚度為1.6 mm。分別采用單向拉伸試驗(yàn)和FLD試驗(yàn)獲得材料性能參數(shù)和破裂準(zhǔn)則，另外分別采用控制拉形變形量和應(yīng)變速率的方法來(lái)避免滑移線和粗晶的出現(xiàn)。

根據(jù)零件數(shù)模信息，毛料為1.6 mm×1845 mm×4795 mm?？紤]到拉形實(shí)際過(guò)程中懸空段的長(zhǎng)度及包覆角對(duì)成形結(jié)果的影響，適當(dāng)延長(zhǎng)拉伸方向長(zhǎng)度，另外為使夾鉗邊緣子夾鉗有效夾持毛料，此處設(shè)置毛料尺寸為1.6 mm×2000 mm×6400 mm。設(shè)置輔助工作臺(tái)高為300 mm。在蒙皮縱拉工藝設(shè)計(jì)與制造系統(tǒng)（XSTRStretch）［4］中建立的零件的工藝分析模型如圖3所示。

圖3 蒙拉工藝分析模型Fig.3 Analysis model of stretch forming process

2 零件加載軌跡優(yōu)化設(shè)計(jì)

根據(jù)零件的工藝分析，以及鋁合金材料2524-T3的材料性能，首先設(shè)計(jì)蒙皮加載的延伸率為2.0%。為了保證零件成形能夠盡量包覆模具且不產(chǎn)生過(guò)多的局部拉伸量，按照以下原則設(shè)計(jì)夾鉗的包覆角度。即成形終了位置，毛料盡量沿模具兩端切線延伸，毛料剛好貼合模具且不產(chǎn)生過(guò)多的包覆。如圖4所示。

圖4 毛料與夾鉗的包覆Fig.4 Blank and wrapping of jaws

因?yàn)榱慵孑^為平緩，因此其成形過(guò)程如下:初始夾持→預(yù)拉伸→曲鉗口→補(bǔ)拉伸。按照以上的描述，設(shè)計(jì)蒙皮零件的加載軌跡參數(shù)如表1和表2所示。

表1 拉形工藝參數(shù)Table 1 Process parameters of stretch forming

表2 曲鉗口角度Table 2 Angles of curving jaws （°）

將上述加載軌跡參數(shù)輸入蒙皮拉形工藝設(shè)計(jì)與制造系統(tǒng)，輸出有限元前置處理文件，導(dǎo)入有限元軟件Pam-Stamp中進(jìn)行模擬，模擬結(jié)果如圖5所示。

從圖5可以看到，毛料在拉伸過(guò)程中，與模具貼合良好，毛料內(nèi)部區(qū)域沒有大變形的發(fā)生，因此加載工藝參數(shù)設(shè)計(jì)合理。但是在毛料與夾鉗的4個(gè)角上，即圖中橢圓位置，發(fā)生了主應(yīng)變?yōu)?.23的變形，此變形量已經(jīng)超過(guò)了2524-T3鋁合金材料的極限拉伸量。通過(guò)調(diào)整模具的高度，增加最終拉伸時(shí)夾鉗的仰角，可以減小毛料兩側(cè)邊與夾鉗交叉的位置的拉伸量。將輔助工作臺(tái)高度設(shè)置為500 mm，再采用上述拉形工藝參數(shù)進(jìn)行有限元模擬，模擬結(jié)果如圖6所示。

圖6 修改模具高度后的拉形模擬主應(yīng)變分布Fig.6 Distribution of major strain after increase the height of tool

從圖6可以看到，增加模具高度后，最終毛料的最大主應(yīng)變減少為0.072，已經(jīng)屬于2524-T3材料合理的變形范圍。由于模具高度與模具頂部貼模情況有一定關(guān)系，如果再增加模具高度，毛料邊緣與夾鉗交叉位置應(yīng)變會(huì)有減少，但會(huì)造成與模具頂部位置接觸的材料發(fā)生起皺，因此500 mm的輔助工作臺(tái)高度是最理想的高度。

3 典型零件生產(chǎn)性試驗(yàn)

按照優(yōu)化的加載軌跡，在數(shù)控拉形機(jī)上進(jìn)行1∶1的典型零件自動(dòng)拉形，驗(yàn)證加載軌跡設(shè)計(jì)的合理性與有限元模擬的準(zhǔn)確性，并積累工藝經(jīng)驗(yàn)和實(shí)驗(yàn)數(shù)據(jù)。

3.1 生產(chǎn)試驗(yàn)過(guò)程

在進(jìn)行生產(chǎn)性試驗(yàn)之前，需要對(duì)毛料進(jìn)行網(wǎng)格印制（如圖7所示），用于試驗(yàn)完成后，測(cè)量成形后的板料內(nèi)的局部變形量，與設(shè)計(jì)的延伸率和有限元模擬結(jié)果進(jìn)行對(duì)比。

圖7 毛料印制網(wǎng)格區(qū)域Fig.7 Areas of grid in blank

根據(jù)優(yōu)化后的零件加載軌跡，應(yīng)用蒙皮縱拉工藝設(shè)計(jì)與制造系統(tǒng)生成設(shè)備的控制代碼，按照數(shù)控代碼自動(dòng)執(zhí)行拉形過(guò)程。成形過(guò)程中的零件如圖8所示。

圖8 試驗(yàn)中的零件Fig.8 The part in stretch forming process

由圖8可以看出，在拉伸過(guò)程中，毛料與模具貼合情況良好。在模具兩端，毛料剛好沿模具兩端切線延伸，沒有產(chǎn)生過(guò)度包覆和不貼模的情況。

3.2 試驗(yàn)結(jié)果測(cè)量

成形結(jié)束后，松開夾鉗，將毛料卸載，通過(guò)測(cè)量零件中間和兩邊位置的毛料長(zhǎng)度，可以對(duì)比優(yōu)化設(shè)計(jì)的加載軌跡與實(shí)際軌跡伸長(zhǎng)量，驗(yàn)證加載軌跡設(shè)計(jì)與數(shù)控代碼的準(zhǔn)確性。測(cè)量結(jié)果見表3。

表3 拉伸長(zhǎng)度測(cè)量結(jié)果Table 3 The length measurement of stretched blank mm

由拉伸長(zhǎng)度測(cè)量情況可以看出（如圖9所示），有限元與加載軌跡設(shè)計(jì)的長(zhǎng)度是材料沒有發(fā)生回彈前的材料長(zhǎng)度，考慮到零件回彈，板料伸長(zhǎng)的計(jì)算精度滿足要求。

圖9 拉形板條試驗(yàn)結(jié)果對(duì)比Fig.9 The length comparison of stretched blank

在成形結(jié)束后，針對(duì)毛料上的印制的網(wǎng)格區(qū)域，采用網(wǎng)格應(yīng)變測(cè)量系統(tǒng)GMAS軟件測(cè)量應(yīng)變［5］，測(cè)量結(jié)果如圖10和11所示。

圖10 網(wǎng)格區(qū)域Fig.10 Areas of grid

圖11 網(wǎng)格測(cè)量結(jié)果Fig.11 Strain measurement

由網(wǎng)格測(cè)量結(jié)果可以看出，有限元模擬的結(jié)果與實(shí)際網(wǎng)格應(yīng)變測(cè)量結(jié)果一致，有限元的計(jì)算能夠?qū)?shí)際生產(chǎn)起到指導(dǎo)作用。

由圖12和圖13可以看出，零件最終成形切邊后，零件與模具的最大間隙約為3 mm，主要發(fā)生在模具兩端最高處和模具中間最低的兩側(cè)位置。零件的貼模情況理想，且與實(shí)際情況相符，滿足低應(yīng)力裝配的要求，不需要對(duì)模具型面進(jìn)行回彈修正。

圖12 試驗(yàn)有限元模擬貼模情況（mm）Fig.12 The fittability of FEM simulation

圖13 實(shí)際零件貼模情況Fig.13 The fittability of actual part

4 結(jié)論

大型飛機(jī)蒙皮零件雖然零件曲面平坦，但是有著尺寸大、曲率分布復(fù)雜等特點(diǎn)，在成形過(guò)程中會(huì)存在不貼模和拉伸量過(guò)大的情況。文中通過(guò)應(yīng)用蒙皮拉形工藝設(shè)計(jì)與制造系統(tǒng)，對(duì)典型蒙皮零件的加載軌跡進(jìn)行了優(yōu)化設(shè)計(jì)，采用有限元仿真對(duì)優(yōu)化結(jié)果進(jìn)行了評(píng)估;將優(yōu)化設(shè)計(jì)后的加載軌跡應(yīng)用于實(shí)際零件的成形，獲得了滿足質(zhì)量要求的零件。通過(guò)研究得到了以下結(jié)論。

1）有限元模擬具有較好的精確度。通過(guò)正交試驗(yàn)方法設(shè)計(jì)有限元模擬試驗(yàn)，可以方便快速地找到較好的工藝參數(shù)組合，實(shí)現(xiàn)工藝參數(shù)優(yōu)化。

2）對(duì)于大型飛機(jī)雙曲率蒙皮零件，模具的高度位置的確定，直接影響成形結(jié)果的好壞;若模具設(shè)置太低，會(huì)使得到毛料兩側(cè)拉伸量過(guò)大，直接造成毛料斷裂;若模具設(shè)置太高，會(huì)導(dǎo)致毛料中間與兩側(cè)產(chǎn)生不均勻的拉伸，導(dǎo)致毛料中間發(fā)生起趨。

3）一般情況下可以找到合適的拉形加載軌跡參數(shù)，使得回彈后的零件在小力作用下能夠貼模，滿足零件質(zhì)量要求，即可不用對(duì)拉形模具進(jìn)行回彈補(bǔ)償，從而降低設(shè)計(jì)和制造成本。

［1］張彥敏，周賢賓.飛機(jī)蒙皮拉伸成形工藝參數(shù)優(yōu)化［J］.航空學(xué)報(bào)，2006，27（6）:1203—1208.

ZHANG Yan-min，ZHOU Xian-bin.Parameter Optimization in AircraftSkin Stretch Forming Process［J］.Acta Aeronautica et Astronautica Sinica，2006，27（6）:1203—1208.

［2］韓金全，萬(wàn)敏，陳雪梅，等.飛機(jī)復(fù)雜蒙皮拉形工藝參數(shù)優(yōu)化設(shè)計(jì)［J］.塑性工程學(xué)報(bào)，2009，16（6）:96—101.

HAN Jin-quan，WAN Min，CHEN Xue-mei，et al.Optimization of Process Parameters in Stretch Forming of Aircraft Complex Skin［J］.Journal of Plasticity Engineering，2009，16（6）:96—101.

［3］張曉靜，周賢賓.板料成形回彈模擬［J］.塑性工程學(xué)報(bào)，1999，6（3）:56—62.

ZHANG Xiao-jing，ZHOU Xian-bin.Simulation of Springback in Sheet Metal Forming［J］.Journal of Plasticity Engineering，1999，6（3）:56—62.

［4］李衛(wèi)東.飛機(jī)蒙皮數(shù)字化拉形技術(shù)研究與系統(tǒng)開發(fā)［D］.北京:北京航空航天大學(xué)，2007.

LI Wei-dong.Technological Research and System Development of Numerical Stretch Forming for Aircraft Skin［D］.Beijing:Beihang University，2007.

［5］吳向東，萬(wàn)敏，李盛.便攜式板料應(yīng)變測(cè)量系統(tǒng)GMAS［J］，鍛壓技術(shù)，2005（z1）:208—210.

WU Xiang-dong，WAN Min，LI Sheng.Potable Sheet Metal Strain Measurement System GMAS［J］.Forging ＆Stamping Technology，2005（z1）:208—210.

［6］彭靜文，李衛(wèi)東，萬(wàn)敏.機(jī)翼前緣蒙皮拉形工藝參數(shù)優(yōu)化與試驗(yàn)［J］.塑性工程學(xué)報(bào)，2011，18（4）:74—78.

PENG Jing-wen，LI Wei-dong，WAN Min.Optimization and Experimental Research of Processing Parameters in Stretch Forming of Wing Leading Edge Skin［J］.Journal of Plasticity Engineering，2011，18（4）:74—78.

［7］韓金全.飛機(jī)蒙皮拉形模具型面與工藝參數(shù)優(yōu)化研究［D］.北京:北京航空航天大學(xué)，2009.

HAN Jing-quan.Design and Optimization of Tool Surface and Processing Parameters in Aircraft Stretching［D］.Beijing:Beihang University，2009.

［8］BAHLOUL R，MKADDEM A，DAL SANTO P，et al.Sheet Metal Bending Optimization UsingResponseSurface Method，Numerical Simulation and Design of Experiments［J］.International Journal of Mechanical Sciences，2006（48）:991—1003.

［9］JANSSON T，ANDERSSON A，NILSSON L.Optimization of Draw-in for an Automotive Sheet Metal Part An Evaluation Using Surrogate Models and ResponseSurfaces［J］.Journal of Materials Processing Technology，2005（159）:426—434.

［10］MKADDEMA A，BAHLOUL R.Experimental and Numerical Optimization of the Sheet Products Geometry Using Response Surface Methodology［J］.Journal of Materials Processing Technology，2007（189）:441—449.

［11］HE De-hua，LI Dong-sheng，LI Xiao-qiang.Optimization on Springback Reduction in Cold Stretch Forming of Titanum-alloy Aircraft Skin［J］.Transaction of Nonferrous Metals Society of China，2010（20）:2350—2357.

［12］HE De-hua，LI Xiao-qiang，LI Dong-sheng.Process Design for Multi-stage Stretch Forming of Aluminum Alloy Aircraft Skin［J］.Transaction of Nonferrous Metals Society of China，2010（20）:1053—1058.

［13］O'Donnell M.Finite Element Modeling of a Multi-stage Stretch-forming Operation Using Aerospace Alloys［D］.Belfast:University of Ulster，2003.

［14］HEUNG K，SEOK H.FEM-based Optimum Design of Multi-stage Deep Drawing Process of Molybdenum Sheet［J］.Journal of Materials Processing Technology，2007，184（1/3）:354—362.

［15］CAI Zhong-yi，WANG Shao-hui，XU Xu-dong.Numerical Simulation for the Multi-point Stretch Forming Process of Sheet Metal［J］.Journal of Materials Processing Technology，2009（209）:396—407.

［16］ESMAEILIZADEH Reza，KHALILI Kourosh，MOHAMMADSADEGHI Bagher，et al.Simulated and Experimental Investigation of Stretch Sheet Forming of Commercial AA1200 Aluminum Alloy［J］.Transactions of Nonferrous Metals Society of China，2014，24（2）:484—490.