王陶+雷君相+彭彥君+高文靜+郭豐偉

摘要： 利用ABAQUS 6.14軟件模擬了在不同壁厚下，材料為20鋼雙直徑圓管吸能元件的液壓脹形、折疊和自由翻轉(zhuǎn)的過程.并分析了不同壁厚下液壓脹形與折疊后的壁厚分布，準(zhǔn)靜態(tài)壓縮翻轉(zhuǎn)變形模式以及準(zhǔn)靜態(tài)壓縮時翻轉(zhuǎn)力與翻轉(zhuǎn)半徑與推導(dǎo)公式之間的關(guān)系.結(jié)果表明：壁厚減薄程度與管材初始壁厚大小有關(guān)；20鋼雙直徑圓管成形后發(fā)生軸向壓縮時其自由翻轉(zhuǎn)的模式多為內(nèi)管外翻；推導(dǎo)的翻轉(zhuǎn)力與翻轉(zhuǎn)半徑理論公式與實際測得翻轉(zhuǎn)力與翻轉(zhuǎn)半徑的結(jié)果比較吻合.

關(guān)鍵詞： 液壓脹形與折疊； 雙直徑圓管； 翻轉(zhuǎn)力； 翻轉(zhuǎn)半徑

中圖分類號： TG 386文獻(xiàn)標(biāo)志碼： A

Influence of Tube Thickness on the Forming and Free

Inversion of Energy Absorption of Double Diameter Tube

WANG Tao，LEI Junxiang，PENG Yanjun，GAO Wenjing，GUO Fengwei

（School of Materials Science & Engineering， University of Shanghai for

Science and Technology， Shanghai 200093， China）

Abstract： Using ABAQUS 6.14 software to simulate the process of hydraulic forming and folding and free inversion of the original 20 steel double diameter tube with different thickness.The relationship between the thickness distribution of the hydraulic forming and folding，the mode of quasistatic compression deformation and the relationship between the inversion force and the inversion radius of the quasistatic compression and the formula are analyzed.The results show that the tube thickness reduction is related to the degree of initial tube thickness.The mode of selfreversal in the axial compression of 20pipe doublediameter pipe is mostly inner tube valgus.The theoretical formula of the reversal force and the turning radius is consistent with the actual measured turning force and the turning radius.

Keywords： hydraulic forming and folding； double diameter tube； inversion force； inversion radius

近年來，隨著公共安全的保護(hù)得到越來越多的關(guān)注，人們加強了關(guān)于結(jié)構(gòu)耐撞性和能量吸收裝置的研究[1-2].雙直徑圓管作為一種有效吸能元件，能夠?qū)崿F(xiàn)載荷最為穩(wěn)定的目的[3]，因而國內(nèi)外對其研究也越來越多[4].

利用內(nèi)高壓成形技術(shù)[5]，制備出合格的雙直徑圓管吸能元件.雙直徑圓管自由翻卷能量吸收結(jié)構(gòu)載荷恒定平穩(wěn)，成為材料和結(jié)構(gòu)的能量吸收中最理想的能量吸收結(jié)構(gòu)[6-7].在對雙直徑圓管施加準(zhǔn)靜態(tài)載荷時，在圓角變形區(qū)自由翻轉(zhuǎn)時的3種主要有效變形模式：雙直徑圓管內(nèi)管外翻、雙直徑圓管外管內(nèi)翻和雙直徑圓管內(nèi)管外翻與外管內(nèi)翻復(fù)合[8-9].翻轉(zhuǎn)力隨著軸向位移的不斷增加而增大，逐漸趨向于穩(wěn)定[10].

本文通過對雙直徑圓管吸能元件的結(jié)構(gòu)設(shè)計，模擬20鋼管材在不同壁厚下液壓脹形與折疊過程，分析影響管材壁厚減薄與初始壁厚的關(guān)系.不同壁厚下，模擬20鋼雙直徑圓管的準(zhǔn)靜態(tài)壓縮過程.對比分析各自載荷位移曲線及其自由翻轉(zhuǎn)模式.最后對準(zhǔn)靜態(tài)壓縮時翻轉(zhuǎn)力、翻轉(zhuǎn)半徑與推導(dǎo)公式之間的關(guān)系進(jìn)行分析.

1有限元模擬

第3期王陶，等：管材壁厚對雙直徑圓管吸能元件成形與自由翻轉(zhuǎn)的影響有 色 金 屬 材 料 與 工 程2017年 第38卷1.1成形過程建模

在ABAQUS 6.14軟件上建立三維有限元，模擬模型見圖1，沖頭和模具型腔幾何形狀為三維殼體，類型為剛體，管坯為三維變形體.管坯劃分采用六面體（Hex）單元，結(jié)構(gòu)化網(wǎng)格劃分技術(shù).為了觀察圓角處的變形，將管坯劃分為3層.沖頭和模具型腔為剛體，不需要劃分網(wǎng)格.管坯和模具型腔之間摩擦接觸，摩擦因數(shù)為0.1；管坯與沖頭之間為無摩擦接觸.利用ABAQUS 6.14軟件分別對壁厚t為1.2，1.4，1.6和1.8 mm的20無縫鋼管的雙直徑圓管進(jìn)行模擬液壓脹形與折疊過程.成形后的雙直徑圓管零件需要滿足壁厚減薄率小、變形區(qū)足長以及成形區(qū)壁厚分布均勻的條件，方能稱作合格的雙直徑圓管零件，見圖2.雙直徑圓管吸能元件液壓脹形與折疊在一副模具中完成且能一次成形出兩個管件.材料參數(shù)見表1.

1.2液壓脹形與折疊和自由翻轉(zhuǎn)過程模擬

根據(jù)圖3中液壓加載路徑以及成形區(qū)間圖對壁厚t為1.2，1.4，1.6和1.8 mm的20無縫鋼管的雙直徑圓管進(jìn)行模擬液壓脹形與折疊過程.對成形后合格的雙直徑圓管，分別模擬不同壁厚下的準(zhǔn)靜態(tài)壓縮過程.

2模擬結(jié)果分析

2.1成形零件壁厚分布比較

管材在液壓脹形與折疊過程中壁厚的減薄，對雙直徑圓管的自由翻轉(zhuǎn)過程影響較大.最大減薄率是判斷該零件成形是否合格的重要標(biāo)準(zhǔn)[11].最大減薄率：

ε=δ0-δδ0（1）

式中：δ0為管材初始壁厚；δ為成形后最小壁厚.

利用ABAQUS 6.14軟件的后處理模塊，模擬不同壁厚下管材液壓脹形與折疊成形后的零件壁厚分布，如圖4所示.從圖4中可以看出，最大壁厚減薄處均發(fā)生在外翻圓角“S”型過渡區(qū)處.4種不同壁厚下最大增厚處在管材兩端，出現(xiàn)的原因是沖頭在軸向進(jìn)給的時候，兩端受到擠壓造成鐓粗所致.初始壁厚為1.8 mm的雙直徑圓管最大減薄率為13.48%，初始壁厚為1.2 mm的雙直徑圓管最大減薄率為19.88%；初始壁厚為1.4 mm與1.6 mm的雙直徑圓管最大減薄率初始壁厚為1.2～1.8 mm.從整體來看，4種不同初始壁厚管材的壁厚減薄的變化相對比較均勻，滿足成形質(zhì)量要求.

