一种飞行试验测点位置的流态判读方法专利检索-飞行器飞机类型专利检索查询-专利查询网

一种飞行试验测点位置的流态判读方法

阅读：0发布：2020-06-16

专利汇可以提供一种飞行试验测点位置的流态判读方法专利检索，专利查询，专利分析的服务。并且本发明涉及一种飞行试验测点位置的流态判读方法。该方法的设计原理为：能量守恒与湍流特性。该方法的实现方式为：首先根据飞行试验测量得到的温度曲线T(t)计算得到温升速率T'(t)；然后根据飞行试验弹道和压力测量值计算获得pu(t)曲线；再根据飞行试验弹道获得雷诺数 Re曲线；最后对比T'(t)与pu(t)曲线，当T'(t)开始增加而pu减小时表明该测点转捩进入湍流状态，Re曲线用于辅助判读。本发明基于热功转换原理，利用飞行试验自身测量数据及其衍生数据，无需额外计算，能够实现测点位置流态的快速判读；本发明大大减小了飞行试验后数据处理的时间，基本上可以达到实时处理。，下面是一种飞行试验测点位置的流态判读方法专利的具体信息内容。

权利要求

1.一种飞行试验测点位置的流态判读方法，其特征在于，包括如下步骤：
(1)由飞行试验数据，获取测点处的温度数据T(t)，表面压力数据p(t)，飞行高度H(t)，飞行马赫数Ma(t)，t为记录时刻；
(2)根据飞行高度H(t)计算大气声速c(t)，运动粘性系数ν(t)；
(3)计算获得飞行速度u(t)＝c(t)·Ma(t)，单位雷诺数Re(t)＝u(t)/ν(t)和组合变量W(t)＝p(t)·u(t)；
(4)根据T(t)计算获得温升速率曲线T'(t)＝dT(t)/dt；绘制W(t)曲线；
(5)对比W(t)和T'(t)曲线，当W(t)减小而T'(t)增加时，测点处流态为湍流；当W(t)增加而T'(t)减小时，测点处流态为层流。
2.根据权利要求1所述的一种飞行试验测点位置的流态判读方法，其特征在于：步骤(3)还包括计算单位雷诺数Re(t)＝u(t)/ν(t)；步骤(4)还包括绘制单位雷诺数Re(t)曲线；
还包括步骤(5)判别流态为湍流时，计算当地雷诺数Re＝Re(t)·xt，如果当地雷诺数不大于设定雷诺数阈值，则判断该湍流为误判读，如果大于设定雷诺数阈值，则判读有效，xt为测点所在位置距飞行器端头距离。
3.根据权利要求1或2所述的一种飞行试验测点位置的流态判读方法，其特征在于：步骤(2)中根据飞行高度H(t)计算大气声速c(t)，运动粘性系数ν(t)采用boops大气参数计算软件。
4.根据权利要求1或2所述的一种飞行试验测点位置的流态判读方法，其特征在于：步骤(4)中在Tecplot 绘图软件中导入数据T(t)，W(t)，根据T(t)计算获得温升速率曲线T'(t)＝dT(t)/dt；绘制W(t)曲线。
5.根据权利要求4所述的一种飞行试验测点位置的流态判读方法，其特征在于：步骤(4)还包括利用Tecplot 数据处理中的光滑功能对T(t)进行滤波，使得曲线光顺，再利用其分析功能计算获得温升速率曲线T'(t)＝dT(t)/dt。

说明书全文

一种飞行试验测点位置的流态判读方法

技术领域

[0001] 本发明属于空气动力学技术领域，用于飞行试验测点所在位置的流态判读。

背景技术

[0002] 边界层流态包括层流、湍流以及层流到湍流的转捩(转捩是指层流到湍流的过渡阶段)，由于层流和湍流对质量、动量和能量的输运性质不同，湍流状态下的摩阻和热流是层流时的3-5倍，因此，流态判读是气动特性预测的基础。

[0003] 目前，关于边界层转捩的飞行试验增多，因此，涉及测点所在位置流态的流态判读问题。比较常用的方法是，沿飞行弹道分别开展全层流和全湍流计算，获得层流热流和湍流热流曲线，飞行试验实际测量的曲线一般位于两者之间，把实测曲线偏离层流预测曲线对应的时间视为进入转捩。这种流态判读方法的主要缺点是，需要额外补充大量弹道点的全层流和全湍流计算，计算量较大，耗时长，判读不及时。

发明内容

