張雪彪，鹿想，劉海涵

(大連理工大學(xué) a.運(yùn)載工程與力學(xué)學(xué)部船舶工程學(xué)院；b.工業(yè)裝備結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧 大連 116024)

在船廠，帆形板的水火彎板成形多采用三角加熱工藝，即火炬在靠近板邊時(shí)進(jìn)行擺動(dòng)加熱，從而增大板邊的加熱面積和溫度，以得到大的板邊收縮量，達(dá)到帆形板的成型要求。對(duì)于感應(yīng)熱源的水火彎板工藝，見圖1，當(dāng)感應(yīng)器在板邊處進(jìn)行擺動(dòng)加熱時(shí)，也可以獲得同樣的收邊加熱效果。

圖1 帆形板變形過程

對(duì)于鋼板多加熱線的變形計(jì)算問題，采用電磁-熱-結(jié)構(gòu)耦合的有限元模型進(jìn)行計(jì)算是困難的，因?yàn)槟Ｐ蛷?fù)雜，計(jì)算量巨大，所以研究者多采用基于固有應(yīng)變法的彈性有限元模型進(jìn)行計(jì)算。有學(xué)者等使用圓盤-彈簧模型作為三角加熱過程的簡(jiǎn)化熱彈塑性力學(xué)模型，在三角加熱的固有應(yīng)變區(qū)域內(nèi)進(jìn)行分區(qū)域的等效載荷加載[1]。有學(xué)者提出計(jì)算三角加熱形狀和位置的算法，建立加熱形狀和殘余變形的函數(shù)關(guān)系。相關(guān)的研究均采用火焰熱源的圓盤-彈簧模型來計(jì)算固有應(yīng)變，力學(xué)模型不適用于除圓形外的其他感應(yīng)器形狀，也不能簡(jiǎn)單地應(yīng)用于鋼板感應(yīng)加熱工藝。為此，考慮采用鋼板電磁-熱耦合的熱彈塑性模型得到加熱線附近的塑性應(yīng)變值，然后對(duì)塑性應(yīng)變進(jìn)行積分，從而獲得固有應(yīng)變沿加熱線的分布規(guī)律，在進(jìn)行彈性有限元分析時(shí)考慮固有應(yīng)變的分布規(guī)律，使之符合鋼板三角感應(yīng)加熱的變形特點(diǎn)。

1 加熱實(shí)驗(yàn)與數(shù)值模擬

針對(duì)鋼板的三角感應(yīng)加熱工藝，文獻(xiàn)[3]在大量實(shí)驗(yàn)的基礎(chǔ)上確定了三角感應(yīng)加熱的擺動(dòng)加熱路徑，并建立了基于輥彎板的鋼板三角感應(yīng)加熱的電磁熱耦合模型，該數(shù)值模型為進(jìn)一步應(yīng)用固有應(yīng)變法預(yù)測(cè)鋼板變形提供了條件。

1.1 實(shí)驗(yàn)布置

在大連理工大學(xué)造船工藝實(shí)驗(yàn)室進(jìn)行加熱實(shí)驗(yàn)，輥彎板的曲率半徑為2 000 mm，鋼板尺寸及加熱的工藝參數(shù)見表1，實(shí)驗(yàn)方案布置見圖2，距離加熱線中間偏右側(cè)20 mm處的P點(diǎn)為溫度測(cè)量點(diǎn)，通過測(cè)量加熱線兩側(cè)11對(duì)測(cè)量點(diǎn)的距離變化得到加熱前后的橫向收縮量。

表1 感應(yīng)加熱的工藝參數(shù)

圖2 實(shí)驗(yàn)布置和測(cè)量方案

三角感應(yīng)加熱路徑見圖3，即先線狀加熱而后折線加熱，從而實(shí)現(xiàn)三角加熱的效果。直線部分210 mm，速度v1為8 mm/s，折線加熱區(qū)域的長(zhǎng)度90 mm，速度v2為9 mm/s。加熱過程中的局部現(xiàn)場(chǎng)見圖4，加熱結(jié)束后的鋼板表面見圖5，在靠近板邊處的加熱線區(qū)域大致呈三角形狀。

