嚴(yán) 磊, 朱樂東,3

(1. 同濟(jì)大學(xué) 橋梁工程系, 上海 200092; 2. 同濟(jì)大學(xué) 土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 上海 200092; 3. 同濟(jì)大學(xué) 橋梁結(jié)構(gòu)抗風(fēng)技術(shù)交通行業(yè)重點(diǎn)實(shí)驗(yàn)室, 上海 200092)

格柵湍流場(chǎng)風(fēng)參數(shù)沿風(fēng)洞軸向變化規(guī)律

嚴(yán) 磊1,2, 朱樂東1,2,3

(1. 同濟(jì)大學(xué) 橋梁工程系, 上海 200092; 2. 同濟(jì)大學(xué) 土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 上海 200092; 3. 同濟(jì)大學(xué) 橋梁結(jié)構(gòu)抗風(fēng)技術(shù)交通行業(yè)重點(diǎn)實(shí)驗(yàn)室, 上海 200092)

格柵湍流場(chǎng)常用于各類風(fēng)洞試驗(yàn)中，橋梁斷面氣動(dòng)參數(shù)均與來流風(fēng)參數(shù)有關(guān)。對(duì)不同格柵湍流場(chǎng)下的風(fēng)參數(shù)沿風(fēng)洞軸向變化的規(guī)律進(jìn)行了研究，得出風(fēng)參數(shù)主要與測(cè)點(diǎn)與格柵斷面間距、格柵板條厚度和單元格柵邊長(zhǎng)有關(guān)。隨著間距增加，平均風(fēng)速短期內(nèi)急劇衰減，然后趨于穩(wěn)定。隨著間距增加，來流湍流的總能量減小而主要能量遷移至高頻區(qū)域：湍流度呈指數(shù)衰減，湍流積分尺度呈增長(zhǎng)趨勢(shì)；約化風(fēng)譜形狀保持相似，峰值對(duì)應(yīng)約化頻率值基本保持不變，但峰值變小。另外根據(jù)風(fēng)參數(shù)的變化規(guī)律，調(diào)試出2類特定的湍流場(chǎng)參數(shù)：湍流度相似而湍流積分尺度不同的湍流場(chǎng)和湍流積分尺度相似而湍流度不同的湍流場(chǎng)，以便于詳盡研究主要風(fēng)參數(shù)對(duì)橋梁斷面氣動(dòng)參數(shù)的影響。

格柵湍流場(chǎng)；多點(diǎn)同步測(cè)量；風(fēng)參數(shù)；湍流度；湍流積分尺度；脈動(dòng)風(fēng)速自功率譜；風(fēng)洞試驗(yàn)

0 引 言

