橫通道設(shè)置對車隧氣動多波峰的影響

王英學(xué), 任文強(qiáng), 駱 陽, 何 俊

(西南交通大學(xué)交通隧道工程教育部重點(diǎn)實驗室, 四川 成都 610031)

列車突入隧道時,會在隧道內(nèi)產(chǎn)生壓縮波,并在隧道出口形成微氣壓波釋放,影響洞口周圍環(huán)境[1-2].

文獻(xiàn)[3-4]中分析了隧道緩沖結(jié)構(gòu)的設(shè)計參數(shù)及降低壓力梯度的理論,得出一些緩沖結(jié)構(gòu)設(shè)計的規(guī)律、方法.文獻(xiàn)[5]中論述了國外有關(guān)治理車隧氣動效應(yīng)的治理措施及微氣壓波的放大現(xiàn)象,并對我國高鐵隧道設(shè)計參數(shù)的選取進(jìn)行了探討.

由于微氣壓波存在放大效應(yīng),僅在洞口設(shè)置緩沖結(jié)構(gòu)有時難以滿足控制微氣壓波要求.實際中,在隧道主洞兩側(cè)開設(shè)支洞,以加快工程進(jìn)度，在不增加工程投資的情況下,可以起到控制微氣壓波的目的.

對于橫通道在緩解微氣壓波中的應(yīng)用,一些學(xué)者已經(jīng)進(jìn)行了研究.文獻(xiàn)[6]中對隧道支洞緩解微氣壓波效果進(jìn)行研究和測試,微氣壓波壓力峰值與隧道出口附近洞內(nèi)壓縮波梯度峰值成正比,表明隧道支洞可以起到緩解微壓波的效果.文獻(xiàn)[7-9]中分析了橫通道緩解微壓波效果，并對結(jié)構(gòu)參數(shù)進(jìn)行了優(yōu)化.

本文將采用理論分析和數(shù)值方法,對隧道設(shè)置橫通道后的壓力曲線進(jìn)行分析,探討分析曲線波動機(jī)理,通過壓力梯度討論橫通道對微氣壓波的影響.以期更好地指導(dǎo)工程設(shè)計.

1 隧道支洞緩解微氣壓波機(jī)理

在實際工程中,為加快施工進(jìn)度所修建的支洞形式多種多樣,可以采用單側(cè)或雙側(cè)布置,橫通道可以是單個或多個.對于橫通道設(shè)施各個支洞的名稱,一般將平行于隧道主洞的通道稱為平導(dǎo)(即平行導(dǎo)洞),與隧道和平導(dǎo)垂直的支洞稱為橫通道.

橫通道及平導(dǎo)使隧道洞口區(qū)段所產(chǎn)生的氣動效應(yīng)復(fù)雜化[10].橫通道的功能與豎井、斜井等的功能類似,只是豎井、斜井距離隧道入口的距離相對較遠(yuǎn)[11].文獻(xiàn)[12]中對壓縮波經(jīng)過豎井的氣動機(jī)理進(jìn)行分析,明確了壓縮波在豎井與隧道交叉點(diǎn)的傳播、反射特性.

基于對豎井的入射波反射規(guī)律的分析,可以確定橫通道壓力波的傳播規(guī)律,將橫通道的氣動壓力變化分為3個階段:

(1) 首波傳播及泄壓(首波標(biāo)記為①,壓縮波為實線,膨脹波為虛線,下同)

該過程是首波到達(dá)橫通道處的壓力釋放過程,與豎井的傳播過程一致.

列車進(jìn)入隧道產(chǎn)生壓縮波,傳播到豎井與隧道連接處時,一部分波以反號波反射,另一部分波則以2個同號波向前傳播,如圖1(a)所示，圖中：lh為橫通道距離隧道洞口；dh為主洞與平導(dǎo)間距.當(dāng)沿橫通道向外傳播的波到達(dá)橫通道出口時,以反號波反射回來,向隧道入口傳播的波到達(dá)隧道入口也產(chǎn)生同樣反射,如此往復(fù)傳播.

(2) 2次波產(chǎn)生與傳播

列車到達(dá)橫通道處,由于空間突變,列車經(jīng)過再次形成壓縮波.該波分為2個分支,隧道主洞部分(標(biāo)記為②)和橫通道部分(標(biāo)記為③),均為壓縮波,如圖1(b)所示.

