付相球，潘旦光，2

（1.北京科技大學(xué)土木工程系，北京100083；2.北京科技大學(xué)城市地下空間工程北京市重點(diǎn)實(shí)驗(yàn)室，北京100083）

引言

非比例阻尼體系廣泛存在于現(xiàn)有的結(jié)構(gòu)體系中，如下部鋼筋混凝土上部鋼結(jié)構(gòu)的豎向混合結(jié)構(gòu)［1-2］、附加阻尼器的振動(dòng)控制結(jié)構(gòu)［3］，土-結(jié)構(gòu)相互作用體系［4-7］等。雖然直接積分法可用于非比例阻尼體系的動(dòng)力反應(yīng)分析，但計(jì)算工作量較大。模態(tài)疊加法計(jì)算效率高，且可識(shí)別對(duì)結(jié)構(gòu)動(dòng)力反應(yīng)具有顯著貢獻(xiàn)的模態(tài)，因而得到廣泛的應(yīng)用［8-9］。

對(duì)于非比例阻尼體系，F(xiàn)oss［10］提出以狀態(tài)空間法求解非比例阻尼系統(tǒng)的復(fù)模態(tài)，該方法可以得到精確的動(dòng)力反應(yīng)。但計(jì)算維數(shù)擴(kuò)大了一倍，且將實(shí)數(shù)域運(yùn)算轉(zhuǎn)換到復(fù)數(shù)域，計(jì)算量約為實(shí)模態(tài)分析的8倍［11］。為提高非比例阻尼體系模態(tài)疊加法的計(jì)算效率，Cronin［12］通過忽略阻尼矩陣中的非對(duì)角元素，強(qiáng)制解耦非比例阻尼體系以應(yīng)用實(shí)模態(tài)疊加法求解體系的近似反應(yīng)，簡(jiǎn)單易行但誤差難以控制［13-14］。Roesset等［15］提出等效阻尼比的概念用于近似解耦非比例阻尼，并常被用于混合結(jié)構(gòu)的地震反應(yīng)分析［16-17］，這種方法實(shí)際上也屬于強(qiáng)制解耦方法，各階模態(tài)的阻尼比在一定程度上考慮了耦合阻尼的影響，提高了計(jì)算精度，但不能解決計(jì)算誤差無法估計(jì)的問題［18］。Ibrahimbegovic等［11］提出了應(yīng)用實(shí)模態(tài)進(jìn)行非比例阻尼體系動(dòng)力反應(yīng)的迭代模態(tài)疊加法，這種方法將耦合阻尼的影響作為非線性荷載，得到高精度的計(jì)算結(jié)果。張靜等［19］針對(duì)大規(guī)模非對(duì)稱實(shí)矩陣的標(biāo)準(zhǔn)特征值問題，發(fā)展了適用于大型復(fù)特征值問題求解的Lanczos-QR方法，提高了計(jì)算效率，同時(shí)避免了范數(shù)誤差。事實(shí)上，非比例阻尼體系復(fù)模態(tài)疊加法的理論是完備的，問題在于狀態(tài)空間復(fù)模態(tài)求解及復(fù)數(shù)域的疊加計(jì)算時(shí)間較長。為提高復(fù)模態(tài)的計(jì)算效率，可將非比例阻尼系統(tǒng)看成由比例阻尼體系攝動(dòng)后得到的新系統(tǒng)，應(yīng)用攝動(dòng)法求解特征 方程［20-21］。Cha［22］提出一階攝動(dòng)方 法求解弱非比例阻尼系統(tǒng)；Tang和Wang［23］提出一種可以處理重根問題的攝動(dòng)方法。但一般的攝動(dòng)法只適用于弱非比例阻尼的體系，且需要用到原系統(tǒng)完整的模態(tài)空間。Lou等［24-25］提出的模態(tài)攝動(dòng)法利用原系統(tǒng)的非完整模態(tài)即可得到二階攝動(dòng)精度，并已應(yīng)用到無阻尼離散和連續(xù)系統(tǒng)特征值問題的分析。Pan等［26］應(yīng)用模態(tài)疊加法基本原理，構(gòu)建等效比例阻尼系統(tǒng)作為原系統(tǒng)，提出一種可求解強(qiáng)非比例體系特征值問題的計(jì)算方法，并以一個(gè)兩自由度體系和一個(gè)單層框架結(jié)構(gòu)驗(yàn)證計(jì)算方法的精度。

模態(tài)攝動(dòng)法本質(zhì)上屬于Ritz法，計(jì)算精度與附加模態(tài)的數(shù)目有關(guān)，文獻(xiàn)［26］中兩個(gè)算例的自由度都比較小，采用完整模態(tài)或接近完整模態(tài)進(jìn)行計(jì)算，結(jié)果的精度都很高。本文在文獻(xiàn)［26］的基礎(chǔ)上，進(jìn)一步研究自由度較多體系非完整模態(tài)下模態(tài)攝動(dòng)法的計(jì)算精度，在此基礎(chǔ)上，建立一種基于無阻尼體系實(shí)模態(tài)的復(fù)模態(tài)疊加法。同時(shí)，以非比例阻尼體系的地震反應(yīng)為例，討論復(fù)模態(tài)疊加法中的模態(tài)截?cái)鄦栴}，分析了基于累積振型參與質(zhì)量、累積振型貢獻(xiàn)系數(shù)及累積振型加速度貢獻(xiàn)系數(shù)所得模態(tài)數(shù)對(duì)計(jì)算精度的影響，并對(duì)比研究本文方法和狀態(tài)空間復(fù)模態(tài)疊加法計(jì)算效率的差異。最后，將本文方法應(yīng)用到一個(gè)帶附加阻尼鋼結(jié)構(gòu)框架地震反應(yīng)分析，驗(yàn)證了本文方法的適用性。

1 基于實(shí)模態(tài)的復(fù)模態(tài)疊加法

一個(gè)N自由度線性黏滯阻尼體系的強(qiáng)迫振動(dòng)方程為

式中u，u?和u?分別為N×1階的位移、速度和加速度向量；f(t)為荷載向量；m，c和k分別為N×N階的質(zhì)量、阻尼和剛度矩陣。對(duì)于非比例阻尼體系，狀態(tài)空間法把N維2階微分方程轉(zhuǎn)變?yōu)?N維1階微分方程［10］

為應(yīng)用模態(tài)疊加法求解式（2）的運(yùn)動(dòng)方程，須求解特征方程

式中γ和ψ分別為復(fù)特征值和相應(yīng)的特征向量。若已知前r（r≤N）對(duì)特征值γj和γj+r=γˉj，以及相應(yīng)的特征向量ψj和ψj+r=ψˉj（j=1，2，…，r），其中“”表示復(fù)數(shù)共軛。則式（2）強(qiáng)迫振動(dòng)的反應(yīng)y可表示為

式中zj為廣義坐標(biāo)，zˉj表示zj的共軛復(fù)數(shù)，r為截?cái)嗟哪B(tài)數(shù)。由模態(tài)疊加法的基本原理，廣義坐標(biāo)zj（j=1，2，…，r）的運(yùn)動(dòng)方程為

從計(jì)算過程看，基于狀態(tài)空間的復(fù)模態(tài)疊加法與常規(guī)模態(tài)疊加法相同，但是，式（3）特征值求解的維數(shù)增加了1倍且為復(fù)數(shù)運(yùn)算，使計(jì)算時(shí)間顯著增加。為提高計(jì)算效率，關(guān)鍵是減少式（3）和式（4）的計(jì)算時(shí)間。為此，本文基于模態(tài)攝動(dòng)法原理，以無阻尼體系的實(shí)模態(tài)為基礎(chǔ)，計(jì)算非比例阻尼的復(fù)模態(tài)，以減少計(jì)算時(shí)間。

非比例阻尼矩陣可分解為

式中cp為相應(yīng)的比例阻尼矩陣，Δc為比例阻尼矩陣與原阻尼矩陣的偏差。若已知無阻尼體系特征方程前n階自振頻率ωj和相應(yīng)的模態(tài)φj（j=1，2，…，n），則比例阻尼矩陣cp滿足正交特性

將m，cp和k組成的系統(tǒng)稱之為等效比例阻尼系統(tǒng)，作為非比例阻尼體系的原系統(tǒng)。則采用狀態(tài)空間法表示的等效比例阻尼體系中的特征值方程為

