埋置塊式基礎地基阻抗函數的簡化計算方法研究

胡燦陽，陳清軍，徐慶陽，熊琴琴

(1.南京審計學院 江蘇省公共工程審計重點實驗室，南京 210029;2.同濟大學 土木工程防災國家重點實驗室，上海 200092)

埋置塊式基礎地基阻抗函數的簡化計算方法研究

胡燦陽1，2，陳清軍2，徐慶陽1，熊琴琴1

(1.南京審計學院 江蘇省公共工程審計重點實驗室，南京 210029;2.同濟大學 土木工程防災國家重點實驗室，上海 200092)

目前關于地基阻抗的分析和計算方法都比較繁瑣，不易為工程設計人員所掌握。基于此，提出了埋置基礎地基阻抗函數的簡化計算方法。首先，采用錐體模型分析均質半空間內埋置圓盤一維振動波的傳播過程，得到了平動和轉動時的地基阻抗。在此基礎上，利用埋置圓盤組來模擬埋置基礎，分析得到了埋置基礎地基阻抗函數。最后，通過算例表明此方法計算效率高，誤差小，適用于工程計算。

錐體模型;塊式基礎;阻抗函數;埋置基礎;土－結構相互作用

土－結構動力相互作用是新近發展起來的與多學科相互交叉和滲透的學科，描述其動力學行為的數學模型(通常為一組控制方程)十分復雜。由于實際問題的復雜性和計算過程的繁瑣性，使得這些模型不便用于實際工程中。近年來，一種采用錐體模型求解土－結構相互作用中地基阻抗函數的方法，因其簡便性和實用性，受到了重視，但這一方法目前尚處于起步階段，需要做進一步的研究。

錐體模型最早是由Ehlers［1］在1942年提出，最初用它來分析明置圓盤的平動，繼而也用它分析轉動。Gazetas［2］和 Gazetas 和 Dobry［3］使用楔形體和錐體來闡述二維和三維情況下的輻射阻尼。Meek和Wolf使用錐體模型得到了均質半空間表面基礎動力響應的簡便計算方法。他們又分別計算了剛性基巖覆蓋層表面基礎的動力響應，分析了柔性基巖覆蓋層表面基礎的動力響應。1993 年，Meek 和 wolf［4］利用 Boussinesq 等理論詳細說明了錐體模型的物理意義，討論了利用錐體模型代替彈性半空間的可行性，從而消除了眾多學者對錐體模型的疑慮。此外，2002年，Jaya和 Prasad［5］使用了同樣的截頭錐段研究了層狀土內埋置基礎的動力剛度。然而，這種截頭錐段法也存在缺陷，它的主要問題是:當低頻時，阻尼系數會變成負值，這在物理上是不可能的。文獻［6］對于錐體模型的初步應用有詳細地總結。2003年，Pradhan等［7］利用波在錐中的傳播，計算了層狀土表面圓形基礎的動力阻抗，這就克服了利用上面截頭錐段分析時的缺陷。2008年，Pradhan等［8］利用錐模型對于層狀土表面機器基礎振動進行了分析，并同時進行了實驗研究，通過對比，驗證了錐體模型對于機器基礎振動分析的有效性。目前，國內對錐體模型的研究較少，張劍寒等［9，10］基于錐體模型概念，得到了層狀半空間地基上明置基礎阻抗函數的計算方法，并與薄層法進行對比來研究它的精度。

1 埋置圓盤

可以把地基－基礎接觸面看作一個無質量剛性板，下面的均質半空間地基近似看作具有不同頂點高度z0的截頭半無限彈性錐體，分為平動和轉動錐體，表面圓盤相對應的錐體模型頂高半徑比等參數以及平動錐模型與轉動錐模型中波反射系數和相應的阻抗函數見文獻［6］。

對于埋置在均質半空間里的半徑為r0，埋深為e的圓盤，利用錐模型可以得到阻抗函數S(ω)的閉合解。先分析如圖1所示豎向自由度的情況。使用雙邊錐模型，每個初始錐的動力剛度都相同。圓盤的初始位移幅值(ω)與圓盤上的作用力幅值 P0(ω)的關系為［11］:

式中:K是半空間表面圓盤的靜剛度系數;c是圓盤振動波傳播速度，對于水平振動和扭轉振動用剪切波速，對于豎向振動和扭轉振動用脹縮波速;ω為頻率。

圖1 埋置在均質半空間里圓盤波傳播過程Fig.1 Wave propagation from disk embedded in homogeneous half-space

波在圓盤下面的初始錐內(頂點為1)向下傳播，沒有反射波產生，錐的張角也即頂點高半徑比z0/r0，可以通過均質半空間表面圓盤的錐模型靜力剛度系數與相應的用三維彈性理論求得的半空間上圓盤的靜剛度精確解相等求出［6］，不同的平截錐具有不同的張角。波在圓盤上面的初始錐內(頂點為2)向上傳播時，在自由表面產生反射，自由表面處入射波幅值為:

