卓德志， 謝延敏， 熊文誠， 唐 維， 黃仁勇

(西南交通大學(xué) 機械工程學(xué)院先進設(shè)計與制造技術(shù)研究所， 四川 成都 610031)

基于偏差調(diào)節(jié)法的蠕變時效成形回彈補償研究

卓德志， 謝延敏， 熊文誠， 唐 維， 黃仁勇

(西南交通大學(xué) 機械工程學(xué)院先進設(shè)計與制造技術(shù)研究所， 四川 成都 610031)

為解決蠕變時效成形后板料回彈較大的問題，基于偏差調(diào)節(jié)法，在不同補償方向建立了考慮水平方向偏差的模具型面修正方案.利用ABAQUS軟件建立了板料成形的有限元模型并對成形過程進行了模擬，并分別用偏差調(diào)節(jié)法和改進后的方法通過反復(fù)迭代對構(gòu)件進行回彈補償.結(jié)果表明，使用偏差調(diào)節(jié)法至少需要5次迭代才能使構(gòu)件滿足成形精度要求，而考慮水平偏差的方法僅需要3次，而且獲得的構(gòu)件的精度更高.研究表明考慮水平偏差的方法具有收斂速度快、成形精度高的特點.

時效成形; 回彈; 模面修正; 補償方向

蠕變時效成形是一種將材料的人工時效和蠕變成形相結(jié)合的成形方法.這種方法是利用材料在一定溫度下應(yīng)力隨時間減小而使部分彈性應(yīng)變轉(zhuǎn)化為塑性應(yīng)變的特性，實現(xiàn)零件的成形并同時完成時效處理，大大改善了合金內(nèi)部的微結(jié)構(gòu)狀態(tài)，提高了合金的強度[1].經(jīng)過蠕變時效成形后的構(gòu)件，往往具有殘余應(yīng)力小、產(chǎn)品精度高、材料機械性能好等特點.這種成形方法適用于可時效強化型合金的整體帶筋和變厚度大曲率復(fù)雜外形和結(jié)構(gòu)的整體壁板結(jié)構(gòu)的成形制造，其被認為是下一代大型民用飛機建造過程中特別重要的金屬成形工藝之一[2].但該方法在成形過程中彈性應(yīng)變并不能完全轉(zhuǎn)化為塑性應(yīng)變，因此板料成形過后會出現(xiàn)回彈.為提高成形精度，需要對模具型面進行修正.目前，以有限元數(shù)值模擬技術(shù)為基礎(chǔ)，蠕變時效成形的回彈預(yù)測和型面補償研究已經(jīng)廣泛開展.

目前主要的回彈補償方法有回彈向前法[3-4]、響應(yīng)面法[5-7]和偏差調(diào)節(jié)法[8-9]等.Gan等[8]提出的偏差調(diào)節(jié)法(displacement adjustment method)將板料成形結(jié)束后的回彈量反向補償?shù)匠跏寄＞呱希⒁源藰?gòu)造出新模具型面，這種方法收斂速度快、適用范圍廣；甘忠等[10-11]通過線性函數(shù)、小波變換閉環(huán)控制等方法實現(xiàn)了蠕變時效成形的回彈補償；劉大海等[12]應(yīng)用有限元軟件ABAQUS實現(xiàn)了筋板件的時效成形型面補償；黃霖等[13]用偏差調(diào)節(jié)法實現(xiàn)了在時效成形中的模面修正.但板料上同一位置在回彈前后存在水平方向的位移[14]，而偏差調(diào)節(jié)法僅將回彈量補償?shù)截Q直方向.因此本文在其基礎(chǔ)上，在蠕變時效成形的模面修正中考慮了水平方向的補償，并用有限元模型進行驗證.

1 考慮水平方向的補償方案

1.1 偏差調(diào)節(jié)法

圖1 偏差調(diào)節(jié)法示意圖Fig.1 Diagram of displacement adjustment method

偏差調(diào)節(jié)法由Wagoner提出，是一種基于有限元的迭代修正方法.其基本含義是在板料上取一定數(shù)量的節(jié)點，以節(jié)點為補償對象，通過向回彈反方向移動節(jié)點的位置來構(gòu)造新的模具型面，即：將成形回彈前后各節(jié)點的位移量加到回彈的反方向來確定新模具型面.在圖1中，以目標(biāo)形狀為模具，板料經(jīng)過成形、卸載后出現(xiàn)回彈，與目標(biāo)形狀間有偏差，記為Δy.在構(gòu)造新的模具時，將Δy加到節(jié)點的位置上，得到新的模具型面，然后板料用新模具進行成形.如果回彈后的形狀仍然不能滿足精度要求，則進行下一次迭代.

(1)

式中K為調(diào)整控制系數(shù)，其值一般取0.8～1.2，根據(jù)回彈率來確定.回彈量較大時，K取較大值，可提高回彈下降速度；回彈量較小時，K取較小值，可避免修正過度.

如圖2所示，設(shè)Y為目標(biāo)型面的拱高，回彈大小的評價方法為

(2)

式中SP為最大回彈量，當(dāng)SP≤ε時，則視為構(gòu)件滿足精度要求.本文中ε=5%.

圖2 板料成形前后狀態(tài)圖Fig.2 The state diagram of sheet before and after forming

1.2 考慮水平方向補償?shù)钠钫{(diào)節(jié)法

Wagnoer提出的偏差調(diào)節(jié)法將回彈量補償?shù)經(jīng)_壓成形的反方向，但未考慮回彈前后節(jié)點位置存在的水平方向偏移量.因此本文在Wagnoer的補償法基礎(chǔ)上，考慮在水平方向進行一定補償，并檢驗其效果.

圖3中，根據(jù)等弧長原則，o，o′，o″分別為構(gòu)件回彈前、回彈后及緊貼補償后模具時同一節(jié)點的位置；oa是豎直方向，也是傳統(tǒng)偏差調(diào)節(jié)法的補償方向；ob是oo′的連線方向；oc是o點所在弧線的法線方向.本文選用連線方向補償和法線方向補償這2種考慮水平方向補償?shù)姆桨?，分別對應(yīng)圖3中ob和oc方向.連線補償是將回彈量補償?shù)焦?jié)點回彈前后位置的連線方向,即ob方向；法向補償則是將回彈量補償?shù)交貜椙霸摴?jié)點的法線方向上，即oc方向，每次補償?shù)拇笮榘辶仙贤还?jié)點回彈前后的距離.由于ob，oc與豎直方向夾角不同，所以2種補償方案在水平方向的分量是不同的.

圖3 補償方向示意圖Fig.3 Diagram of compensation direction

2 研究對象及有限元模型

2.1 研究對象

選用2324鋁合金板作為研究對象，其泊松比為0.33，彈性模量為68 GPa.

蠕變本構(gòu)模型是反映金屬蠕變行為的表達模型.目前描述金屬蠕變行為的本構(gòu)模型一般分為冪指數(shù)函數(shù)和雙曲正弦函數(shù)兩種基本形式.Kowalewski等[15]提出了冪指數(shù)函數(shù)形式的基于機理的單向受力條件下蠕變本構(gòu)方程，其能夠更好地描述金屬蠕變行為的3個階段：瞬時蠕變階段、穩(wěn)態(tài)蠕變階段和加速蠕變階段.其中第3階段為金屬蠕變速率快速升高直到材料發(fā)生斷裂破壞.而蠕變時效成形過程中應(yīng)力小、時間短，一般認為時效成形沒有進入第3階段.因此Huang等[16]在不考慮第3階段的情況下，針對Kowalewski的蠕變損傷本構(gòu)方程提出了更簡化的本構(gòu)模型.其本構(gòu)形式如下：

(3)

通過標(biāo)準蠕變試驗獲得2324鋁合金在190 ℃下不同應(yīng)力水平的蠕變應(yīng)變數(shù)曲線.通過擬合，得到本構(gòu)模型的參數(shù)如表1[17].

表1 本構(gòu)模型參數(shù)

2.2 蠕變時效成形的有限元模型

在ABAQUS中建立板料和模具的有限元模型.板料采用可變形實體建模，尺寸為400 mm×80 mm×1.5 mm；模具的設(shè)計型面為半徑為800 mm的圓柱面，并以此作為初始模具型面，設(shè)置為離散剛體，如圖4.

圖4 板料和模具的有限元模型Fig.4 Finite element model of sheet and die

分析過程為：

1)加載.在板料上施加大小為1 MPa的均布載荷，使板料變形并貼緊模具表面，分析時間為1 s.

2)保壓.保持載荷不變，進行蠕變時效模擬.在此階段，材料內(nèi)部發(fā)生應(yīng)力松弛，彈性應(yīng)變轉(zhuǎn)化為塑性應(yīng)變，回彈量減小，分析時間為18 000 s (5 h).

3)卸載.將施加在板料上的載荷作用解除，使板料自由回彈.

使用Newton-Cotes積分方法，采用Fortran語言將2324鋁合金的蠕變本構(gòu)方程編寫成ABAQUS/Standard蠕變子程序(user subroutines， creep)，在提交分析任務(wù)時調(diào)用，實現(xiàn)對蠕變時效過程的模擬.

3 不同方案的模面修正

為對比不同補償方案的效果，本文中取調(diào)整控制系數(shù)K=1.

3.1 第1次成形模擬

以設(shè)計型面為初始模具進行第1次成形模擬.圖5為二分之一截面的分析結(jié)果，計算得出最大回彈量SP=41.8%,需要進行模具型面修正.

圖5 第1次成形前后板料狀態(tài)Fig.5 The states of sheet before and after the first forming

3.2 偏差調(diào)節(jié)法修正

使用偏差調(diào)節(jié)法進行模面修正，修正過程中僅將補償量補償?shù)截Q直方向.圖6為模面迭代修正的計算結(jié)果.可以看出經(jīng)過5次修正，板料的最大回彈量為3.21%，滿足精度要求.

圖6 豎直方向補償時的最大回彈量下降趨勢Fig.6 The downward trend of maximum SP by compensation in vertical direction

3.3 考慮水平方向補償?shù)钠钫{(diào)節(jié)法修正

使用考慮水平方向補償?shù)姆桨高M行補償.分別將第1次成形回彈前后各個節(jié)點在軟件中的坐標(biāo)輸入計算程序中，得出修正后模具型面的節(jié)點坐標(biāo)，構(gòu)造出下一次成形的模具型面，然后進行成形模擬.表2為使用第1次補償后，3種補償方向的最大回彈量，可以看出考慮了水平補償?shù)倪B線和法線方向補償使回彈量下降更快.

表2 第1次不同方向補償后板料最大回彈量

Table 2 The maximumSPof sheet after the first compensation in different direction

補償方向最大回彈量豎直方向24.83%連線方向19.81%法線方向12.17%

在第1次修正后繼續(xù)進行修正，直到成形后構(gòu)件最大回彈量在允許范圍內(nèi).圖7為豎直方向、連線方向及法線方向補償?shù)淖畲蠡貜椓肯陆第厔輬D.

圖7 不同方向補償?shù)淖畲蠡貜椓肯陆第厔軫ig.7 The downward trends of maximum SP by compensation in different directions

圖7表明：連線方向補償經(jīng)過4次修正后構(gòu)件就滿足精度要求；而法線方向補償只需要3次；第1次補償時，法線方向補償在水平方向的補償分量最大，豎直方向補償分量最小，經(jīng)過第1次補償后法線方向補償?shù)臉?gòu)件回彈量下降最為明顯.

表3表明，當(dāng)構(gòu)件經(jīng)過最少次數(shù)的補償均滿足精度要求時，使用考慮水平方向補償?shù)姆椒ǖ玫降臉?gòu)件最大回彈更小，精度更高.

表3 滿足精度要求后構(gòu)件的最大回彈量

Table 3 The maximumSPof components after meeting the precision requirements

補償方向最大回彈量豎直方向3.21%連線方向3.07%法線方向0.58%

4 結(jié) 論

本文將有限元方法、蠕變本構(gòu)方程、數(shù)值模擬等方法運用到了蠕變時效成形的數(shù)值模擬及型面修正中.

1)對偏差調(diào)節(jié)法進行了改進，考慮了水平方向補償在模具型面修正中的作用.

2)研究結(jié)果表明，考慮了水平方向補償?shù)姆椒ū葍H在豎直方向補償具有更快的收斂速度；而且在構(gòu)件滿足精度要求時，利用法向補償獲得的構(gòu)件的最大回彈量最小，精度較高.

3)3種補償方案的結(jié)果對比表明，水平方向補償對回彈補償效果有較大影響.

[1] 湛利華，陽凌.時效蠕變與時效應(yīng)力松弛行為轉(zhuǎn)換關(guān)系[J].塑性工程學(xué)報，2013，20(3):126-131． ZHAN Li-hua， YANG Ling. Research on conversion relationship between aging reep and aging stress relaxation[J]. Journal of Plasticity Engineering， 2013, 20(3):126-131．

[2] 曾元松，黃遐，黃碩.蠕變成形技術(shù)研究現(xiàn)狀及發(fā)展趨勢[J] .塑性工程學(xué)報，2008，15(3):1-8． ZENG Yuan-song， HUANG Xia， HUANG Shuo. The research situation and the developing tendency of creep age forming technology[J]. Journal of Plasticity Engin-eering， 2008, 15(3):1-8．

[3] KARAFILLIS A P， BOYCE M C. Tooling design in sheet metal forming using springback calculation[J]. International Journal of Mechanical Sciences， 1992，34 (2): 113-131．

[4] KARAFILLIS A P， BOYCE M C. Tooling design accommodating springback error[J]. Journal of Materials Processing Technology， 1992， 32(1/2): 499-508．

[5] GUO Y Q， GATI W， NACEUR H， et al. An efficient DKT rotation free shell element for springback simulation in sheet metal forming[J]. Computer & Structur-es， 2002， 80(27/30): 2299-2312．

[6] NACEUR H， GUO Y Q， BEN-ELECHI S. Response surface methodology for design of sheet forming parameters to control springback effects[J]. Computers and Structures， 2006， 84(26/27): 1651-1663．

[7] 阮光明. 基于響應(yīng)面法的DP800高強鋼沖壓回彈工藝參數(shù)優(yōu)化[J]. 鍛壓技術(shù)，2015， 40(4):40-44． RUAN Guang-ming. Technological parameters optimization of springback for DP800 high strength steel based on response surface method[J]. Forming & Stamping Technology， 2015， 40(4):40-44．

[8] GAN W， WAGONER R H. Die design method for sheet springback[J]. International Journal of Mechanical Sciences, 2004, 46(7):1097-1113．

[9] LINGBEEK R A， GAN W， WAGONER R H, et al. Theoretical verification of the displacement adjustment and springforward algorithms for springback compensation[J]. International Journal of Material Forming， 2008,1(3):159-168．

[10] 甘忠，莊振民，黃官平.基于傳遞函數(shù)的時效成形模具修正算法[J].塑性工程學(xué)報，2013，20(5):38-86． GAN Zhong， ZHUANG Zhen-min， HUANG Guan-ping. Springback correction algorithm of age forming experiments based on transfer function[J]. Journal of Plasticity Engineering， 2013， 20(5):38-86．

