CN111487440A - 一种五孔探针的标定方法 - Google Patents

一种五孔探针的标定方法 Download PDF

Info

Publication number
CN111487440A
CN111487440A CN202010291024.2A CN202010291024A CN111487440A CN 111487440 A CN111487440 A CN 111487440A CN 202010291024 A CN202010291024 A CN 202010291024A CN 111487440 A CN111487440 A CN 111487440A
Authority
CN
China
Prior art keywords
mach number
target
characteristic
position information
mach
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.)
Withdrawn
Application number
CN202010291024.2A
Other languages
English (en)
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.)
Dalian Linghai Huawei Technology Service Co ltd
Original Assignee
Dalian Linghai Huawei Technology Service Co ltd
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 Dalian Linghai Huawei Technology Service Co ltd filed Critical Dalian Linghai Huawei Technology Service Co ltd
Priority to CN202010291024.2A priority Critical patent/CN111487440A/zh
Publication of CN111487440A publication Critical patent/CN111487440A/zh
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P21/00—Testing or calibrating of apparatus or devices covered by the preceding groups
    • G01P21/02—Testing or calibrating of apparatus or devices covered by the preceding groups of speedometers
    • G01P21/025—Testing or calibrating of apparatus or devices covered by the preceding groups of speedometers for measuring speed of fluids; for measuring speed of bodies relative to fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Measuring Fluid Pressure (AREA)

Abstract

本发明公开了一种五孔探针的标定方法，所述方法基于若干标定马赫数的特征曲线，分别确定马赫数与各特征参数的函数关系；基于所述马赫数与各特征参数的函数关系，确定所述待标定范围内各马赫数的特征曲线。本发明通过在每个马赫数的若干角度位置信息下进行标定，对于每个角度位置信息，基于各马赫数在该角度位置信息下的特征参数，确定该角度信息下马赫数与特征参数对应的函数关系，并基于该函数关系确定待标定马赫数在各角度位置信息的特征参数，测量时马赫数由初始值通过振荡逼近真实值，直到满足精度要求。这样在保证精度的前提下，实现了对一个马赫区间的标定，减少了五孔探针在使用前的标定工作量，扩大了五孔探针的使用范围。

Description

一种五孔探针的标定方法
技术领域
本发明涉及工程测量技术领域，特别涉及一种五孔探针的标定方法。
背景技术
五孔探针是测量复杂三维流场标量和矢量特性的常用手段，因其精度高，设备简单，成本低等优点被广泛应用。现有的基于五孔探针的流场测量技术采用线性插值法，但在实际的流场测量中，线性插值方法对于变化的流场没有很好的适应性，Ma超出标定Ma一定范围将存在较大误差，实际测量中往往需进行大量的标定实验，人力时间成本高。
发明内容
本发明要解决的技术问题在于，针对现有技术的不足，提供一种五孔探针的标定方法。
为了解决上述技术问题，本发明所采用的技术方案如下：
一种五孔探针的标定方法，所述标定方法包括：
基于若干标定马赫数的特征曲线，分别确定马赫数与各特征参数的函数关系，其中，若干标定马赫数均处于待标定范围，各特征参数包括俯仰方向校正系数、偏转方向校正系数、总压校准系数以及静压校准系数；
基于所述马赫数与各特征参数的函数关系，确定所述待标定范围内各马赫数的特征曲线。
所述五孔探针的标定方法，其中，所述若干标定马赫数中各标定马赫数包括若干角度位置信息的特征参数，各标定马赫数包括的若干角度位置信息一一对应，并且若干角度位置信息中的每个角度位置信息均对应有马赫数与该各角度位置信息的各特征参数的函数关系，其中，所述角度位置信息包括俯仰角和偏转角。
所述五孔探针的标定方法，其中，所述马赫数与各特征参数的函数关系均为一元多次函数，并且所述各特征参数对应的一元多次函数的拟合阶数相同。
所述五孔探针的标定方法，其中，所述基于所述马赫数与各特征参数的函数关系，确定所述待标定范围内各马赫数的特征曲线具体包括：
对于预设范围内的一待标定马赫数，获取角度位置信息对应的所述马赫数与各特征参数的函数关系，其中，所述角度位置信息包括俯仰角和偏转角；
对于每个角度位置信息，根据该角度位置信息对应的马赫数与各特征参数的函数关系，确定角度位置信息对应的各特征参数；
基于获取到所有角度位置信息对应的各特征参数，确定该待定马赫数的特征曲线。
一种五孔探针数据处理方法，所述方法包括：
对于测量点的测量数据，分别记录五孔探针的五孔压力值，并所述五孔压力值计算该测量点对应的各特征参数，其中，各特征参数包括俯仰方向校正系数、偏转方向校正系数、总压校准系数以及静压校准系数；
设定初始马赫数，根据所述初始马赫数以及特征参数对应的特征曲线确定目标总压和目标静压，并利用气体关系转换成目标马赫数；
若目标马赫数与初始马赫数的差值大于预设差值阈值，则将目标马赫数作为初始马赫数，根据所述俯仰角以及偏转角选取马赫数与各特性参数的函数关系，并基于选取到的函数关系确定初始马赫数对应的俯仰方向校正系数以及偏转方向校正系数；
继续执行根据所述初始马赫数以及各特性参数确定目标总压和目标静压的步骤，直至所述差值小于或等于预设差值阈值。
所述五孔探针数据处理方法，其中，若目标马赫数与初始马赫数的差值小于或等于预设差值阈值，则将所述目标马赫数作为该测量数据对应的马赫数。
所述五孔探针数据处理方法，其中，所述方法包括：
基于该测量数据对应的马赫数，确定该测量点对应的场流参数。
所述五孔探针数据处理方法，其中，所述所述设定初始马赫数，根据所述初始马赫数以及特征参数对应的特征曲线确定目标总压和目标静压，并利用气体关系转换成目标马赫数具体包括：
设定初始马赫数，根据所述初始马赫数、俯仰方向校正系数以及偏转方向校正系数确定俯仰角以及偏转角；
根据所述俯仰角以及所述偏转角确定目标总压特性系数以及目标静压特性系数；
根据目标总压特性系数和目标静压特性系数计算目标总压和目标静压，并利用气体关系转换成目标马赫数。
一种流场测试系统，其包括上位机以及测试组件；所述测试组件与所述上位机相连接，所述测试组件用于采集五孔压力数据；所述上位机用于执行如上任一所述的五孔探针数据处理方法的步骤。
一种流场测试系统，所述上位机用于执行如上任一所述的五孔探针的标定方法。
有益效果：与现有技术相比，本发明公开了一种五孔探针的标定方法，所述方法基于若干标定马赫数的特征曲线，分别确定马赫数与各特征参数的函数关系；基于所述马赫数与各特征参数的函数关系，确定所述待标定范围内各马赫数的特征曲线。本发明通过在每个马赫数的若干角度位置信息下进行标定，对于每个角度位置信息，基于各马赫数在该角度位置信息下的特征参数，确定该角度信息下马赫数与特征参数对应的函数关系，并基于该函数关系确定待标定马赫数在各角度位置信息的特征参数，测量时马赫数由初始值通过振荡逼近真实值，直到满足精度要求。这样在保证精度的前提下，实现了对一个马赫区间的标定，减少了五孔探针在使用前的标定工作量，扩大了五孔探针的使用范围。
附图说明
图1为本发明提供的五孔探针的标定方法的流程图。
图2为本发明提供的五孔探针的结构原理图。
图3为本发明提供的五孔探针的标定方法中一阶、三阶以及五阶对若干标定马赫数对应的特征参数的拟合曲线。
图4为本发明提供的五孔探针的标定方法中二阶、四阶以及六阶对若干标定马赫数对应的特征参数的拟合曲线。
图5为本发明提供的五孔探针的标定方法中马赫数对应的方向校准系数特征曲线的示意图。
图6为本发明提供的五孔探针的标定方法中方向校准系数特征曲线的拟合插值图。
图7为本发明提供的五孔探针数据处理方法中总压校准系数特征曲线以及静压校准系数特征曲线的示意图。
图8为本发明提供的五孔探针数据处理方法的流程示意图。
具体实施方式
本发明提供一种五孔探针的标定方法，为使本发明的目的、技术方案及效果更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
本技术领域技术人员可以理解，除非特意声明，这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是，本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件，但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解，当我们称元件被“连接”或“耦接”到另一元件时，它可以直接连接或耦接到其他元件，或者也可以存在中间元件。此外，这里使用的“连接”或“耦接”可以包括无线连接或无线耦接。这里使用的措辞“和/或”包括一个或更多个相关联的列出项的全部或任一单元和全部组合。
本技术领域技术人员可以理解，除非另外定义，这里使用的所有术语(包括技术术语和科学术语)，具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是，诸如通用字典中定义的那些术语，应该被理解为具有与现有技术的上下文中的意义一致的意义，并且除非像这里一样被特定定义，否则不会用理想化或过于正式的含义来解释。
如图1和图2所示，本实施提供了一种五孔探针的标定方法，所述方法可以包括以下步骤：
S10、基于若干标定马赫数的特征曲线，分别确定马赫数与各特征参数的函数关系，其中，若干标定马赫数均处于待标定范围，各特征参数包括俯仰方向校正系数、偏转方向校正系数、总压校准系数以及静压校准系数。
具体地，所述若干标定马赫数中各马赫数互不相同，并且各标定马赫数均处于待标定范围内。其中，所述待标定范围为马赫数范围，例如，0-1等。在一个实现方式中，所述若干标定马赫数中包括所述马赫数范围的起始马赫数以及结束马赫数，以便于通过插值法确定待标定范围内各马赫数对应的特征曲线时，可以减少插值和实际值的误差，提高拟合效果。
所述若干标定马赫数均若干角度位置信息的特征参数，其中，所述角度位置信息包括俯仰角和偏转角。可以理解的是，对于每个马赫数，该马赫数对应一个特征参数集，特征参数集中包括若干组特征参数，每组特征参数确定五孔探针的一个角度位置信息，即每个角度位置信息对应一个组特征参数。此外，各标定马赫数包括的若干角度位置信息一一对应，例如，若干标定马赫数包括0.2马赫和0.5马赫，0.2马赫对应有角度位置信息(α1，β1)和角度位置信息(α2，β2)，其中，α1，α2表示偏转角，β1、β2表示俯仰角，那么0.5马赫对应有角度位置信息(α1，β1)和角度位置信息(α2，β2)。相应的，0.2马赫对应有特征参数(Kα1，Kβ1，CPt1，CPs1)和特征参数(Kα2，Kβ2，CPt2，CPs2)，其中，Kα1，Kα2表示偏转方向校正系数，Kβ1、Kβ2表示俯仰方向校正系数，CPt1、CPt2表示总压校准系数，CPs1、CPs2表示静压校准系数，那么0.5马赫对应有特征参数(Kα1，Kβ1，CPt1，CPs1)和特征参数(Kα2，Kβ2，CPt2，CPs2)，当然，值得说明的是，0.2马对应的特征参数与0.5马赫对应的特征参数可以不同，这里就是用符号表示说明，不并表示0.2马赫对应的特征参数和0.5马赫对应的特征参数相同。
进一步，所述马赫数与各特征参数的函数关系均为一元多次函数，并且所述各特征参数对应的一元多次函数的拟合阶数相同。可以理解的是，对于每个角度位置信息，马赫数与所述偏转方向校正系数确定一元多次函数，马赫数与所述俯仰方向校正系数确定一元多次函数，以形成一组一元多次函数，其中，偏转方向校正系数对应的一元多次函数、俯仰方向校正系数对应的一元多次函数、总压校准系数对应的一元多次函数以及静压校准系数对应的一元多次函数的拟合阶数相同。可以理解的是，对于每个角度位置信息，将该角度位置信息对应的一元多次函数组是根据各标定马赫数中该角度位置信息对应的特征参数进行拟合得到。由此，基于若干标定马赫数确定的一元多次函数为一元多次函数集，一元多次函数集中包括若干一元多次函数组，每组一元多次函数组对应一个角度位置信息。
进一步，在本实施例的一个实现方式中，所述一元多次函数可以通过最小二乘法拟合得到，即对于每个角度位置信息下的若干标定马赫以及各标定马赫的特征参数(Kα，Kβ，CPt，CPs)，利用最小二乘法寻找特征参数与马赫数之间的对应关系，以建立马赫数与特征参数之间的函数关系，其中，Kα表示偏转方向校正系数，Kβ表示俯仰方向校正系数，CPt表示总压校准系数，CPs表示静压校准系数。
进一步，所述一元多次函数可以表示为：
Y＝a0+a1n+a2n2+a3n3+...+aknk
其中，a0，a1，a2，a3，...，ak为多项式系数，k为拟合阶数，n为马赫数，Y为俯仰方向校正系数Kβ、偏转方向校正系数Kα，总压校准系数CPt以及静压校准系数CPs。
进一步，由于一元多次函数的拟合阶数k对拟合精度影响较大，从而拟合阶数k需要满足如下原则：1、不选高阶，因为高阶多项式虽然更贴合数据点，但在数据点之间存在较大振荡，数据的拟合效果并不好，插值结果不准确；2、不选低阶，因为阶数太低，数据点离曲线较远，插值结果与实际值也偏差较大。由此，在本实施例的一个实现方式中，所述拟合阶数可以为2阶、3阶或4阶中一种。例如，如图3和4所示，对一组六个马赫数下同Kα、Kβ点的特征值进行一阶、二阶、三阶、四阶、五阶以及六阶拟合，可以得到一阶和六阶的拟合效果不好。
此外，在实际应用中，所述阶数可以根据确定拟合函数的马赫数的数量来确定，以使得拟合得到的一元多次函数精度更高，并且拟合阶数可以随着标定马赫数的增加而增加。例如，对0.2、0.8两个马赫数下所对应特性曲线进行一阶拟合；对于0.2、0.5、0.8三个马赫数下所对应特性曲线进行二阶拟合；对于0.2、0.4、0.7、0.8四个马赫数下所对应特性曲线进行三阶拟合；对于0.2、0.4、0.5、0.7、0.8五个马赫数下所对应特性曲线进行三阶拟合；对于0.2、0.3、0.4、0.5、0.7、0.8六个马赫数下所对应特性曲线进行四阶拟合。在本实施例的一个具体实现方式中，权衡考虑节约时间与保证精度时，在亚声速流场测量区间，对三个马赫数下的特性曲线进行三维非线性拟合便可完成整个区间的标定。
S20、基于所述马赫数与各特征参数的函数关系，确定所述待标定范围内各马赫数的特征曲线。
具体地，对于待标定范围内的各马赫数，基于马赫数与各特征参数的函数关系，确定标定马赫数中各角度位置信息对应的各特征参数。可以理解的是，对于每个待标定范围内的待标定马赫数，分别获取各角度位置信息对应的一元多次函数组，基于获取到的一元多次函数组确定待标定马赫数的特征参数，在获取到所有角度位置信息对应的特征参数后，可以确定待标定马赫数对应的特征曲线，其中，所述特征曲线包括方向校准系数特征曲线，总压校准系数特征曲线以及静压校准系数曲线。例如，对于俯仰角校准系数和偏转角校准系数，将各相邻两个特征参数点相邻即得到待标定马赫数对应的如图6所示的方向校准系数特征曲线。对于总压校准系数以及静压校准系数，总压校准系数特征曲线以及静压校准系数曲线可以如图7所示。当然，值得说明的是，待标定马赫数对应的各角度位置信息是根据基于若干标定马赫数建立的一元多次函数集中各一元多次函数组对应的角度位置信息确定。也就是说，分别将待标定马赫带入一元多次函数集中的各一元多次函数组，即可以得到所有角度位置信息对应的特征参数。例如，如图5所示，0.6马赫、0.8马赫和1马赫为特征曲线为标定马赫的特征区域；0.9马赫的特征曲线为基于马赫数与各特征参数的函数关系生成的方向校准系数特征曲线。
基于上述五孔探针的标定方法，如图8所示，本实施例还提供了一种五孔探针数据处理方法，所述方法还包括：
H10、对于测量点的测量数据，分别记录五孔探针的五孔压力值，并所述五孔压力值计算该测量点对应的各特征参数，其中，各特征参数包括俯仰方向校正系数、偏转方向校正系数、总压校准系数以及静压校准系数；
H20、设定初始马赫数，根据所述初始马赫数以及特征参数对应的特征曲线确定目标总压和目标静压，并利用气体关系转换成目标马赫数；
H30、若目标马赫数与初始马赫数的差值大于预设差值阈值，则将目标马赫数作为初始马赫数，根据所述俯仰角以及偏转角选取马赫数与各特性参数的函数关系，并基于选取到的函数关系确定初始马赫数对应的俯仰方向校正系数以及偏转方向校正系数；
H40、继续执行根据所述初始马赫数以及各特性参数确定目标总压和目标静压的步骤，直至所述差值小于或等于预设差值阈值。
具体地，所述五孔压力值分别记为P1、P2、P3、P4以及P5，在获取到五孔压力值后，可以根据五孔压力值与俯仰方向校正系数的对应关系，以及五孔压力值与偏转方向校正系数的对应关系，计算得到偏转方向校正系数以及俯仰方向校正系数，其中，偏转方向校正系数与五孔压力值的对应关系以及俯仰方向校正系数与五孔压力值的对应关系的表达式可以分别为：
Figure BDA0002450395230000081
其中，
Figure BDA0002450395230000082
为1、3、4以及5孔的平均压力，
Figure BDA0002450395230000083
进一步，所述初始马赫数可以为预先设定，所述初始马赫数可以为若干标定马赫数中，也可为预设标定范围内的任一马赫数。可以理解的是，可以在若干标定马赫数中选取一个标定马赫数作为预先设的初始马赫数，或者在预设标定范围内选取一马赫数作为初始马赫数，当需要对测量点的测量数据进行处理时，采用预先设的初始马赫数作为初始马赫数；也可以当需要对测量点的测量数据进行处理时，在若干标定马赫数中随机选取一个标定马赫数作为初始马赫数。此外，在确定初始马赫数之后，可以根据初始马赫数对应的特性曲线作为初次计算的校准文件，利用线性差值确定来流的俯仰角以及偏转角；并基于俯仰角、偏转角、总压校准系数以及静压校准系数的目标总压和目标静压，以便于根据目标总压和目标静压计算得到目标马赫数。
在本实施例的一个具体实现方式中，所述设定初始马赫数，根据所述初始马赫数以及各特性参数确定俯仰角、偏转角、目标总压和目标静压，并利用气体关系转换成目标马赫数具体包括：
设定初始马赫数，根据所述初始马赫数、俯仰方向校正系数以及偏转方向校正系数确定俯仰角以及偏转角；
根据所述俯仰角以及所述偏转角确定目标总压特性系数以及目标静压特性系数；
根据目标总压特性系数和目标静压特性系数计算目标总压和目标静压，并利用气体关系转换成目标马赫数。
具体地，在设定初始马赫数后，可以确定初始马赫数对应的俯仰方向校正系数以及偏转方向校正系数对应的方向校准系数曲线，以确定偏转角以俯仰角。在确定偏转角和俯仰角后，可以根据初始马赫数对应的总压校准系数曲线以及静压校准系数曲线确定总压校准系数以及静压校准系数。此外，在确定目标总压校准系数和目标静压校准系数后，基于目标总压校准系数和目标静压校准系数的关系式可以确定目标总压和目标静压，其中，目标总压校准系数CPt与目标总压Pt的关系式，以及目标静压校准系数CPs与目标静压Ps的关系式分别为：
Figure BDA0002450395230000091
进一步，在确定目标总压和目标静压后，根据目标总压和目标静压确定目标马赫数的关系式可以为：
Figure BDA0002450395230000092
其中，Ma为马赫数，a、b为系数，例如，a＝5，b＝3.5等；Pa为标准大气压。
进一步，在确定目标总压和目标静压后，根据得到目标总压和目标静压计算得到目标马赫数，基于马赫数与特征参数的函数关系得到该马赫数下的特性曲线，并将该目标马赫数与初始马赫数的差值与预设差值阈值进行比较，若差值大于预设差值阈值，则将该目标马赫数作为初始马赫数，并基于马赫数与特征参数的对应关系更新偏转方向校正系数以及俯仰方向校正系统，并重新执行计算根据所述初始马赫数、俯仰方向校正系数以及偏转方向校正系数确定目标总压特性系数、目标静压特性系数、俯仰角以及偏转角的步骤，以不断震荡逐渐逼近真实的来流马赫数，直至所述差值小于或等于预设差值阈值。
进一步，当偏转角、目标总压若目标马赫数与初始马赫数的差值小于或等于预设差值阈值，则将所述目标马赫数作为该测量数据对应的马赫数；并基于该测量数据对应的马赫数，确定该测量点对应的场流参数。其中，场流参数包括场流速度，所述场流速度的计算公式可以为：
Figure BDA0002450395230000101
其中，
Figure BDA0002450395230000102
t为温度。
进一步，根据间接测量的误差传递公式可以得到场流速度误差σv：
Figure BDA0002450395230000103
其中，σt表示温度误差，
Figure BDA0002450395230000104
为总压误差，
Figure BDA0002450395230000105
为静压误差。
所述总压误差、静压误差、偏转角误差以及俯仰角误差的表达式可以分别为：
Figure BDA0002450395230000106
Figure BDA0002450395230000107
Figure BDA0002450395230000108
Figure BDA0002450395230000109
基于上述实施例，本实施例还提供了一种流场测试系统，其包括上位机以及测试组件；所述测试组件与所述上位机相连接，所述测试组件用于采集五孔压力数据；所述上位机用于执行如上述实施例所述的五孔探针数据处理方法的步骤和/或执行如上述实施例所述的五孔探针的标定方法。
此外，上述上位机的具体工作过程在上述方法中已经详细说明，在这里就不再一一陈述。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种五孔探针的标定方法，其特征在于，所述标定方法包括：
基于若干标定马赫数的特征曲线，分别确定马赫数与各特征参数的函数关系，其中，若干标定马赫数均处于待标定范围，各特征参数包括俯仰方向校正系数、偏转方向校正系数、总压校准系数以及静压校准系数；
基于所述马赫数与各特征参数的函数关系，确定所述待标定范围内各马赫数的特征曲线。
2.根据权利要求1所述五孔探针的标定方法，其特征在于，所述若干标定马赫数中各标定马赫数包括若干角度位置信息的特征参数，各标定马赫数包括的若干角度位置信息一一对应，并且若干角度位置信息中的每个角度位置信息均对应于马赫数与该各角度位置信息的各特征参数的函数关系上，其中，所述角度位置信息包括俯仰角和偏转角。
3.根据权利要求2所述五孔探针的标定方法，其特征在于，所述马赫数与各特征参数的函数关系均为一元多次函数，并且所述各特征参数对应的一元多次函数的拟合阶数相同。
4.根据权利要求1所述五孔探针的标定方法，其特征在于，所述基于所述马赫数与各特征参数的函数关系，确定所述待标定范围内各马赫数的特征曲线具体包括：
对于预设范围内的一待标定马赫数，获取角度位置信息对应的所述马赫数与各特征参数的函数关系，其中，所述角度位置信息包括俯仰角和偏转角；
对于每个角度位置信息，根据该角度位置信息对应的马赫数与各特征参数的函数关系，确定角度位置信息对应的各特征参数；
基于获取到所有角度位置信息对应的各特征参数，确定该待定马赫数的特征曲线。
5.一种五孔探针数据处理方法，其特征在于，所述方法包括：
对于测量点的测量数据，分别记录五孔探针的五孔压力值，并所述五孔压力值计算该测量点对应的各特征参数，其中，各特征参数包括俯仰方向校正系数、偏转方向校正系数、总压校准系数以及静压校准系数；
设定初始马赫数，根据所述初始马赫数以及特征参数对应的特征曲线确定目标总压和目标静压，并利用气体关系转换成目标马赫数；
若目标马赫数与初始马赫数的差值大于预设差值阈值，则将目标马赫数作为初始马赫数，根据所述俯仰角以及偏转角选取马赫数与各特性参数的函数关系，并基于选取到的函数关系确定初始马赫数对应的俯仰方向校正系数以及偏转方向校正系数；
继续执行根据所述初始马赫数以及各特性参数确定目标总压和目标静压的步骤，直至所述差值小于或等于预设差值阈值。
6.根据权利要求5所述五孔探针数据处理方法，其特征在于，若目标马赫数与初始马赫数的差值小于或等于预设差值阈值，则将所述目标马赫数作为该测量数据对应的马赫数。
7.根据权利要求6所述五孔探针数据处理方法，其特征在于，所述方法包括：
基于该测量数据对应的马赫数，确定该测量点对应的场流参数。
8.根据权利要求5所述五孔探针数据处理方法，其特征在于，所述设定初始马赫数，根据所述初始马赫数以及特征参数对应的特征曲线确定目标总压和目标静压，并利用气体关系转换成目标马赫数具体包括：
设定初始马赫数，根据所述初始马赫数、俯仰方向校正系数以及偏转方向校正系数确定俯仰角以及偏转角；
根据所述俯仰角以及所述偏转角确定目标总压特性系数以及目标静压特性系数；
根据目标总压特性系数和目标静压特性系数计算目标总压和目标静压，并利用气体关系转换成目标马赫数。
9.一种流场测试系统，其特征在于，其包括上位机以及测试组件；所述测试组件与所述上位机相连接，所述测试组件用于采集五孔压力数据；所述上位机用于执行如权利要求5-8任一所述的五孔探针数据处理方法的步骤。
10.一种流场测试系统，其特征在于，所述上位机用于执行如权利要求1-4任一所述的五孔探针的标定方法。
CN202010291024.2A 2020-04-14 2020-04-14 一种五孔探针的标定方法 Withdrawn CN111487440A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010291024.2A CN111487440A (zh) 2020-04-14 2020-04-14 一种五孔探针的标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010291024.2A CN111487440A (zh) 2020-04-14 2020-04-14 一种五孔探针的标定方法

Publications (1)

Publication Number Publication Date
CN111487440A true CN111487440A (zh) 2020-08-04

Family

ID=71794885

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010291024.2A Withdrawn CN111487440A (zh) 2020-04-14 2020-04-14 一种五孔探针的标定方法

Country Status (1)

Country Link
CN (1) CN111487440A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114876842A (zh) * 2022-04-29 2022-08-09 沈阳鼓风机集团股份有限公司 一种三维流场测量方法
CN116007884A (zh) * 2023-03-21 2023-04-25 中国空气动力研究与发展中心高速空气动力研究所 一种考虑压缩性的五孔探针球坐标系数据处理方法

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6176130B1 (en) * 1997-06-03 2001-01-23 National Aerospace Laboratory Of Science & Technology Agency Flight velocity vector measuring system in wide velocity region using truncated pyramid-shape probe
US6336060B1 (en) * 2000-02-25 2002-01-01 National Aerospace Laboratory Of Science And Technology Agency Arithmetic processing method and system in a wide velocity range flight velocity vector measurement system using a square truncated pyramid-shape five-hole pitot probe
CN103177192A (zh) * 2013-04-12 2013-06-26 国家电网公司 五孔探针角度特性曲线数据处理方法
CN103868670A (zh) * 2014-03-21 2014-06-18 西北工业大学 一种连续式跨声速风洞实验段流场马赫数控制方法
CN104346522A (zh) * 2013-08-02 2015-02-11 霍尼韦尔国际公司 用于计算马赫数和真空速的系统和方法
CN106096223A (zh) * 2016-05-10 2016-11-09 中国科学院工程热物理研究所 一种五孔探针数据处理方法
CN106989896A (zh) * 2017-04-17 2017-07-28 北京航空航天大学 一种测量亚音三维非定常流场的动态温度压力组合探针
CN108362911A (zh) * 2018-02-09 2018-08-03 西京学院 一种基于最小二乘法和插值法的流速仪标定方法
CN108562420A (zh) * 2017-12-26 2018-09-21 北京航天益森风洞工程技术有限公司 一种多孔探针的自动标定系统及方法
CN108593964A (zh) * 2018-03-22 2018-09-28 中国科学院工程热物理研究所 一种七孔气流探针的标定方法
CN110514390A (zh) * 2019-08-12 2019-11-29 北京航空航天大学 一种三孔压力探针测量二维流场不确定度评定方法
CN110514391A (zh) * 2019-08-12 2019-11-29 北京航空航天大学 一种五孔压力探针测量三维流场不确定度评定方法
CN110568869A (zh) * 2019-09-10 2019-12-13 中国航发沈阳发动机研究所 一种提高控制探针自动跟踪测试精度的控制方法

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6176130B1 (en) * 1997-06-03 2001-01-23 National Aerospace Laboratory Of Science & Technology Agency Flight velocity vector measuring system in wide velocity region using truncated pyramid-shape probe
US6336060B1 (en) * 2000-02-25 2002-01-01 National Aerospace Laboratory Of Science And Technology Agency Arithmetic processing method and system in a wide velocity range flight velocity vector measurement system using a square truncated pyramid-shape five-hole pitot probe
CN103177192A (zh) * 2013-04-12 2013-06-26 国家电网公司 五孔探针角度特性曲线数据处理方法
CN104346522A (zh) * 2013-08-02 2015-02-11 霍尼韦尔国际公司 用于计算马赫数和真空速的系统和方法
CN103868670A (zh) * 2014-03-21 2014-06-18 西北工业大学 一种连续式跨声速风洞实验段流场马赫数控制方法
CN106096223A (zh) * 2016-05-10 2016-11-09 中国科学院工程热物理研究所 一种五孔探针数据处理方法
CN106989896A (zh) * 2017-04-17 2017-07-28 北京航空航天大学 一种测量亚音三维非定常流场的动态温度压力组合探针
CN108562420A (zh) * 2017-12-26 2018-09-21 北京航天益森风洞工程技术有限公司 一种多孔探针的自动标定系统及方法
CN108362911A (zh) * 2018-02-09 2018-08-03 西京学院 一种基于最小二乘法和插值法的流速仪标定方法
CN108593964A (zh) * 2018-03-22 2018-09-28 中国科学院工程热物理研究所 一种七孔气流探针的标定方法
CN110514390A (zh) * 2019-08-12 2019-11-29 北京航空航天大学 一种三孔压力探针测量二维流场不确定度评定方法
CN110514391A (zh) * 2019-08-12 2019-11-29 北京航空航天大学 一种五孔压力探针测量三维流场不确定度评定方法
CN110568869A (zh) * 2019-09-10 2019-12-13 中国航发沈阳发动机研究所 一种提高控制探针自动跟踪测试精度的控制方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HOMAM NIKPEY SOMEHSARAEI: "《A novel approach based on artificial neural network for calibration of multi-hole pressure probes》", 《FLOW MEASUREMENT AND INSTRUMENTATION》 *
TOLGA YASA: "《Robust procedure for multi-hole probe data processing》", 《FLOW MEASUREMENT AND INSTRUMENTATION》 *
廖安文: "《基于改进最小二乘法的五孔探针数据处理方法》", 《东方汽轮机》 *
张晓东: "《五孔探针在涡轮导向器出口流场测量中的应用》", 《燃气涡轮试验与研究》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114876842A (zh) * 2022-04-29 2022-08-09 沈阳鼓风机集团股份有限公司 一种三维流场测量方法
CN116007884A (zh) * 2023-03-21 2023-04-25 中国空气动力研究与发展中心高速空气动力研究所 一种考虑压缩性的五孔探针球坐标系数据处理方法
CN116007884B (zh) * 2023-03-21 2023-05-30 中国空气动力研究与发展中心高速空气动力研究所 一种考虑压缩性的五孔探针球坐标系数据处理方法

