兩維地震激勵(lì)下格構(gòu)拱動(dòng)力響應(yīng)簡(jiǎn)化評(píng)估方法

曲 揚(yáng)，羅永峰，朱釗辰，黃青隆

(同濟(jì)大學(xué) 土木工程學(xué)院， 上海 200092)

自汶川大地震以來(lái)，我國(guó)地震災(zāi)害的機(jī)理研究、監(jiān)測(cè)預(yù)警、結(jié)構(gòu)抗震設(shè)防和鑒定加固體系深度革新[1]，在此背景下，地震反應(yīng)分析對(duì)象從原有的一維、單體結(jié)構(gòu)逐漸向多維化、區(qū)域化的建筑群過(guò)渡，其中大跨度網(wǎng)格結(jié)構(gòu)作為最常見(jiàn)的公共建筑結(jié)構(gòu)形式，其多維性能化抗震設(shè)計(jì)的研究狀況滯后于工程應(yīng)用的現(xiàn)象越發(fā)明顯.目前，大跨度網(wǎng)格結(jié)構(gòu)在強(qiáng)震下彈塑性反應(yīng)多采用時(shí)程分析方法計(jì)算，但該方法依賴于地震波的選擇，且計(jì)算代價(jià)較高[2].鑒于此，脫胎于反應(yīng)譜理論發(fā)展而來(lái)的推覆分析方法，由于概念清晰、計(jì)算簡(jiǎn)便，常用于推定結(jié)構(gòu)在罕遇地震作用下的彈塑性行為，該方法計(jì)算精度能夠滿足抗震性能化設(shè)計(jì)需求[3].

源于經(jīng)典Pushover方法基本假定，推覆分析方法在具體實(shí)踐中需面臨荷載模式的確定這一關(guān)鍵問(wèn)題.侯愛(ài)波等[4]在對(duì)不同荷載模式的適用性進(jìn)行深入研究后認(rèn)為，對(duì)于振型豐富且復(fù)雜的結(jié)構(gòu)形式，應(yīng)采用考慮結(jié)構(gòu)多階振型的荷載模式進(jìn)行推覆分析.相陽(yáng)等[5]采用SRSS(平方和開(kāi)平方根)方法對(duì)振型位移進(jìn)行組合，給出了適用于網(wǎng)格結(jié)構(gòu)的推覆荷載模式，但該荷載模式僅考慮了結(jié)構(gòu)的自振特性，并未考慮地震動(dòng)輸入的特性.Ohsaki等[6]針對(duì)大跨度格構(gòu)拱結(jié)構(gòu)提出了一種多模態(tài)線性組合的推覆分析方法，雖然考慮了位移耦合和多階振型，但多次組合模態(tài)進(jìn)行推覆分析的方法，既缺乏物理意義解釋，也喪失了簡(jiǎn)潔性.Poursha等[7]基于對(duì)不同振型坐標(biāo)相關(guān)性的分析，給出了適用于多高層結(jié)構(gòu)的振型向量線性組合的荷載模式，但該方法需事先根據(jù)彈性譜計(jì)算振型坐標(biāo)的相關(guān)性，且針對(duì)網(wǎng)格結(jié)構(gòu)及多維地震激勵(lì)的適用性尚未被驗(yàn)證.

自Chopra提出模態(tài)推覆分析方法后，推覆分析的熱點(diǎn)逐漸由多振型問(wèn)題轉(zhuǎn)向多維化問(wèn)題，已有國(guó)內(nèi)外學(xué)者對(duì)此進(jìn)行了深入研究.Reyes等[8]基于模態(tài)推覆分析方法和主軸模型，率先提出了三維模態(tài)推覆分析方法，該方法對(duì)多振型響應(yīng)進(jìn)行CQC(完全二次項(xiàng)平方根)組合、對(duì)多維響應(yīng)進(jìn)行SRSS組合，從而得到結(jié)構(gòu)總響應(yīng).但該方法僅討論了結(jié)構(gòu)在水平兩維激勵(lì)下的動(dòng)力響應(yīng)，并未涉及對(duì)大跨度網(wǎng)格結(jié)構(gòu)尤為重要的豎向地震作用.為擺脫多維模態(tài)推覆分析方法對(duì)主軸模型的依賴，Manoukas等[9]引入了兩維地震動(dòng)輸入呈線性關(guān)系的基本假定，基于一個(gè)能夠同時(shí)考慮兩維激勵(lì)和響應(yīng)的等效單自由度體系，提出了兩維多模態(tài)推覆分析方法，通過(guò)一次模態(tài)推覆即可同時(shí)考慮兩維響應(yīng)，但該方法相比三維模態(tài)推覆分析方法實(shí)際上并未減少所需選取的振型階數(shù)，且計(jì)算精度并無(wú)顯著提升.王豐等[10]引入了兩維彈性反應(yīng)譜呈線性相關(guān)關(guān)系的假定，針對(duì)多高層結(jié)構(gòu)提出了兩階段多維推覆分析方法，通過(guò)改進(jìn)的需求譜求解目標(biāo)性能點(diǎn)，但該方法確定荷載模式時(shí)僅考慮了基底剪力，顯然不適用于多維耦合的網(wǎng)格結(jié)構(gòu).針對(duì)已有推覆分析成果尚未涉及網(wǎng)格結(jié)構(gòu)豎向地震作用的現(xiàn)狀，楊木旺[11]沿用Pushover方法的思路，利用基于延性系數(shù)的方法建立了彈塑性豎向需求譜，從而提出了評(píng)估網(wǎng)格結(jié)構(gòu)豎向抗震性能的Pushdown方法，但該方法仍無(wú)法考慮多維位移耦合效應(yīng).

針對(duì)上述問(wèn)題，本文通過(guò)對(duì)振型響應(yīng)進(jìn)行合理的線性組合，構(gòu)造了兩維地震激勵(lì)下的格構(gòu)拱結(jié)構(gòu)的變形模式，在對(duì)兩維地震動(dòng)輸入進(jìn)行適當(dāng)簡(jiǎn)化的基礎(chǔ)上，建立了基于給定變形模式的等效單自由度體系，并推導(dǎo)了其動(dòng)力方程，進(jìn)而給出了等效單自由度體系的A-D格式能力曲線，從而提出了大跨度格構(gòu)拱結(jié)構(gòu)在兩維地震激勵(lì)下動(dòng)力響應(yīng)的簡(jiǎn)化評(píng)估方法，并總結(jié)了簡(jiǎn)化方法的計(jì)算步驟.選取一個(gè)大跨度柱承格構(gòu)拱結(jié)構(gòu)為分析對(duì)象，對(duì)結(jié)構(gòu)分別在硬土和軟土場(chǎng)地若干條地震波兩維分量激勵(lì)下的動(dòng)力響應(yīng)進(jìn)行計(jì)算，并與時(shí)程法計(jì)算結(jié)果進(jìn)行對(duì)比，以驗(yàn)證簡(jiǎn)化方法的適用性與實(shí)用性.

1 簡(jiǎn)化評(píng)估方法基本理論

1.1 變形模式的構(gòu)造

根據(jù)振型疊加法，結(jié)構(gòu)體系的總響應(yīng)常采用CQC或SRSS振型響應(yīng)組合法進(jìn)行計(jì)算，但CQC方法過(guò)于復(fù)雜，實(shí)用性低，SRSS方法忽略振型間相互影響，且采用SRSS組合時(shí)，式中僅出現(xiàn)被不規(guī)則放大的正向位移響應(yīng)，而無(wú)負(fù)向位移響應(yīng)，造成變形模式進(jìn)一步失真，尤其不適用于振型復(fù)雜的大跨度網(wǎng)格結(jié)構(gòu).此外，基于振型疊加法的模態(tài)推覆分析方法，需對(duì)選取的每一階振型均進(jìn)行一次推覆分析，對(duì)于有效質(zhì)量參與系數(shù)累積緩慢的網(wǎng)格結(jié)構(gòu)，對(duì)多階振型進(jìn)行模態(tài)推覆分析在一定程度上犧牲了推覆分析方法的簡(jiǎn)潔性.基于此，本節(jié)給出線性組合振型響應(yīng)進(jìn)而構(gòu)造變形模式的方法，從而兼顧地震動(dòng)輸入特性和多階振型的影響，將多階振型推覆分析縮減為一次.

對(duì)于線彈性體系，結(jié)構(gòu)在地震作用下的動(dòng)力響應(yīng)u(t)可擴(kuò)展表示為

(1)