(3) 2次波反射

分配到橫通道中的壓縮波,到達(dá)橫通道外開闊空氣域轉(zhuǎn)變?yōu)榕蛎洸?向隧道方向傳播,在橫通道與隧道連接處分為3支:隧道入口方向(膨脹波)、隧道出口方向(膨脹波)和反射回橫通道方向(轉(zhuǎn)換為壓縮波).

分配到隧道入口方向傳播的膨脹波到達(dá)隧道入口處,釋放部分能量后,轉(zhuǎn)為壓縮波向隧道出口方向傳播.當(dāng)傳播經(jīng)過橫通道時,其傳播、反射與首波的特征一致.由于該波為壓縮波,因此,可引起隧道內(nèi)壓力升高,使壓力梯度波形曲線出現(xiàn)正壓波峰.

(a) 首波

(b) 2次波產(chǎn)生、傳播與反射圖1 橫通道氣動壓力變化Fig.1 Compression wave transmission and reflection at cross adit

2 數(shù)值分析

2.1 紊流模型的選取及列車進(jìn)入隧道過程的實現(xiàn)

在流體分析時,紊流模型的選取對計算結(jié)果有很大影響,可選的紊流模式也較多.直接數(shù)值模擬(direct numerical simulation, DNS)不用任何湍流模型，直接求解完整的三維非定常的N-S方程組.雷諾平均模型(Reynolds-averaged Navier-Stokes equations， RANS)應(yīng)用紊流統(tǒng)計理論,將非穩(wěn)態(tài)的N-S方程對時間做平均,其工程應(yīng)用的紊流模型有雷諾應(yīng)力模型(Reynolds stress model，RSM),該模型把所有尺度的雷諾應(yīng)力都使用模型進(jìn)行模擬.大渦模擬(large eddy simulation, LES)中，動量、能量、質(zhì)量及其標(biāo)量主要由大渦輸送,對大渦及小渦采取不同處理策略,進(jìn)行流場的分析計算.在k-ε模型中,主要求解紊動能和耗散率方程模擬紊流狀態(tài).

許多學(xué)者在車隧氣動分析中,所選取的紊流模型也略有差別,比較常用的有k-ε模型[13]、Les模型[14]等.大渦模型紊流模型被認(rèn)為是具有應(yīng)用前景的紊流模式.

本文利用商用工程軟件Ansys中的Fluent模塊,開展列車進(jìn)入隧道氣動過程的模擬分析.紊流模型選用大渦模型(LES),計算中介質(zhì)采用理想氣體.考慮到大渦模型對網(wǎng)格的質(zhì)量要求較高,在建模及網(wǎng)格劃分時,對車頭、隧道入口等一些關(guān)鍵部位進(jìn)行了網(wǎng)格加密處理.

采用Fluent模塊中的動網(wǎng)格法[15]可以模擬列車-隧道之間的相對運(yùn)動,文獻(xiàn)[7,11-12,16-17]中采用該方法進(jìn)行模擬分析,取得了良好效果.本文也將采用該方法,模擬列車通過隧道過程.研究結(jié)果表明，本文的計算結(jié)果與現(xiàn)場測試數(shù)據(jù)[18]的壓縮波壓力峰值接近,表明選取的模型是正確的.

2.2 計算參數(shù)

在計算中,列車采用CHR3型動車,機(jī)車橫截面積為11.624 m2.車長為100 m,車速vc=350 km/h.隧道面積為100 m2,時速350 km/h雙線隧道標(biāo)準(zhǔn)設(shè)計斷面.隧道洞口設(shè)置1個橫通道,lh=20 m,dh=15 m,橫通道面積為29 m2,平導(dǎo)段隧道面積為35 m2.模型的隧道長度為500 m,隧道兩端空氣域長度為300 m,列車距離隧道入口為100 m,隧道及橫通道布置如圖2所示.

2.3 數(shù)值模型及網(wǎng)格劃分

在計算中,模型的最小網(wǎng)格尺寸為0.02 m×0.02 m,并在某些關(guān)鍵區(qū)域進(jìn)行了網(wǎng)格細(xì)化.隧道出入口區(qū)域采用無反射邊界,車頭、隧道、橫通道模型網(wǎng)格如圖3所示.計算中的迭代次數(shù)為7次,計算步長為0.000 5 s,最大殘差0.001.

