梯形翼风洞试验模型数值模拟技术

基于雷诺平均Navier-Stokes (RANS)方程和结构网格技术,采用二阶空间离散精度的MUSCL格式,并结合k-ω剪切应力输运(SST)两方程湍流模型和γ-Re_θ转捩模型,研究了梯形翼风洞试验模型中前缘缝翼、后缘襟翼连接装置对气动特性的影响。简要介绍了本文采用的计算方法;介绍了梯形翼的风洞试验模型及风洞试验结果;在网格收敛性研究的基础上,采用"全湍流"方式和转捩模型研究了梯形翼试验模型连接装置对气动特性的影响。通过与不带连接装置的计算结果的对比,采用"全湍流"模拟方式,计算模型中考虑试验模型的连接装置引起升力系数下降、阻力系数下降、低头力矩减小以及失速迎角提前;通过与试验数据的对比,进一步考虑转捩影响可以提高梯形翼风洞试验模型气动特性的计算结果与试验结果的吻合程度,梯形翼风洞试验模型失速迎角附近的气动特性数值模拟技术还需要进一步的研究。     

第一届航空CFD可信度研讨会总结

总结了第一届航空可信度研讨会(AeCW-1)的数值模拟结果。AeCW-1会议由中国空气动力学会计算空气动力学专委会主办,并在第一届空气动力学大会上做了专题研讨。本次研讨会吸引了14家单位参与,共提交了26组计算结果。活动的主要目的是针对典型运输机巡航构型,评估国内CFD当前技术状态,探索CFD进一步发展的研究方向,促进CFD验证和确认工作的开展,为大型客机的研制提供技术支撑。AeCW-1选择的研究模型为中国空气动力研究与发展中心(CARDC)设计的单通道运输机模型(CHN-T1),CHN-T1包括机身/机翼/平尾/立尾等部件,分别在CRADC的2.4m×2.4m风洞(FL-26)、欧洲的2.4m×2.0m风洞中(European Transonic Wind Tunnel,ETW)开展了风洞试验。AeCW-1组委会提供了系列结构网格和非结构网格,研究工况主要包括网格收敛性研究、抖振特性研究、雷诺数影响研究,风洞模型支撑装置和静气动弹性变形对数值模拟结果的影响是重要的研究内容。本文对提交组委会的数值模拟结果进行统计分析,并与FL-26及部分ETW的风洞试验结果进行了对比,对未来的研究工作提出了意见和建议。     

飞行仿真气动力数据机器学习建模方法

基于机器学习思想,提出了一种大空域、宽速域的气动力建模方法。该方法利用飞行仿真弹道数据辨识的气动力数据,采用人工神经网络技术,实现了对高度、速度、姿态和舵偏角等多维度强非线性特性的全弹道气动力数据的高精度逼近。首先,分析了神经网络层数、隐含层神经元个数等对建模误差的影响,通过对典型弹道气动数据的神经网络建模计算,确定了较合适的神经网络层数和较优的隐层神经元个数。进而,利用飞行仿真的弹道数据辨识出沿弹道的气动力,采用神经网络建立了包含多个弹道融合的气动力模型,输出量分别为三轴气动力系数和力矩系数。最后通过气动模型输出量与原样本数据的对比,以及4条未参与训练弹道气动数据的预测,验证了该气动力建模方法具有较高的精度。建模结果表明:采用神经网络方法建立的飞行器气动力模型,对拟合多源耦合输入全弹道非线性气动力是可行的和有效的,在样本覆盖的高度、速度、姿态和控制舵偏角范围内,气动力拟合能力较强,并具有一定的外推性。该项研究可以为基于飞行试验数据的气动建模提供新的方法,并且能为飞行器气动力数据挖掘、飞行仿真和总体性能分析提供参考。     

基于成分的燃烧室出口燃气温度计算方法

针对燃烧室不同温度下燃气组成的不同,提出建立不同燃气模型以推算燃气温度的两种计算方法。基于物质守恒、化学平衡和能量守恒原理建立温度计算方程,根据能量守恒方程的不同表达形式给出平均比定压热容法和焓值法,同时给出两种求解非线性方程组的数学解法:三变量迭代法和牛顿法,并采用VB(Visual Basic)编制燃气分析测试程序实现温度算法。列举计算实例比较两种算法的区别,结果显示:随着油气比的增加,两种计算方法得到的燃气温度相差越大,考虑热离解的高温算法结果更小,原因在于燃气中离解组分的摩尔分数随着温度升高而增加,燃烧实际放热量减小,导致燃气温度减小,所以当燃气温度较高时建议采用高温算法更符合实际燃烧情况。      

基于UKF的共轴式无人直升机模型辨识

建立了共轴式无人直升机系统非线性模型,并针对其非线性强,不同飞行模态下气动参数差异等问题,将无迹卡尔曼滤波(UKF)引入共轴式直升机系统非线性模型辨识,不但避免了直升机线性模型仅仅适用于悬停模态的局限性,同时为直升机系统在线自适应控制提供了基础条件,使得共轴式无人直升机自主全包线飞行成为可能.以北京航空航天大学FH-1共轴式无人直升机为例进行了仿真辨识实验.实验结果表明基于该方法的共轴式直升机在线非线性模型辨识不依赖于参数初值的选取,模型参数能在10s内收敛,各状态量辨识精度达到80%以上,明显高于传统的预报误差法(PEM),具有一定的实用性.      

基于PD-MSM-Lyot 滤波器的多波长掺铒光纤激光器

提出了一种基于偏振依赖多模- 单模- 多模光纤Lyot (PD-MSM-Lyot) 滤波器的 多波长掺铒光纤激光器。PD-MSM-Lyot 光纤滤波器由一个PD-MSM 滤波器和一段保偏光 纤构成,其相邻滤波通带具有相同波长间隔。充分利用偏振依赖MSM 光纤滤波器的起 偏作用并结合腔内色散位移光纤中的非线性偏振旋转效应, 构成了强度相关损耗机 制, 从而抑制掺铒光纤中的增益竞争作用, 获得了室温下19 个波长间隔为0.2 nm 的密 集多波长激光输出。通过对比可得, 输出激光光谱随泵浦功率增大而变的更平坦, 且 波长数更多。     

适用于倾转旋翼机气动分析的新型嵌套网格方法

为克服围绕复杂倾转旋翼机生成高质量一体化结构网格的困难,将分块网格策略与嵌套网格方法融合,提出并建立了一套适用于倾转旋翼机气动干扰分析的新型嵌套网格生成方法。该方法中,将整体网格按照倾转旋翼机结构特性进行了分块生成,在保证插值精度的基础上降低了网格规模。为高效连接不同部件网格,首先,引入Inverse Map方法来辅助部件网格定位,采用新颖的多向投影方法保证了部件间洞边界的连续性;其次,基于独立挖洞嵌套策略完成机身与背景网格的组合,并使用了相应的辅助计算坐标进行背景边界网格标识。在贴体网格区域采用基于S A（Spalart Allmaras）湍流模型可压缩RANS(雷诺平均Navier Stokes）方程求解流场,背景网格区域采用欧拉方程,并运用了SPMD(single program, multiple data)模式的并行加速技术,建立了适合倾转旋翼机气动干扰特性分析的高效混合CFD方法。采用Robin直升机机身作为数值算例,验证了CFD方法的有效性。在此基础上着重对倾转旋翼机在悬停、过渡和前飞巡航模式下的干扰流场进行了数值分析研究,分别得出了悬停状态下倾转旋翼机典型的“喷泉效应”干扰现象、过渡状态下的非定常气动干扰特性以及巡航状态下的部件气动变化特性,计算结果表明建立的新型嵌套网格方法能够较好地表征倾转旋翼机外形特性,并能够有效地用于倾转旋翼机的气动特性分析,加速比能超过5.0。      

高重合度齿轮传动系统的振动参数稳定性

研究了高重合度齿轮转子-滚动轴承传动系统的振动强度稳定性以及振动轨道稳定性随系统参数的变化规律。通过振动强度参数稳定域的计算,量化了高重合度齿轮副对齿轮转子 滚动轴承传动系统稳定性提升的贡献大小。结果发现在转子质量偏心较小的条件下,高重合度齿轮系统的振动强度稳定性要比普通重合度齿轮系统提升很多;但在转子质量偏心较大时,高重合度齿轮副的采用对传动系统参数稳定域的扩张并没有显著贡献。通过振动轨道稳定性全局分岔图的计算,明确了滚动轴承游隙、转子质量偏心等参数对高重合度齿轮转子-滚动轴承系统振动轨道稳定性的影响规律,获得了系统各种稳定周期轨道及非周期轨道与对应参数区间的映射关系。      

热声载荷下薄壁结构振动响应试验验证与疲劳分析

由于热声环境下金属薄壁结构表现出复杂的大挠度强非线性振动响应特性,影响结构的疲劳性能与寿命,结合有限元法与降阶模态法对四边固支高温合金矩形薄壁结构的热声响应进行计算。结果研究发现:屈曲后结构出现跳变运动且应力循环呈三角状分布,热声载荷的相对强弱决定了跳变形式。采用改进雨流计数法、Morrow平均应力模型、Miner线性损伤累积理论计算热声疲劳寿命,屈曲前到临界屈曲时应力循环损伤量级显著增大,由10-5增大到10-4,寿命随温度增加呈先减小后增加趋势。开展薄壁结构热声试验,并将仿真计算结果与试验结果进行对比,结果表明结构的模态频率偏差不超过1Hz,动态应变响应结果的量值相当,验证了薄壁结构热声响应计算方法与模型的有效性。      

基于灰色补偿反步法的机动指令跟踪控制

基于三自由度飞行动力学模型和运动学模型,提出了一种带有灰色补偿的反步(Back-stepping)控制方法,以实现某些预定的战术机动动作。首先,将飞行动力学模型和运动学模型分解成3个轴方向的相对简单的子系统;然后,针对3个子系统采用反步法分别设计出模型中没有干扰时的控制律,并采用GM(0,N)预测模型对不确定部分的模型参数进行辨识,进一步根据估计出的参数设计灰色补偿控制律;最后,根据3个子系统设计出的控制律解算出同时满足3个子系统渐进稳定的控制输入,这样使得整个系统渐进稳定。仿真结果表明,GM(0,N)预测模型能精确预测不确定模型参数,该算法能够高效实现预定的战术机动动作。