A dynamic calibration method for a dynamometric system in impulse combustion facilities
-
摘要: 针对脉冲燃烧风洞中的测力系统,提出了一种动态标定方法。利用力锤在模型表面上不同位置,沿不同方向施加一系列集中载荷,由输入载荷和天平输出辨识出该表面对应的单位脉冲响应函数(UIRF),再将各表面对应的UIRF加权得到系统的UIRF,加权系数由试验状态下各表面的压力分布确定。辨识某表面对应的UIRF时,通过将其参数化使反卷积问题转化为参数优化问题以回避问题的病态特性。求解参数优化问题时,先用遗传算法搜索到参数全局最优解的近似值,再以此作为单纯形方法的初值继续优化得到参数最优值。在ANSYS中模拟了动态标定过程,考虑了实际试验中输出应变含有较大噪声的情况,验证了这种动态标定方法的准确性。Abstract: A new dynamic calibration method for a dynamometric system in impulse combustion facilities is proposed. The calibration involved uses an instrumented impact hammer to apply individual calibration forces in different directions at different positions on a face of the model and calculates the unit impulse response function (UIRF) of the face from input loads and output strains. UIRFs of different faces are weighted to obtain the UIRF of the dynamometric system and the weighting coefficients are determined by the pressure on each face under the experimental condition. By parameterization, the problem is converted into a parameter optimization problem to solve the UIRF. Using a genetic algorithm to obtain the approximation of the global optimal solution of the parameters and setting it as the initial value of a simplex algorithm, the exact solution is obtained by the simplex algorithm then. ANSYS simulation of the dynamic calibration is presented. Input loads and output strains are recorded and noises are added to the output strains to simulate the actual experimental situation. The simulation validates the accuracy and feasibility of the dynamic calibration method.
-
表 1 试验系统的前6阶模态
Table 1. The first six order modes of the testing system
Order 1 2 3 4 5 6 f/Hz 5.91 9.55 12.48 16.13 17.03 23.30 Mode Description around x axis around y axis around z axis along z axis along x axis along y axis 表 2 前10个振动成分
Table 2. The first ten frequencies and amplitudes
f/Hz A/με 12.29 15.41 5.77 6.87 16.92 6.03 9.39 -4.50 23.01 3.17 23.20 0.12 16.76 0.12 12.42 0.11 22.85 0.11 5.83 -0.08 表 3 优化初值
Table 3. Initial values of optimization
f/Hz A/με γ 5.77 6.80 0 9.40 -4.50 0 12.29 15.52 0 16.92 6.14 0 23.01 3.40 0 表 4 优化结果
Table 4. Results of optimization
f/Hz A/με γ 5.77 6.89 2×10-4 9.39 -4.50 4×10-4 12.29 15.48 5×10-5 16.92 6.19 2×10-4 23.01 3.46 4×10-4 表 5 优化结果
Table 5. Results of optimization
f/Hz A/με γ 5.77 6.86 1×10-4 9.34 -4.50 5×10-4 12.29 15.47 4×10-5 16.92 6.17 2×10-4 23.01 3.46 4×10-4 -
[1] 乐嘉陵, 刘伟雄, 贺伟, 等.脉冲燃烧风洞及其在火箭和超燃发动机研究中的应用[J].实验流体力学, 2005, 19(1):1-10. http://www.syltlx.com/CN/abstract/abstract9335.shtmlLe J L, Liu W X, He W, et al. Impulse combustion wind tunnel and its application in rocket and scramjet research[J]. Journal of Experiments in Fluid Mechanics, 2005, 19(1):1-10. http://www.syltlx.com/CN/abstract/abstract9335.shtml [2] 刘伟雄, 谭宇, 毛雄兵, 等.一种新运行方式脉冲燃烧风洞研制及初步应用[J].实验流体力学, 2007, 21(4):59-64. http://www.syltlx.com/CN/abstract/abstract9590.shtmlLiu W X, Tan Y, Mao X B, et al. The development and preliminary application of a pulse combustion wind tunnel with new running way[J]. Journal of Experiments in Fluid Mechanics, 2007, 21(4):59-64. http://www.syltlx.com/CN/abstract/abstract9590.shtml [3] Robinson M J, Mee D J, Tsai C Y, et al. Three-component force measurements on a large scramjet in a shock tunnel[J]. Journal of Spacecraft and Rockets, 2004, 41(3):416-425. doi: 10.2514/1.10699 [4] Robinson M J, Hannemann K, Schramm J M. Design and implementation of an internal stress wave force balance in a shock tunnel[J]. CEAS Space Journal, 2011, 1(1):45-57. doi: 10.1007/s12567-010-0003-5/fulltext.html [5] 贺伟, 于时恩, 李宏斌.高超声速一体化飞行器推阻特性测量研究[J].实验流体力学, 2010, 24(2):65-68. http://www.syltlx.com/CN/abstract/abstract9811.shtmlHe W, Yu S E, Li H B. Experimental investigation on thrust drag performance of hypersonic integrative vehicle[J]. Journal of Experiments in Fluid Mechanics, 2010, 24(2):65-68. http://www.syltlx.com/CN/abstract/abstract9811.shtml [6] 贺伟, 童泽润, 李宏斌.单模块超燃发动机推力测量天平研制[J].航空动力学报, 2010, 25(10):2285-2289. http://www.cnki.com.cn/Article/CJFDTOTAL-HKDI201010022.htmHe W, Tong Z R, Li H B. Investigation of thrust balance for the single module scramjet[J]. Journal of Aerospace Power, 2010, 25(10):2285-2289. http://www.cnki.com.cn/Article/CJFDTOTAL-HKDI201010022.htm [7] Abdel-Jawad M M, Mee D J, Morgan R G. New calibration technique for multiple-component stress wave force balances[J]. Review of Scientific Instruments, 2007, 78(6):065101-1-065101-7. doi: 10.1063/1.2744235 [8] 王锋, 任虎, 周正, 等.载荷辨识方法用于脉冲风洞模型阻力测量研究[J].振动与冲击, 2015, 34(24):202-208. http://www.cnki.com.cn/Article/CJFDTOTAL-ZDCJ201524034.htmWang F, Ren H, Zhou Z, et al. Drag force measurement in impulse facilities by using load identification method[J]. Journal of Vibration and Shock, 2015, 34(24):202-208. http://www.cnki.com.cn/Article/CJFDTOTAL-ZDCJ201524034.htm [9] 李东旭.高等结构动力学[M].北京:科学出版社, 2010:273-298.Li D X. Advanced structural dynamics[M]. Beijing:The Science Publishing Company, 2010:273-298. [10] Doherty L J, Smart M K, Mee D J. Measurement of three-components of force on an airframe integrated scramjet at Mach 10[R]. AIAA-2015-3523, 2015. [11] 刘国春, 费强, 赵武云, 等.基于Matlab遗传算法优化工具箱的应用[J].机械研究与应用, 2014, 27(2):71-73. http://www.cnki.com.cn/Article/CJFDTOTAL-JXYJ201402023.htmLiu G C, Fei Q, Zhao W Y, et al. Application of genetic algorithm optimization toolbox based on matlab[J]. Mechanical Research and Application, 2014, 27(2):71-73. http://www.cnki.com.cn/Article/CJFDTOTAL-JXYJ201402023.htm [12] 林鸿彬.基于遗传算法的数据拟合在MATLAB环境中的实现[J].湖南农机, 2010, 37(3):92-97. http://www.cnki.com.cn/Article/CJFDTOTAL-HNNJ201005046.htmLin H B. Data fitting based on genetic algorithm implementation in MATLAB environment[J]. Hunan Agricultural Machinery, 2010, 37(3):92-97. http://www.cnki.com.cn/Article/CJFDTOTAL-HNNJ201005046.htm [13] 罗述全.传统优化算法与遗传算法的比较[J].湖北工业大学学报, 2007, 22(3):32-35. http://www.cnki.com.cn/Article/CJFDTOTAL-HBGX200703011.htmLuo S Q. Comparison between traditional optimized algorithm and heredity algorithm[J]. Hubei University of Technology Journal, 2007, 22(3):32-35. http://www.cnki.com.cn/Article/CJFDTOTAL-HBGX200703011.htm [14] 郭海双, 梁佳雯, 张劭昀. MATLAB遗传算法工具箱GADS优化及应用[J].电子设计工程, 2015, 23(10):27-30. doi: 10.3969/j.issn.1674-6236.2015.10.009Guo H S, Liang J W, Zhang S Y. Optimization and examples in Matlab GA toolbox GADS[J]. Electronic Design Engineering, 2015, 23(10):27-30. doi: 10.3969/j.issn.1674-6236.2015.10.009 [15] 杨改强, 霍丽娟, 杨国义, 等.利用MATLAB拟合van Genuchten方程参数的研究[J].土壤, 2010, 42(2):268-274. http://www.cnki.com.cn/Article/CJFDTOTAL-TURA201002020.htmYang G Q, Huo L J, Yang G Y, et al. Research on fitting van genuchten equation parameter with MATLAB software[J]. Soils, 2010, 42(2):268-274. http://www.cnki.com.cn/Article/CJFDTOTAL-TURA201002020.htm