魯 正， 陳筱一， 王佃超， 呂西林

(1.同濟(jì)大學(xué) 土木工程防災(zāi)國家重點實驗室，上海 200092； 2.同濟(jì)大學(xué) 結(jié)構(gòu)工程與防災(zāi)研究所，上海 200092)

顆粒調(diào)諧質(zhì)量阻尼器減震控制的數(shù)值模擬

魯 正1,2， 陳筱一2， 王佃超2， 呂西林1,2

(1.同濟(jì)大學(xué) 土木工程防災(zāi)國家重點實驗室，上海 200092； 2.同濟(jì)大學(xué) 結(jié)構(gòu)工程與防災(zāi)研究所，上海 200092)

提出一種顆粒調(diào)諧質(zhì)量阻尼器，即通過擺繩將顆粒阻尼器懸掛于主體結(jié)構(gòu)上，將目前廣泛應(yīng)用的調(diào)諧質(zhì)量阻尼器與減震效能優(yōu)越的顆粒阻尼器結(jié)合起來，以擴(kuò)大減振頻帶，增加減振魯棒性?；谝欢ǖ牡刃г瓌t將多顆粒阻尼器等效為單顆粒阻尼器，提出一種近似的顆粒調(diào)諧質(zhì)量阻尼器數(shù)值模擬方法，并介紹了系統(tǒng)參數(shù)的取值方法。進(jìn)行了附加顆粒調(diào)諧質(zhì)量阻尼器的單自由度結(jié)構(gòu)實際地震輸入下的振動臺試驗，通過對比試驗結(jié)果與數(shù)值模擬結(jié)果，發(fā)現(xiàn)兩者吻合較好，說明該簡化方法可以較好地模擬顆粒調(diào)諧質(zhì)量阻尼器的減震控制效果。

顆粒調(diào)諧質(zhì)量阻尼器；調(diào)諧質(zhì)量阻尼器；顆粒阻尼器；數(shù)值模擬；結(jié)構(gòu)振動控制

調(diào)諧質(zhì)量阻尼器(Tuned Mass Damper,TMD)由質(zhì)量塊、彈簧和阻尼系統(tǒng)組成，通過調(diào)整TMD的自振頻率，使其接近主體結(jié)構(gòu)的被控頻率，當(dāng)主體結(jié)構(gòu)受迫振動時，TMD隨之產(chǎn)生與主體結(jié)構(gòu)相反的運動并施加反作用力于主體結(jié)構(gòu)，輸入結(jié)構(gòu)的振動能量通過調(diào)諧質(zhì)量阻尼器逐步消散，從而控制結(jié)構(gòu)的振動[1]。近年來，通過顆粒-結(jié)構(gòu)、顆粒-顆粒之間的非彈性碰撞與摩擦來實現(xiàn)動量交換及結(jié)構(gòu)能量耗散的顆粒阻尼器，克服了調(diào)諧質(zhì)量阻尼器減振頻帶窄的缺點，具有對原系統(tǒng)改動小、可靠度高、耐久性好、對溫度變化不敏感(在顆粒金屬熔點以下均可正常使用，鎢粉能承受近2 000 ℃高溫)、易于用在惡劣環(huán)境等優(yōu)點[2]，受到了廣大土木工程研究人員的重視，但目前在實際建筑結(jié)構(gòu)工程中的應(yīng)用甚少。本文嘗試將目前已成熟應(yīng)用于工程實際的TMD與顆粒阻尼技術(shù)結(jié)合起來，通過擺繩將填充顆粒群的容器懸掛于主體結(jié)構(gòu)上，形成一種新型的顆粒調(diào)諧質(zhì)量阻尼器，以增強(qiáng)其在地震激勵下的能量耗散，并擴(kuò)大其減振頻帶。

對于顆粒阻尼器在地震、風(fēng)振等荷載下顆粒行為的數(shù)值模擬始終是一個難點。考慮到多顆粒碰撞引起運動量的突變，使其動力學(xué)行為表現(xiàn)出很強(qiáng)的非線性，系統(tǒng)很難得到解析解，學(xué)者們研究出一些簡化方法和數(shù)值方法。PAPALOU等[3]將多顆粒阻尼器通過一定的等效原則簡化為單顆粒阻尼器；FRIEND等[4]提出一種解析方法，即把多顆粒模擬為一個凝聚的質(zhì)量塊，并將各種機(jī)理引起的能量耗散由一個通過實驗擬合得到的“有效恢復(fù)系數(shù)”表示；魯正等[5]基于離散單元法，建立了顆粒阻尼器對多自由度結(jié)構(gòu)進(jìn)行減震控制的數(shù)值模擬方法；閆維明等[6]建立了調(diào)諧型顆粒阻尼器的數(shù)值分析簡化模型，并提出了該模型等效阻尼比的能量估算方法。此外，學(xué)者們還提出了回歸模型法[7]和恢復(fù)力曲面法[8]等應(yīng)用于顆粒阻尼技術(shù)的數(shù)值模擬方法。

本文提出一種新型的顆粒調(diào)諧質(zhì)量阻尼器，基于等效前后顆粒阻尼器中腔體空隙體積相等以及顆粒質(zhì)量相等的原則，將多顆粒阻尼器等效為單顆粒阻尼器，建立顆粒調(diào)諧質(zhì)量阻尼器的數(shù)值模型，并進(jìn)行試驗驗證，為今后工程應(yīng)用提供一種簡單高效的分析方法。

1 數(shù)值模擬方法

本節(jié)采用等效的簡化方法來分析顆粒調(diào)諧質(zhì)量阻尼器的減震效果。顆粒調(diào)諧質(zhì)量阻尼器將顆粒阻尼器腔體通過擺繩與主體結(jié)構(gòu)相連，對于腔體可按照單擺進(jìn)行簡化，在傳統(tǒng)顆粒阻尼器的基礎(chǔ)上增加一個自由度即可。以下重點討論對多顆粒阻尼器的簡化方法。

1.1 多顆粒阻尼器等效方法

由于多顆粒在碰撞時表現(xiàn)出很強(qiáng)的非線性，因此要進(jìn)行帶有多顆粒阻尼器的動力系統(tǒng)的精確分析是相當(dāng)困難的。鑒于單顆粒阻尼器的理論分析比較完善，并考慮阻尼器的運動特點，本節(jié)采用以單顆粒阻尼器等效替代多顆粒阻尼器的方法來研究其振動響應(yīng)[3]。等效原則如下：

(1)等效前后阻尼器腔體中的空隙體積相等。

(2)等效前多顆粒總質(zhì)量與等效后單顆粒質(zhì)量相等。

(3)等效前后顆粒均為球形，且密度均為ρ。

(4)等效前的多顆粒阻尼器為長方體，等效后的單顆粒阻尼器為圓柱體，且圓柱體的直徑與等效后單顆粒的直徑相等。

依據(jù)以上原則等效后得到的單顆粒阻尼器示意圖，如圖1所示。等效前后顆粒調(diào)諧質(zhì)量阻尼器的參數(shù)及符號見表1。

圖1 多顆粒阻尼器和等效的單顆粒阻尼器Fig.1 Diagrams of particle damper and equivalent single unit impact damper

類型物理意義符號表示多顆粒阻尼器腔體尺寸(長，寬，高)(dx，dy，dz)單個顆粒質(zhì)量m1p顆?？傎|(zhì)量m顆粒直徑Dp單顆粒阻尼器顆粒質(zhì)量m3顆粒直徑D顆粒自由運動距離d

簡化方法如下：

等效前顆粒在阻尼器中占據(jù)的體積Vspd由式(1)得出，式中N是顆粒的總數(shù)，即N=m/m1p：

(1)

多顆粒阻尼器中的空隙體積Vepd由式(2)得出，式中Vpd是顆粒阻尼器的體積，并假設(shè)盒子顆粒與腔體的體積比，即顆粒填充率ρp=Vspd/Vpd：

