• <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鍛件熱壓后的熱處理工藝
    АⅤ资源中文在线天堂| 成人综合一区亚洲| 看十八女毛片水多多多| 26uuu在线亚洲综合色| 亚洲av电影在线观看一区二区三区 | 欧美区成人在线视频| 日本色播在线视频| 视频中文字幕在线观看| 女人被狂操c到高潮| 看黄色毛片网站| 欧美日韩一区二区视频在线观看视频在线 | 国产又黄又爽又无遮挡在线| 免费av不卡在线播放| 国产真实伦视频高清在线观看| 久久久久久久午夜电影| 日韩 亚洲 欧美在线| 国产综合懂色| 午夜福利在线观看免费完整高清在| 日本wwww免费看| 日本免费a在线| 国产麻豆成人av免费视频| 啦啦啦观看免费观看视频高清| 国产老妇女一区| 搡女人真爽免费视频火全软件| 超碰av人人做人人爽久久| 久久人人爽人人片av| 秋霞伦理黄片| 国产av在哪里看| 国产精品日韩av在线免费观看| 久久久色成人| 日韩一区二区视频免费看| 国产成人精品久久久久久| 亚洲国产精品专区欧美| 99久久精品国产国产毛片| 国产在视频线在精品| 国产精品1区2区在线观看.| 国产不卡一卡二| 日本av手机在线免费观看| 久久人人爽人人片av| 小蜜桃在线观看免费完整版高清| 国产av在哪里看| 噜噜噜噜噜久久久久久91| 国产麻豆成人av免费视频| 成人欧美大片| 日韩人妻高清精品专区| 欧美性猛交╳xxx乱大交人| 麻豆精品久久久久久蜜桃| 国产午夜精品一二区理论片| 人人妻人人澡欧美一区二区| 最近中文字幕高清免费大全6| 在线免费观看不下载黄p国产| 午夜爱爱视频在线播放| 国产高清视频在线观看网站| 一个人看的www免费观看视频| 成年女人看的毛片在线观看| 久久久久精品久久久久真实原创| 国产午夜福利久久久久久| 久久精品国产自在天天线| 一级二级三级毛片免费看| 亚洲人与动物交配视频| 免费观看精品视频网站| 精品不卡国产一区二区三区| 午夜a级毛片| 最近中文字幕高清免费大全6| 韩国av在线不卡| 啦啦啦韩国在线观看视频| 网址你懂的国产日韩在线| 亚洲av日韩在线播放| 小说图片视频综合网站| 久久精品91蜜桃| 婷婷色麻豆天堂久久 | 男人的好看免费观看在线视频| 男人和女人高潮做爰伦理| 久久99热6这里只有精品| 十八禁国产超污无遮挡网站| 亚洲成av人片在线播放无| av福利片在线观看| 亚洲精品自拍成人| 亚洲欧美成人综合另类久久久 | 国产激情偷乱视频一区二区| av在线亚洲专区| av天堂中文字幕网| 热99在线观看视频| 精品一区二区三区视频在线| 亚洲成av人片在线播放无| 亚洲成人精品中文字幕电影| 亚洲熟妇中文字幕五十中出| 五月伊人婷婷丁香| 看免费成人av毛片| 伦精品一区二区三区| 日本黄色视频三级网站网址| 亚洲国产高清在线一区二区三| 精品酒店卫生间| 两个人的视频大全免费| 日本与韩国留学比较| 久久这里只有精品中国| or卡值多少钱| 国产亚洲av嫩草精品影院| 啦啦啦啦在线视频资源| 国产精品麻豆人妻色哟哟久久 | 日韩视频在线欧美| 亚洲国产精品合色在线| 中文字幕av成人在线电影| 国产精品国产三级国产av玫瑰| 久久久久九九精品影院| 中文字幕av在线有码专区| 国产极品天堂在线| 男女视频在线观看网站免费| 可以在线观看毛片的网站| 日韩av不卡免费在线播放| 亚洲精品国产成人久久av| 韩国av在线不卡| 日韩国内少妇激情av| 久久精品久久久久久噜噜老黄 | 91久久精品电影网| 综合色丁香网| 超碰97精品在线观看| 国产色婷婷99| 91aial.com中文字幕在线观看| 免费观看人在逋| 午夜福利成人在线免费观看| 日本黄大片高清| 国产国拍精品亚洲av在线观看| 国产黄色视频一区二区在线观看 | 国模一区二区三区四区视频| 日韩一区二区视频免费看| 免费看a级黄色片| 两个人视频免费观看高清| 亚洲不卡免费看| 中文亚洲av片在线观看爽| 高清av免费在线| 国产中年淑女户外野战色| 天天一区二区日本电影三级| 亚洲熟妇中文字幕五十中出| 又爽又黄无遮挡网站| 精品酒店卫生间| 最近最新中文字幕免费大全7| 一本一本综合久久| 亚洲,欧美,日韩| 欧美bdsm另类| 亚洲国产色片| 久久鲁丝午夜福利片| 亚洲精品日韩av片在线观看| 日本一二三区视频观看| 亚洲,欧美,日韩| 草草在线视频免费看| 99久久精品国产国产毛片| 国产在线男女| 久久久成人免费电影| 99久国产av精品国产电影| 午夜视频国产福利| 精品久久久久久电影网 | 91aial.com中文字幕在线观看| 卡戴珊不雅视频在线播放| 老师上课跳d突然被开到最大视频| 亚洲精品aⅴ在线观看| 国产精品精品国产色婷婷| 午夜福利成人在线免费观看| 最近中文字幕2019免费版| 亚洲真实伦在线观看| 一个人免费在线观看电影| 久99久视频精品免费| 国产亚洲av片在线观看秒播厂 | 免费观看a级毛片全部| 色播亚洲综合网| 一级爰片在线观看| 久久国产乱子免费精品| 亚洲精品国产av成人精品| 成人一区二区视频在线观看| 色噜噜av男人的天堂激情| 一二三四中文在线观看免费高清| 日韩一区二区视频免费看| 亚洲精品日韩av片在线观看| 国产片特级美女逼逼视频| 18+在线观看网站| 成人三级黄色视频| 亚洲成人久久爱视频| 欧美激情国产日韩精品一区| 久久婷婷人人爽人人干人人爱| 美女cb高潮喷水在线观看| 国产精品不卡视频一区二区| 高清日韩中文字幕在线| 成人鲁丝片一二三区免费| 白带黄色成豆腐渣| 亚洲国产精品专区欧美| 日本一二三区视频观看| 欧美激情久久久久久爽电影| 变态另类丝袜制服| 亚洲美女搞黄在线观看| 汤姆久久久久久久影院中文字幕 | 性色avwww在线观看| 看片在线看免费视频| 在线观看一区二区三区| 国产精品国产三级国产专区5o | 欧美xxxx性猛交bbbb| 久久鲁丝午夜福利片| 国产男人的电影天堂91| 亚洲国产成人一精品久久久| 日韩av不卡免费在线播放| 91精品伊人久久大香线蕉| 欧美潮喷喷水| 亚洲av.av天堂| 人妻系列 视频| .国产精品久久| 91久久精品国产一区二区三区| 欧美极品一区二区三区四区| 免费电影在线观看免费观看| 嘟嘟电影网在线观看| 欧美三级亚洲精品| 男女视频在线观看网站免费| 国产成人一区二区在线| 亚洲熟妇中文字幕五十中出| 午夜福利成人在线免费观看| 九九在线视频观看精品| 一个人免费在线观看电影| 在线播放国产精品三级| 国产三级在线视频| 久久久国产成人免费| 久久鲁丝午夜福利片| 变态另类丝袜制服| 国产 一区 欧美 日韩| ponron亚洲| 欧美激情在线99| av专区在线播放| 久久婷婷人人爽人人干人人爱| 秋霞在线观看毛片| 精品一区二区三区人妻视频| 成人午夜高清在线视频| 欧美日本视频| 天堂√8在线中文| 午夜精品国产一区二区电影 | 成人一区二区视频在线观看| 我的女老师完整版在线观看| 国产一级毛片在线| 亚洲激情五月婷婷啪啪| 久久午夜福利片| 精品久久久久久久久久久久久| 亚洲精品日韩在线中文字幕| 在现免费观看毛片| 白带黄色成豆腐渣| 听说在线观看完整版免费高清| 男女那种视频在线观看| www.色视频.com| 国产视频内射| 看十八女毛片水多多多| 超碰av人人做人人爽久久| 成人美女网站在线观看视频| 最新中文字幕久久久久| 久久草成人影院| 成人毛片a级毛片在线播放| 亚洲欧美日韩卡通动漫| 亚洲在久久综合| 精品久久国产蜜桃| 国产色爽女视频免费观看| 国产精品无大码| 好男人在线观看高清免费视频| 人妻制服诱惑在线中文字幕| 深爱激情五月婷婷| 亚洲精华国产精华液的使用体验| 五月玫瑰六月丁香| 国产在视频线在精品| 日韩欧美精品v在线| 老司机影院成人| 波野结衣二区三区在线| 搡老妇女老女人老熟妇| 国产乱人偷精品视频| 亚洲第一区二区三区不卡| 精品久久久久久电影网 | eeuss影院久久| 成人av在线播放网站| 免费无遮挡裸体视频| 欧美性猛交黑人性爽| 尤物成人国产欧美一区二区三区| 深夜a级毛片| 水蜜桃什么品种好| 三级国产精品欧美在线观看| 性插视频无遮挡在线免费观看| 日韩在线高清观看一区二区三区| 九色成人免费人妻av| 在线播放国产精品三级| 国产精品av视频在线免费观看| 人妻系列 视频| 全区人妻精品视频| 久久韩国三级中文字幕| 天天一区二区日本电影三级| 亚洲av一区综合| 天堂av国产一区二区熟女人妻| 99久久精品热视频| 青春草国产在线视频| 我的老师免费观看完整版| 精品久久久久久久末码| 大香蕉久久网| 深爱激情五月婷婷| 91狼人影院| 有码 亚洲区| 精品久久久久久久久久久久久| 日本午夜av视频| 久久久久性生活片| 精品人妻偷拍中文字幕| 能在线免费观看的黄片| 久久久久久九九精品二区国产| 亚洲综合精品二区| 在线观看66精品国产| 日本爱情动作片www.在线观看| 亚洲av中文字字幕乱码综合| 久久久久久国产a免费观看| 观看免费一级毛片| 国产精品,欧美在线| 欧美成人午夜免费资源| 欧美bdsm另类| 国产精品久久电影中文字幕| 精品国产三级普通话版| 99热这里只有精品一区| 26uuu在线亚洲综合色| 免费看av在线观看网站| 国产极品天堂在线| 99热这里只有是精品在线观看| 中文天堂在线官网| 乱系列少妇在线播放| 亚洲欧美精品自产自拍| 国产亚洲av嫩草精品影院| 国产黄片美女视频| 久久精品夜夜夜夜夜久久蜜豆| 久热久热在线精品观看| 久久国产乱子免费精品| 晚上一个人看的免费电影| 国产三级在线视频| 婷婷色av中文字幕| 国产一区二区三区av在线| 在线免费观看不下载黄p国产| 国产高潮美女av| 最近手机中文字幕大全| 色网站视频免费| 91久久精品国产一区二区成人| 久久久久久久久中文| 国产不卡一卡二| 国产片特级美女逼逼视频| 国产精品伦人一区二区| 99久久精品国产国产毛片| 亚洲欧美成人综合另类久久久 | 日本一本二区三区精品| 国产高清不卡午夜福利| 免费不卡的大黄色大毛片视频在线观看 | 国产av码专区亚洲av| 日韩三级伦理在线观看| 国产中年淑女户外野战色| 国产毛片a区久久久久| 我的女老师完整版在线观看| 在线免费观看不下载黄p国产| 亚洲在久久综合| 日韩大片免费观看网站 | 免费一级毛片在线播放高清视频| av免费观看日本| 2021少妇久久久久久久久久久| 成人欧美大片| 五月伊人婷婷丁香| 日韩亚洲欧美综合| 亚洲在久久综合| 亚洲av福利一区| 亚洲欧美一区二区三区国产| 91精品一卡2卡3卡4卡| 亚洲国产高清在线一区二区三| 最近视频中文字幕2019在线8| 免费看日本二区| 成人亚洲欧美一区二区av| 九色成人免费人妻av| 男人狂女人下面高潮的视频| 欧美三级亚洲精品| 一夜夜www| 成人美女网站在线观看视频| 国产真实伦视频高清在线观看| 日韩欧美精品免费久久| 欧美高清成人免费视频www| www日本黄色视频网| 午夜精品在线福利| 日韩国内少妇激情av| av又黄又爽大尺度在线免费看 | 国产成人一区二区在线| 国产成人免费观看mmmm| 国产麻豆成人av免费视频| h日本视频在线播放| 中文字幕熟女人妻在线| av在线亚洲专区| 亚洲av一区综合| 国产成人freesex在线| 69人妻影院| 黄色一级大片看看| 日韩欧美精品免费久久| 三级国产精品欧美在线观看| av女优亚洲男人天堂| 欧美日韩一区二区视频在线观看视频在线 | 深爱激情五月婷婷| 久久精品国产亚洲av天美| 日本免费在线观看一区| 97在线视频观看| 高清在线视频一区二区三区 | 免费不卡的大黄色大毛片视频在线观看 | 日韩成人av中文字幕在线观看| 黄色欧美视频在线观看| www日本黄色视频网| 久99久视频精品免费| 欧美日韩综合久久久久久| 欧美性感艳星| 久久久久久久国产电影| 日韩av在线免费看完整版不卡| 99久久无色码亚洲精品果冻| 日韩在线高清观看一区二区三区| 国产精品国产三级国产av玫瑰| 国产精品伦人一区二区| 在线天堂最新版资源| 日本猛色少妇xxxxx猛交久久| 久久人妻av系列| 亚洲av福利一区| 久久这里有精品视频免费| 狂野欧美激情性xxxx在线观看| 国产午夜精品论理片| 少妇裸体淫交视频免费看高清| 国内精品一区二区在线观看| 一级毛片我不卡| 欧美三级亚洲精品| 51国产日韩欧美| 人人妻人人澡欧美一区二区| 欧美极品一区二区三区四区| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲成人中文字幕在线播放| 老司机福利观看| 搡女人真爽免费视频火全软件| 国产大屁股一区二区在线视频| 国产亚洲91精品色在线| 日韩一区二区视频免费看| 国产精品精品国产色婷婷| 国产精品乱码一区二三区的特点| 卡戴珊不雅视频在线播放| 人体艺术视频欧美日本| 在现免费观看毛片| 亚洲精品乱码久久久久久按摩| 国产精品国产三级国产专区5o | 亚洲在线自拍视频| 九色成人免费人妻av| 国产亚洲精品av在线| av在线播放精品| h日本视频在线播放| 日本一二三区视频观看| 中文天堂在线官网| 建设人人有责人人尽责人人享有的 | 啦啦啦韩国在线观看视频| 99在线人妻在线中文字幕| 久久久欧美国产精品| 在线观看美女被高潮喷水网站| 亚洲国产精品合色在线| 日本一二三区视频观看| 婷婷色麻豆天堂久久 | 国产一区亚洲一区在线观看| 国产精品人妻久久久久久| 亚洲av一区综合| 欧美+日韩+精品| av卡一久久| 亚洲电影在线观看av| 最近2019中文字幕mv第一页| 亚洲精品日韩av片在线观看| 国产综合懂色| 男女边吃奶边做爰视频| 免费看光身美女| 一区二区三区乱码不卡18| 我的老师免费观看完整版| 三级国产精品片| 26uuu在线亚洲综合色| 真实男女啪啪啪动态图| 91久久精品国产一区二区三区| 天天一区二区日本电影三级| 美女cb高潮喷水在线观看| 国产三级中文精品| 亚洲av不卡在线观看| 亚洲精品乱久久久久久| 在线播放无遮挡| 一个人观看的视频www高清免费观看| 亚洲国产高清在线一区二区三| 国产成人aa在线观看| 亚洲成人精品中文字幕电影| 中文字幕人妻熟人妻熟丝袜美| 男女边吃奶边做爰视频| www.av在线官网国产| 插阴视频在线观看视频| 91午夜精品亚洲一区二区三区| 一级av片app| 亚洲伊人久久精品综合 | 久久久久久久久久黄片| 国产又色又爽无遮挡免| 91精品伊人久久大香线蕉| 日本黄大片高清| 麻豆av噜噜一区二区三区| 亚洲伊人久久精品综合 | 黑人高潮一二区| 国产精品电影一区二区三区| 两个人视频免费观看高清| 欧美三级亚洲精品| 好男人视频免费观看在线| 欧美精品一区二区大全| 麻豆乱淫一区二区| 免费av毛片视频| 国产伦精品一区二区三区四那| 搡老妇女老女人老熟妇| 亚洲欧美精品自产自拍| 国产69精品久久久久777片| 国产成人91sexporn| 国产在线男女| 亚洲高清免费不卡视频| 亚洲激情五月婷婷啪啪| 日本三级黄在线观看| 色哟哟·www| 欧美+日韩+精品| 丝袜美腿在线中文| 小说图片视频综合网站| 国产亚洲av嫩草精品影院| 免费观看在线日韩| 精品欧美国产一区二区三| 最近中文字幕2019免费版| 久久精品久久精品一区二区三区| 日韩成人伦理影院| 久久久精品94久久精品| 国产欧美日韩精品一区二区| 亚洲内射少妇av| 波多野结衣巨乳人妻| 噜噜噜噜噜久久久久久91| 日日摸夜夜添夜夜添av毛片| 最近视频中文字幕2019在线8| 亚洲国产日韩欧美精品在线观看| 超碰av人人做人人爽久久| 搞女人的毛片| 99国产精品一区二区蜜桃av| 网址你懂的国产日韩在线| 日韩av在线大香蕉| h日本视频在线播放| 日韩高清综合在线| 在线天堂最新版资源| 国产老妇伦熟女老妇高清| 国产三级中文精品| 成人鲁丝片一二三区免费| 五月伊人婷婷丁香| 内地一区二区视频在线| 婷婷六月久久综合丁香| 久久人妻av系列| 一边亲一边摸免费视频| 亚洲av中文字字幕乱码综合| 中国国产av一级| 国产精品乱码一区二三区的特点| 69av精品久久久久久| 一级黄色大片毛片| 国产精品熟女久久久久浪| 日韩一本色道免费dvd| 麻豆精品久久久久久蜜桃| 国产一区二区在线观看日韩| 久久人妻av系列| 一个人免费在线观看电影| 亚洲av日韩在线播放| 干丝袜人妻中文字幕| 99热网站在线观看| 亚洲性久久影院| 午夜福利视频1000在线观看| 久久国产乱子免费精品| 久久精品熟女亚洲av麻豆精品 | 麻豆成人av视频| 午夜福利在线在线| 精品人妻一区二区三区麻豆| 成人漫画全彩无遮挡| 天堂中文最新版在线下载 | 搡老妇女老女人老熟妇| 午夜福利网站1000一区二区三区| 亚洲人成网站在线播| 淫秽高清视频在线观看| 全区人妻精品视频| 99久国产av精品| 久久久久久国产a免费观看| 在线观看66精品国产| 国产精品一二三区在线看| 免费av毛片视频| 啦啦啦观看免费观看视频高清| 国产精品日韩av在线免费观看| 久久精品国产亚洲av涩爱| 搡老妇女老女人老熟妇| 国产精品国产高清国产av| 午夜老司机福利剧场| 日本-黄色视频高清免费观看| 三级毛片av免费| 男女下面进入的视频免费午夜| 高清在线视频一区二区三区 | 国产精品熟女久久久久浪| 亚洲精品自拍成人| 国产v大片淫在线免费观看| 亚洲综合色惰| 久久这里有精品视频免费| 综合色丁香网| 精品不卡国产一区二区三区| 国产淫片久久久久久久久| 亚洲精品久久久久久婷婷小说 | 欧美激情久久久久久爽电影| 性色avwww在线观看| 日韩一本色道免费dvd| 欧美成人免费av一区二区三区| 七月丁香在线播放| 天堂av国产一区二区熟女人妻| 国产黄片视频在线免费观看| 熟女电影av网| 一级黄片播放器| 在线观看一区二区三区| 国产日韩欧美在线精品| 国产精品精品国产色婷婷| 一个人看的www免费观看视频| 国产精品嫩草影院av在线观看| 中文字幕av成人在线电影| 亚洲欧美精品综合久久99| 国产精品国产三级国产av玫瑰| 赤兔流量卡办理| 在线观看美女被高潮喷水网站| av在线亚洲专区|