Similar Documents

Publication Publication Date Title
CN115950557B (zh) 一种基于压力变送器的温度智能补偿方法
CN106096223B (zh) 一种五孔探针数据处理方法
CN110514391B (zh) 一种五孔压力探针测量三维流场不确定度评定方法
CN108362911B (zh) 一种基于最小二乘法和插值法的流速仪标定方法
CN107271081B (zh) 基于两阶段最小二乘拟合的硅压阻式压力变送器温度补偿方法及装置
CN111487440A (zh) 一种五孔探针的标定方法
CN109142792B (zh) 一种石英挠性加速度计温度误差标定补偿方法
CN105043511B (zh) 基于雷诺数的涡街流量计系数修正方法
CN114323536B (zh) 一种提高五孔探针测量精度的插值方法
CN110987068B (zh) 一种多传感器综合控制系统数据融合方法
CN114486164B (zh) 一种基于神经网络的三孔探针梯度流场测试标定方法
WO2024011779A1 (zh) 群孔测量摆角规划方法、可读介质和设备
CN115358026A (zh) 一种基于多元线性回归与曲面拟合的五孔探针数据处理方法
CN117129179B (zh) 一种连续式风洞翼下双支撑试验的马赫数修正方法
CN114611420A (zh) 非定常气动力计算精度评估及修正方法
CN110553667A (zh) 一种对星敏感器进行热变形补偿的方法
CN108090253B (zh) 一种直升机数字式大气系统空速修正方法
CN110069842B (zh) 舵效快速估算方法
CN108489602B (zh) 一种基于模态测试技术的惯性参数识别方法
CN106645794B (zh) 一种超声速飞行器表面压力测量故障的软冗余方法
CN107860346B (zh) 一种测量坐标系配准方法
CN112965966B (zh) 一种基于实测飞参数据的快速预处理方法、系统及计算机相关产品
CN108151932A (zh) 一种电动执行机构测试系统温度漂移补偿方法
CN110222382B (zh) 一种飞机坐标系优化拟合方法
CN110988787A (zh) 无线信号测向监测中基于聚类分析实现最优测向角测量的方法

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
WW01 Invention patent application withdrawn after publication

Application publication date: 20200804

WW01 Invention patent application withdrawn after publication