不同氣動(dòng)措施對(duì)特大型冷卻塔風(fēng)致響應(yīng)及穩(wěn)定性能影響分析*1

柯世堂， 杜凌云

(南京航空航天大學(xué) 土木工程系，江蘇 南京　210016)

?

不同氣動(dòng)措施對(duì)特大型冷卻塔風(fēng)致響應(yīng)及穩(wěn)定性能影響分析*1

柯世堂?， 杜凌云

(南京航空航天大學(xué) 土木工程系，江蘇 南京210016)

為研究不同氣動(dòng)措施對(duì)特大型冷卻塔結(jié)構(gòu)風(fēng)致強(qiáng)度及穩(wěn)定性能的影響，以?xún)?nèi)陸某核電特大型冷卻塔為例，對(duì)無(wú)氣動(dòng)措施和增設(shè)3種氣動(dòng)措施冷卻塔進(jìn)行剛體測(cè)壓風(fēng)洞試驗(yàn).基于試驗(yàn)結(jié)果對(duì)比分析了不同氣動(dòng)措施下冷卻塔表面平均和脈動(dòng)風(fēng)壓特性，然后采用有限元方法進(jìn)行不同氣動(dòng)措施下特大型冷卻塔的動(dòng)力特性、風(fēng)致響應(yīng)、局部和整體穩(wěn)定性能研究，最終提煉出不同氣動(dòng)措施對(duì)特大型冷卻塔結(jié)構(gòu)抗風(fēng)性能的影響規(guī)律.

特大型冷卻塔；氣動(dòng)措施；風(fēng)洞試驗(yàn)；風(fēng)壓特性；風(fēng)致響應(yīng)；穩(wěn)定性

隨著能源產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整，作為火/核發(fā)電廠重要構(gòu)筑物之一的冷卻塔的規(guī)模日趨高大化，國(guó)內(nèi)規(guī)范條款僅針對(duì)高度在165m以下的冷卻塔，其已無(wú)法滿(mǎn)足當(dāng)前特大型冷卻塔建設(shè)的需求，同時(shí)塔高和直徑的增大使特大型冷卻塔在風(fēng)荷載作用下的承載能力和穩(wěn)定性能成為制約其發(fā)展的瓶頸之一[1-2].

國(guó)內(nèi)外學(xué)者采用風(fēng)洞試驗(yàn)和有限元方法對(duì)大型冷卻塔的風(fēng)致穩(wěn)定性能進(jìn)行了大量研究[3-6]，為其抗風(fēng)設(shè)計(jì)提供了很好地技術(shù)支持.然而，國(guó)內(nèi)外已有研究成果均未涉及不同氣動(dòng)措施[7-8]下特大型冷卻塔的風(fēng)壓分布特性，更缺乏不同氣動(dòng)措施對(duì)其風(fēng)致響應(yīng)和穩(wěn)定性影響的定性和定量分析.

鑒于此，本文以?xún)?nèi)陸某核電特大型冷卻塔工程為背景，對(duì)無(wú)氣動(dòng)措施和增設(shè)3種不同氣動(dòng)措施的冷卻塔進(jìn)行剛體模型測(cè)壓風(fēng)洞試驗(yàn)，基于試驗(yàn)結(jié)果對(duì)比分析表面平均和脈動(dòng)風(fēng)壓特性；再采用有限元方法進(jìn)行不同氣動(dòng)措施下特大型冷卻塔的風(fēng)致響應(yīng)及整體和局部穩(wěn)定性研究，并與無(wú)氣動(dòng)措施下冷卻塔的計(jì)算結(jié)果進(jìn)行對(duì)比分析，最終提煉出不同氣動(dòng)措施對(duì)特大型冷卻塔結(jié)構(gòu)抗風(fēng)設(shè)計(jì)的影響規(guī)律，主要結(jié)論可為此類(lèi)特大型冷卻塔氣動(dòng)措施的選取提供依據(jù).

1　風(fēng)洞試驗(yàn)及結(jié)果分析

1.1剛體測(cè)壓風(fēng)洞試驗(yàn)

本試驗(yàn)結(jié)構(gòu)原型采用江西某核電特大型自然通風(fēng)冷卻塔[1]，塔高215m，淋水面積18 300m2，喉部高度160m，中面半徑49.64m，塔筒分段等厚，最小厚度0.26m，最大厚度1.8m.表1給出了該工程冷卻塔的主要結(jié)構(gòu)尺寸.

測(cè)壓試驗(yàn)所用風(fēng)洞為全鋼結(jié)構(gòu)閉口回流式低速大氣邊界層風(fēng)洞，剛體模型采用1∶500縮尺比，沿環(huán)向和子午向共布置36×12個(gè)測(cè)點(diǎn).同時(shí)在來(lái)流前部放置三角尖劈和地面粗糙元來(lái)模擬B類(lèi)地貌的大氣邊界層風(fēng)場(chǎng).由于物理風(fēng)洞本身的局限性，可通過(guò)適當(dāng)改變模型表面粗糙度和調(diào)整試驗(yàn)風(fēng)速來(lái)近似模擬冷卻塔高雷諾數(shù)時(shí)的繞流特性[9-10].通過(guò)比較確定采用二三層紙帶間隔分布的形式沿圓周均勻布置寬5mm，厚0.1mm共計(jì)36條豎向通長(zhǎng)粗糙紙帶和來(lái)流風(fēng)速10m/s手段來(lái)模擬雷諾數(shù)效應(yīng)(模型雷諾數(shù)為1.51×105).圖1給出了在10m/s試驗(yàn)風(fēng)速下冷卻塔中間斷面平均表面壓力系數(shù)分布與規(guī)范值[11]的對(duì)比曲線，由圖比較可知二者吻合較好，故后續(xù)不同氣動(dòng)措施下的冷卻塔模型均采用此雷諾數(shù)模擬方法.

表1　特大型冷卻塔主要結(jié)構(gòu)尺寸Tab．1　The size characteristics of main componentfor super large cooling tower

3種氣動(dòng)措施分別為：在進(jìn)風(fēng)口上部設(shè)置外部進(jìn)水槽、矩形導(dǎo)風(fēng)板和弧形導(dǎo)風(fēng)板，每種導(dǎo)風(fēng)裝置的尺寸如圖2所示.相應(yīng)計(jì)算模型簡(jiǎn)稱(chēng)無(wú)措施、措施1、措施2和措施3.其中不同氣動(dòng)措施模型如圖3所示.

角度/(°)圖1　試驗(yàn)與規(guī)范壓力分布比較曲線圖Fig.1　Specifications for test results and comparisonof surface pressure distribution

圖2　不同氣動(dòng)措施下模型尺寸示意圖(mm)Fig.2　Detail sizes of four conditions with different aerodynamic measures(mm)

圖3　不同氣動(dòng)措施下冷卻塔剛體測(cè)壓模型示意圖Fig.3　The sketches of cooling towers with different aerodynamic measures

1.2結(jié)果分析

圖4給出了子午向0～50m，50～100m，100～150m和150～215m高度區(qū)間內(nèi)4種氣動(dòng)措施下冷卻塔平均風(fēng)壓隨環(huán)向角度變化曲線.由圖可知，隨著高度區(qū)間的增大，不同氣動(dòng)措施下冷卻塔外表面的平均風(fēng)壓系數(shù)與無(wú)氣動(dòng)措施的分布差別越來(lái)越小，最大影響量從27.69%降至18.02%，在150～215m高度區(qū)域內(nèi)基本與無(wú)氣動(dòng)措施風(fēng)壓曲線分布一致，局部點(diǎn)差異較大；不同氣動(dòng)措施對(duì)120°～240°角度范圍內(nèi)即背風(fēng)區(qū)域的風(fēng)壓系數(shù)影響較大，對(duì)側(cè)風(fēng)區(qū)域和迎風(fēng)區(qū)域的平均風(fēng)壓影響相對(duì)較小.

環(huán)向角度/(°) (a)0～50 m高度

環(huán)向角度/(°) (b)50～100 m高度

環(huán)向角度/(°) (c)100～150 m高度

環(huán)向角度/(°) (d)150～215 m高度圖4　4種冷卻塔不同高度區(qū)域平均風(fēng)壓對(duì)比曲線Fig.4　Average wind pressure coefficientfor four types of cooling tower

圖5給出了子午向0～50m，50～100m，100～150m和150～215m高度區(qū)間內(nèi)不同氣動(dòng)措施下冷卻脈動(dòng)風(fēng)壓系數(shù)均方根對(duì)比曲線，對(duì)比可知在0～50m范圍內(nèi)，4種冷卻塔脈動(dòng)風(fēng)壓系數(shù)均方根差別較大，最大百分比為33.87%，但隨著塔高的增大不同氣動(dòng)措施下的脈動(dòng)風(fēng)壓均方根逐漸接近無(wú)氣動(dòng)措施的分布曲線，且沿環(huán)向角度的變化規(guī)律趨于一致；在120°～240°角度范圍內(nèi)即背風(fēng)區(qū)域不同氣動(dòng)措施對(duì)脈動(dòng)風(fēng)壓的影響最為顯著，在其他角度范圍內(nèi)區(qū)別相對(duì)較小.

