CN104657559B - 基于圆柱型速度参数截面的地月自由返回轨道设计方法 - Google Patents
基于圆柱型速度参数截面的地月自由返回轨道设计方法 Download PDFInfo
- Publication number
- CN104657559B CN104657559B CN201510094547.7A CN201510094547A CN104657559B CN 104657559 B CN104657559 B CN 104657559B CN 201510094547 A CN201510094547 A CN 201510094547A CN 104657559 B CN104657559 B CN 104657559B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- moon
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明涉及基于圆柱型速度参数截面的地月自由返回轨道设计方法，属于航空宇航科学与技术的航天器动力学与控制领域。本方法首先给出地月转移轨道圆柱型速度参数截面定义及特点说明，然后基于双曲剩余速度环给出自由返回轨道解存在性分析方法，最后给出了自由返回轨道解求解策略。本发明解析地给出了自由返回轨道解存在性判定条件，方法快速高效；同时，给出了采用解析方法求解初值，高精度动力学模型迭代进行精确轨道计算的求解策略，方法简单可靠、计算速度快，具有重要的应用前景。
Description
技术领域
本发明涉及一种设计地月自由返回轨道的简单解析方法，属于航空宇航科学与技术的航天器轨道动力学与控制领域。
背景技术
地月自由返回轨道是载人登月任务初期，为提高安全性而采用的飞行轨道，载人探月飞行器若采用该类轨道飞行，能在地月转移入轨后，滑行至月球附近，以较低近月点高度绕月，并以一定再入条件返回地球大气，期间不需任何额外轨道控制。传统设计方法包括解析方法和数值方法两大类。
在解析方法中，首先引入月球影响球概念，将地月转移轨道分为地心段和月心段两段二体轨道，即双二体假设，然后通过寻找月球影响球入口点参数来达到设计整条轨道的目的。采用这种方法设计自由返回轨道通常需要提前大量打靶仿真，得到入口点参数解集范围及参数基本特性，该方法常见于任务总体设计决策时使用。该方法中月球影响球拼接参数并不完全具有实际物理意义，双二体假设动力学模型引起的模型误差也比较大。
数值方法是目前为止公认的与实际飞行轨迹最接近的方法，动力学模型中一般采用JPL(Jet Propulsion Laboratory)星历表求解日地月等相对位置关系，采用勒让德多项式准确计算非球形摄动项，一般采用变步长积分器有效提升数值积分速度。由于载人登月自由返回轨道解的高敏感性和解集域的非光滑性，通常采用数值优化算法求解轨道，初值一般通过解析方法求解。这种串行优化求解的方法虽然能有效求解一条轨道，但计算所需时间和内存资源等不适宜多条轨道快速求解。
无论是传统解析方法还是数值方法，在某些复杂约束条件下，都存在载人登月自由返回轨道解不收敛情况，即无法对自由返回轨道解存在性给予明确判断。本发明提出一种基于圆柱型速度参数截面的地月自由返回轨道设计方法，既能给出地月自由返回轨道解存在性条件，也能采用该方法方便求解地月自由返回轨道。目前尚无公开文献涉及此类方法。论述中提到一些名词概念及坐标系等是本领域人员公知的知识，不一一进行详述。
发明内容
本发明的目的是提供一种解析计算载人登月自由返回轨道的简单方法。本方法首先给出地月转移轨道圆柱型速度参数截面定义及特点说明，然后基于双曲剩余速度环给出自由返回轨道解存在性分析方法，最后给出了自由返回轨道解求解策略。
本发明的技术方案如下:
第一步:白道相关坐标系定义
白道是月球绕地心公转的轨道面。为清晰描述本发明方法步骤，先给出与白道面相关的三个坐标系定义：
(1)白道惯性系
包括地心白道系oExEyEzE和月心白道惯性系oLxLyLzL，原点分别取在地心或者月心；坐标系方向定义参考某基准时刻t0；x轴指向月球升交点方向，即图1中PM方向，对应t0记为z轴指向轨道角动量方向，即图1中HM方向，对应t0记为y轴与另外两轴垂直，构成右手系，即图中QM方向，对应t0记为白道面与赤道面夹角记为εM。
(2)白道瞬时坐标系
包括地心白道瞬时系oExyz和月心白道瞬时系oLxyz，原点分别取在地心和月心；坐标系方向定义参考当前时刻t，x轴由地心指向月心方向，即图1中rM方向，位于PM、QM构成的平面内；z轴指向轨道角动量方向，即图1中HM方向；y轴与另外两周垂直，构成右手系。近似分析中可取PM、QM、HM分别与等同。rM方向相对PM的相位角度记为uM。
(3)月心轨道惯性系
在绕月阶段飞行器相对月球从双曲轨道减速为环月轨道，这时的轨道面调整近似为共面调整，可以在环月轨道坐标系内进行分析。月心轨道惯性系oLx1y1z1，原点取在月心；坐标系方向定义参考某基准时刻t0；x1方向指向绕月轨道相对白道的升交点，即P∞方向；z1轴指向绕月轨道角动量方向，即图1、2中H∞方向；y1轴与另外两轴垂直，构成右手系，即图1、2中Q∞反向。P∞与夹角记为Ω∞M，轨道面与白道夹角记为i∞M，如图1、2所示。
第二步:地月转移轨道圆柱形速度参数截面选取
地球出发的地月转移轨道为大椭圆轨道，一般近地点会约束在给定高度，远地点有较大的设计空间，对于给定的月球位置，飞行器到达月球影响球入口点的速度可以在地球惯性系内沿地心径向和横向分解为vr和vτ。月地返回轨道为实现大气捕获，近地点高度会限定在大气层高度范围内，远地点参数为大于月球影响球出口点的值。
下面基于地心二体轨道理论给出仿真分析。图3给出了地月转移轨道远地点-横向速度关系曲线。此算例取近地点高度为200km，入口点地心距33.5万km，横坐标代表远地点相对入口点高度，单位万km；纵坐标为入口点处vτ，单位m/s。图4给出了地月转移轨道远地点-径向速度关系曲线。
此算例取出口点地心距33.5万km；横坐标代表远地点相对出口点高度，单位万km；纵坐标为近地点高度，单位km；等高线标识出口点处vτ，单位m/s。
由图3、4和图5可以看出，对于给定的地月转移问题，当近地点区间有限，远地点有较大设计空间的情况下，月球影响球出口点(入口点)处的横向速度变化不大，可以把它近似为与远地点无关的常数，这样在月球影响球出口点(入口点)处的速度取值空间可以近似为一个圆柱面，本报告将其称为圆柱形速度参数截面。此参数截面实际为近地点固定约束在出口点(入口点)速度参数空间的体现。从严格的意义上此界面为一近似圆柱面回转面。
为便于对不同工况分析，图6给出了出口点(入口点)地心距与横向速度的对应关系。图中，横坐标为出口点(入口点)地心距，单位m；纵坐标为横向速度vτ，单位m/s。实线对应近地点高度200km的情况，虚线对应近地点高度50km的情况。横向速度取值对应远地点高度相对出口点(入口点)40万km的情况。
图7以返回轨道为例给出了出口点近似与月球同向情况下出口点圆柱形速度参数截面示意图。图中，横向速度vτ取为常数，径向速度给定上下限vr1和vr2。
第三步:双曲剩余速度环计算
在oLx1y1z1坐标系中，自由返回绕月轨道为双曲线，相对月球进入速度v∞in和v∞out的大小相同，双曲剩余速度v∞与近月点参数rp、vp关系为
即
c＝ae
双曲线的渐进角δ和焦点到渐近线的距离d可以表示为
近似分析中，在月球影响球处的v∞in、v∞out可以看作就是沿着渐近线方向。这样，当相对月球进入速度v∞in给定，近月点距离rp给定的情况下，所有可能的v∞out矢量分布构成一个以v∞in为轴线的圆锥，圆锥半锥角为2δ，矢量端点构成一个圆环，圆环半径为v∞sin2δ；
设v∞in沿x轴方向，则
其中，θ∞代表圆环上的幅角；
图10-15给出了双曲剩余速度随v∞的变化关系。图10为渐近角变化曲线，单位度；图11为焦点到渐近线的距离变化曲线，单位m；图12为影响球内总飞行时间变化曲线，单位s；图13为v∞＝1200m/s时的轨道高度与渐近角的关系曲线；图14为v∞＝1200m/s时的hp-δ关系曲线，进一步计算可知rp在影响球附近δ≈87°；图15为v∞-Δv关系曲线,Δv为双曲轨道到环月轨道变轨速度。
第四步:自由返回轨道解
在月心白道瞬时系oLxyz中，入口点位置矢量可以表示为rMin，|rMin|＝Rim，Rim为月球影响球半径；在地心白道瞬时系oExyz中，入口点位置rEin＝rM+rMin，入口点方向可以用rEin与oxy平面的夹角βin、和rEin在oxy平面投影与x轴的夹角αin来表征。近似分析中可以取βin＝0、αin＝7°，对应入口点在月球影响球前方45度左右的情况；
在地心白道瞬时系oExyz中，月球速度vM在y轴附近，入口点速度在圆柱形速度参数截面上，由柱面坐标vr，θ可以确定入口点处相对地球的速度
入口点处相对月球的入口速度
v∞in＝vEin-vM (5)
令v∞in的方位角为β∞in和α∞in，则有(3)式中
或者因此有出口速度环方程
类似地，可以给出出口点处的位置关系rEout＝r* M+rMout，r* M表示出口时刻月球的位置，出口点方向角可以定义为βout和αout。出口点处速度关系
v∞out＝vEout-vM (8)
如果忽略绕月双曲轨道中焦点到渐近线的距离d，即取d≈0，这时
v* M相对vM绕z轴转动角度可以表示为θMv＝nMT rp,v∞，T rp,v∞为图12所示的绕月飞行时间，nM为绕月期间月球的平均轨道角速度，θMv也可以直接通过月球星历求取；
自由返回解的存在条件即要求(7)、(8)式同时成立，即要求出口速度环与出口速度圆柱面存在交点。为形象表示，把图16、17和图18、19合并，使其x轴重合，如图20、21所示，注意速度环的轴相对v∞in绕z轴旋转了θMv，v* M和vM近似重合。可以看出，这时出口交点可能有至多2个。空间圆环和柱面的交线一般解难以解析给出，但很方便通过一维搜索算法求解v∞out。v∞out确定后，vEout可以通过(8)式确定；
绕月轨道面法线为
可以求得，入口点和出口点相对(9)式位置线的距离
入口点和出口点位置
由以上过程可以得到圆锥曲线拼接意义下的自由返回轨道近似解析解。
第五步:自由返回轨道解存在性分析
利用上部分解可以对自由返回轨道解的存在性进行分析。定性分析时，不妨设入口点和出口点的圆柱形速度参数截面都沿月球位置方向，月球速度与月球位置矢量垂直，则可以得到图21形式的投影图，这时入口柱面和出口柱面重合。以vM顶端为圆心，v∞为半径作球面，交圆柱于两个空间环(v∞较小时，也可能是一个环)。设vM＝1000m/s，vτ＝200m/s，则可知存在两个相交环时，单环对球心的最大张角为v∞＝1200m/s情况时的48度。由此可以推论，自由返回轨道可以存在以下两大类，每类有可以分为两小类：
(1)出入口反向解
即出口解和入口解vr符号相反的情况。这类解一般v∞in和v∞out有较大的夹角，是比较实际的地月转移轨道。
v∞in和v∞out的夹角极大值和极小值可以利用白道面投影近似分析。夹角近似为
θio＝θMv+π-2δ (13)
图22给出了θio随v∞、hp的分布等高线图。图中，横坐标代表v∞，取值从800m/s到1800m/s变化；纵坐标代表近地点高度hp，取值从0km到50000km变化。可以看出θio随v∞、hp的增大而减小；
因为hp最小为0，而v∞增大时两个相交环之间的距离会拉大，所以自由返回解的v∞存在上限。当v∞＞vM+vτ时，球面与柱面存在两个相交环，两个相交环之间点的最小夹角可以近似表达为图23给出了存在出入口反向解的最大v∞近似分析。图中设vM＝1000m/s，vτ＝200m/s，横坐标为v∞，纵坐标为角度；实线对应两个相交环之间点的最小夹角；虚线对应近地点为月球半径时的最大θio。可以看出，当v∞≈1540m/s时两线相交，即存在出入口反向自由返回解的最大v∞为1540m/s。实际问题中，由于柱面偏斜、月球速度倾角非零等因素，该最大数值需要修正；
实用的自由返回解要求近月点不能离月球过远，这时的θio有较大的取值，由图22、23可看出θio取值在80-130度之间。可以推论，v∞存在最小值，对于前面的参数取值，有最小值约1200m/s。进一步可以推论，这种情况下自由返回轨道的轨道面与白道面夹角小于考虑前面的参数取值，有最大夹角约13度；
出入口反向解又可以分为两大类
(a)升轨进入降轨返回解
这是常见载人登月轨道，具有飞行时间短的特点；
(b)降轨进入升轨返回解
这时的轨道类似心型，具有较长的飞行时间；
(2)出入口同向解
即出口解和入口解符号相同的情况。这类解一般v∞in和v∞out的夹角很小，最大也就对应v∞＝vM+vτ的情况，约48度。这类轨道具有较高的近月点，并不适合用作地月转移轨道，而适合用于辅助轨道面调整等飞行器轨道控制；
球面与柱面相交单边环的张角下图12给出了单边环张角随的变化关系，对照图22可知，v∞过大时单边环张角很小，这时的出入口同向解趋近影响球边界。当v∞＝vM+vτ时，有出入口同向解hp最小值情况，约7000km；
出入口同向解的特点是轨道面的大角度调整乃至反向；
出入口同向解又可以分为两大类：
(a)升轨解
出入口vr都为正；
(b)降轨解
出入口vr都为负。
第六步:自由返回轨道解求解策略
前面的计算模型用于圆锥曲线拼接模型时，圆柱面轴线指向没有显式给出，需要叠代求解，另外v∞的方向在d较大时存在一定的模型误差，因此需要设计数值求解方法，保证解的收敛性。另外以前面模型为出发点，可以应用高精度模型，进行叠代求解；
首先，在任意的时间都存在自由返回解，可以近似认为该解具有随月球公转的不变性。可以设进入月球影响球的时刻为初始时刻t0，这时对应有月球位置rM t0、速度vM t0可以由月球星历得到。
以入口点的速度为基本设计参数，在圆柱形速度参数截面上描述参数为柱面坐标vrin,θin，横向速度vτin初始值可以采用图6计算给出，初步设定αin、βin，从而可由(5)求得v∞in；
对于给定的近月点高度hp，可以求得δ、T、θMv。通过图20、21中的分析，可以数值求得圆环和出口速度柱面轴线的最短距离；若最短距离小于vτout，则说明存在两个交点，分别数值求出；若最短距离大于vτout则说明无解，为防止漏解，在大于vτout不多的情况，可以取最小距离点解算，叠代后该点如果仍然柱面外则证实无解。根据解得的v∞out，可以确定rMin、rMout，从而得到rEin、rEout，从而类似前面叠代求解出各参数。在这种策略中，t0、vrin、θin、hp四个参数为给定值，其他参数都需要叠代求解；
高精度轨道求解时，vτin、rEin的特性和前面一致，这时对din的切向和径向调整目标可以选择hp和vτout，其他参数通过积分导出即可。
本发明的优点在于：
1、以实际物理量描述了自由返回轨道求解过程，并明确给出自由返回轨道存在性判据；
2、通过简单解析方式计算自由返回轨道，简单可靠且计算速度快。
附图说明
图1为白道与赤道关系图；
图2为白道与环月轨道面关系图；
图3为地月转移轨道远地点-横向速度关系图；
图4为地月转移轨道远地点-径向速度关系图；
图5为月地转移轨道远地点-近地点-横向速度关系图；
图6为地心距与横向速度的对应关系图；
图7返回轨道圆柱形速度参数截面图；
图8-9双曲剩余速度锥面图；
图10-15双曲剩余速度特性图；
图16为入口位置关系图；
图17为入口速度关系图；
图18为出口位置关系图；
图19为出口速度关系图；
图20为圆环与柱面图；
图21为解存在关系图；
图22为v∞in、v∞out夹角随v∞、hp分布图；
图23为存在出入口反向解的最大v∞图；
图24为单边环的张角图。
具体实施方式
本发明的具体实施流程如下：
第一步:根据实际工程背景给出约束条件；
第二步:采用圆柱型速度参数截面方法判定自由返回轨道存在性；
第三步:采用圆柱型速度参数截面方法快速求解自由返回轨道。
采用本发明可以达到以下技术效果：
1、本发明较传统月球影响球参数拼接法，选取描述轨道的参数具有实际物理意义；
2、本发明明确给出了自由返回轨道存在性判定条件，第一次理论上阐述了以往求解方法不收敛的真实原因；
3、本发明不需要进行非线性方程的求解、数值积分或插值，方法简单、可靠、计算速度快。
总之，本发明具有简单、可靠、精度较高和计算速度快等的优点，具有较好的应用前景。
Claims (1)
1.基于圆柱型速度参数截面的地月自由返回轨道设计方法，首先给出地月转移轨道圆柱型速度参数截面定义，然后基于双曲剩余速度环给出自由返回轨道解存在性分析方法，最后给出自由返回轨道解求解方法，其特征在于，该方法包括如下步骤：
第一步:白道相关坐标系定义
白道是月球绕地心公转的轨道面，先给出与白道面相关的三个坐标系定义：
(1)白道惯性系
包括地心白道系oExEyEzE和月心白道惯性系oLxLyLzL，地心白道系原点取在地心，月心白道惯性系原点取在月心；坐标系方向定义参考某基准时刻t0；x轴指向月球升交点方向，即PM方向，对应t0记为z轴指向轨道角动量方向，即HM方向，对应t0记为y轴与另外两轴垂直，构成右手系，即QM方向，对应t0记为白道面与赤道面夹角记为εM；
(2)白道瞬时坐标系
包括地心白道瞬时系oExyz和月心白道瞬时系oLxyz，原点分别取在地心和月心；坐标系方向定义参考当前时刻t，x轴由地心指向月心方向，rM方向，位于PM、QM构成的平面内；z轴指向轨道角动量方向，HM方向；y轴与另外两轴垂直，构成右手系，近似分析中可取PM、QM、HM分别与等同，rM方向相对PM的相位角度记为uM；
(3)月心轨道惯性系
在绕月阶段飞行器相对月球从双曲轨道减速为环月轨道，这时的轨道面调整近似为共面调整，在环月轨道坐标系内进行分析，月心轨道惯性系oLx1y1z1，原点取在月心；坐标系方向定义参考基准时刻t0；x1方向指向绕月轨道相对白道的升交点，即P∞方向；z1轴指向绕月轨道角动量方向，即H∞方向；y1轴与另外两轴垂直，构成右手系，即Q∞方向，P∞与夹角记为Ω∞M，轨道面与白道夹角记为i∞M，
第二步:地月转移轨道圆柱形速度参数截面选取
地球出发的地月转移轨道为大椭圆轨道，对于给定的月球位置，飞行器到达月球影响球入口点的速度在地球惯性系内沿地心径向和横向分解为vr和vτ，月地返回轨道为实现大气捕获，近地点高度会限定在大气层高度范围内，远地点参数为大于月球影响球出口点的值；
对于给定的地月转移问题，当近地点区间有限，远地点有较大设计空间的情况下，月球影响球出口点或入口点处的横向速度变化不大，把它近似为与远地点无关的常数，这样在月球影响球出口点或入口点处的速度取值空间近似为一个圆柱面，将其称为圆柱形速度参数截面，此参数截面实际为近地点固定约束在出口点或入口点速度参数空间的体现，
第三步:双曲剩余速度环计算
在oLx1y1z1坐标系中，自由返回绕月轨道为双曲线，相对月球进入速度v∞in和v∞out的大小相同，双曲剩余速度v∞与近月点参数rp、vp关系为
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>&infin;</mi>
</msub>
<mo>=</mo>
<msqrt>
<mfrac>
<msub>
<mi>&mu;</mi>
<mi>M</mi>
</msub>
<mi>a</mi>
</mfrac>
</msqrt>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>a</mi>
<mo>=</mo>
<mfrac>
<mi>p</mi>
<mrow>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>e</mi>
<mo>=</mo>
<mfrac>
<mi>p</mi>
<msub>
<mi>r</mi>
<mi>p</mi>
</msub>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>r</mi>
<mi>p</mi>
</msub>
<msup>
<msub>
<mi>v</mi>
<mi>p</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>M</mi>
</msub>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
即
c＝ae
双曲线的渐进角δ和焦点到渐近线的距离d可以表示为
d＝ccosδ (2)
在月球影响球处的v∞in、v∞out看作是沿着渐近线方向，这样，当相对月球进入速度v∞in给定，近月点距离rp给定的情况下，所有可能的v∞out矢量分布构成一个以v∞in为轴线的圆锥，圆锥半锥角为2δ，矢量端点构成一个圆环，圆环半径为v∞sin2δ，
设v∞in沿x轴方向，则
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>v</mi>
<mi>&infin;</mi>
</msub>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>&infin;</mi>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>v</mi>
<mi>&infin;</mi>
</msub>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>&infin;</mi>
</msub>
<mo>)</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mi>&pi;</mi>
<mo>-</mo>
<mn>2</mn>
<mi>&delta;</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>&infin;</mi>
</msub>
<mo>)</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mo>(</mo>
<mi>&pi;</mi>
<mo>-</mo>
<mn>2</mn>
<mi>&delta;</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>sin</mi>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>&infin;</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，θ∞代表圆环上的幅角，
第四步:自由返回轨道解
在月心白道瞬时系oLxyz中，入口点位置矢量可以表示为rMin，|rMin|＝Rim，Rim为月球影响球半径；在地心白道瞬时系oExyz中，入口点位置rEin＝rM+rMin，入口点方向可以用rEin与oxy平面的夹角βin、和rEin在oxy平面投影与x轴的夹角αin来表征，近似分析中可以取βin＝0、αin＝7°，对应入口点在月球影响球前方45度的情况；
在地心白道瞬时系oExyz中，月球速度vM在y轴附近，入口点速度在圆柱形速度参数截面上，由柱面坐标vr,θ可以确定入口点处相对地球的速度
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>E</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<msub>
<mi>v</mi>
<mi>r</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>&tau;</mi>
</msub>
<mi>cos</mi>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>&tau;</mi>
</msub>
<mi>sin</mi>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<msub>
<mi>v</mi>
<mi>r</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>&tau;</mi>
</msub>
<mi>cos</mi>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>&tau;</mi>
</msub>
<mi>sin</mi>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
入口点处相对月球的入口速度
v∞in＝vEin-vM (5)
令v∞in的方位角为β∞in和α∞in，则有(3)式中
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>v</mi>
<mi>&infin;</mi>
</msub>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
或者因此有出口速度环方程
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>&infin;</mi>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
给出出口点处的位置关系rEout＝r* M+rMout，r* M表示出口时刻月球的位置，出口点方向角可以定义为βout和αout，出口点处速度关系
v∞out＝vEout-v* M (8)
忽略绕月双曲轨道中焦点到渐近线的距离d，即取d≈0，这时
<mrow>
<mfrac>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mi>&infin;</mi>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<msub>
<mi>r</mi>
<mrow>
<mi>M</mi>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
</mfrac>
<mo>,</mo>
<mfrac>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mi>&infin;</mi>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<msub>
<mi>r</mi>
<mrow>
<mi>M</mi>
<mi>i</mi>
<mi>n</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
v* M相对vM绕z轴转动角度可以表示为θMv＝nMT(rp,v∞)，T(rp,v∞)为绕月飞行时间，nM为绕月期间月球的平均轨道角速度，θMv也可以直接通过月球星历求取，
自由返回解的存在条件即要求(7)、(8)式同时成立，即要求出口速度环与出口速度圆柱面存在交点，速度环的轴相对v∞in绕z轴旋转了θMv，v* M和vM近似重合，这时出口交点有至多2个，通过一维搜索算法求解v∞out，v∞out确定后，vEout通过(8)式确定，
绕月轨道面法线为
<mrow>
<msub>
<mi>H</mi>
<mi>&infin;</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
求得，入口点和出口点相对(9)式位置线的距离
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>H</mi>
<mi>&infin;</mi>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<mi>d</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>H</mi>
<mi>&infin;</mi>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<mi>d</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
入口点和出口点位置
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mi>M</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mi>&infin;</mi>
</msub>
</mfrac>
<msqrt>
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>+</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mi>M</mi>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>v</mi>
<mrow>
<mi>&infin;</mi>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mi>&infin;</mi>
</msub>
</mfrac>
<msqrt>
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>+</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
由以上过程得到圆锥曲线拼接意义下的自由返回轨道近似解析解；
第五步:自由返回轨道解存在性分析
利用上部分解对自由返回轨道解的存在性进行分析，定性分析时，设入口点和出口点的圆柱形速度参数截面都沿月球位置方向，月球速度与月球位置矢量垂直，得到投影图，这时入口柱面和出口柱面重合，以vM顶端为圆心，v∞为半径作球面，交圆柱于两个空间环，设vM＝1000m/s，vτ＝200m/s，则知存在两个相交环时，单环对球心的最大张角为v∞＝1200m/s情况时的48度，由此推论，自由返回轨道存在以下两大类，每类分为两小类：
(1)出入口反向解
即出口解和入口解vr符号相反的情况，这类解v∞in和v∞out有较大的夹角，是实际的地月转移轨道，
v∞in和v∞out的夹角极大值和极小值利用白道面投影近似分析，夹角近似为
θio＝θMv+π-2δ (13)
看出θio随v∞、hp的增大而减小；
当v∞＞vM+vτ时，球面与柱面存在两个相交环，两个相交环之间点的最小夹角近似表达为设vM＝1000m/s，vτ＝200m/s，横坐标为v∞，纵坐标为角度；当v∞≈1540m/s时，即存在出入口反向自由返回解的最大v∞为1540m/s，
推论，v∞存在最小值，对于前面的参数取值，有最小值1200m/s，进一步推论，这种情况下自由返回轨道的轨道面与白道面夹角小于考虑前面的参数取值，有最大夹角13度，
出入口反向解又分为两大类
(a)升轨进入降轨返回解
常见载人登月轨道，
(b)降轨进入升轨返回解
轨道类似心型，
(2)出入口同向解
即出口解和入口解符号相同的情况，这类解v∞in和v∞out的夹角很小，最大也就对应v∞＝vM+vτ的情况，48度，
球面与柱面相交单边环的张角v∞过大时单边环张角很小，这时的出入口同向解趋近影响球边界，当v∞＝vM+vτ时，有出入口同向解hp最小值情况，7000km；
出入口同向解的特点是轨道面的大角度调整乃至反向；
出入口同向解又分为两大类：
(a)升轨解
出入口vr都为正，
(b)降轨解
出入口vr都为负，
第六步:自由返回轨道解求解方法
首先，在任意的时间都存在自由返回解，近似认为该解具有随月球公转的不变性，设进入月球影响球的时刻为初始时刻t0，这时对应有月球位置rM t0、速度vM t0由月球星历得到，
以入口点的速度为基本设计参数，在圆柱形速度参数截面上描述参数为柱面坐标vrin,θin，横向速度vτin初始值计算给出，初步设定αin、βin，由公式(5)求得v∞in；
对于给定的近月点高度hp，求得δ、T、θMv，求得圆环和出口速度柱面轴线的最短距离；若最短距离小于vτout，则说明存在两个交点，分别求出数值；若最短距离大于vτout，则说明无解，为防止漏解，在大于vτout的情况，取最小距离点解算，叠代后该点如果仍然在柱面外则证实无解，根据解得的v∞out，确定rMin、rMout，从而得到rEin、rEout，从而类似前面叠代求解出各参数，t0、vrin、θin、hp四个参数为给定值，其他参数都需要叠代求解；
高精度轨道求解时，vτin、rEin的特性和前面一致，这时对din的切向和径向调整目标选择hp和vτout，其他参数通过积分导出即可。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510094547.7A CN104657559B (zh) | 2015-03-04 | 2015-03-04 | 基于圆柱型速度参数截面的地月自由返回轨道设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510094547.7A CN104657559B (zh) | 2015-03-04 | 2015-03-04 | 基于圆柱型速度参数截面的地月自由返回轨道设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104657559A CN104657559A (zh) | 2015-05-27 |
| CN104657559B true CN104657559B (zh) | 2018-02-23 |
Family
ID=53248679
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510094547.7A Active CN104657559B (zh) | 2015-03-04 | 2015-03-04 | 基于圆柱型速度参数截面的地月自由返回轨道设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104657559B (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105574261B (zh) * | 2015-12-15 | 2018-09-21 | 北京理工大学 | 一种月球借力约束的地月平动点转移轨道设计方法 |
| CN105631095B (zh) * | 2015-12-18 | 2017-02-15 | 中国人民解放军国防科学技术大学 | 一种等间隔发射的多约束地月转移轨道簇搜索方法 |
| CN105718659B (zh) * | 2016-01-21 | 2019-01-18 | 西北工业大学 | 一种高面质比航天器轨道动力学分析方法 |
| CN107609267B (zh) * | 2017-09-12 | 2019-07-16 | 北京理工大学 | 一种月球有限推力多次捕获轨道实现方法 |
| CN110704952B (zh) * | 2019-09-30 | 2022-09-09 | 中国人民解放军国防科技大学 | 一种月地三脉冲返回轨道速度增量分析方法 |
| CN110909461B (zh) * | 2019-11-13 | 2021-09-17 | 清华大学 | 基于可达集概念的地月/月地直接转移轨道设计方法 |
| CN111324954B (zh) * | 2020-02-18 | 2022-02-11 | 西安交通大学 | 一种月基探测平台对地可见性的判别方法 |
| CN113310496B (zh) * | 2021-05-08 | 2024-01-09 | 北京航天飞行控制中心 | 一种确定月地转移轨道的方法及装置 |
| CN113978768B (zh) * | 2021-10-28 | 2024-01-05 | 中国西安卫星测控中心 | 利用月球引力辅助部署Retro-GEO轨道设计方法 |
| CN116626386B (zh) * | 2023-05-11 | 2023-12-15 | 中国科学院国家空间科学中心 | 一种基于地月锥时空剖分的频谱信息表征方法及系统 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103728980A (zh) * | 2014-01-08 | 2014-04-16 | 哈尔滨工业大学 | 航天器相对轨道的控制方法 |
-
2015
- 2015-03-04 CN CN201510094547.7A patent/CN104657559B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103728980A (zh) * | 2014-01-08 | 2014-04-16 | 哈尔滨工业大学 | 航天器相对轨道的控制方法 |
Non-Patent Citations (4)
| Title |
|---|
| Launch window for manned Moon-to-Earth trajectories;Jingyang Li 等;《Aircraft Engineering and Aerospace Technology》;20121231;第84卷(第5期);第344-356页 * |
| 地月转移轨道设计的改进微分校正方法;贺波勇 等;《国防科技大学学报》;20141231;第36卷(第6期);第60-64页 * |
| 载人登月地月转移轨道快速设计及特性分析;彭祺擘 等;《载人航天》;20120530;第18卷(第3期);第18-24页 * |
| 载人登月转移轨道偏差传播机理分析与稳健性设计;贺波勇 等;《物理学报》;20131231;第62卷(第19期);第190505-1-8页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104657559A (zh) | 2015-05-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104657559B (zh) | 基于圆柱型速度参数截面的地月自由返回轨道设计方法 | |
| CN104035335B (zh) | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导方法 | |
| CN105631095B (zh) | 一种等间隔发射的多约束地月转移轨道簇搜索方法 | |
| CN109344449B (zh) | 航天器月地转移轨道逆向设计方法 | |
| CN103728980A (zh) | 航天器相对轨道的控制方法 | |
| Zhang et al. | Optimal two-impulse rendezvous using constrained multiple-revolution Lambert solutions | |
| CN104076819B (zh) | 一种圆参考轨道下卫星有界伴飞的边界控制方法 | |
| CN105574261B (zh) | 一种月球借力约束的地月平动点转移轨道设计方法 | |
| CN106383994A (zh) | 一种基于脉冲和气动辅助结合的低轨轨道面转移方法 | |
| CN106021784B (zh) | 一种基于两层优化策略的全轨迹优化设计方法 | |
| Wolf et al. | Systems for pinpoint landing at Mars | |
| CN105511493B (zh) | 一种基于火星大气辅助的低轨星座部署方法 | |
| CN108387140A (zh) | 一种考虑多个禁飞区约束的解析再入制导方法 | |
| Wolf et al. | Performance trades for Mars pinpoint landing | |
| CN104484493B (zh) | 一种飞船返回预定落点回归轨道设计方法 | |
| CN104019818B (zh) | 一种基于预测轨迹的行星导航轨道器布局优化方法 | |
| CN105573337A (zh) | 一种满足再入角和航程约束的离轨制动闭路制导方法 | |
| Wolf et al. | Toward improved landing precision on Mars | |
| CN104309822A (zh) | 一种基于参数优化的航天器单脉冲水滴形绕飞轨迹悬停控制方法 | |
| Valenzuela et al. | Conflict resolution in converging air traffic using trajectory patterns | |
| CN108082538B (zh) | 一种考虑始末约束的多体系统低能量捕获轨道方法 | |
| CN114019995B (zh) | 一种燃料最优的异面多目标快速飞掠轨迹规划方法 | |
| CN105718660A (zh) | 临近空间大范围机动弹道三维包络计算方法 | |
| Shen et al. | A 6DOF mathematical model of parachute in Mars EDL | |
| CN113093246A (zh) | 地面多目标点成像快速判定及任务参数计算方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |