李蘇杰, 楊　敏, 劉　昆, 李　星

(江蘇科技大學(xué) 船舶與海洋工程學(xué)院，江蘇 鎮(zhèn)江 212003)

?

基于LS-DYNA的船舶型材冷彎回彈研究

李蘇杰, 楊敏, 劉昆, 李星

(江蘇科技大學(xué) 船舶與海洋工程學(xué)院，江蘇 鎮(zhèn)江 212003)

摘要：冷彎型材成型過(guò)程非常復(fù)雜，回彈又是冷彎成型研究中重要部分，準(zhǔn)確預(yù)測(cè)回彈對(duì)提高加工效率和降低加工成本具有重要意義。以上海某船廠肋骨冷彎?rùn)C(jī)為模型，采用顯式動(dòng)力學(xué)有限元法對(duì)型材加載過(guò)程進(jìn)行了模擬計(jì)算，利用隱式靜力學(xué)有限元法對(duì)后續(xù)回彈進(jìn)行分析。計(jì)算表明，彎曲曲率與回彈曲率近似成線性關(guān)系，繼而對(duì)線性關(guān)系的系數(shù)進(jìn)行回歸分析，得到角鋼的回彈預(yù)測(cè)公式；最后將仿真結(jié)果與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行比較，結(jié)果顯示仿真結(jié)果具有一定的可信度。

關(guān)鍵詞：船舶型材；冷彎；回彈；LS-DYNA

0引言

回彈是船舶型材冷加工過(guò)程中存在的一種普遍現(xiàn)象，它直接影響著型材加工精度。獲得符合尺寸的型材的困難問(wèn)題之一就是回彈，如果不充分了解彎曲回彈規(guī)律，要想獲得預(yù)期精度的彎曲件，需要反復(fù)試彎，或者在成形結(jié)束后還要對(duì)彎曲件進(jìn)行校形,這樣將會(huì)導(dǎo)致生產(chǎn)周期的延長(zhǎng)和生產(chǎn)成本的提高，因此需要深入研究回彈問(wèn)題,從而很好地控制彎曲件的回彈。目前，關(guān)于型材彎曲回彈的有限元方法并不鮮見(jiàn)[1—3]，但專門(mén)針對(duì)船舶型材加工工藝特點(diǎn)并不多。

1塑性成型的有限元計(jì)算

上海某船廠冷加工車(chē)間數(shù)控肋骨冷彎?rùn)C(jī)如圖1所示。工作時(shí)，兩側(cè)夾頭模具上下施加載荷夾住型材腹板，約束腹板的6個(gè)方向的自由度，中間油壓沖頭模具相對(duì)肋骨向前、后運(yùn)動(dòng)，從而向肋骨上施加集中載荷，實(shí)現(xiàn)型材的正向或是反向彎曲。通過(guò)工人師傅的經(jīng)驗(yàn)，適當(dāng)過(guò)量彎曲肋骨，最后松開(kāi)夾頭模具和沖頭模具與型材之間的接觸，可以使肋骨回彈到指定的肋骨型線上。

根據(jù)上述分析的型材成型方式，應(yīng)用ANSYS/LS-DYNA建立了型材塑性加工成型有限元計(jì)算模型，如圖2所示。這里型材選用隱式-顯式可以相互轉(zhuǎn)換的單元shell163單元，且使用能容忍彎曲能力Belytschko-wong-chiang算法，沿厚道方向設(shè)置4個(gè)積分點(diǎn)，單元大小為10 mm×10 mm，采用映射劃分網(wǎng)格。肋骨冷彎?rùn)C(jī)模具變型可以忽略，選用剛體、隱式-顯式可以相互轉(zhuǎn)換的體單元solid164，單元大小為20 mm×20 mm×20 mm，采用掃掠劃分網(wǎng)格，全積分控制沙漏。回彈分析時(shí)，將彎曲階段的應(yīng)力場(chǎng)和位移場(chǎng)加到型材上，同時(shí)約束型材兩端面Y、Z方向位移，使型材在X方向回彈，這樣既保證回彈分析收斂，又可以使型材內(nèi)能充分釋放出來(lái)。

圖1　某船廠肋骨冷彎?rùn)C(jī)

Q235A的腹板高度×厚度為330 m×12 mm，面板寬度×厚度為120 mm×18 mm，在肋骨冷彎?rùn)C(jī)兩側(cè)夾模夾持力為400 kN、位移載荷為6 cm作用后，型材內(nèi)的應(yīng)力云圖如圖3(a)所示，卸載回彈后型材內(nèi)的殘余應(yīng)力云圖如圖3(b)所示。

