卫星动量轮在轨姿态无偏差起旋和消旋的控制方法

摘要

本发明提供了卫星动量轮在轨姿态无偏差起旋和消旋的控制方法，包括如下步骤：遥控注数设置卫星最终需要的目标角动量；管理反作用飞轮的动量，控制所述卫星在期望中心转速范围内工作；计算磁前馈力矩将卫星控制器计算的输出力矩磁前馈力矩常规的解耦力矩和干扰补偿力矩合成得到三轴指令力矩控制周期内积分得到三轴卫星的指令角动量，根据卫星的安装矩阵计算得到正在使用的卫星的指令转速。通过磁控力矩作用在星体上，同时通过前馈力矩的方式合成在飞轮的指令转速中抵消磁力矩对星体的影响，由此达到偏置角动量增加和减少过程中姿态无偏差起旋和消旋的的效果。

权利要求

1.卫星动量轮在轨姿态无偏差起旋和消旋的控制方法，其特征在于，包括如下步骤：A、遥控注数设置卫星最终需要的目标角动量；

B、管理飞轮的动量，控制所述卫星在期望中心转速范围内工作，当设置了起旋和消旋指令和俯仰轴目标转速后，增加至磁卸载目标角动量；

C、根据磁卸载力矩 和俯仰轴干扰补偿力矩，计算磁前馈力矩D、将卫星控制器计算的输出力矩 磁前馈力矩 常规的解耦力矩 和干扰补偿力矩 合成得到三轴指令力矩 控制周期内积分得到三轴卫星的指令角动量，根据卫星的安装矩阵计算得到正在使用的卫星的指令转速；

步骤C中所述磁前馈力矩 的计算公式为：其中，Tdyconst为俯仰轴常值干扰力矩，根据卫星在轨实际需要进行补偿；

所述磁卸载力矩 的计算公式为：其中， 为磁场强度在卫星本体的分量，kx,ky,kz为卸载系数，为三轴卸载的角动量，正常情况下 当软件判断处于偏置动量轮起旋模式时，自主将卸载角动量调整为

其中，hQX＝JΩ0为卫星俯仰轴希望达到的目标角动量；

三轴指令角动量

其中， 为切换到起/消旋控制时三轴角动量；

将需要增加或消除的偏置角动量合成在所述磁卸载力矩 中，通过磁控力矩作用在星体上。

2.根据权利要求1所述的控制方法，其特征在于：步骤A中所述目标角动量包括起旋时的俯仰轴偏置转速和/或角动量。

3.根据权利要求1所述的控制方法，其特征在于：步骤A中所述目标角动量包括消旋时的0转速和/或零动量。

4.根据权利要求1所述的控制方法，其特征在于：步骤A中所述卫星目标角动量为一次性设置。

5.根据权利要求1所述的控制方法，其特征在于：步骤D中所述三轴指令力矩 的计算公式为：

6.根据权利要求5所述的控制方法，其特征在于：所述输出力矩 的计算公式为：其中，KPi、KIi和KDi(i＝x,y,z)为控制器参数， θc、ψc为控制用姿态角，ωox、ωoy、ωoz为根据惯性系角速度测量值、轨道及姿态角解算的星体轨道系角速度。

说明书

卫星动量轮在轨姿态无偏差起旋和消旋的控制方法

技术领域

[0001] 本发明属于卫星姿态控制技术领域，特别涉及卫星动量轮在轨姿态无偏差起旋和消旋的控制方法。

背景技术

[0002] 随着民用航天技术与应用需求的发展，国家对于具备对地观测功能的民用卫星寿命的要求越来越高，中低轨的对地观测卫星也逐步向八年以上在轨寿命期发展。

[0003] 为确保完成卫星在轨应用任务，需要对其运动进行精确控制，卫星姿态控制系统是卫星重要的组成部分，随着现代卫星技术的发展，对它的姿态控制系统提出了更高的要
求，要求姿态控制系统简单、可靠。易于实现。高精度、较好的机动性和灵活性成为评价卫星
姿态控制系统优劣的重要标准。

[0004] 偏置动量卫星在寿命、可靠性、容错设计上，具有明显的优势，如可以在偏航轴没有测量信息的情况下正常对地控制，可以仅通过俯仰轴飞轮实现对地控制等。因此，现在很
多长寿命卫星都考虑设计偏置动量控制方式。但是，一般观测卫星的载荷为了同时兼顾重
访周期、幅宽和分辨率等功能，需要卫星定期进行姿态机动控制，这又需要卫星工作于零动
量控制方式。因此，最佳解决办法是卫星同时具备两种控制方式，并且二者之间能够安全可
靠地切换。

[0005] 由于观测卫星载荷对轨道高度的严格要求和卫星本身的安全性要求，喷气系统一般在初始入轨、轨道控制和姿态抢救时才选择使用，飞轮和磁力矩器成为卫星姿态控制的
主要执行机构。也有方法提出通过地磁力矩实现飞轮起旋，但是在起旋过程还是会对卫星
姿态产生较为明显的干扰，在此期间载荷一般不能正常工作。此处提出一种改进方法，同样
通过地磁力矩实现偏置动量轮起旋和消旋，通过对磁力矩前馈的方法，不仅利用了磁控的
最大能力，还能保证在起旋过程三轴继续保持正常飞轮控制，载荷正常成像。

发明内容

[0006] 针对现有技术存在的技术缺陷，本发明提供卫星动量轮在轨姿态无偏起旋和消旋的控制方法，实现卫星达到偏置角动量增加、减少过程中姿态无偏的效果。

[0007] 为解决上述技术问题，本发明是通过以下的技术方案实现卫星动量轮在轨姿态无偏差起旋和消旋的控制方法，包括如下步骤：

[0008] A、遥控注数设置卫星最终需要的目标角动量；

[0009] B、管理反作用飞轮的动量，控制所述卫星在期望中心转速范围内工作，当设置了起旋和消旋指令和俯仰轴目标转速后，增加至磁卸载目标角动量；

[0010] C、根据磁卸载力矩 和俯仰轴干扰补偿力矩，计算磁前馈力矩

[0011] D、将卫星控制器计算的输出力矩 磁前馈力矩 常规的解耦力矩和干扰补偿力矩 合成得到三轴指令力矩 控制周期内积分得到三轴卫星的指令
角动量，根据卫星的安装矩阵计算得到正在使用的卫星的指令转速。

[0012] 优选地，步骤A中所述目标角动量包括起旋时的俯仰轴偏置转速和/或角动量。

[0013] 优选地，步骤A中所述目标角动量包括消旋时的0转速和/或零动量。

[0014] 优选地，步骤A中所述卫星目标角动量为一次性设置。

[0015] 优选地，步骤C中所述磁前馈力矩 的计算公式为：

[0016]

[0017] 其中，Tdyconst为俯仰轴常值干扰力矩，根据卫星在轨实际需要进行补偿。

[0018] 优选地，所述磁卸载力矩 的计算公式为：

[0019]

[0020] 其中， 为磁场强度在卫星本体的分量，kx ,ky ,kz为卸载系数，为三轴卸载的角动量，正常情况下 当软件判断处于偏置
动量轮起旋模式时，自主将卸载角动量调整为

[0021]

[0022] 其中，hQX＝JΩ0为卫星俯仰轴希望达到的目标角动量。

[0023] 优选地，所述三轴指令角动量

[0024]

[0025] 其中， 为切换到起/消旋控制时三轴角动量。

[0026] 优 选 地 ，步 骤 D 中 所 述 三 轴 指 令 力 矩 的 计 算 公 式 为 ：

[0027] 优选地，所述输出力矩 的计算公式为：

[0028]

[0029] 其中，KPi、KIi和KDi(i＝x,y,z)为控制器参数， θc、ψc为控制用姿态角，ωox、ωoy、ωoz为根据惯性系角速度测量值、轨道及姿态角解算的星体轨道系角速度。

[0030] 优选地，将需要增加或消除的偏置角动量合成在所述磁卸载力矩 中，通过磁控力矩作用在星体上。

[0031] 本发明提供了卫星动量轮在轨姿态无偏差起旋和消旋的控制方法，可不使用喷气对大惯量卫星进行零动量控制与偏置动量控制方式平稳切换，通过磁控力矩作用在星体
上，同时通过前馈力矩的方式合成在飞轮的指令转速中抵消磁力矩对星体的影响，由此达
到偏置角动量增加和减少过程中姿态无偏差起旋和消旋的的效果。本发明基于更安全的轮
控模式下进行卫星偏置动量起、消旋控制，同时还可以在卫星动量变化过程中保证光学载
荷正常工作需要的姿态。

