CN112711923A - 多级压气机气动匹配设计方法 - Google Patents

多级压气机气动匹配设计方法 Download PDF

Info

Publication number
CN112711923A
CN112711923A CN202110323072.XA CN202110323072A CN112711923A CN 112711923 A CN112711923 A CN 112711923A CN 202110323072 A CN202110323072 A CN 202110323072A CN 112711923 A CN112711923 A CN 112711923A
Authority
CN
China
Prior art keywords
design
stage
compressor
dimensional
flow
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110323072.XA
Other languages
English (en)
Other versions
CN112711923B (zh
Inventor
姜逸轩
曹传军
吴帆
尹泽勇
李继保
吴志青
李游
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
AECC Commercial Aircraft Engine Co Ltd
Original Assignee
AECC Commercial Aircraft Engine Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by AECC Commercial Aircraft Engine Co Ltd filed Critical AECC Commercial Aircraft Engine Co Ltd
Priority to CN202110323072.XA priority Critical patent/CN112711923B/zh
Publication of CN112711923A publication Critical patent/CN112711923A/zh
Application granted granted Critical
Publication of CN112711923B publication Critical patent/CN112711923B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12—Simultaneous equations, e.g. systems of linear equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Geometry (AREA)
  • Computational Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Algebra (AREA)
  • Data Mining & Analysis (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Automation & Control Theory (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

一种多级压气机气动匹配设计方法，其基于等出口换算流量的原理并以保持“等熵效率不变”为条件，计算获得经过单级设计工作点的等出口换算流量线；在获取单级的等出口换算流量线之后，通过设计系统偏差修正，沿等出口换算流量线，获取修正后的单级设计目标，进而开展单级迭代设计；具体开展某一特定级三维设计工作时，该级进口条件使用S2设计环节获得的二维总温、总压、气流角剖面，使三维设计与一维、二维设计结果保持一致，即压气机各设计环节保持一致。本发明在满足压气机整机气动布局设计目标前提下能有效地将单级设计从多级设计工作中剥离出来，并在不影响压气机整体气动布局的前提下完成单级压气机的气动设计，从而优化了压气机设计流程。

Description

多级压气机气动匹配设计方法
技术领域
本发明涉及航空发动机技术领域，尤其涉及航空发动机高压压气机气动领域，更具体地涉及一种基于等出口换算流量原理的多级压气机气动匹配设计方法。
背景技术
燃气轮机是以连续流动的气体为工质带动叶轮高速旋转，从而将燃料的能量转变为有用功的内燃式动力机械，其实质是一种旋转叶轮式热力发动机，其基本工作原理是：压气机从外界大气环境吸入空气，并经过轴流式压气机逐级压缩使之增压，同时空气温度也相应提高；压缩空气被压送到燃烧室与喷入的燃料混合燃烧生成高温高压的气体；然后再进入到涡轮中膨胀做功，推动涡轮带动压气机和外负荷转子一起高速旋转，实现了气体或液体燃料的化学能部分转化为机械功，并输出电功。
在航空发动机中，压气机是最核心的部件之一，发动机的性能通常取决于压气机的气动性能。压气机是用于压缩气体的部件，一般为多级设计，其各级之间的气动匹配则是影响气动性能的关键因素。因而，在压气机初始设计阶段，需要通过多方面考量获取一个最优的多级压气机气动布局，在这之后再对各级进行单独设计。而多级压气机作为一个整体，其各级之间具有非常强的关联性，各级之间气动性能相互影响。因此，在进行某一级的单独设计时，该级上下游各级的性能乃至整个压气机的气动布局都有可能随之变化，由此可能会对压气机整体气动布局造成不利地影响。因此，如何在不影响压气机整体气动布局的前提下完成单级压气机的气动设计，已然成为了压气机设计过程中极具挑战性的问题。
例如，现有文献CN109165440A（公布日期为2019年1月8日）公开了一种轴流压气机全三维级间气动匹配优化方法，该方法对给定的轴流压气机叶片几何，采用RANS方法，进行全三维粘性内流场数值模拟，根据数值模拟结果，构造静叶各设计截面设计目标载荷，对于转子，计算得到动叶设计截面目标载荷分布，给定出口流量边界条件，保证优化前后流量不变，通过实际计算流量和规定设计流量的差距，在初始背压的基础上，自动完成出口背压的调节，最终满足流量设计要求。
显然，该现有文献CN109165440A采用的方法是根据出口流量来调节出口背压，而非基于等出口换算流量进行单级设计的方法。
再比如，现有文献CN111079236A（公布日期为2020年4月28日）公开了一种用于舰船燃气轮机压气机增容的母型机匹配点获取方法，该方法包括：结合设定的改型压气机的流量和总压比，得到母型压气机的流量和总压比，将得到的母型压气机的流量和总压比绘制在母型机通用特性图上作为母型压气机上的两个匹配点，连接两个匹配点，得到一条直线，满足选定改型机设计参数（压比、流量）的所有可用匹配点都在这条直线及其延伸线上；然后，再改变改型机设计目标（流量、压比），重复计算步骤，可以得到一系列直线，每个直线对应一个设计流量和压比，根据一系列直线簇，可以根据理论上的母型机匹配点，快速分析改型机设计点性能，确定最佳设计参数。
显然，该现有文献CN111079236A的方法是在分布直线两端的两个匹配点对应不同的转速，而非计算得到等出口换算流量线的两个端点的方法。
然而，多级压气机各级之间耦合度极高，各级之间相互影响，增加了压气机的设计难度。如何有效地将单级设计从多级设计中剥离出来，是压气机设计工作中至关重要的问题。显然，上述现有文献都无法有效地解决目前存在的技术难题。
鉴于上述内容可知，目前在本技术领域中还尚不存在这样一种多级压气机气动匹配设计方法，该方法基于等出口换算流量的原理，并以保持“等熵效率不变”为条件，计算获得经过单级设计工作点的等出口换算流量线；在获取单级的等出口换算流量线之后，通过设计系统偏差修正，沿等出口换算流量线，获取修正后的单级设计目标，进而开展单级迭代设计；具体开展某一特定级三维设计工作时，该级进口条件使用S2设计环节获得的二维总温、总压、气流角剖面，使三维设计与一维、二维设计结果保持一致，即压气机各设计环节保持一致。因此，如何能够设计一种实现上述效果的多级压气机气动匹配设计方法俨然成为了亟待解决的技术问题。
发明内容
本发明正是为了解决上述技术问题而作，其目的在于提供一种多级压气机气动匹配设计方法，该方法在满足压气机整机气动布局设计目标前提下，能够有效地将单级设计从多级设计工作中剥离出来，并且在不影响压气机整体气动布局的前提下完成单级压气机的气动设计，从而优化了压气机设计流程。
为了解决上述技术问题，本发明提供了一种多级压气机气动匹配设计方法，其特征在于，所述方法包括以下步骤：
步骤一：开展第i级单级设计，基于一维设计结果，获取单级等出口换算流量线。
步骤二：通过设计系统偏差修正，沿等出口换算流量线，获取修正后的单级设计目标，进而开展第i级单级迭代设计。
较佳地，在本发明的多级压气机气动匹配设计方法中，在步骤一中，通过压气机一维设计，可以获得压气机进口物理流量
Figure 708969DEST_PATH_IMAGE001
，以及第i级单级设计工作点进口换算流量
Figure 777173DEST_PATH_IMAGE002
、总压比
Figure 172383DEST_PATH_IMAGE003
、总温比
Figure 658859DEST_PATH_IMAGE004
、单级效率
Figure 746900DEST_PATH_IMAGE005
、单级平均比热容比
Figure 200884DEST_PATH_IMAGE006
。
较佳地，在本发明的多级压气机气动匹配设计方法中，在步骤一中，求解所述等出口换算流量线上的其中一个端点的坐标，取进口换算流量偏置量
Figure 134205DEST_PATH_IMAGE007
（
Figure 740767DEST_PATH_IMAGE008
为定值），即
Figure 999710DEST_PATH_IMAGE009
………（1）
所述其中一个端点的进口换算流量计算公式为
Figure 675411DEST_PATH_IMAGE010
………（2）
所述其中一个端点的出口换算流量计算公式为
Figure 678002DEST_PATH_IMAGE011
………（3）
式（3）/（2）得到
Figure 139070DEST_PATH_IMAGE012
………（4）
由于已知条件
Figure 303336DEST_PATH_IMAGE013
Figure 544961DEST_PATH_IMAGE014
Figure 272614DEST_PATH_IMAGE015
Figure 650506DEST_PATH_IMAGE016
………（5）
将（5）中各式代入式（4），式（4）可变为
Figure 188935DEST_PATH_IMAGE017
………（4）
第i级单级效率EFF计算公式为
Figure 652277DEST_PATH_IMAGE018
………（5）
由“保持等熵效率不变”的条件可得
Figure 183622DEST_PATH_IMAGE019
………（6）
因此式（4）变为
Figure 681599DEST_PATH_IMAGE020
………（7）
联立式（4）、（7）可得关于
Figure 125350DEST_PATH_IMAGE021
、
Figure 341568DEST_PATH_IMAGE022
两个未知量的二元方程组，求解该方程组可获得
Figure 411024DEST_PATH_IMAGE021
，从而可以获得所述等出口换算流量线上的其中一个端点的坐标位置（
Figure 763508DEST_PATH_IMAGE023
，
Figure 378160DEST_PATH_IMAGE021
）。
较佳地，在本发明的多级压气机气动匹配设计方法中，在步骤一中，采用获得所述等出口换算流量线上的其中一个端点的坐标位置的相同的方法，求得当
Figure 816094DEST_PATH_IMAGE024
时，所述等出口换算流量线上的另一个端点的坐标位置（
Figure 689241DEST_PATH_IMAGE025
，
Figure 896232DEST_PATH_IMAGE026
），从而连接这两个端点即可获得所述等出口换算流量线。
较佳地，在本发明的多级压气机气动匹配设计方法中，在步骤二中，即使所述第i级设计工作点出现偏离，该级的出口换算流量仍然保持不变，因而不会影响该级与下游各级的流量匹配关系，从而保持多级压气机气动设计布局。
较佳地，在本发明的多级压气机气动匹配设计方法中，在步骤二中，开展第i级单级三维设计工作，该级进口条件采用二维S2设计获得的二维总温、总压、气流角剖面设计结果，使三维设计与一维、二维设计结果保持一致，即压气机各设计保持一致。
如上所述，多级压气机的匹配需要考虑多种因素，其中流量匹配是最重要的设计思路。由此，本发明多级压气机气动匹配设计方法的核心技术在于能够在确保单级出口换算流量不变的条件下保持上下游各级之间的气动布局。在完成压气机整体气动布局设计的前提下，单级压气机设计基于等出口换算流量原理，通过设计系统偏差修正，沿等出口换算流量线，获取修正后的单级设计目标，进而开展单级迭代设计。沿等出口换算流量线进行迭代设计，该级的设计工作点发生变化时，其出口换算流量保持不变，即下一级的进口换算流量保持不变，因而下一级流量特性不会发生变化。由于从前往后各级均采用该思路设计，因而压气机从前往后的流量特性保持不变，从而确保了压气机整体气动布局的稳定。
本发明的多级压气机气动匹配设计方法基于等出口换算流量的原理，并以保持“等熵效率不变”为条件，计算获得经过单级设计工作点的等出口换算流量线；在获取单级的等出口换算流量线之后，通过设计系统偏差修正，沿等出口换算流量线，获取修正后的单级设计目标，进而开展单级迭代设计；具体开展某一特定级三维设计工作时，该级进口条件使用S2设计环节获得的二维总温、总压、气流角剖面，使三维设计与一维、二维设计结果保持一致，即压气机各设计环节保持一致。
鉴于上述内容，与现有技术相比，本发明多级压气机气动匹配设计方法具备以下显著的有益效果：
第一、基于等出口换算流量原理，在单级压气机设计工作点发生变化时，保持该级出口换算流量不变，减弱对上下游各级的影响，从而确保了压气机整体气动布局不变；
第二、在满足压气机整体气动布局不变的前提下，为压气机单级设计提供了合理的设计目标，使得耦合度极高的压气机多级设计，简化为单级设计；
第三、实现了压气机单级设计从多级设计剥离出来的目的，简化设计工作；以及
第四、某一特定级三维设计时，进口条件使用S2设计环节获得的二维总温、总压、气流角剖面，使三维设计与一维、二维设计结果保持一致，即压气机各设计环节保持一致。
附图说明
为了更清楚地说明本发明提供的技术方案，下面将对附图作简单地介绍。显而易见地，下面描述的附图仅仅是本发明的一较佳实施例。
图1为本发明多级压气机气动匹配设计方法的一较佳实施例的流程图；
图2为本发明多级压气机气动匹配设计方法的设计系统偏差示意图；
图3为本发明多级压气机气动匹配设计方法的单级压气机进口总压、总温、气流角二维剖面示意图。
具体实施例
为了使得本发明实施的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行更为详细地描述。
对此，首先需要指出的是，在这些实施例的具体描述过程中，为了进行简明扼要的描述，本说明书不可能对实际的实施例的所有特征均作详尽的描述。应当可以理解的是，在任意一种实施例的实际实施过程中，正如在任意一个工程项目或者设计项目的过程中，为了实现开发者的具体目标，为了满足系统相关的或者商业相关的限制，常常会做出各种各样的具体决策，而这也会从一种实施例到另一种实施例之间发生改变。此外，还可以理解的是，虽然这种开发过程中所作出的努力可能是复杂并且冗长的，然而对于与本发明公开的内容相关的本领域的普通技术人员而言，在本公开揭露的技术内容的基础上进行的一些设计、制造或者生产等变更只是常规的技术手段，不应当理解为本公开的内容不充分。
另外，需要说明的是，除非另作定义，权利要求书和说明书中使用的技术术语或者科学术语应当为本发明所属技术领域内具有一般技能的人士所理解的通常意义。“一个”或者“一”等类似词语并不表示数量限制，而是表示存在至少一个。“包括”或者“包含”等类似的词语意指出现在“包括”或者“包含”前面的元件或者物件涵盖出现在“包括”或者“包含”后面列举的元件或者物件及其等同元件，并不排除其它元件或者物件。“连接”或者“相连”等类似的词语并非限定于物理的或者机械的连接，也不限于是直接的还是间接的连接。
总体而言，本发明的多级压气机气动匹配设计方法基于等出口换算流量的原理，并以保持“等熵效率不变”为条件，计算获得经过单级设计工作点的等出口换算流量线；在获取单级的等出口换算流量线之后，通过设计系统偏差修正，沿等出口换算流量线，获取修正后的单级设计目标，进而开展单级迭代设计；具体开展某一特定级三维设计工作时，该级进口条件使用S2设计环节获得的二维总温、总压、气流角剖面，使三维设计与一维、二维设计结果保持一致，即压气机各设计环节保持一致。
下面，首先对以下技术术语的定义先进行解释：
“换算流量”是指：单位时间内从压气机或其中某一级内流过的气体的质量，根据压气机进口或单级进口的总温、总压，换算到标况进口（总压101325Pa、总温288.15K）条件下的得到的流量，单位一般为kg/s；
“匹配”是指：多级压气机，综合考虑多方面因素，对各级之间的的负荷、效率等参数进行的气动布局设计；
“等熵效率”是指：压气机等熵功率与实际压缩工质所需功率之比。
接着，结合图1－3对本发明多级压气机气动匹配设计方法的一具体实施例进行详细阐述，以使本发明的优点和特征能易于被本领域的技术人员理解，从而对本发明范围做出更为清晰的界定。
如图1所示，图中示出了本发明多级压气机气动匹配设计方法的一较佳实施例的流程图。本发明多级压气机气动匹配设计方法的具体方案如下：
首先，如图2所示，图中示出了本发明多级压气机气动匹配设计方法的设计系统偏差示意图。在图2所示的单级“换算流量-压比”图中，获得了一条经过单级设计工作点T的等出口换算流量线，即图2中所示的线段AB。在求解端点A、B坐标的过程中，增加“等熵效率不变”的条件，保证方程组可解。“等熵效率不变”的条件，使得单级设计工作点变化时，其单级等熵效率不变，进一步减小了对上下游各级的气动影响，保证了压气机整体气动布局的一致性。
接着，在如上所述获取单级的等出口换算流量线之后，通过设计系统偏差△修正，沿等出口换算流量线，获取修正后的单级设计目标T’，进而开展单级迭代设计，从而实现了将单级设计从多级设计剥离出来的目的，并且在不影响压气机整体气动布局的前提下完成单级压气机的气动设计，从而优化了压气机设计流程。
最后，在具体开展某一特定级三维设计工作时，该级进口条件使用S2设计环节获得的二维总温、总压、气流角剖面，使三维设计与一维、二维设计结果保持一致，即压气机各设计环节保持一致。
本发明多级压气机气动匹配设计方法的具体实施步骤如下：
步骤一：开展第i级单级设计，基于一维设计结果，获取单级等出口换算流量线。
步骤二：通过设计系统偏差修正，沿等出口换算流量线，获取修正后的单级设计目标，进而开展第i级单级迭代设计。
根据本发明的一较佳实施例，步骤一包括如下具体步骤：
步骤1.1：通过压气机一维设计，可以获得压气机进口物理流量
Figure 947364DEST_PATH_IMAGE001
，以及第i级单级设计工作点进口换算流量
Figure 872595DEST_PATH_IMAGE002
、总压比
Figure 283854DEST_PATH_IMAGE003
、总温比
Figure 610930DEST_PATH_IMAGE004
、单级效率
Figure 832964DEST_PATH_IMAGE005
、单级平均比热容比
Figure 979911DEST_PATH_IMAGE006
。
步骤1.2：求解图2所示的等出口换算流量线AB上的端点A的坐标。取进口换算流量偏置量
Figure 7910DEST_PATH_IMAGE007
（
Figure 376443DEST_PATH_IMAGE008
为定值，一般取1% ~ 2%），即
Figure 831695DEST_PATH_IMAGE009
………（1）
端点A的进口换算流量计算公式为
Figure 669201DEST_PATH_IMAGE027
………（2）
端点A的出口换算流量计算公式为
Figure 235312DEST_PATH_IMAGE028
………（3）
式（3）/（2）得到
Figure 458352DEST_PATH_IMAGE029
………（4）
由于已知条件
Figure 84505DEST_PATH_IMAGE013
Figure 143728DEST_PATH_IMAGE014
Figure 513530DEST_PATH_IMAGE015
Figure 669704DEST_PATH_IMAGE016
………（5）
将（5）中各式代入式（4），式（4）可变为
Figure 653710DEST_PATH_IMAGE017
………（4）
第i级单级效率EFF计算公式为
Figure 262546DEST_PATH_IMAGE030
………（5）
由“保持等熵效率不变”的条件可得
Figure 108142DEST_PATH_IMAGE031
………（6）
因此式（4）变为
Figure 118823DEST_PATH_IMAGE032
………（7）
联立式（4）、（7）可得关于
Figure 211413DEST_PATH_IMAGE021
、
Figure 494496DEST_PATH_IMAGE022
两个未知量的二元方程组，求解该方程组可获得
Figure 206100DEST_PATH_IMAGE021
，从而可以获得图2中等出口换算流量线的端点A的坐标位置（
Figure 274550DEST_PATH_IMAGE023
，
Figure 413407DEST_PATH_IMAGE021
）。
步骤1.3：采用步骤1.2相同的方法，求得当
Figure 731256DEST_PATH_IMAGE024
时，对应的等出口换算流量线的端点B的坐标位置（
Figure 433502DEST_PATH_IMAGE025
，
Figure 418775DEST_PATH_IMAGE026
），连接这两个端点A、B两点即可获得等出口换算流量线。
根据本发明的另一较佳实施例，步骤二包括如下具体步骤：
步骤2.1：第i级一维设计目标工作点为图2中的点T，但是压气机设计过程中由于各种原因，第i级的三维设计最终设计结果很难恰好位于T点。通过本方法，允许第i级设计工作点略偏离T点，但必须沿等出口换算流量线AB变化。沿该线变化，即使设计工作点出现偏离，该级的出口换算流量仍然保持不变，因而不会影响该级与下游各级的流量匹配关系，可以保持压气机一维气动设计的意图，保持多级压气机气动设计布局。通过该方法，实现了压气机单级设计从多级设计中剥离出来的目的，简化设计工作。
步骤2.2：具体开展第i级单级三维设计工作时，如图3所示，该级进口条件采用二维S2设计获得的二维总温、总压、气流角剖面设计结果，使三维设计与一维、二维设计结果保持一致。
具体地讲，在步骤2.1中，如图2所示，在压气机设计过程中，压气机的单级特性线（本文中“单级特性线”是指流量－压比单级特性线）通常由CFD（计算流体力学）软件计算获得的。该CFD软件计算获得的单级特性线与真实结果（“真实结果”是指最终通过压气机试验获得的压气机单级特性线）通常会存在一定偏差。对于某一种特定的CFD计算软件，该偏差通常为设计系统偏差（总是会偏大或偏小）
Figure 666217DEST_PATH_IMAGE033
。该设计系统偏差
Figure 471362DEST_PATH_IMAGE033
产生的原因，通常包括以下两个方面：
Figure 711719DEST_PATH_IMAGE034
（1）CFD计算软件产生的影响
Figure 551499DEST_PATH_IMAGE035
：
具体地讲，不同的CFD计算软件所使用的CFD算法、计算模型设置等是不同的，由此造成计算结果也会产生不同的偏差。通常，对于某一确定的CFD计算软件、确定的计算模型设置而言，其偏差影响
Figure 969842DEST_PATH_IMAGE035
是固定不变的。
（2）CFD计算所用几何模型与真实几何模型差异带来的影响
Figure 527863DEST_PATH_IMAGE036
：
具体地讲，真实的压气机几何模型通常是极为复杂的，在实际进行的压气机设计CFD计算过程中，几何模型是经过一定程度简化的，删减了某些几何特征（比如删减了转子封严篦齿几何特征、压气机叶片叶根倒圆几何特征等），这些被删减的几何特征对计算结果会产生一定影响，这就会造成计算结果与真实结果（试验结果）之间存在一定偏差，该偏差
Figure 306332DEST_PATH_IMAGE036
通常也就是系统偏差。
总地来说，在步骤2.1中，如图2所示，在压气机设计过程中，通过一定的技术积累，事先获取设计系统偏差
Figure 266197DEST_PATH_IMAGE033
。在进行单级压气机设计时，考虑该设计系统偏差
Figure 855442DEST_PATH_IMAGE033
后，将该
Figure 635179DEST_PATH_IMAGE033
沿着等出口换算流量线AB叠加到设计目标点T上，以便获得修正后的设计目标T’，也就是在进行压气机单级设计时，将CFD计算软件计算得到的单级特性匹配到T’点上，由此使得最终的单级压气机真实特性（即通过压气机试验获取的压气机单级特性）匹配到最终想要达到的设计目标T上，从而更好地实现设计期望。
综上所述，相较于现有技术而言，本发明多级压气机气动匹配设计方法基于等出口换算流量的原理，并以保持“等熵效率不变”为条件，计算获得经过单级设计工作点T的等出口换算流量线，从而在确保单级出口换算流量不变的条件下保持上下游各级之间的气动布局。在完成压气机整体气动布局设计的前提下，单级压气机设计基于等出口换算流量原理，每一级沿等出口换算流量线进行迭代设计，该级的设计工作点发生变化时，其出口换算流量保持不变，即下一级的进口换算流量保持不变，因而下一级流量特性不会发生变化。由于从前往后各级均采用该思路设计，因而压气机从前往后的流量特性保持不变，由此确保了压气机整体气动布局的稳定。
以上已详细描述了本发明的较佳实施例，但应理解到，在阅读了本发明的上述讲授内容之后，熟悉本领域的技术人员易于想到其它的优点和修改。因此，在其更宽泛的方面上来说，本发明并不局限于这里所示和所描述的具体细节和代表性实施例。因此，本领域技术人员能够将上述实施例的要素进行合理的组合或者改动，以便在不脱离如所附权利要求书及其等价物所限定的本发明总的发明概念的精神或范围的前提下作出各种修改。
本发明具体实施例中的符号说明：
Mass 压气机进口物理流量
Min 单级进口换算流量
Mout 单级出口换算流量
PR 单级进出口总压比
TR 单级进出口总温比
Ptin 单级进口总压
Ptout 单级出口总压
Ttin 单级进口总温
Ttout 单级出口总温
γ 比热容比
k 进口换算流量偏置量
target 设计工作点T对应的参数
A 端点A对应的参数
B 端点B对应的参数
T 单级设计工作点
T’ 修正后的单级设计工作点
∆ 设计系统偏差
Figure DEST_PATH_IMAGE037
由CFD计算软件产生的系统偏差
Figure 217339DEST_PATH_IMAGE036
由FD计算所用几何模型与真实几何模型差异带来的系统偏差。

Claims (6)

1.一种多级压气机气动匹配设计方法，其特征在于，所述方法包括以下步骤：
步骤一：开展第i级单级设计，基于一维设计结果，获取单级等出口换算流量线；
步骤二：通过设计系统偏差修正，沿等出口换算流量线，获取修正后的单级设计目标，进而开展第i级单级迭代设计。
2.根据权利要求1所述的多级压气机气动匹配设计方法，其特征在于，
在步骤一中，通过压气机一维设计，可以获得压气机进口物理流量
Figure 228486DEST_PATH_IMAGE001
，以及第i级单级设计工作点进口换算流量
Figure 524469DEST_PATH_IMAGE002
、总压比
Figure 872405DEST_PATH_IMAGE003
、总温比
Figure 985855DEST_PATH_IMAGE004
、单级效率
Figure 786451DEST_PATH_IMAGE005
、单级平均比热容比
Figure 10759DEST_PATH_IMAGE006
。
3.根据权利要求2所述的多级压气机气动匹配设计方法，其特征在于，
在步骤一中，求解所述等出口换算流量线上的其中一个端点的坐标，取进口换算流量偏置量
Figure 947622DEST_PATH_IMAGE007
，其中 k为定值，即
Figure 231973DEST_PATH_IMAGE008
………（1）
所述其中一个端点的进口换算流量计算公式为
Figure 519866DEST_PATH_IMAGE009
………（2）
所述其中一个端点的出口换算流量计算公式为
Figure 423231DEST_PATH_IMAGE010
………（3）
式（3）/（2）得到
Figure 604814DEST_PATH_IMAGE011
………（4）
由于已知条件
Figure 669853DEST_PATH_IMAGE012
Figure 569676DEST_PATH_IMAGE013
Figure 11153DEST_PATH_IMAGE014
Figure 188187DEST_PATH_IMAGE015
………（5）
将（5）中各式代入式（4），式（4）可变为
Figure 548761DEST_PATH_IMAGE016
………（4）
第i级单级效率EFF计算公式为
Figure 834684DEST_PATH_IMAGE017
………（5）
由“保持等熵效率不变”的条件可得
Figure 204486DEST_PATH_IMAGE018
………（6）
因此式（4）变为
Figure 970447DEST_PATH_IMAGE019
………（7）
联立式（4）、（7）可得关于
Figure 767502DEST_PATH_IMAGE020
、
Figure 251704DEST_PATH_IMAGE021
两个未知量的二元方程组，求解该方程组可获得
Figure 425197DEST_PATH_IMAGE020
，从而可以获得所述等出口换算流量线上的其中一个端点的坐标位置（
Figure 311244DEST_PATH_IMAGE022
，
Figure 154566DEST_PATH_IMAGE020
）。
4.根据权利要求3所述的多级压气机气动匹配设计方法，其特征在于，
在步骤一中，采用获得所述等出口换算流量线上的其中一个端点的坐标位置的相同的方法，求得当
Figure 250698DEST_PATH_IMAGE023
时，所述等出口换算流量线上的另一个端点的坐标位置（
Figure 837669DEST_PATH_IMAGE024
，
Figure 702856DEST_PATH_IMAGE025
），从而连接这两个端点即可获得所述等出口换算流量线。
5.根据权利要求1－4中任一项所述的多级压气机气动匹配设计方法，其特征在于，
在步骤二中，即使所述第i级设计工作点出现偏离，该级的出口换算流量仍然保持不变，因而不会影响该级与下游各级的流量匹配关系，从而保持多级压气机气动设计布局。
6.根据权利要求5所述的多级压气机气动匹配设计方法，其特征在于，
在步骤二中，开展第i级单级三维设计工作，该级进口条件采用二维S2设计获得的二维总温、总压、气流角剖面设计结果，使三维设计与一维、二维设计结果保持一致，即压气机各设计保持一致。
CN202110323072.XA 2021-03-26 2021-03-26 多级压气机气动匹配设计方法 Active CN112711923B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110323072.XA CN112711923B (zh) 2021-03-26 2021-03-26 多级压气机气动匹配设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110323072.XA CN112711923B (zh) 2021-03-26 2021-03-26 多级压气机气动匹配设计方法

Publications (2)

Publication Number Publication Date
CN112711923A true CN112711923A (zh) 2021-04-27
CN112711923B CN112711923B (zh) 2021-07-27

Family

ID=75550288

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110323072.XA Active CN112711923B (zh) 2021-03-26 2021-03-26 多级压气机气动匹配设计方法

Country Status (1)

Country Link
CN (1) CN112711923B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114444331A (zh) * 2022-04-07 2022-05-06 中国航发四川燃气涡轮研究院 一种多级轴流压气机的级特性匹配方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140223918A1 (en) * 2012-10-04 2014-08-14 Lightsail Energy, Inc. Compressed air energy system integrated with gas turbine
CN105930662A (zh) * 2016-04-26 2016-09-07 中国科学院工程热物理研究所 一种压气机低转速特性外推方法
CN106156436A (zh) * 2016-07-12 2016-11-23 中国航空工业集团公司沈阳发动机设计研究所 一种可调叶片角度分级调控的压气机建模方法
CN108229015A (zh) * 2017-12-30 2018-06-29 中国科学院工程热物理研究所 一种高空两级涡轮增压器变工况匹配设计方法
CN109165440A (zh) * 2018-08-22 2019-01-08 西北工业大学 一种轴流压气机全三维级间气动匹配优化方法
CN110489829A (zh) * 2019-07-31 2019-11-22 中国航发沈阳发动机研究所 一种基于流量特性的空气系统元件设计方法
CN110866312A (zh) * 2019-10-16 2020-03-06 西北工业大学 一种航空燃气涡轮发动机涡轮特性的格式转换方法
US10592622B1 (en) * 2014-10-24 2020-03-17 United States Of America As Represented By The Administrator Of National Aeronautics And Space Administration Assessment of conservatism level in an aircraft turbine engine design
CN111914362A (zh) * 2020-07-22 2020-11-10 中国航发沈阳发动机研究所 一种研发阶段涡扇发动机模型自适应方法
CN112081679A (zh) * 2020-06-24 2020-12-15 中国人民解放军陆军军事交通学院 变海拔柴油机二级可调增压系统与共轨燃油系统优化匹配方法
CN112528575A (zh) * 2020-12-22 2021-03-19 中国船舶重工集团公司第七0三研究所 一种船舶燃气轮机多级轴流压气机多维度协同设计方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140223918A1 (en) * 2012-10-04 2014-08-14 Lightsail Energy, Inc. Compressed air energy system integrated with gas turbine
US10592622B1 (en) * 2014-10-24 2020-03-17 United States Of America As Represented By The Administrator Of National Aeronautics And Space Administration Assessment of conservatism level in an aircraft turbine engine design
CN105930662A (zh) * 2016-04-26 2016-09-07 中国科学院工程热物理研究所 一种压气机低转速特性外推方法
CN106156436A (zh) * 2016-07-12 2016-11-23 中国航空工业集团公司沈阳发动机设计研究所 一种可调叶片角度分级调控的压气机建模方法
CN108229015A (zh) * 2017-12-30 2018-06-29 中国科学院工程热物理研究所 一种高空两级涡轮增压器变工况匹配设计方法
CN109165440A (zh) * 2018-08-22 2019-01-08 西北工业大学 一种轴流压气机全三维级间气动匹配优化方法
CN110489829A (zh) * 2019-07-31 2019-11-22 中国航发沈阳发动机研究所 一种基于流量特性的空气系统元件设计方法
CN110866312A (zh) * 2019-10-16 2020-03-06 西北工业大学 一种航空燃气涡轮发动机涡轮特性的格式转换方法
CN112081679A (zh) * 2020-06-24 2020-12-15 中国人民解放军陆军军事交通学院 变海拔柴油机二级可调增压系统与共轨燃油系统优化匹配方法
CN111914362A (zh) * 2020-07-22 2020-11-10 中国航发沈阳发动机研究所 一种研发阶段涡扇发动机模型自适应方法
CN112528575A (zh) * 2020-12-22 2021-03-19 中国船舶重工集团公司第七0三研究所 一种船舶燃气轮机多级轴流压气机多维度协同设计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YUDONG XU ET AL: "Influence of the Number of Blades of Miniature Centrifugal Impeller on Compressor Performance", 《2020 9TH INTERNATIONAL CONFERENCE ON POWER SCIENCE AND ENGINEERING (ICPSE) 》 *
吴学辉 等: "涡轴发动机总体改型的方法研究", 《机械科学与技术》 *
谭旭刚: "某型多级轴流压气机试验件结构设计及试验研究", 《中国优秀硕士学位论文全文数据库电子期刊 工程科技II辑》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114444331A (zh) * 2022-04-07 2022-05-06 中国航发四川燃气涡轮研究院 一种多级轴流压气机的级特性匹配方法

Also Published As

Publication number Publication date
CN112711923B (zh) 2021-07-27

Similar Documents

Publication Publication Date Title
CN109684597B (zh) 一种用于压气机全工况特性拓展的方法
CN112287580B (zh) 一种基于全三维数值模拟的轴流压气机喘振边界计算方法
US20080301939A1 (en) Gas turbine and manufacturing process of gas turbine
CN112685968B (zh) 基于空间负荷定制思想的轴流压气机气动设计方法
CN113513373A (zh) 一种变几何涡轮一维气动设计方法
US7941300B1 (en) Process for the design of an airfoil
CN112711923B (zh) 多级压气机气动匹配设计方法
CN111027148B (zh) 损失落后角模型自动标定及工业轴流压缩机性能计算方法
Marchukov et al. Improving of the working process of axial compressors of gas turbine engines by using an optimization method
CN113283198A (zh) 优化压气机机匣处理提高稳定裕度的方法、系统、终端
Mojaddam et al. Design and optimization of meridional profiles for the impeller of centrifugal compressors
Sadagopan et al. A design strategy for a 6: 1 supersonic mixed-flow compressor stage
CN113883093B (zh) 一种低反力度压气机叶片设计方法、动叶及压气机
Goryachkin et al. Three-stage low pressure compressor modernization by means of optimization methods
CN117494566A (zh) 一种轴流压气机一维模型建模及修正方法
Wang et al. A forward-curved blade centrifugal compressor for anode recirculation in proton exchange membrane fuel cells
Liang et al. The recent progresses in industrial centrifugal compressor designs
CN112395702B (zh) 轴流压气机叶片的优化设计方法和系统
Li et al. The optimization of a centrifugal impeller based on a new multi-objective evolutionary strategy
CN114239185A (zh) 一种燃气轮机多级轴流压气机的设计方法及装置
CN114491417A (zh) 一种基于cdfs模态变化性能的一维输入修正方法
Avşar et al. Aerodynamic Design and Performance Optimization of a Centrifugal Fan Impeller
Sedunin et al. Redesign of an axial compressor with mass flow reduction of 30%
Krivosheev et al. Optimization of the Stage Number and Parameter Distribution in the Flow Passage at GTE Compressor and Turbine Design
Boyko et al. Devising a method for calculating the turboshaft gas turbine engine performance involving a blade-by-blade description of the multi-stage compressor in a two-dimensional setting

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant