TA的每日心情 | 奋斗 2024-7-1 22:21 |
|---|
签到天数: 2014 天 [LV.Master]伴坛终老
- 自我介绍
- 数学中国站长
 群组: 数学建模培训课堂1 群组: 数学中国美赛辅助报名 群组: Matlab讨论组 群组: 2013认证赛A题讨论群组 群组: 2013认证赛C题讨论群组 |
摘 要:
- E& s8 ^! j' ^5 q本文研究的是变循环发动机部件法建模及优化问题。* g$ b' \, Z5 T! D/ J+ l
对于第一问,首先要求分析风扇特性,通过部件法建模,风扇属于压气机: a+ M, ]0 ]( V& j; H0 m9 F3 a
模型,根据压气机压比函数的定义,利用matlab 线性插值将附录4 中增压比转
2 Z8 ]3 Q/ L* C3 }换成压比值,作出在不同转换速度下的流量随压比函数值的曲线图,如图2 所6 A1 ~6 l) x2 l9 b( b
示。为计算风扇和CDFS 出口的总温、总压和流量,根据附录1 的计算公式,构
9 Z# x3 g5 Q: B' z q! D建进气道模型和风扇、CDFS 的压气机模型。根据已知数据,在亚音速状态下代4 g5 _, F6 V# S& n% |
入构建的进气道、风扇和CDFS 计算模型即可得到所求变量值:风扇出口总压为
/ g% J! H8 S5 G+ y& Q% l4 o% h1.2883、总温为378.3326、流量分别为19.0477;CDFS 的出口总温为419.3767,& |( O- ?" U1 R# ?" ^
总压为1.7745,出口流量为16.9401。其余中间变量见表1 至表3 所示。& J- M5 _ {$ e
第二问,要求解7 个平衡方程组成的非线性方程。首先,分析发动机的整# ?$ K( G' Q6 p
机模型,根据部件法构建发动机的各部分计算模型。在发动机各部件匹配工作
2 ?& {- Q- l) O' t4 J' @时,受7 个平衡方程制约。根据已知的参数和发动机模型,设定七个适当的未
K+ [5 _! h* u# @+ W2 w知数,分别为:高压转速" w e# H7 \6 W0 L* f+ B1 g+ K
H n 、风扇、CDFS、高压涡轮、低压涡轮的压比函数值:6 N# |' }! l% t4 b) d3 B8 n
CL Z 、
6 ~+ N* D5 q2 v0 Y# p- ^, _& |CDFS Z 、
6 T' P! P) D: C& Z: G% F4 g! PCH Z 、TH Z 、
, |3 i) B+ ~9 s' E1 N& t1 CTL Z 以及主燃烧室的出口温度*6 U2 v" p$ Y2 A A6 M# X; J
4 T 。由构建的发动机模
4 X. w% }5 K) l& _ d& `8 L型与其工作时各部件匹配的7 个平衡方程可得到一组非线性方程组。为求解此% A S& `# C" u/ \2 i0 |% p
方程组,采用粒子群算法。将非线性方程组的潜在解视为粒子,设定的未知量可
8 r6 U% X# Q" d6 r看作粒子的位置,粒子群中的粒子个体根据适应性函数引导进行迭代,向非线9 ], o3 c6 H" ~9 z/ r
性方程组的最优解不断逼近 ,通过迭代法,更新粒子的位置,以搜索到全局极
/ j3 B! G- [5 F8 r; P3 Q值,即获取非线性方程组的最优解。经过多次试验,粒子群算法均能搜寻到能
( k/ C7 I& a% Y得到比较准确的非线性方程的解。在给定不同的群体数和初值时,可得到以下 ]2 O* F2 U7 ]3 k% I* C* [
比较理想的解:+ N$ l8 m# q4 A/ O; B
- 2 -. s! L3 G4 G0 f
变量
# A8 _( s6 \2 \8 e& [0 b! ^2 x" \: CH n *( x* W; w2 V" G7 _/ H5 s& c
4 T
' D, O, ~3 Y2 o2 q3 m- U. s: ZCL Z
9 K7 l- J. [* k( ?/ R+ DCDFS Z. C) }( @4 z! [7 P
含义 高压转速0 _6 [. |9 G9 h( q" j, N$ z) m
主燃烧室出口
5 j! F( W" U( m3 N温度5 p/ Y# U+ L0 i6 m
风扇压比函数值% r2 U! Z, L6 J4 r
CDFS压比函7 I$ z/ U' @* J$ a. @ l! G
数值4 ]: l- G U. c" J/ C; D
最优解 0.78 1369.9999 0.3394 1
( U# q' r) ~- n变量
% A: l% K8 I, N. KCH Z TH Z
% ~3 l) _& M ~, \( p XTL Z" c- q2 [0 p7 y/ Y3 p
含义
J+ {+ O! I/ D+ x. z+ k高压压气机压8 m! p8 H T+ a& n# s
比函数值% L2 b6 m& V! [/ Y& W
高压涡轮压比函数
+ Y4 e( Y3 u) S# I* L& _值7 x' d; d8 R z( M0 n; \& b$ r+ @
低压涡轮压比函数# s9 i/ s) @7 |* O
值 G/ r0 |$ ]1 L6 ]
最优解 0.2899 0.246 0.9112
" z! F; B! I6 l# \; s, a* C对于第三问,求解发动机性能最优时,发动机CDFS导叶角度、低压涡轮导叶( a2 _9 {3 U+ Q H. Y
角度和喷管喉道面积3个量的取值。本题的求解为单涵道模式下最优解的查找问
( Q# {7 y6 y$ X题,可按照与问题二类似的步骤,首先根据发动机的整体模型和部件法构建发动
% o- t0 F9 M. Y- `& ?/ j机的各个部分计算模型;在单涵道工作模式下,在发动机各部件匹配工作时,受4 e! ]" K- @5 Y) \/ B
8个平衡方程制约。设定相应的未知量,建立关于发动机CDFS导叶角度CD 、低) w( v5 y5 S) ]$ _
压涡轮导叶角度l
) _( R6 _9 \9 e ]4 ~5 J以及喷管喉道面积A 这三个位置参量的非线性方程组,求出
) R( \/ Z+ X# M/ r, M7 U- i其他未知量与CD , l ' D) ~- i3 D! `) F W3 V. q
, A 的关系。进而表示出发动机性能参数推力F,耗油率
' B) U9 ~5 c( @* v# R3 }5 ]4 K9 {sfc,单位推力Fs的函数关系,另外根据性能最优要求,建立性能评价函数,得到
$ `" ^ O' ^: s% o* _# O7 a关于CD , l ) J8 A9 F; x+ n6 d% S. Q
, A 的三元函数,通过粒子群算法,求解函数的极值,便是性能最
& j- y8 P7 h- B* l! \3 B$ o- Z$ d优时的参数设定。但鉴于时间有限,本文中只建立模型,分析了求解算法,后续
, X: _$ z! T3 B: @0 C! ^的求解过程和结果仍在研究过程中。
/ O8 C; g0 B1 c2 T- g- Z2 U: J5 f- @( s. ^, P( R/ ]8 e
|
zan
|