椭圆轨道交会、悬停与绕飞的全状态反馈控制

面向大椭圆轨道航天器交会对接、编队伴飞以及在轨操控等空间应用的需求,对大椭圆轨道上航天器间的相对运动进行了分析与建模,采用幂级数法分别在脉冲推力和常值推力作用两种情况下对系统进行了近似求解。通过对系统解的变换以及对系统状态的重构,给出了大椭圆轨道上的三种交会制导律。脉冲推力作用假设下的脉冲制导类似近圆轨道的Hill制导方法。常值推力作用假设下的全状态反馈制导律则在交会制导、相对悬停和循迹绕飞控制的过程中实现了对相对位置和相对速度的同步控制。通过构造新的系统状态,改进的变系数全状态反馈制导律提高了相对速度的制导精度,降低了相对制导过程中的最大轨控加速度。三种制导律的制导效果通过数学仿真进行了校验和比较,文中给出的方法实现了椭圆轨道上相对交会制导、悬停保持和循迹绕飞控制。     

Aerodynamic Design Methodology for Blended Wing Body Transport

This paper puts forward a design idea for blended wing body(BWB).The idea is described as that cruise point,maximum lift to drag point and pitch trim point are in the same flight attitude.According to this design idea,design objectives and constraints are defined.By applying low and high fidelity aerodynamic analysis tools,BWB aerodynamic design methodology is established by the combination of optimization design and inverse design methods.High lift to drag ratio,pitch trim and acceptable buffet margin can be achieved by this design methodology.For 300-passenger BWB configuration based on static stability design,as compared with initial configuration,the maximum lift to drag ratio and pitch trim are achieved at cruise condition,zero lift pitching moment is positive,and buffet characteristics is well.Fuel burn of 300-passenger BWB configuration is also significantly reduced as compared with conventional civil transports.Because aerodynamic design is carried out under the constraints of BWB design requirements,the design configuration fulfills the demands for interior layout and provides a solid foundation for continuous work.     

四阶微分方程两点边值问题解的存在唯一性

讨论了非线性四阶常微分方程y（4）=f（x,y,y’,y”,y”’）在混合两点边值条件y’（a）=0,y”（a）＋y”（6）=0,y（b）=0,y”’（b）=0或y’（a）=0,y”＇（a）＋y＇”（b）=0,y（b）=0,y”（6）=0下,解的存在唯一性。其中f在[a,b]×R4上连续且满足Lipschitz条件。并在推广后的Lipschitz条件与Banach压缩映射原理基础上,得到一些新的存在唯一性结果。     

随机事元及其在解决矛盾问题中的应用

通过建立随机事件的可拓模型,给出了随机事元、随机事元的概率、随机事元集和可拓事件集的概念。在对随机事元可拓性研究的基础上给出了随机事元的拓展概率的概念。利用随机事元的多特征性,讨论了当随机事元的某一个或若干个特征及其量值改变时,由可拓变换的传导性,将在随机事元的各要素之间产生传导效应,从而导致随机事元概率的改变,进而为涉及随机事件发生的概率的矛盾问题的解决提供了有效途径。最后举例说明了利用概率的改变处理矛盾问题的方法。     

大型运输机机翼气动外形设计研究

以大型运输机机翼为研究对象,针对其机翼高低速气动外形设计目标,采用改进的Takanashi正反迭代、余量修正方法,进行超临界翼型的大型运输机机翼气动外形反设计.采用优化方法进行三维机翼低速增升装置的优化设计,通过将气动外形优化设计方法与反设计方法相结合以及分布式并行计算等手段来提高设计效率.结果表明,该设计方法是合理可行的.     

基于矩阵奇异值理论的颤振分析新方法

根据颤振分析的基本概念,提出了一类基于矩阵奇异值理论的颤振分析新方法。该方法的特点是,以计算颤振矩阵最小奇异值或条件数的倒数来直接搜索颤振临界点。证明了这两个指标在颤振临界点处的等价性。根据指标在颤振临界点附近取极小值的特点,编制了相应的算法,在确定颤振临界点时无需计算颤振特征根,避免了“窜支”问题,从而减少了人工干预,提高了计算自动化程度。数值算例结果表明,采用该方法计算得到的颤振临界速度和颤振频率与p-k法计算结果的精度相当,且两个指标对应的计算结果一致,验证了其等价性。     

