镁铝贫氧推进剂的能量分析

利用冲压发动机热力计算程序，对镁铝中能贫氧推进剂的能量特性进行了系统研究，并探讨了高能级分对该类贫氧推进剂能量性能的影响，研究结果表明，在镁铝中能贫氧推进剂中增加镁粉含量（或减少铝粉含量）贫氧推进剂的比冲下降，在一定空燃比范围内，增加空燃比有助于提高冲压发动机的比冲。提高贫氧推进剂中ＣＬ－２０和硼粉含量，可以显著提高共能量，而采用叠氮类含能粘合剂取代惰性粘合剂时贫氧推进剂的比冲降低。     

Analysis of the Accuracy of Ballistic Descent from a Circular Circumterrestrial Orbit

The problem of the transportation of the results of experiments and observations to Earth every so often appears in space research. Its simplest and low-cost solution is the employment of a small ballistic reentry spacecraft. Such a spacecraft has no system of control of the descent trajectory in the atmosphere. This can result in a large spread of landing points, which make it difficult to search for the spacecraft and very often a safe landing. In this work, a choice of a compromise scheme of the flight is considered, which includes the optimum braking maneuver, adequate conditions of the entry into the atmosphere with limited heating and overload, and also the possibility of landing within the limits of a circle with a radius of 12.5 km. The following disturbing factors were taken into account in the analysis of the accuracy of landing: the errors of the braking impulse execution, the variations of the atmosphere density and the wind, the error of the specification of the ballistic coefficient of the reentry spacecraft, and a displacement of its center of mass from the symmetry axis. It is demonstrated that the optimum maneuver assures the maximum absolute value of the reentry angle and the insensitivity of the trajectory of descent with respect to small errors of orientation of the braking engine in the plane of the orbit. It is also demonstrated that the possible error of the landing point due to the error of specification of the ballistic coefficient does not depend (in the linear approximation) upon its value and depends only upon the reentry angle and the accuracy of specification of this coefficient. A guided parachute with an aerodynamic efficiency of about two should be used at the last leg of the reentry trajectory. This will allow one to land in a prescribed range and to produce adequate conditions for the interception of the reentry spacecraft by a helicopter in order to prevent a rough landing.     

Optimization of the Flight of a Spacecraft with a Solar Sail from the Earth to Mars with a Perturbation Maneuver near Venus

The trajectories of the fastest flight of a spacecraft (SC) with a solar sail from the Earth's sphere of activity to the Martian sphere of activity including the section of a perturbation maneuver near Venus are investigated. The planetary spheres of activity are assumed to be point-like; i.e., the maneuver section and the initial and final positions of the SC coincide with the corresponding positions of the planets. The initial velocity of the SC is assumed to be equal to the Earth's velocity, so that no leveling of the velocities of the SC and Mars in the final point of the flight is required. The perturbation maneuver is considered as a jump of the heliocentric velocity of the SC at the point of its contact with Venus, which does not change the magnitude of its Venus-centric velocity. The orbits of planets are assumed to be circular and coplanar; the SC trajectory lies at the plane of these orbits. The sail is planar with a specularly reflecting surface. The trajectories of optimum flights are determined as a result of solving the boundary value problem of the Pontryagin maximum principle. The families of solutions to this problem depending on the initial angular positions of Venus and Mars are constructed by the method of continuation over a parameter.     

The Use of Quaternions in the Optimal Control Problems of Motion of the Center of Mass of a Spacecraft in a Newtonian Gravitational Field: I

The problem of optimal control is considered for the motion of the center of mass of a spacecraft in a central Newtonian gravitational field. For solving the problem, two variants of the equations of motion for the spacecraft center of mass are used, written in rotating coordinate systems. Both the variants have a quaternion variable among the phase variables. In the first variant this variable characterizes the orientation of an instantaneous orbit of the spacecraft and (simultaneously) the spacecraft location in this orbit, while in the second variant only the instantaneous orbit orientation is specified by it. The suggested equations are convenient in the respect that they allow the general three-dimensional problem of optimal control by the motion of the spacecraft center of mass to be considered as a composition of two interrelated problems. In the first variant these problems are (1) the problem of control of the shape and size of the spacecraft orbit and (2) the problem of control of the orientation of a spacecraft orbit and the spacecraft location in this orbit. The second variant treats (1) the problem of control of the shape and size of the spacecraft orbit and the orbit location of the spacecraft and (2) the problem of control of the orientation of the spacecraft orbit. The use of quaternion variables makes this consideration most efficient. The problem of optimal control is solved on the basis of the maximum principle. Several first integrals of the systems of equations of the boundary value problems of the maximum principle are found. Transformations are suggested that reduce the dimensions of the systems of differential equations of boundary value problems (without complicating them). Geometrical interpretations are given to the transformations and first integrals. The relation of the vectorial first integral of one of the derived systems of equations (which is an analog of the well-known vectorial first integral of the studied problem of optimal control) with the found quaternion first integral is considered. In this paper, which is the first part of the work, we consider the models of motion of the spacecraft center of mass that employ quaternion variables. The problem of optimal control by the motion of the spacecraft center of mass is investigated on the basis of the first variant of equations of motion. An example of a numerical solution of the problem is given.     

固体火箭超燃冲压发动机补燃室构型的影响分析

针对不同补燃室结构参数对固体火箭超燃冲压发动机补燃室掺混燃烧性能的影响进行研究,分析各级燃烧室的长度与扩张角度对补燃室性能的影响。采用基于密度的二阶迎风格式对补燃室掺混燃烧进行模拟,湍流模型和燃烧模型分别采用SST k-ω模型和涡团耗散模型。结果表明,提高燃烧效率与降低总压损失是相互矛盾的;燃烧效率随燃烧室长度的增加而增大,随燃烧室扩张角度的增加而减小;总压恢复系数随燃烧室长度的增加而减小,随燃烧室扩张角度的增加而增大;一级燃烧室的结构参数对燃烧效率与总压恢复系数的影响最大。当补燃室的总长与出口面积一定时,以发动机的总体性能参数作为补燃室构型的优化目标,对一、二级燃烧室长度与一、三级燃烧室扩张角度进行优化。     

Determination of kinematic state of an orbiting multibody using GNSS signals

Precise attitude determination of the members of a free-flying multibody system is a not so immediate task, due essentially to the large motion of its appendages coupled with their relevant flexibility effects. In fact, sensors used to this aim in current projects, such as optical encoders usually positioned near the joints of each arm, are almost blind to these effects, and clusters of specific redundant sensors should, therefore, be required in order to reconstruct both elastic deformations and rigid motion.Satellite navigation systems (GNSS) offer a suitable and reliable solution to this problem. To exploit the phase of the signal, instead of the traditional pseudo random code, ensures a very high accuracy of the order of magnitude of centimeter. Such a process requires the solution of an initial ambiguity problem, related to the number of integer wavelength included in the length of the member.The aim of the paper is to investigate the capability of this GNSS based technique to reconstruct the kinematics of a flexible multibody system orbiting around the Earth. This analysis requires a simulation including both the multibody dynamics and the navigation system constellation to define the satellites lines-of-sight at each time step.Concerning multibody equations of motion, a Newtonian formulation is adopted in this work. A special attention is required about the choice of the state variables. As the internal forces are associated to the relative displacements between the bodies, which are small fractions of the distance of the multibody spacecraft from the center of the Earth, the task of obtaining these forces from inertial coordinates could be impossible from a numerical point of view. So, the problem is reformulated in such a way that the equation of motion of the system contains global equations, with no internal forces, and local equations, with internal forces. In the latter, only quantities of the same order of the spacecraft dimensions are present.Accuracies achievable in LEO orbit with current GPS and upcoming Galileo systems are evaluated to show the interest of the proposed technique.     

Generalized restricted circular three-body problem as a model for dynamics of binary asteroids

Exploration of the Solar System has recently revealed the existence of a large number of asteroids with satellites, which has stimulated interest in studying the dynamics of such systems. This paper is dedicated to the analysis of the relative motion of a binary asteroid. The conditions of existence of such a system (i.e., when its components do not run away) are derived in the Introduction. Then it is assumed that the satellite has no significant effect on the motion of the main asteroid, the latter being modeled as a dumbbell-like precessing solid body. The equations of motion of this system are a two-parameter generalization of the corresponding equations of the restricted circular three-body problem. It is demonstrated that in the system under consideration there exist steady-state motions in which the small asteroid is equidistant from attracting centers at the ends of the dumbbell (an analog to triangle libration points). The conditions of existence of such motions are derived, and the positions with respect to the dumbbell are analyzed in detail. Examination of the stability of the triangle libration points is reduced to investigation of a characteristic equation of the sixth degree. The stability conditions are derived in the case when the main asteroid executes near-planar motion.     

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.     

空间太阳能电站姿态运动-结构振动耦合建模与分析

针对多旋转关节空间太阳能电站构型，利用基于能量等效原理的连续体等效方法将其等效为柔性梁模型，并考虑重力梯度影响，建立了姿态运动与结构振动的耦合动力学模型；结合Runge-Kutta 法和Newmark法的优点，提出了适用于求解姿态运动与结构振动耦合动力学方程的改进算法，相比于经典Runge-Kutta 法大幅提高了效率；利用改进算法得到了不同参数下的动力学响应。在此基础上，推导了结构振动量级随结构尺寸的六次方量级增加的规律，仿真结果表明尺寸过大引发不稳定现象；分析了姿态运动和重力梯度对结构振动频率和振幅的影响；发现了姿态运动周期受结构柔性影响而增大的现象，这种现象在低轨以及大初始姿态角下影响更为明显。     

