董鵬鵬 賀啟林 王海洲 王儒文

平板錐形膜盒應(yīng)力強度因子計算

董鵬鵬 賀啟林 王海洲 王儒文

(北京宇航系統(tǒng)工程研究所 深低溫技術(shù)研究北京市重點實驗室，北京 100076)

基于梁理論和薄板理論分別建立了平板錐形焊接膜盒的理論計算模型，推導(dǎo)了膜片焊菇處裂紋尖端的應(yīng)力強度因子在軸向外載和壓力下的計算公式。為評估理論計算公式的準(zhǔn)確性，基于ABAQUS的圍線積分方法計算了裂紋尖端的應(yīng)力強度因子。結(jié)果表明薄板模型的應(yīng)力強度因子計算結(jié)果與有限元解吻合更好，誤差在10%以內(nèi)。

平板錐形；焊接膜盒；理論模型；應(yīng)力強度因子；圍線積分

0 引言

目前我國運載火箭的蓄壓器和保險閥多使用平板錐形金屬焊接膜盒。蓄壓器中的膜盒（圖1）是抑制火箭POGO振動的關(guān)鍵部件[1]，保險閥中的膜盒（圖2）是感受壓力的敏感元件[2]，膜盒在工作過程中要經(jīng)歷反復(fù)充壓和振動環(huán)境，因此必須具備足夠的強度和疲勞壽命。閥門和蓄壓器等產(chǎn)品的膜盒是由多個膜片利用精密焊接技術(shù)沿其內(nèi)外邊緣焊接而成的管狀元件，每件膜盒的焊縫數(shù)量從幾十到幾百條不等，其結(jié)構(gòu)特點決定了焊縫應(yīng)力具有典型的集中分布特征[2,3]，在振動載荷作用下，疲勞裂紋擴展是典型的失效模式。目前有關(guān)焊接膜盒的研究主要采用傳統(tǒng)分析方法研究其強度、剛度、穩(wěn)定性和疲勞壽命。樊大鈞[4]基于旋轉(zhuǎn)殼方程，推導(dǎo)了焊接錐形波紋管（焊接膜盒）在受集中力和均布壓力的情況下，其殼體的基本微分方程。王亞軍等[5]針對平板錐形金屬膜盒，得出適用于拉伸、壓縮軸向剛度不一致的膜盒抗彎剛度計算方法。劉錦凡[6]指出蓄壓器金屬膜盒具有一定的機械剛度，其會影響蓄壓器對POGO的抑制效果。余峰等[2]對保險閥膜盒疲勞斷裂問題進行了機理分析，表明膜盒焊縫存在明顯的應(yīng)力集中，容易發(fā)生疲勞破壞。這些傳統(tǒng)分析方法未能充分考慮焊接膜盒應(yīng)力集中的特點，而膜盒的破壞是典型的裂紋擴展導(dǎo)致的破壞，為更好評估焊縫應(yīng)力集中的程度，需要計算焊縫處的應(yīng)力強度因子。國內(nèi)外學(xué)者在應(yīng)力強度因子的求解方面已有較多研究，以前的應(yīng)力強度因子求解主要有圍線積分[7,8]、擴展有限元[9,10]、虛擬裂紋閉合技術(shù)[11]以及理論求解。焊接膜盒作為一類典型的軸對稱結(jié)構(gòu)，在受軸對稱載荷下的應(yīng)力強度因子理論分析具有較普遍的應(yīng)用價值，國內(nèi)張淳源[12]針對雙梁模型開展了理論推導(dǎo)，除此以外未看到對類似結(jié)構(gòu)的理論求解。為方便工程中快速計算平板錐形膜盒的應(yīng)力強度因子，滿足工程設(shè)計，本文建立了膜盒的梁模型和薄板模型，推導(dǎo)了平板錐形膜盒應(yīng)力強度因子計算公式，并將理論計算結(jié)果與ABAQUS仿真結(jié)果進行了對比和分析。

圖1 蓄壓器的膜盒結(jié)構(gòu)

圖2 保險閥結(jié)構(gòu)

1 應(yīng)力強度因子理論計算

1.1 應(yīng)力強度因子與外力功的關(guān)系

圖3 平板錐形膜盒膜片結(jié)構(gòu)圖

表1 膜片參數(shù)尺寸

根據(jù)線彈性斷裂力學(xué)理論，I型裂紋的應(yīng)力強度因子與能量釋放率的關(guān)系為

式中代表I型裂紋的能量釋放率，為I型裂紋的應(yīng)力強度因子，對于平面應(yīng)力來說，對于平面應(yīng)變，為材料的彈性模量，為泊松比。取圖3中的一對膜片，并沿膜片周向截取單位長度，如圖4所示，其受到軸向集中力和均布載荷的作用，對于絕熱系統(tǒng)，能量釋放率的計算可以根據(jù)系統(tǒng)的能量變化關(guān)系確定，如式（2）所示。

