超声速平板边界层旁路转捩直接数值模拟

采用具有七阶精度的有限差分格式和基于Euler方程组的特征边界条件,通过直接数值模拟研究了由固定相位周期性吹吸扰动诱导的、Ma∞=2.25、Re∞=2.5×107/m(635 000/inch)空间发展平板边界层的旁路转捩过程,并建立了相应的直接数值模拟数据库。数值结果表明:扰动相位的随机性并不是转捩发生的必要条件,固定相位周期性吹吸扰动即可诱发边界层的转捩;当采用固定相位时,随着吹吸频率的增加,转捩位置向下游移动;当扰动频率足够大时,周期性吹吸无法诱导出边界层的转捩。随后,在所得数据库中提取了速度及热力学变量的平均值,与已有的数据库和理论分析吻合良好;速度脉动均方根与公认的不可压缩平板边界层实验基本吻合;热力学变量脉动均方根和雷诺应力与已有数据库的统计结果及实验数据具有相同的分布特征。对比发现,不同转捩的诱导因素对平均转捩位置和热力学变量均方根、雷诺应力等高阶统计量的影响较大。     

硅微陀螺平板结构的压膜阻尼研究

压膜阻尼对硅微陀螺的精度具有重要影响。平板结构是硅微陀螺中一种典型结构,本文通过对雷诺方程的推导,得出平板在正弦规律小振幅刚性运动时的阻尼力解析式,并给出低频时阻尼系数的近似算法。利用Matlab和CoventorWare软件对运动平板的阻尼进行仿真计算,得出运动平板的阻尼系数和刚度系数,并分析平板参数对阻尼系数及刚度系数的影响。建立不同类型的带孔平板模型,通过比较可以看出,当阻尼孔按一定方式排列时,可以降低结构的阻尼。     

基于轮廓法测试铝合金拉弯成形残余应力研究

以2026-T3511铝合金T形结构件的拉弯成形件为研究对象,首先采用线切割将T形结构件分割成两块具有T形截面的试样,然后基于轮廓法测试应力原理应用激光位移传感器分别测试两个分割截面的轮廓变化,进而将两个截面的测量值拟合平均后与有限元模型结合以计算该试样的初始残余应力,从而获得该T形结构件的横截面法向应力并分析其分布特征。研究结果表明,基于轮廓法测得的应力分布与有限元仿真结果分布规律相近:峰值压应力为60MPa,峰值拉应力为40MPa;轮廓法在航空结构件应力测试中能够准确地获得其他测试方法无法获得的横截面完整的应力分布。     

连续尺度变换方程在稳态及非稳态流动模拟中的应用

为了提高涡粘性假设的湍流模型对于非稳态流动的求解精度,同时兼顾其对于稳态流动的求解性能,将雷诺应力项与连续变换方程(CSSE)结合而形成新的应力项,使其根据流场尺度、网格尺度及Kolmogorov尺度来自动调节当地的应力雷诺应力模化水平,避免网格因素在流场模拟中产生不利影响,改正了混合RANS/LES方法的速度型偏离对数率问题;同时,该方法并未引入显式亚格子模型(SGS),因此回避了亚格子系数确定对于流场模拟精度产生的影响,改善了湍流模型对于流动不稳定性的辨识精度。在湍流平板算例中,CSSE方法计算的边界层速度型精度与雷诺平均方法(RANS)相当,而对于圆柱尾迹的模拟则证明了CSSE方法具有混合RANS/LES方法的优点,即能够准确模拟流动的不稳定性特征。     

7085铝合金残余应力及加工变形的数值仿真与试验

毛坯残余应力的测量以及零件的加工变形分析是航空整体结构件数控加工工艺研究的难题。应用MSC.Marc软件对某型飞机主起支撑接头的7085铝合金锻件毛坯的淬火-压缩工艺进行了数值仿真,获得了毛坯残余应力分布趋势;应用X射线衍射法与钻孔法测量了毛坯表层残余应力分布,提出了基于系数修正法的毛坯残余应力分析方法,通过最小二乘法（LSM）拟合仿真值与试验值以获得应力修正系数,并对毛坯残余应力分布的仿真结果进行修正,获得了实际毛坯的残余应力分布;由此进行了主起支撑接头缩比零件加工变形仿真与验证试验。研究结果表明：7085铝合金块状毛坯残余应力数值仿真结果与测试结果的趋势吻合,为应力修正提供了物理基础;相对于未修正的应力分布,基于修正后的应力分布仿真获得的变形结果精度提高了50%;7085-T7452铝合金的毛坯残余应力是导致主起支撑接头加工变形的主要因素。     

基于CST方法的高空低雷诺数吸附式叶型耦合优化设计

研究了一种基于类别形状函数变换（CST）方法的吸附式叶型优化设计方法,该方法可以在高空低雷诺数条件下对叶型和抽吸方案耦合优化.结果表明:优化之后在20km高空低雷诺数条件下总压损失降低了65%,静压升提高了0.02,气动性能得到较大提升.而且由于优化过程中罚函数的引入使得优化后吸附式叶型在地面条件下性能也有所提高.对于高空低雷诺数条件下吸附式叶型在抽吸位置之前适当的增加叶型负荷,再通过抽吸来控制附面层,效果最优.并且最佳抽吸位置位于层流分离泡作用区域内.在层流分离泡作用区域内抽吸可以完全消除层流分离泡对叶型性能的影响,并且可以较好控制附面层位移厚度和动量厚度的增加,有效地减小附面层内的动量损失.      

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。     

涡轮导向叶片综合冷却效率实验研究

为了研究流动参数对涡轮导向叶片综合冷却效率的影响,采用红外热像仪对叶片表面的温度分布进行了测量,得到了叶片的综合冷却效率随流量比、温比、主流进口雷诺数和湍流度的变化规律。实验过程中,次流与主流的流量比分别为0.15,0.18,0.20,0.22和0.24;主次流温比分别为1.4,1.7,1.93和2.2;主流进口雷诺数分别为1.0×10~5,1.1×10~5,1.2×10~5,1.3×10~5和1.4×10~5;主流进口湍流度分别为0.506%,8.156%,14.92%。结果表明,综合冷却效率在前缘处最低,沿流向逐渐升高;增大流量比会显著提高叶片的综合冷却效率,在温比为1.93时,流量比由0.15增大至0.24,综合冷却效率平均增加29.3%;温比和主流进口湍流度的增大均不利于综合冷却效率的提升,流量比为0.20时,温比由1.4增大至2.2,综合冷却效率平均下降46.5%,湍流度由0.506%增大至14.92%,综合冷却效率平均降低15.5%;主流进口雷诺数对叶片综合冷却效率的影响很小。     

叶片式预旋喷嘴出口流场实验研究

叶片式预旋喷嘴具有尺寸小,落后角大的特点。为了详细研究小尺寸预旋喷嘴的预旋性能,采用五孔探针对叶片式预旋喷嘴的出口流场进行了实验研究。测量了Ma=0.2,0.3时喷嘴出口的压力分布、速度分布和出口气流角度分布,实验获得了喷嘴的落后角和预旋效率,并进行了与实验工况相同的数值计算。通过实验获得的总压云图以及速度云图,可以发现叶片式预旋喷嘴的端壁二次流损失、尾迹损失严重,有明显的边界层分离现象。Ma=0.2时,喷嘴Re数为5.76×104,落后角2.84°,实验测得的预旋效率为0.73;Ma=0.3时,喷嘴Re数为1.06×105,预旋效率提高至0.77。实验模型端壁的影响使预旋效率实验结果偏低6.5%左右。数值结果与实验测得各参数符合较好:数值结果与测得的喷嘴出口截面平均总压、静压偏差在1%以内;出气速度、周向速度以及出气角度与实验结果偏差在4%以内。数值计算表明,叶片式预旋喷嘴的预旋效率基本不受压比影响,随Re数增大先增大后基本不变,最后基本稳定在0.85。     

大型水陆两栖飞机侧风起降的滑流影响分析

大型水陆两栖飞机AG600的动力装置为安装在机翼上的4台同向旋转涡轮螺旋桨发动机,针对1：15缩比模型带动力风洞试验显示的螺旋桨滑流对侧风起降状态的偏航力矩不稳定影响,对全机带动力风洞试验模型进行了大规模并行非定常数值计算,再现了风洞试验现象,通过流动机理分析明确其产生原因主要是左侧滑时右外翼分离和垂尾背鳍涡破裂,这些原因和数值模拟的准确性也为后期的风洞试验所证实。考虑到模型风洞试验中尺度限制造成的低雷诺数和高螺旋桨转速,为保证飞行安全,继续采用该非定常方法对全尺寸飞机真实侧风起降状态进行了详细数值分析和偏航稳定性评估。研究结果显示,在飞行雷诺数和螺旋桨转速下,相同侧风范围内风洞试验显示的流动不稳定因素基本消失,偏航稳定性允许的侧风范围明显增加。本研究实现了四发螺旋桨飞机起降状态横向气动特性的滑流影响非定常数值分析,建立了基于计算流体力学的风洞与飞行雷诺数效应的相互关系,进行了偏航稳定性的虚拟试飞评估,研究成果也为AG600飞机的首飞和飞行试验所验证。