基于MCIS技术的相位差变化率单站无源定位

针对最大似然估计（ML）方法求解测相位差变化率单站无源定位问题计算量大、定位慢的问题,本文提出一种利用蒙特卡洛重要性抽样技术（MCIS）高精度、低复杂度的估计方法。根据Pincus定理推导出ML问题的近似全局解,利用重要性抽样（IS）技术构建符合高斯分布概率密度（PDF）的重要性函数,作为样本选取的依据,通过逆变换采样获得样本集,统计样本均值直接得到辐射源位置估计结果。MCIS方法简单易实现且运算量低,能够克服传统ML估计多维网格搜索耗时较长的缺陷,而且对目标位置初始估计误差有较低的敏感性。实验结果表明,MCIS算法在相同测量噪声水平下,定位精度优于EKF、NLS算法,有效减小了初始化估计误差对算法定位精度的影响,也进一步讨论分析了算法参数和不同观测条件对定位性能的影响。     

基于逆最优问题的最优制导律及特性分析

在典型的能量最优制导律基础上,将制导律的2个特征根从有限的点/线区域扩展到所有可能的正实根区域,进而提出制导律中的逆最优问题。详细讨论了逆最优问题中性能指标加权矩阵的构造过程,给出了加权矩阵和Riccati矩阵的计算公式;将控制权矩阵选为time-to-go的负n次幂的形式,对加权矩阵的求解进行了举例说明。对8组不同的特征根研究结果表明,尽管每一对可能的特征根取值都能找到最优解释,但这并不能保证与其对应的制导律都能达到与典型能量最优制导律类似的制导性能,特征根取值越靠近典型能量最优制导律,则相对应的制导特性也越接近。     

Design and analysis of an extended mission of CE-2: From lunar orbit to Sun–Earth L2 region

Chang’E-2 (CE-2) has firstly successfully achieved the exploring mission from lunar orbit to Sun–Earth L2 region. In this paper, we discuss the design problem of transfer trajectory and at the same time analyze the visible segment of Tracking, Telemetry & Control (TT&C) system for this mission. Firstly, the four-body problem of Sun–Earth–Moon and Spacecraft can be decoupled in two different three-body problems (Sun–Earth + Moon Restricted Three-Body Problems (RTBPs) and Earth–Moon ephemeris model). Then, the transfer trajectory segments in different model are computed, respectively, and patched by Poincaré sections. The full-flight trajectory including transfer trajectory from lunar orbit to Sun–Earth L2 region and target Lissajous orbit is obtained by the differential correction method. Finally, the visibility of TT&C system at the key time is analyzed. Actual execution of CE-2 extended mission shows that the trajectory design of CE-2 mission is feasible.     

Broad-search algorithms for the spacecraft trajectory design of Callisto–Ganymede–Io triple flyby sequences from 2024 to 2040, Part II: Lambert pathfinding and trajectory solutions

Triple-satellite-aided capture employs gravity-assist flybys of three of the Galilean moons of Jupiter in order to decrease the amount of ΔV$\mathrm{\Delta }V$ required to capture a spacecraft into Jupiter orbit. Similarly, triple flybys can be used within a Jupiter satellite tour to rapidly modify the orbital parameters of a Jovicentric orbit, or to increase the number of science flybys. In order to provide a nearly comprehensive search of the solution space of Callisto–Ganymede–Io triple flybys from 2024 to 2040, a third-order, Chebyshev's method variant of the p-iteration solution to Lambert's problem is paired with a second-order, Newton–Raphson method, time of flight iteration solution to the _V∞$V_{\infty }$-matching problem. The iterative solutions of these problems provide the orbital parameters of the Callisto–Ganymede transfer, the Ganymede flyby, and the Ganymede–Io transfer, but the characteristics of the Callisto and Io flybys are unconstrained, so they are permitted to vary in order to produce an even larger number of trajectory solutions. The vast amount of solution data is searched to find the best triple-satellite-aided capture window between 2024 and 2040.