一种基于空间插值的多因素敏感性分析方法技术

本发明专利技术涉及一种基于空间插值的多因素敏感性分析方法，属于空气动力技术领域。首先选定目标参量，筛选出影响因素。然后进行确定性方案试验。根据试验方案确定影响因素变化范围，选择空间插值方法补全试验数据。最后，基于试验数据，利用方差分析进行多因素敏感性分析。本方法与现有技术相比，能够以较少的实验次数刻画出敏感性分析所需数据，同时填充气动敏感性分析中存在的失效点，不仅节省了计算与实验资源，也能够快速、准确的实现气动敏感性分析，进一步提升敏感性分析效率。此外，针对试验方案中存在失效点的问题，本发明专利技术能够对实验方案中的失效点进行填充。方案中的失效点进行填充。方案中的失效点进行填充。

【技术实现步骤摘要】
一种基于空间插值的多因素敏感性分析方法


[0001]本专利技术涉及一种基于空间插值的多因素敏感性分析方法，属于空气动力


技术介绍

[0002]得益于空气动力学的发展，包括飞行器、压气机等在内的外流、内流领域都处于快速发展阶段。在气动研究中，影响气动性能与计算方法可靠性的因素众多。在这些因素中，不仅包括像攻角、速度、压力等常规影响气动性能的因素，也包括网格量、湍流模型等影响计算方法的因素。众多影响因素对于气动参量的影响程度不同，为了掌握这些影响因素的影响情况，进而为气动优化、气动设计提供指导，迫切需要进行多因素敏感性分析。
[0003]敏感性分析，是指从定量分析的角度研究影响因素发生变化对于研究目标影响程度的一种技术手段。根据作用范围的不同，敏感性分析通常可分为局部与全局敏感性分析两种。其中，局部法适用于局部参数对输出响应的影响，它的使用范围受到参数的变动范围、模型空间等影响。全局敏感性分析是综合整体考虑，分析参数在整个参数空间的对输出的变化影响。在实际工程应用中，相比局部敏感性分析，全局敏感性分析能够更好地从多方位考虑，操作简便，因此，全局敏感性分析受用范围广泛而有效。
[0004]由于全局敏感性分析需要大量试验方案来研究不同影响因素的敏感性，随着影响因素的增加，实验成本也在快速增加。尽管可以通过相应的实验设计方法减小实验次数，减小实验成本，但是其气动实验的成本仍十分巨大，气动仿真实验的耗时也极其漫长。此外，由于空气动力学的复杂性，在实验设计的方案中不可避免的存在一些实验无法测到的方案抑或是仿真无法实现的方案。因此，对于气动敏感性分析而言，急需一种既能够减少实验成本，也能够对失效点进行填充的新的敏感性分析方法。

技术实现思路

[0005]本专利技术的目的是为了解决在气动参量影响因素敏感性分析过程中面临的实验次数过多、实验资源巨大，同时存在失效点等问题，提出一种基于空间插值的多因素敏感性分析方法。
[0006]一种基于空间插值的多因素敏感性分析方法，包括以下步骤：
[0007]首先，选定目标参量，筛选影响因素。
[0008]对于空气动力学中的内外流问题，首先要明确研究对象与研究类型。其中，研究对象是指包括飞行器、压气机等在内的具体模型，研究类型包括实验与数值模拟。
[0009]目标参量，是指研究对象的内外性能参数，对于外流目标参量，包含升力、阻力、压力、速度等外部的性能参数。对于内流目标参量，包含流量、效率、压比、损失系数等。
[0010]以实验为研究类型的研究影响目标参量因素包括：来流马赫数、攻角、附面层特性、表面特性、气动噪声、加工不确定、燃烧、间隙、积冰、非定常流动，等等。以数值模拟方法进行研究的影响因素包括：粗糙度、转捩模型、网格独立性、传热模型、湍流模型、结构/非结
构网格、雷诺平均、大涡模拟，等等。
[0011]影响目标参量的因素有很多，但受到试验条件限制，不可能全面进行考察，因此，技术人员需要通过对实际具体问题的具体分析并结合实验目的，选出主要因素，略去次要因素。同时，要根据具体空气动力学问题确定影响因素的变化范围。
[0012]之后，进行确定性方案试验。
[0013]在明确影响目标参量的因素后，要进行试验设计，目的是能够合理的安排试验因素，控制实验成本，并且能够有效分析实验数据。试验设计方法包括：正交实验设计、回归正交试验设计、均匀设计，等等。
[0014]确定性方案是指：全部试验减去失效方案后的剩余方案，以及受客观条件限制进行的有限个方案。
[0015]其中，失效方案是指：由于气动参量与影响因素相互作用的复杂性，在试验设计的诸多方案中，必然存在实验条件无法测到的方案，或者数值模拟无法计算的方案。实验条件无法测到的方案是指流动由于激波、非定常效应的存在导致无法测量或数据失效的方案，数值模拟无法计算的方案是指计算无法收敛等无法计算的方案。
[0016]受客观条件限制进行的有限方案是指：客观条件包括实验条件、计算资源、时间、人力，等等。由于实验条件有限，对于试验设计的全部方案无法一一进行试验，在有限的实验成本内只能完成数量有限的试验方案。
[0017]然后，选择空间插值方法并补全试验数据。
[0018]空间插值，是指将离散的测量数据转换为连续的数据曲面。空间插值法包括克里金插值法、反距离加权法。空间插值的目的是根据有限个点的数值，通过距离或其它空间关系，推测出其它位置。空间插值方法的选择，具体要根据实际问题进行具体分析选择。利用空间插值方法补全试验数据的步骤如下：
[0019]根据试验方案确定影响因素变化范围，选择空间插值方法，补全试验数据。
[0020]最后，基于试验数据，利用方差分析，进行多因素敏感性分析。
[0021]敏感性分析能够在众多影响因素中指出对于目标有重要影响的因素，并评估各个因素对于目标的影响程度。
[0022]在实践中，方差分析是一种常用的敏感性分析方法。通常，目标函数的总方差分解为各因素引起的方差以及因素之间耦合作用引起的方差，而敏感度的大小则通过各个因素引起的方差占总方差的比例大小来确定：
[0023][0024]其中，表示总体样本的均值；x
i
表示设计变量；是设计变量x
i
的均值；y(x1,
…
,x
n
)表示目标参量；σ2表示总体样本的方差；表示设计变量x
i
的方差。
[0025]有益效果
[0026]本专利技术方法，与现有技术相比，能够以较少的实验次数刻画出敏感性分析所需数据，同时填充气动敏感性分析中存在的失效点，不仅节省了计算与实验资源，也能够快速、
准确的实现气动敏感性分析，进一步提升敏感性分析效率。此外，针对试验方案中存在失效点的问题，本专利技术能够对实验方案中的失效点进行填充。
附图说明
[0027]图1为本专利技术的流程图；
[0028]图2为对照组的流量敏感性分析；
[0029]图3为对照组的压比敏感性分析；
[0030]图4为对照组的效率敏感性分析；
[0031]图5为对照组的气动参数敏感性；
[0032]图6为本专利技术的流量敏感性分析；
[0033]图7为本专利技术的压比敏感性分析；
[0034]图8为本专利技术的效率敏感性分析；
[0035]图9为本专利技术的气动参数敏感性分析。
具体实施方式
[0036]为了更好地说明本专利技术的目的和优点，下面结合附图对本专利技术方法的具体实施方式作详细说明。
[0037]实施例
[0038]本实例针对一高负荷离心压气机气动敏感性进行分析，实验方法为数值模拟。离心压气机主要由叶轮、径向扩压器、轴向扩压器部件组成。其中，决定离心压气机整机性能的气动参数包括流量、压比、效率。
[0039]如图1所述，一种基于空间插值的多因素敏感性分析方法，包括以下步骤：
[0040]步骤1本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于空间插值的多因素敏感性分析方法，其特征在于，包括以下步骤：步骤1：选定目标参量，筛选影响因素；目标参量，是指研究对象的内外性能参数，其中，外流目标参量，至少包括升力、阻力、压力和速度，内流目标参量至少包括流量、效率、压比和损失系数；以实验为研究类型的研究影响目标参量因素包括：来流马赫数、攻角、附面层特性、表面特性、气动噪声、加工不确定、燃烧、间隙、积冰、非定常流动；以数值模拟方法进行研究的影响因素包括粗糙度、转捩模型、网格独立性、传热模型、湍流模型、结构/非结构网格、雷诺平均大涡模拟；步骤2：进行确定性方案试验，完成数值模拟计算；方案试验设计方法包括：正交实验设计、回归正交试验设计、均匀设计；其中，确定性方案是指：全部试验减去失效方案后的剩余方案，以及受客观条件限制进行的有限个方案；失效方案是指：由于气动参量与影响因素相互作用的复杂性，在试验设计的诸多方案中，必然存在实验条件无法测到的方案，或者数值模拟无法计算的方案；实验条件无法测到的方案是指流动由于激波、非定常效应的存在导致无法测量或数据失效的方案，数值模拟无法计算的方案是指计算无法收敛等无法计算的方案；受客观条件限制进行的有限方案是指：客观条件包括实验条件、计算资源、时间、人力，由于条件有限，对于试验设计的全部方案无法全部进行试验，在有限的实验成本内，只能完成数量有限的试验方案；步骤3：根据方案确定影响因素变化范围，利用空间插值方法补全试...

【专利技术属性】
技术研发人员：安宇飞，刘艳明，张希，唐耀璇，王大磊，
申请(专利权)人：北京动力机械研究所，
类型：发明
国别省市：