宋波++殷炳帥++曹謙

摘要：以某電廠660 MW吸收塔結(jié)構(gòu)為研究對象，利用ADINA及FLUENT有限元軟件對吸收塔結(jié)構(gòu)開展了數(shù)值風(fēng)洞分析。研究結(jié)果表明：對于該吸收塔結(jié)構(gòu)，NE風(fēng)向與180°風(fēng)向?qū)儆谧畈焕L(fēng)向，0°風(fēng)向時結(jié)構(gòu)相對較安全；隨著風(fēng)速增大，加速度響應(yīng)隨塔高呈增大趨勢，應(yīng)力呈減小趨勢，塔體下部開口位置有應(yīng)力集中現(xiàn)象，且風(fēng)速對結(jié)構(gòu)應(yīng)力和位移的影響較??；多塔結(jié)構(gòu)之間有較強(qiáng)的相互干擾現(xiàn)象，風(fēng)場上游結(jié)構(gòu)的周圍風(fēng)場受下游結(jié)構(gòu)影響較小，而下游結(jié)構(gòu)周圍風(fēng)場受上游結(jié)構(gòu)的影響較大。

關(guān)鍵詞：吸收塔；最不利風(fēng)向；風(fēng)致動力響應(yīng)；多塔結(jié)構(gòu)

中圖分類號：TU473.1 文獻(xiàn)標(biāo)志碼：A

0 引 言

各國統(tǒng)計數(shù)據(jù)表明，風(fēng)災(zāi)是造成重大生命財產(chǎn)安全損失最嚴(yán)重的自然災(zāi)害之一，故建筑結(jié)構(gòu)在風(fēng)荷載作用下的響應(yīng)問題研究不可忽略。

結(jié)構(gòu)風(fēng)荷載響應(yīng)的研究實質(zhì)上是分析結(jié)構(gòu)的風(fēng)振，是一個空氣彈性力學(xué)問題，實際工程中主要集中在3個方面：顫振、抖振和拉索的風(fēng)振[1]。風(fēng)對結(jié)構(gòu)的作用一般就是考慮順風(fēng)力、橫風(fēng)力及扭力矩3種，所以結(jié)構(gòu)風(fēng)荷載響應(yīng)研究也就主要集中在這3個方面。結(jié)構(gòu)順風(fēng)向響應(yīng)計算的理論己經(jīng)比較成熟。Davenport[2]在他自己提出的近地風(fēng)紊流模型的基礎(chǔ)上建立了計算高層建筑順風(fēng)向響應(yīng)的振風(fēng)因子法[3]，為計算房屋建筑的紊流風(fēng)效應(yīng)（抖振）奠定了基礎(chǔ)。此后的研究工作都是在Davenport提出的方法基礎(chǔ)上考慮一些其他因素建立起來的，如Simiu[45]考慮了迎風(fēng)面與背風(fēng)面風(fēng)壓相關(guān)以及風(fēng)速譜隨高度變化，Solari[6]也考慮了風(fēng)速譜隨高度的變化。

當(dāng)建筑的結(jié)構(gòu)、質(zhì)量及受力關(guān)于中心軸對稱時，結(jié)構(gòu)的順風(fēng)向、橫風(fēng)向及扭轉(zhuǎn)響應(yīng)可以視為相互獨立而分別計算[7]。實際中的結(jié)構(gòu)多是不對稱的，因而結(jié)構(gòu)順風(fēng)向、橫風(fēng)向與扭轉(zhuǎn)風(fēng)振響應(yīng)相互藕合，這個三維藕合振動響應(yīng)研究是風(fēng)工程界尚在研究的一個課題。由于風(fēng)洞試驗的長周期性及條件的局限性，風(fēng)荷載數(shù)值模擬開始作為風(fēng)洞試驗的補(bǔ)充手段或單獨出現(xiàn)在實際工程設(shè)計中，人們已經(jīng)開始利用數(shù)值風(fēng)洞對各種復(fù)雜建筑進(jìn)行表面風(fēng)荷載計算。1997年Stelvam[8]基于LES湍流模型，考慮3種不同的來流條件，利用數(shù)值風(fēng)洞法對德克薩斯科學(xué)研究建筑進(jìn)行了數(shù)值模擬，得到的風(fēng)壓平均值與試驗得到的結(jié)果吻合良好，但峰值風(fēng)壓只有根據(jù)實測數(shù)據(jù)生成脈動來流的數(shù)值模擬結(jié)果才和實測結(jié)果吻合較好。劉若斐等[9]分析雙塔之間的干擾作用，得出根據(jù)不同風(fēng)向角可分為尾流直接作用、尾流和臨近干擾共同作用、臨近干擾作用3個區(qū)域。曹輝[10]采用數(shù)值風(fēng)洞與風(fēng)洞試驗相結(jié)合的方法，對尚信國際大廈進(jìn)行了風(fēng)荷載特性研究，全面分析了典型超高層建筑表面局部體型系數(shù)的分布規(guī)律及其隨風(fēng)向的變化規(guī)律，探討了典型的矩形超高層建筑在不同風(fēng)向下的風(fēng)荷載確定方法。

總體上看，目前各國對風(fēng)荷載導(dǎo)致結(jié)構(gòu)響應(yīng)的研究主要集中在高層、超高層建筑物，而關(guān)于薄壁類高聳殼體吸收塔結(jié)構(gòu)在風(fēng)荷載作用下的動力響應(yīng)研究較少。本文中依托的脫硫塔結(jié)構(gòu)塔體下部具有大型開口，頂部布置有濕式電除塵裝置，結(jié)構(gòu)、荷載不對稱，對風(fēng)荷載的敏感度更高。鑒于此，本文采用ADINA及FLUENT有限元軟件對頂部布置有濕式電除塵裝置的脫硫吸收塔結(jié)構(gòu)開展了數(shù)值風(fēng)洞分析，進(jìn)行了不同風(fēng)向作用下風(fēng)荷載對吸收塔結(jié)構(gòu)的影響研究并明確最不利輸入風(fēng)向。同時為充分研究吸收塔結(jié)構(gòu)與其周邊結(jié)構(gòu)的相互作用，更進(jìn)一步分析實際環(huán)境的風(fēng)場情況，利用CFD數(shù)值模擬軟件FLUENT對吸收塔結(jié)構(gòu)及周邊結(jié)構(gòu)的風(fēng)場進(jìn)行了分析研究，為實際工程中吸收塔類結(jié)構(gòu)的安裝布置提供參考。

1 平均風(fēng)作用下最不利風(fēng)向的確定

本文中所研究的脫硫吸收塔結(jié)構(gòu)為貴州省畢節(jié)地區(qū)某電廠2×660 MW級新建煙氣脫硫塔，結(jié)構(gòu)簡圖如圖1所示。根據(jù)《建筑結(jié)構(gòu)荷載規(guī)范》（GB 50009—2001）[11]可知，畢節(jié)地區(qū)50年重現(xiàn)期基本風(fēng)壓為0.3 kPa，由風(fēng)速風(fēng)壓換算公式計算可知當(dāng)?shù)?0 m高度處10 min內(nèi)平均風(fēng)速為22.4 m·s-1。

吸收塔結(jié)構(gòu)塔身、進(jìn)出口煙道的材料為Q235鋼，密度為7 850 kg·m-3，泊松比為0.3，彈性模量為2.06×1011 Pa，線性膨脹系數(shù)為12.5×10-6 ℃-1，半徑r=9.7 m。吸收塔結(jié)構(gòu)頂部濕式電除塵裝置外殼部分的直徑為13.3 m，材料為玻璃鋼，彈性模量為3.0×1010 Pa，泊松比為0.3，密度為2 100 kg·m-3，熱膨脹系數(shù)為5.0×10-6 ℃-1。

當(dāng)?shù)氐娘L(fēng)向是結(jié)構(gòu)合理布局的重要參考數(shù)據(jù)，依照當(dāng)?shù)氐娘L(fēng)玫瑰圖確定合理的結(jié)構(gòu)布局是工業(yè)設(shè)計的主要方法。風(fēng)玫瑰圖一般可反映2個重要信息：風(fēng)向、風(fēng)頻。吸收塔結(jié)構(gòu)所處貴州省畢節(jié)地區(qū)織金縣的風(fēng)玫瑰圖如圖2所示。按照結(jié)構(gòu)的前期設(shè)計，吸收塔結(jié)構(gòu)的方向布置如圖3所示。

