CN104590587A - 基于时间配比的卫星三轴控制力矩解耦的实现方法 - Google Patents

基于时间配比的卫星三轴控制力矩解耦的实现方法 Download PDF

Info

Publication number
CN104590587A
CN104590587A CN201410706026.8A CN201410706026A CN104590587A CN 104590587 A CN104590587 A CN 104590587A CN 201410706026 A CN201410706026 A CN 201410706026A CN 104590587 A CN104590587 A CN 104590587A
Authority
CN
China
Prior art keywords
satellite
feasible
control
feasible program
control moment
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
CN201410706026.8A
Other languages
English (en)
Other versions
CN104590587B (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.)
Harbin University of Technology Satellite Technology Co.,Ltd.
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201410706026.8A priority Critical patent/CN104590587B/zh
Publication of CN104590587A publication Critical patent/CN104590587A/zh
Application granted granted Critical
Publication of CN104590587B publication Critical patent/CN104590587B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

基于时间配比的卫星三轴控制力矩解耦的实现方法,属于航天器姿态控制技术领域。为了解决卫星推进器布局存在耦合或推进器故障而无法提供三轴解耦控制力矩的问题。它包括:根据实际推进发动机安装位置,求得各个推进器控制力矩;根据时间配比方法,引入推进器控制力矩的工作时间系数通过获得的各个推进器控制力矩与不同时间系数的组合得到实现控制力矩三轴解耦的所有可行方案;根据确定的可行方案的燃料用于姿态控制的效率,确定三轴正负向控制力矩所采用可行方案的优先级,效率越高,对应可行方案优先级越高;根据实际控制信号和确定的优先级最高的可行方案,确定推进器工作时间向量。用于采用推进器对卫星进行姿态控制。

Description

基于时间配比的卫星三轴控制力矩解耦的实现方法
技术领域
本发明属于航天器姿态控制技术领域。
背景技术
在卫星姿态控制领域,常采用飞轮、磁力矩器、推进器作为卫星姿态控制的执行机构。其中,推进器具有控制力矩大、响应速度快等优点,可作为大型卫星姿态稳定以及小卫星快速机动的执行机构。相比其他两种执行机构,采用推进器进行姿态控制时,卫星姿态控制精度往往较低,而且若推进器布局存在耦合,采用推进器对卫星单轴姿态机动时极易引起其他轴姿态变化,更加降低了卫星姿态控制的精度。因此,实际应用中,大型卫星上常采用推进器进行姿态控制,小卫星往往因为对控制精度要求较高而极少仅使用推进器进行姿态控制。
为实现卫星三轴姿态解耦控制,一般需要配置至少12台推进器。但是,受卫星成本、结构、重量等因素的限制,卫星上往往无法部署足够多的推进器。此外,在卫星实际运行时,推进器故障也是不容忽视的,如何在故障情况下提供解耦的三轴控制力矩也成为姿态控制所必须解决的问题。因此,有必要设计一种适当的方法为推进器布局存在耦合的卫星提供三轴解耦的控制力矩,从而保证卫星的姿态控制精度。合理的推进器布局及力矩解耦方法使得将推进器广泛应用到小卫星姿态快速机动成为可能。
发明内容
本发明的目的是为了解决卫星推进器布局存在耦合或推进器故障而无法提供三轴解耦控制力矩的问题,本发明提供一种基于时间配比的卫星三轴控制力矩解耦的实现方法。
本发明的基于时间配比的卫星三轴控制力矩解耦的实现方法,
所述方法包括如下步骤:
步骤一、根据实际推进发动机安装位置,求得各个推进器的控制力矩Tj
Tj=Rj×Fj,Rj为第j个推进器安装位置矢量,Fj为第j个推进器输出推力;
步骤二、根据时间配比方法,引入推进器控制力矩的工作时间系数,并通过步骤一获得的各个推进器的控制力矩与不同时间系数的组合得到能实现控制力矩三轴解耦的所有可行方案:
V={A1,A2,…,Al,…};
Al=(αl1 αl2 … αlj …)T
T L = T A l = Σ j α lj T j ;
T=(T1 T2 … Tj…);
其中,V为实现控制力矩三轴解耦所有可行方案的集合;Al为第l个可行方案,用各个推进器工作时间向量表示;TL为第l个可行方案的合力矩,TL方向应与卫星本体坐标系某一主轴平行,且满足||TL||=1;αlj为可行方案l的第j个推进器的工作时间系数,均为非负值;
步骤三、根据步骤二确定的各个可行方案的燃料用于姿态控制的效率ηl,确定三轴正负向控制力矩所采用可行方案的优先级,效率ηl越高,对应的可行方案优先级越高;
步骤四、根据实际控制信号Uc和步骤三中确定的优先级最高的可行方案,确定推进器工作时间向量。
步骤三中,各个可行方案的燃料用于姿态控制的效率ηl为:
η l = Σ j α lj | | F j | | | | T l | | = Σ j α lj | | F j | | .
步骤四中,根据实际控制信号Uc和步骤三中确定的优先级最高的可行方案,确定推进器工作时间向量:
T o = Σ i = x , y , z u ci A i = u cx A x + u cy A y + u cz A z
其中,uci为控制信号Uc在卫星本体坐标系下i主轴的分量,Ai为在卫星本体坐标系下i主轴对应正向或负向的优先级最高的可行方案。
本发明的有益效果在于,本发明实现了卫星推进器布局存在耦合情况下为三轴姿态提供解耦控制力矩与控制信号到各发动机工作时间的转化,具有如下优点:
(1)推进器耦合布局下,提出时间系数配比的思想,通过推进器组合实现三轴控制力矩的解耦,使卫星各轴姿态控制互不影响;
(2)可得到任一耦合布局时推进器控制力矩解耦的所有可行方案,方便寻优并降低燃料消耗;
(3)推进器故障时,可寻找到燃耗最优的替补推进器与组合方案而不影响控制精度。
附图说明
图1为具体实施方式一所述的基于时间配比的卫星三轴控制力矩解耦的实现方法的原理示意图。
图2为四姿控推进器倾斜布局结构图。
具体实施方式
结合图1说明本实施方式,本实施方式所述的基于时间配比的卫星三轴控制力矩解耦的实现方法,所述方法包括如下步骤:
步骤一、根据实际推进发动机安装位置,求得各个推进器的控制力矩Tj
Tj=Rj×Fj,Rj为第j个推进器安装位置矢量,Fj为第j个推进器输出推力;这些矢量均在本体坐标系下描述。
步骤二、根据时间配比方法,引入推进器控制力矩的工作时间系数,并通过步骤一获得的各个推进器的控制力矩与不同时间系数的组合得到能实现控制力矩三轴解耦的所有可行方案:
V={A1,A2,…,Al,…};
Al=(αl1 αl2 … αlj …)T
T L = T A l = Σ j α lj T j ;
T=(T1 T2 … Tj…);
其中,V为实现控制力矩三轴解耦所有可行方案的集合;Al为第l个可行方案,用各个推进器工作时间向量表示;TL为第l个可行方案的合力矩,TL方向应与卫星本体坐标系某一主轴平行,且满足||TL||=1;αlj为可行方案l的第j个推进器的工作时间系数,均为非负值。
该步骤所述时间配比思想为:由于推进器控制力矩耦合时要实现卫星姿态三轴控制必然满足推进器控制力矩可作为空间向量的基底,因而通过引入时间系数必然可得到卫星各个主轴方向的控制力矩,再对该力矩归一化即可得到各推进器的时间系数。
步骤三、根据步骤二确定的各个可行方案的燃料用于姿态控制的效率ηl,确定三轴正负向控制力矩所采用可行方案的优先级,效率ηl越高,对应的可行方案优先级越高;
η l = Σ j α lj | | F j | | | | T l | | = Σ j α lj | | F j | | .
步骤四、根据实际控制信号Uc和步骤三中确定的优先级最高的可行方案,确定推进器工作时间向量:
T o = Σ i = x , y , z u ci A i = u cx A x + u cy A y + u cz A z
其中,uci为控制信号Uc在卫星本体坐标系下i主轴的分量,Ai为在卫星本体坐标系下i主轴对应正向(uci>0)或负向(uci≤0)的优先级最高的可行方案。
x、y和z分别表示卫星本体坐标系下的三轴。
具体实施例:
以卫星常用的四姿控推进器倾斜布局形式为例,如图2,说明所设计方法的合理性,具体过程为:
各推进器安装位置矢量:
R1=R4=(-0.1 0.1 0)Tm    R2=R3=(-0.1 -0.1 0)Tm
各推进器推力矢量:
F1=F2=(-5 0 -5)TN    F3=F4=(-5 0 5)TN
可得推进器控制力矩为:
T1=(-0.5 -0.5 0.5)TNm T2=(0.5 -0.5 -0.5)TNm
T3=(-0.5 0.5 -0.5)TNm T4=(0.5 0.5 0.5)TNm
通过推进器组合与时间系数加权可得各轴正负向控制推进器组合:
X轴正向组合“2+4”,X轴负向组合“1+3”;
Y轴正向组合“3+4”,Y轴负向组合“1+2”;
Z轴正向组合“1+4”,Z轴负向组合“2+3”;
所有组合推进器权值之比均为1:1。
则所用方法能寻找到所有可行组合,且仅提供单轴控制力矩。

Claims (3)

1.基于时间配比的卫星三轴控制力矩解耦的实现方法,其特征在于,所述方法包括如下步骤:
步骤一、根据实际推进发动机安装位置,求得各个推进器的控制力矩Tj
Tj=Rj×Fj,Rj为第j个推进器安装位置矢量,Fj为第j个推进器输出推力;
步骤二、根据时间配比方法,引入推进器控制力矩的工作时间系数,并通过步骤一获得的各个推进器的控制力矩与不同时间系数的组合得到能实现控制力矩三轴解耦的所有可行方案:
V={A1,A2,…,Al,…};Al=(αl1 αl2 … αlj …)T
T L = TA l = Σ j α lj T j ; T = T 1 T 2 · · · T j · · · ;
其中,V为实现控制力矩三轴解耦所有可行方案的集合;Al为第l个可行方案,用各个推进器工作时间向量表示;TL为第l个可行方案的合力矩,TL方向应与卫星本体坐标系某一主轴平行,且满足||TL||=1;αlj为可行方案l的第j个推进器的工作时间系数,均为非负值;
步骤三、根据步骤二确定的各个可行方案的燃料用于姿态控制的效率ηl,确定三轴正负向控制力矩所采用可行方案的优先级,效率ηl越高,对应的可行方案优先级越高;
步骤四、根据实际控制信号Uc和步骤三中确定的优先级最高的可行方案,确定推进器工作时间向量。
2.根据权利要求1所述的基于时间配比的卫星三轴控制力矩解耦的实现方法,其特征在于,
步骤三中,各个可行方案的燃料用于姿态控制的效率ηl为:
η l = Σ j α lj | | F j | | | | T l | | = Σ j α lj | | F j | | .
3.根据权利要求1或2所述的基于时间配比的卫星三轴控制力矩解耦的实现方法,其特征在于,步骤四中,根据实际控制信号Uc和步骤三中确定的优先级最高的可行方案,确定推进器工作时间向量: T o = Σ i = x , y , z u ci A i = u cx A x + u cy A y + u cz A z
其中,uci为控制信号Uc在卫星本体坐标系下i主轴的分量,Ai为在卫星本体坐标系下i主轴对应正向或负向的优先级最高的可行方案。
CN201410706026.8A 2014-11-27 2014-11-27 基于时间配比的卫星三轴控制力矩解耦的实现方法 Active CN104590587B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410706026.8A CN104590587B (zh) 2014-11-27 2014-11-27 基于时间配比的卫星三轴控制力矩解耦的实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410706026.8A CN104590587B (zh) 2014-11-27 2014-11-27 基于时间配比的卫星三轴控制力矩解耦的实现方法

Publications (2)

Publication Number Publication Date
CN104590587A true CN104590587A (zh) 2015-05-06
CN104590587B CN104590587B (zh) 2016-08-31

Family

ID=53116730

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410706026.8A Active CN104590587B (zh) 2014-11-27 2014-11-27 基于时间配比的卫星三轴控制力矩解耦的实现方法

Country Status (1)

Country Link
CN (1) CN104590587B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104909159A (zh) * 2015-05-29 2015-09-16 楚天科技股份有限公司 一种冻干机联线
CN108820257A (zh) * 2018-06-25 2018-11-16 上海卫星工程研究所 一种卫星用步进电机机构的力矩校核方法
CN109901600A (zh) * 2019-03-08 2019-06-18 宁波天擎航天科技有限公司 一种航天器飞行控制方法、系统及装置
CN111232248A (zh) * 2020-02-14 2020-06-05 哈尔滨工业大学 一种基于轨控推力器脉宽调制的姿轨一体化控制方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2115008C1 (ru) * 1995-02-23 1998-07-10 Рылов Юрий Павлович Устройство управления движением космического аппарата
US5816538A (en) * 1994-10-13 1998-10-06 Hughes Electronics Corporation Dynamic decoupler for improved attitude control
EP1616790A1 (en) * 2004-07-13 2006-01-18 Korea Advanced Institute of Science and Technology Method for attitude control of satellites in elliptic orbits using solar radiation pressure
CN102778891A (zh) * 2012-08-03 2012-11-14 北京理工大学 一种星上控制力矩陀螺群隔振平台的参数选择方法
CN103010486A (zh) * 2012-12-11 2013-04-03 北京控制工程研究所 一种多分支、多路推力器交叉组合喷气分配方法
CN104020778A (zh) * 2014-06-18 2014-09-03 哈尔滨工业大学 基于跟踪时间-能耗最优轨线的挠性卫星姿态机动控制方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5816538A (en) * 1994-10-13 1998-10-06 Hughes Electronics Corporation Dynamic decoupler for improved attitude control
RU2115008C1 (ru) * 1995-02-23 1998-07-10 Рылов Юрий Павлович Устройство управления движением космического аппарата
EP1616790A1 (en) * 2004-07-13 2006-01-18 Korea Advanced Institute of Science and Technology Method for attitude control of satellites in elliptic orbits using solar radiation pressure
CN102778891A (zh) * 2012-08-03 2012-11-14 北京理工大学 一种星上控制力矩陀螺群隔振平台的参数选择方法
CN103010486A (zh) * 2012-12-11 2013-04-03 北京控制工程研究所 一种多分支、多路推力器交叉组合喷气分配方法
CN104020778A (zh) * 2014-06-18 2014-09-03 哈尔滨工业大学 基于跟踪时间-能耗最优轨线的挠性卫星姿态机动控制方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104909159A (zh) * 2015-05-29 2015-09-16 楚天科技股份有限公司 一种冻干机联线
CN108820257A (zh) * 2018-06-25 2018-11-16 上海卫星工程研究所 一种卫星用步进电机机构的力矩校核方法
CN108820257B (zh) * 2018-06-25 2020-05-29 上海卫星工程研究所 一种卫星用步进电机机构的力矩校核方法
CN109901600A (zh) * 2019-03-08 2019-06-18 宁波天擎航天科技有限公司 一种航天器飞行控制方法、系统及装置
CN111232248A (zh) * 2020-02-14 2020-06-05 哈尔滨工业大学 一种基于轨控推力器脉宽调制的姿轨一体化控制方法
CN111232248B (zh) * 2020-02-14 2021-07-27 哈尔滨工业大学 一种基于轨控推力器脉宽调制的姿轨一体化控制方法

