(華南理工大學(xué)土木與交通學(xué)院 廣東 廣州 510641)

一、引言

目前，雖然反射波法能夠較準(zhǔn)確地判斷缺陷的位置與基本類型，但是由于樁側(cè)土的存在，對(duì)基樁缺陷程度仍然只能進(jìn)行定性分析[1]，即還不能夠確定缺陷的大小，因此在定量分析方面的有限元仿真研究，仍然存在一些問(wèn)題。問(wèn)題體現(xiàn)在例如樁側(cè)土阻尼如何確定，樁側(cè)土與樁身的相互作用定義等。本文探討了有限元模擬的實(shí)現(xiàn)方式，并對(duì)此給出了適當(dāng)?shù)慕ㄗh。

二、三種有限元模擬方式定義

為了嘗試定量分析，探究缺陷處相對(duì)入射波的波峰比值與缺陷程度的關(guān)系，這里將樁土相互作用簡(jiǎn)化為一維波動(dòng)問(wèn)題，樁視為一維桿件。我們定義了三種數(shù)值模擬方式，其區(qū)別在于定義樁身與樁周土的互相作用。

第一種方法將樁身和樁周土之間定義為接觸關(guān)系，樁身為主面，土為從面，同時(shí)定義切向行為是摩擦作用和法向行為是硬接觸；

第二種方法直接將樁身和樁周土之間采用綁定約束，互相之間無(wú)相對(duì)位移；方法一二均將樁身的材料設(shè)為均勻，此處其阻尼不考慮。

阻尼的因素主要考慮為樁周土的材料阻尼。

第三種方法將土對(duì)樁的作用簡(jiǎn)化為粘滯阻尼器;

通過(guò)將土壤分層，各層土對(duì)樁身作用看做一對(duì)一對(duì)的彈簧并聯(lián)阻尼。

三、有限元模擬方式存在的問(wèn)題

對(duì)于方法一二，在abaqus中，土體材料阻尼只能通過(guò)Rayleigh阻尼定義，僅用質(zhì)量陣和剛度陣系數(shù)α和β來(lái)表示阻尼的大小。同時(shí)有研究指出，α對(duì)樁身反射波衰減幾乎無(wú)影響，因此基本只有β在控制土體材料阻尼[2]，β的確定暫未有合理科學(xué)的公式來(lái)反映其與土體剪切波速、土體密度、樁身尺寸等關(guān)鍵參數(shù)的聯(lián)系。因此此處我們重點(diǎn)對(duì)方法三進(jìn)行了研究。對(duì)于方法三，相對(duì)方法一二，雖然模型進(jìn)行了簡(jiǎn)化，但是對(duì)于樁側(cè)土阻尼參數(shù)的確定，可以按下列公式[3]：

ci=ρicsi
cd=ρdc
kd=2DG/(1-υ)

(1)

式中：ci表示第i層土對(duì)樁側(cè)面的阻尼系數(shù)，Ns/m3；

cd表示樁底土對(duì)樁端的阻尼系數(shù)，Ns/m3；

p均表示土的密度；

csi和c表示第i層土的剪切波波速和樁底土的縱波波速，m/s;

kd表示樁底土的剛度系數(shù)，N/m；

D表示樁直徑，G表示樁底土剪切模量，v表示樁底土泊松比；

四、仿真計(jì)算

(一)算例概況與計(jì)算參數(shù)

所采用的鋼筋混凝土預(yù)制圓樁，樁長(zhǎng) 20m,泥面以上樁段長(zhǎng)5m,缺陷考慮為縮頸，最窄處等效半徑為0.3m,缺陷段距樁頂9m,缺陷段長(zhǎng)1m。計(jì)算參數(shù)如下：

表1 土體參數(shù)

(二)有限元計(jì)算結(jié)果

加載方式均為在整個(gè)面上的分布加載，縱波波速取1500m/s，剪切波速為100m/s.注：下圖橫坐標(biāo)為傳播距離L，與時(shí)間的換算關(guān)系為

L=T*c/2

(2)

L為長(zhǎng)度，單位m;T為時(shí)間，單位s;c為彈性波樁身內(nèi)波速,m/s。

某點(diǎn)離樁頂距離

D=L(X)-L0

(3)

L(X)為任意一點(diǎn)橫坐標(biāo)，L0為首波峰傳播距離，L0此處約為1.60m。

觀察發(fā)現(xiàn)兩類解法首波峰和樁底反射幅值基本吻合，但有限元解缺陷處速度相對(duì)解析解[4]衰減過(guò)大，經(jīng)過(guò)試算確定，樁側(cè)阻尼降低到原來(lái)一半時(shí)，有限元解與解析解較為接近，此時(shí)對(duì)比圖如下所示：

圖3 解析解對(duì)比1/2原樁側(cè)阻尼下有限元解法

五、結(jié)論

綜上所述，對(duì)于定量分析確定缺陷程度，反射波法有限元模擬方法存在問(wèn)題主要集中在樁側(cè)土阻尼無(wú)法確定；采用樁土耦合模型(樁土設(shè)置接觸單元、樁土綁定)的問(wèn)題在于Rayleigh阻尼的剛度阻尼系數(shù)無(wú)法準(zhǔn)確反映土體參數(shù)對(duì)阻尼的影響，缺乏科學(xué)合理的公式確定阻尼的大小。采用簡(jiǎn)化模型，即土彈簧阻尼模型，有經(jīng)驗(yàn)公式可供指導(dǎo)。有研究表明,理論阻尼系數(shù)應(yīng)用于實(shí)際工程時(shí)需要進(jìn)行修正,修正系數(shù)是剪切波速的函數(shù)[5]，可通過(guò)場(chǎng)地實(shí)驗(yàn)測(cè)定剪切波速與阻尼系數(shù)的若干組實(shí)際數(shù)據(jù)，擬合出兩者的函數(shù)關(guān)系，從而得出修正系數(shù)。

[1]黃理興．動(dòng)測(cè)樁身完整性的新理念[J][巖石力學(xué)與工程學(xué)報(bào)]2002·21(3)

[2]基樁缺陷反射波法檢測(cè)數(shù)值模擬分析_李晴[D]成都，西南交通大學(xué)2012

[3]樁基工程手冊(cè)第三版 [S](人民交通出版社)

[4]單排高樁基礎(chǔ)豎向振動(dòng)理論及應(yīng)用研究[D]杭州，呂述暉，浙江大學(xué)

[5]時(shí)域內(nèi)低應(yīng)變基樁缺陷定量分析方法研究-[D]天津，智勝英，天津大學(xué)