TIG焊熱源特性的反演分析

魏紅莉，于成奎，金偉紅

（徐州徐工施維英機械有限公司，江蘇徐州221004）

TIG焊熱源特性的反演分析

魏紅莉，于成奎，金偉紅

（徐州徐工施維英機械有限公司，江蘇徐州221004）

采用反演分析法，結(jié)合軟件ANSYS的參數(shù)化程序設(shè)計語言APDL，對TIG焊進行三維動態(tài)、非線性溫度場模擬，得到熔池長度、熔寬和熔深等參數(shù)，與焊接實驗得到的實際參數(shù)進行比較分析。利用直接反演分析方法，對比模擬得到的TIG焊接中熱源集中系數(shù)的取值范圍與由實驗確定的數(shù)值，來驗證實驗數(shù)據(jù)的準確性，為提高數(shù)值模擬精度，深入了解焊接過程熱、應(yīng)力本質(zhì)提供依據(jù)。

數(shù)值模擬；反演分析；焊接溫度場

0 前言

TIG焊接是近年來廣泛應(yīng)用的焊接方法，由于其焊接過程穩(wěn)定易于控制、可以自動清除工件氧化膜、適合全位置焊接以及焊接質(zhì)量高等特點在生產(chǎn)中大量應(yīng)用。焊接過程中熔池內(nèi)發(fā)生的各種復(fù)雜的物理化學(xué)變化對于焊縫成形和焊接接頭力學(xué)性能有著重要的影響，因此進一步了解焊接熔池的物理過程，能夠更好地實現(xiàn)焊接過程控制以及改進焊接方法。

TIG焊熱源特性的反演分析充分發(fā)揮了焊接過程數(shù)值模擬節(jié)省人力、物力和財力的巨大優(yōu)勢。TIG焊接過程其電弧加熱的熱效率η應(yīng)用范圍在文獻[3]有明確記載，本研究利用軟件ANSYS建立了TIG焊熔池的三維數(shù)學(xué)模型，利用APDL參數(shù)化程序設(shè)計語言進行數(shù)值模擬，模擬結(jié)果與實操實驗檢測結(jié)果進行對比性驗證，反演分析出TIG焊熱源半徑r及熱源集中系數(shù)K的取值范圍[1]。

1 焊接試驗

試樣材料為低碳鋼Q345，試件尺寸0.27 m× 0.10m×0.006m。TIG焊接參數(shù)為：焊接電流200 A，焊接電壓16V，焊接速度0.004m/s，材料的物理參數(shù)選擇為：板厚0.006 m，熔點1 530℃，密度7 750 kg/m3，熱導(dǎo)率25 W/m·k，比熱容460 J/（kg·k）[2]。按照上述焊接參數(shù)施焊，得到熔池長度、熔寬和熔深參數(shù)如表1所示。

表1 試驗實測數(shù)據(jù)

2 焊接過程數(shù)值模擬

根據(jù)焊接試驗中設(shè)定的相關(guān)參數(shù)及模擬條件，

結(jié)合實測所得數(shù)據(jù)，采用反演分析方法進行成因反推。在電弧電壓、焊接電流、焊接速度等參數(shù)保持相同的情況下，熱源效率選定在0.78～0.85范圍時進行數(shù)值模擬，通過不斷改變熱源有效半徑值r，得到熔池長、熔寬和熔深三個參數(shù)的模擬值，將模擬結(jié)果與實驗測得數(shù)值進行對比分析，逆向反推在誤差允許的范圍內(nèi)確定熱源有效半徑r，帶入公式K=3/r2，計算出熱源集中系數(shù)K[3]。

采用焊接試驗給定的數(shù)據(jù)進行APDL程序設(shè)計，模擬過程中采取非線性瞬態(tài)熱傳導(dǎo)控制工程進行焊接溫度場分析。分析前提中設(shè)定此焊接板材為各向同性材料，溫度T（x，y，z，τ）作為空間坐標（x，y，z）與時間τ的函數(shù)，在區(qū)域Ω中的任何點應(yīng)滿足導(dǎo)熱微分方程[4]

式中k為熱導(dǎo)率[單位：W/（m·k）]；c為比熱容[單位：J/（kg·k）]；ρ為材料密度（單位：kg/m3）；q為區(qū)域Ω中的熱源（單位：W/m3）；?2為拉普拉斯運算符號。

利用Galerkin法，將上式寫成有限元法的矩陣表達式，單元總體合成為

式中K為熱導(dǎo)矩陣；C為熱容矩陣；P為熱流列向量；T為溫度列向量。

熱源模型采用高斯分布，焊接電流、電弧電壓和熱效率已知，電弧中心處最大比熱流為

距電弧中心r處比熱流為

式中K為能量集中系數(shù)（熱源集中系數(shù)）（單位：cm-2）。

在單位時間步長l內(nèi)長線熱源上的比熱流

式中f為熱源作用面積；l為線熱源長度。l=焊接速度×?xí)r間步長=0.004 m/s×0.2 s=0.0008 m，焊接段狀熱源模型如圖1所示[5]。

圖1 焊接段狀熱源模型

將焊接對稱面考慮為絕熱邊界條件，試件上下表面及周圍邊界作換熱邊界條件處理，即

3 試驗分析

熱源效率η在0.78～0.85范圍內(nèi)取值，應(yīng)用軟件ANSYS的APDL語言編輯程序反復(fù)驗證得到熱源有效半徑r的區(qū)間范圍。以η=0.82，r=0.006 m處模擬得到的圖形為例進行分析，如圖2所示。

實驗材料為低碳鋼Q235，其熔點為1 530℃，將圖2b、圖2c與圖2d上溫度為1 530℃處的熔池長、熔深與熔寬計算列表如表2所示。

表2 ANSYS軟件模擬數(shù)據(jù)

由表2可知，在η=0.82，r=0.006 m時模擬得到的數(shù)值與實驗測量的數(shù)值誤差均在5％以內(nèi)，說明此種情況下熱源有效半徑的選擇符合實際要求。當

熱源效率在0.78～0.85之間不斷變化，得到熱源有效半徑的取值范圍為0.005～0.007 m。利用公式轉(zhuǎn)換得到熱源集中系數(shù)K的取值范圍為6.1～12.0。

圖2 η=0.82時的模擬計算結(jié)果

4 結(jié)論

（1）依據(jù)焊接試驗中設(shè)定的相關(guān)參數(shù)及模擬條件，結(jié)合實驗實測所得數(shù)據(jù)結(jié)果，采用反演分析方法進行成因反推。在熱源效率選定在0.78～0.85范圍內(nèi)進行數(shù)值模擬，通過不斷改變熱源有效半徑值，得到熔池長、熔寬和熔深三個參數(shù)的模擬值，將模擬結(jié)果與實驗測得的數(shù)值進行對比分析，逆向反推在誤差允許的范圍內(nèi)確定熱源有效半徑取值范圍為0.005～0.007 m。

（2）采用計算機焊接數(shù)值模擬手段支撐的反演分析方法，可以反算出許多焊接過程熱、力學(xué)本質(zhì)中涉及的不確定參數(shù)，是一種解決問題的有效手段。

[1]陶軍，李冬青，張紹娟，等.焊接過程參數(shù)反演分析進展[J].焊接，2005（9）：13-16.

[2]飯?zhí)镄抟?，大野和郎，澤田正三，?物理學(xué)常用數(shù)表[M].北京：科學(xué)出版社，1987.

[3]陳楚，汪建華，楊洪慶，等.數(shù)值分析在焊接中的應(yīng)用[M].上海：上海交通大學(xué)出版社，1985：140-144.

[4]汪建華.焊接數(shù)值模擬技術(shù)及其應(yīng)用[M].上海：上海交通大學(xué)出版社，2003.

[5]陳家權(quán)，肖順湖，楊新彥，等.焊接過程數(shù)值模擬熱源模型的研究進展[J].裝備制造技術(shù)，2005（3）：10-13.

Inversion analysis of TIG welding heat source characteristic

WEI Hongli，YU Chengkui，JIN Weihong
（Xuzhou XCMG Schwing Machinery Co.，Ltd，Xuzhou 221004，China）

Using inverse analysis,and combined with the software ANSYS parametric programming language APDL，three dimensional dynamic and nonlinear temperature field simulation of TIG welding is carried out，and the parameters such as length，width and depth of the weld pool are obtained，compared with the actual parameters of welding experiment.Using direct inverse analysis method，compare the scale of heat source focus modulus data in simulation count with the data in experiment，validate the experiment data's veracity，offer a foundation todeeplyresearch heat and stress ofwelding.

numerical simulation；inversion analysis；welding temperature field

TG409

A

1001－2303（2016）07－0134-03

10.7512/j.issn.1001-2303.2016.07.31

2016-03-31

魏紅莉（1979—），女，湖北襄陽人，碩士，主要從事材料學(xué)的相關(guān)研究工作。