崔春義 孟坤 武亞軍

摘要： 基于復(fù)剛度傳遞徑向多圈層平面應(yīng)變模型，采用黏性阻尼模型描述樁周土材料阻尼，考慮地基土縱向成層特性，建立徑向非均質(zhì)、縱向成層土中變截面和變模量管樁縱向振動(dòng)簡(jiǎn)化分析模型，采用Laplace變換和復(fù)剛度傳遞方法，遞推得出樁周、樁芯土體與管樁樁體界面處復(fù)剛度，進(jìn)而利用樁土完全耦合條件和阻抗函數(shù)傳遞性，推導(dǎo)得到變模量和變截面管樁樁頂動(dòng)力阻抗解析解答，并與已有相關(guān)解析解進(jìn)行退化驗(yàn)證其合理性；通過(guò)進(jìn)一步參數(shù)化分析探討了樁身變模量段、縮頸、擴(kuò)徑及變模量、變截面位置對(duì)管樁縱向振動(dòng)動(dòng)力響應(yīng)特性的影響規(guī)律，可為具體工程實(shí)踐提供理論指導(dǎo)和參考作用。

關(guān)鍵詞： 管樁； 動(dòng)力響應(yīng)； 復(fù)剛度傳遞模型； 徑向非均質(zhì)； 黏性阻尼

中圖分類號(hào)：TU473.1

文獻(xiàn)標(biāo)志碼： A

文章編號(hào)： 1004-4523（2018）04-0707 11

DOI：10.16385/j.cnki.issn.1004-4523.2018.04.019

引 言

縱向成層土中的樁基相互作用體系耦合振動(dòng)特性研究是樁基抗震、防震設(shè)計(jì)及樁基動(dòng)力檢測(cè)等工程技術(shù)領(lǐng)域的理論基礎(chǔ)，一直以來(lái)亦是土動(dòng)力學(xué)、巖土工程和結(jié)構(gòu)-地基相互作用領(lǐng)域的交叉熱點(diǎn)問(wèn)題[1-2]。在樁基施工過(guò)程中，由于擠土、松弛以及其他擾動(dòng)因素的影響，使得樁周土體沿樁基徑向存在一定不均勻性，即徑向非均質(zhì)效應(yīng)[3-4]。為考慮此種徑向非均質(zhì)效應(yīng)，國(guó)內(nèi)外諸多學(xué)者展開了較多的相關(guān)研究工作。首先，Novak等[5-6]基于滯回阻尼平面應(yīng)變模型將樁周土體分為單層內(nèi)部擾動(dòng)區(qū)域和外部半無(wú)限大未擾動(dòng)區(qū)域，假設(shè)內(nèi)部區(qū)域土體質(zhì)量為零，首次推導(dǎo)得出了能簡(jiǎn)化考慮單層內(nèi)部區(qū)域土體軟化效應(yīng)的樁基縱向和扭轉(zhuǎn)阻抗解析解答，初步分析了樁周土體內(nèi)外區(qū)域非均質(zhì)性對(duì)樁-土耦合振動(dòng)特性的影響規(guī)律。在此基礎(chǔ)上，Veletsos等[7-8]通過(guò)進(jìn)一步考慮樁周土單層內(nèi)部擾動(dòng)區(qū)域土體質(zhì)量慣性效應(yīng)，發(fā)展得出了徑向非均質(zhì)性樁周土中樁基縱向阻抗解析解答，并分析了單層內(nèi)部擾動(dòng)區(qū)域土體軟化程度對(duì)樁基振動(dòng)特性的影響規(guī)律。Nogami等[9-10]基于平面應(yīng)變模型計(jì)算得出遠(yuǎn)場(chǎng)區(qū)域土體復(fù)剛度對(duì)應(yīng)的彈簧和阻尼系數(shù)，近場(chǎng)采用非線性Winkler模型，兩者結(jié)合建立了樁-土耦合系統(tǒng)，在此基礎(chǔ)上進(jìn)行了考慮樁周土徑向非均質(zhì)性的樁基縱向振動(dòng)特性分析。Doston等[11]通過(guò)單純定義樁周土內(nèi)部區(qū)域（單層）土體剪切模量呈指數(shù)變化函數(shù)，利用平面應(yīng)變理論，給出了樁周徑向非均質(zhì)土體中樁基縱向和扭轉(zhuǎn)阻抗的解析表達(dá)式。Vaziri等[12]通過(guò)定義內(nèi)部區(qū)域土體剪切模量為隨徑向位移呈拋物線變化函數(shù)，且與外部區(qū)域剪切模量連續(xù)，進(jìn)行了考慮樁周土徑向非均質(zhì)性的樁基縱向振動(dòng)特性分析。Han等[13-14]基于土體剪切模量為隨徑向位移拋物線函數(shù)來(lái)描述土體軟化和硬化，分析了樁周土體徑向軟化和硬化非均質(zhì)性對(duì)樁土耦合振動(dòng)特性的影響規(guī)律。EI Naggar[15]基于平面應(yīng)變模型求得遠(yuǎn)場(chǎng)區(qū)域復(fù)剛度（即看成彈簧和阻尼器組成），將近場(chǎng)區(qū)域簡(jiǎn)化為彈簧、阻尼器及滑移塊體系，建立了能考慮樁周的非均勻性的樁土耦合振動(dòng)模型。進(jìn)一步地，EI Naggar[16]將內(nèi)部區(qū)域沿徑向分為N個(gè)圈層，通過(guò)平面應(yīng)變模型求得的每圈層復(fù)剛度串聯(lián)建立內(nèi)部區(qū)域復(fù)剛度體系再與外部區(qū)域采用平面應(yīng)變模型得到的復(fù)剛度串聯(lián)，來(lái)考慮樁周土的徑向不均勻性，進(jìn)行了考慮樁周土徑向非均質(zhì)性的樁土耦合縱向振動(dòng)特性分析。王奎華和楊冬英等[17-18]指出EI Naggar模型受樁周土剪切模量影響很小，與實(shí)測(cè)曲線和已有理論結(jié)果不符，有必要對(duì)其進(jìn)行修正。同時(shí)，Wang和Yang等[19-20]提出了基于滯回阻尼模型的徑向復(fù)剛度傳遞的多圈層平面應(yīng)變模型，通過(guò)圈層間復(fù)剛度遞推求得樁周土體對(duì)樁體作用的復(fù)剛度，在此基礎(chǔ)上建立了能考慮樁周土徑向非均質(zhì)性的樁土耦合振動(dòng)系統(tǒng)，進(jìn)行了樁頂振動(dòng)特性的頻域和時(shí)域響應(yīng)分析。樁基縱向振動(dòng)理論應(yīng)用到工程檢測(cè)中重點(diǎn)是針對(duì)于非完整樁，基于此，王奎華等[21]、王宏志等[22]、劉東甲等[23]、馮世進(jìn)等[24]研究了縱向成層土中變模量、變截面實(shí)心樁的動(dòng)力響應(yīng)，使得低應(yīng)變動(dòng)測(cè)方法具有更完善的理論指導(dǎo)。

