胡繩蓀，蔣秀曄，申俊琦，朱麗娜

(1. 天津大學材料科學與工程學院，天津 300072；2. 天津市現(xiàn)代連接技術(shù)重點實驗室，天津 300072)

基于ANSYS二次開發(fā)的焊接溫度場前處理系統(tǒng)

胡繩蓀1,2，蔣秀曄1,2，申俊琦1,2，朱麗娜1,2

(1. 天津大學材料科學與工程學院，天津 300072；2. 天津市現(xiàn)代連接技術(shù)重點實驗室，天津 300072)

要：針對大型應用軟件ANSYS在進行有限元分析時存在重復建模缺點，基于VC與ANSYS的參數(shù)化設計語言APDL，對有限元前處理部分進行二次開發(fā)，編制出了管道焊接結(jié)構(gòu)溫度場的前處理系統(tǒng)．該系統(tǒng)封裝了多種管道焊接結(jié)構(gòu)建模以及不同焊接工藝的溫度場分析過程的參數(shù)化命令流模板．通過友好直觀的程序界面的操作，即可實現(xiàn)多種管道焊接結(jié)構(gòu)的參數(shù)化建模、溫度場數(shù)值模擬，從而提高了計算效率，保證了計算精度．并通過管道多層多道定位焊實例驗證了該前處理系統(tǒng)的準確性及可用性．

ANSYS二次開發(fā)；溫度場數(shù)值模擬；前處理系統(tǒng)

隨著社會經(jīng)濟的發(fā)展，地下管道已經(jīng)成為現(xiàn)代工業(yè)生產(chǎn)和城鎮(zhèn)生活的大動脈，不僅應用在水、油、氣(汽)、煤的輸送，還應用在通信、供電、交通和排水等方面．焊接作為管道結(jié)構(gòu)的重要連接技術(shù)，在管道結(jié)構(gòu)工程建設中發(fā)揮著重要作用．為保證結(jié)構(gòu)的尺寸要求及焊接接頭的質(zhì)量，需要通過數(shù)值模擬技術(shù)準確模擬實際焊接過程中焊接構(gòu)件的殘余應力和變形[1-3].傳統(tǒng)數(shù)值模擬技術(shù)在高校及科研機構(gòu)發(fā)展較為成熟，但由于操作復雜、難以掌握、對使用人員的專業(yè)性要求較高、焊接模擬分析時重復性差等特點，在工業(yè)工程中難以推廣[4-6]．為了解決以上問題，一些企業(yè)和研究院所開展了關(guān)于ANSYS二次開發(fā)的研究．其二次開發(fā)主要依賴于ANSYS提供的參數(shù)化程序設計語言APDL[7-10]，不僅可以實現(xiàn)參數(shù)化分析，還可以根據(jù)需要改變程序以滿足特定的設計和分析．此外，焊接結(jié)構(gòu)的溫度場、變形及殘余應力分析的技術(shù)路線已經(jīng)相對比較成熟，可以形成較專業(yè)的分析模塊．基于以上研究現(xiàn)狀，筆者通過研究參數(shù)化的分析過程，并運用VC和ANSYS二次開發(fā)工具APDL建立管道焊接結(jié)構(gòu)溫度場數(shù)值模擬的專業(yè)前處理分析系統(tǒng)．通過友好的菜單操作，實現(xiàn)常見焊接結(jié)構(gòu)的建模、熱分析過程．

1 前處理系統(tǒng)總體設計

該前處理系統(tǒng)采用了可視化編程和結(jié)構(gòu)化設計方法，主要包括兩大模塊：計算模塊和界面交互模塊．計算模塊是保證系統(tǒng)精度的核心，包括各種焊接結(jié)構(gòu)的APDL命令流模板文件；界面交互模塊利用了VC對ANSYS前處理功能進行二次開發(fā)，生成可視化界面，包括幾何建模、網(wǎng)格劃分、材料性能、焊接參數(shù)和提交數(shù)據(jù)5大功能界面，如圖1所示.

圖1 前處理系統(tǒng)總體設計示意Fig.1 Layout of pretreatment system

該系統(tǒng)的開發(fā)環(huán)境包括：基于Visual C++6.0和ANSYS12.1的系統(tǒng)開發(fā)平臺以及基于Microsoft Access2.3的數(shù)據(jù)庫管理系統(tǒng)．

1.1 界面交互模塊

界面交互模塊負責搜索計算模塊中的匹配命令流模板，并將接收到的參數(shù)化建模尺寸參數(shù)、熱分析參數(shù)填寫在模板中的相應位置．此外，通過VC調(diào)用ANSYS，將新生成的建模命令流傳遞到ANSYS進程中，從而完成自動建模和焊接溫度場的數(shù)值模擬．

為了滿足對不同管道焊接結(jié)構(gòu)模型的有限元分析，該前處理系統(tǒng)設計了多種焊接結(jié)構(gòu)，包括平板對接、管對接、管道法蘭接頭等結(jié)構(gòu)．并對不同的結(jié)構(gòu)模型設計了多種坡口形式的焊接參數(shù)的選擇，包括I、V、X．此外，該用戶界面具有參數(shù)修改和保存、圖形顯示等功能，如圖2(a)所示．在保證計算精度的同時盡量降低計算時間，本系統(tǒng)采用過渡網(wǎng)格劃分形式，如圖2(b)所示．通過調(diào)用ANSYS進行建模和網(wǎng)格劃分，將模型圖片顯示在ANSYS計算界面中．用戶通過觀察網(wǎng)格劃分的疏密，調(diào)整參數(shù)設置，從而實現(xiàn)重復建模．

在焊接過程中，涉及到不同母材、焊縫等不同焊接材料的多種性能，而且焊接材料的性能是一個與溫度有關(guān)的參數(shù)，故其參數(shù)量較大．為此，材料性能界面建立隨溫度變化的動態(tài)材料庫，可以對其進行增刪查改功能，如圖2(c)所示．為了直觀地反映材料性能隨溫度的變化情況，增加了函數(shù)曲線繪制功能．

由于熱源模型在溫度場數(shù)值模擬時起著重要的作用，通過對常用焊接溫度場熱源模型的研究，該系統(tǒng)選取了高斯熱源和雙橢球熱源2種熱源形式，如圖2(d)所示．高斯熱源形式適用于二維溫度場數(shù)值模擬，雙橢球體常作為體熱源應用于三維溫度場數(shù)值模擬．

為了準確進行不同焊接工藝的有限元分析，焊接工藝參數(shù)模塊考慮了多層多道的焊接參數(shù)，并利用Access數(shù)據(jù)庫建立了相應的焊接參數(shù)數(shù)據(jù)庫，如圖2(e)所示．由于定位焊和焊接順序?qū)附討ψ冃蔚挠绊懞艽?，所以該模塊還考慮了定位焊的個數(shù)、冷卻時間、定位焊的長度、不同焊接結(jié)構(gòu)的焊接順序，如圖2(f)所示．

1.2 計算模塊

計算模塊是參數(shù)化命令流模板的集合，負責向系統(tǒng)提供命令流模板，包括管對接、平板對接、管道法蘭焊接的參數(shù)化建模命令流模板、熱分析命令流模板．它是保證系統(tǒng)精度的核心，直接影響著系統(tǒng)的準確性及普遍使用性．參數(shù)化程度越高，準確性、普遍使用性也相應提高；模板的算法越貼近實際，準確性越高．管道焊接結(jié)構(gòu)有限元分析APDL命令流文件是經(jīng)過實驗驗證核實的可靠文件，可以保證其準確性的要求．

計算模塊運用了多種模擬技巧，除考慮材料的熱物理性能、過渡網(wǎng)格劃分方式、不同的熱源形式以及常用的生死單元技術(shù)，還根據(jù)實際焊接經(jīng)驗增加了一些技巧．例如，增大熔池金屬處單元的導熱系數(shù)，模擬焊接過程中熔滴對熔池的沖刷作用，在模擬過程中考慮定位焊的工藝，以貼合實際．

圖2 典型界面模塊Fig.2 Typical interface modules

2 接口技術(shù)的實現(xiàn)

2.1 VC與ANSYS接口技術(shù)

