Compensation of gravitational attraction disturbance to pure gravity orbit for Inner Formation Flying System

The Inner Formation Flying System (IFFS) consisting of an outer satellite and an inner satellite which is a solid sphere proof mass freely flying in the shield cavity can construct a pure gravity orbit to precisely measure the earth gravity field. The gravitational attraction on the inner satellite due to the outer satellite is a significant disturbance source to the pure gravity orbit and is required to be limited to 10⁻¹¹ m s⁻² order. However, the gravitational disturbance force was on 10⁻⁹ m s⁻² order actually and must be reduced by dedicated compensation mass blocks. The region of relative motion of the inner satellite about its nominal position is within 1 cm in dimension, which raises the complexity of the compensation blocks design. The iterative design strategy of the compensation blocks based on reducing the gravitational attraction at the nominal position of the inner satellite is presented, aiming to guarantee the gravitational force in the relative motion region within requirements after the compensation. The compensation blocks are designed according to the current status of IFFS, and the gravitational disturbance force in the region is reduced to 10⁻¹¹ ms⁻² order with minimized adding mass.     

Calculation of quasisynchronous orbits of a spacecraft around Phobos for solving the problem of landing on its surface

The problem of selecting quasi-synchronous orbits of a spacecraft around Phobos is considered. These quasi-synchronous orbits are far (with respect to the Hill’s sphere) quasi-satellite orbits with retrograde rotation in the restricted three body problem. The orbit should pass through a given point at a specified time instant. It should also possess a property of minimum distance from the Phobos surface at every passage above the region of planned landing. The equations of dynamics are represented in the form describing the orbit as a combination of motions in two drifting ellipses, inner and outer ellipses. The center of the outer ellipse is located on the inner ellipse. A formula is derived that relates averaged values of half-axes of the inner and outer ellipses. It is used for construction of the first approximation of numerically designed orbit, which makes it possible to simplify and speed up the computing process. The tables of initial conditions obtained as a result of calculations are presented.     

载人航天器生活舱内湿度场的稳态数值模拟

研究载人航天器生活舱内部的湿度情况具有重要意义。由于生活舱设备外形和结构相当复杂,难以建模,因此按照速度梯度、湿度梯度较大的优先原则,合理选择简化舱体结构及舱内设备。采用基于有限元网格的控制容积法对计算区域离散,对离散方程组求解得出水汽质量浓度分布后,结合舱内温度场通过一定的转换关系最终得到了生活舱内部的湿度分布。数值分析结果表明生活舱内处于相对湿度较高的状态下,提示还需进一步采取措施以改善舱内湿度水平。     

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.     

磁悬浮飞轮优化设计与实验

针对在轨磁悬浮飞轮转子系统,考虑发射主动段振动工况,提出一种基于灵敏度分析的多学科优化设计方法。通过比较分析内、外锁紧方案对锁紧状态飞轮转子系统共振频率的影响,得到了较优的抱式外锁紧方案。在此基础上,对飞轮转子进行灵敏度分析,并选择高灵敏度结构参数作为优化设计变量。考虑其结构强度、转动惯量、控制性能、自由状态和锁紧状态一阶共振频率,以飞轮转子质量最小为优化目标,利用序列二次规划法对其进行多学科优化设计。优化结果表明,在满足发射振动工况下,飞轮转子极转动惯量与质量的数值比达到最大为0.0070,比初始值0.0064提高了9.4%。该优化方法提高了飞轮转子设计的可靠性和效率,其首次在轨实验成功对我国磁悬浮飞轮空间应用具有重要意义。     

Venus transfer design by combining invariant manifolds and low-thrust arcs