以上研究均是針對(duì)實(shí)心樁展開，而對(duì)于大直徑管樁，由于樁芯土的存在，必然使得其與實(shí)心樁的振動(dòng)特性存在差異。丁選明等[25]和鄭長(zhǎng)杰等[26]同時(shí)考慮樁周土和樁芯土，對(duì)徑向均質(zhì)土中完整管樁振動(dòng)特性進(jìn)行求解，并與實(shí)心樁結(jié)果進(jìn)行對(duì)比，說(shuō)明了在豎向荷載作用下管樁表現(xiàn)出與實(shí)心樁動(dòng)力特性的不同。同樣的管樁在施工過(guò)程中同樣會(huì)引起土體擾動(dòng)，Li等[27]同時(shí)考慮土層縱向和徑向非均質(zhì)，土體采用滯回阻尼模型，基于多圈層復(fù)剛度傳遞模型，對(duì)完整管樁-土耦合縱向振動(dòng)響應(yīng)特性進(jìn)行求解分析。

不難看出，考慮樁身非完整性的研究成果大多集中于實(shí)體樁，同時(shí)在考慮土體徑向非均質(zhì)效應(yīng)時(shí)均假定土體材料阻尼為滯回阻尼。而對(duì)非諧和激振問(wèn)題特別是瞬態(tài)激振條件下樁體時(shí)域振動(dòng)響應(yīng)問(wèn)題，土阻尼力與振幅有關(guān)也與應(yīng)變速率有關(guān)，采用滯回阻尼模型在概念上會(huì)引起矛盾，此時(shí)用黏性阻尼模型更為合適 [28-32]。因此本文考慮樁周土體施工擾動(dòng)，土體采用黏性阻尼模型，基于復(fù)剛度傳遞多圈層平面應(yīng)變模型，對(duì)任意激振力作用下徑向非均質(zhì)、縱向成層黏性阻尼土中變截面、變模量管樁縱向振動(dòng)特性進(jìn)行解析理論研究。

1 定解問(wèn)題數(shù)學(xué)模型的建立

1.1 計(jì)算簡(jiǎn)圖及基本假定

將管樁-土體耦合振動(dòng)系統(tǒng)沿縱向由地基土層界面及樁身缺陷截面分成m段（如果截面變化處不與土層分界線重合，每個(gè)截面變化處對(duì)應(yīng)一條虛擬的土層分界線），將樁長(zhǎng)為H管樁自樁身底部由下往上依次編號(hào)為1，2，…，i，…，m層段，各層段厚度分別為l1，l2，…，li，…，lm，各層段頂部埋深分別為h1，h2，…，hi，…，hm。第i層段管樁內(nèi)徑、外徑、樁段截面積、密度和彈性模量分別為ri0，ri1，APi，ρPi和EPi，樁底黏彈性支承剛度系數(shù)為δp，kp。在縱向第i層段中，樁芯土體剪切模量、黏性阻尼系數(shù)和密度分別為Gi0，ηi0，ρi0。

同時(shí)，將縱向第i層段樁周土體沿徑向劃分為內(nèi)部擾動(dòng)區(qū)域和外部區(qū)域，樁周土體內(nèi)部擾動(dòng)區(qū)域徑向厚度為bi，并將內(nèi)部擾動(dòng)區(qū)域沿徑向劃分n個(gè)圈層，第j圈層土體剪切模量、黏性阻尼系數(shù)和密度分別為Gij，ηij，ρij，第j-1個(gè)圈層與第j圈層的界面處半徑為rij。特別地，內(nèi)部區(qū)域和外部區(qū)域界面處的半徑為ri（n+1），外部區(qū)域則為徑向半無(wú)限均勻黏彈性介質(zhì)。管樁樁頂作用任意激振力p（t），第i層段樁芯土和樁周土對(duì)樁身產(chǎn)生的切應(yīng)力分別為fS0i和fS1i，樁土耦合振動(dòng)體系力學(xué)簡(jiǎn)化模型如圖1所示。

基本假定如下：

（1）管樁第i段樁身假定為均質(zhì)等截面彈性體，樁體底部為黏彈性支承，樁身相鄰層段之間滿足力平衡和位移連續(xù)條件。

（2）第i層段樁周土體內(nèi)部擾動(dòng)區(qū)域沿徑向所劃分的n個(gè)圈層和樁芯土體都為均質(zhì)、各向同性黏彈性體；外部區(qū)域?yàn)閺较虬霟o(wú)限均勻黏彈性介質(zhì)。

（3）管樁-土體耦合振動(dòng)系統(tǒng)滿足線彈性和小變形條件。

（4）在各層段中，樁周土、樁芯土與管樁內(nèi)外樁壁土界面上產(chǎn)生的剪應(yīng)力，分別通過(guò)各自樁土界面剪切復(fù)剛度傳遞給樁身，樁土之間完全接觸。

（5）各層段中樁周土剪切波速?gòu)耐獠繀^(qū)域至內(nèi)部擾動(dòng)區(qū)域最內(nèi)圈層呈現(xiàn)線性變化，即剪切模量呈現(xiàn)二次函數(shù)變化規(guī)律，樁周土體黏性阻尼系數(shù)與剪切模量相同。

1.2 定解問(wèn)題

3 算例分析

本文算例基于圖1所示黏彈性土體中管樁耦合振動(dòng)力學(xué)模型，采用前述推導(dǎo)求解所得基于黏性阻尼模型徑向非均質(zhì)、縱向成層樁周土中變模量和變截面管樁樁頂縱向振動(dòng)動(dòng)力阻抗和時(shí)域速度響應(yīng)函數(shù)解答。

已有研究結(jié)果表明，當(dāng)樁周土體徑向圈層數(shù)量n=20以上時(shí)[16]，可以忽略劃分圈層數(shù)對(duì)計(jì)算結(jié)果的影響。這樣，在后續(xù)分析中將管樁-土體耦合振動(dòng)系統(tǒng)沿縱向分成5個(gè)層段，各層段中的樁周土體內(nèi)部區(qū)域則統(tǒng)一沿徑向劃分為20個(gè)圈層。

3.1 解答合理性驗(yàn)證

為了驗(yàn)證本文推導(dǎo)所得徑向非均質(zhì)、縱向成層黏性阻尼土體中變模量和變截面管樁樁頂縱向振動(dòng)動(dòng)力阻抗解析解答合理性，將徑向非均質(zhì)、縱向成層樁周土退化成均質(zhì)土，變模量和變截面管樁退化成均勻樁，并與文獻(xiàn)[33]已有均質(zhì)土中均勻管樁縱向動(dòng)力阻抗解析解進(jìn)行對(duì)比驗(yàn)證。如圖2所示為本文推導(dǎo)所得管樁縱向振動(dòng)動(dòng)力阻抗退化解隨無(wú)量綱頻率θ變化曲線，及其與文獻(xiàn)[33]已有阻抗解的對(duì)比情況。從圖中不難看出，本文所推導(dǎo)的雙向非均質(zhì)土體中變模量和變截面管樁縱向振動(dòng)動(dòng)力阻抗退化解答曲線與文獻(xiàn)[33]中對(duì)應(yīng)結(jié)果吻合。

3.2 樁身缺陷對(duì)樁頂動(dòng)力響應(yīng)影響分析

夾泥、空洞和離析等缺陷均會(huì)引起樁身模量突變，為分析此類缺陷對(duì)樁頂動(dòng)力響應(yīng)的影響，參照文獻(xiàn)[21]的方法將樁身缺陷概化為一模量突變段，并定義λ為缺陷段樁身模量與正常段樁身模量之比。假設(shè)在埋深3.8 m處存在一長(zhǎng)度為1 m的模量突變段，圖3所示為樁身模量突變對(duì)管樁樁頂動(dòng)力響應(yīng)特性的影響。由圖可見(jiàn)，均質(zhì)管樁（λ=1.0時(shí)）樁頂速度導(dǎo)納曲線振幅逐漸衰減，而變模量管樁樁身存在反射界面，其樁頂速度導(dǎo)納曲線振幅存在疊加現(xiàn)象，即具有大、小峰值交替循環(huán)特征。相較于均勻樁而言，當(dāng)λ<1時(shí)，樁頂速度導(dǎo)納曲線的大、小峰幅值差隨λ減小而增大，共振頻率隨λ減小而減小，且均小于均質(zhì)樁共振頻率。（當(dāng)頻率接近樁土系統(tǒng)的固有頻率時(shí)產(chǎn)生共振，共振頻率則對(duì)應(yīng)速度導(dǎo)納曲線上振幅極大值處）。當(dāng)λ>1時(shí)，樁頂速度導(dǎo)納曲線上的大、小峰幅值差和共振頻率均隨λ增加而增大，且對(duì)應(yīng)共振頻率均大于均質(zhì)樁情況。樁頂反射波曲線中顯著呈現(xiàn)出模量突變界面處的反射信號(hào)。具體地，當(dāng)λ>1時(shí)，由于樁身模量在突變段上界面處增大，從而呈現(xiàn)出與樁尖處反相的反射信號(hào)，在突變段下界面處樁身彈性模量減小則呈現(xiàn)出與樁尖處同相的反射信號(hào)，即當(dāng)λ>1時(shí)，模量突變界面反射信號(hào)特征為先反相后同相。相反地，當(dāng)λ<1時(shí)，模量突變界面反射信號(hào)特征為先同相后反相。

圖4所示為管樁樁身模量突變段深度位置對(duì)樁頂動(dòng)力響應(yīng)特性的影響。由圖可見(jiàn)，管樁樁身模量突變段位置對(duì)樁頂速度導(dǎo)納曲線影響顯著。具體地，較均質(zhì)樁而言，管樁樁身模量突變段位置深度越大，其樁頂速度導(dǎo)納曲線中的大、小峰幅值差越小，其對(duì)速度導(dǎo)納曲線的峰值水平影響亦越小。由反射波曲線可以看出，管樁樁身模量突變段位置深度越小，樁身模量突變界面處反射波信號(hào)出現(xiàn)的時(shí)間越早，反射信號(hào)幅值水平也越高。

為分析樁身縮頸對(duì)樁頂動(dòng)力響應(yīng)特性的影響規(guī)律，假設(shè)在埋深3.8 m處存在一長(zhǎng)度為1 m的縮頸段，管樁樁身等壁厚縮頸的內(nèi)、外徑工況變化如表1所列，其中工況Case1為均勻截面樁工況。圖5所示為所列內(nèi)、外徑工況下管樁樁頂速度導(dǎo)納和反射波曲線變化情況。

綜合圖5和表1可知，相對(duì)于均勻截面樁而言，縮頸樁速度導(dǎo)納曲線呈現(xiàn)大、小峰值交替現(xiàn)象。在所列各工況中，僅外徑減小引起的速度導(dǎo)納曲線大、小峰幅值差最大，而外徑減小且內(nèi)徑增大工況次之，僅內(nèi)徑增大引起的大小峰值幅值差則最小。管樁樁頂反射波曲線顯著呈現(xiàn)出樁身縮頸界面處的反射信號(hào)。具體地，由于樁身截面在縮頸段上界面處減小，從而呈現(xiàn)同相反射信號(hào)，在縮頸段下界面處樁身截面增大則呈現(xiàn)反相反射信號(hào)。即樁身出現(xiàn)縮頸段時(shí)，縮頸段界面反射信號(hào)特征為先同相后反相。此外，較均勻截面樁而言，樁身縮頸段的存在，使得樁尖反射波信號(hào)幅值水平降低。

