环月轨道器运动中的月球物理天平动问题

采用类似对地球岁差章动的处理方法,讨论月球物理天平动对月心赤道坐标系 以及月球卫星轨道的影响。在对月球物理天平动的分析表达式与高精度数值历表进行比较的 基础上,给出了相应的月球引力位的变化及相应的坐标系附加摄动解,清楚地表明了对月球 卫星轨道影响的规律。所获结果与数值解进行了比对,证实了从定性和定量两个方面来看本 文的讨论都是有意义的,从而表明月球物理天平动分析解的简单表达式在某些问题中（特别 是定轨和预报）是有实用价值的,而且在建立轨道摄动分析解时,无需像地球卫星那样,去 引进混合形式的轨道坐标系,采用历元月心平赤道坐标系即可。最终表明：无论是采用数值 法定轨和预报还是分析法定轨和预报,均可采用统一的月心平赤道坐标系,这可避免一些不 必要的坐标转换。     

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.     

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°.     

月球卫星的自主轨道确定

分别利用紫外三轴姿态敏感器测出的月心方向和测距仪测出的月心距作为观测量,采用推广卡尔曼滤波方法实时确定出月球卫星环月期间的轨道,考虑到月球卫星在环月期间太阳对月球光照条件的限制,提出分段虑波策略,并给出了仿真分析.     

Determining Spacecraft Motion from Four Star Sensor Measurements

The BOKZ-M60 star sensor (a module that measures the coordinates of stars) has been designed for determining the parameters of the orientation of the intrinsic coordinate system relative to the inertial system from observations of stellar sky sections. The methods and results of processing of measurements by a set of four BOKZ-M60 sensors on the Resurs-P satellite no. 2 have been described. The time interval at which the satellite was in orbital orientation exceeds three orbital revolutions (19003 s). The joint processing of measurements by the four sensors conducted at the same time instants allowed the sensors to be associated with the universal coordinate system. With a root-mean-square error of less than 0.4′′ for each angle of rotation around its axes, this system is consistent with the model of the satellite’s rotational motion. The position of the universal system with respect to the instrumental coordinate system of the satellite was determined from the angular velocity measurements. Here, the root-mean-square errors for the values determined by the angles of rotation of the universal system around its axes were 0.044°, 0.051°, and 0.18°. The low-frequency (with frequencies less than 0.05 Hz) variations in the positions of intrinsic sensor coordinate systems relative to the universal system do not exceed 10′′. These are periodic variations with a fundamental frequency equal to the orbital frequency. The root-mean-square values of high-frequency components of these variations do not exceed 18′′.     

Estimate of accuracy of determining the orientation of the star sensor system according to the experimental data

The BOKZ-M60 star sensor (Unit for Measuring Star Coordinates) is intended for determining the parameters of the orientation of the axes of the intrinsic coordinate system relative to the axes of the inertial system by observations of the regions of the stellar sky. It is convenient to characterize an error of the single determination of the orientation of the intrinsic coordinate system of the sensor by the vector of an infinitesimal turn of this system relative to its found position. Full-scale ground-based tests have shown that, for a resting sensor the root-mean-square values of the components of this vector along the axes of the intrinsic coordinate system lying in the plane of the sensor CCD matrix are less than 2″ and the component along the axis perpendicular to the matrix plane is characterized by the root-mean-square value of 15″. The joint processing of one-stage readings of several sensors installed on the same platform allows us to improve the indicated accuracy characteristics. In this paper, estimates of the accuracy of systems from BOKZ-M60 with two and four sensors performed from measurements carried out during the normal operation of these sensors on the Resurs-P satellite are given. Processing the measurements of the sensor system allowed us to increase the accuracy of determining the each of their orientations and to study random and systematic errors in these measurements.     

Determination of the Attitude Motion of a Satellite Based on a Kinematical Model of Motion and Using Data of Measurements of Its Angular Velocity and the Earth’s Magnetic Field Strength

An integral statistical procedure of determination of the attitude motion of a satellite using the data of onboard measurements of angular velocity vectors and the strength of the Earth’s magnetic field (EMF) is suggested. The procedure uses only the equations of kinematics of a solid body and is applicable to determining both controlled and uncontrollable motions of a satellite at any external mechanical moments acting upon it. When applying this procedure, the data of measurements of both types, accumulated during a certain interval of time, are processed jointly. The data of measuring the angular velocity are smoothed by discrete Fourier series, and these series are substituted into kinematical Poisson equations for elements of the matrix of transition from a satellite-fixed coordinate system to the orbital coordinate system. The equations thus obtained represent a kinematical model of the satellite motion. The solution to these equations (which approximate the actual motion of a satellite) is found from the condition of the best (in the sense of the least squares method) fit of the data of measuring the EMF strength vector to its calculated values. The results of testing the suggested procedure using the data of measurements of the angular velocity vectors onboard the Foton-12 satellite and measurements of EMF strengths are presented.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 4, 2005, pp. 295–305.Original Russian Text Copyright © 2005 by Abrashkin, Volkov, Voronov, Egorov, Kazakova, Pankratov, Sazonov, Semkin.     

