楊瑞康，張勤練，周 舟，董 飛

(1.西安航天發(fā)動(dòng)機(jī)有限公司，陜西 西安 710100; 2.93956部隊(duì)，甘肅 張掖 361023)

0 引言

國(guó)內(nèi)外在汽車、航空、航天領(lǐng)域逐漸實(shí)現(xiàn)了機(jī)器人自動(dòng)焊接，但國(guó)內(nèi)在液體火箭發(fā)動(dòng)機(jī)推力室上機(jī)器人焊接技術(shù)還應(yīng)用得很少，焊接過(guò)程中參數(shù)控制技術(shù)應(yīng)用得更少。焊接控制技術(shù)應(yīng)用難度一般較大，主要因?yàn)楹附邮悄：茖W(xué)，焊接參數(shù)一般在一定范圍內(nèi)變動(dòng)，同時(shí)焊接質(zhì)量精準(zhǔn)評(píng)價(jià)困難[1]?，F(xiàn)有產(chǎn)品自動(dòng)焊接時(shí)因?yàn)楫a(chǎn)品坡口狀態(tài)的差異性容易產(chǎn)生局部咬邊，余高不均勻，或者未溶透氣孔等問(wèn)題，需要手工再補(bǔ)焊[2]。歐洲、美國(guó)NASA、SPACEX等航天企業(yè)逐漸應(yīng)用了焊接過(guò)程控制技術(shù)，國(guó)內(nèi)航天發(fā)動(dòng)機(jī)領(lǐng)域焊接過(guò)程自適應(yīng)控制還沒(méi)有實(shí)際的應(yīng)用[3-5]。

航天發(fā)動(dòng)機(jī)產(chǎn)品焊接質(zhì)量要求高，由于鈑金件熱加工帶來(lái)了坡口狀態(tài)不穩(wěn)定、存在間隙等問(wèn)題。本文結(jié)合試驗(yàn)與數(shù)值模擬技術(shù)，建立了推力室外壁坡口自適應(yīng)焊接參數(shù)模型，輸入為對(duì)接間隙、坡口寬度，輸出為焊接電流、送絲速度與擺動(dòng)幅度等工藝參數(shù)。并在SIMULINK平臺(tái)下建立了機(jī)器人自適應(yīng)仿真系統(tǒng)，以坡口間隙與坡口寬度為輸入量、焊接電流及焊接速度為輸出控制量，進(jìn)行了機(jī)器人焊接過(guò)程自適應(yīng)控制系統(tǒng)仿真驗(yàn)證。

1 試驗(yàn)設(shè)備及軟件

試驗(yàn)設(shè)備為機(jī)器人熱絲TIG焊接系統(tǒng)及SERVE激光焊縫坡口識(shí)別系統(tǒng)。相對(duì)于傳統(tǒng)TIG焊，熱絲TIG焊是將填充焊絲進(jìn)行焊前加熱，大大加快了焊絲的熔化速度(增加達(dá)20～50 g/min)[6]。

本文采用的SYSWELD焊接有限元軟件可用于熱溫度場(chǎng)模擬、電磁分析、力學(xué)分析和氫擴(kuò)散，同時(shí)焊接過(guò)程數(shù)值模擬采用了生死單元技術(shù)[7-8]。焊接工藝參數(shù)的調(diào)整過(guò)程與電弧特性、焊接電源特性等多方面因素相關(guān)，為了控制精確且響應(yīng)快速，本文通過(guò)MATLAB SIMULINK平臺(tái)搭建仿真系統(tǒng)進(jìn)行研究與驗(yàn)證。

圖1 機(jī)器人熱絲TIG焊接系統(tǒng)Fig.1 Robot heat TIG welding system

2 推力室外壁焊接參數(shù)模型建立

2.1 試板焊接試驗(yàn)

進(jìn)行步長(zhǎng)為5 A的變焊接參數(shù)熱絲TIG焊接試驗(yàn)，記錄焊接開(kāi)始出現(xiàn)缺陷的焊接參數(shù)。試驗(yàn)直到未熔透缺陷如圖2(a)所示，記錄此時(shí)的焊接電流Imin，然后將焊接電流不斷調(diào)大，直到焊漏缺陷如圖2(c)所示，記錄此時(shí)的焊接電流Imax。定義初選參數(shù)Im=0.5(Imin+Imax)。

圖2 焊接熔合狀態(tài)Fig.2 Fusion state of welding

同時(shí)根據(jù)焊接背部熔寬及焊接熱輸入最終確定優(yōu)選焊接參數(shù)。優(yōu)化參數(shù)需要滿足焊縫背面熔寬t最大、熱輸入W較小，即i={i|ti=maxt}，i={i|Wi=minW}。

2.2 焊接參數(shù)模型建立及數(shù)值模擬

本文通過(guò)數(shù)值模擬的方式得到熔池面積與電流參數(shù)、板厚及焊接速度的關(guān)系，在此之前需要調(diào)整數(shù)值模擬模型，使得數(shù)值模擬的結(jié)果與試驗(yàn)實(shí)測(cè)的結(jié)果保持一致性，從而準(zhǔn)確地進(jìn)行焊接過(guò)程數(shù)值模擬計(jì)算[9-10]。推力室材料有QCr0.8合金與S-06不銹鋼，通過(guò)對(duì)SYSWELD進(jìn)行二次開(kāi)發(fā)，建立材料QCr0.8與材料S-06的材料文件。

熱絲TIG焊過(guò)程中熔池尾端隨著兩側(cè)液體金屬的凝固逐漸變細(xì)，熔池頭部則是半圓形向前擴(kuò)展，根據(jù)對(duì)于熱源模型的研究，其中雙橢球移動(dòng)熱源模型用來(lái)描述這種狀態(tài)[11]，如圖3所示。

圖3 雙橢球熱源模型Fig.3 Double ellipsoid heat source model

三維非線性瞬態(tài)傳熱問(wèn)題的控制方程為：

(1)

(2)

焊接數(shù)值模擬過(guò)程包括以下步驟：幾何模型建立、網(wǎng)格劃分、材料物性參數(shù)整理與設(shè)置、熱源校核、邊界條件及約束條件設(shè)置、求解方法與參數(shù)設(shè)置、計(jì)算溫度場(chǎng)、計(jì)算應(yīng)力場(chǎng)及后處理[14]。提取計(jì)算結(jié)果，如圖4所示。

圖4 焊接試驗(yàn)試板溫度場(chǎng)云圖Fig.4 Temperature field of welding test plate

本文為了簡(jiǎn)化對(duì)焊接焊縫尺寸形狀的描述，約定：焊縫在垂直于焊接方向的穩(wěn)定的橫截面簡(jiǎn)稱為焊縫截面，其對(duì)應(yīng)的截面積為焊縫截面積。

設(shè)計(jì)三因素四水平回歸試驗(yàn)，探究焊縫截面積與焊接電流、焊接壁厚及焊接速度的關(guān)系模型。針對(duì)V型60°坡口數(shù)值模擬試驗(yàn)，變量為焊接電流100～400 A(打底焊接電流為50～300 A)、板厚5～12 mm與焊接速度10～25 cm/min。數(shù)值模擬試驗(yàn)結(jié)束后，設(shè)定材料熔化溫度線(固相線、液相線)后進(jìn)行后處理。并對(duì)熔池區(qū)域面積進(jìn)行計(jì)算測(cè)量，本文對(duì)熔池區(qū)域面積測(cè)量時(shí)定義3個(gè)測(cè)量值：高度X1，寬度X2及深度X3，如圖5所示。

圖5 焊縫截面區(qū)域測(cè)量圖Fig.5 Measurement of welding section

通過(guò)板件焊接數(shù)值模擬進(jìn)行各個(gè)參數(shù)的焊接過(guò)程計(jì)算，得到數(shù)據(jù)結(jié)果如表1所示。

查閱國(guó)內(nèi)外文獻(xiàn)在描述多變量單響應(yīng)量時(shí)，通常的方法是建立式(3)所示模型，并根據(jù)相關(guān)檢驗(yàn)法對(duì)回歸模型進(jìn)一步進(jìn)行檢驗(yàn)調(diào)整：

(3)

式中：w為響應(yīng)變量(焊縫截面積)；β為回歸系數(shù)；v為焊接工藝參數(shù)；ε為殘差;k=2。

表1 填充層數(shù)值模擬回歸試驗(yàn)結(jié)果

將表1數(shù)據(jù)處理求出定義面積后，進(jìn)行偏最小二乘回歸分析，經(jīng)過(guò)多次回歸，剔除影響過(guò)小的二次項(xiàng)、交叉項(xiàng)后得到填充層焊縫截面積與3個(gè)參數(shù)的回歸公式:

