适用于多迎角结冰试验的混合翼型设计

为突破结冰风洞对翼型模型尺寸的限制,提出了一种新的混合翼型设计方法,可使用一套混合翼型模拟原始翼型在不同迎角下的结冰试验,弥补了以往混合翼型只能在单个设计迎角下使用的缺陷。方法采用多段翼的形式设计混合翼型,以多目标迎角等结冰试验条件作为设计输入,优化设计主翼外形和襟翼的位置、偏转角度,利用襟翼位置和偏转角度的变化实现混合翼型和原始翼型前缘表面的压力系数在不同迎角下能够保持一致,继而保证前缘结冰外形的一致性。设计的混合翼型较原始翼型弦长减小接近40%,在冰风洞中对混合翼型和原始翼型在目标结冰条件下进行试验,对比结果显示,在选取的多个目标迎角下,混合翼型和原始翼型二维截面的结冰外形基本一致,提出的混合翼型设计方法能够有效的模拟原始翼型在不同迎角下的结冰外形。     

基于自适应差分进化算法的变循环发动机模型求解方法研究

为了降低传统迭代算法在求解变循环发动机非线性模型时对初值的依赖性,将模型的求解问题转换为求最小值的优化问题,引入差分进化算法进行模型的求解,并提出一种自适应差分进化算法(ADE)。ADE借助轮盘赌选择法,利用种群的进化经验可以自适应的选择最适合当前种群的差分策略与算法控制参数。针对变循环发动机四个典型工作点的模型求解问题,研究了标准差分进化算法(SDE)的控制参数对其性能的影响,获取了SDE在求解四个典型工作点时的最优控制参数组合,对比分析了ADE与SDE的性能差异,最后研究了种群规模对ADE性能的影响。结果表明：SDE在求解发动机模型时具有较好的鲁棒性,在求解不同工作点时算法的最优控制参数并不完全相同；相比于使用最优控制参数的SDE,ADE可以在不影响算法鲁棒性的情况下提升效率50%以上；减少ADE的种群规模会在提升算法效率的同时破坏鲁棒性。     

共轴刚性旋翼气动干扰数值计算方法

建立了一个适用于共轴刚性旋翼气动特性分析的数值模拟方法。该方法采用任意拉格朗日欧拉方法(Arbitrary Lagrange Euler,ALE)描述的可压缩Navier-Stokes(N-S)方程求解流场,采用低数值耗散的Roe格式进行空间离散;使用多重嵌套网格方法以模拟双旋翼的运动。针对共轴刚性旋翼配平,引入"差量修正"策略解决了传统配平中雅克比矩阵计算复杂的问题。首先,对Harrington-2共轴双旋翼的悬停气动性能进行了计算,然后,对某2 m直径共轴双旋翼的悬停及前飞状态进行了计算,并与试验值进行了对比。结果表明:在典型状态下拉力系数的计算结果与试验值误差在3%以内,扭矩系数的计算结果与试验值误差基本在5%以内;所采用的数值计算方法对旋翼涡尾迹特征具有较高的捕捉精度,可以有效模拟共轴刚性旋翼悬停和小速度前飞下的复杂流场及其细节特征。     

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

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

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

基于Deform软件二次开发和BP神经网络的TA15多向锻造微观组织预报

通过热压缩实验构建了TA15合金的热变形真应力-应变曲线，以此为基础分别建立了合金双相区及单相区温度区间的热变形本构方程；基于热压缩试样动态再结晶的统计数据建立了TA15合金的动态再结晶模型。借助Deform提供的二次开发功能实现相关数学模型的程序化，制定正交实验方案，实现了TA15合金多向锻造变形的微观组织仿真。通过正交实验分析得出各项因子的影响对象及强弱差异，提出了双相区及单相区温度区间内的多向锻造最佳因子组合。建立了TA15合金多向锻造变形微观组织的BP （Back Propagation）神经网络预测模型，将预测结果与有限元仿真结果进行比较，结果表明两种方法的预测结果基本一致，但神经网络具备有限元仿真难以实现的良好细节预测能力，能更为细致地实现对微观组织分布状态的划分。     

高雷诺数槽道湍流的壁面模化大涡模拟研究

选取基准壁湍流的槽道流动,研究了多种模型的壁面模化大涡模拟.模型包括经典的大涡模拟、Spalart-Allmaras、分离涡模拟和一种动态混合模型.基于摩擦速度的雷诺数范围为395~12 000,采用3组粗糙网格,流向和展向维数分别同取37,49和65,法向维数保证y+(1)~1.主要研究平均速度、雷诺切应力分布、详细分析了各模型的特性差异并展示了相应的湍流结构.研究表明:在高雷诺数粗糙网格下,大涡模拟失去求解精度,分离涡模拟出现对数律不匹配,动态混合模型的计算接近直接数值模拟,其对数率区可解应力约占雷诺切应力的93%,边界层外层可解应力约占99%.这说明合适的混合模型可以在经济成本下保证计算精度,具有解决实际问题的潜力.     

旋转机翼悬停气动特性研究

鸭式旋转机翼(CRW)是一种先进的高速直升机方案,旋翼同时也是机翼,是其最关键气动部件之一。采用结构分区拼接网格技术进行空间离散,分区建立参考系,通过求解多重参考系下的N-S方程来计算旋翼流场,首先以传统的Caradonna-Tung实验旋翼的亚、跨声速悬停流场分析为例验证该方法的可靠性,进而采用该方法对旋转机翼悬停流场进行了数值计算,旋翼拉力计算值和地面实验值吻合较好,结果分析表明旋转机翼的悬停流场有着不同于传统旋翼的流场特性。     

战损飞机动力学特性分析及杀伤研究

现代战争中,防空武器性能不断提升,对飞机造成的威胁越来越大,确定飞机受到打击时的杀伤状态及杀伤等级成为飞机设计时需要考虑的重要问题。飞机被防空武器击中时,可能发生机翼脱离损伤,其质量、面积产生明显变化,将导致飞机的重心、气动特性的变化,严重危害飞机的稳定性。针对上述情况,以某型作战飞机为例,建立损伤情况下的飞机模型,综合分析机翼脱离损伤对飞机重量、重心和机翼面积带来的影响;建立损伤飞机的气动模型,计算损伤飞机气动特性的变化情况。基于这些变化,分析损伤对飞机动力学特性的影响,建立损伤状态下的动力学方程,计算损伤飞机的失控时间,依此确定飞机的杀伤等级。本文所做研究可用于在飞机设计过程中分析其杀伤状态及杀伤等级,为飞机设计时飞机杀伤的分析评估提供一种有效的解法。     

A coupling VWM/CFD/CSD method for rotor airload prediction

A coupling fluid-structure method with a combination of viscous wake model(VWM),computational fluid dynamics(CFD) and comprehensive structural dynamics(CSD) modules is developed in this paper for rotor unsteady airload prediction. The hybrid VWM/CFD solver is employed to model the nonlinear aerodynamic phenomena and complicated rotor wake dynamics;the moderate deflection beam theory is implemented to predict the blade structural deformation; the loose coupling strategy based on the ‘delt method' is used to couple the fluid and structure solvers.Several cases of Helishape 7A rotor are performed first to investigate the effect of elastic deformation on airloads. Then, two challenging forward flight conditions of UH-60 A helicopter rotor are investigated, and the simulated results of wake geometry, chordwise pressure distribution and sectional normal force show excellent agreement with available test data; a comparison with traditional CFD/CSD method is also presented to illustrate the efficiency of the developed method.