劉闖，張銳，滕煜晟，任天锃，楊博

(1. 長沙理工大學(xué) 交通運輸工程學(xué)院，長沙 410114；2. 海南省交通運輸廳，海口 570204；3.重慶交通大學(xué) 土木工程學(xué)院，重慶 400074)

高液限土在中國南方地區(qū)分布廣泛，天然含水率高，水穩(wěn)性差，用于公路路基填筑時，通常用作下路堤填料，并在濕法重型擊實標準控制下壓實。由于高液限土下路堤頂面回彈模量(以下簡稱模量)較低，需在其上部進行剛度補償[1]，以達到路基頂面設(shè)計模量和驗收彎沉值的要求。高液限土下路堤頂面模量的快速準確測定對于高液限土路基施工質(zhì)量控制和剛度補償設(shè)計至關(guān)重要。

現(xiàn)場測定路基模量的方法主要有承載板法、貝克曼梁法、落錘彎沉儀FWD法和便攜式落錘彎沉儀PFWD法等。承載板法和貝克曼梁法耗時費力，所需測試人員多；FWD需專用車牽引，偏僻地方不宜到達，這3種方法多用于成型路基，難以用于路基施工質(zhì)量控制[2]，雖然PFWD測試簡便，但以往反算模量時僅取實測荷載和位移時程曲線的峰值，并假設(shè)路基為線彈性體，忽略了其他的時程數(shù)據(jù)和荷載位移的非線性[3-5]。對于沖擊荷載作用下變形響應(yīng)滯后且荷載位移非線性顯著的路基，基于線彈性模型和靜力分析方法反算的模量(簡稱常規(guī)反算模量)與承載板法直接實測的靜模量(簡稱實測模量)相差較大[6]。有學(xué)者通過現(xiàn)場試驗，建立了常規(guī)反算模量和實測模量的回歸關(guān)系[7-8]，但現(xiàn)場檢測工作量大，適用范圍有限。學(xué)者們曾對FWD和PFWD沖擊荷載作用下路基、路面非線性動力響應(yīng)開展研究，有學(xué)者將變形滯后歸因于彈性波傳遞時間，基于彈性動力學(xué)改進了常規(guī)反算方法，但僅將反算結(jié)果與基于線彈性模型的反算模量進行了對比[9-11]。也有學(xué)者將變形滯后歸因于材料的粘彈性，基于粘彈性理論，改進了反算方法[12-18]，并提高了模量反算精度[19-23]，但主要針對路面多層模量反算，且多采用多參數(shù)復(fù)雜理論模型，若用于高液限土路基模量反算，將增加計算難度，也可能導(dǎo)致結(jié)果不收斂。

為解決PFWD常規(guī)反算中因理論假設(shè)與實際情況不符導(dǎo)致的反算高液限土路基模量與實測模量相差較大的問題，筆者基于粘彈性模型，采用準靜態(tài)動力響應(yīng)分析方法，提出并驗證了既能反映高液限土路基非線性動力響應(yīng)特征，又能用動載數(shù)據(jù)反算靜力參數(shù)的高液限土路基模量反算新方法，編制了相應(yīng)的反算程序。同時，將其應(yīng)用于路基施工質(zhì)量控制與剛度補償設(shè)計。

1 高液限土路基模量現(xiàn)場試驗

海南萬寧至洋浦高速公路橫貫海南島，全長163 km，沿線開挖并用于路基填筑的高液限土約240萬m3。在高液限土填方量較大的第六、第九和第十合同段，分別選取K67+190～K67+290 (WY-6)、K110+180～K110+280 (WY-9)和K121+460～K121+560 (WY-10)作為試驗段。對試驗段填料進行基本性質(zhì)指標測試，結(jié)果如表1所示。由于高液限土天然含水率均遠大于其最佳含水率，在海南多雨氣候條件下，很難通過翻曬將高液限土填料含水率大幅降至最佳含水率壓實。因此，在滿足下路堤填料強度和變形要求前提下，在3個試驗段均只將高液限土直接用作下路堤填料，在大于最佳含水率8%～10%狀態(tài)下壓實，并適當(dāng)降低壓實度控制標準[24]。

表1 試驗段高液限土土樣基本性質(zhì)Table 1 Basic properties of soil samples of trial sections

利用高液限土填筑下路堤后，分別在3個試驗段下路堤頂面沿縱向每10 m一個測點，依次進行PFWD檢測和承載板試驗。使用的PFWD設(shè)備主要參數(shù)為：載荷板半徑0.15 m，載荷范圍1～15 kN，載荷脈沖形狀為半正弦波；載荷傳感器相對誤差小于1%，絕對誤差±0.1 kN；位移傳感器有效測試精度1 μm，最大測試范圍2 200 μm；數(shù)據(jù)采集周期0～60 ms。PFWD測試方法和步驟參照文獻[2]進行?，F(xiàn)場承載板試驗參照《公路路基路面現(xiàn)場測試規(guī)程》(JTG E 60—2008)相關(guān)規(guī)定進行。

限于篇幅，僅列出3個試驗段部分代表性PFWD試驗結(jié)果，如圖1所示。從圖1(a)～(c)可以看出，3個測點PFWD荷載與豎向位移時程曲線均有相似性，當(dāng)施加的荷載達到峰值時，該時刻對應(yīng)的豎向位移并未達到峰值，變形未達到最大；當(dāng)位移達到峰值時，荷載已處于減小階段。從圖1(d)可以看出，3個測點的荷載-位移曲線表現(xiàn)出明顯的非線性和滯回圈。因此，沖擊荷載作用下高液限土填筑層變形滯后效應(yīng)顯著。

圖1 高液限土下路堤PFWD測試結(jié)果Fig.1 PFWD results of lower embankment of

PFWD常規(guī)模量反算法是假設(shè)荷載與變形關(guān)系近似為線彈性，與承載板法一樣，均采用圓形垂直剛性分布靜荷載作用下的彈性半空間布辛尼斯克解(見式(1))計算路基模量。

(1)

式中：E為路基模量值，MPa；p為載荷板壓力峰值，MPa；R為載荷板半徑，m；μ為泊松比，取0.35；w為載荷板中心彎沉峰值，mm；π為圓周率。

由于沖擊動荷載作用下變形存在滯后，且線彈性靜力模型無法反映高液限土路基的非線性動力響應(yīng)特征，若仍采用式(1)反算模量會造成結(jié)果的偏差。現(xiàn)場試驗結(jié)果表明，圖1(a)～(c)3個測點的常規(guī)反算模量分別為27.6、24.4、31.5 MPa，而實測模量分別為18.5、15.1、21.2 MPa。3個試驗段30個測點的常規(guī)反算模量與實測模量如圖2所示，統(tǒng)計的常規(guī)反算模量平均絕對誤差10.9 MPa，平均相對誤差52.5%。

圖2 高液限土下路堤常規(guī)反算與實測模量Fig.2 Back-calculated and measured modulus of lower embankment of high liquid limit

2 基于粘彈性模型的模量反算方法

由于以往PFWD反算路基模量時采用的是線彈性靜力學(xué)模型，而沖擊荷載作用下高液限土路基具有顯著變形滯后效應(yīng)和荷載-位移非線性特征，常規(guī)反算方法已不適用，因此，有必要建立一種可較好地反映沖擊荷載作用下高液限土路基非線性動力響應(yīng)特征的簡便模量反算方法。

2.1 選用的粘彈性理論模型

眾所周知，在正弦周期動荷載作用下，粘彈性材料的應(yīng)變也呈周期變化，而且應(yīng)變響應(yīng)滯后于應(yīng)力響應(yīng)。粘彈性理論模型也普遍應(yīng)用于行車荷載作用下瀝青路面動力響應(yīng)分析?？紤]到?jīng)_擊荷載作用下高液限土路基變形滯后效應(yīng)與上述粘彈性材料力學(xué)行為相似，因此，基于粘彈性本構(gòu)模型建立了高液限土模量反算新方法。

粘彈性模型中，Maxwell模型是由彈簧-阻尼器串聯(lián)而成，主要用于模擬粘彈性材料應(yīng)力松弛過程。三參量固體模型、Burgers模型和廣義Maxwell模型等是用彈簧-阻尼器并聯(lián)和串聯(lián)多種形式構(gòu)造出的更加復(fù)雜的力學(xué)模型，可以描述動荷載作用下材料的粘彈性以及殘余變形性質(zhì)。而PFWD測試時，會通過多次落錘消除測點的塑性變形，因此，反算高液限土路基模量不需考慮殘余變形，不必采用形式復(fù)雜的多參數(shù)粘彈性模型增加反算難度。

Kelvin模型是由彈簧-阻尼器并聯(lián)構(gòu)成的粘彈性模型[11]。該模型形式簡潔，便于高效地迭代運算。其數(shù)學(xué)描述見式(2)。

(2)

(3)

根據(jù)式(3)可得到加卸載作用下的應(yīng)力-應(yīng)變曲線，如圖3所示。由圖3可以看出，當(dāng)引入粘滯系數(shù)后，加載時應(yīng)力首先增至最大值，應(yīng)變隨后增至最大值，應(yīng)變具有滯后性，應(yīng)力-應(yīng)變曲線形式與圖1(d)所示的PFWD實測荷載-位移曲線形式極為相似。因此，選用Kelvin模型用于沖擊荷載作用下高液限土路基動力響應(yīng)分析。

圖3 Kelvin模型應(yīng)力-應(yīng)變曲線[25]Fig.3 Stress and strain curve of Kelvin

2.2 準靜態(tài)反算方法

動態(tài)法和準靜態(tài)法[15,26]是動荷載作用下彈性體系動力響應(yīng)分析的兩種方法。

動態(tài)法反算[27]是將Navier動力方程中引入粘彈性理論模型，通過Laplace和Hankel反算導(dǎo)出豎向位移的顯式解，并考慮復(fù)雜混合邊值問題反算模量。由于反算時隱式偏微分方程結(jié)構(gòu)較復(fù)雜，需采用諸多簡化數(shù)值算法進行，比較耗時。

采用準靜態(tài)法反算思路是將靜荷載作用下彈性半空間體表面的豎向位移顯式解作Laplace變換，并將變換式中的彈性參數(shù)代之以粘彈性算子再作反算。這種方法可用動載數(shù)據(jù)反算靜力參數(shù)。后續(xù)的結(jié)果也驗證了對于高液限土路基采用準靜態(tài)法進行反算，其結(jié)果精度能夠滿足要求。因此，采用該方法和PFWD動載數(shù)據(jù)反算的模量屬于靜模量，從而與承載板法測得的靜模量在類型上保持一致。

軸對稱荷載作用下的彈性半空間體表面的豎向位移計算式為

(4)

(5)

由式(5)可知，當(dāng)t=0時，豎向位移為0；當(dāng)t趨近于無窮時，豎向位移等于傳統(tǒng)彈性半空間理論計算得到的豎向位移。粘滯系數(shù)η所在的項，隨著時間的變化是一個始終大于0小于1的數(shù)，限制了整體豎向變形，使理論計算位移時程的峰值縮小和相位向右延展。

2.3 模量反算程序的編制

由式(5)的形式可知，反算參數(shù)為模量值E和粘滯系數(shù)η。因式(5)比較復(fù)雜，難以得到反算參數(shù)的解析解，故采用最小二乘法，即式(6)迭代求極值的優(yōu)化算法進行反算。

(6)

式中：min(Δ)為實測值與理論值的誤差平方和；K為精度控制參數(shù)；W為豎向位移實測值，mm；w為豎向位移理論計算值，mm；β=(E,η)為反算參數(shù)向量；i和n分別為PFWD實測位移時程數(shù)據(jù)編號和數(shù)據(jù)總個數(shù)。

迭代的過程采用修正Gauss-Newton法，即令min(Δ)=f(β)，將式(6)變換為

f(β)=KrT(β)r(β)

(7)

式中：K的意義同式(6)；r(β)為實測與理論值的殘差列向量；T為向量轉(zhuǎn)置符。

在計算時，首先設(shè)定初始計算參數(shù)向量β0，一般按經(jīng)驗取一個初始模量和粘滯系數(shù)值，并按式(8)確定下一個迭代參數(shù)向量。

βk+1=βk+λkAk

(8)

迭代循環(huán)過程采用式(9)進行控制。

(9)

式中：T為向量轉(zhuǎn)置符；ε為迭代控制參數(shù)。

按照式(10)計算得到最終迭代值βg，即模量和粘滯系數(shù)最優(yōu)反算值。

βg=(βk+1+βk)/2

(10)

由于基于最小二乘優(yōu)化算法得到的是采用函數(shù)定義域內(nèi)全局求極值的反算結(jié)果，該反算結(jié)果可能與實際情況相差較大，因此，在編制反算程序時，針對PFWD檢測高液限土路基模量的具體情況作了如下處理：由于數(shù)據(jù)采集周期0～60 ms內(nèi)落錘會在橡膠塊上振動，進而帶動承載板產(chǎn)生小幅跳動，故只取0～25 ms內(nèi)峰值波形的數(shù)據(jù)；計算點位r=0 m，即荷載作用面表面的中心處；精度控制參數(shù)K=0.5和迭代控制系數(shù)ε=0.01，兩個控制值可按照不同精度要求進行調(diào)整。

進行反算時還需要將式(8)進行變形得到式(11)，以控制反算模量值的權(quán)重。當(dāng)計算位移時程的相位大于實測位移時程或計算位移峰值小于實測值，但迭代還未滿足式(9)的要求時，停止參數(shù)η的迭代，而繼續(xù)進行參數(shù)E的迭代。

βk+1=βk+λkAkB

(11)

式中：B為迭代權(quán)重控制向量，按照式(12)計算。

(12)

式中：E(βk)為實測位移和計算位移的峰值誤差限制條件，按E(βk)=|Wmax-wmax(βk)|計算；Wmax為實測位移峰值，mm；wmax(βk)為計算位移峰值，mm。γ為峰值誤差控制量，mm，取γ=0.002 mm，可根據(jù)具體要求適當(dāng)調(diào)整，以改變最終反算精度。

得到反算的最優(yōu)解后，正算位移，并用式(13)評價計算精度。

(13)

式中：R2為精度系數(shù)；W、w、i和n的意義同式(6)；βk為參數(shù)向量。

按照上述方法，利用Matlab編寫了考慮沖擊荷載作用下變形滯回效應(yīng)的高液限土路基模量反算程序[28]，既適用于FWD反算，也適用于PFWD反算。反算流程如圖4所示。

圖4 高液限土路基模量反算流程Fig.4 Back-calculation process for modulus of high

3 模量反算新方法的驗證與應(yīng)用

根據(jù)圖1所示3個測點的PFWD荷載和位移時程曲線，利用本文提出的基于粘彈性模型的高液限土路基模量反算方法(簡稱新方法)進行了模量反算，并與常規(guī)方法的反算結(jié)果和實測結(jié)果進行了對比，如表2所示。3個測點的反算結(jié)果表明，新方法反算的模量僅比相應(yīng)的實測模量大2 MPa左右，而常規(guī)反算模量比實測模量要大10 MPa左右，新方法的模量反算精度提高了40%左右。

表2 高液限土下路堤部分測點模量反算結(jié)果對比Table 2 Comparison of several back-calculation modulus of lower embankment of high liquid limit soil

分別利用新方法和常規(guī)方法反算得到的模量正算位移時程，并與實測位移時程進行對比。其中，圖5為測點K110+200位移時程曲線對比結(jié)果。由圖5可以看出，基于粘彈性模型的新方法計算得到的位移時程與實測位移時程非常吻合。根據(jù)式(13)，新方法的計算精度系數(shù)則為0.935，常規(guī)方法得到的精度系數(shù)僅為0.693。K110+200和K121+500精度系數(shù)對比結(jié)果相似，一同列于表2。以上位移時程與精度系數(shù)對比結(jié)果表明，引入粘彈性模型后，新方法能較為準確地反映沖擊荷載作用下高液限土路基變形響應(yīng)滯后這一特征。

圖5 根據(jù)常規(guī)方法與新方法計算的K110+200位移時程曲線對比Fig.5 Comparison of displacement-time curves of K110+200 calculated according to conventional and new