圖6所示為管樁樁身縮頸段深度位置對(duì)樁頂動(dòng)力響應(yīng)特性的影響情況。由圖可見(jiàn)，管樁樁身縮頸段深度位置對(duì)樁頂速度導(dǎo)納曲線影響顯著。具體地，較均勻截面樁而言，管樁樁身縮頸段位置深度越大，其樁頂速度導(dǎo)納曲線中的大、小峰幅值差越小，其對(duì)速度導(dǎo)納曲線的峰值水平影響亦越小。由反射波曲線可以看出，管樁樁身縮頸段位置深度越小，樁身縮頸段界面處反射波信號(hào)出現(xiàn)的時(shí)間越早，反射信號(hào)幅值水平也越高。

為分析樁身擴(kuò)頸對(duì)樁頂動(dòng)力響應(yīng)特性的影響規(guī)律，假設(shè)在埋深3.8 m處存在一長(zhǎng)度為1 m的擴(kuò)頸段，管樁樁身等壁厚擴(kuò)頸的內(nèi)、外徑工況變化如表2所列，其中工況Case5為均勻截面樁工況。圖7所示為所列內(nèi)、外徑工況下管樁樁頂速度導(dǎo)納和反射波曲線變化情況。綜合圖7和表2可知，相對(duì)于均勻截面管樁而言，擴(kuò)頸管樁樁頂速度導(dǎo)納曲線呈現(xiàn)大、小峰值交替現(xiàn)象。在所列各工況中，僅內(nèi)徑減小引起的速度導(dǎo)納曲線大、小峰幅值差最大，而僅外徑增大工況次之，外徑增大且內(nèi)徑減小工況引起的大小峰值幅值差則最小。管樁樁頂反射波曲線顯著呈現(xiàn)出樁身擴(kuò)頸界面處的反射信號(hào)。具體地，由于樁身截面在擴(kuò)頸段上界面處增大，從而呈現(xiàn)出與樁尖處反相的反射信號(hào)，在擴(kuò)頸段下界面處樁身截面減小則呈現(xiàn)出與樁尖處同相的反射信號(hào)。即樁身出現(xiàn)擴(kuò)頸段時(shí)，擴(kuò)徑段界面反射信號(hào)特征為先反相后同相。

圖8所示為管樁樁身擴(kuò)頸段深度位置對(duì)樁頂動(dòng)力響應(yīng)特性的影響情況。由圖可見(jiàn)，管樁樁身擴(kuò)頸段深度位置對(duì)樁頂速度導(dǎo)納曲線影響顯著。具體地，較均勻截面樁而言，管樁樁身擴(kuò)頸段位置深度越大，其樁頂速度導(dǎo)納曲線中的大、小峰幅值差越小，其對(duì)速度導(dǎo)納曲線的峰值水平影響亦越小。由反射波曲線可以看出，管樁樁身擴(kuò)頸段位置深度越小，樁身擴(kuò)頸段界面處反射波信號(hào)出現(xiàn)的時(shí)間越早，反射信號(hào)幅值水平也越高。

4 結(jié) 論

本文考慮樁周土體徑向施工擾動(dòng)效應(yīng)，基于黏性阻尼土體平面應(yīng)變模型和復(fù)剛度傳遞法，推導(dǎo)得出了任意激振力作用下徑向非均質(zhì)、縱向成層的黏性阻尼土中管樁縱向振動(dòng)動(dòng)力阻抗和速度導(dǎo)納頻域解析解答，以及速度時(shí)域響應(yīng)半解析解答，并將本文解與均質(zhì)土中管樁已有解答進(jìn)行了退化對(duì)比驗(yàn)證，在此基礎(chǔ)上，分別探討了樁身變模量、縮頸和擴(kuò)徑對(duì)樁頂縱向動(dòng)力響應(yīng)特性的影響規(guī)律，計(jì)算分析結(jié)果表明：

