CN106289641B - 挠性航天器质心位置和转动惯量参数联合辨识方法 - Google Patents

挠性航天器质心位置和转动惯量参数联合辨识方法 Download PDF

Info

Publication number
CN106289641B
CN106289641B CN201610784812.9A CN201610784812A CN106289641B CN 106289641 B CN106289641 B CN 106289641B CN 201610784812 A CN201610784812 A CN 201610784812A CN 106289641 B CN106289641 B CN 106289641B
Authority
CN
China
Prior art keywords
spacecraft
flexible
identification
equation
matrix
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
Application number
CN201610784812.9A
Other languages
English (en)
Other versions
CN106289641A (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.)
Shanghai Aerospace Control Technology Institute
Original Assignee
Shanghai Aerospace Control Technology Institute
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 Shanghai Aerospace Control Technology Institute filed Critical Shanghai Aerospace Control Technology Institute
Priority to CN201610784812.9A priority Critical patent/CN106289641B/zh
Publication of CN106289641A publication Critical patent/CN106289641A/zh
Application granted granted Critical
Publication of CN106289641B publication Critical patent/CN106289641B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M1/00—Testing static or dynamic balance of machines or structures
    • G01M1/10—Determining the moment of inertia
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M1/00—Testing static or dynamic balance of machines or structures
    • G01M1/12—Static balancing; Determining position of centre of gravity
    • G01M1/122—Determining position of centre of gravity
    • G01M1/125—Determining position of centre of gravity of aircraft

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Navigation (AREA)

Abstract

本发明涉及一种挠性航天器质心位置和转动惯量参数联合辨识方法，其中将航天器动力学方程转化为辨识观测方程，对航天器质量特性参数进行辨识。本发明提出的方法将柔性附件的挠性振动干扰力矩引入辨识方程，增加辨识观测方程中观测量信息，进而提高柔性航天器的质心位置和转动惯量参数辨识精度。与现有技术相比，解决了柔性航天器质心位置和转动惯量参数高精度联合辨识问题，可广泛应用于带大柔性附件卫星的质心位置和转动惯量参数的辨识。

Description