圖3 三角感應(yīng)加熱路徑示意

圖4 三角感應(yīng)加熱圖

圖5 三角感應(yīng)加熱后鋼板表面

1.2 數(shù)值模擬

利用文獻(xiàn)[3-4]的輥彎板移動(dòng)式感應(yīng)加熱電磁-熱耦合有限元模型對(duì)加熱過程進(jìn)行模擬。輥彎板的移動(dòng)式感應(yīng)加熱示意見圖6，移除空氣包后的橫截面模型見圖7。

圖6 輥彎板的擺動(dòng)感應(yīng)加熱示意

圖7 鋼板-感應(yīng)器-空氣局部網(wǎng)格

整個(gè)加熱系統(tǒng)分為兩部分：①由感應(yīng)器、遠(yuǎn)場(chǎng)空氣包，以及上層運(yùn)動(dòng)空氣組成的移動(dòng)組件；②由鋼板和下層靜空氣組成的靜止組件。該模型通過節(jié)點(diǎn)約束方程技術(shù)，實(shí)現(xiàn)了感應(yīng)器、空氣間隙和鋼板之間相對(duì)移動(dòng)的計(jì)算域求解過程，計(jì)算鋼板的溫度和變形，揭示鋼板和感應(yīng)器之間的電磁耦合效應(yīng)。

P點(diǎn)溫度實(shí)驗(yàn)值和數(shù)值計(jì)算值的比較見圖8，橫向收縮的實(shí)驗(yàn)值與數(shù)值計(jì)算結(jié)果的對(duì)比見圖9，兩者的平均相對(duì)誤差為15.3%。橫向收縮和溫度的變化情況都與實(shí)驗(yàn)值一致，橫向收縮量的分布從板里到板邊呈逐漸增大的趨勢(shì)，表明所設(shè)計(jì)的三角感應(yīng)加熱的路徑是合理的。

圖8 P點(diǎn)溫度的數(shù)值實(shí)驗(yàn)對(duì)比

圖9 橫向收縮量比較

2 固有應(yīng)變法應(yīng)用

根據(jù)固有應(yīng)變法，焊接過程中引起的無彈性的應(yīng)變是固有應(yīng)變?chǔ)?[5]，即

ε*=ε-εe

(1)

式中：ε為總應(yīng)變量；εe為彈性應(yīng)變。

固有應(yīng)變是塑性應(yīng)變?chǔ)舙、熱應(yīng)變?chǔ)舤和相應(yīng)變?chǔ)纽热咧停?/p>

ε*=εp+εt+εθ

(2)

在實(shí)際加熱過程結(jié)束后，熱應(yīng)變?yōu)榱?，如不考慮相變應(yīng)變，固有應(yīng)變就是塑性應(yīng)變值。為了得到鋼板在加熱過程中的固有應(yīng)變，將塑性應(yīng)變?cè)跈M截面上積分得到固有應(yīng)變之和，從而得到固有應(yīng)變值[6]。

2.1 鋼板三角感應(yīng)加熱固有應(yīng)變的計(jì)算

對(duì)于輥彎板的三角感應(yīng)加熱過程，對(duì)數(shù)值模擬得到的各個(gè)截面的塑性應(yīng)變?chǔ)舙進(jìn)行積分，得到橫向固有應(yīng)變之和Wx和縱向固有應(yīng)變之和Wy。

(3)

(4)

式中：εx為橫向固有應(yīng)變；εy為縱向固有應(yīng)變；A為固有應(yīng)變區(qū)域沿垂直于加熱線方向的截面積。

對(duì)于殼單元或?qū)嶓w單元模型，在加載區(qū)域?qū)挾葹閎，加載深度為h的情況下，定義加載區(qū)域的材料具有x、y方向的熱膨脹系數(shù)，熱膨脹系數(shù)值αx、αy分別為

(5)

(6)

