【文章內容簡介】 接本身各種因素的非線性的特點產生的，而不管是幾維問題都要急需攻克的一個關鍵點，就是收斂和精度效率問題，由于大的非線性使得整個計算過程收斂十分困難，計算時間過長，王建華(上海交通)和日本大阪大學進行了三維問題的研究探討，就三維分析中的關鍵問題，收斂和精度效率問題加已分析，提出了一些解決的方法，而且對一些實例進行了實際的模擬分析，對三維分析還有待進一步的研究，以能在保證又快又準確的條件下，更好的模擬復雜構件 [2022]。熱源的選取是焊接模擬過程中的關鍵的問題，在以往的二維問題中應用高斯熱源就可以很好的模擬，但對于在三維條件下，高斯熱源有些時候就不能很好的進行模擬，蔡洪能等人提出了熱源模型雙橢球運動電弧作用，應用節(jié)點熱燴的判斷方法，進行低碳鋼(A3鋼)板試驗中的焊接溫度場中的溫度分布進行分析，最終和試驗中的結果的數(shù)據(jù)符合[16,18]。 國內外的焊接殘余應力研究的現(xiàn)狀起始于20世紀30年代的一些簡單的試驗的測量和數(shù)據(jù)的整理，開始了人們對焊接應力應變進行分析和研究，然后50～60年代通過研究人員理論經驗和數(shù)據(jù)的不斷積累逐漸的形成了一些在理論方面的權威理論作品例如梅蘭和帕爾庫斯的《由于定常溫度場而產生的熱應力》[2223]和帕爾庫斯單獨寫的《非定常熱應力》[24]，全面的總結當時人們在焊接應力應變方面取得的一些進展。起初對焊接應力應變的分析從一維的問題開始，應用圖解法分析焊接過程，其中奧凱爾布洛母等人(前蘇聯(lián)學者)在分析中加進了溫度變化對于材料屬性的影響 [13]。對于尺寸稍大的構件，僅僅的解析方法無法滿足分析的需求，人們開始向計算機程序編寫的方面推進，首次應用編寫的程序模擬一維板中堆焊由Tall等人完成 [14]。隨著一維焊接應力應變的發(fā)展和完善，人們把分析逐步的向二維領域發(fā)展，70年代初，對接焊和平板堆焊的二維應力應變分析程序就由Iwkai和Muraki編制完善，二維分析成為了可能[13]。這一年代的又一個突破就是對焊接過程的另一個關鍵因素的考慮，就是焊接過程中金屬的熔敷產生的相變組織變化 [27]。隨著焊接應力應變理論基礎的不斷的完善，80年代，人們開始注重更深層次的研究，開始研究更加準確的焊接應力應變在分布上的趨勢，通過一些計算數(shù)據(jù)，如Josefson等人定位焊和薄壁管件等焊接應力過程研究數(shù)據(jù)的分析，提出了一些精度更高的焊接應力分布趨勢和一些消減焊接應力應變的方法和措施 [28]。焊接應力應變的研究與計算機的發(fā)展密切結合，向著更加精確和細致的方面拓寬，對影響焊接應力應變的因素考慮越來越多，90年代，焊接應力應變開始考慮焊接過程當中熱源的輻射問題、焊接過程中金屬熔化產生的熔敷現(xiàn)象和焊接構件和空氣之間的熱傳遞的問題，Mahin等人在考慮了以上的因素的條件下，還考慮了溫度場和應力場的耦合關系，選用實驗的方法來矯正熱源的分布趨勢，計算出的焊接應力和用衍射中子得到的實驗結果很好的得到了吻合；接下來TTnouce，考慮了一些對焊接過程產生影響的溫度因素和焊接溫度場、應力場及其焊接過程出現(xiàn)的金屬相變潛熱耦合作用，并且構件在上述考慮條件下的焊接過程熱的本構方程，為以后在這方面的研究的各位研究人士提供了很好的可以參考的資料；焊接應力應變的分析由薄板向著厚板，單層焊接向著雙層焊接方向發(fā)展，Shim等人(美國)就在熱彈塑性平面應變有限元理論下對多層焊接的厚板件實行了計算，并且比較了不同的焊縫形式(坡口)下的焊接殘余應力，以此為依據(jù)，提取了關于厚板的焊接殘余應力的分布的趨勢 [29]?？萍嫉陌l(fā)展總是由簡單到復雜，二維焊接應力應變的完善，使得三維分析初見端倪，Chidiac(加拿大)首次在厚板的應力應變中建立了三維的有限元的溫度場分析模型，考慮了焊接過程中，由于高溫熔化，而導致的材料本身顯微晶體組織的變化和生長。近年來，，得到隨著焊接時間加長，焊縫周圍的拉伸殘余應力有減小的趨勢，而隨著焊接電壓的升高，拉伸殘余應力有加強的趨勢[30]。Cleition Carvalho Silva,Jesualdo Pereira Farias 等人用手工電弧焊對小管進行了對接焊，并用X射線測量儀器進行測量，發(fā)現(xiàn)在相似的焊接參數(shù)條件下，殘余應力得到兩種結果，即可能很大，也可能很小，但是他們還沒有確定到底是哪種參數(shù)的小改變導致焊接殘余應力的不一致性[31]。，得到焊接速度對焊接殘余應力的影響最大，提高焊接的速度有助于降低焊接殘余應力[32]，此外，激光能量大小對焊接殘余應力有巨大的影響。Dean Deng和Hidekazu Murakawa的結論是在薄板焊接過程中厚度方向幾乎不存在溫度梯度，利用大變形理論數(shù)值模擬結果與實驗值比較吻合，而利用小變形理論數(shù)值模擬與實驗結果相差較大[33]。二十世紀九十年代，汪建華等人對焊接過程應用三維數(shù)值模擬，認為焊接溫度場均屬于典型的非線性瞬態(tài)熱傳導，探討了影響收斂進度的因素算，并提出了若干解決計算精度的問題。1996年，顧強陳和紹蕃，應用熱彈塑性應力分析理論，有限元法計算了厚板焊接箱型截面的殘余應力分布[34]。2003年，陳麗敏和陳思作，根據(jù)熱彈塑性應力理論、有限元理論,用大型有限元軟件ANSYS對焊接工字型截面梁進行殘余應力分析，分析結果表明，焊接殘余壓應力的分布與截面幾何參數(shù)有關,為用有限元分析焊接工字型截面梁殘余應力提供了一種方法[35]。2006年，北京大學的楊娜，龍麗華等人以有限元分析軟件ANSYS為工作平臺,基于非線性板殼有限元理論,采用殼單元對輕型門式剛架中H型鋼楔形薄壁梁進行考慮雙重非線性的全過程分析,分析了殘余應力，對變截面和等截面H型鋼梁的相關屈曲性能的影響,殘余應力的存在與否對構件相關屈曲性能影響很大,因此,在H型鋼楔形梁的相關屈曲分析中必須考慮縱向殘余應力的影響，殘余應力峰值越大,構件延性越好,但是同時極限承載力越低[36,37]。2007年清華大學的楊文等人針對鋼板對接焊縫及腹板與翼緣角焊縫連接的工字型截面梁，研究了焊接的溫度場，殘余應力分布及殘余變形，并且通過有限元計算進一步研究了由于焊接殘余應力的存在，熱影響區(qū)內鋼材受力性能的變化[38]。同年哈爾濱工業(yè)大學的張壯南等人通過比較試驗和有限元分析,研究了焊接殘余應力對單軸對稱工字形懸伸梁和連續(xù)梁的整體穩(wěn)定承載力影響，并在其博士論文[39]中采用截面法對焊接單軸對稱工字形截面的殘余應力分布進行了測量，根據(jù)實測值，建立了接近試驗情況的單軸對稱工字形截面殘余應力分布模型[11]。同濟大學的吳蕓和張其林焊接鋁合金構件殘余應力試驗研究,通過對測試數(shù)據(jù)的分析整理,得出縱向焊接工字型截面構件殘余應力的分布情況,為進一步總結焊接鋁合金構件殘余應力的分布規(guī)律及研究殘余應力對構件承載力的影響提供了基礎。2008年，西安建筑大學的馮艷輝在其碩士論文中分析了殘余應力對焊接梁的影響[9]。2010年清華大學的班慧勇等在殘余應力實驗研究中提出了適用于Q420高強等邊角鋼的較為準確和安全的殘余應力分布模型和計算公式[40]。 焊接殘余應力數(shù)值模擬的發(fā)展趨勢利用計算機研究焊接應力和變形問題始于20世紀60年代，主要研究一維焊接應力的產生機制[41]。20世紀70年代以來由于有限元技術的發(fā)展，數(shù)值模擬方法在焊接應力變形中的研究和應用日益廣泛，但基本上是針對二維問題[4243]。在大多數(shù)情況下，焊接應力變形是三維問題，特別是現(xiàn)代焊接結構越來越大型復雜化，而且還存在許多不確定因素,雖然有些可以簡化為二維問題分析[4445]，但是在實際工程中真正可以簡化的例子并不多，因此三維焊接應力變形模擬是必然趨勢。由于焊接過程的復雜性以及焊接結構三維數(shù)值模擬中自由度大、計算效率低、計算精度難以保證等特征。通過近年來三維焊接應力變形數(shù)值模擬研究現(xiàn)狀，取得的進展和現(xiàn)有關鍵問題的分析可知，提高模擬計算的精度和效率是今后焊機模擬的關鍵問題[40]：(1)采用合適的熱源模型和動態(tài)可逆自適應網(wǎng)格生成技術可以減少網(wǎng)格和結點數(shù)目，在保證計算精度的同時，可大大提高運算效率,但其技術細節(jié)及準確性等迄今尚無報道；(2)通過適當調整材料高溫性能參數(shù)有利于有限元解的收斂性，提高計算效率。利用單元“死活”技術描述多層焊及焊縫金屬的熔敷；(3)利用并行計算技術可提高焊接數(shù)值模擬的計算效率,開發(fā)高性能的并行程序和分布處理系統(tǒng)，是今后發(fā)展的趨勢。并行操作和高性能數(shù)據(jù)交換開關是分布式并行系統(tǒng)的開中重要的兩個方面；(4)在焊接物理模擬過程中采用相似理論可以有效減小模擬件的幾何尺寸，減少結構的自由度和計算工作量，盡管結構的尺寸被縮小，位移、溫度應力的結果也很好吻合，節(jié)省計算時間的目的沒有達到，所以還需要進一步的簡化模型。應用有限元進行三維模擬還存在著這樣或那樣的問題,將三維焊接過程數(shù)值模擬技術廣泛運用于工程實際，仍然需要進一步提高計算效率和保證計算精度。 課題的研究內容本課題是單軸對稱焊接工字梁殘余應力數(shù)值模擬及整體穩(wěn)定性分析，采用理論和有限元分析相結合,循序漸進的技術路線和步驟，主要研究內容如下：(1)建立可行的焊接工字形梁焊接處的三維溫度場及殘余應力的動態(tài)模擬分析方法。(2)通過變化參數(shù)對單軸對稱焊接工字梁進行數(shù)值分析，分析參數(shù)變化對焊后殘余應力的分布規(guī)律的影響。 課題研究技術路線本課題是單軸對稱焊接工字梁殘余應力數(shù)值模擬，采用理論和有限元分析相結合,循序漸進的技術路線和步驟：首先收集資料進行理論分析，結合已有實驗模型對構件作初步設計，勾勒出具體的模型，應用Ansys中的熱分析和靜力分析,對所建模型進行三維有限元模擬，求出模型焊接殘余應力的分布與實驗結果相比較，驗證分析的正確性。然后建立單軸對稱焊接工字梁模型，計算模型在不同參數(shù)下的焊接殘余應力的分布規(guī)律。第二章 焊接過程有限元分析理論 焊接過程有限元分析特點采用時間和空間有限元法，相同程度構件細節(jié)采用彈性構件分析，探討焊接過程中，焊料和焊件之間的熱力關系、焊接殘余應力應變，以下為焊接過程中應用有限元分析的特點:(1)建立三維數(shù)值模擬模型，至少焊縫和焊縫附件是如此，用來考慮表面和內部由于冷卻產生的不用影響；(2)溫度場模擬是典型的瞬態(tài)非線性，材料的熱物理屬性與溫度有關，產生的溫度梯度與時間和位置相關，但又是極不相同； (3)局部熱產生的歷史和力學的應力應變歷史將決定局部焊接過程中材料的瞬態(tài)行為；(4)焊接過程中材料將發(fā)生熔敷以及凝固，將改變材料的顯微組織而且在焊后將會改變構件的連接；(5)連續(xù)介質的概念將受到懷疑，因為可能在臨界點發(fā)生缺陷和裂紋；即使是現(xiàn)在，也無法很好的解決收斂和誤差估計困難的難題。另一個工業(yè)和加工過程中存在的問題是，在分析前需要確定的很多材料隨溫度變化的性能參數(shù)，現(xiàn)在高溫下的性能參數(shù)數(shù)據(jù)非常少，基本屬于空白狀態(tài)。 焊接有限元模型的簡化焊接過程是一個涉及到多個學科共同作用和多個物理現(xiàn)象的復雜過程。，在焊接有限元分析中所要輸入和輸出的基本的參數(shù)。并在分析中突出了相變對焊接過程的影響，焊接溫度場應力場及顯微組織之間的影響。，強烈影響用實箭頭表示，較弱影響用虛線箭頭表示。，在做焊接過程有限元分析時，由于焊接溫度場和顯微組織變化對焊接應力場的影響要大于后者對后者的影響，所有只需要考慮前者對后者的影響，就可以滿足要求即考慮單向耦合。本文僅研究了單向耦合問題。此外，本文還考慮了相變潛熱對溫度場的影響。熱源模型，熔池模型，熱力學性能化學成分，焊接參數(shù)，晶粒尺寸瞬態(tài)溫度場相變行為力學特性應力與變形顯微組織的形成和發(fā)展熱 應 力變 形 熱相變潛熱相變相變應力動力學相變 顯微組織變化、溫度場、應力場互相影響 Changes in microstructure、the temperature field、stress field interaction金屬學顯微組織狀態(tài)場由熱應力引起的相變力 學應 力 和 變 形熱 力 學溫 度 場相變潛熱顯微組織轉變熱應力 焊接熱應力、相變潛熱和溫度的相互耦合關系 Welding heat stress、latent heat and temperature relationship of mutual coupling 焊接溫度場分析理論 焊接過程傳熱的基本形式和方程焊接構件的局部受熱，造成焊接中產生，根據(jù)研究結果顯示，焊接熱傳導主要研究的問題是隨時間變化的溫度分布，加以考慮對流和輻射對焊接過程的影響，隨著焊接過程的局部快速的加熱和冷卻，焊接構件的溫度和隨溫度變化的材料屬性在空間和時間上都發(fā)生著強烈的改變，相變潛熱金屬熔敷等現(xiàn)象都伴隨著發(fā)生，所以說焊接過程是典型的瞬態(tài)非線性傳導熱分析[8,48]。溫度場三維熱傳導控制方程如下： ()式中，—內熱源；—材料比熱容；—溫度場溫度分布函數(shù)；—材料密度；—傳熱時間；—材料導熱系數(shù)。在給定初始條件的情況下可以對上述的泛定方程求得定解。過程開始時整個區(qū)域已經具有的已知溫度，即：或 ()式中，—表示均勻初始溫度的已知常數(shù)；—表示不均勻初始溫度的已知函數(shù)以下幾類為焊接溫度場中常用的邊界條件：第一類邊界條件，已知邊界上的分度分布 ()第二類邊界條件，已知邊界上熱流密度函數(shù)分布 ()第三類邊界條件，已知邊界上物體與周圍介質熱交換函數(shù)分布 ()式中：—單位面積上外熱源；—表面熱交換系數(shù)；—周圍介質溫度分布函數(shù)；,—各軸邊界外法線方向余弦。 非線性瞬態(tài)