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

    真空管道磁浮交通車體熱壓載荷分布特征及其非定常特性

    2023-03-21 01:41:02胡嘯馬天昊王瀟飛鄧自剛張繼旺張衛(wèi)華
    實(shí)驗(yàn)流體力學(xué) 2023年1期
    關(guān)鍵詞:尾車熱壓激波

    胡嘯,馬天昊,王瀟飛,鄧自剛,*,張繼旺,張衛(wèi)華

    1.西南交通大學(xué) 牽引動(dòng)力國家重點(diǎn)實(shí)驗(yàn)室,成都 610031

    2.西南交通大學(xué) 超高速真空管道磁浮交通研究中心,成都 610031

    3.西南交通大學(xué) 力學(xué)與航空航天學(xué)院,成都 610031

    0 引言

    超高速磁浮列車具有綠色、高速、輕量化等優(yōu)點(diǎn),是未來軌道交通的重要發(fā)展方向之一,相關(guān)技術(shù)研究也是當(dāng)下熱議話題[1]。若將磁浮列車置于低真空管道內(nèi)運(yùn)行,將大幅度降低氣動(dòng)阻力,減少氣動(dòng)噪聲污染,使得磁浮列車速度有望達(dá)到跨聲速甚至超聲速[2-3]。從2019年開始,國家陸續(xù)出臺(tái)了《交通強(qiáng)國建設(shè)綱要》和《國家綜合立體交通網(wǎng)規(guī)劃綱要》等系列文件,支持超高速磁浮列車的發(fā)展。600 km/h常導(dǎo)磁浮交通系統(tǒng)以及高溫超導(dǎo)高速磁浮工程化樣車及試驗(yàn)線的下線啟用,標(biāo)志著我國高速磁浮領(lǐng)域取得了跨越式進(jìn)展。為了開展高速試驗(yàn),2020年西南交通大學(xué)開始建設(shè)大科學(xué)裝置“多態(tài)耦合軌道交通動(dòng)模型試驗(yàn)平臺(tái)”(最高試驗(yàn)速度1 500 km/h 的超高速真空管道磁浮交通試驗(yàn)系統(tǒng))。2021年,中北大學(xué)與中國航天科工集團(tuán)第三研究院采用超導(dǎo)電動(dòng)磁浮技術(shù)聯(lián)合建設(shè)高速飛車大同(陽高)試驗(yàn)線,目標(biāo)速度為1 000 km/h。

    國內(nèi)外學(xué)者對(duì)真空管道磁浮交通相關(guān)的基礎(chǔ)科學(xué)問題也開展了大量的研究。盡管磁浮列車在低真空管道內(nèi)運(yùn)行,減阻效果明顯,氣動(dòng)阻力值與管內(nèi)氣壓呈線性關(guān)系[4-7],但與空氣動(dòng)力學(xué)相關(guān)的問題仍然存在。由于列車兩端為流線型頭型,運(yùn)動(dòng)的列車與管道系統(tǒng)構(gòu)成動(dòng)態(tài)的拉伐爾噴管[8-10](C–D 噴管)。根據(jù)等熵極限與Kantrowitz 極限,給定某個(gè)列車速度,存在著一個(gè)臨界車/管阻塞比,當(dāng)實(shí)際阻塞比大于該值時(shí),管內(nèi)氣流將進(jìn)入雍塞狀態(tài),在尾流區(qū)產(chǎn)生激波,并在管道壁面反射傳播,惡化列車運(yùn)行環(huán)境[7-8,11-12],這與進(jìn)氣道中啟動(dòng)/不啟動(dòng)問題相似[13]。

    Li 等[14]基于一維無黏等熵理論,推導(dǎo)了亞聲速管道列車在雍塞/非雍塞狀態(tài)下的氣流參數(shù),與仿真結(jié)果吻合較好,同時(shí)Jang 等[15]對(duì)比了這兩種方法對(duì)氣動(dòng)阻力的預(yù)測結(jié)果。侯自豪等[16-17]采用準(zhǔn)一維數(shù)值模擬方法計(jì)算了3 種流動(dòng)狀態(tài)(亞聲速通流、壅塞和超聲速通流)下特長管道內(nèi)部流場特性,研究了管壁摩擦和列車加速對(duì)激波運(yùn)動(dòng)的影響。Yu 等[18]以質(zhì)量流量為指標(biāo),理論分析了管道雍塞的相關(guān)規(guī)律和影響因素,同時(shí)基于二維軸對(duì)稱模型,利用CFD 仿真分析了管道列車在雍塞/非雍塞狀態(tài)下的氣動(dòng)熱特性[19-20]。Zhong 等[21]考慮了車軌間隙,建立了三維車/管模型,重點(diǎn)比較了管道列車在雍塞/非雍塞狀態(tài)下的尾流演化,展示了尾渦結(jié)構(gòu)和三維激波形態(tài),與二維軸對(duì)稱模型的計(jì)算結(jié)果有著顯著差異,Hu 等[22]也得出類似結(jié)論。為了緩解雍塞效應(yīng),提高列車運(yùn)行效率,國內(nèi)外學(xué)者從列車外形[23-24]和管道結(jié)構(gòu)[9,25]兩方面開展了優(yōu)化設(shè)計(jì)研究。

    綜上,管內(nèi)雍塞狀態(tài)對(duì)列車運(yùn)行環(huán)境帶來了嚴(yán)峻挑戰(zhàn),目前的研究主要基于低維模型以及對(duì)稱模型對(duì)管內(nèi)流場的預(yù)測,鮮有研究關(guān)注三維列車幾何表面熱壓載荷的分布特性[21,26]及其非定常特性,而車軌間隙流場[11]、尾渦脫落[22]和激波與邊界層干擾[27]均會(huì)造成車體載荷波動(dòng),帶來列車壁板氣動(dòng)彈性問題[28]以及蒙皮材料的氣動(dòng)疲勞損傷[29],從而威脅列車運(yùn)行安全性、影響乘客舒適性。

    鑒于此,本文建立高溫超導(dǎo)磁浮列車三維幾何模型,利用IDDES 模型求解管道列車周圍流場,探明雍塞/非雍塞狀態(tài)下列車表面熱壓載荷的分布特征;此外,由于通過傳統(tǒng)的CFD 結(jié)果難以厘清列車表面載荷的全局非定常特性[30],本文利用本征正交分解(Proper Orthogonal Decomposition,POD)提取流場重要相干結(jié)構(gòu),識(shí)別列車表面載荷非定常較強(qiáng)區(qū)域,揭示列車表面載荷時(shí)空演化規(guī)律。研究結(jié)果可為管道列車車體結(jié)構(gòu)設(shè)計(jì)提供參考,助力真空管道交通模型試驗(yàn)平臺(tái)的建設(shè)。

    1 管內(nèi)流動(dòng)狀態(tài)判定

    1.1 等熵極限與Kantrowitz 極限

    車/管阻塞比和列車運(yùn)行速度共同決定管內(nèi)流動(dòng)狀態(tài)(雍塞與非雍塞狀態(tài)),一般來說,給定某個(gè)列車速度,可通過等熵極限(式(1))與Kantrowitz 極限(式(2))求出臨界阻塞比[31-32]。

    式中:β為車/管阻塞比,即列車車體橫截面積與管道凈空面積之比,下標(biāo)IL 和KL 分別代表等熵極限和Kantrowitz 極限計(jì)算出的臨界值;vtr為列車運(yùn)行速度;c 為聲速;γ為空氣比熱比,一般取1.4。

    根據(jù)式(1)和(2)繪制了管內(nèi)不同狀態(tài)下阻塞比與列車速度的對(duì)應(yīng)關(guān)系圖(圖1)。可以看到,等熵極限與Kantrowitz 極限將管內(nèi)流動(dòng)狀態(tài)分成3 個(gè)區(qū)域,分別是非雍塞區(qū)域、雍塞區(qū)域和雙解區(qū)域[16-17]。雙解區(qū)域位于超聲速等熵極限和 Kantrowitz 極限之間,當(dāng)列車運(yùn)行條件位于此區(qū)域時(shí),則管內(nèi)流動(dòng)存在兩種可能的解。

    1.2 計(jì)算工況

    目前真空管道磁浮交通的設(shè)計(jì)目標(biāo)時(shí)速為高亞聲速和跨聲速,本文固定列車速度為800 km/h。為了比較管道雍塞和非雍塞狀態(tài)下車體熱壓載荷的非定常特性,這里根據(jù)圖1 設(shè)計(jì)了A、B 和C 這3 個(gè)計(jì)算工況,其中工況A 和B 屬于雍塞類型,而工況C 為非雍塞類型,具體參數(shù)如表1所示。

    圖1 阻塞比與列車速度的臨界關(guān)系Fig.1 The critical relationship between blockage ratio and train speed

    表1 計(jì)算工況Table 1 Calculation case

    2 數(shù)值模型

    2.1 幾何模型

    如圖2所示,本研究中使用的列車模型是基于高溫超導(dǎo)磁浮系統(tǒng)特點(diǎn)建立的全尺寸高速磁浮樣車。安裝在懸浮架上的低溫容器(也稱杜瓦),提供保溫隔熱作用。置于杜瓦底部的超導(dǎo)體在液氮的冷卻作用下進(jìn)入超導(dǎo)狀態(tài),并與永磁軌道提供的磁場環(huán)境相互作用,在宏觀上產(chǎn)生與磁浮列車本身重力平衡的懸浮力,并提供橫向穩(wěn)定所需的導(dǎo)向力。高溫超導(dǎo)磁浮列車具有自懸浮、自導(dǎo)向、自穩(wěn)定特性[1]。

    圖2 真空管道磁浮交通系統(tǒng)幾何模型Fig.2 Geometric model of the evacuated tube maglev transportation system

    列車編組方式為三車編組,包括頭車、中間車和尾車,各車廂之間采用風(fēng)擋連接,每節(jié)車安裝兩個(gè)懸浮架,編號(hào)如圖2(a)所示,即頭車處為B1-1、B1-2,中間車處為B2-1、B2-2,尾車處為B3-2、B3-1。列車高度(Htr)為3.8 m,作為流場特征長度。列車流線型鼻長(Lsn)、車寬(Wtr)和列車總長(Ltr)分別為2.82Htr、0.89Htr和21.51Htr。此外,懸浮間隙(Hsg),即杜瓦底部和永磁軌道之間的垂向距離,為5.26×10?3Htr。低真空管道截面形狀參考高速鐵路的隧道斷面形狀,根據(jù)阻塞比大小確定工況A、B 和C 的管道凈空面積,其值分別為39.43、59.15 和118.29 m2。

    2.2 計(jì)算區(qū)域與邊界條件

    管道列車流動(dòng)具有內(nèi)外流場耦合特征,這與傳統(tǒng)高速列車的空氣動(dòng)力學(xué)特性有較大區(qū)別[33-34]。具體來說,在流向方向上需要足夠長的計(jì)算域來應(yīng)對(duì)由激波傳播引起的較大范圍干擾。根據(jù)初步的迭代計(jì)算結(jié)果,將列車放置在圖3所示的位置處(頭車距離入口130 Htr,尾車距離出口80 Htr),保證在采樣時(shí)間內(nèi)激波不會(huì)從計(jì)算域的兩端溢出,從而避免邊界處能量損失。

    模型中應(yīng)用了三種邊界條件。綠色虛線表示自由流邊界,也被稱為壓力遠(yuǎn)場,用于計(jì)算域的入口和出口,其中氣流的馬赫數(shù)、表壓和溫度分別設(shè)定為0.65、0 Pa 和288 K。為了模擬地面效應(yīng)和管道效應(yīng),地面、軌道和管道使用移動(dòng)壁面,即圖3 中的黑色實(shí)線,移動(dòng)壁面的切向速度大小指定為來流馬赫數(shù)。紅色實(shí)線表示靜止壁面,應(yīng)用于列車表面。

    圖3 計(jì)算區(qū)域與邊界條件示意圖Fig.3 Schematic diagram of the computational domain and boundary conditions

    此外,管內(nèi)的初始溫度(T0)和壓力(p0)分別為288 K、0.01 atm(1 013.25 Pa)。為了便于后續(xù)討論,坐標(biāo)系的原點(diǎn)O 固定在尾車鼻尖下的地面上,坐標(biāo)軸x、y、z 分別表示流向、展向和垂直方向。

    2.3 計(jì)算設(shè)置

    本文重點(diǎn)關(guān)注車體表面熱壓載荷的非定常特性,因此湍流的建模至關(guān)重要。本文采用IDDES 混合方法對(duì)湍流進(jìn)行求解,指定SSTk ?ω模型對(duì)邊界層處進(jìn)行建模,這種方法被廣泛應(yīng)用于高速列車[35-36]和真空管道交通[21-22]的氣動(dòng)特性研究中。

    考慮管內(nèi)氣流的強(qiáng)可壓縮性,采用STAR–CCM+軟件中的耦合隱式流體求解器。選擇AUSM+格式處理對(duì)流通量,可以準(zhǔn)確捕捉激波的不連續(xù)性。對(duì)于對(duì)流項(xiàng),選擇混合二階迎風(fēng)/有界中心方案,混合系數(shù)為0.15。二階隱式時(shí)間步進(jìn)格式用于離散時(shí)間項(xiàng),時(shí)間步長為0.013 tref(tref= Htr/ vtr),保證大部分流場單元的庫朗數(shù)小于1。當(dāng)計(jì)算時(shí)間到達(dá)46.78 tref后開始取樣計(jì)算平均量,此時(shí)流場的物理量沒有明顯的波動(dòng),除了由瞬態(tài)流場結(jié)構(gòu)引起的非定常振蕩。整個(gè)采樣時(shí)間持續(xù)到112.58 tref。

    2.4 網(wǎng)格劃分和獨(dú)立性驗(yàn)證

    利用STAR–CCM+軟件中Trimmer 網(wǎng)格和Prism layer 網(wǎng)格對(duì)計(jì)算區(qū)域進(jìn)行空間離散。為了捕捉車體表面熱壓載荷波動(dòng),對(duì)車體表面和尾流區(qū)域進(jìn)行多級(jí)網(wǎng)格加密,如圖4所示,列車表面最小網(wǎng)格尺寸為0.013 2Htr。為了求解列車壁面速度邊界層和熱邊界層,對(duì)列車表面劃分22 層棱柱層網(wǎng)格,拉伸比為1.2,第一層棱柱層網(wǎng)格厚度為5.26×10?5Htr。工況A、B 和C 的網(wǎng)格總數(shù)分別為2.44×107、2.62×107和2.91×107。

    圖4 計(jì)算網(wǎng)格加密方案Fig.4 Refinement scheme of the calculation grid

    雍塞狀態(tài)存在氣流阻塞以及尾部激波現(xiàn)象,流動(dòng)較為復(fù)雜。這里挑選工況A(β=0.3,雍塞)來驗(yàn)證網(wǎng)格獨(dú)立性,劃分3 套具有相同網(wǎng)格策略但密度不同的網(wǎng)格,詳細(xì)的網(wǎng)格參數(shù)如表2所示,其中y+為壁面第一層網(wǎng)格的無量綱高度。

    表2 用于網(wǎng)格獨(dú)立性研究的3 套網(wǎng)格分辨率Table 2 Three sets of grid resolutions for grid independence studies

    壓力系數(shù)Cp定義如下:

    式中:p為靜壓,ρ為空氣密度。

    圖5 比較了3 套網(wǎng)格下列車側(cè)表面平均熱壓載荷的計(jì)算結(jié)果。可以看到,3 套網(wǎng)格計(jì)算的熱壓載荷分布趨勢一致,主要差異在于局部流動(dòng)分離區(qū)域(如風(fēng)擋處、尾車流線型處),而這些局部流動(dòng)分離區(qū)往往是非定常特性較強(qiáng)區(qū)域,此外中網(wǎng)格和細(xì)網(wǎng)格的結(jié)果差異比粗網(wǎng)格和細(xì)網(wǎng)格的結(jié)果差異小得多。對(duì)比圖5(a)和(b)可以發(fā)現(xiàn),相較于壓力,溫度對(duì)網(wǎng)格尺寸更加敏感。

    圖5 不同網(wǎng)格密度對(duì)列車側(cè)面(z=0.53Htr)平均熱壓載荷的影響Fig.5 Effect of different grid densities on the time-averaged temperature and pressure load on the side of the train(z=0.53Htr)

    綜合考慮計(jì)算資源和計(jì)算準(zhǔn)確性,本文使用中等網(wǎng)格密度來評(píng)估真空管道磁浮交通車體熱壓載荷分布及其非定常特性。

    2.5 風(fēng)洞驗(yàn)證

    截至目前,尚未有公開的真空管道交通系統(tǒng)氣動(dòng)試驗(yàn)數(shù)據(jù),相關(guān)高速試驗(yàn)平臺(tái)正在建設(shè)中。本節(jié)利用Inger 等[37]在克蘭菲爾德大學(xué)航空學(xué)院跨聲速風(fēng)洞測試的試驗(yàn)數(shù)據(jù)來驗(yàn)證數(shù)值方法在模擬管內(nèi)雍塞狀態(tài)氣動(dòng)特性的可靠性。如圖6所示,風(fēng)洞試驗(yàn)底部平順,上表面安裝了半徑(R)0.58 m、高度0.02 m的圓弧凸塊,其截面變化構(gòu)成了一個(gè)拉瓦爾噴管,與本研究中的列車截面積變化非常類似。為與風(fēng)洞試驗(yàn)進(jìn)行對(duì)比,設(shè)置驗(yàn)證工況參數(shù)與試驗(yàn)一致(來流馬赫數(shù)0.73,初始溫度310 K),坐標(biāo)系為x'O'y',其他參數(shù)可參見文獻(xiàn)[37]。

    圖6 跨聲速風(fēng)洞試驗(yàn)段模型設(shè)置示意圖Fig.6 Schematic drawing of model setup in transonic wind tunnel test section

    圖7 比較了模型上表面壓力系數(shù)(Cp)的試驗(yàn)數(shù)據(jù)與數(shù)值仿真結(jié)果。從結(jié)果來看,兩者吻合良好,并且數(shù)值仿真對(duì)激波位置和強(qiáng)度的捕捉也較為準(zhǔn)確。觀察壓力系數(shù)隨截面變化的趨勢可以發(fā)現(xiàn):在收縮段(0 m

    圖7 上表面壓力系數(shù)分布對(duì)比Fig.7 Comparison of upper surface pressure coefficient distribution

    3 本征正交分解(POD)

    本征正交分解(POD)也被稱為主成分分析,是數(shù)據(jù)降維和模態(tài)分解中最普遍使用的方法之一,被廣泛用于識(shí)別列車尾流相干結(jié)構(gòu)[30,38-39]??偟膩碚f,POD 算法將流場矩陣分解為一系列相互正交的基函數(shù)和對(duì)應(yīng)的時(shí)間系數(shù),保證基函數(shù)在最小二乘意義上是最優(yōu)的,以捕獲盡可能多的能量?;瘮?shù)代表流場的空間模態(tài),時(shí)間系數(shù)則表示時(shí)間模態(tài),即POD 算法主要分析空間模態(tài)隨時(shí)間的演化過程,其具體分解過程如下:

    1)對(duì)一段時(shí)間區(qū)間[t1,tn]的空間流場進(jìn)行采樣,在ti(i=1,2,···,n)時(shí)刻,m 個(gè)空間網(wǎng)格點(diǎn)流場數(shù)據(jù)構(gòu)成的向量u(ti)=[u(x1,ti),u(x2,ti),···,u(xm,ti)]稱為一個(gè)流場快照??臻g網(wǎng)格點(diǎn)流場數(shù)據(jù)可以是瞬態(tài)數(shù)據(jù),也可以是瞬態(tài)數(shù)據(jù)減去平均值后的脈動(dòng)值,本文使用脈動(dòng)值進(jìn)行降階分析。

    2)將n 個(gè)流場快照數(shù)據(jù)合并為矩陣U:

    m × n 維矩陣U集合了流場的時(shí)空信息。按照POD 定義,將矩陣U分解為如下形式:

    式中:?j(x)為一個(gè)基函數(shù)(空間模態(tài)),aj(t)為對(duì)應(yīng)的時(shí)間系數(shù)。

    3)為保證基函數(shù)在最小二乘意義上是最優(yōu)的,求解矩陣U的協(xié)方差矩陣R的特征向量和特征值。協(xié)方差矩陣R定義如下:

    在工程應(yīng)用中,通常降階區(qū)域的網(wǎng)格數(shù)量m 很大,矩陣R的維度(m × m)也會(huì)非常大,使得求解困難。本文采用Sirovich[40]提出的快照POD 算法,將R簡化如下:

    求矩陣R的特征向量ψj和特征值λj的方法如下:

    4)在確定了上述較小維度特征向量ψj后,通過以下方式擴(kuò)展到原始 POD 模態(tài):

    時(shí)間系數(shù)aj(t)為空間向量u(t)在基上的投影:

    當(dāng)?j為單位化的基,上式可寫成:

    本文對(duì)列車表面溫度和壓力標(biāo)量進(jìn)行采樣,采樣時(shí)間間隔為0.065tref,快照總數(shù)為1 000 個(gè),對(duì)頭車、中間車以及尾車熱壓載荷分別做POD 降階分析。

    4 結(jié)果與討論

    4.1 車體熱壓載荷時(shí)均分布特征

    圖8 比較了管內(nèi)雍塞/非雍塞狀態(tài)下列車表面壓力與溫度分布。根據(jù)列車截面變化,將列車表面分成3 個(gè)區(qū)域:頭車流線型區(qū)域,即收縮段(?21.51

    從圖8(a)和(b)可以看出,雍塞狀態(tài)(工況A 和B)和非雍塞狀態(tài)(工況C)下的列車上表面熱壓載荷分布特征對(duì)應(yīng)著拉瓦爾噴管的兩種模式[41]。結(jié)合圖9(a)和(b)可以看到,對(duì)于工況A 和B,氣流在收縮段加速,車體上表面壓力和溫度快速下降;在平直車身區(qū)域,除了風(fēng)擋截面變化引起的局部流動(dòng)分離,車體上表面熱壓載荷繼續(xù)下降,但下降幅度小于收縮段,這是由于邊界層發(fā)展造成的氣流流動(dòng)面積進(jìn)一步減小,本質(zhì)上頭車流線型區(qū)域和平直車身區(qū)域均屬于收縮段,前者是幾何截面變化,后者屬于物理流動(dòng)面積變化;在擴(kuò)張段,氣流繼續(xù)加速至超聲速,直到遇到激波,氣流減速,溫度和壓力上升。同時(shí)比較工況A 和B,除了激波產(chǎn)生位置和激波角度(圖9(d)和(e)),兩者分布規(guī)律基本相同。增大阻塞比,激波位置向車尾鼻尖移動(dòng),易發(fā)生激波脫離[26]。與雍塞狀態(tài)相比,非雍塞狀態(tài)(工況C)下車體上表面熱壓載荷分布最大差異體現(xiàn)在擴(kuò)張段,由于喉部速度未加速到馬赫數(shù)Ma=1.0(圖9(c)),氣流在擴(kuò)張段減速,壓力和溫度升高,車尾也無激波產(chǎn)生。其次在平直車身區(qū)域,非雍塞狀態(tài)氣流順利通過,車體上表面壓力和溫度下降速度遠(yuǎn)小于雍塞狀態(tài)。

    圖9 y=0 平面時(shí)均馬赫數(shù)分布Fig.9 Mean Mach number distribution projected on a plane at y=0

    如圖8(c)所示,由于狹小的車軌間隙和懸浮架腔體的截面突變,列車下表面熱壓載荷分布較為復(fù)雜,使得列車核心部件—懸浮器處于惡劣的環(huán)境,非常值得關(guān)注,而這部分在已公開發(fā)表的文獻(xiàn)中鮮有涉及[22]。從壓力分布來看,沿流向(車頭→頭尾)總體呈下降趨勢,在懸浮架腔體處壓力出現(xiàn)大幅度波動(dòng),特別是頭車第一個(gè)B1-1 處,工況A 在此處的壓力系數(shù)波動(dòng)幅值達(dá)到3.36,減小阻塞比能夠緩解波動(dòng)強(qiáng)度。氣流從頭車鼻尖分流,一部分沿著頭車流線型向管道上方發(fā)展,另外一部分突入車軌間隙處,當(dāng)遇到第一個(gè)懸浮架腔體,會(huì)形成局部滯止點(diǎn),這也是壓力大幅度波動(dòng)的原因。

    圖8 不同流動(dòng)狀態(tài)下列車表面熱壓載荷時(shí)均分布對(duì)比Fig.8 Comparison of the time-averaged distribution of train surface temperature and pressure load under different flow conditions

    從溫度分布(圖8(d))來看,不論管內(nèi)是雍塞還是非雍塞狀態(tài),車底溫度都顯著高于車頂,這對(duì)車底蒙皮材料強(qiáng)度提出了挑戰(zhàn)。對(duì)于雍塞工況A 和B,由于氣流在平直車身區(qū)域阻塞,熱量無法及時(shí)耗散,在尾車第二個(gè)懸浮架B3-2 附近達(dá)到峰值,工況A 和B 的峰值分別為382 和375 K。在較小阻塞比(工況C)下,溫度峰值反而比雍塞工況大,高達(dá)420 K,位于尾車第一個(gè)懸浮架B3-1 腔內(nèi)。

    為探究其產(chǎn)生原因,繪制了尾車底部附近時(shí)均流線分布,如圖10所示。在工況A 中,激波使得尾車邊界層分離,隨后重新附著,在列車頂部和側(cè)墻下洗氣流與底部上洗氣流相互作用下再度分離,在尾車鼻尖附近形成了兩個(gè)連貫的分離泡V1 和V2,兩種氣流在鞍點(diǎn)S 處混合為一股氣流,向尾流后方發(fā)展。在工況B 中,列車頂部和側(cè)墻下洗氣流增強(qiáng),形成回流突入尾車流線型底部,并與底部上洗氣流相互作用,形成一條分割線l。在工況C 中,頂部和側(cè)墻下洗氣流繼續(xù)增強(qiáng),在尾流中占據(jù)主導(dǎo)地位,以回流方式進(jìn)入尾車底部,將底部上洗氣流“堵”在尾車第一個(gè)懸浮架B3-1 腔內(nèi),造成局部阻塞,產(chǎn)生局部高溫。

    圖10 y=0 平面時(shí)均流線與溫度分布Fig.10 Mean flow velocity streamlines projected on a plane at y=0 colored by temperature

    4.2 車體熱壓載荷瞬態(tài)波動(dòng)特征

    為了研究車體表面壓力和溫度載荷隨時(shí)間的演變規(guī)律,在列車上下表面典型位置布置一系列氣動(dòng)參數(shù)監(jiān)控點(diǎn),如圖11所示。其中序號(hào)P1-1和P3-1分別位于頭、尾車鼻尖,P1-2、P2-1和P3-2分別在頭、中、尾車平直車身頂部,P1-3、P1-4、P2-2、P2-3、P3-3和P3-4分別在懸浮架腔體內(nèi)部中心處。

    圖11 列車上下表面監(jiān)控點(diǎn)分布Fig.11 Distribution of monitoring points on the top and bottom surfaces of the train

    圖12 和13 分別給出了尾車鼻尖監(jiān)控點(diǎn)P3-1、中間車懸浮架腔內(nèi)監(jiān)控點(diǎn)P2-3的熱壓載荷時(shí)間歷程曲線和頻域分布。限于篇幅,其余監(jiān)控點(diǎn)的相關(guān)統(tǒng)計(jì)信息整理在表3 中,不再單獨(dú)贅述。圖中時(shí)間用無量綱時(shí)間t*表示:

    如圖12所示,隨著阻塞比增大,尾車鼻尖流場波動(dòng)加劇,特別是溫度載荷。工況A 的溫度波動(dòng)最大幅值(最高溫度–最低溫度)達(dá)到75 K。流動(dòng)分離和尾渦脫落是引起載荷波動(dòng)的主要原因,結(jié)合圖10可以發(fā)現(xiàn),工況A 尾車鼻尖附近存在兩個(gè)連貫的分離泡,工況B 僅有一個(gè)分離泡,而工況C 僅可觀察到一個(gè)較小的回流區(qū),因此相較于工況A 和B,工況C 載荷波動(dòng)較小。從頻域來看,對(duì)于阻塞比為0.3 的工況A,壓力和溫度載荷波動(dòng)的主頻分別為0.89 和13.35 Hz,同時(shí)還存在多個(gè)峰值的次頻,說明存在著多尺度非定常流場結(jié)構(gòu)。對(duì)于阻塞比分別為0.2 和0.1 的工況B 和C,壓力和溫度載荷波動(dòng)的主頻均為小于5 Hz 的低頻。

    如圖13所示,由于列車底部流動(dòng)非常復(fù)雜,監(jiān)控點(diǎn)熱壓載荷均值與阻塞比沒有對(duì)應(yīng)關(guān)系,主要由底部流動(dòng)狀態(tài)決定。3 個(gè)工況下壓力載荷均呈現(xiàn)準(zhǔn)周期波動(dòng),波動(dòng)頻率在15 Hz 附近。溫度載荷從總體看呈現(xiàn)上升趨勢,這是因?yàn)榈撞繜崃恳驘o法及時(shí)耗散而累積,使得溫度波動(dòng)主頻為低頻(0.89 Hz),同時(shí)溫度載荷在上升過程中,由于底部流動(dòng)分離和尾渦脫落也存在著小范圍的準(zhǔn)周期波動(dòng)。

    圖13 中間車懸浮架腔內(nèi)監(jiān)控點(diǎn)P2-3 瞬時(shí)載荷波動(dòng)和頻域分布Fig.13 Transient load fluctuations and frequency domain distribution of monitoring point P2-3 at the bogie cavity of the middle car

    4.3 POD 降階分析

    基于POD 降階方法,提取3 個(gè)工況頭車、中間車以及尾車的熱壓載荷主要振蕩模態(tài)和頻率,識(shí)別非定常強(qiáng)度較強(qiáng)區(qū)域。圖14 給出了尾車載荷前25 階模態(tài)能量占比。由圖可見,能量占比隨著模態(tài)階數(shù)增加衰減很快,前幾個(gè)高階模態(tài)對(duì)載荷非定常貢獻(xiàn)較大,捕獲了大部分振蕩特性,特別是溫度載荷,工況A、B 和C 的第一階模態(tài)能量占比分別為74%、68%和68%。因此本文重點(diǎn)分析前兩階模態(tài)模式與頻率特性。

    圖14 尾車熱壓載荷各階模態(tài)能量占比對(duì)比Fig.14 Comparison of the energy contribution of each mode in the pressure and temperature loads of the tail car

    圖15~17 分別給出了頭車、中間車和尾車的壓力載荷前兩階POD 模態(tài),其中為降階之后的無量綱壓力系數(shù)。如圖15所示,對(duì)頭車壓力載荷波動(dòng)貢獻(xiàn)較大的區(qū)域位于車底懸浮架腔內(nèi),從模態(tài)云圖來看,頭車懸浮架底部存在多種尺度的非定常流動(dòng)結(jié)構(gòu)。從時(shí)間系數(shù)的頻率分布亦可看出,除了主頻峰值外,還存在多個(gè)峰值的次頻。

    圖15 頭車壓力載荷前兩階POD 模態(tài)Fig.15 The first two POD modes of head car pressure load

    如圖16所示,與頭車壓力載荷波動(dòng)相比,中間車波動(dòng)較為規(guī)律。3 個(gè)工況下占主導(dǎo)地位的均是頻率在14 Hz 附近的單一流動(dòng)結(jié)構(gòu),再結(jié)合表3 和圖13(c)可以發(fā)現(xiàn),這與監(jiān)控點(diǎn)P2-2和P2-3的壓力主頻或次頻較為一致。從云圖分布來看,除了相位差導(dǎo)致的正負(fù)差異,3 個(gè)工況的模態(tài)結(jié)構(gòu)非常相似。

    圖16 中間車壓力載荷前兩階POD 模態(tài)Fig.16 The first two POD modes of middle car pressure load

    如圖17所示,尾車壓力波動(dòng)模態(tài)主要位于兩個(gè)區(qū)域。其一是列車上表面和側(cè)面,只存在雍塞工況(工況A 和B)。對(duì)于雍塞工況,尾車表面產(chǎn)生斜激波面,激波與邊界層相互作用,造成流動(dòng)分離,是一個(gè)非定常源。由于激波強(qiáng)度的差異,工況A 激波導(dǎo)致的非定常波動(dòng)出現(xiàn)在一階模態(tài),而工況B 則出現(xiàn)在二階模態(tài)。其二是列車底部懸浮架腔內(nèi),從模態(tài)云圖來看,波動(dòng)較強(qiáng)區(qū)域主要位于尾車第二個(gè)懸浮架B3-2 處,并且流動(dòng)結(jié)構(gòu)與中間車底部較為類似,主頻也在14 Hz 附近,這說明14 Hz 為該懸浮架的一個(gè)特征頻率。同時(shí)對(duì)于工況C,二階模態(tài)的振動(dòng)頻率除了主頻16.02 Hz 外,還存在29.37 Hz 次頻,該頻率下的流動(dòng)結(jié)構(gòu)可能與二階模態(tài)云圖中區(qū)域Ⅰ相關(guān)。這也與圖10 分析結(jié)果吻合,在工況C 下,頂部和側(cè)墻下洗氣流將底部上洗氣流“堵”在尾車第一個(gè)懸浮架B3-1 腔內(nèi),造成非定常脈動(dòng)。

    圖17 尾車壓力載荷前兩階POD 模態(tài)Fig.17 The first two POD modes of tail car pressure load

    圖18~20 分別給出了頭車、中間車和尾車的溫度載荷前兩階POD 模態(tài),其中T*為降階之后的無量綱溫度。如圖18所示,與壓力載荷類似,對(duì)頭車溫度載荷波動(dòng)貢獻(xiàn)較大區(qū)域也位于車底懸浮架腔內(nèi),特別是兩個(gè)懸浮架后端杜瓦表面。一階模態(tài)和二階模態(tài)空間結(jié)構(gòu)相似,且波動(dòng)頻率相同,因此一、二階模態(tài)為一對(duì)模態(tài),且存在相位差,這種關(guān)系對(duì)應(yīng)于流動(dòng)結(jié)構(gòu)中的對(duì)流[30,38]。

    圖18 頭車溫度載荷前兩階POD 模態(tài)Fig.18 The first two POD modes of head car temperature load

    如圖19所示,中間車溫度載荷一階模態(tài)與壓力載荷截然不同,從時(shí)間系數(shù)來看,溫度載荷一階模態(tài)呈現(xiàn)增長趨勢,這是因?yàn)橹虚g車底部熱量在密閉低氣壓環(huán)境下無法及時(shí)傳遞與耗散,隨著列車運(yùn)行累積在底部,這也是表3 中監(jiān)控點(diǎn)P2-2和P2-3溫度載荷波動(dòng)主頻為0.89 Hz 的原因。本質(zhì)上該現(xiàn)象不是非定常流動(dòng)結(jié)構(gòu)造成的,而是熱量累積導(dǎo)致的超低頻頻率。由于采樣時(shí)長等因素的影響,本文最低識(shí)別頻率為0.89 Hz。從模態(tài)云圖(圖19(a))看,熱量累積較強(qiáng)區(qū)域位于中間車第二個(gè)懸浮架B2-2前端區(qū)域Ⅱ,需要重點(diǎn)關(guān)注此處的溫度極限。溫度載荷二階模態(tài)則為非定常流動(dòng)結(jié)構(gòu)導(dǎo)致的周期波動(dòng),波動(dòng)主頻在14 Hz 附近,與壓力載荷波動(dòng)頻率一致。

    圖19 中間車溫度載荷前兩階POD 模態(tài)Fig.19 The first two POD modes of middle car temperature load

    如圖20所示,對(duì)于尾車溫度波動(dòng)一階模態(tài),工況A 與工況B、C 有較大差異,從時(shí)間系數(shù)來看,工況B、C 呈現(xiàn)上升趨勢,與中間車溫度載荷一階模態(tài)類似,而工況A 則無明顯增長趨勢。這主要是因?yàn)? 種工況中尾車底部也均存在著熱量累積,與中間車類似,熱量累積較強(qiáng)區(qū)域位于尾車懸浮架前端區(qū)域Ⅲ和Ⅳ處。但是工況A 尾車上表面還存在激波引發(fā)的局部非定常區(qū)域,與底部熱量累積效應(yīng)“中和”,因此在時(shí)間系數(shù)上表現(xiàn)為無明顯增長趨勢。這點(diǎn)在尾車鼻尖監(jiān)控點(diǎn)P3-1的瞬時(shí)溫度載荷波動(dòng)(圖12(b))中也有所體現(xiàn)。尾溫度載荷二階模態(tài)為非定常流動(dòng)結(jié)構(gòu)導(dǎo)致的周期波動(dòng),3 個(gè)工況波動(dòng)主頻均在14 Hz附近,與壓力載荷波動(dòng)頻率一致。同時(shí)觀察到工況C 的波動(dòng)幅值最大,這是由于下洗氣流回流導(dǎo)致的。

    圖20 尾車溫度載荷前兩階POD 模態(tài)Fig.20 The first two POD modes of tail car temperature load

    5 結(jié)論

    本文從時(shí)均分布、瞬時(shí)波動(dòng)以及POD 降階3 方面研究了真空管道磁浮交通雍塞狀態(tài)(β=0.3,0.2)與非雍塞狀態(tài)(β=0.1)下車體熱壓載荷分布特征及其非定常特性,得到結(jié)論如下:

    1)列車上表面載荷分布特征與拉瓦爾噴管相似。兩種狀態(tài)下載荷分布差異主要位于擴(kuò)張段,雍塞狀態(tài)下氣流加速至超聲速,載荷降低,直至激波產(chǎn)生,氣流減速。非雍塞狀態(tài)則相反,且不會(huì)產(chǎn)生激波。

    2)列車下表面熱壓載荷分布較為復(fù)雜。氣流在第一個(gè)懸浮架腔體形成局部滯止點(diǎn),造成壓力大幅度振蕩。同時(shí)熱量在列車底部無法及時(shí)耗散,雍塞狀態(tài)下在尾車第二個(gè)懸浮架達(dá)到峰值,而非雍塞工況的峰值則在尾車第一個(gè)懸浮架腔內(nèi),且大于雍塞工況的峰值。

    3)頭車壓力載荷脈動(dòng)較強(qiáng)區(qū)域位于車底懸浮架腔內(nèi),存在多種尺度的非定常流動(dòng)結(jié)構(gòu),造成中間車壓力波動(dòng)為14 Hz 附近的單一流動(dòng)結(jié)構(gòu)。雍塞狀態(tài)下尾車壓力波動(dòng)位于激波處和底部懸浮架腔內(nèi),非雍塞狀態(tài)則只位于懸浮架腔內(nèi)。

    4)頭車溫度載荷脈動(dòng)較強(qiáng)區(qū)域位于懸浮架后端杜瓦表面。中間車溫度載荷一階模態(tài)體現(xiàn)為熱量累積過程,二階模態(tài)為非定常流動(dòng)結(jié)構(gòu)導(dǎo)致的周期脈動(dòng)。不同阻塞比的尾車溫度載荷一階模態(tài)存在差異,工況B 和C 為熱量累積過程,而工況A 則是激波引發(fā)的局部非定常效應(yīng)與熱量累積耦合作用。

    在面向真空管道交通運(yùn)行的車輛設(shè)計(jì)過程中,需特別關(guān)注列車底部熱量累積,避免蒙皮材料損傷,同時(shí)需進(jìn)一步優(yōu)化平順懸浮架結(jié)構(gòu),降低載荷波動(dòng)。此外,雍塞狀態(tài)的尾車激波作用處也需特別關(guān)注。

    致謝:感謝國家超級(jí)計(jì)算鄭州中心支持!

    猜你喜歡
    尾車熱壓激波
    混勻堆料機(jī)尾車調(diào)偏研究
    一起熱壓三通開裂失效機(jī)理研究
    單堆單取斗輪機(jī)尾車防脫鉤保護(hù)裝置及連鎖
    一種基于聚類分析的二維激波模式識(shí)別算法
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    斗輪堆取料機(jī)尾車改造研究
    基于FPGA熱壓罐的溫度和氣壓控制研究
    陶瓷纖維擺塊式熱壓卷圓模設(shè)計(jì)
    模具制造(2019年3期)2019-06-06 02:10:52
    斜激波入射V形鈍前緣溢流口激波干擾研究
    SA266 Gr.2鍛件熱壓后的熱處理工藝
    欧美日韩中文字幕国产精品一区二区三区 | 啦啦啦免费观看视频1| 1024香蕉在线观看| 久久精品成人免费网站| 欧美精品av麻豆av| 97人妻天天添夜夜摸| 丝袜喷水一区| 国产av精品麻豆| 99热国产这里只有精品6| 老汉色∧v一级毛片| 精品高清国产在线一区| 少妇粗大呻吟视频| 国产一区二区三区综合在线观看| 别揉我奶头~嗯~啊~动态视频 | 国产一区二区三区综合在线观看| 亚洲熟女精品中文字幕| 免费在线观看日本一区| 久久久久久人人人人人| 亚洲精品国产精品久久久不卡| 国产成人影院久久av| 亚洲avbb在线观看| 在线十欧美十亚洲十日本专区| 亚洲少妇的诱惑av| 美女脱内裤让男人舔精品视频| 人人妻人人澡人人看| 亚洲精品国产av蜜桃| 国产精品二区激情视频| 后天国语完整版免费观看| 国产片内射在线| 亚洲一区中文字幕在线| 久热这里只有精品99| 在线看a的网站| 日韩大片免费观看网站| 夜夜骑夜夜射夜夜干| 人人妻,人人澡人人爽秒播| 亚洲精品粉嫩美女一区| 一区二区av电影网| 国产99久久九九免费精品| 亚洲精品中文字幕在线视频| 国产在线免费精品| 女警被强在线播放| 亚洲国产欧美网| 久久久久久久国产电影| 国产男女超爽视频在线观看| e午夜精品久久久久久久| 热re99久久国产66热| 老熟妇仑乱视频hdxx| 国产一区二区三区综合在线观看| 日韩一卡2卡3卡4卡2021年| 一区二区三区四区激情视频| 国产深夜福利视频在线观看| 久久人妻福利社区极品人妻图片| 成年女人毛片免费观看观看9 | 国产97色在线日韩免费| 午夜视频精品福利| 日韩,欧美,国产一区二区三区| 十八禁人妻一区二区| av一本久久久久| 丝袜脚勾引网站| 欧美在线黄色| 免费看十八禁软件| 中文字幕av电影在线播放| 国产高清视频在线播放一区 | 99久久99久久久精品蜜桃| 男人舔女人的私密视频| 丝袜美腿诱惑在线| 欧美另类亚洲清纯唯美| 欧美日韩av久久| 美女中出高潮动态图| 最近最新中文字幕大全免费视频| 91大片在线观看| 精品少妇一区二区三区视频日本电影| 99国产精品免费福利视频| tocl精华| 搡老熟女国产l中国老女人| 精品亚洲乱码少妇综合久久| 欧美精品一区二区大全| 国产一区二区 视频在线| 婷婷色av中文字幕| 91精品三级在线观看| 亚洲天堂av无毛| 脱女人内裤的视频| 精品亚洲乱码少妇综合久久| cao死你这个sao货| 国产高清videossex| 又紧又爽又黄一区二区| 久久中文看片网| 亚洲五月婷婷丁香| 91精品伊人久久大香线蕉| 成人国产一区最新在线观看| 国产伦人伦偷精品视频| 国产成人啪精品午夜网站| 久热爱精品视频在线9| 热99久久久久精品小说推荐| 91国产中文字幕| 久久青草综合色| 久久人人97超碰香蕉20202| 视频区欧美日本亚洲| 国产成人av教育| 色播在线永久视频| 精品一品国产午夜福利视频| 亚洲天堂av无毛| 中文字幕高清在线视频| 亚洲国产毛片av蜜桃av| 18在线观看网站| 亚洲欧美日韩高清在线视频 | 岛国毛片在线播放| 亚洲精品乱久久久久久| 日韩中文字幕欧美一区二区| 亚洲精品美女久久久久99蜜臀| 在线观看免费午夜福利视频| 国产免费一区二区三区四区乱码| 人妻 亚洲 视频| 视频在线观看一区二区三区| 窝窝影院91人妻| 日韩欧美国产一区二区入口| 天天操日日干夜夜撸| 国产人伦9x9x在线观看| 婷婷色av中文字幕| 成年美女黄网站色视频大全免费| 青春草亚洲视频在线观看| 妹子高潮喷水视频| 91成年电影在线观看| 桃花免费在线播放| 亚洲精品第二区| 亚洲男人天堂网一区| 国产熟女午夜一区二区三区| 免费在线观看视频国产中文字幕亚洲 | 亚洲熟女精品中文字幕| 不卡av一区二区三区| 考比视频在线观看| 免费在线观看黄色视频的| 一区二区三区乱码不卡18| 日韩 欧美 亚洲 中文字幕| 久久天躁狠狠躁夜夜2o2o| 国产日韩一区二区三区精品不卡| 成人av一区二区三区在线看 | 国产精品1区2区在线观看. | 国产成人欧美| 丝袜美腿诱惑在线| 9色porny在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 汤姆久久久久久久影院中文字幕| 男男h啪啪无遮挡| 精品亚洲成a人片在线观看| 国产成人精品久久二区二区91| 黑人猛操日本美女一级片| 国产亚洲一区二区精品| 国产av精品麻豆| 国产激情久久老熟女| 一级a爱视频在线免费观看| 高清视频免费观看一区二区| av又黄又爽大尺度在线免费看| 九色亚洲精品在线播放| 国产在线视频一区二区| 777米奇影视久久| 亚洲精品乱久久久久久| 另类精品久久| 成年美女黄网站色视频大全免费| 久久午夜综合久久蜜桃| 精品亚洲成国产av| 国产精品成人在线| 无限看片的www在线观看| 一区二区日韩欧美中文字幕| 一区二区av电影网| 曰老女人黄片| 亚洲国产日韩一区二区| 中文精品一卡2卡3卡4更新| 侵犯人妻中文字幕一二三四区| 十八禁网站网址无遮挡| 中文精品一卡2卡3卡4更新| 黑人巨大精品欧美一区二区mp4| 欧美激情极品国产一区二区三区| 熟女少妇亚洲综合色aaa.| 建设人人有责人人尽责人人享有的| 亚洲精品国产精品久久久不卡| 狂野欧美激情性xxxx| 侵犯人妻中文字幕一二三四区| 19禁男女啪啪无遮挡网站| 久久久久视频综合| 久久久久久人人人人人| 一级,二级,三级黄色视频| 日日摸夜夜添夜夜添小说| 国产一区有黄有色的免费视频| 亚洲人成电影观看| 色播在线永久视频| 国产欧美日韩一区二区三 | 亚洲欧美激情在线| 中文字幕另类日韩欧美亚洲嫩草| 青青草视频在线视频观看| 日韩欧美国产一区二区入口| 国产精品成人在线| 热99国产精品久久久久久7| 宅男免费午夜| 可以免费在线观看a视频的电影网站| 国产精品国产三级国产专区5o| 成人亚洲精品一区在线观看| 欧美久久黑人一区二区| 国产精品香港三级国产av潘金莲| 国产国语露脸激情在线看| 成年人午夜在线观看视频| 久久久久精品人妻al黑| 丰满人妻熟妇乱又伦精品不卡| 久久久久久久国产电影| 欧美亚洲 丝袜 人妻 在线| 午夜日韩欧美国产| 欧美黑人欧美精品刺激| 夫妻午夜视频| 国产亚洲精品第一综合不卡| 波多野结衣一区麻豆| 欧美性长视频在线观看| 女性生殖器流出的白浆| 国产亚洲欧美精品永久| 新久久久久国产一级毛片| 这个男人来自地球电影免费观看| 热re99久久国产66热| 日韩欧美一区二区三区在线观看 | 性高湖久久久久久久久免费观看| 国产日韩一区二区三区精品不卡| 青春草视频在线免费观看| 久久香蕉激情| 国产老妇伦熟女老妇高清| 精品少妇黑人巨大在线播放| 69av精品久久久久久 | 免费在线观看视频国产中文字幕亚洲 | 精品欧美一区二区三区在线| 国产精品久久久久久精品古装| 少妇裸体淫交视频免费看高清 | 男人添女人高潮全过程视频| 亚洲国产av新网站| 亚洲精品中文字幕在线视频| 老鸭窝网址在线观看| 精品卡一卡二卡四卡免费| 美女中出高潮动态图| 这个男人来自地球电影免费观看| 日本vs欧美在线观看视频| 久久久水蜜桃国产精品网| 丝袜美腿诱惑在线| 亚洲欧美成人综合另类久久久| 99国产精品99久久久久| 12—13女人毛片做爰片一| 欧美黄色片欧美黄色片| 国产一区二区在线观看av| 亚洲五月色婷婷综合| 精品久久久久久久毛片微露脸 | 9热在线视频观看99| 亚洲中文av在线| 国产高清videossex| 日本91视频免费播放| 久久九九热精品免费| 色老头精品视频在线观看| 久久久水蜜桃国产精品网| 多毛熟女@视频| 女人高潮潮喷娇喘18禁视频| 99热国产这里只有精品6| 久久久久国产精品人妻一区二区| 久久精品国产亚洲av高清一级| 久久国产亚洲av麻豆专区| 欧美午夜高清在线| 午夜久久久在线观看| 亚洲精品一区蜜桃| 中文字幕人妻熟女乱码| av免费在线观看网站| 久久人妻福利社区极品人妻图片| 99久久综合免费| 国产伦理片在线播放av一区| 午夜成年电影在线免费观看| 国产精品影院久久| 99香蕉大伊视频| 九色亚洲精品在线播放| 丝袜美腿诱惑在线| netflix在线观看网站| 亚洲,欧美精品.| 大码成人一级视频| 久久精品国产a三级三级三级| 久久九九热精品免费| 下体分泌物呈黄色| 满18在线观看网站| av不卡在线播放| 性高湖久久久久久久久免费观看| 脱女人内裤的视频| 在线精品无人区一区二区三| 国产日韩一区二区三区精品不卡| 国产欧美日韩一区二区三 | 亚洲av成人不卡在线观看播放网 | 午夜免费成人在线视频| 欧美日韩一级在线毛片| 丰满饥渴人妻一区二区三| 男女免费视频国产| 热99国产精品久久久久久7| 91成年电影在线观看| 午夜免费鲁丝| 亚洲精品一二三| 大码成人一级视频| 国产日韩欧美在线精品| 久久久久久亚洲精品国产蜜桃av| 99久久人妻综合| 亚洲色图 男人天堂 中文字幕| 一级a爱视频在线免费观看| 日本a在线网址| 热re99久久国产66热| 一二三四社区在线视频社区8| 欧美另类一区| 久久久精品免费免费高清| 欧美精品一区二区免费开放| 亚洲专区国产一区二区| 麻豆av在线久日| 两性午夜刺激爽爽歪歪视频在线观看 | 国产男人的电影天堂91| 巨乳人妻的诱惑在线观看| 国产精品秋霞免费鲁丝片| 日韩一区二区三区影片| 国产精品久久久久久人妻精品电影 | 国产精品熟女久久久久浪| 老司机影院毛片| 中文字幕av电影在线播放| 国产精品二区激情视频| 色视频在线一区二区三区| 欧美变态另类bdsm刘玥| 黑人猛操日本美女一级片| 欧美在线黄色| 成人三级做爰电影| 男女午夜视频在线观看| av片东京热男人的天堂| 丁香六月天网| 成人国产av品久久久| 精品一品国产午夜福利视频| 亚洲国产中文字幕在线视频| 男女无遮挡免费网站观看| 国产区一区二久久| 国产在线一区二区三区精| 在线观看免费午夜福利视频| 999久久久国产精品视频| 免费观看a级毛片全部| 一区二区三区激情视频| 欧美精品一区二区免费开放| 亚洲中文字幕日韩| 亚洲美女黄色视频免费看| 国产日韩欧美视频二区| 在线观看免费日韩欧美大片| 久久久精品区二区三区| 国产精品成人在线| 国产成人欧美| 国产极品粉嫩免费观看在线| 天天躁夜夜躁狠狠躁躁| 人人妻人人澡人人看| 亚洲国产成人一精品久久久| 国产欧美日韩综合在线一区二区| 999久久久国产精品视频| 久久久国产成人免费| 亚洲第一青青草原| www.精华液| 精品一品国产午夜福利视频| 欧美另类一区| 18禁国产床啪视频网站| 国产成人av激情在线播放| 久久亚洲国产成人精品v| 久久ye,这里只有精品| 国产精品久久久人人做人人爽| 熟女少妇亚洲综合色aaa.| 久久狼人影院| www.熟女人妻精品国产| 亚洲视频免费观看视频| 久久久久久久久免费视频了| 欧美日本中文国产一区发布| 91精品伊人久久大香线蕉| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品乱久久久久久| 美女午夜性视频免费| 可以免费在线观看a视频的电影网站| 免费黄频网站在线观看国产| 午夜激情久久久久久久| 国产免费av片在线观看野外av| √禁漫天堂资源中文www| 人妻 亚洲 视频| 欧美另类一区| 午夜福利,免费看| 青春草亚洲视频在线观看| 久久中文看片网| 国产一区二区 视频在线| 美女高潮喷水抽搐中文字幕| 日本av手机在线免费观看| 亚洲男人天堂网一区| 91老司机精品| 大香蕉久久网| 大片免费播放器 马上看| 久久女婷五月综合色啪小说| 无限看片的www在线观看| 免费观看av网站的网址| 下体分泌物呈黄色| 欧美亚洲 丝袜 人妻 在线| 激情视频va一区二区三区| 成年人黄色毛片网站| 免费在线观看影片大全网站| 午夜老司机福利片| 午夜福利视频在线观看免费| 亚洲精品国产一区二区精华液| 久久精品熟女亚洲av麻豆精品| av网站在线播放免费| 成年动漫av网址| 黑人猛操日本美女一级片| 欧美午夜高清在线| 久久99热这里只频精品6学生| 岛国毛片在线播放| 人妻人人澡人人爽人人| 交换朋友夫妻互换小说| 黄网站色视频无遮挡免费观看| 香蕉国产在线看| 黄色 视频免费看| 成人黄色视频免费在线看| 国产高清视频在线播放一区 | 午夜福利乱码中文字幕| 丰满迷人的少妇在线观看| 男女床上黄色一级片免费看| 免费一级毛片在线播放高清视频 | 亚洲免费av在线视频| 国产精品麻豆人妻色哟哟久久| 大片电影免费在线观看免费| 1024香蕉在线观看| 窝窝影院91人妻| 国产一区二区三区在线臀色熟女 | 亚洲一区二区三区欧美精品| 999精品在线视频| 一进一出抽搐动态| 欧美av亚洲av综合av国产av| 丰满人妻熟妇乱又伦精品不卡| 黄片小视频在线播放| 脱女人内裤的视频| 叶爱在线成人免费视频播放| 欧美日韩av久久| 嫁个100分男人电影在线观看| 国产成人欧美| 夜夜夜夜夜久久久久| 无遮挡黄片免费观看| 亚洲精品国产一区二区精华液| 国产不卡av网站在线观看| 日韩欧美免费精品| 免费在线观看视频国产中文字幕亚洲 | 国产免费视频播放在线视频| 叶爱在线成人免费视频播放| 国产精品99久久99久久久不卡| 国产黄频视频在线观看| 亚洲 国产 在线| 久久久久久亚洲精品国产蜜桃av| 两个人免费观看高清视频| 国产男女内射视频| 欧美黑人欧美精品刺激| 亚洲av成人一区二区三| 人成视频在线观看免费观看| 亚洲中文字幕日韩| 99久久国产精品久久久| 欧美97在线视频| 高清在线国产一区| 亚洲精品国产精品久久久不卡| 美国免费a级毛片| 国产成人精品无人区| 精品少妇一区二区三区视频日本电影| 欧美日韩福利视频一区二区| 久久久国产欧美日韩av| 久久久久精品国产欧美久久久 | 亚洲国产中文字幕在线视频| 精品少妇黑人巨大在线播放| av电影中文网址| 美女大奶头黄色视频| 在线观看www视频免费| 热99re8久久精品国产| 欧美av亚洲av综合av国产av| 制服人妻中文乱码| 人人妻人人澡人人爽人人夜夜| 老司机深夜福利视频在线观看 | 嫩草影视91久久| 中文字幕另类日韩欧美亚洲嫩草| 狂野欧美激情性bbbbbb| 69精品国产乱码久久久| videos熟女内射| 国产精品 欧美亚洲| 十八禁网站免费在线| 在线观看舔阴道视频| 99国产精品一区二区蜜桃av | 女性被躁到高潮视频| 国产成人影院久久av| 男人舔女人的私密视频| 一级片'在线观看视频| 欧美激情 高清一区二区三区| 窝窝影院91人妻| 免费女性裸体啪啪无遮挡网站| 亚洲五月色婷婷综合| 免费在线观看黄色视频的| 欧美精品高潮呻吟av久久| 国产高清国产精品国产三级| 一个人免费在线观看的高清视频 | 黄色视频,在线免费观看| 欧美日韩精品网址| 精品熟女少妇八av免费久了| 国产真人三级小视频在线观看| 久久热在线av| 亚洲性夜色夜夜综合| av免费在线观看网站| 亚洲一卡2卡3卡4卡5卡精品中文| 波多野结衣av一区二区av| 中文字幕最新亚洲高清| 国产精品一区二区在线观看99| 不卡av一区二区三区| 这个男人来自地球电影免费观看| 99久久人妻综合| 午夜成年电影在线免费观看| 五月开心婷婷网| 不卡av一区二区三区| 亚洲精品自拍成人| 国产深夜福利视频在线观看| 欧美国产精品一级二级三级| 国产亚洲欧美在线一区二区| 亚洲av片天天在线观看| 久久女婷五月综合色啪小说| 国产高清videossex| 9191精品国产免费久久| 蜜桃在线观看..| 亚洲欧洲精品一区二区精品久久久| 午夜福利乱码中文字幕| 日本一区二区免费在线视频| 欧美成人午夜精品| 亚洲人成电影观看| 韩国高清视频一区二区三区| 两性夫妻黄色片| 国产精品久久久久久精品古装| 国产av又大| 悠悠久久av| 最黄视频免费看| 久久人人97超碰香蕉20202| 久久这里只有精品19| 日本五十路高清| 女人精品久久久久毛片| av在线老鸭窝| 2018国产大陆天天弄谢| 精品一品国产午夜福利视频| 男女高潮啪啪啪动态图| 久久精品国产亚洲av香蕉五月 | 一区福利在线观看| 伦理电影免费视频| 亚洲男人天堂网一区| 黄色a级毛片大全视频| 日韩三级视频一区二区三区| 国产成人av激情在线播放| 久久国产精品大桥未久av| 少妇 在线观看| 色婷婷av一区二区三区视频| kizo精华| 国产成人精品久久二区二区免费| 亚洲精品自拍成人| 国产亚洲精品第一综合不卡| 久久人人爽av亚洲精品天堂| 免费女性裸体啪啪无遮挡网站| 超碰97精品在线观看| 亚洲国产欧美在线一区| 黄色视频不卡| av免费在线观看网站| 国产一区二区在线观看av| 1024香蕉在线观看| 亚洲欧美清纯卡通| 制服诱惑二区| 国产淫语在线视频| 亚洲成人手机| 黄片大片在线免费观看| 亚洲欧美色中文字幕在线| 久久99热这里只频精品6学生| 永久免费av网站大全| 91字幕亚洲| 伦理电影免费视频| 久久人妻福利社区极品人妻图片| 国产精品久久久久久精品古装| 午夜激情久久久久久久| videosex国产| 亚洲三区欧美一区| 爱豆传媒免费全集在线观看| 久久久水蜜桃国产精品网| 777久久人妻少妇嫩草av网站| 日韩制服骚丝袜av| 亚洲第一青青草原| 老司机亚洲免费影院| 美女主播在线视频| 国产精品久久久人人做人人爽| 精品久久久精品久久久| 亚洲色图 男人天堂 中文字幕| 国产淫语在线视频| 亚洲国产欧美网| 嫁个100分男人电影在线观看| 午夜影院在线不卡| 国产免费视频播放在线视频| 汤姆久久久久久久影院中文字幕| 午夜影院在线不卡| 少妇猛男粗大的猛烈进出视频| 18在线观看网站| 午夜影院在线不卡| 真人做人爱边吃奶动态| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久av美女十八| 丁香六月欧美| 亚洲人成电影观看| 国产真人三级小视频在线观看| 中文字幕高清在线视频| 青草久久国产| 在线观看免费视频网站a站| 亚洲av成人一区二区三| 成年人黄色毛片网站| 午夜福利乱码中文字幕| 一区二区三区乱码不卡18| 国产熟女午夜一区二区三区| 日韩精品免费视频一区二区三区| 99九九在线精品视频| 亚洲熟女精品中文字幕| 欧美变态另类bdsm刘玥| 老司机深夜福利视频在线观看 | 日本欧美视频一区| 啦啦啦视频在线资源免费观看| 汤姆久久久久久久影院中文字幕|