楊 芳，龐洪臣※，王玉興

（1.廣東海洋大學(xué)機(jī)械與動(dòng)力工程學(xué)院，廣東湛江 524088；2.華南農(nóng)業(yè)大學(xué)工程學(xué)院，廣州 510642）

0 引言

農(nóng)田水利防護(hù)和洪澇災(zāi)害治理是世界各地廣泛面臨的重大問(wèn)題，防汛搶險(xiǎn)螺旋樁是根據(jù)消防要求專門設(shè)計(jì)的一種便攜式快速防汛搶險(xiǎn)裝備，該裝備可解決傳統(tǒng)擊打沉樁耗時(shí)耗力和大型沉樁設(shè)備易產(chǎn)生次生災(zāi)害等問(wèn)題。沉樁過(guò)程中樁土受力及土體大變形決定了單純運(yùn)用理論分析求解極其復(fù)雜，同時(shí)土體特性也決定了利用試驗(yàn)研究樁土作用也十分困難，因此仿真分析便成了研究沉樁過(guò)程中土體變化的重要手段。邵康等[1]對(duì)螺旋鋼樁進(jìn)行了數(shù)值模擬分析，研究了在豎向載荷作用下對(duì)螺旋樁的擾動(dòng)問(wèn)題；Danda Shi 等[2]利用DEM 對(duì)螺旋樁的沉樁速比進(jìn)行了研究；吳萍等[3]利用FLAC3D對(duì)螺旋樁的承載力進(jìn)行了數(shù)值模擬；韋謝恩等[4]建立了防汛搶險(xiǎn)螺旋樁的基本幾何參數(shù)，利用有網(wǎng)格的有限元法建立了樁體和土體模型，運(yùn)用流固耦合算法對(duì)防汛搶險(xiǎn)螺旋樁的沉樁過(guò)程進(jìn)行了仿真；龐洪臣等[5]利用無(wú)網(wǎng)格法對(duì)沉樁扭矩進(jìn)行了仿真對(duì)比分析，對(duì)樁體結(jié)構(gòu)進(jìn)行了優(yōu)化，降低了沉樁扭矩。以上對(duì)防汛搶險(xiǎn)螺旋樁的研究提供了寶貴經(jīng)驗(yàn)。為進(jìn)一步明確防汛搶險(xiǎn)螺旋樁的沉樁機(jī)理，設(shè)計(jì)省時(shí)省力的樁體結(jié)構(gòu)，本文利用SPH（Smoothed Particle Hydrodynamics）粒子可模擬破裂、大變形的優(yōu)勢(shì)構(gòu)建了土體仿真模型[6-11]，著重對(duì)沉樁過(guò)程中的土體受力情況進(jìn)行了仿真分析和研究，并利用高頻動(dòng)態(tài)壓力傳感器對(duì)樁周的土體壓力進(jìn)行了初步試驗(yàn)研究。

1 樁土建模及主要參數(shù)設(shè)置

1.1 螺旋樁建模

螺旋樁主體結(jié)構(gòu)主要包括樁體和螺旋葉，實(shí)物如圖1 所示，整個(gè)樁體為空心結(jié)構(gòu)。螺旋葉則首先采用APDL（ANSYS Parametric Design Language）建立模型外徑和內(nèi)徑的關(guān)鍵點(diǎn)，再利用樣條曲線將關(guān)鍵點(diǎn)連接起來(lái)生成螺旋葉內(nèi)側(cè)和外側(cè)螺旋線，最后沿螺旋線掃描生成螺旋葉模型。整個(gè)模型采用Belytschko-Lin-Tsay 殼單元對(duì)樁體和螺旋葉網(wǎng)格進(jìn)行劃分，有限元模型及各部分名稱如圖2所示，幾何參數(shù)如表1所示，并選用QT235作為設(shè)計(jì)材料。

圖1 螺旋樁實(shí)物

圖2 螺旋樁有限元模型

表1 螺旋樁幾何參數(shù) mm

1.2 土體建模

