考慮橫向慣性效應(yīng)時飽和土中大直徑樁的縱向振動研究

龔志超， 楊冬英

(蘇州科技大學(xué) 土木工程學(xué)院，江蘇 蘇州 215011)

樁-土耦合振動問題的難點之一就是如何建立合理的土介質(zhì)模型，縱觀現(xiàn)有的研究成果，有將樁周土視為單相介質(zhì)，不考慮液相(土中水)的影響，比如Keshavarz、Wang、吳文兵等[1-3]；有將土視為固液兩相介質(zhì)的，即樁周土體為飽和土比如Zeng、Wang、李強(qiáng)、楊驍、劉林超等[4-10]，對飽和土中單樁縱向振動特性進(jìn)行了研究；也有將土視為非飽和土來研究樁-土耦合振動問題，比如張智卿、王星華等[11-12]，對非飽和土中單樁的縱向振動進(jìn)行了理論與試驗研究，取得了豐碩的成果。其中，現(xiàn)代工程中，由于工程規(guī)模日趨變大，大直徑樁運用范圍極廣，而大直徑樁的振動特性與普通直徑樁的動力特性有所區(qū)別，李強(qiáng)等基于Biot飽和土理論并將樁等效為Rayleigh-Love桿研究了飽和土中大直徑嵌巖樁的縱向振動問題，指出了大直徑樁不同于普通直徑樁的振動特性。劉林超等基于飽和多孔介質(zhì)理論，利用Novak薄層法[13]研究了飽和土中樁基的縱向振動問題，得出了樁頂導(dǎo)納和樁頂復(fù)剛度的表達(dá)式，并將研究結(jié)果與三維模型的結(jié)果對比驗證了結(jié)果的正確性。楊驍?shù)然陲柡投嗫捉橘|(zhì)理論，研究了考慮橫向慣性效應(yīng)時成層飽和土中單樁縱向振動特性，得出飽和土中樁頂復(fù)剛度的解析解。

在樁受縱向動力荷載作用時，不僅產(chǎn)生縱向變形，橫向變形也會隨之產(chǎn)生。在小應(yīng)變情況下橫向變形是否能夠忽略不計跟樁的尺寸有關(guān)，已有研究結(jié)果[14]表明對于大直徑樁來說不能忽略樁身的橫向變形對樁基縱向振動的影響。為了既能合理反映飽和土中大直徑樁振動問題，又能使此耦合問題方便解耦獲取大直徑樁-土耦合振動解，本文基于飽和多孔介質(zhì)理論和Novak薄層法，考慮樁的橫向慣性效應(yīng)，建立飽和土中樁基的縱向振動模型，并分析水的體積分?jǐn)?shù)、樁徑等參數(shù)對樁頂動力響應(yīng)的影響，以豐富飽和土中樁基的振動理論。

1 數(shù)學(xué)模型與基本假定

本文研究的是均質(zhì)飽和土中樁基的縱向振動特性，樁-土耦合作用模型如圖1所示。各向同性的均質(zhì)飽和土中有一樁長為H、樁半徑為r0的黏彈性摩擦樁，樁頂受到豎向振動力q(t) 的作用，樁周飽和土體對樁身單位長度上的摩擦力f(z)=2πr0τ。為簡化計算，作出以下假定：①樁為等截面的粘彈性樁，彈性模量和材料阻尼分別為Ep和ηp；②樁底為黏彈性支承，樁底分布式(單位面積上)彈性系數(shù)和阻尼系數(shù)分別為kb、ηb；③樁周土體為固液兩相飽和多孔介質(zhì)，并且徑向無限延伸；④樁與樁側(cè)土完全接觸，沒有相對滑移和脫離，樁側(cè)土體對樁的作用力通過樁側(cè)土復(fù)剛度傳遞。

圖1 樁-飽和土相互作用模型Fig.1 Model of saturated soil-pile interaction

