雙層地基中傾斜偏心荷載下基樁非線性模型試驗(yàn)研究

李微哲+婁平

摘 要：引入初始地基比例系數(shù)、初始地基系數(shù)、水平位移特征值，導(dǎo)出一種新型雙曲線型p-y曲線，并給出了參數(shù)經(jīng)驗(yàn)選取參考值范圍。在上層粘土下層砂土地基中進(jìn)行了2組模型樁試驗(yàn)，含2根水平加載的鋼樁和6根傾斜偏心荷載下的木樁。試驗(yàn)表明，水平荷載、豎向荷載和偏心彎矩荷載下柔性木質(zhì)模型樁非線性特征顯著。按雙曲線型p-y曲線計(jì)算的模型樁地面處和樁頂水平位移與實(shí)測值吻合良好，驗(yàn)證了雙曲線型p-y曲線。水平位移特征值增大時(shí)，樁身位移將減小，樁身最大剪力將增大，樁身最大彎矩將減小，樁側(cè)最大土壓力將增大。雙層地基中，下層砂土位移特征值變化對試樁受力影響極小，從簡化計(jì)算參數(shù)和應(yīng)用角度出發(fā)，實(shí)際應(yīng)用時(shí)可將下層砂土水平位移特征取值與上層粘土的相同。

關(guān)鍵詞：雙層地基；木質(zhì)模型樁；初始地基比例系數(shù)；初始地基系數(shù)；水平位移特征值

中圖分類號(hào)：TU470 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1674-4764（2018）01-0068-10

Nonlinear behavior of pilesunder inclined and eccentric loads indouble layered soils by model test

Li Weizhe1，2， Lou Ping2

（1.School of Civil Engineering， Central South University， Changsha 410000，P.R.China；2. China Coal Technology & Engineering Group， Chongqing Design & Research Institute Co. Ltd.， Chongqing 400016，P.R.China）

Abstract：A new hyperbolic p-y curve was derived to generally simulate then onlinearity of piles in layered soils， and the new p-y curve was presented containing parameters such as the initial ratio of elastic foundation coefficient， the initial elastic foundation coefficient，characteristic value of lateral displacementfor soils， and etc.. Model test of six timber piles and two steel piles in two groups settled in layer soils with upper clay and lower sand were done. The calculated lateral displacements of model piles both at the top and at the ground are consistent with the test values by p-y curve method， and the new hyperbolic p-y curve is verified. While the characteristic value of lateral displacement for layered soils increases， the displacement and the maximum bending moment of the pile decrease， the maximum shear and the maximum soil pressure around the pile increase. The variety of characteristic value of lateral displacement for the lower sand has extremely small effect on mechanical property of the pile， then the characteristic value of lateral displacement for lower soils can be defined as the same as that of the top soil during engineering application.

Keywords：double layered soils； tiber model piles；the initial ratio of elastic foundation coefficient；the initial elastic foundation coefficient；characteristic value of lateral displacement

模擬樁側(cè)土非線性一般采用p-y曲線法。學(xué)者們對此進(jìn)行了大量研究，并基于單一土層中的模型樁水平承載力試驗(yàn)提出了很多不同形式的p-y曲線。但基于多層地基中的水平荷載、豎向荷載、偏心彎矩荷載下的基樁非線性室內(nèi)模型試驗(yàn)研究則較少。工程基樁大多位于成層地基中，不僅承受水平荷載，還有豎向荷載和偏心彎矩荷載，因此，基于多層地基中復(fù)雜荷載下的模型試驗(yàn)研究基樁非線性特性具有實(shí)際意義。

Matlock[1]、Reese等[2]、Stevens等[3]提出粘性土p-y曲線，其樁側(cè)土壓力與樁身水平位移成冪函數(shù)關(guān)系；王惠初等[4]、田平等[5]基于鋼管樁模型試驗(yàn)提出的河海大學(xué)新統(tǒng)一法中樁側(cè)粘土壓力與樁身水平位移為雙曲線函數(shù)關(guān)系。王騰等[6]基于鋼管模型樁試驗(yàn)提出的粉土p-y曲線、戚春香等[7]提出的弱飽和土p-y曲線，其樁側(cè)土壓力與水平位移為雙曲正切函數(shù)關(guān)系。王國粹等[8]認(rèn)為砂土p-y曲線亦可用雙曲線函數(shù)模擬。李雨潤等[9]提出了液化土中p-y曲線的修正算法，凌賢長等[10]、劉紅軍等[11-12]較為全面地介紹了液化土p-y曲線研究進(jìn)展。樓曉明等[13]提出了一種根據(jù)p-y曲線計(jì)算地基比例系數(shù)的方法，首次將p-y曲線與m法結(jié)合了起來。劉紅軍等根據(jù)玻璃管模型樁試驗(yàn)提出了黃河水下三角洲快速沉積粉土層的p-y曲線，樁側(cè)土壓力與樁身水平位移為雙曲線函數(shù)關(guān)系。上述模型樁試驗(yàn)提出p-y曲線側(cè)重單一類型土非線性研究，且模型樁試驗(yàn)荷載僅限于水平荷載。endprint

用上層粘土、下層砂土模擬成層地基，進(jìn)行了兩組木質(zhì)模型樁試驗(yàn)，一組模型樁按水平荷載→豎向荷載→偏心彎矩順序加載，另一組按豎向荷載→水平荷載→偏心彎矩順序加載，以此模擬樁頂復(fù)雜荷載工況，進(jìn)而研究復(fù)雜荷載下基樁非線性特性。根據(jù)王惠初等[14]提出的粘土p-y曲線和王國粹等[8]提出的砂土p-y曲線，推導(dǎo)了一種新的雙曲線型p-y曲線，并用于模型樁試驗(yàn)結(jié)果分析。

1 雙曲線型p-y曲線

1.1 粘土p-y曲線

根據(jù)王惠初等[4]提出的粘土p-y曲線為

式中：y為水平位移；p為發(fā)生水平位移y時(shí)的樁側(cè)土壓力；pu為極限土壓力；D為樁徑（在應(yīng)用于計(jì)算系數(shù)A時(shí)，單位為m）；ε50為土三軸試驗(yàn)中最大主應(yīng)力達(dá)到極限主應(yīng)力一半時(shí)的應(yīng)變值；β=9～11，軟土取小值，硬土取大值。A為經(jīng)驗(yàn)系數(shù)，軟土和大變位樁取大值，硬土和小變位樁取小值。Matlock建議取0.5～2.5。

根據(jù)API規(guī)范，粘土極限抗力計(jì)算為

式中：J為與土類別有關(guān)的經(jīng)驗(yàn)系數(shù)；其余符號(hào)意義同式（1）。

整理式（1）和式（2）可得淺層粘土p-y曲線

式中：yL為粘土水平位移特征值，m0為粘土初始地基比例系數(shù)，C0為初始地基系數(shù)。當(dāng)經(jīng)驗(yàn)選定了ε50值后，即可按式（1）和式（3）計(jì)算yL、m0、C0。

工程實(shí)際土層較多，建議m0、C0值經(jīng)驗(yàn)選取或根據(jù)水平承載力試驗(yàn)確定。則式（3）可寫成

式中：c為粘土粘聚力；φ為內(nèi)摩擦角；J為大于1的系數(shù)，一般可取1.5～2.0（模型試驗(yàn)中取值為1.8）；初始地基比例系數(shù)m0可按規(guī)范經(jīng)驗(yàn)取值的2～4倍選取，或根據(jù)基樁水平承載力試驗(yàn)，取其1～3級(jí)水平荷載位移值反算初始地基比例系數(shù)m0值再乘以大于1的經(jīng)驗(yàn)系數(shù)確定，一般取1.1～2.0，非線性特征小取小值，非線性特征大則取大值；確定了初始地基比例系數(shù)后則可按式（4）計(jì)算確定初始地基系數(shù)。

1.2 砂土p-y曲線

王國粹等[8]提出的砂土p-y曲線為

式中：Ki為計(jì)算點(diǎn)初始地基系數(shù)，kN·m-3；pu為計(jì)算點(diǎn)樁側(cè)土極限土壓力，kPa；y為計(jì)算點(diǎn)基樁水平位移，m。

根據(jù)《海上固定平臺(tái)的入級(jí)與建造規(guī)范》（1992），水面以上樁側(cè)砂土極限土壓力（kPa）計(jì)算式為

式中：φ為砂土內(nèi)摩擦角；γ為砂土容重。