采用SPH 建立土體模型時(shí)，要求SPH 粒子排列應(yīng)盡可能規(guī)則和均勻，且粒子質(zhì)量差異不能太大。據(jù)此，將土體幾何模型設(shè)為圓柱，圓柱模型橫截面上采用如圖3所示的粒子排列方法，并對(duì)該圓柱模型的周邊粒子和底面粒子的所有自由度進(jìn)行約束，如圖4所示。根據(jù)多次試探性運(yùn)算，為節(jié)省運(yùn)算時(shí)間并盡可能真實(shí)模擬沉樁過(guò)程中樁土之間的作用關(guān)系，選用了高度為1 000 mm、半徑為200 mm的圓柱SPH粒子組作為土體模型，總粒子數(shù)為76 895個(gè)。土體參數(shù)設(shè)置如表2所示，樁土接觸仿真模型如圖5所示。

圖3 SPH粒子排列

圖4 SPH粒子約束

表2 土體參數(shù)

2 仿真結(jié)果分析

利用前述建模和參數(shù)設(shè)置進(jìn)行仿真計(jì)算，根據(jù)仿真結(jié)果分析如下。

2.1 樁土作用分析

圖5 樁土接觸仿真模型（SPH粒子）

圖6 仿真土體應(yīng)力分布

仿真后得出的土體變形圖如圖5 所示，土體應(yīng)力分布圖如圖6 所示，沉樁試驗(yàn)現(xiàn)場(chǎng)如圖7所示，由圖中可以直觀看出，該仿真不僅模擬出了防汛搶險(xiǎn)螺旋樁在沉樁過(guò)程中樁體對(duì)土體的擠壓效應(yīng)，螺旋葉對(duì)土體剪切使土體破裂及地表土體上移的效應(yīng)也有良好的仿真。樁周土體的應(yīng)力分布隨著沉樁深度的增加而逐漸增大，在樁體下端錐部周邊一定區(qū)域內(nèi)形成半球形塑性區(qū)。當(dāng)樁體沉入到極限位置時(shí)下端螺旋葉處的土體應(yīng)力達(dá)到屈服強(qiáng)度，這與圖1 所示的下端螺旋葉和錐體處有明顯的光亮度相符，說(shuō)明土體與螺旋樁的相互作用力主要集中在下端。

為了節(jié)省計(jì)算量，仿真所用的土體幾何模型尺寸比實(shí)際土體影響域小的多，并且如前所述對(duì)周邊土體和底面土體粒子的所有自由度進(jìn)行了約束，這從圖6顯示的土體模型邊界處仍然有較大土體應(yīng)力可以得到驗(yàn)證。

2.2 土體壓力和土體位移對(duì)比分析

在仿真模型中選取距地表300 mm處，即埋樁深度的中部位置且距樁體軸線分別為0 mm、50 mm、100 mm、150 mm的同一徑向處具有代表性的土體粒子壓力曲線和位移曲線，分別對(duì)應(yīng)如圖8、圖9 所示的A、B、C、D 曲線，進(jìn)行對(duì)比分析如下。

圖7 沉樁試驗(yàn)現(xiàn)場(chǎng)

圖8 土體粒子壓力曲線

圖9 土體粒子位移曲線

（1）從圖8和圖9中均可以看出，土體粒子壓力和位移的波動(dòng)變化主要發(fā)生在300～500 mm 深度，而該段的長(zhǎng)度恰好等于樁體前部錐體的高度，進(jìn)一步說(shuō)明沉樁過(guò)程中樁體對(duì)土的作用主要集中在螺旋樁的下端。

（2）圖8 中顯示當(dāng)錐體部分經(jīng)過(guò)該位置時(shí)，距螺旋樁軸線最近的土體粒子即A 壓力曲線波動(dòng)十分劇烈，在樁體下沉的動(dòng)態(tài)過(guò)程中表現(xiàn)出了時(shí)壓時(shí)拉的狀態(tài)，但是當(dāng)錐體部分全部通過(guò)該位置之后壓力曲線則變的比較穩(wěn)定，但其壓力值仍然比距離軸線較遠(yuǎn)處的土體粒子（B、C、D 曲線）大一些。再對(duì)比圖6中A曲線的位移情況發(fā)現(xiàn)以下情況。

①在土體壓力劇烈波動(dòng)的階段，土體位移則近似呈線性增加。也就是說(shuō)，土體粒子位移的線性增大并沒(méi)有表現(xiàn)出壓力的線性升高，因此該處土體的變化特征不是彈性體，而更多表現(xiàn)出的是土體的塑性流動(dòng)特征。

