毕设中期报告

毕业设计（论文）中期报告

题目：基于飞轮摩擦观测器设计的卫星姿态控制方法研究

专业飞行器设计与工程

学生解延浩

学号1111820131

指导教师吴宝林

日期2015年3月13日

哈尔滨工业大学教务处制

1.绪论

对于要求三轴稳定的小卫星姿态控制系统，当前主要有喷气、反作用飞轮和力矩陀螺三种控制方式。而其中，反作用飞轮由于具有不消耗推进剂、控制精度高、系统简单灵活可靠并且可以实现整星零动量等优点，而得到广泛的运用。但在转速较低时，由于摩擦力矩的存在，尤其是在转速过零时摩擦力的大小、方向均产生突变，导致反作用轮的输入信号和输出力矩间的线性关系被严重破坏，大大限制了该方案的实际应用。因此，如何改善反作用飞轮的低速摩擦性能成为了提高飞轮姿态控制精度的关键课题。

等小卫星姿态控制系统，虽然在精对比当前的国内外主要的诸如自适应、变结构、H

∞

度、稳定性和机动能力上都已达到了较佳的性能，但针对于反作用飞轮的低速摩擦补偿却研究不多，也尚未有较为有效可靠的方案。当前对于低速摩擦补偿，国内外学者主要提出了以下方式：（1）直接补偿：将摩擦力都视作库伦摩擦，并根据库伦摩擦力模型，在飞轮输入信号上叠加一个与转速同向的偏置信号，用以抵消摩擦力。这样虽然可以改善一部分飞轮性能，但由于模型的不精确而且所受摩擦力不全是库伦摩擦，导致在反作用飞轮转速过零时补偿效果减弱甚至失效；（2）变增益：在低速过零时提高飞轮转速的反馈增益，即提高了系统转速跟踪的精度，但过高的反馈增益会导致系统能耗增加而且削弱系统稳定性甚至产生极限环振荡现象；（3）高频线性化：当反作用轮转速低于某个阈值时，将高频正弦震颤信号叠加进输入信号中，使原本不连续的摩擦特性得到了较好的线性化，但此方法受限于附加的高频震颤信号和飞轮的摩擦频率特性。（4）变结构控制：将摩擦力在飞轮过零时的突变视作模型的不确定性，并假设不确定上界已知。这种方式虽然有效，但可能导致较大的震颤和控制量，故仍具有一定的保守型和改进空间。

2.研究内容

本文对比了现有的摩擦补偿方案，决定设计摩擦补偿观测器，并采用变结构控制的方式，使系统按照预定“滑动模态”的状态轨迹运动，改善飞轮的低速摩擦性能。结合现有的知识技能，需要进一步研究的内容主要包括：

1)反作用飞轮的低速摩擦与润滑机理和较准确的低速摩擦模型、收敛的飞轮摩擦补偿状态观测器设计及其应用；

2)变结构控制的基础理论与应用以及包含摩擦补偿观测器的变结构姿态控制系统的设计；

3)考虑摩擦的变结构姿态控制系统的Simulink仿真与结果评价及进一步减小控制量和震颤的优化方案，并在仿真中验证优化的效果。

3.研究方案及进度安排

1)查阅有关飞轮低速摩擦和润滑的相关文献或者进行实验测定数据，给出较为精确

的低速过零时的摩擦模型，并将其合理的连续化，再由现代控制理论得出观测器状态方程，并证明其收敛性。

2)学习滑模变结构控制的相关理论，并研究将之运用于卫星姿态控制的方式。

3)根据卫星姿态动力学模型，首先建立不失真实性的三通道姿态计控制系统（包含不确定的干扰因素），并在Simulink中搭建模块，进行计算机仿真以观察低速摩擦对姿态控制系统的影响。

4)在已有基础上，以减小姿态误差和扰动为目标，设计基于状态观测器的变结构卫星姿态控制系统，并在同等条件下进行仿真，并与之前的仿真结果进行对比。

