CN114510798B - 一种航空发动机气动轴向力防错分析方法 - Google Patents
一种航空发动机气动轴向力防错分析方法 Download PDFInfo
- Publication number
- CN114510798B CN114510798B CN202210089680.3A CN202210089680A CN114510798B CN 114510798 B CN114510798 B CN 114510798B CN 202210089680 A CN202210089680 A CN 202210089680A CN 114510798 B CN114510798 B CN 114510798B
- Authority
- CN
- China
- Prior art keywords
- axial force
- pneumatic
- unit body
- interface
- rotor
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
- 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)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本申请提供了一种航空发动机气动轴向力防错分析方法,所述方法包括:步骤一、零组件级气动轴向力计算;步骤二、零组件级界面参数一致性确认;步骤三、单元体级气动轴向力计算;步骤四、单元体气动轴向力一致性确认;步骤五、单元体间界面参数一致性确认;步骤六、整机级气动轴向力计算;步骤七、整机级气动轴向力一致性确认。本申请的方法通过采用全程双向闭环工作模式,以零组件/单元体界面两侧结构尺寸和气动参数取值一致性为目标,以两种算法的气动轴向力一致性为检验准则,对整机气动轴向力进行逐级计算、逐级校验,确保消除整机轴向力计算过程中的人为失误。
Description
技术领域
本申请属于航空发动机设计技术领域,特别涉及一种航空发动机气动轴向力防错分析方法。
背景技术
航空发动机工作过程中,从风扇前进气端到喷管排气端,气流压力、流量、流速和结构特征都在不断变化,产生了各零部件气动轴向力的分布。气动轴向力可分为外力和内力,其中,外力即发动机推力,产生于发动机进口、出口和外壁面气流的动量和压力变化;内力来源于发动机零部件之间,以及零部件与气体之间的相互作用力。各零部件所受气动轴向力是外力和内力共同作用的结果。从强度及寿命分析角度,需要考虑气动轴向力载荷,其计算结果的正确性至关重要。
如0所示,在发动机整机气动轴向力计算过程中,由于发动机工作过程的复杂性,需要处理多种数据,包括典型装配剖面(转子叶尖、篦齿封严等)和转静子安装边位置的结构尺寸数据、部件性能参数、空气***腔压数据等。这些数据又涵盖总体结构、部件结构和性能、空气***、滑油***多个专业。由于输入数据来源广泛、数据格式多样、数据量庞大,人为环节多、处理流程繁琐,在现有技术中,经常出现各专业计算得到的气动轴向力相加后与整体气动轴向力存在较大误差的情况。为避免上述问题,需要建立一种符合发动机工作原理、简便可行的操作流程或方法,确保发动机气动轴向力计算的准确性,提高工作效率。
发明内容
本申请的目的是提供了一种航空发动机气动轴向力防错分析方法,以解决或减轻背景技术中的至少一个问题。
本申请的技术方案是:一种航空发动机气动轴向力防错分析方法,所述方法包括:
步骤一、确定零组件轴向力防错分析所需中间过程数据,根据发动机各零组件的气动轴向力计算方法进行零组件级气动轴向力计算;
步骤二:对发动机零组件中各转子-静子、转子-转子、静子-静子界面进行检查,确保界面两侧选用的结构和气动参数符合一致性要求;
步骤三、将发动机按照部件划分为若干个由一个或多个零组件构成的单元体,根据各单元体的结构及气动参数分别得到各单元体轴向力基准值;
步骤四、将各单元体内一个或多个零组件的轴向力进行叠加得到各单元体轴向力合成值,判断各单元体轴向力基准值与合成值的差值是否满足要求,若满足要求,则判断单元体轴向力计算中人为失误消除;若不满足要求,则返回步骤二、步骤三和步骤四,复查相应单元体轴向力基准值、合成值的计算过程,排除人为失误;
步骤五、完成各单元体轴向力基准值与单元体轴向力合成值校验后,确保各单元体之间的界面参数符合一致性要求:
步骤六:将发动机整体视为一个气动控制体,根据发动机进气端、排气端、外壁面气动和结构参数计算发动机推力基准值;
步骤七:按照发动机传力路径,将步骤一中各零组件气动轴向力叠加至主安装节所在的机匣,并叠加该机匣自身的轴向力,得到发动机推力合成值;
若发动机推力合成值与基准值的误差满足要求,则判断所有过程已消除人为失误;若发动机推力合成值与基准值的误差不满足要求,则返回步骤四重新检查,定位并修正相应人为失误后,重新完成步骤五、步骤六、步骤七。
在本申请优选实施方式中,所述零组件中各转子-静子、转子-转子、静子-静子界面两侧选用的结构和气动参数一致性要求为:
对于转子-静子界面:
1)转子侧与静子侧选用的结构尺寸一致,均为转子侧结构尺寸;
2)转子侧与静子侧气动参数一致,均为转子侧气动参数;
3)转子侧与静子侧轴向力受力大小相等、方向相反。
对于转子-转子界面:
1)两侧结构尺寸一致,均为其中一侧的结构尺寸;
2)两侧气动参数一致,均为其中一侧的气动参数;
3)两侧轴向力受力大小相等、方向相反。
对于静子-静子界面:
1)两侧结构尺寸一致,均为其中一侧的结构尺寸;
2)两侧气动参数一致,均为其中一侧的气动参数;
3)两侧轴向力受力大小相等、方向相反。
在本申请优选实施方式中,各单元体轴向力基准值与合成值的差值与基准值的比值小于1%,则满足要求;反之,则不满足。
在本申请优选实施方式中,当单元体轴向力基准值与合成值的差值不满足要求,且单元体内多个零组件所形成的安装边数量多于3个时,对单元体进行再次细化,直到所述单元体内的安装边数量为2~3个。
在本申请优选实施方式中,各单元体之间的界面参数一致性要求判据包括:
1)两个单元体界面两侧的结构尺寸一致;
2)两个单元体界面两侧的气动参数一致;
3)两个单元体界面两侧的轴向力方向相反;
如果满足上述3项判据,则判断相邻单元体界面参数一致;如果存在不满足的情况,则修正相应位置结构或气动参数,重新按顺序开展防错分析。
在本申请优选实施方式中,所述发动机推力合成值与基准值的误差满足要求为:
发动机推力合成值与基准值的误差与基准值的比例小于1%。
本申请的方法能够完全杜绝整机轴向力计算过程中的人为失误,并且大大缩短误差源定位时间周期,大幅提高工作效率。
附图说明
为了更清楚地说明本申请提供的技术方案,下面将对附图作简单地介绍。显而易见地,下面描述的附图仅仅是本申请的一些实施例。
图1为现有技术中典型的航空发动机结构示意图。
图2为本申请的航空发动机气动轴向力防错分析流程图。
图3为本申请中的零组件级“转子-静子”界面参数选取方法示意图。
图4为本申请中的单元体划分过程示意图。
图5为本申请中的单元体轴向力合成值计算过程示意图。
图6为本申请中的中介机匣受力示意图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行更加详细的描述。
为了克服现有技术中,航空发动机整机气动轴向力计算中,输入数据量、,来源庞杂、处理环节多,容易产生人为失误,且缺乏有效的误差源定位,使得误差源排查工作陷入无休止的重复计算过程,最终也难以达到理想的计算精度到问题,本申请中提出了一种将轴向力计算与误差源排查全程融合的防错分析方法,以消除人为失误,提高发动机气动轴向力工作效率和计算精度。
本申请的防错分析方法流程图如图2所示,本申请中按照零组件级→单元体级→整机级的顺序,采用全程双向闭环工作模式,以零组件/单元体界面两侧结构尺寸和气动参数取值一致性为目标,以两种算法的气动轴向力一致性为检验准则,对整机气动轴向力进行逐级计算、逐级校验,确保消除整机轴向力计算过程中的人为失误。
具体的,本申请的航空发动机气动轴向力防错分析方法包括如下步骤:
一、零组件级轴向力防错分析
步骤一:梳理设计输入,确定各零组件气动轴向力防错分析中所需的中间过程数据,如表1所示。
表1 零组件轴向力防错分析所需中间过程数据
序号 | 数据项 |
1 | 叶片排进口叶尖直径、出口叶尖直径 |
2 | 叶片排进口叶根直径、出口叶根直径 |
3 | 叶片排进口静压、出口静压 |
4 | 气体流量 |
5 | 叶片排进口气流轴向速度、出口气流轴向速度 |
6 | 内腔外环直径和内环直径 |
7 | 内腔静压 |
8 | 轴向力 |
根据各零组件的相关计算方法计算零组件气动轴向力,各零组件的气动轴向力计算方法为常规方法,本处不再赘述。
步骤二:对发动机零组件级中各转子-静子、转子-转子、静子-静子界面进行检查,确保界面两侧选用的结构和气动参数符合以下条件:
对于转子-静子界面:
1)转子侧与静子侧选用的结构尺寸(环形作用面的内、外环直径)一致,均为转子侧结构尺寸;
2)转子侧与静子侧气动参数(静压、流量、流速)一致,均为转子侧气动参数;
3)转子侧与静子侧轴向力受力大小相等、方向相反。
对于转子-转子界面:
1)两侧结构尺寸(环形作用面的内、外环直径)一致,均为其中一侧的结构尺寸;
2)两侧气动参数(静压、流量、流速)一致,均为其中一侧的气动参数;
3)两侧轴向力受力大小相等、方向相反。
对于静子-静子界面:
1)两侧结构尺寸(环形作用面内、外环直径)一致,均为其中一侧的结构尺寸;
2)两侧气动参数(静压、流量、流速)一致,均为其中一侧的气动参数;
3)两侧轴向力受力大小相等、方向相反。
如图3所示,以某发动机中左侧的转子叶片与右侧的导向叶片(静子叶片)形成的界面为例,对零组件级“转子-静子”界面参数的防错分析及实施步骤进行举例说明:
首先,根据发动机各零组件的气动轴向力计算方法进行零组件级气动轴向力计算。
1)转子侧轴向力计算输入原始参数及轴向力检查,原始参数项目包括:转子叶片排出口叶尖直径(DRBK)、转子叶片排出口叶根直径(DRBH)、转子叶片排出口静压(PRB)、转子叶片排出口流量(WRB)、转子叶片排出口流速(VRB)。
按照发动机原理,转子叶片排出口轴向力为:
FRB=π·PRB·(DRBK^2-DRBH^2)/4+WRB·VRB (公式1)。
2)静子侧轴向力计算输入原始参数及轴向力检查,原始参数项目包括:静子叶片排进口叶尖直径(DSAK)、静子叶片排进口叶根直径(DSAH)、静子叶片排进口静压(PSA)、静子叶片排进口流量(WSA)、静子叶片排进口流速(VSA)。
按照发动机原理,静子叶片排出口轴向力为:
FSA=π·PSA·(DSAK^2-DSAH^2)/4+WSA·VSA (公式2)。
静子叶片排出口轴向力FSA受力对象为静子叶片排,受力方向为逆航向向后。
之后,对发动机零组件中各转子-静子界面进行检查,确保界面两侧选用的结构和气动参数符合一致性要求。
如本申请中零组件级轴向力防错分析方法所阐述的,应对静子侧(此处为“静子叶片排进口”)轴向力计算输入参数进行调整,即:静子叶片排进口叶尖直径应由“DSAK”更正为“DRBK”,静子叶片排进口叶根直径应由“DSAH”更正为“DRBH”,静子叶片排进口静压应由“PSA”更正为“PRB”,静子叶片排进口流量应由“WSA”更正为“WRB”,静子叶片排进口流速应由“VSA”更正为“VRB”。
因此,静子叶片排出口轴向力应由公式2中“FSA=π·PSA·(DSAK^2-DSAH^2)/4+WSA·VSA”,更正为静子叶片排出口轴向力为:
FSA’=π·PRB·(DRBK^2-DRBH^2)/4+WRB·VRB (公式3)
受力对象和受力方向不变:受力对象为静子叶片排,受力方向为逆航向向后。
在数值上,FRB=FSA’。
3)受力对象及受力方向检查,
转子叶片排出口轴向力FRB受力对象为转子叶片排,受力方向为顺航向向前;
更正后的静子叶片排出口轴向力FSA’受力对象为静子叶片排,受力方向为逆航向向后。
二、单元体级轴向力防错分析
步骤三:将发动机按部件划分为若干个较大的气动控制体(又称单元体),分别计算各单元体轴向力基准值,其中单元体进口、外壁面和出口的结构尺寸、气动参数及数据处理方式应分别与步骤一中相应零组件轴向力计算用数据保持一致。
如图4所示为本申请该实施例中将发动机划分为六个单元体,单元体划分、单元体轴向力基准值的确定过程:
1)“单元体1”为图4中“D(No.1轴承)→A(进气机匣进口安装边)→B(风扇机匣出口安装边)→C(No.2轴承)→D(No.1轴承)”各点连线(对于安装边,为相连接两个零件之间的结构界面;对于其余部分,为零件与气体之间的气动界面)所形成的虚拟封闭区域(风扇区域气动控制体)。
在发动机工作条件下(不同计算工况),在步骤一的中间过程数据中,沿着该气动控制体的轮廓线(构成虚拟封闭区域的各连线)选取结构参数(直径)和气动参数(包括静压、流量、流速),计算“单元体1”轴向力。
“单元体1”的气动轴向力基准值由以下4项组成:风扇进口(包括流路、帽罩,轴向力计为FFAN_INT)、风扇机匣外壁面(机匣暴露于环境大气部分的壁面,轴向力计为FFAN_CAS)、风扇出口(包括流路、空气***腔、No.2轴承腔,轴向力计为FFAN_OUT)、风扇盘心腔(包括空气***腔、No.1和No.2轴承腔,轴向力计为FFAN_COR),“单元体1”的气动轴向力为:
FITG_FAN=FFAN_INT+FFAN_CAS+FFAN_OUT+FFAN_COR (公式4)
2)“单元体2”为图4中“C(No.2轴承)→B(中介机匣进口安装边)→H(中介机匣出口外环安装边)→F(中介机匣出口内环安装边)→E(No.3轴承)→C(No.2轴承)”各点“连线”(对于安装边,为相连接两个零件之间的结构界面;对于其余部分,为零件与气体之间的气动界面)所形成的虚拟封闭区域(中介机匣气动控制体)。
在发动机工作条件下(不同计算工况),在步骤一的中间过程数据中,沿着该气动控制体的轮廓线(构成虚拟封闭区域的各连线)选取结构参数(直径)和气动参数(包括静压、流量、流速),计算“单元体2”轴向力。
“单元体2”的气动轴向力基准值由以下4项组成:中介机匣进口(包括流路、空气***腔、No.2轴承腔,轴向力计为FINC_INT)、中介机匣外壁面(机匣暴露于环境大气部分的壁面,轴向力计为FINC_CAS)、中介机匣出口(包括外涵道进口流路、压气机进口流路、压气机转子前端空气***腔,轴向力计为FINC_OUT)、中介机匣内环(No.2和No.3轴承腔,轴向力计为FINC_COR),“单元体2”的气动轴向力为:
FINC=FINC_INT+FINC_CAS+FINC_OUT+FINC_COR (公式5)
3)“单元体3”为图4中“E(No.3轴承)→F(压气机机匣进口安装边)→L(高压涡轮机匣出口安装边)→M(No.4轴承)→E(No.3轴承)”各点连线(对于安装边,为相连接两个零件之间的结构界面;对于其余部分,为零件与气体之间的气动界面)所形成的虚拟封闭区域(核心机气动控制体)。
在发动机工作条件下(不同计算工况),在步骤一的中间过程数据中,沿着该气动控制体的轮廓线(构成虚拟封闭区域的各连线)选取结构参数(直径)和气动参数(包括静压、流量、流速),计算“单元体3”轴向力。
“单元体3”的气动轴向力基准值由以下4项组成:压气机进口(包括流路、空气***腔、No.3轴承腔,轴向力计为FEGC_INT)、压气机机匣-主燃烧室机匣-高压涡轮机匣外壁面(机匣暴露于外涵流路部分的壁面,轴向力计为FEGC_CAS)、高压涡轮出口(包括高压涡轮出口流路、高压涡轮转子后端空气***腔、No.4轴承腔,轴向力计为FEGC_OUT)、核心机内环(No.3和No.4轴承腔,以及压气机盘、主燃烧室内环、高压涡轮盘各空气***腔,轴向力计为FEGC_COR),“单元体3”的气动轴向力为:
FITG_EGC=FEGC_INT+FEGC_CAS+FEGC_OUT+FEGC_COR (公式6)
4)“单元体4”为图4中“M(No.4轴承)→L(低压涡轮机匣进口安装边)→K(涡轮后机匣出口安装边)→N(内锥体尖部)→M(No.4轴承)”各点连线(对于安装边,为相连接两个零件之间的结构界面;对于其余部分,为零件与气体之间的气动界面)所形成的虚拟封闭区域(低压涡轮气动控制体)。
在发动机工作条件下(不同计算工况),在步骤一的中间过程数据中,沿着该控制体的轮廓线(构成虚拟封闭区域的各连线)选取结构参数(直径)和气动参数(包括静压、流量、流速),计算“单元体4”轴向力。
“单元体4”的气动轴向力基准值由以下3项组成:低压涡轮进口(包括流路、空气***腔、No.4轴承腔,轴向力计为FLPT_INT)、低压涡轮机匣-涡轮后机匣外壁面(机匣暴露于外涵流路部分的壁面,轴向力计为FLPT_CAS)、低压涡轮出口(包括涡轮后机匣出口流路、内锥体外壁面,轴向力计为FLPT_OUT),“单元体4”的气动轴向力为:
FITG_LPT=FLPT_INT+FLPT_CAS+FLPT_OUT (公式7)
5)“单元体5”为图4中“F(压气机机匣进口安装边)→H(外涵机匣进口安装边)→J(外涵机匣出口安装边)→K(低压涡轮机匣出口安装边)→L(高压涡轮机匣出口安装边)→F(压气机机匣进口安装边)”各点连线(对于安装边,为相连接两个零件之间的结构界面;对于其余部分,为零件与气体之间的气动界面)所形成的虚拟封闭区域(外涵道气动控制体)。
在发动机工作条件下(不同计算工况),在步骤一的中间过程数据中,沿着该控制体的轮廓线(构成虚拟封闭区域的各连线)选取结构参数(直径)和气动参数(包括静压、流量、流速),计算“单元体5”轴向力。
“单元体5”的气动轴向力基准值由以下4项组成:外涵流路进口(流路,轴向力计为FBYP_INT)、外涵机匣外壁面(机匣暴露于环境大气部分的壁面,轴向力计为FBYP_CAS)、外涵流路出口(流路,轴向力计为FBYP_OUT)、外涵流路内壁面(压气机机匣、主燃烧室机匣、高压涡轮机匣暴露于外涵流路部分的壁面,轴向力计为FBYP_COR),“单元体5”的气动轴向力为:
FITG_BYP=FBYP_INT+FBYP_CAS+FBYP_OUT+FBYP_COR (公式8)
6)“单元体6”为图4中“N(内锥体尖部)→K(涡轮后机匣出口安装边)→J(加力燃烧室进口安装边)→P(喷管出口外环)→O(发动机轴线)”各点连线(对于安装边,为相连接两个零件之间的结构界面;对于其余部分,为零件与气体之间的气动界面)所形成的虚拟封闭区域(出口区域气动控制体)。
在发动机工作条件下(不同计算工况),在步骤一的中间过程数据中,沿着该气动控制体的轮廓线(构成虚拟封闭区域的各连线)选取结构参数(直径)和气动参数(包括静压、流量、流速),计算“单元体6”轴向力。
“单元体6”的气动轴向力基准值由以下4项组成:加力部件内涵进口(包括流路、内锥体外壁面,轴向力计为FABE_COR_INT)、加力部件外涵进口(流路,轴向力计为FABE_BYP_INT)、加力-喷管外壁面(暴露于环境大气部分的壁面,轴向力计为FABE_CAS)、喷管出口流路(流路,轴向力计为FABE_OUT),计为:
FITG_ABE=FABE_COR_INT+FABE_BYP_INT+FABE_CAS+FABE_OUT (公式9)
上述过程中的第1~6项,即公式4~公式9的气动轴向力计算结果,均为单元体单独受到的来自于环境气体/工作气体的气动轴向力基准值,未叠加相邻单元体通过安装边传递的轴向力。
步骤四:将步骤三中单元体内部的各转、静子段轴向力叠加,得到相应单元体轴向力合成值。对比单元体轴向力合成值与步骤三中单元体轴向力基准值的差异,两者误差接近于零(两者误差与基准值的比值小于1%),则判断该过程人为失误消除;若存在明显误差,则返回步骤二、步骤三和步骤四,复查相应单元体轴向力基准值、合成值详细计算过程,排除人为失误。
如果存在人为失误的单元体内安装边较多,例如安装边数量多于3个,可参照步骤三中的方式,对单元体进行再次细化,直到每个气动控制体安装边的数量为2~3个,以便于更高效的定位并消除人为失误源。
如图5所示,以“单元体1”为例,说明单元体轴向力合成值的计算过程:
1)“零组件1”为图5中“D(No.1轴承)→A(进气机匣进口安装边)→a1(进气机匣后安装边)→a2(风扇第一级转子叶片进口叶尖)→D(No.1轴承)”各点连线(对于安装边,为相连接两个零件之间的结构界面;对于其余部分,为零件与气体之间的气动界面)所形成的虚拟封闭区域(气动控制体)。在经过步骤二校验和修正后的步骤一中间过程数据中,提取“零组件1”气动轴向力,计为F01_FAN;
2)“零组件2”为图5中“a1(风扇前机匣前安装边)→b1(风扇前机匣后安装边)→a6(风扇第二级转子叶片进口叶尖)→a5(风扇第二级转子叶片进口叶根)→a4(风扇第一级转子叶片出口叶根)→a3(风扇第一级转子叶片出口叶尖)→a2(风扇第一级转子叶片进口叶尖)→a1’(风扇前机匣前安装边)”各点连线(对于安装边,为相连接两个零件之间的结构界面;对于其余部分,为零件与气体之间的气动界面)所形成的虚拟封闭区域(气动控制体)。在经过步骤二校验和修正后的步骤一中间过程数据中,提取“零组件2”气动轴向力,计为F02_FAN;
3)“零组件3”为0中“b1(风扇后机匣前安装边)→B(风扇机匣出口安装边)→b5(风扇第三级转子叶片进口叶尖)→b4(风扇第三级转子叶片进口叶根)→b3(风扇第二级转子叶片出口叶根)→b2(风扇第二级转子叶片出口叶尖)→a6(风扇第二级转子叶片进口叶尖)→b1(风扇后机匣前安装边)”各点连线(对于安装边,为相连接两个零件之间的结构界面;对于其余部分,为零件与气体之间的气动界面)所形成的虚拟封闭区域(气动控制体)。在经过步骤二校验和修正后的步骤一中间过程数据中,提取“零组件3”气动轴向力,计为F03_FAN;
4)“零组件4”为0中“D(No.1轴承)→a2(风扇第一级转子叶片进口叶尖)→a3(风扇第一级转子叶片出口叶尖)→a4(风扇第一级转子叶片出口叶根)→a5(风扇第二级转子叶片进口叶根)→a6(风扇第二级转子叶片进口叶尖)→b2(风扇第二级转子叶片出口叶尖)→b3(风扇第二级转子叶片出口叶根)→b4(风扇第三级转子叶片进口叶根)→b5(风扇第三级转子叶片进口叶尖)→b6(风扇第三级转子叶片出口叶尖)→C(No.2轴承)→D(No.1轴承)”各点连线(对于安装边,为相连接两个零件之间的结构界面;对于其余部分,为零件与气体之间的气动界面)所形成的虚拟封闭区域(气动控制体)。在经过步骤二校验和修正后的步骤一中间过程数据中,提取“零组件4”气动轴向力,计为F04_FAN;
5)则“单元体1”气动轴向力合成值为,
FSUM_FAN=F01_FAN+F02_FAN+F03_FAN+F04_FAN (公式9)
步骤五:完成各单元体轴向力基准值、合成值校验后,按以下判断条件检查各单元体之间的界面参数:
1)两个单元体界面两侧的结构尺寸(环形作用面内、外环直径)一致;
2)两个单元体界面两侧的气动参数(静压、流量、流速)一致;
3)两个单元体界面两侧的轴向力方向相反。
如果满足上述3项判据,则判断相邻单元体界面参数一致;如果存在不满足的情况,则返回步骤二、①,修正相应位置结构或气动参数,重新按顺序开展步骤三、④、⑤防错分析。
以图4为例,说明单元体之间的“界面”的含义为:
1)对于单元体1和单元体2,两单元体之间的界面指的是单元体1的C→B连线与单元体2的C→B连线;
2)对于单元体2和单元体3,两单元体之间的界面指的是单元体2的E→F连线与单元体3的E→F连线;
3)对于单元体2和单元体5,两单元体之间的界面指的是单元体2的F→H连线与单元体5的F→H连线;
4)对于单元体3和单元体4,两单元体之间的界面指的是单元体3的M→L连线与单元体4的M→L连线;
5)对于单元体3和单元体5,两单元体之间的界面指的是单元体3的F→L连线与单元体5的F→L连线;
6)对于单元体4和单元体5,两单元体之间的界面指的是单元体4的L→K连线与单元体5的L→K连线;
7)对于单元体4和单元体6,两单元体之间的界面指的是单元体4的N→K连线与单元体6的N→K连线;
8)对于单元体5和单元体6,两单元体之间的界面指的是单元体5的K→J连线与单元体6的K→J连线。
三、整机级轴向力防错分析
步骤六:将整个发动机视为一个气动控制体,根据发动机进气端、排气端、外壁面气动和结构参数计算发动机推力基准值。
以图4的发动机为例,该发动机整体指的是“D(No.1轴承)→A(进气机匣进口安装边)→P(喷管出口外环)→O(发动机轴线)”各点连线(为发动机与环境大气之间的气动界面)所形成的虚拟封闭区域(整机气动控制体)。
在发动机工作条件下(不同计算工况),在步骤一等中间过程数据(经过步骤二、步骤三、步骤四和步骤五校验、修正后的数据)中,沿着发动机的轮廓线(构成虚拟封闭区域的各连线)选取结构参数(直径)和气动参数(包括静压、流量、流速),计算“发动机”推力基准值。
发动机推力基准值由以下3项组成:发动机进口(包括流路、帽罩,轴向力计为FENG_INT)、发动机外壁面(暴露于环境大气部分的壁面,轴向力计为FENG_CAS)、喷管出口流路(流路,轴向力计为FENG_OUT),该整机气动轴向力为:
FENG_ITG=FENG_INT+FENG_CAS+FENG_OUT (公式10)
步骤七:按照发动机传力路径,将步骤一中各零组件轴向力(经过步骤二、步骤三、步骤四和步骤五校验、修正后的数据)叠加至主安装节所在的机匣,并叠加该机匣自身的轴向力,给出发动机推力合成值。若发动机推力合成值与基准值的误差接近于零(在本申请中,该误差小于1%),则判断前述所有过程已消除人为失误;若发动机推力合成值与基准值仍然存在明显差值,则返回步骤四重新检查,定位并修正相应人为失误后,重新完成步骤五、步骤六、步骤七。
以图4为例,按照发动机传力路径,传递给中介机匣的气动轴向力(外力)包括:
1)作用于风扇转子和低压涡轮转子的气动轴向力(即“单元体1”、“单元体4”中转子件受到的轴向力),通过No.1轴承传递至进气机匣,连同进气机匣、风扇机匣受到的气动轴向力(即“单元体1”中静子件受到的轴向力),共同传递至“B(中介机匣进口安装边)”,轴向力计为FINC_FWD,如图6所示;
2)作用于外涵机匣、加力筒体、喷管的气动轴向力(即单元体5的“H→J”连线和单元体6的“J→P”连线内、外壁面的气动轴向力),传递至“H(中介机匣出口外环安装边)”,轴向力计为FINC_BWD_BPS,如0所示;
3)作用于压气机、主燃烧室、高压涡轮、低压涡轮各机匣的气动轴向力(即“单元体3”、“单元体4”中静子件所受到的气动轴向力),传递至“F(中介机匣出口内环安装边)”,轴向力计为FINC_BWD_EGC,如图6所示;
4)作用于压气机转子、高压涡轮转子的气动轴向力(即“单元体3”中转子件所受到的气动轴向力),传递至“E(No.3轴承)”,轴向力计为FINC_HPR,如图5所示。
中介机匣自身轴向力,即为步骤三中“单元体2”的轴向力“FINC”(经过步骤四和步骤五校验、修正后的数据),连同1)~4)项外力,也会传递给主安装节。因此,传递给主安装节的发动机推力合成值为,
FENG_SUM=FINC_FWD+FINC_BWD_BPS+FINC_BWD_EGC+FINC_HPR+FINC (公式11)
若根据公式10、公式11计算的结果(FENG_SUM-FENG_ITG)/FENG_ITG<1%,则判断人为失误已消除;否则,返回步骤四重新检查,定位并修正相应人为失误后,重新完成步骤五、步骤六、步骤七。
本申请的方法能够完全杜绝整机轴向力计算过程中的人为失误,并且大大缩短误差源定位时间周期,大幅提高工作效率。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。
Claims (6)
1.一种航空发动机气动轴向力防错分析方法,其特征在于,所述方法包括:
步骤一、确定零组件轴向力防错分析所需中间过程数据,根据发动机各零组件的气动轴向力计算方法进行零组件级气动轴向力计算;
步骤二:对发动机零组件中各转子-静子、转子-转子、静子-静子界面进行检查,确保界面两侧选用的结构和气动参数符合一致性要求;
步骤三、将发动机按照部件划分为若干个由一个或多个零组件构成的单元体,根据各单元体的结构及气动参数分别得到各单元体轴向力基准值;
步骤四、将各单元体内一个或多个零组件的轴向力进行叠加得到各单元体轴向力合成值,判断各单元体轴向力基准值与合成值的差值是否满足要求,若满足要求,则判断单元体轴向力计算中人为失误消除;若不满足要求,则返回步骤二、步骤三和步骤四,复查相应单元体轴向力基准值、合成值的计算过程,排除人为失误;
步骤五、完成各单元体轴向力基准值与单元体轴向力合成值校验后,确保各单元体之间的界面参数符合一致性要求:
步骤六:将发动机整体视为一个气动控制体,根据发动机进气端、排气端、外壁面气动和结构参数计算发动机推力基准值;
步骤七:按照发动机传力路径,将步骤一中各零组件气动轴向力叠加至主安装节所在的机匣,并叠加该机匣自身的轴向力,得到发动机推力合成值;
若发动机推力合成值与基准值的误差满足要求,则判断所有过程已消除人为失误;若发动机推力合成值与基准值的误差不满足要求,则返回步骤四重新检查,定位并修正相应人为失误后,重新完成步骤五、步骤六、步骤七。
2.如权利要求1所述的航空发动机气动轴向力防错分析方法,其特征在于,所述零组件中各转子-静子、转子-转子、静子-静子界面两侧选用的结构和气动参数一致性要求为:
对于转子-静子界面:
1)转子侧与静子侧选用的结构尺寸一致,均为转子侧结构尺寸;
2)转子侧与静子侧气动参数一致,均为转子侧气动参数;
3)转子侧与静子侧轴向力受力大小相等、方向相反;
对于转子-转子界面:
1)两侧结构尺寸一致,均为其中一侧的结构尺寸;
2)两侧气动参数一致,均为其中一侧的气动参数;
3)两侧轴向力受力大小相等、方向相反;
对于静子-静子界面:
1)两侧结构尺寸一致,均为其中一侧的结构尺寸;
2)两侧气动参数一致,均为其中一侧的气动参数;
3)两侧轴向力受力大小相等、方向相反。
3.如权利要求1所述的航空发动机气动轴向力防错分析方法,其特征在于,各单元体轴向力基准值与合成值的差值与基准值的比值小于1%,则满足要求;反之,则不满足。
4.如权利要求3所述的航空发动机气动轴向力防错分析方法,其特征在于,当单元体轴向力基准值与合成值的差值不满足要求,且单元体内多个零组件所形成的安装边数量多于3个时,对单元体进行再次细化,直到所述单元体内的安装边数量为2~3个。
5.如权利要求1所述的航空发动机气动轴向力防错分析方法,其特征在于,各单元体之间的界面参数一致性要求判据包括:
1)两个单元体界面两侧的结构尺寸一致;
2)两个单元体界面两侧的气动参数一致;
3)两个单元体界面两侧的轴向力方向相反;
如果满足上述3项判据,则判断相邻单元体界面参数一致;如果存在不满足的情况,则修正相应位置结构或气动参数,重新按顺序开展防错分析。
6.如权利要求1所述的航空发动机气动轴向力防错分析方法,其特征在于,所述发动机推力合成值与基准值的误差满足要求为:
发动机推力合成值与基准值的误差与基准值的比例小于1%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210089680.3A CN114510798B (zh) | 2022-01-25 | 2022-01-25 | 一种航空发动机气动轴向力防错分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210089680.3A CN114510798B (zh) | 2022-01-25 | 2022-01-25 | 一种航空发动机气动轴向力防错分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114510798A CN114510798A (zh) | 2022-05-17 |
CN114510798B true CN114510798B (zh) | 2023-07-21 |
Family
ID=81550327
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210089680.3A Active CN114510798B (zh) | 2022-01-25 | 2022-01-25 | 一种航空发动机气动轴向力防错分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114510798B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115982894B (zh) * | 2023-03-20 | 2023-06-30 | 中国航发四川燃气涡轮研究院 | 带螺纹主安装节安装***与推力销间隙设计方法 |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1474038A (zh) * | 2003-07-16 | 2004-02-11 | 沈阳黎明航空发动机(集团)有限责任 | 燃气轮机用可倾瓦弹簧片组件 |
EP1847457A2 (en) * | 2006-04-22 | 2007-10-24 | Rolls-Royce plc | Aeroengine mounting |
CN105784380A (zh) * | 2016-04-19 | 2016-07-20 | 上海交通大学 | 航空发动机压气机转子叶片***性故障检测方法 |
EP3404209A1 (en) * | 2017-05-18 | 2018-11-21 | United Technologies Corporation | Turbine cooling arrangement |
CN109209641A (zh) * | 2018-10-31 | 2019-01-15 | 中国科学院工程热物理研究所 | 一种发动机转子组件的连接结构 |
CN208749801U (zh) * | 2018-08-01 | 2019-04-16 | 北京福田康明斯发动机有限公司 | 曲轴止推片、曲轴止推片组件以及发动机 |
CN110222401A (zh) * | 2019-05-30 | 2019-09-10 | 复旦大学 | 航空发动机非线性模型建模方法 |
CN110543694A (zh) * | 2019-08-13 | 2019-12-06 | 中国航发贵阳发动机设计研究所 | 一种航空发动机辅助吊挂拉杆振动计算方法 |
CN111241609A (zh) * | 2020-01-20 | 2020-06-05 | 西北工业大学 | 一种航空发动机转静子装配叶尖间隙的预测方法 |
CN111946460A (zh) * | 2020-08-10 | 2020-11-17 | 中国人民解放军第五七一九工厂 | 用于航空发动机风扇组合件装配的导正装置及装配方法 |
CN112307630A (zh) * | 2020-11-02 | 2021-02-02 | 中国航发沈阳发动机研究所 | 一种航空发动机主轴承寿命分析载荷谱编制方法 |
CN212931909U (zh) * | 2020-09-23 | 2021-04-09 | 中国航发商用航空发动机有限责任公司 | 一种工艺进气道及航空发动机实验平台 |
CN113177276A (zh) * | 2021-04-27 | 2021-07-27 | 中国航发沈阳发动机研究所 | 一种航空发动机高压轴止推轴承的载荷确定方法 |
CN113343357A (zh) * | 2021-06-23 | 2021-09-03 | 中国航发沈阳发动机研究所 | 一种航空发动机高低压转子间径向间隙计算方法 |
CN113536192A (zh) * | 2021-06-23 | 2021-10-22 | 中国航发沈阳发动机研究所 | 一种基于加力燃烧室油气比计算加力总燃油流量的方法 |
-
2022
- 2022-01-25 CN CN202210089680.3A patent/CN114510798B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1474038A (zh) * | 2003-07-16 | 2004-02-11 | 沈阳黎明航空发动机(集团)有限责任 | 燃气轮机用可倾瓦弹簧片组件 |
EP1847457A2 (en) * | 2006-04-22 | 2007-10-24 | Rolls-Royce plc | Aeroengine mounting |
CN105784380A (zh) * | 2016-04-19 | 2016-07-20 | 上海交通大学 | 航空发动机压气机转子叶片***性故障检测方法 |
EP3404209A1 (en) * | 2017-05-18 | 2018-11-21 | United Technologies Corporation | Turbine cooling arrangement |
CN208749801U (zh) * | 2018-08-01 | 2019-04-16 | 北京福田康明斯发动机有限公司 | 曲轴止推片、曲轴止推片组件以及发动机 |
CN109209641A (zh) * | 2018-10-31 | 2019-01-15 | 中国科学院工程热物理研究所 | 一种发动机转子组件的连接结构 |
CN110222401A (zh) * | 2019-05-30 | 2019-09-10 | 复旦大学 | 航空发动机非线性模型建模方法 |
CN110543694A (zh) * | 2019-08-13 | 2019-12-06 | 中国航发贵阳发动机设计研究所 | 一种航空发动机辅助吊挂拉杆振动计算方法 |
CN111241609A (zh) * | 2020-01-20 | 2020-06-05 | 西北工业大学 | 一种航空发动机转静子装配叶尖间隙的预测方法 |
CN111946460A (zh) * | 2020-08-10 | 2020-11-17 | 中国人民解放军第五七一九工厂 | 用于航空发动机风扇组合件装配的导正装置及装配方法 |
CN212931909U (zh) * | 2020-09-23 | 2021-04-09 | 中国航发商用航空发动机有限责任公司 | 一种工艺进气道及航空发动机实验平台 |
CN112307630A (zh) * | 2020-11-02 | 2021-02-02 | 中国航发沈阳发动机研究所 | 一种航空发动机主轴承寿命分析载荷谱编制方法 |
CN113177276A (zh) * | 2021-04-27 | 2021-07-27 | 中国航发沈阳发动机研究所 | 一种航空发动机高压轴止推轴承的载荷确定方法 |
CN113343357A (zh) * | 2021-06-23 | 2021-09-03 | 中国航发沈阳发动机研究所 | 一种航空发动机高低压转子间径向间隙计算方法 |
CN113536192A (zh) * | 2021-06-23 | 2021-10-22 | 中国航发沈阳发动机研究所 | 一种基于加力燃烧室油气比计算加力总燃油流量的方法 |
Non-Patent Citations (1)
Title |
---|
基于航空发动机脉动装配的智能管控技术研究;魏小红 等;《智能生产线》;第63卷(第6期);43-50 * |
Also Published As
Publication number | Publication date |
---|---|
CN114510798A (zh) | 2022-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20180306041A1 (en) | Multiple turbine vane frame | |
CN114510798B (zh) | 一种航空发动机气动轴向力防错分析方法 | |
CN112594064B (zh) | 一种基于轴流压气机级间测量参数的s2流场诊断方法 | |
US7941300B1 (en) | Process for the design of an airfoil | |
Wellborn et al. | Redesign of a 12-stage axial-flow compressor using multistage CFD | |
Ju et al. | Design optimization and experimental study of tandem impeller for centrifugal compressor | |
CN110321586B (zh) | 一种航空发动机偏离设计点工作状态迭代求解的取值方法 | |
US20040030666A1 (en) | Method of designing a multi-stage compressor rotor | |
US20240159154A1 (en) | Turbine blade tip shroud with axially offset cutter teeth, and related surface profiles and method | |
EP2806104B1 (en) | Aerofoil recambering | |
EP3819805A1 (en) | Parametric component design process | |
CN113468788B (zh) | 一种发动机整体叶盘变形协调设计方法 | |
US11480073B2 (en) | Gas turbine engine nacelle and method of designing same | |
CN114065572B (zh) | 一种厚度偏差与叶片厚度分布相关的叶片模型构建方法 | |
Salnikov et al. | Multidisciplinary design optimization of a bladed disc for small-size gas-turbine engines | |
US11105206B1 (en) | Turbine airfoil | |
US20210264072A1 (en) | System and process for designing internal components for a gas turbine engine | |
Immery et al. | Design of the Compression System of a Geared Turbofan | |
Lück et al. | Accuracy assessment of steady and unsteady multistage high pressure compressor simulations | |
Li et al. | The optimization of a centrifugal impeller based on a new multi-objective evolutionary strategy | |
CN111486126A (zh) | 一种涡轮增压器叶轮及其五轴加工方法 | |
US11566525B1 (en) | Turbine blade airfoil profile | |
Muchowski et al. | Influence of axial compressor model simplification and mesh density on surge margin evaluation | |
Becker et al. | Multi-objective optimization in axial compressor design using a linked cfd-solver | |
US11098591B1 (en) | Turbine blade with contoured fillet |
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 |