姜博龍，劉維寧，馬 蒙，李明航

(北京交通大學(xué) 土木建筑工程學(xué)院，北京 100044)

隨著軌道交通路網(wǎng)的加密，運行列車引起的環(huán)境振動問題也愈加凸顯。它對沿途居民的日常生活、精密儀器的正常運轉(zhuǎn)和文物古建的保護(hù)都造成了一定的影響。當(dāng)軌道減振和敏感建筑被動隔振無法徹底解決振動問題時，作為綜合減隔振措施中重要一環(huán)的傳播路徑隔振具有不可替代性。

排樁作為典型的非連續(xù)波屏障，應(yīng)用在傳播路徑隔振中，憑借其造價低、耐久不易損壞、隔振效果好等優(yōu)勢，受到國內(nèi)外學(xué)者廣泛關(guān)注，并針對其展開了大量研究?，F(xiàn)有的針對隔振排樁的研究集中在波場求解方法方面，通常采用二維簡化模型對體波進(jìn)行求解而采用三維數(shù)值模型對表面波進(jìn)行求解[1-5]。這些研究考慮到的樁數(shù)目較少或僅僅是單排，且忽視了其周期性。直到石志飛等[6]充分考慮了群樁的周期性，將周期結(jié)構(gòu)帶隙理論引入到工程隔振領(lǐng)域，為解決環(huán)境振動問題提供了新的思路和依據(jù)。周期結(jié)構(gòu)對通過它的振動波具有選擇透過性，無法透過其傳播的振動波所處的頻帶被稱為禁帶或帶隙[7]。黃建坤[8-10]采用平面波展開法求解周期排樁帶隙，以首階完全禁帶的起、止頻率以及對應(yīng)的帶寬作為評價指標(biāo)，對周期排樁首階帶隙分布影響因素進(jìn)行了分析；劉心男[11]對比了二維和三維模型中周期排樁對帶隙范圍內(nèi)振動的衰減結(jié)果，發(fā)現(xiàn)隨著樁長的增加，三維數(shù)值結(jié)果逐步收斂于二維數(shù)值結(jié)果；蒲興波[12]提出了一種新的表面波識別方法，用于獲取二維土—周期排樁系統(tǒng)的面波帶隙。

目前對于周期排樁的帶隙描述、影響因素分析和作用效果評價主要是采用首階禁帶的起、止頻率以及對應(yīng)的帶寬3個指標(biāo)。周期排樁的隔振效果不僅與初始帶隙相關(guān)，與各階帶隙的分布也有關(guān)系。另外，目前周期排樁帶隙影響因素分析和結(jié)構(gòu)優(yōu)化選型都沒有充分考慮被阻隔對象特征，造成了結(jié)構(gòu)優(yōu)化選型針對性不強。因此，針對軌道交通環(huán)境傳播路徑隔振，本文提出以周期排樁帶隙性能評價函數(shù)Φ為單一指標(biāo)、并考慮帶隙分布的選型方法。

1 隔振周期排樁帶隙計算

對周期排樁帶隙的求解可根據(jù)固體物理晶格理論[13]進(jìn)行簡化：周期排樁可由某一根單樁通過周期性和對稱性的拓?fù)渥儞Q獲得，因此，計算時僅需以某一單樁作為基本單元，并對其施加周期邊界進(jìn)行研究。根據(jù)周期結(jié)構(gòu)帶隙產(chǎn)生機(jī)理的不同，隔振周期排樁可分為散射型排樁和局域共振型排樁，二者對應(yīng)的基本單元如圖1所示，圖中：A，B，C分別代表土體、樁和填充體；R為樁半徑；r為填充體半徑；a為周期常數(shù)。

圖1 排樁基本單元平面示意圖

下面利用平面波展開法[7]進(jìn)行周期排樁帶隙計算。土—周期排樁體系彈性動力學(xué)基本(維納)方程為

(1)

其中，

u(d,t)=(uxuyuz)

