CN102288340B - 反作用飞轮输出力矩测量电路及其测量方法 - Google Patents
反作用飞轮输出力矩测量电路及其测量方法 Download PDFInfo
- Publication number
- CN102288340B CN102288340B CN 201110120115 CN201110120115A CN102288340B CN 102288340 B CN102288340 B CN 102288340B CN 201110120115 CN201110120115 CN 201110120115 CN 201110120115 A CN201110120115 A CN 201110120115A CN 102288340 B CN102288340 B CN 102288340B
- Authority
- CN
- China
- Prior art keywords
- speed
- microprocessor
- pulse
- reaction wheels
- reaction
- 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
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
反作用飞轮输出力矩测量电路及其测量方法。它涉及反作用飞轮输出力矩测量领域,它解决了动力学模型计算过程中反作用飞轮输出力矩测量的问题,它的微处理器连通信接口、计数器和整型模块,并采集飞轮转速方向信号和50ms脉冲信号,整型模块采集飞轮转速脉冲信号,计数器还连晶振,通信接口输出动力学计算模块的信号,方法如下:微处理器初始化,计数器初始化,50ms脉冲中断优于飞轮转速脉冲中断;50ms脉冲中断:记录T50,设置flag=1,退出中断;飞轮转速脉冲中断:一、count值加1,若flag=1,则进行二,否退出中断;二、读取Tflag,计算当前计算周期内的反作用飞轮转速,三、计算力矩,四、将力矩发送给动力学模型,退出中断。它应用卫星姿态控制系统动态测试。
Description
技术领域
本发明涉及反作用飞轮输出力矩测量领域。
背景技术
反作用飞轮是小卫星,尤其是高精度三轴稳定小卫星普遍采用的控制部件,通过不同的输入电压输入得到不同的输出力矩,它具有控制精度高、寿命长、连续控制等特点。采用了反作用飞轮的卫星在进行电测时,尤其在进行姿态控制系统闭环测试或者模飞测试时,需要精确测量反作用飞轮的输出力矩,用于闭环测试的动力学模型计算。国内常用的反作用飞轮一般输出转速和转速方向信号,转速采用单位时间内的脉冲数表示,为TTL电平的方波信号,转速方向采用TTL电平信号表示(比如高电平表示正转,低电平表示反转)。反作用飞轮的转速不太高,每分钟最大不超过几千转,每转采用特定个数(比如18个)脉冲表示。
在卫星姿态控制系统闭环测试时,需要实时、精确测量反作用飞轮输出力矩,实时是指在一次动力学计算步长(几十毫秒)内至少提供一次测量数据。如果采用脉冲计数法,假设飞轮转速信号频率小于2000Hz,动力学计算步长50ms,那么一个计算步长内的转速脉冲个数为小于100,由于计数为整数,计数误差为1,输出力矩的计算误差至少为1%,反作用飞轮转速越低误差越大,不能满足使用要求。如果采用测量脉冲周期的方法,一次周期测量的误差一般要十几微秒(由处理器的中断时间决定),对于2000Hz的信号来说,输出力矩的计算最大误差接近2%(10μs/500μs)。且一个动力学计算步长内有多个脉冲,选择哪个脉冲计算该时间内的转速最合适,能够提高反作用飞轮输出力矩的精度,从而满足卫星测试使用要求,成为现在的一个研究方向。
发明内容
本发明的目的是解决卫星姿态控制系统动态测试时动力学模型计算过程中反作用飞轮输出力矩测量的问题,提供一种反作用飞轮输出力矩测量电路及其测量方法。
本发明反作用飞轮输出力矩测量电路包括微处理器、通信接口、计数器、晶振和整形模块;微处理器的第一输入输出端与通信接口的第一输出输入端连接,通信接口的第二输入输出端为动力学计算模块的信号输入输出端,微处理器的第二输入输出端与计数器的输出输入端连接,微处理器的第三输入端与整形模块的输出端连接,整形模块的输入端为飞轮转速脉冲信号输入端,微处理器的第四输入端为微处理器的飞轮转速方向信号输入端,微处理器的第五输入端为微处理器的50ms脉冲信号输入端,计数器的输入端与晶振的输出端连接。
反作用飞轮输出力矩测量方法如下:
微处理器初始化,设置微处理器的飞轮转速脉冲信号输入端的中断口的中断触发方式为上升沿或者下降沿触发方式,即单个边沿触发,设置微处理器的50ms脉冲信号输入端的中断口的中断触发方式为上升沿触发;微处理器定义飞轮转速脉冲中断计数变量count并清零;设置微处理器内部定时器,定时精度1ms,启动定时;
计数器初始化,设置计数器3的计数增加1代表1μs,启动计数器计数;初始化50ms中断更新标志flag=0;
设置50ms脉冲中断处理程序的优先级高于飞轮转速脉冲中断处理程序的优先级;
50ms脉冲中断处理程序为:
记录微处理器内部定时器当前计数值T50,设置50ms中断更新标志flag=1,退出50ms脉冲中断处理程序;
飞轮转速脉冲中断处理程序为:
步骤一:飞轮转速脉冲中断计数变量count值加1,读取飞轮转速方向信号Direct,若50ms中断更新标志flag=1,则进行步骤二,否则退出飞轮转速脉冲中断处理程序;
步骤二:读取微处理器内部定时器当前计数值Tflag,计算当前计算周期内的反作用飞轮转速:
若Tflag-T50>25ms,则当前计算周期内的反作用飞轮转速ω=0,否则读取并记录计数器的数值T3,按下述公式计算当前计算周期内的反作用飞轮转速:
t=T3k-T3k-1,其中下标k代表当前计算周期,k-1代表上一计算周期,
其中,ω表示当前计算周期内的飞轮转速,一周期50ms,Direct表示转速方向,由飞轮转速方向信号决定,如果飞轮转速方向信号为高电平,则Direct=1,否则Direct=-1;
T3k=T3k-1;设置50ms中断更新标志flag=0,
步骤三:计算反作用飞轮输出力矩T:
其中,T代表飞轮输出力矩,J代表飞轮转动惯量,为常值,ωk表示当前计算周期的飞轮转速,ωk-1表示上一计算周期飞轮转速;ωk=ωk-1;
步骤四:将计算得到的飞轮输出力矩T通过通信接口发送给动力学模型,退出飞轮转速脉冲中断处理程序。
本发明通过采集一个动力学计算步长内的反作用飞轮的转速脉冲数,以及该时间内首尾脉冲边沿之间的时间间隔,计算反作用飞轮输出力矩,提高了计算精度,满足卫星测试使用要求。本发明所测量的飞轮输出力矩,由于同时记录相邻两次50ms脉冲上升沿后的第一个飞轮转速脉冲中断之间的时间间隔以及中断次数,脉冲计数是准确的,误差主要为时间测量误差,由于计数器精度为1μs,因此时间误差主要来自微处理器中断处理引起的误差,大约10μs,飞轮输出力矩的计算误差约10μs/50ms=0.0002,精度提高了近百倍。该电路设计简单、精度高,满足卫星系统电测试的使用要求。
附图说明
图1是本发明的反作用飞轮输出力矩测量电路的结构示意图;图2是整形模块5的结构示意图;图3是倍频处理电路52的信号逻辑图,其中输入信号为光电隔离整形芯片51的输出端的输出信号,延时信号为延时芯片521的输出端的输出信号,输出信号为异或模块522的输出端的输出信号;图4是飞轮转速脉冲采样点示意图。
具体实施方式
具体实施方式一:结合图1说明本实施方式,本实施方式包括微处理器1、通信接口2、计数器3、晶振4和整形模块5;微处理器1的第一输入输出端与通信接口2的第一输出输入端连接,微处理器1的第一输入输出端为通信输入输出端,通信接口2的第二输入输出端为动力学计算模块的信号输入输出端,微处理器1的第二输入输出端与计数器3的输出输入端连接,微处理器1的第三输入端与整形模块5的输出端连接,微处理器1的第一输入端为微处理器1的飞轮转速脉冲信号输入端,整形模块5的输入端为飞轮转速脉冲信号输入端,微处理器1的第四输入端为微处理器1的飞轮转速方向信号输入端,微处理器1的第五输入端为微处理器1的50ms脉冲信号输入端,计数器3的输入端与晶振4的输出端连接。
微处理器1是该测量电路的核心,用于发送飞轮输出力矩的计算数据,接收动力学计算的指令;用于控制计数器的读写;用于测量飞轮转速脉冲信号;用于测量飞轮转速方向信号;用于接收50ms脉冲信号;用于飞轮输出力矩的计算;
通信接口2用于微处理器1与动力学计算模型之间的通信;
计数器3通过计数晶振4脉冲数用于记录时间;
整形模块5用于对飞轮转速脉冲信号进行隔离整形后进行倍频处理。
微处理器1采用C8051F040;通信接口2采用RS232接口。
具体实施方式二:结合图2说明本实施方式,本实施方式与具体实施方式一不同点在于整形模块5包括光电隔离整形芯片51和倍频处理电路52,
光电隔离整形芯片51的输入端为飞轮转速脉冲信号输入端,光电隔离整形芯片51的输出端与倍频处理电路52的输入端连接,
倍频处理电路52包括延时芯片521和异或模块522,倍频处理电路52的输入端同时为延时芯片521的输入端和异或模块522的一个输入端,延时芯片521的输出端与异或模块522的另一个输入端连接,异或模块522的输出端为倍频处理电路52的输出端,并且为整形模块5的输出端。
首先通过光电隔离整形芯片51对整形模块5的输入的飞轮转速脉冲信号进行整形与隔离,然后通过倍频处理电路52对输入的飞轮转速脉冲信号进行倍频。
倍频处理电路52将输入信号与输入信号延时后的信号进行异或运算,其运算结果如图3所示,在输入信号产生跳变开始到信号延时时间内输出信号为高电平,其他时间输出信号为低电平,从而将输入信号的频率提高了一倍。被测信号倍频处理后,与微处理器中断输入口连接。
其它组成和连接方式与具体实施方式一相同。
具体实施方式三:本实施方式的飞轮输出力矩的测量方法:
微处理器1初始化,设置微处理器1的飞轮转速脉冲信号输入端的中断口的中断触发方式为上升沿或者下降沿触发方式,即单个边沿触发,设置微处理器1的50ms脉冲信号输入端的中断口的中断触发方式为上升沿触发;微处理器1定义飞轮转速脉冲中断计数变量count并清零;设置微处理器1内部定时器,定时精度1ms,启动定时;
计数器3初始化,设置计数器3的计数增加1代表1μs,启动计数器3计数;初始化50ms中断更新标志flag=0;
设置50ms脉冲中断处理程序的优先级高于飞轮转速脉冲中断处理程序的优先级;
50ms脉冲中断处理程序为:
记录微处理器1内部定时器当前计数值T50,设置50ms中断更新标志flag=1,退出50ms脉冲中断处理程序;
飞轮转速脉冲中断处理程序为:
步骤一:飞轮转速脉冲中断计数变量count值加1,读取飞轮转速方向信号Direct,若50ms中断更新标志flag=1,则进行步骤二,否则退出飞轮转速脉冲中断处理程序;
步骤二:读取微处理器1内部定时器当前计数值Tflag,计算当前计算周期内的反作用飞轮转速:
若Tflag-T50>25ms,则当前计算周期内的反作用飞轮转速ω=0,否则读取并记录计数器3的数值T3,按下述公式计算当前计算周期内的反作用飞轮转速:
t=T3k-T3k-1,其中下标k代表当前计算周期,k-1代表上一计算周期,
其中,ω表示当前计算周期内的飞轮转速,一周期50ms,Direct表示转速方向,由飞轮转速方向信号决定,如果飞轮转速方向信号为高电平,则Direct=1,否则Direct=-1;
T3k=T3k-1;设置50ms中断更新标志flag=0,
步骤三:计算反作用飞轮输出力矩T:
其中,T代表飞轮输出力矩,J代表飞轮转动惯量,为常值,ωk表示当前计算周期的飞轮转速,ωk-1表示上一计算周期飞轮转速;ωk=ωk-1;
步骤四:将计算得到的飞轮输出力矩T通过通信接口2发送给动力学模型,退出飞轮转速脉冲中断处理程序。
如图4所示,飞轮正转,A、C点为相邻两次50ms脉冲的上升沿,计算反作用飞轮输出力矩T,需要记录A点后的第一次飞轮转速脉冲中断B点时刻tb,C点后的第一次飞轮转速脉冲中断D点时刻td,B点与D点之间中断次数为5,如果tb-ta<25ms,且td-tc<25ms,则反作用飞轮输出力矩T为:
本发明内容不仅限于上述各实施方式的内容,其中一个或几个具体实施方式的组合同样也可以实现发明的目的。
Claims (3)
1.反作用飞轮输出力矩测量电路,其特征在于它包括微处理器(1)、通信接口(2)、计数器(3)、晶振(4)和整形模块(5);微处理器(1)的第一输入输出端与通信接口(2)的第一输出输入端连接,通信接口(2)的第二输入输出端为动力学计算模块的信号输入输出端,微处理器(1)的第二输入输出端与计数器(3)的输出输入端连接,微处理器(1)的第三输入端与整形模块(5)的输出端连接,整形模块(5)的输入端为飞轮转速脉冲信号输入端,微处理器(1)的第四输入端为微处理器(1)的飞轮转速方向信号输入端,微处理器(1)的第五输入端为微处理器(1)的50ms脉冲信号输入端,计数器(3)的输入端与晶振(4)的输出端连接。
2.根据权利要求1所述的反作用飞轮输出力矩测量电路,其特征在于整形模块(5)包括光电隔离整形芯片(51)和倍频处理电路(52),光电隔离整形芯片(51)的输入端为飞轮转速脉冲信号输入端,光电隔离整形芯片(51)的输出端与倍频处理电路(52)的输入端连接,倍频处理电路(52)包括延时芯片(521)和异或模块(522),倍频处理电路(52)的输入端同时为延时芯片(521)的输入端和异或模块(522)的一个输入端,延时芯片(521)的输出端与异或模块(522)的另一个输入端连接,异或模块(522)的输出端为倍频处理电路(52)的输出端,并且为整形模块(5)的输出端。
3.反作用飞轮输出力矩测量方法,其特征在于步骤如下:
微处理器(1)初始化,设置微处理器(1)的飞轮转速脉冲信号输入端的中断口的中断触发方式为在每个上升沿和每个下降沿触发,即单个边沿触发,设置微处理器(1)的50ms脉冲信号输入端的中断口的中断触发方式为上升沿触发;微处理器(1)定义飞轮转速脉冲中断计数变量count并清零;设置微处理器(1)内部定时器,定时精度1ms,启动定时;
计数器(3)初始化,设置计数器(3)的计数增加1代表1μs,启动计数器(3)计数;初始化50ms中断更新标志flag=0;
设置50ms脉冲中断处理程序的优先级高于飞轮转速脉冲中断处理程序的优先级;
50ms脉冲中断处理程序为:
记录微处理器(1)内部定时器当前计数值T50,设置50ms中断更新标志flag=1,退出50ms脉冲中断处理程序;
飞轮转速脉冲中断处理程序为:
步骤一:飞轮转速脉冲中断计数变量count值加1,读取飞轮转速方向信号Direct,若50ms中断更新标志flag=1,则进行步骤二,否则退出飞轮转速脉冲中断处理程序;
步骤二:读取微处理器(1)内部定时器当前计数值Tflag,计算当前计算周期内的反作用飞轮转速:
若Tflag-T50>25ms,则当前计算周期内的反作用飞轮转速ω=0,否则读取并记录计数器(3)的数值T3,按下述公式计算当前计算周期内的反作用飞轮转速:
t=T3k-T3k-1,其中下标k代表当前计算周期,k-1代表上一计算周期,
其中,ω表示当前计算周期内的飞轮转速,一周期50ms,Direct表示转速方向,由飞轮转速方向信号决定,如果飞轮转速方向信号为高电平,则Direct=1,否则Direct=-1;
T3k=T3k-1;设置50ms中断更新标志flag=0;步骤三:计算反作用飞轮输出力矩T:
其中,T代表飞轮输出力矩,J代表飞轮转动惯量,为常值,ωk表示当前计算周期的飞轮转速,ωk-1表示上一计算周期飞轮转速;ωk=ωk-1;
步骤四:将计算得到的飞轮输出力矩T通过通信接口(2)发送给动力学模型,退出飞轮转速脉冲中断处理程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110120115 CN102288340B (zh) | 2011-05-10 | 2011-05-10 | 反作用飞轮输出力矩测量电路及其测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110120115 CN102288340B (zh) | 2011-05-10 | 2011-05-10 | 反作用飞轮输出力矩测量电路及其测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102288340A CN102288340A (zh) | 2011-12-21 |
CN102288340B true CN102288340B (zh) | 2013-01-16 |
Family
ID=45334946
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110120115 Expired - Fee Related CN102288340B (zh) | 2011-05-10 | 2011-05-10 | 反作用飞轮输出力矩测量电路及其测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102288340B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102735399B (zh) * | 2012-07-12 | 2014-07-09 | 杭州电子科技大学 | 直流电机飞轮惯量检测电路 |
CN102866259B (zh) * | 2012-08-21 | 2015-06-10 | 浙江宇视科技有限公司 | 一种风扇速度检测装置 |
CN103345275B (zh) * | 2013-06-07 | 2015-08-19 | 哈尔滨工业大学 | 基于角动量裕度的单轴批量反作用飞轮力矩优化分配方法 |
CN105738690B (zh) * | 2016-02-05 | 2018-11-02 | 上海微小卫星工程中心 | 一种复杂机动卫星的反作用飞轮动态功耗评估方法 |
CN106681175A (zh) * | 2017-02-24 | 2017-05-17 | 上海航天控制技术研究所 | 减小反作用飞轮接入半物理系统时延迟的离散化方法 |
CN114279467B (zh) * | 2021-12-24 | 2023-06-06 | 深圳航天科技创新研究院 | 一种反作用飞轮性能参数智能评估系统及其评估方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2474995Y (zh) * | 2000-12-18 | 2002-01-30 | 中国科学院光电技术研究所 | 一种光电传感器 |
CN100335877C (zh) * | 2006-04-30 | 2007-09-05 | 北京航空航天大学 | 控制力矩陀螺用力矩测试气浮转台 |
JP2009190506A (ja) * | 2008-02-13 | 2009-08-27 | Mitsubishi Electric Corp | 人工衛星用姿勢制御装置および人工衛星の姿勢制御方法 |
KR20100075073A (ko) * | 2008-12-24 | 2010-07-02 | 한국항공우주연구원 | 보정기구, 보정방법 및 이를 이용한 3축 힘토크센서 |
-
2011
- 2011-05-10 CN CN 201110120115 patent/CN102288340B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102288340A (zh) | 2011-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102288340B (zh) | 反作用飞轮输出力矩测量电路及其测量方法 | |
CN102495912B (zh) | 一种具有同步校正功能的多通道高速数据采集系统 | |
CN104443427B (zh) | 飞行器颤振预测系统及方法 | |
CN101685135B (zh) | 用于测试电路板的测试装置和方法 | |
CN103248356B (zh) | 一种基于采用锁相环脉冲插值技术的计数器及实现方法 | |
CN104461972B (zh) | 一种数据信号采样的方法和设备 | |
CN105300508B (zh) | 一种超声波信号传播时间的双阈值检测电路及检测方法 | |
CN102680726B (zh) | 一种用于电机转速测量的高精度自适应装置 | |
CN103954903B (zh) | 一种可实时解算的多模式输出电路测试系统 | |
CN109283354A (zh) | 一种基于增量式光电编码器的变m/t测速方法 | |
CN104569472A (zh) | 一种基于光电脉冲编码器的速度获取方法 | |
CN107065688A (zh) | 一种高品质卫星平台时间同步误差的控制方法 | |
CN104406642A (zh) | 一种时差法超声波流量计精确测量方法 | |
CN100514068C (zh) | 滑动时间窗口法测量旋转轴转速的方法 | |
CN100520419C (zh) | 改进的频率确定 | |
CN104111481B (zh) | 同步时钟相位差测量系统和方法 | |
CN103743923B (zh) | 一种高频车速采集系统及采集方法 | |
CN106771326B (zh) | 一种基于增量式光电编码器测量速度的方法 | |
CN103528739A (zh) | 冲量测试系统 | |
CN105467396A (zh) | 一种脉冲激光测距控制系统及测距运算方法 | |
CN202393904U (zh) | 一种用于电能计量装置的误差检验装置 | |
CN101881784B (zh) | 一种基于感应同步器或旋转变压器的位置及速度测量装置 | |
CN106508018B (zh) | 基于飞机运动方程的飞行数据“野值”检测与剔除方法 | |
CN201387481Y (zh) | 用于数字式电能表的误差检测装置 | |
RU113012U1 (ru) | Ультразвуковой измеритель скоростей потока |
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: 20130116 Termination date: 20130510 |