不同重力条件下水平管内气液两相流流型及空隙率分布的数值模拟

基于VOF方法建立了不同重力条件下水平管内气液两相流动的三维非稳态数学模型并进行了数值求解, 研究了10-4g₀, 0.17g₀, 0.38g₀, 1g₀ (g₀=9.8m·s-2)四种重力条件下水平管内气液两流型及变化规律, 比较了不同重力条件下管内截面空隙率的分布和波动规律. 研究结果表明, 该模型能够正确预测不同重力条件下水平管内气液两相流的流型、截面空隙率和滑速比等重要参数; 同一气液两相表观速度工况下, 随着重力水平的升高, 气相更容易在水平管的上部积聚合并, 致使流型发生变化, 同时, 气液两相滑速比增大, 截面空隙率波动峰值的平均值下降, 波动频率降低; 而随着气液两相表观速度的增大, 两相混合工质内惯性力作用也随之增强, 这将削弱重力变化的影响.      

基于径向基函数响应面的机翼有限元模型修正

用ANSYS三维实体单元SOLID45建立机翼基准有限元模型并计算其自由振动的前6阶模态频率.用均匀设计方法将结构参数分组并分别计算各组结构参数对应的模态频率,建立高斯径向基函数响应面模型.用最小二乘法则拟合系数并检验响应面拟合精度,对基准模型的结构参数施加摄动量建立待修正有限元模型.用响应面模型和基准模型计算所得模态频率的相对误差建立适应度函数的表达式,将混沌搜索机制引入粒子群算法对结构参数的摄动量进行寻优计算,搜索所得优化解代入即得修正后模型,将修正后模型与基准模型在测试频段内段外的模态频率近似度进行比较,证实了修正后模型的有效性.      

双三角翼及其翼身组合的滚转运动特性比较研究

采用Navier-Stokes方程与滚转运动方程耦合计算方法,比较研究了不同后掠角的双三角翼和翼身组合体的滚转运动特性,分析了机翼前缘后掠角及细长机身对非定常滚转力矩时滞环、动态流场结构和物面瞬时压力分布的影响。研究结果表明：主翼迎风面上的融合涡能量在80°/60°双三角翼上耗散较小,而在76°/40°双三角翼上耗散严重,这是造成两模型滚转力矩稳定性与时滞特性差异的主要因素;机身对气流的扰动作用,大幅增强了滚转力矩的线性分量;机身对气流的上洗作用,增强了边条涡与融合涡吸力及其时滞性,同时加剧了主翼背风面的两涡干扰;大滚转角时机身对横流流动的干扰,使得主翼背风面压力分布的时滞差异显著增加。该研究结果有助于认识后掠角与细长机身影响双三角翼滚转运动特性的物理机理。     

飞翼布局隐身翼型优化设计

针对飞翼布局设计中气动与隐身设计矛盾更为突出的问题,采用高精度气动和隐身计算方法,建立了基于Parsec参数化方法、径向基函数（RBF）神经网络、Pareto遗传算法和松散式代理模型管理方法的翼型多目标优化设计平台。根据飞翼布局内外翼不同功能和特点,确定了内外翼翼型不同的优化设计目标和约束条件,开展了兼顾气动与隐身性能要求的翼型综合优化设计研究。结果表明：对兼顾气动与隐身性能要求的飞翼布局,内翼段翼型主要通过弯度、前缘半径、尾缘角及厚度等设计,减小低头力矩和重点方位角的雷达散射截面（RCS）均值。外翼段翼型上表面的几何形状对跨声速气动效率的影响很大,应通过上表面设计提高跨声速气动效率,重点方位角RCS均值的减小则通过下表面设计实现。某些翼型参数对气动和隐身性能均有较大影响,但作用相反,应作为综合优化设计的主要设计参数,并采用不同的优化设计策略。Pareto方法给出的前沿阵面可为飞翼布局的三维设计提供更丰富的信息。     

Integrated optimization design of smart morphing wing for accurate shape control

Smart morphing wing, which is equipped with smart materials and able to change structural geometry adaptively, can further improve aerodynamic efficiency of aircraft. This paper presents a new integrated layout and topology optimization design for morphing wing driven by shape memory alloys (SMAs). By simultaneously optimizing the layout of smart actuators and topology of wing substrate, the ultimately determined configuration can achieve smooth, continuous and accurate geometric shape changes. In addition, aerodynamic analysis is carried out to compare smart morphing wing with traditional hinged airfoil. Finally, the optimized smart wing structure is constructed and tested to demonstrate and verify the morphing functionality. Application setbacks are also pointed out for further investigation.     

概率盒全局灵敏度和活跃子空间跨层降维

航空航天仿真系统中的不确定性通常是多源的、混合的,并且系统参数的维数众多。针对高维混合不确定性量化问题,提出一种结合概率盒全局灵敏度和活跃子空间的跨层降维方法。在随机和认知不确定的概率盒表征基础上,使用不确定性缩减法分析参数的全局灵敏度继而进行参数筛选;基于输出梯度协方差矩阵的特征分解,使用活跃子空间法对参数进行降维;构造出一种概率盒表征下的参数筛选和跨层降维方法。最后以NASA多学科不确定性量化挑战问题为例,通过概率盒全局灵敏度分析进行第1层次的参数筛选,原有的21维输入参数减为13维;随后采用活跃子空间进行第2层次的参数降维,维数进一步降至一维。研究结果表明,所提出的方法能够对混合不确定性参数进行灵敏度排序,还能够有效降低模型输入参数的维度,为高维系统混合不确定性量化和进一步的优化工作奠定了基础。     

三角翼机翼摇滚主动控制多学科耦合数值模拟

细长机身和大后掠机翼气动构型的飞行器大攻角飞行时,由于缺少横向阻尼,易发生以绕体轴滚转振动为主的摇滚运动,飞行安全受到严重威胁。针对三角翼摇滚问题,采用动网格技术,建立了气动、运动和控制多学科耦合的数值模拟方法。通过耦合非定常Navier-Stokes方程、刚体运动方程和经典控制律,采用控制面差动偏转的方式对三角翼摇滚主动控制过程进行了数值模拟,并分析了不同控制状态下三角翼受控滚转的运动特性。在来流马赫数为0.3的条件下,实现了80°后掠三角翼摇滚现象的有效控制。     

小展弦比飞翼大迎角状态空间气动力建模

现代军机大迎角区域性能对空中作战优势的建立有着重要影响,针对大迎角区域建立合适的数学模型对于飞行仿真、稳定性分析和控制律设计都有着重要作用,并且对于解决飞行安全问题和研究飞机失速和尾旋问题具有重要意义。针对小展弦比飞翼标模的迟滞特性,对状态空间模型进行改进,应用大振幅强迫振荡试验数据建立了非定常气动力模型,采用风洞典型机动模拟试验验证状态空间模型的有效性和适用性。结果表明：本文发展并改进的状态空间模型能够准确预测小展弦比飞翼标模不同机动下的非定常气动力特性,具有较强的工程实用性。     

基于CFD方法的大型客机高速气动设计

大型客机高速气动设计需融入型号的设计经验、准确的数值分析方法以及高效的全局优化流程。将型号研制积累的设计经验及准则与现有数值计算工具、优化算法和计算机硬件资源相结合,探索发展了基于CFD的大型客机气动优化设计综合方法,该方法系统综合全局优化与局部寻优、人工经验与数值优化、参数化方法和参数控制以及自动化网格生成等方法和技术,大幅提升了气动设计效率。同时,完善了工程中实用的大型客机高速气动设计方法和流程,设计过程中融入了气动、结冰、静气动弹性等多专业的综合约束,反映了机翼设计多学科综合的本质特征,有助于形成综合最优的设计方案。以大型客机的超临界机翼优化设计为例,叙述了其在高速气动设计工作中的应用。     

机翼结构重量预测的多学科分析优化方法

为了克服现有机翼结构重量计算方法的局限性,提出一种基于多学科分析优化的机翼结构重量计算方法。以客机机翼为例,阐述整个计算流程。计算流程的关键步骤包括机翼外形和结构参数化建模、气动分析模型自动生成与外形优化、结构有限元模型的自动生成和结构优化。应用CAD软件CATIA的二次开发方法,实现机翼外形几何模型、结构布置几何模型和气动分析模型的自动生成;应用MSC.Patran的PCL编程技术,实现结构有限元模型的自动生成;应用等效刚度和等效强度方法,提高结构有限元模型自动生成的稳健性,缩短结构分析和优化的计算时间;应用多学科集成和优化技术,建立机翼结构重量预测的计算平台,实现整个计算过程的自动化。算例表明这种方法稳健、有效,可快速地分析机翼外形参数与结构重量之间的关系,分析不同展向载荷分布和不同选材方案对机翼结构重量的影响。