航天器气动力辅助轨道转移轨迹优化问题研究

针对航天器共面及异面气动辅助轨道转移的轨迹优化问题,利用高斯拟谱法将原始连续两点边值问题离散化并转化为等价的非线性规划问题,应用SNOPT算法对此非线性规划问题进行求解．同时分析了对应于不同变轨角度时大气内飞行轨迹、飞行速度以及热流峰值的变化规律．     

基于飞行品质的飞机控制增稳系统参数估计

提出了一种控制增稳系统参数的设计方法,即利用低阶等效系统和飞行品质边界,通过搜索算法确定满足对应飞行品质的参数集合,再利用多种飞行品质确定的参数集合的交集和数据的插值算法,确定控制增稳系统的控制律参数。通过纵向增稳系统设计仿真表明,所提方法简化了传统设计方法通过性能指标反复选择参数和验证系统飞行品质的过程,为控制增稳系统的设计提供一种简单有效的方法。     

一种新的等效拟配初值计算方法

在简要分析传统等效拟配初值选取方法不足的基础上,提出了一种基于Matlab信号处理工具箱提供的invfreqs函数的初值计算方法。整个计算过程不需要人工参与,所得拟配初值接近优化的终值,与传统初值选取方法相比,提高了拟配效率和成功率。仿真结果表明,该初值计算方法具有良好的工程应用价值。     

结构可靠性分析的拟蒙特卡罗方法

由于蒙特卡罗(MC)方法具有程序结构简单,收敛速度与问题维数无关等优点,故其在结构可靠性分析中得到了广泛应用。但是对于小失效概率等问题,计算效率低这一主要缺点限制了该方法的应用范围。为了解决MC方法存在的问题,通过引进单位超立方体中不同的低偏差点集代替伪随机数序列,并结合了重要抽样技术建立了结构可靠性分析的拟蒙特卡罗(QMC)方法。该方法不但可以大幅度减少抽样点数目,还能够得到确定性的估计值, 避免传统MC方法只能得到概率意义下误差的缺陷。通过数值算例可以看出本文方法具有较高的计算精度和效率,因此适用于结构可靠性分析问题。     

FY-3气象卫星中分辨率光谱成像仪结构仿真分析

用ANSYS软件仿真分析了风云三号（FY-3）气象卫星中分辨率光谱成像仪的结构。根据建立的成像仪有限元模型,进行了静力学、模态和谱分析,给出了X、Y、Z向加速度载荷下消旋机构的应力和变形、1-20阶模态频率,以及各向正弦振动谱作用下各部件的应力和变形。结果表明：有限元法与实验数据基本吻合,可用于成像仪设计的优化。     

稀疏拟谱最优控制法求解Goddard火箭问题

提出了一种新的基于直接转化法的求解基于常微分方程(ODE)和微分代数方程(DAE)的最优控制问题的数值方法.该方法通过Legendre-Gauss拟谱法同时离散化状态变量和控制变量,把最优控制问题转化为一个非线性规划问题,并利用改进的多相处理方法避免优化无控段,同时基于稀疏矩阵探索其一阶导数信息.数值结果表明,与传统的直接转换法相比,该方法是一种通用高效的精度较高的ODE/DAE最优控制直接数值求解法.最后,从工程观点出发,应用该方法成功求解了终端自由有路径约束的奇异最优控制问题Goddard火箭问题.     

马赫数4二维直管内拟似冲击波实验和数值研究

通过计算流体力学数值仿真和高速彩色纹影照片显示技术对马赫数4拟似冲击波进行研究。实验是在日本室兰工业大学的压力—真空型超声速风洞中进行的。数值计算是采用三阶精度的QUICK格式和Spalart-Allmaras湍流模型进行的。用实验验证CFD模型,计算结果与实验显示出较好的一致性。可以得出结论,从目前的CFD模型获得的计算结果是精确的。因而,实验很难获得的拟似冲击波流场内部一些流动量可以通过数值仿真结果进行分析。另外,推测了马赫数4拟似冲击波非对称的原因。     

基于主动约束层阻尼控制的大型挠性航天器动力学模型

以大型挠性航天器简化模型——中心刚体和固接板混合系统为对象,同时考虑板表面粘贴的主动约束层阻尼结构,根据拟坐标拉格朗日方程,建立了包含所有高阶小量的完整动力学模型。建模时,采用混合坐标描述中心刚体和主动约束阻尼板的运动,用Golla-Hughes-Mctavish模型描述阻尼材料的复模量。基于该动力学模型可进行空间柔性板振动的主动约束阻尼控制器设计和大型挠性航天器的姿态动力学研究。     

交会对接远距离导引精度分析

为确定交会对接任务地面远距离导引控制可达到的精度,分析了影响精度的误差因素;提出了利用简化动力学模型、采用拟平均根数法和协方差法进行误差传播计算的方法;推导了误差传播和误差分析数学模型;进行了仿真计算.结果表明:远距离导引终点精度主要由传播误差决定,在一定初始条件下,当外推时间小于6 h时,终点位置精度可达到2~5 km,速度精度可达到2~3 m/s;误差传播中,沿迹和径向误差因素占主要成份,且随外推时间增加沿迹误差影响逐渐增大;终点精度的提高应从抑制沿迹、径向误差着手.提出的地面导引控制精度分析方法综合考虑了各种误差因素,计算便捷,适用于工程方案设计阶段的精度分析与估计.      

扑动翼型的低雷诺数气动特性分析

通过求解引入拟压缩项的不可压Navier-Stokes方程,数值模拟了绕扑动翼型的低雷诺数非定常流动。针对厚度在4%-12%之间的NACA对称翼型,分析了翼型厚度等参数对扑动翼型气动特性的影响。在低雷诺数条件下,对于纯俯仰运动,随着翼型厚度的减小,平均阻力系数也变小。而对于纯沉浮运动,发现翼型厚度对气动特性的影响和俯仰运动有很大的差别,平均阻力系数随着翼型厚度的减小而变大。通过对沉浮运动一个周期流线图的分析,认为这是翼型前缘涡的影响造成的。由于前缘涡的影响,翼型厚度增加,平均压差阻力系数变小,甚至会出现负值。雷诺数的影响研究表明,随着雷诺数的增加,扑动翼型的阻力系数减小的趋势越缓慢。