串列雙方柱的風(fēng)壓特性及其流場機(jī)理

杜曉慶　陳麗萍　董浩天　田新新　許漢林

摘 ? 要：為了研究干擾條件下方柱的風(fēng)壓特性及其流場機(jī)理，以串列雙方柱為研究對象，采用大渦模擬方法，在雷諾數(shù)Re = 8 × 104、間距比P/B = 1.1 ～ 5的條件下，研究了兩個方柱的風(fēng)壓系數(shù)、氣動力系數(shù)、風(fēng)壓非高斯特性、風(fēng)壓相關(guān)性隨間距比的變化規(guī)律，重點(diǎn)探討了雙方柱流場特性及其與風(fēng)壓非高斯特性的內(nèi)在聯(lián)系. 研究結(jié)果表明：隨著間距比的增大，串列雙方柱依次表現(xiàn)為3種流態(tài)，即單一鈍體、剪切層再附和雙渦脫流態(tài)，風(fēng)壓特性與流態(tài)密切相關(guān). 風(fēng)壓的非高斯特性和風(fēng)壓相關(guān)性隨流態(tài)變化呈現(xiàn)為3種類型：在單一鈍體流態(tài)下，柱間回流區(qū)附近的表面風(fēng)壓呈現(xiàn)明顯的非高斯特性且風(fēng)壓相關(guān)性較強(qiáng);在剪切層再附流態(tài)下，方柱尾流的渦脫強(qiáng)度低，風(fēng)壓相關(guān)性弱，但風(fēng)壓非高斯區(qū)域大;在雙渦脫流態(tài)下，受上游方柱尾流旋渦的作用，方柱側(cè)風(fēng)面的風(fēng)壓相關(guān)性較強(qiáng)，下游方柱的側(cè)風(fēng)面和背風(fēng)面出現(xiàn)大范圍的風(fēng)壓非高斯區(qū)域.

關(guān)鍵詞：串列雙方柱;氣動干擾;大渦模擬;風(fēng)壓非高斯特性;流場機(jī)理

中圖分類號：TU 311 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼：A

Wind Pressure Characteristics and Flow Mechanism

of Two Tandem Square Columns

DU Xiaoqing1，2，CHEN Liping1，DONG Haotian3?，TIAN Xinxin1，XU Hanlin1

（1. Department of Civil Engineering，Shanghai University，Shanghai 200444，China;

2. Wind Engineering and Aerodynamic Flow Control Research Center，Shanghai University，Shanghai 200444，China;

3. School of Mechanics and Engineering Science，Shanghai University，Shanghai 200444，China）

Abstract：To study the wind pressure characteristics and flow mechanism of square columns due to flow interference，taking the tandem square columns as the research object，large eddy simulation method is applied under the condition that Reynolds number Re = 8 × 104 and spacing ratio P/B=1.1～5. The influences of spacing ration on wind pressure coefficient，aerodynamic coefficient，non-Gaussian feature of wind pressure，and wind pressure correlation of the columns are discussed. The relationship between flow pattern and non-Gaussian characteristics of wind pressure is revealed. The results show that three major regimes are distinguished along with various spacing ratios，namely single bluff-body regime，shear layer reattachment regime and co-shedding regime，which are closely related with the specific wind pressure patterns. In the single bluff-body regime，the wind pressure near the recirculation zone between two columns shows significant non-Gaussian feature，whose correlation of wind pressure is strong. The shear layer reattachment regime，however，shows a weak vortex shedding in near wake，accompanied by a small wind pressure correlation but a large region of non-Gaussian wind pressure. Finally，in the co-shedding regime，due to vortex shedding of the upstream cylinder，the wind pressure on the side surface of two columns has a strong correlation and there is a large region of non-Gaussian wind pressure on the side and leeward surfaces of the downstream cylinder.

Key words：tandem square columns;aerodynamic interference;large eddy simulation;non-Gaussian features of wind pressure;flow mechanism

