【摘要】本文介紹了非線性振型疊加法的研究成果，并建議采用指定阻尼比和振型疊加法相結(jié)合的方法處理橋梁結(jié)構(gòu)的阻尼以及橋梁局部非線性地震反應(yīng)分析的問(wèn)題。
關(guān)鍵詞 局部非線性 阻尼 振型疊加法
一、引言
近幾十年來(lái)，為促進(jìn)經(jīng)濟(jì)的持續(xù)發(fā)展，我國(guó)政府投入巨資修建交通基礎(chǔ)設(shè)施，橋梁建設(shè)突飛猛進(jìn)，相繼建成了一批大跨斜拉橋、懸索橋等重要工程。同時(shí)，據(jù)專家預(yù)測(cè)，我國(guó)已進(jìn)入一個(gè)新的地震活躍期。本世紀(jì)幾次慘重的地震災(zāi)害給了我們深刻的教訓(xùn)：對(duì)這些重大建設(shè)工程，必須進(jìn)行結(jié)構(gòu)的地震反應(yīng)分析，以評(píng)估結(jié)構(gòu)的地震作用下的安全程度。
目前，地震反應(yīng)分析主要有反應(yīng)譜分析方法和動(dòng)態(tài)時(shí)程分析方法兩種方法，對(duì)于大跨、復(fù)雜的橋梁一般以動(dòng)態(tài)時(shí)程分析方法為主。動(dòng)態(tài)時(shí)程分析方法按照不同的解方程思路分為振型疊加法和直接積分法兩種。無(wú)論采用哪種計(jì)算方法，都涉及到三個(gè)物理參數(shù)的確定，即結(jié)構(gòu)的質(zhì)量測(cè)度和阻尼。質(zhì)量和剛度已經(jīng)有了很成熟的計(jì)算方法，結(jié)果也較穩(wěn)定、可靠，但阻尼的計(jì)算則是一個(gè)比較復(fù)雜的迄今為止未得到很好解決的問(wèn)題。橋梁結(jié)構(gòu)中不同材料的構(gòu)件具有不同的阻尼，大量減隔震裝置的使用再增添了阻尼的多樣化和復(fù)雜性。由于阻尼參數(shù)確定的不同，常導(dǎo)致結(jié)構(gòu)反應(yīng)的計(jì)算結(jié)果有顯著的差異，有時(shí)甚至是幾倍甚至是數(shù)量級(jí)上的差別。因此如何處理阻尼問(wèn)題是正確求解大跨、復(fù)雜橋梁結(jié)構(gòu)地震反應(yīng)的關(guān)鍵問(wèn)題之一。
根據(jù)已發(fā)生歷次地震的橋梁震害記錄，我們發(fā)現(xiàn)震害區(qū)斜拉橋、懸索橋等大跨橋梁基本上未出現(xiàn)破壞，筆者所知的只有臺(tái)灣"9.21大地震"震中區(qū)域一未完工的獨(dú)塔斜拉橋主塔塔柱與橫梁交界處出現(xiàn)破裂，但這次地震很特殊，震源非常淺，加速度達(dá)到蛇，像這樣大的地震如何設(shè)防目前還處在研究階段。此外，根據(jù)同濟(jì)大學(xué)土木防災(zāi)重點(diǎn)實(shí)驗(yàn)室所做的幾十座國(guó)內(nèi)大跨橋的抗震分析，發(fā)現(xiàn)大跨斜拉橋、懸索橋的抗震薄弱環(huán)節(jié)是邊墩、輔助墩等部位，這些部位在地震作用下有可能進(jìn)入彈塑性階段，而主塔、主梁等部位基本上處于彈性階段。因此，我們把這種在地震荷載作用下，主體部位基本處在彈性階段，局部進(jìn)入彈塑性階段的情況稱之為局部非線性的地震反應(yīng)。
對(duì)于非線性問(wèn)題的時(shí)程計(jì)算，我們一般不采用振型疊加法。這是因?yàn)楫?dāng)結(jié)構(gòu)處在非線性階段時(shí)，其剛度矩陣隨時(shí)間不斷變化，這就需要我們不斷地計(jì)算剛度矩陣和求特征向量，尤其是求特征向量將會(huì)耗去大量的機(jī)時(shí)，極不經(jīng)濟(jì)。然而，振型疊加法也有其優(yōu)點(diǎn)，它可以將運(yùn)動(dòng)方程化為振型坐標(biāo)系中一組單自由度運(yùn)動(dòng)方程來(lái)求解，并且不必列出結(jié)構(gòu)所有自由度的卑自由度方程，而只列出對(duì)結(jié)構(gòu)反應(yīng)起主要作用的前若干階單自由度方程進(jìn)行求解，這種把結(jié)構(gòu)總反應(yīng)看作是若干個(gè)振型反應(yīng)的疊加有助于人們對(duì)地震反應(yīng)實(shí)質(zhì)的理解。對(duì)于局部非線性問(wèn)題，結(jié)構(gòu)剛度矩陣只有某些位置的元素發(fā)生變化，而大部分位置的元素并未發(fā)生變化。因此如果我們能夠通過(guò)某種計(jì)算技巧，利用振型疊加法求結(jié)構(gòu)的時(shí)程反應(yīng)，并避免剛度矩陣的重新求特征值，從而解出非線性運(yùn)動(dòng)方程，此時(shí)振型疊加法就不失為一種比較好的計(jì)算方法。
本文試圖在研究阻尼選取問(wèn)題的基礎(chǔ)上，介紹非線性振型疊加法一些已有的研究成果，希望對(duì)抗震研究者求解橋梁局部非線性問(wèn)題有所幫助。
二、阻尼問(wèn)題
1．阻尼模型
結(jié)構(gòu)阻尼是對(duì)振動(dòng)結(jié)構(gòu)所耗散的能量的測(cè)量，通常用振動(dòng)一次的能量耗散率來(lái)表示結(jié)構(gòu)阻尼的強(qiáng)弱。近幾十年來(lái)，人們提出了多種阻尼理論假設(shè)，在眾多的阻尼理論假設(shè)中，用得較多的是兩種線性阻尼理論：粘滯阻尼理論和復(fù)阻尼理論（滯變阻尼理論）。
復(fù)阻尼理論認(rèn)為結(jié)構(gòu)具有復(fù)剛度，在考慮阻尼時(shí)在彈性模量或剛度系數(shù)項(xiàng)前乘以復(fù)常數(shù) 即可，v為復(fù)阻尼系數(shù)。復(fù)阻尼理論對(duì)于一般的結(jié)構(gòu)動(dòng)力響應(yīng)來(lái)說(shuō)，計(jì)算過(guò)程非常復(fù)雜，因此，在動(dòng)力響應(yīng)分析中，復(fù)阻尼理論應(yīng)用不多，本文限于篇幅，也就不再展開了。
粘滯阻尼理論假定阻尼力與運(yùn)動(dòng)速度成正比，通常是用不同頻率的阻尼比ζ來(lái)表征系統(tǒng)的阻尼：
粘滯阻尼理論最顯著的特點(diǎn)在于其阻尼力是直接根據(jù)與相對(duì)速度成正比的關(guān)系給出的，不論是簡(jiǎn)諧振動(dòng)或是非簡(jiǎn)諧振動(dòng)，都可直接寫出系統(tǒng)的運(yùn)動(dòng)方程，而且均為線性微分方程，給理論分析帶來(lái)了很大的方便。
在多自由度系統(tǒng)中采用等效粘滯模態(tài)阻尼，阻尼力向量的表達(dá)式為
若［C」可以通過(guò)模態(tài)向量正交化為對(duì)角矩陣時(shí)，則稱為正交阻尼或比例阻尼。反之，則稱之為非正交阻尼。正交阻尼原則上適用于阻尼特性分布比較均勻的工程結(jié)構(gòu)，但由于其使用方便，分析人員對(duì)大部分橋梁都傾向于使用正交阻尼，非正交阻尼因?yàn)橛?jì)算較為麻煩用得較少。
Rayleigh阻尼模型是廣泛采用的一種正交阻尼模型，其數(shù)學(xué)表達(dá)式如下：
C＝a0M＋a1K （2）
式中， a0和a1稱為Rayleigh阻尼常數(shù)。
在Rayleigh阻尼模型下，各階阻尼比可表示為
式中ζi稱為第i階振型的模態(tài)阻尼比，因此若已知任意兩階振型的阻尼比ζi和ζj，則可定出阻尼常數(shù)
確定了a0和al之后，即可確定出各階振型的模態(tài)阻尼比，并確定阻尼矩陣。
2．實(shí)際抗震分析中由于阻尼選取不同所產(chǎn)生的問(wèn)題
目前，橋梁地震反應(yīng)分析一般以直接積分的時(shí)程分析方法為主。其阻尼模型取Rayleigh阻尼模型，并以主塔或主梁的兩個(gè)較低階振型頻率ωi和ωj對(duì)應(yīng)的阻尼比作為ζi和ζj，接式（3）和式(4)
求出其余各階頻率的阻尼比，并求出阻尼矩陣代人動(dòng)力方程，用直接積分的方法求解動(dòng)力方程。這樣處理阻尼雖然非常簡(jiǎn)單，但也產(chǎn)生了以下兩個(gè)不可忽視的問(wèn)題：
（1）如前所述，Rayleigh阻尼作為一種正交阻尼，適用于阻尼特性分布非常均勻的工程結(jié)構(gòu)。但是大跨橋梁一般來(lái)說(shuō)都不能算作非常均勻的結(jié)構(gòu)。例如，為了提高橋梁的跨越能力，主梁一般采用鋼箱梁或鋼混疊合梁，而主塔和邊墩則采用鋼筋混凝土材料，兩者的阻尼特性相差比較大。即使主梁材料特性與主塔差不多，大跨橋梁由于抗風(fēng)和抗震的要求，經(jīng)常會(huì)在橋梁結(jié)構(gòu)的某些部位加有人工阻尼裝置，比如橋墩上安放高阻尼的抗震支座、橋塔上安放控制振動(dòng)的裝置TMD等，這都會(huì)產(chǎn)生摩擦阻尼或集中阻尼從而造成阻尼特性的不均勻分布。這樣的阻尼均勻性前提得不到滿足的情況下，仍按照
Rayleigh阻尼模型去計(jì)算各階振型對(duì)應(yīng)的阻尼比勢(shì)必會(huì)造成除ωi和ωj兩階之外其他各階振型阻尼比與真實(shí)值有或多或少的差別。
（2）根據(jù)同濟(jì)大學(xué)土木防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室對(duì)國(guó)內(nèi)幾十座大跨橋梁進(jìn)行抗震分析后總結(jié)的經(jīng)驗(yàn)，邊墩。輔助墩等部位是大跨橋梁抗震設(shè)施的重點(diǎn)。但是采用Rayleigh阻尼模型時(shí)，用于計(jì)算其他各階振型阻尼比的ωi和ωj一般取的是較低階的振型，而邊墩輔助墩的振動(dòng)一般都發(fā)生在高階振型。根據(jù)Rayleigh阻尼模型圖，可以看出離ωi和ωj越遠(yuǎn)的振型，其阻尼比就越不準(zhǔn)，而且隨著圖上阻尼比按頻率增加的速度越來(lái)越快，邊墩部分振動(dòng)頻率對(duì)應(yīng)的阻尼比比實(shí)際值往往偏大，從這一點(diǎn)講會(huì)導(dǎo)致邊墩部分反應(yīng)的計(jì)算結(jié)果偏于不安全。
一些橋梁抗震研究人員已經(jīng)注意到了以上兩個(gè)問(wèn)題，他們采取的措施是根據(jù)分析的部位不斷變換所選擇的ωi和ωj，比如計(jì)算橋塔的縱向地震反應(yīng)時(shí)就選擇對(duì)橋塔的縱向反應(yīng)起主要作用的兩階頻率作為ωi和ωj，來(lái)計(jì)算其它各階阻尼比，計(jì)算其它地震反應(yīng)時(shí)也依此類推。這樣就需要分析人員不斷的重復(fù)選擇。和約和進(jìn)行時(shí)程計(jì)算，十分繁瑣。
3．解決方法
由以上論述，我們已經(jīng)了解到阻尼是一個(gè)非常復(fù)雜的問(wèn)題，僅僅依靠Rayleigh阻尼模型，會(huì)對(duì)大跨橋梁尤其是邊墩輔助墩等部位的地震反應(yīng)分析出現(xiàn)不應(yīng)有的誤差。因此，我們嘗試尋找一種既不過(guò)分繁瑣又比較準(zhǔn)確的方法。
在前面的論述中，我們發(fā)現(xiàn)阻尼比是反應(yīng)阻尼的一個(gè)方便而有效的量，它把阻尼特性和振型頻率聯(lián)系起來(lái)，使得動(dòng)力方程分析起來(lái)更為簡(jiǎn)單，而且阻尼比可以通過(guò)橋梁實(shí)測(cè)測(cè)出。
如果我們直接指定對(duì)橋塔。主梁、邊墩等重要部位反應(yīng)起主要作用的一些振型頻率的阻尼比，而對(duì)其余各階振型頻率的阻尼比采用線性內(nèi)插的方法確定，這樣做也可以形成阻尼比矩陣。由于我們通過(guò)以前的工程實(shí)例發(fā)現(xiàn)結(jié)構(gòu)各部位的反應(yīng)來(lái)說(shuō)少數(shù)幾階振型的貢獻(xiàn)最為顯著（這些振型的貢獻(xiàn)占到70％～
80％，甚至更多），因此，這樣做能夠保證計(jì)算的正確性，而且并不繁瑣，此對(duì)，以實(shí)測(cè)試驗(yàn)數(shù)據(jù)作為基礎(chǔ)，更增加了其準(zhǔn)確性。同濟(jì)大學(xué)橋梁系近十幾年來(lái)，通過(guò)為國(guó)內(nèi)幾十座大型橋梁進(jìn)行竣工檢測(cè)、成橋檢測(cè)積累了大量的阻尼實(shí)測(cè)資料，并有研究人員準(zhǔn)備把這些阻尼資料整理形成橋梁阻尼數(shù)據(jù)庫(kù)。有了這些數(shù)據(jù)資料為基礎(chǔ)，通過(guò)指定主要振型頻率阻尼比，來(lái)計(jì)算結(jié)構(gòu)動(dòng)力反應(yīng)是行得通的，并且結(jié)合下面的振型疊加法，會(huì)使計(jì)算更加簡(jiǎn)便。
三、非線性問(wèn)題的振型疊加法
1．線性問(wèn)題的振型疊加法
振型疊加法是利用多自由度系統(tǒng)的固有頻率和振型的特性，將結(jié)構(gòu)動(dòng)力響應(yīng)分解為各個(gè)振型分量，對(duì)各個(gè)振型分量分別求解后疊加得到實(shí)際的響應(yīng)。它假定阻尼滿足Rayleigh阻尼模型，利用正交變換將線性動(dòng)力分析中多自由度體系相互耦合的N個(gè)方程
轉(zhuǎn)換成為N個(gè)不耦合的單自由度方程
這些單自由度運(yùn)動(dòng)方程可以通過(guò)杜哈美積分的數(shù)值算法解出。所以，總反應(yīng)為
其中，[Φ]為 N×N的振型向量矩陣。
如果對(duì)結(jié)構(gòu)反應(yīng)起主要作用的是較低的P階振型，那我們就可以只計(jì)算較低的P階單自由度方程，根據(jù)計(jì)算機(jī)數(shù)值計(jì)算的特點(diǎn)并考慮截去的高階振型的影響后，最后實(shí)用的計(jì)算式為[2]
或
其中［Kp」--對(duì)應(yīng)于P個(gè)特征向量的結(jié)構(gòu)剛度矩陣：
{Fp(t)}--對(duì)應(yīng)于P個(gè)特征向量的荷載向量
2．非線性問(wèn)題的振型疊加法
非線性振型疊加法是我們著重要討論的問(wèn)題。在線性問(wèn)題的振型疊加法中，因?yàn)閯偠染仃嚥浑S時(shí)間變化，所以只需要時(shí)程計(jì)算的開始進(jìn)行一次求特征向量的計(jì)算，以后就利用求出的這組特征向量通過(guò)選代算法求解結(jié)構(gòu)各個(gè)時(shí)刻的反應(yīng)。而在非線性問(wèn)題中，結(jié)構(gòu)的剛度矩陣是隨時(shí)間不斷變化的，這就需要在整個(gè)時(shí)程分析過(guò)程中不斷地計(jì)算剛度矩陣并求解特征向量。求矩陣的特征向量是非常耗費(fèi)機(jī)時(shí)的計(jì)算，因而這樣運(yùn)用振型疊加法將會(huì)耗費(fèi)大量的機(jī)時(shí)，極不經(jīng)濟(jì)。但是如果我們能夠避開多次求解特征向量，而只在計(jì)算的開始求一次，然后利用某種計(jì)算技巧，求出各時(shí)刻的反應(yīng)，這樣，振型疊加法對(duì)求解非線性問(wèn)題特別是對(duì)大跨橋梁地震反應(yīng)這種剛度矩陣在整個(gè)時(shí)間過(guò)程中只是部分元素發(fā)生變化而大部分元素不變的局部非線性問(wèn)題還是有應(yīng)用價(jià)值的。下面本文以材料非線性問(wèn)題為例來(lái)說(shuō)明這種分析方法的過(guò)程［1］
受地震加速度作用的非線性體系的運(yùn)動(dòng)方程為
為便于迭代求解，將上式進(jìn)行轉(zhuǎn)化，得到
其中，[k]為0時(shí)刻結(jié)構(gòu)對(duì)應(yīng)的剛度矩陣。{R(u(t))}是當(dāng)考慮材料非線性時(shí)與材料為彈性所產(chǎn)生的恢復(fù)力的差。
與線性問(wèn)題的振型疊加法類似，將相互耦合的N×N動(dòng)力方程化為在正規(guī)坐標(biāo)系下的N個(gè)單自由度方程：
則單自由度運(yùn)動(dòng)方程變?yōu)?br>
其中
上面的方程根據(jù)Nau的研究成果[3]，解可寫為
其中，[A]和[B]兩矩陣的計(jì)算方法可從文獻(xiàn)「3」查得。
我們只需求解對(duì)結(jié)構(gòu)反應(yīng)起主要影響的較低的P階單自由度方程，在考慮截去的高階振型的影響后，結(jié)構(gòu)總的反應(yīng)可表示為
其中
Fr(t,u)為廣義力：
式（20）～（22）即為結(jié)構(gòu)t時(shí)刻總反應(yīng)的計(jì)算式，該式中{R}與待求量結(jié)構(gòu)變形u(t)有關(guān)，
S狀態(tài)代表t時(shí)刻， s＋l狀態(tài)代表 t十Δt時(shí)刻。因此，求解過(guò)程需進(jìn)行選代。具體迭代算法如下面所示。
3．計(jì)算步驟
（l）設(shè)置初始條件：
(5)根據(jù)結(jié)構(gòu)位移向量計(jì)算該時(shí)刻結(jié)構(gòu)單元應(yīng)力。
(7)重復(fù)以上步驟，計(jì)算下一個(gè)時(shí)刻的位移向量{u}（返回第（4）步）。
根據(jù)文獻(xiàn)[1]的研究結(jié)果，使用非線性疊加法的程序與DRAIN2d程序的直接積分法分析同樣的平面框架的地震反應(yīng)，在誤差為6％的情況下，后者所花時(shí)間平均為前者的兩倍。
四、結(jié)語(yǔ)
本文介紹了阻尼和非線性振型疊加法兩個(gè)方面，目的是針對(duì)大跨橋梁局部非線性問(wèn)題探索一種簡(jiǎn)單前效名用的地震反應(yīng)計(jì)算方法。振型疊加法用于結(jié)構(gòu)局部非線性地震反應(yīng)計(jì)算還不很成熟，還有一些具體問(wèn)題，比如阻尼數(shù)據(jù)編輯整理、迭代方法的穩(wěn)定性研究、編制應(yīng)的分析程序并與直接積分法程序比較計(jì)算效果等有待于進(jìn)一步地探索和解決。
參考文獻(xiàn)
[1] Melod Michel hanna, "An Efficient Mode Superposition
method for the Numerical Dynamic Analysis of Bilinear
Systems",Dissertation of Ph. D, University of
California, Irvine, 1989
[2] Wildson, E, L., and Leger, P., "Modal Summation
methods for Structural Dynamic computations,"
Earthquake Engineering and Structral Dynamics, Vol.
16, 23 ~ 27, 1988
[3] Nau, J. M., "Computation of Inelastic Response
Spectra," J. Of Engineering mechanics, ASCE,
Vol.109, No. l, 279~ 288, Feb., 1983
