CN114021253A - 一种基于Modelica语言的液体火箭发动机动态仿真方法 - Google Patents

一种基于Modelica语言的液体火箭发动机动态仿真方法 Download PDF

Info

Publication number
CN114021253A
CN114021253A CN202111295236.9A CN202111295236A CN114021253A CN 114021253 A CN114021253 A CN 114021253A CN 202111295236 A CN202111295236 A CN 202111295236A CN 114021253 A CN114021253 A CN 114021253A
Authority
CN
China
Prior art keywords
liquid
gas
interface
rocket engine
phase
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.)
Pending
Application number
CN202111295236.9A
Other languages
English (en)
Inventor
王丹
陈宏玉
任孝文
刘占一
周晨初
李舒欣
周康
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Aerospace Propulsion Institute
Original Assignee
Xian Aerospace Propulsion Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Aerospace Propulsion Institute filed Critical Xian Aerospace Propulsion Institute
Priority to CN202111295236.9A priority Critical patent/CN114021253A/zh
Publication of CN114021253A publication Critical patent/CN114021253A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及液体火箭发动机系统仿真,主要涉及液体火箭发动机动态特性仿真方法。本发明的目的是解决现有液体火箭发动机的动态特性仿真方法中缺乏能够涵盖各类组件模型传输需求的接口定义方式的技术问题,提供一种基于Modelica语言的液体火箭发动机动态仿真方法。该方法采用模块化建模思想建模,对不同组件的模型进行切分,使得各组件具有各自独立的组件仿真模型,采用基于Modelica语言的液体火箭发动机动态仿真接口定义,然后通过基于Modelica语言的接口将组件相连,由于接口的定义满足组件间物质和能量传输的需求,从而实现了液体火箭发动机动态特性仿真。

Description

一种基于Modelica语言的液体火箭发动机动态仿真方法
技术领域
本发明涉及液体火箭发动机系统仿真,主要涉及一种基于Modelica语言的液体火箭发动机动态仿真方法。
背景技术
液体火箭发动机起动关机过程涉及机械、燃烧和传热等复杂的物理化学过程,采用模块化建模思想,构建液体火箭发动机动态仿真模型库,以拖拽式建模方式实现系统模型搭建,进行全系统仿真,是发动机型号方案改进和新型号发动机研制的重要手段之一。然而,采用模块化建模思想建模,需要对不同组件的模型进行切分,并通过接口将组件相连,接口的定义需满足组件间物质和能量传输的需求。因此,着眼于液体火箭发动机的动态特性仿真方法,亟需一种能够涵盖各类组件模型传输需求的接口定义方式。
发明内容
本发明的目的是解决现有液体火箭发动机的动态特性仿真方法中缺乏能够涵盖各类组件模型传输需求的接口定义方式的技术问题,提供一种基于Modelica语言的液体火箭发动机动态仿真方法。
为解决上述技术问题,本发明提供的技术解决方案如下:
一种基于Modelica语言的液体火箭发动机动态仿真方法,其特殊之处在于,包括以下步骤:
1)通过模块化建模方式,将液体火箭发动机系统中主要组件的动态仿真模型各自封装,使得各组件具有各自独立的组件仿真模型;
2)根据液体火箭发动机中使用的介质,基于Modelica语言定义液体火箭发动机系统各类组件仿真模型中存在数据传输需求的液体接口、气体接口、燃气接口、两相接口、机械转动接口和控制接口六类通用接口;所述介质包括液体和气体;所述通用接口采用标准接口;
其中,
A)气体接口用于流经纯气体介质的组件,共传递三个变量,包括质量流量m_flow、绝对压力p和温度T;
B)液体接口用于流经单相液体介质的组件,共传递两个变量,包括质量流量m_flow和绝对压力p;
C)两相接口用于流经气液混合介质的组件,共传递八个变量,包括绝对压力p、液相温度Tl、气相温度Tg、液相体积份数Vlbar、推进剂气相占比Vgbar、液相质量流量m_flow_l、气相质量流量m_flow_g和吹除气体质量流量m_flow_ig;
D)燃气接口用于热力组件,共传递七个变量,包括混合比Km、液体燃料流量m_flowlf、液体氧化剂流量m_flowlo、液相体积份数Vlbar,以及与气体接口相同的燃气质量流量m_flowg、绝对压力p和温度T;
E)机械转动接口位于涡轮与泵之间,用于涡轮与泵之间的机械力传输,共传递两个变量,包括旋转角度phi和扭矩tau;
F)控制接口传递的变量为实数信号,用于阀门和/或调节器控制信号的输入,以控制阀门的开启关闭和/或调节器的调节状态;
3)将步骤1)的各组件仿真模型,利用步骤2)定义的相应接口相连,通过接口传递数据;
4)对步骤3)通过接口连接各组件后的液体火箭发动机系统,进行基于Modelica语言的液体火箭发动机动态仿真。
进一步地,步骤A)中,所述流经纯气体介质的组件包括气瓶、气体管路、气体节流孔和/或气体阀门。
进一步地,步骤B)中,所述流经单相液体介质的组件包括液体管路、液体节流孔、液体阀门、离心泵和/或液涡轮。
进一步地,步骤C)中,所述流经气液混合介质的组件包括低温介质充填管路、头腔和/或低温充填冷却套。
进一步地,步骤D)中,所述热力组件包括发生器、燃气导管、燃气涡轮和/或推力室。
进一步地,步骤B)中,所述单相液体介质为氧化剂和/或燃料。
进一步地,步骤1)中,所述头腔为乳化吹除头腔。
本发明相比现有技术具有的有益效果如下:
本发明提供的一种基于Modelica语言的液体火箭发动机动态仿真方法,采用模块化建模思想建模,对不同组件的模型进行切分,使得各组件具有各自独立的组件仿真模型,采用基于Modelica语言的液体火箭发动机动态仿真接口定义,然后通过基于Modelica语言的接口将组件相连,由于接口的定义满足组件间物质和能量传输的需求,从而实现了液体火箭发动机动态特性仿真。
附图说明
图1为本发明实施例中液体火箭发动机系统示意图。
具体实施方式
下面结合实施例对本发明作进一步地说明。
一种基于Modelica语言的液体火箭发动机动态仿真方法,包括以下步骤:
1)通过模块化建模方式,将液体火箭发动机系统中主要组件气瓶、气孔板、贮箱、泵、管路、阀门、头腔、发生器、涡轮、燃气管路和燃气孔板等的动态仿真模型各自封装,使得各组件具有各自独立的组件仿真模型
2)根据液体火箭发动机中使用的介质(液体、气体),涉及的领域(流体、机械、控制),基于Modelica语言定义完全覆盖液体火箭发动机系统各类组件仿真模型(组件模型)中存在数据传输需求的液体接口、气体接口、燃气接口、两相接口、机械转动接口和控制接口六类通用接口,覆盖液体火箭发动机系统仿真各类组件仿真模型的数据传输需求;组件仿真模型均依托标准接口(即通用接口)传递数据,以实现不同组件仿真模型在统一平台下的数据贯通;
其中,
A)气体接口用于流经纯气体介质的组件,共传递三个变量,包括质量流量m_flow、绝对压力p和温度T,增加温度是因为动态仿真过程中需考虑气体的可压缩性;所述流经纯气体介质的组件包括气瓶、气体管路、气体节流孔、气体阀门等;比如气瓶、气孔板、贮箱入口的工作介质均为氮气,通过气体接口连接,传递氮气的质量流量m_flow、绝对压力p和温度T;
B)液体接口用于流经单相液体介质的组件,共传递两个变量,包括质量流量m_flow和绝对压力p;所述流经单相液体介质的组件包括液体管路、液体节流孔、液体阀门、离心泵、液涡轮等;所述单相液体介质为氧化剂和/或燃料;比如贮箱出口、泵、管路、阀门、头腔入口的工作介质均为单液相推进剂(氧化剂/燃料),通过液体接口连接,传递推进剂的质量流量m_flow和绝对压力p;
C)两相接口用于流经气液混合介质的组件,共传递八个变量,包括绝对压力p、液相温度Tl、气相温度Tg、液相体积份数Vlbar、推进剂气相占比Vgbar、液相质量流量m_flow_l、气相质量流量m_flow_g和吹除气体质量流量m_flow_ig;所述流经气液混合介质的组件包括低温介质充填管路、带吹除头腔、低温充填冷却套等(具体而言:由于组件介质为混合介质,因此需传递液相质量流量m_flow_l、气相质量流量m_flow_g、吹除气体质量流量m_flow_ig;与其他接口类似,也需要传递绝对压力p,而温度则分为液相温度Tl和气相温度Tg;考虑充填过程,需要传递液相体积份数Vlbar、推进剂气相占比Vgbar,共传递八个变量);比如燃料头腔为乳化吹除头腔,需同时通过液体接口与上游阀门连接,气体接口与吹除气孔板连接,在头腔内形成气液两相掺混,因此出口为两相接口,传递液相质量流量m_flow_l、气相质量流量m_flow_g、吹除气体质量流量m_flow_ig、绝对压力p、液相温度Tl、气相温度Tg、液相体积份数Vlbar和推进剂气相占比Vgbar;
D)燃气接口用于热力组件,共传递七个变量,包括混合比Km、液体燃料流量m_flowlf、液体氧化剂流量m_flowlo、液相体积份数Vlbar,以及与气体接口相同的燃气质量流量m_flowg、绝对压力p和温度T(具体而言:燃气接口首先需传递与气体接口相同的燃气质量流量m_flowg、绝对压力p和温度T;另外,由于燃气接口主要用于热力组件,如发生器、燃气导管、燃气涡轮和推力室,在热力组件中存在燃气、液态燃料、气态燃料、吹除惰性气体多种组份掺混和反应的过程,因此增加了混合比Km、液体燃料流量m_flowlf、液体氧化剂流量m_flowlo、液相体积份数Vlbar,共传递七个变量);比如发生器入口通过两相接口与头腔相连,内部进行燃烧过程模拟,因此出口为燃气接口,后续连接的涡轮、燃气管路、燃气孔板均通过燃气接口传递数据;燃气接口传递燃气质量流量m_flowg、绝对压力p、温度T、混合比Km、液体燃料流量m_flowlf、液体氧化剂流量m_flowlo和液相体积份数Vlbar;
E)机械转动接口位于涡轮与泵之间,用于涡轮与泵之间的机械力传输,共传递两个变量,包括旋转角度phi和扭矩tau;
F)控制接口传递的变量为实数信号,用于阀门和/或调节器控制信号的输入,以控制阀门的开启关闭和/或调节器的调节状态;比如阀门的开启关闭由控制信号输入,通过控制接口与输入信号连接;
3)将步骤1)的各组件仿真模型,利用步骤2)定义的相应接口相连,通过接口传递数据;
4)对步骤3)通过接口将各组件相连的液体火箭发动机系统,进行基于Modelica语言的液体火箭发动机动态仿真。
一种基于Modelica语言的液体火箭发动机动态仿真方法,基于文中所示的接口构建规则,可以满足火箭发动机动态仿真的各类组件数据传输需求。图1所示的液体火箭发动机系统示意图(发生器循环系统图)为一个实施例,并不代表只能采用这种连接方式。比如说:液体接口可以连管路,也可以连泵、阀门、头腔等,只要组件的工作介质是单相流体,就能够通过液体接口相连,并且只需要传递压力和流量两个变量即可。
列表中变量类型带有flow前缀的代表流变量,其特征为在多个接口连接处各接口该变量值之和为零,典型流变量为流量,其余变量为势变量,其特征为在多个接口连接处各接口该变量值相等,典型势变量为压力。通过流变量和势变量的划分,使得多个组件在同一节点处汇集时,能够自动匹配出各自变量的关系,拓展了模型对于火箭发动机中复杂管网仿真的支持功能。下表为接口变量列表:
表1接口变量列表
Figure BDA0003336303410000051
Figure BDA0003336303410000061
本发明的基于Modelica语言的液体火箭发动机动态仿真方法不仅适用于全系统,还适用于半系统,以及不同循环方式。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制,对于本领域的普通专业技术人员来说,可以对前述实施例所记载的具体技术方案进行修改,或者对其中部分技术特征进行等同替换,而这些修改或者替换,并不使相应技术方案的本质脱离本发明所保护技术方案的范围。