影響壁厚減薄的因素有很多，除了與管材自身的材料力學(xué)性能有關(guān)，還與管材與模具之間的摩擦、模具圓角半徑以及液壓脹形與折疊過程中材料的補充速度有關(guān).本次試驗中，模具的圓角半徑對壁厚的模擬過程的影響較大.模具圓角半徑的尺寸直接影響管材在液壓脹形與折疊過程中的材料的流動性，進(jìn)而影響壁厚的減薄率以及壁厚的均勻度，同時合適的圓角半徑能夠減小加工硬化帶來的不利影響.此次模擬過程中，模具圓角半徑為3 mm，通過合理控制液壓力與軸向進(jìn)給量來控制材料的補充速度，可以控制雙直徑圓管在成形過程中避免變形區(qū)出現(xiàn)起皺、屈曲和破裂的失效現(xiàn)象，進(jìn)而控制壁厚減薄率，管壁厚更為均勻，提高成形合格率.

2.2自由翻轉(zhuǎn)模式與翻轉(zhuǎn)載荷模擬結(jié)果分析

液壓脹形與折疊后，對4種不同壁厚的雙直徑圓管進(jìn)行準(zhǔn)靜態(tài)壓縮變形模擬，壓縮變形的翻轉(zhuǎn)載荷位移曲線如圖5所示.

從圖5中可以看出，4種不同初始壁厚的雙直徑圓管在發(fā)生壓縮變形時，4種雙直徑圓管均能達(dá)到載荷穩(wěn)定階段.初始壁厚為1.8 mm的雙直徑圓管自由翻轉(zhuǎn)載荷上升幅度較快，載荷相對平穩(wěn)，承受的翻轉(zhuǎn)載荷相對較大.初始壁厚為1.2 mm的雙直徑圓管自由翻轉(zhuǎn)載荷在趨于平穩(wěn)過程中存在一定的波動，達(dá)到載荷穩(wěn)定時的行程相對較長.初始壁厚為1.4 mm與1.6 mm的雙直徑圓管自由翻轉(zhuǎn)載荷狀態(tài)介于初始壁厚1.2～1.8 mm.

影響翻轉(zhuǎn)載荷的穩(wěn)定因素有很多，除了材料的性能，還與成形后雙直徑圓管內(nèi)外管的壁厚有關(guān).由于內(nèi)外管壁厚直接影響雙直徑圓管軸向壓縮變形模式，進(jìn)而影響翻轉(zhuǎn)載荷的平穩(wěn).初始壁厚為1.8 mm的雙直徑圓管外管減薄相對較少，受材料加工硬化的影響較大，外管不容易變形，屈服應(yīng)力增大，內(nèi)管向外翻轉(zhuǎn)趨勢比較大，故翻轉(zhuǎn)模式以內(nèi)管外翻為主，載荷較為穩(wěn)定.初始壁厚為1.2 mm的雙直徑圓管外管壁厚相對較大，外管有向內(nèi)翻轉(zhuǎn)的趨勢，但加工硬化的影響仍存在，故在到達(dá)載荷穩(wěn)定階段會出現(xiàn)一段波動，翻轉(zhuǎn)模式會出現(xiàn)小段內(nèi)管外翻的階段，但仍然以內(nèi)管外翻為主要翻轉(zhuǎn)模式.初始壁厚為1.4 mm與1.6 mm的雙直徑圓管自由翻轉(zhuǎn)的模式也以內(nèi)管外翻為主，承受載荷穩(wěn)定.

由文獻(xiàn)[12]推導(dǎo)的雙直徑圓管自由翻轉(zhuǎn)吸能元件發(fā)生內(nèi)管外翻時的翻轉(zhuǎn)力公式以及翻轉(zhuǎn)半徑公式[12]為：

P=πdtβS—ln1+2rdd+tt+2rdKd（2）

r=d2KdKdtd+34td2-t2d（3）

式中：β與材料的本構(gòu)關(guān)系有關(guān)，各向異性時，β=2（1+rd）1+2rd；S—為材料的冪指強化系數(shù)；d為內(nèi)管外徑；Kd為內(nèi)管外翻系數(shù)，Kd=ln1+4rdd.

