柳 偉，張斌偉，柳德龍

（1.蘭州新區(qū)城市發(fā)展投資集團(tuán)有限公司，蘭州730087；2.隴東學(xué)院 土木工程學(xué)院，甘肅 慶陽745000；3.甘肅省高?！包S土工程性質(zhì)及工程應(yīng)用”省級重點(diǎn)實驗室，甘肅 慶陽745000；4.中國鐵路烏魯木齊局集團(tuán)有限公司，烏魯木齊830011）

樁基作為眾多工程領(lǐng)域中的一種深基礎(chǔ)結(jié)構(gòu)形式，以其簡單、適用而廣泛應(yīng)用于鐵路、公路、機(jī)場、海洋平臺、高層建筑地基基礎(chǔ)等，其振動特性一直在科學(xué)以及工程應(yīng)用中備受關(guān)注[1-3]。

土與結(jié)構(gòu)相互作用下這一問題的復(fù)雜性，歷來都是巖土工程、結(jié)構(gòu)工程、固體力學(xué)接觸問題中的棘手難題。蒲育、滕兆春[4]基于Hamilton 原理和線彈性理論，應(yīng)用微分求積法研究彈性地基梁自由振動的無量綱頻率特性，但沒有考慮實際工程中土體阻尼的影響。張阿舟[5]在考慮樁周土剛度及阻尼的情況下采用分離變量法求解了全埋置一維樁的自振頻率及振型，若運(yùn)用該方法在求解二維埋置結(jié)構(gòu)自振特性時，得到的頻率方程以超越方程形式表達(dá)，計算會更加復(fù)雜。彭麗等[6-7]采用復(fù)模態(tài)方法分析了黏彈性三參數(shù)地基上梁的橫向振動特性，并用微分求積方法加以驗證，但沒有分析衰減系數(shù)對自振特性的影響。陳興沖[8]采用瑞利法推導(dǎo)計算淺平基橋墩、樁基橋墩及沉井基礎(chǔ)橋墩基頻的近似公式，但該方法不易推廣至高階自振頻率的計算。楊驍?shù)萚9]將樁等效為Reyleigh 梁，利用精確有限元法求解成層液化土單樁-土-結(jié)構(gòu)系統(tǒng)的固有頻率，但運(yùn)用該方法求解Timoshenko梁理論模型時，它的形函數(shù)不再精確，必須以增加結(jié)構(gòu)單元數(shù)量為代價提高計算精度，不易推廣至現(xiàn)實的工程中去，而波動理論對以上問題的解決提供了一條有效的途徑。

回傳射線矩陣自1998年首次引入結(jié)構(gòu)計算以來，已被成功應(yīng)用于框架結(jié)構(gòu)及各向同性層狀介質(zhì)的瞬態(tài)響應(yīng)及振動分析[10-13]，其物理意義明確、列式統(tǒng)一、易于編程，且具有高精度、低耗時、結(jié)果可讀性強(qiáng)等優(yōu)點(diǎn)，在結(jié)構(gòu)動力學(xué)分析中具有獨(dú)特的優(yōu)勢。截止目前，對黏彈性地基中結(jié)構(gòu)自振特性的相關(guān)研究工作還很少。

本文將回傳射線矩陣法推廣至樁土系統(tǒng)的振動分析中，運(yùn)用回傳射線矩陣法及求根法，利用MATLAB語言編程，通過具體算例，分析了外露長度、埋置深度、樁端約束情況對埋置結(jié)構(gòu)自振特性的影響，所得結(jié)論不僅對結(jié)構(gòu)的質(zhì)量檢測具有指導(dǎo)意義及工程應(yīng)用價值，而且可以為結(jié)構(gòu)設(shè)計和施工計算提供理論基礎(chǔ)。

1 回傳射線法及求根法的基本原理

本文基于Timoshenko 梁理論的Winkler 地基模型，將樁基劃分為2 個單元3 個節(jié)點(diǎn)，節(jié)點(diǎn)的編號如圖1（a），建立整體坐標(biāo)系(X,Y)，引入2 個對偶局部坐標(biāo)系(x,y)JK和(x,y)KJ，如圖1（b）所示。

在局部坐標(biāo)系下，12、23單元的波動控制方程為

圖1 與黏彈性地基相互作用單樁的計算模型及局部坐標(biāo)系

式中：u(x,t)為軸向位移；vb(x,t)為彎矩引起的撓度，vs(x,t)為剪力引起的撓度；l、l12、l23分別為單樁的總長、外露長度、埋置長度；E、G、ρ、h、A、k'、Iz分別為單樁的彈性模量、剪切模量、密度、橫截面高度、橫截面面積、截面剪切系數(shù)、橫截面慣性矩；kuj、kvj、βuj、βvj、γj分別為土體縱向彈簧系數(shù)、土體橫向彈簧系數(shù)、土體縱向阻尼系數(shù)、土體橫向阻尼系數(shù)和樁身截面轉(zhuǎn)動時的土體摩擦阻尼系數(shù)；j=1 時，ku1=βu1=kv1=βu1=γ1=0；j=2 時，ku2=ku，kv2=kv，βu2=βu，βv2=βv，γ2=γ。

將式（1）式（2）直接進(jìn)行Fourier變換并整理得

求解式（3）－式（4）得位移在頻域中的表達(dá)式為

式中：a1(ω)、a2(ω)、a3(ω)為待定的入射波波幅；d1(ω)、d2(ω)、d3(ω)為待定的出射波波幅；k1、k2、k3為波數(shù)，滿足

式中：j=1 時，k11(ω)、k21(ω)、k31(ω)對應(yīng)12 單元的波數(shù)；當(dāng)j=2 時，k12(ω)、k22(ω)、k32(ω)對應(yīng)23 單元的波數(shù)；

對應(yīng)于波數(shù)k2j,3j與的比值為

軸力、彎矩、剪力和轉(zhuǎn)角在頻域中的表達(dá)式為

如圖2所示的受力分析圖，對節(jié)點(diǎn)1建立力平衡和位移協(xié)調(diào)條件

根據(jù)圖3所示的受力分析圖，對節(jié)點(diǎn)2建立力平衡和位移協(xié)調(diào)條件。

圖2 節(jié)點(diǎn)1受力分析圖

圖3 節(jié)點(diǎn)2受力分析圖

根據(jù)圖4所示的受力分析圖，對節(jié)點(diǎn)3建立力平衡和位移協(xié)調(diào)條件

圖4 節(jié)點(diǎn)3受力分析圖

將式（10）－式（13）代入式（14）－式（16），并將局部坐標(biāo)系下的矩陣形式組集成整體坐標(biāo)系下的矩陣形式為

式中：d、a為總體出射波和入射波波幅向量；S、s為整體散射矩陣和整體源矢量。

從局部坐標(biāo)系的角度看，對于任一個單元JK，其中一端的入射波對另一端而言就是出射波，因此，入射波的波幅向量和出射波的波幅向量滿足以下相位關(guān)系

式中：aJK(ω)和dJK(ω)為同一個單元兩個節(jié)點(diǎn)處的入射波波幅和出射波波幅；PJK(lJK,ω)稱為傳播矩陣表示單元JK的長度。

將所有桿件單元的入射波波幅向量aJK(ω)和出射波波幅向量(ω)組集到總體入射波波幅向量a和總體出射波波幅向量d～中，寫成緊湊形式為

由于～與d中各元素相同，而排列順序有所變化，因此，引入置換矩陣U，以調(diào)節(jié)～中各元素在總體坐標(biāo)系中的相對位置，有

將式（20）代入式（17），得

式中：R=SPU為回傳射線矩陣；I為單位矩陣。

黏彈性地基中樁基自由振動時，其波源矩陣s=0，即

式中：d(ω) 有非零解的條件為系數(shù)行列式|I-R(ω) |必須為零，由于結(jié)構(gòu)在黏彈性地基中自由振動時需考慮土體阻尼對振動的衰減作用，令圓頻率ωn=+iδn。其中，實部為所求的自振頻率，虛部δn為對應(yīng)的衰減系數(shù)。

