羅 杰，賈 靜

（河北省水利水電勘測設(shè)計(jì)研究院集團(tuán)有限公司，天津 300220）

目前，在巖土工程項(xiàng)目檢測試驗(yàn)中，盡管現(xiàn)行規(guī)程規(guī)范條文大力提倡通過三軸壓縮試驗(yàn)獲得土的抗剪強(qiáng)度指標(biāo)，但由于直接剪切試驗(yàn)（以下簡稱為直剪試驗(yàn)）的方法簡單、高效和大量經(jīng)驗(yàn)積累，其仍是目前巖土工程中獲得抗剪強(qiáng)度指標(biāo)的重要手段。

在直剪試驗(yàn)數(shù)據(jù)處理方法上，土工試驗(yàn)規(guī)程［1-2］或土工試驗(yàn)方法標(biāo)準(zhǔn)［3］中均未指定試驗(yàn)數(shù)據(jù)線性關(guān)系式的擬合方法， 只在數(shù)據(jù)處理部分的條文中規(guī)定采用圖解法［1-3］確定試樣的抗剪強(qiáng)度指標(biāo)。但隨著計(jì)算機(jī)的廣泛應(yīng)用，越來越多的數(shù)據(jù)處理通過計(jì)算機(jī)來完成。 本文通過實(shí)際工程項(xiàng)目的直剪試驗(yàn)數(shù)據(jù)為例， 分別采用最小二乘法、梯度降減法擬合試驗(yàn)數(shù)據(jù)的線性關(guān)系式，根據(jù)各自關(guān)系式確定抗剪強(qiáng)度指標(biāo)C，φ 值。同時(shí),以實(shí)際試驗(yàn)數(shù)據(jù)為標(biāo)準(zhǔn), 分析采用不同方法擬合的線性關(guān)系式可靠性； 以最小二乘法確定的抗剪強(qiáng)度指標(biāo)為參照值， 比較其他方法確定抗剪強(qiáng)度指標(biāo)與其之間的差異。

1 試樣基本物理性質(zhì)指標(biāo)

試樣的基本物理性質(zhì)指標(biāo)如表1。

表1 試樣基本物理性質(zhì)指標(biāo)

從表1 試樣的基本物理性質(zhì)指標(biāo)可看出：3 個(gè)試樣均屬黏性土。

2 直剪試驗(yàn)數(shù)據(jù)

為更好地比較分析不同線性擬合方法對直剪試驗(yàn)數(shù)據(jù)擬合結(jié)果的影響， 在實(shí)際工程工程項(xiàng)目的土工試驗(yàn)中，選取離散性差別較大的3 組試驗(yàn)數(shù)據(jù)，如表2。

表2 直剪切試驗(yàn)數(shù)據(jù)

從表2 可知：在3 組數(shù)據(jù)中，第1 組數(shù)據(jù)的離散性最大，第2 組數(shù)據(jù)的離散性次之，第3 組數(shù)據(jù)的離散性最小。

3 最小二乘法［4-5］

為確定試樣的抗剪強(qiáng)度指標(biāo)（C，φ），需根據(jù)不同垂直壓力下的抗剪強(qiáng)度繪制關(guān)系曲線或擬合關(guān)系式，一般情況下關(guān)系曲線為直線，關(guān)系式為線性表達(dá)式。 通過線性關(guān)系式中的參數(shù)達(dá)到確定試樣抗剪強(qiáng)度指標(biāo)（C，φ）的目的。目前，為提高工作效率和保證成果質(zhì)量， 基本采用辦公軟件或?qū)I(yè)軟件進(jìn)行數(shù)據(jù)處理。

由于最小二乘法簡便， 是實(shí)際工作中較常用的線性關(guān)系式擬合方法。最小二乘法基本原理：采用偏差平方和最小的原則擬合線性關(guān)系式。如設(shè)：線性表達(dá)式的型式為y=a+bx，針對表2 數(shù)據(jù)，其矩陣方程組為：

解方程組得：

