CN109033515B - 一种微型燃气涡轮发动机起动过程建模方法 - Google Patents

Abstract

本发明属于发动机建模技术领域，提供了一种微型燃气涡轮发动机起动过程建模方法，步骤如下：微型燃气涡轮发动机转速的建模；微型燃气涡轮发动机的性能参数与转速的关系建模；误差分析。现有微型燃气涡轮发动机建模方法多是利用纯机理或机器学习的方法，难以准确刻画微型燃气涡轮发动机的起动过程，而机器学习需要大量的试验数据。在工程实践中，目前普遍采用的纯机理操作繁复，效率低下，建模准确度不高。本发明提供了一种基于机理和辨识方法相结合的微型燃气涡轮发动机起动过程建模方法，弥补了现有技术的不足。本发明操作简单，准确率高，并可实现整个微型燃气涡轮发动机的建模。该方法具有一定的拓展性，可推广到其它领域。

Description

一种微型燃气涡轮发动机起动过程建模方法

技术领域

本发明属于发动机建模技术领域，具体涉及一种微型燃气涡轮发动机起动过程建模方法。

背景技术

微型燃气涡轮发动机是一个复杂的热力学系统，由压气机、燃烧室和涡轮组成。在对微型燃气涡轮发动机起动的研究过程中，如何建立微型燃气涡轮发动机慢车转速以下的数学模型并进行起动性能的数值计算，能为微型燃气涡轮发动机起动过程的控制提供参考。因此，对微型燃气涡轮发动机的起动过程进行建模是非常有价值的。

在微型燃气涡轮发动机建模技术方面，国内外公开发表的文献较少，亦未有相关专利涉及此方面的内容。在已有文献中，多是利用支持向量机等方法对发动机起动过程进行建模，但是支持向量机需要大量的发动机起动过程数据，通过对输入输出数据的训练得到发动机的起动模型，然而对于新研发的微型燃气涡轮发动机而言，一般通过燃气涡轮发动机的设计参数得到起动过程的燃油计划，为了保证燃气涡轮发动机能成功起动，通过人工手动的方法完成微型燃气涡轮发动机的起动过程。通过人工起动得到的实验数据建立燃气涡轮发动机的起动模型，然后再通过起动模型进一步优化发动机起动过程的燃油计划，存在一定的局限性。

对于微型燃气涡轮发动机而言，当前的起动过程建模技术已经渐渐满足不了工程的实际需要。因此，探索出一种有效且适用于工程实践的高效建模方法成为一个亟需解决的问题。微型燃气涡轮发动机起动建模技术也因此具备了广阔的研究和应用前景。

发明内容

本发明是为了解决现有技术存在的燃气涡轮发动机起动过程建模数据要求量大，效率低下，建模精度不高的问题，而提出的一种微型燃气涡轮发动机起动过程建模方法。

一种微型燃气涡轮发动机起动过程建模方法，包括以下步骤：

第一步，微型燃气涡轮发动机转速的建模；

第二步，微型燃气涡轮发动机的性能参数与转速的关系建模；

第三步，误差分析。

现具体阐述如下：

第一步，微型燃气涡轮发动机转速的建模；

步骤1：分析微型燃气涡轮发动机起动过程机理，通过发动机转子效应建模方法，具体是通过微型燃气涡轮发动机起动电机的电流值和微型燃气涡轮发动机的转速值，建立燃气涡轮发动机剩余转矩模型，在此基础上计算得到稳态燃料量，从而通过线性插值的方法得到微型燃气涡轮发动机燃料量和微型燃气涡轮发动机转速之间的关系，运用多项式拟合的方法得到微型燃气涡轮发动机燃料量和转速之间的表达式；

步骤2：根据微型燃气涡轮发动机起动过程的燃料量和发动机转速的试验数据，并与步骤1所建的稳态燃料量和微型燃气涡轮发动机转速的关系式进行对比，得到起动过程的燃料量与稳态燃料量之间的差值；

步骤3：步骤2得到的燃料量差值能使微型燃气涡轮发动机产生加速效应，由于燃料量差值与微型燃气涡轮发动机转速之间的关系是高度非线性的，利用分段线性化的方法去拟合燃料量差值与微型燃气涡轮发动机转速之间的关系；

步骤4：燃料量差值和起动机电流共同作用，使微型燃气涡轮发动机产生加速效应，通过对微型燃气涡轮发动机的加速度进行积分，从而得到微型燃气涡轮发动机的速度；

第二步，微型燃气涡轮发动机的性能参数与转速的关系建模；

步骤5：由于微型燃气涡轮发动机的转速与燃气性能参数之间的关系十分密切，采用多项式拟合的方法得到微型燃气涡轮发动机性能参数(压气机出口压力、压气机出口温度、燃烧室出口温度、燃烧室出口压力、涡轮后出口温度和涡轮后出口压力等)与微型燃气涡轮发动机转速之间的关系；

步骤6：由于通过生成的多项式计算燃气涡轮发动机性能参数会减慢模型的速度，同时也可能引起数值计算不稳定的现象，所以将多项式进行离散化，生成一维插值表的形式，有利于提高起动模型的实时性。

第三步，误差分析。

所建立的微型燃气涡轮发动机起动模型的各个性能参数与试验数据的误差均在5％以内，各个性能参数此处指压气机出口压力、压气机出口温度、燃烧室出口压力、涡轮出口压力、涡轮出口温度、微型燃气涡轮发动机转速。

本发明的有益效果：

现有微型燃气涡轮发动机建模方法多是利用纯机理或机器学习的方法，难以准确刻画微型燃气涡轮发动机的起动过程，而机器学习需要大量的试验数据。在工程实践中，目前普遍采用的纯机理操作繁复，效率低下，建模准确度不高。本发明提供了一种基于机理和辨识方法相结合的微型燃气涡轮发动机起动过程建模方法，弥补了现有技术的不足。本发明操作简单，准确率高，并可实现整个微型燃气涡轮发动机的建模。该方法具有一定的拓展性，可推广到其它领域。

附图说明

图1为微型燃气涡轮发动机建模流程图。

图2为剩余转矩与燃料量的关系图。

图3为燃料量差值与剩余扭矩关系图。

图4为起动过程建模误差对比图，其中，(a)压气机出口温度误差；(b)压气机出口压力误差；(c)燃烧室出口压力误差；(d)涡轮出口压力误差；(e)涡轮出口温度误差；(f)燃机转速误差。

具体实施方式

为了使本发明的目的、技术及优点更加清楚明白，以下结合附图及实例，对本发明进行进一步详细说明。

一种微型燃气涡轮发动机起动过程建模方法，包括以下步骤：

步骤1：微型燃气涡轮发动机的起动过程分为3个阶段，第一阶段，发动机转速由零到涡轮开始产生功率的转速，发动机完全由起动机带动加速；第二阶段，由涡轮开始产生功率的转速到起动机脱开；第三阶段，由起动机脱开转速至慢车转速；

按照发动机起动过程的各个阶段，对其进行建模；

第一阶段，发动机转速由零到涡轮开始产生功率的转速(即从零转速到点火转速)，发动机完全由起动机带动加速，满足发动机转子运动方程，即满足公式(1)。

式中M_st＝K_st×I_st，其中I_st表示起动机电流值，K_st为起动机转矩常数，针对不同类型的起动机，其值不同。此处所指的起动机为电起动机。ω为微型燃气涡轮发动机的角速度，J为归算到微型燃气涡轮发动机轴上的转动惯量，即等效转动惯量，针对不同的微型燃气涡轮发动机的结构，利用机械能守恒进行归算。具体的归算公式参照(2)

式中，J₁为微型燃气涡轮发动机轴带动的负载1的转动惯量，ω₁为微型燃气涡轮发动机轴带动的负载1的角速度；J₂为微型燃气涡轮发动机轴带动的负载2的转动惯量，ω₂为微型燃气涡轮发动机轴带动的负载2的角速度；J_n为微型燃气涡轮发动机轴带动的负载n的转动惯量，ω_n为微型燃气涡轮发动机轴带动的负载n的角速度。

第二阶段，发动机转速由点火转速到起动机脱开转速，发动机在起动机和涡轮的共同带动下进行加速，满足发动机转子运动方程，即满足公式(3)。

式中，M_gas就是点火后微型燃气涡轮发动机的剩余扭矩(M_gas＝M_T-M_C-M_f)，M_T代表涡轮产生的转矩，M_C代表压气机消耗的转矩，一般与压气机转速的平方成正比，M_f表示摩擦引起的转矩消耗，一般与转速的平方成正比；M_st是电起动机提供的转矩，与(1)式中相同。J与公式(1)计算过程相同，ω公式(1)相同。

第三阶段，发动机在剩余转速的作用下由发动机在起动机脱开转速加速慢车转速的过程，满足公式(4)。

