燃气涡轮发动机起动建模方法与流程

本发明涉及燃气涡轮发动机建模仿真技术领域，具体涉及一种燃气涡轮发动机起动建模方法。

背景技术：

在航空燃气涡轮发动机起动控制规律和控制算法的设计开发过程中，需要较为准确的发动机起动模型。如果采用传统的机理建模方法，需要通过试验精确测定出起动过程中的气动部件特性、燃烧效率特性、风车特性和内外涵匹配特性等，需要专门的试验设备、较高的试验成本和较长的周期，对于许多小型发动机而言根本不具备以上试验条件。

技术实现要素：

(一)要解决的技术问题

本发明要解决的技术问题是：如何设计一种简便易行，计算量小，对建模前提条件要求低，模型精度可满足工程使用的需要的燃气涡轮发动机起动建模方法。

(二)技术方案

为了解决上述技术问题，本发明提供了一种燃气涡轮发动机起动建模方法，包括以下步骤：

步骤1、根据发动机试验数据提取高转折合转加速度变化规律：

高转折合转加速度是高转折合转速n2c和折合燃油流量wfc的函数，即有对于i＝1,2,3...n次试验，其中第i次试验中记录了j＝1,2,3...mi个采样时刻的数据，mi为第i次试验中采样时刻总数，采样间隔为dt，对于第i次试验中的第j个采样时刻，进行以下计算：

根据发动机风扇入口总温tt1i,j计算温度折合系数θi,j

θi,j＝tt1i,j/tstd(1)

其中tstd为海平面标准大气温度288.15k；

根据发动机风扇入口总压pt1i,j计算压力折合系数δi,j

δi,j＝pt1i,j/pstd(2)

其中pstd为海平面标准大气压力10132.5pa；

使用温度折合系数θi,j、高转转速n2i,j计算高转折合转速n2ci,j

使用温度折合系数θi,j、压力折合系数δi,j、燃油流量wfi,j计算折合燃油流量wfci,j为：

计算高转转加速度

ni,j+1为第i次试验中的第j+1个采样时刻的发动机转速，ni,j-1为第i次试验中的第j-1个采样时刻的发动机转速；

使用高转转加速度压力折合系数δi,j计算高转折合转加速度

提取风车起动过程高转折合转加速度规律时，对使用风车起动试验数据计算得到的n2ci,j和wfci,j，进行二元多项式拟合，即得到

根据起动机瞬时有效功率psi,j和高转转子转动力矩j计算由起动机贡献的高转转加速度

从中扣除得到高转净转加速度

则高转折合净转加速度为

使用发动机地面起动试验数据计算得到的n2ci,j和wfci,j，进行二元多项式拟合，即得到

步骤2、根据发动机试验数据提取低转折合转速n1c变化规律：

低转折合转速n1c是高转折合转速n2c的函数，即有n1c＝g(n2c)；

使用第i次试验中的第j时刻的低转转速n1i,j计算高转折合转速n1ci,j

使用发动机地面起动试验数据计算得到的n2ci,j和n1ci,j，进行多项式拟合，即可得到n1c＝g(n2c)；

步骤3、根据发动机试验数据提取高压压气机出口折合压力pt3c变化规律：

高压压气机出口折合压力pt3c是高转折合转速n2c的函数，即有pt3c＝h(n2c)；

使用第i次试验中的第j时刻的高压压气机出口压力pt3i,j计算高压压气机出口折合压力pt3ci,j

pt3ci,j＝pt3i,j/δi,j(11)

使用发动机地面起动试验数据计算得到的n2ci,j和pt3ci,j，进行多项式拟合，即可得到pt3c＝h(n2c)；

然后利用步骤1至3的规律建立发动机起动过程建模。

优选地，慢车过程结束时仿真结束。

优选地，发动机转速达到慢车转速时仿真结束。

优选地，然后利用步骤1至3的规律建立发动机起动过程建模的计算过程包括以下步骤：