Vepd=Vpd-Vspd=(1/ρp-1)Vspd

(2)

等效后單顆粒阻尼器的空隙體積Veid為：

(3)

根據(jù)等效原則式(1)，令式(2)和式(3)相等，可得：

(4)

HALES[9]的研究表明，半徑相等的球在發(fā)生球密堆積時體積比不能超過0.74，即ρp≤0.74。根據(jù)等效原則式(2)，m=m3，至此，建立了等效的單顆粒阻尼器中顆粒自由運動長度的表達(dá)式，即式(4)。

1.2 顆粒調(diào)諧質(zhì)量阻尼器系統(tǒng)控制方程

圖2 計算模型簡圖Fig.2 Schematic diagram of computational models

其系統(tǒng)的控制方程用矩陣表示為

[K]{X(t)}={Fp(t)}

(5)

設(shè)振動臺的位移為xg(t)，主體結(jié)構(gòu)的位移為x1(t)，質(zhì)量為m1；腔體的位移為x2(t)，質(zhì)量為m2；顆粒的位移為x3(t)，質(zhì)量為m3。則系統(tǒng)控制方程展開如

圖3 G(y)和的形式Fig.

1.3 數(shù)值模擬流程

本文對于顆粒調(diào)諧質(zhì)量阻尼器的數(shù)值模擬采用如下的流程：①確定主體結(jié)構(gòu)、顆粒調(diào)諧質(zhì)量阻尼器和荷載的相關(guān)參數(shù)，并確定結(jié)構(gòu)的初始運動狀態(tài)；②計算等效的單顆粒阻尼器的自由運動距離d；③計算主體結(jié)構(gòu)動力響應(yīng)，繪制數(shù)值模擬結(jié)果的響應(yīng)時程曲線并與實際測量曲線對比，并對比響應(yīng)的峰值結(jié)果。參數(shù)選取詳見2.2。

該數(shù)值模擬過程運用MATLAB軟件實現(xiàn)，調(diào)用程序中內(nèi)置的求常微分方程數(shù)值解的方法ode45，即四階-五階龍格-庫塔(Runge-Kutta)算法。輸入主體結(jié)構(gòu)和顆粒調(diào)諧質(zhì)量阻尼器的相關(guān)參數(shù)，以及1.2中的控制方程，即可進(jìn)行迭代求解。

2 數(shù)值模擬結(jié)果及試驗驗證

為了驗證本文提出的簡化方法用于顆粒調(diào)諧質(zhì)量阻尼器減震控制模擬的可行性和合理性，開展了一個單自由度結(jié)構(gòu)附加顆粒調(diào)諧質(zhì)量阻尼器的模擬地震振動臺試驗。

1.4 試驗簡介

試驗主體結(jié)構(gòu)模型為單自由度鋼排架，為了模擬常見的高層結(jié)構(gòu)，調(diào)整結(jié)構(gòu)的基頻至1.0 Hz左右，在結(jié)構(gòu)頂部附加6 kg質(zhì)量塊，最終其自振頻率為1.37 Hz。填充顆粒的盒子通過四根細(xì)線垂直懸吊在結(jié)構(gòu)頂部。盒子橫截面尺寸相同，均為60 mm×60 mm，長度為80 mm。盒子內(nèi)置26顆直徑為10 mm的鋼球，顆粒密度為ρ=7 644 kg/m3，單個顆粒質(zhì)量為m1p=0.032 kg。

試驗中采用El-Centro波和Kobe波作為結(jié)構(gòu)底部的振動輸入。在進(jìn)行數(shù)值模擬時，從采集的臺面加速度時程曲線中截取一個完整波作為計算輸入波。為方便起見，將各質(zhì)點的位移和速度初值賦為0。本文將計算得到的結(jié)構(gòu)頂層加速度與位移時程曲線與試驗采集得到的結(jié)果進(jìn)行比較，通過二者的對比來驗證該簡化方法的準(zhǔn)確性。關(guān)于試驗過程及結(jié)果在文獻(xiàn)[10]中有詳細(xì)介紹，本文僅介紹數(shù)值模擬參數(shù)的選取和結(jié)果。

2.2 數(shù)值模擬參數(shù)取值

主體結(jié)構(gòu)的質(zhì)量可以通過測量模型的實際質(zhì)量得到，其圓頻率和阻尼比可以通過試驗前的白噪聲掃頻，并考慮模型加工誤差，采用半功率法得到。同理，可得阻尼器腔體的質(zhì)量，考慮到調(diào)諧的因素，其圓頻率與主體結(jié)構(gòu)的自振圓頻率取值相同。

腔體的阻尼比較難測得，在數(shù)值模擬時通過試算的方式確定，計算發(fā)現(xiàn)當(dāng)腔體阻尼比低于0.05時，主體結(jié)構(gòu)的響應(yīng)幾乎不變，即主體結(jié)構(gòu)的響應(yīng)對于腔體阻尼比并不敏感。本文取腔體阻尼比為ξ2=0.05。

根據(jù)MASRI[11]的研究，當(dāng)ω3≥20ω1時，可以較合適地模擬顆粒與容器壁的相互作用。本文取顆粒圓頻率ω3=20ω1，即ω3=172.16 rad/s。

顆粒的阻尼比ξ3與顆粒的恢復(fù)系數(shù)e存在如圖4所示的關(guān)系[5]。在本文的試驗中，根據(jù)文獻(xiàn)[12]，顆粒的恢復(fù)系數(shù)e一般取為0.25，由圖4關(guān)系可得ξ3=0.375。

數(shù)值模擬中用到的系統(tǒng)參數(shù)匯,總?cè)绫?所示。

由于振動臺面的數(shù)據(jù)采集頻率為500 Hz，因此上述數(shù)值模擬過程均以0.002 s為計算步長。

圖4 顆粒阻尼比與恢復(fù)系數(shù)關(guān)系圖Fig.4 Relationship between damping ratio ξ3 and coefficient of restitution e

Tab.2 Values of system parameters

質(zhì)量/kg圓頻率/(rad·s-1)阻尼比主體結(jié)構(gòu)7．68．6080．032阻尼器腔體0．1288．6080．05等效后的單顆粒0．832172．160．375

2.3 數(shù)值模擬結(jié)果

圖5顯示了附加顆粒調(diào)諧質(zhì)量阻尼器的單層鋼排架模型頂層在El-Centro波地震激勵下加速度和位移響

應(yīng)的計算值和試驗值對比曲線?？傮w上看，兩者符合較好，但有部分區(qū)段計算值和試驗值存在一定誤差。一方面是由于為方便起見，數(shù)值模擬時將位移和速度初值賦為零所致，隨著計算的逐步進(jìn)行，初值的選取對結(jié)構(gòu)在地震激勵下響應(yīng)的影響逐漸減小；另一方面是由于顆粒阻尼器在地震激勵下表現(xiàn)出很強(qiáng)的非線性，本文提出的數(shù)值模擬方法對此有一定程度的簡化，比如并未考慮顆粒之間的碰撞的影響。

此外，表3列出了附加不同顆粒形式的顆粒調(diào)諧質(zhì)量阻尼器的結(jié)構(gòu)模型在不同激勵下的頂層峰值和均方根加速度響應(yīng)計算值與試驗值的對比，發(fā)現(xiàn)峰值加速度兩者吻合較好，而均方根加速度也可以將誤差控制在可接受的范圍內(nèi)。事實上，在高層建筑結(jié)構(gòu)抗震設(shè)計的時候，對于主體結(jié)構(gòu)舒適度的要求，主要就反映在最大加速度的數(shù)值上面，設(shè)計人員對該數(shù)值十分關(guān)心。而本簡化方法能對該數(shù)值給出一定程度的合理預(yù)測。

圖5 附加顆粒調(diào)諧質(zhì)量阻尼器模型頂層在地震激勵(El-Centro波)下的響應(yīng)Fig.5 Response at roof level of the test frame with particle tuned mass damper under earthquake excitations (El-Centro Wave)

地震波類型顆粒形式峰值加速度計算值/(m·s-2)試驗值/(m·s-2)誤差/%均方根加速度計算值/(m·s-2)試驗值/(m·s-2)誤差/%4個20mm鋼球2．60242．56451．460．56960．53785．58El?Centro5個20mm鋼球2．60492．6949-3．460．55820．55380．7919個10mm鋼球1．96741．92991．900．57710．54944．8126個10mm鋼球1．88331．82593．050．57630．56072．714個20mm鋼球3．52153．5518-0．860．50560．5377-6．34Kobe5個20mm鋼球3．74003．73540．120．53450．5799-8．4919個10mm鋼球3．88693．85620．790．52450．5751-9．6426個10mm鋼球3．64103．7484-2．950．49830．5435-9．07

總體而言，本文提出的數(shù)值模擬方法能夠較合理地計算出附加顆粒調(diào)諧質(zhì)量阻尼器系統(tǒng)在實際地震激勵下的響應(yīng)，并且由于其原理簡單、操作方便，具有一定的實用價值。然而，要更精確的模擬附加顆粒調(diào)諧質(zhì)量阻尼系統(tǒng)在地震激勵下的行為，尤其是考慮顆粒之間碰撞對減振性能的影響，還需要進(jìn)一步的研究。

3 結(jié) 論

本文在顆粒阻尼器和調(diào)諧質(zhì)量阻尼器的基礎(chǔ)上，結(jié)合各自的優(yōu)點，提出一種新型的顆粒調(diào)諧質(zhì)量阻尼器。基于等效前后顆粒阻尼器內(nèi)空隙與腔體體積比一致以及多顆?？傎|(zhì)量與等效后單顆粒質(zhì)量相等的原則，將多顆粒阻尼器簡化為單顆粒阻尼器，建立顆粒調(diào)諧質(zhì)量阻尼器的數(shù)值模型，并介紹了系統(tǒng)參數(shù)的取值方法。為了驗證該方法的可行性和準(zhǔn)確性，進(jìn)行了附加顆粒調(diào)諧質(zhì)量阻尼器的單自由度結(jié)構(gòu)實際地震輸入下的振動臺試驗，通過數(shù)值模擬結(jié)果與試驗結(jié)果的對比發(fā)現(xiàn)兩者整體吻合較好，并且加速度峰值和均方根也吻合良好，說明該簡化方法可以較好地模擬顆粒調(diào)諧質(zhì)量阻尼器的減震控制效果，有助于進(jìn)一步了解其物理本質(zhì)，從而為進(jìn)一步推廣其在土木工程中的應(yīng)用提供有效的分析途徑。

[ 1 ] 李曉瑋,施衛(wèi)星. 人行天橋MTMD減振控制的魯棒性研究[J]. 結(jié)構(gòu)工程師,2012，28(4)：7-12. LI Xiaowei，SHI Weixing. Research on robustness of vibration control of pedestrian bridges using multi degree of freedom passive tuned mass-dampers [J].Structure Engineers, 2012，28(4)：7-12.

[ 2 ] 魯正,呂西林,閆維明. 顆粒阻尼技術(shù)研究綜述[J]. 振動與沖擊,2013,32(7):1-7. LU Zheng，Lü Xilin，YAN Weiming. A survey of particle damper technology [J]. Journal of Vibration and Shock, 2013, 32(7):1-7.

[ 3 ] PAPALOU A，MASRI S F. Performance of particle dampers under random excitation [J].Journal of Vibration and Acoustics,1996，118(4):614-621.

[ 4 ] FRIEND R D，KINRA V K. Particle impacting damping [J]. Journal of Sound and Vibration,2000,233(1):93-118.

[ 5 ] 魯正,呂西林. 顆粒阻尼器減震控制的數(shù)值模擬[J]. 同濟(jì)大學(xué)學(xué)報(自然科學(xué)版),2013, 41(8):1140-1144. LU Zheng，Lü Xilin. Numerical simulation of vibration control effects of particle dampers [J]. Journal of Tongji University (Natural Science),2013,41(8):1140-1144.

[ 6 ] 閆維明,許維炳,王瑾,等. 調(diào)諧型顆粒阻尼器簡化力學(xué)模型及其參數(shù)計算方法研究與減震橋梁試驗[J]. 工程力學(xué),2014,31(6):79-84. YAN Weiming，XU Weibing，WANG Jin，et al. Experimental and theoretical research on the simplified mechanical model of a tuned particle damper，its parameter determination method and earthquake-induced vibration control of bridge [J]. Engineering Mechanics,2014, 31(6):79-84.

[ 7 ] 胡溧,黃其柏,許智生. 顆粒阻尼的回歸分析研究[J]. 中國機(jī)械工程,2008,19(23):2834-2837. HU Li,HUANG Qibai,XU Zhisheng. Regression analysis of particle damping [J]. China Mechanical Engineering,2008, 19(23):2834-2837.

[ 8 ] 蔣華,陳前. 恢復(fù)力曲面法在顆粒阻尼器研究中的應(yīng)用[J]. 振動、測試與診斷,2007,27(3):228-231. JIANG Hua,CHEN Qian. Application of restoring force surface method in particle damping research [J]. Journal of Vibration, Measurement & Diagnosis, 2007, 27(3):228-231.

[ 9 ] HALES T C. The sphere packing problem [J]. Journal of Computational and Applied Mathematics,1992，164(1): 53-65.

[10] 施衛(wèi)星,何斌,李曉瑋,等. 一種新型調(diào)諧質(zhì)量阻尼器的試驗研究[J]. 振動與沖擊, 2015,34(12):207-211. SHI Weixing，HE Bin，LI Xiaowei，et a. Experimental study on a new type of tuned mass damper [J]. Journal of Vibration and Shock, 2015, 34(12):207-211.

[11] MASRI S F,IBRAHIM A M. Response of the impact damper to stationary random excitation [J]. The Journal of the Acoustical Society of America, 1973, 53(1):200-211.

[12] 葛藤,賈智宏,周克棟. 鋼球和剛性平面彈塑性正碰撞恢復(fù)系數(shù)研究[J]. 工程力學(xué),2008,25(6):209-213. GE Teng,JIA Zhihong,ZHOU Kedong. Research on elastoplastic normal impact of steel spheres against a rigid plane [J]. Engineering Mechanics, 2008, 25(6):209-213.

Numerical simulation for vibration reduction control of particle tuned mass dampers

LU Zheng1,2, CHEN Xiaoyi2, WANG Dianchao2, Lü Xilin1,2

(1. State Key Laboratory of Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China;2.Research Institute of Structural Engineering and Disaster Reduction, Tongji University, Shanghai 200092, China)

A particle tuned mass damper was proposed by hanging a particle damper on a main structure through a string. It combines a widely used tuned mass damper with an efficient particle damper in vibration control area. Thus it is able to increase the robustness of damping performance within a wider frequency band. Here, an approximate numerical simulation method for particle tuned mass dampers was proposed based on an equivalent priciple to take a multi-particle damper as an equivalents single-particle damper. The method to determine the corresponding system parameters was also introduced. A series of shaking table tests of a single DOF structure attached with a particle tuned mass damper under seismic imputs were performed. A comparative study showed that the simulation results agree well with the test ones, so the proposed simplified method can be used to simulate vibration reduction control effects of particle tuned mass dampers.Key words: particle tuned mass damper; tuned mass damper; particle damper; numerical simulation; structural vibration control

國家自然科學(xué)基金(51478361);中央高?；究蒲袠I(yè)務(wù)費專項資金項目

2015-11-20 修改稿收到日期：2016-01-24

魯正 男，博士，副教授，1982年生

TU352；TB535

A

10.13465/j.cnki.jvs.2017.03.008