張佳文，李明超，韓 帥，閆文鈺

(1.水利工程仿真與安全國家重點(diǎn)實(shí)驗(yàn)室，天津大學(xué)，天津 300350；2.香港理工大學(xué)建筑與房地產(chǎn)學(xué)系智能建造實(shí)驗(yàn)室，香港 999077)

不規(guī)則地形條件對地震波傳播路徑和地震烈度異常區(qū)的發(fā)生均有重大影響，主要體現(xiàn)在局部地形對地震波的散射和波型轉(zhuǎn)換[1-3]，國內(nèi)外震害資料也反映了其對地表位移和加速度幅值放大效應(yīng)等的巨大影響[4-6]。因此，場地反應(yīng)分析一直以來都是地震學(xué)及工程抗震領(lǐng)域關(guān)注的重點(diǎn)。此外，近年來大量臺陣記錄均證實(shí)了地震波入射角度的存在[7]，斜入射地震波對結(jié)構(gòu)物振動響應(yīng)的危害性成為共識[8-9]，地震波的入射角度與行波效應(yīng)等問題的引入使得不規(guī)則地形場地地震波場的求解難度進(jìn)一步升級。

目前，場地反應(yīng)的求解主要有解析法和數(shù)值模擬法[10]，解析方法包括波函數(shù)展開法、分離變量法、幾何射線理論等十幾種類型[11-14]，已成功應(yīng)用在了規(guī)則局部地形中，例如：均質(zhì)長方形場地[15]、半圓和半橢圓形峽谷[16]、任意圓弧凹陷場地[17]等。然而解析方法限制條件較多，針對場地介質(zhì)本構(gòu)關(guān)系復(fù)雜或形狀不規(guī)則地形條件，全域波場的解析求解難以實(shí)現(xiàn)。相比而言，數(shù)值模擬方法適用范圍更加廣泛，且能更好地應(yīng)用到抗震分析中[18]，各類基于有限元法、有限差分法和邊界元法等的研究也正在蓬勃發(fā)展[19-21]。斜入射地震波作用下不規(guī)則地形地震波場的求解作為一個較難的課題，可通過解析推導(dǎo)與數(shù)值模擬相結(jié)合來實(shí)現(xiàn)。

ASHFORD 和SITAR[22]基于廣義一致透射邊界等分析了陡坎地形在傾斜SH 波作用下位移等的放大效應(yīng)，將不規(guī)則場地分區(qū)域進(jìn)行求解，并首次引入了此類問題的時域化解法，但波場計(jì)算的精度不高且并未研究平面內(nèi)的波場；李山有和廖振鵬[23]應(yīng)用顯式有限元法對波動方程進(jìn)行離散化求解，研究了斷層臺階場地的反應(yīng)特點(diǎn)，并指出了地表運(yùn)動的差異性，有效考慮了不規(guī)則地形對地震波型的轉(zhuǎn)換問題，然而波場計(jì)算依賴于邊界條件且計(jì)算繁瑣；顧亮等[24]基于有限元數(shù)值模擬與透射邊界理論，研究了SV 波斜入射下陡坎地形地表峰值加速度的響應(yīng)問題，文中對自由波場的計(jì)算較為簡便，但對局部地形的散射問題考慮不充分，因此精度不高；趙密等[25]基于黏性邊界條件推導(dǎo)了邊界上的等效節(jié)點(diǎn)力，研究了SV 波傾斜入射至階梯型場地時的地震動輸入方法，可將地表位移的峰值誤差控制在5%以內(nèi)，但并未詳細(xì)介紹垂直入射與斜入射時地震波場計(jì)算的不同之處；BAO 等[26]提出了混合波場的思路，結(jié)合有限元法和波動輸入方法研究了SV 波垂直和傾斜入射下階梯型地形的地震動輸入方法，其精確度和效率有了較大的提升，然而地形條件和波型較為單一。以上幾類典型的方法不斷提升了計(jì)算精度與適用范圍，但是對局部地形的散射問題考慮還不夠充分，所采用的地震波波型、地震入射角度、場地地形條件等也有限，且在地震動輸入時普遍對人工邊界的限制條件較高，計(jì)算過程十分復(fù)雜；因此，亟需提出更為通用和高效的地震波場計(jì)算與輸入方法。

本文提出了一種不規(guī)則地形條件下的地震波輸入方法。依據(jù)地震波傳播規(guī)律進(jìn)行波場分離，再將非外行波場轉(zhuǎn)換為人工邊界上的等效節(jié)點(diǎn)力形式以完成計(jì)算。該方法有效考慮了局部地形中散射波場的影響，計(jì)算過程明晰簡便，易于推廣至復(fù)雜場地條件。驗(yàn)證部分選用了不同的地形條件，將計(jì)算結(jié)果分別與精確解、參考解、成熟方法解進(jìn)行對比；并應(yīng)用到多組斜入射P 波和SV 波作用下場地的振動的分析中。

1 方法原理

1.1 地震波場分離

不規(guī)則地形條件下場地的全域自由波場無法統(tǒng)一求解，需要在各個邊界上針對不同地形特征與輸入條件分開計(jì)算。李山有和廖振鵬[23]提出地震波的輸入應(yīng)滿足非外行波的條件(非外行波即通過邊界向有限域內(nèi)傳播的波和平行邊界的地震波)，并證明了地震波垂直入射情況下，底邊界的非外行波即入射波，而側(cè)邊界的非外行波可看作與遠(yuǎn)場介質(zhì)對應(yīng)的自由波場。而本實(shí)驗(yàn)進(jìn)一步考慮到傾斜入射的情況，認(rèn)為斜入射情況下(以左下角入射為例)，右邊界區(qū)域在初始階段便存在反射波場，因此非外行波場也只有入射波場。此外，杜修力和趙密[27]曾也指出，目前所有的人工邊界都只能模擬傳遞到無限域的外行波，而由結(jié)構(gòu)物存在、不規(guī)則地形條件、局部地質(zhì)等引發(fā)的散射問題，都需在邊界處進(jìn)行波場分離。

因此，根據(jù)地震波入射、反射和散射等原理，首先在不規(guī)則地形場地的各個邊界面上進(jìn)行波場分離。當(dāng)?shù)卣鸩ù怪比肷鋾r，在底面邊界上將波場分為入射波場和邊界外行場，左右邊界面上分離為自由波場和散射場；當(dāng)?shù)卣鸩▋A斜入射時(以從左下角入射為例)，在底面邊界和右面邊界將波場分為入射波場和邊界外行場，而左面界面上則分離為自由波場和散射場即可。進(jìn)而輸入非外行波場，如圖1 所示：1) 垂直入射情況下在底部邊界輸入入射波場，兩側(cè)邊界輸入各自高度對應(yīng)的自由波場；2) 傾斜入射情況下，在底部邊界和右側(cè)邊界輸入入射波場，在左邊界輸入自由波場。

圖1 不規(guī)則地形下各邊界地震波場輸入Fig.1 Input of seismic wave field at each boundary under irregular terrain

1.2 地震波場求解

入射波場和自由波場是決定場地響應(yīng)的主要因素，需要分別確定；由局部地形引起的散射波場則通過有限元計(jì)算體現(xiàn)。

1) 入射波場：地震的入射波場求解可以直接通過幾何推導(dǎo)實(shí)現(xiàn)。如圖2 所示，利用人工邊界截斷形成一定的有限域空間，假設(shè)地震P 波從左下角(x0,y0)點(diǎn)進(jìn)行入射且振動位移時程為uP(t)。由于P 波入射時質(zhì)點(diǎn)振動方向與傳播方向一致，其入射波場的水平和豎向位移計(jì)算公式如下：

圖2 斜入射P 波作用下的地震波場Fig.2 Seismic wave field under obliquely incident P waves

式中：θ1為P 波入射的角度；Δt1即為P 波傳播到邊界上任一節(jié)點(diǎn)(x,y)的時間；cP為地震P 波的波速。

假設(shè)地震SV 波從左下角(x0,y0)點(diǎn)進(jìn)行入射且振動位移時程為uSV(t)，由于SV 波入射時質(zhì)點(diǎn)振動方向與傳播方向相互垂直，其入射波場的水平和豎向位移計(jì)算如下：

式中：θ2為SV 波入射的角度；Δt2即為SV 波傳播到邊界上任一節(jié)點(diǎn)(x,y)的時間；cSV為地震SV 波的波速。

2) 自由波場：如圖3 所示，自由波場包括入射場和反射場，例如P 波斜入射時又會在界面分層處同時產(chǎn)生反射P 波和反射SV 波，且幅值系數(shù)發(fā)生了變化，因此不能通過簡單的幾何推導(dǎo)求得。

圖3 斜入射地震波傳播規(guī)律Fig.3 Propagation law of obliquely incident seismic waves

目前在研究斜入射地震波下場地的自由波場時，主要有頻域和時域兩種方法，后者更適用于復(fù)雜的場地條件且與現(xiàn)代大型計(jì)算軟件匹配良好，因而廣泛使用。目前比較流行的有一維時域化方法[28]，其結(jié)合了集中質(zhì)量有限元法和中心差分法以建立二維波動方程，計(jì)算出單一垂線上節(jié)點(diǎn)的位移，再通過幾何推導(dǎo)擴(kuò)展到全域波場。

主要過程如下：圖4 為離散化后的半空間模型，橫向和豎向網(wǎng)格尺寸為Δx和Δy，地震P 波或SV 波從模型底部的人工邊界進(jìn)行輸入。將y軸上的三類節(jié)點(diǎn)包括中間節(jié)點(diǎn)(0,n)、自由表面節(jié)點(diǎn)(0, 0)和人工邊界上的節(jié)點(diǎn)(0,N)在PΔt時刻的運(yùn)動方程表示如下：

式中：M、C、K分別為節(jié)點(diǎn)的質(zhì)量、阻尼和剛度矩陣，下標(biāo)為節(jié)點(diǎn)的坐標(biāo)信息；CB和FB分別為邊界的阻尼系數(shù)和需要輸入的集中力。在P 波和SV 波入射時分別有不同的FB計(jì)算方法，求解上述方程組，即可確定半空間內(nèi)y軸上所有節(jié)點(diǎn)在(p+1)Δt時刻的位移，再結(jié)合幾何關(guān)系推廣到全域即可計(jì)算出整個半無限空間的自由波場。

1.3 地震波動輸入

地震屬于振動中的外源問題，在進(jìn)行波場推求后還需進(jìn)行地震動輸入才能準(zhǔn)確計(jì)算場地的動力響應(yīng)，人工邊界和地震動輸入方式的選擇是計(jì)算中的兩個關(guān)鍵環(huán)節(jié)[29]。數(shù)值模擬過程中通常在無限域地基中截取感興趣的部分作為有限域，并在截斷處賦予各類人工邊界，其主要作用是吸收外行波以達(dá)到對無限地基輻射阻尼的模擬。地震動輸入則主要分為振動法和波動法，波動法最早由劉晶波和呂彥東[30]于1998 年提出，可有效克服振動法中對于場地材料、激勵方式和邊界處理上的簡化造成的誤差，可以處理地震波斜入射問題和復(fù)雜場地的輸入問題，也因其時空解耦性的特點(diǎn)而與各類大型通用軟件適配性極高，是目前模擬精度最高且最為流行的地震動輸入方法。波動法的核心原理是將地震動輸入轉(zhuǎn)換為人工邊界上的等效荷載力，需要引入黏性邊界和黏彈性邊界等人工邊界；因此需要同時計(jì)算由自由場產(chǎn)生的應(yīng)力和引入邊界后的附加應(yīng)力，針對不同網(wǎng)格尺寸和方向不同的應(yīng)力條件也需逐一判斷，計(jì)算過程十分復(fù)雜；尤其是考慮復(fù)雜場地條件時，等效荷載力計(jì)算的難度顯著增大，尤其是所要探討的不規(guī)則地形條件，傳統(tǒng)的波動方法已不具備高效處理的能力，亟需進(jìn)行優(yōu)化。

劉晶波等[31]和寶鑫等[19]對波動法進(jìn)行了改進(jìn)，從一維土柱模型入手證明了在邊界節(jié)點(diǎn)附近的所有相鄰節(jié)點(diǎn)同時施加入射波場位移或自由波場位移，再將輸入位移后計(jì)算得到的荷載力反向輸入回模型中可完成地震波輸入，進(jìn)而提出了基于人工邊界子結(jié)構(gòu)的地震波動輸入方法，如圖5所示。圖5(a)為不規(guī)則地形的有限元模型；圖5(b)中將含有邊界節(jié)點(diǎn)的所有單元稱為邊界子結(jié)構(gòu)，將地震波場位移時程輸入到子結(jié)構(gòu)中可得到邊界節(jié)點(diǎn)的反力，經(jīng)過證明此反力在數(shù)值上恰好等于需要的等效荷載力；圖5(c)表示重新在邊界節(jié)點(diǎn)上輸回此反力即可實(shí)現(xiàn)波動輸入。該方法保留了波動法的諸多優(yōu)勢并簡化了等效荷載力的計(jì)算流程，避免了復(fù)雜等效節(jié)點(diǎn)力的過程，可成功地將不規(guī)則地形場地的地震波輸入問題聚焦為地震波場求解問題。

圖5 基于人工邊界子結(jié)構(gòu)的波動輸入方法Fig.5 Wave input method based on substructure of artificial boundaries

2 算例驗(yàn)證

2.1 有限元模型

為了驗(yàn)證所提出的方法在各類不規(guī)則地形條件下的適配性，分別以階梯型、V 字河谷型和梯形河谷型這三類工程中常見的不規(guī)則地形條件為例進(jìn)行驗(yàn)證。圖6 展示了三類地形有限元模型尺寸信息，在模型的左側(cè)、右側(cè)和底面邊界分別采用了黏彈性邊界條件進(jìn)行約束。由于斜入射地震波場在不規(guī)則地形條件下的頻域解析解不能求得，可擴(kuò)大截取的有限域場地范圍，采用遠(yuǎn)置邊界的方法保證地震波不會通過人工邊界反射后再回到需要觀測的部位，并將遠(yuǎn)置邊界法得到的結(jié)果作為有限元近似精確解(簡稱參考解)。圖7 展示了上述三類場地相對應(yīng)的遠(yuǎn)置邊界有限元模型信息，邊界均無任何約束。圖6 和圖7 中的規(guī)則部分為正方體網(wǎng)格，尺寸為1 m×1 m，A點(diǎn)、B點(diǎn)、C點(diǎn)為自由表面觀測點(diǎn)，利用觀測點(diǎn)的位移信息來進(jìn)行對比驗(yàn)證。三類場地的材料參數(shù)及地震波速信息見表1。

表1 場地模型參數(shù)Table 1 Model parameters of the site

圖6 不規(guī)則地形場地的有限元模型 /mFig.6 Finite element model of irregular terrain site

以平面P 波和SV 波為例，傾斜入射時假設(shè)地震波統(tǒng)一從場地的左下角進(jìn)入，若考慮從右下角只需將左右邊界的波場分離方式互換即可，此處不再詳細(xì)展開。選用脈沖地震波的時程曲線及傅里葉譜曲線如圖8 所示。兩種場地模型參數(shù)條件下滿足一維時域化方法自由場計(jì)算精度的時間間隔分別為0.0029 s 和0.0020 s，統(tǒng)一選取0.002 s 作為地震波的輸入以及后續(xù)有限元計(jì)算的時間間隔。不規(guī)則地形條件下各個邊界輸入波場的計(jì)算方法參考第1 節(jié)；而在遠(yuǎn)置邊界法中，各邊界均采用對應(yīng)高度的自由波場即可。

圖8 入射脈沖波的位移時程和傅里葉譜Fig.8 Displacement and Fourier spectrum of the incident seismic wave

共進(jìn)行了4 組驗(yàn)證實(shí)驗(yàn)，選取了入射角度為0°和30°的平面脈沖P 波和入射角度為0°和15°的平面脈沖SV 波。其中，P 波和SV 波垂直入射計(jì)算時選用階梯型場地，P 波30°入射選用梯形河谷型場地，SV 波15°入射時選用V 字河谷型場地，可充分驗(yàn)證方法對不同條件下的適配性。波場計(jì)算過程均利用python 編程，有限元計(jì)算過程在大型通用軟件ABAQUS 中進(jìn)行。

2.2 P 波入射情況驗(yàn)證

如圖9 所示，將節(jié)點(diǎn)等效荷載力輸入到有限元模型的人工邊界后，得到了P 波垂直和傾斜入射下不規(guī)則場地的自由表面點(diǎn)A、點(diǎn)B、點(diǎn)C的橫向與縱向位移情況，圖例中的Ux和Uy分別代表觀測點(diǎn)橫向和縱向位移，為將結(jié)果顯示清楚模擬解部分的位移時間間隔擴(kuò)大了15 倍?？梢钥闯霰疚姆椒ǖ慕Y(jié)果(模擬解)與參考解吻合很好。其中，需要注意的是由于遠(yuǎn)置邊界模型尺寸和原有的不規(guī)則地形尺寸不同，在比較位移結(jié)果時需要去除前面一段時間(地震波從遠(yuǎn)置邊界到原有邊界處的時間)。此外，觀察結(jié)果發(fā)現(xiàn)：1) 觀測點(diǎn)的豎向高度越高，其位移幅值越大；2) 雖然兩個案例的場地參數(shù)不同，但由于輸入的脈沖波位移一樣，可以在位移幅值上做出比較，P 波傾斜入射相較于垂直入射時，橫向位移幅值增大，豎向位移幅值略有減小。這兩點(diǎn)與實(shí)際情況相符。

圖9 P 波入射下自由表面點(diǎn)位移Fig.9 Displacement of free surface under incident P waves

2.3 SV 波入射情況驗(yàn)證

圖10 為SV 波垂直和傾斜入射下自由表面點(diǎn)A、點(diǎn)B、點(diǎn)C的水平與豎向位移情況，結(jié)果同樣與參考解幾乎吻合，說明方法對P 波和SV 波都是通用的。結(jié)果發(fā)現(xiàn)：1) 觀測點(diǎn)的豎向高度越高，其位移幅值越大，規(guī)律同P 波相同；2) SV 波傾斜入射相較于垂直入射時，橫向位移幅值減小，豎向位移幅值增大。這兩點(diǎn)也符合實(shí)際規(guī)律。

圖10 SV 波入射下自由表面點(diǎn)位移Fig.10 Displacement of free surface under incident SV waves

2.4 結(jié)果與應(yīng)用討論

2.4.1 驗(yàn)證結(jié)果對比

表2 展示了本文方法與遠(yuǎn)置邊界法的對比，表2 中的計(jì)算時長和文件均指第2 次節(jié)點(diǎn)力輸入計(jì)算；實(shí)驗(yàn)使用的電腦型號為Intel(R) Core(TM) i7-8700 CPU @ 3.20 GHz，機(jī)帶RAM 為8 G，ABAQUS為2017 版本?？梢钥闯鍪褂帽疚姆椒稍诒ＷC高精度的情況下有效縮短計(jì)算時長并減小結(jié)果文件存儲空間，當(dāng)研究范圍擴(kuò)展到三維空間后優(yōu)勢將更加明顯。