步骤4、模型初始化：高转转速的初始值n20＝n2ignit，其中，n2ignit为发动机点火转速；

对于第j时刻根据该时刻模型的输入参数：发动机风扇入口总温tt1j、发动机风扇入口总压pt1j、燃油流量wfj和起动机有效输出功率psj进行步骤5的计算；

步骤5、计算j时刻温度折合系数θj和压力折合系数δj

根据tt1j计算θj为：

θj＝tt1j/tstd(12)

根据pt1j计算δj为：

δj＝pt1j/pstd(13)

根据j时刻高转转速n2j计算j时刻高转折合转速n2cj为：

利用拟合得到的低转折合转速特性n1c＝g(n2c)计算j时刻低转折合转速为：

n1cj＝g(n2cj)(15)

计算j时刻低转转速n1j为：

利用拟合得到的高压压气机出口折合压力特性pt3c＝h(n2c)计算j时刻高压压气机折合压力pt3cj为：

pt3cj＝h(n2cj)(17)

计算j时刻高压压气机压力pt3j为：

pt3j＝pt3cjδj(18)

计算j时刻折合燃油流量wfcj为：

利用拟合得到的高转折合转加速度变化规律计算j时刻高转折合转加速度，在风车起动时为：

计算j时刻高转转加速度为：

在地面起动时，先计算未考虑起动机作用的高转净折合转加速度

计算j时刻高转净转加速度为

再补偿起动机拖转作用，得到高转转加速度

其中psj为j时刻的起动机有效输出功率；

计算j+1时刻高转转速n2j+1

之后重复步骤5的计算过程，直至仿真结束。

优选地，所述燃气涡轮发动机为航空燃气涡轮发动机。

优选地，所述燃气涡轮发动机为涡扇发动机。

(三)有益效果

本发明提供一种基于相似理论的燃气涡轮发动机起动建模方法，该方法包括：基于相似理论的高转转加速度、高转转速、低转转速和高压压气机出口压力计算方法；以及使用上述计算结果根据核心机气动热力学过程计算涡轮级间温度的方法。该方法能够以较高精度模拟发动机起动过程的主要参数，仅基于发动机试验数据建立发动机起动模型，具有方法简便，计算量小，对建模前提条件要求低的优点，模型精度可满足工程使用的需要。

附图说明

图1是采用本发明的方法得到的3km高度风车起动过程仿真结果；

图2是采用本发明的方法得到的4.5km高度风车起动过程仿真结果。

具体实施方式

为使本发明的目的、内容、和优点更加清楚，下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。

本发明实施例以涡扇发动机为例，提出一种基于相似理论的燃气涡轮起动建模方法，包括以下步骤：

步骤1、根据发动机试验数据提取高转折合转加速度变化规(折合到标准大气条件下)。

高转折合转加速度是高转折合转速n2c和折合燃油流量wfc的函数，即有对于i＝1,2,3...n次试验，其中第i次试验中记录了j＝1,2,3...mi个采样时刻的数据(mi为第i次试验中采样时刻总数)，采样间隔为dt。对于第i次试验中的第j个采样时刻，进行以下计算：

根据发动机风扇入口总温tt1i,j计算温度折合系数θi,j

θi,j＝tt1i,j/tstd(1)

其中tstd为海平面标准大气温度288.15k。

根据发动机风扇入口总压pt1i,j计算压力折合系数δi,j

δi,j＝pt1i,j/pstd(2)

其中pstd为海平面标准大气压力10132.5pa。

使用温度折合系数θi,j、高转转速n2i,j计算高转折合转速n2ci,j

使用温度折合系数θi,j、压力折合系数δi,j、燃油流量wfi,j计算折合燃油流量wfci,j为：

计算高转转加速度

ni,j+1为第i次试验中的第j+1个采样时刻的发动机转速，ni,j-1为第i次试验中的第j-1个采样时刻的发动机转速；