2 飽和土體的動力響應(yīng)

基于飽和多孔介質(zhì)理論和Novak薄層法的假定可知飽和土體的動力控制方程[9]為

(1)

(2)

式中：

解方程式(2)得

(3)

式中：I0和K0分別為零階的第一類和零階的第二類虛宗量的Bessel函數(shù)，A和B分別為待定復(fù)系數(shù)。由于無窮遠(yuǎn)處土體位移為零，可得B=0，所以

(4)

(5)

3 樁的縱向振動方程與求解

通過Hamilton變分原理[16]可知考慮橫向慣性效應(yīng)時樁的任一截面z處所受的截面力為

(6)

式中：Ep、Ap、ηp、υp、ρp分別為樁的彈性模量、截面面積、材料阻尼、泊松比、樁身密度，u=u(z，t)。在樁的任一截面z處取一微元段隔離體進(jìn)行受力分析，可得考慮橫向慣性效應(yīng)時樁的縱向振動方程：

(7)

樁頂和樁底處的邊界條件分別為

(8)

初始條件為

(9)

結(jié)合初始條件(9)對縱向振動方程(6)進(jìn)行Laplace變換并進(jìn)行整理可得

(10)

U(z,ω)=D1cos(λz)+D2sin(λz)

(11)

式中：D1和D2為待定系數(shù)，由邊界條件決定。

由位移阻抗函數(shù)的定義(力除以位移)并結(jié)合樁底邊界條件可得樁頂位移阻抗函數(shù)，即樁頂復(fù)剛度為

(12)

Zp0=Kr+iKi

其中Kr為動剛度，Ki為動阻尼。

由式(12)可得無量綱化樁頂位移阻抗函數(shù)為

(13)

樁頂速度頻域響應(yīng)函數(shù)為

(14)

對Hv(ω)進(jìn)行無量綱化可得

(15)

(16)

4 參數(shù)影響與分析

4.1 樁徑不同時橫向慣性效應(yīng)對樁頂動力響應(yīng)的影響

當(dāng)令樁的泊松比為0時即可得到樁頂動力響應(yīng)的一維解，圖2～圖5反映的是在其他參數(shù)不變的情況下(nL=0.33，nS=0.67)，樁半徑分別為r0=1.0 m和r0=0.4 m時樁的橫向慣性效應(yīng)對大直徑樁頂動力響應(yīng)的影響。

從圖2(a)樁頂速度時域響應(yīng)曲線可以看出，由于考慮橫向慣性效應(yīng)時波的傳播過程中產(chǎn)生彌散作用，使樁底反射信號的波谷值明顯減小，反射峰值的位置明顯延遲，且反射波曲線出現(xiàn)劇烈震蕩，這將會影響對樁底支承情況的判斷；從圖2(b)樁頂導(dǎo)納曲線可以看出，在低頻段橫向慣性效應(yīng)對樁頂導(dǎo)納的影響幾乎可以忽略，但隨著無量綱角頻率的增大，共振頻率有所減小。這是由于樁的橫向變形使波在樁中的傳播是彌散的，較高的激振頻率會引起較大的能量損失。同時也說明，對于飽和土中大直徑樁基動力檢測，應(yīng)對反射信號在樁底附近振蕩的原因加以區(qū)別。通過圖2、圖4的比較可以看出，當(dāng)樁徑明顯減小時，橫向慣性效應(yīng)對樁頂速度時域響應(yīng)和樁頂導(dǎo)納的影響明顯減小，幾乎可以忽略。

圖3、圖5反映了大直徑樁的橫向慣性效應(yīng)對樁頂復(fù)剛度的影響。從圖3中可以看出在低頻段，橫向慣性效應(yīng)對樁頂復(fù)剛度的影響可以忽略，但在高頻段，是否考慮橫向慣性效應(yīng)將對大直徑樁頂復(fù)剛度產(chǎn)生明顯的影響：考慮橫向慣性效應(yīng)會使樁頂復(fù)剛度的峰值減小，并且頻率越大，峰值差越大；通過圖3與圖5的比較可以看出，當(dāng)樁徑明顯減小時，無論是低頻段還是高頻段，橫向慣性效應(yīng)對樁頂復(fù)剛度的影響均明顯減小，甚至可以忽略。綜合圖2～圖5的結(jié)論可以說明對于大直徑樁而言，有必要考慮橫向慣性效應(yīng)對樁頂動力響應(yīng)的影響，而忽略橫向慣性效應(yīng)的一維解只適用于細(xì)長樁。

圖2 樁的橫向慣性效應(yīng)對樁頂速度響應(yīng)的影響(r0=1.0 m)Fig.2 Influence of transversal inertia effect on velocity response at top of pile(r0=1.0 m)

圖3 樁的橫向慣性效應(yīng)對樁頂復(fù)剛度的影響(r0=1.0 m)Fig.3 Influence of transversal inertia effect on complex stiffness at top of pile(r0=1.0 m)

圖4 樁的橫向慣性效應(yīng)對樁頂速度響應(yīng)的影響(r0=0.4 m)Fig.4 Influence of transversal inertia effect on velocity response at top of pile(r0=0.4 m)

4.2 樁徑對樁頂動力響應(yīng)的影響

圖6、圖7反映的是在其他條件不變的情況下(nL=0.33，nS=0.67)，樁徑對樁頂動力響應(yīng)的影響。

從圖6(a)中的樁頂速度時域響應(yīng)曲線可得，隨著樁徑的增大，樁底反射信號幅值明顯增大，同時速度時域響應(yīng)曲線的振蕩也明顯加劇。從圖6(b)中可得，樁徑越大，導(dǎo)納曲線的幅值越大，共振頻率越小。而在一維解中，樁徑的增大并不會引起速度時域響應(yīng)曲線的振蕩，也不會使共振頻率減小[17]，這是由于受樁的橫向慣性效應(yīng)的影響，且樁徑越大，影響越大。

圖7反映的是在其他條件不變的情況下，在低頻范圍內(nèi)樁徑的變化對樁頂復(fù)剛度的影響。從樁頂動剛度曲線可得，在低頻范圍內(nèi)，樁半徑越大，樁頂動剛度越??；隨著激振頻率的增大，樁頂動剛度先增大后減小，且樁半徑越大衰減的越快。從樁頂動阻尼曲線可得，當(dāng)樁半徑不變時，激振頻率越大，樁頂動阻尼越大；對于同一激振頻率，樁半徑越大，樁頂動阻尼越小。

4.3 水的體積分?jǐn)?shù)對樁頂動力響應(yīng)的影響

圖8反映的是在其他參數(shù)都不變的情況下(其中r0=0.8 m)，水的體積分?jǐn)?shù)對樁頂速度響應(yīng)的影響(水的體積分?jǐn)?shù)是指飽和土中水的體積與土體總體積的比值)。從圖中可以看出，飽和土中水的體積分?jǐn)?shù)越大，樁底反射信號幅值越大，樁頂導(dǎo)納的幅值也越大，這是由于水的體積分?jǐn)?shù)提高導(dǎo)致土體剪切模量降低的原因。

圖5 樁的橫向慣性效應(yīng)對樁頂復(fù)剛度的影響(r0=0.4 m)Fig.5 Influence of transversal inertia effect on complex stiffness at top of pile(r0=0.4 m)

圖6 樁的半徑對樁頂速度響應(yīng)的影響Fig.6 Influence of radius on velocity response at top of pile

圖7 樁的半徑對樁頂復(fù)剛度的影響Fig.7 Influence of radius on complex stiffness at top of pile

圖8 水的體積分?jǐn)?shù)對樁頂速度響應(yīng)的影響Fig.8 Influence of water volume fraction on velocity response at top of pile

