Measurement of CO concentration in flat flame based on mid-infrared absorption spectroscopy
-
摘要: CO是碳氢化合物燃烧的主要产物之一，准确测量超燃冲压发动机出口的CO浓度是评估碳氢燃料燃烧效率的重要依据。中红外波段的CO谱线相较近红外而言，具有吸收更强、谱线丰富且谱线对相对孤立、不受其他气体干扰等明显优势。本文基于中红外吸收光谱技术，计算研究了CO中红外光谱特性，选择了适用于高温流场CO测量的特征谱线，设计并搭建了高温流场CO浓度检测系统，开展了气体池浓度标定和不同当量比下平面火焰CO测量验证，实现了某超燃冲压发动机出口高温流场CO测量，反映了航空煤油燃烧过程中CO浓度和温度的变化情况，为超燃冲压发动机的燃烧和流动机理研究提供了有力的研究手段和丰富的实验数据。Abstract: CO is one of the main products of hydrocarbon combustion. Accurate measurement of CO concentration at the scramjet outlet is an important basis for evaluating the combustion efficiency of the hydrocarbon fuel. Compared with the near infrared band, the CO absorption spectrum in the mid infrared band has the advantages of stronger absorption, rich spectral lines, relatively isolated spectral lines and no interference from other gases. Based on the mid infrared absorption spectrum technology, this paper studies the mid infrared spectrum characteristics of CO, selects the characteristic spectrum which is suitable for the measurement of CO in the high temperature flow field, designs and builds the CO concentration detection system in the high temperature flow field, carries out the CO measurement verification under different equivalence ratios of the plane flame, realizes the CO measurement of the high temperature flow field at the exit of a scramjet, and reflects the changes of concentration and temperature of CO during the combustion of aviation kerosene, thus providing powerful research means and abundant experimental data for the study of combustion and flow mechanism of the scramjet.
-
Keywords:
- mid-infrared absorption spectroscopy /
- TDLAS /
- CO /
- flat flame /
- scramjet
-
0 引言
翼型是飞机机翼、尾翼等空气动力部件的核心要素, 在研制新翼型的过程中必须采用风洞实验验证其各项气动力特性。翼型的升力和阻力是其气动力特性至关重要的气动参数, 一般来说, 翼型阻力比升力小一个数量级, 精确测量阻力相对更加困难, 发展实用可靠的翼型阻力精细测量技术具有十分重要的工程应用价值。
在风洞实验中, 一般采用尾迹积分法测量翼型阻力, 即在翼型下游某截面安装尾耙, 通过测量该截面的总压、静压获得该截面的动量损失, 再对动量损失进行积分获得翼型阻力。尾耙的高度和安装位置等因素对翼型阻力测量存在一定影响:(1)尾耙的安装位置会直接影响阻力测量的精度。如果尾耙安装位置离翼型太近, 尾耙会对翼型绕流产生干扰, 改变翼型压力分布, 而且尾迹区动量损失曲线形状会过于"尖瘦", 对其进行积分计算阻力时会带来较大误差; 如果尾耙安装位置离翼型过远, 虽然动量损失曲线足够"肥胖"光滑, 但其一部分可能会超出尾耙测点区域, 也会给积分带来较大误差, 影响阻力测量精度。对于低速风洞, 尾耙一般安装于翼型下游0.5~1倍弦长处; 而对于高速风洞, 由于尾耙对翼型的扰动增大, 位置一般更靠后一些。(2)尾耙的高度以能捕获整个翼型尾迹为原则, 一般都大于1倍翼型弦长。对于翼型大迎角实验和高速实验, 其尾迹区的宽度较大, 尾耙高度宜选择大一些(可能需3倍弦长以上)。而在满足强度和刚度的条件下, 尾耙的其他尺寸应尽量小一些, 以减少对流场的干扰。
Taylor[1]、Bollay[2]和Goett[3]等最早研究发现尾迹截面的位置对阻力测量的影响, 当尾迹截面距离翼型后缘过近会导致计算误差较大; Vinh[4]通过数值模拟研究在大气中采用尾迹积分法计算单个翼型阻力时发现, 一般在翼型下游4~5倍弦长以后, 尾迹积分的位置对阻力计算结果几乎无影响。周瑞兴等[5]进行了尾迹法测量阻力的研究, 设计了头部15°倒角的总压探头, 保证了在流线有很大偏角时仍能测得比较准确的总压值; 焦予秦、支真莉[6-7]改进了低速风洞中以尾迹动量损失积分法测量翼型阻力的积分公式, 加入了压差影响和气流偏角影响, 得到了更好的高升力构型的阻力测量结果。以上工作均是针对低速实验进行研究, 对于跨声速实验, 周伟等[8]以数值模拟方法研究了风洞中尾耙的轴向安装位置对阻力测量的影响, 得出了"尾耙安装于翼型下游2倍弦长之后的测量值比较可靠"的结论。
在西北工业大学NF-6增压连续式跨声速翼型风洞中, 大多采用安装于翼型下游1.35倍弦长处的尾耙测量阻力。研究结果表明, 此处尾耙离翼型较近, 影响阻力测量结果。在文献的基础上, 根据风洞实际结构, 在翼型下游2.90倍弦长处的风洞侧壁上加工槽缝、安装尾耙进行实验, 对比分析了离模型较近及离模型足够远的尾耙阻力测量结果。另外, 对常用尾耙的多次实验结果分析发现, 在较大迎角或较大马赫数下, 尾耙高度不足以完全覆盖动量损失范围。需基于测得的动量损失曲线得到合适的尾耙高度, 且考虑尾耙位置后移需要更大的尾耙高度, 设计新的尾耙进行实验。
本文采用2个不同高度的尾耙, 在西北工业大学NF-6跨声速风洞实验段的2个轴向安装位置进行同一翼型实验, 通过对比分析同一尾耙高度、不同安装位置以及同一安装位置、不同尾耙高度的阻力实验数据, 获得不同实验条件下合适的尾耙高度及安装位置的组合, 改进NF-6风洞阻力测量精度。
1 实验风洞和模型
实验在西北工业大学NF-6增压连续式跨声速翼型风洞二元实验段进行。实验段尺寸为0.8 m×0.4 m×3 m (高×宽×长), 上下壁为开孔壁, 左右壁为实壁[9]。实验模型为OA309旋翼翼型(弦长c=200 mm), 展长为400 mm。图 1为翼型和尾耙在风洞中的安装照片。
2 阻力测量方法和尾耙
实验中, 以尾耙测量尾迹横截面的动量损失。尾耙上布置总压探针和静压探针。总压探针测量其前端所在的风洞横截面位置的总压(一般总压沿垂直于上下壁面的方向是变化的); 静压探针数量较少, 一般取其测得的平均值作为该横截面上的静压值。
实验中使用2种不同高度的尾耙。A尾耙竖向高度为270 mm(总压探针测量范围为260 mm), B尾耙竖向高度为500 mm(总压探针测量范围为480 mm); A、B尾耙的其他外形尺寸一致。尾耙的2个轴向安装位置分别为1位置和2位置。尾耙安装在1位置时, 总压探针前端距翼型前缘距离为470 mm, 即总压探针前端距离翼型后缘1.35 c(弦长c＝200 mm); 尾耙安装在2位置时, 总压探针前端距离翼型后缘2.90 c。
A尾耙的几何数据如图 2所示。尾耙上共有53根总压探针, 编号为1~53, 上下两端探针分布较疏、中间较密, 探针之间的距离为5 mm×17+4 mm×20+5 mm×17;4根静压探针均布, 编号为54~57。B尾耙探针分布和A尾耙类似:总压探针共103根(5 mm×32+4 mm×40+5 mm×32); 4根静压探针均布。
如图 3所示, 取控制体0截面-Ⅰ截面。0截面在翼型远前方, 此处气流未受翼型扰动, 来流静压和总压分别为p∞、p0; Ⅰ截面在翼型下游无限远, 此处气流静压已恢复至来流静压p∞, 总压为p0, 1, 根据动量定理可得翼型阻力Q[10]:
![]()
图 3 风洞中动量法测翼型阻力原理图Fig. 3 Principles for drag measurement with momentum method in wind tunnel
(1) 式中, ρ为气体密度, v∞为来流速度, v1为Ⅰ截面气流速度, y为纵坐标方向。
由于Ⅰ截面在翼型下游无限远处, 尾迹范围无限大, 不可能在该处安装尾耙测量阻力, 因此, 在翼型下游适当位置取Ⅱ截面安装尾耙。假设尾耙测量截面尾迹内的静压已保持为常数, 并且气流在Ⅰ截面-Ⅱ截面之间没有能量损失, 则尾耙测量截面尾迹内的总压也为p0, 1; 而在尾迹之外, 认为气流没有能量损失, 因此总压p0保持不变。定义Ⅱ截面尾耙处气流总压、静压分别为p0, w、pw。
对于高速风洞实验, 风洞中各截面的气流密度并不相同。在Ⅱ截面处, 对式(1)进行推导可得翼型阻力系数CD为[6, 10-12]:
(2) 式中, CD, wt为空风洞时风洞本身流动损失引起的阻力系数, h为尾耙高度, lm为Ⅱ截面处的无量纲单位高度上的动量损失(Momentum loss), 即:
(3) 式中, γ为气体比热比(空气中γ＝1.4)。
3 实验结果分析
3.1 尾耙不同安装位置对阻力测量的影响
将A尾耙分别安装在1和2位置进行实验。图 4给出了马赫数Ma＝0.50和0.80条件下A尾耙在1、2位置的阻力系数。可以看出:迎角α较小时, A尾耙在1、2位置的阻力系数差别不大; 随着迎角α增大(无论是较大的负迎角还是正迎角), 1位置的阻力系数大于2位置的阻力系数。
![]()
图 4 A尾耙在1、2位置处的阻力系数Fig. 4 Drag coefficient for wake rake A at installation position 1 and 2在Ma=0.50和0.80条件下, A尾耙在1、2位置的动量损失曲线如图 5和6所示。从图 5可以看出, 当迎角大于±10°, 2位置的尾迹动量损失范围超出了尾耙测量范围, 尾耙高度明显不够, 这样积分得到的阻力系数是不真实的; 从图 6可以看出, 2位置的尾迹动量损失范围超出尾耙测量范围的正迎角提前到9°, 负迎角提前到-8°, 说明随着马赫数增大, 同一迎角状态下尾迹动量损失范围扩大, 所需的尾耙高度更大。A尾耙安装在2位置, 其高度不足以覆盖尾迹范围, 测得的阻力系数值小于真实值, 这是2位置的阻力系数小于1位置的主要原因; 另外, 由于A尾耙在1位置时离翼型较近, 翼型扰动相对较大, 对阻力测量结果也会产生一定影响。
![]()
图 5 A尾耙在1、2位置的动量损失曲线(Ma=0.50)Fig. 5 The momentum loss curve for wake rake A at installation position 1 and 2 (Ma=0.50)![]()
图 6 A尾耙在1、2位置的动量损失曲线(Ma=0.80)Fig. 6 The momentum loss curve for wake rake A at installation position 1 and 2 (Ma=0.80)从动量损失曲线图还可以看出, 在小迎角条件下, 对尾迹动量损失曲线积分求阻力时, 2位置的积分计算结果更为精确, 因为1位置的动量损失曲线比2位置更为尖瘦。对比图 5和6的-2°和0°的动量损失曲线, 可以看出:1位置的尾迹中心区动量损失变化十分剧烈, 竖向范围较小, 尾迹中心区总压探针相对较少, 积分求阻力可能会产生一定误差; 而尾耙后移至2位置时, 相同状态下动量损失的峰值减小、竖向范围扩大, 尾迹中心区的总压探针增多, 积分结果更为精确。
3.2 不同尾耙高度对阻力测量的影响
上节实验结果验证了尾耙安装在1位置时过于靠近翼型, 阻力系数测量结果不精确; 当尾耙后移至2位置时, A尾耙的高度不够。本节以尾耙高度更高的B尾耙在2位置进行实验。
图 7给出了Ma=0.50、0.70条件下, A、B尾耙在2位置的阻力系数, 图 8和9为Ma=0.50、0.70条件下A、B尾耙在2位置的动量损失曲线。
![]()
图 7 A、B尾耙在2位置的阻力系数Fig. 7 Drag coefficients for wake rake A and B in installation position 2![]()
图 8 A、B尾耙在2位置的动量损失曲线(Ma=0.50)Fig. 8 The momentum loss curve for wake rake A and B in installation position 2 (Ma=0.50)![]()
图 9 A、B尾耙在2位置的动量损失曲线(Ma=0.70)Fig. 9 The momentum loss curve for wake rake A and B at installation position 2 (Ma=0.70)从图 7可以看出, 整体上, B尾耙的阻力系数略大于A尾耙, 这是因为B尾耙竖向高度较大, 洞壁干扰和堵塞干扰等相应增大, 尾耙本身对阻力测量结果造成较大影响。在大迎角时, 2位置的B尾耙的阻力系数比A尾耙大得多, 除因B尾耙高度较大造成的影响外, 根据图 8和9可以看出, 主要原因是此时A尾耙的高度并未完全覆盖尾迹动量损失范围, 计算所得的阻力系数是不正确的, 比B尾耙测得的阻力系数小是合理的。
另外, 从图 8和9还可以看出:B尾耙的动量损失曲线光滑性较差(尤其是迎角较小时), 积分求得的阻力系数可能存在一定误差; 随着迎角和马赫数增大, 曲线变得更为光滑, 积分所得的阻力系数误差随之减小。
4 结论
本文对西北工业大学NF-6增压连续式跨声速风洞翼型阻力测量精度改进方法进行了研究, 对比分析了2种尾耙高度和安装位置对阻力测量结果的影响, 得出以下结论:
(1) 当实验迎角和马赫数较小时, 采用较小高度的A尾耙(高270 mm)安装于2位置(尾耙总压探针前端距离翼型后缘2.90倍弦长)进行实验比较合适。
(2) 当实验迎角或实验马赫数较大时, 应考虑采用B尾耙(高500 mm)安装于2位置(尾耙总压探针前端距离翼型后缘2.90倍弦长)进行实验。
-
![]()
图 2 CO、CO2和H2O 2000~2250 cm-1内吸收谱线强度
Fig. 2 The absorption line -strength of CO、CO2 and H2O at the range of 2000~2250 cm-1
![]()
图 3 温度在1300 K时2050~2100 cm-1波段内CO的吸收谱线
Fig. 3 The absorption spectrum of CO in 2050~2100 cm-1 band at 1300 K
![]()
图 4 候选谱线对3和4在800~1300 K下的温度敏感性
Fig. 4 The temperature sensitivity of candidate lines for pair3 and pair4 at the temperature range of 800~1300 K
![]()
图 5 测量浓度与标准气体浓度的比较
Fig. 5 Comparison of measured concentration and standard gas concentration
![]()
图 6 采集信号的吸收光度拟合及残差
Fig. 6 Absorption spectrophotometric fitting and residual of acquisition signal
![]()
图 7 平面火焰稳定工况下的CO浓度测量结果
Fig. 7 CO concentration measurement under flame stabilization conditions
![]()
图 8 平面火焰炉不同位置的CO温度测量结果
Fig. 8 CO temperature measurement at different locations of the flat flame furnace
![]()
图 11 降噪处理后发动机加热过程和煤油点火过程的信号
Fig. 11 Signals of engine heating process and kerosene ignition process
![]()
图 14 发动机燃烧室压力和发动机出口CO温度随时间变化
Fig. 14 Change of engine combustion chamber pressure and engine outlet CO temperature with time
![]()
图 15 发动机燃烧室压力和发动机出口CO浓度随时间变化
Fig. 15 Change of engine combustion chamber pressure and engine outlet CO concentration with time
表 1 候选谱线参数
Table 1 The parameters of candidate lines
Spectral
pairCenter
frequency
/cm-1Line strength
/(10-20 cm-1·
(mol·cm-2)-1)Low state-
energy
/cm-1Temperature
/Kpair3 2059.9147
2060.332213.40
2.75806.38280
2543.056701300 pair4 2064.3969
2064.583913.90
2.71729.67740
2489.783101300 pair3 2059.9147
2060.332210.10
3.51806.38280
2543.056701800 pair4 2064.3969
2064.583910.20
3.41729.67740
2489.783101800 
下载: 导出CSV
表 2 不同工况下CO温度测量结果
Table 2 CO temperature measurement results under different working conditions
Order Air
/(L·min-1)CH4
/(L·min-1)ϕ Tm
/KTc
/KRelative
error1 18.9 2.39 1.2 1290.2 1237 4.3% 2 18.8 2.55 1.3 1276.7 1252 1.9% 3 18.5 2.72 1.4 1272.3 1248 1.9%
下载: 导出CSV
-
[1] 许琳, 王高, 吕国义, 等. 超声测温技术在模拟航空发动机燃烧室温度测量中的应用[J]. 测试技术学报, 2019, 33(2): 178-184. DOI: 10.3969/j.issn.1671-7449.2019.02.016 XU L, WANG G, LYU G Y, et al. Ultrasonic temperature measurement technology for simulated aero-engine combustion Chambers[J]. Journal of Test and Measurement Technology, 2019, 33(2): 178-184. doi: 10.3969/j.issn.1671-7449.2019.02.016
[2] 周蕊燕, 贾译钧. 浅析航空燃烧室及其发展趋势[J]. 山东工业技术, 2016(5): 214. doi:10.16640/j.cnki.37-1222/t.2016.05. 201 [3] 胡志云, 叶景峰, 张振荣, 等. 航空发动机地面试验激光燃烧诊断技术研究进展[J]. 实验流体力学, 2018, 32(1): 33-42. DOI: 10.11729/syltlx20170135 HU Z Y, YE J F, ZHANG Z R, et al. Development of laser combustion diagnostic techniques for ground aero-engine testing[J]. Journal of Experiments in Fluid Mechanics, 2018, 32(1): 33-42. doi: 10.11729/syltlx20170135
[4] SPEARRIN R M, JEFFRIES J B, HANSON R K. Mid-infrared absorption sensor for measurements of CO and CO2 in propulsion flows[R]. AIAA 2014-0390, 2014. doi: 10.2514/6.2014-0390
[5] WAKEFIELD S, TOM B A, SELL B C, et al. Development of a mid-IR laser diagnostic for combustion efficiency[R]. AIAA 2014-1356, 2014. doi: 10.2514/6.2014-1356
[6] BUSA K M, BROWN M S, GRUBER M, et al. Common-path measurement of H2O, CO, and CO2 via TDLAS for combustion progress in a hydrocarbon-fueled scramjet[R]. AIAA 2016-0659, 2016. doi: 10.2514/6.2016-0659
[7] 张伯昆. 基于TDLAS的环境CO红外激光光谱监测方法研究[D]. 合肥: 合肥工业大学, 2012. ZHANG B K. Study on infrared laser spectrum monitoring method of environmental CO based on TDLAS[D]. Hefei: Hefei University of technology, 2012.
[8] 邢昆明. 基于TDLAS技术检测卷烟主流烟气中CO和CO2及抽吸耗氧量[D]. 合肥: 中国科学技术大学, 2017. XING K M. Detection of CO and CO2 in mainstream cigarette smoke based on TDLAS technology and oxygen consumption by suction[D]. Hefei: University of Science and Technology of China, 2017.
[9] 张保龙. 基于TDLAS的CO浓度检测系统的研究[D]. 西安: 西安科技大学, 2019. ZHANG B L. Research on CO concentration monitoring system based on TDLAS[D]. Xi'an: Xi'an University of Science and Technology, 2019.
[10] SHAO L G, FANG B, ZHENG F, et al. Simultaneous detection of atmospheric CO and CH4 based on TDLAS using a single 2.3μm DFB laser[J]. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2019, 222: 117118. doi: 10.1016/j.saa.2019.05.023
[11] 许婷. 利用可调谐半导体激光吸收光谱测量高温火焰中的CO浓度[D]. 杭州: 浙江大学, 2011. XU T. Measurement of CO concentration in high temperature flame by tunable semiconductor laser absorption spectroscopy[D]. Hangzhou: Zhejiang University, 2011.
[12] 彭于权. 中红外激光光谱燃烧场CO和NO污染物浓度测量[D]. 合肥: 中国科学技术大学, 2018. PENG Y Q. Measurement of CO and NO concentration in combustion field of mid infrared laser spectrum[D]. Hefei: University of Science and Technology of China, 2018.
[13] 陈允魁. 红外吸收光谱法及其应用[M]. 上海: 上海交通大学出版社, 1993. [14] 吴静静. 基于可调谐激光吸收光谱技术的温度测量方法研究[D]. 哈尔滨: 哈尔滨理工大学, 2018. WU J J. Research on temperature measurement method base on tunable diode laser absorption spectroscopy technology[D]. Harbin: Harbin University of Science and Technology, 2018.
[15] GORDON I E, ROTHMAN L S, HILL C. The HITRAN2016 molecular spectroscopic database[J]. Journal of Quantitative Spectroscopy and Radiative Transfer, 2017, 203: 3-69. doi: 10.1016/j.jqsrt.2010.01.027
[16] 聂伟, 叶擎昊, 许振宇, 等. Voigt线型两翼拟合非均匀流场吸光度的方法研究[J]. 光谱学与光谱分析, 2017, 37(3): 816-821. DOI: 10.3964/j.issn.1000-0593(2017)03-0816-06 NIE W, YE Q H, XU Z Y, et al. Study on the method of Voigt profiles two wings fitting non-uniform flow field absorbance[J]. Spectroscopy and Spectral Analysis, 2017, 37(3): 816-821. doi: 10.3964/j.issn.1000-0593(2017)03-0816-06
[17] 符鹏飞, 超星, 侯凌云, 等. 基于TDLAS技术的煤油燃烧温度与组分分布检测[J]. 工程热物理学报, 2019, 40(9): 2176-2182. https://www.cnki.com.cn/Article/CJFDTOTAL-GCRB201909033.htm FU P F, CHAO X, HOU L Y, et al. Measurement of kerosene combustion temperature and component distribution based on TDLAS technology[J]. Journal of Engineering Thermophysics, 2019, 40(9): 2176-2182. https://www.cnki.com.cn/Article/CJFDTOTAL-GCRB201909033.htm
[18] ZHOU X, LIU X, JEFFRIES J B, et al. Development of a sensor for temperature and water concentration in combustion gases using a single tunable diode laser[J]. Measurement Science and Technology, 2003, 14(8): 1459-1468. doi: 10.1088/0957-0233/14/8/335
-
期刊类型引用(18)
1. 韩笑，高创，郑怡彤，刘小兵. 切角凹槽矩形截面桥塔的气动特性试验研究. 石家庄铁道大学学报(自然科学版). 2024(01): 9-15+23 . 
百度学术
2. 白桦，刘博祥，姬乃川，李加武. 节段模型二元端板合理尺寸估算方法. 振动与冲击. 2023(02): 312-320 . 百度学术
3. 王仰雪，刘庆宽，靖洪淼，李震，孙一飞. 倾斜栏杆对流线型箱梁涡激振动性能影响的试验研究. 振动与冲击. 2023(06): 232-239+254 . 百度学术
4. 于畅，仇法梅，田学东，杨群，刘小兵. 圆角矩形柱平均气动力特性的雷诺数效应试验研究. 石家庄铁道大学学报(自然科学版). 2023(02): 29-34 . 百度学术
5. 杨群，于畅，刘小兵，刘庆宽. 不同圆角率的方形断面斯特罗哈数的雷诺数效应研究. 振动与冲击. 2023(11): 223-231 . 百度学术
6. 刘庆宽，王仰雪，孙一飞，李震，韩原，靖洪淼. 栏杆高度对流线型箱梁涡振性能影响的试验研究. 湖南大学学报(自然科学版). 2023(07): 140-150 . 百度学术
7. 刘路路，杨皓然，邹云峰，何旭辉，韩艳，陈志强. 公铁同层双幅非对称主梁气动干扰特性研究. 铁道科学与工程学报. 2023(10): 3861-3872 . 百度学术
8. 杨群，于畅，于文文，刘小兵. 圆角方柱气动特性的风洞试验研究. 振动与冲击. 2023(24): 59-68 . 百度学术
9. 韩振，李波，甄伟，杨庆山，田玉基. 宽厚比为5的超高层建筑风荷载特性研究. 哈尔滨工程大学学报. 2022(02): 196-202 . 百度学术
10. 张庆华，马文勇，杨杰，张彦，周帅伟. 不同风倾角下典型等边角钢静风力风洞试验研究. 振动工程学报. 2022(02): 277-283 . 百度学术
11. 白桦，王涵，姬乃川，李加武. 节段模型长宽比对风洞测力试验及计算分析的影响. 中国公路学报. 2022(08): 202-212 . 百度学术
12. 杨群，刘庆宽，韩瑞，刘小兵. 不同圆角率的方形断面气动特性的雷诺数效应. 振动与冲击. 2020(04): 150-156 . 百度学术
13. 沈国辉，姚剑锋，郭勇，邢月龙，楼文娟. 直径30 cm圆柱的气动力参数和绕流特性研究. 振动与冲击. 2020(06): 22-28 . 百度学术
14. 温青，池俊豪，华旭刚，王修勇，孙洪鑫. 端部条件和展弦比对矩形断面节段模型气动力特征的影响. 实验流体力学. 2020(04): 36-43 . 本站查看
15. 任若松，梁新华，刘小兵，马文勇，刘庆宽. 准流线型桥梁断面气动力特性的雷诺数效应研究. 工程力学. 2020(S1): 139-144+167 . 百度学术
16. 李海飞，梁新华，孙一飞，崔会敏，刘庆宽. 流线型桥梁断面表面脉动风荷载特性研究. 工程力学. 2020(S1): 242-248+260 . 百度学术
17. 马文勇，汪冠亚，郑熙，陈铁，李智，张程远，方平治. 端部状态对斜置圆柱气动力分布的影响. 实验流体力学. 2019(02): 43-50 . 本站查看
34. 杨群，刘庆宽，孙亚松，刘小兵. 圆角方形断面气动特性试验. 振动.测试与诊断. 2020(01): 140-147+208 . 百度学术
其他类型引用(16)
计量
- 文章访问数: 389
- HTML全文浏览量: 137
- PDF下载量: 31
- 被引次数: 34
