CN113190927A - 一种蒸汽发生器u型传热管流弹激励磨损分析计算方法 - Google Patents

一种蒸汽发生器u型传热管流弹激励磨损分析计算方法 Download PDF

Info

Publication number
CN113190927A
CN113190927A CN202110460323.9A CN202110460323A CN113190927A CN 113190927 A CN113190927 A CN 113190927A CN 202110460323 A CN202110460323 A CN 202110460323A CN 113190927 A CN113190927 A CN 113190927A
Authority
CN
China
Prior art keywords
heat transfer
transfer pipe
analysis
wear
fluid
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110460323.9A
Other languages
English (en)
Other versions
CN113190927B (zh
Inventor
张锴
唐力晨
李晨
张伟
刘畅
矫明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Nuclear Engineering Research and Design Institute Co Ltd
Original Assignee
Shanghai Nuclear Engineering Research and Design Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Nuclear Engineering Research and Design Institute Co Ltd filed Critical Shanghai Nuclear Engineering Research and Design Institute Co Ltd
Priority to CN202110460323.9A priority Critical patent/CN113190927B/zh
Publication of CN113190927A publication Critical patent/CN113190927A/zh
Application granted granted Critical
Publication of CN113190927B publication Critical patent/CN113190927B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明的目的在于公开一种蒸汽发生器U型传热管流弹激励磨损分析计算方法,与现有技术相比,将流体弹性刚度引入到传热管磨损分析中,能够在流速接近但未达到流弹性失稳时,较为准确的预测传热管的磨损情况,在设计阶段及运行服务阶段用于蒸汽发生器设计及寿命管理,实现本发明的目的。

Description

一种蒸汽发生器U型传热管流弹激励磨损分析计算方法
技术领域
本发明涉及一种磨损分析计算方法,特别涉及一种蒸汽发生器U型传热管流弹激励磨损分析计算方法。
背景技术
流致振动和磨损是核电站蒸汽发生器传热管失效的主要因素。在设计阶段,准确预测蒸汽发生器传热管在寿期内的流致振动和磨损情况非常重要,是保证蒸汽发生器寿期内正常工作的重要保障。对于三代核电站蒸汽发生器而言,蒸汽发生器传热管的设计寿命应该与设备本身同寿,即60年。
在传统计算分析中,对于两相流中的管束流弹失稳问题已有大量研究,其中以H.J.Connors在相关领域的工作最有代表性(Fluid-elastic vibration of tube arraysexcited by cross-flow.Flow Induced Vibration of Heat Exchangers,pages 42-56,1970),该论文给出了直管在流弹失稳的一般分析方法。对于管束中的一根管而言,在Connors公式的运用上需要一些改进,这在ASME规范附录N中有具体介绍。而在传热管磨损预测中,采用了流弹失稳系数作为输入,一般认为在未出现流弹失稳时,传热管的磨损可忽略不计。
在实际电站运行经验中发现,部分传热管在流速接近但未达到流弹失稳流速时,传热管出现了较为明显的磨损现象。说明传统磨损预测方法不能准确预测磨损现象的发生。
因此,特别需要一种蒸汽发生器U型传热管流弹激励磨损分析计算方法,以解决上述现有存在的问题。
发明内容
本发明的目的在于提供一种蒸汽发生器U型传热管流弹激励磨损分析计算方法,针对现有技术的不足,将传热管流致振动和微动磨损的计算分析耦合起来,判断考虑管壁减薄抗振条支撑变弱条件下是否会出现流弹失稳现象。同时能够预测流致振动情况变化条件下长时间内的传热管磨损情况。
本发明所解决的技术问题可以采用以下技术方案来实现:
一种蒸汽发生器U型传热管流弹激励磨损分析计算方法,其特征在于,包括如下步骤:
(1)支撑刚度测定试验:
通过支撑刚度测定试验,在不同的间隙或夹持力条件下,测定其结构频率,通过数值拟合,获得其接触条件和支撑刚度K的函数关系,并以此作为后续磨损分析模型的边界条件;
(2)流体弹性刚度计算:
采用以下公式计算流体弹性刚度,
Figure BDA0003042012270000021
其中,D为传热管直径,L为传热管无支撑跨距,ξ为流体载荷系数,ρ为流体密度,v0为L范围内的模态加权管间流速,δ为传热管与抗振条或支撑板的间隙;
(3)传热管随机激励分析和传热管磨损分析:
采用流体功率谱密度,对传热管结构开展随机激励分析,在支撑板和传热管接触点、抗振条和传热管接触点上设置弹簧边界条件,弹簧刚度为K-K’;计算获得支撑点的支撑反力和支撑点的相对位移的RSM值;根据Archard公式计算磨损。
本发明的蒸汽发生器U型传热管流弹激励磨损分析计算方法,与现有技术相比,将流体弹性刚度引入到传热管磨损分析中,能够在流速接近但未达到流弹性失稳时,较为准确的预测传热管的磨损情况,在设计阶段及运行服务阶段用于蒸汽发生器设计及寿命管理,实现本发明的目的。
本发明的特点可参阅本案图式及以下较好实施方式的详细说明而获得清楚地了解。
附图说明
图1为本发明的蒸汽发生器U型传热管流弹激励磨损分析计算方法的流程示意图;
图2为本发明的U型传热管有限元分析模型的结构示意图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体图示,进一步阐述本发明。
如图1和图2所示,本发明的蒸汽发生器U型传热管流弹激励磨损分析计算方法,包括支撑刚度测定试验、流体弹性刚度计算、传热管随机激励分析和传热管磨损分析四个部分,分析流程如图1所示。
支撑刚度测定试验:
在传统分析中,一般认为抗振条和支撑板对传热管的支撑条件可近似假设为简支边界条件。但实际结构中,抗振条、支撑板与传热管之间存在设计间隙和装配公差,使得其实际接触提供的边界条件弱于简支边界条件。可以通过支撑刚度测定试验,在不同的间隙或夹持力条件下,测定其结构频率。通过数值拟合,获得其接触条件和支撑刚度K的函数关系。并以此作为后续磨损分析模型的边界条件。
流体弹性刚度计算:
采用以下公式计算流体弹性刚度,
Figure BDA0003042012270000031
其中,D为传热管直径,L为传热管无支撑跨距,ξ为流体载荷系数,ρ为流体密度,v0为L范围内的模态加权管间流速,δ为传热管与抗振条或支撑板的间隙。
传热管随机激励分析和传热管磨损分析:
采用流体功率谱密度,对传热管结构开展随机激励分析,在支撑板和传热管接触点、抗振条和传热管接触点上设置弹簧边界条件,弹簧刚度为K-K’。计算获得支撑点的支撑反力和支撑点的相对位移的RSM值。根据Archard公式计算磨损。
实施例:
在核电厂蒸汽发生器中,一般有几千根至上万根的传热管。针对其中任意一根传热管,建立如图2所示的有限元模型。其中传热管与管板连接处简化为固支边界条件;与支撑板接触处简化为支撑板接触弹簧边界条件;与抗振条接触处简化为抗振条接触弹簧边界条件。弹簧边界条件的刚度设置为K-K’,其中K由结构决定,K’由流场参数决定。采用流体功率谱密度,对传热管结构开展随机激励分析,计算获得支撑点的支撑反力和支撑点的相对位移的RSM值。根据Archard公式计算磨损。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内,本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (1)

1.一种蒸汽发生器U型传热管流弹激励磨损分析计算方法,其特征在于,包括如下步骤:
(1)支撑刚度测定试验:
通过支撑刚度测定试验,在不同的间隙或夹持力条件下,测定其结构频率,通过数值拟合,获得其接触条件和支撑刚度K的函数关系,并以此作为后续磨损分析模型的边界条件;
(2)流体弹性刚度计算:
采用以下公式计算流体弹性刚度,
Figure FDA0003042012260000011
其中,D为传热管直径,L为传热管无支撑跨距,ξ为流体载荷系数,ρ为流体密度,v0为L范围内的模态加权管间流速,δ为传热管与抗振条或支撑板的间隙;
(3)传热管随机激励分析和传热管磨损分析:
采用流体功率谱密度,对传热管结构开展随机激励分析,在支撑板和传热管接触点、抗振条和传热管接触点上设置弹簧边界条件,弹簧刚度为K-K’;计算获得支撑点的支撑反力和支撑点的相对位移的RSM值;根据Archard公式计算磨损。
CN202110460323.9A 2021-04-27 2021-04-27 一种蒸汽发生器u型传热管流弹激励磨损分析计算方法 Active CN113190927B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110460323.9A CN113190927B (zh) 2021-04-27 2021-04-27 一种蒸汽发生器u型传热管流弹激励磨损分析计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110460323.9A CN113190927B (zh) 2021-04-27 2021-04-27 一种蒸汽发生器u型传热管流弹激励磨损分析计算方法

Publications (2)

Publication Number Publication Date
CN113190927A true CN113190927A (zh) 2021-07-30
CN113190927B CN113190927B (zh) 2023-08-04

Family

ID=76979482

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110460323.9A Active CN113190927B (zh) 2021-04-27 2021-04-27 一种蒸汽发生器u型传热管流弹激励磨损分析计算方法

Country Status (1)

Country Link
CN (1) CN113190927B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107025315A (zh) * 2016-02-02 2017-08-08 上海核工程研究设计院 一种核电站蒸汽发生器u型传热管流致振动及微动磨损耦合分析计算方法
WO2019223282A1 (zh) * 2018-05-25 2019-11-28 深圳中广核工程设计有限公司 核电厂蒸汽发生器u形传热管用支撑板及其蒸汽发生器
CN111859752A (zh) * 2020-07-18 2020-10-30 西安交通大学 一种核反应堆蒸汽发生器流致振动计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107025315A (zh) * 2016-02-02 2017-08-08 上海核工程研究设计院 一种核电站蒸汽发生器u型传热管流致振动及微动磨损耦合分析计算方法
WO2019223282A1 (zh) * 2018-05-25 2019-11-28 深圳中广核工程设计有限公司 核电厂蒸汽发生器u形传热管用支撑板及其蒸汽发生器
CN111859752A (zh) * 2020-07-18 2020-10-30 西安交通大学 一种核反应堆蒸汽发生器流致振动计算方法

Also Published As

Publication number Publication date
CN113190927B (zh) 2023-08-04

Similar Documents

Publication Publication Date Title
CN107025315B (zh) 核电站u型传热管流致振动及微动磨损耦合分析计算方法
Green et al. PWR steam generators
Ko Heat exchanger tube fretting wear: review and application to design
CN109085814A (zh) 一种火电汽轮机组整体设备系统延寿评估方法
CN110909505B (zh) 核电厂疲劳监测和寿命评估系统的瞬态温度场计算方法
Hassan et al. The effects of support geometry on the turbulence response of loosely supported heat exchanger tubes
CN111931319B (zh) 非线性支承管束在横向流中的振动特性分析方法
CN113190927A (zh) 一种蒸汽发生器u型传热管流弹激励磨损分析计算方法
Au-Yang Flow-induced wear in steam generator tubes—Prediction versus operational experience
CN113011111A (zh) 一种堆内构件棒状结构湍流激振评估方法及装置
Park et al. Wear depth prediction for steam generator tubes using an estimated in situ wear coefficient
Jo et al. Fluidelastic instability analysis of operating nuclear steam generator U-tubes
CN111861011B (zh) 一种超临界压力流体对流换热性能预测方法及系统
Jendrzejczyk Dynamic characteristics of heat exchanger tubes vibrating in a tube support plate inactive mode
Kim et al. Structural and vibration analysis considering the flow velocity of the heat exchanger
CN113673130A (zh) 蒸发器管束两相流功率谱密度相关长度的获取方法及系统
Hassan et al. A numerical characterization of flow-induced vibration and fretting wear potential in nuclear steam generator tube bundles
Cai et al. A theory for fluidelastic instability of tube-support-plate-inactive modes
Low Fretting problems and some solutions in power plant machinery
CN111855799B (zh) 一种横向流作用下非线性约束管束的损伤识别方法
Kreider et al. Thermal performance degradation and heat-transfer fouling
Keogh et al. Streamwise Fluidelastic Instability in a Deformed Rotated Triangular Tube Array With Pitch-to-Diameter Ratio of 1.375
CN115712998B (zh) 调峰机组金属通流部件内壁开裂的可修复性判定方法
Payen et al. A probabilistic approach for the computation of non-linear vibrations of tubes under cross-flow
Fesenko et al. Forced oscillations of tube bundles in liquid cross-flow

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: No. 29 Hong Cao Road, Xuhui District, Shanghai

Applicant after: Shanghai Nuclear Engineering Research and Design Institute Co.,Ltd.

Address before: No. 29 Hong Cao Road, Xuhui District, Shanghai

Applicant before: SHANGHAI NUCLEAR ENGINEERING RESEARCH & DESIGN INSTITUTE Co.,Ltd.

GR01 Patent grant
GR01 Patent grant