CN105644810B - A kind of opened loop control reaction wheel rotation speed change and discharging method - Google Patents

A kind of opened loop control reaction wheel rotation speed change and discharging method Download PDF

Info

Publication number
CN105644810B
CN105644810B CN201610081918.2A CN201610081918A CN105644810B CN 105644810 B CN105644810 B CN 105644810B CN 201610081918 A CN201610081918 A CN 201610081918A CN 105644810 B CN105644810 B CN 105644810B
Authority
CN
China
Prior art keywords
flywheel
torque
unloading
control
reaction wheel
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.)
Expired - Fee Related
Application number
CN201610081918.2A
Other languages
Chinese (zh)
Other versions
CN105644810A (en
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 Engineering Center for Microsatellites
Original Assignee
Shanghai Engineering Center for Microsatellites
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 Engineering Center for Microsatellites filed Critical Shanghai Engineering Center for Microsatellites
Priority to CN201610081918.2A priority Critical patent/CN105644810B/en
Publication of CN105644810A publication Critical patent/CN105644810A/en
Application granted granted Critical
Publication of CN105644810B publication Critical patent/CN105644810B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/244Spacecraft control systems

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)
  • Magnetic Bearings And Hydrostatic Bearings (AREA)
  • Electric Propulsion And Braking For Vehicles (AREA)

Abstract

本发明提供了一种开环控制反作用轮(简称飞轮)转速变化及卸载方法,具体步骤如下:步骤一,采用四个飞轮进行控制,并将第四飞轮的反作用轮控制力矩指令大小调整至与理论标称转速摩擦力矩大小一致;步骤二,将控制器输出控制指令限幅在最大输出力矩,并按照飞轮分配矩阵分配每个飞轮的控制力矩,分配后将输出力矩最大的飞轮限幅在最大力矩,并将其它两个飞轮的控制力矩成比例缩放;步骤三,计算四个飞轮合成角动量,任意轴角动量达到卸载阈值范围内,并符合卸载时磁场强度和角动量夹角关系时,启动磁力矩器卸载;步骤四,选择卸载能力最大的方向上的磁力矩器工作,采用开关控制方式,对飞轮进行卸载。

The invention provides a method for open-loop control of the speed change and unloading of the reaction wheel (flywheel for short), the specific steps are as follows: step 1, using four flywheels for control, and adjusting the control torque command of the reaction wheel of the fourth flywheel to the same size as The friction torque of the theoretical nominal speed is the same; step 2, limit the output control command of the controller to the maximum output torque, and distribute the control torque of each flywheel according to the flywheel distribution matrix, and limit the flywheel with the maximum output torque to the maximum output torque after distribution. torque, and scale the control torques of the other two flywheels proportionally; Step 3, calculate the synthetic angular momentum of the four flywheels, when the angular momentum of any axis reaches the unloading threshold range, and conforms to the relationship between the magnetic field strength and the angular momentum when unloading, Start the unloading of the magnetic torque device; Step 4, select the magnetic torque device in the direction with the largest unloading capacity to work, and use the switch control method to unload the flywheel.

Description

一种开环控制反作用轮转速变化及卸载方法An open-loop control reaction wheel speed change and unloading method

技术领域technical field

本发明涉及航天控制技术领域,尤其涉及一种开环控制反作用轮转速变化及卸载方法。The invention relates to the technical field of aerospace control, in particular to an open-loop control method for changing the rotational speed of a reaction wheel and unloading it.

背景技术Background technique

反作用轮是卫星常用执行机构,由于卫星姿态精度要求越来越高,卫星寿命要求也越来越高,因此反作用轮在进行姿态控制时,将其控制在最佳工作转速附近工作,避免转速过零产生干扰影响姿态精度,同时提高寿命。The reaction wheel is a common actuator for satellites. Due to the higher and higher requirements for satellite attitude accuracy and satellite life, the reaction wheel should be controlled to work near the optimum operating speed when performing attitude control to avoid excessive speed. Zero interference affects attitude accuracy while improving lifespan.

采用4个飞轮进行姿态控制,可将飞轮偏置在标称转速附近形成整星零动量,对于力矩控制模式反作用飞轮在摩擦力矩的作用下转速会下降,当四个反作用轮转速同时下降时,也可以保证姿态稳定和整星零动量,即无法通过姿态闭环反馈的方式保持飞轮转速,因此要保持飞轮在标称转速附近需要对每个飞轮进行摩擦补偿,但由于每个飞轮的摩擦力矩不同,补偿力矩导致飞轮标称转速不同,影响飞轮的姿态控制。另外,空间环境干扰力矩会导致飞轮偏离标称转速,通常通过磁力矩器卸载抵消干扰,确保飞轮工作在标称转速附近。Using 4 flywheels for attitude control, the flywheels can be biased near the nominal speed to form zero momentum. For the torque control mode, the speed of the reaction flywheel will decrease under the action of friction torque. When the speed of the four reaction wheels decreases at the same time, It can also ensure stable attitude and zero momentum throughout the star, that is, it is impossible to maintain the flywheel speed through attitude closed-loop feedback. Therefore, to keep the flywheel near the nominal speed requires friction compensation for each flywheel, but because the friction torque of each flywheel is different , the compensation torque leads to the difference in the nominal speed of the flywheel, which affects the attitude control of the flywheel. In addition, the space environment disturbance torque will cause the flywheel to deviate from the nominal speed, and the disturbance is usually offset by unloading the magnetic torque device to ensure that the flywheel works near the nominal speed.

发明内容Contents of the invention

为了解决上述问题,本发明提供了一种开环控制反作用轮转速变化及磁力矩器卸载方法,包括如下步骤:步骤一,采用包括第一飞轮、第二飞轮、第三飞轮和第四飞轮的四个飞轮进行控制,并将第四飞轮的反作用轮控制力矩指令大小调整至与理论标称转速摩擦力矩大小一致;步骤二,将控制器输出控制指令限幅在最大输出力矩,并按照第一飞轮、第二飞轮、第三飞轮顺序构成的飞轮分配矩阵分配每个飞轮的控制力矩,分配后将输出力矩最大的飞轮限幅在最大力矩,并将其它两个飞轮的控制力矩成比例缩放,得到三个飞轮控制力矩;步骤三,计算四个反作用轮合成角动量,任意轴角动量达到卸载阈值范围内,并符合卸载时磁场强度和角动量夹角关系时,启动磁力矩器卸载;步骤四,选择卸载能力最大的方向上的磁力矩器工作,采用开关控制方式,对飞轮进行卸载。In order to solve the above problems, the present invention provides a method for open-loop control of reaction wheel speed variation and magnetic torque device unloading, comprising the following steps: Step 1, adopting a method comprising the first flywheel, the second flywheel, the third flywheel and the fourth flywheel The four flywheels are controlled, and the reaction wheel control torque command of the fourth flywheel is adjusted to be consistent with the friction torque of the theoretical nominal speed; step 2, the controller output control command is limited to the maximum output torque, and according to the first The flywheel allocation matrix composed of flywheel, second flywheel, and third flywheel in sequence distributes the control torque of each flywheel. After distribution, the flywheel with the largest output torque is limited to the maximum torque, and the control torque of the other two flywheels is scaled proportionally. Get three flywheel control torques; step 3, calculate the synthetic angular momentum of the four reaction wheels, and when the angular momentum of any shaft reaches the unloading threshold range, and conforms to the angular relationship between the magnetic field strength and the angular momentum during unloading, start the magnetic torquer unloading; step Fourth, select the magnetic torque device in the direction with the largest unloading capacity to work, and use the switch control method to unload the flywheel.

本发明的优点在于,利用开环方式仅对一个飞轮进行摩擦补偿,该飞轮不参与姿态控制,其它三个飞轮进行姿态控制,由于角动量守恒,全部飞轮均保持在标称转速附近。该方法控制简单,通过调整补偿摩擦力矩的大小还可以控制飞轮工作标称转速,使用灵活方便。The advantage of the present invention is that only one flywheel is used for friction compensation in an open-loop manner, the flywheel does not participate in attitude control, and the other three flywheels perform attitude control. Due to the conservation of angular momentum, all flywheels are kept near the nominal speed. The method is simple to control, and the nominal rotational speed of the flywheel can also be controlled by adjusting the size of the compensation friction torque, which is flexible and convenient to use.

空间环境干扰力矩产生角动量累积导致飞轮转速偏离标称值,通过磁力矩器卸载可以抵消干扰力矩影响。本发明采用内源场作为磁卸载的参考磁场,无需安装磁强计,简化了卫星的硬件配置;同时在充分分析内、外源场的基础上设计卸载控制器的阈值保证卸载控制的正确性;另外采用开关控制方式控制磁力矩器工作,大大简化了驱动控制电路的设计。The accumulation of angular momentum generated by the disturbance torque of the space environment causes the flywheel speed to deviate from the nominal value, and the unloading of the magnetic torque device can offset the influence of the disturbance torque. The present invention adopts the internal source field as the reference magnetic field for magnetic unloading, and does not need to install a magnetometer, which simplifies the hardware configuration of the satellite; at the same time, the threshold value of the unloading controller is designed on the basis of fully analyzing the internal and external source fields to ensure the correctness of unloading control ; In addition, the switch control method is used to control the work of the magnetic torque device, which greatly simplifies the design of the drive control circuit.

附图说明Description of drawings

附图1是本发明具体实施方式中飞轮的卸载判断条件的判断逻辑图;Accompanying drawing 1 is the judging logic diagram of the unloading judging condition of flywheel in the specific embodiment of the present invention;

附图2是本发明具体实施方式的方案在高轨卫星的应用实施例中卸载策略计算磁力矩器控制指令流程图;Accompanying drawing 2 is the scheme of the specific embodiment of the present invention in the application example of the high-orbit satellite unloading strategy calculation magnetic torque device control instruction flow chart;

附图3是本发明具体实施方式的方案在高轨卫星的应用实施例中卫星在轨反作用轮转速曲线图;Accompanying drawing 3 is the scheme of the specific embodiment of the present invention in the application example of high orbit satellite;

附图4是本发明具体实施方式的方案在高轨卫星的应用实施例中卫星姿态曲线图。Accompanying drawing 4 is the satellite attitude graph in the application example of the high-orbit satellite of the scheme of the specific embodiment of the present invention.

具体实施方式detailed description

下面结合附图对本发明提供的方法的具体实施方式做详细说明。The specific implementation of the method provided by the present invention will be described in detail below in conjunction with the accompanying drawings.

步骤一,采用包括第一飞轮、第二飞轮、第三飞轮和第四飞轮的四个飞轮进行控制,并将第四飞轮的控制力矩指令大小调整至与理论标称转速摩擦力矩大小一致。Step 1: Four flywheels including the first flywheel, the second flywheel, the third flywheel and the fourth flywheel are used for control, and the control torque command of the fourth flywheel is adjusted to be consistent with the friction torque of the theoretical nominal speed.

本步骤中,根据第四飞轮标称转速得到摩擦力矩Tf4,对其进行摩擦力矩补偿,得到反作用轮控制指令Twc4In this step, the friction torque T f4 is obtained according to the nominal speed of the fourth flywheel, and the friction torque is compensated to obtain the reaction wheel control command T wc4 :

Twc4=-Tf4 T wc4 = -T f4

以及反作用轮力矩控制电压:and the reaction wheel torque control voltage:

Uwc4=Twc4/ki4 U wc4 =T wc4 /k i4

其中,ki4为第四飞轮力矩系数。Wherein, k i4 is the fourth flywheel moment coefficient.

在轨时可以根据第四飞轮转速变化增大或减小摩擦力矩,调整飞轮的标称转速范围。When in orbit, the friction torque can be increased or decreased according to the change of the speed of the fourth flywheel, and the nominal speed range of the flywheel can be adjusted.

步骤二,将控制器输出控制指令限幅在最大输出力矩,并按照第一飞轮、第二飞轮、第三飞轮顺序构成的飞轮分配矩阵分配每个飞轮的控制力矩,分配后将输出力矩最大的飞轮限幅在最大力矩,并将其它两个飞轮的控制力矩成比例缩放,得到三个飞轮控制力矩。Step 2: limit the output control command of the controller to the maximum output torque, and distribute the control torque of each flywheel according to the flywheel distribution matrix formed by the order of the first flywheel, the second flywheel and the third flywheel. The flywheel is limited at the maximum torque, and the control torques of the other two flywheels are scaled proportionally to obtain the control torques of the three flywheels.

步骤二中,第一飞轮、第二飞轮、第三飞轮安装矩阵Cu,将反作用轮的安装矩阵Cu求逆,可得到对应的反作用轮分配矩阵CPIn step 2, the installation matrix C u of the first flywheel, the second flywheel, and the third flywheel is inverted to the installation matrix C u of the reaction wheel, and the corresponding distribution matrix C P of the reaction wheel can be obtained.

首先生成反作用轮控制指令Twc。由反作用轮控制器控制指令Tc,限幅在±0.075Nm(飞轮最大输出力矩)内。First generate the reaction wheel control command T wc . The command T c is controlled by the reaction wheel controller, and the limit is within ±0.075Nm (the maximum output torque of the flywheel).

经过上步求得的分配矩阵CP,可得到用于指令分配的3个飞轮的控制分配指令。Through the distribution matrix C P obtained in the previous step, the control distribution instructions of the three flywheels for instruction distribution can be obtained.

Twc3=CPTc T wc3 = C P T c

