一种火箭导流槽用耐火材料烧蚀量的确定方法和装置

技术领域

本发明涉及飞行器设计技术领域，尤其涉及一种火箭导流槽用耐火材料烧蚀量的确定方法和装置。

背景技术

航天器点火后，发动机喷出的燃气流具有强大的冲击力与侵蚀性。导流槽作为火箭发射平台的关键基础性设施，不仅需承受发射时产生的高温、高压、高速燃气流的冲刷，还需将高温有毒燃气分流，快速导向远离发射塔的空间，为航天器周围的建筑设施、发动机及相关设备提供安全保障。为了提升导流槽的耐高温性能，通常在导流槽表面浇筑专用耐火混凝土，以保证试验发射安全进行。然而，随着航天技术的飞速发展，对航天器运载能力提出了新的更高的要求，导流槽的烧蚀问题日益严重，迫切需要新技术新方法对其烧蚀行为进行研究与预测，为航天器发射台的搭建提供技术支持。

目前，导流槽耐火混凝土的烧蚀性能一般通过试片级试验预测，该方法只能对材料是否符合国军标规范进行验证，无法得出在实际发射和试车过程中导流槽的烧蚀情况。而且，试验需采用喷管出口燃气温度、速度、压力满足要求的煤油发动机，且需提前制备多个立方体试样，试验成本较高，周期较长。

发明内容

有鉴于此，本发明实施例的目的在于提供一种火箭导流槽用耐火材料烧蚀量的确定方法和装置，以解决现有技术中通过试片级试验预测导流槽耐火混凝土的烧蚀性能，无法得出在实际发射和试车过程中导流槽的烧蚀情况的技术问题。

为达上述目的，第一方面，本发明实施例提供了一种火箭导流槽用耐火材料烧蚀量的确定方法，所述确定方法包括：

获取耐火材料的烧蚀试验数据，对所述烧蚀试验数据进行拟合获得烧蚀量计算公式；

根据所述烧蚀量计算公式中的剥蚀因子确定耐火材料的线烧蚀速率；

根据所述线烧蚀速率计算所述耐火材料在不同时间和不同位置处的烧蚀后退量，根据所述烧蚀后退量创建导流槽模型并绘制不同时刻的烧蚀型面；

根据所述导流槽模型和所述烧蚀型面模拟所述耐火材料在烧蚀过程中的外形变化情况，获得实际烧蚀过程；

根据所述实际烧蚀过程计算所述耐火材料在烧蚀过程中的温度场和应力场的耦合关系；

根据所述温度场和应力场的耦合关系通过所述烧蚀后退量的计算公式计算所述耐火材料的烧蚀量。

在一些可能的实施方式中，所述的根据所述线烧蚀速率计算所述耐火材料在不同时间和不同位置处的烧蚀后退量，具体包括：

通过所述线烧蚀速率在所述不同时间和不同位置处对应的时间的积分计算所述所述烧蚀后退量，所述烧蚀后退量的计算公式为：δ＝f(T,v,P,t)，其中，

δ为烧蚀后退量，T为温度，v为速度，P为压力，t为燃气冲刷时间。

在一些可能的实施方式中，所述的根据所述导流槽模型和所述烧蚀型面模拟所述耐火材料在烧蚀过程中的外形变化情况，获得实际烧蚀过程，具体包括：

将所述导流槽模型导入有限元仿真软件，根据所述烧蚀型面分段模拟所述耐火材料在烧蚀过程中的外形变化情况，获得实际烧蚀过程。

在一些可能的实施方式中，所述的根据所述实际烧蚀过程计算所述耐火材料在烧蚀过程中的温度场和应力场的耦合关系，具体包括：

根据所述实际烧蚀过程创建有限元分析模型，并对所述有限元分析模型进行划分网格；

采用生死单元法将所述有限元分析模型中达到烧蚀条件的网格进行逐层剥离；

获取所述耐火材料的物性参数，根据所述物性参数计算所述逐层剥离的过程中不同时刻的温度场和应力场的耦合关系。

在一些可能的实施方式中，所述的采用生死单元法将所述有限元分析模型中达到烧蚀条件的网格进行逐层剥离，具体包括：

根据所述耐火材料的临界温度设置生死单元，根据所述生死单元将所述有限元分析模型中达到烧蚀条件的网格进行逐层剥离。

在一些可能的实施方式中，所述的对所述有限元分析模型进行划分网格，具体包括：

通过布置节点数量、节点分布、网格尺寸和网格类型对所述有限元分析模型进行划分网格。

第二方面，本发明实施例提供了一种火箭导流槽用耐火材料烧蚀量的确定装置，所述确定装置包括：

获取与拟合模块，用于获取耐火材料的烧蚀试验数据，对所述烧蚀试验数据进行拟合获得烧蚀量计算公式；

第一计算模块，用于根据所述烧蚀量计算公式中的剥蚀因子计算耐火材料的线烧蚀速率，根据所述线烧蚀速率计算所述耐火材料在不同时间和不同位置处的烧蚀后退量；

创建与绘制模块，用于根据所述烧蚀后退量创建导流槽模型并绘制不同时刻的烧蚀型面；

模拟模块，用于根据所述导流槽模型和所述烧蚀型面模拟所述耐火材料在烧蚀过程中的外形变化情况，获得实际烧蚀过程；

第二计算模块，用于根据所述实际烧蚀过程计算所述耐火材料在烧蚀过程中的温度场和应力场的耦合关系；

第三计算模块，用于根据所述温度场和应力场的耦合关系通过所述烧蚀后退量的计算公式计算所述耐火材料的烧蚀量。

在一些可能的实施方式中，所述第二计算模块具体包括：

创建与划分子模块，用于根据所述实际烧蚀过程创建有限元分析模型，并对所述有限元分析模型进行划分网格；

剥离子模块，用于采用生死单元法将所述有限元分析模型中达到烧蚀条件的网格进行逐层剥离；

获取与计算子模块，用于获取所述耐火材料的物性参数，根据所述物性参数计算所述逐层剥离的过程中不同时刻的温度场和应力场的耦合关系。

第三方面，本发明实施例提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现一种火箭导流槽用耐火材料烧蚀量的确定方法。

第四方面，本发明实施例提供了一种计算机设备，其包括：

处理器；

用于存储所述处理器可执行指令的存储器；

其中，所述处理器被配置为执行所述指令，以实现一种火箭导流槽用耐火材料烧蚀量的确定方法。

上述技术方案的有益技术效果在于：

本发明实施例提供的一种火箭导流槽用耐火材料烧蚀量的确定方法和装置，该控制方法包括：获取烧蚀试验数据，对烧蚀试验数据进行拟合获得烧蚀量计算公式；根据烧蚀量计算公式中的剥蚀因子确定耐火材料的线烧蚀速率；根据线烧蚀速率计算耐火材料在不同时间和不同位置处的烧蚀后退量，根据烧蚀后退量创建导流槽模型并绘制不同时刻的烧蚀型面；根据导流槽模型和烧蚀型面模拟耐火材料在烧蚀过程中的外形变化情况，获得实际烧蚀过程；根据实际烧蚀过程计算耐火材料在烧蚀过程中的温度场和应力场的耦合关系；根据温度场和应力场的耦合关系通过烧蚀量计算公式计算耐火材料的烧蚀量。本发明实施例通过拟合算法获得导流槽耐火混凝土在高温高速燃气流条件下的烧蚀量计算公式和线烧蚀速率，并在此基础上基于有限元法与生死单元法，对实际导流槽模型的烧蚀形貌进行预测，并给出各个部位的烧蚀量大小，在材料烧蚀剥离过程中获得温度场/应力场的耦合关系，根据温度场和应力场的耦合关系可以计算出不同时刻不同位置处的烧蚀量，可为未来大推力运载火箭导流槽与发射塔的工程设计提供烧蚀预测手段，为航天器发射任务的顺利进行提供技术支撑。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明实施例的一种火箭导流槽用耐火材料烧蚀量的确定方法的流程图；

