常天根 高月華 劉其鵬

摘要：為高效高精度地預(yù)測(cè)大尺寸結(jié)構(gòu)FSW 焊接變形，基于考慮夾具接觸約束的熱彈塑性仿真結(jié)果，在分析塑性應(yīng)變的區(qū)域分布特點(diǎn)的基礎(chǔ)上運(yùn)用固有應(yīng)變法進(jìn)行分區(qū)映射，發(fā)展了一種有效的焊接變形預(yù)測(cè)方法。針對(duì)小尺寸6061-T6鋁合金板材FSW 焊接，將接觸約束狀態(tài)和固定約束狀態(tài)下的焊后變形預(yù)測(cè)值與實(shí)驗(yàn)值進(jìn)行了對(duì)比分析，結(jié)果表明接觸約束狀態(tài)的預(yù)測(cè)結(jié)果更接近實(shí)際的非對(duì)稱變形。將該方法應(yīng)用于大尺寸6061-T6鋁合金板材焊接，并將所得結(jié)果與大尺寸板材接觸約束狀態(tài)下熱彈塑性法結(jié)果進(jìn)行對(duì)比，結(jié)果表明，該方法可有效地預(yù)測(cè)大尺寸板材攪拌摩擦焊接變形。

關(guān)鍵詞：攪拌摩擦焊;固有應(yīng)變;分區(qū)映射;接觸約束;焊接變形

中圖分類號(hào)： TG453+.9????? 文獻(xiàn)標(biāo)識(shí)碼： A文章編號(hào)：1001-2303（2022）02-0017-09

Prediction of FSW Wwelding Deformation by Inherent Strain Method Considering Fixture Contact Constraints

CHANG Tiangen1， GAO Yuehua1，2， LIU Qipeng2

1. School of Locomotive and Vehicle Engineering， Dalian Jiaotong University， Dalian 116028， China

2. School of Aeronautical Manufacturing Engineering， Nanchang Hangkong University， Nanchang 330063， China

Abstract： In this paper， an effective prediction method is proposed to predict the FSW weld deformation for the large-size structures using partition mapping of inherent strain under the local-global mapping frame. The contact state between fix‐ tures and weldment is considered by contact constraints in the thermoelastic-plastic simulation of local structure. The de‐ tailed partition scheme is determined by the distribution characteristics of residual plastic strain on the basis of the prediction results of local structure， and the partition scheme is implemented to the global structure for the partition mapping of inher‐ ent strain. For the welding of small-size 6061-T6 aluminum alloy sheet， the weld deformations are predicted under the con‐ tact and fixed constraint states by thermoelastic-plastic simulation， and these results are analyzed and compared with the ex‐ perimental data. The results show that the prediction results considering contact constraints are closer to the experimental data and can capture actual asymmetric deformation， thus it is necessary that the contact constraints between fixtures and weldment are modeled in the local structure. The proposed method is applied to the welding of large-size 6061-T6 aluminum alloy plates， and the prediction results are compared with those obtained by thermoelastic-plastic method considering the contact constraints between fixtures and weldment. The results show that the proposed method can effectively predict the welding FSW deformation of large-size plates.

Keywords： friction stir welding; inherent strain; partition mapping; contact constraint; welding deformation

引用格式：常天根，高月華，劉其鵬.考慮夾具接觸約束的FSW焊接變形的固有應(yīng)變法預(yù)測(cè)[J].電焊機(jī)，2022，52（2）：17-25.

Citation：CHANGTiangen， GAO Yuehua， LIU Qipeng. Prediction of FSW Wwelding Deformation by Inherent Strain Method Considering Fixture Contact Constraints[J]. Electric Welding Machine， 2022， 52（2）：17-25.

0? 前言

近年來，軌道交通行業(yè)大力發(fā)展，軌道車輛輕量化成為備受關(guān)注的問題。由于鋁合金材料具有質(zhì)量輕、強(qiáng)度高等優(yōu)點(diǎn)，被廣泛應(yīng)用于軌道交通制造業(yè)[1]。傳統(tǒng)熔焊在焊接鋁合金材料時(shí)會(huì)出現(xiàn)裂紋、氣孔等缺陷。攪拌摩擦焊接（FSW）是一種新型固相連接技術(shù)[2]，能夠避免熔焊易產(chǎn)生的缺陷。雖然FSW技術(shù)具有諸多優(yōu)點(diǎn)，但焊后殘余應(yīng)力、殘余變形等問題依然存在[3-4]。大量的試驗(yàn)不僅耗費(fèi)人力、物力，而且難以再現(xiàn)攪拌摩擦焊接過程。熱彈塑性法是目前應(yīng)用最為普遍且被國內(nèi)外學(xué)者廣泛認(rèn)可的攪拌摩擦焊接數(shù)值模擬預(yù)測(cè)方法，但該方法計(jì)算耗時(shí)長，難以預(yù)測(cè)大尺寸構(gòu)件的殘余狀態(tài)。

日本學(xué)者Ueda[5-8]在20世紀(jì)70年代提出了固有應(yīng)變的概念，并采用固有應(yīng)變法研究了T型和I型接頭的焊接變形。Deng等[9-10]利用熱彈塑性法分析了熱輸入、焊接方法、焊接順序及加強(qiáng)筋間距對(duì)焊接變形的影響。國內(nèi)學(xué)者李鴻[11]首次提出將有限元模型的每層網(wǎng)格看作是一個(gè)桿，把桿收縮看作是固有應(yīng)變的和，并采用了一套完整的公式，將殘余塑性應(yīng)變經(jīng)過公式計(jì)算處理為焊縫附近的收縮力和收縮力矩，并將其理論應(yīng)用于船體焊接變形的預(yù)測(cè)。周晶[12]等首次將固有應(yīng)變法應(yīng)用到鋁合金焊接中，將計(jì)算得到的收縮力和收縮力矩施加到殼單元上，探討了將固有應(yīng)變法應(yīng)用于鋁合金結(jié)構(gòu)焊接變形預(yù)測(cè)與控制的可行性。高博[13]等以散貨船的雙層底分段為例，采用等效載荷法對(duì)其焊接變形進(jìn)行預(yù)測(cè)，驗(yàn)證了固有應(yīng)變等效載荷法預(yù)測(cè)船體雙層底結(jié)構(gòu)焊接變形的可行性。汪建華[14]等利用固有應(yīng)變溫度載荷法分析了大型焊接船體結(jié)構(gòu)對(duì)接環(huán)縫的變形和轎車副車架的焊接變形。侯志剛[15]根據(jù)熱彈塑性法分析的結(jié)果，對(duì)端部焊縫中的固有應(yīng)變值進(jìn)行了修正，首次將固有應(yīng)變法應(yīng)用于點(diǎn)焊結(jié)構(gòu)的變形分析。梁偉[16]等采用熱彈塑性法和固有應(yīng)變法相結(jié)合的方式來計(jì)算復(fù)雜薄板結(jié)構(gòu)的總體變形，考察了外部拘束對(duì)板材焊接變形的影響。固有應(yīng)變法的發(fā)展解決了工程中預(yù)測(cè)焊接變形耗時(shí)長的問題。

本文針對(duì)6061-T6鋁合金板材的FSW焊接，基于局部-全局映射思想，以考慮夾具接觸約束的熱彈塑性有限元法仿真結(jié)果為依據(jù)，分析塑性應(yīng)變的區(qū)域分布特點(diǎn)，采用固有應(yīng)變溫度載荷法進(jìn)行分區(qū)映射，進(jìn)行了FSW焊接變形的快速仿真。在試驗(yàn)驗(yàn)證接觸約束狀態(tài)下所得結(jié)果的有效性基礎(chǔ)上，將本方法應(yīng)用于大尺寸板材焊接，并與大尺寸板材接觸約束狀態(tài)下熱彈塑性法所得結(jié)果進(jìn)行對(duì)比分析，驗(yàn)證了本方法的有效性。

1? 分區(qū)域映射的固有應(yīng)變法

固有應(yīng)變溫度載荷法是通過提取小尺寸板材（局部結(jié)構(gòu)）熱彈塑性有限元法所得到的焊縫及其附近位置產(chǎn)生的焊后殘余塑性應(yīng)變仿真結(jié)果，轉(zhuǎn)化為固有應(yīng)變并施加到大尺寸板材（整體結(jié)構(gòu)）上，進(jìn)而獲得焊接變形[14，17-18]。板材攪拌摩擦焊后沿厚度方向的應(yīng)變幾乎為零，故焊接變形主要由縱向固有應(yīng)變和橫向固有應(yīng)變引起，可不考慮垂向變形的影響。

在FSW數(shù)值模擬工作中板材與墊板、板材與夾具之間的約束方式分為兩種：一種是接觸約束方式，即在板材與墊板、板材與夾具之間通過設(shè)置接觸面間的摩擦力和接觸屬性進(jìn)行約束;另一種是固定約束方式，采用力學(xué)等效約束條件建立有限元模型，對(duì)夾具和墊板區(qū)域處的節(jié)點(diǎn)進(jìn)行剛性固定約束（即x、y、z三個(gè)方向）。一般來講，焊接起始段、焊接中段以及焊接結(jié)束段的塑性應(yīng)變會(huì)存在差別。為了提高預(yù)測(cè)精度，本文進(jìn)行小尺寸板材熱彈塑性分析時(shí)考慮實(shí)際夾具及墊板的接觸約束，以更有效地捕獲這種差異。將殘余塑性應(yīng)變分布沿焊接方向分為3段，如圖1所示。

1.1? 固有應(yīng)變的分區(qū)域提取和計(jì)算

縱向收縮可近似理解為由焊縫及其附近區(qū)域的縱向收縮力造成的。對(duì)于分區(qū)后每一區(qū)域，引用 Tendon Force概念[19]來表示縱向收縮狀況，其公式定義為[20]：

式中E為常溫時(shí)材料的彈性模量;εxp為穩(wěn)定焊接階段橫截面的縱向塑性應(yīng)變;A 為橫截面上產(chǎn)生塑性應(yīng)變的面積。

由此，縱向固有應(yīng)變可表示為：

式中B 為彈性計(jì)算時(shí)施加固有應(yīng)變區(qū)域的寬度; H為施加固有應(yīng)變區(qū)域的厚度。

傳統(tǒng)的橫向固有應(yīng)變提取方法是將橫向塑性應(yīng)變對(duì)其面積進(jìn)行積分，計(jì)算過程復(fù)雜。本文采用遲哲[18]的計(jì)算方法，即在小尺寸板材結(jié)構(gòu)中性面上選擇焊縫兩側(cè)包含焊縫區(qū)和熱影響區(qū)的兩條路徑，通過獲取兩條路徑之間距離的變化量S來計(jì)算橫向收縮量。由此，橫向固有應(yīng)變可表示為：

考慮FTendon所產(chǎn)生的橫向收縮，得到最終橫向固有應(yīng)變?chǔ)舩*計(jì)算公式為：

式中ν為泊松比。

1.2? 固有應(yīng)變的分區(qū)域施加

關(guān)于固有應(yīng)變?cè)诖蟪叽绨宀纳希ㄕw結(jié)構(gòu)）施加區(qū)域的確定方法，本文采用Mun和Seo[21]的方法，將Al6061-T6屈服應(yīng)力快速變化的區(qū)域指定為固有應(yīng)變的施加區(qū)域。圖2給出了Al6061-T6屈服應(yīng)力隨溫度變化曲線，由此確定溫度在380℃及其以上的區(qū)域?yàn)楣逃袘?yīng)變的施加區(qū)域。

由于在大尺寸板材上難以直接施加固有應(yīng)變，通過將其轉(zhuǎn)換為材料的各向異性熱膨脹系數(shù)進(jìn)行等效加載，等效公式為：

式中 α為熱膨脹系數(shù);T為單位溫度載荷。針對(duì)每一區(qū)域，應(yīng)用式（5）計(jì)算出等效的熱膨脹系數(shù)，并進(jìn)行對(duì)應(yīng)區(qū)域施加，經(jīng)過一次彈性計(jì)算即可預(yù)測(cè)焊接變形。

2? 考慮夾具接觸約束的小尺寸板材焊接變形的熱彈塑性法預(yù)測(cè)

2.1? 考慮夾具接觸約束的小尺寸板材焊接的有限元模型