将反作用轮的控制指令Twc中取最大值限制在±0.075Nm以内,Limit the maximum value of the control command T wc of the reaction wheel within ±0.075Nm,

再按照比例缩放计算每个飞轮力矩控制指令,即:Then calculate each flywheel torque control command according to the scaling, that is:

Twci=Twci*Twcmax/Twci(max) T wci =T wci *T wcmax /T wci(max)

反作用轮力矩控制电压:Reaction wheel torque control voltage:

Uwci=Twci/kii U wci = T wci /k ii

其中,kii为反作用飞轮力矩系数。Among them, k ii is the reaction flywheel moment coefficient.

步骤三,计算四个反作用轮合成角动量,任意轴角动量达到卸载阈值范围内,并符合卸载时磁场强度和角动量夹角关系时,启动磁力矩器卸载。Step 3: Calculate the synthetic angular momentum of the four reaction wheels. When the angular momentum of any axis reaches the unloading threshold range and conforms to the angular relationship between the magnetic field strength and the angular momentum during unloading, the magnetic torquer is started to unload.

本步骤中,计算四个反作用飞轮合成角动量后,根据角动量大小判断是否需要卸载,判断条件,以其中一轴为例:hx为飞轮合成X轴角动量,卸载判断条件如附图1所示。In this step, after calculating the synthesized angular momentum of the four reaction flywheels, judge whether unloading is required according to the magnitude of the angular momentum, and the judging conditions, taking one of the axes as an example: h x is the synthesized angular momentum of the X-axis of the flywheel, and the judging conditions for unloading are shown in Figure 1 shown.

a.条件:|hi-Hi|≥k1 a. Condition: |h i -H i |≥k 1

记录:IsDumpi=onRecord: IsDump i = on

b.条件:IsDumpi=on且k1>|hi-Hi|≥k2 b. Condition: IsDump i = on and k 1 >|h i -H i |≥k 2

记录:IsDumpi=onRecord: IsDump i = on

c.条件:IsDimpi=off且k1>|hi-Hi|≥k2 c. Condition: IsDimp i =off and k 1 >|h i -H i |≥k 2

记录:IsDimpi=offRecord: IsDimp i = off

d.条件:|hi-Hi|<k2 d. Condition: |h i -H i |<k 2

记录:IsDimpi=offRecord: IsDimp i = off

其中,反作用飞轮的标称角动量Hi,整星零动量飞轮合成标称角动量为0Nms,k1飞轮卸载角动量开阈值,k2飞轮卸载角动量关阈值,IsDumpi,i=X,Y,Z飞轮卸载状态。Among them, the nominal angular momentum H i of the reaction flywheel, the composite nominal angular momentum of the whole star zero momentum flywheel is 0Nms, k 1 flywheel unloading angular momentum opening threshold, k 2 flywheel unloading angular momentum closing threshold, IsDump i , i=X, Y, Z flywheel unloading status.

并且,在本步骤中,高轨弱磁环境下,内源场的磁场强度已经衰减到一百nT左右,与外源场的磁场强度接近。由于外源场模型复杂,受环境影响大,因此采用内源场计算磁场强度,为了避免外源场的影响,导致卸载力矩反向,夹角大于45°小于135°,其中为ΔH和本体系地磁场强度Bb的单位矢量,卸载方向:Moreover, in this step, under the high-orbit weak magnetic environment, the magnetic field strength of the internal source field has decayed to about 100 nT, which is close to the magnetic field strength of the external source field. Due to the complex model of the external source field, which is greatly affected by the environment, the internal source field is used to calculate the magnetic field strength. In order to avoid the influence of the external source field, the unloading torque is reversed. with The included angle is greater than 45° and less than 135°, where with is the unit vector of ΔH and the local geomagnetic field strength B b , unloading direction:

否则,Mc=0,不进行卸载。Otherwise, M c =0, no unloading is performed.

步骤四,选择卸载能力最大的方向上磁力矩器工作,采用开关控制方式,对其它的飞轮进行卸载。Step 4, select the direction with the largest unloading capacity for the magnetic torque device to work, and use the switch control method to unload other flywheels.

本步骤中,卫星在三个方向安装磁力矩器,卸载过程中每个磁力矩器都会产生卸载力矩,但是三个轴其作用不一样,其主要作用占60%以上,最小占10%左右,选用起主要作用两个磁力矩器。磁力矩器产生的力矩Tmi在ΔH上投影产生力矩ΔTmi In this step, the satellite installs magnetic torque devices in three directions. During the unloading process, each magnetic torque device will generate unloading torque, but the functions of the three axes are different. The main function accounts for more than 60%, and the minimum accounts for about 10%. Two magnetic torque devices are selected to play the main role. The torque T mi generated by the magnetic torque device is projected on ΔH to generate the torque ΔT mi