結(jié)合雙直徑圓管的尺寸，得出4種壁厚雙直徑圓管準(zhǔn)靜態(tài)壓縮下的模擬結(jié)果和理論結(jié)果，如表2所示.

初始壁厚/mm翻轉(zhuǎn)形式穩(wěn)定翻轉(zhuǎn)力/kN理論翻轉(zhuǎn)力/kN外翻圓角半徑/mm理論圓角半徑/mm1.2內(nèi)管外翻22.3421.151.281.341.4內(nèi)管外翻25.1424.901.361.451.6內(nèi)管外翻31.4729.221.491.581.8內(nèi)管外翻37.1236.871.621.66

同時將上述結(jié)果與Guist & Marble提出的翻轉(zhuǎn)力與翻轉(zhuǎn)半徑公式[13-14]和Reddy提出的翻轉(zhuǎn)力與翻轉(zhuǎn)半徑公式[15-16]所得出的結(jié)果進(jìn)行對比，如圖6所示.

從圖6中可以看出，文獻(xiàn)[12]推導(dǎo)的雙直徑圓管自由翻轉(zhuǎn)吸能元件發(fā)生內(nèi)管外翻時的翻轉(zhuǎn)力公式及翻轉(zhuǎn)半徑公式得出的結(jié)果，較Guist & Marble和Reddy提出的翻轉(zhuǎn)力與翻轉(zhuǎn)半徑公式所得出的結(jié)果更接近模擬值.說明該公式使用在雙直徑圓管發(fā)生內(nèi)管外翻時預(yù)測翻轉(zhuǎn)力與翻轉(zhuǎn)半徑較為準(zhǔn)確.

3試驗分析

對成形后的雙直徑圓管在Zwick試驗機上進(jìn)行軸向壓縮試驗.成形后的雙直徑圓管零件和準(zhǔn)靜態(tài)壓縮后剖視的結(jié)果如圖7所示.

準(zhǔn)靜態(tài)壓縮后載荷位移曲線如圖8所示.

準(zhǔn)靜態(tài)壓縮試驗中測得的翻轉(zhuǎn)力和翻轉(zhuǎn)半徑的大小與理論計算的結(jié)果對比如圖9所示.

從圖7中可以看出，在準(zhǔn)靜態(tài)壓縮試驗過程中，雙直徑圓管在發(fā)生自由翻轉(zhuǎn)時，4種不同壁厚的雙直徑圓管翻轉(zhuǎn)模式并非剛開始就發(fā)生內(nèi)管外翻，而是開始翻轉(zhuǎn)階段發(fā)生了小部分的外管內(nèi)翻，隨后發(fā)生內(nèi)管外翻，總體趨勢以內(nèi)管外翻為主.出現(xiàn)此類情況的原因，可能與管材液壓脹形成形后內(nèi)外壁的壁厚以及成形過程中加工硬化有關(guān)，導(dǎo)致不能以單一的翻轉(zhuǎn)模式進(jìn)行.

從圖8中的準(zhǔn)靜態(tài)壓縮試驗結(jié)果和模擬結(jié)果可以看出，載荷位移曲線的變化趨勢大致相同，試驗中承受的載荷值比模擬結(jié)果略大.但整體來看，模擬結(jié)果與試驗結(jié)果基本吻合.

從圖9中可明顯看出，實際測得的翻轉(zhuǎn)力比理論值大，實際測得的翻轉(zhuǎn)半徑比理論值小.但是較Guist & Marble提出的翻轉(zhuǎn)力和翻轉(zhuǎn)半徑公式與Reddy提出的翻轉(zhuǎn)力與翻轉(zhuǎn)半徑公式得出的結(jié)果更接近實際值.

4結(jié)論

（1） 4種不同壁厚的管材均在合適的加載路徑下成形出合格的雙直徑圓管吸能元件，不同壁厚的管材的減薄率各不相同，且與壁厚有關(guān).初始壁厚為1.8 mm的雙直徑圓管的最大減薄率為13.14%.

（2） 相同條件下，4種不同壁厚的20鋼雙直徑圓管在發(fā)生自由翻轉(zhuǎn)時均以內(nèi)管外翻的翻轉(zhuǎn)模式為主，翻轉(zhuǎn)載荷穩(wěn)定可控，可以用作吸能元件和各種碰撞吸能場合.

（3） 從模擬與試驗相結(jié)合的角度來看，當(dāng)雙直徑圓管發(fā)生內(nèi)管外翻時，由文獻(xiàn)[12]推導(dǎo)的雙直徑圓管吸能元件翻轉(zhuǎn)力公式以及翻轉(zhuǎn)半徑公式得出的結(jié)果更為精確可靠.

參考文獻(xiàn)：

[1]殷之平，李玉龍，黃其青.薄壁圓筒結(jié)構(gòu)耐撞性優(yōu)化設(shè)計[J].機械科學(xué)與技術(shù)，2010，29（2）：220-224.

[2]ZHANG X，HU H H.Energy absorption of longitudinally grooved square tubes under axial compression[J].ThinWalled Structures，2009，47（12）：1469-1477.

[3]余同希，盧國興.材料與結(jié)構(gòu)的能量吸收[M].華云龍，譯.北京：化學(xué)工業(yè)出版社，2006.

[4]張立玲，高峰，杜發(fā)榮.不等厚拼焊管軸向壓縮變形的數(shù)值模擬研究[J].塑性工程學(xué)報，2006，13（4）：6-9，19.

[5]ALHASSANI S T S，JOHNSON W，LOWE W T.Characteristics of inversion tubes under axial loading[J].Journal of Mechanical Engineering Science，1972，14（6）：370-381.

[6]苑世劍，王小松.內(nèi)高壓成形技術(shù)研究與應(yīng)用新進(jìn)展[J].塑性工程學(xué)報，2008，15（2）：22-30.

[7]汪奇超，雷君相，駱協(xié)海.鋁合金防碰撞吸能管液壓成形加載路徑研究[J].鍛壓技術(shù)，2011，36（6）：55-58.

[8]邱信明，賀良鴻.圓管準(zhǔn)靜態(tài)翻轉(zhuǎn)的三維模型分析[J].失效分析與預(yù)防，2011，6（1）：1-7.

[9]劉鋼，董哲，苑世劍.加載路徑對TAl5板材熱拉深氣脹成形筒形件減薄率的影響[C]∥第十二屆全國塑性工程學(xué)術(shù)年會暨第四屆全球華人塑性加工技術(shù)研討會.重慶：中國機械工程學(xué)會，2011.

[10]鄧天根，雷君相.變徑管自由翻轉(zhuǎn)變形的趨向性研究[J].鍛壓技術(shù)，2014，39（2）：144-149，153.

[11]岳永保，楊合，詹梅，等.薄壁管小彎曲半徑數(shù)控彎曲壁厚減薄實驗研究[J].鍛壓技術(shù)，2007，32（5）：58-62.

[12]徐龍江，雷君相，高貴杰.變徑管與薄壁圓管軸向壓縮過程研究[J].機械工程與自動化，2015（3）：137-138，144.

[13]GUIST L R，MARBLE D P.Prediction of the inversion load of a circular tube[R].Washington，D.C.：NASA，1966.

[14]TOMESANI L.Analysis of a tensiondriven outsidein tube inversion[J].Journal of Materials Processing Technology，1997，64（1/3）：379-386.

[15]REDDY T Y.Tube inversionan experiment in plasticity[J].International Journal Mechanical Engineering Education，1989，17（4）：277-291.

[16]REDDY T Y.Guist and Marble revisitedon the natural knuckle radius in tube inversion[J].International Journal of Mechanical Sciences，1992，34（10）：761-768.有色金屬材料與工程第38卷第3期NONFERROUS METAL MATERIALS AND ENGINEERINGVol.38No.32017