三维物体表面摩阻的有限元算法

在有限元法的框架内，给出一种从动量守恒方程出发计算物面摩阻的新方法，并用数值试验进行了验证。     

导叶伺服动力学系统通用建模方法

针对现有导叶伺服系统模型跨平台联合仿真速度慢、通用性差的问题，提出了一种基于动力学特性方程的、部件通用性 强的导叶伺服系统通用建模方法。采用部件建模法从底层方程到部件级模型，再到顶层系统，形成层次化、模块化的建模架构，使 原理与工程应用有机融合；通过对导叶伺服系统各组成部件进行工作机理分析，建立导叶伺服系统的AMESIM模型，推导各组成 部件的动力学特性方程，开发Simulink部件模型库；与AMESIM模型进行对比验证，构建Simulink导叶伺服动力学系统机理模型。 结果表明：系统模型稳态误差不大于0.12%，斜坡响应跟踪误差为3.7%，对于幅值为作动筒位移5%的不同频率正弦指令信号，导 叶伺服闭环系统带宽不小于8 Hz。该导叶伺服动力学系统模型具有较好的准确性、稳定性和动态特性，是可行和有效的，具有明 显的推广应用价值和跨平台多学科联合仿真优势。     

基于弹载捷联惯性导航系统精确导航的双欧拉全姿态方法

垂直发射的地空导弹在大机动飞行时,其四元数法姿态角输出存在奇异性。为解决该问题,提出了正、反欧拉角微分方程切换的双欧拉全姿态解算方法。对正、反欧拉角微分方程进行了详细推导,论述了正、反欧拉角微分方程精华区与奇异区的倒置关系。在此基础上,建立了正、反欧拉角的转换公式。根据正欧拉角的俯仰角对正、反欧拉角无缝切换,有效解决了俯仰角在±90°附近时姿态奇异性问题。仿真和半物理试验的结果表明:双欧拉全姿态算法解算的实际姿态与理论姿态保持一致,在奇异点附近的姿态角仍能正常输出,未出现突变现象,该方法得出的导航精度比四元数法得出的导航精度提高了9.81%。     

超高速撞击中靶后碎片云的内外边界方程研究

研究超高速撞击防护靶形成的碎片云的形状模型对航天器防护结构设计具有重要意义。文章用数值仿真的方法研究靶后碎片云的内外边界模型:将碎片云根据其材料来源分为靶板碎片云和弹丸碎片云,分别根据各碎片云的分布特点建立碎片云形状的内外边界方程;在双纽线边界方程的基础上,修正为由双纽线方程和圆弧方程组成的边界方程。对比整体碎片云的适用率表明,由双纽线方程和圆弧方程组成的边界方程比单一的双纽线边界方程描述碎片云实际分布更为贴合。     

中空纤维膜机载制氮装置的数学建模分析

建立了考虑浓差极化现象的微分方程数学模型,并用正交配置法求解,对中空纤维膜机载制氮装置(OBIGGS)进行了分析,并对部分状态点的计算结果进行了实验验证.结果表明:富氮气体(NEA)中氧气的质量分数随着进气温度的升高而降低,在达到最小值后又呈上升趋势;在进气与排气压力差保持不变的情况下,随着中空纤维膜排气压力的下降,中空纤维膜的富氮气体质量流量逐渐增加;随着中空纤维膜丝长度的增加,丝内气体的压降和富氮气体质量流量均有所增加;中空纤维膜空气分离制得的富氮气体质量流量越大,则所需进气的质量流量越大,且富氮气体氧气质量分数越高.      

基于三维非结构网格的连续伴随优化方法

基于非结构网格,采用Euler方程、连续伴随方法以及自由形面变形技术(Free-form deformation,FFD)结合序列二次规划(Sequential quadratic programming,SQP)优化算法,开发出了三维飞行器的气动外形优化代码。应用Fortran 90语言,动态存储数组和指针,节省了内存。采用Roe形式的二阶迎风格式和LU-SGS隐式时间离散格式求解连续伴随方程,同时采用一种新的参数化坐标求解方法,使得FFD技术更加鲁棒。和有限差分方法求得的梯度进行比较,验证了连续伴随方法所求梯度的正确性,实现了ONERA M6机翼的减阻优化,证明本文开发的优化代码能够实现三维气动外形的减阻优化。     

DLR-F6翼身组合体的高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6翼身组合体的高阶精度数值模拟研究。主要目的是确认WCNS模拟跨声速典型运输机构型的能力。采用粗、中、细3套网格开展了网格收敛性研究,从气动特性、压力分布、表面流态等方面研究了网格密度对DLR-F6翼身组合体气动特性的影响;采用中等网格开展了来流迎角对气动特性的影响研究。通过与试验数据、CFL3D软件和TRIP软件计算结果的对比,表明网格密度主要影响激波位置和压差阻力系数,同时对翼身结合部分离区大小有一定影响;采用高阶精度计算方法显著提高了气动力系数的模拟精度,力矩系数数值模拟结果与试验的差异有待进一步分析。     

基于风洞试验的非定常气动力微分方程建模方法

为准确描述大攻角非定常运动中非定常气动力的气动特性,研究了以微分方程为基础的非定常气动建模方法。在某战斗机大振幅偏航、滚转单自由度及偏航-滚转耦合运动风洞试验的基础上,对偏航力矩、滚转力矩建立微分方程模型,分析研究不同频率、同一侧滑角情况下,气动力的迟滞时间;然后对高频运动下的气动力进行拟合,进而建立高频运动下的气动力微分方程模型。研究表明:模型可精确预测不同机动下非定常气动特性;建模方法具有较强的工程可行性。     

一组新的四元数轨道要素建模方法

针对传统以欧拉角为参数的轨道要素的奇异性、不确定性以及计算效率低等问题,提出了一种新的四元数轨道要素。建立了新轨道要素与经典轨道要素,以及新轨道要素与惯性系下位置、速度的相互转换关系,推导了基于新轨道要素的高斯摄动方程。以J₂项摄动下的轨道推演为例进行仿真验证,结果表明新轨道要素不仅在圆轨道与赤道平面轨道处不存在奇异性和不确定性,而且由于新轨道要素不涉及三角函数运算,新高斯摄动方程积分效率和计算精度明显提高。     

基于傅立叶平均法下的连续小推力动力学分析

通过将小推力展开为偏近点角的傅立叶级数,并对高斯摄动方程在一个轨道周期上的平均,将原方程的推力转化为仅由14个傅立叶系数表示的控制变量。仿真计算表明,平均化后的高斯方程使计算量与牛顿积分相比显著减少,且对小推力而言有足够的精度。对利用平均化后的高斯方程计算轨道根数时产生误差的原因进行了研究,并进一步分析小推力的范围和小推力近似表达式对上述误差的影响,为今后小推力下非开普勒轨道动力学分析提供了理论依据和参数。