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

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

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

中圖分類號：TU 311 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼：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

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

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

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

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

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

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

圖1為計算模型示意圖，來流為均勻來流，基于來流風速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為亞格子應力張量. 采用WALE亞格子尺度模型， WALE常數(shù)為0.325. 采用SIMPLEC算法求解壓力-速度場耦合方程，空間離散方法和時間離散方法分別采用中心插分格式和二階隱式格式.

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

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

由表1可知，本文計算結(jié)果與文獻數(shù)值模擬結(jié)果相近，所有結(jié)果都落在文獻值范圍內(nèi);與文獻試驗值相比，平均阻力系數(shù)、脈動升力系數(shù)和Strouhal吻合較好，而本文Case C1和Case C2的脈動阻力系數(shù)略大，這可能是展向長度偏小造成. 從總體上看，本文計算的7種工況的結(jié)果相近，但不同展向長度下脈動阻力系數(shù)有偏差. 綜合考慮計算精度和計算效率，本文采用Case A2的計算參數(shù)進一步分析串列雙方柱的繞流問題. 即如圖3所示，每個方柱的周向布置300個單元，在角部適當加密;近壁面最小網(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）下串列雙方柱的平均流線圖（坐標的長度單位為B），分別代表了本文計算得到的3種流場結(jié)構(gòu)，即單一鈍體（P/B=1.1～1.5）、剪切層再附（P/B=2～3.5）和雙渦脫（P/B=4～5）流態(tài)，這與文獻[13]的結(jié)果相似.

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

2.2 ? 氣動力系數(shù)

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

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

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

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

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

2.4 ? 風壓非高斯特性

圖8～圖10為3種流態(tài)下（P/B=1.25、3、4）兩個方柱的風壓偏度（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）為測點i的風壓時程根方差;N為風壓系數(shù)時程樣本數(shù)，本文的計算工況所采樣的風壓系數(shù)時程樣本數(shù)為6×104;采樣起始時刻和結(jié)束時刻對應的無量綱時間t*（t*=tU0 /B）分別為330和630. 考慮到方柱的對稱性，本文僅給出上半部分的分布情況. 整體來看，串列方柱的偏度和峰度均與單方柱有明顯差異，上游方柱的風壓大部分呈負偏度，而下游方柱的峰度波動較為劇烈.

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

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

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

參照文獻[11]對風壓非高斯特性的劃分標準：定義偏度|S|>0.2且峰度|K|>3.5的區(qū)域為風壓非高斯區(qū). 本文基于上述風壓偏度和峰度的分布，來劃分串列雙方柱的風壓非高斯區(qū)域.

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

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

3 ? 流場特性

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

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

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

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

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

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

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

4 ? 結(jié) ? 論

本文在雷諾數(shù)8 × 104的均勻來流條件下，對串列雙方柱繞流進行了大渦模擬研究，分析了不同間距比（P/B=1.1～5）情況下的流態(tài)、風壓變化，重點探討了風壓的非高斯特性，并從流場結(jié)構(gòu)角度闡釋了非高斯特性的作用機理. 主要結(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)下的方柱的氣動性能、風壓特性和繞流場特征有明顯差異.

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

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

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

參考文獻

[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] ? ?柯世堂，徐璐，朱鵬. 基于大渦模擬超大型冷卻塔施工期風荷載時域特性分析[J]. 湖南大學學報（自然科學版），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)，余蜜，吳紅華，等. 某低矮模型房屋實測風場和風壓的相關(guān)性研究[J]. 湖南大學學報（自然科學版），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] ? ?韓寧，顧明. 方形高層建筑風壓脈動非高斯特性分析[J]. 同濟大學學報（自然科學版），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] ?莊翔，董欣，丁潔民，等. 矩形高層建筑表面風壓脈動的非高斯特性研究[J]. 建筑結(jié)構(gòu)學報，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]. 浙江大學學報（工學版），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] ?王浩，柯世堂. 不同四塔組合形式對特大型冷卻塔局部非高斯風壓分布影響研究[J]. 工程力學，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] ?杜曉慶，許漢林，馬文勇，等. 串列雙方柱氣動干擾效應試驗研究[J]. 建筑結(jié)構(gòu)學報，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.