一种使用非接触测量的涡轮叶片动应力计算方法技术

本发明专利技术提供了一种使用非接触测量的涡轮叶片动应力计算方法。包括：1、通过叶片非接触测量系统获取旋转态叶片的叶尖振动信息；2、通过有限元计算方法(FEM)获取叶片谐调叶盘模态信息；3、结合基本失谐模型(FMM)方法获得叶盘失谐模态加权因子和失谐模态振型；4、结合叶尖振动位移和失谐模态振型获得叶片响应计算模态叠加的主坐标向量；5、结合FEM计算的谐调模态应力和FMM计算的失谐模态振型获得叶盘失谐模态应力；6、根据失谐模态应力和步骤4得到的主坐标向量计算获得叶片动应力。本方法可通过叶尖计时实验直接在线计算获得叶片动应力。本方法减少了有限元计算时间；不需要获取实际叶盘精确的几何外形数据，节约了经济成本和时间。同时本文方法使用失谐模态叠加得到计算动应力，在原理上更加合理，不单纯依赖于数据库或者经验公式。

【技术实现步骤摘要】
一种使用非接触测量的涡轮叶片动应力计算方法
本申请涉及旋转叶片振动监测
，尤其是涉及一种涡轮机叶片动应力非接触测量方法。
技术介绍
涡轮机是国防和工业生产的关键设备，而叶片是涡轮机的核心部件，其健康状况对于涡轮机的安全、平稳、高效运行具有重大意义。涡轮机运行过程中叶片承受复杂的冲击载荷，容易产生振动问题，而叶片振动导致的疲劳断裂故障和事故尤为突出。因此，对叶片动应力(Dynamicstress)进行振动监测是很有必要的。目前对叶片进行动应力监测主要有两种方式：接触式测量方式和非接触式测量方式。接触式测量方式主要使用应变片直接粘贴于叶片表面，对旋转态叶片进行振动监测，这种方式会对叶盘叶片运行造成影响，同时也不能适应高温高速的叶盘运行环境。而非接触式测量方式成功避免了这一问题，叶尖计时(BladeTipTiming)方法就是一种典型的叶片非接触测量方式。如何通过非接触测量系统监测到的叶尖振动数据计算出叶片的动应力一直是一个难题。目前，基于叶尖计时方法的涡轮机叶片动应力的非接触测量方法，大多是先通过传感器监测叶片的叶尖振动信息，从而根据数据库信息或者经验公式推演计算出叶片的动应力数值。这种方式缺少必要的理论支撑，受叶片运行条件影响而具有较大的误差波动。对比文件：[1]王琰，郭定文，王秋蓉，黄文超.一种用于获得发动机转子叶片振动应力的方法[P].中国专利：CN104483119A,2015.04.01[2]X.关，J.张，S.K.周，N.H.乌莱里希，N.E.金，N.R.特夫斯.使用非接触测量和动态响应重构技术的涡轮机叶片疲劳寿命分析[P].中国专利：CN105008887A,2015.10.28[3]D.M.Feiner；J.H.Griffin.AFundamentalModelofMistuningforaSingleFamilyofModes[J].JournalofTurbomachinery,2002,Vol.124:597-605
技术实现思路
有鉴于此，本专利技术提供一种使用非接触测量的涡轮叶片动应力计算方法，以对涡轮机叶片的振动应力进行实时监测，提高动应力的计算精度，降低成本。本专利技术提出一种基于模态叠加计算理论、基本失谐模型方法(FMM)、有限元计算方法和涡轮机叶片非接触测量方法相结合的动应力在线计算方法，本方法主要解决如何由叶片振动的位移信息推演出叶片振动动应力的问题。该方法主要应用于透平机械中叶片动应力的在线监测中。为实现上述目的，本方案采用如下技术方案：1)通过叶尖计时系统获取旋转态叶片的叶尖振动信息；2)通过有限元计算方法(FEM)获取谐调叶盘叶片的模态振型、模态应力和模态频率，获取单个叶片悬臂梁状态下的模态频率；通过叶盘叶片的模态振型、模态应力，获得叶片谐调模态应力对叶片谐调模态位移的比值；3)结合步骤2)计算得到结果和基本失谐模型(FMM)方法，获得叶盘失谐模态加权因子和失谐模态振型；4)结合步骤1)获得的叶尖振动位移和步骤3)获得的失谐模态振型，获得叶片振动响应的模态叠加形式的主坐标向量；5)由步骤2)有限元计算的谐调模态应力，谐调模态振型和步骤3)获得的FMM计算的失谐模态加权因子获得叶盘失谐模态应力；6)根据步骤5)得到的失谐模态应力和步骤4)得到的主坐标向量计算得到叶片动应力值。进一步的，第1)步中通过叶尖计时系统获取的旋转态叶片的叶尖振动信息包括：叶片叶尖位移响应向量x(t)和每个叶片的振动固有频率其中x(t)＝[x1(t)x2(t)...xm(t)]T，向量中每个元素为同时刻每个叶片的叶尖振动位移值。x(t)中每个元素xi(t)可以是向量也可以是单个数值，当xi(t)为向量时，xi(t)为包含一个叶片多个点的振动信息，x(t)为矩阵形式；当xi(t)为单个数值时，xi(t)为包含一个叶片1个点的振动信息，x(t)为向量形式。进一步的，第2)步中通过有限元计算方法如下：叶盘模型为基于步骤(1)中叶尖计时监测的叶盘实物获得的，为叶盘实物的理想的谐调的叶盘模型，即计算中每个扇区是相同的。可以是每个叶片扇区均具有相同的结构参数和网格参数；也可以是建立单个扇区有限元模型，再进行基于循环对称方式的有限元计算。进而获取谐调叶盘叶片的模态振型、模态应力和模态频率获取单个叶片悬臂梁状态下的模态频率通过谐调叶盘叶片的模态振型、模态应力，获得叶片模态应力对叶片模态位移的比值进一步的，第3)步中获得叶盘失谐模态加权因子和失谐模态振型的计算方式如下：本步骤提到的基本失谐模型方法为D.M.Feiner等发表的论文《AFundamentalModelofMistuningforaSingleFamilyofModes》中介绍的方法。叶片悬臂梁状态下的模态频率谐调叶盘叶片的模态频率每个叶片的振动固有频率三者作为基本失谐模型方法的输入值，根据《AFundamentalModelofMistuningforaSingleFamilyofModes》文中方法构建特征矩阵求得特征向量即为叶盘失谐模态加权因子βj，其中第j阶振型的第s个叶片的模态位移为进一步的，第4)步中，由βj计算得到的组建失谐模态振型B，B为步骤3)中计算得到的失谐模态振型的全部或者一部分。有主坐标向量η(t)计算方式如下：η(t)＝B-1x(t)当x(t)包含全部叶片的位移响应信息时，B为叶片某一族模态的所有模态，B-1为B的逆矩阵。当x(t)只包含部分叶片的位移响应信息时，B取为叶片某一族模态的部分模态，B-1为B的广义逆矩阵。以上可求得模态叠加形式的主坐标向量η(t)。进一步的，第5)步中计算叶盘失谐模态应力方式如下：定义第j阶第s个叶片的失谐模态应力为σjs：其中为由步骤2)计算得到叶片谐调模态应力对叶片模态位移的比值，βjm为步骤3)获得的失谐模态加权因子βj中元素。由所有σjs元素即可组成所有叶片的第j阶失谐模态应力σj，由所有阶失谐模态应力组成失谐模态应力矩阵。进一步的，第6)步中计算叶片动应力值σD计算方式如下：其中σj为步骤5)中所得的第j阶失谐模态应力，ηj为步骤4)所得的主坐标向量η(t)中的第j个元素。与传统有限元计算动应力方法相比，本计算方法考虑了实际叶盘的必然存在的失谐情况，同时不需要建立精确的失谐有限元模型，大大减小了有限元计算的工作量。又可以同叶尖计时方法结合，实现实时的在线叶片动应力计算与监测。与传统基于叶尖计时动应力计算方式相比，本计算考虑了叶片振动响应时的多模态振动情况；而不像传统计算方式一样仅仅将叶片振动当做某一阶模态的振动，将叶尖振动位移与叶片应力作为简单的线性处理；综上本方法将获得更高的动应力计算精度。综上所述，本专利技术旨在提供一种基于叶尖计时方法的叶片动应力计算方法。一方面可以结合叶尖计时方法等叶片振动非接触测量方法实现叶片动应力实时在线监测，另一方面可以减小叶片动应力计算误差，提高叶片剩余寿命预测精度。本方法对涡轮叶片的健康监测具有重要意义。附图说明为了更清楚地说明本申请实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本申请的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获本文档来自技高网...

【技术保护点】
1.一种使用非接触测量的涡轮叶片动应力计算方法，其特征在于，所述监测计算动应力方法包括以下步骤：(1)通过叶尖计时系统获取旋转态叶片的叶尖振动信息；(2)通过有限元计算方法(FEM)获取谐调叶盘叶片的模态振型、模态应力和模态频率，获取单个叶片悬臂梁状态下的模态频率；通过谐调叶盘叶片的模态振型、模态应力，获得叶片模态应力对叶片模态位移的比值。(3)结合步骤(2)计算得到结果和基本失谐模型(FMM)方法，获得叶盘失谐模态加权因子和失谐模态振型；(4)结合步骤(1)获得的叶尖振动位移和步骤(3)获得的失谐模态振型，获得叶片振动响应计算的模态叠加形式的主坐标向量；(5)由FEM计算的谐调模态应力，谐调模态振型和步骤(3)获得的FMM计算的失谐模态加权因子获得叶盘失谐模态应力；(6)根据失谐模态应力和步骤4得到的主坐标向量计算得到叶片动应力值。

【技术特征摘要】
1.一种使用非接触测量的涡轮叶片动应力计算方法，其特征在于，所述监测计算动应力方法包括以下步骤：(1)通过叶尖计时系统获取旋转态叶片的叶尖振动信息；(2)通过有限元计算方法(FEM)获取谐调叶盘叶片的模态振型、模态应力和模态频率，获取单个叶片悬臂梁状态下的模态频率；通过谐调叶盘叶片的模态振型、模态应力，获得叶片模态应力对叶片模态位移的比值。(3)结合步骤(2)计算得到结果和基本失谐模型(FMM)方法，获得叶盘失谐模态加权因子和失谐模态振型；(4)结合步骤(1)获得的叶尖振动位移和步骤(3)获得的失谐模态振型，获得叶片振动响应计算的模态叠加形式的主坐标向量；(5)由FEM计算的谐调模态应力，谐调模态振型和步骤(3)获得的FMM计算的失谐模态加权因子获得叶盘失谐模态应力；(6)根据失谐模态应力和步骤4得到的主坐标向量计算得到叶片动应力值。2.如权利要求1所述的方法，其特征在于步骤(1)中的叶尖计时系统对旋转状态的叶盘...

【专利技术属性】
技术研发人员：王维民，户东方，张旭龙，李维博，
申请(专利权)人：北京化工大学，
类型：发明
国别省市：北京,11