CN109883389B - 一种旋转叶片动应变场测量方法及其系统 - Google Patents
一种旋转叶片动应变场测量方法及其系统 Download PDFInfo
- Publication number
- CN109883389B CN109883389B CN201910225544.0A CN201910225544A CN109883389B CN 109883389 B CN109883389 B CN 109883389B CN 201910225544 A CN201910225544 A CN 201910225544A CN 109883389 B CN109883389 B CN 109883389B
- Authority
- CN
- China
- Prior art keywords
- strain
- blade
- rotating blade
- matrix
- finite element
- 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.)
- Active
Links
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种旋转叶片动应变场测量方法及其系统,所述方法包括以下步骤:建立待测量旋转叶片的三维有限元模型,提取所述三维有限元模型的模态参数;在所述旋转叶片上布置多个应变片,确定所述应变片测点数目、位置和方向;构造应变片测点动应变与全场动应变的转换矩阵;基于所述应变片测量所述旋转叶片相应位置动应变,所述动应变基于所述转换矩阵经由模态处理得到所述旋转叶片任意时刻、任意位置及任意方向的动应变。本发明提供的方法仅利用极少测点便可实现旋转叶片整体动应变场的测量,不仅可实现叶片表面动应变的测量,还可以实现叶片内部节点正应变、剪应变的测量,计算过程简单,测量精度高,易于在线测量,节约大量应变片。
Description
技术领域
本发明属于航空发动机、燃气轮机旋转叶片振动测试技术领域,特别是一种旋转叶片动应变场测量方法及其系统。
背景技术
高速旋转叶片的完整性直接影响航空发动机整体结构的安全运行,受工作环境苛刻、载荷强交变等因素的影响,其在服役过程中极易产生振动疲劳裂纹而导致严重事故。叶片振动过大导致的高周疲劳是航空发动机叶片主要失效模式。叶片高周疲劳主要由各种气动载荷、机械载荷导致的动应力引起,在短时间内便可累计大量循环产生疲劳裂纹,特别是当叶片发生共振时动应力极易导致叶片疲劳失效。在航空发动机研制、生产过程中,为了掌握叶片振动特性,需要对叶片动应变/应力进行测量;在发动机使用过程中,需要对发动机叶片健康状况进行实时监测、预警和诊断。通过对叶片动应变进行测量,可确定叶片在发动机工作转速范围内是否存在共振的危险;测量叶片动应力与发动机工作状态之间的特性关系,可确定叶片在不同状态下能否可靠工作。通过在旋转叶片表面粘贴应变片,仅能测量有限叶片有限位置动应变,其可靠性和持续工作时间较低,特别是高温环境下在涡轮叶片布置大量应变片常常只有很少的应变片可以获取有效信息,存活率极低。
在背景技术部分中公开的上述信息仅仅用于增强对本发明背景的理解,因此可能包含不构成在本国中本领域普通技术人员公知的现有技术的信息。
发明内容
针对现有技术中存在的问题,本发明提出一种旋转叶片动应变场测量方法及其系统,解决了应变片仅能测量叶片有限位置动应变、不能测量叶片整体应变场的难题,并具备同时测量旋转叶片表面与内部正应变、剪应变的优势。
不同的振动模态对叶片响应的贡献量有较大差别,一般前若干阶模态对系统振动贡献较大,其他模态(主要是高阶模态)可以忽略不计。旋转叶片动应变场测量技术是基于有限测点的应变信息反演重构整个应变场,模态降阶与扩张理论基于各阶模态振型关于质量矩阵的正交性,利用振动系统有限的模态来表示多自由度系统振动。模态降阶将物理空间中关于有限测点动应变的动力学方程转化到模态空间关于模态响应的多自由度解耦系统;模态扩展基于应变模态振型将模态空间的模态响应转换为物理空间的动应变。为此,本发明基于模态降阶与扩展理论通过有限测点的振动反演重构旋转叶片整体动应变场,核心是建立叶片有限测点动应变与全场动应变的转换矩阵。
本发明的目的是通过以下技术方案予以实现,一种旋转叶片动应变场测量方法包括以下步骤:
第一步骤中,建立待测量旋转叶片的三维有限元模型,提取所述三维有限元模型的模态参数;
第二步骤中,在所述旋转叶片上布置多个应变片,确定所述应变片测点数目、位置和方向;
第三步骤中,构造应变片测点动应变与全场动应变的转换矩阵;
第四步骤中,基于所述应变片测量所述旋转叶片相应位置动应变;
第五步骤中,所述动应变基于所述转换矩阵经由模态处理得到所述旋转叶片任意时刻、任意位置及任意方向的动应变。
所述的方法中,第一步骤中,通过模态分析提取所述三维有限元模型前nm阶模态参数、模态频率fi和大小为ndof×1的应变模态振型ψi,构造旋转叶片全场应变模态振型矩阵大小为2ndof×nm,其中,nm表示模态数目,i表示模态阶次,ndof表示旋转叶片有限元模型的自由度数目,ndof=3nn,nn表示旋转叶片有限元模型节点。
所述的方法中,第一步骤中,每个旋转叶片有限元模型节点的应变包含3个正应变εx、εy、εz与3个剪应变γxy、γyz、γxz共6个应变分量。
所述的方法中,第二步骤中,旋转叶片应变片测点数目nd大于等于模态数目nm。
所述的方法中,第二步骤中,基于旋转叶片应变模态振型构造应变模态振型的测点选择矩阵大小为2nc×nm,其中nc表示旋转叶片表面有限元网格节点数目;2表示每个节点有2个方向可安装应变片;从测点选择矩阵ψp中随机选择nd个测点构造大小为nd×nm的测点应变模态振型矩阵ψd,并计算其矩阵条件数κ;随机过程重复R次并从中选择矩阵条件数κ最小时的测点布置。
所述的方法中,第五步骤中,基于转换矩阵T和应变实测的位移信号s(t),旋转叶片t时刻叶片表面和内部所有节点应变经由公式S(t)=Ts(t)计算得出,
应变
,所述应变s(t)包括正应变和剪应变;其中,εi,x表示叶片有限元模型第i个节点x方向的正应变,εi,y表示叶片有限元模型第i个节点y方向的正应变,εi,z表示叶片有限元模型第i个节点z方向的正应变,γi,xy表示叶片有限元模型第i个节点x-y方向的剪应变,γi,yz表示叶片有限元模型第i个节点y-z方向的剪应变,γi,xz表示叶片有限元模型第i个节点x-z方向的剪应变。
所述的方法中,第五步骤中,所述模态处理基于所述应变模态振型模态降阶与扩展。
根据本发明的另一方面,一种实施所述方法的测量系统包括,
多个应变片,其布置在旋转叶片上;
计算单元,其连接所述动应变测量模块,所述计算单元包括,
测点优选模块,其配置成优化布置在所述旋转叶片上的应变片的测点数目、位置和方向,其中,基于旋转叶片可安装应变片测点的应变模态振型构造旋转叶片全场应变模态振型矩阵从测点选择矩阵ψp中随机选择nd个测点构造大小为nd×nm的测点应变模态振型矩阵ψd并计算其矩阵条件数κ,随机过程重复R次并从中选择矩阵条件数κ最小时的测点布置,
转换矩阵计算模块,其配置成构造测点动应变与全场动应变的转换矩阵,
动应变场重构模块,用于计算旋转叶片t时刻叶片表面和内部所有节点应变经由公式S(t)=Ts(t)计算得出,
应变
所述应变S(t)包括正应变和剪应变。
有益效果
本发明提供的旋转叶片动应变场重构方法仅利用极少测点动应变便可实现旋转叶片整体动应变场的测量。不仅可实现叶片表面动应变的测量,还可以实现叶片内部节点正应变、剪应变的测量,本发明提供的方法可实现利用有限测点正应变反演重构任意节点剪应变,本发明构造的应变-应变转换矩阵,与频率、时间、甚至边界条件无关。本发明提供的方法计算过程简单,易于在线测量,可节约大量应变片。本发明提供的方法测量精度高,可评估应变片实测数据是否有效。本发明提供的旋转叶片动应变场重构系统过程简单,易于实现。
附图说明
通过阅读下文优选的具体实施方式中的详细描述,本发明各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。说明书附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。显而易见地,下面描述的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。而且在整个附图中,用相同的附图标记表示相同的部件。
在附图中:
图1是本发明提供的一种旋转叶片动应变场重构方法优先实例的流程示意图;
图2(a)至图2(b)是本发明提供的一种旋转叶片动应变场重构系统的结构示意图,其中,图2(a)旋转叶片动应变场重构系统组成;图2(b)动应变测量模块;
图3是一个实施例中模拟转子叶片动载荷激励位置与应变片测点位置示意图;
图4(a)至图4(c)是一个实施例中旋转叶片的应变模态振型,其中,图4(a)一弯振型;图4(b)一扭振型;图4(c)二弯振型;
图5是一个实施例中旋转叶片三个应变片实测信号;
图6(a)至图6(b)是一个实施例中重构的旋转叶片应变场中叶身748号节点动应变与真实动应变比较结果,其中,图6(a)三个正应变分量与真实动应变比较结果;图6(b)三个剪应变分量与真实动应变比较结果;
图7(a)至图7(b)是一个实施例中重构的旋转叶片应变场中叶根1100号节点动应变与真实动应变比较结果,其中,图7(a)三个正应变分量与真实动应变比较结果;图7(b)三个剪应变分量与真实动应变比较结果。
以下结合附图和实施例对本发明作进一步的解释。
具体实施方式
下面将参照附图1至图7(b)更详细地描述本发明的具体实施例。虽然附图中显示了本发明的具体实施例,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。
需要说明的是,在说明书及权利要求当中使用了某些词汇来指称特定组件。本领域技术人员应可以理解,技术人员可能会用不同名词来称呼同一个组件。本说明书及权利要求并不以名词的差异来作为区分组件的方式,而是以组件在功能上的差异来作为区分的准则。如在通篇说明书及权利要求当中所提及的“包含”或“包括”为一开放式用语,故应解释成“包含但不限定于”。说明书后续描述为实施本发明的较佳实施方式,然所述描述乃以说明书的一般原则为目的,并非用以限定本发明的范围。本发明的保护范围当视所附权利要求所界定者为准。
为便于对本发明实施例的理解,下面将结合附图以具体实施例为例做进一步的解释说明,且各个附图并不构成对本发明实施例的限定。
为了更好地理解,图1为一个旋转叶片动应变场测量方法工作流程图,如图1所示,一种旋转叶片动应变场测量方法包括以下步骤:
第一步骤S1中,建立待测量旋转叶片的三维有限元模型,提取所述三维有限元模型的模态参数;
第二步骤S2中,在所述旋转叶片上布置多个应变片,确定所述应变片测点数目、位置和方向;
第三步骤S3中,构造应变片测点动应变与全场动应变的转换矩阵;
第四步骤S4中,基于所述应变片测量所述旋转叶片相应位置动应变;
第五步骤S5中,所述动应变基于所述转换矩阵经由模态处理得到所述旋转叶片任意时刻、任意位置及任意方向的动应变。
所述的方法的一个实施方式中,第一步骤S1中,通过模态分析提取所述三维有限元模型前nm阶模态参数、模态频率fi和大小为ndof×1的应变模态振型ψi,构造旋转叶片全场应变模态振型矩阵大小为2ndof×nm,其中,nm表示模态数目,i表示模态阶次,ndof表示旋转叶片有限元模型的自由度数目,ndof=3nn,nn表示旋转叶片有限元模型节点。
所述的方法的一个实施方式中,第一步骤S1中,每个旋转叶片有限元模型节点的应变包含3个正应变εx、εy、εz与3个剪应变γxy、γyz、γxz共6个应变分量。
所述的方法的一个实施方式中,第二步骤S2中,旋转叶片应变片测点数目nd大于等于模态数目nm。
所述的方法的一个实施方式中,第二步骤S2中,基于旋转叶片应变模态振型构造应变模态振型的测点选择矩阵大小为2nc×nm,其中nc表示旋转叶片表面有限元网格节点数目;2表示每个节点有2个方向可安装应变片;从测点选择矩阵ψp中随机选择nd个测点构造大小为nd×nm的测点应变模态振型矩阵ψd,并计算其矩阵条件数κ;随机过程重复R次并从中选择矩阵条件数κ最小时的测点布置。
所述的方法的一个实施方式中,第五步骤S5中,旋转叶片t时刻叶片表面和内部所有节点应变经由公式S(t)=Ts(t)计算得出,
应变
所述应变S(t)包括正应变和剪应变。
所述的方法的一个实施方式中,第五步骤S5中,所述模态处理基于所述应变模态振型模态降阶与扩展。
为了进一步理解本发明,下面结合附图1至图7(b)及具体实施例对本发明作进一步描述,应该强调的是,下述说明仅仅是示例性的,而本发明的应用对象不局限下述示例。
图1是本发明完成的一种旋转叶片动应变场重构方法的流程示意图,该方法基于模态降阶与扩展理论构造旋转叶片有限动应变测点与全场节点动应变的转换关系,实现叶片动应变场重构;图2(a)至图2(b)是本发明提供的一种旋转叶片动应变场重构系统的结构示意图,其中,1-转子机匣;2-应变片;3-转子叶片;4-轮盘;5-转子;6-转速传感器,方法具体步骤如下:
1)提取叶片三维有限元模型的模态参数:参见图3,利用ANSYS有限元分析软件建立模拟转子直板叶片的三维有限元模型,其中材料为铝,密度2700kg/m3,泊松比0.33,弹性模量72000MPa;叶片长48mm,厚度1mm,宽20mm;有限元单元类型为实体单元SOLID185,节点总数为3153;叶根两侧面固定约束,模拟转子叶片实际工作状态;
利用ANSYS模态分析模态提取前3阶模态参数,即nm=3:模态频率fi、大小为ndof×1的位移模态振型φi、大小为2ndof×1的应变模态振型ψi,其中,前三阶模态频率分别为f1=333.08Hz、f2=1806.03Hz、f3=2076.52Hz;构造旋转叶片全场应变模态振型矩阵大小为2ndof×nm,应变模态振型见图4(a)至图4(c);i表示模态阶次,ndof=9459表示叶片有限元模型的自由度数目,则ndof=3nn,nn=3153表示叶片有限元模型节点的数目;应变包含3个正应变εx、εy、εz与3个剪应变γxy、γyz、γxz共6个应变分量,即每个节点有6个应变模态振型。
2)确定叶片应变片测点数目、位置、方向:旋转叶片应变片测点数目nd不得小于关注的模态数目nm,即nd≥nm;本案例中,关注模拟转子叶片前三阶振动模态,取nm=3;应变片数目取最少,即nd=3;
提取叶片可供安装应变片的有限元模型节点[200,500]区间的应变模态振型,构造关于应变模态振型的测点选择矩阵大小为2nc×nm;其中nc=301表示叶片表面可供安装应变片的节点数目;2表示每个节点有Y和Z两个方向可供安装应变片;
从测点选择矩阵ψp中随机选择nd=3个测点;构造大小为nd×nm=3×3的测点应变模态振型矩阵ψd,并计算其矩阵条件数κ;这个随机过程重复R=1000次,并从中选择条件数κ最小时的测点布局方案;测点优选结果见图3:选中的三个测点分别为335号节点、345号节点、460号节点且均在Y向安装应变片,对应的测点应变模态振型矩阵ψd的条件数为10.50。
4)利用有限应变片获取旋转叶片有限位置动应变:在ANSYS有限元软件中对模拟叶片进行瞬态分析,质量阻尼系数设定为α=12.1380,刚度阻尼系数设定为β=8.1986×10-8,模拟气动载荷对旋转叶片的多模态振动,对叶片叶端1117号节点X方向施加多频简谐激励f(t)=cos(2πf1t)+10cos(2πf2t)+20cos(2πf3t),获得叶片真实的应变场,作为重构结果的参考;图3中三个应变片获取旋转叶片t个时刻的动应变时域信号s(t)=[s1(t),s2(t),s3(t)]T,实测结果参见图5;其中,采样频率fs=25000Hz,信号的数据长度为N=7500,采样时间为t=N/fs=0.3s。
基于模态降阶与扩展理论实现旋转叶片任意时刻、任意位置、任意方向的动应变测量:计算转子叶片t∈[0,0.3]s时刻叶片表面和内部所有节点应变经由公式S(t)=Ts(t)计算得出,
应变
所述应变S(t)包括正应变和剪应变。
取旋转叶片叶身748号节点和叶根1100号节点作为动应变场高精度重构的典型代表(见图3),结论同样适用于其他节点。其中,图6(a)至图6(b)是一个实施例中重构的叶片应变场中叶身748号节点动应变与真实动应变比较结果;图7(a)至图7(b)是一个实施例中重构的叶片应变场中叶根1100号节点动应变与真实动应变比较结果。从图6(a)至图6(b)和图7(a)至图7(b)中可知,重构的动应变信号与真实动应变高度吻合,几乎难以区分;为了定量评价本发明的旋转叶片动应变场测量方法的性能,在t∈[0,0.3]s区间计算重构信号与真实应变的相对误差,图6(a)中叶身748号节点εx、εy、εz三个正应变的相对误差分别仅为2.02%、1.66%和0.74%,图6(b)中叶身748号节点γxy、γyz、γxz三个剪应变的相对误差分别仅为1.05%、1.17%和3.02%;图7(a)中叶根1100号节点εx、εy、εz三个正应变的相对误差分别仅为0.62%、0.64%和0.64%,图7(b)中叶根1100号节点γxy、γyz、γxz三个剪应变的相对误差分别仅为0.67%、1.42%和1.57%。因此,本发明提供的一种旋转叶片动应变场测量方法,可以高精度地重构叶片动应变场。
本发明提供的方法仅利用极少测点便可实现旋转叶片整体动应变场的测量,不仅可实现叶片表面动应变的测量,还可以实现叶片内部节点正应变、剪应变的测量,计算过程简单,测量精度高,易于在线测量,节约大量应变片。以上所述仅为本发明的较佳实施例而已,可应用在航空发动机、燃气轮机、汽轮机等旋转机械的风扇/压气机/涡轮叶片振动测试中,并不用以限制本发明。
在另一个实施例中,方法包括以下步骤:
1)提取叶片三维有限元模型的模态参数;
2)确定叶片应变片测点数目、位置、方向;
3)构造有限测点动应变与全场动应变的转换矩阵;
4)利用有限应变片获取旋转叶片有限位置动应变;
5)基于模态降阶与扩展理论实现旋转叶片任意时刻、任意位置、任意方向的动应变测量。
进一步地,步骤1)建立旋转叶片的三维有限元模型,通过模态分析提取前nm阶模态参数:模态频率fi和大小为ndof×1的应变模态振型ψi;构造旋转叶片全场应变模态振型矩阵大小为2ndof×nm;i表示模态阶次,ndof表示叶片有限元模型的自由度数目,则ndof=3nn,nn表示叶片有限元模型节点的数目;每个节点的应变包含3个正应变εx、εy、εz与3个剪应变γxy、γyz、γxz共6个应变分量,即每个节点有6个应变模态振型。
进一步地,步骤2)旋转叶片应变片测点数目nd不得小于关注的模态数目nm,即nd≥nm;
从测点选择矩阵ψp中随机选择nd个测点;构造大小为nd×nm的测点应变模态振型矩阵ψd,并计算其矩阵条件数κ;这个随机过程重复R次,并从中选择条件数κ最小时的测点布局方案。
进一步地,步骤5)基于模态降阶与扩展理论实现旋转叶片任意时刻、任意位置、任意方向的动应变S(t)测量:旋转叶片t时刻叶片表面和内部所有节点应变经由公式S(t)=Ts(t)计算得出,
应变
所述应变s(t)包括正应变和剪应变。
本发明的另一方面还提供了一种旋转叶片动应变场测量方法用系统,包括:
模态分析模块:利用商业有限元分析软件对叶片的三维有限元模型进行模态分析,用于获取旋转叶片前nm阶模态频率fi和应变模态振型ψi;用于构造应变模态振型矩阵大小为2ndof×nm;ndof表示叶片有限元模型的自由度数目。
测点优选模块:用于确定旋转叶片应变片测点数目nd不得小于关注的模态数目nm,即nd≥nm;用于构造关于应变模态振型的测点选择矩阵大小为2nc×nm;其中nc表示叶片表面有限元网格节点数目;2表示每个节点有2个方向可供安装应变片;也就是说,提取叶片可供安装应变片测点的应变模态振型;从测点选择矩阵ψp中随机选择nd个测点;构造大小为nd×nm的测点应变模态振型矩阵ψd,并计算其矩阵条件数κ;这个随机过程重复R次,并从中选择条件数κ最小时的测点布局方案。
动应变场重构模块:用于计算旋转叶片任意时刻、任意位置、任意方向的动应变S(t):旋转叶片t时刻叶片表面和内部所有节点应变经由公式S(t)=Ts(t)计算得出,
应变
所述应变S(t)包括正应变和剪应变。
一种实施所述方法的测量系统包括,
多个应变片,其布置在旋转叶片上;
计算单元,其连接所述动应变测量模块,所述计算单元包括,
测点优选模块,其配置成优化布置在所述旋转叶片上的应变片的测点数目、位置和方向,其中,基于旋转叶片可安装应变片测点的应变模态振型构造旋转叶片全场应变模态振型矩阵从测点选择矩阵ψp中随机选择nd个测点构造大小为nd×nm的测点应变模态振型矩阵ψd并计算其矩阵条件数κ,随机过程重复R次并从中选择矩阵条件数κ最小时的测点布置,
转换矩阵计算模块,其配置成构造测点动应变与全场动应变的转换矩阵,
动应变场重构模块,用于计算旋转叶片t时刻叶片表面和内部所有节点应变经由公式s(t)=Ts(t)计算得出,
应变
所述应变s(t)包括正应变和剪应变。
在一个实施例中,测量系统还包括显示单元和无线通信设备,无线通信设备包括4G/GPRS或互联网通信模块。
在一个实施例中,模态分析模块、测点优选模块、转换矩阵计算模块或动应变场重构模块为通用处理器、数字信号处理器、专用集成电路ASIC或现场可编程门阵列FPGA,
在一个实施例中模态分析模块、测点优选模块、转换矩阵计算模块或动应变场重构模块包括存储器,所述存储器包括一个或多个只读存储器ROM、随机存取存储器RAM、快闪存储器或电子可擦除可编程只读存储器EEPROM。
尽管以上结合附图对本发明的实施方案进行了描述,但本发明并不局限于上述的具体实施方案和应用领域,上述的具体实施方案仅仅是示意性的、指导性的,而不是限制性的。本领域的普通技术人员在本说明书的启示下和在不脱离本发明权利要求所保护的范围的情况下,还可以做出很多种的形式,这些均属于本发明保护之列。
Claims (5)
1.一种旋转叶片动应变场测量方法,所述方法包括以下步骤:
第一步骤(S1)中,建立待测量旋转叶片的三维有限元模型,提取所述三维有限元模型的模态参数;
第二步骤(S2)中,在所述旋转叶片上布置多个应变片,确定所述应变片测点数目、位置和方向;
第三步骤(S3)中,构造应变片测点动应变与全场动应变的转换矩阵;
第四步骤(S4)中,基于所述应变片测量所述旋转叶片相应位置动应变;
第五步骤(S5)中,所述动应变基于所述转换矩阵经由模态处理得到所述旋转叶片任意时刻、任意位置及任意方向的动应变;
其中,第一步骤(S1)中,通过模态分析提取所述三维有限元模型前nm阶模态参数、模态频率fi和大小为ndof×1的应变模态振型ψi,构造旋转叶片全场应变模态振型矩阵大小为2ndof×nm,其中,nm表示模态数目,i表示模态阶次,ndof表示旋转叶片有限元模型的自由度数目,ndof=3nn,nn表示旋转叶片有限元模型节点;
第二步骤(S2)中,基于旋转叶片应变模态振型构造应变模态振型的测点选择矩阵大小为2nc×nm,其中nc表示旋转叶片表面有限元网格节点数目;2表示每个节点有2个方向可安装应变片;从测点选择矩阵ψp中随机选择nd个测点构造大小为nd×nm的测点应变模态振型矩阵ψd,并计算其矩阵条件数κ;随机过程重复R次并从中选择矩阵条件数κ最小时的测点布置;
第五步骤(S5)中,基于转换矩阵T和应变实测的位移信号s(t),旋转叶片t时刻叶片表面和内部所有节点应变经由公式S(t)=Ts(t)计算得出,
应变
2.根据权利要求1所述的方法,其中,第一步骤(S1)中,每个旋转叶片有限元模型节点的应变包含3个正应变εx、εv、εz与3个剪应变γxy、γvz、γxz共6个应变分量。
3.根据权利要求1所述的方法,其中,第二步骤(S2)中,旋转叶片应变片测点数目nd大于等于模态数目nm。
4.根据权利要求1所述的方法,其中,第五步骤(S5)中,所述模态处理基于所述应变模态振型模态降阶与扩展。
5.一种实施权利要求1-4中任一项所述方法的测量系统,所述测量系统包括,
多个应变片,其布置在旋转叶片上;
计算单元,其连接所述动应变测量模块,所述计算单元包括,
测点优选模块,其配置成优化布置在所述旋转叶片上的应变片的测点数目、位置和方向,其中,基于旋转叶片可安装应变片测点的应变模态振型构造旋转叶片全场应变模态振型矩阵从测点选择矩阵ψp中随机选择nd个测点构造大小为nd×nm的测点应变模态振型矩阵ψd并计算其矩阵条件数κ,随机过程重复R次并从中选择矩阵条件数κ最小时的测点布置,
转换矩阵计算模块,其配置成构造测点动应变与全场动应变的转换矩阵,
动应变场重构模块,用于计算旋转叶片t时刻叶片表面和内部所有节点应变,经由公式S(t)=Ts(t)计算得出,
应变
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910225544.0A CN109883389B (zh) | 2019-03-22 | 2019-03-22 | 一种旋转叶片动应变场测量方法及其系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910225544.0A CN109883389B (zh) | 2019-03-22 | 2019-03-22 | 一种旋转叶片动应变场测量方法及其系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109883389A CN109883389A (zh) | 2019-06-14 |
CN109883389B true CN109883389B (zh) | 2020-08-25 |
Family
ID=66933886
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910225544.0A Active CN109883389B (zh) | 2019-03-22 | 2019-03-22 | 一种旋转叶片动应变场测量方法及其系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109883389B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110426208B (zh) * | 2019-08-13 | 2021-06-22 | 中国航发贵阳发动机设计研究所 | 一种高温应变计安装工艺可靠性验证装置 |
CN111060271A (zh) * | 2019-10-29 | 2020-04-24 | 北京动力机械研究所 | 一种小型涡轮发动机涡轮叶片振动应力动测试验方法 |
CN111238961B (zh) * | 2020-01-16 | 2020-12-08 | 西安交通大学 | 基于多粒子追踪技术的软物质非连续应变场测量方法 |
CN111325718B (zh) * | 2020-01-23 | 2023-07-25 | 深圳大学 | 应变模态分析方法及相关装置 |
CN111239432B (zh) * | 2020-02-10 | 2021-05-28 | 西安交通大学 | 基于叶端定时的叶片到达时间测量方法及其采集系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2136189A2 (en) * | 2008-06-18 | 2009-12-23 | Rolls-Royce plc | Method for analysing vibration in rotor blades |
EP2818647A1 (en) * | 2013-05-30 | 2014-12-31 | Rolls-Royce plc | Blade tip timing system and method of analysing blade tip displacement |
CN104501728A (zh) * | 2014-12-12 | 2015-04-08 | 天津大学 | 一种基于全光纤叶尖定时的叶尖间隙测量方法 |
CN104697798A (zh) * | 2015-03-05 | 2015-06-10 | 西安交通大学 | 航空发动机叶尖间隙三维表征与光纤动态检测系统及方法 |
CN106529049A (zh) * | 2016-11-18 | 2017-03-22 | 南京航空航天大学 | 一种用于模型修正的应变模态相关性评价方法 |
CN108827567A (zh) * | 2018-03-17 | 2018-11-16 | 北京化工大学 | 一种传感器安装角度连续可调的叶片振动监测试验装置 |
CN109101769A (zh) * | 2018-09-20 | 2018-12-28 | 中国人民解放军国防科技大学 | 一种基于压缩感知的叶端定时传感器数目确定方法 |
-
2019
- 2019-03-22 CN CN201910225544.0A patent/CN109883389B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2136189A2 (en) * | 2008-06-18 | 2009-12-23 | Rolls-Royce plc | Method for analysing vibration in rotor blades |
EP2818647A1 (en) * | 2013-05-30 | 2014-12-31 | Rolls-Royce plc | Blade tip timing system and method of analysing blade tip displacement |
CN104501728A (zh) * | 2014-12-12 | 2015-04-08 | 天津大学 | 一种基于全光纤叶尖定时的叶尖间隙测量方法 |
CN104697798A (zh) * | 2015-03-05 | 2015-06-10 | 西安交通大学 | 航空发动机叶尖间隙三维表征与光纤动态检测系统及方法 |
CN106529049A (zh) * | 2016-11-18 | 2017-03-22 | 南京航空航天大学 | 一种用于模型修正的应变模态相关性评价方法 |
CN108827567A (zh) * | 2018-03-17 | 2018-11-16 | 北京化工大学 | 一种传感器安装角度连续可调的叶片振动监测试验装置 |
CN109101769A (zh) * | 2018-09-20 | 2018-12-28 | 中国人民解放军国防科技大学 | 一种基于压缩感知的叶端定时传感器数目确定方法 |
Non-Patent Citations (2)
Title |
---|
利用动态位移信息的全场动应变;陈帅;《万方学位论文》;20181202;2.3 模态缩减/扩展 * |
旋转叶片振动性能参数测试技术研究;欧阳涛;《中国优秀学位论文全文数据库工程科技II辑》;20090715;3.4 叶片固有频率和模态分析,4.1 叶尖定时测振原理,4.2 叶尖定时测振系统构成,5.1 某型号设备旋转叶片振动测试实验 * |
Also Published As
Publication number | Publication date |
---|---|
CN109883389A (zh) | 2019-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110608710B (zh) | 一种基于叶端定时的转子叶片动应变场测量方法及其系统 | |
CN109883380B (zh) | 一种基于叶端定时的转子叶片位移场测量方法及其系统 | |
CN109870134B (zh) | 一种旋转叶片非接触式动应变场测量方法及其系统 | |
CN109883389B (zh) | 一种旋转叶片动应变场测量方法及其系统 | |
CN109885976B (zh) | 一种旋转叶片位移场反演重构方法及其系统 | |
CN108895974B (zh) | 一种结构变形光纤监测与重构方法及系统 | |
CN110375690B (zh) | 一种旋转叶片非接触式位移场测量方法及其系统 | |
CN110032776B (zh) | 基于模态解耦的叶片动应变场与位移场重构方法及其系统 | |
EP2904368B1 (en) | Turbine blade fatigue life analysis using non-contact measurement and dynamical response reconstruction techniques | |
CN111507042A (zh) | 基于叶端定时的旋转叶片动应力测量方法及其系统 | |
Kim et al. | Shape estimation with distributed fiber Bragg grating sensors for rotating structures | |
CN111507043A (zh) | 一种基于叶端定时的转子叶片动应力场测量方法及其系统 | |
CN101122541A (zh) | 汽轮机叶片振动试验方法及装置 | |
CN111563340A (zh) | 一种转子叶片动应力重构方法及其系统 | |
CN104296924B (zh) | 动力总成惯性参数测试方法 | |
Ao et al. | Blade dynamic strain non-intrusive measurement using L1/2-norm regularization and transmissibility | |
CN111563342A (zh) | 一种转子叶片的应变片测点布局优化方法 | |
CN104462862A (zh) | 一种基于三次b样条尺度函数的机械结构动载荷识别方法 | |
CN112464535B (zh) | 一种转子叶片动应变测量数据一致性评估方法 | |
CN113029481A (zh) | 一种针对叶片扭转振动的测量方法 | |
CN109883379B (zh) | 基于模态振型的叶片位移应变测量方法 | |
CN115062410A (zh) | 单传感器转子叶片多模态应变场重构方法及测量系统 | |
CN115081271B (zh) | 基于数字模拟器的叶端定时系统校核方法及校核系统 | |
US12007309B2 (en) | Non-contact dynamic strain field measuring method and system for rotating blade | |
JP3114358B2 (ja) | 構造物の振動試験装置及び方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |