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.     

Determination of attitude motion of the <Emphasis Type="Italic">Foton M-3</Emphasis> satellite according to the data of onboard measurements of the Earth’s magnetic field

The results of reconstruction of rotational motion of the Foton M-3 satellite during its uncontrolled flight in September 2007 are presented. The reconstruction was performed by processing the data of onboard measurements of the Earth’s magnetic field obtained by the DIMAC instruments. The measurements were carried out continuously throughout the flight, but the processing technique dealt with the data portions covering time intervals of a few orbital revolutions. The data obtained on each such interval were processed jointly by the least squares method with using integration of the equations of satellite motion relative to its center of mass. When processing, the initial conditions of motion and the used mathematical model’s parameters were estimated. The results of processing 16 data sets gave us complete information about the satellite motion. This motion, which began at a low angular velocity, had gradually accelerated and in five days became close to the regular Euler precession of an axisymmetric solid body. At the end of uncontrolled flight the angular velocity of the satellite relative to its lengthwise axis was 0.5 deg/s; the angular velocity projection onto the plane perpendicular to this axis had a magnitude of about 0.18 deg/s.     

Uncontrolled Attitude Motion of the Orbital Station MIR in the Last Months of Its Flight

The results of determination of the uncontrolled attitude motion of the orbital station MIR on four prolonged segments of its unmanned flight in 2000 and 2001 are presented. The determination was carried out on the basis of the data of onboard measurements of the Earth's magnetic field. The data obtained on a time interval of several hours were processed jointly by the least squares method by integration of the equations of motion of the station with respect to its center of mass. The processing resulted in the estimation of the initial conditions of the motion and of the parameters of the mathematical model used. Several types of regular motion were observed on sufficiently prolonged time intervals on the studied segments. Some of these motions were planned; others were established spontaneously.     

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.     

Uncontrollable Rotational Motion of the Mir Orbital Station

The results of determining the rotational motion of the Mir orbital station are presented for four long segments of its unmanned uncontrolled flight in 1999–2000. The determination was carried out using the data of onboard measurements of the Earth's magnetic field intensity. These data, taken for a time interval of several hours, were jointly processed by the least squares method with the help of integration of the equations of station motion relative to its center of mass. As a result of this processing, the initial conditions of motion and the parameters of the mathematical model used were evaluated. The technique of processing is verified using the telemetry data on angular velocity of the station and its attitude parameters. Two types of motion were applied on the investigated segments. One of them (three segments) presents a rotation around the axis of the minimum moment of inertia. This axis executes small oscillations with respect to a normal to the orbit plane. Such a motion was used for the first time on domestic manned orbital complexes. The second type of motion begins with a biaxial rotation which, in a few weeks, goes over into a motion very similar to the rotation around the normal to the orbit plane, but around the axis of the maximum moment of inertia.     

Determination of attitude motion of the Foton M-2 satellite using the data of onboard measurements of its angular velocity

We present the resutls of a prompt determination of the uncontrolled attitude motion of the Foton M-2 satellite, which was in orbit from May 31 to June 16, 2005. The data of onboard measurements of the angular velocity vector were used for this determination. The measurement sessions were carried out once a day, each lasting 83 min. Upon terminating a session, the data were transmitted to the ground to be processed using the least squares method and integrating the equations of motion of the satellite with respect to its center of mass. As a result of processing, the initial conditions of motion during a session were estimated, as well as parameters of the mathematical model used. The satellite’s actual motion is determined for 12 such sessions. The results obtained in flight completely described the satellite’s motion. This motion, having begun with a small angular velocity, gradually became faster, and in two days became close to the regular Euler precession of an axisymmetric solid body. On June 14, 2005 the angular velocity of the satellite with respect to its longitudinal axis was approximately 1.3 degrees per second, and the angular velocity projection onto a plane perpendicular to this axis had a magnitude of about 0.11 degrees per second. The results obtained are consistent with more precise results obtained later by processing the data on the Earth’s magnetic field measured on the same satellite, and they complement the latter in determination of the motion in the concluding segment of the flight, when no magnetic measurements were performed.     

Determination of the Microacceleration Quasisteady Component onboard the International Space Station

A comparison of two methods of determination of the microacceleration quasisteady component arising onboard the International Space Station was performed. In the first method the acceleration was calculated using the relative motion of the station reconstructed on the basis of telemetry data. The second method was a direct measurement of the microacceleration by a low-frequency accelerometer and a smoothing of the data obtained. The used measurements were made by the American accelerometer MAMS. The above comparison can theoretically be used to refine the position of the station center of mass relative to its body.     

Uncontrolled motion of the <Emphasis Type="Italic">Foton M-2</Emphasis> satellite and quasistatic microaccelerations on its board

The results of determining the uncontrolled rotational motion of the Foton M-2 satellite (in orbit from May 31 to June 16, 2005) are presented. The determination was made using the data of onboard measurements of the Earth’s magnetic field strength. Segments 270 min long (three orbits) were selected from these data covering the first two thirds of the flight. On each such segment the data were processed jointly by the least squares method using integration of the equations of attitude motion of the satellite. In processing, the initial conditions of motion and parameters of the used mathematical model were estimated. The thus obtained results gave a complete overview of the satellite motion. This motion, having started with a small angular velocity, gradually accelerated, and in two days became close to the regular Euler precession of an axisymmetric solid body. On June 09, 2005 (the last day of measurements) the angular velocity of the satellite relative to its lengthwise axis was about 1.1 deg/s, while the projection of the angular velocity onto a plane perpendicular to this axis had an absolute value of about 0.11 deg/s. Deviations of the lengthwise axis from a normal to the orbit plane did not exceed 60°. Based on the results of determination of the rotational motion of the satellite, calculations of quasi-static microaccelerations on its board are performed.     

Estimating the accuracy of the technique of reconstructing the rotational motion of a satellite based on the measurements of its angular velocity and the magnetic field of the Earth

The paper has studied the accuracy of the technique that allows the rotational motion of the Earth artificial satellites (AES) to be reconstructed based on the data of onboard measurements of angular velocity vectors and the strength of the Earth magnetic field (EMF). The technique is based on kinematic equations of the rotational motion of a rigid body. Both types of measurement data collected over some time interval have been processed jointly. The angular velocity measurements have been approximated using convenient formulas, which are substituted into the kinematic differential equations for the quaternion that specifies the transition from the body-fixed coordinate system of a satellite to the inertial coordinate system. Thus obtained equations represent a kinematic model of the rotational motion of a satellite. The solution of these equations, which approximate real motion, has been found by the least-square method from the condition of best fitting between the data of measurements of the EMF strength vector and its calculated values. The accuracy of the technique has been estimated by processing the data obtained from the board of the service module of the International Space Station (ISS). The reconstruction of station motion using the aforementioned technique has been compared with the telemetry data on the actual motion of the station. The technique has allowed us to reconstruct the station motion in the orbital orientation mode with a maximum error less than 0.6° and the turns with a maximal error of less than 1.2°.     

Determination of the Quasistatic Component of Microaccelerations at the MirStation

We describe the determination of the quasistatic component of microaccelerations, as it was done during the space experiments with the instruments DACON and ALICE-2. This component was calculated using the telemetric information related to the motion of the station with respect to its center of mass. The information consists of the values of the station angular velocity vector and the quaternion which specifies its orientation, determined at discrete instants of time. The quaternion is determined with a step of about 1 min, the angular velocity with a step of about 10 s. The information is used in the following manner. At first, the quaternion components corresponding to some time interval are smoothed by splines. Then, employing the obtained splines and kinematic equations, the angular velocity and acceleration of the station are calculated on this interval. Finally, the microacceleration is calculated as a function of time at the point of the location of the instrument. The data of measurements of the angular velocity are used for the purpose of control. As a rule, the available telemetric information allows one to find the quasi-static component for the entire time interval of carrying out the experiment. Examples of determination of this component for some experiments are presented. A comparison with the results of calculating the microacceleration quasistatic component by other methods is made.     

Uncontrolled Attitude Motion of the Foton-12 Satellite and Quasi-Steady Microaccelerations onboard It

The results of determination of the uncontrolled attitude motion of the Foton-12 satellite (placed in orbit on September 9, 1999, terminated its flight on September 24, 1999) are presented. The determination was carried out by the onboard measurement data of the Earth's magnetic field strength vector. Intervals with a duration of several hours were selected from data covering almost the entire flight. On each such interval the data were processed simultaneously using the least squares method by integrating the satellite's equations of motion with respect to the center of mass. The initial conditions of motion and the parameters of the mathematical model employed were estimated in processing. The results obtained provided for a complete representation of the satellite's motion during the flight. This motion, beginning with a small angular velocity, gradually sped up. The growth of the component of the angular velocity with respect to the longitudinal axis of the satellite was particularly strong. During the first several days of the flight this component increased virtually after every passage through the orbit's perigee. As the satellite's angular velocity increased, its motion became more and more similar to the regular Euler precession of an axisymmetric rigid body. In the last several days of flight the satellite's angular velocity with respect to its longitudinal axis was about 1 deg/s and the projection of the angular velocity onto the plane perpendicular to this axis had a magnitude of approximately 0.15 deg/s. The deviation of the longitudinal axis from the normal to the orbit plane did not exceed 60°. The knowledge of the attitude motion of the satellite allowed us to determine the quasi-steady microacceleration component onboard it at the locations of the technological and scientific equipment.     

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.     

Internal reference frames for representation and storage of visual information: the role of gravity.