图2是本发明实施例的一种耐火混凝土高温烧蚀过程中的形成的物质图；

图3是本发明实施例的一种高温烧蚀过程中耐火混凝土表面的能量平衡关系图；

图4是本发明实施例的一种火箭导流槽用耐火材料烧蚀量的确定装置的结构框图；

图5是本发明实施例的一种第二计算模块的结构框图；

图6是本发明实施例的一种计算机可读存储介质的功能框图；

图7是本发明实施例的一种计算机设备的功能框图。

具体实施方式

下面将详细描述本发明的各个方面的特征和示例性实施例。在下面的详细描述中，提出了许多具体细节，以便提供对本发明的全面理解。但是，对于本领域技术人员来说很明显的是，本发明可以在不需要这些具体细节中的一些细节的情况下实施。下面对实施例的描述仅仅是为了通过示出本发明的示例来提供对本发明的更好的理解。在附图和下面的描述中，至少部分的公知结构和技术没有被示出，以便避免对本发明造成不必要的模糊；并且，为了清晰，可能夸大了部分结构的尺寸。此外，下文中所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施例中。

实施例一

图1是本发明一种火箭导流槽用耐火材料烧蚀量的确定方法的流程图，如图1所示，该确定方法包括如下步骤：

步骤S1，获取耐火材料的烧蚀试验数据，对烧蚀试验数据进行拟合获得烧蚀量计算公式。

本实施例中，耐火材料可以包括耐火混凝土，通过高温射流试验平台对耐火混凝土试件进行烧蚀获得耐火混凝土的烧蚀试验数据，该烧蚀试验数据包括烧蚀前质量、烧蚀后质量、烧蚀后退量和温度等，本实施例中，试验时会在耐火混凝土试验件上安装温度测量传感器，可以在试验过程中测量每一时刻的温度数值，可以在试验过程中通过测量获得，根据这些烧蚀试验数据通过拟合编程的方式获得该耐火混凝土的烧蚀量计算公式，该烧蚀量计算公式中包含耐火混凝土、温度场与应力场三者之间的耦合计算关系。

步骤S2，根据烧蚀量计算公式中的剥蚀因子确定耐火材料的线烧蚀速率。

本实施例中，航天器点火后，发动机喷出的燃气流具有强大的冲击力与侵蚀性，产生的高温、高压、高速燃气流的冲刷导流器上的耐火混凝土时，耐火混凝土会减少，减少的这部分也可以叫做烧蚀退后量，烧蚀退后量中包括被高温烧掉的和被气流冲刷掉的，被气流给冲刷掉的这一部分可以称作剥蚀掉的，本实施例中的剥蚀因子为剥蚀掉的耐火材料在整个烧蚀退后量里面的质量的占比。根据能量守恒、动量守恒和质量守恒，通过剥蚀因子计算耐火混凝土烧蚀过程中的线烧蚀速率。

步骤S3，根据线烧蚀速率计算耐火材料在不同时间和不同位置处的烧蚀后退量，根据烧蚀后退量创建导流槽模型并绘制不同时刻的烧蚀型面。

本实施例中，通过线烧蚀速率在不同时间和不同位置处对应的时间的积分计算烧蚀后退量，烧蚀后退量的计算公式为：δ＝f(T,v,P,t)，其中，δ为烧蚀后退量，T为温度，v为速度，P为压力，t为燃气冲刷时间；所以，可以根据线烧蚀速率实时地计算不同时刻不同位置的耐火混凝土的烧蚀后退量，根据烧蚀后退量在CAD软件中创建导流槽模型并绘制烧蚀不同时刻烧蚀型面。

步骤S4，根据导流槽模型和烧蚀型面模拟耐火材料在烧蚀过程中的外形变化情况，获得实际烧蚀过程。

具体的，将导流槽模型导入有限元计算程序中，通过有限元软件对烧蚀型面进行仿真，通过模拟耐火混凝土在烧蚀过程中烧蚀量、温度场和应力场三者之间的变化关系，获得实际烧蚀过程，即耐火混凝土在实际烧蚀过程中外形变化的模型。

