CN106202688A - 六轴隔振平台的动力学完全各向同性设计方法 - Google Patents

六轴隔振平台的动力学完全各向同性设计方法 Download PDF

Info

Publication number
CN106202688A
CN106202688A CN201610523893.7A CN201610523893A CN106202688A CN 106202688 A CN106202688 A CN 106202688A CN 201610523893 A CN201610523893 A CN 201610523893A CN 106202688 A CN106202688 A CN 106202688A
Authority
CN
China
Prior art keywords
platform
vibration
mounting plate
upper mounting
kinetics
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.)
Granted
Application number
CN201610523893.7A
Other languages
English (en)
Other versions
CN106202688B (zh
Inventor
杨小龙
任薇
吴洪涛
李耀
陈柏
刘艳梨
程世利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201610523893.7A priority Critical patent/CN106202688B/zh
Publication of CN106202688A publication Critical patent/CN106202688A/zh
Application granted granted Critical
Publication of CN106202688B publication Critical patent/CN106202688B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Landscapes

  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Vibration Prevention Devices (AREA)

Abstract

本发明公开一种六轴隔振平台的动力学完全各向同性设计方法,给出了动力学完全各向同性条件,包括:一个几何协调方程、两个惯量方程和三个质量几何约束方程。根据给出的完全各向同性条件设计隔振平台的几何、质/惯量参数,可以使隔振平台六个非零自然频率全部相等,在任意方向上获得相同的隔振性能。

Description

六轴隔振平台的动力学完全各向同性设计方法
技术领域
本发明属于机械装备的动力学、优化设计研究领域,尤其是基于并联机构的隔振平台的动力学性能优化设计。
背景技术
高精度空间系统的精确指向需要超静环境,空间干涉仪、激光通信设备等高精度航天器需达到纳米级的运动稳定性。然而,航天器在轨运行期间,星上反作用飞轮、控制力矩陀螺、低温制冷机、太阳电池阵驱动机构等部件的运动将产生频带较宽的低幅机械振动,称为航天器微振动。为了使航天器在六个自由度上达到极高的光学稳定水平,可以设计一种基于Stewart机构的六轴隔振平台。它由上、下两个平台和六个并联的腿组成,每条腿形成一个单轴隔振器,通过两个球铰与上、下两平台相连。近二十年已有不少研究者开展了这方面的研究。Geng和Haynes最早提出了一种基于特殊构型Stewart机构的六轴隔振方法和系统,这种特殊构型称为立方Stewart平台,成为设计六轴隔振系统的典型形式。采用立方构型的隔振平台的优势是能够建立简化的解耦动力学模型,每条支腿对平台的作用力相对独立,可以采用分散控制策略实现主动隔振。之后,Spanos和Rahman,Hauge和Campbell以及Preumont等也基于立方Stewart平台研究了六轴隔振系统的动力学建模和主动控制。然而,上述隔振平台的六个非零自然频率不同,而每条腿机械结构、作动器、传感器和控制增益是相同的,这就不能针对所有的模态得到一致的隔振性能。因此,隔振平台的设计应使其最大自然频率与最小自然频率的比值最小。这个问题称为隔振平台动力学各向同性设计,目的是使其六个非零自然频率相等,以便在任意方向上获得相同的隔振性能。在动力学各向同性设计方面,Ma和Angeles定义了基于广义质量矩阵的动力学各向同性指标,并以Stewart平台为例进行了动力学各向同性设计研究,得到的最优指标值(等价于最大最小频率比)为1.1242。Jiang和He等人针对Stewart平台提出了基于自然频率矩阵的动力学各向同性设计,得到的最优指标值(定义为最大最小频率比)为1.2508。基于自然频率矩阵的动力学各向同性指标的优化会使得系统的各阶自然频率更加接近,当各向同性指标为最优(值为1)时,系统的各阶自然频率相同,称之为动力学完全各向同性。目前,动力学完全各向同性设计问题仍然没有解决。
故,需要一种新的技术方案以解决上述问题。
发明内容
本发明的目的是针对现有技术存在的不足,提供一种六轴隔振平台的动力学完全各向同性设计方法。
为解决上述问题,本发明的六轴隔振平台的动力学完全各向同性设计可采用如下技术方案:
1、一种六轴隔振平台的动力学完全各向同性设计方法,所述的六轴隔振平台包括上平台、下平台和连接上平台及下平台的六条腿,每一条腿是一个单轴隔振器,每一条腿的两端用球铰分别与上平台、下平台相连;其中,与上平台连接的球铰为上球铰;与下平台连接的球铰为下球铰;其特征在于,该方法包括如下步骤:
(1)确定动力学各向同性设计变量
连接上平台的六个上球铰安装在一个平面内,且相邻两个上球铰点的夹角为α;连接下平台的六个下球铰同样安装在一个平面内,相邻两个下球铰点的夹角为β;定义夹角差γ=(α-β)/2,上平台半径为ra,下平台半径为rb,每条腿上安装一个作动器,作动器的等效刚度为k;六轴隔振平台处于平衡位置时高度为h;定义上平台的质量ma,上平台相对于质心的转动惯量为Ia=diag(Iax Iay Iaz)3×3,下平台的质量是mb,下平台相对于质心的转动惯量为Ib=diag(Ibx Iby Ibz)3×3;六轴隔振平台的设计变量为上、下平台的质/惯量参数和隔振平台的机构几何参数,共十二个,满足如下约束条件:
h>0,ζ>0,μ>0,0<γ<π/2
ma>0,mb>0,Iax>0,Iay>0,Iaz>0,Ibx>0,Iby>0,Ibz>0
上式中,ζ=ra/h,μ=rb/h,ζ为上平台的半径和平台高度比;μ为下平台的半径和平台高度比;
(2)按照以下动力学各向同性条件确定设计参数值
一个几何协调方程:
c o s γ = ζ 2 μ 2 + 1 - 1 ζ μ ;
两个惯量方程:
I a x = I a y I b x = I b y ;
三个质量几何约束方程:
( 1 m a + 1 m b ) / ( 1 I a z + 1 I b z ) = 2 h 2 ( 1 + ζ 2 μ 2 - 1 ) 1 I a x / ( 1 I a z + 1 I b z ) = 2 + 2 μ 2 - 2 1 + ζ 2 μ 2 1 I b x / ( 1 I a z + 1 I b z ) = 2 + 2 ζ 2 - 2 1 + ζ 2 μ 2 .
附图说明
图1是本发明中六轴隔振平台的几何描述图。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
其中,如图1所示,本发明中所述的六轴隔振平台包括上平台(载荷平台)、下平台(基平台)和连接上平台及下平台的六条腿,每一条腿是一个单轴隔振器,设计有作动器和传感器,每一条腿的两端用球铰分别与上平台、下平台相连;其中,与上平台连接的球铰为上球铰;与下平台连接的球铰为下球铰。而针对该六轴隔振平台包,本发明提供的六轴隔振平台的动力学完全各向同性设计方法包括如下步骤。
1求解六轴隔振平台的自然频率
连接上平台的六个上球铰安装在一个平面内,且相邻两个上球铰点的夹角为α;连接下平台的六个下球铰同样安装在一个平面内,相邻两个下球铰点的夹角为β;α和β按圆周循环对称分布。定义夹角差γ=(α-β)/2,上平台半径为ra,下平台半径为rb,每条腿上安装一个压电作动器,作动器的等效刚度为k。隔振平台处于平衡位置时高度为h。定义上平台的质量ma,上平台相对于质心的转动惯量为Ia=diag(Iax Iay Iaz)3×3,下平台的质量是mb,下平台相对于质心的转动惯量为Ib=diag(Ibx Iby Ibz)3×3
建立三个笛卡尔直角坐标系,分别为惯性系G-xyz,下平台质心坐标系O-xyz和上平台质心坐标系O'-xyz。O'-xyz和O-xyz相对于G-xyz的姿态分别用转动矩阵Ra和Rb表示。pa0和pb0分别表示上、下平台的初始位置,pa和pb分别表示上、下平台的位移。每支腿与上平台的连接点在O'-xyz中的矢径记为:a1...a6,与下平台的连接点在O-xyz中的矢径记为:b1...b6。ai和bi在G-xyz中用qai和qbi表示,那么存在转换关系:qai=Raai,qbi=Rbbi。仅考虑隔振平台的平衡位置处的自然频率,那么Ra=Rb=E,E是3×3单位矩阵。
隔振平台的位置方程为
liei=pa+pa0+qai-pb-pb0-qbi(i=1,2,…,6) (1)
li是每条腿的长度,ei是沿着每条腿的单位矢量,表示腿的方向。
隔振平台的质量矩阵为:
M=diag(maE Ia mbE Ib)12×12 (2)
隔振平台的刚度矩阵为:
K=JTdiag(k)6×6J (3)
J是Stewart机构雅可比矩阵:
J l u = e 1 T ( q a 1 × e 1 ) T - e 1 T - ( q b 1 × e 1 ) T . . . . . . . . . . . . e 6 T ( q a 6 × e 6 ) T - e 6 T - ( q b 6 × e 6 ) T - - - ( 4 )
最后求解K关于M的广义特征值,即自然频率的平方:
ω 1 2 = ω 2 2 = ω 3 2 = ω 4 2 = ω 5 2 = ω 6 2 = 0
ω 7 2 = 6 kh 2 l 0 2 ( 1 m a + 1 m b ) , ω 8 2 = 3 kr a 2 r b 2 ( 1 - cos 2 γ ) l 0 2 ( 1 I a z + 1 I b z ) - - - ( 5 )
ω 9 2 = kh 2 ( A + B ) / l 0 2 , ω 10 2 = kh 2 ( A - B ) / l 0 2
ω 11 2 = kh 2 ( C + D ) / l 0 2 , ω 12 2 = kh 2 ( C - D ) / l 0 2
上式中,
A = 3 2 h 2 I a x I b x m a m b ( I a x I b x ( l 0 2 - h 2 ) ( m a + m b ) + h 2 m a m b ( I b x r a 2 + I a x r b 2 ) ) B = 3 2 h 2 I a x I b x m a m b ( I a x I b x ( l 0 2 - h 2 ) ( m a + m b ) + h 2 m a m b ( I b x r a 2 + I a x r b 2 ) ) 2 - 4 h 2 I a x I b x m a m b ( h 2 m a m b + ( I a x + I b x ) ( m a + m b ) ) r a 2 r b 2 sin 2 γ C = 3 2 h 2 I a y I b y m a m b ( I a y I b y ( l 0 2 - h 2 ) ( m a + m b ) + h 2 m a m b ( I b y r a 2 + I a y r b 2 ) ) D = 3 2 h 2 I a y I b y m a m b ( I a y I b y ( l 0 2 - h 2 ) ( m a + m b ) + h 2 m a m b ( I b y r a 2 + I a y r b 2 ) ) 2 - 4 h 2 I a y I b y m a m b ( h 2 m a m b + ( I a y + I b y ) ( m a + m b ) ) r a 2 r b 2 sin 2 γ - - - ( 6 )
六个为零的自然频率表示隔振平台的自由漂浮状态。分别是两个二次代数方程组的根。动力学完全各向同性设计的目标是使得ω7,…,ω12相同。
2动力学完全各向同性条件
隔振平台的设计变量确定为上、下平台的质/惯量参数和隔振平台的机构几何参数,共十二个,满足如下约束条件:
h > 0 , &zeta; > 0 , &mu; > 0 , 0 < &gamma; < &pi; / 2 m a > 0 , m b > 0 , I a x > 0 , I a y > 0 , I a z > 0 , I b x > 0 , I b y > 0 , I b z > 0 - - - ( 7 )
上式中,ζ=ra/h,μ=rb/h。
令隔振平台的六阶非零自然频率相等:ω7=ω8=…=ω12,可以获得隔振平台获得动力学各向同性条件,包括
一个几何协调方程:
c o s &gamma; = &zeta; 2 &mu; 2 + 1 - 1 &zeta; &mu; - - - ( 8 )
两个惯量方程:
I a x = I a y I b x = I b y - - - ( 9 )
三个质量几何约束方程:
( 1 m 2 + 1 m b ) / ( 1 I a z + 1 I b z ) = 2 h 2 ( 1 + &zeta; 2 &mu; 2 - 1 ) 1 I a x / ( 1 I a z + 1 I b z ) = 2 + 2 &mu; 2 - 2 1 + &zeta; 2 &mu; 2 1 I b x / ( 1 I a z + 1 I b z ) = 2 + 2 &zeta; 2 - 2 1 + &zeta; 2 &mu; 2 - - - ( 10 )
3应用示范
在进行隔振平台动力学各向同性设计时,根据(8)、(9)和(10)共6个约束方程,无法唯一地确定12个设计变量,此时可以根据其他设计条件确定一部分参数后再求解其他参数,完成设计。这里将给出已知质/惯量参数情形下的应用示范。
当质/惯量参数确定后,由方程(10)可以确定h,ζ,μ,γ再由方程(8)解出。
表1给出了三组设计条件,表2是按照各向同性条件得到的三组设计结果,数据均采用国际单位制(SI)。尽管表2中的γ大于π/3,在结构设计时可以改变腿的设计,使相邻两腿不相交。
表1三组设计实例的质/惯量参数
表2三组设计实例的结果

Claims (1)

1.一种六轴隔振平台的动力学完全各向同性设计方法,所述的六轴隔振平台包括上平台、下平台和连接上平台及下平台的六条腿,每一条腿是一个单轴隔振器,每一条腿的两端用球铰分别与上平台、下平台相连;其中,与上平台连接的球铰为上球铰;与下平台连接的球铰为下球铰;其特征在于,该方法包括如下步骤:
(1)确定动力学各向同性设计变量
连接上平台的六个上球铰安装在一个平面内,且相邻两个上球铰点的夹角为α;连接下平台的六个下球铰同样安装在一个平面内,相邻两个下球铰点的夹角为β;定义夹角差γ=(α-β)/2,上平台半径为ra,下平台半径为rb,每条腿上安装一个作动器,作动器的等效刚度为k;六轴隔振平台处于平衡位置时高度为h;定义上平台的质量ma,上平台相对于质心的转动惯量为Ia=diag(Iax Iay Iaz)3×3,下平台的质量是mb,下平台相对于质心的转动惯量为Ib=diag(Ibx Iby Ibz)3×3;六轴隔振平台的设计变量为上、下平台的质/惯量参数和隔振平台的机构几何参数,共十二个,满足如下约束条件:
h>0,ζ>0,μ>0,0<γ<π/2
ma>0,mb>0,Iax>0,Iay>0,Iaz>0,Ibx>0,Iby>0,Ibz>0
上式中,ζ=ra/h,μ=rb/h;ζ为上平台的半径和平台高度比;μ为下平台的半径和平台高度比;
(2)按照以下动力学各向同性条件确定设计参数值
一个几何协调方程:
c o s &gamma; = &zeta; 2 &mu; 2 + 1 - 1 &zeta; &mu; ;
两个惯量方程:
I a x = I a y I b x = I b y ;
三个质量几何约束方程:
( 1 m a + 1 m b ) / ( 1 I a z + 1 I b z ) = 2 h 2 ( 1 + &zeta; 2 &mu; 2 - 1 ) 1 I a x / ( 1 I a z + 1 I b z ) = 2 + 2 &mu; 2 - 2 1 + &zeta; 2 &mu; 2 1 I b x / ( 1 I a z + 1 I b z ) = 2 + 2 &zeta; 2 - 2 1 + &zeta; 2 &mu; 2 .
CN201610523893.7A 2016-07-04 2016-07-04 六轴隔振平台的动力学完全各向同性设计方法 Active CN106202688B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610523893.7A CN106202688B (zh) 2016-07-04 2016-07-04 六轴隔振平台的动力学完全各向同性设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610523893.7A CN106202688B (zh) 2016-07-04 2016-07-04 六轴隔振平台的动力学完全各向同性设计方法

Publications (2)

Publication Number Publication Date
CN106202688A true CN106202688A (zh) 2016-12-07
CN106202688B CN106202688B (zh) 2019-01-25

Family

ID=57466248

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610523893.7A Active CN106202688B (zh) 2016-07-04 2016-07-04 六轴隔振平台的动力学完全各向同性设计方法

Country Status (1)

Country Link
CN (1) CN106202688B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108897239A (zh) * 2018-06-29 2018-11-27 北京控制工程研究所 一种航天器两级姿态控制模拟系统
CN109555812A (zh) * 2018-07-18 2019-04-02 中国航空工业集团公司洛阳电光设备研究所 一种压电驱动的各向同性多自由度隔振平台

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
佟志忠 等: "基于加速度各向同性的空间光学仪器主动隔振Stewart平台设计", 《振动与冲击》 *
吴迎: "基于Stewart平台的卫星微振动主动控制方法研究及装置优化设计", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *
陈登硕: "基于电磁作动器的Stewart隔振平台振动主动控制研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108897239A (zh) * 2018-06-29 2018-11-27 北京控制工程研究所 一种航天器两级姿态控制模拟系统
CN108897239B (zh) * 2018-06-29 2021-03-26 北京控制工程研究所 一种航天器两级姿态控制模拟系统
CN109555812A (zh) * 2018-07-18 2019-04-02 中国航空工业集团公司洛阳电光设备研究所 一种压电驱动的各向同性多自由度隔振平台

Also Published As

Publication number Publication date
CN106202688B (zh) 2019-01-25

Similar Documents

Publication Publication Date Title
CN106094528B (zh) 一种空间柔性机械臂振动抑制算法
Koop Global gravity field modelling using satellite gravity gradiometry
CN111099045B (zh) 双超卫星动力学与控制气浮平台全物理仿真方法
Chen et al. Dynamical modelling and control of space tethers: a review of space tether research
CN106950853B (zh) 一种球形贮箱微重力环境下液体晃动的建模方法
CN110329550B (zh) 用于敏捷卫星应用的姿态控制
CN103807353B (zh) 无角位移隔振装置及其支链
Hu et al. Dynamics and vibration suppression of space structures with control moment gyroscopes
WO2019191294A1 (en) Methods and apparatus for in-situ measurements of atmospheric density
Lin et al. Recent advances in precision measurement & pointing control of spacecraft
CN104238563B (zh) 可变面倾角的控制力矩陀螺群设计方法
Gutierrez Performance assessment and enhancement of precision controlled structures during conceptual design
CN106202688A (zh) 六轴隔振平台的动力学完全各向同性设计方法
CN109683628A (zh) 一种基于有限时间分布式速度观测器的航天器相对位置控制方法
CN105912020B (zh) 一种弱稳定边界区域内航天器编队飞行方法
Bainum et al. Motion and stability of a dual-spin satellite with nutation damping
Robb et al. Distributed magnetic attitude control for large space structures
CN107270905A (zh) 一种静基座下惯性平台连续翻滚自标定自对准方法
Wasfy et al. Multibody dynamic simulation of the next generation space telescope using finite elements and fuzzy sets
CN111605733B (zh) 航天器自治协同粗精分层主被一体三超控制参数确定方法
Khoroshylov et al. Dynamics and attitude control of space-based synthetic aperture radar
CN104484498B (zh) 确定有效推进加速度折损的太阳帆瞬态动力学分析方法
Dodds et al. Sliding-mode control system for the three-axis attitude control of rigid-body spacecraft with unknown dynamics parameters
CN108614578B (zh) 一种小推力悬浮轨道上航天器编队飞行方法
Chen et al. Balancing control of a cubical robot balancing on its corner

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant