对一个管壳式换热器进行温度和压力载荷作用下的有限元强度分析,并对这两种载荷作用 下的结构应力响应做分类研究,然后对结构进行改进,作同样的分析。比较这些结果,得出考虑温度载荷作用下的换热器强度校核的规律和结构设计的特点。

管壳式换热器是化工、石油、轻工、能源等工业应用广泛的过程设备之一，它具有选材范围广，换热表面清洗较方便，适用性较强，处理能力大，能承受高温和高压等特点。管壳式换热器的结构设计主要依据是GB151［1］，GB151 中关于换热器管板强度校核是根据弹性基础上薄板理论，在轴对称结构的条件下，将薄板的三维变形简化为二维梁式变形，由此来计算其强度的。而换热器壳体厚度的选择，主要是根据壳体所受到的壳程压力来确定。

换热器由于其工作特点，不仅有管程压力和壳程压力等载荷作用，而且还要受到工作介质的温度载荷作用。在GB151 中对压力载荷，给出了管板和壳体的尺寸选择，及固定管板兼作法兰的管板和壳体的连接方式。然而，对于在温度载荷作用下，这些尺寸却没有具体的说明要求。

本文通过一个管壳式换热器的强度校核，将载荷分类为压力载荷和温度载荷，来说明结构在这些载荷作用下的应力响应特点，进而提出该结构改进的意见。本文采用三维有限元的分析方法，来研究其内在规律。

换热器结构尺寸及载荷工况

BEM型换热器结构如图1 所示，管板上共有500 根换热管，分布在管板的上半部分，左右对称。结构尺寸和材料：管板内径：1300 mm；管板厚度：80mm；法兰外径：1460 mm；管板材料：00Cr19Ni10；壳体厚度：24mm。

为简化计算，在建立有限元模型时，只考虑换热器的管板、壳体、管束和膨胀节等主要结构，法兰垫片用等效的均布比压来代替。

由于整台换热器结构是前后左右对称，所以只取组合体的四分之一，换热管长度取一半。换热管是细长形状，所以用杆单元来模拟。这样的作法对管子附近的管板应力计算是不准确的，事实上，如果考虑换热管和管板的胀焊连接，该处的真实应力也很难计算。根据圣维南原理，这样处理对远处的非布管区管板的计算影响不大。

壳体和膨胀节用20节点六面体单元划分网格，管板用10 节点四面体单元来划分网格，在管板和壳体过渡的区域是13节点的五面体锥型单元。管束用2节点杆单元来划分网格。表示换热管的杆单元节点与表示管板的实体单元节点在对应位置上重合。整个分析过程?**用ANSYS有限元软件来完成。

在结构的前后对称面和左右对称面上加上对称边界条件，即这些面上的法向位移为零。换热管的一端固定位移，另一端与管板连接。此外，还需要限制整体结构的刚体位移。

取正常操作工况为校核的工况，其具体数值为：壳程压力：Ps .= -0.1 MPa；壳程温度：Ts=230℃；管程压力：Pt =0.2MPa；管程温度：Tt=50℃；管板温度：T=140℃；法兰垫片压力：Pc=69MPa；法兰螺栓预紧力：116.5kN。

由于管板上的开孔面积只占其总面积的约5%，所以忽略当量压力和管程压力之间的差别。

根据换热器操作工况下温度载荷的经验数据，其分布如下：在结构的下半部，管板的中截面处的温度为140 ℃，膨胀节处的温度为230 ℃，在结构的上半部，管板的中截面处的温度为120℃，膨胀节处的温度为200 ℃，从管板到膨胀节，根据空间位置的不同，进行双线性插值。筋板的温度，也按照这一规律进行插值。