三体模型下二维平面地月转移轨道设计与特性分析

为提高地月转移轨道搜索速度和提供精确的初值猜测,在三体模型下研究了双脉冲地月转移轨道的快速设计方法及其轨道特性。首先基于三体模型建立了地心旋转系下地月转移轨道系统模型,选择物理意义明显的地月转移加速速度增量和地心旋转系对准角作为控制变量;其次在二体模型假设下给出控制变量的初值估计方法,提出了不同运行方向地月转移轨道的搜索策略,并采用序列二次规划算法对4种类型地月转移轨道进行求解,分析其控制变量差异;最后分析了转移时间、近地点高度和近月点高度对地月转移轨道的影响,以及月球引力对二体模型下初值估计的修正量。仿真结果表明,本文提出的设计方法能快速准确地搜索出指定类型的地月转移轨道,同时其轨道特性分析结论可为后续月球探测任务地月转移轨道选择提供依据。     

地月L2点空间站转移轨道设计与特性分析

针对低运输成本的未来地月L2点Halo轨道空间站补给任务轨道转移问题,研究了基于不变流形和月球借力的间接转移轨道设计方法及轨道参数特性。首先建立了三脉冲从近地停泊轨道出发到Halo轨道的间接转移轨道设计模型;然后提出了二体解析初值搜索和局部梯度优化相结合的串行优化方法,并采用圆型限制性三体问题模型验证了方法的正确性与快速性;最后,通过大量仿真算例分析了间接转移轨道参数特性。结果表明不变流形插入点相位和Halo轨道幅值对间接转移轨道燃料消耗影响较大,而月球借力高度影响则比较小。     

基于不变流形的地月L2点Halo轨道转移轨道设计

针对Halo轨道转移轨道设计不易收敛的问题,结合月球探测背景,分析了地月L2点Halo轨道及其不变流形可靠近月球的运动轨迹特性,给出基于二分法的先粗选后精选的零消耗转移轨道搜索方法并研究得出Halo轨道全相位点入轨的零消耗转移轨道不超过2条,提出一种最小x轴约束的近月点终止条件和自适应退步搜索的改进微分修正算法,对Halo轨道全相位点入轨的转移轨道设计进行求解。仿真结果表明,该改进微分修正算法收敛速度快,能有效避免奇异,且适应性强,能搜索出全相位点入轨的所有转移轨道。     

地球低轨道空间站往返月球的转移轨道设计

围绕地球低轨道空间站约束下的地月往返转移任务,基于圆锥曲线拼接模型和兰伯特问题提出了该约束下的地月/月地转移轨道解析设计方法,并利用高精度动力学模型对初始轨道设计方法进行了验证;对中国空间站和国际空间站往返于月球的转移轨道进行了特征分析,包括轨道转移窗口、速度增量、转移时间以及月固坐标系下的可达月球低轨道和出发月球低轨道。结果表明:中国空间站或国际空间站与月球低轨道之间的往返轨道转移窗口每月至少存在3个;空间站轨道倾角不影响速度增量与转移时间,且不影响可达月球低轨道或出发月球低轨道在月固坐标系下的轨道倾角-升交点经度分布;利用高精度动力学模型可有效针对轨道面约束下的初始轨道设计结果进行修正。     

基于弹道逃逸和小推力捕获的地月转移轨道设计

针对地月空间货运任务和环月轨道空间设施建设任务,提出一种弹道逃逸和小推力捕获相结合的新型地月轨道转移模式,并建立了一整套该类型轨道设计方法。首先,在三体模型假设下分别建立地心弹道逃逸轨道和月心小推力捕获轨道的二维极坐标动力学模型。对于弹道逃逸轨道,将地心旋转系对准角和地月转移加速速度增量作为控制变量,提出初值估计解析公式,并应用序列二次规划算法进行快速求解。对于小推力捕获轨道,以月心距为参考量设置与弹道逃逸轨道的拼接点约束,提出能量匹配方法预估飞行时间,采用最优螺旋轨道的初始伴随状态解析式预估近月点伴随变量初值。基于混合法和轨道逆推思想,采用人工免疫算法进行小推力捕获轨道求解。仿真结果表明,基于弹道逃逸和小推力捕获的地月轨道转移方式大幅降低了近月制动燃料消耗,能快速穿越地球辐射带,且飞行时间适中;同时,提出的轨道设计方法能快速搜索到基于弹道逃逸和小推力捕获的地月转移轨道,验证了该方法的有效性。     

地月L2点Halo轨道支持的登月轨道优化设计

针对低成本月球探测任务轨道设计问题,基于不变流形提出了地月L2点Halo轨道支持的登月方案,并研究了轨道优化设计方法。建立了3脉冲从近地停泊轨道(LEO)到环月轨道(LLO)的不变流形转移轨道数学模型;设计了全局遍历初值搜索和局部梯度优化的串行优化设计方法,并采用圆型限制性三体模型(CR3BP)优化了一条从LEO出发,通过稳定流形轨道、地月L2点Halo轨道、不稳定流形轨道,最终抵达LLO的燃料最优转移轨道;给出了仿真实例,仿真结果验证了文章所述方法的正确性。研究结论可为未来低成本地月转移或Halo轨道空间站支持的登月方案提供参考。     

月地转移轨道精确轨道设计

以基于Lambert算法的快速轨道设计结果为初值,开展精确轨道设计研究.通过对月地返回飞行阶段的摄动项和量级分析,建立了月地转移轨道的动力学方程,提出了一种双向嵌套循环搜索算法,采用该算法求解同时满足两端约束条件的精确月地转移轨道.该算法以出月球影响球的时刻和位置、速度为中间变量,一方面采用前向数值积分和微分改正法搜索满足地球再入端的轨道,另一方面采用后向数值积分并进行倾角和近月距修正得到满足月球端的轨道,通过这种双向嵌套循环,使得两段轨道在月球影响球边界处的位置和速度连续,从而获得一条完整的满足两端约束条件的月地转移精确轨道.最后以2017年1月26日出月球影响球作为返回窗口,给出了具体的设计算例,并通过STK软件仿真验证了程序的设计结果.     

月球返回轨道再入角变化特征

再入角是航天器返回大气层时在再入点处速度方向与"地平面"之间的夹角。若忽略地球的非球形因素,则可近似的看做轨道切向与横向之间的夹角。为了避免探测器过热问题,一般再入角不宜太大,在3°～8°之间。文章以只在近月点进行一次制动的月球探测器的霍曼转移型的返回轨道为例,通过对轨道性质的分析和数值计算,说明地月相对位置和地球自转对月球返回轨道再入角的影响。分析和计算得到以下结论:1)对于相同的转移时间和固定的再入点,当月球位于南纬最高点时,则再入角的绝对值可以取到最小值;2)对于相同的转移时间和固定的再入角,当月球位于南纬最高点时,再入点的纬度可以取到最大值;3)转移时间越短,再入角的绝对值可以取到更小值,而再入点纬度可以取到更大值。以上这些极值对应的都是极轨轨道。     

月地低能返回轨道最优控制策略研究

为解决月球货运飞船采用的低能返回轨道对初值极为敏感的问题,研究了月地低能返回轨道的最优控制策略。首先基于椭圆四体动力学模型,分析了月地低能返回轨道的动力学特性;进而引入协方差分析法分析了轨道初始飞行状态的误差传播特性,确定了保证终端再入点高度约束要求的飞行器入轨点位置、速度以及入轨时刻控制精度需求;根据轨道对不同时期施加控制的敏感性不同,设计了一种三脉冲轨道控制策略,以实现既精确控制落点约束,又节约控制燃料消耗的目的。从仿真结果可知,该策略可有效控制月地低能返回轨道终端再入点精度,降低初始敏感度。该控制方案用于月地转移可显著降低对推进控制系统的精度需求,提高转移方案的工程可实现性。     

小推力航天器的地月低能转移轨道

 在限制性四体模型下研究基于小推力方式的地月低能转移问题,通过借助于平动点轨道的相空间结构来揭示小推力转移的机理。重点研究了小推力转移自由飞行段的构造：经由LL₁点穿越获得最小能量的低能转移;而经由LL₁点Halo轨道穿越,得到(M,N)圈穿越轨道;由于Halo轨道相对于平动点增加了一维度的选择,根据(2,2)圈穿越轨道构造该转移的自由飞行段。在地球势阱逃逸和月球势阱捕获段,分别设计了合适的小推力的控制律及发动机开/关机时间,成功实施近地球段的小推力加速和近月球段的减速。尽管未对所得到的结果进行优化,所得转移轨道的燃料消耗也与类似边界条件的SMART-1轨道基本一致。     

环形燃烧室两相燃烧数值研究

根据环形燃烧室设计需要,运用计算机模拟技术研究模型环形燃烧室三维有旋液雾燃烧流场,研究两相流动与混合,液雾蒸发对燃烧过程的影响。两相流动模型采用颗粒群轨道模型,气相紊流输运方程求解采用双方程k-ε模型和EBU-Arrhenius紊流燃烧模型,热辐射采用热通量法模型,本文提供计算方法可以用来研究不同性能参数对环形燃烧室两相燃烧特性影响。     

Halo轨道转移及中途修正问题研究(英文)

This article addresses the design of the trajectory transferring from Earth to Halo orbit, and proposes a timing closed-loop strategy of correction maneuver during the transfer in the frame of circular restricted three body problem （CR3BP）. The relation between the Floquet multipliers and the magnitudes of Halo orbit is established, so that the suitable magnitude for the aerospace mission is chosen in terms of the stability of Halo orbit. The stable manifold is investigated from the Poincar6 mapping defined which is different from the previous researches, and six types of single-impulse transfer trajectories are attained from the geometry of the invariant manifolds. Based on one of the trajectories of indirect transfer which are ignored in the most of literatures, the stochastic control theory for imperfect information of the discrete linear stochastic system is applied to design the trajectory correction maneuver. The statistical dispersion analysis is performed by Monte-Carlo simulation,     

相位延迟边界条件在叶轮机械颤振分析中的应用

基于带相位延迟的周期边界条件,建立了某跨声速转子的双通道高效气动阻尼计算模型.数值计算了该转子的气动性能、颤振边界和叶片模态,和实验数据吻合较好.通过传统的多通道能量法以及双通道方法计算了叶片在一弯模态,不同叶片间相位角条件下的气动阻尼,获得了基本一致的计算结果,而双通道方法相比于传统的多通道能量法计算效率提升约7.7倍,内存需求约为后者的0.45倍.不同叶片振幅对气动阻尼结果的影响研究表明,对于较小的叶片振幅,流动非线性对气动阻尼计算结果仍然有显著的影响.不同工况的计算结果表明:叶片间相位角对转子叶片的气动阻尼有显著的影响,对于该转子最小的气动阻尼均在叶片间相位角为-42.4°时得到;同时,在近颤振状态,不同叶片间相位角对应的气动阻尼均小于近设计状态.     

RF航段衔接ILS进近程序的最大下降角设计方法

RNP 程序的RF 航段与传统ILS 精密进近程序相结合能充分利用两类程序的优势,但两种程序的过渡阶段需要考虑最大下降角度的设计。本文分析过渡阶段程序的平面与剖面几何设计参数限制条件,以及从星基导航方式转为陆基无线电导航方式需要稳定截获航向信号与下滑台信号的运行限制条件;结合两种限制条件,设计一种可以从不同最后进近定位点（FAP）高度计算下滑道的偏离量与到最后航向截获点（FACF）的弧距离关系进行判断,从而确定最大下降角度的方法,并进行仿真验证。结果表明：该方法可以在不同过渡航渡长度与不同最后进近下降角度的限制下,根据给定最高程序假设温度计算得到符合两种限制条件的最大下降角,为RNP 程序与ILS 程序衔接过渡的下降剖面设计提供了一种基于理论基础且具有可行性的设计方法。     

环形燃烧室两相喷雾燃烧的大涡模拟

在三维任意曲线坐标系下采用欧拉-拉格朗日两相大涡模拟方法对环形燃烧室火焰筒气液两相紊流瞬态反应流进行数值模拟.采用椭圆偏微分方程生成三维贴体网格,计算中所采用的数学度模型有:k方程亚网格尺度模型估算亚网格紊流粘性;亚网格EBU燃烧模型估算化学反应速率;热通量辐射模型估算辐射换热.并在非交错网格体系下气相采用SIMPLE算法和混合差分格式求解,液相采用随机离散模型(Stachasttc Separated Flow,简称SSF),在拉格朗日坐标系下追踪各油珠群沿各自轨道运动、质量损失及能量变化.通过计算结果与实验数据相比较,表明在三维贴体坐标系下对燃烧室火焰筒两相紊流油雾燃烧流场进行大涡模拟,采用欧拉-拉格朗日两相大涡模拟方法能反映两相紊流化学反应流流动及实际燃烧过程.     

变截面U形通道内肋高对换热特性影响的实验

采用实验方法研究了在定常状态下, 肋的高度对横肋变截面U形通道内的换热特性的影响。通道内肋间距为15 mm, 肋的宽度均为2 mm, 肋的高度分别为1, 1.5, 2 mm和2.5 mm。雷诺数在12 000到56 000之间变化, 采用单元分析法, 以相同形状的无肋通道作为基准, 研究了肋的高度对肋间距为15 mm通道的换热效果及阻力影响, 得到了包括几何参数在内的经验关联式。结果表明, 有肋的通道换热都得到了增强;肋高度的变化对变截面蛇形通道换热的影响不是单调的, 而是呈多峰分布;肋的高度对阻力系数的影响却是随着肋高的增大而增大的。      

多级轴流压气机非设计性能的数值预估

采用LU-SGS隐式解法以及叶排间掺混面模型数值模拟压气机三维粘性流动,研究了多级轴流压气机的设计与非设计状态性能的预估问题。对一台五级高压压气机进行了详细的流场数值计算,包括了设计转速的100%,90%及80%等多个转速及各种工作状态。结果表明,本文所建立的数值方法与计算机程序对多级轴流压气机各种非设计状态的计算都具有收敛容易和快速的优点,计算得到的压气机总性能和各分级性能与设计参数进行了对比,证明该方法和程序的预估结果是比较准确的。根据设计转速下的计算结果,预估了压气机设计点的喘振裕度。此外,还对不同状态下叶片进口气流角的变化进行了分析。      

弹道导弹中段防御的拦截窗口分析

根据弹道导弹中段飞行近似为二体运动的特点，应用二体运动理论分析了弹道导弹中段防御的拦截窗口。假设最低拦截点由拦截弹末段寻的飞行正常工作的最低地面高度决定，最高拦截点由拦截弹投射能力决定，给出了根据当前时刻弹道导弹位置和速度求解拦截窗口的算法。对采用最小能量弹道飞行的弹道导弹，给出了根据射程求解拦截窗口的算法，并对不同射程的最小能量弹道，给出了求解拦截窗口的算例。