Determination of parameters of rotational motion of the <Emphasis Type="Italic">Monitor-E</Emphasis> spacecraft from the telemetry data on solar batteries’ current

The Monitor-E spacecraft executed uncontrolled flight due to emergency situation, no telemetry information on parameters of the spacecraft’s attitude motion being available. So, the problem arose to determine the spacecraft’s rotational motion from the accessible indirect information—the electric current provided by solar batteries. In this paper the integrated statistical technique is described, that allows one to solve this problem. The values of current, obtained over the time interval some tens of minutes long, have been processed simultaneously by the least squares method using the integration of the equations of spacecraft’s rotational motion. As a result of processing, the initial conditions of motion were estimated, and the spacecraft’s moments of inertia were updated, as well as the angles, specifying solar batteries position in the spacecraft-fixed coordinate system. The results of processing of 12 data sets are presented, which allowed us to reconstruct the actual rotational motion of the spacecraft.     

The Use of a Solar Electrojet Propulsion System for Jupiter Satellite Injection

The problem of optimization of the interplanetary trajectory of flight for a multistage spacecraft with high- and low-thrust engines into the Jupiter satellite orbit is considered. Low-thrust engines (stationary plasma engines) are used on a heliocentric flight segment. Their operation is maintained with electric power supply from solar batteries. The principal feasibility of the realization of such a project is shown, and estimations of the mass of a spacecraft placed into Jupiter's satellite orbit are presented.     

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.     

Realization of modes of satellite attitude motion with a small level of microaccelerations using electromechanical executive devices

Quasi-static microaccelerations are estimated for a satellite specially designed to perform space experiments in the field of microgravity. Three modes of attitude motion of the spacecraft are considered: passive gravitational orientation, orbital orientation, and semi-passive gravitational orientation. In these modes the lengthwise axis of the satellite is directed along the local vertical, while solar arrays lie in the orbit plane. The second and third modes are maintained using electromechanical executive devices: flywheel engines or gyrodynes. Estimations of residual microaccelerations are performed with the help of mathematical modeling of satellite’s attitude motion under the action of gravitational and aerodynamic moments, as well as the moment produced by the gyro system. It is demonstrated that all modes ensure rather low level of quasi-static microaccelerations on the satellite and provide for a fairly narrow region of variation for the vector of residual microacceleration. The semi-passive gravitational orientation ensures also a limited proper angular momentum of the gyro system.     

Experimental investigation of the modes of operation of uncontrolled attitude motion of the Progress spacecraft

Results of in-flight tests of three modes of uncontrolled attitude motion of the Progress spacecraft are described. These proposed modes of experiments related to microgravity are as follows: (1) triaxial gravitational orientation, (2) gravitational orientation of the rotating satellite, and (3) spin-up in the plane of the orbit around the axis of the maximum moment of inertia. The tests were carried out from May 24 to June 1, 2004 onboard the spacecraft Progress M1-11. The actual motion of this spacecraft with respect to its center of mass, in the above-mentioned modes, was determined by telemetric information about an electric current tapped off from solar batteries. The values of the current obtained during a time interval of several hours were processed jointly using the least squares method by integration of the equations of the spacecraft’s attitude motion. The processing resulted in estimation of the initial conditions of motion and of the parameters of mathematical models used. For the obtained motions the quasi-static component of microaccelerations was computed at a point onboard, where installation of experimental equipment is possible.     

Analysis of the orbital motion of the asteroid Apophis’ satellite

We have analyzed the orbital disturbed spacecraft motion near an asteroid. The equations of the asteroidocentric spacecraft motion have been used with regard to three perturbations from celestial bodies, the asteroid’s nonsphericity, and solar radiation pressure. It has been shown that the orbital parameters of the main spacecraft and a small satellite with a radio beacon can be selected such that the orbits are rather stable for a fairly long period of time, i.e., a few weeks for the main spacecraft with an orbit initial radius of ~0.5 km and a few years before approaching Apophis with the Earth in 2029, for a small satellite at an orbit initial radius of ~1.5 km. The initial orientation of the spacecraft orbital plane perpendicular to the sunward direction is optimal from the point of view of the stability of the spacecraft flight near an asteroid.     

Optimal control of the deployment process of solar wings on spacecraft

The optimal attitude control problem of spacecraft during the stretching process of solar wings is investigated in this paper. The dynamical equations of the nonholonomic system are derived from the conservation principle of the angular momentum of the multibody system. Attitude control of the spacecraft with internal motion is reduced to a nonholonomic motion planning problem. The spacecraft attitude control is transformed into the steering problem for a drift free control system. The optimal solution for steering a spacecraft with solar wings is presented. The controlled motion of spacecraft is simulated for two cases. The numerical results demonstrate the effectiveness of the optimal control approach.     

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.     

Reconstruction of spacecraft rotational motion using a Kalman filter

Quasi-static microaccelerations of four satellites of the Foton series (nos. 11, 12, M-2, M-3) were monitored as follows. First, according to measurements of onboard sensors obtained in a certain time interval, spacecraft rotational motion was reconstructed in this interval. Then, along the found motion, microacceleration at a given onboard point was calculated according to the known formula as a function of time. The motion was reconstructed by the least squares method using the solutions to the equations of satellite rotational motion. The time intervals in which these equations make reconstruction possible were from one to five orbital revolutions. This length is increased with the modulus of the satellite angular velocity. To get an idea on microaccelerations and satellite motion during an entire flight, the motion was reconstructed in several tens of such intervals. This paper proposes a method for motion reconstruction suitable for an interval of arbitrary length. The method is based on the Kalman filter. We preliminary describe a new version of the method for reconstructing uncontrolled satellite rotational motion from magnetic measurements using the least squares method, which is essentially used to construct the Kalman filter. The results of comparison of both methods are presented using the data obtained on a flight of the Foton M-3.     

Uncontrolled rotational motion of the <Emphasis Type="Italic">Aist</Emphasis> small spacecraft prototype

The results of reconstructing the uncontrolled rotational motion of the Aist small spacecraft prototype during its flight in early 2014 have been presented. The reconstruction was carried out by processing data from onboard measurements of the Earth’s magnetic field. The processing procedure used portions of data covering intervals of time with durations ranging from a few dozen minutes to three hours. Data obtained in each such interval were processed jointly by the least-squares method by integrating the equations of the satellite motion relative to the center of mass. The initial conditions of the motion and the parameters of the used mathematical model during processing have been estimated. The results of processing for several data intervals have provided a fairly complete picture of the satellite motion. This was the weakly disturbed Euler–Poinsot motion.     

共面快速受控绕飞轨迹设计与控制

绕飞运动在航天器在轨服务与在轨支援、辅助航天员舱外活动、航天器编队飞行、空间交会对接等空间活动中具有重要应用。分析了快速受控绕飞的可行性和主要过程，建立了适用于目标航天器运行在圆轨道上的共面快速绕飞和进入绕飞与退出绕飞的轨迹设计模型，采用多速度脉冲控制方法和等角度，等时间控制方式对绕飞轨迹进行控制。仿真计算结果表明所提出的快速受控绕飞轨迹设计模型和控制方法可以实现对圆轨道目标航天器的共面快速受控绕飞。     

The use of the Earth’s gravitational potential for interplanetary flights

An efficient scheme of the use of the Earth’s gravity in interplanetary flights is suggested, which opens up new opportunities for exploration of the solar system. The scheme of the gravitational maneuver allows one to considerably reduce the spacecraft mass consumption for a flight and the time of flight. An algorithm of the gravitational maneuver is suggested that takes into account the restriction on the altitude of a planet flyby. Estimates are made of transport capabilities for delivery of a spacecraft to the orbits of Jupiter, Saturn, and Uranus. The spacecraft is based on a middle-class carrier launcher of the Soyuz type and includes chemical and electric plasma jet engines of the SPD-140 type, which use solar energy.     

Trajectory correction of the Spektr-R spacecraft motion

The results of refining the parameters of the Spektr-R spacecraft (RadioAstron project) motion after it was launched into the orbit of the Earth’s artificial satellite in July 2011 showed that, at the beginning of 2013, the condition of staying in the Earth’s shadow was violated. The duration of shading of the spacecraft exceeds the acceptable value (about 2 h). At the end of 2013 to the beginning of 2014, the ballistic lifetime of the spacecraft completed. Therefore, the question arose of how to correct the trajectory of the motion of the Spektr-R satellite using its onboard propulsion system. In this paper, the ballistic parameters that define the operation of onboard propulsion system when implementing the correction, and the ballistic characteristics of the orbital spacecraft motion before and after correction are presented.     

Effect of Microaccelerations on the Distribution of a Dopant in a Semiconductor Melt during Space Flight

The effect of residual microaccelerations on the distribution of a dopant in a semiconductor melt located in a heated closed cavity onboard an Earth-orbiting satellite is considered in the context of a model problem of thermal convection. The amplitude–frequency characteristics of the response of this distribution to the perturbing microaccelerations are obtained. It is demonstrated that the effect of low-frequency microaccelerations decreases when the frequency increases. A comparison is made of the macroscopic inhomogeneities of the dopant concentration due to the actual low-frequency (quasi-static) component of microaccelerations onboard different spacecraft: the orbital station Mir, the satellite Foton-11, a Space Shuttle orbiter, and the International Space Station. A substantial effect of the rotational motion of the spacecraft on the character of the time behavior of a macroscopic inhomogeneity is demonstrated.     

USA-193 decay predictions using public domain trajectory data and assessment of the post-intercept orbital debris cloud

