周清橋

0 引言

結(jié)構(gòu)的抗震計(jì)算分析方法有很多，本文根據(jù)高層建筑結(jié)構(gòu)的特點(diǎn)，對(duì)其進(jìn)行抗震計(jì)算分析。

1 計(jì)算模型

多層結(jié)構(gòu)往往采用層剪切模型來(lái)計(jì)算。而像更高層數(shù)的結(jié)構(gòu)，其結(jié)構(gòu)往往是筒體或框架剪力墻結(jié)構(gòu)，則適合采用彎剪層模型。而桿模型以梁、柱等單根構(gòu)件為基礎(chǔ)單元，將樓層質(zhì)量分別集中于結(jié)構(gòu)各節(jié)點(diǎn)，形成質(zhì)點(diǎn)。桿模型是較為精確的計(jì)算模型，我們采用的計(jì)算模型是層剪切模型，其方程如下［1］:

2 數(shù)值計(jì)算方法

2． 1 結(jié)構(gòu)非線性計(jì)算方法

式(1)是非線性方程，需要用數(shù)值方法求解。但這里我們采用比較有代表性的Newton-Raphson法，在進(jìn)行非線性求解時(shí)具體采用增量Newton-Raphson方法，此方法是先把荷載分成若干增量(或荷載步)，然后在每個(gè)增量步內(nèi)用Newton-Raphson方法進(jìn)行迭代求解。

圖1是增量Newton-Raphson方法求解非線性方程F=Kδ的圖解。圖1中荷載 F 分成了 i個(gè)增量步(F1，F(xiàn)2，…，F(xiàn)i-1，F(xiàn)i)，現(xiàn)以第一增量步為例說(shuō)明Newton-Raphson方法的迭代過(guò)程。在直角坐標(biāo)系F—δ中，第一增量步內(nèi)直線F1=R和曲線F=Kδ的交點(diǎn)A的橫坐標(biāo)δA就是方程F=Kδ在第一增量步內(nèi)的解。

圖1 增量Newton-Raphson方法

在實(shí)際使用中，Newton-Raphson法開(kāi)始需先按線性理論求出位移δ1作為第一次迭代近似值，即圖1中A1點(diǎn)的橫坐標(biāo)，如果荷載R并不因變形而改變它的數(shù)值和方向(假設(shè)為小變形)，則得［2］:

其中，KT為曲線F=Kδ的斜率(或切線剛度)。

第二次迭代從B1點(diǎn)做曲線F=Kδ的切線交直線F1=R于A2點(diǎn)，取A2點(diǎn)的橫坐標(biāo)δ2，從圖1中可見(jiàn):

由圖1看出δ2就是位移的第二次迭代近似值。如此重復(fù)可得以下迭代公式:

不斷重復(fù)迭代直至滿足收斂條件為止，便可得到最終結(jié)果。其他增量步內(nèi)的求解過(guò)程與第一增量步內(nèi)相同，實(shí)踐證明該法收斂性較好。

2． 2 矩陣特征值與特征向量的求解

對(duì)結(jié)構(gòu)進(jìn)行振型分析，求特征值:

注意到式(1)并不是標(biāo)準(zhǔn)的特征值表達(dá)式，可采用矩陣運(yùn)算的方法把式(1)化為標(biāo)準(zhǔn)的特征值表達(dá)式。

在式(1)左右同時(shí)乘以［M-1］，即可化為:

令 A=［K］［M-1］，λ =ω2，式(1)即可變成標(biāo)準(zhǔn)的特征值表達(dá)式。這里［E］為單位剛度矩陣。

當(dāng)然也可以左右同時(shí)乘以柔度矩陣［K-1］，兩者是等效的。

上式問(wèn)題常用的解法有雅可比法、迭代法以及近年來(lái)發(fā)展的子空間迭代法，此處采用的方法屬于迭代法求解。由于總剛矩陣［A］并不是實(shí)對(duì)稱矩陣，不適用雅可比法，而須先把其化解為赫申伯格矩陣，然后用QR方法來(lái)求解其特征值。

1)約化赫申伯格步驟，對(duì)k=2，3…作如下變換:

a．從第k-1列的第一個(gè)以下的元素中選出絕對(duì)值最大的元素 al，k-1。

b．交換第l1行與第k行，交換第l1列與第k列。

c．對(duì)于 i=k+1，…，n 作變換:

2)求赫申伯格矩陣全部特征值的QR方法步驟:

設(shè)不可約的上H矩陣為A。

a．確定一個(gè)初等正交矩陣Q0:

對(duì)A作相似變換A1=Q0AQ0。

b．利用同樣的方法確定正交對(duì)稱矩陣 Q1，Q2，…，Qn-2，對(duì)A1，A2，…，An-2作相似變換:

反復(fù)進(jìn)行以上迭代，直到將上H矩陣A變?yōu)閷?duì)角塊全部是一階塊或二階塊為止，此時(shí)就可以直接從各一階塊或二階塊解出全部特征值。

3)求特征矩陣的特征向量:

利用求得的特征值代入特征方程，采用全選主元消去法的步驟:

設(shè)方程為AX=B。

消去過(guò)程:

a．全選主元。即從系數(shù)矩陣的第k行、第k列開(kāi)始的右下子陣中選取絕對(duì)值最大的元素，并將它交換到主元素的位置上。

b．歸一化。即:

c．消去。即:

回代過(guò)程:

2． 3 結(jié)構(gòu)動(dòng)力非線性分析的Newmark-β方法

在結(jié)構(gòu)動(dòng)力非線性分析中，本文使用Newmark-β時(shí)間積分方法，在離散的時(shí)間點(diǎn)上求解這些方程，在tn+Δt時(shí)刻，動(dòng)力方程式為:

其中，［M］，［C］，［K］分別為質(zhì)量矩陣、阻尼矩陣和剛度矩陣;{}，}，{y}分別為加速度矢量、速度矢量和位移n+1矢量;{fn+1}為作用力矢量。

主要計(jì)算步驟:

1)初始計(jì)算。

b．選擇時(shí)間步長(zhǎng) Δt，參數(shù) β，γ，計(jì)算積分常數(shù)。

c．形成等效剛度矩陣［⌒K］=［K］+a0［M］+a1［C］。

2)對(duì)每一時(shí)刻計(jì)算。

a．計(jì)算t+Δt時(shí)刻等效荷載:

c．計(jì)算t+Δt時(shí)刻的加速度和速度:

在Newmark-β法中，參數(shù)β，γ的取值影響著算法的精度和穩(wěn)定性，且迭代的步長(zhǎng)決定迭代的收性。為保證算法具有不低于二階精度，要求參數(shù)β取值為1/2;一般情況下，參數(shù)γ取1/8～1/4即可獲得穩(wěn)定的解。本文中 β取1/2，γ取1/4，為有條件穩(wěn)定［3，4］。

3 計(jì)算結(jié)果

以某12層框架結(jié)構(gòu)為例，因?yàn)槠錁菍虞^低，其抗側(cè)力構(gòu)件薄弱，且層數(shù)較低，因此結(jié)構(gòu)的計(jì)算模型采用層間剪切模型。鑒于計(jì)算結(jié)果的工程判斷以模型的層間剪力和層間位移為主，通常以等效層間模型為主要的分析模型，其結(jié)果與用桿單元的計(jì)算模型精度相差不超過(guò)5%。在選用場(chǎng)地波形方面，我們選用最不利波El-Centro波。而El-Centro波計(jì)算時(shí)原結(jié)構(gòu)的位移反應(yīng)最大。El-Centro波為振動(dòng)型破壞地震波。Taft波為沖擊型破壞地震波。人工隨機(jī)波為能量型地震波。地震波對(duì)結(jié)構(gòu)的作用自然應(yīng)該首先體現(xiàn)在振動(dòng)型破壞作用上。為了分析共振效應(yīng)，我們進(jìn)一步分析了El-Centro波，下面我們列出這El-Centro波的傅立葉譜，見(jiàn)圖2。

圖2 El-Centro E-W傅立葉譜

經(jīng)過(guò)分析各條波的傅立葉譜，我們發(fā)現(xiàn)各條波的卓越圓頻率段為:El-Centro E-W波—0 s-1～3．75 s-1。以上四條波都將波峰集中在結(jié)構(gòu)的圓頻率附近，所以易發(fā)生共振［5］。

［1］ 張 敏．高層框架動(dòng)力分析剛度矩陣的彎剪層模型［J］．華東交通大學(xué)學(xué)報(bào)，2005，22(5):34-36．

［2］ 鐘萬(wàn)勰，何 窮，劉正興．?dāng)?shù)值計(jì)算方法［M］．北京:中國(guó)建筑工業(yè)出版社，1991:130-138．

［3］ 沈聚敏，周錫元．抗震工程學(xué)［M］．北京:中國(guó)建筑工業(yè)出版社，2000:1-299．

［4］ 呂西林．高層建筑結(jié)構(gòu)［M］．武漢:武漢工業(yè)大學(xué)出版社，2001:44-50．

［5］ 李 輝，曹 亮．在役結(jié)構(gòu)的抗震分析［J］．浙江工業(yè)大學(xué)學(xué)報(bào)，2005，33(2):223-226．