轴流压气机上限流量测试方法技术

一种轴流压气机上限流量测试方法，首先根据压气机的几何参数和边界条件经前处理构造出适于流线曲率法计算的流道，而后设定上限流量的初始迭代值及其递增量、迭代次数，每个迭代步中进行流场计算和上限流量判别，直至判断出上限流量，终止时的迭代值作为上限流量，最后记录上限流量并对其评估，本发明专利技术预测跨声速工况上限流量精度高，不依赖压气机经验数据库，节省了设计方案所需的成本，计算流程模块化，与流场计算结果耦合，可靠性高。

【技术实现步骤摘要】

本专利技术涉及的是一种航空发动机领域的技术，具体是一种轴流压气机上限流量测试方法。
技术介绍
轴流压气机采用基于流线曲率法的设计方案时，流量作为初始推导参数，由于不存在上限流量的理论解，需要通过计算经验模型的方式得到方案设计中的上限流量，即确定压气机的最大工作流量。现有方法是根据设计方案模块化以往压气机实验数据，依照经验库提供的结果求解出整机设计状态下的流量，再由设计状态流量来预测上限流量。但该方法受经验库内压气机设计水平的限制，难以预测出新型压气机设计方案中的流动特点。现有的轴流压气机转子叶尖相对马赫数可达到1.5，轴流压气机的超声速区域流动范围扩大。超声速流动越显著，压气机的上限流量受到超声速流动理论限制愈专利技术显，使得通过经验库预测的精度偏低。
技术实现思路
本专利技术针对现有技术无法得到非设计状态的上限流量等不足，提出一种轴流压气机上限流量测试方法。本专利技术是通过以下技术方案实现的：本专利技术首先根据压气机的几何参数和边界条件经前处理构造出适于流线曲率法计算的流道，而后设定上限流量的初始迭代值及其递增量、迭代次数，每个迭代步中进行流场计算和上限流量判别，直至判断出上限流量为止，终止时的迭代值作为上限流量，最后记录上限流量并对其评估。本专利技术包括以下具体步骤：1)对设计方案确定的压气机几何参数和边界条件进行前处理，构造出适于流线曲率法计算的流道；2)设定上限流量的初始迭代值、迭代次数和递增量；3)进行本迭代步的流场计算；4)对流场计算结果进行上限流量判别，判断出上限流量则终止，否则，若未达到迭代次数，按递增量得到新的迭代值并回到步骤3)，若达到迭代次数，重新设定迭代次数并回到步骤2)；5)将终止时的迭代值作为上限流量并进行记录评估。所述的步骤1)具体包括以下步骤：1.1)依次构造外机匣、轮毂、叶片前缘和叶片尾缘的型线函数；1.2)构造无叶片区内计算站的型线函数；1.3)依照流量径向分配需求确定流线分布规律，并应用于所有计算站，从而构造出流道的流线型线函数。所述的流线分布规律为：其中：i＝1,2,...,N-1，S(i)为第i条流线距轮毂的距离，N为流线总数，A、B、C和D选取不同值则流量径向分配需求不同，满足A+B＝C+D＝1，且A、B、C和D皆为正值。所述的步骤3)具体包括以下步骤：3.1)确定流场计算次数；3.2)计算初始流场；3.3)逐一求解流场，求解各个计算站的流线曲率法控制方程组；3.4)若流线曲率法控制方程组收敛或达到流场计算次数，则进行下一步，否则重新构造流线型线函数，并回到步骤3.2)。3.5)储存并输出流场计算结果。所述的流线曲率法控制方程组包括：熵方程、连续性方程和动量方程。所述的连续性方程为：其中：ρ为密度，W为相对速度，r为径向，z为轴向，n为流面法向量。所述的动量方程为：其中：p为压力，E为粘性应力，ω为角速度。所述的熵方程为：其中：s为熵，cp为等压比热容，为相对滞止温度，R为气体常数，为相对滞止压力。所述的步骤4)具体包括以下步骤：4.1)如果流场计算结果发散时，到达迭代次数则重新设置更大的迭代次数，重新执行步骤2)，未达迭代次数则按递增量得到新的迭代值并回到步骤3)，如果流场计算结果收敛则执行下一步；4.2)调用流场计算结果中的部分参数；4.3)应用气动拥塞公式和上限流量判断准则判断上限流量，若迭代值为上限流量则终止整个迭代，若迭代值不为上限流量且未达迭代次数则按递增量得到新的迭代值并回到步骤3)，若迭代值不为上限流量且达到迭代次数则重新设置更大的迭代次数，重新执行步骤2)。所述的气动拥塞公式为：其中：临界声速临界密度流线所在截面的无量纲临界面积为：所述的上限流量判断准则指达到气动拥塞的流线流量之和大于当前迭代值的一半时的迭代值即为上限流量。技术效果与现有技术相比，本专利技术预测跨声速工况上限流量精度高，不依赖压气机经验数据库，节省了设计方案所需的成本，计算流程模块化，与流场计算结果耦合，可靠性高。附图说明图1为本专利技术流程示意图；图2为三维几何与子午面投影示意图；图3为流道示意图；图4为流场计算结果图；图5为上限流量判断结果与试验结果对比图。具体实施方式下面对本专利技术的实施例作详细说明，本实施例在以本专利技术技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本专利技术的保护范围不限于下述的实施例。实施例1本实施例通过轴流压气机上限流量测试方法来预测一台跨声速轴流压气机的上限流量，即首先根据压气机的几何参数和边界条件经前处理构造出适于流线曲率法计算的流道，而后设定上限流量的初始迭代值及其递增量、迭代次数，每个迭代步中进行流场计算和上限流量判别，直至判断出上限流量，终止时的迭代值作为上限流量，最后记录上限流量并对其评估。所述的压气机的几何参数包括：前缘半径、前缘轴向位置、尾缘半径、尾缘轴向位置、轮毂半径、轮毂轴向位置、外机匣半径、外机匣轴向位置、叶顶间隙、前缘几何角、尾缘几何角、前缘小圆厚度、叶片最大厚度、尾缘小圆厚度、叶片最大厚度处轴向位置、罩量角、弦长、安装角和叶片数。所需的边界条件包括：进口滞止温度、滞止压力和周向速度。压气机的主要设计参数如表1所示。表1压气机设计参数如图1所示，压气机上限流量的预测包括以下步骤：1)对设计方案确定的如图2所示的压气机几何参数和边界条件进行前处理，构造出适于流线曲率法计算的流道：1.1)依次构造外机匣、轮毂、叶片前缘和叶片尾缘的型线函数。对输入的几何参数进行三次自然样条插值获得几何参数的型线函数即几何型线，把已知的离散几何点转化为连续函数。1.2)构造无叶片区内计算站的型线函数。如图3所示，按照相同的插值方法，在进气段和出气段的无叶片区构造了径向计算站，且进气段的计算站之间的轴向距离按照一定比例衰减。构造了叶片之间无叶片区的计算站，使计算站的轴向位置位于两排叶片的中点处，总共构造了18处计算站。1.3)依照流量径向分配需求确定流线分布规律，并应用于无叶片区内计算站、叶片前缘和叶片尾缘，从而构造出流线型线函数，如图3所示构造了36条流线，形成流道。所述的流线分布规律为：其中：i＝1,2,...,N-1，S(i)为第i条流线距轮毂的距离，N为流线总数，A、B、C和D选取不同值则流量径向分配需求不同。但需满足A+B＝C+D＝1，且A、B、C和D皆取正值。2)设定上限流量的初始迭代值，迭代次数和递增量。初始迭代值可为小于上限流量的任意正值，递增量按照需求的测试精度决定，迭代次数的设置仅为终止异常值引起的无限循环，故应使得达到迭代次数时的迭代值大于上限流量，若因迭代次数过少引起计算终止，仅需将迭代次数调大即可，选取初始迭代值为20kg/s，递增量为0.01kg/s，迭代次数为1000次。3)进行本迭代步的流场计算：3.1)确定流场计算次数，流场计算次数过大会显著延长计算时间，流场计算次数过少则流场难以收敛，这里选择流场计算次数为1000次。3.2)计算初始流场，气体动力学管道流动规律即忽略叶片对流动的影响，起到构造初始流场的作用。在初始流量即初始迭代值为20kg/s，管道进口流速为165m/s，由于流道不断收缩，在亚声速条件下流动持续加速，获得初始流场。3.3)逐一求解各个计算站的流线曲率法控制方程组，各个计算站首先迭代求解动量方程，本文档来自技高网...

【技术保护点】
一种轴流压气机上限流量测试方法，其特征在于，首先根据压气机的几何参数和边界条件经前处理构造出适于流线曲率法计算的流道，而后设定上限流量的初始迭代值及其递增量、迭代次数，每个迭代步中进行流场计算和上限流量判别，直至判断出上限流量，终止时的迭代值作为上限流量，最后记录上限流量并对其评估。

【技术特征摘要】
1.一种轴流压气机上限流量测试方法，其特征在于，首先根据压气机的几何参数和边界条件经前处理构造出适于流线曲率法计算的流道，而后设定上限流量的初始迭代值及其递增量、迭代次数，每个迭代步中进行流场计算和上限流量判别，直至判断出上限流量，终止时的迭代值作为上限流量，最后记录上限流量并对其评估。2.根据权利要求1所述的轴流压气机上限流量测试方法，其特征是，具体包括以下步骤：1)对设计方案确定的压气机几何参数和边界条件进行前处理，构造出适于流线曲率法计算的流道；2)设定上限流量的初始迭代值、迭代次数和递增量；3)进行本迭代步的流场计算；4)对流场计算结果进行判别，判断出上限流量则终止，否则，若未达到迭代次数，按递增量得到新的迭代值并回到步骤3)，若达到迭代次数，重新设定迭代次数并回到步骤2)；5)将终止时的迭代值作为上限流量并进行记录评估。3.根据权利要求2所述的轴流压气机上限流量测试方法，其特征是，所述的步骤1)具体包括以下步骤：1.1)依次构造外机匣、轮毂、叶片前缘和叶片尾缘的型线函数；1.2)构造无叶片区内计算站的型线函数；1.3)依照流量径向分配需求确定流线分布规律，并应用于无叶片区内计算站、叶片前缘和叶片尾缘，从而构造出流道的流线型线函数。4.根据权利要求3所述的轴流压气机上限流量测试方法，其特征是，所述的流线分布规律为：其中：S(i)为第i条流线距轮毂的距离，N为流线总数，A、B、C和D选取不同值则流量径向分配需求不同，A+B＝C+D＝1，且A、B、C和D皆为正值。5.根据权利要求3所述的轴流压气机上限流量测试方法，其特征是，所述的步骤3)具体包括以下步骤：3.1)确定流场计算次数；3.2)计算初始流场；3....

【专利技术属性】
技术研发人员：吴东润，斯夏依，滕金芳，羌晓青，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：上海;31