1、不锈钢舱化学品船温度应力分析房玉慧 马 骏大连理工大学摘要：本文以某不锈钢舱化学品船为研究对象，借助于大型工程软件 ANSYS，建立了 2 舱段有限元模型，然后进行温度场求解，算出等效节点温度载荷后，将热分析单元转换成结构分析单元，将计算出的温度载荷作为体载荷施加到船体结构上，再施加外力载荷，进行热与结构应力耦合计算。最后，分析结果。结果表明：加强筋由于温度载荷的影响，变形很大，并且在端部产生了应力集中，因此应适当优化加强筋形状。关键词：不锈钢舱；温度场；热应力；有限元分析1 引 言化学品船作为一种运载危险液体货物的特殊船舶，经常运输一些对人体或环境有害的液体或气体，这些物质在运输途中一般需要
2、保持一定的特殊存贮环境，例如运输沥青、木馏油、硫黄等化学物质时，通常需要持续加热并使其维持在高温状态(90-250 )，以防止货物凝固。这就要求化学品船能够经受住长时间高热源作用下引起的温度应力，保证船体结构不发生屈服或开裂破损等问题。因此对化学品船进行专门的温度场及温度应力分析，确保其符合强度要求，对保证其在运输过程中的安全具有重要意义和价值。2 计算原理2.1 温度场计算原理热分析遵循热力学第一定律。即能量守恒定律。稳态热分析可以通过有限元计算确定由于稳定的热载荷引起的温度、热梯度、热流率、热流密度等参数。如果系统的净热流率为 0，即流入系统的热量加上系统自身产生的热量等于流出系统的热量：
3、，则系统处于热稳态。稳态热分析的能量平衡方程为：q流 入 生 成 流 出 KTQ式中： 为传导矩阵； 为节点温度向量； 为节点热流率向量。KT、 、 由 ANSYS 自动生成，前提条件是输入模型几何参数、材料热性能参数以及所TQ施加的边界条件。2.2 温度应力计算原理有限元法的分析过程，概括起来可分为下列六个步骤：(1) 结构离散化(2) 选择位移模式(3) 分析单元的力学特性(4) 计算等效节点力(5) 集合所有单元的刚度方程建立结构平衡方程(6) 将边界条件代入，求解未知节点位移并进行单元分析3 计算模型及计算条件3.1 有限元模型计算模型均包括 1 个整舱和 2 个半舱（1/2+1+1/
4、2 模型）内的全部结构构件（包括不锈钢舱壁） 。计算模型见图 3-1、图 3-2 和图 3-3 所示。图 5.3 有限元计算模型图 3 温度场物理模型3.2 材料参数计算采用的材料特性及主要参数见表 1 和表 2表 1 船体材料特性20 C S31803 双相不锈钢 普通钢密 度 8000 kg/m3 7850 kg/m3弹性模量 2.01011 2.11011泊松比 0.3 0.3热传导系数 20 W/(m2 C) 60 W/(m2 C)热膨胀系数 1.310-6 C-1 1.310-6 C-1表 2 温度场计算的主要参数项 目 数 值货物对内底板、舷侧内板的对流换热系数 100 W/(m2
5、 C)海水对外底板、舷侧外板的对流换热系数 8606 W/(m2 C)空气对甲板、舷侧外板的对流换热系数 23.3 W/(m2 C)舷边舱内空气对舷侧外板的对流换热系数 10 W/(m2 C)舷边舱内空气对甲板的对流换热系数 10 W/(m2 C)货物温度 130 C，恒温海水、空气温度 0 C，恒温3.3 计算条件本文在求解温度场时，忽略辐射的影响，考虑的热传递过程具体如下：1. 高温货物以自然对流方式将热量传递给甲板、舱壁和内底板等；2. 舷边舱内空气以对流方式与舷侧内板、舷侧外板换热；3. 外界空气与海水以强迫对流方式与甲板、舷侧外板和外底板进行热交换；4. 结构构件之间以热传导方式进行
6、热传递；5. 对于双层底内空气与内底板、外底板之间不发生对流换热。在进行结构应力分析时，计算模型的边界条件采用 CCS 规范规定的边界条件3.4 计算工况根据规范要求的一道纵舱壁油船所规定 9 种计算工况（见表 3.1）和所考虑附加的危险计算工况以及强度校核所得结论，在不锈钢舱热应力分析中，选取 1、5、9 等三种主要结构构件合成应力最大的计算工况为分析工况。这三种工况对应不同的装载要求，使两个不锈钢舱分别为满载和部分装载；同时不锈钢舱内设定最高液货温度为 130 度，相邻普通液货舱温度为 0 度，这时不锈钢舱内和相邻液货舱船体结构构件承受的温度应力和结构耦合应力均达到最大值。所以按这三种工况
7、校核所得到结构耦合应力应为船体构件在营运中所承受载荷的峰值应力，这样得到的分析结果对实际工程是偏于安全的。表 3.1 一道纵舱壁油船计算工况4 计算成果本文建立两舱段的有限元计算模型，采用 CCS 船级社的计算处理方法，用 ANSYS 对该有限元模型进行了详细的计算。计算得到大量的数据结果。4.1 温度场与温度应力计算情况图 4.1 船体温度场云图（满载） 图 4.2 船体温度场云图（部分装载）Fig.4.1 Nephogram of Hull Temperature Field Fig.4.2.Nephogram of Hull Temperature Field(Full Load) (P
8、artial Load)图 4.3 船体结构应力云图， （Seqv,满载） 图 4.4 船体结构应力云图（Seqv，部分装载）Fig.4.3.Nephogram of Structural Stress Fig.4.4.Nephogram of Structural Stress(Seqv，F ull Load） （Seqv，Partial Load ）图 4.5 船体结构应力云图（Sx，满载） 图 4.6 船体结构应力云图（Sy，部分装载）Fig.4.5. Nephogram of Structural Stress Fig.4.6. Nephogram of Structural Stre
9、ss(Sx，Full Load） ） （Sy，Partial Load）将船体构件的温度场及温度应力绘出曲线如下所示：-200204060801001200 6.66 13.32 19.98 26.64纵 向 坐 标 （ m)温度（）内 底 板 纵 向 温 度 分 布甲 板 纵 向 温 度 分 布图 4.7 温度场纵向分布曲线Fig.4.7 Longitudinal distribution of temperature field -20020406080100120140-20 -10 0 10 20横 向 位 置 (m)温度（）内 底 横 向 温 度 场 分 布甲 板 横 向 温 度 场
10、 分 布图 4.8 温度场沿横向分布曲线Fig.4.8 Transversal distribution of temperature field 0501001502002503003504000 6.66 13.32 19.98 26.64纵 向 位 置 (m)温度应力(MPa)内 底 纵 骨 温 度 应 力甲 板 纵 骨 温 度 应 力图 4.9 温度应力沿纵向分布曲线Fig.4.9 Longitudinal distribution of thermal stress 050100150200250300350400450-20 -10 0 10 20横 向 位 置 (m)温度应力(M
11、Pa)内 底 温 度 应 力 沿 船 宽 分 布甲 板 温 度 应 力 沿 船 宽 分 布图 4.10 温度应力沿横向方向分布曲线Fig.4.10Transversal distribution of thermal stress 从图中可以看出，温度场在纵向、横向分布各不相同，但是温度场关于模型中横剖面和中纵剖面呈对称分布。从总体来看，船体内部靠近高温液货的构件温度较高，受到温度载荷的影响很大。下面列出温度和结构载荷的耦合应力表 4.11 耦合应力结果（工况一）Tab.4.11Results of coupling stress (loading condition one)最大应力（MPa
12、）计算值 许用值工况一 位置Sx Sy Sz S甲板 289 180 206甲板纵骨 342 243内底板 369 359 206内底纵骨 299 243纵舱壁 133 54 90 206横舱壁 100 139 158 206满载不锈钢舱壁 165 27 127 206甲板 172 120 206甲板纵骨 208 243内底板 341 309 206部分装载内底纵骨 302 243纵舱壁 133 54 90 206横舱壁 100 139 158 206不锈钢舱壁 151 18 73 206另外，本文专门对加强筋进行了应力分析，结果如下：图 4.12 应力云图（Seqv） 图 4.13 位移云图
13、（Usum）Fig.4.12. Nephogram of Stress（ Seqv） Fig.4.13. Nephogram of Displacement（Usum）分析结果表明：加强筋左部靠近高温液货的地方，受到较大的温度应力，甚至两端出现了温度应力集中，从位移云图中可以得到，由于温度载荷作用，加强筋的变形很大，但是，加强筋的存在，减小了温度载荷对其他连接构件的影响。5 结 论本文利用大型通用有限元分析软件 ANSYS 对其进行不同温度场下的热应力分析计算。研究了同一不锈钢液货舱内装载不同温度的高温货物和相邻液货舱内装载不同温度高温货物两种条件下产生的不均匀温度场的分布，以及由此两种不均匀
14、温度场引起的船体结构内部热应力分布及变化状况，探讨由此引起的船体局部强度，给出了船体结构所能承受的热载荷极限值。在一定假设条件下，根据实际模型，进行了有限元计算，得到了大量计算数据，完成了不锈钢舱化学品船的温度场和温度应力分析，计算结果表明：接触高温液货的船体构件或靠近高温液货的构件，不仅温度较高，而且温度应力也较大，某些构件受到的耦合应力甚至超过许用应力，因此温度应力不能忽视。槽型舱壁对减小温度应力作用明显，起缓冲作用。对于与不锈钢舱壁连接的加强筋，受温度影响也较大，产生较大的变形，因此造成在内底板上的应力集中，应考虑适当优化加强筋形状。参考文献:1 陈伯真，胡毓仁.船体温度分布及温度应力计算J. 上海交通大学学报， 1995，29(3):33-41.2 陈辉，袁海庆等.温度作用对钢筒仓结构影响的有限元分析J.四川建筑科学研究，2006，32(6)：47-49.3 滕晓青，顾永宁.双壳型船体结构稳态温度场和温度应力 J.中国造船，2000，41(2):59-64.4 Det Norske Veritas. Hull structural design ships with length 100 metres and above , Part 3 Chapter 1, 1997