環(huán)向角度/(°) (a)0～50 m高度

環(huán)向角度/(°) (b)50～100 m高度

環(huán)向角度/(°) (c)100～150 m高度

環(huán)向角度/(°) (d)150～215 m高度圖5　4種冷卻塔不同高度區(qū)域脈動(dòng)風(fēng)壓系數(shù)均方根對(duì)比曲線Fig.5　 RMS fluctuating wind pressure coefficient for four types of cooling tower

2　動(dòng)力特性分析

采用有限元方法分別對(duì)4種氣動(dòng)措施下的特大型冷卻塔進(jìn)行動(dòng)力特性分析，圖6給出了對(duì)應(yīng)的有限元整體及局部模型，塔筒采用Shell63單元，子午向和環(huán)向分別劃分為135和160個(gè)單元，支柱和環(huán)基采用Beam188單元，環(huán)基下部采用空間彈簧單元模擬彈性地基，每根樁基均采用3個(gè)力和力矩彈簧單元分別模擬樁沿豎向、環(huán)向、徑向、繞豎向、繞環(huán)向和繞徑向的作用，彈簧單元一端與環(huán)基剛性連接，另一端與地面固接約束，環(huán)基與支柱下部剛接，支柱上部與塔筒下部節(jié)點(diǎn)耦合.

圖6　 4種冷卻塔的整體和局部有限元模型Fig.6　Global and local finite element model of four types of cooling tower

圖7給出了4種塔型前100階頻率對(duì)比曲線.圖8給出了4種冷卻塔第一階振型對(duì)比圖．由圖可看出：無(wú)氣動(dòng)措施冷卻塔首階振型為環(huán)向諧波的有3個(gè)，豎向諧波的有2個(gè)，而3種氣動(dòng)措施下冷卻塔第一階振型特性完全一致，環(huán)向諧波有4個(gè)，豎向諧波2個(gè)且下部豎向諧波形狀完整.由圖可看出，四者基頻分別為：氣動(dòng)措施2(0.657 0)> 氣動(dòng)措施1(0.656 6)> 氣動(dòng)措施3(0.656 2)> 無(wú)氣動(dòng)措施(0.643 8)；措施2冷卻塔的各階頻率相比其他三者為最大，在0～60階范圍內(nèi)，措施1和措施3兩種冷卻塔的頻率十分接近.

振型階數(shù)圖7　4種冷卻塔前100階頻率對(duì)比圖Fig.7　The first 100 order frequency comparisonchart of four types of cooling tower

圖8　第一階振型對(duì)比圖Fig.8　The first order modal contrast figure

3　同組風(fēng)壓下冷卻塔的受力性能分析

本節(jié)均采用無(wú)氣動(dòng)措施冷卻塔外表面的平均風(fēng)壓系數(shù)對(duì)4種冷卻塔模型進(jìn)行加載分析.

3.1環(huán)基與支柱響應(yīng)

圖9和圖10給出了同一組風(fēng)荷載作用下不同氣動(dòng)措施冷卻塔環(huán)基徑向位移曲線和支柱軸力曲線，由圖分析可知：1)氣動(dòng)措施的設(shè)置對(duì)于環(huán)基側(cè)風(fēng)區(qū)域的徑向位移影響最大，對(duì)背風(fēng)區(qū)域影響相對(duì)略小，對(duì)迎風(fēng)區(qū)域的環(huán)基徑向位移影響最小，4者幾乎吻合；2)不同氣動(dòng)措施下冷卻塔的支柱軸力分布趨勢(shì)一致，在側(cè)風(fēng)區(qū)和背風(fēng)區(qū)局部支柱處略有差別.

環(huán)向角度/(°)圖9　同組風(fēng)荷載作用下冷卻塔環(huán)基徑向位移曲線Fig.9　The radial displacement curve of ring basecooling tower under the same wind load

支柱編號(hào)圖10　 同組風(fēng)荷載作用下冷卻塔支柱頂部軸向力曲線Fig.10　The axial force curve of pillar top of coolingtower under the same wind load

3.2塔筒響應(yīng)

