CN109165414B - 一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法 - Google Patents
一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法 Download PDFInfo
- Publication number
- CN109165414B CN109165414B CN201810849950.XA CN201810849950A CN109165414B CN 109165414 B CN109165414 B CN 109165414B CN 201810849950 A CN201810849950 A CN 201810849950A CN 109165414 B CN109165414 B CN 109165414B
- Authority
- CN
- China
- Prior art keywords
- parameter
- vibration
- vibration damping
- follows
- equation
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Vibration Prevention Devices (AREA)
Abstract
本发明公开了一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法,该方法针对空间非合作目标抓捕过程中的碰撞冲击问题,通过在抓捕机构以及抓捕航天器之间装置仿生非线性结构方法实现抓捕被动减振,减小非合作目标抓捕过程中的冲击碰撞对基座的扰动;具体来说,针对空间非合作目标抓捕过程的冲击减振问题,考虑到仿生结构阻尼的质量和减振效果,对非线性结构参数进行优化设计,求解最优冲击碰撞扰动被动抑制结构设计。
Description
技术领域
本发明属于空间非合作目标抓捕控制领域,具体涉及一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法。
背景技术
随着空间开发利用活动的增加,空间碎片的数量急剧增多。而空间轨道的资源是有限的,如果不清理空间碎片,待其数量增大到一定程度后,已有碎片之间的碰撞会使产生碎片的速度大于清理的速度,从而引起“凯斯勒”效应,空间将不能再被利用。因此空间碎片的清理是亟待解决的任务。由于大量空间碎片具有旋转特性,在抓捕过程中,不可避免会出现碰撞冲击等。因此需要针对空间碎片抓捕过程中的振动冲击控制问题,提出高效可行的解决方法。
发明内容
本发明的目的在于提供一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法。
本发明采用如下技术方案来实现的:
一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法,该方法针对空间非合作目标抓捕过程中的碰撞冲击问题,通过在抓捕机构以及抓捕航天器之间装置仿生非线性结构方法实现抓捕被动减振,减小非合作目标抓捕过程中的冲击碰撞对基座的扰动;具体来说,针对空间非合作目标抓捕过程的冲击减振问题,考虑到仿生结构阻尼的质量和减振效果,对非线性结构参数进行优化设计,求解最优冲击碰撞扰动被动抑制结构设计。
本发明进一步的改进在于,根据生物界普遍存在的X构型结构特征,设计一种具有非线性阻尼特征的减振机构,该减振机构由多段X型结构组成,每个X型结构由两根杆组成;杆件之间的摩擦以及结构形变可有效减缓抓捕过程扰动碰撞冲击对抓捕航天器的影响。
本发明进一步的改进在于,包括以下步骤:
步骤一:减振机构结构设计
步骤二:减振机构动力学建模
设四自由度圆盘中“X”构型数目为n′,中间单自由度部分“X”构型数目为n,θ为X构型与底面相连的部分与底盘之间的初始夹角;为X构型工作之后的位置与其初始位置之间的夹角;qi为X构型总体的纵向变化量;si为X构型总体的横向变化量;L为X构型每一部分的杆的长度,由几何关系可得到如下关系式:
由公式(1)可求得:
而参数si可表示为
系统势能表示如下,其中Ωmp(T)为上盘的角速度矩阵;IT为上盘的转动惯量矩阵:
参数I(T)的表达式如下,其中Ix,Iy与Iz分别为X,Y,Z方向的转动惯量;rT为顶盘的半径:
而上层板的旋转角速度表达式如下:
假定V=[Px,Py,Pz,α,β,γ],带入式(10)到式(6),则整个系统的动能为:
其中:
式中M44=Ixcos2γ+Iysin2γ;M55=cos2α(Ixcos2γ+Iysin2γ)+Izsin2α;M66=Iz;M45=(Ix-Iy)cosαcosγsinγ;M56=-Izsinα;
其中(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:给定优化方法的优化目标;
其中
式中(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可表示为
系统势能表示如下,其中Ωmp(T)为上盘的角速度矩阵;IT为上盘的转动惯量矩阵:
系统势能为:
参数I(T)的表达式如下,其中Ix,Iy与Iz分别为X,Y,Z方向的转动惯量;rT为顶盘的半径:
而上层板的旋转角速度为:
假定V=[Px,Py,Pz,α,β,γ],带入式(10)到式(6),则整个系统的动能为:
其中:
式中M44=Ixcos2γ+Iysin2γ;M55=cos2α(Ixcos2γ+Iysin2γ)+Izsin2α;M66=Iz;M45=(Ix-Iy)cosαcosγsinγ;M56=-Izsinα;
其中(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:给定优化方法的优化目标;
其中
式中(Px,Py,Pz)分别为x轴y轴和z轴的位置偏移量。
Step3:在给定外界干扰下,设定减振系统的减振比为优化目标。采用粒子群算法,在结构参数区间内,选定方程(18)中的结构参数值,求解最优减振结构参数。
Claims (1)
1.一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法,其特征在于,包括以下步骤:
步骤一:减振机构结构设计,根据生物界普遍存在的X构型结构特征,设计一种具有非线性阻尼特征的减振机构,该减振机构由多段X型结构组成,每个X型结构由两根杆组成;杆件之间的摩擦以及结构形变可有效减缓抓捕过程扰动碰撞冲击对抓捕航天器的影响;
步骤二:减振机构动力学建模
设四自由度圆盘中“X”构型数目为n′,中间单自由度部分“X”构型数目为n,θ为X构型与底面相连的部分与底盘之间的初始夹角;为X构型工作之后的位置与其初始位置之间的夹角;qi为X构型总体的纵向变化量;si为X构型总体的横向变化量;L为X构型每一部分的杆的长度,由几何关系可得到如下关系式:
由公式(1)可求得:
而参数si可表示为
系统势能表示如下,其中Ωmp(T)为上盘的角速度矩阵;IT为上盘的转动惯量矩阵:
参数I(T)的表达式如下,其中Ix,Iy与Iz分别为X,Y,Z方向的转动惯量;rT为顶盘的半径:
而上层板的旋转角速度表达式如下:
假定V=[Px,Py,Pz,α,β,γ],带入式(10)到式(6),则整个系统的动能为:
其中:
式中M44=Ixcos2γ+Iysin2γ;M55=cos2α(Ixcos2γ+Iysin2γ)+Izsin2α;M66=Iz;M45=(Ix-Iy)cosαcosγsinγ;M56=-Izsinα;
其中(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:给定优化方法的优化目标;
其中
式中(Px,Py,Pz)分别为x轴y轴和z轴的位置偏移量;
Step3:在给定外界干扰下,设定减振系统的减振比为优化目标,采用粒子群算法,在结构参数区间内,选定方程(18)中的结构参数值,求解最优减振结构参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810849950.XA CN109165414B (zh) | 2018-07-28 | 2018-07-28 | 一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810849950.XA CN109165414B (zh) | 2018-07-28 | 2018-07-28 | 一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109165414A CN109165414A (zh) | 2019-01-08 |
CN109165414B true CN109165414B (zh) | 2022-10-04 |
Family
ID=64898603
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810849950.XA Active CN109165414B (zh) | 2018-07-28 | 2018-07-28 | 一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109165414B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111077896B (zh) * | 2019-12-27 | 2022-07-19 | 南京理工大学 | 一种基于改进分层算法的充液挠性航天器参数优化方法 |
CN111814277B (zh) * | 2020-08-24 | 2022-03-15 | 西北工业大学 | 一种隔振平台设备及其动力学模型构建方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103863583A (zh) * | 2014-02-27 | 2014-06-18 | 中国空间技术研究院 | 一种飞舌抓捕机构及空间目标捕获方法 |
US9663252B1 (en) * | 2016-12-07 | 2017-05-30 | Beihang University | Method for attitude controlling based on finite time friction estimation for flexible spacecraft |
CN107520844A (zh) * | 2017-09-21 | 2017-12-29 | 西北工业大学 | 一种空间机械臂抓捕非合作目标的多面体碰撞动力学分析方法 |
-
2018
- 2018-07-28 CN CN201810849950.XA patent/CN109165414B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103863583A (zh) * | 2014-02-27 | 2014-06-18 | 中国空间技术研究院 | 一种飞舌抓捕机构及空间目标捕获方法 |
US9663252B1 (en) * | 2016-12-07 | 2017-05-30 | Beihang University | Method for attitude controlling based on finite time friction estimation for flexible spacecraft |
CN107520844A (zh) * | 2017-09-21 | 2017-12-29 | 西北工业大学 | 一种空间机械臂抓捕非合作目标的多面体碰撞动力学分析方法 |
Non-Patent Citations (2)
Title |
---|
空间微重力环境地面模拟系统的控制器设计;陈三风等;《机器人》;20080515(第03期);全文 * |
空间机器人目标捕获过程的碰撞冲击分析;谢冉等;《计算机仿真》;20160215(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109165414A (zh) | 2019-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104267732B (zh) | 基于频域分析的挠性卫星高稳定度姿态控制方法 | |
CN109165414B (zh) | 一种空间非合作目标抓捕过程中振动冲击抑制机构优化设计方法 | |
Zhang et al. | Attitude control of a quadrotor aircraft subject to a class of time-varying disturbances | |
CN102778891B (zh) | 一种星上控制力矩陀螺群隔振平台的参数选择方法 | |
Elias | Seismic energy assessment of buildings with tuned vibration absorbers | |
CN110466808B (zh) | 基于多等碰撞概率线法的凸多边形航天器安全控制方法 | |
Pfeiffer et al. | Contacts in multibody systems | |
Zappulla et al. | Experiments on autonomous spacecraft rendezvous and docking using an adaptive artificial potential field approach | |
CN108628345B (zh) | 一种电磁航天器编队悬停协同控制方法及系统 | |
Chu et al. | Modeling and stabilization control for space-borne series-wound capturing mechanism with multi-stage damping | |
CN112276952B (zh) | 一种多机器人系统的鲁棒同时镇定方法及系统 | |
Xie et al. | Acceleration-level trajectory planning for a dual-arm space robot | |
CN106096206A (zh) | 一种柔性航天器分布式执行机构和敏感器优化配置方法 | |
Fu et al. | Multi-UAV formation control method based on modified artificial physics | |
CN112389679B (zh) | 一种考虑多种摄动力的小卫星常值推力轨道递推方法 | |
CN111783285B (zh) | 多点支撑结构的载荷传递路径优化方法 | |
CN109188900A (zh) | 一种外界扰动引起的航天器柔性部件振动被动抑制方法 | |
Chen et al. | Post-capture angular momentum management of space robot with controllable damping joints | |
CN111814277A (zh) | 一种隔振平台及其动力学模型构建方法 | |
CN109987258A (zh) | 一种空间机器人捕获非合作目标后的消旋方法 | |
NGUYEN | An analytical solution of dynamic vibration absorber to suppress the vibration of a pendulum structure subjected to moving loads caused by a hanging point | |
Zhong et al. | Minimum time-jerk trajectory generation for a mobile articulated manipulator | |
Wang et al. | Multi-Objective Path Planning Based on Parallel Elite Particle Swarm Optimization Algorithm | |
CN117666604A (zh) | 一种基于分布式的多无人机编队重构控制方法 | |
Minowa et al. | Semi-active vibration control of double-layer cylindrical lattice shell with variable damping mechanism |
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 |