目前僅考慮線彈性，因此不存在塑性功。而裂紋擴展單位面積釋放的能量為

對于圖4中的單個膜片，在彈性范圍內(nèi)，外力功與應(yīng)變能為

將（7）、（8）代入式（3）和（4），系統(tǒng)能量釋放率可以表示為

式中為裂紋長度，為厚度（本文取單位厚度=1），將式（9）代入式（1）即可求出應(yīng)力強度因子為

1.2 梁模型的應(yīng)力強度因子

膜盒膜片是軸對稱旋轉(zhuǎn)殼體，如圖3所示，為求解其應(yīng)力強度因子，可取用半個波，由于波高比內(nèi)外徑小得多，可取周向截取寬度單位為1的小段，近似為梁[13]，在軸向力的作用下，其力與位移的關(guān)系如圖5所示。

圖5 梁模型

圖5中，C點力矩平衡有

梁的應(yīng)變能為

將式（14）代入（13），結(jié)合式（6），外力功為

將式（15）代入式（10）可得應(yīng)力強度因子為

裂紋長度（a=13.9mm）下，應(yīng)力強度因子隨軸向力P和分布力q變化的曲線如圖6所示。

1.3 薄板模型的應(yīng)力強度因子

圖7 膜片受力圖

邊界條件如式（23）所示

式中

將式（26）代入式（10），求出平面應(yīng)變下應(yīng)力強度因子表達(dá)式，如下式所示

應(yīng)力強度因子隨軸向力和分布力q變化的曲線如圖8所示。

2 應(yīng)力強度因子有限元計算

對于應(yīng)力強度因子的求解，ABAQUS軟件提供了圍線積分法、擴展有限元法、虛擬裂紋閉合技術(shù)三種求解方法，對于靜態(tài)裂紋，使用圍線積分的方式較為簡單，設(shè)置材料屬性時不需要對破壞準(zhǔn)則與裂紋擴展條件相關(guān)參數(shù)進行定義。求解過程如圖9所示。

表2 應(yīng)力強度因子計算結(jié)果

圖9 有限元求解步驟

圖10 膜片結(jié)構(gòu)及其剖分圖

圖11 裂紋尖端奇異單元

圖12 應(yīng)力應(yīng)變云圖

表3 應(yīng)力強度因子計算結(jié)果

圖13 應(yīng)力強度因子變化曲線

3 計算結(jié)果分析

將梁模型、薄板模型的求解結(jié)果與基于圍線積分的有限元解進行比較，如圖14所示，可以看出基于梁模型與薄板模型求解的焊接膜盒裂紋尖端的應(yīng)力強度因子與有限元結(jié)果比較接近。驗證了理論推導(dǎo)公式的相對正確性。

圖14 應(yīng)力強度因子變化曲線

各模型應(yīng)力強度因子與有限元解的相對誤差如表4所示。表中的數(shù)據(jù)為給定軸向力時，不同內(nèi)壓下的平均值，從表中可以看出，梁模型的應(yīng)力強度因子平均誤差約為17.46%，而薄板模型的平均誤差為9.32%，因此薄板模型的應(yīng)力強度因子公式計算結(jié)果更接近于有限元解。

表4 平均相對誤差

4 結(jié)論

1）本文基于梁理論、薄板理論分別建立了平板錐形焊接膜盒的梁模型、薄板模型，并推導(dǎo)出膜片焊菇處的應(yīng)力強度因子表達(dá)式?；贏BAQUS的圍線積分方法，計算了軸向集中載荷以及內(nèi)壓作用下焊菇處裂紋尖端的應(yīng)力強度因子。

2）通過對比應(yīng)力強度因子的理論計算解與有限元解，梁模型的平均相對誤差在17%左右，薄板模型平均相對誤差在9%左右，因此可使用薄板模型開展膜盒焊菇處裂紋尖端的應(yīng)力強度因子的工程計算。

[1] 張智, 王楠, 劉竹生.中國載人運載火箭POGO抑制技術(shù)研究[J].中國科學(xué):技術(shù)科學(xué),2014,44:504-509.[Zhang Zhi, Wang Nan, Liu Zhusheng. POGO reduction technology of Chinese manned launch vehicles[J]. Sci Sin Tech, 2014, 44:504-509.]

[2] 余峰, 武園浩, 孫法國, 等.保險閥膜盒疲勞斷裂分析及結(jié)構(gòu)改進[J].強度與環(huán)境, 2019, 46(2): 42-47. [Yu Feng, Wu Yuanhao, Sun Faguo, et al. Fracture analysis and structure improvement of welded metal bellows of safety valve[J].Structure & environment engineering, 2019,46(2):42-47.]