5)结合工程实际方法，从多角度尝试对这种控制策略进行进一步优化，目前尚在论证中的方案是在低速时增加系统阻尼或采用极点连续变化的方式，以改善其震颤特性，并进行仿真验证。

4.进度完成情况

5. 已完成的研究工作及成果

1)

卫星姿态动力学模型建模（含飞轮）

卫星的姿态动力学模型即卫星在内外力矩的共同作用下绕质心转动的规律，包括整星运动和星体各部分的相对运动。出于简化考虑，本文采用单刚体假设，即不考虑星体挠性和液体晃动对姿态运动的影响，将卫星视作简单刚体进行计算。

卫星通常由星体和转子组成，以卫星质心O 为原点，并以某需要的特征方向为坐标轴建立卫星本体坐标系Oxyz 。I s 为卫星本体（不含飞轮）在Oxyz 系内的转动惯量，I R 为飞轮转动惯量。转动角速度ω在本体系中的分量依次为ωx 、ωy 、ωz ，则可记作ω=(ωx ,ωy ,ωz )T

，飞轮相对于主体的角速度记作Ω=(Ω1,Ω2,Ω3)，则克得到卫星相对其质心O 的动量矩：

S R =I I +H ωΩ (1) 设卫星受环境干扰力矩T e ，则由动量矩定理有：

+⨯=e H ωH T

(2) 代公式(2)入公式(1)，可得方程：

()S R S R e I I I I ω⨯+++=ωωωΩT

(3)

其中，00

0z y z

x y x

ωωωωωωω⨯⎡⎤-⎢

⎥

=-⎢⎥⎢⎥-⎣

⎦

为角速度的叉乘矩阵，R I ω为飞轮对本体的实际

输出力矩，记作T u 故得到卫星姿态动力学方程：

1I (I I )s S R u e ωω-⨯⨯

=---+ωωΩT T

(4) 同理，对于飞轮，其动力学方程为：

R I ==-u C f ΩT T T

(5)

其中T C 为根据控制方程求出的应有控制力矩，T f 为摩擦力矩，主要包括轴承固体间摩擦力矩和润滑剂的粘滞摩擦力矩。

对于卫星姿态运动的描述，考虑到避免奇异，采用姿态四元数q 描述卫星姿态。 以四元数形式表示的卫星姿态运动学方程：

0331

2T q q I q ⨯

⨯⨯⎡⎤-⎢

⎥=⎢⎥+⎣⎦

q ω (6)

其中：00

123q q q q q ⎡⎤

⎢⎥⎡⎤

⎢⎥==⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦q q 为卫星姿态四元数；

2)

摩擦模型建立及状态观测器设计

反作用飞轮在轴承的支撑下转动，不仅仅受到电机的驱动力矩，也受到轴承的摩擦力矩，二者合力矩才是飞轮的和外力矩，即姿态控制系统真正的控制力矩：

-=r c f T T T

(7)

其中，Tr 为卫星实际控制力矩，Tc 为根据控制律得出的应受控制力矩，T f 主要包括轴承固体间摩擦力矩、库伦摩擦力矩和润滑剂的粘滞摩擦力矩，粘性摩擦力矩与转速成正比，比例系数为K v ，而固体摩擦则参考Dahl 于1968年提出的固体间摩擦模型，即Dahl 模型。该模型与实际摩擦情况较为相似，已被广泛应用，而库伦摩擦则可视作定值。

Dahl 模型如下：

0≠Ω 时：()21v K sign ⎛⎫ ⎪

⎪=++⨯ ⎪⎛⎫ ⎪+ ⎪

⎪⎝⎭⎝

⎭S f k

s T T ΩT ΩΩΩ (8)

0=Ω 时：⎧⎪⎨⎪⎩m

c m f c m m c m c m m

T T >T T =T -T -T

(9)