張玥 劉子琦 石慧慧

摘?要：為得到山區(qū)峽谷典型Ⅴ型地形下的橋梁抗風(fēng)設(shè)計的關(guān)鍵風(fēng)特性參數(shù)，通過“數(shù)值風(fēng)洞試驗”建立模型，以實橋為工程背景構(gòu)建三維地形網(wǎng)格，基于3種典型湍流模型和流體壁面粗糙度對橋位的風(fēng)場特性進行分析，選取最優(yōu)湍流模型和流體壁面粗糙度對該地形橋位處的風(fēng)場特性進行分析；圍繞兩個影響橋址區(qū)風(fēng)特性的地形參數(shù)（山體高度和夾角），闡明Ⅴ型峽谷地形下風(fēng)特性參數(shù)的變化規(guī)律，并與現(xiàn)行規(guī)范進行對比，推算橋面設(shè)計基準風(fēng)速。結(jié)果表明：該Ⅴ型峽谷地形的最佳湍流模型和流體壁面粗糙度分別為RNG ?k-ε和20 m；隨著山體高度及山體夾角兩參數(shù)的變化，風(fēng)速增加幅度分別為2.37%～12.56%和1.24%～6.98%；山體高度及山體夾角兩變參數(shù)下，梯度風(fēng)高度分別為700 m和800 m左右；實際湍流強度大于規(guī)范中四類地表規(guī)范值，在離地高度40 m范圍內(nèi)，湍流強度更接近于D類地表。實際Ⅴ型峽谷橋址區(qū)不能簡單按規(guī)范歸為C類或D類地形，在設(shè)計風(fēng)參數(shù)計算時應(yīng)綜合考慮規(guī)范和經(jīng)驗公式，還可借助模擬手段來提高其準確性。

關(guān)鍵詞：風(fēng)特性；Ⅴ型峽谷；數(shù)值風(fēng)洞試驗；湍流模型；地形參數(shù)；設(shè)計基準風(fēng)速

中圖分類號：U 442

文獻標志碼：A

文章編號：1672-9315（2024）03-0575-12

DOI：10.13800/j.cnki.xakjdxxb.2024.0317開放科學(xué)（資源服務(wù)）標識碼（OSID）：

Wind characteristics of bridge site under typical canyon

V-shaped landform

ZHANG Yue，LIU Ziqi，SHI Huihui

（College of Architecture and Civil Engineering，Xian University of Science and Technology，Xian 710054，China）

Abstract：In order to obtain the key wind characteristic parameters for wind resistance design of Bridges in typical V-shaped terrain in mountain canyons，an ?model was established by the method of “numerical wind tunnel test”，a three-dimensional terrain grid was constructed with real bridges as the engineering object，and the wind field characteristics of bridge sites were analyzed based on three typical turbulence models and fluid wall roughness.The optimal turbulence model and fluid wall roughness were selected to explore the wind field characteristics at the bridge site.From these two topographic parameters，mountain height and angle，which affect the wind characteristics of the bridge site，the variation law of the downwind characteristic parameters of the V-shaped canyon terrain is summarized，and the design reference wind speed of the bridge deck is calculated by comparing with the current norms.The results show that the optimal turbulence model and fluid wall roughness of the V-shaped canyon terrain are RNG and 20m respectively.With the change of mountain height and mountain angle，the increase of wind speed is 2.37%～12.56% and 1.24%～6.98%，respectively.Under the two variable parameters of mountain height and mountain angle，the gradient wind height is about 700 m ?and 800 m respectively.The actual turbulence intensity is greater than the standard values of the four types of surface in the specification，and the turbulence intensity is closer to the D type surface in the range of 40 m above the ground.The actual V-shaped canyon bridge site can not be classified as C or D terrain

simply according to the code.The code and empirical formula must be considered comprehensively when designing the wind parameter calculation，and the accuracy could be improved by means of simulation.

Key words：wind characteristics；V-shaped canyon；numerical wind tunnel test；turbulence model；terrain parameters；basic design wind speed

0?引?言

隨著經(jīng)濟和科技的快速發(fā)展，在西部山區(qū)建造的橋梁日益增多，當風(fēng)遇到障礙物山體和峽谷時，會出現(xiàn)峽谷風(fēng)的放大效應(yīng)和山體繞流、附著效應(yīng)，尤其在復(fù)雜山區(qū)地形，山體對橋址區(qū)的風(fēng)場分布規(guī)律和流場作用機理的影響極大。為了保證橋梁的正常運營和后期維護，需探究山區(qū)典型地形橋址區(qū)風(fēng)場特性的分布規(guī)律和典型風(fēng)向角影響下的分布機理[1-3]。

