波紋管的極限載荷分析

李志平　段成紅

波紋管的極限載荷分析

李志平*段成紅

（北京化工大學機電工程學院）

采用極限載荷分析法對單層和多層波紋管進行極限載荷分析，不僅研究了波紋管拉伸位移及相關參數(shù)對波紋管極限承載能力的影響，也研究了同等厚度單層波紋管和多層波紋管極限載荷之間的關系。結果表明：波紋管極限載荷不隨拉伸位移和層間接觸狀態(tài)的變化而變化，但隨著波紋管波形參數(shù)增大而減??；多層波紋管極限載荷小于同等厚度的單層波紋管。

波紋管極限載荷層間接觸有限元模型摩擦

0　前言

波紋管作為補償管線，也作為容器熱變形以及吸震、減震的柔性元件，已經越來越廣泛應用于各種工程實踐中。關于波紋管力學性能的研究前人做了大量工作，但大都采用理論方法［1-3］或有限元應力分析法進行研究［4-6］。目前雖有波紋管極限載荷分析方面的研究，但其分析主要集中在波紋管極限承載能力方面，并未考慮相關因素對波紋管極限載荷的影響［7］。不僅如此，工程上多層波紋管仍然按照單層設計，未考慮層間摩擦和層間間隙問題［8-9］。為精確計算波紋管極限載荷，研究波紋管拉伸位移、結構參數(shù)、層間摩擦和層間間隙對波紋管極限載荷的影響，本文運用極限載荷分析法對U形單層、多層波紋管進行極限載荷分析，其結果對工程實踐具有一定的參考價值。

1　波紋管結構參數(shù)和材料參數(shù)

1.1波紋管結構及其參數(shù)

波紋管是母線呈波紋形的回轉薄壁殼體，有多種形式，如U形、V形和Ω形等。本文所研究的U形波紋管具體結構如圖1所示，結構參數(shù)取自某工程實際膨脹節(jié)波紋部分，具體如表1所示。

圖1　波紋管結構

表1　單（多）層波紋管結構尺寸

1.2波紋管材料參數(shù)

本文研究的U形波紋管材料為S30408不銹鋼，材料參數(shù)如表2所示。

2　有限元模型

2.1單層波紋管有限元模型

由于單層波紋管厚度方向尺寸遠小于其他方向尺寸，因此可以當做殼體考慮。又由于分析涉及非線性，為減少計算成本，提高計算效率，本文采用支持塑性分析的殼單元Shell 181建立單層雙波波紋管全模型。

單層模型厚度為實際成形厚度tp，計算公式如式（1）所示：

式中t——單層名義厚度，mm；

Db——波紋管內徑，mm；

Dm——波紋管波紋的平均直徑，mm；

w——波高，mm；

n——層數(shù)。

單層波紋管有限元模型采用規(guī)則四邊形進行網(wǎng)格劃分，網(wǎng)格尺寸2×2 mm2，如圖2所示。

圖2　單層波紋管有限元模型

2.2多層波紋管有限元模型

考慮到極限載荷分析計及幾何非線性問題以及多層波紋管層與層之間的相互影響，為減少計算成本、提高計算效率，本文采用支持塑性及接觸計算的平面單元Plane 182，建立三層雙波波紋管軸對稱模型。波紋管尺寸如表1所示，模型厚度為實際成形厚度tp，計算公式如式（1）所示。

三層波紋管有限元模型采用規(guī)則四邊形進行網(wǎng)格劃分，厚度方向4份。為保證計算精度，網(wǎng)格尺寸長寬比一般不超過2∶1，故模型經向方向網(wǎng)格尺寸為0.5 mm。

根據(jù)接觸算法特點及波紋管結構特性，經過多次試算，接觸計算所用接觸剛度FKN取0.5；侵入容差FTOL取0.1，且采用內層波紋管外表面為目標面、外層波紋管內表面為接觸面的柔-柔接觸。三層波紋管有限元模型如圖3所示。

圖3　多層（3層）波紋管有限元模型

3　單層波紋管極限載荷分析

3.1拉伸位移對波紋管極限載荷的影響

為研究拉伸位移對波紋管極限載荷的影響，本節(jié)計算波紋管一端固定，一端分別施加5 mm、10 mm、15 mm、20 mm、25 mm和30 mm拉伸位移時的極限載荷。

由圖4可知，在材料為理想彈塑性的情況下，波紋管拉伸時極限載荷比兩端固定（圖中位移為0）時的極限載荷低，但隨拉伸位移的增大，波紋管極限載荷并沒有下降，而是基本保持不變。

