航天器零推进剂姿态机动启发式规划方法

针对零推进剂姿态机动约束耦合,导致规划效率低的问题,提出一种基于约束评价启发式的差分进化算法。基于零推进剂姿态机动约束分析,在规划空间中建立约束评价函数,对航天器约束和控制力矩陀螺约束进行解耦分类处理。设计基于约束评价函数的启发式策略,用于改进差分进化算法的种群变异,主动控制规划进程,引导零推进剂姿态机动规划序列满足多种约束,得到安全的姿态机动路径和控制力矩陀螺框架角路径。最后,数值仿真结果验证了该方法能够实现航天器零推进剂姿态机动规划。     

解析壁面函数的可压缩效应修正研究

基于适用于不可压缩流动的解析壁面函数,针对可压缩湍流边界层特征,考虑壁面网格内对流项变化和能量方程中黏性耗散项的影响,提出一种适用于可压缩流动的解析壁面函数。基于二维超声速和高超声速激波湍流边界层干扰流动,完成了粗网格高雷诺数k-ε模型加标准壁面函数、原始解析壁面函数、可压缩修正解析壁面函数和密网格低雷诺数Lauder-Sharma k-ε模型的对比计算。结果表明:四种湍流效应模拟策略都可以准确预测壁面压力和摩擦系数,而本文发展的考虑对流项分布和黏性耗散项影响的解析壁面函数不仅消除了原始解析壁面函数的非物理振荡,而且大幅提升了壁面函数壁面热流的预测精度,与密网格解最大差异不超过40%,预测结果接近于密网格低雷诺数模型结果,而标准壁面函数和原始解析壁面函数预测的壁面热流符号相反,且数值最大差异达500%。对于Ma5斜激波边界层干扰算例,本文构造的壁面函数计算时间仅为低雷诺数模型的5%左右,相较于其它壁面模型,计算时间仅增加了1%。     

A variational Bayesian-based robust adaptive filtering for precise point positioning using undifferenced and uncombined observations

In the application of precise point positioning (PPP), especially in the dynamic mode, the classical Kalman filter (KF) usually produces a large number of estimation errors or diverges when there are gross errors in the observation data or unexpected turbulences occur in target motion state or both of them. For such problem, a variational Bayesian (VB)-based robust adaptive Kalman filtering (VB-RAKF) is proposed in this paper. This filter introduces a classification robust equivalent weight function to resist observation gross error and the inverse Wishart prior to model inaccurate process noise covariance matrix (PNCM). To improve the instantaneous accuracy of state estimation, the VB approach is used to obtain better estimations of inaccurate PNCM. Several sets of observation data collected by IGS reference stations and vehicles are employed to check the robustness and positioning accuracy of the VB-RAKF model. The results show that the VB-RAKF algorithm is more robust than the KF, and can effectively resist the gross error in observation data and control state disturbance. In the IGS reference station tests, when compared to the KF, the static positioning accuracies of the VB-RAKF in the north, east and up directions are improved by 13%, 8% and 22%, respectively, and the simulated dynamic positioning accuracies of the VB-RAKF in the north, east and up directions are improved by 19%, 9% and 21%, respectively. The in-vehicle dynamic test verifies that the VB-RAKF outperforms the KF, and shows that the VB-RAKF has better performance than the KF when dealing with observation data which has obvious gross errors, and similar performance as the KF when gross errors are small.     

Reliability estimation of rotary lip seal in aircraft utility system based on time-varying dependence degradation model and its experimental validation

With several attractive properties, rotary lip seals are widely used in aircraft utility system, and their reliability estimation has drawn more and more attention. This work proposes a reliability estimation approach based on time-varying dependence analysis. The dependence between the two performance indicators of rotary lip seals, namely leakage rate and friction torque, is modeled by time-varying copula function with polynomial to denote the time-varying parameters, and an efficient copula selection approach is utilized to select the optimal copula function. Parameter estimation is carried out based on a Bayesian method and the reliability during the whole lifetime is calculated based on a Monte Carlo method. Degradation test for rotary lip seal is conducted and the proposed model is validated by test data. The optimal copula function and optimal order of polynomial are determined based on test data. Results show that this model is effective in estimating the reliability of rotary lip seals and can achieve a better goodness of fit.     

A global empirical orthogonal function model of plasmaspheric electron content

Based on the method for establishing a global plasmaspheric model using observations from COSMIC and MetOp-A orbit determination GNSS receivers, Chen et al. (2017) obtained a global plasmaspheric total electron content product with a spatial resolution of 2.5° × 5° and a time resolution of 4 h. In this paper, we use those global plasmaspheric electron content product in 2008, 2010, 2011, 2012, 2014 for 1446 days to establish a global plasmaspheric empirical model based on empirical orthogonal function (EOF). The model can well characterize the spatiotemporal variation of plasmaspheric electron content (PEC) and the influence of solar radiation on it. Only the first four orders of EOF sequences can characterize the 98.43% features of the original PEC dataset. The principal component coefficient Pk is decomposed twice during modeling, and the combination of trigonometric function and linear function is used to model Pk to characterize the solar cycle, annual cycle, semi-annual cycle and quarter-cycle variation. We compare the PEC model values with the actual observation data, the results show that the empirical PEC model values are highly correlated with the actual observations. The correlation between the two is above 0.96, and the RMS maximum of the difference between the PEC model values and the observed values are 0.70 TECU, and the average of the difference between the PEC model values and the observed values are −0.18 TECU, respectively. In addition, we validate the reliability of the global plasmaspheric model established by two empirical orthogonal function decomposition method using actual observation data, according to the global distribution of the differences between the PEC model values and the observed values in low solar activity and high solar activity, it can be seen that under low solar activity and high solar activity conditions, the model has good adaptability.     

