TA的每日心情 | 奋斗 2024-7-1 22:21 |
|---|
签到天数: 2014 天 [LV.Master]伴坛终老
- 自我介绍
- 数学中国站长
 群组: 数学建模培训课堂1 群组: 数学中国美赛辅助报名 群组: Matlab讨论组 群组: 2013认证赛A题讨论群组 群组: 2013认证赛C题讨论群组 |
摘 要:( o5 @( z$ e1 t8 ] \. p
本文对变循环发动机部件法建模,非线性方程与方程组及多变量性能寻优7 F; s: p& v [ M E$ G
问题进行了研究。! d( } A9 t4 \: I. l
在问题一中,本文首先利用附录3 提供的压比函数值定义式计算即可得到
* X$ S0 |# t4 M9 n0 _! r4 v: x各换算转速、各增压比值对应的压比函数值,然后画出流量随压比函数值变化
. E3 D# G5 t* o- t的图形(图2)。接着利用附录1 及附录3 的公式在Matlab 中编程求解,其中
- T' ?, n9 T2 N% O. z8 V增压比和换算流量涉及到了插值计算。求解结果如下:风扇的出口总温3 K/ D) s0 a8 v& |+ r
378.34852K 、总压1.28834bar 和流量19.04771,CDFS 的出口总温419.37667K 、2 l* r8 S/ ~' V. i9 c
出口总压1.77447bar 和流量16.93963。
v/ I% \' g D针对问题二首先通过部件法建模法,将进气道、风扇、CDFS 等12 个部件0 K+ L# E- z5 F( q- B! a( y, k
的特性计算组成涉及7 个变量的发动机全流程参数计算,而7 个平衡方程组成
* w) Q7 ], p8 c1 `6 Q8 @) D了非线性方程组。该非线性方程组无明显的数学表达形式,因此不可能得到解) ? |: I" [, W% _
析解,需要采用数值计算方法来求解。问题的实质是求解一组隐式非线性方程
; F% {6 S4 @8 z. h& W组,满足残差最小的优化目标。本文首先分别应用牛顿-拉弗森算法与遗传算法4 {! `0 x& k. k. C" l6 A
分别对变循环发动机模型的非线性方程组求解。考虑到牛顿-拉弗森算法与遗传
! w0 Y2 Y. U" f/ X算法的各自优点,本文设计了遗传算法与牛顿-拉弗森算法相结合的混合算法。- A, J8 E) l! E9 F
在所设计的算法流程中,首先由牛顿-拉弗森算法求解,对于牛顿-拉弗森算法
. A- L k# r6 ~" r; |2 h$ p3 R求解不收敛的工作点,则调用遗传算法重新确定一点,作为牛顿-拉弗森算法新- X1 k( ^, K0 P7 z& g* K) I
的迭代初值,如此反复。既保持了牛顿-拉弗森算法的高计算效率,又吸收了遗; }3 E9 q: G: B% N" l# J1 A
传算法全局收敛的优点,在保持与现有模型计算效率一致的情况下,彻底突破 T+ u! P' u( g' T
了发动机非线性数学模型中平衡方程与非线性方程常规解法的局限,实现了模2 B5 x% c) B2 f8 E6 Z
2
5 ]7 s8 N( y0 a# N. a B$ I5 F/ N型在整个包线范围内的收敛。经过综合求解,最后出 7个变量的值为高速 转速 0.879 、主燃烧出口温度 、主燃烧出口温度 1520 、 风扇、 风扇CDFSCDFS CDFS、压气机高涡 轮低、压气机高涡 轮低、压气机高涡 轮低、压气机高涡 轮低轮的压比函数分别为 0.554 、0.634 、0.793 、0.213 、0.0254 。- A: ]! [6 m1 b3 ~( E
问题三 的实质是单涵道模式下,变循环发动机 CDFSCDFS CDFS导叶角度、低压涡轮 导叶角度和喷管喉道面积 三个 可调 变量在各自的范围内,需要对可调 变量进行性能寻优, 使得设 定的性能目标达到最优 。本文 采用遗传算法对变循 环发动机的变几何量进行编码, 执遗传操作对某一性能目标寻优分别 对变循环发动机的在 固定和变化 马赫数下进行最大推力、小耗油率 两种优化 目标的 性能寻 优控制 ,结果见表 3、表 4。
* q4 |, K) B/ e本文综合运用数理统计、 牛顿 -拉弗森 法 、遗传算法 、遗传算等数学 方法, 利用 Matlab Matlab软件, 结合多种模型对所提出问题进行了研究,具有很好的实用性与推广。 最后,总结了模型的 最后,总结了模型的 优点与不足 ,并 提出了基于 牛顿 -拉弗森 法和 BP 神经网络 的混合模型,为后续研究此类问题学者提供了一个新思路 。8 X8 e5 L( r$ l, G4 }% |
关键词 关键词 :变循环 发动机、 发动机、 部件 法建模、 牛顿 -拉弗森 法、 遗传算BP 神经网 络) l1 B* O" J9 J( b$ l) O t+ Y
, C! Z0 `: a$ ?' d! K
|
zan
|