郄祿文 許璐璐 陶佳驥

(河北大學(xué)建筑工程學(xué)院，河北 保定 071002)

地震作用下粘彈性邊界邊坡模型處理方法★

郄祿文 許璐璐 陶佳驥

(河北大學(xué)建筑工程學(xué)院，河北 保定 071002)

針對(duì)兩側(cè)高度存在差異的巖土邊坡問題，采取Plane82標(biāo)準(zhǔn)單元體來模擬巖體邊坡，并在所截取矩形區(qū)域外圍施加粘彈性邊界，隨后在對(duì)稱模型的計(jì)算條件下，截取地震記錄包含峰值在內(nèi)的前10 s歷程，將應(yīng)力時(shí)程作為荷載輸入，分析地震波在巖石邊坡中的傳播。結(jié)果顯示，針對(duì)擁有粘彈性邊界的邊坡穩(wěn)定性提出的模型處理方法簡單易行且計(jì)算結(jié)果合理可靠。

地震作用，邊坡，對(duì)稱性，粘彈性邊界

在固定邊界條件下，當(dāng)?shù)卣鸩ㄓ谶吰聝蓚?cè)傳入、傳出時(shí)，在介質(zhì)邊界上將出現(xiàn)反射波，在地震波輸入界面上會(huì)阻礙地震波的傳播作用，在地震波輸出界面上地震反射波將在邊坡內(nèi)部出現(xiàn)慣性力，將對(duì)穩(wěn)定分析的最終結(jié)果產(chǎn)生持續(xù)的影響。而在粘彈性界面條件下，由于巖石邊坡模型穩(wěn)定性的低頻特性，會(huì)對(duì)巖石邊坡的穩(wěn)定性分析結(jié)果產(chǎn)生明顯影響。LYSMER等[1]在粘彈性界面配置了粘滯阻尼器，考慮了介質(zhì)的蓄能作用卻忽視了其外力取消后的恢復(fù)作用。DEEKS等[2]建議利用粘彈性邊界針對(duì)所有頻率的巖土邊坡都具有較好的穩(wěn)定性，模型計(jì)算分析將會(huì)得到較為精確的結(jié)果。

在工程的應(yīng)用中，粘彈性邊界對(duì)于形狀規(guī)則如矩形、弧形的地基有著較準(zhǔn)確的計(jì)算結(jié)果，而對(duì)于擁有不對(duì)稱邊界的邊坡模型卻將得到錯(cuò)誤的結(jié)果。因此，本文針對(duì)邊坡模型提出一種簡單易行且計(jì)算結(jié)果可靠的計(jì)算方式。

1 粘彈性邊界單元參數(shù)選取

當(dāng)垂直人工邊界底邊界入射地震波后，經(jīng)過一段時(shí)間傳播后到達(dá)地表，會(huì)以反射波的形式在地球表面或巖體底界面處再次返射回巖土地基，并引起地表結(jié)構(gòu)的振動(dòng)，該振動(dòng)也會(huì)反射回地基中進(jìn)行傳播。在這種情況下，地基中3種形式的波同時(shí)存在，即自由場入射波、反射波以及次生場散射波[3]。

為模擬地震波在實(shí)際巖石邊坡中的波動(dòng)傳播規(guī)律，要保證波在傳入、傳出人工邊界時(shí)與原連續(xù)介質(zhì)擁有一致的傳播特性，不存在反射效應(yīng)，而被人工邊界全部吸收，地震波在巖石邊坡中發(fā)生完全透射。為消除地震波的反射，將彈簧與阻尼器設(shè)置在人工邊界上。在該邊界完全吸收散射波的前提下，這種構(gòu)造的粘彈性人工邊界的彈簧剛度與阻尼系數(shù)，可基于柱面或球面的散射波場來確定[4]。

在無限介質(zhì)中模擬地震波的傳播，為了消除邊界上散射波的反射，并保證波的傳播方式在人工邊界處和不存在人為截取邊界的真實(shí)情況非常接近?；谡硰椥苑ㄏ颉⑶邢蛉斯み吔鐟?yīng)力，假設(shè)在某半徑處截?cái)嘟橘|(zhì)，在粘彈性邊界切向和法向分別設(shè)置線性彈簧和粘性阻尼器。

