一种发动机部件特性拓展方法技术

本发明专利技术公开了一种发动机部件特性拓展方法，向已有部件特性数据中添加参数化的零转速线特性数据，基于径向基函数近似建模技术建立部件特性近似模型，根据近似模型得到待拓展特性并计算特性拓展误差，再利用差分进化算法对零转速线参数进行优化，最终得到符合误差限制需求的特性结果。应用本发明专利技术技术方案的发动机部件特性拓展方法，能够获取发动机部件拓展特性，同时保证所获取的特性光顺、精确且不产生有违实际物理意义的特性点。解决了现有常用发动机部件特性拓展方法，需要人为额外设置参数，并且特性拓展结果曲线奇异、不精确，有违实际物理意义的问题。

A Method for Expanding the Characteristics of Engine Components

【技术实现步骤摘要】
一种发动机部件特性拓展方法
本专利技术涉及航空发动机领域，具体而言，涉及一种发动机部件特性拓展方法。
技术介绍
对于航空发动机而言，对部件特性进行拓展对研究航空发动机在风车、起动工况下的整机特性具有重要意义。无论是数值仿真或部件试验都只能获取有限转速范围的精确部件特性，因此在工程实践中往往通过现有的精确部件特性图，通过相似原理或特性变化规律对部件特性图进行拓展。然而对于现有的部件特性拓展技术，往往需要工程人员根据自身经验给定一些拓展系数，且往往对部件特性的拓展仅仅依靠转速较低的几条特性线，无法充分利用已有的全部特性数据，并且拓展得到的特性结果有可能会出现曲线奇异的现象，往往需要人为对特性进行修正、光顺。不仅如此，许多方法由于自身限制导致拓展的特性结果会受限于某些特性区域。这些缺陷会导致特性拓展方法不通用，拓展结果不精确，受人为因素影响大等问题。进而会导致对发动机在风车、起动工况下整机性能参数数值模拟精度差甚至基本工作特性错误。因此，发展新的部件特性拓展方法以获取精确、光顺、不受人为因素影响的部件特性具有至关重要的意义。
技术实现思路
本专利技术解决的技术问题是：本专利技术的主要目的在于提供一种发动机部件特性拓展方法，以至少解决现有技术中特性拓展结果奇异、范围受限、精度低的问题。本专利技术的技术方案是：一种发动机部件特性拓展方法，包括以下步骤：步骤一，发动机部件已有特性数据为压比、换算流程和效率，将已有特性数据中的效率转换为换算能量参数；得到新的特性数据参数，即零转速线参考点，表示为(Wc_p,Pr_p,Tc_p)；Pr_p为参考点压比，Wc_p为参考点换算流量，Tc_p为参考点换算能量参数；其中转换方式如下式所示：Tc＝(πk-1/k-1)/ηTc为换算能量参数，k为绝热指数，η为效率；步骤二，建立进化种群，种群大小为m；即进化种群中含有m个父代个体，且每个父代个体均包含(Wc_p,Pr_p,Tc_p)三个特性参数；初始化m个父代个体，初始化的参数取值需要保证Wc_p＞0,0＜Pr_p＜1,Tc_p＞0。步骤三，对步骤二中的m个父代个体采取进化操作；根据当前的零转速线参考点进化种群中的m个父代个体，采用差分进化算法的变异、交叉算子，得到m个子代个体。步骤四，初始化进化种群适应度计算进程。取i＝1，其中i代表当前计算适应度的个体序号，此时进化种群中包含父代个体及步骤三中产生的子代个体，因此种群大小为2*m；步骤五，向已有特性数据中加入以进化种群中第i个个体的参考点数值参数化表达的零转速线特性数据，包括以下子步骤：子步骤一，给定参数化表达的零转速线。根据进化种群中第i个个体的参考点(Wc_p,Pr_p,Tc_p)的三个参数，得到当前零转速线表达式。零转速线可参数化表达为：其中，Pr为压比，Wc为换算流量，Tc为换算能量参数，Kp和KT为待定系数。待定系数可用下列公式表达：子步骤二，将参数化的零转速线特性数据加入已有特性数据中。根据已有特性数据中各个转速线上的特性点数，利用参数化的零转速线计算等量的特性数据点，并将这些特性数据点加入已有特性数据中。步骤六，根据步骤五得到的特性数据建立近似模型，并计算近似模型总误差，包括以下子步骤：子步骤一，初始化近似模型总误差计算过程。设置n为2，近似模型总误差为0。子步骤二，获得当前近似模型建立数据。在步骤五中获得的特性数据中，从低至高选取第n条转速线特性数据并将其去除；子步骤三，初始化径向基函数近似模型。建立O·M个径向基函数，其中O为包括零转速特性线、待拓展转速特性线以及已有特性线的特性线总条数，M则为一条转速特性线上放置的径向基函数个数，O·M需要小于用于近似建模的特性数据点数。径向基函数的中心位置向量由转速和Beta值确定。其中转速即为该径向基函数被分配到的转速特性线的转速，而同一条转速特性线上的径向基函数中心位置向量中的Beta值则按照取值范围为0至1的等差数列赋值。子步骤四，计算当前近似模型误差。根据子步骤二中得到的特性数据，以及子步骤三中初始化的径向基函数近似模型，建立径向基函数近似模型。并根据近似模型计算子步骤一中被去除的转速线特性数据，并与被去除的特性数据进行对比并计算均方根误差。子步骤五，计算近似模型总误差。子步骤一中的n值加1，将子步骤四中得到的均方根误差累加至子步骤一中的近似模型总误差。若n值小于已有特性线数据总条数，则转至子子步骤二，否则转至步骤七。步骤七，计算进化种群中第i个个体的适应度。设置步骤五中选取的第i个个体的适应度为步骤六中计算得到的近似模型总误差。步骤八，判断是否所有个体适应度计算完毕。i值加一，若i值大于2*m，则转至步骤九，否则转至步骤五。步骤九，对进化种群进行筛选。根据步骤七中计算得到的适应度，对进化种群中m个父代个体和m个子代个体进行从小至大排序，选取前m个个体替换步骤二中的父代个体，并舍弃其余个体。步骤十，判断优化是否结束。若优化进程达到结束条件，则转至步骤十，否则转至步骤三。优化结束条件一般可选取为总优化步数。步骤十一，计算拓展结果，包括以下子步骤：子步骤一，得到优化后的零转速线参考点参数。选定步骤九中排序第一的个体，并得到其参考点参数(Wc_p,Pr_p,Tc_p)；子步骤二，得到近似建模所用的数据。向已有特性数据中加入以子步骤一中得到的参考点参数(Wc_p,Pr_p,Tr_p)参数化表达的零转速线特性数据。子步骤三，建立径向基函数近似模型，根据子步骤二中得到的特性数据建立径向基函数近似模型。子步骤四，得到待拓展结果。根据子步骤三中得到的径向基函数近似模型计算待拓展转速线的特性结果，此结果即是最终拓展结果。专利技术效果本专利技术的技术效果在于：应用本专利技术技术方案的发动机部件特性拓展方法，通过零转速特性线的引入，以向已有特性数据中加入隐含的低转速区域特性物理规律，并利用径向基函数近似建模技术挖掘这一规律以形成能够较好地拟合部件特性的近似模型。从而将基于数学拟合的特性拓展方法和基于物理规律的特性拓展方法有机结合。从而使利用本专利技术技术方案的发动机部件特性拓展方法可以综合部件物理规律和全部已有特性的数学规律，克服利用传统部件特性拓展方法时参数选取没有理论指导、拓展结果受人为因素影响大且不能正确反映部件在低转速下的物理性质的缺陷。特别地，对于压缩部件低转速特性拓展而言，应用本方法可以获取压缩部件在低转速下的“压气机模式”、“搅拌器模式”、“涡轮模式”下的部件特性，同时拓展特性结果光顺且无物理意义不明确的点。且近似模型的精度，即拓展特性的误差估计可降低至流量均方根误差在1.0*10-3以下，压比和换算能量参数的均方根误差在1.5*10-3以下。附图说明图1是根据本专利技术实施例可选的一种发动机特性拓展方法流程示意图；图2是根据本专利技术实施例可选的一种发动机特性拓展方法中近似模型总误差计算方法；图3是根据本专利技术实施例可选的一种发动机特性拓展方法中近似建模径向基函数中心参数选取方法；其中，在上述附图中包括以下参数；n：当前转速线序号；N：转速线总条数；m：每条转速线上径向基函数序号；M：每条转速线上径向基函数总数；O：包含零转速线、待拓展转速线以及已有转速线的条数；BetaRbf(n,m)：第n条转速线上第m个径向基函数的beta值即该径向基函数中心参数本文档来自技高网...