格柵湍流場(chǎng)常用于湍流場(chǎng)節(jié)段模型測(cè)力、測(cè)壓和測(cè)振試驗(yàn)中，用來測(cè)量橋梁斷面抖振力系數(shù)譜、氣動(dòng)導(dǎo)納函數(shù)、跨向相關(guān)性參數(shù)、湍流場(chǎng)中氣動(dòng)導(dǎo)數(shù)和顫振臨界風(fēng)速，并研究湍流度和湍流積分尺度等湍流場(chǎng)風(fēng)參數(shù)對(duì)上述橋梁斷面氣動(dòng)參數(shù)的影響[1-8]。在這類試驗(yàn)中，一方面要求格柵湍流場(chǎng)是空間均勻的，即在模型所在位置，平均風(fēng)速、湍流度、湍流積分尺度和脈動(dòng)風(fēng)速自功率譜等湍流場(chǎng)風(fēng)參數(shù)在風(fēng)洞截面內(nèi)的離散度必須充分小。然而，由于格柵湍流場(chǎng)的湍流是格柵對(duì)來流的干擾和旋渦脫落所產(chǎn)生，而格柵的桿件是離散布置的，因此這種擾動(dòng)在風(fēng)洞截面內(nèi)從一開始就是不均勻的，但由于擾動(dòng)和旋渦的橫向擴(kuò)散與混合，隨著下游截面離格柵的距離的增加，湍流場(chǎng)風(fēng)參數(shù)在截面內(nèi)的均勻性逐步改善；另一方面要求湍流場(chǎng)風(fēng)參數(shù)，湍流強(qiáng)度和湍流積分尺度是分離的，即要保證流場(chǎng)的湍流強(qiáng)度相互近似但湍流積分尺度不同或者流場(chǎng)的湍流積分尺度相互近似但湍流強(qiáng)度不同，以便研究湍流場(chǎng)風(fēng)參數(shù)對(duì)橋梁斷面氣動(dòng)參數(shù)的影響。然而，目前湍流理論的研究尚處于數(shù)學(xué)模型階段，其實(shí)際應(yīng)用缺乏可靠的理論指導(dǎo)和可移植的試驗(yàn)結(jié)果?，F(xiàn)階段研究人員可以對(duì)不同的湍流發(fā)生裝置(格柵，尖劈和擋板等)進(jìn)行組合[9]以確定合適布置，但是組合后湍流場(chǎng)脈動(dòng)風(fēng)速約化風(fēng)譜差異較大，影響后面的參數(shù)分析；另外也可以僅對(duì)某一種湍流發(fā)生裝置進(jìn)行調(diào)試，例如在格柵湍流場(chǎng)下通過改變格柵斷面間距、格柵板條厚度和單元格柵邊長(zhǎng)來確定合適布置。由于試驗(yàn)條件的限制和調(diào)試工作的繁瑣，就目前作者的認(rèn)識(shí)水平，僅僅文獻(xiàn)[10-11]在格柵湍流場(chǎng)下調(diào)試出了湍流度相似而湍流積分尺度不同的湍流場(chǎng)，遺憾的是研究者并未對(duì)風(fēng)參數(shù)沿風(fēng)洞軸向變化的規(guī)律進(jìn)行總結(jié)，也未調(diào)試出湍流積分尺度相似而湍流度不同的湍流場(chǎng)。所以本文通過在同濟(jì)大學(xué)土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室TJ-2風(fēng)洞中進(jìn)行的不同被動(dòng)格柵湍流場(chǎng)測(cè)試工作，對(duì)湍流場(chǎng)風(fēng)參數(shù)變化規(guī)律進(jìn)行了研究。然后依據(jù)該規(guī)律，調(diào)試出2類特定的湍流場(chǎng)參數(shù)：湍流度相似湍流積分尺度不同的湍流場(chǎng)和湍流積分尺度相似湍流度不同的湍流場(chǎng)，為下一步研究湍流場(chǎng)風(fēng)參數(shù)對(duì)橋梁斷面氣動(dòng)參數(shù)的影響規(guī)律打下基礎(chǔ)。

1 湍流場(chǎng)風(fēng)參數(shù)

湍流場(chǎng)風(fēng)參數(shù)主要包括平均風(fēng)速、湍流度、湍流積分尺度以及脈動(dòng)風(fēng)速自功率譜。其中湍流度是風(fēng)速脈動(dòng)強(qiáng)度的一個(gè)指標(biāo)，其被定義為脈動(dòng)風(fēng)速的標(biāo)準(zhǔn)差與平均風(fēng)速的比值，即：

(1)

(2)

(3)

式中：λ為擬合參數(shù)；r為y或z方向兩個(gè)測(cè)點(diǎn)間的間距。

脈動(dòng)風(fēng)速a沿r方向的湍流積分尺度則為：

(4)

由于沿x方向同時(shí)布置多個(gè)測(cè)風(fēng)速探針時(shí)，上游探針會(huì)對(duì)下游探針產(chǎn)生較明顯的干擾，因此在實(shí)際應(yīng)用中湍流風(fēng)速a沿x方向的湍流積分尺度采用泰勒的“渦流凍結(jié)傳輸”假說[12]進(jìn)行計(jì)算，即：假設(shè)湍流中的漩渦是不衰減地以平均風(fēng)速向下游傳輸，由此湍流積分尺度的公式可以轉(zhuǎn)化為時(shí)間尺度的積分：

(5)

(6)

式中λ為擬合參數(shù)。

這樣，便可以得到風(fēng)場(chǎng)中脈動(dòng)風(fēng)速a沿順風(fēng)向x軸的湍流積分尺度：

(7)

圖1和2分別為典型歸一化協(xié)方差曲線和典型歸一化自協(xié)方差曲線。

圖1 典型歸一化協(xié)方差曲線

圖2 典型歸一化自協(xié)方差曲線

脈動(dòng)風(fēng)速自功率譜(簡(jiǎn)稱風(fēng)譜)描述的是湍流運(yùn)動(dòng)能量隨頻率或波長(zhǎng)的分布情況，并且其一般都被表示為約化形式，因?yàn)榧s化風(fēng)譜的曲線至少有一個(gè)與湍流典型尺度相對(duì)應(yīng)的峰值。直接使用風(fēng)洞來流實(shí)測(cè)風(fēng)譜這樣離散數(shù)據(jù)不僅麻煩，而且可能會(huì)帶來顯著的誤差。因此，一般需要對(duì)其進(jìn)行帶約束條件的最小二乘曲線擬合[13]，而3個(gè)方向的脈動(dòng)風(fēng)速約化風(fēng)譜均可采用以下的目標(biāo)函數(shù)來擬合：

(8)

(9)

其中，fz為一種無量綱約化頻率，稱為莫寧坐標(biāo)。a,b,m,γ為4個(gè)擬合參數(shù)，一般令γ為1。c一般為常系數(shù)5/3，而待擬合系數(shù)a,b,m必須滿足下列約束條件：

(10)

(11)

這里，Γ為伽馬函數(shù)。

2 試驗(yàn)概況

TJ-2風(fēng)洞的試驗(yàn)段為3.0m×2.5m×15m(寬×高×長(zhǎng))，可調(diào)風(fēng)速最大可達(dá)65m/s，流場(chǎng)不均勻指標(biāo)δU/U≤1.0%，湍流度Iu≤0.5%，氣流水平向和豎向偏角均小于0.5°。

試驗(yàn)中采用澳大利亞Turbulence Flow公司生產(chǎn)的Series100眼鏡蛇探頭測(cè)量各點(diǎn)處的風(fēng)速時(shí)程數(shù)據(jù)。眼鏡蛇探頭設(shè)置參數(shù)分別為采集頻率5kHz，輸出頻率1kHz，采樣點(diǎn)數(shù)65536。坐標(biāo)系定義如下：沿著風(fēng)洞來流方向?yàn)閤方向，垂直向上為z方向，y為水平方向，按右手法則確定(指向TJ-2風(fēng)洞操作臺(tái))，坐標(biāo)系原點(diǎn)設(shè)在格柵平面與風(fēng)洞底板交線的中點(diǎn)；此坐標(biāo)系與眼鏡蛇探頭約定脈動(dòng)風(fēng)速方向一致。

試驗(yàn)中采用沿著y或z方向一次布置7個(gè)眼鏡蛇探頭進(jìn)行多點(diǎn)同步采樣，而探頭之間的最小間距一般需要滿足2個(gè)條件，一方面小于y或z方向的湍流積分尺度，另一方面大于探頭尺寸的10倍距離。因此非參考探頭相對(duì)參考探頭交錯(cuò)布置，且距離分別設(shè)置為3.35、9.0、11.6、17.6、25.5和35cm，其布置情況如圖3所示。對(duì)于脈動(dòng)風(fēng)速a沿y、z方向的6個(gè)湍流積分尺度，采用多點(diǎn)同步測(cè)量法計(jì)算；對(duì)于脈動(dòng)風(fēng)速a沿x方向的3個(gè)湍流積分尺度采用泰勒的“渦流凍結(jié)傳輸”假說計(jì)算。

圖3 眼鏡蛇探頭布置情況

3 試驗(yàn)結(jié)果

3.1 典型格柵湍流場(chǎng)風(fēng)特性參數(shù)沿風(fēng)洞軸向變化

(12)

式中：n為頻率，Lax為x方向湍流積分尺度，U為平均風(fēng)速。

平均風(fēng)速在短期內(nèi)衰減非?？?，這是因?yàn)闅饬鹘?jīng)格柵斷面時(shí)能通過的斷面面積急劇減少，而風(fēng)機(jī)短時(shí)間內(nèi)提供的能量不變，單位時(shí)間內(nèi)某斷面通過的氣體體積一定，故此時(shí)會(huì)有明顯的加速過程，而氣團(tuán)離開格柵斷面時(shí)，同理會(huì)明顯地減小；之后則衰減很慢，甚至不變。

圖4 編號(hào)6×5×05格柵

圖5 平均風(fēng)速、湍流度和湍流積分尺度隨歸一化距離變化情況

Fig.5 Mean wind velocity, turbulence intensity and integral length scale variation with the non-dimensional distance

圖6 順風(fēng)向脈動(dòng)擬合風(fēng)速譜隨歸一化距離變化情況

Fig.6 Fitting non-dimensional power-spectral density for the longitudinal turbulence component variation with the non-dimensional distance

圖7 水平風(fēng)向脈動(dòng)擬合風(fēng)速譜隨歸一化距離變化情況

Fig.7 Fitting non-dimensional power-spectral density of the lateral turbulence component variation with non-dimensional distance

圖8 豎直風(fēng)向脈動(dòng)擬合風(fēng)速譜隨歸一化距離變化情況

Fig.8 Fitting non-dimensional power-spectral density of the vertical turbulence component variation with non-dimensional distance

圖9 編號(hào)6×5×05+6×5×05+6×5×05格柵

約化風(fēng)譜形狀保持相似，但峰值變小，而峰值對(duì)應(yīng)約化頻率值基本保持不變。此表明沿離格柵越遠(yuǎn)，風(fēng)速脈動(dòng)能量越小，這與湍流度衰減機(jī)理一致。

3.2 兩類特定湍流場(chǎng)組

本試驗(yàn)還進(jìn)行了多組格柵湍流的調(diào)試，甚至包括3層格柵的調(diào)試。3層格柵裝置如圖9所示，中間格柵為參考格柵斷面，沿著來流方向的前后格柵斷面至中間參考斷面的間距分別為410cm和200cm。根據(jù)湍流風(fēng)參數(shù)對(duì)影響變量的不同依賴程度，給出2類特定的湍流場(chǎng)參數(shù)，即湍流度大體相同，湍流積分尺度相差較大；湍流積分尺度大體相同，湍流度相差較大。表1～6給出了相應(yīng)的調(diào)試結(jié)果，而表中湍流場(chǎng)類型項(xiàng)中6×5×05D200的D200表示測(cè)點(diǎn)距參考格柵斷面200cm，而6×5×05與上文含義一致。

表1 湍流度Iu約為11%而湍流積分尺度不同的湍流場(chǎng)組

表2 湍流度Iu約為8%而湍流積分尺度不同的湍流場(chǎng)組

表3 湍流度Iu約為6%而湍流積分尺度不同的湍流場(chǎng)組

表4 湍流度Iu約為4%而湍流積分尺度不同的湍流場(chǎng)組

表5 湍流積分尺度約為0.16 m而湍流度不同的湍流場(chǎng)組

表6 湍流積分尺度約為0.13 m而湍流度不同的湍流場(chǎng)組

將測(cè)試結(jié)果中Iu接近于11%、8%、6%、4%的格柵類型分別歸為一組，每組篩選出湍流積分尺度相差較大的3種格柵類型，從而獲得4組湍流度大致相同而紊流積分尺度不同的湍流場(chǎng)類型，如表1～4所示。

從中可以看出，在湍流度大致相同的基礎(chǔ)上獲得較大湍流積分尺度可以通過以下3種途徑：(1) 通過增大板條寬度B并且增大測(cè)點(diǎn)至參考格柵斷面的距離D；(2) 通過增加格柵的層數(shù)并且增大測(cè)點(diǎn)至參考格柵斷面的距離D；(3) 通過增大單元格子的邊長(zhǎng)A并且增大測(cè)點(diǎn)至參考格柵斷面的距離D。

從中可以看出，在湍流積分尺度大致相同的基礎(chǔ)上獲得較大湍流度可以通過以下3種途徑：(1) 通過增大板條寬度B但減小測(cè)點(diǎn)至參考格柵斷面的距離D；(2) 通過增加格柵的層數(shù)但減小測(cè)點(diǎn)至參考格柵斷面的距離D；(3) 通過增大單元格子的邊長(zhǎng)A但減小測(cè)點(diǎn)至參考格柵斷面的距離D。

4 結(jié) 論

從試驗(yàn)結(jié)果可以得到以下格柵湍流場(chǎng)的風(fēng)參數(shù)隨截面離格柵距離增加的變化規(guī)律：

(5) 根據(jù)湍流風(fēng)參數(shù)對(duì)影響變量的不同依賴程度，可以調(diào)試出給出2類特定的湍流場(chǎng)參數(shù)，即湍流度大體相同，湍流積分尺度相差較大；湍流積分尺度大體相同，湍流度相差較大。為下一步研究湍流場(chǎng)風(fēng)參數(shù)對(duì)橋梁抖振力參數(shù)的影響規(guī)律打下基礎(chǔ)。

[1] 文水兵, 朱樂東, 丁泉順. 扁平雙邊肋主梁橋梁斷面抖振力及氣動(dòng)導(dǎo)納頻譜特性研究[J]. 結(jié)構(gòu)工程師, 2008, 24(3): 96-101.

Wen S B, Zhu L D, Ding Q S. Research on buffeting force and aerodynamic admittance spectra characteristics for flat double-grid bridge decks[J]. Structural Engineers, 2008, 24(3): 96-101.

[2] 朱樂東, 趙傳亮, 丁泉順. 扁平箱主梁斷面抖振力頻譜特性以及橋跨方向相關(guān)性研究[J]. 結(jié)構(gòu)工程師, 2009, 25(4): 96-101.

Zhu L D, Zhao C L, Ding Q S. Spectral characteristics and spanwise coherence of buffeting forces on a flat box-girder bridge deck[J]. Structural Engineers, 2009, 25(3): 96-101.

[3] 周奇. 考慮特征湍流效應(yīng)和力空間相關(guān)性的橋梁抖振分析[D]. 上海: 同濟(jì)大學(xué), 2011.

Zhou Q. Buffeting analysis of bridges with signature turbulence effect and spatial coherence of buffeting forces [D]. Shanghai: Tongji University, 2011.

[4] 徐自然, 周奇, 朱樂東. 考慮模型抖振力跨向不完全相關(guān)性效應(yīng)的氣動(dòng)導(dǎo)納識(shí)別[J] . 實(shí)驗(yàn)流體力學(xué), 2014, 28(4): 1-8.

Xu Z R, Zhou Q, Zhu L D. Identification of aerodynamic admittances by considering the effect of incomplete span-wise correclation of buffeting forces on sectional model[J]. Journal of Experiments in Fluid Mechanics, 2014, 28(4): 1-8.

[5] Kimura K, Fujino Y, Nakato S. Characteristics of buffeting forces on flat cylinders[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1997, 69: 365-374.

[6] Sarkar P P, Jones N P, Scanlan R H. Identification of aeroelastic parameters of flexible bridges[J]. Engineering Mechanics, ASCE, 1994, 120(8): 1718-1741.

[7] Scanlan R H, Lin W H. Effects of turbulence on bridge flutter derivatives[J]. Journal of the Engineering Mechanics Division, 1978, 104(4): 719-733.

[8] Bucher C G, Lin Y K. Effects of wind turbulence on motion stability of long-span bridges[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1988, 36: 1355-1364.

[9] Larose G L, Tanaka H, Gimsing N J, et al. Direct measurements of buffeting wind forces on bridge decks[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1998, 74-76: 809-818.

[10] Haan Jr F L, Kareem A. Anatomy of turbulence effects on the aerodynamics of an oscillating prism[J]. Journal of Engineering Mechanics, 2009, 135(9): 987-999.

[11] Hui M C. Turbulent wind action on long span bridges with separated twin-girder decks[D]. Shanghai: Tongji University, 2002.

[12] Batchelor G K. The theory of homogeneous turbulence[M]. London: Cambridge University Press, 1953.

[13] Zhu L D. Buffeting response of long span cable-supported bridges under skew winds: field measurement and analysis[D]. HongKong: The Hong Kong Polytechnic University, 2002.

(編輯：李金勇)

Wind characteristics of grid-generated wind field along the wind tunnel

Yan Lei1,2, Zhu Ledong1,2,3

(1. Department of Bridge Engineering, Tongji University, Shanghai 200092, China;2. State Key Laboratory of Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China; 3. Key Laboratory for Wind Resistance Technology of Bridges of Ministry of Transport, Tongji University, Shanghai 200092, China)

Grid-generated wind fields are commonly applied in various wind tunnel tests to obtain aerodynamic parameters of flexible bridges. In this work, wind characteristics of grid-generated wind field were investigated at different distances from the grid section along the wind tunnel. The wind characteristics principally depend on the distance, the thickness of slat and the length of the unit grid side. The results indicate that the mean wind velocity descends steeply and soon afterwards steadies with the distance. The total turbulent energy declines and the dominant energy migrates to the higher frequency region with the distance: the turbulence intensity declines and the integral length scale increases; non-dimensional power-spectral density forms and non-dimensional frequencies corresponding to the peak value remain unchanged, but the peak values decline. Two specific turbulent conditions were adjusted according to these contributing factors: similar turbulence intensities with different integral length scales; similar integral length scales with different turbulence intensities, founding the basis to investigate the effects of main wind characteristics on the aeroelastic parameters.

grid-generated wind filed;multipoint simultaneously measurements;wind characteristics;turbulence intensities;integral length scales;power spectral density of turbulence;wind tunnel test

1672-9897(2015)01-0049-06

10.11729/syltlx20140075

2014-07-05;

2014-08-12

科技部國(guó)家重點(diǎn)實(shí)驗(yàn)室基礎(chǔ)研究資助項(xiàng)目(SLDRCE08-A-02);國(guó)家自然科學(xué)基金項(xiàng)目(91215302)

YanL,ZhuLD.Windcharacteristicsofgrid-generatedwindfieldalongthewindtunnel.JournalofExperimentsinFluidMechanics, 2015, 29(1): 49-54. 嚴(yán) 磊, 朱樂東. 格柵湍流場(chǎng)風(fēng)參數(shù)沿風(fēng)洞軸向變化規(guī)律. 實(shí)驗(yàn)流體力學(xué), 2015, 29(1): 49-54.

P404; TU317.1

A

嚴(yán) 磊(1988-)，男，湖南岳陽人，博士研究生。研究方向：橋梁和建筑結(jié)構(gòu)抗風(fēng)研究。通信地址：上海市四平路1239號(hào)同濟(jì)大學(xué)橋梁系橋梁館309室(200092)。E-mail: ledong@#edu.cn