脉动流技术在管壳式换热器振动分析中的应用（修改）..docVIP

脉动流技术在管壳式换热器振动分析中的应用 管壳式换热器是一种广泛应用于石油、化工、机械、电力、食品等行业的工艺设备[1]。随着过程工业生产规模的逐渐扩大，管壳式换热器越来越趋向于大型化，管束无支撑跨度增大，刚性变差，为了提高传热效率，流体流速不断提高，并且新型传热技术[2-3]（如脉动流）的广泛研究，使换热器流体诱导振动问题越来越突出[4]。可见，为了实现脉动流作用下提高换热器的换热效率，并使该设备能正常工作，则对脉动流作用下换热器的振动分析具有重要的意义。 本文以脉动流换热器为研究对象进行振动分析， （-1） 式中，— 结构的质量矩阵；— 结构的阻尼矩阵；— 结构的刚度矩阵；— 结构的位移向量；— 结构的载荷向量. 当结构系统无外载作用若结构在简谐载荷作用下谐响应分析若结构在任意随时间变化的载荷作用下瞬态动力学分析 (1-2) 式中，—流体的质量阵；—流体的阻尼阵；—流体的刚度阵；—的形函数； —流体边界的法向向量； —位移分量、、的形函数。 —耦合质量阵；—耦合质量阵；—耦合刚度阵. 2.模态分析 2.1 结构尺寸与工艺参数 该换热器筒体长度为1800mm，换热管规格为Ф25x0.8mm，换热管弹性模量为2.1x1011Pa，泊松比为0.3，密度为7850kg/m3。其结构图和实验室设备,如下图1和图2所示： 图1 结构示意图 图2 实验室试验设备 换热器管程上加入旁路，并在旁路上安装脉动流发生器，无脉动流条件下，壳程入口速度为1.5m/s，管程入口速度为0.5m/s，模拟条件与实验未加脉动流时条件一致；脉动流条件下，换热管管程入口速度随时间作周期性变化，脉动流函数如式（1-）1.5m/s，模拟条件与实验加入脉动流时条件一致。 （1-） 根据换热器的实际模型，建立如图3所示的流体模型和换热管模型。为了保证计算精度的前提下节约计算机资源，选用非结构化网格对流体模型进行网格划分，划分后生成的网格数约为20万个，经验证该网格已收敛，其网格模型如图4所示[9]。 图3 实验换热器流固耦合模型 图4 流体网格模型 2.3 加载与求解 首先在Fluent中进行流场分析，模型选用模型，定义物性参数为水，整个区域设置为流体域。再将管、壳程入口定义为速度入口边界条件，并指定入口速度大小；对于加入脉动流后的流场分析需要通过UDF的方式改变管程入口速度。同时设置回流条件（由湍流强度和水力直径控制），壳程出口为压力出口边界条件，壳壁设置为壁面边界。 流体分析结束后，打开结构分析模块，将模型中流体模型去掉，为换热管模型划分网格，得到换热管有限元模型。由于换热管两端固定，则应对换热管两端所有节点施加全约束，再将流体分析的结果通过流固耦合交界面加载到换热管模型中，设置分析类型为模态，提取前五阶模态，采用分块兰索斯法对模型进行求解。 2.4 模态结果分析 图5与图6分别为有无脉动流作用下的换热管前两阶模态振型图，从图5中可以看出，换热管中间段相对位移最大，振型显示一个波峰；从图6中可以看出，换热管振型左右两侧分别出现一个波峰和一个波谷。 a.无脉动流 b.有脉动流 图5 换热管第一阶模态振型图 a.无脉动流 b.有脉动流 图6 换热管第二阶模态振型图 通过对换热管在无脉动流和有脉动流作用下进行模态分析，得到两种工况下的固有频率计算结果如表1所示。从表1中可以看出，两种工况下固有频率相差很小，说明外界流场的变化对换热管在流体环境中的固有振动特性的影响很小。 表1 有/无脉动流作用下换热管固有频率表 无脉动流 有脉动流 一阶固有频率/Hz 47.365 47.372 二阶固有频率/Hz 130.10 130.12 三阶固有频率/Hz 253.85 253.88 四阶固有频率/Hz 417.18 417.22 五阶固有频率/Hz 618.87 618.93 3.动态响应分析 3.1无脉动流时换热管谐响应分析 在无脉动流作用下，参照相关公式[10]，计算得到加载在换热管入口和出口处的载荷分别为10.8N和7.7N从图中可以看出，换热管在流体激振力作用下，其第一、三阶固有频率处振动响应十分强烈，共振能量很大。距固定端0.15m处的节点在第一、三阶固有频率处振幅较大，最大振幅值达到0.43mm，而在第二阶固有频率处产生的振幅相对较小，随着节点越靠近中间位置处，振动幅值将达到最大，其值为2.1mm。 图7 距固定端0.15m处的频率振幅曲线 图8 换热管中间位置的频率振幅曲线 对比图9和图10可以看出，换热管端部在第三阶固有频率处应力最大，其值达到95MPa，可见，在换热管端部最易出现强度破坏