群體超高層建筑之間存在復(fù)雜的氣動干擾現(xiàn)象，這使得群體建筑的風(fēng)荷載與單體建筑有很大差異. 影響群體超高層建筑風(fēng)荷載的因素眾多，如來流風(fēng)特性、建筑體型、建筑數(shù)量、建筑間距、高寬比、風(fēng)向角等，情況非常復(fù)雜[1-5]，其流場干擾機(jī)理有待研究.

風(fēng)壓的非高斯特性是超高層建筑圍護(hù)結(jié)構(gòu)抗風(fēng)設(shè)計需重點(diǎn)關(guān)注的問題. 已有研究[6-7]表明，單體超高層建筑的表面風(fēng)壓存在明顯的非高斯特性，若按照傳統(tǒng)的具有高斯特性的峰值因子法估計極值風(fēng)壓，會使結(jié)構(gòu)設(shè)計偏于不安全. Ko等[8]通過模擬風(fēng)壓，研究了高斯風(fēng)荷載與非高斯風(fēng)荷載對方形建筑側(cè)面的影響，結(jié)果表明：在非高斯風(fēng)荷載作用下，建筑表面易出現(xiàn)較大的負(fù)壓極值. 不少學(xué)者通過大氣邊界層風(fēng)洞試驗(yàn)研究了單個柱體的風(fēng)壓非高斯特性. 韓寧等[9]基于風(fēng)洞試驗(yàn)分析了不同風(fēng)向角下方形截面建筑的風(fēng)壓非高斯特性，發(fā)現(xiàn)來流直接作用面主要為高斯區(qū)域，而非高斯區(qū)域主要在分離流和尾流作用處. 莊翔等[10]對矩形截面高層建筑表面風(fēng)壓脈動的非高斯特性進(jìn)行研究，劃分了兩種風(fēng)向下建筑表面風(fēng)壓脈動的高斯區(qū)和非高斯區(qū). 樓文娟等[11]研究了帶切角的菱形超高層建筑的風(fēng)壓非高斯特性，給出了劃分風(fēng)壓非高斯區(qū)域的依據(jù)，并發(fā)現(xiàn)側(cè)面前緣分離區(qū)、背風(fēng)面及迎風(fēng)面切角區(qū)域存在較強(qiáng)的風(fēng)壓非高斯特性. 上述研究均未考慮周圍建筑的干擾效應(yīng). 王浩等[12]對5種四塔方案進(jìn)行剛體測壓風(fēng)洞試驗(yàn)，結(jié)果表明：與單塔工況相比，串列、矩形、菱形、L形和斜L形方案的風(fēng)壓非高斯區(qū)域均有所增加，但并未給出其流場干擾機(jī)理.

均勻來流作用下的串列雙方柱是群體超高層建筑的簡化模型. Sohankar[13]基于數(shù)值模擬結(jié)果，將串列雙方柱繞流場分為三種流態(tài)：單一鈍體、剪切層再附和雙渦脫流態(tài). Yen等[14]研究了雷諾數(shù)、間距比及風(fēng)向角對串列方柱的影響，結(jié)果表明：其流態(tài)、Strouhal數(shù)和平均阻力系數(shù)受到較大影響. 杜曉慶等[10]通過風(fēng)洞試驗(yàn)研究了串列雙方柱間的氣動干擾，發(fā)現(xiàn)兩方柱的風(fēng)壓非高斯區(qū)域與間距有密切聯(lián)系，但尚未從流場角度解釋風(fēng)壓非高斯特性的機(jī)理. 理解干擾條件下風(fēng)壓特性的流場機(jī)理，有助于提出更好的氣動優(yōu)化措施.

本文以串列雙方柱繞流為研究對象，采用大渦模擬方法，在雷諾數(shù)8 × 104條件下（易與相同雷諾數(shù)下的風(fēng)洞試驗(yàn)結(jié)果進(jìn)行比較和驗(yàn)證），考慮9種不同間距比（P/B = 1.1 ～ 5），研究了雙方柱的風(fēng)壓系數(shù)、氣動力系數(shù)、風(fēng)壓非高斯特性、風(fēng)壓相關(guān)性隨間距比的變化規(guī)律，重點(diǎn)探討了串列雙方柱的流場干擾機(jī)理及其與風(fēng)壓非高斯特性的內(nèi)在聯(lián)系.

1 ? 計算模型和結(jié)果驗(yàn)證

1.1 ? 計算模型和計算網(wǎng)格

圖1為計算模型示意圖，來流為均勻來流，基于來流風(fēng)速U0和方柱邊長B計算得到的雷諾數(shù)為Re = 8 × 104，兩方柱的中心間距為P，共計算了9種不同間距比P/B分別為1.1、1.25、1.5、2、2.5、3、3.5、4和5.

圖2所示為計算域模型及邊界條件，本文采用半徑為30B的O型計算域，阻塞率為1.67%. 計算域入口和出口分別采用速度入口邊界條件（Velocity-inlet）和自由出口邊界條件（Outflow），方柱展向兩端采用周期性邊界條件（Periodic），方柱表面為無滑移壁面（Wall）. 采用大渦模擬（LES）湍流模型，大尺度渦通過濾波后的不可壓縮Navier-Stokes方程直接求解：

= 0 ? ? ? （1）

+=-+ν- ? ? （2）

式中：[u]i為濾波后速度分量;xi為空間位置分量;t為時間;ρ為流體密度;[p]為濾波后壓力;ν為流體動力黏度;τij為亞格子應(yīng)力張量. 采用WALE亞格子尺度模型， WALE常數(shù)為0.325. 采用SIMPLEC算法求解壓力-速度場耦合方程，空間離散方法和時間離散方法分別采用中心插分格式和二階隱式格式.

1.2 ? 計算參數(shù)和結(jié)果驗(yàn)證

為驗(yàn)證本文采用的計算方法及計算參數(shù)的正確性，以單方柱為對象，分別研究了周向網(wǎng)格數(shù)量、無量綱時間步長Δt*（Δt*=ΔtU0 /B，其中Δt為有量綱時間步）及展向長度對結(jié)果的影響. 表1列出了平均阻力系數(shù)、脈動阻力系數(shù)、脈動升力系數(shù)和Strouhal數(shù)等計算結(jié)果，并與文獻(xiàn)的試驗(yàn)值[15-17]和數(shù)值模擬結(jié)果[18-19]進(jìn)行了對比.

由表1可知，本文計算結(jié)果與文獻(xiàn)數(shù)值模擬結(jié)果相近，所有結(jié)果都落在文獻(xiàn)值范圍內(nèi);與文獻(xiàn)試驗(yàn)值相比，平均阻力系數(shù)、脈動升力系數(shù)和Strouhal吻合較好，而本文Case C1和Case C2的脈動阻力系數(shù)略大，這可能是展向長度偏小造成. 從總體上看，本文計算的7種工況的結(jié)果相近，但不同展向長度下脈動阻力系數(shù)有偏差. 綜合考慮計算精度和計算效率，本文采用Case A2的計算參數(shù)進(jìn)一步分析串列雙方柱的繞流問題. 即如圖3所示，每個方柱的周向布置300個單元，在角部適當(dāng)加密;近壁面最小網(wǎng)格厚度為0.001B，近壁面的y+≈1;展向長度為4B，共劃分40層單元. 串列雙方柱的計算模型總單元數(shù)為320萬至420萬不等，無量綱時間步長Δt*=0.005.

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

2.1 ? 流場結(jié)構(gòu)

圖 4 給出了典型間距比（P/B=1.25、3、4）下串列雙方柱的平均流線圖（坐標(biāo)的長度單位為B），分別代表了本文計算得到的3種流場結(jié)構(gòu)，即單一鈍體（P/B=1.1～1.5）、剪切層再附（P/B=2～3.5）和雙渦脫（P/B=4～5）流態(tài)，這與文獻(xiàn)[13]的結(jié)果相似.

從圖4（a）可以看出：當(dāng)P/B=1.25時，上游方柱形成的剪切層將兩個方柱完整的包裹起來，并在柱間形成較小的對稱回流，兩方柱就像一個鈍體;隨著間距比增大（P/B=3），如圖4（b）所示，上游的剪切層會在下游方柱的側(cè)風(fēng)面形成再附，進(jìn)而尾流回流區(qū)也變窄;當(dāng)間距比繼續(xù)增大（P/B=4），如圖4（c）所示，上游方柱的剪切層不再形成再附，上游方柱接近單方柱流場形態(tài)，且兩個方柱的尾流處均形成對稱回流并伴隨旋渦脫落.

2.2 ? 氣動力系數(shù)

圖5為上、下游方柱的氣動力系數(shù)隨間距比的變化曲線. 由圖5（a）可知，本文平均阻力系數(shù)與文獻(xiàn)值變化趨勢一致，吻合情況良好. 當(dāng)P/B≤3.5時，兩方柱的平均阻力系數(shù)變化較平穩(wěn)，當(dāng)P/B為3.5～4時出現(xiàn)氣動力跳躍，這是由繞流場結(jié)構(gòu)從剪切層再附流態(tài)變?yōu)殡p渦脫流態(tài)造成的，發(fā)生該氣動力跳躍的間距比一般稱為臨界間距比. 值得注意的是，不同文獻(xiàn)中的臨界間距比有一定差異，可能受到來流湍流度、雷諾數(shù)等的影響.

圖5（b）和圖5（c）為串列雙方柱的脈動氣動力系數(shù)隨間距比的變化. 可以看出，兩個方柱的脈動氣動力值隨間距比變化的趨勢一致，同樣在臨界間距比處氣動力出現(xiàn)跳躍，且與文獻(xiàn)值吻合良好. 當(dāng)P/B小于臨界間距比時，脈動氣動力值均小于單方柱，發(fā)生氣動力跳躍后脈動阻力和脈動升力都大于單方柱，這種變化是伴隨串列雙方柱流態(tài)的轉(zhuǎn)變發(fā)生的. 此外，在雙渦脫流態(tài)下，上游方柱的脈動阻力系數(shù)小于下游方柱，而脈動升力系數(shù)卻相反. 還需注意的是，當(dāng)P/B≤2時，脈動升力系數(shù)隨間距比增大有逐漸減小的趨勢.

2.3 ? 表面風(fēng)壓系數(shù)

圖6為間距比P/B=1.25、3、4時雙方柱的平均風(fēng)壓系數(shù)分布圖，方柱表面第i個測點(diǎn)t時刻的風(fēng)壓系數(shù)時程Cp（i，t） = （式中：p（i，t）為i測點(diǎn)的總風(fēng)壓時程，p0為來流靜壓，ρ為空氣密度，U0為來流風(fēng)速）. 由圖6可知，3種不同間距的上、下游方柱平均風(fēng)壓系數(shù)分布與單方柱有所不同. 對于上游方柱，平均風(fēng)壓系數(shù)在迎風(fēng)面處的分布一致，但在背風(fēng)面和側(cè)風(fēng)面處存在兩種分布趨勢，當(dāng)P/B=4時上游方柱的平均風(fēng)壓系數(shù)分布與單方柱的趨勢一致. 對于下游方柱，由于受上游方柱的干擾作用，其平均風(fēng)壓分布與單方柱的差異更大，整體都呈現(xiàn)負(fù)壓;當(dāng)P/B小于臨界間距比時（P/B≤3.5）迎風(fēng)面處負(fù)壓最強(qiáng)，而P/B=4時背風(fēng)面處負(fù)壓最強(qiáng).

圖7為雙方柱的脈動風(fēng)壓系數(shù)分布圖，可以看出不同間距下的脈動風(fēng)壓分布有顯著差異. P/B小于臨界間距比時（P/B≤3.5），即在單一鈍體和剪切層再附流態(tài)下，上、下游方柱的風(fēng)壓脈動值較小，均顯著低于單方柱;而在雙渦脫流態(tài)下（P/B = 4～5），兩個方柱的風(fēng)壓脈動值明顯增大，上游方柱的脈動風(fēng)壓與單方柱相當(dāng)，而下游方柱的風(fēng)壓脈動值總體較單方柱大.

2.4 ? 風(fēng)壓非高斯特性

圖8～圖10為3種流態(tài)下（P/B=1.25、3、4）兩個方柱的風(fēng)壓偏度（S=（Cp（i，t）-Cp（i））3 σ-3

cp（i，t））和峰度（K = （Cp（i，t） - Cp（i））4 σ-4

cp（i，t））分布. 其中，σ

cp（i，t）為測點(diǎn)i的風(fēng)壓時程根方差;N為風(fēng)壓系數(shù)時程樣本數(shù)，本文的計算工況所采樣的風(fēng)壓系數(shù)時程樣本數(shù)為6×104;采樣起始時刻和結(jié)束時刻對應(yīng)的無量綱時間t*（t*=tU0 /B）分別為330和630. 考慮到方柱的對稱性，本文僅給出上半部分的分布情況. 整體來看，串列方柱的偏度和峰度均與單方柱有明顯差異，上游方柱的風(fēng)壓大部分呈負(fù)偏度，而下游方柱的峰度波動較為劇烈.

圖8為單一鈍體流態(tài)（P/B=1.25）時風(fēng)壓偏度和峰度的分布. 對于上游方柱，如圖8（a）所示，偏度和峰度在迎風(fēng)面和側(cè)風(fēng)面處都較穩(wěn)定，只在后角點(diǎn)處有明顯波動，而背風(fēng)面處存在較大的負(fù)偏度和峰度. 由圖8（b）可知，下游方柱的峰度在迎風(fēng)面和后角點(diǎn)處有明顯波動，并對應(yīng)較大的負(fù)偏度和峰度.

圖9為剪切層再附流態(tài)（P/B=3）時風(fēng)壓偏度和峰度的分布. 上游方柱在迎風(fēng)面和側(cè)風(fēng)面處的分布較平緩，如圖9（a）所示，僅在背風(fēng)面呈現(xiàn)較大的正偏度和峰度;對于下游方柱，由圖9（b）可知，峰度的分布在迎風(fēng)面和背風(fēng)面處均出現(xiàn)極值并對應(yīng)較大的負(fù)偏度.

圖10為雙渦脫流態(tài)（P/B=4）時的風(fēng)壓偏度和峰度的分布. 由圖10（a）可知，上游方柱偏度和峰度分布與單方柱接近;而下游方柱的分布趨勢與單方柱明顯不同，如圖10（b）所示，在側(cè)風(fēng)面和背風(fēng)面處均出現(xiàn)較大的負(fù)偏度和峰度.

參照文獻(xiàn)[11]對風(fēng)壓非高斯特性的劃分標(biāo)準(zhǔn)：定義偏度|S|>0.2且峰度|K|>3.5的區(qū)域?yàn)轱L(fēng)壓非高斯區(qū). 本文基于上述風(fēng)壓偏度和峰度的分布，來劃分串列雙方柱的風(fēng)壓非高斯區(qū)域.

2.5 ? 表面風(fēng)壓的相關(guān)性

圖11為3種流態(tài)下（P/B=1.25、3、4）上、下游方柱的表面風(fēng)壓的相關(guān)系數(shù)，測點(diǎn)j和測點(diǎn)k的風(fēng)壓相關(guān)系數(shù)R（Cp j，Cpk） = ，cov（Cp j，Cp k）為測點(diǎn)j和k風(fēng)壓時程的協(xié)方差. 考慮到方柱的對稱性，本文僅給出上游方柱典型測點(diǎn)與下游方柱所有測點(diǎn)的風(fēng)壓相關(guān)系數(shù). 結(jié)果表明，3種流態(tài)下風(fēng)壓相關(guān)性分布明顯不同，隨間距增大相關(guān)系數(shù)整體呈現(xiàn)減小又增大的趨勢. 由圖11（a）可知，對于單一鈍體流態(tài)（P/B=1.25），上游方柱背風(fēng)面與下游方柱迎風(fēng)面的風(fēng)壓相關(guān)性最強(qiáng)，最大值為0.9左右;此外，兩方柱側(cè)風(fēng)面的風(fēng)壓相關(guān)性也較強(qiáng). 相比之下，對于剪切層再附流態(tài)（P/B=3），如圖11（b）所示，上游方柱與下游方柱的風(fēng)壓相關(guān)性整體較差，相關(guān)性曲線分布在0附近. 而對于雙渦脫流態(tài)（P/B=4），如圖11（c）所示，兩方柱又表現(xiàn)出較強(qiáng)的風(fēng)壓相關(guān)性，上、下游方柱側(cè)風(fēng)面風(fēng)壓相關(guān)系數(shù)接近0.8.