2 地震荷載的輸入方式

對(duì)于地震等來自無窮遠(yuǎn)處的外源波動(dòng)，采取波場分離法。該方法假設(shè)地震產(chǎn)生的自由波場與連續(xù)介質(zhì)產(chǎn)生的散射波場共同疊加形成邊坡模型兩側(cè)邊界的振動(dòng)波場，而邊坡模型底邊界的波場則是由入射波場與散射波場疊加而成，并將地震波輸入問題轉(zhuǎn)化為波源問題。這里作用在人工邊界上的地震荷載將采用等效荷載法計(jì)算[5]。

假設(shè)，對(duì)于來自任意方向的入射波u0(x,y,t)，人工施加在邊界上的位移與應(yīng)力分別為u(xB,yB,t),σ(xB,yB,t)，其在模型人工邊界上產(chǎn)生的位移為u0(xB,yB,t)，應(yīng)滿足：

u(xB,yB,t)=u0(xB,yB,t)

(1)

σ(xB,yB,t)=σ0(xB,yB,t)

(2)

現(xiàn)將應(yīng)力fB(t)作為地震荷載于人工邊界輸入，分離彈簧—阻尼器與邊界的聯(lián)系，可得到人工邊界上點(diǎn)的應(yīng)力關(guān)系：

fB(t)=σ(xB,yB,t)+fB′(t)

(3)

代入彈簧和粘性阻尼的物理關(guān)系可得：

fB(t)=σ0(xB,yB,t)+CBu0′(xB,yB,t)+KBu0(xB,yB,t)

(4)

上式表明，作業(yè)在粘彈性邊界上的荷載分解成兩部分：在邊界處自由場引發(fā)的振動(dòng)力和邊界克服彈簧變形與阻尼力。因此，設(shè)置了粘彈性邊界條件下的位移時(shí)程更符合實(shí)際情況。

3 算例

如圖1所示的計(jì)算模型，由于兩側(cè)面積不相等，模型兩側(cè)所受的合力不等，在0.29 s以后高低兩側(cè)面的合力出現(xiàn)偏差，這種差異隨著地震波向上傳播并在頂面發(fā)生反射后逐漸增加，最終模型的計(jì)算結(jié)果偏向一側(cè)[5]。如圖2中A為低側(cè)面所受的合力變化，B為高側(cè)面所受的合力變化。

因此，在計(jì)算中采用形式上對(duì)稱的計(jì)算模型，在相同的地震荷載作用下得到兩側(cè)立面的合力變化。對(duì)于如圖1所示邊坡模型，選擇Plane82標(biāo)準(zhǔn)單元來模擬巖石邊坡，其單元長度尺寸為5 m，阻尼比系數(shù)為0.05。通過模態(tài)計(jì)算分析，巖石邊坡第一階和第二階振動(dòng)頻率分別是0.112 2和0.091。將粘彈性人工邊界設(shè)置在所截取的矩形區(qū)域邊緣，人工邊界參數(shù)列于表1。對(duì)于DP型的巖石邊坡模型，其主體材料為材料1，邊坡上部覆蓋巖塊為材料2，兩種材料的參數(shù)值見表2。

表1 人工邊界參數(shù)表

表2 巖體材料參數(shù)值

地震過程的全程加速度時(shí)程曲線如圖3所示。地震記錄為包含峰值在內(nèi)的前10 s歷程，以0.01 s為時(shí)間間隔，共采集1 001個(gè)加速度記錄。邊坡模型設(shè)置的粘彈性人工邊界屬于應(yīng)力邊界條件，在荷載輸入之前，對(duì)加速度記錄中各點(diǎn)進(jìn)行拉格朗日二次插值，并對(duì)插值公式進(jìn)行2次積分，將地震過程中的加速度記錄需要轉(zhuǎn)變?yōu)橄鄳?yīng)速度與位移時(shí)程，如圖4所示。

