无反作用力矩空间机器人轨迹跟踪控制

针对目前关节驱动方式的不足,提出一种采用无反作用力矩驱动方式的空间机器人设计概念。系统平台与各节机械臂上均安装一组金字塔构型的控制力矩陀螺（CMGs）作为力矩执行机构,各臂间由自由球铰连接。采用Rodriguez参数描述平台姿态和机械臂角位移,利用Kane方程建立系统动力学模型。在此基础上,设计了渐近稳定的轨迹跟踪控制律,使得平台姿态和机械臂角位移跟踪期望运动轨迹;并设计了CMGs操纵律,使之准确输出期望控制力矩。此外,研究了机械臂工作空间到关节空间的轨迹规划算法,使得所设计的控制律也可应用于工作空间的轨迹跟踪控制。由三关节系统的仿真结果,验证了无反作用力矩设计概念的可行性和轨迹规划算法的有效性。     

基于HLA的航天飞行任务联合仿真系统设计

简要说明了HLA的基本特点,描述了一个已经实现的航天测控仿真系统,并将其结构思想与HLA进行了对照。在此基础上,提出了一个依据HLA标准,对航天员座舱仿真、火箭和飞船仿真、测控网仿真等进行集成的航天飞行任务联合仿真系统框架,分析了各个仿真联邦成员的功能需求和信息接口,并着重介绍了针对任务要求的时间管理、数据通信等关键RTI机制的设计和实现方法,以及参照FEDEP模型和借助商业成品软件工具的系统开发策略。     

基于小波变换的单脉冲雷达目标RCS测量方法研究

随着现代信号处理技术的发展,对非平稳信号分析和处理的小波分析技术已成功应用于雷达目标特性分析领域,大功率单脉冲雷达作为我国航天测控网的主干设备,具有一定的目标特性识别能力。本文主要针对脉冲雷达RCS测量原理,讨论了基于小波变换的单脉冲雷达空间目标RCS测量方法,提出应发挥窄带低分辨率雷达在未来空间目标识别中的作用。     

空间碎片材料检测技术研究

为了对空间碎片材料的成分进行快速检测,分析了空间碎片的来源和构成,并对空间碎片主要成分的国内外检测技术手段进行了详细对比,采用了一种基于LIBS的网络化空间碎片检测方案。该方案具有分析速度快、测量精度和灵敏度高、可多元素同时测量、空间适应性好且无需样品预处理等特点,同时可实现检测流程全自动化和多线程化。     

空间太阳能电站微波能量传输验证方案设计

微波能量传输技术作为空间太阳能电站（SpaceSolarPowerStation,SSPS）的关键技术之一,目前的研究和验证工作均集中在各单项技术的突破和验证,缺乏针对SSPS系统特点的全面优化设计。文章根据SSPS的工作模式给出了全面验证空间太阳能电站微波能量传输的验证系统方案设计,对收发天线进行了一体化设计,利用了幅度近似高斯分布的发射阵列场分布设计和低反射的接收整流阵列设计,以高精度来波方向测量和高精度移相控制为波束指向控制的技术途径。对验证系统的波束收集效率进行了分析,收集效率可达94.2%,比传统均匀分布系统高出17.6%。验证系统可从系统规模缩比、波束扫描范围、发射天线口径场分布、整流天线处功率密度、反向波束控制方法等方面模拟SSPS微波能量传输工作模式,推动SSPS系统技术的发展。     

微波能量传输系统设计及试验

面向空间太阳能电站应用,进行了固态体制微波能量传输技术研究。针对能量传输波束扩散导致收集效率低的问题,研究了基于人工媒质理论设计的完美匹配层的能量接收整流表面,通过调节人工媒质单元的结构参数实现天线输出阻抗与整流电路输入阻抗的共轭匹配,同时抑制整流电路高次谐波,省去原有匹配及低通滤波器,简化电路结构、实现高效微波能量吸收与转换。以空间太阳能电站规定的波束中心传输微波功率密度限制作为能量接收整流表面设计的约束条件,设计能量接收整流表面,结合固态体制微波能量发射端,搭建5.8GHz小规模微波能量传输系统开展了地面试验验证,实测结果显示整流表面能量转换效率最高为57.7%。此次试验验证了先进的固态能量传输试验系统,为空间太阳能地面缩比试验及未来空间太阳能电站的建设提供了技术支撑。     

Yarkovsky-Schach effect on space debris motion

