1、实 验 流 体 力 学 风洞静态压力分布测量实验实验设计及数据处理程序编制数据采集与分析姓 名 学 号 实验日期 指导老师 1一、实验目的风洞测压试验是一种在风洞中测量模型表面压力分布的试验。目的是通过测量飞行器及其部件，如机翼、机身、尾翼、操纵面、外挂物等表面上的压力分布，为飞行器及其部件结构强度计算提供气动载荷分布的原始数据。通过测压实验，能够给出定量化的结果，获得模型上的压力分布信息。因此，测压试验是研究模型气动特性、验证数值计算方法的一种重要手段。本次实验内容是测定标准模型在不同实验状态下各截面测压点的压力值，并进行数据处理，最后得到各截面的压力分布曲线随风速及迎角的变化规律。 二实验
2、设备1、风洞风洞是产生人工气流的设备，本次实验所用风洞为开口回流式低速风洞，如图 1 所示。其主要组成部分为实验段、扩压段、拐角和导流片、稳定段、收缩段以及动力段。 D4 风 洞 实 验 段风 扇图 1 D4 风洞示意图实验段尺寸：长度3.5m，宽度 1.5m，高度1.5m，收缩比9；实验段风速：闭口最高风速为80m / s,开口最高风速为60m / s；实验段湍流度为0.08% 。2、风速控制系统D4 风洞采用可控硅控制无级调速； 风速控制系统组成如图 1 所示。本文的工作是在系统的外层增加了稳风速的闭环控制系统。风洞风速的控制采用直接数字式闭环控制。首先通过 PCL727 进行 D/A 转
3、换，将数字量转换成模拟量 420mA 驱动电流，经过西门子驱动器来控制可控硅的输出电流，从而控制电机转速，电机拖动风扇，产生气流，使试验段获得所需的速压 (其中 为气0P0p体总压， 为气体静压)。由差压变送器将压差转换成 1-5V 电压，再由压差风速转换公式计算得出p风速值，经 PCI1716 进行 A/D 转换，将数字量输入计算机，通过数字 PID 控制器输出控制量，从而改变输出的驱动电流，达到控制风速的目的。三、实验模型：1、实验模型：非圆截面机身标准测压模型，如图2所示。2图2 实验模型2、实验条件：1、固定模型迎角= 30 ；调节实验风速： 10 米/秒- 25 米/ 秒 2、固定实
4、验风速V= 25 米/秒；调节模型迎角： = 10 - 503、测压截面测压点分布及坐标定义为绘制压力系数沿截面周向的分布曲线，给出压力分布曲线 Cp=f(y/b)中横坐标的定义，如图3 所示。其中，b=D/2，y 为 Y 轴坐标。图 3 模型测压点分布及压力分布曲线横坐标定义为方便测压实验的进行，首先将测压点进行了编号，如图 4 所示，按照顺时针方向从 1 到22。实验中将模型上的 22 个测压点，通过测压软管，按照编号顺序连接到 DTC 压力采集模块的 1到 22 通道上。另外，DTC 压力采集模块的 23 通道测量实验中自由来流的总压，24 通道测量自由来流的静压。这样就可以通过测得的压
5、力数据计算出模型表面的压力系数。图 4 模型测压点编号3四、测压实验数据处理数据处理就是将数据转换为信息的过程。其基本目的是从大量的数据中根据数据自身的规律和及其相互联系，抽取并推导出对于某些特定的人们来说是有价值、有意义的数据，为进一步的分析提供依据。1、实验数据说明采集得到的原始数据文件以 .txt 形式保存。数据第 1 列为采集时间，第 2 列到第 25 列为传感器采集到的压力，其中，第 2 列到第 23 列对应测压点从 1 到 22 的压力，第 24 列为总压，第 25 列为静压。采集得到的压力为相对于大气压的差压，单位为 Pa。由于测量流场为定常流场，实验中对所有测点进行了多次测量，
6、数据处理时可通过平均得到每个测点的测量值。2、根据实验数据，编写程序，计算每个测压点测得的压力系数，并绘制出模型机身截面在不同实验状态下的 Cp=f(y/b)曲线，如下图 5 所示。压力系数要通过测得的测压点压力及总压、静压值计算。并对不同实验状态下的结果进行对比分析。公式： 21/iiU图 5 绘制完成的压力系数沿模型周向的分布曲线3、根据模型截面的 Cp=f(y/b)曲线，编写程序，积分该截面的截面法向力系数 CN 和截面侧向力系数 CY。并对不同实验状态下的结果进行对比分析。计算截面力系数时，截面轴向宽度取单位宽度 1mm，则每个测压点对应的面积为 ds1 mm2。积分采用线性叠加的方式
7、。C N 的参考面积为 D1 mm2，C Y 的参考面积为 H1 mm2。图 6 积分过程中相关量示意注：为计算压力系数和截面法向力、侧向力系数，所需相关数据见表 1。4表 1 测压点相关参数测压点编号 1 2 3 4 5 6 7 8Y 坐标/mm -15 -30 -45 -60 -60 -50 -40 -30Z 坐标 /mm -66.75 -57.17 -41.07 -18.57 13.93 25.714 35.359 42.882法向 Z 轴分量 0.921 0.758 0.614 0.504 -0.614 -0.682 -0.758 -0.841法向 Y 轴分量 -0.39 -0.652
8、 -0.789 -0.864 -0.789 -0.731 -0.652 -0.541积分长度 ds/mm 9.52 11.403 14.03 19.787 14.218 8.389 7.55 6.819测压点编号 9 10 11 12 13 14 15 16Y 坐标/mm -20 -10 0 10 20 30 40 50Z 坐标 /mm 48.203 51.396 52.5 51.396 48.203 42.882 35.359 25.714法向 Z 轴分量 -0.921 -0.977 -1 -0.977 -0.921 -0.841 -0.758 -0.682法向 Y 轴分量 -0.39 -0
9、.212 0 0.212 0.39 0.541 0.652 0.731积分长度 ds/mm 6.244 5.884 5.76 5.884 6.244 6.819 7.55 8.389测压点编号 17 18 19 20 21 22 Y 坐标/mm 60 60 45 30 15 0 Z 坐标 /mm 13.93 -18.57 -41.07 -57.17 -66.75 -70 法向 Z 轴分量 -0.614 0.504 0.614 0.758 0.921 1 法向 Y 轴分量 0.789 0.864 0.789 0.652 0.39 0 积分长度 ds/mm 14.218 19.787 14.03
10、11.403 9.52 8.829 4、根据实验所得数据编写程序，得到各个实验条件下得 Cp 值。 （程序见附件 1，计算结果见附件2-9）根据数据处理结果，作出各个条件下的 Cp-y/b 图，如下所示：图 7 v10-a305图 8 v15-a30图 9 v20-a306图 10 v25-a30图 11 v25-a107图 12 v25-a20图 13 v25-a408图 14 v25-a50从图中可见：气流经过该机身截面的流动图画，有类似圆柱扰流的特征。最大正压点为 22 点；最大负压点为 8、14 点附近，该点为速度最大点。通过图 710 可以看出，当攻角确定为 30速度从 10 m/s
11、25 m/s变化时，有如下规律：随着速度的增加，流动现象基本不变，特征点位置基本没有变化。只有最大负压增加，升力增加通过图 1013 可以看出，当速度确定为 25 m/s，攻角从 1050变化时，有如下规律：攻角为 10时，截面上全部都为负压区随着攻角增加，从截面尾部逐渐出现正压区5、根据模型截面的 Cp=f(y/b)曲线，编写程序，积分得到了该截面的截面法向力系数 CN 和截面侧向力系数 CY，将结果列表如下（程序见附录 11）v a CN CY10 30 .26841 .0049115 30 .27535 .0104120 30 .28086 .0096725 30 .28539 .0074425 10 .02728 -.0006825 20 .11880 -.0016425 40 .39287 .0058425 50 .47948 -.00242攻角不变时，随着速度的增加，CN 即升力系数增加，CY 是个小量，先增加后减小。9速度不变时，随着攻角的增加，CN 增加，CY 的变化规律不太明显。通过比较可知，对于机身模型，升力随着速度和攻角的增加而增大，失速攻角很大，超过 50。