非比例阻尼體系新系統(tǒng)的特征值γj和模態(tài)ψj可以利用原系統(tǒng)的特征值sj和模態(tài)ηj進(jìn)行簡(jiǎn)單的攝動(dòng)分析而近似地求得［26］，即設(shè)：

式 中D11，D12，D21，D22，E11，E12，E21和E22都 為n×n階 的 方 陣；R1和R2為n×1階向 量；x為2n×1階向量。各矩陣和向量的元素具體表達(dá)見附錄。

式（12）將復(fù)特征方程（3）轉(zhuǎn)換成了2n維的非線性代數(shù)方程組，求解代數(shù)方程通常比求解特征方程要簡(jiǎn)單得多，本文采用Newton-Raphson迭代方法求解。在得到未知向量x后，非比例阻尼體系的第j個(gè)特征值為

采用有阻尼體系的頻率、阻尼比和特征值的表示方法

則：

式中Re代表取元素的實(shí)部。為了避免與相應(yīng)的等效比例阻尼的自振頻率ωj和阻尼比ζj混淆，和表示第j階的擬無阻尼自振頻率和阻尼比。當(dāng)式（10）中的n=N時(shí)，2N個(gè)線性無關(guān)的特征向量ηj組成2N維的完備空間，因此，非比例阻尼體系的特征向量Ψj可以在這個(gè)2N維的空間中精確地展開，此時(shí)攝動(dòng)解也將收斂到狀態(tài)空間下的精確解。而通常高階模態(tài)相比于低階模態(tài)影響很小，因此方程（12）的維數(shù)n通常比N小得多，這表明本文所采用的攝動(dòng)法無需求解原系統(tǒng)所有階模態(tài)即可得到高精度的攝動(dòng)解，從而減少計(jì)算時(shí)間。

式（12）中的j從1到r依次取值可得非比例阻尼體系前r階復(fù)特征值和復(fù)模態(tài)。由此，式（5）中pj(t)和aj用實(shí)模態(tài)表示為：

方程（5）得到廣義坐標(biāo)zj的解。令其中quj與qlj分別為qj的前n和后n個(gè)元素組成的向量。由式（4）和（10）可得體系的位移為

式中qu和ql分別為r個(gè)quj與qlj向量組成的n×r維的矩陣，z={z1z2...zr}T為r×1的廣義坐標(biāo)。由式（16），式（17）和式（18）可知，基于實(shí)模態(tài)的復(fù)模態(tài)疊加法將復(fù)模態(tài)、廣義單自由度體系的部分系數(shù)和位移表示為實(shí)模態(tài)的線性組合，簡(jiǎn)化了計(jì)算。

注重農(nóng)業(yè)技術(shù)創(chuàng)新，加快完善和填補(bǔ)優(yōu)勢(shì)農(nóng)畜產(chǎn)品和特色農(nóng)畜產(chǎn)品標(biāo)準(zhǔn)化生產(chǎn)操作規(guī)程，加強(qiáng)標(biāo)準(zhǔn)推廣和使用指導(dǎo)，大力宣傳培訓(xùn)農(nóng)牧業(yè)投入品使用規(guī)范，督促生產(chǎn)經(jīng)營主體嚴(yán)格落實(shí)間隔期休藥期規(guī)定。發(fā)展農(nóng)牧業(yè)適度規(guī)模經(jīng)營，推進(jìn)農(nóng)畜產(chǎn)品標(biāo)準(zhǔn)化生產(chǎn)基地建設(shè)，提高農(nóng)畜產(chǎn)品生產(chǎn)經(jīng)營主體的專業(yè)化、組織化程度，先行對(duì)農(nóng)畜產(chǎn)品新型經(jīng)營主體開展農(nóng)畜產(chǎn)品質(zhì)量安全綜合指數(shù)評(píng)價(jià)，將是否按標(biāo)生產(chǎn)作為政策支持的重要條件，從源頭上增強(qiáng)農(nóng)畜產(chǎn)品經(jīng)營主體的安全意識(shí)，推進(jìn)農(nóng)畜產(chǎn)品的安全生產(chǎn)。

2 算例分析

為研究地震作用下基于實(shí)模態(tài)的復(fù)模態(tài)疊加法適用性，借鑒文獻(xiàn)［11］的算例構(gòu)建圖1所示的框架結(jié)構(gòu)，討論模態(tài)攝動(dòng)法的精度、模態(tài)截?cái)嘀笜?biāo)以及復(fù)模態(tài)疊加法的計(jì)算精度和效率。該框架有限元模型包括27個(gè)梁?jiǎn)卧?2個(gè)自由度，材料的彈性模量和密度分別為40 N/m2，1 kg/m3，梁截面慣性矩和面積分別為1 m4，1 m2，單元長度為1 m。采用集中質(zhì)量模型，由靜力凝聚消除轉(zhuǎn)動(dòng)自由度后，體系包含48個(gè)有效無阻尼實(shí)模態(tài)?？蚣芙Y(jié)構(gòu)在節(jié)點(diǎn)3，8，12，17，20與25的水平方向有附加阻尼器，阻尼器系數(shù)分別為c1=3 N·s/m，c2=5 N·s/m。以耦合系數(shù)α表示非比例阻尼矩陣耦合的程度［27］

圖1 附加阻尼器的三層框架結(jié)構(gòu)Fig.1 Three-layer frame with concentrated dampers

式中Clk表示阻尼矩陣C中的第l行第k列元素。體系的耦合系數(shù)為1，為強(qiáng)非比例阻尼體系。輸入地震時(shí)程為1940年5月在加州地震中記錄到的El Centro波的N-S方向的加速度分量，如圖2所示。

圖2 El Centro波加速度時(shí)程Fig.2 El Centro ground acceleration history

2.1 模態(tài)攝動(dòng)法的精度

當(dāng)模態(tài)攝動(dòng)法采用完整的無阻尼體系特征向量空間計(jì)算時(shí)，其結(jié)果收斂到精確解，但大型復(fù)雜結(jié)構(gòu)通常僅計(jì)算有限階模態(tài)，為保證計(jì)算精度，計(jì)算非比例阻尼體系第j階的模態(tài)時(shí)，式（10）中所需的無阻尼體系的模態(tài)數(shù)n為式中Δn為計(jì)算所需的附加模態(tài)數(shù)。以狀態(tài)空間法所得的擬無阻尼頻率與阻尼比為精確解，n階模態(tài)下由模態(tài)攝動(dòng)法得到的和為近似解，則擬無阻尼頻率與阻尼比的相對(duì)誤差為：

圖3所示為模態(tài)攝動(dòng)法所得前3階模態(tài)相對(duì)誤差隨附加模態(tài)數(shù)的變化情況。表1列出前3階采用不同方法所得擬無阻尼頻率與阻尼比的計(jì)算結(jié)果，表2為不同方法的計(jì)算誤差。其中n=48表示完整模態(tài)空間，強(qiáng)制解耦法為直接去除模態(tài)阻尼矩陣的非對(duì)角元素的解耦方法。

表1 框架結(jié)構(gòu)的擬無阻尼自振頻率與阻尼比Tab.1 The quasi-undamped natural frequencies and damping ratios of the frame system

表2 框架結(jié)構(gòu)的擬無阻尼自振頻率與阻尼比的相對(duì)誤差/%Tab.2 The errors of quasi-undamped natural frequencies and damping ratios of the frame system/%

圖3 前3階頻率與阻尼比的相對(duì)誤差Fig.3 Relative errors of the first three natural frequencies and damping ratios

計(jì)算結(jié)果表明：（1）對(duì)于模態(tài)攝動(dòng)法，隨著附加模態(tài)數(shù)的增加，前3階頻率和阻尼比趨近于精確解；采用完整實(shí)模態(tài)的模態(tài)攝動(dòng)法計(jì)算結(jié)果與精確解一致。附加模態(tài)數(shù)為8時(shí)，前3階頻率的相對(duì)誤差小于1%，前3階模態(tài)阻尼比的相對(duì)誤差小于2%。這表明對(duì)于大型復(fù)雜結(jié)構(gòu)，當(dāng)采用附加模態(tài)數(shù)大于8的非完整模態(tài)進(jìn)行模態(tài)攝動(dòng)法計(jì)算時(shí)，即可得到很高的計(jì)算精度。（2）強(qiáng)制解耦方法計(jì)算所得的頻率誤差小于5%，但阻尼比誤差可達(dá)9.596%。這是因?yàn)閺?qiáng)制解耦過程改變了阻尼矩陣，直接對(duì)模態(tài)阻尼比產(chǎn)生影響，因此阻尼比與精確解存在明顯的差別；而由單自由度體系有阻尼頻率可知，當(dāng)阻尼比小于0.2時(shí)，阻尼引起體系有阻尼頻率的變化很小，因此強(qiáng)制解耦方法計(jì)算時(shí)頻率與精確解相差不大。

2.2 復(fù)模態(tài)疊加法地震反應(yīng)分析

對(duì)于一般結(jié)構(gòu)而言，少數(shù)低階模態(tài)的疊加即可得到滿足工程要求的解，因此存在模態(tài)截?cái)嗟膯栴}［28］。然而對(duì)于非比例阻尼體系的復(fù)模態(tài)疊加法，能否采用實(shí)模態(tài)疊加法中的模態(tài)截?cái)喾椒ㄐ枰M(jìn)一步的研究。針對(duì)本文的三層框架結(jié)構(gòu)，三類實(shí)模態(tài)的模態(tài)截?cái)嘀笜?biāo)：累積振型參與質(zhì)量、首層平動(dòng)位移的累積振型貢獻(xiàn)系數(shù)以及首層平動(dòng)的累積振型加速度貢獻(xiàn)系數(shù)如表3所示。以模態(tài)截?cái)嘀笜?biāo)超過90%作為模態(tài)選取的依據(jù)，則由累積振型參與質(zhì)量、首層平動(dòng)位移的累積振型貢獻(xiàn)系數(shù)和首層平動(dòng)的累積振型加速度貢獻(xiàn)系數(shù)選取的模態(tài)數(shù)分別為5階，3階和20階。

表3 不同模態(tài)截?cái)嘀笜?biāo)計(jì)算結(jié)果/%Tab.3 Results of different modal truncation indexes/%

式中ei-max表示第i個(gè)自由度位移與精確解的峰值誤差，ei-sum表示第i個(gè)自由度位移與精確解的累積誤差表示第i個(gè)自由度位移的精確解，ui(t)表示第i個(gè)自由度位移的近似解，T為積分時(shí)長。在El Centro波作用下，前3階、前5階、前20階和全部48階復(fù)模態(tài)采用式（18）計(jì)算所得頂層、中間層、底層水平位移時(shí)程如圖4所示，在采用模態(tài)攝動(dòng)法計(jì)算特征值時(shí)，附加實(shí)模態(tài)的個(gè)數(shù)均取8。同時(shí)為了進(jìn)行比較，將狀態(tài)空間法精確解也一并繪出。各自由度的峰值誤差、累積誤差如表4所示。

由圖4與表4可以看出：1）當(dāng)采用48階復(fù)模態(tài)疊加計(jì)算時(shí)，計(jì)算結(jié)果與狀態(tài)空間法完全重合，說明完整模態(tài)的復(fù)模態(tài)疊加法可以得到精確解；2）底層水平位移的模態(tài)截?cái)嗾`差最大，頂層水平位移的模態(tài)截?cái)嗾`差最小，說明底層水平位移反應(yīng)對(duì)高階模態(tài)更為敏感；3）當(dāng)采用前3階、前5階復(fù)模態(tài)疊加計(jì)算時(shí)，各層水平位移的誤差較大，其中底層水平位移峰值誤差達(dá)到30%，累積誤差超過了50%，說明對(duì)于該框架體系，采用累積振型參與質(zhì)量和首層平動(dòng)位移的累積振型貢獻(xiàn)系數(shù)作為模態(tài)截?cái)嘀笜?biāo)難以滿足計(jì)算精度要求；4）當(dāng)采用前20階復(fù)模態(tài)疊加計(jì)算時(shí)，頂層和中間層的水平位移峰值誤差及累積誤差均在10%以內(nèi)，底層水平位移累積誤差偏大，略大于10%，但峰值誤差只有5.8%，結(jié)合圖4（c）的位移時(shí)程結(jié)果可以看出前20階復(fù)模態(tài)疊加時(shí)底層水平位移與精確解結(jié)果較吻合。總體而言，累積振型參與質(zhì)量和首層平動(dòng)位移的累積振型貢獻(xiàn)系數(shù)所取的復(fù)模態(tài)數(shù)偏少，基于首層平動(dòng)的累積振型加速度貢獻(xiàn)系數(shù)的模態(tài)截?cái)喾椒梢詽M足計(jì)算精度要求。

表4 三層框架結(jié)構(gòu)水平位移誤差/%Tab.4 Errors of horizontal displacements of the threelayer frame/%

圖4 地震作用下水平位移計(jì)算結(jié)果Fig.4 Horizontal displacement under earthquake excitation

為對(duì)比不同方法的計(jì)算精度，圖5給出了前20階模態(tài)下本文方法與狀態(tài)空間法和強(qiáng)制解耦非比例阻尼矩陣后的實(shí)模態(tài)疊加法（簡(jiǎn)稱實(shí)模態(tài)疊加法）的水平位移計(jì)算結(jié)果，并與精確解相比較。不同方法相應(yīng)的峰值誤差及累積誤差如表5所示。

由圖5、表5結(jié)果可以看出：1）本文方法與狀態(tài)空間法的前20階計(jì)算結(jié)果接近，兩種方法的時(shí)程曲線與精確解吻合很好；三層水平位移的峰值誤差均小于10%，除底層位移累積誤差偏大高于10%以外，頂層及中間層位移累積誤差也均小于10%，顯示了良好的精度；2）實(shí)模態(tài)疊加法計(jì)算結(jié)果與精確解存在明顯誤差，三層水平位移的累積誤差均大于20%，其中底層位移的峰值誤差也高于20%，這說明對(duì)于強(qiáng)非比例阻尼的體系，實(shí)模態(tài)疊加方法忽略了非對(duì)角線阻尼的影響，導(dǎo)致較大的計(jì)算誤差。

表5 三層框架結(jié)構(gòu)水平位移誤差/%Tab.5 Errors of horizontal displacements of the threelayer frame/%

圖5 地震作用下三種方法的水平位移計(jì)算結(jié)果Fig.5 Horizontal displacement of three methods under earthquake excitation

為對(duì)比本文方法與狀態(tài)空間法的計(jì)算效率，前20階模態(tài)疊加下計(jì)算時(shí)間如表6所示，計(jì)算過程統(tǒng)一采用2.19 GHz處理器的筆記本，Matlab計(jì)算軟件?？梢钥闯觯B(tài)攝動(dòng)法求解特征方程所耗費(fèi)的時(shí)間小于狀態(tài)空間法直接求解復(fù)特征方程的時(shí)間，說明模態(tài)攝動(dòng)法提高了求解特征方程的效率；本文方法由于計(jì)算過程中部分采用實(shí)模態(tài)線性組合，減少了復(fù)數(shù)運(yùn)算，計(jì)算效率也明顯高于狀態(tài)空間方法。綜合計(jì)算精度結(jié)果，本文方法兼顧了計(jì)算效率與精度，適用于強(qiáng)非比例阻尼體系的地震反應(yīng)分析。

表6 計(jì)算時(shí)間/msTab.6 Calculation time/ms

3 五層鋼結(jié)構(gòu)地震反應(yīng)分析

為驗(yàn)證本文方法對(duì)真實(shí)結(jié)構(gòu)的適用性，下面采用不同方法分析北京市通州區(qū)某公司總部辦公樓的地震反應(yīng)。該辦公樓主體為鋼框架結(jié)構(gòu)，層高17.7 m，建筑長50.4 m，分為7跨，其中一榀橫向框架鋼結(jié)構(gòu)計(jì)算簡(jiǎn)圖如圖6所示。鋼結(jié)構(gòu)的梁和柱采用Q345鋼，樓板選用C25混凝土。橫向7.2 m，6 m跨度框架梁截面為HN450 mm×200 mm×9 mm×14 mm，橫向3 m跨度框架梁截面為HN300 mm×150 mm×6.5 mm×9 mm（數(shù)字代表工字鋼的截面高度，截面寬度，腹板厚度，翼緣厚度）；底層柱為400 mm×400 mm×16 mm×16 mm箱型截面，其他層柱為400 mm×400 mm×14 mm×14 mm（數(shù)字代表箱型柱（方管）截面高度，截面寬度，上下厚度，左右厚度）箱型截面。辦公樓重力荷載如樓板、墻體等以等效質(zhì)量的方式施加在框架結(jié)構(gòu)上，其中屋面重力荷載取6.38 kN/m2，樓面重力荷載取7 kN/m2。鋼結(jié)構(gòu)阻尼比取0.02，采用瑞利阻尼分析，瑞利阻尼系數(shù)的參考頻率選用基頻與荷載卓越頻率。同時(shí)為提高該框架結(jié)構(gòu)的抗震性能，各樓層橫向A-B跨和CD跨之間各設(shè)置一個(gè)黏滯阻尼器，采用斜對(duì)角布置如圖6所示，附加阻尼器的阻尼系數(shù)為c=1500 kN·s/m。

