張克偉

摘 要：1stOpt的計算核心是基于一種“通用全局優(yōu)化算法”，該算法之最大特點是克服了使用迭代法必須給出合適初始值的難題。本文首先采用《超聲回彈綜合法檢測混凝土強度技術(shù)規(guī)程》（CECS02：88）提供的三十組超聲回彈綜合法數(shù)據(jù)，利用1stOpt提供的包維爾法+通用全局優(yōu)化算法等多種算法，建立超聲回彈綜合法三參數(shù)和五參數(shù)冪函數(shù)曲面的的擬合方程，并在“合肥某科學研究院行政服務中心綜合樓”在建工程混凝土強度檢測的基礎(chǔ)上對比分析各算法的優(yōu)劣，為各方法在混凝土檢測中的應用提供參考依據(jù)。

關(guān)鍵詞：1stOpt；超聲回彈綜合法；三參數(shù)；五參數(shù)；混凝土強度

1.基于麥夸特法的數(shù)學模型

1.1建立綜合法三參數(shù)模型

采用《規(guī)程》提供的實驗數(shù)據(jù)，編寫代碼本如下：

Parameters a，b，c；

Variable x1，x2，y；

Function y=a*x1^b*x2^c；

Data；

4.8 31.0 25.3

4.7 30.8 26.0

4.6 30.5 27.1

（以上每行數(shù)據(jù)從左到右依次為為構(gòu)件混凝土中聲速值Vi、回彈值Ri和混凝土抗壓強度值fi）

優(yōu)化算法采用麥夸特法+通用全局優(yōu)化法，收斂判斷指標為1.00E-10，最大迭代次數(shù)為1000次。麥夸特法+通用全局優(yōu)化法建立的綜合法三參數(shù)冪函數(shù)曲面可采用如下方程：

-結(jié)構(gòu)或構(gòu)件第i個測區(qū)混凝土抗壓強度換算值，MPa；

-修正后測區(qū)回彈代表值；

-測區(qū)聲速代表值。

1.2三參數(shù)模型評價

由三參數(shù)冪函數(shù)曲面方程檢測混凝土抗壓強度的主要計算公式如下：

構(gòu)件混凝土強度推定值主要取決于 和 ，進而可建立該方法的統(tǒng)一數(shù)學模型如下：

由模型 可看出混凝土抗壓強度推定值 的不確定度來源于修正后測區(qū)混凝土中聲速代表值 及修正后測區(qū)回彈代表值 。由于兩種測強曲線的形式相同而參數(shù)不同，故不確定度分析中，回彈值和聲速值的靈敏系數(shù)不同，而其他因素引起的不確定度均相同。

（1）聲速代表值的靈敏系數(shù)

在所研究點聲速代表值的靈敏系數(shù)即為改點強度換算值對 的偏導數(shù)，即 ，回彈代表值取平均值 ，聲速代表值取 。

（2）回彈代表值的靈敏系數(shù)

求出回彈代表值的靈敏系數(shù) ?；貜棿碇等∑骄岛吐曀俅碇狄簿∑淦骄担M一步計算可得各構(gòu)件的靈敏系數(shù)。合成標準不確定度uc（y）確定后，與一個包含因子k相乘，可以期望在[y-U，y+U]包含了測量結(jié)果可能值的較大部分，這里k=2。

1.3建立綜合法五參數(shù)模型

麥夸特法+通用全局優(yōu)化法建立綜合法五參數(shù)冪函數(shù)曲面可采用如下方程：

2.結(jié)論

本節(jié)采用《超聲回彈綜合法檢測混凝土強度技術(shù)規(guī)程》（CECS02：88）提供的三十組超聲回彈綜合法數(shù)據(jù)，利用1stOpt提供的麥夸特法+通用全局優(yōu)化算法，建立超聲回彈綜合法三參數(shù)冪函數(shù)曲面和五參數(shù)冪函數(shù)曲面的的擬合方程，并在“合肥某科學研究院行政服務中心綜合樓”在建工程混凝土強度檢測的基礎(chǔ)上對比分析各算法的優(yōu)劣，總體來看，五參數(shù)冪函數(shù)曲面比三參數(shù)曲面較準，不確定度較小。以上結(jié)論可為各方法在混凝土檢測中的應用提供依據(jù)。

參考文獻

[1]JJF1059-1999，測量不確定度評定與表示[S].

[2]CECS 02：2005，超聲回彈綜合法檢測混凝土強度技術(shù)規(guī)程[S].

[3]GB/T4883-2008，數(shù)據(jù)的統(tǒng)計處理和解釋正態(tài)樣本離群值的判斷和處理[S].

[4]呂列民等.鉆芯法檢測結(jié)構(gòu)混凝土合格性評定分析[J].水利水電技術(shù)，2011，42（8）：21-23.

[5]林洪樺.測量誤差與不確定度評估[M].北京：機械工業(yè)出版社，2009.8.

（作者單位：安徽省建筑設(shè)計研究總院股份有限公司）