式中：φl(shuí)、ql(t)分別為第l階振型的振型向量和廣義振型坐標(biāo).

根據(jù)結(jié)構(gòu)動(dòng)力學(xué)，第l階振型的響應(yīng)ul(t)為

ul=ΓlDl(t)φl(shuí)

(2)

式中：Γl、Dl(t)分別為第l階振型的振型參與系數(shù)和譜位移.從而式(1)可寫(xiě)為

(3)

對(duì)于非線性體系，地震作用下的位移響應(yīng)仍可按線彈性體系的模態(tài)分解方法(式(3))進(jìn)行擴(kuò)展分解，但隨著結(jié)構(gòu)進(jìn)入彈塑性階段，由于剛度退化，振型已不同于線彈性階段，若仍按照式(3)形式采用線彈性階段的振型組合得到結(jié)構(gòu)位移響應(yīng)，可通過(guò)振型組合系數(shù)的形式考慮非線性體系的彈塑性發(fā)展以及各階振型的相關(guān)性(如文獻(xiàn)[5-7]).

有效質(zhì)量參與系數(shù)反映振型空間分布與動(dòng)力輸入空間分布間的相似程度，是衡量結(jié)構(gòu)振型反應(yīng)對(duì)結(jié)構(gòu)整體振動(dòng)反應(yīng)貢獻(xiàn)的合理指標(biāo).一般而言，振型的有效質(zhì)量參與系數(shù)越高，則該階振型對(duì)于結(jié)構(gòu)整體響應(yīng)的控制作用往往越顯著，其振型組合系數(shù)越高[12-13].因此，針對(duì)格構(gòu)拱結(jié)構(gòu)，為考慮橫向和豎向位移的耦合，同時(shí)兼顧地震動(dòng)輸入特性和多階振型的相互影響，本文以有效質(zhì)量參與系數(shù)作為振型組合系數(shù)，構(gòu)造兩維地震激勵(lì)下目標(biāo)性能點(diǎn)處的變形模式Φr為

Φr=∑rnΓnDnφn

(4)

式中：rn為第n階振型的有效質(zhì)量參與系數(shù).計(jì)算時(shí)，可依據(jù)有效質(zhì)量參與系數(shù)遴選結(jié)構(gòu)主振型，并使其累積值達(dá)到抗震規(guī)范的要求.

基于式(4)構(gòu)造的變形模式，僅需對(duì)結(jié)構(gòu)實(shí)施一次非線性推覆分析，計(jì)算結(jié)果即可較為全面地反映地震動(dòng)輸入特性及各主振型對(duì)結(jié)構(gòu)總響應(yīng)的貢獻(xiàn).

1.2 動(dòng)力方程的轉(zhuǎn)化

(5)

式中：u(t)為由于兩維地震激勵(lì)產(chǎn)生的結(jié)構(gòu)體系動(dòng)力位移向量；M和C分別為結(jié)構(gòu)體系的質(zhì)量矩陣和阻尼矩陣；lx和lz分別為x向和z向影響因子；Fs(t)為結(jié)構(gòu)體系恢復(fù)力向量，與當(dāng)前變形狀態(tài)和加載歷史相關(guān).

現(xiàn)假定兩維地震動(dòng)輸入呈相關(guān)關(guān)系，即

(6)

式中：κ為比例系數(shù).于是，式(5)右邊的有效地震力peff(t)可簡(jiǎn)化為

(7)

式中：lxz=lx+κlz為兩維等效影響因子.

假定結(jié)構(gòu)體系在兩維地震激勵(lì)下的變形模式為Φr，則動(dòng)力位移向量為

u(t)=Φrqr(t)

(8)

式中：qr(t)為對(duì)應(yīng)變形模式Φr的廣義坐標(biāo).

(9)

(10)

(11)

(12)

取等效振型參與系數(shù)Γr=Lr/Mr，令qr(t)=ΓrDr(t)，代入式(12)整理得到

(13)

式(13)即為與原MDF體系相對(duì)應(yīng)的等效單自由度(Equivalent Single Degree of Freedom, ESDF)體系的基本動(dòng)力方程，其中等效位移響應(yīng)Dr(t)與結(jié)構(gòu)整體動(dòng)力位移響應(yīng)u(t)滿足如下關(guān)系

u(t)=ΓrΦrDr(t)

(14)

(15)

