CN111523182B - 一种基于流固耦合的叶片振动响应分析方法 - Google Patents

一种基于流固耦合的叶片振动响应分析方法 Download PDF

Info

Publication number
CN111523182B
CN111523182B CN202010157238.0A CN202010157238A CN111523182B CN 111523182 B CN111523182 B CN 111523182B CN 202010157238 A CN202010157238 A CN 202010157238A CN 111523182 B CN111523182 B CN 111523182B
Authority
CN
China
Prior art keywords
blade
fluid
analysis
solid coupling
pneumatic
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
Application number
CN202010157238.0A
Other languages
English (en)
Other versions
CN111523182A (zh
Inventor
漆文凯
许正华
田少杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202010157238.0A priority Critical patent/CN111523182B/zh
Publication of CN111523182A publication Critical patent/CN111523182A/zh
Application granted granted Critical
Publication of CN111523182B publication Critical patent/CN111523182B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

本发明公开了一种基于流固耦合的叶片振动响应分析方法;首先,采用非线性谐波法对叶片流道内的流场进行非定常分析,得到叶片表面的非定常载荷分布,采用傅里叶变换提取关键频率所对应的载荷,获得叶片所受气动激振力;然后,对叶片及其流道内的流场进行流固耦合分析,分析叶片振动时气动力对叶片所做的气动功,根据能量等效原理获得叶片在流场中振动的各阶模态气动阻尼比;最后,将气动激振力与各阶模态气动阻尼比加载到叶片上,采用模态叠加法对叶片进行振动响应分析;该分析方法可以高效、准确地分析叶片在非定常气动力作用下的振动响应,为发动机的安全保障和设计工作提供支持。

Description

一种基于流固耦合的叶片振动响应分析方法
技术领域
本发明涉及航空系统技术领域,主要涉及一种基于流固耦合的叶片振动响应分析方法。
背景技术
随着现代飞机的发展,人们对航空发动机有着更高的要求,一方面要求其推力大、推重比、高和耗油率低,另一方面要求其可靠性高、寿命长和重量轻等特点。目前各种新技术的使用,提高了压气机的增压比和效率,但使得压气机叶片表面气动载荷大大增加,叶片振动问题严重,进而引起叶片的高周疲劳失效(HCF)。叶片一旦发生断裂故障,断片可能将其它级叶片打坏甚至打穿机匣,从而严重损坏发动机,甚至造成飞行事故。
为了保证发动机的安全,需要在发动机设计之初对叶片的振动应力进行预估,因此,建立一套准确的叶片振动响应分析方法成为发动机设计工作的关键点之一。但是,发动机叶片的工作环境是复杂的非定常流场,分析叶片的振动响应是一个流固耦合问题,完整的求解固体和流体方程来计算叶片振动响应的精确解是一项计算量巨大、计算成本高昂的工作。将气体力对叶片的作用分解为激励作用和阻尼作用,分别分析叶片在流场中的气动激振力和气动阻尼,并在各个部分中考虑流固耦合的作用,同时保证了计算的精度。
发明内容
发明目的:本发明提供了一种基于流固耦合的叶片振动响应分析方法,将气体力对叶片的作用分解为激励作用和阻尼作用,分别分析叶片在流场中的气动激振力和气动阻尼,并在各个部分中考虑流固耦合的作用,同时保证了计算的精度。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种基于流固耦合的叶片振动响应分析方法,包括气动激振力分析、气动阻尼分析以及叶片振动响应分析三个部分,具体包括以下步骤:
步骤S1、采用非线性谐波法对叶片流道内的流场进行非定常分析,获得叶片表面所受的非定常载荷,通过傅里叶变换对非定常结果进行频谱分析,根据不同频率对应的幅值,根据振幅从大到小顺序提取若干阶关键频率及其对应载荷来模拟原流场,通过流固耦合接口程序,获得叶片表面所受的气动激振力;
步骤S2、基于能量法和弱流固耦合分析法,对叶片以及其流道的流场进行流固耦合分析,对叶片进行模态分析,获得各阶固有频率与振型,编写动网格控制程序使得流体域内的叶型边界面以叶片的各阶固有模态形式做周期运动,通过气动阻尼积分计算程序,获得气动力对叶片所做的气动功,根据能量等效原理获得叶片的各阶模态气动阻尼比;
步骤S3、基于步骤S1获得的叶片表面的气动力激振力和步骤S2获得的叶片各阶模态气动阻尼比,采用模态叠加法对叶片进行振动响应分析。
进一步地,所述步骤S1中对叶片流道内流场进行非定常分析过程中,所述流场非定常分析的对象为整圈流道,流场进口输入条件包括轮盘转速、总温、总压边界条件,出口输入条件包括静压边界条件,并且进口总温、总压和出口静压边界条件都以特定的函数形式随时间变化。
进一步地,所述流场非定常分析采用的计算软件为NUMECA。
进一步地,所述步骤S1中所述的流固耦合接口程序,采用插值法将流场计算得到的气动激振力加载到固体叶片上。
进一步地,所述步骤S1中所述的流固耦合接口程序,采用MATLAB以及ANSYS APDL二次开发语言编写。
进一步地,所述步骤S2中的动网格控制程序,由MATLAB语言编写,将叶片固体域网格各个节点的运动幅值按照振型向量设置后,采用插值法加载到流体域的流固耦合边界面的网格节点上,从而控制流体域的流固耦合边界面运动。
进一步地,步骤S2所述的气动阻尼积分计算程序,由MATLAB语言编写,将叶片振动的一个周期分为100个时间步,每个时间步内各个节点的位移与气动力的乘积即为气动力所做的气动功;将流固耦合边界面上所有节点、所有时间步的气动功累加,得到叶片振动一个周期内气动力对整个叶片所做的气动功。
进一步地,所述步骤S2中计算叶片各阶气动阻尼比时,阶数不低于6阶。
有益效果:本方法具备以下优点:
(1)将气体力的作用分解为激励作用与阻尼作用的,并且在两个部分中都考虑了流固耦合效应的影响,所以本发明在保证求解精度的前提下,减小了求解叶片振动响应过程中的计算量,提高了计算效率。
(2)本发明在求解过程中,无需编写针对特定问题专用的流固耦合求解程序,仅在现有商用计算软件基础上进行分析就可以得到结果,所以本发明有着较低的计算成本。
(3)本发明在求解过程中,多处采用了插值法进行流体域与固体域之间的数据传递,所以本发明对所求解叶片的叶型有更广泛的适用性。
(4)本发明解决了非定常气动力作用下叶片振动响应的准确、高效求解问题,所以本发明具有实际的工程意义,可供工程人员参考,具有可观的应用前景。
附图说明
图1是本发明提供的分析方法流程图;
图2是本发明实施方式中的整圈叶片模型;
图3是本发明实施方式中的气动激振力分布云图;
图4是本发明实施方式中的单一叶片固体模型;
图5是本发明实施方式中的单一叶片流道的流场模型;
图6是本发明实施方式中的叶盆振动应力分布云图;
图7是本发明实施方式中的叶背振动应力分布云图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
针对某型发动机压气机中的叶片模型,建立1.5级三排(前排静子、转子和后排静子)单通道扇区流场计算模型,如图2所示,各排叶片数目如表1所示。
表1叶排叶片数目
首先,采用非线性谐波法对叶片流道内的流场进行非定常分析,获得叶片表面所受的非定常载荷,采用傅里叶变换对非定常结果进行频谱分析,根据不同频率对应的幅值,提取振幅较大的前几阶关键频率和对应载荷来模拟原流场,通过流固耦合接口程序,获得叶片表面所受的气动激振力。
非定常计算时,设置进口总温、进口总压、出口静压以及轮盘转速。采用HarmonicMethod进行非定常计算,以相同进、出口边界条件下的定常计算结果为初场,采用非线性谐波法求解。为提高计算效率,每个扰动的谐波数设置为3。
流体域的计算得到了叶片表面非定常压力分布,编写流固耦合接口程序,由于流体域网格与固体域网格不匹配,采用MATLAB语言编写流固网格节点插值程序,采用ANSYSAPDL二次开发语言,将插值过后的非定常压力完成向气动激振力的转化,各面单元的压力与该单元的面积相乘,得到该单元的气动力,气动力分解,分别加载到面单元各节点,如图3所示。
接着,基于能量法和弱流固耦合分析法,对叶片以及其流道的流场进行流固耦合分析,对叶片进行模态分析,获得各阶固有频率与振型,编写动网格控制程序使得流体域内的叶型边界面以叶片的各阶固有模态形式做周期运动,通过气动阻尼积分计算程序,获得气动力对叶片所做的气动功,根据能量等效原理获得叶片的各阶模态气动阻尼比。
建立单一叶片固体模型以及其流道的流场模型,如图4和图5所示,在固体域对叶片施加根部固支约束,设置轮盘转速,进行模态分析,获得叶片各阶固有频率与振型。采用ANSYS APDL二次开发语言将叶片各阶模态振型导出,获得叶片固体域网格节点的位移,通过MATLAB语言编写插值程序把叶片固体域网格节点位移插值到流体域网格节点中,然后将流体域网格节点的位移写入到CFX的动网格控制程序中,并定义运动频率为叶片固有频率。
非定常计算时,导入动网格控制程序,设置进口总温、进口总压、出口静压以及轮盘转速。设置非定常计算的计算时间步长和总时间步长,为了提高计算的精度,叶片振动一个周期设置为含有100个时间步。计算得到每个时间步内每个节点的位移和气动力,通过MATLAB语言编写计算程序计算叶片振动一个周期内,气动力对整个叶片所做总的气动功。基于能量等效原理,获得前八阶模态的气动阻尼比,如表2所示。
表2叶片各阶模态气动阻尼比
最后,将所得气动力激振力和叶片各阶模态气动阻尼比,采用模态叠加法对叶片进行振动响应分析。
将气动激振力加载到转子叶片上,各阶模态阻尼比设置为计算获得的叶片各阶模态气动阻尼比,通过ANSYS中的模态叠加法进行谐响应计算,叶盆和叶背的振动应力分布云图如图6和图7所示,除去叶片固支带来的在叶片根部应力集中,叶身最大振动应力约为100MPa。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (8)

1.一种基于流固耦合的叶片振动响应分析方法,其特征在于,包括气动激振力分析、气动阻尼分析以及叶片振动响应分析三个部分,具体包括以下步骤:
步骤S1、采用非线性谐波法对叶片流道内的流场进行非定常分析,获得叶片表面所受的非定常载荷,通过傅里叶变换对非定常结果进行频谱分析,根据不同频率对应的幅值,根据振幅从大到小顺序提取若干阶关键频率及其对应载荷来模拟原流场,通过流固耦合接口程序,获得叶片表面所受的气动激振力;
步骤S2、基于能量法和弱流固耦合分析法,对叶片以及其流道的流场进行流固耦合分析,对叶片进行模态分析,获得各阶固有频率与振型,编写动网格控制程序使得流体域内的叶型边界面以叶片的各阶固有模态形式做周期运动,通过气动阻尼积分计算程序,获得气动力对叶片所做的气动功,根据能量等效原理获得叶片的各阶模态气动阻尼比;
步骤S3、基于步骤S1获得的叶片表面的气动激振力和步骤S2获得的叶片各阶模态气动阻尼比,采用模态叠加法对叶片进行振动响应分析。
2.根据权利要求1所述的一种基于流固耦合的叶片振动响应分析方法,其特征在于,所述步骤S1中对叶片流道内流场进行非定常分析过程中,所述流场非定常分析的对象为整圈流道,流场进口输入条件包括轮盘转速、总温、总压边界条件,出口输入条件包括静压边界条件,并且进口总温、总压和出口静压边界条件设定为随时间变化的函数形式。
3.根据权利要求2所述的一种基于流固耦合的叶片振动响应分析方法,其特征在于,所述流场非定常分析采用的计算软件为NUMECA。
4.根据权利要求1所述的一种基于流固耦合的叶片振动响应分析方法,其特征在于,所述步骤S1中所述的流固耦合接口程序,采用插值法将流场计算得到的气动激振力加载到固体叶片上。
5.根据权利要求4所述的一种基于流固耦合的叶片振动响应分析方法,其特征在于,所述步骤S1中所述的流固耦合接口程序,采用MATLAB以及ANSYS APDL二次开发语言编写。
6.根据权利要求1所述的一种基于流固耦合的叶片振动响应分析方法,其特征在于,所述步骤S2中的动网格控制程序,由MATLAB语言编写,将叶片固体域网格各个节点的运动幅值按照振型向量设置后,采用插值法加载到流体域的流固耦合边界面的网格节点上,从而控制流体域的流固耦合边界面运动。
7.根据权利要求1所述的一种基于流固耦合的叶片振动响应分析方法,其特征在于:步骤S2所述的气动阻尼积分计算程序,由MATLAB语言编写,将叶片振动的一个周期分为100个时间步,每个时间步内各个节点的位移与气动力的乘积即为气动力所做的气动功;将流固耦合边界面上所有节点、所有时间步的气动功累加,得到叶片振动一个周期内气动力对整个叶片所做的气动功。
8.根据权利要求1所述的一种基于流固耦合的叶片振动响应分析方法,其特征在于:所述步骤S2中计算叶片各阶气动阻尼比时,阶数不低于6阶。
CN202010157238.0A 2020-03-09 2020-03-09 一种基于流固耦合的叶片振动响应分析方法 Active CN111523182B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010157238.0A CN111523182B (zh) 2020-03-09 2020-03-09 一种基于流固耦合的叶片振动响应分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010157238.0A CN111523182B (zh) 2020-03-09 2020-03-09 一种基于流固耦合的叶片振动响应分析方法

Publications (2)

Publication Number Publication Date
CN111523182A CN111523182A (zh) 2020-08-11
CN111523182B true CN111523182B (zh) 2024-02-13

Family

ID=71901512

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010157238.0A Active CN111523182B (zh) 2020-03-09 2020-03-09 一种基于流固耦合的叶片振动响应分析方法

Country Status (1)

Country Link
CN (1) CN111523182B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112859592B (zh) * 2020-12-29 2022-08-09 中国航空工业集团公司西安飞机设计研究所 一种用于涡桨飞机结构模态频率控制的方法
CN112861288B (zh) * 2021-03-08 2021-11-23 浙江水泵总厂有限公司 一种流体机械校核方法
CN114151146B (zh) * 2021-10-20 2023-05-05 中国航发四川燃气涡轮研究院 多联带冠涡轮转子叶片气流激振力参数的获取方法
CN114893442B (zh) * 2022-05-09 2023-05-23 北京航空航天大学 一种导叶、压气机及压气机的气动布局设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2003524A2 (en) * 2007-06-15 2008-12-17 United Technologies Corporation Aeroelastic model
CN101599104A (zh) * 2009-07-16 2009-12-09 北京航空航天大学 一种航空涡轮发动机叶片颤振边界的模拟方法
CN105699074A (zh) * 2016-04-01 2016-06-22 山东大学 一种叶片-转子系统流固耦合动态特性试验装置
CN110162813A (zh) * 2018-06-08 2019-08-23 北京机电工程研究所 一种叶轮机械气动阻尼计算方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7082371B2 (en) * 2003-05-29 2006-07-25 Carnegie Mellon University Fundamental mistuning model for determining system properties and predicting vibratory response of bladed disks
US9739167B2 (en) * 2012-07-25 2017-08-22 Siemens Energy, Inc. Method and system for monitoring rotating blade health

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2003524A2 (en) * 2007-06-15 2008-12-17 United Technologies Corporation Aeroelastic model
CN101599104A (zh) * 2009-07-16 2009-12-09 北京航空航天大学 一种航空涡轮发动机叶片颤振边界的模拟方法
CN105699074A (zh) * 2016-04-01 2016-06-22 山东大学 一种叶片-转子系统流固耦合动态特性试验装置
CN110162813A (zh) * 2018-06-08 2019-08-23 北京机电工程研究所 一种叶轮机械气动阻尼计算方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Comparison of Fluid-Structure Coupling Methods for Blade Forced Response Prediction;Florent Payer 等;ASME 2012 Gas Turbine India Conference;第1-8页 *
Evaluation of the Principle of Aerodynamic Superposition in Forced Response Calculations;Stefan Schmitt 等;UNSTEADY AERODYNAMICS, AEROACOUSTICS AND AEROELASTICITY OF TURBOMACHINES;第133-144页 *
带分流叶片叶轮的气动阻尼仿真;王延忠 等;计算机仿真;第36卷(第6期);第77-82页 *
航空发动机叶片气动弹性动力响应的数值方法研究进展;张章 等;力学进展;第42卷(第5期);第572-582页 *
航空发动机叶片气动阻尼的研究;郭雪莲 等;第十一届全国空气弹性学术交流会会议论文集;第470-474页 *

Also Published As

Publication number Publication date
CN111523182A (zh) 2020-08-11

Similar Documents

Publication Publication Date Title
CN111523182B (zh) 一种基于流固耦合的叶片振动响应分析方法
Hall et al. Calculation of three-dimensional unsteady flows in turbomachinery using the linearized harmonic Euler equations
Jennions et al. Three-Dimensional Navier–Stokes Computations of Transonic Fan Flow Using an Explicit Flow Solver and an Implicit κ–ε Solver
Ekici et al. Harmonic balance analysis of blade row interactions in a transonic compressor
Liu et al. Investigation of unsteady impeller-diffuser interaction in a transonic centrifugal compressor stage
Valkov et al. Effect of Upstream Rotor Vortical Disturbances on the Time-Averaged Performance of Axial Compressor Stators: Part 1—Framework of Technical Approach and Wake–Stator Blade Interactions
Chen et al. Comparison of computations using Navier-Stokes equations in rotating and fixed coordinates for flow through turbomachinery
Sun et al. Numerical investigation of surge prediction in a transonic axial compressor with a hybrid BDF/Harmonic Balance Method
Winhart et al. Application of the time transformation method for a detailed analysis of multistage blade row interactions in a shrouded turbine
Yang et al. A novel three-dimensional inverse method for axial compressor blade surface design
Song et al. Aerodynamic and structural multidisciplinary optimization design method of fan rotors based on blade curvature constraints
Stapelfeldt Advanced methods for multi-row forced response and flutter computations
Wang et al. Numerical investigation of the effects of system volume and average mass flow on the surge characteristics of an axial compressor
Vanti et al. An integrated numerical procedure for flutter and forced response assessment of turbomachinery blade-rows
Sondak et al. Simulation of vortex shedding in a turbine stage
Castillon Evaluation of a multiple frequency phase lagged method for unsteady numerical simulations of multistage turbomachinery
CN112632726B (zh) 一种面向叶轮机械叶片气动弹性模拟的流场重构方法
Davis et al. Computational approach for predicting stall inception in multistage axial compressors
Wang et al. Critical endwall blockage attenuation-based automatic optimization of casing treatment design for transonic axial flow compressor
Baturin et al. Application of 3D Navier-Stokes Equations and Mathematical Optimization Techniques to Improve the Efficiency of Seven-Stage Axial Compressor.
Su et al. An efficient coupled method of cascade flutter investigation based on reduced order model
Li et al. Investigations on Unsteady Flow Excitation and Mechanical Performance of Last Turbine Stage Long Blade Using Fluid-Structure Interaction Method
Dalbanjan et al. Sensitivity Study of Stagger Angle on the Aerodynamic Performance of Transonic Axial Flow Compressors
Gao et al. Forced Response Analysis of a Radial Turbine With Different Modelling Methods
Sun et al. Stall Inception Analysis of Transonic Compressors With Chordwise and Axial Sweep

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