陳高翔， 王亞鋒

(西安西電高壓開關有限責任公司， 陜西 西安 710018)

0 引言

型材彎曲件在汽車、家用電器、船舶、建材等民用領域、導彈及其發(fā)射架制造、飛行器制造等軍用領域的應用日益廣泛.飛機制造中滾彎工藝主要用于成形框肋緣條，機身前后段和發(fā)動機短艙的長桁，這些都是組成飛機骨架的受力零件.滾彎加工過程十分復雜，它涉及到幾何非線性、材料非線性和接觸非線性的多重非線性問題，成型規(guī)律不易把握[1-4].

由于滾彎成形過程影響因素很多，傳統的物理實驗法成本高、周期長，而且存在著各種誤差.近年來，有限元分析被認為是對型材彎曲回彈預測的最有效的方法[5,6].本文將運用顯式動力學有限元分析軟件ANSYS/LS-DYNA對型材滾彎過程進行模擬[5,6],并分析Z型型材滾彎過程腹板中性層的變化，以及成形與回彈過程型材腹板切向應力、應變分布.

1 滾彎成形過程顯式有限元模擬

本文模擬數控四軸型材滾彎機的二維成形過程，型材所用材料為7075-O,滾輪材料為45鋼.在ANSYS/LS-DYNA中建立的有限元模型[7]如圖1所示.型材在上下滾輪的夾持和傳導作用下(靠摩擦力進給)向Z軸正向移動，左右滾輪成角度向上運動，上下滾輪為主轉動輪,且采用分片組裝式結構[8]，左右滾輪在型材的帶動下旋轉.通過改變左右滾輪的上升位置即可獲得不同的滾彎半徑的型材彎曲件.在有限元模擬中，材料模型為經典雙線隨動硬化(BKIN)彈塑性材料模型，型材采用SHELL163單元，滾輪設為剛體.對型材、左右滾輪采用掃略網格，對上下滾輪采用自由網格.各部分材料性能參數見表1.

表1 各部件材料性能參數

圖1 滾彎成形有限元模型

由于在滾彎過程中，型材的進給依靠上下滾輪的摩擦力驅動，故設定上下滾輪繞X軸自轉(轉向不同)，通過設定摩擦系數即可實現型材的進給運動.由于左右滾輪作用在型材上的壓力很難測定，而其移動位移易于測量，因此在文中采取對左右滾輪施加位移載荷來模擬滾彎過程.最終滾彎成形模擬如圖2所示.通過對模擬結果進行測量可知，滾彎后型材變形區(qū)腹板高度減小而其長度卻增加.

圖2 滾彎成形過程模擬圖

2 滾彎成形應力與應變中性層分析

對滾彎區(qū)域沿型材中心位置的腹板高度方向(逐步向內)取9個單元，對比分析滾彎過程中腹板高度方向應力與應變隨時間變化的分布規(guī)律，如圖3所示.從圖3可以知道：

(1) 沿型材腹板高度方向由外側到內側應力、應變變化依次為正值、零值、負值(單元A-I)，則最外層受拉應力，最內層受壓應力.而根據中性層理論，腹板高度的h/2處為中性層，由圖3可知E單元處于中性層位置.

(2) 隨著滾彎過程的進行，應力和應變中性：層不斷波動，并且隨著曲率的增大逐漸向曲率中心方向移動(壓應力方向).

(3) 滾彎開始階段，應力中性層與應變中性層重合，并通過腹板斷面重心，隨著滾彎進行，應力中性層與應變中性層不再重合，而且應力中性層內移量大于應變中性層內移量.

圖3 腹板高度方向應力與應變隨時間變化分布圖

3 回彈前后應力、應變分析

3.1 回彈前后腹板長度方向應力與應變分布

在滾彎加工過程中，型材一部分發(fā)生彎曲，另一部分卻發(fā)生回彈，即回彈也是一個動態(tài)過程.所以，型材滾彎過程中的回彈分析比較復雜[9].本文對回彈前后型材腹板長度方向等效Mises彈性應變和應力進行對比分析，結果如圖4所示.

圖4 回彈前后腹板長度方向應力與應變分布圖

從圖示中可以看出在滾彎的過程中，型材逐步經歷彎曲-再彎曲-回彈變化過程.結合圖a、圖b可知，型材滾彎一段后，左滾輪與上滾輪之間的應力和應變大于右滾輪與上滾輪之間的應力和應變，這是因為型材在滾彎的過程中還在向z軸正向進給，因而左側應力、應變應大于右側.型材在滾彎開始與上滾輪接觸部位應力最大，隨著逐步進給，由于摩擦的產生使得應力、應變峰值降低.脫離接觸后，由于滾彎過程加深，使得這部分達到新的應力、應變峰值.而在左右滾輪接觸的位置,由于滾輪對型材有一向上的力,使得該處在塑性變形區(qū)域內應力較小.而應力趨近于零的部分為型材直邊部分.

從圖示中還可以看出，回彈前后型材與滾輪動態(tài)接觸過程中存儲的彈性能在成形壓力消失后釋放，使得滾彎區(qū)域彈性應變回復，等效Mises應力大大降低.然而，在與上滾輪接觸位置，回彈前后等效Mises彈性應變變化最小.

3.2 回彈前后腹板高度方向應力分布

回彈過程起始于卸載階段，因此，為了更好了解回彈規(guī)律，必須結合彎曲成形最終應力變化和卸載過程應力變化[10].由于型材采用彈塑性材料，滾彎成形后，型材在左右滾輪的作用下，型材腹板高度方向上的切向應力分布如圖5(a)所示.從圖示中可知，成形后，腹板在滾彎拉伸區(qū)應力為正，在壓縮區(qū)應力為負，中性層向內側移動，且型材在成形過程中經歷了從彈性變形到塑性變形的轉變.卸載時，根據卸載定理，相當于施加反向彈性彎矩，故其應力呈線性分布，如圖5(b)所示.卸載后，將卸載前與卸載時切向應力疊加，可得到回彈后型材腹板高度切向應力變化，其應力如圖5(c)所示.

圖5 彈塑性彎曲卸載過程中型材腹板高度方向切向應力變化

圖6 回彈前后腹板高度方向切向應力對比

本文對有限元模擬結果進行了數值分析，得到了型材腹板高度方向在滾彎成形、回彈兩個階段的切向應力分布狀況，如圖6所示.從圖示中可知，數值模擬結果與上文理論分析結果一致.

4 結論

該文運用顯式動力分析軟件ANSYS/LS-DYNA成功實現了對飛機Z型材彎曲件滾彎及其回彈過程的數值模擬.通過比較型材切向應力模擬分布與理論變化，發(fā)現兩者一致，說明采用顯、隱式相結合的方法計算彎曲回彈問題是比較成功的.

模擬結果表明：型材變形區(qū)腹板高度減少及長度增加，腹板應力中性層和應變中性層內移，回彈后，彈性能釋放，等效Mises彈性應變和等效Mises應力的峰值變小，但在型材的起始滾彎位置仍有一個較小的殘余應力集中.

[1] 楊建國.大尺寸板材滾彎過程應力場特點分析[J]．材料科學與工藝，2008，16(5)：684-687．

[2] 孫衛(wèi)東.數控四軸型材滾彎機應用研究[D].西安：西北工業(yè)大學，2001.

[3] 周養(yǎng)萍，亓江文.型材變曲率滾彎過程有限元模擬[J].機械科學與技術，2008，27(8)：1113-1120．

[4] 鄭軍興.冷彎成形有限元數值模擬及回彈分析[D]．昆明：昆明理工大學，2006．

[5] MeiZhan,HeYang,LiangHuang,et al. Springback analysis of numerical control bending of thin-walled tube using numerical-analytic method[J]. Journal of Materials Processing Technology,2006,177:197-201．

[6] S.K.Panthi,N.Ramakrishnan ,K.K.Pathak,et al. An analysis of springback in sheet metal bending using fnite element method(FEM)[J]. Journal of Materials Processing Technology,2007,186:120-124.

[7] 何 濤． ANSYS10.0/LS-DYNA非線性有限元分析實例指導教程[M]．北京：機械工業(yè)出版社，2006．

[8] 周養(yǎng)萍．型材滾彎成形研究[D]．西安：西安理工大學，2007．

[9] 胡軍峰,楊建國，方洪淵，等.滾彎過程的三維動態(tài)仿真模擬[J]．塑性工程學報，2005，12(3)：51-55．

[10] 李文平．彎曲回彈變分原理及其數值模擬研究[D]．秦皇島：燕山大學，2006．