（1）相對(duì)于均質(zhì)樁而言，存在模量突變段管樁的樁頂速度導(dǎo)納曲線呈現(xiàn)大、小峰交替現(xiàn)象。當(dāng)變模量系數(shù)λ<1時(shí)，樁頂速度導(dǎo)納曲線的大、小峰幅值差隨λ減小而增大，共振頻率隨λ減小而減小，且均小于均質(zhì)樁共振頻率。相對(duì)于樁尖處反射信號(hào)而言，模量突變界面反射信號(hào)特征為先反相后同相。當(dāng)變模量系數(shù)λ>1時(shí)，樁頂速度導(dǎo)納曲線上的大、小峰幅值差和共振頻率均隨λ增加而增大，且對(duì)應(yīng)共振頻率均大于均質(zhì)樁情況。而模量突變界面反射信號(hào)特征則為先同相后反相。

（2）相對(duì)于均勻截面樁而言，縮頸管樁樁頂速度導(dǎo)納曲線呈現(xiàn)大、小峰值交替現(xiàn)象。在所列各工況中，僅外徑減小引起的速度導(dǎo)納曲線的大、小峰幅值差最大，而僅內(nèi)徑增大引起的大小峰值幅值差則最小。管樁樁頂反射波曲線顯著呈現(xiàn)出樁身縮頸界面處的反射波信號(hào)。相對(duì)于樁尖處反射信號(hào)而言，縮頸段界面反射信號(hào)特征為先同相后反相。此外，較均勻截面樁而言，樁身縮頸段的存在，使得樁尖反射波信號(hào)幅值水平降低。

（3）相對(duì)于均勻截面管樁而言，擴(kuò)頸管樁樁頂速度導(dǎo)納曲線呈現(xiàn)大、小峰值交替現(xiàn)象。在所列各工況中，僅內(nèi)徑減小引起的樁頂速度導(dǎo)納曲線大、小峰幅值差最大，而外徑增大且內(nèi)徑減小工況引起的大小峰值幅值差則最小。管樁樁頂反射波曲線顯著呈現(xiàn)出樁身擴(kuò)頸界面處的反射波信號(hào)。相對(duì)于樁尖處反射信號(hào)而言，擴(kuò)頸段界面反射信號(hào)特征為先反相后同相。

（4）樁身模量突變段、縮頸段和擴(kuò)徑段位置深度越大，其對(duì)樁頂速度導(dǎo)納曲線的峰值影響越小，且對(duì)應(yīng)的大、小峰幅值差亦越小。對(duì)于樁頂反射波曲線而言，模量突變段、縮頸段和擴(kuò)徑段位置深度越小，各自對(duì)應(yīng)的樁身缺陷段界面處反射波信號(hào)出現(xiàn)時(shí)間越早，反射信號(hào)幅值水平也越高。

參考文獻(xiàn)：

[1] 崔春義， 張石平， 楊 剛， 等.考慮樁底土層波動(dòng)效應(yīng)的飽和黏彈性半空間中摩擦樁豎向振動(dòng)[J]. 巖土工程學(xué)報(bào)，2015，37（5）：878—892.

CUI Chun-yi， ZHANG Shi-ping， YANG Gang， et al.Vertical vibration of floating piles in saturated viscoelastic half-spaceconsidering wave effect of subsoil under pile bottom [J]. Chinese Journal of Geotechnical Engineering， 2015，37（5）：878—892.

[2] Cui Chun-yi， Zhang Shi-ping， Yang Gang， et al. Vertical vibration of a floating pile in a saturated viscoelastic soil layer overlaying bedrock[J]. Journal of Central South University，2016，23（1）： 220—232.

[3] Han Y， Vaziri H. Dynamic response of pile groups under lateral loading[J]. Soil Dynamics & Earthquake Engineering， 1992， 11（2）：87—99.

[4] Wu W B， Jiang G S， Dou B， et al. Vertical dynamic impedance of tapered pile considering compacting effect[J]. Mathematical Problems in Engineering， 2013，（2）：1—9.

[5] Novak M，Sheta M. Approximate approach to contact problems of piles [C]. Proceedings of the Geotechnical Engineering Division， American Society of Civil Engineering National Convention， Florida， 1980： 53—79.

[6] Novak M， Han Y C. Impedances of soil Layer with boundary zone[J]. Journal of Geotechnical Engineering， 1990， 116（6）：1008—1014.

[7] Veletsos A S， Dotson K W. Impedances of soil layer with disturbed boundary zone[J]. Journal of Geotechnical Engineering， 1986， 112（3）：363—368.

[8] Veletsos A S， Dotson K W. Vertical and torsional vibration of foundations in inhomogeneous media[J]. Journal of Geotechnical Engineering， 1988， 114（9）：1002—1021.

[9] Nogami T， Konagai K. Dynamic response of vertically loaded nonlinear pile foundations[J]. Journal of Geotechnical Engineering， 1987， 113（2）：147—160.

