Rayleigh波場的數值模擬及其應用1

白建方 馬立龍

Rayleigh波場的數值模擬及其應用1

白建方 馬立龍

（石家莊鐵道大學，土木工程學院，石家莊 050043）

本文總結了Rayleigh波場數值模擬的4種思路和各自特點，根據Lamb問題的理論分析成果提出了基于地表激振的Rayleigh波場數值模擬技術，論述了地表集中震源作用下引起地表波動場的特點及主要影響因素。在此基礎上，利用Plaxis 2D有限元軟件實現了Rayleigh波場的數值模擬，并結合算例驗證了所提方法的可行性和結果的合理性。最后，針對多層建筑結構，研究了不同地震動輸入模式下結構動力反應的特點。結果表明，Rayleigh波作用下結構的動力反應特性明顯區別于在底部輸入剪切波時的結果，不同地震動輸入模式對結構的振動形態和破壞模式有著明顯的影響。

Rayleigh波 Lamb問題 Plaxis 2D軟件 數值模擬 地震動輸入 多層結構

引言

地震發生時，震源釋放的巨大能量以體波的形式向外傳播。隨著傳播距離的增加，由于受輻射阻尼和材料阻尼的影響，體波在距震中較遠處已衰減得較弱，此時地表附近的建筑物或構筑物主要受地震面波的影響，其中又以瑞利（Rayleigh）面波為主要成分。Rayleigh波并非1種新的波型，它是由2類不同的體波（縱波和橫波）傳播到地表后，在一定條件下發生反射而產生的非均勻縱波和非均勻橫波相互干涉和疊加而成。其能量分布一般僅限于距離半空間自由表面1.5—2倍Rayleigh波波長范圍的薄層內（薩瓦林斯基，1981；楊桂通等，1988）。在傳播時，其質點在波的傳播方向與表面層法向組成的平面內作逆進的橢圓運動，波速與介質的物理特性有關，一般接近但稍小于橫波波速（約為橫波波速的0.92倍），且振幅大，傳播距離遠（吳世明，1997）。已有的震害調查表明，對于淺源地震，在距震中約0.6—5倍震源深度，或20—50km以外的區域，將能觀察到明顯的面波震害特征（崔杰等，2008）。因此，對于絕大多數的工程結構，尤其位于遠場的建筑結構，實際上多處在Rayleigh波起主要作用的地帶，故其動力響應分析應該考慮Rayleigh波的影響。

Rayleigh波自1887年被發現以來，已經在地質勘查、隧道施工超前預報等眾多領域得到了廣泛應用（佘德平，2004；王朝令等，2014；Yu等，2015）。但對Rayleigh波作用下建筑結構的動力反應特性方面的研究還相對較少，且一般多限于對特殊形狀的地下結構或特殊地形所做的理論分析，如梁建文等（2009）和巴振寧等（2014）采用直接剛度法等理論工具分析了Rayleigh波入射情況的各類特殊場地地形、地貌情況下的動力特征。對于實際地震作用下土體進入塑性之后，或者對于復雜結構的地震反應以及土-結構動力相互作用的問題，理論方法很難得到解析解。因此，研究如何利用現有成熟的有限元技術實現Rayleigh波場的數值模擬，對于推動Rayleigh波作用下工程結構的抗震性能研究具有現實意義。

本文首先總結了目前進行Rayleigh波場數值模擬的4種思路，其次根據Lamb問題的理論研究成果提出了基于地表激振的Rayleigh波場數值模擬技術，結合Plaxis 2D軟件，通過算例驗證了所提方法的可行性以及結果的合理性。最后，針對某多層建筑結構，采用所提方法分析了Rayleigh波作用下結構的動力反應特點，并與傳統的剪切波輸入以及剪切波和Rayleigh波聯合輸入下的結果進行了對比，論述了不同地震動輸入模式對建筑結構動力反應特性的影響。

1 Rayleigh波場數值模擬的4種思路

