宮林東 韓志仁

（沈陽航空航天大學(xué)航空航天學(xué)部，遼寧 沈陽 110136）

引言

橡皮成形的有限元數(shù)值模擬技術(shù),隨著計(jì)算機(jī)技術(shù)、有限元技術(shù)和塑性加工理論的發(fā)展日趨成熟,將數(shù)值模擬的分析方法應(yīng)用于橡皮成形,是一種新的研究途徑。目前沖壓過程的計(jì)算機(jī)分析仿真（CAE）已經(jīng)能夠在工程實(shí)際中幫助解決傳統(tǒng)方法難以解決的模具設(shè)計(jì)和沖壓工藝的制定等難題，如拉裂的預(yù)測(cè)，回彈的計(jì)算和起皺預(yù)測(cè)等等。

在成形過程中板料與橡皮均會(huì)發(fā)生很大的變形，橡皮直接與板料相作用，因此對(duì)橡皮變形的精確模擬對(duì)于橡皮成形的模擬是非常重要的。

1 橡皮成形中的回彈問題

回彈是板料成形卸載后,彈性變形部分產(chǎn)生回復(fù)的一種現(xiàn)象,目前解決回彈問題的有效途徑,是通過有限元軟件,采用數(shù)值模擬的方法進(jìn)行回彈預(yù)測(cè)、回彈補(bǔ)償,并進(jìn)行合理修模的方法。因此,回彈的數(shù)值模擬預(yù)測(cè)精度的準(zhǔn)確性,將不僅影響成形零件的質(zhì)量,并對(duì)改善航空企業(yè)傳統(tǒng)的不修?；蛲ㄟ^人工校形的生產(chǎn)方式產(chǎn)生影響。

本文利用ANSYS軟件,對(duì)橡皮成形中飛機(jī)機(jī)翼前緣蒙皮零件進(jìn)行有限元數(shù)值模擬,建立有限元數(shù)值模型,對(duì)單元積分類型、網(wǎng)格尺寸方案,以及接觸與摩擦等影響因素進(jìn)行量化研究,旨在提高橡皮成形數(shù)值模擬的回彈精度預(yù)測(cè),通過模擬結(jié)果,獲得在成形與回彈預(yù)測(cè)的精度上更精準(zhǔn)的模型及改進(jìn)方案。

2 橡皮成形的有限元模擬

2.1 有限元模型

顯式算法基于動(dòng)力學(xué)方程，因此無需迭代；而靜態(tài)隱式算法基于虛功原理，一般需要迭代計(jì)算。此成形模擬過程采用顯式算法,回彈模擬過程采用隱式算法。

模擬分析的材料參數(shù)如下：鈑料為AL2B06-T3，厚度1.5mm，密度為2.87e-6g/mm3，屈服強(qiáng)度為 400MPa,切線模量為100Mpa,泊松比為0.3,楊氏模量為73000MPa。剛性凸模的密度為7.58e-6g/mm3，泊松比為0.3,楊氏模量為2.07 e5 MPa。橡皮的密度為1.15e-6g/mm3，剪切模量為 1040MPa。

鈑料的橡皮成形采用ANSYS/LS-DYNA進(jìn)行數(shù)值模擬。

2.2 網(wǎng)格尺寸對(duì)成形的影響

合理的網(wǎng)格劃分是保證獲得準(zhǔn)確模擬結(jié)果的前提，單元尺寸的大小會(huì)直接影響算精度。網(wǎng)格尺寸可以小一些以保精確的模具表面描述，能充分表示板料的變形特征。Mattiasson發(fā)現(xiàn)回彈值與單元尺寸成反比關(guān)系,單元尺寸越小,計(jì)算的回彈量越大；反之，回彈量越小,甚至可能出現(xiàn)反向現(xiàn)象。

通過模擬分析得到：回彈角與網(wǎng)格尺寸成反比,即網(wǎng)格尺寸越小,計(jì)算的回彈角越大，反之，回彈量越小，同時(shí)對(duì)橡皮成形應(yīng)力應(yīng)變的影響極大，導(dǎo)致模擬結(jié)果不準(zhǔn)確。因此，在保證回彈預(yù)測(cè)精度的前提下，對(duì)模具幾何形狀變化劇烈的地方應(yīng)劃分尺寸較小的單元，其他地方可采用較大的單元尺寸。

2.3 摩擦系數(shù)大小對(duì)成形的影響

接觸問題是塑性成形模擬領(lǐng)域的難點(diǎn)問題之一。在板料成形過程中板料上的作用力是通過板料表面與模具表面的接觸來傳遞的。模具表面為主約束表面，板料表面為從表面。如何正確判定主從表面的接觸狀態(tài)并選擇合適的位移約束條件和邊界摩擦條件，將直接影響到模擬運(yùn)算的運(yùn)行。摩擦對(duì)零件回彈的影響較大，隨著摩擦系數(shù)的增大，板材回彈減小。但是過大的摩擦系數(shù)又會(huì)讓鈑料表面質(zhì)量變壞，網(wǎng)格會(huì)撕裂。

2.4 ANSYS回彈求解

板料彎曲及回彈的數(shù)值模擬過程是:凸模向下移動(dòng)壓迫板料彎曲變形,在彎曲成形過程結(jié)束后,把顯式單元類型轉(zhuǎn)變?yōu)殡[式,刪除材料模型中所有的非線性材料性質(zhì),將彎曲后的幾何形狀作為回彈求解的初始構(gòu)形,從顯式分析的結(jié)果中輸入單元應(yīng)力數(shù)據(jù),去除模具,對(duì)鈑料施加相關(guān)載荷,然后即可求解。圖1為成形零件的mises屈服應(yīng)力分析云圖,圖2為鈑料的回彈示圖。

圖1 mises屈服應(yīng)力分析云圖

然后，在y=275面截鈑料所得的曲線上選取有限個(gè)節(jié)點(diǎn)，根據(jù)節(jié)點(diǎn)的坐標(biāo)前后的變化得到其回彈數(shù)據(jù)圖表。圖3回彈前、后的節(jié)點(diǎn)坐標(biāo)和回彈示圖

圖2 回彈示圖

圖3 回彈前、后的節(jié)點(diǎn)坐標(biāo)和回彈示圖

利用航空手冊(cè)上的經(jīng)驗(yàn)數(shù)據(jù)對(duì)比模擬回彈數(shù)值：

1、模擬的回彈角?；貜椙暗慕嵌葹?25°，回彈后的角度為131°回彈角為6°。

2、手冊(cè)中的經(jīng)驗(yàn)回彈角?；貜椙暗慕嵌葹?125°，回彈后的角度為 131°30′，回彈角為7°10′，回彈角比經(jīng)驗(yàn)值小 1°10′。

3 結(jié)論

(1)網(wǎng)格尺寸與回彈角成反比,即網(wǎng)格尺寸越小,計(jì)算的回彈角越大，要使回彈角小，網(wǎng)格尺寸應(yīng)該大，但同時(shí)網(wǎng)格尺寸越大，對(duì)橡皮成形應(yīng)力應(yīng)變的影響極大，導(dǎo)致模擬結(jié)果不準(zhǔn)確。通過模擬我們選擇網(wǎng)格大小為長2mm、寬0.25mm。(2)摩擦對(duì)零件回彈的影響較大，隨著摩擦系數(shù)的增大，板材回彈減小。但是摩擦系數(shù)過大又會(huì)讓鈑料表面質(zhì)量變壞，還有網(wǎng)格會(huì)撕裂。經(jīng)過多次模擬和經(jīng)驗(yàn)數(shù)據(jù)對(duì)比得出，選用動(dòng)摩擦系數(shù):0.15，靜摩擦系數(shù):0.1。(3)回彈角度比手冊(cè)經(jīng)驗(yàn)值略小，有限元模擬值與試驗(yàn)經(jīng)驗(yàn)值的趨勢(shì)基本一致,回彈模擬值比試驗(yàn)值偏小?？梢宰鳛樵O(shè)計(jì)模具時(shí)的參考依據(jù)。

[1]陳磊,王俊彪等.橡皮液壓成形直彎邊回彈試驗(yàn)與分析研究[J].塑性工程學(xué)報(bào),2008，15(3):47-48

[2]Leu Daw Kwei.A simplified approach for evaluating bendability and springback in plastic bending of aniso-tropic sheet metals[J].Journal of Materials Processing Technology,1997，（66）:9-17

[3]Flavio Cimolin,Roberto Vadori,Claudio Canuto.Spri-ngback compensation in deep drawing applications[J].International Journal of Material Forming,2008， （43）:101-113

[4]趙國偉,王元?jiǎng)?陳建橋,李榮峰,祝洪川.板料成形回彈的數(shù)值模擬與影響因素[J].鍛壓裝備與制造技術(shù),2005（3）:55-56

[5]李瀧杲,王書恒等.金屬板料成形有限元模擬基礎(chǔ)[M].北京:北京航空航天大學(xué)出版社,2008.55-60