三维制导的几何方法与鲁棒控制方法

 在考虑导弹速度和目标速度均为时变的情况下,将微分几何方法与李雅普诺夫稳定理论结合起来,提出了一种导弹三维导引规律设计新方法。对机动目标,PPNG方法的主要问题是在导引末段,它要求导弹在俯仰和偏航通道上具有较大的加速度,为了解决这一问题,本文应用基于李雅普诺夫稳定理论的鲁棒控制方法,提出了一种三维鲁棒制导算法。这种方法在导引末段不需要过大的加速度命令,不需要知道目标精确的加速度和速度方位信息。     

一种基于并联关系模型的TC4合金本构方程

 以 TC4合金流动应力σ为目标函数,通过热态物理模拟试验,系统地研究了主要锻造热力参数 (变形温度、变形速度和变形程度等 )对目标函数的影响,基于数理统计学原理,科学地分析并回归出了能综合反映锻造热力参数对材料成形性能影响的本构方程,得出了 TC4合金变形工艺参数 (变形温度 T、等效应变速率ε、等效应变ε等 )对目标函数的影响规律,从而,为 TC4合金锻造过程的数值模拟和锻造热力参数的合理制订与控制提供了基础。     

槽道流动中湍斑的数值模拟

采用Navier-Stokes方程直接数值模拟槽道流动中湍斑的形成和发展.建立高精度、高分辨率的三维耦合中心差分格式,并和时间分裂法、高阶迎风差分格式一起组成高效的高阶紧致差分方法;克服一般高阶中心差分格式不适用于边界邻域的困难,可用于包括边界邻域在内的整个计算区域;在边界上满足Navier-Stokes方程的无反射和压力条件,保证了格式的精度和一致性.通过模拟槽道流动中从壁面脉冲到湍斑的演化过程,显示出湍斑内雷诺应力和速度、压力脉动等湍流的基本特征,表明了湍斑的生成与壁面脉冲的强度和Re数大小直接相关.     

星孔型装药发动机三维两相流场的数值模拟

为了研究颗粒在星孔型装药固体火箭发动机燃烧室和喷管中的运动轨迹以及颗粒与发动机壁面的碰撞情况,针对可压两相流动,采用了高雷诺数下的k ε湍流模型和欧拉 拉格朗日两相流模型,用全速度SIMPLE方法对方程组进行求解,并用PSIC方法进行气固耦合计算。计算得出了流场内两相的速度、温度等参数的分布及多种情况下固体颗粒的运动轨迹。在燃气生成量确定的情况下,从距离喷管较近的某些位置进入流场的颗粒比较容易撞击壁面;颗粒的尺寸和局部产生的旋涡对颗粒的轨迹和碰撞也会产生较大的影响。     

遗传算法在航空发动机非线性数学模型中的应用

收敛性是航空发动机非线性数学模型的重要指标。现有发动机平衡方程迭代解法还不能保证模型大范围收敛。以某涡扇发动机为对象,采用遗传算法求解发动机非线性数学模型,将模型中的发动机平衡方程求解转换为极小值优化问题,建立了遗传算法计算模型,重点分析了采用遗传算法求解的适应度函数设计方法。数值仿真结果表明,与牛顿-拉夫逊解法相比,采用遗传算法方法可实现模型大范围收敛。     

翼型积冰的数值模拟

本文采用计算流体力学的方法对NACA0012翼型前沿的Rimeice积冰进行了预测。在流场计算中,用有限体积(FV)法对二维定常不可压粘流的时均N-S方程进行离散求解;采用四阶龙戈-库塔(RK)法求解水滴运动方程;假定水滴在与翼型相碰撞的点处完全凝结,并且冰沿着与当地翼型表面法向一致的方向增长,以此来预测积冰的形状。文中同时对迎角为4°时结冰翼型的气动特性进行了分析。     

多级火箭级间分离流动特性的数值模拟

针对某型号研制提出一种级间分离设计方案，从三维薄层近似N－S方程出发，对多种火箭级间分离发动机喷流与超声速外流形成的复杂流动进行了数值模拟，比较了栅栏构架和排焰孔两种设计方案间段流场，从流动形态分析三维排焰孔模型更为稳定，研究给出了不同分离距离下火箭的轴向力系数，发现在较近分离距离内，后体受到的作用力随距离不是单调减小，与分离喷流参数有密切关系，计算结果可作为分离段设计的参考。     

大小叶片轴流压气机转子流动特性分析

采用全三维粘性数值模拟手段，通过与按常规设计转子对比的方法，分析研究了大小叶片轴流压气机转子流动特性，分析结果表明，大小叶片转子流场在叶根亚声速区的流通能力增强，在叶尖跨声速区可以产生更加有利的激波体系，小叶片可以有效地控制叶栅槽道中气流扩散，在较高负荷和相同喘振裕度的条件下，大小叶片转子可以比按常规设计的转子在更高的压比，效率和流量下工作。     

导弹运动方程的逆解析解

针对在采用反馈线性化理论设计导弹飞行控制系统时控制器算法对导弹逆模型的要求 ,推导了导弹运动方程的逆解析解 ,包括导弹质心运动的运动学方程及其动力学方程的逆解析解、姿态运动的运动学方程及其动力学方程的逆解析解 ,详细给出了逆模型的求解过程 ,为一类飞行控制系统运用全局线性化设计奠定了基础。最后 ,结合工程应用进行了仿真计算 ,仿真结果表明逆解析解的推导是正确的     

疲劳应力变幅的断口反推研究

通过对18Cr2Ni4WA钢在不同应力比和不同载荷下的疲劳裂纹扩展速率试验,拟合得Paris关系式,然后对试验断口进行疲劳条带宽度S(微观裂纹扩展速率)测量,将疲劳条带宽度值代入拟合的Paris关系式中,反推计算断口试样的疲劳应力变幅.结果表明,疲劳条带宽度只在疲劳裂纹扩展的一定阶段与宏观裂纹扩展速率相等,选择该阶段的数据点进行疲劳应力变幅反推计算,相对误差在10%以下;数据点是否位于宏、微观裂纹扩展速率相等的阶段,可由S～yπa双对数座标S形曲线判定,并发现由于裂纹扩展速率的误差传递系数小于0.5,计算结果对裂纹扩展速率测量的相对误差不敏感.