温度对聚丙烯环管反应器热膨胀的影响研究_王政文.pdf

2023年05月|125工艺管控套之间用换热夹套连接管相连，同时也考虑了换热夹套连接管因受温度影响而发生热膨胀的情况，因此在换热夹套连接管上也设置了换热夹套连接管膨胀节，以消除两端换热夹套因受温差作用发生的横向膨胀位移以及换热夹套连接管因受温差作用自身发生的轴向膨胀位移(工作时换热夹套连接管膨胀节处于压缩状态)1-3。1 温度对换热夹套上端膨胀节的影响研究图 1 所示为对环管反应器换热夹套上端膨胀节的应力分析结果，内容包括网格划分、施加载荷、位移约束、总变形云图、轴向变形云图以及等效应力云图。首先，在软件 ANSYS Static Structural 里的 Mesh Mechanical 下对换热夹套上端膨胀节进行结构化网格划分，网格划分总数为 42 480，并经过网格无关性0 引言本文对于环管反应器因受温度影响而发生热膨胀的研究，是基于反应器内筒体与换热夹套因受温差作用而发生膨胀(轴向与横向)的过程，但内筒体所受的温差大于换热夹套所受的温差，因此内筒体的轴向膨胀位移大于换热夹套的轴向膨胀位移，导致换热夹套对内筒体轴向膨胀上有一定的约束作用，而为了消除内筒体由于换热夹套约束而产生的较大温差应力，决定在换热夹套上设置膨胀节。而环管反应器在环形支座位置处将换热夹套与内筒体通过径向支撑板相连且固定，并将换热夹套分为上下两段，因此为了能够更好的消除内筒体上产生的较大温差应力，经过研究决定在换热夹套上下两端都设置了膨胀节(即换热夹套上端膨胀节与换热夹套下端膨胀节，工作时换热夹套上下两端膨胀节都处于拉伸状态)。两个换热夹温度对聚丙烯环管反应器热膨胀的影响研究王政文1，张万尧1，余金锭2(1.天华化工机械及自动化研究设计院有限公司 国家干燥技术及装备工程技术研究中心，甘肃 兰州 730060；2.福建中景石化有限公司，福建 福清 350309)摘要：文章针对温度对聚丙烯环管反应器热膨胀的影响进行了研究，首先通过运用三维软件、ANSYS 软件等，对聚丙烯环管反应器膨胀节进行三维模型的建立以及网格划分；其次通过运用 ANSYS 软件对膨胀节进行应力分析，通过分析计算可以得到受外力与温差作用下膨胀节膨胀情况的应力云图；最后通过对应力分析结果的对比研究，进一步验证受温度影响膨胀节膨胀的位移与筒体膨胀位移是否一致，减小不同步膨胀产生的应力约束；并且采用第三强度理论与第四强度理论分别对环管反应器温差叠加后的总应力水平进行有效评价。关键词：聚丙烯环管反应器；膨胀节；ANSYS；应力分析中图分类号：TB3文献标志码：A 文章编号：1008-4800(2023)14-0125-06DOI:10.19900/ki.ISSN1008-4800.2023.14.038Study on the Effect of Temperature on the Thermal Expansion of Polypropylene Loop ReactorWANG Zheng-wen1,ZHANG Wan-yao1,YU Jin-ding2(1.National Drying Technology and Equipment Engineering Technology Research Center of Tianhua Chemical Machinery and Automation Research and Design Institute Co.,Ltd.,Lanzhou 730060,China;2.Fujian Zhongjing Petrochemical Co.,Ltd.,Fuqing 350309,China)Abstract:In this paper,the influence of temperature on the thermal expansion of polypropylene loop reactor is studied.Firstly,through the use of three-dimensional software,ANSYS software,the three-dimensional model of the expansion joint of polypropylene loop reactor is established and the mesh is divided.Secondly,through the use of ANSYS software to analyze the stress of the expansion joint,through analysis and calculation can be obtained under the action of external force and temperature difference of the expansion of the stress cloud diagram.Finally,through the comparative study of the stress analysis results,it is further verified that the expansion displacement of the expansion joint affected by temperature is consistent with the expansion displacement of the cylinder,so as to reduce the stress constraint caused by the asynchronous expansion.The third strength theory and the fourth strength theory are used to evaluate the total stress level after the superposition of temperature differences in the ring-tube reactor.Keywords:polypropylene loop reactor;expansion joint;ANSYS;stress analysis工艺管控工艺管控126|2023年05月与应力分布情况。通过以上对应力分析前的描述与分析，决定开始对环管反应器换热夹套上端膨胀节进行位移约束与施加载荷，并对求解计算后得到的换热夹套上端膨胀节总变形云图、轴向变形云图、等效应力云图进行研究与分析。计算结果表明，换热夹套上端膨胀节在受载荷作用下的总变形最大值为 12.619 mm；轴向位移最大值为 4.598 mm，最小值为 1.622 5 mm，平均值为 3.110 25 mm；等效应力最大值为 141.7 MPa，最小值为0.729 48 MPa。通过计算结果得出的数据可知，膨胀节在内部温度 45(室外年平均气温为 20，换热夹套内部循环冷却水的温度为 45，温差则为 25)、内部压力 1.1 MPa、压力 F=6.16105 N 下的膨胀情况。而环管反应器长度为 55 000 mm 的内筒体在内部工作温度 75(换热夹套内部循环冷却水的温度为45，温差则为30)下的膨胀情况可由公式(1)得出：P1=L1 T1=55 0001.2 10-530=19.8 mm (1）将换热夹套长度也同样视为 55 000 mm，则换热夹套在内部温度 45(温差为 25)下的膨胀情况可由公式(2)得出：P2=L2 T2=55 0001.2 10-525=16.5 mm (2)通过公式(1)与(2)可知，环管反应器在正常工作时内筒体的膨胀位移为 19.8 mm、换热夹套的膨胀位移为 16.5 mm，膨胀位移差值为 3.3 mm。由于环管反应器在环形支座位置处将换热夹套分为上下两段，通过查阅图纸可知，换热夹套上段为 50 300 mm、换热夹套下段为 4 700 mm。假定将内筒体也同样分为上下两段且长度与换热夹套相同，通过计算与分析可以确定换热夹套上段与内筒体上段的理论膨胀位移分别为 15.09 mm、18.108 mm；换热夹套下段与内筒体下段的理论膨胀位移分别为 1.41 mm、1.692 mm，因此换热夹套上端膨胀节的理论膨胀位移为 3.018 mm、换热夹套下端膨胀节的理论膨胀位移为 0.282 mm。而换热夹套上端膨胀节的理论膨胀位移为 3.018 mm 与膨胀节在受载荷作用下的轴向位移平均值 3.110 25 mm 几乎相等，验证了换热夹套上端膨胀节应力分析的准确性。以上对于反应器内筒体与换热夹套因受温差作用而发生膨胀位移的分析，结论并不能对膨胀节的有效性进行整体评估，而侧重点主要还在于对膨胀节因受外力与温差作用下发生膨胀时的应力进行评估，因此为了能够更准确的对膨胀节进行应力评估，决定采用两种方法，即通过采用第三强度理论(也称“最大剪验证，且网格质量均大于 0.8。其次，为了能够更准确的模拟出内筒体与换热夹套因受温差作用发生膨胀时，换热夹套对内筒体轴向膨胀上的约束作用，决定对膨胀节下端施加位移约束，对膨胀节上端施加垂直向下的压力 N，施加内部压力 1.1 MPa。最后，在进行求解计算前，选取总变形、轴向变形、等效应力对换热夹套上端膨胀节进行应力分析。应力分析时，膨胀节下端与反应器直管段连接部位被视为固定位置，是为了给膨胀节添加位移约束；膨胀节上端施加垂直向下的压力，是为了简单还原内筒体与换热夹套在整体法兰处相连对膨胀节不受外力影响而发生自由膨胀时的约束作用所带来的约束力；膨胀节施加的内部压力即环管反应器在正常工作时换热夹套的内压，此处考虑到工况下膨胀节的变形 图1 换热夹套上端膨胀节应力分析结果2023年05月|127()5253bpdESC eh C+=(14)组合应力，分别由公式(15)、(16)、(17)表示：p=2+3(15)d=4+5(16)R=0.7p+d(17)通过查阅 GB/T 1502011 压力容器 可知，膨胀节材料在设计温度下的屈服强 st=345 MPa、许用应力 t=189 MPa。通过应力计算，内压引起的直边段周向薄膜应力 z=115.33 MPat；内压引起的加强环周向薄膜应力 c=127.53 MPat；内压引起的膨胀节周向薄膜应力 1=66.41 MPat；内压引起的膨胀节经向薄膜应力 2=140.62 MPat；内压引起的膨胀节经向弯曲应力 3=6.25 MPa；轴向位移引起的膨胀节经向薄膜应力 4=42.02 MPa；轴向位移引起的膨胀节经向弯曲应力 5=277.02 MPa；组合应力p=146.87 MPa1.5st=517.5 MPa、d=319.04 MPa、R=421.849 MPa2st=690 MPa，因此通过采用第三强度理论对换热夹套上端膨胀节膨胀时内部应力水平的校核符合要求。通过 ANSYS 软件对换热夹套上端膨胀节进行应力分析后得到的应力结果最大值为 141.7 MPa，小于膨胀节材料的许用应力 t=189 MPa。验证了第三强度理论计算结果的准确性，说明了换热夹套上端膨胀节的应力水平低于理论应力允许值，在应力允许范围内，既提供了膨胀节的设计依据，同时又通过采用第四强度理论对环管反应器内外管温差叠加后的总应力水平进行了有效评价。2 温度对换热夹套下端膨胀节的影响研究图 2 所示为对环管反应器换热夹套下端膨胀节的应力分析结果，内容包括网格划分、施加载荷、位移约束、总变形云图、轴向变形云图以及等效应力云图。首先，在软件 ANSYS Static Structural 里的 Mesh Mechanical 下对换热夹套下端膨胀节进行结构化网格划分，网格划分总数为 105 12，并经过网格无关性验证，且网格质量均大于 0.8。其次，为了能够更准确的模拟出内筒体与换热夹套因受温差作用发生膨胀时，换热夹套对内筒体轴向膨胀上的约束作用，决定对膨胀节上端施加位移约束，对膨胀节下端施加垂直向上的力 F=9.5105 N(此处对换热夹套下端膨胀节施加的力大于对换热夹套上端膨胀节施加的力，是因为考虑了重力等因素)，施加内部压力 1.1 MPa。最后，应力理论”)对膨胀节膨胀时内压与轴向位移作用下产生的薄膜应力、弯曲应力进行研究；通过采用第四强度理论(也称“畸变能理论”)，并运用 ANSYS 软件模拟计算得到的等效应力云图，再次对膨胀节内部应力的分布情况进行校核。通过查阅 GB/T 167492018 压力容器波形膨胀节 可知，内压引起的直边段周向薄膜应力，由公式(3)表示，其中 D0、Dc、k 分别由公式(4)、(5)、(6)表示：()204402tbzttbccccpD L E kmSE L DS kE L D=+(3)0 bDDmS=+(4)2cbcDDmSS=+(5)41.5bLkD S=，当 k1 时，取 k=1 (6)内压引起的加强环周向薄膜应力，由公式(7)表示：()24402tcccttbccccpD L E kmSE L DS kE L D=+(7)内压引起的膨胀节周向薄膜应力，由公式(8)表示，其中 Dm由公式(9)表示：11220.571mppDhmSW|=|+|(8)0mDDh=+(9)内压引起的膨胀节经向薄膜应力，由公式(10)表示，其中 Sp由公式(11)表示：22pphmS=(10)0.5bpmDSSD=|(11)内压引起的膨胀节经向弯曲应力，由公式(12)表示：232ppphCm S=|(12)轴向位移引起的膨胀节经向薄膜应力，由公式(13)表示：()2432bpfESCeh C+=(13)轴向位移引起的膨胀节经向弯曲应力，由公式(14)表示：工艺管控128|2023年05月套下端膨胀节总变形云图、轴向变形云图、等效应力云图进行研究与分析。计算结果表明，换热夹套下端膨胀节在受载荷作用下的总变形最大值为 8.561 6 mm；轴向位移最大值为 2.031 1 mm，最小值为-1.265 5 mm，平均值为0.382 8 mm；等效应力最大值为 49.585 MPa，最小值为 0.748 31 MPa。并结合上节内容可知，换热夹套下端膨胀节的理论膨胀位移为 0.282 mm，与膨胀节在受载荷作用下的轴向位移平均值 0.382 8 mm 几乎相等，验证了换热夹套下端膨胀节应力分析的准确性。在对膨胀节膨胀时内压与位移作用下产生的薄膜应力、弯曲应力的分析研究中，通过运用上节公式(1.3)(1.17)分别计算出内压引起的直边段周向薄膜应力 z=102.23 MPat；内压引起的加强环周向薄膜应力 c=115.56 MPat；内压引起的膨胀节周向薄膜应力 1=57.84 MPat；内压引起的膨胀节经向薄膜应力 2=123.42 MPat；内压引起的膨胀节经向弯曲应力 3=28.22 MPa；轴向位移引起的膨胀节经向薄膜应力 4=10.25 MPa；轴向位移引起的膨胀节经向弯曲应力 5=232.72 MPa；组合应力 p=181.26 MPa1.5st=517.5 MPa、d=242.97 MPa、R=369.852 MPa2st=690 MPa，因此通过采用第三强度理论对换热夹套下端膨胀节膨胀时内部应力水平的校核符合要求。由 ANSYS 软件对换热夹套下端膨胀节进行应力分析后得到的应力结果最大值为 49.585 MPa，小于膨胀节材料的许用应力 t=189 MPa。验证了第三强度理论计算结果的准确性，说明了换热夹套下端膨胀节的应力水平低于理论应力允许值，在应力允许范围内，既提供了膨胀节的设计依据，同时又通过采用第四强度理论对环管反应器内外管温差叠加后的总应力水平进行了有效评价。3 温度对换热夹套连接管膨胀节的影响研究图 3 所示为对环管反应器换热夹套连接管膨胀节的应力分析结果，内容包括网格划分、施加载荷、位移约束、总变形云图、轴向变形云图以及等效应力云图。首先在软件 ANSYS Static Structural 里的 Mesh Mechanical 下对换热夹套连接管膨胀节进行结构化网格划分，网格划分总数为 35 760，并经过网格无关性验证，且网格质量均大于 0.8。其次对膨胀节中间位置施加位移约束，对膨胀节左端施加水平向右的压力F=1.38105 N，对膨胀节右端施加水平向左的压力F=1.38105 N，施加内部压力 1.1 MPa。最后在进行在进行求解计算前，选取总变形、轴向变形、等效应力对换热夹套下端膨胀节进行应力分析。应力分析时，膨胀节上端与反应器直管段连接部位被视为固定位置，是为了给膨胀节添加位移约束；膨胀节下端施加垂直向上的力，是为了简单还原内筒体与换热夹套在整体法兰处相连对膨胀节不受外力影响而发生自由膨胀时的约束作用所带来的约束力；膨胀节施加的内部压力即环管反应器在正常工作时换热夹套的内压，此处考虑到工况下膨胀节的变形与应力分布情况。通过以上对应力分析前的描述与分析，决定开始对环管反应器换热夹套下端膨胀节进行位移约束与施加载荷，并对求解计算后得到的换热夹 图2 换热夹套下端膨胀节应力分析结果2023年05月|129的内压，此处考虑到工况下膨胀节的变形与应力分布情况。通过以上对应力分析前的描述与分析，决定开始对环管反应器换热夹套连接管膨胀节进行位移约束与施加载荷，并对求解计算后得到的换热夹套连接管膨胀节总变形云图、轴向变形云图、等效应力云图进行研究与分析。计算结果表明，换热夹套连接管膨胀节在受载荷作用下的总变形最大值为 0.758 99 mm；膨胀节左端轴向位移最大值为-0.379 48 mm，最小值为-0.497 75 mm，平均值为-0.438 615 mm；膨胀节右端轴向位移最大值为 0.687 27 mm，最小值为 0.552 74 mm，平均值为 0.620 005 mm；膨胀节轴向总位移为 1.058 62 mm(因膨胀节左端轴向位移与膨胀节右端轴向位移方向相反，因此膨胀节轴向总位移等于膨胀节左右两端轴向位移平均值的绝对值之和)；等效应力最大值为118.93 MPa，最小值为 1.817 8 MPa。通过计算结果得出的数据可知，膨胀节在内部压力 1.1 MPa、左右两端压力 F=1.38105 N 下的膨胀情况。而筒体长度为3 500 mm 的换热夹套连接管在内部温度 45 (室外年平均气温为 20，温差为 25)下的膨胀情况可由公式(18)得出：P3=L3 T3=35 0001.2 10-525=1.05 mm (18)通过公式(18)可知，换热夹套连接管在正常工作时筒体的膨胀位移 1.05 mm 与膨胀节在受载荷作用下的轴向总变形量 1.058 62 mm 几乎相等，验证了换热夹套连接管膨胀节应力分析的准确性。在对膨胀节膨胀时内压与位移作用下产生的薄膜应力、弯曲应力的分析研究中，通过运用公式(3)(17)分别计算出内压引起的直边段周向薄膜应力 c=108.72 MPat；内压引起的膨胀节周向薄膜应力 1=42.26 MPat；内压引起的膨胀节经向薄膜应力 2=123.42 MPat；内压引起的膨胀节经向弯曲应力 3=3.56 MPa；轴向位移引起的膨胀节经向薄膜应力 4=22.72 MPa；轴向位移引起的膨胀节经向弯曲应力 5=195.77 MPa；组 合 应 力 p=126.98 MPa 1.5st=517.5 MPa、d=218.49 MPa、R=307.376 MPa 2st=690 MPa，因此通过采用第三强度理论对换热夹套连接管膨胀节膨胀时内部应力水平的校核符合要求。由 ANSYS 软件对换热夹套连接管膨胀节进行应力分析后得到的应力结果最大值为 118.93 MPa，小于膨胀节材料的许用应力 t=189 MPa。验证了第三强度理论计算结果的准确性，说明了换热夹套连接求解计算前，选取总变形、轴向变形、等效应力对换热夹套连接管膨胀节进行应力分析。应力分析时，膨胀节两端向里压缩，因此决定对膨胀节中间位置施加位移约束；膨胀节两端施加向里的压力，是为了简单还原筒体受热发生膨胀时对膨胀节的压缩作用所产生的压力以及膨胀节不受外力影响而发生自由膨胀时所受的约束力；膨胀节施加的内部压力即环管反应器在正常工作时换热夹套连接管 图3 换热夹套连接管膨胀节应力分析结果工艺管控130|2023年05月性进行整体评估，首先决定采用第三强度理论对膨胀节膨胀时内压引起的直边段周向薄膜应力、内压引起的加强环周向薄膜应力、内压引起的膨胀节周向薄膜应力、内压引起的膨胀节经向薄膜应力、内压引起的膨胀节经向弯曲应力、轴向位移引起的膨胀节经向薄膜应力、轴向位移引起的膨胀节经向弯曲应力、组合应力进行研究与校核。其次决定采用第四强度理论，并通过运用 ANSYS 软件对膨胀节应力分析后得到的应力结果再一次进行评估与校核。两次应力水平的校核结果均符合要求，说明了膨胀节的应力水平低于理论应力允许值，在应力允许范围内，既提供了膨胀节的设计依据，同时又对温差叠加后的总应力水平进行了有效评价。管膨胀节的应力水平低于理论应力允许值，在应力允许范围内，既提供了膨胀节的设计依据，同时又通过采用第四强度理论对环管反应器内外管温差叠加后的总应力水平进行了有效评价。4 结语本文对于环管反应器因受温度影响而发生热膨胀的研究中，通过对换热夹套上端膨胀节、换热夹套下端膨胀节与换热夹套连接管膨胀节的应力分析，不仅模拟出了膨胀节工作时的膨胀状态，还验证了膨胀节在外力与内压下发生的轴向变形量与筒体因受温度作用下的膨胀量几乎相等，同时也验证了应力分析结果的准确性。而为了能够对膨胀节的有效符号说明P1反应器内筒体膨胀位移mmL1反应器内筒体长度mmT1反应器内筒体内外温差P2换热夹套膨胀位移mmL2换热夹套长度mmT2换热夹套内外温差P3换热夹套连接管膨胀位移mmL3换热夹套连接管长度mmT3换热夹套连接管内外温差p设计压力MPam膨胀节的层数D0直边段平均直径mmDb膨胀节圆环内径mmL4直边段长度mmLc直边段加强环长度mmEb室温下膨胀节材料的弹性模量MPaEbt设计温度下膨胀节材料的弹性模量MPaEct设计温度下加强环材料的弹性模量MPaSp考虑到成形减薄后膨胀节一层的有效厚度mmS膨胀节一层的名义厚度mmSc直边段加强环的有效厚度mme一个波的轴向位移mmDc直边段加强环平均直径mmDm膨胀节平均直径mmDo换热夹套筒体外直径mmh膨胀节的波高mmW膨胀节一个波的波长mmC腐蚀裕量mmk、Cp、Cf、Cd系数z内压引起的直边段周向薄膜应力MPac内压引起的加强环周向薄膜应力MPa1内压引起的膨胀节周向薄膜应力MPa2内压引起的膨胀节经向薄膜应力MPa3内压引起的膨胀节经向弯曲应力MPa4轴向位移引起的膨胀节经向薄膜应力MPa5轴向位移引起的膨胀节经向弯曲应力MPat许用应力MPap、d、R组合应力MPast屈服强度MPa参考文献：1 孙胜仁.大型膨胀节设计方法研究D.北京：北京化工大学，2015.2 李进楠.超大型U 形膨胀节成形模拟及强度与刚度分析D.北京：北京化工大学，2021.3 李进楠，于洪杰，钱才富，等.超大型U 形膨胀节成形模拟与分析C.大连：第十六届全国膨胀节学术会议，2021.作者简介：王政文(1997-)，男，天华化工机械及自动化研究设计院有限公司，从事压力容器的研发与设计等工作，。