专栏内容(赛前预售价99,比赛期间299): 2024数学建模国赛期间会发布思路、代码和优秀论文。(本专栏达不到国一的水平,适用于有一点点基础冲击省奖的同学,近两年有二十几个国二,但是达不到国一,普遍获得省奖,请勿盲目订阅)
python全套教程(一百篇博客):从新手到掌握使用python,可以对数学建模问题进行建模分析。
35套模型算法(优秀论文示例):马尔科夫模型、遗传算法、逻辑回归、逐步回归、蚁群算法、蒙特卡洛模拟、聚类模型、线性规划、粒子群算法、神经网络、相关系数、灰色预测、灰色关联分析、模糊综合评价、模拟退火、时间序列ARMA、方差分析支持向量机、插值、排队论、拟合模型、微分方程、层次分析法、小波分析、多元回归、图论floyd算法、图论Dijkstra模型、因子分析、动态规划、博弈论、决策树、典型相关分析、元胞自动机、主成分分析、TOPSIS法。
- 1 -
参赛密码
(由组委会填写)
第十届华为杯全国研究生数学建模竞赛
题 目
《变循环发动机部件法建模及优化》
摘 要:
本文研究的是变循环发动机部件法建模及优化问题。
对于第一问,首先要求分析风扇特性,通过部件法建模,风扇属于压气机
模型,根据压气机压比函数的定义,利用 matlab 线性插值将附录 4 中增压比转
换成压比值,作出在不同转换速度下的流量随压比函数值的曲线图,如图 2 所
示。为计算风扇和 CDFS 出口的总温、总压和流量,根据附录 1 的计算公式,构
建进气道模型和风扇、CDFS 的压气机模型。根据已知数据,在亚音速状态下代
入构建的进气道、风扇和 CDFS 计算模型即可得到所求变量值:风扇出口总压为
1.2883、总温为
378.3326
、流量分别为
19.0477
;
CDFS
的出口总温为
419.3767
,
总压为
1.7745
,出口流量为
16.9401
。其余中间变量见表
1
至表
3
所示。
第二问,要求解 7 个平衡方程组成的非线性方程。首先,分析发动机的整
机模型,根据部件法构建发动机的各部分计算模型。在发动机各部件匹配工作
时,受 7 个平衡方程制约。根据已知的参数和发动机模型,设定七个适当的未
知数,分别为:高压转速
n
H
、风扇、CDFS、高压涡轮、低压涡轮的压比函数值:
Z
CL
、
Z
CDFS
、
Z
CH
、
Z
TH
、
Z
TL
以及主燃烧室的出口温度
T
4 *
。由构建的发动机模
型与其工作时各部件匹配的 7 个平衡方程可得到一组非线性方程组。为求解此
方程组,采用粒子群算法。将非线性方程组的潜在解视为粒子,设定的未知量可
看作粒子的位置,粒子群中的粒子个体根据适应性函数引导进行迭代,向非线
性方程组的最优解不断逼近 ,通过迭代法,更新粒子的位置,以搜索到全局极
值,即获取非线性方程组的最优解。经过多次试验,粒子群算法均能搜寻到能
得到比较准确的非线性方程的解。在给定不同的群体数和初值时,可得到以下
比较理想的解:
- 2 -
变量
n
H
T
4 *
Z
CL
Z
CDFS
含义
高压转速
主燃烧室出口
温度
风扇压比函数值
CDFS
压比函
数值
最优解
0.78
1369.9999
0.3394
1
变量
Z
CH
Z
TH
Z
TL
含义
高压压气机压
比函数值
高压涡轮压比函数
值
低压涡轮压比函数
值
最优解
0.2899
0.246
0.9112
对于第三问,求解发动机性能最优时,发动机CDFS导叶角度、低压涡轮导叶
角度和喷管喉道面积3个量的取值。本题的求解为单涵道模式下最优解的查找问
题,可按照与问题二类似的步骤,首先根据发动机的整体模型和部件法构建发动
机的各个部分计算模型;在单涵道工作模式下,在发动机各部件匹配工作时,受
8个平衡方程制约。设定相应的未知量,建立关于发动机
CDFS
导叶角度
CD
、低
压涡轮导叶角度
l
以及喷管喉道面积
A
这三个位置参量的非线性方程组,求出
其他未知量与
CD
,
l
,
A
的关系。进而表示出发动机性能参数推力
F
,耗油率
sfc
,单位推力
Fs
的函数关系,另外根据性能最优要求,建立性能评价函数,得到
关于
CD
,
l
,
A
的三元函数,通过粒子群算法,求解函数的极值,便是性能最
优时的参数设定。但鉴于时间有限,本文中只建立模型,分析了求解算法,后续
的求解过程和结果仍在研究过程中。
- 3 -
一、问题重述
1.1 问题背景
由飞机
/
发动机设计原理可知,对于持续高马赫数飞行任务,需要高单位推
力的涡喷循环,反之,如果任务强调低马赫数和长航程,就需要低耗油率的涡
扇循环。双涵道变循环发动机可以同时具备高速时的大推力与低速时的低油耗。
变循环发动机的内在性能优势,受到了各航空
s
强国的重视,是目前航空发动
机的重要研究方向。
双涵道变循环发动机的基本构造如图 1 所示:
进气道
风扇
副外涵道
前混合器
核心驱动风扇级(
CDFS
)
主外涵道
高压压气机
主燃烧室
高压涡轮
低压涡轮
后混合器
加力燃烧室
尾喷管
模式转换活门
涡
扇工
作模式(
图上半部分)
涡喷工作模式(图下半部分)
图
1
双涵道变循环发动机基本构造图
1.2 问题的提出
根据变循环发动机的构造及基本原理,建立数学模型完成以下任务:
1. 1)根据附录 4 中风扇特性数据表中流量随压比函数值变化的图形;
2)设在发动机飞行高度
11
H km
,飞行马赫数
Ma
0.8
的亚音速巡航
点,导叶角度均设置为
0°
,风扇和
CDFS
的物理转速都为
0.95
,风扇和
CDFS
的压比函数值都为
0.5
,
求风扇和
CDFS
的出口总温、总压和流量。
2. 设在发动机飞行高度
11
H km
,飞行马赫数
Ma
0.8
的亚音速巡航点,
采用双涵道模式,导叶角度均设置为
0°
,选择活门完全打开,副外涵道
面积设为
1.8935e+003
,后混合器出口总面积设置为
2.8518e+004
,尾喷
管喉道面积
A
8
9.5544e+003
,
n
L
=0.85
。请运用或设计适当的算法求解
由发动机
7
个平衡方程组成的非线性方程组。要求陈述算法的关键步骤
及其解释,尽可能讨论算法的有效性。
3.
1
)设在发动机飞行高度
11
H km
,飞行马赫数
Ma
1.5
的超音速巡航
点,发动机采用单涵道模式,将选择活门面积设置为
0
,风扇导叶角度、
高压压气机导叶角度、高压涡轮导叶角度均设置为
0°
,后混合器面积设
置为
2.8518e+004
。请问发动机
CDFS
导叶角度、低压涡轮导叶角度和
喷管喉道面积
3
个量为多少时,发动机的性能最优
?
2)试研究发动机飞行高度
11
H km
,飞行马赫数从
Ma
1.1
变化到
Ma
1.6
,
发动机特性最优时,CDFS 导叶角度、低压涡轮导叶角度,尾喷
管喉道面积随飞行马赫数的变化规律。此时发动机采用单涵道模式,将
选择活门面积设置为
0
,风扇导叶角度、高压压气机导叶角度、高压涡
轮导叶角度均设置为
0°
,后混合器出口总面积设置为
2.8518e+004
,后
混合器内、外涵道面积可调。
更多数学建模资料请关注微店店铺“数学建模学习交流”
https://k.weidian.com/RHO6PSpA
- 4 -
二、模型的基本假设
2.1
系统处于稳定状态,不考虑单/双涵道模式转换时的动态过程;
2.2
变循环发动机是十分复杂的气动热力学系统,为简化模型,对部分部件进
行适当的简化;
2.3
忽略燃油在燃烧室中燃烧延迟的影响;
2.4
气流在各部件的流动都按准一维流动处理,即各个截面气流参数作平均处
理[1];
2.5
外界环境气压为标准大气压;
三、符号说明
注释表
zz
:
压比
pr
:
增压比
H
:
发动机所在高度
p
0
:
静压
T
0
:
静温
M
a
:
飞行马赫数
:
气体绝热指数
i
:
总压恢复系数
n
cor
:
转换速度
n
:
物理速度
pr
:
增压比
W
:
换算流量
k
:
修正系数
R
:
气体常数
:
熵
h
:
焓
T
:
总温
p
:
总压
f
:
油气比
:
导叶角
N
:
功
l
:
功率
A
:
面积
( )
j
F X
:
适应度函数
下标
map
:
插值
in
:
进口
out
:
出口
pr
:
增压
/
落压比
ei
:
理想
max
:
最大
min
:
最小
air
:
空气
gas
:
燃气
- 1 -
四、问题分析
本文研究的是变循环发动机部件法的建模及优化问题,根据各部件的已知
特性,模拟整个发动机的性能。
问题一,首先要求分析风扇特性。在变循环发动机的工作过程中,风扇相当
于低压压气机。根据附录 3 中压气机压比函数的定义式,利用 matlab 对附录 4
中原有数据进行处理,将增压比转换成压比值,作出在不同转换速度下的流量
随压比函数值的曲线图。为计算风扇和 CDFS 出口的总温、总压和流量,采用部
件法,根据附录 1 中给出的发动机部件计算公式,构建进气道和压气机模型。
由于风扇和 CDFS 均相当于压气机,在建模时可直接采用压气机模型。根据已知
数据,在亚音速状态下代入构建的三个部件的计算模型即可得到所求变量值。
问题二要求解 7 个平衡方程组成的非线性方程。首先,分析发动机的整机
模型,根据部件法构建发动机的各部分计算模型。在发动机各部件匹配工作时,
受 7 个平衡方程制约。设定适当的未知参数,由构建的发动机模型可将平衡方
程等价于关于未知参数的非线性方程组。为求解此方程组,采用粒子群算法,将
非线性方程组的潜在解视为粒子,设定的未知量可看作粒子的位置,粒子群中
的粒子个体根据适应性函数引导进行迭代,向非线性方程组的最优解不断逼
近 ,通过迭代法,更新粒子的位置,以搜索到全局极值,即获取非线性方程组
的最优解。
问题三要求解发动机性能最优时,发动机
CDFS
导叶角度、低压涡轮导叶
角度和喷管喉道面积
3
个量的取值。这个问题涉及到最优解的问题,根据题目
意思,可知,最优解时,发动机有较强的动力性和燃油经济性。并且,发动机
的新能参数推力
F
,耗油率
sfc
,单位推力
Fs
可以用发动机
CDFS
导叶角度、
低压涡轮导叶角度和喷管喉道面积这三个变量表示,问题转化成关于这三个变
量的发动机性能参数的最优解问题。所以首先需要做的是将性能参数用这三个
变量表示出,同问题二一样的做法,要根据发动机的整体模型和部件法构建发
动机的个个部分计算模型,注意此时是超音速飞行,飞机需要有较强的动力,
所以此时飞机发动机工作在单涵道模式,部件法建模计算时,同样受到平衡方
程制约,在单涵道模式下,没有混合器流量制约,但是单涵道模式中,可以认
为,所有气体流量全部走内涵道,那么有前后级部件流量是相等的。设定相应
的未知量以及中间变量,构建未知参数的非线性方程组,最终可以消除所有中
间变量,将发动机性能参数用题目中的三个变量表示,然后再求解函数的最优
解。对于函数最优解问题,仍然可以用粒子群算法。
五、模型的建立与问题求解
5.1
问题
1
的求解
5.1.1
风扇特性分析
在变循环发动机工作中,气体进入进气道后,经过风扇,气压升高。风扇
的工作原理相当于低压压气机。根据附录 3 可知,压气机的压比的物理含义为:
在某一换算转速下,压气机增压比的标准变化量。由此定义该换算转速下,风
- 2 -
扇的增压比与压比函数关系式:
min
max min
zz=
pr pr
pr pr
(1)
根据式(1),利用 excel 直接计算出附录 4 风扇特性数据表中各个换算转速
下每个增压比对应的压比值。再将所得的压比值及对应的流量数据导入matlab,
拟合出不同换算转速下流量随压比函数值的曲线,如图 2 所示:
图 2 各换算转速下流量随压比函数值的曲线
5.1.2
风扇与
CDFS
的出口总温、总压及流量计算
气体进入进气道后,流经风扇,由风扇对气流进行初步增压;部分气流再
经由 CDFS 对其进一步做功,使得气流的压力和温度进一步升高。
采用部件法,根据附录 1 中给出的发动机各部件计算公式,构建进气道及
压气机模型,而风扇和 CDFS 均相当于压气机。根据飞行高度,飞行马赫数,导
叶角和 CDFS 的物理转速,风扇和 CDFS 的压比函数值等已知数据,按照构建的
模型分成进气道、风扇和 CDFS 三步进行计算,即可得到风扇和 CDFS 出口的总
温、总压和流量。
5.1.2.1
进气道模型
进气道的作用是提供与发动机相匹配的空气流量。根据附录 1 进气道的计
算模型,在已知
H km
11
时,可得到气道出口的总温、总压。
1)根据已知的高度 H,可计算标准大气条件下环境压力
p
0
(静压),环境
- 3 -
温度
T
0
(静温):
图 5 3 双涵道变循环发动机模型气流流动示意图
常规混排涡扇发动机不同的是,双涵道变循环发动机具有两种典型的工作
模式:单涵模式和双涵模式。
双涵道模式
当处于起飞或亚音速巡航时,发动机以此模式工作,气流流动过程为图 4,
发动机在亚音速巡航的低功率工作状态,风扇后的模式转换活门因为副外涵与
风扇后的压差打开,使更多空气进入副外涵,同时前混合器面积开大,打开后
混合器,增大涵道比,降低油耗,此时为发动机的涡扇模式。
单涵道模式
在飞机加速、爬升以及超音速飞行时,发动机将以此模式工作,气流流动
过程为图 4 虚线框部分。关闭模式选择活门,关小前、后 混合器 出口面积,
使前段风扇出口空气几乎全部通过核心驱动风扇级和高压压气机,仅少量空气
经外涵道流出以冷却尾喷管,这时发动机以涡喷的工作方式工作,产生高单位
推力,以满足上述状态对推力的需求。[1]
5.2.1.2
构建发动机其余各部件模型
根据图 5 中的整机模型,采用部件法对各个模块建立计算模型。其中进气
道与压气机模型在 5.1 中已建立,这里不再赘述。
5.2.1.2.1
主燃烧室
主燃烧室的作用是通过燃油的燃烧使气流的温度升高。在假设理想的情况
下,主燃烧室的进口总温、总压、流量即为高压压气机的出口总温、总压以及
流量值。在主燃烧室进口总温
T
in
、进口总压
p
in
、空气流量
W
ain
以及出口温度
T
4
已知的情况下,进行如下过程,可搭建主燃烧室的计算模型。
1)出口油气比
流气
气流
气流
- 8 -
根据主燃烧室的出口温度
T
4
由焓函数:式
(10)
计算出口焓
h
out
,而进口焓
h
in
已知;因此根据式(16)可得到主燃烧室的出口油气比。
out in
b
b u in
f h h
H h
(16)
其中,燃烧效率
b
=0.99
,燃油热值
H
u
=42900000
。
2)计算燃油流量
f ain b
W W f
(17)
3)出口总压
* *
p p
out in b
(18)
其中,
b
=0.98
为主燃烧室总压恢复系数。
5.2.1.2 .2
涡轮
变循环发动机中有高压与低压涡轮两部分涡轮结构,两者使用相同的涡轮
模型。其中高压涡轮与高压压气机和 CDFS 同轴,其主要作用是使高温气流在其
中膨胀做功,以带动 CDFS 和高压压气机;而低压压气机与风扇同轴,类似高压
涡轮,在正常运转时使高温气流在其中膨胀,带动前段风扇。
在涡轮进口总温
T
in
,总压
p
in
、压比函数值
zz
、物理转速
n
、以及导叶角已知的情
况下,进行如下计算,构建涡轮模型。
1)求涡轮换算转速
*
,
*
in d
cor
in
T
n n
T
(19)
其中,高压涡轮:
*
,
T
in d
=1850
,低压涡轮:
*
,
T
in d
= 1.5405e+003
。
2
)线性插值法计算涡轮落压比、效率和换算流量
根据 1)中所得到的换算速度,以及附录 4 给出的涡轮特性,使用线性插值
法,利用 matlab 先根据换算速度,对落压比、流量和效率进行线性插值;再根
据一次插值后所得的增压比计算各压比函数值;最后根据给出的压比函数值,
对落压比、换算流量及效率进行第二次线性插值,可得到对应的落压比
,
pr
c map
、
效率
,
c map
和换算流量
,
W
c map
。
3)修正落压比、效率和换算流量
将 2)中线性插值法得到涡轮的落压比、效率和换算流量进行如下修正:
- 9 -
,
,
2
,
( 1)(1 ) 1
100
(1 )
100
(1 )
100
pr
c pr c map
w
c W c map
c c map
k
pr C pr
k
W C W
k
C
(20)
高压、低压涡轮的修正系数
pr
w
k k k
、 、
均分别取
1
,
1
,
0.01
;系数
C
的取
值为:
高压涡轮:
1.5342
13.2121
1.0121
W pc
C
C
C
,低压涡轮
0.7902
0.3881
1.0061
pr
W
C
C
C
。
4)计算涡轮流量
涡轮出口流量可由式(18)计算所得:
*
*
,
* *
,
in d
in
g
in in d
T
p
W W
T p
c
(21)
其中,高压涡轮:
*
,
p
in d
=28.7297
,低压涡轮:
*
,
p
in d
=11.3371
。
5)涡轮出口总压
*
*
in
out
c
p
p
pr
(22)
6)涡轮出口总温
由计算得的落压比
pr
c
、效率
c
以及已知的入口总温,可根据式
(20)
计算得
口总压:
/
*
*
1 1 /
Cp R
out
c
in
T
pr
T
c
(23)
其中,高压涡轮平方等压比热
CP
1.2988e+003
,低压涡轮平均等压比热
CP
1.2745e+003
,
R
为气体常数。
7)求涡轮的混合气体的进出口焓
涡轮的功和功率的计算过程主要为:首先根据进口总温
T
in
*
分别计算纯燃气
的进口焓
,
h
st in
(式
24
)以及此时空气的进口焓
,
h
air in
(式
10
),再计算在此油气
比下的燃气的总进口焓
h
in
(式 25);按照同样的步骤,根据出口总温
*
T
out
计算燃
- 10 -
气的总出口焓
h
out
。
6 3 2
3 5 4 9
5 12 6 16 7
0.11152575 10 0.31020206 10 2.9961197
0.27934788 2 0.18746407 10 0.73499597 10
0.15062602 10 0.12510984 10
Hst T T
e T T
T T T
(24)
1
f
Hgas Hair Hst
f
(25)
其中
f
为油气比。
8)求涡轮功和功率
T in out
T g T m
l h h
N W l
(26)
其中,
m
=0.99
为涡轮机械效率。
5.2.1.2.3
涵道
发动机中涵道的模型结构很简单,主要作用为允许部分气体流通。在涵道
进口总温
*
T
in
、总压
*
p
in
、流量
,
W
a in
、总压恢复系数
duct
已知的情况下,可由以
下过程直接获得涵道的出口总温、总压、流量。
* *
* *
, ,
out in
out in duct
a out a in
T T
p p
W W
(27)
其中
duct
0.98
为涵道的总压恢复系数。
5.2.1.2.4
涵道前混合器,选择活门,副外涵道建模
与常规混排涡扇发动机不同的是,双外涵变循环发动机在副外涵进口处,
安装有改变发动机工作模式的选择活门,模式选择活门的开闭,将改变变循环
发动机内部气动热力学循环。为了准确描述发动机在亚音速巡航点和超音速巡
航点的运转状态,建立了包含模式选择活门、副外涵及 CDFS 涵道的气流掺混
模型[1]。如图 4 所示。
- 11 -
图 4 变循环发动机局部简图
(
图中数字序号表示发动机各截面参数定义的下脚标
)
1) 求气动函数
125
q
( )
和
125
在副外涵道、
CDFS
涵道出口面积及外涵道流量其流量已知的情况下,
根据式(
27
)求得
125
q
( )
和
125
。
*
125
125 125 125
*
125
1 1
2
1 1
1
1
( )
1 1
( ) ( ) (1 )
2 1
2
1
g m
m
p
W k A q
T
q
k
R
(27)
其中,中
CDFS
涵道出口面积
A
125
= 608.4252
,
*
p
125
为
CDFS
涵道出口总压,
*
T
125
为
CDFS
涵道出口总温。
为气体绝热指数,
k
m
为流量系数,R 为气体常数。
纯空气:
=1.4
287
m
0.0404
R
k
,燃气:
=1.33
287.31
m
0.0394
R
k
。
2)计算 CDFS 涵道出口静压
p
125
*
p p
125 125 125
( )
(28)
其中
125
(
)
满足气动函数:
- 12 -
2 1
1
1
1
( )
(
)
(29)
3)求
225
(
)
和
225
由前混合器静压平衡可直接求出
225
( )
和
225
,
*
p
225
为副外涵道出口总压;
225 125
*
225 225 225
p p
p p
(
)
(30)
4)计算副外涵道出口流量
*
225
225 225 225
*
225
( )
g m
p
W k A q
T
(31)
其中副外涵道面积(选择活门面积):
A
225
=1.8395e+003
,
*
T
225
为副外涵出
口总温。
5)计算前混合器出口总温
*
T
15
、总压
*
p
15
、流量
W
g
15
225 125 15
225 225 125 125 15 15
* * *
225 225 225 125 125 125 15 15 225 125
15
( ) ( ) ( )( )
g g g
g g g
g
W W W
W h W h W h
p f A P f A p f A A
W
*
15
15 15
*
15
( )
m
p
k A q
T
(32)
其中:
15 225 125
A A A
,
h
15
为前混合器出口焓,
h
125
为前混合器
CDFS
涵
道出口焓,
h
225
为前混合器副外涵出口焓,
W
g
225
为前混合器副外涵出口流量,
W
g
125
为前混合器
CDFS
涵道出口流量。
*
T
15
是由式(
32
)的第二个公式求出的。
且
f
( )
满足气动函数:
1
-1
2
( ) ( ) ( )
1
f q z
(33)
5.2.1.2.5
后混合室
后混合室的主要作用是将各气流均匀混合,提高燃气的利用率。在已知混
合气两股参混前的总温、总压、流量及面积,即可根据以下建立的后混合室模
型计算出口总温、总压、流量等参数。
1)
计算
61
q
( )
和
61
*
61
61 61 61
*
61
( )
g m
p
W k A q
T
(34)
其中内涵出口面积
A
61
= 5.3061e+003
,
*
p
61
为内涵出口总压,
*
T
61
为内涵出口总温。
- 13 -
2)计算
62
q
( )
和
62
*
62
62 62 62
*
62
( )
g m
p
W k A q
T
(35)
其中外涵出口面积
A
62
= 2.3212e+004
,
*
p
62
为外涵出口总压,
*
T
62
为外涵出口总温。
3)计算内涵道及外涵道静压
*
61 61 61
*
62 62 62
( )
( )
p p
p p
(36)
4
)混合器出口总温
T
6 *
、总压
*
p
6
、流量
W
g
6
61 62 6
61 61 62 62 6 6
* * *
61 61 61 62 62 62 6 6 61 62
*
6
6 6 6
6 *
( ) ( ) ( )( )
( )
g g g
g g g
g m
W W W
W h W h W h
p f A P f A p f A A
p
W k A q
T
(37)
其中,
A A A
6 61 62
,
h
6
为后混合器出口焓,
h
61
为后混合器内涵出口焓,
h
62
为后混
合器外涵出口焓,
W
g
61
为后混合器内涵出口流量,
W
g
62
为后混合器外涵出口流量。
5.2.1.2.6
加力燃烧室
加力燃烧室的主要作用是将混合的燃气和空气进行燃烧,其构建的模型也
十分简单。当进口总温、总压、流量已知时,其出口总温、总压及流量可由式
(38)得到:
* *
6 6 6
* *
6 6
6 6
out in
out in
g out g in
P P
T T
W W
(38)
其中,
P
6 *
in
、
P
6 *
out
分别为进出口总压,
T
6 *
in
、
T
6 *
out
分别为进出口总温,
*
6
W
g in
、
*
6
W
g out
分别为进出口流量,
6
=1
为总压恢复系数。
5.2.1.2.7
喷尾管
喷尾管的作用是使燃气在其中膨胀,将燃气部分热能转化动能,从而产生
推力。本文采用拉瓦尔(收敛-扩张)喷尾管(如图 5 所示)建立模型。
- 14 -
7
8
9
图 5 拉瓦尔喷尾管示意图
在拉瓦尔尾喷管中,任意截面总温、总压、流量均不变,则由流量公式
可以得到:
( )
q A const
(38)
因此在已知任意截面的面积或者
,就可以求出该截面的参数。
拉瓦尔尾喷管有三种工作状态:临界、亚临界和超临界。当处于临界时,尾
喷管喉部
Ma
1
,喉部之后气流变为超音速气流,尾喷管出口静压与大气压相等
(完全膨胀);处于超临界时,喉部
Ma
1
,此时尾喷管出口面积会自动改变(增
大)使尾喷管出口静压与大气压相等,使尾喷管变为临界状态,但尾喷管出口面
积有最大限制,当到达最大限制值时,尾喷管出口静压不能与大气压相等,则通
过
( )
q A const
重新计算出口参数;处于亚临界时,喉部
Ma
1
,喉部之后不能
加速到超音速。
在已知尾喷管进口总温
*
T
in
、总压
*
p
in
、流量
W
g
、大气环境压力
p
0
的条件下,
建立拉瓦阿尔喷尾管模型。
1
) 计算尾喷管喉道面积
A
8
,出口面积
A
9
假设尾喷管始终处于临界或超临界状态,即速度系数
8
=1
。
由流量公式
8 *
8 8 8
8 *
( )
g m
p
W k A q
T
计算出
A
8
;
p p
9 0
,计算出
( ) /
9 9 9 *
p p
,并求出
9
;
由流量公式
9 *
9 9 9
9 *
( )
g m
p
W k A q
T
计算出
A
9
;
判断
9
8
A ar
A
(这里
ar
3
),如果是,则
9 8
3
A A
,利用流量公式重新计
算
9
。
2
) 计算尾喷管出口静温
*
T T
9 7 9
( )
;
(39)
- 15 -
其中
( )
为气动函数,可由式
(40)
求出:
2
( ) 1 1
1
(40)
3
) 计算尾喷管出口气流速度
*
9 9 9
nz
2
c h h
(41)
其中:
nz
0.98
,焓
h
9 *
和焓
h
9
分别由尾喷管出口总温和静温求出。
5.2.2
发动机平衡方程
假设发动机处于稳定运行在单/双涵道时,不考虑起飞、加速、爬升、单双
涵道切换时的动态情况,发动机各部件匹配工作时,受如下 7 个平衡方程制约。
1) 低压轴功率平衡
N N
CL TL mL
0
(42)
其中
N
CL
是风扇消耗功率,
N
TL
是低压涡轮发出功率,
mL
=0.99
为中间轴机
械效率。
2) 高压轴功率平衡
N N N
CH CDFS TH mH
0
(43)
N
CH
和
N
CDFS
分别是高压压气机和CDFS的消耗功率,
N
TH
是高压涡轮发出功
率,
mH
=0.99
是高速轴的机械效率。
3) 高压涡轮进口截面流量平衡
41 41
0
g g
W W
(44)
W
g
41
是高压涡轮进口截面气体流量,即主燃烧室出口气体流量和冷却空气流
量,
W
g
41
是通过高压涡轮特性数据线性插值得到的高压涡轮流量,这里忽略冷却
的空气流量。
4) 低压涡轮进口截面流量平衡
45 45
0
g g
W W
(45)
W
g
45
是低压涡轮进口截面流量,由主燃烧室出口气体流量和冷却空气流量计
算得到,
W
g
45
是通过低压涡轮特性数据线性插值得到的低压涡轮流量,这里忽略
冷却的空气流量。
- 16 -
5) 后混合器静压平衡
p p
61 62
0
(46)
p
61
和
p
62
分别是后混合器内、外涵道(主外涵道)的静压,二者应平衡。后
混合器见图2的6截面.
6) 尾喷管面积平衡
8 8
A A
0
(47)
'
A
8
为给定的尾喷管8截面的面积,这里
'
A
8
9.4575e+003
、
A
8
为按附录1尾喷
管的有关公式计算出的尾喷管8截面的面积,二者应相等。
7) 风扇出口流量平衡
2 21 13
0
a a a
W W W
(48)
其中风扇出口的流量
W
a
2
分流为副外涵流量
W
a
13
和CDFS进口流量
W
a
21
,三者之
间应存在平衡关系。副外涵流量
W
a
13
是按附录1中前混合器的有关公式计算出的。
5.2.3
模型的求解
发动各部件匹配工作时,受 5.2.3 中的七个平衡方程制约,若设定适当的
未知参量,各个平衡方程中的参数均可根据 5.1 及 5.2.1 中构建的发动机各部
件模型得到其关于已知量、假定参数或中间变量的表达方式。由此平衡方程组
转换成一组关于假定参数非线性方组,求解非线性方程组常用的方法有数值迭
代法和进化计算方法等。Newton-Raphson 迭代法虽然收敛较快,但是每步迭代
都要计算 Jacobi 矩阵并求解线性方程组,工作量较大。当 Jacobi 矩阵奇异或
接近奇异时,NR 法的计算将无法进行。Newton-Raphson 迭代法对初值的选择
比较敏感,不恰当的初值使得迭代难以收敛。在发动机的非线性方程组中,由
于不知道方程组最终解的范围,初值的选取具有随机性,难以将七个变量的初
值选在解的领域中,所以不宜使用 Newton-Raphson 迭代法。粒子群算法
(Particle Swarm Optimization,PSO)为一种基于随机性的非数值算法,是
一种新颖的、简便的智能算法,整个过程中毋须导数,同传统的优化方法相比,
PSO 具有更强的全局优化能力,也能较快地收敛于可接受解,对初值的选择比
较不敏感。所以,PSO 能够比较准确地求解发动机的非线性方程组。
5.2.3.1
假设参量的选定
由 5.1 及 5.2.1 中所构建的各部件的模型可知,若已知发动机飞行高度
H
、
飞行马赫数、低压转速、高压转速、各压气机和涡轮的压比函数值和导叶角、
主燃烧室的出口温度、各涵道与混合器的面积等参量时,即可获取发动机的所
有热动数据。
在假设理想的情况下,每一级部件的进口总温、总压等于上一级部件的出
口总温、总压。根据题目及附录总所给的数据,可知低压转速
n
L
、各压气机及
- 17 -
涡轮模型的导叶角、飞行高度、飞行马赫数、各涵道与混合器面积均是已知量,
则若假设高压转速、各压气机与涡轮的压比函数值以及主燃烧室的出口温度也
已知,便可根据部件法所构建的发动机模型与七大平衡方程得到发动机的稳态
运行过程。
选定的假设参数极其说明如下表 5 所示:
表 5 选定的七个假设参数说明
试给参数
说明
试给参数
说明
n
H
高压转速(高压压气机、
CDFS、高压涡轮物理转速)
T
4 *
主燃烧室出口温度
Z
CL
风扇压比函数值
Z
TH
高压涡轮压比函数值
Z
CDFS
CDFS
压比函数值
Z
TL
低压涡轮压比函数值
Z
CH
高压压气机压比函数值
5.2.3.2
粒子群算法求解非线性方程
粒子群优化算法(PSO)起源对简单社会系统的模拟,在对粒子群算法数学化
后,PSO 是一种很好的非线性方程的求解方式。
5.2.3.2.1
粒子群算法原理
在粒子群算法(
Particle Swarm Optimization, PSO
)中,每个优化问题的潜在
解都是搜索空间中的一只鸟,称为粒子。每个优化问题的潜在解都是搜索空间
中的一只鸟,称之为粒子。所有的粒子都有一个由被优化的函数决定的适值
( fitness value)
,每个粒子还有一个速度决定它们飞翔的方向和距离。粒子们就
追随当前的最优粒子在解空间中搜索最佳位置[6]。
假设在一个群体中有 N 个粒子,不同的个体有不同的位置
X
j
,它在控件中
以某个速度
V
j
飞行,不同的位置对应于不同的与优化目标函数值相关的个体适
应度函数值
( )
j
F X
。
POS 初始化为一群随机粒子(随机解),然后通过迭代找到最优解。在每一
次的迭代中,粒子通过跟踪两个极值来更新自己:第一个就是粒子本身所经历
的最好位置,为
1 2 7
[ , ]
T
j j j j
P p p p
…
,也称为
P
best
;第二个就是群体中所有粒子
经历的最好位置,为
1 2 7
[ , ]
T
g g g g
P p p p
…
,也成为
g
best
。若已知第 k 代第 j 各粒
子的速度
( )
k
V
j
及位置
( )
k
j
X
,则第(k+1)代第 j 各粒子的速度及位置为:
( 1) ( ) ( ) ( ) ( )
1 1 2 2
( ) ( )
k k k k k
j j j j j j g j
V w V c r P X c r P X
(49)
( 1) ( ) ( 1)
k k k
j j j
X X V
(50)
- 18 -
( )
max max min max
( ) /
k
w w k w w k
(51)
其中:
( )
k
w
:惯性权重系数,是迭代次数的函数,且随迭代次数线性减少;
w
max
:初始惯性权重;
w
min
:终止惯性权重;
k
max
:最大迭代次数。
c
1
:粒子自身加速度权重系数,一般在
0~2
之间取值;
c
2
:全局加速度权重系数,一般也在
0~2
之间取值;
1 2
,
r r
:
[0,1]
范围内两个相互独立的、均匀分布的随机数。
[4]
5.2.3.2.2
粒子群算法流程
基本粒子群算法的流程如下:
Step1
:初始化粒子群,包括群体规模
N
,每个粒子的位置
X
j
和速度
V
j
Step2
:计算每个粒子的适应度值
( )
j
F X
;
Step3
:对每个粒子,用它的适应度值
( )
j
F X
和个体极值
P
best
比较,如果
( )
j
F X
优于
P
best
,
( )
j
F X
替换掉
P
best
;
Step4
:对每个粒子,用它的适应度值
( )
j
F X
和全局极值
g
best
比较,如果
( )
j
F X
优于
g
best
则用
( )
j
F X
替
g
best
;
Step5
:根据公式(
49
)
-
(
51
)更新粒子的速度
V
j
和位置
X
j
;
Step6
:如果满足结束条件
(
误差足够好或到达最大循环次数
)
退出,否则返回
Step2
。
5.2.3.2.3
粒子群算法本题的实例化
在本文中,有七个假定的未知量
1 2 7
,
x x x
…
(对应关系为:x1:
n
H
高速物理
转速(高压压气机、CDFS、高压涡轮物理转速) x2:
Z
CL
风扇压比函数值; x3:
Z
CDFS
CDFS 压比函数;4;
Z
CH
高压压气机压比函数; x5:
T
4 *
主燃烧室出口温
度; x6:
Z
TH
高压涡轮压比函数; x7:
Z
TL
低压涡轮压比函数。)每个粒子在
某一时刻都有由这七个未知量组成的位置
1 2 7
[ , ]
T
j j j j
X x x x
…
。
- 19 -
设定函数
1
( )
j
f X
…
7
( )
j
f X
为平衡方程 1-8(式 42-式 48)在位置
X
j
时的计
算值。若
X
j
为非线性方程组的解时,满足平衡方程,即
1
( )
j
f X
…
7
( )
j
f X
均等于
0。令个体适应度函数为
( )
j
F X
:
1 2 7
( ) ( ) ( ) ( )
j j j j
F X f X f X f X
…
(49)
当
( ) 0
j
F X
时,位置
X
j
为最佳位置,即为非线性方程组的解。
在寻找最佳位置(即最优解)时,
( )
j
F X
的值越小,位置越佳。粒子根据
式(49)-(51)来更新自己的速度和位置。
对于本题的粒子群算法的流程图如图 6 所示:
- 20 -
开始
初始化每个粒子的速度和位置
计算每个粒子的适应度值
求出每个粒子的个体最优
求出整个群体的全局最优
粒子当前位置的适应度值小
于个体最优适应度值?
将当个体最优设为当前粒子位置
粒子当前位置的适应度值小
于全局最优适应度值?
根据式
49
对粒子的速度进行优化
根据式
50
对粒子的位置进行优化
适应度值对于
0
的误差够小?
达到最大迭代次数?
输出结果
结束
将当全局最优设为当前粒子位置
是
是
否
否
是
否
否
是
图
6
粒子群算法流程图
5.2.3.2.4
粒子群算法求解结果
根据部件法建立的模型,粒子群算法计算过程中,选定不同数目的群体个
数,经多次试验,发现在给定的不同的初值位置
X
j
情况下,都能有较好的收敛
效果,使得群体中各个粒子的适应度函数
1 2 7
( ) ( ) ( ) ( )
j j j j
F X f X f X f X
…
均能收敛到相近的数值,精度设为:eps=10^(-4),且该数值在允许的范围内可
近似认为等于 0。
如根据表 6 中各参数的设定,经过粒子群算法计算后,可得到一组比较理
- 21 -
想的非线性方程的解。
表 6 设定的初始参数列表
附加
参量
N
T
c
1
c
2
含义
群体粒子数
最大迭代次
数
粒子自身加速度
权重系数
全局加速度权
重系数
数值
30
100
2
2
初始
位置
n
H
T
4 *
Z
CL
Z
CDFS
含义
高压转速
主燃烧室出
口温度
风扇压比函数值
CDFS
压比函数
值
数值
0.6
1000
0.5
0.5
初始
位置
Z
CH
Z
TH
Z
TL
含义
高压压气机压
比函数值
高压涡轮压比函
数值
低压涡轮压比函数值
数值
0.5
0.5
0.5
在表 6 的给定的参数下,最终可得到一组最优解,如表 7 所示。粒子群算
法中选取的 30 组粒子在迭代过程中的趋势如图 7 所示。
表 7 线性方程的一组最优解
变量
n
H
T
4 *
Z
CL
Z
CDFS
含义
高压转速
主燃烧室出口
温度
风扇压比函数值
CDFS
压比函
数值
最优解
0.78
1369.9999
0.3394
1
变量
Z
CH
Z
TH
Z
TL
含义
高压压气机压
比函数值
高压涡轮压比函数
值
低压涡轮压比函数
值
最优解
0.2899
0.246
0.9112
由于群体粒子数 N=30,迭代一次需要的时间比较长。在计算时,为较好地
保存结果,按迭代次数分批次进行统计计算。图 7-图 9 分别为 PSO 算法三次迭
代的结果。
图 7 所示的结果中迭代了 20 次,适应值 F(x)从 1189208.9387 降低至
481109.0440。图 8 所示的结果中迭代了 20 次,适应值 F(x)从 481109.0440 降
低至 53218.4441。图 9 所示的结果中迭代了 20 次,适应值 F(x)从 53218.4441
降低至 9434.49741。此时 F(x)相对来说已经很小了,比较接近方程的解。由于
时间紧迫,未能进行更多次数的迭代。因此,以第三次迭代的结果作为一组最
优解。迭代过程如图 7 至图 9 所示。
图
7
第一次迭代结果
图
8
第二次迭代结果
- 22 -
- 23 -
图
9
第三次迭代结果
获取最优解的过程中的过程中在迭代过程中计算的各个模型的中间量结果
如表 8 所示。
表 8 中间变量计算结果
参数
T
out
P
out
W
out
l
含义
出口总温
出口总压
出口流量
功率
进气道
244.3812
0.3447
-
-
风扇
333.2381
0.9083
15.8659
1.41215*10
6
副外涵道
333.29
0.89014
1.9769
-
CDFS 涵道
368.4658
1.2273
3.0950
-
前混合器
355.3631
0.3168
5.0719
-
CDFS
369.4658
1.2524
9.7639
-
主外涵道
355.3631
0.31042
5.0719
-
高压压气机
527.75539
2.9973
6.6688
1.07850*10
6
主燃烧室
-
2.9373
6.8234
1.454718*10
6
高压涡轮
1124.9190
1.0766
4.8727
1.426426*10
6
低压涡轮
894.6882
0.3572
5.2770
-
- 24 -
加力燃烧室
894.6882
0.3572
5.2770
-
后混合器
637.8398
0.8096
10.3490
-
从上面数据可以看出表格 7 得出的一组最优解是存在一定误差的,由于时
间问题,未能继续求解,但是从图 7-图 9 的计算过程中可以看出,粒子群算法
对于发动机的非线性方程组求解是有效的,粒子个体逐步逼近方程组的最优解。
5.3
问题
3
的求解
问题三要求在发动机处于单涵道模式时,若使得发动机性能最优,发动机
CDFS
导叶角度、低压涡轮导叶角度和喷管喉道面积
3
个量的取值。这个问题
涉及到最优解的问题。由题可知,在最优解时,发动机有较强的动力性和燃油
经济性。并且,发动机的新能参数推力
F
,耗油率
sfc
,单位推力
Fs
可以用发
动机
CDFS
导叶角度、低压涡轮导叶角度和喷管喉道面积这三个变量表示,问
题转化成关于这三个变量的发动机性能参数的最优解问题。因此首先需要将性
能参数用这三个变量表示出,同
5.2
的求解过程类似,根据发动机的整体模型
和部件法构建发动机的各个部分计算模型。由于此时是超音速飞行,飞机需要
有较强的动力,所以此时飞机发动机工作在单涵道模式,部件法建模计算时,
同样受到平衡方程制约,在单涵道模式下,没有混合器流量制约,但是单涵道
模式中,可以认为,所有气体流量全部走内涵道,那么有前后级部件流量是相
等的。设定相应的未知量以及中间变量,构建未知参数的非线性方程组,最终
可以消除所有中间变量,将发动机性能参数用题目中的三个变量表示,然后再
求解函数的最优解。对于本题中函数最优解问题,仍然可以用粒子群算法。
5.3.1
模型的建立
5.3.1.1
单涵道模式下的整机模型
根据 5.2 中各部件的计算模型以及单涵道模式的运行过程,可建立如图
10 所示的在发动机模型在单涵道下的气流流动示意图。
CDFS
进
气
道
模
型
气流
压
气
机
模
型
气流 气流
压
气
机
模
型
风扇
压
气
机
模
型
高压压气机
主
燃
烧
室
模
型
气流 气流 气流 涡
轮
模
型
主燃烧室 高压涡轮
涡
轮
模
型
气流
低压涡轮
加
力
燃
烧
室
模
型
气流 气流 尾
喷
管
模
型
气流
加力燃烧室 尾喷管
燃油
燃油
模式转换活门关闭
图
8
双涵道变循环发动机单涵道工作时气流流动示意图
5.3.1.2
单涵道模式下各部件工作状态
根据
5.2
中建立的各部件计算模型以及单涵道模式下的各个部件的工作状
- 25 -
态,可进行如下模型的建立。
5.3.1.2.1
进气道模型,
根据公式(
2
)
~
(
5
),可以求出,进气道出口的总温
T
1
、总压
p
1
5.3.1.2.2
风扇模型
假设风扇转速
n
L
,风扇压比函数值
Z
CL
,已知风扇导叶角
L
,
进气道出口
的总温
T
1
、总压
p
1
(进气道模型求出),根据公式(
6
)
~
(
13
),可以求出风扇
出口总温
T
out
,1
,总压
*
p
out
,1
,出口流量
W
a
,1
,功率
N
c
,1
。
5.3.1.2.3
核心驱动风扇机
CDFS
模型
假设
CDFS
转速为
n
H
,
CDFS
压比函数值
Z
CDFS
,CDFS 导叶角
CDFS
,已
知入口总温
T
out
,1
,总压
*
p
out
,1
(风扇模型已经求出)。
根据公式(
6
)
~
(
13
),可以求出
CDFS
出口总温
T
out
,2
,
CDFS
出口总压
*
p
out
,2
,
出口流量
W
a
,2
,功率
N
c
,2
。
5.3.1.2.4
高压压气机模型
假设高压压气机压比函数值
Z
CH
,已知高压压气机转速
n
H
,高压压气机导
叶角
H
,(前面已经假设),高压压气机入口总温
T
out
,2
,入口总压
*
p
out
,2
(
CDFS
模型求出)。
根据公式(
6
)
~
(
13
),可以求出高压压气机出口总温
T
out
,3
,出口总压
*
p
out
,3
,
出口流量
W
a
,3
,功率
N
c
,3
。
5.3.1.2.5
主燃烧室模型
假设出口温度
T
4 *
,已知主燃烧室入口温度
T
out
,3
,入口总压
*
p
out
,3
,入口流量
W
a
,3
(高压压气机模型求出)。
根据公式(
16
)
~
(
18
)可以求出主燃烧室出口总压
*
p
out
,4
,出口流量
W
g
,4
,
出口油气比
f
b
。
- 26 -
5.3.1.2.6
高压涡轮模型
假设高压涡轮压比函数值
Z
TH
,已知高压转速
n
H
,高压涡轮导向角,入口
总温
T
4 *
、总压
*
p
out
,4
。
根据公式(
19
)
~
(
26
)可以求出高压涡轮出口总压
*
p
out
,5
,出口总温
*
T
out
,5
,
出口流量
W
g
,5
,功率
N
c
,5
。
5.3.1.2.7
低压涡轮模型
假设低压涡轮压比函数值
Z
TL
,低压涡轮导叶角
CH
,
已知 入口总温
*
T
out
,5
,
入口总压
*
p
out
,5
(高压涡轮求得)。
根据公式(
19
)
~
(
26
),可以求得低压涡轮出口总压
*
p
out
,6
,出口总温
*
T
out
,6
,
出口流量
W
g
,41
,功率
N
c
,6
。
5.3.1.2.8
后混合器模型
+
加力燃烧室模型
已知内涵道截面积,入口总温,总压,流量,根据公式(
34
)
~
(
38
)可以
求出出口流量
W
g
,7
,出口总温
*
T
out
,7
,
*
p
out
,7
。
5.3.1.2.9
尾喷管模型
假设
A8
截面积,已知入口总温,总压,流量。
根据公式(
39
)
~
(
41
),可以求出尾喷管出口气流速度
c
9
,尾喷管出口静温
T
9
,
出口截面积
A
9
。
5.3.2
平衡方程构建
与 5.2 中类似,也是适当选取一组位置参数,根据一些列平衡方程,通过
5.3.1 中建立的模型,可得到一组非线性方程组,再通过粒子群算法求解该线
性方程组。
5.3.2.1
平衡方程组
同样各个部件之间也受如下一些平衡方程影响和制约:
1) 低压轴功率平衡:
N N
CL TL mL
0
2) 高压轴功率平衡:
N N N
CH CDFS TH mH
0
- 27 -
3) 高压涡轮进口截面流量平衡:
41 41
0
g g
W W
4) 低压涡轮进口截面流量平衡:
45 45
0
g g
W W
5) 尾喷管面积平衡:
8 8
A A
0
6
) 各个部件之间由于是单涵道,所以还有流量平衡:
W
a
,1
=
W
a
,2
7
)
W
a
,2
=
W
a
,3
8
) 混合器面积公式
61 62 6
61 61 62 62 6 6
* * *
61 61 61 62 62 62 6 6 61 62
*
6
6 6 6
6 *
( ) ( ) ( )( )
( )
g g g
g g g
g m
W W W
W h W h W h
p f A P f A p f A A
p
W k A q
T
5.3.2.2
未知参量选取
根据各个部件的模型计算,可以得知假设变量有:风扇转速
n
L
,风扇压比
函数值
Z
CL
,
CDFS
转速为
n
H
,
CDFS
压比函数值
Z
CDFS
,CDFS 导叶角
CDFS
,
高压压气机压比函数值
Z
CH
,
出口温度
T
4 *
,高压涡轮压比函数值
Z
TH
,低压涡轮
压比函数值
Z
TL
,低压涡轮导叶角
CH
,尾喷管截面积
A
8
一共
11
个变量。
结合以上 8 个平衡方程和公式,完全可以消去变量定义中 8 个未知量。这
样可以用发动机
CDFS
导叶角度、低压涡轮导叶角度和喷管喉道面积
3
个量表
示出其他的变量。
那么对于发动机性能参数推力
F
,耗油率
sfc
,单位推力
Fs
便也随之可以
由这三个变量表示出。
选定的三个未知参量极其说明如下表
9
所示:
表 9 选定的 3 个假设参数说明
试给参数
说明
CD
发动机
CDFS
导叶角度
l
低压涡轮导叶角度
A
喷管喉道面积
- 28 -
5.3.3
单涵道模型的求解未解决的难点:
1
、在求解上述
8
个方程时,由于变量数多于方程数,所以必然有三个变量
无法解除,所以需要找出一种方法,可以快速变量求解出方程。
2
、假设已经找出三个变量和另外其他变量之间的关系后,还需要建立一个
性能评比函数,通过计算评比函数的最优值,如何选择这样的一个性能评比函
数,能够较可信的确定函数他的最值就是能够代表发动机的最优性能。需要进
行选择函数以及相应的计算。
3
假设已经确定性能评比函数,那么找到这个函数的最大值,就能确定
CDFS
导叶角度、低压涡轮导叶角度和喷管喉道面积三个量的最佳值,由于这是一个
三元的函数,求解一个三元函数的最大值问题可以用牛顿迭代法,遗传算法,
粒子群算法等等,鉴于前面分析各种方法优缺点,考虑用粒子群算法求解最优。
但鉴于时间有限,本文中只建立模型,以及分析求解算法,后续的求解过程和
结果仍在研究过程中。
六、模型评价与改进
6.1
模型的优缺点
6.1.1
部件法
本文采用部件法构建双涵道变循环发动机的模型。部件法的主要优势在于
将发动机的整个模型分成各个模块,各模块根据上级参数可独立分析,简化了
模型分析过程。
6.1.2
粒子群算法
1)同传统的优化方法相比 , PSO 具有更强的全局优化能力 ,也能较快地收敛
于可接受解;
2)PSO 采用实数编码,不需要像遗传算法一样采用二进制编码(或者采用针对
实数的遗传操作)。凡能用遗传算法求解的问题均能用粒子群算法, 不需要诸如
选择、交叉、变异和保留等遗传操作 ,因此控制参数少 ,算法简单效率较高[4];
3)PSO 也是基于群体的非数值并行算法 ,整个过程毋须导数信息;
4)对于多解的非线性方程组 ,可以通过独立计算多次来得到多个完全数值解;
6.2
模型的缺点
1)基本粒子群算法的性能依赖算法的参数选择;
2)由于 PSO 中粒子向自身历史最佳位置和邻域或群体历史最佳位置聚集,形成
粒子种群的快速趋同效应,容易出现陷入局部极值、早熟收敛或停滞现象;
3)基本粒子群算法容易陷入局部最优,计算量较大,求解时间比较长,如果参
数设置不但,求解速度比较慢。
6.3
改进方向
考虑到基本粒子群算法存在局限性,结合粒子群算法(PSO)和遗传算法(GA)
的优点,采用一种新的算法:遗传粒子群算法(GAPSO)。该算法首先计算适应度
函数值,对群体中的 体按适应度函数值的大小进行排序,再采用 PSO 算法对群
体中的优秀个体进行提高,提高以后的个体被保留并进入下一代,其他适应度
函数值差的个体则被淘汰。对已经提高的优秀个体再通过选择交叉和变异步骤
得到剩下的个体。采用新方法得到的下一代个体比采用遗传和粒子群方法得到
的下一代个体要优秀,因此 GAPSO 算法的性能有很大的提高[6]。
参考文献
[1]
苟学中
.
变循环发动机建模及控制规律研究
[D].
南京航空航天大学
, 2012.
[2]
苟学中
,
周文祥
,
黄金泉
.
变循环发动机部件级建模技术
[J].
[3]
陈长忆
,
叶永春
.
基于粒子群算法的非线性方程组求解
[J].
计算机应用与
软件
, 2006, 23(5): 137-139.
[4]
陈长忆
,
叶永春
.
基于粒子群算法的非线性方程组求解
[J][J].
计算机应用
与软件
, 2006, 23(5): 137-139.
[5] Schneider F. B., Andrews G. R., Concepts for Concurrent Programming [J],
LNCS224, Springer2Verlag, 1985. 669
~
716
[6] Kennedy J, Eberhart R
.
Particle swarm optimization [A]
.
in: Proceedings of the
4th IEEE International Conference on Neural Networks [C]
,
Piscataway: IEEE
Service Center, 1995, pp.1942 -1948
.
[6]
彭晓波
,
桂卫华
,
黄志武
,
等
. GAPSO:
一种高效的遗传粒子混合算法及其应
用
[J].
系统仿真学报
, 2008, 20(18): 5025-5027.
- 29 -