采用材料為6061-T6鋁合金、尺寸為620 mm× 200 mm×5 mm的兩塊小尺寸板材FSW進(jìn)行焊接變形仿真。攪拌針直徑7 mm，攪拌針長度4.9 mm，軸肩直徑22 mm，軸肩壓力8 000 N，焊速60 mm/min，攪拌頭轉(zhuǎn)速400 r/min 。焊接從焊縫前端25 mm處開始，預(yù)熱15 s，至焊縫末端25 mm處終止。將小尺寸模型中焊縫區(qū)域的網(wǎng)格細(xì)化，最小網(wǎng)格尺寸2 mm× 2 mm×1.67 mm，得到的有限元模型如圖3所示。有限元模擬過程中將接觸屬性設(shè)置分為切向和法向兩方面，切向接觸屬性設(shè)置為庫倫摩擦，取摩擦系數(shù) 0.3。法向的接觸屬性設(shè)置為硬接觸，采用罰方法，其中硬接觸表示當(dāng)兩個(gè)接觸面發(fā)生接觸，接觸面間傳遞的法向壓力其大小不受限制;當(dāng)接觸面間的接觸壓力變?yōu)榱慊蜇?fù)值時(shí)，表示接觸面分離，軟件根據(jù)接觸材料屬性自動(dòng)計(jì)算罰剛度。在硬接觸中，法向采用罰方法，ABAQUS軟件根據(jù)接觸材料屬性自動(dòng)計(jì)算罰剛度，既能有效避免數(shù)值奇異問題，也能充分考慮接觸材料的屬性，適合攪拌摩擦焊接工況。

考慮環(huán)境溫度為20℃，薄板與空氣接觸處對(duì)流換熱系數(shù)為30 W/（m2·℃）。薄板上、下分別由夾具和墊板固定，其對(duì)流換熱系數(shù)為1 000 W/（m2·℃）。夾具放置在距離焊縫兩側(cè)各40 mm處，寬度30 mm。6061-T6鋁合金的熱物理性能和熱力學(xué)性能參數(shù)如表 1所示[22]。

在熱彈塑性仿真數(shù)值模擬過程中，將攪拌頭產(chǎn)熱以熱源的形式進(jìn)行施加，包括軸肩和焊件的接觸區(qū)域以及攪拌針的作用區(qū)域。利用攪拌頭的扭矩來計(jì)算焊接時(shí)總的熱量輸入，總熱量分配到軸肩和焊件接觸區(qū)域的比例為75%，分配到攪拌針作用區(qū)域的比例為25%。熱源功率計(jì)算公式可表示為[23]：

式中Qtotal為焊接總輸入功率，Qs 0.75Qtotal，Qp0.25Qtotal;η為焊接過程的熱利用率，為95%;ω為攪拌頭旋轉(zhuǎn)角速度;n為攪拌頭旋轉(zhuǎn)速度;R0為攪拌頭軸肩半徑;R 1為攪拌針半徑;H為攪拌針高度;Mz為攪拌頭扭矩，對(duì)于6061-T6鋁合金材料的FSW其攪拌頭扭矩可用式（9）進(jìn)行計(jì)算[22]。此外，模型中考慮了攪拌頭下壓力的影響，其計(jì)算公式為[24]：

式中FZ 為攪拌頭的下壓力;P為添加到模型上的壓力。

2.2? 小尺寸板材FSW焊接變形的熱彈塑性法預(yù)測(cè)

基于上述有限元數(shù)值仿真分析方案進(jìn)行計(jì)算，得到考慮夾具接觸約束的小尺寸板材6061-T6鋁合金FSW焊后變形情況，如圖4a所示。垂向變形圖呈中間下凹、兩端上凸的反馬鞍形狀，垂向相對(duì)位移量為 10.02 mm 。為進(jìn)行有效對(duì)比，也對(duì)固定約束（墊板和夾具接觸處為位移固定約束）進(jìn)行熱彈塑性仿真，如圖4b所示，垂向相對(duì)位移量為8.67 mm 。針對(duì)本節(jié)研究的6061-T6鋁合金板材FSW焊接變形，文獻(xiàn)[25]給出了試驗(yàn)測(cè)量的變形結(jié)果。提取圖5中的 A-A、B-B、C-C、D-D 四個(gè)方向變形量進(jìn)行對(duì)比分析（圖中試驗(yàn)值引自文獻(xiàn)[25]），如圖6所示。

由圖6可知，在A-A和B-B兩個(gè)方向上，接觸約束狀態(tài)下曲線峰值位置較固定約束狀態(tài)下曲線峰值位置左移，預(yù)測(cè)出了A-A和B-B兩個(gè)方向上垂向位移呈不對(duì)稱的特征，且綜合對(duì)比兩種約束方式在各方向上的相對(duì)變形量時(shí)發(fā)現(xiàn)：接觸約束狀態(tài)的相對(duì)變形量更接近實(shí)驗(yàn)值。而固定約束情況下的變形結(jié)果在焊接方向上基本呈對(duì)稱趨勢(shì)，與實(shí)際存在偏差。對(duì)比實(shí)驗(yàn)數(shù)據(jù)和兩種約束方式的相對(duì)變形量，帶墊板和夾具模型的接觸約束狀態(tài)下的焊接變形情況更符合實(shí)際變形，因此基于該方法所得的應(yīng)變結(jié)果進(jìn)行大尺寸板材的焊接變形預(yù)測(cè)將更加準(zhǔn)確。

