CN107065547A - 一种基于零空间方法的非合作目标自主交会策略 - Google Patents

一种基于零空间方法的非合作目标自主交会策略 Download PDF

Info

Publication number
CN107065547A
CN107065547A CN201710225648.2A CN201710225648A CN107065547A CN 107065547 A CN107065547 A CN 107065547A CN 201710225648 A CN201710225648 A CN 201710225648A CN 107065547 A CN107065547 A CN 107065547A
Authority
CN
China
Prior art keywords
task
msub
mrow
noncooperative target
kernel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710225648.2A
Other languages
English (en)
Other versions
CN107065547B (zh
Inventor
袁建平
李琪
朱战霞
袁源
马卫华
何立文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201710225648.2A priority Critical patent/CN107065547B/zh
Publication of CN107065547A publication Critical patent/CN107065547A/zh
Application granted granted Critical
Publication of CN107065547B publication Critical patent/CN107065547B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Abstract

本发明公开了一种基于零空间方法的非合作目标自主交会策略,针对非合作目标在轨服务任务,首先对自主交会过程进行任务分析,将整个交会过程分为非合作目标对接抓捕点跟踪任务和碰撞规避任务,对上述两种任务进行优先级排序,并对每种基本任务所对应的零空间进行求解,最后再根据零空间方法对这两种任务进行综合,将低优先级任务投影到高优先级任务的零空间,从而保证了交会轨迹的安全性。本发明所提供的基于零空间方法的自主交会策略是一种分层递阶控制的方法,它能够确保高优先级任务最先得到完成,从而避免不同优先级任务之间发生冲突。基于零空间方法的非合作目标自主交会策略能够同时实现对非合作目标对接抓捕点的跟踪以及碰撞规避要求。

Description

一种基于零空间方法的非合作目标自主交会策略
技术领域
本发明属于航天器动力学领域,具体涉及一种基于零空间方法的非合作目标自主交会策略。
背景技术
随着空间科学技术的快速发展,人类对于空间的探索活动也在近几十年急剧增加。迄今为止,人类已经将数千颗人造卫星成功送入轨道,用于执行各种各样的空间任务。然而就在全人类正在为自己所取得的巨大进步暗自高兴时,人类也逐渐认识到了空间环境正在快速恶化。目前,空间碎片正在以一种“凯斯勒”效应在快速增长,同时现有的卫星都不具备可维修的能力,一旦出现故障就只能报废。在这种情况下,空间在轨服务概念广泛提及,与此同时在轨服务技术也成为了当前以及未来空间技术的研究热点。空间在轨服务技术具有广泛的应用前景,可以对空间故障卫星在轨维修、部件升级,可以对未成功入轨的卫星进行辅助入轨、任务恢复,也可以用于空间碎片清除与空间攻防对抗等领域。鉴于其对在空间态势及国民经济等方面的重大意义,世界各国纷纷展开了针对在轨服务技术的相关研究,包括轨道快车技术、工程试验卫星-VII以及试验卫星系统在内的多项任务都是对在轨服务相关技术的试验与验证。
在进行空间在轨操作时,通常需要经过以下几个阶段:自主安全交会,抓捕对接操作以及抓捕后稳定操作。其中,自主安全交会作为后续其他操作的前提,直接关系到了整个在轨服务的成败。现有的针对合作目标的自主交会对接技术已经十分的成熟,并且被广泛地应用于空间站补给、火星采样返回等任务中。然而,对于非合作目标的自主交会技术从未进行过相关的在轨演示验证,非合作目标自主交会技术仍是当前国内外研究的热点。非合作目标是指那些不能通过星间链路传递自身信息,没有合作标识器的一类空间物体。非合作目标在自身残余角动量以及空间复杂摄动影响下通常会处于翻滚状态,同时非合作目标的运动信息不能被服务航天器完全准确地获得,这使得对于非合作的自主安全交会变得十分的困难。本发明对非合作目标的自主交会问题进行了研究,采用基于零空间方法的轨迹规划算法,设计出一条安全的自主交会轨迹。零空间方法是一种能够同时处理多种任务要求的方法,它根据任务的紧迫性不同,对不同任务进行优先级排序,将低优先级的任务向高优先级任务的零空间投影,从而保证低优先级的任务不会影响到高优先级任务的执行。只有当高优先级任务完成后才会完成低优先级任务,避免了任务之间的相互冲突。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足,提供一种基于零空间方法的非合作目标自主交会策略,对交会过程中的碰撞规避要求和对接点跟踪要求进行优先级排序,从而规划出一条安全的交会轨迹。
本发明采用以下技术方案:
一种基于零空间方法的非合作目标自主交会策略,首先对自主交会过程进行任务分析,将整个交会过程为非合作目标对接抓捕点跟踪任务和碰撞规避任务,对所述非合作目标对接抓捕点跟踪任务和碰撞规避任务进行优先级排序,并对每种基本任务所对应的零空间进行求解,最后再根据零空间方法对所述非合作目标对接抓捕点跟踪任务和碰撞规避任务进行综合,将低优先级任务投影到高优先级任务的零空间,从而保证了交会轨迹的安全性。
优选的,具体包括以下步骤:
S1、对空间非合作目标的翻滚运动进行建模,确定翻滚过程中抓捕对接点的运动规律;
S2、对步骤S1所述非合作目标自主交会任务进行分解,并分别对所有基本任务进行分析描述;
S3、对所述基本任务进行优先级排序,同时对其零空间进行求解;
S4、对步骤S3所述基本任务的速度输出进行求解,并对整个自主交会任务的速度输出进行综合,得到整个任务的速度输出;
S5、基于零空间方法生成自主交会轨迹,并进行仿真分析与验证。
优选的,步骤S1中,所述非合作目标的翻滚运动可以通过姿态运动方程描述如下:
其中,Jt表示非合作目标的惯量张量,ωt表示非合作目标本体系相对于惯性坐标系的角速度矢量,表示非合作目标本体坐标系相对于惯性坐标系的姿态四元数,τ表示外界扰动力矩。
优选的,所述非合作目标对接抓捕点在目标轨道坐标系中的表示为:
其中,rd表示对接抓捕点在目标本体坐标系中的分量,ρd表示对接抓捕点在目标轨道坐标系中的分量,表示由惯性坐标系到目标轨道系的方向余弦矩阵,表示由惯性坐标系到目标本体坐标系的方向余弦矩阵,表示对方向余弦矩阵的转置。
优选的,所述非合作目标对接抓捕点在轨道坐标系中的分量可以表示为:
其中,vd表示非合作目标抓捕点在轨道坐标系中的矢量表示,n表示非合作目标的轨道角速率。
优选的,步骤S2中,所述非合作目标自主交会任务包括碰撞规避任务和非合作目标对接抓捕点跟踪任务;
所述碰撞规避任务为服务航天器和非合作目标之间的距离:
m1=min{||ρ-ρo||-d,0}
其中,m1为碰撞规避任务中的可控变量,ρ表示当前服务航天器的位置矢量,ρo表示非合作目标的位置矢量,d表示交会过程中所设定的服务航天器和非合作目标之间的安全距离;
所述对接抓捕点跟踪任务为:
m2=ρ
其中,m2表示对对接抓捕点跟踪任务中的可控变量。
优选的,步骤S3中,所述碰撞规避任务具有最高的优先级,所述对接抓捕点跟踪任务为最低优先级。
优选的,所述可控变量m1的零空间为:
其中,表示可控变量m1所对应的零空间,J1(ρ)表示对应于可控变量m1的雅可比矩阵,表示矩阵J1(ρ)的广义拟矩阵;
所述可控变量m2的零空间为:
其中,03表示3×3的零矩阵,表示可控变量m2所对应的零空间,J2(ρ)表示对应于可控变量m2的雅可比矩阵,表示矩阵J2(ρ)的广义拟矩阵。
优选的,步骤S4中,根据零空间方法,对所述碰撞规避任务和对接抓捕点跟踪任务的输出v进行综合如下:
其中,v1为碰撞规避任务的速度输出,v2为对接抓捕点任务的速度输出,v表示任务综合后的最终速度输出。
优选的,步骤S5中,对积分求解并进行离散化,得到如下:
ρsafe(tk)=ρsafe(tk-1)+v(tk)Δt
其中,tk表示第k个采样时间,ρsafe(tk)表示tk时刻由零空间方法得到的安全位置矢量,Δt表示采样周期。
与现有技术相比,本发明至少具有以下有益效果:
本发明首先对自主交会过程进行任务分析,将整个交会过程为两个基本任务,分别是非合作目标对接抓捕点跟踪任务和碰撞规避任务。在此基础上,本发明对上述两种任务进行优先级排序,并对每种基本任务所对应的零空间进行求解,最后再根据零空间方法对这两种任务进行综合,将低优先级任务投影到高优先级任务的零空间,从而保证了交会轨迹的安全性。
进一步的,本发明先对空间非合作目标的翻滚运动进行建模,确定翻滚过程中抓捕对接点的运动规律,再对非合作目标自主交会任务分解,并对所有基本任务分别进行分析描述,然后对上述基本任务进行优先级排序,同时对其零空间进行求解,再对上述基本任务的速度输出进行求解,并对整个自主交会任务的速度输出进行综合,最后基于零空间方法生成自主交会轨迹,并进行仿真分析与验证,将一个复杂的任务分解成几个基本的任务单元,然后对基本任务进行优先级排序。通过将低优先级的任务投影到高优先级任务的零空间,从而保证高优先级任务得到优先处理,避免任务之间的冲突,基于零空间方法的非合作目标自主交会策略,可以使得服务航天器在逼近非合作目标的过程中避免发生碰撞,从而确保交会任务的顺利完成。
本发明所提供的基于零空间方法的自主交会策略是一种分层递阶控制的方法,它能够确保高优先级任务最先得到完成,从而避免不同优先级任务之间发生冲突。基于零空间方法的非合作目标自主交会策略能够同时实现对非合作目标对接抓捕点的跟踪以及碰撞规避要求。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明的流程图;
图2为本发明自主交会过程中相对距离的变化曲线示意图;
图3为本发明自主交会过程中位置跟踪误差的变化曲线示意图;
图4为本发明自主交会过程中跟踪速度的变化曲线示意图。
具体实施方式
本发明提供了一种基于零空间方法的非合作目标自主交会策略,对空间非合作目标自主交会过程中的安全交会轨迹进行设计,确保服务航天器在逼近非合作目标的过程避免与其发生碰撞。
请参阅图1,本发明公开了一种基于零空间方法的非合作目标自主交会策略,首先对自主交会过程进行任务分析,将整个交会过程为非合作目标对接抓捕点跟踪任务和碰撞规避任务。对所述非合作目标对接抓捕点跟踪任务和碰撞规避任务进行优先级排序,并对每种基本任务所对应的零空间进行求解,最后再根据零空间方法对所述非合作目标对接抓捕点跟踪任务和碰撞规避任务进行综合,将低优先级任务投影到高优先级任务的零空间,从而保证了交会轨迹的安全性。
具体步骤如下:
S1、对空间非合作目标的翻滚运动进行建模,确定翻滚过程中抓捕对接点的运动规律;
非合作目标在自身残余角动量及空间复杂摄动环境影响下通常处于自由翻滚状态,因此非合作目标的翻滚运动可以通过姿态运动方程描述如下:
其中,Jt表示非合作目标的惯量张量,ωt表示非合作目标本体系相对于惯性坐标系的角速度矢量,表示非合作目标本体坐标系相对于惯性坐标系的姿态四元数,τ表示外界扰动力矩。对于近距离自主交会任务,通常任务的持续时间很短,因此外界扰动力矩可以忽略即τ=0。
由于非合作目标处于自由翻滚状态,因此导致固连在非合作目标本体上的对接抓捕点也是随时间发生动态变化。为了确保下一步抓捕对接操作的顺利进行,必须要求服务航天器在自主交会的末端时刻完成对对接抓捕点的精确跟踪,因此需要将对接抓捕点的运动转换到目标轨道坐标系中,转化过程如下:
由非合作目标的轨道六要素可以方便地得到由惯性坐标系到目标轨道坐标系的方向余弦矩阵如下:
其中,表示由惯性坐标系到目标轨道系的方向余弦矩阵,Ω,i,θ分别表示非合作目标的升交点赤经,轨道倾角以及纬度幅角。
同时,由惯性坐标系到目标本体坐标系的方向余弦矩阵可以表示为
其中,表示由惯性坐标系到目标本体坐标系的方向余弦矩阵。
因此,非合作目标对接抓捕点在目标轨道坐标系中的表示为
其中,rd表示对接抓捕点在目标本体坐标系中的分量,ρd表示对接抓捕点在目标轨道坐标系中的分量,表示由惯性坐标系到目标轨道系的方向余弦矩阵,表示由惯性坐标系到目标本体坐标系的方向余弦矩阵,表示对方向余弦矩阵的转置。
此外,非合作目标对接抓捕点在轨道坐标系中的分量可以表示为
其中,vd表示非合作目标抓捕点在轨道坐标系中的矢量表示,n表示非合作目标的轨道角速率。
S2、对非合作目标自主交会任务分解,并对所有基本任务分别进行分析描述;
对于整个自主交会任务,可以将其分为两个基本任务,分别是碰撞规避任务和非合作目标对接抓捕点跟踪任务。这两个基本任务具体描述如下:
在近距离交会过程中,碰撞规避任务能够保证服务航天器的安全和完整,为后续任务的实现提供基本的保障。通常避障任务可以表述为服务航天器和非合作目标之间的距离
m1=min{||ρ-ρo||-d,0}
其中,m1为碰撞规避任务中的可控变量,ρ表示当前服务航天器的位置矢量,ρo表示非合作目标的位置矢量,d表示交会过程中所设定的服务航天器和非合作目标之间的安全距离。
除了避障任务外,服务航天器还需要跟踪非合作目标的对接抓捕点,以便进行下一步的抓捕操作。根据上述描述,对接抓捕点跟踪任务可以方便地描述为
m2=ρ
其中,m2表示对接点跟踪任务中的可控变量。
S3、对上述基本任务进行优先级排序,同时对其零空间进行求解;
根据步骤S2可知,自主交会过程中涉及到两个基本任务,分别是碰撞规避任务和对接抓捕点跟踪任务。从这两个任务中可以轻易地分析出其优先性,碰撞规避任务是整个后续任务进行的前提,故而其具有最高的优先级,所以对接抓捕点跟踪任务的优先级最低。
由于零空间方法需要将低优先级的任务投影到高优先级任务的零空间,因此就需要对两种基本任务的任务零空间进行求解。根据零空间方法,需要对可控变量m1和m2分别求雅可比矩阵,具体求解如下:
对可控变量m1求雅可比矩阵可得
其中,J1(ρ)表示对应于可控变量m1的雅可比矩阵。由J1(ρ)可以得到其逆矩阵表示为
其中,表示矩阵J1(ρ)的广义拟矩阵。
因此根据矩阵论可知,可控变量m1的零空间为
其中,表示可控变量m1所对应的零空间,J1(ρ)表示对应于可控变量m1的雅可比矩阵,表示矩阵J1(ρ)的广义拟矩阵。
对可控变量m2求雅可比矩阵可得
其中,J2(ρ)表示对应于可控变量m2的雅可比矩阵,I3表示3×3的单位矩阵。于是可以方便地得到矩阵J2(ρ)的逆矩阵
于是可控变量m2的零空间为
其中,03表示3×3的零矩阵,表示可控变量m2所对应的零空间,J2(ρ)表示对应于可控变量m2的雅可比矩阵,表示矩阵J2(ρ)的广义拟矩阵。
S4、对上述基本任务的速度输出进行求解,并对整个自主交会任务的速度输出进行综合;
根据零空间方法,可以得到碰撞规避任务的速度输出为:
其中,λ1为任务增益,m1d表示碰撞规避任务的期望值,这里可以将其设置为m1d=0。
同理,可以得到对接抓捕点任务的速度输出为
其中,λ2为任务增益,m2d表示碰撞规避任务的期望值,这里可以将其设置为m1d=ρd,v2d表示对接抓捕点跟踪任务的期望速度,这里设置为v2d=vd
根据零空间方法,对上述两个任务的输出进行综合,可得整个任务的速度输出
其中,v1为碰撞规避任务的速度输出,v2为对接抓捕点任务的速度输出,v表示任务综合后的最终速度输出
S5、基于零空间方法生成自主交会轨迹,并进行仿真分析与验证。
由于零空间方法是一种速度级的方法,因此要想获得最终的安全交会轨迹,需要进行积分求解,即
在数值求解过程中必须对上述积分进行离散化,于是可得
ρsafe(tk)=ρsafe(tk-1)+v(tk)Δt
其中tk表示第k个采样时间,ρsafe(tk)表示tk时刻由零空间方法得到的安全位置矢量,Δt表示采样周期。
仿真实例,为了说明本方法的可行性及有效性,假设非合作目标的惯量矩阵为
Jt=diag(10,15,10)kg/m2
其中,初始角速度为ωt=[0.4,0.5,0.8]Trad/s,初始的姿态四元数为q(0)=[1,0,0,0]T,假定非合作目标位于圆轨道上轨道半径为r=7178160(m),轨道倾角i、经地点幅角Ω、以及纬度幅角都取为0rad。同时服务航天器相对于目标的初始位置设置为ρ(0)=[15;8;5](m),交会过程中的安全距离d=2(m),对抓捕点在目标本体坐标系中的位置矢量为rd=[0,0,-2]T(m),任务增益λ1和λ2都取为0.8。
基于零空间方法的非合作目标自主交会仿真结果如图2-4所示。图2表示在接近非合作目标的过程中相对距离随时间的变化曲线,图3表示交会过程中位置跟踪误差随时间的变化曲线,图4表示交会过程中跟踪速度随时间的变化曲线。
从图2可以发现服务航天器在接近非合作目标的过程中能够很好的避免和非合作目标发生碰撞,从而确保碰撞规避任务的完成。图3表明服务航天器能够有效地跟踪非合作目标抓捕点的运动,进而确保对接点跟踪任务的顺利完成。图4表示服务航天器在跟踪非合作目标对接点运动过程中所需要的跟踪速度,从图中可以发现跟踪速度连续有界,因此很容易实现。
由图2-4可以发现,基于零空间方法的非合作目标自主交会策略能够保证碰撞规避和对接点跟踪任务同时完成,并且避免两种任务之间发生冲突,此外优先级最高的碰撞规避任务最先得到保证,这样可以提高服务航天器在接近过程中的安全性,同时又能在此前提下实现最终的安全交会。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。

Claims (10)

1.一种基于零空间方法的非合作目标自主交会策略,其特征在于,首先对自主交会过程进行任务分析,将整个交会过程为非合作目标对接抓捕点跟踪任务和碰撞规避任务,对所述非合作目标对接抓捕点跟踪任务和碰撞规避任务进行优先级排序,并对每种基本任务所对应的零空间进行求解,最后再根据零空间方法对所述非合作目标对接抓捕点跟踪任务和碰撞规避任务进行综合,将低优先级任务投影到高优先级任务的零空间,从而保证了交会轨迹的安全性。
2.根据权利要求1所述的一种基于零空间方法的非合作目标自主交会策略,其特征在于,具体包括以下步骤:
S1、对空间非合作目标的翻滚运动进行建模,确定翻滚过程中抓捕对接点的运动规律;
S2、对步骤S1所述非合作目标自主交会任务进行分解,并分别对所有基本任务进行分析描述;
S3、对所述基本任务进行优先级排序,同时对其零空间进行求解;
S4、对步骤S3所述基本任务的速度输出进行求解,并对整个自主交会任务的速度输出进行综合,得到整个任务的速度输出;
S5、基于零空间方法生成自主交会轨迹,并进行仿真分析与验证。
3.根据权利要求2所述的一种基于零空间方法的非合作目标自主交会策略,其特征在于,步骤S1中,所述非合作目标的翻滚运动可以通过姿态运动方程描述如下:
<mrow> <msub> <mi>J</mi> <mi>t</mi> </msub> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>t</mi> </msub> <mo>+</mo> <msup> <msub> <mi>&amp;omega;</mi> <mi>t</mi> </msub> <mo>&amp;times;</mo> </msup> <msub> <mi>J</mi> <mi>t</mi> </msub> <msub> <mi>&amp;omega;</mi> <mi>t</mi> </msub> <mo>=</mo> <mi>&amp;tau;</mi> </mrow>
<mrow> <msub> <mover> <mi>q</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>v</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <msup> <msub> <mi>q</mi> <mi>v</mi> </msub> <mo>&amp;times;</mo> </msup> <msub> <mi>&amp;omega;</mi> <mi>t</mi> </msub> <mo>+</mo> <msub> <mi>q</mi> <mn>0</mn> </msub> <msub> <mi>&amp;omega;</mi> <mi>t</mi> </msub> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mover> <mi>q</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <msub> <mi>q</mi> <mi>v</mi> </msub> <mi>T</mi> </msup> <msub> <mi>&amp;omega;</mi> <mi>t</mi> </msub> </mrow>
其中,Jt表示非合作目标的惯量张量,ωt表示非合作目标本体系相对于惯性坐标系的角速度矢量,表示非合作目标本体坐标系相对于惯性坐标系的姿态四元数,τ表示外界扰动力矩。
4.根据权利要求3所述的一种基于零空间方法的非合作目标自主交会策略,其特征在于,所述非合作目标对接抓捕点在目标轨道坐标系中的表示为:
<mrow> <msub> <mi>&amp;rho;</mi> <mi>d</mi> </msub> <mo>=</mo> <msubsup> <mi>C</mi> <mi>i</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> </msubsup> <msup> <mrow> <mo>(</mo> <msubsup> <mi>C</mi> <mi>i</mi> <mrow> <mi>t</mi> <mi>b</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mi>T</mi> </msup> <msub> <mi>r</mi> <mi>d</mi> </msub> </mrow>
其中,rd表示对接抓捕点在目标本体坐标系中的分量,表示由惯性坐标系到目标轨道系的方向余弦矩阵,表示由惯性坐标系到目标本体坐标系的方向余弦矩阵,ρd表示对接抓捕点在目标轨道坐标系中的分量,表示对方向余弦矩阵的转置。
5.根据权利要求4所述的一种基于零空间方法的非合作目标自主交会策略,其特征在于,所述非合作目标对接抓捕点在轨道坐标系中的分量可以表示为:
<mrow> <msub> <mi>v</mi> <mi>d</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msubsup> <mi>C</mi> <mi>i</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> </msubsup> <msup> <mrow> <mo>(</mo> <msubsup> <mi>C</mi> <mi>i</mi> <mrow> <mi>t</mi> <mi>b</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mi>T</mi> </msup> <msub> <mi>&amp;omega;</mi> <mi>t</mi> </msub> <mo>-</mo> <mo>&amp;lsqb;</mo> <mn>0</mn> <mo>;</mo> <mn>0</mn> <mo>;</mo> <mi>n</mi> <mo>&amp;rsqb;</mo> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <msub> <mi>&amp;rho;</mi> <mi>d</mi> </msub> </mrow>
其中,vd表示非合作目标抓捕点在轨道坐标系中的矢量表示,n表示非合作目标的轨道角速率。
6.根据权利要求2所述的一种基于零空间方法的非合作目标自主交会策略,其特征在于,步骤S2中,所述非合作目标自主交会任务包括碰撞规避任务和非合作目标对接抓捕点跟踪任务;
所述碰撞规避任务为服务航天器和非合作目标之间的距离:
m1=min{||ρ-ρo||-d,0}
其中,m1为碰撞规避任务中的可控变量,ρ表示当前服务航天器的位置矢量,ρo表示非合作目标的位置矢量,d表示交会过程中所设定的服务航天器和非合作目标之间的安全距离;
所述对接抓捕点跟踪任务为:
m2=ρ
其中,m2表示对对接抓捕点跟踪任务中的可控变量。
7.根据权利要求6所述的一种基于零空间方法的非合作目标自主交会策略,其特征在于,步骤S3中,所述碰撞规避任务具有最高的优先级,所述对接抓捕点跟踪任务为最低优先级。
8.根据权利要求7所述的一种基于零空间方法的非合作目标自主交会策略,其特征在于,所述可控变量m1的零空间为:
其中,表示可控变量m1所对应的零空间,J1(ρ)表示对应于可控变量m1的雅可比矩阵,表示矩阵J1(ρ)的广义拟矩阵;
所述可控变量m2的零空间为:
其中,03表示3×3的零矩阵,表示可控变量m2所对应的零空间,J2(ρ)表示对应于可控变量m2的雅可比矩阵,表示矩阵J2(ρ)的广义拟矩阵。
9.根据权利要求8所述的一种基于零空间方法的非合作目标自主交会策略,其特征在于,步骤S4中,根据零空间方法,对所述碰撞规避任务和对接抓捕点跟踪任务的输出v进行综合如下:
<mrow> <mi>v</mi> <mo>=</mo> <msub> <mi>v</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mover> <mi>m</mi> <mo>&amp;OverBar;</mo> </mover> <mn>1</mn> </msub> <msub> <mi>v</mi> <mn>2</mn> </msub> </mrow>
其中,v1为碰撞规避任务的速度输出,v2为对接抓捕点任务的速度输出,v表示任务综合后的最终速度输出。
10.根据权利要求2所述的一种基于零空间方法的非合作目标自主交会策略,其特征在于,步骤S5中,对积分求解并进行离散化,得到如下:
ρsafe(tk)=ρsafe(tk-1)+v(tk)Δt
其中,tk表示第k个采样时间,ρsafe(tk)表示tk时刻由零空间方法得到的安全位置矢量,Δt表示采样周期。
CN201710225648.2A 2017-04-07 2017-04-07 一种基于零空间方法的非合作目标自主交会策略 Active CN107065547B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710225648.2A CN107065547B (zh) 2017-04-07 2017-04-07 一种基于零空间方法的非合作目标自主交会策略

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710225648.2A CN107065547B (zh) 2017-04-07 2017-04-07 一种基于零空间方法的非合作目标自主交会策略

Publications (2)

Publication Number Publication Date
CN107065547A true CN107065547A (zh) 2017-08-18
CN107065547B CN107065547B (zh) 2019-12-03

Family

ID=59602302

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710225648.2A Active CN107065547B (zh) 2017-04-07 2017-04-07 一种基于零空间方法的非合作目标自主交会策略

Country Status (1)

Country Link
CN (1) CN107065547B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109062044A (zh) * 2018-08-03 2018-12-21 北京航空航天大学 一种终端迭代学习对接控制方法
CN111882184A (zh) * 2020-07-14 2020-11-03 福州大学 一种多智能体系统零空间行为控制动态任务优先级规划方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070228219A1 (en) * 2006-03-31 2007-10-04 Behrens John W Two part spacecraft servicing vehicle system with universal docking adaptor
CN101423121A (zh) * 2008-12-08 2009-05-06 北京航空航天大学 应用双视线测量的非合作目标航天器相对导航制导方法
CN102236825A (zh) * 2010-04-26 2011-11-09 兰州交通大学 一种机场货运物流系统的输送机任务调度方法
CN102692572A (zh) * 2012-05-14 2012-09-26 北京航空航天大学 一种基于时间可用度的飞机射频设备间电磁兼容性分析方法
CN103438888A (zh) * 2013-07-24 2013-12-11 西北工业大学 一种对空间非合作目标自主交会的相对导航方法
CN104316060A (zh) * 2014-06-06 2015-01-28 清华大学深圳研究生院 空间非合作目标的交会对接方法与装置
CN104407611A (zh) * 2014-09-30 2015-03-11 同济大学 一种仿人机器人稳定行走的控制方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070228219A1 (en) * 2006-03-31 2007-10-04 Behrens John W Two part spacecraft servicing vehicle system with universal docking adaptor
CN101423121A (zh) * 2008-12-08 2009-05-06 北京航空航天大学 应用双视线测量的非合作目标航天器相对导航制导方法
CN102236825A (zh) * 2010-04-26 2011-11-09 兰州交通大学 一种机场货运物流系统的输送机任务调度方法
CN102692572A (zh) * 2012-05-14 2012-09-26 北京航空航天大学 一种基于时间可用度的飞机射频设备间电磁兼容性分析方法
CN103438888A (zh) * 2013-07-24 2013-12-11 西北工业大学 一种对空间非合作目标自主交会的相对导航方法
CN104316060A (zh) * 2014-06-06 2015-01-28 清华大学深圳研究生院 空间非合作目标的交会对接方法与装置
CN104407611A (zh) * 2014-09-30 2015-03-11 同济大学 一种仿人机器人稳定行走的控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张业鑫: "基于TS模糊理论的微小型航天器群的协同控制研究", 《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》 *
高登巍: "接近和跟踪非合作机动目标的非线性最优控制", 《宇航学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109062044A (zh) * 2018-08-03 2018-12-21 北京航空航天大学 一种终端迭代学习对接控制方法
CN109062044B (zh) * 2018-08-03 2020-10-09 北京航空航天大学 一种终端迭代学习对接控制方法
CN111882184A (zh) * 2020-07-14 2020-11-03 福州大学 一种多智能体系统零空间行为控制动态任务优先级规划方法
CN111882184B (zh) * 2020-07-14 2022-10-14 福州大学 多智能体系统零空间行为控制动态任务优先级规划方法

Also Published As

Publication number Publication date
CN107065547B (zh) 2019-12-03

Similar Documents

Publication Publication Date Title
CN104898691B (zh) 编队卫星有限时间构型包含控制方法
Breger et al. Safe trajectories for autonomous rendezvous of spacecraft
Izzo et al. Autonomous and distributed motion planning for satellite swarm
CN106707751A (zh) 航天器终端接近的有限时间饱和避碰控制方法
CN106945020A (zh) 一种空间双机械臂系统运动协调控制方法
Meng et al. Space robots with flexible appendages: dynamic modeling, coupling measurement, and vibration suppression
Li et al. Disturbance observer based control for spacecraft proximity operations with path constraint
CN105138010A (zh) 一种编队卫星分布式有限时间跟踪控制方法
Wilde et al. Experimental characterization of inverse dynamics guidance in docking with a rotating target
CN105159304A (zh) 接近并跟踪空间非合作目标的有限时间容错控制方法
CN103728980A (zh) 航天器相对轨道的控制方法
Chan et al. Verifying safety of an autonomous spacecraft rendezvous mission
Chu et al. Hybrid tension control method for tethered satellite systems during large tumbling space debris removal
CN104656447A (zh) 一种航天器抗干扰姿态跟踪的微分几何非线性控制方法
CN106054613A (zh) 一种自主安全逼近翻滚故障卫星的轨迹规划方法
Morgan et al. Spacecraft swarm guidance using a sequence of decentralized convex optimizations
Yuan et al. Probability-based hazard avoidance guidance for planetary landing
Liu et al. Dynamic boundary layer based neural network quasi-sliding mode control for soft touching down on asteroid
CN107065547B (zh) 一种基于零空间方法的非合作目标自主交会策略
Gong et al. Barrier Lyapunov function-based planetary landing guidance for hazardous terrains
Chen et al. Integrated guidance and control for microsatellite real-time automated proximity operations
Bucchioni et al. Open loop safe trajectory design for cislunar NRHO rendezvous
Zhan et al. A novel collision-free detumbling strategy for a space robot with a 7-DOF manipulator in postcapturing phase
CN105678061A (zh) 一种满足被动安全要求的空间站安全区域设计方法
CN111596678B (zh) 一种航天器逼近空间翻滚目标的相对位置自适应控制方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant