空气动力学实验报告

1、实验一边界层流动测量实验摘要：边界层，又称为流动边界、附面层，它是流体流动过程中，紧贴壁而的粘性阻力不可 忽略的一层薄薄的流体，它对主要流体运动的影响很大。口普朗特捉出该概念起，边界层研 究就一直是流体力学研究中一个焦点和难点课题°本实验通过热线风速仪测量距离凹口平板 前缘不同位置点流体的速度分布情况，并对实验数据加以分析处理，从而确定出在不同工况 屮的边界层的厚度、位移厚度，以及避免粘性力等参数，最终分析边界层的特性。关键词：边界层，热线风速仪，粘性力，雷诺数，拟合，标定1. 实验简介此次实验是在一个开口式风洞中进行的，该风洞试验段截面尺寸为：500mm*500mmo 设置风洞风机

2、的运行频率为20hz和30hz、，利用热线风速仪测量凹槽分离点20mm的边界 层上的速度分布。然后用两种不同的方法拟合热线风速仪实验前后标定曲线，得出标定误差 值，从而分析比较这两种拟合方法的优缺点，并分析出实验中热线性能的稳定性。2. 实验步骤1）将皮托管固定在风洞试验段，轴线和来流速度方向平行。记录皮托管标定系数k。皮托 管静压连接到压力传感器负压接口，丿支托管总压连接到压力传感器通道1；2）热线风速仪探头安装在二位坐标架上，连接热线探头与恒温控制器输入、输出。此时热 线恒温控制器切勿通电！将热线探头移至和皮托管同一高度；3）热线输出连接到数据采集卡ai0,皮托管输出连接到数据采集卡all

3、；4）将热线ik温控制器通电，打开matlab热线风速仪标定程序“hw calibnnion.m”，改 变文件名运行程序；5）将热线移动至测量点（距离i叫腔分离点x=20mm）上方口由来流中，调整风洞风速， 风机运彳亍频率f=30hz, matlab运行热线速度分布测量程序"hw measurement.m n改变 文件存储名称。改变风洞风速，风机运行频率f=2（）hz,重复步骤4；6）打开matlab热线风速仪标定程序'hw calibration.m，,改变标定参数存储文件名， 重新运行标定程序。7）按要求处理实验数据。3. 实验设备1）开口式风洞，实验段截面5（）（）m

4、m*5（）（）mim实验段风洞壁面装有凹腔模型，测量模型分 离点上游边界层速度分布。2）程控二位坐标架：计算机软件可精确控制探针在边界层内移动。3）压力传感器：将压力信号转化为电信号，灵敏度为0.0245 v/mm.4）皮托管：测量來流静压几和总压几，用于测量來流速度u和标定热线风速仪。5）热线风速仪:测虽气体流场（尤其是湍流场）中速度大小，具冇较高的频率响应（10khz）o6）数据采集卡：ni-usb6212bnc （16位，8差分通道）；7）计算机：控制二维处标架，控制和执行数据采集；其中，热线风速仪的工作原理如下：热线风速仪前段的热线探头作为后端控制电路中电桥的一个桥端电阻。电桥工作时把

5、探 头上的电阻丝加热，当流体以速度u流过探头时带走一部分热量q电桥设计是为了保持探 头热丝温度恒定，探头热丝损失热量q后，两端电压e发生变化从而提供相应的热量以弥 补损失掉的那部分热暈，最终保证探头热丝的温度怛定。通过标准速度源可以将速度u与 电桥电压e之间的关系拟合出来。这样在测最时，只要测最到电桥电压e,就可以知道当地 风速u。热线风速仪的标定：川风洞自由來流速度作为标准速度，自由來流速度urti皮托管测量得到。拟合关系式 用：1）king" s law： e = a + bun,n = 0.45 ；2）u =a + be + ce2 + de3 + fe4. 实验数据处理1）分