Claims (7)

1.一种基于Modelica语言的液体火箭发动机动态仿真方法,其特征在于,包括以下步骤:
1)通过模块化建模方式,将液体火箭发动机系统中主要组件的动态仿真模型各自封装,使得各组件具有各自独立的组件仿真模型;
2)根据液体火箭发动机中使用的介质,基于Modelica语言定义液体火箭发动机系统各类组件仿真模型中存在数据传输需求的液体接口、气体接口、燃气接口、两相接口、机械转动接口和控制接口六类通用接口;所述介质包括液体和气体;所述通用接口采用标准接口;
其中,
A)气体接口用于流经纯气体介质的组件,共传递三个变量,包括质量流量m_flow、绝对压力p和温度T;
B)液体接口用于流经单相液体介质的组件,共传递两个变量,包括质量流量m_flow和绝对压力p;
C)两相接口用于流经气液混合介质的组件,共传递八个变量,包括绝对压力p、液相温度Tl、气相温度Tg、液相体积份数Vlbar、推进剂气相占比Vgbar、液相质量流量m_flow_l、气相质量流量m_flow_g和吹除气体质量流量m_flow_ig;
D)燃气接口用于热力组件,共传递七个变量,包括混合比Km、液体燃料流量m_flowlf、液体氧化剂流量m_flowlo、液相体积份数Vlbar,以及与气体接口相同的燃气质量流量m_flowg、绝对压力p和温度T;
E)机械转动接口位于涡轮与泵之间,用于涡轮与泵之间的机械力传输,共传递两个变量,包括旋转角度phi和扭矩tau;
F)控制接口传递的变量为实数信号,用于阀门和/或调节器控制信号的输入,以控制阀门的开启关闭和/或调节器的调节状态;
3)将步骤1)的各组件仿真模型,利用步骤2)定义的相应接口相连,通过接口传递数据;
4)对步骤3)通过接口连接各组件后的液体火箭发动机系统,进行基于Modelica语言的液体火箭发动机动态仿真。
2.根据权利要求1所述的基于Modelica语言的液体火箭发动机动态仿真方法,其特征在于:
步骤A)中,所述流经纯气体介质的组件包括气瓶、气体管路、气体节流孔和/或气体阀门。
3.根据权利要求1所述的基于Modelica语言的液体火箭发动机动态仿真方法,其特征在于:
步骤B)中,所述流经单相液体介质的组件包括液体管路、液体节流孔、液体阀门、离心泵和/或液涡轮。
4.根据权利要求1所述的基于Modelica语言的液体火箭发动机动态仿真方法,其特征在于:
步骤C)中,所述流经气液混合介质的组件包括低温介质充填管路、头腔和/或低温充填冷却套。
5.根据权利要求1所述的基于Modelica语言的液体火箭发动机动态仿真方法,其特征在于:
步骤D)中,所述热力组件包括发生器、燃气导管、燃气涡轮和/或推力室。
6.根据权利要求1至5任一项所述的基于Modelica语言的液体火箭发动机动态仿真方法,其特征在于:
步骤B)中,所述单相液体介质为氧化剂和/或燃料。
7.根据权利要求6所述的基于Modelica语言的液体火箭发动机动态仿真方法,其特征在于:
步骤1)中,所述头腔为乳化吹除头腔。
CN202111295236.9A 2021-11-03 2021-11-03 一种基于Modelica语言的液体火箭发动机动态仿真方法 Pending CN114021253A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111295236.9A CN114021253A (zh) 2021-11-03 2021-11-03 一种基于Modelica语言的液体火箭发动机动态仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111295236.9A CN114021253A (zh) 2021-11-03 2021-11-03 一种基于Modelica语言的液体火箭发动机动态仿真方法

