基于后處理再造型的連續(xù)制造工藝過(guò)程仿真技術(shù)

胡廣旭　孟梅　劉冰

摘要：制造工藝仿真技術(shù)目前僅應(yīng)用于單點(diǎn)工藝分析，缺乏在連續(xù)工藝實(shí)施中考慮應(yīng)力和變形等變化對(duì)前后工序的影響，不能實(shí)現(xiàn)連續(xù)的制造工藝仿真，因此提出一種基于后處理再造型（Post-Processing Re-Modeling，PPRM）的仿真技術(shù).該技術(shù)利用前道工序仿真結(jié)果的后處理數(shù)據(jù)重構(gòu)后道工序的仿真模型，從而可以系統(tǒng)、連貫地考慮各制造階段的工藝變形、應(yīng)力和應(yīng)變等，實(shí)現(xiàn)連續(xù)制造工藝過(guò)程仿真.將這項(xiàng)技術(shù)應(yīng)用于船體結(jié)構(gòu)部件裝焊變形的仿真計(jì)算中，計(jì)算多個(gè)零件連續(xù)裝焊工藝過(guò)程的焊接變形積累，從而預(yù)測(cè)并監(jiān)控多道裝焊工藝前后工序之間的互相影響干涉.該技術(shù)在船體制造中的肋板彎曲工藝、分段堆放和分段吊裝等工藝中應(yīng)用前景較好.

關(guān)鍵詞：制造工藝仿真； 熱彈塑性有限元； 船體結(jié)構(gòu)； 焊接變形； 船舶制造； 分段吊裝

中圖分類號(hào)： TH164；TG404

文獻(xiàn)標(biāo)志碼：B

Abstract：The current manufacturing process simulation technology is only applied in single point process analysis without considering the effect of the change of stress and deformation on former and later processes in the implementation of continuous process， so the continuous manufacturing process can not be implemented. Hence， a simulation technology is proposed on the basis of Post-Processing Re-Modeling（PPRM）. With the technology， the post-processing data of former process simulation is used to re-model the subsequent process， and the process deformation， stress and strain in each manufacturing stage can be systematically and continuously considered and the contiuous manufacturing process simulation can be implemented. The technology is applied in the simulation and calculation of welding deformation of components of ship hull structure， and the welding deformation accumulation are calculated in the continuous welding process of multi-components， and the mutual interference before and after the process of multi-welding steps is predicted and monitored. The application prospect of the technology is good for ship hull manufacturing processes such as rib bending， block laying， block lifting， and so on.

Key words：manufacture process simulation； thermal-elastic-plastic finite element； ship hull structure； welding deformation； shipbuilding； block lifting

0引言

在大型工業(yè)產(chǎn)品，如飛機(jī)及船舶與海洋工程裝備等大型結(jié)構(gòu)的制造過(guò)程中，從材料切割到最終成品，往往歷經(jīng)多道制造工藝.在每道工藝操作中，設(shè)備都會(huì)因工藝設(shè)備操作時(shí)產(chǎn)生的物理或化學(xué)作用而引起變形、內(nèi)應(yīng)力、組織轉(zhuǎn)變和裂紋等，影響產(chǎn)品制造精度和疲勞壽命.制造質(zhì)量是決定產(chǎn)品能否滿足用戶需求的關(guān)鍵.不論多么優(yōu)秀的產(chǎn)品設(shè)計(jì)，如果缺少優(yōu)質(zhì)的制造工藝，都難以成為優(yōu)質(zhì)的產(chǎn)品.

針對(duì)制造工藝質(zhì)量問(wèn)題，很多研究均采用制造工藝仿真技術(shù)，在工藝設(shè)計(jì)之初進(jìn)行計(jì)算、預(yù)測(cè)和分析，一方面提前預(yù)判工藝實(shí)施方案、優(yōu)化工藝參數(shù)，另一方面可以提前準(zhǔn)備修補(bǔ)方案.目前，制造工藝仿真技術(shù)在很多領(lǐng)域得到廣泛應(yīng)用，包括熱處理、鑄造、焊接、沖壓和冷彎加工等.[1-5]

工藝仿真技術(shù)主要應(yīng)用于單點(diǎn)工藝仿真計(jì)算，僅有少數(shù)研究在焊接工藝仿真后進(jìn)行焊后熱處理仿真，且這些仿真技術(shù)僅限于幾何模型不變的情況.[6-8]因此，目前仍缺乏產(chǎn)品制造全過(guò)程的仿真研究，如在船體外板制造過(guò)程中，冷彎加工后的鋼材內(nèi)部會(huì)產(chǎn)生殘余變形和應(yīng)力，這些材料變化對(duì)后序水火彎板工藝或焊接裝配工藝的影響很難用單點(diǎn)工藝仿真進(jìn)行預(yù)測(cè)分析.由于缺少多道連續(xù)工序的仿真技術(shù)手段，導(dǎo)致前道工序的變形、應(yīng)力和組織轉(zhuǎn)變等缺陷很難在后道工序中考慮到，使得工藝仿真計(jì)算脫離實(shí)際制造過(guò)程.

針對(duì)上述問(wèn)題，提出一種基于后處理再造型（Post-Processing Re-Modeling， PPRM）的仿真技術(shù)，利用該技術(shù)可實(shí)現(xiàn)連續(xù)工藝仿真，從而可以考慮連續(xù)工藝各階段缺陷的積累.

1PPRM仿真技術(shù)

在多道連續(xù)工序的仿真計(jì)算中，利用前一道工序的仿真后處理數(shù)據(jù)重構(gòu)后一道工序的仿真模型，從而在后道工序仿真計(jì)算中考慮前道工序的工藝影響，以達(dá)到全過(guò)程的仿真計(jì)算目的，即為PPRM技術(shù)，其具體實(shí)現(xiàn)流程見(jiàn)圖1.首先在前道工序建模過(guò)程中分別導(dǎo)出幾何或有限元網(wǎng)格的標(biāo)準(zhǔn)模板數(shù)據(jù)，之后定義前道工序初始條件和邊界條件，進(jìn)行仿真計(jì)算.在計(jì)算完畢后，采用后處理開(kāi)發(fā)技術(shù)輸出仿真后處理結(jié)果數(shù)據(jù)，并依據(jù)標(biāo)準(zhǔn)模板幾何或網(wǎng)格進(jìn)行數(shù)據(jù)轉(zhuǎn)換，進(jìn)一步與前道工序標(biāo)準(zhǔn)模板數(shù)據(jù)進(jìn)行數(shù)據(jù)疊加和重構(gòu)，實(shí)現(xiàn)幾何或網(wǎng)格再造型，獲得再造型模型.以此為基礎(chǔ)，在后道工序中加載再造型模型的仿真邊界條件，同時(shí)導(dǎo)入前道工序的物理場(chǎng)后處理數(shù)據(jù)結(jié)果，實(shí)現(xiàn)PPRM技術(shù)全過(guò)程.反復(fù)按照PPRM技術(shù)流程建立連續(xù)工藝過(guò)程仿真模型，最終實(shí)現(xiàn)全過(guò)程連續(xù)工藝仿真.

實(shí)現(xiàn)PPRM技術(shù)有以下要點(diǎn).

1）PPRM規(guī)劃.針對(duì)不同的制造工藝仿真問(wèn)題，PPRM關(guān)注的重點(diǎn)不同，所需的后處理數(shù)據(jù)也不同，例如：研究制造精度時(shí)往往著重于變形數(shù)據(jù)再造，研究產(chǎn)品制造對(duì)強(qiáng)度的影響時(shí)則更關(guān)注制造應(yīng)力和應(yīng)變甚至裂紋等缺陷，研究熱加工情況下的金屬性能變化時(shí)則需模擬加工過(guò)程中組織的變化.在確定研究關(guān)鍵點(diǎn)后即可規(guī)劃PPRM技術(shù)提取的后處理數(shù)據(jù)，確定模型重構(gòu)所需算法.

2）標(biāo)準(zhǔn)模板幾何或網(wǎng)格的數(shù)據(jù)操作技術(shù).目前造型和仿真軟件大多采用國(guó)外商用軟件，其模型數(shù)據(jù)格式各異.為實(shí)現(xiàn)PPRM技術(shù)，首先將各種數(shù)據(jù)格式轉(zhuǎn)換為標(biāo)準(zhǔn)模板數(shù)據(jù)格式（如iges，stp和inp等），然后基于標(biāo)準(zhǔn)模板數(shù)據(jù)格式開(kāi)發(fā)數(shù)據(jù)操作模塊.針對(duì)不同的模板數(shù)據(jù)類型，采用的開(kāi)發(fā)函數(shù)不同，可根據(jù)具體需要進(jìn)行選取.

3）后處理結(jié)果數(shù)據(jù)自動(dòng)輸出技術(shù).在有限元仿真計(jì)算中，輸出結(jié)果數(shù)據(jù)以單元或節(jié)點(diǎn)為基本數(shù)據(jù)存儲(chǔ)單元，由于仿真計(jì)算網(wǎng)格數(shù)據(jù)量巨大，因此難以手動(dòng)輸出.目前，國(guó)外商用仿真計(jì)算軟件一般提供相應(yīng)后處理二次開(kāi)發(fā)語(yǔ)言和開(kāi)發(fā)函數(shù)庫(kù)，可按照需求靈活地開(kāi)發(fā)自動(dòng)數(shù)據(jù)輸出功能.較為常用的后處理語(yǔ)言有Python和C等.

4）后處理數(shù)據(jù)與再造型模型數(shù)據(jù)交互技術(shù).實(shí)施PPRM的模型再造過(guò)程是利用仿真結(jié)果后處理數(shù)據(jù)對(duì)初始模型數(shù)據(jù)進(jìn)行疊加、修改和重構(gòu)的交互過(guò)程，涉及有限元節(jié)點(diǎn)的幾何、力學(xué)和物理場(chǎng)等信息加載，具體包括后處理數(shù)據(jù)格式轉(zhuǎn)換、標(biāo)準(zhǔn)模板數(shù)據(jù)格式識(shí)別以及后處理數(shù)據(jù)添加和替換等操作程序，可依據(jù)仿真后處理數(shù)據(jù)類型和存儲(chǔ)格式開(kāi)發(fā)對(duì)應(yīng)的數(shù)據(jù)交互功能.

2PPRM技術(shù)仿真實(shí)例

由于PPRM技術(shù)的應(yīng)用不局限于某種特定工藝的仿真計(jì)算，本文將該技術(shù)應(yīng)用于船體某部件上多個(gè)肋板的連續(xù)裝焊仿真中，通過(guò)仿真觀察前道肋板裝焊后產(chǎn)生的變形對(duì)后道肋板裝焊的影響.

2.1多肋板船體部件焊接變形的PPRM建模

焊接變形是影響船體結(jié)構(gòu)制造精度的關(guān)鍵因素之一.采用PPRM技術(shù)預(yù)測(cè)船體某部件肋板裝焊變形過(guò)程，明確每條肋板焊接對(duì)結(jié)構(gòu)的影響，揭示焊接變形積累過(guò)程，輔助研究設(shè)計(jì)精度控制對(duì)策.

某船體外板結(jié)構(gòu)模型見(jiàn)圖2.圖2中虛線為部件與曲面外板連接的裝配線.該裝配線制造精度直接影響其后續(xù)與外板的裝配，因此通過(guò)仿真揭示肋板裝焊過(guò)程引起的裝配線上的變形.

在部件現(xiàn)場(chǎng)裝焊過(guò)程中，為避免操作空間干涉，肋板按照一定順序裝焊（見(jiàn)圖3），依次為Rib1，Rib2，Rib3和Rib4，具體仿真計(jì)算時(shí)在每條肋板裝焊變形仿真后均進(jìn)行PPRM.圖3位移提取點(diǎn)D1， D2， D3， D4和D5均為裝配線上的節(jié)點(diǎn).由于本例主要關(guān)注裝焊精度，因此在PPRM技術(shù)應(yīng)用中輸出標(biāo)準(zhǔn)模板三維幾何或網(wǎng)格數(shù)據(jù)，此處選取inp網(wǎng)格模板文件.由于肋板間距較遠(yuǎn)，可忽略應(yīng)力的相互影響，后處理主要提取空間位移作為模型重構(gòu)數(shù)據(jù)，與初始模型的inp網(wǎng)格模板文件進(jìn)行數(shù)據(jù)交互，重構(gòu)其數(shù)據(jù)實(shí)現(xiàn)再造型，然后在再造型模型上建立預(yù)裝焊的肋板模型，如此反復(fù)循環(huán)，最終完成所有肋板裝焊仿真計(jì)算.具體流程見(jiàn)圖4.

