Experimental research on the effect of processing technology of cascade test article on the performance measurement
-
摘要: 针对近年来平面叶栅试验研究过程中出现的试验件加工问题，对直叶片加工工艺、有机玻璃栅板加工工艺以及成套叶栅试验件关键参数的检测方法进行了工艺试验和叶栅试验器试验验证。结果表明，综合考虑加工成本和粗糙度对性能的影响规律，建议高速高负荷叶栅叶片表面粗糙度应不低于Ra=1.6。航空有机玻璃的选材、加工工序、工艺以及装配使用不当会导致试验过程中出现应力纹和黑圈，从而对高速状态下的可视化测量造成不利影响。通过改进成套叶栅试验件关键参数的检测方法，试验件加工质量有了较大提升，有助于获得更为准确可靠的性能试验数据。Abstract: In order to solve the processing problems of the test article occurred in the cascade experimental research, investigation and study of the processing technologies of the test blade and Aero-Plexiglass sidewall have been carried out through processing test and experimental validation. And then the detection methods of the key parameters are put forward based on trilinear coordinates measuring instrument. It can be seen from the results that considering the processing cost and the effects of the blade surface roughness on the cascade performance, the suitable surface roughness of the blade for high speed and high load cascade is suggested to be not more than Ra=1.6. Stress lines and black areas may appear on the Aero-Plexiglass sidewall due to some factors, such as material selections, manufacturing procedures, processing technologies and improper assembly or use, etc., which would cause adverse effect on schlieren visual measurement in high speed state. The processing quality of the test article is greatly improved by developing the detection methods of the key parameters, which will be helpful to get more accurate and reliable performance data.
-
Keywords:
- cascade test /
- processing technology /
- surface roughness /
- glass cascades /
- detecting technology
-
0 引言
目前，要求地面试验设备完全模拟高超声速飞行环境是非常困难的，在一个试验设备上进行所有环境的模拟试验更是不可能的。因此，需要研制多种试验设备，以满足高超声速试验需求[1-5]。
常规高超声速风洞主要模拟参数是马赫数和雷诺数[6]，其他一些参数则主要通过高焓试验设备来模拟[7-8]。
在高超声速风洞运行中，气源不仅要维持风洞所必须达到的压力比，还需要满足雷诺数模拟的要求。一般来说，风洞的总压要足够高。如果风洞连续工作，高的压力比和高总压将使风洞消耗很大的动力，因而高超声速风洞多为暂冲式。
暂冲式常规高超声速风洞的运行方式是：气罐中的压缩空气经过加热系统达到所需温度，然后通过型面喷管，在试验段形成所需的高超声速流场，最后经由超声速扩散段升压后进入引射排气系统，排入大气，或经过冷却器进入真空系统，由真空泵系统排入大气。常规高超声速风洞一般加热温度在800K以下，主要是为了防止气流冷凝，不能真实模拟实际总温(马赫数6时达到1800K)，因此，开展发动机及飞行器一体化带动力试验还需要满足总温要求的高焓设备。
燃烧加热风洞是目前高超声速飞行器地面模拟试验尤其是带动力试验的主要设备。在过去的十几年间，中国空气动力研究与发展中心研制了不同尺度的脉冲燃烧风洞，探索了一体化飞行器设计、计算与试验技术[7-9]。在此基础上，发展了基于脉冲燃烧风洞的大尺度飞行器带动力一体化试验技术[10-11]，提出了一种基于一体化试验直接测力结果的飞行器和发动机性能评估方法。
由于采用燃烧加热方式，无论是采用氢燃料还是碳氢燃料，都会产生水蒸气(H2O)、二氧化碳(CO2)等燃烧产物，即形成所谓“污染”。污染组分将造成风洞试验气体物理化学属性与真实空气存在一定差异，导致地面试验难以完全模拟真实飞行状态下的所有来流参数，且给地面试验结果向真实飞行状态的外推带来不确定性。为了尽可能降低这种不确定性，地面试验一般需要慎重考虑试验来流与模拟对象环境之间的参数匹配问题：即通过有选择地调整污染来流的某些状态参数，使之逼近对应的真实飞行环境，同时放弃一些无法兼顾的非关键参数，以尽可能达到飞行器气动与推进性能的可靠模拟。
目前主要有两种模拟方式：对于气动力试验，一般采用静温、静压、马赫数模拟；对于发动机性能试验，也可以采用总焓、动压、马赫数模拟。
中国科技大学杨基明、罗喜胜等[12-13]通过典型升力体飞行器的试验和计算研究，证明污染组分使得斜激波波后参数相对纯净空气有一定变化，从而造成模型表面压力及模型的气动力有一定变化，但变化量较小。相对于纯净空气来流，污染空气来流时发动机推力、单位推力、比冲均有所下降。脉冲燃烧风洞采用氢氧燃烧加热，污染空气来流时推力性能下降相对较小，与纯净空气来流时较为接近。
西北工业大学宋文艳等[14]研究了H2O/CO2污染煤油燃料对超声速燃烧室的影响，认为污染对点火和超声速燃烧具有一定的抑制作用，污染组分的存在会导致燃烧室模态转换点发生变化。
谭宇等[15]在酒精燃烧和氢气燃烧两种加热方式的风洞设备上开展了匹配方案对超燃冲压发动机性能影响的试验研究，比较了两种目前较常用的气流参数匹配方案，结果表明：对于采用氢气燃烧加热方式的风洞设备，总焓动压马赫数(h0QM)匹配比静温静压马赫数(TPM)匹配能够获得更高的壁面静压和推力收益；对于采用酒精燃烧加热方式的风洞设备，两种匹配方案表现相当。
壁温比通常定义为模型壁面温度与恢复温度的比值。常规高超声速风洞来流温度较低，气流壁面恢复温度与壁面温度相差不大，接近于绝热状态；而燃烧加热脉冲风洞来流温度较高，风洞试验时间很短(300~600ms)，壁面温度基本保持常温，接近于等温状态，气流壁面恢复温度与试验模型表面温差较大，气流对壁面的加热效应明显。壁面温度条件会影响到高超声速边界层内的流动参数，进而影响到流场的波系结构和气动性能。壁温降低引起摩擦系数的增加和边界层厚度减小，使得激波边界层作用区域变小；同时，冷壁使得边界层内亚声速部分声速较低，马赫数更高。
为了评估真实飞行条件下的飞行器气动性能，需重点考虑壁温对燃烧加热脉冲风洞试验结果的影响。
本文采用不通气标模，在常规高超声速风洞以及两个不同尺度(Φ2.4m和Φ600mm)的脉冲燃烧风洞中开展对比测力试验，结合数值计算，研究脉冲燃烧风洞水凝结、雷诺数及壁温比对模型气动性能的影响规律。
1 试验模型
试验模型采用不通气标模。对应Φ2.4m脉冲燃烧风洞，试验模型为大不通气标模，采用背部支撑，如图 1所示；对应常规高超声速风洞和Φ600mm脉冲燃烧风洞，试验模型为1/5缩比的小不通气标模，采用尾部支撑，如图 2所示。
2 脉冲燃烧风洞水凝结影响研究
为了获得水凝结对脉冲燃烧风洞试验数据的影响，采用小不通气标模，在Φ600mm脉冲燃烧风洞开展不同总温条件下的对比试验研究，分析测力结果。
2.1 对比试验参数
对Φ600mm脉冲燃烧风洞喷管进行了试验配套改造，在来流总温、组分不同的条件下，保持风洞出口马赫数一致。以现有的Ma6喷管为基础(总温1500K)，设计制造总温为1200和1800K的喉道段，喷管出口马赫数与现有喷管一致。表 1为不通气高超标模试验参数(p0、T0分别为总压、总温，p、T分别为静压、静温)。
表 1 不同总温高超标模试验参数Table 1 Flow parameters of different T0 testsMa p0/MPa T0/K p/Pa T/K Re 6.0 4.5 1200 2368 162 7.08×106 6.0 5.1 1500 2398 210 4.98×106 6.0 5.8 1800 2374 263 3.61×106 2.2 试验结果
图 3、4给出了不通气高超标模的轴向力系数CA、法向力系数CN随迎角α变化的试验结果曲线。
![]()
图 3 高超标模不同总温试验CA~α曲线Fig. 3 CA~α graph of different T0 tests for the typical hypersonic model![]()
图 4 高超标模不同总温试验CN~α曲线Fig. 4 CN~α graph of different T0 tests for the typical hypersonic model试验表明：在总温1200~1800K试验状态下，最大轴向力系数差别为5.4%，说明水凝结对不通气高超标模的轴向力影响很小，约为5%左右。在-2°迎角时，总温1200K(有水凝结)试验状态的法向力系数较总温1500K减小了约40%；在0°~6°迎角下，法向力系数减小12%~5%，说明随着迎角的增加，模型下表面的凝结水蒸发，导致水凝结对法向力的影响减小。在-2°迎角时，总温1800K试验状态的法向力系数较总温1500K增加了约30%，说明总温1500K、-2°迎角状态在模型的下表面也发生了水凝结，导致模型下表面流场压力升高，法向力系数减小；而在其他迎角状态，总温1500K时并没有发生水凝结现象，因此法向力系数差别较小。
3 脉冲燃烧风洞雷诺数影响研究
3.1 对比试验参数
为了获得雷诺数对试验结果的影响，用不同尺度的不通气标模在Φ2.4m脉冲燃烧风洞和Φ600mm脉冲燃烧风洞进行对比试验，试验参数见表 2(CHIF：脉冲燃烧风洞)。
表 2 试验参数Table 2 Test parametersFacility Ma T0/K p0/MPa Φ2.4m CHIF 6 1575 4.76 Φ600mm CHIF 6 1500 4.61 3.2 试验和计算结果
图 5和6给出了大、小不通气标模在两个脉冲燃烧风洞的轴向力系数CA、法向力系数CN试验结果，同时也给出了两个模型按照对应风洞来流条件计算获得的气动力系数结果。
![]()
图 5 不通气标模轴向力系数计算和试验结果Fig. 5 Numerical and experimental axis force coefficients of the test model![]()
图 6 不通气标模法向力系数计算和试验结果Fig. 6 Numerical and experimental normal force coefficients of the test model计算结果表明：模型尺度对不通气标模的轴向力有一定影响，小不通气标模的轴向力略大于大不通气标模，相差约2%~5%。由于两个脉冲燃烧风洞模拟条件基本相同，试验时间也接近，因此这个差异可以解释为模型尺度(雷诺数)造成的：小不通气标模的雷诺数小于大不通气标模，因此轴向力略大。试验并未能够反映出雷诺数的影响，这可能是由于两个模型支撑方式不同以及试验测量误差掩盖了差异。
两个模型的法向力结果不一致，主要是由于支撑方式不同造成的。前期研究表明，腹部支撑对轴向力影响不大，但对法向力影响显著。
计算与试验结果符合较好，使计算和试验得到了相互验证。
4 脉冲燃烧风洞壁温比影响研究
采用小不通气标模，通过Φ600mm脉冲燃烧风洞及常规高超声速风洞的试验对比，结合数值计算，研究壁温比对脉冲燃烧风洞试验结果的影响。
4.1 对比试验参数
Φ600mm脉冲燃烧风洞与Φ1m常规高超声速风洞的试验流场存在总温、雷诺数、壁温比等方面的差异。假设壁温300K，风洞流场参数对比如表 3所示(HWT：常规高超声速风洞；Tratio为壁温比)。
表 3 Φ1m常规高超声速风洞与Φ600mm脉冲燃烧风洞流场参数Table 3 Test parameters of Φ1m hypersonic wind tunnel and Φ600mm combustion heated impulse facilityFacility Ma p0/MPa T0/K γ Re Tratio Φ1m HWT 6.0 3.00 500 1.40 2.62×107 5.17 Φ600mm CHIF 6.0 4.65 1560 1.36 4.10×106 1.49 4.2 试验和计算结果
脉冲燃烧风洞试验与计算结果表明：采用等温壁300K条件计算获得的阻力系数与试验数据比较接近，误差在6%以内；绝热壁计算阻力系数较等温壁300K计算小5%(迎角6°)~15%(迎角-2°)；不同壁温条件对升力系数的计算结果影响不大，与试验结果的误差分别在6%和8%以内；采用等温壁300K条件计算获得的结果与试验基本一致。图 7和8给出了脉冲燃烧风洞试验和计算的阻力系数CD及升力系数CL的对比曲线。
![]()
图 7 CD~α试验与计算结果对比(脉冲风洞，总温1500K)Fig. 7 CD~α comparison of numerical and experimental results from Φ600mm combustion heated impulse facility![]()
图 8 CL~α试验与计算结果对比(脉冲风洞，总温1500K)Fig. 8 CL~α comparison of numerical and experimental results from Φ600mm combustion heated impulse facility常规高超声速风洞试验和计算结果表明：不同壁温条件对计算结果影响不大，计算结果与试验符合较好，阻力误差在4%以内，升力误差在3%以内。图 9和10给出了相应的阻力系数和升力系数比较曲线。
![]()
图 9 CD~α试验与计算结果对比(常规高超声速风洞)Fig. 9 CD~α comparison of numerical and experimental results from Φ1m hypersonic wind tunnel![]()
图 10 CL~α试验与计算结果对比(常规高超声速风洞)Fig. 10 CL~α comparison of numerical and experimental results from Φ1m hypersonic wind tunnel通过试验验证，表明数值计算具有较高的可信度。为了进一步分析规律，采用数值计算分析了脉冲燃烧风洞与常规高超声速风洞阻力系数的数据相关性，如图 11所示。
![]()
图 11 CD~α计算与试验结果对比(脉冲燃烧风洞和常规高超声速风洞)Fig. 11 CD~α comparison of numerical and experimental results from combustion heated impulse facility and hypersonic wind tunnel(1) 脉冲燃烧风洞等温壁计算得到的阻力系数比常规高超声速风洞大7%(迎角6°)~17%(迎角-2°)，迎角越小，差异越大。这个差异由多方面因素导致，主要有雷诺数、比热比、壁温比等。
(2) 针对同一个模型，分别采用脉冲燃烧风洞和常规高超声速风洞的模拟参数计算(壁面都取绝热壁条件)，此时得到的结果差异应该排除了壁温比的影响，而反映了两个风洞模拟条件(如雷诺数、比热比)不同引起的变化量，计算获得的这一部分的影响量约为3%(迎角6°)~6%(迎角-2°)。
(3) 除去两个风洞模拟条件(如雷诺数、比热比)不同引起的阻力变化量，可得到壁温比的影响量约为4%(迎角6°)~12%(迎角-2°)，迎角越小，差异越大。
5 结论
通过脉冲燃烧风洞与常规高超声速风洞不通气高超标模对比试验和计算，得到如下结论：
(1) 脉冲燃烧风洞获得的气动力系数变化规律与常规高超声速风洞一致；
(2) 雷诺数、壁温比对阻力系数均有影响，其中壁温比影响显著，脉冲燃烧风洞获得的阻力系数明显大于常规高超声速风洞；
(3) 水凝结对不通气高超标模的轴向力影响较小，对法向力影响较大，且不同迎角状态影响量不同。
-
![]()
图 2 不同试验雷诺数状态下叶片表面等熵马赫数分布
Fig. 2 Isentropic Mach number distribution at different Reynolds number states
![]()
图 3 i=0°、Re=1.29×106时吸力面油流图片
Fig. 3 Suction side oil flow pictures at i=0°, Re=1.29×106
![]()
图 4 不同试验马赫数状态下叶片表面等熵马赫数分布
Fig. 4 Isentropic Mach number distribution at different Mach number states
![]()
图 5 总压损失系数随进口马赫数的变化曲线
Fig. 5 Total pressure loss variation with inlet Mach number at different blade surface roughness states
![]()
图 7 存在抛光精磨问题的航空有机玻璃栅板
Fig. 7 Cascade models with Aero-Plexiglass sidewall which have polishing and fine grinding problems
![]()
图 8 孔槽塌边造成的黑圈示意图
Fig. 8 Schilieren picture with black areas due to the chamfer edge collapse
![]()
图 11 试验件尾缘直线度及栅距检测位置示意图
Fig. 11 Schematic diagram of detection methods for trailing edge straightness and pitch
![]()
图 12 喉宽及安装角测量示意图
Fig. 12 Schematic diagram of detection methods for the throat width and stagger angle
![]()
图 13 叶片表面等熵马赫数分布对比
Fig. 13 Comparison of isentropic Mach number distribution between adjacent channels
表 1 不同线切割方式对比
Table 1 Comparison of different wire-electrode cutting methods
走丝方式 走丝速度/(mm·s-1) 加工粗糙度(Ra) 特点 快走丝 10~20 2.4~3.2 电极丝高速往返运动, 精度低 中走丝 粗加工8~12
精加工1~31.2~1.6 多次切割，精度提高 慢走丝 ＜0.25 0.8及以上 电极丝低速单向运动，精度高 注：获得的加工粗糙度范围基于工艺试验 
下载: 导出CSV
-
[1] 金秀杰, 马艳玲, 高见. 压气机叶片型面精密数控加工技术[C]//第十五届中国科协年会第13分会场: 航空发动机设计、制造与应用技术研讨会论文集. 2013. [2] 刘维伟, 李杰光, 赵明, 等. 航空发动机薄壁叶片加工变形误差补偿技术研究[J]. 机械设计与制造, 2009(10): 175-177. DOI: 10.3969/j.issn.1001-3997.2009.10.070 LIU W W, LI J G, ZHAO M, et al. Research on the compensation of deformation error in NC machining of thin-walled blades[J]. Machinery Design & Manufacture, 2009(10): 175-177. doi: 10.3969/j.issn.1001-3997.2009.10.070
[3] 高丽敏, 蔡宇桐, 曾瑞慧, 等. 叶片加工误差对压气机叶栅气动性能的影响[J]. 推进技术, 2017, 38(3): 525-531. DOI: 10.13675/j.cnki.tjjs.2017.03.007 GAO L M, CAI Y T, ZENG R H, et al. Effects of blade machining error on compressor cascade aerodynamic performance[J]. Journal of Propulsion Technology, 2017, 38(3): 525-531. doi: 10.13675/j.cnki.tjjs.2017.03.007
[4] 李正, 余华蔚, 尹红顺, 等. 椭圆前缘锐化度对亚声速压气机叶片性能的影响[J]. 燃气涡轮试验与研究, 2018, 31(3): 14-17, 42. DOI: 10.3969/j.issn.1672-2620.2018.03.003 LI Z, YU H W, YIN H S, et al. The influence of ellipse leading edge sharpness on the performance of subsonic compressorblade[J]. Gas Turbine Experiment and Research, 2018, 31(3): 14-17, 42. DOI: 10.3969/j.issn.1672-2620.2018.03.003
[5] 曹传军, 邱毅, 李斌. 亚声叶型前缘形状对压气机气动性能的影响[J]. 燃气涡轮试验与研究, 2018, 31(6): 1-7. DOI: 10.3969/j.issn.1672-2620.2018.06.001 CAO C J, QIU Y, LI B. Influence of leading edge shape of subsonic airfoil on compressor aerodynamicperformance[J]. Gas Turbine Experiment and Research, 2018, 31(6): 1-7. DOI: 10.3969/j.issn.1672-2620.2018.06.001
[6] 张伟昊, 邹正平, 李维, 等. 叶型偏差对涡轮性能影响的非定常数值模拟研究[J]. 航空学报, 2010, 31(11): 2130-2138. https://www.cnki.com.cn/Article/CJFDTOTAL-HKXB201011005.htm ZHANG W H, ZOU Z P, LI W, et al. Unsteady numerical simulation investigation of effect of blade profile deviation on turbine performance[J]. Acta Aeronautica et Astronautica Sinica, 2010, 31(11): 2130-2138. https://www.cnki.com.cn/Article/CJFDTOTAL-HKXB201011005.htm
[7] 张伟昊, 邹正平, 刘火星, 等. 叶型偏差对整机环境中涡轮性能的影响[J]. 工程热物理学报, 2010, 31(11): 1830-1834. https://www.cnki.com.cn/Article/CJFDTOTAL-GCRB201011010.htm ZHANG W H, ZOU Z P, LIU H X, et al. Effect of profile deviation on turbine performance in whole engine environment[J]. Journal of Engineering Thermophysics, 2010, 31(11): 1830-1834. https://www.cnki.com.cn/Article/CJFDTOTAL-GCRB201011010.htm
[8] BACK S C, HOBSON G V, SONG S J, et al. Effects of Reynolds number and surface roughness magnitude and location on compressor cascade performance[J]. Journal of Turbomachinery, 2012, 134(5): 051013. doi: 10.1115/1.4003821
[9] SCHLICHTING H, GERSTEN K. Boundary-layer theory[M]. New York: McGraw-Hill, 1979.
[10] MONTIS M, NIEHUIS R, FIALA A. Aerodynamic measurements on a low pressure turbine cascade with different levels of distributed roughness[C]//Proceedings of ASME 2011 Turbo Expo: Turbine Technical Conference and Exposition. 2012. doi: 10.1115/GT2011-45015
[11] 余心明, 戴青. 影响线切割加工表面粗糙度的因素及应对措施[J]. 科技与企业, 2013(22): 392. doi: 10.13751/j.cnki.kjyqy.2013.22.365 [12] 杨蕾, 蔡安江. 电火花线切割加工参数对加工速度和表面粗糙度影响的研究[J]. 机床与液压, 2011, 39(15): 45-47. DOI: 10.3969/j.issn.1001-3881.2011.15.014 YANG L, CAI A J. Research on influence of processing parameters on processing speed and surface roughness in WEDM[J]. Machine Tool & Hydraulics, 2011, 39(15): 45-47. DOI: 10.3969/j.issn.1001-3881.2011.15.014
[13] 周桂莲, 张迪. 汽轮机叶片材料1Cr13钢线切割表面粗糙度正交试验研究[J]. 工具技术, 2010, 44(5): 33-35. DOI: 10.16567/j.cnki.1000-7008.2010.05.008 ZHOU G L, ZHANG D. Orthogonal experiment investigation on surface roughness of WEDM in machining 1Cr13 steel turbineblade[J]. Tool Engineering, 2010, 44(5): 33-35. doi: 10.16567/j.cnki.1000-7008.2010.05.008
[14] 万枝铭. 有机玻璃零件的加工[J]. 机械工程师, 2007(5): 147-148. DOI: 10.3969/j.issn.1002-2333.2007.05.076 WAN Z M. The processing of PMMA components[J]. Mechanical Engineer, 2007(5): 147-148. DOI: 10.3969/j.issn.1002-2333.2007.05.076
[15] 万庆, 李晓波, 李鹏飞, 等. 有机玻璃零件的数控加工[J]. 金属加工(冷加工), 2018(11): 28-29. DOI: 10.3969/j.issn.1674-1641.2018.11.014 WAN Q, LI X B, LI P F, et al. The CNC machining of PMMA components[J]. Metal Working (Metal Cutting), 2018(11): 28-29. DOI: 10.3969/j.issn.1674-1641.2018.11.014
[16] 唐凯, 刘志刚, 王晖, 等. 叶片表面粗糙度对前弯压气机叶片流动特性影响的试验研究[J]. 推进技术, 2020, 41(8): 1710-1719. DOI: 10.13675/j.cnki.tjjs.190486 TANG K, LIU Z G, WANG H, et al. Experimental investigation of surface roughness on flow around forward-skewed compressor blade[J]. Journal of Propulsion Technology, 2020, 41(8): 1710-1719. doi: 10.13675/j.cnki.tjjs.190486
[17] 孙玉福. 实用工程材料手册[M]. 北京: 机械工业出版社, 2014: 353. SUN Y F. Practical engineering material manual[M]. Beijing: China Machine Press, 2014: 353.
[18] 中国航空材料手册编辑委员会. 中国航空材料手册6[M]. 北京: 中国标准出版社, 1989.
计量
- 文章访问数: 261
- HTML全文浏览量: 179
- PDF下载量: 10
