TA的每日心情 | 奋斗 2024-7-1 22:21 |
|---|
签到天数: 2014 天 [LV.Master]伴坛终老
- 自我介绍
- 数学中国站长
 群组: 数学建模培训课堂1 群组: 数学中国美赛辅助报名 群组: Matlab讨论组 群组: 2013认证赛A题讨论群组 群组: 2013认证赛C题讨论群组 |
摘 要:7 D! E; u4 d! H' L5 F! u
本文采用部件级建模法精确模拟发动机的各个部件,依据各部件匹配工作时) m% [6 J5 C0 C; ~- V0 r6 T
的7 个平衡方程,对发动机的性能进行模拟。
3 b% k4 l% g; Y# g针对问题一,为了求解风扇和CDFS 的出口总温、总压和流量,建立模型对1 t- [, L+ Y( k# G' c' @6 a, K* _6 M
这两个部件的特性进行精确模拟,利用给定的发动机飞行高度和飞行马赫数,求% y7 h R* c4 g$ g$ i& G1 Z5 T
解出风扇的出口总温、总压、流量分别为379.4985、1.3087、19.0483,CDFS 的$ ]( \& X2 h/ q6 g* J! x, ~
出口总温、总压、流量分别为420.5365、1.8012、17.164。分析得出,气流在进
1 f' [9 h. o, k" X入风扇和CDFS 两个压气机部件至流出过程中,总温、总压增大,而气体从风扇
1 ?* N1 [8 }/ G3 R4 b! ]6 `流入到CDFS 的过程中,总温、总压亦增大,流量减小。此结论符合压气机压缩" p& Q9 x, [" p/ A
气体导致温度升高、压强增大、流量减小的功能特点。; ~' n7 Y3 w1 l2 l* c0 x8 s) O% \
针对问题二,根据发动机整机模型,由七个参数值可计算出平衡残差量。以% X' Q3 K0 I; e% S$ a, K+ Q
平衡残差量最小为原则,对离散化的待估参数进行变域、变步长的搜索,根据当' f+ Q4 Q' n; T/ I* [: q
前的最优解与次优解确定下一步的搜索域与搜索步长,逐步缩小搜索范围、减小7 W: D/ j. A& s6 |# D4 C! r
搜索步长,搜索的终止条件设为:(1)高压转速、压比函数值的搜索步长减小至
9 L) q3 S3 E# X+ N: s0.01,主燃烧室出口温度的搜索步长减小至10;(2)最优解与次优解相同。搜
& D, [* s& n, t4 F: t+ \索的终止条件保证了解的精度与收敛性。依此算法搜索得到高压转速、压比函数
* ~+ d8 W3 F0 ~2 F: t) K值(风扇、CDFS、高压压气机、高压涡轮、低压涡轮)、主燃烧室出口温度的最
8 x2 E3 f u: U0 }0 Z优解分别为1.00,0.33,0.43,0.53,0.14,0.12,1520,此时平衡方程残差量为
" ]& H0 r9 K* Y* C( Q% B0.2550。逐步搜索过程中参数的解与平衡方程的残差趋于固定值,参数的解为模
" C- ]/ o K2 P5 g9 b: D9 H2 t型的收敛解。" C; X; w) Z/ m6 W C6 u, A1 l; {
针对问题三(1),为了保证发动机性能最优,求解CDFS 导叶角度、低压涡- X4 H) x4 ?: B
轮导叶角度和喷管喉道面积3 个变量,实质上是一个优化的问题。本文建立优化3 b) A! B( y4 A
模型, 采用单位推力和耗油率的线性组合构建一个新的性能评价指标
& {% | J+ a- {% g1 26 v" V8 V0 I* A, d, O. c: |
A Fˆs sfˆc。( 1
( V+ Z+ X, x! I( X' I7 o F. b、2 为比例系数),以其最小值作为目标函数,同时借鉴5 M$ b* C* w$ D- j$ J1 h7 A
2
& f- n& }% g2 z; {. v问题二中求解非线性方程组的方法,利用参数遍历法对模型进行解算。最终得到; i% y: ?) ~" X* L2 ~ y/ w, x( B* t
CDFS 导叶角度、低压涡轮导叶角度和喷管喉道面积3 个量分别为3,15,6952.496
! Y2 m3 H; L: c* q# e+ |6 I: u* m时,发动机的性能最优,此时单位推力和耗油率分别为1293.092,0.000239 。同
/ k5 B3 c% c5 j1 Y3 n5 s" D时,通过对遍历过程中部分参数对应的发动机的性能大小,分析得出规律:低压
# ]) B2 j3 R+ X& ~5 y7 ?转速对发动机的性能无太大影响,提高主燃烧室的出口温度可以有效降低耗油率/ P6 [( f1 v$ ]+ d) A u
,增大风扇的压比函数值则能有效地增强单位推动力、降低耗油率。 d7 h# |5 u" k3 X j
针对问题三(2),探索CDFS 导叶角、低压涡轮导叶角和尾喷管喉部面积在# k8 E; N* P9 h, b/ ?2 `) n
发动机性能最优条件下随飞行马赫数的变化规律,基于工作点的变步长的搜索方* h( t: ?7 n% y, F$ f, @9 U7 D8 \
法,以发动机性能局部最优作为约束条件,以马赫数、CDFS 导叶角、低压涡轮# L7 {* R8 y: R c
导叶角为输入值,以尾喷管喉部面积、局部最优时对应的马赫数、CDFS 导叶角; R2 `6 Y- c! x9 y7 ^- n9 f* B3 H
、低压涡轮导叶角为输出值,建立了变步长最优化模型。得出的结果显示,在某
1 L8 c2 W* s8 [* m1 R; [4 j) Z个具体工作点时发动机性能最优的条件下,当马赫数增加时,CDFS 导叶角、低% |2 L+ N, q, v" v4 O& u) V0 o
压涡轮导叶角为恒定值,相关系数为0;而尾喷管喉部面积随马赫数的增大呈现
7 f3 M! B2 t0 u8 O, g阶梯性递减的情况,当马赫数增加到某个具体的值时,面积保持恒定。这与整机 r+ g- _( |7 w' w
模型中尾喷管喉部面积的规律描述相符。在本文给出的工作点1 下,压比函数值6 X5 u7 i6 [7 m. ~9 E
处在中位,CDFS 导叶角的值恒为35,低压涡轮导叶角恒为15,尾喷管喉部面积
6 i" S0 P& H* c; e/ F2 R" W z从4109.696 递减到4087.818 后保持恒定;工作点2 下,压比函数值处在高位,
0 X' O# T5 M. _4 X5 ~CDFS 导叶角的值恒为29,低压涡轮导叶角恒为15,尾喷管喉部面积从3336.678/ I/ |% l9 @) m$ E& q
递减到3283.023 后保持恒定;工作点3 下,压比函数值处在低位,CDFS 导叶角
* H% m% ^2 q: H4 k9 a4 @的值恒为29,低压涡轮导叶角恒为14,尾喷管喉部面积不变,为3369.63。0 C4 y7 Q3 [) p/ Y7 N z l
关键词:变循环发动机;部件法建模;平衡方程;变域变步长搜索
7 K1 p7 t) F: t" ~
3 V5 ]) w/ Z9 U6 \) ~/ h |
zan
|