基于CFD/CSD耦合的叶轮机叶片失速颤振计算

为了研究叶轮机叶片的失速颤振特性,发展了一种计算流体力学与计算结构力学(CFD/CSD)时域耦合方法。该方法通过每一物理时刻CFD和CSD的循环迭代实现了耦合计算。在CFD分析中,采用鲁棒性较好的空间离散格式AUSM+-UP,并基于延迟脱体涡模型(DDES)模拟了带分离流动。在结构分析中,通过模态法构建了旋转叶片动力学方程并运用杂交多步方法进行求解。以孤立转子Rotor37为例,计算了不同工况下流场总体与细节参数,与实验结果的对比验证了CFD算法的精度。对某转子叶片进行了颤振特性研究,计算所得的广义位移时间响应曲线表明该叶片在近失速工况下会发生失速颤振,其表现形式为一阶弯曲模态发散且各阶模态之间不耦合。分析表明,流场不稳定和非定常效应是引起失速颤振的关键因素,同时折合频率的降低也会导致原本气动弹性稳定的叶片发生失速颤振。     

考虑空腔的高超声速多流动区域同步数值模拟

先进的高超声速飞行器具有薄壁空腔结构,在飞行过程中受热会产生空腔内气体流动现象,从而影响流场和结构的温度分布。采用数值方法准确模拟高超声速流场、结构温度场和空腔内流动对热结构分析是很有必要的。以研究空腔流动对结构温度分布影响为目的,发展了一种适用于多流动区域流场/结构温度场耦合问题的同步计算方法,并以高超声速带空腔结构物体为例,数值研究了其外部气动热/结构热传导引起的空腔热对流问题。以已发展的高超声速外流场/结构温度场同步计算方法为基础,为了进一步考虑空腔内低速流场,采用了预处理矩阵方法。在流场与结构温度场的交界面两侧分别引入虚拟单元,从而高效地实现相邻场之间物理信息交换。首先通过标准算例验证了方法在求解单独气动热/结构热传导问题以及空腔自然对流问题中的准确性。进而对封闭和带有开孔的两种高超声速运动圆环分别进行多流动区域同步数值模拟。计算结果表明,由于结构温度不均匀引起的空腔内热对流反之也会对结构温度场分布产生轻微的影响。在空腔内气体流动的影响下,封闭圆环的前缘温度在35s内最多下降0.8%左右。对于带开孔空腔的圆环,其孔壁周边温度在0.5s内能够超过外流前缘驻点温度。     

高超声速全动翼面全时域耦合分析方法及应用

在流场-结构温度场同步计算方法的基础上,建立了多物理场全时域耦合分析方法,将方法应用于沿轨道运动的高超声速全动翼面热气动弹性稳定性分析。采用基于有限体积模型的CFD同步计算方法求解高超声速流场和结构温度场,建立映射关系实现结构有限元模型气动载荷加载和温度场赋值。采用移动坐标系和动网格相结合的方式描述变速度飞行和翼面偏转过程。通过坐标系变换将翼面偏转过程和振动过程的网格变形量叠加,考虑翼面振动和偏转的耦合非定常效应。针对沿轨道运动的高超声速飞行器,建立了同步计算方法与全时域耦合分析方法相结合的热气动弹性稳定性分析流程。研究发现,与同步计算方法相比,全时域耦合分析方法能够模拟结构振动对流场和结构温度场的影响,计算得到的监测点热流密度波动幅值占热流峰值的10%左右,而温度变化并不明显,相比于刚体模型,监测点温度只下降了0.3%左右。全时域耦合热气动弹性分析方法得到的颤振临界点在4-5号状态点之间,颤振形式为铰链扭转模态与一阶弯曲模态的耦合颤振,与"冻结"模态的热颤振方法结果一致。