由于R(ω)中的各項是關(guān)于自振頻率ωˉn及衰減系數(shù)δn的二維復(fù)數(shù)超越方程。在數(shù)學(xué)上此類問題只有數(shù)值解而沒有解析解?？紤]到以上因素，根據(jù)回傳射線矩陣法的列式特點(diǎn)提出將二分法和黃金分割法(簡稱求根法)結(jié)合起來進(jìn)行迭代求解復(fù)雜復(fù)數(shù)方程，即：分別對ωn的實部和虛部δn進(jìn)行循環(huán)，當(dāng)[I-R(ω) ]的模小于預(yù)先給定的誤差時，取出對應(yīng)的ωn，則ωn的實部ωˉn為所求的自振頻率，虛部δn為對應(yīng)的衰減系數(shù)。

求得結(jié)構(gòu)的自振頻率及衰減系數(shù)后即可求得自由振動時的振型曲線，由線性代數(shù)的知識可知，矩陣與其行列式的值有如下關(guān)系

式中：adj[I-R]為矩陣[I-R]的伴隨矩陣；det[I-R]為矩陣[I-R]的模。

當(dāng)det[I-R]中ω的實部取自振頻率ωˉn，虛部δn取衰減系數(shù)時，式（23）為

設(shè)N是出射波波幅向量dk的維數(shù)，m為1≤m≤N的任意正整數(shù)，矩陣adj[I-R]的第m列記為dk，則

當(dāng)結(jié)構(gòu)振動頻率等于其自振頻率時，adj[I-R(ωk) ]的每一列均可看作是式（25）的解。即adj[I-R(ωk) ]的每一個非零列都可以做為黏彈性地基中單樁在自由振動時非零出射波的波幅向量，求得dk以后，代入式(22)可求得ak，將ak、dk代入式(5)－式(7)中，可求解結(jié)構(gòu)任意點(diǎn)處的位移，將各節(jié)點(diǎn)位移歸一化處理后即可得到黏彈性地基中單樁的振型曲線。

2 數(shù)值結(jié)果與分析

與黏彈性地基相互作用的單樁計算模型如圖1（a）所示，樁土相互作用用并聯(lián)的彈簧和黏壺來模擬，樁底為黏彈性支撐，樁徑D=1.25 m，泊松比μ=0.2337 ，彈 性 模 量E=4.321×1010Pa ，密 度ρ=267 6 kg m3，截面剪切系數(shù)k′=0.82；土體計算參數(shù)參考文獻(xiàn)[14]中的數(shù)值，如表1所列。

表1 黏彈性地基中土體的計算參數(shù)

2.1 埋置長度、外露長度對單樁自振特性的影響

以黏彈性地基中樁底彈性支撐、樁頂自由的單樁為例，根據(jù)回傳射線矩陣法及求根法，編寫計算任意階數(shù)的自振頻率和振動模態(tài)的MATLAB程序，如表2所列。

表2 3種工況下單樁的埋置長度、外露長度及樁長

圖5為埋置長度及外露長度對單樁自振頻率的影響。

圖5 埋置長度及外露長度對單樁自振頻率的影響

由圖可知，1 階自振頻率趨于重合，隨著黏彈性地基中單樁自振階數(shù)的增大，其自振頻率逐漸增大，尤其第4 階之后，這種增幅很明顯；埋置深度相同時，隨著外露長度越長，結(jié)構(gòu)的各階自振頻率越??；外露長度相同時，隨著埋置深度越深，結(jié)構(gòu)的各階自振頻率越小。

圖6 埋置長度及外露長度對單樁衰減系數(shù)的影響

圖6為埋置長度及外露長度對單樁衰減系數(shù)的影響。由圖可知，隨著黏彈性地基中單樁自振階數(shù)的增大，其衰減系數(shù)逐漸增大；埋置深度相同時，隨著外露長度越長，結(jié)構(gòu)的各階衰減系數(shù)越?。煌饴堕L度相同時，隨著埋置深度越深，結(jié)構(gòu)的各階衰減系數(shù)越大。