已有的關于Rayleigh波場數值模擬技術的研究，可以歸結為4種思路：

（1）直接模擬方法。即利用有限元技術直接建立包含震源在內的大尺度場地模型，結合現有震源模型模擬地震的產生、能量的釋放以及波在大范圍場地條件下的傳播特性。只要震源參數合理、模型尺寸足夠大，傳播介質的物理和幾何參數足夠準確，則理論上該方法可以完整地再現體波由震源到場地的傳播過程以及在地表經折射疊加進而形成Rayleigh波后的傳播情況。羅超（2017）曾利用該思路實現了某河谷地形的地震波場模擬。在實際應用中，該方法面臨的主要問題是如何平衡計算效率和計算精度的問題。為了提高計算精度，需要盡可能準確地建立反應實際傳播介質幾何、物理特征的場地模型，盡管已經提出了旨在提高計算效率的算法，但對于實際工程，如此大的計算規模其耗時也是巨大的。因此，在實際應用中多采用“兩步法”，對于結構所在局部場地采用精細化有限元分析，而對于遠域則假設為均勻場地或分層場地近似計算，近域和遠域通過人工邊界進行耦合（Kontoe等，2009；劉晶波等，2018）。

（2）等效人工邊界波動輸入法。劉晶波等（2006）在分析地下結構動力反應時，首先建立了包含地下結構和周圍場地的小范圍有限元模型，模型底部和側向邊界采用粘彈性人工邊界。為得到自由Rayleigh波場，首先假定某一基準點的水平位移時程，再利用Rayleigh波在地表附近任一質點所滿足的運動方程，即水平運動和豎向運動的關系，結合傅立葉變換技術得到相應的豎向位移時程。以此Rayleigh波任意角度入射到粘彈性人工邊界上，求得人工邊界節點上的波動輸入等效節點荷載，進一步以此荷載為輸入，完成地下結構的動力反應計算。

（3）近似Rayleigh波場模擬方法。岳慶霞等（2008）在研究地下綜合管廊抗震特性時利用已有的地震記錄，將其視為水平方向的波動，仍然參照Rayleigh波在地表附近任一質點所需滿足的運動方程，即水平運動和豎向運動的關系，利用傅立葉變換得到豎直方向的波動，并根據Rayleigh波在土層內部沿豎向的衰減規律，生成不同深度處的位移場，從而得到整個區域Rayleigh波位移場，在此基礎上利用有限元方法分析位于其中的地下結構的動力反應。

（4）地表激振法。在地質勘探以及隧道施工超前地質預報領域，王朝令等（2012）、蘇傳行（2017）曾以Lamb問題為原型，研究得到自由表面豎向集中脈沖激振時產生的振動在遠離振源處將以Rayleigh波為主要成分，從中詳細探討了用有限元軟件模擬Lamb問題時的網格劃分、單元選取、人工邊界的設置等問題。證實了當遠離激振位置一定距離后地面運動將主要受Rayleigh波的影響，從而提供了1種生成Rayleigh波場的數值模擬思路。施有志等（2017a，2017b）利用該思路，采用三角脈沖激振方式，通過合理設置激振力的振幅和激振時間獲得了遠場地表的Rayleigh波場，進而分析了地下綜合管廊在此波場下的動力反應。

綜上所述，在實際應用中4種思路均能在一定程度上定性地模擬Rayleigh波的波場特性，其精度無明顯差別。第1種思路理論上最完備，但計算效率和精度難以兼顧，且尚無商業程序；第2、3種思路的物理本質一樣，均以Rayleigh方程為基礎，結合Rayleigh波沿水平和豎向的衰減特征近似實現整個波場的模擬，區別在于第2種思路將外域波動場的影響等效為近場有限元模型的邊界條件，第3種思路則是直接由Rayleigh波的運動學特征生成整個波動場。這2種思路均需在現有商業程序中進行二次開發，以計算Rayleigh波場或人工邊界等效節點力。相對而言，第4種思路最為簡單易用，一方面，該方法可歸結為Lamb問題的數值模擬，可利用Lamb問題的理論研究成果；另一方面，可充分利用現有商業軟件簡單高效的前后處理技術和豐富的單元庫，無需編程或二次開發。因此，本文選擇以第4種思路為基礎展開進一步的研究。

2 Lamb問題及其研究成果

Lamb（1904）研究了作用在彈性半無限體表面的1條線或1點處以及內部的1條線或1點處等4種沖擊載荷條件下彈性波傳播問題的解析解，這些問題被統稱為Lamb問題。由于該問題的解答在動力機器基礎設計方面有一定應用價值，因此學者們對不同激振形式下Lamb問題的解答進行了廣泛的研究。王貽蓀（1979，1980，1982）系統總結了該領域的研究成果，并采用圓心位移影響函數法推導出了半空間豎向集中諧和力和水平集中諧和力作用下產生的表面位移的解析表達式；劉凱欣等（2004）和劉廣裕等（2007）進一步給出了脈沖載荷作用于彈性半空間的垂直點源問題的1個代數形式精確解，并采用積分變換方法，求得了區域脈沖載荷下1個二維Lamb問題的代數形式的精確解。本文僅給出與Rayleigh波場數值模擬相關的主要研究結論。

2.1 面源激振作用下地表振幅衰減規律

對于面源激振，Woods（1968）對均質各向同性彈性體半空間表面在豎向諧和振動的圓形基礎作用下的波動場特征進行了研究，結果表明在地表處，體波以1/2衰減，而Rayleigh波以1/衰減（為地表不同位置到震源的距離），與點源激振情況下的衰減規律一致，所以在遠場處可獲得由Rayleigh波起控制作用的波動場。

2.2 遠場和近場的劃分

研究表明縱波和橫波在近場的效應不可忽略，隨著距離的增加，在遠場體波的影響很小，Rayleigh波起控制作用。對于遠、近場的劃分，基于現場監測數據和基于理論解分析所給出的公式不完全相同。

令特征因數a=s，其中代表激振頻率，代表震源距，s代表場地剪切波速，王貽蓀等（1979）根據集中諧和力作用下地表豎向位移的解析公式，認為a＜2.0時可視為近場，至少在a＞3.5后可視為遠場。

根據Rayleigh波長和震源距的關系，高廣運（1998）認為對于一般土場地（泊松比≥0.35）、＞2.5時（為Rayleigh波波長）以及對于較硬巖石場地（泊松比≤0.3）、＞5.0時，可視為遠場，小于該距離時為近場。

于文福（2017）利用Rayleigh波運動方程推導出了Rayleigh波引起的質點橢圓軌跡的橢圓極化率（橢圓長軸與短軸之比）理論公式，以此作為參照標準，通過對比點源激振方式下地表不同距離質點運動軌跡的橢圓極化率的數值結果與理論值之間的差異，得出了對于不同激振頻率的信號，在震源距小于5倍Rayleigh波波長區間屬于近場區，此時體波干擾比較明顯；而當震源距大于5倍Rayleigh波波長后，則是Rayleigh波占優的遠場區，體波干擾較小。這一結論與高廣運（1998）在較硬巖石場地下的研究結論一致。

楊先健等（2013）對于動力機器基礎的研究結果認為，當滿足＜[15，20]0時為近場波動（0為機器基礎當量半徑），＞[15，20]0時為遠場波動。

2.3 振源半徑和介質剪切波速的影響

在一定激振頻率下，振源半徑越大，面波所占比重越高，因而幾何衰減也趨緩慢。反之振源半徑越小，面波所占比重相應減小，體波所占比重增加，幾何衰減趨于加快。

當介質的剪切模量增大時，相應的剪切波速增大。在振源半徑不變時，近源面波所占比重減小，近源幾何衰減加快（楊先健等，2013）。

2.4 激振頻率的影響

按照理想彈性半空間的理論分析結果，激振頻率越大，地面豎向振幅的衰減越慢；并且隨著激振頻率（撓頻）的增大，地面豎向振幅的衰減曲線隨距離的起伏波動越明顯（即豎向位移的衰減并非單調遞減）（王貽蓀，1982），因此在波動場的數值模擬時應合理設置激振頻率。但需要指出的是，上述結論是根據理想彈性半空間的假設得到的，對于實際半空間場地，大量實測資料顯示的則是隨著撓頻增大，衰減越快，這可能是由于實際場地的非彈性性質所致（楊桂通等，1988）。

2.5 場地剪切波速的影響

對于均勻的彈性半空間，Rayleigh波波速與激振頻率無關，只與介質物理特性（泊松比和剪切波速等）有關，即在數值模擬Rayleigh波場時若采用均質彈性半空間假設，無論激振頻率多大，選定場地材料特性后其Rayleigh波速為確定值。此外，根據Rayleigh波的運動方程可知，在地表處其水平位移與豎向位移相差90°相位，因此在Rayleigh波作用下將會先感受到豎向振動，之后再感受到水平運動，這與實際地震相一致。

3 Rayleigh波場的數值模擬

采用Plaxis 2D程序（劉志祥等，2017），模擬Lamb遠場問題，并分析自由場表面質點的運動規律，以校驗利用現有有限元軟件模擬Rayleigh波場的可行性和結果的合理性。

3.1 模型參數的確定原則

根據Lamb問題的理論研究成果，選取具備一定激振面積的面源激振方式，荷載形式為簡諧力荷載。根據在場地某處擬獲得的地面運動，諧振力的振幅可由理論衰減公式初步確定，也可通過簡單試算來確定；激振頻率由擬獲取的主要頻率范圍來確定；建模時由于問題的空間軸對稱性質，采用二維軸對稱模型。設置有限元的邊界距離時，為避免波在人工邊界處反射的干擾研究區域的波動場，盡可能使波達不到邊界處，或即使達到邊界處、其反射波尚未傳播到擬研究區域，其可通過調整場地剪切波速、激振時間以及場地尺寸來實現，同時在側向邊界處施加Plaxis中的粘性邊界作為吸收邊界，以減少波的反射效應（Brinkgreve等，2016）。

根據上述原則，初步設計出模擬Rayleigh波場的有限元模型，如圖1所示。該模型采用15節點三角形單元建立軸對稱模型，場地土層水平尺寸50m、深30m，不考慮地下水的影響，共包含3539個節點。振源假設為安置在厚0.2m、直徑1.0m的混凝土基礎上的振動裝置。動力荷載通過簡諧荷載模擬，頻率10Hz、周期0.1s，振幅經試算后確定為5×104kN/m2，方向向下，此時位于遠場地表的豎向位移約0.01m。除了基礎的重量，振動裝置本身的重量簡化為8kN/m2的均布荷載。地基土為砂質黏土，假設為線彈性，參數如表1所示，其中地基土的彈性模量較高，原因為動荷載較快使得地基土的動力剛度一般大于靜力剛度。振動基礎由Plaxis 2D中的板單元模擬，其物理參數如表2所示。

根據上述參數算出該場地剪切波速為97.07m/s，Rayleigh波波速約為90m/s，波長約為9m。在遠、近場的劃分上，按照王貽蓀等（1979）的判斷標準，約為6m，可視為遠場；按照高廣運（1998）或于文福（2017）的標準，約45m，可視為遠場；按照楊先健等（2013）的經驗公式，約為7.5—10m，視為遠場。可見高光遠、于文福的標準最為嚴格，王貽蓀的標準最寬松，而楊先健的經驗公式介于兩者之間。

圖1 Rayleigh波場數值模擬的有限元網格劃分

表1 Plaxis 2D中土層材料參數

表2 Plaxis 2D中振動基礎材料參數

3.2 數值模擬結果

動力分析設為2個階段，首先簡諧荷載強迫振動0.5s（根據Rayleigh波速可知此時傳播距離約為45m，尚未達到側向人工邊界），其次自由振動0.5s。為了與上述理論研究成果對比，算例中未考慮土體阻尼的影響。為研究地表距振源不同距離質點的振動情況，分別在地表距振源2m、4m、6m、10m、15m、20m、25m、30m和40m處設置檢測點A—I（圖1）。

在設計的激振荷載下，自由振動期間不同時段的有限元網格變形情況（放大50倍后）如圖2所示。

圖2 不同時段有限元網格變形示意圖

由圖2可以看出，在面源簡諧荷載激振作用下，地表運動呈明顯的滾動特征，在波動傳播到某一位置時，該處網格將產生橢圓形扭曲，與Rayleigh波運動特性一致。

為對比遠、近場的不同劃分標準，分別給出圖1中地表不同位置處質點的運動軌跡，如圖3所示。

圖3 地表不同位置處質點運動軌跡

圖3中的曲線代表地表不同位置處質點的水平位移與豎向位移關系，即質點的運動軌跡。由圖可以看出，由于每個點距離震源的距離不同，橢圓的形態也體現出較大差別，靠近震源處的橢圓均存在較明顯的傾斜，且距離不同，橢圓的傾斜角度也不盡相同，隨著震源距的增大（大于15m），橢圓形態逐漸趨于穩定，并且最終形成1個長軸垂直于自由表面的直立橢圓。造成這種現象的原因是靠近震源處，體波干擾較強，且面波的轉換還沒有完全成型，橢圓形態不穩定，而在遠離震源的位置，由于傳播速度的不同以及衰減速度的差異，體波的影響在面波主能量到達之前已經很微弱甚至消失，所以質點的運動開始趨于理想的Rayleigh面波橢圓運動。地表I處（距離震源40m）質點的運動未完成1個橢圓周期，說明此時Rayleigh波剛傳到此處，也證明了該模型邊界距離的選取符合上述設置原則。

對比上述遠、近場的判斷標準可以看出，按照王貽蓀（1979）的標準，當＞6m視為遠場，但從數值結果看，盡管此時形成了橢圓運動，但橢圓形狀距理想Rayleigh波引起的橢圓形狀尚有明顯差異，說明此時Rayleigh波已開始占主導成分，可作為遠、近場判斷標準的下限值；按照高廣運（1998）和于文福（2017）的標準，當＞45m可視為遠場，此標準稍顯嚴格，可作為判斷條件上限值；按楊先健等（2013）的經驗公式，當＞10m可視為遠場，該經驗公式由實際勘測資料確定，計入了實際場地的材料阻尼，考慮到本算例尚無計入土體阻尼影響的因素，但此標準與數值結果符合較好，可作為Rayleigh波場數值模擬中初步判斷遠、近場的標準。在實際應用中，還應考慮計算的目的、對精度的要求以及場地模型的大小和計算效率等因素，在上、下限之間合理選擇。

為細致觀察該激振條件下地表動力反應的特點，以地表C處為例，給出其質點運動軌跡，如圖4所示。

圖4 地表C處質點運動軌跡

由圖4可以看出，在加載之初，運動雜亂無章且振幅非常小，這是由于縱波的傳播速度最快，質點先接收到縱波信息，隨后接收到橫波信息，除此之外還有縱波、橫波在地表干涉而成的首波效應，在諸多因素作用下質點運動無明顯規律。隨著加載繼續，后期質點的運動呈明顯的橢圓軌跡，并且振幅在整個振動過程中最大，表明質點在這一時期內主要受Rayleigh面波影響而發生橢圓形的運動，但由于此時距離震源位置尚不遠，體波的干擾比較明顯，使得橢圓的長軸有明顯的傾角；當Rayleigh面波到達后，質點的振動又逐步回到平衡位置。

圖5給出了地表C處水平運動和豎向運動的位移時程，從中可以看出最早接收到體波信息，隨后大部分時間受Rayleigh波控制，在Rayleigh波作用下會同時產生水平運動和豎向運動，且一般豎向位移分量大于水平位移分量，兩者之間有一定的相位差。在自由振動階段，即使沒有材料阻尼，但由于幾何阻尼的存在，波動也會逐漸衰減。

圖5 地表C處質點不同方向位移時程

4 不同地震動輸入模式對多層建筑結構動力反應的影響

土體參數、激振裝置和震源參數的設置均與上述算例相同。在距離震源14—24m處有一多層建筑結構（此時建筑結構基本位于遠場區域，地面運動以Rayleigh波為主），共兩跨，總寬度10m，地上5層，每層3m，地下1層，為2m高的地下室，地下室墻體與周圍土體間設置界面單元以模擬土-結構動力相互作用。建筑物的墻體、樓板和地下室底板用Plaxis 2D中的板單元模擬，參數如表3所示。中間立柱用點對點錨桿單元模擬，參數如表4所示。由于建筑物的存在，軸對稱假設不再適用，需要建立完整的場地模型，含有建筑物的有限元模型如圖6所示。

圖6 含建筑結構的場地有限元模型

表3 建筑物材料屬性（板單元）

表4 點對點錨桿材料屬性

4.1 Rayleigh波作用下多層建筑結構的動力反應結果

采用地表激振法實現Rayleigh波場的數值模擬，激振參數與上述算例保持一致，為了便于觀察多層結構的后期地震反應特點，將自由振動時間增加到4.5s，強迫振動時間仍為0.5s。在此激振條件下，不同時刻建筑結構的振動形態（放大50倍）如圖7所示，頂層兩側角點處的水平和豎向位移時程如圖8所示。

圖7 Rayleigh波作用下不同時刻建筑物的振動形態

圖8 Rayleigh波作用下建筑物頂層兩側角點的水平和豎向位移時程

4.2 底部剪切波作用下多層建筑結構的動力反應結果

傳統的抗震反應分析通常在有限元模型的底部輸入垂直向上傳播的剪切波，本文用簡諧荷載進行描述，基底水平運動振幅設為0.01m（近似等于Rayleigh波作用下自由場分析時建筑結構所在位置的豎向位移大小），頻率10Hz，激振時間與Rayleigh波作用情況一致，即0.5s的強迫振動和4.5s的自由振動，其它參數與上述算例相同。不同時刻建筑物的振動形態（放大200倍）如圖9所示，頂層兩側角點處的水平和豎向位移時程如圖10所示。

圖9 剪切波作用下不同時刻建筑物的振動形態

圖10 剪切波作用下建筑物頂層兩側角點的水平和豎向位移時程

4.3 剪切波和Rayleigh波共同作用下建筑結構的動力反應結果

實際地震作用時，由于不同地震波傳播速度的差異，在遠場的建筑結構首先接收到體波，考慮到縱波的衰減較快，此時主要以剪切波為主，隨后將接收到以Rayleigh波為主的面波。為模擬實際地震作用的特點，針對上述場地模型，首先考慮在底部輸入水平運動時程以模擬剪切波作用，激振時間0.5s，隨后留出0.5s自由振動時間，然后由地表激振方式產生Rayleigh波，荷載大小、強迫振動時間與上述算例相同，自由振動時間設為3.5s。在該剪切波和Rayleigh波的聯合作用下，建筑物的前半段振動形態與圖9類似，后半段則與圖7類似，頂層兩側角點處的水平和豎向位移時程如圖11所示。

4.4 數值試驗結果分析

對比圖7和圖9可以看出，Rayleigh波作用下建筑結構的振動形態與剪切波作用下存在明顯差異，在剪切波作用下，建筑結構主要表現為左右搖晃，以剪切變形為主；而Rayleigh波由于會同時引起地面的豎向運動和水平運動，所以導致建筑結構整體運動形態由左右擺動和上下晃動耦合而成。振動形態的不同進而將導致不同的破壞位置和破壞模式。

圖11 剪切波與Rayleigh波聯合作用下建筑物頂層兩側角點的水平和豎向位移時程

