[发明专利]一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法

权利要求书

1.一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法，其特征在于，包括以下步骤：

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

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

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

由公式(1)可求得：

而参数s_i可表示为

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

系统势能表示如下，其中Ω_mp(T)为上盘的角速度矩阵；I_T为上盘的转动惯量矩阵：

参数I_(T)的表达式如下，其中I_x，I_y与I_z分别为X，Y，Z方向的转动惯量；r_T为顶盘的半径：

减振机构的上层板和下层板的关系如下，其中Ω_mp(T)为上盘的角速度矩阵；Ω_mp(B)为下盘的角速度矩阵；为坐标转移矩阵：

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

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

假定V＝[P_x,P_y,P_z,α,β,γ]，带入式(10)到式(6),则整个系统的动能为：

其中：

式中M₄₄＝I_xcos²γ+I_ysin²γ；M₅₅＝cos²α(I_xcos²γ+I_ysin²γ)+I_zsin²α；M₆₆＝I_z；M₄₅＝(I_x-I_y)cosαcosγsinγ；M₅₆＝-I_zsinα；

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

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

其中(c_α,c_β,c_γ,c_x,c_y,c_z)为减振系统在六个自由度上的阻尼参数，而(k_h,k_v)为横向和纵向的弹簧弹性系数；

公式(13)又可表示为：

其中

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

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

其中，M₄₄,M₄₅,M₅₅,M₅₆,M₆₆为相应的力矩；为相应方向的加速度；为三个方向的角速度；α,β,γ为相应的转动角度；为相应的虚位移；为相应的虚速度；J_i(q_i)为提出的参数；c_zi,c_αi,c_βi,c_γi为相应的阻尼系数；

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

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

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

Step1：给定X构型的动力学方程如公式(18)所示；给定X构型的结构参数(c_α,c_β,c_γ,c_x,c_y,c_z,k_h,k_v)的范围；

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

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

其中

式中(P_x,P_y,P_z)分别为x轴y轴和z轴的位置偏移量；

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