目前，由于現(xiàn)行抗風(fēng)規(guī)范中對地表的分類較粗略，并未詳細劃分典型的山區(qū)峽谷類型而是統(tǒng)一歸為D類，導(dǎo)致在具體的典型峽谷地貌中，規(guī)范確定值與數(shù)值風(fēng)洞模擬所得結(jié)果產(chǎn)生差異。周繼等以湘江特大橋為例研究復(fù)雜地形橋位下的風(fēng)特性，發(fā)現(xiàn)經(jīng)過CFD數(shù)值模擬分析的主梁設(shè)計基準風(fēng)速與抗風(fēng)規(guī)范確定的基準風(fēng)速相比偏大[4]。所以對于山區(qū)復(fù)雜地形尤其是橫跨兩山峽谷地貌下的橋梁結(jié)構(gòu)，其抗風(fēng)設(shè)計所需風(fēng)荷載，仍需借助物理風(fēng)洞試驗和現(xiàn)場實測等方法進行相關(guān)風(fēng)參數(shù)的修正。針對上述問題，當風(fēng)洞試驗或現(xiàn)場實測的數(shù)據(jù)支撐難以得到時，數(shù)值模擬快速準確地確定設(shè)計所需風(fēng)特性參數(shù)的優(yōu)勢就體現(xiàn)出來了。借助“數(shù)值風(fēng)洞”保證抗風(fēng)設(shè)計的準確性和可靠性就成了這一方法的關(guān)鍵，一般來說計算域中湍流模型的選取很大程度會導(dǎo)致計算結(jié)果的準確性。目前針對湍流模型的研究很多，

張向旭、李永樂、于艦涵等選取SST k-ε模型對復(fù)雜山體地形進行風(fēng)場模擬計算分析[5-7]；薛亞飛、孫亭亭、王鐵強等采用了RNG k-ε模型對山區(qū)橋位處風(fēng)場進行了模擬研究[8-10]。其中，王鐵強基于6種湍流模型進行了分析，得出Standard k-ω模型和SST ?k-ω模型用于高雷諾數(shù)流場的近地面風(fēng)場計算時會出現(xiàn)較大的誤差；而RNG ?k-ε模型和Realizable ??k-ε、RSM模型的模擬結(jié)果在總體誤差比較相近且均較??；Realizable ??k-ε模型在主要計算旋轉(zhuǎn)和靜態(tài)流動區(qū)域時不能提供自然的湍流黏度；RSM用于復(fù)雜流動時有潛力達到更高的精度，但計算量大。最終可知RNG ???k-ε模型較適合在風(fēng)場模擬中使用[10]；于濤，唐昂等[11-18]采用了Realizable ??k-ε模型對山區(qū)風(fēng)場進行了模擬運算，其中，馮林以黃賢村山體風(fēng)場為研究對象，選取了3種典型湍流模型與實測數(shù)據(jù)比較，得出RNG ??k-ε模型結(jié)果誤差最大，平均誤差約達5%，而Realizable ??k-ε模型的結(jié)果誤差最小，誤差值約為3%，Realizable ??k-ε模型模擬精度較高；在地形參數(shù)方面，樓文娟、陳平、李正昊等通過數(shù)值模擬和風(fēng)洞試驗較全面的研究了山體之間的間距、山體形狀以及峽谷長度等地形參數(shù)對峽谷內(nèi)部風(fēng)場的影響[19-21]；BOWEN、CAO、LUBITZ等分別研究了山體坡度、形狀和來流風(fēng)向角等流場影響因素對山地風(fēng)場的影響[22-24]。

針對峽谷地形參數(shù)對風(fēng)場影響的研究多以理想二維模型或由單峰山丘輪廓模型組成的三維模型為研究對象，二維模型無法全面地考慮實際山體情況，而單峰山丘組成的三維模型無法在其基礎(chǔ)上改變山體高度和山體夾角來進行研究?；诖耍缘湫廷跣蛵{谷實際地形為研究對象，以矮寨大橋為工程背景，對山體地形參數(shù)改變時的風(fēng)場變化規(guī)律進行分析，并提出梯度風(fēng)高度范圍和計算出該峽谷地形下橋址區(qū)的基準風(fēng)速，以及這種地形下適應(yīng)的地表規(guī)范歸屬類型，為后續(xù)相似地形模型的建立及抗風(fēng)參數(shù)的研究提供參考。

1?Ⅴ型谷三維地形模型

1.1?實際地形模型

矮寨大橋位于湖南湘西州吉首市矮寨鎮(zhèn)，海拔高度680～900 m，跨中橋面標高與谷底高差達340 m。根據(jù)三維數(shù)字地形圖LSV與下載軟件定位矮寨大橋并選取計算域，橋位坐標（N28°19′58.48″，E109°35′53.75″）?？紤]到地形邊界對風(fēng)場回流的影響和保證流場的充分發(fā)展，以及為了更接近真實的風(fēng)場，確定將計算域選為2 200 m×2 200 m×2 000 m，通過提取該地形高程數(shù)據(jù)點，將高程點換算為經(jīng)緯度坐標后，利用萬能坐標轉(zhuǎn)換器轉(zhuǎn)換為直角坐標系即XYZ形式，導(dǎo)入Rhinoceros生成三維地形曲面，再通過Fluent前處理軟件Gambit對犀牛軟件Rhinoceros導(dǎo)出的sat文件進行建模，利用Gambit生成計算模型。矮寨大橋的橋位地形如圖1所示，由犀牛軟件Rhinoceros生成的地形曲面如圖2所示。

1.2?網(wǎng)格劃分以及邊界條件劃分

運用Gambit進行非結(jié)構(gòu)性網(wǎng)格劃分、設(shè)定邊界條件，采用三角形非結(jié)構(gòu)化網(wǎng)格，近地面最小網(wǎng)格尺度30 mm。速度入口邊界參數(shù)時一般采用均勻流或?qū)L(fēng)速設(shè)置成風(fēng)速剖面計算，以此模擬實際峽谷地形中風(fēng)速沿高度變化而變化的特性，試驗?zāi)M采用均勻流。數(shù)值模擬計算僅針對最不利均勻流場進行分析，由唐金旺針對地形做的風(fēng)洞試驗[14]可知均勻流場的最不利風(fēng)速為10 m/s，故入口風(fēng)速值設(shè)置為10 m/s。其中，針對計算域邊界條件和網(wǎng)格劃分如圖3～圖5所示。

1.3?湍流模型及計算工況

為了確保此類地形風(fēng)場特性計算的準確性，選取兩方程帶旋修正Realizable ?k-ε、重整化群RNG ?k-ε和剪應(yīng)力傳輸SST ?k-ε 3種湍流模型進行對比分析，湍流模型求解過程中的各項具體參數(shù)設(shè)置見表1。童祎巍對矮寨大橋安裝設(shè)置觀測點進行現(xiàn)場實測[25]，由實測所統(tǒng)計得出的強風(fēng)樣本可知矮寨大橋所處峽谷的大風(fēng)大都來自偏北風(fēng)，其次是偏西方向和偏東方向?？紤]矮寨大橋?qū)崪y強風(fēng)風(fēng)向及橋址處峽谷兩側(cè)山體對風(fēng)的遮擋效應(yīng)，按照標準逆時針設(shè)置以下3種典型風(fēng)向角作為計算工況進行數(shù)值模擬研究：風(fēng)向角為0°時，

記為工況1，風(fēng)向角為90°時記為工況2，風(fēng)向角

為270°時記為工況3，3種工況均采用圖4中設(shè)定的計算域邊界條件。為了分析橋位處風(fēng)速和湍流強度沿高程的變化，擬選取矮寨大橋上的5個特定位置作為測點研究其豎方向上的風(fēng)速和湍流強度大小。測點0～4分別對應(yīng)跨中點、吉首岸橋跨四分點、吉首岸橋塔位置、茶洞岸跨四分點、茶洞岸橋塔位置。3種工況的風(fēng)向角與測點布置如圖6、圖7所示。

2?實際橋址區(qū)的風(fēng)特性分析

2.1?風(fēng)速和湍流強度通過提?、跣蛵{谷3種湍流模型下3種工況和5個橋位測點風(fēng)速的計算結(jié)果，分析風(fēng)參數(shù)隨著離地高度變化的分布規(guī)律。橋梁跨中0號測點風(fēng)速和湍流強度高程剖面分布如圖8、圖9所示。根據(jù)峽谷地形下橋址區(qū)各工況典型位置處的豎向風(fēng)速剖面和湍流強度分析得出，各湍流模型的變化趨勢基本一致。跨中0號測點，在工況1和工況2的情況下，變化曲線走向一致但波動性較大，由于橋面兩側(cè)復(fù)雜山體的干擾作用，高程風(fēng)速剖面分布呈Z型變化，在300～500 m范圍出現(xiàn)兩次拐點：340～500 m高程段內(nèi)風(fēng)速隨高程增加而減小，高程達到500 m之后又隨高度增加而增加，其中風(fēng)速減小段可能是由橋面與兩側(cè)山體的相互作用影響，高程達到500 m之后風(fēng)速大致符合規(guī)范中冪指數(shù)的增長規(guī)律；工況3在高程500 m以下，風(fēng)剖面斜率較大，風(fēng)速增長速度明顯大于工況1和工況2?？赡茉蚴牵汗r3的情況下山體的遮擋效應(yīng)被削弱和峽谷兩側(cè)山體的擠壓導(dǎo)致風(fēng)速迅速增大，最大風(fēng)速與最小風(fēng)速的差值達到近0.93，梯度風(fēng)高度為700 m；橋梁跨中位置的最大風(fēng)速出現(xiàn)在RNG湍流模型上，達到10.2 m/s。

1號、2號和3號測點的風(fēng)速和湍流強度高程剖面如圖10、圖11所示。1號測點和2號測點分別為茶洞岸和吉首岸四分之一跨處，3種工況下的風(fēng)速剖面基本符合常規(guī)峽谷地形橋位風(fēng)剖面的特征，湍流模型風(fēng)速基本隨高度增加而增加，在高度達到700～800 m時風(fēng)速趨于不變。3號測點位于茶洞岸塔底山體的山頂處，此處有明顯的山頂風(fēng)速增大效應(yīng)，在600～640 m迅速增大，最大風(fēng)速為工況2下達到了11.8 m/s，比入口平均風(fēng)速增大了18%，之后趨于平緩。

3種典型湍流模型下各監(jiān)測點的湍流強度隨高程遞減的變化趨勢基本一致，其中尤以跨中0號測點下工況1和工況2兩種工況更為接近，變化趨勢較為復(fù)雜，在300～700 m發(fā)生了3次突轉(zhuǎn)，其他測點基本在高程范圍內(nèi)湍流強度隨高度增加而減少，且在700 m后趨于平緩，這是由不同風(fēng)向角下山體遮擋效應(yīng)引起的；各測點在高程方向下湍流強度最大與最小差值分別為：0.808、0.882、0.833、0.756；橋面處湍流強度削弱基本一致，以0號測點為例：Realizable ?k-ε模型下3種工況橋面高度處湍流強度與高程方向下最大湍流強度相比，分別降低了12.8%、14.5%和8.7%，RNG k-ε模型下分別降低了12.0%、10.7%和5.6%，SST ??k-ε模型下分別降低了19.4%、16.9%和9.3%；橋塔處有增強，最大達到了12.6%。由以上分析矮寨大橋Ⅴ型峽谷梯度風(fēng)高度大致為700 m左右，而抗風(fēng)規(guī)范中矮寨大橋規(guī)定為D類地形，由規(guī)范查得梯度風(fēng)高度為450 m，這與模擬所得結(jié)果有很大差異。這一差異說明規(guī)范中對地形粗略分類得出的風(fēng)特性與實際山區(qū)橋址區(qū)風(fēng)特性有所不同，在實際工程中應(yīng)借助現(xiàn)場實測及風(fēng)洞試驗、模擬等手段分析具體風(fēng)場特性。

2.2?最優(yōu)數(shù)值模型選取

通過對Ⅴ型峽谷不同方向角和典型橋址位置處各湍流模型風(fēng)特性隨高程變化情況的分析，并借助已有學(xué)者的風(fēng)洞試驗結(jié)果[26]，并設(shè)置不同粗糙高度（0，10，20 m）選出適合分析本地形風(fēng)場特性的最優(yōu)湍流模型。

2.2.1?流體壁面粗糙度

流體壁面粗糙度是Fluent對地形曲面計算域WALL的一種邊界模擬條件，刻畫了山體表面的黏滯阻力。通過改變他可以得出該底面邊界參數(shù)對流場風(fēng)特性的影響，以此達到模擬不同復(fù)雜地形下相適應(yīng)的計算情況。研究所選的橋址區(qū)為起伏較大的丘陵地區(qū)，粗糙高度選擇1 m，在利用Fluent設(shè)置壁面問題時應(yīng)基于山體地表植被的覆蓋情況，設(shè)置相應(yīng)的粗糙度高度h，為了選出最適宜的粗糙度高度，此處選取0，10和20 m進行研究。

2.2.2?模擬結(jié)果與試驗結(jié)果的比較

為得到最佳湍流模型，首先選取典型橋位跨中0號測點的平均風(fēng)速與丁冬的風(fēng)洞試驗結(jié)果[26]進行比較，該風(fēng)洞試驗使用2塊主體模型和4塊緩和導(dǎo)風(fēng)板設(shè)計出橋位地形對稱的Ⅴ型峽谷地形模型，試驗?zāi)Ｐ偷目s尺比是1∶1 000，結(jié)果見表2。其次選擇順峽谷方向工況1下各測點風(fēng)速計算結(jié)果進行對比，最后將該工況下不同壁面粗糙度下的湍流模型風(fēng)速進行對比，分析數(shù)值模擬與風(fēng)洞試驗的差異，以此選取最適合本地形的湍流模型，工況1下各測點風(fēng)速對比見表3。

從表2、表3可以看出，對于Ⅴ型峽谷來說，3種湍流模型計算出的風(fēng)速差別不大，并且和風(fēng)洞試驗數(shù)據(jù)相差不大。在用風(fēng)洞試驗進行結(jié)果驗證時，風(fēng)速相對誤差以不超過10%為宜，以上3種湍流模型都可適用，但RNG ?k-ε湍流模型計算出來的風(fēng)速與風(fēng)洞試驗數(shù)據(jù)最為接近，故最終選擇RNG ?k-ε湍流模型作為最佳方案。

將順峽谷方向下典型測點處3種不同的壁面粗糙度的RNG ?k-ε湍流模型風(fēng)速與風(fēng)洞試驗結(jié)果進行對比如圖12所示。3種粗糙度下風(fēng)速與風(fēng)洞試驗更為相近的是20 m，故選擇該粗糙高度下的RNG ?k-ε湍流模型作為最佳湍流模型計算。

3?峽谷地形參數(shù)的影響

為研究地形參數(shù)對Ⅴ型峽谷地形橋位風(fēng)特性的影響，在原Ⅴ型峽谷地形的基礎(chǔ)上，通過改變離散點高程的縱坐標和經(jīng)度來改變?nèi)S實際山體地形模型的高度和夾角，基于矮寨大橋數(shù)值模擬優(yōu)選出來的湍流模型和粗糙度高度對峽谷周圍風(fēng)場特性開展研究。該三維山體地形模型的測點布置如圖13所示，模擬結(jié)果選取典型的谷底0號測點和谷頂1號測點進行分析。

3.1?山體高度

由于在原地形山體高度的基礎(chǔ)上，200 m以內(nèi)風(fēng)速的變化值比較明顯，故建立山體高度增高120，80，40 m與減小40，80，120 m 6種不同山體的三維模型進行分析，將山體隨高度變化的工況分別列為H1、H2、H3、H4、H5和H6。對Ⅴ型谷6種不同山體高度在風(fēng)向角為0°和90°兩種工況下各測點至計算域頂部風(fēng)速的計算結(jié)果進行分析，提取0號測點與1號測點的山體風(fēng)特性分布規(guī)律，以便于分析風(fēng)速和湍流強度在不同山體高度的情況下隨高度變化的分布規(guī)律。由此初步得出風(fēng)速和湍流強度在不同山體高度各工況（風(fēng)向角）下隨高度的變化規(guī)律，如圖14～圖17所示。

從圖14、圖15可以看出，山谷谷底的近地面風(fēng)速隨著高度的增加會有一個過渡減緩期，減緩期高度在350～400 m，其影響程度可能受到峽谷兩側(cè)山體高度的增加而有所削弱，正是峽谷風(fēng)緊縮升高的一個放大效應(yīng)；山頂風(fēng)速則會隨著山體高度的增加有明顯的增強效應(yīng)；梯度風(fēng)高度大約在700 m左右，在該高度處，山體高度越高，風(fēng)速增強效果越大，最大達到12.56%。由規(guī)范可知在C類、D類地形中，梯度風(fēng)高度為400 m和450 m，這與本節(jié)計算結(jié)果差距較大，說明實際工程中山體不能按規(guī)范簡單歸為4類。

從圖16、圖17可以看出，山谷的湍流強度在近地面400～500 m有一個停滯期，整體上看其強度會隨著山體高度的增加有所下降；相反山頂?shù)耐牧鲝姸葎t隨山體高度的增加有所增強。

3.2?山體夾角為分析山體不同角度下橋址處的風(fēng)特性變化規(guī)律，建立山體夾角分別為30°、60°、90°和120°的4種山體模型，分別由θ1～4表示，如圖18所示。

提?、跣凸?種不同角度在風(fēng)向角為0°和90°兩種工況下谷底和谷頂至計算域頂部風(fēng)速的計算結(jié)果，分析風(fēng)速和湍流強度在不同山體角度下隨山體高度變化的規(guī)律，如圖19～圖22所示。

從圖19、圖20可以看出，兩個測點在風(fēng)向角為0°和90°2種工況下的不同山體夾角風(fēng)速隨高程變化規(guī)律基本一致，都隨著離地高度的增加而增加；山谷和山頂測點的風(fēng)速都隨著山體夾角的減小而增大，即山體越陡峭，風(fēng)速值越大；梯度風(fēng)高度大約在800 m左右，在該高度處，山體夾角越小，山體越陡峭，風(fēng)速增強效果越大，最大達到了6.98%。由規(guī)范可知在C類、D類地形中，梯度風(fēng)高度為400 m和450 m，這與分析結(jié)果差距較大，說明實際工程中山體不能按規(guī)范簡單歸為4類。

從圖21、圖22可以看出，山谷與山頂湍流強度均隨著夾角變大，山體越平緩，其湍流強度越小，且山谷的湍流強度有一個200 m的過渡期，這種過渡期隨著山體的陡峭程度越大越明顯，可能原因是受到兩側(cè)山體的遮擋擠壓影響。

3.3?Ⅴ型峽谷湍流強度與現(xiàn)行抗風(fēng)規(guī)范的對比

通過矮寨大橋在0號測點工況一的數(shù)值計算結(jié)果與抗風(fēng)規(guī)范離地高度160 m以下設(shè)計湍流強度的對比，如圖23所示，峽谷地形下的數(shù)值湍流強度較抗風(fēng)規(guī)范中4類地表有明顯的差異，在40 m以下，更接近于D類地表規(guī)范，在20 m高度處兩值最為接近，僅相差2.59%；而在40 m以上高度，峽谷地形的湍流強度較規(guī)范中4類地形的差異均較大，最大在離地150 m處相差32.06%。這些結(jié)果進一步說明不同峽谷地形下地形參數(shù)變化值幅度大，在工程設(shè)計中應(yīng)進一步借鑒實測數(shù)據(jù)、風(fēng)洞試驗或數(shù)值風(fēng)洞的結(jié)果，同時結(jié)合規(guī)范綜合考量，確保選取的湍流強度是正確的。

綜合分析可知，隨山體高度的增高，山谷處湍流強度呈現(xiàn)下降趨勢，山頂處湍流強度隨高度增加增強；隨山體夾角變大，山谷和山頂處湍流強度下降。受山體兩側(cè)高度和夾角的影響，山體高度變化對湍流強度的影響效應(yīng)明顯大于山體夾角變化的影響。

3.4?山區(qū)峽谷橋面設(shè)計基準風(fēng)速的推算與比較

依據(jù)橋梁地理位置，選取距離橋址較近的吉首市氣象站作為基本風(fēng)速參考站點，查《公路橋梁抗風(fēng)設(shè)計規(guī)范》（JTG/T 3360-01—2018）附錄A，Ⅴ型峽谷地形所處的吉首地區(qū)不考慮風(fēng)向地面以上10 m高度處100 a重現(xiàn)期、10 min平均的年最大風(fēng)速為24.3 m/s。

