廖志華 梁騰龍

(廣西大學(xué)土木建筑工程學(xué)院 廣西 南寧 530001)

一、緒論

利用有限單元法進(jìn)行地震動力分析時，所截取有限元模型尺寸比實際場地尺寸要小得多，地震波在有限元計算區(qū)域傳播時，波動很快就傳播到計算區(qū)域邊界而發(fā)生反射，影響計算結(jié)果。若要保證在計算時間內(nèi)，波動傳播不到計算區(qū)域邊界上，需要截取較大計算區(qū)域進(jìn)行計算，處理辦法為采用遠(yuǎn)置邊界。所謂遠(yuǎn)置邊界就是將人工邊界取與廣義結(jié)構(gòu)距離足夠遠(yuǎn)，其人工邊界采用固定處理，上述處理避開了反射波動對監(jiān)測點或廣義結(jié)構(gòu)的二次作用，保證結(jié)果精確。當(dāng)?shù)卣饎恿Ψ治鲇嬎銜r間較長，計算區(qū)域介質(zhì)波速較大時，采用遠(yuǎn)置邊界需截取較大計算區(qū)域。有限元計算是將有限計算區(qū)域劃分為有限個單元進(jìn)行分析，在單元大小確定的情況下，若計算區(qū)域過大，就會增加計算自由度而大量地消耗計算機資源。

建立合理的人工邊界來模擬有限計算區(qū)域截斷前的狀態(tài)，同時能夠合理模擬波動在截斷介質(zhì)中傳播是十分有必要的[1]。動力學(xué)一般采用人工邊界來模擬土體輻射阻尼，使得波動在截取有限計算區(qū)域內(nèi)順利透過人工邊界而不產(chǎn)生反射，邊界應(yīng)力、位移在介質(zhì)截取邊界處與原來相同，以達(dá)到還原地震場的情況，保證模擬的精度。

粘彈性邊界由Deeks于1994年基于柱面波動理論的基礎(chǔ)上推導(dǎo)，彌補粘性邊界存在的問題，保證截取計算區(qū)域邊界既能模擬波動穿過邊界被吸收的情況，同時模擬地基彈性恢復(fù)力情況。其基本思路是在截取有限計算區(qū)域邊界上，通過設(shè)置阻尼來吸收波動，通過彈簧來模擬地基彈性恢復(fù)力。粘彈性邊界是目前國內(nèi)應(yīng)用于波動動力分析的主要邊界之一，諸多學(xué)者對粘彈性人工邊界進(jìn)行研究和改進(jìn)并運用到實際問題的研究中，驗證其適用性并取得許多有意義的研究成果[2]。

二、有限元模型的建立

本文通過建立模型驗證設(shè)置的粘彈性人工邊界的精確度。分別對遠(yuǎn)置邊界模型和粘彈性邊界模型輸入正弦波，在中心設(shè)置觀測點，對比觀測結(jié)果。其中參數(shù)如下表1所示。

表1 模型的物理參數(shù)Table 1 Physical Parameters of Model

遠(yuǎn)置邊界模型中，左側(cè)長度L1=400，中間長度L0=20，右側(cè)長度L2=400，高度為H=20，寬為10。

粘彈性模型長L=820，高H=20,寬為10，剛度方程及阻尼方程分別為：

Kbt1=0.5*G4/R cbt1=P4*cs1

Kbn1=1*G4/R cbn1=P4*cp1

圖1 遠(yuǎn)置邊界模型示意圖Figure 1 Model Schematic Diagram of Distant Boundary

圖2 粘彈性邊界模型示意圖Fig.2 schematic diagram of viscoelastic boundary model

三、結(jié)果分析

分別對建立的遠(yuǎn)置邊界模型及粘彈性邊界模型輸入幅值為0.45的正弦波，得到觀測點的位移時程曲線如圖3(b)～(c)所示。

(a)輸入的正弦波 (b)遠(yuǎn)置邊界下的水平位移

(c)粘彈性邊界下的水平位移 (d)兩者對比 圖3 兩種邊界下的檢測位移曲線Fig.3 detection displacement curves under two boundaries

由圖3可知，截斷區(qū)域邊界處采用粘彈性人工邊界模擬計算結(jié)果與遠(yuǎn)置邊界模擬計算結(jié)果在峰值上完全吻合，在非峰值處曲線也具有高吻合度，而采用固定邊界的計算結(jié)果發(fā)生持續(xù)震蕩。遠(yuǎn)置邊界雖能得到逼近真實解的計算結(jié)果，但計算時間卻是采用粘彈性人工邊界的模型計算時間的50倍左右，這對于計算機資源產(chǎn)生極大消耗和浪費。表明采用彈簧-阻尼器組成的粘彈性人工邊界能較好地模擬截斷區(qū)域邊界特性，在進(jìn)行動分析時具有相當(dāng)高的計算精度且有利于節(jié)省計算機資源。

四、結(jié)論

通過分別建立模型、輸入地震波得到檢測點的位移時程曲線，得出采用彈簧-阻尼器組成的粘彈性人工邊界能較好地模擬截斷區(qū)域邊界特性，在進(jìn)行動分析時具有相當(dāng)高的計算精度且有利于節(jié)省計算機資源。