TA的每日心情 | 奋斗 2024-7-1 22:21 |
|---|
签到天数: 2014 天 [LV.Master]伴坛终老
- 自我介绍
- 数学中国站长
 群组: 数学建模培训课堂1 群组: 数学中国美赛辅助报名 群组: Matlab讨论组 群组: 2013认证赛A题讨论群组 群组: 2013认证赛C题讨论群组 |
摘 要:% c; b& }! E! J# ]* p
本文研究了题目中的所有问题,对问题一、二中的双涵道(涡扇)模式和
F' b3 Y. r5 ?- c% ?问题三中的单涵道(涡喷)模式分别建立了相应模式下的部件模型和整机气动" p5 v- k' ^ S) U D
热力学模型,利用遗传算法能解决高复杂性、多变量非线性方程问题的优点,
& L) Z+ K. R5 D8 ~+ r) o本文基于遗传算法对各个模型进行求解与优化,并最终获得了相应结果。
' t ~# W \! V: p# `' g$ J" E针对问题一:) T; q. g, [- T# p ~
首先,作出风扇的特性曲线图;对图进行分析可知:对于某一特定换算转速3 M3 ~, A( v" h# n- X0 q; o
下,当压比函数值不断减小时,流量逐渐趋近平稳;当压比函数值一定时,换算6 Y# h$ y0 L+ w' H' B& F
转速越大流量就越大,当压比函数值增大到一定程度后流量会减小。4 V7 o5 H6 P8 M1 _" W' K( N
其次,分别对进气道、风扇和CDFS 进行建模,运用二次线性插值获得对应( M/ i8 C Y- k7 A2 ^
的流量和效率代入模型,用牛顿迭代法反解出口总温;结果如下表1 所示:
. K( I% {8 O z7 ?2 r/ j8 j E* H/ |表1 问题一的数值结果
5 a* D; o5 P3 u. h2 l参数名称/ q7 B7 y5 p4 w$ R( L
部件名称( }2 B+ Z5 [' p/ A; P9 T6 ?1 K
出口总温 出口总压 出口流量 功 功率
b. D8 f. x! i9 t; X风扇 379.2843 1.3057 19.0477 135391.2 2578893.0
& K+ H9 o6 J" a+ \% [! _CDFS 420.3160 1.7973 17.1330 41573.3 712275.63 \; g( W6 L6 K
针对问题二:8 n6 |1 R0 b% m
此问题是发动机处于双涵(涡扇)工作模式的研究,首先,对此模式下发动
. A% I- v& {$ H8 H. H5 {: M6 I9 i+ d& D机各个部件进行了建模,之后根据气动热力学过程对各部件模型进行整合建成整/ i* E/ o3 M& M d0 @& r+ i
机模型,确定非线性方程组。$ g0 @' P+ ?& p4 [$ Y
2
* \- x& H" i7 @! P7 J! G其次,根据发动机的平衡方程确立非线性方程组的求解约束条件:3 h c1 G1 U2 L6 A# Q
1 | ( - * ERR NCL NTL mL ) / NTL *mL| (1)
3 o3 y ? [7 E# I: O2 | ( - * ) / * | CH CDFS TH mH TH mH ERR N N N N (2)
0 u* w9 S* a: c1 I l3 Y' '" }5 f' u7 y/ H' t ]: ^9 P
41 41 41 ERR3 | (Wg -Wg ) /Wg | (3)6 f. u# O* x( n6 y1 `2 G$ ~
' '
9 R6 n0 K" X+ f' H& A45 45 45 ERR4 | (Wg -Wg ) /Wg | (4)
1 P* }. e& W& `* D61 62 62 ERR5 | ( p - p ) / p | (5)! Y9 }8 i# g& g3 C m9 V1 t0 l) }
' '
0 n& g4 L" A S: f5 c8 8 ERR6 | (A8 - A ) / A | (6), D- \& N# x- G
2 21 13 21 13 7 | ( - - ) / ( ) | a a a a a ERR W W W W W (7)" x. w% T2 Z" @* a. o
最后,采用遗传算法解非线性方程。但由于遗传算法为仿生算法,每次运行9 F& U, I7 E1 E& ?
程序的结果可能不尽相同,但是每次运算的误差是保证在10%以内的,这保证了+ l* N; r8 w8 n. ~
模型及结果的有效性。求解结果如下表2 所示:
( f- E+ W0 Q, \; W表2 非线性方程组的求解结果; z" C/ ~4 m, ~' a R, s
变量名称 H n CL z CDFS z CH z 4 T TH z TL z- g8 @% e; ?+ _8 m4 E m2 g: a+ a
求解结果 0.9157 0.3332 0.6367 0.5835 1551.490 0.6620 0.89398 j5 @( m4 p! |) r0 F0 ?
针对问题三:' F* ^7 _7 ?8 B9 ^
此问题中发动机处于单涵模式工作状态,本问题其实是在发动机最优性能的1 r/ M( n9 b4 H+ R% {: G
条件下,通过调节规律t 4 T const,求出三个导叶角的最优组合值。首先,在性
?0 A! i6 G" z% q% X) U) T能最优的条件下,从发动机设计全局出发,建立单涵模式下发动机性能最优模型。
( |# q! U, ~- T* _9 R$ J# K通过对实际情况的考虑,把发动机的单位推力最大化作为评价标准,模型求解结+ [8 R9 P5 a7 ?( F6 V# {2 @/ X
果如下表3 所示:
2 d8 c- }) z& o1 ^* [* B表3 发动机性能最优时各变量取值
5 X& s) ]% `% K, T, ~& K. M变量
( N- }0 O1 _6 N. S9 _名称 H n CL z CDFS z CH z 4 T TH z TL z C TL NL
; z# J& Q3 B* q l求解/ ~- r* [# p- F. n0 q( m
结果 0.997 0.958 0.045 0.032 1700 0.025 0.020 8.431 3.296 0.888
/ V) w" p, @) V! }8 Y R对于第三大问中的问题(2)由于能力有限,没能达到较好的收敛效果,变
: \! q% u) o' J: Q, a+ h量较多导致结果跳变较大,最后决定每次合理的确定三个所求量中的两个来对
, u0 E: s$ ^) F: M+ J3 G6 r第三个进行求值,以获得不同马赫数下的变量,进而得出变化规律。如下图所. v' ~9 j. `0 c/ a) F1 g
示:
$ l% c$ s- W" m; [% r/ Z(a)推力随a M 变化规律 (b)CDFS 导叶角度变化
3 }7 e) P2 t5 n4 \+ N(c)低压涡轮导叶角变化 (d)喉管面积随a M 变化规律
. s( ]$ m& w4 l- x% O( E7 |图1 单涵模式发动机性能最优几何变量随马赫数的变化规律
% U C2 ]5 w0 R* \0 I关键词:变循环航空发动机,部件级建模,遗传算法,双涵模式,单涵模; O. D9 s0 H! @9 H6 P8 E+ [
/ b6 ?3 S* d& `1 f- f* d
|
zan
|