若|Tmi|≠0If |T mi |≠0

否则,ΔTmi=0Otherwise, ΔT mi =0

选择三轴中卸载作用最大两个轴的磁力矩器工作Select the magnetic torque device with the largest unloading effect on the two axes among the three axes

通过比较|ΔTm1|,|ΔTm2|,|ΔTm3|,ΔTmi最大2个轴卸载控制器控制指令为Mci,其他轴控制指令为0。By comparing |ΔT m1 |, |ΔT m2 |, |ΔT m3 |, ΔT mi maximum 2 axis unloading controller control command is M ci , other axis control command is 0.

下面给出上述方案在某高轨卫星上的实施例。Provide the embodiment of above-mentioned scheme on a certain high-orbit satellite below.

卫星上采用三正交一斜装反作用飞轮,飞轮最大角动量为45Nms,转速范围±6000rpm,最大输出力矩0.075Nm,控制电压±10V。四个反作用飞轮安装矩阵为:The satellite adopts a three-orthogonal and one-oblique reaction flywheel with a maximum angular momentum of 45Nms, a speed range of ±6000rpm, a maximum output torque of 0.075Nm, and a control voltage of ±10V. The four reaction flywheel installation matrices are:

第一步:将斜装反作用轮4控制在标称转速附近,其他正交反作用轮控制姿态,根据角动量守恒原理,在斜装轮作用下,四个飞轮都工作在标称转速附近。Step 1: Control the obliquely mounted reaction wheel 4 near the nominal speed, and the other orthogonal reaction wheels control the attitude. According to the principle of conservation of angular momentum, under the action of the obliquely mounted wheel, all four flywheels work near the nominal speed.

反作用轮在轨工作最佳转速在2000rpm左右,设置4个飞轮(代号X,Y,Z,S1)标称角动量分别为12Nms,12Nms,12Nms,20.784Nms,对应转速1596rpm,1596rpm,1596rpm,2764rpm,斜装反作用飞轮S1根据在2764rpm附近摩擦力矩对应控制电压控制飞轮S1工作,如果反作用轮S1转速低于2764rpm,提高控制电压,反之降低控制电压,将反作用飞轮S2标称转速控制在1800rpm~3000rpm附近。其他三个反作用根据控制器输出控制力矩转化成控制电压控制反作用轮X,Y,Z工作,具体计算方法:The optimum rotational speed of the reaction wheel on orbit is around 2000rpm, and the nominal angular momentums of 4 flywheels (code X, Y, Z, S1) are respectively 12Nms, 12Nms, 12Nms, 20.784Nms, corresponding to rotational speeds of 1596rpm, 1596rpm, 1596rpm, 2764rpm , the obliquely installed reaction flywheel S1 controls the work of the flywheel S1 according to the control voltage corresponding to the friction torque near 2764rpm. If the speed of the reaction wheel S1 is lower than 2764rpm, increase the control voltage; nearby. The other three reactions are converted into control voltage according to the output control torque of the controller to control the work of the reaction wheels X, Y, and Z. The specific calculation method:

将正交安装的3个反作用轮的安装矩阵Cu求逆,可得到对应的反作用轮分配矩阵CPThe installation matrix C u of the three reaction wheels installed orthogonally is inverted, and the corresponding distribution matrix C P of the reaction wheels can be obtained.

由反作用轮控制器控制指令Tc,限幅在±0.075Nm内。The command T c is controlled by the reaction wheel controller, and the limit is within ±0.075Nm.

经过上步求得的分配矩阵CP,可得到用于指令分配的3个飞轮的控制分配指令。Through the distribution matrix C P obtained in the previous step, the control distribution instructions of the three flywheels for instruction distribution can be obtained.

Twc=CPTc T wc = C P T c

将反作用轮X/Y/Z的控制指令Twc中取最大值限制在±0.075Nm以内,Limit the maximum value of the control command T wc of the reaction wheel X/Y/Z within ±0.075Nm,

再计算每个飞轮力矩控制指令,即:Twci=Twci*Twcmax/Twci(max) Then calculate each flywheel torque control command, namely: T wci =T wci *T wcmax /T wci(max)

将反作用轮X/Y/Z控制指令Twc转换为力矩指令形式Uwc,用于驱动反作用飞轮工作:Convert the reaction wheel X/Y/Z control command T wc into the torque command form U wc for driving the reaction flywheel:

Uwc=Twc/0.0075U wc = T wc /0.0075

