曼型干式煤氣柜動(dòng)力特性分析及風(fēng)振響應(yīng)主要貢獻(xiàn)模態(tài)識(shí)別

李正良1，2， 劉欣鵬1， 晏致濤1，2， 焦紅偉1， 俞登科1

(1.重慶大學(xué)土木工程學(xué)院，重慶400045；2.重慶大學(xué)山地城鎮(zhèn)建設(shè)與新技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，重慶400045)

摘要：曼型干式煤氣柜不同于一般空間結(jié)構(gòu)，其動(dòng)力特性十分復(fù)雜，作為研究其頻域內(nèi)風(fēng)振響應(yīng)的基礎(chǔ)，首先通過(guò)建立精細(xì)化有限元模型研究其在不同工況下的動(dòng)力特性；其次，分別計(jì)算背景、共振響應(yīng)振型能量參與系數(shù)，并按降序排列，選取結(jié)構(gòu)風(fēng)振響應(yīng)的主要模態(tài);最后，通過(guò)選取模態(tài)進(jìn)行頻域計(jì)算并與傳統(tǒng)頻域計(jì)算方法(完全二次振型疊加，CQC法)對(duì)比。結(jié)果表明：結(jié)構(gòu)自振頻率分布密集；動(dòng)力特性受活塞位置及內(nèi)部氣壓影響較大；背景、共振響應(yīng)振型能量參與系數(shù)能準(zhǔn)確反映各階模態(tài)能量對(duì)響應(yīng)總能量的貢獻(xiàn)程度；通過(guò)振型能量參與系數(shù)法確定的主要模態(tài)進(jìn)行頻域分析具有計(jì)算效率高，計(jì)算精度可控的優(yōu)點(diǎn)。

關(guān)鍵詞：動(dòng)力特性；背景響應(yīng)；共振響應(yīng)；脈動(dòng)風(fēng)荷載；振型能量參與系數(shù)；頻域分析

中圖分類號(hào):TU973+.213；TU33+3

文獻(xiàn)標(biāo)志碼：A

DOI:10.13465/j.cnki.jvs.2015.19.012

Abstract:The dynamic characteristics of a MAN type dry gas storage tank are very complex, they are different from those of ordinary spatial structures. Here, as the basis of its wind-induced response analysis in frequency domain, its dynamic characteristics were analyzed under different working conditions by using the finite element method. Its modes were sorted with order-descending according to the calculation results of mode energy participation factors of background response and resonant response, respectively. Then, the modes with significant contributions to its wind-induced vibration were selected. Finally, the results calculated in frequency domain with a selected modes were compared with those calculated with CQC method. The results showed that the modal frequency distribution of the structure is dense; the location of piston and internal pressure have larger influences on the dynamic characteristics of the gas storage tank; the mode energy participation factors of background response and resonant one can reflect the contributions of different order modes’ energy to the total energy of the response; the analysis in frequency domain with the main modes determined using the mode energy participation factors has advantages of higher calculation efficiency and controllable computation accuracy.

基金項(xiàng)目：國(guó)家自然科學(xué)基金(61004067)資助項(xiàng)目

收稿日期：2014-08-01修改稿收到日期：2014-09-18

Dynamic characteristics of a MAN type dry gas storage tank and recognition of main modes with contributions to its wind-induced vibration

LIZheng-liang1,2,LIUXin-peng1,YANZhi-tao1,2,JIAOHong-wei1,YUDeng-ke1(1. School of Civil Engineering, Chongqing University, Chongqing 400045, China；2. Key Laboratory of New Technology for Construction of Cities in Mountain Area,Ministry of Education, Chongqing University, Chongqing 400045, China)

Key words:dynamic characteristics; background response; resonant response; fluctuating wind loads; mode energy participation factor; analysis in frequency domain

曼型干式煤氣柜(下文簡(jiǎn)稱新型煤氣柜)是一種儲(chǔ)存可燃?xì)怏w的大型重要建筑物，廣泛用于冶金、石化和市政等行業(yè)，其在節(jié)約能源和保護(hù)環(huán)境上發(fā)揮著重大作用。近年來(lái)，隨著我國(guó)鋼鐵行業(yè)的發(fā)展，國(guó)內(nèi)開(kāi)始大量建造該類煤氣柜。

新型煤氣柜作為空間結(jié)構(gòu)的代表，與傳統(tǒng)高層建筑結(jié)構(gòu)多側(cè)重抗震性能的研究不同[1]，其風(fēng)振響應(yīng)已成為控制結(jié)構(gòu)設(shè)計(jì)的關(guān)鍵因素之一。因此對(duì)此類結(jié)構(gòu)風(fēng)振響應(yīng)分析成為了亟待解決的重要問(wèn)題。從目前的研究狀況來(lái)看，在模態(tài)分解法的基礎(chǔ)上，利用隨機(jī)振動(dòng)理論來(lái)求解結(jié)構(gòu)在脈動(dòng)風(fēng)作用下的動(dòng)力風(fēng)致響應(yīng)即頻域法，是分析空間結(jié)構(gòu)的首選方法。

然而，與一般高層建筑結(jié)構(gòu)形式不同，新型煤氣柜具有高跨比較??；柜壁較薄，易產(chǎn)生局部振動(dòng)；內(nèi)部活塞質(zhì)量巨大，對(duì)煤氣柜動(dòng)力特性影響顯著等特點(diǎn)。國(guó)內(nèi)對(duì)于此類新型結(jié)構(gòu)形式的動(dòng)力特性分析尚屬空白。因此作為頻域分析法的基礎(chǔ)，對(duì)其動(dòng)力特性的精細(xì)化分析尤為重要。

另外，在頻域內(nèi)計(jì)算大跨度空間結(jié)構(gòu)的脈動(dòng)風(fēng)振響應(yīng)一般采用兩種方法：基于荷載空間分布模式的Ritz向量疊加法[2]和完全二次振型疊加法(CQC)[3-4]。Ritz向量疊加法基本特點(diǎn)是在確定結(jié)構(gòu)主要模態(tài)時(shí)，考慮了外荷載的空間分布模式對(duì)激發(fā)結(jié)構(gòu)模態(tài)特性的影響，過(guò)濾掉與荷載分布向量正交的振型，從而在進(jìn)行頻域計(jì)算是大大減少了計(jì)算量，提高了計(jì)算效率，然而對(duì)于時(shí)空分布較復(fù)雜的脈動(dòng)風(fēng)荷載來(lái)說(shuō)，確定荷載分布模式較為困難，限制了Ritz向量疊加法在風(fēng)振響應(yīng)分析中的應(yīng)用。對(duì)于傳統(tǒng)的頻域計(jì)算方法，即完全二次振型疊加法，該方法自動(dòng)包含了所有主要模態(tài)的耦合項(xiàng)，計(jì)算精度較高。工程分析中，往往取前10～20階模態(tài)進(jìn)行計(jì)算，但這種做法，對(duì)于大型空間結(jié)構(gòu)來(lái)說(shuō)，往往會(huì)忽略某些對(duì)風(fēng)振響應(yīng)貢獻(xiàn)較大的高階模態(tài)[5-10]，因此如何合理考慮高階振型對(duì)結(jié)構(gòu)風(fēng)振響應(yīng)的影響是采用頻域分析法計(jì)算結(jié)構(gòu)風(fēng)振響應(yīng)的重要問(wèn)題。對(duì)于曼型煤氣柜這種新型結(jié)構(gòu)形式，國(guó)內(nèi)對(duì)其動(dòng)力特性及主要貢獻(xiàn)模態(tài)的研究尚屬空白。本文總結(jié)了振型疊加法計(jì)算結(jié)構(gòu)風(fēng)振響應(yīng)的公式[11]，以及背景響應(yīng)與共振響應(yīng)的計(jì)算方法[12-14]，并將背景—共振響應(yīng)振型能量參與系數(shù)法引入該類復(fù)雜空間結(jié)構(gòu)中，從能量的角度定量的選取風(fēng)振響應(yīng)中的主要貢獻(xiàn)模態(tài)，解決了該類結(jié)構(gòu)風(fēng)振響應(yīng)計(jì)算精度的評(píng)價(jià)問(wèn)題。

1煤氣柜動(dòng)力特性分析

1.1煤氣柜簡(jiǎn)介

圖1　新型煤氣柜幾何模型圖 Fig.1 The geometry model of new type dry gas tank

圖2　活塞幾何模型圖 Fig.2 The geometry model of new type dry gas tank piston

曼型干式煤氣柜由外筒、柜頂、內(nèi)部活塞以及附屬結(jié)構(gòu)(通風(fēng)塔樓、外部電梯、外部樓梯等)組成，本文研究的煤氣柜的儲(chǔ)氣規(guī)模為30萬(wàn)m3，其中該類煤氣柜的主要外形尺寸為：煤氣柜總高120m，柜體側(cè)壁高108m，外筒半徑66m，活塞總行程96m，柜頂與活塞的拱頂曲率半徑為88m，所有型鋼及鋼板為Q235鋼，見(jiàn)圖1～圖2。

1.2無(wú)內(nèi)壓作用煤氣柜模態(tài)分析

不考慮柜體內(nèi)外氣體壓力差，利用ANSYS有限元分析軟件分別對(duì)煤氣柜在三種工況(活塞高位、中位、低位)下的動(dòng)力特性進(jìn)行分析。由于煤氣柜是軸對(duì)稱結(jié)構(gòu)，因此在模態(tài)分析時(shí)均將振動(dòng)形式相同的模態(tài)合并。

1.2.1柜體在活塞高位狀態(tài)下的模態(tài)和頻率

煤氣柜在活塞位于高位時(shí)(活塞底部距離地面97.2m)，通過(guò)模態(tài)分析得到結(jié)構(gòu)的頻率和模態(tài)如表1和圖3(由于篇幅限制僅給出典型模態(tài)，下同)。

表1　 新型煤氣柜結(jié)構(gòu)的自振頻率

圖3　煤氣柜振型圖 Fig.3 Vibration chart of new type gas tank

分析可知：當(dāng)活塞位于高位時(shí)，其自振頻率的分布比較密集，前40階模態(tài)頻率分布在1.3～7Hz；柜頂?shù)恼駝?dòng)模態(tài)僅出現(xiàn)在第14和25階，這一定程度說(shuō)明柜頂?shù)膭偠纫黠@大于柜體；結(jié)構(gòu)前25階振型中，一、二階振型為x,y方向的懸臂梁式振動(dòng)，第三階振型為結(jié)構(gòu)的扭轉(zhuǎn)振動(dòng)，其余各階均為殼體的局部振動(dòng)。

1.2.2柜體在活塞中位狀態(tài)下的模態(tài)和頻率

煤氣柜在活塞位于中位時(shí)(活塞底部距離地面54m)，通過(guò)模態(tài)分析得到結(jié)構(gòu)的頻率和振型如表2和圖4。

表2　新型煤氣柜結(jié)構(gòu)的自振頻率

圖4　煤氣柜振型圖 Fig.4 Vibration chart of new type gas tank

分析可知：當(dāng)活塞位于中位時(shí)，煤氣柜模態(tài)頻率比高位模態(tài)頻率大，柜體整體剛度增大；與高位時(shí)模態(tài)一樣，自振頻率的分布比較密集，前40階振型的頻率分布在2～7Hz；結(jié)構(gòu)前25階振型中，一、二階振型為x,y方向的懸臂梁式振動(dòng)，第三階振型為結(jié)構(gòu)的扭轉(zhuǎn)振動(dòng)，這與活塞高位時(shí)的振動(dòng)形態(tài)相似，但其對(duì)應(yīng)的自振頻率與活塞高位差別較大。對(duì)于高階模態(tài)，盡管兩種工況下(活塞高、中位)均為殼體的局部振動(dòng)，但相應(yīng)階數(shù)的模態(tài)振動(dòng)形式及其對(duì)應(yīng)自振頻率均變化較大，這是因?yàn)榛钊奈恢米兓瘜?duì)結(jié)構(gòu)的動(dòng)力特性影響很大。

1.2.3柜體在活塞低位狀態(tài)下的模態(tài)和頻率

煤氣柜在活塞低位狀態(tài)下(活塞底部距離地面3.6m)，通過(guò)模態(tài)分析得到結(jié)構(gòu)的頻率和振型如表3和圖5。

表3　新型煤氣柜結(jié)構(gòu)的自振頻率

圖5　煤氣柜振型圖 Fig.5 Vibration chart of new type gas tank

分析可知：結(jié)構(gòu)模態(tài)頻率比前兩種工況相應(yīng)模態(tài)頻率大；前40階振型的頻率分布在2.7～7Hz，自振頻率的分布相對(duì)高位和中位更密集；隨著活塞位置的繼續(xù)下降，煤氣柜前25階振型中，不同于高位和中位振型，低位的一階振型為環(huán)向四個(gè)波形的殼體振動(dòng)，第二、三階振型為懸臂梁式振動(dòng)，柜體的扭轉(zhuǎn)振動(dòng)出現(xiàn)在第18階，其余為殼體的局部模態(tài)。

1.3內(nèi)壓作用下煤氣柜結(jié)構(gòu)模態(tài)分析

煤氣柜處于工作狀態(tài)時(shí)，內(nèi)部封閉有一定壓力的煤氣，30萬(wàn)m3的煤氣柜設(shè)計(jì)內(nèi)壓為15kPa。內(nèi)壓的存在將會(huì)影響結(jié)構(gòu)的剛度，并且對(duì)結(jié)構(gòu)的動(dòng)力特性產(chǎn)生影響。故本文主要針對(duì)活塞處于高位狀態(tài)時(shí)，內(nèi)壓作用與否對(duì)煤氣柜模態(tài)分析。

表4　新型煤氣柜結(jié)構(gòu)的自振頻率

表4為活塞在高位狀態(tài)下，考慮內(nèi)壓與未考慮內(nèi)壓煤氣柜結(jié)構(gòu)的頻率與振型特點(diǎn)描述(由于振型相似及篇幅原因內(nèi)壓作用下的振型截圖不再給出)。從表4可以看出，考慮內(nèi)壓對(duì)煤氣柜的平動(dòng)振型和扭轉(zhuǎn)振型幾乎沒(méi)有影響，有內(nèi)壓和無(wú)內(nèi)壓的平動(dòng)頻率和扭轉(zhuǎn)頻率均相同；而內(nèi)壓對(duì)柜體殼體局部振動(dòng)影響較顯著，考慮內(nèi)壓的影響，殼體振動(dòng)的頻率要比不考慮內(nèi)壓最大增加12%左右。由此可知，柜體內(nèi)部氣壓對(duì)其局部剛度影響顯著，而對(duì)整體剛度幾乎沒(méi)有影響。

2背景-共振響應(yīng)振型能量參與系數(shù)

2.1脈動(dòng)風(fēng)振響應(yīng)振型疊加原理

對(duì)于新型煤氣柜結(jié)構(gòu)在外施動(dòng)力荷載作用下的動(dòng)力方程為：

(1)

結(jié)構(gòu)振型矩陣Φ=[φ1,φ2,…,φn]，且ΦTM Φ=I(其中I為單位矩陣)，則ΦTK Φ=Ω2(其中Ω為結(jié)構(gòu)前n階圓頻率組成的對(duì)角矩陣，即diag[ω1,ω2,…,ωn]結(jié)構(gòu)在脈動(dòng)風(fēng)作用下的位移響應(yīng)方差矩陣為[15]：

LTΦH*T(ω)ΦTdω

(2)

式中，SPP(ω)測(cè)壓點(diǎn)脈動(dòng)風(fēng)壓時(shí)程互譜矩陣；H(ω)為頻響函數(shù)矩陣，且H(ω)=[H1(ω),H2(ω),…,Hn(ω)]，其中，Hj(ω)為第j階模態(tài)的頻響函數(shù)。

2.2背景響應(yīng)及其能量參與系數(shù)的計(jì)算

將脈動(dòng)風(fēng)作用下的結(jié)構(gòu)振動(dòng)看作靜力過(guò)程，即不考慮結(jié)構(gòu)的動(dòng)力效應(yīng)，可得到結(jié)構(gòu)總的背景響應(yīng)及各振型的背景響應(yīng)。

(1)總背景響應(yīng)

脈動(dòng)風(fēng)作用下，總的背景響應(yīng)Xb可表示為：

KXb=LPd(t)

(3)

Xb的方差可表示為：

(4)

式中，XP為靜力平衡方程KXb=L的解；b為背景響應(yīng)。

(2)振型背景響應(yīng)

設(shè)Xb=ΦQb，將其代入式(3)，并對(duì)等式左右同時(shí)左乘ΦT，那么式(3)可改寫(xiě)為：

Ω2Qb=ΦTLPd(t)

(5)

得到：

Qb=H0ΦTLPd(t)

(6)

則第j階模態(tài)背景響應(yīng)為：

qb,j=Hj,0φTjLPd(t)

(7)

第j階模態(tài)對(duì)總背景響應(yīng)的貢獻(xiàn)為：

Xb,j=φjHj,0φTjLPd(t)

(8)

振型位移Qb的方差為：

(9)

式中，第j階模態(tài)背景響應(yīng)的方差為：

(10)

(3)背景響應(yīng)振型能量參與系數(shù)

脈動(dòng)風(fēng)荷載在背景響應(yīng)上做的功為：

(11)

式中，∑diag[]表示矩陣的對(duì)角元素之和。

對(duì)背景響應(yīng)功Wb取數(shù)學(xué)期望，得到

(12)

脈動(dòng)風(fēng)荷載在第j階模態(tài)背景響應(yīng)上做的功為：

(13)

對(duì)在第j階模態(tài)背景響應(yīng)功Wb,j取數(shù)學(xué)期望，得到：

(14)

第j階振型背景響應(yīng)的能量參與系數(shù)是：

(15)

由式(15)可知，背景響應(yīng)的振型能量參與系數(shù)反映了振型位移對(duì)總背景響應(yīng)位移的平均貢獻(xiàn)率。背景響應(yīng)位移計(jì)算值的精度可通過(guò)振型累積參與系數(shù)以及振型參與系數(shù)來(lái)控制。

2.4共振響應(yīng)及其能量參與系數(shù)的計(jì)算

共振響應(yīng)是指脈動(dòng)風(fēng)荷載在共振頻率處激起的結(jié)構(gòu)振動(dòng)，對(duì)第j階振型的運(yùn)動(dòng)方程兩側(cè)進(jìn)行傅里葉變換，得到第j階振型響應(yīng)的頻域解為：

qj(ω)=Hj(ω)φTjLPd(ω)

(16)

第j階振型的共振響應(yīng)為：

qr,j(ω)=Hj(ω)φTjLPd(ωj)

(17)

式中，下標(biāo)“r”表示共振響應(yīng)；Pd(ωj)為脈動(dòng)風(fēng)壓傅里葉譜值在第j階振型頻率ωj處的取值。

第j階振型共振響應(yīng)的位移方差為：

(18)

式中，下標(biāo)“d”表示位移響應(yīng)；ξj為第j階模態(tài)阻尼比；SPP(ωj)為脈動(dòng)風(fēng)壓功率譜在第j階模態(tài)頻率ωj處的取值。

第j階振型共振響應(yīng)的速度方差為：

(19)

式中，下標(biāo)“v”表示位移響應(yīng)。

脈動(dòng)風(fēng)荷載在第j階振型共振響應(yīng)上的功等于第j階振型的動(dòng)能和彈性能之和：

(20)

對(duì)Wr,j取數(shù)學(xué)期望，得到：

(21)

共振響應(yīng)總能量，從理論上來(lái)說(shuō)等于全部模態(tài)共振響應(yīng)的能量之和，但當(dāng)結(jié)構(gòu)自由度十分龐大時(shí)，共振總能量很難得到。因此，為了方便起見(jiàn)，將共振能量參與系數(shù)定義為各階共振響應(yīng)能量與背景響應(yīng)總能量之比，如下式：

(22)

3煤氣柜主要貢獻(xiàn)模態(tài)選取

新型煤氣柜為軸對(duì)稱結(jié)構(gòu)體系，因此對(duì)其進(jìn)行風(fēng)振響應(yīng)頻域分析時(shí)僅考慮0°風(fēng)向角情況，本文運(yùn)用線性濾波AR模型法[16]對(duì)脈動(dòng)風(fēng)速進(jìn)行模擬，模擬結(jié)果見(jiàn)圖6～圖7。

圖6　脈動(dòng)風(fēng)時(shí)程曲線 Fig.6 Turbulence wind velocity time-history curve

圖7　模擬風(fēng)譜與目標(biāo)譜的比較 Fig.7 The target spectrum and simulation spectrum

3.1煤氣柜背景響應(yīng)主要貢獻(xiàn)模態(tài)

由于煤氣柜內(nèi)部活塞上下移動(dòng)而造成結(jié)構(gòu)體系動(dòng)力特性的時(shí)變性，因此本文首先給出活塞三個(gè)典型位置即高位、中位、低位前100階背景響應(yīng)能量參與系數(shù)，見(jiàn)圖8。前文分析可知，煤氣柜內(nèi)部氣壓對(duì)結(jié)構(gòu)模態(tài)影響較大，因此以下工況分析均與實(shí)際情況相同即考慮內(nèi)部氣壓對(duì)結(jié)構(gòu)的影響。

接下來(lái)分別對(duì)三種工況的背景響應(yīng)能量參與系數(shù)按照降序排列，并進(jìn)行累加計(jì)算。當(dāng)背景響應(yīng)振型能量參與系數(shù)<2%，并且累積系數(shù)使之>90%時(shí)，對(duì)背景響應(yīng)主要貢獻(xiàn)模態(tài)選取結(jié)束，篩選結(jié)果見(jiàn)圖9。

圖8　背景響應(yīng)能量參與系數(shù) Fig.8 The mode energy participation factorof background response

圖9　背景響應(yīng)主要貢獻(xiàn)模態(tài) Fig.9 The natural modes with significant contribution to background response

3.2煤氣柜共振響應(yīng)主要貢獻(xiàn)模態(tài)

分別給出三種工況下的共振能量參與系數(shù)，見(jiàn)圖10。與背景響應(yīng)主要模態(tài)選取原則不同，對(duì)共振響應(yīng)主要貢獻(xiàn)模態(tài)選取時(shí)，當(dāng)共振能量參與系數(shù)<2%時(shí)即可結(jié)束篩選，這主要是因?yàn)樵谟?jì)算共振能量參與系數(shù)時(shí)，所用歸一化系數(shù)為背景響應(yīng)總能量值，這將造成共振累積能量很難達(dá)到90%以上，篩選結(jié)果見(jiàn)圖11。

需要注意的是，本文風(fēng)荷載的采樣頻率為16Hz，所以在對(duì)共振主要模態(tài)的選取時(shí)，僅能考慮頻率低于16Hz的模態(tài)，對(duì)于頻率大于16Hz的模態(tài)，隨著模態(tài)階數(shù)的增加，其頻率越來(lái)越大，相對(duì)低頻的荷載，越來(lái)越難以激起高階模態(tài)的共振，因此高階模態(tài)共振響應(yīng)越來(lái)越小，可以忽略。

圖10　共振響應(yīng)能量參與系數(shù) Fig.10 The mode energy participation factor of resonant response

圖11　共振響應(yīng)主要貢獻(xiàn)模態(tài) Fig.11 The natural modes with significant contribution to resonant response

對(duì)于共振響應(yīng)能量，由圖11分析可知，篩選出的主要模態(tài)對(duì)共振響應(yīng)總能量的貢獻(xiàn)都不大，這主要是兩方面造成的：結(jié)構(gòu)模態(tài)分布密集，脈動(dòng)風(fēng)壓能夠激起多階模態(tài)的振動(dòng)，這可以理解為多階模態(tài)均攤了共振響應(yīng)總能量；由于結(jié)構(gòu)局部模態(tài)剛度較大，各局部模態(tài)自振頻率與脈動(dòng)風(fēng)荷載頻率相差較遠(yuǎn)，使得脈動(dòng)風(fēng)荷載無(wú)法激起某階模態(tài)的較大共振。

為了驗(yàn)證上述模態(tài)選取方法的準(zhǔn)確性，本文首先采用完全二次振型疊加法(CQC法)按式(2)分別計(jì)算三種工況下煤氣柜在脈動(dòng)風(fēng)作用下的總響應(yīng)，并將前100階模態(tài)計(jì)算得到的位移均方根作為準(zhǔn)確值。接著按上述方法識(shí)別煤氣柜主要貢獻(xiàn)模態(tài)，并通過(guò)其進(jìn)行頻域分析。最后，給出按上述兩種方法計(jì)算得到的結(jié)構(gòu)典型位置測(cè)點(diǎn)(高度100.8m)的脈動(dòng)響應(yīng)均方根值，見(jiàn)圖12。

