宋飛飛，吳金平，楊英麗，趙 彬，羅媛媛，郭荻子

(西北有色金屬研究院，陜西 西安 710016)

Ti35合金管材數(shù)控彎曲成形規(guī)律有限元模擬研究

宋飛飛，吳金平，楊英麗，趙 彬，羅媛媛，郭荻子

(西北有色金屬研究院，陜西 西安 710016)

利用Abaqus有限元軟件對(duì)新型耐蝕Ti35合金管材的數(shù)控彎曲過(guò)程進(jìn)行了模擬研究。研究了彎曲角度、芯棒伸出量、壓塊相對(duì)助推速度和相對(duì)彎曲半徑對(duì)Ti35合金管材成形結(jié)果的影響。結(jié)果表明，Ti35合金管材數(shù)控彎曲截面扁化率和回彈角隨彎曲角度的增大而增大;彎曲變形越劇烈(如減小彎曲半徑、壓塊相對(duì)助推速度，或增大芯棒伸出量)，壁厚減薄率越大，回彈角越小。截面扁化率隨芯棒伸出量、相對(duì)彎曲半徑的增大而減小。

Ti35合金;鈦合金管;數(shù)控彎曲;有限元模擬

0 引言

在石油、煤、天然氣等石化能源日趨緊張的今天，核能以其能量密集、功率高、清潔等優(yōu)點(diǎn)正逐漸成為傳統(tǒng)石化能源的替代品。核乏燃料后處理在核工業(yè)燃料循環(huán)體系中占有重要地位，其主要目的是回收寶貴的裂變?nèi)剂希?］，由于核乏燃料后處理關(guān)鍵設(shè)備長(zhǎng)期在強(qiáng)酸性與放射性環(huán)境中工作，其關(guān)鍵設(shè)備用材料的研發(fā)受到許多學(xué)者的關(guān)注［2－3］。鈦及鈦合金很容易與氧結(jié)合生成氧化膜，在硝酸、濃硫酸等存在氧化劑的溶液中可以保持鈍化，具有優(yōu)良的耐蝕性。Ti35合金是一種Ti－Ta系二元α型鈦合金，具有良好的加工性能及耐腐蝕性能，特別在具有氧化性陽(yáng)離子的核乏燃料模擬料液環(huán)境中，其腐蝕速率＜0.01 mm/a，耐腐蝕性能優(yōu)異，作為超低碳不銹鋼的替代材料，具有在核乏燃料后處理關(guān)鍵設(shè)備上應(yīng)用的巨大前景［4］。

作為一種新型的設(shè)備用材，Ti35合金管材主要應(yīng)用于換熱管和料液管路或管道等部位，其中管路或管道必然存在彎曲的部分。數(shù)控彎管技術(shù)以其能夠?qū)崿F(xiàn)管材批量、快速生產(chǎn)的特點(diǎn)，被廣泛應(yīng)用于管材的彎曲成形制造當(dāng)中［5］。然而，管材彎曲后會(huì)發(fā)生管壁減薄、截面扁化及回彈等現(xiàn)象，這些現(xiàn)象如果過(guò)度顯著，超過(guò)使用要求，就會(huì)成為成形缺陷，如彎曲外側(cè)斷裂、截面扁化嚴(yán)重、彎曲內(nèi)側(cè)起皺嚴(yán)重及過(guò)度回彈等。相比于不銹鋼，鈦合金屈強(qiáng)比高，造成塑性加工區(qū)間小，屈彈比高，造成回彈更加顯著［6］，其管材彎曲性能及成形質(zhì)量嚴(yán)重影響管路或管道的使用性能。因此，有必要對(duì)Ti35合金管材的彎曲性能進(jìn)行系統(tǒng)研究，從而為工業(yè)化應(yīng)用Ti35合金設(shè)備的設(shè)計(jì)提供相關(guān)依據(jù)。

1 研究方法

1.1 三維彈塑性有限元建模