VC與ANSYS的接口技術(shù)，主要通過Create Process函數(shù)創(chuàng)建進程去實現(xiàn)ANSYS的調(diào)用過程．程序啟動后，系統(tǒng)調(diào)用CreateProcess創(chuàng)建管理進程的內(nèi)核對象，將系統(tǒng)的程序的代碼和界面輸入的數(shù)據(jù)加載到新開辟的虛擬地址空間，從而成功創(chuàng)建新進程和主線程．lpApplicationName和lpzCommand Line參數(shù)是命令行字符串．通過該字符串可設置待打開的可執(zhí)行文件名和傳遞給新進程的命令行．當系統(tǒng)找到可執(zhí)行文件后將界面數(shù)據(jù)賦值給計算模塊的APDL命令流模板，并將賦值后的可執(zhí)行文件的代碼和數(shù)據(jù)映射到新進程的地址空間中，生成新進程的命令行，從而調(diào)用ANSYS讀取命令行進行數(shù)值模擬．圖3所示為CreateProcess函數(shù)進程控制界面數(shù)據(jù)傳遞以及ANSYS的計算分析過程．

圖3 CreateProcess進程控制示意Fig.3Block diagram of process control by CreateProcess function

在調(diào)用ANSYS時，通過命令行字符串實現(xiàn)與它的控制和通信．ANSYS識別的命令字符有：-b批處理方式，-g圖形交互方式，-dir設置工作目錄，-i輸入文件的全路徑和名稱，-o指定輸出文件的全路徑和名稱．

2.2 VC與Access接口技術(shù)

軟件中使用數(shù)據(jù)庫來存儲材料的物理性能、多層多道焊的焊接工藝參數(shù)．針對材料庫中每種物理性能建立一個表，每個表中包含材料編號、溫度、性能值3字段．此外，還建立了索引表，包含材料名稱和材料編號的對應關(guān)系，如表1所示．

表1 材料庫索引Tab.1 Index of material database

程序中使用ADO方式操作Access數(shù)據(jù)庫．VC與Access的連接方法如下：引入msado15.dll的動態(tài)連接庫，創(chuàng)建ADO智能指針，初始化COM環(huán)境，連接數(shù)據(jù)庫，實現(xiàn)對表的增刪改查，釋放相應資源．

3 溫度場數(shù)值模擬實例驗證

3.1 實例概況

為了驗證系統(tǒng)的準確性及可用性，以X形坡口的管道對接焊為例，管道長為150,mm，外徑為150,mm，壁厚為12,mm．為了保證溫度場的計算精度和計算速度，單元類型選用節(jié)點數(shù)較少的八節(jié)點單元，選取焊縫區(qū)與母材區(qū)之比為1∶9的大模型過渡網(wǎng)格劃分方式．本例所采用的材料為Q345，從材料性能模塊的材料庫中調(diào)用．管道焊接熱源模型選用雙橢球熱源，熱源尺寸參數(shù)根據(jù)焊縫的熔合線初步選取[11-12]．本例采用多層多道定位焊，焊接層數(shù)及道數(shù)如圖4所示，焊接工藝條件及參數(shù)如圖2(e)和圖2(f)所示．

3.2 模擬結(jié)果及驗證分析

將管道對接焊的幾何尺寸參數(shù)、材料熱物理性能和焊接工藝參數(shù)等輸入該前處理系統(tǒng)界面中，系統(tǒng)通過給命令流模板賦值生成可執(zhí)行的命令流文件，并調(diào)用ANSYS進行后臺計算．計算完畢后，通過ANSYS的后處理程序查看其計算結(jié)果．圖5為管道網(wǎng)格劃分模型，圖6為管道外表面最后一道焊縫的溫度場云圖．從圖6可以看出，熔池處的最高溫度為1,671,℃，溫度場的形貌為一個橢圓狀．焊接過程中溫度場基本上屬于準穩(wěn)態(tài)，整個過程溫度場變化也比較符合實際焊接變化情況．

圖4 X型坡口多層多道焊示意Fig.4Schematic plot of multi-layer and multi-pass welding of X groove

圖5 管道網(wǎng)格劃分模型Fig.5 Pipe grid partition model

圖6 溫度場云圖Fig.6 Nephogram of temperature field

進一步分析驗證該溫度場的合理性，通過POST1繪制了節(jié)點A、B的全過程熱循環(huán)曲線，并運用POST26繪制了垂直焊縫方向的溫度分布曲線．如圖7所示，直線L1為管道外徑上垂直于焊縫的路徑線段，點A為管道內(nèi)徑上第一道焊縫處的節(jié)點，點B為管道外徑上最后一道焊縫處的節(jié)點．

圖8為最后一道焊時垂直焊縫方向的溫度分布，溫度最高點所對應的位置是最后一道焊縫的熱源中心，熱源中心溫度為1,671,℃；在距熱源中心3,mm的范圍內(nèi)的焊縫溫度在1,450,℃以上，處于熔化狀態(tài)，這與所定義的熱源尺寸寬度相吻合；在前7.9,mm內(nèi)溫度下降很快，溫度梯度很大；在7.9,mm以后，溫度下降緩慢，維持在240,℃以下．由此可知，該曲線與文獻[13-15]中所標注的軸向溫度距離曲線的趨勢基本保持一致．

圖7 路徑和節(jié)點示意Fig.7 Schematic plot of paths and nodes

圖8 垂直焊縫方向的溫度分布Fig.8 Temperature distribution of vertical weld direction

由圖9(a)可知，該曲線是典型的多層多道焊熱循環(huán)曲線．點A所處焊道為第1層第1道，最高峰之后的3個小峰即由于后焊的3道對點A的熱作用而產(chǎn)生．點B所處焊道為最后1層最后1道，對比圖9(a)和圖9(b)中2條曲線的時間軸可知，點B所對應的時間正好是點A最后1個小峰所在位置，兩者相吻合，并且與本例中設置的X型3層4道焊的焊接工藝保持了高度一致性．

圖9 點A和點B全過程熱循環(huán)曲線Fig.9 Thermal cycling curves of points A and B during whole welding process

4 結(jié) 論

(1) 基于VC強大的界面交互功能，該前處理系統(tǒng)界面友好、交互性強，降低了對分析者的要求，便于操作和分析．此外，VC強大的數(shù)據(jù)庫連接功能，便于后續(xù)焊接工藝參數(shù)、材料庫等大量數(shù)據(jù)的保存和分析．并且，運用其中的CreateProcess函數(shù)的進程控制可以實現(xiàn)VC與ANSYS的調(diào)用過程．

(2) 該前處理系統(tǒng)采用參數(shù)化建模和參數(shù)化分析的方法，通過ANSYS的APDL語言建立不同焊接結(jié)構(gòu)和工藝的管道焊接溫度場的命令流文件，將命令流文件中的參數(shù)變量化，建立相應的模板文件庫．該模板文件庫都是通過工程實踐驗證的命令流文件，從而保證了計算的可用性和準確性．

(3) X型管道多層多道焊溫度模擬的實驗驗證結(jié)果表明，該前處理系統(tǒng)模擬結(jié)果符合多層多道焊的實際情況．

［1］ Habib L M，Kleen U，Otremba F. Numerical simulation of weld residual stresses and countermeasures in austenitic steel piping[C] // International Conference on Nuclear Engineering. Nice，F(xiàn)rance，1997：25-29.

［2］ Basavaraju Chakrapani. Simplified analysis of shrinkage in pipe to pipe butt welds[J]. Nuclear Engineering and Design，2000，197(3)：239-247.

［3］ Zhang Y M，Kovacevic R，Wu L. Dynamic analysis and identification of gas tungsten arc welding process for weld penetration control[J]. Journal of Engineering for Industry，1996，118(1)：123-136.

［4］ Zhang Tao. Ansys second development method and application in transient thermal simulation[C] // International Conference on Computer Science and Software Engineering. Wuhan，China，2008：643-645.

［5］ Sun Yuantao，Wang Shaomei. ANSYS second development method and application in computer aided engineering[C] // International Conference on Computer- AidedIndustrial Design and Conceptual Design. Hangzhou，China，2006：17-19.

［6］ Li Bin，Zhu Bo. Parameterized computing system on structure of lift scaffolding attached to the wall[C] // International Conference on Mechanical and Electronics Engineering. Hefei，China，2011：3972-3977.

［7］ 吳建軍. 基于VC++與ANSYS單纖維拔出實驗模擬分析系統(tǒng)的研究與開發(fā) [D]. 上海：上海交通大學機械工程學院，2009.

Wu Jianjun. Development of Simulation and Analysis System of Single Fibre Pull-Out Test Based on VC++ and ANSYS[D]. Shanghai：School of Mechanical Engineering，Shanghai Jiaotong University，2009(in Chinese).

［8］ 朱 彬. 基于Ansys的鋁合金TIG焊溫度場應力場數(shù)值模擬及二次開發(fā) [D]. 鎮(zhèn)江：江蘇大學材料科學與工程學院，2007.

