风洞静态压力分布测量实验

实实验验流流体体力力学学 风洞静态压力分布测量实验风洞静态压力分布测量实验 实验设计及数据处理程序编制 数据采集与分析 姓姓 名名 学学 号号 实验日期实验日期 指导老师指导老师 1 一 实验目的 风洞测压试验是一种在风洞中测量模型表面压力分布的试验 目的是通过测量飞行器及其部件 如机翼 机身 尾翼 操纵面 外挂物等表面上的压力分布 为飞行器及其部件结构强度计算提供 气动载荷分布的原始数据 通过测压实验 能够给出定量化的结果 获得模型上的压力分布信息 因此 测压试验是研究模型气动特性 验证数值计算方法的一种重要手段 本次实验内容是测定标准模型在不同实验状态下各截面测压点的压力值 并进行数据处理 最 后得到各截面的压力分布曲线随风速及迎角的变化规律 二 实验设备 1 风洞 风洞是产生人工气流的设备 本次实验所用风洞为开口回流式低速风洞 如图 1 所示 其主要 组成部分为实验段 扩压段 拐角和导流片 稳定段 收缩段以及动力段 D4 风洞 实验段 风 扇 图 1 D4 风洞示意图 实验段尺寸 长度3 5m 宽度1 5m 高度1 5m 收缩比9 实验段风速 闭口最高风速为80m s 开口最高风速为60m s 实验段湍流度为0 08 2 风速控制系统 D4 风洞采用可控硅控制无级调速 风速控制系统组成如图 1 所示 本文的工作是在系统的外 层增加了稳风速的闭环控制系统 风洞风速的控制采用直接数字式闭环控制 首先通过 PCL727 进 行 D A 转换 将数字量转换成模拟量 4 20mA 驱动电流 经过西门子驱动器来控制可控硅的输出电 流 从而控制电机转速 电机拖动风扇 产生气流 使试验段获得所需的速压 其中为气 0 PP 0 p 体总压 为气体静压 由差压变送器将压差转换成 1 5V 电压 再由压差风速转换公式计算得出p 风速值 经 PCI1716 进行 A D 转换 将数字量输入计算机 通过数字 PID 控制器输出控制量 从 而改变输出的驱动电流 达到控制风速的目的 三 实验模型 1 实验模型 非圆截面机身标准测压模型 如图2所示 2 图2 实验模型 2 实验条件 1 固定模型迎角 30 调节实验风速 10 米 秒 25 米 秒 2 固定实验风速V 25 米 秒 调节模型迎角 10 50 3 测压截面测压点分布及坐标定义 为绘制压力系数沿截面周向的分布曲线 给出压力分布曲线 Cp f y b 中横坐标的定义 如图 3 所示 其中 b D 2 y 为 Y 轴坐标 图 3 模型测压点分布及压力分布曲线横坐标定义 为方便测压实验的进行 首先将测压点进行了编号 如图 4 所示 按照顺时针方向从 1 到 22 实验中将模型上的 22 个测压点 通过测压软管 按照编号顺序连接到 DTC 压力采集模块的 1 到 22 通道上 另外 DTC 压力采集模块的 23 通道测量实验中自由来流的总压 24 通道测量自由 来流的静压 这样就可以通过测得的压力数据计算出模型表面的压力系数 图 4 模型测压点编号 3 四 四 测压实验数据处理 数据处理就是将数据转换为信息的过程 其基本目的是从大量的数据中根据数据自身的规律和 及其相互联系 抽取并推导出对于某些特定的人们来说是有价值 有意义的数据 为进一步的分析 提供依据 1 实验数据说明 采集得到的原始数据文件以 txt 形式保存 数据第 1 列为采集时间 第 2 列到第 25 列为传感 器采集到的压力 其中 第 2 列到第 23 列对应测压点从 1 到 22 的压力 第 24 列为总压 第 25 列 为静压 采集得到的压力为相对于大气压的差压 单位为 Pa 由于测量流场为定常流场 实验中 对所有测点进行了多次测量 数据处理时可通过平均得到每个测点的测量值 2 根据实验数据 编写程序 计算每个测压点测得的压力系数 并绘制出模型机身截面在不同实 验状态下的 Cp f y b 曲线 如下图 5 所示 压力系数要通过测得的测压点压力及总压 静压值计 算 并对不同实验状态下的结果进行对比分析 公式 2 1 2 i i pp Cp U 图 5 绘制完成的压力系数沿模型周向的分布曲线 3 根据模型截面的 Cp f y b 曲线 编写程序 积分该截面的截面法向力系数CN和截面侧向力系 数CY 并对不同实验状态下的结果进行对比分析 计算截面力系数时 截面轴向宽度取单位宽度 1mm 则每个测压点对应的面积为 ds 1 mm2 积分采用线性叠加的方式 CN的参考面积为 D 1 mm2 CY的参考面积为 H 1 mm2 图 6 积分过程中相关量示意 注 为计算压力系数和截面法向力 侧向力系数 所需相关数据见表 1 4 表 1 测压点相关参数 测压点编号12345678 Y 坐标 mm 15 30 45 60 60 50 40 30 Z 坐标 mm 66 75 57 17 41 07 18 5713 9325 71435 35942 882 法向 Z 轴分量0 9210 7580 6140 504 0 614 0 682 0 758 0 841 法向 Y 轴分量 0 39 0 652 0 789 0 864 0 789 0 731 0 652 0 541 积分长度 ds mm9 5211 40314 0319 78714 2188 3897 556 819 测压点编号910111213141516 Y 坐标 mm 20 1001020304050 Z 坐标 mm48 20351 39652 551 39648 20342 88235 35925 714 法向 Z 轴分量 0 921 0 977 1 0 977 0 921 0 841 0 758 0 682 法向 Y 轴分量 0 39 0 21200 2120 390 5410 6520 731 积分长度 ds mm6 2445 8845 765 8846 2446 8197 558 389 测压点编号171819202122 Y 坐标 mm60604530150 Z 坐标 mm13 93 18 57 41 07 57 17 66 75 70 法向 Z 轴分量 0 6140 5040 6140 7580 9211 法向 Y 轴分量0 7890 8640 7890 6520 390 积分长度 ds mm14 21819 78714 0311 4039 528 829 4 根据实验所得数据编写程序 得到各个实验条件下得 Cp 值 程序见附件 1 计算结果见附件 2 9 根据数据处理结果 作出各个条件下的 Cp y b 图 如下所示 图 7 v10 a30 5 图 8 v15 a30 图 9 v20 a30 6 图 10 v25 a30 图 11 v25 a10 7 图 12 v25 a20 图 13 v25 a40 8 图 14 v25 a50 从图中可见 气流经过该机身截面的流动图画 有类似圆柱扰流的特征 最大正压点为 22 点 最 大负压点为 8 14 点附近 该点为速度最大点 通过图 7 10 可以看出 当攻角确定为 30 速度从 10 m s 25 m s 变化时 有如下规律 随着速度的增加 流动现象基本不变 特征点位置基本没有变化 只有最大负压增加 升力增加 通过图 10 13 可以看出 当速度确定为 25 m s 攻角从 10 50 变化时 有如下规律 攻角为 10 时 截面上全部都为负压区 随着攻角增加 从截面尾部逐渐出现正压区 5 根据模型截面的 Cp f y b 曲线 编写程序 积分得到了该截面的截面法向力系数 CN和截面侧 向力系数 CY 将结果列表如下 程序见附录 11 v a CN CY 10 30 26841 00491 15 30 27535 01041 20 30 28086 00967 25 30 28539 00744 25 10 02728 00068