[0004] 本发明提供一种飞行试验测点位置的流态判读方法，使得工程人员能够根据飞行试验数据快速判读出测点位置的流态，快速获取测点进入转捩和退出转捩的时刻。

[0005] 本发明的技术解决方案是：

[0006] 提供一种飞行试验测点位置的流态判读方法，包括如下步骤：

[0007] (1)由飞行试验数据，获取测点处的温度数据T(t)，表面压力数据p(t)，飞行高度H(t)，飞行马赫数Ma(t)，t为记录时刻；

[0008] (2)根据飞行高度H(t)计算大气声速c(t)，运动粘性系数ν(t)；

[0009] (3)计算获得飞行速度u(t)＝c(t)·Ma(t)，单位雷诺数Re(t)＝u(t)/ν(t)和组合变量W(t)＝p(t)·u(t)；

[0010] (4)根据T(t)计算获得温升速率曲线T'(t)＝dT(t)/dt；绘制W(t)曲线；

[0011] (5)对比W(t)和T'(t)曲线，当W(t)减小而T'(t)增加时，测点处流态为湍流；当W(t)增加而T'(t)减小时，测点处流态为层流。

[0012] 优选的，步骤(3)还包括计算单位雷诺数Re(t)＝u(t)/ν(t)；步骤(4)还包括绘制雷诺数曲线Re(t)；还包括步骤(5)判别流态为湍流时，计算当地雷诺数Re＝Re(t)·xt，如果当地雷诺数不大于设定雷诺数阈值，则判断该湍流为误判读，如果大于设定雷诺数阈值，则判读有效，xt为测点所在位置距飞行器端头距离。

[0013] 优选的，步骤(2)中根据飞行高度H(t)计算大气声速c(t)，运动粘性系数ν(t)采用boops大气参数计算软件。

[0014] 优选的，步骤(4)中在Tecplot 绘图软件中导入数据T(t)，W(t)，根据T(t)计算获得温升速率曲线T'(t)＝dT(t)/dt；绘制W(t)曲线。

[0015] 优选的，步骤(4)还包括利用Tecplot 数据处理中的光滑功能对T(t)进行滤波，使得曲线光顺，再利用其分析功能计算获得温升速率曲线T'(t)＝dT(t)/dt。

[0016] 本发明与现有技术相比的优点在于：

[0017] (1)本发明利用成熟的商业软件进行数据后处理，简单易用；

[0018] (2)本发明基于热功转换原理，利用飞行试验自身测量数据及其衍生数据，无需额外计算，能够实现测点位置流态的快速判读；本发明大大减小了飞行试验后数据处理的时间，基本上可以达到实时处理。

[0019] (3)本发明通过雷诺数进行辅助判断，可以提高流态判断的准确性。附图说明

[0020] 图1为T(t)，T'(t)，Re(t)和W(t)曲线；

[0021] 图2为本发明测点位置流态判读流程图。

具体实施方式

[0022] 对于飞行器表面流态，通常基于热流q的变化进行判读，壁面热流为温度法向梯度，k为导热系数。对于边界层流动，与W(t)正相关，即W(t)越大，越大；T'(t)与q正相关，即q越大，T'(t)越大；因此，如果W(t)减小(即减小)，而T'(t)增加(q增加)，说明只可能是导热系数k增加引起的，代表流态由层流变为湍流，因此，该时刻是测点所在位置层流到湍流的转捩时刻。

[0023] 同样如果W(t)增加，而T'(t)开始减小代表流态由湍流变为层流，对应时刻为测点所在位置湍流到层流的转捩时刻。

[0024] 结合图2，下面以钝楔飞行试验为例对本发明进行更详细的说明：

[0025] (1)从飞行试验获得测点处的温度数据T(t)，表面压力数据p(t)，飞行高度H(t)，飞行马赫数Ma(t)。

[0026] (2)利用boops软件(大气参数计算软件)根据飞行高度H(t)计算给出大气声速c(t)，运动粘性系数ν(t)；

[0027] (3)计算获得飞行速度u(t)＝c(t)·Ma(t)，单位雷诺数Re(t)＝u(t)/ν(t)和组合变量W(t)＝p(t)·u(t)；

[0028] (4)在Tecplot绘图软件中导入数据T(t)，W(t)和Re(t)，利用Tecplot软件的数据处理功能对T(t)数据进行滤波，使数据光滑；再利用该软件的分析功能计算获得温升速率曲线T'(t)＝dT(t)/dt；

[0029] (5)各曲线图如图1所示，对比W(t)和T'(t)曲线，可获得以下结论：

[0030] (a)在t＝21s时，W(t)开始减小，而T'(t)反而增加，这一趋势持续到t＝47s，表明在[21s，47s]区间内，测点处的流态处于湍流状态。

[0031] (b)在t＝47s时T'(t)开始减小，而W(t)先减小后增加，这一趋势持续到t＝59s，表明在[47s，59s]区间内，测点处的流态处于层流状态。

[0032] (c)在t＝59s时T'(t)再次开始增加，而W(t)在减小，这一趋势持续到t＝71s，表明在[59s，71s]区间内，测点处的流态处于湍流状态。

[0033] (d)根据Re(t)曲线进行辅助判读。对于钝楔体，层流到湍流的转捩雷诺数阈值可5 5
取5.0×10 ，即Re(t)·xt<5.0×10时只可能是层流流态，其中xt为测点所在位置距飞行器端头距离。图1中Re(t)的单位为105，分析的测点位置为xt＝0.8m，可以看到，在判读的两个湍流区，均有Re(t)·xt>5.0×105，因此，此次判读是有效的。

[0034] 飞行试验测点位置的流态判读方法，包括如下步骤：

[0035] (1)由飞行试验数据，获取测点处的温度数据T(t)，表面压力数据p(t)，飞行高度H(t)，飞行马赫数Ma(t)，t为记录时刻；

[0036] (2)根据飞行高度H(t)计算大气声速c(t)，运动粘性系数ν(t)；

[0037] (3)计算获得飞行速度u(t)＝c(t)·Ma(t)，单位雷诺数Re(t)＝u(t)/ν(t)和组合变量W(t)＝p(t)·u(t)；

[0038] (4)根据T(t)计算获得温升速率曲线T'(t)＝dT(t)/dt；绘制W(t)曲线；

[0039] (5)对比W(t)和T'(t)曲线，当W(t)减小而T'(t)增加时，测点处流态为湍流；当W(t)增加而T'(t)减小时，测点处流态为层流。

[0040] 步骤(2)中根据飞行高度H(t)计算大气声速c(t)，运动粘性系数ν(t)采用boops大气参数计算软件。

[0041] 步骤(3)还包括计算单位雷诺数Re(t)＝u(t)/ν(t)；步骤(4)还包括绘制单位雷诺数Re(t)曲线；还包括步骤(5)判别流态为湍流时，计算当地雷诺数Re＝Re(t)·xt，如果当地雷诺数不大于设定雷诺数阈值，则判断该湍流为误判读，如果大于设定雷诺数阈值，则判读有效，xt为测点所在位置距飞行器端头距离。

[0042] 步骤(4)中在Tecplot绘图软件中导入数据T(t)，W(t)，根据T(t)计算获得温升速率曲线T'(t)＝dT(t)/dt；绘制W(t)曲线。

[0043] 步骤(4)还包括利用Tecplot数据处理中的光滑功能对T(t)进行滤波，使得曲线光顺，再利用其分析功能计算获得温升速率曲线T'(t)＝dT(t)/dt。

[0044] 本发明大大减小了飞行试验后数据处理的时间，基本上可以达到实时处理。

[0045] 本发明未详细描述内容为本领域技术人员公知技术。

[0046] 以上通过钝楔飞行试验为例描述了本发明的具体实施方式，但并非是对本发明的限制，实际上，只要飞行试验给出了弹道，压力和温度测量数据，均可用本发明所介绍的方法判读测点所在位置流态。还需说明的是，凡是依据本发明的技术实质对以上实例所做的任何简单修改均属于本发明的技术范围。

标题	发布/更新时间	阅读量
一种可调节式扑翼飞行器测力平台及其安装使用方法	2020-10-26	0
用于飞行器的混合网络的模块化元件和飞行器	2021-01-29	1
虚拟应答器	2021-08-27	2
一种高速飞行器强制转捩装置及设计方法	2020-10-12	2
一种灭火飞行器	2022-12-04	0
用于在飞行器机身的区段中安装系统部件的方法	2020-05-25	2
一种时空可变流场下飞行器的球面轨道编队跟踪控制方法	2020-08-06	1
石墨烯型超导无线传输电能飞行器与传感器阵列系统	2021-04-02	0
桥梁检测系统	2020-09-30	2
AIRCRAFT WING SHIFT DEVICE	2023-08-01	0

一种飞行试验测点位置的流态判读方法

一种飞行试验测点位置的流态判读方法

技术领域

背景技术

发明内容

具体实施方式

该功能需要专业版企业版VIP权限，您可以：