CN104375512A - 一种基于频谱分析的航天器机动路径优化方法 - Google Patents
一种基于频谱分析的航天器机动路径优化方法 Download PDFInfo
- Publication number
- CN104375512A CN104375512A CN201410602042.2A CN201410602042A CN104375512A CN 104375512 A CN104375512 A CN 104375512A CN 201410602042 A CN201410602042 A CN 201410602042A CN 104375512 A CN104375512 A CN 104375512A
- Authority
- CN
- China
- Prior art keywords
- tau
- max
- omega
- motor
- driven
- 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
Landscapes
- Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
一种基于频谱分析的航天器机动路径优化方法,步骤为:(1)获取航天器某单个姿态轴的机动角度和执行机构能力,基于标准正弦机动路径规划方法,获取初始正弦机动路径;(2)计算初始正弦机动路径的频谱函数及航天器帆板基频处的振幅;(3)利用振幅频谱上界约束条件,解算得到航天器帆板基频与机动加速时长的关系,得到优化后的加速时间和优化后的最大角加速度;(4)基于优化后的最大角加速度及加速时间,在机动角度约束下,根据正弦机动路径设计方法,重新获取机动滑行时间由此得到最终的优化后的正弦机动路径。本发明方法适用于带有大型挠性帆板航天器姿态机动过程的姿态角加速度设计。
Description
技术领域
本发明属于控制领域,涉及一种航天器轨道的优化方法。
背景技术
根据挠性帆板的振动方程以及帆板的模态频率特性可知,帆板振动模态容易在其固有频率点处被外激励信号所激发。星体的角加速度作为挠性帆板的外激励,其在模态频率处的幅值大小直接决定了帆板模态被激发的程度。从设计角度考虑,为了缩短航天器机动完成后的稳定时间,应尽量减小机动过程对帆板产生的激励。这就需要减小星体角加速度在帆板模态频率处的幅值。
机动路径优化的主要目的是从避免激起挠性帆板振动的角度,对已设计好的机动路径进一步优化。以往的机动路径设计及优化主要考虑执行机构能力及机动时间的约束,而执行机构能力与机动时间不能反映航天器挠性帆板的物理特性,从而在设计机动路径时不能主动规避帆板基频,存在机动过程激励帆板大幅振动的可能。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种基于频谱分析的航天器机动路径优化方法,考虑了帆板基频的影响,从帆板挠性频率分析的角度,显式的给出了航天器机动角加速度需满足的与帆板基频相关的约束条件,并在此基础上实现了对航天器机动路径的优化,适用于带有大型挠性帆板航天器姿态机动过程的姿态角加速度设计。
本发明的技术解决方案是:一种基于频谱分析的航天器机动路径优化方法,包括如下步骤:
(1)获取航天器某单个姿态轴的机动角度θd、执行机构的能力τmax,基于标准正弦机动路径规划方法,获取初始正弦机动最大角加速度amax、最大角速度ωmax、机动加速时间τ、滑行时间τ1,由此得到初始正弦机动路径为:
其中,机动加速时间 滑行时间 如果计算得到的滑行时间τ1≤0,则确定滑行时间τ1=0,并根据角度约束条件θd调整机动加速时间为
(2)计算所获得的初始正弦机动路径的频谱函数F(ω),
由此得到振幅频谱的上界约束条件
(3)计算航天器帆板基频ω1处的振幅|F(ω1)|,
(4)利用步骤(2)得到的振幅频谱上界约束条件,对步骤(3)中获取的|F(ω1)|提出约束,要求由此解算得到航天器帆板基频ω1与机动加速时长τ的关系,作为对加速时间的约束,即在此基础上结合步骤(1)确定的加速时间,得到优化后的加速时间
(5)根据步骤(4)中获取的对加速时间的约束计算得到对最大角加速度的约束,综合初始的最大角加速度amax,得到优化后的最大角加速度
(6)基于优化后的最大角加速度及加速时间在机动角度约束下,根据正弦机动路径设计方法,重新获取机动滑行时间当τ1>0时,优化 当τ1=0,优化后的
(7)利用步骤(6)的计算结果,得到最终的优化后的正弦机动路径为:
本发明与现有技术相比的优点在于:
(1)本发明方法通过研究机动过程姿态角加速度的频谱特性,给出了机动角加速度需满足的与帆板基频相关的约束条件,并在此基础上对设计的机动角加速度进行优化,主动规避了在帆板固有频率点处激发帆板振动的可能性,尤其适用于带有大型挠性帆板航天器姿态机动过程的姿态角加速度设计;
(2)本发明方法显式的给出了航天器机动角加速度需满足的与帆板基频相关的约束条件,设计人员可以根据帆板特性选择机动过程需避开的频域范围,设计思想简单,易于工程实现。
附图说明
图1为本发明方法的流程框图;
图2为标准正弦机动路径。
具体实施方式
对于带有大型挠性附件的复杂航天器,对机动路径进行优化以避免机动过程对挠性帆板的激励,是提高航天器机动到位后姿态控制精度的重要手段。
由分析可知,正弦型路径机动速度的幅值谱比较复杂,但主要的频谱成分均分布在低频区域。本发明根据工程需求,设计星体角加速度的主要频谱成分分布在小于帆板的第一阶模态频率区间内;并根据设计约束,明确给出最大角加速度应满足的与帆板模态频率相关的条件,进而对机动角加速度进行优化。本发明不同于以往的机动路径优化方法,结合工程需求和设计人员的实际经验,从物理本质优化机动路径,实现了机动过程尽量避免帆板挠性振动的目的。
如图1所示,为本发明的路径优化设计流程图,具体实施步骤如下:
(1)根据执行机构能力获取航天器的最大机动角加速度amax,根据执行机构角动量容量和敏感器量程获取航天器的最大受限角速度ωmax;根据航天器某姿态轴机动角度θd(单个姿态轴)及航天器转动惯量J的约束,基于标准正弦机动路径规划方法,获取机动加速时间滑行时间如果计算得到的滑行时间τ1≤0,则确定卫星的滑行时间τ1=0,并根据角度约束条件θd调整机动加速时间即正弦函数的周期
根据标准正弦机动路径设计方法,获取最大角速度ωmax、机动加速时间τ以及滑行时间τ1之后,可以构造出标准正弦机动路径,如图2所示,其角加速度曲线f(t)的表达式可以写为:
(2)对所获得的正弦机动路径进行分析,计算其频谱函数,为 其振幅频谱满足经放缩,可以获得振幅频谱的上界,即
(3)假设航天器帆板基频为ω1,则该频率处的幅值为
(4)根据工程需求,对步骤(3)中获取的|F(ω1)|与振幅频谱上界的比值提出约束,要求为了避免设计的机动路径对柔性帆板的激励,通常选取0<α<1。
(5)根据该约束,解算帆板基频ω1与机动加速时长τ的关系,得到对加速段时长的约束,即
(6)将步骤(5)中获取的对加速时间的约束代入步骤(1),得到对最大角加速度的约束,即(当τ1>0时),(当τ1=0时);优化最大角加速度为(当τ1>0时),或(当τ1=0时);结合步骤(1)设计的加速时间,选取优化后的加速时间
(7)基于优化的最大角加速度及加速时间对机动滑行时间τ1进行优化,即当τ1>0时,优化当τ1=0,优化后的从而最终实现航天器基于频谱分析的机动路径优化,优化机动路径为:
实施例
(1)以某典型复杂航天器滚动轴机动70°过程的路径设计为例,假设航天器转动惯量2962kg.m2,执行机构为控制力矩陀螺,能提供的最大力矩为20Nm,角动量容量为75Nms,敏感器为陀螺,最大量测角速率为2.5°/s。则根据正弦机动路径获取卫星的最大机动角加速度为0.3869°/s2,最大角速度为2.5°/s。考虑到工程实现中需要为陀螺测量和执行机构留有一定的余量,因此设计路径时选取最大机动角加速度为amax=0.3831°/s2,最大角速度ωmax=2.4389°/s,正弦三角函数周期加速时间滑行时间
(2)计算获得当前设计的正弦机动路径的振幅频谱
(3)获取航天器帆板基频ω1=0.2Hz,根据步骤(2)中获取的振幅频谱,计算该频率处的幅值为
(4)根据工程需求,对步骤(3)中获取的F(ω1)与振幅频谱上界的比值提出约束,要求本实例选取α=0.1;
(5)根据该约束,解算帆板基频ω1与机动加速时长τ的关系,得 从而得到对加速段时长的约束,即
(6)将步骤(5)中获取的对加速时间的约束显然根据标准正弦机动路径得到的加速时间10s不满足此约束,选取新的加速时间并对最大角加速度进行优化。将加速时间约束代入步骤(1),得到对最大角加速度的约束,即优化最大角加速度为
(7)根据优化后的最大角加速度及加速时间重新设计正弦机动路径,匀速时间 完成航天器基于频谱分析的机动路径优化。
最终得到的机动路径角加速度函数表达形式为:
其中
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (1)
1.一种基于频谱分析的航天器机动路径优化方法,其特征在于包括如下步骤:
(1)获取航天器某单个姿态轴的机动角度θd、执行机构的能力τmax,基于标准正弦机动路径规划方法,获取初始正弦机动最大角加速度amax、最大角速度ωmax、机动加速时间τ、滑行时间τ1,由此得到初始正弦机动路径为:
其中,机动加速时间 滑行时间 如果计算得到的滑行时间τ1≤0,则确定滑行时间τ1=0,并根据角度约束条件θd调整机动加速时间为
(2)计算所获得的初始正弦机动路径的频谱函数F(ω),
由此得到振幅频谱的上界约束条件
(3)计算航天器帆板基频ω1处的振幅|F(ω1)|,
(4)利用步骤(2)得到的振幅频谱上界约束条件,对步骤(3)中获取的|F(ω1)|提出约束,要求0<α<1,由此解算得到航天器帆板基频ω1与机动加速时长τ的关系,作为对加速时间的约束,即在此基础上结合步骤(1)确定的加速时间,得到优化后的加速时间
(5)根据步骤(4)中获取的对加速时间的约束计算得到对最大角加速度的约束,综合初始的最大角加速度amax,得到优化后的最大角加速度
(6)基于优化后的最大角加速度及加速时间在机动角度约束下,根据正弦机动路径设计方法,重新获取机动滑行时间当τ1>0时,优化 当τ1=0,优化后的
(7)利用步骤(6)的计算结果,得到最终的优化后的正弦机动路径为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410602042.2A CN104375512B (zh) | 2014-10-30 | 2014-10-30 | 一种基于频谱分析的航天器机动路径优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410602042.2A CN104375512B (zh) | 2014-10-30 | 2014-10-30 | 一种基于频谱分析的航天器机动路径优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104375512A true CN104375512A (zh) | 2015-02-25 |
CN104375512B CN104375512B (zh) | 2017-01-11 |
Family
ID=52554507
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410602042.2A Active CN104375512B (zh) | 2014-10-30 | 2014-10-30 | 一种基于频谱分析的航天器机动路径优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104375512B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109269504A (zh) * | 2018-10-24 | 2019-01-25 | 北京控制工程研究所 | 一种具有末端约束的姿态机动路径规划方法 |
CN110647151A (zh) * | 2019-10-16 | 2020-01-03 | 北京京东乾石科技有限公司 | 坐标转换方法及装置、计算机可读存储介质、电子设备 |
CN110673617A (zh) * | 2019-09-23 | 2020-01-10 | 北京控制工程研究所 | 一种点对点姿态机动智能轨迹规划与调整方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102023051B (zh) * | 2010-09-29 | 2012-06-13 | 北京控制工程研究所 | 一种测量星上有效载荷三轴角位移高频微振动的方法 |
CN102519473B (zh) * | 2011-12-08 | 2014-05-28 | 北京控制工程研究所 | 一种适于帆板基频较高卫星的混合正弦机动路径引导方法 |
CN103235597B (zh) * | 2013-04-09 | 2015-05-20 | 北京理工大学 | 一种航天器的快速姿态机动快速稳定联合控制方法 |
CN104062976B (zh) * | 2014-06-10 | 2016-08-24 | 北京控制工程研究所 | 一种基于角加速度导数为正弦曲线的飞行器姿态快速机动方法 |
CN104020778B (zh) * | 2014-06-18 | 2017-07-28 | 哈尔滨工业大学 | 基于跟踪时间‑能耗最优轨线的挠性卫星姿态机动控制方法 |
-
2014
- 2014-10-30 CN CN201410602042.2A patent/CN104375512B/zh active Active
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109269504A (zh) * | 2018-10-24 | 2019-01-25 | 北京控制工程研究所 | 一种具有末端约束的姿态机动路径规划方法 |
CN109269504B (zh) * | 2018-10-24 | 2021-04-13 | 北京控制工程研究所 | 一种具有末端约束的姿态机动路径规划方法 |
CN110673617A (zh) * | 2019-09-23 | 2020-01-10 | 北京控制工程研究所 | 一种点对点姿态机动智能轨迹规划与调整方法 |
CN110673617B (zh) * | 2019-09-23 | 2020-10-23 | 北京控制工程研究所 | 一种点对点姿态机动智能轨迹规划与调整方法 |
CN110647151A (zh) * | 2019-10-16 | 2020-01-03 | 北京京东乾石科技有限公司 | 坐标转换方法及装置、计算机可读存储介质、电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN104375512B (zh) | 2017-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102538821B (zh) | 一种快速、参数分段式捷联惯性导航***自对准方法 | |
CN104062976B (zh) | 一种基于角加速度导数为正弦曲线的飞行器姿态快速机动方法 | |
CN103033186A (zh) | 一种用于水下滑翔器的高精度组合导航定位方法 | |
CN103777635A (zh) | 动力定位船舶鲁棒自适应航迹跟踪控制*** | |
CN102915036B (zh) | 一种参数不确定性飞行器倾斜角控制***极限环抑制方法 | |
JP2014527176A5 (zh) | ||
CN104375512A (zh) | 一种基于频谱分析的航天器机动路径优化方法 | |
US9045240B2 (en) | Flight control device, spacecraft, and reference trajectory correcting method | |
CN104536442A (zh) | 一种基于动态规划的水下航行器路径规划方法 | |
CN105547274A (zh) | 一种用于机抖激光陀螺的主动式抖动消振控制方法 | |
CN103454662B (zh) | 一种基于ckf的sins/北斗/dvl组合对准方法 | |
CN102162733A (zh) | 一种基于svm的auv舰位推算导航误差实时修正方法 | |
CN103217699A (zh) | 一种基于偏振信息的组合导航***递推优化初始对准方法 | |
CN103389113A (zh) | 基于惯性传感器输出误差补偿的动基座惯导初始对准方法 | |
CN103954282B (zh) | 基于加速度计输出增量的捷联惯性导航方法 | |
Park et al. | Vehicle positioning based on velocity and heading angle observer using low-cost sensor fusion | |
CN104819710A (zh) | 一种具有温度补偿结构的谐振式硅微机械陀螺 | |
CN105759827A (zh) | 一种抑制不期望柔性振动的航天器姿态控制*** | |
CN104467612A (zh) | 一种音圈电机控制方法及镜头对焦*** | |
CN103940429A (zh) | 一种惯性导航***横坐标系下载体姿态的实时测量方法 | |
CN101124456B (zh) | 产生导航参数和竖直位置的方法 | |
CN107270899A (zh) | 基于切换控制的长航时惯性导航***阻尼切换方法 | |
CN107664511A (zh) | 基于速度信息的摇摆基座粗对准方法 | |
CN102880049A (zh) | 一种基于帆板挠性形变测量的自适应振动控制方法 | |
CN103344251B (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 |