[10]Nogami T， Konagai K. Time domain flexural response of dynamically loaded single piles[J]. Journal Engineering Mechanics， ASCE， 1988， 114（9）： 1512—1525.

[11]Dotson K W， Veletsos A S. Vertical and torsional impedances for radially inhomogeneous viscoelastic soil layers[J]. Soil Dynamics & Earthquake Engineering， 1990， 9（3）：110—119.

[12]Vaziri H， Han Y. Impedance functions of piles in inhomogeneous media[J]. Journal of Geotechnical Engineering， 1993， 119（9）：1414—1430.

[13]Han Y C， Sabin G C W. Impedances for radially inhomogeneous viscoelastic soil media[J]. Journal of Engineering Mechanics， 1995， 121（9）：939—947.

[14]Han Y C. Dynamic vertical response of piles in nonlinear soil[J]. Journal of Geotechnical & Geoenvironmental Engineering， 1997， 123（123）：710—716.

[15]EI Naggar M H， Novak M. Nonlinear axial interaction in pile dynamics[J]. Journal of Geotechnical Engineering， 1994， 120（4）： 678—696.

[16]EI Naggar M H. Vertical and torsional soil reactions for radially inhomogeneous soil layer[J]. Structural Engineering & Mechanics， 2000， 10（4）：299—312.

[17]王奎華， 楊冬英， 張智卿，等. 基于復(fù)剛度傳遞多圈層平面應(yīng)變模型的樁動(dòng)力響應(yīng)研究[J]. 巖石力學(xué)與工程學(xué)報(bào)， 2008， 27（4）： 825—831.

WANG Kui-hua， YANG Dong-ying， ZHANG Zhi-qing， et al. Study on dynamic response of pile based on complex stiffness transfer model of radial multizone plane strain[J]. Chinese Journal of Rock Mechanics & Engineering， 2008， 27（4）：825—831.

[18]楊冬英， 王奎華. 任意圈層徑向非均質(zhì)土中樁的縱向振動(dòng)特性[J]. 力學(xué)學(xué)報(bào)， 2009， 41（2）：243—252.

YANG Dong-ying， WANG Kui-hua. Vertical vibration of pile in radially inhomogeneous soil layers[J]. Chinese Journal of Theoretical & Applied Mechanics， 2009， 41（2）：243—252.

[19]Wang K H， Yang D Y， Zhang Z Q， et al.A new approach for vertical impedance in radially inhomogeneous soil layer[J]. International Journal for Numerical & Analytical Methods in Geomechanics， 2012， 36（6）：697—707.

[20]Yang D Y， Wang K H， Zhang Z Q， et al.Vertical dynamic response of pile in a radially heterogeneous soil layer[J]. International Journal for Numerical & Analytical Methods in Geomechanics， 2009， 33（8）：1039—1054.

[21]王奎華， 謝康和， 曾國(guó)熙. 變截面阻抗樁受迫振動(dòng)問(wèn)題解析解及應(yīng)用[J]. 土木工程學(xué)報(bào)， 1998， 31（6）， 56—67.

WANG Kui-hua， XIE Kang-he， ZENG Guo-xi. An analytical solution to forced vibration of foundation pile with variable section impedance and its application[J]. China Civil Engineering Journal， 1998， 31（6）：56—67.

[22]王宏志， 陳云敏， 陳仁朋. 分層土中變截面樁振動(dòng)問(wèn)題半解析解的研究[J]. 振動(dòng)與沖擊， 2001， 20（1）：55—58.

WANG Hong-zhi， CHEN Yun-min， CHEN Ren-peng. Semi-analytical solution for vibration of pile with variable section in layered soil[J]. Journal of Vibration & Shock， 2001， 20（1）： 55—58.

[23]劉東甲. 不均勻土中多缺陷樁的軸向動(dòng)力響應(yīng)[J]. 巖土工程學(xué)報(bào)， 2000， 22（4）：391—395.

LIU Dong-jia. Dynamic axial response of multi-defective piles in nonhomogeneous soil[J]. Chinese Jounal of Geotechnical Engineering， 2000， 22（4）： 391—395.

[24]馮世進(jìn)， 柯 瀚， 陳云敏，等. 成層土中粘彈性變截面樁縱向振動(dòng)分析及應(yīng)用[J]. 巖石力學(xué)與工程學(xué)報(bào)， 2004， 23（16）：2798—2803.

FENG Shi-jin， KE Han， CHEN Yun-min， et al. Study on vertical vibration of viscoelastic pile with variable sections in layered soils[J]. Chinese Journal of Rock Mechanics & Engineering， 2004， 23（16）：2798—2803.

[25]丁選明， 劉漢龍. 大直徑管樁在瞬態(tài)集中荷載作用下的振動(dòng)響應(yīng)時(shí)域解析解[J]. 巖土工程學(xué)報(bào)， 2013， 35（6）：1010—1017.

DING Xuan-ming， LIU Han-long. Time-domain analytical solution of the vibration response of a large-diameter pipe pile subjected to transient concentrated load[J]. Chinese Journal of Geotechnical Engineering， 2013， 35（6）：1010—1017.

[26]鄭長(zhǎng)杰， 丁選明， 劉漢龍，等. 考慮土體三維波動(dòng)效應(yīng)的現(xiàn)澆大直徑管樁縱向振動(dòng)響應(yīng)解析解[J]. 巖土工程學(xué)報(bào)， 2013， 35（12）：2247—2254.

ZHENG Chang-jie， DING Xuan-ming， LIU Han-long， et al. Analytical solution to vertical vibration of cast-in-place concrete large-diameter pipe piles by considering 3D wave effect of soils[J]. Chinese Journal of Geotechnical Engineering， 2013， 35（12）： 2247—2254.

[27]Li Z Y， Wang K H， Wu W B， et al. Vertical vibration of a large-diameter pipe pile considering the radial inhomogeneity of soil caused by the construction disturbance effect[J]. Computers & Geotechnics， 2017： 90—102.

[28]胡昌斌， 王奎華， 謝康和. 樁與黏性阻尼土耦合縱向振動(dòng)時(shí)樁頂時(shí)域響應(yīng)研究[J]. 振動(dòng)工程學(xué)報(bào)， 2004， 17（1）：72—77.

HU Chang-bin， WANG Kui-hua， XIE Kang-he. Time domain axial response of dynamically loaded pile in viscous damping soil layer[J]. Journal of Vibration Engineering， 2004， 17（1）：72—77.

[29]Nogami T， Novak M. Soil-pile interaction in vertical vibration[J]. Earthquake Engineering & Structural Dynamics， 1976， 4（3）：277—293.

[30]Militano G， Rajapakse R K N D. Dynamic response of a pile in a multi-layered soil transient torsional and axial loading[J]. Géotechnique， 1999， 49（1）：91—109.

[31]胡海巖. 結(jié)構(gòu)阻尼模型及系統(tǒng)時(shí)域動(dòng)響應(yīng)[J]. 應(yīng)用力學(xué)學(xué)報(bào)， 1993，（1）：8—16.

HU Hai-yan. Structural damping model and system dynamic response at time domain[J]. Chinese Journal of Applied Mechanics， 1993，（1）：8—16.

[32]廖振鵬. 工程波動(dòng)理論導(dǎo)論[M]. 北京：科學(xué)出版社， 2002.

LIAO Zhen-peng.Introduction to Engineering Wave Theory[M]. Beijing： Science Press， 2002.

[33]丁選明， 陳 磊. 粘彈性地基中管樁的縱向振動(dòng)特性研究[J]. 長(zhǎng)江科學(xué)院院報(bào)， 2009， 26（3）：32—35.

DING Xuan-ming， CHEN Lei. Study on vertical vibration characteristics of pipe piles in viscoelastic soil[J]. Journal of Yangtze River Scientific Research Institute， 2009， 26（3）：32—35.

Abstract： A simplified mechanical model of vertical vibration of the pipe pile embedded in radially inhomogeneous and longitudinally layered viscoelastic soil considering both the variations of cross-section and modulus along the pile shaft is proposed， by employing the viscous damping model and the plane stain model of annular complex stiffness transfer. Firstly， the complex stiffness at the interface between soil and pile is derived based on Laplace transform and complex stiffness transfer method. Secondly， the analytical solutions for dynamic impedance， velocity admittance and reflected wave signal at the pile head are obtained by using the pile-soil compatibility condition of pipe pile and radial inhomogeneous surrounding soil. Furthermore， the obtained analytical solution for dynamic impedance at the pile head is also reduced to verify its validity with existing solution. Finally，an extensive parametric analysis is conducted to investigated the effects of the defective pile segments with variable modulus， necking and expanding， and the corresponding defect locations of pile shaft on the dynamic response at the pile head of the pipe pile embedded in radially inhomogeneous and longitudinal layered viscoelastic soil， which can provide reference and guide for specific engineering practice.

Key words： pipe pile； dynamic response； complex stiffness transfer model； radial heterogeneity； viscous damping

作者簡(jiǎn)介：

崔春義（1978—），男，副教授，博士后。E-mail： cuichunyi@dlmu.edu.cn

通訊作者：王科盛（1978—），男，博士，副教授。E-mail： keshengwang@uestc.edu.cn