申请/专利权人：西安交通大学

申请日：2020-07-24

公开（公告）日：2024-04-05

公开（公告）号：CN112084583B

主分类号：G06F30/17

分类号：G06F30/17;G06F30/20;G06F119/04

优先权：

专利状态码：有效-授权

法律状态：2024.04.05#授权;2021.01.01#实质审查的生效;2020.12.15#公开

摘要：本发明公开了基于数字孪生的转子叶片寿命预测方法和预测系统，所述方法包括：基于转子叶片尺寸和材料参数建立转子叶片的初始三维有限元模型；叶端定时测量获取转子叶片振动参数和叶端振动位移数据，实时更新初始三维有限元模型获得数字孪生模型及转子叶片坎贝尔图；基于应力重构和叶端振动位移数据计算所述数字孪生模型的转子叶片应力；基于所述转子叶片应力，获取应力时间历程以计算转子叶片寿命。

主权项：1.一种基于数字孪生的转子叶片寿命预测方法，所述方法包括以下步骤：第一步骤中，基于转子叶片尺寸和材料参数建立转子叶片的初始三维有限元模型；第二步骤中，叶端定时测量获取转子叶片振动参数和叶端振动位移数据，实时更新初始三维有限元模型获得数字孪生模型及转子叶片坎贝尔图，叶端定时测量包括，在转子叶片机匣周向布置多个叶端定时传感器测量转子叶片旋转状态下到达传感器的时间，转子叶片在无振动下的到达时间为ts：式中，R为叶端点到回转中心的旋转半径，ωR为转子转频，NB为叶片数目，转子叶片实际到达某叶端定时传感器为te，则所述转子叶片在该转测得的叶尖位移为d，d＝ωRRte-ts＝ωRRΔt，式中，Δt为转子叶片旋转状态下无振动和产生振动而到达叶端定时传感器的时间差，获得各转子叶片到达各叶端定时传感器的时间差以计算叶端振动位移，基于各传感器测得的转子叶片每圈振动位移，采用周向傅里叶拟合算法计算转子叶片振动参数，令某转子叶片的振动共有m阶振动模态振型，其振动信号Xt为： 式中，Xit表示解耦后的第i阶模态振动信号，c表示叶片静变形，ωi表示叶片多模态振动圆周频率，Ai表示第i阶振动最大振幅，表示第i阶振动初始相位，t表示叶片振动时刻，由此得到转子叶片振动信号以拟合叶片振动曲线，基于叶片坎贝尔图得到转子叶片多振动模态的振动参数，基于求转子叶片多模态振动参数，基于有限元修正和更新初始三维有限元模型获得转子叶片的数字孪生模型；第三步骤中，基于应力重构和叶端振动位移数据计算所述数字孪生模型的转子叶片应力，转子叶片叶端定时测点对应数字孪生模型中自由度a，共有m阶模态参与转子叶片振动响应，数字孪生模型中测点对应的位移模态振型Ua为：Ua＝[Ua1,Ua2,…,Uam]，其中，Ua大小为1×m，当有m阶模态参与转子叶片振动响应时，n自由度转子叶片全场等效应力模态振型矩阵Vn为：Vn＝[Vn1,Vn2,…,Vnm]，其中，Vn大小为n×m，构造转子叶片叶端测点位移响应与叶片全场应力的转换矩阵： 其中，转换矩阵T的大小为n×m，Uai表示测点第i阶位移模态振型，由换算矩阵T和所述解耦的多模态振动信号Xt，计算转子叶片任意时刻t的等效应力：σt＝TXt＝T[X1t,X2t,…,Xmt]T；第四步骤中，基于所述转子叶片应力，获取应力时间历程以计算转子叶片寿命，基于所述转子叶片应力获取应力时间历程包括，采用雨流循环计数将不规则应力时间历程转化成循环预定载荷历程的应力加载历程，其中，预定循环载荷历程中包含k个不同应力水平的应力幅，按照线性疲劳累积损伤Palmgren-Miner模型假设，在k种应力水平作用下，转子叶片到出现裂纹时的疲劳重复次数Bf为： 其中，Ni表示第i个应力水平下的循环次数，Nfi表示第i个应力水平下叶片达到破坏时的循环次数，由此计算转子叶片运行到出现裂纹时的疲劳重复次数Bf为疲劳寿命，第i个应力水平下叶片达到破坏时的循环次数Nfi由平均应力的S-N曲线获得： 其中，σai表示第i个应力水平下的应力幅；σm来自数字孪生模型计算的平均应力；σ'f和b分别表示转子叶片材料疲劳强度系数参数和疲劳强度指数。

全文数据：

权利要求：

百度查询： 西安交通大学 基于数字孪生的转子叶片寿命预测方法和预测系统

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD