半潜入喷管收敛段碳/酚醛材料烧蚀特性

采用流固热耦合数值方法,以及Abaqus的ALE（arbitrary Lagrangian-Eulerian,任意拉格朗日欧拉方法）自适应网格技术,对位于半潜入喷管收敛段的碳/酚醛的烧蚀现象进行了预估,与试验结果进行了对比分析。结果表明:在收敛段,Al/Al₂O₃液滴或颗粒对材料的传热烧蚀起到了关键作用,当气流与材料表面平行时,液滴或颗粒的影响很小;从碳/酚醛热分解角度出发,基于完全碳化即被剥蚀的假设,基本能够预估碳/酚醛材料的烧蚀特征,烧蚀速率大约为0.3 mm/s;后效碳化对材料碳化过程的影响明显,发动机工作期间,分解层的厚度大约为2.0 mm,考虑后效作用,分解层厚度可能会增加至4.0 mm左右;与喉衬接触的材料区域,喉衬的传热会明显加剧材料的碳化过程。      

适用于涡激振荡问题研究的并行高精度方法

提出了一种用于研究大振幅及物体间小间距等极端情况下涡激振荡（VIV）问题的高精度方法。该方法针对三角形/四边形非结构混合网格,采用高精度谱差分（SD）格式对Navier-Stokes方程进行空间离散。通过引入非均匀滑移网格方法,将计算域分割成多个互相不重合的子区域,从而实现了子区域网格的独立变形,子区域之间通过粘接元进行信息传递。采用全耦合方法精确求解VIV问题中流体和固体间的相互作用。通过研究并行策略以及使用消息传递接口,实现了求解器的并行计算能力。数值模拟表明：对于无黏和有黏流动问题,该方法均能保持SD方法的高精度特性;对于动网格下的定常均匀来流,求解器能够做到自由流保持;单个圆柱VIV问题仿真与现有文献结果符合较好,验证了方法的可靠性;对于流场变形情况比较复杂的涡激振荡问题,求解器可以有效实现网格变形,并具有理想的并行效率。     

基于致动线方法的5 MW海上风力机气动弹性分析

随着风力机功率的不断增大以及新型复合材料的应用,叶片的柔性和几何非线性变形成为风力机设计中不可忽略的因素,结构和气动弹性的分析也随之变得更加复杂。然而,风能行业中传统的叶片分析方法无法准确预测现代复杂叶片的气动弹性特性,从而导致风力机性能预测出现较大误差。本文基于柔性多体动力学,建立了一种新型双向流固耦合模型,在结合致动线方法和大涡模拟的基础上,考虑了结构和气动弹性对风力机性能的影响,可用于动态结构载荷预测及流固耦合分析。对5 MW基准风力机进行建模,验证了计算模型的准确性,并讨论了叶片的瞬时结构响应,分析了叶片变形对风力机功率、尾迹的影响。研究结果表明,叶片的柔性在风力机气动弹性设计中不可忽略,同时本文模型可以准确捕捉风力机的尾迹结构（包括叶尖涡和叶根涡）,更适用于现代兆瓦级复合材料风力机气动弹性和尾迹分析。     

水下环境爆炸对方形水池冲击载荷数值模拟研究

水下环境爆炸问题属于较为复杂、较难进行试验的科学研究,理论公式对水下环境爆炸相关问题分析能力有限,随着数值仿真技术准确性和精确性的提高,利用数值仿真技术对水下环境爆炸问题进行分析显得尤为重要。开展了爆炸冲击对水池壁面作用力数值仿真,采用流固耦合算法和remap技术,将一维楔形爆炸模型映射到三维球体模型中,以较小的计算规模模拟3kg TNT药包在较大水池（9000mm×5000mm×5000mm）内爆炸冲击过程,揭示了爆炸冲击波在水池中发生多次反射机制,得出了水池壁面受力规律,仿真计算与分析验证了本文提出的模型方法的有效性。     

飞机液压系统出油管路振动抑制方法研究

随着飞机液压系统的高压、高功率化发展,液压管路的振动问题开始凸显出来,降低液压管路的振动,对于提高飞机航行安全性具有重要意义。本文针对航空液压管路进行振动特性研究,获得管路振动固有频率随内部压力和流速的变化规律,进一步根据液压管路动力学参数可变的特点,设计刚度可调的半主动吸振装置,建立适合复杂系统动力学响应分析的多软件联合仿真方法,验证吸振装置和半主动控制算法的振动抑制效果。仿真结果表明,半主动振动控制方法对液压管路系统的振动能够起到有效的抑制作用,在吸振器有效频率范围内管路振动衰减能够达45dB以上。     

基于流固耦合的斜激波冲击作用下曲壁板气动弹性分析

采用自主发展的双向流固耦合求解器,研究了斜激波冲击作用下曲壁板的气动弹性响应特性。曲壁板的几何非线性大变形运动方程采用有限差分法求解,流体控制方程基于有限体积法求解,双向流固耦合采用交错迭代算法。计算结果表明:当动压小于临界颤振动压时,曲壁板表现出静平衡状态,且随着动压的增大,壁板变形的非对称性越明显。当动压大于临界颤振动压时,壁板振动位移先增大后减小,最终达到稳定颤振状态,且该极限环颤振并不关于初始位置正负对称的。同时,随着动压的增大,壁板颤振的正向峰谷值、负向峰谷值和振幅均逐渐增大,颤振频率则逐渐减小。壁板振动响应规律并不随着壁板弯曲高度的改变而单调递增或递减,较小的弯曲高度可以降低壁板颤振临界动压值,但是当弯曲高度进一步增大后,由于气动非线性特性增强,准周期无规则运动状态被激发了出来,临界颤振动压迅速升高。     

非线性静气动弹性分析方法在阻力方向舵上的应用

对于无尾飞翼布局飞机,阻力方向舵是一种十分有效的航向操纵装置,在开裂状态下利用其产生的阻力形成的偏航力矩进行偏航控制.在这一过程中,阻力方向舵的气动分布具有较强的非线性特性,由此导致的舵面弹性变形会对其操纵效率产生较大影响.采用基于CFD/CSD的流固耦合方法,完成了非线性静气动弹性设计建模及仿真计算方法研究,开展了阻...     

一种屏蔽式总温探头设计及其测量误差估计

在忽略辐射误差的情况下，使用流固耦合传热数值模拟的方法预测了一种总温探头在不同工况下的速度误差和导热误差。结果表明：在马赫数为0.2～0.6的范围内，导热误差都保持在较小值；总误差的90%以上由速度误差贡献；最大测量误差为1.13 K，比结构A、结构B、结构C以及结构D分别低了331.6%、119.4%、61.6%以及59.5%；速度误差和导热误差存在互相影响的关系，存在最优解可使总误差最小；适用于等温来流和马赫数大于0.5的高速非等温来流的总温测量。最后探讨了热电偶节点位置对测量精度的影响，节点在从设计点远离支持体的过程中，测量误差遵循着先增大后减小的变化规律。     

跨音速风扇全环叶片颤振特性的流固耦合分析

发展了求解叶片颤振问题的流固耦合计算方法和全环叶片振动的气动弹性模型,在每一时间步同步求解流体运动方程和叶片振动方程并交换边界信息;流体域求解了非定常雷诺平均N-S方程,得到每一步由于叶片变形而引起的流场变化;叶片变形则由积分叶片表面受到的气动力并求解结构动力学方程得到.颤振分析是在全环叶片模型上进行的,并解除了预先设定叶片间相位角的限制.此方法的显著特征是在一次气动弹性计算过程中,可同时分析叶片多个固有模态、多个节径下的气动弹性稳定性,大大提高了使用时域法进行叶片排气弹分析的计算效率.考察了NASA rotor 67风扇全环模型在堵塞点、最高效率点和近喘点3个气动工况下,节径变化对叶片气动弹性稳定性的影响,给出了不同模态下气弹最不稳定状态对应的叶片振动节径形式.结果表明,振动形式对于叶片气动弹性稳定性的影响很大.