步骤S5，根据实际烧蚀过程计算耐火材料在烧蚀过程中的温度场和应力场的耦合关系。

具体的，根据耐火混凝土在实际烧蚀过程中外形变化的模型，计算出温度场和应力场的耦合关系，本实施例中，在有限元软件中可选择仿真类型为应力场与温度场耦合分析。

步骤S6，根据温度场和应力场的耦合关系，通过烧蚀后退量的计算公式计算耐火材料的烧蚀量。

具体的，在步骤S1中获得的烧蚀量计算公式中包含了烧蚀量、温度场和应力场三者之间的变化关系，因此，在通过步骤S5中计算出温度场和应力场的耦合关系后，再通过烧蚀后退量的计算公式δ＝f(T,v,P,t)计算耐火材料的烧蚀量，其中，δ为烧蚀后退量，T为温度，v为速度，P为压力，t为燃气冲刷时间。

本实施例中，步骤S1至步骤S4中通过仿真获得了实际烧蚀过程中外形的变化情况，即得到了实际烧蚀过程的模型，而在步骤S5和步骤S6中通过这个模型就可以计算出实际烧蚀过程中的温度场和应力场的耦合关系，再通过烧蚀量计算公式即可实时地计算出耐火混凝土在不同时刻的烧蚀量。

本发明实施例通过拟合的算法获得导流槽耐火混凝土在高温高速燃气流条件下的烧蚀量计算公式和线烧蚀速率，并在此基础上基于有限元法，对实际导流槽模型的烧蚀形貌进行仿真，并给出各个部位的烧蚀量大小，在材料烧蚀剥离过程中获得温度场和应力场的耦合关系，根据温度场/应力场的耦合关系可以计算出不同时刻不同位置处的烧蚀量，可为未来大推力运载火箭导流槽与发射塔的工程设计提供烧蚀预测手段，为航天器发射任务的顺利进行提供技术支撑。

在一些实施例中，步骤S4中根据导流槽模型和烧蚀型面模拟耐火材料在烧蚀过程中的外形变化情况，获得实际烧蚀过程，具体包括：

将导流槽模型导入有限元仿真软件，根据烧蚀型面依次模拟耐火材料在烧蚀过程中的外形变化情况，获得实际烧蚀过程。本发明实施例中，采用依次模拟的方法对烧蚀过程中外形的变化进行仿真获得实际烧蚀过程，即一步接着一步地模拟，极大地减少了运算量，提高数值模拟速度。

在一些实施例中，步骤S5中根据实际烧蚀过程计算耐火材料在烧蚀过程中的温度场和应力场的耦合关系，具体包括：

步骤S51，根据实际烧蚀过程创建有限元分析模型，并对有限元分析模型进行划分网格；

本实施例中，根据实际烧蚀过程创建有限元分析模型之后，为了便于后续的计算，可以通过布置节点数量、节点分布(节点是连接相邻单元的网格点)、网格尺寸和网格类型对有限元分析模型进行划分网格。

步骤S52，采用生死单元法将有限元分析模型中达到烧蚀条件的网格进行逐层剥离。

图2是本发明实施例的一种耐火混凝土高温烧蚀过程中的形成的物质图，图3是本发明实施例一种高温烧蚀过程中耐火混凝土表面的能量平衡关系图，如图2和图3所示，耐火混凝土物质层(D)受高温气流烧蚀时，首先混凝土中的水合化合物分解发生脱水反应，形成分解层(C)，随着温度的进一步升高，脱水后的产物之间以及产物与骨料之间进行了一系列的化学反应，发生转化、烧结现象，形成烧结层(B)；之后一部分脱水后的产物达到了熔点变成液相，形成液态层(A)，并附着于耐火混凝土的烧结层(B)表面(耐火混凝土)，烧蚀过程中，表面烧蚀后退产生一定的退后量(E)，并由此建立烧蚀表面能量平衡方程为：

q

式中，q

