首页/文章/ 详情

管道热应力分析原理详解(理论计算与workbench实现)

2月前浏览1307
正文共: 1649字 14图     预计阅读时间: 5分钟

1 前言

管道热应力分析是高温系统设计的一个重要工作,特别是对于SOFC系统,设备&管道的运行温度可达800℃以上。这种情况下,设备&管道的热应力分析很重要,因为蠕变将会是一个重要问题,设备管路集成设计时,需要尽可能降低热应力水平。热应力分析可能对于大多数工程设计人员而言,就是在专业软件上熟悉一下操作,因为整个热应力分析流程已经被设计成流水线模式。但是一个合格的设计人员,对于这种分析更为底层的信息笔者觉得有必要去研究研究。今天,我们用一个简单的管道模型,完整地演示一下热应力分析的底层逻辑。

2 模型问题描述

考虑如下一根钢管,规格为φ73×3,长度L=1000mm,不考虑内部压力,假设管道温度由22℃升高到900℃,分别计算不同的力学边界条件下管道的热应力情况。我们考虑三种力学边界条件:1)两端固定;2)一端固定,一端自由;3)一端固定,另一端连接设备,并且设备放置于滑动支座上,滑动摩擦系数μ=0.3,设备质量60kg。
假设钢管的在900℃的弹性模量E为2e5MPa,22℃到900℃范围内的线膨胀系数α为1.2e-5mm/mm.℃

3 理论计算

我们先采用理论方法计算上述几种条件的热应力。
首先,对于边界条件2),很显然,管道是自由的。因此,轴向上的应力为0MPa。同时,管道受热需要膨胀,其自由膨胀长度△L=αL△t=1.2e-5×1000×(900-22)=10.536mm。当然,管道两端对支座的作用力为0。
其次,对于边界条件1),很显然,由于两端固定,因此管道的膨胀长度△L=0。但是,管道所受的热应力,相当于L+△L长度的管道被压缩成L长度而造成的轴向应力。所以,管道的应变为ε=△L/L,而应力可以根据胡克定律计算得到,即σ=εE=0.010536×2e5=2107.2MPa。管道两端对支座的作用力就等于管道的轴向力,亦即F=σA=2107.2×π/4×(73^2-67^2)= 1390192.448N
再次,对于边界条件3),相当于边界条件1)的弱化情况,介于边界条件2)和边界条件1)之间。我们首先对管道进行受力平衡分析:管道受热膨胀后会对设备形成推力F,这个推力起始阶段会大于设备的静摩擦力f,设备发生位移,直至推力F刚好等于静摩擦力f。此时管道所受的轴向力(推力)F=f=μmg=0.3×60×9.81=176.58N,因此管道的轴向应力即热应力σ=F/A=176.58/[π/4×(73^2-67^2)]= 0.2676MPa。管道的应变ε=σ/E=1.34e-6,因此由于设备约束推力形成的位移为εL=1.34e-3mm,管道总变形的自由膨胀的位移减去设备约束推力位移:△L=αL△t-εL=10.535mm。
可见,要减小管道的热应力,减弱其约束是根本。采用滑动支座是有效的弱化约束的方法,另外减小支座摩擦系数也很重要。

4 workbench实现

接下来,我们在workbench平台分别计算上述三种工况。注意,我们创建管道1/4对称模型,便于边界条件设置。
边界条件2),管道一端设置位移约束,轴向位移为0,其余方向自由;1/4剖面设置无摩擦约束,亦即对称边界;管道本体温度900℃。
管道的变形和应力云图如下,另一端位移10.536mm,管道的轴向应力几乎为0,与理论计算一致。
边界条件1),管道两端设置位移约束,轴向位移为0,其余方向自由;1/4剖面设置无摩擦约束,亦即对称边界;管道本体温度900℃。
管道的变形和应力云图如下,管道轴向位移几乎为0,管道的轴向应力2107.2MPa,管端对支座的推力为1390200N(由于1/4模型,因此总推力需要乘以4),与理论计算一致。
边界条件3),管道一端设置位移约束,轴向位移为0,其余方向自由;管道另一端施加44.145N压力(由于1/4模型,因此推力需要除以4);1/4剖面设置无摩擦约束,亦即对称边界;管道本体温度900℃。
管道的变形和应力云图如下,另一端位移10.535mm,管道的轴向应力约为0.267MPa,与理论计算一致。

来源:仿真与工程
Workbench理论管道
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-09-15
最近编辑:2月前
余花生
签名征集中
获赞 195粉丝 372文章 306课程 0
点赞
收藏
作者推荐

FLUENT换热器模型相似试验模拟

正文共: 1083字 7图 预计阅读时间: 3分钟1 前言基于相似原理开展模型试验是装备开发的一个非常高效、经济的手段,之前我们做过一个相似原理的模拟(2022年6月19日推文,“一个多此一举的模拟——相似原理”)。今天,我们继续探讨一下相似原理,今天的案例相对复杂一点,以换热器为对象。2 建模与网格我们创建如下的套管换热器,逆流形式,长度0.6m,内外径分别为Φ40mm和Φ60mm,忽略壁厚的影响。划分六面体结构化网格,节点数约15.3玩,最小正交质量0.8。3 求解与讨论3.1原型假设换热器的性能是按下列参数设置的计算结果,物性参数是基于定性温度的值(进出口平均温度)。注意由于雷诺数很小,按层流计算。3.2 模型1假设由于某种原因,需要按5:1的比例尺开展缩小模型的试验,介质还是原来的介质不变。对于传热问题,由于努塞尔数(换热系数)通常是雷诺数和普朗特数的关联式,因此这种情况相对而言是简单的,因为物性参数不变,那么普朗特数自然是相似的,只需要调整流速即可实现雷诺数相似。雷诺数与特征长度成比例,因此为了使得雷诺数相似,模型缩小5倍后,流速应当放大5倍。因此,我们在FLUENT对模型进行笛卡尔坐标系下三个坐标轴按比例缩小5倍,同时修改流速至5倍大小,后重新计算。模型1的计算结果如下,可用看出各温度结果和原型完全一致,因此对数温差LMTD也一致,而流量随流速成比例变化,因此功率也缩小5倍,而换热面积缩小25倍,因此传热系数U放大了5倍,两侧的膜系数也按比例变化,这些变化规律可通过量纲分析来阐明。又由于努塞尔数和水力直径成比例(Nu=h*D/k),因此模型1的努塞尔数和原型是一致的。模型1和原型的阻力系数完全一致。综上,可用看出模型1和原型是相似的。3.3 模型2假设由于某种原因,无法采用原型的高温介质开展试验,但是尺寸不缩放,这种情况难度相对模型1较高,假设有如下的介质,恰好很好地满足普朗特数相似。这种其实是很难的,因为物性参数本身就与温度相关,因此需要找到定性温度下的普朗特数和原型基本相等的介质。而雷诺数相似却较好实现,因为只需要调整流速即可。我们看一下这种情况下的计算结果,很显然,换热系数、膜系数、努塞尔数和原型有较大差距,摩擦系数可以认为相似,因此模型2和原型不相似。我们分析一下,由于摩擦系数和雷诺数相关,因此雷诺数相似则摩擦系数相似,结果的偏差来自于计算误差。综上可知,对于换热器的整体模型试验,采用原型介质开展模型缩比试验很容易实现;但是采用原型尺寸,换其他介质开展试验是不现实的。不过,当研究对象不是换热器整体,而是研究冷侧或热侧的膜系数的相似试验,又是另外的情况了,我们以后讨论。来源:仿真与工程

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