[11] 甘忠，馮爽，錢志偉,等.基于小波變換閉環(huán)控制的時效成形回彈補償[J].塑性工程學(xué)報，2015，22(1):111-118． GAN Zhong， FENG Shuang， QIAN Zhi-wei, et al. Age forming springback compensation using wavelet transform based closed loop control[J]. Journal of Plasticity Engineering， 2015, 22(1):111-118．

[12] 劉大海，黎俊初，常春.基于ABAQUS的筋板件時效成形型面回彈補償算法[J]. 塑性工程學(xué)報，2012，19(6):74-78. LIU Da-hai， LI Jun-chu， CHANG Chun. An ABAQUS-based springback compensation algorithm of die profile for creep age forming of ribweb parts[J]. Journal of Plasticity Engineering， 2012， 19(6): 74-78．

[13] 黃霖，萬敏.鋁合金厚板時效成形回彈補償算法[J].航空學(xué)報，2008，29(5):1406-1410. HUANG Lin， WAN Min. Compensation algorithm for springback in age forming for aluminum alloy thick plate[J]. Acta Aeronautica et Astronautica Sinica， 2008， 29(5):1406-1410．

[14] 陽湘安，阮峰.模面幾何修正的回彈補償方向分析[J].塑性工程學(xué)報，2010，17(2)2:6-10. YANG Xiang-an， RUAN Feng. Analysis on springback compensation direction for die-face adjustment[J]. Journal of Plasticity Engineering， 2010, 17(2): 6-10．

[15] KOWALEWSKI Z L， HAUHURST D R， DYSON B F. Mechansims-based creep constitutive equations for an aluminum alloy[J]. Journal of Strain Analysis， 1994， 29(4):309-316．

[16] HUANG L， WAN M， CHI C L， et al. FEM analysis of springback in age forming for aluminum alloy plate[J]. Chinese Journal of Aeronautics， 2007， 20(6):564-569．

[17] 黃碩，曾元松，黃遐.2324鋁合金蠕變時效成形有限元分析[J].塑性工程學(xué)報，2009，16(4):129-133． HUANG Shuo， ZENG Yuan-song， HUANG Xia. FE analysis of creep age forming for Aluminum alloy 2324[J]. Journal of Plasticity Engineering, 2009， 16(4):129-133.

Research on springback compensation of creep age forming based on displacement adjustment method

ZHUO De-zhi， XIE Yan-min， XIONG Wen-cheng， TANG Wei， HUANG Ren-yong

(Institute of Advanced Design and Manufacturing， School of Mechanical Engineering,Southwest Jiaotong University， Chengdu 610031， China)

To solve the large springback problems of sheet after the creep age forming， based on the displacement adjustment (DA) method， the die surface adjustment methods which considered horizontal direction displacement in different compensation directions were proposed. The finite element model of sheet forming was established by using the software ABAQUS， and the forming processes were simulated. The processes of springback compensation were carried out separately by using the original displacement adjustment method and the improved ones. The results showed that the forming parts could meet the precision requirement with 5 times of iteration by using the original method， while only 3 times were needed by using the improved methods and forming parts reached a higher precision. The advantages of improved methods which are fast convergence rate and high precision are confirmed．

creep age forming; springback; die surface adjustment; compensation direction

2016-02-24.

本刊網(wǎng)址·在線期刊：http://www.zjujournals.com/gcsjxb

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

卓德志(1991—)，男，山東滕州人，碩士生，從事金屬蠕變成形型面修正及回彈研究,E-mail：zhuodezhi@163.com. http://orcid.org//0000-0003-3248-7354 通信聯(lián)系人：謝延敏(1975—)，男，四川安岳人，副教授，碩士生導(dǎo)師，博士，從事先進塑性加工技術(shù)仿真和板料成形的穩(wěn)健設(shè)計等研究，E-mail:xie_yanmin@swjtu.edu.cn.

10.3785/j.issn. 1006-754X.2016.06.005

TG 146.21

A

1006-754X(2016)06-0553-05