《航天器的姿态与轨道最优控制_1.docx》由会员分享，可在线阅读，更多相关《航天器的姿态与轨道最优控制_1.docx（14页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、航天器的姿态与轨道最优控制航天器的姿态与轨道最优控制董丽娜唐晓华吴朝俊司渭滨(第八小组)西安交通大学电气工程学院，陕西省，西安市710049【摘要】从航天器的轨道运动学方程出发,运用线性离散系统最优控制理论,提出了一种用于航天器轨道维持与轨道机动的最优控制方法,建立了相关的最优控制模型并给出了求解该模型的算法。仿真计算结果表明,本文提出的最优控制方法是正确和可行的。【关键词】航天器轨道保持轨道机动最佳控制OptimalControlofSpacecraftStateandOrbitDongLiNa，TangXiaoHua，WuChaoJun，SiWeiBinEESchoolofXianJiao
2、tonguniversity，Xian,Shannxiprovince,710049【Abstract】Thispaperprovidesanewoptimalcontrolmethodfororbitalmaintenanceandmaneuver,whichbeginswiththekineticsequationofspacecraftandisbasedonthelineardiscreteoptimalcontroltheory,establishestherelativeoptimalcontrolmodelandgivesitssolution.Thesimulationresu
3、ltsshowthatthegivenoptimalcontrolmethodinthispaperiscorrectandfeasible.【Keyword】Spacecraft，Orbitalkeeping，Orbitalmaneuver，Optimalcontrol1引言一般地,常见的航天器有:运载火箭、人造卫星、载人飞船、宇宙飞船、空间站等。宇宙飞)。对于无人航天器由一下几部分组成：(1)有效载荷(2)构造(3)动力系统(4)温度控制系统(5)姿态和轨道控制(6)电源系统和遥测、遥控、跟踪系统星载部分其中姿态和轨道控制是航天器运行经过中比拟重要的技术。包括被动式和主动式。(1)被动式：
4、利用航天器本身的动力特性和自然环境力矩来控制姿态。a.不需消耗航天器上的能源b.构造简单，适用于较长寿命的航天器c.控制精度不高,常用的有自旋稳定，重力梯度稳定，磁稳定等。(2)主动式：当控制器送来信号后则启动执行机构，控制卫星姿态的执行机构。a.以喷气三轴控制卫星上装有轴向，横向和切向喷嘴b.以飞轮为主的三轴控制即角动量控制姿态和轨道控制必要性:(1)卫星进入轨道后，长时间的运行经过中会碰到各种干扰，即便极稀薄的空气也会构成空气阻力，阻力固然很小但长久的作用也会使卫星速度变慢，偏离预定的轨道。(2)卫星的天线和观察设备要卫星有一定的指向，当偏离后就需要调整它的姿态。(3)返回式卫星在轨工作完
5、成后，也要调整卫星姿态。使其发动机喷口对准预定方向，点火产生推力，使卫星进入轨道，然后返回地面。由此可见，航天器轨道维持和轨道机动是轨道力学研究的重要问题。工程实践中,轨道维持通常是按轨道平面维持和轨道形状维持等多个方面分别进行的,这不但繁琐而且不易实现燃料消耗的最小化。如何才能以最小的能耗实现航天器的轨道维持呢?对于已偏离其标称轨道(由于入轨误差和空间各种扰动因素的影响)、同时又肩负着变轨责任的在轨航天器而言,又该怎样调整预先设定好的各变轨点速度脉冲,最终以最小的能耗代价飞抵原定目的轨道或临时变更的新目的轨道呢?针对这些问题,本文运用线性离散系统最优控制理论对轨道维持与机动问题进行了研究,提
6、出了轨道维持与机动的最优控制模型及其求解算法。该模型的提出,从最优控制的角度出发对轨道维持与机动问题展开研究。2最优控制模型的建立与求解设计目的：当航天器在运行工作中，其中某一动量飞轮发生故障或失效，航天器能保持有效的姿态控制和任意的定位。设计要求：利用最优控制方法提出求解带有2个动量飞轮航天器的姿态控制算法,通过数值仿真,表明该方法对航天器姿态控制是有效的分析思路：本课题考虑航天器在2个动量飞轮作用下，当系统角动量为0时，将航天器系统的姿态控制问题转化为无漂移系统的运动规划问题,利用最优控制方法确定动量飞轮控制输入规律，以到达航天器主体的期望姿态。2.1系统建模如图1所示以系统质点O为原点建
7、立相对惯性空间平动的坐标系。i分别是iO到系统总质心O的矢径两个飞轮相对系统总质心的位置:其中：又由牛顿定律,航天器系统相对点的动量矩可表示为:当航天器系统起始动量矩为零。上式可化作为:以上各式中：为航天器的绝对角速度矢量。iI(i=0,1,2)分别为主刚体和2个飞轮的惯性张量。10112022,dbdb=+=+121122,00xxyybbbbbb?=?2122000R(,)JJ()JII(IJ)JbbiiiiiiiiiiiiiHbj=T=+?=+-?=22111(JJ)biiiii-=-+i(i=1,2)分别为飞轮绕的转动角。Ji(i=1,2)分别为动量飞轮绕的惯量矩。又由于能够用卡尔丹角
8、及其导数表示为:代入上式,则有:其中：(,)q=2.2控制方法将动量飞轮相对转动角速度i取作输入变量,记作u,定义航天器主刚体的位形(,)qT=状态变量,则系统的状态方程为：B(x)uq=其中：22-111211B(x)=L(),u=(,)iiiiJJb-T=-+根据最小能量控制原理,选择航天器动量飞轮转动耗散能作为最优控制指标,性能指标函数为：0(),dTJuuut=?在式中()ut为Hilbert空间2L的可测向量函数。实际计算时,只需考虑有限维的情况,则u可表示为Fourier基向量1Niie=的线性组合1NiiiUe=其中()1,2,iiN=为函数在1Niie=基上的投影。将视作新的控
9、制变量,考虑终端约束条件,指标函数()Ju写为：221(,)()NifiJxTx=+-为罚因子,能够证实,当,时,式()Ju与(,)J有一样的最优解。()XT是由控制输入u给定在tT=时的解,为RN的函数,记作()()XTf=。coscossin0coscoscos0sin01?=-?221111L(JJ)iiiiiqb-=-+给定N和,则上式可变为：(),()fJfx=+-因而寻找控制输入使式()Ju为最小值的问题转化为寻找式指标函数()J为最小值的问题。能够用无约束算法对变换后的系统求最优控制,将带入即得到原问题的最优控制u。利用牛顿法,在极小点附近用指标函数()Ju的Taylor多项式展
10、开,取二阶近似并令9()/90nJ+=,得到牛顿法的迭代公式1212nnnnJJ-+?=-?其中：2()nnfnJfxT?=+A-?()24212()infiiinJfxT=?=E+AA+-H?此式中A为f的Jacobi矩阵,iH为f的分量if的Hesse矩阵,E为单位阵。因()EAAT+恒为正定,指标函数()J又具有平方和形式,上式可用修正牛顿法,即Gauss-Newton迭代公式表示11()nnnnffx+-TT?=-E+AA+A-?式中1/=,为步长因子,013模型仿真求解考虑航天器利用2个动量飞轮转动控制姿态运动,航天器系统质量几何参数为：121212021202122b100b(0,
11、1,0)0.2,0.5Kgm500,5Kg,Idiag(86.215,85.07,113.565)KgmIdiag(0.5,0.25,0.25)KgmIdiag(0.25,0.5,0.25)KgmddmjjmKgmmTT=?=?=?=?2，设航天器初始和终端位形分别为：000q0,q00/6f?=?仿真试验中选取10个Fourier正交基矢量,其中61()iiet=为：1234560.5sincossincossine,e,e,e,e,e000000ttttt?=?107()iiet=分别由上式各基矢量行轮换得到。上例通过19次迭代到达最优指标值,19=61.67,误差精度为210-。图2为航
12、天器动量飞轮相对转动的最优控制输入规律,图3为航天器主刚体姿态运动优化时间历程。图2动量飞轮相对转动的最优控制输入规律文档视界航天器的姿态与轨道最优控制航天器的姿态与轨道最优控制0q0/6f?=?的有效控制，进而讲明航天器在系统角动量为零的情况下，能够利用2个动量飞轮对刚体航天器的姿态进行控制。非线性最优控制算法在解决航天器姿态控制问题中是有效的。从仿真算例可看出刚体航天器可控并能完全到达指定终端位形。参考文献：1刘延柱.航天器姿态动力学.北京:国防工业出版社,1995.2黄圳圭.航天器姿态动力学.长沙,国防科技大学出版社,19973杨嘉墀.航天器轨道动力学与控制.北京,宇航出版社,20014
13、中国人民当兵的人总装备部军事训练教材编辑工作委员会.航天器飞行控制与仿真.北京,国防工业出版社2004.25CrouchPE.Spacecraftattitudecontrolandstabilization:applicationofgeometriccontroltheorytorigidbodymodelsJ.IEEETransactionsonAutomaticControl,1984,29(4):87-956ChrisFernandes,leonidGurvits,ZexiangLi,Member,-OptimalNonholonomicMotionPlanningForaSystemofCoujpledRigidBodies.450-463第八组人员名单：