CN110316402A - 一种编队控制模式下的卫星姿态控制方法 - Google Patents
一种编队控制模式下的卫星姿态控制方法 Download PDFInfo
- Publication number
- CN110316402A CN110316402A CN201910480432.XA CN201910480432A CN110316402A CN 110316402 A CN110316402 A CN 110316402A CN 201910480432 A CN201910480432 A CN 201910480432A CN 110316402 A CN110316402 A CN 110316402A
- Authority
- CN
- China
- Prior art keywords
- formation
- thruster
- control
- satellite
- disturbance torque
- 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.)
- Pending
Links
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 117
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000005484 gravity Effects 0.000 claims abstract description 20
- 230000009471 action Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims 1
- 238000013461 design Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000007774 longterm Effects 0.000 description 3
- 230000035508 accumulation Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000003032 molecular docking Methods 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
- B64G1/245—Attitude control algorithms for spacecraft attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/28—Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect
- B64G1/283—Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect using reaction wheels
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
一种编队控制模式下的卫星姿态控制方法,在由编队推力器干扰力矩引起的星体三轴角动量的积累值小于等于飞轮可吸收的角动量值时,采用反作用飞轮实现卫星姿态控制,通过反作用飞轮PID控制对编队推力器干扰力矩进行补偿,在由编队推力器干扰力矩引起的星体三轴角动量的积累值大于飞轮可吸收的角动量值时,采用推力器实现卫星姿态控制,通过推力器PD控制对编队推力器干扰力矩进行补偿。本发明实现了满足编队飞行的姿态控制精度和卫星快速姿态稳定控制要求,提高了卫星在轨编队飞行工作寿命和在轨可靠性,减少了卫星编队飞行姿态控制对地面测控资源的依赖。
Description
技术领域
本发明涉及卫星编队控制技术领域,尤其涉及一种编队控制模式下的卫星姿态控制方法。
背景技术
目前,卫星的编队任务大致可以分为合作目标的编队飞行、空间交会对接和空间服务几大类。空间交会对接和空间服务等短期任务性质的卫星编队已实现工程应用,而合作目标卫星编队按构形设计要求长期高可靠飞行尚处于工程研制阶段。
卫星编队飞行控制需要根据任务要求进行构形初始化、构形保持、构形切换、构形重构等编队控制,与传统的单星航天器轨道飞行控制不同,编队卫星除基于精准的绝对轨道测量,同时需根据星间相对运动状态精准测量,采用更精确的推力脉冲控制实现在轨自主运行。
编队卫星长期绕飞,因星上燃料携带资源有限,地面测控站资源受限,卫星编队控制需重点进行可靠性设计,从卫星寿命保证及编队任务需求覆盖等,提高编队卫星在轨自主性水平。编队卫星在轨自主性管理主要反映在两个方面:一个是在轨自主姿态控制,一个是在轨自主轨道控制,而轨道控制必定产生姿态干扰。编队卫星姿态控制,需结合卫星具体编队任务,选择合适的卫星姿态控制执行机构,减小卫星姿态控制长期依靠推力器的依赖,保证卫星在轨运行寿命。
发明内容
本发明提供一种编队控制模式下的卫星姿态控制方法,实现了满足编队飞行的姿态控制精度和卫星快速姿态稳定控制要求,提高了卫星在轨编队飞行工作寿命和在轨可靠性,减少了卫星编队飞行姿态控制对地面测控资源的依赖。
为了达到上述,本发明提供一种编队控制模式下的卫星姿态控制方法,在由编队推力器干扰力矩引起的星体三轴角动量的积累值小于等于飞轮可吸收的角动量值时,采用反作用飞轮实现卫星姿态控制,通过反作用飞轮PID控制对编队推力器干扰力矩进行补偿,在由编队推力器干扰力矩引起的星体三轴角动量的积累值大于飞轮可吸收的角动量值时,采用推力器实现卫星姿态控制,通过推力器PD控制对编队推力器干扰力矩进行补偿。
所述的编队推力器干扰力矩∑ΔM包含因质心变化引起的干扰力矩ΔM1、因编队推力器引起的干扰力矩ΔM2、以及因推力不对称误差引起的干扰力矩ΔM3:∑ΔM=ΔM1+ΔM2+ΔM3。
因质心变化引起的干扰力矩ΔM1:
展开后为:
其中,质心变化引起的力臂变化为ΔRi=[ΔRx ΔRy ΔRz]T,F为根据编队推力器安装布局对应到等效控制方向的推力大小。
因编队推力器引起的干扰力矩ΔM2:
ΔM2=R×ΔF;
其中,ΔF为引起作用方向的干扰推力:
Fx的偏斜角为α。
因推力不对称误差引起的干扰力矩ΔM3:
其中,ΔF为不对称推力:
δ为成对工作的推力器推力之差的绝对值与平均推力的百分比。
在编队构形初始化,或构形切换,或构形重构时,采用推力器实现卫星姿态控制,在编队构形保持时,采用反作用飞轮实现卫星姿态控制。
采用推力器实现卫星姿态控制时,对于每个采样周期,姿控推力器产生的冲量矩应等效于控制指令力矩产生的冲量矩,即Mi·ton=Tci·ts,Mi为推力器冲量,ts为控制周期,通过调制脉宽ton获取所需的控制力矩Tci,控制力矩计算公式如下:
其中,为姿态角估计值,为姿态角速度估计值,其中KP、Kd分别为控制参数。
采用反作用飞轮实现卫星姿态控制时,反作用飞轮的控制力矩为:
其中,KP、Ki、Kd分别为控制参数,TC为控制力矩,为误差姿态角,
ωe为误差姿态角速度,e,a代表任意数的限幅处理。
本发明具有以下有益效果:通过星体动力学方程引入编队推力器干扰力矩项,设计不同执行机构的控制算法,实现了满足编队飞行的姿态控制精度和卫星快速姿态稳定控制要求,提高了卫星在轨编队飞行工作寿命和在轨可靠性,减少了卫星编队飞行姿态控制对地面测控资源的依赖。
附图说明
图1是本发明提供的一种编队控制模式下的卫星姿态控制方法的流程图。
图2是编队控制推力器控制框图。
图3是编队控制执行机构采用推力器控制的姿态角控制曲线。
图4是编队控制执行机构采用推力器控制的姿态角速度控制曲线。
图5是编队控制反作用飞轮控制框图。
图6是编队控制执行机构采用反作用飞轮控制的姿态角控制曲线。
图7是编队控制执行机构采用反作用飞轮控制的姿态角速度控制曲线。
具体实施方式
以下根据图1~图7,具体说明本发明的较佳实施例。
目前,空间合作目标的编队卫星长期绕飞控制任务尚处于工程研制阶段,卫星编队控制长期绕飞国内还未有可参考的工程应用。通过对编队卫星任务分析,同时考虑卫星长期编队绕飞控制需求,本发明提出采用反作用飞轮、推力器两种不同执行机构引入姿态控制的设计方法。
如图1所示,本发明提供一种编队控制模式下的卫星姿态控制方法,包含以下步骤:
步骤S1、计算编队控制模式下的编队推力器干扰力矩;
步骤S2、判断由编队推力器干扰力矩引起的星体三轴角动量的积累值是否小于等于飞轮可吸收的角动量值,若否,进行步骤S3,若是,进行步骤S4;
步骤S3、采用推力器实现卫星姿态控制,通过推力器PD控制对编队推力器干扰力矩进行补偿;
步骤S4、采用反作用飞轮实现卫星姿态控制,通过反作用飞轮PID控制对编队推力器干扰力矩进行补偿。
具体来说,所述的步骤S1中,分析计算编队控制模式下的编队推力器干扰力矩主要考虑卫星特性,开展动力学分析,并建立动力学方程。
根据卫星运行轨道,建立卫星环境模型。环境模型中主要考虑的外干扰力矩含气动力矩、重力梯度力矩、太阳光压、常值轨道陀螺力矩、卫星剩磁引起的地磁力矩等。外干扰力矩中除地磁力矩为周期项外,其余为常值力矩。在编队控制过程中,干扰力矩变化表现为编队推力器因质心、推力变化等引起的干扰力矩。
编队推力器干扰力矩∑ΔM包含因质心变化引起的干扰力矩ΔM1、因编队推力器引起的干扰力矩ΔM2、以及因推力不对称误差引起的干扰力矩ΔM3。
∑ΔM=ΔM1+ΔM2+ΔM3 (公式1)。
根据编队控制推力器布局,考虑推力器推力大小和质心变化,以卫星本体控制坐标系为基准,计算干扰力矩变化量,F为根据编队推力器安装布局对应到等效控制方向的推力大小。
设质心变化引起的力臂变化为ΔRi=[ΔRx ΔRy ΔRz]T,则干扰力矩展开后为公式2:
推力器误差考虑因制造及安装误差引起的推力方向偏斜导致的干扰误差,以星体X向编队控制的推力器为例,若Fx的偏斜角为α,引起作用方向的干扰推力为公式3,干扰力矩为公式4:
ΔM2=R×ΔF (公式4)
其余方向编队推力器引起的干扰力矩类此计算,获取各向编队推力器引起的干扰推力。
编队控制下,因推力不对称误差引起的干扰力矩,设推力不对称为δ(δ为成对工作的推力器推力之差的绝对值与平均推力的百分比),不对称推力为公式5,形成的干扰力矩为公式6:
所述的步骤S2中,卫星编队控制包含编队构形初始化、构形保持、构形切换、构形重构等。编队构形初始化、构形切换、构形重构时,编队推力器喷气时间过长,由编队推力器干扰力矩带来的角动量累积量比较大,若在飞轮控制下进行将引起飞轮角动量的输出饱和,卫星姿态将失控,处于这种编队时长控制状态时,主要采用推力器进行姿态控制。当编队构形已经建立和构形重构完成后,编队进入保持控制时,星体三轴角动量积累角动量在飞轮角动量输出的正常范围之内,完全可以被飞轮吸收,采用反作用飞轮进行星体姿态控制。
其中飞轮输出力矩以飞轮最大力矩输出设计值为约束,充分利用飞轮一阶惯性的微分特性,通过对卫星动力学分析,并隔离大部件的振动干扰,简化动力学方程,如公式7所示。
上式中(X,Y,Z代表星体三轴干扰力矩分量)为作用在星体上的干扰力矩;(X,Y,Z代表星体三轴控制力矩分量)为作用在星体上的控制力矩,由反作用飞轮角动量变化产生;I为卫星惯量(X,Y,Z代表星体三轴主惯量,YZ为星体YZ向交惯量)。当外干扰力矩使星体三轴角动量积累并超过飞轮角动量可吸收的能力时,卫星姿态将会失控。
所述的步骤S3中,卫星编队姿态控制采用推力器时,卫星姿态动力学需将编队推力器工作时引起的干扰力矩作为常值项进行控制,并对编队推力器进行角度解耦偏置补偿,对反作用飞轮转速进行处理,避免反作用飞轮产生干扰力矩。
推力器控制采用脉宽调制设计,将控制力矩指令转化为输给推进系统可以执行的脉冲指令。对于每个采样周期,姿控推力器产生的冲量矩应等效于控制指令力矩产生的冲量矩,即Mi·ton=Tci·ts,Mi为推力器冲量,ts为控制周期,通过调制脉宽ton获取所需的控制力矩Tci,Tci计算见公式8。
推力器控制方案采用星敏+陀螺或轨道罗盘确定三轴姿态和姿态角速度信息,以(0°,0°,0°)为姿态控制目标值,采用PD(比例微分)控制,输出相应脉宽的喷气量。为了限制各轴的角速度,对姿态角加入限幅环节。
控制力矩计算如公式7:
其中,为姿态角估计值;为姿态角速度估计值,其中KP、Kd分别为控制参数。
按照图2推力器控制框图,进行星体动力学分析,根据公式1,计算编队推力器干扰力矩,并作为常值项干扰力矩。考虑编队推力器的耦合影响,根据相应方向编队推力器布局产生的耦合力矩进行偏置角度补偿。取卫星初始姿态角为[2 2 2]°,初始姿态角速度为[0.3 0.3 0.3]°/s,以编队保持控制2000s时间长度作用下,采用推力器进行姿态控制,系统的三轴姿态角速度估值在±0.1°/s以内,其三轴姿态角度估值为0.5°以内,达到系统控制的精度和稳定度要求,仿真曲线如图3和图4所示,其中编队推力器工作时间为1000s~3000s。
所述的步骤S4中,卫星编队姿态控制采用反作用飞轮时,需消除因编队推力器工作产生的干扰力矩,并根据飞轮配置和安装构形,分配控制矩阵,发送飞轮控制指令。为保证卫星姿态控制裕度,对飞轮输出进行力矩限幅,并进行卫星姿态限幅,使卫星保持一定的姿态精度和稳定度。
计算编队推力器干扰力矩,根据飞轮配置和安装构形,确定反作用飞轮三轴的最大输出力矩,并保证反作用飞轮三轴合成力矩足够克服编队推力器干扰力矩的裕度,建议取2以上,从而消除因编队推力器干扰力矩对星体的姿态稳定度影响。
由于飞轮输出的角动量有限,一定时间内编队推力器干扰力矩对星体角动量会带来累积效应。因此,编队控制时一般仅在编队保持下选择反作用飞轮作为执行机构,进行卫星姿态控制。
在编队保持控制下引入角速度误差和姿态误差,设计PID(比例积分微分)控制器,如公式9:
其中KP、Ki、Kd分别为控制参数,TC为控制力矩,为误差姿态角,ωe为误差姿态角速度,定义如公式10:
公式中e,a代表任意数的限幅处理。
控制器中引入积分环节用以消除静态误差。为了实现存在各种受限及约束的前提下满足快速性的要求,对姿态偏差角进行限幅。控制回路中根据卫星姿态任务具体要求可以引入导引控制的角度和角速度信息,克服导引控制带来的耦合影响,从而实现系统的高稳定度控制。同时,根据接入系统控制的飞轮配置进行星体控制系下力矩指令限幅,避免因飞轮角动量饱和导致星体姿态出现非期望控制。
按照图5飞轮控制框图,进行星体动力学分析,根据公式1计算编队推力器干扰力矩,并将其作为常值项干扰力矩。在编队保持控制下,以编队控制100s进行仿真分析,在2100s处进行编队控制,控制时长100s,编队控制过程中滚动角引起的最大偏差为0.33°,俯仰角为0.21°,偏航角为0.17°,对应各轴各姿态角速度的变化范围为(0.015,0.008,0.007)°/s,仿真曲线如图6和图7所示。编队控制结束后,星体可以快速实现高精度高稳定度姿态控制,达到设计目标。
本发明具有以下有益效果:通过星体动力学方程引入编队推力器干扰力矩项,设计不同执行机构的控制算法,实现了满足编队飞行的姿态控制精度和卫星快速姿态稳定控制要求,提高了卫星在轨编队飞行工作寿命和在轨可靠性,减少了卫星编队飞行姿态控制对地面测控资源的依赖。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (8)
1.一种编队控制模式下的卫星姿态控制方法,其特征在于,在由编队推力器干扰力矩引起的星体三轴角动量的积累值小于等于飞轮可吸收的角动量值时,采用反作用飞轮实现卫星姿态控制,通过反作用飞轮PID控制对编队推力器干扰力矩进行补偿,在由编队推力器干扰力矩引起的星体三轴角动量的积累值大于飞轮可吸收的角动量值时,采用推力器实现卫星姿态控制,通过推力器PD控制对编队推力器干扰力矩进行补偿。
2.如权利要求1所述的编队控制模式下的卫星姿态控制方法,其特征在于,所述的编队推力器干扰力矩∑△M包含因质心变化引起的干扰力矩△M1、因编队推力器引起的干扰力矩△M2、以及因推力不对称误差引起的干扰力矩△M3:∑△M=△M1+△M2+△M3。
3.如权利要求2所述的编队控制模式下的卫星姿态控制方法,其特征在于,因质心变化引起的干扰力矩△M1:
展开后为:
其中,质心变化引起的力臂变化为△Ri=[△Rx △Ry △Rz]T,F为根据编队推力器安装布局对应到等效控制方向的推力大小。
4.如权利要求2所述的编队控制模式下的卫星姿态控制方法,其特征在于,因编队推力器引起的干扰力矩△M2:
△M2=R×△F;
其中,R代表力臂,△F为引起作用方向的干扰推力:
Fx的偏斜角为α。
5.如权利要求2所述的编队控制模式下的卫星姿态控制方法,其特征在于,因推力不对称误差引起的干扰力矩△M3:
其中,R代表力臂,△F为不对称推力:
δ为成对工作的推力器推力之差的绝对值与平均推力的百分比。
6.如权利要求1所述的编队控制模式下的卫星姿态控制方法,其特征在于,在编队构形初始化,或构形切换,或构形重构时,采用推力器实现卫星姿态控制,在编队构形保持时,采用反作用飞轮实现卫星姿态控制。
7.如权利要求1所述的编队控制模式下的卫星姿态控制方法,其特征在于,采用推力器实现卫星姿态控制时,对于每个采样周期,姿控推力器产生的冲量矩应等效于控制指令力矩产生的冲量矩,即Mi·ton=Tci·ts,Mi为推力器冲量,ts为控制周期,通过调制脉宽ton获取所需的控制力矩Tci,控制力矩计算公式如下:
其中,为姿态角估计值,为姿态角速度估计值,其中KP、Kd分别为控制参数。
8.如权利要求1所述的编队控制模式下的卫星姿态控制方法,其特征在于,采用反作用飞轮实现卫星姿态控制时,反作用飞轮的控制力矩为:
其中,KP、Ki、Kd分别为控制参数,TC为控制力矩,为误差姿态角,ωe为误差姿态角速度,e,a代表任意数的限幅处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910480432.XA CN110316402A (zh) | 2019-06-03 | 2019-06-03 | 一种编队控制模式下的卫星姿态控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910480432.XA CN110316402A (zh) | 2019-06-03 | 2019-06-03 | 一种编队控制模式下的卫星姿态控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110316402A true CN110316402A (zh) | 2019-10-11 |
Family
ID=68119459
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910480432.XA Pending CN110316402A (zh) | 2019-06-03 | 2019-06-03 | 一种编队控制模式下的卫星姿态控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110316402A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111319799A (zh) * | 2020-03-02 | 2020-06-23 | 上海航天控制技术研究所 | 一种火星探测环绕器分离安全自主推力系统及方法 |
CN111422379A (zh) * | 2020-03-20 | 2020-07-17 | 上海航天控制技术研究所 | 一种编队卫星协同轨道控制方法 |
CN111443722A (zh) * | 2020-03-23 | 2020-07-24 | 上海航天控制技术研究所 | 一种编队卫星定时段自主保持方法 |
CN112193438A (zh) * | 2020-09-14 | 2021-01-08 | 中国空间技术研究院 | 一种卫星入轨初期放气扰动的评估系统与评估方法 |
CN112591149A (zh) * | 2020-12-01 | 2021-04-02 | 上海航天控制技术研究所 | 周期性干扰力矩辨识及补偿方法 |
CN112607063A (zh) * | 2020-12-29 | 2021-04-06 | 中国人民解放军32039部队 | 一种卫星位置保持的控制方法和装置 |
CN114019992A (zh) * | 2021-10-09 | 2022-02-08 | 北京控制工程研究所 | 一种面向推力器复用的三轴解耦姿态控制方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6089507A (en) * | 1996-12-05 | 2000-07-18 | Parvez; Shabbir Ahmed | Autonomous orbit control with position and velocity feedback using modern control theory |
CN101056079A (zh) * | 2007-04-26 | 2007-10-17 | 北京航空航天大学 | 一种磁悬浮反作用飞轮 |
CN102636818A (zh) * | 2012-04-27 | 2012-08-15 | 清华大学 | 星星跟踪复合编队地球重力场测量系统及其方法 |
CN103072702A (zh) * | 2013-01-30 | 2013-05-01 | 北京控制工程研究所 | 卫星轨道和姿态控制方法 |
CN103112602A (zh) * | 2013-01-30 | 2013-05-22 | 北京控制工程研究所 | 一种基于推力器和动量轮联合的三轴角速度阻尼控制方法 |
CN105955283A (zh) * | 2016-05-30 | 2016-09-21 | 上海航天控制技术研究所 | 多轴快速姿态机动喷气控制方法 |
CN106094854A (zh) * | 2016-07-15 | 2016-11-09 | 中国人民解放军装备学院 | 电磁编队卫星姿态和轨道相对控制方法 |
CN108639384A (zh) * | 2018-04-03 | 2018-10-12 | 上海航天控制技术研究所 | 一种基于姿控推力器复用的卫星推进剂管理方法 |
CN108639385A (zh) * | 2018-05-15 | 2018-10-12 | 浙江大学 | 一种无基准最简姿控系统快速稳定控制的实现方法 |
CN109507875A (zh) * | 2019-01-08 | 2019-03-22 | 哈尔滨工业大学 | 一种欧拉旋转卫星姿态机动递阶饱和pid控制方法 |
-
2019
- 2019-06-03 CN CN201910480432.XA patent/CN110316402A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6089507A (en) * | 1996-12-05 | 2000-07-18 | Parvez; Shabbir Ahmed | Autonomous orbit control with position and velocity feedback using modern control theory |
CN101056079A (zh) * | 2007-04-26 | 2007-10-17 | 北京航空航天大学 | 一种磁悬浮反作用飞轮 |
CN102636818A (zh) * | 2012-04-27 | 2012-08-15 | 清华大学 | 星星跟踪复合编队地球重力场测量系统及其方法 |
CN103072702A (zh) * | 2013-01-30 | 2013-05-01 | 北京控制工程研究所 | 卫星轨道和姿态控制方法 |
CN103112602A (zh) * | 2013-01-30 | 2013-05-22 | 北京控制工程研究所 | 一种基于推力器和动量轮联合的三轴角速度阻尼控制方法 |
CN105955283A (zh) * | 2016-05-30 | 2016-09-21 | 上海航天控制技术研究所 | 多轴快速姿态机动喷气控制方法 |
CN106094854A (zh) * | 2016-07-15 | 2016-11-09 | 中国人民解放军装备学院 | 电磁编队卫星姿态和轨道相对控制方法 |
CN108639384A (zh) * | 2018-04-03 | 2018-10-12 | 上海航天控制技术研究所 | 一种基于姿控推力器复用的卫星推进剂管理方法 |
CN108639385A (zh) * | 2018-05-15 | 2018-10-12 | 浙江大学 | 一种无基准最简姿控系统快速稳定控制的实现方法 |
CN109507875A (zh) * | 2019-01-08 | 2019-03-22 | 哈尔滨工业大学 | 一种欧拉旋转卫星姿态机动递阶饱和pid控制方法 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111319799A (zh) * | 2020-03-02 | 2020-06-23 | 上海航天控制技术研究所 | 一种火星探测环绕器分离安全自主推力系统及方法 |
CN111422379A (zh) * | 2020-03-20 | 2020-07-17 | 上海航天控制技术研究所 | 一种编队卫星协同轨道控制方法 |
CN111422379B (zh) * | 2020-03-20 | 2021-08-10 | 上海航天控制技术研究所 | 一种编队卫星协同轨道控制方法 |
CN111443722A (zh) * | 2020-03-23 | 2020-07-24 | 上海航天控制技术研究所 | 一种编队卫星定时段自主保持方法 |
CN111443722B (zh) * | 2020-03-23 | 2022-09-23 | 上海航天控制技术研究所 | 一种编队卫星定时段自主保持方法 |
CN112193438A (zh) * | 2020-09-14 | 2021-01-08 | 中国空间技术研究院 | 一种卫星入轨初期放气扰动的评估系统与评估方法 |
CN112591149A (zh) * | 2020-12-01 | 2021-04-02 | 上海航天控制技术研究所 | 周期性干扰力矩辨识及补偿方法 |
CN112607063A (zh) * | 2020-12-29 | 2021-04-06 | 中国人民解放军32039部队 | 一种卫星位置保持的控制方法和装置 |
CN112607063B (zh) * | 2020-12-29 | 2021-10-08 | 中国人民解放军32039部队 | 一种卫星位置保持的控制方法和装置 |
CN114019992A (zh) * | 2021-10-09 | 2022-02-08 | 北京控制工程研究所 | 一种面向推力器复用的三轴解耦姿态控制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110316402A (zh) | 一种编队控制模式下的卫星姿态控制方法 | |
CN106873611B (zh) | 一种多通道线性自抗扰控制器的设计方法 | |
CN104527994B (zh) | 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法 | |
US20080315039A1 (en) | System and methods for space vehicle torque balancing | |
CN103991559A (zh) | 一种洛伦兹航天器悬停控制方法 | |
MacKunis et al. | Adaptive satellite attitude control in the presence of inertia and CMG gimbal friction uncertainties | |
CN104309822A (zh) | 一种基于参数优化的航天器单脉冲水滴形绕飞轨迹悬停控制方法 | |
Chung et al. | Propellant-free control of tethered formation flight, part 1: Linear control and experimentation | |
Jang et al. | Design of launch vehicle flight control systems using ascent vehicle stability analysis tool | |
Kusuda et al. | Feedback control with nominal inputs for agile satellites using control moment gyros | |
Subramanian | Nonlinear control strategies for quadrotors and CubeSats | |
Zhang et al. | Piece-wise affine MPC-based attitude control for a CubeSat during orbital manoeuvres | |
Varma | Control of satellites using environmental forces: aerodynamic drag/solar radiation pressure | |
CN105068425A (zh) | 一种适用于敏捷卫星姿态确定的状态反馈鲁棒非脆弱控制方法 | |
CN105739511B (zh) | 一种缺失迹向控制的欠驱动航天器悬停渐近控制方法 | |
Marchetti et al. | Electric propulsion and controller design for drag-free spacecraft operation | |
Cacciatore et al. | The Design of the GNC of the Re-entry Module of Space Rider | |
Xu et al. | PD-type magnetic-based sun-pointing attitude control scheme | |
CN112393835B (zh) | 一种基于扩展卡尔曼滤波的小卫星在轨推力标定方法 | |
Bayat et al. | A heuristic design method for attitude stabilization of magnetic actuated satellites | |
Daitx et al. | Development of a combined attitude and position controller for a satellite simulator | |
Reijneveld et al. | Attitude control system of the Delfi-n3Xt satellite | |
Lee et al. | Developing an automatic control system of unmanned aircrafts with a single-antenna GPS receiver | |
Xie et al. | Autonomous guidance, navigation, and control of spacecraft | |
Hall et al. | Ares i flight control system overview |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191011 |