基于CFD和CSM耦合的通用静气弹分析方法

提出了一种适用于有限元精细化建模的流固耦合插值点选择方法,通过RBF(径向基函数)方法实现流固耦合面的数据交换,实现了基于CFD/CSM(computational fluid dynamics/computational structural mechanics)耦合的通用非线性静气弹分析方法。以HIRENASD(high Reynolds number aero-structural dynamics)风洞试验模型为验证对象,数值结果很好地与风洞试验结构变形、气动压力分布吻合,验证了所发展非线性CFD/CSM耦合静气弹求解器的精度。详细研究了HIRENASD模型在大迎角(AOA)流动下的静气动弹性特性,以及该模型弹性变形对机翼气动特性影响规律。研究表明:HIRENASD弹性模型变形后其升力小于刚性模型;在小迎角范围内刚性、弹性模型升力差随迎角增大呈线性增长;当迎角大于4°后,升力差先减小后基本保持不变,呈非线性关系。      

基于代理模型的大型民机机翼气动优化设计

先进的气动优化设计思想与方法,对于提升大型民机气动与综合性能具有至关重要的意义。探讨了大型民机超临界机翼气动优化设计的基本准则和要点,并结合代理优化算法,提出了一种面向工程应用的多轮次高效全局气动优化设计方法。首先,通过一系列解析函数/翼型/机翼优化测试算例进行了验证。其次,将所提出的方法与人工修型相结合,开展了针对宽体客机超临界机翼的两轮气动优化设计,使其气动性能得到显著改善。最后,采用不同的雷诺平均Navier-Stokes（RANS）方程求解器对安装优化机翼的全机巡航构型进行了典型状态的气动性能综合评估。研究结果表明,所提出的代理优化算法具有很高的优化效率、较强的约束处理和全局优化能力;将所发展的基于代理模型的多轮次气动优化设计方法与人工修型相结合,能够获得满足设计要求的气动外形,验证了该方法在大型民机超临界机翼气动设计中的有效性和工程实用性。     

基于参数化分析的柔性后缘优化设计

介绍了一种自适应柔性机翼后缘变形控制优化设计方法,该方法采用了由单个驱动器驱动的单块式分布柔性结构。该方法的关键在于设计一种合适的柔性后缘结构拓扑形式,使该自适应柔性机翼后缘能够达到精确的变形。阐明了柔性结构拓扑优化设计的数学模型,针对自适应柔性机翼后缘概念设计提出了基于参数化分析的拓扑优化设计方案,并进行了工程化圆整和尺寸优化。最终,通过非线性有限元分析和功能测试,验证了本文提出的柔性后缘优化设计方法的合理性。     

含有SMA弹簧驱动器的可变倾斜角翼梢小翼研究

 针对传统翼梢小翼在非设计状态减阻效果不佳的缺点,提出一种含有形状记忆合金(SMA)弹簧驱动器的变体翼梢小翼结构,它能根据飞行状态主动调整小翼的倾斜角,实时优化飞机的阻力特性.采用力-热-应变耦合法设计了所需的SMA弹簧驱动器,并通过有限元仿真与风洞试验验证了变体翼梢小翼的变形能力,最后初步研究了变体翼梢小翼的闭环控制方法.研究结果表明,在飞机的起飞阶段(自由来流流速为26 m/s,迎角为3°),变体翼梢小翼的倾斜角能在1 min内自主完成预定变化过程,倾斜角的最大变化量为23°,控制精度的最大误差为12%,各项指标均符合设计要求.     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

几何非线性机翼本征梁元素模型的高效化改进

 采用Hodges等提出的时间-空间离散化的几何精确非线性本征梁通用模型处理柔性机翼结构动力学问题时,当离散化的节点数增大时,该方法的未知数数量成倍地增长,而且方程组是严重病态的,因此数值模拟计算的速度非常缓慢。针对机翼中最常见的悬臂梁结构,根据空间离散化的边界条件,提出了空间缩聚法把空间离散差分方程缩聚为常系数矩阵格式,得到了只与时间相关的微分方程组,进一步推导得到了该方程组的雅可比矩阵,因而大大减少了方程组的数量以及求解过程的循环和迭代步数。采用Gear方法分别求解了原始的本征梁元素模型和本文提出的缩聚模型,结果表明空间缩聚模型在相同条件下可提高运算速度约5.1倍,而且对不同类型的外载荷都具有较好的通用性、稳定性和高效性。     

一种柔性蒙皮设计技术及其在后缘变弯度机翼结构中的应用

变体飞机可在飞行包线内的不同飞行状态下获得最优气动效能,增强执行多种任务目标的能力,实现变体飞机结构的关键技术之一是柔性蒙皮的设计。本文介绍一种基于可变形蜂窝的柔性蒙皮结构,并研究了这种柔性蒙皮在后缘变弯度机翼结构上的应用,通过仿真分析和全尺寸原理样件地面试验验证,该结构满足真实飞行载荷条件下对变形能力和承载能力的要求。     

跨声速机翼流场的N—S方程计算

本文给出了三维机翼粘性流场的数值模拟过程,控制方程为时均化的三维可压缩薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,湍流模型采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数湍流模型,空间离散采用中心有限体积格式,时间采用Ｒｕｎｇｅ－Ｋｕｔｔａ多齿格式进行式积分。     

待机状态下机翼结冰的快速计算方法

飞机在飞行过程中遭遇结冰气象,尤其是过冷大水滴结冰条件所导致的机翼结冰,将严重影响飞机的气动性能与操纵品质,从而导致飞行故障或飞行事故,是飞机飞行安全的重要影响因素,是飞机研制中必须重点解决的难题。虽然针对结冰已建立有多种计算仿真程序,但快速获得结冰模拟冰形,包括过冷大水滴条件下的结冰模型是工程师们一直追求的解决方案。利用前向多层神经网络,尝试建立了一种针对待机状态下的机翼结冰模型。此方法基于坐标转换原理,将获取的采用直角坐标表达的标准翼型的机翼结冰数据转换为对应的极坐标表达的数据,以飞行参数、气象参数和极坐标角度作为输入,极坐标模值作为输出训练构建出此种机翼结冰模型。使用所建立的模型,仿真预测结果表明该计算方法具备快速性和精确性,其计算精度可满足实际要求。