CN103136444B - 一种导航卫星不同控制模式下太阳帆板的受照分析方法 - Google Patents
一种导航卫星不同控制模式下太阳帆板的受照分析方法 Download PDFInfo
- Publication number
- CN103136444B CN103136444B CN201310032558.3A CN201310032558A CN103136444B CN 103136444 B CN103136444 B CN 103136444B CN 201310032558 A CN201310032558 A CN 201310032558A CN 103136444 B CN103136444 B CN 103136444B
- Authority
- CN
- China
- Prior art keywords
- satellite
- sun
- coordinate system
- unit vector
- orbit
- 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)
- Navigation (AREA)
Abstract
一种导航卫星不同控制模式下太阳帆板的受照分析方法,根据太阳星历,计算目标时刻太阳在地心赤道惯性坐标系的单位矢量;根据卫星的初始轨道参数,计算目标时刻卫星的轨道参数;根据目标时刻卫星的轨道参数,计算目标时刻太阳在轨道坐标系的单位矢量;根据卫星的主要飞行事件及时序确定目标时刻卫星的工作模式;根据目标时刻卫星的工作模式及太阳在轨道坐标系的单位矢量,得到目标时刻帆板法向和太阳在卫星本体坐标系的单位矢量;根据目标时刻帆板法向和太阳在卫星本体系的单位矢量,得到帆板法向和太阳光线的夹角,完成太阳帆板受照分析。本发明对帆板法线与太阳的夹角进行事前预报,为飞控期间电源状态监视、功率计算以及热分析提供依据。
Description
技术领域
本发明涉及一种导航卫星不同控制模式下太阳帆板的受照分析方法,特别涉及一种转移轨道阶段和工作轨道段的太阳帆板受照分析方法。
背景技术
太阳帆板是卫星航天器在轨运行时的主要能量来源。太阳帆板的输出功率,是卫星在轨运行期间关注的一项重要指标。对于面积恒定的太阳帆板,太阳帆板法线与太阳光线的夹角决定了太阳帆板的有效受照面积,直接影响其输出功率。
对于导航系列卫星,星箭分离后通常要经过三到四次远地点变轨后才能进入准工作轨道。卫星在变轨期间,可能会出现太阳帆板法线与太阳光线的夹角较大,帆板输出功率较小,蓄电池供电的现象。然而目前已有方法主要是通过在轨遥测参数来估算帆板输出功率和太阳帆板法线与太阳光线的夹角。该方法属于事后处理,不利于卫星飞控期间的电源和热控分系统的状态监视。
从以上分析可知,有必要提供一种几何方法对卫星在轨运行期间不同控制模式下尤其是飞控期间转移轨道阶段的太阳帆板受照情况进行事前预报,可为卫星在轨运行期间制定卫星监视手册、电源和热控分系统的功率计算和热分析提供依据。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提出一种适用于导航卫星不同控制模式下太阳帆板的受照分析方法。
本发明的技术解决方案是:
一种导航卫星不同控制模式下太阳帆板的受照分析方法,步骤如下:
(1)根据太阳星历,计算目标时刻太阳在地心赤道惯性坐标系中的单位矢量;
(2)根据卫星的初始轨道参数,利用卫星预报模型计算得到目标时刻卫星的轨道参数;
(3)根据目标时刻卫星的轨道参数,计算地心赤道惯性坐标系到卫星轨道坐标系的转换矩阵,得到目标时刻太阳在卫星轨道坐标系的单位矢量;
(4)根据卫星的主要飞行事件及时序确定目标时刻卫星的工作模式;
(5)根据目标时刻卫星的工作模式及太阳在卫星轨道坐标系的单位矢量,得到目标时刻帆板法向和太阳在卫星本体坐标系的单位矢量;
(6)根据目标时刻太阳在卫星本体坐标系中的单位矢量Sunbody和帆板法向在卫星本体中的单位矢量FBbody,得到帆板法向和太阳光线的夹角θS=cos-1(Sunbody*FBbody),完成太阳帆板的受照分析。
所述步骤(3)中根据目标时刻卫星的轨道参数,计算地心赤道惯性坐标系到卫星轨道坐标系的转换矩阵,具体为:
其中,Loi为地心赤道惯性坐标系到卫星轨道坐标系的转换矩阵,i为目标时刻卫星轨道参数中的倾角;Ω为目标时刻卫星轨道参数中的升交点赤经;ω为目标时刻卫星轨道参数中的近地点幅角;u为目标时刻卫星轨道参数中的轨道幅角。
所述步骤(5)中根据卫星的工作模式及太阳在卫星轨道坐标系中的单位矢量,得到帆板法向和太阳在卫星本体坐标系的单位矢量,具体如下:
a.当卫星工作在太阳捕获模式时,帆板法向和太阳在卫星本体坐标系的单位矢量均为(0,0,-1);
b.当卫星工作在地球指向模式时,若卫星帆板法向指向星本体-Z轴,则帆板法向在卫星本体坐标系的单位矢量为(0,0,-1);若卫星帆板法向指向星本体-X轴,则帆板法向在卫星本体坐标系的单位矢量为(-1,0,0);若卫星姿态为地球指向姿态,则太阳在卫星本体坐标系的单位矢量等于太阳在卫星轨道坐标系中的单位矢量;若卫星姿态为点火姿态,则太阳在卫星本体坐标系的单位矢量SunBody可以通过公式SunBody=LboSunOrbit得到,其中
SunOrbit为太阳在卫星轨道坐标系的单位矢量,为点火姿态偏航角;
c.当卫星工作在远地点点火模式时,帆板法向在卫星本体坐标系的单位矢量为(-1,0,0),太阳在卫星本体坐标系的单位矢量SunBody可以通过公式SunBody=LboSunOrbit得到,其中SunOrbit为太阳在卫星轨道坐标系的单位矢量,为点火姿态偏航角;
d.当卫星工作在正常模式时,若卫星姿态处于动态偏置,则帆板法向和
太阳在卫星本体坐标系的单位矢量均为若卫星姿态处于零偏置,则帆板法向在卫星本体坐标系的单位矢量为太阳在卫星本体坐标系的单位矢量等于太阳在卫星轨道坐标系中的单位矢量(SOx,SOy,SOz)。
本发明与现有技术相比的优点在于:
(1)本发明方法属于事前预报,相比现有技术的事后处理,可为卫星在轨运行期间尤其是飞控期间电源和热控分系统的状态监视、功率计算以及热分析提供依据。
(2)本发明方法结合卫星的姿态要求,通过分析太阳、卫星以及帆板指向的关系,采用几何分析的方法进行预报。该方法可以推广应用到其他已知姿态和帆板控制信息的卫星中。
附图说明
图1为本发明的控制流程框图;
图2为本发明导航卫星的主要飞行模式流程示意图;
具体实施方式
本发明针对导航卫星在轨运行期间不同控制模式下卫星的姿态要求,通过分析太阳、卫星以及帆板指向三者之间的关系,采用几何分析的方法对太阳帆板法线与太阳光线的夹角进行的事前预报,可用于分析太阳帆板的输出功率。
如图1所示,为本发明控制方法的流程框图,具体步骤如下:
导航卫星在轨运行期间设置不同的控制模式,通过以下步骤分析不同模式下太阳帆板法线与太阳光线的夹角,为太阳帆板的功率分析提供参考。
(1)根据太阳星历,计算目标时刻太阳在地心赤道惯性坐标系中的单位矢量;
(2)根据卫星的初始轨道参数,利用卫星预报模型计算得到目标时刻卫星的轨道参数;
(3)根据目标时刻卫星的轨道参数,计算地心赤道惯性坐标系到卫星轨道坐标系的转换矩阵,得到目标时刻太阳在卫星轨道坐标系的单位矢量;
计算方法如下:
SunOrbit=LoiSunECI
其中,SunOrbit为太阳在卫星轨道坐标系中的单位矢量;SunECI为太阳在卫星轨道坐标系中的单位矢量;
转换矩阵为:
其中,Loi为地心赤道惯性坐标系到卫星轨道坐标系的转换矩阵,i为目标时刻卫星轨道参数中的倾角;Ω为目标时刻卫星轨道参数中的升交点赤经;ω为目标时刻卫星轨道参数中的近地点幅角;u为目标时刻卫星轨道参数中的轨道幅角.
(4)根据卫星的主要飞行事件及时序确定目标时刻卫星的工作模式,所述目标时刻卫星的工作模式包括卫星的姿态和帆板法向指向;
如图2所示,为导航卫星的主要飞行模式流程示意图。导航卫星星箭分离后首先进入对日定向的太阳捕获模式,之后要经过三到四次远地点变轨进入准工作轨道,最后在位置保持模式下通过相位捕获进入工作轨道,运行在正常工作模式下。对于前两次变轨,卫星首先从太阳捕获模式转入地球指向模式进行变轨前的状态设置,接着进入远地点点火模式进行变轨,变轨结束后转入地球指向模式,待状态恢复完毕后转入太阳捕获模式。对于最后几次变轨,卫星变轨结束入地球指向模式,不再转入太阳捕获模式。
导航卫星共有的主要工作模式包括:太阳捕获模式、地球指向模式、远地点点火模式、正常模式。
(5)根据目标时刻卫星的工作模式及太阳在卫星轨道坐标系的单位矢量,得到目标时刻帆板法向和太阳在卫星本体坐标系的单位矢量;具体如下:
a、当卫星工作在太阳捕获模式时,帆板法向和太阳在卫星本体坐标系的单位矢量均为(0,0,-1);
b、当卫星工作在地球指向模式时,若卫星帆板法向指向星本体-Z轴,则帆板法向在卫星本体坐标系的单位矢量为(0,0,-1);若卫星帆板法向指向星本体-X轴,则帆板法向在卫星本体坐标系的单位矢量为(-1,0,0);若卫星姿态为地球指向姿态,则太阳在卫星本体坐标系的单位矢量等于太阳在卫星轨道坐标系中的单位矢量;若卫星姿态为点火姿态,则太阳在卫星本体坐标系的单位矢量SunBody可以通过公式SunBody=LboSunOrbit得到,其中
SunOrbit为太阳在卫星轨道坐标系的单位矢量,为点火姿态偏航角;
c、当卫星工作在远地点点火模式时,帆板法向在卫星本体坐标系的单位矢量为(-1,0,0),太阳在卫星本体坐标系的单位矢量SunBody可以通过公式SunBody=LboSunOrbit得到,其中
SunOrbit为太阳在卫星轨道坐标系的单位矢量,为点火姿态偏航角;
d、当卫星工作在正常模式时,若卫星姿态处于动态偏置,则帆板法向和太阳在卫星本体坐标系的单位矢量均为若卫星姿态处于零偏置,则帆板法向在卫星本体坐标系的单位矢量为太阳在卫星本体坐标系的单位矢量等于太阳在卫星轨道坐标系中的单位矢量(SOx,SOy,SOz)。
(6)根据目标时刻太阳在卫星本体坐标系中的单位矢量Sunbody和帆板法向在卫星本体中的单位矢量FBbody,得到帆板法向和太阳光线的夹角θS=cos-1(Sunbody*FBbody),完成太阳帆板的受照分析。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (2)
1.一种导航卫星不同控制模式下太阳帆板的受照分析方法,其特征在于步骤如下:
(1)根据太阳星历,计算目标时刻太阳在地心赤道惯性坐标系中的单位矢量;
(2)根据卫星的初始轨道参数,利用卫星预报模型计算得到目标时刻卫星的轨道参数;
(3)根据目标时刻卫星的轨道参数,计算地心赤道惯性坐标系到卫星轨道坐标系的转换矩阵,得到目标时刻太阳在卫星轨道坐标系的单位矢量;
(4)根据卫星的主要飞行事件及时序确定目标时刻卫星的工作模式;
(5)根据目标时刻卫星的工作模式及太阳在卫星轨道坐标系的单位矢量,得到目标时刻帆板法向和太阳在卫星本体坐标系的单位矢量;
(6)根据目标时刻太阳在卫星本体坐标系中的单位矢量Sunbody和帆板法向在卫星本体坐标系中的单位矢量FBbody,得到帆板法向和太阳光线的夹角θS=cos-1(Sunbody*FBbody),完成太阳帆板的受照分析;
所述步骤(3)中根据目标时刻卫星的轨道参数,计算地心赤道惯性坐标系到卫星轨道坐标系的转换矩阵,具体为:
其中,Loi为地心赤道惯性坐标系到卫星轨道坐标系的转换矩阵,i为目标时刻卫星轨道参数中的倾角;Ω为目标时刻卫星轨道参数中的升交点赤经;ω为目标时刻卫星轨道参数中的近地点幅角;u为目标时刻卫星轨道参数中的轨道幅角。
2.根据权利要求1所述的一种导航卫星不同控制模式下太阳帆板的受照分析方法,其特征在于:所述步骤(5)中根据卫星的工作模式及太阳在卫星轨道坐标系中的单位矢量,得到帆板法向和太阳在卫星本体坐标系的单位矢量,具体如下:
a.当卫星工作在太阳捕获模式时,帆板法向和太阳在卫星本体坐标系的单位矢量均为(0,0,-1);
b.当卫星工作在地球指向模式时,若卫星帆板法向指向卫星本体坐标系的Z轴,则帆板法向在卫星本体坐标系的单位矢量为(0,0,-1);若卫星帆板法向指向卫星本体坐标系的X轴,则帆板法向在卫星本体坐标系的单位矢量为(-1,0,0);若卫星姿态为地球指向姿态,则太阳在卫星本体坐标系的单位矢量等于太阳在卫星轨道坐标系中的单位矢量;若卫星姿态为点火姿态,则太阳在卫星本体坐标系的单位矢量SunBody可以通过公式SunBody=LboSunOrbit得到,其中
SunOrbit为太阳在卫星轨道坐标系的单位矢量,为点火姿态偏航角;
c.当卫星工作在远地点点火模式时,帆板法向在卫星本体坐标系的单位矢量为(-1,0,0),太阳在卫星本体坐标系的单位矢量SunBody可以通过公式SunBody=LboSunOrbit得到,其中SunOrbit为太阳在卫星轨道坐标系的单位矢量,为点火姿态偏航角;
d.当卫星工作在正常模式时,若卫星姿态处于动态偏置,则帆板法向和太阳在卫星本体坐标系的单位矢量均为若卫星姿态处于零偏置,则帆板法向在卫星本体坐标系的单位矢量为太阳在卫星本体坐标系的单位矢量等于太阳在卫星轨道坐标系中的单位矢量(SOx,SOy,SOz)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310032558.3A CN103136444B (zh) | 2013-01-29 | 2013-01-29 | 一种导航卫星不同控制模式下太阳帆板的受照分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310032558.3A CN103136444B (zh) | 2013-01-29 | 2013-01-29 | 一种导航卫星不同控制模式下太阳帆板的受照分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103136444A CN103136444A (zh) | 2013-06-05 |
CN103136444B true CN103136444B (zh) | 2015-10-21 |
Family
ID=48496263
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310032558.3A Active CN103136444B (zh) | 2013-01-29 | 2013-01-29 | 一种导航卫星不同控制模式下太阳帆板的受照分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103136444B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103407581A (zh) * | 2013-07-12 | 2013-11-27 | 西北工业大学 | 一体化薄膜太阳能电池太阳帆 |
CN103488833B (zh) * | 2013-09-24 | 2017-01-04 | 北京空间飞行器总体设计部 | 一种导航卫星复杂模型的热辐射力建模方法 |
CN106096148B (zh) * | 2016-06-14 | 2019-05-24 | 中国空间技术研究院 | 一种简单姿态控制下的大倾角轨道卫星太阳帆板指向方法 |
CN109063259B (zh) * | 2018-07-02 | 2023-03-31 | 上海卫星工程研究所 | 地球静止轨道大口径光学成像卫星载荷遮阳板设计方法 |
CN110466807B (zh) * | 2019-08-12 | 2020-11-20 | 北京控制工程研究所 | 一种惯性定向小角度旋转卫星的帆板控制方法 |
CN111846289B (zh) * | 2020-08-06 | 2021-11-16 | 上海航天控制技术研究所 | 一种太阳帆板偏置安装时的卫星对日定向控制方法及卫星 |
CN112613197B (zh) * | 2021-01-08 | 2021-08-24 | 中国科学院微小卫星创新研究院 | 一种星体对转动太阳帆板遮挡面积的分析方法 |
CN113624238B (zh) * | 2021-08-11 | 2024-05-14 | 中国科学院微小卫星创新研究院 | 星间链路指向恢复方法 |
CN114252151B (zh) * | 2021-11-29 | 2024-08-23 | 中国科学院国家空间科学中心 | 一种星载光学仪器在轨工作模式自主调整方法及其装置 |
CN114313318B (zh) * | 2021-12-31 | 2024-02-02 | 浙江时空道宇科技有限公司 | 卫星太阳能帆板转角处理方法、装置、卫星及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102105818A (zh) * | 2009-03-18 | 2011-06-22 | 株式会社博思科 | 日照量的评价方法和评价装置 |
CN102567586A (zh) * | 2012-01-04 | 2012-07-11 | 中国电子科技集团公司第二十八研究所 | 一种航天任务三维综合态势的仿真系统及仿真方法 |
CN102565797A (zh) * | 2011-12-21 | 2012-07-11 | 北京航空航天大学 | 一种针对聚束模式星载sar图像的几何校正方法 |
-
2013
- 2013-01-29 CN CN201310032558.3A patent/CN103136444B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102105818A (zh) * | 2009-03-18 | 2011-06-22 | 株式会社博思科 | 日照量的评价方法和评价装置 |
CN102565797A (zh) * | 2011-12-21 | 2012-07-11 | 北京航空航天大学 | 一种针对聚束模式星载sar图像的几何校正方法 |
CN102567586A (zh) * | 2012-01-04 | 2012-07-11 | 中国电子科技集团公司第二十八研究所 | 一种航天任务三维综合态势的仿真系统及仿真方法 |
Non-Patent Citations (4)
Title |
---|
Variable structure control for satellite attitude stabilization in elliptic orbits using solar radiation pressure;T.R.Patel et al.;《ACTA ASTRONAUTICA》;20090228;第64卷(第2-3期);第359~373页 * |
三维稳定卫星姿态确定和姿态控制系统研究;柏林;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20020615(第01期);第11~23页 * |
双向反射分布函数在空间目标可见光反射特性建模中的应用;孙成明 等;《光学技术》;20080915;第34卷(第5期);第750~753,757页 * |
组合体太阳帆板遮挡分析_;黄海兵 等;《系统仿真学报》;20090605;第21卷(第11期);第3215~3218页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103136444A (zh) | 2013-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103136444B (zh) | 一种导航卫星不同控制模式下太阳帆板的受照分析方法 | |
CN104142686B (zh) | 一种卫星自主编队飞行控制方法 | |
CN103112603B (zh) | 欠驱动高速自旋卫星建立正常姿态的方法 | |
CN103488166B (zh) | 一种整星全动量管理功能测试系统及测试方法 | |
CN103676941B (zh) | 基于运动学和动力学模型的卫星控制系统故障诊断方法 | |
CN102759358B (zh) | 基于失效卫星表面参考点的相对位姿动力学建模方法 | |
CN104176275B (zh) | 一种使用动量轮与磁力矩器联合的速率阻尼方法 | |
CN105716615A (zh) | 火星探测器自主导航方法 | |
CN104181930B (zh) | 一种倾斜轨道卫星偏航机动的自主控制方法 | |
CN103955224B (zh) | 一种用于相对运动视线跟踪的姿态控制方法 | |
CN103676918B (zh) | 一种基于未知输入观测器的卫星执行机构故障诊断方法 | |
CN108897023B (zh) | 一种星上自主的非合作机动目标跟踪保持变轨方法 | |
CN103019247A (zh) | 一种火星探测器无陀螺自主空间姿态机动控制方法 | |
CN104298128A (zh) | 空间飞行器导航制导技术地面仿真方法 | |
CN103264776A (zh) | 一种基于信息融合的控制系统工作模式设置及切换方法 | |
Wolf et al. | Toward improved landing precision on Mars | |
CN102004491A (zh) | 一种卫星初入轨段的初始太阳捕获方法 | |
CN105947239B (zh) | 一种基于矢量投影的航天器光学敏感器视场分析方法 | |
CN104021241A (zh) | 一种导航卫星复杂模型的地球反照光压摄动建模方法 | |
Qiao et al. | Asteroid approaching orbit optimization considering optical navigation observability | |
CN102514734A (zh) | 基于日地系统Halo轨道探测器构型与姿态指向的姿态递推方法 | |
CN104765373A (zh) | 一种星上相对运动状态获取方法 | |
CN104950668A (zh) | 卫星编队解析式燃料优化控制方法及系统 | |
CN103274066B (zh) | 从Halo轨道出发探测深空目标的逃逸轨道设计方法 | |
CN105573332A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |