大型冷卻塔結構風致氣動和氣彈效應綜合研究

趙 林,展艷艷,王志男,梁譽文,劉曉鵬,程霄翔,張軍鋒,柯世堂,王小松,曹曙陽,葛耀君

(同濟大學 土木工程防災國家重點實驗室,上海 200092)

大型冷卻塔結構風致氣動和氣彈效應綜合研究

趙 林,展艷艷*,王志男,梁譽文,劉曉鵬,程霄翔,張軍鋒,柯世堂,王小松,曹曙陽,葛耀君

(同濟大學 土木工程防災國家重點實驗室,上海 200092)

大型冷卻塔作為空間薄壁高聳結構,結構風效應突出,在設計中風荷載屬重要的控制因素。本文系統地介紹了同濟大學結構風工程研究團隊在冷卻塔結構抗風方向取得的研究進展,涉及風荷載特征與統計、結構風效應試驗與分析、整體結構多目標優化、災害氣候條件荷載與效應分析和全過程集成軟件平臺開發等5個方面,在冷卻塔結構抗風問題中的風荷載、結構響應和試驗模型設計等方面得出相應結論。

大型冷卻塔;風荷載;結構風效應;多目標優化;災害氣候;全過程分析

0 引 言

大型冷卻塔為空間薄壁高聳殼體結構,對風荷載作用敏感,在設計中風荷載是一個重要的控制因素,結構風荷載的準確評價是其抗風研究的基礎。結合風洞試驗和現場實測手段,研究了塔筒表面內外壓分布規律[1-6],發現了超高雷諾數條件結構表面動態繞流發展規律[7-10],建立了基于計算流體力學方法的冷卻塔結構風荷載計算方法,揭示了復雜群塔干擾條件風荷載分布模式和結構靜、動力荷載響應規律[11],并對原型冷卻塔開展了長期動態風壓實測以驗證試驗和計算結果[12-14]。

大型冷卻塔振型復雜,風振效應突出[15-17]。為從風洞試驗的角度研究結構風效應,提出了等效梁格氣動彈性模型設計理論和方法[18]并開展同步測壓和測振試驗[19],用一致耦合分析法分析冷卻塔風致響應[20-21];建立了計入結構荷載行為[22]的多目標一致原則的等效風荷載分析方法[2,23-25];考慮不同評價指標[26-29],研究雙塔[17,30-31]、四塔[27]、六塔[26]和八塔[28,32]群塔組合干擾效應[33-34];開展了考慮非線性效應的冷卻塔結構屈曲穩定分析。

基于大型冷卻塔全壽命經濟性,提出了考慮各種荷載模式下[35-36]結構穩定和強度安全等綜合效應的塔筒選型多目標優化方法[37-39],優化過程同時考慮復雜荷載組合作用下的結構安全性和工程經濟性。

傳統冷卻塔研究工作中的風洞試驗和計算分析完全基于良態氣候模式(季風等),大型冷卻塔對臺風、龍卷風等特殊極端風作用非常敏感,此類風環境特征不同于良態氣候模式。為此,開展了臺風和龍卷風2類罕遇災害模式作用下的理論、試驗和實測研究[40-42],在臺風脈動特性[43-44]、極值風特性隨機模擬[42]等方面取得進展;分別利用多風扇主動控制風洞[45-48]和自主開發的龍卷風模擬裝置[49]對臺風和龍卷風風場特性進行試驗研究工作,保證了災害氣候條件下冷卻塔環向風壓分布、風振系數及內吸力影響等因素相關技術參數取值正確,滿足在我國沿海風災多發區建設大型冷卻塔的需求。

基于冷卻塔設計預研、施工建造和整體運營等全壽命各個環節,研發了涵蓋結構建模、內力計算、雙向配筋、整體優化和自動繪圖[50]等冷卻塔設計的軟件集成平臺,并可針對多種災害氣候模式(臺風和龍卷風)分析風效應規律和抗風穩定性。

1 風荷載特征與統計

1.1 動態風壓現場實測

冷卻塔風洞試驗模型無法完整模擬雷諾數效應,難于準確再現冷卻塔表面動態風荷載與來流條件、塔群組合等參數間的關系,全尺寸現場實測是對風洞試驗的有力補充和驗證。

研發了全天候動態風壓采集設備,對徐州某電廠冷卻塔(高約167m)進行了塔筒表面動、靜態風壓長期現場觀測[12-13]。實測現場地處廣闊平原(可近似看作B類地貌),在正北與正東方向無任何高大建筑,在正西方向有少部分高大建筑(見圖1),考慮風速風向現場實測數據并結合冷卻塔剛性模型測壓風洞試驗對比,規避了特定風向臨近建筑物的干擾影響[13]。為確定風壓測量時的來流狀態,在測量塔正東向、距地面高度約10m處設置超聲波風速儀。如圖2所示,在冷卻塔130m高度處沿塔筒外表面環向均勻布置36個壓阻式風壓傳感器,在90m高度處外表面環向均勻布置18個壓阻式風壓傳感器。測試期間記錄到了6次風速較大(130m高度處風速介于11～17m/s)且不受周邊高大建筑物干擾的風壓時程。圖3給出了6次不同風速下130m高度處沿環向表面脈動壓力分布曲線及相應的最小二乘平均值擬合曲線,同時與國內外歷次冷卻塔脈動風壓的實測結果進行對比。環向脈動壓力分布曲線可以分為3個區域:迎風區(0～40°)、側風區(40°～120°)和背風區(120°～180°)。迎風區在駐點處的脈動風壓最大,在40°附近達到最小值,此點為平均壓力系數為0的角度;環向脈動壓力的最大值出現在側風區70°～90°的區域內,該點同時為平均壓力分布中最小負壓出現的角度;110°～120°為側風區與尾流區的過渡區域,脈動風壓急劇減小,來流在此區域脫離;背風區的風壓脈動平穩且數值較小。

1.2 塔筒表面風荷載非高斯特性

為了說明不同斷面不同區域測點間風壓的空間相關特性,圖4[25]給出了單體冷卻塔塔底、喉部和塔頂3個典型斷面迎風區、負壓極值區和分離區測點風壓與環向所有測點的相關性結果。由于塔底斷面存在明顯的三維端部效應,測點間的相關性較弱,主要呈現的是高斯分布特性;在中部和塔頂區域,迎風面互相關性衰減極快,測點表現出高斯分布特性,在負壓極值至分離區域,測點間的相關性較強,并存在明顯的大偏斜和高峰態現象,風壓信號表現出非高斯特性[56],進入背風區其測點的相關性衰減,高斯特性又逐漸明顯。

圖5給出了冷卻塔中上部斷面的環向高斯與非高斯區域的劃分圖,簡單直接地反應不同區域的風壓分布特性[25]。圖中冷卻塔表面劃分高斯及非高斯區域的標準為:偏斜值|Cpisk|＞0.2且峰態值|Cpiku|＞3.5的風壓信號為非高斯分布。

1.3 統一脈動風壓公式

根據脈動風壓現場實測和風洞試驗數據所表明的脈動風壓與紊流度之間的相關性,建立了不同紊流度條件下統一的脈動風壓環向分布模式,提出了考慮迎風點脈動風壓系數修正擬合的統一脈動風壓系數均方值環向分布公式:

式中:αk(Pu)為與迎風點脈動風壓系數相關的八項式參數;k=0～7;Pu是迎風點脈動風壓系數;Iu是來流順風向紊流度。擬合參數如表1所示。

表1 脈動風壓分布統一擬合式及擬合參數Table 1 Uniform fitting formula and fitting parameters of fluctuating wind pressure distribution

為了進行直觀的對比,將現有現場實測曲線按照迎風點脈動風壓系數數值大小分為3組:第1組:徐州現場實測,迎風點脈動系數約0.115,來流紊流度約7.32%;第2組:Ruscheweyh[51]/Davenport[52]/Sageau[53]現場實測結果[57],迎風點脈動系數約0.2,來流紊流度約12.52%;第3組:周良茂現場實測,迎風點脈動系數約0.245,來流紊流度約15.2%。分別采用迎風點脈動風壓系數,利用修正脈動風壓公式計算得到的結果,與上述3組現場實測脈動風壓系數進行對比。圖6為徐州現場實測值與統一脈動風壓公式計算結果的對比,可以看出,徐州實測脈動風壓曲線與統一脈動風壓分布類似。圖7為已有現場實測值與統一脈動風壓系數的對比,Ruscheweyh和Davenport實測與統一脈動風壓系數在各個區域較為接近,而Sageau實測值在各個區域均大于其余測量值;周良茂實測脈動風壓曲線與統一脈動風壓系數較為接近,除了90°附近略有差異。綜上所述,除了Sageau實測脈動風壓系數與統一脈動風壓系數差別較大外,徐州冷卻塔脈動風壓實測結果與Davenport[52]、Ruscheweyh[51]、周良茂[54]實測曲線當考慮來流紊流度效應后與統一脈動風壓擬合結果在數值上差別不大,故可以定性地說明來流紊流度的不同是造成實測脈動風壓系數產生差異的主要原因。結合本文推薦的脈動風壓曲線統一擬合式(即公式(1)和(2))解釋了圖3中所示歷史上多次脈動風壓分布存在差異的原因。

2 結構風效應試驗與分析

2.1 等效梁格氣彈模型設計

氣彈模型風洞試驗是研究冷卻塔風振響應的有效手段,氣彈模型等效梁格設計方法針對傳統“連續介質氣彈模型”在物理參數和氣動力參數相似比模擬等方面存在的不足做出相應的改進(見表2),在模型頻率、振型、氣動力參數模擬、風振系數分布和試驗易操作性方面均取得良好的效果[18]。等效梁格方法采用空間縱橫垂直交叉桁梁網狀結構近似模擬連續殼筒體結構動力特性,解決了連續介質模型彎扭剛度與軸向剛度不協調的問題;按照幾何相似比的要求,采用具有可張拉性能的彈性、輕質薄膜整體張貼在鋼骨架外表面模擬實際冷卻塔結構的外形;為實現質量系統的模擬,用銅鉛塊為配重補充不足部分的質量;模型環向和子午向構件采用A3鋼線切割加工,考慮外置彈性蒙皮附加阻尼比效應后中,總阻尼比仍可以控制在3.5%左右,基本滿足鋼筋混凝土冷卻塔5%的阻尼比要求。根據等效梁格方法設計的冷卻塔氣彈模型如圖8所示。由表3可見,等效梁格冷卻塔模型與原型結構振型具有較好的相似性且表面壓力分布與規范值較好吻合。

表2 冷卻塔氣彈模型等效梁格設計方法與傳統方法比較Table 2 Comparison between equivalent beam-net method and traditional method

表3 冷卻塔氣彈模型設計參數與實測值Table 3 Design parameters and measurement values of cooling tower model

2.2 共振響應

冷卻塔結構的脈動風振響應以共振分量為主[16],背景和交叉項分量的貢獻比較接近,一般均在10%以下,最大可達到20%左右。圖9給出了超大塔結構總的脈動風振響應三維分布圖。由圖可見,隨著子午向高度的增大,總脈動風振響應逐漸變大,在80m高度處出現第1個峰值,然后再減小到140m高度后又逐步增大至塔頂,其在子午向上出現2個明顯的峰值;在環向斷面上其響應也出現多個峰值,并且沒有對稱特性;結構的最大脈動風振響應數值達到6.6mm。

2.3 群塔干擾評價準則

為量化風荷載條件冷卻塔群塔干擾效應,形成了眾多干擾效應定義準則,派生出數值結果差異性較大的荷載比例放大因子——群塔比例系數,認識和評判準則存在明顯的差異。為了系統地研究冷卻塔復雜群塔條件荷載干擾效應,本項研究以某超大型冷卻塔為例,基于風洞試驗、結構有限元分析和結構設計配筋方案,系統分析了六塔典型布置條件的群塔干擾效應,詳細比較了不同塔距條件矩陣、菱形布置方案的基于荷載、受力和設計配筋3個比較準則層面的群塔干擾效應。圖10定義和比較了25種指標下的群塔比例系數,可見多種群塔比例系數準則在數值大小和離散性方面存在明顯差異,在描述風荷載最不利來流角度方面具有較好的一致性,其中配筋層面的群塔比例系數離散性顯著小于荷載層面結果[58]。

根據配筋包絡,采取分項群塔比例系數放大規范風壓(見圖11),配筋曲線與按全局統一的群塔比例系數的配筋曲線的對比如圖12所示。由圖可知,按照分項群塔比例系數計算的塔筒配筋曲線在實現配筋包絡的同時減少其它模板位置的配筋富余[59]。推薦多風向角多設計荷載組合獲得的塔筒配筋包絡作為冷卻塔設計參考依據,基于配筋包絡比選的在塔筒高度范圍內變化的分項干擾系數作為工程應用群塔比例系數。

2.4 多目標等效風荷載

多目標等效靜風荷載的實質是通過一定的數值優化方法找到一個“誤差最小”的等效靜風荷載分布模式,同時接近多個響應目標的極值。以喉部斷面環向18個節點(平均間隔20°取點)位移響應極值作為等效目標,得出了對應的多目標等效靜風荷載三維分布圖和等值線圖,如圖13和14所示[25]。由圖可知:多目標等效靜風荷載沒有具體的物理含義,在數值上沒有出現分布很不合理的區域或現象,均在-15～18k N之間;在頂部斷面出現的峰值數量減小,并且數值也更加合理,與實際的等效靜風荷載分布特征接近。

為對比多目標等效靜風荷載和規范給出的單塔表面風荷載分布特征的差別,圖15給出了冷卻塔喉部環向斷面上以典型結點位移響應極值作為等效目標得到的等效風荷載和規范分布曲線圖,并給出了各個典型斷面上隨環向角度變化的等效風荷載擬合多項式。由圖可知,考慮多目標的等效風荷載曲線在分布模式和數值大小上均與規范數值有較大的差別;在不同斷面分布模式不同,不存在絕對的正、負壓區域,并且其峰值也不一定出現在迎風面,因此分斷面表示更為合理。

為驗證多目標等效靜風荷載的精度,將其作為靜力荷載加載到大型冷卻塔結構上,求解相應的風振響應,并與結構實際響應極值進行對比,對比結果如圖16所示[25]。由圖可見,采用多目標等效靜風荷載計算出來的結果與實際響應極值相差很小,誤差控制滿足要求。

3 整體結構多目標優化

塔筒選型多目標優化方法考慮各種荷載組合下結構穩定、強度安全和總體造價等的綜合效應,結合響應面法與梯度搜索法,得到了適用于多種風荷載分布模式的塔型。結構優化選型過程分為2個階段[38]:第一階段確定在不同風荷載分布模式下關于荷載加載方向角的最不利工況,優化目標為PT→PTmax;第二階段在對應風荷載分布模式下對典型結構參量取值構成的可行域內進行整體結構優化選出經濟結構尺寸,優化目標為PT→PTmin,待優化參數如表4和圖17所示。

表4 主要優化變量Table 4 Main optimization variables

選取規范對稱風壓、試驗迎風非對稱風壓、試驗側風非對稱風壓和試驗對稱放大風壓4種風荷載分布模式,在4種風荷載分布模式下分別進行冷卻塔結構優化,得到各荷載模式對應的優化塔型子午線型。優化塔型子午線型的半徑及壁厚隨高度的變化如圖18所示[38]。

從塔筒配筋的角度考察優化塔型1的強度安全性,圖19[38]給出其在4種風荷載分布模式下的塔筒理論配筋量。由圖可知,在規范對稱風壓下,優化塔型1的子午向外側、子午向內側理論配筋量在絕大多數子午向模板都大于其他風荷載分布模式的相應值,環向外側、環向內側理論配筋量至少70%的子午向模板理論配筋量都不小于其他風荷載分布模式的相應值,其余曲線交叉情況大多出現在第80階及以上或第5階及以下的子午向模板。在結構設計中可通過實際局部配筋加強保證每一層子午向模板的強度安全。綜上所述,優化后的塔筒子午向線形1適用于多種風荷載分布模式。

在規范對稱風壓下,優化塔型1的穩定系數相較于初始塔型有所提高,鋼筋造價比降低5.4%,總造價降低19.0%。具體優化效果隨優化迭代步數的變化趨勢如圖20所示。

表5 優化塔型1與初始塔型性能對比Table 5 Comparison between optimized and initial tower

4 災害氣候條件荷載與效應分析

4.1 強臺風環境主動風洞模擬

由于氣象站在進行風速統計時無法區分極值風速是良態風引起還是臺風引起,當某地區良態風的發生權重大于臺風時,就會降低規范所統計的基本風速,若在該地區僅采用規范風荷載進行結構抗風設計而不考慮臺風氣候的影響可能會導致結構潛在的風險。

為了模擬臺風場極值風環境特征,選取了上海地區臺風極值風環境作為主動風洞試驗模擬對象,本次試驗利用了日本宮崎大學的三維多風扇主動控制風洞(見圖21)。臺風氣候下B類工程地貌的平均風、紊流度、紊流積分尺度剖面模擬如圖22所示。

以上海地區登陸的臺風作為例,收集和整理了1949～2005年間該地臺風年鑒統計資料和上海崇明島侯家鎮氣象站的逐時(連續)氣象觀測記錄[41],通過敏感度分析方法選出了對臺風極值風環境貢獻程度較大的參數,采用Monte-Carlo臺風風場隨機模型以越界峰值法和廣義Pareto分布探討了幾類工程場地目標重現期內極值風速預測過程,給出了3類典型工程場地臺風風環境的梯度風高度、平均風剖面、極值風速和陣風因子取值(見表6)。可以看出:臺風氣候A類、B類場地下,梯度風高度普遍低于規范值,平均風速剖面也較規范值陡;同時臺風氣候下的10min平均風速及瞬時風速一般低于規范值。

表6 3種典型工程場地臺風與良態氣候模式風環境參數比較(百年設計風速)Table 6 Parameter comparison between typhoon and normal wind at 3 typical engineering sites

基于某電廠大型冷卻塔現場實測和風洞試驗,發現臺風氣候下的紊流度Iu(TC)與規范風荷載下的紊流度Iu(Code)存在如公式(3)所示的關系:

式中:α值對應A、B、C 3類工程場地分別取1.60、1.48和1.36;3類工程場地下的Iu(Code)分別 為IuA(z)=12% × (z/10)-0.12、IuB(z)=14% ×(z/10)-0.15、IuC(z)=23%×(z/10)-0.22。圖23給出了臺風氣候與規范荷載下3類工程場地的紊流度剖面。圖24(a)和(b)分別給出了冷卻塔喉部區域與迎風駐點的相關系數及與風壓最小點的相關系數,可以看出:臺風氣候下紊流度高于規范風荷載的紊流度,因此相關系數也較大,即紊流度越大,相關系數越強。

4.2 龍卷風條件

同濟大學設計開發的龍卷風模擬器[49]可以通過改變導流板角度得到不同渦流比(不同尺度、不同流動形態)的龍卷風,通過改變風機轉速得到不同轉速的龍卷風,同時可模擬龍卷風水平移動。龍卷風風洞及示意圖如圖25和26所示。由龍卷風風洞模擬器得到的龍卷風三維風速和風壓場[61]如圖27所示。

改變龍卷風渦核與冷卻塔塔筒中心的相對位置關系,測量塔筒內外表面風壓分布。如圖28[62]所示,試驗中龍卷風和冷卻塔模型的相對位置關系通過r/rc量化,其中r代表冷卻塔中心與龍卷風中心之間的距離,rc代表冷卻塔塔筒高度處的龍卷風半徑;沿冷卻塔高度布置3層測壓點,環向每層均勻布置12個,測壓點的布置和編號如圖29所示。定義渦流比,變化范圍為1～12。

以上部測點的平均風壓系數分布為例,給出不同渦流比S、不同相對位置條件下的試驗結果[62]如圖30、31和32所示。圖中橫坐標為測壓點編號,縱坐標為平均壓力系數。由圖可見,在龍卷風特異氣流影響下冷卻塔內、外表面均呈現負壓,與規范風荷載與臺風氣候下的冷卻塔表面風荷載分布存在本質區別;冷卻塔距離龍卷風越近,塔筒內、外表面吸力越大,反之,塔筒內、外表面吸力越小。

5 全過程集成軟件平臺開發

目前國內采用的冷卻塔計算分析及設計程序難于滿足超出規范要求的超大型冷卻塔結構抗風和抗震設計要求,缺少靈活的風洞試驗與現場實測風荷載與國外相關規范規定的荷載輸入條件,不具備整體結構考慮多種荷載組合條件全局優化設計功能,分析結果未與通用商業軟件進行系統的校核,不方便設計人員完成冷卻塔結構設計工作。為適應我國新形式下超大型冷卻塔建設工作的順利展開,基于設計預研、施工建造和整體運營等全壽命各個環節,研發了涵蓋結構建模、內力計算、雙向配筋、整體優化和自動繪圖等冷卻塔設計的軟件集成平臺——同濟風向標Wind Lock軟件,并可針對多種災害氣候模式(臺風和龍卷風)分析風效應規律和抗風穩定性。

Wind Lock功能模塊主要包括:(1)良態與臺風氣候極值風環境模擬與預測;(2)大跨空間結構(冷卻塔)動力及等效風荷載分析;(3)大跨空間結構(冷卻塔)多種荷載組合與設計分析;(4)輔助模塊 (氣動力參數數據庫、復雜群樁特性分析等)。具體分項功能實現示例如圖33所示。

6 結 語

本文在概括介紹同濟大學結構風工程研究團隊在冷卻塔結構抗風方向取得的研究進展的基礎上,重點說明了在風荷載特征與統計、結構風效應試驗與分析、整體結構多目標優化、災害氣候條件荷載與效應分析和全過程集成軟件平臺開發等5個方面取得的代表性研究成果。主要結論如下:

(1)進行了長期的原形冷卻塔動態風壓現場實測,得到了環向脈動風壓分布曲線;發現了塔筒表面風荷載的非高斯分布特性,給出了塔筒環向非高斯區域劃分;建立了考慮紊流度影響的統一脈動風壓擬合公式,對歷次現場實測值均有較好的擬合效果。

(2)提出并應用等效梁格氣彈模型研究冷卻塔結構動力響應,量化了沿塔筒高度脈動響應分布;比較不同指標對群塔干擾效應的評價異同;研究了多目標等效靜風荷載的分布規律及其合理性。

(3)提出針對最不利工況確定、冷卻塔結構參數選取和子午線性改進方面的整體結構多目標優化方法。

(4)研究了臺風、龍卷風等災害氣候條件下冷卻塔結構的風效應,對比良態風條件及災害氣候條件下冷卻塔內、外表面風荷載分布差異。

(5)自主開發了全過程集成軟件平臺Wind-Lock,實現大型冷卻塔結構參數化建模、動靜力分析、整體結構優化分析等。

[1] 李鵬飛,趙林,葛耀君,等.超大型冷卻塔風荷載特性風洞試驗研究[J].工程力學,2008,25(6):60-67.Li P F,Zhao L,Ge Y J,et al.Wind tunnel investigation on wind load characteristics for super large cooling towers[J].Engineering Mechanics,2008,25(6):60-67.

[2] 趙林,李鵬飛,葛耀君.等效靜風荷載下超大型冷卻塔受力性能分析[J].工程力學,2008,25(7):79-86.Zhao L,Li P F,Ge Y J.Numerical investigation on performance of super large cooling towers under equivalent static wind load[J].Engineering Mechanics,2008,25(7):79-86.

[3]Zhao L,Ge Y J.Wind loading characteristics of super-large cooling towers[J].Wind and Structures,2010,13(3):257-273.

[4]Cao S,Wang J,Cao J,et al.Experimental study of wind pressures acting on a cooling tower exposed to stationary tornadolike vortices[J].Journal of Wind Engineering and Industrial Aerodynamics,2015,145:75-86.

[5]Cheng X X,Zhao L,Ge Y J.Multiple loading effects on windinduced static performance of super-large cooling towers[J].International Journal of Structural Stability and Dynamics,2013,13(8):1-13.

[6] 李鵬飛.脈動風特性及其對橋梁主梁斷面的抖振作用研究[D].上海:同濟大學,2007.Li P F.Research on turbulent wind characteristics and its effects on buffeting responses of bridge girder sections[D].Shanghai:Tongji University,2007.

[7] 操金鑫,趙林,葛耀君,等.雙曲線圓截面建筑結構雷諾數效應模擬實踐[J].實驗流體力學,2009,23(4):46-50,55.Cao J X,Zhao L,Ge Y J,et al.Practice on simulation of Reynolds number effects for the hyperbolic circular section structures[J].Journal of Experiments in Fluid Mechanics,2009,23(4):46-50,55.

[8] 董銳,趙林,葛耀君,等.雙曲圓截面冷卻塔壁面粗糙度對其繞流動態特性影響[J].空氣動力學學報,2013,31(2):250-259.Dong R,Zhao L,Ge Y J,et al.Investigation of surface roughness and its influence to flow dynamic characteristics of hyperbolic cooling tower[J].Acta Aerodynamica Sinica,2013,31(2):250-259.

[9] 劉天成,趙林,丁志斌.圓形截面冷卻塔不同表面粗糙度時繞流特性的試驗研究[J].工業建筑,2006,36:301-304.Liu T C,Zhao L,Ding Z B.Test research of flow feature for hyperbolic circular section cooling tower with different superfacial roughness[J].Industrial Construction,2006,36:301-304.

[10]Liu X P,Zhao L,Ge Y J.Extreme value distribution of surface aerodynamic pressure of hyperbolic cooling tower[C].The 7thinternational colloquium on bluff body areodynamics and applications,Shanghai,2012.

[11]Zhang J F,Ge Y J,Zhao L.Influence of latitude wind pressure distribution on the responses of hyperbolodial cooling tower shell[J].Wind and Structures,2013,16(6):579-601.

[12]Zhao L,Ge Y J,Ahsan K.Fluctuating Wind Pressure Distribution around Full-scale Cooling Tower[J].Journal of Wind Engineering&Industrial Aerodynamics,2017,165:34-45.

[13]劉曉鵬.大型冷卻塔動態風壓現場實測及其效應綜合評價[D].上海:同濟大學,2013.Liu X P.On-spot measurement of stochastic wind pressure and its evaluation of multiple effects[D].Shanghai:Tongji University,2013.

[14]Zhao L,Ge Y J.Dynamic wind pressure distribution around full size cooling tower[C].Proceedings of the International Association for Shell and Spatial Structures(IASS)Symposium 2015,2015.

[15]Ke S T,Ge Y J.The influence of self-excited forces on wind loads and wind effects for super-large cooling towers[J].Journal of Wind Engineering and Industrial Aerodynamics,2014,132:125-135.

