一种利用磁力矩器卸载卫星合成角动量的双滞环方法技术

本发明专利技术涉及一种利用磁力矩器卸载卫星合成角动量的双滞环方法，属于卫星姿态控制技术领域。该方法针对磁卸载在地磁场方向刚满足卸载条件的轨道弧段附近磁控电压的频繁切换的问题，针对角动量误差以及地磁场矢量与角动量误差矢量夹角引入双滞环卸载策略，最后对卸载电压进行滤波平滑后输出，有效降低磁控力矩对卫星姿态稳定度的影响，实现了卫星的高稳定度姿态控制。相比传统磁卸载方法，采用本发明专利技术的方法可降低磁控力矩引起的卫星姿态角和角速度的波动，有效克服磁卸载对卫星的姿态稳定度产生的恶劣影响。

A Double Hysteresis Method for Synthesizing Angular Momentum of Unloading Satellite Using Magneto-Torque

【技术实现步骤摘要】
一种利用磁力矩器卸载卫星合成角动量的双滞环方法
本专利技术涉及一种利用磁力矩器卸载卫星合成角动量的双滞环方法，属于卫星姿态控制

技术介绍
卫星在轨飞行过程中，会受到各种环境干扰的影响，使卫星姿态角和姿态角速度出现偏差。为了保持卫星三轴姿态的稳定，当代卫星一般采用动量轮等角动量交换装置，将外界环境对卫星姿态的干扰转化为角动量存储在动量轮中，当角动量达到一定程度时，将启动磁力矩器产生磁卸载力矩对动量轮进行角动量卸载，这就是所谓的角动量的“磁卸载”。传统磁卸载一般根据卫星的角动量积累情况进行磁控电压计算，但是在刚满足卸载条件时，会出现磁控电压频繁切换的情况，此时磁控力矩会引起卫星姿态角和角速度的波动，对卫星的高稳定度控制产生恶劣影响。
技术实现思路
本专利技术的技术解决问题是：克服现有技术的不足，提出一种利用磁力矩器卸载卫星合成角动量的双滞环方法，该方法首先在卫星角动量误差超过Δhlimit1且地磁场矢量与角动量误差矢量夹角大于ξlimit1度时进行磁卸载电压计算，直到角动量误差小于Δhlimit0或地磁场矢量与角动量误差矢量夹角小于ξlimit0度时为止。计算出的卸载电压经过滤波平滑后输出，进一步降低磁控力矩对卫星姿态稳定度的影响，实现了卫星的高稳定度姿态控制。本专利技术的技术解决方案是：一种利用磁力矩器卸载卫星合成角动量的双滞环方法，使用三台沿本体坐标系正装的磁力矩器卸载卫星在本体坐标系下三个方向的合成角动量，该方法包括如下步骤：(1)根据卫星在轨道上的位置，计算卫星所处空间环境的磁场Bb；Bb＝[BbxBbyBbz]T，其中Bbx，Bby，Bbz为磁场强度在卫星本体坐标系的三轴分量；(2)计算卫星的合成角动量ΔHb；初始t0时刻卫星的合成角动量ΔHb0＝0；t1时刻卫星的合成角动量t2时刻卫星的合成角动量…ti时刻卫星的合成角动量ti+1时刻卫星的合成角动量其中，Js为卫星的转动惯量，ωOI为卫星轨道坐标系相对于惯性坐标系的角速度；ωt1为t1时刻卫星相对于惯性坐标系的角速度，CBO1为t1时刻卫星从轨道坐标系到本体坐标系的转换矩阵，为t1时刻卫星动量轮组合的合成角动量；ωt2为t2时刻卫星相对于惯性坐标系的角速度，CBO2为t2时刻卫星从轨道坐标系到本体坐标系的转换矩阵，为t2时刻卫星动量轮组合的合成角动量；ωti为ti时刻卫星相对于惯性坐标系的角速度，CBOi为ti时刻卫星从轨道坐标系到本体坐标系的转换矩阵，为ti时刻卫星动量轮组合的合成角动量；ωti+1为ti+1时刻卫星相对于惯性坐标系的角速度，CBOi+1为ti+1时刻卫星从轨道坐标系到本体坐标系的转换矩阵，为ti+1时刻卫星动量轮组合的合成角动量；(3)计算当前时刻步骤(1)得到的卫星所处空间环境的磁场Bb的方向与步骤(2)得到的当前时刻卫星的合成角动量ΔHb的夹角ε的余弦值cosε的绝对值|cosε|；(4)计算当前时刻卫星的合成角动量ΔHb的范数||ΔHb||；(5)根据步骤(3)得到的余弦值cosε的绝对值|cosε|和步骤(4)得到的范数||ΔHb||判断当前时刻卫星的合成角动量是否需要进行卸载，初始时刻时不对卫星的合成角动量进行卸载；第一种情况，当当前时刻|cosε|<a时且当前时刻||ΔHb||＞b时对卫星的合成角动量进行卸载；第二种情况，当当前时刻|cosε|＞c时或当前时刻||ΔHb||<d时不对卫星的合成角动量卸载；a、b、c、d均为设定阈值，且a<c，b>d；第三种情况，当当前时刻|cosε|<a时且当前时刻d≤||ΔHb||≤b时，根据当前时刻的前一时刻的卸载情况进行卸载或不卸载；第四种情况，当当前时刻||ΔHb||＞b时且当前时刻a≤|cosε|≤c时，根据当前时刻的前一时刻的卸载情况进行卸载或不卸载；第五种情况，当当前时刻a≤|cosε|≤c时且当前时刻d≤||ΔHb||≤b时，根据当前时刻的前一时刻的卸载情况进行卸载或不卸载；(6)根据步骤(5)的判断结果，如果需要进行卸载则进入步骤(6)，如果不需要进行卸载则输出给磁力矩器的卸载电压为零；当判断结果为进行卸载时，则卫星本体坐标系下x方向的卸载电压Vm1*为：卫星本体坐标系下y方向的卸载电压Vm2*为：卫星本体坐标系下z方向的卸载电压Vm3*为：其中，kPMU为磁卸载系数，为一设定值；Mmlf为磁力矩器的最大磁矩；Vmmlf为电压的限幅值，为一设定值；ΔHbx、ΔHby、ΔHbz为卫星的合成角动量ΔHb在卫星本体坐标系下的三个分量；Bbx、Bby、Bbz为卫星所处空间环境的磁场Bb在卫星本体坐标系下的三个分量；(7)对步骤(6)得到的卸载电压进行滤波，输出电压给磁力矩器，使磁力矩器对卫星在本体坐标系下三个方向的合成角动量进行卸载；滤波公式如下：其中kvm为滤波系数，为一设定值，为上一时刻输出的滤波后的电压。本专利技术与现有技术相比的有益效果是：(1)本专利技术利用合成角动量及其与磁场强度夹角的双滞环设计，完全解决了磁控电压频繁切换引起的航天器姿态角速度波动的问题；(2)本专利技术对磁控电压进行了滤波，使用该算法能够有效降低磁卸载对航天器姿态稳定度的影响；(3)本专利技术提出的算法所要求的计算量小，无需增加额外的计算资源，适合星载计算机的实现。(4)本专利技术涉及一种利用磁力矩器卸载的双滞环的方法，属于卫星姿态机动控制
该方法针对磁卸载在地磁场方向刚满足卸载条件的轨道弧段附近磁控电压的频繁切换的问题，针对角动量误差以及地磁场矢量与角动量误差矢量夹角引入双滞环卸载策略，最后对卸载电压进行滤波平滑后输出，有效降低磁控力矩对卫星姿态稳定度的影响，实现了卫星的高稳定度姿态控制。相比传统磁卸载方法，采用本专利技术的方法可降低磁控力矩引起的卫星姿态角和角速度的波动，有效克服磁卸载对卫星的姿态稳定度产生的恶劣影响。附图说明图1为夹角卸载条件判断滞环示意图；图2为角动量卸载条件判断滞环示意图；图3为本体坐标系下的磁场强度示意图；图4为本体坐标系下的卫星合成角动量示意图；图5为磁力矩器输出电压示意图。具体实施方式以下结合附图和具体实施例对本专利技术进行详细说明。实施例下面对本专利技术进行具体说明：假设为运行在轨道高度约500km，倾角97°的太阳同步轨道上。卫星配置了三台最大磁矩为Mmlf＝200Am2的磁力矩器，如图1-5所示；(1)根据卫星的地理位置，计算卫星所处空间环境的磁场Bb；假设卫星当前所在地理位置为东经94°，北纬0.03°，高度492km，其环境磁场强度在本体坐标系下的三轴分量为t1时刻：Bb＝[-3.0084979825743828e-005,5.9109881461567932e-006,-8.4364482182156017e-006]nT。t2时刻：Bb＝[-6.8490665271550280e-007,1.0552072686590731e-005,-4.2079269982111809e-005]nT。(2)计算卫星系统的合成角动量ΔHb；假设卫星转动惯量为JS＝[14621.6,0.61,-267.410.61,12589.8,-365.91-267.41,-365.91,9798.9]km2，t1时刻：惯性姿态角速度ωt＝[-2.214本文档来自技高网...

【技术保护点】
1.一种利用磁力矩器卸载卫星合成角动量的双滞环方法，其特征在于：使用三台沿本体坐标系正装的磁力矩器卸载卫星在本体坐标系下三个方向的合成角动量，该方法包括如下步骤：(1)根据卫星在轨道上的位置，计算卫星所处空间环境的磁场Bb；(2)计算卫星的合成角动量ΔHb；初始t0时刻卫星的合成角动量ΔHb0＝0；t1时刻卫星的合成角动量ΔHb1＝Js·(ωt1–CBO1·ωOI)+ΔHΣmw1；t2时刻卫星的合成角动量ΔHb2＝Js·(ωt2–CBO2·ωOI)+ΔHΣmw2；…ti时刻卫星的合成角动量ΔHbi＝Js·(ωti–CBOi·ωOI)+ΔHΣmwi；ti+1时刻卫星的合成角动量ΔHbi+1＝Js·(ωti+1–CBoi+1·ωOI)+ΔHΣmwi+1；(3)计算当前时刻步骤(1)得到的卫星所处空间环境的磁场Bb的方向与步骤(2)得到的当前时刻卫星的合成角动量ΔHb的夹角ε的余弦值cosε的绝对值|cosε|；(4)计算当前时刻卫星的合成角动量ΔHb的范数||ΔHb||；(5)根据步骤(3)得到的余弦值cosε的绝对值|cosε|和步骤(4)得到的范数||ΔHb||判断当前时刻卫星的合成角动量是否需要进行卸载，如果需要进行卸载则进入步骤(6)，如果不需要进行卸载则输出给磁力矩器的卸载电压为零；(6)根据步骤(5)的判断结果，当判断结果为进行卸载时，则卫星本体坐标系下x方向的卸载电压Vm1...

【技术特征摘要】
1.一种利用磁力矩器卸载卫星合成角动量的双滞环方法，其特征在于：使用三台沿本体坐标系正装的磁力矩器卸载卫星在本体坐标系下三个方向的合成角动量，该方法包括如下步骤：(1)根据卫星在轨道上的位置，计算卫星所处空间环境的磁场Bb；(2)计算卫星的合成角动量ΔHb；初始t0时刻卫星的合成角动量ΔHb0＝0；t1时刻卫星的合成角动量ΔHb1＝Js·(ωt1–CBO1·ωOI)+ΔHΣmw1；t2时刻卫星的合成角动量ΔHb2＝Js·(ωt2–CBO2·ωOI)+ΔHΣmw2；…ti时刻卫星的合成角动量ΔHbi＝Js·(ωti–CBOi·ωOI)+ΔHΣmwi；ti+1时刻卫星的合成角动量ΔHbi+1＝Js·(ωti+1–CBoi+1·ωOI)+ΔHΣmwi+1；(3)计算当前时刻步骤(1)得到的卫星所处空间环境的磁场Bb的方向与步骤(2)得到的当前时刻卫星的合成角动量ΔHb的夹角ε的余弦值cosε的绝对值|cosε|；(4)计算当前时刻卫星的合成角动量ΔHb的范数||ΔHb||；(5)根据步骤(3)得到的余弦值cosε的绝对值|cosε|和步骤(4)得到的范数||ΔHb||判断当前时刻卫星的合成角动量是否需要进行卸载，如果需要进行卸载则进入步骤(6)，如果不需要进行卸载则输出给磁力矩器的卸载电压为零；(6)根据步骤(5)的判断结果，当判断结果为进行卸载时，则卫星本体坐标系下x方向的卸载电压Vm1*为：卫星本体坐标系下y方向的卸载电压Vm2*为：卫星本体坐标系下z方向的卸载电压Vm3*为：其中，kPMU为磁卸载系数，为一设定值；Mmlf为磁力矩器的最大磁矩；Vmmlf为电压的限幅值，为一设定值；ΔHbx、ΔHby、ΔHbz为卫星的合成角动量ΔHb在卫星本体坐标系下的三个分量；Bbx、Bby、Bbz为卫星所处空间环境的磁场Bb在卫星本体坐标系下的三个分量；(7)对步骤(6)得到的卸载电压进行滤波，输出电压给磁力矩器，使磁力矩器对卫星在本体坐标系下三个方向的合成角动量进行卸载；滤波公式如下：其中kvm为滤波系数，为上一时刻输出的滤波后的电压。2.根据权利要求1所述的一种利用磁力矩器卸载卫星合成角动量的双滞环方法，其特征在于：所述的步骤(1)中，Bb＝[BbxBbyBbz]T，其中Bbx，Bby，Bbz为磁场强度在卫星本体坐标系的三轴分量。3.根据权利要求1所述的一种利用磁力矩器卸载卫星合成角动量的双滞环方法，其特征在于：所述的步骤(2)中，Js为卫星的转动惯量，ωOI为卫星轨道坐标系相...

【专利技术属性】
技术研发人员：陆栋宁，王淑一，雷拥军，顾斌，刘洁，田科丰，傅秀涛，
申请(专利权)人：北京控制工程研究所，
类型：发明
国别省市：北京,11