3? 考慮夾具接觸約束的大尺寸板材焊接變形的固有應(yīng)變法預(yù)測(cè)

3.1? 分區(qū)固有應(yīng)變的提取和計(jì)算

基于熱彈塑性方法得到的小尺寸板材縱向殘余塑性應(yīng)變?cè)茍D如圖7所示。可以看出，縱向塑性應(yīng)變分布不均勻。為了提高預(yù)測(cè)精度，沿焊接方向?qū)⑵浞譃?個(gè)區(qū)域（見圖7），分別提取4個(gè)區(qū)域?qū)?yīng)單元的縱向塑性應(yīng)變值及其面積，求出縱向塑性應(yīng)變平均值。依據(jù)1.1節(jié)方法進(jìn)行固有應(yīng)變的計(jì)算和施加。

圖8為小尺寸板材穩(wěn)態(tài)焊接過程中的溫度分布云圖。由于焊速較低，受溫度影響區(qū)域呈現(xiàn)出不明顯的橢圓形，但可以看出攪拌頭前方溫度梯度較大，后方溫度梯度較小。在1.2節(jié)中，確定溫度在380℃及其以上的區(qū)域?yàn)楣逃袘?yīng)變的施加區(qū)域，其溫度分布云圖如圖9所示。由圖9可知，在380℃及其以上時(shí)，板材上、下表面的溫度分布大致呈圓形。由此可確定固有應(yīng)變施加區(qū)域的斷面近似為倒梯形，如圖10所示。依據(jù)1.2節(jié)方法，各區(qū)域計(jì)算所得的熱膨脹系數(shù)值如表2所示。

3.2? 考慮夾具接觸約束的大尺寸板材焊接變形的固有應(yīng)變法預(yù)測(cè)

針對(duì)兩塊尺寸為1 240 mm×200 mm×5 mm 的6061-T6鋁合金板材進(jìn)行FSW對(duì)接焊的焊接變形預(yù)測(cè)試驗(yàn)。將該模型依據(jù)圖7進(jìn)行同比例的區(qū)域劃分，將2.2節(jié)計(jì)算得到的等效材料膨脹系數(shù)分區(qū)域施加，并進(jìn)行彈性計(jì)算，預(yù)測(cè)結(jié)果如圖11a所示，最大變形量約24.7 mm 。為驗(yàn)證結(jié)果的有效性，對(duì)該焊接結(jié)構(gòu)進(jìn)行考慮夾具接觸約束的熱彈塑性分析，變形結(jié)果如圖11b所示，最大變形量約為24.6 mm。

按圖5所示提取A-A、B-B、C-C、D-D 四個(gè)方向變形量進(jìn)行對(duì)比分析，如圖12所示。本文方法與熱彈塑性法的變形曲線趨勢(shì)基本一致。兩種方法所得的相對(duì)變形量的數(shù)據(jù)如表3所示，可以看出，文中方法有較高的預(yù)測(cè)精度。

4? 結(jié)論

本文針對(duì)小尺寸6061-T6鋁合金板材FSW 焊接，探討了不同約束方式對(duì)焊接變形預(yù)測(cè)的影響，結(jié)果表明考慮夾具、墊板與焊板之間的接觸約束所得的變形沿焊接方向是不對(duì)稱的，與實(shí)際的焊接變形相吻合且整體預(yù)測(cè)值精度較高。該結(jié)果說明在進(jìn)行大尺寸結(jié)構(gòu)FSW焊后變形的快速預(yù)測(cè)時(shí)，局部模型中考慮夾具、墊板與焊板之間的接觸約束是必要的。

依據(jù)接觸約束狀態(tài)下局部結(jié)構(gòu)殘余塑性應(yīng)變的區(qū)域分布特點(diǎn)，采用固有應(yīng)變法對(duì)整體結(jié)構(gòu)進(jìn)行分區(qū)映射，發(fā)展出一種有效的焊接變形預(yù)測(cè)方法。將此方法應(yīng)用于大尺寸板材FSW焊接變形預(yù)測(cè)，其結(jié)果與熱彈塑方法所得結(jié)果較為一致，具有較高的預(yù)測(cè)精度。