3 ? 流場特性

為了進(jìn)一步探討風(fēng)壓特性與流場特性之間的內(nèi)在聯(lián)系，圖12～圖14給出了三種流態(tài)下的典型流場圖（平均風(fēng)壓場、瞬時渦量圖及瞬時流線圖），對下游方柱升力系數(shù)達(dá)到最大值時刻的流場特征進(jìn)行分析，進(jìn)而說明不同流態(tài)與風(fēng)壓特性之間的內(nèi)在聯(lián)系. 如上文所述，平均風(fēng)壓場中方柱側(cè)面的加粗實(shí)線代表風(fēng)壓呈現(xiàn)非高斯特性的分布區(qū)域.

3.1 ? 單一鈍體流態(tài)

圖12為單一鈍體流態(tài)的流場圖（P/B=1.25）. 由平均風(fēng)壓場可以看到，兩方柱側(cè)風(fēng)面及方柱間有較強(qiáng)負(fù)壓，柱間的負(fù)壓最強(qiáng)且呈現(xiàn)明顯的風(fēng)壓非高斯特性，風(fēng)壓非高斯區(qū)域主要集中在柱間及下游方柱后角點(diǎn)附近. 由瞬時流場可知，單一鈍體流態(tài)下對應(yīng)的瞬時渦量較大，兩個方柱之間存在變化的回流，受同一回流區(qū)影響，該回流區(qū)附近的風(fēng)壓相關(guān)性很強(qiáng)且呈現(xiàn)明顯的風(fēng)壓非高斯特性. 此外，受同一剪切層包裹作用，兩方柱側(cè)風(fēng)面也存在較強(qiáng)風(fēng)壓相關(guān)性.

3.2 ? 剪切層再附流態(tài)

圖13為剪切層再附流態(tài)的流場圖（P/B=3）.由平均風(fēng)壓場可知，兩方柱的負(fù)壓區(qū)向下游方柱前角點(diǎn)移動，下游方柱側(cè)風(fēng)面負(fù)壓減弱，而風(fēng)壓非高斯區(qū)域卻較大，兩方柱的背風(fēng)面及下游迎風(fēng)面前角處為主要風(fēng)壓非高斯區(qū)域. 由瞬時流場可知，剪切層再附流態(tài)下方柱尾流渦脫強(qiáng)度低，故而兩方柱風(fēng)壓相關(guān)性整體很弱;上游方柱的剪切層會撞擊下游方柱，并導(dǎo)致下游方柱迎風(fēng)面呈風(fēng)壓非高斯特性.

3.3 ? 雙渦脫流態(tài)

圖14為雙渦脫流態(tài)的流場圖（P/B=4）.由平均風(fēng)壓場可知，兩方柱背風(fēng)面均有較強(qiáng)的負(fù)壓;下游方柱的風(fēng)壓非高斯區(qū)域大于上游方柱，主要集中在背風(fēng)面及側(cè)風(fēng)面前角處. 由瞬時流場可知，雙渦脫流態(tài)下的上游方柱尾流形成強(qiáng)烈的旋渦脫落，脫落的旋渦會撞擊下游方柱前角點(diǎn)附近的迎風(fēng)面和側(cè)風(fēng)面，故兩方柱側(cè)風(fēng)面存在較強(qiáng)的風(fēng)壓相關(guān)性且下游方柱側(cè)風(fēng)面前角處呈顯著的風(fēng)壓非高斯特性.

4 ? 結(jié) ? 論

本文在雷諾數(shù)8 × 104的均勻來流條件下，對串列雙方柱繞流進(jìn)行了大渦模擬研究，分析了不同間距比（P/B=1.1～5）情況下的流態(tài)、風(fēng)壓變化，重點(diǎn)探討了風(fēng)壓的非高斯特性，并從流場結(jié)構(gòu)角度闡釋了非高斯特性的作用機(jī)理. 主要結(jié)論如下：

