摘要：填方路基的不均勻沉降或過大的工后沉降，都會引起路面出現(xiàn)破損甚至破壞，致使公路的安全系數(shù)減小、通行能力降低和使用壽命縮短，因此，相關(guān)規(guī)范作出了必須進行路基沉降監(jiān)測的規(guī)定，以保障道路始終處于穩(wěn)定的良性的狀態(tài)，最大限度地發(fā)揮效益。路基沉降，是運行車輛的動荷載、道路本身的靜荷載與含水等因素持續(xù)、同時、再三作用的結(jié)果，隨時間發(fā)生、發(fā)育、發(fā)展之特征是明顯的，因此，監(jiān)測須定期、反復(fù)進行，并且每一次完整的監(jiān)測都不單單是外業(yè)測量，還包括沉降分析、評價和變化趨勢預(yù)測預(yù)報等項，其數(shù)據(jù)處理工作勢必繁復(fù)、大量。若采用計算機編程方法完成上述任務(wù)，一則能夠極大減輕計算量度；二則能夠讓處理程序程式化，從而嚴謹化；三則能夠?qū)崿F(xiàn)數(shù)據(jù)分析圖像化、檢評結(jié)論可視化；四則能夠讓數(shù)學(xué)功底稍差、計算機知識略遜的人員正常實施作業(yè)。

關(guān)鍵詞：路基；沉降；VB語言；點擊；界面

中圖分類號：TU433 文獻標志碼：Ａ文章編號：１６７３－２９１Ｘ（２０11）17－０312－02

引言

山區(qū)道路建設(shè)中常見的高填方路基，在行車動荷載的沖擊力、大量土石方靜荷載的垂直壓力、堤體含水量反復(fù)發(fā)生變化及不可能盡善盡美的施工質(zhì)量等因素持續(xù)、同時、再三作用下，易產(chǎn)生過大的工后沉降或不均勻沉降，引起路面出現(xiàn)破損甚至破壞，導(dǎo)致道路安全系數(shù)減小、通行能力降低和使用壽命縮短。因此，《公路路基設(shè)計規(guī)范》、《公路路基施工技術(shù)規(guī)范》、《公路工程質(zhì)量檢驗評定標準》都明確規(guī)定，無論是施工期間還是運營階段，都必須開展路基沉降監(jiān)測，以控制沉降程度與變化趨向，采取有效措施，確保路基始終處于穩(wěn)定的良性的狀態(tài)，支持道路整體，最大限度地發(fā)揮其經(jīng)濟效益和社會效益。由于路基沉降具有隨時間發(fā)生、發(fā)育、發(fā)展之顯著特征，監(jiān)測需定期、反復(fù)進行，并且每一次完整的監(jiān)測都不單單是外業(yè)測量，還包括沉降分析、評價和變化趨勢預(yù)測預(yù)報等項，其數(shù)據(jù)處理工作勢必繁復(fù)、大量，此時，若采用計算機編程方法完成上述任務(wù)，一則能夠極大減輕計算量度；二則能夠讓處理程序程式化，從而嚴謹化；三則能夠?qū)崿F(xiàn)數(shù)據(jù)分析圖像化，檢評結(jié)論可視化；四則能夠讓數(shù)學(xué)功底稍差、計算機知識略遜的人員正常實施作業(yè)。