第二步:计算反作用飞轮合成角动量,根据卸载策略计算磁力矩器控制指令如图2所示,驱动磁力矩器工作。The second step: Calculate the synthetic angular momentum of the reaction flywheel, calculate the control command of the magnetic torque device according to the unloading strategy, as shown in Figure 2, and drive the magnetic torque device to work.

依据上述方案,卫星在轨反作用轮转速曲线如图3所示,反作用轮工作在标称转速±500rpm以内,卫星姿态曲线如图4所示,姿态控制精度在0.03度,一个月累计卸载时间29477s,在该控制方案下,整星保持零动量,在高轨空间环境下,干扰力矩小,在该卸载作用下能很好保持飞轮转速。According to the above scheme, the rotation speed curve of the satellite’s in-orbit reaction wheel is shown in Figure 3, the reaction wheel works within the nominal speed of ±500rpm, the satellite attitude curve is shown in Figure 4, the attitude control accuracy is 0.03 degrees, and the cumulative unloading time of one month is 29477s , under this control scheme, the entire satellite maintains zero momentum, and in the high-orbit space environment, the disturbance torque is small, and the flywheel speed can be well maintained under the unloading effect.

本实施例已成功应用于型号卫星,在轨数据表明,采用该飞轮控制策略和卸载方法简单有效,能保证反作用飞轮工作在标称转速附近,且姿态控制精度高,卸载能力足够,且能避开外源场的影响。This embodiment has been successfully applied to model satellites. The in-orbit data show that the flywheel control strategy and unloading method are simple and effective, and can ensure that the reaction flywheel works near the nominal speed, and the attitude control accuracy is high, the unloading capacity is sufficient, and it can avoid Open the influence of the external source field.

以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。The above is only a preferred embodiment of the present invention, it should be pointed out that for those of ordinary skill in the art, without departing from the principle of the present invention, some improvements and modifications can also be made, and these improvements and modifications should also be considered Be the protection scope of the present invention.

Claims (7)

