CN104143021B - 一种空间绳网机器人的逼近动力学建模方法 - Google Patents

Abstract

本发明公开了一种空间绳网机器人的逼近动力学建模方法，针对空间绳网机器人的逼近动力学问题，研究其动力学建模及解算方法，首先建立建模参考坐标系并提出建模假设条件，再建立柔性网模型以及空间绳网机器人的逼近动力学模型，最后进行系绳运动的速度跳变建模；由于空间绳网机器人在逼近目标过程中，四个单片网机构均处于绷紧状态，以平面有限元理论中的T3单元为基础，利用位置矢量的矩阵形式可以有效描述单片网机构的动力学。本发明可以有效降低单片网结构动力学模型复杂度，提高模型解算速度。

Description

一种空间绳网机器人的逼近动力学建模方法

技术领域

本发明属于新型航天器动力学建模研究的领域，具体涉及一种空间绳网机器人的逼近动力学建模方法。

背景技术

空间绳网机器人是一种新型的“空间平台+连接系绳+柔性网(含自主机动单元)”结构的空间机器人系统，可针对各种目标进行捕获，且具有极大的容差能力，主要应用于空间垃圾清理任务。其主要工作流程为：首先由搭载空间绳网机器人的空间平台靠近目标至空间绳网机器人的作用距离范围内；然后由空间平台发射空间绳网机器人的柔性网，再由自主机动单元在自身的位姿控制机构和导航设备以及空间平台测控系统的支持下，控制柔性网按预定的轨迹与速度向目标进行逼近飞行；最后由柔性网碰撞并包裹目标，从而完成对于目标的捕获，并由空间平台利用连接系绳将目标拖曳至大气层或坟墓轨道，然后切断连接系绳。空间绳网机器人的结构如附图1所示。1表示空间平台，2表示连接系绳，3表示柔性网，4为连接系绳与柔性网的连接点，5,6,7,8表示安装于柔性网上的四个自主机动单元。

但是，柔性网使得空间绳网机器人的逼近动力学变得十分复杂，动力学建模非常困难。文献(于洋，宝音贺西，李俊锋.空间柔性网抛射展开过程动力学建模与仿真[J].宇航学报，2010，31(5)：1289-1296)和文献(敬忠良，袁建平等.航天器自主操作的测量与控制[M].北京：中国宇航出版社，2011：493-518)将柔性网离散化为质点和弹性杆单元的结构，在此基础上建立了柔性网的动力学模型；文献(Provot X.Deformation constraints in amass-spring model to describe rigid cloth behavior[C].Proceedings of GraphicsInterface，Quebec，Canada，1995)采用质点弹簧模型建立了网结构的动力学模型。但是这类质点弹簧模型虽然建模原理简单，但建立的模型十分复杂，且以各质点为对象的数学模型很难用于控制系统设计。另外，由于组成柔性网的系绳杨氏模量极高，模型仿真计算缓慢。研究设计一种空间绳网机器人的新型建模方法势在必行。

发明内容

本发明的目的是针对空间绳网机器人的逼近动力学问题，研究其动力学建模及解算方法，提供一种空间绳网机器人的逼近动力学建模方法，该方法为空间绳网机器人的研究奠定了基础。

为了实现上述目的，本发明所采用的技术方案包括以下步骤：

1)建立建模参考坐标系并提出建模假设条件；

2)建立单片柔性网的模型；

3)建立空间绳网机器人的逼近动力学模型；

4)建立系绳连接点与自主机动单元间系绳运动的速度跳变模型。

所述的步骤1)中，建立建模参考坐标系的具体方法是：

首先建立地心惯性系OXYZ，其坐标原点O位于地球中心，OX轴指向地球的春分点，OZ轴指向地球北极，OY轴在赤道平面内垂直于OX轴；然后建立轨道坐标系O_oX_oY_oZ_o，其坐标原点O_o位于空间平台的质心，O_oX_o轴沿机动空间平台轨道的切向，指向空间平台运动方向，O_oZ_o轴与原点O_o和地心O的连线重合，指向地心，O_oY_o轴沿轨道平面负法线方向。

所述的步骤1)中，建模假设条件具体如下：

假设1.空间平台运行于圆轨道，且质量要远大于系绳、柔性网和自主机动单元质量的总和；

假设2.忽略空间平台与柔性网的连接系绳的弹性与质量，忽略自主机动单元体积，逼近任务中，连接系绳处于绷紧状态；

假设3.柔性网的质量分布均匀，网孔很小，在逼近目标过程中不大幅变形；

假设4.由于系绳杨氏模量极大，假设系绳连接点与自主机动单元间、自主机动单元间的系绳不可伸长。

所述的步骤2)中，建立单片柔性网模型的具体方法是：

系绳连接点与柔性网的四个自主机动单元将柔性网分为四块，利用假设3，将四片柔性网分别建模为三角形薄壳；以系绳连接点与柔性网的第一自主机动单元和第二自主机动单元组成的单片柔性网为例说明，A表示系绳连接点，B表示第一自主机动单元，C表示第二自主机动单元，R₁，R₂，R₃分别表示A,B,C在地心惯性系下的位置矢量；

针对单片柔性网，采用平面有限元理论中的T3单元来描述；对于单片柔性网上任一点D，其在地心惯性系下的位置矢量为：

R≈s₁R₁+s₂R₂+s₃R₃ (1)

式中，s₁、s₂、s₃表示点D在薄壳上的面积坐标，它们满足：

其中，表示单片柔性网在未发生任何变形条件下三角形的面积；

建立系绳连接点与第一自主机动单元和第二自主机动单元组成的单片柔性网的拉格朗日函数：

式中，dΣ表示三角形薄壳上的面积微元，ρ表示单片柔性网的平均面密度，m_W表示整个柔性网的总质量，G表示万有引力常数，M表示地球的质量；

对单片柔性网的拉格朗日函数求变分，得：

其中，δ为变分符号，t₁,t₂表示积分时间；

利用C-W方程将上式转换到轨道坐标系下，得到：

其中，r为单片柔性网上点D在轨道坐标系O_oX_oY_oZ_o下的位置矢量；ω为空间平台轨道运动的平均角速度；

同理，写出其他三片柔性网的表达式；设r₄,r₅,r₆,r₇,r₈分别表示系绳连接点、第一自主机动单元、第二自主机动单元、第三自主机动单元以及第四自主机动单元在轨道坐标系O_oX_oY_oZ_o下的位置矢量，将其写为矩阵形式：

则：

其中，L₄₆₇L₄₇₈L₄₈₅分别表示系绳连接点与第二自主机动单元、第三自主机动单元，系绳连接点与第三自主机动单元、第四自主机动单元，系绳连接点与第四自主机动单元、第一自主机动单元分别组成的三个单片柔性网的拉格朗日函数；表示矩阵的直积运算；

(lmn＝456,467,478,485)，M^lmn为5×5矩阵，其任意元素满足：

所述的步骤3)中，建立空间绳网机器人的逼近动力学模型的具体方法是：

首先分析系绳连接点与四个自主机动单元间、四个自主机动单元间系绳的不可伸长假设对逼近动力学建模带来的影响；

以系绳连接点和第一自主机动单元间系绳为例，设L₄₅为连接点和第一自主机动单元之间系绳的标称长度，则，||r₄-r₅||≤L₄₅；引入间隙函数g₄₅对这一现象进行描述：

g₄₅＝L₄₅-||r₄-r₅||≥0 (7)

对应的约束反力λ₄₅满足：

后两种形式的约束反力与系统状态的二次导数相关，将其写为隐式方程：

同理，写出系绳连接点与第二自主机动单元，连接点与第三自主机动单元，连接点与第四自主机动单元，第一自主机动单元与第二自主机动单元，第二自主机动单元与第三自主机动单元，第三自主机动单元与第四自主机动单元，第四自主机动单元与第一自主机动单元之间系绳的间隙函数(g₄₆,g₄₇,g₄₈,g₅₆,g₆₇,g₇₈,g₈₅)和约束反力λ₄₆,λ₄₇,λ₄₈,λ₅₆,λ₆₇,λ₇₈,λ₈₅表达；

将间隙函数和约束反力写为矩阵，得：

g_N＝[g₄₅,g₄₆,g₄₇,g₄₈,g₅₆,g₆₇,g₇₈,g₈₅]^T

λ_N＝[λ₄₅,λ₄₆,λ₄₇,λ₄₈,λ₅₆,λ₆₇,λ₇₈,λ₈₅]^T

利用拉格朗日法对空间绳网机器人的逼近动力学进行建模

整个系统的拉格朗日函数满足：

式中，m_P表示空间平台的质量，m_M表示自主机动单元的质量；R_o表示空间平台在地心惯性系下的坐标，R_i(i＝5,…,8)表示第一自主机动单元、第二自主机动单元、第三自主机动单元以及第四自主机动单元在地心惯性系下的坐标；

另外，系统中非保守力所做的功满足：

式中，r₄表示系绳连接点在轨道坐标系O_oX_oY_oZ_o下的位置矢量，F_T表示空间平台和柔性网的连接系绳拉力，F_i(i＝5,…,8)表示作用在自主机动单元上的推力；

由广义Hamilton原理可知，空间绳网机器人的逼近动力学满足：

将式(10)和式(11)代入方程(12)并使用分步积分进行化简可得：

式中，

对(13)式进行整理，同时考虑变分的任意性得系统的动力学方程为：

式中，

整个系统的解算通过构建隐式方程，利用Lemke算法进行求解。

所述的步骤4)中，系绳运动的速度跳变建模的具体方法是：

假设在t^-时刻，节点i和节点j之间的距离达到标称长度L_ij，而两个节点的相对速度不为0，且有相互远离的趋势；于是，在t时刻连接两个节点的系绳绷紧，由于系绳的刚度和阻尼都比较高，所以类似于碰撞过程，系绳会在极短时间内产生了比较大的张紧力，从而使得两个节点在t⁺时刻具有了沿系绳方向一致的速度；

为了描述瞬间的速度跳变，将动力学方程(14)在[t^-,t⁺]上积分，得：

由于在极短的时间内，位置矢量r_N还来不及发生变化，同时由于控制力F是有限的，因此方程左侧的后两项和右侧的第一项为0；对于那些满足系绳绷直条件的约束，即满足g_i＝0且其约束力在无限小的时间内趋向于无穷大，因此其积分不为0，而对于那些不满足绷直条件的约束，其约束反力为有限值，因而积分也为0；

为了便于表达，将满足绷直条件的约束取出，将其编号从小到大排列得到列向量按照同样的顺序排列对应的间隙函数和约束反力能够得到列向量和另外，分别将约束反力向量λ_NT和在时间段[t^-,t⁺]记作Λ_NT和两者之间满足：

其中，

于是，方程(15)能够化简为：

另外，对于间隙函数向量它在t⁺时刻必须满足：

则：绷紧系绳的速度跳变模型为：

与现有技术相比，本发明具有以下有益效果：

本发明针对空间绳网机器人的逼近动力学问题，研究其动力学建模及解算方法，首先建立建模参考坐标系并提出建模假设条件，再建立柔性网模型以及空间绳网机器人的逼近动力学模型，最后进行系绳运动的速度跳变建模；由于空间绳网机器人在逼近目标过程中，四个单片网机构均处于绷紧状态，以平面有限元理论中的T3单元为基础，利用位置矢量的矩阵形式可以有效描述单片网机构的动力学。本发明的主要优势是可以有效降低单片网结构动力学模型复杂度，提高模型解算速度。

进一步的，本发明步骤3)，利用系绳杨氏模量极高的特点，将系绳建模为不可伸长；并利用间隙函数和约束反力建模系绳松弛时张力为零的特点。模型稳定度高，计算速度快，可以避免系绳的珠式模型(质点弹簧模型)仿真步长小，计算缓慢，易于发散的问题。

进一步的，本发明步骤4)，有效反映了系绳绷紧时的动力学特性，将珠式模型中极短时间发生的动力学用瞬态形式进行描述，简化建模复杂度。同步骤3)一起，从根本上解决了珠式模型计算步长必须很小，否则极易发散的问题。

附图说明

图1为空间绳网机器人的结构示意图；

图2为空间绳网机器人的单片柔性网机构示意图；

图3为系绳绷紧瞬间，两端点的速度跳变示意图。

其中，1为空间平台；2为连接系绳；3为柔性网；4为连接点；5为第一自主机动单元；6为第二自主机动单元；7为第三自主机动单元；8为第四自主机动单元。

具体实施方式

下面结合附图对本发明作进一步详细的说明：

参见图1至图3，本发明一种空间绳网机器人的逼近动力学建模方法，包括以下步骤：

一：建立建模参考坐标系并提出建模假设条件

首先建立地心惯性系OXYZ，其坐标原点O位于地球中心，OX轴指向地球的春分点，OZ轴指向地球北极，OY轴在赤道平面内垂直于OX轴；然后建立轨道坐标系O_oX_oY_oZ_o，其坐标原点O_o位于空间平台的质心，O_oX_o轴沿机动平台轨道的切向，指向平台运动方向，O_oZ_o轴与原点O_o和地心O的连线重合，指向地心，O_oY_o轴沿轨道平面负法线方向。坐标系示意图如附图1所示。

空间绳网机器人的逼近动力学建模基于如下假设：

(1)空间平台运行于圆轨道，且质量要远大于连接系绳、柔性网和自主机动单元质量的总和；

(2)忽略柔性网与平台连接系绳的弹性与质量，忽略自主机动单元体积，逼近任务中，连接系绳处于绷紧状态；

(3)柔性网的质量分布均匀，网孔很小，在逼近目标过程中不大幅变形；

(4)由于系绳杨氏模量极大，假设系绳连接点4与自主机动单元间、自主机动单元间的系绳不可伸长；

二：建立单片柔性网的模型

由于系绳连接节点与四个自主机动单元将柔性网分为四块，利用假设3，可将四块网结构建模为三角形薄壳。以系绳连接点与自主机动单元5,6组成的单片网结构为例，其建模示意如附图2所示。A表示系绳连接点4，B表示自主机动单元5，C表示自主机动单元6。R₁，R₂，R₃分别表示A,B,C在地心惯性系下的位置矢量。

针对单片网机构，采用平面有限元理论中的T3单元来描述。对于网结构上任一点D，其在地心惯性系下的位置矢量为：

R≈s₁R₁+s₂R₂+s₃R₃ (1)

式中，s₁、s₂、s₃表示点D在薄壳上的面积坐标，它们满足：

其中，表示柔性网在未发生任何变形条件下三角形的面积。

于是，建立连接点4与自主机动单元5，6组成的单片网结构的拉格朗日函数：

式中，dΣ表示三角形薄壳上的面积微元，ρ表示柔性网的平均面密度，m_W表示柔性网的总质量，G表示万有引力常数，M表示地球的质量。

对单片网结构的拉格朗日函数求变分，得：

其中，δ为变分符号，t₁,t₂表示积分时间。

利用C-W方程将上式转换到轨道坐标系下。

其中，r为网结构上点D在轨道坐标系O_oX_oY_oZ_o下的位置矢量。ω为空间平台轨道运动的平均角速度。

同理，可写出其他三片网结构的表达式。设r₄,r₅,r₆,r₇,r₈分别表示系绳连接点4、自主机动单元在轨道坐标系O_oX_oY_oZ_o下的位置矢量，将其写为矩阵形式：

则：

其中，L₄₆₇L₄₇₈L₄₈₅分别表示系绳连接点4与第二自主机动单元6和第三自主机动单元7，系绳连接点4与第三自主机动单元7和第四自主机动单元8，系绳连接点4与第四自主机动单元8和第一自主机动单元5组成的单片网结构的拉格朗日函数。表示矩阵的直积运算。

(lmn＝456,467,478,485)，M^lmn为5×5矩阵，其任意元素满足：

由于空间绳网机器人在逼近目标过程中，四个单片网机构均处于绷紧状态，以利用平面有限元理论中的T3单元为基础，利用式(5)可以有效描述单片网机构的动力学。本发明的主要优势是可以有效降低单片网结构动力学模型复杂度，提高模型解算速度。

三：建立空间绳网机器人的逼近动力学模型

首先分析系绳连接点与四个自主机动单元间、四个自主机动单元间系绳的不可伸长假设对逼近动力学建模带来的影响。

以系绳连接点4与自主机动单元5间系绳为例,设L₄₅为系绳连接点4与自主机动单元5之间系绳的标称长度，则，||r₄-r₅||≤L₄₅。引入间隙函数g₄₅对这一现象进行描述：

g₄₅＝L₄₅-||r₄-r₅||≥0 (7)

对应的约束反力λ₄₅满足：

后两种形式的约束反力与系统状态的二次导数相关，可将其写为隐式方程：

同理，相应写出系绳连接点4与自主机动单元6，系绳连接点4与自主机动单元7，系绳连接点4与自主机动单元8，自主机动单元5与自主机动单元6，自主机动单元6与自主机动单元7，自主机动单元7与自主机动单元8，自主机动单元8与自主机动单元5间系绳的间隙函数g₄₆,g₄₇,g₄₈,g₅₆,g₆₇,g₇₈,g₈₅和约束反力λ₄₆,λ₄₇,λ₄₈,λ₅₆,λ₆₇,λ₇₈,λ₈₅表达。

将间隙函数和约束反力写为矩阵，得：

g_N＝[g₄₅,g₄₆,g₄₇,g₄₈,g₅₆,g₆₇,g₇₈,g₈₅]^T

λ_N＝[λ₄₅,λ₄₆,λ₄₇,λ₄₈,λ₅₆,λ₆₇,λ₇₈,λ₈₅]^T

下面利用拉格朗日法对空间绳网机器人的逼近动力学进行建模。

整个系统的拉格朗日函数满足：

式中，m_P表示空间机动平台的质量，m_M表示自主机动单元的质量。R_o表示空间平台在地心惯性系下的坐标，R_i(i＝5,…,8)表示自主机动单元在地心惯性系下的坐标。

另外，系统中非保守力(包括约束反力)所做的功满足：

式中，r₄表示系绳连接点4在轨道坐标系O_oX_oY_oZ_o下的位置矢量，F_T表示空间平台和柔性网的连接系绳拉力，F_i(i＝5,…,8)表示作用在自主机动单元上的推力。

由广义Hamilton原理可知，空间绳网机器人的逼近动力学满足：

将式(10)和式(11)代入方程(12)并使用分步积分进行化简可得：

式中，

对(13)式进行整理，同时考虑变分的任意性可得系统的动力学方程为：

式中，

系统的解算可通过构建隐式方程，利用Lemke算法进行求解。

利用系绳杨氏模量极高的特点，将系绳建模为不可伸长；并利用间隙函数和约束反力建模系绳松弛时张力为零的特点。模型稳定度高，计算速度快，可以避免系绳的珠式模型(质点弹簧模型)仿真步长小，计算缓慢，易于发散的问题。

四、建立系绳连接点4与自主机动单元间系绳运动的速度跳变模型

如附图3所示，假设在t^-时刻，节点i和节点j之间的距离达到标称长度L_ij，而两个节点的相对速度不为0，且有相互远离的趋势。于是，在t时刻连接两个节点的系绳绷紧，由于系绳的刚度和阻尼都比较高，所以类似于碰撞过程，系绳会在极短时间内产生了比较大的张紧力，从而使得两个节点在t⁺时刻具有了沿系绳方向一致的速度。

为了描述瞬间的速度跳变，将动力学方程(14)在[t^-,t⁺]上积分，得：

由于在极短的时间内，位置矢量r_N还来不及发生变化，同时由于控制力F是有限的，因此方程左侧的后两项和右侧的第一项为0。对于那些满足系绳绷直条件的约束，即满足g_i＝0且其约束力在无限小的时间内趋向于无穷大，因此其积分不为0，而对于那些不满足绷直条件的约束，其约束反力为有限值，因而积分也为0。

为了便于表达，我们将满足绷直条件的约束取出，将他们的编号从小到大排列得到列向量按照同样的顺序排列对应的间隙函数和约束反力可以得到列向量和另外，分别将约束反力向量λ_NT和在时间段[t^-,t⁺]记作Λ_NT和两者之间满足：

其中，

于是，方程(15)可以化简为：

另外，对于间隙函数向量它在t⁺时刻必须满足：

则：绷紧系绳的速度跳变模型为：

本步骤的速度跳变建模有效反映了系绳绷紧时的动力学特性，将珠式模型中极短时间发生的动力学用瞬态形式进行描述，简化建模复杂度。同第三步一起，从根本上解决了珠式模型计算步长必须很小，否则极易发散的问题。

以上内容仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明权利要求书的保护范围之内。

Claims (5)

1.一种空间绳网机器人的逼近动力学建模方法，其特征在于，包括以下步骤：

1)建立建模参考坐标系并提出建模假设条

件；

建模假设条件具体如下：

假设1.空间平台(1)运行于圆轨道，且质量大于系绳、柔性网(3)和自主机动单元质量的总和；

假设2.忽略空间平台与柔性网的连接系绳(2)的弹性与质量，忽略自主机动单元体积，逼近任务中，连接系绳(2)处于绷紧状态；

假设3.柔性网的质量分布均匀，网孔很小，在逼近目标过程中不大幅变形；

假设4.由于系绳杨氏模量极大，假设系绳连接点(4)与自主机动单元间、自主机动单元间的系绳不可伸长；

2)建立单片柔性网的模型，具体方法是：

系绳连接点(4)与柔性网(3)的四个自主机动单元将柔性网分为四块，利用假设3，将四片柔性网分别建模为三角形薄壳；以系绳连接点(4)与柔性网(3)的第一自主机动单元(5)和第二自主机动单元(6)组成的单片柔性网为例说明，A表示系绳连接点(4)，B表示第一自主机动单元(5)，C表示第二自主机动单元(6)，R₁，R₂，R₃分别表示A,B,C在地心惯性系下的位置矢量；

针对单片柔性网，采用平面有限元理论中的T3单元来描述；对于单片柔性网上任一点D，其在地心惯性系下的位置矢量为：

R≈s₁R₁+s₂R₂+s₃R₃ (1)

式中，s₁、s₂、s₃表示点D在薄壳上的面积坐标，它们满足：

$s_1=\frac{{\stackrel{\‾}{S}}_{\mathrm{\Δ}BCD}}{{\stackrel{\‾}{S}}_{\mathrm{\Δ}ABC}},s_2=\frac{{\stackrel{\‾}{S}}_{\mathrm{\Δ}CAD}}{{\stackrel{\‾}{S}}_{\mathrm{\Δ}ABC}},s_3=\frac{{\stackrel{\‾}{S}}_{\mathrm{\Δ}ABD}}{{\stackrel{\‾}{S}}_{\mathrm{\Δ}ABC}}$

其中，表示单片柔性网在未发生任何变形条件下三角形的面积；

建立系绳连接点(4)与第一自主机动单元(5)和第二自主机动单元(6)组成的单片柔性网的拉格朗日函数：

$\begin{array}{c}{L}_{456}=\∫{\∫}_{\mathrm{\Δ}ABC}\frac{1}{2}\mathrm{\ρ}{\stackrel{\·}{R}}^T\stackrel{\·}{R}d\Σ-\∫{\∫}_{\mathrm{\Δ}ABC}(-\mathrm{\ρ}\frac{GM}{\left|\right|R\left|\right|})d\Σ\\ =\frac{m_W}{4}{\∫}_0^1{\mathrm{ds}}_1{\∫}_0^{1-s_1}(\frac{1}{2}{\stackrel{\·}{R}}^T\stackrel{\·}{R}+\mathrm{\ρ}\frac{GM}{\left|\right|R\left|\right|}){\mathrm{ds}}_2\end{array}---\left(2\right)$

式中，dΣ表示三角形薄壳上的面积微元，ρ表示单片柔性网的平均面密度，m_W表示整个柔性网的总质量，G表示万有引力常数，M表示地球的质量；

对单片柔性网的拉格朗日函数求变分，得：

$\mathrm{\δ}{\∫}_{t_1}^{t_2}{L}_{456}dt={\∫}_{t_1}^{t_2}\[\frac{m_W}{4}{\∫}_0^1{\mathrm{ds}}_1{\∫}_0^{1-s_1}{\mathrm{\δR}}^T(-\stackrel{\·\·}{R}-\mathrm{\ρ}\frac{GM}{\left|\right|R|{|}^3}R){\mathrm{ds}}_2\]dt---\left(3\right)$

其中，δ为变分符号，t₁,t₂表示积分时间；

利用C-W方程将上式转换到轨道坐标系下，得到：

${\∫}_{t_1}^{t_2}{\mathrm{\δL}}_{456}dt\≈{\∫}_{t_1}^{t_2}\[\frac{m_W}{4}{\∫}_0^1{\mathrm{ds}}_1{\∫}_0^{1-s_1}-{\mathrm{\δr}}^T(\stackrel{\·\·}{r}+M_{\stackrel{\·}{r}}\stackrel{\·}{r}+M_{r}r){\mathrm{ds}}_2\]dt---\left(4\right)$

其中，r为单片柔性网上点D在轨道坐标系O_oX_oY_oZ_o下的位置矢量；ω为空间平台(1)轨道运动的平均角速度；

同理，写出其他三片柔性网的表达式；设r₄,r₅,r₆,r₇,r₈分别表示系绳连接点(4)、第一自主机动单元(5)、第二自主机动单元(6)、第三自主机动单元(7)以及第四自主机动单元(8)在轨道坐标系O_oX_oY_oZ_o下的位置矢量，将其写为矩阵形式：

$r_N={\[r_4^T,r_5^T,r_6^T,r_7^T,r_8^T\]}^T$

则：

$\left\{\begin{array}{c}{\∫}_{t_1}^{t_2}{\mathrm{\δL}}_{456}dt\≈{\∫}_{t_1}^{t_2}-{\mathrm{\δr}}_N^T\[\frac{m_W}{4}(M_1^{456}{\stackrel{\·\·}{r}}_N+M_2^{456}{\stackrel{\·}{r}}_N+M_3^{456}{r}_N)\]dt\\ {\∫}_{t_1}^{t_2}{\mathrm{\δL}}_{467}dt\≈{\∫}_{t_1}^{t_2}-{\mathrm{\δr}}_N^T\[\frac{m_W}{4}(M_1^{467}{\stackrel{\·\·}{r}}_N+M_2^{467}{\stackrel{\·}{r}}_N+M_3^{467}{r}_N)\]dt\\ {\∫}_{t_1}^{t_2}{\mathrm{\δL}}_{478}dt\≈{\∫}_{t_1}^{t_2}-{\mathrm{\δr}}_N^T\[\frac{m_W}{4}(M_1^{478}{\stackrel{\·\·}{r}}_N+M_2^{478}{\stackrel{\·}{r}}_N+M_3^{478}{r}_N)\]dt\\ {\∫}_{t_1}^{t_2}{\mathrm{\δL}}_{485}dt\≈{\∫}_{t_1}^{t_2}-{\mathrm{\δr}}_N^T\[\frac{m_W}{4}(M_1^{485}{\stackrel{\·\·}{r}}_N+M_2^{485}{\stackrel{\·}{r}}_N+M_3^{485}{r}_N)\]dt\end{array}\right.---\left(5\right)$

其中，L₄₆₇、L₄₇₈、L₄₈₅分别表示系绳连接点(4)与第二自主机动单元(6)、第三自主机动单元(7)，系绳连接点(4)与第三自主机动单元(7)、第四自主机动单元(8)，系绳连接点(4)与第四自主机动单元(8)、第一自主机动单元(5)分别组成的三个单片柔性网的拉格朗日函数；表示矩阵的直积运算；

(lmn＝456,467,478,485)，M^lmn为5×5矩阵，其任意元素满足：

3)建立空间绳网机器人的逼近动力学模型；

4)建立系绳连接点与自主机动单元间系绳运动的速度跳变模型。

2.根据权利要求1所述的空间绳网机器人的逼近动力学建模方法，其特征在于：所述的步骤1)中，建立建模参考坐标系的具体方法是：

首先建立地心惯性系OXYZ，其坐标原点O位于地球中心，OX轴指向地球的春分点，OZ轴指向地球北极，OY轴在赤道平面内垂直于OX轴；然后建立轨道坐标系O_oX_oY_oZ_o，其坐标原点O_o位于空间平台(1)的质心，O_oX_o轴沿机动空间平台(1)轨道的切向，指向空间平台(1)运动方向，O_oZ_o轴与原点O_o和地心O的连线重合，指向地心，O_oY_o轴沿轨道平面负法线方向。

3.根据权利要求1所述的空间绳网机器人的逼近动力学建模方法，其特征在于：所述的步骤3)中，建立空间绳网机器人的逼近动力学模型的具体方法是：

首先分析系绳连接点与四个自主机动单元间、四个自主机动单元间系绳的不可伸长假设对逼近动力学建模带来的影响；