由圖2可以看出，該地區(qū)東北風(fēng)（NE）的頻率較高，其次是北偏東風(fēng)（NNE）。吸收塔結(jié)構(gòu)的方向布置形式是否合理，需要通過平均風(fēng)作用下結(jié)構(gòu)的響應(yīng)進(jìn)行分析，因而本文展開了不同風(fēng)向作用下吸收塔結(jié)構(gòu)的數(shù)值風(fēng)洞試驗，探求在實際風(fēng)向頻率最高輸入角度情況下結(jié)構(gòu)的風(fēng)致動力響應(yīng)，以及不同方向角輸入下結(jié)構(gòu)的最不利風(fēng)向輸入，并明確風(fēng)向輸入角度與結(jié)構(gòu)響應(yīng)的規(guī)律，為后續(xù)展開的脈動風(fēng)作用下結(jié)構(gòu)的動力響應(yīng)研究提供參考依據(jù)。

1.1 吸收塔結(jié)構(gòu)數(shù)值風(fēng)洞模型

根據(jù)《建筑結(jié)構(gòu)荷載規(guī)范》（GB 50009—2001），吸收塔所在畢節(jié)地區(qū)50年重現(xiàn)期基本風(fēng)壓為0.3 kPa，畢節(jié)地區(qū)屬于B類地區(qū)，其地面粗糙度系數(shù)取為0.16。風(fēng)壓風(fēng)速關(guān)系式為

ω0=12ρv20

（1）

式中：ω0為風(fēng)壓；ρ為大氣密度；v0為平均風(fēng)速。

根據(jù)式（1）可計算得到當(dāng)?shù)?0 m高度處10 min內(nèi)的平均風(fēng)速為22.4 m·s-1。同樣參照該規(guī)范可知畢節(jié)地區(qū)10年重現(xiàn)期基本風(fēng)速為18.3 m·s-1，100年重現(xiàn)期基本風(fēng)速為24.2 m·s-1。

由于實際工程結(jié)構(gòu)完成后，將在結(jié)構(gòu)外層附加保溫層，結(jié)構(gòu)的外表面趨于光滑。另外由于濕式電除塵裝置、煙道及AFT漿液收集盤的存在，導(dǎo)致結(jié)構(gòu)的形式為高聳異形結(jié)構(gòu)，在數(shù)值風(fēng)洞分析過程中對流場與結(jié)構(gòu)面的耦合造成一定困難，因此對結(jié)構(gòu)進(jìn)行合理的簡化是有必要的。由于煙道開口，將導(dǎo)致結(jié)構(gòu)與風(fēng)場在煙道進(jìn)出氣口處無法耦合，因此，在這兩處進(jìn)行了封頭處理，采用低強(qiáng)度、輕質(zhì)量的薄板封頭，防止風(fēng)進(jìn)入塔內(nèi)。

風(fēng)場模型采用了不壓縮性流體來模擬，以此建立了包圍著吸收塔結(jié)構(gòu)的三維流體單元模型，同時在風(fēng)入口邊界處施加了速度條件，風(fēng)速的大小根據(jù)前文求解得出的風(fēng)速確定，并在風(fēng)場的側(cè)面設(shè)置了固壁（Wall）條件。將風(fēng)場與結(jié)構(gòu)模型的接觸面設(shè)置為流固耦合界面，繼而展開后續(xù)相應(yīng)的數(shù)值分析。根據(jù)上述分析，建立的結(jié)構(gòu)模型及數(shù)值風(fēng)洞模型如圖4所示。

1.2 風(fēng)向輸入變化對結(jié)構(gòu)的影響

貴州省畢節(jié)地區(qū)50年重現(xiàn)期的基本風(fēng)速為22.4 m·s-1，以此風(fēng)速分別從NE風(fēng)向（該地實際頻率最高的風(fēng)向）、0°風(fēng)向、90°風(fēng)向及180°風(fēng)向輸入平均風(fēng)，分析4個風(fēng)向工況下結(jié)構(gòu)最大應(yīng)力位移以及風(fēng)場的風(fēng)速響應(yīng)，研究不同風(fēng)向輸入對結(jié)構(gòu)的影響，同時為實際結(jié)構(gòu)的方向布置提供一定的參考依據(jù)，并為后續(xù)的脈動風(fēng)研究提供基礎(chǔ)。

NE風(fēng)向是畢節(jié)地區(qū)最高頻率風(fēng)向，要展開靜風(fēng)作用下結(jié)構(gòu)的安全性能分析首先需要進(jìn)行NE風(fēng)向的分析。在該風(fēng)向輸入下，風(fēng)場Wall邊界的速度云圖如圖5（a）所示，風(fēng)速矢量圖如圖5（b）所示。沿平行于風(fēng)向的塔體中心做豎向剖面，其速度云圖如圖6（a）所示；沿濕式電除塵器與下部塔體連接部位的水平面做剖面，其風(fēng)速矢量圖如圖6（b）所示。

由圖5和圖6可以看出，NE風(fēng)向時，結(jié)構(gòu)的迎風(fēng)面風(fēng)速較大，結(jié)構(gòu)的背風(fēng)面?zhèn)蕊L(fēng)速減小，結(jié)構(gòu)遠(yuǎn)處風(fēng)場的風(fēng)速受結(jié)構(gòu)的干擾相對較小。

為分析不同風(fēng)向作用下結(jié)構(gòu)的動力響應(yīng)，分別從0°風(fēng)向、90°風(fēng)向及NE風(fēng)向輸入平均風(fēng)進(jìn)行分析。為方便對比分析結(jié)構(gòu)，將不同風(fēng)向作用下結(jié)構(gòu)應(yīng)力、位移和加速度隨高度變化趨勢表示出來，如圖7所示，并將不同風(fēng)向作用下結(jié)構(gòu)的最大動力響應(yīng)匯總，如表1所示。

由圖7可以看出，在不同風(fēng)向輸入下，結(jié)構(gòu)的最大位移響應(yīng)出現(xiàn)在180°風(fēng)向（與NE風(fēng)向動力響應(yīng)較接近，為使模擬結(jié)果可視性更強(qiáng)，用NE風(fēng)向代替180°風(fēng)向），最大應(yīng)力響應(yīng)出現(xiàn)在NE風(fēng)向，結(jié)構(gòu)的最小應(yīng)力及位移響應(yīng)均出現(xiàn)在0°風(fēng)向。因此可以判斷，NE風(fēng)向與180°風(fēng)向?qū)τ诮Y(jié)構(gòu)而言為相對較不利風(fēng)向，0°風(fēng)向時結(jié)構(gòu)相對安全，這為吸收塔結(jié)構(gòu)方向布置提供了一定參考。

由表1可以看出：結(jié)構(gòu)的最大應(yīng)力為60.08 MPa，出現(xiàn)在進(jìn)氣口煙道下端，該值遠(yuǎn)小于Q235鋼材允許應(yīng)力113 MPa；結(jié)構(gòu)的最大位移為4.6 cm，遠(yuǎn)小于《高聳結(jié)構(gòu)設(shè)計規(guī)范》允許值[12]。2 脈動風(fēng)作用下結(jié)構(gòu)的動力響應(yīng)

通常把風(fēng)作用的風(fēng)壓平均值當(dāng)做穩(wěn)定風(fēng)壓或者平均風(fēng)壓，而實際風(fēng)壓是在平均風(fēng)壓值上下波動的。在平均風(fēng)壓作用下，建筑物產(chǎn)生一定的側(cè)向位移，而在脈動風(fēng)壓作用下，建筑物在該側(cè)向位移值附近振動。尤其對高度大、剛度小的高層建筑物，脈動風(fēng)壓將引起不可忽視的動力效應(yīng)。當(dāng)建筑物高度大于30 m或高寬比大于1.5時，包括塔架、鋼構(gòu)架、煙囪等高聳結(jié)構(gòu)均需考慮風(fēng)振。本文中吸收塔結(jié)構(gòu)整體高度達(dá)73 m，高徑比大于1.5，必須進(jìn)行風(fēng)振分析。參考柳美玉[13]基于風(fēng)洞試驗研究計算的50年、100年周期設(shè)計風(fēng)荷載，進(jìn)行了10年一遇、50年一遇及100年一遇脈動風(fēng)分析，同時基于第1.2節(jié)中結(jié)論，從結(jié)構(gòu)的不利風(fēng)向180°風(fēng)向輸入脈動風(fēng)。

2.1 脈動風(fēng)速譜模擬

展開吸收塔結(jié)構(gòu)的風(fēng)致動力響應(yīng)分析，具備較強(qiáng)適用性及代表性的風(fēng)速時程曲線是基本前提。本文應(yīng)用MATLAB軟件編譯適用于吸收塔結(jié)構(gòu)的風(fēng)速時程曲線以模擬風(fēng)速，根據(jù)前文所述畢節(jié)地區(qū)10年一遇、50年一遇及100年一遇基本風(fēng)速分別為18.3，22.4，24.2 m·s-1，頻率的取值范圍為0.001～1 Hz，根據(jù)規(guī)范采用Davenport譜模擬當(dāng)?shù)氐拿}動風(fēng)速譜，得到的貴州畢節(jié)地區(qū)10 m高度處10年一遇、50年一遇及100年一遇脈動風(fēng)時程曲線如圖8所示。同時為驗證應(yīng)用AR自回歸法模擬脈動風(fēng)速譜的準(zhǔn)確性，本文將模擬譜和目標(biāo)譜（Davenport譜）進(jìn)行對比，如圖9所示。

由圖8可以看出，在3種不同基本風(fēng)速下，脈動風(fēng)均圍繞著基本風(fēng)速波動，且隨著風(fēng)速增加，波形變化越發(fā)明顯。同時由圖9可以看出，應(yīng)用AR自回歸法模擬的脈動風(fēng)速譜與目標(biāo)譜（Davenport譜）的擬合情況良好，說明應(yīng)用AR自回歸法模擬所得的脈動風(fēng)速譜具備準(zhǔn)確性，且未發(fā)生脈動風(fēng)速時程失真的現(xiàn)象，與實際工況比較接近，模擬所得的脈動風(fēng)速譜可用于下一步的研究中。

2.2 不同風(fēng)速下吸收塔結(jié)構(gòu)的動力響應(yīng)分析

根據(jù)不同種類設(shè)備使用年限的相關(guān)規(guī)定，電廠設(shè)備的使用年限以40年為基準(zhǔn)，所以分析在吸收塔服役期間10年一遇罕見風(fēng)速對結(jié)構(gòu)的影響具有重要的工程意義。參照第2.1節(jié)模擬得到的10年一遇脈動風(fēng)速譜，由180°風(fēng)向輸入至結(jié)構(gòu)中，風(fēng)場的風(fēng)速云圖如圖10所示。

由圖10可以看出，由于在結(jié)構(gòu)的180°風(fēng)向施加初始風(fēng)速邊界條件，風(fēng)場中空氣流動逐漸由施加風(fēng)速一側(cè)向另一側(cè)移動，最終風(fēng)場中風(fēng)速變化趨于與施加風(fēng)速變化一致的形式，由于風(fēng)場的總長度為200 m，該過程大約需要10～20 s。同時由于風(fēng)場的前10～20 s風(fēng)速較小，結(jié)構(gòu)的動力響應(yīng)相對也較小，為了更加直觀明顯地表示出時程曲線的變化過程，特將時程曲線前10 s的曲線段刪除，則吸收塔結(jié)構(gòu)在脈動風(fēng)作用下位移、加速度及應(yīng)力的動力響應(yīng)最大位置的時程曲線如圖11所示。

由圖11可以看出：結(jié)構(gòu)的最大位移為45.99 mm，出現(xiàn)在下部塔體與上部濕式電除塵結(jié)構(gòu)連接的變截面處；結(jié)構(gòu)的最大加速度響應(yīng)為10.6 mm·s-2；結(jié)構(gòu)的最大應(yīng)力響應(yīng)為48.13 MPa。說明在10年一遇風(fēng)速下，結(jié)構(gòu)的最大位移響應(yīng)未超過74 cm，最大應(yīng)力響應(yīng)均未超過允許值113 MP，吸收塔結(jié)構(gòu)是安全的。

按照上述10年一遇罕見風(fēng)速對結(jié)構(gòu)影響的分析過程，分別用50年一遇及100年一遇風(fēng)速譜進(jìn)行分析（過程不做贅述）。為了對比分析不同脈動風(fēng)風(fēng)速作用下結(jié)構(gòu)的動力響應(yīng)，分別沿結(jié)構(gòu)的高度方向提取了位移、加速度以及應(yīng)力響應(yīng)，如圖12所示。