在進(jìn)行固有應(yīng)變加載時(shí)使用溫度載荷法進(jìn)行加載，在ANSYS軟件中設(shè)置材料的各向異性熱膨脹系數(shù)，然后施加單位溫度載荷，進(jìn)行一次彈性計(jì)算獲得變形結(jié)果。溫度與應(yīng)變的關(guān)系為

ε=α·Δt

(7)

式中:α為材料熱膨脹系數(shù);Δt為單位溫度載荷，設(shè)置為1 ℃;ε為產(chǎn)生的熱應(yīng)變數(shù)值。

沿著加熱線每隔15 mm選取一個(gè)截面，共20個(gè)截面，提取每個(gè)截面的節(jié)點(diǎn)塑性應(yīng)變值，積分后得到橫向和縱向固有應(yīng)變積分曲線見圖10。

圖10 固有應(yīng)變積分

可以看出橫向和縱向固有應(yīng)變積分曲線沒有明顯的變化規(guī)律，因此加載固有應(yīng)變時(shí)必須考慮其實(shí)際的分布規(guī)律。這一點(diǎn)不同于焊接產(chǎn)生的固有應(yīng)變，因?yàn)楹缚p長(zhǎng)度足夠長(zhǎng)，超過臨界焊接長(zhǎng)度后固有應(yīng)變有穩(wěn)定值[7]，可以用平均固有應(yīng)變來進(jìn)行加載。

首先加熱線長(zhǎng)度是有限的；另外，三角加熱的加熱路徑?jīng)Q定了鋼板加熱后產(chǎn)生非均勻的固有變形分布。鑒于此，在彈性有限元計(jì)算時(shí)須加載實(shí)際的固有應(yīng)變分布，而不能采用平均固有應(yīng)變。

2.2 固有應(yīng)變的加載

各截面處的塑性應(yīng)變?cè)茍D見圖11。在板邊Y=0 mm處，固有應(yīng)變集中分布區(qū)域?qū)挾茸畲?09.0 mm，Y=300 mm處寬度為44.2 mm，寬度值從板里向板邊逐漸增大。

圖11 各截面橫向塑性應(yīng)變?cè)茍D

鋼板彈性有限元分析模型見圖12。

圖12 鋼板整體網(wǎng)格圖

單元類型是SOLID 45單元。在加熱線附近的單元進(jìn)行了網(wǎng)格細(xì)化，網(wǎng)格尺寸為1 mm，其他非加熱區(qū)域的網(wǎng)格為5 mm。在彈性有限元計(jì)算時(shí)采取圖13所示的加載方法，在每個(gè)加載段內(nèi)分別加載不同的熱膨脹系數(shù)，經(jīng)過多次計(jì)算，發(fā)現(xiàn)梯形加載和矩形加載的變形結(jié)果差別很小，而矩形加載卻更方便，這與文獻(xiàn)[8]結(jié)論一致。因此，將固有應(yīng)變區(qū)寬度設(shè)置為50 mm。根據(jù)文獻(xiàn)[9]推導(dǎo)的考慮角變形的固有應(yīng)變加載深度公式，此處加載深度為14 mm。

圖13 固有應(yīng)變加載方法

圖14a)為橫向變形云圖，圖14b)為垂向變形云圖。彈性有限元計(jì)算的結(jié)果與熱彈塑性有限元模型的結(jié)果十分相近，計(jì)算時(shí)間約1 min，變形形狀也符合三角加熱的變形特點(diǎn)。

圖14 彈性有限元分析的位移云圖

橫向收縮量的對(duì)比見圖15，與實(shí)驗(yàn)值相比較彈性有限元分析結(jié)果相對(duì)誤差為13.84%。加熱線中心線處的垂向位移對(duì)比見圖16，兩者的相對(duì)誤差為3.75%。面內(nèi)變形和面外變形結(jié)果都符合工程要求，表明計(jì)算方法是可靠的。

圖15 橫向收縮量比較

圖16 垂向位移比較

3 輥彎板三角感應(yīng)加熱變形計(jì)算

