復(fù)雜高層結(jié)構(gòu)隨機地震反應(yīng)中阻尼模型的影響分析

2014-09-27 13:50:19楊永勝陳清軍

湖南大學(xué)學(xué)報·自然科學(xué)版 2014年6期

關(guān)鍵詞：抗震設(shè)計阻尼

楊永勝+陳清軍

文章編號：16742974(2014)06002807

收稿日期：20130822

基金項目：國家自然科學(xué)基金資助項目(50978198)

作者簡介：楊永勝（1984-），男，山西忻州人，同濟大學(xué)博士研究生

通訊聯(lián)系人，E-mail:chengj@#edu.cn

摘要：扼要比較了Rayleigh阻尼和Caughey阻尼兩種阻尼模型，利用ANSYS軟件中的APDL編程語言進行二次開發(fā)，將Caughey阻尼模型應(yīng)用到復(fù)雜高層建筑結(jié)構(gòu)的隨機地震反應(yīng)分析中；以某實際復(fù)雜高層建筑為例，采用現(xiàn)場測試和數(shù)值模擬相結(jié)合的方法，分析了不同阻尼模型對復(fù)雜高層結(jié)構(gòu)隨機地震反應(yīng)計算結(jié)果的影響，討論了Rayleigh阻尼模型和Caughey阻尼模型的適用范圍.結(jié)果表明，Rayleigh阻尼模型的曲線形狀與所選取的控制頻率相關(guān)，兩控制頻率間的計算阻尼比小于實測結(jié)果，兩控制頻率以外部分的計算阻尼比則大于實測阻尼比；Caughey阻尼模型能在較大頻率范圍內(nèi)較準確地反映結(jié)構(gòu)的阻尼特性.



關(guān)鍵詞：阻尼；抗震設(shè)計；隨機地震反應(yīng)；ANSYS二次開發(fā)；振動測試

中圖分類號：TU311.3 文獻標識碼：A

Influence of the Damping Model of Complex Highrise

Building Structure on Its Random Seismic Response



YANG Yongsheng1,2，CHEN Qingjun1,2

(1.State Key Laboratory of Disaster Reduction in Civil Engineering, Tongji Univ, Shanghai 200092, China;

2.Research Institute of Structural Engineering and Disaster Reduction, Tongji Univ, Shanghai 200092, China)

Abstract:A comparative analysis between Rayleigh damping model and Caughey damping model was first conducted, and then, by ANSYS secondary developing based on APDL programming language, the Caughey damping model was successfully applied to the random seismic response analysis of complex highrise building structure. An actual complex highrise building structure was selected as the research object. By using the combination method of field vibration test and numerical simulation, the influence of different damping models on the random seismic response of complex highrise building was analyzed, and the applicability and effectiveness of the two damping models were also discussed. The results have indicated that the curve shape of Rayleigh damping model is closely related to the control frequencies, the calculated damping ratio is lower than the actual test damping ratio between two control frequencies, and it is contrary among other frequency ranges. The Caughey damping model can more accurately meet the actual damping in a wider frequency range.

Key words:damping; seismic design; random seismic response; ANSYS secondary development; vibration test



近年來，國內(nèi)外相繼興建了許多復(fù)雜高層及超高層建筑，這些特殊的高層及超高層建筑結(jié)構(gòu)的組成成分、耗能機理和振動特性都比較復(fù)雜，其阻尼矩陣的形成也有其特殊性[1].目前，在高層建筑地震反應(yīng)分析中往往采用Rayleigh阻尼模型來構(gòu)建阻尼矩陣，對于一般結(jié)構(gòu)，往往選取結(jié)構(gòu)的第1和第2階振型參數(shù)來形成Rayleigh阻尼矩陣.但對于復(fù)雜高層及超高層結(jié)構(gòu)而言，高階振型對結(jié)構(gòu)地震反應(yīng)的影響將不容忽視，而且對于復(fù)雜高層結(jié)構(gòu)的某些構(gòu)件，在地震反應(yīng)分析中起控制作用的可能不再是低階振型，因此如果采用只考慮2個低階振型的Rayleigh阻尼模型，將會導(dǎo)致結(jié)構(gòu)的地震反應(yīng)分析結(jié)果產(chǎn)生較大的偏差[2].