由圖12（a）可以看出：隨著結(jié)構(gòu)高度的增加，在風(fēng)荷載作用下結(jié)構(gòu)的位移響應(yīng)逐漸增大，但在下部塔體與上部濕除結(jié)構(gòu)連接的變截面處位移出現(xiàn)突增，說明在風(fēng)荷載作用下塔體整體將呈現(xiàn)上部位移大于下部位移的趨勢，但由于吸收塔結(jié)構(gòu)是“頭重腳輕”類結(jié)構(gòu)，濕式電除塵裝置的中部有質(zhì)量較大的支撐梁，風(fēng)荷載作用下該部位相對其下部塔段要更穩(wěn)定，位移響應(yīng)要相對較小，從而導(dǎo)致了下部塔體與上部濕除結(jié)構(gòu)連接的變截面處位移相對較大；隨著風(fēng)速增大，塔體頂部位移響應(yīng)的差異較塔底明顯。由圖12（b）可以看出，隨著結(jié)構(gòu)的增高，加速度響應(yīng)逐漸增大。圖12（c）可以看出：塔體底部應(yīng)力相對較大，隨結(jié)構(gòu)高度增加，應(yīng)力整體呈減小趨勢，但在煙氣入口部位及濕除結(jié)構(gòu)連接的變截面處應(yīng)力集中；在不同風(fēng)速下結(jié)構(gòu)應(yīng)力響應(yīng)的差異在塔體底部較為明顯。3 多塔結(jié)構(gòu)數(shù)值風(fēng)洞試驗

由于多塔結(jié)構(gòu)高度大，且相距較近，宜考慮群體效應(yīng)[14]。為了深入研究吸收塔結(jié)構(gòu)與其周邊結(jié)構(gòu)的相互作用，進(jìn)一步分析實際環(huán)境的風(fēng)場情況，利用CFD數(shù)值模擬[15]軟件FLUENT對吸收塔結(jié)構(gòu)及其周邊結(jié)構(gòu)的風(fēng)場進(jìn)行分析研究。根據(jù)實際工程資料，吸收塔分為1#塔體與2#塔體，中間布置煙囪，其平面布置情況如圖13所示。

利用GAMBIT軟件建立了整體結(jié)構(gòu)的有限元模型，如圖14所示。根據(jù)畢節(jié)地區(qū)風(fēng)速玫瑰圖可知當(dāng)?shù)仡l率較高風(fēng)向為東北風(fēng)，因此風(fēng)場的長度方向沿東北西南（NESW）方向，其總長度為800 m，寬度為600 m，高度為400 m。計算風(fēng)場入口采用的是速度進(jìn)流邊界條件，并采用完全發(fā)展出流邊界作為出口邊界條件；結(jié)構(gòu)的表面及地面采用的是無滑移壁面條件；風(fēng)場的頂部及兩側(cè)采用的是對稱邊界條件。同時根據(jù)前文計算分析可知，吸收塔所在地50年重現(xiàn)期的基本風(fēng)壓為0.3 kPa，則基本風(fēng)速為22.4 m·s-1，該地區(qū)屬于B類地區(qū)，其地面粗糙度系數(shù)取為0.16。

本文數(shù)值風(fēng)洞分析模型的水力直徑為480 m，水力直徑是在管內(nèi)流動中引入的，其目的是給非圓管流動的分析模型取一個合適的特征長度，以此來計算其雷諾數(shù)。水力直徑為4倍濕橫截面面積與濕圓周長的比值，即水力直徑D=4A/P，A為濕橫截面面積，P為濕圓周長。模型的入口邊長分別為400，600 m，因此A=240 000 m2，P=1 000 m，則所求水力直徑為480 m；根據(jù)經(jīng)驗，本文模型的湍流強(qiáng)度缺省參數(shù)設(shè)置為5。湍流強(qiáng)度表示為湍流強(qiáng)度漲落標(biāo)準(zhǔn)差與平均速度的比值，其值是衡量湍流強(qiáng)弱的指標(biāo)。本文模型的求解精度為0.001，求解過程至第62步開始收斂。

經(jīng)計算分析，得到模型的風(fēng)壓云圖如圖15所示，同時為了顯示不同高度處結(jié)構(gòu)的各類響應(yīng)值，定義了的不同高度面，分別為0，14（進(jìn)氣口中部位置），49（濕除結(jié)構(gòu)與下部塔體連接變截面處），100 m（煙囪中部）。對上述3個高度處的風(fēng)壓云圖進(jìn)行分析可知：3個結(jié)構(gòu)在迎風(fēng)面上的最大風(fēng)壓基本相近，最大風(fēng)壓達(dá)426 Pa，而在結(jié)構(gòu)的背風(fēng)側(cè)由于2#吸收塔結(jié)構(gòu)處于迎風(fēng)層的最前端，其背風(fēng)側(cè)風(fēng)壓值較煙囪及1#吸收塔結(jié)構(gòu)的背風(fēng)側(cè)風(fēng)壓值更大，最小風(fēng)壓為-838 Pa，出現(xiàn)在煙氣入口煙道下方。同時，可以看到由于3個結(jié)構(gòu)對風(fēng)流向的相互影響，負(fù)壓區(qū)域不是直接分布于結(jié)構(gòu)的背風(fēng)側(cè)，為了進(jìn)一步分析其原因及效果，分別提取了14 m及49 m高度處結(jié)構(gòu)附近的風(fēng)速矢量圖，如圖16所示。

可以看出，在NE風(fēng)向下，3個結(jié)構(gòu)對風(fēng)場流動存在相互干擾的現(xiàn)象。氣流在流經(jīng)2#吸收塔之后，由結(jié)構(gòu)兩側(cè)流向背風(fēng)側(cè)，未形成明顯渦流現(xiàn)象，但由圖16可以看出，風(fēng)速在背風(fēng)側(cè)明顯減小，由貝努力定理可知氣流對2#吸收塔的壓力增大；之后該氣流經(jīng)由煙囪結(jié)構(gòu)的一側(cè)流過，同樣與經(jīng)由煙囪的氣流在其背風(fēng)側(cè)相聚，依次經(jīng)由1#吸收塔結(jié)構(gòu)流過，形成的流線圖如圖17所示（14 m高度處）。

通過上述分析可以看出，3個結(jié)構(gòu)之間存在較強(qiáng)的相互干擾現(xiàn)象，但是上游結(jié)構(gòu)的周圍風(fēng)場受下游結(jié)構(gòu)影響不大，在單獨計算上游結(jié)構(gòu)時可以不考慮干擾效應(yīng)。下游結(jié)構(gòu)周圍風(fēng)場受上游結(jié)構(gòu)的影響較大，在設(shè)計分析下游結(jié)構(gòu)時必須考慮上游結(jié)構(gòu)的影響。

4 結(jié) 語

（1）通過對比分析不同風(fēng)向作用下結(jié)構(gòu)的風(fēng)致動力響應(yīng)得出：本文結(jié)構(gòu)的最大位移響應(yīng)出現(xiàn)在180°風(fēng)向，最大應(yīng)力響應(yīng)出現(xiàn)在NE風(fēng)向，結(jié)構(gòu)的最小應(yīng)力及位移響應(yīng)均出現(xiàn)在0°風(fēng)向。因此可以判斷，NE風(fēng)向與180°風(fēng)向?qū)τ诮Y(jié)構(gòu)而言為相對不利風(fēng)向，0°風(fēng)向時結(jié)構(gòu)相對安全，此項計算對實際工程中確定不對稱吸收塔結(jié)構(gòu)的安裝方向具有一定借鑒意義。

（2）隨著結(jié)構(gòu)高度的增加，在風(fēng)荷載作用下結(jié)構(gòu)的位移響應(yīng)逐漸增大，在濕除結(jié)構(gòu)與塔體連接處由于材料的不連續(xù)性使得結(jié)構(gòu)位移產(chǎn)生突變，所以實際工程中要注意變截面、不同材料銜接處的加強(qiáng)。

（3）風(fēng)速的增大引起結(jié)構(gòu)應(yīng)力、加速度及頂部位移的增大，但是影響較小。相比較而言，在實際工程中應(yīng)該重點考慮風(fēng)向的影響。

（4）多塔結(jié)構(gòu)的數(shù)值風(fēng)洞試驗分析表明，3個結(jié)構(gòu)之間存在較強(qiáng)的相互干擾現(xiàn)象，但是上游結(jié)構(gòu)的周圍風(fēng)場受下游結(jié)構(gòu)影響不大，在單獨計算上游結(jié)構(gòu)時可以不考慮干擾效應(yīng)。下游結(jié)構(gòu)周圍風(fēng)場受上游結(jié)構(gòu)的影響較大，在設(shè)計分析下游結(jié)構(gòu)時必須考慮上游結(jié)構(gòu)的影響。

（5）對于吸收塔結(jié)構(gòu)動力響應(yīng)規(guī)律的研究，本文只考慮了風(fēng)荷載，但是實際工程中地震荷載的影響不可忽略，還需對風(fēng)荷載與地震荷載共同作用下吸收塔的動力特性進(jìn)行分析。

參考文獻(xiàn)：

References：

[1] 項海帆.結(jié)構(gòu)風(fēng)工程研究的現(xiàn)狀和展望[J].振動工程學(xué)報，1997，10（3）：258263.

XIANG Haifan.State of the Art and Prospect in Studies on Structural Wind Engineering[J].Journal of Vibration Engineering，1997，10（3）：258263.

[2]DAVENPORT A G.The Application of Statistical Concepts to the Wind Loads of Structures[J].Proceedings of the Institute of Civil Engineering，1961，19（4）：449472.

[3]DAVENPORT A G.Gust Loading Factors[J].Journal of the Structural Division，1967，93（3）：1134.

[4]SIMIU E.Gust Factors and Alongwind Pressure Correlation[J].Journal of the Structural Division，1973，99（4）：773783.

[5]SIMIU E.Equivalent Static Wind Loads for Tall Building Design[J].Journal of the Structural Division，1976，102（4）：719737.

[6]SOLARI G.Alongwind Response Estimation：Closed Form Solution[J].Journal of the Structural Division，1982，108（1）：225244.

[7]BENEKE D L，KWOK K C S.Aerodynamic Effect of Wind Induced Torsion on Tall Buildings[J].Journal of Wind Engineering and Industrial Aerodynamics，1993，50：271280.

[8]STELVAM R P.Computation of Pressures on Texas Tech University Building Using Large Eddy Simulation[J].Journal of Wind Engineering and Industrial Aerodynamics，1997，6768：647657.

[9]劉若斐，沈國輝，孫炳楠.大型冷卻塔風(fēng)荷載的數(shù)值模擬研究[J].工程力學(xué)，2006，23（增1）：177183.

LIU Ruofei，SHEN Guohui，SUN Bingnan.Numerical Simulation Study of Wind Load on Large Hyperbolic Cooling Tower[J].Engineering Mechanics，2006，23（S1）：177183.

[10]曹 輝.尚信國際大廈風(fēng)洞試驗及風(fēng)荷載特性研究[D].重慶：重慶大學(xué)，2012.

CAO Hui.Research on Wind Tunnel Test and Wind Loads of Shangxin International Building[D].Chongqing：Chongqing University，2012.

[11]GB 50009—2001，建筑結(jié)構(gòu)荷載規(guī)范[S].

GB 50009—2001，Load Code for the Design of Building Structures[S].

[12]GB 50135—2006，高聳結(jié)構(gòu)設(shè)計規(guī)范[S].

GB 50135—2006，Code for Design of Highrising Structures[S].

[13]柳美玉.高層雙塔建筑體型系數(shù)的風(fēng)洞試驗研究及數(shù)值模擬[D].北京：清華大學(xué)，2008.

LIU Meiyu.Study on Wind Load Shape Coefficient of Highrise Doubletower Building[D].Beijing：Tsinghua University，2008.

[14]虞終軍，丁潔民，林禎彬，等.綠地中央廣場雙塔結(jié)構(gòu)抗風(fēng)設(shè)計[J].建筑結(jié)構(gòu)，2013，43（9）：4246.

YU Zhongjun，DING Jiemin.LIN Zhenshan，et al.Windresistant Design on Doubletower of Greenland Central Plaza[J].Building Structure，2013，43（9）：4246.

[15]唐 意，金新陽，楊立國.錯列布置超高層建筑群的干擾效應(yīng)研究[J].土木工程學(xué)報，2012，45（8）：97103.

TANG Yi，JIN Xinyang，YANG Liguo.Study of the Interference Effects of Wind Loads on Tall Buildings in Staggered Arrangement[J].China Civil Engineering Journal，2012，45（8）：97103.