圖7為埋置長度及外露長度對單樁模態(tài)的影響。

由圖7可知，單樁基礎(chǔ)發(fā)生自由振動時，埋置長度和外露長度對結(jié)構(gòu)模態(tài)并無明顯影響；3 種工況下，隨著樁長的增大，結(jié)構(gòu)的1 階振型曲線趨于重合，隨著黏彈性地基中單樁自振階數(shù)的增大，其橫向位移增大，表明樁長對振型中橫向位移的影響較大；埋置部分的橫向位移小于外露部分的橫向位移，表明黏彈性地基對結(jié)構(gòu)的振動衰減作用明顯。

2.2 樁端約束情況對單樁自振特性的影響

圖8為樁端約束情況對單樁自振頻率的影響。

由圖可知，隨著黏彈性地基中單樁自振階數(shù)的增大，其自振頻率逐漸增大；樁頂固定工況下單樁各階自振頻率最大，樁頂鉸接工況下單樁各階自振頻率最小，樁頂自由工況下單樁各階自振頻率略小于

圖7 埋置長度及外露長度對單樁模態(tài)的影響

圖8 樁端約束情況對單樁自振頻率的影響

其在樁頂固定工況下的各階自振頻率，又大于其在樁頂鉸接工況下的各階自振頻率，這一結(jié)果與文獻(xiàn)［7］的結(jié)論一致。

圖9 樁端約束情況對單樁衰減系數(shù)的影響

圖9為樁端約束情況對單樁衰減系數(shù)的影響。由圖可知，隨著黏彈性地基中單樁自振階數(shù)的增大，其衰減系數(shù)逐漸增大；樁頂自由工況下單樁各階衰減系數(shù)最大，樁頂固定工況下單樁各階衰減系數(shù)最小，樁頂鉸接工況下單樁各階衰減系數(shù)介于其在樁頂自由和樁頂固定工況下的各階衰減系數(shù)之間。

圖10 樁端約束情況對單樁模態(tài)的影響

圖10為樁端約束情況對黏彈性地基中單樁基礎(chǔ)模態(tài)的影響。由圖可知，樁端相對樁底，其橫向位移和軸向位移均增大，表明不同的樁端約束情況對模態(tài)的影響較大；樁頂自由工況下單樁各階振型峰值相對其在樁頂鉸接和樁頂固定工況下的各階振型峰值較小，并且樁頂自由時，隨著黏彈性地基中單樁自振階數(shù)的增大，樁端的橫向位移逐漸減小。

3 結(jié)語

將回傳射線矩陣法推廣至樁土系統(tǒng)的振動分析中，對比分析了黏彈性地基中樁的外露長度、埋置深度、樁端約束情況對單樁自振頻率、衰減系數(shù)和模態(tài)的影響。得出以下結(jié)論：

（1）隨著黏彈性地基中單樁自振階數(shù)的增大，其自振頻率增幅、衰減系數(shù)逐漸增大，尤其是第四階之后，各階自振頻率增幅很明顯。

（2）隨著外露長度越長，結(jié)構(gòu)的各階自振頻率和衰減系數(shù)越?。浑S著埋置深度越深，結(jié)構(gòu)的各階自振頻率越小，衰減系數(shù)越大。

（3）埋置長度和外露長度對結(jié)構(gòu)的模態(tài)并無明顯影響，但樁長對振型中橫向位移的影響較大；埋置部分的橫向位移小于外露部分的橫向位移。

（4）樁頂固定工況下單樁各階自振頻率大于其在樁頂鉸接和樁頂自由工況下的各階自振頻率；樁頂自由工況下單樁各階衰減系數(shù)大于其在樁頂固定和樁頂鉸接工況下的各階衰減系數(shù)；樁頂自由工況下單樁各階振型峰值相對其在樁頂鉸接和樁頂固定工況下的各階振型峰值較小。