于地基底邊界與兩側(cè)邊界輸入應(yīng)力時(shí)程模擬地震波作用在巖石邊坡中的傳播，并選取危巖體與地表上兩點(diǎn)，在前文截取的10 s地震波作用下分析兩測點(diǎn)的位移時(shí)程，最后兩點(diǎn)位移時(shí)程與邊坡兩側(cè)面合力時(shí)程分析結(jié)果如圖5，圖6所示。

可見，危巖體頂點(diǎn)的最大位移約為11.0 cm，坡體對(duì)地震波具有明顯的放大作用。值得注意的是，最大地表位移絕對(duì)值將比底部地震入射波幅值大一倍以上(見圖5)，其主要是入射波與反射波在自由邊界處共同疊加的影響。由于這種入射波在地表的放大作用，因此在對(duì)巖石邊坡模型進(jìn)行時(shí)程分析時(shí)，應(yīng)將作為響應(yīng)輸入的地震波減小至原來的1/2，獲得的模型地表震動(dòng)波形和地震產(chǎn)生的波形大致相同。結(jié)果表明，在粘彈性人工邊界前提下，采取地震作用以應(yīng)力時(shí)程形式輸入模型，模擬結(jié)果與真實(shí)的地表運(yùn)動(dòng)基本一致。并且由合力時(shí)程曲線可知，兩側(cè)立面在運(yùn)動(dòng)過程中受力大致相等。

4 結(jié)語

本文針對(duì)兩側(cè)立面高度存在差異的巖石邊坡模型，采取Plane82標(biāo)準(zhǔn)單元來模擬巖體邊坡，并在所截取矩形區(qū)域外圍施加粘彈性邊界，隨后在對(duì)稱模型的計(jì)算條件下，截取地震記錄包含峰值在內(nèi)的前10 s歷程，將應(yīng)力時(shí)程作為荷載于兩側(cè)邊界與地基底邊界輸入，數(shù)值模擬地震波在巖體中的傳播。最后結(jié)合力的時(shí)程曲線可知，在對(duì)稱模型的計(jì)算條件下，兩側(cè)立面合力時(shí)程趨于一致，在運(yùn)動(dòng)過程中模型整體受力平衡，只發(fā)生很小的偏移。因此，對(duì)于立面高度不等的邊坡模型，采用對(duì)稱的計(jì)算方式是合理的。

[1] LYSMER J,KULEMYER R L.Finite dynamic model for infinite media[J].Journal of the Engineering Mechanics Division,1969,95(4):859-878.

[2] DEEKS,ANDREW J,RANDOLPH,et al.Axisymmetric time-domain transmitting boundaries[J].Journal of the Engineering Mechanics Division,1994,1(25):733-739.

[3] 杜修力.工程波動(dòng)理論與方法[M].北京：科學(xué)出版社,2009.

[4] 黃連娣.高混凝土拱壩的數(shù)值模擬及可靠度研究[D].大連:大連理工大學(xué),2011.

[5] 谷 音.結(jié)構(gòu)—地基動(dòng)力相互作用問題高效數(shù)值方法研究及工程應(yīng)用[D].北京:清華大學(xué),2005.

[6] 蔣新新.復(fù)雜地基條件下粘彈性人工邊界模型及其應(yīng)用[D].大連:大連理工大學(xué),2014.

Approachforslopemodelwithviscoelasticboundaryunderearthquakeaction★

QieLuwenXuLuluTaoJiaji

(CollegeofCivilEngineeringandArchitecture,HebeiUniversity,Baoding071002,China)

For a model with different heights on both sides, applying a viscoelastic boundary to the boundary， rock mass with Plane82 unit was simulated. Under the condition of the symmetric model, seismic records in the bottom boundary and both sides were inputted to simulate the earthquake effect. By analyzing the force-time curve, it was showed that the model is mechanically balanced during motion. Therefore, for the asymmetric slope model, the symmetrical calculation method is more reasonable.

earthquake action, slope, asymmetry, viscoelastic boundary

1009-6825(2017)29-0069-02

2017-08-04 ★:國家自然科學(xué)基金項(xiàng)目(61374184)資助

郄祿文(1966- )，男，博士，教授

TU457

A