數(shù)控彎管成形是多模具耦合作用的復(fù)雜非線性過(guò)程(如圖1所示)，理論解析法因存在假設(shè)而導(dǎo)致研究精度有限，實(shí)驗(yàn)法又耗費(fèi)大量人力物力，而有限元法因可以處理復(fù)雜非線性過(guò)程，成為處理該問(wèn)題的最佳選擇?；贏baqus軟件，兼顧計(jì)算精度和計(jì)算效率，對(duì)管材彎曲和回彈過(guò)程分別采用動(dòng)態(tài)顯式(Dynamic，Explicit)和靜態(tài)隱式(Static，General)算法求解。管材選擇六面體8節(jié)點(diǎn)減縮積分實(shí)體單元(C3D8R)，該單元具有沙漏控制和減縮積分功能，厚度方向選擇4個(gè)積分點(diǎn)來(lái)滿足精度要求。模具選擇離散剛體單元(R3D4)。對(duì)不同模具與管材間的摩擦采用庫(kù)倫摩擦模型描述接觸面間的摩擦行為，對(duì)管材內(nèi)表面與芯棒、外表面與模具的接觸分別采用罰函數(shù)法和動(dòng)態(tài)約束算法?？紤]到管材彎曲對(duì)稱性，采用半管模型進(jìn)行計(jì)算。在網(wǎng)格劃分問(wèn)題上，為了提高計(jì)算精度，對(duì)管材的彎曲變形區(qū)劃分較細(xì)的網(wǎng)格，而對(duì)各模具以及非變形區(qū)的管材采用較粗的網(wǎng)格進(jìn)行劃分［7］。

圖1 管材數(shù)控彎曲示意圖Fig.1 Principle of tube NC bending

依據(jù) GB/T 228—2010，對(duì)規(guī)格為 φ30 mm×4 mm的Ti35合金管進(jìn)行整管室溫拉伸試驗(yàn)，試驗(yàn)設(shè)備為INSTRON 5985萬(wàn)能材料試驗(yàn)機(jī)。獲得管材拉伸曲線后，采用冪指數(shù)化加工硬化模型σ=k(ε－a)n來(lái)擬合管材屈服后的應(yīng)力應(yīng)變行為，試驗(yàn)及擬合結(jié)果如表1所示，該數(shù)據(jù)將錄入軟件材料模型中來(lái)描述變形行為。

表1 Ti35合金管材料參數(shù)Table 1 Material parameters of Ti35 tube

1.2 試驗(yàn)方案

采用單因素法進(jìn)行研究，選取的工藝參數(shù)均對(duì)彎曲結(jié)果有顯著影響，包括彎曲角度α、芯棒伸出量e、壓塊相對(duì)助推速度Vp/V和相對(duì)彎曲半徑R/D，每因素至少采用3水平，如表2所示，其中粗體為默認(rèn)設(shè)置。考察的指標(biāo)包括Ti35合金管材彎曲過(guò)程中的缺陷形式，是否有起皺、凸包、破裂及其它缺陷出現(xiàn)。如果可以順利彎曲成形，則要考察彎曲后管材的形狀尺寸精度和性能，包括外脊線的壁厚、圓度、回彈等。

表2 Ti35合金管材數(shù)控彎曲工藝參數(shù)模擬方案Table 2 Simulation scheme of forming parameters for Ti35 alloy tube numerical-control bending

1.3 結(jié)果表征方法

Ti35合金管材彎曲后的圓度可以用截面扁化率表征，外脊線的壁厚可以用壁厚減薄率表征。采樣頻率為從彎曲起始面至彎曲終止面，每隔幾度的測(cè)量結(jié)果，其計(jì)算方法分別如式(1)、(2)所示。彎曲回彈可以用回彈角來(lái)表征，即彎曲角度與回彈后的最終成形角度的差值，如式(3)所示。

式中，D0和D1分別為Ti35合金管材彎曲前外徑和彎曲后短軸長(zhǎng)度，t0和t1分別為彎曲前壁厚及彎曲后外脊線處壁厚，α0和α1分別為彎曲角度和成形角度。

2 結(jié)果與討論

2.1 不同彎曲角度的影響

圖2 彎曲角度對(duì)壁厚減薄率、截面扁化率和回彈角的影響Fig.2 Effect of bending angle on wall thinning ratio，ovalization ratio and springback angle

管材在彎曲變形時(shí)，外側(cè)脊線承受最大的拉應(yīng)力，因此，選擇對(duì)管材外脊線處的壁厚進(jìn)行測(cè)量。圖2a為不同彎曲角度下管材外脊線處的壁厚減薄率變化規(guī)律?？梢钥闯觯煌瑥澢嵌认?，壁厚減薄率的變化趨勢(shì)是一致的，最大壁厚減薄發(fā)生部位均靠近彎曲后端。因?yàn)閺澢蠖说玫叫景舻闹危谛景舻淖饔孟?，此處管材?yīng)力狀態(tài)更為復(fù)雜，減薄率較大。圖2b是彎曲角對(duì)截面扁化率的影響圖，可以看出，截面最大扁化率隨彎曲角的增加而增大，這是因?yàn)閺澢窃酱?，變形區(qū)域越大，芯棒對(duì)管材的支撐作用就相對(duì)減弱，引起扁化率增大。通過(guò)隱式算法求解回彈過(guò)程，圖2c所示為不同彎曲角度下的回彈角變化規(guī)律，可以看出回彈角隨彎曲角的增大而增大，且和彎曲角基本呈線性關(guān)系。實(shí)踐中常通過(guò)該規(guī)律，對(duì)實(shí)際控制成形角進(jìn)行指導(dǎo)［8］，即首先進(jìn)行20°、120°彎曲，擬合出回彈角－彎曲角關(guān)系曲線，再根據(jù)所要成形的角度，反推出彎曲角度。

