NUMECA颤振计算方法介绍.ppt

NUMECA 颤振计算方法介绍 非耦合法：能量法 弱耦合法：模态叠加法、与CSD迭代的流固耦合法 强耦合法 颤振----叶轮机 弹性系统在均匀来流中的自激振动，属于气动弹性稳定性范畴 大量实验研究发现，叶片颤振发作时的频率与低阶的弯曲或扭转固有频率颇为一致 假设叶片颤振发作时的频率与振型同叶片或叶片-轮盘组件的固有频率与振型相一致 除旋转离心刚性效应影响动频外，频率与颤振振型不受周围非定常流场的影响 分析工作大为简化 节径/节圆振动 叶片-轮盘转子组件系统→节径振动 →不同振型 →叶片间相位角（体现叶片与叶片间的相互作用） 同一节径，两重根 问题简化 不同振型→不同叶片间相位角→不同气动阻尼→最小气动阻尼→是否稳定 颤振发生时，所有叶片频率相同，叶片间存在固定的相位差 计算不同叶片间相位角下的气动阻尼，判断稳定性 能量法 叶片以固有频率作初始微幅振动，与周围气体发生能量交换 一个振动周期内系统与外界能量交换的正负来判别叶片颤振是否发生 气动阻尼 = N*单个叶片气动功/(4*系统平均振动动能) = 对数衰减率 略去机械阻尼（远小于负气动阻尼，偏安全） 准则：气动阻尼≥0 结构动力学分析与非定常气动分析独立，非耦合 实现方式 1、给定振动方式--微幅简谐振动（小振幅+频率+相位角，实验常用） 固有频率→振动频率 机理研究、参数化研究 2、给定振型，不耦合每个叶片给定同样大小的振型，按不同时刻给定一个脉冲激励 一次只能算一个叶片间相位角，特定的相位角需要特定的叶片数 An example of a design process: For a given blade motion, check the flow conditions if the flutter is at the onset, damped or amplified using aerodynamic damping parameter 给定确定的位移方式 求解非定常流场 确定非定常累计功和叶片负荷 确定气动阻尼系数 计算流程 刚体运动模块 刚体运动模块下的三种自定义刚体运动方式 Oscillating Body Off-Centered Rotor 标4算例，涡轮624工况结果 刚体运动计算方式，与实验方式吻合 时间步长按一个周期100~200个点选取，弯曲振动 相位延迟边条 气动功折合成气动阻尼系数 趋势一致，能确定最不稳定叶片间相位角 后行波（涡轮）最不稳定 压力面对吸力面影响强于吸力面对压力面 数值不准，特别是不稳定区域 实验相位角范围更宽，绝对值更大 局部稳定性判别标准 与振动方向定义，初始相位角，扭转中心有关，并非唯一 总体一致，耦合方法得到的气动阻偏小，非耦合方法计算结果偏保守 耦合和非耦合计算结果对比 结论 绝对值相同的前后行波，同等扰动，其幅值基本一致，差异主要体现在相位，最不稳定的是后行波 180°相位角，实验和计算结果稳定最好，0°最差 180°相位角下，相邻叶片反向振动，通道面积的变化最显著，是流场中主要的大扰动，相位角误差造成的扰动影响则减弱。 结论 吸力面整体压力幅值大于压力面，压力面压力高，波动小 尾缘压力幅值高于前缘 0-180之间，180相位下压力幅值最大，相邻通道面积差最大，扰动最剧烈 不稳定区域集中在吸力面尾缘，幅值高，是正气动功的主要来源 Numeca中的流固耦合计算方法 Aerodynamic Effects + Structure Deformation MpCCI With FINE?/Turbo or FINE?/Open Modal Approach Prescribed Deformation (One-Way Coupling) FINE?/FSI-Oofelie Fluid-Structure Interaction Aerodynamic Effects + Structure Deformation + Structure Stresses Modal Approach Resolved Deformation (Two-Way Coupling) Sloshing Cold-to-Hot Analysis for turbomachinery Forced Response Analysis Strongly-coupled configurations Flutter Analysis 模态叠加法 对于大多数叶轮机械，气动力对结构域的特征属性的影响可忽略不计 结构域的特征属性由自然频率和振型给定，其作为模态计算的输入 结构域的特征属性可由商用有限元软件或实验测量数据确定 求解刚度方程，刚度方程整合至一个flow solver中，不