梁玉旺，周鳳璽，2，郝磊超

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

引言

我國工業(yè)、交通等方面飛速發(fā)展，在提高人們生活質(zhì)量的同時，各種工業(yè)活動及交通運(yùn)行誘發(fā)的人工振動問題對建筑物、精密儀器以及居民正常生活等產(chǎn)生了嚴(yán)重影響。甚至國際上已將其列為七大公害之一。因此，人工振動的治理已經(jīng)成為巖土工程領(lǐng)域的研究熱點(diǎn)。

基巖上單一土層的振動存在截止頻率，當(dāng)激振頻率低于截止頻率時，土層中沒有波的傳播?；诖嗽?，Chouw等［1-2］首先提出在地基中人工設(shè)置一個硬夾層，以形成有限尺寸的人工基巖來達(dá)到隔振目的，將其命名為波阻板（Wave Impeding Block），并通過二維頻域邊界元法對其主動和被動隔振效果進(jìn)行了分析。Peplow等［3］采用邊界積分方程法研究了雙層地基中波阻板隔振效果，由于此方法比較繁瑣，并不適用于分析多層地基振動問題。因此，高廣運(yùn)等［4-7］結(jié)合薄層法—邊界元法（半解析邊界元法）分析了二維、三維層狀地基中波阻板隔振效果，之后，通過現(xiàn)場試驗(yàn)對半解析邊界元法進(jìn)行了驗(yàn)證，得到了較為滿意的結(jié)果［8］。

此外，還有部分學(xué)者研究了不同類型波阻板隔振效果：Takemiya［9］提出了一種蜂窩狀波阻板，并通過現(xiàn)場試驗(yàn)分析了其對高架鐵路樁基振動的隔振效果?；谟邢拊椒?，李志華等［10］驗(yàn)證了蜂窩狀波阻板相對于實(shí)體波阻板具有更好的隔振效果。周鳳璽等［11-12］研究了含液飽和多孔波阻板和梯度波阻板的地基振動控制研究，結(jié)果表明，這2種類型的波阻板更具優(yōu)越性。

然而，波阻板邊界多為不規(guī)則形狀，邊值問題比較復(fù)雜。劉殿魁等［13］提出“域函數(shù)”的概念，即基于復(fù)變函數(shù)理論的保角映射方法將不規(guī)則邊界轉(zhuǎn)換為單位圓邊界，實(shí)現(xiàn)了邊值問題的簡化，為求解彈性波散射問題［14-16］提供了一種有效方法。

文中基于保角映射方法建立幾何模型和邊界條件，考慮穩(wěn)態(tài)平面Rayleigh波入射，運(yùn)用波函數(shù)展開法求解二維地基中埋設(shè)波阻板時對Rayleigh 波的散射解答。通過數(shù)值算例，對埋置深度、剪切模量比、入射頻率等參數(shù)對波阻板隔振效果的影響進(jìn)行了分析。

1 保角映射

文中采用復(fù)變函數(shù)理論中的保角映射方法，引入保角映射函數(shù)z=ω(η)，將z 平面中矩形截面波阻板邊界和外域映射為η平面中單位圓邊界和外域，保角映射坐標(biāo)系如圖1 所示，其中，坐標(biāo)系z=x+iy，η=ξ+iζ，取單位圓η=eiθ。

圖1 波阻板保角映射坐標(biāo)系Fig.1 Conformal mapping coordinate system of wave impeding block

根據(jù)Schwarz-Christoffel 積分定理，矩形波阻板保角映射函數(shù)的表達(dá)式為［17］：

式中：R為由波阻板尺寸決定的實(shí)數(shù)；c1～c7為波阻板形狀決定的復(fù)常數(shù)；系數(shù)k與波阻板的長寬比有關(guān)。

2 波阻板對Rayleigh波的散射解答

2.1 入射波場

