一种改进阻尼倒数的空间机器人避奇异算法

针对自由漂浮空间机器人运动规划过程中存在的动力学奇异问题,提出一种动力学奇异转化加改进阻尼倒数的奇异回避算法。首先,通过运动学分解,间接求解机器人逆运动学方程,避免直接求解逆广义雅可比矩阵,以实现奇异参数的分离;然后,采用改进的高斯分布阻尼因子的阻尼倒数代替奇异参数的普通倒数,回避运动学奇异的影响。与现有阻尼倒数方法相比,提高了算法在奇异点附近的连续性,进一步减小了算法带来的跟踪误差。在PUMA类型空间机器人模型上进行仿真试验,校验算法的有效性。     

载体姿态无扰的自由漂浮空间机器人运动学特性研究

针对载体姿态无扰的自由漂浮空间机器人的运动学问题,引入具有非完整特性的载体 姿态无扰约束方程,推导了载体姿态无扰的自由漂浮空间机器人广义雅可比矩阵。采用等效 杆的概念,分析了普通状态的自由漂浮空间机器人和载体姿态无扰的自由漂浮空间机器人可 达工作空间;基于载体姿态无扰的自由漂浮空间机器人广义雅可比矩阵所得到的奇异方程, 研究了载体姿态无扰的自由漂浮空间机器人动力学奇异性。仿真结果验证了普通状态的自由 漂浮空间机器人与载体姿态无扰的自由漂浮空间机器人运动学问题的差异性。研究结果为载 体姿态无扰的自由漂浮空间机器人运动规划及控制研究奠定了理论基础。
     

Impulsive control for angular momentum management of tumbling spacecraft

We discuss an angular momentum control of a tumbling spacecraft. The proposed control method is to apply an impulse by a space robot arm, to measure and control the relative position and attitude between the target spacecraft, and then to apply another impulse until the rotational motion of the target spacecraft is well damped. A discrete controller is designed using the simplified equations of rotational motion through appropriate coordinate transformation. The stationary response under contact model uncertainty is investigated and stability condition is analytically derived. Numerical simulations are given to validate the proposed approach.     

一种新型并联微动机器人运动学分析

提出了一种新型的少自由度微动并联机器人,利用螺旋理论对于微动并联机器人的自由度进行分析,得到机构正常运动时所需要的驱动个数。利用矢量方法得到微动并联机器人的运动学反解。基于雅可比矩阵的可逆性,研究了微动机器人的奇异问题,从而确定该微动并联机器人的输入和输出之间的关系。     

自由飞行空间机器人的运动学建模研究

空间机器人将在未来的空间项目中发挥重要的作用。但在空间环境中，由于缺乏一个固定的基座，机械臂的运动会对机器人本体的位置和姿态产生扰动，其运动学将变得很复杂。为建立自由飞行空间机器人的运动学模型，通常是将系统的动量表达为某个预先选定连杆的速度的函数，进而得到具有广义雅可比矩阵形式的运动学模型。这种方法的缺点是系统的动量无法以通用表达式描述，并且对如何最优选取连杆也缺乏必要的说明和讨论。本文对此进行了研究，得到了自由飞行空间机器人系统动量的通用表达式；通过选取相应的连杆，得到了具有广义雅可比矩阵形式的运动学模型；也分析和讨论了连杆优选问题。     

空间机械臂关节故障引起速度突变的抑制

针对空间机械臂关节发生故障后的末端速度突变问题,提出一种基于退化雅克比矩阵的突变抑制方法。首先利用旋量理论建立空间机械臂的通用运动学模型,然后通过重构数学模型和调整控制模型与运行参数,引入补偿项,实现模型重构和在线调整中的平滑特性。并且开展了机械臂正常运行、故障运行和故障后抑制运行这三种工况下的控制仿真。结果表明,此方法对于实现单关节故障状态下机械臂末端速度突变抑制是有效的,对提高任务完成的连续性能力和降低故障影响具有重要的在轨应用意义。     

空间机器人动力学奇异回避的笛卡尔轨迹规划

针对空间机器人动力学奇异的回避问题,提出一种基于组合函数的笛卡尔轨迹规划方法。将梯形规划与正弦函数相结合,对机械臂末端的位姿进行参数化。机械臂沿着某一轨迹运动时,根据阻尼最小方差法(DLS)的特点,提出一种判断是否发生动力学奇异的方法,并据此进行轨迹规划。该方法可以保证空间机械臂运动过程中不会遇到动力学奇异。此外,将基座姿态扰动和机械臂运动时间作为目标函数的一部分,最终,轨迹规划问题转化为多目标优化问题,并利用混合整数规划的混沌粒子群优化算法(CPSO)进行求解。该优化算法能够改善标准粒子群算法(PSO)的“早熟”现象。仿真结果表明,新方法能够有效处理动力学奇异问题,减小基座姿态扰动及运动时间。     

轨道交会、悬停及绕飞控制的解析解方法

面向航天器交会对接、编队伴飞以及在轨操控等空间应用的需求,分别对近圆、椭圆轨道上航天器间的相对运动进行了分析与建模,在常值推力作用假设下进行了相对运动的解析求解。采用模型预测的方法获得航天器相对位置和相对速度的预期偏差。通过广义逆变换构造关于预期偏差的最小范数、最小二乘全状态反馈控制器。提出了一种普遍适用于近圆、椭圆轨道,可以实现轨道交会、相对悬停保持和循迹绕飞,对相对位置和相对速度进行同步控制的高精度、高稳定度相对制导律。仿真结果校验了方法的可行性和有效性。     

Modes of uncontrolled rotational motion of the <Emphasis Type="Italic">Progress M-29M</Emphasis> spacecraft

We have reconstructed the uncontrolled rotational motion of the Progress M-29M transport cargo spacecraft in the single-axis solar orientation mode (the so-called sunward spin) and in the mode of the gravitational orientation of a rotating satellite. The modes were implemented on April 3–7, 2016 as a part of preparation for experiments with the DAKON convection sensor onboard the Progress spacecraft. The reconstruction was performed by integral statistical techniques using the measurements of the spacecraft’s angular velocity and electric current from its solar arrays. The measurement data obtained in a certain time interval have been jointly processed using the least-squares method by integrating the equations of the spacecraft’s motion relative to the center of mass. As a result of processing, the initial conditions of motion and parameters of the mathematical model have been estimated. The motion in the sunward spin mode is the rotation of the spacecraft with an angular velocity of 2.2 deg/s about the normal to the plane of solar arrays; the normal is oriented toward the Sun or forms a small angle with this direction. The duration of the mode is several orbit passes. The reconstruction has been performed over time intervals of up to 1 h. As a result, the actual rotational motion of the spacecraft relative to the Earth–Sun direction was obtained. In the gravitational orientation mode, the spacecraft was rotated about its longitudinal axis with an angular velocity of 0.1–0.2 deg/s; the longitudinal axis executed small oscillated relative to the local vertical. The reconstruction of motion relative to the orbital coordinate system was performed in time intervals of up to 7 h using only the angularvelocity measurements. The measurements of the electric current from solar arrays were used for verification.     

Ballistics, navigation, and motion control for a spacecraft during its landing on the surface of Phobos

Basic concepts and algorithms laid as foundations of the scheme of landing on the Martian moon Phobos (developed for the Phobos-Grunt project) are presented. The conditions ensuring the landing are discussed. Algorithms of onboard navigation and control are described. The equations of spacecraft motion with respect to Phobos are considered, as well as their use for correction of the spacecraft motion. The algorithm of estimation of the spacecraft’s state vector using measurements with a laser altimeter and Doppler meter of velocity and distance is presented. A system for modeling the landing with a firmware complex including a prototype of the onboard computer is described.     

Gravitation effect on a flux of sporadic micrometeoroids in the vicinity of near-Earth orbits

The new approach to gravitation effect determination in calculating the flux of sporadic micrometeoroids in the near-Earth space is proposed. The technique is based on integration of the equations of motion of sporadic micrometeoroids with accounting for bending their trajectories when particles are approaching the Earth. The technique and results of calculation of the gravitational focusing factor k_g for various conditions are presented. The feature of the proposed technique for calculating coefficient k_g consists in the fact that this coefficient does not explicitly depend on the values of particles velocity at the last point. The results of investigation of coefficient k_g have shown that, for the given initial velocity of micrometeoroids, the values of this coefficient depend on deflection of its direction from the direction to the Earth center. It is shown that for low-altitude orbits the flux density can increase up to 60%. The distribution of probabilities of various directions of particles flying to spacecraft structural elements is found to be non-uniform.     

Singularity avoidance of control moment gyros by one-step ahead singularity index

A steering law of control moment gyros for spacecraft attitude control by using one-step ahead singularity index is addressed in this paper. In some recent studies, the null motion approaches or singularity robustness steering laws have been extensively investigated to avoid singular configurations for a control momentum gyro (CMG) cluster. As a novel approach different from them, the proposed approach is based on optimization technique by minimizing the one-step ahead singularity index. Modified approaches are also presented in this paper. The proposed one-step prediction method ultimately gives an optimized solution of gimbal rates with advanced ability to avoid a singularity. A singularity index for reliable computation of a gradient vector is also introduced. Finally, performance of the proposed algorithm is demonstrated by numerical simulations.     