本文將首先扼要比較Rayleigh阻尼和Caughey阻尼2種阻尼模型，利用ANSYS軟件中的APDL編程語言進行二次開發(fā)，以將Caughey阻尼模型應(yīng)用到復(fù)雜高層建筑結(jié)構(gòu)的隨機地震反應(yīng)分析中，并以某實際復(fù)雜高層建筑為例，采用現(xiàn)場測試和數(shù)值模擬相結(jié)合的方法，分析不同阻尼模型對復(fù)雜高層結(jié)構(gòu)隨機地震反應(yīng)計算結(jié)果的影響.

1 兩種阻尼模型比較

在結(jié)構(gòu)動力分析中，一般假定阻尼力的大小與質(zhì)點振動的速度成正比，即采用黏滯阻尼模型[3-4].為了應(yīng)用振型疊加法，一般假定振型關(guān)于阻尼矩陣C正交.目前，多自由度體系的阻尼矩陣主要有Rayleigh阻尼和Caughey阻尼2種模型.

Rayleigh阻尼模型作為目前應(yīng)用最為廣泛的阻尼模型，該模型假定結(jié)構(gòu)的阻尼矩陣是質(zhì)量矩陣和剛度矩陣的線性組合，通常稱為比例阻尼：

C=a0M+a1K．(1)

式中：a0和a1為比例常數(shù)．利用振型正交性，可得：

a0a1=2ωmωnω2n－ω2mωn－ωm－1ωn1ωmξmξn.(2)

式中：ωm,ωn, ξm,ξn分別為體系第m和n階振型的自振頻率及對應(yīng)的振型阻尼比．振型阻尼比是結(jié)構(gòu)體系耗能特性的宏觀描述，易于從實驗和實測試驗中獲?。敻麟A阻尼比相同時，即ξm=ξn=ξ，式(2)可簡化為：

a0a1=2ξωn+ωmωmωm1.(3)

對于一般工程結(jié)構(gòu)而言，因為低階振型在結(jié)構(gòu)的動力反應(yīng)中起主導(dǎo)作用，通常取低階振型來確定比例系數(shù)a0和a1，其他各階振型的阻尼比可通過式(4)求解：

ξj=12a0ωj+a1ωj.(4)

式中：ωj為第j階的自振頻率; ξj為相應(yīng)于第j階的阻尼比．

Caughey阻尼又稱擴展的Rayleigh阻尼．若要求更多階振型的阻尼比滿足實測阻尼比，則必須構(gòu)造更多項的線性組合，Caughey阻尼能較好地滿足這個要求．Caughey阻尼可寫成如下形式[3]：

C=M∑b∈Sab(M-1K)b=∑b∈SCb.(5)

式中:ab是Caughey阻尼模型中的比例系數(shù)，集合S包含n個任意整數(shù)值，n的大小和S的具體元素取決于阻尼比曲線擬合結(jié)果．根據(jù)黏滯阻尼假定，推導(dǎo)得出項b對于廣義阻尼值的貢獻為：

c*nb=φnTCbφn=abω2bnm*n. (6)

在此基礎(chǔ)上，與任意振型n相關(guān)的廣義阻尼為：

c*n=∑b∈Sc*nb=∑b∈Sabω2bnm*n. (7)

同時，利用振型關(guān)于阻尼矩陣正交的假定，亦可推導(dǎo)得出第n階振型的廣義阻尼為：

c*n=φnTCφn=2ξnωnm*n.(8)

從中可得出第n階自振頻率ωn對應(yīng)的振型阻尼比為：

ξn=c*n2ωnm*n=12ωn∑b∈Sabω2bn. (9)

上式提供了根據(jù)指定振型數(shù)目、自振頻率和阻尼比計算系數(shù)ab的途徑：將已知的n個自振頻率以及實測的n階阻尼比代入該式，得到n階線性代數(shù)方程組，從中可以解得n個待定常數(shù)ab值；然后按式(5)形成阻尼矩陣，這樣使得參與計算的n階振型的阻尼比與實際量測結(jié)果吻合.理論上S可包含-∞到+∞之間的任意整數(shù)值，但實踐上，S中元素值選擇越接近于零越好，這樣有助于保證數(shù)值計算的穩(wěn)定，當S={0，l}時，即為式(1)所示的Rayleigh阻尼模型.本文利用ANSYS軟件中的APDL編程語言進行二次開發(fā)，從而借助ANSYS軟件將Caughey阻尼模型應(yīng)用到復(fù)雜高層建筑結(jié)構(gòu)的隨機地震反應(yīng)分析中.

2 某高層建筑結(jié)構(gòu)動力特性測試

本文采用同濟大學(xué)某復(fù)雜高層建筑為例，進行阻尼模型對復(fù)雜高層結(jié)構(gòu)隨機地震反應(yīng)結(jié)果的影響分析.該建筑上部結(jié)構(gòu)為21層巨型鋼框架結(jié)構(gòu)體系.框架的層高為4 m，每邊9跨，跨度為5.4 m，結(jié)構(gòu)平面尺寸為48.6 m×48.6 m，主體結(jié)構(gòu)采用鋼框架加外圍支撐的形式.主樓樓層中部樓板缺失，形成L型3平面，往上順時針旋轉(zhuǎn)，每3層旋轉(zhuǎn)90°，到頂層共旋轉(zhuǎn)6次.該建筑較多樓層的樓板開洞面積大于樓板面積的35%，有效樓板寬小于典型樓板寬度的40%；結(jié)構(gòu)沿高度每隔3層設(shè)2 m高設(shè)備層，豎向剛度突變，該建筑屬于典型的平面不規(guī)則和豎向剛度不規(guī)則的復(fù)雜高層結(jié)構(gòu).柱子采用鋼管混凝土柱，截面尺寸900 mm×900 mm，梁采用H型鋼(700 mm ×300 mm×24 mm×13 mm)(高度×寬度×腹板厚度×翼緣厚度)，弦桿采用H型鋼(500 mm×200 mm×16 mm×10 mm)(高度×寬度×腹板厚度×翼緣厚度)，腹桿采用H型鋼(200 mm×200 mm×12 mm×8 mm)(高度×寬度×腹板厚度×翼緣厚度)，支撐圓管(φ300 mm×16 mm)，表1為結(jié)構(gòu)主要構(gòu)件的材料參數(shù).

表1 結(jié)構(gòu)部件材料參數(shù)

Tab.1 Material parameters of structural parts

結(jié)構(gòu)名稱

材料類型

彈性模量

/MPa

密度

/(kg?m-3)

泊松比

承臺

混凝土

3.0×104

2.5×103

0.2

柱

鋼筋混凝土

4.7×104

3.2×103

0.2

梁

鋼材

2.1×105

7.8×103

0.3

弦桿

鋼材

2.1×105

7.8×103

0.3

腹桿

鋼材

2.1×105

7.8×103

0.3

支承

鋼材

2.1×105

7.8×103

0.3