挠性航天器质心位置和转动惯量参数联合辨识方法
技术领域
本发明涉及航天器的动力学参数辨识技术，尤其是带大型柔性附件航天器的质心位置和转动惯量的联合辨识方法。
背景技术
由空间环境和任务的变化带来不确定因素，以及模型的建立误差致使航天器的质量特性参数不能通过地面实验手段计算得到，只能通过在轨辨识技术获取。柔性航天器由于安装了大型柔性附件，而航天器在进行姿态调整的过程中易激发柔性附件的挠性振动，并对航天器的中心刚体产生耦合挠性干扰力矩，直接影响挠性航天器质量特性参数辨识精度，不利于挠性航天器的高精度姿态控制实现。对于柔性航天器的质心位置和转动惯量参数辨识，随着柔性附件与卫星平台耦合作用越强，产生的挠性振动干扰力矩对卫星姿态的影响作用越强，在质量特性参数辨识观测方程的测量中若仅考虑执行机构的控制作用，则无法实现对质量特性参数辨识的有效性和准确性。
发明内容
本发明的目的在于提供一种挠性航天器的质心卫星和转动惯量参数的联合辨识方法，实现对挠性航天器动力学模型参数修正。
本发明的技术方案是提供一种挠性航天器质心位置和转动惯量参数联合辨识方法：
定义
Jc＝[Jx,Jy,Jz,Jyz,Jxz,Jxy]T
Jx、Jy、Jz为航天器的转动惯量矩阵J中的主惯量；Jxy、Jxz、Jyz为转动惯量矩阵J中的惯量积；ΔR为推力器在本体坐标系中的位置误差矢量；
Zi为第i次测量的测量方程
Zi＝HiX+Vi i＝1,2,…,k (10)
其中，在将航天器动力学方程转化为辨识观测方程后可以求得，Vi为零均值的高斯白噪声；
则前k+1次测量为
其中，
Zk+1为第k+1次测量的测量方程
Zk+1＝Hk+1X+Vk+1 (12)
由前k次测量确定的加权最小二乘估计为
其中，
W1，W2，…，Wk为加权矩阵；
令
基于递推最小二乘估计算法进行质量特性参数辨识：
定义柔性航天器的本体坐标系为O-XbYbZb，质心在本体坐标系中的位置矢量为Ro，推力器在本体坐标系中的位置矢量为Rc；卫星本体相对于惯性坐标系的转速为ω，各轴分别为ωx、ωy、ωz，柔性附件的模态坐标向量为η，模态阻尼比矩阵为ξ，模态刚度矩阵为Λ，耦合矩阵为B；
在对挠性航天器的质心位置和转动惯量参数进行联合辨识之前，将航天器动力学方程转化为辨识观测方程：
挠性航天器的动力学方程为
其中，F为推力器推力矢量；
航天器的转动惯量矩阵J为
转化得到航天器的辨识观测方程为
其中，为F的斜对称矩阵；
为了解决航天器姿态机动过程中，柔性附件产生挠性振动耦合干扰力矩，直接影响航天器的姿态控制性能的问题，本发明提出的方法将柔性附件的挠性振动干扰力矩引入辨识方程，增加辨识观测方程中观测量信息，进而提高柔性航天器的质心位置和转动惯量参数辨识精度。与现有技术相比，解决了柔性航天器质心位置和转动惯量参数高精度联合辨识问题，可广泛应用于带大柔性附件卫星的质心位置和转动惯量参数的辨识。
附图说明
图1为挠性航天器结构示意图；
图2为本发明所述方法的流程图。
具体实施方式
本发明对挠性航天器的质心位置和转动惯量参数进行辨识，实现对挠性航天器的动力学模型参数修正，主要包括建立观测方程，进行参数辨识，将柔性附件的挠性振动干扰力矩引入观测方程，保证观测量信息全面，实现挠性航天器质量特性参数的高精度辨识。图1为挠性航天器结构示意图。
本发明的挠性航天器质心位置和转动惯量参数辨识方法，定义柔性航天器的本体坐标系为O-XbYbZb，质心在本体坐标系中的位置矢量为Ro，推力器在本体坐标系中的位置矢量为Rc，航天器的转动惯量矩阵为J，其中主惯量为Jx、Jy、Jz，惯量积为Jxy、Jxz、Jyz，卫星本体相对于惯性坐标系的转速为ω，各轴分别为ωx、ωy、ωz，N个推力器的推力矢量为fi(i＝1,2,…,N)，柔性附件的模态坐标向量为η，模态阻尼比矩阵为ξ，模态刚度矩阵为Λ，耦合矩阵为B。
如图2所示，所述挠性航天器质心位置和转动惯量参数联合辨识方法，包括如下步骤：
1)将航天器动力学方程转化为辨识观测方程
航天器的动力学方程为
其中，和分别是η的一阶导数和二阶导数，为ω的一阶导数，ΔR为推力器在本体坐标系中的位置误差矢量，F为推力器推力矢量。
航天器的辨识方程为
其中，为F的斜对称矩阵。
Jc＝[Jx,Jy,Jz,Jyz,Jxz,Jxy]T
2)对挠性航天器的质心位置和转动惯量参数进行辨识：
定义
根据公式(18)，令Zi为第i次测量，测量方程为
Zi＝HiX+Vi i＝1,2,…,k (19)
其中，Vi为零均值的高斯白噪声。
则前k+1次测量为
其中，
Zk+1为第k+1次测量，测量方程为
Zk+1＝Hk+1X+Vk+1 (21)
由前k次测量确定的加权最小二乘估计为
其中，W1，W2，…，Wk为加权矩阵。
令
则
由前k+1次测量确定的加权最小二乘估计为
其中，
则基于递推最小二乘估计算法进行质心位置和转动惯量参数辨识算法为
本发明采用的方法与现有技术相比，其优点和有益效果是：解决了挠性航天器的质心位置和转动惯量的参数辨识问题，实现了对挠性航天器的动力学模型参数修正。