圖2 數(shù)值計算中隧道及橫通道布置示意(單位:m)Fig.2 Tunnel and cross adit setting dimension in numeral simulation (unit: m)

(a) 列車車頭

(b) 隧道入口及橫通道圖3 數(shù)值計算模型網(wǎng)格劃分Fig.3 Numeral simulation mesh

2.4 計算結(jié)果

通過計算,得到在無橫通道時,距離隧道洞口10、15、20 m時隧道內(nèi)測點(diǎn)壓力及壓力梯度監(jiān)測曲線，如圖4所示.橫通道距離隧道入口20 m,設(shè)置橫通道后,隧道內(nèi)距隧道入口200 m處測點(diǎn)壓力及壓力梯度曲線如圖5所示,各曲線的壓縮波及壓力梯度峰值如表1所示.

(a) 壓力(b) 壓力梯度圖4 無橫通道時隧道內(nèi)測點(diǎn)壓力及壓力梯度變化曲線Fig.4 Variation of pressure and pressure gradient in tunnel without cross adit

(a) 壓力(b) 壓力梯度圖5 隧道內(nèi)測點(diǎn)壓力及壓力梯度曲線Fig.5 Variation of pressure and pressure gradient in tunnel

表1 隧道內(nèi)測點(diǎn)壓力及壓力梯度峰值Tab.1 Pressure and pressure gradient peak value in tunnel

2.5 計算結(jié)果分析

2.5.1波動曲線

從有、無橫通道的壓力及壓力梯度曲線可知:

(1) 壓力梯度在列車進(jìn)入隧道較短距離即達(dá)到峰值,壓縮波在列車完全進(jìn)入隧道前一直增加.

列車進(jìn)入隧道10～20 m過程中,壓縮波處于上升期,曲線峰值是逐漸增加的(如圖4(a)),而壓力梯度很快上升,在10 m左右已經(jīng)達(dá)到峰值(如圖4(b)).

(2) 隧道設(shè)置橫通道與無橫通道相比,壓力曲線波動相差較大,會出現(xiàn)次級波峰、波谷.

在無橫通道時,隧道內(nèi)200 m處壓力梯度可分為3段:迅速上升區(qū)段、平緩上升區(qū)段、迅速下降區(qū)段.迅速上升區(qū)段為列車頭部段突入隧道所產(chǎn)生;平緩上升區(qū)段為列車車身進(jìn)入隧道階段所產(chǎn)生;迅速下降區(qū)段為列車車尾進(jìn)入隧道所產(chǎn)生的膨脹波所引起.設(shè)置橫通道以后,隧道內(nèi)200 m處壓縮波曲線的迅速上升區(qū)段和迅速下降區(qū)段與無橫通道的情況基本一致;在平緩上升區(qū)段則出現(xiàn)多個波峰、波谷.

2.5.2壓力曲線多波出現(xiàn)機(jī)理

(1) 首波、2次波峰值壓力梯度峰值形成距離差異無論是壓縮波還是壓力梯度,要到達(dá)峰值都需要一定時間,為更好分析壓力梯度多峰值出現(xiàn)機(jī)理,假定列車進(jìn)入隧道形成首波壓力梯度峰值所需距離為x1,列車經(jīng)過橫通道形成2次波峰值所需距離為x2.根據(jù)前面對洞口多測點(diǎn)壓力梯度對比,取x1=10 m.圖6為壓力波在隧道內(nèi)傳播到達(dá)時刻計算簡圖.

圖6 隧道內(nèi)測點(diǎn)壓力波動特性分析簡圖Fig.6 Calculation graph of pressure variation at measure point in tunnel

計算中,橫通道距離隧道入口lh=20 m,橫通道在隧道方向?qū)挾葁h=5.9 m,壓縮波傳播速度為vb=340 m/s.以隧道入口為時間起始點(diǎn),則首波和2次波分別到達(dá)隧道內(nèi)P處距離所用時間為

t1=x1/vc+(lh-x1+wh+x2+x)/vb，

(1)

t2=(lh+wh+x2)/vc+x/vb，

(2)

Δt1=t2-t1=(lh+wh+x2)/vb+x/vb-

x1/vc-(lh-x1+wh+x2+x)/vb=

(lh+wh+x2-x1)/vc-

(lh+wh+x2-x1)/vb.

記Δx=x1-x2,則有

Δt1=(lh+wh-Δx)/vc-(lh+wh-Δx)/vb，

Δx=((lh+wh)(1/vc-1/vb)-Δt1)/

(1/vc-1/vb).

經(jīng)計算得到,隧道內(nèi)200 m測點(diǎn)處首波、2次波壓力梯度峰值出現(xiàn)時刻分別為t1=1.648 s，t2=1.826 s，則，Δt1=0.178 s.計算得到Δx=1.7 m,x1>x2.

(2) 2次波波谷

2次波波谷產(chǎn)生的原因可能有兩種:

① 列車進(jìn)入隧道所形成的首波分配到橫通道的壓縮波分支,傳播至橫通道出口后,反射成為膨脹波向隧道方向傳播.該膨脹波到達(dá)橫通道與隧道交叉處,又分為3支:第1支轉(zhuǎn)換為壓縮波又向橫通道出口傳播;第2支作為膨脹波向隧道入口方向傳播;第3支作為膨脹波向隧道出口方向傳播,該分支正是產(chǎn)生2次波波谷的原因.

在該情況下,引起2次波波谷的波源與隧道內(nèi)首波是一致的,只是行走路徑不同,在橫通道內(nèi)有1個折返的過程,于是得到橫通道反射膨脹波與首波歷程時間差.

以橫通道的中心線作為波的傳播路徑,隧道入口為時間起算點(diǎn),則有

t1=x1/vc+(lh-x1+wh+x2+x)/vb，

(3)

t2=x1/vc+(lh-x1)/vb+2(wh+dh+lh)/vb+

(x2+x)/vb.

(4)

于是有

Δt2=2(dh+lh)/vb+wh/vb=0.223 s.

② 列車到達(dá)橫通道時,產(chǎn)生2次波在隧道內(nèi)傳播的同時,在橫通道內(nèi)也產(chǎn)生壓縮波.該橫通道壓縮波傳播至橫通道出口后,反射成為膨脹波向隧道方向傳播.該膨脹波到達(dá)橫通道與隧道交叉處,又分為3支:第1支轉(zhuǎn)換為壓縮波又向橫通道出口傳播;第2支作為膨脹波向隧道入口方向傳播;第3支作為膨脹波向隧道出口方向傳播,該分支正是產(chǎn)生2次波波谷的原因.

在該情況下,引起2次波波谷的波源是由列車到達(dá)橫通道時誘發(fā)引起的壓縮波,在橫通道內(nèi)傳播并折返的過程.于是得到橫通道反射膨脹波與首波歷程時間差.

以橫通道的中心線作為波的傳播路徑，隧道入口為時間起算點(diǎn),則有

t1=x1/vc+(lh-x1+wh+x2+x)/vb，

(5)

t2=(lh+wh)/vc+2(dh+lh)/vb+

(x2+x)/vb.

(6)

于是有

Δt2=(lh+wh-x1)/vc+

(2dh+lh+x1-wh)/vb=0.323 s.

由計算曲線上檢測到的Δt2=0.314 s.所以,第2種誘發(fā)原因更為合理.

(3) 橫通道產(chǎn)生壓縮波的2次反射波(壓縮波).

對2次波波谷的分析表明，從橫通道傳播到隧道都是膨脹波.在監(jiān)測曲線中看到有第3個波峰,而波峰一定是壓縮波所引起的,其產(chǎn)生原因為分配到橫通道中的壓縮波,經(jīng)橫通道出口及隧道入口反射,再次向隧道出口方向傳播的壓縮波.

參考圖6,首波及該壓縮波到達(dá)隧道內(nèi)測點(diǎn)P的時刻為

t1=x1/vc+(lh-x1+wh+x2+x)/vb，

(7)

t2=(lh+wh)/vc+2(dh+lh)/vb+

(2lh+wh)/vb+(x2+x)/vb.

(8)

于是有

Δt3=(lh+wh-x1)/vc+

(2dh+3lh+x1)/vb=0.457 s.

實際檢測曲線的時間差為Δt3=0.45 s.

上述計算,隧道內(nèi)測點(diǎn)壓力梯度峰值出現(xiàn)時間差如表2所示.

表2 壓力梯度多波峰出現(xiàn)時刻對照表Tab.2 Pressure gradient multiple peak appearing period comparison table s

2.6 壓力波傳播特性應(yīng)用

通過分析可知,列車在經(jīng)過的設(shè)置橫通道后,隧道內(nèi)壓力出現(xiàn)多波峰、波谷的原因是在橫通道處壓縮波、膨脹波的傳播、反射.利用這些特性,可以更好理解壓力梯度曲線中波峰、波谷的產(chǎn)生原因.

針對第3壓力梯度峰值,可知橫通道的不利位置.分析思路如下.

當(dāng)壓縮波首波分配到橫通道內(nèi)的分支,傳播到達(dá)橫通道出口,反射到達(dá)橫通道與隧道交叉點(diǎn)處時刻與列車經(jīng)過該點(diǎn)時刻重合,會使分配到橫通道的壓縮波疊加，則有

(9)

算例中,x1=10 m,dh=15 m,則lh=49 m.

lh=50 m和lh=20 m時的壓力梯度對比分析曲線如圖7所示.

(a) 壓力

(b) 壓力梯度圖7 不同橫通道設(shè)置位置時隧道內(nèi)測點(diǎn)壓力及壓力梯度曲線(測點(diǎn)距隧道入口200 m)Fig.7 Pressure and pressure gradient curve at measure point in tunnel with different cross adit setting position( the measure point is 200 m to tunnel entrance )

P2=6.6 kPa/s,P3=5.3 kPa/s.

3 結(jié) 論

在設(shè)置的橫通道條件下,通過對車隧氣動特性的理論分析和數(shù)值模擬,得到如下結(jié)論:

(1) 在列車進(jìn)入隧道過程中,壓縮波壓力達(dá)到峰值需要較長時間,在列車車體進(jìn)入隧道過程中一直是增加的,而壓力梯度達(dá)到峰值所需時間較短,在列車進(jìn)入隧道10 m左右,即達(dá)到峰值.

(2) 列車進(jìn)入隧道形成首波壓力梯度峰值所需運(yùn)行距離要略大于列車經(jīng)過橫通道形成2次壓力梯度波峰值所需運(yùn)行距離;

(3) 列車經(jīng)過橫通道時,在橫通道內(nèi)也產(chǎn)生壓縮波,其反射、傳播影響隧道內(nèi)壓力梯度多個波峰、波谷的出現(xiàn)和峰值的大小.

(4) 隨著橫通道距隧道入口距離增大,列車經(jīng)過橫通道所產(chǎn)生的次級壓力波動效應(yīng)會有所降低.橫通道位置設(shè)置在50 m與設(shè)置在20 m處相比,首波壓力梯度峰值基本不變,但次級波的壓力梯度峰值更小.

[1] 趙文誠,高波,漆泰岳. 高速鐵路隧道出口微壓波及其主被動減緩措施[J]. 石家莊鐵道學(xué)院學(xué)報,2004,17(3): 5-9.

ZHAO Wencheng, GAO Bo, QI Taiyue. Micro-pressure wave around high-speed railway tunnel exit and active and passive method of reduction[J]. Journal of Shijiazhuang Railway Institute, 2004, 17(3): 5-9.

[2] 趙文成,高波,王英學(xué). 高速鐵路隧道口微壓波減緩措施分析[J]. 地下空間與工程學(xué)報,2005(2): 274-277.

ZHAO Wencheng, GAO Bo, WANG Yingxue. Analysis of the mitigation methods of micro-pressure wave generated around high-speed railway tunnel exits[J]. Chinese Journal of Underground Space and Engineering, 2005(2): 274-277.

[3] HOWE M S, LIDA M, MAEDA T T， et al. Aero-acoustics of a tunnel-entrance hood with a rectangular window[J]. Journal of Fluid Mechanics, 2003, 4(8): 211-243.

[4] HOWE M S. Design of a tunnel-entrance hood with multiple windows and variable cross-section[J]. Journal of Fluids and Structures, 2003, 17(8): 1111-1121.

[5] 王建宇. 列車通過隧道時誘發(fā)的空氣動力學(xué)問題和高速鐵路隧道設(shè)計參數(shù)[J]. 世界隧道,1995(1): 3-13.

[6] YAMAMOTO A, OZAWA S, MAEDA T. Reduction of micro-pressure wave radiated from tunnel exit by side branches in tunnel[J]. Quarterly Report of Rail way Technical Research Institute of Japanese National Rail Ways, 1984, 25(3): 102-105.