目前結(jié)構(gòu)自振特性主要測試方法有穩(wěn)態(tài)正弦激振法、傳遞函數(shù)法、自由振動法和脈動測試法.其中脈動測試法是利用環(huán)境隨機振動作為結(jié)構(gòu)物激振的振源來測定并分析結(jié)構(gòu)物固有特性，是近年來隨著計算機技術(shù)及FFT理論的普及而發(fā)展起來的，現(xiàn)已被廣泛應(yīng)用于建筑物的動力分析研究中.本文采用脈動測試法對該復(fù)雜高層結(jié)構(gòu)進行了動力特性測試，所用儀器為美國凱尼公司生產(chǎn)的Basalt型數(shù)字振動測試儀．為了排除高頻背景噪聲的干擾，使得測試結(jié)果的低頻成分更真實，本次測試將儀器采樣頻率調(diào)整到50 Hz.地脈動激勵可以看作為各態(tài)歷經(jīng)的平穩(wěn)隨機過程，只要測試時間足夠長，可以用單個樣本函數(shù)上的時間平均來描述這個過程的所有樣本的平均特性．本次測試記錄時間設(shè)定為20 min以上，本次測試的測點選擇在該高層結(jié)構(gòu)21層樓梯間拐角處，其中某測試樣本的時程如圖1所示.

脈動激勵的頻譜通常情況下是比較平坦的，可

時間/min(a) 東西方向

時間/min(b) 南北方向

圖1測試得到的某樣本時程

Fig.1 A time history sample of the field test results

以看成是有限帶寬的白噪聲，建筑結(jié)構(gòu)振動信號的功率譜在自振頻率處會出現(xiàn)較大的峰值.將測試所得的信號求其自功率譜，通過峰值所在頻率點來識別結(jié)構(gòu)的自振頻率.對各階自振頻率對應(yīng)的阻尼比，可采用半功率法求取，由半功率點幅值Hω/2確定半功率帶寬Δω，然后根據(jù)式(10)計算：

ξk=Δω2ωk. (10)

式中: ωk為第k階自振頻率; H(ωk)為k階自振頻率處的峰值.圖1所示測試樣本的功率譜曲線見圖2，表2給出了本次測試的結(jié)構(gòu)自振頻率和阻尼比.

頻率/Hz(a) 東西方向

頻率/Hz(b) 南北方向

圖2 自功率譜曲線

Fig.2 Autopower spectrum curves 

表2 結(jié)構(gòu)自振頻率和阻尼比的測試結(jié)果

Tab.2 Natural frequency and damping ratio of the field test

方向

自振頻率/Hz

阻尼比/%

1階

2階

3階

1階

2階

3階

東西方向

0.47

1.62

2.82

3.15

2.11

2.46

南北方向

0.50

1.75

2.78

3.38

1.47

1.98

3 某復(fù)雜高層結(jié)構(gòu)有限元模型及隨機地震

動輸入

3.1 結(jié)構(gòu)有限元模型的建立

為進行復(fù)雜高層結(jié)構(gòu)的隨機地震反應(yīng)分析，采用ANSYS軟件建立三維有限元模型[5]，其中梁、柱、斜撐和腹桿采用Beam188單元來模擬，樓板用Shell63單元來模擬，建立的某復(fù)雜高層建筑結(jié)構(gòu)三維有限元模型如圖3所示.

由于地脈動激勵強度較弱，一般低階自振頻率和對應(yīng)的阻尼比精度較高，與實際比較相符.本文采用東西方向和南北方向第一階自振頻率對應(yīng)阻尼比的平均值作為該復(fù)雜高層結(jié)構(gòu)的參考阻尼比，具體取值為3.3%.

數(shù)值模型的自振特性見表3，將前三階自振頻率與實測結(jié)果比較，第一階的相對誤差為2%，第二

階的相對誤差為5.2%，第三階的相對誤差為7.4%，表明本次有限元模型具有較高精度.表中同時也給出了振型對應(yīng)的振型參與系數(shù)和累積有效質(zhì)量參與系數(shù).復(fù)雜高層結(jié)構(gòu)進行抗震設(shè)計時，要求各振型累積參與質(zhì)量之和不低于90%[6-7]，為了更精確地比較2種阻尼模型下復(fù)雜高層結(jié)構(gòu)地震反應(yīng)的不同以及考慮到高階振型對復(fù)雜高層結(jié)構(gòu)地震反應(yīng)的影響，本文取前150階振型參與計算，輸入激勵的頻率取0~20 Hz.



圖3 結(jié)構(gòu)三維有限元模型

Fig.3 The three dimension finite element model

表3 結(jié)構(gòu)自振特性表

Tab.3 Natural vibration characteristics of the finite element model

階數(shù)

100

150

頻率/Hz

0.48

0.50

0.76

1.54

1.73

2.43

2.61

2.85

4.04

4.88

6.39

9.91

13.29

15.10

17.64

19.14

周期/s

2.08

1.96

1.32

0.65

0.58

0.41

0.38

0.34

0.25

0.21

0.16

0.10

0.08

0.07

0.06

0.05

東西方向振型

參與系數(shù)

1.00

0.26

0.18

0.38

0.09

0.04

0.19

0.01

0.02

0.01

0.04

0.01

南北方向振型

參與系數(shù)

0.27

1.00

0.17

0.10

0.44

0.05

0.01

0.25

0.03

0.08

0.01

東西方向累積有

效質(zhì)量參與系數(shù)

0.73

0.78

0.88

0.89

0.92

0.96

0.97

0.98

0.99

1.00

南北方向累積有

效質(zhì)量參與系數(shù)

0.05

0.73

0.75

0.76

0.89

0.93

0.95

0.97

0.98

0.99

1.00

3.2 隨機地震動輸入模型

本文采用杜修力于1998年提出的修正地表隨機地震動模型作為輸入激勵，其功率譜密度表達式為：[8]

S(ω)=1+4ξ2gωωg21－ωωg22+4ξ2gωωg2×

11+(Dω)2×ω4(ω2+ω20)2S0.(11)

式中:ωg和ξg分別為場地土的卓越頻率和阻尼比;ω0≈1.83為低頻拐角頻率;D=1/28π為與震源有關(guān)的參數(shù).杜修力修正模型假定基巖地震動為有色譜，在大量實際地震記錄的統(tǒng)計基礎(chǔ)上，加入高通濾波器ω4/ω2+ω202和低通濾波器1/［1+Dω2］，合理修正了高頻能量和低頻能量，與實際地震動特征比較相符.

4 兩種不同阻尼模型下的計算結(jié)果分析

4.1 Raleigh阻尼模型下的計算結(jié)果分析

首先，取結(jié)構(gòu)實測的參考阻尼比3.3%作為各階振型阻尼比，根據(jù)振型疊加原理，利用隨機振動CQC(complete quadratic combination)法進行計算，其計算結(jié)果作為參考真實解.選取5組Rayleigh阻尼模型進行計算，第一組取結(jié)構(gòu)1階、4階自振頻率作為控制頻率，記作R1；第二組取結(jié)構(gòu)1階、10階自振頻率作為控制頻率，記作R2；第三組取結(jié)構(gòu)1階、25階自振頻率作為控制頻率，記作R3；第四組取結(jié)構(gòu)1階、40階自振頻率作為控制頻率，記作R4；最后一組取結(jié)構(gòu)的1階、150階自振頻率作為控制頻率，記作R5．各組阻尼模型的擬合曲線見圖4, 從圖中可看出,阻尼曲線的形狀主要與控制頻率點及兩控制頻率之間的間距有關(guān)：兩控制頻率之間的部分計算阻尼小于實際的結(jié)構(gòu)阻尼,兩控制頻率以

外的部分計算阻尼則大于實際的結(jié)構(gòu)阻尼.所以,控制頻率的合理選取是計算結(jié)果準確與否的關(guān)鍵.

