一种油膜干涉表面摩擦应力测量方法技术

本发明专利技术公开了一种油膜干涉表面摩擦应力测量方法，其具体包括以下的步骤：步骤一、在t1时刻和t2时刻之间，获取油膜粘性系数                                               随时间的变化关系，并测量得到t1时刻和t2时刻的油膜厚度；步骤二、设定初始摩擦应力，并根据步骤一得到的油膜粘性系数随时间的变化关系，t1时刻的油膜厚度，数值计算薄油膜方程，并采用显式时间步推进到t2时刻，从而计算得到t2时刻的油膜厚度；步骤三、比较计算得到的t2时刻的油膜厚度与测量得到的t2时刻的油膜厚度，并根据比较结果对摩擦应力进行修正，重新从t1时刻开始显式时间步推进到t2时刻，当计算所得的油膜厚度与测量得到的油膜厚度一致时循环结束，此时的即为求解的摩擦应力。通过上述方法得到油膜干涉表面摩擦应力。

【技术实现步骤摘要】
一种油膜干涉表面摩擦应力测量方法
本专利技术涉及试验流体力学
，本专利技术公开了一种油膜干涉表面摩擦应力测量方法。
技术介绍
摩擦阻力在飞行器的设计中占有非常重要的地位。对于商业飞机而言，在巡航状态时摩擦阻力几乎占到总阻的40％～50％，有评估认为减少1％的阻力大约可以增加10名乘客[1]，无论是从经济效益还是节能环保方面来说，减阻都是非常重要的；对于高超声速带动力飞行器而言，摩擦阻力占到了总阻的30％，直接关系到飞行器是否可以获得正推力以及设计方案是否成功；对于超高声速再入滑翔类武器而言，摩擦阻力的大小直接影响其射程，关系到武器设计的成败。传统的模型表面摩擦应力测量技术有摩阻天平、Preston管等方法，这些方法需要在实验前进行仔细的校准，在实验后对测量数据进行修正，因此整个实验非常繁复，测量误差较大。尤其是在高超声速风洞实验中，由于来流速压较小，模型表面的摩擦应力只有几帕斯卡到几十帕斯卡，基于应变原理的传感器很难测量准确。而Preston管等一些方法需要用探针测量靠近壁面的压力分布，探针直接插入流场，对探针附近的流动不可避免的产生干扰。因此需要发展非插入、更准确的测量方法。近二十年来得到大力发展的测量方法主要有MEMS、剪切液晶和油膜干涉(OilFilmInterferometry)。其中油膜干涉方法能适应较大速度范围的流动情况。它是一种非插入、不需要校准、能应用于曲面的摩擦应力测量方法，分辨率较高。除此之外，它对模型的加工要求不高。油膜干涉方法的主要原理为，覆盖在模型表面的油膜，受边界层的剪切作用会形成楔形形状，其厚度随时间的变化率主要与表面摩擦应力有关。油膜的厚度通常在微米量级，对流场基本没有干扰，可以通过光学干涉的方法准确测量。当入射光(波长为λ)照射模型表面时，分别在空气/油膜界面和油膜/模型界面上发生反射。两道反射光在无穷远或相机透镜的聚焦下相互干涉，光强相长或相消，表现为黑白相间的条纹。1962年Squire提出的薄油膜公式给出了油膜厚度与摩擦应力的关系(见公式1)。Tanner等人首次将其应用于干涉法摩擦应力测量。Zilliac提出用油滴代替油线，将之应用于三维流动，用前缘附近两三道条纹求解摩擦应力大小，用油滴流动方向指示摩擦应力的方向。传统的油膜干涉摩擦应力测量技术主要应用于低速和亚跨声速范围，在整个吹风过程中油膜的粘性系数变化较小，通常被当成一个常数进行处理。只需要精确测量模型表面目标区域内一点在某时刻的温度，即可获得用于求解的粘性系数。有些解算方法虽然考虑了油膜的粘性系数随时间的变化情况，但在实际应用中需要连续测量油膜的温度并进行积分。这一过程非常繁琐，且实现起来很困难。如Garrison和Ackman提出的迭代求解方法，计算公式如下：其中Cf是摩擦应力系数，n是流线间距离，q是来流动压，i是迭代步数，迭代的初始值使用下式。(2)式使用的是单时间步求解方法。对于常规直加热式高超声速风洞而言，风洞启动时间较长，难以准确确定油膜变形的开始时间，且由于启动和关车阶段气流冲击载荷较大，使用单时间步求解方法误差较大。因此，常规高超声速风洞宜采用双时间步求解方法，即在吹风过程中，在t1和t2时刻分别记录油膜干涉图像。另外使用红外热像仪记录t1到t2时刻油膜温度的变化情况。由于油膜的厚度在微米量级，可以认为模型表面的温度就是油膜的温度。在高超声速风洞实验的吹风过程中，油膜的温度随空间和时间均是变化的。对于使用油滴法的求解区域(5mm×5mm)，其随时间的变化更为明显。例如，来流M∞＝7.94的平板模型，在实验时间内，油膜温度变化约为0.55℃/s。显然，当双时间步的时间间隔较长时，由于温度的变化引入的误差将非常明显。高超声速风洞实验的摩擦应力系数较小，以平板模型为例，距前缘100mm处约为0.00036，在较短的时间内(<1s)作用在油膜上的效果(油膜厚度的变化)不明显，由相应的干涉图像计算油膜厚度时相对误差较大。减小油膜的粘性系数似乎是一个可行的方法，当油膜粘性系数足够小时，短时间内油膜的厚度也会发生明显变化。但在实际的风洞实验中，当油膜粘性系数较小时，很容易受各种扰动的影响，导致干涉条纹扭曲变形，因而不能获得准确的油膜厚度。具体到数据处理方面，油膜干涉摩擦应力测量技术主要由两部分构成，一是由光学干涉图像获得油膜的厚度，二是由油膜的厚度随时间的变化率获得摩擦应力。前者显然与温度的变化无关，其处理方法可参考任何一本有关光学干涉方面的书籍。后者在求解薄油膜方程时需要考虑油膜粘性系数随时间的变化。以一维情况为例，见公式(4)。双时间步解法一种是将(4)式在t1到t2时间段积分，积分过程中将μ当作常数；一种是做定常摩擦应力假设(如(3)式)，分别测得t1时刻的条纹宽度为d1，t2时刻的条纹宽度为d2，条纹宽度对应的油膜厚度差为其中λ为入射光波长，nf、na分别为硅油和空气的折射率，θi为入射角，则摩擦应力为上述两种方法要么根本没有考虑μ的变化，要么只是简单的取μ的平均值来代替粘性系数的变化。这样的处理带来的误差比较大，尤其是在高超声速实验中。以双时间步解法的第二种方法为例，t1时刻的不同、t1与t2之间间隔的不同会导致不同的结果。
技术实现思路
针对上述的摩擦应力测量方法存在没有考虑粘性系数变化或者对粘性系数变化采用简单平均的方式导致测量结果误差较大的技术问题，本专利技术提出了一种新的油膜干涉表面摩擦应力测量方法。本专利技术的目的通过下述技术方案来实现：一种油膜干涉表面摩擦应力测量方法，其具体包括以下的步骤：步骤一、在t1时刻和t2时刻之间，获取油膜粘性系数μ随时间的变化关系，并测量得到t1时刻和t2时刻的油膜厚度；步骤二、设定初始摩擦应力τ0，并根据步骤一得到的油膜粘性系数μ随时间的变化关系，t1时刻的油膜厚度，数值计算薄油膜方程，并采用显式时间步推进到t2时刻，从而计算得到t2时刻的油膜厚度；步骤三、比较计算得到的t2时刻的油膜厚度与测量得到的t2时刻的油膜厚度，并根据比较结果对摩擦应力τ0进行修正，重新从t1时刻开始显式时间步推进到t2时刻，当计算所得的油膜厚度与测量得到的油膜厚度一致时循环结束，此时的τ0即为求解的摩擦应力τ。通过上述方法得到油膜干涉表面摩擦应力，这样的方法充分结合粘性系数的变化，同时又不需要连续测量油膜的温度，就能够得到准确的油膜干涉表面摩擦应力。更进一步地，上述采用显式时间步推进到t2时刻，从而计算得到t2时刻的油膜厚度具体包括以下的步骤：将一维情况下的薄油膜公式进行离散，时间上采用向前差分，空间上采用中心差分，得到其中n为时间步数，i为空间步数，μn为第n个时间步时油膜的粘性系数，h表示油膜厚度，t表示时间，μ表示硅油的粘性系数，n＝(t2-t1)/Δt，x是物面坐标。更进一步地，上述获取油膜粘性系数μ随时间的变化关系具体步骤为：首先多次测量油膜温度，获得此时间段内硅油油膜温度随时间的变化关系，并根据预先标定的硅油粘性系数和温度曲线获得油膜粘性系数随时间的变化关系。更进一步地，上述步骤二中初始摩擦应力τ0为：其中：μ1、μ2分别为t1和t2时刻油膜的粘性系数，d1、d2分别为t1和t2时刻油膜干涉条纹的间距(即相邻亮条纹中心线间距离)，Δh为一倍条纹间距对应的油本文档来自技高网...