基于参数辨识的冗余自由飞行空间机器人多臂协调运动规划

根据空间机器人姿态干扰特性，本文提出了利用多层间向神经网络来辨识自由飞行空间机器人本体质量的算法。其次，提出了基于姿态受限广义雅可比矩阵的多臂协调运动规划方法，较好地解决了多臂自由飞行空间机器人进行空间作业时的运动规划和姿态控制问题。     

机器人动态受力感知及零重力运动模拟技术

针对航天器在轨服务任务的地面零重力模拟需求,研究基于工业机器人的零重力运动模拟技术。建立基于BP神经网络的受力感知预测模型,该模型采用机器人末端姿态、加速度、角速度和角加速度作为输入层参数,采用机器人末端六维力传感器数据作为输出层参数,实现了对机器人末端负载的高精度动态受力感知。设计正交试验方法确定机器人的运动路径点进行样本数据采集,实现了受力感知预测模型对机器人全工作空间的覆盖。进一步,基于对机器人末端负载的受力感知数据,应用动力学理论计算负载在失重状态下的运动速度,并控制机器人执行相应的运动,实现了对机器人末端负载的零重力运动模拟。     

Designing corrections for the trajectory of the <Emphasis Type="Italic">Spektr-R</Emphasis> spacecraft in the event of immersions into the Moon’s sphere of influence

The results of updating the parameters of motion of the Spektr-R spacecraft at the end of 2016 have shown that, in January 2018, with a probability close to unity, the condition that a spacecraft stay in the Earth’s shadow is violated; however, in May of the same year, the ballistic life of the spacecraft will be terminated. Thus, in 2017, the question arose of how to design the correction of flight of this spacecraft using its onboard propulsion system. The correction was designed with allowance for the fact that, for the first time since it was launched, the spacecraft in the course of several years, beginning with 2017, repeatedly approaches the Moon, deeply immersing into its sphere of influence. This paper presents the technologically and organizationally convenient, allowable versions of upcoming correction of the Spektr-R spacecraft trajectory and justifies the particular scheme of its implementation.     

Rotational motion of <Emphasis Type="Italic">Foton M-4</Emphasis>

The actual controlled rotational motion of the Foton M-4 satellite is reconstructed for the mode of single-axis solar orientation. The reconstruction was carried out using data of onboard measurements of vectors of angular velocity and the strength of the Earth’s magnetic field. The reconstruction method is based on the reconstruction of the kinematic equations of the rotational motion of a solid body. According to the method, measurement data of both types collected at a certain time interval are processed together. Measurements of the angular velocity are interpolated by piecewise-linear functions, which are substituted in kinematic differential equations for a quaternion that defines the transition from the satellite instrument coordinate system to the inertial coordinate system. The obtained equations represent the kinematic model of the satellite rotational motion. A solution of these equations that approximates the actual motion is derived from the condition of the best (in the sense of the least squares method) match between the measurement data of the strength vector of the Earth’s magnetic field and its calculated values. The described method makes it possible to reconstruct the actual rotational satellite motion using one solution of kinematic equations over time intervals longer than 10 h. The found reconstructions have been used to calculate the residual microaccelerations.     

基于机动平台序列图像的空间飞行器定轨技术研究

针对机动观测平台单目光学成像系统的特点,在不能测定目标飞行器位置和速度的前提下,通过对成像系统与空间飞行器空间关系的分析,提出了视平均运动角速度与真平均运动角速度的概念,并构建了关于二者的约束方程,实现了基于测角数据的观测斜距的估计,从而解算出定轨所需的初始状态参数。基于观测斜距估计的轨道确定方法把对空间飞行器的定轨问题,归结为根据图像序列计算目标测角和根据测角数据确定观测斜距,解决了利用空间单目光学成像数据的定轨问题,并以高轨卫星为实例对定轨精度进行了仿真验证。     

Manifold dynamics in the Earth–Moon system via isomorphic mapping with application to spacecraft end-of-life strategies