Estimation of differential code biases with Jason-2/3 onboard GPS observations

With the continuous deployment of Low Earth Orbit (LEO) satellites, the estimation of differential code biases (DCBs) based on GNSS observations from LEO has gained increasing attention. Previous studies on LEO-based DCB estimation are usually using the spherical symmetry ionosphere assumption (SSIA), in which a uniform electron density is assumed in a thick shell. In this study, we propose an approach (named the SHLEO method) to simultaneously estimate the satellite and LEO onboard receiver DCBs by modeling the distribution of the global plasmaspheric total electron content (PTEC) above the satellite orbit with a spherical harmonic (SH) function. Compared to the commonly used SSIA method, the SHLEO model improves the GPS satellite DCB estimation accuracy by 13.46% and the stability by 22.34%, respectively. Compared to the GPS satellite DCBs estimated based on the Jason-3-only observations, the accuracy and monthly stability of the satellite DCBs can be improved by 14.42% and 26.8% when both Jason-2 and Jason-3 onboard observations are jointly processed. Compared with the Jason-2 solutions, the GPS satellite DCB estimates based on the fusion of Jason-2 and Jason-3 observations have an improved consistency of better than 18.26% and 9.71% with the products provided by the Center for Orbit Determination in Europe (CODE) and Chinese Academy of Sciences (CAS). Taking the DCB products provided by the German Aerospace Center (DLR) as references, there is no improvement in accuracy of the GPS satellite DCB estimates based on the fusion of Jason-2 and Jason-3 observations than the Jason-2 solutions alone. A periodic variation is found in the time series of both the Jason-3 and Jason-2 onboard receiver DCB estimates. Preliminary analysis of the PTEC distribution based on the estimated SH coefficients are also presented.     

考虑跟踪性能约束的近空间飞行器控制器设计

针对一类可变后掠角的近空间飞行器(NSV)指令跟踪问题，考虑其受到外界扰动及参数不确定的影响，同时考虑系统的跟踪性能约束及姿态角速度约束，提出了一种基于模糊系统的切换控制器设计方法，确保系统在扰动影响以及给定约束下能够对给定信号进行稳定跟踪。建立了包含未知扰动及不确定项的近空间飞行器非线性切换模型；通过设计模糊系统对系统所受的总干扰进行实时估计，并基于反步法进行了切换控制器设计，在控制器中对干扰进行补偿；通过公共Barrier Lyapunov方法对系统稳定性及动态性能进行了分析。数值仿真算例校验了所提出方法的有效性及优越性。     

漂浮基空间机器人捕获卫星过程动力学模拟及捕获后混合体运动的RBF神经网络控制

 讨论了漂浮基空间机器人在轨捕获目标卫星过程的碰撞动力学建模,以及捕获操作结束后空间机器人与卫星混合体的稳定控制问题。首先采用多刚体动力学建模方法并结合空间机器人捕获目标卫星过程中的碰撞动力学特性,建立了漂浮基空间机器人在轨捕获漂浮卫星过程的动力学模型,并在此基础上计算出完成捕获操作后空间机器人与目标卫星混合体关节的运动速度。然后针对卫星及空间机器人系统惯性参数均是未知的复杂情况,应用上述模型、神经网络控制理论和Lyapunov稳定性理论,设计了空间机器人与卫星混合体在捕获过程碰撞冲击影响下稳定运动的高斯径向基函数神经网络控制方案,以达到对捕获卫星的有效控制。此外,高斯径向基函数神经网络控制方案具有不需要测量和反馈载体位置、移动速度与加速度的显著优点。系统数值仿真证实了上述控制方案的有效性。     

0-1线性规划的连续化求解方法

针对0-1线性规划的优化问题,提出一种惩罚函数方法。考虑到0-1线性规划的最优值特征,通过在目标函数中加上惩罚函数,将0-1离散线性规划模型连续化成非线性规划模型,并使用Matlab的Fmincon函数进行求解。经对多个算例的计算,并和其他算法比较,结果表明惩罚函数法的可行性和有效性。将该方法应用于实际的飞机排班问题上,取得比较满意的结果。     

一种求解Euler方程的简化气体动理学模型

基于一种简化的平衡分布函数,提出一种新的动理学通量矢量分裂（KineticFluxVectorSplitting,KFVS）格式来求解Euler方程。算例证实,新的KFVS格式继承了原始KFVS格式中捕捉激波准确、无数值震荡的优点。与原始的KFVS格式相比,由于简化的平衡分布函数与Maxwellian平衡分布函数相比复杂程度大大降低,计算中不需要再计算误差函数,新的KFVS格式的计算效率得到了很大提高。在同样的计算精度下,新KFVS格式比原始KFVS格式计算效率提高了一倍以上,显示了新格式具有良好的数值模拟性能。