CN113221267A - 一种基于在轨数据的发动机性能参数修正方法 - Google Patents
一种基于在轨数据的发动机性能参数修正方法 Download PDFInfo
- Publication number
- CN113221267A CN113221267A CN202110477761.6A CN202110477761A CN113221267A CN 113221267 A CN113221267 A CN 113221267A CN 202110477761 A CN202110477761 A CN 202110477761A CN 113221267 A CN113221267 A CN 113221267A
- Authority
- CN
- China
- Prior art keywords
- engine
- orbital transfer
- orbit
- flow
- rail
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/12—Timing analysis or timing optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Plasma Technology (AREA)
Abstract
本发明提出的一种基于在轨数据的发动机性能参数修正方法,利用在轨实际表现数据,对发动机的性能参数进行修正,有益于提高在轨数据的分析水平,更加深入了解发动机的工作性能;本发明的方法无需多余的硬件,即可实现对发动机性能参数的修正,不仅有益于提高在轨剩余推进剂的估算精度,进一步提高在轨推力的预测精度和变轨效率,避免推进剂的浪费,对后续卫星的长期在轨管理具有重大的意义。
Description
技术领域
本发明涉及一种发动机性能参数修正方法。
背景技术
我国地球静止轨道卫星一般采用双组元统一推进系统,并配置远地方发动机和姿控推力器。远地点发动机主要用于为航天器轨道机动入轨提供推力。
在卫星轨道机动期间,制定变轨策略需要用到整星质量、发动机的比冲及推力等参数。其中发动机比冲可以通过发动机地面热标试验得到,发动机推力和整星质量则需要通过引入发动机的小偏差方程。因此,整星的变轨控制效果及长期在轨管理均与发动机参数的精确性密切相关。如果发动机的实际性能偏差过大,不仅会影响变轨次数,耗费不必要的人力、物力,更会造成推进剂的浪费,从而导致卫星的寿命减少,带来不可避免的损失。
由于,发动机性能参数一般由额定点附近的小偏差方程给出,因此,在参数处理的方法本身存在一定的系统误差,使用受限。另外,发动机在轨工作的工况与地面热标工况不尽相同,且地面热标精度有待提高,冷试与热试之间的参数转换也存在一定的方法误差。因此,发动机的实际在轨性能与地面相比,往往存在不可忽略的差异性。
发明内容
本发明要解决的技术问题是:提供了一种发动机性能在轨修正方法,利用在轨数据,结合参数寻优方法,实现了对发动机性能参数进行修正的目的,适用于双组元推进系统卫星发动机的在轨性能修正。
本发明所采用的技术方案是:一种基于在轨数据的发动机性能参数修正方法,包括如下步骤:
(1)每次卫星变轨前,打开卫星配置的超声波流量计,测量发动机工作过程中的推进剂流量,轨控发动机点火结束后,关闭超声波流量计,流量测量停止;每次变轨期间,对超声波流量计的输出进行积分得到该次变轨过程中的推进剂消耗量Δmui;
其中,Qo、Qf分别为氧化剂和燃烧剂贮箱下游对应的超声波流量计输出流量,ti为第i次轨道机动的时间;i为正整数;
(2)记录每次卫星变轨的沉底推力器点火时间,结合沉底推力器的流量,由BK法(薄记法)获取本次变轨期间沉底推力器的推进剂消耗量Δmci;
Δmci=ωc×tci
tci为第i次变轨对应的沉底推力器的点火时间,ωc为沉底推力器的流量。
(3)记录每次卫星变轨期间姿控推力器的点火时间,结合姿控推力器的流量,由BK法获取本次变轨期间的姿控推力器的推进剂消耗量Δmzi;
Δmzi=ωz×tzi
tzi为第i次变轨对应的姿控推力器的点火时间,ωz为姿控推力器的流量。
(4)结合步骤(1)、步骤(2)和步骤(3)的计算结果,获得第i次变轨期间的轨控发动机实际消耗量Δmei;
Δmei=Δmui-Δmci-Δmzi;
(5)通过第i次变轨结束后由轨道标定得到的推进剂消耗量数据Δmoi,获得第i次变轨的速度增量Δvi;
Ispt为发动机在地面标定的额定比冲,Moi为第i次变轨前的卫星质量。
(6)将第i次变轨前的卫星质量Moi及步骤(5)计算得到的速度增量Δvi,代入火箭方程,获得卫星变轨发动机比冲Isp对应的轨控发动机推进剂消耗量ΔMei;
(7)设置目标函数J=|ΔMei-Δmei|,采用单变量寻优法,求解J(Isp)=0的最优解Isp0,Isp0为修正的卫星变轨发动机比冲;
(8)由固定比冲法,结合步骤(7)得到的标定比冲Isp0和步骤(4)得到的轨控发动机消耗量Δmei,获得第i次变轨修正后的推力Fi;
(9)由步骤(1)中超声波流量计的输出流量对变轨期间的发动机平均流量进行标定,即
(10)根据第i次变轨过程中的平均压力和温度遥测数据,使用轨控发动机的小偏差方程计算该次变轨的氧化剂流量ωo和燃烧剂流量ωf;
氧化剂输送管路和燃烧剂输送管路流量的计算方程为:
ωo=a11(Po-ΔPlo)+b11(Pf-ΔPlf)+c11To+d11Tf+e11,
ωf=a12(Po-ΔPlo)+b12(Pf-ΔPlf)+c12To+d12Tf+e12,
式中,系数a11、b11、c11、d11、e11、a12、b12、c12、d12、e12均通过轨控发动机地面试验获得;Po、To分别为氧化剂贮箱的遥测平均压力、平均温度;Pf、Tf分别为燃烧剂贮箱的遥测平均压力、平均温度;ΔPlo、ΔPlf分别为待修正的氧化剂输送管路和燃烧剂输送管路流阻。
ΔPlo0、ΔPlf0分别为修正的氧化剂输送管路和燃烧剂输送管路流阻;
将ΔPlo0、ΔPlf0分别带入步骤(10)的氧化剂输送管路和燃烧剂输送管路流量的计算方程中,得到修正后的轨控发动机的小偏差方程,用于变轨过程中的发动机推力及推进剂消耗量计算。
本发明与现有技术相比的优点在于:
(1)本发明提出的基于在轨数据的发动机性能参数修正方法,利用在轨实际表现数据,对发动机的性能参数进行修正,有益于提高在轨数据的分析水平,更加深入了解发动机的工作性能;
(2)本发明的方法无需多余的硬件,即可实现对发动机性能参数的修正,不仅有益于提高在轨剩余推进剂的估算精度,进一步提高在轨推力的预测精度和变轨效率,避免推进剂的浪费,对后续卫星的长期在轨管理具有重大的意义;
(3)本发明的方法通用性强,适用于常规双组元推进系统单台轨控发动机的参数修正,也可推广至多台轨控发动机的在轨参数修正,具有广泛的应用价值和推广前景。
附图说明
图1为本发明的流程图。
具体实施方式
下面结合附图对本发明进行说明。
如图1所示,一种基于在轨数据的发动机性能参数修正方法,包括如下步骤:
(1)每次卫星变轨前,打开卫星配置的超声波流量计,测量发动机工作过程中的推进剂流量,轨控发动机点火结束后,关闭超声波流量计,流量测量停止;每次变轨期间,对超声波流量计的输出进行积分得到该次变轨过程中的推进剂消耗量Δmui;
其中,Qo、Qf分别为氧化剂和燃烧剂贮箱下游对应的超声波流量计输出流量,ti为第i次轨道机动的时间;i为正整数;
(2)记录每次卫星变轨的沉底推力器点火时间,结合沉底推力器的流量,由BK法(薄记法)获取本次变轨期间沉底推力器的推进剂消耗量Δmci;
Δmci=ωc×tci
tci为第i次变轨对应的沉底推力器的点火时间,ωc为沉底推力器的流量。
(3)记录每次卫星变轨期间姿控推力器的点火时间,结合姿控推力器的流量,由BK法获取本次变轨期间的姿控推力器的推进剂消耗量Δmzi;
Δmzi=ωz×tzi
tzi为第i次变轨对应的姿控推力器的点火时间,ωz为姿控推力器的流量。
(4)结合步骤(1)、步骤(2)和步骤(3)的计算结果,获得第i次变轨期间的轨控发动机实际消耗量Δmei;
Δmei=Δmui-Δmci-Δmzi;
(5)通过第i次变轨结束后由轨道标定得到的推进剂消耗量数据Δmoi,获得第i次变轨的速度增量Δvi;
Ispt为发动机在地面标定的额定比冲,Moi为第i次变轨前的卫星质量。
(6)将第i次变轨前的卫星质量Moi及步骤(5)计算得到的速度增量Δvi,代入火箭方程,获得卫星变轨发动机比冲Isp对应的轨控发动机推进剂消耗量ΔMei;
(7)设置目标函数J=|ΔMei-Δmei|,采用单变量寻优法,求解J(Isp)=0的最优解Isp0,Isp0为修正的卫星变轨发动机比冲;
(8)由固定比冲法,结合步骤(7)得到的标定比冲Isp0和步骤(4)得到的轨控发动机消耗量Δmei,获得第i次变轨修正后的推力Fi;
(9)由步骤(1)中超声波流量计的输出流量对变轨期间的发动机平均流量进行标定,即
(10)根据第i次变轨过程中的平均压力和温度遥测数据,使用轨控发动机的小偏差方程计算该次变轨的氧化剂流量ωo和燃烧剂流量ωf;
氧化剂输送管路和燃烧剂输送管路流量的计算方程为:
ωo=a11(Po-ΔPlo)+b11(Pf-ΔPlf)+c11To+d11Tf+e11,
ωf=a12(Po-ΔPlo)+b12(Pf-ΔPlf)+c12To+d12Tf+e12,
式中,系数a11、b11、c11、d11、e11、a12、b12、c12、d12、e12均通过轨控发动机地面试验获得;Po、To分别为氧化剂贮箱的遥测平均压力、平均温度;Pf、Tf分别为燃烧剂贮箱的遥测平均压力、平均温度;ΔPlo、ΔPlf分别为待修正的氧化剂输送管路和燃烧剂输送管路流阻。
ΔPlo0、ΔPlf0分别为修正的氧化剂输送管路和燃烧剂输送管路流阻;
将ΔPlo0、ΔPlf0分别带入步骤(10)的氧化剂输送管路和燃烧剂输送管路流量的计算方程中,得到修正后的轨控发动机的小偏差方程,用于变轨过程中的发动机推进剂消耗量计算。
本发明未详细说明部分属于本领域技术人员公知技术。
Claims (10)
1.一种基于在轨数据的发动机性能参数修正方法,其特征在于,包括如下步骤:
(1)每次卫星变轨前,打开卫星配置的超声波流量计,测量发动机工作过程中的推进剂流量,轨控发动机点火结束后,关闭超声波流量计,流量测量停止;每次变轨期间,对超声波流量计的输出进行积分得到该次变轨过程中的推进剂消耗量Δmui;
(2)记录每次卫星变轨的沉底推力器点火时间,结合沉底推力器的流量,由BK法获取本次变轨期间沉底推力器的推进剂消耗量Δmci;
(3)记录每次卫星变轨期间姿控推力器的点火时间,结合姿控推力器的流量,由BK法获取本次变轨期间的姿控推力器的推进剂消耗量Δmzi;
(4)结合步骤(1)、步骤(2)和步骤(3)的计算结果,计算获得第i次变轨期间的轨控发动机实际消耗量Δmei;
(5)通过第i次变轨结束后由轨道标定得到的推进剂消耗量数据Δmoi,获得第i次变轨的速度增量Δvi;
(6)将第i次变轨前的卫星质量Moi及步骤(5)计算得到的速度增量Δvi,代入火箭方程,获得卫星变轨发动机比冲Isp对应的轨控发动机推进剂消耗量ΔMei;
(7)设置目标函数J=|ΔMei-Δmei|,采用单变量寻优法,求解J(Isp)=0的最优解Isp0,Isp0为修正的卫星变轨发动机比冲;
(8)由固定比冲法,结合步骤(7)得到的标定比冲Isp0和步骤(4)得到的轨控发动机消耗量Δmei,获得第i次变轨修正后的推力Fi;
(10)根据第i次变轨过程中的平均压力和温度遥测数据,使用轨控发动机的小偏差方程计算该次变轨的氧化剂流量ωo和燃烧剂流量ωf;
ΔPlo、ΔPlf分别为待修正的氧化剂输送管路和燃烧剂输送管路流阻,ΔPlo0、ΔPlf0分别为修正的氧化剂输送管路和燃烧剂输送管路流阻;
将ΔPlo0、ΔPlf0分别带入步骤(10)中,得到修正后的轨控发动机的小偏差方程,用于变轨过程中的发动机推力及推进剂消耗量计算。
3.根据权利要求1或2所述的一种基于在轨数据的发动机性能参数修正方法,其特征在于,步骤(2)中,Δmci=ωc×tci;
tci为第i次变轨对应的沉底推力器的点火时间,ωc为沉底推力器的流量。
4.根据权利要求3所述的一种基于在轨数据的发动机性能参数修正方法,其特征在于,步骤(3)中,Δmzi=ωz×tzi;
tzi为第i次变轨对应的姿控推力器的点火时间,ωz为姿控推力器的流量。
5.根据权利要求4所述的一种基于在轨数据的发动机性能参数修正方法,其特征在于,步骤(4)中,Δmei=Δmui-Δmci-Δmzi。
10.根据权利要求9所述的一种基于在轨数据的发动机性能参数修正方法,其特征在于,步骤(10)中,氧化剂输送管路和燃烧剂输送管路流量的计算方程为:
ωo=a11(Po-ΔPlo)+b11(Pf-ΔPlf)+c11To+d11Tf+e11,
ωf=a12(Po-ΔPlo)+b12(Pf-ΔPlf)+c12To+d12Tf+e12,
式中,系数a11、b11、c11、d11、e11、a12、b12、c12、d12、e12均通过轨控发动机地面试验获得;Po、To分别为氧化剂贮箱的遥测平均压力、平均温度;Pf、Tf分别为燃烧剂贮箱的遥测平均压力、平均温度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110477761.6A CN113221267B (zh) | 2021-04-30 | 2021-04-30 | 一种基于在轨数据的发动机性能参数修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110477761.6A CN113221267B (zh) | 2021-04-30 | 2021-04-30 | 一种基于在轨数据的发动机性能参数修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113221267A true CN113221267A (zh) | 2021-08-06 |
CN113221267B CN113221267B (zh) | 2023-06-16 |
Family
ID=77090101
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110477761.6A Active CN113221267B (zh) | 2021-04-30 | 2021-04-30 | 一种基于在轨数据的发动机性能参数修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113221267B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113989954A (zh) * | 2021-10-26 | 2022-01-28 | 上海空间推进研究所 | 上面级主发动机在轨推力和推进剂剩余量计算方法及系统 |
CN114117750A (zh) * | 2021-11-08 | 2022-03-01 | 中国长城工业集团有限公司 | 一种化学推进与电推进混合配置的geo卫星推进剂预算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5934621A (en) * | 1997-10-08 | 1999-08-10 | Kent; Stan | Batched optimized method for transfer orbit construction |
CN106114911A (zh) * | 2016-06-23 | 2016-11-16 | 中国空间技术研究院 | 一种卫星推进系统在轨飞行混合比调整方法 |
CN106198033A (zh) * | 2016-06-23 | 2016-12-07 | 中国空间技术研究院 | 一种卫星发动机比冲在轨计算方法 |
CN111688954A (zh) * | 2020-07-02 | 2020-09-22 | 北京空间技术研制试验中心 | 一种航天器贮箱系统的排放不平衡系数在轨估计方法 |
-
2021
- 2021-04-30 CN CN202110477761.6A patent/CN113221267B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5934621A (en) * | 1997-10-08 | 1999-08-10 | Kent; Stan | Batched optimized method for transfer orbit construction |
CN106114911A (zh) * | 2016-06-23 | 2016-11-16 | 中国空间技术研究院 | 一种卫星推进系统在轨飞行混合比调整方法 |
CN106198033A (zh) * | 2016-06-23 | 2016-12-07 | 中国空间技术研究院 | 一种卫星发动机比冲在轨计算方法 |
CN111688954A (zh) * | 2020-07-02 | 2020-09-22 | 北京空间技术研制试验中心 | 一种航天器贮箱系统的排放不平衡系数在轨估计方法 |
Non-Patent Citations (2)
Title |
---|
刘格军: "基于视加速度的固体火箭发动机飞行比冲计算", 固体火箭技术, vol. 31, no. 5, pages 464 - 465 * |
曹伟 等: "嫦娥三号推进系统在轨推进剂耗量计算方法", 火箭推进, vol. 41, no. 01, pages 63 - 69 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113989954A (zh) * | 2021-10-26 | 2022-01-28 | 上海空间推进研究所 | 上面级主发动机在轨推力和推进剂剩余量计算方法及系统 |
CN114117750A (zh) * | 2021-11-08 | 2022-03-01 | 中国长城工业集团有限公司 | 一种化学推进与电推进混合配置的geo卫星推进剂预算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113221267B (zh) | 2023-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113221267A (zh) | 一种基于在轨数据的发动机性能参数修正方法 | |
CN111613818B (zh) | 一种燃料电池氢系统的剩余氢气计算方法 | |
CN110525694B (zh) | 一种考虑并联贮箱压力差异的推进系统混合比控制方法 | |
CN111426364B (zh) | 一种基于多源数据融合的推进剂剩余量测量方法 | |
GB2500755B (en) | Fluid mass flow measurement apparatus and method | |
CN112109923B (zh) | 一种基于自主标定的卫星半长轴保持自主轨控方法 | |
CN115562325A (zh) | 一种通过姿态偏置利用面外推力器实现离轨控制方法 | |
CN108791957B (zh) | 一种卫星变轨期间横向质心高精度实时辨识方法 | |
US20170328567A1 (en) | Multivariable fuel control and estimator (mfce) for preventing combustor blowout | |
CN111470075B (zh) | 一种基于人工智能算法的航天器在轨推力预测方法 | |
CN106198033B (zh) | 一种卫星发动机比冲在轨计算方法 | |
US11598289B1 (en) | Propellant volume and mixture ratio control | |
CN108563914B (zh) | 基于夏式最小二乘的轨道控制推力拟合系数计算方法 | |
US5934621A (en) | Batched optimized method for transfer orbit construction | |
US10415484B2 (en) | Engine control system | |
CN111290433B (zh) | 一种长期自主编队联合管道保持方法 | |
CN113218660B (zh) | 一种电推力矢量在轨标定方法 | |
CN114035611A (zh) | 可重复使用高超声速飞行器上升段轨迹优化与制导方法 | |
CN111409864B (zh) | 一种编队卫星系统星间气动外形差异在轨标定方法 | |
JPH09287550A (ja) | イオンエンジン装置 | |
Matthijssen et al. | State-of-the-Art Gauging Components for Improved Propellant Management on 3-Axis Stabilized Spacecraft | |
RU2388665C1 (ru) | Способ калибровки импульса тяги маршевого двигателя | |
CN113968361B (zh) | 一种适用地球同步卫星定点控制规划的解析计算方法 | |
CN113944570B (zh) | 一种基于电动泵发动机的轨道控制方法及计算设备 | |
RU2714286C1 (ru) | Способ приведения космического аппарата к долготе стояния на геостационарной орбите |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |