界限含水率程序設(shè)計

1、界限含水率試驗數(shù)據(jù)處理技術(shù)10界限含水率概述界限含水率試驗是土工試驗中最常見、應(yīng)用最廣的試驗之一。土的液塑限指標(biāo)是細(xì)粒土進(jìn)行分類和定名的最基本指標(biāo)，在巖土工程中，液塑限指標(biāo)的準(zhǔn)確性涉及土壤定名的正確性，同時影響土樣狀態(tài)的確定，進(jìn)而影響到土的承載能力的確定；在公路鐵路路基工程中，還直接影響細(xì)粒土的填料分組，所以正確地確定土壤的液塑限指標(biāo)對工程具有很重要的意義。21 界限含水率試驗數(shù)據(jù)處理原理 根據(jù)極限平衡理論，沿圓錐與土的接觸面的極限剪應(yīng)力（見圖1）為：Prlh3 （1）P圓錐重力錐角度數(shù)A圓錐與土的接觸面積，其計算方法如下： r為圓錐入土的半徑h入土深度l圓錐入土的斜邊長1 界限含水率試驗數(shù)據(jù)

2、處理原理 4一般的，由于為定值，令 則 （2） 對式(2)兩邊取對數(shù)得： （3）其中l(wèi)gKP為常數(shù)。由于重塑土的無側(cè)限抗壓強(qiáng)度與含水率存在雙對數(shù)關(guān)系，即 （4）1 界限含水率試驗數(shù)據(jù)處理原理 5將上式繪成雙對數(shù)線，是一條直線。對多種土進(jìn)行不同含水率和抗剪強(qiáng)度與對應(yīng)圓錐入土深度試驗，結(jié)果表明，理論曲線與試驗曲線在特性上相一致。經(jīng)整理得： （5） 其中a和b是由試驗擬合確定的常數(shù)，根據(jù)以上原理，界限含水率試驗應(yīng)采用雙對數(shù)坐標(biāo)系。1 界限含水率試驗數(shù)據(jù)處理原理 6目前界限含水率數(shù)據(jù)處理常用的方法是規(guī)范作圖法，此外還有直接線性最小二乘法和公式法。下面列舉某次試驗中的兩組數(shù)據(jù)(見表1)，并且用以上各種方

3、法計算界限含水率偏差。錐沉量h(mm)含水率(%)錐沉量h(mm)含水率(%)第一組數(shù)據(jù)第二組數(shù)據(jù)3.4022.595.6036.709.2030.279.4543.2815.9036.0315.6753.212 試驗數(shù)據(jù)處理方法及偏差計算7直接線性擬合函數(shù)簡單，在Excel表中處理方便，有些試驗人員在試驗時直接用線性擬合（采用了最小二乘法），且在Excel表中有時擬合度達(dá)到0.99以上。但根據(jù)以上原理，顯然直接線性擬合是不合理的，且不能估計塑限差距是否滿足要求，造成數(shù)據(jù)計算偏差較大。表2是直線擬合法的計算結(jié)果。2.1直接線性擬合法第一組數(shù)據(jù)第二組數(shù)據(jù)擬合方程y = 1.0697x + 19.

4、468擬合方程y = 1.636x + 27.644擬合度R2=0.9847擬合度0.9997p (%)21.61p (%)30.9210mmL(%)30.1710mmL(%)44.0017mmL(%)37.6517mmL(%)55.468目前界限含水率試驗數(shù)據(jù)處理一般采用TB10102-2004鐵路工程土工試驗規(guī)程法，以含水率為橫坐標(biāo)，圓錐下沉深度為縱坐標(biāo)，在雙對數(shù)坐標(biāo)紙上繪制如圖2所示的關(guān)系曲線，三點(diǎn)應(yīng)連成一條直線。當(dāng)三點(diǎn)不在一條直線上，則通過高含水率這一點(diǎn)與其余兩點(diǎn)連成兩條直線，在圓錐下沉深度為2mm處，可查得相應(yīng)的兩個含水率，當(dāng)這兩個含水率的差值小于2%時，應(yīng)以這兩點(diǎn)的含水率平均值與高

5、含水率的點(diǎn)連成一條直線。當(dāng)這兩個含水率之差值大于或等于2%時，則應(yīng)重做試驗。 2.2 規(guī)范作圖法第一組數(shù)據(jù)11010010100w(%)h(mm)9規(guī)范作圖法是土工試驗規(guī)范中所建議的方法，但該法有如下缺點(diǎn)：作圖在雙對數(shù)紙上進(jìn)行，由于對數(shù)坐標(biāo)不均勻分格，1、2之間和3、4之間的間距不相同，作圖時一般是按線性讀數(shù)，顯然會造成讀數(shù)偏差；讀數(shù)受人為因素影響較大；對2%的要求無法計算；當(dāng)塑性指數(shù)處于10和17附近時，直接影響了土的定名。所以在試驗中使用此法時求解不確定性較大。表3 作圖法計算第一組數(shù)據(jù)第二組數(shù)據(jù)塑限(%)18.50塑限(%)25.0010mmL(%)32.0010mmL(%)45.001