砂土內(nèi)摩擦角φ=20～45°時(shí)，zr=8.6～29.7D，即使假定深層砂土的極限土壓力隨埋深線性增長，其對水平荷載下基樁受力影響也極其有限。地面處樁側(cè)砂土極限土壓力不為0，與砂土弱粘聚力實(shí)際情況吻合。

有學(xué)者認(rèn)為砂土初始地基系數(shù)隨埋深線性增加，也有認(rèn)為與埋深呈非線性關(guān)系，本文亦認(rèn)為砂土初始地基系數(shù)與埋深關(guān)系如下呈線性關(guān)系，且與樁側(cè)極限土壓力呈線性關(guān)系。

式中：m0為砂土初始地基比例系數(shù)，可參照現(xiàn)行樁基設(shè)計(jì)規(guī)范選?。籆0為初始地基系數(shù)；yL為砂土水平位移特征值。

由式（5）、式（7）和式（8）可得砂土p-y曲線的另一種形式

當(dāng)經(jīng)驗(yàn)選定了初始地基比例系數(shù)，則可按式（8）計(jì)算初始地基系數(shù)C0和砂土位移特征值yL。

如忽略砂土弱粘聚力，則式（8）簡化為

式中初始地基比例系數(shù)、水平位移特征值可經(jīng)驗(yàn)選取或根據(jù)水平承載力試驗(yàn)確定。

1.3 水平位移特征值的物理意義

式（4）和式（9）中的水平位移特征值yL反映了樁側(cè)土非線性特征，通過雙曲線型非線性修正因子yL（yL+y）對初始地基系數(shù)和初始地基比例系數(shù)進(jìn)行修正的方式計(jì)入樁側(cè)土非線性。當(dāng)樁側(cè)土水平位移增大時(shí)，非線性修正因子yL（yL+y）減小。當(dāng)水平荷載很小時(shí)，水平位移也很小，地基系數(shù)非線性修正因子趨于1，即yL/（yL+y→0）=1，因此，可根據(jù)基樁水平承載力試驗(yàn)初始1～2級(jí)水平位移反算得的地基系數(shù)和地基比例系數(shù)，再乘以大于1的修正系數(shù)即得初始地基系數(shù)或初始地基比例系數(shù)。

2 模型樁試驗(yàn)

進(jìn)行了兩組上層粘土下層砂土中的模型樁試驗(yàn)（圖1），每組模型試驗(yàn)包含2根水平加載的鋼模型樁和3根傾斜荷載下的木制模型樁，模型樁直徑、樁長、彈性模量、抗彎剛度等參數(shù)均由試驗(yàn)測定，見表1。鋼模型樁主要用于確定初始地基系數(shù)。

試驗(yàn)土分層壓實(shí)。上覆粘土層厚30 cm。02組試驗(yàn)粘土含水量w=24.18%，容重γ=19.05 kN·m-3，液性指數(shù)IL=0.115，不排水快剪試驗(yàn)測得c=44.91，φ=9.16°。03組試驗(yàn)粘土含水量w=23.71%，容重γ=18.81 kN·m-3，IL=0.081，不排水快剪試驗(yàn)測得c=38.35，φ=10.13°。

下臥中砂層厚70 cm，顆粒容重26.49 kN/m3。02組試驗(yàn)砂土含水量w=4.34%，γ=14.81 kN·m-3，φ=32.79°。03組試驗(yàn)砂土含水量w=4.03%，γ=14.43 kN·m-3，φ=33.42°。

加載裝置由自平衡豎向加載系統(tǒng)、重物懸臂式偏心加載系統(tǒng)、重物滑輪式水平加載系統(tǒng)組成。荷載穩(wěn)定標(biāo)準(zhǔn)為0.01 mm/5 min。02組模型樁按水平荷載→豎向荷載→偏心彎矩順序加載，本組03#～05#木質(zhì)模型樁樁頂分級(jí)加載，其加載順序和樁頂各級(jí)水平荷載、豎向荷載、偏心彎矩?cái)?shù)值見表2。03組則按豎向荷載→水平荷載→偏心彎矩順序加載。本組08#～10#木質(zhì)模型樁樁頂分級(jí)加載，其加載順序和樁頂各級(jí)水平荷載、豎向荷載、偏心彎矩?cái)?shù)值見表2。采用不同加載順序，可獲得基樁以下多種不同的受力工況：1）水平受荷樁，如03#～05#模型樁在未施加豎向力之前的狀態(tài)；2）考慮施工誤差下的豎向受荷樁，如08#～10#模型樁僅施加了豎向力而未施加水平力的狀態(tài)；3）安裝了固定支座或其他原因無偏心彎矩荷載的傾斜荷載樁，如03#～05#、08#～10#模型樁施加了水平力和豎向力，但尚未施加偏心彎矩荷載時(shí)的狀態(tài)；4）安裝了活動(dòng)支座的傾斜偏心荷載樁，如03#～05#、08#～10#模型樁施加了水平力和豎向力和偏心彎矩荷載時(shí)的狀態(tài)。endprint

3 試驗(yàn)結(jié)果分析

對水平荷載、豎向荷載、偏心彎矩荷載下模型樁非線性分析采用文獻(xiàn)[14-15]中的有限桿單元法，并自編了Matlab有限元分析程序。

02組、03組模型試驗(yàn)的初始地基系數(shù)、初始地基比例系數(shù)如表3所示。模型樁計(jì)算分析采用的水平位移特征值yL如表4所示。

水平荷載、豎向荷載、彎矩荷載作用下模型樁地面處水平位移和樁頂水平位移如表2和圖2～圖13所示。

由圖2～圖7可知，按水平荷載→豎向荷載→偏心彎矩的順序?qū)δ举|(zhì)模型樁頂加載試驗(yàn)，模型樁非線性特征顯著。比較圖2、圖4、圖6可知，地基土非線性越顯著，其水平位移特征值越小，反之亦然；比較圖3、圖5、圖7結(jié)論亦同。

由圖8～圖13可知，按豎向荷載→水平荷載→偏心彎矩的順序?qū)δ举|(zhì)模型樁頂加載試驗(yàn)，模型樁非線性特征顯著。比較圖8、圖10、圖12可知，地基土非線性越顯著，其水平位移特征值越小，反之亦然；比較圖9、圖11、圖13結(jié)論亦同。

6根模型樁共有16個(gè)水平荷載下試樁實(shí)測位移、17個(gè)豎向荷載下計(jì)入試樁傾斜實(shí)測位移、26個(gè)傾斜荷載下試樁實(shí)測位移、23個(gè)傾斜偏心荷載下試樁實(shí)測位移結(jié)果與按本文理論分析結(jié)果基本一致，從多工況不同受力狀態(tài)驗(yàn)證了本文p-y曲線。

02組6根模型樁對應(yīng)的粘土水平位移特征值均不同，最大值為0.85 mm，最小值為0.48 mm，平均值0.617 mm，樣本標(biāo)準(zhǔn)差0.13 mm，變異系數(shù)為0.21<0.25，因此，試驗(yàn)中粘土水平位移特征值取平均值是合理的估計(jì)，即為0.617 mm。

5#模型樁計(jì)算結(jié)果圖6和圖7所示，計(jì)算位移值1為按本文p-y曲線法計(jì)算，其計(jì)算參數(shù)如表3和表4。計(jì)算位移值2為按文獻(xiàn)[17]中宏觀下基于m法下的p-y曲線法計(jì)算結(jié)果，其計(jì)算關(guān)鍵參數(shù)初始地基比例系數(shù)為27 N·cm-4，水平位移特征值為0.4 mm。其位移計(jì)算結(jié)果基本一致，但因宏觀下基于m法的p-y曲線法忽略了粘土地面處初始地基系數(shù)，從而使得初始地基比例系值大大增加，其水平位移特征值則減小。

4 水平位移特征值對基樁受力影響

粘土非線性特征一般比下砂土更為顯著，其下層中砂水平位移特征值yLs比上層粘土的水平位移特征值yLc更大。以5#模型樁為例，僅改變下層砂土位移特征值的情況下，分別假定yLs=yLc和yLs=100yLc的試樁內(nèi)力位移計(jì)算結(jié)果如表5?？梢?，下層砂土水平位移特征值放大100倍后，樁身最大彎矩、剪力以及樁側(cè)最大土壓力變化不超過0.5%，樁頂最大位變化不超過1.67%。可見，下層砂土位移特征值變化對試樁受力影響極小，因此，從簡化計(jì)算參和應(yīng)用角度出發(fā)，假定下層砂土的水平位移特征值與上層粘土的相同是合理可行的。

以5#模型樁為例，保持初始地基系數(shù)、初始地基比例系數(shù)不變，僅改變其雙層土水平位移移特征值時(shí)基樁計(jì)算結(jié)果如表5和圖14～圖16。可見，水平位移特征值增大后，樁身位移將減小，樁身最大剪力將增大，樁身最大彎矩將減小，樁側(cè)最大土壓力將增大。

5 結(jié)論

1）提出一種新的雙曲線型p-y曲線，即式（4）和式（9），并給出了其關(guān)鍵參數(shù)初始地基系數(shù)、初始地基比例系數(shù)、水平位移特征值的選取確定方法。

2）水平位移特征值在本次模型試驗(yàn)中表現(xiàn)出一定的離散型，最大值為0.85 mm，最小值為0.48 mm，平均值0.617 mm，樣本標(biāo)準(zhǔn)差0.13 mm，變異系數(shù)為0.21<0.25，因此，試驗(yàn)中水平位移特征值取平均值是合理的估計(jì)。

3）分別按按水平荷載→豎向荷載→偏心彎矩的順序和按豎向荷載→水平荷載→偏心彎矩的順序?qū)δ举|(zhì)模型樁頂加載試驗(yàn)， 6根模型樁共有16個(gè)水平荷載下試樁實(shí)測位移、17個(gè)豎向荷載下計(jì)入試樁傾斜實(shí)測位移、26個(gè)傾斜荷載下試樁實(shí)測位移、23個(gè)傾斜偏心荷載下試樁實(shí)測位移結(jié)果與理論分析結(jié)果基本一致，從多工況不同受力狀態(tài)驗(yàn)證了本文p-y曲線。

4）下層砂土位移特征值變化對試樁受力影響極小，因此，從簡化計(jì)算參和應(yīng)用角度出發(fā)，假定下層砂土的水平位移特征值與上層粘土的相同是合理可行的。從而，為多層地基選取統(tǒng)一的水平位移特征值提供理論支持。

5）水平位移特征值增大后，樁身位移將減小，樁身最大剪力將增大，樁身最大彎矩將減小，樁側(cè)最大土壓力將增大。

6）文獻(xiàn)[16]中宏觀下基于m法的p-y曲線應(yīng)用在上層粘土的多層地基中基樁非線性分析時(shí)存在不足，因其忽略了地面處粘土初始地基系數(shù)，使得初始地基比例系數(shù)大大增加，可以引入初始地基系數(shù)加以修正。參考文獻(xiàn)：

[1] MATLOCK H S. Correlation for design of laterally loaded pile in soft clay [C]//Proceedings of 2nd Offshore Technology Conference， Houston：[s.n.]，1970.

[2] REESE L C， COX W R， KOOP F D. Field testing and analysis of laterally loaded piles in stiff clay [C]//Proceedings of 7th Offshore Technology Conference， Houston： [s.n.]， 1975.

[3] STEVENS J， AUDIBERT J. Re-examination of p-y curve formulations [C]//Proceedings of 11th Offshore Technology Conference， Houston： [s.n.]， 1979.

[4] 王惠初，武冬青，田平.黏土中橫向靜載樁p-y曲線的一種新的統(tǒng)一法[J].河海大學(xué)學(xué)報(bào)，1991，19（1）：9-17.endprint

WANG H C， WU D Q， TIAN P. A new unified method of p-y curve for laterally loaded pile groups in plasticity soil [J]. Journal of Hohai University，1991， 19（1）： 9-17.（in Chinese）

[5] 田平，王惠初.黏土中橫向周期性荷載樁的p-y曲線統(tǒng)一法[J].河海大學(xué)學(xué)報(bào)，1993，21（1）：9-14.

TIAN P， WANG H C. A unified method of p-y curve for laterally and circularly loaded pile in plasticity soil [J]. Journal of Hohai University， 1993， 21（1）：9-14. （in Chinese）

[6] 王騰，王天霖.粉土p-y 曲線的試驗(yàn)研究[J].巖土力學(xué)，2009，30（5）：1343-1346.

WANG T， WANG T L. Experimental research on silt p-y curves [J]. Rock and Soil Mechanics， 2009，30（5）：1343-1346. （in Chinese）

[7] 戚春香，王建華.弱化飽和土層中樁的雙曲型p-y 曲線[J]. 天津大學(xué)學(xué)報(bào)， 2010，20（3）：215-220.

QI C X， WANG J H.Hyperbolic p-y curve of pile in saturated degradation strata [J]. Journal of Tianjin University， 2010，20（3）：215-220. （in Chinese）

[8] 王國粹，楊敏.砂土中水平受荷樁非線性分析[J].巖土力學(xué)，2011，32（2）：261-267.

WANG G C， YANG M. Nonlinear analysis of laterally loaded piles in sand [J]. Rock and Soil Mechanics， 2011，32（2）：261-267. （in Chinese）

[9] 李雨潤，袁曉銘，梁艷.樁-液化土相互作用p-y曲線修正計(jì)算方法研究[J].巖土工程學(xué)報(bào)，2009，31（4）： 595-599.

LI Y R， YUAN X M， LIANG Y. Modified calculation method of p-y curves for liquefied soil-pile interaction [J]. Chinese Journal of Geotechnical Engineering， 2009，31（4）： 595-599. （in Chinese）

[10] 凌賢長，唐亮.液化場地樁基側(cè)向響應(yīng)分析中p-y曲線模型研究進(jìn)展[J].力學(xué)進(jìn)展，2010，40（3）：250-261.

LING X C， TANG L.Recent advance of p-y curve to model lateral response of pile foundation on liquefied ground [J].Advances in Mechanics， 2010， 40（3）：250-261. （in Chinese）

[11] 劉紅軍，馬明泊，呂小輝，等.黃河水下三角洲快速沉積粉土層樁基p-y曲線試驗(yàn)研究[J]. 中國海洋大學(xué)學(xué)報(bào)，2015，45（10）：81-87

LIU H J， MA M P， LYU X H， et al.Experimental study on p-y curve of rapidly deposited silt in the subaqueous Yellow River Delta [J]. Periodical of Ocean University of China，2015，45（10）：81-87. （in Chinese）

[12] 劉紅軍，西國庚，馬明泊.可液化土p-y曲線模型折減方法的研究進(jìn)展[J]. 中國海洋大學(xué)學(xué)報(bào)，2015，45（7）：107-112.

LIU H J， XI G G， MA M P.Recent advance of reduction methods on p-y curves in liquefiable soil [J].Periodical of Ocean University of China， 2015，45（7）：107-112. （in Chinese）

[13] 樓曉明，吳昊，黃江楓.基于p-y曲線確定飽和黏性土的地基比例系數(shù)[J].巖土工程學(xué)報(bào)，2012，34（12）： 2206-2212.

LOU X M， WU H， HUANG J F.Determination of slope coefficient of subgrade reaction of saturated clay based on p-y curve [J]. Chinese Journal of Geotechnical Engineering， 2012， 34（12）： 2206-2212. （in Chinese）

[14] 趙明華，李微哲，單遠(yuǎn)銘，等.成層地基中傾斜荷載樁改進(jìn)有限桿單元法研究[J].工程力學(xué)，2008，25（5）：79-85.

ZHAO M H，LI W Z， SHAN Y M， et al. Behavior analysis of piles in layered clays under eccentric and inclined loads by improved finite pole element method [J]. Engineering Mechanics，2008， 28（5）： 79-85. （in Chinese）

[15] 趙明華，李微哲，曹文貴.復(fù)雜荷載及邊界條件下基樁有限桿單元方法研究[J].巖土工程學(xué)報(bào)，2006，28（9）： 1059-1064.

ZHAO M H，LI W Z， CAO W G.Study on applying finite pole element method to analysis of piles under complex loads with different boundary restraints[J]. Chinese Journal of Geotechnical Engineering. 2006， 28（9）： 1059-1064. （in Chinese）

[16] 李微哲，婁平.m法下的雙曲線型p-y曲線[J].土木建筑與環(huán)境工程，2017，39（2）：115-122.

LI W Z， LOU P. Analysis of pile by p-y curve method based on m-method [J]. Journal of Civil， Aachitectural & Eneironmental Engineering，2012，39（2）：115-122.（in Chinese）endprint