CN101929872B - 应用于单轴气浮台的星敏感器的简易模拟装置及模拟方法 - Google Patents
应用于单轴气浮台的星敏感器的简易模拟装置及模拟方法 Download PDFInfo
- Publication number
- CN101929872B CN101929872B CN2010102832273A CN201010283227A CN101929872B CN 101929872 B CN101929872 B CN 101929872B CN 2010102832273 A CN2010102832273 A CN 2010102832273A CN 201010283227 A CN201010283227 A CN 201010283227A CN 101929872 B CN101929872 B CN 101929872B
- Authority
- CN
- China
- Prior art keywords
- angular velocity
- air bearing
- star sensor
- bearing table
- chip microcomputer
- 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
Images
Landscapes
- Navigation (AREA)
Abstract
应用于单轴气浮台的星敏感器的简易模拟装置及模拟方法,涉及一种应用于单轴气浮台的星敏感器的简易模拟装置及模拟方法。它解决了现有地面航天器半物理仿真中的星敏感器由于无法接入仿真系统中导致无法对控制算法进行性能考核的问题。其装置:光电码盘的转动部分套在单轴气浮台的转轴上并固定,其信号通过无线通信模块发送到单轴气浮台的台面上的单片机中。其方法:采用光电码盘测量单轴气浮台的角速度,并将角速度数据通过无线通信模块发送至单轴气浮台台面上的单片机中。本发明适用于单轴气浮台星敏感器的功能模拟。
Description
技术领域
本发明涉及一种应用于单轴气浮台的星敏感器的简易模拟装置及模拟方法。
背景技术
航天器控制系统半物理仿真是研制航天器过程中特有的一种仿真方法,它利用气浮转台作为运动模拟器并结合部分实物搭建的半物理仿真环境,进行控制系统方案的仿真验证。单轴气浮转台通过压缩空气在气浮轴承与轴承座之间形成气膜使台面浮于空中,实现气浮转台台面与台体之间近似无摩擦的相对转动,从而模拟航天器在外层空间所受干扰力矩很小的力学环境,进而在地面上利用气浮转台台面的转动模拟刚体航天器的姿态运动。
在基于单轴气浮转台的航天器半物理仿真中,为了考验关键部件对航天器控制算法控制性能的影响,常需要将实物接入到实际控制回路中。这些实物既包括飞轮、陀螺等必要器件,也包括决定航天器高精度控制不可或缺的星敏感器。目前,受星敏感器工作原理的限制,地面半物理仿真时多采用星空模拟器配合方法实现对其性能的考核。然而驱动星空模拟器需要用户提供航天器精确的轨道和姿态信息,这显然有悖于采用星敏感器进行高精度姿态测量的初衷。因此在地面航天器半物理仿真中,星敏感器无法接入闭环仿真系统中对算法进行性能考核。
发明内容
本发明是为了解决现有地面航天器半物理仿真中的星敏感器由于无法接入仿真系统中导致无法对控制算法进行性能考核的问题,从而提出一种应用于单轴气浮台的星敏感器的简易模拟装置及模拟方法。
应用于单轴气浮台的星敏感器的简易模拟装置,它包括光电码盘、转速数字显示仪和无线通信模块,所述无线通信模块包括无线信号发射单元和无线信号接收单元;所述光电码盘的转动部分套在单轴气浮台的转轴上并固定,光电码盘的接收部分固定在单轴气浮台的底座上,所述光电码盘的信号输出端与转速数字显示仪的信号输入端连接;所述转速数字显示仪的信号输出端与无线信号发射单元的信号输入端连接;无线信号接收单元固定在单轴气浮台的台面上。
基于上述装置的应用于单轴气浮台的星敏感器的简易模拟方法,它由以下步骤实现:
步骤一、采用光电码盘测量单轴气浮台的角速度,获得角速度数据,并将所述角速度数据传送至转速数字显示仪;
步骤二、转速数字显示仪显示步骤一获得的角速度数据,并将所述角速度数据通过无线信号发射单元发射;
步骤三、无线信号接收单元接收步骤二发射的角速度数据,并将所述角速度数据传送至单片机;
步骤四、仿真计算机计算航天器另外两个轴的角速度数据和轨道信息并传送至单片机;
步骤五、单片机根据步骤三中接收的角速度数据和步骤四接收的航天器另外两个轴的角速度数据和轨道信息进行计算,获得星敏感器的模拟输出数据。
步骤五所述的根据步骤三中接收的角速度数据和步骤四接收的航天器另外两个轴的角速度数据和轨道信息进行计算,具体步骤是:将接收到的角速度数据利用四阶龙格库塔法进行积分,获得到航天器体坐标系相对于轨道坐标系的姿态四元数Qob;
步骤三中接收到的气浮台的角速度为一维矢量wob,与步骤四接收的航天器另外两个轴的角速度组成航天器三轴的角速度三维矢量w;
航天器上积分时间常数dTime;
航天器初始姿态角速度为w0;
航天器初始姿态四元数为Q0;
四元数计算公式为:
根据姿态四元数运动学方程表示的函数:
计算如下四个系数:
K1=f(Q0,w0)
根据上述获得的四个系数,计算得到姿态四元数Qob:
仿真计算机向单片机中实时输入的航天器运行轨道信息:瞬时轨道六根数,计算轨道坐标系相对惯性坐标系的姿态四元数Qio;
所述瞬时轨道六根数为:升交点赤经Ω,近地点幅角ω,偏心率e,平近点角M,轨道倾角i;
通过公式:
计算:
其中:
获得的轨道坐标系相对惯性坐标系的姿态四元数Qio即实际星敏感器测量的航天器姿态,即星敏感器模拟装置的输出数据。
无线信号发射单元和无线信号接收单元均采用蓝牙模块实现。
转速数字显示仪的信号输出端与无线信号发射单元的信号输入端之间的信号传输、无线信号接收单元的信号输出端和单片机的信号输入端之间的信号传输均通过RS485总线实现。
有益效果:本发明采用光电码盘测量单轴气浮台高精度转动角速度数据,并将所述角速度数据通过无线信号发射单元发射至单轴气浮台的台面上的无线信号接收单元,然后采用单片机计算出对应的航天器某个轴的转动角度,通过将单片机接入单轴气浮台的台面上的仿真系统,从而实现对控制算法进行性能的考核,进而实现星敏感器高精度测量功能和性能的模拟。本发明具有体积小、延迟时间小、模拟数据真实且精度高等特点。
附图说明
图1是本发明的结构示意图。
具体实施方式
具体实施方式一、结合图1说明本具体实施方式,应用于单轴气浮台的星敏感器的简易模拟装置,它包括光电码盘1、转速数字显示仪2和无线通信模块3,所述无线通信模块3包括无线信号发射单元31和无线信号接收单元32;所述光电码盘1的转动部分套在单轴气浮台的转轴上并固定,光电码盘1的接收部分固定在单轴气浮台的底座上,所述光电码盘1的信号输出端与转速数字显示仪2的信号输入端连接;所述转速数字显示仪2的信号输出端与无线信号发射单元31的信号输入端连接;无线信号接收单元32固定在单轴气浮台的台面上。光电码盘1的接收部分位于其转动部分的工作面中,保证能够准确接收到转动部分的光信号。
本具体实施方式中,转速数字显示仪2和无线信号发射单元31均设置在单轴气浮台的底座上或地面上。
本实施方式采用光电码盘1测量单轴气浮台高精度转动角速度数据,并将所述角速度数据通过无线信号发射单元发射至单轴气浮台的台面上的无线信号接收单元,然后采用单片机4计算出航天器的姿态四元数并发送给仿真计算机5,从而实现对控制算法进行性能的考核,进而实现星敏感器高精度测量功能和性能的模拟。
具体实施方式二、本具体实施方式与具体实施方式一所述的应用于单轴气浮台的星敏感器的简易模拟装置的区别在于,它还包括单片机4,单片机4固定在单轴气浮台的台面上,所述无线信号接收单元32的信号输出端与单片机4的信号输入端连接。
具体实施方式三、本具体实施方式与具体实施方式二所述的应用于单轴气浮台的星敏感器的简易模拟装置的区别在于,它还包括仿真计算机5,所述仿真计算机5固定在单轴气浮台的台面上,所述仿真计算机5的信号输入或输出端与单片机4的信号输出或输入端连接。
本具体实施方式中,单片机4通过一个串行接口与无线信号接收单元32的串口连接,用于接收数据,通过另一串行接口作为模拟器信号输出端,用于与仿真计算机5连接。
本实施方式中,单片机4和无线通信模块3的工作电源电压均为+5v,由单轴气浮转台台面上的仿真系统电源提供。
本实施方式具有体积小、延迟时间小、原理简单、研制成本低、模拟数据真实且精度高的优点。
具体实施方式四、本具体实施方式与具体实施方式三所述的应用于单轴气浮台的星敏感器的简易模拟装置的区别在于,它还包括RS485总线,所述光电码盘1的信号输出端与转速数字显示仪2的信号输入端之间、转速数字显示仪2的信号输出端与无线信号发射单元31的信号输入端之间、无线信号接收单元32的信号输出端和单片机4的信号输入端之间、仿真计算机5的信号输入或输出端与单片机4的信号输出或输入端之间均通过RS485总线连接。
具体实施方式五、本具体实施方式与具体实施方式四所述的应用于单轴气浮台的星敏感器的简易模拟装置的区别在于,无线通信模块3为蓝牙无线通信模块。
具体实施方式六、本具体实施方式与具体实施方式五所述的应用于单轴气浮台的星敏感器的简易模拟装置的区别在于,光电码盘1的型号为EAS425绝对值光电码盘(上海精浦机电有限公司),无线通信模块3的型号为BTS5504C2P(重庆金瓯科技发展有限责任公司);单片机4的型号为C8051F040;转速数字显示仪2的型号为XJP-48E(上海上自仪转速表仪表电机有限公司),仿真计算机5的型号为MSM855(瑞士Digital-Logic公司)。
本实施方式中,光电码盘1的工作电源电压为+10v直流电,转速数字显示仪2的工作电源电压为220v交流电,由地面供电电源提供。
本实施方式中各器件的型号选取均不是唯一的。
具体实施方式七、基于具体实施方式一所述的应用于单轴气浮台的星敏感器的简易模拟装置的应用于单轴气浮台的星敏感器的简易模拟方法,它由以下步骤实现:
步骤一、采用光电码盘1测量单轴气浮台的角速度,获得角速度数据,并将所述角速度数据传送至转速数字显示仪2;
步骤二、转速数字显示仪2显示步骤一获得的角速度数据,并将所述角速度数据通过无线信号发射单元31发射;
步骤三、无线信号接收单元32接收步骤二发射的角速度数据,并将所述角速度数据传送至单片机4;
步骤四、仿真计算机5计算航天器另外两个轴的角速度数据和轨道信息并传送至单片机4;
步骤五、单片机4根据步骤三中接收的角速度数据和步骤四接收的航天器另外两个轴的角速度数据和轨道信息进行计算,获得星敏感器的模拟输出数据。
步骤五所述的根据步骤三中接收的角速度数据和步骤四接收的航天器另外两个轴的角速度数据和轨道信息进行计算,具体步骤是:将接收到的角速度数据利用四阶龙格库塔法进行积分,获得到航天器体坐标系相对于轨道坐标系的姿态四元数Qob;
步骤三中接收到的气浮台的角速度为一维矢量wob,与步骤四接收的航天器另外两个轴的角速度组成航天器三轴的角速度三维矢量w;
航天器上积分时间常数dTime;
航天器初始姿态角速度为w0;
航天器初始姿态四元数为Q0;
四元数计算公式为:
根据姿态四元数运动学方程表示的函数:
计算如下四个系数:
K1=f(Q0,w0)
根据上述获得的四个系数,计算得到姿态四元数Qob:
仿真计算机(5)向单片机(4)中实时输入的航天器运行轨道信息:瞬时轨道六根数,计算轨道坐标系相对惯性坐标系的姿态四元数Qio;
所述瞬时轨道六根数为:升交点赤经Ω,近地点幅角ω,偏心率e,平近点角M,轨道倾角i;
通过公式:
计算:
其中:
获得的轨道坐标系相对惯性坐标系的姿态四元数Qio即实际星敏感器测量的航天器姿态,即星敏感器模拟装置的输出数据。
无线信号发射单元31和无线信号接收单元32均采用蓝牙模块实现。
转速数字显示仪2的信号输出端与无线信号发射单元31的信号输入端之间的信号传输、无线信号接收单元32的信号输出端和单片机4的信号输入端之间的信号传输均通过RS485总线实现。
获得的轨道坐标系相对惯性坐标系的姿态四元数Qio即实际星敏感器测量的航天器姿态,即星敏感器模拟装置的输出数据。
本实施方式中,蓝牙模块的传输速率设为19200bps,数据发送频率为10Hz。光电码盘1测量精度为20角秒,因此模拟器的测量精度最高可达20角秒,与目前国际上典型的星敏感器的指标一致。
本实施方式中,单片机4将计算得到的航天器姿态四元数Qio输出至仿真计算机5,为仿真计算机5提供姿态信息,进而实现航天器姿态控制系统仿真。
Claims (7)
1.应用于单轴气浮台的星敏感器的简易模拟装置,其特征是:它包括光电码盘(1)、转速数字显示仪(2)和无线通信模块(3),所述无线通信模块(3)包括无线信号发射单元(31)和无线信号接收单元(32);所述光电码盘(1)的转动部分套在单轴气浮台的转轴上并固定,光电码盘(1)的接收部分固定在单轴气浮台的底座上,所述光电码盘(1)的信号输出端与转速数字显示仪(2)的信号输入端连接;所述转速数字显示仪(2)的信号输出端与无线信号发射单元(31)的信号输入端连接;无线信号接收单元(32)固定在单轴气浮台的台面上;
光电码盘(1)的接收部分位于其转动部分的工作面中,保证能够准确接收到转动部分的光信号;
它还包括单片机(4),单片机(4)固定在单轴气浮台的台面上,所述无线信号接收单元(32)的信号输出端与单片机(4)的信号输入端连接;
它还包括仿真计算机(5),所述仿真计算机(5)固定在单轴气浮台的台面上,所述仿真计算机(5)的信号输入或输出端与单片机(4)的信号输出或输入端连接;
光电码盘(1)用于测量单轴气浮台高精度转动角速度数据,并将所述角速度数据通过无线信号发射单元发射至单轴气浮台的台面上的无线信号接收单元,单片机(4)用于计算出航天器的姿态四元数并发送给仿真计算机(5),从而实现对控制算法进行性能的考核,进而实现星敏感器高精度测量功能和性能的模拟。
2.根据权利要求1所述的应用于单轴气浮台的星敏感器的简易模拟装置,其特征在于它还包括RS485总线,所述光电码盘(1)的信号输出端与转速数字显示仪(2)的信号输入端之间、转速数字显示仪(2)的信号输出端与无线信号发射单元(31)的信号输入端之间、无线信号接收单元(32)的信号输出端和单片机(4)的信号输入端之间、仿真计算机(5)的信号输入或输出端与单片机(4)的信号输出或输入端之间均通过RS485总线连接。
3.根据权利要求2所述的应用于单轴气浮台的星敏感器的简易模拟装置,其特征在于无线通信模块(3)为蓝牙无线通信模块。
4.根据权利要求3所述的应用于单轴气浮台的星敏感器的简易模拟装置,其特征在于光电码盘(1)的型号为EAS425绝对值光电码盘;无线通信模块(3)的型号为BTS5504C2P;单片机(4)的型号为C8051F040;转速数字显示仪(2)的型号为XJP-48E;仿真计算机(5)的型号为MSM855。
5.基于权利要求1所述的应用于单轴气浮台的星敏感器的简易模拟装置的应用于单轴气浮台的星敏感器的简易模拟方法,其特征是:它由以下步骤实现:
步骤一、采用光电码盘(1)测量单轴气浮台的角速度,获得角速度数据,并将所述角速度数据传送至转速数字显示仪(2);
步骤二、转速数字显示仪(2)显示步骤一获得的角速度数据,并将所述角速度数据通过无线信号发射单元(31)发射;
步骤三、无线信号接收单元(32)接收步骤二发射的角速度数据,并将所述角速度数据传送至单片机(4);
步骤四、仿真计算机(5)计算航天器另外两个轴的角速度数据和轨道信息并传送至单片机(4);
步骤五、单片机(4)根据步骤三中接收的角速度数据和步骤四接收的航天器另外两个轴的角速度数据和轨道信息进行计算,获得星敏感器的模拟输出数据;具体过程是:将接收到的角速度数据利用四阶龙格库塔法进行积分,获得到航天器体坐标系相对于轨道坐标系的姿态四元数Qob;
步骤三中接收到的气浮台的角速度为一维矢量wob,与步骤四接收的航天器另外两个轴的角速度组成航天器三轴的角速度三维矢量w;
航天器上积分时间常数dTime;
航天器初始姿态角速度为w0;
航天器初始姿态四元数为Q0;
四元数计算公式为:
根据姿态四元数运动学方程表示的函数:
计算如下四个系数:
K1=f(Q0,w0)
根据上述获得的四个系数,计算得到姿态四元数Qob:
仿真计算机(5)向单片机(4)中实时输入的航天器运行轨道信息:瞬时轨道六根数,计算轨道坐标系相对惯性坐标系的姿态四元数Qio;
所述瞬时轨道六根数为:升交点赤经Ω,近地点幅角ω,偏心率e,平近点角M,轨道倾角i;
通过公式:
计算:
其中:
获得的轨道坐标系相对惯性坐标系的姿态四元数Qio即实际星敏感器测量的航天器姿态,即星敏感器模拟装置的输出数据。
6.根据权利要求5所述的应用于单轴气浮台的星敏感器的简易模拟方法,其特征在于无线信号发射单元(31)和无线信号接收单元(32)均采用蓝牙模块实现。
7.根据权利要求5所述的应用于单轴气浮台的星敏感器的简易模拟方法,其特征在于转速数字显示仪(2)的信号输出端与无线信号发射单元(31)的信号输入端之间的信号传输、无线信号接收单元(32)的信号输出端和单片机(4)的信号输入端之间的信号传输均通过RS485总线实现。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102832273A CN101929872B (zh) | 2010-09-16 | 2010-09-16 | 应用于单轴气浮台的星敏感器的简易模拟装置及模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102832273A CN101929872B (zh) | 2010-09-16 | 2010-09-16 | 应用于单轴气浮台的星敏感器的简易模拟装置及模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101929872A CN101929872A (zh) | 2010-12-29 |
CN101929872B true CN101929872B (zh) | 2012-07-18 |
Family
ID=43369180
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102832273A Expired - Fee Related CN101929872B (zh) | 2010-09-16 | 2010-09-16 | 应用于单轴气浮台的星敏感器的简易模拟装置及模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101929872B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102520719A (zh) * | 2011-12-06 | 2012-06-27 | 北京邮电大学 | 具有五自由度的微重力气浮目标卫星模拟器系统 |
CN103206966B (zh) * | 2013-04-12 | 2015-04-29 | 哈尔滨工业大学 | 一种单轴气浮台精密测量误差校正方法 |
CN105259793A (zh) * | 2015-11-19 | 2016-01-20 | 中国人民解放军国防科学技术大学 | 基于aodv协议无线智能网络的卫星姿态控制仿真系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101503116B (zh) * | 2009-02-17 | 2010-12-08 | 哈尔滨工业大学 | 一种分布式航天器地面仿真系统及其实现方法 |
CN101487698A (zh) * | 2009-02-25 | 2009-07-22 | 哈尔滨工业大学 | 单轴气浮台高精度姿态角测量方法及其装置 |
CN101706512B (zh) * | 2009-11-25 | 2011-06-15 | 哈尔滨工业大学 | 基于星敏感器的姿态测量信息和飞轮的角动量测量信息的航天器伪速率的估计方法 |
CN101813542B (zh) * | 2010-04-23 | 2012-02-01 | 哈尔滨工业大学 | 单轴气浮转台台面外加载荷后自动调整质心的装置及方法 |
CN101833285A (zh) * | 2010-06-12 | 2010-09-15 | 哈尔滨工业大学 | 编队卫星姿态协同控制仿真装置及其仿真验证方法 |
-
2010
- 2010-09-16 CN CN2010102832273A patent/CN101929872B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101929872A (zh) | 2010-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102879014B (zh) | 深空探测接近过程的光学成像自主导航半物理仿真试验系统 | |
CN106200614B (zh) | 一种采用控制力矩陀螺真实力矩的航天器姿控测试系统及方法 | |
CN201865663U (zh) | 一种惯性测量单元以及动力调谐陀螺连续测斜系统 | |
CN103575299B (zh) | 利用外观测信息的双轴旋转惯导系统对准及误差修正方法 | |
CN104898642B (zh) | 一种用于航天器姿态控制算法的集成测试仿真系统 | |
CN101424534B (zh) | 惯性/重力组合导航半实物模拟装置 | |
CN104296908B (zh) | 三自由度气浮台干扰力矩组成测量装置 | |
CN101706284B (zh) | 提高船用光纤陀螺捷联惯导系统定位精度的方法 | |
CN102538819B (zh) | 基于双圆锥红外和星敏感器的自主导航半物理仿真试验系统 | |
CN102003968B (zh) | 光纤陀螺捷联惯性导航系统的单轴转台标定方法 | |
CN101876244A (zh) | 一种惯性测量单元以及动力调谐陀螺连续测斜仪 | |
CN102278108B (zh) | 小口径定向陀螺测斜仪连续测量模式标定方法 | |
CN101995824A (zh) | 星箭一体化航天器姿态控制半物理仿真系统 | |
CN103344226A (zh) | 一种基于mems旋转技术的寻北系统及方法 | |
CN103868648A (zh) | 三轴气浮仿真实验平台的质心测量方法 | |
CN102889076A (zh) | 陀螺测斜仪标定方法 | |
CN101706287A (zh) | 一种基于数字高通滤波的旋转捷联系统现场标定方法 | |
CN102323759B (zh) | 翼伞自主归航半实物仿真系统 | |
CN104848818B (zh) | Stewart平台姿态测量装置及测量方法 | |
CN102636183A (zh) | 基于光纤监测和双轴转台离心机的挠性陀螺二次过载项测试方法 | |
CN105157723A (zh) | 一种基于光纤陀螺的捷联惯导系统的标定方法 | |
CN101929872B (zh) | 应用于单轴气浮台的星敏感器的简易模拟装置及模拟方法 | |
CN100491204C (zh) | 一种利用定轨数据标定加速度计的方法 | |
CN102636184B (zh) | 无角运动环境下基于离心机的挠性陀螺比力敏感项标定方法 | |
CN103674068A (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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120718 Termination date: 20120916 |