圖3　加載應(yīng)力云圖與殘余應(yīng)力云圖

型材彎曲后的應(yīng)力分布、位移分布作為回彈的初始應(yīng)力載荷和初始位移載荷，直接決定回彈量的大小。從圖3(a)可以發(fā)現(xiàn)，在加載過(guò)程中，型材腹板下部受到最大的應(yīng)力，在加工過(guò)程中最容易受到損壞；型材兩端由于受到兩側(cè)夾模的影響，應(yīng)力分布則較為雜亂。型材回彈成型后留在型材內(nèi)的應(yīng)力稱為殘余應(yīng)力，殘余應(yīng)力的大小、分布直接影響回彈量的大小。殘余應(yīng)力的分布可以理解為回彈后留在型材內(nèi)的沒(méi)有釋放出來(lái)的內(nèi)能，分布的越是復(fù)雜，應(yīng)力越是較大，越能說(shuō)明在回彈過(guò)程中釋放的能量越小，回彈量也就越小。從圖3(b)可以發(fā)現(xiàn)，在型材中部和兩側(cè)型材內(nèi)的殘余應(yīng)力較大。

通過(guò)Ansys軟件處理后，可以提取被彎曲型材腹板外緣(因?yàn)楦拱逋饩夒x中間沖模相對(duì)較遠(yuǎn)，受集中載荷影響較小，回彈規(guī)律明顯便于研究)節(jié)點(diǎn)回彈前、后的位移信息，應(yīng)用MATLAB的gradient函數(shù)，求解型材腹板外緣離散節(jié)點(diǎn)曲率與曲率半徑。角鋼腹板外緣回彈前、后曲率曲線如圖4所示。

從圖4可以發(fā)現(xiàn)，在型材兩側(cè)相當(dāng)一段長(zhǎng)度，曲率幾乎為0，隨著向型材中部靠近型材所受彎矩慢慢變大，曲率出現(xiàn)了突然增大；另外，由于受到模具形狀和集中載荷的影響，型材中部附近出現(xiàn)了較大的曲率。從回彈前、后曲率比較來(lái)看，2條曲線在型材中部相差較大，說(shuō)明此位置回彈量也相對(duì)較大。

圖4　角鋼回彈前后X軸向位置-曲率圖

2回彈的回歸分析

圖5　不同鋼彎曲曲率-回彈曲率關(guān)系曲線

表1　Q235A的3組規(guī)格

從圖5可以發(fā)現(xiàn)，隨著彎曲曲率的增加，型材的回彈量也隨之增加。另外，越是高強(qiáng)度鋼，型材的回彈量越大。從圖6可以發(fā)現(xiàn)，隨著型材腹板高度與面板寬度比值的增加，型材的回彈量反而會(huì)減小，這是因?yàn)樾筒牡母拱宓幕貜椥∮谛筒牡拿姘?對(duì)整個(gè)型材的回彈起到限制作用。

圖6　不同h/b彎曲曲率-回彈曲率關(guān)系曲線

由數(shù)值分析的結(jié)果可得知，相對(duì)彎曲曲率與相對(duì)回彈曲率近似為線性關(guān)系,所以可以用一元線性建立其函數(shù)表達(dá)式，并根據(jù)不同模型得到不同的線性關(guān)系，最后將線性關(guān)系的系數(shù)進(jìn)行回歸分析。

預(yù)建立的一元線性方程：

kd=Ak+B

(1)式中：kb為相對(duì)回彈曲率；k為相對(duì)彎曲曲率；A、B為關(guān)系式系數(shù)，其值可采用最小二乘法對(duì)仿真得到的kb與k關(guān)系進(jìn)行線性擬合得到。

針對(duì)仿真回彈分析結(jié)果進(jìn)行回歸分析，需要建立多元回歸模型。建模時(shí)，需要選取對(duì)回彈影響較大的因子作為回歸分析因素，根據(jù)經(jīng)驗(yàn)公式[4]與仿真分析綜合考慮，選取型材的材料力學(xué)性能σS/E(屈服強(qiáng)度與彈性模量比值)、D/E(切線模量與彈性模量比值)，型材幾何參數(shù)h/b(型材高度與寬度比值)、η(型材橫截面形狀系數(shù))作為回歸分析因素。為了提高回歸分析的準(zhǔn)確性，選用了3組不同腹板高，3組不同面板寬，3組不同的材料的角鋼進(jìn)行交叉，共計(jì)15組不同的角鋼，數(shù)據(jù)見(jiàn)表2。

表2　角鋼回歸分析數(shù)據(jù)

將表1數(shù)據(jù)應(yīng)用回歸分析程序[5]，可得到系數(shù)A、系數(shù)B的回歸方程：

A= 161.233 2σS/E-4.025 1D/E-

0.006 1h/b-0.027 7η

B=-4.797 3σS/E+0.275 4D/E+0.000 3h/b+0.000 9η

將A、B帶入式(1)可得角鋼回彈的預(yù)測(cè)公式：

kb=(161.233 2)σS/E-4.025D/E-0.006 1h/b-0.027 7η)k+(-4.797 3σS/E+0.275 4D/E+

0.000 3h/b+0.000 9η)

應(yīng)用回歸方法得到的方程在實(shí)際應(yīng)用中有它特定的適用范圍?，F(xiàn)在所得到的回彈預(yù)測(cè)模型是比較簡(jiǎn)單化的數(shù)學(xué)模型，需要在以后的實(shí)際生產(chǎn)中去驗(yàn)證，并且根據(jù)實(shí)際情況來(lái)修正。

3仿真結(jié)果與實(shí)驗(yàn)結(jié)果的比較

本文將仿真的結(jié)果與薛江山[6]在武漢理工大學(xué)造船工藝研究所內(nèi)用100 t數(shù)控肋骨冷彎?rùn)C(jī)加工的肋骨試驗(yàn)測(cè)量值進(jìn)行比較，所用型材的腹板高度厚度×厚度為70 mm×5 mm，面板寬度×厚度為50 mm×5 mm，試驗(yàn)結(jié)果見(jiàn)表3。文獻(xiàn)中測(cè)量的是彎曲角和回彈角，由文獻(xiàn)[6]得知，彎曲前后曲率和彎曲角的變化規(guī)律是相同的。

從表3可以發(fā)現(xiàn)，仿真和試驗(yàn)測(cè)得的結(jié)果較為一致，但當(dāng)型材變形量較小時(shí)，計(jì)算值與試驗(yàn)值誤差較大，且前者較后者大。隨著型材變形量的不斷增大，仿真值也逐漸向試驗(yàn)值靠近，兩者之間誤差逐漸減小。

表3　試驗(yàn)與仿真回彈量對(duì)比

4結(jié)論

本文通過(guò)ANSYS/LS-DYNA建立肋骨冷彎?rùn)C(jī)有限元模型,分析了型材彎曲后的應(yīng)力分布、殘余應(yīng)力分布，研究了多因素對(duì)回彈的影響。根據(jù)計(jì)算結(jié)果可以得出以下結(jié)論：

(1)同樣情況下，高強(qiáng)度鋼回彈量要高于低強(qiáng)度鋼，h/d越大的型材回彈量越小。

(2)對(duì)回彈因素進(jìn)行回歸分析，得到角鋼的回彈預(yù)測(cè)公式。

(3)將仿真結(jié)果與實(shí)驗(yàn)值對(duì)比，結(jié)果表明兩者基本一致。但當(dāng)型材變形量較小時(shí)，計(jì)算值與試驗(yàn)值誤差相對(duì)較大。隨著型材變形量的不斷增大，仿真值也逐漸向試驗(yàn)值靠近。

參考文獻(xiàn)：

[1]時(shí)陽(yáng).L截面型鋼壓彎的有限元分析與回彈預(yù)測(cè)[D]. 大連：大連交通大學(xué),2013.

[2]谷瑞杰.彎曲成形回彈研究進(jìn)展[J].機(jī)械科學(xué)與技術(shù),2005(6)：653-658.

[3]Eggerlsen P A,Mattiasson K. Material modeling for accurate spring back prediction[J]. International Journal of Material Forming, 2009 (1):793-796.

[4]紀(jì)卓尚.船舶制造工藝力學(xué)[M].北京：國(guó)防工業(yè)出版社,2005.

[5]王巖,隋思漣,王愛(ài)青.數(shù)理統(tǒng)計(jì)與MATLAB工程數(shù)據(jù)分析[M].北京:清華大學(xué)出版社,2006.

[6]薛江山.肋骨冷彎自動(dòng)加工回彈問(wèn)題及成形精度的實(shí)驗(yàn)研究[D]. 武漢：武漢理工大學(xué),2007.

中圖分類號(hào)：U661.43

文獻(xiàn)標(biāo)志碼：A

作者簡(jiǎn)介：李蘇杰(1986—)，男，碩士研究生，研究方向?yàn)榇芭c海洋結(jié)構(gòu)物設(shè)計(jì)制造；楊敏(1963—)，女，副教授，研究方向?yàn)榇霸O(shè)計(jì)與先進(jìn)制造。

收稿日期：2015-11-10