许多商业有限元分析软件,如ＡＮＳＹＳ、ＡＢＡＱＵＳ、ＭＡＲＫ、ＮＡＳＴＲＡＮ,都已经具备了单元自动剖分的功能。但是,到目前为止,对于一般空间问题只能自动生成8节点单元的蜕化元,即四面体或五面体单元,而不能自动生成所有单元都是较规则的8节点六面体单元的有限元模型,四面体或五面体单元是常应变单元,不能较好模拟应力变化较为剧烈的应力分布状态。由于换热器所用膨胀节是弯曲应力变化非常剧烈的部位,因此,四面体或五面体单元用于换热器膨胀节的有限元应力分析效果不好。非金属膨胀节厂家广泛生产波纹膨胀节应用于各行业，许多研究者在对膨胀节进行有限元分析时,为了回避8节点单元模型建模烦琐这一问题,常常由有限元分析软件自动生成单元网格。这样剖分的结果是在膨胀节及其附近高应力变化区形成大量极不规则的四面体单元或五面体单元,而该区域正是应力分布最为复杂的区域。文中提出在对膨胀节进行有限元分析时,保证膨胀节及其附近区域内单元都是比较规整、均匀的8节点单元时才能得到较好的计算结果,这种单元可以使用不协调位移模式,能够较好地模拟膨胀节、开孔接管、弯曲梁及拱坝等受弯曲为主的工程结构。为了验证这一结论,以多乙二醇塔再沸器膨胀节为例进行了有限元分析对比的计算。

多乙二醇塔Ｅ 710再沸器是乙二醇装置的重要设备之一,此再沸器为壳体上带有两组4个膨胀节的固定管板式换热器,型号为ＢＥＭ450 2.5/0.2 6.4 6/19 3Ｉ,其管程设计压力0.2ＭＰａ,设计温度200℃;壳程设计压力2.5ＭＰａ,设计温度为250℃。换热管规格为 19ｍｍ×1.6ｍｍ×6000ｍｍ,材料为ＳＵＳ304,壳体规格为 450ｍｍ×8ｍｍ×5990ｍｍ,壳体材料为ＳＵＳ304,膨胀节按ＧＢ16749—1997《压力容器波形膨胀节》选型,其型号为ＨＺＷＣ450 2.5 1×8×2,材料为ＳＵＳ304。

有限元对比分析

分别采用比较规整的8节点六面体ｓｏｌｉｄ45单元和软件自动剖分功能生成的四面体单元或五面体单元建立膨胀节三维有限元模型,见图3。考虑到结构和载荷都是轴对称的,故取包含两个膨胀节的部分壳体的1/4作为计算模型。整个计算模型共划分为约65280个单元,沿壁厚方向取4层可以较好地模拟径向弯曲应力的剧烈变化。

(1)位移边界条件　在1/4膨胀节壳体的两个纵向轴对称侧面内,沿对称面法线方向的自由度为0;在膨胀节壳体一侧横断面内轴向自由度为0。

(2)力边界条件　在膨胀节壳体另一侧横断面内,由内压和温度载荷的共同作用产生的轴向应力,可由下述方法确定。即由圆柱壳体、膨胀节、换热管束和管板组成的管箱,在内压和温度载荷的作用下是一个静不定系统,解这个静不定系统即可得到壳体横断面内的轴向应力。在内压和温度作用下,管箱处于平衡状态,管箱内力的平衡方程为:ｐｉＡｉ=Ｆｓｈ+Ｆｔｕ,壳体与管束变形协调方程为Δｌｓｈ=Δｌｔｕ,其中,ｐｉ=2.1ＭＰａ,为壳程操作内压;Ａｉ为管板有效受压面积,ｍ2;Ｆｓｈ、Ｆｔｕ分别为壳体、管束内力,Ｎ;Δｌｓｈ、Δｌｔｕ分别为壳体、管束的轴向变形量,ｍ。解内力平衡方程与变形协调方程组成的方程组可得Ｆｓｈ=1 118×104Ｎ,故在膨胀节壳体横断面内的力边界条件为ｐｘ=Ｆｓｈ/Ａｓｈ=-0 97126ＭＰａ,其中Ａｓｈ=0 01151ｍ2,为壳体横截面积。