式中：t為時間；d為位置矢量；u(d,t)為位移矢量；ux,uy,uz為沿x,y,z方向的位移分量；ρ(d)為密度函數(shù)；λ(d)和μ(d)為與位置矢量d相關(guān)的Lamé常數(shù)。

由于密度函數(shù)ρ(d)，與位置矢量d相關(guān)的Lamé常數(shù)λ(d)和μ(d)在空間分布呈現(xiàn)周期性，故其可在倒格矢空間按傅里葉級數(shù)展開，為敘述方便，統(tǒng)一采用f來表示這3個材料參量，有

(2)

其中，

式中：i為虛數(shù)單位；G1為對應(yīng)各組元材料參數(shù)分布場內(nèi)的倒格矢；f(G1)為傅里葉系數(shù)；S為基本單元面積。

對于散射型周期排樁，有

(3)

其中，

式中：η為樁占整個基本單元的填充比；P(G)為結(jié)構(gòu)函數(shù)； J1為第一類第1階的貝塞爾函數(shù)。

對于局域共振型周期排樁，有

(4)

其中，

式中：ηB，ηC，ηf分別表示樁、填充體以及散射體(填充體+樁)占整個基本單元的填充比；P(G1)，Ph(G1)為結(jié)構(gòu)函數(shù)。

式(1)的解為

(5)

式中：ω為角頻率；G2為位移場中的倒格矢；K為限制在第一Brillouin區(qū)內(nèi)的Bloch波矢；uK(G2)為與各材料參數(shù)具有相同周期的Bloch調(diào)幅函數(shù)。

將式(2)和式(5)代入式(1)，可得本征方程為

(6)

G2)(K+G2)(K+G3)uz(K+G2)

(7)

其中，

μG3-G2[(K+G2)j(K+G3)l]}

j=x,y；l=x,y

G3=G1+G2

式中：ux(K+G2)，uy(K+G2)和uz(K+G2)為沿x，y和z方向的位移分量；δjl為克羅內(nèi)克函數(shù)。

令波矢K掃略第一不可約布里淵區(qū)便可得到整個結(jié)構(gòu)的頻散曲線和帶隙特性[7]。

下面用1個具體算例驗證該計算方法，選取文獻(xiàn)[14]中的周期排樁。該排樁為圓截面實心樁結(jié)構(gòu)，拓?fù)湫问綖檎骄Ц?，排樁類型為散射型，周期常?shù)a=2 m，樁的半徑R=0.65 m，具體材料參數(shù)參見表1。

表1 土—樁系統(tǒng)的材料組成及物理力學(xué)參數(shù)

圖2給出了本文和文獻(xiàn)[14]計算帶隙的結(jié)果。二者均在33.3～37.3 Hz存在1個完全帶隙。二者計算結(jié)果的吻合證明了本文方法的正確性。

圖2 本文和文獻(xiàn)[14]計算的帶隙結(jié)果

2 能帶結(jié)構(gòu)帶隙性能評價函數(shù)

圖3給出了2種典型周期排樁的帶隙計算結(jié)果，圖3(a)為按照六角晶格布置的散射型排樁，周期常數(shù)為4 m，混凝土樁半徑為1.2 m；圖3(b)為按照六角晶格布置的局域共振型排樁，周期常數(shù)為3 m，混凝土樁半徑為1 m，橡膠填充體半徑為0.4 m。

目前對于周期排樁的帶隙描述、影響因素分析和作用效果評價主要是采用首階禁帶的起、止頻率以及對應(yīng)的帶寬3個指標(biāo)。由圖3可見，散射型排樁的首階帶隙起始頻率低于局域共振型排樁，且首階完全帶隙的帶寬要寬于局域共振型排樁；但局域共振型排樁80 Hz以下帶寬總和卻大于散射型排樁。因此，對于周期排樁的帶隙描述、影響因素分析和作用效果評價不能只關(guān)注于首階完全帶隙，也要綜合考慮關(guān)心頻段內(nèi)各階帶隙的分布情況。

圖3 典型六角晶格布置排樁的帶隙分布對比

帶隙具備2個典型特征，一個是中心頻率，另一個是帶寬。圖4為被阻隔對象顯著頻率與周期隔振結(jié)構(gòu)帶隙分布關(guān)系示意圖。被阻隔對象的顯著頻率為63 Hz，圖中帶隙1和帶隙2分別為2種周期隔振結(jié)構(gòu)的帶隙分布。從圖4可以看出，與帶隙2相比，具有帶隙1的周期隔振結(jié)構(gòu)的帶隙中心頻率更加接近被阻隔對象振動的顯著頻率，但是帶寬卻較窄。因此若想取得良好的隔振效果，需要使周期隔振結(jié)構(gòu)的帶隙中心頻率貼近被阻隔對象的顯著頻率，且?guī)捿^寬。

圖4 被阻隔對象顯著頻率與周期隔振結(jié)構(gòu)帶隙分布關(guān)系示意圖

基于上述分析，提出一種綜合考慮被阻隔對象振動特征、隔振排樁各階帶隙分布的周期排樁帶隙性能評價函數(shù)Φ為

(8)

式中：n為帶隙數(shù)目；ωj為出現(xiàn)帶隙的下限頻率；ωj+1為同一帶隙的上限頻率；(minωj+1-maxωj)為該帶隙的帶寬；(minωj+1+maxωj)/2為該帶隙的中心頻率；ωD為需要阻隔的關(guān)心頻率。

隨著帶隙中心頻率逐步逼近需要阻隔的關(guān)心頻率ωD，|ωD-(minωj+1+maxωj)/2|將取得極小值，同時隨著帶隙(minωj+1-maxωj)的增大，周期排樁帶隙性能評價函數(shù)Φ將取得極大值，Φ值越大表征該結(jié)構(gòu)形式越適用于被阻隔對象的顯著頻率的振動阻隔。

3 軌道交通特征頻率隔振周期排樁選型

本文以東單地鐵站到建國門地鐵站區(qū)間隧道列車運行引起環(huán)境振動為研究背景[16]進(jìn)行傳播路徑隔振周期排樁選型。

該段地鐵隧道初襯厚度為0.25 m，所采用的混凝土的彈性模量為28.5 GPa，泊松比為0.2，密度為2 450 kg·m-3，材料阻尼比為0.02。二襯厚度為0.3 m，二襯混凝土的彈性模量為30 GPa，泊松比為0.2，密度為2 500 kg·m-3，材料阻尼比為0.02。隧道基底混凝土參數(shù)與二襯一致，列車采用常見的6節(jié)編組，經(jīng)過測試斷面的速度為60 km·h-1，扣件為普通DTⅥ2扣件。隧道埋深與地層相對位置關(guān)系如圖5所示，地層信息參考表2。

圖5 隧道與地層相對位置[16]

土體類型厚度/m剪切波速/(m·s-1)壓縮波速/(m·s-1)動彈性模量/(MN·m-2)動泊松比密度/(kg·m-3)阻尼比雜填土 2.5173.3356.5138.80.3501 6500.05粉質(zhì)黏土 18.0235.0472.6308.00.3272 0100.04卵石和圓礫332.0620.0670.00.2862 0500.03

對該斷面進(jìn)行地表環(huán)境振動監(jiān)測時測得的距離隧道中心線40 m位置處的地表豎向振動響應(yīng)如圖6所示。由圖6可以看出，DTVI2扣件軌道型式下地鐵列車以60 km·h-1速度運行引起的地表豎直方向振動峰值出現(xiàn)在63 Hz，此頻率為軌道和扣件系統(tǒng)的自振頻率。因此以該豎直方向振動的峰值頻率63 Hz作為關(guān)心頻率，即ωD=63 Hz。

隔振周期排樁選取按照正方晶格和六角晶格2種典型晶格類型布置的散射型排樁和局域共振型排樁，其布置如圖7所示。

根據(jù)文獻(xiàn)[15]可知，排樁埋深一般到粉質(zhì)黏土層底部可取得較好效果，鑒于隧道也埋置于粉質(zhì)黏土層，且排樁所在地層中粉質(zhì)黏土的含量占比為主(約90%)，本文計算以粉質(zhì)黏土—周期排樁的帶隙性能評價函數(shù)Φ作為選型依據(jù)。相應(yīng)材料參數(shù)見表3。

圖6 距離隧道中線40 m處地表振動加速度響應(yīng)測試值

組成材料彈性模量/MPa密度/(kg·m-3)泊松比A土體308.01 9000.327B混凝土30 000.02 3000.180C填充體(橡膠)7.81 3000.470

選取14種周期排樁進(jìn)行優(yōu)化選型，14種周期排樁的具體情況及對應(yīng)的帶隙性能評價函數(shù)Φ值見表4。由表4可見， 5號周期排樁的Φ值最大，說明5號周期排樁最適宜用于隔離以豎直方向63 Hz為關(guān)心頻率的環(huán)境振動。

圖7 不同晶格類型布置的排樁

表4 各周期排樁的具體情況及對應(yīng)的帶隙性能評價函數(shù)Φ值

圖8給出了5號周期排樁分別在雜填土和粉質(zhì)黏土中的帶隙分布圖。從圖8可以看出，5號周期排樁在其穿過的2層土體中都有帶隙覆蓋63 Hz的關(guān)心頻率，且在80 Hz以下帶寬分布范圍較大。

圖8 5號周期排樁豎直方向帶隙分布

4 結(jié)果驗證

運用ABAQUS軟件，以前文提到的工程資料建立隧道—地層—周期排樁耦合的三維動力有限元分析模型(圖9)進(jìn)行軌道交通環(huán)境振動傳播路徑隔振周期排樁的隔振效果分析，驗證第3節(jié)中選型結(jié)果的有效性。地層參數(shù)見表2，模型尺寸為100m×60 m×24 m，土體、隧道及排樁均采用實體單元模擬，在隧道基底上輸入的激勵為采用周期性車輛—軌道耦合解析模型[16]計算得到的扣件反力時程，邊界設(shè)置為彈簧阻尼吸收邊界。

圖9 三維動力有限元模型

采用加速度振幅衰減系數(shù)Ar[17]評價排樁隔振性能為

(9)

式中：a1和a2分別為隔振排樁后某點的加速度幅值與同一位置處無隔振排樁的加速度幅值。

Ar→0表明結(jié)構(gòu)具有良好的隔振性能。

圖10為5號周期排樁后距離隧道中心線40 m位置(與實測點對應(yīng))處的振幅衰減系數(shù)及其帶隙分布。

圖10 距離隧道中心線40 m處振動加速度振幅衰減系數(shù)及其帶隙分布

由圖10可以看出，帶隙范圍內(nèi)2層土體的振動都有較好的衰減，說明5號排樁隔振效果明顯，也間接證明了帶隙計算結(jié)果的正確性與選型方法的有效性。

5 結(jié) 語

本文提出了基于帶隙性能評價函數(shù)的軌道交通特征頻段隔振周期排樁的選型方法，性能評價函數(shù)作為單一評價指標(biāo)綜合考慮了帶隙帶寬、中心頻率2個因素，規(guī)避了多指標(biāo)分析的不確定性；其次，將軌道交通環(huán)境振動的顯著頻率作為關(guān)心頻率引入該函數(shù)，充分考慮了被阻隔對象的振動特征。該方法不再單獨以首階帶隙為選型依據(jù)，綜合考慮了各階帶隙的分布情況。本文通過建立三維動力有限元模型對選型結(jié)果進(jìn)行了驗證，選出的隔振周期排樁在帶隙分布范圍內(nèi)有較好的隔振效果。