既然路基沉降是隨時間變化的，所以沉降變形分析須考慮兩個因素：時間和沉降量。又由于“一門科學(xué)只有成功運用了數(shù)學(xué)，才算達到了真正完善的地步?！保R克思語）沉降分析理論依據(jù)數(shù)理統(tǒng)計方法。為構(gòu)筑以時間為自變量、沉降量為函數(shù)的沉降變化數(shù)模，首先繪制沉降數(shù)據(jù)曲線圖（沉降監(jiān)測值為縱標、時間為橫標），閱圖判斷模型，以模型的可能類型為藍本進行方程擬合（若是線性模型，求出回歸系數(shù)a、b嵌入模型藍本；若是非線性模型，先將其線性化，再求回歸系數(shù)a、b，然后嵌入模型藍本），得到初步模型，然后通過統(tǒng)計檢驗證明模型正確性（數(shù)模顯著性檢驗方法有很多，但相關(guān)性檢驗因其計算簡便、結(jié)論可信，成為最常用的檢驗方式，本文亦用此檢驗法），檢驗若證實模型成立，說明所擬方程就是沉降實況的正確表達式。運用回歸方程，代入時間變量，可求出沉降回歸分析（函數(shù)）值，從而從理論上掌握沉降變化。將同一時間的沉降監(jiān)測值與回歸分析值比較，其差值（殘差）愈小，說明實際沉降愈趨于穩(wěn)定；若其差較大，表明沉降正在加深，須加密監(jiān)測并查證原因，采取恰當補救措施控制局面；若差值突然大增，并伴監(jiān)測值互差變化較大，表明路基異常，是發(fā)生大沉陷的先兆，須及時搶救消除事故隱患。

須要說明的兩個問題是，編程所依據(jù)的數(shù)學(xué)原理公式，請參見參考文獻[1]，為精煉篇幅，本文未一一列出。若僅擬合形成一個初步模型，只須驗證其具有足夠的回歸精度即可；如試擬合形成的不止一個初步模型，則在滿足回歸精度的模型中，擇回歸精度最高者保留之。

一、 源數(shù)據(jù)

我國南方某高速公路試驗路段，土方填高28m，從路基填筑完工開始，選擇K27+700斷面上的一點，進行了歷時13個多月的沉降監(jiān)測，數(shù)據(jù)見表1（如果以每天為時間單位的話，橫坐標太長，為了繪圖的方便，以60天為時間單位，即0.45 60天=27天的時候進行了第一次測量，13個多月共監(jiān)測14次）。

二、運行及程序代碼

1. 運行

進入數(shù)據(jù)處理界面→點擊“導(dǎo)入數(shù)據(jù)”，選擇K27+700斷面右沉降數(shù)據(jù)的txt文件，輸入監(jiān)測數(shù)據(jù)→點擊“繪制折線圖”得沉降數(shù)據(jù)的折線圖→點擊“查看圖形”得五種模型的曲線→選擇指數(shù)模型點擊“建立模型” ，得擬合模型→點擊“模型檢驗” ，得相關(guān)系數(shù)、自由度→點擊“保存結(jié)果”，輸出回歸分析值、殘差等（分別見圖1、圖2）。

2. 代碼

Dim day1（1 To 1000） As Single

Dim data1（1 To 1000） As Single

Dim xiyi， x1， y1， s， c， e， b As Single

For i = 1 To Num

If data（i） = 0 Then data1（i） = data（i）

If data（i） <> 0 Then data1（i） = Log（-data（i））

day1（i） = 1 / day（i）

Next i

For i = 1 To Num

s = day1（i） * data1（i） + s

c = day1（i） ^ 2 + c

e = day1（i） + e

b = data1（i） + b

Next i

x1 = e / Num

y1 = b / Num

zb1 = Abs（（s - Num * x1 * y1） / （c - Num * （x1 ^ 2）））

zb0 = （y1 + zb1 * x1）

zb0 = Exp（zb0）

zb1 = Format（zb1， “0.00000”）

zb0 = Format（zb0， “0.00000”）

Text1.Text = “u=” -zb0 “×” “e^” “（” -zb1 “/t）”

“檢驗?zāi)Ｐ汀?/p>

Dim s1， s2， s3， s4， b0， b1， b2， r As Single

Dim day1（1 To 1000） As Single ‘Variant

Dim data1（1 To 1000） As Single ’Variant

For i = 1 To Num

If data（i） = 0 Then data1（i） = data（i）

If data（i） <> 0 Then data1（i） = Log（-data（i））

day1（i） = 1 / day（i）

Next i

For i = 1 To Num

s2 = data1（i） + s2

s1 = day1（i） + s1

Next i

For i = 1 To Num

s3 = （day1（i） - s1 / Num） * （data1（i） - s2 / Num） + s3

Next i

For i = 1 To Num

b1 = （day1（i） - s1 / Num） ^ 2 + b1

b2 = （data1（i） - s2 / Num） ^ 2 + b2

Next i

r = s3 / Sqr（b1 * b2）

Text3 = Format（r， “0.000”）

For i = 1 To Num

s4 = （Format（-zb0 * Exp（-zb1 / day（i）） - data（i）， “0.0”）） ^ 2 + s4

Next i

Text2 = Num - 2

Text4 = Format（s4 / （Num - 2）， “0.00”）

Text5 = Format（Sqr（s4 / （Num - 2））， “0.00”）

List2.Clear

List2.AddItem “時間” ““ “累計形變量（mm）” ““ “回歸分析值（mm）” ““ “V（mm）”

For i = 1 To Num

List2.AddItem Format（day（i）， “0.00”） “ “ Format（data（i）， “0.0”） “ “ Format（-zb0 * Exp（-zb1 / day（i））， “0.0”） ““ Format（-zb0 * Exp（-zb1 / day（i）） - data（i）， “0.0”）

Next i

三、說明與結(jié)論

1. 經(jīng)驗證明，路基沉降變化模型，指數(shù)與對數(shù)形式最多，由于篇幅限制，本文僅列出了程序的核心部分之一；指數(shù)模型回歸及其檢驗程序代碼，與對數(shù)模型代碼的編寫思路基本相似。

2. 圖2顯示，所建模型相關(guān)系數(shù)r＝-0.975，以14-2=12為引數(shù)、α＝0.01為保證率，查相關(guān)系數(shù)表，得相關(guān)系數(shù)臨界值r臨＝0.661，顯然|r|≥|r臨|， 證實本例整體回歸精度滿足要求，所擬沉降模型有效且顯著，回歸分析值均可靠；回歸擬合度分析算得擬合中誤差σ＝5.9 mm，小于容許誤差（2倍觀測誤差）10 mm，再一次確認了上述結(jié)論的正確性；殘差取值范圍[-9.9，7.7]，未超本例路基設(shè)計書規(guī)定的 20 mm 容許范圍。

3. 通過討論，可有把握地斷定，路基K27+700斷面的沉降整體無異常，而且隨時間的推移，前后監(jiān)測值互差愈來愈小，說明路基日趨穩(wěn)定，表明路基施工質(zhì)量情況良好。

4. 筆者之所以使用VB語言編程，除了它直觀易學(xué)、功能強大、開發(fā)空間廣闊外，還因為該語言是我國計算機等級考試的專門語言，在等級考試的有力推動下，VB成為我國使用普及率最高的語言之一； 應(yīng)用本程序開展路基沉降質(zhì)量檢評，操作人員毋須具備完備的數(shù)理統(tǒng)計理論知識，毋須具備計算機程序編寫能力，只須簡單的輸入和運行點擊即可。

參考文獻：

[1] 徐兮.最小二乘原理在隧道施工安全質(zhì)量預(yù)測預(yù)報中的應(yīng)用[J].數(shù)學(xué)的實踐與認識，2008，（6）：113-118.

[2] 白鳳山，么煥民，李春玲，沈繼紅.數(shù)學(xué)建模[M].哈爾濱：哈爾濱工業(yè)大學(xué)出版社，2002.

[3] 陳謙應(yīng)，蔣樹屏，柴賀軍，楊建國.山區(qū)公路路基穩(wěn)定理論與實踐[M].北京：人民交通出版社，2005.

[4] 何國斌.VisualBasic程序設(shè)計教程[M].重慶：重慶大學(xué)出版社，2006.

[5] 馮力.回歸分析方法原理及SPSS實際操作[M].北京：中國金融出版社，2004.

注：本文中所涉及到的圖表、注解、公式等內(nèi)容請以PDF格式閱讀原文