CN105667838B - 一种皮纳卫星的模块化姿态确定与控制装置及其方法 - Google Patents
一种皮纳卫星的模块化姿态确定与控制装置及其方法 Download PDFInfo
- Publication number
- CN105667838B CN105667838B CN201610143290.4A CN201610143290A CN105667838B CN 105667838 B CN105667838 B CN 105667838B CN 201610143290 A CN201610143290 A CN 201610143290A CN 105667838 B CN105667838 B CN 105667838B
- Authority
- CN
- China
- Prior art keywords
- max
- sign
- yes
- satellite
- judge whether
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000006340 racemization Effects 0.000 claims abstract description 16
- 238000012545 processing Methods 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 30
- 230000008569 process Effects 0.000 claims description 10
- 230000000087 stabilizing effect Effects 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 3
- 230000002441 reversible effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 230000005611 electricity Effects 0.000 claims 1
- 238000013461 design Methods 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 108091092878 Microsatellite Proteins 0.000 description 1
- 206010034972 Photosensitivity reaction Diseases 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 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
- B64G1/36—Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors
- B64G1/369—Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors using gyroscopes as attitude sensors
-
- 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/32—Guiding or controlling apparatus, e.g. for attitude control using earth's magnetic field
-
- 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/36—Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors
- B64G1/363—Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors using sun sensors
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
Abstract
本发明属于皮纳卫星技术领域,公开了一种皮纳卫星的模块化姿态确定与控制装置及其方法。装置包括外壳,外壳内设置有用于测定皮纳卫星姿态的姿态敏感单元,用于根据皮纳卫星姿态发出命令信号的姿控计算机,用于执行姿控计算机发出的命令的姿态控制执行单元;方法为:卫星本体坐标系为Obxyz,轨道坐标系Ooxyz,太阳矢量和地磁场矢量在Ooxyz中分别为S1、S2,在Obxyz中分别为在姿控计算机内预存有皮纳卫星的旋转角度的设置阈值Δw,陀螺仪实时测量皮纳卫星的旋转角度w,旋转角度w大于姿控计算机内的设置阈值Δw,启动姿态消旋控制对皮纳卫星进行消旋处理,直至旋转角度w小于等于姿控计算机内的设置阈值Δw后,进入姿态确定和姿态捕获控制的循环阶段。
Description
技术领域
本发明属于皮纳卫星技术领域,特别涉及一种皮纳卫星的模块化姿态确定与控制装置及其方法。
背景技术
皮纳卫星领域飞速发展,各单位、高校争相进入皮纳卫星领域,作为一个新兴产业,它拥有广阔的市场前景。在单位企业里,微小卫星以其实现成本低,设计时间周期短等的快速响应能力应用在应急通信、空间组网等领域。在高校教育中,其低成本、高集成的特性使得皮纳卫星可以很好的融入创新与学习,对学生的发展很有帮助,所以被很多大学所重视。当下已有多所大学的小卫星发射成功。
皮纳卫星的姿态确定与控制***(Attitude Determination and ControlSystem简写为ADCS)为小卫星提供姿态控制能力,它由姿态敏感器、姿态控制执行器与姿控计算机组成。姿态敏感器配合姿控计算机可以提供卫星姿态信息,通过姿控计算机对陀螺仪、卫星导航接收机(GNSS)和磁强计、太阳敏感器、星敏感器等姿态敏感器采集的信息进行处理,解算出当前卫星的姿态,为姿态控制提供保障。姿态控制执行器是姿态控制算法的执行机构,主要包括反作用飞轮、磁力矩器和微推进器等,通过姿控计算机中运行的姿态控制算法真正实现小卫星的姿态控制。姿控计算机提供姿态确定、姿态控制以及数据存储功能,皮纳卫星种类多样,外形多样,结构多样,根据用户的需求设计不同的皮纳卫星上。
现在已经研制成功并试飞成功的一些皮纳卫星,由于体积和功耗相比于100kg级小卫星又降低了1到2个数量级,所以只有很少部分拥有姿态确定与控制功能,拥有ADCS的,也都是针对自己研制的小卫星,量身定做,集成在***内部,重新研制小卫星时需要再次研制其ADCS,没有移植性,并且无法接入其他单位制作的小卫星中。即使有可以应用在其他卫星上的ADCS,也必须是经过硬件修改之后才可以。归结就是没有通用的接口设计,或者没有自己的处理器,需要其他使用卫星的处理器,这样对星务软件的独立设计又提出了巨大要求。这些ADCS很难在着眼于皮纳卫星其他领域应用的企业和高校中推广,无法为其制造便捷的解决方案。这就使得皮纳卫星的研制周期加长,影响研究课题的进度。现在国内没有针对皮纳卫星的模块化通用ADCS产品。
发明内容
本发明的目的是提供一种皮纳卫星的模块化姿态确定与控制装置及其方法,该装置为模块化设计提高了其通用性。
为达到以上目的,本发明采用以下技术方案予以实现。
方案一:
一种皮纳卫星的模块化姿态确定与控制装置,其特征在于:包括用于测定皮纳卫星姿态的姿态敏感单元,用于根据皮纳卫星姿态发出命令信号的姿控计算机,用于执行姿控计算机发出的命令的姿态控制执行单元;
所述姿态敏感单元包含用于测量皮纳卫星在惯性系下的角速度的陀螺仪,用于测量卫星本体坐标系下地磁场矢量的磁强计,用于测量卫星本体坐标系下太阳矢量的太阳敏感器,用于测量皮纳卫星的位置信息的卫星导航接收机;所述姿控计算机的四个输入端分别对应电连接所述陀螺仪的输出端、卫星导航接收机的输出端、磁强计的输出端、太阳敏感器的输出端;
所述姿态控制执行单元包含磁力矩器驱动单元,用于完成对皮纳卫星实时姿态控制的磁力矩器;所述姿控计算机的输出端电连接磁力矩器驱动单元的输入端,所述磁力矩器驱动单元的输出端电连接磁力矩器的控制端。
上述技术方案的特点和进一步改进:
进一步的,所述姿态敏感单元还包含温度传感器和模拟量电流传感器,所述温度传感器的输出端和模拟量电流传感器的输出端对应电连接所述姿控计算机的两个输入端。
进一步的,所述姿态敏感单元、姿控计算机和姿态控制执行单元设置在外壳内,所述外壳的底部设置为开口状,所述外壳的下底面设置有向外延伸的折弯面,所述折弯面上设置有多个通孔,所述外壳的侧壁上设置有供电接口。
进一步的,所述外壳内还设置有5V升压模块和3.3V稳压模块,所述供电接口与5V升压模块电连接,所述3.3V稳压模块与5V升压模块电连接。
进一步的,所述供电接口为J30J-9ZKW-J接口,所述磁力矩器为三轴磁力矩器,所述姿控计算机采用单片微型计算机。
方案二:
一种皮纳卫星的姿态确定与控制方法,基于上述皮纳卫星的模块化姿态确定与控制装置,其特征在于:包括以下步骤,
步骤一,定义卫星本体坐标系为Obxyz,参考坐标系为轨道坐标系Ooxyz,太阳矢量和地磁场矢量在Ooxyz中分别为S1、S2,太阳矢量和地磁场矢量在Obxyz中分别为
步骤二,皮纳卫星进入任务设定轨道后,在姿控计算机内预存有皮纳卫星的旋转角度的设置阈值Δw,陀螺仪实时测量皮纳卫星的旋转角度w,旋转角度w大于姿控计算机内的设置阈值Δw,启动姿态消旋控制对皮纳卫星进行消旋处理,直至旋转角度w小于等于姿控计算机内的设置阈值Δw后,进入姿态确定和姿态捕获控制的循环阶段。
上述技术方案的特点和进一步改进:
(1)进一步的,所述消旋处理通过B-dot方法对皮纳卫星进行速率阻尼,B-dot方法的计算式为:
式中,Mb为卫星输出磁矩,
K为消旋控制率,
Bb为地磁场矢量在卫星本体坐标系中的投影,通过磁强计测量得到,为Bb的变化率,
则,磁力矩器输出消旋力矩对皮纳卫星进行消旋处理,磁力矩器输出的消旋力矩为:
(2)进一步的,所述姿态确定的过程为,
设s3=s1×s2,令s1、s2、s3三个矢量组成矩阵s的三列S=(S1 S2S3),
令矩阵s*为
则有:
Ao→b为轨道坐标系到本体坐标系的姿态旋转矩阵,
S、S*不共线,矩阵S可逆,Ao→b=S*S-1,
根据姿态动力学中姿态旋转矩阵Ao→b与欧拉角的转换关系反推欧拉角的值为:
式中,a21代表姿态旋转矩阵Ao→b当中2行1列的元素,a11代表姿态旋转矩阵Ao→b当中1行1列的元素,a31代表姿态旋转矩阵Ao→b当中3行1列的元素,a32代表姿态旋转矩阵Ao→b当中3行2列的元素,a33代表姿态旋转矩阵Ao→b当中3行3列的元素;从而完成姿态确定。
(3)进一步的,所述姿态捕获控制采用PD控制,
设期望磁控力矩Tc:
式中:Kp为比例参数,
Kd为微分参数,
θ为余仰角,θ=90°-δ,δ为卫星的赤纬;
为余仰角速度,
期望磁控力矩Tc的三轴分量为Tcx、Tcy、Tcz,将期望磁控力矩Tc的三轴分量Tcx、Tcy、Tcz大小按照顺序排列有如下六种排列方式:
Tcx>Tcy>Tcz、Tcx>Tcz>Tcy、Tcy>Tcx>Tcz、Tcy>Tcz>Tcx、Tcz>Tcx>Tcy、Tcz>Tcy>Tcx,
通过公式Tc=cross(M(i),Bb),cross表示叉乘,在不同排列方式下可以分别解得一组卫星输出磁矩M(i):
M(i)=[mxi、myi、mzi],
其中,i=1、2、3、4、5、6,
mxi为卫星输出磁矩M(i)在卫星本体坐标系中x轴上的分量,
myi为卫星输出磁矩M(i)在卫星本体坐标系中y轴上的分量,
mzi为卫星输出磁矩M(i)在卫星本体坐标系中z轴上的分量,
则:Tb(i)=M(i)×Bb,
式中,Tb(i)为每一种可能性下的实际能输出的力矩;
设dis(i)=||Tc|-|Tb(i)||,
使得dis(i)最小的卫星输出磁矩M(i)为最优卫星输出磁矩M,
则,磁力矩器输出姿态捕获控制力矩对皮纳卫星进行姿态捕获控制处理,磁力矩器输出的姿态捕获控制力矩Tb为:
Tb=cross(M,B'b),
其中,B’b为卫星本体坐标系中地磁场矢量量测值。
(4)进一步的,在六种不同排列方式下,所述卫星输出磁矩M(i)的计算过程如下:
Bb的三轴分量为Bbx、Bby、Bbz,
Tcy>Tcx>Tcz时:
在|Bbx|>|Bbz|的状态下,设mz=Tcy/Bbx,然后判断是否|mz|>mmax,若|mz|>mmax为否,则mz=mz且mx=0,若|mz|>mmax为是,则设mz=sign(mz)mmax且mx=-(Tcx-mzBbx)/Bbz;再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
在|Bbx|≤|Bbz|的状态下,设mx=-Tcy/Bbz,然后判断是否|mx|>mmax,若|mx|>mmax为否,则mx=mx且mz=0,若|mx|>mmax为是,则设mx=sign(mx)mmax且mz=(Tcy-mxBbz)/Bbx,再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
然后计算my,再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
Tcz>Tcy>Tcx时:
在|Bby|>|Bbx|的状态下,设mx=Tcz/Bby,然后判断是否|mx|>mmax,若|mx|>mmax为否,则mx=mx且my=0,若|mx|>mmax为是,则设mx=sign(mx)mmax且my=-(Tcz-mxBby)/Bbx;再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
在|Bby|≤|Bbx|的状态下,设my=-Tcz/Bbx,然后判断是否|my|>mmax,若|my|>mmax为否,则my=my且mx=0,若|my|>mmax为是,则设
my=sign(my)mmax且mx=(Tcz-myBbx)/Bby,再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
然后计算mz,再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
Tcy>Tcz>Tcx时:
在|Bbx|>|Bbz|的状态下,设mz=Tcy/Bbx,然后判断是否|mz|>mmax,若|mz|>mmax为否,则mz=mz且mx=0,若|mz|>mmax为是,则设mz=sign(mz)mmax且mx=-(Tcx-mzBbx)/Bbz;再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
在|Bbx|≤|Bbz|的状态下,设mx=-Tcy/Bbz,然后判断是否|mx|>mmax,若|mx|>mmax为否,则mx=mx且mz=0,若|mx|>mmax为是,则设mx=sign(mx)mmax且mz=(Tcy+mxBbz)/Bbx,再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
然后计算my,再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
Tcz>Tcx>Tcy时:
在|Bby|>|Bbx|的状态下,设mx=Tcz/Bby,然后判断是否|mx|>mmax,若|mx|>mmax为否,则mx=mx且my=0,若|mx|>mmax为是,则设mx=sign(mx)mmax且my=-(Tcz-mxBby)/Bbx;再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
在|Bby|≤|Bbx|的状态下,设my=-Tcz/Bbx,然后判断是否|my|>mmax,若|my|>mmax为否,则my=my且mx=0,若|my|>mmax为是,则设my=sign(my)mmax且mx=(Tcz-myBbx)/Bby,再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
然后计算mz,再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
Tcx>Tcy>Tcz时:
在|Bby|>|Bbz|的状态下,设my=Tcx/Bbz,然后判断是否|my|>mmax,若|my|>mmax为否,则my=my且mz=0,若|my|>mmax为是,则设my=sign(my)mmax且mz=-(Tcx-myBbz)/Bby;再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
在|Bby|≤|Bbz|的状态下,设mz=-Tcx/Bby,然后判断是否|mz|>mmax,若|mz|>mmax为否,则mz=mz且my=0,若|mz|>mmax为是,则设mz=sign(mz)mmax且my=(Tcx+mzBby)/Bbz,再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
然后计算mx,再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
Tcx>Tcz>Tcy时:
在|Bby|>|Bbz|的状态下,设my=Tcx/Bbz,然后判断是否|my|>mmax,若|my|>mmax为否,则my=my且mz=0,若|my|>mmax为是,则设my=sign(my)mmax且mz=-(Tcx-myBbz)/Bby;再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
在|Bby|≤|Bbz|的状态下,设mz=-Tcx/Bby,然后判断是否|mz|>mmax,若|mz|>mmax为否,则mz=mz且my=0,若|mz|>mmax为是,则设mz=sign(mz)mmax且my=(Tcx+mzBby)/Bbz,再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
然后计算mx,再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
其中,mmax=NIS,
式中,N表示线圈匝数,I表示线圈内通过的电流,S表示线圈围成的面积。
本发明的有益效果是:
1、高集成、模块化,可移植性:由于本模块完全封装,集成度高,为皮纳卫星提供最快速的姿态确定于控制能力,模块化、可移植,通过预留的打孔位置,可以满足很多皮纳卫星的需要。
2、在皮纳卫星***中相对独立,仅通过供电即可完成对皮纳卫星的控制。
3、接口简单,本产品能够实时将有关姿态的控制指令与量测数据通过RS232串口传输,方便用户使用,可市场化。
4、本发明的姿态确定与控制方法针对皮纳卫星纯磁控的高非线性特点,采用了专家控制这样一种先进的解决非线性控制问题的方法。并且在5000秒内能完成皮纳卫星的姿态确定、姿态控制与姿态捕获。
附图说明
图1为本发明的一种皮纳卫星的模块化姿态确定与控制装置的电气连接结构示意图;
图2为本发明的一种皮纳卫星的模块化姿态确定与控制装置中的外壳的结构示意图;
图3为图2的俯视示意图;
图4为计算最优卫星输出磁矩M-best的流程图;
图5a、图5b、图5c、图5d、图5e、图5f为Tcx>Tcy>Tcz、Tcx>Tcz>Tcy、Tcy>Tcx>Tcz、Tcy>Tcz>Tcx、Tcz>Tcx>Tcy、Tcz>Tcy>Tcx六种不同排列方式下卫星输出磁矩M(i)的计算流程图;
图中:1、外壳;2、折弯面、3、供电接口。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
参照图1,为本发明的一种皮纳卫星的模块化姿态确定与控制装置的电气连接结构示意图;该皮纳卫星的模块化姿态确定与控制装置包括外壳1,外壳1内设置有用于测定皮纳卫星姿态的姿态敏感单元,用于根据皮纳卫星姿态发出命令信号的姿控计算机,用于执行姿控计算机发出的命令的姿态控制执行单元。
姿态敏感单元包含用于测量皮纳卫星在惯性系下的角速度的陀螺仪,用于测量卫星本体坐标系下地磁场矢量的磁强计,用于测量卫星本体坐标系下太阳矢量的太阳敏感器,用于测量皮纳卫星的位置信息的卫星导航接收机;姿控计算机的四个输入端分别对应电连接陀螺仪的输出端、卫星导航接收机的输出端、磁强计的输出端、太阳敏感器的输出端。姿态控制执行单元包含磁力矩器驱动单元,用于完成对皮纳卫星实时姿态控制的磁力矩器;姿控计算机的输出端电连接磁力矩器驱动单元的输入端,磁力矩器驱动单元的输出端电连接磁力矩器的控制端。
参照图2、图3,为外壳的结构示意图;姿态敏感单元、姿控计算机和姿态控制执行单元设置在外壳内,将所有的器件均设置在外壳1内,达到模块化设计,使用时直接安装该模块就可以,为研发团队节省了人力物力;外壳1的底部设置为开口状,外壳1的下底面设置有向外延伸的折弯面2,折弯面2上设置有多个通孔,通孔用于通过螺钉安装该装置,方便安装,外壳1的侧壁上设置有供电接口3。
姿态敏感单元还包含温度传感器和模拟量电流传感器,温度传感器的输出端和模拟量电流传感器的输出端对应电连接姿控计算机的两个输入端。温度传感器、电流传感器为辅助敏感器,在算法中不起决定性作用,只是用来修饰数据。
磁力矩器为三轴磁力矩器。姿控计算机采用单片微型计算机。供电接口为J30J-9ZKW-J接口。
外壳1内还设置有5V升压模块和3.3V稳压模块,供电接口与5V升压模块电连接,3.3V稳压模块与5V升压模块电连接。
该皮纳卫星的姿态确定与控制装置的姿态确定与控制方法,包括以下步骤:步骤一,定义卫星本体坐标系为Obxyz,参考坐标系为轨道坐标系Ooxyz,太阳矢量和地磁场矢量在Ooxyz中分别为S1、S2,太阳矢量和地磁场矢量在Obxyz中分别为
步骤二,皮纳卫星进入任务设定轨道后,在姿控计算机内预存有皮纳卫星的旋转角度的设置阈值Δw,陀螺仪实时测量皮纳卫星的旋转角度w,旋转角度w大于姿控计算机内的设置阈值Δw,启动姿态消旋控制对皮纳卫星进行消旋处理,直至旋转角度w小于等于姿控计算机内的设置阈值Δw后,进入姿态确定和姿态捕获控制的循环阶段。
消旋处理通过B-dot方法对皮纳卫星进行速率阻尼。
B-dot方法的计算式为:
式中,Bb为地磁场矢量在卫星本体坐标系中的投影,通过磁强计测量得到,
为Bb的变化率,
Mb为卫星输出磁矩,
K为消旋控制率,B-dot算法控制率的取值范围很广,可以在1000-1000000范围内取值,取值大小主要影响算法收敛速度,如果控制率取值过大,产生的磁距大于线圈实际可输出的最大磁距,影响算法效率。
则,磁力矩器输出消旋力矩对皮纳卫星进行消旋处理,磁力矩器输出的消旋力矩为:
姿态确定的过程为:
设s3=s1×s2,令s1、s2、s3三个矢量组成矩阵s的三列S=(S1 S2S3),
令矩阵s*为
则有:
Ao→b为轨道坐标系到本体坐标系的姿态旋转矩阵,
S、S*不共线,矩阵S可逆,Ao→b=S*S-1,
根据姿态动力学中姿态旋转矩阵Ao→b与欧拉角的转换关系反推欧拉角的值为:
式中,a21代表姿态旋转矩阵Ao→b当中2行1列的元素,a11代表姿态旋转矩阵Ao→b当中1行1列的元素,a31代表姿态旋转矩阵Ao→b当中3行1列的元素,a32代表姿态旋转矩阵Ao→b当中3行2列的元素,a33代表姿态旋转矩阵Ao→b当中3行3列的元素;从而完成姿态确定。
姿态捕获控制采用PD控制,具体方法如下:
设期望磁控力矩Tc:
式中:θ为余仰角,θ=90°-δ,δ为卫星的赤纬,
为余仰角速度,
Kp为比例参数,
Kd为微分参数,PD算法中参数取值的依据是最终控制产生的磁距是否在00.2以内,同时角速度控制率和角度控制率存在比例关系,一般角速度控制率为角度控制率的10的三次方倍。
期望磁控力矩Tc的三轴分量为Tcx、Tcy、Tcz,将期望磁控力矩Tc的三轴分量Tcx、Tcy、Tcz大小按照顺序排列有如下六种排列方式:
Tcx>Tcy>Tcz、Tcx>Tcz>Tcy、Tcy>Tcx>Tcz、Tcy>Tcz>Tcx、Tcz>Tcx>Tcy、Tcz>Tcy>Tcx,
通过公式Tc=cross(M(i),Bb),cross表示叉乘,在不同排列方式下可以分别解得一组卫星输出磁矩M(i):
M(i)=[mxi、myi、mzi],
其中,i=1、2、3、4、5、6,
mxi为卫星输出磁矩M(i)在卫星本体坐标系中x轴上的分量,
myi为卫星输出磁矩M(i)在卫星本体坐标系中y轴上的分量,
mzi为卫星输出磁矩M(i)在卫星本体坐标系中z轴上的分量,
如图5a、图5b、图5c、图5d、图5e、图5f所示,在六种不同排列方式下,卫星输出磁矩M(i)的计算过程如下:
Bb的三轴分量为Bbx、Bby、Bbz,
Tcy>Tcx>Tcz时:
在|Bbx|>|Bbz|的状态下,设mz=Tcy/Bbx,然后判断是否|mz|>mmax,若|mz|>mmax为否,则mz=mz且mx=0,若|mz|>mmax为是,则设mz=sign(mz)mmax且mx=-(Tcx-mzBbx)/Bbz;再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
在|Bbx|≤|Bbz|的状态下,设mx=-Tcy/Bbz,然后判断是否|mx|>mmax,若|mx|>mmax为否,则mx=mx且mz=0,若|mx|>mmax为是,则设mx=sign(mx)mmax且mz=(Tcy-mxBbz)/Bbx,再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
然后计算my,再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
Tcz>Tcy>Tcx时:
在|Bby|>|Bbx|的状态下,设mx=Tcz/Bby,然后判断是否|mx|>mmax,若|mx|>mmax为否,则mx=mx且my=0,若|mx|>mmax为是,则设mx=sign(mx)mmax且my=-(Tcz-mxBby)/Bbx;再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
在|Bby|≤|Bbx|的状态下,设my=-Tcz/Bbx,然后判断是否|my|>mmax,若|my|>mmax为否,则my=my且mx=0,若|my|>mmax为是,则设my=sign(my)mmax且mx=(Tcz-myBbx)/Bby,再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
然后计算mz,再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
Tcy>Tcz>Tcx时:
在|Bbx|>|Bbz|的状态下,设mz=Tcy/Bbx,然后判断是否|mz|>mmax,若|mz|>mmax为否,则mz=mz且mx=0,若|mz|>mmax为是,则设mz=sign(mz)mmax且mx=-(Tcx-mzBbx)/Bbz;再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
在|Bbx|≤|Bbz|的状态下,设mx=-Tcy/Bbz,然后判断是否|mx|>mmax,若|mx|>mmax为否,则mx=mx且mz=0,若|mx|>mmax为是,则设mx=sign(mx)mmax且mz=(Tcy+mxBbz)/Bbx,再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
然后计算my,再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
Tcz>Tcx>Tcy时:
在|Bby|>|Bbx|的状态下,设mx=Tcz/Bby,然后判断是否|mx|>mmax,若|mx|>mmax为否,则mx=mx且my=0,若|mx|>mmax为是,则设mx=sign(mx)mmax且my=-(Tcz-mxBby)/Bbx;再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
在|Bby|≤|Bbx|的状态下,设my=-Tcz/Bbx,然后判断是否|my|>mmax,若|my|>mmax为否,则my=my且mx=0,若|my|>mmax为是,则设my=sign(my)mmax且mx=(Tcz-myBbx)/Bby,再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
然后计算mz,再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
Tcx>Tcy>Tcz时:
在|Bby|>|Bbz|的状态下,设my=Tcx/Bbz,然后判断是否|my|>mmax,若|my|>mmax为否,则my=my且mz=0,若|my|>mmax为是,则设my=sign(my)mmax且mz=-(Tcx-myBbz)/Bby;再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
在|Bby|≤|Bbz|的状态下,设mz=-Tcx/Bby,然后判断是否|mz|>mmax,若|mz|>mmax为否,则mz=mz且my=0,若|mz|>mmax为是,则设mz=sign(mz)mmax且my=(Tcx+mzBby)/Bbz,再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
然后计算mx,再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
Tcx>Tcz>Tcy时:
在|Bby|>|Bbz|的状态下,设my=Tcx/Bbz,然后判断是否|my|>mmax,若|my|>mmax为否,则my=my且mz=0,若|my|>mmax为是,则设my=sign(my)mmax且mz=-(Tcx-myBbz)/Bby;再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
在|Bby|≤|Bbz|的状态下,设mz=-Tcx/Bby,然后判断是否|mz|>mmax,若|mz|>mmax为否,则mz=mz且my=0,若|mz|>mmax为是,则设mz=sign(mz)mmax且my=(Tcx+mzBby)/Bbz,再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
然后计算mx,再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
其中,mmax=NIS,指的是线圈可以产生的最大磁矩,取决于供电电流、线圈围成的面积和线圈匝数,
式中,N表示线圈匝数,I表示线圈内通过的电流,S表示线圈围成的面积。
则:Tb(i)=M(i)×Bb,
式中,Tb(i)为每一种可能性下的实际能输出的力矩;
设dis(i)=||Tc|-|Tb(i)||
使得dis(i)最小的卫星输出磁矩M(i)为最优卫星输出磁矩M-best,参照图图4,则,磁力矩器输出姿态捕获控制力矩对皮纳卫星进行姿态捕获控制处理,磁力矩器输出的姿态捕获控制力矩Tb为:
Tb=cross(M-best,B'b)
其中,B’b为卫星本体坐标系中地磁场矢量量测值。
本发明的有益效果是:高集成、模块化,可移植性:由于本模块完全封装,集成度高,为皮纳卫星提供最快速的姿态确定于控制能力,模块化、可移植,通过预留的打孔位置,可以满足很多皮纳卫星的需要;在皮纳卫星***中相对独立,仅通过供电即可完成对皮纳卫星的控制;接口简单,本产品能够实时将有关姿态的控制指令与量测数据通过RS232串口传输,方便用户使用,可市场化。
尽管以上结合附图对本发明的实施方案进行了描述,但是本发明并不局限于上述的具体实施方案和应用领域,上述的具体实施方案仅仅是示意性的、指导性的,而不是限制性的。本领域的普通技术人员在说明书的启示下,在不脱离本发明权利要求所保护的范围的情况下,还可以做出很多种的形式,这些均属于本发明保护之列。
Claims (2)
1.一种皮纳卫星的模块化姿态确定与控制装置,其特征在于:包括用于测定皮纳卫星姿态的姿态敏感单元,用于根据皮纳卫星姿态发出命令信号的姿控计算机,用于执行姿控计算机发出的命令的姿态控制执行单元;
所述姿态敏感单元包含用于测量皮纳卫星在惯性系下的角速度的陀螺仪,用于测量卫星本体坐标系下地磁场矢量的磁强计,用于测量卫星本体坐标系下太阳矢量的太阳敏感器,用于测量皮纳卫星的位置信息的卫星导航接收机;所述姿控计算机的四个输入端分别对应电连接所述陀螺仪的输出端、卫星导航接收机的输出端、磁强计的输出端、太阳敏感器的输出端;
所述姿态控制执行单元包含磁力矩器驱动单元,用于完成对皮纳卫星实时姿态控制的磁力矩器;所述姿控计算机的输出端电连接磁力矩器驱动单元的输入端,所述磁力矩器驱动单元的输出端电连接磁力矩器的控制端;
所述姿态敏感单元还包含温度传感器和模拟量电流传感器,所述温度传感器的输出端和模拟量电流传感器的输出端对应电连接所述姿控计算机的两个输入端;
所述姿态敏感单元、姿控计算机和姿态控制执行单元设置在外壳(1)内,所述外壳(1)的底部设置为开口状,所述外壳(1)的下底面设置有向外延伸的折弯面(2),所述折弯面(2)上设置有多个通孔,所述外壳(1)的侧壁上设置有供电接口(3);
所述外壳(1)内还设置有5V升压模块和3.3V稳压模块,所述供电接口与5V升压模块电连接,所述3.3V稳压模块与5V升压模块电连接;
所述供电接口为J30J-9ZKW-J接口,所述磁力矩器为三轴磁力矩器;所述姿控计算机采用单片微型计算机;
基于上述皮纳卫星的模块化姿态确定与控制装置的皮纳卫星姿态确定与控制方法为:包括以下步骤,
步骤一,定义卫星本体坐标系为Obxyz,参考坐标系为轨道坐标系Ooxyz,太阳矢量和地磁场矢量在Ooxyz中分别为S1、S2,太阳矢量和地磁场矢量在Obxyz中分别为
步骤二,皮纳卫星进入任务设定轨道后,在姿控计算机内预存有皮纳卫星的旋转角度的设置阈值Δw,陀螺仪实时测量皮纳卫星的旋转角度w,旋转角度w大于姿控计算机内的设置阈值Δw,启动姿态消旋控制对皮纳卫星进行消旋处理,直至旋转角度w小于等于姿控计算机内的设置阈值Δw后,进入姿态确定和姿态捕获控制的循环阶段;
所述消旋处理通过B-dot方法对皮纳卫星进行速率阻尼,
B-dot方法的计算式为:
式中,Mb为卫星输出磁矩,
K为消旋控制率,
Bb为地磁场矢量在卫星本体坐标系中的投影,通过磁强计测量得到,
为Bb的变化率,
则,磁力矩器输出消旋力矩对皮纳卫星进行消旋处理,磁力矩器输出的消旋力矩为:
所述姿态确定的过程为:
设s3=s1×s2,令s1、s2、s3三个矢量组成矩阵s的三列
S=(S1 S2 S3),
令矩阵s*为
则有:
Ao→b为轨道坐标系到本体坐标系的姿态旋转矩阵,
S、S*不共线,矩阵S可逆,Ao→b=S*S-1,
根据姿态动力学中姿态旋转矩阵Ao→b与欧拉角的转换关系反推欧拉角的值为:
式中,a21代表姿态旋转矩阵Ao→b当中2行1列的元素,a11代表姿态旋转矩阵Ao→b当中1行1列的元素,a31代表姿态旋转矩阵Ao→b当中3行1列的元素,a32代表姿态旋转矩阵Ao→b当中3行2列的元素,a33代表姿态旋转矩阵Ao→b当中3行3列的元素;从而完成姿态确定;
所述姿态捕获控制采用PD控制,
设期望磁控力矩Tc:
式中:Kp为比例参数,
Kd为微分参数,
θ为余仰角,θ=90°-δ,δ为卫星的赤纬;
为余仰角速度,
期望磁控力矩Tc的三轴分量为Tcx、Tcy、Tcz,将期望磁控力矩Tc的三轴分量Tcx、Tcy、Tcz大小按照顺序排列有如下六种排列方式:
Tcx>Tcy>Tcz、Tcx>Tcz>Tcy、Tcy>Tcx>Tcz、Tcy>Tcz>Tcx、Tcz>Tcx>Tcy、Tcz>Tcy>Tcx,
通过公式Tc=cross(M(i),Bb),cross表示叉乘,在不同排列方式下可以分别解得一组卫星输出磁矩M(i):
M(i)=[mxi、myi、mzi],
其中,i=1、2、3、4、5、6,
mxi为卫星输出磁矩M(i)在卫星本体坐标系中x轴上的分量,
myi为卫星输出磁矩M(i)在卫星本体坐标系中y轴上的分量,
mzi为卫星输出磁矩M(i)在卫星本体坐标系中z轴上的分量,
则:Tb(i)=M(i)×Bb,
式中,Tb(i)为每一种可能性下的实际能输出的力矩;
设dis(i)=||Tc|-|Tb(i)||,
使得dis(i)最小的卫星输出磁矩M(i)为最优卫星输出磁矩M,
则,磁力矩器输出姿态捕获控制力矩对皮纳卫星进行姿态捕获控制处理,磁力矩器输出的姿态捕获控制力矩Tb为:
Tb=cross(M,B'b),
其中,B′b为卫星本体坐标系中地磁场矢量量测值。
2.如权利要求1所述的皮纳卫星的模块化姿态确定与控制装置,其特征在于:在六种不同排列方式下,所述卫星输出磁矩M(i)的计算过程如下:
Bb的三轴分量为Bbx、Bby、Bbz,
Tcy>Tcx>Tcz时:
在|Bbx|>|Bbz|的状态下,设mz=Tcy/Bbx,然后判断是否|mz|>mmax,若|mz|>mmax为否,则mz=mz且mx=0,若|mz|>mmax为是,则设mz=sign(mz)mmax且mx=-(Tcx-mzBbx)/Bbz;再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
在|Bbx|≤|Bbz|的状态下,设mx=-Tcy/Bbz,然后判断是否|mx|>mmax,若|mx|>mmax为否,则mx=mx且mz=0,若|mx|>mmax为是,则设mx=sign(mx)mmax且mz=(Tcy-mxBbz)/Bbx,再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
然后计算my,再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
Tcz>Tcy>Tcx时:
在|Bby|>|Bbx|的状态下,设mx=Tcz/Bby,然后判断是否|mx|>mmax,若|mx|>mmax为否,则mx=mx且my=0,若|mx|>mmax为是,则设mx=sign(mx)mmax且my=-(Tcz-mxBby)/Bbx;再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
在|Bby|≤|Bbx|的状态下,设my=-Tcz/Bbx,然后判断是否|my|>mmax,若|my|>mmax为否,则my=my且mx=0,若|my|>mmax为是,则设my=sign(my)mmax且mx=(Tcz-myBbx)/Bby,再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
然后计算mz,再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
Tcy>Tcz>Tcx时:
在|Bbx|>|Bbz|的状态下,设mz=Tcy/Bbx,然后判断是否|mz|>mmax,若|mz|>mmax为否,则mz=mz且mx=0,若|mz|>mmax为是,则设mz=sign(mz)mmax且mx=-(Tcx-mzBbx)/Bbz;再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
在|Bbx|≤|Bbz|的状态下,设mx=-Tcy/Bbz,然后判断是否|mx|>mmax,若|mx|>mmax为否,则mx=mx且mz=0,若|mx|>mmax为是,则设mx=sign(mx)mmax且mz=(Tcy+mxBbz)/Bbx,再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
然后计算my,再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
Tcz>Tcx>Tcy时:
在|Bby|>|Bbx|的状态下,设mx=Tcz/Bby,然后判断是否|mx|>mmax,若|mx|>mmax为否,则mx=mx且my=0,若|mx|>mmax为是,则设mx=sign(mx)mmax且my=-(Tcz-mxBby)/Bbx;再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
在|Bby|≤|Bbx|的状态下,设my=-Tcz/Bbx,然后判断是否|my|>mmax,若|my|>mmax为否,则my=my且mx=0,若|my|>mmax为是,则设my=sign(my)mmax且mx=(Tcz-myBbx)/Bby,再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
然后计算mz,再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
Tcx>Tcy>Tcz时:
在|Bby|>|Bbz|的状态下,设my=Tcx/Bbz,然后判断是否|my|>mmax,若|my|>mmax为否,则my=my且mz=0,若|my|>mmax为是,则设my=sign(my)mmax且mz=-(Tcx-myBbz)/Bby;再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
在|Bby|≤|Bbz|的状态下,设mz=-Tcx/Bby,然后判断是否|mz|>mmax,若|mz|>mmax为否,则mz=mz且my=0,若|mz|>mmax为是,则设mz=sign(mz)mmax且my=(Tcx+mzBby)/Bbz,再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
然后计算mx,再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
Tcx>Tcz>Tcy时:
在|Bby|>|Bbz|的状态下,设my=Tcx/Bbz,然后判断是否|my|>mmax,若|my|>mmax为否,则my=my且mz=0,若|my|>mmax为是,则设my=sign(my)mmax且mz=-(Tcx-myBbz)/Bby;再判断是否|mz|>mmax,若|mz|>mmax为是,则mz=sign(mz)mmax,若|mz|>mmax为否,则mz=mz;
在|Bby|≤|Bbz|的状态下,设mz=-Tcx/Bby,然后判断是否|mz|>mmax,若|mz|>mmax为否,则mz=mz且my=0,若|mz|>mmax为是,则设mz=sign(mz)mmax且my=(Tcx+mzBby)/Bbz,再判断是否|my|>mmax,若|my|>mmax为是,则my=sign(my)mmax,若|my|>mmax为否,则my=my;
然后计算mx,再判断是否|mx|>mmax,若|mx|>mmax为是,则mx=sign(mx)mmax,若|mx|>mmax为否,则mx=mx;
其中,mmax=NIS,
式中,N表示线圈匝数,I表示线圈内通过的电流,S表示线圈围成的面积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610143290.4A CN105667838B (zh) | 2016-03-14 | 2016-03-14 | 一种皮纳卫星的模块化姿态确定与控制装置及其方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610143290.4A CN105667838B (zh) | 2016-03-14 | 2016-03-14 | 一种皮纳卫星的模块化姿态确定与控制装置及其方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105667838A CN105667838A (zh) | 2016-06-15 |
CN105667838B true CN105667838B (zh) | 2017-08-11 |
Family
ID=56307664
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610143290.4A Expired - Fee Related CN105667838B (zh) | 2016-03-14 | 2016-03-14 | 一种皮纳卫星的模块化姿态确定与控制装置及其方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105667838B (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106092106B (zh) * | 2016-06-20 | 2018-11-23 | 哈尔滨工业大学 | 星敏感器与磁传感器间欧拉角标定方法 |
CN106005483B (zh) * | 2016-07-06 | 2018-07-24 | 西北工业大学 | 一种模块化手机星的主动姿态控制方法 |
CN107054697B (zh) * | 2017-03-09 | 2019-03-12 | 西北工业大学 | 一种纳卫星磁力矩器空间温度补偿姿态控制方法 |
CN107444675B (zh) * | 2017-08-18 | 2019-10-01 | 上海微小卫星工程中心 | 一种航天器速率阻尼控制方法 |
CN107491060A (zh) * | 2017-08-20 | 2017-12-19 | 南京理工大学 | 一种立方星姿控分***调试方法 |
CN108072372A (zh) * | 2017-11-28 | 2018-05-25 | 西北工业大学 | 一种用于微小卫星的集成姿态测量*** |
CN108168757B (zh) * | 2017-12-14 | 2019-10-22 | 哈尔滨工业大学 | 空间翻滚目标电磁消旋微小消旋力测试平台及其测试方法 |
CN110803304B (zh) * | 2018-05-02 | 2021-08-10 | 上海微小卫星工程中心 | 一种卫星姿态控制*** |
CN109823571B (zh) * | 2019-01-23 | 2020-07-03 | 清华大学 | 一种遥感微纳卫星的多阶段姿态控制方法 |
CN110775302B (zh) * | 2019-09-24 | 2021-02-09 | 北京控制工程研究所 | 一种基于太阳帆板输出电流信息的应急对日方法 |
CN110589028B (zh) * | 2019-09-29 | 2021-07-06 | 上海航天控制技术研究所 | 一种卫星姿态机动异常时的自主模式切换方法 |
CN110667892B (zh) * | 2019-10-31 | 2021-03-19 | 北京国电高科科技有限公司 | 基于地磁测量的卫星消旋控制方法 |
CN110849391B (zh) * | 2019-11-01 | 2021-07-13 | 上海卫星工程研究所 | 矢量磁力仪与星敏感器非正交安装矩阵标定***及方法 |
CN111056045B (zh) * | 2019-12-26 | 2021-03-23 | 航天行云科技有限公司 | 一种三轴磁力矩器的控制方法以及装置 |
CN111532456B (zh) * | 2020-05-15 | 2021-11-12 | 清华大学 | 一种控制方法和装置 |
CN111874269B (zh) * | 2020-08-10 | 2022-02-01 | 吉林大学 | 一种磁控小卫星的低功耗对日捕获及定向姿态控制方法 |
CN112607064B (zh) * | 2020-12-25 | 2021-08-06 | 上海交通大学 | 避免进入自旋状态的微纳卫星磁阻尼控制方法 |
CN115140318B (zh) * | 2022-07-20 | 2023-05-05 | 南京理工大学 | 一种适用于微纳卫星大角速率消旋的磁阻尼控制方法 |
CN115108047B (zh) * | 2022-07-20 | 2023-04-18 | 南京理工大学 | 一种面向空间在轨服务的失稳航天器附着式消旋方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
BE797511A (fr) * | 1973-03-29 | 1973-07-16 | Org Europeene De Rech | Systeme de mesure d'attitude de satellite |
CN103838242A (zh) * | 2014-03-18 | 2014-06-04 | 西北工业大学 | 一种acs*** |
-
2016
- 2016-03-14 CN CN201610143290.4A patent/CN105667838B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105667838A (zh) | 2016-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105667838B (zh) | 一种皮纳卫星的模块化姿态确定与控制装置及其方法 | |
CN105116910B (zh) | 一种对地面点凝视成像的卫星姿态控制方法 | |
CN1323282C (zh) | 补偿地磁传感器姿态误差的设备和方法 | |
CN106094887B (zh) | 一种集成三轴转台功能及稳定平台功能的*** | |
CN106200614B (zh) | 一种采用控制力矩陀螺真实力矩的航天器姿控测试***及方法 | |
CN109823571A (zh) | 一种遥感微纳卫星的多阶段姿态控制方法 | |
CN104898642A (zh) | 一种用于航天器姿态控制算法的集成测试仿真*** | |
CN102354123A (zh) | 一种跨平台可扩展的卫星动态仿真测试*** | |
CN107054697B (zh) | 一种纳卫星磁力矩器空间温度补偿姿态控制方法 | |
CN105966639A (zh) | 一种卫星对日自旋稳定控制***及方法 | |
Hoang et al. | Pre-processing technique for compass-less madgwick in heading estimation for industry 4.0 | |
Hoang et al. | Measurement optimization for orientation tracking based on no motion no integration technique | |
CN102081360A (zh) | 惯性天文组合导航半物理实验*** | |
CN106599376B (zh) | 一种绳系卫星拉力方向估计方法 | |
CN111007865A (zh) | 以对日指向偏差为约束的卫星平稳对地定向方法 | |
Babcock | CubeSat attitude determination via Kalman filtering of magnetometer and solar cell data | |
CN106649947B (zh) | 基于李群谱算法的卫星姿态数值仿真方法 | |
Blachuta et al. | Attitude and heading reference system based on 3D complementary filter | |
CN107588786B (zh) | 一种用于星敏感器仿真测试的多用途恒星模拟器驱动方法 | |
CN110779554A (zh) | 一种机械臂、基于imu的初始位姿的校准***与方法 | |
Hong et al. | Application of EKF for missile attitude estimation based on “SINS/CNS” integrated guidance system | |
Qiu et al. | Two-wheeled self-balancing robot modeling and nonlinear control | |
CN111649738B (zh) | 微重力场下的加速度计初始姿态解算方法 | |
CN110568387B (zh) | 一种基于磁梯度张量的航天器磁矩测试方法 | |
Ziying et al. | Design and implementation of a MEMS-based attitude angle measuring system for moving objects |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170811 |