Publications (1)

Publication Number Publication Date
CN114021253A true CN114021253A (zh) 2022-02-08

Family

ID=80060472

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111295236.9A Pending CN114021253A (zh) 2021-11-03 2021-11-03 一种基于Modelica语言的液体火箭发动机动态仿真方法

Country Status (1)

Country Link
CN (1) CN114021253A (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114239438A (zh) * 2022-02-18 2022-03-25 中国汽车技术研究中心有限公司 一种氢气循环设备仿真方法和系统
CN115408875A (zh) * 2022-09-16 2022-11-29 苏州清动碳零信息科技有限公司 一种基于Modelica语言的溴化锂吸收式制冷系统动态仿真方法
CN116257942A (zh) * 2023-05-16 2023-06-13 东方空间技术(山东)有限公司 一种火箭仿真模型的确定方法及装置
CN117216899A (zh) * 2023-09-11 2023-12-12 中国人民解放军国防科技大学 一种基于irc方法的液氧煤油发动机预燃室建模方法
CN117216900A (zh) * 2023-09-11 2023-12-12 中国人民解放军国防科技大学 一种基于irc方法的液氧煤油发动机燃烧室建模方法
CN117875095A (zh) * 2024-03-13 2024-04-12 西安航天动力研究所 一种发动机频域仿真模型的建模方法、装置及电子设备

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114239438A (zh) * 2022-02-18 2022-03-25 中国汽车技术研究中心有限公司 一种氢气循环设备仿真方法和系统
CN114239438B (zh) * 2022-02-18 2022-06-17 中国汽车技术研究中心有限公司 一种氢气循环设备仿真方法和系统
CN115408875A (zh) * 2022-09-16 2022-11-29 苏州清动碳零信息科技有限公司 一种基于Modelica语言的溴化锂吸收式制冷系统动态仿真方法
CN115408875B (zh) * 2022-09-16 2024-02-02 苏州清动碳零信息科技有限公司 基于Modelica语言的溴化锂吸收式制冷系统动态仿真方法
CN116257942A (zh) * 2023-05-16 2023-06-13 东方空间技术(山东)有限公司 一种火箭仿真模型的确定方法及装置
CN116257942B (zh) * 2023-05-16 2023-08-15 东方空间技术(山东)有限公司 一种火箭仿真模型的确定方法及装置
CN117216899A (zh) * 2023-09-11 2023-12-12 中国人民解放军国防科技大学 一种基于irc方法的液氧煤油发动机预燃室建模方法
CN117216900A (zh) * 2023-09-11 2023-12-12 中国人民解放军国防科技大学 一种基于irc方法的液氧煤油发动机燃烧室建模方法
CN117216900B (zh) * 2023-09-11 2024-04-02 中国人民解放军国防科技大学 一种基于irc方法的液氧煤油发动机燃烧室建模方法
CN117216899B (zh) * 2023-09-11 2024-05-17 中国人民解放军国防科技大学 一种基于irc方法的液氧煤油发动机预燃室建模方法
CN117875095A (zh) * 2024-03-13 2024-04-12 西安航天动力研究所 一种发动机频域仿真模型的建模方法、装置及电子设备
CN117875095B (zh) * 2024-03-13 2024-06-21 西安航天动力研究所 一种发动机频域仿真模型的建模方法、装置及电子设备

Similar Documents

Publication Publication Date Title
CN114021253A (zh) 一种基于Modelica语言的液体火箭发动机动态仿真方法
Gravdahl et al. Compressor surge and rotating stall: modeling and control
US6748349B1 (en) Generalized fluid system simulation program
Whitmore et al. Closed-loop precision throttling of a hybrid rocket motor
Di Matteo et al. Start-up transient simulation of a liquid rocket engine
CA2966108A1 (en) System and method for stabilizing transcritical air-to-fuel heat exchange
Chesse et al. Performance simulation of sequentially turbocharged marine diesel engines with applications to compressor surge
CN112800562A (zh) 推力控制方法、装置、电子设备及存储介质
CN114021252A (zh) 一种基于Modelica语言的液体火箭发动机动态仿真模型库架构
CN115422772A (zh) 基于布雷顿循环容积设备模型的流量-压力耦合瞬态计算方法
McCarthy et al. A multi-domain component based modeling toolset for dynamic integrated power and thermal system modeling
Rockwell et al. Simulated scramjet shock train control using an all-coefficient adaptive control approach
Yamanishi et al. Transient analysis of the LE-7A rocket engine using the rocket engine dynamic simulator (REDS)
Zhang et al. Control law synthetizing for an innovative indirect precooled airbreathing engine under off-design operation conditions
Su et al. Dynamic characteristics of lox/kerosene variable thrust liquid rocket engine test system based on general modular simulation method
Jia et al. Parametric analysis of variable stator vane system in gas turbines based on cosimulation of its refined model and system dynamic performance model
Moral et al. Espss model of a simplified combined-cycle engine for supersonic cruise
Di Matteo et al. Steady state library for liquid rocket engine cycle design
Murray et al. Design and analysis of a high speed, high pressure peroxide/RP-1 turbopump
Leonardi et al. A methodology to study the possible occurrence of chugging in liquid rocket engines during transient start-up
Wang et al. Modeling of pogo vibration in rockets with staged combustion cycle engine
Lungu et al. Aircraft double-spool single-jet engine's model
Sheikhbahaei et al. Optimal tuner selection using Kalman filter for a real-time modular gas turbine model
Park et al. Monte-Carlo Simulation for Analyzing the Performance Variation of a Liquid Rocket Engine Using Gas-Generator Cycle
Krishnaswami et al. Fuel-efficient operation of compressor stations using simulation-based optimization

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