月球坑内空间环境热流的分布研究

文章推导获得了月面阳光矢量和月球坑内探测器表面关系的表达式,采用蒙特卡罗法数值计算探测器各表面的辐射热流,分析了月球坑内环境辐射热流的分布特征。讨论了月面纬度、月球坑尺寸、月面发射率等参数对探测器表面热流的影响。计算表明:月面纬度增加及月球坑坑口半径的减小均会导致某些时刻阳光无法照射到坑底,从而与其他工况有较大的计算差别,两参数对探测器各面热流的影响规律不同;月面发射率对探测器表面热流的影响作用较大,增加月面发射率能明显降低探测器某些表面的辐射热流。     

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.     

An Analysis of the Low-Frequency Component in Measurements of Angular Velocity and Microacceleration onboard the Foton-12 Satellite

We compare the results of two methods used to determine the angular velocity of the Foton-12 satellite and the low-frequency component of microaccelerations onboard it. The first method is based on reconstruction of the satellite's rotational motion using the data of onboard measurements of the strength of the Earth's magnetic field. The motion (time dependence of the orientation parameters and angular velocity) was found from the condition of best approximation of the measurement data by the functions calculated along the solutions to equations of attitude motion of the satellite. The solutions found were used to calculate the quasistatic component of microaccelerations at certain points of the satellite, in particular, at the point of location of an accelerometer of the QSAM system. Filtration of the low-frequency component of the angular velocity and microacceleration from the data of measurements by a sensor of angular velocity and by the accelerometer of this system served as a second method. The filtration was made using the discrete Fourier series. A spectral analysis of the functions representing the results of determining the angular velocity and microacceleration by both methods is performed. Comparing the frequencies and amplitudes of the harmonic component of these functions allowed us to estimate the accuracy of measurements made by the QSAM system in the low-frequency range.     

Further analyses of the slide lander and of drop delivery systems for improved lunar surface access

The paper discusses two methods for lunar surface access. One method is characterized by very little propellant consumption for landing (lunar slide lander, LSL); the other (drop delivery method, DDM) avoids release of propellant gases at the lunar surface. The LSL is of particular promise and importance. Its analysis and development introduces a new field of cosmic dynamics—harenodynamics, the science and technology of interaction of surfaces with dust and sand at near-orbital speeds down to low velocities. Although the data base is small so far and needs to be enlarged, analyses of the LSL so far indicate promise as to feasibility. Its realization will revolutionize present conventional concepts of lunar development. The DDM offers cost-effective alternatives to conventional lunar landing by retrothrust, but on a more selective basis, because propellant-savings are secondary to avoiding gas releases into the lunar high-vacuum environment. It was found that stationary energy absorption systems (EAS) are required, because vehicle-attached absorption systems are entirely inadequate. This requires large structures which can be built on the lunar surface only after a somewhat higher degree of industrial capability has been established. However, they can be built entirely of lunar materials. Altogether, the LSL is the more significant of the two both in terms of economy and of operational scope.     

The Mode of Gravitational Orientation for the International Space Station

The results of the determination of the uncontrolled attitude motion of the International Space Station during its unmanned flight in 1999 are presented. The data of onboard measurements of three components of the angular velocity are used for this determination. These data covering an interval of slightly less than one orbit were jointly processed by the least squares method, by integrating the equations of motion of the station relative to its center of mass. As a result of this processing, the initial conditions of the motion and the parameters of the mathematical model used were estimated. The actual motion of the station has been determined for 20 such intervals during April–November. Throughout these intervals, the station rotated about the axis of the minimum moment of inertia, the latter executing small oscillations relative to the local vertical. Such a mode, known as the mode of gravitational orientation of a rotating satellite or the mode of generalized gravitational orientation, was planned before the flight. The measurements were made to verify it. The quasistatic component of the microaccelerations aboard the station is estimated for this mode.     

Functional restrictions on the orientation of onboard and ground methods in the RadioAstron project

In this paper, functional restrictions on the orientation of the onboard systems of the space radio telescope, stations of scientific data receiving and ground radio telescopes in the RadioAstron project are considered. Restrictions important for practical problems of scheduling observations with ground-space radio interferometer are discussed in detail. An algorithm for calculating the angles of the narrow-beam antenna drive, which takes into account technological restrictions on the capability of its spatial orientation, is presented.     

北斗系统地固系自主定轨算法

为避免传统惯性系自主定轨算法中的多次坐标转换，针对地固系分布式自主定轨算法原理及其中的关键技术开展研究。给出在地固系中进行自主星历生成的总体框图，推导了地固系自主定轨算法的基本方程，提出一种新的以地固系位置速度为状态量的状态转移矩阵解析计算方法并分析其计算量，证明了地固系自主定轨算法可以进一步降低定轨算法复杂度，且无需上注地球定向参数进行坐标转换。使用仿真数据和真实在轨测量数据的测试结果表明，提出的地固系分布式自主定轨算法与传统算法定轨精度相当，说明了所提算法的可行性，验证了星载自主定轨原型软件的有效性。     

Integrating advanced mobility into lunar surface exploration

With growing knowledge of the lunar surface environment from recent robotic missions, further assessment of human lunar infrastructures and operational aspects for surface exploration become possible. This is of particular interest for the integration of advanced mobility assets, where path planning, balanced energy provision and consumption as well as communication coverage grow in importance with the excursion distance. The existing modeling and simulation tools for the lunar surface environment have therefore been revisited and extended to incorporate aspects of mobile exploration. An extended analysis of the lunar topographic models from past and ongoing lunar orbital missions has resulted in the creation of a tool to calculate and visualize slope angles in selected lunar regions. This allows for the identification of traversable terrain with respect to the mobile system capabilities. In a next step, it is combined with the analysis of the solar illumination conditions throughout this terrain to inform system energy budgets in terms of electrical power availability and thermal control requirements. The combination of the traversability analysis together with a time distributed energy budget assessment then allows for a path planning and optimization for long range lunar surface mobility assets, including manned excursions as well as un-crewed relocation activities. The above mentioned tools are used for a conceptual analysis of the international lunar reference architecture, developed in the frame of the International Architecture Working Group (IAWG) of the International Space Exploration Coordination Group (ISECG). Its systems capabilities are evaluated together with the planned surface exploration range and paths in order to analyze feasibility of the architecture and to identify potential areas of optimization with respect to time-based and location-based integration of activities.     

The use of three-impulse transfer to insert the spacecraft into the high Moon Artificial Satellite orbits

The problem of the optimal spacecraft’s insertion from the Earth into the high circular polar Moon Artificial Satellite’s orbit (MAS) with a radius of 4000–8000 km has been investigated. A comparison of single- and three-impulse insertion schemes has been performed. The analysis was made taking into account the disturbances from the lunar gravity field harmonics and the gravity fields of the Earth and the Sun, as well as the engine’s limited thrust. It has been shown that the three-impulse transfer from the initial selenocentric hyperbola of the approach into the considered final high MAS orbit is noticeably better with respect to the final mass than the ordinary single-impulse deceleration. The control parameters that implement this maneuver and provide nearly the same energy expenses as in the Keplerian case have been presented. It was found that, in contrast to the Keplerian case, in the considered case of the real gravity field, there is the optimal maximum distance of the maneuver. Recently, the Moon exploration problem became actual again.     

Strategies of developing advanced lunar power systems

Electric and thermal power have to be available at the base site on the lunar surface before the first lunar crew arrives. Unlimited solar energy is available during the lunar day, but this must be stored for use during the lunar night unless nuclear energy systems are available. State-of-the-art candidate systems are reviewed and the production of solar cells on the moon is discussed. Various options for developing a lunar power plant are proposed. These must be simulated and optimized in a real lifecycle systems scenario to provide operations and cost data essential for choosing a strategy.     

Method based on artificial excitation of characteristic radiation by an electron beam for remote X-ray spectral elemental analysis of surface rocks on atmosphereless celestial bodies

This article, like our previous one [1], is devoted to advanced space technology concepts. It evaluates the potential for developing active systems to conduct a remote elemental analysis of surface rocks on an atmosphereless celestial body. The analysis is based on the spectrometry of characteristic X-rays (CXR) artificially excited in the surface soil layer. It has been proposed to use an electron beam injected from aboard a spacecraft orbiting the celestial body (or moving in a flyby trajectory) to excite the CXR elements contained in surface rocks. The focus is on specifying technical requirements to the parameters of payloads for a global mapping of the composition of lunar rocks from aboard of a low-orbiting lunar satellite. This article uses the results obtained in [2], our first study that shows the potential to develop an active system for a remote elemental analysis of lunar surface rocks using the above method. Although there has been interest in our research on the part of leading national academic institutions and space technology developers in the Soviet Union, the studies were discontinued because of the termination of the Soviet lunar program and the completion of the American Apollo program.     

Strategies of developing advanced lunar power systems

Electric and thermal power have to be available at the base site on the lunar surface before the first lunar crew arrives. Unlimited solar energy is available during the lunar day, but this must be stored for use during the lunar night unless nuclear energy systems are available. State-of-the-art candidate systems are reviewed and the production of solar cells on the moon is discussed. Various options for developing a lunar power plant are proposed. These must be simulated and optimized in a real life-cycle systems scenario to provide operations and cost data essential for choosing a strategy.     

月表地形环境对着陆器热控系统的影响分析

月球着陆器着陆月表之后的外热流环境与地球轨道卫星外热流环境的最大区别在于月表红外热流的影响.分析月表地形特征是确定月表红外热流影响范围的基础.文章根据调研分析选择了月海地区可能存在的平坦月表、斜坡、石块和撞击坑4种典型地形特征以及月表的热物性参数,并通过计算分析得出了月表地形对着陆器热控系统的一些影响规律--其中平坦月表和斜坡的影响最大,石块和小型撞击坑的影响可以忽略,大型撞击坑可以等效为斜坡的影响.