变循环发动机部件法建模及优化研究生数学建模竞赛参赛作品.doc
《变循环发动机部件法建模及优化研究生数学建模竞赛参赛作品.doc》由会员分享,可在线阅读,更多相关《变循环发动机部件法建模及优化研究生数学建模竞赛参赛作品.doc(28页珍藏版)》请在三一文库上搜索。
1、参赛密码 (由组委会填写)全第十届华为杯全国研究生数学建模竞赛学 校参赛队号队员姓名参赛密码 (由组委会填写) 第十届华为杯全国研究生数学建模竞赛题 目 变循环发动机部件法建模及优化 摘 要:本文利用附录1、2、3、4给出的特性数据以及计算公式得出了风扇特性表中流量随风扇压比函数值的变化图形规律,求出了给定条件下风扇和CDFS出口的总温、总压和流量,并进一步建立了发动机非线性平衡方程组求解的遗传算法模型,并得到了给定条件下非线性方程组求解结果,进而了解了给定条件下变循环发动机双涵道模式的工作状况;同时还给出了计算发动机性能最优对应的发动机CDFS导叶角度、低压涡轮导叶角度和喷管喉道面积的求法。
2、针对问题一。第一小问,借用附录3中的压气机压比函数值定义式,求得附录4中风扇特性数据表中各换算转速下增压比对应的压比函数值(见附表一),然后由所得压比函数值及表中流量数据画出了流量随压比函数值变化的图形(见图5-1、5-2);第二小问,首先借助题中给定的物理转速及压比函数值利用附录4中风扇和CDFS特性数据表插值并利用附录一中公式(2.7)进行修正得到二者对应的增压比、效率和换算流量,由于风扇进口总温、总压=进气道出口总温、总压,CDFS进口总温、总压=风扇出口总温、总压。借助修正的增压比、效率、换算流量及题中所给初始条件并利用附录中的相关公式可依次得到了风扇和CDFS的出口总温、总压及流量分
3、别为378.333、1.288、19.048;431.803、1.774、16.940。针对问题二。首先利用附录1、2、3、4中的信息得到了题中7个非线性方程的基本参数表达形式。将题中的初始条件代入得到了仅含、七个基本参数的具体形式非线性方程组。本文采用了遗传算法对所得非线性方程组进行了求解,依次进行随机生成初始化群体(即解的初始值)、选择、交叉、变异计算、停止判断操作。这些流程通过C语言编码实现。通过计算得到所求非线性方程组解为指明了问题二给定条件下的变循环发动机双涵道模式运行工作状况。针对问题三。根据本问给定初始条件,可调未知参数共有11个,按照飞机发动机常用的调节规律,选择其中的4个参数
4、赋予初值。实际发动机非线性数学模型中,通常采用数值计算方法,首先给出一组7个独立变量的初值,代入发动机模型,计算各平衡方程的误差,通过一定的非线性方程组迭代解法,不断地修正独立变量的试给值,直至各平衡方程的误差满足设定的迭代精度要求,从而得出各个参数的值,然后由所得参数值,结合附录1中计算发动机性能参数中的推力、单位推力、耗油率公式对变循环发动机的性能进行取优评价,此时对应的发动机CDFS导叶角度、低压涡轮导叶角度和喷管喉道面积即为所求。关键词:变循环发动机;非线性方程组;遗传算法;单/双涵道模式一、问题重述1.1 问题背景由飞机/发动机设计原理可知,对于持续高马赫数飞行任务,需要高单位推力的
5、涡喷循环,反之,如果任务强调低马赫数和长航程,就需要低耗油率的涡扇循环。双涵道变循环发动机可以同时具备高速时的大推力与低速时的低油耗。变循环发动机的内在性能优势,受到了各航空强国的重视,因此有必要对变循环发动机部件法建模问题进行深入研究,进而推动航空事业的发展。1.2 现有信息和数据变循环发动机有两种工作模式,分别为涡喷模式和涡扇模式。发动机在亚音速巡航的低功率工作状态,风扇后的模式转换活门因为副外涵与风扇后的压差打开,使更多空气进入副外涵,同时前混合器面积开大,打开后混合器,增大涵道比,降低油耗,此时为发动机的涡扇模式。发动机在超音速巡航、加速、爬升状态时,前混合器面积关小,副外涵压力增大,
6、选择活门关闭,迫使绝大部分气体进入核心机,产生高的推力,此时为发动机的涡喷模式。燃气涡轮发动机是由进气道、压气机、主燃烧室、涡轮、喷管等部件组成的。如果计算机能够对这些部件的性能进行准确的模拟,那么也就能准确地模拟整个发动机的性能。这种建立在准确模拟发动机各部件性能基础上的发动机性能计算方法,称为部件法。该方法是建立在发动机各部件特性已知的基础上的,因此是计算精度较高的一种方法。发动机各部件匹配工作时,受低压轴功率平衡方程、高压轴功率平衡方程、高压涡轮进口截面流量平衡方程、低压涡轮进口截面流量平衡方程、后混合器静压平衡方程、尾喷管面积平衡方程和风扇出口流量平衡方程7个平衡方程的制约,同时该问题
7、给出了发动机部件计算公式、工质热物理性质参数、气动函数和压气机及涡轮特性等数据。1.3 需要解决的问题1.3.1 流量随压比函数值变化图形的求解以及风扇和CDFS的出口总温、总压和流量的求解。(1)、请画出附录4中风扇特性数据表中流量随压比函数值变化的图形。(2)、设在发动机飞行高度,飞行马赫数的亚音速巡航点,导叶角度均设置为0,风扇和CDFS的物理转速都为0.95,风扇和CDFS的压比函数值都为,求风扇和CDFS的出口总温、总压和流量。1.3.2 请运用或设计适当的算法求解由发动机7个平衡方程组成的非线性方程组。要求陈述算法的关键步骤及其解释,尽可能讨论算法的有效性。设在发动机飞行高度,飞行
8、马赫数的亚音速巡航点,采用双涵道模式,导叶角度均设置为0,选择活门完全打开,副外涵道面积设为,后混合器出口总面积设置为2.8518e+004,尾喷管喉道面积,。请运用或设计适当的算法求解由发动机7个平衡方程组成的非线性方程组。要求陈述算法的关键步骤及其解释,尽可能讨论算法的有效性。1.3.3 发动机性能最优的求解以及尾喷管喉道面积随飞行马赫数变化规律的求解。(1)、设在发动机飞行高度,飞行马赫数的超音速巡航点,发动机采用单涵道模式,将选择活门面积设置为0,风扇导叶角度、高压压气机导叶角度、高压涡轮导叶角度均设置为,后混合器面积设置为2.8518e+004。请问发动机CDFS导叶角度、低压涡轮导
9、叶角度和喷管喉道面积3个量为多少时,发动机的性能最优?(2)、试研究发动机飞行高度,飞行马赫数从变化到,发动机特性最优时,CDFS导叶角度、低压涡轮导叶角度,尾喷管喉道面积随飞行马赫数的变化规律。此时发动机采用单涵道模式,将选择活门面积设置为0,风扇导叶角度、高压压气机导叶角度、高压涡轮导叶角度均设置为,后混合器出口总面积设置为2.8518e+004,后混合器内、外涵道面积可调(即不受附录1后混合器给定的内、外涵道面积值的约束)。二、问题分析2.1 问题一针对本问题。第一小问可借用附录3中的压比函数值定义式,求得附录4中风扇特性数据表中增压比对应的压比函数值,然后由所得压比函数值及流量数据画出
10、流量随压比函数值变化的图形。第二小问采用附录1中发动机进气道、压气机的计算公式,依据假定1按进气道-风扇-CDFS的顺序求解1。其中风扇和CDFS均为压气机部件,在计算中特征流量、增压比、效率均为核心转数、压比函数值、导叶角角度的函数,需采用附录4中风扇特性数据表和CDFS特性数据表插值求解。最后即可得到总温、总压和流量。2.2 问题二针对本问题。分两步求解,首先依据附录1中的发动机部件计算公式再代入附录2、3中的辅助公式列出题目中给出的7个发动机平衡方程,建立仅含、七个变量的待求的非线性方程组。然后针对得到的非线性方程组的特点,选用标准遗传算法,按照算法的步骤编程求解。2.3 问题三当给定发
11、动机部件特性时,同时即描述了发动机的工作过程,但由已建立的方程组唯一确定满足发动机共同工作方程的解,还需给出发动机的调节规律7。针对本问题。本文采用2组具有合理初值的参数分别计算得到发动机各参数的最优解,从而得到优化发动机CDFS导叶角度、低压涡轮导叶角度和喷管喉道面积3个量的方法,使发动机性能最优。在马赫数从变化到过程中,若算法计算较慢,则采用几个离散的马赫数求得对于的最优发动机CDFS导叶角度、低压涡轮导叶角度和喷管喉道面积,然后寻找各自的变化规律;若能得到快速收敛的算法,则取较小马赫数增量步长连续得到一系列最优发动机CDFS导叶角度、低压涡轮导叶角度和喷管喉道面积,得到随马赫数变化的规律
12、,分别画出三个参数随马赫数变化的图形。三、基本假设通过阅读有关参考资料和对数据进行简单的分析,我们对该变循环发动机模型做如下假设:3.1 高压压气机后不经主燃烧室的分流气流为冷却气流,在本题中忽略不计。3.2 如果计算机能够对变循环发动机部件的性能进行准确的模拟,那么也就能够准确地模拟整个发动机的性能。3.3 在前混合器计算中,假定CDFS涵出口流量,即为通过CDFS特性数据线性插值得到的CDFS计算流量减去高压压气机特性数据线性插值得到的高压压气机计算流量。3.4 如无特殊要求,发动机中直接相连的下一级主要部件进口总温、总压假定为上一级主要部件出口处的总温、总压(给出总压恢复系数的另外考虑)
13、,直接相连的发动机主要部件顺序如下:风扇-CDFS-高压压气机-主燃烧室-高压涡轮-低压涡轮。例如:风扇计算得到风扇出口的总温、总压,这一总温、总压即为CDFS的进口总温、总压,以此类推。四、符号说明五、模型建立和求解5.1 求解问题一5.1.1 风扇特性数据表中流量随压比函数值变化的图形由附录3中的压气机压比函数值的定义可知,压比函数值zz可以如下算得:设压气机某换算转速所对应的增压比数据(见附录4)的最大值为,最小值为,则定义该换算转速对应的压气机增压比的压比函数值为: (5.1.1-1)由附录4中风扇特性数据,利用公式(5.1.1-1)可分别算出换算转速为0.4、0.5、0.6、0.7、
14、0.81、0.9、0.95、1和1.075对应的压比函数值。然后利用得到的压比函数值数据与流量数据可以得风扇流量随压比函数值变化的图形。结果如下所示:图5-1 风扇流量随压比函数值变化图形(一)图5-2 风扇流量随压比函数值变化图形(二)5.1.2 给定条件下风扇及CDFS出口的总温、总压和流量(1)、计算思路由附录1中公式(2.1)(2.4)算出进气道的出口总温总压,由“风扇进口总温总压=进气道出口总温总压”结合风扇的特性数据算出风扇出口总温总压,再利用“CDFS进口总温总压=风扇出口总温总压”结合CDFS特性数据及相关公式算出CDFS的出口总温总压;风扇和CDFS出口的流量则可以直接利用附
15、录4中的特性数据进行插值得到。(2)、给定条件题中给定的初始条件有:飞行高度H=11km;飞行马赫数为Ma=0.8的亚音速巡航点;各导叶角度均设置为0;风扇和CDFS物理转速均为0.95;风扇和CDFS的压比函数值均为0.5。(3)、算法流程进气道出口总温总压的计算利用附录1中的公式(2.1)(2.4)结合上面给定条件算出进气口的总温、总压分别为:244.38、0.345。风扇出口总温总压及流量的计算A、总温总压计算首先,计算修正后的增压比、效率和换算流量。利用给定的物理转速n=0.95及附录1中的公式(2.5)计算出换算转速;利用及给定的压比函数值zz=0.5从附录4中风扇特性数据表线性插值
16、得到增压比、效率和换算流量;再结合附录1公式(2.7)算出修正后的增压比、效率和换算流量。之后,计算出口总温和总压。由假设有:风扇进口总压=进气道出口总压,则有:由假设有:风扇进口总温=进气道出口总温,再结合附录1中“2.2.2 计算过程3) 计算压气机出口参数”的公式算法步骤算出风扇出口总温。B、出口流量计算利用附录4中的风扇特性数据线性插值并修正求出来的流量结合附录1中的公式(2.8)求出风扇出口流量。CDFS出口总温总压及流量的计算按与风扇同样的方法步骤算出修正后的增压比、效率和换算流量;由于CDFS进口总压=风扇出口总压,CDFS进口总温=风扇出口总温。再利用与风扇总温总压同样的求解方
17、法可以得出CDFS的出口总温总压;同样采用附录4中CDFS特性数据线性插值并修正获得的流量结合附录1中的公式(2.8)求出CDFS出口流量。(4)、计算结果经过计算得到风扇和CDFS的出口总温、总压及流量结果如下表:表5-1 风扇和CDFS的出口总温、总压及流量出口总温出口总压出口流量风扇378.3331.28819.048CDFS431.8031.77416.9405.2 求解问题二如无特殊要求,发动机中直接相连的下一级主要部件进口总温、总压假定为上一级主要部件出口处的总温、总压(给出总压恢复系数的另外考虑),直接相连的发动机主要部件顺序如下:风扇-CDFS-高压压气机-主燃烧室-高压涡轮-
18、低压涡轮。例如:风扇计算得到风扇出口的总温、总压,这一总温、总压即为CDFS的进口总温、总压,以此类推。5.2.1 展开发动机平衡方程组(1)、展开发动机平衡方程组一 (5.2.1-1) (2)、展开发动机平衡方程二 (5.2.1-2) (3)、展开发动机平衡方程三 (5.2.1-3) (4)、展开发动机平衡方程四 (5.2.1-4)(5)、展开发动机平衡方程五 (5.2.1-5) (6)、展开发动机平衡方程六 (5.2.1-6) (7)、展开发动机平衡方程七 (5.2.1-7)对于第二问所给条件,本方程组7个一共有7个待定变量,他们分别是、。5.2.2 遗传算法模型的建立与求解由于发动机的7
19、个平衡方程组成的非线性方程组,利用传统的非线性方程组求解,无论是从算法的选择还是算法本身的构造都与索要解决的问题的特性有很大的关系,很多情况下算法中算子的构造及有效性会成为我们解决问题的巨大障碍。而遗传算法作为一种灵活的自适应算法无需过多的考虑问题的具体形式,具备全局收敛性,突破了传统算法的诸多限制与不足,因而本文采用遗传算法进行发动机非线性方程组的求解。(1)、非线性方程组求解转化成遗传算法的优化问题设有限空间内的非线性方程组:其中:常量、为相应自变量的上下限。记,上述方程组可表示为:构造非线性函数为n维欧式空间的一个有界区域。在整个自变量的有界区域,当上述非线性方程组有解时,函数满足;即,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 循环 发动机 部件 建模 优化 研究生 数学 竞赛 参赛 作品
链接地址:https://www.31doc.com/p-2081077.html