吸收塔塔体封头应力分析.docxVIP

第二章吸收塔塔底封头应力分析 2.1吸收塔结构与工艺 吸收塔的工艺，操作工况，主要结构陈述一下，把图纸放上去。 2.2有限元建模与计算2.2.1模型的建立 由于塔底封头处的应力状态为本文研究的重点，因此，为简化计算，不考虑塔体上接管等局部结构对塔底应力状态的影响，将整个塔体作为轴对称结构处理，建立二维轴对称模型，如图2-4所示。其中，将塔板与格栅组合成的复杂结构等效为刚度和强度都较大（设置较高的弹性模量和屈服强度）的钢板置于塔内，将填料、吸附剂等塔内内件及工艺介质对各层钢板的作用力等效为均布压力P，塔体及钢板的自重设为G。  2.2.2网格划分塔体的网格划分如图2-5所示，采用CAX4R或CAX3单元，共划分21598个节点，18647个单元，其中局部可能应力集中的区域采用较密集的网格。 2.2.3材料参数 冷却塔塔体材料为16MnR钢，其性能参数见表2-1,此外，塔板与格栅的等效钢板的弹性模量E取为600GPa，屈服强度八取为600MPa。 表2-1塔体材料16MnR刚的性能参数 材料 温度T °C 密度p (g/mm3) 弹性模量E (GPa) 屈服强度。 y (MPa) 16MnR 200 7.86 200 319 2.2.3载荷与边界条件 根据工厂提供的吸附塔操作数据可知，增加塔的产量后，塔体内部填料、吸附剂及介质等的总重量约为2550t，将此重量均匀分布给各层塔板与格栅等效钢板，可得各层钢板承受的质量为2550/13=196t，等效钢板的面积为;[30002-(425+30)2]=2.76x107mm2，因此每层等效钢板所承受的均布压力户=0.07MPa。另外，取重力加速度g=9.8m/s2。 设置图2-4所示模型的左边界为轴对称边界，限制裙坐支撑的塔体部分的各个方向移动与转动。 2.3结果与讨论 2.3.1整体应力分布 利用有限元ABAQUS软件，根据上述模型及载荷与边界条件，计算得到塔体的应力分布。图2-6给出了塔体的整体应力分布图，由图2-6可见，塔体在底部封头支撑处存在明显的应力集中，应力最大位置位于封头外表面支撑所在位置的两侧，最大值达到319MPa，超过了塔体封头材料16MnR钢在200°C时的屈服强度。 图2-6塔体整体应力分布 2.3.2塔体底部封头局部应力分布 为详细考察塔底封头局部区域的应力分布情况，取图2-7中所示的路径P1,图2-8?2-12分别给出了此局部区域的Mises应力、最大主应力、径向应力S11、轴向应力S22、周向应力S33的分布云图及沿路径P1的应力分布规律。取路径P1即是为了更直观的表示此局域的应力大小与分布情况。 由图2-8可见，封头在此支撑对应位置的外表面存在明显应力集中，有较大区域的应力超过了材料的屈服强度，即发生了屈服。而由图2-9可知，封头在此局部区域的最大主应力甚至达到了365MPa，超过了16MnR钢的屈服强度。 P1图2-7路径P1(a)Mises应力分布云图 0100200300400500600700800距离d/(mm0100200300400500600700800 距离d/(mm) 图2-8塔体封头局部区域的Mises应力分布 SfMax.Principal(Avg:75%)365.003326.913288.823250.733212.643174.553136.46398.37360.28322.193-15.897-53.987-92.07705050504332211X匕X力应主大最L(a)最大主应力分布云图00100200300400500600700800距离d/(mm) (b)沿路径P1的最大应力分布 图2-9塔体封头局部区域的最大主应力分布S,Sil(Avg:75%)357.678301.629245.581~k189.532~k133.483_k77.434M-21.385_k-34.664—k-90.712M--146.761-202.810-258.859■■-314.908(a)径向应力S11分布云图 400350300250200150100500100200300400500600700800距离d/(mm) 400 (b)沿路径P1的径向应力分布 图2-10塔体封头局部区域的径向应力分布 (a)轴向应力S11分布云图S,S22(Avg:75%)221.236179.525.-137.813——96.102——54.390——12.679--29.032■--70.744112.455.--154.166--195.878 (a)轴向应力S11分布云图 -237.589■■-279.301 应向轴 (b)沿路径P1的轴向应力分布 图2-11塔体封头局部区域的轴向应力分布 (a)周向