周鳳璽， 鄭 琦

(1. 蘭州理工大學(xué) 土木工程學(xué)院，蘭州 730050； 2. 西部土木工程防災(zāi)減災(zāi)教育部工程研究中心，蘭州 730050)

隨著城鎮(zhèn)化建設(shè)和現(xiàn)代工業(yè)的迅速發(fā)展，各種振動(dòng)問題日益頻繁，環(huán)境振動(dòng)已經(jīng)成為一種新形式的環(huán)境污染，并被列為當(dāng)前世界的七大環(huán)境公害之一。而隨著人們生活水平的提高，對(duì)生活環(huán)境和工作環(huán)境的質(zhì)量要求也越來越高，對(duì)振動(dòng)的限制越來越嚴(yán)格，環(huán)境振動(dòng)及其治理已經(jīng)成為巖土工程的一個(gè)研究熱點(diǎn)問題。

自Woods[1]通過一些現(xiàn)場原位試驗(yàn)，研究了近場主動(dòng)隔振和遠(yuǎn)場被動(dòng)隔振的隔振效果，并且提出了一個(gè)評(píng)判屏障隔振效果的重要參數(shù)振幅衰減系數(shù)以來，國內(nèi)外學(xué)者對(duì)連續(xù)屏障和非連續(xù)屏障的減振隔振效果進(jìn)行了一系列的研究工作。比如在振源或需保護(hù)結(jié)構(gòu)地基中設(shè)置隔振屏障(比如空溝、多排樁、空心管樁等)來降低其振動(dòng)水平是目前進(jìn)行地基減振隔振的主要措施[2-5]。由于受表面激振的層狀地基中波的傳播存在截止頻率，當(dāng)激振頻率低于截止頻率時(shí)，地基中不存在波的傳播現(xiàn)象?；诖爽F(xiàn)象，Schmid等[6]建議使用一個(gè)剛性層來形成有限尺寸的人工基巖，并將這個(gè)人工基巖稱為波阻板(Wave Impeding Block,WIB)。隨后他們采用二維頻域邊界元法，分析了WIB主動(dòng)隔振和被動(dòng)隔振，結(jié)果表明WIB的被動(dòng)隔振效果好于填充溝。采用邊界積分方程法，Peplow等[7]研究了二維雙層地基波阻板主動(dòng)隔振的隔振效果。利用邊界元和有限元法，Takemiya[8-9]研究了WIB的隔振效果?；趯訝罱橘|(zhì)中土與結(jié)構(gòu)動(dòng)力相互作用的半解析邊界元法模型，高廣運(yùn)等[10-11]分別分析了二維和三維層狀地基中波阻板的隔振效果。隨后，針對(duì)三維豎向非均勻和飽和地基模型，高廣運(yùn)等[12-13]對(duì)軌道交通荷載作用下豎向非均勻地基和飽和地基中波阻板的隔振性能進(jìn)行了研究。李寧等[14]通過現(xiàn)場試驗(yàn)，先研究了水平激振下混凝土材料的WIB對(duì)層狀地基的隔振效果，又分析了WIB的埋深對(duì)隔振效果的影響，并對(duì)比了實(shí)測的位移結(jié)果與半解析邊界元方法計(jì)算的結(jié)果?；贐iot多孔介質(zhì)波動(dòng)理論，結(jié)合PML邊界條件，高貝貝[15]用時(shí)域有限元方法研究了二維地基中儲(chǔ)液多孔波阻板的隔振性能?；贐iot多孔介質(zhì)波動(dòng)理論和功能梯度材料，周鳳璽等[16-17]研究了含液飽和多孔波阻板和梯度波阻板的地基振動(dòng)控制研究，結(jié)果表明含液飽和多孔波阻板和梯度波阻板的地基隔振體系更加具有優(yōu)越性。利用傅里葉級(jí)數(shù)展開的方法，引入輔助函數(shù)和波函數(shù)展開技術(shù)，張海等[18]研究了彈性半空間中含直邊界半圓形襯砌隧道對(duì)SH波的散射解析解。結(jié)合Green函數(shù)法及裂紋“切割”技術(shù)，采用保角映射的方法，楊在林等[19]研究了SH波下方入射時(shí)半無限空間界面附近橢圓形彈性夾雜及任意位置直線裂紋對(duì)SH波的散射問題。

由于WIB隔振設(shè)計(jì)屬于土與結(jié)構(gòu)動(dòng)力相互作用問題，由于其復(fù)雜性目前常用的分析方法主要是數(shù)值模擬(邊界元法、有限元法)和試驗(yàn)驗(yàn)證。但是，對(duì)于WIB隔振性能的解析分析較少。因此，本文借助復(fù)變函數(shù)的積分變換法將有限長度的波阻板變換為單位圓，建立了二維彈性地基中矩形波阻板對(duì)波場散射問題的控制方程。考慮平面SH波入射，運(yùn)用波函數(shù)展開法得到了二維地基波阻板隔振性能的解析解答。通過數(shù)值算例，分析了波阻板剪切模量、波阻板埋深、彈性波入射角度對(duì)隔振效果的影響。

1 WIB的保角映射變換

根據(jù)復(fù)變函數(shù)理論，通過映射函數(shù)z=ω(η)可將波阻板所在z平面上所占的區(qū)域變換為復(fù)平面η上的單位圓，保角映射前后的坐標(biāo)系統(tǒng)如圖1所示，其中H表示W(wǎng)IB埋置深度，α表示SH波入射角度，相應(yīng)的保角映射公式為[20]。

(1)

式中：

c1=cos(2χπ)

(2a)

(2b)

(2c)

(2d)

式中：χ為取決于WIB長寬比a/b的系數(shù)；系數(shù)c1∶c7僅與χ有關(guān)；R為與WIB尺寸有關(guān)的實(shí)數(shù)，可由以下兩式中的任意一式計(jì)算：

(3a)

(3b)

對(duì)于正方形，即長寬比a/b取1時(shí)，χ=0.25；而其他矩形，可采用試算法求χ，即估計(jì)一個(gè)χ值，由式(2)求得c1∶c7，代入式(1)可得映射函數(shù)，它是一個(gè)微曲的矩形。以寬度0.4 m，長度4.0 m的WIB為例，經(jīng)試算得到χ=0.106 5，圖2給出映射波阻板以及與單位圓的對(duì)應(yīng)關(guān)系，將矩形波阻板所在z平面上所占的區(qū)域變換為復(fù)平面η上的單位圓。

圖2 映射波阻板與單位圓的對(duì)應(yīng)關(guān)系

2 控制方程

穩(wěn)態(tài)平面SH波入射含WIB的二維彈性地基中，引起了散射波位移場為垂直xy平面的位移分量w(x,y)(省略了時(shí)間因子e-iωt)，且應(yīng)滿足Helmholtz方程：

(4)

應(yīng)力和位移的關(guān)系為：

(5a)

(5b)

(6)

(7a)

(7b)

而在復(fù)坐標(biāo)系中，應(yīng)力表達(dá)式為：

(8a)

(8b)

式中：θ為曲線坐標(biāo)系(r,θ,z)中r軸與直角坐標(biāo)系中x軸之間的夾角。

3 WIB對(duì)SH波的散射

3.1 入射波場

(9a)

(9b)

式中：w0為入射波的波幅；α為入射角；kS為波數(shù)。

由(7)可知w(i)和w(r)相對(duì)應(yīng)的應(yīng)力場為：

(10a)

(10b)

(10c)

(10d)

3.2 散射波場

(11)

由式(7)可得到散射波對(duì)應(yīng)的散射應(yīng)力場為：

(12a)

(12b)

同理，WIB結(jié)構(gòu)內(nèi)的散射波為：

(13)

所對(duì)應(yīng)的應(yīng)力場為：

Jn+1(kB|z|)(z/|z|)n+1e-iθ

(14a)

Jn+1(KB|z|)(z/|z|)n+1e-iθ}

(14b)

3.3 邊界條件

(15a)

(15b)

3.4 問題的解答

將式(9)，式(11)和(13)代入式(15a)，整理后可得

(16a)

將式(10)，式(12)和式(14)代入式(15b)，可得

(16b)

用e-imθ乘方程(16)的兩邊，并在區(qū)間(-π,π)上積分，則得到關(guān)于待定復(fù)系數(shù)An和Bn的理論解的無窮線性方程組：