在自由表面處，反射系數 －α(ω)=1，反射波幅值為:

反射波在頂點為3的新錐內傳播。當它傳播到圓盤所在深度時，幅值衰減為:

因為半空間是均質的，這個波會繼續在截頭錐內向無限遠處傳播，而不再產生反射。圓盤處的總位移u0(ω)是由這個反射波［式(4)］和原始激勵波［式(1)］疊加而成的。于是有幅值:

引入幾何參數k=2e/z0和傳播的時間參數T=2e/c，后得到:

式(7)也即為平動錐的阻抗函數，對于水平振動和豎向振動，分別取相應的K、T和k進行計算。

對于轉動自由度，分析過程是類似的。埋置在均質半空間里的圓盤地基阻抗函數:

同樣，對于搖擺振動和扭轉振動，取相應的K?、T和k進行計算。

2 埋置圓盤組

如圖2所示，為了表示在半空間里半徑為r0，深度為e的圓柱形基礎，可以在基礎位置處放置m個圓盤。也就是說，在以后將要開挖的地方可以看作是m個圓盤和其中的 m－1個分層。第一個圓盤在層狀半空間的表面，第 m個圓盤在埋置基礎的剛性基底處。文獻［5］建議兩個相鄰圓盤的最大豎向間距Δe應該滿足以下要求:每個波長度內至少要分為6層。

圖2 半空間內埋置圓盤組Fig.2 Stack of disks to represent cylindrical foundation embedded in half-space

為了精確，本文采用每個波長度內要分為10層，波長為λ=c/(ω/2π):

式中ω代表了動力模型的最高頻率，c是波速(通常是剪切波速cs)。當使用無量綱頻率a0=ωr0/cs時，有:

對于豎向自由度(圖3)，有:

式(10)給出了埋置在未開挖的層狀半空間內圓盤的力－位移關系，這里［Sf(ω)］是指與自由場圓盤組有關的阻抗矩陣。在自由場中，幅值為{u(ω)}的豎向位移分量之間是相互獨立的。在剛性基礎中，各點的豎向位移和基底的豎向位移相同，幅值都為u0(ω)。這種剛體位移表示為:

運動約束矢量{A}={111…1}T。力Q0(ω)等于各個荷載Pi(ω)之和:

則，由式(10)、式(11)、式(12)得:

扭轉自由度的分析過程類似，埋置基礎地基阻抗是:

這里m?是開挖區域的質量極慣性矩，{A}和上面定義相同。

圖3 豎向自由度時剛體位移和相互作用力Fig.3 Enforcement of rigid-body displacement for verticaldegree of freedom

分析耦合水平和搖擺自由度(圖4)。埋置圓盤的水平位移幅值{u(ω)}和轉動幅值{?(ω)}都和基底的剛體運動相關，基底的水平位移幅值為{u0(ω)}，轉動幅值為{?0(ω)}，則有:

于是有埋置基礎的阻抗函數矩陣:

對上面的過程很容易進行一定的擴展。比如，當基礎墻沒有剛度和剛體約束時，矩陣［A］不再需要了。在這種情況下，位移和轉動幅值［u(ω)］和［?(ω)］都是未知的，由這些自由度來確定埋置基礎的地基阻抗函數。圓盤中間的區域仍然通過對內部域的動力剛度矩陣進行抽取進行分析。后者是(上標e表示開挖):

其中，［K］和［M］分別為開挖部分相應于{u(ω)}和{?(ω)}的靜力剛度矩陣和質量矩陣，可以通過將內部域模型化為具有豎軸的桿或者梁得到。修改后的相互作用力－位移關系為:

埋置基礎的地基阻抗函數矩陣等于自由場阻抗矩陣與開挖部分剛度矩陣的差:

3 算例

如圖5所示，在泊松比為υ=0.25的均質半空間里，有半徑為r0埋置深度為e的剛性圓柱基礎，計算它的所有自由度的阻抗函數。埋深率e/r0分別為0.5、1和2(這個范圍包含了大多數實際工程的情況)。結果根據無量綱頻率a0=ωr0/cs范圍從0到4來確定。

由式(9b)可知兩個相鄰圓盤的最大豎向距離為Δe=(π/20)r0，這時a0=4。對于最大埋深率e/r0=2時，至少要分為13層，也就是要14個圓盤。本文使用16層，因此 Δe==0.125r0。同理，當 e/r0=1 時，模型分為8層，當e/r0=0.5時，模型分為4層。

使用錐模型時，埋置基礎靜力剛度系數的計算也可以用同樣的過程，這時ω=0。它們除以相同運動下明置基礎的相應靜力剛度系數，得到無量綱的靜力剛度因子，并且在圖6中和文獻［12］的數值解進行對比。圖6中的曲線相應于文獻［13］中的經驗公式，它適用于e/r0≤2，由圖上可見錐模型的精度對于所有自由度都是很好的。