[3] 王升林, 馮震, 劉士杰, 等.焊接結(jié)構(gòu)對某型金屬波紋管使用壽命的影響分析[J].強度與環(huán)境, 2021, 48(1):32-39.[Wang Shenglin, Feng Zhen, Liu Shijie. et al. Analysis of the influence of welding structure on the service life of a certain metal bellows[J].Structure & Environment Engineering, 2021, 48(1):32-39.]

[4] 樊大鈞, 著.波紋管設(shè)計學(xué)[M].北京:北京理工大學(xué)出版社, 1988: 286-289.

[5] 王亞軍, 陳牧野, 周浩洋. 平板錐形金屬膜盒內(nèi)壓柱失穩(wěn)理論研究[J].導(dǎo)彈與航天運載技術(shù), 2019(6):11-16.[Wang Yajun, Chen Muye, Zhou Haoyang. Research on column buckling calculation formula under internal pressure of plane-cone shaped bellows[J]. Missiles and Space Vehicles, 2019(6):11-16.]

[6] 劉錦凡, 孫丹, 陳雪巍, 等.蓄壓器膜盒機械剛度對液體火箭POGO振動影響研究[J]. 振動與沖擊, 2016, 35(19): 168-171.[Liu Jinfan, Sun Dan, Chen Xuewei, et al. Influences of mechanical stiffness of accumulator on POGO vibration of liquid rockets [J].Journal of Vibration and Shock,2016, 35(19): 168-171.]

[7] Rice J R. A path-independent integral and the approxi-mate analysis of strain concentration by notches and cracks[J]. ASME Journal of Applied Mechanics, 1968, 35(2): 379-386.

[8] 楊曉翔, 范家齊, 匡震邦. 求解混合型裂紋應(yīng)力強度因子的圍線積分法[J]. 計算結(jié)構(gòu)力學(xué)及其應(yīng)用, 1996, 13(1):84-89. [Yang Xiaoxiang, Fan Jiaqi, Kuang Zhenbang.A contour integral method for stress intensity factors of mixed-mode crack[J]. Computational Structural Mech-anics and Applications, 1996, 13(1): 84-89. ]

[9] Huang R, Sukumar N, Prevost J H. Modeling quasi-static crack growth with the extended finite element method. part ii:numerical applications[J]. International Journal of Solids and Structures, 2003, 40(26): 7539-7552.

[10] 郭歷倫, 陳忠富, 羅景潤, 等. 擴展有限元法及應(yīng)用綜述[J].力學(xué)季刊, 2011, 32(4): 612-625. [Guo Lilun, Chen Zhongfu, Luo Jingrun, et al. A review of the extended finite element method and its applications[J]. Chinese Quarterly of Mechanics, 2011, 32(4): 612-625.]

[11] Raju I D. Calculation of strain-energy release rates with high-order and singular finite elements[J]. Engineering Fracture Mechanics, 1987, 28: 251-274.

[12] 張淳源. 雙梁、雙板模型裂紋的能量釋放率[J]. 固體力學(xué)學(xué)報,1980(1): 117-122.[Zhang Chunyuan. Energy release rate of cracks in double beam and double plate model[J]. Chinese Journal of Solid Mechanics, 1980(1): 117-122.]

[13] 李添祥, 羅小平, 郭志, 等.聚四氟乙烯膨脹節(jié)軸向剛度計算與性能研究[J]. 華南理工大學(xué)學(xué)報(自然科學(xué)版), 1999,27(3):3-5.[Li Tianxiang, Luo Xiaoping, Guo Zhi, et al. Research on axial spring rate calculation and behavior of polytetrafluoroethylene expansion joint[J]. Journal of South China University of Technology (Natural Science Edition), 1999, 27 (3): 3-5.]

[14] 陸明萬, 羅學(xué)富, 編著. 彈性理論基礎(chǔ)[M]. 北京: 清華大學(xué)出版社, 2001: 143-147.

Calculation of Stress Intensity Factors for Plane-cone Shaped Bellows

DONG Peng-peng HE Qi-lin WANG Hai-zhou WANG Ru-wen

(Beijing Institute of Astronautical Systems Engineering, Beijing Key Laboratory of Cryogenic Technology Research, Beijing 100076, China )

Based on beam and thin plate theory, the theoretical model of plane-cone shaped welded bellows is established, and the formula of the stress intensity factor of the crack tip is derived for the axial load and pressure condition. The stress intensity factor is also calculated by the use of contour integral method of ABAQUS code to evaluate the accuracy of theoretical formula. The results show that the thin plate model coincides better with the finite element method with the error less than 10%.

Plane-cone shaped; Welded bellows; Theoretical model; Stress intensity factor; Contour integral

O346.1, V421.3+3

A

1006-3919(2021)04-0030-07

10.19447/j.cnki.11-1773/v.2021.04.006

2021-03-17；

2021-05-20

董鵬鵬(1995—)，男，碩士研究生，研究方向：焊接膜盒的疲勞與斷裂行為研究；（100076）北京9200信箱10分箱12號.