超声回弹综合法不同模型的对比分析0824.docVIP

超声回弹综合法不同模型的对比分析 引言 1stOpt是七维高科有限公司（7D-Soft High Technology Inc.）开发并拥有完全自主知识产权的一套数学优化分析综合工具软件包，在非线性回归，曲线拟合，非线性复杂模型参数估算求解，线性/非线性规划等领域功能强大。其计算核心是基于一种“通用全局优化算法”(Universal GlobalOptimization- UGO)，该算法之最大特点是克服了使用迭代法必须给出合适初始值的难题，即由 1stOpt随机给出参数初始值，通过其独特的全局优化算法，最终找出最优解。以非线性回归为例，该领域常用的软件工具包诸如 Matlab, OriginPro, SAS, SPSS,DataFit, GraphPad 等，均需提供适当的参数初始值以便计算能够收敛并找到最优解，如果设定的参数初始值不当则计算难以收敛，其结果是无法求得正确结果。而在实际应用当中，对大多数用户来说，给出恰当的初始值是件相当困难的事，而1stOpt凭借其超强的寻优、容错能力，在大多数情况下（大于 90%)，从任一随机初始值开始，都能求得正确结果。 本节首先采用《超声回弹综合法检测混凝土强度技术规程》（CECS02:88）提供的三十组超声回弹综合法数据，利用1stOpt提供的包维尔法(Powell Optimization-PO)+通用全局优化算法(Universal GlobalOptimization - UGO)等多种算法，分别建立超声回弹综合法（以下简称综合法）三参数幂函数曲面、五参数幂函数曲面的的拟合方程，并在“安徽省·水利部淮委水利科学研究院行政服务中心综合楼”在建工程混凝土强度检测的基础上对比分析各算法的优劣，为各方法在混凝土检测中的应用提供参考依据。 1.基于麦夸特法的数学模型 1.1建立综合法三参数模型 表1是引自《超声回弹综合法检测混凝土强度技术规程》（CECS02:88）提供的超声、回弹、抗压强度实验数据。 序号i 声速Vi 回弹值Ri 强度值fi 序号i 声速Vi 回弹值Ri 强度值fi 1 4.80 31.0 25.3 16 4.41 28.2 19.0 2 4.75 30.8 26.0 17 4.37 26.1 18.7 3 4.66 30.5 27.1 18 4.50 27.0 19.6 4 4.87 38.6 39.0 19 4.63 31.6 23.8 5 4.85 36.6 38.8 20 4.65 27.2 20.8 6 4.91 38.0 40.7 21 4.58 30.1 23.9 7 4.07 20.7 10.0 22 4.62 30.5 24.9 8 4.08 18.8 10.0 23 4.70 30.7 25.5 9 4.23 17.2 10.2 24 4.60 29.9 25.0 10 4.40 20.9 14.8 25 4.77 38.6 41.8 11 4.56 21.5 15.9 26 4.79 40.4 47.9 12 4.45 20.0 14.6 27 4.75 36.8 39.0 13 4.28 24.0 13.0 28 4.78 41.0 46.8 14 4.20 25.2 13.3 29 4.70 35.0 36.3 15 4.25 24.6 13.5 30 4.75 36.3 36.7 利用表1数据编写代码本如下： Parameters a,b,c; Variable x1,x2,y; Function y=a*x1^b*x2^c; Data; 4.80 31.0 25.3 4.75 30.8 26.0 4.66 30.5 27.1 。。。。。。 4.78 41.0 46.8 4.70 35.0 36.3 4.75 36.3 36.7 （以上每行数据从左到右依次为为构件混凝土中声速值Vi、回弹值Ri和混凝土抗压强度值fi） 优化算法 收敛判断指标 最大迭代数 麦夸特法+通用全局优化法 1.00E-10 1000 达到收敛判断指标是计算终止，计算结果如下。 迭代数 计算用时 相关系数R 22 4秒484微秒 0.9865 a b c 0.00526 1.481 1.820 由以上计算结果可以认为，麦夸特法+通用全局优化法建立的综合法三参数幂函数曲面可采用如下方程： 式中： ——结构或构件第i个测区混凝土抗压强度换算值，MPa； ——修正后测区回弹代表值； ——测区声速代表值。 1.2三参数模型评价 由麦夸特法+通用全局优化法建立的综合法三参数幂函数曲面方程检测混凝土抗压强度的主要计算公式如下： 式中： —— 结构或构件第i个测区混凝土抗