對比圖8、圖10和圖11可以發現，剪切波作用下結構頂層的反應以水平運動為主，豎向位移很小；而在Rayleigh波作用下，整個動力分析時段內的豎向位移基本大于水平位移；在剪切波和Rayleigh波聯合作用下，頂層位置前半段以水平運動為主，在Rayleigh波到達后豎向運動與水平運動大致在同一量級，此時建筑結構的動力反應最大。此外，對于抗側剛度不大的多層框架結構類型，無論哪種地震動輸入模式，其水平方向的動力反應持續時間均遠大于豎向反應，如在Rayleigh波作用下，激振時間設為0.5s時，頂層水平位移約在8s后才趨于0，而該層的豎向震動約在3s后就已經結束。由此可見，不同地震動輸入模式的影響主要反映在地震波的持時內，對結構自由振動期間的影響不大，原因是無論哪種激振方式，結構后期的振動都將以水平反應為主。最后，無論哪種地震動輸入模式，頂層左、右角點的水平運動均一致，而其豎向運動則有明顯差異，主要是由于樓板單元平面內的剛度遠大于其平面外的剛度。

數值試驗表明，本文所提出的地表激振Rayleigh波數值模擬技術，可用于各類復雜結構在不同地震動輸入模式下的動力特性研究。

5 結論

本文總結了Rayleigh波場數值模擬的4種思路，結合Lamb問題的研究成果，采用地表激振的方法實現了Rayleigh波場的數值模擬，結合算例驗證了該方法的可行性和結果的合理性。在此基礎上，針對某多層建筑結構對比了Rayleigh波輸入、剪切波輸入以及剪切波和Rayleigh波聯合輸入等不同地震動輸入模式對多層建筑結構動力反應的影響。

結果表明，在不同的地震動輸入方式下，多層建筑結構的振動形態和破壞模式會產生明顯的差異，因此在結構抗震分析時，有必要根據建筑場地所處位置選擇更符合實際的地震動輸入方式。基于地表激振的Rayleigh波數值模擬技術實現方便、計算效率高，同時還可充分利用現有商業軟件豐富的本構模型和強大的處理功能，為今后研究土-結構動力相互作用問題、不同地震動輸入模式、材料非線性或復雜地形條件對Rayleigh波場的影響等問題奠定了基礎。

巴振寧，梁建文，2014．瑞雷波斜入射下層狀半空間中沉積谷地周圍的三維散射研究．地震學報，36（4）：571—583.

崔杰，張為，張建國，2008．Rayleigh波對淺地表地基震害的影響．華南地震，28（2）：10—18．

高廣運，1998．非連續屏障地面隔振理論與應用．杭州：浙江大學．

梁建文，魏新磊，Lee V．W.，2009．圓弧形沉積谷地對Rayleigh波三維散射解析解．天津大學學報，42（1）：24—34．

劉廣裕，劉凱欣，2007．區域脈沖載荷下二維Lamb問題的精確求解．固體力學學報，28（4）：341—346．

劉晶波，李彬，2006．Rayleigh波作用下地下結構的動力反應分析．工程力學，23（10）：132—135，131．

劉晶波，譚輝，寶鑫等，2018．土-結構動力相互作用分析中基于人工邊界子結構的地震波動輸入方法．力學學報，50（1）：32—43．

劉凱欣，劉廣裕，2004．垂直點源問題的一個精確解．科學通報，49（5）：419—423．

劉志祥，張海清，2017．PLAXIS 2D基礎教程．北京：機械工業出版社．

羅超，2017．河谷地形及土-結構相互作用對大跨度橋梁地震反應的影響．上海：同濟大學，78—94.

薩瓦林斯基，1981．地震波．段星北，譯．北京：科學出版社．

佘德平，2004．波場數值模擬技術．勘探地球物理進展，27（1）：16—21．

施有志，林樹枝，趙花麗，2017a．瑞利阻尼參數對瑞利波作用下場地動力響應的影響．人民長江，48（3）：76—80．

施有志，孫愛琴，林樹枝等，2017b．Rayleigh波作用下地下綜合管廊動力響應三維數值分析．世界地震工程，33（4）：196—210．

蘇傳行，2017．二維起伏地形瑞利面波數值模擬研究，成都：西南交通大學，18—22．

王朝令，劉爭平，2012．黏彈性邊界條件在ANSYS有限元波場模擬中的實現．大地測量與地球動力學，32（2）：28—31．

王朝令，劉爭平，黃云艷等，2014．隧道地震預報波場的有限元數值模擬．吉林大學學報（地球科學版），44（4）：1369—1381．

王貽蓀，1979．關于拉姆問題的精確解．湖南大學學報，6（1）：61—76．

王貽蓀，1980．半無限體表面在豎向集中諧和力作用下表面豎向位移的精確解．力學學報，12（4）：386—391．

王貽蓀，1982．地面波動分析若干問題．建筑結構學報，3（2）：56—67．

吳世明，1997．土介質中的波．北京：科學出版社．

楊桂通，張善元，1988．彈性動力學．北京：中國鐵道出版社．

楊先健，徐建，張翠紅，2013．土-基礎的振動與隔振．北京：中國建筑工業出版社．

于文福，2017．點源橫向非均勻介質瑞利面波波場數值模擬研究．成都：西南交通大學，66—72．

岳慶霞，李杰，2008．近似Rayleigh地震波作用下地下綜合管廊響應分析．防災減災工程學報，28（4）：409—416．

Brinkgreve R. B. J., Kumarswamy S., Swolfs, W. M., 2016. PLAXIS 2D reference manual. Delft, The Netherlands: PLAXIS BV, 286—290.

Kontoe S., Zdravkovic L., Potts D. M., 2009. An assessment of the domain reduction method as an advanced boundary condition and some pitfalls in the use of conventional absorbing boundaries. International Journal for Numerical and Analytical Methods in Geomechanics, 33(3): 309—330.

Lamb H., 1904. On the propagation of tremors over the surface of an elastic solid. Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character, 203: 1—42.

Woods R. D., 1968. Screening of surface waves in soils. Journal of Soil Mechanics and Foundations Division ASCE, 94(SM4): 951—979.

Yu W. F., Liu Z. P., 2015. A numerical study of the Rayleigh wave particle motions excited by a point source and Poisson’s ratio for lateral inhomogeneous half-spaces. Journal of Applied Geophysics, 123: 242—255.

Numerical Modeling Techniques of Rayleigh Wave Field and Its Application

Bai Jianfang and Ma Lilong

(School of Civil Engineering, Shijiazhuang Railway University, Shijiazhuang 050043, China)

Four different approaches about numerical modeling techniques of Rayleigh wave field were summarized. The techniques base on vertical harmonic point load on the free surface has been proposed in this paper. The main research achievements of Lamb issues, such as the characteristics of Rayleigh wave field under point load and its main influencing factors were discussed. On this basis the numerical modeling of Rayleigh wave field with Plaxis 2D was conducted and the effectiveness of the approach was verified with an example. The dynamic response of a multi-story building structure during propagation of Rayleigh wave was analyzed using the proposed method, and the results were compared to that from the traditional methods. Numerical analysis results show that different wave motion input methods have great effect on the dynamic response of building structures. The vibration shape and failure modes of the structure were greatly influenced by seismic input modes.

Rayleigh wave; Lamb issues; Plaxis 2D software; Numerical modeling; Wave motion input; Multi-story building

10.11899/zzfy20190207

河北省大型基礎設施防災減災協同創新中心項目，河北省重點學科建設項目（橋梁與隧道工程）

2018-07-23

白建方，男，生于1976年。講師。主要從事工程系統抗震方面的研究。E-mail：bjf2004@126.com

白建方，馬立龍，2019．Rayleigh波場的數值模擬及其應用．震災防御技術，14（2）：328—340．