本实施例中，步骤S52采用生死单元法将有限元分析模型中达到烧蚀条件的网格进行逐层剥离，即耐火材料(如耐火混凝土)达到一定的温度，然后它就会被剥蚀，会被气流给吹走，剥蚀和吹走后就会烧蚀往后退，网格生成了之后，就会有很多层，可以根据这个温度来定义生死单元，选择要剥蚀的网格，根据生死单元来控制剥离过程，即一层一层的剥离。

步骤S53，获取耐火材料的物性参数，根据物性参数计算逐层剥离的过程中不同时刻的温度场和应力场的耦合关系。

耐火材料的物性参数是表征该材料固有属性的参数，物性参数一般包括该耐火材料的热导率和比热等，用于计算温度场和应力场之间的耦合关系。本实施例中，可以通过力热耦合的方法根据热流等边界条件，在材料烧蚀剥离过程中计算得出各个时刻的温度场和应力耦合计算结果，其中，热流是通过热环境计算获得，是一个输入条件。

在一些实施例中，步骤S52采用生死单元法将有限元分析模型中达到烧蚀条件的网格进行逐层剥离，具体包括：根据耐火材料的临界温度设置生死单元，根据生死单元将有限元分析模型中达到烧蚀条件的网格进行逐层剥离。

本实施例中，每一种耐火材料都会有一个临界温度，生死单元考虑要不要逐层剥离，生死单元法的核心思想是通过单元的“生”(激活)或者“死”(抑制)来模拟一些特殊的物理工况，即通过生死单元来控制这个剥离过程。

本发明实施例通过拟合算法获得导流槽耐火混凝土在高温高速燃气流条件下的烧蚀量计算公式和线烧蚀速率，并在此基础上基于有限元法与生死单元法，对实际导流槽模型的烧蚀形貌进行预测，并给出各个部位的烧蚀量大小，在材料烧蚀剥离过程中获得温度场/应力场的耦合关系，根据温度场/应力场的耦合关系可以计算出不同时刻不同位置处的烧蚀量，可为未来大推力运载火箭导流槽与发射塔的工程设计提供烧蚀预测手段，为航天器发射任务的顺利进行提供技术支撑。

图4是本发明实施例的一种火箭导流槽用耐火材料烧蚀量的确定装置的结构框图，如图4所示，该确定装置100包括：

获取与拟合模块101，用于获取耐火材料的烧蚀试验数据，对烧蚀试验数据进行拟合获得烧蚀量计算公式。

本实施例中，耐火材料可以包括耐火混凝土，通过高温射流试验平台对耐火混凝土试件进行烧蚀获得耐火混凝土的烧蚀试验数据，该烧蚀试验数据包括烧蚀前质量、烧蚀后质量、烧蚀后退量和温度等，本实施例中，试验时会在耐火混凝土试验件上安装温度测量传感器，可以在试验过程中测量每一时刻的温度数值，可以在试验过程中通过测量获得，根据这些烧蚀试验数据通过拟合编程的方式获得该耐火混凝土的烧蚀量计算公式，该烧蚀量计算公式中包含耐火混凝土、温度场与应力场三者之间的耦合计算关系。

第一计算模块102，用于根据烧蚀量计算公式中的剥蚀因子计算耐火材料的线烧蚀速率，根据线烧蚀速率计算耐火材料在不同时间和不同位置处的烧蚀后退量。

本实施例中，航天器点火后，发动机喷出的燃气流具有强大的冲击力与侵蚀性，产生的高温、高压、高速燃气流的冲刷导流器上的耐火混凝土时，耐火混凝土会减少，减少的这部分也可以叫做烧蚀退后量，烧蚀退后量中包括被高温烧掉的和被气流冲刷掉的，被气流给冲刷掉的这一部分可以称作剥蚀掉的，本实施例中的剥蚀因子为剥蚀掉的耐火材料在整个烧蚀退后量里面的质量的占比。根据能量守恒、动量守恒和质量守恒，通过剥蚀因子计算耐火混凝土烧蚀过程中的线烧蚀速率。

创建与绘制模块103，用于根据烧蚀后退量创建导流槽模型并绘制不同时刻的烧蚀型面。