Architectural design for space tourism

The paper describes the main issues for the design of an appropriately planned habitat for tourists in space.Due study and analysis of the environment of space stations (ISS, MIR, Skylab) delineate positive and negative aspects of architectonical design. Analysis of the features of architectonical design for touristic needs and verification of suitability with design for space habitat.Space tourism environment must offer a high degree of comfort and suggest correct behavior of the tourists. This is intended for the single person as well as for the group. Two main aspects of architectural planning will be needed: the design of the private sphere and the design of the public sphere.To define the appearance of environment there should be paid attention to some main elements like the materiality of surfaces used; the main shapes of areas and the degree of flexibility and adaptability of the environment to specific needs.     

飞行过载下燃烧室凝相粒子沉积特征数值研究

固体发动机燃烧室部分凝相粒子在飞行过载作用下产生沉积,严重影响发动机工作性能。通过确定燃烧室粒子参数和建立燃烧室内两相流场数值方法,获得了发动机不同轴向过载下粒子运动及沉积规律。与试验数据进行对比分析,验证了计算方法的准确性。数值结果表明,随着轴向过载增大,后封头及喷管潜入段粒子沉积质量逐渐增大。沉积粒子粒径大于凝相粒子平均值,即粒子粒径越大,越容易沉积。轴向过载增大,减小了粒子在发动机内部的驻留时间,凝相粒子平均驻留时间均大于0.13 s。     

带机械臂的月面巡视探测器的路径规划方法研究

针对巡视探测器路径规划和机械臂路径规划的不同和两者需要依次执行的连接需求,文章将月面巡视探测器的路径规划过程分为三个阶段：巡视探测器的路径规划、机械臂的路径规划、器臂动态联合的路径规划,针对不同的规划分别提出了不同的规划方法,并进行了仿真验证。研究结果表明,巡视探测器的整个路径规划是一个复杂的运算规划过程,在非结构化月面环境下,基于月面三维数字高程图采用改进的启发式搜索（A＊）算法,可以比较高效地完成巡视探测器的路径规划;在采用蒙特卡罗法建立机械臂可达工作空间的基础上,可以比较简单、准确地获取机械臂的规划路径;巡视探测器整个就位探测过程的实现需要两种规划的动态联合。     

复合材料层合板热-力作用下的失效研究

碳纤维复合材料力学性能优异,在航空航天等领域广泛使用,其在热-力联合作用下的损伤失效研究对于结构的损伤破坏和强度预测具有重要意义。发展了热力耦合条件下复合材料结构渐进损伤分析方法,建立了三维有限元热烧蚀模型,并验证了计算模型的可靠性;采用三维Hashin失效准则,结合材料刚度突然退化模式,建立了失效分析模型,仿真分析了热-力联合作用下复合材料层合板损伤演化全过程。结果表明,该方法不仅能够较好地模拟复合材料层合板从局部失效的萌生、扩展直至结构完全失效的全过程,而且可以直观地显示结构的损伤失效模式,预测结构在不同条件下的承载能力。     

Evolution of the Orbit of an Earth Satellite with a Solar Sail