[16]Ke S T,Ge Y J,Zhao L,et al.A new methodology for analysis of equivalent static wind loads on super-large cooling towers[J].Journal of Wind Engineering and Industrial Aerodynamics,2012,111:30-39.

[17]許林汕,趙林,葛耀君.超大型冷卻塔隨機風振響應分析[J].振動與沖擊,2009,28(4):180-184.Xu L S,Zhao L,Ge Y J.Wind-excited stochastic responses of super large cooling towers[J].Journal of Vibration and Shock,2009,28(4):180-184.

[18]趙林,葛耀君,曹豐產.雙曲薄殼冷卻塔氣彈模型的等效梁格方法和實驗研究[J].振動工程學報,2008,21(1):31-37.Zhao L,Ge Y J,Cao F C.Equivalent beam-net design theory of aero-elastic model about hyperbolic thin-shell cooling towers and its experimental investigation[J].Journal of Vibration Engineering,2008,21(1):31-37.

[19]柯世堂,侯憲安,趙林,等.超大型冷卻塔風荷載和風振響應參數分析:自激力效應[J].土木工程學報,2012,45(12):45-53.Ke S T,Hou X A,Zhao L,et al.Parameter analysis of wind loads and wind induced responses forsuper-large cooling towers:self-excited force effect[J].China Civil Engineering Journal,2012,45(12):45-53.

[20]Zhao L,Ke S T,Ge Y J,et al.Wind-induced performance of super large cooling tower considering coupling effects between resonant and background components[C].International Symposium on Industrial Chimneys and Cooling Towers.Prague,2014.

[21]葛耀君,趙林.大型冷卻塔風荷載共振效應的耦合分析[C].全國結構振動與動力學學術研討會暨第四屆結構動力學專業委員會會議,中國 蘇州,2011.

[22]Ding Z B,Tamura Y,Yoshida A.Contributions to member stresses due to overall wind-induced behaviors of thin-walled cylindrical shell[J].Journal of Wind Engineering and Industrial Aerodynamics,2012,107:192-201.

[23]柯世堂,陳少林,趙林,等.超大型冷卻塔等效靜風荷載精細化計算及應用[J].振動測試與診斷,2013,33(5):824-830.Ke S T,Chen S L,Zhao L,et al.Fine calculation and application of equivalent static wind load of super large cooling tower[J].Journal of Vibration,Measurement&Diagnosis,2013,33(5):824-830.

[24]Ke S T,Ge Y J,Zhao L,et al.Wind-induced responses of super-large cooling towers[J].Journal of Central South University,2013,20(11):3216-3228.

[25]柯世堂.大型冷卻塔結構風效應和等效靜風荷載研究[D].上海:同濟大學,2011.Ke S T.Investigation of structrual wind effects and equivalent static wind loads of large cooling towers[D].Shanghai:Tongji University,2011.

[26]程霄翔.大型冷卻塔表面風壓及其雷諾數效應的風洞試驗和現場實測[D].上海:同濟大學,2015.Cheng X X.Study of wind pressures and their Reynold’s number effects on large cooling towers using wind tunnel tests and field measurements[D].Shanghai:Tongji University,2015.

[27]Zhao L,Chen X,Ke S T,et al.Aerodynamic and aero-elastic performances of super-large cooling towers[J].Wind and Structures,2014,19(4):443-465.

[28]程霄翔,趙林,葛耀君.典型矩形八塔超大型冷卻塔塔群風致干擾效應試驗[J].中南大學學報(自然科學版),2013,44(1):372-380.Cheng X X,Zhao L,Ge Y J.Wind tunnel investigation on interference effect of eight grouped super large cooling towers with rectangular arrangement[J].Journal of Central South U-niversity(Science and Technology),2013,44(1):372-380.

[29]張軍鋒,葛耀君,趙林.群塔布置對冷卻塔整體風荷載和風致響應的不同干擾效應[J].工程力學,2016,(8):15-23,44.Zhang J F,Ge Y J,Zhao L.Interference effects on global wind loads and wind induced responses for group hyperboloidal cooling towers[J].Engineering Mechanics,2016,(8):15-23,44.

[30]趙林,葛耀君,許林汕,等.超大型冷卻塔風致干擾效應試驗研究[J].工程力學,2009,(1):149-154.Zhao L,Ge Y J,Xu L S.et al.Wind tunnel investigation on wind-induced interference effects for super large cooling towers[J].Engineering Mechanics,2009,(1):149-154.

[31]張軍鋒,趙林,柯世堂,等.大型雙冷卻塔雙塔組合表面風壓干擾效應試驗研究[J].哈爾濱工業大學學報,2011,43(4):81-87.Zhang J F,Zhao L,Ke S T,et al.Wind tunnel investigation on wind-pressure interference effects for two large hyperbolodial cooling towers[J].Journal of Harbin Institute of Technology,2011,43(4):81-87.

[32]Zhao L,Zhan Y Y,Ge Y J.Unified wind-induced interference criteria based on the reinforcement envelope of cooling towers[C].International Symposium on Industrial Chimneys and Cooling Towers,Rotterdam,2016.

[33]張軍鋒.大型冷卻塔結構特性與風致干擾效應研究[D].上海:同濟大學,2012.Zhang J F.Structure behavior and wind-induced interference effects of large cooling towers[D].Shanghai:Tongji University,2012.

[34]Zhao L,Chen X,Ge Y J.Investigations of adverse wind loads on a large cooling tower for the six-tower combination[J].Applied Thermal Engineering,2016,105:988-999.

[35]張軍鋒,葛耀君,柯世堂,等.中美日三國規范高層結構風荷載標準值對比[J].湖南大學學報(自然科學版),2011,38(10):18-25.Zhang J F,Ge Y J,Ke S T,et al.Comparative study on the nominal value of wind loads on tall buildings among the codes of China,America and Japan[J].Journal of Hunan University(Natural Sciences),2011,38(10):18-25.

[36]張軍鋒,葛耀君,趙林.雙曲冷卻塔結構規范風荷載標準值對比[J].建筑科學與工程學報,2011,28(2):52-61.Zhang J F,Ge Y J,Zhao L.Comparisons of nominal values of wind loads on hyperboloidal cooling towers[J].Journal of Architecture and Civil Engineering,2011,28(2):52-61.

[37]張軍鋒,葛耀君,趙林.雙曲冷卻塔下部子午線形對結構性能的影響[J].建筑結構,2012,(2):118-121.Zhang J F,Ge Y J,Zhao L.Effect of bottom meridian curve on the mechanics behavior of hyperboloidal cooling tower[J].Building Structure,2012,(2):118-121.

[38]梁譽文.考慮多種風荷載分布模式的冷卻塔結構優化選型[D].上海:同濟大學,2016.Liang Y W.Structural optimization of large cooling towers considering various distribution patterns of wind loadings[D].Shanghai:Tongji University,2016.

[39]Zhao L,Liang Y W,Ge Y J,et al.Multi-object optimazation analysis of structural design for large cooling towers[J].Heat Transfer engineering,2017,38(11-12):1135-1145.

[40]趙林,葛耀君,項海帆.臺風風場隨機參數敏感性分析[J].同濟大學學報(自然科學版),2005,33(6):727-731.Zhao L,Ge Y J,Xiang H F.Stochastic parameter sensitivity analysis of typhoon wind field[J].Journal of Tongji University(Natural Science),2005,33(6):727-731.

[41]趙林,朱樂東,葛耀君.上海地區臺風風特性Monte-Carlo隨機模擬研究[J].空氣動力學學報,2009,(1):25-31.Zhao L,Zhu L D,Ge Y J.Monte-Carlo simulation about typhoon extreme value wind characteristics in Shanghai region[J].Acta Aerodynamica Sinica,2009,(1):25-31.

[42]趙林,葛耀君,宋麗莉,等.廣州地區臺風極值風特性蒙特卡羅隨機模擬[J].同濟大學學報(自然科學版),2007,08:1034-1038,1068.Zhao L,Ge Y J,Song L L,et al.Monte-Carlo simulation analysis of typhoon extreme value wind characteristics in Guangzhou[J].Journal of Tongji University(Natural Science),2007,08:1034-1038,1068.

[43]武占科,趙林,朱樂東,等.“羅莎”(0716)臺風高空實測脈動風特性分析[J].空氣動力學學報,2010,28(3):291-296.Wu Z K,Zhao L,Zhu L D,et al.High-altitude observation about turbulence characteristics for“Krosa”(0716)strong typhoon[J].Acta Aerodynamica Sinica,2010,28(3:)291-296.

[44]潘晶晶,趙林,冀春曉,等.東南沿海登陸臺風近地脈動特性分析[J].建筑結構學報,2016,(1):85-90.Pan J J,Zhao L,Ji C X,et al.Fluctuating wind characteristics near flat ground during typhoon landfall in southeast coastal areas of China[J].Journal of Building Structures,2016,(1):85-90.

[45]潘韜,趙林,曹曙陽,等.多風扇主動控制風洞類平板斷面抖振力識別研究[J].振動與沖擊,2010,29(6):178-183.Pan T,Zhao L,Cao S Y,et al.Buffeting force analysis of thin plate section in multiple fans active control wind tunnel[J].Journal of Vibration and Shock,2010,29(6):178-183.

[46]許坤,趙林,葛耀君,等.主動來流閉口箱梁氣動力時頻相關性[J].同濟大學學報自然科學版,2014,42(3):338-345.Xu K,Zhao L,Ge Y J,et al.Time-frequency coherences of aerodynamic loads on a streamlined bridge girder in an actively controlled wind tunnel[J].Journal of Tongji University(Natural Science),2014,42(3):338-345.

[47]潘韜,趙林,曹曙陽,等.主動來流條件類平板斷面氣動力荷載效應分析[J].實驗流體力學,2010,24(6):32-37.Pan T,Zhao L,Cao S Y,et al.Analysis of aerodynamic load effects on thin plat section under active control flow condition[J].Journal of Experiments in Fluid Mechanics,2010,24(6):32-37.

[48]趙楊,曹曙陽,Yukio Tamura,等.雷暴沖擊風模擬及其荷載的風洞試驗研究[J].振動與沖擊,2009,(4):1-3,9.Zhao Y,Cao S Y,Yukio T,et al.Simulation of downburst and its loads with wind tunnel test[J].Journal of Vibration and Shock,2009,(4):1-3,9.

[49]王錦,周強,曹曙陽,等.龍卷風風場的試驗模擬[J].同濟大學學報(自然科學版),2014,11:1654-1659.Wang J,Zhou Q,Cao S Y,et al.Physical study on tornadolike flow based on tornado vortex simulator[J].Journal of Tongji University(Natural Science),2014,11:1654-1659.

[50]鄭鴻飛,王小松,毛國輝.CAD二次開發在冷卻塔軟件開發中的應用[J].計算機技術與發展,2013,23(10):188-192.Zheng H F,Wang X S,Mao G H.Application of CAD secondary development in cooling tower software development[J].Computer Technology and Development,2013,23(10):188-192.