以系绳连接点(4)和第一自主机动单元(5)间系绳为例，设L₄₅为连接点(4)和第一自主机动单元(5)之间系绳的标称长度，则，||r₄-r₅||≤L₄₅；引入间隙函数g₄₅对这一现象进行描述：

g₄₅＝L₄₅-||r₄-r₅||≥0 (7)

对应的约束反力λ₄₅满足：

${\mathrm{\λ}}_{45}=\left\{\begin{array}{cc}0& g_{45}>0\\ 0& g_{45}=0,{\stackrel{\·}{g}}_{45}>0\\ 0& g_{45}=0,{\stackrel{\·}{g}}_{45}=0,{\stackrel{\·\·}{g}}_{45}>0\\ \≥0& g_{45}=0,{\stackrel{\·}{g}}_{45}=0,{\stackrel{\·\·}{g}}_{45}=0\end{array}\right.---\left(8\right)$

后两种形式的约束反力与系统状态的二次导数相关，将其写为隐式方程：

$\left\{\begin{array}{c}{g}_{45}=0,{\stackrel{\·}{g}}_{45}=0\\ {\stackrel{\·\·}{g}}_{45}{\mathrm{\λ}}_{45}=0,{\stackrel{\·\·}{g}}_{45}\≥0,{\mathrm{\λ}}_{45}\≥0\end{array}\right.---\left(9\right)$

同理，写出系绳连接点(4)与第二自主机动单元(6)，连接点(4)与第三自主机动单元(7)，连接点(4)与第四自主机动单元(8)，第一自主机动单元(5)与第二自主机动单元(6)，第二自主机动单元(6)与第三自主机动单元(7)，第三自主机动单元(7)与第四自主机动单元(8)，第四自主机动单元(8)与第一自主机动单元(5)之间系绳的间隙函数(g₄₆,g₄₇,g₄₈,g₅₆,g₆₇,g₇₈,g₈₅)和约束反力λ₄₆,λ₄₇,λ₄₈,λ₅₆,λ₆₇,λ₇₈,λ₈₅表达；

将间隙函数和约束反力写为矩阵，得：

g_N＝[g₄₅,g₄₆,g₄₇,g₄₈,g₅₆,g₆₇,g₇₈,g₈₅]^T

λ_N＝[λ₄₅,λ₄₆,λ₄₇,λ₄₈,λ₅₆,λ₆₇,λ₇₈,λ₈₅]^T

利用拉格朗日法对空间绳网机器人的逼近动力学进行建模

整个系统的拉格朗日函数满足：

$L=\frac{1}{2}{m}_P{\stackrel{\·}{R}}_o^T{\stackrel{\·}{R}}_o-\frac{{\mathrm{GMm}}_P}{\left|\right|R_o\left|\right|}+L_{456}+L_{467}+L_{478}+L_{485}+\underset{i=5}{\overset{8}{\Σ}}{m}_M\[\frac{1}{2}{R}_i^T{R}_i-\frac{GM}{\left|\right|R_i\left|\right|}\]---\left(10\right)$

式中，m_P表示空间平台(1)的质量，m_M表示自主机动单元的质量；R_o表示空间平台(1)在地心惯性系下的坐标，R_i(i＝5,…,8)表示第一自主机动单元(5)、第二自主机动单元(6)、第三自主机动单元(7)以及第四自主机动单元(8)在地心惯性系下的坐标；

另外，系统中非保守力所做的功满足：

${\mathrm{\δ}}^{\′}W=-{\mathrm{\δr}}_4^T\frac{r_4}{\left|\right|r_4\left|\right|}{F}_T+\underset{i=5}{\overset{8}{\Σ}}{\mathrm{\δR}}_i^T{F}_i+{\mathrm{\δg}}_N{\mathrm{\λ}}_N---\left(11\right)$

式中，r₄表示系绳连接点(4)在轨道坐标系O_oX_oY_oZ_o下的位置矢量，F_T表示空间平台和柔性网的连接系绳拉力，F_i(i＝5,…,8)表示作用在自主机动单元上的推力；

由广义Hamilton原理可知，空间绳网机器人的逼近动力学满足：

${\∫}_{t_0}^{t_1}(\mathrm{\δ}L+{\mathrm{\δ}}^{\′}W)dt=0---\left(12\right)$

将式(10)和式(11)代入方程(12)并使用分步积分进行化简可得：

${\∫}_{t_0}^{t_1}\{-{\mathrm{\δr}}_N^T\[\frac{m_W}{4}(M_1^W{\stackrel{\·\·}{r}}_N+M_2^W{\stackrel{\·}{r}}_N+M_3^W{r}_N)\]-\underset{i=5}{\overset{8}{\Σ}}{\mathrm{\δr}}_i^T\[m_M({\stackrel{\·\·}{r}}_i+M_{\stackrel{\·}{r}}{\stackrel{\·}{r}}_i+M_r{r}_i)\]+{\mathrm{\δr}}_N^{T}F+{\mathrm{\δr}}_N^T{\left(\frac{\∂g_N}{\∂r_N}\right)}^T{\mathrm{\λ}}_N\}dt=0---\left(13\right)$

式中，

$M_i^W=M_i^{456}+M_i^{467}+M_i^{478}+M_i^{485},(i=1,2,3),$

$F={\[-F_T{r}_4^T/\left|\right|r_4\left|\right|,F_5^T,F_6^T,F_7^T,F_8^T\]}^T,$

对(13)式进行整理，同时考虑变分的任意性得系统的动力学方程为：

$M{\stackrel{\·\·}{r}}_N+C{\stackrel{\·}{r}}_N+{\mathrm{Kr}}_N=F+{\left(\frac{\∂g_N}{\∂r_N}\right)}^T{\mathrm{\λ}}_N---\left(14\right)$

式中，

$M=\frac{m_W}{4}{M}_1^W+m_{M}diag(0,1,1,1,1)\⊗diag(1,1,1),$

$C=\frac{m_W}{4}{M}_2^W+m_{M}diag(0,1,1,1,1)\⊗M_{\stackrel{\·}{r}},$

$K=\frac{m_W}{4}{M}_2^W+m_{M}diag(0,1,1,1,1)\⊗M_r.$

4.根据权利要求3所述的空间绳网机器人的逼近动力学建模方法，其特征在于：整个系统的解算通过构建隐式方程，利用Lemke算法进行求解。

5.根据权利要求3所述的空间绳网机器人的逼近动力学建模方法，其特征在于：所述的步骤4)中，系绳运动的速度跳变建模的具体方法是：

假设在t^-时刻，节点i和节点j之间的距离达到标称长度L_ij，而两个节点的相对速度不为0，且有相互远离的趋势；于是，在t时刻连接两个节点的系绳绷紧，由于系绳的刚度和阻尼都比较高，所以类似于碰撞过程，系绳会在极短时间内产生了比较大的张紧力，从而使得两个节点在t⁺时刻具有了沿系绳方向一致的速度；

为了描述瞬间的速度跳变，将动力学方程(14)在[t^-,t⁺]上积分，得：

$\underset{t^{-}}{\overset{t^{+}}{\∫}}M{\stackrel{\·\·}{r}}_{N}dt+\underset{t^{-}}{\overset{t^{+}}{\∫}}C{\stackrel{\·}{r}}_{N}dt+\underset{t^{-}}{\overset{t^{+}}{\∫}}{\mathrm{Kr}}_{N}dt=\underset{t^{-}}{\overset{t^{+}}{\∫}}Fdt+{\left(\frac{\∂g_N}{\∂r_N}\right)}^T\underset{t^{-}}{\overset{t^{+}}{\∫}}{\mathrm{\λ}}_{N}dt---\left(15\right)$

由于在极短的时间内，位置矢量r_N还来不及发生变化，同时由于控制力F是有限的，因此方程左侧的后两项和右侧的第一项为0；对于那些满足系绳绷直条件的约束，即满足g_i＝0且其约束力在无限小的时间内趋向于无穷大，因此其积分不为0，而对于那些不满足绷直条件的约束，其约束反力为有限值，因而积分也为0；

为了便于表达，将满足绷直条件的约束取出，将其编号从小到大排列得到列向量按照同样的顺序排列对应的间隙函数和约束反力能够得到列向量和另外，分别将约束反力向量λ_NT和在时间段[t^-,t⁺]记作Λ_NT和两者之间满足：

${\mathrm{\Λ}}_{NT}=M_{RT}{\stackrel{\‾}{\mathrm{\Λ}}}_{NT}---\left(16\right)$

其中，

$M_{RT}(i,j)=\left\{\begin{array}{cc}1& i={\stackrel{\‾}{I}}_{NI}\left(j\right)\\ 0& others\end{array}\right.$

于是，方程(15)能够化简为：

$M({\stackrel{\·}{r}}_N{|}_{t^{+}}-{\stackrel{\·}{r}}_N{|}_{t^{-}})={\left(\frac{\∂g_N}{\∂r_N}\right)}^T{M}_{RT}{\stackrel{\‾}{\mathrm{\Λ}}}_{NT}---\left(17\right)$

另外，对于间隙函数向量它在t⁺时刻必须满足：

${\stackrel{\·}{\stackrel{\‾}{g}}}_{NT}{|}_{t^{+}}=0---\left(18\right)$

则：绷紧系绳的速度跳变模型为：

$\left\{\begin{array}{c}M({\stackrel{\·}{r}}_N{|}_{t^{+}}-{\stackrel{\·}{r}}_N{|}_{t^{-}})={\left(\frac{\∂g_N}{\∂r_N}\right)}^T{M}_{RT}{\stackrel{\‾}{\mathrm{\Λ}}}_{NT}\\ {\stackrel{\·}{\stackrel{\‾}{g}}}_{NT}{|}_{t^{+}}=0\end{array}\right.---\left(19\right).$