將土體視為各向同性的均質(zhì)線彈性體，在半空間地基中有限深度處埋設(shè)矩形截面WIB，幾何模型如圖2所示。a、b對應(yīng)表示W(wǎng)IB 的長度和寬度，H表示W(wǎng)IB 幾何中心與水平地表之間的距離。坐標(biāo)系xoy和x1o1y1分別對應(yīng)復(fù)數(shù)坐標(biāo)，且復(fù)數(shù)坐標(biāo)量z=z1+iH，。

圖2 波阻板幾何模型Fig.2 Geometric model of wave impeding block

Rayleigh 波是縱波和橫波在自由表面處耦合而形成的面波，沿半空間表面?zhèn)鞑?，不考慮時間諧和作用（省去時間因子e-iωt），在直角坐標(biāo)系中入射波勢函數(shù)可表達(dá)為：

式中，Aφ和Bψ分別代表縱波和橫波的幅值，且滿足，kp、ks、kr分別為P波、SV波、Rayleigh波波數(shù)。

在z平面，入射波勢函數(shù)為：

這件事后來我們嶺北人都知道了，大家都覺得馱子說得有道理，這個周小羽性格孤僻，目無尊長，讀書一點(diǎn)讀不來，亂七八糟的涂幾張畫有什么用，還不如彈棉花來得實(shí)在，直接掙錢了。如果再這樣放任周小羽亂涂亂畫下去，指不定什么時候再讓李老師上門來，而且還指不定弄出什么大事來。

式中，

根據(jù)保角映射函數(shù)z=w(η)，在η平面中，式（4）可進(jìn)一步表示為：

2.2 散射波場

平面Rayleigh波滿足波動控制方程：

式中：φ、ψ為P 波和SV 波的標(biāo)量位移勢函數(shù)；cp為介質(zhì)P 波波速，cp=；cs為介質(zhì)剪切波速，cs=；l、m為Lamé常數(shù)。

引入復(fù)變量z=x+iy和共軛復(fù)變量=x-iy，式（7）在復(fù)平面上的穩(wěn)態(tài)波動控制方程為：

Rayleigh 波在半空間邊界處的散射解答較為復(fù)雜，不像求解SH 波散射問題時可以直接構(gòu)造出滿足水平邊界應(yīng)力自由條件的波場表達(dá)式。文中采用大圓弧假定法［18-19］處理邊界問題，即以一個半徑很大的圓（?。﹣肀平肟臻g水平邊界，并給出波函數(shù)的一般形式解，如圖3 所示。曲線Г表示大圓弧邊界，O2代表大圓弧的圓心，l表示大圓弧的半徑。

圖3 大圓弧假設(shè)示意圖Fig.3 Schematic diagram of large arc assumption

當(dāng)入射波傳播至波阻板和地基界面處，會發(fā)生波的散射。根據(jù)穩(wěn)態(tài)波動控制方程式（8）得到波函數(shù)的通解，Rayleigh 波在波阻板界面處的散射波勢函數(shù)，邊界Г 處散射波勢函數(shù)，波阻板內(nèi)部的散射波勢函數(shù)表達(dá)式如下：

2.3 邊值問題

在極坐標(biāo)系（r，θ）下，徑向位移和環(huán)向位移勢函數(shù)表達(dá)式為：

在極坐標(biāo)系（r，θ）下，正應(yīng)力分量和切應(yīng)力分量表達(dá)式如下：

在波阻板表面處滿足應(yīng)力和位移連續(xù)條件，在半空間表面滿足應(yīng)力自由條件，所以，邊值問題可以表述為：

將式（6）、式（9）～式（14）代入式（17）中，化簡整理可得如下方程組：

式中：i=j=（1，2，3，4，5，6）；Xn=（An，Bn，Cn，Dn，En，F(xiàn)n）。

在式（18）兩邊同乘以e-imθ，并以(-π,π)為積分區(qū)間對θ積分，即可得到關(guān)于待定復(fù)系數(shù)理論解的無窮線性方程組，如式（19）所示，

式中m可取(0,±1,±2,···)，對上述方程組采用截?cái)嘤邢揄?xiàng)的方法求解，確定波函數(shù)方程中待定復(fù)系數(shù)的值。

3 驗(yàn)證分析

為驗(yàn)證文中解答的正確性，考慮取波阻板密度ρB=ρS=1 950 kg/m3，彈性模量EB=ES=35 MPa，剪切模量μB=μS=13.3 MPa。即將文中模型退化為無波阻板的半空間自由場問題，并將計(jì)算結(jié)果與Rayleigh 波入射半空間自由表面的動力響應(yīng)解答進(jìn)行對比［21］。如圖4所示，文中的位移幅值與自由場位移幅值吻合，說明文中計(jì)算方法的正確性和有效性。

圖4 地表位移幅值對比Fig.4 The contrast of displacement for surface amplitude

在無窮線性方程組求解待定系數(shù)時需要采通過取有限項(xiàng)的方法實(shí)現(xiàn)，圖5 給出了n取不同項(xiàng)時的計(jì)算結(jié)果，當(dāng)截?cái)囗?xiàng)數(shù)取n=10時滿足計(jì)算精度要求，因此文中取n=10進(jìn)行計(jì)算。

圖5 不同n值時地表豎向位移變化曲線Fig.5 The variation curves of surface vertical displacement with different n

在參數(shù)分析時，首先引入位移振幅衰減系數(shù)Ar評價(jià)隔振屏障的隔振效果：

式中：U表示有波阻板時某點(diǎn)入射波場、反射波場、散射波場的位移之和；U0表示半空間入射波場的位移。

波阻板的寬度對隔振效果的影響較小，參數(shù)分析時，忽略波阻板寬度（取波阻板寬度為1）對隔振效果的影響，將入射波頻率進(jìn)行無量綱化處理為：

另外將波阻板的幾何參數(shù)歸一化到λR（Rayleigh 波波長），歸一化后的波阻板埋深為，波阻板長度為。

圖6 不同埋深時Ar變化曲線Fig.6 The variation curves of Ar with different depths

為了分析波阻板剪切模量比對隔振效果的影響，選取ηR=0.26，ρB/ρS=1.2，H/λr=0.75。不同剪切模量比的振幅衰減系數(shù)曲線如圖7所示。當(dāng)波阻板剪切模量與剪切模量相等時，波阻板后側(cè)依然能達(dá)到20%以上的隔振效果，并且隨著剪切模量增大，隔振效果明顯增加。因此，提高剪切模量比是提高隔振效果的有效措施。

圖7 不同剪切模量比時Ar變化曲線Fig.7 The variation curves of Ar with different shear modulus ratio

圖8 不同入射頻率時Ar變化曲線Fig.8 The variation curves of Ar with different incident frequency

圖9 不同長度時Ar變化曲線Fig.9 The variation curves of Ar with different length

4 結(jié)論

在二維彈性地基中有限深度處埋設(shè)矩形波阻板，采用復(fù)變函數(shù)理論和波函數(shù)展開法得到波阻板對平面Rayleigh 波散射的理論解答，并引入振幅衰減系數(shù)分析了波阻板埋置深度、模量比、彈性波頻率和波阻板長度等對隔振效果的影響。結(jié)果表明：

（1）隨著埋置深度增加，波阻板隔振效果逐漸降低，尤其是深度大于1 倍Rayleigh 波長時，隔振程度很小，因此，波阻板埋設(shè)深度不宜過大。

（2）波阻板的隔振效果隨著剪切模量比的增大而增大，在減隔振工程中，在控制成本的情況下，可以通過提高剪切模量實(shí)現(xiàn)更好的隔振效果。

（3）波阻板對低頻波的隔振效果更具優(yōu)勢。

（4）波阻板長度取值一定范圍內(nèi)，隨著長度的增大，波阻板隔振效果越好