CN108549412A - 一种考虑太阳角变化率和控制因子权重可变的磁控太阳捕获方法 - Google Patents
一种考虑太阳角变化率和控制因子权重可变的磁控太阳捕获方法 Download PDFInfo
- Publication number
- CN108549412A CN108549412A CN201810304610.9A CN201810304610A CN108549412A CN 108549412 A CN108549412 A CN 108549412A CN 201810304610 A CN201810304610 A CN 201810304610A CN 108549412 A CN108549412 A CN 108549412A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- solar
- control
- angle
- vector
- 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 47
- 238000013016 damping Methods 0.000 claims description 40
- 238000004422 calculation algorithm Methods 0.000 claims description 16
- 238000000205 computational method Methods 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 8
- 239000000463 material Substances 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000012163 sequencing technique Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000006340 racemization Effects 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/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/28—Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect
- B64G1/283—Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect using reaction wheels
-
- 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)
- Remote Sensing (AREA)
- Aviation & Aerospace Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种考虑太阳角变化率和控制因子权重可变的磁控太阳捕获方法,包括:进行太阳角、太阳角速度以及磁矢量变化率的模的计算;进行太阳矢量计算有效性判断;当太阳矢量计算有效时,基于飞轮启旋偏置完成情况进行施加力矩方向计算;进行磁矩施加方向计算,并计算磁矩与力矩方向的夹角;进行磁矩与力矩方向的夹角是否在上下限之间的范围内判断;当磁矩与力矩的夹角处于上下限之间的范围内,采用速率加太阳角控制加太阳角速度补偿算法计算磁控输出;以及进行磁控输出磁矩限幅。
Description
技术领域
本发明涉及卫星姿态控制领域,特别涉及一种考虑太阳角变化率和控制因 子权重可变的磁控太阳捕获方法。
背景技术
卫星在某些工作情况约束下需要保持帆板对太阳持续定向的任务,如卫星 单机故障等姿控异常情况下进入安全模式,需要保持帆板对太阳定向以保障卫 星的能源安全。
通常情况下卫星进入安全模式下,一般采用推进控制或基于磁控控制的方 式保持对日定向。使用推进控制控制的卫星受限于推进可靠性、燃料等限制不 宜长时间进行对日定向控制,且对于部分卫星而言甚至没有配置推进组件。因 而,使用磁控保持对日定向不消耗工质,是一种适用于长时间、可靠性相对较 高且工程上易于实现的一种对日定向控制方式。
基于磁控的对日定向通常采用卫星以一定角速度自旋而对日定向、卫星反 作用轮启旋偏置或三轴纯磁控方式进行稳定对日定向。以卫星自旋+磁控方式 实现对日定向需要将卫星首先旋转起来,再进行磁控对日定向控制;基于三轴 纯磁控方式进行稳定对日定向的控制方式存在控制器参数设计难点和设计约 束;而基于卫星反作用轮启旋偏置+磁控的方式可以避免需要将启旋卫星的问 题,且较三轴纯磁控方式在工程上易于实现。
基于卫星反作用轮启旋偏置加磁控的方式,现阶段磁控算法存在长时间控 制收敛不稳定,导致卫星帆板太阳能电池面与太阳矢量夹角(该夹角简称太阳 角)起伏变化较大,另一方面磁控算法也存在稳定控制太阳角精度不理想,太 阳角控制精度低于20度,甚至40度,导致卫星太阳能供电受到影响。
发明内容
针对现有技术中的基于卫星反作用轮启旋偏置加磁控方式对日定向磁控 算法存在长时间控制收敛不稳定,控制太阳角精度不理想导致卫星太阳能供电 受到影响的问题,根据本发明的一个实施例,提供一种考虑太阳角变化率和控 制因子权重可变的磁控太阳捕获方法,包括:进行太阳角、太阳角速度以及磁 矢量变化率的模的计算;进行太阳矢量计算有效性判断;当太阳矢量计算有效 时,基于飞轮启旋偏置完成情况进行施加力矩方向计算;进行磁矩施加方向计 算,并计算磁矩与力矩方向的夹角;进行磁矩与力矩方向的夹角是否在上下限 之间的范围内判断;当磁矩与力矩的夹角处于上下限之间的范围内,采用速率 加太阳角控制加太阳角速度补偿算法计算磁控输出;以及进行磁控输出磁矩限 幅。
在本发明的一个实施例中,太阳角、太阳角速度以及磁矢量变化率的模的 计算的公式分别为:
其中gama为当前周期计算获得的太阳角,Lb是卫星 本体系中帆板产生太阳能面法线矢量,Sb是卫星本体系中的太阳矢量;
其中gamadot为太阳角速度,ΔT为姿控软件周期;
磁矢量变化率的模Bdotnorm计算分为两步,先进行磁场矢量变化速度 Bdot计算:
其中10ΔT为磁強计采样周期,Bb为本次进行 磁強计采样时的磁场矢量,Bb -为上一次进行磁強计采样时的磁场矢量,
再进行磁矢量变化率的模Bdotnorm计算:
在本发明的一个实施例中,进行太阳矢量计算有效性判断即判断太阳矢量 是否太敏无效或递推无效。
在本发明的一个实施例中,基于飞轮启旋偏置完成情况进行施加力矩方向 计算方法为:
若飞轮偏置完成,则计算方法为Tn=(Lb×Sb)×Lb;
若飞轮(轮Y)偏置未完成,则计算方法为Tn=Lb×Sb;以及
对力矩方向Tn进行单位化,单位化方法为其中Tn为力矩方向, Lb为帆板法线矢量,Sb是太阳矢量。
在本发明的一个实施例中,进行磁矩施加方向计算,并计算磁矩与力矩方 向的夹角的方法为:
磁矩施加方向由力矩方向计算,计算方法为Ang_TnBb=Tn×Bb,磁矩与力 矩方向的夹角的计算方法为其中Ang_TnBb是非归一化的, Bb为磁矩方向,Tn为力矩方向。
在本发明的一个实施例中,速率加太阳角控制加太阳角速度补偿算法的计 算方法为:
Pb=-k2_Bdot-(k6gama+k6k3_.gamadot).Ang_TnBb
其中Pb表示磁矩,Bdot为磁场矢量的变化率,k2为Bdot阻尼因子,gama 为太阳角,k6为太阳角控制因子,gamadot为太阳角速度k6k3为太阳角速度控 制因子,Ang_TnBb为磁矩施加方向。
在本发明的一个实施例中,在进行所述的采用速率加太阳角控制加太阳角 速度补偿算法计算磁控输出前,还包括判断卫星角速度是否大于第一阈值 BdotLimt(120nT)或者太阳角是否大于第二阈值GamaLimt。
在本发明的一个实施例中,如果卫星角速度大于第一阈值BdotLimt(120nT) 或者太阳角大于第二阈值GamaLimt,则采用高权重速率阻尼控制系数;如果 卫星角速度小于第一阈值BdotLimt(120nT)且太阳角小于第二阈值GamaLimt, 则采用高权重太阳角控制系数。
在本发明的一个实施例中,高权重速率阻尼控制系数为k2_=k2.k4,k3_= k3.k5,其中k4是速率阻尼修正控制系数,k5是太阳角速度修正控制系数,其值 为均大于1的正实数;所述的高权重太阳角控制系数为k2_=k2,k3_=k3。
在本发明的一个实施例中,当太阳矢量计算失效时,卫星仅进行速率阻尼 控制,并进行速率阻尼控制磁控输出的磁矩限幅。
在本发明的一个实施例中,当磁矩与力矩的夹角超出上下限之间的范围时, 卫星仅进行速率阻尼控制,并进行速率阻尼控制磁控输出的磁矩限幅。
在本发明的一个实施例中,速率阻尼控制的技术公式为Pb=-k1Bdot,其中 Pb表示磁矩,k1表示角速度阻尼因子,Bdot为磁场矢量的变化率。
在本发明的一个实施例中,该方法还包括在每次进行算法执行的同时,进 行飞轮启旋条件判断,若磁矢量变化率的模小于第三阈值(例如,200nT), 且当前太阳角小于等于第四阈值(例如,50度)时,飞轮启旋;否则飞轮转 速保持。
本发明提供一种考虑太阳角变化率和控制因子权重可变的磁控太阳捕获 方法。采用速率加太阳角控制加太阳角速度补偿控制算法,并引入了卫星速率 大小和太阳角大小改变磁控系数,从而改变卫星速率控制、太阳角控制和太阳 角速度控制在磁控计算中的比重,结合磁控+卫星反作用轮启旋偏置实现太阳 捕获。避免了先单纯进行速率阻尼控制、再单纯进行太阳角控制易出现角速度 变化较大引起太阳角变化起伏问题;同时使得在卫星速率较大或太阳角较大时 增大速率阻尼或太阳角速度速率阻尼的控制比重,在卫星速率较小且太阳角较 小时减小速率阻尼或太阳角速度速率阻尼的控制比重,而增大太阳角的控制比 重,提高太阳角的控制精度;使用反作用轮偏置避免了需要卫星启旋偏置的问题。
附图说明
为了进一步阐明本发明的各实施例的以上和其它优点和特征,将参考附图 来呈现本发明的各实施例的更具体的描述。可以理解,这些附图只描绘本发明 的典型实施例,因此将不被认为是对其范围的限制。在附图中,为了清楚明了, 相同或相应的部件将用相同或类似的标记表示。
图1示出根据本发明的一个实施例提供的一种考虑太阳角变化率和控制 因子权重可变的磁控太阳捕获方法的算法流程示意图。
图2示出根据本发明的具体实施例提供的一种考虑太阳角变化率和控制 因子权重可变的磁控太阳捕获方法仿真获得的帆板法线与太阳矢量夹角变化 示意图。
图3示出根据本发明的具体实施例提供的一种考虑太阳角变化率和控制 因子权重可变的磁控太阳捕获方法仿真获得的卫星三轴对地姿态角和角速度 变化示意图。
图4示出根据本发明的具体实施例提供的一种考虑太阳角变化率和控制 因子权重可变的磁控太阳捕获方法仿真获得的飞轮转速变化示意图。
具体实施方式
在以下的描述中,参考各实施例对本发明进行描述。然而,本领域的 技术人员将认识到可在没有一个或多个特定细节的情况下或者与其它替换 和/或附加方法、材料或组件一起实施各实施例。在其它情形中,未示出或 未详细描述公知的结构、材料或操作以免使本发明的各实施例的诸方面晦 涩。类似地,为了解释的目的,阐述了特定数量、材料和配置,以便提供 对本发明的实施例的全面理解。然而,本发明可在没有特定细节的情况下 实施。此外,应理解附图中示出的各实施例是说明性表示且不一定按比例 绘制。
在本说明书中,对“一个实施例”或“该实施例”的引用意味着结合 该实施例描述的特定特征、结构或特性被包括在本发明的至少一个实施例 中。在本说明书各处中出现的短语“在一个实施例中”并不一定全部指代 同一实施例。
需要说明的是,本发明的实施例以特定顺序对工艺步骤进行描述,然而这 只是为了方便区分各步骤,而并不是限定各步骤的先后顺序,在本发明的不同 实施例中,可根据工艺的调节来调整各步骤的先后顺序。
本发明提供一种考虑太阳角变化率和控制因子权重可变的磁控太阳捕获 方法。采用速率加太阳角控制加太阳角速度补偿控制算法,并引入了卫星速率 大小和太阳角大小改变磁控系数,从而改变卫星速率控制、太阳角控制和太阳 角速度控制在磁控计算中的比重,结合磁控+卫星反作用轮启旋偏置实现太阳 捕获。避免了先单纯进行速率阻尼控制、再单纯进行太阳角控制易出现角速度 变化较大引起太阳角变化起伏问题;同时使得在卫星速率较大或太阳角较大时 增大速率阻尼或太阳角速度速率阻尼的控制比重,在卫星速率较小且太阳角较 小时减小速率阻尼或太阳角速度速率阻尼的控制比重,而增大太阳角的控制比 重,提高太阳角的控制精度;使用反作用轮偏置避免了需要卫星启旋偏置的问题。
在介绍具体实施方式之前,先设定基本前提及参数:
1)卫星轨道为太阳同步轨道;
2)卫星的磁控周期为5秒;
3)卫星的姿控软件周期为0.5秒;
4)卫星磁力矩器控制方式为开关控制;
5)卫星对太阳定向控制方式为三轴磁控结合飞轮启旋偏置,并假设为Y 轮偏置控制,卫星飞轮轮组安装方式为三正交一斜装的安装方式。
下面结合图1来详细描述基于本发明的一个实施例提供的一种考虑太阳 角变化率和控制因子权重可变的磁控太阳捕获方法的具体实施方案。图1示出 根据本发明的一个实施例提供的一种考虑太阳角变化率和控制因子权重可变 的磁控太阳捕获方法的算法流程示意图。
首先,如图1所示,进行太阳角、太阳角速度以及磁矢量变化率的模的计 算。具体计算方式如下:
A)进行太阳角gama计算
太阳角即太阳矢量与帆板法线夹角,由于前述卫星的姿控软件周期为0.5 秒,设定每个姿控软件周期进行太阳角计算并转化为角度所示,计算公式如下:
其中Lb是卫星本体系中帆板产生太阳能面法线矢量在卫星本体系中的表 示,不同的安装方式,该矢量不同;Sb是卫星本体系中的太阳矢量。连续存储 最近十个软件周期内的太阳矢量与帆板法线夹角。gama为当前周期计算获得 的太阳角,gama-为前1软件周期太阳角,以此类推,gama10-为前10姿控软 件周期太阳角。
B)进行太阳角速度gamadot计算
在本发明的该实施例中,选取10个姿控软件周期10ΔT的太阳角变化来计 算太阳角速度,计算公式如下:
其中gamadot为太阳角速度,ΔT为姿控软件周期。
C)进行磁矢量变化率的模Bdotnorm计算
磁矢量变化率的模Bdotnorm计算进一步分为两步:
第一步,进行磁场矢量变化速度Bdot计算,计算公式如下:
单位nT/s,其中ΔT为姿控软件周期,在本发明的该实施例中,ΔT=0.5秒, 为上一次进行磁強计采样时的磁场矢量,这里假设磁控周期为5秒,磁強计 采样周期为5秒。
第二步,进行磁矢量变化率的模Bdotnorm计算,计算公式如下:
接下来,如图1所示,进行太阳矢量计算有效性判断。
当太阳矢量计算失效时,也即太敏无效及太阳矢量递推无效时,卫星仅进 行速率阻尼控制,磁矩计算公式如下:
Pb=-k1Bdot (5)
其中,Pb表示磁矩,k1表示角速度阻尼因子,Bdot为磁场矢量的变化率。
最后磁控输出磁矩限幅,即判断磁矩Pbi计算是否超过最大值Pbmax,若超 过最大值按最大值输出,即:如果|Pbi|≥Pbmax i=x,y,z,则Pbi=sign(Pbi)·Pbmax。
当太阳矢量计算有效时,如图1所示,基于飞轮启旋偏置完成情况进行施 加力矩方向Tn计算。
根据帆板的法线矢量和太阳矢量计算控制太阳角需要的作用力矩。该力矩 的计算需考虑卫星处于偏置情况和非偏置情况,也即:
若飞轮(轮Y)偏置完成,则计算公式如下:
Tn=(Lb×Sb)×Lb (6)
力矩方向Tn,Lb是帆板法线矢量,Sb是太阳矢量。
若飞轮(轮Y)偏置未完成,则计算公式如下:
Tn=Lb×Sb (7)
然后,对力矩方向Tn进行单位化,单位化方式如下:
接下来,再如图1所示,计算磁矩施加方向Ang_TnBb,并计算磁矩与力 矩方向的夹角aa。
由力矩方向计算磁矩施加方向的计算方式如下:
Ang_TnBb=Tn×Bb (9)
这里Ang_TnBb是非归一化的,Bb为磁矩方向。
力矩与磁矩夹角aa的计算公式如下:
然后,再如图1所示,判断磁矩与力矩的夹角aa是否处于上限aalimtup 和下限aalimtdown之间。
当磁矩与力矩的夹角aa超出上限或者低于下限时,仅进行速率阻尼控制, 磁矩计算公式如上面公式(5)所示,同样也进行磁控输出磁矩限幅控制。
当磁矩与力矩的夹角aa处于上限aalimtup和下限aalimtdown之间的时, 即当aa满足aalimtdown≤aa≤aalimtup时,采用速率加太阳角控制加太阳角速 度补偿控制的卫星帆板对日定向的太阳捕获控制算法和基于卫星速度和太阳 角而优化控制系数的控制算法,也即:
Pb1=-k2_Bdot (11)
Pb2=-(k6gma+k6k3_·gamadot)·Ang_TnBb (12)
Pb=Pb1+Pb2 (13)
这里k2_为磁场矢量的变化率Bdot的阻尼因子,k6为太阳角控制因子,k6k3为太阳角速度控制因子。
采用速率加太阳角控制加太阳角速度补偿控制算法是本发明技术方案的 一个创新点,而本发明技术方案的另一个创新点是引入了卫星磁矢量变化率和 太阳角大小改变磁控因子,从而改变卫星速率控制、太阳角控制和太阳角速度 控制因子在磁控计算中的权重。具体实现如下:
如果Bdotnorm>BdotLimt(120nT)或gama>GamaLimt(40deg),则满足此 条件时卫星角速度较大或者太阳角大于阈值GamaLimt,此时速率阻尼或太阳 角速度控制因子的权重较大,以达到快速减小卫星角速度和太阳角速度的作用, 通过修正控制系数来调整k2和k3的值,计算公式如下所示。
k2_=k2·k4 (14)
k3_=k3·k5
这里k4是速率阻尼修正控制系数,k5是太阳角速度修正控制系数,其值为 均大于1的正实数。
如果Bdotnorm≤BdotLimt(120nT)且gama≤GamaLimt(40deg)条件满足, 也即此时Bdotnorm较小,卫星角速度较小且帆板与太阳矢量的夹角小于等于 GamaLimt时,提高太阳角控制因子权重,保障太阳角帆板对太阳定向的控制 精度。
k2_=k2 (15)
k3_=k3
最后,再如图1所示,进行磁控输出磁矩限幅,判断磁矩计算是否超过最 大值,若超过最大值按最大值输出,即:如果|Pbi|≥Pbmax i=x,y,z,则 Pbi=sign(Pbi)·Pbmax。
同时,在每次进行算法执行的同时,进行飞轮启旋偏置,具体的飞轮启旋 偏置方法如下:
判断Bdotnorm<200nT/s,且当前太阳角小于等于50度时,Y轮飞轮按每 0.5秒变化-0.2rpm启旋偏置到进入安全模式的初始转速-300rpm。卫星偏置的 飞轮转速和偏置转速变化率应根据卫星实际的惯量、飞轮的惯量以及帆板的安 装方向进行合理的设置。
启旋偏置飞轮的选择与帆板进行太阳能电池面的法线有关,这里假设帆板 太阳能电池面的法线矢量在卫星本体系中表示为[0-10],则选择卫星的Y 轴飞轮进行启旋。
下面结合图2-图4来描述基于本发明的考虑太阳角变化率和控制因子权 重可变的磁控太阳捕获方法进行的具体仿真实例。
首先设定仿真的初始条件,具体设置如下:
轨道设置:太阳同步轨道,轨道高度600km,降交点地方时6:00;
初始对地角度:[170,40,20]度;
对地角速度:[1.5 1.5 1.5]度/秒;
初始太阳矢量与帆板夹角:167.14度;
帆板太阳能电池面法线矢量:[0 -1 0];
卫星惯量设置:[500 600 500 0 0 0]kgm2;
最终偏置状态:;卫星Y轴3Nms偏置;
偏置方式:Y飞轮偏置-300rpm;
飞轮初始转速:[0 0 0 0]rpm;
偏置目标转速:[0-300 0 0]rpm。
仿真结果如图2-图4所示,图2示出根据本发明的具体实施例提供的一种 考虑太阳角变化率和控制因子权重可变的磁控太阳捕获方法仿真获得的帆板 法线与太阳矢量夹角变化示意图;图3示出根据本发明的具体实施例提供的一 种考虑太阳角变化率和控制因子权重可变的磁控太阳捕获方法仿真获得的卫 星三轴对地姿态角和角速度变化示意图;图4示出根据本发明的具体实施例提 供的一种考虑太阳角变化率和控制因子权重可变的磁控太阳捕获方法仿真获 得的飞轮转速变化示意图。
仿真结果分析:可以看出本发明的技术方案可以将卫星角速度消旋并完成 大角速度大角度太阳捕获,并最终将太阳帆板法线与太阳矢量夹角稳定控制在 10度以内。
本发明提供的该种考虑太阳角变化率和控制因子权重可变的磁控太阳捕 获方法。采用速率加太阳角控制加太阳角速度补偿控制算法,并引入了卫星速 率大小和太阳角大小改变磁控系数,从而改变卫星速率控制、太阳角控制和太 阳角速度控制在磁控计算中的比重,结合磁控+卫星反作用轮启旋偏置实现太 阳捕获。避免了先单纯进行速率阻尼控制、再单纯进行太阳角控制易出现角速 度变化较大引起太阳角变化起伏问题;同时使得在卫星速率较大或太阳角较大 时增大速率阻尼或太阳角速度速率阻尼的控制比重,在卫星速率较小且太阳角 较小时减小速率阻尼或太阳角速度速率阻尼的控制比重,而增大太阳角的控制 比重,提高太阳角的控制精度;使用反作用轮偏置避免了需要卫星启旋偏置的问题。
尽管上文描述了本发明的各实施例,但是,应该理解,它们只是作为 示例来呈现的,而不作为限制。对于相关领域的技术人员显而易见的是, 可以对其做出各种组合、变型和改变而不背离本发明的精神和范围。因此, 此处所公开的本发明的宽度和范围不应被上述所公开的示例性实施例所限 制,而应当仅根据所附权利要求书及其等同替换来定义。
Claims (13)
1.一种考虑太阳角变化率和控制因子权重可变的磁控太阳捕获方法,包括:
进行太阳角、太阳角速度以及磁矢量变化率的模的计算;
进行太阳矢量计算有效性判断;
当太阳矢量计算有效时,基于飞轮启旋偏置完成情况进行施加力矩方向计算;
进行磁矩施加方向计算,并计算磁矩与力矩方向的夹角;
进行磁矩与力矩方向的夹角是否在上下限之间的范围内判断;
当磁矩与力矩的夹角处于上下限之间的范围内,采用速率加太阳角控制加太阳角速度补偿算法计算磁控输出;以及
进行磁控输出磁矩限幅。
2.如权利要求1所述的方法,其特征在于,所述的太阳角、太阳角速度以及磁矢量变化率的模的计算的公式分别为:
其中gama为当前周期计算获得的太阳角,Lb是卫星本体系中帆板产生太阳能面法线矢量,Sb是卫星本体系中的太阳矢量;
其中gamadot为太阳角速度,ΔT为姿控软件周期;
磁矢量变化率的模Bdotnorm计算分为两步,先进行磁场矢量变化速度Bdot计算:
其中10ΔT为磁強计采样周期,Bb为本次进行磁強计采样时的磁场矢量,Bb -为上一次进行磁強计采样时的磁场矢量,
再进行磁矢量变化率的模Bdotnorm计算:
3.如权利要求1所述的方法,其特征在于,所述的进行太阳矢量计算有效性判断即判断太阳矢量是否太敏无效或递推无效。
4.如权利要求1所述的方法,其特征在于,所述的基于飞轮启旋偏置完成情况进行施加力矩方向计算方法为:
若飞轮偏置完成,则计算方法为Tn=(Lb×Sb)×Lb;
若飞轮(轮Y)偏置未完成,则计算方法为Tn=Lb×Sb;以及
对力矩方向Tn进行单位化,单位化方法为其中Tn为力矩方向,Lb为帆板法线矢量,Sb是太阳矢量。
5.如权利要求1所述的方法,其特征在于,所述的进行磁矩施加方向计算,并计算磁矩与力矩方向的夹角的方法为:
磁矩施加方向由力矩方向计算,计算方法为Ang_TnBb=Tn×Bb,磁矩与力矩方向的夹角的计算方法为其中Ang_TnBb是非归一化的,Bb为磁矢量,Tn为力矩方向。
6.如权利要求1所述的方法,其特征在于,所述的速率加太阳角控制加太阳角速度补偿算法的计算方法为:
Pb=-k2_Bdot-(k6gama+k6k3_.gamadot).Ang_TnBb
其中Pb表示磁矩,Bdot为磁场矢量的变化率,k2_为Bdot阻尼因子,gama为太阳角,k6为太阳角控制因子,gamadot为太阳角速度,k6k3为太阳角速度控制因子,Ang_TnBb为磁矩施加方向。
7.如权利要求1所述的方法,其特征在于,在进行所述的采用速率加太阳角控制加太阳角速度补偿算法计算磁控输出前,还包括判断卫星角速度是否大于第一阈值或者太阳角是否大于第二阈值。
8.如权利要求7所述的方法,其特征在于,如果卫星角速度大于第一阈值或者太阳角大于第二阈值,则采用高权重速率阻尼控制系数;如果卫星角速度小于第一阈值且太阳角小于第二阈值,则采用高权重太阳角控制系数。
9.如权利要求7所述的方法,其特征在于,所述的高权重速率阻尼控制系数为k2_=k2.k4,k3_=k3.k5,其中k4是速率阻尼修正控制系数,k5是太阳角速度修正控制系数,其值为均大于1的正实数;所述的高权重太阳角控制系数为k2_=k2,k3_=k3。
10.如权利要求1所述的方法,其特征在于,当太阳矢量计算失效时,卫星仅进行速率阻尼控制,并进行速率阻尼控制磁控输出的磁矩限幅。
11.如权利要求1所述的方法,其特征在于,当磁矩与力矩的夹角超出上下限之间的范围时,卫星仅进行速率阻尼控制,并进行速率阻尼控制磁控输出的磁矩限幅。
12.如权利要求10或11所述的方法,其特征在于,所述速率阻尼控制的技术公式为Pb=-k1Bdot,其中Pb表示磁矩,k1表示角速度阻尼因子,Bdot为磁场矢量的变化率。
13.如权利要求1所述的方法,其特征在于,还包括在每次进行算法执行的同时,进行飞轮启旋条件判断,若磁矢量变化率的模小于第三阈值,且当前太阳角小于等于第四阈值时,飞轮启旋;否则飞轮转速保持。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810304610.9A CN108549412B (zh) | 2018-04-08 | 2018-04-08 | 一种考虑太阳角变化率和控制因子的磁控太阳捕获方法 |
CN202011395040.2A CN112389681B (zh) | 2018-04-08 | 2018-04-08 | 一种用于卫星对日定向的磁控方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810304610.9A CN108549412B (zh) | 2018-04-08 | 2018-04-08 | 一种考虑太阳角变化率和控制因子的磁控太阳捕获方法 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011395040.2A Division CN112389681B (zh) | 2018-04-08 | 2018-04-08 | 一种用于卫星对日定向的磁控方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108549412A true CN108549412A (zh) | 2018-09-18 |
CN108549412B CN108549412B (zh) | 2020-11-24 |
Family
ID=63514002
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011395040.2A Active CN112389681B (zh) | 2018-04-08 | 2018-04-08 | 一种用于卫星对日定向的磁控方法 |
CN201810304610.9A Active CN108549412B (zh) | 2018-04-08 | 2018-04-08 | 一种考虑太阳角变化率和控制因子的磁控太阳捕获方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011395040.2A Active CN112389681B (zh) | 2018-04-08 | 2018-04-08 | 一种用于卫星对日定向的磁控方法 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN112389681B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109649693A (zh) * | 2019-01-21 | 2019-04-19 | 上海微小卫星工程中心 | 一种纯磁控自旋对日定向方法 |
CN111874269A (zh) * | 2020-08-10 | 2020-11-03 | 吉林大学 | 一种磁控小卫星的低功耗对日捕获及定向姿态控制方法 |
CN113353292A (zh) * | 2021-06-26 | 2021-09-07 | 山东航天电子技术研究所 | 磁控非自旋对日定向方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0372434B1 (en) * | 1988-12-09 | 1993-09-01 | Hughes Aircraft Company | Satellite control system |
US6285927B1 (en) * | 1999-05-26 | 2001-09-04 | Hughes Electronics Corporation | Spacecraft attitude determination system and method |
CN101381004A (zh) * | 2008-08-20 | 2009-03-11 | 南京航空航天大学 | 基于大气阻力的微小卫星编队飞行控制方法及控制装置 |
US20100001142A1 (en) * | 2004-12-07 | 2010-01-07 | Lockheed Martin Corporation | Optimized land mobile satellite configuration and steering |
CN101934863A (zh) * | 2010-09-29 | 2011-01-05 | 哈尔滨工业大学 | 基于磁力矩器和飞轮的卫星姿态全方位控制方法 |
CN102591349A (zh) * | 2012-03-12 | 2012-07-18 | 北京控制工程研究所 | 高轨道卫星大初始角速率情况的无陀螺太阳捕获控制方法 |
CN102582850A (zh) * | 2012-03-16 | 2012-07-18 | 上海微小卫星工程中心 | 提高卫星磁控精度的方法 |
CN104097791A (zh) * | 2014-06-24 | 2014-10-15 | 上海微小卫星工程中心 | 一种基于磁敏感器和星敏感器的全姿态捕获方法及其装置 |
CN104097793A (zh) * | 2014-06-24 | 2014-10-15 | 上海微小卫星工程中心 | 一种卫星零动量磁控太阳捕获装置及方法 |
CN105966639A (zh) * | 2016-05-11 | 2016-09-28 | 上海微小卫星工程中心 | 一种卫星对日自旋稳定控制系统及方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB593446A (en) * | 1943-12-31 | 1947-10-16 | Bendix Aviat Corp | Step-by-step follow-up and anti-hunt mechanism |
US6231011B1 (en) * | 1998-11-02 | 2001-05-15 | University Of Houston System | Satellite angular momentum control system using magnet-superconductor flywheels |
CN103274060B (zh) * | 2013-04-27 | 2015-04-22 | 中国空间技术研究院 | 一种基于日光反射的航天器补能系统 |
CN103365302B (zh) * | 2013-06-18 | 2016-03-09 | 陕西理工学院 | 三相限太阳光线传感器与光斑图像传感器融合的太阳光线跟踪算法 |
US9981446B2 (en) * | 2013-09-03 | 2018-05-29 | The Boeing Company | Structural inserts for honeycomb structures |
DE102015114819B3 (de) * | 2015-09-04 | 2016-12-22 | Rockwell Collins Deutschland Gmbh | Drallradvorrichtung zur Lagestabilisierung eines Raumflugkörpers |
US9663252B1 (en) * | 2016-12-07 | 2017-05-30 | Beihang University | Method for attitude controlling based on finite time friction estimation for flexible spacecraft |
CN106809406B (zh) * | 2017-01-19 | 2019-04-30 | 上海航天控制技术研究所 | 一种基于地磁力矩的飞轮起旋控制方法 |
CN107600464B (zh) * | 2017-09-18 | 2019-08-23 | 上海航天控制技术研究所 | 利用星敏感器信息的飞轮控制捕获太阳及对日定向方法 |
-
2018
- 2018-04-08 CN CN202011395040.2A patent/CN112389681B/zh active Active
- 2018-04-08 CN CN201810304610.9A patent/CN108549412B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0372434B1 (en) * | 1988-12-09 | 1993-09-01 | Hughes Aircraft Company | Satellite control system |
US6285927B1 (en) * | 1999-05-26 | 2001-09-04 | Hughes Electronics Corporation | Spacecraft attitude determination system and method |
US20100001142A1 (en) * | 2004-12-07 | 2010-01-07 | Lockheed Martin Corporation | Optimized land mobile satellite configuration and steering |
CN101381004A (zh) * | 2008-08-20 | 2009-03-11 | 南京航空航天大学 | 基于大气阻力的微小卫星编队飞行控制方法及控制装置 |
CN101934863A (zh) * | 2010-09-29 | 2011-01-05 | 哈尔滨工业大学 | 基于磁力矩器和飞轮的卫星姿态全方位控制方法 |
CN102591349A (zh) * | 2012-03-12 | 2012-07-18 | 北京控制工程研究所 | 高轨道卫星大初始角速率情况的无陀螺太阳捕获控制方法 |
CN102582850A (zh) * | 2012-03-16 | 2012-07-18 | 上海微小卫星工程中心 | 提高卫星磁控精度的方法 |
CN104097791A (zh) * | 2014-06-24 | 2014-10-15 | 上海微小卫星工程中心 | 一种基于磁敏感器和星敏感器的全姿态捕获方法及其装置 |
CN104097793A (zh) * | 2014-06-24 | 2014-10-15 | 上海微小卫星工程中心 | 一种卫星零动量磁控太阳捕获装置及方法 |
CN105966639A (zh) * | 2016-05-11 | 2016-09-28 | 上海微小卫星工程中心 | 一种卫星对日自旋稳定控制系统及方法 |
Non-Patent Citations (1)
Title |
---|
侯建文: "FY-1C卫星姿态控制系统", 《上海航天》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109649693A (zh) * | 2019-01-21 | 2019-04-19 | 上海微小卫星工程中心 | 一种纯磁控自旋对日定向方法 |
CN113955152A (zh) * | 2019-01-21 | 2022-01-21 | 上海微小卫星工程中心 | 一种星体对日定向控制方法 |
CN113955152B (zh) * | 2019-01-21 | 2024-03-01 | 上海微小卫星工程中心 | 一种星体对日定向控制方法 |
CN111874269A (zh) * | 2020-08-10 | 2020-11-03 | 吉林大学 | 一种磁控小卫星的低功耗对日捕获及定向姿态控制方法 |
CN111874269B (zh) * | 2020-08-10 | 2022-02-01 | 吉林大学 | 一种磁控小卫星的低功耗对日捕获及定向姿态控制方法 |
CN113353292A (zh) * | 2021-06-26 | 2021-09-07 | 山东航天电子技术研究所 | 磁控非自旋对日定向方法 |
CN113353292B (zh) * | 2021-06-26 | 2022-06-07 | 山东航天电子技术研究所 | 磁控非自旋对日定向方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112389681B (zh) | 2022-05-17 |
CN108549412B (zh) | 2020-11-24 |
CN112389681A (zh) | 2021-02-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108549412A (zh) | 一种考虑太阳角变化率和控制因子权重可变的磁控太阳捕获方法 | |
Carrasco et al. | Aperture-coupled reflectarray element with wide range of phase delay | |
CN106096148B (zh) | 一种简单姿态控制下的大倾角轨道卫星太阳帆板指向方法 | |
CN104527994A (zh) | 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法 | |
CN103072701B (zh) | 一种欠驱动的卫星消旋控制方法 | |
CN102923317B (zh) | 一种适用于卫星姿态角速度阻尼的欠驱动控制方法 | |
CN103034237B (zh) | 使用两个单框架控制力矩陀螺的航天器姿态机动控制方法 | |
CN108583938B (zh) | 一种可应用于运行于太阳同步晨昏轨道的全向天线通信卫星姿态控制系统及其方法 | |
CN108255062B (zh) | 基于改进差分进化机理的动力定位节能推力分配方法 | |
CN109625329A (zh) | 一种基于离散喷气的飞轮角动量自主卸载方法 | |
CN109649693B (zh) | 一种纯磁控自旋对日定向方法 | |
Zhang et al. | Hybrid threshold event-triggered control for sail-assisted USV via the nonlinear modified LVS guidance | |
CN106292297A (zh) | 基于pid控制器和l1自适应控制器的姿态控制方法 | |
CN109977587A (zh) | 一种直流电机推进无人水面艇速度控制器设计方法 | |
CN109445448A (zh) | 一种轮控小卫星自适应积分滑模姿态控制器 | |
CN109144085A (zh) | 基于特征结构配置的鲁棒Hinf航天器指向控制方法及系统 | |
CN110498063A (zh) | 一种利用太阳敏感器的全姿态序列对日定向方法 | |
US6439511B1 (en) | Thermal shock avoidance for satellite solar panels | |
Wang et al. | Solar sail spacecraft trajectory optimization based on improved imperialist competitive algorithm | |
CN110712769A (zh) | 一种基于cmg的无陀螺太阳定向控制方法 | |
CN114661056B (zh) | 一种考虑推进器伺服控制的差动水面无人船轨迹跟踪方法 | |
Zhang et al. | An adjustable control for inertia momentum wheel with disturbance compensation | |
Dong et al. | Attitude control design for aerospace vehicle with blended reaction control systems and aerodynamic surfaces | |
CN107757951A (zh) | 一种微小卫星零动量实现太阳能捕获的装置及方法 | |
Zheng et al. | RBF adaptive sliding control for five-axis flexible satellite |
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 |