CN108052694B - 隐身性加力燃烧室内锥体设计方法 - Google Patents

隐身性加力燃烧室内锥体设计方法 Download PDF

Info

Publication number
CN108052694B
CN108052694B CN201711130079.XA CN201711130079A CN108052694B CN 108052694 B CN108052694 B CN 108052694B CN 201711130079 A CN201711130079 A CN 201711130079A CN 108052694 B CN108052694 B CN 108052694B
Authority
CN
China
Prior art keywords
section
afterburner
area
calculating
cone
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.)
Active
Application number
CN201711130079.XA
Other languages
English (en)
Other versions
CN108052694A (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 Shenyang Engine Research Institute
Original Assignee
AECC Shenyang Engine Research Institute
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 Shenyang Engine Research Institute filed Critical AECC Shenyang Engine Research Institute
Priority to CN201711130079.XA priority Critical patent/CN108052694B/zh
Publication of CN108052694A publication Critical patent/CN108052694A/zh
Application granted granted Critical
Publication of CN108052694B publication Critical patent/CN108052694B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Landscapes

  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Ignition Installations For Internal Combustion Engines (AREA)
  • Testing Of Engines (AREA)

Abstract

本发明涉及航空发动机技术领域,具体提供了隐身性加力燃烧室内锥体设计方法,首先获得内锥体的加力燃烧室内涵进口面积、出口面积、进口总压、总温、流量以及扩张段长度,然后计算当量扩张半角和气动参数,得到进口截面速度系数和气流静压,再计算出气流速度系数减小到设定值时对应的流道面积,并根据流道面积、内涵进口截面面积及当量扩张半角确定首段扩张长度,确定分段点,之后计算内涵进口截面和分段点截面各自的静压并计算出每等份的截面面积,进而得到首段型面,最后由分段点和内涵出口截面内边界点确定直锥型面,得到内椎体型面。该内锥体设计方法在继承了传统方法包含低损失的减速扩压设计的同时,又可实现低雷达可探测性要求。

Description

隐身性加力燃烧室内锥体设计方法
技术领域
本发明涉及航空发动机技术领域,特别涉及隐身性加力燃烧室内锥体设计方法。
背景技术
加力燃烧室进口气流速度大,必须先进行扩压减速,才能开展组织燃烧设计。内锥体是对加力燃烧室内气流扩压的主要构件之一,其设计优劣将直接影响加力燃烧室性能。
为了实现低损失的扩压效果,传统的设计方法得到的加力内锥体常采用如图1所示的全锥的形式,或者如图2所示的截锥形式。全锥尤其是截锥的镜面反射效应,使后视雷达可探测性大幅增强,无法满足战斗机、发动机的隐身要求。根据隐身系统研究表明,低雷达可探测性的内锥体的型面,应避免采用类似镜面或球面的截锥及全锥,而且对直锥与加力燃烧室气流轴心线的夹角也有一定要求。因此迫切需要开发一款新的加力燃烧室内锥体设计方法,在内锥体设计过程,可同时满足扩压要求和隐身要求。
发明内容
为克服上述现有技术存在的至少一种缺陷,本发明提供了隐身性加力燃烧室内锥体设计方法,包括如下步骤:
步骤一,设计输入准备:获得设计内锥体的边界条件和设计指标,所述边界条件和设计指标包括加力燃烧室内涵进口面积、出口面积、进口总压、总温、流量,以及扩张段长度L;
步骤二,当量扩张半角计算:根据步骤一中获得的加力燃烧室内涵进口面积、出口面积和扩张段长度L计算当量扩张半角;
步骤三,气动参数计算:根据步骤一中获得的加力燃烧室内涵进口总压、总温、流量,采用一维定常连续方程计算进口截面速度系数λ,如公式(1)所示:
Figure GDA0002822148190000021
其中Qm为质量流量,C为常数,Pt为总压,Tt为总温,A为流道面积,q(λ)为流量函数,ε(λ)为静密度和总密度的比值,k为气体比热比,R为气体常数且R=287J/(kg﹒K);
并采用总静压换算方程计算在速度系数λ下的气流静压P,如公式(2)所示:
Figure GDA0002822148190000022
步骤四,分段点确定:根据步骤三中的公式(1),计算出气流速度系数减小到设定值时对应的流道面积A1,根据流道面积A1、内涵进口截面面积及当量扩张半角确定首段扩张长度L1,与首段扩张长度L1对应的点为分段点;
步骤五,首段设计:将首段扩张长度L1分成若干等份,计算内涵进口截面和分段点截面各自的静压,然后按每等份静压增长梯度相同的原则计算出每等份的截面面积,再根据截面面积折算对应的内锥体径向坐标,最后以平滑曲线连接每等份的径向坐标点得到首段型面;
步骤六,直锥设计:由分段点和内涵出口截面内边界点确定直锥型面,判断直锥与加力燃烧室气流轴心线的夹角是否满足隐身要求的角度,若不符合要求则返回步骤四并重新选择分段点,重新选取后判断对应的当量扩张半角是否在气动要求范围内;若符合要求,则将分段点与内边界点确定的直线,延伸至加力燃烧室轴心线,形成直锥型面,得到内椎体型面。
优选的,步骤一中的所述边界条件和设计指标还包括合流环型面和前安装边几何位置。
优选的,步骤二中计算当量扩张半角时将所述加力燃烧室内涵进口面积和出口面积分别折算成完整的圆形。
优选的,步骤四中的设定值为0.27。
本发明提供的隐身性加力燃烧室内锥体设计方法,在继承了传统方法包含低损失的减速扩压设计的同时,又可实现低雷达可探测性要求。
附图说明
以下参考附图描述的实施例是示例性的,旨在用于解释和说明本发明,而不能理解为对本发明的保护范围的限制。
图1是现有技术中加力燃烧室内全锥形式的示意图;
图2是现有技术中加力燃烧室内截锥形式的示意图;
图3是内椎体的几何边界示意图;
图4是计算当量扩张半角的示意图;
图5是确定分段点的示意图;
图6是获得的隐身性加力燃烧室内锥体的型面示意图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行更加详细的描述。
需要说明的是:在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明一部分实施例,而不是全部的实施例,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,均仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明保护范围的限制。
在本文中,“示意性”表示“充当实例、例子或说明”,不应将在本文中被描述为“示意性”的任何图示、实施方式解释为一种更优选的或更具优点的技术方案。
本发明提供了隐身性加力燃烧室内锥体设计方法,包括如下步骤:
步骤一,设计输入准备:获得设计内锥体的边界条件和设计指标,所述边界条件和设计指标包括加力燃烧室内涵进口面积、出口面积、进口总压、总温、流量,以及扩张段长度L。本实施例中,边界条件和设计指标还包括合流环型面和前安装边几何位置。
步骤二,当量扩张半角计算:如图3及图4所示,根据步骤一中获得的加力燃烧室内涵进口面积、出口面积(即混合截面内涵面积)和扩张段长度L计算当量扩张半角。计算当量扩张半角时将所述加力燃烧室内涵进口面积和出口面积分别折算成完整的圆形,两个圆之间的距离与扩压段长度相同,即内涵进、出口间距离。如图4所示,本实施例中当量扩张半角为17°。
步骤三,气动参数计算:根据步骤一中获得的加力燃烧室内涵进口总压、总温、流量,采用一维定常(等熵)连续方程计算进口截面速度系数λ,如公式(1)所示:
Figure GDA0002822148190000051
其中Qm为质量流量(kg/s),C为常数,对于空气,C=0.04042,燃气k=1.33时C=0.03969,Pt为总压(N/m2),Tt为总温(K),A为流道面积(m2),q(λ)为流量函数,ε(λ)为静密度和总密度的比值,k为气体比热比,空气为1.4,燃气1.33~1.36左右,R为气体常数且R=287J/(kg﹒K);
并采用总静压换算方程计算在速度系数λ下对应流动截面的气流静压P(N/m2),如公式(2)所示:
Figure GDA0002822148190000052
步骤四,分段点确定:如图5所示,根据步骤三中的公式(1),计算出气流速度系数减小到设定值时对应的流道面积A1,本实施例中,设定值为0.27,根据流道面积A1、内涵进口截面面积及当量扩张半角确定首段扩张长度L1,与首段扩张长度L1对应的点为分段点;此处所选取的当量扩张半角,可参照内涵进口截面速度系数的大小适当调整,当内涵进口截面速度系数偏大时,所选当量扩张半角向偏小方向修正(当量扩张半角根据内涵进口截面速度系数决定),当量扩张半角选取越小,则首段扩张长度L1将越长。
步骤五,首段设计:如图5所示,将首段扩张长度L1分成若干等份,计算内涵进口截面和分段点截面各自的静压,然后按每等份静压增长梯度相同的原则计算出每等份的截面面积,再根据截面面积折算对应的内锥体径向坐标,最后以平滑曲线连接每等份的径向坐标点得到首段型面。
步骤六,直锥设计:如图6所示,由分段点和内涵出口截面内边界点确定直锥型面,判断直锥与加力燃烧室气流轴心线的夹角是否满足隐身要求的角度,若不符合要求则返回步骤四并重新选择分段点,重新选取后判断对应的当量扩张半角是否在气动要求范围内,因为重新选择分段点后,当量扩张半角会随之变化,因此需要重新判断当量扩张半角是否在气动要求范围内;若符合要求,则将分段点与内边界点确定的直线,延伸至加力燃烧室轴心线,形成直锥型面,得到内椎体型面。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (4)

1.隐身性加力燃烧室内锥体设计方法,其特征在于,包括如下步骤:
步骤一,设计输入准备:获得设计内锥体的边界条件和设计指标,所述边界条件和设计指标包括加力燃烧室内涵进口面积、出口面积、进口总压、总温、流量,以及扩张段长度L;
步骤二,当量扩张半角计算:根据步骤一中获得的加力燃烧室内涵进口面积、出口面积和扩张段长度L计算当量扩张半角;
步骤三,气动参数计算:根据步骤一中获得的加力燃烧室内涵进口总压、总温、流量,采用一维定常连续方程计算进口截面速度系数λ,如公式(1)所示:
Figure FDA0002822148180000011
其中Qm为质量流量,C为常数,Pt为总压,Tt为总温,A为流道面积,q(λ)为流量函数,ε(λ)为静密度和总密度的比值,k为气体比热比,R为气体常数且R=287J/(kg﹒K);
并采用总静压换算方程计算在速度系数λ下的气流静压P,如公式(2)所示:
Figure FDA0002822148180000012
步骤四,分段点确定:根据步骤三中的公式(1),计算出气流速度系数减小到设定值时对应的流道面积A1,根据流道面积A1、内涵进口截面面积及当量扩张半角确定首段扩张长度L1,与首段扩张长度L1对应的点为分段点;
步骤五,首段设计:将首段扩张长度L1分成若干等份,计算内涵进口截面和分段点截面各自的静压,然后按每等份静压增长梯度相同的原则计算出每等份的截面面积,再根据截面面积折算对应的内锥体径向坐标,最后以平滑曲线连接每等份的径向坐标点得到首段型面;
步骤六,直锥设计:由分段点和内涵出口截面内边界点确定直锥型面,判断直锥与加力燃烧室气流轴心线的夹角是否满足隐身要求的角度,若不符合要求则返回步骤四并重新选择分段点,重新选取后判断对应的当量扩张半角是否在气动要求范围内;若符合要求,则将分段点与内边界点确定的直线,延伸至加力燃烧室轴心线,形成直锥型面,得到内椎体型面。
2.根据权利要求1所述的隐身性加力燃烧室内锥体设计方法,其特征在于,步骤一中的所述边界条件和设计指标还包括合流环型面和前安装边几何位置。
3.根据权利要求1所述的隐身性加力燃烧室内锥体设计方法,其特征在于,步骤二中计算当量扩张半角时将所述加力燃烧室内涵进口面积和出口面积分别折算成完整的圆形。
4.根据权利要求1所述的隐身性加力燃烧室内锥体设计方法,其特征在于,步骤四中的设定值为0.27。
CN201711130079.XA 2017-11-15 2017-11-15 隐身性加力燃烧室内锥体设计方法 Active CN108052694B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711130079.XA CN108052694B (zh) 2017-11-15 2017-11-15 隐身性加力燃烧室内锥体设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711130079.XA CN108052694B (zh) 2017-11-15 2017-11-15 隐身性加力燃烧室内锥体设计方法

Publications (2)

Publication Number Publication Date
CN108052694A CN108052694A (zh) 2018-05-18
CN108052694B true CN108052694B (zh) 2021-01-19

Family

ID=62120215

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711130079.XA Active CN108052694B (zh) 2017-11-15 2017-11-15 隐身性加力燃烧室内锥体设计方法

Country Status (1)

Country Link
CN (1) CN108052694B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113361081B (zh) * 2021-05-23 2022-07-01 中国航发沈阳发动机研究所 一种带导流支板的加力燃烧室外涵气动面积确定方法
CN114021262B (zh) * 2021-11-15 2022-09-20 中国航发沈阳发动机研究所 一种三通道前置扩压器设计方法
CN115614173B (zh) * 2022-12-19 2023-02-28 中国航发沈阳发动机研究所 一种径向搭接轴向插接内锥体安装边结构

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015073498A1 (en) * 2013-11-12 2015-05-21 Aero Systems Engineering, Inc. Systems and methods for improved accuracy
CN105373637A (zh) * 2014-08-26 2016-03-02 中国航空工业集团公司沈阳发动机设计研究所 一种波瓣混合器流量系数的数值计算方法
CN206514768U (zh) * 2017-02-16 2017-09-22 中国航发沈阳发动机研究所 一种加力燃烧室隔热屏冷却通道高度测量工具

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015073498A1 (en) * 2013-11-12 2015-05-21 Aero Systems Engineering, Inc. Systems and methods for improved accuracy
CN105373637A (zh) * 2014-08-26 2016-03-02 中国航空工业集团公司沈阳发动机设计研究所 一种波瓣混合器流量系数的数值计算方法
CN206514768U (zh) * 2017-02-16 2017-09-22 中国航发沈阳发动机研究所 一种加力燃烧室隔热屏冷却通道高度测量工具

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
The Analysis and Diagnosis of Vibration Fault of a Certain Type of Aero-Engine;Zhengyi Jiang等;《Applied Mechanics and Materials》;20140930;全文 *
整流支板和火焰稳定器的一体化设计加力燃烧室性能的数值模拟;李锋 等;《航空发动机》;20121030;全文 *

Also Published As

Publication number Publication date
CN108052694A (zh) 2018-05-18

Similar Documents

Publication Publication Date Title
CN108052694B (zh) 隐身性加力燃烧室内锥体设计方法
US10731556B2 (en) Isentropic compression inlet for supersonic aircraft
CA2379091C (en) Supersonic external-compression diffuser and method for designing same
EP1643083B1 (en) Coupled parametric design of flow control and duct shape
Lucas et al. Effect of BLI–type inlet distortion on turbofan engine performance
Otto et al. Inward-turning streamline-traced inlet design method for low-boom, low-drag applications
US20220145777A1 (en) Gas turbine engine nacelle and method of designing same
CN103987949A (zh) 一种涡轮发动机的缩放喷嘴
Vázquez et al. The effect of airfoil clocking on efficiency and noise of low pressure turbines
Migdal Supersonic annular nozzles
US20210277802A1 (en) Fan case for a gas turbine engine
CA2898663C (en) Gas turbine engine exhaust fluid passage duct
Tomita et al. Investigation on characteristics of conventional-nozzle-based altitude compensating nozzles by cold-flow tests
Slater Design factors for two-dimensional, external-compression supersonic inlets
US20150101338A1 (en) Exhaust centerbody for a turbine engine
US20160290645A1 (en) Axisymmetric offset of three-dimensional contoured endwalls
Conners et al. Wind tunnel testing of an axisymmetric isentropic relaxed external compression inlet at Mach 1.97 design speed
Flegel-McVetta et al. Aerodynamic measurements of a variable-speed power-turbine blade section in a transonic turbine cascade at low inlet turbulence
Slater Supersonic External-Compression Inlets for Mach 1.4 to 2.0
Abbott et al. Computational study of the aerodynamic performance of three-dimensional subsonic inlets
C¸ evik et al. A robust design optimization of a mixed-flow compressor impeller
Jegannath et al. Comprehensive Studies on the Aerodynamic Characteristics of Subsonic Scarf Inlets
Mikheev et al. Shaping a SPRE multinozzle cluster with out-of-round nozzles
Colville et al. Inverse Design of Hypersonic Star Inlets

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