The design of interplanetary trajectories based on patched circular restricted three body models is gradually becoming a valuable alternative to the classical patched conic approach. The main advantage offered by such a model is the possibility to exploit the manifold dynamics to move naturally far from or toward a body. Generally, propulsive maneuvers are required to match these structures. Low-thrust arcs offer the possibility to have a significant propellant mass reduction when moving from manifold to manifold. The aim of this paper is to present a methodology to design low-thrust trajectories between two planetary orbits connecting the manifolds of two circular three body systems. The approach is based on a grid search on the main parameters governing the solution to identify those trajectories moving within the manifold images on given Poincarè sections. The value of the Jacoby constant of the target libration point periodic orbit is chosen as stop condition for the thrusting phases. Ballistic arcs follow up to the proper Poincarè section intersection. A grid search for an Earth to Venus transfer is presented as test case.     

Collisions of stars by oscillating orbits of the second kind

Ordinary estimations of the number of star collisions in our galaxy—by simple kinematic considerations—lead to a very small number of such collisions: about one or even less every millions of years. However star collisions can occur through the following indirect way which has a much higher probability. (a) Binary stars are very common in our galaxy, about 30–50% of the stars. (b) If two binary stars meet a triple system can be formed by an ordinary exchange type motion. (c) A triple system is generally decomposed into the “inner orbit” (i.e. the relative orbit of the two nearest stars) and the “outer orbit” (i.e. the relative orbit of the third star with respect to the center of mass of the two nearest stars). The major axes of these two orbits have generally small perturbations and it is the same for the eccentricity of the outer orbit. On the contrary, if the relative inclination of the two orbits is large, the perturbations of the eccentricity of the inner orbit are important and can even in some cases lead to an eccentricity equal to one, that is to a collision of the two stars of the inner orbit.Such orbits can be called “oscillating orbits of the second kind”, indeed the first oscillating orbits—conceived by Khilmi and described for the first time in an example by Sitnikov—have unbounded mutual distances r_ij, but the system always come back to small sizes, it has an infinite number of very large expansions followed by strong contractions and, in the three-body case, an upper bound of lim inf (r_1.2 + r_1.3 + r_2.3) can be given in terms of the three masses and the integrals of motion. For the oscillating orbits of the second kind the mutual distances r_ij are bounded, but the velocities are unbounded (i.e. lim inf r_ij = 0 for at least one r_ij) and the system goes to a collision if the bodies have non-zero radius even small. The analytical study of the oscillating orbits of the second kind is a part of the general analytical study of the three-body problem, a part which must be valid for large eccentricities and large inclinations. The use of Delaunay's variables and of a Von Zeipel transformation lead to a first order integrable approximation, valid for any eccentricities and any inclinations, and giving the following results: (a) The oscillating orbits of the second kind occur when the angular momentum of the outer orbit has a modulus sufficiently close to the modulus of the total angular momentum of the three-body system. Hence these orbits occur for inclinations in the vicinity of 90°. (b) The oscillating orbits represent a set of positive measure of phase space and the first order study allows to give a rough estimation of the probability of collisions—even for stars of infinitely small radius. This probability, for given initial major axes and eccentricities and for isotropic arbitrary initial orientations, is generally of the order of m₃$\text{R}\text{M}$ (m₃ being the mass of the outer star, M the total mass and R the ratio of the period of the inner orbit to that of outer orbit).One question remains to be solved: how many collisions of stars are due to that phenomenon? That question is difficult because the probability of formation of a triple system by a random meeting of two binaries is very uneasy to estimate. However it seems that, compared to the usual evaluations based on pure kinematic considerations without gravitational effects, the number of collisions must be multiplied by a factor between one thousand and one million.     

非结构网格技术应用于固体火箭发动机内流场数值模拟

从二维欧拉方程出发，采用非结构网格技术对固体火箭发动机内流场进行数值模拟，研究了点火初期、燃烧过程中和燃烧结束时固体火箭发动机燃烧室不规则物理区域内形成的复杂流动，分析了流场结构及特性，并对在这一领域应用非结构网格技术存在的问题进行分析。     

超高速撞击中靶后碎片云的内外边界方程研究

