一种高精度非接触气动摩擦阻力测量方法及测量装置的制造方法
【技术领域】
[0001] 本发明涉及激光测速技术领域中的粒子图像测速（ParticleImage Velocimetry，简称PIV)技术，具体来说，是一种高精度非接触式气动摩擦阻力测量方法及 测量装置。
【背景技术】
[0002] 对于航空航天飞行器而言，由气流相对运动引起的壁面摩擦阻力是影响飞行器空 气动力学特性的重要力学指标，是在地面风洞试验中需要准确测量的物理量。测量壁面摩 擦阻力的方法包括：牛顿内摩擦定律方法，边界层动量积分方法，Preston管法，Stanton管 法，热比拟方法，油膜干涉法和液晶涂层法等。这些方法各有利弊，简述如下：直接测量法基 于牛顿内摩擦定律，通过测量壁面处平均流向速度的法向梯度来计算壁面摩擦阻力，虽然 容易实现，但受测速仪器分辨率的限制，测量精度较低；边界层动量积分方法通过对速度型 (平均流向速度的法向分布）积分可以有效消除部分测量误差，从而提高精度，但是测量工 作量较大，至少需要测量两个流向站点的边界层速度型；Preston管和Stanton管实际上是 Pitot管的一种发展和延伸，通过测量壁面外某点的平均流速来估计壁面处的速度梯度，虽 然操作容易，但是测速管道的引入会对流场产生一定的干扰，测量精度不高；基于MEMS技 术的微型传感器尺寸小、惯性小，一定程度上提高了测量的时间和空间分辨率，但传感器加 工困难，成本较高；油膜干涉法利用薄油膜在剪切应力作用下厚度发生变化的特性实现对 剪切应力的测量，但标定复杂，且对光路系统有很高要求；液晶涂层法虽然操作方便，但测 试精度很低。此外，上述方法都引入了测试探头或涂层与流场进行接触测量，不可避免的对 测试流场产生干扰，从而影响测试精度。近年发展起来的PIV技术以非接触光学测量的方 法测量靠近壁面的流场速度，能够解决探头/涂层干扰问题。在PIV测速的基础上应用牛 顿内摩擦定律，即可获得壁面摩擦阻力。由于壁面摩擦阻力计算依赖于壁面处流向速度的 法向梯度，因此提高摩擦阻力测量精度的关键是得到空间分辨率较高的速度场。
[0003] 传统的PIV图像处理算法因为查询窗口的存在而面临空间分辨率不足的难题，具 体方式为：
[0004] 应用传统PIV互相关算法计算跨帧图像对中示踪粒子的运动速度时，通常采用窗 口互相关方法，查询窗口大小一般选择32像素X32像素。假设查询窗口大小为2MX2N，则 相关函数为：
[0005]
：￡1>
[0006] 式中，I为t时刻图像灰度，I'为t+At时刻图像灰度。查询窗口是一个2MX2N 大小的矩形，（i，j)表示在矩形中位置坐标为（i，j)处的像素，假设坐标原点设在矩形的中 心处；（Ar,As)表示t+At时刻图像上与（i,j)坐标点相距（Ar,As)处的像素。
[0007] 通过在相关区（根据来流速度，实验的视场范围以及微距镜头与测试平面的距离 等条件确定）寻找使得相关函数取得最大值时的（X，y)，得到示踪粒子组在At时间内的位 A.y 移As，根据速度等于位移除时间，即便可得到查询窗口区域内速度的平均值，此 时的空间分辨率为32像素。
[0008] 但由于在近壁区速度梯度急剧变化，查询窗口的存在会导致速度估值偏高，从而 造成壁面处速度梯度计算不准。为了将PIV测速原理应用到非接触气动摩擦阻力测量上， 必须解决PIV图像处理算法空间分辨率不高的问题。
【发明内容】
[0009] 本发明提供一种高精度非接触气动摩擦阻力测量方法及装置，用于在空气动力学 风洞试验中通过非接触光学方法，快速测量模型表面气动摩擦阻力的一维分布。
[0010] 本发明一种高精度非接触气动摩擦阻力测量方法，通过下述步骤实现：
[0011] 步骤1 :在流场中播撒示踪粒子。
[0012] 步骤2 :照明流场。
[0013] 步骤3 :待流场稳定。；
[0014] 步骤4 :米集不踪粒子跨帧图像对。
[0015] 对于典型风洞试验，采样时间需要超过2分30秒，采集的粒子跨帧图像对需要超 过2000对；每对示踪粒子跨帧图像对为t与t+At时刻记录下的粒子图像，At为跨帧时 间。
[0016] 步骤5 :对步骤4得到示踪粒子跨帧图像对进行处理，采用单像素系综互相关算法 得到模型表面近壁区平均流向速度场及速度梯度。
[0017] A、对粒子跨帧图像对进行分组；
[0018] 将步骤4中采集的全部示踪粒子跨帧图像对分开成两个组，第一组图像为每对示 踪粒子跨帧图像对中在t时刻记录下的示踪粒子图像，记为setl。第二组图像为每对示踪 粒子跨帧图像对中在t+At时刻记录下的示踪粒子图像，记为set2。
[0019] B、setl与set2中的图像匹配；
[0020] 以灰度强度显示setl中每幅图像中的每一个示踪粒子，并采用单像素系综互相 关方法，在set2中找到与setl中图像对应的图像中的示踪粒子进行匹配，具体方式为：
[0021] 令示踪粒子跨帧图像对中，setl中的示踪粒子跨帧图像中像素点坐标为（i，j)， set2中对应的像素点邻域内相距坐标（Ar，As)的像素点坐标为（i+Ar，j+As)。
[0022] 则Nf个示踪粒子跨帧图像对的系综相关函数用公式表示为：
[0023]
[0024] 式中，/；：,?和分别为第n个示踪粒子跨帧图像对中位于setl中图像与set2 中图像像素点的灰度值；U和是Nf个示踪粒子跨帧图像对中位于setl中图像与set2中图像像素点的灰度强度系综平均值，〇表示灰度强度的标准差：
[0025]
[0026] 得到Nf个示踪粒子跨帧图像对的系综相关函数之后，利用互相关算法找到系综相 关函数R&As(i，j)的峰值，进而得到示踪粒子在setl和set2间的位移，最终得到示踪粒 子的运动速度，且利用中心差分法得到某一法向高度y处的速度梯度。
[0027] 步骤6 :根据牛顿内摩擦定律计算出飞行器模型表面的摩擦阻力。
[0028] 针对上述高精度非接触气动摩擦阻力测量方法的测量装置，包括风洞、激光器、片 光系统、CCD相机、同步控制器与计算机；
[0029] 其中，风洞用来产生并且控制气流；激光器的发射端安装有片光系统，片光系统用 来将激光器发射的激光束扩散成片光；片光与流场的速度方向平行且垂直于飞行器模型需 要进行摩擦阻力测量的表面；CCD相机的光轴与片光平面相垂直；通过CCD相机采集一定时 间内沿飞行器模型需要进行摩擦阻力测量的表面外法线方向一定高度内的示踪粒子的跨 帧图像对，转变成数字信号并存入控制计算机；同步控制器具有三个接口，分别通过连接线 与激光器、CCD相机以及计算机相连接；通过同步控制器接收计算机发送的数字信号，同时 用来触发CCD相机和激光器工作，使得CCD相机可在双曝光模式下采集示踪粒子的跨帧图 像对。
[0030] 1、本发明高精度非接触气动摩擦阻力测量方法及测量装置，在图像采集过程中采 用长焦微距镜头能够以较高的分辨率记录粒子图像，利用单像素系综相关算法进行图像处 理可以得到每一个像素上的时间平均速度，因此具有很高的空间分辨率，为准确获得壁面 摩擦阻力奠定了基础；
[0031] 2、本发明高精度非接触气动摩擦阻力测量方法及测量装置，能够获得更高精度的 壁面摩擦阻力一维分布，并且对测量设备要求较低，测量时间快，对温度、湿度等外界条件 不敏感。
【附图说明】
[0032] 图1为本发明粒子图像测速方法流程图。
[0033] 图2为单像素精度算法示意图。
[0034] 图3为根据牛顿内摩擦定律计算壁面摩擦阻力的原理示意图。
[0035] 图4为本发明在单像素精度算法实施例中计算得出的零压力梯度下典型二维湍 流边界层的流向平均速度分布U+ (y+)和流向速度脉动强度分布u"s+ (y+)。
[0036]图5为本发明在单像素精度算法实施例中于低速风洞试验中测量得到的典型粒 子图像。
[0037] 图6a为本发明在单像素精度算法实施例中采用单像素系综相关算法和传统窗口 相关算法计算的速度绝对误差|%|的壁面法向分布。
[0038] 图6b为本发明在单像素精度算法实施例中采用单像素系综相关算法计算的平均 速度随机误差〇 (%)的壁面法向分布。
[0039]图7为本发明在单像素精度算法实施例中采用单像素系综相关算法和传统窗口 相关算法计算的壁面摩擦系数Cf相对误差lij的对比。
[0040]图8为本发明高精度非接触气动摩擦阻力测量仪结构示意图。
[0041 ]图中：
[0042] 1-风洞 2-激光器 3-片光系统
[0043] 4-C⑶相机5-同步控制器6-计算机
[0044] 7-示踪粒子
【具体实施方式】
[0045] 下面结合附图对本发明作进一步详细说明。
[0046] 本发明高精度非接触气动摩擦阻力测量方法，如图1所示，通过下述步骤实现：
[0047] 步骤1 :在流场中播撒示踪粒子；
[0048] 在PIV实验中，示踪粒子大小及浓度的选择一定程度上影响着测量精度，故在进 行实验前需要根据实验条件及要求选取适当大小的示踪粒子及浓度，通常情况下气体环境 中粒径在1-5ym时综合性能较好；而由本发明通过验证实验已经得知，利用单像素系综相 关算法的最优粒子直径为3~4像素，具体选择还要根据具体实验的视场范围、微距镜头与 测试平面的距离、激光器输出光强等参数进一步确定。
[0049] 步骤2 :采用激光器照明流场。
[0050] 步骤3 :开启风洞待流场稳定。
[0051] 步骤4 :通过C⑶相机，采集示踪粒子跨帧图像对；
[0052] 对于典型风洞试验，C⑶相机采样时间需要超过2分30秒，采集的粒子跨帧图像 对需要超过2000对；每对示踪粒子跨帧图像对为CCD相机曝光t与t+At时刻记录下的粒 子图像，At为跨帧时间。上述CCD相机曝光的间隔时间和跨帧时间的选择要根据流场的 流速以及视野范围来选择，选择的依据就是使得示踪粒子在跨帧时间内移动的像素要大于 粒子在图像中所占据的像素个数，小于10个像素大小左右。
[0053] 步骤5:将由步骤4得到示踪粒子跨帧图像对，导入计算机进行数据处理，采用单 像素系综互相关算法进行计算，得到模型表面近壁区平均流向速度场及速度梯度；
[0054] A、对粒子跨帧图像对进行分组；
[0055] 将由步骤4中由CCD相机采集的全部示踪粒子跨帧图像对分开成两个组，如图3 所示，第一组图像为每对示踪粒子跨帧图像对中CCD相机在曝光t时刻记录下的示踪粒子 图像，记为setl;而第二组图像为每对示踪粒