1.3 兩維地震動(dòng)輸入的簡(jiǎn)化

文獻(xiàn)[9]和[10]均引入了兩維地震動(dòng)輸入呈線性相關(guān)關(guān)系的基本假定進(jìn)行研究，雖大大簡(jiǎn)化了理論推導(dǎo)過(guò)程，但其基本假定中線性相關(guān)系數(shù)的取值均未考慮地震動(dòng)輸入特性，即對(duì)于不同的地震動(dòng)，兩維輸入的線性相關(guān)關(guān)系不變，與實(shí)際不符.本節(jié)基于地震動(dòng)輸入、彈性反應(yīng)譜、結(jié)構(gòu)峰值位移響應(yīng)量三者之間的對(duì)應(yīng)關(guān)系，提出改進(jìn)的兩維地震動(dòng)輸入簡(jiǎn)化方法.

(19)

對(duì)于兩維地震動(dòng)輸入呈相關(guān)關(guān)系這一基本假定(式(6))，其關(guān)鍵是確定比例系數(shù)κ，結(jié)合式(17)和(19)可知：

(20)

(21)

1.4 ESDF體系的建立

本節(jié)采用非線性推覆分析轉(zhuǎn)化得到Ar-Dr關(guān)系，從而建立ESDF體系.

取結(jié)構(gòu)動(dòng)力變形模式為Φr，則結(jié)構(gòu)體系的位移響應(yīng)ur為

ur=ΓrDrΦr

(22)

設(shè)線彈性體系的恢復(fù)力增量ΔFs為

ΔFs=ΚΔur=ΓrΔDrΚΦr=

(23)

從而推覆荷載增量ΔFr可表示為

ΔFr=ΔFs=ΔχrMΦr

(24)

根據(jù)式(15)，偽加速度增量ΔAr和等效位移增量ΔDr分別為

(25)

(26)

進(jìn)入非線性階段，由于剛度退化，振型已不同于線彈性階段，但考慮到振型相互之間耦聯(lián)性仍較弱，推覆分析的推覆荷載增量ΔFr仍可采用式(24)表示，因而，此時(shí)結(jié)構(gòu)體系的恢復(fù)力增量ΔFs可以表示為

ΔFs=ΔFr=ΔχrMΦr

(27)

從而偽加速度增量ΔAr為

(28)

(29)

至此，通過(guò)對(duì)結(jié)構(gòu)進(jìn)行非線性推覆分析，得到了非線性體系的Ar-Dr關(guān)系.

1.5 簡(jiǎn)化方法計(jì)算步驟

本文簡(jiǎn)化方法的計(jì)算步驟如下：

(1) 對(duì)結(jié)構(gòu)進(jìn)行模態(tài)分析，選取主振型，根據(jù)式(4)構(gòu)造變形模式Φr；

(2) 針對(duì)選取的地震波兩維分量，按1.3節(jié)對(duì)兩維地震動(dòng)輸入進(jìn)行簡(jiǎn)化；

(3) 由變形模式導(dǎo)出荷載模式(式(24))，對(duì)結(jié)構(gòu)進(jìn)行非線性推覆分析，得到Ar-Dr關(guān)系；

(4) 基于Ar-Dr關(guān)系，建立對(duì)應(yīng)變形模式Φr的ESDF體系；

(5) 對(duì)ESDF體系進(jìn)行時(shí)程分析，得到其目標(biāo)性能點(diǎn)；

(6) 根據(jù)目標(biāo)性能點(diǎn)處的結(jié)構(gòu)彈塑性響應(yīng)，對(duì)結(jié)構(gòu)的抗震性能進(jìn)行評(píng)估.

2 數(shù)值算例

為了驗(yàn)證本文所提出簡(jiǎn)化方法的適用性，基于 ANSYS 和MATLAB平臺(tái)，選取一系列地震波兩維分量，采用本文簡(jiǎn)化方法和時(shí)程分析方法對(duì)一個(gè)大跨度柱承格構(gòu)拱結(jié)構(gòu)在兩維地震激勵(lì)下的動(dòng)力響應(yīng)進(jìn)行計(jì)算分析，并以兩向節(jié)點(diǎn)位移、單元最大應(yīng)力以及計(jì)算精度和耗時(shí)作為對(duì)比參數(shù)，分析簡(jiǎn)化方法的誤差和效率.