使用高转转加速度压力折合系数δi,j计算高转折合转加速度

提取风车起动过程高转折合转加速度规律时，对使用风车起动试验数据按照上述公式(3)、(4)、(6)计算得到的n2ci,j和wfci,j，进行二元多项式拟合，即可得到

地面起动过程中使用起动机带转高转转子，因此需要扣除起动机的影响，提取地面起动过程中的高转净折合转加速度变化规律根据起动机瞬时有效功率psi,j和高转转子转动力矩j计算由起动机贡献的高转转加速度

从中扣除得到高转净转加速度

则高转折合净转加速度为

使用发动机地面起动试验数据计算得到的n2ci,j和wfci,j，进行二元多项式拟合，即可得到

步骤2、根据发动机试验数据提取低转折合转速n1c变化规律

低转折合转速n1c是高转折合转速n2c的函数，即有n1c＝g(n2c)。

使用第i次试验中的第j时刻的低转转速n1i,j计算高转折合转速n1ci,j

使用发动机地面起动试验数据计算得到的n2ci,j和n1ci,j，进行多项式拟合，即可得到n1c＝g(n2c)。

步骤3、根据发动机试验数据提取高压压气机出口折合压力pt3c变化规律。

高压压气机出口折合压力pt3c是高转折合转速n2c的函数，即有pt3c＝h(n2c)。

使用第i次试验中的第j时刻的高压压气机出口压力pt3i,j计算高压压气机出口折合压力pt3ci,j

pt3ci,j＝pt3i,j/δi,j(11)

使用发动机地面起动试验数据计算得到的n2ci,j和pt3ci,j，进行多项式拟合，即可得到pt3c＝h(n2c)。

在完成上述规律的提取后，使用上述规律建立发动机起动过程建模，模型计算过程包括以下步骤：

步骤4、模型初始化：高转转速的初始值n20＝n2ignit,其中，n2ignit为发动机点火转速。

对于第j时刻根据该时刻模型的输入参数：发动机风扇入口总温tt1j、发动机风扇入口总压pt1j、燃油流量wfj和起动机有效输出功率psj进行步骤五的计算。

步骤5、计算j时刻温度折合系数θj和压力折合系数δj

根据tt1j计算θj为：

θj＝tt1j/tstd(12)

根据pt1j计算δj为：

δj＝pt1j/pstd(13)

根据j时刻高转转速n2j计算j时刻高转折合转速n2cj为：

利用拟合得到的低转折合转速特性n1c＝g(n2c)计算j时刻低转折合转速为：

n1cj＝g(n2cj)(15)

计算j时刻低转转速n1j为：

利用拟合得到的高压压气机出口折合压力特性pt3c＝h(n2c)计算j时刻高压压气机折合压力pt3cj为：

pt3cj＝h(n2cj)(17)

计算j时刻高压压气机压力pt3j为：

pt3j＝pt3cjδj(18)

计算j时刻折合燃油流量wfcj为：

利用拟合得到的高转折合转加速度变化规律计算j时刻高转折合转加速度。在风车起动时为：

计算j时刻高转转加速度为：

在地面起动时，先计算未考虑起动机作用的高转净折合转加速度

计算j时刻高转净转加速度为

再补偿起动机拖转作用，得到高转转加速度

其中psj为j时刻的起动机有效输出功率。

计算j+1时刻高转转速n2j+1

之后重复步骤5的计算过程，直至仿真结束，本实施例中，发动机转速达到慢车转速时，即慢车过程结束时仿真结束。

采用本发明的方法进行3km高度风车起动的仿真，仿真与试验对比结果如图1所示。采用本发明的方法还进行了4.5km高度地面起动的仿真，仿真与试验对比结果如图2所示。仿真从点火时刻开始运行，运行周期间隔为0.01s。可以看出，n2和pt3的模型示值误差可以控制在10％以内；n1起动初期误差较大，在第8s后示值误差也可以降低到10％以下。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。