考慮橫向慣性效應的非飽和土中單樁的豎向動力響應1)

郭鵬飛周順華 楊龍才 肖軍華

(同濟大學道路與交通工程教育部重點實驗室,上海 201804)

考慮橫向慣性效應的非飽和土中單樁的豎向動力響應1)

郭鵬飛2)周順華 楊龍才 肖軍華

(同濟大學道路與交通工程教育部重點實驗室,上海 201804)

基于非飽和土的動力控制方程，考慮橫向慣性效應,建立了三相非飽和介質中嵌巖樁的豎向動力響應連續(xù)介質模型，對樁側非飽和土的動力控制方程進行Laplace變換，在頻域內(nèi)，通過引入勢函數(shù)、算子分解等手段對控制方程進行解析，得到了樁側土體剪應力及豎向振動位移的表達式.結合樁基的豎向振動方程及樁–土接觸面的連續(xù)性條件，使樁土耦合振動系統(tǒng)得以解答，最終在頻域內(nèi)得到了樁頂復剛度、導納、樁–土系統(tǒng)振動位移及應力的解析解，借助Laplace逆變換得到了半正弦激勵載荷下樁頂?shù)乃俣葧r程曲線.最后，通過算例分析驗證了計算結果的準確性，分析了橫向慣性、泊松比、飽和度、長徑比、樁土模量比等因素對樁基動力響應的影響.結果表明：(1)單樁動剛度、阻尼、導納等變量隨頻率變化發(fā)生周期性振蕩，在樁基各階固有頻率處發(fā)生共振；(2)泊松比、飽和度、長徑比、樁土模量比等因素對樁基的動力響應有較大影響，且頻率越大，影響越明顯；(3)泊松比越大，單樁動剛度、阻尼、導納的波動幅值及對應的頻率越小，樁頂時程曲線中的樁底反射信號越弱；(4)飽和度越大，對應各動力響應的波動幅值越大，且樁底反射信號的波峰越大.

非飽和土，算子分解，橫向慣性，樁基動力響應

引言

作為常見的地基結構，樁基具有承載力大、穩(wěn)定性好、沉降量小、工程適應性強等特點，在大跨橋梁、動力廠房、風力發(fā)電、高層建筑等眾多領域得到廣泛應用.樁基在使用期間，除靜力載荷外，樁基往往還要承受上部結構帶來的交通、風、海浪、重型吊車等動力載荷作用[1].對于靜載作用下的樁基，研究多關注于樁基的承載及沉降特性，成果較多且理論完備.而對于動力載荷下的樁基，仍有較多問題亟待解決，如抗震設計中樁基的動剛度、動阻尼取值，動力廠房及大跨度橋梁引起的周邊環(huán)境振動，樁基的完整性檢測及長期沉降計算等問題.樁基的振動理論即求解動力載荷作用下樁基的動力響應，包括樁基的剛度、阻尼、位移幅頻及導納等，其核心是樁基與樁周土的動力相互作用，是解決樁基動力問題的基礎[2].

對此，許多學者對單相彈性介質或飽和兩相土中樁基的動力響應進行了研究.Berger[3]用ALUSH有限元程序進行了地震作用下的樁–土相互作用數(shù)值分析，對樁頂振動位移進行了分析.Hu等[4-5]建立了樁周介質的黏彈性流體多孔介質模型，模型考慮為軸對稱問題，在連續(xù)介質力學框架內(nèi)構建了飽和土中樁–土的耦合系統(tǒng)控制方程，在空間域內(nèi)采用DQM方法對控制方程進行離散，在時域內(nèi)得到了樁及樁周介質振動的數(shù)值解.王奎華等[6]基于連續(xù)介質理論，將樁周土視為飽和兩相介質，構建了單樁豎向振動的虛土樁模型，通過Laplace變換對勢函數(shù)進行求解，得到了有關樁基振動的解析解.

在數(shù)值解方面，常用方法有有限元、邊界元及兩者的耦合等.Blaney[7]采用一致邊界矩陣的有限元公式代表自由場，利用軸對稱有限元對樁的振動特性進行了三維有限元分析，Sen等[8-9]采用邊界元法(BM)，將樁體劃分為許多圓柱形邊界元，得到了群樁豎直振動的位移場，對非均質土中群樁基礎的橫向和縱向動力響應進行分析.Mamoon等[10]則采用有限元、邊界元耦合算法，計算了單相半空間彈性介質中單樁、群樁的豎向動力阻抗.有限元及邊界元方法在處理復雜邊界問題方面具有優(yōu)勢，但其屬于數(shù)值解或半解析解，需要較大的計算量和較長的計算時間，且計算誤差較大.

