一种基于有限元模型的汽轮机转子蠕变损伤实时计算方法与流程

本发明涉及汽轮机转子蠕变损伤计算与评估领域，特别是涉及有限元数值计算方法和多元回归实时计算方法。

背景技术：

转子是汽轮机机组的核心部件，通常运行在高温高压的恶劣工作条件，容易产生裂纹，金属材料会出现损伤。转子的损伤主要包括低周疲劳损伤和高温蠕变损伤，低周疲劳损伤和高温蠕变损伤对转子产生的作用是耦合的，在机组实际运行中有必要考虑到两类损伤的耦合作用。金属材料的蠕变损伤主要是指长期高温环境下，塑性应变逐渐变大的现象。转子金属材料一般是低碳合金钢，当汽轮机机组进入平稳状态后，各类热工参数逐渐稳定，热应力变化较小，其损伤形式由低周疲劳损伤过渡到高温蠕变损伤。危害主要体现在，由于金属蠕变损伤的累积，转子部件发生过量的塑性变形而无法使用；或者蠕变进入到加速发展阶段，材料会发生蠕变破裂，导致严重的安全事故。因此，对于汽轮机机组，尤其是最昂贵的转子部件，需要进行连续的、实时的、高精度的蠕变损伤计算。

现有技术中，既有单独计算高温蠕变损伤的方法，也有考虑疲劳-蠕变耦合情况的方法，例如非线性力学方法、数值计算法等，这类方法理论基础好，计算精度高，但是算法复杂，计算速度慢。另外一类是解析法，这类算法强调计算速度，对损伤计算过程和模型做了一定的简化，而牺牲了计算精度。这两类方法对热能装备的损伤计算都具有积极的指导作用。但是，如何将二者相结合，既能满足热能装备的蠕变损伤实时计算，又可以保证计算精度，是现在急需解决的一个问题。

技术实现要素：

鉴于以上所述现有技术的缺点，本发明的目的在于提供一种基于有限元模型的汽轮机转子蠕变损伤实时计算方法，用于解决现有技术对汽轮机转子蠕变损伤计算中实时性与高精度难以共融的问题。

为实现上述目的及其他相关目的，提供一种基于有限元模型的汽轮机转子蠕变损伤实时计算方法，其特征在于，所述方法包括：通过有限元分析软件计算转子的温度场和应力场；提取关键部位的温度和应力数据；将蠕变损伤计算相关的温度和应力数据进行归一化处理；采用多元回归方法建立蠕变应力实时计算模型；解析并拟合应力-损伤函数关系，实时计算蠕变损伤；建立蠕变损伤实时计算系统结构。

优选地，通过有限元分析软件计算汽轮机转子的温度场和应力场。

(1)温度场计算可采用二维轴对称模型，考虑温度函数为t(z，r，t)，则有：

转子外叶轮外表面与蒸汽对流换热，为第三类边界条件：

其中，r为转子径向参量；r0为表面半径；h为蒸汽与转子换热系数；tf为蒸汽温度。转子与轴承润滑油接触部位可认为是第一类边界条件，即温度与润滑油温相等:

转子两端忽略转子热流密度。边界条件表示如下；

其中，γ表示转子中心及两端边界，q为热流密度。

(2)应力场计算时对转子固定点部位进行零位移约束具体为：

u|z＝0＝0

优选地，提取关键部位的温度和应力数据。一般情况下转子调节级根部是应力较大的部位，可提取调节级根部应力最大处的温度和应力曲线，再通过数据导出模块将所需的温度数据和应力数据分别导出保存。

优选地，将蠕变损伤计算相关的温度、压力、转速以及应力数据进行归一化处理。由于不同评价指标往往具有不同的量纲和量纲单位，这样的情况会影响到数据分析的结果。为了消除指标之间的量纲影响，数据集需要进行数据标准化处理，具体为：

其中max为原始样本数据的最大值，min为原始样本数据的最小值，x^*为归一化后的数值，x为归一化前的数值。

优选地，采用多元回归方法建立蠕变应力实时计算模型。多元线性回归模型：

y＝b0+b1x1+l+bmxm+εε～n(0,σ²)

m元线性回归方程估计值为：

目标函数为最小化残差平方和：

其中x为自变量；y为因变量；b0～bm为回归系数；ε为随机变量；n为满足期望为0，方差为σ²的正态分布；x为自变量x组成的矩阵；y为n个因变量y组成的矩阵；x^t为x的转置矩阵；为b0～bm的估计值；为y的估计值；为yi的估计值；σ²为方差；为σ²的无偏估计值；qe为残差平方和。

优选地，解析并拟合应力-损伤函数关系，实时计算蠕变损伤。已知最大应力，根据以下三式计算高温蠕变损伤。

(1)循环应力-应变关系式为：

其中e是弹性模量，k'是循环强度系数，n'是循环应变硬化指数。

(2)蠕变损伤发生在平稳状态，此时主蒸汽温度为540℃的时候，应变与循环致裂周次nf关系为：

ε＝0.00332(2nf)^-0.0697+0.6264(2nf)^-0.7553，

(3)损伤d和应变ε的关系为：

其中的多项式拟合系数p1＝9.7904×10¹²，p2＝-3.8943×10¹¹，p3＝6.011×10⁹，p4＝-4.6037×10⁷，p5＝1.8292×10⁵，p6＝-3.3086×10²，p7＝0.27029，p8＝0。

优选地，建立蠕变损伤实时计算系统架构。将训练温度输入到adina有限元软件中来求得热应力数据，再将压力和分析得到的离心力和求得的热应力一起输入到多元回归模型中，同时将训练的温度、压力、转速以及实测的温度、压力、转速也输入到多元回归模型中进行拟合对比分析，将拟合后得到的应力进行应力-应变的关系式得到应变，再通过应变-损伤的关系模型进行多项式拟合，计算得到蠕变损伤。

附图说明

图1显示为一种基于有限元模型的汽轮机转子蠕变损伤实时计算方法流程示意图。

图2显示为某国产300mw汽轮机转子冷态启动温度曲线。

图3显示为冷态启动转子蠕变应力曲线。

图4显示为冷态启动蠕变损伤应力-应变曲线。

图5显示为蠕变应力计算训练数据曲线。

图6显示为蠕变应力计算测试数据曲线。

图7显示为高温蠕变损伤实时计算曲线图。

图8显示为一种基于有限元模型的汽轮机转子蠕变损伤实时计算系统架构。

具体实施方式

以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易的了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。

请参阅图1至图8。需要说明的是，本实施例中所提供的图示仅以示意方式说明本发明的基本构想，遂图示中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制，其实际实施时各组件的形态、数量及比例可为一种随意的改变，且其组件布局形态也可能更为复杂。

在汽轮机转子高温蠕变损伤计算领域，计算方法主要可以分为理论分析法和解析法。前者理论基础好，计算精度高，但是算法复杂，计算速度慢。后者强调计算速度，对损伤计算过程和模型做了简化，牺牲了计算精度。目前需要寻求一种可以具备一定精度的转子高温蠕变损伤实时计算方法，本发明就是基于这些理念而形成的。

本发明的目的在于提供一种转子蠕变损伤实时计算方法，用以解决现有方法中转子蠕变损伤实时计算和高精度计算难以共融的问题。以下将详细描述本发明的一种基于有限元模型的汽轮机转子蠕变损伤实时计算方法的原理和实施方式，使本领域技术人员不需要创造性劳动即可理解本发明的一种基于有限元模型的汽轮机转子蠕变损伤实时计算方法。

如图1所示，本发明一种基于有限元模型的汽轮机转子蠕变损伤实时计算方法，所述方法步骤包括：

s1，通过有限元分析软件计算转子的温度场和应力场；

s2，提取关键部位的温度和应力数据；

s3，将蠕变损伤计算相关的温度、压力、转速以及应力数据进行归一化处理；

s4，采用多元回归方法建立蠕变应力实时计算模型；

s5，解析并拟合应力-损伤函数关系，实时计算蠕变损伤；

s6，建立蠕变损伤实时计算系统架构。

下面结合具体实施方式对本发明进行详细说明。本发明实施在adina8.5软件和matlabr2016a环境下完成。具体方法如下：使用的汽轮机转子是以某电厂的1号机组为原型，该汽轮机是单轴、凝汽式汽轮机。转子的材料为30cr1mo1v。高压转子的长度约为4800mm，具有单列调节级和11级的压力级。通过adina建立该转子二维模型，根据《国产300mw汽轮机运行导则》，在调节级和高压级施加工况模拟。

首先执行步骤s1，使用有限元分析软件计算汽轮机转子温度场和应力场。转子呈对称形状，建立适当的几何模型，同时施加荷载、确定边界条件及材料参数，通过有限元分析软件进行分析求解。

在步骤s2中，提取关键部位的温度和应力数据。根据汽轮机转子的温度场和应力场云图可以发现，转子前几级温度较大，且温度梯度较为明显；从应力场云图也可看到，转子调节级根部是应力较大的部位。提取调节级根部应力最大处的温度曲线和应力曲线以及温度和应力的等值线，再通过数据导出模块将所需的温度数据、疲劳损伤实测应力数据、蠕变损伤实测应力数据和蠕变损伤应力-应变数据导出为excel文档并保存。某国产300mw汽轮机转子冷态启动温度曲线如图2所示，冷态启动转子蠕变应力曲线如图3所示，冷态启动蠕变损伤应力-应变曲线如图4所示。

在步骤s3中，将蠕变损伤计算相关的温度、压力、转速以及应力数据进行归一化处理。原始数据经过数据归一化处理后，各指标处于同一数量级，本发明采用min-max标准化对原始数据的线性变换，使结果值映射到[min,max]之间。转换函数为：在matlab中调用mapminmax函数，对温度、压力、转速以及应力数据进行归一化处理。

在步骤s4中，采用多元回归方法建立蠕变应力实时计算模型，根据多元回归模型y＝b0+b1x1+l+bmxm+εε～n(0,σ²)，采用最小二乘法和求偏导数，整理后得到方程组，引入矩阵：

则方程组可写成：x^txb＝x^ty，于是可解得：

的估计值：

方程为m元线性回归方程。在matlab软件中调用regress多元回归函数、mse均方误差函数对预测数据和实测数据进行拟合分析。蠕变应力计算训练数据图像如图5所示，蠕变应力计算测试数据图像如图6所示。

在步骤s5中，解析并拟合应力-损伤函数关系，实时计算蠕变损伤。如果已知最大应力值，那么可以通过过程量的数值计算得到高温蠕变损伤。循环应力-应变关系式为：其中e是弹性模量，k'是循环强度系数，n'是循环应变硬化指数。最大应力发生在主蒸汽温度达到540℃的时候，此时汽轮机转子高温蠕变损伤可以通过下式计算得到：ε＝0.00332(2nf)^-0.0697+0.6264(2nf)^-0.7553，循环致裂次数nf与应变ε的关系可以通过多项式拟合得到，并可以进一步得到以金属材料应变作为输入，以循环致裂次数nf作为输出的关系式。因此，通过多项式拟合重新建立起应力-损伤关系模型，解决了损伤的实时计算问题。其中损伤d和应变ε的关系可通过多项式拟合得到，高温蠕变损伤实时计算曲线如图7所示。

在步骤s6中，建立蠕变损伤实时计算系统架构。将训练温度输入到adina有限元软件中来求得热应力数据，再将压力和分析得到的离心力和求得的热应力一起输入到多元回归模型中，同时将训练的温度、压力、转速以及实测的温度、压力、转速也输入到多元回归模型中进行拟合对比分析，将拟合后得到的应力进行应力-应变的关系式得到应变，再通过应变-损伤的关系模型进行多项式拟合，从中即可得到蠕变损伤，如图8所示。