游溢 晏致濤 陳俊帆 李文勝 何成 閆渤文

摘 要：輸電塔結(jié)構(gòu)具有高柔、輕質(zhì)和小阻尼的特性，因此風荷載是輸電塔結(jié)構(gòu)的主要控制荷載.以蘇通大跨越輸電塔結(jié)構(gòu)作為工程背景，開展了一系列的剛性節(jié)段模型風洞試驗，研究了在光滑均勻流作用下的圓鋼管格構(gòu)式塔架氣動力特性.并基于計算流體力學（Computational Fluid Dynamics，簡稱CFD）技術(shù)，采用大渦模擬（Large-eddy simulation，簡稱LES）方法分別對圓鋼管格構(gòu)式塔架節(jié)段的順風向、橫風向和扭轉(zhuǎn)向風荷載特性進行計算，并將計算結(jié)果與實驗值以及各國規(guī)范進行了對比.此外，討論了湍流度分別為5%、10%、15%和20%四種情況下的三分力系數(shù).結(jié)果表明，本文的數(shù)值模擬結(jié)果與各國規(guī)范及風洞試驗結(jié)果較為一致，且湍流度的變化對結(jié)構(gòu)三分力系數(shù)的脈動成分有顯著影響.此外，從節(jié)段模型的表面風壓力分布情況來看，最大正壓出現(xiàn)在迎風面的主管正面和輔材相交節(jié)點的正面，最大負壓出現(xiàn)在主管的側(cè)面.通過渦量圖對節(jié)段模型流場結(jié)構(gòu)分析發(fā)現(xiàn)，x向的渦量離散程度更高，而y和z向上渦量分布更加均勻連續(xù).本文的研究成果對實際工程中的圓管截面格構(gòu)式結(jié)構(gòu)風荷載評估具有重要的參考價值.

關(guān)鍵詞：格構(gòu)式塔架；大渦模擬；氣動力；輸電塔

中圖分類號：TM75，TU312.1 文獻標志碼：A

文章編號：1674—2974（2018）07—0054—07

Abstract： Due to the high flexibility， light self-weight and low damping ratio of the power transmission tower， wind loads are the dominant lateral loads at its design stage. This study conducted a series of rigid model tests on the aerodynamic forces of circular steel tubular lattice structures under smoothly uniform flow conditions in the wind tunnel laboratory， and the models were fabricated based on a super high-rise power transmission tower in Huainan-Nanjing-Shanghai UHVAC transmission line. The Large-eddy simulation（LES） of lattice structures were performed to simulate the aerodynamic forces in the longitudinal， transverse and lateral directions， and the results were compared to the corresponding specifications in the available wind codes or standards. In addition， the aerodynamic forces of lattice structures were simulated by using LES under the flow conditions with turbulence intensities of 5%， 10%， 15% and 20%. The results showed that the numerical results were in relatively reasonable agreement with the experimental measurements， and the incoming turbulence intensity had indispensable influences on the fluctuating components of the aerodynamic forces. Moreover， the maximum positive pressures were observed on the windward surface at the junction of the strut and the diagonal member， while the maximum negative pressures were recorded on the lateral surface of the strut. Furthermore， the vorticity analysis showed that the vorticity in the x-direction was more scattered while those in the y- and z- directions were more uniform and continuous. The findings of this study are of great use for evaluating the aerodynamic forces of circular steel tubular lattice structures in the practice constructions.

Key words： lattice towers；LES；aerodynamic force；transmission towers

格構(gòu)式塔架由于具有受力合理、節(jié)省材料、風荷載體型系數(shù)小、回轉(zhuǎn)半徑大等結(jié)構(gòu)優(yōu)點，應用遍及輸送電、通訊信息以及廣播、電視信號傳送等各個領(lǐng)域.然而，格構(gòu)式塔架輕質(zhì)、高柔、小阻尼的結(jié)構(gòu)特性，決定了風荷載是其設(shè)計中的主要控制性荷載[1-2].因此，研究格構(gòu)式塔架的風荷載特性以及作用機理，可為設(shè)計提供荷載取值依據(jù)，對保證結(jié)構(gòu)的安全具有重要意義.

國內(nèi)外學者通過風洞試驗的方法對格構(gòu)式塔架結(jié)構(gòu)的風荷載特性進行了大量的研究.Holmes等[3-5]基于準定常假定，深入研究了順風向格構(gòu)式塔架的風振響應.Bayar等[6]通過風洞試驗的方法，采用剛性模型高頻天平測力試驗，開展了對格構(gòu)式塔架靜態(tài)風效應的研究.Celio等[7]基于試驗的塔架平均和脈動阻力系數(shù)，對比研究了風向角、擋風系數(shù)、遮蔽效應以及紊流度與阻力系數(shù)的關(guān)系.鄒良浩等[8]則開展了對輸電塔、通訊塔以及電視塔三種不同形式格構(gòu)式塔架，在動力風荷載下三分力系數(shù)的研究.樓文娟、程志軍等[9-10]通過氣動彈性模型的風洞試驗，對格構(gòu)式塔架的風振響應進行了研究.風洞試驗雖然能夠提供準確的氣動力系數(shù)，但參數(shù)化分析成本較高，并且無法給出流場結(jié)構(gòu)特征.除采用風洞試驗研究外，數(shù)值模擬計算也成為了研究的重要手段之一.然而，目前采用數(shù)值風洞來模擬格構(gòu)式塔架風荷載特性的文獻比較少見，謝華平等[11]通過CFD技術(shù)對角鋼格構(gòu)式塔架進行數(shù)值模擬，分析了計算域、網(wǎng)格劃分、湍流模型和湍流強度等因素對塔架平均風壓的影響，但模擬的結(jié)果缺乏試驗的驗證.周志勇等[12]基于CFD技術(shù)采用雷諾時均法（RANS）對河南電視塔的體型系數(shù)進行了數(shù)值模擬.RANS方法的湍流模型中包含了人為的假定，不能完全反映實際的湍流運動特征.

為更好地適應目前計算機的計算能力，較真實和精確地反映湍流的運動、發(fā)展和耗散，本文采用大渦模擬法（LES）對格構(gòu)式塔架進行數(shù)值模擬，研究圓管截面格構(gòu)式塔架節(jié)段的靜風三分力系數(shù)，通過與相關(guān)實驗數(shù)據(jù)進行對比[13]，為實際工程提供一定的參考.

1 數(shù)值計算模型及工況

1.1 數(shù)值模型

計算域及模型尺寸如圖1所示，節(jié)段高度方向與向上的z方向一致，計算域高度與節(jié)段高度一致，取474 mm.其中，主材圓管直徑為D = 57 mm，橫桿和斜材的直徑分別為17 mm和20 mm.計算域入口邊界為速度入口，設(shè)定常流速U0 = 18.4 m/s，出口為壓力出口，上下邊界以及左右邊界為對稱邊界，桿件表面為無滑移墻面邊界（no-slip wall）.根據(jù)國內(nèi)外相關(guān)圓柱繞流大渦模擬研究[14-15]，工程結(jié)構(gòu)風荷載的大渦模擬宜采用一方程亞格子模型，例如動力Smagorinsky-Lilly模型，計算時間步長為5 × 10-4 s.

為在滿足阻塞比的情況下，盡可能地減少計算網(wǎng)格數(shù)，因此以節(jié)段模型底面中心為原點，計算域上下兩側(cè)各取20D，來流方向取20D，出流方向上為保證尾流的充分發(fā)展取為40D.由于橫桿、斜桿和主材（豎桿）之間交叉接觸的地方夾角較小，且三者又是不同直徑的圓截面桿件相交，要實現(xiàn)外流場結(jié)構(gòu)化網(wǎng)格的劃分十分困難.另一方面，即使采用結(jié)構(gòu)化網(wǎng)格也無法保證足夠高的網(wǎng)格質(zhì)量，且可能存在網(wǎng)格分布方向與流體流向不一致時產(chǎn)生的偽耗散，故核心區(qū)域的網(wǎng)格采用的是非結(jié)構(gòu)化網(wǎng)格，而核心區(qū)外的網(wǎng)格則采用的是結(jié)構(gòu)化網(wǎng)格.為了保證計算結(jié)果的精度，通過加密核心區(qū)近壁面處的網(wǎng)格并增設(shè)棱柱型邊界層網(wǎng)格，來保證近壁面第一層網(wǎng)格y+≤1，從而確保取得誤差較小精度較高的計算結(jié)果.計算域網(wǎng)格劃分如圖2所示：

為了準確模擬來流風場的脈動風速譜特性，采用Huang等提出的Discretizing and Synthetic Random Flow Field Generation method （簡稱DSRFG）方法生成滿足卡門譜的入口湍流.該方法基于嚴格的理論推導，能產(chǎn)生滿足任意形式功率譜及各向異性的湍流脈動風速場，可調(diào)空間關(guān)聯(lián)性、通用性， 能夠嚴格滿足流體連續(xù)性條件，從而保證了大渦模擬計算的穩(wěn)定性，易于并行化處理[16-17].模型前方的脈動風速譜如圖3所示，結(jié)果表明DSRFG生成的脈動風速譜和目標卡門譜吻合良好，準確地再現(xiàn)了風洞試驗的來流脈動風速譜特性.