②該土體粒子位移達(dá)到最大值后又產(chǎn)生了較小的反向移動(dòng)，即回松狀態(tài)，而后位移保持穩(wěn)定。

③該處位置的土體粒子位移遠(yuǎn)遠(yuǎn)大于其他各處的土體粒子。

（3）在圖8中，B壓力曲線相比A壓力曲線波動(dòng)沒(méi)有那么明顯。但是在沉樁過(guò)程中表現(xiàn)出先呈近似線性增大，達(dá)到最大壓力值時(shí)表現(xiàn)出明顯的波動(dòng)性，而后又呈近似線性下降最后再穩(wěn)定在一個(gè)定值附近，且下降的斜率要比上升的斜率大的多。這說(shuō)明該處的土體壓力在螺旋樁前部錐體斜面的擠壓下線性增大，而當(dāng)錐體部分通過(guò)后土體粒子出現(xiàn)了回松現(xiàn)象使得壓力急劇下降。再對(duì)比圖9中B位移曲線的位移情況發(fā)現(xiàn)：

①在深度為300～400 mm 之間的壓力上升階段土體粒子的位移卻基本保持不變，具有一定的滯后性；

②在深度為400 mm處的最大壓力階段土體粒子才開(kāi)始產(chǎn)生近似線性的位移；

③在深度為450～500 mm 之間土體粒子的壓力和位移基本成同步線性下降且該處土體粒子反向移動(dòng)的距離要比A 曲線大的多；

④距離樁體較遠(yuǎn)的C 曲線和D 曲線的壓力變化較為緩慢，沒(méi)有A 曲線和B 曲線那么大的波動(dòng)，且距離軸線較近粒子對(duì)應(yīng)的C 曲線在壓力增長(zhǎng)的時(shí)間段內(nèi)斜率較D 曲線大，即壓力增長(zhǎng)快些，而后兩者壓力值皆保持恒定，但查看他們對(duì)應(yīng)的位移曲線可知兩者所對(duì)應(yīng)的位移卻基本沒(méi)有變化，既沒(méi)有表現(xiàn)出類似A曲線和B曲線的大位移，也沒(méi)有表現(xiàn)出回松的特性。

（4）圖8 中顯示在錐體經(jīng)過(guò)該位置期間距樁體軸線較遠(yuǎn)的B曲線比A曲線表現(xiàn)出更大的土體壓力。

3 土體壓力試驗(yàn)

3.1 試驗(yàn)器材

試驗(yàn)器材包括：CYY9 高頻響動(dòng)態(tài)土體壓力傳感器2 個(gè)；數(shù)字示波記錄儀（NEC RA2300 數(shù)據(jù)采集系統(tǒng)）1 套；錘頭1件；自制擴(kuò)孔鋼管等。

3.2 土體壓力傳感器主要參數(shù)

測(cè)量土體壓力采用CYY9 高頻響動(dòng)態(tài)土體壓力傳感器，如圖10 所示，頻率響應(yīng)約500 kHz，壓力滿量程為50 MPa，傳感器的厚度為10 mm、直徑為20 mm，輸出電壓信號(hào)為0～5 V，可測(cè)試土體中的靜態(tài)和動(dòng)態(tài)壓力變化。傳感器已由生產(chǎn)廠家在出廠前進(jìn)行了標(biāo)定，壓力值和輸出的電壓信號(hào)之間呈線性關(guān)系如公式（1）所示，通過(guò)采集傳感器電壓信號(hào)的變化可以獲知傳感器實(shí)際所測(cè)壓力的動(dòng)態(tài)變化。

圖10 CYY9高頻動(dòng)態(tài)土體壓力傳感器

式中：Pt為與傳感器輸出電壓對(duì)應(yīng)的動(dòng)態(tài)壓力值；VF為傳感器滿量程輸出電壓值；V0為傳感器零位電壓值；Vt為傳感器動(dòng)態(tài)輸出電壓值。

3.3 試驗(yàn)過(guò)程

at為測(cè)量防汛搶險(xiǎn)螺旋樁沉樁過(guò)程中的土體壓力變化，分別在深300 mm、距離樁體軸線150 mm 位置處埋入 2 個(gè)壓力傳感器且傳感器埋入點(diǎn)與軸線垂線成90h，如圖11所示。

根據(jù)傳感器尺寸，設(shè)計(jì)1 個(gè)圓柱形鐵管，將其前端加工成與傳感器厚度、寬度一致的扁平結(jié)構(gòu)。將此鐵管與樁體軸線平行的插入地下300 mm 深度、徑向距離樁體軸線150 mm的指定位置。拔出后將傳感器沿著鐵管擠壓土體形成的空腔放入最底部，再利用土體將空腔填滿壓實(shí)。然后，將螺旋樁沉入土體中，利用NEC RA2300數(shù)據(jù)采集系統(tǒng)記錄傳感器數(shù)輸出電壓隨時(shí)間的變化值，根據(jù)測(cè)得的數(shù)據(jù)和標(biāo)定函數(shù)公式（1），并對(duì)數(shù)據(jù)濾波后可以得到沉樁過(guò)程中2 個(gè)傳感器所測(cè)得的試驗(yàn)土體壓力曲線如圖12所示。

圖11 傳感器埋入方式

圖12 試驗(yàn)土體壓力曲線

3.4 試驗(yàn)分析

（1）由于埋入傳感器時(shí)，一方面對(duì)鐵管的錘擊將導(dǎo)致空腔周邊土體產(chǎn)生很大硬化，另一方面雖然將埋入傳感器的空腔壓實(shí)但土體已非原狀土，所以測(cè)得的土體壓力值將與真實(shí)土體狀態(tài)壓力值具有很大的誤差。

（2）圖12 中2 個(gè)傳感器所測(cè)的最大土體壓力分別約為0.03 106Pa 和0.02 106Pa，而圖8中仿真土體壓力曲線D 中顯示的最大土體壓力約為0.1 106Pa，兩者相差約一個(gè)量級(jí)，但與仿真曲線D在變化趨勢(shì)上相比兩者仍然表出現(xiàn)出了較好的一致性。

（3）根據(jù)2 條試驗(yàn)土體壓力曲線可知，雖然埋入深度相同，在與螺旋樁軸線徑向距離相同的圓周上的，不同角度位置的土體壓力仍然會(huì)不盡相同，即螺旋樁樁周同一半徑圓周上的土體壓力并非均勻分布。

（4）根據(jù)2條試驗(yàn)土體壓力曲線并結(jié)合仿真圖8曲線D均可以看出，隨著沉樁逐漸深入，在所測(cè)位置處，土體壓力初始階段近似呈線性增加，在壓力值達(dá)到頂峰之后將會(huì)出現(xiàn)一個(gè)泄壓的過(guò)程，最后壓力值將基本保持恒定。

4 結(jié)束語(yǔ)

（1）利用SPH 構(gòu)建仿真土體對(duì)沉樁過(guò)程進(jìn)行仿真計(jì)算，能較真實(shí)地仿真出沉樁過(guò)程中土體的擠壓及剪切破裂、土體的應(yīng)力分布和土體位移變化等工作情形。

（2）防汛搶險(xiǎn)螺旋樁在沉樁過(guò)程中，樁尖錐體部位會(huì)形成半球形土體塑性區(qū)，土體與樁體的相互作用力主要集中在下端錐體部位。

（3）在沉樁過(guò)程中樁周土體的壓力變化和位移變化并不同步，樁周同一圓周處的壓力不均布。

（4）緊靠樁體軸線的土體粒子表現(xiàn)出時(shí)壓時(shí)拉的受力狀態(tài)，同時(shí)該處土體粒子的大位移表現(xiàn)出土體的塑性流動(dòng)的特征。

（5）最大的土體壓力不一定發(fā)生在緊靠樁體軸線的位置，而是距樁體軸線一定的距離。

（6）試驗(yàn)和仿真均難以獲得深層土體壓力的準(zhǔn)確值，但通過(guò)兩者的相互印證仍可獲得沉樁過(guò)程中土體壓力較好的變化規(guī)律。

本文結(jié)合仿真和土體壓力試驗(yàn)分析了沉樁過(guò)程中的樁土作用狀態(tài)及土體受力變化情況，在一定程度上進(jìn)一步揭示了防汛搶險(xiǎn)螺旋樁的沉樁機(jī)理，為設(shè)計(jì)減阻降耗的螺旋樁樁體結(jié)構(gòu)提供了一定借鑒。