表2 兩種方法的結(jié)果對比Table 2 Comparison of the two methods

2.4.2 應(yīng)用分析

由于本文方法涉及的理論推導(dǎo)和步驟較多，可以通過每次計(jì)算后的關(guān)鍵點(diǎn)位移和波場快照等進(jìn)行分步檢驗(yàn)，確保其合理性后再繼續(xù)進(jìn)行。主要流程分兩次輸入與計(jì)算，統(tǒng)一以P 波入射為例進(jìn)行說明。

1) 非外行波場計(jì)算及第1 次位移輸入：圖11(a)為P 波垂直入射時采用波場分離技術(shù)計(jì)算后，三個邊界上典型位置點(diǎn)的輸入位移時程(具體位置如圖5，G1為底部邊界中點(diǎn))。兩側(cè)自由表面的E1、F1點(diǎn)的橫向位移幅值分別為底部G1點(diǎn)的兩倍，這是因?yàn)樽杂刹▓霭巳肷洳▓龊头瓷洳▓觯v向位移均為0，結(jié)果是符合規(guī)律的。圖11(b)顯示P 波入射角度為30°時，左邊界E3點(diǎn)的橫、豎向位移均為F3點(diǎn)和G3點(diǎn)的兩倍，也是由于自由波場同時包括了入射場和反射場造成的位移幅值加倍；F3點(diǎn)的位移曲線同底部G3點(diǎn)變化趨勢一樣，但稍有延遲，這是由行波效應(yīng)引起的，整體符合波場傳播的規(guī)律。

圖11 場地關(guān)鍵點(diǎn)的位移輸入情況Fig.11 Input displacement of key points in the site

2) 第2 次等效荷載力輸入：即第1 次輸入波場位移后進(jìn)行有限元計(jì)算，提出各邊界上的等效節(jié)點(diǎn)力再反向輸回的過程。圖12 和圖13 為第2 次地震波輸入后的波場快照，可以看出：開始階段波場基本沿平面形狀進(jìn)行入射，證明各個邊界等效節(jié)點(diǎn)力計(jì)算及輸入的相對時間關(guān)系準(zhǔn)確；在傳播到不規(guī)則地表后波陣面形狀會發(fā)生變化，可看出局部地形的影響；最后隨著時間增加波場位移逐漸消失，說明黏彈性邊界有效吸收了散射波。因此，不需要將邊界進(jìn)行遠(yuǎn)置即可獲得較準(zhǔn)確的結(jié)果。

圖12 P 波垂直入射時的波場快照Fig.12 Wave field snapshot of site under vertical incident P waves

圖13 P 波入射角度30°時的波場快照Fig.13 Wave field snapshot of site under P waves with the incident angle of 30°

2.5 方法對比驗(yàn)證與拓展

2.5.1 與頻域精確解的對比

由于不規(guī)則地形的地震波場無頻域精確解，故前文結(jié)果均與參考解進(jìn)行對比。實(shí)際上，本方法針對大部分半空間場地是通用的，為進(jìn)一步證明其正確性與精度，建立規(guī)則地形模型進(jìn)行驗(yàn)證，規(guī)則地形地震波場的建立基于較多的假設(shè)與簡化，因而有精確的解析解。

參考了文獻(xiàn)[12]選取相應(yīng)的解析方法，圖14為規(guī)則地形場地的有限元模型，網(wǎng)格劃分方式、場地模型參數(shù)、邊界約束、荷載施加方式和地震波形式均與2.1 節(jié)的模型一致。選取自由表面中點(diǎn)C點(diǎn)為觀測點(diǎn)，圖15 為P 波30°入射和SV 波垂直入射下C點(diǎn)的位移?？梢钥闯霰疚姆椒ㄋ贸龅慕Y(jié)果與精確解幾乎吻合，驗(yàn)證了所提出方法的可靠性。

圖14 規(guī)則地形場地的有限元模型 /mFig.14 Finite element model of regular terrain site

圖15 地震波入射下規(guī)則場地自由表面點(diǎn)位移Fig.15 Displacement of free surface of regular site under seismic waves

2.5.2 與已有方法的對比

不規(guī)則地形場地的波場求解已有一些成熟的方法，選取文獻(xiàn)[25]中的方法進(jìn)行進(jìn)一步對比。文獻(xiàn)中將波場在人工邊界處進(jìn)行分解，用黏性邊界模擬輸入波場外的其余波場，并根據(jù)地基參數(shù)確定邊界阻尼參數(shù)等。將此方法應(yīng)用到圖6 中的階梯型場地與梯形河谷型場地模型中，分別選取了自由表面C1和C3點(diǎn)為觀測點(diǎn)，圖16 展示了其在P 波30°入射和SV 波垂直入射下的位移。

圖16 地震波入射下階梯型場地自由表面點(diǎn)位移Fig.16 Displacement of free surface of stepped site under seismic waves

從圖16 中可以看出，本方法和文獻(xiàn)[25]中方法得出的位移結(jié)果基本一致，證實(shí)了其在垂直入射和斜入射時波場計(jì)算均較為高效，同時本文方法相比于成熟方法而言計(jì)算節(jié)點(diǎn)力過程更加簡單，對人工邊界的限制也較低。

2.5.3 拓展應(yīng)用

波場分離方法同樣適用復(fù)雜的場地條件，只是當(dāng)介質(zhì)參數(shù)較為復(fù)雜時，入射波場和自由波場的求解難度會顯著增大。為使方法推廣性更強(qiáng)，這里以最為常見的成層場地為例，補(bǔ)充一下其計(jì)算思路。

成層介質(zhì)自由波場的推求已有較多頻域和時域的成熟方法[32-33]。而大部分文獻(xiàn)中沒有詳細(xì)介紹成層地基內(nèi)入射波場求解，這里推薦2 種方案：1) 頻域法求解可結(jié)合傳遞函數(shù)法和剛度矩陣而推求；2)時域方法求解可根據(jù)幾何關(guān)系[34-36]，截取地震波還未入射到自由表面前的位移時程，即可從自由波場中去除反射波場的影響，得到入射波場。

此外，更為復(fù)雜的三維真實(shí)地形條件下地震波場的構(gòu)建也是后續(xù)研究的重點(diǎn)。

3 響應(yīng)結(jié)果分析

3.1 入射角度對結(jié)果的影響

以圖5(a)中的階梯型模型為例，分別計(jì)算P 波和SV 波不同入射角度下自由表面A1點(diǎn)的位移響應(yīng)情況。因?yàn)槠矫鍿V 波入射時存在臨界入射角，一般在35°左右，因此選取SV 波入射角度范圍為0°～30°，P 波入射角度為0°～90°。如圖17 和圖18所示，隨著P 波入射角度的增大，橫向位移幅值先增大后減小，大約在60°左右時達(dá)到最大；豎向位移呈現(xiàn)逐漸減小的趨勢。但當(dāng)P 波入射角度為90°時，位移呈現(xiàn)波動變化趨勢并未完全收斂，說明本文方法對于超大角度的波動輸入問題仍有一定問題需要處理。隨著SV 波入射角度的增大，橫向位移呈現(xiàn)逐漸減小的趨勢但變化幅度較小，縱向位移幅值呈現(xiàn)逐漸增大的趨勢。對于不規(guī)則地形來說，入射角度對振動響應(yīng)的影響也是符合基本變化規(guī)律，且不同角度情況下位移變化幅度較大，說明對于地震波入射角度的考慮是必要的。

圖17 P 波不同入射角度下階梯型場地A1 點(diǎn)位移Fig.17 Displacement of A1 on the stepped site under P waves with different incident angles

