CN107991080A - 一种基于非接触测振与仿真计算的高频叶片模态分析方法 - Google Patents
一种基于非接触测振与仿真计算的高频叶片模态分析方法 Download PDFInfo
- Publication number
- CN107991080A CN107991080A CN201711270610.3A CN201711270610A CN107991080A CN 107991080 A CN107991080 A CN 107991080A CN 201711270610 A CN201711270610 A CN 201711270610A CN 107991080 A CN107991080 A CN 107991080A
- Authority
- CN
- China
- Prior art keywords
- blade
- modal
- frequency
- high frequency
- vibration measuring
- 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.)
- Withdrawn
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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]
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于非接触测振与仿真计算的高频叶片模态分析方法,可满足1000Hz‑15000Hz高频范围内的模态分析,包括以下步骤:(1)采用有限元计算方法,计算自由模态与振型;(2)根据自由模态计算结果,选择测点,利用激光多普勒测振技术,测试叶片原点频响函数,获取叶片真实的固有频率;(3)修正仿真模型的弹性模量、泊松比参数,使仿真模型的计算结果与真实固有频率一致;对修正后的仿真模型施加工作边界条件,计算工作状态下的叶片模态;本发明可针对高频叶片,获得精确的模态频率与振型,为发动机叶片共振、疲劳失效等故障,提供可靠的分析基准。
Description
技术领域
本发明涉及一种高频叶片的模态分析方法,具体的说是一种基于非接触测振与仿真计算的高频叶片模态分析方法。
背景技术
离心压气机叶盘与向心涡轮叶盘的叶片共振疲劳、掉块等故障是微小型航空发动机的典型故障,从共振角度出发,为研究并解决该故障,对叶片进行精准的模态分析是十分必要的,离心压气机叶盘和向心涡轮叶盘因结构小,型面复杂,固有频率高等特征,很难通过常规的模态测试手段对其完成精确的模态分析;同时,因参数误差、结构误差、离散误差的客观存在,仿真模型的计算精度也有限,因此,很难利用上述方法获得精确的叶片模态。
发明内容
本发明所要解决的技术问题是,针对以上现有的技术缺陷,提出一种基于非接触测振与仿真计算的高频叶片模态分析方法,该方法是针对高频叶片,基于非接触式频响函数测试,结合模态仿真计算,通过修正弹性模量和泊松比参数,获取精确的叶片模态,为发动机叶片共振、疲劳时效等故障,提供可靠的分析基准,该方法克服了因叶片型面复杂、固有频率高等带来的测试难度高、精度差的缺点,也避免了因结构、参数、离散引入的仿真计算误差。
为了解决上述技术问题,本发明提供一种基于非接触测振与仿真计算的高频叶片模态分析方法,该方法是针对高频叶片,基于非接触式频响函数测试,结合模态仿真计算,通过修正弹性模量和泊松比参数,获取精确的叶片模态,为发动机叶片共振、疲劳时效等故障,提供可靠的分析基准,具体包括以下步骤:
(1) 采用有限元计算方法,计算高频叶片的自由模态,得到自由模态的固有频率与模态振型;
(2) 根据计算的模态振型,选择测点,利用激光多普勒测振方式测试叶片原点频响函数,获取叶片真实的固有频率;
(3) 修正仿真模型的弹性模量、泊松比参数,使得仿真模型的计算结果与真实固有频率保持一致,对修正后的仿真模型施加工作边界条件,计算得到工作状态下的叶片模态。
本发明进一步限定的技术方案是:
前述基于非接触测振与仿真计算的高频叶片模态分析方法中,步骤(1)中高频叶片是包括微小型航空发动机的离心压气机叶盘、向心涡轮叶盘,其叶片模态频率在1000Hz-15000Hz高频范围内。
前述基于非接触测振与仿真计算的高频叶片模态分析方法中,步骤(2)中选择测点时测点在同一个叶片上,根据叶片模态振型,选择振幅最大的点作为测试点。
前述基于非接触测振与仿真计算的高频叶片模态分析方法中,步骤(2)中测试频响函数的激励方式包括:声场激励、电磁激励和力锤激励,其中:对于钢件叶片选择电磁激励和力锤激励;对于非钢件选择声场激励和力锤激励。
前述基于非接触测振与仿真计算的高频叶片模态分析方法中,步骤(2)中利用激光多普勒测振方式测试叶片原点频响函数时,测试仪器与被测点的距离不小于0.5米,测试激光束与被测平面法向的角度不大于20°。
前述基于非接触测振与仿真计算的高频叶片模态分析方法中,步骤(3)中修正的目标是被测件的固有频率计算值与步骤(2)中的测试值之间误差在0.3%以内,修正的参数变量包括被测件的弹性模量、泊松比参数。
前述基于非接触测振与仿真计算的高频叶片模态分析方法中,步骤(3)施加被测件的工作边界条件具体为:对于离心压气机叶盘,施加发动机实际的工作转速对应的离心载荷场;对于向心涡轮叶盘,施加发动机实际的工作转速对应的离心载荷场和温度场。
本发明的有益效果是:
由于采用了非接触式激光多普勒测振技术,避免了因传感器附加质量引入的误差,测试的固有频率值精确无误;基于几何建模准确的前提下,叶片的自由模态振型是准确无误的;
本发明充分利用了二者方法的各自优势,以测试固有频率为基准,修正仿真模型的弹性模量、泊松比参数,获取精确的叶片模态频率与振型,基于该模型下的模态分析,频率与振型准确无误,这将对叶片共振、疲劳、掉块等故障产生直接的工程指导价值。
附图说明
图1为本发明的流程图;
图2 离心叶轮的模态振型;
图3 离心叶轮的自由模态测试图;
图4 离心叶轮的叶片频响函数。
具体实施方式
下面对本发明做进一步的详细说明:
实施例1
本实施例提供一种基于非接触测振与仿真计算的高频叶片模态分析方法,该方法是针对高频叶片,基于非接触式频响函数测试,结合模态仿真计算,通过修正弹性模量和泊松比参数,获取精确的叶片模态,为发动机叶片共振、疲劳时效等故障,提供可靠的分析基准,具体包括以下步骤,参照附图1:
(1) 采用有限元计算方法,计算高频叶片的自由模态,得到自由模态的固有频率与模态振型;
(2) 根据计算的模态振型,选择测点,利用激光多普勒测振方式测试叶片原点频响函数,获取叶片真实的固有频率;
(3) 修正仿真模型的弹性模量、泊松比参数,使得仿真模型的计算结果与真实固有频率保持一致,对修正后的仿真模型施加工作边界条件,计算得到工作状态下的叶片模态。
参照附图2,案例对象为小型涡喷发动机的离心压气机叶轮,对其建立几何模型,设置材料参数,采用有限元计算方法,计算自由模态,得到固有频率与振型,根据固有频率,选取频率测试范围1000Hz-15000Hz。
参照附图3,将叶轮自由悬挂,根据模态振型,选取叶尖区域作为测量点,利用激光多普勒测振技术,采用力锤激励,测试叶片原点频响函数,如附图4所示,测试时保证激光光线路径能够达到测点位置。
分析测试频率与仿真分析频率的差异,并以测试频率为基准,修正计算模型的弹性模量和泊松比,得到准确的高频叶片模型,对高频叶片模型进行工作状态(即:旋转状态)边界条件下的模态计算,获取工作状态下,高频叶片的固有频率与振型。
在本实施例中,步骤(2)中利用激光多普勒测振方式测试叶片原点频响函数时,测试仪器与被测点的距离不小于0.5米,测试激光束与被测平面法向的角度不大于20°;步骤(3)中修正的目标是被测件的固有频率计算值与步骤(2)中的测试值之间误差在0.3%以内,修正的参数变量包括被测件的弹性模量、泊松比参数;步骤(3)施加被测件的工作边界条件具体为:对于离心压气机叶盘,施加发动机实际的工作转速对应的离心载荷场;对于向心涡轮叶盘,施加发动机实际的工作转速对应的离心载荷场和温度场。
这样通过案例的技术方案,该离心叶轮的叶片频率便能够准确获得,这对叶片的共振、疲劳、掉块等故障排除的工作开展,具有切实的工程指导意义。本发明具有分析结论精确,操作相对简单,便于测试等特点。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (7)
1.一种基于非接触测振与仿真计算的高频叶片模态分析方法,其特征在于,该方法是针对高频叶片,基于非接触式频响函数测试,结合模态仿真计算,通过修正弹性模量和泊松比参数,获取精确的叶片模态,为发动机叶片共振、疲劳时效等故障,提供可靠的分析基准,具体包括以下步骤:
(1) 采用有限元计算方法,计算高频叶片的自由模态,得到自由模态的固有频率与模态振型;
(2) 根据计算的模态振型,选择测点,利用激光多普勒测振方式测试叶片原点频响函数,获取叶片真实的固有频率;
(3) 修正仿真模型的弹性模量、泊松比参数,使得仿真模型的计算结果与真实固有频率保持一致,对修正后的仿真模型施加工作边界条件,计算得到工作状态下的叶片模态。
2.根据权利要求1所述的基于非接触测振与仿真计算的高频叶片模态分析方法,其特征在于,所述步骤(1)中高频叶片是包括微小型航空发动机的离心压气机叶盘、向心涡轮叶盘,其叶片模态频率在1000Hz-15000Hz高频范围内。
3.根据权利要求1所述的基于非接触测振与仿真计算的高频叶片模态分析方法,其特征在于,所述步骤(2)中选择测点时测点在同一个叶片上,根据叶片模态振型,选择振幅最大的点作为测试点。
4.根据权利要求1所述的基于非接触测振与仿真计算的高频叶片模态分析方法,其特征在于,所述步骤(2)中测试频响函数的激励方式包括:声场激励、电磁激励和力锤激励,其中:对于钢件叶片选择电磁激励和力锤激励;对于非钢件选择声场激励和力锤激励。
5.根据权利要求1所述的基于非接触测振与仿真计算的高频叶片模态分析方法,其特征在于,所述步骤(2)中利用激光多普勒测振方式测试叶片原点频响函数时,测试仪器与被测点的距离不小于0.5米,测试激光束与被测平面法向的角度不大于20°。
6.根据权利要求1所述的基于非接触测振与仿真计算的高频叶片模态分析方法,其特征在于,所述步骤(3)中修正的目标是被测件的固有频率计算值与步骤(2)中的测试值之间误差在0.3%以内,修正的参数变量包括被测件的弹性模量、泊松比参数。
7.根据权利要求2所述的基于非接触测振与仿真计算的高频叶片模态分析方法,其特征在于,所述步骤(3)施加被测件的工作边界条件具体为:对于离心压气机叶盘,施加发动机实际的工作转速对应的离心载荷场;对于向心涡轮叶盘,施加发动机实际的工作转速对应的离心载荷场和温度场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711270610.3A CN107991080A (zh) | 2017-12-05 | 2017-12-05 | 一种基于非接触测振与仿真计算的高频叶片模态分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711270610.3A CN107991080A (zh) | 2017-12-05 | 2017-12-05 | 一种基于非接触测振与仿真计算的高频叶片模态分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107991080A true CN107991080A (zh) | 2018-05-04 |
Family
ID=62036055
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711270610.3A Withdrawn CN107991080A (zh) | 2017-12-05 | 2017-12-05 | 一种基于非接触测振与仿真计算的高频叶片模态分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107991080A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109033574A (zh) * | 2018-07-10 | 2018-12-18 | 西安工程大学 | 考虑运行温度变化的输电导线模态分析方法 |
CN109060284A (zh) * | 2018-08-07 | 2018-12-21 | 广东工业大学 | 一种基于dic技术的试验模态分析方法 |
CN109883379A (zh) * | 2019-03-22 | 2019-06-14 | 西安交通大学 | 基于模态振型的叶片位移应变测量方法 |
CN109932151A (zh) * | 2019-03-28 | 2019-06-25 | 东北大学 | 一种行波激励作用下整体叶盘节径运动测试装置及方法 |
CN110134990A (zh) * | 2019-04-03 | 2019-08-16 | 南京航空航天大学 | 航空发动机叶片前缘危险部位的检测方法 |
CN110375690A (zh) * | 2019-03-22 | 2019-10-25 | 西安交通大学 | 一种旋转叶片非接触式位移场测量方法及其系统 |
CN111024214A (zh) * | 2019-12-25 | 2020-04-17 | 华中科技大学 | 一种实时获取声共振混合机运行过程中固有频率的方法 |
CN111850442A (zh) * | 2020-07-06 | 2020-10-30 | 中国人民解放军空军工程大学 | 一种钛合金整体叶盘叶片预防高阶振型诱发叶尖掉块的强化方法 |
CN112444365A (zh) * | 2020-11-30 | 2021-03-05 | 哈尔滨工业大学 | 一种基于力锤摆动法及激光多普勒法的卫星太阳翼基板展开低频模态测试方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009085788A (ja) * | 2007-09-28 | 2009-04-23 | Takenaka Komuten Co Ltd | スラブの振動同定方法、制振装置、制振装置配置方法、建築床構造、及び振動測定装置 |
CN103835979A (zh) * | 2012-11-27 | 2014-06-04 | 重庆长安汽车股份有限公司 | 一种非接触式的风扇模态测试系统及方法 |
CN104776963A (zh) * | 2015-04-14 | 2015-07-15 | 北京强度环境研究所 | 一种声激励非接触式模态试验系统及方法 |
JP2016024134A (ja) * | 2014-07-23 | 2016-02-08 | 三井精機工業株式会社 | モーダル解析支援装置及び同様の支援機構を備えた実稼働解析支援装置 |
CN105912772A (zh) * | 2016-04-08 | 2016-08-31 | 北京宇航系统工程研究所 | 一种集频率、振型和频响函数为一体的模型修正方法 |
CN106529055A (zh) * | 2016-11-18 | 2017-03-22 | 南京航空航天大学 | 一种基于应变模态振型相关性的模型修正方法 |
CN107169172A (zh) * | 2017-04-18 | 2017-09-15 | 上海交通大学 | 一种复杂装配体的振动分析仿真方法 |
-
2017
- 2017-12-05 CN CN201711270610.3A patent/CN107991080A/zh not_active Withdrawn
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009085788A (ja) * | 2007-09-28 | 2009-04-23 | Takenaka Komuten Co Ltd | スラブの振動同定方法、制振装置、制振装置配置方法、建築床構造、及び振動測定装置 |
CN103835979A (zh) * | 2012-11-27 | 2014-06-04 | 重庆长安汽车股份有限公司 | 一种非接触式的风扇模态测试系统及方法 |
JP2016024134A (ja) * | 2014-07-23 | 2016-02-08 | 三井精機工業株式会社 | モーダル解析支援装置及び同様の支援機構を備えた実稼働解析支援装置 |
CN104776963A (zh) * | 2015-04-14 | 2015-07-15 | 北京强度环境研究所 | 一种声激励非接触式模态试验系统及方法 |
CN105912772A (zh) * | 2016-04-08 | 2016-08-31 | 北京宇航系统工程研究所 | 一种集频率、振型和频响函数为一体的模型修正方法 |
CN106529055A (zh) * | 2016-11-18 | 2017-03-22 | 南京航空航天大学 | 一种基于应变模态振型相关性的模型修正方法 |
CN107169172A (zh) * | 2017-04-18 | 2017-09-15 | 上海交通大学 | 一种复杂装配体的振动分析仿真方法 |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109033574A (zh) * | 2018-07-10 | 2018-12-18 | 西安工程大学 | 考虑运行温度变化的输电导线模态分析方法 |
CN109033574B (zh) * | 2018-07-10 | 2023-06-06 | 西安工程大学 | 考虑运行温度变化的输电导线模态分析方法 |
CN109060284A (zh) * | 2018-08-07 | 2018-12-21 | 广东工业大学 | 一种基于dic技术的试验模态分析方法 |
CN110375690B (zh) * | 2019-03-22 | 2021-02-26 | 西安交通大学 | 一种旋转叶片非接触式位移场测量方法及其系统 |
CN109883379A (zh) * | 2019-03-22 | 2019-06-14 | 西安交通大学 | 基于模态振型的叶片位移应变测量方法 |
CN110375690A (zh) * | 2019-03-22 | 2019-10-25 | 西安交通大学 | 一种旋转叶片非接触式位移场测量方法及其系统 |
CN109883379B (zh) * | 2019-03-22 | 2020-07-10 | 西安交通大学 | 基于模态振型的叶片位移应变测量方法 |
CN109932151A (zh) * | 2019-03-28 | 2019-06-25 | 东北大学 | 一种行波激励作用下整体叶盘节径运动测试装置及方法 |
CN110134990A (zh) * | 2019-04-03 | 2019-08-16 | 南京航空航天大学 | 航空发动机叶片前缘危险部位的检测方法 |
CN110134990B (zh) * | 2019-04-03 | 2020-09-25 | 南京航空航天大学 | 航空发动机叶片前缘危险部位的检测方法 |
CN111024214A (zh) * | 2019-12-25 | 2020-04-17 | 华中科技大学 | 一种实时获取声共振混合机运行过程中固有频率的方法 |
CN111850442B (zh) * | 2020-07-06 | 2022-06-10 | 中国人民解放军空军工程大学 | 一种钛合金整体叶盘叶片预防高阶振型诱发叶尖掉块的强化方法 |
CN111850442A (zh) * | 2020-07-06 | 2020-10-30 | 中国人民解放军空军工程大学 | 一种钛合金整体叶盘叶片预防高阶振型诱发叶尖掉块的强化方法 |
CN112444365A (zh) * | 2020-11-30 | 2021-03-05 | 哈尔滨工业大学 | 一种基于力锤摆动法及激光多普勒法的卫星太阳翼基板展开低频模态测试方法 |
CN112444365B (zh) * | 2020-11-30 | 2023-08-29 | 哈尔滨工业大学 | 一种卫星太阳翼基板展开低频模态测试方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107991080A (zh) | 一种基于非接触测振与仿真计算的高频叶片模态分析方法 | |
Rigosi et al. | Synchronous vibration parameters identification by tip timing measurements | |
Ewins | The effects of detuning upon the forced vibrations of bladed disks | |
Chen et al. | New step to improve the accuracy of blade tip timing method without once per revolution | |
CN109238718B (zh) | 一种发动机风扇叶片高周疲劳性能测试方法 | |
CN105352713B (zh) | 通过检测叶片静位移实现检测叶片疲劳裂纹的方法 | |
Zhao et al. | A transient characteristic-based balancing method of rotor system without trail weights | |
Mangardich et al. | A fracture mechanics based approach for the fretting fatigue of aircraft engine fan dovetail attachments | |
CN109063287A (zh) | 一种离心叶轮中心孔裂纹扩展特征模拟件设计方法 | |
CN109029884A (zh) | 一种汽车悬臂结构件振动疲劳分析的方法 | |
Hemberger et al. | Investigations on maximum amplitude amplification factor of real mistuned bladed structures | |
CN110319922A (zh) | 基于叶尖定时原理的叶片扭转振动位移测量方法 | |
Wang et al. | A novel method to improve the precision of BTT under rapid speed fluctuation conditions | |
Wu et al. | Blade crack detection using blade tip timing | |
Mohamed et al. | The determination of steady-state movements using blade tip timing data | |
CN115114721A (zh) | 基于非接触测量的叶片多模态最大应力预测方法及系统 | |
Stadtmüller et al. | Experimental and numerical investigation of wake-induced transition on a highly loaded LP turbine at low Reynolds numbers | |
Zhang et al. | Five dimensional movement measurement method for rotating blade based on blade tip timing measuring point position tracking | |
Kaszynski et al. | Accurate blade tip timing limits through geometry mistuning modeling | |
Kharyton et al. | Simulation of tip-timing measurements of a cracked bladed disk forced response | |
Hemberger et al. | Mistuning modeling and its validation for small turbine wheels | |
Przysowa | The analysis of synchronous blade vibration using linear sine fitting | |
Kielb et al. | Experimental and numerical study of forced response in a full-scale rotating turbine | |
Gillaugh et al. | Mistuning evaluation comparison via as-manufactured models, traveling wave excitation, and compressor rigs | |
Hartung et al. | A Practical Approach for Evaluation of Equivalent Linear Damping From Measurements of Mistuned and/or Non-Linear Stages and Forced Response Validation |
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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20180504 |
|
WW01 | Invention patent application withdrawn after publication |