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

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

基于纯比例导引的拦截碰撞角约束制导策略

拦截碰撞角约束制导是当前导弹制导研究的关键问题之一。首先基于理想比例导引（IPN）律拦截非机动目标的解析解,推导了纯比例导引律（PPN）拦截固定目标的解析解,得到了弹目相对距离、制导指令加速度和导弹前置角的显示表达式,并进一步得到了拦截碰撞角与弹目相对运动状态和比例导引系数之间的解析表达式。其次,基于该解析表达式,提出了基于PPN的拦截碰撞角约束制导策略（PPNIACG）,并探讨了在铅垂面内进行落角约束打击和水平面内进行拦截碰撞角约束打击的2种实现方式。最后,以弹道成型制导律（TSG）和最优碰撞角约束制导律（OIACG）为参考,通过数值仿真算例,对PPNIAC的拦截性能进行了对比分析,验证了所提出制导策略的有效性和正确性。     

Optimization and verification of free flight separation similarity law in high-speed wind tunnel

Based on the similarity of separation time, a similarity law optimization method for high-speed weapon delivery test is derived. The typical separation state under wind load is simulated by the numerical method. The real separation data of aircraft, separation data of previous test methods, separation data of ideal wind tunnel test of previous methods, and simulation data of the proposed optimization method are obtained. A comparison of the data shows that the method proposed can improve the performance of tracking. Similarity law optimization starts with the development of motion equations and dynamic equations in the windless state to address the problems of mismatching between vertical and horizontal displacement, and to address the problems of separation trajectory distortion caused by insufficient gravity acceleration of the scaling model of existing light model. The ejection velocity of the model is taken as a factor/vector, and is adjusted reasonably to compensate the linear displacement insufficiency caused by the insufficient vertical acceleration of the light model method, so as to ensure the matching of the vertical and horizontal displacement of the projectile, and to improve the consistency between the test results of high-speed projection and the actual separation trajectory. The optimized similarity law is applicable to many existing free-throwing modes of high-speed wind tunnels. The optimized similarity law is not affected by the ejection velocity and hanging mode of the projectile. The optimized similarity law is suitable not only for the launching of the buried ammunition compartment and external stores, but also for the test design of projectile launching and gravity separation.     

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

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

基于目标距离的末制导雷达最优搜索图模型

深入研究了反舰导弹末制导雷达目标捕捉模型,以末制导雷达第一方位搜索周期捕捉概率和全程搜索捕捉概率为约束条件,建立了以末制导雷达第一方位搜索周期扫描面积为目标函数的末制导雷达最优搜索图求取模型,解决了基于目标距离的末制导雷达最优搜索图参数求取问题,并利用多项式拟合,得到了最优搜索图参数工程实现方法,该方法易编程、计算量小、占用内存少,是一种简单、实用的中远程反舰导弹末制导雷达最优搜索图参数实现方法。     

基于鲁棒变结构控制方法的非线性导弹自动驾驶仪设计

对一类非线性不确定导弹控制系统,采用自适应方法进行了基于滑模理论的变结构控制器设计。该方案使得采用反演手段产生的微分爆炸问题得到较好的回避,同时构造自适应律,消除了不确定性对系统的影响,最后通过构造Lyapunov函数方法,确保了系统的鲁棒稳定性。文末仿真结果表明了该方案的有效性。     

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

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

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

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

二维翼型结冰数值计算

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

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

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