先進(jìn)行初始肋板裝焊的幾何和網(wǎng)格建模，即Ribn，輸出標(biāo)準(zhǔn)模板網(wǎng)格數(shù)據(jù)，同時(shí)加載初始條件和邊界條件，完成Ribn的焊接變形仿真計(jì)算.進(jìn)一步提取Ribn焊接仿真產(chǎn)生的變形位移信息，輸出并與Ribn標(biāo)準(zhǔn)模板進(jìn)行數(shù)據(jù)交互，從而將Ribn產(chǎn)生的焊接變形數(shù)據(jù)重構(gòu)至后道工序Ribn+1網(wǎng)格模型中，使Ribn+1仿真模型考慮Ribn的殘余焊接變形.如此循環(huán)就可實(shí)現(xiàn)連續(xù)工序焊接過(guò)程的集成仿真.

2.2三維焊接變形有限元模型

變形仿真主要采用熱彈塑性有限元理論[9]，船體部件三維有限元模型見(jiàn)圖5，采用全六面體網(wǎng)格，且應(yīng)用過(guò)渡網(wǎng)格技術(shù)控制網(wǎng)格數(shù)量.焊接熱源模型采用雙橢球熱源模型模擬CO2氣體保護(hù)焊[10]，材料為C15鋼，材料參數(shù)隨溫度變化見(jiàn)文獻(xiàn)[11]，熱交換采用空氣對(duì)流邊界條件.

2.3全裝焊過(guò)程仿真結(jié)果

有效顯示為焊接溫度場(chǎng)結(jié)果，采用后處理技術(shù)將超過(guò)2 000 ℃的部分統(tǒng)一歸為熔化金屬.肋板Rib1和Rib4焊接瞬態(tài)溫度場(chǎng)仿真結(jié)果見(jiàn)圖6.

再造型前后網(wǎng)格模型對(duì)比見(jiàn)圖7，可清晰地觀察到由焊接工藝引起的結(jié)構(gòu)變形.

由圖8可知，每條肋板焊后都引起肋板附近區(qū)域的變形，尤其是在裝配線上.一方面裝配線邊緣缺少約束，自由變形較大；另一方面由于肋板的焊縫端部距離邊緣裝配線較近，導(dǎo)致變形較大.肋板Rib1，

Rib2，Rib3和Rib4焊后的再造型模型與初始模型在相同位置節(jié)點(diǎn)位移差值見(jiàn)表1.通過(guò)相同節(jié)點(diǎn)在不同肋板裝焊階段的位移差值可以觀察全過(guò)程變形的積累.Rib4焊后差值是所有肋板焊接變形的最終積累結(jié)果，即為最終變形.由圖3可知，Rib1，Rib2和Rib3方向趨近于與裝配線垂直，其肋板端部與邊緣接近，且焊縫較長(zhǎng)，導(dǎo)致焊后變形較大.由表1可知，Rib1，Rib2和Rib3焊后分別在焊縫較近的D2，D4和D3節(jié)點(diǎn)位移較大，而Rib4焊縫由于距離裝配線邊緣較遠(yuǎn)且焊縫較短，其焊后變形較小，幾乎不影響節(jié)點(diǎn)位移.

基于PPRM技術(shù)實(shí)現(xiàn)焊接變形仿真，考慮前道裝焊變形對(duì)后道裝焊的影響，使制造工藝仿真更接近現(xiàn)場(chǎng)工藝實(shí)施操作，有助于焊接變形控制方案的制定.

3PPRM在船舶制造領(lǐng)域應(yīng)用探討

船舶制造過(guò)程包括多個(gè)工序環(huán)節(jié)，例如在將鋼板制成船體曲面外板過(guò)程中，需經(jīng)過(guò)預(yù)處理、等離子切割、初步冷彎成型、水火彎板成型和肋板裝焊等過(guò)程.在各工序?qū)嵤┻^(guò)程中，不可避免地會(huì)產(chǎn)生熱變形、應(yīng)力集中、精度偏差和材質(zhì)缺陷等，這些缺陷會(huì)影響后道工序的成型和裝配等，影響產(chǎn)品整體質(zhì)量.探索采用基于PPRM技術(shù)的集成仿真技術(shù)監(jiān)控、預(yù)測(cè)和發(fā)現(xiàn)全過(guò)程工藝缺陷，進(jìn)而輔助工藝方案制定和實(shí)施，提高工藝水平.PPRM技術(shù)在船舶制造工藝仿真中的應(yīng)用方向主要包括以下方面：

1）考慮初始缺陷的船體結(jié)構(gòu)強(qiáng)度計(jì)算.目前，船體結(jié)構(gòu)強(qiáng)度設(shè)計(jì)中大多采用理想三維模型進(jìn)行校核計(jì)算，而實(shí)際船舶產(chǎn)品由于各階段建造工藝產(chǎn)生的缺陷，其結(jié)構(gòu)內(nèi)部存在大范圍的殘余應(yīng)力集中，這些應(yīng)力缺陷在強(qiáng)度計(jì)算模型中缺乏考慮，難以保證設(shè)計(jì)強(qiáng)度和疲勞的有效性.[9]采用PPRM技術(shù)，通過(guò)應(yīng)力后處理數(shù)據(jù)進(jìn)行模型重構(gòu)，將全工藝過(guò)程應(yīng)力考慮到船體結(jié)構(gòu)模型中，有助于提高強(qiáng)度計(jì)算的有效性.

2）分段堆放與分段吊裝工藝變形連續(xù)計(jì)算.船體分段建造完畢，往往堆放在分段堆場(chǎng)等待后續(xù)吊裝和總組.由于分段質(zhì)量一般為幾十至幾百噸，且體積大，在其堆放過(guò)程中由于重力產(chǎn)生的變形較大；在吊裝過(guò)程中同樣產(chǎn)生明顯的重力變形，影響制造精度和后續(xù)總組裝焊工藝實(shí)施.采用PPRM技術(shù)，首先進(jìn)行分段堆放重力變形仿真計(jì)算，然后利用變形后處理數(shù)據(jù)進(jìn)行模型重構(gòu)，在此基礎(chǔ)上建立吊裝變形仿真計(jì)算模型，最終實(shí)現(xiàn)分段堆放和吊裝工藝的連續(xù)變形計(jì)算，預(yù)測(cè)連續(xù)工藝變形，輔助工藝方案制定，提前探索避免變形的工藝措施，保障后序分段總組裝焊工藝順利實(shí)施.

3）大型船體分段結(jié)構(gòu)裝焊的熱彈塑性仿真計(jì)算.由于船體結(jié)構(gòu)較大，采用熱彈塑性有限元法計(jì)算焊接變形時(shí)網(wǎng)格數(shù)量多、計(jì)算量巨大，難以一次實(shí)現(xiàn)分段結(jié)構(gòu)的焊接熱彈塑性計(jì)算.[12]采用PPRM技術(shù)，在建模仿真過(guò)程中對(duì)當(dāng)前焊接工序的焊縫附近網(wǎng)格進(jìn)行網(wǎng)格細(xì)化，對(duì)后續(xù)工序的焊縫網(wǎng)格可粗化處理或忽略；而在后續(xù)焊縫仿真階段，可將之前仿真完畢的焊縫附近網(wǎng)格粗化，進(jìn)而減少各工序的計(jì)算網(wǎng)格數(shù)量，縮短計(jì)算周期，保證大型結(jié)構(gòu)計(jì)算仿真的實(shí)現(xiàn).采用該方法，一方面可以按照焊接工序分階段建立大型船體結(jié)構(gòu)的焊接仿真模型，大量減少單個(gè)仿真階段的網(wǎng)格數(shù)量，保證計(jì)算順利實(shí)施；另一方面可以考慮前道焊接工序?qū)笮蜓b焊工藝的影響，實(shí)現(xiàn)大型船體分段裝焊的全過(guò)程仿真，有利于分析結(jié)構(gòu)全過(guò)程裝焊的應(yīng)力變形狀態(tài).

4結(jié)束語(yǔ)

提出一種可用于連續(xù)工藝過(guò)程集成仿真的PPRM技術(shù)，用于預(yù)測(cè)連續(xù)工藝中前道工序?qū)嵤?duì)后道工序的影響，輔助工作人員在工序制定時(shí)從全過(guò)程工藝角度考慮方案，有助于提高制造工藝水平.

通過(guò)PPRM技術(shù)實(shí)例應(yīng)用，實(shí)現(xiàn)基于熱彈塑性理論的連續(xù)瞬態(tài)焊接工藝變形仿真，得到連續(xù)工藝變形積累結(jié)果.結(jié)果表明由于結(jié)構(gòu)肋板距離部件裝配線邊緣較近，邊緣引起的焊接變形明顯.

針對(duì)船舶制造工藝特征，探討PPRM技術(shù)在船體制造領(lǐng)域地應(yīng)用前景，明確該技術(shù)在船舶制造工藝中的應(yīng)用主研究方向，有助于進(jìn)一步深入研究PPRM.

參考文獻(xiàn)：

[1]JAROSLAV M. Finite element analysis and simulation of quenching and other heat treatment processes[J]. Comput Mat Sci， 2003， 27（3）： 313-332.

[2]張貴寶， 陳軍， 王曉方. 基于板料成型數(shù)值模擬的沖壓模具結(jié)構(gòu)分析方法[J]. 機(jī)械工程學(xué)報(bào)， 2008， 44（8）： 174-179.

ZHANG Guibao， CHEN Jun， WANG Xiaofang. Methodology of stamping die structure analysis based on sheet metal forming numerical simulation[J]. Chin J Mech Eng， 2008， 44（8）： 174-179.

[3]HU Z， LABUDOVIC M， WANG H， et al. Computer simulation and experimental investigation of sheet metal bending using laser beam scanning[J]. Int J Machine Tools & Manufacture， 2001， 41（4）： 589-607.

[4]張響， 童水光， 閆勝昝， 等. 集成鑄造數(shù)值模擬的車輪疲勞分析[J]. 機(jī)械工程學(xué)報(bào)， 2008， 44（8）： 46-49.

ZHANG Xiang， TONG Shuiguang， YAN Shengzan， et al. Integrating casting simulation into automotive wheel fatigue analysis[J]. Chin J Mech Eng， 2008， 44（8）： 46-49.

[5]王蘋， 方洪淵， 韓德成， 等. 基于熱-力耦合法的高速車地板焊接變形數(shù)值模擬[J]. 焊接學(xué)報(bào)， 2012， 33（6）： 15-18.

WANG Ping， FANG Hongyuan， HAN Decheng， et al. Numerical simulation on as-welded deflection of high speed rail floor based on thermo-mechanical coupling[J]. Trans China Welding Institution， 2012， 33（6）： 15-8.

[6]DESHPANDE A A， TANNER D W J， SUN W， et al. Combined butt joint welding and post weld heat treatment simulation using SYSWELD and Abaqus[EB/OL]. （2014-04-01）[2011-01-01]. http：//pil.sagepub.com/content/225/1/1.full.pdf+html.

[7]FISK M， LUNDBACK A. Simulation and validation of repair welding and heat treatment of an alloy 718 plate[J]. Finite Elements Anal & Des， 2012， 58（1）： 66-73.

[8]BERGLUNDA D， ALBERB H， RUNNEMALMA H. Simulation of welding and stress relief heat treatment of an aero engine component[J]. Finite Elements Anal & Des， 2003， 39（1）： 865-881.

[9]胡廣旭， 孟梅， 劉冰， 等. 基于連續(xù)瞬態(tài)計(jì)算模型的焊接殘余應(yīng)力對(duì)船體結(jié)構(gòu)強(qiáng)度影響研究[J]. 船舶力學(xué)， 2013， 17（12）： 1439-1446.

HU Guangxu， MENG Mei， LIU Bing， et al. Research on the influence of welding residual stresses on the strength of ship structures based on continuous transient model[J]. J Ship Mech， 2013， 17（12）： 1439-1446.

[10]GOLDAK J， CHAKRAVARTI A， BIBBY M. A new finite element model for welding heat sources[J]. Metallurgical Trans： B， 1984， 15（2）： 299-305.

[11]DENG D. FEM prediction of welding residual stress and distortion in carbon steel considering phase transformation effects[J]. Mat & Des， 2009， 30（2）： 359-366.

[12]陳建波， 羅宇， 龍哲. 大型復(fù)雜結(jié)構(gòu)焊接變形熱彈塑性有限元分析[J]. 焊接學(xué)報(bào)， 2008， 29（4）： 69-72.

CHEN Jianbo， LUO Yu， LONG Zhe. Analysis on welding distortion of large complicated structure by thermal elastieplastic finite element method[J]. Transaction China Welding Institution， 2008， 29（4）： 69-72.

（編輯武曉英）