In early 2008, the need arose to predict the orbital decay of the American spacecraft USA-193, whose characteristics, function and orbit were classified information. With no orbit data and independent Italian tracking capability available, we turned our attention on the orbits determined by a worldwide network of about 20 visual satellite observers. The orbits of USA-193 obtained from their visual observations were therefore used as the sole source of orbit information. Contrary to our expectations, this exercise was extremely successful and we learned a lot in the process. The orbits provided by the visual observers were very accurate for such a low satellite (although the minimum and very stable level of solar activity helped considerably); however, data gaps of a few days were sometimes possible, due to unfavorable pass geometry or weather and light conditions. In any case, the orbital period and the semimajor axis were so accurate that it was possible for us to obtain very good decay fits using special perturbation software, including various atmospheric density models together with all the other relevant perturbing accelerations. We were therefore able to estimate accurate values of the ballistic parameter and the resulting decay and reentry predictions were extremely stable. Amateur optical observations and images of USA-193 had also led to a rough estimation of the shape and sizes of the satellite, revealing that the solar arrays had never been deployed. With this information, and taking into account our estimates of the ballistic parameter, we obtained reasonable and consistent values of the spacecraft mass. Based on previous reentry fragmentation analyses, we were then able to guess the expected USA-193 casualty area, casualty expectancy, debris ground footprint and probability of impact in Italy. Lastly, after the decision by the US Government to destroy the satellite, we independently predicted the interception time windows and the post-event ground tracks. Following the successful spacecraft breakup, we analyzed the evolution of the resulting debris cloud and assessed its (very limited) adverse impact on the circumterrestrial environment.     

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.     

一种精确计算航天器本体对太阳电池阵遮挡的方法

针对航天器太阳电池阵电设计和热设计中需要准确考虑阴影影响的问题,提出了一种可以精确计算太阳电池阵阴影的方法。首先使用三角网格来建立3D模型,其次考虑了航天器本体构件间的相对运动,然后用一个“包围盒”去截取模型上的三角网格点,再把这些点投射到太阳电池阵上,最后把太阳电池阵分成小方格,分别使用逐点比较法和最小矩形法来生成阴影图形。给出了月球车遮挡计算的实例,仿真分析表明生成阴影图形时,最小矩形法具有更高的速度。使用本文所述的方法,可计算结构复杂、构件间有相对运动的航天器本体对太阳电池阵造成的遮挡,能生成精确的阴影图形,为后续计算受遮挡的太阳电池阵的输出特性奠定了基础。     

Variation of the Satellite Orbit Altitude by the Light Pressure Force

The possibility of the uncontrolled increase of the altitude of an almost circular satellite orbit by the force of the light pressure is investigated. The satellite is equipped with a damper and a system of mirrors (solar batteries can serve as such a system). The flight of the satellite takes place in the mode of a single-axis gravitational orientation, the axis of its minimum principal central moment of inertia makes a small angle with the local vertical and the motion of the satellite around this axis constitutes forced oscillations under the impact of the moment of force of the light pressure. The form of the oscillations and the initial orbit are chosen so that the transverse component of the force of the light pressure acting upon the satellite be positive and the semimajor axis of the orbit would continuously increase. As this takes place, the orbit remains almost circular. We investigate the evolution of the orbit over an extended time interval by the method which employs separate integration of the equations of the orbital and rotational motions of the satellite. The method includes outer and inner cycles. The outer cycle involves the numerical integration of the averaged equations of motion of the satellite center of mass. The inner cycle serves to calculate the right-hand sides of these equations. It amounts to constructing an asymptotically stable periodic motion of the satellite in the mode of a single-axis gravitational orientation for current values of the orbit elements and to averaging the equations of the orbital motion along it. It is demonstrated that the monotone increase of the semimajor axis takes place during the first 15 years of motion. In actuality, the semimajor axis oscillates with a period of about 60 years. The eccentricity and inclination of the orbit remain close to their initial values.     

基于脉冲星的空间飞行器自主导航技术研究

研究了利用X射线脉冲星实现空间飞行器自主导航的原理,指出通过测量脉冲星发射的脉冲到达航天器与到达太阳系质心的时间差可以确定该航天器的位置。给出了两种适合用于处理脉冲星观测量的方法,即基于最小二乘(LS)的方法和基于卡尔曼滤波(KF)的方法。通过数值仿真验证了X射线脉冲星自主导航技术的可行性,并且说明采用KF算法能够有效利用卫星轨道运动方程中包含的导航信息,从而获得更高的导航精度。     

Study of microaccelerations onboard the <Emphasis Type="Italic">International Space Station</Emphasis> with the DAKON-M convection sensor

The results of experiments with the DAKON-M convection sensor onboard the Russian orbital segment of the International Space Station are described. A comparison of the sensor measurements with the results of calculation of the quasistatic microacceleration component at the point of installation is made. For this comparison we have used three measurement intervals of the experiments in 2009, during which spacecraft were docked with the station, undocked from it, and actuation of jet engines of the attitude control system took place. When calculating microacceleration, we use the measurement data of the low-frequency MAMS accelerometer, installed on the American segment, and the telemetry data on the ISS rotational motion. This information allowed one to convert the MAMS measurements to the point of installation of the DAKON-M convection sensor. A comparison of sensor measurements with calculated microaccelerations showed sufficiently accurate coincidence between the calculated and measured data.