Quadratic programming equivalent mapping method for external aerodynamic force in flight load analysis
-
摘要:
传统的气动力映射方法中,守恒型方法能保证总积分载荷的守恒,但计算复杂;非守恒型方法实现简单,能保证气动压力的分布,但不能保证总积分载荷的守恒。利用等式约束保证积分载荷相等,并基于非守恒型加入线性基的径向基函数(RPIM)构造最小二乘保持原气动压力的分布特征,将气动力等效传递问题转化为一个带等式约束的二次规划问题。通过在非守恒方法的基础上构造出新的守恒型方法,兼顾了2种方法的优点。对某一大展弦比机翼上表面的三维气动力等效映射至二维平板面元气动力的结果表明,所提方法在牺牲了一定气动压力分布精度的基础上,严格保证了总积分载荷的相等,同时也继承了非守恒插值方法的高效性。所提方法可用于各类复杂外部气动力数据的快速等效映射,以及其他学科类似的数据映射问题。
Abstract:In conventional aerodynamic force mapping methods, conservative interpolation method can ensure the equivalence of the integral load, but the calculation is complex. While the non-conservative interpolation method is simple to implement with consistency of the distribution of aerodynamic pressure, but no conservation of the overall load. In this paper, equality constraints are used to ensure that the integral loads are equal, with the least squares based on point interpolation method using radical basis function (RPIM) to maintain the distribution characteristics of the original aerodynamic pressure, thus transforming the aerodynamic force mapping problem into a quadratic programming problem with equality constraints. This operation builds up a new conservative method on an original non-conservative method, which has preserved the advantages of each method. Mapping the aerodynamic force from a three-dimensional mesh to low-order panel mesh of a high-aspect-ratio wing is taken as an example. The results show that the method proposed in this paper can not only strictly guarantee the equality of integral loads, but also inherit the high efficiency of non-conservative interpolation method, despite the loss of some accuracy in the distribution of aerodynamic force. Besides, it is suitable for fast equivalent mapping of various complex external aerodynamic data, and similar data mapping problem in other disciplines.
-
表 1 径向基函数
Table 1. Radial basis function
类型 R(x)定义式 高斯型(G) e-α||x||2 对数型(TPS) ||x||2ln||x|| 复合2次型 (c2+||x||2)1/2 逆复合2次型 1/(c2+||x||2)1/2 Wendland’s C0 (1-||x||)+2 Wendland’s C2 (1-||x||)+4(4||x||+1) Wendland’s C4 (1-||x||)+6(35||x||2+18||x||+3) 表 2 映射前后总积分载荷系数
Table 2. Total integral load coefficient before and after mapping
系数 原气动
网格RPIM插值方法
(误差)二次规划方法
(误差)Cz 0.149 08 0.147 0(1.4%) 0.149 08(0%) CMx 0.027 657 5 0.027 081(2.08%) 0.027 657 5(0%) CMy -0.060 124 -0.059 259(1.44%) -0.060 124(0%) -
[1] 《飞机设计手册》总编委会.飞机设计手册第6册:气动设计[M].北京:航空工业出版社, 2001:49-186.General Editorial Board of Aircraft Design Manual.Aircraft design manual(6th volume):Aerodynamic design[M].Beijing:Aviation Industry Press, 2001:49-186(in Chinese). [2] 万志强, 邓立东, 杨超, 等.基于非线性试验气动力的飞机静气动弹性响应分析[J].航空学报, 2005, 26(4):439-445. doi: 10.3321/j.issn:1000-6893.2005.04.012WAN Z Q, DENG L D, YANG C, et al.Aircraft static aeroelastic response analysis based on nonlinear experimental aerodynamic data[J].Acta Aeronautica et Astronautica Sinica, 2005, 26(4):439-445(in Chinese). doi: 10.3321/j.issn:1000-6893.2005.04.012 [3] 安伟刚, 梁生云, 陈殿宇.一种局部动态数据交换方法在流固耦合分析中的应用[J].航空学报, 2013, 34(3):541-546. http://d.old.wanfangdata.com.cn/Periodical/hkxb201303010AN W G, LIANG S Y, CHEN D Y.Local dynamic data exchange in fluid structure interaction analysis[J].Acta Aeronautica et Astronautica Sinica, 2013, 34(3):541-546(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/hkxb201303010 [4] 黄宇, 阎超, 王文, 等.混合重叠网格插值方法的改进及应用[J].北京航空航天大学学报, 2017, 43(2):285-292. http://d.old.wanfangdata.com.cn/Periodical/bjhkhtdxxb201702010HUANG Y, YAN C, WANG W, et al.An improved interpolation method for hybrid overset grid and its application[J].Journal of Beijing University of Aeronautics and Astronautics, 2017, 43(2):285-292(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/bjhkhtdxxb201702010 [5] YAO Z, LIOU M S.Progress in the three-dimensional dragon grid scheme: AIAA-2001-2540[R].Reston: AIAA, 2001. [6] XU K L, SUN G, CAI J Y.On interface conservative attributions for computations of the complex flows of high-lift system based on chimera technique[C]//International Conference on Computer and Automation Engineering.Piscataway, NJ: IEEE Press, 2010: 279-283. [7] 张宇飞, 陈海昕, 符松.基于高阶守恒重映对窗口嵌入技术的改进[J].计算物理, 2011, 28(2):167-173. doi: 10.3969/j.issn.1001-246X.2011.02.002ZHANG Y F, CHEN H X, FU S.Improvement on window embedment technology with high order conservative remapping[J].Chinese Journal of Computational Physics, 2011, 28(2):167-173(in Chinese). doi: 10.3969/j.issn.1001-246X.2011.02.002 [8] XIANG Z.An implicit and globally conservative unstructured chimera grid method: AIAA-2011-777[R].Reston: AIAA, 2011. [9] FARRELL P E, MADDISON J R.Conservative interpolation between volume meshes by local galerkin projection[J].Computer Methods in Applied Mechanic and Engineering, 2011, 200(1):89-100. http://cn.bing.com/academic/profile?id=aac95ff70f0ddb440fb234e1c139f5d9&encoded=0&v=paper_preview&mkt=zh-cn [10] 徐春光, 董海波, 刘君.基于单元相交的混合网格精确守恒插值方法[J].爆炸与冲击, 2016, 36(3):305-312. http://d.old.wanfangdata.com.cn/Periodical/bzycj201603003XU C G, DONG H B, LIU J.An accurate conservative interpolation method for the mixed grid based on the intersection of grid cells[J].Explosion and Shock Waves, 2016, 36(3):305-312(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/bzycj201603003 [11] 田书玲.基于非结构网格方法的重叠网格算法研究[D].南京: 南京航空航天大学, 2008: 63.TIAN S L.Investigation of overset unstructured grids algorithm[D].Nanjing: Nanjing University of Aeronautics and Astronautics, 2008: 63(in Chinese). [12] ALLEN C B, RENDALL T C S.Unified fluid-structure interpolation and mesh motion using radial basis functions[J].International Journal for Numerical Methods in Engineering, 2008, 74:1519-1559. doi: 10.1002/nme.2219 [13] 林言中, 陈兵, 徐旭.基于径向基函数插值的气动弹性计算方法[J].北京航空航天大学学报, 2014, 40(7):953-958. http://d.old.wanfangdata.com.cn/Periodical/bjhkhtdxxb201407016LIN Y Z, CHEN B, XU X.Numerical method of aeroelasticity based on radial basis function interpolation[J].Journal of Beijing University of Aeronautics and Astronautics, 2014, 40(7):953-958(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/bjhkhtdxxb201407016 [14] 孙岩, 孟德虹, 王运涛, 等.基于径向基函数与混合背景网格的动态网格变形方法[J].航空学报, 2016, 37(5):1462-1472. http://d.old.wanfangdata.com.cn/Periodical/hkxb201605007SUN Y, MENG D H, WANG Y T, et al.Dynamic grid deformation method based on radial basis function and hybrid background grid[J].Acta Aeronautica et Astronautica Sinic, 2016, 37(5):1462-1472(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/hkxb201605007 [15] LIU G R, GU Y T.无网格法理论及程序设计[M].王建明, 周学军, 译.济南: 山东大学出版社, 2007: 45-48.LIU G R, GU Y T.An Introduction to meshfree methods and their programming[M].WANG J M, ZHOU X J, translated.Jinan: Shandong University Press, 2007: 45-48(in Chinese). [16] 陈文, 傅卓佳, 魏星.科学与工程计算中的径向基函数方法[M].北京:科学出版社, 2014:11-13.CHEN W, FU Z J, WEI X.Radial basis function method in scientific and engineering computation[M].Beijing:Science Press, 2014:11-13(in Chinese). [17] WENDLAND H.Piecewise polynomial, positive definite and compactly supported radial basis functions of minimal degree[J].Advances Computational Mathematics, 1995, 4(1):389-396. doi: 10.1007/BF02123482 [18] SMITH M J, HODGES D H, CESNIK C E S.Evaluation of computational algorithms suitable for fluid-structure interactions[J].Journal of Aircraft, 2000, 37(2):282-294. doi: 10.2514/2.2592 [19] 高立.数值最优化方法[M].北京:北京大学出版社, 2014:217-226.GAO L.Numerical optimization method[M].Beijing:Peking University Press, 2014:217-226(in Chinese). [20] BOER A D, ZUIJLEN A H V, BIJL H.Review of coupling methods for non-matching meshes[J].Computer Methods in Applied Mechanics and Engineering, 2007, 196(8):1515-1525. doi: 10.1016/j.cma.2006.03.017