贮箱内推进剂微重力下晃动的有限元模拟及试验

采用有限体积(VOF)法 任意拉格朗日．欧拉(ALE)法，建立分析液体火箭贮箱内推进剂大幅晃动的数学模型，给出了计算的基本方程。对常重力和低重力下不同形状贮箱内推进剂的非线性晃动进行了数值模拟，并用落塔法进行了微重力下贮箱内液体晃动试验。结果表明，数值模拟法与落体试验的结果相近，在工程上有一定的实用价值。     

微重力下球形贮箱内液体晃动有限元模拟

以卫星球形贮箱为研究对象,用任意拉格朗日-欧拉(ALE)法建立了数学模型,给出了计算的基本方程,并基于MSC-Dytran有限元程序,对常重力和微重力环境中贮箱内推进剂的晃动进行了数值模拟。结果表明,微重力环境中液体更易发生大幅晃动,甚至产生飞溅,工程中常用的等效力学模型已不再适用,而有限元数值模拟方法可较好地解决该问题,为卫星及大充液复杂结构航天器控制系统的设计和分析提供理论参考。     

航天器贮箱大幅液体晃动三维质心面等效模型研究

针对航天器贮箱内的大幅液体晃动建立了三维质心面等效模型,此模型将贮箱内液体等效为位于液体质心处的质点,该质心点可以在质心约束面内任意运动,同时利用质心点动力学方程及接触碰撞模型近似计算液体晃动对贮箱的作用力。质心面等效模型的数值仿真结果与基于VOF方法的流体动力学数值计算结果基本吻合,从而验证了质心面等效模型用于航天器贮箱大幅液体晃动分析的有效性。     

用多维模态理论分析航天器贮箱液体有限幅晃动力

首次将多维模态理论应用到求解航天工程中常见的圆柱贮箱液体非线性晃动问题中。针对贮箱作水平横向运动，根据Narirrmnov-Moiseev的三阶渐近假设关系，通过选取主导模态以及确定它们的阶次关系，将一般形式的无穷维模态系统降为5维渐近模态系统。在模态系统基础之上，推导出了液体作用于贮箱壁的力的简洁表达式。数值仿真结果表明，相对于线性晃动理论结果，有限幅晃动液体将对贮箱壁产生更大的横向作用力，同时对贮箱底部还作用有明显的纵向力。可将此公式应用于工程实际问题中部分充液贮箱液体有限幅晃动力的估算。     

非轴对称贮箱液体晃动的等效力学模型

针对非轴对称贮箱液体晃动问题,从一般意义上推导了作三维运动的任意三维贮箱内液体晃动的等效力学模型,采用有限元方法建立计算等效力学模型参数的数值算法,利用矩阵相似变换的性质分离液体零频以消除液体刚度矩阵的奇异性,从而将等效力学模型参数的计算归结为广义特征值问题和线性代数方程组的求解,具有较高的计算效率。通过非轴对称贮箱液体晃动的算例说明,不同的等效晃动方向可能不正交,晃动平面可能不过液体质心,因此非轴对称贮箱液体晃动的等效力学模型比轴对称贮箱液体晃动的更加复杂。最后以三棱柱形充液贮箱为例,采用商用CFD（计算流体动力学）分析软件FLOW-3D验证了在贮箱作典型三维运动时等效力学模型的有效性。
     

滑行段低温推进剂流动及换热特性对气枕压力的影响研究

运载火箭在飞行过程中需要进行姿态调整以满足入轨要求,贮箱内推进剂在外界干扰力的作用下将发生晃动,由此引入了诸如气液接触面积、蒸发、冷凝过程及推进剂流动变化等不确定影响因素。实际飞行过程尤其是进入滑行段的初始推进剂晃动对贮箱内气枕压力及推进剂流动行为具有重要影响。在调研国内外运载火箭末级飞行过程中低温贮箱压力及推进剂流动特性的基础上,建立仿真模型,采用流体体积函数方法(VOF)分析滑行段推进剂流动特性变化对贮箱气枕压力的影响。     

航天器贮箱内液体晃动阻尼研究(二):数值计算

针对航天器贮箱中燃料晃动阻尼进行估算.在液体晃动阻尼理论研究的基础上,考虑了容器壁面、液体内部以及受污染的自由液面处产生的耗散,建立了一套计算刚性容器中液体小幅晃动阻尼的有限元计算方案,并将计算值与实验值以及经验公式值进行了比较.结果发现大部分情况下的实验值落在两种极端情况的计算值之间,此外如果将液面污染的Stokes边界层阻尼以一定的加权因子引入,可以获得与由实验数据拟合的经验公式一致的结果.     

微重力下板式贮箱内液体晃动性能研究

板式贮箱内推进剂管理装置(PMD)以板式结构部件为主,其抑制液体晃动性能应能够满足卫星平台高稳定度和快速机动的需求。针对微重力下板式贮箱内液体晃动性能,建立基于有限体积(VOF)方法的气液两相流流动模型,通过数值模拟计算手段,比较分析研究了微重力环境下有防晃叶片和无防晃叶片的板式PMD晃动性能。计算分析结果表明,有防晃叶片的板式PMD抑制液体晃动性能明显优于无防晃叶片的板式PMD,使液面快速趋于稳定。同时,将计算分析结果与液体晃动微重力试验结果进行了分析比较,微重力试验结果与计算分析结果非常相近。研究结果表明,有防晃叶片的板式PMD具有更好的抑制液体晃动功能,能够显著地抑制贮箱内推进剂晃动,以满足卫星平台的高稳定度和快速机动的需求。     