圖4　波紋管極限載荷隨拉伸位移變化曲線

3.2波形參數(shù)對波紋管極限載荷的影響

波形參數(shù)波高w和波距q是波紋管的關鍵參數(shù)，波高對波紋管經向彎曲應力及剛度有非常大的影響。波高w和波距q存在內在聯(lián)系，通常薄壁波紋管波高w與波距q的比值w/q為0.8～1.2。

（1）波高對波紋管極限載荷的影響

本節(jié)研究波距q為64 mm情況下，w/q處于0.8～1.2之間時波高對波紋管極限載荷的影響。根據(jù)q＝64 mm、w/q＝0.8～1.2，分別選取w為52 mm、56 mm、60 mm、64 mm、68 mm、72 mm和76 mm進行極限載荷分析。計算時波紋管兩端分別施加固定約束。

由圖5可知，在材料為理想彈塑性的情況下，波紋管極限載荷隨著波高的增大而明顯下降。因此，在波距確定和w/q＝0.8～1.2的情況下，應盡量選擇小波高波紋管。

（2）波距對波紋管極限載荷的影響

本節(jié)研究波高w為52 mm情況下，w/q處于0.8～1.2之間時波距對波紋管極限載荷的影響。根據(jù)w＝52 mm、w/q＝0.8～1.2，分別選取q為44 mm、48 mm、52 mm、56 mm、60 mm和64 mm進行極限載荷分析。計算時波紋管兩端分別施加固定約束。

由圖6可知，在材料為理想彈塑性的情況下，波紋管極限載荷隨著波距的增大而明顯下降。因此，在波高確定和w/q＝0.8～1.2的情況下，應盡量選擇小波距波紋管。

圖5　波紋管極限載荷隨波高變化曲線

圖6　波紋管極限載荷隨波距變化曲線

4　多層波紋管極限載荷分析

4.1摩擦系數(shù)對波紋管極限載荷的影響

對于不銹鋼材料，層間靜摩擦系數(shù)為0.15，動摩擦系數(shù)為0.1［10］。由于實際應用中波紋管層與層之間處于動靜摩擦之間的某個狀態(tài)，故摩擦系數(shù)在0.1～0.15之間波動，本節(jié)研究摩擦系數(shù)為0.1、0.11、0.12、0.13、0.14和0.15時波紋管的極限載荷。計算時波紋管兩端分別施加固定約束。

由圖7可知，在材料為理想彈塑性的情況下，波紋管極限載荷隨波紋管層間摩擦系數(shù)變化很?。ㄗ畲笞兓考s3%），因此可以認為波紋管層間摩擦系數(shù)對波紋管極限載荷幾乎沒有影響。

4.2層間間隙對波紋管極限載荷的影響

多層波紋管是由多層薄圓筒套合沖壓而成，加工時層間間隙不可能消除，一般間隙為0.1～0.3 mm［11］。由本文4.1節(jié)可知，摩擦系數(shù)對波紋管極限載荷幾乎沒有影響，故本節(jié)研究摩擦系數(shù)為0.1時，層間間隙分別為0.1 mm、0.15 mm、0.2 mm、0.25 mm和0.3 mm時波紋管的極限載荷。

圖7　波紋管極限載荷隨摩擦系數(shù)變化曲線

由圖8和圖9可知，在材料為理想彈塑性的情況下，波紋管極限載荷幾乎不隨層間間隙變化（最大變化量約1%），所以可以認為波紋管層間間隙對波紋管極限載荷幾乎沒有影響。不過需要注意的是，不同層間間隙下波紋管剛度會發(fā)生變化?？紤]到波紋管間隙中存在的氣體在溫度升高時會膨脹而導致層與層間分離、膨脹，從而影響其使用，因此，應盡可能減小波紋管層間間隙。

圖8　不同層間間隙波紋管載荷位移曲線

5　單層、多層波紋管極限載荷對比

雖然上文分別討論了拉伸位移、結構參數(shù)以及層間接觸參數(shù)對波紋管極限載荷的影響，但并未進行同等厚度的單層、多層波紋管極限載荷的比較，故本節(jié)對此進行研究。采用單層3 mm波紋管和每層1 mm的三層波紋管進行計算和對比分析，波紋管兩端分別固定約束。圖10為單層和三層波紋管載荷位移曲線。

圖9　波紋管極限載荷隨層間間隙變化曲線

