吳 邊 張 敏

(四川建筑職業(yè)技術(shù)學(xué)院，四川 德陽 618000)

?

·結(jié)構(gòu)·抗震·

大跨超長類結(jié)構(gòu)多點(diǎn)激勵輸入地震動應(yīng)用研究★

吳 邊 張 敏

(四川建筑職業(yè)技術(shù)學(xué)院，四川 德陽 618000)

圍繞合理確定大跨、超長類結(jié)構(gòu)的多點(diǎn)激勵輸入地震動這一問題，對當(dāng)前在應(yīng)用相干函數(shù)模型及多點(diǎn)地震動合成中遇到的一些問題作了研究，分析了相干函數(shù)模型的適用性，并給出了H-V模型和Qu模型在相干程度下的參數(shù)取值，基于Hao方法,歸納了空間地震動的人工合成方法。

多點(diǎn)激勵，相干函數(shù)模型，地震動合成，大跨結(jié)構(gòu)

0 引言

近年來我國建成了大量的大跨，超長類結(jié)構(gòu)工程，該類結(jié)構(gòu)大多建設(shè)在地震區(qū)，有的甚至是在高烈度區(qū)，結(jié)構(gòu)的抗震問題突出；而且，該類結(jié)構(gòu)往往擔(dān)負(fù)著重要功能，如大跨度橋梁、航站樓、輸電塔線等，若出現(xiàn)地震后結(jié)構(gòu)失效會造成不可估量的損失；因此，對大跨、超長類結(jié)構(gòu)的抗震設(shè)計就顯得尤為重要[1]。

由于大跨、超長類結(jié)構(gòu)屬于空間延展結(jié)構(gòu)，結(jié)構(gòu)所處區(qū)域范圍廣，場地條件變化較大。在對該類結(jié)構(gòu)進(jìn)行動力時程分析時，若仍采用一致地震動輸入會產(chǎn)生較大誤差，更為合理的方法是采用能在一定程度上反映地震動的空間變化特征的多點(diǎn)激勵地震動輸入。當(dāng)前對行波效應(yīng)考慮較成熟，也有規(guī)范可循；然而在實際工程中對相干效應(yīng)和局部場地效應(yīng)的考慮很少。本文主要就針對如何確定合理的多點(diǎn)激勵輸入地震動加以研究，并給出了一套相應(yīng)的方法，能夠為大跨、超長類結(jié)構(gòu)在確定輸入地震動的問題上提供參考。

1 多點(diǎn)激勵輸入地震動應(yīng)用現(xiàn)狀

目前對大跨、超長類結(jié)構(gòu)考慮空間效應(yīng)的地震響應(yīng)分析方法有修正反應(yīng)譜法，隨機(jī)振動分析方法，時程分析法。由于前兩類分析方法在考慮地震動空間效應(yīng)時有較大的局限性，而時程分析方法計算原理成熟，商用有限元分析軟件(如SAP2000，MIDAS等)中可直接使用，在實際工程中應(yīng)用廣泛。為了在時程分析中考慮地震動的空間變化，逐漸發(fā)展出一套多點(diǎn)激勵的時程分析方法；多點(diǎn)激勵即在結(jié)構(gòu)的各個支承點(diǎn)輸入具有空間相關(guān)性的一組地震波，而非同一條地震波。該類方法較傳統(tǒng)的一致輸入的分析方法，可在一定程度內(nèi)考慮地震動的空間變化，其考慮程度取決于采用的地震動空間變化模型。

從當(dāng)前我國在大跨類工程的實踐來看，在選擇輸入地震動時，若要完整地考慮地震動的空間變化全部因素，還不成熟；其中對相干函數(shù)模型的選擇問題尤為明顯。其中一個突出的矛盾是抗震分析需要計算多組“相干程度”工況，然而，除極個別研究者[6]外，絕大多數(shù)相干函數(shù)模型沒有給出多工況下相應(yīng)的參數(shù)取值。

為解決以上問題，本文以地震動臺陣記錄為依據(jù)研究了多點(diǎn)激勵輸入地震動合成中的有關(guān)問題，并給出了一套合成方法(見圖1)，最后通過算例加以驗證。

2 多點(diǎn)激勵輸入地震動合成

大跨、超長類結(jié)構(gòu)時程分析的第一步就是確定合理的多點(diǎn)輸入地震動。雖然是采用真實的差動臺陣記錄較為可靠，然而在工程應(yīng)用很難實現(xiàn)。對多點(diǎn)輸入地震動而言，最佳的解決辦法是人工合成多點(diǎn)地震動的方法。

2.1 自譜密度估計

已知實際地震動記錄，可采用Welch修正周期圖法得出自譜密度。若未給出實際地震動，則基于過濾白噪聲過程，采用功率譜(自譜密度)模型，當(dāng)前較多采用的是改進(jìn)K-T功率譜模型(杜修力、陳厚群[7])。

2.2 相干函數(shù)模型分析

考慮到當(dāng)前提出的相干函數(shù)模型較多(見表1)，就存在相干函數(shù)模型的選擇性問題。本文首先針對實際臺陣記錄，對不同的相干函數(shù)模型參數(shù)加以擬合，根據(jù)擬合效果(均方差、殘差平方、擬合優(yōu)度系數(shù))判定各模型的通用性，最后選取了2種相干函數(shù)模型。

表1 相干函數(shù)模型統(tǒng)計

1)Harichandran-Vanmarcke模型(以下稱H-V Model)。

參數(shù)取值：A=0.736;k=766;f0=1.09;b=2.78。

2)屈鐵軍模型(以下稱QuModel)。

|γ(ξ，ω)|=exp[-a(ω)ξb(ω)]。

其中，a(ω)=a1ω2+a2；b(ω)=b1ω+b2。

參數(shù)取值：a1=0.000 016 78;a2=0.001 219;b1=-0.005 5;b2=0.767 4。

2.2.1 模型參數(shù)擬合

以SMART1臺陣在1986年11月14日記錄的Event45為例(見圖2)，給出模型參數(shù)的擬合過程。

任意選取某一參考臺站(本文選擇I06)其他非參考臺數(shù)據(jù)站與其共同構(gòu)成樣本點(diǎn)。選擇方法：臺站基本處于一個方向；各個臺站間的距離應(yīng)該差別明顯，因此不宜在相似的距離選擇多個臺站；總的選擇數(shù)量不能太多，選擇多了會使擬合結(jié)果出現(xiàn)較大的離散性。根據(jù)上述原則，確定了4組臺站擬合樣本(I06～I(xiàn)07，I06～C00，I06～I(xiàn)12，I06～M06)，樣本點(diǎn)取平滑過后的值。參數(shù)估計采用UniversalGlobalOptimization算法，回歸計算采用Levenberg-Marquardt算法，達(dá)到收斂標(biāo)準(zhǔn)后，得出最終參數(shù)取值。

2.2.2 模型參數(shù)統(tǒng)計分析

借鑒抗震規(guī)范中統(tǒng)計設(shè)計反應(yīng)譜的思路，采用兩次擬合的方法進(jìn)行參數(shù)統(tǒng)計分析(總共統(tǒng)計了1 254組地震動記錄)。表2給出了H-VModel和QuModel模型3組代表三種“相干度”的參數(shù)取值建議。

表2 相干函數(shù)模型參數(shù)取值建議

2.3 空間地震動合成

空間地震動的時空序列描述可以用譜密度矩陣來表示：

其中，[S(iω)]為對稱正定陣Hermite矩陣，對[S(iω)]進(jìn)行Cholesky分解得到：

[S(iω)]=[L(iω)][LH(iω)]。

則有：

其中，第n點(diǎn)的地震動時程可表示為n項三角級數(shù)之和(Hao[14])；

(1)

Anm(ωk)和θnm(ωk)與譜密度矩陣之關(guān)系，可表達(dá)為(屈鐵軍等[12])：

由于式(1)無法考慮和后面生成點(diǎn)的相關(guān)性，將式(1)寫成離散的形式：

(2)

其中，l為時間;m為地點(diǎn);k為頻率。對式(2)進(jìn)行離散Fourier變換：

(3)

在計算機(jī)編程時，通過式(3)先生成Fourierspectrum，然后使用InverseFourierTransform得到最終的地震動時程。

3 算例分析

根據(jù)上節(jié)所述方法，在Matlab中編制計算程序MulPoint，根據(jù)不同的合成初始條件來合成(圖3中1為已知天然波，2為無天然波，3為已知功率譜模型)。

1)條件1。初始地震動選擇SMART1臺陣的Event45事件中C00臺站記錄到的EW水平地震波分量。假定視波速500m/s。采用H-V模型(中相干程度)。分段強(qiáng)度包線模型：

2)條件2。Ⅰ類場地，8度設(shè)防烈度，設(shè)計地震分組第3組，結(jié)構(gòu)的阻尼比ξ=0.05。杜修力—陳厚群模型。其他同條件1。

圖4,圖5分別給出了合成地震動的功率譜和相干函數(shù)值，從圖中可以看出，無論是采用哪一種初始條件，從地震動的功率譜來看，S1,S2，S3差異很小，與初始條件(C00或功率譜模型)差異也不大，是比較合理的；從吻合相干函數(shù)模型的情況來看，也是吻合很好的；總的來說，最終合成效果理想。