帆形板成形時(shí)需要在鋼板上布置加熱線，基于固有應(yīng)變法的彈性有限元分析方法計(jì)算鋼板的變形。因?yàn)殇摪灞砻嬗卸鄺l加熱線，在忽略加熱線彼此之間影響的前提下進(jìn)行變形預(yù)測(cè)，實(shí)驗(yàn)時(shí)設(shè)置加熱線間距為200 mm。

實(shí)驗(yàn)用鋼板尺寸為1 200×800×20 mm，輥彎半徑2 000 mm，實(shí)驗(yàn)方案見圖17。

圖17 輥彎板三角感應(yīng)實(shí)驗(yàn)布置

采用靠近板邊處進(jìn)行擺動(dòng)加熱的三角加熱方式，共10條加熱線，加熱線的速度信息見表2，加熱線長(zhǎng)度300 mm，空氣間隙9 mm，加熱功率50 kW。鋼板在全部加熱線完成之后的加熱效果見圖18。

表2 加熱線速度信息 mm/s

圖18 加熱后的鋼板

輥彎板的彈性有限元模型見圖19，單元類型為SOLID45，均勻劃分鋼板網(wǎng)格，長(zhǎng)度、寬度方向的網(wǎng)格尺寸均為10 mm，板厚方向的網(wǎng)格尺寸為0.5 mm。使用溫度載荷法進(jìn)行固有應(yīng)變的加載，在各個(gè)加熱線處的單元施加X、Y方向的熱膨脹系數(shù)，限制剛體位移約束四個(gè)頂點(diǎn)的X、Y、Z方向的自由度。

圖19 鋼板有限元模型

測(cè)量線示意圖見圖20。鋼板曲面的測(cè)量采用徠卡掌上全站儀S910進(jìn)行測(cè)量，測(cè)量精度1.0 mm，測(cè)量現(xiàn)場(chǎng)見圖21，測(cè)量曲板坐標(biāo)點(diǎn)數(shù)據(jù)后進(jìn)行數(shù)據(jù)處理得到撓度值。鋼板的垂向位移云圖見圖22，鋼板中間區(qū)域的垂向位移較大。橫向變形云圖見圖23，收縮量在板邊處較大，鋼板整體呈現(xiàn)為帆形。

圖20 測(cè)量線布置示意

圖21 徠卡掌上全站儀S910測(cè)量鋼板曲面現(xiàn)場(chǎng)

圖22 垂向變形云圖

圖23 橫向變形云圖

各測(cè)量線的撓度測(cè)量值和計(jì)算值的對(duì)比見圖24，各計(jì)算值和實(shí)驗(yàn)值一致。其中A1線與A5線位于縱向板邊位置，因?yàn)榘暹吿幍墓逃袘?yīng)變值較大，計(jì)算出的垂向位移有波浪型的變化。A1～A5線的實(shí)驗(yàn)值和計(jì)算值的平均相對(duì)誤差分別為：14.2%，10.7%，7.8%，13.4%，13.1%，誤差均在工程允許的范圍內(nèi)。

圖24 測(cè)量線的撓度值比較

對(duì)于多條加熱線的輥彎板變形計(jì)算，在忽略加熱線之間影響的前提下，基于固有應(yīng)變法的彈性有限分析可以實(shí)現(xiàn)鋼板的整體變形計(jì)算，計(jì)算時(shí)間短，可以用于水火彎板工藝的整體變形預(yù)測(cè)。

4 結(jié)論

鋼板的三角感應(yīng)加熱工藝產(chǎn)生復(fù)雜的固有應(yīng)變分布規(guī)律，在進(jìn)行固有應(yīng)變的彈性有限元計(jì)算時(shí)在加載區(qū)域內(nèi)需要加載實(shí)際的固有應(yīng)變值，而不是平均固有應(yīng)變。同時(shí)，進(jìn)行輥彎板的多條加熱線的感應(yīng)加熱實(shí)驗(yàn)，使用建立的固有應(yīng)變加載模型進(jìn)行彈性有限元計(jì)算，計(jì)算的鋼板撓度值與實(shí)驗(yàn)值相符，整體變形效果與實(shí)驗(yàn)一致。