魏紅莉，于成奎，金偉紅

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

TIG焊熱源特性的反演分析

魏紅莉，于成奎，金偉紅

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

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

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

0 前言

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

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

1 焊接試驗(yàn)

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

表1 試驗(yàn)實(shí)測數(shù)據(jù)

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

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

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

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

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

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

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

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

距電弧中心r處比熱流為

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

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

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

圖1 焊接段狀熱源模型

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

3 試驗(yàn)分析

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

實(shí)驗(yàn)材料為低碳鋼Q235，其熔點(diǎn)為1 530℃，將圖2b、圖2c與圖2d上溫度為1 530℃處的熔池長、熔深與熔寬計(jì)算列表如表2所示。

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

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

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

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

4 結(jié)論

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

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

[1]陶軍，李冬青，張紹娟，等.焊接過程參數(shù)反演分析進(jìn)展[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ìn)展[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)研究工作。