本实施例中，创建与绘制模块103通过线烧蚀速率在不同时间和不同位置处对应的时间的积分计算烧蚀后退量，烧蚀后退量的计算公式为：δ＝f(T,v,P,t)，其中，δ为烧蚀后退量，T为温度，v为速度，P为压力，t为燃气冲刷时间；所以，可以根据线烧蚀速率实时地计算不同时刻不同位置的耐火混凝土的烧蚀后退量，根据烧蚀后退量在CAD软件中创建导流槽模型并绘制烧蚀不同时刻烧蚀型面。

模拟模块104，用于根据导流槽模型和烧蚀型面模拟耐火材料在烧蚀过程中的外形变化情况，获得实际烧蚀过程。

具体的，模拟模块104将导流槽模型导入有限元计算程序中，通过有限元软件对烧蚀型面进行仿真，通过模拟耐火混凝土在烧蚀过程中烧蚀量、温度场和应力场三者之间的变化关系，获得实际烧蚀过程，即耐火混凝土在实际烧蚀过程中外形变化的模型。

第二计算模块105，用于根据实际烧蚀过程计算耐火材料在烧蚀过程中的温度场和应力场的耦合关系。

具体的，第二计算模块105根据耐火混凝土在实际烧蚀过程中外形变化的模型，计算出温度场和应力场的耦合关系，本实施例中，在有限元软件中可选择仿真类型为应力场与温度场耦合分析。

第三计算模块106，用于根据温度场和应力场的耦合关系，通过烧蚀后退量的计算公式计算耐火材料的烧蚀量。

具体的，在获取与拟合模块101中获得的烧蚀量计算公式中包含了烧蚀量、温度场和应力场三者之间的变化关系，因此，在通过第二计算模块105计算出温度场和应力场的耦合关系后，再通过计算获得烧蚀量。

具体的，在获取与拟合模块101中获得的烧蚀量计算公式中包含了烧蚀量、温度场和应力场三者之间的变化关系，因此，在通过二计算模块105中计算出温度场和应力场的耦合关系后，最后通过烧蚀后退量的计算公式δ＝f(T,v,P,t)计算耐火材料的烧蚀量，其中，δ为烧蚀后退量，T为温度，v为速度，P为压力，t为燃气冲刷时间。

本实施例中，先通过仿真获得了实际烧蚀过程中外形的变化情况，即得到了实际烧蚀过程的模型，再通过这个模型就可以计算出实际烧蚀过程中的温度场和应力场的耦合关系，再通过烧蚀量计算公式即可实时地计算出耐火混凝土的烧蚀量。

图5是本发明实施例的一种第二计算模块的结构框图，如图5所示，该第二计算模块105具体包括：

创建与划分子模块1051，用于根据实际烧蚀过程创建有限元分析模型，并对有限元分析模型进行划分网格。

本实施例中，根据实际烧蚀过程创建有限元分析模型之后，为了便于后续的计算，可以通过布置节点数量、节点分布、网格尺寸和网格类型对有限元分析模型进行划分网格。

剥离子模块1052，用于采用生死单元法将有限元分析模型中达到烧蚀条件的网格进行逐层剥离。

本实施例中，采用生死单元法将有限元分析模型中达到烧蚀条件的网格进行逐层剥离，即耐火材料(如耐火混凝土)达到一定的温度，然后它就会被剥蚀，会被气流给吹走，剥蚀和吹走后就会烧蚀往后退，网格生成了之后，就会有很多层，可以根据这个温度来定义生死单元，选择要剥蚀的网格，根据生死单元来控制剥离过程，即一层一层的剥离。

获取与计算子模块1053，用于获取耐火材料的物性参数，根据物性参数计算逐层剥离的过程中不同时刻的温度场和应力场的耦合关系。

耐火材料的物性参数是表征该材料固有属性的参数，物性参数一般包括该耐火材料的热导率和比热等，用于计算温度场和应力场之间的耦合关系。本实施例中，可以通过力热耦合的方法根据热流等边界条件，在材料烧蚀剥离过程中计算得出各个时刻的温度场/应力耦合计算结果，其中热流是通过热环境计算获得，是一个输入条件。

本发明实施例通过拟合算法获得导流槽耐火混凝土在高温高速燃气流条件下的烧蚀量计算公式和线烧蚀速率，并在此基础上基于有限元法与生死单元法，对实际导流槽模型的烧蚀形貌进行预测，并给出各个部位的烧蚀量大小，在材料烧蚀剥离过程中获得温度场/应力场的耦合关系，根据温度场/应力场的耦合关系可以计算出不同时刻不同位置处的烧蚀量，可为未来大推力运载火箭导流槽与发射塔的工程设计提供烧蚀预测手段，为航天器发射任务的顺利进行提供技术支撑。

实施例三

图6是本发明实施例的一种计算机可读存储介质的功能框图。如图6所示，本发明实施例还提供了一种计算机可读存储介质200，计算机可读存储介质200内存储有计算机程序210，该计算机程序210被处理器执行时，实现：

获取烧蚀试验数据，对烧蚀试验数据进行拟合获得烧蚀量计算公式；

根据烧蚀量计算公式中的剥蚀因子确定耐火材料的线烧蚀速率；

根据线烧蚀速率计算耐火材料在不同时间和不同位置处的烧蚀后退量，根据烧蚀后退量创建导流槽模型并绘制不同时刻的烧蚀型面；

根据导流槽模型和烧蚀型面模拟耐火材料在烧蚀过程中的外形变化情况，获得实际烧蚀过程；

根据实际烧蚀过程计算耐火材料在烧蚀过程中的温度场和应力场的耦合关系；

根据温度场和应力场的耦合关系通过烧蚀后退量的计算公式计算耐火材料的烧蚀量。

所述集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、电载波信号、电信信号以及软件分发介质等。当然，还有其他方式的可读存储介质，例如量子存储器、石墨烯存储器等等。需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括电载波信号和电信信号。

实施例四

图7是本发明实施例的一种计算机设备的功能框图。本发明实施例还提供一种计算机设备，请参考图7，在硬件层面，该计算机设备包括处理器，可选地还包括内部总线、网络接口、存储器。其中，存储器可能包含内存，例如高速随机存取存储器(Random-AccessMemory，RAM)，也可能还包括非易失性存储器(non-volatile memory)，例如至少1个磁盘存储器等。当然，该电子设备还可能包括其他业务所需要的硬件。

处理器、网络接口和存储器可以通过内部总线相互连接，该内部总线可以是工业标准体系结构ISA总线、外设部件互连标准PCI总线或扩展工业标准结构EISA总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示，图7中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。

存储器，用于存放程序。具体地，程序可以包括程序代码，所述程序代码包括计算机操作指令。存储器可以包括内存和非易失性存储器，并向处理器提供指令和数据。处理器从非易失性存储器中读取对应的计算机程序到内存中然后运行，在逻辑层面上形成一种基于三摆伺服发动机火箭姿态的控制方法。处理器，执行存储器所存放的程序，并具体用于执行图1所示实施例揭示的一种基于三摆伺服发动机火箭姿态的控制方法。

上述如图1所示实施例揭示的一种基于三摆伺服发动机火箭姿态的控制方法可以应用于处理器中，或者由处理器实现。处理器可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器可以是通用处理器，包括中央处理器(CentralProcessingUnit，CPU)、网络处理器(Network Processor，NP)等；还可以是数字信号处理器(DigitalSignal Processor，DSP)、专用集成电路(Application Specific IntegratedCircuit，ASIC)、现场可编程门阵列(Field－Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。

当然，除了软件实现方式之外，本发明的电子设备并不排除其他实现方式，比如逻辑器件抑或软硬件结合的方式等等，也就是说以下处理流程的执行主体并不限定于各个逻辑单元，也可以是硬件或逻辑器件。上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。

虽然本发明提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。

本发明是参照根据本发明实施例的方法、设备、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

本说明书中的各个实施例均采用相关的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于装置、电子设备及可读存储介质实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

以上所述仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等，均包含在本发明的保护范围内。