在解析解方面，對半空間介質中樁基的豎向振動，可歸結為軸對稱問題，即將半空間介質沿徑向進行離散，通過適當?shù)暮喕斑吔鐥l件處理，可以獲得較為完整的解析解.Novak[11]將土視為線性黏彈性體，考慮樁土界面的應力、位移連續(xù)，首先運用連續(xù)介質力學的方法，建立了彈性樁–土作用的二維模型，計算了單樁的豎向動力響應并通過實驗驗證.隨后，Nogami等[12-13]從三維波動方程出發(fā)，通過構造勢函數(shù)，獲得了自由場水平動力阻抗，并在此基礎上結合邊界條件對樁的水平振動方程進行了求解.李強等[14-16]將Novak的連續(xù)介質力學方法拓展到飽和土領域，引入Biot飽和土的動力方程，樁周視為飽和土，樁假定為一維Bernoulli彈性桿，通過邊界條件的處理求得了飽和土層的豎向動力阻抗，并得到了單樁豎向動力響應的解析解.此外,吳文兵等[17-18]在連續(xù)介質力學的基礎上發(fā)展了相應的虛土樁模型，分析了變截面樁的豎向動力響應，胡昌斌等[19-20]采用該方法對樁基的扭轉問題進行了研究，陸建飛[21]基于Muki等[22]提出的虛土樁法，求解了飽和地基中單樁的水平動力響應.

通過上述分析，早期的研究一般將樁周土假定為單相的彈性連續(xù)介質或飽和的兩相多孔介質.但由于非飽和土中樁基的動力響應問題的復雜性，有關研究成果尚不多見.而非飽和土在地球上分布廣泛，且性質與飽和土存在差異，因此有關非飽和土的研究得到了高度重視.目前考慮各相相互運動與作用的三相非飽和土的動力控制方程已經(jīng)建立[23]，文獻[24]對非飽和土中樁基的動力響應進行了初步研究，但其忽略了振動中樁基的徑向振動位移，與實際情況存在偏差.本文結合非飽和土的動力控制方程，將樁周介質視為非飽和土，樁視為Rayleigh-Love桿，建立相應的連續(xù)介質模型，通過三角函數(shù)和Bessel函數(shù)對非飽和土的控制方程進行離散，結合樁土界面的連續(xù)性條件，得到樁土系統(tǒng)的動力響應解.

1 樁土相互作用的連續(xù)介質模型

如圖1所示，非飽和土中的樁頂作用有激勵載荷F(t)，樁長為H，樁基半徑r0，樁周土為彈性均勻的三相非飽和介質，土層表面自由，正應力和剪應力均為零，土層底部為基巖剛性支撐，樁土緊密接觸，樁與樁側土應力與位移連續(xù)，樁土變形為小變形.

圖1 樁–土相互作用模型Fig.1 The pile-soil model

傳統(tǒng)的樁基縱向振動理論將樁視為一維的彈性桿，即Bernoulli桿，不考慮樁身材料的側向振動，這與實際情況存在明顯偏差.若將樁視為嚴格的三維柱體則又太過復雜，Rayleigh和Love從能量原理出發(fā)，對三維桿進行簡化修正，得到了較為實用的Rayleigh-Love理論[25-26]，根據(jù)此理論,樁的豎向振動方程為

其中，Ep為樁身的彈性模量，ρp為樁身材料密度，υp為樁身材料的泊松比，f(z)為樁側阻力,wp為樁身豎向振動位移.

地表透水透氣、正應力為零

其中，下標r代表徑向，z代表豎向，τzr,σz則分別為樁側土體的剪應力和正應力，e,ε,ξ分別代表土骨架、氣相、液相的體應變，pl為樁側孔隙液壓力，pg為孔隙氣壓力.

剛性基底底部豎向位移為零

樁土完全接觸，無相對滑移脫離且邊界不透水

接觸面土體剪應力與樁身摩阻力相等

樁頂?shù)倪吔?/p>

樁端豎向振動位移為零

式中，ui(i=r,z)為固相位移，wi,vi(i=r,z)分別為液相、氣相與土骨架間的相對位移.

2 非飽和土的動力控制方程求解

引入徐明江[23]建立的非飽和土動力控制方程，基于連續(xù)介質力學理論，結合空間平均化方法，采用非飽和土的水土特征曲線V-G模型推導了非飽和土的實用波動方程.該方程基于Mualem理論考慮孔隙流體的滲透系數(shù).引入方程并將其轉化為u,v,w的矢量形式有

滲流連續(xù)性方程為

方程中ρ為非飽和多孔介質的平均密度

v為孔隙水與固體骨架的相對位移

w氣體與固體骨架間的相對位移

其余各變量及意義見附錄.

將矢量控制方程(10)展開為柱坐標系下的形式，并將樁–土耦合作用模型考慮為軸對稱問題并進行無量綱化，得到的控制方程為

其中,拉普拉斯算子為

固體骨架、液體及氣體的體應變分別為

王星華等[24]曾采用直接對方程(12)進行變換處理的方法對非飽和土的動力控制方程進行解析，這里將采用勢函數(shù)分解的方法進行解析.首先引入相應的勢函數(shù)，將固體骨架位移和液體、氣體與固體骨架的相對位移寫成勢函數(shù)的表達式

式中，φs,ψs,φl,ψl,φg,ψg分別為固、液、氣三相的位移勢，且都是坐標r,z與時間t的函數(shù)，將勢函數(shù)(15)代入方程(12)且對方程進行Laplace變換后得到

將上述方程進行無量綱化，并寫為矩陣的形式后有

若上述微分方程組(17a)有非零解，則需滿足方程組的系數(shù)行列式矩陣等于零的條件，則有

方程(18)中的各系數(shù)分別為

對于方程(18a)，先將方程分解為

由于R1(r),Z1(z)互不相關，則方程(21)左右兩側均等于一常數(shù).因此其可以轉化為一個零階虛宗量Bessel函數(shù)(22a)，與一個二階常系數(shù)微分方程(22b)

其中，I0(cr)和K0(cr)分別為第一類和第二類零階虛宗量Bessel函數(shù)，E1,F1,G1,H1為待定系數(shù).

基于同樣的求解方法，則由方程 (20b)、方程(20c)可得

由徑向無限遠處土體的位移、應力及應變?yōu)榱?以及零階第一類Bessel函數(shù)的性質知E1,E2,E3,E4均等于零.

因此有

根據(jù)同樣的方法可以得到

將方程(24a)～方程(24c)代回方程組(17a)可得到關系式

將方程(24d)～方程(24f)代回方程組(17b)可得到關系式

根據(jù)方程(25)、方程(24)、方程(14)可以得到各相介質的豎向和徑向位移表達式以及土體骨架正應力、剪應力的表達式

根據(jù)邊界條件方程(2)可知地表正應力為零且K0(c11r),K0(c21r),K0(c31r)線性不相關，得到

為了進一步對表達式 (26)進行化簡，將方程(26d)與方程(26f)代入樁–土邊界條件方程(6)可以得到

方程中的G11,G12,G13,G14均為待求變量，通過推導可以將G13,G14表達為

其中θ=H/r0為樁基的長徑比.

為了表述方便，令

則進一步可以得到

至此，通過對非飽和土動力控制方程的解析，得到了樁側各相介質的振動位移及應力表達式(見式(32)～式(34)，表達式中仍有G11,G12兩個未知量尚未確定，下面將借助樁–土界面的連續(xù)性條件，通過對樁基豎向振動方程的求解，最終確定各未知量的表達式并求得系統(tǒng)的動力響應.

3 樁–土耦合系統(tǒng)求解

將樁基視為Rayleigh-Love桿，對方程(1)進行Laplace變換并通過無量綱化后得到

根據(jù)邊界條件(6)，將(32)代入上述二階常系數(shù)微分方程，得到樁基豎向振動位移的解為

此時，樁基的橫向振動位移為

根據(jù)樁土接觸面位置的位移連續(xù)性條件，將方程(33)、式(34)、式(36)、式(37)代入邊界條件(4)可以得到如下等式

此外將方程(36)、方程(37)代入樁頂與樁底的邊界條件方程(8)、方程(9)可以得到

令s=iω，則根據(jù)樁基的樁頂復剛度的定義有

Z′(ω)為無量綱復阻抗

樁頂無量綱復剛度為

式中，Kd代表樁基的復剛度，K為動剛度的實部，反映樁基抵抗變形的能力，虛部C為樁基的阻尼，反映樁–土系統(tǒng)的能量耗散特性.

樁頂位移頻域響應函數(shù)為

位移頻域響應函數(shù)的相位差為

樁頂導納(速度頻域響應函數(shù))為

此外，通過傅里葉逆變換可以得到樁基的單位脈沖時域響應函數(shù)

當在樁頂施加半正弦的脈沖激振力時，即樁頂載荷F(t)=Fmaxsin(πt/t0)，其中t∈(0,t0)，t0為載荷持續(xù)時間.通過傅里葉逆變換和卷積定律可以得到樁頂?shù)臒o量綱時域響應函數(shù)

方程中的a0為無量綱頻率，a0=ωTp，Tp為波從樁頂傳播到樁底所用時間，為無量綱脈沖寬度因子，無量綱時間

本節(jié)求得的各物理量表達式均為無量綱化后的結果，實際應用時只需通過相應的轉化關系將無量綱量轉化為實際結果即可.

4 豎向振動特性分析

樁基的豎向振動特性主要包括動剛度、阻尼及導納等[25-27]，上文已根據(jù)連續(xù)介質模型得到了各動力響應的解析表達式.通過推導可以發(fā)現(xiàn)，樁基的動力響應與載荷特性、材料性質及幾何尺寸等多種因素有關，例如飽和度Sr、樁土模量比Ep/μ、長徑比H/r0等因素都會對樁基的振動特性產(chǎn)生較大影響[26]，本節(jié)結合文獻[24]，分別計算在考慮和不考慮橫向慣性效應兩種情況下的樁基動力響應，分析各樁土參數(shù)對單樁動力響應的影響.

選用文獻[23]中給出的有關非飽和土的材料參數(shù)，樁身材料的物理力學參數(shù)按文獻[24]進行選取，如表1所示.由于結果較多，本文僅對單樁的剛度、阻尼、導納及樁頂速度時程進行分析.結果中的剛度、阻尼在樁基的抗震及動力機器基礎設計中具有重要作用，而導納和速度時程則是樁基小應變及樁身完整性檢測的理論基礎.圖中的橫坐標a0=ωr0/cp為載荷的無量綱頻率，其中代表樁身材料的等效振動波速[28-29].

表1 樁及樁周土物理力學性質Table 1 Parameters of pile and soil

首先分析樁身泊松比對動力響應的影響并驗證結果的準確性，作為對照，計算時取Sr=0.6,κ=1.0×106，H=20 m，其中不考慮橫向慣性效應為文獻[24]的結果，圖2對其進行了進一步的無量綱化處理，從圖2可以發(fā)現(xiàn)，υp=0時，本文結果能完整地退化為不考慮橫向慣性效應時的樁基振動解，進而驗證了本文結果的準確性.圖2(a)～圖2(c)分別為樁基的無量綱動剛度、阻尼及速度幅頻(導納)曲線，圖中阻尼由徑向輻射阻尼與各相介質間相互運動產(chǎn)生的阻尼組成，可以發(fā)現(xiàn)，樁基各變量隨頻率變化發(fā)生周期性振蕩，在各固有頻率處發(fā)生共振.頻率較小時，樁基動剛度近似等于靜剛度，動阻尼近似為零，當頻率增大時，動剛度波動幅度不斷增大，動阻尼的波峰和波谷也不斷增大，而導納波動幅度則基本不變.圖2(d)代表半正弦激勵載荷下樁頂?shù)乃俣葧r程曲線，計算時無量綱脈沖寬度因子取0.25，可以發(fā)現(xiàn)，載荷作用下，樁頂無量綱速度近似為1，一定間隔后樁頂收到樁底的反射信號，且隨著能量耗散，信號幅值不斷減小.

混凝土泊松比對樁基動力響應有明顯影響，圖2(a)～圖2(c)可以發(fā)現(xiàn)，當考慮橫向慣性效應時，隨頻率的增大，剛度、阻尼及導納的波動幅度不斷減小，同時各幅值對應的載荷頻率(自振頻率)也不斷減小，且泊松比越大，規(guī)律越明顯；圖2(c)中，隨頻率的增加，速度幅值不斷減小，即頻率越大，樁身振動的能量耗散越大，考慮泊松比能更真實的反應樁身振動過程中的能量耗散.在時程曲線圖2(d)中，考慮橫向慣性效應時的樁底反射值明顯小于不考慮橫向慣性效應的反射值，且樁底反射對應的時間也有所偏差，此外，考慮橫向慣性效應時，樁底反射值出現(xiàn)了明顯的振蕩現(xiàn)象，且泊松比越大振蕩越明顯.

圖3為飽和度對樁頂動力響應的影響，分析時Sr分別取0.3,0.9，υp=0.25，從圖3(a)～圖3(c)中可以發(fā)現(xiàn)，與不考慮橫向慣性效應規(guī)律相似，考慮橫向慣性效應時，飽和度對樁基系統(tǒng)的固有頻率影響不大，即飽和度發(fā)生變化時，樁基剛度、阻尼、導納達到波動幅值時的頻率沒有變化，這主要是因為系統(tǒng)的固有振型主要由樁基的物理性質決定，飽和度變化對系統(tǒng)的參振質量影響不大，因此主頻變化不大.飽和度對各動力響應的振動幅值有明顯影響，且頻率越大，飽和度對其影響程度也越大.圖3(d)為飽和度對樁頂速度時程的影響，可以發(fā)現(xiàn)，飽和度對樁頂速度幅值有明顯影響，飽和度越大，樁底反射峰值越小，飽和度對樁底反射波出現(xiàn)的時間沒有影響.

圖2 樁身泊松比對動力響應的影響Fig.2 The e ff ect of pile Poisson ratio on dynamic response

圖2 樁身泊松比對動力響應的影響(續(xù))Fig.2 The e ff ect of pile Poisson ratio on dynamic response(continued)

圖3 飽和度對單樁動力響應的影響Fig.3 The e ff ect of saturation on dynamic response

圖4為長徑比對樁基動力響應的影響.計算時Sr=0.6，保持樁基截面不變，θ分別取10,20，可以發(fā)現(xiàn)無論考慮橫向慣性與否，長徑比都是影響樁基動力響應的最主要因素，長徑比越大，剛度、阻尼及導納的震蕩幅值越小，但對應的各階主頻卻逐漸減小.長徑比越大，樁基的剛度幅值越小，其主要原因是，本文的研究對象為嵌巖樁，樁基越短，基巖埋深越淺，則樁基的剛度越大；長徑比越大，阻尼的波動幅度減小，但樁基的整體阻尼卻不斷增大，樁基的阻尼主要由樁側土提供，即樁側土越厚，樁基的阻尼越大.圖4(d)為樁頂速度時程曲線，可以發(fā)現(xiàn)，樁基越長，樁底反射越不明顯，長徑比差別越大，樁底反射幅值的差距越大，長徑比為60時，樁底第二次反射已非常微弱，即長徑比越大，能量在樁身傳播的能量耗散越大.

圖5為樁–土模量比對樁基動力響應的影響，計算時θ=20，Sr=0.6，保持樁身模量不變改變土的剪切模量，參照工程實際，模量比分別取1 000,500，可見，隨著模量比的增大，剛度、阻尼等的波動幅值不斷增大，而波長(固有頻率)卻逐漸減小，模量比對各變量的波動幅度有較大影響，但對波長影響較小.可見當樁側土的模量較大時，有利于提高樁基的抗震性能.從圖5(d)中可以發(fā)現(xiàn)，模量比越大，能量在樁身傳播工程中耗散越多，輻射出去的能量越大，樁底反射越微弱.

圖4 長徑比對單樁動力響應的影響Fig.4 The e ff ect of pile length-diameter ratio on dynamic response

圖5 樁土模量比對單樁動力響應的影響Fig.5 The e ff ect of pile-soil modulus ratio on dynamic response

圖5 樁土模量比對單樁動力響應的影響(續(xù))Fig.5 The e ff ect of pile-soil modulus ratio on dynamic response(continued)

5 結論

本文將樁周介質視為三相非飽和土，假定樁為Rayleigh-Love桿，考慮樁基的橫向慣性效應，根據(jù)樁土接觸的連續(xù)性條件，建立了非飽和土中嵌巖樁豎向動力響應的理論求解模型.通過勢函數(shù)、算子分解等手段對模型進行解析，得到了周期載荷作用下，單樁的動剛度、阻尼、導納等解析表達式.借助Laplace逆變換得到了樁頂?shù)乃俣葧r程曲線.通過案例分析得到如下結論：

(1)單樁動剛度、阻尼、導納等變量隨頻率發(fā)生周期性振蕩，在樁基各階固有頻率處發(fā)生共振.當頻率較小時，動剛度近似等于靜剛度，動阻尼近似為零，當載荷頻率增大時，動剛度波動幅度不斷增大，動阻尼整體不斷增大，而導納的波動幅度則基本穩(wěn)定.

(2)混凝土泊松比為零時，本文計算結果與文獻[24]相吻合，驗證了本文計算結果的準確性；考慮樁基的橫向慣性效應時，隨頻率的增大，剛度、阻尼及導納的波動幅度不斷減小，各幅值的對應頻率(自振頻率)也不斷減小，泊松比越大，規(guī)律越明顯；泊松比對速度時程曲線影響較大，泊松比越大，樁底反射越弱，波形的振蕩越明顯，反射波的對應時間也會有所偏差.

(3)飽和度對樁基的動力響應幅值有明顯影響，且頻率越大，影響越明顯；飽和度對樁基的固有頻率影響較??；對樁頂?shù)乃俣葧r程影響較大，飽和度越小，速度時程的樁底反射越弱.

(4)與不考慮橫向慣性效應時規(guī)律相似，長徑比和樁–土模量比對單樁的動力響應影響較大，尤其長徑比，不僅對動力響應的波動幅值產(chǎn)生明顯影響，當長徑比改變時，系統(tǒng)的固有振動頻率也會發(fā)生明顯變化.

1 郭鵬飛,王旭,楊龍才等.長期豎向循環(huán)載荷作用下黃土中單樁沉降特性模型試驗研究.巖土工程學報，2015,37(3):551-558(Guo Pengfei,Wang Xu,Yang Longcai,et al.Model tests on settlement of single pile in loess under long-term axial cyclic loading.Chinese Journal of Geotechnical Engineering，2015,37(3):551-558 (in Chinese))

2 楊龍才,郭慶海,周順華等.高速鐵路橋樁在軸向循環(huán)載荷長期作用下的承載和變形特性試驗研究.巖石力學與工程學報,2005, 24(13):2362-2368(Yang Longcai,Guo Qinghai,Zhou Shunhua, et al.Dynamic behaviors of pile foundation of high-speed railway bridge under long-term cyclic loading in soft soil.Chinese Journal of Rock Mechanics and Engineering,2005,24(13):2362-2368(in Chinese))

3 Berger E.Simplifie method for evaluating soil-pile-structure interaction e ff ects//Proceedings of 9th O ff-shore Technology Conference,1977:589-601

4 Hu YJ,Zhu YY,Cheng CJ.DQEM for large deformation analysis of structures with discontinuity conditions and initial displacements.Engineering Structures,2008,30:1473-1487

5 胡育佳.樁基非線性靜動力學特性研究.[博士論文].上海：上海大學,2008(Hu Yujia.The research on nonlinear static and dynamic characteristics of piles.[PhD Thesis].Shanghai:Shanghai University,2008(in Chinese))

6 王奎華,王寧,劉凱等.三維軸對稱條件下基于虛土樁法的單樁縱向振動分析.巖土工程學報,2012,34(5):885-892(Wang Kuihua,Wang Ning,Liu Kai,et al.Longitudinal vibration of piles in 3D axisymmetric soil based on fictitiou soil pile method.Chinese Journal of Geotechnical Engineering,2012,34(5):885-892(in Chinese))

7 Blaney GW.Dynamic sti ff ness of piles//Proceeding of 2nd International Conference on Numerical Method in Geomechanics,1976,Ⅱ:1001-1012

8 Sen R,Kausel E,Banerjee RK.Dynamic analysis of piles and pile groups embedded in homogenous soils.Earthquake Engineering and Structural Dynamics,1985:13(1):53-65

9 Sen R,Kausel E,Banerjee RK.Dynamic analysis of piles and pile groups embedded in non-homogenous soils.Int.J.Num.Anal. Meth.Geomech,1985,9:507-524

10 Mamoon SM,Kaynia AM,Banerjee PK.Frequency domain dynamic analysis of piles and pile groups.J Eng Mech—ASCE,1990: 116(10):2237-2257

11 Novak M.Dynamic sti ff ness and damping of piles.Canadian Geotechnical Journal,1974,11:574-598

12 Nogami T,Novak M.Resistance of soil to a horizontally vibrating pile.Earthquake Engineering and Structural Dynamics,1977,5: 249-261

13 Novak M,Nogami T.Soil-pile interaction in horizontal vibration.Earthquake Engineering and Structural Dynamics,1977,5:263-281

14 李強,王奎華,謝康和.飽和土樁縱向振動引起土層復阻抗分析研究.巖土工程學報,2004,26(5):679-683(Li Qiang,Wang Kuihua,Xe Kanghe.Study on resistance factor of saturated soil caused by longitudinal vibration of pile.Chinese Journal of Geotechnical Engineering,2004,26(5):679-683(in Chinese))

15 李強,王奎華,謝康和.飽和土中端承樁縱向振動特性研究.力學學報,2004,36(4):435-442(Li Qiang,Wang Kuihua,Xie Kanghe. Vertical vibration of an end bearing pile embedded in saturated soil.Acta Mech Sinica,2004,36(4):435-442(in Chinese))

16 李強,鄭輝,王奎華.飽和土中摩擦樁豎向振動解析解及應用,工程力學,2011,28(1):157-170(Li Qiang,Zheng Hui,Wang Kuihua.Analytical solution and its application of vertical vibration of a friction pile in saturated soil.Engineering Mechanics,2011,28(1): 157-170(in Chinese))

17 吳文兵,王奎華,張智卿等.半空間地基中虛土樁模型的精度分析及應用.應用基礎與工程科學學報,2012,20(1):121-129(Wu Wenbing,Wang Kuihua,Zhang Zhiqing,et al.Accuracy and application of virtual soil pile model in half space foundation.Journal of Basic Science and Engineering,2012,20(1):121-129(in Chinese))

18 吳文兵,蔣國盛,王奎華等.土塞效應對管樁縱向動力特性的影響研究.巖土工程學報,2014,36(6):1129-1141(Wu Wenbing,Jiang Guosheng,Wang Kuihua,et al.Influenc of soil plug e ff ect on vertical dynamic response of pipe piles.ChineseJournalof Geotechnical Engineering,2014,36(6):1129-1141(in Chinese))

19 胡昌斌,張濤.考慮樁土耦合作用時樁基扭轉振動特性研究.工程力學,2007,24(3):147-153(Hu Changbin,Zhang Tao.Soil-pile interaction in torsional vibrations of a pile in viscious damping soil layer.Engineering Mechanics,2007,24(3):147-153(in Chinese))

20 胡昌斌,張濤.樁土耦合扭轉振動理論研究與特性分析.巖土工程學報,2007,29(2):184-190(Hu Changbin,Zhang Tao.Study on soil-pile interaction in torsional vibrations.Chinese Journal of Geotechnical Engineering,2007,29(2):184-190(in Chinese))

21 陸建飛.頻域內(nèi)半空間飽和土中水平受荷樁的動力分析.巖石力學與工程學報,2002,21(4):577-581(Lu Jianfei.Dynamic frequency domain analysis on laterally loaded pile in saturated soil.Chinese Journal of Rock Mechanics and Engineering,2002,21(4): 577-581(in Chinese))

22 Muki R,Sternberg E.Elastostatic load-transfer to a half-space from apartially embedded axially loaded rod.J Eng Mech,1970,125(9): 1048-1055

23 徐明江.非飽和土地基與基礎的動力響應研究.[博士論文].廣州：華南理工大學,2010(Xu Mingjiang.Investigation on dynamic response of unsaturated soils and foundation.[PhD Thesis]. Guangzhou:South China University of Technology,2010(in Chinese))

24 王星華,章敏.非飽和土半空間中單樁豎向振動特性研究.巖土工程學報,2013,35(7):1210-1220(Wang Xinghua,Zhang Min. Dynamical response of vertically loaded pile in half-space of unsaturated porous media.Chinese Journal of Geotechnical Engineering, 2013,35(7):1210-1220(in Chinese))

25 王寧,王奎華.樁底土的成層性對樁體縱向剛度的影響.巖石力學與工程學報,2013,32(5):1042-1048(Wang Ning,Wang Kuihua.Influenc of layering of stratum under pile tip on pile longitudinal sti ff ness.Chinese Journal of Rock Mechanics and Engineering, 2013,32(5):1042-1048(in Chinese))

26 李強.考慮三維波動的飽和土中樁縱向耦合振動理論.[博士論文].杭州：浙江大學，2004(Li Qiang.Study on pile longitudinal vibrationinsaturatedsoilconsideringthreedimensionalwavee ff ect. [PhD Thesis].Hangzhou:Zhejiang University,2004(in Chinese))

27 呂述暉,王奎華,張鵬.梁–樁–土豎向耦合振動特性分析.力學學報,2015,47(1):169-173(LShuhui,Wang kuihua,Zhang Peng. Vertical coupling virbration of beam-pile-soil system.Chinese Journal of Theoretical and Applied Mechanics,2015,47(1):169-173(in Chinese))

28 王桂敏,李強,王奎華.單層飽和土中樁豎向振動簡化模型及其解析解.巖石力學與工程學報,2006,25(2):4233-4240(Wang Guimin,Li Qiang,Wang Kuihua.Simplifie model for vertical vibration of pile in single-layer saturated soil and its analytical solution.Chinese Journal of Rock Mechanics and Engineering,2006, 25(2):4233-4240(in Chinese))

29 周香蓮,周光明,王建華.水平簡諧載荷作用下飽和土中群樁的動力反應.巖石力學與工程學報,2005,24(8):1433-1438(Zhou Xianglian,Zhou Guangming,Wang Jianhua.Dynamic response of pile group in saturated soil subjected to horizontal loads.Chinese Journal of Rock Mechanics and Engineering,2005,24(8):1433-1438(in Chinese))

附錄

Sr和Sw0分別代表流體飽和度和束縛飽和度；γ為有效應力系數(shù)，通常等于Sr，n為土體孔隙率，λ,μ分別為土骨架的拉梅常數(shù)，ηl和ηg分別代表流體和空氣的黏滯系數(shù)；krl和krg分別代表流體和空氣的滲透系數(shù)；κ為土的固有滲透系數(shù)，Ks,Kl,Kg,Kb分別是土顆粒、流體、空氣及土骨架的體積壓縮模量；α,m,d分別為擬合參數(shù)，符號s,l,g分別表示非飽和土體中的固、液、氣三相.

ANALYTICAL SOLUTION OF THE VERTICAL DYNAMIC RESPONSE OF ROCK-SOCKED PILE CONSIDERING TRANSVERSE INERTIAL EFFECT IN UNSATURATED SOIL1)

Guo Pengfei2)Zhou Shunhua Yang Longcai Xiao Junhua
(Key Laboratory of Road and Traffic Engineering of the Ministry of Education,Tongji University,Shanghai201804,China)

Based on the dynamic governing equations of unsaturated soil,a continuous medium model regarding vertical dynamic response of rock-socked pile in three-phase unsaturated medium was established in this paper,taking the transverse inertial e ff ect into consideration.The Laplace transformation is then used to solve the dynamic governing equations. In frequency domain,potential function and operator decomposition methods are also used to resolve the governing equations,thus obtained the expressions of soil vertical vibration displacement and shear stress.To solve the soil-pile coupling system,the vertical vibration equation of pile foundation and the continuity conditions are both adopted.Eventually,complex sti ff ness and admittance of the pile butt,vertical vibration displacement and shear stress of the model are obtained in frequency domain.And the travel curve of vibration velocity under half sine excitation load has also been deductedwith the help of Laplace inverse transformation.The model is verifie through the comparison with the saturated models. Finally,a case study of dynamic response of a pile in unsaturated soil is presented.The influence of transverse inertial effect,Poisson’s ratio,saturation,pile length-diameter ratio and pile-soil modulus ratio on response of pile are investigated. The results show that:(1)Dynamic sti ff ness,damping and admittance oscillate with the frequency,and the pile resonance occur at the natural frequency of pile.(2)Dynamic response of the pile is sensitive to Poisson’s ratio,saturation,pile length-diameter ratio and pile-soil modulus,while the sensitivity with higher frequency.(3)Larger Poisson’s ratio causes smaller amplitude of dynamic sti ff ness,damping and admittance,as well as smaller reflectio signal at bottom of the pile in its speed history curve.(4)Larger saturation leads to larger amplitude of pile response and peak value of reflectio signal from the pile bottom.

unsaturated soil,operator decomposition,transverse inertial e ff ect,response of pile

TU45

A

10.6052/0459-1879-16-286

2016–10–17收稿，2017–01–03錄用,2017–01–04網(wǎng)絡版發(fā)表.

1)國家自然科學基金資助項目(51478353,41472247).

2)郭鵬飛，博士研究生，主要研究方向：樁基動力學.E-mail：tjdxgpf@163.com

郭鵬飛，周順華，楊龍才，肖軍華.考慮橫向慣性效應的非飽和土中單樁的豎向動力響應.力學學報,2017,49(2):344-358

Guo Pengfei,Zhou Shunhua,Yang Longcai,Xiao Junhua.Analytical solution of the vertical dynamic response of rock-socked pile considering transverse inertial e ff ect in unsaturated soil.Chinese Journal of Theoretical and Applied Mechanics,2017,49(2):344-358