CN113190927A - 一种蒸汽发生器u型传热管流弹激励磨损分析计算方法 - Google Patents
一种蒸汽发生器u型传热管流弹激励磨损分析计算方法 Download PDFInfo
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force 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型传热管流弹激励磨损分析计算方法。
背景技术
流致振动和磨损是核电站蒸汽发生器传热管失效的主要因素。在设计阶段,准确预测蒸汽发生器传热管在寿期内的流致振动和磨损情况非常重要,是保证蒸汽发生器寿期内正常工作的重要保障。对于三代核电站蒸汽发生器而言,蒸汽发生器传热管的设计寿命应该与设备本身同寿,即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)流体弹性刚度计算:
其中,D为传热管直径,L为传热管无支撑跨距,ξ为流体载荷系数,ρ为流体密度,v0为L范围内的模态加权管间流速,δ为传热管与抗振条或支撑板的间隙;
(3)传热管随机激励分析和传热管磨损分析:
采用流体功率谱密度,对传热管结构开展随机激励分析,在支撑板和传热管接触点、抗振条和传热管接触点上设置弹簧边界条件,弹簧刚度为K-K’;计算获得支撑点的支撑反力和支撑点的相对位移的RSM值;根据Archard公式计算磨损。
本发明的蒸汽发生器U型传热管流弹激励磨损分析计算方法,与现有技术相比,将流体弹性刚度引入到传热管磨损分析中,能够在流速接近但未达到流弹性失稳时,较为准确的预测传热管的磨损情况,在设计阶段及运行服务阶段用于蒸汽发生器设计及寿命管理,实现本发明的目的。
本发明的特点可参阅本案图式及以下较好实施方式的详细说明而获得清楚地了解。
附图说明
图1为本发明的蒸汽发生器U型传热管流弹激励磨损分析计算方法的流程示意图;
图2为本发明的U型传热管有限元分析模型的结构示意图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体图示,进一步阐述本发明。
如图1和图2所示,本发明的蒸汽发生器U型传热管流弹激励磨损分析计算方法,包括支撑刚度测定试验、流体弹性刚度计算、传热管随机激励分析和传热管磨损分析四个部分,分析流程如图1所示。
支撑刚度测定试验:
在传统分析中,一般认为抗振条和支撑板对传热管的支撑条件可近似假设为简支边界条件。但实际结构中,抗振条、支撑板与传热管之间存在设计间隙和装配公差,使得其实际接触提供的边界条件弱于简支边界条件。可以通过支撑刚度测定试验,在不同的间隙或夹持力条件下,测定其结构频率。通过数值拟合,获得其接触条件和支撑刚度K的函数关系。并以此作为后续磨损分析模型的边界条件。
流体弹性刚度计算:
其中,D为传热管直径,L为传热管无支撑跨距,ξ为流体载荷系数,ρ为流体密度,v0为L范围内的模态加权管间流速,δ为传热管与抗振条或支撑板的间隙。
传热管随机激励分析和传热管磨损分析:
采用流体功率谱密度,对传热管结构开展随机激励分析,在支撑板和传热管接触点、抗振条和传热管接触点上设置弹簧边界条件,弹簧刚度为K-K’。计算获得支撑点的支撑反力和支撑点的相对位移的RSM值。根据Archard公式计算磨损。
实施例:
在核电厂蒸汽发生器中,一般有几千根至上万根的传热管。针对其中任意一根传热管,建立如图2所示的有限元模型。其中传热管与管板连接处简化为固支边界条件;与支撑板接触处简化为支撑板接触弹簧边界条件;与抗振条接触处简化为抗振条接触弹簧边界条件。弹簧边界条件的刚度设置为K-K’,其中K由结构决定,K’由流场参数决定。采用流体功率谱密度,对传热管结构开展随机激励分析,计算获得支撑点的支撑反力和支撑点的相对位移的RSM值。根据Archard公式计算磨损。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内,本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (1)
1.一种蒸汽发生器U型传热管流弹激励磨损分析计算方法,其特征在于,包括如下步骤:
(1)支撑刚度测定试验:
通过支撑刚度测定试验,在不同的间隙或夹持力条件下,测定其结构频率,通过数值拟合,获得其接触条件和支撑刚度K的函数关系,并以此作为后续磨损分析模型的边界条件;
(2)流体弹性刚度计算:
其中,D为传热管直径,L为传热管无支撑跨距,ξ为流体载荷系数,ρ为流体密度,v0为L范围内的模态加权管间流速,δ为传热管与抗振条或支撑板的间隙;
(3)传热管随机激励分析和传热管磨损分析:
采用流体功率谱密度,对传热管结构开展随机激励分析,在支撑板和传热管接触点、抗振条和传热管接触点上设置弹簧边界条件,弹簧刚度为K-K’;计算获得支撑点的支撑反力和支撑点的相对位移的RSM值;根据Archard公式计算磨损。
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)
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 | 西安交通大学 | 一种核反应堆蒸汽发生器流致振动计算方法 |
-
2021
- 2021-04-27 CN CN202110460323.9A patent/CN113190927B/zh active Active
Patent Citations (3)
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 |