王大勇，時陽，黃永紅

(1.大連交通大學材料科學與工程學院，遼寧大連116028;2.中國北車集團大連機車車輛有限公司工藝部，遼寧大連116028)*

Q235的L截面型鋼壓彎的有限元分析及回彈預測

王大勇1，時陽1，黃永紅2

(1.大連交通大學材料科學與工程學院，遼寧大連116028;2.中國北車集團大連機車車輛有限公司工藝部，遼寧大連116028)*

分析發(fā)電機風道前托架的彎曲工藝，引入相對彎曲半徑系數(shù)C，對傳統(tǒng)板材壓彎回彈預測公式進行修正，得到計算L截面型鋼回彈的公式.利用有限元軟件DYNAFORM對Q235材料的L截面型鋼壓彎模擬，對回彈進行預測，結合試驗反推出C值，總結出C值和相對彎曲半徑r0/H關系曲線.利用該曲線可以計算出所需的凸模圓角半徑值，優(yōu)化了模具設計和壓彎工藝，為實際生產中的壓彎模的設計提供參考，減少了試模次數(shù)，降低生產成本.

回彈;壓彎L截面型鋼;有限元;相對彎曲半徑系數(shù)C

0 引言

回彈是板材沖壓成形過程中產生的一種缺陷，回彈的存在是不可避免的，影響回彈的因素多種多樣，關于回彈問題人們也做了大量工作，主要是通過對回彈的預測和控制兩方面進行研究，首先是對回彈進行預測，方法主要包括理論解析法，數(shù)值模擬法，實驗法，然后通過改善成形工藝和模具補償法來實現(xiàn)對加工零件回彈的控制，近年來，板材回彈的研究逐漸成熟［1－7］.

回彈缺陷同樣影響著型材的加工成形，但相對于板材型材回彈預測缺少準確的數(shù)據(jù)做指導，傳統(tǒng)的沖壓手冊中關于回彈的經(jīng)驗公式及模具設計數(shù)據(jù)也主要針對平面板材，應用于型材有很大的誤差，使得型材彎曲模具設計費時費力，嚴重影響生產效率.壓彎成形作為一種經(jīng)濟，高效的加工方法在實際的生產中得到普遍應用，但是其精度較低，而且關于型材壓彎回彈的資料數(shù)量更是有限.因此，參考板材回彈研究方法對型材壓彎法進行回彈研究對實際生產具有一定的指導意義.

1 工藝及理論分析

圖1所示為機車上使用的一個構件—發(fā)電機風道前托架，由材料為Q235的L截面型材通過壓彎成形.從尺寸上看零件的彎曲半徑較大，很難進行一次性壓彎得到，必須通過多道壓彎才能實現(xiàn)成形.

圖1 零件尺寸

根據(jù)板材的回彈預測公式［8］:

式中，rp為凸模圓角半徑(mm);r0為工件回彈后的半徑(mm);A為材料的彈性模量簡化系數(shù)(A=3σs/E);C為相對彎曲半徑系數(shù);H為型材翼板的高度(mm).

將式(2)作為預測L截面型鋼回彈的公式，確定C值的大小是要解決的主要問題，需要通過對具有不同的相對彎曲半徑r0/H的L截面型材進行試驗得到，試驗時多次試模和修模工作量和費用巨大，因此，需要借助有限元模擬軟件的輔助.

2 模擬分析

2.1 建立有限元模型

隨著計算機技術的快速發(fā)展數(shù)值模擬法被越來越多的應用于解決塑性成形問題，其中DYNAFORM軟件在分析板材的回彈問題方面被廣泛的使用，通過數(shù)值模擬的方法對公式進行檢驗首先，建立有限元模型并劃分網(wǎng)格如圖2所示.單元類型采用BT殼單元，抽取坯料中性層來作為模型.劃分模具和坯料的最小網(wǎng)格尺寸均為為0.5 mm，最大為5 mm.設置凸模運行速度為300 mm/s.利用動力顯式算法計算壓彎過程，再通過靜力隱式計算卸載后的回彈.將每一步完成回彈計算坯料的文件(后綴名為dynain)重新導入到DYNAFORM前處理中，調整坯料與模具的相對位置，再進行下一道壓彎和回彈計算［9］.這樣做的目的是由于前一道工序的模擬結果直接影響下一道工序的模擬精度.通常在模擬板料的多步?jīng)_壓過程時，都忽略了每道工序后板料的回彈，即忽略了板料在每道工序后的應力變化，而板料的回彈計算精度又與其變形歷史密切相關［10］.在多步?jīng)_壓模擬中間加進回彈計算，相對傳統(tǒng)的多步?jīng)_壓模擬過程而然，避免了因忽略每道工序后的回彈而造成的應力變化模擬不正確，進而影響最終回彈精度［11］.成形計算過程采用全積分單元，厚度方向5個積分點，LS－DYNA中提供了12種迭代方法，其中BFGS法是缺省方法，這是目前最成功的一個變尺度算法，是一種擬牛頓法，收斂速度快，且具有較好的數(shù)值穩(wěn)定性.隱式參數(shù)的設置對回彈精度和收斂有重要影響［12-13］.對于初步設計方案，進行每一道壓彎時坯料的送進量為1倍的凸模寬度，發(fā)現(xiàn)在翼板上存在一定的起皺現(xiàn)象，如果減小送料量會改善這種現(xiàn)象，并調整凸模和擋板的間隙為1.1倍的板厚，經(jīng)過多次模擬將每次的送料量確定為2/3倍的凸模寬度，修改后的模型如圖3所示，修改模型后既改善了翼板起皺現(xiàn)象又保證了給料量不至于影響效率，如圖4所示.

圖2 有限元模型的建立及網(wǎng)劃分

圖3 調整坯料給料量后的模型

圖4 兩種方案下的起皺現(xiàn)象分析圖

對截面30 mm×30 mm，材料為Q235的L截面型鋼，在不同的相對彎曲半徑r0/H情況下進行模擬壓彎，并結合試驗測量其回彈后的半徑值r0，再計算回彈半徑Δr=r0－rP，最終得到相對彎曲半徑r0/H與回彈半徑Δr的關系如圖5所示.模擬結果和試驗結果得到較好的吻合，驗證了模擬的可靠性.

圖5 相對彎曲半徑r0/H和回彈半徑Δr的關系

圖5 相對彎曲半徑r0/H與相對彎曲半徑系數(shù)C的關系

2.2 試驗驗證

對圖1所示的發(fā)電機風道前托架進行回彈的預測，相對彎曲半徑r0/H=725/30≈24，彈性模量簡化系數(shù)A=0.005，由圖5可以查得相對彎曲半徑系數(shù)C≈0.7，根據(jù)式(2)設計凸模圓角半徑rp≈619 mm，實際設計凸模圓角半徑為640 mm，與實際吻合的較好.

3 結論

(1)借助有限元模擬軟件DYNAFORM可以對回彈進行預測，并可以預先知道模具設計的不足;

(2)模擬軟件的參數(shù)設置的合理性直接影響模擬結果的準確性，需要結合一定的試驗來對預測精度加以保證;

(3)將板材回彈預測公式中的板厚t換成翼板高度H，并引入相對彎曲半徑系數(shù)C值，修改后的公式可以比較準確的預測L截面型鋼壓彎回彈值，為實際生產中凸模的設計提供參考，減少試模次數(shù)，降低成本.

［1］劉桂花，馮再新，賈青云，等.板料沖壓成形過程中回彈預測及控制的研究進展［J］.金屬鑄鍛焊技術，2011(17):122-125.

［2］余同希，章亮熾.塑性彎曲理論及其應用［M］.北京:科學出版社，1992.

［3］HILL R.The Mathematical Theory of Plasticity［M］.London:［s.n.］，1950.

［4］JOHNSON W，YU T X.Springback after the biaxial elastic-plastic pure bending of a rectangular plate-I［J］.International journal of mechanical sciences，1981，23:619-630.

［5］JOHNSON W，YU T X.On the range of applicability of results for the springback of an elastic－perfectly plastic rectangular plate after subjecting it to biaxial pure bending II［J］.International journal of mechanical sciences，1981，23:631-637.

［6］林忠欽，劉罡，李淑慧，等.應用正交試驗設計提高U型件的成形精度［J］.機械工程學報，2002，38:83-89.

［7］楊雪春，董懿瓊，彭偉，等.U型件回彈控制模具補償法的研究［J］.鍛壓技術，2009，2(1):48-52.

［8］楊玉英，崔令江.實用沖壓工藝及模具設計手冊［M］.北京:機械工業(yè)出版社，2005.

［9］陳文亮.板料成型CAE分析教程［M］.北京:機械工業(yè)出版社，2005.

［10］徐丙坤，施法中.板料沖壓成形回彈的數(shù)值模擬［J］.北京航空航天大學學報，2001(2):194-197.

［11］陳煒，王曉路，高霖.板料多步?jīng)_壓回彈的數(shù)值模擬研究［J］.塑性工程學報，2005(5):12-15.

［12］龔紅英.板料沖壓成形CAE實用教程［M］.北京:化學工業(yè)出版社，2009.

［13］BRADLEY N，MAKER，XINHAI ZHU.Input Parameters for Metal Forming Simulation using LS-DYNA［M］.Livermore:Software Technology Corporation，2000:1-10.

FEA and Prediction Study of Q235 L-Section Steel in Press Bending Springback Process

WANG Da-yong1，SHI Yang1，HUANG Yong-hong2

(1.School of Materical Science Engineering，Dalian Jiaotong University，Dalian 116028，China;2.Dalian Locomotive and Rolling Stock Co.，Ltd，CNR Group，Dalian 116022，China)

Bending process of generator air duct bracket was analysed，and the relative bending radius coefficient C was introduced to revise the traditional plate bending springback prediction formula to obtain springback prediction formula of L-section steel.The bending process of the L-section steel was simulated by DYNAFORM，and the springback was predicted.Value C was

by experiment，and the relationship between value C and relative bending radiusr0/Hwas summed up.The die radius needed can be obtained by the relation curve，then the design of the die and the bending craft was optimized，which is reference to the die design to reduce the experiment and the production cost.

springback;press bending L-section;FEA;relative bending radius coefficient C

A

式中，rp為凸模圓角半徑;r0為板材回彈后的半*

2013-03-28

10.13291/j.cnki.djdxac.2014.03.021

1673-9590(2014)03-0088-04

王大勇(1971－)，男，副教授，學士，主要從事金屬塑性成形技術的研究

E-mail:wangdy@djtu.edu.cn.徑;σS/E為屈服強度和彈性模量的比值;t為板材厚度，而對于同種材料，L截面型鋼相對于等厚的板材只增加了一個高度為H的翼板，翼板對回彈起到限制的作用，板厚的作用不如翼板，所以，假設將式(1)中的t換成H，型材的相對彎曲半徑為r0/H，并且引入相對彎曲半徑系數(shù)C，因此將式(1)修改為