[51]Ruscheweyh H.Wind loadings on hyperbolic natural draft cooling towers[J].Journal of Wind Engineering and Industrial Aerodynamics,1976,1:335-340.

[52]Davenport A G,Isyumov N.The dynamic and static action of wind on hyperbolic cooling towers[R].Rresearch report No.BLWT1-66,Univ of Western Ontario,London,Ontario,Canada,1966.

[53]Sageau J F.Caractérisation des champs de pression moyens et fluctuantsàla surface des grands aérorefrigérants,Electricitéde France[C].Direction desétudes et Recherches,6 quai Watier,Chatou,France,1979.

[54]周良茂,李培華.兩個鄰近全尺寸雙曲型冷卻塔風壓分布的測量[J].氣動實驗與測量控制,1992,6(3):37-44.Zhou L M,Li P H.A measurement of the wind pressure distribution on two neighboring full-scale hyperbolic cooling towers[J].Aerodynamic Experiment and Measurement&Control,1992,6(3):37-44.

[55]趙林,劉曉鵬,高玲,等.大型冷卻塔表面脈動風壓原型實測與分布準則[J].土木工程學報,2017,50(1):1-11.Zhao L,Liu X P,Gao L,et al.Full scale measurement and distribution rules of surface fluctuating wind pressure of a large cooling tower[J].China Civil Engineering Journal,2017,50(1):1-11.

[56]柯世堂,趙林,葛耀君.大型雙曲冷卻塔表面脈動風壓隨機特性——風壓極值探討[J].實驗流體力學,2010,24(4):7-12.Ke S T,Zhao L,Ge Y J.Features of fluctuating wind pressure on large hyperbolic cooling tower:discussions on extreme wind pressure[J].Journal of Experiments in Fluid Mechanics,2010,24(4):7-12.

[57]Simiu E,Scanlan R H,Simiu E,et al.Wind effects on structures:funddementals and applications to design[M].New York:John Wiley,1996.

[58]Zhao L,Zhan Y Y,Liang Y W,et al.Comparison of wind-induced interference criteria and effects under typical towersgroup arrangements[C]//Spatial Structures in the 21st Century,Proceedings of the IASS Annual Symposium 2016,Tokyo,Japan,2016.

[59]Zhao L,Zhan Y Y,Ge Y J.Unified wind-induced interference criteria based on the reinforcement envelope of cooling towers[C].International Symposium on Industrial Chimneys and Cooling Towers,Rotterdam,Holland,2016.

[60]Nishi A,Kikugawa H,Matsuda Y,et al.Turbulence control in multiple-fan wind tunnels[J].Journal of Wind Engineering&Industrial Aerodynamics,1997,67-68(4):861-872.

[61]Wang J,Cao S,Pang W,et al.Experimental study on effects of ground roughness on flow characteristics of tornado-like vortices[J].Boundary-Layer Meteorology,2016:1-21.

[62]Cao S,Wang J,Cao J,et al.Experimental study of wind pressures acting on a cooling tower exposed to stationary tornadolike vortices[J].Journal of Wind Engineering&Industrial Aerodynamics,2015,145:75-86.

Comprehensive investigation into wind-induced aerodynamic and aeroelastic effects of large cooling towers

Zhao Lin,Zhan Yanyan*,Wang Zhinan,Liang Yuwen,Liu Xiaopeng,Cheng Xiaoxiang,Zhang Junfeng,Ke Shitang,Wang Xiaosong,Cao Shuyang,Ge Yaojun
(State Key Laboratory of Disaster Reduction in Civil Engineering,Tongji University,Shanghai 200092,China)

As typical thin-wall flexible structures,large cooling towers are of critical sensitivity to wind loads during its life-cycle.It's introduced systematically the research development of large cooling towers made by Tongji wind-resistance group during the past decades,including the study on statistic characteristics of the wind load from the full-scale observation,experiment and analysis of wind-induced structural responses,multi-objective structural optimization,the study on wind-induced loads and effects under disaster climates,the development of the self-developed whole process integration software platform,and so on.Some key problems in engineering application and fundamental investigation discovered during the development of large cooling towers in China have been resolved successfully.

large cooling towers;wind loads;wind-induced structural effects;multi-objective optimization;disaster climate;full-process analysis

TU973.2+13

:A

(編輯:李金勇)

2016-12-21;

:2017-03-13

國家自然科學基金項目(51678451和50978203);科技部重大科技項目(核電超大型冷卻塔結構研究及技術支持,2009ZX06004-010-HYJY);國家自然科學基金委優秀青年基金項目(51222809);教育部“新世紀優秀人才支持計劃”(NCET-13-0429)聯合資助項目

*通信作者E-mail:Zhan-yanyan@126.com

Zhao L,Zhan Y Y,Wang Z N,et al.Comprehensive investigation into wind-induced aerodynamic and aeroelastic effects of large cooling towers.Journal of Experiments in Fluid Mechanics,2017,31(3):1-15.趙 林,展艷艷,王志男,等.大型冷卻塔結構風致氣動和氣彈效應綜合研究.實驗流體力學,2017,31(3):1-15.

1672-9897(2017)03-0001-15

10.11729/syltlx20160201

趙 林(1974-),黑龍江牡丹江人,教授,博導。研究方向:橋梁和建筑結構風工程。通信地址:上海市四平路1239號同濟大學風工程館305室(200092)。E-mail:zhaolin@tongji.edu.cn