TA的每日心情 | 奋斗 2024-7-1 22:21 |
|---|
签到天数: 2014 天 [LV.Master]伴坛终老
- 自我介绍
- 数学中国站长
 群组: 数学建模培训课堂1 群组: 数学中国美赛辅助报名 群组: Matlab讨论组 群组: 2013认证赛A题讨论群组 群组: 2013认证赛C题讨论群组 |
摘 要:
7 O8 `: p, z/ @# I% x9 j: A5 I$ T" I本文采用部件级建模法精确模拟发动机的各个部件,依据各部件匹配工作时
2 t/ q- F7 E6 s0 x的7 个平衡方程,对发动机的性能进行模拟。) f2 E9 _+ T% E/ T9 v9 t
针对问题一,为了求解风扇和CDFS 的出口总温、总压和流量,建立模型对2 c- }+ q9 s2 e) }
这两个部件的特性进行精确模拟,利用给定的发动机飞行高度和飞行马赫数,求. R( d' c. U. ^$ h) V
解出风扇的出口总温、总压、流量分别为379.4985、1.3087、19.0483,CDFS 的/ ]- x; v4 A F6 b, V3 V2 e
出口总温、总压、流量分别为420.5365、1.8012、17.164。分析得出,气流在进5 r* Y/ Y6 c' ?
入风扇和CDFS 两个压气机部件至流出过程中,总温、总压增大,而气体从风扇
/ f/ r2 {5 j( T4 N7 f, o* _0 {流入到CDFS 的过程中,总温、总压亦增大,流量减小。此结论符合压气机压缩! _3 U0 X6 q& P6 i2 ?7 F, l9 o
气体导致温度升高、压强增大、流量减小的功能特点。
1 F8 ?! W5 i$ w针对问题二,根据发动机整机模型,由七个参数值可计算出平衡残差量。以8 l3 u8 {, U: Y5 z: h, ^' Q, o
平衡残差量最小为原则,对离散化的待估参数进行变域、变步长的搜索,根据当
% Y1 M, H% e& m, i7 B2 R2 ~前的最优解与次优解确定下一步的搜索域与搜索步长,逐步缩小搜索范围、减小
. p" G9 P$ W' t( [搜索步长,搜索的终止条件设为:(1)高压转速、压比函数值的搜索步长减小至9 } k* t: n: D" K& F6 U
0.01,主燃烧室出口温度的搜索步长减小至10;(2)最优解与次优解相同。搜
, o% J+ F& n6 }4 p, z* ?, { h索的终止条件保证了解的精度与收敛性。依此算法搜索得到高压转速、压比函数4 A6 c2 l6 ]+ L, `! p& B- D
值(风扇、CDFS、高压压气机、高压涡轮、低压涡轮)、主燃烧室出口温度的最. f& [ q ~/ r7 I& F, ?
优解分别为1.00,0.33,0.43,0.53,0.14,0.12,1520,此时平衡方程残差量为- P/ ?! K9 E% }1 `
0.2550。逐步搜索过程中参数的解与平衡方程的残差趋于固定值,参数的解为模
( ~9 B! t0 E5 c/ {: T% }! s型的收敛解。
7 {, r6 D# q( B) {+ P8 L* t针对问题三(1),为了保证发动机性能最优,求解CDFS 导叶角度、低压涡
; |% N, e- C# u0 Q轮导叶角度和喷管喉道面积3 个变量,实质上是一个优化的问题。本文建立优化6 [% T% V3 D6 ?' r
模型, 采用单位推力和耗油率的线性组合构建一个新的性能评价指标
- t1 Z& P' @/ c, s1 2
: Q" Y, e: N9 c4 f- FA Fˆs sfˆc。( 1 & C0 a+ K$ B& q+ T' c0 k. O
、2 为比例系数),以其最小值作为目标函数,同时借鉴9 K! ]8 B& B- W6 v
2
T; o8 r( I4 `' Q* D问题二中求解非线性方程组的方法,利用参数遍历法对模型进行解算。最终得到
: t4 J- x0 Q0 g( r6 OCDFS 导叶角度、低压涡轮导叶角度和喷管喉道面积3 个量分别为3,15,6952.496& H0 \/ q" C- G2 i- n
时,发动机的性能最优,此时单位推力和耗油率分别为1293.092,0.000239 。同
$ ]' @* E5 v2 R2 C% |/ A+ v时,通过对遍历过程中部分参数对应的发动机的性能大小,分析得出规律:低压
+ D" q& ^3 E, N转速对发动机的性能无太大影响,提高主燃烧室的出口温度可以有效降低耗油率
2 d. a8 O6 i0 D. U% H,增大风扇的压比函数值则能有效地增强单位推动力、降低耗油率。
2 F, S0 M9 C/ p1 r* e" }针对问题三(2),探索CDFS 导叶角、低压涡轮导叶角和尾喷管喉部面积在 V( s" t6 ]* z' m
发动机性能最优条件下随飞行马赫数的变化规律,基于工作点的变步长的搜索方1 m8 V i! r F0 F
法,以发动机性能局部最优作为约束条件,以马赫数、CDFS 导叶角、低压涡轮. g* {$ U) Q1 @$ {% S8 ^
导叶角为输入值,以尾喷管喉部面积、局部最优时对应的马赫数、CDFS 导叶角. x( b* C' o" I- a, ]. Y# }
、低压涡轮导叶角为输出值,建立了变步长最优化模型。得出的结果显示,在某
l' S6 F, u3 H# ^个具体工作点时发动机性能最优的条件下,当马赫数增加时,CDFS 导叶角、低
3 T0 q+ b3 `% n& |5 X压涡轮导叶角为恒定值,相关系数为0;而尾喷管喉部面积随马赫数的增大呈现9 k: J: ~- l8 [6 O; k
阶梯性递减的情况,当马赫数增加到某个具体的值时,面积保持恒定。这与整机5 ^) B" a v7 n( N2 r( f' \
模型中尾喷管喉部面积的规律描述相符。在本文给出的工作点1 下,压比函数值
( ?& w- s3 }5 M$ F; a处在中位,CDFS 导叶角的值恒为35,低压涡轮导叶角恒为15,尾喷管喉部面积
7 l; n9 q: M/ v( y( S7 I T从4109.696 递减到4087.818 后保持恒定;工作点2 下,压比函数值处在高位,% B# X" j3 h8 ]& N+ q# D
CDFS 导叶角的值恒为29,低压涡轮导叶角恒为15,尾喷管喉部面积从3336.678! T$ C% D5 V3 y7 e+ H9 W
递减到3283.023 后保持恒定;工作点3 下,压比函数值处在低位,CDFS 导叶角9 `" M; N* r2 u, `8 G
的值恒为29,低压涡轮导叶角恒为14,尾喷管喉部面积不变,为3369.63。' k2 o: b7 X: b6 w* E
关键词:变循环发动机;部件法建模;平衡方程;变域变步长搜索, K+ l5 {, ~/ c' z& |
- L/ Q( e( B5 d! h# s2 b |
zan
|