6、7mmL(%)37.517mmL(%)55.002.2 規(guī)范作圖法10設(shè)試驗參數(shù)如圖3所示。為便于計算，我們將lg(h)作為自變量，lg()作為因變量。AB的直線方程為同理AC的方程為 2.3 公式法(6)(8)11分別把h2mm代入式(6)(8)可得AB2和AC2，則由兩條直線引起的塑限偏差值為：p=AB2-AC2 (10) 當(dāng)p2%公式法的優(yōu)點(diǎn)是能明確p的大小，控制試驗是否需要重做。但該法以錐沉量最大者為準(zhǔn)確值，調(diào)整B、C兩點(diǎn)顯然也不盡合理，且不能保證三點(diǎn)按規(guī)范所擬合的誤差最小，特別是塑限偏差值接近2%時。第一組數(shù)據(jù)第二組數(shù)據(jù)塑限(%)18.93塑限(%)24.1410mmL(%)31.1

7、210mmL(%)44.7817mmL(%)36.7917mmL(%)54.352.3 公式法13公式法計算Dim Kab As Double, Kac As Double, Wab2 As Double, Wac2 As Double, Wpcha As DoubleDim ha As Double, hb As Double, hc As Double, wa As Double, wb As Double, wc As Doubleha = Val(Text48)hb = Val(Text47)hc = Val(Text46)wa = Val(Text51) / 100wb = Val(

8、Text50) / 100wc = Val(Text49) / 100If ha = 0 Or hb = 0 Or hc = 0 Or wa = 0 Or wb = 0 Or wc = 0 Or ha = hb Or hb = hc Or wa = wb Or wb = wc Then MsgBox 原始數(shù)據(jù)輸入有誤，請重輸 Exit SubEnd If輸入ha,hb,hcWa,wb,wc判別輸入數(shù)據(jù)有沒有空值14Kab = (Lg(wa) - Lg(wb) / (Lg(ha) - Lg(hb)Kac = (Lg(wa) - Lg(wc) / (Lg(ha) - Lg(hc) 兩段直線的斜率W

9、ab2 = 10 (Kab * Lg(2) + Lg(wa) - Kab * Lg(ha)Wac2 = 10 (Kac * Lg(2) + Lg(wa) - Kac * Lg(ha) 計算兩段直線所對應(yīng)的塑限Wpcha = (Wab2 - Wac2)塑限差值wText60(0) = Format(Wpcha * 100, #0.000)Dim tmp As Double, Wp As Double, Wl10 As Double, Wl17 As Doubletmp = (Lg(wa) - Lg(Wab2 + Wac2) / 2) / (Lg(ha) - Lg(2)tmp = tmp * (L

10、g(2) - Lg(ha) + Lg(wa)Wp = 10 tmpText60(1) = Format(Wp * 100, #0.000) 塑限求解Public Function Lg(ByVal x As Double) As DoubleLg = Log(x) / Log(10)End Function公式法的優(yōu)點(diǎn)，可以判定符合規(guī)范與否15tmp = (Lg(wa) - Lg(Wab2 + Wac2) / 2) / (Lg(ha) - Lg(2)tmp = tmp * (Lg(10) - Lg(ha) + Lg(wa)Wl10 = 10 tmpText60(2) = Format(Wl10

11、 * 100, #0.000) 10mm液限求解tmp = (Lg(wa) - Lg(Wab2 + Wac2) / 2) / (Lg(ha) - Lg(2)tmp = tmp * (Lg(17) - Lg(ha) + Lg(wa)Wl17 = 10 tmpText60(3) = Format(Wl17 * 100, #0.000) 17mm液限求解Text60(4) = Format(Wl10 - Wp) * 100, #0.000)Text60(5) = Format(Wl17 - Wp) * 100, #0.000) 塑性指數(shù)的求解特點(diǎn)：變量隨時可以定義，但在一個過程中不能兩次聲明一個變量

12、162.4先求對數(shù)再最小二乘擬合法 最小二乘法是德國數(shù)學(xué)家C.F.高斯在1794年解決行星軌道預(yù)測問題時首先提出的。線性擬合的最小二乘法基本原理為：設(shè)經(jīng)驗方程是y=f(x)，方程中含有一些待定系數(shù)an，給出真實值(xi,yi)|i=1,2,.n，將這些x，y值代入方程然后作差，可以描述誤差：yi- f (xi)，為了考慮整體的誤差，可以取平方和，之所以要平方是考慮到誤差可正可負(fù)直接相加可以相互抵消，所以記誤差為： 172.4先求對數(shù)再最小二乘擬合法 如果經(jīng)驗方程是線性的，形如y=ax+b，就是線性回歸。按上面的分析，誤差函數(shù)為： 于是得到關(guān)于a ，b的線性方程組： 要令e為最小，則e對a和b求

13、偏導(dǎo)數(shù)為零得：n為數(shù)據(jù)點(diǎn)數(shù)設(shè) (13)(14)(15)(16)(17)(18)182.4先求對數(shù)再最小二乘擬合法 則方程化簡為： (19) 解出上述方程組得：根據(jù)以上原理，利用Visual Basic編寫的最小二乘法線性擬合的關(guān)鍵函數(shù)程序段如下：192.4先求對數(shù)再最小二乘擬合法 Private Sub Nihe(x() As Double, y() As Double, a As Double, b As Double, e As Double)Dim n As IntegerDim aa As Double, bb As Double, cc As Double, dd As Double

14、, delta As Doubleaa = 0: bb = 0: cc = 0: dd = 0n = UBound(x) + 1For i = 0 To UBound(x) aa = aa + x(i) * x(i) bb = bb + x(i) cc = cc + x(i) * y(i) dd = dd + y(i)Nextdelta = aa * n - bb * bbIf Abs(delta) 0.0000000001 Then MsgBox errorElse a = (cc * n - bb * dd) / delta b = (aa * dd - cc * bb) / deltaE

15、nd Ife = 0For i = 0 To UBound(x)求解誤差 e = e + (y(i) - a * x(i) - b) 2Next iEnd Sub根據(jù)原理編寫如何調(diào)用過程（函數(shù)）20對上述函數(shù)調(diào)用Private Sub Command18_Click()Dim x(2) As Double, Y(2) As DoubleDim a As Double, b As Double, e As DoubleDim Wl10 As Double, Wl17 As Double, Wp2 As Double, Ip As Doublex(0) = Lg(Val(Text48)x(1) =

16、 Lg(Val(Text47)x(2) = Lg(Val(Text46)Y(0) = Lg(Val(Text51)Y(1) = Lg(Val(Text50)Y(2) = Lg(Val(Text49) 先求對數(shù)再擬合不求對數(shù)就是直接線性擬合法21過程調(diào)用續(xù) Nihe x(), Y(), a, b, e或者call nihe(x(),y(),a,b,e)過程調(diào)用的兩種寫法22反求液塑限續(xù)Wp2 = 10 (a * Lg(2) + b)Wl10 = 10 (a * Lg(10) + b)Wl17 = 10 (a * Lg(17) + b) 反求液塑限Ercheng(1) = Format(Wp2,

17、#0.000)Ercheng(2) = Format(Wl10, #0.000)Ercheng(3) = Format(Wl17, #0.000)Ercheng(4) = Format(Wl10 - Wp2, #0.000)Ercheng(5) = Format(Wl17 - Wp2, #0.000)Ercheng(0) = Format(1 - e, 擬合度 & #0.0000)23任意沉降對應(yīng)的含水量Dim Wren As DoubleIf Val(Text55) = 30 Thentext56.text= 出錯Exit SubElsetmp = 10 (a * Lg(Val(Text55

18、) + b)text56.text= Format(tmp, 0.#) 任一沉降對應(yīng)的含水率求解End IfEnd Sub242.4先求對數(shù)再最小二乘擬合法 注：x()數(shù)組是A、B、C三點(diǎn)的橫坐標(biāo)的對數(shù)，y()數(shù)組是相應(yīng)的縱坐標(biāo)的對數(shù)，a、b返回要求解的線性方程系數(shù)，e返回誤差的平方和。把A(lgA,lghA)，B(lgB,lghB)，C(lgC,lghC)三點(diǎn)坐標(biāo)組成的數(shù)組x(0)=lgA，x(1)=lgB，x(2)=lgC和y(0)=lghA，y(1)=lghB，x(2)=lghC兩個數(shù)組調(diào)用上述函數(shù)，即可得出擬合后的線性方程參數(shù)a，b。把y=lg(2)及y=lg(10)代入擬合公式可得塑限和液限。計算結(jié)果見表5。但該法也有自身的缺點(diǎn)，即擬合度很高時，也無法保證p2%的規(guī)范要求。25第一組數(shù)據(jù)第二組數(shù)據(jù)塑限(%)19.21塑限(%)25.1210m