Also Published As

Publication number Publication date
CN104590587B (zh) 2016-08-31

Similar Documents

Publication Publication Date Title
Jiang et al. Fixed-time rendezvous control of spacecraft with a tumbling target under loss of actuator effectiveness
CN104590587A (zh) 基于时间配比的卫星三轴控制力矩解耦的实现方法
CN104536452B (zh) 基于时间‑燃料最优控制的航天器相对轨道转移轨迹优化方法
CN102566578B (zh) 基于奇异值分解的单框架控制力矩陀螺群协调控制方法
CN103076807B (zh) 一种欠驱动柔性航天器姿态稳定的控制方法
CN102073280A (zh) 一种复杂挠性航天器模糊奇异摄动建模与姿态控制方法
CN104898686B (zh) 一种基于反作用轮摩擦特性的抗干扰姿态控制方法
CN103699069A (zh) 微小卫星先进电子集成系统
CN104656447A (zh) 一种航天器抗干扰姿态跟踪的微分几何非线性控制方法
CN102923317B (zh) 一种适用于卫星姿态角速度阻尼的欠驱动控制方法
Yang Direct optimization of low-thrust many-revolution earth-orbit transfers
CN104794311A (zh) 一种时间自由条件下的全局最优双脉冲轨道转移方法
CN104914873B (zh) 一种姿轨控发动机的耦合方法
CN103336528B (zh) 一种欠驱动航天器三轴姿态稳定控制方法
Shen et al. Point return orbit design and characteristics analysis for manned lunar mission
Wang et al. Six-DOF trajectory optimization for reusable launch vehicles via Gauss pseudospectral method
CN107187617A (zh) 一种用于航天器接管控制的多单元力矩分配方法
EP3457230B1 (en) Single-axis pointing pure magnetic control algorithm for spacecraft based on geometrical analysis
CN107491082A (zh) 航天器姿态控制混合执行机构优化控制方法
CN103010486A (zh) 一种多分支、多路推力器交叉组合喷气分配方法
CN102426020B (zh) 一种航姿系统的地球自转误差补偿方法
CN103345275B (zh) 基于角动量裕度的单轴批量反作用飞轮力矩优化分配方法
CN109799704A (zh) 一种挠性航天器姿态与结构振动的耦合控制方法
Wu et al. Study on multi-parameter coupling mechanism of a dual-channel controlled ejection seat
CN206421232U (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
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20211228

Address after: Room 1107, 11 / F, National University Science Park, Harbin Institute of technology, No. 434, youyou street, Nangang District, Harbin City, Heilongjiang Province

Patentee after: Harbin Institute of Technology Asset Management Co.,Ltd.

Address before: 150001 No. 92 West straight street, Nangang District, Heilongjiang, Harbin

Patentee before: HARBIN INSTITUTE OF TECHNOLOGY

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220119

Address after: 150028 6 / F, building 1, innovation and entrepreneurship Plaza, science and technology innovation city, high tech Industrial Development Zone, Harbin City, Heilongjiang Province

Patentee after: Harbin University of Technology Satellite Technology Co.,Ltd.

Address before: Room 1107, 11 / F, National University Science Park, Harbin Institute of technology, No. 434, youyou street, Nangang District, Harbin City, Heilongjiang Province

Patentee before: Harbin Institute of Technology Asset Management Co.,Ltd.