TA的每日心情 | 奋斗 2024-7-1 22:21 |
|---|
签到天数: 2014 天 [LV.Master]伴坛终老
- 自我介绍
- 数学中国站长
 群组: 数学建模培训课堂1 群组: 数学中国美赛辅助报名 群组: Matlab讨论组 群组: 2013认证赛A题讨论群组 群组: 2013认证赛C题讨论群组 |
摘 要:
$ S+ L0 k! |, @本文对变循环发动机部件法建模,非线性方程与方程组及多变量性能寻优
2 z2 d0 {% P1 N$ ]# T% y9 ^; p# R问题进行了研究。
) j5 f7 ]: N8 ]在问题一中,本文首先利用附录3 提供的压比函数值定义式计算即可得到
; D5 [' D0 W3 f A$ J6 |4 ?4 v各换算转速、各增压比值对应的压比函数值,然后画出流量随压比函数值变化7 H/ P) `" B4 x5 z4 h" F$ |* Z1 o
的图形(图2)。接着利用附录1 及附录3 的公式在Matlab 中编程求解,其中+ D6 g# p$ |. {% B# A! C& U
增压比和换算流量涉及到了插值计算。求解结果如下:风扇的出口总温
* V5 a" A- u! f, Y4 [. k' c0 r3 B7 c$ O378.34852K 、总压1.28834bar 和流量19.04771,CDFS 的出口总温419.37667K 、
/ ^9 O. f5 Z- U9 T$ y出口总压1.77447bar 和流量16.93963。2 Z0 r$ q* K$ @% m
针对问题二首先通过部件法建模法,将进气道、风扇、CDFS 等12 个部件
! P- }# Y, Z' E& E0 U1 k的特性计算组成涉及7 个变量的发动机全流程参数计算,而7 个平衡方程组成
% e$ S" P! ?* H% r9 Q3 M# e4 y了非线性方程组。该非线性方程组无明显的数学表达形式,因此不可能得到解6 G' u9 Z- `0 n* |3 b0 M. F
析解,需要采用数值计算方法来求解。问题的实质是求解一组隐式非线性方程
/ z' @4 Z7 c2 G组,满足残差最小的优化目标。本文首先分别应用牛顿-拉弗森算法与遗传算法
0 e5 Q+ I+ k4 z! Y分别对变循环发动机模型的非线性方程组求解。考虑到牛顿-拉弗森算法与遗传3 S/ J9 N2 g4 }1 U
算法的各自优点,本文设计了遗传算法与牛顿-拉弗森算法相结合的混合算法。7 c# w# v* l: o& m
在所设计的算法流程中,首先由牛顿-拉弗森算法求解,对于牛顿-拉弗森算法
" t2 U$ F v& M/ W求解不收敛的工作点,则调用遗传算法重新确定一点,作为牛顿-拉弗森算法新; |% i6 g$ z" l+ o- L
的迭代初值,如此反复。既保持了牛顿-拉弗森算法的高计算效率,又吸收了遗
7 p3 C, x, |" i8 |: B传算法全局收敛的优点,在保持与现有模型计算效率一致的情况下,彻底突破
o& B- w* e5 [- }8 B/ R了发动机非线性数学模型中平衡方程与非线性方程常规解法的局限,实现了模
0 C% k$ j; n+ m; y5 x7 p2. M" g3 n) R& O! ?
型在整个包线范围内的收敛。经过综合求解,最后出 7个变量的值为高速 转速 0.879 、主燃烧出口温度 、主燃烧出口温度 1520 、 风扇、 风扇CDFSCDFS CDFS、压气机高涡 轮低、压气机高涡 轮低、压气机高涡 轮低、压气机高涡 轮低轮的压比函数分别为 0.554 、0.634 、0.793 、0.213 、0.0254 。3 Z2 q, z& u. L F' J: B
问题三 的实质是单涵道模式下,变循环发动机 CDFSCDFS CDFS导叶角度、低压涡轮 导叶角度和喷管喉道面积 三个 可调 变量在各自的范围内,需要对可调 变量进行性能寻优, 使得设 定的性能目标达到最优 。本文 采用遗传算法对变循 环发动机的变几何量进行编码, 执遗传操作对某一性能目标寻优分别 对变循环发动机的在 固定和变化 马赫数下进行最大推力、小耗油率 两种优化 目标的 性能寻 优控制 ,结果见表 3、表 4。
6 r- _* @8 s$ S1 x+ J9 P本文综合运用数理统计、 牛顿 -拉弗森 法 、遗传算法 、遗传算等数学 方法, 利用 Matlab Matlab软件, 结合多种模型对所提出问题进行了研究,具有很好的实用性与推广。 最后,总结了模型的 最后,总结了模型的 优点与不足 ,并 提出了基于 牛顿 -拉弗森 法和 BP 神经网络 的混合模型,为后续研究此类问题学者提供了一个新思路 。
: h& d' b. [. |关键词 关键词 :变循环 发动机、 发动机、 部件 法建模、 牛顿 -拉弗森 法、 遗传算BP 神经网 络: X8 R9 B6 g' ^) ?( [& I" a+ ~
" B5 S/ u/ F. ]5 m' y |
zan
|