一种高精度非接触气动摩擦阻力测量方法及测量装置制造方法及图纸

本发明专利技术公开了一种基于粒子图像测速(PIV)原理的高精度非接触气动摩擦阻力测量方法及测量装置，用于在空气动力学风洞试验中通过非接触光学方法，快速测量模型表面气动摩擦阻力的一维分布，测量仪包括：双脉冲激光器和片光系统、接收和记录粒子散射光的CCD相机、同步控制器与计算机。在风洞试验中使用该测量仪测量得到模型近壁面的粒子图像对序列，其后采用单像素系综互相关算法对其进行图像处理计算，能够得到空间分辨率较高的近壁区平均流向速度场，根据牛顿内摩擦定理，由壁面处的速度梯度计算出精度较高的摩擦力及其一维分布。本发明专利技术具有精度高，对设备要求较低，测量时间快，对温度、湿度等外界条件不敏感等优点。

【技术实现步骤摘要】
一种高精度非接触气动摩擦阻力测量方法及测量装置
本专利技术涉及激光测速
中的粒子图像测速(ParticleImageVelocimetry，简称PIV)技术，具体来说，是一种高精度非接触式气动摩擦阻力测量方法及测量装置。
技术介绍
对于航空航天飞行器而言，由气流相对运动引起的壁面摩擦阻力是影响飞行器空气动力学特性的重要力学指标，是在地面风洞试验中需要准确测量的物理量。测量壁面摩擦阻力的方法包括：牛顿内摩擦定律方法，边界层动量积分方法，Preston管法，Stanton管法，热比拟方法，油膜干涉法和液晶涂层法等。这些方法各有利弊，简述如下：直接测量法基于牛顿内摩擦定律，通过测量壁面处平均流向速度的法向梯度来计算壁面摩擦阻力，虽然容易实现，但受测速仪器分辨率的限制，测量精度较低；边界层动量积分方法通过对速度型(平均流向速度的法向分布)积分可以有效消除部分测量误差，从而提高精度，但是测量工作量较大，至少需要测量两个流向站点的边界层速度型；Preston管和Stanton管实际上是Pitot管的一种发展和延伸，通过测量壁面外某点的平均流速来估计壁面处的速度梯度，虽然操作容易，但是测速管道的引入会对流场产生一定的干扰，测量精度不高；基于MEMS技术的微型传感器尺寸小、惯性小，一定程度上提高了测量的时间和空间分辨率，但传感器加工困难，成本较高；油膜干涉法利用薄油膜在剪切应力作用下厚度发生变化的特性实现对剪切应力的测量，但标定复杂，且对光路系统有很高要求；液晶涂层法虽然操作方便，但测试精度很低。此外，上述方法都引入了测试探头或涂层与流场进行接触测量，不可避免的对测试流场产生干扰，从而影响测试精度。近年发展起来的PIV技术以非接触光学测量的方法测量靠近壁面的流场速度，能够解决探头/涂层干扰问题。在PIV测速的基础上应用牛顿内摩擦定律，即可获得壁面摩擦阻力。由于壁面摩擦阻力计算依赖于壁面处流向速度的法向梯度，因此提高摩擦阻力测量精度的关键是得到空间分辨率较高的速度场。传统的PIV图像处理算法因为查询窗口的存在而面临空间分辨率不足的难题，具体方式为：应用传统PIV互相关算法计算跨帧图像对中示踪粒子的运动速度时，通常采用窗口互相关方法，查询窗口大小一般选择32像素×32像素。假设查询窗口大小为2M×2N，则相关函数为：式中，I为t时刻图像灰度，I'为t+Δt时刻图像灰度。查询窗口是一个2M×2N大小的矩形，(i，j)表示在矩形中位置坐标为(i,j)处的像素，假设坐标原点设在矩形的中心处；(Δr,Δs)表示t+△t时刻图像上与(i,j)坐标点相距(Δr,Δs)处的像素。通过在相关区(根据来流速度，实验的视场范围以及微距镜头与测试平面的距离等条件确定)寻找使得相关函数取得最大值时的(x,y)，得到示踪粒子组在Δt时间内的位移Δs，根据速度等于位移除时间，即便可得到查询窗口区域内速度的平均值，此时的空间分辨率为32像素。但由于在近壁区速度梯度急剧变化，查询窗口的存在会导致速度估值偏高，从而造成壁面处速度梯度计算不准。为了将PIV测速原理应用到非接触气动摩擦阻力测量上，必须解决PIV图像处理算法空间分辨率不高的问题。
技术实现思路
本专利技术提供一种高精度非接触气动摩擦阻力测量方法及装置，用于在空气动力学风洞试验中通过非接触光学方法，快速测量模型表面气动摩擦阻力的一维分布。本专利技术一种高精度非接触气动摩擦阻力测量方法，通过下述步骤实现：步骤1：在流场中播撒示踪粒子。步骤2：照明流场。步骤3：待流场稳定。；步骤4：采集示踪粒子跨帧图像对。对于典型风洞试验，采样时间需要超过2分30秒，采集的粒子跨帧图像对需要超过2000对；每对示踪粒子跨帧图像对为t与t+Δt时刻记录下的粒子图像，Δt为跨帧时间。步骤5：对步骤4得到示踪粒子跨帧图像对进行处理，采用单像素系综互相关算法得到模型表面近壁区平均流向速度场及速度梯度。A、对粒子跨帧图像对进行分组；将步骤4中采集的全部示踪粒子跨帧图像对分开成两个组，第一组图像为每对示踪粒子跨帧图像对中在t时刻记录下的示踪粒子图像，记为set1。第二组图像为每对示踪粒子跨帧图像对中在t+Δt时刻记录下的示踪粒子图像，记为set2。B、set1与set2中的图像匹配；以灰度强度显示set1中每幅图像中的每一个示踪粒子，并采用单像素系综互相关方法，在set2中找到与set1中图像对应的图像中的示踪粒子进行匹配，具体方式为：令示踪粒子跨帧图像对中，set1中的示踪粒子跨帧图像中像素点坐标为(i，j)，set2中对应的像素点邻域内相距坐标(Δr,Δs)的像素点坐标为(i+Δr，j+Δs)。则Nf个示踪粒子跨帧图像对的系综相关函数用公式表示为：式中，和分别为第n个示踪粒子跨帧图像对中位于set1中图像与set2中图像像素点的灰度值；和是Nf个示踪粒子跨帧图像对中位于set1中图像与set2中图像像素点的灰度强度系综平均值，σ表示灰度强度的标准差：得到Nf个示踪粒子跨帧图像对的系综相关函数之后，利用互相关算法找到系综相关函数RΔr,Δs(i,j)的峰值，进而得到示踪粒子在set1和set2间的位移，最终得到示踪粒子的运动速度，且利用中心差分法得到某一法向高度y处的速度梯度。步骤6：根据牛顿内摩擦定律计算出飞行器模型表面的摩擦阻力。针对上述高精度非接触气动摩擦阻力测量方法的测量装置，包括风洞、激光器、片光系统、CCD相机、同步控制器与计算机；其中，风洞用来产生并且控制气流；激光器的发射端安装有片光系统，片光系统用来将激光器发射的激光束扩散成片光；片光与流场的速度方向平行且垂直于飞行器模型需要进行摩擦阻力测量的表面；CCD相机的光轴与片光平面相垂直；通过CCD相机采集一定时间内沿飞行器模型需要进行摩擦阻力测量的表面外法线方向一定高度内的示踪粒子的跨帧图像对，转变成数字信号并存入控制计算机；同步控制器具有三个接口，分别通过连接线与激光器、CCD相机以及计算机相连接；通过同步控制器接收计算机发送的数字信号，同时用来触发CCD相机和激光器工作，使得CCD相机可在双曝光模式下采集示踪粒子的跨帧图像对。1、本专利技术高精度非接触气动摩擦阻力测量方法及测量装置，在图像采集过程中采用长焦微距镜头能够以较高的分辨率记录粒子图像，利用单像素系综相关算法进行图像处理可以得到每一个像素上的时间平均速度，因此具有很高的空间分辨率，为准确获得壁面摩擦阻力奠定了基础；2、本专利技术高精度非接触气动摩擦阻力测量方法及测量装置，能够获得更高精度的壁面摩擦阻力一维分布，并且对测量设备要求较低，测量时间快，对温度、湿度等外界条件不敏感。附图说明图1为本专利技术粒子图像测速方法流程图。图2为单像素精度算法示意图。图3为根据牛顿内摩擦定律计算壁面摩擦阻力的原理示意图。图4为本专利技术在单像素精度算法实施例中计算得出的零压力梯度下典型二维湍流边界层的流向平均速度分布U+(y+)和流向速度脉动强度分布urms+(y+)。图5为本专利技术在单像素精度算法实施例中于低速风洞试验中测量得到的典型粒子图像。图6a为本专利技术在单像素精度算法实施例中采用单像素系综相关算法和传统窗口相关算法计算的速度绝对误差的壁面法向分布。图6b为本专利技术在单像素精度算法实施例中采用单像素系综相关算法计算的平均速度随本文档来自技高网...

【技术保护点】
一种高精度非接触气动摩擦阻力测量方法，其特征在于：通过下述步骤实现：步骤1：在流场中播撒示踪粒子；步骤2：照明流场；步骤3：待流场稳定；步骤4：采集示踪粒子跨帧图像对；对于典型风洞试验，采样时间需要超过2分30秒，采集的粒子跨帧图像对需要超过2000对；每对示踪粒子跨帧图像对为t与t+Δt时刻记录下的粒子图像，Δt为跨帧时间；步骤5：对步骤4得到示踪粒子跨帧图像对进行处理，采用单像素系综互相关算法得到模型表面近壁区平均流向速度场及速度梯度；A、对粒子跨帧图像对进行分组；将步骤4中采集的全部示踪粒子跨帧图像对分开成两个组，第一组图像为每对示踪粒子跨帧图像对中在t时刻记录下的示踪粒子图像，记为set1；第二组图像为每对示踪粒子跨帧图像对中在t+Δt时刻记录下的示踪粒子图像，记为set2；B、set1与set2中的图像匹配；以灰度强度显示set1中每幅图像中的每一个示踪粒子，并采用单像素系综互相关方法，在set2中找到与set1中图像对应的图像中的示踪粒子进行匹配，具体方式为：令示踪粒子跨帧图像对中，set1中的示踪粒子跨帧图像中像素点坐标为(i，j)，set2中对应的像素点邻域内相距坐标(Δr,Δs)的像素点坐标为(i+Δr，j+Δs)；则Nf个示踪粒子跨帧图像对的系综相关函数用公式表示为：RΔr,Δs(i,j)=1Nf-1×Σn=1Nf(Ii,jn-Ii,j‾)(Ii+Δr,j+Δsn-Ii+Δr,j+Δs‾)σi,jσi+Δr,j+Δs]]>式中，和分别为第n个示踪粒子跨帧图像对中位于set1中图像与set2中图像像素点的灰度值；和是Nf个示踪粒子跨帧图像对中位于set1中图像与set2中图像像素点的灰度强度系综平均值，σ表示灰度强度的标准差：σi,j=1Nf-1Σn=1Nf(Ii,jn-Ii,j‾).]]>得到Nf个示踪粒子跨帧图像对的系综相关函数之后，利用互相关算法找到系综相关函数RΔr,Δs(i,j)的峰值，进而得到示踪粒子在set1和set2间的位移，最终得到示踪粒子的运动速度，且利用中心差分法得到某一法向高度y处的速度梯度；步骤6：根据牛顿内摩擦定律计算出飞行器模型表面的摩擦阻力。...

【技术特征摘要】
1.一种高精度非接触气动摩擦阻力测量方法，其特征在于：通过下述步骤实现：步骤1：在流场中播撒示踪粒子；步骤2：照明流场；步骤3：待流场稳定；步骤4：采集示踪粒子跨帧图像对；对于典型风洞试验，采样时间需要超过2分30秒，采集的粒子跨帧图像对需要超过2000对；每对示踪粒子跨帧图像对为t与t+Δt时刻记录下的粒子图像，Δt为跨帧时间；步骤5：对步骤4得到示踪粒子跨帧图像对进行处理，采用单像素系综互相关算法得到模型表面近壁区平均流向速度场及速度梯度；A、对粒子跨帧图像对进行分组；将步骤4中采集的全部示踪粒子跨帧图像对分开成两个组，第一组图像为每对示踪粒子跨帧图像对中在t时刻记录下的示踪粒子图像，记为set1；第二组图像为每对示踪粒子跨帧图像对中在t+Δt时刻记录下的示踪粒子图像，记为set2；B、set1与set2中的图像匹配；以灰度强度显示set1中每幅图像中的每一个示踪粒子，并采用单像素系综互相关方法，在set2中找到与set1中图像对应的图像中的示踪粒子进行匹配，具体方式为：令示踪粒子跨帧图像对中，set1中的示踪粒子跨帧图像中像素点坐标为(i，j)，set2中对应的像素点邻域内相距坐标(Δr,Δs)的像素点坐标为(i+Δr，j+Δs)；则Nf个示踪粒子跨帧图像对的系综相关函数用公式表示为：式中，和分别为第n个示踪粒子跨帧图像对中位于set1中图像与set2中图像像素点的灰度值；和是Nf个示踪粒子跨帧图像对中位于set1中图像与set2中图像像素点的灰度强度系综平均值，σ表示灰度强度的标准差：得到Nf个示踪粒子跨帧图像对的系综相关函数之后，利...

【专利技术属性】
技术研发人员：潘翀，王建杰，申俊琦，王晋军，李志波，李磊，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11