李小康 謝壯寧 王湛

(華南理工大學(xué)亞熱帶建筑科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，廣東廣州510640)

作為一種能量耗散裝置，調(diào)諧質(zhì)量阻尼器(Tuned Mass Damper，TMD)可用來(lái)增加結(jié)構(gòu)的總阻尼從而抑制結(jié)構(gòu)振動(dòng)，目前在大型復(fù)雜土木結(jié)構(gòu)的風(fēng)振控制中獲得較為廣泛的應(yīng)用.TMD的設(shè)計(jì)制作有賴(lài)于對(duì) TMD風(fēng)振控制效果的精確計(jì)算，因此TMD風(fēng)振控制計(jì)算引起了結(jié)構(gòu)風(fēng)工程界的廣泛關(guān)注［1-6］.目前，對(duì)高層建筑TMD風(fēng)振控制機(jī)理與計(jì)算方法的研究已趨成熟，但相關(guān)工作主要集中在順風(fēng)向控制，針對(duì)橫風(fēng)向風(fēng)振控制的研究較少.盡管就控制機(jī)理本身，橫、順風(fēng)向并無(wú)本質(zhì)差別，但由于荷載形成機(jī)理的差異導(dǎo)致橫風(fēng)向氣動(dòng)荷載要遠(yuǎn)比順風(fēng)向復(fù)雜，目前尚無(wú)法類(lèi)似于順風(fēng)向給出統(tǒng)一的表達(dá)形式.風(fēng)洞試驗(yàn)作為唯一可精確獲取結(jié)構(gòu)橫風(fēng)向氣動(dòng)荷載的有效手段被廣泛應(yīng)用，但針對(duì)試驗(yàn)數(shù)據(jù)的風(fēng)振控制計(jì)算分析則少有研究.隨著結(jié)構(gòu)高度超過(guò)200m，結(jié)構(gòu)橫風(fēng)向響應(yīng)逐漸超越順風(fēng)向成為控制風(fēng)向，因此對(duì)橫風(fēng)向風(fēng)振控制的計(jì)算研究更具實(shí)際意義和工程價(jià)值.

針對(duì)任意復(fù)雜形式的大型空間結(jié)構(gòu)，文中提出一種多調(diào)諧質(zhì)量阻尼器(MTMD)風(fēng)振控制的快速算法，該方法計(jì)算效率高、內(nèi)存消耗少，屬于一種精確算法，使得需要大量計(jì)算的單調(diào)諧質(zhì)量阻尼器(STMD)參數(shù)分析成為可能，同時(shí)也可進(jìn)行更為復(fù)雜的 MTMD風(fēng)振控制的精確計(jì)算分析.文中以439m高的深圳京基金融中心(KFT)為背景進(jìn)行風(fēng)致響應(yīng)及控制的計(jì)算，分析了STMD的頻率、阻尼比和質(zhì)量等參數(shù)對(duì)風(fēng)振控制效果的影響，并對(duì)不同形式MTMD的控制效果進(jìn)行詳細(xì)的分析計(jì)算和討論.

1 方法概述

受l個(gè)TMD控制的高層建筑，其風(fēng)振反應(yīng)的運(yùn)動(dòng)方程可表示為

式中:M、C、K分別是結(jié)構(gòu)系統(tǒng)的質(zhì)量、阻尼及剛度矩陣;MT、CT、KT分別表示TMD系統(tǒng)的質(zhì)量、阻尼及剛度矩陣;B為N×l階TMD系統(tǒng)的位置指示矩陣，bji=δjk(i=1，2，…，N);x為結(jié)構(gòu)位移向量;w為T(mén)MD相對(duì)位移向量;p為氣動(dòng)力向量.

令x=Φq，其中Φ表示結(jié)構(gòu)模態(tài)矩陣(N×r階)，q為前r階模態(tài)響應(yīng)向量，由振型分解法可得受控結(jié)構(gòu)在模態(tài)坐標(biāo)下的運(yùn)動(dòng)方程:

式中，ζs、ωs分別為前r階結(jié)構(gòu)模態(tài)阻尼比和模態(tài)頻率矩陣，ζT、ωT分別為l階TMD阻尼比和頻率矩陣(均為對(duì)角陣)，g=ΦTp為廣義力向量，，Mp=ΦTMΦ為模態(tài)質(zhì)量矩陣.

對(duì)式(2)進(jìn)行簡(jiǎn)單處理，TMD運(yùn)動(dòng)方程可寫(xiě)為

對(duì)于式(3)，由于方程右邊第二項(xiàng)的存在，使得受TMD控制下的結(jié)構(gòu)運(yùn)動(dòng)方程不能簡(jiǎn)單地轉(zhuǎn)化至頻域內(nèi)進(jìn)行求解，常用的辦法是將式(3)中方程右邊的第二項(xiàng)移動(dòng)至方程的左邊，將變化后的式(3)與式(4)共同組成另一方程組進(jìn)行求解，但當(dāng)式(3)中結(jié)構(gòu)多階模態(tài)參與響應(yīng)貢獻(xiàn)時(shí)，組裝后方程組再次出現(xiàn)需要解耦的問(wèn)題.盡管目前對(duì)于高層建筑的風(fēng)振控制計(jì)算可進(jìn)行適當(dāng)簡(jiǎn)化，采用一階模態(tài)響應(yīng)近似估算結(jié)構(gòu)總響應(yīng)［4］，但這種處理會(huì)影響其普適性，尤其是對(duì)于大型復(fù)雜空間結(jié)構(gòu)(也包括考慮橫、順、扭耦合作用的超高層建筑)，僅考慮一階模態(tài)顯然是不夠的.

采用文獻(xiàn)［7-8］中的基本思想推導(dǎo)結(jié)構(gòu)在受TMD控制下的結(jié)構(gòu)運(yùn)動(dòng)方程.假定pu、qu、gu、xu和wu分別為氣動(dòng)力p、模態(tài)響應(yīng)q、廣義力g、響應(yīng)x和w的截尾函數(shù)矢量，同樣xu=Φqu，則有

由于截尾函數(shù)矢量pu滿(mǎn)足絕對(duì)可積條件，可對(duì)其進(jìn)行傅里葉變換，有

式中，i為虛數(shù)單位;t、ω分別表示時(shí)間和圓頻率變量，Pu是ω和截尾函數(shù)的截取時(shí)間長(zhǎng)度T0的函數(shù).同樣，Gu、Qu和Wu分別表示gu、qu和wu的傅里葉變換，其表達(dá)式類(lèi)似于式(7).可以證明［9］，當(dāng)p為平穩(wěn)隨機(jī)過(guò)程時(shí)，其功率譜密度矩陣Sp為

式中，上標(biāo)*表示共軛.由此，式(5)、(6)在頻域內(nèi)可表示為

記

hs=(-ω2+2ζsωsωi)-1(為r階對(duì)角陣，表示無(wú)TMD控制的結(jié)構(gòu)廣義坐標(biāo)的頻響函數(shù)矩陣)，

由式(10)得

代入式(9)，有

經(jīng)分析可知，式(12)中括號(hào)項(xiàng)的逆陣是存在的，記

它實(shí)際為T(mén)MD控制下的結(jié)構(gòu)主坐標(biāo)頻響函數(shù)矩陣，由此式(12)可簡(jiǎn)述為

由式(7)及gu定義可知:

同時(shí)令

則TMD控制下的結(jié)構(gòu)風(fēng)致響應(yīng)x的功率譜函數(shù)矩陣Sx可采用下式計(jì)算:

對(duì)于TMD系統(tǒng)本身，將式(15)代入式(11)，并定義

則TMD系統(tǒng)響應(yīng)表示為

分別對(duì)結(jié)構(gòu)系統(tǒng)和TMD系統(tǒng)的功率譜函數(shù)進(jìn)行積分，即可獲得響應(yīng)協(xié)方差矩陣:同樣其功率譜函數(shù)矩陣可采用下式進(jìn)行計(jì)算:

由上可以看出，在無(wú)TMD控制時(shí)，方程(1)中位置矩陣B為零矩陣，方程(1)退化為無(wú)TMD控制的運(yùn)動(dòng)方程，可見(jiàn)文中算法適用于任意復(fù)雜形式的結(jié)構(gòu)風(fēng)致響應(yīng)及TMD控制計(jì)算，屬于對(duì)HEM算法［7］的擴(kuò)展.同時(shí)注意到有無(wú)TMD控制的結(jié)構(gòu)響應(yīng)計(jì)算僅在于系統(tǒng)頻響函數(shù)之間的區(qū)別，與無(wú)TMD控制的結(jié)構(gòu)頻響函數(shù)矩陣hs相比，受TMD控制的結(jié)構(gòu)頻響函數(shù)矩陣hT并不是對(duì)角陣，但作為一確定性矩陣而言對(duì)計(jì)算效率基本不存在影響.

特別需要指出的是，如果采用傳統(tǒng)完全二次型組合方法進(jìn)行風(fēng)振控制計(jì)算，仍然會(huì)出現(xiàn)3個(gè)3維矩陣相乘的問(wèn)題，其譜矩陣的存儲(chǔ)需要非常大的內(nèi)存空間，計(jì)算也頗為耗時(shí)，即便可勉強(qiáng)計(jì)算下去，計(jì)算效率也令人難以忍受［8］.

2 實(shí)例驗(yàn)證及分析

2.1 KFT結(jié)構(gòu)特征及風(fēng)洞試驗(yàn)簡(jiǎn)介

KFT工程的基本建筑外形如圖1所示，塔樓主立面(南北立面)基本為等寬度圓弧收縮曲面，兩邊有1.8m懸挑翼板.圖1(b)中西側(cè)立面(東西兩側(cè)立面相同)為圓錐形平面，立面上設(shè)置遮陽(yáng)構(gòu)件，塔樓頂部高度為439m，從上至下依次設(shè)置5個(gè)設(shè)備層和避難層，其編號(hào)及高度如圖中標(biāo)注.圖1(c)給出了塔樓基本平面布置及風(fēng)洞試驗(yàn)坐標(biāo)系.KFT模型采用有機(jī)玻璃制作，幾何縮尺比為1∶450，模型總高度為0.976m，試驗(yàn)詳細(xì)信息可見(jiàn)文獻(xiàn)［10］.KFT風(fēng)致響應(yīng)計(jì)算中取結(jié)構(gòu)前12階模態(tài)進(jìn)行響應(yīng)計(jì)算，結(jié)構(gòu)前3階模態(tài)周期分別為7.52、7.19和3.92s，分別對(duì)應(yīng)x、y和扭轉(zhuǎn)方向，結(jié)構(gòu)的模態(tài)阻尼比分別取為3.5%和1.5%，用于位移和結(jié)構(gòu)頂部加速度計(jì)算.

2.2 STMD參數(shù)分析

在結(jié)果分析之前，定義以下基于頂部加速度的減振系數(shù)βa和頂部位移的減振系數(shù)βd:

圖1 建筑示意圖及參考坐標(biāo)系(單位:m)Fig.1 Building drawing and definition of coordinate system (Unit:m)

式中，σa，T和σa分別為10年重現(xiàn)期受控和不受控的結(jié)構(gòu)頂部加速度標(biāo)準(zhǔn)差，σd，T和σd為100年重現(xiàn)期受控和不受控的結(jié)構(gòu)頂部位移響應(yīng)標(biāo)準(zhǔn)差.由定義可知，βa和βd越小，減振效果越好.為表述簡(jiǎn)潔，下文出現(xiàn)的加速度響應(yīng)均為10年重現(xiàn)期結(jié)果，位移響應(yīng)均為100年重現(xiàn)期結(jié)果.

首先考慮STMD情況，對(duì)受控情況下的結(jié)構(gòu)風(fēng)致響應(yīng)進(jìn)行計(jì)算，分析STMD各參數(shù)對(duì)控制效果的影響.結(jié)果包含結(jié)構(gòu)系統(tǒng)各自由度位移、加速度及等效靜風(fēng)荷載響應(yīng)和TMD系統(tǒng)的位移、加速度響應(yīng)等，限于篇幅，文中僅給出橫風(fēng)向響應(yīng)結(jié)果.