Claims (1)

1.一种挠性航天器质心位置和转动惯量参数联合辨识方法，其特征在于，
对挠性航天器的质心位置和转动惯量参数进行联合辨识；
定义
Jc＝[Jx,Jy,Jz,Jyz,Jxz,Jxy]T
Jx、Jy、Jz为航天器的转动惯量矩阵J中的主惯量；Jxy、Jxz、Jyz为转动惯量矩阵J中的惯量积；ΔR为推力器在本体坐标系中的位置误差矢量；
Zi为第i次测量的测量方程
Zi＝HiX+Vi i＝1,2,…,k (1)
其中，Hi是通过将航天器动力学方程转化为辨识观测方程后求取的参数，Vi为零均值的高斯白噪声；
则前k+1次测量为
其中，
Zk+1为第k+1次测量的测量方程
Zk+1＝Hk+1X+Vk+1 (3)
由前k次测量确定的加权最小二乘估计为
其中，
W1，W2，…，Wk为加权矩阵；
令
基于递推最小二乘估计算法进行质量特性参数辨识：
定义柔性航天器的本体坐标系为O-XbYbZb，质心在本体坐标系中的位置矢量为Ro，推力器在本体坐标系中的位置矢量为Rc；卫星本体相对于惯性坐标系的转速为ω，各轴分别为ωx、ωy、ωz，柔性附件的模态坐标向量为η，模态阻尼比矩阵为ξ，模态刚度矩阵为Λ，耦合矩阵为B；
在对挠性航天器的质心位置和转动惯量参数进行联合辨识之前，将航天器动力学方程转化为辨识观测方程：
挠性航天器的动力学方程为
其中，F为推力器推力矢量；
航天器的转动惯量矩阵J为
转化得到航天器的辨识观测方程为
其中，为F的斜对称矩阵；
则，第i次测量的测量方程中的参数
CN201610784812.9A 2016-08-31 2016-08-31 挠性航天器质心位置和转动惯量参数联合辨识方法 Active CN106289641B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610784812.9A CN106289641B (zh) 2016-08-31 2016-08-31 挠性航天器质心位置和转动惯量参数联合辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610784812.9A CN106289641B (zh) 2016-08-31 2016-08-31 挠性航天器质心位置和转动惯量参数联合辨识方法

Publications (2)

Publication Number Publication Date
CN106289641A CN106289641A (zh) 2017-01-04
CN106289641B true CN106289641B (zh) 2018-12-25

Family

ID=57674111

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610784812.9A Active CN106289641B (zh) 2016-08-31 2016-08-31 挠性航天器质心位置和转动惯量参数联合辨识方法

Country Status (1)

Country Link
CN (1) CN106289641B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107389266B (zh) * 2017-08-31 2019-08-23 中国航空工业集团公司沈阳飞机设计研究所 一种飞行器重心柔性测量方法
CN107703748B (zh) * 2017-10-09 2020-01-07 东南大学 一种基于偏置板设计的重载机器人静态刚度辨识方法
CN109870271A (zh) * 2019-01-31 2019-06-11 西北工业大学 大尺度挠性航天器转动惯量卡尔曼滤波辨识方法
CN109870273A (zh) * 2019-01-31 2019-06-11 西北工业大学 基于动量守恒的航天器在轨质心辨识方法
CN109974933A (zh) * 2019-01-31 2019-07-05 西北工业大学 考虑卫星挠性耦合的转动惯量递推最小二乘辨识方法
CN111707414B (zh) * 2020-06-23 2022-05-27 重庆大学 地磅式大型重型飞行器质量质心测量装置及方法
CN111707412B (zh) * 2020-07-10 2022-06-07 哈尔滨工业大学 一种大型三轴气浮台转动惯量在线辨识系统及惯量算法
CN112326117B (zh) * 2020-10-23 2022-07-05 北京控制工程研究所 一种基于cmg的航天器惯量在轨辨识方法
CN112577669B (zh) * 2020-12-15 2022-08-12 上海卫星工程研究所 分体式卫星在轨两舱质心辨识方法及系统
CN113359431B (zh) * 2021-06-17 2023-05-12 北京控制工程研究所 一种针对航天器挠性振动的在线辨识与抑制方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AT381171B (de) * 1985-03-18 1986-09-10 Steyr Daimler Puch Ag Verfahren zum messen des massentraegheitsmomentes des der fahrkupplung nachgeordneten, die triebraeder mitumfassenden antriebsstranges eines kraftfahrzeuges
KR0149369B1 (ko) * 1992-03-17 1998-12-01 전성원 무게중심 및 관성 모멘트 측정방법과 그 장치
US7623534B1 (en) * 2005-09-09 2009-11-24 At&T Intellectual Property I, Lp Method and systems for content access and distribution
CN102620886B (zh) * 2012-03-27 2014-10-15 南京航空航天大学 两步在轨辨识组合航天器转动惯量估计方法
CN103970142A (zh) * 2013-02-01 2014-08-06 上海新跃仪表厂 在轨拖曳的组合体航天器姿轨复合控制方法
CN105259786B (zh) * 2015-10-29 2018-09-25 中国科学院力学研究所 待辨识目标的惯性参数辨识方法和装置

Also Published As

Publication number Publication date
CN106289641A (zh) 2017-01-04

Similar Documents

Publication Publication Date Title
CN106289641B (zh) 挠性航天器质心位置和转动惯量参数联合辨识方法
CN110174899B (zh) 一种基于敏捷卫星的高精度成像姿态指向控制方法
CN112558621A (zh) 一种基于解耦控制的飞行机械臂系统
CN104526695A (zh) 一种最小化基座碰撞扰动的空间机械臂轨迹规划方法
CN109426147B (zh) 捕获卫星后组合航天器的自适应增益调整控制方法
CN105182801B (zh) 一种基于扩张状态观测器的Stewart平台主动隔振PD控制方法
CN114291295B (zh) 一种单磁悬浮控制敏感陀螺的卫星双轴姿态测控一体方法
CN109901402B (zh) 一种基于航向平滑技术的自主水下机器人路径跟踪方法
CN105022402A (zh) 一种双刚体航天器快速机动的最短时间确定方法
CN106502261A (zh) 纳星重构故障卫星姿控功能的辨识控制一体化方法
CN110058603B (zh) 一种形变再入飞行器变形前指令力矩确定方法
CN105912007A (zh) 空间机械臂抗干扰姿态稳定的微分几何非线性控制方法
CN110567462B (zh) 一种近似自旋非合作航天器三轴转动惯量比的辨识方法
CN115366109A (zh) 一种旋翼飞行机械臂复合分层抗干扰方法
CN110509276B (zh) 一种机场跑道检测机器人的运动建模及参数辨识方法
CN103514362B (zh) 基于模型误差补偿的两行根数生成方法
CN117742156A (zh) 一种基于rbf神经网络的四旋翼无人机控制方法及系统
CN110955255B (zh) 基于cmg的高精度轨控姿态维持方法、系统及介质
CN118111430A (zh) 基于最小误差熵卡尔曼的交互多模型auv组合导航方法
CN102749846A (zh) 一种两平行构型VSDGCMGs奇异回避操纵律设计方法
CN104850128B (zh) 一种用于具有大惯量积航天器的动量轮布局配置方法
CN111546344A (zh) 一种用于对准的机械臂控制方法
Diao et al. An output feedback attitude tracking controller design for quadrotor unmanned aerial vehicles using quaternion
CN107870063B (zh) 基于动量守恒的航天器转动惯量在轨测量方法
CN114671050B (zh) 基于一体化线性算子和抗饱和技术的航天器跟踪控制方法

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