圖6 辦公樓計(jì)算簡(jiǎn)圖Fig.6 Calculation diagram of an office building

該框架結(jié)構(gòu)的基本自振周期為1 s，耦合系數(shù)α=0.57。累積振型加速度貢獻(xiàn)系數(shù)超過90%時(shí)，所需的最少模態(tài)數(shù)為7階。在El Centro波作用下，該框架結(jié)構(gòu)的底層和頂層位移反應(yīng)時(shí)程如圖7所示，峰值與累積誤差如表7所示。其中模態(tài)攝動(dòng)法計(jì)算特征值時(shí)，附加實(shí)模態(tài)的個(gè)數(shù)仍取8。

圖7 El Centro波作用下三種方法的水平位移計(jì)算結(jié)果Fig.7 Horizontal displacement of three methods under El Centro wave

由圖7及表7可以看出，本文方法所得底層與頂層位移的峰值誤差均在1%以內(nèi)，累積誤差不超過5%，具有很高的計(jì)算精度。而實(shí)模態(tài)疊加法底層位移峰值誤差大于10%，累積誤差達(dá)29.40%?？傮w而言，本文方法計(jì)算精度高，可用于強(qiáng)非比例阻尼體系的地震反應(yīng)分析，而強(qiáng)制解耦的實(shí)模態(tài)疊加法易造成無法預(yù)計(jì)的誤差。

表7 結(jié)構(gòu)水平位移誤差/%Tab.7 Errors of horizontal displacements of the frame/%

4 結(jié)論

工程實(shí)際結(jié)構(gòu)中的阻尼常具有非比例阻尼的特征，本文基于模態(tài)攝動(dòng)法提出基于實(shí)模態(tài)的復(fù)模態(tài)疊加法，用于地震荷載下的強(qiáng)迫振動(dòng)分析。通過理論分析和數(shù)值計(jì)算可以得出以下結(jié)論：

1）當(dāng)采用無阻尼體系所有模態(tài)構(gòu)成的完備空間進(jìn)行模態(tài)攝動(dòng)法計(jì)算時(shí)，可以得到非比例阻尼體系精確的復(fù)特征值和特征向量；對(duì)于非完備空間，當(dāng)附加模態(tài)數(shù)大于8時(shí)，模態(tài)攝動(dòng)法所得的復(fù)特征值誤差小于2%；

2）非比例阻尼體系的復(fù)模態(tài)疊加法計(jì)算時(shí)，基于累積振型參與質(zhì)量及累積振型貢獻(xiàn)系數(shù)進(jìn)行模態(tài)截?cái)嗨玫哪B(tài)數(shù)偏少，由此導(dǎo)致結(jié)構(gòu)地震反應(yīng)的誤差較大，建議采用累積振型加速度貢獻(xiàn)系數(shù)作為復(fù)模態(tài)疊加法的模態(tài)截?cái)嘀笜?biāo)；

3）本文方法將非比例阻尼體系的特征向量和強(qiáng)迫振動(dòng)的解采用實(shí)模態(tài)線性組合進(jìn)行計(jì)算，從而減少了復(fù)數(shù)運(yùn)算，提高了計(jì)算效率；

4）對(duì)于強(qiáng)非比例阻尼體系的地震反應(yīng)分析，本文方法兼顧了計(jì)算效率與精度，而強(qiáng)制解耦非比例阻尼矩陣的實(shí)模態(tài)疊加法忽略了耦合阻尼的影響，計(jì)算誤差難以控制。

附錄：