【技术保护点】
1.一种发动机部件特性拓展方法，其特征在于，包括以下步骤：步骤一，发动机部件已有特性数据为压比、换算流程和效率，将已有特性数据中的效率转换为换算能量参数；得到新的特性数据参数，即零转速线参考点，表示为(Wc_p,Pr_p,Tc_p)；Pr_p为参考点压比，Wc_p为参考点换算流量，Tc_p为参考点换算能量参数；其中转换方式如下式所示：Tc＝(π

【技术特征摘要】
1.一种发动机部件特性拓展方法，其特征在于，包括以下步骤：步骤一，发动机部件已有特性数据为压比、换算流程和效率，将已有特性数据中的效率转换为换算能量参数；得到新的特性数据参数，即零转速线参考点，表示为(Wc_p,Pr_p,Tc_p)；Pr_p为参考点压比，Wc_p为参考点换算流量，Tc_p为参考点换算能量参数；其中转换方式如下式所示：Tc＝(πk-1/k-1)/ηTc为换算能量参数，k为绝热指数，η为效率；步骤二，建立进化种群，种群大小为m；即进化种群中含有m个父代个体，且每个父代个体均包含(Wc_p,Pr_p,Tc_p)三个特性参数；初始化m个父代个体，初始化的参数取值需要保证Wc_p＞0,0＜Pr_p＜1,Tc_p＜0。步骤三，对步骤二中的m个父代个体采取进化操作；根据当前的零转速线参考点进化种群中的m个父代个体，采用差分进化算法的变异、交叉算子，得到m个子代个体。步骤四，初始化进化种群适应度计算进程。取i＝1，其中i代表当前计算适应度的个体序号，此时进化种群中包含父代个体及步骤三中产生的子代个体，因此种群大小为2*m；步骤五，向已有特性数据中加入以进化种群中第i个个体的参考点数值参数化表达的零转速线特性数据，包括以下子步骤：子步骤一，给定参数化表达的零转速线。根据进化种群中第i个个体的参考点(Wc_p,Pr_p,Tc_p)的三个参数，得到当前零转速线表达式。零转速线可参数化表达为：其中，Pr为压比，Wc为换算流量，Tc为换算能量参数，Kp和KT为待定系数。待定系数可用下列公式表达：子步骤二，将参数化的零转速线特性数据加入已有特性数据中。根据已有特性数据中各个转速线上的特性点数，利用参数化的零转速线计算等量的特性数据点，并将这些特性数据点加入已有特性数据中。步骤六，根据步骤五得到的特性数据建立近似模型，并计算近似模型总误差，包括以下子步骤：子步骤一，初始化近似模型总误差计算过程。设置n为2，近似模型总误差为0。子步骤二，获得当前近似模型建立数据。在步骤五中获得的特性数据中，从低至高选取第n条转速线特性数据并将其去除；子步骤三，初始化径向基函数近似模型。...

【专利技术属性】
技术研发人员：叶一帆，王占学，张晓博，张明阳，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西,61