2.2 芯棒伸出量的影響

圖3a為芯棒伸出量對(duì)壁厚減薄率的影響圖。從圖中可以看出壁厚減薄率在整個(gè)彎曲段變化不大，隨著芯棒伸出量的增大，最大減薄率逐漸增大。這是因?yàn)殡S著芯棒伸出量的增加，阻礙了材料順利越過(guò)彎曲切點(diǎn)向前流動(dòng)，這樣在后續(xù)的彎曲過(guò)程中，由于沒(méi)有后續(xù)材料的補(bǔ)償，外側(cè)因受拉應(yīng)力而產(chǎn)生的壁厚減薄就會(huì)加劇。圖3b為芯棒伸出量對(duì)截面扁化率的影響圖。從圖中可以看出，不同芯棒伸出量下的管材截面扁化率的變化規(guī)律相同，即截面扁化率從夾塊端到壓塊端都是逐漸增加。隨著芯棒伸出量的增加，截面扁化率整體上呈減小的趨勢(shì)。這是因?yàn)樵黾有景羯斐隽烤驮黾恿斯懿膬?nèi)壁受到的支撐，從而有助于減小彎管截面扁化。

圖3 芯棒伸出量對(duì)壁厚減薄率、截面扁化率和回彈角的影響Fig.3 Effect of mandrel protrusion length on wall thinning ratio，ovalization ratio and springback angle

圖3c是芯棒伸出量對(duì)回彈角的影響圖。從圖中可以看出，隨著芯棒伸出量的增加回彈角呈減小趨勢(shì)。這是因?yàn)樾景羯斐隽吭龃螅箯澢瓮鈧?cè)受拉應(yīng)力增大，拉應(yīng)力可以減小回彈［8］。

2.3 壓塊相對(duì)助推速度的影響

圖4a為壓塊相對(duì)助推速度對(duì)壁厚減薄率的影響圖，可以看出，當(dāng)壓塊相對(duì)助推速度為0.9時(shí)，壁厚減薄率明顯大于壓塊相對(duì)助推速度大于等于1時(shí)的壁厚減薄率。這是因?yàn)楫?dāng)壓塊速度小于彎曲速度時(shí)，壓塊對(duì)管材的摩擦力與管材移動(dòng)方向相反，這樣就增大了管材外壁受到的拉應(yīng)力，使得外側(cè)壁厚減薄率增大。當(dāng)壓塊速度與彎曲速度相等或大于彎曲速度時(shí)，壓塊對(duì)管材向后的摩擦力消失，改善管材外壁受所受拉應(yīng)力起到側(cè)推作用，從而減小壁厚減薄。當(dāng)壓塊相對(duì)助推速度大于1時(shí)壁厚減薄率略有減小。圖4b為壓塊相對(duì)助推速度對(duì)截面扁化率的影響圖，從圖中可以看出當(dāng)壓塊相對(duì)助推速度為0.9時(shí)截面扁化率明顯增大，當(dāng)壓塊相對(duì)助推速度大于1時(shí)截面扁化率略微減小，且當(dāng)壓塊相對(duì)助推速度大于等于1時(shí)，截面扁化率變化不大。圖4c為壓塊相對(duì)助推速度對(duì)回彈角影響圖，從圖中可以看出，回彈角在壓塊相對(duì)助推速度為0.9時(shí)，回彈角最小，這是因?yàn)樾∮?的壓塊相對(duì)助推速度在管材彎曲過(guò)程中實(shí)際施加了向后的拉力，彎曲切向拉力可以減小回彈角。當(dāng)壓塊相對(duì)助推速度大于1時(shí)，回彈角隨壓塊相對(duì)助推速度的增大略微增大。

圖4 壓塊相對(duì)助推速度對(duì)壁厚減薄率、截面扁化率和回彈角的影響Fig.4 Effect of relative boosting velocity on wall thinning ratio，ovalization ratio and springback angle

2.4 相對(duì)彎曲半徑的影響