圖10　單層、三層波紋管載荷位移曲線

由圖10可知，在材料為理想彈塑性的情況下，當載荷與約束情況相同，單層波紋管極限載荷為1.58 MPa，三層波紋管極限載荷為0.99 MPa，三層波紋管極限載荷明顯小于總厚度相同的單層波紋管。經分析，這是因為在總厚度相同情況下，相對于單層波紋管而言，三層波紋管每層厚度較小而導致在較小載荷下單層就發(fā)生屈服，從而使整體極限載荷下降。

6　結論

本文運用極限載荷分析法對波紋管進行極限載荷分析，研究波紋管結構參數(shù)及層間摩擦和層間間隙對波紋管極限載荷的影響，并在此基礎上對比總厚度相同的單層和三層波紋管極限載荷。在材料為理想彈塑性時主要有以下結論：

（1）在波紋管許用拉伸范圍內，一端固定一端拉伸時的極限載荷比兩端固定時的極限載荷低；但隨著拉伸位移的增大，其極限載荷并沒有下降，而是基本保持不變。由于實際應用中，波紋管兩端不可能完全固定，所以實際極限載荷應較兩端固定時的極限載荷低。

（2）波紋管極限載荷隨著波紋管波形參數(shù)波高和波距的增大明顯降低，故實際應用中，應盡量選用滿足w/q＝0.8～1.2的較小波高和波距。

（3）多層波紋管層間摩擦系數(shù)和層間間隙對波紋管極限載荷影響很小，但為防止層間間隙氣體在高溫下膨脹而導致層與層間分離、膨脹，影響波紋管使用，應盡可能減小波紋管層間間隙。

（4）同等厚度情況下，多層波紋管極限載荷較單層波紋管小。

［1］錢偉長.半圓弧波紋管的計算——細環(huán)殼理論的應用［J］.清華大學學報，1979，19（1）：84-99.

［2］朱衛(wèi)平，黃黔.中細柔性圓環(huán)殼整體屈曲的一般解及在波紋管計算中的應用（Ⅰ）——基本方程和一般解［J］.應用數(shù)學和力學，2002，23（8）：790-797.

［3］朱衛(wèi)平.U型波紋管及相關結構環(huán)向屈曲的有限元分析（Ⅰ）——基本方程及環(huán)板的屈曲［J］.應用力學學報，2002，19（4）：19-25.

［4］Stephen O Pierce，Jeffrey L Evans.Failure analysis of a metal bellows flexible hose subjected to multiple pressure cycles［J］.Engineering Failure Analysis，2012，22：11-20.

［5］楊建良，楊湖.應用TEMA-9th進行膨脹節(jié)分析設計［J］.化工設備與管道，2015，52（3）：24-27.

［6］Charles Becht.An evaluation of EJMA stress calculations for unreinforced bellows［J］.Journal of Pressure Vessel Technology，2001，124（1）：124-129.

［7］錢才富，陳朝暉，姚琳.Ω型膨脹節(jié)在內壓作用下極限承載能力分析［C］.第二屆中國CAE工程分析技術年會暨2006全國計算機輔助工程（CAE）技術與應用高級研討會論文集.青島，2006：262-268.

［8］StandardsofTheExpansionJointManufactures Association［S］.9th Edition.2008.

［9］GB 16749—1997壓力容器波形膨脹節(jié)［S］.

［10］成大先.機械設計手冊（第一卷：摩擦因數(shù)）［M］.第5版.北京：化學工業(yè)出版社，2010.

［11］李杰，段玫.多層波紋管接觸分析及穩(wěn)定性屈曲分析［J］.材料開發(fā)與應用，2011，26（6）：53-57.

Limit-Load Analysis for Bellows

Li Zhiping Duan Chenghong

The limit-load analysis method was applied to study the limit load of monolayer and multilayer bellows in this paper.Not only the influence of tensile displacement and related parameters on the ultimate load carrying capacity of bellows was studied，but also the relationship between the limit load of single-layer bellows and multilayer bellows was studied.The analysis results showed that the limit load of bellows did not change along with the change of tensile displacement and interlayer contact state，but it decreased as the bellows waveform parameters increased；the multilayer bellows'limit load was less than single-layer bellows in the same thickness.

Bellows；Limit-load；Interlayer contact；Finite element model；Friction

TQ 050.3

10.16759/j.cnki.issn.1007-7251.2016.10.009

2016-03-15）

*李志平，男，1988年生，碩士研究生。北京市，100029。