圖9 水的體積分?jǐn)?shù)對樁頂復(fù)剛度的影響Fig.9 Influence of water volume fraction on complex stiffness at top of pile

圖9反映的是在其他參數(shù)不變時，在低頻階段，飽和土中水的體積分?jǐn)?shù)對樁頂復(fù)剛度的影響。在激振頻率不變時，隨著水的體積分?jǐn)?shù)的增加，樁頂動剛度明顯減小，樁底動阻尼亦減??；在水體積分?jǐn)?shù)不變時，隨著激振頻率的增加，樁頂動剛度先增加后降低，且水的體積分?jǐn)?shù)越大衰減的越快，而樁底動阻尼隨著激振頻率的增加而增加。

5 結(jié) 論

本文分析了均質(zhì)飽和土中考慮橫向慣性效應(yīng)時樁基的縱向振動特性，得出了樁頂導(dǎo)納的解析解和樁頂速度時域響應(yīng)的半解析解，并分析了相關(guān)參數(shù)對樁頂動力響應(yīng)的影響，得到如下結(jié)論：

(1)對于大直徑樁基，橫向慣性效應(yīng)對樁頂動力響應(yīng)有著明顯的影響：橫向慣性效應(yīng)會使樁頂速度時域響應(yīng)曲線的振蕩加劇，樁底反射信號的波谷值明顯減??；在高頻范圍內(nèi)，橫向慣性效應(yīng)會使樁頂復(fù)剛度的峰值減小，但在低頻范圍內(nèi)，可以忽略橫向慣性效應(yīng)對樁頂復(fù)剛度的影響。

(2)忽略橫向慣性效應(yīng)的一維解只適用于細(xì)長樁。

(3)其他條件不變時，樁徑的增大使樁底反射信號明顯增強(qiáng)，反射波曲線的振蕩明顯加劇，樁頂導(dǎo)納曲線的峰值增大，共振頻率提前。在低頻范圍內(nèi)，樁徑的增大會使樁頂動剛度的衰減加快，會使樁頂動阻尼減小。

(4)其他參數(shù)不變時，隨著飽和土中水的體積分?jǐn)?shù)的增大，樁底反射信號幅值和樁頂導(dǎo)納幅值均增大。在低頻范圍內(nèi)，水的體積分?jǐn)?shù)越大，樁頂動剛度越小，且衰減越快，樁頂動阻尼也越小。

[ 1 ] KESHAVARZ A, MALEKZADEH P, HOSSEINI A. Time domain dynamic analysis of floating piles under impact loads[J]. International Journal of Geomechanics, 2016, 17(2): 04016051.

[ 2 ] WANG K, WU W, ZHANG Z, et al. Vertical dynamic response of an inhomogeneous viscoelastic pile[J]. Computers and Geotechnics, 2010, 37(4): 536-544.

[ 3 ] 吳文兵, 王奎華, 竇斌. 任意層地基中黏彈性楔形樁縱向振動響應(yīng)研究[J]. 振動與沖擊, 2013, 32(8): 120-127.

WU Wenbing, WANG Kuihua, DOU Bin. Vertical dynamic response of a tapered pile embedded in layered foundation[J]. Journal of Vibration and Shock, 2013, 32(8): 120-127.

[ 4 ] ZENG X, RAJAPAKSE R. Dynamic axial load transfer from elastic bar to poroelastic medium[J]. Journal of Engineering Mechanics, 1999, 125(9): 1048-1055.

[ 5 ] WANG J H, LI Y L, QING B H, et al. Integral equation method for the dynamic responses of the pile-saturated-soil system to a moving load[J]. Geotechnical Special Publication, 2009, 27(186): 592-599.

[ 6 ] 李強(qiáng), 王奎華, 謝康和. 飽和土中大直徑嵌巖樁縱向振動特性研究[J]. 振動工程學(xué)報, 2006, 18(4): 500-505.

LI Qiang, WANG Kuihua, XIE Kanghe. Dynamic response for vertical vibration of large diameter pile in saturated soil[J]. Journal of Vibration Engineering, 2006, 18(4): 500-505.

[ 7 ] 楊驍, 蔡雪瓊. 考慮橫向效應(yīng)飽和黏彈性土層中樁的縱向振動[J]. 巖土力學(xué), 2011, 32(6): 1857-1863.

YANG Xiao, CAI Xueqiong. Vertical vibration of pile in saturated visco elastic soil layer considering transversal effects[J]. Rock and Soil Mechanics, 2011, 32(6): 1857-1863.

[ 8 ] 楊驍, 唐潔. 成層飽和土中考慮橫向慣性的單樁縱向振動[J]. 巖土力學(xué), 2013, 34(6): 1560-1566.

YANG Xiao, TANG Jie. Vertical vibration of single pile with transversal inertia effect in stratified saturated soil[J]. Rock and Soil Mechanics, 2013, 34(6): 1560-1566.

[ 9 ] 劉林超, 楊驍. 基于薄層法的飽和土樁縱向振動研究[J]. 巖土力學(xué), 2010, 31(1): 92-98.

LIU Linchao, YANG Xiao. Study of longitudinal vibrations of pile in saturated soil based on layer method[J]. Rock and Soil Mechanics, 2010, 31(1): 92-98.

[10] 劉林超, 閆啟方. 飽和土中管樁的縱向振動特性[J]. 水利學(xué)報, 2011, 42(3): 366-372.

LIU Linchao, YAN Qifang. Vertical vibration characteristics of pipe pile in saturated soil[J].Journal of Hydraulic Engineering, 2011, 42(3): 366-372.

[11] 張智卿, 王奎華, 李強(qiáng), 等. 非飽和土中端承樁縱向振動問題簡化解[J]. 工程力學(xué), 2010,27 (5): 159-165.

ZHANG Zhiqing, WANG Kuihua, LI Qiang, et al. Simplified solution for vertical vibration of an end bearing pile embedded in unsaturated soil[J]. Engineering Mechanics, 2010,27 (5): 159-165.

[12] 王星華, 章敏. 非飽和土半空間中單樁豎向振動特性研究[J]. 巖土工程學(xué)報, 2013, 35(7): 1210-1220.

WANG Xinghua, ZHANG Min. Dynamical response of vertically loaded pile in an unsaturated porous media half-space[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(7): 1210-1220.

[13] NOVAK M, ABOUL-ELLA F, NOGAMI T. Dynamic soil reactions for plane strain case[J]. Journal of the Engineering Mechanics Division, 1978, 104(4): 953-959.

[14] 黎正根, 龔育齡. 波在大直徑樁中傳播的三維效應(yīng)現(xiàn)象[J]. 巖石力學(xué)與工程學(xué)報, 1998, 17(4): 434-439.

LI Zhenggen，GONG Yuling. Discussion on three-dimensional effect of wave propagating in the pile with large diameter[J]. Chinese Journal of Rock Mechanics and Engineering，1998，17(4) : 434-439.

[15] LYSMER J, RICHART F E. Dynamic response of footings to vertical loading[J]. Journal of Soil Mechanics & Foundations Div, 1966,92(1): 65-91.

[16] 楊桂通, 張善元. 彈性動力學(xué)[M]. 北京: 中國鐵道出版社, 1988： 240-242.

[17] 王奎華, 應(yīng)宏偉. 廣義 Voigt 土模型條件下樁的縱向振動響應(yīng)與應(yīng)用[J]. 固體力學(xué)學(xué)報, 2003, 24(3): 293-303.

WANG Kuihua, YANG Hongwei. Vertical vibration response and application of pile with general Voigt model[J]. Journal of Solid Mechanics, 2003, 24(3): 293-303.