(17)

其中：

利用截?cái)嘤邢揄?xiàng)的方法解出方程組(17)。就可以確定波函數(shù)中未知系數(shù)An和Bn，并將未知系數(shù)An和Bn代入式(9)、(11)、(13)，即求得彈性半空間介質(zhì)的位移場和WIB的位移場。

4 分析討論

為驗(yàn)證本文解答的正確性，選取王貽蓀[22]對(duì)Lamb問題的解析解答。利用突加力問題的閉合解求得了Lamb問題的解答，避免了對(duì)難以積分形式解的求積。取彈性地基密度ρs=2×103kg/m3，剪切波速為cs=200 m/s，固有圓頻率為ω=10 Hz，剪切模量為μS=80 MPa。對(duì)比結(jié)果見圖3。由圖3看出，本文計(jì)算值與文獻(xiàn)[22]結(jié)果相一致，說明本文的算法有效。

圖3 SH荷載激振地表振幅解答對(duì)比

對(duì)于WIB隔振效果的評(píng)價(jià)，Woods提出用振幅衰減系數(shù)AR來衡量[1]：

AR=有屏障時(shí)振幅/無屏障時(shí)振幅

且AR≤0.25時(shí)隔振效果較好。

為了分析級(jí)數(shù)項(xiàng)系數(shù)對(duì)數(shù)值計(jì)算結(jié)果的影響。選取參數(shù)a=4 m，b=0.4 m，ρS=2×103kg/m3，ω=10 Hz，μS=80 MPa，kS=0.05，kS/kB=1/3，μB/μS=2，H=0.5b，α=30°。圖4給出了不同n的情況下地表豎向位移振幅的變化曲線。由圖4知，當(dāng)散射波系數(shù)n取10時(shí)可以滿足計(jì)算精度要求。

分析剪切模量比對(duì)隔振效果的影響，選取a=4 m，b=0.4 m，ρS=2×103kg/m3，ω=10 Hz，μS=80 MPa，kS=0.05，kS/kB=1/3，H=0.5b，α=30°。圖5繪出了地表振幅衰減系數(shù)的變化曲線。由圖5知，當(dāng)WIB的剪切模量與地基的剪切模量比為μB/μS=0.5時(shí)，采用WIB隔振體系隔振效果不是很明顯，當(dāng)剪切模量比為μB/μS=1時(shí)，可以看出WIB有顯著的隔振效果；當(dāng)μB/μS=2時(shí)，隔振效果顯著增強(qiáng)。所以，應(yīng)根據(jù)隔振設(shè)計(jì)實(shí)踐的需求，合理地選擇WIB的剪切模量可以有效提高隔振效果。

圖4 豎向位移振幅的變化曲線

圖5 不同WIB剪切模量時(shí)振幅衰減系數(shù)AR變化曲線

其余參數(shù)不變，取模量比μB/μS=2時(shí)，圖6給出了WIB的埋深對(duì)地表振幅衰減系數(shù)的變化曲線。從圖6可以發(fā)現(xiàn)，隨著埋深的增大對(duì)地表振幅的隔振效果減小。

圖6 不同WIB埋深時(shí)振幅衰減系數(shù)AR變化曲線

為了討論了不同的入射角度變化對(duì)地表幅值的影響。圖7給出了kS/kB=1/3，H=0.5b時(shí)地表振幅衰減系數(shù)的變化曲線。由圖7知，地震波入射角度α對(duì)隔振效果有著明顯的影響，隨著入射角度的增大隔振效果明顯減小，當(dāng)α=90°時(shí)，隔振效果減小尤為明顯。

圖7 不同入射角度時(shí)振幅衰減系數(shù)AR變化曲線

5 結(jié) 論

運(yùn)用復(fù)變函數(shù)的保角映射方法和波場勢函數(shù)展開法，得到二維彈性地基WIB隔振的解析解答，通過參數(shù)分析討論了WIB剪切模量、埋深、彈性波入射角度對(duì)隔振效果的影響，研究表明：隨著WIB剪切模量的增大其隔振效果明顯增加，隨著WIB埋置深度增加隔振效果逐漸減小，并且WIB的隔振效果隨彈性波入射角度的增大而減小。