CN105819004B - 一种卫星的太阳翼控制方法、系统及卫星 - Google Patents
一种卫星的太阳翼控制方法、系统及卫星 Download PDFInfo
- Publication number
- CN105819004B CN105819004B CN201610251520.9A CN201610251520A CN105819004B CN 105819004 B CN105819004 B CN 105819004B CN 201610251520 A CN201610251520 A CN 201610251520A CN 105819004 B CN105819004 B CN 105819004B
- Authority
- CN
- China
- Prior art keywords
- angle
- satellite
- solar wing
- sun
- solar
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000005286 illumination Methods 0.000 claims description 24
- 108091092878 Microsatellite Proteins 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 13
- 230000007246 mechanism Effects 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 241000197727 Euscorpius alpha Species 0.000 description 2
- 238000012369 In process control Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010965 in-process control Methods 0.000 description 2
- 230000000903 blocking effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000003306 harvesting Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
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/42—Arrangements or adaptations of power supply systems
- B64G1/44—Arrangements or adaptations of power supply systems using radiation, e.g. deployable solar arrays
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Remote Sensing (AREA)
- Aviation & Aerospace Engineering (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供一种卫星的太阳翼控制方法、系统及卫星。所述卫星的太阳翼控制方法包括：获取所述卫星的太阳翼的最大转动角度，所述卫星的地球遮挡角；当设定的控制周期到达时，获取当前时间点t的轨道坐标系下的归一化太阳矢量，计算当前时间点t的太阳翼主轴指向太阳的转动角；确定所述太阳翼主轴指向太阳时的转动角落入的区域；根据所述太阳翼主轴指向太阳时的转动角落入的区域求解太阳翼的目标姿态角；所述目标姿态角包括偏航角、滚转角以及俯仰角，所述偏航角恒为0。本发明实现尽最大可能提高能源获取效率。
Description
技术领域
本发明涉及一种航天技术领域，特别是涉及一种卫星的太阳翼控制方法、系统及卫星。
背景技术
现代微小卫星具有重量轻、体积小等特点，逐步广泛应用于遥感、通信、空间科学试验等诸多领域，能够适用于的载荷也从以前的小功率逐步提高为大功率。然而微小卫星要想满足大功率载荷应用，就必须实现获取大功率能源获取。微小卫星的太阳翼在选取相同的单体太阳电池片情况下，主要能够从两方面获取大功率能源。一方面需要增大太阳翼的面积提高有效面积，另一方面需要通过调节太阳翼姿态指向太阳提高供电效率。然而，微小卫星由于自身体积和重量受限，通过太阳翼获取大功率能源时可能面临很多困难。比如：1)若太阳翼固联安装不对日跟踪指向，则存在着某段时间内无法获取能源问题(低倾角轨道时更明显)；2)安装相对更大的太阳翼在对日跟踪指向时会对微小卫星的星敏、天线、星间通信等单机造成遮挡；3)太阳翼运动时候会对卫星产生干扰力矩，而微小卫星通常安装较小型飞轮克服干扰力矩能力较差；4)微小卫星研制成本受限，其太阳翼跟踪指向机构较难采用可360度连续转向的滑环机构；5)微小卫星特别是由其构成的星座网络卫星数目较多，地面运控支援较少，需要星上自主完成不同太阳受照情况下的能源获取。
鉴于此，如何找到一种适用于微小卫星获得大功率能源的太阳翼控制方案就成了本领域技术人员亟待解决的问题。
发明内容
鉴于以上所述现有技术的缺点，本发明的目的在于提供一种卫星的太阳翼控制方法、系统及卫星，用于解决现有技术中卫星由于自身体积和重量受限，通过太阳翼获取大功率能源时面临很多困难的问题。
为实现上述目的及其他相关目的，本发明提供一种卫星的太阳翼控制方法，所述卫星的太阳翼控制方法包括：获取所述卫星的太阳翼的最大转动角度αmax，所述卫星的地球遮挡角γ；当设定的控制周期Δt到达时，获取当前时间点t的轨道坐标系下的归一化太阳矢量Rso(t)＝[xso(t),yso(t),zso(t)]，计算当前时间点t的太阳翼主轴指向太阳的转动角αd(t)＝arccos(zso(t))；确定所述太阳翼主轴指向太阳时的转动角αd(t)落入的区域，所述区域包括：αd(t)≤αmax、αmax<αd(t)<π-γ，αd(t)≥π-γ中的任一个；根据所述太阳翼主轴指向太阳时的转动角αd(t)落入的区域求解太阳翼的目标姿态角；所述目标姿态角包括偏航角、滚转角φd(t)以及俯仰角θd(t)，所述偏航角恒为0。
可选地，所述轨道坐标系包括VVLH坐标系。
可选地，当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≤αmax时，则当前卫星处于光照区且太阳翼主轴指向太阳的转动角αd(t)在太阳翼的最大转动角度αmax内，此时按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
可选地，当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αmax<αd(t)<π-γ时，则当前卫星处于光照区且太阳翼主轴指向太阳的转动角αd(t)在太阳翼的最大转动角度αmax以外，按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
其中，
可选地，当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≥π-γ时，判断是否同时满足αmax<|αd(t-Δt)|<π-γ，如果满足，则t0＝t；此时对姿态导引律参数进行修正计算：|β|＝arcsin(|yso(t)|)；T为卫星的轨道周期；aθ＝θd(t0)； aφ＝φd(t0)；
可选地，当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≥π-γ时，按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
可选地，所述卫星的太阳翼控制方法还包括：运用PID控制器实时获与所述目标姿态角对应的驱动信号。
可选地，所述PID控制器采用的PID控制算法包括：
其中，Oθ(t)和Oφ(t)分别为对应θ、φ角度的控制驱动信号；Kpθ、Kiθ、Kdθ、Kpφ、Kiφ、Kdφ为常规PID算法的控制参数，θ(t)、φ(t)为实时测量得到的θ、φ角度值。
可选地，所述卫星的太阳翼控制方法还包括：根据所述驱动信号驱动所述卫星的二维转动机构调整所述太阳翼的主轴指向。
本发明还提供一种卫星的太阳翼控制系统，所述卫星的太阳翼控制系统包括：基本参数获取模块，用于获取所述卫星的太阳翼的最大转动角度αmax，所述卫星的地球遮挡角γ；转动角计算模块，用于当设定的控制周期Δt到达时，获取当前时间点t的轨道坐标系下的归一化太阳矢量Rso(t)＝[xso(t),yso(t),zso(t)]，计算当前时间点t的太阳翼主轴指向太阳的转动角αd(t)＝arccos(zso(t))；确定所述太阳翼主轴指向太阳时的转动角αd(t)落入的区域，所述区域包括：αd(t)≤αmax、αmax<αd(t)<π-γ，αd(t)≥π-γ中的任一个；目标姿态角计算模块，用于根据所述太阳翼主轴指向太阳时的转动角αd(t)落入的区域求解太阳翼的目标姿态角；所述目标姿态角包括偏航角、滚转角φd(t)以及俯仰角θd(t)，所述偏航角恒为0。
可选地，当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≤αmax时，则当前卫星处于光照区且太阳翼主轴指向太阳的转动角αd(t)在太阳翼的最大转动角度αmax内，此时按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
可选地，当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αmax<αd(t)<π-γ时，则当前卫星处于光照区且太阳翼主轴指向太阳的转动角αd(t)在太阳翼的最大转动角度αmax以外，按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：其中，
可选地，当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≥π-γ时，判断是否同时满足αmax<|αd(t-Δt)|<π-γ，如果满足，则t0＝t；此时对姿态导引律参数进行修正计算：|β|＝arcsin(|yso(t)|)；T为卫星的轨道周期；aθ＝θd(t0)； aφ＝φd(t0)；
可选地，当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≥π-γ时，按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
可选地，所述卫星的太阳翼控制系统还包括驱动信号获取模块，用于运用PID控制器获与所述目标姿态角对应的驱动信号。
可选地，所述PID控制器采用的PID控制算法包括：
其中，Oθ(t)和Oφ(t)分别为对应θ、φ角度的控制驱动信号；Kpθ、Kiθ、Kdθ、Kpφ、Kiφ、Kdφ为常规PID算法的控制参数，θ(t)、φ(t)为实时测量得到的θ、φ角度值。
可选地，所述轨道坐标系包括VVLH坐标系。可选地，所述卫星的太阳翼控制系统还包括太阳翼调整模块，用于根据所述驱动信号驱动所述卫星的二维转动机构以调整所述太阳翼的主轴指向。
可选地，所述二维转动机构安装于所述卫星的顶部、底部或内部。
本发明还提供一种卫星，所述卫星包括如上所述的卫星的太阳翼控制系统。
可选地，所述卫星为微小卫星。
如上所述，本发明的一种卫星的太阳翼控制方法、系统及卫星，具有以下有益效果：能够实现卫星的太阳翼在既不对卫星正常工作造成遮挡，也不对卫星平台产生较大的干扰力矩情况下自主二自由度跟踪太阳，从而实现尽最大可能提高能源获取效率，确保微小卫星的大功率能源获取。
附图说明
图1显示为本发明的卫星的太阳翼控制方法的一实施例的流程示意图。
图2显示为本发明的卫星的太阳翼控制方法的一实施例的地球遮挡角示意图。
图3显示为本发明的卫星的太阳翼控制方法的一实施例的控制原理示意图。
图4显示为本发明的卫星的太阳翼控制方法的另一实施例的流程示意图。
图5显示为本发明的卫星的太阳翼控制系统的一实施例的模块示意图。
元件标号说明
1 卫星的太阳翼控制系统
11 基本参数获取模块
12 转动角计算模块
13 目标姿态角计算模块
22 卫星
23 最大转动角形成的锥面
24 地球遮挡角形成的锥面
25 太阳绕卫星轨道面法线形成的锥面
26 卫星轨道面
S1～S3 步骤
具体实施方式
以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。
需要说明的是，本实施例中所提供的图示仅以示意方式说明本发明的基本构想，遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制，其实际实施时各组件的型态、数量及比例可为一种随意的改变，且其组件布局型态也可能更为复杂。
本发明提供一种卫星的太阳翼控制方法。所述卫星的太阳翼控制方法用于调整卫星的太阳翼的指向，让卫星的太阳翼持续跟踪指向太阳，以得到最大程度的太阳光照，从而获得大功率的太阳能。在一个实施例中，如图1所示，所述卫星的太阳翼控制方法包括：
步骤S1，获取所述卫星的太阳翼的最大转动角度αmax，所述卫星的地球遮挡角γ。太阳翼最大转动角αmax指太阳翼主轴的最大转角，即太阳翼主轴只能在以卫星矢径为轴的αmax半锥角的圆锥内转动。地球遮挡角γ定义如图2所示，即卫星垂直与地球球心的矢量(天地矢量)与卫星到地球切线的连线矢量的夹角。计算公式为
其中Re为地球半径，h为卫星轨道高度。当卫星与太阳的连线位于地球遮挡角的范围时，太阳将被地球所遮挡，此时，卫星处于阴影区，即不能得到太阳的光照。
步骤S2，当设定的控制周期Δt到达时，获取当前时间点t的轨道坐标系下的归一化太阳矢量Rso(t)＝[xso(t),yso(t),zso(t)]，计算当前时间点t的太阳翼主轴指向太阳的转动角αd(t)＝arccos(zso(t))；确定所述太阳翼主轴指向太阳时的转动角αd(t)落入的区域，所述区域包括：αd(t)≤αmax、αmax<αd(t)<π-γ，αd(t)≥π-γ中的任一个。控制周期Δt为本太阳翼控制方法中的控制周期，每经过一个控制周期Δt，太阳翼都将计算太阳翼主轴指向太阳时的转动角αd(t)，并求解太阳翼的目标姿态角，进而控制太阳翼跟踪太阳转动，使得太阳翼得到最大程度的太阳光照，获得最大的太阳能。轨道坐标系下的归一化太阳矢量Rso(t)＝[xso(t),yso(t),zso(t)]为当前时间点t的太阳在轨道坐标系下的坐标，xso(t)为轨道坐标系下的归一化太阳矢量中的x方向的值，yso(t)为轨道坐标系下的归一化太阳矢量中的y方向的值，zso(t)为轨道坐标系下的归一化太阳矢量中的z方向的值。控制原理如图3所示，其中，以卫星轨道面26为参照坐标，卫星22太阳翼的最大转动角形成的锥面23，地球遮挡角形成的锥面24，太阳绕卫星轨道面法线形成的锥面25。当太阳的位置位于太阳翼的最大转动角形成的锥面23之上时，表明卫星的转动角在最大转动角的范围内，即αd(t)≤αmax。当太阳的位置位于地球遮挡角形成的锥面24之下时，表明卫星处于阴影区，此时，由于地球的遮挡，卫星得不到太阳光，即αd(t)≥π-γ。当太阳的位置位于太阳翼的最大转动角形成的锥面23与地球遮挡角形成的锥面24之间时，表明卫星处于光照区，但卫星的转动角超出了最大转动角的范围，即αmax<αd(t)<π-γ。在一个实施例中，所述轨道坐标系包括VVLH轨道坐标系。VVLH(Vehicle Velocity Local Horizontal)轨道坐标系，原点在飞行器质心；Z轴指向地心；X轴指在轨道面内飞行方向，垂直Z轴；Y按右手法则确定。
步骤S3，根据所述太阳翼主轴指向太阳时的转动角αd(t)落入的区域求解太阳翼的目标姿态角；所述目标姿态角包括偏航角、滚转角φd(t)以及俯仰角θd(t)，所述偏航角恒为0。通常，确定导弹(或飞机)在空间中的方向需要用三个角度，分别为偏航角、俯仰角和滚转角，这三个角度通常称为欧拉角，或弹体的姿态角。滚转角定义为弹体的Oy轴与惯性坐标系的夹角。偏航角即为导弹(或飞机)的飞行偏离目标方向(沿oy坐标轴)的角度。俯仰角就是导弹(或飞机)相对于惯性坐标系的XOY平面“俯仰”的角度。针对所述转动角αd(t)落入的区域不同，计算所述目标姿态角的算法也不同。
在具体的实施例中，当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≤αmax时，则当前卫星处于光照区且太阳翼主轴指向太阳的转动角αd(t)在太阳翼的最大转动角度αmax内，此时按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αmax<αd(t)<π-γ时，则当前卫星处于光照区且太阳翼主轴指向太阳的转动角αd(t)在太阳翼的最大转动角度αmax以外，按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
其中，
当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≥π-γ时，由于此时卫星在阴影区内，太阳翼主轴对日定向没有意义，转回太快又会出现干扰力矩较大等问题，故采用平滑回转导引律生成方法，将主轴指向从进阴影时刻的指向平滑转到出阴影时刻的指向。具体地，按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
其中，姿态导引律参数包括aθ,bθ,cθ,dθ以及aφ,bφ,cφ。当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≥π-γ时，判断是否同时满足αmax<|αd(t-Δt)|<π-γ，如果满足，则将此时的当前时间点作为本周期T中卫星由光照区进入阴影区的分界点t0＝t；此时对姿态导引律参数进行修正计算。在本周期T内，计算得到的姿态导引律参数不再修改，直接用于计算目标姿态角。阴影时间计算方法如下：
cosγ＝cosα·cos|β|
其中γ为地球遮挡角，β为轨道β角，即太阳矢量与轨道面的夹角，T为轨道周期，ts为阴影区时长。参数调节边界条件如下：
其中，采用多项式拟合下：
fθ(t)＝aθ+bθ×(t-t0)+cθ×(t-t0)2+dθ×(t-t0)3t∈[t0,ts]
fφ(t)＝aφ+bφ×(t-t0)+cφ×(t-t0)2t∈[t0,ts]
其中，|β|＝arcsin(|yso(t)|)；T为卫星的轨道周期，Δt为控制周期。
修改得到的姿态导引律参数如下：aθ＝θd(t0)； aφ＝φd(t0)；
在一个实施例中，所述卫星的太阳翼控制方法还包括：运用PID控制器实时获与所述目标姿态角对应的驱动信号。PID控制器(亦称PID调节器)是在过程控制中，按偏差的比例(P)、积分(I)和微分(D)进行控制的一种自动控制器。本发明中的PID控制器包括二自由度PID控制，就是使目标值跟踪特性为最优的PID参数和使外扰抑制特性最优的PID参数，能分别独立地进行整定，使两特性同时达到最优。这样，既提高了控制系统设计的自由度，又改善了系统的品质。在一个实施例中，所述PID控制器采用的PID控制算法包括：
其中，Oθ(t)和Oφ(t)分别为对应θ、φ角度的控制驱动信号；Kpθ、Kiθ、Kdθ、Kpφ、Kiφ、Kdφ为常规PID算法的控制参数，θ(t)、φ(t)为实时测量得到的θ、φ角度值。
在一个实施例中，所述卫星的太阳翼控制方法还包括：根据所述驱动信号驱动所述卫星的二维转动机构调整所述太阳翼的主轴指向。让卫星的太阳翼持续跟踪指向太阳，以得到最大程度的太阳光照，从而获得大功率的太阳能。
在一个实施例中，对所述姿态导引律参数的修正可以先行执行，如图4所示，所述卫星的太阳翼控制方法包括：
1，先确定是否处于阴影光照分界点，太阳翼主轴指向太阳时的转动角αd(t)是否同时满足
2，当太阳翼主轴指向太阳时的转动角αd(t)是否同时满足时，修改导引律参数。具体包括：赋值分界点时刻t0＝t，计算轨道β角|β|＝arcsin(|yso(t)|)，计算阴影区时间进而根据参数调节边界条件得到得到的姿态导引律参数如下：aθ＝θd(t0)； aφ＝φd(t0)；
3，判别导引模式，即确定所述太阳翼主轴指向太阳时的转动角αd(t)落入的区域，所述区域包括：αd(t)≤αmax、αmax<αd(t)<π-γ，αd(t)≥π-γ中的任一个。
4，当αd(t)≤αmax时，采用对日导引律生成目标姿态角。
5，当αmax<αd(t)<π-γ时，采用准对日导引律生成目标姿态角：
其中，
6，当αd(t)≥π-γ时，采用平滑回转导引律生成目标姿态角：
7，运用PID控制器实时获与所述目标姿态角对应的驱动信号，并驱动太阳翼转动。
本发明还提供一种卫星的太阳翼控制系统，所述卫星的太阳翼控制系统可以采用如上所述的卫星的太阳翼控制方法工作。在一个实施例中，如图5所示，所述卫星的太阳翼控制系统1包括基本参数获取模块11，转动角计算模块12以及目标姿态角计算模块13，其中：
基本参数获取模块11用于获取所述卫星的太阳翼的最大转动角度αmax，所述卫星的地球遮挡角γ。太阳翼最大转动角αmax指太阳翼主轴的最大转角，即太阳翼主轴只能在以卫星矢径为轴的αmax半锥角的圆锥内转动。地球遮挡角γ定义如图2所示，即卫星垂直与地球球心的矢量(天地矢量)与卫星到地球切线的连线矢量的夹角。计算公式为
其中Re为地球半径，h为卫星轨道高度。当卫星与太阳的连线位于地球遮挡角的范围时，太阳将被地球所遮挡，此时，卫星处于阴影区，即不能得到太阳的光照。
转动角计算模块12与基本参数获取模块11相连，用于当设定的控制周期Δt到达时，获取当前时间点t的轨道坐标系下的归一化太阳矢量Rso(t)＝[xso(t),yso(t),zso(t)]，计算当前时间点t的太阳翼主轴指向太阳的转动角αd(t)＝arccos(zso(t))；确定所述太阳翼主轴指向太阳时的转动角αd(t)落入的区域，所述区域包括：αd(t)≤αmax、αmax<αd(t)<π-γ，αd(t)≥π-γ中的任一个。控制周期Δt为本太阳翼控制方法中的控制周期，每经过一个控制周期Δt，太阳翼都将计算太阳翼主轴指向太阳时的转动角αd(t)，并求解太阳翼的目标姿态角，进而控制太阳翼跟踪太阳转动，使得太阳翼得到最大程度的太阳光照，获得最大的太阳能。轨道坐标系下的归一化太阳矢量Rso(t)＝[xso(t),yso(t),zso(t)]为当前时间点t的太阳在轨道坐标系下的坐标，xso(t)为轨道坐标系下的归一化太阳矢量中的x方向的值，yso(t)为轨道坐标系下的归一化太阳矢量中的y方向的值，zso(t)为轨道坐标系下的归一化太阳矢量中的z方向的值。控制原理如图3所示，其中，以卫星轨道面26为参照坐标，卫星22太阳翼的最大转动角形成的锥面23，地球遮挡角形成的锥面24，太阳绕卫星轨道面法线形成的锥面25。当太阳的位置位于太阳翼的最大转动角形成的锥面23之上时，表明卫星的转动角在最大转动角的范围内，即αd(t)≤αmax。当太阳的位置位于地球遮挡角形成的锥面24之下时，表明卫星处于阴影区，此时，由于地球的遮挡，卫星得不到太阳光，即αd(t)≥π-γ。当太阳的位置位于太阳翼的最大转动角形成的锥面23与地球遮挡角形成的锥面24之间时，表明卫星处于光照区，但卫星的转动角超出了最大转动角的范围，即αmax<αd(t)<π-γ。在一个实施例中，所述轨道坐标系包括VVLH轨道坐标系。VVLH(Vehicle Velocity LocalHorizontal)轨道坐标系，原点在飞行器质心；Z轴指向地心；X轴指在轨道面内飞行方向，垂直Z轴；Y按右手法则确定。
目标姿态角计算模块13与转动角计算模块12相连，用于根据所述太阳翼主轴指向太阳时的转动角αd(t)落入的区域求解太阳翼的目标姿态角；所述目标姿态角包括偏航角、滚转角φd(t)以及俯仰角θd(t)，所述偏航角恒为0。
在具体的实施例中，当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≤αmax时，则当前卫星处于光照区且太阳翼主轴指向太阳的转动角αd(t)在太阳翼的最大转动角度αmax内，此时按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αmax<αd(t)<π-γ时，则当前卫星处于光照区且太阳翼主轴指向太阳的转动角αd(t)在太阳翼的最大转动角度αmax以外，按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
其中，
当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≥π-γ时，由于此时卫星在阴影区内，太阳翼主轴对日定向没有意义，转回太快又会出现干扰力矩较大等问题，故采用平滑回转导引律生成方法，将主轴指向从进阴影时刻的指向平滑转到出阴影时刻的指向。具体地，按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
其中，姿态导引律参数包括aθ,bθ,cθ,dθ以及aφ,bφ,cφ。当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≥π-γ时，判断是否同时满足αmax<|αd(t-Δt)|<π-γ，如果满足，则将此时的当前时间点作为本周期T中卫星由光照区进入阴影区的分界点t0＝t；此时对姿态导引律参数进行修正计算。在本周期T内，计算得到的姿态导引律参数不再修改，直接用于计算目标姿态角。阴影时间计算方法如下：
cosγ＝cosα·cos|β|
其中γ为地球遮挡角，β为轨道β角，即太阳矢量与轨道面的夹角，T为轨道周期，ts为阴影区时长。参数调节边界条件如下：
其中，采用多项式拟合下：
fθ(t)＝aθ+bθ×(t-t0)+cθ×(t-t0)2+dθ×(t-t0)3 t∈[t0,ts]
fφ(t)＝aφ+bφ×(t-t0)+cφ×(t-t0)2 t∈[t0,ts]
其中，|β|＝arcsin(|yso(t)|)；T为卫星的轨道周期，Δt为控制周期。
修改得到的姿态导引律参数如下：aθ＝θd(t0)； aφ＝φd(t0)；
在一个实施例中，所述卫星的太阳翼控制系统还包括驱动信号获取模块，用于运用PID控制器获得与所述目标姿态角对应的驱动信号。PID控制器(亦称PID调节器)是在过程控制中，按偏差的比例(P)、积分(I)和微分(D)进行控制的一种自动控制器。本发明中的PID控制器包括二自由度PID控制，就是使目标值跟踪特性为最优的PID参数和使外扰抑制特性最优的PID参数，能分别独立地进行整定，使两特性同时达到最优。这样，既提高了控制系统设计的自由度，又改善了系统的品质。在一个实施例中，所述PID控制器采用的PID控制算法包括：
其中，Oθ(t)和Oφ(t)分别为对应θ、φ角度的控制驱动信号；Kpθ、Kiθ、Kdθ、Kpφ、Kiφ、Kdφ为常规PID算法的控制参数，θ(t)、φ(t)为实时测量得到的θ、φ角度值。
在一个实施例中，所述卫星的太阳翼控制系统还包括太阳翼调整模块：用于根据所述驱动信号驱动所述卫星的二维转动机构调整所述太阳翼的主轴指向。让卫星的太阳翼持续跟踪指向太阳，以得到最大程度的太阳光照，从而获得大功率的太阳能。所述二维转动机构可以安装于所述卫星的顶部、底部或内部。
本发明还提供一种卫星，所述卫星包括如上所述的卫星的太阳翼控制系统。所述卫星为微小卫星。重量在1000千克以下的人造卫星统称为"微小卫星"，微小卫星具有种种优势和特点：高新技术含量高，研制周期短一年左右，研制经费低人民币数千万元，且可以进一步组网形成分布式星座形“虚拟大卫星”，以代替目前的大卫星。
综上所述，本发明的一种卫星的太阳翼控制方法、系统及卫星能够实现卫星的太阳翼在既不对卫星正常工作造成遮挡，也不对卫星平台产生较大的干扰力矩情况下自主二自由度跟踪太阳，从而实现尽最大可能提高能源获取效率，确保微小卫星的大功率能源获取。所以，本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。
上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。
Claims (19)
1.一种卫星的太阳翼控制方法，其特征在于，所述卫星的太阳翼控制方法包括：
获取所述卫星的太阳翼的最大转动角度αmax，所述卫星的地球遮挡角γ；
当设定的控制周期Δt到达时，获取当前时间点t的轨道坐标系下的归一化太阳矢量Rso(t)＝[xso(t),yso(t),zso(t)]，计算当前时间点t的太阳翼主轴指向太阳的转动角αd(t)＝arccos(zso(t))；确定所述太阳翼主轴指向太阳时的转动角αd(t)落入的区域，所述区域包括：αd(t)≤αmax、αmax＜αd(t)＜π-γ，αd(t)≥π-γ中的任一个；
根据所述太阳翼主轴指向太阳时的转动角αd(t)落入的区域求解太阳翼的目标姿态角；所述目标姿态角包括偏航角、滚转角φd(t)以及俯仰角θd(t)，所述偏航角恒为0。
2.根据权利要求1所述的卫星的太阳翼控制方法，其特征在于：当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≤αmax时，则当前卫星处于光照区且太阳翼主轴指向太阳的转动角αd(t)在太阳翼的最大转动角度αmax内，此时按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
3.根据权利要求1所述的卫星的太阳翼控制方法，其特征在于：当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αmax＜αd(t)＜π-γ时，则当前卫星处于光照区且太阳翼主轴指向太阳的转动角αd(t)在太阳翼的最大转动角度αmax以外，按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
其中，
4.根据权利要求1所述的卫星的太阳翼控制方法，其特征在于：当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≥π-γ时，判断是否同时满足αmax＜|αd(t-Δt)|＜π-γ，如果满足，则t0＝t；此时对姿态导引律参数aθ,bθ,cθ,dθ,aφ,bφ,cφ进行修正计算：|β|＝arcsin(|yso(t)|)；T为卫星的轨道周期；aθ＝θd(t0)； aφ＝φd(t0)；
5.根据权利要求4所述的卫星的太阳翼控制方法，其特征在于：当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≥π-γ时，按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
6.根据权利要求1所述的卫星的太阳翼控制方法，其特征在于：所述卫星的太阳翼控制方法还包括：运用PID控制器实时获取与所述目标姿态角对应的驱动信号。
7.根据权利要求6所述的卫星的太阳翼控制方法，其特征在于：所述PID控制器采用的PID控制算法包括：
其中，Oθ(t)和Oφ(t)分别为对应θd(t)、φd(t)角度的控制驱动信号；Kpθ、Kiθ、Kdθ、Kpφ、Kiφ、Kdφ为常规PID算法的控制参数，θ(t)、φ(t)为实时测量得到的θd(t)、φd(t)角度值。
8.根据权利要求6所述的卫星的太阳翼控制方法，其特征在于：所述卫星的太阳翼控制方法还包括：根据所述驱动信号驱动所述卫星的二维转动机构调整所述太阳翼的主轴指向。
9.一种卫星的太阳翼控制系统，其特征在于：所述卫星的太阳翼控制系统包括：
基本参数获取模块，用于获取所述卫星的太阳翼的最大转动角度αmax，所述卫星的地球遮挡角γ；
转动角计算模块，用于当设定的控制周期Δt到达时，获取当前时间点t的轨道坐标系下的归一化太阳矢量Rso(t)＝[xso(t),yso(t),zso(t)]，计算当前时间点t的太阳翼主轴指向太阳的转动角αd(t)＝arccos(zso(t))；确定所述太阳翼主轴指向太阳时的转动角αd(t)落入的区域，所述区域包括：αd(t)≤αmax、αmax＜αd(t)＜π-γ，αd(t)≥π-γ中的任一个；
目标姿态角计算模块，用于根据所述太阳翼主轴指向太阳时的转动角αd(t)落入的区域求解太阳翼的目标姿态角；所述目标姿态角包括偏航角、滚转角φd(t)以及俯仰角θd(t)，所述偏航角恒为0。
10.根据权利要求9所述的卫星的太阳翼控制系统，其特征在于：当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≤αmax时，则当前卫星处于光照区且太阳翼主轴指向太阳的转动角αd(t)在太阳翼的最大转动角度αmax内，此时按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
11.根据权利要求9所述的卫星的太阳翼控制系统，其特征在于：当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αmax＜αd(t)＜π-γ时，则当前卫星处于光照区且太阳翼主轴指向太阳的转动角αd(t)在太阳翼的最大转动角度αmax以外，按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
其中，
12.根据权利要求9所述的卫星的太阳翼控制系统，其特征在于：当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≥π-γ时，判断是否同时满足αmax＜|αd(t-Δt)|＜π-γ，如果满足，则t0＝t；此时对姿态导引律参数aθ,bθ,cθ,dθ,aφ,bφ,cφ进行修正计算：|β|＝arcsin(|yso(t)|)；T为卫星的轨道周期；aθ＝θd(t0)； aφ＝φd(t0)；
13.根据权利要求12所述的卫星的太阳翼控制系统，其特征在于：当所述太阳翼主轴指向太阳时的转动角αd(t)落入区域αd(t)≥π-γ时，按照以下公式求解太阳翼目标姿态角的滚转角φd(t)以及俯仰角θd(t)：
14.根据权利要求9所述的卫星的太阳翼控制系统，其特征在于：所述卫星的太阳翼控制系统还包括驱动信号获取模块，用于运用PID控制器获取与所述目标姿态角对应的驱动信号。
15.根据权利要求14所述的卫星的太阳翼控制系统，其特征在于：所述PID控制器采用的PID控制算法包括：
其中，Oθ(t)和Oφ(t)分别为对应θd(t)、φd(t)角度的控制驱动信号；Kpθ、Kiθ、Kdθ、Kpφ、Kiφ、Kdφ为常规PID算法的控制参数，θ(t)、φ(t)为实时测量得到的θd(t)、φd(t)角度值。
16.根据权利要求14所述的卫星的太阳翼控制系统，其特征在于：所述卫星的太阳翼控制系统还包括太阳翼调整模块，用于根据所述驱动信号驱动所述卫星的二维转动机构以调整所述太阳翼的主轴指向。
17.根据权利要求16所述的卫星的太阳翼控制系统，其特征在于：所述二维转动机构安装于所述卫星的顶部、底部或内部。
18.一种卫星，其特征在于：所述卫星包括如权利要求9至17中的任一权利要求所述的卫星的太阳翼控制系统。
19.根据权利要求18所述的卫星，其特征在于：所述卫星为微小卫星。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610251520.9A CN105819004B (zh) | 2016-04-21 | 2016-04-21 | 一种卫星的太阳翼控制方法、系统及卫星 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610251520.9A CN105819004B (zh) | 2016-04-21 | 2016-04-21 | 一种卫星的太阳翼控制方法、系统及卫星 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105819004A CN105819004A (zh) | 2016-08-03 |
| CN105819004B true CN105819004B (zh) | 2018-01-16 |
Family
ID=56527233
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610251520.9A Expired - Fee Related CN105819004B (zh) | 2016-04-21 | 2016-04-21 | 一种卫星的太阳翼控制方法、系统及卫星 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105819004B (zh) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106570614A (zh) * | 2016-10-14 | 2017-04-19 | 上海微小卫星工程中心 | 星上自主分布式任务调度方法 |
| CN106843249B (zh) * | 2017-01-24 | 2019-12-17 | 上海航天控制技术研究所 | 一种二维导引姿态控制方法 |
| CN106970573B (zh) * | 2017-05-22 | 2019-03-15 | 上海航天控制技术研究所 | 一种航天器太阳翼永磁同步电机驱动控制系统 |
| CN107839903B (zh) * | 2017-09-25 | 2020-07-14 | 上海卫星工程研究所 | 单翼卫星转移轨道段帆板放气时间估计方法 |
| CN108657467B (zh) * | 2018-05-16 | 2019-05-21 | 清华大学 | 一种利用虚拟太阳矢量的航天器偏航机动控制方法及系统 |
| CN109032158B (zh) * | 2018-07-24 | 2021-07-09 | 北京控制工程研究所 | 一种直线规划的头对日—尾对日交替连续偏航姿态控制方法 |
| CN109677637B (zh) * | 2019-02-22 | 2021-05-18 | 北京空间技术研制试验中心 | 基于光学测角相机的空间非合作目标跟踪指向方法 |
| CN110641741B (zh) * | 2019-10-23 | 2020-12-18 | 北京国电高科科技有限公司 | 双自由度太阳帆板控制方法及其控制系统 |
| CN110963087B (zh) * | 2019-11-11 | 2021-10-29 | 上海航天控制技术研究所 | 一种多模式复杂流程的空间站太阳翼的控制方法 |
| CN111717415B (zh) * | 2020-06-22 | 2022-05-17 | 中国科学院微小卫星创新研究院 | 一种星座卫星太阳电池阵对日跟踪方法 |
| CN112777003B (zh) * | 2020-12-31 | 2021-09-21 | 中国人民解放军63923部队 | 一种防止sada电弧放电对卫星造成损伤的方法及装置 |
| CN113110612B (zh) * | 2021-04-13 | 2023-11-28 | 深圳航天东方红卫星有限公司 | 一种双轴太阳电池阵驱动装置及其运动规划方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5697582A (en) * | 1991-05-07 | 1997-12-16 | Deutsche Aerospace Ag | Method of adjusting the position of satellites by means of solar pressure torques |
| US5794891A (en) * | 1995-01-06 | 1998-08-18 | Matra Marconi Space France | Method of controlling the attitude control for satellites on an orbit inclined relative to the equator |
| US6164597A (en) * | 1999-08-05 | 2000-12-26 | Space Systems/Loral, Inc. | Using reflected solar radiation pressure to provide stationkeeping maneuvering of a spacecraft |
| CN101858747A (zh) * | 2010-03-26 | 2010-10-13 | 航天东方红卫星有限公司 | 一种有效利用地球辐照能的卫星帆板对日定向目标姿态的解析确定方法 |
| CN103448922A (zh) * | 2013-08-08 | 2013-12-18 | 上海卫星工程研究所 | 静止轨道卫星太阳电池阵在轨被遮挡的分析计算方法 |
-
2016
- 2016-04-21 CN CN201610251520.9A patent/CN105819004B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5697582A (en) * | 1991-05-07 | 1997-12-16 | Deutsche Aerospace Ag | Method of adjusting the position of satellites by means of solar pressure torques |
| US5794891A (en) * | 1995-01-06 | 1998-08-18 | Matra Marconi Space France | Method of controlling the attitude control for satellites on an orbit inclined relative to the equator |
| US6164597A (en) * | 1999-08-05 | 2000-12-26 | Space Systems/Loral, Inc. | Using reflected solar radiation pressure to provide stationkeeping maneuvering of a spacecraft |
| CN101858747A (zh) * | 2010-03-26 | 2010-10-13 | 航天东方红卫星有限公司 | 一种有效利用地球辐照能的卫星帆板对日定向目标姿态的解析确定方法 |
| CN103448922A (zh) * | 2013-08-08 | 2013-12-18 | 上海卫星工程研究所 | 静止轨道卫星太阳电池阵在轨被遮挡的分析计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105819004A (zh) | 2016-08-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105819004B (zh) | 一种卫星的太阳翼控制方法、系统及卫星 | |
| CN107600464B (zh) | 利用星敏感器信息的飞轮控制捕获太阳及对日定向方法 | |
| CN105468009B (zh) | 应用于微小型飞行器的多动力融合飞控系统以及方法 | |
| CN105698762B (zh) | 一种单机航迹上基于不同时刻观测点的目标快速定位方法 | |
| CN1039301C (zh) | 卫星姿态控制的测试设备及其应用 | |
| CN104369877B (zh) | 一种深空探测器天线指向的设计方法 | |
| CN111609857A (zh) | 一种空间碎片定轨遍历观测方法及系统 | |
| CN105891821A (zh) | 一种机载下视测量目标的自动跟踪方法 | |
| JPH0133400B2 (zh) | ||
| US10386857B2 (en) | Sensor-centric path planning and control for robotic vehicles | |
| CN106483466B (zh) | 一种卫星入轨阶段太阳电池阵输出电流的估算方法 | |
| CN112572835B (zh) | 一种具有姿态切换的卫星在轨角动量管理及控制方法 | |
| CN103955224B (zh) | 一种用于相对运动视线跟踪的姿态控制方法 | |
| CN112198885B (zh) | 一种满足机动平台自主降落需求的无人机控制方法 | |
| CN110697086B (zh) | 以单控制力矩陀螺和两飞轮实现卫星三轴稳定控制的方法 | |
| CN111319791B (zh) | 一种适用于火星探测的姿态控制模式管理方法 | |
| CN110775302A (zh) | 一种基于太阳帆板输出电流信息的应急对日方法 | |
| CN110641741A (zh) | 双自由度太阳帆板控制方法及其控制系统 | |
| CN109657417A (zh) | 一种用于惯性空间观测卫星的能源保障设计方法 | |
| CN108427429B (zh) | 一种考虑动态指向约束的航天器视轴机动控制方法 | |
| CN109625329A (zh) | 一种基于离散喷气的飞轮角动量自主卸载方法 | |
| CN111007865A (zh) | 以对日指向偏差为约束的卫星平稳对地定向方法 | |
| CN111483618A (zh) | 一种适用于火星探测的多指向驱动任务自主管理控制方法 | |
| CN113772130A (zh) | 一种确定太阳电池阵法线矢量的方法 | |
| RU2414392C1 (ru) | Способ ориентации осей космического аппарата в солнечно-орбитальную систему координат |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180116 Termination date: 20180421 |