Zhu Bin. Numerical Simulation of TIG Welding Temperature and Stress Field for Aluminum Alloy Based on Ansys and Secondary Development[D]. Zhenjiang：School of Material Science and Engineering，Jiangsu University，2007(in Chinese).

［9］ 朱麗娜. 海洋工程焊接結(jié)構(gòu)應力及變形分析系統(tǒng)的研究與開發(fā) [D]. 天津：天津大學材料科學與工程學院，2010.

Zhu Lina. The Research of Welding Residual Stress and Deformation Analysis System for Welding Structure Used on Ocean Engineering[D]. Tianjin：School of Materials Science and Engineering，Tianjin University，2010(in Chinese).

［10］ 邢艷洪. 基于ANSYS二次開發(fā)的電梯參數(shù)化有限元分析系統(tǒng) [D]. 沈陽：東北大學機械工程學院，2006.

Xing Yanhong. Parameterized Lift Finite Element Analysis System Based on Second Development of ANSYS[D]. Shenyang：School of Mechanical Engineering，Dongbei University，2006(in Chinese).

［11］ 王 煜，趙海燕，吳 甦，等. 高能束焊接雙橢球熱源模型參數(shù)的確定 [J]. 焊接學報，2003，24(2)：67-70.

Wang Yu，Zhao Haiyan，Wu Su，et al. Shape parameter determination of double ellipsoid heat source model in numerical simulation of high energy beam welding[J]. Transactions of the China Welding Institution，2003，24(2)：67-70(in Chinese).

［12］ 曾 志，王立君，王 月，等. 5A06鋁合金間斷焊溫度場的數(shù)值模擬 [J]. 天津大學報，2008，41(7)：849-853.

Zeng Zhi，Wang Lijun，Wang Yue，et al. Numerical simulation of temperature distribution during discontinuous welding of 5A06 aluminum alloy[J]. Journal of Tianjin University，2008，41(7)：849-953(in Chinese).

［13］ 宋立新，王 勇，韓 濤，等. 管線鋼在役焊接多道焊的數(shù)值模擬 [J]. 壓力容器，2007，24(11)：18-21.

Song Lixin，Wang Yong，Han Tao，et al. Numerical modeling of the multi-pass welding of in-service welding onto the pipeline steel[J]. Pressure Vessel Technology，2007，24(11)：18-21(in Chinese).

［14］ Deng D， Kiyoshima S. FEM prediction of welding residual stresses in a SUS304 girth-welded pipe with emphasis on stress distribution near weld start/end location[J]. Computational Materials Science，2010，50(2)：612-621.

［15］ Deng D， Murakawa H. Finite element analysis of temperature field，microstructure and residual stress in multi-pass butt-welded 2.25Cr-1Mo steel pipes[J]. Computational Materials Science，2008，43(4)：681-695.

Pretreatment System of Welding Temperature Field Based on Secondary Development of ANSYS

Hu Shengsun1,2，Jiang Xiuye1,2，Shen Junqi1,2，Zhu Lina1,2
(1. School of Materials Science and Engineering，Tianjin University，Tianjin 300072，China；2. Tianjin Key Laboratory of Advanced Joining Technology，Tianjin University，Tianjin 300072，China)

When carrying out numerical analysis by large application software ANSYS，the programs usually need to be reworked or modified if the working conditions have changed. In order to solve the problem，a friendly pretreatment system of finite element，based on the VC platform and APDL language，was built by the secondary development of ANSYS. The pretreatment system was used to simulate the temperature field in pipelines，which encapsulated the modeling process of many different kinds of pipe joint forms，as well as parameterized command flow templates of temperature field simulation with different welding processes. Through friendly and intuitive interface operation，parameterized modeling and temperature field simulation of various welded structures can be accomplished，thus improving the computational efficiency and ensuring the calculation precision. The analytical methods turns out to be accurate and reliable as compared with the experimental results of multi-layer and multi-pass butt-welded pipes for positioning welding.

secondary development of ANSYS；temperature field numerical simulation；pretreatment system

TG456.7

A

0493-2137(2013)11-1039-06

DOI 10.11784/tdxb20131115

2012-04-19；

2012-08-15.

胡繩蓀（1956— ），男，教授，huss@tju.edu.cn.

申俊琦，shenjunqi@tju.edu.cn.