第一作者 樓夢(mèng)麟 男，博士，教授，博士生導(dǎo)師，1947年生

應(yīng)用通用程序計(jì)算深覆蓋土層地震反應(yīng)的幾個(gè)問題

樓夢(mèng)麟，邵新剛

(同濟(jì)大學(xué) 土木工程防災(zāi)國家重點(diǎn)實(shí)驗(yàn)室，上海200092)

摘要：討論了利用ANSYS和LS-DYNA程序進(jìn)行深覆蓋土層地震反應(yīng)時(shí)域計(jì)算時(shí)所涉及的幾個(gè)問題。綜合前期研究成果，通過對(duì)比通用程序的數(shù)值計(jì)算結(jié)果，分析了數(shù)值計(jì)算時(shí)土層計(jì)算范圍的選取及單元網(wǎng)格劃分、阻尼矩陣的形成等對(duì)計(jì)算精度的影響。分析結(jié)果表明：在確定深覆蓋土層阻尼矩陣的比例系數(shù)時(shí)，土層基頻與地震激勵(lì)主要特征頻率之間的關(guān)系是一個(gè)重要影響因素，當(dāng)前者小于后者，特別是遠(yuǎn)小于后者時(shí)，應(yīng)基于土層的基頻與地震激勵(lì)的主要特征頻率來形成阻尼矩陣的比例系數(shù)；豎向、橫向網(wǎng)格尺寸及人工邊界的選取會(huì)對(duì)計(jì)算的結(jié)果造成顯著影響。

關(guān)鍵詞：深覆蓋土層；地震反應(yīng)；土層范圍；有限元網(wǎng)格；Rayleigh阻尼矩陣

基金項(xiàng)目：國家科技支撐項(xiàng)目(2011BAG07B01);國家自然科學(xué)

收稿日期：2013-09-30修改稿收到日期：2014-06-19

中圖分類號(hào):TU449文獻(xiàn)標(biāo)志碼： A

Several problems in seismic response calculation of soil layer with deep deposit using general software

LOUMeng-lin,SHAOXin-gang(State Key Laboratory of Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China)

Abstract:Several problems in seismic response analysis of soil layer with deep deposit were discussed using ANSYS and LS-DYNA. The effects of some factors, such as, selection of calculation scope of soil layer, division of mesh elements, formation of damping matrix, etc. calculation accuracy were analyzed by comparing the numerical results obtained with the general programs. Numerical examples demonstrated that the relationship between the fundamental frequency of soil layer and the predominant frequency of earthquake excitation is an important factor to affect calculation accuracy; when the former is smaller than the latter, especially, much smaller than the latter, the proportional coefficients in damping matrix should be formed with the two frequencies; the vertical and horizontal grid sizes and the selection of artificial boundary have a significant influence on calculation results.

Key words: soil layer with deep deposit; seismic response; scope of soil layer; finite element mesh; Rayleigh damping matrix

復(fù)雜場地的土層地震反應(yīng)分析對(duì)合理確定上部結(jié)構(gòu)的地震動(dòng)輸入?yún)?shù)有重要意義[1]，是進(jìn)行重要工程場地地震安全性評(píng)價(jià)工作的主要環(huán)節(jié)之一[2]。在進(jìn)行土層地震反應(yīng)的確定性分析時(shí)，通常有兩種分析方法，即時(shí)域分析方法和頻域分析方法。在進(jìn)行時(shí)域分析時(shí)，一般采用粘滯阻尼假定，即阻尼力與質(zhì)點(diǎn)的運(yùn)動(dòng)速度成正比，通過求解實(shí)系數(shù)微分方程組，可求得各土層各自由度地震反應(yīng)的時(shí)程。而在頻域分析時(shí)，通常采用Fourier變換，通過求解復(fù)代數(shù)方程組，可求得土層地震反應(yīng)的穩(wěn)態(tài)解，此時(shí)一般采用滯后阻尼假定。這兩種分析方法各有不同的適用范圍，采用頻域法時(shí)，可以方便地對(duì)線性體系的土層地震反應(yīng)進(jìn)行計(jì)算，但對(duì)于非線性土層地震反應(yīng)的分析只能采用等效線性化方法[3]，無法進(jìn)行真非線性分析。而時(shí)域分析方法能夠適應(yīng)非線性體系地震反應(yīng)分析的不同要求，既可以采用等效線性化的分析方法，也可以采用真非線性分析方法[4]。除了簡單體系如梁、柱等外，一般工程動(dòng)力體系地震反應(yīng)的時(shí)域計(jì)算多采用直接積分法，一般又分為隱式算法和顯式算法。不同特點(diǎn)的程序采用不同方法，如在當(dāng)前應(yīng)用廣泛的ANSYS程序中采用了隱式直接積分的算法，而在LS-DYNA程序中采用的是顯式直接積分的算法[5]。深覆蓋土層(泛指土層厚度50 m以上或基本周期1 s以上的土層)地震反應(yīng)分析時(shí)地震工程中經(jīng)常遇到的問題，針對(duì)應(yīng)用ANSYS和LS-DNYA兩個(gè)通用程序進(jìn)行深覆蓋土層地震反應(yīng)分析時(shí)所涉及的幾個(gè)問題進(jìn)行討論。

1有限土層計(jì)算范圍

在實(shí)際工程應(yīng)用中，常把無限域中的地震波傳播問題轉(zhuǎn)化為離散體系的振動(dòng)問題，因此需要合理地把半無限域轉(zhuǎn)化為有限土層的計(jì)算體系，本節(jié)重點(diǎn)討論土層地震反應(yīng)計(jì)算中首先遇到的土層計(jì)算范圍有限化的問題，一般需解決以下幾個(gè)問題。

1.1地震輸入面的確定

在土層的地震反應(yīng)分析中，有時(shí)并不在基巖面輸入地震波，采用在土層的假想基巖面處輸入地震波的方式計(jì)算土層的動(dòng)力反應(yīng)。工程中常按以下原則確定假想基巖面[6]：假想基巖面下臥各層巖土的剪切波速均不小于500 m/s；當(dāng)?shù)孛? m以下存在剪切波速大于其上部各土層剪切波速2.5 m/s的土層，且該層及下臥各層巖土的剪切波速均不小于400 m/s時(shí)，可按地面至該土層頂面的距離確定。文獻(xiàn)[7]中以上海地區(qū)的4個(gè)不同深度鉆孔為例計(jì)算了不同假想基巖面時(shí)土層的自振周期和地震反應(yīng)。數(shù)值結(jié)果表明：以土層剪切波速500 m/s作為假想基巖面的判別依據(jù)并不科學(xué)，有時(shí)會(huì)產(chǎn)生50%左右的計(jì)算誤差且無規(guī)律可言。這是因?yàn)橥翆由疃鹊母淖儠?huì)使土層的自振頻率隨之變化，使得土層自振頻率與輸入地震波的頻譜特性之間的關(guān)系發(fā)生變化而致使地震反應(yīng)產(chǎn)生不同變化；同時(shí)從本文“4”的分析中可知，土層自振頻率不同也對(duì)土層阻尼矩陣的形成產(chǎn)生影響，也會(huì)影響地震反應(yīng)計(jì)算的準(zhǔn)確性。在當(dāng)前計(jì)算機(jī)性能大幅提升的條件下，取深覆蓋土層下臥基巖面作為土層地震反應(yīng)計(jì)算的激勵(lì)面是很有必要的，下臥基巖的彈性變形對(duì)深覆蓋土層地震反應(yīng)的影響可以忽略[8]。

1.2側(cè)向豎向人工邊界的選擇

在土層地震反應(yīng)計(jì)算和土-結(jié)構(gòu)動(dòng)力相互作用問題的計(jì)算中，土層兩側(cè)的邊界(或稱豎向邊界)問題曾是廣泛關(guān)注的課題?；谟?jì)算機(jī)能力有限的前提下，人們研究了多種人工邊界，以期大幅減少計(jì)算規(guī)模，實(shí)現(xiàn)深覆蓋土層的動(dòng)力計(jì)算。目前較常用的人工邊界有：粘性邊界、疊加邊界、旁軸邊界、一致邊界，粘彈性邊界、透射邊界等。由于粘性邊界物理意義清晰，處理方法簡單，雖然只有一階精度，但易于實(shí)現(xiàn)。因而一些通用程序得到應(yīng)用，如在LS-DYNA程序中，有一無反射邊界(non-reflecting boundary)即為粘性阻尼邊界，其基本思想就是利用邊界阻尼力模擬無限域的影響，邊界力定義為：

