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

    真空管道列車流固耦合研究進展及關鍵技術分析

    2023-07-12 01:12:46寇杰符澄高興龍孫運強
    實驗流體力學 2023年3期

    寇杰,符澄,高興龍,孫運強

    中國空氣動力研究與發(fā)展中心 設備設計與測試技術研究所,綿陽 621000

    0 引 言

    作為一種新型軌道交通工具,磁浮列車不存在傳統(tǒng)列車的輪軌黏著限制、輪軌摩擦阻力等問題,理論上可以實現(xiàn)更高速度的商業(yè)運行,近年來受到越來越多研究者的關注,已經(jīng)成為軌道交通領域的熱點研究對象[1-2]。

    磁浮列車依靠懸浮力運行,與運行軌道間有一定間隙,高速運行時存在明顯的地面效應;同時,高速磁浮列車運行的目標速度為600~1 000 km/h(圖1 為我國正在研發(fā)的600 km/h 速度等級高速磁浮列車),遠高于現(xiàn)有高速列車。高速運行條件下的氣動阻力等氣動載荷更大,對磁浮列車的外形、結(jié)構及控制系統(tǒng)設計都帶來了巨大挑戰(zhàn)[3]。為實現(xiàn)磁浮列車高速運行,研究者提出了真空管道列車概念,其中以2013 年埃隆·馬斯克提出的Hyperloop 膠囊列車最為著名[4-5]。真空管道列車的原理是通過抽吸空氣實現(xiàn)封閉管道內(nèi)部的低壓狀態(tài),降低磁浮列車周圍流場空氣密度,從而大幅降低列車運行中的氣動阻力,實現(xiàn)超高速運行[6]。

    圖1 600 km/h 速度等級磁浮列車[1]Fig.1 The maglev train at speed of 600 km/h [1]

    低壓環(huán)境及管道的封閉邊界條件,給真空管道列車帶來了一系列氣動問題。列車在管道內(nèi)超高速運動,其前方形成壓縮波,后方產(chǎn)生膨脹波[7];壓縮波和膨脹波在管道內(nèi)傳播,進一步形成各種斜激波;各種激波反射疊加,引起管道內(nèi)壓力波動;再疊加車頭、車身和車尾通過時的壓力擾動,整個管道內(nèi)的流場環(huán)境十分復雜。真空管道列車前方和后方產(chǎn)生的壓縮波包括斜激波、正激波、反射激波、Lamda 激波、菱形激波等復雜的激波簇結(jié)構,這些復雜結(jié)構與列車相互作用,形成高度非穩(wěn)定的列車運行環(huán)境[8-9]。懸浮運行的真空管道列車和管道之間存在明顯間隙,列車運行中的不穩(wěn)定將導致流場環(huán)境邊界條件發(fā)生變化,進一步影響真空管道內(nèi)部流場。因此,真空管道列車周圍的氣體壓力和速度均受到列車速度和姿態(tài)的影響,存在嚴重的流固耦合效應。

    本文對真空管道列車氣動研究進展及軌道列車流固耦合特性研究進展等進行了綜述,分析了真空管道列車流固耦合關鍵技術,并對下一步研究提出了建議。

    1 真空管道列車氣動研究進展

    1.1 真空管道列車氣動特性

    采用真空管道的目的是降低列車運行中的氣動阻力,真空管道列車的氣動阻力特性是研究者重點關注的問題。周曉等[10]建立了二維不可壓縮真空管道列車數(shù)值計算模型,研究了真空管道氣體壓力、阻塞比、列車速度對氣動阻力的影響,發(fā)現(xiàn)阻塞比增大會導致列車氣動阻力急劇增大,管道內(nèi)部壓力降低則會降低氣動阻力,表明真空環(huán)境對于列車高速運行具有基礎性作用。Kim 等[11]采用相似方法開展研究,得出的真空管道列車氣動阻力與列車速度、阻塞比和管道壓力的關系與周曉等的研究具有較好的一致性;Kim 等還關注了真空管道列車前方產(chǎn)生激波的臨界速度,該速度主要與阻塞比相關。劉加利等[12-13]將真空管道列車數(shù)值計算模型拓展至三維,建立了3 節(jié)車列車模型和真空管道內(nèi)部流場的三維計算模型,如圖2 所示。

    圖2 真空管道列車三維計算模型[12]Fig.2 A three-dimensional calculation model for a vacuum pipe train[12]

    該研究將列車氣動阻力D 分為兩部分,即壓差阻力Dp和摩擦阻力Dτ。為分析各阻力的關系,定義了無量綱的氣動阻力系數(shù)CD、壓差阻力系數(shù)CDp和摩擦阻力系數(shù)CDτ:

    研究發(fā)現(xiàn),在不同運行狀態(tài)下,摩擦阻力系數(shù)CDτ均遠小于壓差阻力系數(shù)CDp,因此,壓差阻力對列車氣動阻力起決定性作用。從圖3 和4 可以看出,當管道壓力和列車速度固定時,氣動阻力隨阻塞比增大而增大。低壓條件下,列車速度和管道壓力對列車壓差阻力和氣動阻力基本無影響,而阻塞比則對氣動阻力影響明顯。

    圖3 不同管道壓力下氣動阻力隨阻塞比的變化[12]Fig.3 Aerodynamic drag changes with the blockage under different pipeline pressures[12]

    圖4 不同列車速度下氣動阻力隨阻塞比的變化[12]Fig.4 Aerodynamic drag changes with the blockage under different speeds[12]

    在真空管道列車氣動載荷特性研究方面,王博[14]采用三維數(shù)值仿真與模型試驗方法對真空管道高溫超導磁浮列車氣動阻力與管道壓力、列車速度、阻塞比的關系進行了研究,提出了一種低壓慣性試驗方法。黃尊地等[15]依據(jù)最小空間尺寸計算克努森數(shù),判斷真空管道內(nèi)的流動狀態(tài),考慮三維定常可壓縮效應,利用滑移網(wǎng)格技術,分析列車速度、管道真空度、阻塞比及環(huán)境溫度對列車氣動阻力的影響。王志飛等[16]提出了一種基于數(shù)值仿真和正交理論相結(jié)合的參數(shù)設計研究方法,分析了影響真空管道列車運行的3 種因素(列車速度、阻塞比和管道壓力)對列車氣動阻力的影響趨勢及各因素影響的顯著性。馮瑞龍等[17]基于RBF(Radial Basis Function)和BP 神經(jīng)網(wǎng)絡(Back Propagation Network),訓練了一種真空管道列車氣動阻力預測模型,分析了阻塞比、列車速度和管道壓力影響下的氣動阻力。

    陳緒勇[18-19]和Zhang[20]等則關注了不同頭車外形對真空管道列車氣動阻力的影響,建立了二維計算模型,對不同頭車外形的真空管道列車氣動阻力進行數(shù)值分析。Pandey 等[21]分析了頭車外形與尾車外形對真空管道列車氣動阻力的影響,如圖5 和6 所示,不同的頭車和尾車流線型設計對車頭和車尾的氣壓差影響有限,因而對列車壓差阻力影響較小。Yang 等[22]也研究了頭車外形和尾車外形對真空管道列車氣動阻力的影響,并在研究中加入了管道壓力、列車速度等影響因素。Zhang[23]研究了真空管道直徑和阻塞比對列車氣動阻力的影響,研究結(jié)果顯示對于半徑1.5 m 的真空管道列車,最佳管道直徑為2~4 m。Bi 等[24]構建了一段550 m 的封閉管道,研究發(fā)現(xiàn)真空管道列車的氣動阻力在激波和反射波共同作用下是一個時變量,且受到激波傳播的時間影響,如圖7 所示(圖中α 為阻塞度、β 為管道內(nèi)氣壓與標準大氣壓比值、vt為列車運行速度)。Ma 等[25]建立了真空管道列車氣動阻力與阻塞比、列車速度和管道壓力之間的理論模型,并采用試驗裝置(圖8)獲得的數(shù)據(jù)進行了驗證,但該研究僅針對真空管道列車低速運行開展,且理論分析結(jié)果與試驗結(jié)果差異明顯。

    圖5 不同頭車形狀條件下的氣壓差[19]Fig.5 Pressure differentials under different header shapes[19]

    圖6 不同尾車形狀條件下的氣壓差[19]Fig.6 Pressure differentials under different trailer shapes[19]

    圖7 封閉管道內(nèi)氣動阻力時域變化曲線[24]Fig.7 A change curve of the aerodynamic drag within the closed pipe[24]

    圖8 真空管道列車試驗裝置[25]Fig.8 Vacuum pipe train test device[25]

    從現(xiàn)有研究可以看出,真空管道列車氣動阻力主要與列車速度、阻塞比、管道壓力和列車外形等因素有關。管道內(nèi)氣體的阻塞效應對列車氣動阻力系數(shù)起決定性作用,需在真空管道列車氣動設計中重點考慮管道的阻塞比[26-28]。真空管道列車實際運行中,在管道內(nèi)還可能受到升力、橫向氣動力和氣動力矩的影響,但現(xiàn)有氣動載荷特性研究基本局限于氣動阻力研究,尚未對其他氣動載荷給予充分關注。

    1.2 真空管道列車流場特性

    列車在具有一定真空度的封閉管道內(nèi)運行,帶動管道內(nèi)氣體流動,在列車周圍形成一定的氣壓波。管道內(nèi)存在壓縮波、膨脹波及各類斜激波組成的復雜波系[29],嚴重影響真空管道列車運行的穩(wěn)定性和安全性,同時還是管道壁面承載壓力的主要影響因素。因此,對真空管道內(nèi)各種壓力波的產(chǎn)生、傳播及演變問題應給予關注。

    真空管道列車壓力波的產(chǎn)生機制和傳播機理復雜。張曉涵等[30]建立了二維亞聲速真空管道列車模型,采用數(shù)值計算方法研究了真空管道列車運行前方氣動壅塞現(xiàn)象與后方非對稱尾部激波現(xiàn)象。圖9為不同列車速度下真空管道氣動壅塞特性的對比(圖中,Ma0、Mat分別為來流馬赫數(shù)、管道喉道馬赫數(shù),m0和m1分別為列車-管道前方自由空間的空氣流量、自由來流進入列車-管道之間的空氣流量)。列車低速運行時,真空管道內(nèi)未發(fā)生壅塞現(xiàn)象;亞聲速運行時,管道內(nèi)發(fā)生臨界壅塞現(xiàn)象,列車頂部出現(xiàn)馬赫數(shù)為1 的臨界點(即“等熵極限”);在高亞聲速、超聲速運行時,管道內(nèi)持續(xù)保持雍塞,列車后方產(chǎn)生弓形激波,影響列車后方的氣體流動。

    圖9 真空管道內(nèi)部等熵流動流場特征分析[30]Fig.9 Analysis of flow field characteristics of isentropic flow in vacuum pipe[30]

    列車在高亞聲速和超聲速運行時,其頭部前方出現(xiàn)了由溢出氣流堆積而成的高壓區(qū)。通過計算,可以擬合出列車前方的雍塞長度和列車運行時間的關系:

    式中:L 為列車前方壅塞長度,t 為列車運行時間,Ma 為列車運行馬赫數(shù)。

    圖10 為600 km/h 運行速度下真空管道列車尾部區(qū)域的馬赫數(shù)分布。此時真空管道列車尾部區(qū)域存在膨脹波和激波。當來流為高亞聲速時,尾部形成局部超聲速區(qū)域,誘發(fā)產(chǎn)生一系列膨脹波,膨脹波經(jīng)過尾部流線型、管道上表面的多次反射與疊加,形成了較短的激波。激波在管道上表面邊界層中與剪切層相互作用,導致列車上方和管道內(nèi)壁之間激波反射不穩(wěn)定,產(chǎn)生激波串并向后耗散。

    圖10 600 km/h 運行速度下列車尾部馬赫數(shù)分布[30]Fig.10 Mach number distribution at train trail at speed of 600 km/h[30]

    Oh[31]、Gillani[32]、Kang[33]、Zhou[34]和Bao[35]等采用穩(wěn)態(tài)分析方法研究了壓力波與列車運行速度、列車長度、管道壓力、管道溫度和阻塞率的關系。列車實際運行產(chǎn)生的壓力波并非完全穩(wěn)定,其前方壓縮波、后方膨脹波及激波的傳播與反射疊加,以及管道的壁面效應,使壓力波呈現(xiàn)非穩(wěn)態(tài)[36]。Kim 等[11]將真空管道列車后方斜激波的產(chǎn)生和傳播機理作為非穩(wěn)態(tài)問題進行研究,Niu[37]、Sui[38]和Yang[22]等也采用非穩(wěn)態(tài)二維數(shù)值模擬方法研究了斜激波與列車速度、管道壓力、阻塞率等的關系。Niu 等[39]還利用非穩(wěn)態(tài)數(shù)值模擬方法研究了列車加速和減速過程中斜激波的特性。Sui[38]和Bi[24]等則研究了壓力波和膨脹波的產(chǎn)生和傳播機理。周鵬等[8]研究了真空管道列車超高速運行中產(chǎn)生的各類波系(圖11 和12 展示了弓形激波、正激波、反射激波、Lamda 激波、菱形激波等激波簇結(jié)構的產(chǎn)生和傳播特性。圖中,P1、P2、P3 分別為激波的反射位置1~3;R2、R3 為反射波2 和3;I1 為弓形激波1;N 為正激波;M3 為馬赫桿3),研究結(jié)果表明,真空管道內(nèi)的各類激波使列車氣動環(huán)境比運行于開放流場的普通高速列車更為復雜,且各類激波的傳播、反射相互疊加,使列車運行流場環(huán)境具有高度的非穩(wěn)定性。

    圖11 車頭處壓力云圖和激波結(jié)構圖[8]Fig.11 Pressure clouds and shock cluster structure maps at the head of the vehicle[8]

    圖12 車尾處壓力云圖和激波結(jié)構圖[8]Fig.12 Pressure clouds and shock cluster structure maps at the end of the vehicle[8]

    不穩(wěn)定的激波簇結(jié)構對真空管道列車產(chǎn)生時變氣動載荷,在時變氣動升力、氣動力矩及其他各種帶有沖擊性的氣動載荷作用下,列車運行特性必將發(fā)生變化。因此,對真空管道列車氣動阻力之外的氣動載荷也亟待開展深入研究。

    2 軌道列車流固耦合研究進展

    作為流體力學和固體力學相互耦合發(fā)展出的力學分支,流固耦合研究在機翼顫振、油箱內(nèi)燃油晃動和高層建筑風載變形等實際工程問題中均體現(xiàn)出重要價值。傳統(tǒng)流固耦合研究主要研究流場和固體場相互作用的結(jié)果,在數(shù)值計算方法方面,需要研究兩種場的交界面上的各種物理參數(shù)的平衡與協(xié)調(diào)。軌道列車的流固耦合研究屬于流體力學和多體動力學耦合的交叉學科研究,與傳統(tǒng)流固耦合研究重點關注的流體力學和固體力學耦合差別較大,但這種交叉學科研究在軌道交通研究領域仍被稱為“流固耦合研究”。

    真空管道列車與傳統(tǒng)軌道列車具有一定的相似性,其流固耦合研究可以借鑒軌道列車流固耦合研究的分析方法、分析手段并參考相關分析結(jié)果。

    2.1 軌道列車流固耦合分析方法

    軌道列車流固耦合研究重點關注外界環(huán)境風作用下的列車運行狀態(tài)。于夢閣[40-41]、Bettle[42]、王永冠[43]、Baker[44]、Thomas[45]、李田[46]等利用氣動計算軟件或理論模型分析得到了外界環(huán)境風作用下高速列車受到的氣動載荷,進而采用車輛動力學計算軟件分析了氣動載荷作用下高速列車的運動特性。該離線方法可以快速計算高速列車氣動載荷及動力學指標,但并未考慮氣動載荷導致的列車運行姿態(tài)變化及變化的列車運行姿態(tài)對流場特性產(chǎn)生的影響。

    楊吉忠等[47]采用有限體積法和任意拉格朗日歐拉法(Arbitrary Lagrangian Eulerian,ALE)進行了氣動計算和車輛動力學計算的耦合分析。與之類似,為分析側(cè)風作用下的列車動力學特性,崔濤基于ALE 及車輛動力學模型對比分析了流固耦合效應對車輛氣動力和氣動力矩的影響,發(fā)現(xiàn)在列車氣動分析中考慮列車運行姿態(tài)變化和列車流固耦合效應,會對車體搖頭力矩(側(cè)向力不作用于質(zhì)心時產(chǎn)生的力矩)產(chǎn)生明顯影響。崔濤等[48]建立了列車流體動力學模型,并采用參數(shù)傳遞、求解控制和動態(tài)網(wǎng)格技術將氣動載荷參數(shù)從流體動力學模型傳遞至車輛動力學模型,還提出針對不同問題可以采用不同耦合方法聯(lián)立流場方程和車輛動力學方程。聯(lián)立耦合方法可分為離線耦合、顯式耦合和隱式耦合:離線耦合僅在車輛動力學計算中引入氣動載荷;顯式耦合利用車輛系統(tǒng)響應子模塊與氣動子模塊相互交換數(shù)據(jù)實現(xiàn)實時交互求解;與前兩者相比,隱式耦合選擇性地進行數(shù)據(jù)交互,需在車輛系統(tǒng)響應子模塊和氣動子模塊都達到穩(wěn)定狀態(tài)的前提下進行數(shù)據(jù)交互。

    李田等[49-52]構建了基于氣動軟件Fluent 和車輛動力學軟件Simpack 的聯(lián)合仿真模型。如圖13 所示,利用一個接口程序?qū)luent 和Simpack 相互連接,以Fluent 計算得到氣動載荷,經(jīng)接口程序傳遞至Simpack,以Simpack 計算列車在氣動載荷下的運行姿態(tài)和運動參數(shù),并將相關參數(shù)經(jīng)接口程序傳遞至Fluent,F(xiàn)luent 根據(jù)更新的輸入?yún)?shù)進行新一步計算。

    圖13 聯(lián)合仿真求解過程Fig.13 Joint simulation resolution process

    在車輛-軌道耦合動力學方程中引入2 個積分參數(shù)(μ和λ),用以構造新的顯式積分格式:

    式中:Xn-1、Xn和Xn+1分別表示第n-1、n 和n +1 個迭代步的位移矩陣;Δt 為時間積分步長。

    基于Fluent 和Simpack 的流固耦合計算需經(jīng)接口程序?qū)崿F(xiàn)兩者之間的數(shù)據(jù)交換,其中一個程序的計算過程中,另一程序基本處于等待參數(shù)的狀態(tài);同時,氣動計算和車輛動力學計算的時間步長存在數(shù)量級差異,導致后者等待時間較長,嚴重影響計算效率。此外,由于是給定時間步長進行迭代,氣動計算比車輛動力學計算的迭代步長更長,導致氣動計算得到的氣動載荷在2 個輸出的迭代步之間變化較大,可能導致車輛動力學計算在突變氣動載荷輸入下不收斂。

    針對聯(lián)合仿真的計算效率和計算收斂問題,李田[53]提出了一種內(nèi)嵌式高速列車流固耦合聯(lián)合仿真方法,如圖14 所示。

    圖14 內(nèi)嵌式聯(lián)合仿真求解過程Fig.14 Embedded joint simulation resolution process

    該方法改進了車輛動力學計算部分,通過編程將車輛動力學計算程序包內(nèi)嵌至氣動計算軟件,無需借助接口程序進行數(shù)據(jù)傳遞,在氣動計算軟件中即可實現(xiàn)氣動計算和車輛動力學計算,避免車輛動力學計算因缺乏輸入長期處于等待狀態(tài),大幅提高了聯(lián)合計算效率。同時,為降低突變氣動載荷對計算收斂性的影響,在車輛動力學計算程序中線性插值當前時刻和上一時刻的氣動載荷,使得加載的氣動載荷連續(xù)。t 時刻作用于列車-軌道耦合動力學模型的氣動載荷Fa表達式為:

    式中:tn、tn+1分別為第n 個和第n + 1 個迭代步時間;Fn、Fn+1分別為tn、tn+1時刻Fluent 計算得到的氣動載荷。

    在常值的橫風作用下,列車的氣動狀態(tài)和車輛動力學狀態(tài)會達到穩(wěn)定?;诖嗽?,李田等[54]提出了一種新的流固耦合計算方法—平衡狀態(tài)法,用以快速計算橫風作用下的高速列車流固耦合特性。在計算中,為判定計算的氣動力是否達到穩(wěn)定狀態(tài),對每一項指標設置一定的收斂判斷準則:

    式中:Fn和Fn+1分別為第n 個和第n + 1 個迭代步的氣動力指標值;ε 為收斂誤差,取ε=0.001。

    在平衡狀態(tài)法中,首先需計算初始狀態(tài)下高速列車的穩(wěn)態(tài)流場,并對計算得到的氣動載荷進行收斂性判定,若達到收斂狀態(tài)則將氣動載荷輸入車輛動力學計算程序,計算氣動載荷作用下的列車姿態(tài)參數(shù),并根據(jù)姿態(tài)參數(shù)更新氣動計算的網(wǎng)格數(shù)據(jù),作新一步氣動計算后再次輸出氣動載荷數(shù)據(jù),如是循環(huán),直至氣動計算結(jié)果和車輛動力學計算結(jié)果均達到收斂,則表明計算進入平衡狀態(tài),即可輸出車輛動力學和氣動計算結(jié)果。如圖15 所示。

    圖15 平衡狀態(tài)法聯(lián)合仿真求解過程[54]Fig.15 Joint simulation resolution process in balanced state method[54]

    根據(jù)李田的研究,可以得到如表1 所示的3 種流固耦合方法計算效率的對比。離線仿真法的氣動計算總迭代步數(shù)和動力學計算總時間均最少,但存在前文所述的準確性和收斂性問題,難以反映實際的流固耦合特性。對比平衡狀態(tài)法和交互式聯(lián)合仿真法可以發(fā)現(xiàn),前者氣動計算每個時間步的迭代步數(shù)(100 步)多于后者(20 步),但前者的時間迭代總步數(shù)(50 步)遠少于后者(7 500 步)。由于平衡狀態(tài)法每次迭代均需使用動力學計算,其動力學計算總時間比交互式聯(lián)合仿真法多250 s,但在整個流固耦合仿真過程中,氣動計算所需時間遠多于動力學計算,因此平衡狀態(tài)法計算效率更高。

    表1 計算效率比較[54]Table 1 The contrast of computational efficiency[54]

    在目前的軌道列車流固耦合計算方法研究中,李田等的研究成果兼顧了計算精度和計算效率,具有較好的創(chuàng)新性和實用性;提出的聯(lián)合仿真法、內(nèi)嵌式聯(lián)合仿真法、平衡狀態(tài)法均可作為真空管道列車流固耦合研究的參考和借鑒(在研究中,仍需考慮真空管道列車流固耦合研究的特殊性)。

    2.2 流固耦合作用下的列車運行特性

    在流固耦合作用下,列車的動力學性能受到氣動載荷的激勵發(fā)生明顯改變。為準確計算氣動載荷作用下的列車運行特性,研究者采用流固耦合計算方法將氣動計算和車輛動力學計算聯(lián)系起來,相關研究包括了第2.1 節(jié)提及的離線式流固耦合計算和交互式流固耦合計算。

    表2 和3 分別為采用不同流固耦合計算方法得到的列車氣動載荷指標和運行姿態(tài)指標[54]。從表中可以看出:平衡狀態(tài)法和交互式聯(lián)合仿真法得到的結(jié)果差異相對較?。浑x線仿真法得到的氣動載荷指標和運行姿態(tài)指標與其他2 種方法差異明顯。離線仿真法得到的結(jié)果不夠準確,難以體現(xiàn)真實的列車流固耦合特性。

    表2 氣動載荷指標比較[54]Table 2 The contrast of aerodynamic load[54]

    表3 列車運行姿態(tài)指標比較[54]Table 3 The contrast of vehicle body attitude[54]

    崔濤等[55-56]采用離線式流固耦合計算方法研究了高速列車通過站臺的安全性,但該方法僅適用于無環(huán)境風條件。崔濤[55]采用顯式耦合方法分析了高速列車交會時的流場特性和列車運行安全性,發(fā)現(xiàn)在考慮流固耦合效應時,交會時的高速列車表面壓力波動、氣動載荷、列車運行的動力學響應及運行安全性指標均有一定程度增大,列車交會的安全性降低。圖16 和17 給出了有/無環(huán)境風條件下高速列車350 km/h 等速交會時的表面壓力分布。崔濤等[57-59]還采用穩(wěn)態(tài)流固耦合計算方法對側(cè)風環(huán)境下高速列車交會安全性開展了研究。在側(cè)風作用下,列車姿態(tài)發(fā)生變化,但當側(cè)風恒定時,列車姿態(tài)最終進入穩(wěn)定狀態(tài)。因此,可以采用隱式耦合方法對側(cè)風作用下列車姿態(tài)變化的最終狀態(tài)進行分析,并采用離線仿真方法計算列車的動力學指標。如圖18 所示,與不考慮列車姿態(tài)變化相比,在考慮列車姿態(tài)變化的情況下,頭車受到的橫向氣動力、垂向氣動力、側(cè)滾氣動力矩和俯仰氣動力矩均明顯增大。由于氣動載荷增大,列車運行脫軌系數(shù)和輪重減載率等均明顯增大,列車運行的安全性降低,如圖19 所示。

    圖16 無風條件下350 km/h 等速交會[57]Fig.16 Same velocity rendezvous at speed of 350 km/h without wind[57]

    圖17 10 m/s 側(cè)風條件下350 km/h 等速交會 [57]Fig.17 Same velocity rendezvous at speed of 350 km/h with sidewind at speed of 10 m/s[57]

    圖18 頭車、中間車和尾車的氣動載荷[57]Fig.18 Aerodynamic load of head,middle and tail vehicle[57]

    圖19 運行安全性指標[57]Fig.19 The operational safety factor[57]

    現(xiàn)有的列車流固耦合計算分析結(jié)果表明,流固耦合效應會導致列車氣動載荷指標和安全性指標增大,列車運行安全性降低,在相關研究中必須考慮流固耦合效應。此外,離線仿真法雖然具有更高的效率,但是由于未實時傳輸數(shù)據(jù),導致計算結(jié)果偏離真實情況。因此,在真空管道列車流固耦合研究中,應采用在線的流固耦合計算方法。

    3 真空管道列車流固耦合關鍵技術

    3.1 真空管道列車流場分析技術

    真空管道列車高速運行過程中,復雜變化的管道壓力環(huán)境導致車外壓力波動及列車壓力載荷變化,激勵列車系統(tǒng)產(chǎn)生動力學響應,使列車運行姿態(tài)產(chǎn)生變化,進而加劇流場變化。與傳統(tǒng)高速列車流場分析相比,真空管道列車流場分析具有特殊性。真空管道形成的封閉環(huán)境使列車運行環(huán)境中的波系更為復雜,懸浮運行的列車姿態(tài)變化也需在流場分析中作為重要輸入。基于現(xiàn)有真空管道流場研究成果,在后續(xù)研究中需重點將真空管道列車姿態(tài)這一因素引入流場分析,建立計算精度較高,能夠準確反映流場特性,體現(xiàn)復雜波系的產(chǎn)生、反射和傳播特性且兼顧計算效率的流場分析模型。針對不同問題,可以采用二維模型或三維模型進行計算,提高研究效率。

    3.2 真空管道列車流固耦合聯(lián)合仿真技術

    真空管道列車流固耦合分析需同時進行氣動計算和車輛動力學計算。針對高速列車流固耦合問題,李田、崔濤等提出了多種流固耦合計算方法,如離線計算、交互式計算、內(nèi)嵌式計算及平衡狀態(tài)法等,主要思路是聯(lián)合氣動計算軟件與動力學計算軟件(或程序),相互傳遞氣動載荷數(shù)據(jù)和列車姿態(tài)參數(shù)作為彼此的輸入條件。

    真空管道內(nèi)波系復雜,列車處于時變載荷狀態(tài),在激波作用下,列車局部可能形成突變載荷,尤其是在列車加速或減速狀態(tài)下。因此,在流固耦合計算中需實時傳輸時變氣動載荷數(shù)據(jù)和列車姿態(tài)參數(shù)。此外,列車姿態(tài)參數(shù)對真空管道內(nèi)波系的產(chǎn)生和傳播條件也存在影響,列車姿態(tài)變化幅度過大,會產(chǎn)生新的壓縮波或膨脹波。因此,針對真空管道列車復雜多變的流場環(huán)境及懸浮運行的車輛動力學特點,需基于現(xiàn)有高速列車流固耦合計算方法和技術,提出一套可準確計算真空管道列車流固耦合問題的方法和技術。在列車動力學計算模型中應施加時變氣動力、氣動力矩及局部氣動力,以體現(xiàn)氣動載荷的特殊性。

    3.3 氣動載荷下的真空管道列車控制技術

    真空管道列車懸浮運行于真空管道中,在流固耦合效應及管道內(nèi)各種激波作用下,列車受到時變氣動載荷及沖擊載荷,其運行姿態(tài)發(fā)生改變。為避免列車姿態(tài)及氣動環(huán)境發(fā)生劇烈變化,導致列車運行安全性降低(甚至與管道發(fā)生碰撞),需在列車動力學計算模型中加入列車控制技術,將氣動載荷下的列車運行姿態(tài)變化控制于合適范圍內(nèi),確保列車安全運行。

    控制技術研究可以從2 個方面開展:一方面,設計合理的真空管道列車氣動外形,減小列車的氣動力變化;另一方面,可以在列車內(nèi)部設置懸掛及控制機構,主動調(diào)整列車相對懸浮軌道的姿態(tài)。

    4 結(jié)論與展望

    真空管道列車利用磁懸浮技術和管道真空運行環(huán)境,理論上可實現(xiàn)超過1 000 km/h 的運行速度。但封閉的真空管道導致氣動環(huán)境復雜,同時列車懸浮運行使列車運行姿態(tài)極易發(fā)生改變,流固耦合效應明顯。本文針對真空管道列車的流固耦合問題,對真空管道列車氣動研究進展、軌道列車流固耦合特性研究進展進行了綜述,分析了真空管道列車流固耦合關鍵技術,得到以下結(jié)論:

    1)真空管道列車管道內(nèi)氣體的阻塞效應對列車氣動阻力系數(shù)起決定性作用?,F(xiàn)有氣動載荷特性研究局限于對真空管道列車氣動阻力的研究,而列車在實際運行中還受到升力、橫向氣動力、氣動力矩的影響,研究者需對這些氣動載荷給予關注。

    2)真空管道內(nèi)存在各類激波,其氣動環(huán)境比開放流場遠為復雜,且各類激波的傳播、反射相互疊加使列車運行流場環(huán)境具有高度的非穩(wěn)定性。不穩(wěn)定的激波簇結(jié)構對列車產(chǎn)生時變氣動載荷,列車在時變的氣動升力、氣動力矩及各種沖擊性氣動載荷作用下,運行特性發(fā)生變化。

    3)在軌道列車流固耦合研究中提出的聯(lián)合仿真法、內(nèi)嵌式聯(lián)合仿真法及平衡狀態(tài)法,均可作為真空管道列車流固耦合研究的參考和借鑒,但真空管道列車流固耦合研究還需考慮真空管道列車的特殊性,重點關注時變氣動載荷和局部沖擊載荷的影響。在動力學計算模型中應施加時變氣動力、氣動力矩及局部氣動力,以體現(xiàn)氣動載荷的特殊性。

    4)流固耦合效應會導致列車受到的氣動載荷指標和列車安全性指標增大,運行安全性降低。離線流固耦合分析方法雖然具有更高的效率,但會導致計算結(jié)果偏離真實情況。在真空管道列車流固耦合研究中,需采用在線流固耦合計算方法。

    在真空管道列車流固耦合問題研究中,需重點研究開發(fā)真空管道列車流場分析技術,關注列車姿態(tài)對流場特性的影響;開發(fā)真空管道列車流固耦合分析技術,實現(xiàn)高效準確的流固雙向耦合計算,準確分析列車氣動特性和動力學特性;開發(fā)真空管道列車控制技術,實現(xiàn)真空管道列車受控安全運行。

    流固耦合問題是真空管道列車技術發(fā)展面臨的重要問題。在未來研究中,需進一步借鑒傳統(tǒng)高速列車和常規(guī)磁浮列車的研究方法和手段,并積極尋求創(chuàng)新與突破。

    亚洲成人中文字幕在线播放| 美女内射精品一级片tv| 22中文网久久字幕| 91精品一卡2卡3卡4卡| 久久久a久久爽久久v久久| 一卡2卡三卡四卡精品乱码亚洲| 干丝袜人妻中文字幕| 久久综合国产亚洲精品| 亚洲一级一片aⅴ在线观看| 91aial.com中文字幕在线观看| 午夜福利在线在线| 91狼人影院| 亚洲国产精品国产精品| 婷婷色av中文字幕| 亚洲欧美日韩高清专用| 午夜精品在线福利| 日韩av在线大香蕉| 狂野欧美激情性xxxx在线观看| 亚洲精品色激情综合| 日韩三级伦理在线观看| 99久久精品热视频| 成人三级黄色视频| www日本黄色视频网| 午夜福利在线在线| 亚洲精品日韩在线中文字幕| 久久久久国产网址| 国产精品乱码一区二三区的特点| 日韩国内少妇激情av| 熟妇人妻久久中文字幕3abv| 美女被艹到高潮喷水动态| 欧美成人午夜免费资源| 欧美成人午夜免费资源| 欧美成人午夜免费资源| 女的被弄到高潮叫床怎么办| 18禁动态无遮挡网站| 欧美丝袜亚洲另类| 又爽又黄无遮挡网站| 超碰97精品在线观看| 桃色一区二区三区在线观看| 一级黄片播放器| 久久精品国产99精品国产亚洲性色| 久久精品国产99精品国产亚洲性色| 亚洲av免费高清在线观看| 性色avwww在线观看| 午夜亚洲福利在线播放| 免费观看性生交大片5| av在线天堂中文字幕| 18禁在线播放成人免费| 国国产精品蜜臀av免费| 亚洲熟妇中文字幕五十中出| 久久久久久国产a免费观看| 久久精品综合一区二区三区| 最近最新中文字幕免费大全7| 国产伦精品一区二区三区视频9| 久久久久久久久久黄片| 99在线人妻在线中文字幕| 99在线人妻在线中文字幕| 日日撸夜夜添| 精品一区二区免费观看| 色吧在线观看| 国产精品国产三级国产av玫瑰| 午夜精品一区二区三区免费看| 免费看av在线观看网站| 国产片特级美女逼逼视频| 中文乱码字字幕精品一区二区三区 | 免费黄色在线免费观看| 大香蕉97超碰在线| 免费黄色在线免费观看| 欧美成人精品欧美一级黄| 欧美成人精品欧美一级黄| 中文在线观看免费www的网站| 国产亚洲av片在线观看秒播厂 | 卡戴珊不雅视频在线播放| 男插女下体视频免费在线播放| 亚洲欧美一区二区三区国产| 国产中年淑女户外野战色| 国产一级毛片在线| 久久草成人影院| 亚洲精品乱码久久久v下载方式| 亚洲三级黄色毛片| 久久久久久久久中文| 波多野结衣巨乳人妻| 高清视频免费观看一区二区 | 国产伦精品一区二区三区四那| 99久久精品一区二区三区| 人人妻人人看人人澡| 小蜜桃在线观看免费完整版高清| 成人综合一区亚洲| 久久这里有精品视频免费| 日韩欧美在线乱码| 免费av观看视频| 免费看光身美女| 国产午夜福利久久久久久| 寂寞人妻少妇视频99o| 亚洲精品久久久久久婷婷小说 | 中文字幕熟女人妻在线| 青青草视频在线视频观看| 91在线精品国自产拍蜜月| av卡一久久| 国产老妇伦熟女老妇高清| 特大巨黑吊av在线直播| 99热这里只有是精品50| 搞女人的毛片| 国产探花极品一区二区| 国产精品.久久久| 男女啪啪激烈高潮av片| 看非洲黑人一级黄片| 亚洲av成人精品一二三区| 国产男人的电影天堂91| 中国国产av一级| 亚洲国产精品成人综合色| 午夜福利在线观看吧| 99久久精品国产国产毛片| 99热全是精品| 高清毛片免费看| 秋霞伦理黄片| 99久久九九国产精品国产免费| 在线观看av片永久免费下载| 精品人妻视频免费看| 日本午夜av视频| 亚洲内射少妇av| 精品久久久久久久末码| 麻豆成人午夜福利视频| 人妻系列 视频| av黄色大香蕉| 国产又色又爽无遮挡免| 插逼视频在线观看| 日韩欧美精品v在线| 爱豆传媒免费全集在线观看| 一级毛片久久久久久久久女| 18禁在线无遮挡免费观看视频| 亚洲精品国产av成人精品| 一边摸一边抽搐一进一小说| 一本一本综合久久| 麻豆精品久久久久久蜜桃| 又黄又爽又刺激的免费视频.| 久久精品国产亚洲av涩爱| 国产av码专区亚洲av| 精品人妻偷拍中文字幕| 婷婷六月久久综合丁香| 久久久久久久久久黄片| 久久综合国产亚洲精品| 插阴视频在线观看视频| 非洲黑人性xxxx精品又粗又长| 精品99又大又爽又粗少妇毛片| 成人综合一区亚洲| 亚洲欧美日韩东京热| 人妻夜夜爽99麻豆av| 少妇人妻一区二区三区视频| 国产精品一及| 日韩成人av中文字幕在线观看| 禁无遮挡网站| 国产精品国产三级国产专区5o | 国产三级在线视频| 亚洲av.av天堂| 综合色丁香网| 中文乱码字字幕精品一区二区三区 | 免费电影在线观看免费观看| 美女xxoo啪啪120秒动态图| 亚洲国产最新在线播放| 色视频www国产| 亚洲图色成人| 国产亚洲精品av在线| 免费看光身美女| 免费无遮挡裸体视频| 亚洲精品,欧美精品| 亚洲av成人精品一二三区| 国产亚洲最大av| 欧美xxxx黑人xx丫x性爽| 水蜜桃什么品种好| 国产乱人偷精品视频| 亚洲av中文字字幕乱码综合| 三级男女做爰猛烈吃奶摸视频| 欧美日本亚洲视频在线播放| 亚州av有码| 欧美zozozo另类| 99热全是精品| 欧美一级a爱片免费观看看| 国产伦一二天堂av在线观看| 精品一区二区三区人妻视频| 美女大奶头视频| 亚洲av成人av| 日日摸夜夜添夜夜爱| 久久精品国产亚洲网站| 国产色爽女视频免费观看| 久久鲁丝午夜福利片| 成人毛片a级毛片在线播放| 成人一区二区视频在线观看| 在线观看av片永久免费下载| 只有这里有精品99| 国产伦精品一区二区三区四那| videos熟女内射| 午夜精品一区二区三区免费看| 免费观看a级毛片全部| 99在线人妻在线中文字幕| 色播亚洲综合网| 欧美丝袜亚洲另类| 久久精品国产鲁丝片午夜精品| 国产探花极品一区二区| 日本-黄色视频高清免费观看| 日韩制服骚丝袜av| 免费看美女性在线毛片视频| 成人性生交大片免费视频hd| 婷婷色av中文字幕| 99热精品在线国产| 国产极品精品免费视频能看的| 欧美日韩在线观看h| 噜噜噜噜噜久久久久久91| 午夜福利视频1000在线观看| 久久精品熟女亚洲av麻豆精品 | 99热6这里只有精品| 亚洲欧美中文字幕日韩二区| 色尼玛亚洲综合影院| 国产成人一区二区在线| 国产精品一区二区三区四区免费观看| 夫妻性生交免费视频一级片| 国产免费福利视频在线观看| 久久热精品热| 精品午夜福利在线看| 卡戴珊不雅视频在线播放| 午夜激情欧美在线| 亚洲在线观看片| 国产av码专区亚洲av| 久久精品国产自在天天线| av在线天堂中文字幕| 中文字幕精品亚洲无线码一区| 日韩欧美国产在线观看| 夫妻性生交免费视频一级片| av国产免费在线观看| 九色成人免费人妻av| 久久久色成人| 国产色爽女视频免费观看| 免费av不卡在线播放| 少妇人妻一区二区三区视频| or卡值多少钱| 午夜激情福利司机影院| 中文字幕久久专区| 99视频精品全部免费 在线| 3wmmmm亚洲av在线观看| 久久精品久久精品一区二区三区| 秋霞在线观看毛片| 黄色欧美视频在线观看| 精品不卡国产一区二区三区| 日本五十路高清| 97在线视频观看| 简卡轻食公司| 日本一本二区三区精品| 国产极品精品免费视频能看的| 亚洲伊人久久精品综合 | 一级av片app| 高清视频免费观看一区二区 | 亚洲自拍偷在线| 超碰av人人做人人爽久久| 非洲黑人性xxxx精品又粗又长| 久久久精品大字幕| 免费观看在线日韩| 我的女老师完整版在线观看| 国产精品久久视频播放| 日本三级黄在线观看| 免费无遮挡裸体视频| 搡女人真爽免费视频火全软件| av在线播放精品| 欧美3d第一页| 久久久久久久久久成人| 日日啪夜夜撸| 亚洲最大成人av| 51国产日韩欧美| 亚洲丝袜综合中文字幕| 午夜福利在线观看免费完整高清在| 亚洲18禁久久av| 在线免费十八禁| av黄色大香蕉| 九九在线视频观看精品| 我的女老师完整版在线观看| 一级毛片我不卡| 欧美另类亚洲清纯唯美| 日韩,欧美,国产一区二区三区 | 干丝袜人妻中文字幕| 看片在线看免费视频| 麻豆久久精品国产亚洲av| 国产国拍精品亚洲av在线观看| 1024手机看黄色片| 日韩视频在线欧美| 永久网站在线| 久久精品影院6| 最近的中文字幕免费完整| 久久久久久久久中文| 日韩精品青青久久久久久| 丰满乱子伦码专区| 非洲黑人性xxxx精品又粗又长| 免费av观看视频| av播播在线观看一区| 欧美激情在线99| 成年av动漫网址| 精品国产露脸久久av麻豆 | 亚洲精品aⅴ在线观看| 变态另类丝袜制服| 三级毛片av免费| 中文资源天堂在线| 精品人妻视频免费看| 男女视频在线观看网站免费| 99久久成人亚洲精品观看| 永久网站在线| 成年女人永久免费观看视频| 在线免费观看不下载黄p国产| 欧美高清成人免费视频www| 久久99蜜桃精品久久| 免费观看人在逋| 中文字幕免费在线视频6| 夫妻性生交免费视频一级片| 久久久久久久久久黄片| 1024手机看黄色片| 69av精品久久久久久| 欧美丝袜亚洲另类| 欧美+日韩+精品| 国产一区二区在线观看日韩| 久久6这里有精品| 精品99又大又爽又粗少妇毛片| 亚洲国产精品合色在线| 亚洲综合色惰| 一本一本综合久久| 国产欧美日韩精品一区二区| 神马国产精品三级电影在线观看| 九九热线精品视视频播放| 久久精品国产亚洲网站| 亚洲国产精品合色在线| 日韩精品有码人妻一区| 一区二区三区乱码不卡18| 男女那种视频在线观看| h日本视频在线播放| 午夜日本视频在线| 欧美一区二区国产精品久久精品| 国产69精品久久久久777片| 国内精品一区二区在线观看| 国产伦精品一区二区三区视频9| 亚洲国产最新在线播放| 午夜精品一区二区三区免费看| 男人舔奶头视频| 成人午夜高清在线视频| 国产综合懂色| 国产女主播在线喷水免费视频网站 | 亚洲欧美日韩东京热| 国内揄拍国产精品人妻在线| 久久午夜福利片| 国产成人a区在线观看| 久久热精品热| 两个人视频免费观看高清| 午夜久久久久精精品| 亚洲精品aⅴ在线观看| 听说在线观看完整版免费高清| 久久久国产成人精品二区| 黄色配什么色好看| 免费av毛片视频| 国产真实乱freesex| 国产午夜精品论理片| 搡女人真爽免费视频火全软件| 最近手机中文字幕大全| 国产激情偷乱视频一区二区| 中文欧美无线码| 18禁在线无遮挡免费观看视频| 亚洲三级黄色毛片| 国产亚洲5aaaaa淫片| 欧美激情在线99| 国产一区有黄有色的免费视频 | 免费av观看视频| av免费观看日本| 日韩高清综合在线| 免费观看精品视频网站| 深夜a级毛片| 久久国内精品自在自线图片| 久久久欧美国产精品| av免费观看日本| 一区二区三区乱码不卡18| 亚洲丝袜综合中文字幕| 亚洲成人精品中文字幕电影| 一个人观看的视频www高清免费观看| 成人美女网站在线观看视频| 2021天堂中文幕一二区在线观| 亚洲aⅴ乱码一区二区在线播放| 国产免费又黄又爽又色| 亚洲精品久久久久久婷婷小说 | 久久久国产成人精品二区| 国产黄片美女视频| 久久人人爽人人爽人人片va| 女人被狂操c到高潮| 一级黄色大片毛片| 成年女人永久免费观看视频| 午夜精品国产一区二区电影 | 美女内射精品一级片tv| 亚洲在线观看片| 日本午夜av视频| 久久久色成人| 精品少妇黑人巨大在线播放 | 国产乱人偷精品视频| 国产精品国产三级专区第一集| 欧美成人一区二区免费高清观看| 午夜免费男女啪啪视频观看| 色噜噜av男人的天堂激情| 99久国产av精品国产电影| 亚洲欧美日韩东京热| 国产精品一区二区三区四区久久| 国产单亲对白刺激| 日本黄色视频三级网站网址| 激情 狠狠 欧美| 国产日韩欧美在线精品| 国产av在哪里看| 天天躁日日操中文字幕| av专区在线播放| 亚洲美女视频黄频| 熟女人妻精品中文字幕| 国产黄片美女视频| 欧美性猛交黑人性爽| 午夜久久久久精精品| 亚洲精品成人久久久久久| 国产片特级美女逼逼视频| 久久久久国产网址| 亚洲精品国产成人久久av| 老女人水多毛片| 美女大奶头视频| 亚洲人与动物交配视频| 夫妻性生交免费视频一级片| 我要看日韩黄色一级片| 最近最新中文字幕大全电影3| 日本与韩国留学比较| 好男人在线观看高清免费视频| 黄色配什么色好看| 精品久久国产蜜桃| 亚洲天堂国产精品一区在线| 亚洲成人中文字幕在线播放| 午夜爱爱视频在线播放| 亚洲欧美日韩东京热| 亚洲18禁久久av| 欧美极品一区二区三区四区| 久久久久久九九精品二区国产| 色播亚洲综合网| 看黄色毛片网站| 色吧在线观看| 99久国产av精品| 麻豆久久精品国产亚洲av| 国产免费视频播放在线视频 | 波多野结衣巨乳人妻| 亚洲美女视频黄频| 久久精品国产亚洲网站| 日本-黄色视频高清免费观看| 一二三四中文在线观看免费高清| 日韩欧美在线乱码| 免费电影在线观看免费观看| 超碰97精品在线观看| 日韩av不卡免费在线播放| 久久午夜福利片| 黄色一级大片看看| 丰满乱子伦码专区| 亚洲av成人精品一区久久| 久久精品国产自在天天线| 免费黄色在线免费观看| 日本黄大片高清| 97热精品久久久久久| 久久99精品国语久久久| 嫩草影院新地址| 99久久精品一区二区三区| 七月丁香在线播放| 蜜桃亚洲精品一区二区三区| 国产爱豆传媒在线观看| 成人午夜高清在线视频| 日韩成人伦理影院| 天堂√8在线中文| 亚洲国产精品久久男人天堂| 成年版毛片免费区| 乱码一卡2卡4卡精品| 免费看a级黄色片| 国产亚洲av嫩草精品影院| 中文字幕av在线有码专区| 七月丁香在线播放| 欧美成人精品欧美一级黄| 日日啪夜夜撸| 亚洲国产色片| 人体艺术视频欧美日本| 一夜夜www| 欧美变态另类bdsm刘玥| 国产精品一及| 亚洲久久久久久中文字幕| 国产亚洲5aaaaa淫片| 麻豆成人午夜福利视频| 99视频精品全部免费 在线| 久久韩国三级中文字幕| 特级一级黄色大片| 国产精品爽爽va在线观看网站| 亚洲在线自拍视频| 久久久成人免费电影| 亚洲国产成人一精品久久久| 亚洲乱码一区二区免费版| 精品人妻熟女av久视频| 成年女人看的毛片在线观看| 91精品伊人久久大香线蕉| 日韩国内少妇激情av| 国产高清国产精品国产三级 | 伦精品一区二区三区| 亚洲性久久影院| 精品国产三级普通话版| 免费看av在线观看网站| 汤姆久久久久久久影院中文字幕 | 男女国产视频网站| 国产成人精品一,二区| 美女国产视频在线观看| 欧美另类亚洲清纯唯美| 亚洲欧美精品专区久久| 午夜爱爱视频在线播放| 成年女人看的毛片在线观看| 欧美成人一区二区免费高清观看| 看十八女毛片水多多多| 99国产精品一区二区蜜桃av| 国产成人freesex在线| 嫩草影院新地址| 欧美+日韩+精品| 亚洲av一区综合| 亚洲国产日韩欧美精品在线观看| 亚洲熟妇中文字幕五十中出| 精品久久久久久久久久久久久| 男女下面进入的视频免费午夜| 三级国产精品片| 日韩精品青青久久久久久| 国产真实乱freesex| 国产极品精品免费视频能看的| 能在线免费看毛片的网站| 夫妻性生交免费视频一级片| 女人十人毛片免费观看3o分钟| 国产精品乱码一区二三区的特点| 偷拍熟女少妇极品色| 欧美日韩国产亚洲二区| 狂野欧美激情性xxxx在线观看| 免费人成在线观看视频色| av又黄又爽大尺度在线免费看 | 国产精品一区二区在线观看99 | 国产精品久久久久久精品电影| 亚洲图色成人| 色噜噜av男人的天堂激情| 精品99又大又爽又粗少妇毛片| 男人舔女人下体高潮全视频| 蜜桃久久精品国产亚洲av| 成人国产麻豆网| 亚洲久久久久久中文字幕| 99久久精品热视频| 欧美成人免费av一区二区三区| 亚洲欧美精品专区久久| 99热全是精品| 免费看a级黄色片| 国产精品一二三区在线看| 蜜臀久久99精品久久宅男| 青春草国产在线视频| 中文欧美无线码| 女人被狂操c到高潮| 村上凉子中文字幕在线| 免费在线观看成人毛片| 国产精品国产三级专区第一集| 国产精品国产三级国产专区5o | 国产老妇伦熟女老妇高清| 永久网站在线| 精品一区二区三区视频在线| 两个人的视频大全免费| 国产成人aa在线观看| 久久欧美精品欧美久久欧美| 一级黄片播放器| 人人妻人人澡欧美一区二区| 久久精品久久久久久久性| 蜜臀久久99精品久久宅男| 啦啦啦啦在线视频资源| av免费观看日本| 国产91av在线免费观看| 婷婷六月久久综合丁香| 高清在线视频一区二区三区 | 青青草视频在线视频观看| 久久精品国产亚洲av涩爱| av视频在线观看入口| 白带黄色成豆腐渣| 一区二区三区免费毛片| 色网站视频免费| 插阴视频在线观看视频| 精品熟女少妇av免费看| 国产片特级美女逼逼视频| 久久久亚洲精品成人影院| 简卡轻食公司| 国产三级中文精品| 乱系列少妇在线播放| 人妻夜夜爽99麻豆av| 在线免费十八禁| 日本-黄色视频高清免费观看| 女人十人毛片免费观看3o分钟| 国产精品一区二区在线观看99 | АⅤ资源中文在线天堂| 美女被艹到高潮喷水动态| 国产精品爽爽va在线观看网站| 欧美日韩国产亚洲二区| 欧美zozozo另类| 夜夜看夜夜爽夜夜摸| 久久国内精品自在自线图片| 麻豆久久精品国产亚洲av| 久久精品夜色国产| 国产精品无大码| 蜜桃久久精品国产亚洲av| 99久久中文字幕三级久久日本| 国产探花在线观看一区二区| 高清日韩中文字幕在线| 老女人水多毛片| 国国产精品蜜臀av免费| 18+在线观看网站| av在线播放精品| 成人亚洲欧美一区二区av| 久久久久网色| 日本五十路高清| 欧美成人免费av一区二区三区| 国产私拍福利视频在线观看| 美女被艹到高潮喷水动态| 视频中文字幕在线观看| 国产白丝娇喘喷水9色精品| 亚洲国产最新在线播放|