按照新方法對3個試驗段30個測點進行模量反算的結(jié)果如圖6所示。對比圖2可看出，新方法反算模量相對常規(guī)反算模量更接近實測模量，平均絕對誤差僅為2.6 MPa，平均相對誤差為9.2%。同時也證明了，由于采用了準靜態(tài)分析方法，可根據(jù)沖擊荷載作用下路基動力響應(yīng)(荷載和位移時程數(shù)據(jù))反算靜力參數(shù)，從而使得新方法反算模量更接近由承載板法實測的靜模量。

圖6 高液限土下路堤頂面的新方法反算模量與實測模量Fig.6 New back-calculation and measured modulus on the top of lower embankment of high liquid limit

良好的路基施工質(zhì)量必然要求在各層位施工完成后路基強度變異性小，均勻一致[29]。為此，應(yīng)用PFWD和模量反算新方法，對海南萬洋高速公路高液限土下路堤頂面模量的變異性進行了評價，快速、準確地確定了模量過低需要補壓處理的點位，減小了下路堤模量的變異性。由于高液限土下路堤模量較低，而萬洋高速公路路基頂面設(shè)計彎沉值為120(0.01 mm)，需要在上路堤和路床填筑合適填料進行剛度補償。為此，采用文獻[2]的方法，繼續(xù)以WY-9為試驗段，進行了路基剛度補償設(shè)計和現(xiàn)場試驗，確定并實施了用高液限土摻30%砂礫填上路堤，用高液限土摻40%砂礫填下路床，用低液限土摻4%水泥填上路床方案。實測路基彎沉的代表值為88(0.01 mm)，滿足驗收要求。

為研究新方法對于高模量路基測試的適用性，在WY-9試驗段路基頂面分別進行了PFWD與承載板試驗。其中，K110+200處PFWD測試結(jié)果如圖7所示。與圖1(b)、(d)相比，由于路基模量提升，沖擊荷載作用下變形滯后效應(yīng)和荷載-位移滯回圈雖顯著減小，但仍然存在，同樣會導(dǎo)致常規(guī)反算誤差。承載板實測與PFWD模量反算結(jié)果表明：實測模量為196.1 MPa，常規(guī)反算模量為221.10 MPa，新方法反算模量為205.6 MPa。

圖7 K110+200處路基頂面PFWD測試結(jié)果Fig.7 PFWD results on the top of subgrade of

圖8為WY-9試驗段路基頂面所有10個測點的模量實測與反算結(jié)果。統(tǒng)計結(jié)果表明：新方法反算模量平均絕對誤差為8.5 MPa，平均相對誤差僅為4.3%；而常規(guī)反算模量平均絕對誤差為24.1 MPa，平均相對誤差為12.2%?？梢钥闯?，當(dāng)路基模量較高時，新方法和常規(guī)方法反算模量的平均相對誤差均有所減小，但新方法反算結(jié)果的精度更高。

根據(jù)高液限土路基試驗段現(xiàn)場試驗結(jié)果，由于在低模量的下路堤和高模量的路床頂面新方法反算模量與實測模量之間的相對誤差均在10%以內(nèi)，因此，可將新方法反算的模量直接作為高液限土路基剛度補償設(shè)計的參數(shù)。

圖8 WY-9試驗段路基頂面的反算模量與實測模量Fig.8 Back-calculation and measured modulus on the top of subgrade of WY-9 test

4 結(jié)論

1)沖擊荷載作用下高液限土路基變形峰值滯后于荷載峰值，荷載-位移曲線彈性滯回明顯。

2)基于線彈性靜力模型的PFWD模量常規(guī)反算方法忽略了路基非線性動力響應(yīng)特征，用于高液限土路基模量測定時，反算模量與實測模量偏差較大。

3)基于粘彈性模型和準靜態(tài)法的PFWD模量反算新方法既可反映路基非線性動力響應(yīng)特征，又可用動載數(shù)據(jù)反算靜模量；在下路堤和路床頂面，采用新方法反算的模量均與實測模量接近，誤差在10%以內(nèi)。

4)模量反算新方法可用高液限土路基各層位頂面模量變異性控制，建議將新方法反算模量作為路基剛度補償?shù)脑O(shè)計參數(shù)。