圖11給出了同一組風(fēng)壓下塔筒70°子午線上的節(jié)點(diǎn)徑向位移和單元環(huán)向應(yīng)力沿高度變化曲線圖.由圖可見(jiàn)當(dāng)對(duì)不同氣動(dòng)措施冷卻塔施加同一組風(fēng)壓時(shí)，70°子午線上的節(jié)點(diǎn)徑向位移和單元環(huán)向應(yīng)力幾乎沒(méi)有差別，僅在氣動(dòng)措施所在處高度及喉部有較小差別.

徑向位移/m(a)節(jié)點(diǎn)徑向位移

環(huán)向應(yīng)力/Pa(b)單元環(huán)向應(yīng)力圖11　同組風(fēng)荷載作用下冷卻塔70°子午線上節(jié)點(diǎn)徑向位移和單元環(huán)向應(yīng)力曲線Fig.11　The node radial displacement and element hoop stress curves of 70° meridian ring under the same wind load

3.3整體穩(wěn)定性驗(yàn)算

進(jìn)行整體穩(wěn)定性[12]驗(yàn)算的輸入荷載組合為自重+K(風(fēng)荷載+內(nèi)吸力)，K為失穩(wěn)特征值，失穩(wěn)臨界風(fēng)速是K與基本風(fēng)速的乘積，此時(shí)的風(fēng)荷載均為無(wú)氣動(dòng)措施下冷卻塔的表面風(fēng)荷載.

計(jì)算得到4種氣動(dòng)措施下冷卻塔的屈曲系數(shù)、臨界風(fēng)速及屈曲模態(tài)如表2所示.由表可知:1)氣動(dòng)措施的設(shè)置可以提高冷卻塔的靜風(fēng)整體穩(wěn)定性；2)同一風(fēng)荷載作用下措施1對(duì)提高冷卻塔整體穩(wěn)定性的影響效果最好，此時(shí)對(duì)應(yīng)的屈曲失穩(wěn)臨界風(fēng)速為217.29m/s.

表2　同組風(fēng)荷載作用下冷卻塔屈曲模態(tài)Tab.2　Buckling mode for four types of cooling tower under the same wind load

4　風(fēng)致響應(yīng)特性

本節(jié)采用不同氣動(dòng)措施下對(duì)應(yīng)各自的冷卻塔外表面平均風(fēng)壓系數(shù)對(duì)四種冷卻塔進(jìn)行靜風(fēng)加載，具體研究不同氣動(dòng)措施及其風(fēng)壓分布對(duì)冷卻塔受力性能和屈曲穩(wěn)定的影響.

4.1環(huán)基位移

圖12給出了不同氣動(dòng)措施下冷卻塔環(huán)基的徑向、環(huán)向和豎向位移曲線．由圖可見(jiàn)：1)對(duì)應(yīng)風(fēng)荷載作用不同氣動(dòng)措施對(duì)冷卻塔環(huán)基的變形影響作用較大，4種塔型位移變化規(guī)律一致； 2)不同氣動(dòng)措施對(duì)冷卻塔的位移影響在側(cè)風(fēng)區(qū)域和背風(fēng)區(qū)域影響較大，在迎風(fēng)區(qū)域影響較?。?)冷卻塔的豎向位移變化劇烈，不同范圍內(nèi)的節(jié)點(diǎn)豎向位移突變嚴(yán)重.

環(huán)向角度/(°) (a)徑向位移

環(huán)向角度/(°) (b)環(huán)向位移

環(huán)向角度/(°) (c)豎向位移圖12　4種冷卻塔環(huán)基位移曲線Fig.12　The displacement curve of ringbase for four types of cooling tower

4.2支柱內(nèi)力

圖13給出了不同氣動(dòng)措施下冷卻塔支柱頂部軸力變化曲線，按支柱傾斜方向分為奇數(shù)支柱和偶數(shù)支柱.由圖看出：1)奇數(shù)支柱與偶數(shù)支柱軸向力呈軸對(duì)稱(chēng)；2)不同氣動(dòng)措施下冷卻塔的支柱軸向力分布趨勢(shì)幾乎相同，無(wú)氣動(dòng)措施奇數(shù)和偶數(shù)支柱軸向力分別在支柱編號(hào)16～24范圍內(nèi)和編號(hào)24～32范圍內(nèi)突然減小，其它范圍內(nèi)冷卻塔支柱軸向力以氣動(dòng)措施1作用下最大，以無(wú)氣動(dòng)措施作用下為最小，氣動(dòng)措施2和3對(duì)支柱軸向力影響相當(dāng).

支柱編號(hào)圖13　4種冷卻塔支柱頂部軸向力曲線Fig.13　 The axial force curve of pillar topof four types of cooling tower

4.3塔筒位移

選擇迎風(fēng)點(diǎn)(0°)、零壓力系數(shù)點(diǎn)(30°)、負(fù)壓極大值點(diǎn)(70°)及背風(fēng)點(diǎn)(180°)4個(gè)代表性區(qū)域進(jìn)行不同氣動(dòng)措施冷卻塔的筒壁位移響應(yīng)分析.圖14給出了不同氣動(dòng)措施冷卻塔在各自風(fēng)荷載作用下的塔筒0°,30°,70°及180°子午線上徑向位移隨高度的變化曲線.對(duì)比分析可得：1)不同氣動(dòng)措施冷卻塔在0°和70°子午線上節(jié)點(diǎn)徑向位移在喉部以下比較接近，在喉部以上數(shù)值稍有差異，最大相差12.67%；2)在30°子午線上的徑向位移差異較大，125m以下無(wú)氣動(dòng)措施冷卻塔位移最大，帶弧形導(dǎo)風(fēng)板冷卻塔位移最小，達(dá)到喉部高度后位移突然減小，其中以無(wú)氣動(dòng)措施冷卻塔減小趨勢(shì)最明顯；3)180°子午線上節(jié)點(diǎn)位移在塔筒中下部以無(wú)氣動(dòng)措施冷卻塔最大，達(dá)到喉部高度后位移均開(kāi)始減小.

喉部壁厚較薄，屬于冷卻塔的薄弱部位，有必要對(duì)其徑向位移分布特性進(jìn)行研究.圖15給出了4種氣動(dòng)措施下冷卻塔的喉部徑向位移隨角度的變化曲線，可將0.00圓環(huán)假定為冷卻塔喉部原形.由圖可見(jiàn)4種冷卻塔的喉部徑向位移大小和變化趨勢(shì)幾乎一致，其中喉部最大徑向負(fù)位移-0.041出現(xiàn)在正迎風(fēng)角0°處，最大正位移0.036出現(xiàn)在±70°附近；在0°～45°范圍內(nèi)，徑向位移為負(fù)，且逐漸減小；45°～70°范圍內(nèi)，徑向位移為正，且逐漸增大；70°～100°范圍內(nèi)，徑向位移為正，且逐漸減小；100°～180°范圍內(nèi)，徑向位移先增大后減小至0繼而增大至0.005左右.

徑向位移/m(a)0°子午線

徑向位移/m(b)30°子午線

徑向位移/m(c)70°子午線

徑向位移/m(d)180°子午線圖14　塔筒典型子午線徑向位移變化曲線Fig.14　The radial displacement curves of typicalmeridian ring for tower drum

圖15　4種冷卻塔喉部徑向位移曲線Fig.15　 The radial displacement curvesof throat pf four types of cooling tower

4.4塔筒應(yīng)力

以無(wú)氣動(dòng)措施冷卻塔為例，圖16提取了塔筒所有節(jié)點(diǎn)的環(huán)向和子午向應(yīng)力等勢(shì)線圖．由圖可發(fā)現(xiàn)環(huán)向和子午向應(yīng)力較大值均出現(xiàn)在±70°左右，故后續(xù)以70°子午線上的應(yīng)力值進(jìn)行4種氣動(dòng)措施下冷卻塔的應(yīng)力對(duì)比研究.

環(huán)向角度/(°) (a)環(huán)向應(yīng)力

環(huán)向角度/(°) (b)子午向應(yīng)力圖16　無(wú)氣動(dòng)措施冷卻塔應(yīng)力等勢(shì)線圖Fig.16　The stress to the equipotential chart ofcooling tower without aerodynamic measure

圖17給出了70°子午向4種氣動(dòng)措施下冷卻塔的環(huán)向和子午向應(yīng)力沿塔高的變化曲線對(duì)比圖.由圖對(duì)比分析可見(jiàn)：1)4種冷卻塔的環(huán)向、子午向應(yīng)力沿塔高變化基本相同；2)由于自重的積累，塔筒底部應(yīng)力最大，在氣動(dòng)措施設(shè)置高度處應(yīng)力驟減；3)隨著冷卻塔標(biāo)高的增大，應(yīng)力逐漸減小，在喉部處略有突變，但總體趨勢(shì)不變，冷卻塔頂部壁厚較薄，平臺(tái)外伸，應(yīng)力突增.

環(huán)向角度/(°) (a)環(huán)向應(yīng)力

環(huán)向角度/(°) (b)子午向應(yīng)力圖17　4種冷卻塔70°子午線環(huán)向、子午向應(yīng)力曲線Fig.17　The circular and meridian stress curve of 70° meridian ring of four types of cooling tower

5　靜風(fēng)穩(wěn)定性分析

5.1計(jì)算特征值法整體穩(wěn)定性驗(yàn)算

采用不同氣動(dòng)措施下對(duì)應(yīng)各自的冷卻塔外表面風(fēng)荷載進(jìn)行穩(wěn)定性驗(yàn)算，計(jì)算得到4種氣動(dòng)措施下冷卻塔的屈曲系數(shù)、臨界風(fēng)速及屈曲模態(tài)如表3所示.由表可知：1)氣動(dòng)措施的設(shè)置可以提高冷卻塔的靜風(fēng)整體穩(wěn)定性；2)措施2和3對(duì)冷卻塔整體穩(wěn)定性的影響效果相當(dāng)，氣動(dòng)措施1對(duì)冷卻塔的整體穩(wěn)定性改善效果最為顯著，對(duì)應(yīng)的屈曲失穩(wěn)臨界風(fēng)速為71.45m/s，小于統(tǒng)一風(fēng)壓分布模式下的臨界風(fēng)速71.7 0m/s.

表3　冷卻塔屈曲模態(tài)Tab．3　Buckling mode for four types of cooling tower

5.2局部穩(wěn)定性驗(yàn)算

采用規(guī)范[13-14]提出的屈曲應(yīng)力計(jì)算方法分析冷卻塔的局部穩(wěn)定性.圖18給出了4種氣動(dòng)措施下塔筒局部穩(wěn)定系數(shù)云圖，由于計(jì)算得到的不同部位局部穩(wěn)定因子[15]數(shù)值相差較大，為便于對(duì)比，圖中數(shù)值由實(shí)際局部穩(wěn)定因子取對(duì)數(shù)給出.由圖可見(jiàn)：1)4種氣動(dòng)措施下冷卻塔的最小局部穩(wěn)定系數(shù)為：無(wú)氣動(dòng)措施(4.34)<氣動(dòng)措施3(4.46)<氣動(dòng)措施2(4.50)<氣動(dòng)措施1(4.56)，出現(xiàn)位置大致相同，均在環(huán)向±75°附近、高度30～80m范圍內(nèi)，故在設(shè)計(jì)時(shí)應(yīng)局部加固此處；2)無(wú)氣動(dòng)措施冷卻塔的局部穩(wěn)定因子最小，氣動(dòng)措施2下冷卻塔的數(shù)值分布范圍比其他三者廣；3)氣動(dòng)措施1與氣動(dòng)措施3下冷卻塔的數(shù)值分布相似，但后者偏于安全，故在局部穩(wěn)定性方面，建議選擇弧形導(dǎo)風(fēng)裝置進(jìn)行氣動(dòng)措施設(shè)置.

環(huán)向角度/(°) (a)無(wú)氣動(dòng)措施

環(huán)向角度/(°) (b)氣動(dòng)措施1

環(huán)向角度/(°) (c)氣動(dòng)措施2

環(huán)向角度/(°) (d)氣動(dòng)措施3圖18　 塔筒單元局部穩(wěn)定安全系數(shù)Fig.18　Local stability safety coefficientof all tower drum elements

6　結(jié)　論

結(jié)合風(fēng)洞試驗(yàn)和有限元分析，對(duì)不同氣動(dòng)措施下冷卻塔表面風(fēng)壓特性和風(fēng)致強(qiáng)度及穩(wěn)定性進(jìn)行了對(duì)比研究，得到主要結(jié)論如下：

1)不同氣動(dòng)措施對(duì)于大型冷卻塔表面平均和脈動(dòng)風(fēng)壓分布有一定影響，主要體現(xiàn)在0～50m高度區(qū)域內(nèi)，最大影響量分別為27.69%和33.87%.隨著塔高增加不同氣動(dòng)措施冷卻塔表面環(huán)向平均和脈動(dòng)風(fēng)壓分布逐漸與無(wú)氣動(dòng)措施的分布曲線接近，在150～215m高度區(qū)域內(nèi)風(fēng)壓分布基本一致；不同氣動(dòng)措施對(duì)環(huán)向120°～240°角度范圍內(nèi)即背風(fēng)區(qū)域的風(fēng)荷載影響較大，對(duì)側(cè)風(fēng)區(qū)域和迎風(fēng)區(qū)域影響相對(duì)較??；

2)針對(duì)不同氣動(dòng)措施對(duì)冷卻塔風(fēng)致位移和內(nèi)力響應(yīng)的影響，采用統(tǒng)一風(fēng)壓分布模式計(jì)算要比采用對(duì)應(yīng)模型表面風(fēng)壓要小，而相應(yīng)的屈曲失穩(wěn)臨界風(fēng)速前者則要偏大為71.70m/s；

3)不同氣動(dòng)措施均可提高冷卻塔結(jié)構(gòu)基頻,分別為：氣動(dòng)措施2(0.657 0)> 氣動(dòng)措施1(0.656 6)> 氣動(dòng)措施3(0.656 2)> 無(wú)氣動(dòng)措施(0.643 8)，且對(duì)冷卻塔的振型影響較為顯著；

4)不同氣動(dòng)措施對(duì)環(huán)基位移、筒壁位移和應(yīng)力分布特性影響較小，但對(duì)支柱軸向力影響較大，帶有外部進(jìn)水槽的冷卻塔軸向力最大，兩種導(dǎo)風(fēng)板對(duì)支柱軸向力的影響相當(dāng)；

5)不同氣動(dòng)措施均可提高冷卻塔的靜風(fēng)整體穩(wěn)定性，其中以外部進(jìn)水槽對(duì)冷卻塔的整體穩(wěn)定性改善效果最優(yōu)，對(duì)應(yīng)屈曲失穩(wěn)臨界風(fēng)速為71.45m/s；增設(shè)外部進(jìn)水槽冷卻塔局部穩(wěn)定因子值較大，在局部穩(wěn)定性方面相對(duì)其他氣動(dòng)措施更為安全.

[1]柯世堂, 侯憲安, 趙林,等. 特大型冷卻塔風(fēng)荷載和風(fēng)振響應(yīng)參數(shù)分析：自激力效應(yīng)[J]. 土木工程學(xué)報(bào), 2012, 45(12): 45-53.

KEShi-tang,HOUXian-an,ZHAOLin, et al.Parameteranalysisofwindloadsandwindinducedresponsesforsuper-largecoolingtowers:self-excitedforceeffects[J].ChinaCivilEngineeringJournal, 2012, 45(12): 45-53. (InChinese)

[2]趙林, 李鵬飛, 葛耀君. 等效靜風(fēng)荷載下特大型冷卻塔受力性能分析[J]. 工程力學(xué), 2008, 25(7): 79-86.

ZHAOLin,LIPeng-fei,GEYao-jun.Numericalinvestigationonperformanceofsuperlargecoolingtowersunderequivalentstaticwindload[J].EngineerngMechanics, 2008, 25(7): 79-86. (InChinese)

[3]KEST,LIANGJ,ZHAOL, et al .InfluenceofventilationrateontheaerodynamicinterferencefortwoIDCTsbyCFD[J].WindandStructures,AnInternationalJournal, 2015, 20(3): 449-468.

[4]柯世堂, 趙林, 葛耀君. 大型雙曲冷卻塔表面脈動(dòng)風(fēng)壓隨機(jī)特性——風(fēng)壓極值探討[J]. 實(shí)驗(yàn)流體力學(xué), 2010, 24(4): 7-12.

KEShi-tang,ZHAOLin,GEYao-jun.Featuresoffluctuatingwindpressureonlargehyperboliccoolingtower:discussionsonextremewindpressure[J].JournalofExperimentsinFluidMechanics, 2010, 24(4): 7-12. (InChinese)

[5]鄒云峰, 牛華偉, 陳政清. 特大型冷卻塔單塔外表面風(fēng)荷載三維效應(yīng)及其設(shè)計(jì)取值[J]. 振動(dòng)與沖擊, 2013, 32(24): 76-82.

ZOUYun-feng,NIUHua-wei,CHENZheng-qing.Three-dimensionaleffectanddesignvaluesofoutersurfacewindloadingforasinglesuper-largecoolingtower[J].JournalofVibrationandShock, 2013, 32(24): 76-82. (InChinese)

[6]沈國(guó)輝, 余關(guān)鵬, 孫炳楠,等. 大型冷卻塔雙塔干擾的風(fēng)洞試驗(yàn)研究[J]. 振動(dòng)與沖擊, 2011, 30(3): 110-114.

SHENGuo-hui,YUGuan-peng,SUNBing-nan, et al.Studyoninterferenceeffectsoftwolargehyperboliccoolingtowersbyusingwindtunneltest[J].JournalofVibrationandShock, 2011, 30(3): 110-114. (InChinese)

[7]SPALDINGB,SINGHAMR.Dieleistungvonkaminkuhlern[J].VergleichTheorieandPraxisChamTechnik, 1996, 18(7): 385-391.

[8]艾·漢佩. 冷卻塔[M]. 胡賢,譯. 北京: 電力工業(yè)出版社, 1978:118-119.

HANMPEE.Coolingtower[M] ．TranslatedbyHUXian．Translate.Beijing:ElectricPowerIndustialPress,1978:118-119. (InChinese)

[9]PIMERM.Windpressurefluctuationsonacoolingtower[J].JournalofWindEngineeringandIndustrialAerodynamics, 1982,10: 343-360.

[10]董國(guó)朝, 陳政清, 羅建輝, 等. 冷卻塔混凝土粗糙度對(duì)平均風(fēng)壓系數(shù)的影響[J]. 湖南大學(xué)學(xué)報(bào):自然科學(xué)版, 2011, 38(7): 6-12.

DONGGuo-chao,CHENZheng-qing,LUOJian-hui, et al.Effectofsurfaceroughnessonthemeanpressurecoefficientofconcretecoolingtower[J].JournalofHunanUniversity:NaturalSciences,2011, 38(7): 6-12. (InChinese)

[11]GB/T50102-2003工業(yè)循環(huán)水冷卻設(shè)計(jì)規(guī)范[S]. 北京: 中國(guó)計(jì)劃出版社, 2003:17-18.

GB/T50102-2003Codefordesignofcoolingforindustrialrecirculatingwater[S].Beijing:ChinaSharpePress, 2003:17-18. (InChinese)

[12]張軍鋒, 葛耀君, 趙林. 基于風(fēng)洞試驗(yàn)的雙曲冷卻塔靜風(fēng)整體穩(wěn)定研究[J]. 工程力學(xué), 2012, 29(5): 68-77.

ZHANGJun-feng,GEYao-jun,ZHAOLin.Studyonglobalaerostaticstabilityofhyperboloidalcoolingtowersbasedonthewindtunneltests[J].EngineerngMechanics, 2012, 29(5): 68-77. (InChinese)

[13]PAULAYT.Couplingbeamsofreinforcedconcreteshearwalls[J].JournaloftheStructuralDivision, 1971, 97(ST3): 843-862.

[14]DL/T5339-2006火力發(fā)電廠水工設(shè)計(jì)規(guī)范[S]. 北京: 中國(guó)電力出版社, 2006:115-116.

DL/T5339-2006Codeforhydraulicdesignoffossilfuelpowerplants[S].Beijing:ChinaPowerPress, 2006: 115-116. (InChinese)

[15]KEST,GEYJ,ZHAOL, et al.Anewmethodologyforanalysisofequivalentstaticwindloadsonsuper-largecoolingtowers[J].JournalofWindEngineeringandIndustrialAerodynamics, 2012，111(3): 30-39.

Impact Study on Wind-induced Response and Stability for Super Large Cooling Tower with Different Aerodynamic Measures

KE Shi-tang?, DU Ling-yun

(Dept of Civil Engineering, Nanjing Univ of Aeronautics and Astronautics, Nanjing, Jiangsu210016, China)

Tostudythewind-inducedstrengthandstabilitypropertiesoflargecoolingtowerswithdifferentaerodynamicmeasures,thepressuresoftherigid-bodymodelswithoutorwiththreedifferentaerodynamicmeasuresweremeasuredbythewindtunneltests.Anuclearsuper-largecoolingtowerinlandwasconsideredasthetestspecimen.Thesurfacewindmeanandfluctuatingpressurecharacteristicswereinvestigatedfromthetestresults.Furthermore,finiteelementanalysiswasconductedtoevaluatethedynamiccharacteristicsofnaturalvibration,thewind-inducedresponse,andtheoverallandlocalstability.Finally,theeffectiverulesofthedifferentaerodynamicmeasuresonwindresistanceforsuper-largecoolingtowerwereproposed.

superlargecoolingtowers;aerodynamicmeasures;windtunneltest;windpressurecharacteristic;wind-inducedresponse;stability

1674-2974(2016)05-0079-11

2015-07-04

國(guó)家自然科學(xué)基金資助項(xiàng)目(51208254)，National Natural Science Foundation of China(51208254)；博士后科學(xué)基金資助項(xiàng)目(2013M530255, 1202006B)作者簡(jiǎn)介：柯世堂(1982-)，男，安徽池州人，南京航空航天大學(xué)副教授

TU279.741

A

?通訊聯(lián)系人，E-mail:keshitang@163.com