圖12　CQC法與振型能量參與系數(shù)法的結(jié)果 Fig.12 The results of CQC method and mode energy participation factor method

圖13　振型能量參與系數(shù)法的計(jì)算效率 Fig.13 The efficiency ofmode energy participation factor method

結(jié)果表明，三種工況下，分別通過(guò)主要模態(tài)與前100階模態(tài)進(jìn)行頻域分析的結(jié)果吻合很好，證明主要貢獻(xiàn)模態(tài)選取合理，準(zhǔn)確。

為了驗(yàn)證通過(guò)主要貢獻(xiàn)模態(tài)進(jìn)行結(jié)構(gòu)頻域分析的計(jì)算效率，本文給出了煤氣柜6個(gè)典型測(cè)點(diǎn)在三種工況下參與計(jì)算的主要貢獻(xiàn)模態(tài)數(shù)不同，其分析結(jié)果與標(biāo)準(zhǔn)值的比值，見(jiàn)圖13。

結(jié)果表明，隨著主要貢獻(xiàn)模態(tài)數(shù)的不斷累加，分析結(jié)果迅速收斂，這比傳統(tǒng)的頻域分析方法效率更高，且計(jì)算精度可以通過(guò)調(diào)整篩選原則控制。對(duì)于這種自振頻率密集，參振模態(tài)多的結(jié)構(gòu)體系，振型能量參與系數(shù)能夠很好的識(shí)別主要模態(tài)。

4結(jié)論

(1)活塞位置對(duì)新型煤氣柜動(dòng)力特性的影響顯著，隨著活塞位置降低，煤氣柜自振頻率逐漸變大。

(2)煤氣柜在活塞高位和中位狀態(tài)下，一、二階振型為x,y方向的懸臂梁式振動(dòng)，第三階振型為結(jié)構(gòu)的扭轉(zhuǎn)振動(dòng)；不同于高位和中位振型，低位的一階振型為環(huán)向四個(gè)波形的殼體振動(dòng)，第二、三階振型為懸臂梁式振動(dòng)，柜體的扭轉(zhuǎn)振動(dòng)出現(xiàn)在第18階，其余各階局部振動(dòng)。

(3)內(nèi)壓對(duì)煤氣柜的平動(dòng)振型和扭轉(zhuǎn)振型幾乎沒(méi)有影響，而內(nèi)壓對(duì)柜體殼體振動(dòng)影響較顯著，考慮內(nèi)壓的影響，殼體振動(dòng)的頻率要比不考慮內(nèi)壓最大增加12%左右。

(4)脈動(dòng)風(fēng)振響應(yīng)的計(jì)算精度與模態(tài)選取密切相關(guān)，模態(tài)的選取取決于各階模態(tài)響應(yīng)對(duì)總響應(yīng)的貢獻(xiàn)大小。背景能量參與系數(shù)法能夠準(zhǔn)確識(shí)別出煤氣柜脈動(dòng)風(fēng)振響應(yīng)的主要貢獻(xiàn)模態(tài)。

(5)通過(guò)調(diào)整篩選原則，可以對(duì)煤氣柜頻域分析計(jì)算精度進(jìn)行控制，進(jìn)而能有效提高頻域分析效率。

參考文獻(xiàn)

[1]沈世釗.大跨空間結(jié)構(gòu)的發(fā)展—回顧與展望[J].土木工程學(xué)報(bào),1998,31(3):5-14.

SHEN Shi-zhao. Development Of Long-Span Structures—A Review And Prospect[J]. China Civil Engineering Journal,1998,31(3):5-14.

[2]Wilson E L, YUAN Ming-wu. Dynamic analysis by direct superposition of Ritz vectors [J]. Earthquake Engineering and Structural Dynamics, 1982(10):813-821.

[3]顧明,黃鵬,周毅,等.北京首都機(jī)場(chǎng)3號(hào)航站樓風(fēng)荷載和響應(yīng)研究[J].土木工程學(xué)報(bào),2005,38(1): 40-44.

GU Ming, HUANG Peng, ZHOU Yi,et al. A study on wind loads and responses of terminal 3 at beijing capital airport [J]. China Civil Engineering Journal, 2005,38(1):40-44.

[4]王國(guó)硯,黃本才,林穎儒,等.基于CQC方法的大跨屋蓋結(jié)構(gòu)隨機(jī)風(fēng)振響應(yīng)計(jì)算[J].空間結(jié)構(gòu),2003,9(4):22-26.

WANG Guo-yan, HUANG Ben-cai, LIN Ying-ru. Solution of wind induced random vibration of large span roof structures based on CQC method[J]. Spatial Structures,2003,9(4):22-26.

[5]Nakayama M, Sasaky Y, Masuda K, et al. An efficient method for selection of vibration modes contributory to wind response on dome like roofs[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1998, 73(1):31-43.

[6]Wu Y, Zhang J S, Chen B, et al. Identification of the dominant vibration modes of single-layer reticulated shells under wind action[J]. International Journal of Space Structures, 2007, 22(2):123- 132.

[7]何艷麗,董石麟,龔景海. 空間網(wǎng)格結(jié)構(gòu)頻域風(fēng)振響應(yīng)分析模態(tài)補(bǔ)償法[J].工程力學(xué),2002,19(4):1-6.

HE Yan-li, DONG Shi-lin, GONG Jing-hai. Wind-induced response of spatial structures with mode compensation in frequency domain[J].Engineering Mechanics,2002,19(4):1-6.

[8]胡繼軍,網(wǎng)殼風(fēng)振及控制研究[D].上海:上海交通大學(xué),2000.

[9]Mataki Y, Iwasa Y, Fukao Y, et al. Wind induced response of low-profile cable-reinforced air-supported structures [J]. Journal of Wind Engineering and Industrial Aerodynamics, 1988, 29:253-262.

[10]鄧華,董石麟,何艷麗,等.深圳游泳跳水館主館屋蓋結(jié)構(gòu)分析及風(fēng)振響應(yīng)計(jì)算[J].建筑結(jié)構(gòu)學(xué)報(bào), 2004, 25(2):72-78.

DENG Hua, DONG Shi-lin, HE Yan-li, et al. Behavior analysis and wind-induced vibration computationof mail hall roof structure of shenzhen swimming and diving hall[J]. Journal of Building Structures,2004,25(2):72-78.

[11]田玉基,楊慶山,范重,等.國(guó)家體育場(chǎng)大跨度屋蓋結(jié)構(gòu)風(fēng)振系數(shù)研究 [J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2007, 28(2): 36-41.

TIAN Yu-ji, YANG Qing-shan, FAN Zhong, et al. Gust factors for large-span roof of the national stadium[J]. Journal of Building Structures,2007,28(2):36-41.

[12]Ashraf Ali M, Gould P L. On the resonant component of the response of single degree of freedom systems under wind loading [J]. Engineering Structures, 1985, 7(4): 280-292.

[13]陳波, 武岳, 沈世釗. 背景響應(yīng)、共振響應(yīng)定義及其相關(guān)性分析方法 [J]. 振動(dòng)工程學(xué)報(bào), 2008, 21(2):140-145.

CHEN Bo, WU Yue, SHEN Shi-zhao. Definitions and correlation analyses for background response and resonant response[J]. Journal of VibrationEngineering, 2008, 21(2):140-145.

[14]田玉基,楊慶山. 大跨度屋蓋結(jié)構(gòu)脈動(dòng)風(fēng)振響應(yīng)的參與振型[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào),2009(6):146-149.

TIAN Yu-ji, YANG Qing-shan. Dominant modes of fluctuating wind-induced response for large-span roof [J]. Journal of Harbin Institute of Technology,2009(6):146-149.

[15]方同.工程隨機(jī)振動(dòng)[M].北京:國(guó)防工業(yè)出版社,1995.

[16]舒新玲, 周岱.風(fēng)速時(shí)程AR模型及其快速實(shí)現(xiàn)[J].空間結(jié)構(gòu), 2003, 12(9): 27-32.

SHU Xin-ling, ZHOU Dai. AR model of wind speed timeseries and its rapid implementation[J].Spatial Structures,2003,12(9):27-32.

第一作者劉均男，副教授，1976年8月生