一种平面管束流体弹性失稳评定方法技术

本发明专利技术提供的一种平面管束流体弹性失稳评定方法包括模态分析和Connors公式计算两部分。模态分析部分负责计算换热管的固有频率和对应的振型；Connors公式计算则根据模态分析结果计算不同模态下的等效流速和临界流速。本发明专利技术具有如下优点：利用有限元本身的离散特点，可以实现不同管外流体参数、非均匀横向流速、复杂换热管形状的模态分析计算；在有限元模态分析结果的基础上，采用高斯积分方法对Connors公式进行计算，实现了振型对等效流速的加权效应。

【技术实现步骤摘要】
一种平面管束流体弹性失稳评定方法
本专利技术涉及一种对平面管束流体弹性失稳进行评定的方法。
技术介绍
管壳式换热器广泛应用于电力、化工领域，主要利用换热管束实现管内外流体介质的大面积热交换。管外流体在流经管束时，其流动方向可分解为纵向和横向。如果横向流速过大，则会产生流体弹性失稳现象，导致换热管大幅振动并最终破坏失效。流体弹性失稳是一种流固强耦合问题。每根换热管的运动会影响流体作用力和其他换热管的运动，并产生自激力。对于失稳现象的发生，一般认为是流固耦合作用对结构质量、阻尼和作用力产生了不利影响。由于其机理复杂，目前没有非常完善的理论，大部分预测失稳发生的公式都是通过试验测量得到的经验公式。热交换器标准GB/T151-2014的附录C对流体诱发振动的计算和校核给出了一系列计算公式和判定准则，但使用范围仅限于管外单一介质均匀流过直管或U形管的问题，对换热管的支撑形式有明确的要求。ASMEBPVC-IIIAppendices中的N-1330章节介绍了流体弹性失稳产生的一些特征，简要给出了推荐公式以及经验常数的选择，但没有对具体问题给出明确的计算方法。使用上述标准往往会产生以下困难：1)无本文档来自技高网...

【技术保护点】
1.一种平面管束流体弹性失稳评定方法，其特征在于，包括以下步骤：步骤1、利用换热管的轴线，将换热管简化为三维空间中的一维曲线模型；将一维曲线模型拆分为若干线段，须保证离散线段模型与原模型在几何上基本一致；在空间中定义三维笛卡儿坐标系，对节点坐标、单元节点编号、节点约束进行列表，以便后续计算，其中，将离散线段定义为单元，将离散线段的端点定义为节点，单元的节点包括离散线段的两个端点以及线段中点，共三个：节点坐标列表LOC表示为：

【技术特征摘要】
1.一种平面管束流体弹性失稳评定方法，其特征在于，包括以下步骤：步骤1、利用换热管的轴线，将换热管简化为三维空间中的一维曲线模型；将一维曲线模型拆分为若干线段，须保证离散线段模型与原模型在几何上基本一致；在空间中定义三维笛卡儿坐标系，对节点坐标、单元节点编号、节点约束进行列表，以便后续计算，其中，将离散线段定义为单元，将离散线段的端点定义为节点，单元的节点包括离散线段的两个端点以及线段中点，共三个：节点坐标列表LOC表示为：式中，xi、yi、zi分别为第i个节点的x轴、y轴、z轴坐标；单元节点编号列表NOD表示为：式中，分别为第e个单元的第1、2、3个节点的编号，第e个单元的第1、2、3个节点为第e个单元的两个端点及一个中间节点；节点约束列表UROT表示为：式中，uxi、uyi、uzi、rotxi、rotyi、rotzi分别为第i个节点的6个自由度约束；步骤2、根据Timoshenko梁理论、有限元理论以及数值积分方法，计算单元质量矩阵Me及单元刚度矩阵Ke：式中，ngaussm和ngaussk分别为质量阵和刚度阵的高斯积分点数量，ξi为高斯点位置，wi为加权；je＝le/2为标准域变换至单元局部坐标的Jacobian行列式，le是单元线段的长度：N(ξ)是关于标准域坐标ξ的矩阵函数：式中，为单元节点1的形函数；为单元节点2的形函数；N3＝N3(ξ)＝1-ξ2为单元节点3的形函数；矩阵C是Me的截面相关矩阵：式中，ρ，A，J，Iy，Iz分别为单元密度、截面积、扭转惯性矩、y向弯曲惯性矩、z向弯曲惯性矩；B(ξ)是关于标准域坐标ξ的矩阵函数：式中，矩阵D是Ke的截面相关矩阵：式中，E、Ky、Gy、Kz、Gz、Kx、G分别为杨氏模量、y向剪切修正系数、y向剪切模量、z向剪切修正系数、z向剪切模量、扭转修正系数、剪切模量；矩阵T是单元局部坐标至全局坐标的变换矩阵：式中，步骤3、将步骤2中计算的单元质量矩阵Me及单元刚度矩阵Ke按照节点自由度的编号依次累加进整体质量矩阵M及整体刚度矩阵K中：式中的并非求和符号，而是代表整体...

【专利技术属性】
技术研发人员：张马骏，杨志刚，陆毓颖，王佳欢，杨洁，
申请(专利权)人：上海电气电站设备有限公司，
类型：发明
国别省市：上海,31