Experimental studies of visual mechanisms suggests that the CNS represents image information with respect to preferred horizontal and vertical axes, as shown by a phenomenon known as the "oblique effect". In the current study we used this effect to evaluate the influence of gravity on the representation and storage of visual orientation information. Subjects performed a psychophysical task in which a visually-presented stimulus line was aligned with the remembered orientation of a reference stimulus line presented moments before. The experiments were made on 5 cosmonauts during orbital space flight and additionally on 13 subjects in conditions of normal gravity with a tilting chair. Data were analyzed with respect to response variability and timing. On earth, these measurements for this task show a distinct preference for horizontally and vertically oriented stimuli when the body and gravitational axes were aligned. This preference was markedly decreased or disappeared when the body axis was tilted with respect to gravity; this effect was not connected with ocular counter-rolling nor could we find a preference of any other intermediate axis between the gravity and body aligned axes. On the other hand, the preference for vertical and horizontal axes was maintained for tests performed in microgravity over the course of a 6 month flight, starting from flight day 6. We concluded that subjects normally process visual orientation information in a multi-modal reference frame that combines both proprioceptive and gravitational cues when both are available, but that a proprioceptive reference frame is sufficient for this task in the absence of gravity after a short period of adaptation. Some of the results from this study have been previously published in a preliminary report. Grant numbers: 99-04-48450.     

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.     

Determination of the inertia tensor of the <Emphasis Type="Italic">International Space Station</Emphasis> on the basis of telemetry data

We describe the method and results of determination of the inertia tensor of the International Space Station using telemetry data related to its attitude motion and the total angular momentum of gyrodines. A linear system of differential equations describing the variation of the total angular momentum of gyrodines on some time interval is derived on the basis of the data related to the station orientation in the same time interval. This linear system represents the theorem related to the variation of the total angular momentum of the station and gyrodines and takes into account the action of gravitational and aerodynamic moments upon the station. The solution to the system depends linearly on the components of the inertia tensor of the station and on the parameters specifying the aerodynamic moment. The estimates of these quantities are carried out by the least squares method on the condition of the best approximation by the solutions to the considered linear system of the telemetry values of the total angular momentum of the gyrodines.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 2, 2005, pp. 135–146.Original Russian Text Copyright © 2005 by Banit, Belyaev, Dobrinskaya, Efimov, Sazonov, Stazhkov.     

Low-Frequency Microaccelerations onboard the Foton-11 Satellite

We analyze the microacceleration measurements carried out onboard the Foton-11 satellite with the three-component accelerometer BETA. The microaccelerations were recorded virtually throughout the entire orbital flight of the Foton-11 satellite. The data obtained were analyzed in the following way. First they were used to determine the actual rotational motion of the satellite for several arbitrarily selected time intervals 4 h long. This problem was solved by constructing the approximation of the microacceleretation low-frequency component (previously determined from the data) by its calculated analog computed along the solutions to differential equations of rotational motion of the satellite. The approximation was made by the least squares method. As a result, those mathematical model parameters and the solutions to equations of motion were found that gave the best consistency of the microacceleretation low-frequency component and its calculated analog. Then the spectral analysis of the low-frequency component and its calculated analog was made. It was shown that, although basic harmonics of these functions coincided sufficiently well, some harmonics of the low-frequency component failed to be interpreted in terms of the satellite's rotational motion.     

Analysis of Low-Frequency Microaccelerations onboard the Foton-11Satellite

The method and the results of investigating the low-frequency component of microaccelerations onboard the Foton-11satellite are presented. The investigation was based on the processing of data of the angular velocity measurements made by the German system QSAM, as well as the data of measurements of microaccelerations performed by the QSAM system and by the French accelerometer BETA. The processing was carried out in the following manner. A low-frequency (frequencies less than 0.01 Hz) component was selected from the data of measurements of each component of the angular velocity vector or of the microacceleration, and an approximation was constructed of the obtained vector function by a similar function that was calculated along the solutions to the differential equations of motion of the satellite with respect to its center of mass. The construction was carried out by the least squares method. The initial conditions of the satellite motion, its aerodynamic parameters, and constant biases in the measurement data were used as fitting parameters. The time intervals on which the approximation was constructed were from one to five hours long. The processing of the measurements performed with three different instruments produced sufficiently close results. It turned out to be that the rotational motion of the satellite during nearly the entire flight was close to the regular Eulerian precession of the axially symmetric rigid body. The angular velocity of the satellite with respect to its longitudinal axis was about 1 deg/s, while the projection of the angular velocity onto the plane perpendicular to this axis had an absolute value of about 0.2 deg/s. The magnitude of the quasistatic component of microaccelerations in the locations of the accelerometers QSAM and BETA did not exceed 5 × 10^–5–10^–4m/s²for the considered motion of the satellite.     

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.     

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.     

相对测量信息不全的航天器相对轨道确定方法研究

针对近圆轨道航天器交会或远距离伴飞相对测量导航过程中,测量信息不全情况下的航天器自主相对轨道确定问题进行了研究.给出适合描述较远距离相对运动的二阶近似模型,并在采用雷达或光学测量的基础上设计了扩展卡尔曼滤波器.数学仿真结果表明,在观测量较少或存在部分区域不可测情况下,通过扩展卡尔曼滤波算法能够以较高精度估计出目标航天器的相对轨道.