CN111776251A - 卫星的长期对日跟踪控制方法 - Google Patents
卫星的长期对日跟踪控制方法 Download PDFInfo
- Publication number
- CN111776251A CN111776251A CN202010705648.4A CN202010705648A CN111776251A CN 111776251 A CN111776251 A CN 111776251A CN 202010705648 A CN202010705648 A CN 202010705648A CN 111776251 A CN111776251 A CN 111776251A
- Authority
- CN
- China
- Prior art keywords
- satellite
- axis
- control
- representing
- coordinate system
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000007774 longterm Effects 0.000 title claims abstract description 16
- 238000005096 rolling process Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 abstract description 5
- 238000009825 accumulation Methods 0.000 abstract description 3
- 230000009471 action Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 241000037831 Polygonatum sibiricum Species 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了的一种卫星的长期对日跟踪控制方法,所述方法包括:控制卫星的X、Y轴实现卫星的‑Z轴对日跟踪;卫星的Z轴机动将卫星的X轴保持在轨道平面内。本发明通过将卫星的X轴控制到轨道平面,建立了卫星一个稳态区域,减小了星体角动量的累积,解决了安装有大型天线、太阳阵等大型展开部件的卫星重力梯度力矩大,无法长期保持对日跟踪的问题。
Description
技术领域
本发明涉及卫星姿态控制领域,特别涉及一种卫星的长期对日跟踪控制方法。
背景技术
卫星一般会设计安全模式,确保在轨运行发生严重故障时可以确保卫星安全。其中一类安全模式通过姿态控制卫星的太阳阵对日,保证卫星能源,等待地面进行故障清除后重新转入正常工作模式。
对于安装有大型天线、太阳阵等大型展开部件的卫星,其惯量积较大,导致卫星受到很大的重力梯度力矩。当卫星进入安全模式,转入对日定向和跟踪时由于重力梯度力矩的作用,导致卫星动量执行机构饱和,失去对日跟踪甚至失控。
一种解决办法是选用更大角动量的执行机构,带来的问题是卫星重量增加,付出的代价较大,甚至是不能满足卫星总体设计、运载能力的要求。
发明内容
本发明的目的是提供一种卫星的长期对日跟踪控制方法,以解决针对安装有大型天线、太阳阵等大型展开部件的卫星安全模式姿态捕获容易导致执行机构饱和的问题。
为了解决上述问题,本发明通过以下技术方案实现:
一种卫星的长期对日跟踪控制方法,包括:控制卫星的X、Y轴实现卫星的-Z轴对日跟踪;卫星的Z轴机动将卫星的X轴保持在轨道平面内。
可选地,所述控制卫星的X、Y轴实现卫星的-Z轴对日跟踪的步骤包括:
可选地,对日定向的所述第一控制力矩采用如下公式进行计算
Tcy_b=qPKPyθs+qDKDyωbiy
Tcz_b=0
式中,Tcx_b表示卫星b坐标系x方向的控制力矩;Tcy_b表示卫星b坐标系y方向的控制力矩;Tcz_b表示卫星b坐标系z方向的控制力矩;和θs表示太阳角,对和θs进行限幅,即|θs|≤π/18rad;ωbix、ωbiy和ωbiz分别为卫星b坐标系的三轴角速度;Kpi和KDi分别表示根据卫星质量特性设计的控制参数,其中i=x,y,z;qP,qD的缺省值均为1。
可选地,所述卫星的Z轴机动将卫星的X轴保持在轨道平面内的步骤包括:获取偏航方位角和指令力矩,根据所述指令力矩,卫星的Z轴机动将卫星的X轴保持在轨道平面内;
所述偏航方位角采用如下公式进行计算:
式中,ψm表示偏航方位角;
所述指令力矩采用如下公式进行计算:
Tcy_b=qPKPyθs+qDKDyωbiy
Tcz_b=qPKPzψm+qDKDzωbiy
式中,Tcx_b表示卫星b坐标系x方向的控制力矩;Tcy_b表示卫星b坐标系y方向的控制力矩;Tcz_b表示卫星b坐标系z方向的控制力矩;ψm表示偏航方位角;和θs表示太阳角,对和θs进行限幅,即|θs|≤π/18rad;ωbix、ωbiy和ωbiz分别为卫星b坐标系的三轴角速度;Kpi和KDi分别表示根据卫星质量特性设计的控制参数,其中i=x,y,z;qP,qD的缺省值均为1。
可选地,还包括:对卫星在阴影区进行姿态控制:当卫星进入阴影区时,则卫星b坐标系中的X轴、Y轴进行惯性保持控制,Z轴利用陀螺积分计算偏航方位角,
卫星b坐标系相对于惯性系的姿态角估值:
第二控制力矩采用如下公式进行计算:
式中,Tcx_b表示卫星b坐标系x方向的控制力矩;Tcy_b表示卫星b坐标系y方向的控制力矩;Tcz_b表示卫星b坐标系z方向的控制力矩;
本发明至少具有以下优点之一:
本发明提供的一种卫星的长期对日跟踪控制方法,通过将卫星的X轴控制到轨道平面,建立了卫星一个稳态区域,减小了星体角动量的累积,解决了安装有大型天线、太阳阵等大型展开部件的卫星重力梯度力矩大,无法长期保持对日跟踪的问题,本发明为一种姿态快速捕获控制技术,实现了使用小角动量执行机构可以实现此类卫星快速太阳捕获,并保持长期对日跟踪的目的。
附图说明
图1为本发明一实施例提供的一种卫星的长期对日跟踪控制方法的流程示意图。
具体实施方式
以下结合附图1和具体实施方式对本发明一种卫星的长期对日跟踪控制方法作进一步详细说明。根据下面说明,本发明的优点和特征将更清楚。需要说明的是,附图采用非常简化的形式且均使用非精准的比例,仅用以方便、明晰地辅助说明本发明实施方式的目的。为了使本发明的目的、特征和优点能够更加明显易懂,请参阅附图。须知,本说明书所附图式所绘示的结构、比例、大小等,均仅用以配合说明书所揭示的内容,以供熟悉此技术的人士了解与阅读,并非用以限定本发明实施的限定条件,故不具技术上的实质意义,任何结构的修饰、比例关系的改变或大小的调整,在不影响本发明所能产生的功效及所能达成的目的下,均应仍落在本发明所揭示的技术内容能涵盖的范围内。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
本实施例提供的一种卫星的长期对日跟踪控制方法,包括:控制卫星的X周、Y轴实现卫星的-Z轴对日跟踪;卫星的Z轴机动将卫星的X轴保持在轨道平面内。
所述控制卫星的X、Y轴实现卫星的-Z轴对日跟踪的步骤包括:根据太阳敏感器得到的太阳角和θs,陀螺测量到的卫星b坐标系的三轴角速度ωbix、ωbiy和ωbiz,计算第一控制力矩,以控制卫星的X、Y轴实现星体的-Z轴对日跟踪。
对日定向的所述第一控制力矩采用如下公式进行计算
式中,Tcx_b表示卫星b坐标系x方向的控制力矩;Tcy_b表示卫星b坐标系y方向的控制力矩;Tcz_b表示卫星b坐标系z方向的控制力矩;和θs表示太阳角,对和θs进行限幅,即|θs|≤π/18rad;ωbix、ωbiy和ωbiz分别为卫星b坐标系的三轴角速度;Kpi和KDi分别表示根据卫星质量特性设计的控制参数,其中i=x,y,z;qP,qD的缺省值均为1。
可选地,所述卫星的Z轴机动将卫星的X轴保持在轨道平面内的步骤包括:获取偏航方位角和指令力矩,根据所述指令力矩,卫星的Z轴机动将卫星的X轴保持在轨道平面内;
所述偏航方位角采用如下公式进行计算:
经过卫星轨道计算得到平黄赤交角ε,黄经Λ,太阳赤经Ωs,太阳赤纬δs,平轨道倾角h,平纬度幅角u、平升交点赤经Ω,则太阳矢量在惯性坐标系下的矢量Si可表示为
式中,Six表示太阳矢量在惯性坐标系下x方向的投影,Siy表示太阳矢量在惯性坐标系下y方向的投影,Siz表示太阳矢量在惯性坐标系下z方向的投影,其中i表示惯性坐标系;
定义两个中间变量M、N:
M=-cosΩssinδssinΩsin(h)+sinΩssinδscosΩsin(h)+cosδscos(h) (3)
N=sinΩssinΩsin(h)+cosΩscosΩsin(h)
计算J2000惯性坐标系到卫星轨道坐标系的坐标转换矩阵Aoi:
Aoi11=-sinucosΩ-cosucosisinΩ
Aoi12=-sinusinΩ+cosucosicosΩ
Aoi13=cosusini
Aoi21=-sinisinΩ
Aoi122=sinicosΩ
Aoi23=-cosi
Aoi31=-cosucosΩ+sinucosisinΩ
Aoi32=-cosusinΩ-sinucosicosΩ
Aoi33=-sinusini
计算太阳矢量在轨道面上的投影:
Sox=Aoi11cosΛ+Aoi12sinΛcosε+Aoi13sinΛsinε
Soy=Aoi21cosΛ+Aoi22sinΛcosε+Aoi23sinΛsinε
Soz=Aoi31cosΛ+Aoi32sinΛcosε+Aoi33sinΛsinε (5)
Sox表示太阳矢量在轨道坐标系下x方向的投影,Soy表示太阳矢量在轨道坐标系下y方向的投影,Soz表示太阳矢量在轨道坐标系下z方向的投影;
太阳矢量在轨道系的坐标分量v1=[Sox Soy Soz]T。
太阳矢量在卫星b坐标系的坐标分量记为r1,计算方法如下:
通过磁场公式得到磁场矢量在轨道系的坐标分量Bo,记为v2
通过磁强计,得到磁场矢量在卫星b坐标系的坐标分量Bb,记为r2。
建立如下两个坐标系
其中,r1,r2为太阳、磁场矢量在卫星b坐标系的坐标分量。
R=[R1 R2 R3]
其中,v1,v2为太阳、磁场矢量在轨道系的坐标分量。
V=[V1 V2 V3] (8)
卫星b坐标系相对轨道坐标系的姿态矩阵Ab←o为
Ab←o=RVT。
将Ab←o转换为四元数,得到qb←o,
式中,qbi表示惯性坐标系到卫星b坐标系的转换四元数,qoi表示惯性坐标系到轨道坐标系的转换四元数;
式中,qs表示卫星b坐标系转换到目标坐标系的转换四元数;qs1、qs2、qs3表示表示四元数qs的矢量部分,qs4表示表示四元数qs的标量部分;
先对qs进行归一化处理,然后进行如下计算
式中,ψm表示偏航方位角;
将ψm转换到[-π,π]的范围内,具体方法如下:
式中,ψm表示偏航方位角。
所述指令力矩采用如下公式进行计算:
式中,Tcx_b表示卫星b坐标系x方向的控制力矩;Tcy_b表示卫星b坐标系y方向的控制力矩;Tcz_b表示卫星b坐标系z方向的控制力矩;ψm表示偏航方位角;和θs表示太阳角,对和θs进行限幅,即|θs|≤π/18rad;ωbix、ωbiy和ωbiz分别为卫星b坐标系的三轴角速度;Kpi和KDi分别表示根据卫星质量特性设计的控制参数,其中i=x,y,z;qP,qD的缺省值均为1。可以理解的是,所述预设时间段可以为10s(20拍),但本发明不以此为限。
本实施例还包括:判断所述卫星是否位于阴影区,若否则重复上述过程,若是,则对卫星在阴影区进行姿态控制:当卫星进入阴影区时,则卫星b坐标系中的X轴、Y轴进行惯性保持控制,Z轴利用陀螺积分计算偏航方位角,
卫星b坐标系相对于惯性系的姿态角估值:
第二控制力矩采用如下公式进行计算:
式中,Tcx_b表示卫星b坐标系x方向的控制力矩;Tcy_b表示卫星b坐标系y方向的控制力矩;Tcz_b表示卫星b坐标系z方向的控制力矩;
由此,本实施例提供的一种卫星的长期对日跟踪控制方法,通过将卫星的X轴控制到轨道平面,建立了卫星一个稳态区域,减小了星体角动量的累积,解决了安装有大型天线、太阳阵等大型展开部件的卫星重力梯度力矩大,无法长期保持对日跟踪的问题,本发明为一种姿态快速捕获控制技术,实现了使用小角动量执行机构可以实现此类卫星快速太阳捕获,并保持长期对日跟踪的目的。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (5)
1.一种卫星的长期对日跟踪控制方法,其特征在于,包括:
控制卫星的X轴、Y轴实现卫星的-Z轴对日跟踪;
卫星的Z轴机动将卫星的X轴保持在轨道平面内。
4.如权利要求1所述的卫星的长期对日跟踪控制方法,其特征在于,
所述卫星的Z轴机动将卫星的X轴保持在轨道平面内的步骤包括:
获取偏航方位角和指令力矩,根据所述指令力矩,卫星的Z轴机动将卫星的X轴保持在轨道平面内;
所述偏航方位角采用如下公式进行计算:
式中,ψm表示偏航方位角;
所述指令力矩采用如下公式进行计算:
Tcy_b=qPKPyθs+qDKDyωbiy
Tcz_b=qPKPzψm+qDKDzωbiy
5.如权利要求1所述的卫星的长期对日跟踪控制方法,其特征在于,
还包括:对卫星在阴影区进行姿态控制:
当卫星进入阴影区时,则卫星b坐标系中的X轴、Y轴进行惯性保持控制,Z轴利用陀螺积分计算偏航方位角,
卫星b坐标系相对于惯性系的姿态角估值:
第二控制力矩采用如下公式进行计算:
式中,Tcx_b表示卫星b坐标系x方向的控制力矩;Tcy_b表示卫星b坐标系y方向的控制力矩;Tcz_b表示卫星b坐标系z方向的控制力矩;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010705648.4A CN111776251B (zh) | 2020-07-21 | 2020-07-21 | 卫星的长期对日跟踪控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010705648.4A CN111776251B (zh) | 2020-07-21 | 2020-07-21 | 卫星的长期对日跟踪控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111776251A true CN111776251A (zh) | 2020-10-16 |
CN111776251B CN111776251B (zh) | 2022-01-28 |
Family
ID=72764727
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010705648.4A Active CN111776251B (zh) | 2020-07-21 | 2020-07-21 | 卫星的长期对日跟踪控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111776251B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114485641A (zh) * | 2022-01-24 | 2022-05-13 | 武汉梦芯科技有限公司 | 一种基于惯导卫导方位融合的姿态解算方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6550721B2 (en) * | 2000-03-09 | 2003-04-22 | The Boeing Company | Safing mode for high momentum states in body stabilized spacecraft |
CN106155074A (zh) * | 2016-08-15 | 2016-11-23 | 上海航天控制技术研究所 | 一种保证星地链路的卫星三轴对日定向控制方法 |
CN106864774A (zh) * | 2017-03-13 | 2017-06-20 | 上海航天控制技术研究所 | 卫星从任意姿态到对日定向的控制方法和*** |
CN108657467A (zh) * | 2018-05-16 | 2018-10-16 | 清华大学 | 一种利用虚拟太阳矢量的航天器偏航机动控制方法及*** |
CN110228605A (zh) * | 2019-06-18 | 2019-09-13 | 北京电子工程总体研究所 | 一种基于太阳敏感器的卫星安全对日控制方法 |
-
2020
- 2020-07-21 CN CN202010705648.4A patent/CN111776251B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6550721B2 (en) * | 2000-03-09 | 2003-04-22 | The Boeing Company | Safing mode for high momentum states in body stabilized spacecraft |
CN106155074A (zh) * | 2016-08-15 | 2016-11-23 | 上海航天控制技术研究所 | 一种保证星地链路的卫星三轴对日定向控制方法 |
CN106864774A (zh) * | 2017-03-13 | 2017-06-20 | 上海航天控制技术研究所 | 卫星从任意姿态到对日定向的控制方法和*** |
CN108657467A (zh) * | 2018-05-16 | 2018-10-16 | 清华大学 | 一种利用虚拟太阳矢量的航天器偏航机动控制方法及*** |
CN110228605A (zh) * | 2019-06-18 | 2019-09-13 | 北京电子工程总体研究所 | 一种基于太阳敏感器的卫星安全对日控制方法 |
Non-Patent Citations (2)
Title |
---|
于松柏,蔡娅雯,徐菁宇: "倾斜轨道光照特性及卫星总体设计对策", 《高分辨率对地观测学术联盟 第六届高分辨率对地观测学术年会论文集(上)》 * |
刘俊泽,杨小锋,朱光明: ""神舟"五号飞船轨道舱留轨长期管理期间的飞行控制", 《载人航天》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114485641A (zh) * | 2022-01-24 | 2022-05-13 | 武汉梦芯科技有限公司 | 一种基于惯导卫导方位融合的姿态解算方法及装置 |
CN114485641B (zh) * | 2022-01-24 | 2024-03-26 | 武汉梦芯科技有限公司 | 一种基于惯导卫导方位融合的姿态解算方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN111776251B (zh) | 2022-01-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106155074B (zh) | 一种保证星地链路的卫星三轴对日定向控制方法 | |
CN111897357B (zh) | 一种卫星对地扫描的姿态跟踪控制方法 | |
CN110174899A (zh) | 一种基于敏捷卫星的高精度成像姿态指向控制方法 | |
CN104848860B (zh) | 一种敏捷卫星成像过程姿态机动规划方法 | |
CN108508918B (zh) | 一种静轨遥感卫星数传天线高精度实时对地指向控制方法 | |
US5080307A (en) | Spacecraft earth-pointing attitude acquisition method | |
US5508932A (en) | Method and device for earth acquisition using the pole star for a three-axis stabilized satellite in a low inclination orbit | |
US8131409B2 (en) | Gyroless transfer orbit sun acquisition using only wing current measurement feedback | |
US7823836B2 (en) | Optimal sun safe attitude for satellite ground tracking | |
CN107380485B (zh) | 一种微小卫星大面阵广域多模凝视成像控制方法 | |
JP2625336B2 (ja) | 太陽センサを使用して人工衛星の太陽翼を制御する方法および装置 | |
JPH06510499A (ja) | 地球を指向する3軸安定化衛星および付属する太陽と地球を捕捉する方法 | |
CN109782787B (zh) | 一种太阳光压辅助下欠驱动航天器姿态的双模mpc控制方法 | |
CN110803304B (zh) | 一种卫星姿态控制*** | |
CN111007865B (zh) | 以对日指向偏差为约束的卫星平稳对地定向方法 | |
CN111776251B (zh) | 卫星的长期对日跟踪控制方法 | |
CN112550767A (zh) | 一种卫星偏航导引下的飞轮组动量管理方法 | |
CN110162069B (zh) | 一种近地轨道航天器阳光反射凝视期望姿态解析求解方法 | |
CN108427429B (zh) | 一种考虑动态指向约束的航天器视轴机动控制方法 | |
CN110955255B (zh) | 基于cmg的高精度轨控姿态维持方法、***及介质 | |
CN110697085B (zh) | 一种双sgcmg与磁力矩器组合的卫星控制方法 | |
CN110502038B (zh) | 一种机动过程中天线预置的高稳定度控制方法 | |
JP2021011257A5 (zh) | ||
US7797085B2 (en) | Process for determining the position of a spacecraft with the aid of a directional vector and a total angular momentum measurement | |
CN115622594A (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 |