4 結(jié)語

圍繞合理確定大跨、超長類結(jié)構(gòu)的多點(diǎn)激勵輸入地震動這一問題。對相干函數(shù)模型的應(yīng)用以及多點(diǎn)地震動合成中遇到的一些問題加以研究。主要工作和結(jié)論可歸納為：

1)在實際地震動差動臺陣記錄的擬合結(jié)果的基礎(chǔ)上，對比了當(dāng)前較普遍的相干函數(shù)模型，給出了其中的H-V Model和Qu Model的考慮不同相干度的參數(shù)取值建議。

2)基于不同的初始條件(已知基準(zhǔn)地震動、功率譜模型)給出了多點(diǎn)激勵地震動的人工合成方法，并通過算例加以初步驗證。

[1] 李英民,劉立平.工程結(jié)構(gòu)的設(shè)計地震動[M].北京：科學(xué)出版社,2011.

[2] Newmark N M, Rosenblueth E. Fundamentals of earthquake engineering[J]. Civil engineering and engineering mechanics series,1971(12)：76.

[3] 趙基達(dá),藍(lán) 天.中國空間結(jié)構(gòu)三十年的進(jìn)展及今后展望[J].工業(yè)建筑,2013(4):131-138.

[4] 全 偉.大跨橋梁多維多點(diǎn)地震反應(yīng)分析研究[D].大連:大連理工大學(xué),2008.

[5] 白鳳龍.空間變化地震動激勵下大跨度結(jié)構(gòu)的反應(yīng)研究[D].大連:大連理工大學(xué),2010.

[6] Kaiming B,Hong H,Nawawi C. Required separation distance between decks and at butments of a bridge crossing a canyon site to avoid seismic pounding[J]. Earthquake Engineering Mechanics,1986,112(2):154-174.

[7] 杜修力,陳厚群.地震動隨機(jī)模擬及其參數(shù)確定方法[J].地震工程與工程振動,1994(4):1-5.

[8] Loh C H. Analysis of the spatial variation of seismic waves and ground movements from smart-1 array data[J]. Earthquake engineering & structural dynamics,1985,13(5):561-581.

[9] Harichandran R S, Vanmarcke E H. Stochastic variation of earthquake ground motion in space and time[J]. Journal of Engineering and Structural Dynamics,2009,39(3):303-329.

[10] Abrahamson.N.A,Schneider.J.F,Stepp.J.C Empirical spatial coherency functions for application to soil-structure interaction analyses[J]. Earthquake spectra,1991,7(1):1-27.

[11] 馮啟民,胡聿賢.空間相關(guān)地面運(yùn)動的數(shù)學(xué)模型[J].地震工程與工程振動,1981,1(2):1-8.

[12] 屈鐵軍,王君杰,王前信.空間變化的地震動功率譜的實用模型[J].地震學(xué)報,1996(1):55-62.

[13] 王君杰,陳 虎.面向設(shè)計應(yīng)用的地震動空間相干函數(shù)模型[J].地震工程與工程振動,2007(1):16-23.

[14] Hao H, Oliveira C S, Penzien J. Multiple-station ground motion processing and simulation based on SMART-1 array data[J]. Nuclear Engineering and Design,1989,111(3):293-310.

[15] Menke W, Lerner-Lam A L, Dubendorff B, et al. Polarization and coherence of 5 to 30 Hz seismic wave fields at a hard-rock site and their relevance to velocity heterogeneities in the crust[J]. Bulletin of the Seismological Society of America,1990,80(2):430-449．

[16] Luco J E, Wong H L. Response of a rigid foundation to a spatially random ground motion[J]. Earthquake engineering & structural dynamics,1986,14(6):891-908.

[17] Somerville P G, Mclaren J P, Saikia C K, et al. Site-specific estimation of spatial incoherence of strong ground motion[C]. ASCE,1988.

Applied study on inputted multi-point ground motions of long-span structures★

Wu Bian Zhang Min

(Sichuan College of Architectural Technology, Deyang 618000, China)

On the issue of identifying the multi-point ground motions of the large-span and super-long structures, the paper researches some problems in the application of the coherence function model and the integration of multi-point, analyzes the adaptability of the coherence functions, undertakes the parameter value with H-V model and Qu model under the coherence, and sums up the manual integration method of the space motion.

multi-point motion, coherence function, motion integration, large-span structure

1009-6825(2017)08-0023-03

2017-01-04★：四川建筑職業(yè)技術(shù)學(xué)院科研項目(2016KJ03)

吳 邊(1987- )，男，碩士，助教; 張 敏(1989- )，女，碩士，助教

TU318.1

A