以橋址附近的“虛擬標準氣象站”基本風(fēng)速作為峽谷邊界入口的基本風(fēng)速，借用現(xiàn)有研究的經(jīng)驗關(guān)系式計算建于Ⅴ型峽谷橋面的設(shè)計基準風(fēng)速，經(jīng)驗關(guān)系為[27]

Ud=U10

1.0+

πH5（B1+B2）

=kU10

（1）

式中?Ud為設(shè)計基準風(fēng)速，m/s；U10為虛擬氣象站的基本風(fēng)速，m/s；H為橋面高度，m；B1為峽谷上口處寬度，一般取橋長，m；B2為峽谷下底面寬度，m；

k為峽谷效應(yīng)修正系數(shù)。

根據(jù)以上修正公式，參考矮寨大橋設(shè)計（圖24）：

U10=24.3 m/s

，B1=1 000 m，B2=250 m，H=340 m，得k=1.171。

從而橋面高度處的設(shè)計基準風(fēng)速：

Ud=1.171×24.3 m/s=28.46 m/s

。

由該經(jīng)驗公式計算得到的橋梁抗風(fēng)設(shè)計基準風(fēng)速與文獻[17]中的擬合計算公式基本一致，與文獻[28]設(shè)計階段的橋面基準風(fēng)速34.9 m/s相差18.5%，說明在現(xiàn)場實測風(fēng)速數(shù)據(jù)難以獲得的山區(qū)地帶，可參考經(jīng)驗公式初步獲得跨越Ⅴ型峽谷橋梁的基準風(fēng)速。

根據(jù)以上初步判斷的地表類型并利用現(xiàn)有最新的抗風(fēng)規(guī)范計算該橋面高度的設(shè)計基準風(fēng)速為

Ud=kf

Z10

α0Us10

（2）

式中?kf為抗風(fēng)風(fēng)險系數(shù)；α0為橋址處的地表粗糙度系數(shù)；

Z為橋面高度，m;

Us10為橋梁設(shè)計基本風(fēng)速，m/s。

橋梁設(shè)計基本風(fēng)速Us10計算式為

Us10=kcU10

（3）

式中?kc為基本風(fēng)速地表類別轉(zhuǎn)換系數(shù)。

根據(jù)《抗風(fēng)規(guī)范》可知，kf=1.0，α0=0.3，kc=0.564，從而橋梁設(shè)計基本風(fēng)速：

Ud=1.0×

34010

0.3×0.564×24.3 m/s=39.48 m/s。

經(jīng)驗公式與現(xiàn)有規(guī)范相差27.9%，因此在進行峽谷地形橋址處的設(shè)計風(fēng)參數(shù)計算時，應(yīng)綜合考慮橋梁抗風(fēng)規(guī)范、相關(guān)經(jīng)驗公式，并可借助橋梁工程研究中常用數(shù)值模擬方法[29]作為參考來進一步提高其準確性。

4?結(jié)?論

1）受峽谷地形的影響，在3種湍流模型綜合對比并考慮地表粗糙度對近地面風(fēng)場干擾的條件下，該橋址地區(qū)范圍內(nèi)地形最優(yōu)模型為流體壁面粗糙度20 m的RNG k-ε模型，此模型可為后續(xù)類似地貌下橋梁結(jié)構(gòu)的風(fēng)場數(shù)值模擬研究提供參考。

2）受山體高度和山體夾角兩種地形參數(shù)的影響，該類峽谷地形在兩參數(shù)變化時的梯度風(fēng)高度分別為700 m和800 m左右，大于規(guī)范中4類地表的取值，且在離地高度40 m范圍內(nèi)，湍流強度更接近于D類地表；在梯度風(fēng)高度處，山體高度和夾角對風(fēng)速增強影響范圍分別為2.37%～12.56%和1.24%～6.98%，明顯山體高度對此峽谷地形的影響效果要大于山體夾角的變化，同時由于兩側(cè)山體的干擾，湍流強度在一定高度處（400～500 m）會有一段停滯期。

3）由于峽谷兩側(cè)山體的遮擋作用，風(fēng)剖面分布呈Z型，在300～500 m出現(xiàn)兩次拐點，之后風(fēng)速大致符合規(guī)范中冪指數(shù)的增長規(guī)律；根據(jù)峽谷地形參數(shù)分析，兩側(cè)山體夾角越小，風(fēng)速值越大，山體高度越高，風(fēng)速值增長越快。

4）依據(jù)Ⅴ型峽谷橋梁的經(jīng)驗公式得出的橋面高度處的設(shè)計基準風(fēng)速與按抗風(fēng)規(guī)范計算得出的風(fēng)速相比，更加接近實際橋址區(qū)。因此，針對此類地形風(fēng)洞試驗或?qū)嶋H工程項目的研究，在沒有現(xiàn)場實測數(shù)據(jù)支撐時，可以通過“數(shù)值風(fēng)洞”計算得出風(fēng)速數(shù)據(jù)及各類風(fēng)場特性作為參考。

參考文獻（References）：

[1]

UCHIDA T，OHYA Y.Large-eddy simulation of turbulent airflow over complex terrain[J].Journal of Wind Engineering and Industrial Aerodynamics，2003，91（1-2）：219-229．

[2]BULLARD J E，WIGGS G F S，NASH D J.Experimental study of wind directional variability in the vicinity of a model valley[J].Geomorphology，2000（35）：127-143.

[3]劉黎陽，張志田，汪志雄，等.不同規(guī)模地形模型對某山區(qū)橋梁設(shè)計風(fēng)特性確定的影響[J].實驗流體力學(xué)，2018，32（6）：49-54.

LIU Liyang，ZHANG Zhitian，WANG Zhixiong，et al.

Scope effects of terrain models on wind properties design of a bridge located at mountainous area[J].Journal of Experiments in Fluid Mechanics，2018，32（6）：49-54.

[4]周繼，王新國，劉志文，等.大跨度山區(qū)橋梁風(fēng)特性數(shù)值模擬及試驗研究[J].鐵道科學(xué)與工程學(xué)報，2022，19（4）：995-1003.

ZHOU Ji，WANG Xinguo，LIU Zhiwen，et al.Numerical simulation and experimental study on wind characteristics of long span bridge in mountainous area[J].Journal of Railway Science and Engineering，2022，19（4）：995-1003.

[5]張向旭，鄭史雄，唐煜，等.山區(qū)峽谷大跨度拱橋橋址區(qū)風(fēng)參數(shù)數(shù)值模擬研究[J].鐵道科學(xué)與工程學(xué)報，2018，15（2）：398-406.

ZHANG Xiangxu，ZHENG Shixiong，TANG Yu，et al.Numerical study on wind parameters for a long span arch bridge in the area of a gorge[J].Journal of Railway Science and Engineering，2018，15（2）：398-406.

[6]李永樂，遆子龍，汪斌，等.山區(qū)Y形河口附近橋址區(qū)地形風(fēng)特性數(shù)值模擬研究[J].西南交通大學(xué)學(xué)報，2016，51（2）：341-348.

LI Yongle，TI Zilong，WANG Bin，et al.Numerical simulation of wind characteristics over bridge site near Y-shaped river junction in mountainous area[J].Journal of Southwest Jiaotong University，2016，51（2）：341-348.

[7]于艦涵，李明水，廖海黎.山區(qū)地形對橋位風(fēng)場影響的數(shù)值模擬[J].西南交通大學(xué)學(xué)報，2016，51（4）：654-662.

YU Jianhan，LI Mingshui，LIAO Haili.Numerical simulation of effect of mountainous topography on wind field at bridge site[J].Journal of Southwest Jiaotong University，2016，51（4）：654-662.

[8]薛亞飛.復(fù)雜地形橋位風(fēng)場實測與CFD數(shù)值模擬[D].長沙：湖南大學(xué)，2015.

XUE Yafei.Field measurement and CFD numerical si-mulation of wind over bridge site with complex terrain[D].Changsha：Hunan University，2015.

[9]孫亭亭，楊吉新，黎建華，等.基于脈動風(fēng)作用下的大跨度橋梁拉索風(fēng)場數(shù)值模擬[J].山東農(nóng)業(yè)大學(xué)學(xué)報（自然科學(xué)版），2018，49（5）：809-814.

SUN Tingting，YANG Jixin，LI Jianhua，et al.Numerical simulation of wind field on large span bridge cable based on pulsating winds[J].Journal of Shandong Agricultural University（Natural Science Edition），2018，49（5）：809-814.

[10]王鐵強.實際地形下風(fēng)場CFD模擬與參數(shù)分析[D].廣州：中山大學(xué)，2009.

WANG Tieqiang.CFD wind field simulation and para-meter analysis over natural terrain[D].Guangzhou：Sun Yat-sen University，2009.

[11]于濤.山區(qū)峽谷風(fēng)特性的數(shù)值模擬與現(xiàn)場實測的對比研究[D].湘潭：湖南科技大學(xué)，2014.

YU Tao.Study on wind-environment numerieal simulation and field measurement about canyon[D].Xiangtan：Hunan University of Science and Technology，2014.

[12]張玥.西部山區(qū)谷口處橋位風(fēng)特性觀測與風(fēng)環(huán)境數(shù)值模擬研究[D].西安：長安大學(xué)，2009.

ZHANG Yue.Field investigation and wind-environment numerical simulation of bridge site in the western gap of valley areas[D].Xian：Changan University，2009.

[13]唐昂.基于實際峽谷地形的近地風(fēng)場特性研究[D].杭州：浙江理工大學(xué)，2015.

TANG Ang.Study on characteristcs of wind field in ground layer on actual canyon terrain[D].Hangzhou：Zhejiang Sci-Tech University，2015.

[14]唐金旺.喇叭口地形橋址強風(fēng)特性風(fēng)洞試驗及CFD模擬研究[D].西安：西安科技大學(xué)，2017.

TANG Jinwang.The strong wind characteristics at bridge site of Trumpet-shaped topography with wind tunnel test and CFD simulation[D].Xian：Xian University of Science and Technology，2017.

[15]姚劍鋒，沈國輝，姚旦，等.峽谷和埡口地形風(fēng)場特征的CFD數(shù)值模擬[J].哈爾濱工業(yè)大學(xué)學(xué)報，2016，48（12）：165-171.