圖18 SV 波不同入射角度下階梯型場地A1 點(diǎn)位移Fig.18 Displacement of A1 on the stepped site under SV waves with different incident angles

3.2 局部地形對結(jié)果的影響

選用圖14 中的規(guī)則地形場地作為參考，即100 m×50 m 的長方形場地，其動力參數(shù)均取為相同。表3 和表4 分別展示了P 波和SV 波入射條件下兩種地形條件下自由表面A1點(diǎn)的位移幅值對比；其中，相差百分比的計(jì)算方法統(tǒng)一為(規(guī)則地形下位移-階梯型地形下位移)/階梯型地形下位移。由表3 和表4 可見，在兩種場地條件下，位移隨著地震波入射角度的變化均呈現(xiàn)相同的變化規(guī)律，與3.1 節(jié)的規(guī)律一致。同時可以看出局部地形條件對場地位移的影響較大，在P 波入射條件下大部分情況階梯型場地下位移較大，SV 波入射下則相反。說明地形對地震波散射等傳播規(guī)律的影響復(fù)雜多變，會隨著波型和入射角度的變化而變化，在地震計(jì)算中對于地形因素的考慮也是必要的。

表3 P 波入射時不同地形條件下自由表面A1 點(diǎn)位移Table 3 Displacement of A1 under different terrain conditions when P waves incident

表4 SV 波入射時不同地形條件下自由表面A1 點(diǎn)位移Table 4 Displacement of A1 under different terrain conditions when SV waves incident

4 結(jié)論

合理的地震動輸入是動力分析的必要基礎(chǔ)，采取了解析推導(dǎo)結(jié)合數(shù)值模擬的手段，提出了一種基于波場分離技術(shù)和應(yīng)力型人工邊界的地震波場計(jì)算與輸入方法，分別詳細(xì)推導(dǎo)了地震P 波和SV 波在垂直入射和傾斜入射下的地震波場，在充分考慮到局部地形產(chǎn)生的散射波場的同時，采用了改進(jìn)的波動輸入方法實(shí)現(xiàn)地震波的高效輸入。此外，對比分析了多組不同入射角度下的規(guī)則和不規(guī)則地形的場地反應(yīng)情況，主要得到了以下結(jié)論：

(1) 通過與遠(yuǎn)置邊界法得到的自由表面位移的參考解進(jìn)行對比，證明了方法的精確度與高效性，且適用于多類不規(guī)則地形條件的場地以及不同的地震波型和入射角度。

(2) 所采用的改進(jìn)的波動輸入法，極大簡化了等效節(jié)點(diǎn)力計(jì)算，將地震波輸入的重難點(diǎn)聚焦于場地自由波場和入射波場的計(jì)算，使得所提出的方法易于推廣至三維場地以及成層地基情況。

(3) 地震波的入射角度變化對不規(guī)則場地地表位移響應(yīng)影響較大，自由表面點(diǎn)隨著P 波入射角度的增大，橫向位移幅值先增大后減小，豎向位移呈現(xiàn)逐漸減小的趨勢；隨著SV 波入射角度增大，橫向位移呈現(xiàn)逐漸減小的趨勢，縱向位移規(guī)律相反。整體變換規(guī)律與規(guī)則場地條件下的相一致且符合實(shí)際。

(4) 不規(guī)則地形條件對場地的地表位移響應(yīng)也有較大影響，地震波通過局部地形時傳播路徑會發(fā)生變化并產(chǎn)生額外的散射波場，因此，其振動反應(yīng)相較于規(guī)則場地有明顯不同，并同時受到地震波型和入射角度等的影響。

所提出的方法和得出的規(guī)律性結(jié)論可為不規(guī)則地形上的振動響應(yīng)分析提供手段，后續(xù)研究將進(jìn)一步聚焦于三維真實(shí)地形條件下地震波場的構(gòu)建與輸入方法改進(jìn)提升。