(4)

對(duì)上述得到的回歸方程模型進(jìn)行顯著性檢驗(yàn)，可采用F檢驗(yàn)法檢驗(yàn)回歸的總體效果。其F值為52.62，遠(yuǎn)大于F0.01(8，9)=5.47，回歸效果顯著；回歸方程相關(guān)系數(shù)R2=98.33%，也表明回歸效果非常顯著。同樣得到打底層焊縫截面積模型與蓋面層焊縫截面積模型：

(5)

(6)

以上3個(gè)模型均是根據(jù)S-06鋼物性參數(shù)所得，1Cr21Ni5Ti的計(jì)算焊縫截面積為S-06鋼的1.091倍。根據(jù)3種近似焊縫截面幾何模型圓、橢圓和梯形，計(jì)算得出圓、橢圓的弦長(zhǎng)等信息，進(jìn)而獲得模型規(guī)劃焊縫截面積S0。如圖6所示，根據(jù)焊接填充入金屬量及形狀(梯形、圓、橢圓)等建立坐標(biāo)系，計(jì)算出橢圓(圓)方程，得出規(guī)劃焊縫截面寬度、深度(最小的熔透所需熔池區(qū)域)所需的焊縫截面積。

圖6 預(yù)估焊縫截面形狀計(jì)算模型示意圖Fig.6 Calculation model of estimating welding section shape

焊縫截面按照橢圓與半圓/梯形計(jì)算如圖6所示。(x1，y1)，(x2，y2)，(x3，y3)及(x4，y4)為填充規(guī)劃面積后的層間高度與坡口側(cè)壁交點(diǎn)。其中打底層的梯形焊縫截面區(qū)域底部寬設(shè)為e。焊縫截面規(guī)劃面積S0只是焊縫截面積區(qū)域規(guī)劃的最小值，為保證焊接參數(shù)的可靠性及穩(wěn)定性，需對(duì)S0公式進(jìn)行修正，通過(guò)乘上一個(gè)大于1修正系數(shù)實(shí)現(xiàn)。此修正系數(shù)通過(guò)試驗(yàn)擬合得出。定義焊縫熔寬、深熔合修正系數(shù)為α，β。通過(guò)數(shù)據(jù)擬合得到熔合修正系數(shù)，填充：α=1.132 12，β=1.095 6；蓋面：α=1.123 8，β=1.061 2；打底：α=1.084 4，β=1。參數(shù)計(jì)算模型流程如圖7所示。

圖7 自適應(yīng)焊接參數(shù)模型框圖Fig.7 Block diagram of adaptive welding parameter

基于之前的試驗(yàn)研究，推力室及其工裝質(zhì)量很大導(dǎo)致轉(zhuǎn)動(dòng)慣性較大，難以精確地控制轉(zhuǎn)動(dòng)速度，故在焊接過(guò)程中固定15 cm/min 為常用焊接速度。確定焊接電流、送絲速度與擺幅三者為焊接自適應(yīng)調(diào)整參數(shù)，其三者分別可以控制熱輸入、坡口填入金屬量與熱源分布，滿足焊接過(guò)程中需要調(diào)整的3個(gè)方面。鎖底坡口送絲量調(diào)整示意圖如圖8所示，圖b中為間隙。

(a)打底層 (b)填充層 (c)蓋面層圖8 鎖底坡口送絲量調(diào)整示意圖Fig.8 Adjustment of wire feeding amount in the groove

打底層焊接送絲速度(cm/min)求解：

(7)

打底層焊接電流模型求解：

(8)

(9)

IW=182+23.532 2b+2.117 9b2

(10)

同樣得到填充層及蓋面層模型：

填充層模型

SW=390+48.412 1b

(11)

IW=283+2.224 4b2+46.089 7b

(12)

蓋面層模型

SW=370+39.125b

(13)

IW=274+51.251 6b+5.914 6b2

(14)

As=±(2.5+0.6b)

(15)

2.3 自適應(yīng)焊接控制系統(tǒng)仿真

傳統(tǒng)上，焊接電源采用簡(jiǎn)單下降的平特性形式， 這種電源簡(jiǎn)化為一階慣性環(huán)節(jié)，即電源電弧系統(tǒng)傳遞函數(shù)與伺服電機(jī)驅(qū)動(dòng)的送絲機(jī)構(gòu)傳遞函數(shù)：

(16)

其中

Thd=0.001 2

Khd=22.5

Tm=0.001s

Ts1=0.006

式中Kss為Ud到vs的靜態(tài)倍率，約為10 mm/S·V。

利用MATLAB SIMULINK仿真平臺(tái)構(gòu)建模糊控制模塊fuzzy-1，將自適應(yīng)參數(shù)模型嵌入到仿真系統(tǒng)中，進(jìn)一步建立了焊接電流與送絲速度的自適應(yīng)控制系統(tǒng)如圖9所示，并進(jìn)行了仿真，得到了優(yōu)化的各個(gè)比例因子的取值。焊接電流控制模糊控制系統(tǒng)的最佳比例因子K3=50；調(diào)整因子α1=3，α2=0.2；送絲速度控制模糊控制系統(tǒng)的最佳比例因子K3=800；調(diào)整因子α1=4，α2=0.5。

進(jìn)一步通過(guò)使用sin(x)函數(shù)作為輸入量模擬坡口間隙的變化，進(jìn)行焊接電流及送絲速度自適應(yīng)模糊控制算法仿真。如圖10所示，該自適應(yīng)控制系統(tǒng)控制準(zhǔn)確性良好，響應(yīng)速度小于80 ms且無(wú)超調(diào)現(xiàn)象，可實(shí)現(xiàn)自適應(yīng)焊接參數(shù)快速響應(yīng)控制。

3 試驗(yàn)驗(yàn)證

本文試驗(yàn)驗(yàn)證了推力室擴(kuò)張段樣機(jī)及其坡口試驗(yàn)件焊接，焊接結(jié)束后，產(chǎn)品焊縫焊接外觀良好，無(wú)咬邊、焊漏等缺陷。

圖9 焊接自適應(yīng)控制系統(tǒng)Fig.9 Welding adaptive control system

圖10 焊接自適應(yīng)系統(tǒng)仿真曲線Fig.10 Simulation curve of welding adaptive system

如圖11所示，焊縫余高分別為1.85 mm，1.89 mm及1.88 mm，誤差為2%。本試驗(yàn)在推力室鎖底坡口寬度波動(dòng)最大值為1.3 mm情況下焊接一致性良好。進(jìn)行了推力室坡口的模擬件試驗(yàn)，焊縫熔合良好，并進(jìn)行X光檢測(cè)合格。對(duì)焊縫拉伸件進(jìn)行了拉伸試驗(yàn)，兩試件斷裂位置均發(fā)生在熱影響區(qū)，強(qiáng)度均達(dá)到母材強(qiáng)度的90%。

圖11 推力室坡口試板自適應(yīng)焊接試驗(yàn)圖Fig.11 Adaptive welding experiment of thrust chamber plate

4 結(jié)論

1)建立的推力室外壁坡口焊接參數(shù)模型，可實(shí)現(xiàn)根據(jù)坡口精確計(jì)算匹配出焊接電流、送絲速度等關(guān)鍵工藝參數(shù)。使得焊接參數(shù)能夠適應(yīng)推力室產(chǎn)品對(duì)接坡口狀態(tài)的差異性，實(shí)現(xiàn)對(duì)不同坡口狀態(tài)的產(chǎn)品焊接參數(shù)根據(jù)坡口狀態(tài)進(jìn)行調(diào)整。

2)針對(duì)液體火箭發(fā)動(dòng)機(jī)推力室專用材料S06，1Cr21Ni5Ti等材料建立了熱絲TIG焊接參數(shù)計(jì)算模型?？梢愿鶕?jù)板厚、焊接速度得出合適的焊接參數(shù)以及焊接參數(shù)窗口。建立板厚、焊接速度、焊接電流與焊縫截面積的參數(shù)模型，可以預(yù)估計(jì)算出焊縫的尺寸。

3)在SIMULINK平臺(tái)下引入模糊控制建立了機(jī)器人焊接自適應(yīng)系統(tǒng)，針對(duì)響應(yīng)速度仿真優(yōu)化了調(diào)節(jié)因子，驗(yàn)證了模型及控制系統(tǒng)的有效性。 采用自適應(yīng)焊接系統(tǒng)進(jìn)行了推力室坡口的模擬件試驗(yàn)，焊縫熔合良好，并進(jìn)行X光檢測(cè)合格。