卫星闭环姿控下的双组元四贮箱液体晃动联合仿真研究

首先介绍了基于UDP通讯协议实现的贮箱液体晃动联合仿真方法.该方法采用FLUENT计算运动贮箱的液体晃动力、力矩、质心结果，星体动力学程序在获取贮箱晃动结果后通过控制算法和策略计算星体运动的线速度和角速度.两者之间的双向数据交换，通过C/C++语言开发的通讯接口来实现.然后借助上述联合仿真方法，开展了贮箱晃动联合仿真分析，并与在轨遥感数据进行了对比，然后分析了俯仰0~25°机动工况的4贮箱液体晃动特性，最后对比分析了带孔和无孔叶片在星体动力学控制下的液体晃动特性.非沿轨迹工况下，2、3号贮箱的晃动力和力矩峰值大于1、4号贮箱，4个贮箱的Y向晃动力峰值在1.05~1.96 N，力矩值在0.05~0.36 N·m；带孔叶片对液体的搅动作用要弱于无孔叶片.无孔叶片搅动的气泡出现了明显的拉伸与离心运动，并导致更大的流体流动速度.     

侧向微重力环境下板式表面张力贮箱内推进剂晃动行为分析

侧向微重力是航天器在轨飞行时在东西位置保持和南北位置保持状态时所处的加速度环境.由于航天器贮箱内推进剂在侧向加速度环境下的重定位过程易产生晃动，因此侧向加速度环境对贮箱内推进剂管理装置(PMD)的管理能力的要求更高.为确保板式贮箱对推进剂的在轨管理能力，需要开展一系列的数值仿真与试验对PMD的管理能力进行验证.本文以板式表面张力贮箱为研究对象，采用VOF两相流模型对侧向微重力环境下重定位过程中贮箱内推进剂的运动进行数值模拟，并对推进剂运动中的晃动行为进行分析.最终，通过数值模拟结果验证了板式表面贮箱的管理能力，为微重力落塔试验与空间站试验提供参考依据.     

基于SPH方法的球形贮箱液体晃动动力学分析

针对发生液体破碎现象的航天器球形贮箱液体晃动问题,基于光滑粒子流体动力学方法（SPH方法）对N S方程进行离散近似，分别对核函数、应力张量、边界问题、邻域粒子搜索等方面进行选择确定或计算处理，最终得到基于SPH方法的流体动力学模型。通过建立包括不同燃料液体体积及不同挡板安装位置的仿真模型，分析了航天器球形贮箱的受力及局部压强。通过分析贮箱受力情况及设置测量点压强情况，得到运动参数对晃动的影响规律，为航天器球形贮箱的设计提供了理论依据。     

大容量推进剂贮箱液体晃动性能试验

对某种内带推进剂管理装置（PMD）的大容量推进剂贮箱内液体晃动性能进行试验.开展不同充液比工况下空壳贮箱液体晃动试验,分析比较晃动试验结果与采用等效动力学模型的液体晃动理论计算结果,两者结果一致吻合,试验系统可靠性和理论模型的正确性得到良好验证.然后开展了不同充液比工况下内带PMD贮箱液体晃动试验,试验结果表明贮箱内液体晃动性能受内部PMD结构影响明显.该试验研究结果为运载火箭和卫星的姿态和轨道控制系统的设计优化提供重要参考和数据支撑.     

微重力板式贮箱推进剂重定位的数值仿真

板式表面张力贮箱内推进剂重定位对确定推进剂分布情况、研究晃动影响、提高控制精度等具有重要意义．为研究板式贮箱内推进剂重定位的规律,对微重力下板式贮箱内液体重定位问题进行数值仿真．计算时使用三维非定常两相流动流体（VOF）模型,对某一板式贮箱寿命末期在不同微重力加速度情况下各种重定位过程进行数值仿真,得到各种工况下重定位的全过程,以及定位后推进剂的分布情况．数值仿真结果为板式贮箱的设计提供有利依据．     

液体非线性晃动试验与等效力学模型仿真

航天器的在轨机动可能导致液体燃料出现包括非线性晃动在内的复杂动力学行为。采用运动脉动球模型（MPBM）对液体非线性旋转晃动问题进行了研究,并通过开展地面试验验证了MPBM仿真结果的准确性。首先,在MPBM的建模过程中以晃动回复力的形式考虑重力的影响,具体是通过虚功原理将重力作用引入到脉动球的能量变化过程。然后,搭建地面晃动试验平台,实现了对球型贮箱晃动现象的准确观测以及对晃动力和晃动力矩的实时测量。最后,以试验实测的加速度激励和晃动力为参照,采用MPBM进行了仿真与对比,仿真结果较好地匹配了试验结果。文中对液体晃动的合理抽象进一步地提升了等效力学模型对复杂晃动行为的仿真能力,为实现航天器更为高效、准确的在轨控制提供基础。     

月球软着陆器动力学建模研究

建立准确、实用的月球着陆探测器动力学模型,是制导、导航和控制系统方案设计中的首要工作,文中月球软着陆探测器的推进剂贮箱采用横向对称分散布局的形式,并分析了布局方式对液体推进剂晃动的影响;将贮箱内液体考虑一阶弹簧质量模型,应用虚功率原理建立包含液体推进剂晃动的月球软着陆探测器从动力下降段至最终着陆段的动力学模型,并对动力学下降段进行了数值仿真,仿真结果可为月球软着陆探测器控制系统设计提供理论参考。     

液体火箭贮箱增压计算公式的分析

<正>一、问题的提出火箭推进剂贮箱增压系统是液体火箭上不可缺少的一个系统,增压系统保证的推进剂贮箱压力值直接关系到液体发动机工作的成败,攸关重要。然而要在地面靠模拟试验确定贮箱压力值,不仅耗资太大,并且边界条件复杂,模拟困难。因此,增压压力的理论计算就显得十分重要了。增压系统理论计算的要点是根据能量守恒原理及气体状态方程式,对进入贮箱的能量与输出的能量进行平衡,从而计算出推进剂贮箱内的增压气体压力Px。以往的增压计算,其基本公式均采用传统的方程式:     

部分充液球形中心贮箱的旋转卫星系统的运动分析

本文分析并计算了半充液球形贮箱内的液体晃动问题。以流体运动的基本方程和系统运动的Euler方程为基础,用特征函数展开法计算出液体的速度场,并在此基础上估算液体的能量耗散率及章动角的时间常数。      

罐车转向时液体晃动的仿真分析

采用VOF模型对罐车转向时液体晃动进行了数值模拟,对单室受力及侧倾力矩进行了计算.计算结果表明,无防波板时,随向心加速度增加,竖直、水平方向受力峰值及侧倾力矩均变大;随充液比增加,竖直方向受力峰值增大,水平方向受力峰值与侧倾力矩峰值先增大后减小;带防波板时,竖直、水平方向受力峰值减小,当防波板面积大于纵剖面积30%时防波板面积的增加能显著改善罐体受力;侧倾力矩峰值先增大后减小,在防波板面积等于纵剖面积30%时达到最大值,当防波板面积大于纵剖面积40%时防波板面积增加能提高防晃效果.     

微重力环境下半球形贮箱液体重定位试验研究

本文利用落塔形成的微重力环境,对半球形贮箱中液体重定位问题进行了试验研究。所用液体与半球形箱壁之间的接触角为零并在加载前形成了失重初始凹液面.本试验分侧向不加干扰及加干扰两种初始条件.通过高速摄影机观察了液体重定位的全过程并测量了重定位流速和重定位时间。试验结果表明,液体重定位呈壁流流型.实测重定位时间与理论估算值一致.      

毛细泵回路蒸发器的数值研究

文章针对毛细泵回路(CPL)蒸发器内的多孔芯结构, 建立了描述换热和流动的非稳态数值模型。该模型采用焓模型来处理多孔介质内工质的相变问题, 通过分区计算既避免了采用耗时复杂的自适应网格技术又能方便地跟踪汽液界面, 为处理多孔介质中的毛细力带来方便; 在处理流动问题时, 通过对源项的处理考虑了毛细抽力的作用, 体现了毛细抽力是CPL的动力源, 从给定的算例来分析, 该数值模型较好地反映了毛细芯的工作特性。     

一种柔性铰链对铰接板振动特性影响试验研究

针对用于连接小型航天器和太阳能帆板的一种弹簧刚度可调的柔性铰链，研究铰链刚度和铰链转动副的摩擦与铰接结构中阻尼对柔性铰接板系统振动响应的影响.首先，建立了压电柔性铰接柔性板试验平台，采用压电片传感器和激光位移传感器检测，压电片驱动器用于持续激励振动.然后，在柔性铰链扭转刚度不同的情况下，分别进行了外部激励和持续压电驱动器激励的柔性铰接板振动模态分析的实验研究.实验结果分析表明，柔性铰链的刚度和旋转副摩擦阻尼因素，会影响柔性铰接板系统的基频振动固有频率和阻尼比.实验研究结果为分析柔性铰接板的铰链影响特性提供参考.     

罐车制动时液体晃动的仿真分析

采用Volume of Fluid (VOF) 模型对罐车制动时液体的晃动进行了数值模拟,对单室受力、受力位置及整车轴荷分配进行了计算,并通过与实验结果的对比验证了计算方法的可靠性.计算结果表明,无防波板时,随减速度增加,单室 x,y方向受力峰值增大,整车轴荷比增大;随充液比增加,单室x方向受力峰值先增大后减小,y 方向受力峰值增大,制动初始与结束时充液比越大轴荷比越大,1s左右充液比越大轴荷比越小;单室带防波板时,随防波板面积增加, x,y 方向受力峰值减小,当防波板面积大于横截面的40%时,增加防波板面积能显著改善罐体受力,且防波板面积越大轴荷比峰值越小.     

Injection mechanism for cosmic ray acceleration

The quasilinear theory of MHD waves excitation by cosmic rays accelerated at a front of supernova shock has been constructed. It is shown that the energetic particles excite the waves propagating from the shock front, the intensity and the spectrum of these waves is obtained. The role of nonlinear Landau damping in the formation of such spectrum has been analysed. The diffusive scattering length of the high energy particles in the preshock region has been calculated and it is shown with the help of these formulae that the effective Fermi acceleration at the shock front is possible upto the values of the relativistic factor = 10⁴ - 10⁵. The injection mechanism for cosmic rays acceleration has been proposed. It is based on stochastic Fermi acceleration of the thermal plasma by MHD waves excited in the preshock region. Different possibilities for wave phase velocity dispersion needed for stochastic Fermi acceleration are analysed, those are the excitation of the oblique magnetosonic waves as well as the excitation of parallel Alfven waves propagating in opposite directions. The distribution function of the suprathermal particles accelerated by MHD waves is obtained, the cosmic rays density as well as the lower boundary of their energy spectrum realised in the proposed mechanism are also calculated.     

星间洛仑兹力编队飞行的平衡点及零速度曲面

提出了一种新的利用星间洛仑兹力控制卫星相对运动的方法:使主星产生自旋磁场,副星带电,通过控制副星所受的星间洛仑兹力进行编队。假设主星产生的人造磁场表现为偶极子并且运动在一条开普勒圆形轨道上,副星恒定带电并在主星附近运动,同时假设星间洛仑兹力只影响副星的运动而不影响主星的运动。推导了副星在主星HCW(当地垂直当地地平坐标系)坐标系下的相对运动动力学方程。针对偶极子与HCW坐标轴X轴重合的情况下,推导了动力系统下的平衡点(并采用稳定性分析方法分析其线性化意义下的稳定性)、积分常数和零速度曲面,证明了有界相对运动的存在。最后用数值仿真验证了上述结论。     

建立在视线方程基础上的空战动力学模型

首先定义并介绍了视线动力学方程,在此基础上讨论视线坐标系上的机动策略问题。介绍视线坐标系上的空战策略的一种决策方法,以及视线坐标系上的空战策略向机体坐标系上的转化问题。最后建立在机体系上的力和力矩方程,用来说明如何操控飞机达到所需要的空战策略。