6、别用两种不同方法拟合热线风速仪实验前后标定曲线，给出标定误差值，分析拟合方 法优缺点。将实验前后热线风速仪标定illi线画在同一张图上，将实验后热线标定加上温 度侨正和前面两组曲线应在同一张图上，分析实验中热线性能是否稳定及温度侨止的作 用。表41实验前标定时采集的数据流 il/u电压/e题/t3.08361.696221.2363.99461.742521.2944.94061.782421.3925.86891.817121.5216.7911.847521.6547.7121.876121.767&62041.901521.8729.55781.925621.98910.4681

7、.947322.173表4-2实验后标定时采集的数据流il/u电压/e温度/t3.05831.692422.0293.99231.738822.0094.93641.778922.0335.84981.813322.0916.78111.844422.1167.69911.872622.1918.61191.89822.2829.54051.921622.402按照king's law： e = a + bun = 0.45关系拟合程序：clc;clear;data=load('e:2014标定实验hw calibration-xxxxxx(u_e)_0527.txt'

8、)；%载入数据 v=data (: 1);w=v. a0.45;e=data (:, 2);f=inline (,a(l)+a(2)*w,；a'；w,);b,c=lsqcurvefit(f,1.35,0.205,w,e);%b 是(a,b)向量组合,c 是目标函数值;b,cfor i=l:xv(i)=data(i,l)h(i)=b*l;v(i).a0.45;endfigure;plot(vzh；-ar')xlabelc流体流速u');ylabel('电桥电压e1);title。电压e与流速u的关系j;按照u = a + be + ce2 + de3 + fe4关

9、系拟合程序：clc;clear;data=load('e:2014标定实验hw calibration-xxxxxx(u_e)_0527.txt')；%载入数据 v 二 data(:j);e=data(:,2);f=inline('a(l)+a (2)*e+a(3)*e.a2+a(4)*e.a3+a(5)*e.a4,/'a',e,);b,c=lsqcurvefit(f,l 111 l,e,v);%b 是(a,b)向量组合,c 是目标函数值;b,cfor i=l:xe(i)=data(iz2);h(i)=b*l e(i) e(i).a2 e(i).a3 e

10、(i).a4'endfigure;plot(e,h；-.g')xlabelc电桥电压e');ylabef流体流速it); title。电压e与流速u的关系j;拟合曲线为：由图1可以看出，未拟合前实验前后的曲线间距明显，这就表明，在实验结束后热线探 头发生了一定距离的偏移。而且使用kings law拟合后，拟合前后的曲线还是存在一定的 谋差。根据0=丄£ j仏(data)一坷(mod，实验前llll线的标定误差为2.618%,实验示1111线的标定谋羞为2.433%。31图2拟合方法二标定曲线使用同样的课差公式计算得到：对实验前的曲线的标定谋差为：0.598%,

11、对实验示的曲 线的标定误差为：0.803%。由上图可以看出，使用拟合方法二拟合后的曲线可以与未拟合的曲线比鮫好地重合，说 明拟合误差很小。通过比较上面两图，可以得出，第二种拟合方法明显比第一种方法得到的拟合效果好, 这就说，第二种拟介方法的标定误差较小。同时，从实验前后的两条曲线可以看出，实验前 后热线探头发牛了一定距离的偏移，这证明该实验使用的热线性能还不够稳定。2) 利用文献方法，预测测量点距离壁而粹确距离，测量点摩擦速度，计算壁而粘性应力， 边界层厚度，边界层位移厚度，边界层动最厚度，基于动最厚度的边界层雷诺数等参数。根据文献指导并参考徐进同学的程序,然后进行必要的校匸完善，得到以下程序

12、及图形: clc;clear;file=*e:2014边界层试验x20f20 profile.txt'data=dlmread(file);jun k=sortrows(data);yi=junk(:,l)./1000;u=ju nk(:,2);umea n 二 mean(data(l:20,4);u99=umean*0.99;% j 由来流风速；yd=-9.9:0.01:-4;for j=l:length(yd) %截 llifflim循环迭代;y=yi+yd(j)./1000;ut=0.20:0.01:l;lp=1.5e-5./ut;k=0.41;b=5.0;for i=l:len

13、gth(ut) %壁面摩擦应力循环迭代;yp=y./lp(i);up=u./ut(i);ye=up+exp(-k*b)*(exp(k*up)-l-k*up-(k*up).a2/2-(k*up).a3/6); err(i,j)=sum(abs(yp(l:35)-ye(l:35)./yp(l:35)/35;% 戈余函数计算;endendm indexi=min(err);%行数;n indexj=min(m);%列数；yoff=yd(indexj)% 截止距离;utf=ut(indexi(indexj)%yg 面摩擦应力;lpf=1.5e-5./utf; yp=(yi+yoff./1000)./l

14、pf;upp=l:0.5:25;up=junk(:,2)./utf; ue=(l./k).*log(yp)+b;ye=upp+exp(-k*b)*(exp(k*upp)-l-k*upp-(k*upp).a2/2-(k*upp).a3/6);figure(l);semilogx(yp,up,'-ob7linewidth',l);hold on;semilogx(ye/upp/-or,/linewidth,/l);hold on;semilogx(yp,ue/-k7linewidth',l);xlabel('ya+');ylabel(,ua+')；l

15、egend('experimentai data'z'spalding profile'/log-law');axis(l 1000 0 26);yabs=yi+yoff./1000;ul=l-u/umea n;ya bs p=t ra pz(ya bs, u i). * 1000% 边界层位移厚度;uul=(u/umea n).*(l u ./umea n);t h eta=t ra pz(ya bs, u u i). * 1000% 边界层动最焊度;re=1.2*umean*theta/1000./(1.5e-5)%基于动量厚度的边界层雷诺数； fi

16、gure(2);yreal=yi.*1000+yabsp;ureal=junk(:,2);plot(ureal,yreal,'-oc')xlabelcurearbylabelcyreal');axis(0 8 0 120);収f=20hz时，其风速近似为6.26nvs,经循环迭代得到，测点距壁血的误差距离为 9.34mm；测量点摩擦速度为().34m/s；计算出壁而的剪切应力为0.149n/m2；边界层厚度从 上图可以近似得到，为12mim边界层位移厚度为2.3755mm；边界层动量厚度为2.0819mm；取f=3()hz时，其风速近似为9.3m/s,经循环迭代得到，测

17、点距壁而的谋差距离为7.99mm；测虽:点摩擦速度为0.46m/s；计算出壁面的剪切应力为0.273n/m2；边界层厚度从上图可以近似得到，为11mm；边界层位移厚度为2.0217mm；边界层动量厚度为1.5098mm；基于动 量厚度的雷诺数为1122。3) 画出边界层速度分布曲线(基于边界层厚度和自山來流无量纲化)，边界层速度rms分 布illi线。画出边界层速度分布illi线(基于边界层内层参数无量纲化)。(图例为20hz) clc;clear;file='e:2014边界层试验x20f20 profile.txt'data=dlmread(file);jun k=sort

18、rows(data);yi=junk(:,l);ui=jun k(:,2);y=yi./12;u=ui./6.26;plot(u,y；-ob');xlabef基于自由来流风速无量纲化-bylabelc基于边界层厚度无量纲化')；axis(0.45 1.10 10)clc;clear;file='e:2014边界层试验x20f20 profile.txf; data=dlmread(file);jun k=sortrows(data);yi=junk(:j);rms=junk(:,3);y=yi-9.34;plot(rms,y,-ob,);xlabelfrms'h

19、ylabelc 边界层');axis(o 0.8 0 100)clc;clear;file='e:2014边界层试验x20f20 profile.txt'data=dlmread(file);jun k=sortrows(data);yi=junk(:,l)；ui=jun k(:,2);y=(yi-9.34)./12;u=ui./0.34;plot(u,y；-ob');xlabelc基于边界层摩擦速度无最纲化lylabef基于边界层厚度无最纲化axis(8 20 0 10)以上图形均为在20hz下的数据得到的，30hz的相应图形用同样的方法即可得到，在 此不再赘

20、述。4）对比分析讨论不同风洞运行速度下边界层速度分布曲线特性，边界层基本参数。10201000607080911 112圧干自由来流风連无蜗化1020020304rms05060708¥卜边連度f;和.墨干吕由来庚风速无列化0 102030405060708rms108 61014边2181项目20hz30hz自由來流速度m/s6.269.3测量点摩擦速度m/s0.340.46边界层厚度/mm1211边界层位移厚度/mm2.37552.0217边界层动量厚度/mm2.08191.5098边界层雷诺数10511122通过mat lab编程将两种频率下的三类图形放在一起进行比较，由数据文

21、件获知，在风 洞运行频率f=20hz时，自由來流速度约为6.26m/s,而在风洞运行频率f=30hz时，自由來 流速度为9.3m/s,从图上可以看到，无虽纲化后的图形儿乎一致，边界层速度分布曲线其基 本走势均为距离壁血较近的区域速度逐渐增加，到一定程度后就不再增加，这段距离即为边 界层厚度。5) 对比分析讨论不同风洞运行速度下剪切层速度大小，rms分布曲线特性。clc;clear;filel='e:2014边界层试验x20f20 profile.txt'file2='e:2014边界层试验x20f30 profile.txt'datal=dlmread(file

22、l);data2=dlmread(file2); junkl=sortrows(datal);ju nk2 二 sortrows(data2); yil=junkl(:,l);uil二junkl(:,2);yi2=junk2(:,l);ui2=junk2(:,2);yl=yil-9.34;ul=uil; y2=yi2-7.99;u2=ui2;figure(l);plot(ul,yl,-og,u2,y2；-*c'); legend('20hz730hz');xlabelcvlyiabelfy');titled不同频率下剪切层速度大小比较j;hold onyil=j

23、unkl(:,l);rmsl=junkl(:,3)； yi2=junk2(:,l)；rms2=junk2(:,3);yl=yil-9.34;y2=yi2-7.99;figure(2);plot(rmsl/yl；-og',rms2/y2/'-*c');legend('20hz730hz');xlabelcrms'bylabelcy');title('不同频率t rms大小比较从上面剪切层速度大小分布比较图和rms分布比较图两图可以看到,距壁面相同距离时,20hz下的剪切层速度小于30hz下的剪切层速度，从rms图上可以看到，基本趋势

24、两种风速下是一致的，但是距壁面相同距离时，30hz h的bps值要人于20hz时的值，而且30hz 时的rms最大值为1.18,而20hz时的最大值为0.704。6）在urmsm大值处，速度的能量谱图，分析能量谱特征。从图上可以看到峰值出现的地方，对应的频率即为该信号主要能量来源，从低到高频率 阶次依次为100hz, 160hz,220hz,280hz,而h.rfl图形趋势看到其主要能量集中于低频处，频 率越髙处能量越低。实验二 数据采集与处理1. 实验仪器1）信号发牛器及示波器:dsox2012a；2）数据采集卡：ni-usb 6212 bnc （16 位，8 差分通道，400 ks/s）；

25、3）计算机：控制和执行数据采集；2. 实验步骤1）连接信号发生器，数据采集卡，电脑等设备；2）设置信号发生器产生正弦波信号，信号频率为f=200hz,信号幅值为a=（）.（）5v；3）信号发生器产生信号连接到数据采集卡ai0；4）打开matlab数据采集程序acquisition.m,设置采样频率fs=100hz,通道输入范围('inputrange')为卜0.2 0.2。更改数据存储文件名,运行程序;5) 检查步骤4屮数据是否成功保持。重复步骤4改变采样频率fs=200hzo更改数据存储 文件名，运行程序；6) 检查步骤5中数据是否成功保持。重复步骤4改变采样频率fs=400

26、hzo更改数据存储文件名，运行程序；刀 检杳步骤6中数据是否成功保持。重复步骤4改变采样频率fs=1000hzo更改数据存储文件名，运行程序；8)检查步骤7中数据是否成功保持。重复步骤7,改变通道输入范围(tnpu(rang&)为卜1010oo更改数据存储文件名，运行程序。3. 实验数据处理1)画出数据随时间变化曲线(带点线)x(t)；不同参数下所得到的曲线画到同一张图上，用不同颜色标示。讨论区别，及采样参数影响。clc;clear;fipl=fopen(*e:2014数据采集试验 glyl00r02.bin7rf)；fip2=fopen(*e:2014数据采集试验 gly200r02

27、.bin','rt')；fip3=fopen(*e:2014数据采集试验 gly400r02.bin7rt')；fip4=fopen('e:2014数据采集试验 glyl000r02.bin'；rt')；fip5=fopen('e:2014数据采集试验 gly3000r02.bin','rt,);fip6=fopen(*e:2014数据采集试验 gly3000rl0.bin7rt');al/numl=fread(fipl/inf/,double')a2/num2=fread(fip2/inf/,do

28、uble')a3znum3=fread(fip3,inf/double')a4,num4=fread(fip4,inf,'double')a5,num5=fread(fip5,inf/double')a6/num6=fread(fip6/inf/,double')tl=(0:l:9)*10a-2;figure(l);plot(tlzal(l:10)；- or');% 第一hold ont2=(0:l:19)/2*10a-2;plot(t2za2(l:20)；-+m,);hold ont3=(0:l:39)/4*10a-2;plot(t3,

29、a3(l:40)；-*g,)；hold ont4=(0:l:99)*10a-3;plot(t4,a4(l:100),'>c)；hold ont5=(0:l:299)/3*10a-3;plot(t5,a5(l:300)；-<b');xlabel(时间 t/s1);ylabelfipm值 v/v*);title。不同采样频率k的曲线比较') figure(2);t6=(0:l:299)/3*10a-3;plot(t6,a5(l:300),yb);hold ont7=(0:l:299)/3*10a-3;plot(t7,a6(l:300)；-+k,);xlabel(

30、'时间 t/s');ylabel('幅值 v/v');title(不同釆样范围下的曲线比较jo o§8o幅值v/v§oooosog吕3o o49回湘求琢佛t3出隣長剂o§o人02«1bvaoc9m虽o>1可刃求落0st3陪帶一一芯从图上可以清址的看到：采样频率取loohz和200hz时完全失真，采到的为两条横线， 当取400hz （原始信号频率为200hz）时，采样信号有一点变化，但还离原始信号差很远, 因为不满足香农釆样定理（釆样频率400hz）,当把频率取到1000hz吋，釆样信号有了原 始信号的轮廓，但述是失真

31、比较大，当把频率取为3000hz时，采样信号己经很大程度的逼 近了原始信号，然后将采样幅值范围取为-10, 10,在峰值处有了平台信号失真，所以幅值 取为1-0.2,0.2非常接近真实信号。因此在采集数据时，必须满足采样定理，同时信号的灵敏 度要求幅值与真实信号幅值接近。2）计算采集到数据的平均值xm,均方差xrmso比较不同采样参数下得到值的区别。讨论 区别，及采样参数影响。clc;clear;fipl=fopen（'e:2014数据采集试验glyl00r02.bin,；rt,）；al,numl=fread(fipl/inf/,double1)； suml=0;for i=l:num

32、lsuml=al(i)+suml;endxl=suml/numl sum2=0;for i=l:numlcha=al(i)-xl;fang=cha.a2;sum2=sum2+fa ng;endxrmsl=sqrt(sum2/numl)依次算出6种情况下的均值和均方差(即标准差)。表2-1不同情况下采集数据的均值和均方差列表比较项目100hz0.2200hz0.2400hz0.21000hz0.23000hz0.23000hz10平均0.0110-0.0168-3.7972e-4-3.5217e-4-3.5032e-4-0.0010值均方0.00209.5665e40.00320.01770.0

33、1770.0177差编程计算求得各自的均值和均方差。从上表可以看出，100hz和200fk时均方差较小， 而均值较大，这与实际情况一致，因为具图形近乎为一直线；400hz时严重失真，图形呈锯 齿波形，均值较小近乎为0；而后面三种高频采集的信号均值都很小，而h均方差值接近, 较好地反映了原始信号，相比z下低的采样频率采集的信号均方差的误差较大。3）用mat lab编程计算所得数据的傅里叶变化x（k）,验证parsovals theorenio画出|x（k） | 随频率变化曲线。不同采样参数得到的曲线画到同一张图上，用不同颜色标示。讨论区 另u，及采样参数影响。clc;clear;fipl=fop

34、en（'e:2014数据釆集试验glylooro.2.bin','rt'）；fip2=fopen（'e:2014数据采集试验gly200r0.2.bin：'卍）;fip3=fopen（'e:2014 数据采集试验 gly400r0.2.bin7rt'）;fip4=fopen（'e:2014 数据采集试验 glyl000r0.2.bin7rt'）;fip5=fopen'e:2014数据采集试验 xglxylooooro.bin'/rt'）;fip6=fopen（'e:2014数据釆集试

35、验glyloooorlo.bin','rt'）;datal=fread（fipl,512/double'）data2=fread（fip2/512/,double1）data3=fread(fip3,512/double') data4=fread(fip4z512,'double') data5=fread(fip5z512/double') data6=fread(fip6,512/double') n=512; %512个采样点 figure ；fs=loo; %采样频率t=(o:n-l)/fs; %时域y=fft(

36、datal,n); %快速傅里叶变换 f=(o:n-l)*fs/n; %频域的点数 plot(f/abs(y)；-sr,) % 频谱图 hold onfs=200;t=(o:n-l)/fs;y=fft(data2?n);f=(o:n-l)*fs/n;plot(f,abs(y)/-dm')hold onfs=400;t=(o:n-l)/fs;y=fft(data3,n);f=(o:n-l)*fs/n;plot(f,abs(y)；-hg')hold onfs=looo;t=(o:n-l)/fs;y=fft(data4,n); f=(o:n-l)*fs/n;plot(f,abs(y)

37、/-*c')hold onfs=loooo;t=(o:n-l)/fs;y=fft(data5,n);f=(o:n-l)*fs/n; plot(fzabs(y)；-oy') hold onfs=loooo;t=(o:n-l)/fs;y=fft(data6,n); f=(o:n-l)*fs/n;plot(f,abs(y)/-pk')204060801004i3<>2105010015020061«'«42i份0 - -'02000400060008000100000200040006000800010000从上图6种采样信号单

38、独的傅立叶变换后的图形，可以很明显地看到，loohz和200hz 两种情况下，只有（）hz处有幅值，属于肓流分量，而其余部分近乎为零，没有能蜃，没有 采集到信号，或者说采集到的信号完全失真。一般来说，满足采样定理的信号频谱图是对称的，后面的儿种情况就是这样，明显看到在200hz处出现峰值，表示能量集屮在此频率下， 说明原始信号的频率为200hz,这与实际情况是相符的。在后两种情况f,采样频率是相同 的(10000hz),但是采样范围是不同的(卜0.2,0.2,-10,10),而从图上看到其傅立叶变换的 频谱图是完全相同的，表明频谱图与采样频率密切相关的。4) 画岀能量谱图，不同采样参数得到的曲线画到同一张图上，用