研究超高速撞击防护靶形成的碎片云的形状模型对航天器防护结构设计具有重要意义。文章用数值仿真的方法研究靶后碎片云的内外边界模型:将碎片云根据其材料来源分为靶板碎片云和弹丸碎片云,分别根据各碎片云的分布特点建立碎片云形状的内外边界方程;在双纽线边界方程的基础上,修正为由双纽线方程和圆弧方程组成的边界方程。对比整体碎片云的适用率表明,由双纽线方程和圆弧方程组成的边界方程比单一的双纽线边界方程描述碎片云实际分布更为贴合。     

发动机药柱和推进剂方坯老化性能相关性研究

通过长期贮存的CTPB推进剂方坯性能变化和发动机中推进剂药柱性能变化比较,研究了发动机药柱和推进剂方坯老化性能的相关性,发现发动机中不同位置的推进剂性能的变化有明显差异,内层推进剂“变软”的速率比外层慢得多。当外层推进剂强度降低较大时,内层推进剂仍有较高的保持率,几乎和推进剂初始性能相同,并且强度由内向外逐渐变化。因此。单用推进剂方坯的老化性能难于推断发动机药柱的寿命,并对这一现象对发动机寿命的影响进行了讨论。     

Did evolution select a nonrandom "alphabet" of amino acids?

The last universal common ancestor of contemporary biology (LUCA) used a precise set of 20 amino acids as a standard alphabet with which to build genetically encoded protein polymers. Considerable evidence indicates that some of these amino acids were present through nonbiological syntheses prior to the origin of life, while the rest evolved as inventions of early metabolism. However, the same evidence indicates that many alternatives were also available, which highlights the question: what factors led biological evolution on our planet to define its standard alphabet? One possibility is that natural selection favored a set of amino acids that exhibits clear, nonrandom properties-a set of especially useful building blocks. However, previous analysis that tested whether the standard alphabet comprises amino acids with unusually high variance in size, charge, and hydrophobicity (properties that govern what protein structures and functions can be constructed) failed to clearly distinguish evolution's choice from a sample of randomly chosen alternatives. Here, we demonstrate unambiguous support for a refined hypothesis: that an optimal set of amino acids would spread evenly across a broad range of values for each fundamental property. Specifically, we show that the standard set of 20 amino acids represents the possible spectra of size, charge, and hydrophobicity more broadly and more evenly than can be explained by chance alone.     

航天器机动时DGMSCMG磁悬浮转子干扰补偿控制

双框架磁悬浮控制力矩陀螺（DGMSCMG）具有寿命长、综合效益好等突出优势,但航天器机动时,航天器及DGMSCMG内、外框架系统的转动均导致磁悬浮高速转子产生一定的耦合运动,影响磁悬浮转子系统的稳定性,同时使输出力矩精度下降,从而严重影响航天器姿态控制的精度。本文建立了基于DGMSCMG的航天器动力学模型,分析航天器、外框架、内框架、磁悬浮转子四者之间的动力学耦合关系。针对磁悬浮转子的非线性耦合干扰,提出一种基于复合控制的补偿方法,通过磁轴承产生相应的电磁力,对陀螺耦合力矩和惯性耦合力矩进行补偿控制。仿真结果表明,干扰补偿控制能有效抑制航天器及框架对磁悬浮转子的耦合干扰,也有效提高了磁悬浮转子系统的稳定性。     

Hypervelocity impact testing of advanced materials and structures for micrometeoroid and orbital debris shielding

A series of 66 hypervelocity impact experiments have been performed to assess the potential of various materials (aluminium, titanium, copper, stainless steel, nickel, nickel/chromium, reticulated vitreous carbon, silver, ceramic, aramid, ceramic glass, and carbon fibre) and structures (monolithic plates, open-cell foam, flexible fabrics, rigid meshes) for micrometeoroid and orbital debris (MMOD) shielding. Arranged in various single-, double-, and triple-bumper configurations, screening tests were performed with 0.3175 cm diameter Al2017-T4 spherical projectiles at nominally 6.8 km/s and normal incidence. The top performing shields were identified through target damage assessments and their respective weight. The top performing candidate shield at the screening test condition was found to be a double-bumper configuration with a 0.25 mm thick Al3003 outer bumper, 6.35 mm thick 40 PPI aluminium foam inner bumper, and 1.016 mm thick Al2024-T3 rear wall (equal spacing between bumpers and rear wall). In general, double-bumper candidates with aluminium plate outer bumpers and foam inner bumpers were consistently found to be amongst the top performers. For this impact condition, potential weight savings of at least 47% over conventional all-aluminium Whipple shields are possible by utilizing the investigated materials and structures. The results of this study identify materials and structures of interest for further, more in-depth, impact investigations.     