依據(jù)上述兩式和表2 數(shù)據(jù),分別計(jì)算各試樣試驗(yàn)數(shù)據(jù)擬合的現(xiàn)行表達(dá)式參數(shù)a，b 值。計(jì)算過程及計(jì)算結(jié)果如表3。各試樣的抗剪強(qiáng)度與垂直壓力關(guān)系曲線如圖1。

圖1 抗剪強(qiáng)度與垂直壓力的關(guān)系曲線

表3 各試樣試驗(yàn)數(shù)據(jù)擬合的現(xiàn)行表達(dá)式參數(shù)計(jì)算結(jié)果

從圖1 可看到：所顯示公式中的a，b 值與表3 中的計(jì)算值完全一致，因此，借助Excel 選取散點(diǎn)圖，其數(shù)據(jù)的線性趨勢線擬合方法采用的是最小二乘法。

4 梯度降減法［4,6］

梯度降減法需進(jìn)行迭代計(jì)算，有時(shí)還需進(jìn)行優(yōu)化后才能確定最后結(jié)果，因此在電算化沒有廣泛普及的年代，其使用受到一定限制。但今天，廣泛普及的計(jì)算機(jī)應(yīng)用已完全化解此問題。

本文應(yīng)用梯度降減法的基本思路：在坐標(biāo)圖中，沿著代價(jià)函數(shù)在某點(diǎn)處切線的降低方向， 通過各待求變量變化的適當(dāng)步長，確定代價(jià)函數(shù)的下一點(diǎn)，并對應(yīng)1 個(gè)預(yù)測函數(shù)。 然后以此點(diǎn)處代價(jià)函數(shù)切線的降低方向，通過相同的變量變化步長，確定代價(jià)函數(shù)的另一個(gè)下降點(diǎn)，同樣對應(yīng)1 個(gè)預(yù)測函數(shù)。如此不斷地找下去，直至找到滿足要求的代價(jià)函數(shù)最小值或約定的變化值所對應(yīng)的點(diǎn)，此時(shí)對應(yīng)的預(yù)測函數(shù)是所需的最終結(jié)果。

設(shè)擬合的線性表達(dá)式（預(yù)測函數(shù)）為：y=a+bx

為評判獲得的預(yù)測函數(shù)是否滿足要求，引入均方差代價(jià)函數(shù)。

設(shè)均方差代價(jià)函數(shù)（或稱平方差代價(jià)函數(shù)、損失函數(shù)）為：

式中 m 為數(shù)據(jù)組數(shù)（樣本數(shù)），本文實(shí)例中m=4；1/2 為常量， 為在求梯度時(shí)能把二次方乘數(shù)加以抵消，這樣方便計(jì)算且不會(huì)影響計(jì)算結(jié)果；yi為試驗(yàn)實(shí)測數(shù)據(jù)值；ya，b為預(yù)測函數(shù)。根據(jù)a，b 值得到預(yù)測的y值。即：ya,b(xi)=a+bxiaj，bj值的迭代計(jì)算公式如式（8）～式（9）：

根據(jù)表2 中試樣編號為003004 的試驗(yàn)數(shù)據(jù)，采用梯度降減法進(jìn)行迭代計(jì)算的過程及計(jì)算結(jié)果如表4。 通過優(yōu)化， 在計(jì)算過程中aj，bj的初始值選取為：a0=20，b0=0.4；步長值選取為：αa=0.3，αb=0.00002。

表4 直剪試驗(yàn)數(shù)據(jù)線性擬合表達(dá)式參數(shù)a，b 值迭代計(jì)算

從表4 中aj，bj的迭代計(jì)算過程和計(jì)算結(jié)果可看出： 當(dāng)以均方差代價(jià)函數(shù)的最小值或約定的變化值為判斷標(biāo)準(zhǔn)時(shí)，a=76.717404，b=0.368912，修約后，a=76.72，b=0.3689。 預(yù) 測 函 數(shù) 的 最 終 關(guān) 系 式 為：y=0.3689x+76.72。 與上述采用最小二乘法獲得的關(guān)系式非常接近。同時(shí)，從表4 還可看出：如果持續(xù)增加計(jì)算步驟， 那么最終的a，b 值與采用最小二乘法獲得的值完全相同。

同理，根據(jù)表2 中試樣編號為003012 的試驗(yàn)數(shù)據(jù)， 通過優(yōu)化， 在計(jì)算過程中aj，bj的初始值a0=20，b0=0.4。步長αa=0.3，αb=0.00002。當(dāng)以均方差代價(jià)函數(shù)的最小值為判斷標(biāo)準(zhǔn)時(shí)，a=42.9775，b=0.341277，修約后，a=42.98，b=0.3413。 預(yù)測函數(shù)的最終關(guān)系式為：y=0.3413x+42.98。

同理，根據(jù)表2 中試樣編號為003032 的試驗(yàn)數(shù)據(jù)， 通過優(yōu)化， 在計(jì)算過程中aj，bj的初始值a0=20，b0=0.4。步長αa=0.3，αbt=0.00002。當(dāng)以均方差代價(jià)函數(shù)的最小值為判斷標(biāo)準(zhǔn)時(shí)，a =27.447728，b =0.399208，修約后，a=27.44，b=0.3992。預(yù)測函數(shù)的最終關(guān)系式為：y=0.3992x+27.44。

5 關(guān)系式可靠性分析

為驗(yàn)證上述關(guān)系式的可靠性， 以試驗(yàn)實(shí)測數(shù)據(jù)為基準(zhǔn)值， 比較采用各方法擬合關(guān)系式的計(jì)算值與基準(zhǔn)值的偏差情況，如表5、表6。從表5、表6 可看出： 采用不同方法擬合的線性關(guān)系式并分別計(jì)算在不同垂直壓力下的抗剪強(qiáng)度，計(jì)算結(jié)果略有不同，相差不大，但與實(shí)測數(shù)據(jù)均存在一定偏差，其偏差絕對值隨各組試驗(yàn)數(shù)據(jù)離散性增大而增大。

表5 各垂直壓力下抗剪強(qiáng)度的實(shí)測值和關(guān)系式計(jì)算值

表6 采用不同計(jì)算方法線性擬合關(guān)系式計(jì)算的抗剪強(qiáng)度與實(shí)測值偏差情況

6 不同方法確定抗剪強(qiáng)度指標(biāo)的差異

最小二乘法：為方便求解，需關(guān)心誤差函數(shù)在迭代值處的局部性質(zhì)， 而不用考慮在迭代值處的全局性質(zhì)，這種方法機(jī)器學(xué)習(xí)領(lǐng)域使用的非常廣泛。

梯度降減法： 是用于找到可微函數(shù)的局部最小值的一階迭代優(yōu)化算法。 為了使用梯度下降找到函數(shù)的局部最小值，采取與該數(shù)在當(dāng)前點(diǎn)的梯度（或近似梯度）的負(fù)值成比例的步驟。梯度下降法每次都以梯度的反方向下降，所以，有可能會(huì)走出鋸齒路線，從而增加迭代次數(shù)。

不論什么方法，只要選好初始參數(shù)，就可減少迭代次數(shù)，從而減少計(jì)算量。在實(shí)際使用中，可根據(jù)具體情況決定用哪種方法優(yōu)化求解。

不同線性擬合方法確定的試樣抗剪強(qiáng)度指標(biāo)如表7。

表7 不同線性方法確定的抗剪強(qiáng)度指標(biāo)

從表7 可看出： 兩種方法確定的抗剪強(qiáng)度指標(biāo)幾乎沒有差別。

7 結(jié)語

（1） 對直剪試驗(yàn)數(shù)據(jù)線性的擬合方法不是單一的。 擬合方法不同， 其結(jié)果也可能不同。為了保證數(shù)據(jù)的連續(xù)性，盡可能采用同一種線性擬合方法。

（2） 在兩種線性擬合方法中，最小二乘法最簡單、最方便。一組試驗(yàn)數(shù)據(jù)采用某種方法擬合線性關(guān)系式后，要驗(yàn)證關(guān)系式的可靠性，當(dāng)發(fā)現(xiàn)綜合誤差較大且不滿足實(shí)際要求時(shí)，應(yīng)考慮采用其他線性擬合方法。