1）隨間距比增大，串列雙方柱流場依次表現(xiàn)為單一鈍體（P/B=1.1～1.5）、剪切層再附（P/B=2～3.5）和雙渦脫（P/B=4～5）3種流態(tài)，3種流態(tài)下的方柱的氣動性能、風(fēng)壓特性和繞流場特征有明顯差異.

2）對于單一鈍體流態(tài)（P/B=1.1～1.5），表面風(fēng)壓在柱間和下游方柱后角點(diǎn)附近呈非高斯特性. 受柱間回流的影響，回流區(qū)附近表面風(fēng)壓具有顯著的非高斯特性，且風(fēng)壓相關(guān)性很強(qiáng);此外，由于受同一剪切層包裹，雙方柱側(cè)風(fēng)面的風(fēng)壓相關(guān)性也較強(qiáng).

3）對于剪切層再附流態(tài)（P/B=2～3.5），風(fēng)壓非高斯區(qū)域主要為雙方柱的背風(fēng)面和下游方柱迎風(fēng)面前角處. 由于渦脫強(qiáng)度較低，風(fēng)壓相關(guān)性整體差;上游方柱的剪切層撞擊下游方柱前角點(diǎn)，導(dǎo)致下游方柱迎風(fēng)面呈現(xiàn)大范圍的風(fēng)壓非高斯區(qū)域.

4）對于雙渦脫流態(tài)（P/B=4～5），雙方柱背風(fēng)面和側(cè)風(fēng)面近角點(diǎn)處為主要風(fēng)壓非高斯區(qū)域. 受上游方柱尾流渦脫的撞擊作用，兩方柱側(cè)風(fēng)面的風(fēng)壓相關(guān)性較強(qiáng)，下游方柱側(cè)風(fēng)面前角處呈顯著的風(fēng)壓非高斯特性.

參考文獻(xiàn)

[1] ? ?AMIN J A， AHUJA A K. Characteristics of wind forces and responses of rectangular tall buildings[J]. International Journal of Advanced Structural Engineering（IJASE），2014，6（3）：1—14.

[2] ? ?柯世堂，徐璐，朱鵬. 基于大渦模擬超大型冷卻塔施工期風(fēng)荷載時域特性分析[J]. 湖南大學(xué)學(xué)報（自然科學(xué)版），2018，45（11）：62—72.

KE S T，XU L，ZHU P. Analysis on time domain characteristics of wind loads during construction period of super large cooling towers based on large eddy simulation[J]. Journal of Hunan University （Natural Sciences），2018，45（11）：62—72. （In Chinese）

[3] ? ?HUI Y，TAMURA Y，YOSHIDA A，et al. Pressure and flow field investigation of interference effects on external pressures between high-rise buildings[J]. Journal of Wind Engineering and Industrial Aerodynamics，2013，115：150—161.

[4] ? ?SADEK F，SIMIU E. Peak non-Gaussian wind effects for database-assisted low-rise building design[J]. Journal of Engineering Mechanics，2002，128（5）：530—539.

[5] ? ?李正農(nóng)，余蜜，吳紅華，等. 某低矮模型房屋實(shí)測風(fēng)場和風(fēng)壓的相關(guān)性研究[J]. 湖南大學(xué)學(xué)報（自然科學(xué)版），2016，43（5）：70—78.

LI Z N，YU M，WU H H，et al. Correlation research of the measured wind field and wind pressure of a low-rise building[J]. Journal of Hunan University （Natural Sciences），2016，43（5）：70—78. （In Chinese）

[6] ? ?SURESH KUMAR K，STATHOPOULOS T. Fatigue analysis of roof cladding under simulated wind loading[J].Journal of Wind Engineering and Industrial Aerodynamics，1998，77/78：171—183.

[7] ? ?KUMAR K S，STATHOPOULOS T. Synthesis of non-Gaussian wind pressure time series on low building roofs[J]. Engineering Structures，1999，21（12）：1086—1100.

[8] ? ?KO N H，YOU K P，KIM Y M. The effect of non-Gaussian local wind pressures on a side face of a square building[J]. Journal of Wind Engineering & Industrial Aerodynamics，2005，93（5）：383—397.

[9] ? ?韓寧，顧明. 方形高層建筑風(fēng)壓脈動非高斯特性分析[J]. 同濟(jì)大學(xué)學(xué)報（自然科學(xué)版），2012，40（7）：971—976.

HAN N，GU M.Analysis on non-Gaussian features of fluctuating wind pressure on square tall buildings[J]. Journal of Tongji University （Natural Science），2012，40（7）：971—976. （In Chinese）

[10] ?莊翔，董欣，丁潔民，等. 矩形高層建筑表面風(fēng)壓脈動的非高斯特性研究[J]. 建筑結(jié)構(gòu)學(xué)報，2016，37（S1）：13—18.

ZHUANG X，DONG X，DING J M，et al. Non-Gaussian features of wind pressure fluctuations on a rectangular high-rise building[J]. Journal of Building Structures，2016，37（S1）：13—18.（In Chinese）

[11] ?樓文娟，李進(jìn)曉，沈國輝，等. 超高層建筑脈動風(fēng)壓的非高斯特性[J]. 浙江大學(xué)學(xué)報（工學(xué)版），2011，45（4）：671—677.

LOU W J，LI J X，SHEN G H，et al. Non-Gaussian feature of wind-induced pressure on super-tall building[J].Journal of Zhejiang University （Engineering Science），2011，45（4）：671—677.（In Chinese）

[12] ?王浩，柯世堂. 不同四塔組合形式對特大型冷卻塔局部非高斯風(fēng)壓分布影響研究[J]. 工程力學(xué)，2018，35（8）：162—171.

WANG H，KE S T. Research on non-Gaussian wind pressure of four super-large cooling towers under different layouts[J]. Engineering Mechanics，2018，35（8）：162—171. （In Chinese）

[13] ?SOHANKAR A. A numerical investigation of the flow over a pair of identical square cylinders in a tandem arrangement[J]. International Journal for Numerical Methods in Fluids，2012，70（10）：1244—1257.

[14] ?YEN S C，SAN K C，CHUANG T H. Interactions of tandem square cylinders at low Reynolds numbers[J]. Experimental Thermal & Fluid Science，2008，32（4）：927—938.

[15] ?杜曉慶，許漢林，馬文勇，等. 串列雙方柱氣動干擾效應(yīng)試驗(yàn)研究[J]. 建筑結(jié)構(gòu)學(xué)報，2019，40（11）：27—34.

DU X Q，XU H L，MA W Y，et al. Experimental study on aerodynamic interference effect of two square cylinders in a tandem arrangement[J]. Journal of Building Structures，2019，40（11）：27—34. （In Chinese）

[16] ?LEE B E. The effect of turbulence on the surface pressure field of a square prism[J]. Journal of Fluid Mechanics，1975，69（2）：263—282.

[17] ?VICKERY B J. Fluctuating lift and drag on a long cylinder of square cross-section in a smooth and in a turbulent stream[J].Journal of Fluid Mechanics，1966，25（3）：481—494.

[18] ?SOHANKAR A，DAVIDSON L，NORBERG C. Large eddy simulation of flow past a square cylinder：comparison of different subgrid scale models[J]. Journal of Fluids Engineering，2000，122（1）：39—47.

[19] ?RODI W. Comparison of LES and RANS calculations of the flow around bluff bodies[J]. Journal of Wind Engineering and Industrial Aerodynamics，1997，69/70/71：55—75.

[20] ?MAHBUB ALAM M D，MORIYA M，TAKAI K，et al. Suppression of fluid forces acting on two square prisms in a tandem arrangement by passive control of flow[J]. Journal of Fluids and Structures，2002，16（8）：1073—1092.

[21] ?SOHANKAR A. A LES study of the flow interference between tandem square cylinder pairs[J]. Theoretical and Computational Fluid Dynamics，2014，28（5）：531—548.

[22] ?LIU C H，CHEN J M. Observations of hysteresis in flow around two square cylinders in a tandem arrangement[J]. Journal of Wind Engineering and Industrial Aerodynamics，2002，90（9）：1019—1050.