基于大气阻力实时辨识的Drag-free卫星最优控制研究

研究了基于大气阻力实时辨识的Drag-free卫星最优控制。将Drag-free卫星和其内部的验证质量等效为两颗内、外编队卫星并建立动力学方程,推导了基于卫星和验证质量的相对运动状态观测值反演大气阻力的算法。建立了状态最优调节器模型,采用动态规划求解经典的二次型最优控制。对低轨圆轨道Drag-free卫星的仿真计算结果表明方法的求解精度较高,计算消耗较小。     

The use of a movable telescoping end mass system for the time-optimal control of spinning spacecraft

The time-optimal control of a spin-stabilized spacecraft with a movable telescoping appendage (boom), is considered analytically and numerically. The motion of a control mass at the end of the boom is determined such that the terminal time will be minimized for two-axis control of a symmetric spacecraft. The equations of rotational motion are linearized about the desired state of spin about the symmetry axis. The equations for the transverse angular velocity components have the form of a coupled two dimensional harmonic oscillator with boom motion as a control force. The control function which brings the system to the desired state is known to be a series of positive and negative pulses. If the initial state is such that the system can be driven to rest in a single switch, the responses, switching and final times, and required boom motion may be determined analytically. Some typical numerical results based on these solutions are discussed.     

液体火箭发动机阻尼机架减振特性研究

对具有橡胶——惯性结构的火箭发动机机架简化模型中动态载荷传播特性进行了研究。机架内动态载荷在依次通过橡胶-惯性结构粘贴位置时，携带的部分动能传递到橡胶材料内，导致机架内应力幅值显著低于能量交换前的应力。通过分析橡胶材料内外表面的速度，发现速度差引起了橡胶阻尼材料沿轴向发生剪切加载和卸载过程，并导致机架传递到橡胶材料的动态能量发生剪切变形的过程中被耗散吸收；此外，机架内动态载荷大小与机架内质点速度的变化量有关，机架内质点动能减小及运动速度降低是引起机架动态载荷幅值降低的根本原因。     

树脂含量对芳纶纤维/环氧复合材料性能的影响

研究了树脂含量对芳纶纤维／环氧复合材料的性能的影响，主要是浸胶纤维强度，压力容器中纤维强度的发挥系数，容器特性系数pV/Wc。给出了按混合定律计算复合材料强度的修正系数，提出了对容器进行表面处理以改善其内外层含胶量的均匀度，并进行了实验，实验证明，容器内外层树脂含量分布更为均匀，容器的性能得到提高。     

电容放电式火工品点火电路参数设计与仿真

文章介绍了火工品电容放电式点火电路,阐述了电容充放电原理和充放电电路的参数计算方法;仿真分析了充电电压和储能电容等参数对点火电路放电性能的影响。通过仿真分析得出,合理选择储能电容和充电电压能够获得较好的点火性能。     

航天器瞬态热平衡试验技术的新探索

文章以某型号航天器热平衡试验为工程应用背景,在正样星中用星内仪器按在轨飞行状态工作的方法实现仪器的内热源,用表面薄膜加热器模拟航天器表面多层部分的外热流,用红外加热笼模拟散热面的外热流,藉此实现了瞬态热平衡试验.     

旋风式精过滤器内流场模拟的若干问题

概述了固体燃气发生器中应用的旋风式精过滤器内流场模拟中的主要问题，包括其网格生成、湍流模拟、两相流模拟、边界条件及数值计算方法等。对于全面展开三维两相旋转湍流的内流场模拟有一定的参考价值。最后，结合固体推进剂燃气过滤的特点进行了简要讨论。