使用GasTurb 13/14开展航空发动机/燃气轮机性能分析：以AMT Lynx涡喷发动机为例，设计点分析与模型修正

1 准备工作

1.1 GasTurb 13/14 航空发动机/燃气轮机性能仿真与分析软件

1.2 荷兰AMT公司的160kgf推力级涡喷发动机Lynx技术参数

山猫（Lynx）发动机是AMT-NL Nike设计的衍生型号，采用单级离心式压气机、直流环形燃烧室和单级轴流涡轮结构。其卓越性能和优异的推重比归功于新型扩压器设计。
该革命性扩散器技术同时实现了发动机的紧凑直径。

1.2.1 初步规格 (设计转速)

所有数据基于标准温压条件±2%

参数	公制单位	英制单位
发动机直径	260毫米	10.2英寸
发动机长度	690毫米	27.1英寸
发动机重量	17010克	37.5磅
系统飞行重量*	19500克	43磅
设计转速推力(46,000 RPM)	1569牛	352.0磅力
最大允许转速	50,000 RPM	50,000 RPM
怠速推力	70牛	15.7磅力
压比	4:1	4:1
设计转速质量流量	2200克/秒	4.84磅/秒
常规排气温度(双探头平均)	800°C	1472°F
最大排气温度(双探头平均)	875°C	1607°F
燃油消耗率	3600克/分	127盎司/分
燃油消耗系数	38.24克/(千牛·秒)	1.35磅/(磅力·小时)
启动方式	直喷煤油启动系统(<25秒启动时间)

*总重量包含：发动机、ECU、油泵、锂聚合物电池、温度传感器、阀门、安装支架

1.2.2 工作条件

参数	公制单位	英制单位
最低工作温度*	-25°C	-13°F
最高工作温度	+50°C	122°F
最大工作海拔	8000米	26250英尺
最大工作速度	1000公里/小时	620英里/小时

*使用锂聚合物电池时最低工作温度限制为-10°C

1.2.3 控制接口

接口类型	功能说明
1 x RS232	节流阀与开关设置(含特殊功能*)
2 x PWM(1-2毫秒)**	节流阀与开关设置(按钮激活特殊功能*)
2 x 模拟输入(0-5V)**	节流阀与开关设置(按钮激活特殊功能*)

*特殊功能包括：预注油、点火器、启动器、快速冷却

NL7001-山猫发动机测试数据
所有数据基于标准温压条件(15°C/1013毫巴)

时间戳	油门(%)	转速(千转)	排气温度(°C)	输出电压(V)	推力(牛)	总燃油(cc/分)	燃烧室(cc/分)	轴承(cc/分)	耗油系数(克/千牛·秒)
2013-06-13 11:52:35	0	17.5	649	1.37	110	923	902	21	105.67
2013-06-13 11:52:07	10	25.0	640	1.83	254	1418	1378	40	70.38
2013-06-13 11:51:47	20	30.0	635	2.29	418	1811	1753	58	54.70
2013-06-13 11:51:28	30	34.0	635	2.71	585	2175	2101	74	46.88
2013-06-13 11:51:09	40	37.0	644	3.14	756	2543	2454	89	42.43
2013-06-13 11:50:50	50	39.0	672	3.56	912	2881	2779	102	39.86
2013-06-13 11:50:30	60	40.5	690	3.99	1062	3229	3114	115	38.36
2013-06-13 11:50:11	70	42.0	718	4.44	1197	3543	3419	124	37.33
2013-06-13 11:49:52	80	43.0	736	4.87	1327	3864	3732	132	36.74
2013-06-13 11:49:33	90	44.0	750	5.29	1441	4154	4016	138	36.37
2013-06-13 11:49:14	100	45.0	778	5.75	1558	4438	4293	145	35.93

AMT公司官网：https://amtjets.com/

2 热力循环设计点性能仿真与分析

2.1 设计点输入参数

Property	Unit	Value	Comment
Inlet Corr. Flow W2Rstd	kg/s	2.5
Intake Pressure Ratio		1
Pressure Ratio		4
Isentr.Compr.Efficiency		0.78
Burner Exit Temperature	K	1219
Burner Design Efficiency		0.91
Burner Partload Constant		1.6	used for off design only
Fuel Heating Value	MJ/kg	42.9
Overboard Bleed	kg/s	0
Isentr.Turbine.Efficiency		0.84
Power Offtake	kW	0
Mechanical Efficiency		0.99
Burner Pressure Ratio		0.94
Turbine Exit Duct Press Ratio		1
Speicified Thrust Coeff CFG		0.985
Isentr.Compr.Efficiency		0.81

2.2 仿真结果

3 模型修正的初步开展

下面将分别初步介绍对发动机耗油率影响的燃烧室效率修正方法和对推力影响的尾喷管推力系数修正方法。

3.1 燃烧室燃烧效率特性修正

3.1.1 GasTurb采用的燃烧室效率计算方法总结

3.1.1.1 核心公式

a. 燃烧室负荷参数

$\Omega =\frac{W_{31}}{P_3^{1.8}\cdot e^{\frac{T_3}{300K}}\cdot Vol}$

参数说明：

• $W_{31}$：空气质量流量（kg/s）
• $P_3$：燃烧室入口压力（bar）
• $T_3$：燃烧室入口温度（K）
• $Vol$：燃烧室体积（m³）

b. 相对燃烧室负荷

$\text{相对负荷}=\frac{\Omega }{{\Omega }_{\text{des}}}$

• ${\Omega }_{\text{des}}$：设计点的燃烧室负荷。

c. 燃烧效率计算公式

$\log \left(1-\eta \right)=a+b\cdot \log \left(\frac{\Omega }{{\Omega }_{\text{des}}}\right)$
其中：

• $a=\log \left(1-{\eta }_{\text{des}}\right)$
• ${\eta }_{\text{des}}$：设计点燃烧效率（输入值）
• $b$：燃烧室部分负荷常数（默认值1.6，建议取值范围1.0~2.2）

d. 燃烧效率反推公式

$\eta =1-10^{a+b\cdot \log \left(\frac{\Omega }{{\Omega }_{\text{des}}}\right)}$

3.1.1.2 关键参数说明

参数	说明
${\eta }_{\text{des}}$	设计点燃烧效率（如输入0.999表示99.9%效率）
$b$	燃烧室部分负荷常数，默认值为1.6，反映效率随负荷变化的敏感性
$\Omega$	燃烧室负荷参数，与流量、压力、温度、体积相关
$\Theta$	燃烧室负荷参数的倒数，$\Theta =1/\Omega$

3.1.1.3 计算步骤

a. 设计点计算

• 输入设计点参数（$W_{31},P_3,T_3,Vol$）计算${\Omega }_{\text{des}}$。
• 输入设计点燃烧效率${\eta }_{\text{des}}$，计算常数$a=\log (1-{\eta }_{\text{des}})$。

b. 非设计点计算

• 根据部分负荷参数计算当前$\Omega$。
• 计算相对负荷$\Omega /{\Omega }_{\text{des}}$。
• 使用公式$\eta =1-10^{a+b\cdot \log (\Omega /{\Omega }_{\text{des}})}$得到非设计点效率。

3.1.1.4 图表辅助分析

• 效率与负荷关系图：通过调整$b$值（部分负荷常数），可模拟不同燃烧室的效率变化趋势。
• 默认值适用性：GasTurb 13默认使用$b=1.6$，若需精确匹配实际燃烧室性能，需通过实验或文献数据校准。

3.1.1.5 注意事项

• 燃烧室体积$Vol$为固定值，设计点与部分负荷计算中保持不变。
• 若默认方法不适用，可通过自定义表格（General Table）结合迭代方法调整燃烧效率。

3.1.2 GasTurb燃烧效率特性修正方法

若GasTurb提供的上述燃烧效率模拟方法不适用于您的燃烧室类型，可采用以下步骤进行燃烧效率计算：

a. 使用组合值与通用表格

通过构建通用表格（General Table）定义燃烧效率与工况参数的映射关系。

b. 迭代燃烧室设计效率

将"Burner Design Efficiency"设为迭代变量，通过循环计算使其满足：
${\eta }_{\text{calculated}}={\eta }_{\text{composed value}}$

重要说明
在非设计点模拟中，"Burner Design Efficiency"需先被显式定义为输入量（Input Quantity），方可激活为可迭代变量。此操作需在GasTurb软件中通过勾选参数输入状态实现。

3.2 实现自定义的喷管推力系数特性

3.2.1 GasTurb采用的推力系数定义

推力系数是用于修正总推力速度项的修正因子，用于计算不同喷管条件下的净推力。以下是两类推力系数的具体定义：

3.2.1.1 实际推力系数

a. 应用场景

适用于分排涡扇发动机的净推力计算，包含核心喷管和外涵喷管的推力系数。

b. 公式

$F_N=C_{FG,8}\cdot \left(\frac{W_8}{W_{ref}}\cdot \sqrt{\frac{P_{s,8}}{P_{amb}}}\cdot A_8\cdot V_8\right)+C_{FG,18}\cdot \left(\frac{W_{18}}{W_{ref}}\cdot \sqrt{\frac{P_{s,18}}{P_{amb}}}\cdot A_{18}\cdot V_{18}\right)$

参数说明：

• $C_{FG,8}$：核心喷管推力系数
• $C_{FG,18}$：外涵喷管推力系数
• $W_8,W_{18}$：核心/外涵空气质量流量（kg/s）
• $P_{s,8},P_{s,18}$：核心/外涵喷管静压（Pa）
• $P_{amb}$：环境压力（Pa）
• $A_8,A_{18}$：喷管截面积（m²）
• $V_8,V_{18}$：喷管出口流速（m/s）

3.2.1.2 理想推力系数

a. 应用场景

描述收敛-扩张喷管在完全膨胀到环境压力时的理想情况。

b. 公式

$C_{FG,id}=\frac{F_G}{F_{G,id}}$
其中：

• $F_{G,id}$：喷管完全膨胀至环境压力时的理想总推力
• $F_G$：实际喷管出口面积$A_9$膨胀后的总推力

3.2.1.3 关键区别

类型	描述
实际推力系数	用于实际工程计算，区分核心与外涵喷管的推力贡献
理想推力系数	理论参考值，反映喷管设计的理想性能，用于评估实际推力损失

3.2.1.4 使用注意事项

1. 符号一致性：避免混淆$C_{FG}$（实际推力系数）与$C_{FG,id}$（理想推力系数）。
2. 单位统一：所有压力参数需统一单位（如均用Pa或bar）。
3. 适用范围：理想推力系数仅适用于收敛-扩张喷管的理论分析，实际计算需结合具体喷管类型调整。

3.2.2 GasTurb自定义喷管推力系数特性方法

a. 使用组合值与通用表格

通过构建通用表格（General Table）定义喷管推力系数与工况参数（例如，喷管落压比或压气机换算转速）的映射关系。

b. 迭代喷管推力系数

将"CFG8"设为迭代变量，通过循环计算使其满足：
$CFG{8}_{\text{input}}=CFG{8}_{\text{composed value}}$

4 模型计算与参考对比

若需进一步咨询，可将需求信息发送至邮箱：pyengine@163.com