2.1 格構(gòu)拱結(jié)構(gòu)模型

大跨度柱承格構(gòu)拱結(jié)構(gòu)跨度45 m，高度15 m，上部拱結(jié)構(gòu)高度5 m，厚度1.2 m，矢跨比1/9，下部支承柱高度10 m，寬度1.5 m，支承柱腳為固定鉸支座.上部拱結(jié)構(gòu)外環(huán)節(jié)點(diǎn)集中質(zhì)量0.5 t，內(nèi)環(huán)節(jié)點(diǎn)集中質(zhì)量0.1 t，結(jié)構(gòu)模型和桿件截面如圖1所示.材料為Q235鋼，采用雙線性隨動(dòng)強(qiáng)化模型，彈性模量206 GPa，屈服強(qiáng)度235 MPa，屈服后彈性模量0.8 GPa，結(jié)構(gòu)阻尼比為0.02.上部拱結(jié)構(gòu)節(jié)點(diǎn)編號(hào)為順時(shí)針序，外環(huán)節(jié)點(diǎn)1～35號(hào)，內(nèi)環(huán)節(jié)點(diǎn)36～68號(hào)，下部支承柱節(jié)點(diǎn)自下而上、自左而右編號(hào)為 69～88號(hào)，如圖1中括號(hào)標(biāo)注所示.上部拱結(jié)構(gòu)單元編號(hào)同樣為順時(shí)針序，上弦桿單元1～34號(hào)，下弦桿單元35～66號(hào)，直腹桿單元67～99號(hào)，斜腹桿單元100～133號(hào)，下部支承柱單元如圖1中局部詳圖標(biāo)注所示.本算例研究格構(gòu)拱結(jié)構(gòu)在兩維地震激勵(lì)下的動(dòng)力響應(yīng)，x向和z向主振型信息見(jiàn)表1和圖2，所列兩向主振型的累積質(zhì)量參與系數(shù)均超過(guò)了70%.

圖1 格構(gòu)拱結(jié)構(gòu)模型及桿件截面Fig.1 Structural layout and member sections of the arch

圖2 格構(gòu)拱結(jié)構(gòu)主振型Fig.2 Fundamental vibration modes of the arch

表1 格構(gòu)拱結(jié)構(gòu)主振型信息Tab.1 Information of fundamental modes of the arch

2.2 結(jié)構(gòu)ESDF體系

為兼顧簡(jiǎn)化方法的準(zhǔn)確性和簡(jiǎn)潔性，本算例選取x向第1階和z向前2階主振型構(gòu)造變形模式，累積質(zhì)量參與系數(shù)達(dá)到70%，進(jìn)而根據(jù)式(24)得到荷載模式進(jìn)行非線性推覆分析，求得ESDF體系的Ar-Dr關(guān)系，對(duì)非線性推覆分析所得Ar-Dr曲線進(jìn)行雙線性化處理[14]，即得到可用于時(shí)程分析求解的ESDF體系A(chǔ)r-Dr關(guān)系曲線.以地震波FS1和SS1為例，ESDF體系推覆分析得到的雙線性化Ar-Dr關(guān)系曲線如圖3所示.

圖3 ESDF在地震波FS1和SS1作用下的Ar-Dr關(guān)系曲線Fig.3 Ar-Dr curves for ESDF under FS1 and SS1

圖4 地震波偽加速度反應(yīng)譜圖Fig.4 Pseudo acceleration response spectra of the seismic waves

2.3 結(jié)構(gòu)動(dòng)力響應(yīng)

2.3.1地震動(dòng)輸入 為驗(yàn)證簡(jiǎn)化方法對(duì)于不同場(chǎng)地周期的適用性，本文從日本地球科學(xué)與防災(zāi)技術(shù)研究所(NIED)的K-NET強(qiáng)震數(shù)據(jù)庫(kù)中選取2003年以來(lái)震級(jí)6.5級(jí)以上的8次地震，提取其硬土和軟土場(chǎng)地各12條地震波的x向和z向分量作為地震動(dòng)輸入，其偽加速度譜見(jiàn)圖4.根據(jù)主軸模型[15-16]，兩維反應(yīng)譜加速度Axz(T)為

(30)

2.3.2結(jié)構(gòu)節(jié)點(diǎn)位移響應(yīng) 采用簡(jiǎn)化方法(Simplified Procedure, SP)和時(shí)程分析(Response History Analysis, RHA)方法得到的格構(gòu)拱結(jié)構(gòu)節(jié)點(diǎn)位移(Unode)如圖5所示.相比時(shí)程分析方法計(jì)算結(jié)果，簡(jiǎn)化方法動(dòng)力響應(yīng)計(jì)算結(jié)果的平均誤差ρ采用下式計(jì)算：

(31)

式中：rSP和rRHA分別為采用簡(jiǎn)化方法和時(shí)程分析方法計(jì)算得到的響應(yīng)量.

圖5 簡(jiǎn)化方法和時(shí)程分析方法節(jié)點(diǎn)位移計(jì)算結(jié)果對(duì)比Fig.5 Nodal displacement responses given by SP and RHA

計(jì)算結(jié)果表明，采用簡(jiǎn)化方法計(jì)算得到的目標(biāo)性能點(diǎn)處的結(jié)構(gòu)位移響應(yīng)與時(shí)程分析方法的計(jì)算結(jié)果比較接近：① 兩者計(jì)算結(jié)果的位移值隨節(jié)點(diǎn)號(hào)的變化趨勢(shì)大致相同，表明兩者在目標(biāo)性能點(diǎn)處的變形模式基本一致，證明了本文構(gòu)造變形模式的合理性；② 與時(shí)程分析方法的計(jì)算結(jié)果相比，簡(jiǎn)化方法在x方向最大位移節(jié)點(diǎn)處位移計(jì)算值的平均誤差為15.1%，z方向最大位移節(jié)點(diǎn)處位移計(jì)算值的平均誤差為12.5%，且絕大部分節(jié)點(diǎn)誤差在25%以內(nèi)，說(shuō)明本文提出的兩維地震動(dòng)輸入簡(jiǎn)化方法具備一定的適用性，能夠有效估計(jì)結(jié)構(gòu)的地震需求，但對(duì)結(jié)構(gòu)體系目標(biāo)性能點(diǎn)的求解精度仍有待進(jìn)一步改進(jìn).此外，由圖5可見(jiàn)，對(duì)于z方向最大位移，簡(jiǎn)化方法計(jì)算結(jié)果均大于時(shí)程分析方法計(jì)算結(jié)果，說(shuō)明簡(jiǎn)化方法對(duì)于豎向位移的計(jì)算評(píng)估偏于安全.

圖6 簡(jiǎn)化方法和時(shí)程分析方法單元最大應(yīng)力云圖對(duì)比Fig.6 The extreme stresses of elements given by SP and RHA

圖7 簡(jiǎn)化方法計(jì)算誤差Fig.7 Calculating errors of SP

2.3.3結(jié)構(gòu)最大單元應(yīng)力響應(yīng) 采用簡(jiǎn)化方法和時(shí)程分析方法得到的格構(gòu)拱結(jié)構(gòu)單元最大應(yīng)力(σmax)云圖如圖6所示.由圖可見(jiàn)，采用簡(jiǎn)化方法求得的目標(biāo)性能點(diǎn)處的單元最大應(yīng)力與時(shí)程分析方法的計(jì)算結(jié)果比較接近，兩者計(jì)算結(jié)果的應(yīng)力分布和變化趨勢(shì)基本相同，表明兩者在目標(biāo)性能點(diǎn)處的單元應(yīng)力狀態(tài)基本一致，說(shuō)明基于變形模式導(dǎo)出的推覆荷載模式具備合理性.采用式(31)計(jì)算簡(jiǎn)化方法所得結(jié)構(gòu)最大單元應(yīng)力響應(yīng)的誤差，各單元應(yīng)力響應(yīng)計(jì)算結(jié)果的誤差繪于圖7(a)，其平均誤差為16.6%，大部分單元應(yīng)力誤差在30%以內(nèi).此外，簡(jiǎn)化方法計(jì)算得到的單元應(yīng)力普遍大于時(shí)程分析方法計(jì)算結(jié)果，說(shuō)明簡(jiǎn)化方法對(duì)于單元應(yīng)力響應(yīng)的計(jì)算評(píng)估偏于安全.將簡(jiǎn)化方法和時(shí)程分析方法預(yù)測(cè)得到的屈服桿件數(shù)量統(tǒng)計(jì)示于圖7(b)，可見(jiàn)簡(jiǎn)化方法對(duì)于結(jié)構(gòu)屈服桿件的預(yù)測(cè)可靠性較高，屈服桿件數(shù)量略高于時(shí)程分析方法，說(shuō)明預(yù)測(cè)結(jié)果偏于安全.

2.4 計(jì)算效率對(duì)比

與對(duì)格構(gòu)拱結(jié)構(gòu)直接進(jìn)行時(shí)程分析相比，采用簡(jiǎn)化方法，依據(jù)累積質(zhì)量參與系數(shù)截取不同振型階數(shù)構(gòu)造變形模式，從而分析得到的節(jié)點(diǎn)位移和單元最大應(yīng)力計(jì)算誤差和計(jì)算耗時(shí)占比如表2所示.綜合圖5～7和表2可以發(fā)現(xiàn)，在保證了一定的計(jì)算精度的前提下，簡(jiǎn)化方法相比時(shí)程分析方法在計(jì)算效率方面優(yōu)勢(shì)顯著，耗時(shí)僅為時(shí)程分析方法的20%～25%.此外，隨著選取主振型累積質(zhì)量參與系數(shù)的增大，簡(jiǎn)化方法對(duì)位移響應(yīng)和單元應(yīng)力響應(yīng)的計(jì)算精度分別提高9.1%和8.4%，而計(jì)算耗時(shí)也隨之增加，達(dá)到了時(shí)程分析方法的30%～35%，但相比之下仍然優(yōu)勢(shì)顯著，說(shuō)明隨著截取振型的增多，本文構(gòu)造的變形模式考慮了更高階振型，故計(jì)算結(jié)果有所提高.綜上，當(dāng)累積質(zhì)量參與系數(shù)達(dá)到70%時(shí)，簡(jiǎn)化方法可較好地兼顧精度和耗時(shí)，計(jì)算效率較高.

表2 簡(jiǎn)化方法與時(shí)程分析方法計(jì)算效率對(duì)比Tab.2 Comparison of computing efficiency by SP and RHA

3 結(jié)論

本文通過(guò)對(duì)各階振型響應(yīng)按有效質(zhì)量參與系數(shù)進(jìn)行線性組合，給出了兩維地震激勵(lì)下的格構(gòu)拱結(jié)構(gòu)的變形模式，并在對(duì)兩維地震動(dòng)輸入進(jìn)行適當(dāng)簡(jiǎn)化的基礎(chǔ)上，提出了大跨度格構(gòu)拱結(jié)構(gòu)在兩維地震激勵(lì)下動(dòng)力響應(yīng)的簡(jiǎn)化評(píng)估方法，同時(shí)給出了簡(jiǎn)化方法的理論推導(dǎo)過(guò)程和計(jì)算步驟，通過(guò)兩種場(chǎng)地條件下兩維地震激勵(lì)的數(shù)值算例與時(shí)程分析方法進(jìn)行對(duì)比分析，得出以下結(jié)論：

(1) 本文提出的變形模式能夠考慮雙向位移響應(yīng)，且避免了SRSS方法對(duì)網(wǎng)格結(jié)構(gòu)位移響應(yīng)的不合理放大，并將推覆分析縮減為一次；

(2)研究表明，本文對(duì)兩維地震動(dòng)輸入的簡(jiǎn)化，在保證了一定的計(jì)算精度的基礎(chǔ)上考慮了地震動(dòng)輸入特性，從而大大簡(jiǎn)化了理論推導(dǎo)；

(3) 數(shù)值分析結(jié)果表明，無(wú)論是硬土還是軟土場(chǎng)地，簡(jiǎn)化方法對(duì)大跨度格構(gòu)拱結(jié)構(gòu)在兩維地震激勵(lì)下的動(dòng)力響應(yīng)均具有較好的適用性，節(jié)點(diǎn)兩向位移平均誤差分別為15.1%和12.5%，單元最大應(yīng)力平均誤差為16.6%，計(jì)算精度較高；

(4) 簡(jiǎn)化方法概念簡(jiǎn)潔、可操作性強(qiáng)，且僅需較小的計(jì)算代價(jià)，即可得到較高的精度，計(jì)算結(jié)果偏于安全，可用于初步評(píng)估大跨度格構(gòu)拱結(jié)構(gòu)在兩維地震激勵(lì)下的抗震性能.