基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

战斗机嵌入式训练系统中的智能虚拟陪练

智能化"实虚"对抗是现代先进战斗机嵌入式训练系统的重要功能需求。自主空战决策控制技术在未来空战装备发展中扮演关键角色。将当前的功能需求和发展中的技术结合起来,得到了空战智能虚拟陪练的概念。先进控制决策技术的引入使得智能虚拟陪练能够帮助飞行员完成复杂的战术训练,而训练中真实的对抗场景为技术的验证提供了理想的环境,大量的训练数据为技术的持续迭代优化提供了保障。作为可学习和进化的空战战术专家,智能陪练在人机对抗和自我对抗中不断优化,当其具备与人相当甚至超越人的战术能力时,可应用于未来的无人空战系统。智能虚拟陪练需要具备4项基本能力：智能决策能力、知识学习能力、对抗自优化能力和参数化表示能力。对其包含的关键技术进行了分析,提出并实现了一个基于模糊推理、神经网络和强化学习的解决方案,展示了其各项基本能力及目前达到的空战水平。未来更多的模型和算法可在智能虚拟陪练的框架中进行验证和优化。     

外挂物投放问题的二维非结构动态网格技术研究

以动态网格守恒型方程为控制方程,在二维非结构网格基础上发展了弹簧式光顺法和局部网格重构法两种动态网格变形技术。结合运动轨迹计算方程,将动态网格法应用于飞机外挂物投放问题中,成功模拟了导弹自由投放的整个过程,获得了详细的包括弹体运动速度、运动轨迹和受力情况等在内的丰富的流场信息。成功的应用不仅验证了动态网格法的可行性;而且也表明,动态网格法的进一步完善与发展将成为解决外挂物投放问题的有力工具。     

超声速变马赫数风洞流场参数线性变化验证

旋转喷管型面使超声速变马赫数风洞在单次运行过程中可连续调节实验区的马赫数,便于研究飞行器的机动过程、进气道起动过程中的气动问题。在控制喷管型面旋转过程中,流场参数能否线性变化是衡量超声速变马赫数风洞性能的一个重要指标。分析变马赫数风洞实验区流场参数的线性变化规律,利用弹簧光顺的动网格技术建立数值仿真模型,验证喷管位于马赫数3.041~3.215 范围所对应的位置时,实验区流场参数是否满足线性变化规律。结果表明：通过对喷管型面旋转的控制实现了风洞实验区流场参数的线性变化,动态计算结果与预期实验区流场参数线性变化规律吻合良好;在不同加速度的流场参数线性变化过程中,各时刻实验区的平均参数与预期参数之间的偏差均小于0.13%。     

基于振动传递符号有向图的齿轮箱嵌入式传感器优化配置模型与算法

针对嵌入式条件下振动信号的采集方式、传递路径与传统方式不同的特点,采用有限元仿真方法,分析了嵌入式传感下齿轮箱故障振动的传递特性,在此基础上建立了基于故障振动传递符号有向图的优化配置模型,并研究了基于粒子群优化的求解算法.案例仿真结果验证了上述方法的有效性,为在齿轮箱内部合理嵌入传感器提供了一种可行途径.      

薄壁结构的加筋布局优化设计

采用拓扑优化方法研究了薄壁结构的加筋布局优化问题.考虑到航空航天领域薄壁加筋结构的大量使用与结构形式的复杂性,提出了一种适用于有限元自由网格剖分的加筋设计新方法--几何背景网格法.该方法一方面通过定义几何背景网格的大小,实现了加筋设计域内任意离散网格沿加筋高度方向的布局参数化定义;另一方面,通过背景网格曲线坐标系下的定义,实现了三维曲面薄壁壳结构的加筋布局设计.采用该方法分别对平面结构与曲面薄壁结构进行了以结构刚度最大为目标的加筋优化设计.数值结果表明,所提出的方法能有效获得合理的加筋布局.     

压力载荷下的结构拓扑-形状协同优化

 压力载荷作用下的结构轻量化设计是工程中的常见问题,由于压力加载面的可设计性,现有以固定载荷为基础的拓扑优化技术不能很好地处理这类问题。直接采用CAD参数化样条或B样条曲线描述压力加载面,通过拓扑和形状变量的联合优化满足了工程实际对结构轻量化与边界的功能性与光滑性设计要求。同时,为了避免结构边界形状变化时有限元网格刷新引起的定义拓扑伪密度变量的困难,用所提出的背景网格和密度点技术实现了每一步单元密度设计迭代结果的自动传递,并采用网格变形技术实现了形状设计变量灵敏度分析。采用4个数值算例验证了方法的有效性,其中发动机承力框架的设计结果充分说明该方法在航空结构设计中的重要应用价值。     

二维翼型结冰数值计算

翼型结冰计算是随着结冰时间步长变化的准定常过程,因此在每个时间步长后都要更新各个模块的计算工作,如此反复直到达总结冰计算时间.在各模块计算中,分别采用了泊松方程法生成网格、求解雷诺时均方程计算流场、拉格朗日追踪法求解水滴轨迹、基于Messinger模型计算冰形、采用变距滤点法重构边界.引入动态链接库技术将各模块集成为结冰计算平台,并将三种冰形(楔形冰、槽状冰及混合冰)计算结果与国外软件及试验数据进行了对比.      

一种三维结构网格并行算法

基于三维结构网格,以MPICH2作为并行开发环境,研究了对等模式下流场的并行计算问题。在三维结构网格已经生成的前提下,针对计算区域的划分,网格输出格式和数据通信,提出了一种并行算法。算法的核心思想是将需要通讯的其他区域的单元按区号从1到n排列,并获得本区域通信单元的I-J-K编号与其他区域通讯单元的排列序号的映射关系。最后,在IBM BladeCenter HS22刀片机平台上进行了算例验证,结果表明算法具有很好的可靠性与通用性。     

F-35战斗机气动及隐身特性分析

以F-35战斗机为研究目标,对目标体进行三维外形重建,并对重建后的模型进行气动及隐身特性的计算和分析,首先采用了基于飞机三视图进行轮廓线提取来重建F-35全机理论外形的方法,引入了基于草图跟踪的CATIA三视图的标定,大大地提高了模型重建的精度。其次采用非结构网格对F-35的计算区域进行网格划分,采用EULER方程完成了F-35在亚声速、跨声速及超声速等飞行条件下的流场计算,分析了不同状态下的升力、阻力和大迎角气动特性。最后利用曲面像素法对F-35全机高频雷达目标特性进行了计算,提出了一种基于IGS数据转换格式的隐身计算网格生成方法,比较了不同俯仰角及方位角下的RCS特性曲线,分析了对RCS影响比较大的部件。