CN105867167B - 跨平台的自旋稳定卫星飞行仿真系统 - Google Patents
跨平台的自旋稳定卫星飞行仿真系统 Download PDFInfo
- Publication number
- CN105867167B CN105867167B CN201610228224.7A CN201610228224A CN105867167B CN 105867167 B CN105867167 B CN 105867167B CN 201610228224 A CN201610228224 A CN 201610228224A CN 105867167 B CN105867167 B CN 105867167B
- Authority
- CN
- China
- Prior art keywords
- module
- control
- data
- processing
- telemetering
- 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
- 238000004088 simulation Methods 0.000 title claims abstract description 46
- 238000012545 processing Methods 0.000 claims abstract description 102
- 238000005094 computer simulation Methods 0.000 claims abstract description 37
- 238000005259 measurement Methods 0.000 claims abstract description 19
- 238000012544 monitoring process Methods 0.000 claims abstract description 19
- 230000005540 biological transmission Effects 0.000 claims abstract description 18
- 238000009825 accumulation Methods 0.000 claims abstract description 11
- 239000000446 fuel Substances 0.000 claims abstract description 11
- 230000033001 locomotion Effects 0.000 claims abstract description 10
- 238000009987 spinning Methods 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 29
- 230000008569 process Effects 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 20
- 230000001133 acceleration Effects 0.000 claims description 10
- 238000000926 separation method Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 9
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims description 8
- 230000006870 function Effects 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 4
- 238000004080 punching Methods 0.000 claims description 4
- 238000002485 combustion reaction Methods 0.000 claims description 2
- 238000005314 correlation function Methods 0.000 claims description 2
- 238000002347 injection Methods 0.000 claims description 2
- 239000007924 injection Substances 0.000 claims description 2
- 230000008450 motivation Effects 0.000 claims 1
- 230000003993 interaction Effects 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 10
- OAKJQQAXSVQMHS-UHFFFAOYSA-N Hydrazine Chemical compound NN OAKJQQAXSVQMHS-UHFFFAOYSA-N 0.000 description 4
- 230000001186 cumulative effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- JEIPFZHSYJVQDO-UHFFFAOYSA-N ferric oxide Chemical compound O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及一种跨平台的自旋稳定卫星飞行仿真系统，主要包括飞行仿真控制系统和动力学仿真系统。其中飞行仿真控制系统包括运行控制、数据发送、数据接收、遥控处理、遥测处理、监视显示和实时记盘模块；动力学仿真系统主要包括数据输入模块、四级发动机点火模拟、燃料量/贮箱压力及发动机推力计算、加速度计章动测量值计算、自主章动控制计算、章控/姿控的选择、控制力及力矩计算、姿态运动动力学、控制总冲累加计算和轨道量测量计算。本系统具有较强的通用性，可以与自旋稳定卫模平台、载人航天卫模平台、探索卫模平台进行数据交互和模拟演练，完成卫星自主章控、转速控制、时延测量、姿态控制、应急章控、应急起旋等状态的模拟仿真。
Description
技术领域
本发明涉及一种跨平台的自旋稳定卫星模拟验证系统，特别是一种跨平台的自旋稳定卫星飞行仿真系统。卫模仿真软件运行于测控系统仿真计算机上，充分模拟卫星飞行时的各种技术状态，完成卫星自主章控、转速控制、时延测量、姿态控制、应急章控、应急起旋等状态的模拟仿真。
背景技术
自旋稳定卫星测控任务准备过程中，需要使用卫星测控模拟器与船载测控系统协同工作，完成星地大回路演练，检验、验证测控系统的软硬件技术状态的正确性和协调性。
卫星测控模拟器正常情况下完成自星箭分离至同步定点各飞行段的卫星运行过程遥测、外测和遥控仿真，主要包括：控制系统敏感器输出的姿态脉冲的动态仿真；卫星各飞行段的姿态(含章动)和转速变化过程仿真，包括星箭分离、转移轨道、准同步轨道、定点后，及地面参与控制(姿态机动、章动控制、轨道控制)时卫星姿态的实时变化过程；完成以上卫星飞行段和变轨段的轨道动力学仿真，并实时输出卫星轨道信息；完成星箭分离至四级发动机分离前卫星章动加速度测量和自主章控回路的仿真；仿真卫星推进系统遥测参数(肼瓶压力与温度、催床温度等)随轨道、姿态控制的实时变化。异常情况下完成卫星不起旋或起旋不正常的仿真；卫星自主章控失效的仿真。
目前测控系统尚未有一套成熟的仿真软件系统，具备与各套卫模平台(自旋稳定卫模平台、探索卫模平台、载人航天卫模平台)进行数据交换、遥测遥控数据仿真和动力学仿真的能力，不同型号试验任务必须单独研发相应的模拟器仿真软件，造成了仿真软件配置项多，技术状态多，管理困难等问题。因此，研制一种跨平台的自旋稳定卫星飞行仿真系统具有重要意义。
发明内容
本发明所要解决的技术问题是针对上述现有技术提供一种跨平台的自旋稳定卫星飞行仿真系统，采用高级抽象通用化的数据接口设计，解决了卫星飞行仿真系统与各套卫模平台之间进行数据交换、模拟仿真的问题，可以作为一种通用性较强的自旋稳定卫星飞行仿真系统。
本发明解决上述问题所采用的技术方案为：一种跨平台的自旋稳定卫星飞行仿真系统，该系统主要包括飞行仿真控制系统和动力学仿真系统两部分，其中飞行仿真控制系统由仿真控制模块、数据发送模块、数据接收模块、遥控处理模块、遥测处理模块、监视显示模块和记盘处理模块组成；动力学仿真系统由数据输入模块、四级发动机点火模拟模块、燃料量/压力/推力计算模块、日/月星历计算模块、章控阀值确定模块、加速度计章动测量值计算模块、自主章动控制计算模块、平均推力系数确定模块、章控/姿控选择模块、控制力及力矩计算模块、星体质量特性计算模块、姿态运动动力学计算模块、姿态脉冲时刻计算模块、控制总冲累加计算模块、轨道量测量计算模块、轨道运动动力学计算模块、遥测帧中断服务模块、遥控指令接收模块和时统中断服务模块组成，其中，
仿真控制模块用于接收并处理人机界面命令标识，负责初始化、启动相关功能模块、创建各功能事件、生成系统参数描述表和加载动力学仿真系统，并在系统退出时结束系统的运行；
数据发送模块用于发送数据接收模块输出的信息，并将数据信息送监视显示模块进行显示处理和实时记盘模块进行记盘处理；
数据接收模块用于接收管道数据，进行正确性检查和判断，并放入相应的缓冲区，处理后通知遥控处理模块、遥测处理模块、监视显示模块和实时记盘模块，或转发给动力学仿真系统；
遥控处理模块用于实时接收和处理遥控数据帧、执行脉宽和控制参数，将相关遥控数据帧、执行脉宽和控制参数发送至动力学仿真系统，由动力学仿真系统模拟卫星的姿态变化过程；并将所有遥控指令、执行脉宽进行编码遥测反演计算，通过遥测信号模拟卫星各系统随遥控指令的变化过程；
遥测处理模块用于接收动力学仿真系统仿真的卫星编码遥测、模拟遥测信息以及遥控处理功能模块仿真的编码遥测信息，组装卫星编码遥测和模拟遥测信息帧，形成网发信息送数据发送模块；
监视显示模块用于监视显示本机软、硬件工作状态；
实时记盘模块用于完成任务执行过程中的收发数据记盘工作；
数据输入模块用于在每个仿真周期的开始时，向系统注入发动机工作脉冲信息，四级发动机点火模拟模块用于根据收到的四级脉冲，计算四级点火起始时刻和终止时刻；燃料量/压力/推力计算模块用于计算发动机系统的剩余燃料量、贮箱压力以及相应发动机推力；日/月星历计算模块用于根据经典的日、月星历的计算公式计算太阳、月亮位置；章控阀值确定模块用于根据遥控指令章控阀值，确定阀值；加速度计章动测量值计算模块用于模拟加速度计测量值经主动章控回路的带通滤波器后的输出，同时组装形成编码遥测信息；自主章动控制计算模块用于进行粗、精控判断，形成控制脉冲前后沿，提供遥测输出及状态字设置；发动机平均推力系数确定模块用于计算发动机平均推力系数；章控/姿控选择模块用于根据遥控指令进行章控和姿控的选择设置；控制力及力矩计算模块用于根据发动机的推力，计算力矩和力矩在惯性系投影，并计算本次工作周期的力冲量的累加值；星体质量特性计算模块用于计算四级点火过程、二次分离后卫星的质量及惯量，计算二次分离后的惯量比和章动频率；姿态运动动力学计算模块用于计算卫星角速度及姿态角；姿态脉冲时刻计算模块用于根据卫星实时姿态、太阳星历和卫星轨道位置，计算姿态脉冲出现的时刻，并按格式组装成模拟遥测数据；控制总冲累加计算模块用于本工作周期内的控制总冲的累加计算，提供给燃料量/压力/推力计算模块进行控制总冲累加值计算；轨道量测量计算模块用于测量船的位置及速度以及轨道量测量计算；轨道运动动力学计算模块用于计算地心距、地球半张角以及卫星轨道位置的赤经赤纬值；遥测帧中断服务模块用于向主程序输送遥测参数，并负责每秒显示刷新；遥控指令接收模块用于接收遥控指令，并存储、显示；时统中断服务模块用于接收时统中断板发送的定时中断信号，并向相应模块发出中断消息。
优选地，所述数据接收模块负责接收自旋稳定卫模平台、载人航天卫模平台、探索卫模平台的接口数据，转换为内部自定义格式的遥控数据，送遥控处理模块进行处理。
优选地，所述数据接收模块根据网收数据源地址，按信息分类在以下情况分别进行相应处理：在计算机演练模式下，对接收到的遥测遥控处理机数据，首先进行三判二处理，成功后，根据信息组装相应的返回信息帧和回答信息帧，并把信息交给遥控处理模块处理，同时将处理信息转发动力学仿真系统；在星地大回路演练模式下，收到卫星测控模拟平台帧中断时刻信息时，立即转发给动力学仿真系统；收到通用卫星测控模拟平台的应答机遥测数据信息时，交给遥测处理模块处理；其他数据信息立即交给遥控处理模块处理，数据接收模块利用管道技术，接收动力学仿真系统的数据，按以下情况处理：在计算机演练模式下，将收到的全帧编码遥测数据，实时发送遥测通知事件，交给遥测处理模块处理；将收到的全帧模拟遥测数据、返回的同控参数、姿章联控/应急章控、非同步控制单脉冲、非同步控制脉冲串等信息转发数据发送模块处理；在星地大回路演练模式下，将收到的编码遥测数据包，实时发送遥测通知事件，交给遥测处理功能模块处理；将收到的模拟遥测数据包，将遥测时间增加一工作周期后，转发数据发送模块处理。
与现有技术相比，本发明的优点在于：
1、本卫星飞行仿真系统具备与自旋稳定卫模平台、载人航天卫模平台、探索卫模平台数据交互的能力，满足了一套卫星飞行仿真系统适应三套卫模平台的仿真需求。
2、本卫星飞行仿真系统满足了大回路模式和计算机模式两种仿真模式下的卫星飞行仿真需求。
3、本卫星飞行仿真系统满足了自旋稳定卫星自主章控、转速控制、时延测量、姿态控制、应急章控、应急起旋等状态的模拟仿真。
附图说明
图1为飞行仿真控制系统与动力学仿真系统的数据流图。
图2为动力学仿真系统流程图。
图3为运行控制模块流程图。
图4为数据发送模块流程图。
图5为数据接收模块流程图。
图6为遥控处理模块流程图。
图7为遥测处理模块流程图。
图8为监视显示模块流程图。
图9为实时记盘模块流程图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
本发明涉及一种跨平台的自旋稳定卫星飞行仿真系统，该系统主要包括飞行仿真控制系统和动力学仿真系统两部分。其中飞行仿真控制系统由仿真控制模块11、数据发送模块12、数据接收模块13、遥控处理模块14、遥测处理模块15、监视显示模块16和记盘处理模块17组成；动力学仿真系统2由数据输入模块、四级发动机点火模拟模块202、燃料量/压力/推力计算模块203、日/月星历计算模块204、章控阀值确定模块205、加速度计章动测量值计算模块206、自主章动控制计算模块207、平均推力系数确定模块208、章控/姿控选择模块209、控制力及力矩计算模块210、星体质量特性计算模块211、姿态运动动力学计算模块212、姿态脉冲时刻计算模块213、控制总冲累加计算模块214、轨道量测量计算模块215、轨道运动动力学计算模块216、遥测帧中断服务模块217、遥控指令接收模块218和时统中断服务模块219组成。
卫星模拟器通用测控平台通过网络接口将星上6个发动机工作执行脉冲及四级点火脉冲及各种遥控指令，传送给飞行仿真系统。嵌在飞行仿真系统内的动力学仿真和控制软件接收到这些信息后，模拟卫星的姿态(包括章动)及轨道动力学，将接收的发动机工作参数，变换成相应的控制力及力矩，作用于姿态及轨道动力学方程，改变卫星的姿态(包括章动)、转速及轨道。并根据太阳、地球、卫星的实时位置及卫星姿态计算出姿态敏感器6个姿态脉冲时刻，根据卫星章动运动计算出加速度计测量信息以及自主章控回路有关参数，与控制有关的肼瓶压力等参数，传送给飞行仿真系统，由它转发给测控平台。测控平台将这些参数经相应变化后送到地面测控站。
如图3所示，仿真控制模块11接收处理人机界面交互命令，获取初始化配置文件完成系统的初始化工作，将资源句柄、同步事件送内存共享缓冲区；将监视显示信息送监视显示模块16显示处理；将控制信息数据发送模块12和数据接收模块13处理；将数据实时记盘控制信息送记盘处理模块17处理。
如图4所示，数据发送模块12从管道读取网发数据，对数据进行核检、判断网发管道标志进行网发处理，并对错误信息进行出错处理；查网发通讯信息描述表，按照各自的信宿标识，获取网发目的组地址及网发帧数据长度，实施网络发送。处理发送遥测处理模块15产生的编码编遥和模拟遥测数据模遥；处理发送动力学仿真系统2返回的同控参数、姿章联控/应急章控、非同步控制单脉冲、非同步控制脉冲串等信息以及动力学仿真软件仿真的轨道参数；处理发送遥控处理模块14产生的小环结果；并将处理发送的相关数据送监视显示模块16和实时记盘模块17进行显示处理和记盘处理。
如图5所示，数据接收模块13利用SOCKET接收网络数据，根据网收数据源地址，按信息分类在以下情况分别进行相应处理：在计算机——计算机演练模式下，对接收到的遥测遥控处理机数据，首先进行3判2处理，成功后，根据信息组装相应的返回信息帧和回答信息帧，并把信息交给遥控处理模块14处理，同时将同控器控制命令、同控参数、姿章联控/应急章控、非同步控制单脉冲、非同步控制脉冲串等信息转发动力学仿真系统2。在星——地大回路演练模式下，收到卫星测控模拟平台帧中断时刻(256ms包中断)信息时，立即转发给动力学仿真系统2；收到通用卫星测控模拟平台的应答机遥测数据信息时，交给遥测处理模块15处理；其他数据信息，立即交给遥控处理模块14处理。
数据接收模块13利用管道技术，接收动力学仿真系统2的数据，按以下情况处理：在计算机——计算机演练模式下，将收到的全帧编码遥测数据，实时发送遥测通知事件，交给遥测处理模块15处理；将收到的全帧模拟遥测数据、返回的同控参数、姿章联控/应急章控、非同步控制单脉冲、非同步控制脉冲串等信息转发数据发送模块12处理；在星——地大回路演练模式下，将收到的编码遥测数据包，实时发送遥测通知事件，交给遥测处理模块15处理；将收到的模拟遥测数据包，将遥测时间增加一工作周期后，转发数据发送模块12处理。
如图6所示，遥控处理模块14负责实时接收和处理遥控数据帧、执行脉宽和控制参数，将相关遥控数据帧、执行脉宽和控制参数发送至动力学仿真系统2，由动力学仿真系统2模拟卫星的姿态变化过程；并将所有遥控指令、执行脉宽进行编码遥测反演计算，通过遥测信号模拟卫星各系统随遥控指令的变化过程。可分为下列两种情况：
在星——地大回路演练模式中，遥控功能模块14实时接收卫星测控模拟平台发送的遥控数据帧、执行脉冲等信息；收到遥控数据帧后，将其转换为指令代号。将与动力学仿真系统2有关的遥控指令代号、执行脉冲实时转发给动力学仿真系统2，并将遥控指令代号和执行脉冲分类加工处理，按要求反演成编码遥测相应参数，发送给遥测处理模块15处理。
在计算机—计算机仿真模式下，遥控处理模块14实时接收遥测遥控处理机发送出的遥控指令代号、执行脉冲；当收到遥控代号时，实时组装小环比对结果送遥测遥控处理机；将与动力学仿真系统2有关的指令代号、执行脉冲等信息实时转发给动力学仿真系统2，并将遥控指令代号以及执行脉冲分类加工处理，按要求反演成编码遥测相应参数，发送给遥测处理模块15处理。
如图7所示，遥测处理模块15主要负责接收动力学仿真系统2仿真的卫星编码遥测、模拟遥测信息以及遥控处理功能模块仿真的编码遥测信息。在星——地大回路演练模式下，按卫星测控模拟平台的要求组装卫星编码遥测、模拟遥测信息帧，形成网发信息送数据发送模块12；在计算机——计算机演练模式下，按卫星遥测地面传输格式组装卫星编码遥测、模拟遥测信息，形成网发信息送数据发送能模块12。
如图8所示，监视显示模块16负责监视本机软、硬件工作状态，并实时显示各模块发送的数据帧。完成数据发送模块12的数据发送信息显示处理；完成数据接收模块13的数据接收信息显示处理；完成遥控处理模块14的遥控类信息显示处理；完成遥测处理模块15的遥测类信息显示处理。
如图9所示，记盘处理模块17完成仿真过程中的收发数据记盘处理，负责处理数据发送模块12和数据接收模块13产生的记盘数据。
数据输入模块：在每个仿真周期的开始，向飞行动力学仿真与控制软件注入发动机工作脉冲信息。根据收到的发动机脉冲的起始时刻tb、收到的发动机脉冲的终止时刻te，工作周期起始时刻t0(可由时统约定给出)，发动机脉冲类别特征数INW＝(1，0，-1)和发动机脉冲数N。计算发动机类别特征数IF(N)，N＝1……7。INW在一般情况下(长脉冲序列和星上自主章控粗控)置1，其它情况下根据需要置0或-1。可通过人机界面进行注入。发动机类别特征数IF(N)如下表所示。
若IF(i)＝1，表明发动机i工作；若IF(i)＝0，表明发动机i不工作。实际上传来的是表示发动机类别的标志码，通过计算机变成便于运算的特征数IF(N)。
四级发动机点火模拟模块：在收到四级脉冲后IF(7)置1，te被传给仿真软件，计算出四级点火起、终时刻t7b和t7e，为了提高精度，将四级发动机推力F(7)的计算放在控制力及力矩计算模块中。
燃料量、压力及推力计算模块：当星上自主章控回路工作时，无论是粗控还是精控，将采用分别累计A、B两发动机系统的发动机工作累积冲量SA、SB的方法。假设发动机的比冲为常值，则可由累积冲量求出燃料消耗量△MA、△MB。从而计算出剩余燃料量MA、MB。再利用气态方程求出贮箱压力PA、PB。最后利用贮箱压力与推力关系求出发动机推力，按长脉冲序工作状态(N>100)且脉宽为150ms，INW＝1计算。
日、月星历计算模块：根据经典的日，月星历的计算公式计算太阳、月亮位置。根据aE表示地球赤道半径。计算相对于历元时平春分点的太阳直角坐标ξS、ηS、ζS，春分点平移角θγ，轨道坐标系的太阳直角坐标XS、YS、ZS，太阳方向赤经、赤纬角αS、δS，月亮直角坐标ξS、ηS、ζS、XS、YS、ZS。算法或流程：先计算太阳轨道半长轴、偏心率、倾角、升交点赤经、近地点辐角以及平近点角，然后通过解开普勒方程给出输出数据。
章控阀值确定模块：根据遥控指令章控阀值A、B、C、D，确定阀值。当没有新遥控指令时，开、关阀值保持不变。输入遥控指令所确定的章控阀值，粗、精控门限值a1A、a3A；a1B、a3B；a1C、a3C；a1D、a3D，输出关阀值a2。算法或流程：关阀值a2＝a1-滞环宽。
加速度计章动测量值计算模块：模拟加速度计测量值经主动章控回路的带通滤波器后的输出，同时，形成编码遥测信息传送出去。INA、INB表示章控回路工作状态标志，由控制章控回路A、B工作的遥控指令确定，ψa(I)、γa(I)、ha(I)表示加速度计安装量；计算出加速度测量值A(I)，遥测加速度测量值Am(I)。算法或流程：仅当A路(或B路)章控回路工作时，加速度计章动测量值才产生，并经遥测传送地面；章控回路不工作时，则不计算也不下传。在编程中用数组A(I)来区别A，B路。I＝1时为A路，I＝2时为B路。考虑到经带通滤波后的放大及相移以及加速度计测量的饱和，加速度测量值可按如下计算：
AMX约为5.25米/秒2。引入INA、INB描述A、B路工作状态(工作时为1，不工作为0)，由遥控指令K165～K170的状态码确定。计算出的加速度章动测量值将通过编码遥测传到地面，每帧传送12个数值，间隔为44ms，记ta(I，K)为遥测值采样时间在一帧内的相对时间，其中I＝1，2代表A，B路，K＝1，……12表12个采样点。遥测加速度测量值按下式计算：
Am(I,K)＝Ca[A(I)-(A(I)-Ab(I))·(tz-ta(I,K)/△t)
其中Ab(I)为前一点加速度值，Ca为遥测因子。
加速度计输出峰值按下式计算：
加峰值＝惯量比×(2－惯量比)×0.929×ωZ 2×tan(章动角)；
自主章动控制计算模块进行粗、精控判断，形成控制脉冲前后沿，提供遥测输出及状态字设置。Ab(I)、A(I)表示前次及本次加速度计测量值。N1、N2、N3表示相位钟参数。输出正过零时刻TC0、负过零时刻TC和双脉冲起终了时刻TN1、TN2、TN3、TN4。算法或流程：该模块首先确定过零检测及双脉冲的起始和终了时刻。然后进行章动的粗、精控或不控的判断，这里引入了控制特征数，值为1表示控，为0表不控。随后，由粗、精控判断结果确定本步是否形成控制，当本步有控制时将控制特征数置为1，否则为零。
发动机平均推力系数确定模块计算发动机平均推力系数。输入发动机工作脉冲起终了时刻tb、te，发动机脉冲推力时延TPW，计算步长△t，本步起始时刻t。输出平均推力系数KF。
算法或流程：发动机平均推力系数是指发动机实际工作时间与计算步长的比，用平均推力在全步长积分，这样控制总量是精确的，保证了控制精度。
章控、姿控选择模块根据章控或姿控遥控指令进行章控和姿控的选择设置。IN表示卫星处于章控或姿控状态的标志。
算法或流程图：引入表示卫星是处于章控还是姿控的特征数IN，由遥控指令章控转姿控(K171)和姿控转章控(K172)决定。当IN为1时表示章控状态，当IN为0时表示姿控状态。
控制力及力矩计算模块：根据四级发动机推力(当其工作时)，以及6个发动机的推力，确定力矩，给出力在惯性系投影，并计算本次工作周期的力冲量的累加值。根据参数XF(N)、YF(N)、ZF(N)(N＝1～6)表示6个发动机的安装点坐标，iF(N)、jF(N)、kF(N)表示喷管推力单位方向在本体系上的投影，DLX(N)、DLY(N)、DLZ(N)表示单位推力下的力矩。计算输出LX、LY、LZ表示力矩在本体系上的投影量，FX、FY、FZ表示力在本体系上的投影量，SFXI、SFYI、SFZI表示一次工作周期内的力冲量累加。
算法或流程：已知了各发动机推力以及发动机在星体上的安装点及角度方向(给定常值)，则可以确定力矩，作为这步姿态动力学积分的力函数项。还要把每次计算步长内的力冲量进行累加，直至一个工作周期计算完成，获得的力的冲量累加输入到轨道动力学方程。
星体质量特性计算模块：计算四级点火过程、二次分离后卫星的质量及惯量，计算二次分离后的惯量比和章动频率。根据输入参数：△JX、△JY、△JZ表示四级点火前后卫星转动惯量的变化量，△MS表示四级点火前后卫星质量的变化量。计算输出：JX、JY、JZ表示四级点火过程、二次分离后卫星的转动惯量，MS表示四级点火过程、二次分离后卫星的质量，μ表示二次分离后卫星的惯量比，Ω表示章动频率。
姿态运动动力学模块：采用龙格——库塔法解姿态运动动力学方程。根据输入参数：τ表示章动发散时间常数，ρ1表示过渡轨道段发散时间常数，ρ2表示准同步轨道段阻尼时间常数，CX、CY、CZ表示变质量动力学阻尼系数，仅当四级点火时考虑。计算输出卫星角速度及姿态角ωX、ωY、ωZ、φ、θ、
姿态脉冲时刻计算模块：根据卫星实时姿态(欧拉角)，太阳星历，卫星轨道位置。计算姿态脉冲出现的时刻(包括标志码)，按照包含全部6个脉冲(太1、太2、南中、南出、北中、北出)的统一格式传送出去。输入参数：φ、θ、表示卫星实时欧拉角，ξS、ηS、ζS、XS、YS、ZS表示太阳星历位置，卫星轨道位置。计算输出姿态脉冲出现时刻。
红外脉冲时刻计算流程：南红外地平仪探头方向与自旋轴夹角95°，北红外地平仪探头方向与自旋轴夹角85°。卫星自旋时，红外探头扫到地球时产生一个扫地脉冲，经线路处理后得到扫入和扫出脉冲，再经取中线路处理得到地中脉冲，地中及地出脉冲经遥测传到地面。则南、北红外地平仪扫过地球时，将产生北中、北出、南中、南出四个红外脉冲，仿真软件将给出这些脉冲时刻，传送给飞行仿真主程序。
地中脉冲产生：由上一自旋周期出脉冲与入脉冲相减取半产生半弦宽并存储，到本自旋周期入脉冲出现后加上半弦宽形成地中脉冲。
控制总冲累加计算模块：计算本工作周期内的控制总冲的累加量。根据输入参数控制冲量AF(N)，计算输出控制总冲累加值S1A、S1B。算法或流程：
S1A＝S1A+AF(N)·△t，A系统控制总冲累加。
S1B＝S1B+AF(N)·△t，B系统控制总冲累加。
轨道量测量计算模块：计算测量船的位置、速度以及轨道量测量。根据输入参数测控站纬度φ(i)、测控站经度λ(i)、测控站高度h(i)和格林威治恒星时S。计算输出
测量船位置RX(i)、RY(i)、RZ(i)和测控船速度以及轨道量测量ρ、A、E。
轨道运动动力学模块：采用COWELL方法直角坐标形式的动力学方程，其优点是计算简单。动力学积分采用龙格—库塔方法。在积分给出卫星位置后，计算地心距、地球半张角以及卫星轨道位置的赤经赤纬值。
遥测帧中断服务模块217
遥测帧中断服务模块217向主程序输送遥测参数，并负责秒显示刷新。算法或流程：在一次工作周期计算完成之后，计算机将全部待输出的数据存贮单元内的数据，等到512ms帧中断信号时送出。每两次中断刷新一次屏幕显示。
遥控指令接收模块218
遥控指令接收模块218接收遥控指令，并存储、显示。算法或流程：每个控制周期都要进行遥控指令的接收和处理。先从飞行仿真主程序接收指令(消息)，各指令位的意义。然后，读到的指令数据放在相应的内存中。
时统中断服务219
接收时统中断板来的定时中断信号，向相应模块发出中断消息。
除上述实施例外，本发明还包括有其他实施方式，凡采用等同变换或者等效替换方式形成的技术方案，均应落入本发明权利要求的保护范围之内。
Claims (2)
1.一种跨平台的自旋稳定卫星飞行仿真系统，其特征在于：该系统包括飞行仿真控制系统(1)和动力学仿真系统(2)两部分，其中飞行仿真控制系统由仿真控制模块(11)、数据发送模块(12)、数据接收模块(13)、遥控处理模块(14)、遥测处理模块(15)、监视显示模块(16)和记盘处理模块(17)组成；动力学仿真系统(2)由数据输入模块、四级发动机点火模拟模块(202)、燃料量/压力/推力计算模块(203)、日/月星历计算模块(204)、章控阀值确定模块(205)、加速度计章动测量值计算模块(206)、自主章动控制计算模块(207)、平均推力系数确定模块(208)、章控/姿控选择模块(209)、控制力及力矩计算模块(210)、星体质量特性计算模块(211)、姿态运动动力学计算模块(212)、姿态脉冲时刻计算模块(213)、控制总冲累加计算模块(214)、轨道量测量计算模块(215)、轨道运动动力学计算模块(216)、遥测帧中断服务模块(217)、遥控指令接收模块(218)和时统中断服务模块(219)组成，其中，
仿真控制模块(11)用于接收并处理人机界面命令标识，负责初始化、启动相关功能模块、创建各功能事件、生成系统参数描述表和加载动力学仿真系统(2)，并在系统退出时结束系统的运行；
数据发送模块(12)用于发送数据接收模块(13)输出的信息，并将数据信息送监视显示模块(16)进行显示处理和记盘处理模块(17)进行记盘处理；
数据接收模块(13)用于接收管道数据，进行正确性检查和判断，并放入相应的缓冲区，处理后通知遥控处理模块(14)、遥测处理模块(15)、监视显示模块(16)和记盘处理模块(17)，或转发给动力学仿真系统(2)；所述数据接收模块(13)负责接收自旋稳定卫模平台、载人航天卫模平台、探索卫模平台的接口数据，转换为内部自定义格式的遥控数据，送遥控处理模块(14)进行处理；
遥控处理模块(14)用于实时接收和处理遥控数据帧、执行脉宽和控制参数，将相关遥控数据帧、执行脉宽和控制参数发送至动力学仿真系统(2)，由动力学仿真系统(2)模拟卫星的姿态变化过程；并将所有遥控指令、执行脉宽进行编码遥测反演计算，通过遥测信号模拟卫星各系统随遥控指令的变化过程；
遥测处理模块(15)用于接收动力学仿真系统(2)仿真的卫星编码遥测、模拟遥测信息以及遥控处理功能模块仿真的编码遥测信息，组装卫星编码遥测和模拟遥测信息帧，形成网发信息送数据发送模块(12)；
监视显示模块(16)用于监视显示本机软、硬件工作状态；
记盘处理模块(17)用于完成任务执行过程中的收发数据记盘工作；
数据输入模块用于在每个仿真周期的开始时，向系统注入发动机工作脉冲信息，四级发动机点火模拟模块(202)用于根据收到的四级脉冲，计算四级点火起始时刻和终止时刻；燃料量/压力/推力计算模块(203)用于计算发动机系统的剩余燃料量、贮箱压力以及相应发动机推力；日/月星历计算模块(204)用于根据经典的日、月星历的计算公式计算太阳、月亮位置；章控阀值确定模块(205)用于根据遥控指令章控阀值，确定阀值；加速度计章动测量值计算模块(206)用于模拟加速度计测量值经主动章控回路的带通滤波器后的输出，同时组装形成编码遥测信息；自主章动控制计算模块(207)用于进行粗、精控判断，形成控制脉冲前后沿，提供遥测输出及状态字设置；发动机平均推力系数确定模块(208)用于计算发动机平均推力系数；章控/姿控选择模块(209)用于根据遥控指令进行章控和姿控的选择设置；控制力及力矩计算模块(210)用于根据发动机的推力，计算力矩和力矩在惯性系投影，并计算本次工作周期的力冲量的累加值；星体质量特性计算模块(211)用于计算四级点火过程、二次分离后卫星的质量及惯量，计算二次分离后的惯量比和章动频率；姿态运动动力学计算模块(212)用于计算卫星角速度及姿态角；姿态脉冲时刻计算模块(213)用于根据卫星实时姿态、太阳星历和卫星轨道位置，计算姿态脉冲出现的时刻，并按格式组装成模拟遥测数据；控制总冲累加计算模块(214)用于本工作周期内的控制总冲的累加计算，提供给燃料量/压力/推力计算模块(203)进行控制总冲累加值计算；轨道量测量计算模块(215)用于测量船的位置及速度以及轨道量测量计算；轨道运动动力学计算模块(216)用于计算地心距、地球半张角以及卫星轨道位置的赤经赤纬值；遥测帧中断服务模块(217)用于向主程序输送遥测参数，并负责每秒显示刷新；遥控指令接收模块(218)用于接收遥控指令，并存储、显示；时统中断服务模块(219)用于接收时统中断板发送的定时中断信号，并向相应模块发出中断消息。
2.根据权利要求1所述的一种跨平台的自旋稳定卫星飞行仿真系统，其特征在于：所述数据接收模块(13)根据网收数据源地址，按信息分类在以下情况分别进行相应处理：在计算机演练模式下，对接收到的遥测遥控处理机数据，首先进行三判二处理，成功后，根据信息组装相应的返回信息帧和回答信息帧，并把信息交给遥控处理模块(14)处理，同时将处理信息转发动力学仿真系统(2)；在星地大回路演练模式下，收到卫星测控模拟平台帧中断时刻信息时，立即转发给动力学仿真系统(2)；收到通用卫星测控模拟平台的应答机遥测数据信息时，交给遥测处理模块(15)处理；其他数据信息立即交给遥控处理模块(14)处理，数据接收模块(13)利用管道技术，接收动力学仿真系统(2)的数据，按以下情况处理：在计算机演练模式下，将收到的全帧编码遥测数据，实时发送遥测通知事件，交给遥测处理模块(15)处理；将收到的全帧模拟遥测数据、返回的同控参数、姿章联控/应急章控、非同步控制单脉冲、非同步控制脉冲串信息转发数据发送模块(12)处理；在星地大回路演练模式下，将收到的编码遥测数据包，实时发送遥测通知事件，交给遥测处理功能模块(15)处理；将收到的模拟遥测数据包，将遥测时间增加一工作周期后，转发数据发送模块(12)处理。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610228224.7A CN105867167B (zh) | 2016-04-13 | 2016-04-13 | 跨平台的自旋稳定卫星飞行仿真系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610228224.7A CN105867167B (zh) | 2016-04-13 | 2016-04-13 | 跨平台的自旋稳定卫星飞行仿真系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105867167A CN105867167A (zh) | 2016-08-17 |
| CN105867167B true CN105867167B (zh) | 2018-12-11 |
Family
ID=56636696
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610228224.7A Expired - Fee Related CN105867167B (zh) | 2016-04-13 | 2016-04-13 | 跨平台的自旋稳定卫星飞行仿真系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105867167B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109087552B (zh) * | 2018-10-19 | 2022-10-04 | 西安基石睿盛信息技术有限公司 | 一种航天器运控模拟训练系统 |
| CN111177891B (zh) * | 2019-12-11 | 2023-08-22 | 上海卫星工程研究所 | 高轨转移段变轨策略仿真验证方法及系统 |
| CN111487065B (zh) * | 2020-04-24 | 2022-03-04 | 中国空间技术研究院 | 一种卫星控制系统及电推进系统闭环联试方法 |
| CN114861320B (zh) * | 2022-05-19 | 2023-02-10 | 北京航天飞行控制中心 | 一种航天器姿控推力建模及定轨解算方法 |
| CN115421543B (zh) * | 2022-11-02 | 2023-05-16 | 北京宇航系统工程研究所 | 一种低温贮箱压力控制方法及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7584085B2 (en) * | 2000-08-25 | 2009-09-01 | The Boeing Company | Method and apparatus for simulating multiple ground station ranging for spacecraft |
| CN102354123A (zh) * | 2011-07-18 | 2012-02-15 | 北京航空航天大学 | 一种跨平台可扩展的卫星动态仿真测试系统 |
| CN103760898A (zh) * | 2014-02-12 | 2014-04-30 | 航天东方红卫星有限公司 | 一种快速建立小卫星控制测试系统的方法 |
| CN104260908A (zh) * | 2014-08-28 | 2015-01-07 | 上海微小卫星工程中心 | 一种跨平台对地观测卫星联合演示验证系统 |
| CN105468817A (zh) * | 2015-11-12 | 2016-04-06 | 上海科梁信息工程股份有限公司 | 一种多模型实时仿真系统 |
-
2016
- 2016-04-13 CN CN201610228224.7A patent/CN105867167B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7584085B2 (en) * | 2000-08-25 | 2009-09-01 | The Boeing Company | Method and apparatus for simulating multiple ground station ranging for spacecraft |
| CN102354123A (zh) * | 2011-07-18 | 2012-02-15 | 北京航空航天大学 | 一种跨平台可扩展的卫星动态仿真测试系统 |
| CN103760898A (zh) * | 2014-02-12 | 2014-04-30 | 航天东方红卫星有限公司 | 一种快速建立小卫星控制测试系统的方法 |
| CN104260908A (zh) * | 2014-08-28 | 2015-01-07 | 上海微小卫星工程中心 | 一种跨平台对地观测卫星联合演示验证系统 |
| CN105468817A (zh) * | 2015-11-12 | 2016-04-06 | 上海科梁信息工程股份有限公司 | 一种多模型实时仿真系统 |
Non-Patent Citations (4)
| Title |
|---|
| 受控航天器的通用仿真框架;郑伟 等;《航天控制》;20010131(第1期);全文 * |
| 基于分层体系结构的卫模仿真软件通用化设计;饶爱水 等;《兵工自动化》;20110831(第8期);全文 * |
| 航天器控制系统通用仿真技术研究;李青松;《中国优秀硕士学位论文全文数据库信息科技辑》;20100515(第05期);全文 * |
| 航天测控系统通用仿真平台;闫高宗;《系统仿真学报》;20060831;第18卷;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105867167A (zh) | 2016-08-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105867167B (zh) | 跨平台的自旋稳定卫星飞行仿真系统 | |
| CN101464935B (zh) | 基于网络的auv智能容错组合导航仿真系统 | |
| CN100575877C (zh) | 基于多信息融合的航天器组合导航方法 | |
| CN101320524B (zh) | 多处理器实时仿真平台 | |
| CN110262280A (zh) | 航天器交会对接飞行控制智能数据分析与决策支持系统 | |
| CN104898642B (zh) | 一种用于航天器姿态控制算法的集成测试仿真系统 | |
| CN101876710B (zh) | 基于实时卫星信号模拟器同步增减可见卫星装置及构建法 | |
| CN106840196A (zh) | 一种捷联惯性导航计算机测试系统及实现方法 | |
| CN100476359C (zh) | 一种基于星光角距的深空探测器upf自主天文导航方法 | |
| CN106338753A (zh) | 一种基于地面站/星间链路/gnss联合测量的地球同步轨道星座定轨方法 | |
| CN101354251B (zh) | 一种深空探测器等效转移轨道确定方法 | |
| CN106019323A (zh) | 卫星导航接收机仿真测试系统 | |
| Acharya | Understanding satellite navigation | |
| CN102591343A (zh) | 基于两行根数的卫星轨道维持控制方法 | |
| Hill | Autonomous navigation in libration point orbits | |
| CN107589756A (zh) | 一种奔月卫星编队初始化方法 | |
| CN104848862A (zh) | 一种环火探测器精密同步定位守时方法及系统 | |
| CN108075826A (zh) | 一种实现卫星半物理仿真试验星地时间同步的方法及系统 | |
| CN111731513B (zh) | 一种基于单脉冲轨控的高精度引力场中回归轨道维持方法 | |
| Cork | Aircraft dynamic navigation for unmanned aerial vehicles | |
| CN105865432B (zh) | 一种针对陀螺仪多源噪声的混合滤波方法与测试平台 | |
| CN105487405A (zh) | 低低跟踪重力测量卫星半物理仿真系统 | |
| US6282467B1 (en) | Three-axis inertial attitude determination for spinning spacecraft | |
| CN107054702B (zh) | 一种半物理仿真测试中地磁场模拟方法 | |
| Di Domenico | Development of a hardware-in-the-loop simulation framework for interplanetary transfers on smaller timescales |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20181026 Address after: 214400 103 mailbox 850, Jiangyin City, Wuxi, Jiangsu Applicant after: CHINA SATELITE MARITIME MEASUREMENT AND CONTROL BUREAU Address before: 214400 103 mailbox 850, Jiangyin City, Wuxi, Jiangsu Applicant before: UNIT 63680 OF PLA |
|
| 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: 20181211 |