TA的每日心情 | 奋斗 2024-7-1 22:21 |
|---|
签到天数: 2014 天 [LV.Master]伴坛终老
- 自我介绍
- 数学中国站长
 群组: 数学建模培训课堂1 群组: 数学中国美赛辅助报名 群组: Matlab讨论组 群组: 2013认证赛A题讨论群组 群组: 2013认证赛C题讨论群组 |
摘 要:
8 d7 j/ t6 L3 E7 l+ i/ k. U本文对变循环发动机部件法建模,非线性方程与方程组及多变量性能寻优
5 t7 X+ ^- }% ~问题进行了研究。
2 _' k' l' {4 E. \! {7 q在问题一中,本文首先利用附录3 提供的压比函数值定义式计算即可得到) z: \6 q& j2 l/ A
各换算转速、各增压比值对应的压比函数值,然后画出流量随压比函数值变化
3 k6 ?4 q/ w, t8 a- k3 p的图形(图2)。接着利用附录1 及附录3 的公式在Matlab 中编程求解,其中. p( v9 h! v' P
增压比和换算流量涉及到了插值计算。求解结果如下:风扇的出口总温
5 d2 G/ a+ f' {! r- L378.34852K 、总压1.28834bar 和流量19.04771,CDFS 的出口总温419.37667K 、
5 D4 [; f! F D- _" [% M2 `出口总压1.77447bar 和流量16.93963。 ^: Q) G$ i9 P
针对问题二首先通过部件法建模法,将进气道、风扇、CDFS 等12 个部件
" |7 D5 I1 E+ w8 C) c8 t) y的特性计算组成涉及7 个变量的发动机全流程参数计算,而7 个平衡方程组成2 W# C }4 {* a/ u$ i
了非线性方程组。该非线性方程组无明显的数学表达形式,因此不可能得到解
% U. c6 z2 M; r析解,需要采用数值计算方法来求解。问题的实质是求解一组隐式非线性方程
6 Q6 q2 G# F& ^) ^/ \组,满足残差最小的优化目标。本文首先分别应用牛顿-拉弗森算法与遗传算法) Q' u+ l- R4 M" o6 p# ^9 g9 I
分别对变循环发动机模型的非线性方程组求解。考虑到牛顿-拉弗森算法与遗传' U6 P. f- M3 H! z8 H C5 V
算法的各自优点,本文设计了遗传算法与牛顿-拉弗森算法相结合的混合算法。
% e/ w" s. F5 p8 S, A" ^在所设计的算法流程中,首先由牛顿-拉弗森算法求解,对于牛顿-拉弗森算法
, u* ?1 C; p B) \! |) j; K求解不收敛的工作点,则调用遗传算法重新确定一点,作为牛顿-拉弗森算法新
3 _- m# ~5 ^% C2 f3 S: S的迭代初值,如此反复。既保持了牛顿-拉弗森算法的高计算效率,又吸收了遗' ~& t* C8 _) _6 U* _" u1 }
传算法全局收敛的优点,在保持与现有模型计算效率一致的情况下,彻底突破
* n; m- ]9 H9 z7 Q1 ]* B了发动机非线性数学模型中平衡方程与非线性方程常规解法的局限,实现了模
/ O1 r0 c+ T# W20 I8 u! d V% o" d" S" m
型在整个包线范围内的收敛。经过综合求解,最后出 7个变量的值为高速 转速 0.879 、主燃烧出口温度 、主燃烧出口温度 1520 、 风扇、 风扇CDFSCDFS CDFS、压气机高涡 轮低、压气机高涡 轮低、压气机高涡 轮低、压气机高涡 轮低轮的压比函数分别为 0.554 、0.634 、0.793 、0.213 、0.0254 。
" \/ w8 ~6 I) W4 M2 C3 b问题三 的实质是单涵道模式下,变循环发动机 CDFSCDFS CDFS导叶角度、低压涡轮 导叶角度和喷管喉道面积 三个 可调 变量在各自的范围内,需要对可调 变量进行性能寻优, 使得设 定的性能目标达到最优 。本文 采用遗传算法对变循 环发动机的变几何量进行编码, 执遗传操作对某一性能目标寻优分别 对变循环发动机的在 固定和变化 马赫数下进行最大推力、小耗油率 两种优化 目标的 性能寻 优控制 ,结果见表 3、表 4。
( K) ^3 v7 u9 ]本文综合运用数理统计、 牛顿 -拉弗森 法 、遗传算法 、遗传算等数学 方法, 利用 Matlab Matlab软件, 结合多种模型对所提出问题进行了研究,具有很好的实用性与推广。 最后,总结了模型的 最后,总结了模型的 优点与不足 ,并 提出了基于 牛顿 -拉弗森 法和 BP 神经网络 的混合模型,为后续研究此类问题学者提供了一个新思路 。5 W# a3 S! r, n; W, E# }& r8 }' ?
关键词 关键词 :变循环 发动机、 发动机、 部件 法建模、 牛顿 -拉弗森 法、 遗传算BP 神经网 络$ M# M! f$ o0 j, ?+ m3 g% M
- y3 k; y$ [! s9 y
|
zan
|