外文翻譯--對塔式起重機臂三維裂紋的有限元斷裂分析 中文版

1、<p>　　對塔式起重機臂三維裂紋的有限元斷裂分析</p><p>　　曹宗杰,匡振邦和李昌平</p><p>　　中國上海交通大學工程力學系， 200240</p><p>　　中國長春航空空軍大學機械工程系，130022</p><p>　　zjcao96@sina.com,zbkiang@mail.sjtu.edu.cn,Li

2、changpingl@sina.com</p><p>　　摘要：本文介紹了在三維裂紋前端位移域的局部解和構造三維裂紋問題的一種新方法，和提出了一種新方法用于計算三維裂紋問題的應力強度因子。用現(xiàn)有的方法，對結構強度與縱向閉合裂紋軸向張力進行了分析。裂紋擴展是可能的，因為軸壓塔式起重機主臂的關節(jié)附近有環(huán)拉伸應力。應力和應力強度因子用于塔式起重機主臂關節(jié)計算。探討了焊接殘余應力和應力強度因子的影響。</p>

3、;<p>　　關鍵詞：應力強度因子；殘余應力奇異有限元；裂紋</p><p><b>　　1引言</b></p><p>　　電力系統(tǒng)塔式起重機是塔式起重機類型中最重的的一個。一些災難性事故發(fā)生時塔式起重機都是在使用工程中的。磁測量，封閉式縱向裂縫在塔式起重機主臂的表面經常發(fā)現(xiàn)。裂縫影響結構的強度，因此對縱向閉合裂紋塔式起重機主臂的強度必須進行分析。因為

4、裂紋尖端，一般方法很難準確地確定裂紋如何影響塔式起重機結構強度。許多研究人員在裂紋尖端的奇異領域提出了各種模擬和處理。每一種方法都有其優(yōu)點和缺點存在。他們中的大多數(shù)是需要密集的網(wǎng)格計算的有限元方法。一些更為精確的元素，如雜交/混合有限元，有更精確的結果，但是它很難控制結果的穩(wěn)定性。因此需要在工程上應用一種特殊的算法。</p><p>　　奇異有限元求解應力強度因子裂紋問題時有特殊的優(yōu)勢。在本文中的三維裂紋的前端位

5、移域解決方案，介紹了一種構建二維裂紋彈性問題和三維裂紋問題的新的奇異單元，并提出了一種新方法用于計算應力強度因子的三維裂紋問題。此方法不僅可以節(jié)省內存和計算時間，也增加了對涉及到應力強度因子有限域問題的計算解的精度。與現(xiàn)有的方法相比，對縱向閉合裂紋的結構強度軸向張力進行了分析。對塔起重機的主臂關節(jié)應力和應力強度因子進行計算。探討了焊接殘余應力和應力強度因子的影響。</p><p><b>　　2奇異有限

6、元方法</b></p><p>　　由于應力裂紋奇異性的存在，常規(guī)有限元應力域精確的表達強度是非常困難的。因此，任意裂紋形狀的三維物體，局部位移場，應力場和應變場可以表示為如下形式：</p><p><b>　?。?）</b></p><p><b>　　（2）</b></p><p>

7、<b>　?。?）</b></p><p>　　這里的、和分別是分布函數(shù)矩陣的局部應力場，應變場和位移場。的應力強度因子（SIF）的陣列在裂紋前沿的任意一點，二次插值函數(shù)模擬元素中裂紋前沿的應力強度因子是該點放置在裂紋前端的功能。然后下面的公式可給出：</p><p><b>　?。?）</b></p><p><b

8、>　　（5）</b></p><p><b>　?。?）</b></p><p>　　其中L是應力強度因子的二次插值函數(shù)矩陣；是應力強度因子在單元節(jié)點的陣列，是在裂紋前沿的應力強度因子的節(jié)點數(shù)組。局部解式（1）-（3）為等參單元，對任意形狀的三維裂紋單元位移場如下：</p><p><b>　?。?）</b&g

9、t;</p><p>　　其中q是節(jié)點位移，并用奇異的元素數(shù)。M是M的節(jié)點單元邊界，L是元素邊界節(jié)點。根據(jù)應變與位移之間的關系，為，可以得到以下方程</p><p><b>　　（8）</b></p><p>　　其中和為應變向量應力矢量；是微分算子；B是應變矩陣；是奇異應變矩陣；D是彈性系數(shù)矩陣。根據(jù)應力和應變之間的關系，可以得到以下方程：&

