A Method for Inverse Solution of Structural Interval Parameters based on DIRECT Algorithm
-
摘要: 提出一种基于DIRECT算法的结构区间参数反求方法。对于结构不确定性参数反求问题,一般转化为不确定性传播和模型参数优化的双层求解问题。首先,区间模型用来描述响应和待识别结构参数的不确定性,并建立了相应区间参数反求模型。其次,在迭代反求过程中自适应更新径向基函数用来近似原系统模型,并利用DIRECT算法来求解内层不确定性传播问题。最后,通过遗传算法来求解外层的优化模型,从而识别结构不确定性参数的区间。数值算例用来验证了该方法的正确性和有效性,并将该方法应用来反求车辆乘员约束系统中的不确定性参数。Abstract: This paper presents a structural interval parameter inverse solution method based on DIRECT algorithm. The inverse solution of structural interval parameters is generally translated into a two-layer solution of uncertainty analysis and model parameter optimization. First, an interval model is applied to describing the uncertainties of measured responses and identified structural parameters, and the corresponding inverse model of structural interval parameters is established. Second, in the process of iterative inversion, the adaptive updating radial basis function is applied to approximating the original system model, and the DIRECT algorithm is used to analyze the uncertainty of the inner layer. Finally, the optimization model of the outer layer is solved with the genetic algorithm to identify the upper and lower bounds of structural interval parameters. Numerical examples are used to verify the correctness and effectiveness of the method presented in the paper, which is applied to identifying the uncertainty parameters in the vehicle occupant constraint system.
-
表 1 不同核函数形式
序号 名称 核函数 形参数 1 高斯函数(Gaussian/GS) e-αr2 α 2 MQ函数(Multiquadric) (r2+R2)q R, q 3 对数路径函数(Logistic) 1/(1+eηr) η 4 立方函数(Cubic) (r2+λ)3 λ 5 IMQ函数(Inverse multiquadric) β 表 2 计算与测量响应结果对比
方法 计算响应 测量响应 下界误差/% 上界误差/% 本文方法 Y1 [19.973, 49.173] [20.000, 49.000] 0.14 0.35 Y2 [28.021, 57.684] [28.000, 57.500] 0.08 0.32 基于一阶泰勒展开的方法 Y1 [22.250, 51.776] [20.000, 49.000] 11.25 5.67 Y2 [30.588, 59.597] [28.000, 57.500] 9.24 3.65 表 3 计算与测量响应结果对比
WIC 计算响应 测量响应 下界误差/% 上界误差/% WIC1 [0.510, 0.544] [0.508, 0.542] 0.39 0.37 WIC2 [0.437, 0.472] [0.435, 0.469] 0.46 0.64 -
[1] Jin B. Conjugate gradient method for the robin inverse problem associated with the Laplace equation[J]. International Journal for Numerical Methods in Engineering, 2007, 71(4):433-453 doi: 10.1002/(ISSN)1097-0207 [2] Liu J, Xu C, Han X, et al. Determination of the state parameters of explosive detonation products by computational inverse method[J]. Inverse Problems in Science and Engineering, 2016, 24(1):22-41 doi: 10.1080/17415977.2014.993981 [3] Tarantola A. Popper, Bayes and the inverse problem[J]. Nature Physics, 2006, 2(8):492-494 doi: 10.1038/nphys375 [4] Fonseca J R, Friswell M I, Mottershead J E, et al. Uncertainty identification by the maximum likelihood method[J]. Journal of Sound and Vibration, 2005, 288(3):587-599 doi: 10.1016/j.jsv.2005.07.006 [5] Nichols J M, Link W A, Murphy K D, et al. A Bayesian approach to identifying structural nonlinearity using free-decay response:application to damage detection in composites[J]. Journal of Sound and Vibration, 2010, 329(15):2995-3007 doi: 10.1016/j.jsv.2010.02.004 [6] Zhang W, Liu J, Cho C, et al. A fast Bayesian approach using adaptive densifying approximation technique accelerated MCMC[J]. Inverse Problems in Science and Engineering, 2016, 24(2):247-264 doi: 10.1080/17415977.2015.1017488 [7] Liu J, Hu Y F, Xu C, et al. Probability assessments of identified parameters for stochastic structures using point estimation method[J]. Reliability Engineering & System Safety, 2016, 156:51-58 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=49c47a1144aa241eb38f59269c6d3d4a [8] 邱志平, 马丽红, 王晓军.不确定非线性结构动力响应的区间分析方法[J].力学学报, 2006, 38(5):645-651 doi: 10.3321/j.issn:0459-1879.2006.05.011Qiu Z P, Ma L H, Wang X J. Interval analysis for dynamic response of nonlinear structures with uncertainties[J]. Chinese Journal of Theoretical and Applied Mechanics, 2006, 38(5):645-651(in Chinese) doi: 10.3321/j.issn:0459-1879.2006.05.011 [9] 王晓军, 邱志平.含不确定参数弹簧质量系统振动反问题的区间分析法[J].固体力学学报, 2004, 25(4):461-466 doi: 10.3969/j.issn.0254-7805.2004.04.017Wang X J, Qiu Z P. Inverse vibration problem of spring-mass system with parameter uncertainty via interval analysis[J]. Acta Mechanica Solida Sinica, 2004, 25(4):461-466(in Chinese) doi: 10.3969/j.issn.0254-7805.2004.04.017 [10] Jiang C, Liu G R, Han X. A novel method for uncertainty inverse problems and application to material characterization of composites[J]. Experimental Mechanics, 2008, 48(4):539-548 doi: 10.1007/s11340-007-9081-5 [11] Liu J, Han X, Jiang C, et al. Dynamic load identification for uncertain structures based on interval analysis and regularization method[J]. International Journal of Computational Methods, 2011, 8(4):667-683 doi: 10.1142/S0219876211002757 [12] 郭红玲, 杨海天, 赵潇.蚁群算法求解弹性本构参数区间反问题[J].工程力学, 2012, 29(1):7-12, 19 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201200066327Guo H L, Yang H T, Zhao X. Solving an inverse problem of intervals of elastic constitutive parameters via ant colony algorithm[J]. Engineering Mechanics, 2012, 29(1):7-12, 19(in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201200066327 [13] Zhang W, Liu J, Han X, et al. A computational inverse technique for determination of encounter condition on projectile penetrating multilayer medium[J]. Inverse Problems in Science and Engineering, 2012, 20(8):1195-1213 doi: 10.1080/17415977.2012.659733 [14] Feng X X, Zhuo K M, Wu J L, et al. A new interval inverse analysis method and its application in vehicle suspension design[J]. SAE International Journal of Materials and Manufacturing, 2016, 9(2):315-320 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=fb4bec3388081bd0856208f187e11e88 [15] Zhu H, Bogy D B. DIRECT algorithm and its application to slider air-bearing surface optimization[J]. IEEE Transactions on Magnetics, 2002, 38(5):2168-2170 doi: 10.1109/TMAG.2002.802794 [16] 陈国栋, 韩旭, 刘桂萍, 等.基于自适应径向基函数的整车耐撞性多目标优化[J].中国机械工程, 2011, 22(4):488-493 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=SciencePaper201303280000000484Chen G D, Han X, Liu G P, et al. Multi-objective design optimization on crashworthiness of full vehicle based on adaptive radial basis function[J]. China Mechanical Engineering, 2011, 22(4):488-493(in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=SciencePaper201303280000000484 [17] Tian L F, Collins C. An effective robot trajectory planning method using a genetic algorithm[J]. Mechatronics, 2004, 14(5):455-470 doi: 10.1016/j.mechatronics.2003.10.001 [18] Wang J G, Liu G R. On the optimal shape parameters of radial basis functions used for 2-D meshless methods[J]. Computer Methods in Applied Mechanics and Engineering, 2002, 191(23-24):2611-2630 doi: 10.1016/S0045-7825(01)00419-4 [19] 陈国栋.基于代理模型的多目标优化方法及其在车身设计中的应用[D].长沙: 湖南大学, 2012 http://cdmd.cnki.com.cn/Article/CDMD-10532-1012482435.htmChen G D. Multi-objective optimization method based on metamodel and its applications in vehicle body design[D]. Changsha: Hunan University, 2012(in Chinese) http://cdmd.cnki.com.cn/Article/CDMD-10532-1012482435.htm [20] 何洪文, 黄贤广, 张亚明, 等.基于DIRECT算法的混合动力汽车参数优化研究[J].汽车工程学报, 2011, 1(2):35-41 http://d.old.wanfangdata.com.cn/Periodical/qcgcxb201101005He H W, Huang X G, Zhang Y M, et al. Study on parameters optimization of hybrid electric vehicle based on DIRECT algorithm[J]. Chinese Journal of Automotive Engineering, 2011, 1(2):35-41(in Chinese) http://d.old.wanfangdata.com.cn/Periodical/qcgcxb201101005