• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    2車道公路隧道射流風(fēng)機(jī)空間布局優(yōu)化的CFD分析

    2016-05-28 07:49:30閆治國
    隧道建設(shè)(中英文) 2016年4期

    趙 黎, 閆治國,2,3

    (1. 同濟(jì)大學(xué)土木工程學(xué)院地下建筑與工程系, 上海 200092; 2. 同濟(jì)大學(xué)土木工程防災(zāi)國家重點(diǎn)實(shí)驗(yàn)室, 上?!?00092; 3. 同濟(jì)大學(xué)巖土及地下工程教育部重點(diǎn)實(shí)驗(yàn)室, 上?!?00092)

    ?

    2車道公路隧道射流風(fēng)機(jī)空間布局優(yōu)化的CFD分析

    趙黎1, 閆治國1,2,3

    (1. 同濟(jì)大學(xué)土木工程學(xué)院地下建筑與工程系, 上海200092; 2. 同濟(jì)大學(xué)土木工程防災(zāi)國家重點(diǎn)實(shí)驗(yàn)室, 上海200092; 3. 同濟(jì)大學(xué)巖土及地下工程教育部重點(diǎn)實(shí)驗(yàn)室, 上海200092)

    摘要:為解決2車道公路隧道射流風(fēng)機(jī)的空間布局優(yōu)化問題,依托明堂山隧道工程,采用以往公路隧道射流風(fēng)機(jī)通風(fēng)效果數(shù)值模擬研究中不同的邊界條件,將隧道出入口邊界條件均設(shè)為大氣壓強(qiáng),并不預(yù)先給定隧道入口風(fēng)速大小。對影響射流風(fēng)機(jī)升壓力大小的因素,如風(fēng)機(jī)縱向間距、風(fēng)機(jī)布置高度、風(fēng)機(jī)橫向間距展開討論,得到: 風(fēng)機(jī)應(yīng)設(shè)在距建筑限界15~30 cm高度處;風(fēng)機(jī)橫向凈距應(yīng)設(shè)為1.5~2倍風(fēng)機(jī)直徑;風(fēng)機(jī)縱向間距應(yīng)設(shè)在150 m以上。明堂山隧道實(shí)際風(fēng)機(jī)布置方式所采用的參數(shù)均在優(yōu)化結(jié)果范圍內(nèi),按隧道實(shí)際長度及設(shè)計(jì)射流風(fēng)機(jī)臺(tái)數(shù)建模,模擬結(jié)果表明風(fēng)機(jī)在進(jìn)行優(yōu)化布局后,隧道通風(fēng)效果能夠達(dá)到設(shè)計(jì)要求。

    關(guān)鍵詞:明堂山公路隧道; 射流風(fēng)機(jī); 空間布局優(yōu)化; 數(shù)值模擬; CFD

    0引言

    伴隨著我國公路隧道的大發(fā)展,隧道縱向通風(fēng)方式得到越來越廣泛的應(yīng)用,而任何形式的縱向通風(fēng)方法都離不開射流風(fēng)機(jī)的作用。射流風(fēng)機(jī)在公路隧道中的布置位置受限于建筑限界與拱頂之間的空間,隧道中的風(fēng)機(jī)射流是一種有限空間射流,具有多股平行并列、多組縱向串聯(lián)、以通風(fēng)氣流為伴隨流動(dòng)等特點(diǎn)[1],因此,射流風(fēng)機(jī)的通風(fēng)效果與其空間布局有很大關(guān)系。

    我國很多學(xué)者都對優(yōu)化射流風(fēng)機(jī)在隧道空間中的布局問題展開過研究。如張林[1]對公路隧道射流風(fēng)機(jī)的安裝參數(shù)優(yōu)化進(jìn)行了數(shù)值模擬研究,分析了風(fēng)機(jī)與拱頂距離、風(fēng)機(jī)左右間距及不同功率風(fēng)機(jī)組合等對流態(tài)的影響;楊秀軍等[2]對公路隧道通風(fēng)中射流風(fēng)機(jī)縱向最小間距進(jìn)行了研究;胡純良等[3]對不同風(fēng)機(jī)軸距下某隧道內(nèi)氣流分布進(jìn)行了研究,得出在射流風(fēng)機(jī)后150 m處風(fēng)速基本不再變化,隧道內(nèi)氣流均勻性較好,處于完全穩(wěn)定的狀態(tài);方勇等[4]對3車道公路隧道射流風(fēng)機(jī)的設(shè)置位置展開了數(shù)值模擬分析,得出為了減小2股射流之間的干擾及隧道壁面影響,2臺(tái)風(fēng)機(jī)的最佳橫向凈距應(yīng)在3~4倍風(fēng)機(jī)直徑。

    在已有的對隧道射流通風(fēng)的數(shù)值模擬研究中,均給定了隧道入口風(fēng)速,但在實(shí)際施工中,隧道入口風(fēng)速是由射流風(fēng)機(jī)、交通風(fēng)和自然風(fēng)綜合影響決定的,這種邊界條件顯然與實(shí)際不符。因此,本文提出了將隧道入口設(shè)為壓力邊界條件(即給定大氣壓強(qiáng)值)的模擬方法,并不強(qiáng)制隧道入口風(fēng)速為定值,而是完全靠射流風(fēng)機(jī)的升壓作用吸進(jìn)風(fēng)量,即隧道入口的風(fēng)速是計(jì)算得到的。同時(shí),為簡化計(jì)算,忽略交通風(fēng)和自然風(fēng)的影響,將其提供的動(dòng)力或阻力用射流風(fēng)機(jī)或沿程阻力損失代替[5]。

    本文依托明堂山公路隧道,利用CFD軟件Fluent按隧道實(shí)際長度及設(shè)計(jì)射流風(fēng)機(jī)臺(tái)數(shù)建模,分析2車

    道公路隧道中射流風(fēng)機(jī)的安裝高度、橫向間距及縱向間距對通風(fēng)效果的影響,并將模擬結(jié)果應(yīng)用于優(yōu)化明堂山隧道中射流風(fēng)機(jī)的空間布局,應(yīng)用效果證明隧道通風(fēng)效果能夠達(dá)到設(shè)計(jì)要求。

    1基于壓力邊界的公路隧道射流通風(fēng)數(shù)值模擬方法

    1.1隧道模型說明

    岳武高速明堂山特長公路隧道長7 500 m,雙向4車道,單洞隧道截面寬11 m,高7 m,面積為65.18 m2。其中,左線采用全射流通風(fēng)方式,擬布置54臺(tái)1120型射流風(fēng)機(jī),見圖1。

    圖1 明堂山隧道左線風(fēng)機(jī)布置平面圖

    按照明堂山隧道的結(jié)構(gòu)參數(shù)建模,模型取500 m長,在隧道內(nèi)布置2組(每組2臺(tái))射流風(fēng)機(jī)。風(fēng)機(jī)出口風(fēng)速為30 m/s,入口風(fēng)量為29.56 m3/s。射流風(fēng)機(jī)用直徑1.12 m、長1.5 m的圓筒模擬[6],如圖2所示。

    圖2 隧道模型網(wǎng)格正視圖

    1.2邊界條件設(shè)定

    數(shù)值計(jì)算是在有限區(qū)域內(nèi)進(jìn)行的,因此,在區(qū)域的邊界上需要給定邊界條件。邊界條件要求在數(shù)學(xué)上滿足適定性,在物理上具有明顯的意義。本文的計(jì)算邊界類型及條件如表1所示。

    隧道進(jìn)出口均為壓力邊界條件,且設(shè)為大氣壓強(qiáng)值,這與隧道的實(shí)際情況相符。風(fēng)機(jī)進(jìn)口為質(zhì)量流量邊界條件,其數(shù)值為所選用風(fēng)機(jī)的吸入風(fēng)量參數(shù);風(fēng)機(jī)出口為速度邊界條件,其數(shù)值為所選用風(fēng)機(jī)的出口流速參數(shù),保證了質(zhì)量守恒定律[7]。隧道底部和頂部為壁面邊界條件,其粗糙高度由隧道實(shí)際粗糙度給定[8]。

    表1 計(jì)算模型邊界條件設(shè)置

    1.3通風(fēng)效果優(yōu)劣判斷標(biāo)準(zhǔn)

    (1)

    式中: vr為隧道斷面平均風(fēng)速; vj為射流風(fēng)機(jī)出口風(fēng)速; Ar為隧道截面積; Aj為風(fēng)機(jī)截面積; n為風(fēng)機(jī)臺(tái)數(shù)。

    (2)

    以往學(xué)者在對射流風(fēng)機(jī)空間布局進(jìn)行數(shù)值模擬分析時(shí),因給定了隧道入口風(fēng)速,并將出口相對靜壓強(qiáng)制為0,所以,往往將得到的入口靜壓與沿程損失之和作為射流提供升壓力的模擬值[9],但這樣得到的升壓力模擬值不能很好地體現(xiàn)射流風(fēng)機(jī)的升壓作用。本文在確定升壓力模擬值時(shí),是從計(jì)算結(jié)果中找出氣流在到達(dá)射流風(fēng)機(jī)前隧道截面靜壓的最小值和氣流經(jīng)射流風(fēng)機(jī)升壓后隧道截面靜壓的最大值,二者之差即為射流風(fēng)機(jī)提供的升壓力。

    2射流風(fēng)機(jī)空間布置位置優(yōu)化分析

    2.1射流風(fēng)機(jī)安裝高度對通風(fēng)效果的影響

    JTG/T D70/2-02—2014《公路隧道通風(fēng)設(shè)計(jì)細(xì)則》[10]中規(guī)定,射流風(fēng)機(jī)在隧道橫斷面上的布置應(yīng)不侵入建筑限界,射流風(fēng)機(jī)的邊沿與隧道建筑限界的凈距不宜小于15 cm。因此,為研究射流風(fēng)機(jī)安裝高度對通風(fēng)效果的影響,設(shè)定風(fēng)機(jī)距建筑限界分別為15、30、45、60 cm 4組工況進(jìn)行計(jì)算。

    第1組風(fēng)機(jī)距隧道入口50 m,第2組風(fēng)機(jī)距隧道入口200 m,氣流在到達(dá)各組射流風(fēng)機(jī)前隧道截面靜壓的最小值和氣流經(jīng)射流風(fēng)機(jī)升壓后隧道截面靜壓的最大值如圖3所示。

    圖3射流風(fēng)機(jī)前后截面最大及最小靜壓值(風(fēng)機(jī)距建筑限界不同時(shí))

    Fig. 3Maximum and minimum static pressures of tunnel cross-section(with different distances between fan and building boundary)

    風(fēng)機(jī)布置高度對升壓力的影響見表2。由表2可知,風(fēng)機(jī)布置得越高,越接近隧道拱頂,風(fēng)機(jī)調(diào)壓綜合影響系數(shù)k值越小,表明風(fēng)機(jī)提供的升壓力越小。因此,建議明堂山隧道在進(jìn)行風(fēng)機(jī)布置時(shí),風(fēng)機(jī)設(shè)在距建筑限界15~30 cm高度處。

    表2 風(fēng)機(jī)布置高度對升壓力的影響

    風(fēng)機(jī)距建筑限界60 cm和15 cm時(shí),風(fēng)機(jī)出口處風(fēng)速矢量圖如圖4和圖5所示。隨風(fēng)機(jī)布置高度增加沿風(fēng)機(jī)軸線剖面速度變化云圖如圖6所示。

    圖4風(fēng)機(jī)距建筑限界60 cm時(shí)風(fēng)機(jī)出口處風(fēng)速矢量圖(單位: m/s)

    Fig. 4Vectogram of air velocity at fan exit when distance between fan and building boundary of 60 cm (m/s)

    圖5風(fēng)機(jī)距建筑限界15 cm時(shí)風(fēng)機(jī)出口處風(fēng)度矢量圖(單位: m/s)

    Fig. 5Vectogram of air velocity at fan exit when distance between fan and building boundary of 15 cm (m/s)

    對比圖4和圖5可知,風(fēng)機(jī)布置在距建筑限界60 cm高度時(shí),其出口附近風(fēng)速明顯要比布置在距建筑限界15 cm高度時(shí)小。由圖6可知,隨風(fēng)機(jī)布置高度的增加,出口附近風(fēng)速越來越小。這是因?yàn)轱L(fēng)機(jī)出口風(fēng)流受到拱頂影響沒有足夠豎向空間充分?jǐn)U展,高速風(fēng)流撞擊到拱頂或受到隧道壁面摩阻影響導(dǎo)致動(dòng)能大量損失,使其帶動(dòng)未被風(fēng)機(jī)吸入的低速風(fēng)流流動(dòng)的效果變差[11],這從流體運(yùn)動(dòng)學(xué)的機(jī)制上解釋了風(fēng)機(jī)布置得越高,其提供的升壓力越小這一結(jié)論。

    圖6隨風(fēng)機(jī)布置高度增加沿風(fēng)機(jī)軸線剖面速度變化云圖(單位: m/s)

    Fig. 6Nephograms of air velocity variation along fan axis vs. installing height of fan (m/s)

    2.2射流風(fēng)機(jī)橫向間距對通風(fēng)效果的影響

    JTG/T D70/2-02—2014《公路隧道通風(fēng)設(shè)計(jì)細(xì)則》[10]中規(guī)定,當(dāng)同一斷面布置2臺(tái)及以上射流風(fēng)機(jī)時(shí),相鄰風(fēng)機(jī)的橫向凈距不宜小于1倍風(fēng)機(jī)直徑。因此,為研究射流風(fēng)機(jī)橫向間距對通風(fēng)效果的影響,設(shè)定2臺(tái)風(fēng)機(jī)凈距1倍直徑、2臺(tái)風(fēng)機(jī)凈距1.5倍直徑、2臺(tái)風(fēng)機(jī)凈距2倍直徑、2臺(tái)風(fēng)機(jī)凈距2.5倍直徑4組工況進(jìn)行計(jì)算。

    第1組風(fēng)機(jī)距隧道入口50 m,第2組風(fēng)機(jī)距隧道入口200 m,氣流在到達(dá)各組射流風(fēng)機(jī)前隧道截面靜壓的最小值和氣流經(jīng)射流風(fēng)機(jī)升壓后隧道截面靜壓的最大值如圖7所示。

    圖7射流風(fēng)機(jī)前后截面最大及最小靜壓值(風(fēng)機(jī)橫向凈距不同時(shí))

    Fig. 7Maximum and minimum static pressures of tunnel cross-section(with different transverse distances among fans)

    風(fēng)機(jī)橫向間距對升壓力的影響見表3。由表3可知,當(dāng)2臺(tái)風(fēng)機(jī)橫向凈距在1~2倍直徑時(shí),風(fēng)機(jī)調(diào)壓綜合影響系數(shù)k值隨橫向凈距的增大而增大;但當(dāng)風(fēng)機(jī)橫向凈距為2.5倍直徑時(shí),k值則大大減小。因此,建議明堂山隧道在進(jìn)行風(fēng)機(jī)布置時(shí),風(fēng)機(jī)橫向凈距設(shè)為1.5~2倍風(fēng)機(jī)直徑。

    表3風(fēng)機(jī)橫向凈距對升壓力的影響

    Table 3Influence of transverse distance among fans on pressure rise

    風(fēng)機(jī)橫向凈距psj/Papj/Pak2臺(tái)風(fēng)機(jī)凈距1倍直徑47.5636358.870.8082臺(tái)風(fēng)機(jī)凈距1.5倍直徑50.7642458.870.8622臺(tái)風(fēng)機(jī)凈距2倍直徑55.7169158.870.9462臺(tái)風(fēng)機(jī)凈距2.5倍直徑41.6124058.870.707

    當(dāng)2臺(tái)風(fēng)機(jī)橫向凈距為2倍直徑時(shí),風(fēng)機(jī)前后縱向靜壓變化云圖如圖8和圖9所示。

    圖8風(fēng)機(jī)橫向凈距為2D時(shí)第1組射流風(fēng)機(jī)前后縱向靜壓變化云圖(單位: Pa)

    Fig. 8Nephogram of static pressure variation of Jet Fan Group No. 1 when clear distance among fans is 2D(Pa)

    圖9風(fēng)機(jī)橫向凈距為2D時(shí)第2組射流風(fēng)機(jī)前后縱向靜壓變化云圖(單位: Pa)

    Fig. 9Nephogram of static pressure variation of Jet Fan Group No. 2 when clear distance among fans is 2D(Pa)

    2臺(tái)風(fēng)機(jī)橫向凈距為1倍和2倍風(fēng)機(jī)直徑時(shí),風(fēng)機(jī)出口處速度矢量圖如圖10和圖11所示。對比圖10和圖11可知,風(fēng)機(jī)橫向凈距為1D時(shí),其出口附近風(fēng)速明顯要比風(fēng)機(jī)橫向凈距為2D時(shí)小。

    圖102臺(tái)風(fēng)機(jī)凈距為1倍風(fēng)機(jī)直徑時(shí)風(fēng)機(jī)出口處風(fēng)速矢量圖 (單位: m/s)

    Fig. 10Vectogram of air velocity at fan exit when clear distance among fans is 1D(m/s)

    隨風(fēng)機(jī)橫向間距增加沿風(fēng)機(jī)軸線剖面速度變化云圖如圖12所示。由圖12可知,當(dāng)風(fēng)機(jī)橫向凈距在1~2倍直徑時(shí),風(fēng)機(jī)出口附近風(fēng)速逐漸增大;而當(dāng)風(fēng)機(jī)橫向凈距為2.5倍直徑時(shí),其出口附近風(fēng)速大大減小,甚至比風(fēng)機(jī)凈距為1倍直徑時(shí)還要小。這是因?yàn)轱L(fēng)機(jī)橫向間距越小,2臺(tái)風(fēng)機(jī)出口風(fēng)流越易相互受到影響,即風(fēng)機(jī)噴射的高速風(fēng)流在風(fēng)機(jī)后方匯聚時(shí)相互卷吸和干擾,使本應(yīng)被射流卷吸進(jìn)來的周圍速度相對小的隧道氣流未能通過射流微團(tuán)產(chǎn)生的橫向脈動(dòng)得到動(dòng)量與能量交換[12],因此通風(fēng)效果變差。當(dāng)風(fēng)機(jī)橫向凈距為2.5倍直徑時(shí),雖然2臺(tái)風(fēng)機(jī)間距離足夠大,但每臺(tái)風(fēng)機(jī)與各自一側(cè)的隧道壁面距離變得很小,而風(fēng)機(jī)與隧道壁面空間越小,射流發(fā)展越受到限制,升壓作用越不明顯,所以此時(shí)k值變小。

    2.3射流風(fēng)機(jī)縱向間距對通風(fēng)效果的影響

    在射流通風(fēng)中,考慮到電纜線路的電力損耗、風(fēng)機(jī)管理和維護(hù)等因素,射流風(fēng)機(jī)一般集中布置。在集中布置中只有每組風(fēng)機(jī)之間保持足夠的距離,噴射氣流才會(huì)在到達(dá)下一組風(fēng)機(jī)位置之前充分減速。如果噴射氣流減速不完全,將會(huì)影響到下一級風(fēng)機(jī)的工作性能[13]。因此,為研究射流風(fēng)機(jī)縱向間距對通風(fēng)效果的影響,設(shè)定2組風(fēng)機(jī)縱向間距50 m、2組風(fēng)機(jī)縱向間距100 m、2組風(fēng)機(jī)縱向間距150 m、2組風(fēng)機(jī)縱向間距200 m 4組工況進(jìn)行計(jì)算。

    圖112臺(tái)風(fēng)機(jī)凈距為2倍風(fēng)機(jī)直徑時(shí)風(fēng)機(jī)出口處風(fēng)速矢量圖 (單位: m/s)

    Fig. 11Vectogram of air velocity at fan exit when clear distance among fans is 2D(m/s)

    圖12隨風(fēng)機(jī)橫向間距增加沿風(fēng)機(jī)軸線剖面速度變化云圖 (單位: m/s)

    Fig. 12Nephograms of air velocity variation along fan axis vs. transverse distance among fans (m/s)

    隨射流風(fēng)機(jī)縱向間距的增加沿風(fēng)機(jī)軸線剖面風(fēng)速變化云圖如圖13所示。由圖13可知,2組風(fēng)機(jī)間隧道部分的速度是先增大后減小,這是由于射流微團(tuán)的卷吸作用,將動(dòng)量和能量交換給隧道氣流,使射流范圍擴(kuò)展,流量增加,壓力上升形成射流發(fā)展過程,這時(shí)隧道內(nèi)速度是增大的;當(dāng)伴隨流動(dòng)范圍逐漸減少,整個(gè)隧道氣流沿著縱向呈現(xiàn)一種漸變的、非均勻的逆壓流動(dòng),直至射流發(fā)展完全,伴隨流消失,斷面形成均勻速度分布[14],這一階段隧道內(nèi)速度是減小的。當(dāng)2組風(fēng)機(jī)縱向間距為50 m和100 m時(shí),第2組風(fēng)機(jī)前的速度還沒有完全下降到隧道斷面平均速度大??;而當(dāng)2組風(fēng)機(jī)縱向間距為150 m和200 m時(shí),第2組風(fēng)機(jī)前的速度已基本下降到隧道斷面平均速度大小。表明隨著射流風(fēng)機(jī)縱向間距的增加,噴射氣流逐漸減速完全,射流得到更充分發(fā)展。因此,建議明堂山隧道在進(jìn)行風(fēng)機(jī)布置時(shí),將風(fēng)機(jī)縱向間距設(shè)在150 m以上。

    圖13隨射流風(fēng)機(jī)縱向間距的增加沿風(fēng)機(jī)軸線剖面風(fēng)速變化云圖 (單位: m/s)

    Fig. 13Nephograms of air velocity variation along fan axis vs. longitudinal distance among fans (m/s)

    3明堂山隧道左線全射流通風(fēng)方案數(shù)值模擬驗(yàn)證

    岳武高速明堂山特長公路隧道長7 500 m,分左右2條線。其中,左線采用全射流通風(fēng)方案,設(shè)計(jì)風(fēng)量為456 m3/s(即斷面平均風(fēng)速可達(dá)7 m/s),共布置了54臺(tái)1120型射流風(fēng)機(jī)。

    隧道實(shí)際風(fēng)機(jī)布置方式為: 風(fēng)機(jī)距建筑限界高度為20 cm,風(fēng)機(jī)橫向凈距為2倍風(fēng)機(jī)直徑,風(fēng)機(jī)縱向間距為170 m。

    現(xiàn)按隧道實(shí)際長度及設(shè)計(jì)射流風(fēng)機(jī)臺(tái)數(shù)建模,風(fēng)機(jī)空間布局采用明堂山隧道實(shí)際布置方案。模型總長5 000 m,第1組風(fēng)機(jī)距隧道入口200 m,每組風(fēng)機(jī)縱向間距170 m,最后一組風(fēng)機(jī)距隧道出口380 m。為確保數(shù)值模擬的準(zhǔn)確性,將距出口200 m處的隧道截面設(shè)置為局部阻力格柵(邊界條件為散熱器,損失系數(shù)經(jīng)大量試算后設(shè)為2.7),以保證5 000 m隧道模型的總阻力損失達(dá)到克服實(shí)際隧道長度阻力所需的升壓力值。隧道出入口和風(fēng)機(jī)出入口的邊界條件與2.2中設(shè)定的一致。

    計(jì)算結(jié)果表明,隧道入口截面的平均風(fēng)速為6.58 m/s,隧道出口截面的平均風(fēng)速為6.61 m/s,與設(shè)計(jì)值7 m/s相差很小,說明在這54臺(tái)射流風(fēng)機(jī)的作用下,隧道吸入的風(fēng)量與設(shè)計(jì)風(fēng)量基本一致。小于7 m/s是因?yàn)樗淼涝O(shè)計(jì)風(fēng)量是由阻滯工況時(shí)的需風(fēng)量決定的,而設(shè)計(jì)風(fēng)機(jī)臺(tái)數(shù)是由隧道通風(fēng)阻抗力、交通風(fēng)動(dòng)力和自然風(fēng)阻力共同決定的[15];但模擬時(shí)未考慮交通風(fēng)動(dòng)力的作用,導(dǎo)致在54臺(tái)射流風(fēng)機(jī)作用下隧道入口未能吸入達(dá)到設(shè)計(jì)值的風(fēng)量,所以數(shù)值模擬計(jì)算得到的隧道出入口風(fēng)速略小于設(shè)計(jì)值。風(fēng)機(jī)軸線高度處截面的平均縱向速度如圖14所示。

    圖14 風(fēng)機(jī)軸線高度處截面的平均縱向速度

    Fig. 14Diagram of average longitudinal velocity of cross-section at the axial height of fan

    4結(jié)論與建議

    本文依托明堂山隧道工程,針對射流風(fēng)機(jī)的空間布局優(yōu)化問題,開展了CFD數(shù)值模擬分析,得到以下結(jié)論。

    1)在2車道公路隧道斷面上布置射流風(fēng)機(jī)時(shí),風(fēng)機(jī)布置得越高,越接近隧道拱頂,風(fēng)機(jī)調(diào)壓綜合影響系數(shù)越小,即提供的升壓力越小。建議明堂山隧道在進(jìn)行風(fēng)機(jī)布置時(shí),將風(fēng)機(jī)設(shè)在距建筑限界15~30 cm高度處。

    2)在2車道公路隧道斷面上布置射流風(fēng)機(jī)時(shí),隨著2臺(tái)風(fēng)機(jī)橫向凈距的增加,風(fēng)機(jī)調(diào)壓綜合影響系數(shù)也在增加,即提供的升壓力也在增大,但離得過遠(yuǎn)又會(huì)使每臺(tái)風(fēng)機(jī)與各自一側(cè)隧道壁面的空間不足以致射流發(fā)展受限,造成升壓力下降。建議明堂山隧道在進(jìn)行風(fēng)機(jī)布置時(shí),將風(fēng)機(jī)橫向凈距設(shè)為1.5~2倍風(fēng)機(jī)直徑。

    3)在2車道公路隧道斷面上布置射流風(fēng)機(jī)時(shí),風(fēng)機(jī)縱向間距應(yīng)設(shè)置在150 m以上,以保證噴射氣流在到達(dá)下一組風(fēng)機(jī)位置之前充分減速,使射流得到充分發(fā)展。建議明堂山隧道在進(jìn)行風(fēng)機(jī)布置時(shí),將風(fēng)機(jī)縱向間距設(shè)在150 m以上。

    隧道實(shí)際風(fēng)機(jī)布置方式所采用的參數(shù)均在優(yōu)化結(jié)果范圍內(nèi),可見本文對2車道公路隧道射流風(fēng)機(jī)空間布局優(yōu)化問題所做的CFD數(shù)值分析具有指導(dǎo)意義。

    本文在對明堂山隧道實(shí)際長度和設(shè)計(jì)射流風(fēng)機(jī)臺(tái)數(shù)建模時(shí),未考慮交通風(fēng)動(dòng)力和自然風(fēng)阻力的影響,僅做了簡化處理,在后續(xù)研究中,期望借助Fluent動(dòng)網(wǎng)格技術(shù)考慮交通風(fēng)的影響,使數(shù)值模擬更加貼合實(shí)際。

    參考文獻(xiàn)(References):

    [1]張林.公路隧道射流通風(fēng)數(shù)值模擬與風(fēng)機(jī)安裝參數(shù)優(yōu)化[D].重慶: 重慶交通大學(xué),2009.(ZHANG Lin. Numerical simulation of jet ventilation and optimization of installation parameters of fan of road tunnel[D].Chongqing: Chongqing Jiaotong University, 2009.(in Chinese))

    [2]楊秀軍,王曉雯,陳建忠.公路隧道通風(fēng)中射流風(fēng)機(jī)縱向最小間距研究[J].重慶交通大學(xué)學(xué)報(bào)(自然科學(xué)版),2008,27(1): 40-43.(YANG Xiujun, WANG Xiaowen, CHEN Jianzhong. Study on minimum longitudinal distance of jet fan in ventilation of highway tunnel[J]. Journal of Chongqing Jiaotong University(Natural Science), 2008, 27(1): 40-43.(in Chinese))

    [3]胡純良,張海琳,王敏,等.不同風(fēng)機(jī)軸距下某隧道內(nèi)氣流分布研究[J].四川建筑科學(xué)研究,2014,40(1): 94-97.(HU Chunliang, ZHANG Hailin, WANG Min, et al . Study on air flow distribution in tunnel under different axial distances among fans [J].Sichuan Building Science,2014,40(1): 94-97. (in Chinese))

    [4]方勇,蘭宇,曾艷華,等. 三車道公路隧道射流風(fēng)機(jī)設(shè)置位置研究[J].現(xiàn)代隧道技術(shù),2009,46(2): 90-93.(FANG Yong, LAN Yu, ZENG Yanhua, et al. Study on the installation of jet fan in three-lane highway tunnel[J].Modern Tunnelling Technology,2009,46(2): 90-93.(in Chinese))

    [5]石平.公路隧道通風(fēng)局部效應(yīng)三維數(shù)值模擬分析與研究[D].西安: 長安大學(xué),2004.(SHI Ping. 3D numerical simulation analysis and research on local effect of road tunnel ventilation[D]. Xi’an: Chang’an University, 2004. (in Chinese))

    [6]王永東.長大公路隧道縱向通風(fēng)數(shù)值模擬研究[D].西安: 長安大學(xué),2000.(WANG Yongdong. Numerical simulation research on longitudinal ventilation of long and large road tunnel[D].Xi’an: Chang’an University,2000.(in Chinese))

    [7]孫詳海.流體力學(xué)[M].上海:上海交通大學(xué)出版社,2000: 23-35.(SUN Xianghai. Fluid mechanics[M]. Shanghai: Shanghai Jiaotong University Publisher,2000: 23-35.(in Chinese))

    [8]王福軍.計(jì)算流體動(dòng)力學(xué)分析:CFD軟件原理與應(yīng)用[M].北京: 清華大學(xué)出版社,2004:135-180.(WANG Fujun. Computational fluid dynamics analysis:Theory and application of CFD software[M]. Beijing: Tsinghua University Publisher,2004: 135-180.(in Chinese))

    [9]趙學(xué)錄.長大隧道通風(fēng)與防災(zāi)技術(shù)研究[J].隧道建設(shè),2007,27(增刊2): 171-173.(ZHAO Xuelu. Research on techniques of ventilation and disaster prevention in long tunnels[J]. Tunnel Construction,2007,27(S2): 171-173.(in Chinese))

    [10]公路隧道通風(fēng)設(shè)計(jì)細(xì)則: JTG/T D70/2-02—2014 [S].北京: 人民交通出版社,2014.(Guidelines for design of ventilation of highway tunnels:JTG/T D70/2-02—2014[S].Beijing: China Communications Press,2014.(in Chinese))

    [11]Axel Bring, Tor-Goran Malmstrom, Carl Axel Boman. Simulation and measurement of road tunnel ventilation[J]. Tunnelling and Underground Space Technology, 1997,12(3): 417-424.

    [12]Bari S, Naser J. Simulation of airflow and pollution levels caused by severe traffic jam in a road tunnel[J]. Tunnelling and Underground Space Technology,2010,25(1): 70-77.

    [13]Wongwises S, Songneti T, Lokathada N, et al. A comparison of the flow characteristics of refrigerants flowing through adiabatic capillary tubes[J]. International Communications in Heat and Mass Transfer,2011,38(3): 398-404.

    [14]Gonzalez J A. A practical approach for evaluating ventilation air duct performance in an existing road tunnel[C]//International Symposium on Aerodynamics and Ventilation of Vehicle Tunnels. Luzern: BHR Group, 2003, 11 (2): 429-448.

    [15]Mizuno A, Azuma T, Lchikawa A. Appropriate JF installation interval obtained by booster capacity measurements[C]//BHR Group Conference Series Publication. Bury St. Edmunds: Professional Engineering Publishing, 2000: 431-448.

    Computational Fluid Dynamics (CFD) Analysis of Spatial Layout Optimization of Jet Ventilation Fan Used in Two-lane Road Tunnel

    ZHAO Li1, YAN Zhiguo1,2,3

    (1.DepartmentofGeotechnicalEngineering,TongjiUniversity,Shanghai200092,China; 2.StateKeyLaboratoryforDisasterReductioninCivilEngineering,TongjiUniversity,Shanghai200092,China; 3.KeyLaboratoryofGeotechnicalandUndergroundEngineeringofMinistryofEducation,TongjiUniversity,Shanghai200092,China)

    Abstract:The jet ventilation of Mingtangshan Two-lane Road Tunnel is numerically simulated so as to optimize the spatial layout of ventilation fans; and then the factors for the pressure rise of ventilation fan, i.e. longitudinal distance among fans, setting height of fan and transverse distance among fans, are discussed. The authors suggest that the jet fans should be installed at the position higher than building boundary about 15-30 cm; the transverse distance among fans should be 1.5-2 times fan diameter; and the longitudinal distance among fans should be larger than 150 m. The computational fluid dynamics (CFD) simulation results show that the jet ventilation fan after optimization can achieve good effect.

    Keywords:Mingtangshan Road Tunnel; jet ventilation fan; spatial layout optimization; numerical simulation; computational fluid dynamics (CFD)

    中圖分類號:U 453.5

    文獻(xiàn)標(biāo)志碼:B

    文章編號:1672-741X(2016)04-0411-07

    DOI:10.3973/j.issn.1672-741X.2016.04.007

    第一作者簡介:趙黎(1991—),男,河北秦皇島人,同濟(jì)大學(xué)隧道與地下建筑工程專業(yè)在讀碩士,研究方向?yàn)樗淼劳L(fēng)。E-mail: 1009448957@qq.com。

    基金項(xiàng)目:山區(qū)隧道防災(zāi)、資源利用及環(huán)保技術(shù)集成研究(2013318J02120)

    收稿日期:2015-09-28; 修回日期: 2015-12-09

    午夜福利在线免费观看网站| 人人妻人人添人人爽欧美一区卜| 亚洲熟女精品中文字幕| 国产熟女午夜一区二区三区| 久久久久精品久久久久真实原创| 国产国语露脸激情在线看| 少妇人妻久久综合中文| 99热网站在线观看| 国产精品免费大片| 黑人欧美特级aaaaaa片| 热re99久久国产66热| 超碰成人久久| 丝袜脚勾引网站| 新久久久久国产一级毛片| 国产一区二区激情短视频 | 咕卡用的链子| 久久久a久久爽久久v久久| 99国产精品免费福利视频| 母亲3免费完整高清在线观看 | 香蕉精品网在线| a级片在线免费高清观看视频| 人成视频在线观看免费观看| 国产精品三级大全| 不卡av一区二区三区| 超色免费av| 久久久久精品久久久久真实原创| 又粗又硬又长又爽又黄的视频| 2018国产大陆天天弄谢| 黄色怎么调成土黄色| 9191精品国产免费久久| 老司机亚洲免费影院| 午夜精品国产一区二区电影| 天天躁狠狠躁夜夜躁狠狠躁| 1024香蕉在线观看| 亚洲欧美精品自产自拍| 亚洲精品久久久久久婷婷小说| 欧美bdsm另类| 日日撸夜夜添| 欧美中文综合在线视频| 午夜免费鲁丝| 亚洲欧美一区二区三区久久| 不卡av一区二区三区| 另类精品久久| 欧美av亚洲av综合av国产av | 黄色配什么色好看| av.在线天堂| 亚洲激情五月婷婷啪啪| 少妇人妻 视频| 男女国产视频网站| 人妻人人澡人人爽人人| 国产精品蜜桃在线观看| videosex国产| 男人操女人黄网站| 日韩制服骚丝袜av| 天堂中文最新版在线下载| 侵犯人妻中文字幕一二三四区| 黑人巨大精品欧美一区二区蜜桃| 美女脱内裤让男人舔精品视频| 亚洲综合色惰| 在现免费观看毛片| 久久毛片免费看一区二区三区| 久久免费观看电影| 欧美精品av麻豆av| 黄色 视频免费看| 久久久久国产网址| 亚洲综合色惰| 亚洲精品日本国产第一区| 国产色婷婷99| freevideosex欧美| 啦啦啦在线免费观看视频4| 中文字幕亚洲精品专区| 欧美激情极品国产一区二区三区| 国产成人91sexporn| 精品久久久精品久久久| av在线app专区| 久久国产精品大桥未久av| 国产麻豆69| 国产1区2区3区精品| 久久久久久人人人人人| 精品国产一区二区三区久久久樱花| 国产一区二区在线观看av| av卡一久久| 久久精品国产a三级三级三级| 大陆偷拍与自拍| 久久久久久伊人网av| 国语对白做爰xxxⅹ性视频网站| 在现免费观看毛片| 亚洲人成电影观看| 丰满迷人的少妇在线观看| 日韩一区二区三区影片| 国产免费现黄频在线看| 又粗又硬又长又爽又黄的视频| a 毛片基地| 亚洲人成网站在线观看播放| 久久精品亚洲av国产电影网| 黄色怎么调成土黄色| 国产爽快片一区二区三区| 久久女婷五月综合色啪小说| 在线观看三级黄色| 韩国av在线不卡| 欧美人与善性xxx| 婷婷色综合大香蕉| 巨乳人妻的诱惑在线观看| 国产成人精品久久二区二区91 | 男人操女人黄网站| 春色校园在线视频观看| 亚洲精华国产精华液的使用体验| 最黄视频免费看| 大片免费播放器 马上看| 三上悠亚av全集在线观看| 日韩,欧美,国产一区二区三区| 欧美人与性动交α欧美精品济南到 | 久久韩国三级中文字幕| 日本91视频免费播放| 久久人妻熟女aⅴ| 国产精品三级大全| 狠狠精品人妻久久久久久综合| 成人亚洲欧美一区二区av| 国产成人免费观看mmmm| 欧美日韩精品网址| 国产欧美亚洲国产| 爱豆传媒免费全集在线观看| 国产精品一二三区在线看| www.av在线官网国产| 最黄视频免费看| 精品久久久精品久久久| h视频一区二区三区| 免费高清在线观看视频在线观看| 国产精品久久久久久av不卡| 97在线人人人人妻| 久久99精品国语久久久| 大片免费播放器 马上看| 国产成人精品一,二区| 国产成人精品无人区| 久久久精品94久久精品| 久久国内精品自在自线图片| 18禁动态无遮挡网站| 18禁国产床啪视频网站| 中文字幕人妻丝袜一区二区 | 天天躁夜夜躁狠狠久久av| 久久久精品国产亚洲av高清涩受| 国产一区亚洲一区在线观看| 免费观看性生交大片5| 国产精品久久久久久精品电影小说| 少妇 在线观看| 国产日韩欧美亚洲二区| 日韩,欧美,国产一区二区三区| 男女啪啪激烈高潮av片| 欧美日韩亚洲国产一区二区在线观看 | 人人澡人人妻人| 久久热在线av| 狂野欧美激情性bbbbbb| 久久精品熟女亚洲av麻豆精品| 免费高清在线观看日韩| 亚洲人成电影观看| 男女高潮啪啪啪动态图| 国产 精品1| 国产乱人偷精品视频| 黄色一级大片看看| 99香蕉大伊视频| 久久人人爽av亚洲精品天堂| 亚洲精品在线美女| 伦理电影大哥的女人| 国产成人精品一,二区| 国产免费福利视频在线观看| 两个人看的免费小视频| 国产野战对白在线观看| av网站免费在线观看视频| 色视频在线一区二区三区| 亚洲美女搞黄在线观看| 视频在线观看一区二区三区| av国产久精品久网站免费入址| 精品视频人人做人人爽| 国产欧美日韩综合在线一区二区| 日韩 亚洲 欧美在线| 亚洲在久久综合| 亚洲综合色惰| 美女视频免费永久观看网站| 男的添女的下面高潮视频| 看免费av毛片| 免费少妇av软件| 婷婷色麻豆天堂久久| 免费观看无遮挡的男女| 日韩欧美一区视频在线观看| 国产精品秋霞免费鲁丝片| 女性被躁到高潮视频| 十八禁网站网址无遮挡| 免费久久久久久久精品成人欧美视频| 亚洲精品中文字幕在线视频| 精品人妻偷拍中文字幕| videossex国产| av视频免费观看在线观看| 国产熟女午夜一区二区三区| 最近的中文字幕免费完整| 国产精品.久久久| 久久久国产一区二区| 热re99久久精品国产66热6| 最近中文字幕2019免费版| 高清av免费在线| 老司机影院成人| 天天躁日日躁夜夜躁夜夜| 两性夫妻黄色片| 国产97色在线日韩免费| 午夜日韩欧美国产| 激情视频va一区二区三区| 秋霞伦理黄片| 边亲边吃奶的免费视频| 亚洲国产欧美在线一区| 国语对白做爰xxxⅹ性视频网站| 国产日韩欧美亚洲二区| 在线天堂最新版资源| 七月丁香在线播放| 一区二区三区四区激情视频| 久久久久久免费高清国产稀缺| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产黄频视频在线观看| 在线观看国产h片| 久久精品夜色国产| 亚洲国产精品999| 日本av免费视频播放| 嫩草影院入口| 午夜影院在线不卡| 午夜福利网站1000一区二区三区| 午夜91福利影院| 考比视频在线观看| 久久久久视频综合| 妹子高潮喷水视频| 日韩电影二区| 亚洲国产av新网站| 丝袜在线中文字幕| 中文天堂在线官网| 国产在视频线精品| 女性生殖器流出的白浆| 人人妻人人添人人爽欧美一区卜| 欧美97在线视频| 99热国产这里只有精品6| av女优亚洲男人天堂| 岛国毛片在线播放| 亚洲国产精品一区三区| 精品一品国产午夜福利视频| 国产精品亚洲av一区麻豆 | 亚洲男人天堂网一区| 日韩中文字幕视频在线看片| 18禁动态无遮挡网站| 久久免费观看电影| 亚洲婷婷狠狠爱综合网| 日本黄色日本黄色录像| 搡老乐熟女国产| 在线观看人妻少妇| 又大又黄又爽视频免费| 一边摸一边做爽爽视频免费| 男女下面插进去视频免费观看| 久久韩国三级中文字幕| 精品人妻熟女毛片av久久网站| 亚洲精品国产一区二区精华液| 久久久久久久国产电影| av免费观看日本| 最近的中文字幕免费完整| 国产精品.久久久| 久久影院123| 欧美av亚洲av综合av国产av | 午夜福利视频精品| 久久毛片免费看一区二区三区| 成人国语在线视频| 亚洲少妇的诱惑av| 亚洲美女视频黄频| 另类精品久久| 亚洲精品,欧美精品| 精品少妇内射三级| 精品亚洲成a人片在线观看| 一本大道久久a久久精品| 久久久久久久国产电影| 久久国产精品大桥未久av| 中国三级夫妇交换| 午夜老司机福利剧场| 国产成人精品一,二区| 天天躁夜夜躁狠狠久久av| 欧美日韩视频高清一区二区三区二| 晚上一个人看的免费电影| 成人免费观看视频高清| 十分钟在线观看高清视频www| 好男人视频免费观看在线| 国产乱人偷精品视频| 国产成人精品久久二区二区91 | 天天操日日干夜夜撸| 日日撸夜夜添| 男人爽女人下面视频在线观看| 午夜免费男女啪啪视频观看| 欧美人与性动交α欧美精品济南到 | 老司机影院成人| 夫妻性生交免费视频一级片| 亚洲,一卡二卡三卡| 熟女av电影| 国产男女内射视频| 久久婷婷青草| 七月丁香在线播放| 精品亚洲成国产av| 七月丁香在线播放| 国产日韩欧美亚洲二区| 亚洲av男天堂| 99久久人妻综合| 五月伊人婷婷丁香| 婷婷成人精品国产| av免费观看日本| 激情视频va一区二区三区| 一级,二级,三级黄色视频| 我要看黄色一级片免费的| 波野结衣二区三区在线| 99久久综合免费| 欧美激情极品国产一区二区三区| 亚洲,一卡二卡三卡| 亚洲男人天堂网一区| 91久久精品国产一区二区三区| 亚洲国产毛片av蜜桃av| 成人国产av品久久久| 蜜桃在线观看..| 亚洲国产欧美在线一区| 精品一区二区三卡| 又粗又硬又长又爽又黄的视频| 精品亚洲成a人片在线观看| 亚洲美女视频黄频| 精品亚洲乱码少妇综合久久| av片东京热男人的天堂| 女人高潮潮喷娇喘18禁视频| 99热网站在线观看| 午夜福利视频在线观看免费| 女人精品久久久久毛片| 女人高潮潮喷娇喘18禁视频| 少妇的逼水好多| 肉色欧美久久久久久久蜜桃| 欧美 亚洲 国产 日韩一| 国产成人一区二区在线| 中国国产av一级| 人妻 亚洲 视频| 久久久久国产一级毛片高清牌| 精品亚洲乱码少妇综合久久| 桃花免费在线播放| 国产一区二区在线观看av| 一级黄片播放器| 最近最新中文字幕大全免费视频 | 成人手机av| 国产在线免费精品| 亚洲欧美一区二区三区黑人 | 色婷婷久久久亚洲欧美| 欧美变态另类bdsm刘玥| 99香蕉大伊视频| 成年人午夜在线观看视频| 人人妻人人爽人人添夜夜欢视频| 最近最新中文字幕免费大全7| 欧美成人午夜精品| 国产精品一区二区在线观看99| www.熟女人妻精品国产| 国产成人欧美| 叶爱在线成人免费视频播放| 亚洲精品乱久久久久久| 麻豆精品久久久久久蜜桃| 90打野战视频偷拍视频| 国产激情久久老熟女| 韩国av在线不卡| 国产精品一区二区在线观看99| 啦啦啦中文免费视频观看日本| 宅男免费午夜| 亚洲情色 制服丝袜| 国产精品不卡视频一区二区| 丰满迷人的少妇在线观看| 18+在线观看网站| 男女边吃奶边做爰视频| 久久精品aⅴ一区二区三区四区 | 久久精品国产鲁丝片午夜精品| 婷婷色综合大香蕉| 国产在线一区二区三区精| 欧美精品一区二区免费开放| 妹子高潮喷水视频| 最新中文字幕久久久久| 久久精品国产a三级三级三级| 国产精品麻豆人妻色哟哟久久| 日韩制服骚丝袜av| av在线老鸭窝| 各种免费的搞黄视频| 免费黄频网站在线观看国产| 高清av免费在线| 午夜福利乱码中文字幕| av一本久久久久| 成人午夜精彩视频在线观看| 国产高清国产精品国产三级| 十八禁网站网址无遮挡| 精品少妇内射三级| 久久人妻熟女aⅴ| 观看av在线不卡| 国产一级毛片在线| 久久久久精品久久久久真实原创| 精品人妻一区二区三区麻豆| 中文精品一卡2卡3卡4更新| 国产精品久久久av美女十八| √禁漫天堂资源中文www| 热99久久久久精品小说推荐| 校园人妻丝袜中文字幕| 波野结衣二区三区在线| 不卡av一区二区三区| 午夜福利网站1000一区二区三区| 亚洲精品成人av观看孕妇| 国产精品国产三级专区第一集| 日韩一区二区三区影片| 欧美日韩一级在线毛片| 丁香六月天网| 国产精品三级大全| 欧美人与善性xxx| 啦啦啦在线免费观看视频4| 国产免费现黄频在线看| 日本欧美国产在线视频| av网站免费在线观看视频| 精品少妇一区二区三区视频日本电影 | 精品国产乱码久久久久久男人| 国产精品久久久久成人av| 制服丝袜香蕉在线| 秋霞在线观看毛片| 最黄视频免费看| av线在线观看网站| 中文精品一卡2卡3卡4更新| 日韩欧美一区视频在线观看| 天堂8中文在线网| 大香蕉久久网| 国产成人免费无遮挡视频| 大陆偷拍与自拍| 亚洲四区av| 国产精品 国内视频| 亚洲精品成人av观看孕妇| 制服诱惑二区| 亚洲图色成人| 中文精品一卡2卡3卡4更新| 午夜福利一区二区在线看| 美女xxoo啪啪120秒动态图| 亚洲在久久综合| 精品少妇内射三级| 国产成人精品久久二区二区91 | 成年av动漫网址| av片东京热男人的天堂| 日韩一区二区视频免费看| 丰满迷人的少妇在线观看| 黄色一级大片看看| 2021少妇久久久久久久久久久| 纵有疾风起免费观看全集完整版| 国产日韩欧美视频二区| 欧美日韩一级在线毛片| 老汉色∧v一级毛片| av网站免费在线观看视频| 国产老妇伦熟女老妇高清| 亚洲人成77777在线视频| 免费观看无遮挡的男女| 国产精品不卡视频一区二区| 欧美中文综合在线视频| 精品亚洲成a人片在线观看| 99香蕉大伊视频| 国产精品久久久久久精品古装| 亚洲av中文av极速乱| 老司机影院成人| 中文字幕人妻丝袜制服| 黄色毛片三级朝国网站| 捣出白浆h1v1| 国产一区二区 视频在线| 另类亚洲欧美激情| 2018国产大陆天天弄谢| 欧美av亚洲av综合av国产av | √禁漫天堂资源中文www| 久久国产精品男人的天堂亚洲| 人人妻人人添人人爽欧美一区卜| 国产深夜福利视频在线观看| 亚洲av中文av极速乱| a 毛片基地| 成人毛片60女人毛片免费| 性色avwww在线观看| 国产男女超爽视频在线观看| 啦啦啦中文免费视频观看日本| 国产毛片在线视频| 黑人猛操日本美女一级片| 成年女人在线观看亚洲视频| 777米奇影视久久| 久久久久久久亚洲中文字幕| 一本—道久久a久久精品蜜桃钙片| 成年av动漫网址| 日韩伦理黄色片| 亚洲五月色婷婷综合| 一边摸一边做爽爽视频免费| 成人黄色视频免费在线看| 国产日韩欧美亚洲二区| 国产视频首页在线观看| 看十八女毛片水多多多| 日韩免费高清中文字幕av| 欧美日韩综合久久久久久| 青春草亚洲视频在线观看| 王馨瑶露胸无遮挡在线观看| 少妇被粗大的猛进出69影院| 久久99蜜桃精品久久| 亚洲欧洲日产国产| 成年动漫av网址| 免费av中文字幕在线| 欧美激情 高清一区二区三区| 国产视频首页在线观看| 久久精品aⅴ一区二区三区四区 | 人人妻人人添人人爽欧美一区卜| 赤兔流量卡办理| √禁漫天堂资源中文www| 秋霞伦理黄片| 男女边吃奶边做爰视频| 日韩一区二区视频免费看| 多毛熟女@视频| 两个人免费观看高清视频| 多毛熟女@视频| 免费看不卡的av| 亚洲精品国产色婷婷电影| 80岁老熟妇乱子伦牲交| 国产片内射在线| 黄片小视频在线播放| 久久综合国产亚洲精品| 国产精品国产av在线观看| 丰满乱子伦码专区| 亚洲人成网站在线观看播放| 国产视频首页在线观看| 一区福利在线观看| 男女高潮啪啪啪动态图| a级片在线免费高清观看视频| 国产精品女同一区二区软件| 高清视频免费观看一区二区| 天美传媒精品一区二区| 捣出白浆h1v1| 午夜福利视频精品| 国产男女超爽视频在线观看| 亚洲精品国产av成人精品| 久久精品aⅴ一区二区三区四区 | 日韩制服骚丝袜av| 亚洲av在线观看美女高潮| 精品一区二区免费观看| 亚洲第一av免费看| 性少妇av在线| 免费av中文字幕在线| 人妻 亚洲 视频| 97在线视频观看| a级片在线免费高清观看视频| 青青草视频在线视频观看| 99re6热这里在线精品视频| 免费观看性生交大片5| 一本—道久久a久久精品蜜桃钙片| 97在线人人人人妻| 各种免费的搞黄视频| 亚洲国产看品久久| 国产精品免费大片| 久久这里有精品视频免费| 一级毛片电影观看| 久热久热在线精品观看| 亚洲成色77777| 一级毛片我不卡| 久久久久久久久久久免费av| 人妻人人澡人人爽人人| 亚洲精品aⅴ在线观看| 国产精品久久久久久精品电影小说| 国产一区二区三区综合在线观看| 久久久久久久久久久久大奶| 建设人人有责人人尽责人人享有的| 视频在线观看一区二区三区| 少妇 在线观看| 国产精品一国产av| 亚洲精品国产av蜜桃| 免费在线观看黄色视频的| 国产野战对白在线观看| 少妇的丰满在线观看| 三上悠亚av全集在线观看| 女人精品久久久久毛片| 男女无遮挡免费网站观看| www.精华液| 久久精品aⅴ一区二区三区四区 | 午夜老司机福利剧场| 天天操日日干夜夜撸| 最近手机中文字幕大全| 在线 av 中文字幕| 久久久久久久亚洲中文字幕| 日韩中文字幕视频在线看片| 在线观看免费视频网站a站| 多毛熟女@视频| 深夜精品福利| 91在线精品国自产拍蜜月| 赤兔流量卡办理| 精品少妇久久久久久888优播| 欧美最新免费一区二区三区| 黄色怎么调成土黄色| 欧美日韩精品网址| 熟女电影av网| 美女主播在线视频| 国产深夜福利视频在线观看| 国产成人免费观看mmmm| 久久久a久久爽久久v久久| 成人国语在线视频| 两个人免费观看高清视频| 亚洲情色 制服丝袜| 成人亚洲精品一区在线观看| 久久国内精品自在自线图片| 1024视频免费在线观看| 如何舔出高潮| 久久久久视频综合| 亚洲精品一二三| 91精品伊人久久大香线蕉| 国产免费福利视频在线观看| 人妻系列 视频| 国产极品粉嫩免费观看在线| 久久97久久精品| 秋霞伦理黄片| 综合色丁香网| 欧美激情高清一区二区三区 | 精品卡一卡二卡四卡免费| 久久久国产精品麻豆| 热99国产精品久久久久久7| 亚洲国产看品久久| 中文欧美无线码| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩亚洲国产一区二区在线观看 | 国产精品久久久久成人av|