板壳式换热器有限元强度计算.docVIP

????特殊板壳式换热器封板的有限元模拟及强度研究????????????????????????????????????孙爱芳1,刘敏珊2,董其伍2(1.河南机电高等专科学校,河南新乡453002;2.郑州大学热能研究中心,河南郑州450002)????摘要:对于板壳式换热器中的特殊开口封板,利用APDL语言(ANSYSparametricdesignlanguage)实现全参数?化驱动的三维有限元模型,通过三维有限元模拟,并对关键区域作路径分析,获得封板的应力分布规律.实验结果?验证了数值模拟的有效性.经过分析,归纳出封板中的应力影响因素,提出工程关联式模型.以不开口区环板上的?最大综合应力为研究对象,对封板进行参数化数值分析,并对数值模拟结果进行数据处理与拟合,得到封板强度计?算工程算式.拟合公式所得数据与原模拟数据相比,平均绝对误差为0·238MPa,平均相对误差为1·8%,这对于工?程应用已经足够精确,可为新产品的工程应用提供参考.????关键词:板壳式换热器;封板;有限元模拟;强度????中图分类号:TQ051·5　文献标识码:A　文章编号:1006-754X(2008)04-0295-05????针对实际生产中众多设备的腐蚀失效,作者研?发了导热复合材料紧凑型板壳式换热器[1],其中封?板采用“薄管板”的贴面式结构[2],是连接壳体、芯体?和封头,并承受压力和热膨胀以及来自此三部件载?荷的主要部件.与管壳式换热器的管板所处的位置?和作用相似,封板是板壳式换热器中举足轻重的关?键元件,但有着与管板开孔形式迥然不同的开口方式(如图1所示),其与强度计算有关的问题虽然可?以参考管板的解决思路,但毕竟是一个全新的课题.?为了开拓其应用范围,保证设备安全运行,需要对封?板进行深入的研究[3].????用解析法求解封板的应力是非常繁琐的,难以?提出一个普遍适用的公式用以计算各种情况下的危?险应力[4].通过建立与实际设备相似的模型,无疑会大大提高计算结果的正确度[5].应用计算机辅助设?计技术、三维有限元结构分析技术和三维温度场分?析技术进行封板结构力学分析[5-7],也是换热器分析?与设计的总的发展趋势.本文利用参数化三维有限?元模拟整个封板连接系统,获得封板强度设计的工?程关联式,为新型换热器设计提供理论和工程应用?的依据.?????????1　封板的三维有限元分析????根据换热器和封板几何结构的特点,采用实体?单元和壳单元相结合模拟,封板和板束材料为短切?碳纤维增强聚四氟乙烯复合材料,研究表明可视为?各向同性的弹性材料[1],模型的建立以及分析过程?采用APDL语言实现完全参数化驱动.????1.1　温度场分析????封板上下表面分别为板程和壳程温度,存在温?差,加之复合材料本身导热性不良,需要考察其温度?应力的分布情况.本文分析重点在于换热器稳定工?作情况下封板的温度分布情况,因此采用稳态热分?析数学模型.????图2为封板连接系统的温度场.其温度载荷在?计算后存储于相应的RTH文件中,结构分析时直?接读入即可.????在封板中心芯体连接区,靠近低温侧管板出现?了较大的温度梯度,这里结构尺寸变化大,热应力的?影响可能比较大.取封板中心区厚度路径上的温度?分布如图3所示,不存在温度急剧变化的区域.可以?预测,不会产生较大的热应力.????1.2　应力场分析????主要考虑封板的应力分布规律,忽略封头等,考?虑到结构和载荷的对称性,沿换热器的纵向对称面?切开取其四分之一作为分析模型体.结构纵向对称面约束了法向位移,壳程筒体横截面约束了轴向位?移.封板两侧分别施加板程和壳程压力.网格剖分如?图4所示.?????????????????为使封板与隔板的胶结接头不致承受大的剥离?应力,换热器开启时要求先开板程,封板强度按只承?受板程的设计压力来进行设计计算.应力强度云图?如图5所示,可见封板上最大应力强度可能发生在?板束外缘处,为此设定四条路径来分析封板的应力.?四条路径上的应力分布见图6和图7.????可见,在机械载荷作用下,封板的应力云图出现?有规律的变化.封板径向(路径1和路径2)上,开口区应力由中心向边缘缓慢地逐渐增大,不开口区应?力有波动,并存在较大的应力;在封板与壳程筒体相?连接的地方出现了较大的应力集中,这与经典理论?相符.????在路径1和路径2中不开口区显示的应力最大?值处的周向上(路径3和路径4),也存在较大应力?区,位于开口区的拐角处,这是应力集中与无板束支承双重原因造成的.?通过路经分析,了解了封板的整体应力分布情?况,其中路径4上既避开了明显的应力集中,又存在?最大应力,后续的分析均以此为研究对象.????1.3　实验研究????作者设计了一台试验装置