σ法向=-ρcdv法向

(1)

σ切向=-ρcsv切向

(2)

式中:ρ為材料密度，cd和cs分別為壓縮波速和剪切波速，v法向和v切向分別為邊界法向及切向質(zhì)點(diǎn)運(yùn)動(dòng)速度，σ法向和σ切向分別為無反射邊界作用的法向應(yīng)力及切向應(yīng)力。

已有研究表明[9-10]，在土層地震反應(yīng)這類外源激勵(lì)問題中，粘性阻尼邊界并不具備縮減土層范圍的有效性，設(shè)置遠(yuǎn)置側(cè)移邊界是十分必要的，建議至少應(yīng)遠(yuǎn)置5倍土層深度以外。本文以水平均勻土層與河谷地形土層為研究對(duì)象，進(jìn)一步通過應(yīng)用LS-DYNA計(jì)算所得的算例結(jié)果來驗(yàn)證上述文獻(xiàn)的研究結(jié)論。

圖1所示水平均勻土層幾何尺寸與物理特性如下：土層厚度h為10 m，土介質(zhì)的彈性模量391 MPa，泊松比0.35，質(zhì)量密度2 000 kg/m3，阻尼比0.05，剪切波速269 m/s。以地表中點(diǎn)A的加速度及位移反應(yīng)為研究對(duì)象，分別采用粘性邊界、側(cè)移邊界和固定邊界三種形式，邊界都在圖1的模型兩側(cè)施加。共計(jì)算了7種工況，分別取L/h=0.5、1、2.5、5、10、20和40。固定邊界是指將兩端人工邊界的節(jié)點(diǎn)三個(gè)方向自由度全部約束住，側(cè)移邊界是指不約束與水平地震激勵(lì)方向一致的水平自由度，約束其他方向自由度。粘性阻尼邊界即在自由人工邊界的基礎(chǔ)上再附加了產(chǎn)生阻尼力的粘性阻尼器。為了排除豎向網(wǎng)格對(duì)于計(jì)算精度的影響，本文考慮豎向網(wǎng)格d=λmin/32來進(jìn)行分析，λmin是指根據(jù)地震波有效截止頻率fmax所對(duì)應(yīng)的正弦波波長，本文輸入ElCentro波(見圖2，峰值調(diào)整為1m/s2)激勵(lì)，fmax取為10Hz，豎向網(wǎng)格均取d=1 m。

圖1　土層有限元網(wǎng)格圖 Fig.1 Finite element mesh of soil layer

圖2　 El Centro波 Fig.2 El Centro wave

河谷地形土層的材料參數(shù)與均勻土層相同，以河谷底部中部位置A點(diǎn)的加速度及位移反應(yīng)為研究對(duì)象，兩側(cè)豎向人工邊界也采用粘性邊界、側(cè)移邊界與固定邊界，河谷左岸、底部和右岸的土層厚度分別為69 m、65 m和71 m，幾何尺寸見圖3。該土層設(shè)計(jì)了6種計(jì)算工況，即L1/h1=L2/h2= 2.5、5、10、20、40、60，h1=69 m，h2=71 m。豎向網(wǎng)格劃分規(guī)則同水平均勻土層。

圖3　河谷土層 Fig.3 River valley soil layer

1.3土層計(jì)算模型的自振頻率

三種不同人工邊界下，計(jì)算所得到的土層基頻表1。由表1可知：不同人工邊界對(duì)水平均勻土層和河谷土層基頻的計(jì)算值影響顯著。固定邊界下土層基頻最大，自由邊界最小。對(duì)于水平均勻土層來說，在側(cè)移邊界條件下，計(jì)算模型所得基頻與解析解是相同的[11]。從中可以看出：土層的兩側(cè)設(shè)為固定邊界時(shí)，對(duì)土層基頻計(jì)算結(jié)果影響最大，自由邊界次之，側(cè)移邊界對(duì)于基頻的影響最小。對(duì)于固定邊界而言，當(dāng)L/h≥5時(shí)，邊界對(duì)于基頻的影響在5%以內(nèi)，隨著L/h的增加，土層基頻趨近于側(cè)移邊界的?？傮w上看，對(duì)于兩種土層，應(yīng)取L/h(L1/h1和L2/h2)≥5，與以前的分析結(jié)論相同。

表1　不同人工邊界時(shí)土層的基頻(單位:Hz)

1.4地震反應(yīng)計(jì)算結(jié)果

7種不同計(jì)算范圍下水平均勻土層頂點(diǎn)A處的加速度峰值和位移峰值見表2。對(duì)于水平均勻土層，可通過求解波動(dòng)方程得到其解析解[11]，其中加速度峰值a=2.15 m/s2，位移峰值u=1.16 mm。不同土層范圍有限元計(jì)算結(jié)果的相對(duì)誤差表示為：

(3)

式中:amax為水平均勻土層地表中點(diǎn)A處加速度(或位移)峰值有限元解，a*為波動(dòng)解得到的地表地震加速度(或位移)峰值。加速度和位移峰值的相對(duì)誤差見表3，從中可看出7種計(jì)算工況下有限元解的近似性的變化情況。

表2　 均勻土層A點(diǎn)加速度峰值和位移峰值

表3　均勻土層A點(diǎn)加速度峰值和位移峰值的相對(duì)誤差(%)

由表2和表3數(shù)據(jù)可知，當(dāng)采用固定邊界時(shí)，均勻土層地震反應(yīng)的計(jì)算結(jié)果并不隨土層范圍的擴(kuò)大而單調(diào)變化，這里涉及到不同的有限元土層范圍的土層自振頻率與輸入地震波頻率關(guān)系的變化因素，至少要取L為20倍土層深度h時(shí)，才能達(dá)到計(jì)算誤差控制在6%以內(nèi)；而采用粘性邊界時(shí)，數(shù)值計(jì)算結(jié)果變化情況與固定邊界相類似，總體上加速度的計(jì)算結(jié)果要優(yōu)于固定邊界，這是由于邊界阻尼器的設(shè)置，能夠消除部分人工邊界上的反射能量，但是效果仍不如側(cè)移邊界。對(duì)于固定邊界及粘滯邊界而言，在L/h很小時(shí)邊界范圍取值對(duì)土體基頻的影響非常顯著，這是造成在這兩種邊界條件下土層地震反應(yīng)計(jì)算結(jié)果誤差很大的主要原因。特別是在自由邊界條件下，當(dāng)土層范圍取得較小(如L/h≤1)時(shí)，土層基頻遠(yuǎn)遠(yuǎn)低于解析解，即使土層范圍取得較廣，但土層基頻仍與解析解有一定差距。我們知道，低頻分量對(duì)土層位移反應(yīng)的貢獻(xiàn)要遠(yuǎn)遠(yuǎn)大于對(duì)土層加速度反應(yīng)的貢獻(xiàn)[12]，因此自由邊界造成土層基頻變小致使在位移反應(yīng)中，粘滯邊界下的位移計(jì)算精度在很多情況下不及固定邊界，特別是在L/h=0.5時(shí)，土層基頻減低1.5倍，致使位移峰值增大3倍多。而側(cè)移邊界下土層基頻并不隨L/h的變化而變化，這樣也保證了計(jì)算的精度。顯然對(duì)于水平均勻土層，無論土層范圍取多大，側(cè)移邊界下的計(jì)算精度都很好，這是因?yàn)樗降卣鹱饔孟路謱泳鶆蛲翆拥卣鸱磻?yīng)分析為反對(duì)稱問題，兩側(cè)的側(cè)移邊界正好適應(yīng)了這一問題的對(duì)稱性條件，因此可用兩側(cè)具有側(cè)移邊界的單排單元的土條來計(jì)算水平地震作用下分層均勻土層的地震反應(yīng)[13]。

對(duì)于河谷地形土層，由于沒有解析解，以6種工況中L/h1=L/h2=60時(shí)的土層地震反應(yīng)為比較基準(zhǔn)，相關(guān)計(jì)算結(jié)果見表4、表5，可近似看出其他5種計(jì)算工況下河谷土層地震反應(yīng)解的近似性。

表4　河谷土層中A點(diǎn)的加速度峰值和位移峰值

表5不同計(jì)算范圍時(shí)A點(diǎn)加速度峰值

和位移峰值的相對(duì)誤差(%)

Tab.5 Relative error of peak acceleration and displacement at point A of River valley soil layer(%)

由表4和表5可以看出，河谷地形的地表A點(diǎn)的加速度或者位移在20倍邊界以上時(shí)基本沒有區(qū)別，誤差很小。對(duì)于固定邊界，當(dāng)L/h1=L/h2≤10時(shí)，計(jì)算結(jié)果誤差很大，粘性邊界相對(duì)于固定邊界，計(jì)算精度總體上要好很多，但是在L/h1=L/h2=2.5時(shí)，粘性邊界下土層位移反應(yīng)的誤差要大于固定邊界，其原因也是自由邊界下，土層基頻下降較多。對(duì)于側(cè)移邊界，當(dāng)2.5≤L/h1=L/h2≤40倍時(shí)，計(jì)算誤差都在2%以內(nèi)，精度較高。

結(jié)合以前和現(xiàn)在的分析結(jié)果，采用ANSYS、LS-DYNA進(jìn)行深覆蓋土層地震反應(yīng)計(jì)算時(shí)，宜采用遠(yuǎn)置側(cè)移人工邊界，一般單側(cè)范圍宜取土層厚度的5倍以上。土層在豎向地震激勵(lì)下的地震反應(yīng)計(jì)算精度對(duì)有限范圍大小取值的敏感性不及水平激勵(lì)時(shí)突出，實(shí)際工程計(jì)算中，如要進(jìn)行豎向地震輸入，必有水平地震輸入，因此按水平基底運(yùn)動(dòng)輸入所確定的土層有限范圍完全適用于豎向基底運(yùn)動(dòng)輸入時(shí)土層地震反應(yīng)的計(jì)算。

2計(jì)算網(wǎng)格劃分的要求

對(duì)于動(dòng)力分析問題，當(dāng)連續(xù)介質(zhì)離散為有限單元時(shí)，就有了濾波和散射作用[14]，動(dòng)力分析中應(yīng)重視單元?jiǎng)澐种械某叨葐栴}，使得在有效激勵(lì)頻率范圍內(nèi)，計(jì)算的精度得到保障。

2.1豎向網(wǎng)格尺寸

2.2橫向網(wǎng)格尺寸

3阻尼矩陣建模方式對(duì)地震反應(yīng)計(jì)算的影響

3.1地震反應(yīng)分析方程中的阻尼矩陣

在進(jìn)行復(fù)雜土層的地震動(dòng)力反應(yīng)時(shí)域分析時(shí)，一般采用集中質(zhì)量的有限元方法。對(duì)計(jì)算體系進(jìn)行離散化處理后，可以得到如下動(dòng)力方程式：

(1)

[C]=a0[M]+a1[K]

(2)

式中：a0和a1為比例常數(shù)。利用振型正交性，可得：

(3)

式中:ωm、ωn和ξm、ξn分別為計(jì)算體系第m、n階振型的自振頻率及對(duì)應(yīng)的阻尼比。振型阻尼比是動(dòng)力體系耗能特性的宏觀描述，對(duì)土層介質(zhì)來說易于從室內(nèi)土樣實(shí)驗(yàn)和實(shí)測試驗(yàn)中獲取。當(dāng)各階振型阻尼比相同時(shí)，即ξm=ξn=ξ，上式可簡化為:

(4)

這樣由式(4)確定的第i階振型阻尼比的計(jì)算值為：

(5)

式(5)表示的第i階振型阻尼比與第i階自振頻率之間的關(guān)系(見圖4)。

圖4　振型阻尼比與頻率的關(guān)系 Fig.4 Relationship of modal damping ratio and frequency

在土動(dòng)力分析中，土層體系的阻尼矩陣也常采用如下只與質(zhì)量矩陣成比例或與剛度矩陣成比例的形式。

(6)

(7)

這兩種單參數(shù)阻尼矩陣所產(chǎn)生的各階振型阻尼比隨頻率變化的趨勢也列入圖4中。很顯然，在質(zhì)量比例的阻尼矩陣下，動(dòng)力體系的各階振型阻尼比隨頻率升高而降低，而在剛度比例阻尼矩陣下，動(dòng)力體系的各階振型阻尼比隨頻率升高而線性增大。

由于深覆蓋土層的自振周期較長，因此其低階自振頻率遠(yuǎn)小于基巖地震波中的主要分量的激振頻率。已有研究表明：對(duì)于這類長周期動(dòng)力體系(指低階振型頻率遠(yuǎn)小于外部動(dòng)力荷載或動(dòng)力作用中主要分量的激勵(lì)頻率的動(dòng)力體系)在建立動(dòng)力體系阻尼矩陣時(shí)，必須重視外部動(dòng)力激勵(lì)中主要頻率成分的影響，對(duì)于深覆蓋土層地震反應(yīng)同樣需加以關(guān)注。下面根據(jù)我們的研究成果，結(jié)合ANSYS和LS-DYNA兩個(gè)計(jì)算軟件的阻尼建模功能，介紹進(jìn)行深覆蓋土層地震反應(yīng)分析時(shí)合理地建立阻尼矩陣應(yīng)關(guān)注的問題。

3.2ANSYS軟件

應(yīng)用ANSYS程序進(jìn)行深覆蓋土層地震反應(yīng)計(jì)算時(shí)，可采用多種形式的阻尼矩陣，其中最常用的是式(2)所示Rayleigh阻尼和式(7)所示的剛度比例阻尼矩陣。前者為雙參數(shù)模式，后者為單參數(shù)模式。在兩種計(jì)算模式并存的情況，雙參數(shù)模式的可調(diào)控性較好，建議在進(jìn)行深覆蓋土層地震反應(yīng)分析時(shí)不要選用式(7)，下面討論如何選雙參數(shù)。對(duì)于一般工程結(jié)構(gòu)而言，因?yàn)榈碗A振型在結(jié)構(gòu)的地震反應(yīng)中起主導(dǎo)作用，一般選取體系基頻及其所對(duì)應(yīng)的振型阻尼比來確定比例系數(shù)，即取ωm=ω1=2πf1；設(shè)ωn=2πf，如何選取f是建立Rayleigh阻尼矩陣應(yīng)加以特別關(guān)注的問題。通常選取動(dòng)力體系的第2階自振頻率f2或第3階自振頻率f3及其對(duì)應(yīng)的振型阻尼比。我們研究[18]表明對(duì)深覆蓋土層而言，這一選取方法是不可取的，應(yīng)選取基巖輸入地震波反應(yīng)譜的峰值頻率fR。這樣可以保障在土層第1階振型和主要激勵(lì)頻率相近的高階振型對(duì)土層地震反應(yīng)的貢獻(xiàn)的正確性，從而使深覆蓋土層地震反應(yīng)計(jì)算更加準(zhǔn)確。

3.3LS-DYNA軟件

應(yīng)用LS-DYNA程序進(jìn)行深覆蓋土層反應(yīng)計(jì)算時(shí)采用顯式算法的質(zhì)量比例阻尼矩陣進(jìn)行直接積分，為此取：[C]=a0[M]，并建議:

a0=2ω1ξ

(8)

式中：ω1為土層的基頻，ξ為土層的振型阻尼比，一般假定土層各階振型阻尼比相同。以ω1來確定比例系數(shù)a0符合人們的常規(guī)的認(rèn)識(shí)，即在結(jié)構(gòu)體系地震反應(yīng)中，第1階振型頻率貢獻(xiàn)最大。當(dāng)土層基頻接近于或高于輸入地震的主要頻率分量時(shí)，顯然由于土層第1階振型的振型參與系數(shù)大，其對(duì)體系地震反應(yīng)的貢獻(xiàn)遠(yuǎn)大于其他高階振型，因此當(dāng)?shù)?階振型阻尼比與實(shí)際阻尼比相等時(shí)，在第1階振型地震反應(yīng)計(jì)算正確得到保證的前提下，其他高階振型阻尼比偏小時(shí)使得土層地震反應(yīng)的計(jì)算誤差非常有限。而當(dāng)計(jì)算體系的基頻遠(yuǎn)低于輸入地震波的主要頻率分量時(shí)，很顯然與地震波中主要頻率分量相近的土層高階振型將存在共振效應(yīng)。眾所周知，共振放大因子為Di=1/(2ξi)，顯然當(dāng)?shù)凸勒裥妥枘岜葧r(shí)，將高估共振放大因子，從而使得該振型的地震反應(yīng)計(jì)算值增大。例如當(dāng)ω1/ωi=1/2時(shí)，則Di增大1倍，依此類推，此時(shí)高階振型的共振效應(yīng)不可忽略。根據(jù)我們的研究成果[19]，在質(zhì)量比例的阻尼模型中，應(yīng)對(duì)式(8)進(jìn)行修正，采用如下計(jì)算公式：

a0=2ωeξ

(9)

式中：ωe為待定的頻率參數(shù)。當(dāng)土層基頻ω1小于輸入基巖地震波反應(yīng)譜峰值頻率ωR(ωR=2πfR)時(shí)，建議采用取為土層基頻ω1和輸入地震波反應(yīng)譜峰值頻率ωR的算術(shù)平均值，基本可保證計(jì)算所得土層加速度峰值略高于真實(shí)解。

4結(jié)論

在深覆蓋土層地震反應(yīng)時(shí)域計(jì)算中，應(yīng)特別關(guān)注土層有效范圍的選取、計(jì)算網(wǎng)格劃分、阻尼矩陣形成方式等問題，處理不當(dāng)會(huì)嚴(yán)重影響計(jì)算結(jié)果的可靠性。

(1) 深覆蓋土層地震反應(yīng)計(jì)算中取下臥基巖面作為土層地震反應(yīng)計(jì)算的地震激勵(lì)面是很有必要的；土層計(jì)算時(shí)側(cè)移邊界的計(jì)算效果很好，宜采用遠(yuǎn)置側(cè)移人工邊界，單側(cè)范圍一般宜取5倍以上土層厚度。

(3) 在深覆蓋土層地震反應(yīng)的隱式計(jì)算中，Rayleigh阻尼矩陣建模時(shí)，除取其中一個(gè)振型頻率ωm為土層基頻外，建議另一個(gè)頻率ωn取為輸入基巖地震波反應(yīng)譜的峰值頻率ωR；在顯式計(jì)算中，形成質(zhì)量比例阻尼矩陣時(shí)，建議ωs取為土層基頻和輸入地震波反應(yīng)譜峰值頻率的算術(shù)平均值，基本可保證計(jì)算所得土層加速度峰值略高于真實(shí)解。

參 考 文 獻(xiàn)

[1] 廖振鵬．近場波動(dòng)問題的有限元解法[J]．地震工程與工程振動(dòng)，1984, 4(2)：1-14．

LIAO Zhen-peng. A finite element method for near-field wave motion in heterogeneous materials[J]. Earthquake Engineering andEngineering Vibration, 1984, 4(2): 1-14.

[2] GB 17741-2005. 工程場地地震安全性評(píng)價(jià)[S]. 北京：地震出版社，2012.

[3] Martin P P, Seed H B. One-dimensional dynamic ground response analyses [J]. Journal of the Geotechnical Engineering Division, 1982, 108(7): 935-952.

[4] 李小軍，廖振鵬，關(guān)慧敏．粘彈性場地地形對(duì)地震動(dòng)影響分析的顯式有限元-有限差分方法[J]．地震學(xué)報(bào)，1995, 17(3)：362-369．

LI Xiao-jun, LIAO Zhen-peng, GUAN Hui-min. An explicit finite element-finite difference method for analyzing the effect of visco-elastic An explicit finite element-finite difference method for analyzing the effect of visco-elastic Local topograph on the earthquake motion [J]. ACTA Seismologica Sinica, 1995, 17(3): 362-369.

[5] 尚曉江，蘇建宇，王化鋒．ANSYS/LS-DYNA 動(dòng)力分析方法與工程實(shí)例[M]．北京：中國水利出版社，2008.

[6] GB50011-2010. 建筑抗震設(shè)計(jì)規(guī)范 [S]. 北京：中國建筑工業(yè)出版社，2010．

[7] 上海市地震局、同濟(jì)大學(xué)．上海市地震動(dòng)參數(shù)區(qū)劃[M]．北京：地震出版社，2004.

[8] 潘旦光，樓夢(mèng)麟．基巖彈性剛度對(duì)土層地震反應(yīng)的影響[J]．地震工程與工程振動(dòng)，2004(4)：158-163．

PAN Dan-guang, LOU Meng-lin. Influence of rock stiffness on seismic response of layered soil [J]. Earthquake Engineering and Engineering Vibration, 2004(4): 158-163.

[9] 樓夢(mèng)麟，潘旦光，范立礎(chǔ)．土層地震反應(yīng)分析中側(cè)向人工邊界的影響[J]．同濟(jì)大學(xué)學(xué)報(bào)，2003, 7(7):757-761．

LOU Meng-lin, PAN Dan-guang, FAN Li-chu. Effect of vertical artificial boundary on seismic response of soil layer[J]. Journal of Tongji University , 2003, 7(7): 757-761.

[10] 喬冠東，樓夢(mèng)麟．大型沉井基礎(chǔ)地震反應(yīng)分析中側(cè)向邊界的影響[J]．工程力學(xué)，2012, 26(S1): 111-114．

QIAO Guan-dong, LOU Meng-lin. Effect of vertical artifical boundary on seismic responses analysis of large-scale open caisson [J]. Engineering Mechanics, 2012, 26(S1): 111-114.

[11] 潘旦光，樓夢(mèng)麟，董聰．一致輸入作用下土層的地震反應(yīng)分析[J]．計(jì)算力學(xué)學(xué)報(bào)，2005, 22(05): 562-567．

PAN Dan-guang, LOU Meng-lin, DONG Cong. Seismic response analysis of soil layer under uniform excitation[J]. Chinese Journal of Computational Mechanics, 2005, 22( 05):562-567.

[12] 樓夢(mèng)麟，董云．加速度反應(yīng)譜規(guī)準(zhǔn)化對(duì)場地位移反應(yīng)的影響[J]．震災(zāi)防御技術(shù)，2014，9(2)：149-158.

LOU Meng-lin, DONG Yun. Effect of normalizing acceleration response spectra on displacement time history[J]. Technology for Earthquake Disaster Prevention, 2014, 9(2): 149-158.

[13] 鄔都，樓夢(mèng)麟．水平成層土層地震反應(yīng)分析的一維有限元方法[J]．震災(zāi)防御技術(shù)，2008( 01): 45-52．

WU Dou, LOU Meng-lin. Seismic response analysis in horizontal layer with 1-D finite element method [J]. Technology for Earthquake Disaster Prevention, 2008(01): 45-52.

[14] 宗福開．波傳播問題中有限元分析的頻散特性及離散化準(zhǔn)則[J]．爆炸與沖擊，1984, 4(4): 16-21．

ZONG Fu-kai. Frequency-dispersion characteristics and discretization of the finite element analysis in wave propagation problems [J]. Explosion and Shock Waves, 1984, 4(4): 16-21.

[15] 白建方，樓夢(mèng)麟．復(fù)雜場地土層有限元建模的兩個(gè)問題研究[J]．巖土工程界，2008, 11(9): 27-30．

BAI Jian-fang, LOU Meng-lin. The two problems research of finite element modeling of complex ground soil [J]. Geotechnical Engineering World, 2008, 11(9): 27-30.

[16] 楊燕，樓夢(mèng)麟．汶川地震中遠(yuǎn)場深覆蓋土層動(dòng)力反應(yīng)分析[J]．防災(zāi)減災(zāi)工程學(xué)報(bào)，2011, 31(4): 462-468．

YANG Yan, LOU Meng-lin. Dynamic response analysis of far-field soil layer with deep deposit during wenchuan earthquake [J]. Journal of Disaster Prevent and Mitigation ENG, 2011, 31(4): 462-468.

[17] 白建方，樓夢(mèng)麟．行波輸入下土層有限元網(wǎng)格的劃分原則[J]．世界地震工程，2009, 25(1): 48-51．

BAI Jian-fang, LOU Meng-lin. Partition principle of the finite element meshes of soil under traveling wave excitations [J]. World Information on Earthquake Engineering, 2009, 25(1): 48-51.

[18] 樓夢(mèng)麟，邵新剛. 深覆蓋土層Rayleigh阻尼矩陣建模問題的討論[J]. 巖土工程學(xué)報(bào),2013, 35(7): 1272-1279.

LOU Meng-lin,SHAO xin-gang. Discussion on modeling issues of Rayleigh damping matrix in soil layers with deep deposit [J]. Chinese Journal of Geotechnical Engineering,2013, 35(7): 1272-1279.

[19] 樓夢(mèng)麟，邵新剛. 土層地震反應(yīng)顯式計(jì)算中阻尼矩陣系數(shù)的選取[J]. 同濟(jì)大學(xué)學(xué)報(bào),2013,41(8):1126-1132.

LOU Meng-lin, SHAO Xin-gang.Coefficient selection of damping matrix in explicit calculation of soil layer seismic response [J]. Journal of Tongji University, 2013,41(8): 1126-1132.