发明名称 一种换热器传热管漩涡脱落诱发振动的分析方法
摘要 本发明公开了一种换热器传热管漩涡脱落诱发振动的分析方法,包括以下步骤:获得换热器一次侧及二次侧流体沿传热管分布的流场参数,基于流场参数和传热管的结构参数判断流场参数属哪一范围,若流场参数属于第一预设范围,则利用流体力系数数据库获得具体流场参数下的流体力系数;若流场参数属于第二预设范围,则利用公式(1)~(4)计算换热器传热管的振动响应;若流场参数属于所述第三预设范围,则利用公式(1)、(2)、(5)、(6)计算换热器传热管的振动响应。本发明的有益效果是:既可以考虑传热管内流体的流动,也可以考虑传热管内流体-传热管外流体-传热管振动的耦合作用,便捷、高效地实现了传热管的漩涡脱落诱发振动计算。
申请公布号 CN105205258A 申请公布日期 2015.12.30
申请号 CN201510605752.5 申请日期 2015.09.22
申请人 中国核动力研究设计院 发明人 冯志鹏;姜乃斌;臧峰刚;张毅雄;齐欢欢;黄旋
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 成都行之专利代理事务所(普通合伙) 51220 代理人 王记明
主权项 一种换热器传热管漩涡脱落诱发振动的分析方法,其特征在于,包括以下步骤:S1、利用热工水力系统分析软件对换热器的一次侧及二次侧流场进行计算,获得所述一次侧及二次侧流体沿传热管分布的流场参数,该流场参数包括:传热管内流体与传热管外流体的密度和流速;S2、基于所述流场参数和所述传热管的结构参数,判断所述流场参数是否属于第一预设范围;S3、若所述流场参数属于所述第一预设范围,则利用流体力系数数据库通过插值获得具体流场参数下的流体力系数,所述的流体力系数数据库为:<img file="dest_path_image001.GIF" wi="507" he="347" />表中:<i>U</i>为横向流速,<i><img file="dest_path_image003.GIF" wi="26" he="24" /></i>为脉动阻力系数,<i><img file="dest_path_image005.GIF" wi="23" he="21" /></i>为脉动升力系数,<img file="dest_path_image007.GIF" wi="29" he="31" />为稳态阻力系数,<i>S<sub>t</sub></i>为Strouhal频率;S4、基于所述流场参数和所述传热管的结构参数,判断所述流场参数是否属于第二预设范围,若所述流场参数属于所述第二预设范围,则利用公式1~公式4计算换热器传热管的漩涡脱落诱发振动响应;S5、基于所述流场参数和所述传热管的结构参数,判断所述流场参数是否属于第三预设范围,若所述流场参数属于所述第三预设范围,则利用公式1、公式2、公式5、公式6计算换热器传热管的漩涡脱落诱发振动响应;其中,公式1~公式6分别指:公式1:<img file="dest_path_image009.GIF" wi="451" he="90" />;公式2:<img file="dest_path_image011.GIF" wi="428" he="110" />;公式3:<img file="dest_path_image013.GIF" wi="367" he="61" />;公式4:<img file="dest_path_image015.GIF" wi="367" he="61" />;公式5:<img file="dest_path_image017.GIF" wi="342" he="60" />;公式6:<img file="dest_path_image019.GIF" wi="346" he="58" />;式中,<img file="dest_path_image021.GIF" wi="129" he="48" />,<img file="dest_path_image023.GIF" wi="104" he="46" />,<i><img file="dest_path_image025.GIF" wi="93" he="25" /></i>,<i><img file="dest_path_image027.GIF" wi="87" he="24" /></i>;<i>σ<sub>x</sub></i>=0.3×4、<i>σ<sub>y</sub></i>=0.3、<i>Q<sub>x</sub></i>=12×4、<i>Q<sub>y</sub></i>=12,<i><img file="dest_path_image028.GIF" wi="95" he="34" /></i>,<i><img file="dest_path_image029.GIF" wi="88" he="29" /></i>;<i><img file="dest_path_image030.GIF" wi="35" he="30" /></i>,<i><img file="dest_path_image031.GIF" wi="31" he="29" /></i>,<img file="71694dest_path_image007.GIF" wi="29" he="31" />利用所述流体力系数数据库通过插值得到;<i>EI</i>为传热管的弯曲刚度,<i>A</i><sub>p</sub>为传热管的横截面积,<i>D</i>为传热管的外径,<i>L</i>为传热管的长度,<i>c</i>为结构阻尼,<i>ρ</i>为流体密度,<i>m</i><sub>f</sub>和<i>m</i><sub>s</sub>为单位长度的传热管内流体的质量和单位长度的传热管的质量;<i>u</i>和<i>w</i>为传热管的流向振动位移和横向振动位移,<i>ω</i><sub>L</sub>和<i>ω</i><sub>D</sub>为升力方向的漩涡脱落频率和阻力方向的漩涡脱落频率,其中,<i>ω</i><sub>D</sub>=2<i>ω</i><sub>L</sub>,<i>ω</i><sub>L</sub>=2π<i>US<sub>t</sub></i>/<i>D</i>,<i>S<sub>t</sub></i>利用所述流体力系数数据库通过插值得到,<i>U</i>为横向流速,<i>V</i>为传热管内流体的速度;<i>x</i>为传热管的流向振动方向,<i>y</i>为传热管的横向振动方向,<i>z</i>为传热管初始横向挠度为零时的轴向中心线,<i>t</i>为时间。
地址 610000 四川省成都市一环路南三段28号