圖5a為相對(duì)彎曲半徑對(duì)壁厚減薄率的影響圖，從圖中可以看出隨著相對(duì)彎曲半徑的增大，壁厚減薄率逐漸減小。這是因?yàn)榇蟮南鄬?duì)彎曲半徑下，同樣彎曲角度下變形量相對(duì)較小，因此壁厚減薄率也較小。圖5b為相對(duì)彎曲半徑對(duì)截面扁化率的影響圖，可以看出隨著相對(duì)彎曲半徑的增大，截面扁化率逐漸減小。規(guī)律與壁厚減薄率變化一致。圖5c為相對(duì)彎曲半徑對(duì)回彈角的影響圖?？梢钥闯觯貜椊请S相對(duì)彎曲半徑的增大而增大，這是因?yàn)橄鄬?duì)彎曲半徑增大，同樣變形中的彈性變形增大，導(dǎo)致卸載后回彈角增大。

圖5 相對(duì)彎曲半徑對(duì)壁厚減薄率、截面扁化率和回彈角的影響Fig.5 Effect of relative bending radius on wall thinning ratio，ovalization ratio and springback angle

3 結(jié)論

(1)彎曲角度對(duì)Ti35合金管數(shù)控彎曲壁厚減薄率影響較小。隨著彎曲角度增大，截面扁化率增大，且回彈角與彎曲角度呈正比例增長(zhǎng)關(guān)系。

(2)Ti35合金管材彎曲變形越劇烈(如減小相對(duì)彎曲半徑、增大芯棒伸出量，或降低壓塊相對(duì)助推速度)，壁厚減薄率越大、回彈角越小。截面扁化率隨芯棒伸出量、相對(duì)彎曲半徑的增大而減小。當(dāng)壓塊相對(duì)助推速度小于1時(shí)，壁厚減薄率及截面扁化率明顯增大，但回彈角減小。

(3)通過(guò)有限元模擬，可以優(yōu)化Ti35合金管材數(shù)控彎曲成形工藝參數(shù)，減小試驗(yàn)次數(shù)，降低人力和物力成本，提高Ti35合金管材彎曲成形精度和效率。

［1］徐瀟瀟，邱邵宇.Ti合金在核燃料后處理設(shè)備中的應(yīng)用研究進(jìn)展［J］.材料導(dǎo)報(bào)，2005，19(11):57－59.

［2］Dey P K，Bansal N K.Spent fuel reprocessing:a vital link in Indian nuclear power program［J］.Nuclear Engineering and Design，2006，236(7/8):723－729.

［3］Baldev Raj，Kamachi Mudali U.Materials development and corrosion problems in nuclear fuel reprocessing plants［J］.Progress in Muclear Energy，2006，48(4):283 －313.

［4］吳金平，楊英麗，奚正平，等.Ti35合金焊接接頭在高溫硝酸中的腐蝕性能研究［J］.鈦工業(yè)進(jìn)展，2012，29(1):22－25.

［5］Jiang Z Q，Yang H，Zhan M.Establishment of a 3D FE model for the bending of a titanium alloy tube［J］.International Journal of Mechnical Sciences，2010，52(9):1115 －1124.

［6］Yang H，Li H，Zhang Z Y，et al.Advances and Trends on Tube bending Forming Technologies［J］.Chinese Journal of Aeronautics，2012，25(1):1 －12.

［7］Song F F，Yang H，Li H，et al.Springback prediction of thick-walled high-strength titanium tube bending［J］.Chinese Journal of Aeronautics，2013，26(5):1336 －1345.

［8］宋飛飛，楊合，楊英麗，等.鈦合金管冷彎回彈及其控制研究現(xiàn)狀［J］.鈦工業(yè)進(jìn)展，2014，31(6):7－11.

Finite Element Simulation of Forming Rules of Ti35 Alloy Tube Numerical Control Bending

Song Feifei，Wu Jinping，Yang Yingli，Zhao Bin，Luo Yuanyuan，Guo Dizi
(Northwest Institute for Nonferrous Metal Research，Xi’an 710016，China)

The numerical-control(NC)bending process of new corrosion-resistant Ti35 alloy tube is numerical analyzed using finite element(FE)software Abaqus.In the simulation，the effects of bending angle，mandrel protrusion length，relative boosting velocity and relative bending radius on forming process of Ti35 alloy tube were investigated.It is observed that the thickness thinning ratio is increasing，while the springback angle is decreasing when intensifying the bend forming，such as decrease the bending radius and relative boosting velocity，or increase the mandrel protrusion length.The mandrel protrusion length and relative bending radius have the negative effect.FE simulation can effectively decrease the manpower and material resource costs，and it also provides some reference for further industrialized application of Ti35 alloy.

Ti35 alloy;titanium alloy tube;NC bending;FE simulation

10.13567/j.cnki.issn1009-9964.2015.04.005

2015－05－07

陜西省重點(diǎn)科技創(chuàng)新團(tuán)隊(duì)計(jì)劃“鈦合金研發(fā)創(chuàng)新團(tuán)隊(duì)”(2012KCT-23)

宋飛飛(1987—)，男，助理工程師。