[7] 李志偉,梁習(xí)鋒,張鍵. 橫通道對緩解隧道瞬變壓力的研究[J]. 鐵道科學(xué)與工程學(xué)報,2010,7(4): 37-41.

LI Zhiwei, LIANG Xifeng, ZHANG Jian. Study of alleviating transient pressure with cross passage in a tunnel[J]. Journal of Railway Science and Engineering, 2010, 7(4): 37-41.

[8] 蒲佳. 豎井、橫通道對高速鐵路隧道出口微壓波影響分析[D]. 成都:西南交通大學(xué),2009.

[9] 戴文,高波,趙文成. 高速鐵路隧道橫通道氣動效應(yīng)研究[J]. 鐵道建筑,2011(12): 57-60.

DAI Wen, GAO Bo, ZHAO Wencheng. Study on aerodynamic effects in cross passage of high speed railway tunnel[J]. Railway Engineering, 2011(12): 57-60.

[10] 岳楹沁,高波. 橫通道對隧道出口微壓波影響的數(shù)值模擬研究[J]. 華東交通大學(xué)學(xué)報,2006,23(5): 20-23.

YUE Yingqin, GAO Bo. Study of the Influence over cross passage on high-speed railway tunnel exit micro-pressure wave[J]. Journal of East China Jiaotong University, 2006, 23(5): 20-23.

[11] 劉佩斯. 豎井對高速鐵路隧道氣動效應(yīng)的影響研究[D]. 成都:西南交通大學(xué),2009.

[12] 駱建軍,高波,王英學(xué),等. 高速列車穿越有豎井隧道流場的數(shù)值模擬[J]. 西南交通大學(xué)學(xué)報,2004,39(4): 442-446.

LUO Jianjun, GAO Bo, WANG Yingxue, et al. Simulation of unsteady three-dimensional flow induced by high-speed train entering tunnel with shaft[J]. Journal of Southwest Jiaotong University, 2004, 39(4): 442-446.

[13] 駱建軍,吳盡,陳鵬飛. 高速鐵路并聯(lián)隧道橫通道對隧道內(nèi)壓力變化的影響[J]. 北京交通大學(xué)學(xué)報,2015,39(1): 8-13.

LUO Jianjun, WU Jin, CHEN Pengfei. Pressure change from the cross aisle when the train passing through parallel tunnel of high-speed railway[J]. Journal of Beijing Jiaotong University, 2015, 39(1): 8-13.

[14] 馬靜,張杰,楊志剛. 橫風(fēng)下高速列車非定?？諝鈩恿μ匦匝芯縖J]. 鐵道學(xué)報,2008,30(6): 109-114.

MA Jing, ZHANG Jie, YANG Zhigang. Study on the unsteady aerodynamic characteristics of a high-speed train under cross wind[J]. Journal of the China Railway Society, 2008, 30(6): 109-114.

[15] 張來平,鄧小剛,張涵信. 動網(wǎng)格生成技術(shù)及非定常計算方法進(jìn)展綜述[J]. 力學(xué)進(jìn)展,2010,40(4): 424-447.

ZHANG Laiping, DENG Xiaogang, ZHANG Hanxin. Reviews of moving grid generation techniques and numerical methods for unsteady flow[J]. Advances in Mechanics, 2010, 40(4): 424-447.

[16] 趙文成,高波,駱建軍. 高速列車通過隧道的三維數(shù)值模擬[J]. 中國鐵道科學(xué),2003,24(5): 96-100.

ZHAO Wencheng, GAO Bo, LUO Jianjun. Three-dimensional numerical simulation of high-speed trains passing through tunnel[J].China Railway Science, 2003, 24(5): 96-100.

[17] 王英學(xué),高波,陸振華,等. 隧道中部存在車站時列車會車氣動效應(yīng)數(shù)值分析[J]. 空氣動力學(xué)學(xué)報,2009,27(3): 369-372.

WANG Yingxue, GAO Bo, LU Zhenhua, et al. Aerodynamics effect numeral simulation of high speed train passing through tunnel with mid-station[J]. Acta Aerodynamica Sinica, 2009, 27(3): 369-372.

[18] 趙有朋,馬偉斌,程愛君,等. 高速鐵路隧道氣動效應(yīng)[M]. 北京:中國鐵道出版社,2012： 139-144.