The efficiency of using the light pressure of solar radiation for increasing the semimajor axis of the orbit of an Earth Satellite carrying a solar sail is estimated. The orbit is nearly circular and has an altitude of about 900 km. The satellite is in the mode of single-axis solar orientation: it rotates at an angular velocity of 1 deg/s around the axis of symmetry, which traces the direction to the Sun. This mode is maintained by the solar sail, which serves in this case as a solar stabilizer. The following method of increasing the semimajor axis of the orbit (which is equivalent to increasing the total energy of the satellite's orbital motion) is considered. On those sections of the orbit, where the angle between the light pressure force acting upon the sail and the vector of geocentric velocity of the satellite does not exceed a specified limit, the sail is functioning as a solar stabilizer. On those sections of the orbit, where the above-indicated angle exceeds this limit, the sail is furled by way of turning the edges of the petals towards the Sun. Such a control increases the semimajor axis by more than 150 km for three months of flight. In this case, the accuracy of solar orientation decreases insignificantly.     

Convection of liquid with internal heat release in a rotating container

The convection of heat-generating fluid in a rotating horizontal cylinder is experimentally investigated. The threshold of convection excitation, the structure of convective flows and the heat transfer in the cylinder depending on the heat release capacity, liquid viscosity and aspect ratio of the cavity are studied. It is found that the average convection is excited by the thermovibrational mechanism —the gravity force, rotating in the cavity frame, produces the oscillations of non-isothermal fluid relative to the wall, which in turn result in excitation of mean convective flows. It is shown that the structure of convective flows depends on the dimensionless velocity of rotation. At relatively low rotation velocity the convection develops in the form of a periodic system of vortices regularly distributed along the cylinder axis. The threshold of excitation (critical value of vibration parameter) of three-dimensional vortex structures grows with rotation velocity. Above some definite rotation velocity the convection develops as two-dimensional rolls parallel to the axis of rotation. The threshold of two-dimensional structures excitation does not depend on the rotation velocity. Besides the structure of thermal convective flows the analysis of the relatively weak currents generated by the inertial waves below the threshold of convection is performed.     

Stochastic motion of geosynchronous satellites

Long-period evolution of the orbital elements of geosynchronous objects moving in the neighborhood of separatrices of the region of libration resonance is studied. At the phase plane “the longitude of the subsatellite point-semi-major axis” the regions of migration of separatrices where the separatrices can move due to the influence of perturbations are isolated. The dimensions of these regions and the distances between them have the same order of magnitude, several kilometers along the semi-major axis. The regions of migration do not overlap at the initial value of the orbit inclination 0°. When inclination increases, gradual overlapping of these regions takes place. On average, the width of the zones of stochasticity increases when the area to mass ratio increases and decreases when inclination of the orbit increases. The zones of stochasticity corresponding to different separatrices are separated from each other. A catalog of the passive geosynchronous objects moving in the neighborhood of separatrices is compiled. The evolution of the objects of the catalog is studied.     

Microscope in orbit calibration procedure for a test of the equivalence principle at 10(-15)

The scientific objectives of the MICROSCOPE space mission impose a very fine calibration of the on-board accelerometers. However the required performance cannot be achieved on ground because of the presence of high disturbing sources. On-board the CHAMP satellite, accelerometers similar in the concept to the MICROSCOPE instrument, have already flown and analysis of the provided data then allowed to characterise the vibration environment at low altitude as well as the fluctuation of the drag. The requirements of the in-orbit calibration procedure for the MICROSCOPE instrument are demonstrated by modelling the expected applied acceleration signals with the developed analytic model of the mission. The proposed approach exploits the drag-free system of the satellite and the sensitivity of the accelerometers. A specific simulator of the attitude control system of the satellite has been developed and tests of the proposed solution are performed using nominal conditions or disturbing conditions as observed during the CHAMP mission.     

多自由度微振动环境时域波形复现的数值仿真

为满足航天器微振动环境模拟的需要,开展了多自由度微振动时域波形复现控制方法研究。首先,介绍了基于时域波形复现的多自由度微振动环境模拟控制理论方法。其次,针对六自由度微振动激励系统,应用MATLAB软件建立了基于实测传递函数矩阵的多输入多输出微振动激励仿真系统,针对微振动时域波形复现闭环控制过程进行了算法编程,并给出了仿真的闭环控制流程图。最后,通过算例对多自由度微振动时域波形复现进行了数值仿真,以给定的白噪声为输入,模拟对实际存在的系统非线性、测量误差等影响因素的控制效果。仿真结果验证了多自由度微振动时域波形复现控制方法的可行性及有效性,所得结论可以为研究多自由度微振动时域波形复现控制系统提供参考。     

适于惯导系统初始对准的神经网络实时算法研究

通常卡尔曼滤波器被用于解决惯导系统的初始对准问题。由于卡尔曼滤波的运算时间与系统阶次的立方成正比 ,所以当系统阶次较高时 ,滤波器会失去实时性。而神经网络具有函数逼近性能 ,实时性又好。为此 ,本文研究了一种基于扩展卡尔曼滤波原理的权值更新多层神经网络学习算法 ,对此算法进行了详细的推证 ,并将该算法运用到惯导系统的初始对准过程。仿真结果表明了这种神经网络结构用于惯导系统初始对准问题的有效性 ,既真正获得了与扩展卡尔曼滤波器相同的对准精度 ,又大大提高了系统的实时性     

积分球在实验室内用于空间遥感器的辐射定标

文中介绍了国外可见光实时传输型航天遥感器在发射前使用积分球对其进行大气条件下辐射定标的技术,得出的定标系数A用于图像的辐射校正。阐述了积分球的设计原理与方法。计算了内壁涂料光谱反射率、光源功率、光源分布色温以及安放位置对积分球出射窗口输出辐亮度的影响。设计球体达到的基本指标为:球体直径Φ=1600mm,开口直径φ=300mm,最低级识分辐亮度(0.45～0.52μm)时L=1.03W/m~2.str,最高级L=12.79W/m~2.str,朗伯特性(中心点)±30°内优于1％,开口面均匀性η=0.9％,辐亮度不确定度δ=6％。用该球对CCD相机定标,结果与用“6S”软件计算的相一致。