附图说明

[0032] 通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：

[0033] 图1示出了本发明的具体实施方式的，卫星动量轮在轨姿态无偏差起旋和消旋的控制方法；

[0034] 图2示出了本发明的实施例一的，同时具备零动量和偏置动量控制能力的飞轮模型；

[0035] 图3示出了本发明的实施例一的，计算磁前馈力矩 的步骤；

[0036] 图4示出了本发明的实施例一的，计算卫星三轴指令力矩 的步骤。

具体实施方式

[0037] 为了更好的使本发明的技术方案清晰的表示出来，下面结合附图对本发明作进一步说明。

[0038] 图1示出了本发明的具体实施方式的，卫星动量轮在轨姿态无偏差起旋和消旋的控制方法。具体地，包括如下步骤：

[0039] 首先，进入步骤S101，遥控注数设置卫星最终需要的目标角动量。

[0040] 然后，进入步骤S102，管理反作用飞轮的动量，控制所述卫星在期望中心转速范围内工作，当设置了起旋和消旋指令和俯仰轴目标转速后，增加至磁卸载目标角动量。

[0041] 再者，进入步骤S103，根据磁卸载力矩 和俯仰轴干扰补偿力矩，计算磁前馈力矩

[0042] 最后，进入步骤S104，将卫星控制器计算的输出力矩 磁前馈力矩 常规的解耦力矩 和干扰补偿力矩 合成得到三轴指令力矩 控制周期内积分得到
三轴卫星的指令角动量，根据卫星的安装矩阵计算得到正在使用的卫星的指令转速。

[0043] 图2示出了本发明的第一实施例的，同时具备零动量和偏置动量控制能力的飞轮模型，ObXbYbZb为卫星本体坐标系，偏置动量轮与反作用飞轮a、c、d中任两个组合可以进行
偏置动量起旋和消旋，实现偏置动量控制反作用飞轮a、b、c、d中任两个组合可以实现偏置
动量控制。具体地，包括如下步骤：

[0044] 某卫星偏置动量控制方式下，俯仰轴额定角动量‑60N.m.s，对应偏置动量轮中心转速‑5300rpm。卫星零动量工作方式下，偏置动量轮不工作或工作与0prm附近，当卫星需要
转至偏置动量控制方式时，

[0045] 首先，上注动量轮起旋指令包，该指令包中包含偏置动量轮目标转速Ω0，该值可以为‑5300rpm，也可以根据飞轮摩擦力矩的实际情况设置为其它值；

[0046] 然后，计算磁前馈力矩 如图3示出，具体步骤如下：

[0047] 1)进入步骤201，计算卫星三轴指令角动量 公式如下：

[0048]

[0049] 其中， 为切换到起/消旋控制时三轴角动量

[0050] 2)进入步骤202，计算磁卸载力矩 公式如下：

[0051]

[0052] 其中， 为磁场强度在卫星本体的分量，kx ,ky ,kz为卸载系数，为三轴卸载的角动量，正常情况下

[0053] 当软件判断处于偏置动量轮起旋模式时，自主将卸载角动量调整为

[0054]

[0055] 其中，hQX＝JΩ0为卫星俯仰轴希望达到的目标角动量，如‑60N.m.s或其它目标值。

[0056] 3)进入步骤203，计算磁前馈力矩 公式如下：

[0057]

[0058] 其中，Tdyconst为俯仰轴常值干扰力矩，根据卫星在轨实际需要进行补偿。

[0059] 最后，计算卫星三轴指令力矩 如图4示出，具体步骤如下：

[0060] 1)进入步骤301，计算卫星PID控制器力矩 公式如下：

[0061]

[0062] 其中，KPi、KIi和KDi(i＝x,y,z)为控制器参数， θc、ψc为控制用姿态角，ωox、ωoy、ωoz为根据惯性系角速度测量值、轨道及姿态角解算的星体轨道系角速度。

[0063] 2)进入步骤302，计算卫星三轴指令力矩

[0064] 三轴指令力矩包括卫星控制计算的输出力矩 磁前馈力矩 以及常规的解耦力矩 干扰补偿力矩 公式如下：

[0065]

[0066] 以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影
响本发明的实质内容。