Rotor67颤振分析案例.docxVIP

采用TBR-FT方法进行Rotor67叶

片颤振分析

SunilPatil

ANSYSInc

主要内容

?简介(气动弹性,TBR-FT)

?测试案例描述

?Workbench工作流程

?网格细节

?稳态计算结果

?瞬态计算模型细节

?气动阻尼系数预测

?结论

2?2011ANSYS,Inc.June8,2020

气弹问题

SurgePeakEff.ChokeSurgePeakEff.Choke

100%rpm100%rpm

85%rpm85%rpm

70%rpm70%rpm

50%rpm50%rpm

Flutter

CorrectedMassFlowCorrectedMassFlow

只考虑气动因素时压气机特性线

图

综合考虑气弹问题和气动因素

时压气机特性线图

3?2011ANSYS,Inc.June8,2020

气弹问题

流固双向耦合(FSI)

CFDFEM

计算成本高

更切合实际的

分析方法

颤振分析

(气动阻尼系数)

?叶片振动已被提前指定(FEM模态分析)

?基于流体瞬态分析过程进行阻尼系数的计算(基于

CFD)

受迫振动分析

?叶片气动激励由气动计算得到

?计算叶片的振动和应力(基于FEM或ROM).

解耦.

解耦→与ROM耦合

可以预测

叶片是否会发生颤振

可以预测

叶片位移结构失效

4?2011ANSYS,Inc.June8,2020

叶片颤振

气动阻尼系数计算

?叶片颤振(气动阻尼系数计算)

●颤振一般发生于叶片自然振动状态下

●取决于气动和结构的共同作用结果

●模态分析→自然频率和模态振型

●计算时只采用预先指定的模态振型

5?2011ANSYS,Inc.June8,2020

TBR–FourierTransformation方法

M

ftAe??

()

=?

j(mt)m

mM

=?

双通道分析

?=

T

ND

?

ft=ft+?T?ft+?T

()()()BBB

122

ft=ft??T?ft??T

()()()BBB

322

测试案例描述

叶片数22

设计转速16043rpm

设计流量34.07kg/s

雷诺数

(基于叶中弦长)

1.0E6

压比1.68

节径数-8to+8

模态1st弯曲

频率534Hz

叶片最大位移变

量

1.5%,3%弦长

NASARotor67跨音风扇

7?2011ANSYS,Inc.June8,2020

Workbench工作流程

8?2011ANSYS,Inc.June8,2020

叶片网格

9?2011ANSYS,Inc.June8,2020

稳态计算结果

网格Pressure

Ratio

Temperature

RatioEfficiency

粗糙1.681.1792.1

一般1.681.1792.8

精细1.681.1792.6

10?2011ANSYS,Inc.June8,2020

稳态计算结果

3DStreamlines10%叶高速度矢量图

11?2011ANSYS,Inc.June8,2020

稳态计算结果

叶尖马赫数分布50%叶高处马赫数分布云图

12?2011ANSYS,Inc.June8,2020

MeshMotion细节

?相对于初始网格位置发生的网格位移

?Meshmotion模型–Displacement

diffusion(网格刚度与controlvolumes

的二次方成反比

?Meshmotion在叶尖处的设置

?Surfaceofrevolution?Paralleltoboundary

?叶片meshmotion–Periodic

displacement

13?2011ANSYS,Inc.June8,2020

叶片MeshDisplacement设置

14?2011ANSYS,Inc.June8,2020

模型细节

?SST模型+HighspeedWallHeat

transfer模型

?TotalEnergyheattransfer模型+

viscousworkterm

?TBR–FT–叶片颤振

?Timeperiod–1/Frequency

?Numb