,280 北京科技大学学报 第32卷 种规格（通过在地板上间隔20am设一凹槽，根据要 求面积不同，壁板安装在不同的凹槽内)高度方向 分别设0.20.30.40.5和0.6m五种高度测试方 n 案（通过在壁板上每间隔10am设凹槽，根据高度要 求，将顶板安置在相应高度凹槽内，通过吊绳与壁板 9 凹槽结合施加约束)，每种测试方案只需要将壁板 和顶板调整到相应位置，相关壁板位置调试好后，进 行密封，尽可能降低漏风，调整好顶板的平整度，用 50450 00 0.5 钢丝绳固定，当测试条件准备好后，剪断吊绳，则1 50 0.4 e300 500.2 0.3m 板首先掉落，2二、3板失去约束相继掉落，然后是4、 图3不盖板数据二次插值后的图像 5板掉落，风速测试采用数字式风表量取最大值 Fig3 mage of uncovered boad data after inserting quadmatic vahes 根据理论公式(12)中风速与高度、面积比呈幂 函数的关系，构造幂函数： v=ah"s. 取对数得线性函数： lgv=ga十blgh十cgs 对表1中不盖板部分实验数据利用二元线性回 归求得参数估计： lga=0.3536=a=2.2574 b=0.3447,c=0.3058 图2风速测定实验装置 所以幂函数为： Fig 2 Experinental device for esting the w ind speed v=2.2574h0.370-3058 (13) 4实验分析 4.2类绕流模型分析（盖板） 4.1类打气筒模型分析（不盖板） 对表1中盖板部分实验数据利用MATLAB作 对表1中不盖板部分实验数据利用MATLAB 二次插值后图像如图4所示 作二次插值后图像如图3所示 表1冲击气浪风速测定实验数据 11 Tab le 1 Experment data of w ind speed of air shock waves 10 m.s-l 9 高度/面积比125面积比250面积比375面积比500 m 盖板不盖板盖板不盖板盖板不盖板盖板不盖板 6 6.18.54.26.46.77.98.48.5 0.26.47.84.56.86.37.78.28.3 0.5 5.88.65.47.26.68.18.48.6 0.4 8.39.35.87.77.28.48.710.4 面积比 30075062 0.37.99.45.97.97.18.58.310.2 图4盖板数据二次插值后的图像 7.58.75.68.07.68.38.69.7 Fig 4 mage of covered boand data after inserting quadratic values 8.611.06.58.68.19.29.59.6 0.48.111.37.78.78.49.48.710.2 构造幂函数： 8.311.87.78.77.99.19.311.3 v-ah's. 9.813.38.89.49.010.39.711.9 取对数得线性函数： 0.58.212.58.210.19.311.110.211.7 lgv=ga十blgh十clgs 9.312.68.29.68.910.79.811.1 6.28.18.39.99.911.311.312.1 对表1中盖板部分实验数据利用二元线性回归 0.66.38.38.510.210.111.611.412.7 求得参数估计： 6.99.78.410.19.710.911.112.3 lga=-0.0086=a=0.9804北 京 科 技 大 学 学 报 第 32卷 种规格 （通过在地板上间隔 20ｃｍ设一凹槽‚根据要 求面积不同‚壁板安装在不同的凹槽内 ）‚高度方向 分别设 0∙2、0∙3、0∙4、0∙5和 0∙6ｍ五种高度测试方 案 （通过在壁板上每间隔 10ｃｍ设凹槽‚根据高度要 求‚将顶板安置在相应高度凹槽内‚通过吊绳与壁板 凹槽结合施加约束 ）．每种测试方案只需要将壁板 和顶板调整到相应位置‚相关壁板位置调试好后‚进 行密封‚尽可能降低漏风‚调整好顶板的平整度‚用 钢丝绳固定．当测试条件准备好后‚剪断吊绳‚则 1 ＃ 板首先掉落‚2 ＃、3 ＃板失去约束相继掉落‚然后是 4 ＃、 5 ＃板掉落．风速测试采用数字式风表量取最大值． 图 2 风速测定实验装置 Ｆｉｇ．2 Ｅｘｐｅｒｉｍｅｎｔａｌｄｅｖｉｃｅｆｏｒｔｅｓｔｉｎｇｔｈｅｗｉｎｄｓｐｅｅｄ 4 实验分析 4∙1 类 “打气筒 ”模型分析 （不盖板 ） 对表 1中不盖板部分实验数据利用 ＭＡＴＬＡＢ 作二次插值后图像如图 3所示． 表 1 冲击气浪风速测定实验数据 Ｔａｂｌｅ1 Ｅｘｐｅｒｉｍｅｎｔｄａｔａｏｆｗｉｎｄｓｐｅｅｄｏｆａｉｒｓｈｏｃｋｗａｖｅｓ ｍ·ｓ—1 高度／ ｍ 面积比 125 面积比 250 面积比 375 面积比 500 盖板 不盖板 盖板 不盖板 盖板 不盖板 盖板 不盖板 6∙1 8∙5 4∙2 6∙4 6∙7 7∙9 8∙4 8∙5 0∙2 6∙4 7∙8 4∙5 6∙8 6∙3 7∙7 8∙2 8∙3 5∙8 8∙6 5∙4 7∙2 6∙6 8∙1 8∙4 8∙6 8∙3 9∙3 5∙8 7∙7 7∙2 8∙4 8∙7 10∙4 0∙3 7∙9 9∙4 5∙9 7∙9 7∙1 8∙5 8∙3 10∙2 7∙5 8∙7 5∙6 8∙0 7∙6 8∙3 8∙6 9∙7 8∙6 11∙0 6∙5 8∙6 8∙1 9∙2 9∙5 9∙6 0∙4 8∙1 11∙3 7∙7 8∙7 8∙4 9∙4 8∙7 10∙2 8∙3 11∙8 7∙7 8∙7 7∙9 9∙1 9∙3 11∙3 9∙8 13∙3 8∙8 9∙4 9∙0 10∙3 9∙7 11∙9 0∙5 8∙2 12∙5 8∙2 10∙1 9∙3 11∙1 10∙2 11∙7 9∙3 12∙6 8∙2 9∙6 8∙9 10∙7 9∙8 11∙1 6∙2 8∙1 8∙3 9∙9 9∙9 11∙3 11∙3 12∙1 0∙6 6∙3 8∙3 8∙5 10∙2 10∙1 11∙6 11∙4 12∙7 6∙9 9∙7 8∙4 10∙1 9∙7 10∙9 11∙1 12∙3 图 3 不盖板数据二次插值后的图像 Ｆｉｇ．3 Ｉｍａｇｅｏｆｕｎｃｏｖｅｒｅｄｂｏａｒｄｄａｔａａｆｔｅｒｉｎｓｅｒｔｉｎｇｑｕａｄｒａｔｉｃｖａｌｕｅｓ 根据理论公式 （12）中风速与高度、面积比呈幂 函数的关系‚构造幂函数： ｖ＝ａｈ ｂｓ ｃ． 取对数得线性函数： ｌｇｖ＝ｌｇａ＋ｂｌｇｈ＋ｃｌｇｓ． 对表 1中不盖板部分实验数据利用二元线性回 归求得参数估计： ｌｇａ＝0∙3536⇒ａ＝2∙2574‚ ｂ＝0∙3447‚ｃ＝0∙3058． 所以幂函数为： ｖ＝2∙2574ｈ 0∙3447ｓ 0∙3058 （13） 4∙2 类 “绕流 ”模型分析 （盖板 ） 对表 1中盖板部分实验数据利用 ＭＡＴＬＡＢ作 二次插值后图像如图 4所示． 图 4 盖板数据二次插值后的图像 Ｆｉｇ．4 Ｉｍａｇｅｏｆｃｏｖｅｒｅｄｂｏａｒｄｄａｔａａｆｔｅｒｉｎｓｅｒｔｉｎｇｑｕａｄｒａｔｉｃｖａｌｕｅｓ 构造幂函数： ｖ＝ａｈ ｂｓ ｃ． 取对数得线性函数： ｌｇｖ＝ｌｇａ＋ｂｌｇｈ＋ｃｌｇｓ． 对表 1中盖板部分实验数据利用二元线性回归 求得参数估计： ｌｇａ＝—0∙0086⇒ａ＝0∙9804‚ ·280·