CN110058603A - 一种形变再入飞行器变形前指令力矩确定方法 - Google Patents
一种形变再入飞行器变形前指令力矩确定方法 Download PDFInfo
- Publication number
- CN110058603A CN110058603A CN201910276877.6A CN201910276877A CN110058603A CN 110058603 A CN110058603 A CN 110058603A CN 201910276877 A CN201910276877 A CN 201910276877A CN 110058603 A CN110058603 A CN 110058603A
- Authority
- CN
- China
- Prior art keywords
- deformation
- axis
- instruction
- aircraft
- torque
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开一种形变再入飞行器变形前指令力矩确定方法,包括:获取飞行器形变后的指令滚转角速度;预估飞行器形变前的转动惯量;确定飞行器形变前的指令角速度;基于导航姿态角和指令姿态角计算误差姿态角;计算飞行器形变前调姿指令力矩;按照姿控推力器开关逻辑开机进行姿态控制,使飞行器指向指令姿态,以满足控制精度要求;及获取飞行器变形前的启旋指令力矩。本发明中通过预估飞行器形变前的转动惯量,并计算误差姿态角和形变前的调姿指令力矩,按照姿态控制推力器开关逻辑开机进行姿态控制,使得形变再入飞行器指向指令姿态达到预设的控制精度要求,能够满足形变再入飞行器形变后角速度的要求。
Description
技术领域
本发明涉及指令力矩确定方法。更具体地,涉及一种形变再入飞行器变形前指令力矩确定方法。
背景技术
在地球大气层以外的宇宙空间,基本上按照天体力学的规律运行的各类飞行器,又称空间飞行器。对于空间飞行器来说,其制导、导航与控制分系统,简称GNC分系统,GNC分系统承担着空间飞行器从起飞到返回的全部运动控制任务。
航天器再入是指完成预定空间任务后,通过制动进入大气层,然后在地面上着陆的过程,包括返回式卫星、飞船、空天飞行器等。其中,形变再入飞行器是一种再入前通过变形来增大迎风面积,提高阻力系数的新型飞行器。该飞行器在再入前,将姿态调整到指令姿态后,通过安装在其尾部的姿控推力器使其旋转到一定的角速度,然后变形,确保其以一定的自旋角速度再入大气层,这样可以使其在大气层内飞行时保持姿态稳定,并可减少再入过程中多种因素造成的落点散布。
传统的再入飞行器往往是刚性飞行器,其指令力矩只需要将指令角速度直接引入姿态控制回路即可计算。然而,对于形变再入飞行器,由于其转动惯量在形变前后存在很大的变化,根据动量守恒定律,其角速度也会发生较大的变化。因此,传统的指令力矩确定方法无法满足形变再入飞行器形变后的角速度要求,在使用上具有一定的局限性。
因此,需要提供一种形变再入飞行器变性前指令力矩确定方法。
发明内容
本发明的目的在于提供一种形变再入飞行器变性前指令力矩确定方法。
为达到上述目的,本发明采用下述技术方案:
一种形变再入飞行器变形前指令力矩确定方法,该方法包括:
获取形变再入飞行器形变后的指令滚转角速度;
预估形变再入飞行器形变前的转动惯量;
确定形变再入飞行器形变前的指令角速度;
基于导航姿态角和指令姿态角计算误差姿态角;
计算形变再入飞行器形变前调姿指令力矩;
按照姿控推力器开关逻辑开机进行姿态控制,使飞行器指向指令姿态,以满足控制精度要求;及
获取形变再入飞行器变形前的启旋指令力矩。
可选的,获取形变再入飞行器形变后的指令滚转角速度包括:
获取形变再入飞行器形变后的最佳滚转角速度;及
将最佳滚转角速度作为指令滚转角速度P_ωr并上注给形变再入飞行器。
可选的,根据以下公式预估形变再入飞行器变形前的转动惯量Jx:
其中,转动惯量系统坐标系包括X轴、Y轴和Z轴;Jx为飞行器形变前绕X轴的转动惯量;Jxempty为飞行器形变前空载绕X轴的转动惯量;Jxfull为飞行器形变前满载绕X轴的转动惯量;mempty为飞行器空载状态下质量;mfull为飞行器满载状态下质量;m为飞行器当前质量。
可选的,根据以下公式确定形变再入飞行器形变前的指令角速度ωxr:
其中,ωxr为飞行器形变前的绕X轴的指令角速度;Jxopen为飞行器形变后下绕X轴的转动惯量。
可选的,根据以下公式计算误差姿态角:
θe=θr-θ
φe=φr-φ
γe=γr-γ
其中,θe为误差俯仰角,φe为误差偏航角,γe为误差滚转角;θr为指令俯仰角,φr为指令偏航角,γr为指令滚转角;θ为导航俯仰角,φ为导航偏航角,γ为导航滚转角。
可选的,根据以下公式计算形变再入飞行器形变前调姿指令力矩:
Mx0=KPIx·(γe+∫γedt)+kdx·(ωx-ωxr)
My0=KPIy·(φe+∫φedt)+kdy·ωy
Mz0=KPIz·(θe+∫θedt)+kdz·ωz
其中,Mx0为X轴的形变前调姿指令力矩,My0为Y轴的形变前调姿指令力矩,Mz0为Z轴的形变前调姿指令力矩;KPIx为X轴的比例项和积分项系数,KPIy为Y轴的比例项和积分项系数,KPIz为Z轴的比例项和积分项系数;kdx为X轴的阻尼系数,kdy为Y轴的阻尼系数,kdz为Z轴的阻尼系数;ωx为导航得到的角速度在飞行器的X轴的分量,ωy为导航得到的角速度在飞行器的Y轴的分量,ωz为导航得到的角速度在飞行器的Z轴的分量。
可选的,在调姿的过程中,若连续N帧满足以下公式,则判断调姿到位:
(|θe|<Lev)&&(|φe|<Lev)&&(|γe|<Lev)
其中,Lev为姿控精度门限。
可选的,获取形变再入飞行器变形前的启旋指令力矩包括:
令X轴、Y轴、Z轴的比例项和积分项系数为零;
按照以下公式计算启旋指令力矩:
Mx1=kdx·(ωx-ωxr)
My1=kdy·ωy
Mz1=kdz·ωz
其中,Mx1为X轴的启旋指令力矩;My1为Y轴的启旋指令力矩;Mz1为Z轴的启旋指令力矩。
本发明的有益效果如下:
本发明的一种形变再入飞行器变性前指令力矩确定方法,通过预估形变再入飞行器形变前的转动惯量,并通过计算误差姿态角和形变前的调姿指令力矩,按照姿态控制推力器开关逻辑开机进行姿态控制,使得形变再入飞行器指向指令姿态达到预设的控制精度要求,能够满足形变再入飞行器形变后角速度的要求。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明。
图1示出本发明中形变再入飞行器变形前指令力矩确定方法的流程图。
图2示出本发明实施例中姿态角的仿真结果。
图3示出本发明实施例中角速度的仿真结果。
图4示出本发明实施例中指令力矩的仿真结果。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
应注意的是,本发明的说明书和权利要求书及上述附图中的属于“第一”、“第二”等是用于区别不同的对象,而不是用于描述特定顺序。此外,术语“包括”和“具有”以及它们任何变形,意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元,而是可选地还包括没有列出的步骤或单元,或可选地还包括对于这些过程、方法或设备固有的气体步骤或单元。
本发明中公开一种形变再入飞行器变形前指令力矩确定方法,该方法包括:获取形变再入飞行器形变后的指令滚转角速度;预估形变再入飞行器形变前的转动惯量;确定形变再入飞行器形变前的指令角速度;基于导航姿态角和指令姿态角计算误差姿态角;计算形变再入飞行器形变前调姿指令力矩;按照姿控推力器开关逻辑开机进行姿态控制,使飞行器指向指令姿态,以满足控制精度要求;及获取形变再入飞行器变形前的启旋指令力矩。
本发明的一种形变再入飞行器变性前指令力矩确定方法,通过预估形变再入飞行器形变前的转动惯量,并通过计算误差姿态角和形变前的调姿指令力矩,按照姿态控制推力器开关逻辑开机进行姿态控制,使得形变再入飞行器指向指令姿态达到预设的控制精度要求,能够满足形变再入飞行器形变后角速度的要求。
下面结合一个具体实施例进行说明,如图1所示
第一步,获取形变后的指令滚转角速度P_ωr。
地面技术人员在飞行器返回前,通过技术分析,获得飞行器形变后的最佳滚转角速度,采用该滚转角速度再入后,飞行器能够最大限度地保持姿态稳定,克服大气多种因素的干扰,落点精度最高。
将该最佳角速度作为指令P_ωr,上注给飞行器。
在本实施例中,形变后的指令滚转角速度P_ωr=30°/s。
第二步,估计飞行器形变前的转动惯量。
飞行器形变前的转动惯量无法直接测量,然而,可以根据质量与转动惯量的对应关系进行估计。一般来说,可以根据以下公式进行估算:
式中:
Jx:飞行器形变前绕X轴的转动惯量;
Jxempty:飞行器形变前空载绕X轴的转动惯量;
Jxfull:飞行器形变前满载绕X轴的转动惯量;
mempty:飞行器空载状态下质量;
mfull:飞行器满载状态下质量;
m:飞行器当前质量。
在本实施例中,飞行器形变前空载绕X轴的转动惯量Jxempty=25kg·m2,飞行器形变前满载绕X轴的转动惯量Jxfull=30kg·m2,飞行器空载状态下质量mempty=350kg,飞行器当前质量m=300kg,经计算,飞行器形变前绕X轴的转动惯量Jx=25kg·m2。
第三步,确定飞行器形变前的指令角速度。
利用下面公式计算
其中:
ωxr——飞行器形变前的绕X轴的指令角速度;
Jxopen——飞行器形变后绕X轴的转动惯量。
在本实施例中,飞行器形变后绕X轴的转动惯量Jxopen=100kg·m2,则经计算飞行器形变前的绕X轴的指令角速度ωxr=120°/s。
第四步,根据导航姿态角和指令姿态角,计算误差姿态角。
θe=θr-θ
φe=φr-φ
γe=γr-γ
其中:
θe,φe,γe——误差俯仰角、误差偏航角、误差滚转角;
θr,φr,γr——指令俯仰角、指令偏航角、指令滚转角;
θ,φ,γ——导航俯仰角、导航偏航角、导航滚转角;
第五步,计算形变前调姿指令力矩。
利用以下公式计算形变前调姿指令力矩。
Mx0=KPIx·(γe+∫γedt)+kdx·(ωx-ωxr)
My0=KPIy·(φe+∫φedt)+kdy·ωy
Mz0=KPIz·(θe+∫θedt)+kdz·ωz
其中:
Mx0、My0、Mz0——形变前调姿指令力矩;
KPIx、KPIy、KPIz——X轴、Y轴、Z轴的比例项和积分项系数;
kdx、kdy、kdz——X轴、Y轴、Z轴的阻尼系数;
ωx、ωy、ωz——导航得到的角速度在飞行器的X轴、Y轴、Z轴的分量。
第六步,按照姿控推力器开关逻辑开机进行姿态控制,使飞行器指向指令姿态,直到达到控制精度要求。
在调姿的过程中,若连续N帧满足以下公式,则认为调姿到位。
(|θe|<Lev)&&(|φe|<Lev)&&(|γe|<Lev)
式中,Lev为姿控精度门限。
在本实施例中,姿控精度门限Lev=0.5°,N=10。
应说明的是,推力器的具体开关逻辑与推力器的具体布局相关,不具有通用性,在此不做研究。
第七步,令X轴、Y轴、Z轴的比例项和积分项系数为零,计算启旋指令力矩。
令KPIx、KPIy、KPIz均为零后,第五步的指令力矩演变为:
Mx1=kdx·(ωx-ωxr)
My1=kdy·ωy
Mz1=kdz·ωz
Mx1、My1、Mz1——启旋指令力矩。
至此,完成了形变再入飞行器变形前指令力矩确定方法。
下面通过仿真验证该方法的可行性
假设某飞行器形变前后的转动惯量分别为:Jx=25kg·m2,Jxopen=100kg·m2,为了确保再入过程姿态稳定,要求形变后的角速度大于30°/s。
形变前的指令姿态分别为θr=-55.25°,φr=-162.1°,γr=64.68°。
仿真结果如下,其中图2为姿态角的仿真结果,图3为角速度的仿真结果,图4为指令力矩的仿真结果。
由以上仿真结果可知,0~46s为飞行器形变前调姿段,按照本申请的第一到第六步完成姿态调整,46s~50s为启旋阶段,按照本申请的第七步确定的指令力矩进行控制,将X轴角速度提升至121°/s。55s~66s为形变阶段,在形变过程中,随着转动惯量的逐渐增大,X轴的角速度逐渐下降,最终,角速度下降到32°/s,满足要求的30°/s的要求。
由以上仿真结果及分析可知,本发明中的形变再入飞行器形变前指令力矩的确定方法能够满足形变再入飞行器形变后角速度的要求,具有可行性。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以对本发明进行若干改进和修饰,这些改进和修饰也落入本发明权利要求的保护范围内。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (8)
1.一种形变再入飞行器变形前指令力矩确定方法,其特征在于,所述方法包括:
获取所述形变再入飞行器形变后的指令滚转角速度;
预估所述形变再入飞行器形变前的转动惯量;
确定所述形变再入飞行器形变前的指令角速度;
基于导航姿态角和指令姿态角计算误差姿态角;
计算所述形变再入飞行器形变前调姿指令力矩;
按照姿控推力器开关逻辑开机进行姿态控制,使飞行器指向指令姿态,以满足控制精度要求;及
获取所述形变再入飞行器变形前的启旋指令力矩。
2.根据权利要求1所述的指令力矩确定方法,其特征在于,所述获取所述形变再入飞行器形变后的指令滚转角速度包括:
获取所述形变再入飞行器形变后的最佳滚转角速度;及
将所述最佳滚转角速度作为指令滚转角速度P_ωr并上注给所述形变再入飞行器。
3.根据权利要求2所述的指令力矩确定方法,其特征在于,所述根据以下公式预估所述形变再入飞行器变形前的转动惯量Jx:
其中,
所述转动惯量系统坐标系包括X轴、Y轴和Z轴;
Jx为飞行器形变前绕X轴的转动惯量;
Jxempty为飞行器形变前空载绕X轴的转动惯量;
Jxfull为飞行器形变前满载绕X轴的转动惯量;
mempty为飞行器空载状态下质量;
mfull为飞行器满载状态下质量;
m为飞行器当前质量。
4.根据权利要求3所述的指令力矩确定方法,其特征在于,根据以下公式确定所述形变再入飞行器形变前的指令角速度ωxr:
其中,
ωxr为飞行器形变前的绕X轴的指令角速度;
Jxopen为飞行器形变后下绕X轴的转动惯量。
5.根据权利要求4所述的指令力矩确定方法,其特征在于,根据以下公式计算所述误差姿态角:
θe=θr-θ
φe=φr-φ
γe=γr-γ
其中,
θe为误差俯仰角,φe为误差偏航角,γe为误差滚转角;
θr为指令俯仰角,φr为指令偏航角,γr为指令滚转角;
θ为导航俯仰角,φ为导航偏航角,γ为导航滚转角。
6.根据权利要求5所述的指令力矩确定方法,其特征在于,根据以下公式计算所述形变再入飞行器形变前调姿指令力矩:
Mx0=KPIx·(γe+∫γedt)+kdx·(ωx-ωxr)
My0=KPIy·(φe+∫φedt)+kdy·ωy
Mz0=KPIz·(θe+∫θedt)+kdz·ωz
其中,
Mx0为X轴的形变前调姿指令力矩,My0为Y轴的形变前调姿指令力矩,Mz0为Z轴的形变前调姿指令力矩;
KPIx为X轴的比例项和积分项系数,KPIy为Y轴的比例项和积分项系数,KPIz为Z轴的比例项和积分项系数;
kdx为X轴的阻尼系数,kdy为Y轴的阻尼系数,kdz为Z轴的阻尼系数;
ωx为导航得到的角速度在飞行器的X轴的分量,ωy为导航得到的角速度在飞行器的Y轴的分量,ωz为导航得到的角速度在飞行器的Z轴的分量。
7.根据权利要求6所述的指令力矩确定方法,其特征在于,在调姿的过程中,若连续N帧满足以下公式,则判断调姿到位:
(|θe|<Lev)&&(|φe|<Lev)&&(|γe|<Lev)
其中,Lev为姿控精度门限。
8.根据权利要求7所述的指令力矩确定方法,其特征在于,所述获取所述形变再入飞行器变形前的启旋指令力矩包括:
令X轴、Y轴、Z轴的比例项和积分项系数为零;
按照以下公式计算启旋指令力矩:
Mx1=kdx·(ωx-ωxr)
My1=kdy·ωy
Mz1=kdz·ωz
其中,
Mx1为X轴的启旋指令力矩;
My1为Y轴的启旋指令力矩;
Mz1为Z轴的启旋指令力矩。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910276877.6A CN110058603B (zh) | 2019-04-08 | 2019-04-08 | 一种形变再入飞行器变形前指令力矩确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910276877.6A CN110058603B (zh) | 2019-04-08 | 2019-04-08 | 一种形变再入飞行器变形前指令力矩确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110058603A true CN110058603A (zh) | 2019-07-26 |
CN110058603B CN110058603B (zh) | 2022-03-04 |
Family
ID=67318568
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910276877.6A Active CN110058603B (zh) | 2019-04-08 | 2019-04-08 | 一种形变再入飞行器变形前指令力矩确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110058603B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110750102A (zh) * | 2019-11-25 | 2020-02-04 | 北京电子工程总体研究所 | 形变再入飞行器变形前指令力矩确定方法 |
CN111159942A (zh) * | 2019-12-26 | 2020-05-15 | 北京电子工程总体研究所 | 一种基于定常模拟的带翼飞行器滚转阻尼力矩计算方法 |
CN112596532A (zh) * | 2020-11-19 | 2021-04-02 | 北京电子工程总体研究所 | 一种h个空间飞行器姿控发动机控制指令动态分配方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003252299A (ja) * | 2002-02-26 | 2003-09-10 | Mitsubishi Electric Corp | 人工衛星の姿勢制御装置 |
US20140263823A1 (en) * | 2013-01-10 | 2014-09-18 | SZ DJI Technology Co., Ltd | Transformable aerial vehicle |
CN107908203A (zh) * | 2017-11-07 | 2018-04-13 | 哈尔滨工业大学 | 一种空间柔性电帆自旋调速控制方法 |
CN110320927A (zh) * | 2019-07-26 | 2019-10-11 | 中国人民解放军国防科技大学 | 智能变形飞行器的飞行控制方法及系统 |
CN112325710A (zh) * | 2020-09-24 | 2021-02-05 | 北京航天自动控制研究所 | 一种运载火箭大推力直接入轨高精度姿态控制方法和系统 |
-
2019
- 2019-04-08 CN CN201910276877.6A patent/CN110058603B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003252299A (ja) * | 2002-02-26 | 2003-09-10 | Mitsubishi Electric Corp | 人工衛星の姿勢制御装置 |
US20140263823A1 (en) * | 2013-01-10 | 2014-09-18 | SZ DJI Technology Co., Ltd | Transformable aerial vehicle |
CN107908203A (zh) * | 2017-11-07 | 2018-04-13 | 哈尔滨工业大学 | 一种空间柔性电帆自旋调速控制方法 |
CN110320927A (zh) * | 2019-07-26 | 2019-10-11 | 中国人民解放军国防科技大学 | 智能变形飞行器的飞行控制方法及系统 |
CN112325710A (zh) * | 2020-09-24 | 2021-02-05 | 北京航天自动控制研究所 | 一种运载火箭大推力直接入轨高精度姿态控制方法和系统 |
Non-Patent Citations (7)
Title |
---|
HAITAO SHUI, SHENGJUN PENG, XUN LI, AND HONGXU MA: "Coordinated Manipulator and Spacecraft Motion Planning for Free-Floating Space Robots", 《IEEE INTERNATIONAL CONFERENCE ON ROBOTICS AND BIOMIMETICS (ROBIO)》 * |
SANG-HOON YOON,HAESEONG CHO,JUNHEE LEE,CHONGAM KIM: "Effects of camber angle on aerodynamic performance of flapping-wing micro air vehicle", 《JOURNAL OF FLUIDS AND STRUCTURES》 * |
万航等: "基于动态逆的空天变体飞行器姿态控制", 《空天防御》 * |
岳彩红,唐胜景,郭 杰,王 肖,张浩强: "高超声速伸缩式变形飞行器再入轨迹快速优化", 《系统工程与电子技术》 * |
徐天富: "柔性自旋飞行器动力学分析与稳定性研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅱ辑》 * |
薛辰: "折叠翼飞行器变形过程飞行控制研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 * |
许来斌等: "M型机翼变体飞机动力学分析", 《北京理工大学学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110750102A (zh) * | 2019-11-25 | 2020-02-04 | 北京电子工程总体研究所 | 形变再入飞行器变形前指令力矩确定方法 |
CN111159942A (zh) * | 2019-12-26 | 2020-05-15 | 北京电子工程总体研究所 | 一种基于定常模拟的带翼飞行器滚转阻尼力矩计算方法 |
CN111159942B (zh) * | 2019-12-26 | 2023-09-15 | 北京电子工程总体研究所 | 一种基于定常模拟的带翼飞行器滚转阻尼力矩计算方法 |
CN112596532A (zh) * | 2020-11-19 | 2021-04-02 | 北京电子工程总体研究所 | 一种h个空间飞行器姿控发动机控制指令动态分配方法 |
CN112596532B (zh) * | 2020-11-19 | 2022-09-23 | 北京电子工程总体研究所 | 一种h个空间飞行器姿控发动机控制指令动态分配方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110058603B (zh) | 2022-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104035335B (zh) | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导方法 | |
CN110058603A (zh) | 一种形变再入飞行器变形前指令力矩确定方法 | |
JP6989252B2 (ja) | エネルギー効率に優れた衛星の操縦 | |
Sarigul-Klijn | Survey of planetary entry guidance algorithms | |
CN105807781B (zh) | 一种基于比例导引的空间近距飞越末制导方法 | |
CN104570742B (zh) | 基于前馈pid控制的异面交叉快变轨道快速高精度相对指向控制方法 | |
CN108803649A (zh) | 一种垂直起降重复使用运载器自抗扰滑模控制方法 | |
Wolf et al. | Systems for pinpoint landing at Mars | |
CN108663936B (zh) | 模型不确定航天器无退绕姿态跟踪有限时间控制方法 | |
CN108820264A (zh) | 一种用于清除空间碎片的绳系拖曳系统及方法 | |
CN104281150A (zh) | 一种姿态机动的轨迹规划方法 | |
CN103942085A (zh) | 一种舰载直升机飞行模拟器的着舰仿真方法 | |
CN106248082B (zh) | 一种飞行器自主导航系统及导航方法 | |
CN112541225A (zh) | 一种面向火箭回收大姿态翻转的虚实结合仿真模拟方法 | |
CN109774977A (zh) | 一种基于四元数的时间最优的卫星姿态快速机动方法 | |
CN111806734B (zh) | 重力补偿系统、航天器相对运动地面悬吊实验系统及方法 | |
CN105799949B (zh) | 一种亚轨道卫星的压心设计方法、姿态控制方法和系统 | |
CN110597274B (zh) | 一种适应姿态重定向的sgcmg动态框架角速度确定方法 | |
Shen et al. | A 6DOF mathematical model of parachute in Mars EDL | |
CN110723316B (zh) | 一种sgcmg的框架角速度确定方法 | |
CN104850128B (zh) | 一种用于具有大惯量积航天器的动量轮布局配置方法 | |
CN104932520B (zh) | 一种存在测量误差下的航天器雨滴形状悬停构型控制方法 | |
CN107870063B (zh) | 基于动量守恒的航天器转动惯量在轨测量方法 | |
Mavrakis et al. | On-orbit robotic grasping of a spent rocket stage: Grasp stability analysis and experimental results | |
CN110466803A (zh) | 基于等倾角姿态控制的自旋稳定卫星姿态预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |