空间机器人接触操作的风险判定方法、系统、介质及设备

1.本发明涉及空间机器人接触碰撞动力学领域，尤其涉及空间机器人接触操作的风险判定方法、系统、介质及设备。


背景技术：

2.空间机器人在轨执行接触性操作任务时，主要关注最大接触碰撞力和基座姿态扰动两个指标。由于传感器误差以及目标飞行器在轨速度的影响，采用零相对速度的软接触在实际工程中很难实现，接触操作不可避免的会产生接触碰撞力，过大接触力可能会对目标物或者机器人造成结构损伤导致任务失败；针对基座惯量较小的空间机器人系统，接触碰撞会引起基座姿态较大的扰动，而基座姿态往往与空间机器人对地通信和能量获取密切相关。因此，空间机器人在轨执行接触性操作任务时，很有必要对上述两个指标同时进行关注，并根据分析结果对潜在风险预警。现有的研究主要针对空间机器人接触碰撞力和基座姿态扰动单一目标进行优化，优化目标以“尽可能小”等定性语句阐述，然而由于这两个重要指标不具有一致性甚至空间机器人某些构型下存在相反的可能，即一种指标得到优化另外一种指标可能会变差，这种定性的优化不能确保通过优化一种指标使得两种指标同时满足要求。本发明提出一种自由漂浮空间机器人接触操作风险评估方法，基于两种指标形成风险预判区域，可客观定量描述空间机器人接触操作风险程度，同时可为风险指标优化提供直观定量的参考目标。
3.目前在机器人碰撞风险评估方面，主要有基于轨迹预测的风险评估方法、基于传感器的风险评估方法和基于障碍物运动姿态信息的风险评估方法，此类机器人碰撞风险评估方法旨在避免机器人与目标物发生不当接触碰撞而产生风险，其核心均为避免接触。而针对空间机器人系统，接触是其执行在轨操作的前提和基础，因此上述风险评估方法并不适用。经过广泛的调研目前尚未见到针对空间机器人接触操作风险评估的方法，大多数学者致力于对空间机器人最大接触碰撞力(接触碰撞脉冲)和基座姿态扰动等单个风险指标进行推导和优化，优化目标以“尽可能小”等定性语句阐述，然而由于这两个重要指标不具有一致性甚至空间机器人某些构型下存在相反的可能，即一种指标得到优化另外一种指标可能会变差，这种定性的优化不能确保通过优化一种指标使得两种指标同时满足要求。此外，当某些指标本就在安全区域内时，没有必要浪费太空中宝贵的资源对其进行优化。


技术实现要素：

4.本发明所要解决的技术问题是提供空间机器人接触操作的风险判定方法、系统、介质及设备。
5.本发明解决上述技术问题的技术方案如下：一种空间机器人接触操作的风险判定方法，包括：
6.步骤1，输入模块获取空间机器人的风险判断参数；
7.步骤2，风险评估模块基于所述风险判断参数计算最大接触碰撞力以及基座姿态
扰动二范数；
8.步骤3，输出模块基于所述最大接触碰撞力以及所述基座姿态扰动二范数得出风险等级判定结果。
9.本发明的有益效果是：通过最大接触碰撞力以及基座姿态扰动二范数的获取可以避免机器人与目标物发生不当接触碰撞而产生风险，且，本技术结合了上述两个条件，相较于现有技术中仅针对一种条件进行风险评估的判断方法根据有准确性，另外还可以有效减少太空中宝贵的资源的浪费情况。
10.在上述技术方案的基础上，本发明还可以做如下改进。
11.进一步，步骤2具体包括：
12.风险评估模块根据接触碰撞力f的计算公式对所述风险判断参数进行计算得到最大接触碰撞力，根据姿态变化量计算公式对所述风险判断参数进行计算得到基座姿态扰动二范数；
13.其中，所述接触碰撞力f的计算公式的建立过程具体为：
14.风险评估模块基于所述空间机器人的关节角度以及基座位姿构建动力学模型，基于所述动力学模型构建所述空间机器人的改进等效质量模型，基于接触碰撞力模型以及所述改进等效质量模型构建接触碰撞力f的计算公式；
15.所述姿态变化量计算公式的建立过程具体为：
16.基于所述动力学模型，得出姿态变化量计算公式。
17.进一步，所述动力学模型具体为：
[0018][0019]
其中，为所述空间机器人的基座的线加速度，为所述空间机器人的基座的角加速度，为所有关节角的角加速度，(
·
)
t
代表矩阵的转置；c＝[cb,cm]
t
，cb和cm分别代表所述空间机器人的基座和所述空间机器人的机械臂相关的非线性速度依赖项；f
b_m
＝[f
bt
,τ
mt
]
t
，fb＝[f
bt
,τ
bt
]
t
，fb和τb分别为基座所受外力和外力矩，τm为所述机械臂所有关节的输出力矩；fe代表施加在机械臂末端的力或力矩；j
b_m
代表空间机器人系统增广雅克比矩阵。
[0020]
进一步，所述基于所述动力学模型构建所述空间机器人的改进等效质量模型具体为：
[0021]
风险评估模块构建空间机器人的机械臂末端的速度计算公式，基于所述速度计算公式以及所述动力学模型构建改进等效质量模型，所述改进等效质量模型具体为：
[0022][0023]
其中，u代表所述机械臂末端施加力的单位方向矢量，j
b_mv
代表空间机器人系统线速度增广雅克比矩阵。
[0024]
进一步，所述f计算公式具体为：
[0025][0026]
其中，为机械臂末端与目标物接接触碰撞初始时刻相对速度，cr为接触环境恢复系数，k为接触环境刚度系数，δ为机械臂末端与目标物接触碰撞压缩量，α为常量，λ为接触环境磁滞阻尼系数，为机械臂末端与目标物接触碰撞压缩速度。
[0027]
进一步，所述步骤3具体为：
[0028]
输出模块基于所述最大接触碰撞力、所述基座姿态扰动二范数以及风险阈值进行风险等级判定；
[0029]
其中，所述风险阈值通过机械设计标准以及任务需求设计划分得出。
[0030]
进一步，所述风险判断参数包括：
[0031]
空间机器人的基座位姿、关节角度、机械臂末端速度，目标物的质量以及速度，接触环境的杨氏模量和泊松比、接触点处曲率半径和接触碰撞力方向。
[0032]
本发明解决上述技术问题的另一种技术方案如下：一种空间机器人接触操作的风险判定系统，包括：
[0033]
输入模块，用于获取空间机器人的风险判断参数；
[0034]
风险评估模块，用于基于所述风险判断参数计算最大接触碰撞力以及基座姿态扰动二范数；
[0035]
输出模块，用于基于所述最大接触碰撞力以及所述基座姿态扰动二范数得出风险等级判定结果。
[0036]
本发明的有益效果是：通过最大接触碰撞力以及基座姿态扰动二范数的获取可以避免机器人与目标物发生不当接触碰撞而产生风险，且，本技术结合了上述两个条件，相较于现有技术中仅针对一种条件进行风险评估的判断方法根据有准确性，另外还可以有效减少太空中宝贵的资源的浪费情况。
[0037]
进一步，风险评估模块具体用于：
[0038]
根据接触碰撞力f的计算公式对所述风险判断参数进行计算得到最大接触碰撞力，根据姿态变化量计算公式对所述风险判断参数进行计算得到基座姿态扰动二范数；
[0039]
其中，所述接触碰撞力f的计算公式的建立过程具体为：
[0040]
风险评估模块基于所述空间机器人的关节角度以及基座位姿构建动力学模型，基于所述动力学模型构建所述空间机器人的改进等效质量模型，基于接触碰撞力模型以及所述改进等效质量模型构建接触碰撞力f的计算公式；
[0041]
所述姿态变化量计算公式的建立过程具体为：
[0042]
基于所述动力学模型，得出姿态变化量计算公式。
[0043]
进一步，所述动力学模型具体为：
[0044][0045]
其中，为所述空间机器人的基座的线加速度，为所述空间机器人的基座的角加速度，为所有关节角的角加速度，(
·
)
t
代表矩阵的转置；c＝[cb,cm]
t
，cb和cm分别代表所述空间机器人的基座和所述空间机器人的机械臂相
关的非线性速度依赖项；f
b_m
＝[f
bt
,τ
mt
]
t
，fb＝[f
bt
,τ
bt
]
t
，fb和τb分别为基座所受外力和外力矩，τm为所述机械臂所有关节的输出力矩；fe代表施加在机械臂末端的力或力矩；j
b_m
代表空间机器人系统增广雅克比矩阵。
[0046]
进一步，所述基于所述动力学模型构建所述空间机器人的改进等效质量模型具体为：
[0047]
风险评估模块构建空间机器人的机械臂末端的速度计算公式，基于所述速度计算公式以及所述动力学模型构建改进等效质量模型，所述改进等效质量模型具体为：
[0048][0049]
其中，u代表所述机械臂末端施加力的单位方向矢量，j
b_mv
代表空间机器人系统线速度增广雅克比矩阵。
[0050]
进一步，所述f计算公式具体为：
[0051][0052]
其中，为机械臂末端与目标物接接触碰撞初始时刻相对速度，cr为接触环境恢复系数，k为接触环境刚度系数，δ为机械臂末端与目标物接触碰撞压缩量，α为常量，λ为接触环境磁滞阻尼系数，为机械臂末端与目标物接触碰撞压缩速度。
[0053]
进一步，输出模块具体用于：
[0054]
基于所述最大接触碰撞力、所述基座姿态扰动二范数以及风险阈值进行风险等级判定；
[0055]
其中，所述风险阈值通过机械设计标准以及任务需求设计划分得出。
[0056]
进一步，所述风险判断参数包括：
[0057]
空间机器人的基座位姿、关节角度、机械臂末端速度，目标物的质量以及速度，接触环境的杨氏模量和泊松比、接触点处曲率半径和接触碰撞力方向
[0058]
本发明解决上述技术问题的另一种技术方案如下：一种存储介质，所述存储介质中存储有指令，当计算机读取所述指令时，使所述计算机执行如上述任一项所述的一种空间机器人接触操作的风险判定方法。
[0059]
本发明的有益效果是：通过最大接触碰撞力以及基座姿态扰动二范数的获取可以避免机器人与目标物发生不当接触碰撞而产生风险，且，本技术结合了上述两个条件，相较于现有技术中仅针对一种条件进行风险评估的判断方法根据有准确性，另外还可以有效减少太空中宝贵的资源的浪费情况。
[0060]
本发明解决上述技术问题的另一种技术方案如下：一种电子设备，包括上述存储介质、执行上述存储介质内的指令的处理器。
[0061]
本发明的有益效果是：通过最大接触碰撞力以及基座姿态扰动二范数的获取可以避免机器人与目标物发生不当接触碰撞而产生风险，且，本技术结合了上述两个条件，相较于现有技术中仅针对一种条件进行风险评估的判断方法根据有准确性，另外还可以有效减少太空中宝贵的资源的浪费情况。
附图说明
[0062]
图1为本发明一种空间机器人接触操作的风险判定方法的实施例提供的流程示意图；
[0063]
图2为本发明一种空间机器人接触操作的风险判定系统的实施例提供的主体结构框图；
[0064]
图3为本发明一种空间机器人接触操作的风险判定方法的实施例提供的完整结构示意图；
[0065]
图4为本发明一种空间机器人接触操作的风险判定方法的实施例提供的二维风险评估空间示意图；
[0066]
图5为本发明一种空间机器人接触操作的风险判定方法的实施例提供的三自由度自由漂浮空间机器人结构示意图；
[0067]
图6为本发明一种空间机器人接触操作的风险判定方法的实施例提供的风险评估结果示意图。
具体实施方式
[0068]
以下对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
[0069]
如图1所示，一种空间机器人接触操作的风险判定方法，包括：
[0070]
步骤1，输入模块100获取空间机器人的风险判断参数；
[0071]
步骤2，风险评估模块200基于所述风险判断参数计算最大接触碰撞力以及基座姿态扰动二范数；
[0072]
步骤3，输出模块300基于所述最大接触碰撞力以及所述基座姿态扰动二范数得出风险等级判定结果。
[0073]
在一些可能的实施方式中，通过最大接触碰撞力以及基座姿态扰动二范数的获取可以避免机器人与目标物发生不当接触碰撞而产生风险，且，本技术结合了上述两个条件，相较于现有技术中仅针对一种条件进行风险评估的判断方法根据有准确性，另外还可以有效减少太空中宝贵的资源的浪费情况。
[0074]
需要说明的是，该风险评估方法由三个模块组成：输入模块100、风险评估模块200和输出模块300，如图3所示。输入模块100包括空间机器人运动学和动力学参数、基座位姿和关节角度、机械臂末端速度，目标物质量和速度(假设目标物为单体，且仅发生对心碰撞)，接触环境的杨氏模量和泊松比、接触点处曲率半径和接触碰撞力方向；风险评估模块200包括最大接触碰撞力和基座姿态扰动两种风险评估指标，依据机械设计标准和任务需求设计的可允许最大接触碰撞力和可允许基座姿态扰动；输出模块300包括：风险等级和操作建议，其中风险等级包括四种情况：绝对安全、次安全、接近危险以及危险，根据风险等级的不同也会发出不同的警报，如声音或警报文件等，操作建议包括四种情况，直接操作无需优化、操作且注意观察某些风险指标、可操作性但某些风险指标建议进行优化以及停止操作但必须进行风险指标优化直至退出危险区域。
[0075]
其中，输入模块100包括的内容为风险判断参数，将上述参数输入至风险评估模块200中进行风险评估，评估中所利用到的公式或模型等均是预先建立好的，该部分可参考实
施例1进行辅助理解。经过风险评估模块200的计算后得出两个参数，即最大接触碰撞力以及所述基座姿态扰动二范数，输出模块300基于风险评估模块200计算的结果进一步进行判断，得出风险等级以及该风险等级对应的操作建议。
[0076]
实施例1，对于空间机器人的动力学模型部分，
[0077]
采用拉格朗日方法建立以关节角度和基座位姿为变量的动力学模型：
[0078][0079]
其中，为基座的线加速度，为基座的角加速度，为关节角加速度，(])
t
代表矩阵的转置；c＝[cb,cm]
t
，cb和cm分别代表基座和机械臂相关的非线性速度依赖项；f
b_m
＝[f
bt
,τ
mt
]
t
，fb＝[f
bt
,τ
bt
]
t
，fb和τb分别为基座所受外力和外力矩，τm为机械臂关节输出力矩；fe代表施加在基座和机械臂末端的力或力矩；j
b_m
代表空间机器人系统增广雅克比矩阵，h为空间机器人系统的惯性矩阵。该动力学模型为描述空间机器人的动力学模型，用以反应空间机器人系统的动力学特性，其中，关节角度是指机械臂各个臂杆的转动角度，具体见图5所标注；基座位姿是指基座的位置和姿态。a为基座质心到第一个关节的距离，b,c,d分别为机械臂臂杆1,2,3的长度，θ1,θ2,θ3分别为关节1，关节2和关节3的转角。以上均为空间机器人运动学参数，在计算动力学模型、改进等效质量、基座姿态扰动二范数都会用到，同时也是输入模块里的参数。
[0080]
对于风险指标1：最大接触碰撞力，
[0081]
由空间机器人速度级运动学知识可以建立：
[0082][0083]
其中，为空间机器人的机械臂末端的速度，v0为基座的线速度，ω0为基座的角速度，为关节角加速度。
[0084]
由于矩阵h可逆，在公式(1)两边同时乘以j
b_m
h-1
并结合公式(2)可得：
[0085][0086]
其中，为的广义逆矩阵，通常用来代表广义逆，或者伪逆。
[0087]
基于公式(3)，建立针对空间机器人的改进等效质量模型：
[0088][0089]
其中，u代表施加力的单位方向矢量，j
b_mv
代表线速度增广雅克比矩阵，me为机械臂末端改进等效质量。
[0090]
两物体之间连续型接触碰撞动力学经典模型为：
[0091]
[0092]
其中，k为刚度系数，为阻尼力项，δ为碰撞压缩量，α为常量一般取值为1.5，λ为磁滞阻尼系数可采用zhang et al模型如下：
[0093][0094]
其中，为接触碰撞初始时刻相对速度，cr为恢复系数。
[0095]
结合公式(4)可得：
[0096][0097]
其中，m
t
为目标物的质量，为机械臂末端与目标物接触碰撞压缩加速度。
[0098]
根据公式(7)进行进一步推导可得：
[0099][0100][0101]
其中，t为整个接触碰过程持续时间。
[0102]
结合公式(5)、(6)、(8)、(9)可得：
[0103][0104]
至此，可以得到空间机器人整个接触碰撞过程数据，通过对比可获得最大接触碰撞力。通过接触碰撞力公式(8)求取在整个接触碰撞过程中的不同位置力值进行比较获得的最大接触碰撞力数值解。
[0105]
对于风险指标2：最大基座姿态扰动
[0106]
将公式(1)拆分可得：
[0107][0108]
其中，m为空间机器人系统总质量，e为单位矩阵，h
bmν
和h
bmω
分别为空间机器人系统的线速度相关和空间机器人系统的角速度相关耦合惯性矩阵，hm为机械臂惯性矩阵，h
ω
为空间机器人的基座角速度相关惯性矩阵，r
0g
为从空间机器人的基座质心到空间机器人系统质心向量，r
0g
×
为r
0g
的反对称矩阵，c
bν
和c
bω
分别为空间机器人的基座线速度相关和角速度相关的非线性项，j
bν
和j
bω
分别为空间机器人的基座线速度和角速度雅克比矩阵，为所有
关节角的角加速度，jm为机械臂雅克比矩阵。n代表关节总个数。
[0109]
由公式(11)消除和可得：
[0110][0111]
其中，其中，其中，其中，
[0112]
针对自由漂浮空间机器人，由于基座不受外力和外力矩的作用，则fb＝[0,0,0]
t
,τb＝[0,0,0]
t
，对公式(12)在整个接触碰撞过程进行积分可得基座由接触碰撞引起的姿态变化量：
[0113][0114]
其中，t0为接触碰撞初始时刻，δt为接触碰撞持续时间，为基座角速度相关惯性矩阵。
[0115]
其中，机械设计标准和任务需求设计的可允许最大接触碰撞力和可允许基座姿态扰动均为认为定义的外界输入条件。
[0116]
实施例2，以最大接触碰撞力和基座姿态扰动形成的二维风险评估空间，x轴代表基座姿态扰动因子，y轴代表最大接触碰撞力因子，f
max
为最大接触碰撞力，||δω0||为基座姿态扰动二范数，s1是末端可允许的最大接触碰撞力，s2是可允许的最大基座扰动的二范数，故f
max
/s1∈[0,1]且||δω0||/s2∈[0,1]是确保接触操作任务安全的必要条件。η1,η2和η
′1,η
′2是根据机械设计标准和具体任务需求而预先设定的风险阈值，据此可以分为以下四个风险等级。s1以及s2均为人为设定的数据。
[0117]
风险等级i：f
max
/s1∈[0,η1]且||δω0||/s2∈[0,η
′1]；
[0118]
对于大多数操作任务，f
max
/s1和||δω0||/s2两个指标都是期望越小越好，因此，风险等级i意味着绝对安全，操作建议是直接操作无需优化。
[0119]
风险等级ii：f
max
/s1∈(η1,η2]且||δω0||/s2∈[0,η
′2]或者f
max
/s1∈[0,η2]且||δω0||/s2∈(η
′1,η
′2]；
[0120]
次安全区域，优化建议是操作，但需关注超出绝对安全区域的风险指标变化。
[0121]
风险等级iii：f
max
/s1∈(η2,1]且||δω0||/s2∈[0,1]或者f
max
/s1∈[0,1]且||δω0||/s2∈(η
′2,1]；
[0122]
当f
max
/s1→
1或者||δω0||/s2→
1时，意味着f
max
或者||δω0||已经趋近于其极限值了。因此，风险等级iii意味着接近风险，操作建议是可操作，但某些风险指标建议进行优化。
[0123]
风险等级iv：f
max
/s1∈(1, ∞)或者||δω0||/s2∈(1, ∞)；
[0124]
在此种条件下，f
max
或者||δω0||已经超出了其可允许范围，因此，风险等级iv意味着危险，操作建议是停止操作，必须进行风险指标优化直至退出危险区域。
[0125]
根据具体的任务需要，风险等级和风险指标还可以进行细分，形成更高维度的操作风险评估空间。
[0126]
优选地，在上述任意实施例中，步骤2具体包括：
[0127]
风险评估模块200根据接触碰撞力f的计算公式对所述风险判断参数进行计算得到最大接触碰撞力，根据姿态变化量计算公式对所述风险判断参数进行计算得到基座姿态扰动二范数；
[0128]
其中，所述接触碰撞力f的计算公式的建立过程具体为：
[0129]
风险评估模块200基于所述空间机器人的关节角度以及基座位姿构建动力学模型，基于所述动力学模型构建所述空间机器人的改进等效质量模型，基于接触碰撞力模型以及所述改进等效质量模型构建接触碰撞力f的计算公式；
[0130]
所述姿态变化量计算公式的建立过程具体为：
[0131]
基于所述动力学模型，得出姿态变化量计算公式。
[0132]
优选地，在上述任意实施例中，所述动力学模型具体为：
[0133][0134]
其中，为所述空间机器人的基座的线加速度，为所述空间机器人的基座的角加速度，为所有关节角的角加速度，(
·
)
t
代表矩阵的转置；c＝[cb,cm]
t
，cb和cm分别代表所述空间机器人的基座和所述空间机器人的机械臂相关的非线性速度依赖项；f
b_m
＝[f
bt
,τ
mt
]
t
，fb＝[f
bt
,τ
bt
]
t
，fb和τb分别为基座所受外力和外力矩，τm为所述机械臂所有关节的输出力矩；fe代表施加在机械臂末端的力或力矩；j
b_m
代表空间机器人系统增广雅克比矩阵。
[0135]
需要说明的是，n代表关节总个数，空间机器人系统为包含一条n自由度的串联机械臂和一个自由漂浮状态的基座组成的。
[0136]
优选地，在上述任意实施例中，所述基于所述动力学模型构建所述空间机器人的改进等效质量模型具体为：
[0137]
风险评估模块200建空间机器人的机械臂末端的速度计算公式，基于所述速度计算公式以及所述动力学模型构建改进等效质量模型，所述改进等效质量模型具体为：
[0138][0139]
其中，u代表所述机械臂末端施加力的单位方向矢量，j
b_mv
代表空间机器人系统线速度增广雅克比矩阵。
[0140]
优选地，在上述任意实施例中，所述接触碰撞力f的计算公式具体为：
[0141][0142]
其中，为机械臂末端与目标物接接触碰撞初始时刻相对速度，cr为接触环境恢复系数，k为接触环境刚度系数，δ为机械臂末端与目标物接触碰撞压缩量，α为常量，λ为接触环境磁滞阻尼系数，为机械臂末端与目标物接触碰撞压缩速度。
[0143]
优选地，在上述任意实施例中，所述步骤3具体为：
[0144]
输出模块300基于所述最大接触碰撞力、所述基座姿态扰动二范数以及风险阈值进行风险等级判定；
[0145]
其中，所述风险阈值通过机械设计标准以及任务需求设计划分得出。
[0146]
优选地，在上述任意实施例中，所述风险判断参数包括：
[0147]
空间机器人的基座位姿、关节角度、机械臂末端速度，目标物的质量以及速度，接触环境的杨氏模量和泊松比、接触点处曲率半径和接触碰撞力方向。
[0148]
如图2所示，一种空间机器人接触操作的风险判定系统，包括：
[0149]
输入模块100，用于获取空间机器人的风险判断参数；
[0150]
风险评估模块200，用于基于所述风险判断参数计算最大接触碰撞力以及基座姿态扰动二范数；
[0151]
输出模块300，用于基于所述最大接触碰撞力以及所述基座姿态扰动二范数得出风险等级判定结果。
[0152]
在一些可能的实施方式中，通过最大接触碰撞力以及基座姿态扰动二范数的获取可以避免机器人与目标物发生不当接触碰撞而产生风险，且，本技术结合了上述两个条件，相较于现有技术中仅针对一种条件进行风险评估的判断方法根据有准确性，另外还可以有效减少太空中宝贵的资源的浪费情况。
[0153]
优选地，在上述任意实施例中，风险评估模块200具体用于：
[0154]
根据接触碰撞力f的计算公式对所述风险判断参数进行计算得到最大接触碰撞力，根据姿态变化量计算公式对所述风险判断参数进行计算得到基座姿态扰动二范数；
[0155]
其中，所述接触碰撞力f的计算公式的建立过程具体为：
[0156]
风险评估模块200基于所述空间机器人的关节角度以及基座位姿构建动力学模型，基于所述动力学模型构建所述空间机器人的改进等效质量模型，基于接触碰撞力模型以及所述改进等效质量模型构建接触碰撞力f的计算公式；
[0157]
所述姿态变化量计算公式的建立过程具体为：
[0158]
基于所述动力学模型，得出姿态变化量计算公式。
[0159]
优选地，在上述任意实施例中，所述动力学模型具体为：
[0160][0161]
其中，为所述空间机器人的基座的线加速度，为所述空间机器人的基座的角加速度，为所有关节角的角加速度，(
·
)
t
代表矩阵的转置；c＝[cb,cm]
t
，cb和cm分别代表所述空间机器人的基座和所述空间机器人的机械臂相关的非线性速度依赖项；f
b_m
＝[f
bt
,τ
mt
]
t
，fb＝[f
bt
,τ
bt
]
t
，fb和τb分别为基座所受外力和外力矩，τm为所述机械臂所有关节的输出力矩；fe代表施加在机械臂末端的力或力矩；j
b_m
代表空间机器人系统增广雅克比矩阵。
[0162]
优选地，在上述任意实施例中，所述基于所述动力学模型构建所述空间机器人的改进等效质量模型具体为：
[0163]
风险评估模块200构建空间机器人的机械臂末端的速度计算公式，基于所述速度计算公式以及所述动力学模型构建改进等效质量模型，所述改进等效质量模型具体为：
[0164][0165]
其中，u代表所述机械臂末端施加力的单位方向矢量，j
b_mv
代表空间机器人系统线速度增广雅克比矩阵。
[0166]
优选地，在上述任意实施例中，所述接触碰撞力f的计算公式具体为：
[0167][0168]
其中，为机械臂末端与目标物接接触碰撞初始时刻相对速度，cr为接触环境恢复系数，k为接触环境刚度系数，δ为机械臂末端与目标物接触碰撞压缩量，α为常量，λ为接触环境磁滞阻尼系数，为机械臂末端与目标物接触碰撞压缩速度。
[0169]
优选地，在上述任意实施例中，输出模块300具体用于：
[0170]
基于所述最大接触碰撞力、所述基座姿态扰动二范数以及风险阈值进行风险等级判定；
[0171]
其中，所述风险阈值通过机械设计标准以及任务需求设计划分得出。
[0172]
优选地，在上述任意实施例中，所述风险判断参数包括：
[0173]
空间机器人的基座位姿、关节角度、机械臂末端速度，目标物的质量以及速度，接触环境的杨氏模量和泊松比、接触点处曲率半径和接触碰撞力方向。
[0174]
需要说明的是，目标物为太空中的待操作物体，包括目标飞行器、失效卫星、空间碎片等。
[0175]
本发明解决上述技术问题的另一种技术方案如下：一种存储介质，所述存储介质中存储有指令，当计算机读取所述指令时，使所述计算机执行如上述任一项所述的一种空间机器人接触操作的风险判定方法。
[0176]
在一些可能的实施方式中，通过最大接触碰撞力以及基座姿态扰动二范数的获取可以避免机器人与目标物发生不当接触碰撞而产生风险，且，本技术结合了上述两个条件，相较于现有技术中仅针对一种条件进行风险评估的判断方法根据有准确性，另外还可以有效减少太空中宝贵的资源的浪费情况。
[0177]
本发明解决上述技术问题的另一种技术方案如下：一种电子设备，包括上述存储介质、执行上述存储介质内的指令的处理器。
[0178]
在一些可能的实施方式中，通过最大接触碰撞力以及基座姿态扰动二范数的获取可以避免机器人与目标物发生不当接触碰撞而产生风险，且，本技术结合了上述两个条件，相较于现有技术中仅针对一种条件进行风险评估的判断方法根据有准确性，另外还可以有效减少太空中宝贵的资源的浪费情况。
[0179]
实施例3，根据本发明提供的一种自由漂浮空间机器人接触操作风险评估方法，以如附图3所示的三自由度自由漂浮空间机器人为例进行说明，其中a＝0.6m,b＝1.3m,c＝1.5m,d＝0.5m，动力学参数见表1，从与基座固连的连杆到机械臂末端连杆依次编号为连杆1，连杆2和连杆3。
[0180]
表1动力学参数
[0181][0182]
假设接触环境刚度参数为k＝109n/m
1.5
，恢复系数为cr＝0.8，初始接触速度为[0,0,-0.1]m/s，机器人捕获构型为θ＝[60,235.68,154.32]deg，经风险评估模型计算接触操作风险指标：最大接触碰撞力为1000.88n，基座姿态扰动量的二范数为1.45
×
10-2
deg/s。根据机械结构要求和任务需求，设置η1＝0.3,η
′1＝0.5,η2＝0.7,η
′2＝0.8，则共存在四个风险等级如附图6所示。当s1＝5000n,s2＝0.4deg/s时，风险落点为p1在绝对安全区域内，则此次接触操作任务可直接执行无需优化；当s1＝2000n,s2＝0.2deg/s时，风险落点为p2在次安全区域内，则此次接触性操作任务可执行，但需关注风险指标基座姿态扰动；当s1＝1200n,s2＝0.4deg/s时，风险落点为p3在接近危险区域内，任务可执行但建议对最大接触碰撞力进行优化；当s1＝5000n,s2＝0.1deg/s时，风险落点为p4在危险区域内，则必须停止操作，优化基座姿态扰动直到使风险落点落在安全区域内。
[0183]
读者应理解，在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
[0184]
在本技术所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的方法实施例仅仅是示意性的，例如，步骤的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个步骤可以结合或者可以集成到另一个步骤，或一些特征可以忽略，或不执行。
[0185]
上述方法如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-onlymemory)、随机存取存储器(ram，randomaccessmemory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0186]
以上，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要
求的保护范围为准。