一种基于自适应超扭滑模控制的平动点轨道保持方法与流程

[0001]
本发明属于航天器飞行动力学与控制领域，具体涉及一种基于自适应超扭滑模控制的平动点轨道保持方法。


背景技术：

[0002]
随着月球探测的新一轮高潮，充分开发地月空间、可持续更经济的利用月球资源成为主要课题之一。由于平动点轨道拥有特殊的动力学特性，其在深空探测领域具有重要的应用价值。然而，平动点轨道固有的不稳定特性要求航天器需要不断进行主动轨道保持，才能在附近长时间工作。轨道保持技术可以分为两大类：基于动力学特性的轨道保持技术和基于控制理论的轨道保持技术。相比于基于动力学特性的轨道保持技术，基于控制理论的轨道保持技术的优势在于它可以保证受控系统受到不确定因素干扰时的性能以及稳定性。但是，基于控制理论的轨道保持技术设计的轨道保持器性能保守、算法复杂、计算效率低，不适合于工程应用。


技术实现要素：

[0003]
针对现有技术中存在的问题，本发明提供了一种基于自适应超扭滑模控制的平动点轨道保持方法，用以解决上述技术问题。
[0004]
为了解决上述技术问题，本发明通过以下技术方案予以实现：
[0005]
一种基于自适应超扭滑模控制的平动点轨道保持方法，包括：
[0006]
在地月系统旋转坐标系下，建立地月系统中航天器平动点轨道运动的相对运动动力学方程；
[0007]
建立航天器平动点轨道运动的干扰加速度模型；
[0008]
根据干扰加速度模型得到干扰加速度，利用自适应超扭滑模控制方法和航天器平动点轨道运动的相对运动动力学方程，设计航天器平动点轨道保持的控制率。
[0009]
进一步地，所述建立地月系统中航天器平动点轨道运动的相对运动动力学方程，具体如下：
[0010]
定义：r
ref
＝[x
ref y
ref z
ref
]
t
表示参考轨迹，δr＝[δx δy δz]
t
表示受控轨迹r
c
与参考轨迹r
ref
之间的相对位置矢量，r
c
＝r
ref
+δr；
[0011]
航天器平动点轨道运动的相对运动动力学方程为：
[0012][0013][0014]
[0015]
其中，μ为地月系统质量参数；a＝[a
x a
y a
z
]
t
为航天器的控制加速度；
[0016][0017]
进一步地，所述建立航天器平动点轨道运动的干扰加速度模型具体包括：
[0018]
建立月球轨道偏心率扰动引起的干扰加速度模型为：
[0019][0020]
其中，d
e
为月球轨道偏心率扰动引起的干扰加速度；r
c
为月球圆轨道圆心与地月系统平动点间的距离；r
e
为月球椭圆轨道近地点与地月系统平动点间的距离；
[0021]
建立太阳引力引起的干扰加速度模型为：
[0022][0023]
其中，d
s
为太阳引力引起的干扰加速度；μ
s
为太阳的无量纲质量；r
se
为太阳到地球的位置矢量；r
sm
为太阳到月球的位置矢量；r
s
为太阳到航天器的位置矢量。
[0024]
进一步地，所述根据干扰加速度模型得到干扰加速度，利用自适应超扭滑模控制方法和航天器平动点轨道运动的相对运动动力学方程，设计航天器平动点轨道保持的控制率，具体为：
[0025]
定义自适应超扭滑模控制方法中的滑模变量为：
[0026][0027]
其中，k为参考轨迹与受控轨迹之间的位置误差的参数；c为参考轨迹与受控轨迹之间的速度误差的参数；
[0028]
根据航天器平动点轨道运动的相对运动动力学方程，并考虑月球轨道偏心率扰动引起的干扰加速度模型和太阳引力引起的干扰加速度模型得到的干扰加速度，滑模变量的一阶导数为：
[0029][0030]
其中，
[0031][0032]
[0033][0034]
其中，d＝d
e
+d
s
＝[d
x d
y d
z
]
t
为航天器平动点轨道运动的干扰加速度；
[0035]
根据自适应超扭滑模控制的基本原理，航天器平动点轨道保持的控制率设计为：
[0036][0037][0038][0039]
其中，a
x
,b
x
，a
y
,b
y
和a
z
,b
z
分别为在地月系统旋转坐标系中x,y和z方向三个分量的自适应增益。
[0040]
进一步地，a
x
,b
x
，a
y
,b
y
和a
z
,b
z
的求取方法为：
[0041][0042][0043][0044]
其中，ε
x
,λ
x
,r
x
,w
x
，ε
y
,λ
y
,r
y
,w
y
和ε
z
,λ
z
,r
z
,w
z
为实数。
[0045]
与现有技术相比，本发明至少具有以下有益效果：本发明相比于其他轨道保持控制器，设计的基于自适应超扭滑模控制器，可以在有干扰的情况下，实现目标轨道的高精度保持且算法收敛速度快；相比于传统滑模方法，本发明所设计的控制律可以实现连续控制输入并消除抖振，并具有算法简单、鲁棒性和自适应性强的优点。
[0046]
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
[0047]
为了更清楚地说明本发明具体实施方式中的技术方案，下面将对具体实施方式描
述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0048]
图1为地月系统旋转坐标系。
具体实施方式
[0049]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0050]
作为本发明的某一具体实施方式，一种基于自适应超扭滑模控制的平动点轨道保持方法，包括：
[0051]
步骤一、在地月系统旋转坐标系下，建立圆型限制性三体问题(cr3bp)地月系统中航天器平动点轨道运动的相对运动动力学方程，具体如下：
[0052]
定义：r
ref
＝[x
ref y
ref z
ref
]
t
表示参考轨迹，δr＝[δx δy δz]
t
表示受控轨迹r
c
与参考轨迹r
ref
之间的相对位置矢量，r
c
＝r
ref
+δr；
[0053]
其中，x
ref
表示r
ref
在旋转坐标系下的x轴方向的位置分量；y
ref
表示r
ref
在旋转坐标系下的y轴方向的位置分量；z
ref
表示r
ref
在旋转坐标系下的z轴方向的位置分量；δx表示δr在旋转坐标系下的x轴方向的位置偏差；δy表示δr在旋转坐标系下的y轴方向的位置偏差；δz表示δr在旋转坐标系下的z轴方向的位置偏差；
[0054]
如图1所示地月系统旋转坐标系中，航天器平动点轨道运动的相对运动动力学方程为：
[0055][0056][0057][0058]
其中，μ为地月系统质量参数；a＝[a
x a
y a
z
]
t
为航天器的控制加速度；
[0059][0060]
步骤二、建立航天器平动点轨道运动的干扰加速度模型具体如下：
[0061]
在深空环境中，航天器面临行星引力，执行机构非线性等各种扰动；在地月系统crtbp模型中，月球轨道偏心率扰动和太阳引力扰动为主要的干扰因素。因此，本发明只考
虑这两种干扰力，即建立月球轨道偏心率扰动引起的干扰加速度模型为：
[0062][0063]
其中，d
e
为月球轨道偏心率扰动引起的干扰加速度；r
c
为月球圆轨道圆心与地月系统平动点间的距离；r
e
为月球椭圆轨道近地点与地月系统平动点间的距离；
[0064]
建立太阳引力引起的干扰加速度模型为：
[0065][0066]
其中，d
s
为太阳引力引起的干扰加速度；μ
s
为太阳的无量纲质量；r
se
为太阳到地球的位置矢量；r
sm
为太阳到月球的位置矢量；r
s
为太阳到航天器的位置矢量。
[0067]
步骤三、根据干扰加速度模型得到干扰加速度，利用自适应超扭滑模控制方法和航天器平动点轨道运动的相对运动动力学方程，设计航天器平动点轨道保持的控制率，具体为：
[0068]
定义自适应超扭滑模控制方法中的滑模变量为：
[0069][0070]
其中，k为参考轨迹与受控轨迹之间的位置误差的参数；c为参考轨迹与受控轨迹之间的速度误差的参数；
[0071]
根据航天器平动点轨道运动的相对运动动力学方程，并考虑月球轨道偏心率扰动引起的干扰加速度模型和太阳引力引起的干扰加速度模型得到的干扰加速度，滑模变量的一阶导数为：
[0072][0073]
其中，
[0074][0075][0076][0077]
其中，d＝d
e
+d
s
＝[d
x d
y d
z
]
t
为航天器平动点轨道运动的干扰加速度；
[0078]
根据自适应超扭滑模控制的基本原理，航天器平动点轨道保持的控制率设计为：
[0079]
[0080][0081][0082]
其中，a
x
,b
x
，a
y
,b
y
和a
z
,b
z
分别为在地月系统旋转坐标系中x,y和z方向三个分量的自适应增益，具体的，a
x
,b
x
，a
y
,b
y
和a
z
,b
z
由下述方程得到：
[0083][0084][0085][0086]
其中，ε
x
,λ
x
,r
x
,w
x
，ε
y
,λ
y
,r
y
,w
y
和ε
z
,λ
z
,r
z
,w
z
为实数。
[0087]
本发明针对问题，基于圆型限制性三体问题模型，并考虑深空中普遍存在干扰力，设计了一种具有鲁棒性的自适应超扭滑模控制器，可以实现轨道保持的高精度快速响应。
[0088]
最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。