1.2 計算工況

將豎直節(jié)段模型在光滑均勻流作用下（無湍流）的數(shù)值模擬結(jié)果與風洞試驗結(jié)果進行對比分析.同時還對同風速下，湍流度分別為5%、10%、15%和20%的工況進行模擬計算，討論研究湍流度對圓管截面格構(gòu)式塔架三分力系數(shù)的影響，詳細計算工況如表1所示.

截斷模型由4根豎桿（主管），8根斜桿以及4根橫桿組成.由于各個桿件的截面尺寸不相同，故各個桿件在同一工況中的雷諾數(shù)也相互不同，各桿的雷諾數(shù)具體情況如表2所示.

2 計算結(jié)果及分析

2.1 數(shù)值模擬與試驗結(jié)果比較

CFD計算結(jié)果與文獻[11]中風洞試驗三分力的對比如表3所示.

從表3中可以看出，豎直節(jié)段模型的三分力系數(shù)與試驗值存在一定的差異，其中阻力系數(shù)CD與試驗值的誤差為9.4%，升力系數(shù)CL和扭轉(zhuǎn)系數(shù)CM均小于試驗值.分析造成二者誤差的原因，有以下幾點：

1）試驗模型在輔材相交的部分采用了方形鋼夾片進行連接固定，使得迎風面積略有增加，經(jīng)計算會造成約3%的誤差；

2）一方面，試驗模型由于加工制作的關(guān)系，所有桿件尺寸不可能做到與設(shè)計完全一致，連接固定的地方也不可能完全嚴絲密縫且完全剛接，而數(shù)值模型卻是完全一致且在幾何上完全對稱的；另一方面，風洞試驗時模型的迎風角度不可能完全準確，可能出現(xiàn)較小的偏角.這兩方面的原因都會造成試驗結(jié)果的升力系數(shù)CL和扭轉(zhuǎn)系數(shù)CM比數(shù)值模擬偏大一些.

CFD計算結(jié)果與各國規(guī)范[18-24]的對比如圖4所示.

其中密實比的定義如下：

按照上述公式得出試驗模型的Φ為0.314.圖4中與CFD計算結(jié)果對比的規(guī)范值，取的是0°風向角下方形塔架體型系數(shù)與密實比Φ的關(guān)系曲線，從圖中能夠看出，本文的計算結(jié)果與規(guī)范IEC60826給定值最為接近.

2.2 湍流度對三分力系數(shù)的影響

圖5（a）（b）給出了模型的三分力系數(shù)隨來流湍流度增加的變化情況，從圖中可以看出，湍流度的變化對模型三分力系數(shù)均值的影響較小.湍流度的增大對模型三分力系數(shù)的脈動分量影響較大，湍流度的增大會導致模型三分力系數(shù)的脈動分量隨之變大，脈動幅值增大.

圖6給出了不同來流湍流度情況下模型三分力系數(shù)的功率譜，通過頻譜上的分析能夠更清楚地看到，三個方向的氣動力功率譜隨著頻率的增大先緩慢減小最后急劇減小.并且隨著湍流度的增加，在低頻部分結(jié)構(gòu)三分力系數(shù)的脈動幅值逐漸增大.

將三分力功率譜S*i（f）按如下表達式進行歸一化處理，結(jié)果如圖6所示.

圖7表明，隨著湍流度的增加，功率譜峰值及其對應卓越頻率逐漸前移，且湍流度對阻力譜高頻部分影響顯著.

2.3 壓力分布云圖

圖8分別給出了模型迎風面和背風面的風壓系數(shù)分布云圖.從圖中可以看出，節(jié)段模型的迎風面受正壓，側(cè)風面和背面受負壓.其中，結(jié)構(gòu)所受最大正風壓系數(shù)約為1.0，出現(xiàn)在迎風面的兩根主管正面和輔材相交的連接處，因為這兩個局部的相對過流截面較小，來流會在這兩個地方產(chǎn)生較大的風壓.而由于上游主管和輔材的遮擋效應，可以明顯地看出下游桿件迎風面的正風壓分布是小于上游桿件的.

結(jié)構(gòu)所受的最大負壓出現(xiàn)在主管的側(cè)面，從圖8（b）中可以較清楚地看出，斜桿和水平橫桿對負壓的分布是有影響的，在主管無輔材的一側(cè)，負壓分布呈長條形，而在有輔材的一側(cè)，負壓分布呈橢圓形，在輔材的位置處負壓變小回收，表明流體的再附點提前了.就背風面而言，主管和輔材上的壓力分布并不十分規(guī)整，這是由于主管和側(cè)面斜桿的間距是不斷變化的，互相會對彼此的尾流渦脫產(chǎn)生影響，導致其壓力分布的不均勻性.

2.4 流場結(jié)構(gòu)分析

在流場結(jié)構(gòu)的分析中，渦量反映了旋渦的能量強度和方向，描述了速度矢量場的特性，因此渦量圖較速度分布云圖能夠更直觀地表達結(jié)構(gòu)的渦脫特征，給出更多的流場結(jié)構(gòu)信息，因此，可以通過渦量圖的分析來描述節(jié)段模型的流場結(jié)構(gòu)特征.

圖9分別給出了模型在三個方向的渦量圖，可以清楚地看到流體在流經(jīng)結(jié)構(gòu)桿件時的附著和分離，并形成隨機脫落的旋渦.結(jié)果表明x向的渦量離散程度更高，y和z向上渦量相對連續(xù)性更高一些.

3 結(jié) 論

采用大渦模擬方法對圓管截面格構(gòu)式塔架節(jié)段模型進行了數(shù)值模擬計算，同時為了研究來流湍流度對結(jié)構(gòu)三分力系數(shù)的影響，對湍流度分別為5%、10%、15%和20%時的工況進行分析.通過數(shù)值模擬的手段，對結(jié)構(gòu)的風壓分布、流場結(jié)構(gòu)進行了更細致的研究，經(jīng)過對比分析得到了以下的結(jié)論：

1）采用大渦模擬方法對格構(gòu)式塔架節(jié)段模型進行數(shù)值模擬，能夠得到較為滿意的計算結(jié)果，與風洞的試驗值對比，阻力系數(shù)CD的誤差能夠控制在10%以內(nèi)，并與規(guī)范IEC60826給定值最為接近，數(shù)值模擬計算的升力系數(shù)CL和扭轉(zhuǎn)系數(shù)CM均略小于試驗值.

2）來流湍流度的增加，對模型三分力系數(shù)的均值基本沒有影響，但脈動分量會隨之增加，即脈動幅值變大.同時，隨著湍流度的增加，三分力系數(shù)功率譜峰值對應的卓越頻率逐漸前移，且湍流度對功率譜高頻部分影響顯著.

3）通過對節(jié)段模型的壓力分布的分析發(fā)現(xiàn)，最大正壓出現(xiàn)在迎風面的主管正面和輔材相交節(jié)點的正面，最大負壓出現(xiàn)在主管的側(cè)面.

4）通過渦量圖對節(jié)段模型流場結(jié)構(gòu)的分析發(fā)現(xiàn)，x向的渦量離散程度更高，y和z向上渦量相對連續(xù)性更高一些.

參考文獻

[1] PARK J H，MOON B W，MIN K W.et al.Cyclic loading test of friction-type reinforcing members upgrading wind-resistant performance of transmission towers[J].Engineering Structures，2007， 29（11）：3185—3196.

[2] 梁樞果，鄒良浩， 趙林，等. 格構(gòu)式塔架動力風荷載解析模型[J].同濟大學學報（自然科學版），2008，36（2）：166—171.

LIANG S G， ZOU L H， ZHAO L，et al.Analytical model of dynamic wind loads on lattice towers [J].Journal of Tongji University（Natural Science），2008，36（2）：166—171.（In Chinese）

[3] HOLMES J D.Along-wind responses of lattice tower： part I—derivation of expressions for gust response factors [J]. Engineering Structures，1994，6（4）：287—292.

[4] HOLMES J D. Along-wind responses of lattice tower—Ⅱ.aerodynamics damping and deflections [J]. Engineering Structures， 1996，18（7）：483—488.

[5] HOLMES J D. Along-wind responses of lattice tower—Ⅲ.effective load distributions[J]. Engineering Structures，1996，18（7）：489—494.

[6] BAYAR D C. Drag coefficients of latticed towers[J]. Journal of Structural Engineering， 1986， 112 （2）： 417—430.

[7] CELIO F C J，ISYUMOV N，BRASIL R M L R F.Experimental study of the wind forces on rectangular latticed communication towers with antennas[J].Journal of Wind Engineering and Industrial Aerodynamics，2003，91（8）：1007—1022.