1.一种开环控制反作用轮转速变化及卸载方法,其特征在于,该方法包括如下步骤:1. an open-loop control reaction wheel speed variation and unloading method, it is characterized in that, the method comprises the steps: 步骤一,采用包括第一飞轮、第二飞轮、第三飞轮和第四飞轮的四个飞轮进行控制,并将第四飞轮的飞轮力矩控制指令大小调整至与理论标称转速对应的摩擦力矩一致,其中,飞轮即为反作用轮;Step 1: Use four flywheels including the first flywheel, the second flywheel, the third flywheel and the fourth flywheel for control, and adjust the size of the flywheel torque control command of the fourth flywheel to be consistent with the friction torque corresponding to the theoretical nominal speed , where the flywheel is the reaction wheel; 步骤二,将控制器输出控制指令限幅在最大输出力矩,并按照第一飞轮、第二飞轮、第三飞轮顺序构成的飞轮分配矩阵分配每个飞轮的控制力矩,分配后将输出力矩最大的飞轮限幅在最大输出力矩,并将其它两个飞轮的控制力矩成比例缩放,得到三个飞轮控制力矩;Step 2: limit the output control command of the controller to the maximum output torque, and distribute the control torque of each flywheel according to the flywheel distribution matrix formed by the order of the first flywheel, the second flywheel and the third flywheel. The flywheel is limited to the maximum output torque, and the control torque of the other two flywheels is scaled proportionally to obtain the control torque of the three flywheels; 步骤三,计算四个飞轮合成角动量,任意轴角动量达到卸载阈值范围内,并符合卸载时磁场强度和角动量夹角关系时,启动磁力矩器卸载;Step 3: Calculating the composite angular momentum of the four flywheels. When the angular momentum of any axis reaches the unloading threshold range and meets the angular relationship between the magnetic field strength and the angular momentum during unloading, the magnetic torquer is started to unload; 步骤四,选择卸载能力最大的方向上的磁力矩器工作,采用开关控制方式,对飞轮进行卸载。Step 4, select the magnetic torque device in the direction with the largest unloading capacity to work, and use the switch control method to unload the flywheel. 2.根据权利要求1所述的开环控制反作用轮转速变化及卸载方法,其特征在于:2. The open-loop control reaction wheel speed change and unloading method according to claim 1, characterized in that: 步骤一中,根据第四飞轮标称转速得到摩擦力矩Tf4,对其进行摩擦力矩补偿,得到飞轮力矩控制指令Twc4 In step 1, the friction torque T f4 is obtained according to the nominal speed of the fourth flywheel, and the friction torque is compensated to obtain the flywheel torque control command T wc4 Twc4=-Tf4 T wc4 = -T f4 和飞轮力矩控制电压:and flywheel torque control voltage: Uwc4=Twc4/ki4 U wc4 =T wc4 /k i4 其中,ki4为第四飞轮力矩系数。Wherein, k i4 is the fourth flywheel moment coefficient. 3.根据权利要求1所述的开环控制反作用轮转速变化及卸载方法,其特征在于:步骤二中,第一飞轮、第二飞轮、和第三飞轮的安装矩阵Cu,将其求逆可得到对应的飞轮分配矩阵CP3. The open-loop control reaction wheel speed change and unloading method according to claim 1, characterized in that: in step 2, the installation matrix C u of the first flywheel, the second flywheel, and the third flywheel is inversely calculated The corresponding flywheel allocation matrix C P can be obtained. 4.根据权利要求1所述的开环控制反作用轮转速变化及卸载方法,其特征在于,步骤三中,计算四个飞轮合成角动量后,根据角动量大小判断是否需要卸载。4. The open-loop control reaction wheel speed change and unloading method according to claim 1, characterized in that in step 3, after calculating the synthetic angular momentum of the four flywheels, it is judged whether unloading is required according to the angular momentum. 5.根据权利要求1所述的开环控制反作用轮转速变化及卸载方法,其特征在于,步骤三中采用内源场计算磁场强度。5. The open-loop control reaction wheel speed change and unloading method according to claim 1, characterized in that in step 3, the internal source field is used to calculate the magnetic field strength. 6.根据权利要求1所述的开环控制反作用轮转速变化及卸载方法,其特征在于,步骤四中选用起主要作用的磁力矩器进行卸载工作。6. The open-loop control reaction wheel speed change and unloading method according to claim 1, characterized in that, in step 4, the magnetic torque device that plays a major role is selected to perform the unloading work. 7.根据权利要求1所述的开环控制反作用轮转速变化及卸载方法,其特征在于,第一飞轮、第二飞轮和第三飞轮为正交安装的飞轮,第四飞轮为斜装飞轮。7. The open-loop control reaction wheel speed change and unloading method according to claim 1, characterized in that the first flywheel, the second flywheel and the third flywheel are orthogonally mounted flywheels, and the fourth flywheel is a diagonally mounted flywheel.
CN201610081918.2A 2016-02-05 2016-02-05 A kind of opened loop control reaction wheel rotation speed change and discharging method Expired - Fee Related CN105644810B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610081918.2A CN105644810B (en) 2016-02-05 2016-02-05 A kind of opened loop control reaction wheel rotation speed change and discharging method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610081918.2A CN105644810B (en) 2016-02-05 2016-02-05 A kind of opened loop control reaction wheel rotation speed change and discharging method

Publications (2)

Publication Number Publication Date
CN105644810A CN105644810A (en) 2016-06-08
CN105644810B true CN105644810B (en) 2017-09-12

Family

ID=56488332

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610081918.2A Expired - Fee Related CN105644810B (en) 2016-02-05 2016-02-05 A kind of opened loop control reaction wheel rotation speed change and discharging method

Country Status (1)

Country Link
CN (1) CN105644810B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106218922B (en) * 2016-07-27 2018-06-15 中国科学院长春光学精密机械与物理研究所 The joint actuating mechanism controls method of flexible agility satellite
CN106542120B (en) * 2016-09-30 2018-11-02 上海航天控制技术研究所 In conjunction with the satellite three-axis attitude control method of magnetic torquer when flywheel drive lacking
CN106542118B (en) * 2016-10-08 2018-08-28 上海航天控制技术研究所 A method of normal attitude being restored to from magnet-wheel joint control state control using flywheel
CN106742071B (en) * 2017-01-20 2019-01-25 上海航天控制技术研究所 A method of zero passage is avoided using the control of orthogonal installation flywheel
CN107600463B (en) * 2017-08-22 2019-05-24 长光卫星技术有限公司 A kind of agility small satellite attitude fast reserve control method
CN108528765B (en) * 2018-04-08 2019-07-05 上海微小卫星工程中心 A method for determining the unloading threshold of low-orbit satellites based on a composite mechanism
CN108983637A (en) * 2018-09-04 2018-12-11 上海微小卫星工程中心 A Semi-Physical Simulation Method of Satellite Control System Using Reaction Wheels
CN109080856B (en) * 2018-10-10 2020-03-17 上海微小卫星工程中心 Method for maintaining rotation speed of reaction wheel
CN109649690B (en) * 2018-12-27 2021-07-09 上海航天控制技术研究所 High-precision hierarchical saturation attitude control method and system for satellite based on flywheel configuration
CN110127088B (en) * 2019-05-05 2020-11-20 北京控制工程研究所 A double-hysteresis loop method for unloading satellite synthetic angular momentum using a magnetic torquer
CN111897355B (en) * 2020-08-06 2022-09-13 中国科学院微小卫星创新研究院 Satellite attitude maneuver trajectory planning method
CN113501144B (en) * 2021-07-27 2022-05-03 北京微纳星空科技有限公司 Flywheel rotating speed control method, device, equipment and storage medium

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5201833A (en) * 1991-07-19 1993-04-13 General Electric Company Attitude control system with reaction wheel friction compensation
JP3281061B2 (en) * 1992-10-14 2002-05-13 株式会社東芝 Spacecraft unloading device
US7198232B1 (en) * 2004-08-26 2007-04-03 Lockheed Martin Corporation Optimal speed management for reaction wheel control system and method
CN100439204C (en) * 2007-07-03 2008-12-03 北京控制工程研究所 Estimation and Compensation Method for Reducing the Influence of Attitude Control Jet on Orbit
CN102343985B (en) * 2011-07-08 2013-07-24 北京航空航天大学 Satellite time optimal posture maneuvering method with reaction flywheel
CN102616386B (en) * 2012-03-29 2014-04-02 哈尔滨工业大学 Uniaxial quick maneuverable spacecraft flywheel configuration and optimization method thereof

Also Published As

Publication number Publication date
CN105644810A (en) 2016-06-08

Similar Documents

Publication Publication Date Title
CN105644810B (en) A kind of opened loop control reaction wheel rotation speed change and discharging method
Yoon et al. Spacecraft adaptive attitude and power tracking with variable speed control moment gyroscopes
Hu et al. Robust finite-time control allocation in spacecraft attitude stabilization under actuator misalignment
Zhang et al. Fuzzy-PID control for the position loop of aerial inertially stabilized platform
CN106933241B (en) Single-gimbal control momentum gyro spacecraft fault tolerant control method based on fault de couple
CN106292681A (en) A kind of satellite Active Fault-tolerant Control Method distributed based on observer and On-line Control
CN107438806B (en) Method and apparatus for controlling the attitude of a spacecraft
CN103034237B (en) Spacecraft attitude maneuver control method using two single-frame control moment gyroscopes
CN109533396B (en) Satellite spin orientation method based on magnetic measurement and control
Huang et al. Optimal Lorentz-augmented spacecraft formation flying in elliptic orbits
US8014911B2 (en) Methods and systems for imposing a momentum boundary while reorienting an agile vehicle with control moment gyroscopes
CN109625329A (en) A kind of autonomous discharging method of flywheel angular momentum based on discrete jet
US9567112B1 (en) Method and apparatus for singularity avoidance for control moment gyroscope (CMG) systems without using null motion
CN109677638B (en) An Improved Pure Magnetron Spin-to-Sun Orientation Method Based on Geomagnetic Field Measurement Parameters
CN113955152B (en) Star sun-to-day directional control method
JP2009173069A (en) Posture control device for artificial satellites, and posture control method
CN111897352B (en) Attitude Maneuvering Dynamic Optimization Method and Medium Based on Control Torque Output Constraints
JP2022121605A (en) Satellite attitude control system using eigen vector, non-linear dynamic inversion, and feedforward control
Zidek et al. Receding horizon drift counteraction and its application to spacecraft attitude control
CN111891401B (en) CMG group return nominal configuration control method, system and medium based on zero-motion optimization
CN109911250B (en) Method for magnetically measuring magnetic control rate damping
CN108528765B (en) A method for determining the unloading threshold of low-orbit satellites based on a composite mechanism
Karami et al. Spacecraft momentum dumping using fewer than three external control torques
JP7416396B2 (en) Attitude control device
Lee et al. Balancing control of a single-wheel robot considering power-efficiency and gyroscopic instability suppression

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
CB03 Change of inventor or designer information

Inventor after: Bai Tao

Inventor after: Xiong Shujie

Inventor after: Wu Guoqiang

Inventor after: Lin Baojun

Inventor after: Qin Guijun

Inventor after: Yan Yanteng

Inventor after: Wang Haoguang

Inventor before: Bai Tao

Inventor before: Xiong Shujie

Inventor before: Wu Guoqiang

Inventor before: Qin Guijun

Inventor before: Yan Yanteng

CB03 Change of inventor or designer information
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170912

Termination date: 20180205

CF01 Termination of patent right due to non-payment of annual fee