气动力辅助椭圆轨道转移技术研究

研究了两种共面椭圆轨道最优转移方法--基于单纯推力的双脉冲对称转移和基于气动力辅助变轨技术的协同机动方法。推导了双脉冲对称转移问题中转移椭圆轨道的求解公式,采用遗传算法求解最优变轨点位置,给出了变轨脉冲的计算方法。将气动力辅助对称转移过程分为3段：真空飞行段、大气飞行段、真空飞行段。以能量最小为性能指标,将大气飞行段轨迹优化问题转化为标准最优控制问题模型,并利用Gauss伪谱法进行求解,得到“巡航+滑行+巡航+滑行”的大气最优飞行轨迹。最后对两种转移方法进行详细比较,指出初始椭圆轨道近地点越接近大气层,气动力辅助对称转移方法比双脉冲对称转移方法越节省能量。     

星间洛仑兹力编队飞行的平衡点及零速度曲面

提出了一种新的利用星间洛仑兹力控制卫星相对运动的方法:使主星产生自旋磁场,副星带电,通过控制副星所受的星间洛仑兹力进行编队。假设主星产生的人造磁场表现为偶极子并且运动在一条开普勒圆形轨道上,副星恒定带电并在主星附近运动,同时假设星间洛仑兹力只影响副星的运动而不影响主星的运动。推导了副星在主星HCW(当地垂直当地地平坐标系)坐标系下的相对运动动力学方程。针对偶极子与HCW坐标轴X轴重合的情况下,推导了动力系统下的平衡点(并采用稳定性分析方法分析其线性化意义下的稳定性)、积分常数和零速度曲面,证明了有界相对运动的存在。最后用数值仿真验证了上述结论。     

天琴计划轨道构型长期漂移特性分析

为了实现天琴计划航天器轨道构型的长期稳定漂移并满足轨道构型的稳定性指标,文章应用优化算法搜索航天器初始最优瞬时轨道要素。建立考虑摄动情况下航天器轨道动力学模型,对导致等边三角形轨道构型发生长周期漂移的轨道要素进行了分析,并且将这些要素作为轨道设计的待优化参量,应用粒子群（ParticleSwarmOptimization,PSO）智能寻优算法进行优化,应用最终优化得到的轨道要素进行轨道构型仿真。仿真结果表明,在不进行轨道控制的情况下,4年内轨道构型呼吸角变化小于1.5°,臂长变化小于1500km,臂长变化速率小于10km/s,初步满足轨道构型的稳定性指标。     

An alternative representation of the relative motion: The local orbital elements

The study of the relative motion between two very-close non-massive objects has been usually done through Hill equations. More recent approaches use the differences of orbital elements to extrapolate relative motion. In this paper, we present another possible representation of the relative motion. This representation takes advantage of the very particular most common relative trajectory: the ellipse. Our representation induces a better temporal comprehension of the motion but it also provides the topology of the natural relative motions. Our approach is useful for circular and slightly eccentric reference orbits.     

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

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

第二类无奇点根数的北斗历书参数设计

为统一北斗三类卫星的历书拟合算法,提出基于第二类无奇点根数进行历书参数设计的方法,设计了以倾角向量变率作为摄动参数的历书参数模型,并给出新的历书模型的用户使用算法。利用覆盖2013年全年的实际在轨卫星的数值轨道进行了历书拟合试验。结果表明,地球静止轨道（GEO）和倾斜地球同步轨道（IGSO)卫星的拟合位置误差为2~4km,拟合用户距离误差（URE）约为1km,中圆地球轨道（MEO）卫星的拟合位置误差为1~2km,拟合URE约为500m。通过分析6个轨道根数和2个摄动参数全年的变化范围,对新历书模型进行量化单位和占用比特位的通信接口设计,定量分析量化单位对历书表达精度的影响。结果表明,参数截断后对位置误差的影响小于50m,对URE误差的影响小于5m。因此,历书量化误差对信号捕获以及首次定位时间带来的影响可以忽略不计。     

箱梁颤振过程中旋涡的演化规律及气动力的变化特性

利用粒子图像测速技术（PIV）观测了箱梁颤振过程中模型周围流场的旋涡特征,以模型扭转振动位移作为参考信号,采用相位平均的方法研究了旋涡规律性演化对模型周期性振动的驱动作用。风速较低时,箱梁振幅很小,其尾部风嘴附近上侧的旋涡尺度也很小,不易观测到,而下侧的旋涡尺度较大,其形状接近于圆形。当风速接近颤振临界风速时,箱梁振幅明显增大,并且模型尾部风嘴上侧的旋涡尺度也显著增大,达到与下侧旋涡尺度相匹配的程度,模型尾部风嘴上下侧旋涡的交替作用主导了结构振动直到模型振动发散。基于流固松耦合的计算策略,采用计算流体动力学（CFD）的数值方法模拟了箱梁颤振临界状态下的绕流特性,结合正交特征分解（POD）的方法研究了模型颤振时刻表面压力的空间分布特征,通过分析发现在颤振过程中箱梁表面波动压力的主控成分向迎风侧风嘴漂移。     

轨道拦截预警算法研究

针对空间轨道拦截作战模式,提出了一种拦截预警的计算方法.推导了求解轨道交点位置的表达式,然后结合拦截卫星的末端威胁距离给出了判断拦截是否发生的地心距判据和时间差判据,最后通过仿真算例验证了交点表达式以及判据的正确性,并讨论了拦截卫星的威胁距离与时间差门限的关系.     

冷轧张应力分布模型及其在辊系变形计算中的应用

以往张应力分布计算多以实验为基础,或理论模型建立过程中采用了过多假设,使张应力分布计算精度受到很大限制,针对这些情况,提出利用变分法建立单元出口张应力计算模型,并与影响函数法相结合建立了冷轧前张力横向分布计算模型及辊系弹性变形计算模型,并开发了相应的计算程序.用该程序计算了某公司单机架可逆冷轧机的单位宽轧制力、辊间压力、轧件厚度以及张应力等的横向分布并进行了分析.实例计算表明,计算结果与生产实测数据吻合较好,证明了理论模型是实用可靠的.     

基于自适应模拟退火遗传算法的最优Lambert转移

主要研究了航天器采用Lambert二脉冲变轨的优化问题。对于初始位置、目标位置和转移时间都不固定的Lambert二脉冲转移,由于多变量以及方程本身的复杂性,采用传统的优化方法效率低甚至无法求解.采用了自适应遗传算法(AGA),寻求多变量的最优解.同时结合模拟退火算法,得到了自适应模拟退火遗传算法(ASAGA),该算法既具有全局搜索能力,又改善了一般遗传算法的局部寻优能力.通过仿真,比较了遗传算法和自适应模拟退火遗传算法的寻优结果,表明两者寻求最优转移的有效性,以及自适应模拟退火算法具有更强的寻优能力.