圖7到圖10給出了水平、豎向、搖擺和扭轉運動阻抗函數對應于無量綱頻率a2的彈性系數和阻尼系數。由圖上可見，水平運動時，動力彈性系數的錐模型結果與文獻［12］的數值解在低頻部分吻合很好，在高頻部分存在很小的偏差，動阻尼系數的錐模型結果與文獻［12］的偏差不超過10%，在豎向、搖擺和扭轉的情況下，錐模型的結果與文獻［12］的值吻合的很好(這里是埋深率e/r0=1的情況，也可以在其它情況下進行比較)。由此可見，使用錐模型分析均質半空間里的埋置基礎阻抗函數具有很高的精度。

圖10 扭轉阻抗函數(e/r0=1，υ=0.25)Fig.10 Torsional impedance function(e/r0=1，υ =0.25)

4 結論

本文首先基于錐體模型，對振動波在錐段內的傳播進行分析，得到了埋置無質量剛性圓盤的地基阻抗函數。進而用一組埋置圓盤來模擬埋置基礎，可以得到地基阻抗函數。最后，通過算例表明，此方法物理意義明確、分析和計算過程簡單，結果具有很高的精度，適于工程人員掌握和應用。

［1］Ehlers G.The effect of soil flexibility on vibrating systems［J］.Beton Eisen，1942，41(21/22):197 －203.

［2］Gazetas G. Simple physical methods for foundation impedance.Dynamic behaviour of foundations and buried structures［M］.In:Banerjee PK，Butterfield R，editors.Developments in soil mechanics and foundation engineering，vol.3.London:Elsevier，1987.［chapter 2］.

［3］Gazetas G，Dobry R.Simple radiation damping model for piles and footings［J］.Journal of Geotechnical Engineering，ASCE，1984，110(6):937 －956.

［4］Meek J W，Wolf J P.Why cone models can represent the elastic half-space［J］.Earthquake Engineering and Structural Dynamics，1993，22(1):759 －771.

［5］Jaya K P，Prasad A M.Embedded foundation in layered soil under dynamic excitations［J］.Soil Dynamics and Earthquake Engineering，2002，22(6):485 －498.

［6］Wolf J P.Foundation vibration analysis using simple physical models［M］.Englewood Cliffs，NJ:Prentice-Hall，1994.

［7］Pradhan P K，Baidya D K，Ghosh D P.Impedance functions of circular foundation resting on layered soil using cone model［J］.Electronic Journal of Geotechnical Engineering，2003，8(B).

［8］ Pradhan P K，Mandal A，Baidya D K，et al.Dynamic response of machine foundation on layered soil:cone model versus experiments［J］. Geotechnicaland Geological Engineering，2008，26(4):453 －468.

［9］張劍寒，蔣 通.層狀地基上塊式基礎地基阻抗函數簡化計算方法［J］.工業建筑，2007，37(8):68－72.

［10］張劍寒，劉志杰，蔣 通.層狀地基對塊體基礎地基剛度的影響［J］.工業建筑，2008，38(s1):664 －666.

［11］ Meek J W，Wolf J P.Cone models for embedded foundation［J］.Journal of Geotechnical Engineering Division，ASCE，1994，120(1):60－80.

［12］ Aspel R J，Luco J E.Impedance functions for foundations embedded in a layered medium: an integralequation approach［J］. Earthquake Engineering and Structural Dynamics，1987，15(2):213 －231.

［13］ Pais A， KauselE. Approximateformulasfordynamic stiffnesses of rigid foundations［J］.Soil Dynamics and Earthquake Engineering，1988，7(4):213 －227.

Simplified calculation method for impedance function of embedded block foundation

HU Can-yang1，2，CHEN Qing-jun2，XU Qing-yang1，XIONG Qin-qin1

(1.Nanjing Audit University，Jiangsu Key Laboratory of Public Project Audit，Nanjing 210029，China;2.State Key laboratory for Disaster Reduction in Civil Engineering，Tongji University，Shanghai 200092，China)

Analysis and calculation methods for impedance function of a foundation are complicated，and it is difficult for engineering designers to grasp and use them.A simplified method was proposed here to calculate impedance function of a foundaton embedded in a half-space.Firstly，based on a cone model，one-dimensional vibration wave propagating of a rigid massless circular disk embedded in a homogeneous half-space was analysed.Impedance functions of translational and rotational cone models were derived.Then，the embedded foundation was modeled with a stack of disks and its impedance function was calculated.The model was validated by comparing the results with the rigorous analysis ones. The numerical results indicated that this method has higher precision and efficiency，it is applicable for engineering calculation.

cone model;block foundation;impedance function;embedded foundation;soil-structural interaction

TU470+.3

A

科技部國家重點實驗室基礎研究資助項目(SLDRCE08－B－03);江蘇省高校自然科學研究項目(10KJB560003);江蘇省高校優勢學科建設工程資助項目(PADD)

2009－12－21 修改稿收到日期:2010－04－06

胡燦陽 男，博士，講師，1974年1月生