針對(duì)復(fù)雜結(jié)構(gòu)，F(xiàn)SW焊接變形的快速有效預(yù)測(cè)仍然需要進(jìn)一步研究，后續(xù)工作將在本研究的基礎(chǔ)上進(jìn)行拓展和改進(jìn)，嘗試應(yīng)用于復(fù)雜結(jié)構(gòu)的FSW焊接變形預(yù)測(cè)。

參考文獻(xiàn)：

[1] 趙衍華.鋁合金的攪拌摩擦焊接[M].北京：中國宇航出版社，2009.

ZHAO Yanhua. Friction Stir Welding of Aluminum Al‐ loy[M]. Beijing：China Aerospace Press，2009.

[2]孫宜華，杜良.攪拌摩擦焊的研究進(jìn)展與應(yīng)用[J].新技術(shù)新工藝，2011，（6）：70-73.

SUN Yihua，DU Liang. Development and Application of Friction Stir Welding[J]. New Technology & New Process，2011（6）：70-73.

[3]王春蕊.海洋平臺(tái)生活模塊鋁合金材料的攪拌摩擦焊接技術(shù)研究[J].化工管理，2015，（18）：95-95.

WANG Chunrui. Friction? Stir Welding Technology of Aluminum Alloy for Offshore Platform Living Module [J].Chemical Industry Management，2015（18）：95-95.

[4]張?zhí)飩}，郭德倫，欒國紅，等.固相連接新技術(shù)-攪拌摩擦焊技術(shù)[J].航空制造技術(shù)，1999（2）：35-36.ZHANG Tiancang，GUO Delun，LUAN Guohong，et al. A New Technology for Solid Phase Bonding Friction Stir Welding[J]. Aeronautical Manufacturing Technol‐ ogy，1999（2）：35-36.

[5] Ueda Y，Nalacho K. A new measuring method of re‐sidual stress with the aid of finiteelement method of reliability of estinated values[J]. Journal of the Society of Naval Architects of Japan，1971，2（2）：499-507.

[6] Ueda Y，Nalacho K，Tanigawa W，et al. New measur‐ing? method? of three? in? dimensional? stress? based? on theory of inherent strain[J]. Tansactions of the JWRI （Japan Welding Research? Institute）（Japan），1979，8（2）：89-96.

[7] Ueda Y，Tanigawa W. Prediction of residual stresses inbutt welded plates using inherent strains[J]. Journal of? Materials? Processing? Technobgy，1993，115（10）：417-423.

[8] Yuan M G，Ueda Y. Prediction of residual stresses inwelded T- and I joints using inherent strains[J]. Journal? of Materials Processing Technology ，1996，118（4）：229-234.

[9] Dean Deng，WeiLiang. Numerical simulation of weld‐ing distortion in large structures[J]. Computer Methods in Applied Mechanics and Engineering，2007，196（45-48）：4613-4627.

[10] Dean? Deng ，WeiLiang. FEM? prediction? of bucklingdistortion induced by welding in the? plate panel struc‐tures[J]. Computational Materials Science，2008，43（4）：591-607.

[11]李鴻.基于固有應(yīng)變的船體分段焊接變形預(yù)測(cè)[D].黑龍江：哈爾濱工程大學(xué)，2005.

LI Hong. Welding Deformation Prediction of Hull Sec‐ tions based on Inherent Strain[D]. Harbin：Harbin En‐ gineering University，2005.

[12]周晶，常保華，張驊，等.采用固有應(yīng)變法預(yù)測(cè)鋁合金焊接變形[J].焊接技術(shù)，2010，39（06）：6-11.

ZHOU Jing，CHANG Baohua，ZHANG Hua，et al. A New? Method? for? Predicting? Welding? Deformation? of Aluminum Alloy[J]. Welding Technology，2010，39（06）：6-11.

[13]高博，陳章蘭，李昆鵬，等.基于等效載荷法的復(fù)雜結(jié)構(gòu)焊接變形預(yù)測(cè)[J].焊接技術(shù)，2013，42（9）：64-67.?? GAO Bo，CHEN Zhanglan，LI Kunpeng，et al. Weld‐ ing? Deformation? Prediction? of? ComplexStructuresbased on Equivalent Load Method[J]. Welding Tech‐ nology，2013，42（9）：64-67.

[14]汪建華，陸?zhàn)?，魏良?固有應(yīng)變有限元法預(yù)測(cè)焊接變形理論及其應(yīng)用[J].焊接學(xué)報(bào)，2002，23（6）：36-40. WANG Jianhua，LU Hao，WEI Liangwu. The Predic‐ tion of Welding Deformation by Inherent Strain Finite Element Method and its Application[J]. Transactions of the China Welding Institution，2002，23（6）：36-40.

[15]侯志剛.薄板結(jié)構(gòu)焊接變形的預(yù)測(cè)與控制[D].湖北：華中科技大學(xué)，2005.

HOU Zhigang. Prediction and Control of Welding De‐ formation? of Thin? Plate? Structure[D]. Hubei：Hua‐ zhong University of Science and Technology，2005.

[16]梁偉，鄭穎，鄧德安，等.外拘束對(duì)鋁合金薄板結(jié)構(gòu)焊接變形的影響[J].機(jī)械工程學(xué)報(bào)，2021，57（6）：71-72. LIANG Wei，ZHENG Ying，DENG Dean，et al. Ef‐ fect of External Constraint on Welding Deformation of Aluminum Alloy? Sheet[J]. Chinese? Journal? of Me‐ chanical Engineering，2021，57（6）：71-72.

[17]徐濟(jì)近，陳立功.基于固有應(yīng)變法筒體對(duì)接多道焊焊接變形的預(yù)測(cè)[J].焊接學(xué)報(bào)，2007，28（1）：77-80.?? ?XU Jijing，CHEN Ligong. Prediction of Welding De‐ formation? in? Multi-pass? Welding? of? Cylinder? Butt based on Natural Strain Method[J]. Transactions of the China Welding Institution，2007，28（1）：77-80.

[18]遲哲，楊鑫華.基于固有應(yīng)變法的地鐵側(cè)墻FSW焊接變形仿真[J].電焊機(jī)，2017，10（47）：1-7.

CHI Zhe，YANG Xinhua. Deformation Simulation of FSW Welding of Subway Side Wall Based on Inherent Strain Method[J]. Electric Welding Machine，2017，10（47）：1-7.

[19]羅宇，魯華益，謝雷，等. Tendon Force的概念及計(jì)算方法[J].造船技術(shù)，2004（04）：38-40.

LUO Yu，LU Huayi，XIE Lei，et al. An Overview of Tendon? Force? and? its? Applications[J]. Shipbuilding Technology，2004（04）：38-40.

[20] Deng D，Murakawa H，Ling W. Numerical Simulationor Welding? Distortion? in? Large? Structures[J]. Com‐ puter? Methods? In Applied? Mechanics And? Engineer‐ ing，2007，196（45-48）：4613-4627.

[21] Hyung-Suk Mun，Sung-Il Seo. Welding strain analysisof friction stir-welded aluminum? alloy structures using? inherent? strain-based? equivalent? loads[J]. Journal? of? Mechanical? Science and Technology，2013，27（9）：2775-2782.

[22]劉玉龍.軌道車輛板材攪拌摩擦焊接殘余狀態(tài)仿真研究[D].遼寧：大連交通大學(xué)，2018.

LIU Yulong. Simulation of Residual? State of Friction Stir Welding for Rail Vehicle Plate[D]. Liaoning：Da‐ lianJiaotong University，2018.

[23]Vuyst T，Dealvise L D. Inverse Analysis using a Ge‐netic Algorithm? for? the? Finite? Element? Modeling? of Friction Stir Welding[C]. Proc.5th Int. SympOnFric‐ tion Stir Welding. USA：ParkCity，2004.

[24]張正偉.攪拌摩擦焊接構(gòu)件殘余狀態(tài)和疲勞壽命研究[D].遼寧：大連理工大學(xué)，2014.

ZHANG Zhengwei. Residual State and Fatigue Life of Friction Stir Welded Components[D]. Liaoning：Da‐ lian University of Technology，2014.

[25]付強(qiáng).鋁合金薄板攪拌摩擦焊殘余應(yīng)力及失穩(wěn)變形的預(yù)測(cè)研究[D].天津：天津大學(xué)，2010.

FU Qiang，Effect of Friction Stir Welding on the Re‐ sidual Stress and Instability Deformation of Aluminum Alloy[D]. Tianjin：Tianjin University，2010.

編輯部網(wǎng)址：http：//www.71dhj.com