分別計算各Raleigh阻尼模型工況下結(jié)構(gòu)的隨機振動響應(yīng)，選取結(jié)構(gòu)頂層節(jié)點作為參考點，圖5為結(jié)構(gòu)頂層的位移功率譜計算結(jié)果；分別在每層樓板處選取一個參考點，計算其隨機振動響應(yīng)，然后根據(jù)隨機振動理論，通過計算各響應(yīng)量的各階譜矩，來估算響應(yīng)最大值的均值，圖6為樓層的位移極值包絡(luò)圖.

頻率/Hz

圖4 不同控制頻率下的Rayleigh阻尼模型

Fig.4 Rayleigh damping model under 

different control frequencies



頻率/Hz(a)東西方向

頻率/Hz(b)南北方向

圖5 頂層位移功率譜密度

Fig.5Displacement power spectral density of roof response



位移/m(a)東西方向

位移/m(b)南北方向

圖6 樓層位移極值包絡(luò)圖

Fig.6Maximum relative displacement response along structural height



從圖5和圖6可看出由結(jié)構(gòu)1階、5階自振頻率為阻尼控制頻率(R2)的計算結(jié)果與真實解最為接近，由結(jié)構(gòu)前兩階自振頻率為阻尼控制頻率(R1)的計算結(jié)果明顯偏小，說明其計算阻尼偏大，低估了

結(jié)構(gòu)的響應(yīng).其余3種情況(R3,R4,R5)的計算結(jié)果比真實值偏大，說明其計算阻尼比實際阻尼要小.

4.2 Caughey阻尼模型下的計算結(jié)果分析

由式(9)可知，當i=2時，Caughey阻尼模型將與Raleigh阻尼模型相同，本文分別選取i=5,7,9三種情況確定的Caughey阻尼模型進行計算.i=5時，分別選取結(jié)構(gòu)的第1,5,10,15,25階自振頻率作為控制頻率；i=7時，分別選取結(jié)構(gòu)的第1, 5, 10, 15, 25, 40, 60階自振頻率作為控制頻率；i=9時，分別選取結(jié)構(gòu)的第1,5,10,15,25,40,60,100,150階自振頻率作為控制頻率；各組阻尼模型的擬合曲線見圖7.

頻率/Hz

圖7 不同控制頻率下的Caughey阻尼模型

Fig.7 Caughey damping model under 

different control frequencies

圖8為各組Caughey阻尼模型下結(jié)構(gòu)頂層的位移功率譜計算結(jié)果，圖9為樓層位移極值包絡(luò)圖.從圖中可以看出9個基本頻率點控制的Caughey阻尼模型計算結(jié)果與參考真實值最為接近，隨著控制頻率點數(shù)的增加，計算結(jié)果越來越靠近真實值.由此可見，控制頻率點的數(shù)目及控制頻率選取的合理性是決定Caughey阻尼模型精確與否的關(guān)鍵.

頻率/Hz(a)東西方向

頻率/Hz(b)南北方向

圖8 頂層位移功率譜密度

Fig.8 Displacement power spectral 

density of roof response



位移/m(a)東西方向

位移/m(b)南北方向

圖9 樓層位移極值包絡(luò)圖

Fig.9 Maximum relative displacement response 

along structural height

對比圖5,圖6與圖8, 圖9可知，采用Caughey阻尼模型的計算結(jié)果比采用Reyleigh阻尼模型的計算結(jié)果更接近真實解，Caughey阻尼在廣義上可以看作是對已知的阻尼比頻率數(shù)據(jù)以冪函數(shù)為底的非線性擬合，能在較大頻域范圍內(nèi)更準確地反映結(jié)構(gòu)的阻尼特性.在進行隨機地震反應(yīng)分析時，采用Caughey阻尼的計算結(jié)果更能反映實際情況.

5 結(jié) 語

本文以某實際復(fù)雜高層建筑為例，采用現(xiàn)場測試和數(shù)值模擬相結(jié)合的方法，分析了不同阻尼模型對復(fù)雜高層結(jié)構(gòu)隨機地震反應(yīng)計算結(jié)果的影響，討論了Rayleigh阻尼模型和Caughey阻尼模型的適用范圍.本文的分析表明：

1) Rayleigh阻尼模型的曲線形狀與所選取的控制頻率相關(guān)，兩控制頻率間的計算阻尼比小于實測結(jié)果，兩控制頻率以外部分的計算阻尼比則大于實測阻尼比；對于高階振型對地震反應(yīng)貢獻較大的復(fù)雜高層建筑而言，選取前兩階振型作為Rayleigh阻尼的控制頻率會產(chǎn)生較大誤差.

2) Caughey阻尼模型在廣義上可以看作是對已知的阻尼比頻率數(shù)據(jù)以冪函數(shù)為底的非線性擬合，能在較大頻域范圍內(nèi)更準確地反映結(jié)構(gòu)的阻尼特性.對于高階振型對地震反應(yīng)貢獻較大的復(fù)雜高層建筑而言，在結(jié)構(gòu)隨機地震反應(yīng)分析中，采用Caughey阻尼模型可獲得更高精度的計算結(jié)果.

3) 本文分析的2種阻尼模型(Rayleigh阻尼模型和Caughey阻尼模型)屬于比例阻尼范疇.對于多于1種材料組成的結(jié)構(gòu)，由于不同材料在結(jié)構(gòu)的不同部分提供的能量損失機制差別較大，將導(dǎo)致阻尼不是成比例的；對此，可將每一個明顯的結(jié)構(gòu)組成部分都建立一個比例阻尼矩陣，然后將其直接集裝形成組合的體系矩陣加以處理.對于更一般的廣義比例阻尼及非比例阻尼問題則需要進一步研究.

參考文獻

［1］段紹偉，向湘林，沈蒲生. 復(fù)雜阻尼結(jié)構(gòu)阻尼模型研究[J]. 振動與沖擊，2011，30(2)：73-76.

DUAN Shaowei, XIANG Xianglin, SHEN Pusheng. Damping model of a complex structure[J]. Journal of Vibration and Shock, 2011, 30(2): 73-76.(In Chinese)

［2］沈飛，樓夢麟. 超高層建筑地震反應(yīng)中高階振型影響分析[J].工程力學(xué)，2012，29(SⅠ)：23-28.

SHEN Fei, LOU Menglin. Influence of high modes of super highrise building on its seismic responses[J]. Engineering Mechanics, 2012, 29(S I): 23-28. (In Chinese)



［3］ CLOUGH R W, PENZIEN J. Dynamics of structures[M]. New York: McGrawHill, 1993:183-193.

［4］ CHOPRA A K. Matrix dynamics of structures: theory and applications to earthquake engineering[M]. 2nd ed. New Jersey: PrenticeHall, 2001:454-462.

［5］王新敏. ANSYS工程結(jié)構(gòu)數(shù)值分析[M].北京：人民交通出版社，2007:542-553.

WANG Xinmin. ANSYS numerical analysis of engineering structure[M]. Beijing: China Communications Press, 2007:542-553. (In Chinese)

［6］ GB 50011—2010建筑抗震設(shè)計規(guī)范[S]. 北京：中國建筑工業(yè)出版社，2010:320-325.

GB 50011—2010 Code for seismic design of buildings [S]. Beijing: China Architecture & Building Press, 2010:320-325. (In Chinese)

［7］ WILSON E L. Threedimensional static and dynamic analysis of structures[M]. Berkeley, California:Computers and Structures, Inc, 2002:12-13.

［8］杜修力. 水工建筑物抗震可靠度設(shè)計和分析用的隨機地震輸入模型[J].地震工程與工程振動,1998,18(14)：76-81.

DU Xiuli. A random seismic input model used in seismic reliability design and analysis of hydraulic buildings[J]. Earthquake Engineering and Engineering Vibration, 1998, 18(14): 76-81. (In Chinese)