-
摘要:
在矩独立重要性分析过程中,重要性指标往往用于衡量结构系统输出不确定性向输入变量不确定性的逆向分配问题。假设输入参数的方差可以减缩一定比例因子,那么矩独立重要性指标可以定义为该缩减因子的函数。同时,假设输入参数的方差缩减因子为一随机变量,那么可以取矩独立重要性指标函数的均值定义一个新的平均矩独立重要性指标。由于使用Sobol方法计算平均矩独立重要性指标的模型需要循环抽样,计算量很高,故引入拒绝抽样(RS)方法,通过重复利用矩独立重要性分析中的一组输入输出样本,就可以额外计算得到矩独立指标函数和平均矩独立重要性指标,这大大节约了计算成本。本文所提指标函数及平均指标的有效性和RS方法的准确性、高效性通过数值和工程算例得以验证。
-
关键词:
- 重要性分析 /
- 矩独立重要性指标函数 /
- 平均矩独立重要性指标 /
- Sobol方法 /
- 拒绝抽样(RS)
Abstract:In the process of moment-independent importance analysis, importance index are always used to quantify the inverse allocation of structural system output uncertainty to input uncertainty.An assumption is given that the variance of a given factor can be reduced by future research, which leads to the development the moment-independent index function. The moment-independent importance index function provides an index for a given factor as a function of the amount of variance of that factor can be reduced. Meanwhile, by assuming the reduction amount of a particular factor variance as a random variable, the average moment-independent importance index is defined by taking average of the moment-independent importance index function. Estimating the average moment-independent importance index involves a large amount of computation using Sobol's method, and thus rejection sampling (RS) method is introduced here with the generated samples used in Sobol's method. Consequently, RS can use the samples generated during Sobol's method to accurately estimate the moment-independent importance index function and the average moment-independent importance index without any further model evaluation, which greatly reduces the computational cost. Numerical and engineering examples are demonstrated to show the effectiveness of the proposed measures and the accuracy and availability of the RS method.
-
表 1 屋架结构模型中输入变量的分布参数
Table 1. Distribution parameters of input variables of roof truss structure model
分布参数 均值 变异系数 q 20 000 N/m 0.07 l 12 m 0.01 AS 9.82×10-4 m2 0.06 AC 0.04 m2 0.12 ES 1×1011 N/m2 0.06 EC 2×1010 N/m2 0.06 表 2 十杆桁架结构模型中输入变量的分布参数
Table 2. Distribution parameters of input variables of ten-bar truss structure model
分布参数 均值 变异系数 L 1 m 0.05 E 100 GPa 0.05 P1 80 kN 0.05 P2 10 kN 0.05 P3 10 kN 0.05 Ai 0.001 m2 0.15 -
[1] 吕震宙, 宋述芳, 李洪双, 等.结构机构可靠性及可靠性灵敏度分析[M].北京:科学出版社, 2009:2-7.LU Z Z, SONG S F, LI H S, et al.Reliability and reliability analysis for structure and mechanics[M].Beijing:Science Press, 2009:2-7(in Chinese). [2] SALTELLI A, TARANTOLA S, CAMPLONGO F, et al.Sensitivity analysis in practice:A guide to assessing scientific models[M].New York:Wiley, 2004:3-10. [3] WANG W, ZHOU C, GAO H, et al.Application of non-probabilistic sensitivity analysis in the optimization of aeronautical hydraulic pipelines[J].Structural and Multidisciplinary Optimization, 2018, 57(6):2177-2191. doi: 10.1007/s00158-017-1848-4 [4] BORGONOVO E.Measuring uncertainty importance:Investigation and comparison of alternative approaches[J].Risk Analysis, 2006, 26(5):1349-1361. doi: 10.1111/risk.2006.26.issue-5 [5] HELTON J C, DAVIS F J.Latin hypercube sampling and the propagation of uncertainty in analyses of complex systems[J].Reliability Engineering & System Safety, 2003, 81:23-69. https://www.sciencedirect.com/science/article/pii/S0951832003000589 [6] SALTELLI A, MARIVOET J.Non-parametric statistics in sensitivity analysis for model output:A comparison of selected techniques[J].Reliability Engineering & System Safety, 1990, 28:229-253. https://www.sciencedirect.com/science/article/pii/095183209090065U [7] MORRIS M.Factorial sampling plans for preliminary computational experiments[J].Technometrics, 1991, 33(2):161-174. doi: 10.1080/00401706.1991.10484804 [8] LAMBONIA M, IOOSS B, POPELIN A L et al.Derivative-based global sensitivity measures:General links with Sobol' indices and numerical tests[J].Mathematics and Computers in Simulation, 2013, 87:45-54. doi: 10.1016/j.matcom.2013.02.002 [9] SOBOL I M.Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates[J].Mathmatics and Computers in Simulation, 2001, 55(1):221-280. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=84bf6510004c1dd0e68e68e1ce4e738f [10] CASTILLO E, MINGUEZ R, CASTILLO C.Sensitivity analysis in optimization and reliability problems[J].Reliability Engineering & System Safety, 2008, 93(12):1788-1800. https://www.sciencedirect.com/science/article/pii/S0951832008000999 [11] BORGONOVO E.A new uncertainty importance measure[J].Reliability Engineering & System Safety, 2007, 92(6):771-784. http://d.old.wanfangdata.com.cn/OAPaper/oai_arXiv.org_0911.1455 [12] CHUN M H, HAN S J, TAK N I.An uncertainty importance measure using a distance metric for the change in a cumulative distribution function[J].Reliability Engineering & System Safety, 2000, 70(3):313-321. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=e0ae54c4ba810fb7a7bac1b6ea0cc94d [13] ARCHER K J, KIMES R V.Empirical characterization of random forest variable importance measures[J].Computational Statistics & Data Analysis, 2008, 52(4):2249-2260. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=29b5837a80a90735efbbede0f85c80ac [14] CUI L J, LU Z Z, ZHAO X P.Moment-independent importance measure of basic random variable and its probability density evolution solution[J].Science China Technological Sciences, 2010, 53:1138-1145. doi: 10.1007/s11431-009-0386-8 [15] 王文选, 高行山, 周长聪.基于点估计的矩独立重要性测度分析方法[J].机械工程学报, 2017, 53(8):16-24. http://d.old.wanfangdata.com.cn/Periodical/jxgcxb201708003WANG W X, GAO H S, ZHOU C C.The moment-independent importance measure analysis method based to point-estimate[J].Chinese Journal of Mechanical Engineering, 2017, 53(8):16-24(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/jxgcxb201708003 [16] ZHOU C, LU Z, LI W.Sparse grid integration based solutions for moment-independent importance measures[J].Probabilistic Engineering Mechanics, 2015, 39:46-55. doi: 10.1016/j.probengmech.2014.12.002 [17] LI L Y, LU Z Z.Moment-independent importance measure of input variable and its state-dependent parameter solution[J].Structrual Safety, 2012, 38:40-47. doi: 10.1016/j.strusafe.2012.04.001 [18] WEI P F, LU Z Z, HAO W R, et al.Efficient sampling methods for global reliability sensitivity analysis[J].Computer Physics Communications, 2012, 183:1728-1743. doi: 10.1016/j.cpc.2012.03.014 [19] ALLAIRE D L, WILLCOX K E.A variance-based sensitivity function for factor prioritization[J].Reliability Engineering & System Safety, 2012, 107:107-114. https://www.sciencedirect.com/science/article/pii/S0951832011001712 [20] DEGROOT M, SCHERVISH M.Probability and statistics[M].3rd ed.Boston:Addison-Wesley, 2002:35-40. [21] WANG X.Improving the rejection sampling method in quasi-Monte Carlo methods[J].Journal of Computational & Applied Mathematics, 2000, 114:231-246. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=d8b5e8216c81012df7a5c9431acc67ec [22] STORLIE C, SWILER L, HELTON J, et al.Implementation and evaluation of nonparametric regression procedures for sensitivity analysis of computationally demanding models[J].Reliability Engineering & System Safety, 2009, 94:1735-1763. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=841bafee53d4545efaf6661e8aa12720 [23] LU Z Z, LIU C L, YUE Z F.Probabilistic safe analysis of the working life of a powder-metallurgy turbine disk[J].Materials Science and Engineering:A, 2005, 395:153-159. doi: 10.1016/j.msea.2004.12.008