TA的每日心情 | 奋斗 2024-7-1 22:21 |
---|
签到天数: 2014 天 [LV.Master]伴坛终老
- 自我介绍
- 数学中国站长
 群组: 数学建模培训课堂1 群组: 数学中国美赛辅助报名 群组: Matlab讨论组 群组: 2013认证赛A题讨论群组 群组: 2013认证赛C题讨论群组 |
摘 要:
4 h9 T1 [, }1 E9 ?/ z本文研究的是变循环发动机部件法建模及优化问题。
8 e- b$ a7 Z6 m' M) m2 G* z对于第一问,首先要求分析风扇特性,通过部件法建模,风扇属于压气机8 c! F" E- j# _* h; ]9 a$ _
模型,根据压气机压比函数的定义,利用matlab 线性插值将附录4 中增压比转9 b' |9 @( l: j3 B
换成压比值,作出在不同转换速度下的流量随压比函数值的曲线图,如图2 所
4 P9 O6 ^* d/ l! z. ~% Z示。为计算风扇和CDFS 出口的总温、总压和流量,根据附录1 的计算公式,构' a- O0 N G% ]1 w* O
建进气道模型和风扇、CDFS 的压气机模型。根据已知数据,在亚音速状态下代
1 L! d& a8 k& E) u入构建的进气道、风扇和CDFS 计算模型即可得到所求变量值:风扇出口总压为/ C5 h/ i- ^0 z( Z7 B5 m! u
1.2883、总温为378.3326、流量分别为19.0477;CDFS 的出口总温为419.3767," e$ e9 I* l9 Q- W) C
总压为1.7745,出口流量为16.9401。其余中间变量见表1 至表3 所示。
* Z- J4 T7 o1 k, x第二问,要求解7 个平衡方程组成的非线性方程。首先,分析发动机的整
$ s; \0 N0 K+ j: g机模型,根据部件法构建发动机的各部分计算模型。在发动机各部件匹配工作2 F0 _$ Y4 r' }+ ?3 T& H& \0 h3 B* o
时,受7 个平衡方程制约。根据已知的参数和发动机模型,设定七个适当的未6 I1 J% K. d1 @' G: b5 R
知数,分别为:高压转速
8 K) `, p( Z7 W. t) e KH n 、风扇、CDFS、高压涡轮、低压涡轮的压比函数值:/ P b2 u9 X1 z9 L' q; ^, D, X/ E
CL Z 、/ I& |6 i5 O g6 w. y
CDFS Z 、
4 a3 k1 S$ c+ n5 ^$ t, M# }CH Z 、TH Z 、
1 d/ b8 \: R Y7 uTL Z 以及主燃烧室的出口温度*
/ f1 S/ U) R f5 x( V B( |4 T 。由构建的发动机模
- `8 A# Q& K3 I& D/ k# Z- }型与其工作时各部件匹配的7 个平衡方程可得到一组非线性方程组。为求解此/ a- u9 {; c* X/ v: i
方程组,采用粒子群算法。将非线性方程组的潜在解视为粒子,设定的未知量可
) u: F- n* o8 E6 i" d* t6 ~看作粒子的位置,粒子群中的粒子个体根据适应性函数引导进行迭代,向非线3 |' R) L/ W$ E: X
性方程组的最优解不断逼近 ,通过迭代法,更新粒子的位置,以搜索到全局极
2 K4 E/ y* w' D) S5 g! d4 L7 W值,即获取非线性方程组的最优解。经过多次试验,粒子群算法均能搜寻到能
" C% k. U8 \# {3 l1 l4 Y6 T% D' m得到比较准确的非线性方程的解。在给定不同的群体数和初值时,可得到以下
- G# W$ K1 e' L/ Y0 P/ Z比较理想的解:
* K$ ?; d: z9 W$ T- 2 -& S! L9 J" O+ A9 T- B) U7 Y7 m2 f
变量
8 \+ n; \, b! E5 M% a' nH n *
0 q# j$ g) [) v. Z( ?4 T
% m& c. Y: v# }2 ?5 U% Z: d4 B% ?CL Z6 d. l: X( P/ l1 b; L
CDFS Z
0 B& |4 Y4 Z) G5 {7 O8 s含义 高压转速
! `" f, P% g( s- z主燃烧室出口
2 j$ ~/ r; b. u+ p2 Y7 _$ o; ]温度- {8 i1 O% O" @' d5 l
风扇压比函数值" `8 Q8 x& @7 D4 n d
CDFS压比函
* \0 S- ]" M, w6 o8 a6 {数值
7 G+ v5 O/ `9 v5 T最优解 0.78 1369.9999 0.3394 1. u5 B j* k0 H% [, D% V( K
变量+ ]- D1 `" t9 B0 v. g4 `
CH Z TH Z$ v1 f/ v* \8 N5 }) ~. z
TL Z
5 G: W! ~# w* L含义) ^# Z: F! C1 V! R( B8 J
高压压气机压
; `* Z+ _3 k& s/ E比函数值
1 P; k J' z+ [7 M7 Y高压涡轮压比函数
' @3 |. H0 ^" d6 O值0 t6 ~8 n- ^, x7 c
低压涡轮压比函数) B9 |4 K$ x6 {/ w' u
值! w( j# {) s- P5 m+ W M2 h
最优解 0.2899 0.246 0.9112; S4 H: S9 m! a0 ]- s3 O
对于第三问,求解发动机性能最优时,发动机CDFS导叶角度、低压涡轮导叶3 ]' v1 H% d; h
角度和喷管喉道面积3个量的取值。本题的求解为单涵道模式下最优解的查找问
. [6 I; l. a1 H+ I+ o! W X" Z2 O题,可按照与问题二类似的步骤,首先根据发动机的整体模型和部件法构建发动
" W/ c$ Y. `# }3 u& F机的各个部分计算模型;在单涵道工作模式下,在发动机各部件匹配工作时,受
6 K) I" c% G) b% g8个平衡方程制约。设定相应的未知量,建立关于发动机CDFS导叶角度CD 、低 x4 F/ y4 }+ W0 [1 O( [* `% \
压涡轮导叶角度l % @; B T( G/ T1 E, H
以及喷管喉道面积A 这三个位置参量的非线性方程组,求出
9 z! [+ ~" [: n, N9 c( `其他未知量与CD , l
# F- x4 v' h3 V8 ~( q" ~, l1 h, A 的关系。进而表示出发动机性能参数推力F,耗油率% F3 L: N- X/ E* m$ f; w& G
sfc,单位推力Fs的函数关系,另外根据性能最优要求,建立性能评价函数,得到
' \" e$ z- {* C& g4 J关于CD , l
5 I% d( e" G4 }! ^" O, A 的三元函数,通过粒子群算法,求解函数的极值,便是性能最& |( _; L% z3 I! M2 f$ ?8 m
优时的参数设定。但鉴于时间有限,本文中只建立模型,分析了求解算法,后续+ d |' Z1 G, W) |
的求解过程和结果仍在研究过程中。
, J1 c; v9 q4 w( D$ y) K8 m% w' {& K
|
zan
|