2022年毕设中期报告 .pdf

《2022年毕设中期报告 .pdf》由会员分享，可在线阅读，更多相关《2022年毕设中期报告 .pdf（12页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、毕 业 设 计 论 文 中 期 报 告题 目：基于飞轮摩擦观测器设计的卫星姿态控制方法研究专业飞行器设计与工程学生解延浩学号1111820131 指导教师吴宝林日期2015 年 3 月 13 日哈尔滨工业大学教务处制精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页，共 12 页1. 绪论对于要求三轴稳定的小卫星姿态控制系统，当前主要有喷气、 反作用飞轮和力矩陀螺三种控制方式。而其中，反作用飞轮由于具有不消耗推进剂、控制精度高、系统简单灵活可靠并且可以实现整星零动量等优点，而得到广泛的运用。但在转速较低时，由于摩擦力矩的存在，尤其是在转速过零时

2、摩擦力的大小、方向均产生突变，导致反作用轮的输入信号和输出力矩间的线性关系被严重破坏，大大限制了该方案的实际应用。因此，如何改善反作用飞轮的低速摩擦性能成为了提高飞轮姿态控制精度的关键课题。比照当前的国内外主要的诸如自适应、变结构、H等小卫星姿态控制系统，虽然在精度、稳定性和机动能力上都已到达了较佳的性能，但针对于反作用飞轮的低速摩擦补偿却研究不多，也尚未有较为有效可靠的方案。当前对于低速摩擦补偿，国内外学者主要提出了以下方式：1直接补偿：将摩擦力都视作库伦摩擦，并根据库伦摩擦力模型，在飞轮输入信号上叠加一个与转速同向的偏置信号，用以抵消摩擦力。 这样虽然可以改善一部分飞轮性能， 但由于模型的

3、不精确而且所受摩擦力不全是库伦摩擦，导致在反作用飞轮转速过零时补偿效果减弱甚至失效； 2变增益：在低速过零时提高飞轮转速的反馈增益，即提高了系统转速跟踪的精度， 但过高的反馈增益会导致系统能耗增加而且削弱系统稳定性甚至产生极限环振荡现象； 3高频线性化：当反作用轮转速低于某个阈值时，将高频正弦震颤信号叠加进输入信号中，使原本不连续的摩擦特性得到了较好的线性化，但此方法受限于附加的高频震颤信号和飞轮的摩擦频率特性。4变结构控制：将摩擦力在飞轮过零时的突变视作模型的不确定性，并假设不确定上界已知。这种方式虽然有效，但可能导致较大的震颤和控制量，故仍具有一定的保守型和改良空间。2. 研究内容本文比照

4、了现有的摩擦补偿方案，决定设计摩擦补偿观测器，并采用变结构控制的方式，使系统按照预定“滑动模态”的状态轨迹运动，改善飞轮的低速摩擦性能。结合现有的知识技能，需要进一步研究的内容主要包括：1)反作用飞轮的低速摩擦与润滑机理和较准确的低速摩擦模型、收敛的飞轮摩擦补偿状态观测器设计及其应用；2)变结构控制的基础理论与应用以及包含摩擦补偿观测器的变结构姿态控制系统的设计；3)考虑摩擦的变结构姿态控制系统的Simulink 仿真与结果评价及进一步减小控制量和震颤的优化方案，并在仿真中验证优化的效果。3. 研究方案及进度安排1)查阅有关飞轮低速摩擦和润滑的相关文献或者进行实验测定数据，给出较为精确精选学习

5、资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页，共 12 页的低速过零时的摩擦模型，并将其合理的连续化，再由现代控制理论得出观测器状态方程，并证明其收敛性。2)学习滑模变结构控制的相关理论，并研究将之运用于卫星姿态控制的方式。3)根据卫星姿态动力学模型，首先建立不失真实性的三通道姿态计控制系统包含不确定的干扰因素，并在Simulink 中搭建模块， 进行电脑仿真以观察低速摩擦对姿态控制系统的影响。4)在已有基础上，以减小姿态误差和扰动为目标，设计基于状态观测器的变结构卫星姿态控制系统，并在同等条件下进行仿真，并与之前的仿真结果进行比照。5)结合工程

6、实际方法，从多角度尝试对这种控制策略进行进一步优化，目前尚在论证中的方案是在低速时增加系统阻尼或采用极点连续变化的方式，以改善其震颤特性，并进行仿真验证。4. 进度完成情况序号阶段及内容起讫日期阶段成果形式完成情况1 通过导师的参考方向和自己的兴趣确定研究课题2015年 3月 3 日至2015年 3月 5 日选题方案已完成2 通过网络等途径搜集相关资料，并进行加工整理，总结研究现状，发现问题2015年 3月 6 日至2015 年 3 月 11 日资料总结已完成3 确定研究课题，梳理资料后撰写开题报告，认真思考导师给的意见和建议，完善开题报告。进行开题报告答辩2015 年 3 月 12 日至20

7、15 年 3 月 19 日开题报告已完成4 摩擦模型的建立及连续化，开始详细学习滑模变结构控制的相关理论，并研究将之运用于卫星姿态控制系统的方案2015 年 3 月 20 日至2015 年 3 月 25 日较为真实的连续化摩擦模型已完成5 状态观测器方程的建立2015 年 3 月 26 日至2015年 4月 5 日收敛的状态观测器已完成6 根据卫星姿态动力学模型，建立不失真实性的三通道卫星姿态计控制系统2015年 4月 6 日至2015 年 4 月 15 日考虑摩擦的三通道姿态计控制系统已完成精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页，

8、共 12 页7 学习卫星姿态动力学模型和 matlab 仿真方法，在电脑中搭建控制系统，并进行仿真包含摩擦等不确定的干扰因素2015 年 4 月 16 日至2015 年 4 月 25 日未经优化的系统仿真曲线已完成8 总结研究结果，撰写中期报告，进行中期答辩2015 年 4 月 25 日至2015 年 4 月 30 日中期报告进行中9 继续学习变结构控制，并开始设计基于状态观测器的变结构卫星姿态控制系统，进行仿真2015年 5月 1 日至2015 年 5 月 20 日仿真曲线未完成10 尝试改良控制策略，改善控制效果2015 年 5 月 21 日至2015年 6月 6 日优化的仿真曲线未完成1

9、1 总结研究结果，撰写论文。2015年 6月 7 日至2015 年 6 月 20 日最终论文未完成表 1 5. 已完成的研究工作及成果1)卫星姿态动力学模型建模含飞轮卫星的姿态动力学模型即卫星在内外力矩的共同作用下绕质心转动的规律，包括整星运动和星体各部分的相对运动。出于简化考虑，本文采用单刚体假设，即不考虑星体挠性和液体晃动对姿态运动的影响，将卫星视作简单刚体进行计算。卫星通常由星体和转子组成，以卫星质心O为原点，并以某需要的特征方向为坐标轴建立卫星本体坐标系Oxyz。Is为卫星本体不含飞轮在Oxyz系内的转动惯量，IR为飞轮转动惯量。转动角速度在本体系中的分量依次为x、y、z，则可记作=(

10、x, y, z)T，飞轮相对于主体的角速度记作=(1, 2, 3) ，则克得到卫星相对其质心 O的动量矩：SR= IIH(1) 设卫星受环境干扰力矩Te，则由动量矩定理有：eH HT(2) 代公式(2) 入公式 (1) ，可得方程：()SRSReIIIIT(3) 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页，共 12 页其中，000zyzxyx为角速度的叉乘矩阵，RI 为飞轮对本体的实际输出力矩，记作 Tu故得到卫星姿态动力学方程：1I (II)sSRueTT(4) 同理，对于飞轮，其动力学方程为：RIuCfTTT(5) 其中 TC为根据

11、控制方程求出的应有控制力矩，Tf为摩擦力矩，主要包括轴承固体间摩擦力矩和润滑剂的粘滞摩擦力矩。对于卫星姿态运动的描述，考虑到防止奇异，采用姿态四元数q 描述卫星姿态。以四元数形式表示的卫星姿态运动学方程：03 312Tqq Iqq(6) 其中：00123qqqqqqq为卫星姿态四元数；2)摩擦模型建立及状态观测器设计反作用飞轮在轴承的支撑下转动，不仅仅受到电机的驱动力矩，也受到轴承的摩擦力矩，二者合力矩才是飞轮的和外力矩，即姿态控制系统真正的控制力矩：-rcfTTT(7) 其中， Tr 为卫星实际控制力矩， Tc 为根据控制律得出的应受控制力矩，Tf主要包括轴承固体间摩擦力矩、 库伦摩擦力矩和

12、润滑剂的粘滞摩擦力矩，粘性摩擦力矩与转速成正比，比例系数为 Kv，而固体摩擦则参考Dahl 于 1968年提出的固体间摩擦模型，即Dahl模型。该模型与实际摩擦情况较为相似，已被广泛应用，而库伦摩擦则可视作定值。Dahl 模型如下：0时：21vKsignSfksTTT(8)0时：mcmfcmmcmcmmTT TT =T -T-T T TT -T-T(9)精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页，共 12 页(8)式中 Kv为粘性摩擦力矩系数， Tm为最大静摩擦力矩， Tc为电机驱动力矩， Tk为库伦摩擦力矩， Ts为 Stribeck

13、摩擦力矩， s为 Stribeck摩擦特征速率。工程中常常取一个极小的正常数 ，定义：停滞区为 。由于粘性部分的建模较为简单，故通常将其与线性系统模型合为一部分，则其余部分为：22tanh111vKsignSSfnfkkssTTTTTT(10) 此处为便于计算，将sign 函数用 tanh 函数代替，其中01。由此可将飞轮动力学模型线性化：112()tanh11cvfcvJKJKSksT TTTT(11) 相应的状态方程为：12()(sign()JDfcfff0TTT TT(12) 可以得出 =0，Tf=0是该式的一个稳定平衡点，在平衡点将方程展开，可得：1120( )0fJDJoTffxTT

14、(13) 这里，x=，TfT，o(x) 表示高阶小量。在这里可以根据上式的线性部分判断系统稳定性：在 =Tf=0点系统是稳定的，即低速摩擦有令飞轮趋于停止的特点。将(13) 上式进一步线性化，即可得到线性系统：10JDfcfTTT(14) 考虑到系统中 可以直接测量，故设计状态观测器：111?JDkkfcfT T T(15) 其中 k1与 k2分别为观测器增益，调节增益可实现对线性系统的状态观测，且经证明此观测器同样可以实现对非线性系统的状态观测。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页，共 12 页3)基于误差四元数的PD控制系统设

15、计a)为观察摩擦对于飞轮姿态控制系统的影响，并进一步熟悉姿态控制系统设计，设计了基于误差四元数和加速度误差的姿态控制系统，系统框图如下图 1 进行控制律设计必须先明确误差信号的表示形式，出于防止奇异的考虑，本文选用姿态四元数表征卫星姿态，而姿态误差则有误差四元数表示，其计算方法为：*eadqqq(16) 其中 qe为误差四元数， qa为实际姿态四元数， qd为目标姿态四元数。 q*为四元数的共轭，四元数乘法计算：0000.,abababababbaq qqqqqq qqqqq(17) 以四元数形式表示的卫星姿态运动学方程：03 312Tqq Iqq(18) 系统采用经典的 PD控制，即以误差四

16、元数和误差角速度作为控制器的输入量，控制方程为：=dssdRpdIIIKKuddeeTq(19) 其中dssdRIIIdd为考虑姿态跟踪的前馈量。且可以证明该控制律可以实现系统的渐进稳定性。4)控制系统的 matlab 实现与仿真matlab 仿真程序设计以 yf =qe1,qe2,qe3, 1, 2, 3, 1, 2, 3 为每一步积分的初值，其中qe1,qe2,qe3为上一控制周期末的姿态与期望姿态间误差四元数的矢量部分，1, 2,3为角速度误差，1, 2, 3为飞轮转速，根据系统动力学与运动学模型：精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -

17、第 7 页，共 12 页103 31I(II)12()sSRufTRfqq IqIuTTqTT(20) 建立 oed45 积分函数。每一个控制周期内控制力矩不变且与上一周期末的姿态误差、姿态角速度相关，积分迭代后即可得出卫星在仿真时间内的姿态运动情况。仿真结果参考“实验三号”小卫星的有关数据，设置：控制周期为 t=250ms 仿真时间设 ts=22000s 迭代步数 n=88000步星体转动惯量矩阵：212.490.670.060.6713.850.060.060.0615.75SIkg m飞轮转动惯量矩阵：20.00230000.00230000.0023RIkg m轨道角速度00.001/

18、rads初始姿态：0.70710.30.40.5Tiq星体初始角速度0.010.0050.0033/Tradsi飞轮初始角速度000/radsri控制系统增益：3pK ,18dK粘性摩擦系数：56.4 10/vKN m s rad库伦摩擦力矩：42.5 10kTN mStribeck摩擦力矩：42.5 10sTN mStribeck摩擦特征速率常数：2.5/Srads仿真曲线及数据：有摩擦时的姿态误差：精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页，共 12 页 轴：0.20.40.60.811.21.41.61.822.2x 104-20

19、020时 间 （ s）飞轮角速度（r/min）飞 轮 角 速 度 曲 线0.20.40.60.811.21.41.61.822.2x 104-202时 间 (s)姿态角误差(deg)姿 态 误 差0.20.40.60.811.21.41.61.822.2x 104-0.0500.05时 间 (s)摩擦摩 擦图 2 轴：0.20.40.60.811.21.41.61.822.2x 104-20020时 间 （ s）飞轮角速度（r/min）飞 轮 角 速 度 曲 线0.20.40.60.811.21.41.61.822.2x 104-202时 间 (s)姿态角误差(deg)姿 态 误 差0.20.

20、40.60.811.21.41.61.822.2x 104-0.0500.05时 间 (s)摩擦摩 擦图 3 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页，共 12 页 轴：0.20.40.60.811.21.41.61.822.2x 104-20020时 间 （ s）飞轮角速度（r/min）飞 轮 角 速 度 曲 线0.20.40.60.811.21.41.61.822.2x 104-202时 间 (s)姿态角误差(deg)姿 态 误 差0.20.40.60.811.21.41.61.822.2x 104-0.0500.05时 间 (s

21、)摩擦摩 擦图 4 无摩擦时的姿态误差 轴：0.20.40.60.811.21.41.61.822.2x 104-20020时 间 （ s）飞轮角速度（r/min）飞 轮 角 速 度 曲 线0.20.40.60.811.21.41.61.822.2x 104-101时 间 (s)姿态角误差(deg)姿 态 误 差0.20.40.60.811.21.41.61.822.2x 104-101时 间 (s)摩擦摩 擦图 5 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页，共 12 页 轴：0.20.40.60.811.21.41.61.822.

22、2x 104-20020时 间 （ s）飞轮角速度（r/min）飞 轮 角 速 度 曲 线0.20.40.60.811.21.41.61.822.2x 104-101时 间 (s)姿态角误差(deg)姿 态 误 差0.20.40.60.811.21.41.61.822.2x 104-101时 间 (s)摩擦摩 擦图 6 轴：0.20.40.60.811.21.41.61.822.2x 104-20-100时 间 （ s）飞轮角速度（r/min）飞 轮 角 速 度 曲 线0.20.40.60.811.21.41.61.822.2x 104-101时 间 (s)姿态角误差(deg)姿 态 误 差0

23、.20.40.60.811.21.41.61.822.2x 104-101时 间 (s)摩擦摩 擦图 7 由仿真结果可见，无摩擦力矩作用时，卫星姿态稳态误差小于0.008deg，而加入摩精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页，共 12 页擦力矩后姿态误差最大到达deg。因此，必须在飞轮的低速区内设法补偿摩擦力矩，尤其是静摩擦力矩，以提高姿态控制性能。6. 后期拟完成的研究工作及进度安排依照原计划，如表格所示：1 总结研究结果， 撰写中期报告，进行中期答辩2015年 4 月 25日至2015年 4 月 30日中期报告2 继续学习变结

24、构控制， 并开始设计基于状态观测器的变结构卫星姿态控制系统， 进行仿真2015 年 5 月 1 日至2015年 5 月 20日仿真曲线3 尝试改良控制策略， 改善控制效果2015年 5 月 21日至2015 年 6 月 6 日优化的仿真曲线4 总结研究结果，撰写论文。2015 年 6 月 7 日至2015年 6 月 20日最终论文表 2 7. 存在的问题与困难对于滑模变结构控制的理论及运用了解不够，需继续学习对于状态观测器在matlab 中的实现及在变结构控制系统中的运用缺乏更多控制系统性能评价手段，对系统的性能评价不够全面8. 论文按时完成的可能性可按时完成精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 12 页，共 12 页