李創(chuàng)第， 楊雪峰， 李宇翔， 葛新廣

(廣西科技大學(xué) 土木建筑工程學(xué)院，柳州 545006)

1 引 言

地震災(zāi)害調(diào)查研究表明[1-3]，相鄰建筑結(jié)構(gòu)在地震激勵(lì)作用下發(fā)生碰撞導(dǎo)致結(jié)構(gòu)構(gòu)件的破壞甚至整體倒塌的問(wèn)題日益凸顯，引發(fā)了專家學(xué)者的廣泛關(guān)注[4-6]。Wu等[7]對(duì)相鄰建筑間設(shè)置粘滯流體阻尼器的組合體系動(dòng)力性能進(jìn)行了研究，研究表明相鄰結(jié)構(gòu)間設(shè)置粘滯流體阻尼器具有良好的減震效果。劉佩等[8]對(duì)考慮碰撞效應(yīng)的相鄰框架結(jié)構(gòu)易損性進(jìn)行了研究，研究表明碰撞對(duì)結(jié)構(gòu)的局部響應(yīng)影響較大。劉紹峰等[9]對(duì)連接粘滯阻尼器的相鄰結(jié)構(gòu)進(jìn)行了振動(dòng)臺(tái)試驗(yàn)研究，研究表明連接阻尼器對(duì)相鄰結(jié)構(gòu)地震反應(yīng)具有一定的控制作用。Hao等[10]對(duì)相鄰結(jié)構(gòu)間連接調(diào)諧粘滯質(zhì)量阻尼器進(jìn)行了研究，發(fā)現(xiàn)調(diào)諧粘滯質(zhì)量阻尼器具有良好的減震性能。在抗震設(shè)計(jì)中層間位移大小決定了結(jié)構(gòu)內(nèi)力的大小，以上研究未能給出層間位移的表達(dá)式，因此，對(duì)結(jié)構(gòu)層間位移的研究具有重要的工程意義。

地震動(dòng)觀測(cè)表明，地震動(dòng)是一種隨機(jī)過(guò)程。工程界已提出了多種隨機(jī)激勵(lì)模型，如Kanai-Tajimi譜[11,12]、Clough-Penzien譜[13]和歐進(jìn)萍譜[14]等。Kanai-Tajimi譜將場(chǎng)地土表層的震動(dòng)看做場(chǎng)地土對(duì)下臥層基巖的震動(dòng)過(guò)濾而來(lái)，而基巖具有白噪聲隨機(jī)過(guò)程的震動(dòng)，是其他更復(fù)雜隨機(jī)地震動(dòng)模型的基礎(chǔ)。結(jié)構(gòu)隨機(jī)地震動(dòng)響應(yīng)分析方法主要有時(shí)域法和頻域法[15,16]。時(shí)域法應(yīng)用的前提[17,18]是隨機(jī)激勵(lì)模型具有協(xié)方差函數(shù)，Kanai-Tajimi雖然具有協(xié)方差表達(dá)式，但是獲得的地震動(dòng)響應(yīng)表達(dá)式復(fù)雜。而頻域法則是利用結(jié)構(gòu)響應(yīng)功率譜密度函數(shù)等于結(jié)構(gòu)的頻響函數(shù)模值的平方與地震動(dòng)激勵(lì)功率譜的代數(shù)乘積，該方法的不足之處是工程應(yīng)用時(shí)需先得到結(jié)構(gòu)的響應(yīng)方差和譜矩，而要獲得這些物理量需要對(duì)響應(yīng)功率譜密度函數(shù)進(jìn)行積分處理，其精度和效率受積分區(qū)間與積分步長(zhǎng)的選取影響較大[19]。XU等[20]研究了在相鄰結(jié)構(gòu)設(shè)置主動(dòng)控制下的組合體系基于Kanai-Tajimi譜隨機(jī)激勵(lì)下的封閉解，但表達(dá)式較為復(fù)雜。葛新廣等[17,18]利用Kanai-Tajimi譜的濾波方程與結(jié)構(gòu)的地震動(dòng)方程重構(gòu)，從而將結(jié)構(gòu)的地震動(dòng)方程轉(zhuǎn)化為易于獲得封閉解的白噪聲激勵(lì)，從時(shí)域法與頻域法聯(lián)合應(yīng)用的角度獲得了結(jié)構(gòu)位移及層間位移的0階～2階譜矩的簡(jiǎn)明封閉解。

本文對(duì)Kelvin型粘彈性阻尼器連接的相鄰結(jié)構(gòu)基于Kanai-Tajimi譜地震動(dòng)激勵(lì)下的結(jié)構(gòu)系列響應(yīng)方差和0階～2階譜矩進(jìn)行了研究，提出了一種簡(jiǎn)明解法。首先，利用Kanai-Tajimi譜的濾波方程與組合體結(jié)構(gòu)的運(yùn)動(dòng)方程聯(lián)立，將復(fù)雜的隨機(jī)地震動(dòng)模型精確轉(zhuǎn)化為白噪聲激勵(lì)模型；其次，利用復(fù)模態(tài)法獲得了組合體結(jié)構(gòu)位移及層間位移等系列響應(yīng)時(shí)域解的統(tǒng)一簡(jiǎn)明表達(dá)式；然后，基于時(shí)域法與頻域法聯(lián)合應(yīng)用獲得了組合體結(jié)構(gòu)相對(duì)于地面的絕對(duì)位移和層間相對(duì)位移的0階～2階譜矩的簡(jiǎn)明封閉解。

2 組合體結(jié)構(gòu)地震動(dòng)方程的重構(gòu)

兩相鄰結(jié)構(gòu)通過(guò)Kelvin型粘彈性阻尼器進(jìn)行連接形成組合體結(jié)構(gòu)，結(jié)構(gòu)計(jì)算簡(jiǎn)圖和Kelvin型粘彈性阻尼器本構(gòu)模型分別如圖1和圖2所示，其左側(cè)及右側(cè)結(jié)構(gòu)的運(yùn)動(dòng)方程分別為

(1)

(2)

圖1 結(jié)構(gòu)計(jì)算簡(jiǎn)圖

圖2 Kelvin型粘彈性阻尼器本構(gòu)模型

式中O1L為JL-1個(gè)元素均為0的行向量，O2 L為nL-JL個(gè)元素均為0的行向量，O1 R為JR-1個(gè)元素均為0的行向量，O2 R為nR-JR個(gè)元素均為0的行向量，JL和JR分別為左右側(cè)結(jié)構(gòu)連接樓層號(hào)，T為向量轉(zhuǎn)置。

將式(1，2)聯(lián)立，則組合體系的地震動(dòng)方程為

(3)

式中

式中o1為nL×nR階矩陣，內(nèi)部各元素均為0，組合體中設(shè)置阻尼器的左側(cè)與右側(cè)樓層的剛度及阻尼系數(shù)分別為kL,JLJL=kL,JL+kL,JL+ 1+kd，kR,JRJR=kR,JR+kR,JR+ 1+kd，其中,kL,JLJL和kR,JRJR分別為左右側(cè)結(jié)構(gòu)剛度矩陣的對(duì)角線元素；cL,JLJL=cL,JL+cL,JL+1+cd，cR,JRJR=cR,JR+cR,JR+ 1+cd，其中，cL,JLJL和cR,JRJR分別為左右側(cè)結(jié)構(gòu)阻尼矩陣中的對(duì)角線元素。組合體中未設(shè)置阻尼器的左側(cè)與右側(cè)樓層的剛度及阻尼系數(shù)分別為kL,i i=kL,i+kL,i + 1，kR,i i=kR,i+kR,i + 1；cL,i i=cL,i+cL,i + 1，cR,i i=cR,i+cR,i + 1。Cd和Kd為輔助矩陣，其中，Cd(JL,JR)=-cd，Kd(JL,JR)=-kd，其余元素均為0。

Kanai-Tajimi模型的濾波方程為[11,12,17,18]

(4a,4b)

CüR(τ)=2πS0δ(τ)

(5)

式中S0為地震動(dòng)強(qiáng)度常數(shù),δ(τ)為Dirac函數(shù)。

聯(lián)立式(3，4)，并用狀態(tài)方程表示

(6,7)

q=2nL+2nR+2

式中o2為(nL+nR)個(gè)元素為0的行向量，o3為元素為0的(nL+nR)階方陣，E為(nL+nR)階單位對(duì)角陣。

3 地震動(dòng)方程的復(fù)模態(tài)解耦

根據(jù)復(fù)模態(tài)理論[15,17-19]，存在特征值矩陣P、右特征向量矩陣U和左特征向量矩陣V，對(duì)式(6)復(fù)模態(tài)解耦，即存在

(8)

式中P為對(duì)角陣。

利用復(fù)模態(tài)變換

y=Uz

(9)

式中z為復(fù)模態(tài)變量。

將式(9)代入式(6)，利用復(fù)模態(tài)理論，式(6)可改寫(xiě)為

(10)

式中

(11)

由于P為對(duì)角陣，故式(10)的分量表達(dá)式為

(j=1,2,…,q)(12)

(j=1,2,…,q)(13)

4 組合體結(jié)構(gòu)地震動(dòng)響應(yīng)的統(tǒng)一解

對(duì)組合體結(jié)構(gòu)各層相對(duì)于地面的地震動(dòng)位移和層間地震動(dòng)變形的分析是工程結(jié)構(gòu)抗震設(shè)計(jì)的重要環(huán)節(jié)；而結(jié)構(gòu)層相對(duì)于地面的地震動(dòng)速度(絕對(duì)速度)和層間速度是結(jié)構(gòu)抗震動(dòng)力可靠度分析的基礎(chǔ)。由此，本文研究了上述參數(shù)的統(tǒng)一解。

4.1 組合體結(jié)構(gòu)地震動(dòng)響應(yīng)的絕對(duì)速度與位移

(14a)

(14b)

式中j=1～nL+nR(當(dāng)j=1～nL時(shí)為組合體系中左側(cè)的結(jié)構(gòu)，當(dāng)j=nL+1～nL+nR時(shí)為組合體系中右側(cè)的結(jié)構(gòu))，uj為右特征向量矩陣U的第j行向量；動(dòng)力響應(yīng)的模態(tài)強(qiáng)度系數(shù)λj,i為

(15)

式中uj,i為uj的第i個(gè)分量。

4.2 組合體結(jié)構(gòu)的層間位移與層間速度

(16a)

ΔxL,j=(uj-uj - 1)z=

(16b)

(16c)

式中j=2～nL。

(17)

(18a)

ΔxR,j=(unL+ j-unL+ j - 1)z=

(18b)

(18c)

式中j=2～nR。

(19)

由式(10，16，18)可得，結(jié)構(gòu)各層位移及其速度、層間位移及其速度可統(tǒng)一表示為

(20)

式中X(t)為地震動(dòng)響應(yīng)量;κi為響應(yīng)量對(duì)應(yīng)的模態(tài)強(qiáng)度系數(shù)，不同的響應(yīng)量分別見(jiàn)式(14，16～19)；Xi(t)分量形式如下,

(i=1～q)(21)

5 地震動(dòng)系列響應(yīng)的簡(jiǎn)明封閉解

5.1 組合體結(jié)構(gòu)地震動(dòng)響應(yīng)的協(xié)方差

由隨機(jī)振動(dòng)理論及式(20)，結(jié)構(gòu)地震動(dòng)響應(yīng)X的協(xié)方差為

CX(τ)=E[X(t)X(t+τ)]=

(22)

式中E[·]為數(shù)學(xué)期望。

由式(21)可知結(jié)構(gòu)地震動(dòng)響應(yīng)的分量的協(xié)方差為

E[Xk(t)Xi(t+τ)]=κkκi×

(23)

式中u和v為積分變量。

把式(5)代入式(23)，得

E[Xk(t)Xi(t+τ)]=2πS0κkκi×

(24)

利用Dirac函數(shù)的性質(zhì)，式(24)簡(jiǎn)化為

(25)

對(duì)式(25)積分可得

(26)

由式(22,26)，組合體結(jié)構(gòu)基于Kanai-Tajimi譜的動(dòng)力響應(yīng)的協(xié)方差為

(27)

由式(27)可知，組合體結(jié)構(gòu)基于Kanai-Tajimi譜隨機(jī)激勵(lì)的動(dòng)力響應(yīng)的協(xié)方差可用耗能體系的震動(dòng)復(fù)特征值函數(shù)的線性組合表示，表達(dá)式簡(jiǎn)潔明了，且為模態(tài)的完全組合解。

(28)

由隨機(jī)振動(dòng)理論，結(jié)構(gòu)動(dòng)力響應(yīng)的單邊功率譜與協(xié)方差存在Wiener-Khinchin關(guān)系[15]，即

(29)

式中SX(ω)為結(jié)構(gòu)響應(yīng)X的功率譜。

把式(27)代入式(29)并積分，可得

(30)

5.2 組合體結(jié)構(gòu)地震動(dòng)響應(yīng)0階～2階譜矩的簡(jiǎn)明封閉解

由平穩(wěn)激勵(lì)的譜矩[15]定義，地震動(dòng)響應(yīng)的0階譜矩αX，0為

(31)

對(duì)式(31)進(jìn)行積分運(yùn)算，可得

(32)

由隨機(jī)振動(dòng)理論，結(jié)構(gòu)隨機(jī)地震動(dòng)響應(yīng)的0階譜矩等于其方差，比較式(28，32)，可知本文方法計(jì)算0階譜矩的正確性。

由隨機(jī)振動(dòng)理論[15,19]，平穩(wěn)地震動(dòng)響應(yīng)變化率的0階譜矩等于其平穩(wěn)響應(yīng)的2階譜矩。因此，組合體系響應(yīng)的2階譜矩可由對(duì)應(yīng)響應(yīng)量的變化率的0階譜矩表示，即

(33)

組合體系響應(yīng)的一階譜矩[15]為

(34)

由文獻(xiàn)[17,18]，譜矩αX,1可表示為

(35)

根據(jù)式(14～19，32，33，35)，組合體結(jié)構(gòu)的地震動(dòng)響應(yīng)的0階～2階譜矩均有簡(jiǎn)明封閉解。

6 算 例

一鋼筋混凝土組合體結(jié)構(gòu)，左側(cè)結(jié)構(gòu)(主結(jié)構(gòu))樓層數(shù)為10層，右側(cè)結(jié)構(gòu)(相鄰結(jié)構(gòu))為6層，左右側(cè)結(jié)構(gòu)層高均為3.6 m，兩建筑結(jié)構(gòu)均在6層設(shè)置粘彈性阻尼器相連接。左側(cè)結(jié)構(gòu)各層質(zhì)量均為 1.5×105kg，各層抗側(cè)剛度均為1.67×107N/m；右側(cè)結(jié)構(gòu)各層質(zhì)量均為0.9×105kg，各層抗側(cè)剛度為2.04×107N/m；結(jié)構(gòu)阻尼采用瑞雷阻尼模型，阻尼比0.05。粘彈性阻尼器采用Kelvin模型，其力學(xué)參數(shù)為kd=1.0×106N/m，cd=1.0×105N·s/m。Kanai-Tajimi隨機(jī)激勵(lì)模型的參數(shù)為ωg=25.13 rad/s，ξg=0.64，S0=23.17×10-4m2/s3。

6.1 功率譜對(duì)比分析

由式(3)可知，基于虛擬激勵(lì)的結(jié)構(gòu)響應(yīng)的頻域解為

(36)

式中

x(ω)={xL,1(ω),…,xL,nL(ω),xR,1(ω),…,xR,nR(ω)}T

結(jié)構(gòu)層間位移的頻響域解Δx(ω)可以表示為

Δx(ω)={xL,1(ω),xL,2(ω)-xL,1(ω),…,

xL,nL(ω)-xL,nL- 1(ω),xR,1(ω),

xR,2(ω)-xR,1(ω),…,

xR,nL(ω)-xR,nL- 1(ω)}

則基于虛擬激勵(lì)法的結(jié)構(gòu)響應(yīng)功率譜為

(37)

為驗(yàn)證本文方法的正確性，分別將本文方法獲得功率譜密度函數(shù)與Kanai-Tajimi譜功率譜密度函數(shù)及虛擬激勵(lì)法獲得的結(jié)構(gòu)位移及結(jié)構(gòu)層間位移功率譜密度函數(shù)進(jìn)行了對(duì)比。

從圖3和圖4可以看出，兩種方法獲得的功率譜密度函數(shù)值完全重合，驗(yàn)證了本文方法計(jì)算功率譜密度函數(shù)的正確性。虛擬激勵(lì)法的功率譜密度函數(shù)為式(37)，每給定圓頻率變量ω，均需要進(jìn)行矩陣運(yùn)算，而本文方法的功率譜密度函數(shù)為式(30)，僅需要在求解組合體系運(yùn)動(dòng)方程的特征值及響應(yīng)系數(shù)時(shí)進(jìn)行一次矩陣運(yùn)算。

圖3 左側(cè)10層結(jié)構(gòu)位移功率譜

structure on the left

圖4 右側(cè)6層層間位移功率譜對(duì)比

6.2 本文方法精度及效率分析

為了驗(yàn)證本文所提組合體系響應(yīng)方差和譜矩封閉解的正確性，與虛擬激勵(lì)法進(jìn)行了對(duì)比分析。由式(37)可知，虛擬激勵(lì)法的結(jié)構(gòu)響應(yīng)功率譜無(wú)顯示表達(dá)式，而隨機(jī)振動(dòng)中的0階～2階譜矩是對(duì)功率譜的積分，故虛擬激勵(lì)法譜矩的計(jì)算只能采用梯形數(shù)值積分的方法進(jìn)行，譜矩的計(jì)算結(jié)果受積分步長(zhǎng)的影響較大。

從圖5～圖7可以看出，隨著虛擬激勵(lì)法積分步長(zhǎng)的逐漸變小，虛擬激勵(lì)法計(jì)算結(jié)果逐漸逼近本文方法，在理論上說(shuō)明本文方法的正確性。

圖5 積分步長(zhǎng)對(duì)左側(cè)結(jié)構(gòu)位移0階譜矩的影響分析

圖6 積分步長(zhǎng)對(duì)左側(cè)結(jié)構(gòu)位移1階譜矩的影響分析

圖7 積分步長(zhǎng)對(duì)左側(cè)結(jié)構(gòu)位移2階譜矩的影響分析

6.3 組合體結(jié)構(gòu)設(shè)置阻尼裝置前后減震效果對(duì)比

為驗(yàn)證相鄰結(jié)構(gòu)設(shè)置粘彈性阻尼器的減震效果，圖8和圖9給出了設(shè)置粘彈性阻尼器前后的位移及層間位移變化情況；圖10和圖11給出了設(shè)置粘彈性阻尼器前后的層間剪力及層間位移角對(duì)比分析。

從圖8可以看出，設(shè)置粘彈性阻尼器的組合體結(jié)構(gòu)能有效降低相鄰結(jié)構(gòu)在地震作用下的相對(duì)應(yīng)地面的絕對(duì)位移，左側(cè)結(jié)構(gòu)頂端位移減少12.7%，右側(cè)結(jié)構(gòu)頂端位移減少15.8%。從圖9可以看出，當(dāng)連接阻尼器設(shè)置在第6層時(shí)，左側(cè)結(jié)構(gòu)連接層以下的層間位移減小幅度較大，但連接層以上的層間位移出現(xiàn)少許增大；右側(cè)結(jié)構(gòu)1層～4層層間位移減小明顯，但6層層間位移出現(xiàn)少許增大。從圖9和圖10可以看出，層間剪力及層間彈性角的變化趨勢(shì)與層間位移一致。因此，相鄰結(jié)構(gòu)設(shè)置連接阻尼裝置時(shí)應(yīng)對(duì)層間位移進(jìn)行綜合評(píng)估來(lái)確定，以確保低層結(jié)構(gòu)的安全。從圖11可以看出，左右側(cè)建筑結(jié)構(gòu)的最大層間位移角發(fā)生在底層。設(shè)置阻尼器之前左側(cè)結(jié)構(gòu)1層的層間位移角為0.00737，未滿足現(xiàn)行高層建筑混凝土結(jié)構(gòu)技術(shù)規(guī)程的規(guī)范值；而設(shè)置阻尼器之后層間位移角最大為0.00626，小于規(guī)范值。右側(cè)結(jié)構(gòu)設(shè)置阻尼器前后層間位移角均滿足現(xiàn)行規(guī)范要求，但設(shè)置阻尼器之后，層間位移角明顯減小。

圖8 左右側(cè)結(jié)構(gòu)減震前后絕對(duì)位移對(duì)比分析

圖9 左右側(cè)結(jié)構(gòu)減震前后層間位移對(duì)比分析

圖10 左右側(cè)結(jié)構(gòu)減震前后層間剪力對(duì)比分析

圖11 左右側(cè)結(jié)構(gòu)減震前后層間位移角對(duì)比分析

7 結(jié) 論

本文針對(duì)相鄰建筑結(jié)構(gòu)設(shè)置Kelvin型粘彈性阻尼器組成的組合體結(jié)構(gòu)基于Kanai-Tajimi譜的相對(duì)于地面的絕對(duì)位移、層間位移等系列響應(yīng)的簡(jiǎn)明封閉解進(jìn)行了研究，獲得如下結(jié)論。

(1) 本文利用Kelvin型粘彈性阻尼器的微分型本構(gòu)關(guān)系及Kanai-Tajimi譜的濾波方程，與相鄰結(jié)構(gòu)的運(yùn)動(dòng)方程聯(lián)合建立全微分型且易于獲得簡(jiǎn)明解的白噪聲激勵(lì)的地震動(dòng)方程；利用復(fù)模態(tài)方法獲得組合體系的位移、層間位移的方差和0階～2階譜矩的簡(jiǎn)明統(tǒng)一解。

(2) 通過(guò)與虛擬激勵(lì)法對(duì)比分析，驗(yàn)證了本文所提的結(jié)構(gòu)系列響應(yīng)的方差和0階～2階譜矩的簡(jiǎn)明封閉解的正確性，避免了虛擬激勵(lì)法通過(guò)試算才能獲得結(jié)構(gòu)響應(yīng)方差和譜矩高精度解的缺點(diǎn)，因此，本文方法更為簡(jiǎn)便且具有較高的精度及效率。

(3) 通過(guò)對(duì)兩相鄰結(jié)構(gòu)設(shè)置連接阻尼裝置前后的結(jié)果對(duì)比，說(shuō)明設(shè)置連接阻尼裝置對(duì)相鄰結(jié)構(gòu)具有良好的減震效果，降低了地震作用下結(jié)構(gòu)碰撞發(fā)生的概率，但某些樓層的層間位移和層間剪力可能增大，故需對(duì)結(jié)構(gòu)層間位移進(jìn)行綜合評(píng)估。