火箭推进剂液体晃动关机响应的数值仿真

本文采用VOF方法对火箭发动机关机阶段所引起的液体晃动进行了数值仿真，以氧化剂贮箱为例，得到了贮箱无阻尼环和带阻尼环时的晃动仿真图和液体质心变化曲线，检验了VOF方法处理碎波的能力。     

超低重力环境下液体推进剂的特性

通讯卫星上液体推进剂的质量比值大会对卫星的姿态产生影响,这与卫星的力学特性和流体运动间的相互作用有关.卫星贮箱内液体燃料的自由表面波动或液面晃动是一个公认的问题.本文提供了作用于贮箱上,且量值相当大的重力或等量加速度力的数值结果,但这种力变小时就应把其它的力考虑进去,如液体燃料与液面上部的气体之间的面际张力.在ALCATEL SPACE和LEMMA的共同努力下,已开发出ANASLOSH软件包,其目的是给任何即便是很小的加速度力提供准确测定晃动模式的方法.首先,考虑表面张力及液-气界面与贮箱壁交点处的接触条件,确定稳态平衡表面的形状,接着,计算出了自由表面的微小振荡.电推进机动飞行过程中通讯卫星晃动模式的确定就是其一项工业应用.     

基于刚体转动四元数模型的航天器再定向方案

介绍了一种新的刚体转动的四元数模型。这些模型可以提高再定向控制的对边界条件的适应能力，在此基础上还介绍了这些模型的一些实际应用。     

On solving the problems of optimization of trajectories of many-revolution orbit transfers of spacecraft

The problem of optimal control over many-revolution spacecraft orbit transfers between circular coplanar orbits of satellites is considered. The spacecraft flight is controlled by a thrust vector of a jet engine with restricted thrust (JERT). The mass expenditure is minimized at a limited time of flight. The optimal control problem is solved based on the maximum principle. The boundary value problem of the maximum principle is solved numerically using the shooting method. A modified computation scheme of the shooting method is suggested (multi-point shooting), as well as a method (correlated with the scheme) of choosing the initial approximation with the use of a solution to the optimization problem in the impulse formulation. The scheme and method allow one to construct many-revolution spacecraft orbit transfers.     

热变形对一种天线指向机构指向精度的影响分析

热变形对天线指向机构的指向精度具有重要影响,但利用传统的设计方法很难实现对指向精度的准确计算,尤其是热变形对指向精度影响的研究工作更是少见。为此建立了机构的三维实体热-结构耦合有限元分析模型,根据有限元分析结果并结合机构的运动链关系,构建了各关键运动副运动误差的计算方法;在此基础上实现了机构指向误差的计算,研究了运动副名义、最大以及最小配合间隙条件下,温度载荷对各关键运动副运动误差和指向精度的影响。研究结果表明,温度载荷对关键运动副的运动误差以及机构的指向精度具有重要影响,随着温度载荷偏离常温数值的增加,指向误差不断增大,名义配合间隙条件下的指向误差最大。     

Determination of parameters of attitude motion of a small communication satellite using the data of measurements of current of solar panels

A communication satellite (small spacecraft) injected into a geosynchronous orbit is considered. Flywheel engines are used to control the rotational spacecraft motion. The spacecraft after the emergency situation has passed into a state of uncontrolled rotation. In this case, no direct telemetric information about parameters of its rotational motion was accessible. As a result, the problem arose to determine the rotational satellite motion according to the available indirect information: current taken from the solar panels. Telemetric measurements of solar panel current obtained on the time interval of a few hours were simultaneously processed by the least squares method integrating the equations of rotational satellite motion. We present the results of processing 10 intervals of the measurement data allowing one to determine the real rotational spacecraft motion and to estimate the total angular momentum of flywheel engines.     

攻角变化对超音速进气道再起动特性的影响

以二维非定常可压缩流的N-S方程为控制方程,采用SST k-ω湍流模型对攻角变化引起的超音速进气道再起动过程进行了数值模拟,研究了攻角变化对超音速进气道再起动特性的影响.结果表明,当超音速进气道不起动时,可通过合理改变攻角实现进气道的再起动工作;超音速进气道的再起动攻角随攻角变化速率的增大近似呈线性增加;攻角变化速率较...     

电缆屏蔽金属网屏蔽效能的工程计算

金属网由于柔韧性好、使用方便等优势而在电缆屏蔽工程上得到广泛应用。本文对电缆用的金属网屏蔽效能进行工程计算,通过计算,对金属网不同材料、不同规格、单层与双层屏蔽及有缝隙情况下的屏蔽效能进行比较与分析,对工程应用有一定的指导意义。     

多通道功率放大器中幅度和相位一致性影响分析

文章利用传输矩阵方法,建立了MPA分析模型,并通过引入幅度和相位随机性,结合各路放大器幅度和相位组合概率,对通道隔离度的影响进行了计算和分析。     

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

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

Determination of an Estimation of Parameters of Motion of a Spacecraft on the Stage of Its Descent on the Basis of Integration of Inertial and Optoelectronic Navigation Systems

The problem of determining an estimation of current spacecraft motion parameters on the stage of its descent is solved on the basis of readings of an optoelectronic navigation system and additional information from the self-contained inertial navigation system at the most general assumptions on the nature of motion of the object and the statistical properties of the noise of the measuring devices.