式中，J与公式(1)相同，ω公式(1)相同。M_gas与公式(3)计算过程相同。

步骤2：第一阶段微型燃气涡轮发动机转速的建模

首先，通过计算得到发动机等效转动惯量(具体参见公式(2))；利用转子运动方程式，结合起动机电流值，根据M_st＝K_st×I_st计算得到起动机加在微型燃气涡轮发动机轴上的转矩值，利用一阶龙格-库塔方法(欧拉方法)对其进行积分，得到该模型的微型燃气涡轮发动机转速，保证该微型燃气涡轮发动机的转速和试验数据中的燃机转速误差在5％以内；

步骤3：第二阶段微型燃气涡轮发动机转速的建模

步骤3.1：通过计算得到发动机等效转动惯量(具体参见公式(2))；利用转子运动方程式，结合起动机电流值，得到燃机剩余转矩值，燃机剩余转矩＝涡轮产生的转矩-压气机消耗的转矩-摩擦引起的转矩，筛选出燃机剩余扭矩绝对值为零的所有点，根据筛选出的燃机剩余扭矩绝对值为零的点，利用线性插值方法确认出各点相对应的燃料量，该燃料量即为对应点的“稳态燃料量”；

步骤3.2：在步骤3.1的基础上，进行多项式拟合，得到燃机转速和稳态燃料量之间的关系表达式，其关系式是一个三项多项式，具体的多项式系数因燃料种类的不同而有所差别，此处采用的燃料为天然气，其具体表达式参见公式(5)：

步骤3.3：微型燃气涡轮发动机的剩余转矩和燃料量差值之间的函数关系十分复杂，呈现高度的非线性，对其进行简化，采用分段线性函数刻画该复杂关系，假设燃机剩余转矩与燃料量差值之间是分段一次函数关系，即采用公式(6)进行刻画。

M_gas＝k(n)×(W_f-W_fss)+b(n) (6)

式中，k(n)、b(n)为系数，不同的燃机转速，其值不相同，即k(n)、b(n)数值是燃机转速n的函数，(W_f—W_fss)为动态燃料量和多项式得到的稳态燃料量之间的差值，燃机转速不同，其值也不同；

步骤3.4：微型燃气涡轮发动机的剩余转矩与起动机作用在燃机轴上的转矩求和，再利用一阶龙格-库塔方法(欧拉方法)对其进行积分，得到该模型计算出的微型燃气涡轮发动机转速，保证该微型燃气涡轮发动机的转速和试验数据中的燃机转速误差在5％以内；

步骤3.5：为了保证步骤3.4中提到的“微型燃气涡轮发动机的转速和试验数据中的燃机转速误差在5％以内”，需要根据微型燃气涡轮发动机的转速来实时调整公式(6)中的k(n)、b(n)，使构建的微型燃气涡轮发动机起动模型的转速最大程度地接近试验转速。

步骤4：第三阶段微型燃气涡轮发动机转速的建模

首先，通过计算得到发动机等效转动惯量(具体参见公式(2))；利用转子运动方程式，结合燃料量，根据公式(6)计算得到起动机加在微型燃气涡轮发动机轴上的转矩值，利用一阶龙格-库塔方法(欧拉方法)对其进行积分，得到该模型的微型燃气涡轮发动机转速，保证该微型燃气涡轮发动机的转速和试验数据中的燃机转速误差在5％以内；需要根据微型燃气涡轮发动机的转速来实时调整公式(6)中的k(n)、b(n)，使构建的微型燃气涡轮发动机起动模型的转速最大程度地接近试验转速。

步骤5：微型燃气涡轮发动机的性能参数与燃气涡轮发动机转速之间的关系十分密切，此处提到的燃气涡轮发动机的性能参数主要包括：压气机出口温度、压气机出口压力、燃烧室出口压力、涡轮出口温度、涡轮出口压力五个参数，利用MATLAB自带的数据拟合工具拟合出燃气涡轮发动机的性能参数和燃气涡轮发动机转速的函数表达式，即：

式(7)中，T₂为压气机出口温度，K_T为压气机出口的温度系数，不同的起动运行过程，其参数不同，N_g表示微型燃气涡轮发动机的转速；

式(8)中，P₂压气机出口压力,K_P为压气机出口的压力系数，不同的起动运行过程，其参数不同，N_g表示微型燃气涡轮发动机的转速；

式(9)中，P_3C为燃烧室出口压力，K_c为燃烧室的压力系数，不同的起动运行过程，其参数不同，N_g表示微型燃气涡轮发动机的转速；

P₄＝constant+k_wp×N_g (10)

式(10)中，故涡轮出口压力P₄建模为常值+k_wp×N_g，常值为标准大气压强的1.03倍，k_wp×N_g与微型涡轮发动机工作转速有关，k_wp随微型涡轮发动机转速的变化而变化；

式(11)中，T₄为涡轮出口温度，K_t为涡轮出口的温度系数，不同的起动运行过程，其参数不同，N_g表示微型燃气涡轮发动机的转速，N_gt表示微型燃气涡轮发动机的额定转速，N_idle表示微型燃气涡轮发动机的慢车转速；

步骤6：将所得的燃气涡轮发动机性能参数和转速之间的多项式函数关系进行线性化离散，即对公式(7)、(8)、(9)、(10)、(11)进行线性离散，为了保证离散精度，微型燃气涡轮发动机转速的步长值取为0.2％×N_gt，将离散化处理后的数据放入相应的一维线性插值表。

误差分析：根据燃气涡轮发动机的试验数据，得到压气机出口压力、压气机出口温度、燃烧室出口压力、涡轮出口压力、涡轮出口温度、微型燃气涡轮发动机转速等参数与燃机起动模型之间的差距均在5％以内，说明本发明所提出的微型燃气涡轮发动机起动过程的建模方法是有效的、可行的。

Claims (1)

1.一种微型燃气涡轮发动机起动过程建模方法，其特征在于，步骤如下：

步骤1：微型燃气涡轮发动机的起动过程分为3个阶段，第一阶段，发动机转速由零到涡轮开始产生功率的转速，发动机完全由起动机带动加速；第二阶段，由涡轮开始产生功率的转速到起动机脱开；第三阶段，由起动机脱开转速至慢车转速；

按照发动机起动过程的各个阶段，对其进行建模；

第一阶段，发动机转速由零到涡轮开始产生功率的转速即从零转速到点火转速，发动机完全由起动机带动加速，满足发动机转子运动方程，即满足公式(1)：

式中M_st＝K_st×I_st，其中I_st表示起动机电流值，K_st为起动机转矩常数；ω为微型燃气涡轮发动机的角速度，J为归算到微型燃气涡轮发动机轴上的转动惯量，即等效转动惯量，针对不同的微型燃气涡轮发动机的结构，利用机械能守恒进行归算；具体的归算公式(2)

式中，J₁为微型燃气涡轮发动机轴带动的负载1的转动惯量，ω₁为微型燃气涡轮发动机轴带动的负载1的角速度；J₂为微型燃气涡轮发动机轴带动的负载2的转动惯量，ω₂为微型燃气涡轮发动机轴带动的负载2的角速度；J_n为微型燃气涡轮发动机轴带动的负载n的转动惯量，ω_n为微型燃气涡轮发动机轴带动的负载n的角速度；

第二阶段，发动机转速由点火转速到起动机脱开转速，发动机在起动机和涡轮的共同带动下进行加速，满足发动机转子运动方程，即满足公式(3)：

式中，M_gas是点火后微型燃气涡轮发动机的剩余扭矩，M_gas＝M_T-M_C-M_f，M_T代表涡轮产生的转矩；M_C代表压气机消耗的转矩，与压气机转速的平方成正比；M_f表示摩擦引起的转矩消耗，与转速的平方成正比；M_st是电起动机提供的转矩，与公式(1)相同；J与公式(1)相同，ω与公式(1)相同；

第三阶段，微型燃气涡轮发动机在剩余转速的作用下，由发动机在起动机脱开转速加速慢车转速的过程，满足公式(4)：

式中，J与公式(1)相同，ω与公式(1)相同；M_gas与公式(3)相同；

步骤2：第一阶段微型燃气涡轮发动机转速的建模

首先，通过计算得到发动机等效转动惯量，见公式(2)；利用转子运动方程式，结合起动机电流值，根据M_st＝K_st×I_st计算得到起动机加在微型燃气涡轮发动机轴上的转矩值，利用一阶龙格-库塔方法对其进行积分，得到公式(1)-公式(4)建立的模型的微型燃气涡轮发动机转速，保证该微型燃气涡轮发动机的转速和试验数据中的燃机转速误差在5％以内；

步骤3：第二阶段微型燃气涡轮发动机转速的建模