在頻率比固定不變的情況下，首先考察ξT對(duì)控制效果的影響，結(jié)果如圖2所示，其中質(zhì)量比UT分別取值為0.005、0.010、0.015，對(duì)應(yīng)實(shí)際質(zhì)量為520、1039、1559t，頻率比f(wàn)T/f0=1(STMD頻率fT與結(jié)構(gòu)受控橫方向第1階固有頻率f0之比，f0=0.139)，ξT取值范圍為0～0.1.

圖2 減振系數(shù)隨STMD阻尼比的變化Fig.2 Vibration reduction of coefficient with.STMD damping ratio

由圖2可見(jiàn):隨ξT增加，βd呈現(xiàn)先減小后增大的趨勢(shì);當(dāng)ξT太小以致接近零時(shí)，控制效果最差;存在一個(gè)ξT最優(yōu)值，當(dāng)ξT超過(guò)最優(yōu)值繼續(xù)增大時(shí)，位移的控制效果反而降低，但變化幅度并不顯著，特別是隨著UT的增加，反映控制效果的參數(shù)βd變化曲線趨于平緩，且由圖2可以看出，隨著UT的增加，ξT最優(yōu)值也逐漸增加，當(dāng)UT取0.015時(shí)，最優(yōu)的ξT接近0.06.

在保證一定控制效果的條件下，實(shí)際TMD設(shè)計(jì)中宜取較大的ξT以保證TMD本身有足夠的能量耗散性能而不至于由于過(guò)度的振蕩產(chǎn)生疲勞破壞而使TMD失效，當(dāng)然TMD本身還會(huì)有限位裝置防止這種情況發(fā)生.因此根據(jù)圖2，在文中后續(xù)的計(jì)算分析中，均按照ξT=0.065進(jìn)行計(jì)算.

圖3給出了STMD控制下結(jié)構(gòu)均方根位移響應(yīng)隨STMD頻率比f(wàn)T/f0的變化情況，fT/f0取值范圍為0.1～2.0.

圖3 減振系數(shù)隨STMD頻率比的變化Fig.3 Vibration of reduction coefficient with STMD frequency ratio

由圖3可見(jiàn):(1)隨UT的增加，βd逐漸減小，STMD的控制效果增加;UT=0.015時(shí)，βd最小為0.723，但在實(shí)際工程中多方面因素限制了TMD質(zhì)量不可能過(guò)大，同時(shí)過(guò)度增加質(zhì)量會(huì)影響整體結(jié)構(gòu)的動(dòng)力參數(shù)指標(biāo)，使得計(jì)算結(jié)果失真.(2)在ξT不變的前提下，隨著fT遠(yuǎn)離結(jié)構(gòu)固有頻率值，βd增加，控制效果減弱;fT在結(jié)構(gòu)固有頻率附近時(shí)，控制效果最佳，但由于阻尼的影響，最佳控制效果的fT略小于結(jié)構(gòu)固有頻率f0.(3)βd對(duì)fT的變化較為敏感，當(dāng)fT/f0小于0.5或大于1.5時(shí)，βd均接近1;STMD的控制效果僅在較小的fT范圍內(nèi)較好.

圖4給出了結(jié)構(gòu)受控前后結(jié)構(gòu)頂部位移響應(yīng)的功率譜密度在不同坐標(biāo)系統(tǒng)下的比較結(jié)果，為了直觀起見(jiàn)，圖中功率譜密度的變量改用頻率f表示，豎線表示結(jié)構(gòu)第2階固有頻率為0.139 Hz(對(duì)應(yīng)結(jié)構(gòu)橫風(fēng)向第1階頻率).由圖可見(jiàn)，橫風(fēng)向響應(yīng)的共振部分遠(yuǎn)大于背景部分，在TMD作用下，響應(yīng)的共振部分被抑制，處于共振區(qū)(結(jié)構(gòu)固有頻率附近范圍)的響應(yīng)功率譜密度函數(shù)值急劇減小，使得結(jié)構(gòu)風(fēng)致響應(yīng)顯著降低.由圖4可以清晰地看出，TMD裝置只能控制其對(duì)應(yīng)模態(tài)頻率附近的功率譜密度值.

圖4 結(jié)構(gòu)頂層位移響應(yīng)功率譜密度函數(shù)Fig.4 PSD of displacement response at the top floor of the building

2.3 MTMD參數(shù)分析

2.3.1 不同樓層MTMD的控制效果

由于超高層建筑的風(fēng)振響應(yīng)一般是一階振型控制的，故在理論上將STMD置于結(jié)構(gòu)頂層可以獲取最佳的控制效果，但實(shí)際應(yīng)用中在某些情況下，受空間、安裝和結(jié)構(gòu)等方面的限制，必須采用每個(gè)TMD相比更小的MTMD方式.

以KFT工程為例，改用3個(gè)質(zhì)量相同的TMD，安裝在在1、2、3號(hào)設(shè)備層(見(jiàn)圖1)上，TMD的總質(zhì)量比UT仍然取為0.01，每個(gè)TMD的fT/f0=1，ξT= 0.065，采用文中方法計(jì)算其控制效果，結(jié)果見(jiàn)表1.注意比較是基于TMD總質(zhì)量一致的原則也即表中的3種控制策略的TMD總質(zhì)量比是一致的，均為0.01.由表1可見(jiàn)，僅在設(shè)備層1安裝TMD的控制效果最佳(即STMD)，βa=0.686，βd=0.778，在3個(gè)設(shè)備層均設(shè)置有TMD時(shí)，βa=0.734，βd=0.820，加速度和位移的控制效果比單個(gè)TMD方案差些，但差別不大，只有0.048和0.042.

表1 不同樓層MTMD作用下結(jié)構(gòu)頂部的橫風(fēng)向響應(yīng)Table 1 Across-wind response at top with MTMD on different fcoors

雖然MTMD的控制效果不如STMD，但MTMD系統(tǒng)中每個(gè)TMD的體積和重量相對(duì)較小，更易于安裝、實(shí)現(xiàn)，從整個(gè)系統(tǒng)上考慮，其可靠度也更好，當(dāng)MTMD中某個(gè)TMD出現(xiàn)故障時(shí)，其它TMD依然可保持工作狀態(tài)且能夠取得一定的控制效果.考慮以上3個(gè)小TMD構(gòu)成MTMD，在1、2、3號(hào)設(shè)備層布置3個(gè)相同的TMD，計(jì)算不同工作狀態(tài)時(shí)的控制效果見(jiàn)表2.其中每個(gè)TMD的參數(shù)為:ξT=0.065，UT= 0.0033，fT/f0=1，工況欄中的數(shù)字表示參加工作的TMD所在的設(shè)備層號(hào)(沒(méi)出現(xiàn)號(hào)碼表示該TMD失效)，如1/2表示只有1、2設(shè)備層的TMD參加工作而3號(hào)設(shè)備層TMD失效.

表2 不同工況MTMD作用下結(jié)構(gòu)頂部橫風(fēng)向響應(yīng)Table 2 Across-wind response at top with MTMD in different experimental cases

計(jì)算結(jié)果表明，盡管MTMD中單個(gè)TMD出現(xiàn)故障導(dǎo)致控制效果有所減弱，但其他層TMD仍可起到較好的控制作用.例如，第3設(shè)備層的TMD失效時(shí)，βa=0.754，βd=0.840;即便是控制效果最好的第1設(shè)備層TMD失效，βa和βd仍有0.819和0.889，較無(wú)控制工況也分別降低了0.181和0.111.

2.3.2 MTMD頻率的影響

MTMD系統(tǒng)中可通過(guò)對(duì)每個(gè)TMD設(shè)置不同的頻率進(jìn)一步擴(kuò)大其控制范圍.在KFT設(shè)備層1、2、3上布置不同頻率的TMD組成MTMD系統(tǒng)，第1、2和3設(shè)備層的TMD的fT/f0分別取為1.10、1.02和0.94，UT=0.0033，ξT=0.065.同時(shí)分析在不同結(jié)構(gòu)周期情況下對(duì)頂部加速度的控制效果并和STMD (fT/f0=1，UT=0.001，ξT=0.065)情況做比較，結(jié)果見(jiàn)圖5.

圖5 減振系數(shù)隨結(jié)構(gòu)周期的變化Fig.5 Vibration of reduction coefficient with structural natural period

由圖5可見(jiàn)，當(dāng)結(jié)構(gòu)周期比在0.90～1.10之間變化時(shí)，總的看來(lái)兩種控制方式的效果隨周期比變化的趨勢(shì)是一樣的，但MTMD的控制效果要略差于STMD的控制效果.

2.3.3 受控結(jié)構(gòu)阻尼的影響

實(shí)際上，結(jié)構(gòu)參數(shù)中阻尼的不確定性會(huì)遠(yuǎn)高于結(jié)構(gòu)頻率，假定結(jié)構(gòu)周期不變，分析以上STMD和MTMD在不同結(jié)構(gòu)阻尼比ζs下的結(jié)構(gòu)頂部加速度的控制效果，結(jié)果見(jiàn)圖6.

圖6 減振系數(shù)隨結(jié)構(gòu)阻尼比的變化Fig.6 Vibration of reduction coefficient with structural damping ratio

由圖6可見(jiàn)，在ζs＞0.003時(shí)，兩種控制方案的βa均隨ζs的增加而單調(diào)增加，控制效果在ζs=0.003時(shí)最好，但由于實(shí)際結(jié)構(gòu)的模態(tài)阻尼比不可能小于0.003，因而圖中ζs的較小值并沒(méi)有太大的實(shí)際借鑒意義，可以認(rèn)為結(jié)構(gòu)阻尼比越小控制效果越好.

2.3.4 不同頻率MTMD的全風(fēng)向控制結(jié)果

以上TMD控制方法的研究主要針對(duì)0°風(fēng)向角(橫風(fēng)向)進(jìn)行，實(shí)際工程中需要對(duì)36個(gè)風(fēng)向角的控制效果進(jìn)行綜合評(píng)判，以圖5中選用的MTMD系統(tǒng)為例，采用文中快速算法計(jì)算KFT工程36個(gè)風(fēng)向角的加速度響應(yīng)，并與無(wú)控制的結(jié)果進(jìn)行比較，結(jié)果見(jiàn)圖7.

圖7 全風(fēng)向角下頂部峰值加速度Fig.7 Peak acceleration of structure at all azimuths

由圖7可知，MTMD作用時(shí)0°和180°風(fēng)向角附近(橫風(fēng)向)結(jié)構(gòu)頂部峰值加速度比無(wú)控制時(shí)顯著降低，無(wú)控時(shí)頂部最大峰值加速度為0.0229g，出現(xiàn)在350°風(fēng)向角，主要是由于上游地王大廈的干擾引起［10］，由《高層建筑混凝土結(jié)構(gòu)技術(shù)規(guī)程》［11］可知，該峰值加速度值不能滿(mǎn)足規(guī)范中有關(guān)高層住宅舒適度的規(guī)定(不大于0.015 g)，僅滿(mǎn)足辦公要求(不大于0.025 g)，而采用MTMD控制之后，可將該風(fēng)向角峰值加速度進(jìn)一步降低至0.0156g，減少約32%，這也是受控結(jié)構(gòu)全風(fēng)向的最大峰值加速度.而在90°和180°風(fēng)向角附近(順風(fēng)向)，由于MTMD并不控制該方向風(fēng)致振動(dòng)，受控前后響應(yīng)差距不大.圖7還表明，在MTMD控制下，所有風(fēng)向角的結(jié)構(gòu)頂部峰值加速度響應(yīng)均比無(wú)控時(shí)要小.

3 結(jié)論

文中提出的多TMD風(fēng)振控制響應(yīng)的快速算法理論上屬于一種精確算法，且計(jì)算效率高，可對(duì)任意復(fù)雜形式的空間結(jié)構(gòu)進(jìn)行風(fēng)振控制計(jì)算.

單TMD的控制效果對(duì)TMD頻率比較敏感，最佳頻率比略小于1;最優(yōu)的TMD阻尼比和TMD質(zhì)量有關(guān)，隨TMD質(zhì)量的增大而增大，考慮到阻尼的耗能作用，建議采用0.065的TMD阻尼比.

在總質(zhì)量相同的條件下，采用在不同樓層設(shè)置TMD方式的減振效果只比在頂層設(shè)置TMD的效果減少5%左右，兩者的控制效果隨結(jié)構(gòu)頻率和阻尼的變化規(guī)律也保持一致.

采用MTMD方案的全風(fēng)向角計(jì)算結(jié)果表明，在MTMD的控制下，結(jié)構(gòu)頂部峰值加速度的最大值下降30%左右，基本滿(mǎn)足規(guī)范中高層住宅和辦公的舒適度要求.

［1］ Xu Y L，Kwok K C S，Samali B.The effect of tuned mass dampers and liquid dampers on cross-wind response of tall/slender structures［J］.Journal of Wind Engineering and Industrial Aerodynamic，1992，40(1):33-54.

［2］ Ricciardelli F.A linear model for structures with tuned mass dampers［J］.Wind and Structures，1999，2(3): 151-171.

［3］ Rana R，Soong T T.Parametric study and simplified design of tuned mass dampers［J］.Engineering Structures，1998，20(3):193-204.

［4］ 瞿偉廉，陶牟華，Chang C C.五種被動(dòng)動(dòng)力減振器對(duì)高層建筑脈動(dòng)風(fēng)振反應(yīng)控制的實(shí)用設(shè)計(jì)方法［J］.建筑結(jié)構(gòu)學(xué)報(bào)，2001，22(2):29-34. Qu Wei-lian，Tao Mu-hua，Chang C C.Practical design method for effect of five kinds of passive dynamic absorbers on fluctuation wind-induced vibration response control of tall buildings［J］.Journal of Building Structures，2001，22(2):29-34.

［5］ Li Chunxiang.Performance of multiple tuned mass dampers for attenuating undesirable oscillations of structures under the ground acceleration［J］.Earthquake Engineering and Structural Dynamics，2000，29(9):1405-1421.

［6］ 唐意，顧明.某超高層建筑TMD風(fēng)振控制分析［J］.振動(dòng)與沖擊，2006，25(1):16-19. Tang Yi，Gu Ming.Analysis on control of wind induced vibration of a super-tall building with TMD［J］.Journal of Vibration and Shock，2006，25(1):16-19.

［7］ 謝壯寧.風(fēng)致復(fù)雜結(jié)構(gòu)隨機(jī)振動(dòng)分析的一種快速算法——諧波激勵(lì)法［J］.應(yīng)用力學(xué)學(xué)報(bào)，2007，24(2): 263-266. Xie Zhuang-ning.New rapid algorithm for wind-induced random vibration of complex structures［J］.Chinese Jour-nal of Applied Mechanics，2007，24(2):263-266.

［8］ 李小康，謝壯寧.大跨度屋蓋結(jié)構(gòu)風(fēng)振響應(yīng)和等效靜風(fēng)荷載的快速算法和應(yīng)用［J］.土木工程學(xué)報(bào)，2010，43(7):29-36. Li Xiao-kang，Xie Zhuang-ning.Efficient algorithm and application for the wind-induced response and equivalent static wind load of large-span roof structures［J］.China Civil Engineering Journal，2010，43(7):29-36.

［9］ 楊位欽，顧嵐.時(shí)間序列分析與動(dòng)態(tài)數(shù)據(jù)建模［M］.北京:北京理工大學(xué)出版社，1988.

［10］ 謝壯寧，石碧青，倪振華，等.深圳京基金融中心氣動(dòng)抗風(fēng)措施試驗(yàn)研究［J］.建筑結(jié)構(gòu)學(xué)報(bào)，2010，31 (10):1-7. Xie Zhuang-ning，Shi Bi-qing，Ni Zhen-hua，et al.Experimental study on reduction of wind loads on Shenzhen Kingkey Financial Tower by aerodynamic strategy［J］. Journal of Building Structures，2010，31(10):1-7.

［11］ GB 50010—2010，高層建筑混凝土結(jié)構(gòu)技術(shù)規(guī)程［S］.