YAO Jianfeng，SHEN Guohui，YAO Dan，et al.CFD-based numerical simulation of wind field characteristics on valley and col terrain[J].Journal of Harbin Institute of Technology，2016，48（12）：165-171.

[16]高亮.內(nèi)陸強風(fēng)特性的現(xiàn)場實測與模擬[D].西安：長安大學(xué)，2012.

GAO Liang.Simulation and field measurement of the characters of strong wind in inner land of China[D].Xian：Changan University，2012.

[17]洪新民，郭文華，熊安平.山區(qū)峽谷風(fēng)場分布特性及地形影響的數(shù)值模擬[J].長安大學(xué)學(xué)報（自然科學(xué)版），2017，37（5）：56-64.

HONG Xinmin，GUO Wenhua，XIONG Anping.Numerical simulation of distribution characteristics of wind fields and terrains influence in mountain canyon[J].Journal of Changan University（Natural Science Edition），2017，37（5）：56-64.

[18]馮林.基于峽谷地形的風(fēng)場數(shù)值模擬研究[D].杭州：浙江理工大學(xué)，2017.

FENG Lin.Numerical simulation of wind field based on canyon terrain[D].Hangzhou：Zhejiang Sci-Tech University，2017.

[19]樓文娟，劉萌萌，李正昊，等.峽谷地形平均風(fēng)速特性與加速效應(yīng)[J].湖南大學(xué)學(xué)報（自然科學(xué)版），2016，43（7）：8-15.

LOU Wenjuan，LIU Mengmeng，LI Zhenghao，et al.Research on mean wind speed characteristics and speed-up effect in canyon terrain[J].Journal of Hunan University（Natural Sciences），2016，43（7）：8-15.

[20]陳平.地形對山地丘陵風(fēng)場影響的數(shù)值研究[D].杭州：浙江大學(xué)，2007.

CHEN Ping.Numerical study of terrain influence on the airflow over hilly land[D].Hangzhou：Zhejiang University，2007.

[21]李正昊，樓文娟，章李剛，等.地貌因素對埡口內(nèi)風(fēng)速影響的數(shù)值模擬[J].浙江大學(xué)學(xué)報（工學(xué)版），2016，50（5）：848-855.

LI Zhenghao，LOU Wenjuan，ZHANG Ligang，et al.Numerical simulation of effects of topographic factors on wind speed in col[J].Journal of Zhejiang University（Engineering Science），2016，50（5）：848-855.

[22]BOWEN A J，LINDLEY D.A wind-tunnel investigation of the wind speed and turbulence characteristics close to the ground over various escarpment shapes[J].Boundary-Layer Meteorology，1977，12（3）：259-271.

[23]CAO S，TAMURA T.Experimental study on roughness effects on turbulent boundary layer flow over a two-dimensional steep hill[J].Journal of Wind Engineering and Industrial Aerodynamics，2006，94（1）：1-19.

[24]LUBITZ W D，WHITE B R.Wind-tunnel and field investigation of the effect of local wind direction on speed-up over hills[J].Journal of Wind Engineering and Industrial Aerodynamics，2007，95（8）：639-661.

[25]童祎巍.矮寨大橋橋址區(qū)脈動風(fēng)特性及高斯條件模擬研究[D].湘潭：湖南科技大學(xué)，2015.

TONG Yiwei.Fluctuating wind characteristics of Aizhai bridge area and of its gaussian condition simulations[D].Xiangtan：Hunan University of Science and Technology，2015.

[26]丁冬.典型山區(qū)地形橋位良態(tài)風(fēng)特性研究[D].長沙：湖南大學(xué)，2017.

DING Dong.Study on normal wind characteristics over bridge sites with typical mountainous terrain[D].Changsha：Hunan University，2017.

[27]

王凱，廖海黎，李明水，等.山區(qū)峽谷橋梁設(shè)計基準風(fēng)速的確定方法[J].西南交通大學(xué)學(xué)報，2013，48（1）：29-35.

WANG Kai，LIAO Haili，LI Mingshui，et al.Determination method for basic design wind speed of mountainous-valley bridge[J].Journal of Southwest Jiaotong University，2013，48（1）：29-35.

[28]陳政清，李春光，張志田，等.山區(qū)峽谷地帶大跨度橋梁風(fēng)場特性試驗[J].實驗流體力學(xué)，2008（3）：54-59，67.

CHEN Zhengqing，LI Chunguang，ZHANG Zhitian，et al.Model test study of wind field characteristics of long-span bridge site in mountainous valley terrain[J].Journal of Experiments in Fluid Mechanics，2008（3）：54-59，67.

[29]谷拴成，姚博語，任翔，等.煤礦輸煤棧橋抗振性能及動力響應(yīng)[J].西安科技大學(xué)學(xué)報，2022，42（1）：8-15.

GU Shuancheng，YAO Boyu，REN Xiang，et al.Research on resonance and dynamic response of coal conveying trestle in coal mine[J].Journal of Xian University of Science and Technology，2022，42（1）：8-15.

（責(zé)任編輯：劉潔）