《武汉工程大学学报》 2014年04期
55-59
出版日期：2014-04-30
ISSN:1674-2869
CN:42-1779/TQ
脉动流条件下弹性管束换热器振动特性仿真
0引言弹性管束换热器是利用流体诱导振动来达到强化传热的效果，并且还有利于降低污垢热阻，实现复合强化传热［1］.通过诱发和控制传热元件自由振动达到流体诱导振动强化传热来防止振动破坏，因此控制换热器内部传热元件的振动至关重要［23］.为了实现换热器的强化传热和流体诱导强化传热技术，研究如何合理地诱发传热元件振动以及实现对传热元件的有效控制有着重要的意义［4］.由于脉动流的脉动参数具有周期性变化的特点，可以将其作为对传热元件进行简谐激励的外部激励源，以实现对传热元件振动的控制［56］.为了控制传热元件振动的效果，采用脉动流发生装置来诱导弹性管束产生一定频率的周期性振动［78］.然后利用得到的周期性振动的脉动流流体流经弹性管束，进而研究弹性管束换热器的振动频率.1弹性管束结构及有限元模型弹性管束由四根弯管和两块不锈钢连接体组成，其三维实体模型如图1所示.连接体5、6两处为固定端，A、B两端可以自由运动.管束从外到内的编号分别为4、3、2、1，换热器工作时管程流体可以从A、B两端的任意一端进入，依次经过四根弹性管束，连接体将相邻管束连通.利用ANSYS软件建立了弹性管束有限元模型，连接体选用六面体实体单元Solid45建模，管束选用板壳单元Shell93建模，网格划分后得到的有限元模型如图2所示. 图1弹性管束结构图Fig.1Schematic drawing of elastic tube bundle structure 图2弹性管束的有限元模型Fig.2Finite element model of elastic tube bundle2模态分析及计算结果利用ANSYS软件中结构力学的模态分析模块对管束进行分析，得到平面弹性管束结构的各阶固有频率和相应振型.管束各部件的几何尺寸和材料属性列于表1和表2中.表3给出了平面弹性管束结构前10阶固有频率和振型，从计算结果可以看出，平面弹性管束的振型存在两种：在管束平面内的面内振型和垂直于管束平面的面外振型.实际工况中，在内外流体的共同作用下，弹性管束呈现出的是一种复合而成的复杂三维运动，并不是简单的处于一种面内振动状态或者面外振动状态.周围流体由于管束这种复杂的三维振动状态会产生扰动，这有利于边界层厚度减薄，进而导致热阻减小，对流传热系数增加，换热性能提高.而且管束振动变形使管束表面污垢不易附着和促进污垢的脱落，导致污垢热阻减小，最终实现复合强化传热.第4期王成刚，等：脉动流条件下弹性管束换热器振动特性仿真武汉工程大学学报第36卷表1弹性管束的几何尺寸Table 1Geometric dimensions of elastic tube bundlemm尺寸名称数值钢管截面半径10钢管厚度1.5钢管1半径75钢管2半径95钢管3半径115钢管4半径135连接体5各边长40×20×20连接体6各边长80×20×20表2弹性管束的材料属性Table 2material properties of elastic tube bundle材料属性名称数值钢管的弹性模量/GPa129钢管的泊松比0.33钢管的密度/kg·m-38 900连接体的弹性模量/GPa210连接体的泊松比0.3连接体的密度/kg·m-37 800表3弹性管束结构模态Table 3The structure mode of elastic tube bundle阶次频率/Hz 主要变形方式122.444面内222.900面外330.179面外434.408面内544.383面外660.107面内789.196面外899.329面内9155.800面内10167.570面外3脉动流发生装置三维仿真计算考虑到是研究无源强化传热形式以及实际生产中安装的可行性，因此产生所需脉动流采用的是在分支出口处安放绕流体的方式.流体流经一定形状的绕流体后内部存在旋涡，会产生旋涡分离，在不同位置流体的速度值和压力值是不同的，导致尾流中流体的流速、压力、密度产生周期性变化，从而出口处流体的速度、压力会随时间而变化，形成脉动流.然而对于不同绕流体的截面形状所产生的脉动流强度和稳定性也不同.斯特罗哈尔数St是绕流体最主要的特性参数，可根据涡脱落频率计算得到：St=fsdu1=fsmdu其中：fs为漩涡脱落频率，Hz；d为绕流体迎流面特征宽度，m;u1为绕流体两侧的平均速度，m/s;m为绕流体两侧流通面积与管道流通面积之比；u为管道内平均流速，m/s.选用三棱柱绕流体和圆柱绕流体，m=0.72.各个结构尺寸：直管部分内径为30 mm、长度为590 mm，分支部分直管管束内径为12 mm、长度为24 mm、三棱柱绕流体横向特征宽度为2.8 mm，三棱柱迎流面距离分支入口为12 mm；圆柱绕流体的直径为4 mm，圆柱迎流面距离分支入口为12 mm.分支部分的有限元模型如图3所示:图3（a）为三棱柱绕流体分支部分有限元模型、图3（b）为圆柱绕流体分支部分有限元模型. （a） （b）图3分支部分的有限元模型Fig.3Finite element model of branch part按照上述的结构尺寸利用FLUENT计算软件，对脉动流发生装置进行仿真计算，在脉动流发生装置入口处流速设为0.4 m/s，在分支出口处检测速度变化情况，得到如图4所示的曲线：图4（a） 为三棱柱绕流体监测点的速度变化曲线、图4（b） 为圆柱绕流体监测点的速度变化曲线. 可以看出两种绕流体监测点的速度变化曲线波形明显并且具有一定的稳定性.三棱柱绕流体的速度变化范围在0.4~1.3 m/s之间，圆柱绕流体的速度变化范围在0.07~0.15 m/s之间，相比较而言三棱柱绕流体产生的波形速度变化范围更大并且速度变化曲线的波形也比较稳定，因此三棱柱绕流体更适合用于脉动流发生装置. （a） （b）图4监测点的速度变化曲线 Fig.4The speed curve of monitoring stations 三棱柱绕流体的速度变化频率由图4（a）可计算得到基本为41.67 Hz.此时的St值为0.21，与相关研究中St所给的范围比较一致.4装置整体的实际频率 利用脉动流发生装置产生的脉动流流体流经弹性管束换热器，对弹性管束的振动频率进行分析计算.通过变频器调节泵的流量来实现脉动流发生装置入口流速的变化，利用FFT(快速傅里叶变换)得到管束的振动频谱图.选用脉动流发生装置入口速度为0.4 m/s时产生的脉动流流体流经弹性管束换热器，模拟了第一排平面管束及换热器的振动频谱图，如图5(a)、5(b)所示.从图5(a)可以看出管束的振动频谱主要包括30 Hz、42 Hz、52 Hz三个频率，从图5（b）可以看出：管束振动频谱中的30 Hz和52 Hz是由外界干扰所产生，由脉动流引起的振动频率为42 Hz，这一数值与之前的结果基本一致. （a） （b）图5振动频率Fig.5Vibration frequency 5结语a.利用ANSYS软件对弹性管束进行模态分析，得到弹性管束的振动是一种复杂的三维运动.这无疑对提高弹性管束的换热性能和抗结垢能力产生重要影响.b.在平面弹性管束分支出口安放绕流体可以产生脉动流，采用安放三棱柱绕流体能够产生波形较好的脉动流.因此三棱柱绕流体更适合用于该脉动流发生装置.c.利用FLUENT软件对脉动流发生装置进行仿真计算得到的频率与实际情况下脉动流流体流经弹性管束换热器引起的振动频率相一致.致谢感谢国家自然科学基金委员会和武汉工程大学研究生创新基金的资助！