【技术保护点】
一种油膜干涉表面摩擦应力测量方法，其具体包括以下的步骤：步骤一、在t1时刻和t2时刻之间，获取油膜粘性系数随时间的变化关系，并测量得到t1时刻和t2时刻的油膜厚度；步骤二、设定初始摩擦应力，并根据步骤一得到的油膜粘性系数随时间的变化关系，t1时刻的油膜厚度，数值计算薄油膜方程，并采用显式时间步推进到t2时刻，从而计算得到t2时刻的油膜厚度；步骤三、比较计算得到的t2时刻的油膜厚度与测量得到的t2时刻的油膜厚度，并根据比较结果对摩擦应力进行修正，重新从t1时刻开始显式时间步推进到t2时刻，当计算所得的油膜厚度与测量得到的油膜厚度一致时循环结束，此时的即为求解的摩擦应力。

【技术特征摘要】
1.一种油膜干涉表面摩擦应力测量方法，其具体包括以下的步骤：步骤一、在t1时刻和t2时刻之间，获取油膜粘性系数μ随时间的变化关系，并测量得到t1时刻和t2时刻的油膜厚度；步骤二、设定初始摩擦应力τ0，并根据步骤一得到的油膜粘性系数μ随时间的变化关系，t1时刻的油膜厚度，数值计算薄油膜方程，并采用显式时间步推进到t2时刻，从而计算得到t2时刻的油膜厚度；步骤三、比较计算得到的t2时刻的油膜厚度与测量得到的t2时刻的油膜厚度，并根据比较结果对摩擦应力τ0进行修正，重新从t1时刻开始显式时间步推进到t2时刻，当计算所得的油膜厚度与测量得到的油膜厚度一致时循环结束，此时的τ0即为求解的摩擦应力τ。2.如权利要求1所述的一种油膜干涉表面摩擦应力测量方法，其特征在于所述采用显式时间步推进到t2时刻，从而计算得到t2时刻的油膜厚度具体包括以下的步骤：将一维情况下的薄油膜公式进行离散，时间上采用向前差分，空间上采用中心差分，得到其中n为时间步数，i为空间步数，μn为第n个时间步时油膜的粘性系数，h表示油膜厚度，t表示时间，μ表示硅油的粘性系数，n＝(t2-t1)/Δt，x是物面坐标，△x为沿流线方向的条纹间距。3.如权利要求1或者2所述的一种油膜干涉表面摩擦应力测量方法，其特征在于所述获取油膜粘性系数μ随时间的变...

【专利技术属性】
技术研发人员：刘志勇，张长丰，
申请(专利权)人：中国空气动力研究与发展中心高速空气动力研究所，
类型：发明
国别省市：四川;51