10、lt;/p><p><b>　　（9）</b></p><p>　　3三維裂紋應力強度因子</p><p>　　根據(jù)最小勢能原理，可以得到以下方程</p><p><b>　　（10）</b></p><p>　　當是整個研究區(qū)時，是奇異的區(qū)。f和分別是表面力的奇異單元。將式（7

11、）-（9）代入式（10），可以得到以下方程：</p><p><b>　　(11)</b></p><p>　　其中Q是全球位移陣列。F和分別是整體區(qū)域研究和奇異區(qū)域邊界力向量。K是整體剛度矩陣。是耦合的整體剛度矩陣元素之間的定期和奇異的元素。是奇異單元的整體剛度矩陣。通過分別改變Q和之間的獨立性，可以得到以下方程</p><p><b&

12、gt;　　(12)</b></p><p><b>　　(13)</b></p><p>　　當時，式（13）可以改寫為下列形式：</p><p><b>　　(14)</b></p><p>　　替代式（14）代入式（12），可以得到以下方程：</p><p>

13、<b>　　(15)</b></p><p>　　當 (16)</p><p><b>　　(17)</b></p><p><b>　　4數(shù)值例子</b></p><p>　　塔式起重機是桁架結構，其主臂是無縫鋼管（）是20鋼材料。裂紋的長

14、度為5～500mm，裂縫的深度為0.5～2.5mm。計算模型的建立應用二十節(jié)點三維有限元素。關節(jié)臂有限元網(wǎng)格與附近的封閉式縱向裂縫的主要危險如圖1所示。裂縫計算長度為500mm。分別計算裂縫深度，為1.5,2.5，3和3.5毫米。計算負荷（）等于。出于焊接殘余應力對主臂與封閉式縱向流的考慮，焊接殘余應力（）相當于等于0.5。泊松分布，材料的楊氏模量的比分別等于0.3和0.07MN /。環(huán)拉伸應力在不同的徑向截面無焊接殘余應力，如圖2所

15、示。環(huán)拉伸應力不同軸向截面無焊接殘余應力如圖3和圖4所示。通過上述計算，塔式起重機主臂的關節(jié)部分有環(huán)的拉伸應力在可能的流環(huán)附近傳播。</p><p>　　從圖2~4可知，主臂縱向閉合裂紋的環(huán)拉伸應力降低時有焊接殘余應力。在不同的情況下，主臂關節(jié)部分應力強度因子的三維裂紋附近無焊接殘余應力和對焊接殘余應力進行了計算，如圖5,圖6所示。從圖5可以看出，最大應力強度因子等于時裂紋的厚度為3.5mm，關節(jié)處發(fā)生較少， 在

16、10cm左右。單邊緣裂紋擴展門檻值，它的長度是0.025到5毫米，是從4.0至的20鋼材料。主臂關節(jié)附近的最大應力強度因子小于單邊緣裂紋擴展門檻值。如圖6，最大應力強度因子等于時，裂紋深度等于3.5mm。最大應力強度因子與焊接殘余應力小于無焊接殘余應力。</p><p><b>　　5結論</b></p><p>　　本文對三維裂紋前端位移域的解決方案進行了介紹和三維

17、裂紋問題的一種新的奇異單元的構造，和提出了一種新方法，計算三維裂紋問題的應力強度因子。與用現(xiàn)有的方法，對結構強度與縱向閉合裂紋軸向張力進行了分析。這是可能的裂紋擴展，因為是有環(huán)的拉伸應力下的軸壓塔式起重機主臂的關節(jié)附近的應力。焊接殘余應力對焊接材料強度有害，但焊接殘余應力使主臂環(huán)拉伸應力下降。在同一時間，應力強度因子的下降，改進了塔式起重機主臂的安全性。</p><p><b>　　參考文獻</b

18、></p><p>　　[1] 王志超，張理蘇,《一種新的準諧調有限元列式方法及在斷裂力學中的應用》1990.37（6）1195-1201</p><p>　　[2] O.C.欽科維奇（Zienkiewicz O.C.），《有限元法》,麥格勞希爾圖書有限公司第三版：669，1977</p><p>　　[3]Barsoum RS.Int J.Numer.Mer

19、h.Engng.1976,10:25~37</p><p>　　[4] Abdel Wahab MM, De Roeck G A Finite Element Solution for Elliptical Cracks Using the ICCI Method^ Engineering Fracture Mechanics, 1996,53:518-526</p><p>　　[5]

20、曹宗杰《三維奇異準諧調單元理論研究及其數(shù)值分析》，1994。</p><p>　　[6] Raju Ls, Newman JC_ Stress-intensiEy Factors for A Wide Range of Semi-eUiptical Surface Cracks in Finite- thickness Plates, Engineering Fracture Mechanics, 1979,11