步骤3.1：通过计算得到发动机等效转动惯量，见公式(2)；利用转子运动方程式，结合起动机电流值，得到燃机剩余转矩值，燃机剩余转矩＝涡轮产生的转矩-压气机消耗的转矩-摩擦引起的转矩，筛选出燃机剩余扭矩绝对值为零的所有点，根据筛选出的燃机剩余扭矩绝对值为零的点，利用线性插值方法确认出各点相对应的燃料量，该燃料量即为对应点的“稳态燃料量”；

步骤3.2：在步骤3.1的基础上，进行多项式拟合，得到燃机转速和稳态燃料量之间的关系表达式，其关系式是三项多项式，表达式(5)：

步骤3.3：微型燃气涡轮发动机的剩余转矩和燃料量差值之间的函数关系十分复杂，呈现高度的非线性，对其进行简化，因此，采用分段线性函数刻画该复杂关系，假设燃机剩余转矩与燃料量差值之间是分段一次函数关系，即采用公式(6)进行简化

M_gas＝k(n)×(W_f-W_fss)+b(n) (6)

式中，k(n)、b(n)为系数，不同的燃机转速，其值不相同，即k(n)、b(n)数值是燃机转速n的函数，(W_f-W_fss)为动态燃料量和多项式得到的稳态燃料量之间的差值，燃机转速不同，其值也不同；

步骤3.4：微型燃气涡轮发动机的剩余转矩与起动机作用在燃机轴上的转矩求和，再利用一阶龙格-库塔方法对其进行积分，得到公式(1)-公式(4)建立的模型计算出的微型燃气涡轮发动机转速，保证该微型燃气涡轮发动机的转速和试验数据中的燃机转速误差在5％以内；

步骤3.5：为了保证步骤3.4中提到的“微型燃气涡轮发动机的转速和试验数据中的燃机转速误差在5％以内”，需根据微型燃气涡轮发动机的转速来实时调整公式(6)中的k(n)、b(n)，使构建的微型燃气涡轮发动机起动模型的转速最大程度地接近实际的试验转速；

步骤4：第三阶段微型燃气涡轮发动机转速的建模

首先，通过计算得到发动机等效转动惯量，见公式(2)；利用转子运动方程式，结合燃料量，根据公式(6)计算得到起动机加在微型燃气涡轮发动机轴上的转矩值，利用一阶龙格-库塔方法对其进行积分，得到公式(1)-公式(4)建立的模型的微型燃气涡轮发动机转速，保证该微型燃气涡轮发动机的转速和试验数据中的燃机转速误差在5％以内；需根据微型燃气涡轮发动机的转速来实时调整公式(6)中的k(n)、b(n)，使构建的微型燃气涡轮发动机起动模型的转速最大程度地接近实际的试验转速；

步骤5：微型燃气涡轮发动机的性能参数与燃气涡轮发动机转速之间的关系十分密切，燃气涡轮发动机的性能参数包括压气机出口温度、压气机出口压力、燃烧室出口压力、涡轮出口温度和涡轮出口压力五个参数，拟合燃气涡轮发动机的性能参数和燃气涡轮发动机转速的函数表达式，即：

式(7)中，T₂为压气机出口温度，K_T为压气机出口的温度系数，不同的起动运行过程，其参数不同，N_g表示微型燃气涡轮发动机的转速；

式(8)中，P₂压气机出口压力,K_P为压气机出口的压力系数，不同的起动运行过程，其参数不同，N_g表示微型燃气涡轮发动机的转速；

式(9)中，P_3C为燃烧室出口压力，K_c为燃烧室的压力系数，不同的起动运行过程，其参数不同，N_g表示微型燃气涡轮发动机的转速；

P₄＝constant+k_wp×N_g (10)

式(10)中，故涡轮出口压力P₄建模为常值+k_wp×N_g，常值为标准大气压强的1.03倍，k_wp×N_g与微型涡轮发动机工作转速有关，k_wp随微型涡轮发动机转速的变化而变化；

式(11)中，T₄为涡轮出口温度，K_t为涡轮出口的温度系数，不同的起动运行过程，其参数不同，N_g表示微型燃气涡轮发动机的转速，N_gt表示微型燃气涡轮发动机的额定转速，N_idle表示微型燃气涡轮发动机的慢车转速；

步骤6：将所得的燃气涡轮发动机性能参数和转速之间的多项式函数关系进行线性化离散，即对公式(7)、(8)、(9)、(10)、(11)进行线性离散，为了保证离散精度，微型燃气涡轮发动机转速的步长值取为0.2％×N_gt，将离散化处理后的数据放入相应的一维线性插值表。

Images