The Yarkovsky-Schach effect is a small perturbation affecting Earth satellites and space debris illuminated by the Sun. It was first applied to the orbit of LAGEOS satellites as an explanation of the residuals in orbital elements. In this work, we carry out several numerical integration tests taking into consideration various orbit and rotation parameters, in order to analyse this effect in a broader context. The semi-major axis variations remain small and depend on the spin axis attitude with respect to the Sun. We show that the force amplitude is maximised for orbits inclined with $i\approx$?20–30°. We also observe the influence on other orbital elements, notably on the orbit inclination. However, these effects are clearly observed only on long timescales; in our simulations, we propagated the orbits for 200?y. The Yarkovsky-Schach effect is thus confirmed to have a minuscule magnitude. It should be taken into account in studies requiring high-precision orbit determination, or on expanded timescales.     

Optimal Earth’s reentry disposal of the Galileo constellation

Nowadays there is international consensus that space activities must be managed to minimize debris generation and risk. The paper presents a method for the end-of-life (EoL) disposal of spacecraft in Medium Earth Orbit (MEO). The problem is formulated as a multiobjective optimisation one, which is solved with an evolutionary algorithm. An impulsive manoeuvre is optimised to reenter the spacecraft in Earth’s atmosphere within 100?years. Pareto optimal solutions are obtained using the manoeuvre $\mathrm{\Delta }v$ and the time-to-reentry as objective functions to be minimised. To explore at the best the search space a semi-analytical orbit propagator, which can propagate an orbit for 100?years in few seconds, is adopted. An in-depth analysis of the results is carried out to understand the conditions leading to a fast reentry with minimum propellant. For this aim a new way of representing the disposal solutions is introduced. With a single 2D plot we are able to fully describe the time evolution of all the relevant orbital parameters as well as identify the conditions that enables the eccentricity build-up. The EoL disposal of the Galileo constellation is used as test case.     

Simulation-based evaluation of a cold atom interferometry gradiometer concept for gravity field recovery

The prospects of future satellite gravimetry missions to sustain a continuous and improved observation of the gravitational field have stimulated studies of new concepts of space inertial sensors with potentially improved precision and stability. This is in particular the case for cold-atom interferometry (CAI) gradiometry which is the object of this paper. The performance of a specific CAI gradiometer design is studied here in terms of quality of the recovered gravity field through a closed-loop numerical simulation of the measurement and processing workflow. First we show that mapping the time-variable field on a monthly basis would require a noise level below $5\text{mE}/\sqrt{\text{Hz}}$. The mission scenarios are therefore focused on the static field, like GOCE. Second, the stringent requirement on the angular velocity of a one-arm gradiometer, which must not exceed ${10}^{-6}$?rad/s, leads to two possible modes of operation of the CAI gradiometer: the nadir and the quasi-inertial mode. In the nadir mode, which corresponds to the usual Earth-pointing satellite attitude, only the gradient $V_{\mathit{yy}}$, along the cross-track direction, is measured. In the quasi-inertial mode, the satellite attitude is approximately constant in the inertial reference frame and the 3 diagonal gradients $V_{\mathit{xx}}\text{,}{V}_{\mathit{yy}}$ and $V_{\mathit{zz}}$ are measured. Both modes are successively simulated for a 239?km altitude orbit and the error on the recovered gravity models eventually compared to GOCE solutions. We conclude that for the specific CAI gradiometer design assumed in this paper, only the quasi-inertial mode scenario would be able to significantly outperform GOCE results at the cost of technically challenging requirements on the orbit and attitude control.     

An estimation of Envisat’s rotational state accounting for the precession of its rotational axis caused by gravity-gradient torque

The rotational state of Envisat is re-estimated using the specular glint times in optical observation data obtained from 2013 to 2015. The model is simplified to a uniaxial symmetric model with the first order variation of its angular momentum subject to a gravity-gradient torque causing precession around the normal of the orbital plane. The sense of Envisat’s rotation can be derived from observational data, and is found to be opposite to the sense of its orbital motion. The rotational period is estimated to be $(120.674\pm 0.068)·\exp \left((4.5095\pm 0.0096)\times {10}^{-4}·t\right)\text{s}$, where t is measured in days from the beginning of 2013. The standard deviation is 0.760?s, making this the best fit obtained for Envisat in the literature to date. The results demonstrate that the angle between the angular momentum vector and the negative normal of the orbital plane librates around a mean value of $8.53°\pm 0.42°$ with an amplitude from about $0.7°$ (in 2013) to $0.5°$ (in 2015), with the libration period equal to the precession period of the angular momentum, from about 4.8?days (in 2013) to 3.4?days (in 2015). The ratio of the minimum to maximum principal moments of inertia is estimated to be $0.0818\pm 0.0011$, and the initial longitude of the angular momentum in the orbital coordinate system is $40.5°\pm 9.3°$. The direction of the rotation axis derived from our results at September 23, 2013, UTC 20:57 is similar to the results obtained from satellite laser ranging data but about $20°$ closer to the negative normal of the orbital plane.