Recently, manifold dynamics has assumed an increasing relevance for analysis and design of low-energy missions, both in the Earth–Moon system and in alternative multibody environments. With regard to lunar missions, exterior and interior transfers, based on the transit through the regions where the collinear libration points L₁ and L₂ are located, have been studied for a long time and some space missions have already taken advantage of the results of these studies. This paper is focused on the definition and use of a special isomorphic mapping for low-energy mission analysis. A convenient set of cylindrical coordinates is employed to describe the spacecraft dynamics (i.e. position and velocity), in the context of the circular restricted three-body problem, used to model the spacecraft motion in the Earth–Moon system. This isomorphic mapping of trajectories allows the identification and intuitive representation of periodic orbits and of the related invariant manifolds, which correspond to tubes that emanate from the curve associated with the periodic orbit. Heteroclinic connections, i.e. the trajectories that belong to both the stable and the unstable manifolds of two distinct periodic orbits, can be easily detected by means of this representation. This paper illustrates the use of isomorphic mapping for finding (a) periodic orbits, (b) heteroclinic connections between trajectories emanating from two Lyapunov orbits, the first at L₁, and the second at L₂, and (c) heteroclinic connections between trajectories emanating from the Lyapunov orbit at L₁ and from a particular unstable lunar orbit. Heteroclinic trajectories are asymptotic trajectories that travels at zero-propellant cost. In practical situations, a modest delta-v budget is required to perform transfers along the manifolds. This circumstance implies the possibility of performing complex missions, by combining different types of trajectory arcs belonging to the manifolds. This work studies also the possible application of manifold dynamics to defining suitable, convenient end-of-life strategies for spacecraft orbiting the Earth. Seven distinct options are identified, and lead to placing the spacecraft into the final disposal orbit, which is either (a) a lunar capture orbit, (b) a lunar impact trajectory, (c) a stable lunar periodic orbit, or (d) an outer orbit, never approaching the Earth or the Moon. Two remarkable properties that relate the velocity variations with the spacecraft energy are employed for the purpose of identifying the optimal locations, magnitudes, and directions of the velocity impulses needed to perform the seven transfer trajectories. The overall performance of each end-of-life strategy is evaluated in terms of time of flight and propellant budget.     

Nonlinear 6-DOF control of spacecraft docking with inter-satellite electromagnetic force

Compared to traditional docking systems, spacecraft docking with inter-satellite electromagnetic mechanism has distinct advantages. However, its 6-DOF control problem has not been adequately investigated. From our knowledge, this paper attempts to study the 6-DOF control problem for the first time. Based on the far-field electromagnetic force model and Hill's model, the dynamic model of translational motion is derived; using tracking control strategy, LQR method and estimate of Extended State Observer (ESO), an optimal and robust translational controller is designed to satisfy relative position/velocity requirements of soft docking. Representing the attitude of the docking spacecraft pair by unit quaternion, the attitude dynamic and kinematic models with quaternion expression are derived; using behavior-based coordinated control approach and ESO, a decentralized attitude controller is designed to simultaneously align one spacecraft with its absolute desired attitude and with the other spacecraft of the docking pair, requiring no angular velocity measurement and exhibiting better robust capability. The feasibility and performance of this proposed 6-DOF controller are validated by theoretical deduction and simulation results.     

Novel reaction control techniques for redundant space manipulators: Theory and simulated microgravity tests

This paper presents two novel redundancy resolution schemes aimed at locally minimizing the reaction torque transferred to the spacecraft during manipulator manoeuvres. The subject is of particular interest in space robotics because reduced reactions result in reduced energy consumption and longer operating life of the attitude control system. The first presented solution is based on a weighted Jacobian pseudoinverse and is derived by using Lagrangian multipliers. The weight matrix is defined by means of the inertia matrix which appears in the spacecraft reaction torque dynamics. The second one is based on a least squares formulation of the minimization problem. In this formulation the linearity of the forward kinematics and of the reaction torque dynamics equations with respect to the joint accelerations is used. A closed-form solution is derived for both the presented methods, and their equivalence is proven analytically. Moreover, the proposed solutions, which are suitable for real-time implementation, are extended in order to take into account the physical limits of the manipulator joints directly inside the solution algorithms. A software simulator has been developed in order to simulate the performance of the presented solutions for the selected test cases. The proposed solutions have then been experimentally tested using a 3D free-flying robot previously tested in an ESA parabolic flight campaign. In the test campaign the 3D robot has been converted in a 2D robot thanks to its modularity in order to perform planar tests, in which the microgravity environment can be simulated without time constraints. Air-bearings are used to sustain the links weight, and a dynamometer is used to measure the reaction torque. The experimental validation of the presented inverse kinematics solutions, with an insight on the effect of joint flexibility on their performance, has been carried out, and the experimental results confirmed the good performance of the proposed methods. In particular, two test cases have been analyzed in order to validate and evaluate the performance of both the unconstrained solution and the solution which takes into account the robot physical limits.