CN116542177A - 一种基于启停工况分析评判的水轮机寿命评估方法及系统 - Google Patents
一种基于启停工况分析评判的水轮机寿命评估方法及系统 Download PDFInfo
- Publication number
- CN116542177A CN116542177A CN202310542959.7A CN202310542959A CN116542177A CN 116542177 A CN116542177 A CN 116542177A CN 202310542959 A CN202310542959 A CN 202310542959A CN 116542177 A CN116542177 A CN 116542177A
- Authority
- CN
- China
- Prior art keywords
- water turbine
- fatigue
- analysis
- finite element
- judgment
- 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.)
- Pending
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 92
- 238000000034 method Methods 0.000 title claims abstract description 84
- 238000004458 analytical method Methods 0.000 title claims abstract description 68
- 238000004088 simulation Methods 0.000 claims abstract description 39
- 230000003068 static effect Effects 0.000 claims abstract description 26
- 239000012530 fluid Substances 0.000 claims abstract description 24
- 239000000463 material Substances 0.000 claims abstract description 20
- 230000007246 mechanism Effects 0.000 claims abstract description 9
- 238000011160 research Methods 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims description 22
- 230000008859 change Effects 0.000 claims description 13
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 3
- 230000006870 function Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 210000000746 body region Anatomy 0.000 claims description 2
- 238000012423 maintenance Methods 0.000 abstract description 7
- 206010016256 fatigue Diseases 0.000 description 63
- 230000035882 stress Effects 0.000 description 26
- 239000000306 component Substances 0.000 description 15
- 230000008569 process Effects 0.000 description 14
- 238000011156 evaluation Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 7
- 230000036541 health Effects 0.000 description 7
- 230000000007 visual effect Effects 0.000 description 7
- 238000004590 computer program Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 239000008358 core component Substances 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000001052 transient effect Effects 0.000 description 3
- 230000004075 alteration Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000013210 evaluation model Methods 0.000 description 2
- 230000000977 initiatory effect Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009960 carding Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
Classifications
-
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/20—Hydro energy
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明属于水轮机寿命评估技术领域,具体提供了一种基于启停工况分析评判的水轮机寿命评估方法及系统,其中方法包括:采用计算流体动力学(CFD)方法,针对水轮机启停工况对水轮机内部场进行三维流场仿真,以确定结构危险部位;针对机构危险部位的三维模型进行有限元网格划分以进行有限元静力分析;基于名义应力法,输入有限元静力分析结果、材料疲劳性能参数及疲劳载荷,利用疲劳分析软件ncode对结构危险部位的疲劳寿命进行仿真研究。通过结构力学与疲劳寿命分析方法对水轮机关键部件进行疲劳寿命预测与风险预警,辅助电站运维采取提前维护与更换的决策,能够有效降低核心组件因疲劳发生意外失效的概率。
Description
技术领域
本发明涉及水轮机寿命评估技术领域,更具体地,涉及一种基于启停工况分析评判的水轮机寿命评估方法及系统。
背景技术
我国在六十年代至八十年代建有大量的中小型水电站,目前多数已接近或达到设计寿命,除少量报废外,大部分仍在使用,设备老化问题突出,延长使用寿命依据不足,风险性很大。如何科学地分析此类电站延续运行的安全性和可行性,准确评估各电站剩余运行寿命,并提出合理的建议,保证机组安全运行就显得十分迫切。水电站设备管理的主要问题在于主设备不宜更换,传统的解决方案一般是对其进行评估检测,主要采用现场评估检查方法,现场评估检查包括查阅资料和实地检查,现场评估结束后汇总到综合寿命评估;或采用现场评估检查和试验检测相结合的综合寿命评估法,根据现场检查评估结果,对需要和必须开展试验检测的设备开展例行和专项试验检测工作,获取各设备设施的性能状况,推算各设备设施的寿命。
但上述两种方法不能对设备运行全流程进行寿命评估和管理,且主要存在以下两个问题:一是水电站的传统寿命评估方法存在着无法对水轮机进行启停工况状态预测和关键部件疲劳风险预报等问题;二是水电站的传统寿命评估方法不具备使水轮机健康状态从数据分析到寿命预测、运维决策的完整评估过程。
发明内容
本发明针对现有技术中存在的一是水电站的传统寿命评估方法存在着无法对水轮机进行启停工况状态预测和关键部件疲劳风险预报等问题;二是水电站的传统寿命评估方法不具备使水轮机健康状态从数据分析到寿命预测、运维决策的完整评估过程的技术问题。
本发明提供了一种基于启停工况分析评判的水轮机寿命评估方法,包括以下步骤:
S1,采用计算流体动力学(CFD)方法,针对水轮机启停工况对水轮机内部场进行三维流场仿真,以确定结构危险部位;
S2,针对机构危险部位的三维模型进行有限元网格划分以进行有限元静力分析;
S3,基于名义应力法,输入有限元静力分析结果、材料疲劳性能参数及疲劳载荷,利用疲劳分析软件ncode对结构危险部位的疲劳寿命进行仿真研究。
优选地,所述S1具体包括:
S11,用有限个离散点上的变量值的集合来代替原始连续的物理量,对流动基本方程进行离散,建立离散点上的变量值之间关系的代数方程组;
S12,通过求解代数方程组获得变量的近似值;
S13,通过数值计算方法,定量地描述流场,并用图像显示的方法展示三维流场内部的变化。
优选地,所述S2具体包括:
S21,将水轮机的结构危险部位的水体区域离散成一组控制体;
S22,在这组控制体上对质量守恒方程组和动量守恒方程组求解,将偏微分方程组离散为代数方程组得到离散方程;
S23,通过一维数值仿真计算迭代求解离散方程,使最后结果逐步贴近流动参数的真实值直到指定精度。
优选地,所述S22具体包括:
S221,将计算区域划分为多个不重复的控制体积,并使每个控制体积都有一个网格节点作为代表,将待求的守恒型微分方程在任一控制体积及一定时间间隔内对空间与时间作积;
S222,对待求解函数及导数对时间及空间的变化型线或插值方式作出假设;
S223,对选定的型线作出积分并整理成一组关于节点上未知量的离散方程。
优选地,所述S23具体包括:
得离散方程后,基于SIMPLE(Semi Implicit Method for Pressure LinkedEquations)算法,通过压力场求解速度场,并对初始压力场进行修正处理,保证迭代计算的连续性。
优选地,所述S23具体包括:根据需要计算的模型假设一个速度场,以便于完善离散方程;
假定压力场的压力为p*,根据p*处理离散方程得到定义为u*、v*的速度;
根据SIMPLE算法以达到满足迭代的联系性要求进行修正,修正后的压力值为p′,速度值为u′及v′;
将(p′+p*)与(u′+u*)、(v′+v*)作为下一步迭代计算的初始量继续进行迭代运算。
优选地,所述S3具体包括:
S31,基于有限元静力分析结果,得到结构危险部位在外载荷作用下的应力分布及危险区域;
S32,输入材料的疲劳性能参数,包括材料的S-N曲线,ε-N曲线或σ-N曲线;
S33,加载疲劳载荷谱,对水轮机的结构危险部位的疲劳寿命进行分析计算,得到危险区域的疲劳寿命。
本发明还提供了一种基于启停工况分析评判的水轮机寿命评估系统,所述系统用于实现基于启停工况分析评判的水轮机寿命评估方法,包括:
流场仿真模块,用于采用计算流体动力学(CFD)方法,针对水轮机启停工况对水轮机内部场进行三维流场仿真,以确定结构危险部位;
有限元静力分析模块,用于针对机构危险部位的三维模型进行有限元网格划分以进行有限元静力分析;
疲劳分析模块,用于基于名义应力法,输入有限元静力分析结果、材料疲劳性能参数及疲劳载荷,利用疲劳分析软件ncode对结构危险部位的疲劳寿命进行仿真研究。
本发明还提供了一种电子设备,包括存储器、处理器,所述处理器用于执行存储器中存储的计算机管理类程序时实现基于启停工况分析评判的水轮机寿命评估方法的步骤。
本发明还提供了一种计算机可读存储介质,其上存储有计算机管理类程序,所述计算机管理类程序被处理器执行时实现基于启停工况分析评判的水轮机寿命评估方法的步骤。
有益效果:本发明提供的一种基于启停工况分析评判的水轮机寿命评估方法及系统,其中方法包括:采用计算流体动力学(CFD)方法,针对水轮机启停工况对水轮机内部场进行三维流场仿真,以确定结构危险部位;针对机构危险部位的三维模型进行有限元网格划分以进行有限元静力分析;基于名义应力法,输入有限元静力分析结果、材料疲劳性能参数及疲劳载荷,利用疲劳分析软件ncode对结构危险部位的疲劳寿命进行仿真研究。通过结构力学与疲劳寿命分析方法对水轮机关键部件进行疲劳寿命预测与风险预警,辅助电站运维采取提前维护与更换的决策,能够有效降低核心组件因疲劳发生意外失效的概率。
附图说明
图1为本发明提供的一种基于启停工况分析评判的水轮机寿命评估方法流程图;
图2为本发明提供的CFD仿真流程示意图;
图3为本发明提供的名义应力法估算构件疲劳寿命的流程图;
图4为本发明提供的典型的S-N曲线图;
图5为本发明提供的ncode疲劳寿命分析流程图;
图6为本发明提供的一种可能的电子设备的硬件结构示意图;
图7为本发明提供的一种可能的计算机可读存储介质的硬件结构示意图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
如图1-图5所示,本发明通过对电站水轮机开展数字孪生建设,结合数值仿真技术,拓展监测数据分析及应用,支撑水轮机状态数字化映射。建立一体化集成状态分析及评估模型,基于过程监测数据,进行价值挖掘,作为寿命评估的支撑数据,保障评估模型构建的真实性与可靠性;
借助CFD软件通过数值仿真技术手段,对整个水轮机全流道1:1的三维数值仿真模型,获取水轮机内部流场及关键部位材料压力载荷信息。基于各类典型工况,在水轮机启动、停机等瞬态过程变化中,模拟出导叶开度变化过程以及流场的变化状态,为水轮机关键部件的寿命评估提供必须的载荷输入。
以水轮机关键部件的图纸以及其它资料为依据,得到水轮机过流部件、水轮机转轮1:1的三维数值仿真模型,实现结构仿真结果的可视化展示,同时作为支撑后续关键部件健康状态评估的数据基础;
利用有限元断裂力学计算结果和合理、可靠的概率分析方法,对水轮机关键部件疲劳裂纹的萌生进行预测,对关键部件当前风险进行预警并提供运维建议。
通过开发水轮机内部流场动态展示、核心部件受力分析等应用。确定关键寿命部件,对部件内部流场和力场进行仿真,梳理反映水电机组典型特征的代表工况,全面覆盖水轮机运行的全周期过程,结合流体仿真结果,针对不同部件的典型问题开展强度计算、模态分析、谐振分析和裂纹分析,全面了解水轮机关键部件的承载状况,从而全面掌握机组的运行状态及关键部件的健康情况。结合水轮机的在线监测数据和实际的有限元结构仿真数据,对启停工况下核心部件(结构危险部位)的剩余寿命评估,对关键部件可能发生的疲劳风险进行及时预警,科学客观评价设备健康状况。
具体方法流程如下:
S1,水轮机内流场仿真。
水轮机启停工况是一个不稳定的瞬态过程,启动和停机时的水力冲击和水压脉动都会对水轮机的部件产生不利影响。由于水轮机内部流场复杂,且无法在水轮机内部安装传感器实现对关键部件的监测,为了获取水轮机内部流场及关键部位材料压力载荷信息。借助计算流体动力学CFD软件通过数值仿真技术手段,针对典型启停工况,复现水轮机内复杂的瞬态流场动态变化过程,全面掌握流体关键参数变化趋势,并实现动态流体的高清晰效果可视化展示,同时获得水轮机整体的受力载荷参数信息。CFD仿真流程图如图2所示。
具体的方案,采用计算流体动力学方法对水轮机内部场进行三维流场仿真。基本思想是用一系列有限个离散点上的变量值的集合来代替原始连续的物理量,通过特定的原则对流动基本方程进行离散,建立起离散点上的变量值之间关系的代数方程组,然后通过求解代数方程组获得变量的近似值。最后通过数值计算方法,定量地描述流场,并用图像显示的方法展示三维流场内部的变化。这样便可以得知水轮机的结构危险部位(即关键部件)以及该危险部位的受力载荷参数信息。
流场仿真具体包括以下几点:
(1)控制方程。CFD仿真计算的过程本质上是求解控制方程。质量守恒定律和动量守恒定律是描述流体运动的基本定律。对于本项目所研究的水轮机内部流动问题,流动介质是低速的水流,可近似看成不可压缩的流体,因此在求解方程时只需对连续性方程和动量方程进行求解。
(2)连续性方程。
其中,t为单位时间,V为控制体体积,p为流体密度,S为控制面积,n为控制面面积,n为微元面积矢量dS法外线的单位向量,U为微元面积dS上的流体速度。
在直角坐标系中可以改写为:
式中,下标i可取值为1,2,3,以表示三个空间坐标。
(3)动量方程。动量守恒定律在流场中的数学表达称为动量方程,即作用在控制体上的外力的合力与每秒钟内通过控制面流入的动量之和,等于每秒钟内控制体中流体动量的增量。
其中,v为控制体体积,p为流体密度,S为控制速度,n为微元面积矢量ds外法线的单位向量,U为微元表面ds上的流体速度,П为微元面积矢量ds的应力张量。
在直角坐标系中可以改写为:
其中,σij为微元面积矢量dS的应力张量П的分量,i、j可取值为1,2,3,以表示三个空间坐标。
(4)控制方程的求解方法。求解控制方程的数值方法有有限差分法、有限元法和有限体积法。本发明拟采用有限体积法对控制方程进行求解,首先将水轮机关键部件的水体区域离散成一组控制体;然后在这组控制体上对质量守恒和动量守恒方程组求解,将偏微分方程组离散为代数方程组;最后通过一维数值仿真计算迭代求解方程组,使最后结果逐步贴近流动参数的真实值。
以水电站水轮机设计图纸以及其它资料为依据,利用Creo 3D建模软件完成对整个水轮机全流道1:1的三维数值仿真模型,采用ICEM构建水轮机的结构化或非结构化的高质量流体仿真网格模型,为全流道可视化仿真提供数字化映射载体。根据水轮机设计工况和运行实际情况,选取具有代表性的若干平稳运行工况、水轮机启动、停止过程变化工况、实际情况,选取具有代表性的若干平稳运行工况、水轮机启动、停止过程变化工况、以及飞逸工况进行流场分析。基于各类典型工况,在水轮机启动、停机等瞬态过程变化中,能够模拟出导叶开度变化过程以及流场的变化状态。同时还可以对整个水轮机内部流场的多维度信息的可视化展示,为水轮机关键部件的材料性能分析与寿命评估提供必须的载荷输入。
S2,核心部件结构仿真。
因水轮机机组在启停过程中常处在振动区内,导致关键部件结构(即结构危险部位)出现共振现象,此外,由调峰、调频引起的功率变化也可能导致大幅度低频动应力。振动和动应力问题可能导致转轮、顶盖等关键部件出现疲劳裂纹问题,因此针对启停工况,对水轮机关键部件的健康状态及疲劳寿命情况进行评估,保障机组安全平稳运行。
以水电站水轮机关键部件的图纸以及其它资料为依据,得到水轮机过流部件、水轮机转轮1:1的三维数值仿真模型。同时根据三维模型的几何特征,进行结构化或非结构化的高质量的有限元网格划分。根据流场仿真结果,结合水轮机实际运行监测数据,完成水轮机关键部件的结构力学有限元静力分析,获取关键部件的模态、位移、应力结果,完成结构仿真结果的可视化展示,同时作为支撑后续关键部件健康状态评估的数据基础。
利用有限元断裂力学计算结果和合理、可靠的概率分析方法,对水轮机关键部件疲劳裂纹的萌生进行预测。获得实际运行工况下的水轮机关键部件的受力状态与断裂损伤评估结果,以可视化模型及三维云图的形式展示水轮机关键部件的应力分布情况及易出现断裂的高风险区域,对关键部件当前风险进行预警。
计算流体力学,基本思路可分为如下三点:
首先,将计算区域划分为一系列不重复的控制体积,并使每个控制体积都有一个网格节点作为代表,将待求的守恒型微分方程在任一控制体积及一定时间间隔内对空间与时间作积。
然后,对待求解函数及导数对时间及空间的变化型线或插值方式作出假设;
最后,对选定的型线作出积分并整理成一组关于节点上未知量的离散方程。
具体操作为:
通过网格体积控制值对流场计算模型进行网格划分处理;离散方程的获得可借助控制方程对不同网格体积积分;根据求解离散方程得到各单元节点处的变量值,经插值处理可获得此控制体积的截面变量值。
得到离散的代数方程组(即离散方程)后,需要对其进行处理。为此,有限体积法为有限元求解问题提供了解决方案,通过选择合适的算法对离散方程进行处理以更加准确便捷地解决实际工程问题。SIMPLE(Semi Implicit Method for Pressure LinkedEquations)算法是目前流体计算领域应用最为普遍的数值计算方法。SIMPLE算法通过压力场求解速度场,但在实际问题中求解出的速度场通常难以满足迭代要求,此时为了保证迭代计算的连续性,则需要对初始压力场进行修正处理。如此反复直到获得收敛解,该算法的基本步骤为:
根据需要计算的模型假设一个速度场,以便于完善离散方程;
假定压力场的压力为p;根据p处理动量离散方程得到定义为u、v*的速度;
根据前面速度对压力场进行修正计算方法,达到满足迭代的联系性要求进行修正,修正后的压力值为p′,速度值为u′及v′;
将(p′+p*)与(u′+u*)、(v′+v*)作为下一步迭代计算的初始量继续进行迭代运算。
S3,疲劳寿命分析。
基于名义应力法,利用疲劳分析软件ncode对水轮机关键部件的疲劳寿命进行仿真研究。
基于名义应力法估算结构疲劳寿命的流程,名义应力法以材料的S-N曲线为基础,对照结构危险部位的名义应力和应力集中系数,并结合疲劳损伤累积理论,对结构疲劳寿命进行计算,流程图如图4所示。
利用名义应力法估算结构的疲劳寿命时,需要建立材料的S-N曲线,一条完整的S-N曲线如图4所示,可分为3段:低周疲劳区(LCF),高周疲劳区(HCF)和亚疲劳区(SF)。N=1/4时,即静拉伸所对应的疲劳强度Smax=Sb,N=106~107时,对应的疲劳强度为疲劳极限,即当循环应力低于此疲劳极限时,结构接近于无限寿命。
搭建疲劳寿命仿真流程。利用有限元方法分析结构的疲劳寿命,采用Ncode软件进行疲劳寿命分析,流程图如图5所示。
利用ncode软件分析结构疲劳寿命时的三个基本条件是:
1.有限元静力分析结果,得到结构危险部位在外载荷作用下的应力分布及危险区域。
2.结构材料的疲劳性能参数,包括材料的S-N曲线,ε-N曲线或σ-N曲线等。
3.疲劳载荷谱。基于以上疲劳寿命分析的基本条件,再设置合适的相关参数就可以对水轮机关键部件的疲劳寿命进行分析计算,得到危险区域的疲劳寿命。而水轮机的结构危险部位的危险区域的疲劳寿命就是水轮机的寿命。
本发明实施例还提供了一种基于启停工况分析评判的水轮机寿命评估系统,所述系统用于实现基于启停工况分析评判的水轮机寿命评估方法,包括:
流场仿真模块,用于采用计算流体动力学(CFD)方法,针对水轮机启停工况对水轮机内部场进行三维流场仿真,以确定结构危险部位;
有限元静力分析模块,用于针对机构危险部位的三维模型进行有限元网格划分以进行有限元静力分析;
疲劳分析模块,用于基于名义应力法,输入有限元静力分析结果、材料疲劳性能参数及疲劳载荷,利用疲劳分析软件ncode对结构危险部位的疲劳寿命进行仿真研究。
请参阅图6为本发明实施例提供的电子设备的实施例示意图。如图6所示,本发明实施例提了一种电子设备,包括存储器1310、处理器1320及存储在存储器1310上并可在处理器1320上运行的计算机程序1311,处理器1320执行计算机程序1311时实现以下步骤:S1,采用计算流体动力学(CFD)方法,针对水轮机启停工况对水轮机内部场进行三维流场仿真,以确定结构危险部位;
S2,针对机构危险部位的三维模型进行有限元网格划分以进行有限元静力分析;
S3,基于名义应力法,输入有限元静力分析结果、材料疲劳性能参数及疲劳载荷,利用疲劳分析软件ncode对结构危险部位的疲劳寿命进行仿真研究。
请参阅图7为本发明提供的一种计算机可读存储介质的实施例示意图。如图7所示,本实施例提供了一种计算机可读存储介质1400,其上存储有计算机程序1411,该计算机程序1411被处理器执行时实现如下步骤:S1,采用计算流体动力学(CFD)方法,针对水轮机启停工况对水轮机内部场进行三维流场仿真,以确定结构危险部位;
S2,针对机构危险部位的三维模型进行有限元网格划分以进行有限元静力分析;
S3,基于名义应力法,输入有限元静力分析结果、材料疲劳性能参数及疲劳载荷,利用疲劳分析软件ncode对结构危险部位的疲劳寿命进行仿真研究。
需要说明的是,在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详细描述的部分,可以参见其它实施例的相关描述。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包括这些改动和变型在内。
Claims (10)
1.一种基于启停工况分析评判的水轮机寿命评估方法,其特征在于,包括以下步骤:
S1,采用计算流体动力学(CFD)方法,针对水轮机启停工况对水轮机内部场进行三维流场仿真,以确定结构危险部位;
S2,针对机构危险部位的三维模型进行有限元网格划分以进行有限元静力分析;
S3,基于名义应力法,输入有限元静力分析结果、材料疲劳性能参数及疲劳载荷,利用疲劳分析软件ncode对结构危险部位的疲劳寿命进行仿真研究。
2.根据权利要求1所述的基于启停工况分析评判的水轮机寿命评估方法,其特征在于,所述S1具体包括:
S11,用有限个离散点上的变量值的集合来代替原始连续的物理量,对流动基本方程进行离散,建立离散点上的变量值之间关系的代数方程组;
S12,通过求解代数方程组获得变量的近似值;
S13,通过数值计算方法,定量地描述流场,并用图像显示的方法展示三维流场内部的变化。
3.根据权利要求1所述的基于启停工况分析评判的水轮机寿命评估方法,其特征在于,所述S2具体包括:
S21,将水轮机的结构危险部位的水体区域离散成一组控制体;
S22,在这组控制体上对质量守恒方程组和动量守恒方程组求解,将偏微分方程组离散为代数方程组得到离散方程;
S23,通过一维数值仿真计算迭代求解离散方程,使最后结果逐步贴近流动参数的真实值直到指定精度。
4.根据权利要求3所述的基于启停工况分析评判的水轮机寿命评估方法,其特征在于,所述S22具体包括:
S221,将计算区域划分为多个不重复的控制体积,并使每个控制体积都有一个网格节点作为代表,将待求的守恒型微分方程在任一控制体积及一定时间间隔内对空间与时间作积;
S222,对待求解函数及导数对时间及空间的变化型线或插值方式作出假设;
S223,对选定的型线作出积分并整理成一组关于节点上未知量的离散方程。
5.根据权利要求4所述的基于启停工况分析评判的水轮机寿命评估方法,其特征在于,所述S23具体包括:
得离散方程后,基于SIMPLE(Semi Implicit Method for Pressure LinkedEquations)算法,通过压力场求解速度场,并对初始压力场进行修正处理,保证迭代计算的连续性。
6.根据权利要求5所述的基于启停工况分析评判的水轮机寿命评估方法,其特征在于,所述S23具体包括:根据需要计算的模型假设一个速度场,以便于完善离散方程;
假定压力场的压力为p*,根据p*处理离散方程得到定义为u*、v*的速度;
根据SIMPLE算法以达到满足迭代的联系性要求进行修正,修正后的压力值为p′,速度值为u′及v′;
将(p′+p*)与(u′+u*)、(v′+v*)作为下一步迭代计算的初始量继续进行迭代运算。
7.根据权利要求1所述的基于启停工况分析评判的水轮机寿命评估方法,其特征在于,所述S3具体包括:
S31,基于有限元静力分析结果,得到结构危险部位在外载荷作用下的应力分布及危险区域;
S32,输入材料的疲劳性能参数,包括材料的S-N曲线,ε-N曲线或σ-N曲线;
S33,加载疲劳载荷谱,对水轮机的结构危险部位的疲劳寿命进行分析计算,得到危险区域的疲劳寿命。
8.一种基于启停工况分析评判的水轮机寿命评估系统,其特征在于,所述系统用于实现如权利要求1-7任一项所述的基于启停工况分析评判的水轮机寿命评估方法,包括:
流场仿真模块,用于采用计算流体动力学(CFD)方法,针对水轮机启停工况对水轮机内部场进行三维流场仿真,以确定结构危险部位;
有限元静力分析模块,用于针对机构危险部位的三维模型进行有限元网格划分以进行有限元静力分析;
疲劳分析模块,用于基于名义应力法,输入有限元静力分析结果、材料疲劳性能参数及疲劳载荷,利用疲劳分析软件ncode对结构危险部位的疲劳寿命进行仿真研究。
9.一种电子设备,其特征在于,包括存储器、处理器,所述处理器用于执行存储器中存储的计算机管理类程序时实现如权利要求1-7任一项所述的基于启停工况分析评判的水轮机寿命评估方法的步骤。
10.一种计算机可读存储介质,其特征在于,其上存储有计算机管理类程序,所述计算机管理类程序被处理器执行时实现如权利要求1-7任一项所述的基于启停工况分析评判的水轮机寿命评估方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310542959.7A CN116542177A (zh) | 2023-05-12 | 2023-05-12 | 一种基于启停工况分析评判的水轮机寿命评估方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310542959.7A CN116542177A (zh) | 2023-05-12 | 2023-05-12 | 一种基于启停工况分析评判的水轮机寿命评估方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116542177A true CN116542177A (zh) | 2023-08-04 |
Family
ID=87444981
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310542959.7A Pending CN116542177A (zh) | 2023-05-12 | 2023-05-12 | 一种基于启停工况分析评判的水轮机寿命评估方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116542177A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117634266A (zh) * | 2024-01-25 | 2024-03-01 | 成都市特种设备检验检测研究院(成都市特种设备应急处置中心) | 用于电梯制动器的紧急制动性能分析方法 |
-
2023
- 2023-05-12 CN CN202310542959.7A patent/CN116542177A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117634266A (zh) * | 2024-01-25 | 2024-03-01 | 成都市特种设备检验检测研究院(成都市特种设备应急处置中心) | 用于电梯制动器的紧急制动性能分析方法 |
CN117634266B (zh) * | 2024-01-25 | 2024-04-16 | 成都市特种设备检验检测研究院(成都市特种设备应急处置中心) | 用于电梯制动器的紧急制动性能分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113221280B (zh) | 一种基于数字孪生的滚动轴承建模与模型更新方法及系统 | |
CN110889238B (zh) | 一种数字孪生驱动的热压罐管控系统 | |
Chengwei et al. | Dynamic parametric modeling-based model updating strategy of aeroengine casings | |
US10436669B2 (en) | Interpolation engine for analysis of time-varying load data signals | |
CN103559418B (zh) | 一种基于故障物理的电子产品任务可靠度计算方法 | |
CN111737811A (zh) | 基于数字孪生的直升机动部件寿命管理方法、设备、介质 | |
CN116380166A (zh) | 一种设备异常监测方法、电子设备及存储介质 | |
CN116542177A (zh) | 一种基于启停工况分析评判的水轮机寿命评估方法及系统 | |
CN113569350B (zh) | 基于数字孪生的离心压缩机叶轮疲劳寿命预测方法 | |
Song et al. | A novel metamodeling approach for probabilistic LCF estimation of turbine disk | |
CN106339541A (zh) | 大型风力发电机组中塔筒门框的焊缝疲劳强度分析方法 | |
CN113420465B (zh) | 一种基于数字孪生模型的液压支架全寿命周期管理方法 | |
CN109492897B (zh) | 基于能量的暂态稳定和振荡动态的统一分析方法、装置及系统 | |
CN109029884A (zh) | 一种汽车悬臂结构件振动疲劳分析的方法 | |
CN115292849A (zh) | 基于相场法和bp神经网络的机械结构剩余寿命预测方法 | |
CN117034800A (zh) | 一种基于Fluent的往复式注水泵流场分析方法 | |
CN105005210B (zh) | 机电一体化仿真系统及使用其的方法 | |
CN110955933B (zh) | 一种基于响应面法的机械结构模糊疲劳可靠度计算方法 | |
KR102245922B1 (ko) | CFD(Computational Fluid Dynamics) 해석 이상 징후 예측 시스템 및 방법 | |
Zhang et al. | An implicit algorithm for high-order DG/FV schemes for compressible flows on 2D arbitrary grids | |
CN106529086A (zh) | 用于机械结构件的静载荷强度有限元仿真精度的确定方法 | |
CN113515848A (zh) | 矿山设备的结构设计优化方法、系统和可读存储介质 | |
US20130179132A1 (en) | Analysis Method, Apparatus and Software for a System With Frequency Dependent Materials | |
CN107180141A (zh) | 基于径向基代理模型的齿轮减速器箱体可靠性优化方法 | |
CN109657285A (zh) | 汽轮机转子瞬态应力的检测方法 |
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 |