一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法与流程

本发明属于空间非合作目标抓捕控制领域，具体涉及一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法。

背景技术：

随着空间开发利用活动的增加，空间碎片的数量急剧增多。而空间轨道的资源是有限的，如果不清理空间碎片，待其数量增大到一定程度后，已有碎片之间的碰撞会使产生碎片的速度大于清理的速度，从而引起“凯斯勒”效应，空间将不能再被利用。因此空间碎片的清理是亟待解决的任务。由于大量空间碎片具有旋转特性，在抓捕过程中，不可避免会出现碰撞冲击等。因此需要针对空间碎片抓捕过程中的振动冲击控制问题，提出高效可行的解决方法。

技术实现要素：

本发明的目的在于提供一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法。

本发明采用如下技术方案来实现的：

一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法，该方法针对空间非合作目标抓捕过程中的碰撞冲击问题，通过在抓捕机构以及抓捕航天器之间装置仿生非线性结构方法实现抓捕被动减振，减小非合作目标抓捕过程中的冲击碰撞对基座的扰动；具体来说，针对空间非合作目标抓捕过程的冲击减振问题，考虑到仿生结构阻尼的质量和减振效果，对非线性结构参数进行优化设计，求解最优冲击碰撞扰动被动抑制结构设计。

本发明进一步的改进在于，根据生物界普遍存在的x构型结构特征，设计一种具有非线性阻尼特征的减振机构，该减振机构由多段x型结构组成，每个x型结构由两根杆组成；杆件之间的摩擦以及结构形变可有效减缓抓捕过程扰动碰撞冲击对抓捕航天器的影响。

本发明进一步的改进在于，包括以下步骤：

步骤一：减振机构结构设计

步骤二：减振机构动力学建模

设四自由度圆盘中“x”构型数目为n′，中间单自由度部分“x”构型数目为n，θ为x构型与底面相连的部分与底盘之间的初始夹角；为x构型工作之后的位置与其初始位置之间的夹角；qi为x构型总体的纵向变化量；si为x构型总体的横向变化量；l为x构型每一部分的杆的长度，由几何关系可得到如下关系式：

由公式(1)可求得：

而参数si可表示为

整个系统的动能表示如下，其中m为整个系统的质量；为z轴方向上的速度：

系统势能表示如下，其中ωmp(t)为上盘的角速度矩阵；it为上盘的转动惯量矩阵：

参数i(t)的表达式如下，其中ix，iy与iz分别为x，y，z方向的转动惯量；rt为顶盘的半径：

减振机构的上层板和下层板的关系如下，其中ωmp(t)为上盘的角速度矩阵；ωmp(b)为下盘的角速度矩阵；为坐标转移矩阵：

假定下层板的在扰动力矩作用下的姿态角度为α,β,γ；则上层板的旋转角速度如下，公式中为坐标轴向量；为三个方向的角速度：

而上层板的旋转角速度表达式如下：

假定v＝[px,py,pz,α,β,γ]，带入式(10)到式(6),则整个系统的动能为：

其中：

式中m44＝ixcos²γ+iysin²γ；m55＝cos²α(ixcos²γ+iysin²γ)+izsin²α；m66＝iz；m45＝(ix-iy)cosαcosγsinγ；m56＝-izsinα；

将式(11)展开，可得下式，其中为z轴方向上的速度；为三个方向的角速度：

整个系统的虚功表达式如下，其中kh为横向弹簧的劲度系数；kv为纵向弹簧的劲度系数；为相应的虚位移；为相应的虚速度；czi,cαi,cβi,cγi为相应的阻尼系数：

其中(cα,cβ,cγ,cx,cy,cz)为减振系统在六个自由度上的阻尼参数，而(kh,kv)为横向和纵向的弹簧弹性系数；

公式(13)又可表示为：

其中

而对公式(12)积分可得：

由公式(16)-(17)可得：

其中，m44,m45,m55,m56,m66为相应的力矩；为相应方向的加速度；为三个方向的角速度；α,β,γ为相应的转动角度；为相应的虚位移；为相应的虚速度；ji(qi)为提出的参数；czi,cαi,cβi,cγi为相应的阻尼系数；

上式即为四自由度减振机构的动力学模型，假定在任何坐标轴有摄动存在，通过设计结构参数能够有效抑制外界干扰；当给定非线性结构参数，采用数值方法求解上方程，即可得到数值解；在给定不同的参数下，x构型的非线性阻尼不同，而减振效果也不同；在考虑减振系统的质量参数和体积参数的同时，优化减振机构的参数，可求解出最优减振机构；

步骤三：采用粒子群算法对被动抑制结构进行设计优化

在给定结构参数范围的条件下，考虑减振机构的质量参数约束、体积参数约束的条件下，以减振比例参数为优化目标，设计最优x结构非线性减振系统，具体步骤为：

step1：给定x构型的动力学方程如公式(18)所示；给定x构型的结构参数(cα,cβ,cγ,cx,cy,cz,kh,kv)的范围；

step2：给定优化方法的优化目标；

在此假定外界扰动、冲击振动在传递过减振机构后，传递的干扰引起的位置偏移、偏移速度以及偏移加速度分别为：c1,c2,c3为相应的系数，则优化指标为：

其中

式中(px,py,pz)分别为x轴y轴和z轴的位置偏移量；

step3：在给定外界干扰下，设定减振系统的减振比为优化目标，采用粒子群算法，在结构参数区间内，选定方程(18)中的结构参数值，求解最优减振结构参数。

本发明具有如下有益的技术效果：

本发明提供的一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法，该方法针对空间非合作目标抓捕过程中的碰撞冲击问题，通过在抓捕机构以及抓捕航天器之间装置仿生非线性结构方法实现抓捕被动减振，减小非合作目标抓捕过程中的冲击碰撞对基座的扰动；具体来说，针对空间非合作目标抓捕过程的冲击减振问题，考虑到仿生结构阻尼的质量和减振效果，对非线性结构参数进行优化设计，求解最优冲击碰撞扰动被动抑制结构设计。

附图说明

图1为x构型形变几何示意图。

具体实施方式

以下结合附图对本发明做出进一步的说明。

本发明提供的一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法，该方法针对空间非合作目标抓捕过程中的碰撞冲击问题，通过在抓捕机构以及抓捕航天器之间装置仿生非线性结构方法实现抓捕被动减振，减小非合作目标抓捕过程中的冲击碰撞对基座的扰动。

根据生物界普遍存在的x构型结构特征，设计的一种具有非线性阻尼特征的减振机构。该减振机构由多段x型结构组成，每个x型结构由两根杆组成。杆件之间的摩擦以及结构形变可有效减缓抓捕过程扰动碰撞冲击对抓捕航天器的影响。

针对空间非合作目标抓捕过程的冲击减振问题，考虑到仿生结构阻尼的质量和减振效果，对非线性结构参数进行优化设计，求解最优冲击碰撞扰动被动抑制结构设计。

空间非合作目标的主动被动混合控制方法的具体步骤包括：

步骤一：减振机构结构设计。

在生物中，x构型普遍存在于丹顶鹤、袋鼠、啄木鸟等生物的结构体中。x构型所具有的非线性阻尼能够有效减缓外界冲击以及振动对这些生物的影响。由此提炼出x构型减振机构，借助其x构型形变产生的非线性阻尼，可以有效减缓外界扰动给本体带来的冲击。

步骤二：减振机构动力学建模

设四自由度圆盘中“x”构型数目为n′，中间单自由度部分“x”构型数目为n。其中参数θ,qi,si,l的定义如图1所示，设四自由度圆盘中“x”构型数目为n′，中间单自由度部分“x”构型数目为n，θ为x构型与底面相连的部分与底盘之间的初始夹角；为x构型工作之后的位置与其初始位置之间的夹角；qi为x构型总体的纵向变化量；si为x构型总体的横向变化量；l为x构型每一部分的杆的长度，由图1的几何关系可得到如下关系式：

由公式(1)可求得：

而参数si可表示为

整个系统的动能表示如下，其中m为整个系统的质量；为z轴方向上的速度：

系统势能表示如下，其中ωmp(t)为上盘的角速度矩阵；it为上盘的转动惯量矩阵：

系统势能为：

参数i(t)的表达式如下，其中ix，iy与iz分别为x，y，z方向的转动惯量；rt为顶盘的半径：

减振机构的上层板和下层板的关系如下，其中ωmp(t)为上盘的角速度矩阵；ωmp(b)为下盘的角速度矩阵；为坐标转移矩阵：

假定下层板的在扰动力矩作用下的姿态角度为α,β,γ；则上层板的旋转角速度如下，公式中为坐标轴向量；为三个方向的角速度：

而上层板的旋转角速度为：

假定v＝[px,py,pz,α,β,γ]，带入式(10)到式(6),则整个系统的动能为：

其中：

式中m44＝ixcos²γ+iysin²γ；m55＝cos²α(ixcos²γ+iysin²γ)+izsin²α；m66＝iz；m45＝(ix-iy)cosαcosγsinγ；m56＝-izsinα；

将式(11)展开，可得下式，其中为z轴方向上的速度；为三个方向的角速度：

整个系统的虚功表达式如下，其中kh为横向弹簧的劲度系数；kv为纵向弹簧的劲度系数；为相应的虚位移；为相应的虚速度；czi,cαi,cβi,cγi为相应的阻尼系数：

其中(cα,cβ,cγ,cx,cy,cz)为减振系统在六个自由度上的阻尼参数，而(kh,kv)为横向和纵向的弹簧弹性系数；

公式(13)又可表示为：

其中

而对公式(12)积分可得：

由公式(16)-(17)可得：

其中，m44,m45,m55,m56,m66为相应的力矩；为相应方向的加速度；为三个方向的角速度；α,β,γ为相应的转动角度；为相应的虚位移；为相应的虚速度；ji(qi)为提出的参数；czi,cαi,cβi,cγi为相应的阻尼系数；

上式即为四自由度减振机构的动力学模型，假定在任何坐标轴有摄动存在，通过设计结构参数可以有效抑制外界干扰。当给定非线性结构参数，采用数值方法求解上方程，即可得到数值解。在给定不同的参数下，x构型的非线性阻尼不同，而减振效果也不同。在考虑减振系统的质量参数和体积参数的同时，优化减振机构的参数，可求解出最优减振机构。

步骤三：采用粒子群算法对被动抑制结构进行设计优化

在给定结构参数范围的条件下，考虑减振机构的质量参数约束、体积参数约束的条件下，以减振比例参数为优化目标，设计最优x结构非线性减振系统。具体步骤为：

step1：给定x构型的动力学方程如公式(18)所示；给定x构型的结构参数(cα,cβ,cγ,cx,cy,cz,kh,kv)的范围。

step2：给定优化方法的优化目标；

在此假定外界扰动、冲击振动在传递过减振机构后，传递的干扰引起的位置偏移、偏移速度以及偏移加速度分别为：c1,c2,c3为相应的系数，则优化指标为：

其中

式中(px,py,pz)分别为x轴y轴和z轴的位置偏移量。

step3：在给定外界干扰下，设定减振系统的减振比为优化目标。采用粒子群算法，在结构参数区间内，选定方程(18)中的结构参数值，求解最优减振结构参数。