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

    橫風(fēng)下高速列車駛?cè)胨淼罆r(shí)瞬態(tài)氣動(dòng)性能研究

    2019-11-22 15:23:03鄧鍔楊偉超尹榮申張平平
    關(guān)鍵詞:橫風(fēng)高速列車隧道

    鄧鍔 楊偉超 尹榮申 張平平

    摘? ?要:針對(duì)列車高速駛?cè)胨淼罆r(shí)流場(chǎng)的三維、非定常及可壓縮湍流等特性,建立了精細(xì)化的隧道-列車-空氣三維CFD數(shù)值模型,對(duì)比分析洞口有無(wú)橫風(fēng)條件下列車駛?cè)胨淼肋^(guò)程中車體周邊的瞬態(tài)流場(chǎng)結(jié)構(gòu)、壓力分布,并研究橫風(fēng)條件下車體的5項(xiàng)氣動(dòng)荷載(氣動(dòng)橫向力、氣動(dòng)升力、傾覆力矩、偏航力矩和點(diǎn)頭力矩)指標(biāo)的瞬變特性以及風(fēng)速和車速變化對(duì)其最大瞬變幅值的影響情況.研究結(jié)果表明:當(dāng)列車在橫風(fēng)環(huán)境下駛?cè)胨淼溃赐獠糠周圀w兩側(cè)流場(chǎng)結(jié)構(gòu)和壓力分布差異顯著,而洞內(nèi)部分差異較小,從而引發(fā)列車進(jìn)洞前后車體壓差突變;列車在進(jìn)洞過(guò)程中,車體的各項(xiàng)氣動(dòng)荷載均存在瞬變效應(yīng),且尾車同時(shí)呈現(xiàn)出傾覆、“上跳”、“蛇形”擺動(dòng)以及“點(diǎn)頭”等行為;風(fēng)速變化對(duì)尾車偏航力矩變化幅值影響較顯著,而車速變化對(duì)頭車偏航力矩變化幅值影響較顯著.

    關(guān)鍵詞:隧道;橫風(fēng);高速列車;氣動(dòng)荷載;瞬變效應(yīng)

    中圖分類號(hào):U25? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1674—2974(2019)09—0069—10

    Abstract: Considering the three dimensional unsteady compressible turbulent characteristics of the flow field when a high-speed train (HST) enters the tunnel, a refined 3D CFD numerical model of tunnel-train-air was built to analyze the transient flow field structure and pressure distribution in the process of HST entering into tunnel under the conditions of crosswind and without crosswind comparatively. The transient change characteristics of the five aerodynamic loads (i.e. aerodynamic side force, lift force, rolling moment, yawing moment and nodding moment) of the train body under crosswind and the influence of wind and vehicle speeds on its maximum variation amplitudes were studied. The main results can be summarized as follow: When a HST enters into tunnel under crosswind environment, flow field structure and pressure distribution on both sides of the train body outside the tunnel have an obvious dissimilarity, and it's smaller inside the tunnel, which causes the sudden change of the train body pressure difference before and after the HST entering the tunnel. In the process of HST entering the tunnel, there are transient change effects in various aerodynamic loads of the train body, and the tail carriage presents capsizing,“jumping”,“snake-like” swing and “nodding” behavior at the same time. The change of the wind speed has a significant effect on the variation amplitude of the tail carriage yawing moment, while the change of the train speed has a significant effect on the variation amplitude of the head carriage yawing moment.

    Key words: tunnel;crosswind;high-speed train (HST);aerodynamic load;transient change effect

    橫風(fēng)是影響高速列車運(yùn)行安全的重要因素,由強(qiáng)橫風(fēng)誘發(fā)的列車失穩(wěn)甚至傾覆等事故在世界各國(guó)均時(shí)有發(fā)生[1-2]. 近年來(lái),國(guó)內(nèi)外諸多學(xué)者以車體的氣動(dòng)荷載為指標(biāo),分別研究了橫風(fēng)條件下,高速列車在平地[3]、路堤[4-5]、橋梁[6-7]以及擋風(fēng)墻[8]等不同基礎(chǔ)設(shè)施上運(yùn)行時(shí)的氣動(dòng)特性. 結(jié)果表明,橫風(fēng)作用下,不同基礎(chǔ)設(shè)施型式對(duì)應(yīng)不同的流場(chǎng)結(jié)構(gòu)和車體氣動(dòng)荷載特性[4]. 當(dāng)列車由一種設(shè)施環(huán)境高速駛?cè)肓硗庖环N設(shè)施環(huán)境時(shí),車體周邊的流場(chǎng)和氣動(dòng)荷載勢(shì)必將發(fā)生突變,而這往往是導(dǎo)致列車發(fā)生安全事故的重要原因.

    相對(duì)于擋風(fēng)墻、路堤和橋梁等結(jié)構(gòu),隧道屬于封閉結(jié)構(gòu),當(dāng)列車由明線橫風(fēng)環(huán)境下的平坦地面高速駛?cè)胨淼罆r(shí),車體氣動(dòng)荷載的變化將更為劇烈,而這又是我國(guó)大風(fēng)地區(qū)列車運(yùn)行常遇的情形. 迄今為止,針對(duì)此類問(wèn)題,國(guó)內(nèi)僅有晉永榮[9]和苗秀娟等[10]學(xué)者對(duì)強(qiáng)橫風(fēng)環(huán)境下高速列車駛出隧道時(shí)的氣動(dòng)性能進(jìn)行了研究,而對(duì)于強(qiáng)橫風(fēng)下列車高速駛?cè)胨淼罆r(shí)車體的氣動(dòng)荷載瞬變特性研究鮮有報(bào)道.

    為此,本文針對(duì)高速列車駛?cè)胨淼肋^(guò)程中流場(chǎng)的三維、非定常及可壓縮湍流等特性,擬采用ANSYS ICEM CFD軟件建立隧道-列車-空氣三維數(shù)值模型,然后導(dǎo)入ANSYS FLUENT進(jìn)行計(jì)算求解,分析了洞口橫風(fēng)條件下列車高速駛?cè)胨淼肋^(guò)程中車體周邊流場(chǎng)結(jié)構(gòu)和壓力的變化特性,以揭示列車在進(jìn)洞過(guò)程中車體的氣動(dòng)力和氣動(dòng)力矩等參數(shù)的瞬變特性及相應(yīng)的氣動(dòng)力學(xué)行為,并研究了風(fēng)速和車速變化對(duì)氣動(dòng)荷載變化幅值的影響規(guī)律,以期為后續(xù)列車運(yùn)行的舒適性和安全性評(píng)估提供參考.

    1? ?數(shù)值模型

    1.1? ?控制方程

    針對(duì)本文所研究的工況特點(diǎn),可知:1)列車高速突入隧道的過(guò)程中車體的位置時(shí)刻變化,屬于瞬態(tài)問(wèn)題;2)馬赫數(shù)(Ma)約等于0.3,可將氣體視為可壓縮流;3)列車周圍流場(chǎng)雷諾數(shù)大于106,處于湍流狀態(tài)。根據(jù)文獻(xiàn)[11]的論述,本文采用RNG k - ε兩方程湍流模型來(lái)模擬列車周邊流場(chǎng)結(jié)構(gòu)的瞬態(tài)演化過(guò)程。其控制方程包括連續(xù)方程、動(dòng)量守恒方程(即Navier-Stokes方程)以及能量守恒方程,分別表示為:

    1.2? ?模型概況

    以我國(guó)高速鐵路常用的CRH3型列車為研究對(duì)象,建立3節(jié)編組的簡(jiǎn)化模型,長(zhǎng)約76 m,如圖1所示,車體表面設(shè)為無(wú)滑移壁面(Wall)邊界條件。計(jì)算模型整體布置如圖2所示(縮尺比例為1 ∶ 1),當(dāng)列車處于初始位置時(shí),車頭鼻尖距隧道洞口90 m;隧道兩端外部大氣場(chǎng)屬于半無(wú)限空間,采用半柱體模擬,柱體直徑約380 m,長(zhǎng)250 m,入口端環(huán)境橫風(fēng)垂直于列車運(yùn)行方向均勻分布(βw = 90°)。大氣場(chǎng)的外邊界設(shè)為壓力遠(yuǎn)場(chǎng)邊界條件(Pressure-far-field),風(fēng)速Vw和風(fēng)向角βw則通過(guò)給定Pressure-far-field邊界中Z和X方向的馬赫數(shù)Ma來(lái)確定。中間為隧道段,隧道長(zhǎng)400 m,隧道橫斷面按照我國(guó)高鐵適用于時(shí)速300~350 km單洞雙線的標(biāo)準(zhǔn)斷面設(shè)置,凈空面積為100 m2。地面以及隧道壁面均設(shè)置為粗糙度0.5的固定壁面邊界條件(Wall).

    計(jì)算區(qū)域均采用結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分,結(jié)合滑移網(wǎng)格法(SMM)和動(dòng)網(wǎng)格法(DMM)[12],將整個(gè)計(jì)算區(qū)域分為靜止網(wǎng)格區(qū)域和動(dòng)網(wǎng)格區(qū)域兩部分(如圖3(a))。其中,靜止區(qū)域保持靜止;動(dòng)網(wǎng)格區(qū)域包含列車及附近的空氣,其橫截面尺寸如圖3(b)所示,運(yùn)用鋪層法實(shí)現(xiàn)列車與隧道間的相對(duì)運(yùn)動(dòng)。具體實(shí)現(xiàn)過(guò)程如下:動(dòng)網(wǎng)格區(qū)域前后端邊界始終保持靜止?fàn)顟B(tài)。在計(jì)算過(guò)程中,區(qū)域最后端緊鄰邊界處的網(wǎng)格層被拉伸,當(dāng)大于設(shè)定的尺寸時(shí),自動(dòng)分裂為兩層;而區(qū)域最前端緊鄰邊界處的網(wǎng)格層被壓縮,當(dāng)小于設(shè)定的尺寸時(shí),自動(dòng)與其緊鄰的網(wǎng)格層合并。列車表面網(wǎng)格及動(dòng)網(wǎng)格區(qū)域內(nèi)部的流體網(wǎng)格類型均設(shè)置為Rigid,即這兩部分網(wǎng)格單元整體向前移動(dòng),不會(huì)出現(xiàn)網(wǎng)格重建及尺寸變化。兩個(gè)區(qū)域的流場(chǎng)信息通過(guò)Interface交換。由于本文主要關(guān)心列車周圍流場(chǎng)情況,為準(zhǔn)確模擬車身周圍渦流的形成、脫落以及附面層效應(yīng),將緊貼車體表面的第1層網(wǎng)格厚度h0設(shè)為0.01 m,并按比率1.1擴(kuò)展至第8層,相應(yīng)的y+值處于50~180的范圍內(nèi),第8層之后繼續(xù)以較大比率擴(kuò)展.列車運(yùn)行方向的網(wǎng)格尺寸范圍為0.01~0.5 m,車身的網(wǎng)格尺寸約為0.4 m,模型總網(wǎng)格單元數(shù)約550萬(wàn),如圖3(c)所示.

    1.3? ?可靠性驗(yàn)證

    為驗(yàn)證本文數(shù)值計(jì)算方法的可靠性,首先依據(jù)文獻(xiàn)[13]的實(shí)車試驗(yàn)場(chǎng)景,另建立8節(jié)編組(長(zhǎng)201.4 m)的列車模型和1 005 m的隧道模型(隧道入口端大氣場(chǎng)長(zhǎng)度取350 m,列車初始位置及隧道斷面仍與圖2保持一致),提取隧道內(nèi)某2個(gè)典型測(cè)點(diǎn)壓力的數(shù)值計(jì)算結(jié)果與文獻(xiàn)[13]中相應(yīng)實(shí)測(cè)結(jié)果進(jìn)行對(duì)比.其中,數(shù)值模型中的測(cè)點(diǎn)位置、隧道長(zhǎng)度、隧道斷面型式以及列車長(zhǎng)度均與實(shí)車試驗(yàn)環(huán)境基本保持一致,洞口無(wú)橫風(fēng)作用.

    如圖4所示,t = 0 s,列車車頭鼻尖剛好抵達(dá)隧道口;t = 14.5 s,列車剛好完全駛離隧道. 其中,測(cè)點(diǎn)1和測(cè)點(diǎn)2分別距隧道入口140 m和200 m,距軌面高度均為1.5 m;列車運(yùn)行速度為300 km/h.

    分析圖4中測(cè)點(diǎn)1和測(cè)點(diǎn)2的數(shù)據(jù)可知,實(shí)測(cè)值與計(jì)算值的波動(dòng)規(guī)律以及峰值大小基本吻合.測(cè)點(diǎn)壓力受列車沖擊波、入口壓縮波、入口膨脹波以及它們的反射波的交替作用而發(fā)生相應(yīng)的變化,即壓縮波引起測(cè)點(diǎn)壓力上升,膨脹波引起測(cè)點(diǎn)壓力下降.對(duì)于測(cè)點(diǎn)1,數(shù)值計(jì)算的測(cè)點(diǎn)壓力最大正峰值為1.26 kPa,相應(yīng)的實(shí)測(cè)值為1.38 kPa,二者相差9%;數(shù)值計(jì)算的測(cè)點(diǎn)壓力最小負(fù)峰值為-1.42 kPa,相應(yīng)的實(shí)測(cè)值為-1.25 kPa,二者相差12%.對(duì)于測(cè)點(diǎn)2,數(shù)值計(jì)算結(jié)果的最大正峰值為1.26 kPa,相應(yīng)的實(shí)測(cè)值為1.45 kPa,二者相差13%;數(shù)值計(jì)算結(jié)果的最小負(fù)峰值為-1.39 kPa,相應(yīng)的實(shí)測(cè)值為-1.35 kPa,二者相差3%.

    此外,為進(jìn)一步驗(yàn)證側(cè)風(fēng)條件下本文數(shù)值模型的可靠性,將Schober等[14]采用ICE3型1 ∶ 15縮尺模型(見(jiàn)圖5)進(jìn)行風(fēng)洞試驗(yàn)所得的尾車氣動(dòng)橫向力和偏航力矩?cái)?shù)據(jù)與本文相應(yīng)的數(shù)值模擬結(jié)果進(jìn)行對(duì)比(CRH3型列車與德國(guó)的ICE3型列車外型基本一致).為確保二者流場(chǎng)雷諾數(shù)一致,對(duì)本文數(shù)值模型另作如下處理:將數(shù)值模型縮尺比例同樣設(shè)置為1 ∶ 15;列車運(yùn)行速度Vt和風(fēng)速Va分別設(shè)定為0 km/h和78 m/s;模型側(cè)面積A和特征高度h分別為10 m2和3 m;最后改變?chǔ)拢?jì)算不同β(10°、 15°、20°、30°、 45° 和 60°)條件下尾車(對(duì)于靜態(tài)列車模型,尾車的研究效果與頭車一致)的氣動(dòng)橫向力系數(shù)Cz和偏航力矩系數(shù)Cmy,如圖6所示.

    由圖6進(jìn)一步分析可知:整體上,本文數(shù)值計(jì)算結(jié)果與Schober的試驗(yàn)結(jié)果變化趨勢(shì)基本一致,即當(dāng)風(fēng)向角β在0°~60°范圍內(nèi)遞增時(shí),二者尾車氣動(dòng)荷載也基本隨之增大.盡管數(shù)值計(jì)算結(jié)果中個(gè)別數(shù)據(jù)與Schober的試驗(yàn)結(jié)果相差較大. 如當(dāng)β = 10°時(shí),數(shù)值計(jì)算結(jié)果Cmy比相應(yīng)的試驗(yàn)結(jié)果大41%.究其原因,可能是由于列車數(shù)值模型中忽略了轉(zhuǎn)向架等細(xì)部構(gòu)造,導(dǎo)致其底面積和側(cè)面積增大. 然而,其余的數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果基本吻合,相差保持在10%左右.

    綜上所述,無(wú)論隧道洞口是否有側(cè)風(fēng)作用,本文所建的數(shù)值模型及其計(jì)算結(jié)果均是可靠的.

    2? ?氣動(dòng)荷載計(jì)算

    高速列車在運(yùn)行過(guò)程中,描述單節(jié)車廂氣動(dòng)荷載指標(biāo)主要有氣動(dòng)阻力Fx、橫向力Fz、氣動(dòng)升力Fy、傾覆力矩Mx、偏航力矩My和點(diǎn)頭力矩Mz . 由于氣動(dòng)阻力Fx主要引起列車運(yùn)行能耗升高,對(duì)行車安全性影響較小,故本文只針對(duì)其余5個(gè)指標(biāo)展開(kāi)研究.圖7給出了列車三節(jié)車廂的氣動(dòng)載荷加載點(diǎn)位置.

    列車在運(yùn)行時(shí),車體表面受到氣動(dòng)壓力和粘性力的共同作用[15],但由于粘性力主要體現(xiàn)在對(duì)氣動(dòng)阻力的影響,對(duì)其余5個(gè)指標(biāo)影響較小,故在計(jì)算這5項(xiàng)氣動(dòng)載荷時(shí),將每節(jié)車廂沿車長(zhǎng)方向劃分為若干段(其中曲線車頭車尾劃分稍密),提取各工況中每段的迎風(fēng)面、背風(fēng)面、頂面和底面平均靜壓(Facet Average Static Pressure)時(shí)程數(shù)據(jù),具體計(jì)算公式如下[16]:

    式中:Fz、Fy分別為單節(jié)車廂所受的氣動(dòng)橫向力和氣動(dòng)升力;Mx、My、Mz分別為單節(jié)車廂所受到的分別繞 x、y、z軸的氣動(dòng)力矩;P uwi? 、P dwi? 、P dsi? 、P usi? 分別為單節(jié)車廂第i段迎風(fēng)面、背風(fēng)面、頂面和底面的平均靜壓; Δxi為單節(jié)車廂第i段的縱向長(zhǎng)度;xi為單節(jié)車廂中第i段中心與該節(jié)車廂力矩中心在x方向上的距離(當(dāng)?shù)趇段中心位于該節(jié)車廂力矩中心的前方,xi取負(fù)值;反之,xi取正值);Thi、Twi分別為單節(jié)車廂第i段對(duì)應(yīng)的車廂高度和寬度;H0為單節(jié)車廂力矩中心距車廂底面的距離,如圖7,本文取1.38 m;Gi為單節(jié)車廂中第i段所受的重力.

    3? ?結(jié)果分析

    3.1? ?流場(chǎng)結(jié)構(gòu)

    當(dāng)列車由明線駛?cè)胨淼罆r(shí),會(huì)形成2種主要的氣流,即:列車風(fēng)和列車進(jìn)入隧道時(shí)的噴射流.當(dāng)洞口存在橫風(fēng)時(shí),橫風(fēng)將改變列車風(fēng)和噴射流的流場(chǎng)結(jié)構(gòu),從而進(jìn)一步影響車體壓力及其氣動(dòng)力學(xué)行為[17].現(xiàn)以軌面上方1.5 m高度處的水平面為例,對(duì)比分析0 m/s和25 m/s橫風(fēng)條件下(βw = 90°),列車以250 km/h的速度駛?cè)胨淼肋^(guò)程中車體附近瞬態(tài)流場(chǎng)結(jié)構(gòu)分布變化特性,如圖8所示,其中t表示列車運(yùn)行時(shí)刻(當(dāng)列車處于初始位置時(shí),t = 0 s).

    由圖8分析可知:

    1)在洞口無(wú)橫風(fēng)條件下,當(dāng)車頭即將進(jìn)洞時(shí),車身兩側(cè)流場(chǎng)結(jié)構(gòu)對(duì)稱,主要表現(xiàn)為:氣流從車頭流向車尾形成以列車為中心的繞流,僅在車尾形成渦旋.隨著列車逐漸進(jìn)入隧道,由于列車對(duì)隧道內(nèi)空氣的排開(kāi)作用,在隧道入口形成與列車運(yùn)行方向相反的噴射流,考慮到雙線隧道中列車與隧道中心并不重合等因素,可認(rèn)為噴射流基本仍呈對(duì)稱結(jié)構(gòu).當(dāng)車尾抵達(dá)隧道入口時(shí),車體繞流主要沿較寬一側(cè)噴出隧道口,而列車行駛一側(cè)相對(duì)狹窄,氣流相對(duì)較少,并對(duì)尾渦流的對(duì)稱性造成局部破壞. 整體上,無(wú)橫風(fēng)條件下,列車駛?cè)胨淼罆r(shí)車體左右兩側(cè)的氣流變化基本對(duì)稱.

    2)在洞口有25 m/s的持續(xù)穩(wěn)定橫風(fēng)條件下,當(dāng)車頭即將進(jìn)洞時(shí),由于列車風(fēng)和橫風(fēng)的聯(lián)合作用,在X方向自車頭至車尾附著一縱向渦流,在洞口也同時(shí)形成一個(gè)豎向的渦流,此渦流是列車進(jìn)洞形成的噴射流和橫風(fēng)共同作用的結(jié)果,可對(duì)車頭壓力造成一定的影響.隨著列車駛?cè)胨淼?,車頭前方的渦流消失,附著車體的渦流不再是從車頭開(kāi)始,而是自隧道入口處開(kāi)始,駛?cè)胨淼啦糠值能圀w附近無(wú)顯著渦流,且車尾后側(cè)渦流向背風(fēng)側(cè)偏轉(zhuǎn).當(dāng)車尾即將抵達(dá)洞口時(shí),車尾渦流區(qū)朝背風(fēng)側(cè)方向偏轉(zhuǎn)更加顯著,并呈現(xiàn)扭轉(zhuǎn)趨勢(shì).

    3.2? ?壓力分布

    車體左右兩側(cè)流場(chǎng)結(jié)構(gòu)的不同必然會(huì)導(dǎo)致車體兩側(cè)附近氣壓的變化. 圖9給出了與圖8對(duì)應(yīng)的Y=1.5 m平面上不同典型時(shí)刻車體周邊瞬態(tài)壓力分布變化.

    由圖9進(jìn)一步分析可知:

    1)當(dāng)洞口無(wú)橫風(fēng)作用時(shí),車體兩側(cè)壓力場(chǎng)分布整體對(duì)稱,且兩側(cè)壓力值變化始終保持同步,不存在顯著壓差.其中明線條件下中間段車身附近壓力基本穩(wěn)定在0 Pa附近;隧道內(nèi)條件下車體附近壓力分布主要由車體位置、隧道壁面的限制、隧道內(nèi)壓縮波和膨脹波等多種因素共同作用,具體表現(xiàn)為:車頭鼻尖處所受正壓較高,頭部近壁側(cè)壓力變化梯度較遠(yuǎn)壁側(cè)稍大,但兩側(cè)壓力分布仍基本對(duì)稱.

    2)在洞口有25 m/s的持續(xù)穩(wěn)定橫風(fēng)作用下,當(dāng)車頭抵達(dá)隧道入口,在背風(fēng)側(cè)車體后約2/3范圍內(nèi)存在負(fù)壓區(qū),此負(fù)壓區(qū)與流場(chǎng)渦流結(jié)構(gòu)對(duì)應(yīng)(見(jiàn)圖8(b)和圖9(b)),負(fù)壓值最低達(dá)到了-3 000 Pa. 隨著列車?yán)^續(xù)駛?cè)胨淼?,進(jìn)入隧道部分車體兩側(cè)壓力趨于平衡;而明線條件下迎風(fēng)側(cè)仍為正壓,背風(fēng)側(cè)仍為負(fù)壓,且車體后方負(fù)壓區(qū)逐漸變短.可見(jiàn),車體兩側(cè)壓差在沿列車運(yùn)行方向上極不均衡.

    3.3? ?氣動(dòng)荷載

    隧道內(nèi)外環(huán)境差異不僅引起車體迎風(fēng)側(cè)和背風(fēng)側(cè)壓力分布的改變,對(duì)整個(gè)車體所受的壓力也將有重要的影響,進(jìn)而會(huì)影響作用于車體的氣動(dòng)力和氣動(dòng)力矩等氣動(dòng)荷載指標(biāo). 圖10給出了在洞口25 m/s的恒定橫風(fēng)條件下(βw = 90°),列車以250 km/h的速度駛?cè)胨淼肋^(guò)程中,三節(jié)車廂的氣動(dòng)橫向力、氣動(dòng)升力、傾覆力矩、偏航力矩和點(diǎn)頭力矩隨時(shí)間的變化過(guò)程,同時(shí)分別給出了對(duì)應(yīng)的氣動(dòng)荷載作用下頭車、中間車和尾車在進(jìn)洞時(shí)的運(yùn)行姿態(tài)變化過(guò)程(圖中t = 1.29 s時(shí)刻列車車頭鼻尖剛好抵達(dá)隧道口,t = 2.39 s時(shí)刻列車剛好完全駛?cè)胨淼溃?

    從圖10中分析可知:

    1)對(duì)于氣動(dòng)橫向力,列車進(jìn)洞前,頭車受到的橫向力最大,中間車次之,尾車最小;列車進(jìn)洞時(shí),頭車和中間車的橫向力直接突降至0 kN,尾車則先迅速降至-15 kN,然后才恢復(fù)至零值附近,對(duì)應(yīng)的尾車將發(fā)生一次整體反向橫移的運(yùn)動(dòng).由于傾覆力矩是通過(guò)橫向力與力臂的乘積求得,故其隨時(shí)間的變化特性與橫向力相同,在此不再贅述,如圖10(d)所示,對(duì)應(yīng)的尾車運(yùn)行姿態(tài)在經(jīng)歷“右傾-左傾”變化后才恢復(fù)平穩(wěn).

    2)對(duì)于氣動(dòng)升力,列車進(jìn)洞前,頭車氣動(dòng)升力在40 kN附近持續(xù)大幅震蕩,而尾車震蕩幅度相對(duì)較小,中間車最為平穩(wěn);進(jìn)洞時(shí),頭車氣動(dòng)升力直接突降,尾車和中間車氣動(dòng)升力將小幅上升后恢復(fù)零值,即中間車和尾車將先后相應(yīng)地產(chǎn)生一個(gè)整體上跳的行為.

    3)對(duì)于偏航力矩,當(dāng)列車在明線運(yùn)行時(shí),中間車偏航力矩最高,頭車次之,尾車最低;列車進(jìn)洞時(shí),三節(jié)車偏航力矩先后發(fā)生突增,且頭車增幅最大,達(dá)344 kN·m. 這是由于列車在進(jìn)洞過(guò)程中,隧道內(nèi)部分車體所受橫向力幾乎為零,而洞外部分車體在橫風(fēng)作用下仍承受巨大橫向力,橫向力沿車廂縱向分布極不均勻,從而使偏航力矩突增,此時(shí)乘客將明顯感覺(jué)車廂劇烈晃動(dòng),列車往往易因偏航力矩過(guò)大而發(fā)生脫軌事故. 隨著車廂完全進(jìn)洞,頭車和中間車的偏航力矩先后直接恢復(fù)至零值,而尾車的尾部將發(fā)生一次向迎風(fēng)側(cè)的反向擺動(dòng)后才恢復(fù)平穩(wěn).

    4)對(duì)于點(diǎn)頭力矩,列車進(jìn)洞前,頭車持續(xù)在100 kN·m附近大幅震蕩,且頭部整體呈上翹姿態(tài),尾車震蕩幅度相對(duì)較小,中間車最平穩(wěn).進(jìn)洞時(shí),頭車點(diǎn)頭力矩經(jīng)歷一個(gè)波谷,相應(yīng)的頭車頭部將產(chǎn)生一個(gè)點(diǎn)頭行為;尾車則經(jīng)歷一個(gè)波峰后才恢復(fù)平穩(wěn),相應(yīng)的尾車尾部也將產(chǎn)生一個(gè)點(diǎn)頭行為.

    綜上所述,列車由明線橫風(fēng)環(huán)境快速駛?cè)胨淼赖倪^(guò)程中,車體的各項(xiàng)氣動(dòng)荷載均表現(xiàn)出顯著的瞬變效應(yīng).在風(fēng)速為25 m/s且車速為250 km/h的運(yùn)行條件下,相對(duì)于中間車和尾車,頭車變化幅值最大,這將嚴(yán)重威脅其運(yùn)行安全性;而尾車由于氣動(dòng)荷載突變作用,引起連續(xù)“右-左”傾覆、“上跳”、“蛇形”擺動(dòng)以及“點(diǎn)頭”等行為同時(shí)發(fā)生,也將在一定程度上降低其乘坐舒適性和運(yùn)行安全性.

    3.4? ?影響因素分析

    由上述分析可知,列車駛?cè)胨淼罆r(shí)氣動(dòng)荷載的瞬變幅度將嚴(yán)重影響列車運(yùn)行的安全性,而影響氣動(dòng)荷載的主要因素有風(fēng)速和車速.下面進(jìn)一步分析風(fēng)速和車速的變化對(duì)氣動(dòng)荷載的影響情況,圖11給出了不同風(fēng)速下(25 m/s、30 m/s、35 m/s、40 m/s)(βw = 90°)列車以250 km/h的速度駛?cè)胨淼肋^(guò)程中各節(jié)車廂氣動(dòng)荷載最大波動(dòng)幅值,圖12給出了在25 m/s的橫風(fēng)環(huán)境下,列車以不同車速(200 km/h、250 km/h、300 km/h、350 km/h)駛?cè)胨淼罆r(shí)各節(jié)車廂的氣動(dòng)荷載最大波動(dòng)幅值.

    由圖11和圖12分析可知:

    1)對(duì)于氣動(dòng)橫向力最大波動(dòng)幅值,隨風(fēng)速增大,三節(jié)車廂均隨之增大,但尾車在風(fēng)速達(dá)到35 m/s以后增大并不顯著;若風(fēng)速不變,車速增大,僅頭車增大,而中間車和和尾車呈微弱減小趨勢(shì).

    2)對(duì)于氣動(dòng)升力最大波動(dòng)幅值,隨風(fēng)速增大,頭車和中間車隨之增大,而尾車增大趨勢(shì)并不明顯;若風(fēng)速不變,車速增大,頭車和中間車微弱增大,尾車仍呈減小趨勢(shì).

    3)對(duì)于偏航力矩最大波動(dòng)幅值,隨著風(fēng)速增大,三節(jié)車均隨之明顯增大,且以尾車增大最為顯著(當(dāng)風(fēng)速為40 km/h時(shí),達(dá)937 kN·m);若風(fēng)速不變,車速增大,三節(jié)車偏航力矩突變幅值變化趨勢(shì)與橫向力相似,但頭車增大尤為顯著.

    4)對(duì)于點(diǎn)頭力矩最大波動(dòng)幅值,隨著風(fēng)速增大,三節(jié)車也均有所增大,且中間車增大較為顯著(當(dāng)風(fēng)速為40 km/h時(shí),達(dá)674 kN·m);若風(fēng)速不變,車速增大,三節(jié)變化趨勢(shì)仍與橫向力相似.

    綜上所述,風(fēng)速變化對(duì)各節(jié)車廂氣動(dòng)荷載變化幅值的影響程度明顯大于車速變化的影響.隨著車速或風(fēng)速的變化,在各項(xiàng)氣動(dòng)荷載中,又以偏航力矩的變化最為顯著.當(dāng)車速不變,隨著風(fēng)速增大,尾車偏航力矩變化幅值增大最為顯著,又因列車運(yùn)行安全性在很大程度上受偏航力矩變化的影響[18],故風(fēng)速變化對(duì)尾車運(yùn)行安全影響較突出;而當(dāng)風(fēng)速不變,隨車速增大,頭車偏航力矩變化幅值增大最為顯著,即車速變化對(duì)頭車運(yùn)行安全影響較突出.

    4? ?結(jié)? ?論

    1)對(duì)于流場(chǎng)結(jié)構(gòu),在洞口橫風(fēng)作用下,洞外部分車體兩側(cè)流場(chǎng)結(jié)構(gòu)差異明顯,背風(fēng)側(cè)附著一縱向渦流.隨著列車駛?cè)胨淼?,附著車體的渦流在隧道入口處開(kāi)始消失,駛?cè)胨淼啦糠值能圀w附近無(wú)顯著渦流,車頭處流場(chǎng)仍呈源流特性.當(dāng)車尾即將抵達(dá)洞口時(shí),車尾渦流區(qū)朝背風(fēng)側(cè)方向偏轉(zhuǎn)更加顯著,并呈現(xiàn)扭轉(zhuǎn)趨勢(shì).

    2)對(duì)于壓力分布,在洞口橫風(fēng)作用下,當(dāng)車頭抵達(dá)隧道入口,在背風(fēng)側(cè)車體后約2/3范圍內(nèi)存在負(fù)壓區(qū).隨著列車?yán)^續(xù)駛?cè)胨淼?,洞?nèi)部分車體兩側(cè)壓力仍趨于平衡,而洞外部分車體兩側(cè)壓差顯著,車體兩側(cè)壓差在縱向上極不均衡.

    3)高速列車由明線橫風(fēng)環(huán)境駛?cè)胨淼赖倪^(guò)程

    中,車體的各項(xiàng)氣動(dòng)荷載均表現(xiàn)出顯著的瞬變效應(yīng).其中尾車在氣動(dòng)荷載大幅瞬變作用下,同時(shí)產(chǎn)生較顯著的連續(xù)“右-左”傾覆、“上跳”、“蛇形”擺動(dòng)以及“點(diǎn)頭”等行為.

    4)風(fēng)速變化對(duì)尾車偏航力矩變化幅值影響較顯著,而車速變化對(duì)頭車偏航力矩變化幅值影響較顯著.

    參考文獻(xiàn)

    [1]? ? BAKER C,CHELI F,ORELLANO A,et al. Cross wind effects on road and rail vehicles[J]. Vehicle System Dynamic,2009,47(8): 983—1022.

    [2]? ? LI X Z,XIAO J,LIU D J,et al. An analytical model for the fluctuating wind velocity spectra of a moving vehicle [J]. Journal of Wind Engineering & Industrial Aerodynamics,2017,164:34—43.

    [3]? ? CHELI C,CORRADI R,ROCCHI D,et al. Wind tunnel tests on train scale models to investigate the effect of infrastructure scenario [J]. Journal of Wind Engineering & Industrial Aerodynamics,2010,98(6/7): 353—362.

    [4]? ? TOMASINI G,GIAPPINO S,CORRADI R,et al. Experimental investigation of the effects of embankment scenario on railway vehicle aerodynamic coefficients [J]. Journal of Wind Engineering & Industrial Aerodynamics,2014,131:59—71.

    [5]? ? ZHOU D,TIAN H Q,THOMPSON M,et al. Numerical and experimental investigations of the flow around a high-speed train on an embankment under side wind conditions [J]. The Aerodynamics of Heavy Vehicles III,2015,79:113—130.

    [6]? ? 李永樂(lè),楊懌,吳夢(mèng)雪,等. 側(cè)向風(fēng)作用下橋上列車交會(huì)過(guò)程的空氣動(dòng)力學(xué)特性[J]. 中國(guó)鐵道科學(xué),2015,36(2):37—44.

    LI Y L,YANG Y,WU M X,et al. Aerodynamic characteristics in the process of two trains passing each other on bridge under cross wind action [J]. China Railway Science,2015,36(2):37—44. (In Chinese)

    [7]? ? BOCCIOLONE M,CHELI F,CORRADI R,et al. Crosswind action on rail vehicles:wind tunnel experimental analyses [J]. Journal of Wind Engineering & Industrial Aerodynamics,2008,96(5):584—610.

    [8]? ? 李田,張繼業(yè),張衛(wèi)華. 橫風(fēng)下高速列車通過(guò)擋風(fēng)墻動(dòng)力學(xué)性能[J]. 鐵道學(xué)報(bào),2012,34(7):30—35.

    LI T,ZHANG J Y,ZHANG W H. Dynamic performance of high-speed train passing windbreak in crosswind [J]. Journal of the China Railway Society,2012,34(7):30—35. (In Chinese)

    [9]? 晉永榮. 大風(fēng)下高速列車單車出隧道空氣動(dòng)力效應(yīng)數(shù)值模擬初探[D]. 蘭州:蘭州交通大學(xué)機(jī)電工程學(xué)院,2014:39—58.

    JIN Y R. A preliminary research on numerical simulation of aerodynamic effect when high-speed train out from tunnel under strong wind condition [D]. Lanzhou:School of Mechanical Engineering,Lanzhou Jiaotong University,2014:39—58. (In Chinese)

    [10]? 苗秀娟,張雷,李志偉,等. 強(qiáng)橫風(fēng)下列車出隧道時(shí)的瞬態(tài)氣動(dòng)性能[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,45(3):958—964.

    MIAO X J,ZHANG L,LI Z W,et al. Aerodynamic performance of train passing out tunnel under cross-wind [J]. Journal of Central South University (Science and Technology),2014,45(3):958—964. (In Chinese)

    [11]? YANG W C,DENG E,LEI M F,et al. Flow structure and aerodynamic behavior evolution during train entering tunnel with entrance in crosswind [J]. Journal of Wind Engineering & Industrial Aerodynamics,2018,175:229—243.

    [12]? 谷正氣,楊濱徽,龔旭,等. 會(huì)車瞬態(tài)氣動(dòng)特性分析與研究[J]. 湖南大學(xué)學(xué)報(bào)(自然科學(xué)版),2010,37(6):27—31.

    GU Z Q,YANG B H,GONG X,et al. Analysis for transient aerodynamic characteristics of two automobiles passing each other[J]. Journal of Hunan University (Natural Sciences),2010,37(6):27—31. (In Chinese)

    [13]? 劉峰,姚松,劉堂紅,等. 高速鐵路隧道壁面氣動(dòng)壓力實(shí)車試驗(yàn)分析[J]. 浙江大學(xué)學(xué)報(bào)(工學(xué)版),2016,50(10):2018—2024.

    LIU F,YAO S,LIU T H,et al. Analysis on aerodynamic pressure of tunnel wall of high-speed railway by full-scale train test [J]. Journal of Zhejiang University (Engineering Science),2016,50(10):2018—2024. (In Chinese)

    [14]? SCHOBER M,WEISE M,ORELLANO A,et al. Wind tunnel investigation of an ICE 3 endcar on three standard ground scenarios [J]. Journal of Wind Engineering & Industrial Aerodynamics,2010,98:345—352.

    [15]? 繆曉郎. 高速列車空氣壓力波作用下橫向振動(dòng)分析[D]. 成都:西南交通大學(xué)機(jī)械工程學(xué)院,2014:37—38.

    MIAO X L. Lateral vibration analysis of high speed train under the air pressure wave[D]. Chengdu:School of Mechanical Engineering,Southwest Jiaotong University,2014:37—38. (In Chinese)

    [16]? 鄧鍔. 橫風(fēng)環(huán)境下高速列車突入隧道時(shí)氣動(dòng)特性及行車安全研究[D]. 長(zhǎng)沙:中南大學(xué)土木工程學(xué)院,2017:45—46.

    DENG E. Study on aerodynamic characteristics and traffic safety of high-speed train burst into tunnel under crosswind environment [D]. Changsha:School of Civil Engineering,Central South University,2017:45—46. (In Chinese)

    [17]? KHIER W,BREUER M,DURST F,et al. Flow structure around trains under side wind conditions:A numerical study [J]. Computers & Fluids,2000,29(2):179—195.

    [18]? 孫建成. 大風(fēng)下高速列車進(jìn)出防風(fēng)明洞及隧道口的氣動(dòng)特性數(shù)值模擬研究[D]. 蘭州:蘭州交通大學(xué)機(jī)電工程學(xué)院,2015:89—92.

    SUN J C. Numerical simulation research on the aerodynamic characteristics of high-speed trains entering and exiting the anti-wind opencut tunnel and the tunnel portal under strong wind [D]. Lanzhou:School of Mechanical Engineering,Lanzhou Jiaotong University,2015:89—92. (In Chinese)

    猜你喜歡
    橫風(fēng)高速列車隧道
    與隧道同行
    神奇的泥巴山隧道
    小讀者(2019年24期)2019-01-10 23:00:37
    橫風(fēng)和軌道不平順聯(lián)合作用下的車輛—軌道系統(tǒng)隨機(jī)分析模型
    氣動(dòng)作用下高速列車響應(yīng)特性研究
    科技資訊(2016年29期)2017-02-28 14:36:58
    黑乎乎的隧道好可怕
    擋風(fēng)墻后動(dòng)車組與普速列車交會(huì)氣動(dòng)特性
    動(dòng)車組橫風(fēng)環(huán)境下的交會(huì)氣動(dòng)效應(yīng)
    高速列車系統(tǒng)集成試驗(yàn)的工藝淺談
    科技傳播(2016年7期)2016-04-28 00:00:02
    動(dòng)車風(fēng)道系統(tǒng)的合理化設(shè)計(jì)
    試論焊后退火工藝對(duì)高速列車轉(zhuǎn)向架焊接接頭組織和性能的影響
    麻豆乱淫一区二区| 精品久久久噜噜| 人妻制服诱惑在线中文字幕| 久久久午夜欧美精品| 久久久久久伊人网av| 成人午夜高清在线视频| 日本黄色片子视频| 看黄色毛片网站| 麻豆久久精品国产亚洲av| 美女国产视频在线观看| 精品久久久精品久久久| 尤物成人国产欧美一区二区三区| 伦理电影大哥的女人| 日本与韩国留学比较| 爱豆传媒免费全集在线观看| 日韩亚洲欧美综合| 亚洲三级黄色毛片| av在线蜜桃| 亚洲av.av天堂| 久久97久久精品| 国产69精品久久久久777片| 亚洲av中文字字幕乱码综合| 在线观看av片永久免费下载| av一本久久久久| 国产精品久久久久久精品电影小说 | 中国美白少妇内射xxxbb| ponron亚洲| 最近最新中文字幕免费大全7| 午夜福利在线观看吧| 午夜福利成人在线免费观看| av天堂中文字幕网| 搡女人真爽免费视频火全软件| 一边亲一边摸免费视频| 2021天堂中文幕一二区在线观| 国产高潮美女av| 日韩欧美 国产精品| 久久久久精品久久久久真实原创| 一级片'在线观看视频| 国产精品爽爽va在线观看网站| 嫩草影院新地址| 在现免费观看毛片| 国产乱来视频区| 直男gayav资源| 午夜免费男女啪啪视频观看| 岛国毛片在线播放| 国产高潮美女av| 精品久久久久久久久亚洲| 日韩av不卡免费在线播放| 久久久亚洲精品成人影院| 美女大奶头视频| 久久久久久伊人网av| a级毛片免费高清观看在线播放| 少妇裸体淫交视频免费看高清| av在线播放精品| 老司机影院成人| 蜜桃亚洲精品一区二区三区| 国产一区有黄有色的免费视频 | 国产精品1区2区在线观看.| 身体一侧抽搐| 男女啪啪激烈高潮av片| 中文天堂在线官网| 汤姆久久久久久久影院中文字幕 | 一级二级三级毛片免费看| 好男人视频免费观看在线| 亚洲国产精品成人久久小说| 午夜精品国产一区二区电影 | 成人鲁丝片一二三区免费| 午夜精品一区二区三区免费看| 国产一区二区三区综合在线观看 | 91久久精品国产一区二区成人| 国产老妇伦熟女老妇高清| 蜜桃亚洲精品一区二区三区| 夜夜爽夜夜爽视频| 99久久九九国产精品国产免费| 欧美三级亚洲精品| 自拍偷自拍亚洲精品老妇| 99久久精品热视频| 国国产精品蜜臀av免费| 亚洲aⅴ乱码一区二区在线播放| 国产精品嫩草影院av在线观看| 熟妇人妻不卡中文字幕| 成人午夜精彩视频在线观看| 91狼人影院| 欧美区成人在线视频| 在线观看av片永久免费下载| 精品久久久久久电影网| 亚洲欧洲日产国产| 国产综合精华液| 禁无遮挡网站| 麻豆av噜噜一区二区三区| 亚洲av在线观看美女高潮| 国产激情偷乱视频一区二区| 精品酒店卫生间| 我的老师免费观看完整版| 中文欧美无线码| 国产精品一及| 精品久久久久久久久av| 视频中文字幕在线观看| 午夜免费激情av| 人妻系列 视频| 极品教师在线视频| 日本三级黄在线观看| 国产欧美日韩精品一区二区| www.av在线官网国产| ponron亚洲| 91aial.com中文字幕在线观看| 亚洲欧美成人综合另类久久久| 国产黄a三级三级三级人| 国产av国产精品国产| 亚洲精品日本国产第一区| 欧美成人精品欧美一级黄| 久久久久久久亚洲中文字幕| 久久久精品免费免费高清| 色综合色国产| 夫妻性生交免费视频一级片| av在线蜜桃| 亚洲精品,欧美精品| 菩萨蛮人人尽说江南好唐韦庄| 床上黄色一级片| 久久97久久精品| 国产黄色视频一区二区在线观看| 亚洲人成网站在线播| 国产激情偷乱视频一区二区| 国产男女超爽视频在线观看| 久久精品人妻少妇| 亚洲美女视频黄频| 久久久久久久午夜电影| 亚洲婷婷狠狠爱综合网| 国产精品一二三区在线看| 老师上课跳d突然被开到最大视频| 人人妻人人看人人澡| 欧美日韩综合久久久久久| 成人毛片a级毛片在线播放| av免费观看日本| 久久久国产一区二区| 久久久精品欧美日韩精品| 国产在视频线精品| 秋霞在线观看毛片| 老女人水多毛片| 日韩,欧美,国产一区二区三区| ponron亚洲| 国产伦精品一区二区三区视频9| 淫秽高清视频在线观看| 18禁裸乳无遮挡免费网站照片| 男人狂女人下面高潮的视频| 精品久久久久久久人妻蜜臀av| 欧美激情久久久久久爽电影| 中文天堂在线官网| 精品一区二区免费观看| 精品一区在线观看国产| 成人无遮挡网站| 建设人人有责人人尽责人人享有的 | 国产人妻一区二区三区在| 蜜桃亚洲精品一区二区三区| 赤兔流量卡办理| 欧美日韩综合久久久久久| 黄片无遮挡物在线观看| 又粗又硬又长又爽又黄的视频| 日韩成人av中文字幕在线观看| 水蜜桃什么品种好| 伊人久久国产一区二区| 国产精品99久久久久久久久| 亚洲电影在线观看av| 在线观看一区二区三区| or卡值多少钱| 亚洲va在线va天堂va国产| 国产69精品久久久久777片| 亚洲在线自拍视频| 色综合色国产| 亚洲真实伦在线观看| av又黄又爽大尺度在线免费看| 欧美zozozo另类| 亚洲欧洲日产国产| 国产黄a三级三级三级人| 国产精品日韩av在线免费观看| 高清av免费在线| 亚洲av福利一区| 亚洲国产欧美人成| 亚洲最大成人av| 最近中文字幕高清免费大全6| 少妇猛男粗大的猛烈进出视频 | 亚洲怡红院男人天堂| 国产 亚洲一区二区三区 | 国产成人精品婷婷| 日本一二三区视频观看| 欧美 日韩 精品 国产| 97超碰精品成人国产| ponron亚洲| 亚洲内射少妇av| 在现免费观看毛片| 99热这里只有是精品50| av福利片在线观看| 久久久久精品久久久久真实原创| 深爱激情五月婷婷| 国产美女午夜福利| 国产免费又黄又爽又色| 肉色欧美久久久久久久蜜桃 | 丰满人妻一区二区三区视频av| 97热精品久久久久久| 免费看日本二区| 免费看不卡的av| 男人舔奶头视频| 国产精品爽爽va在线观看网站| 97精品久久久久久久久久精品| 国产探花在线观看一区二区| 精品国产露脸久久av麻豆 | 午夜福利在线观看吧| 久久久久久久国产电影| 国产成人一区二区在线| 亚洲精品aⅴ在线观看| 爱豆传媒免费全集在线观看| 插阴视频在线观看视频| 亚洲人成网站在线播| 少妇熟女aⅴ在线视频| 亚洲精品成人久久久久久| 舔av片在线| 男人爽女人下面视频在线观看| 亚洲在线自拍视频| 国产精品久久久久久久电影| 亚洲综合色惰| 亚洲婷婷狠狠爱综合网| 建设人人有责人人尽责人人享有的 | 91久久精品国产一区二区成人| 蜜桃亚洲精品一区二区三区| av在线亚洲专区| 亚洲欧美精品自产自拍| 日韩制服骚丝袜av| 毛片一级片免费看久久久久| 男人狂女人下面高潮的视频| 国产午夜精品论理片| 亚洲自拍偷在线| 欧美一区二区亚洲| 久久人人爽人人片av| 欧美高清成人免费视频www| 日韩一本色道免费dvd| 久久久精品94久久精品| 简卡轻食公司| 国产欧美日韩精品一区二区| 男女边摸边吃奶| 丰满乱子伦码专区| 狂野欧美白嫩少妇大欣赏| 尾随美女入室| 2022亚洲国产成人精品| 激情 狠狠 欧美| 最近的中文字幕免费完整| 亚洲人成网站高清观看| 亚洲av日韩在线播放| 国产精品三级大全| 最近2019中文字幕mv第一页| 99久久精品一区二区三区| 国产爱豆传媒在线观看| 夫妻午夜视频| 国产成人aa在线观看| 男女那种视频在线观看| 伦理电影大哥的女人| 热99在线观看视频| 大片免费播放器 马上看| 亚洲av在线观看美女高潮| 精品国产一区二区三区久久久樱花 | 青春草亚洲视频在线观看| 秋霞在线观看毛片| 91精品一卡2卡3卡4卡| 成人毛片60女人毛片免费| 亚洲av成人av| 国产片特级美女逼逼视频| 欧美潮喷喷水| 午夜福利网站1000一区二区三区| 国产高清有码在线观看视频| 嫩草影院入口| 熟妇人妻不卡中文字幕| 精品不卡国产一区二区三区| 亚洲欧美成人综合另类久久久| 国产精品久久久久久av不卡| 亚洲美女搞黄在线观看| 久久久久久久久久人人人人人人| 久久久久精品久久久久真实原创| 一个人观看的视频www高清免费观看| 亚洲av不卡在线观看| h日本视频在线播放| 久久久久久久久大av| 国内精品宾馆在线| 亚洲av中文字字幕乱码综合| 麻豆久久精品国产亚洲av| 精品久久久久久久久av| 久久久久精品性色| 成人一区二区视频在线观看| 高清午夜精品一区二区三区| 国产精品一及| 男人爽女人下面视频在线观看| 少妇的逼水好多| 2021天堂中文幕一二区在线观| 亚洲av一区综合| 午夜精品国产一区二区电影 | 亚洲av不卡在线观看| 国产男女超爽视频在线观看| 免费电影在线观看免费观看| 男女视频在线观看网站免费| 国产精品蜜桃在线观看| 色网站视频免费| 我的女老师完整版在线观看| 国产精品熟女久久久久浪| 久久久精品免费免费高清| 男的添女的下面高潮视频| 欧美精品国产亚洲| 亚洲欧美一区二区三区国产| 人妻一区二区av| 亚洲欧美日韩无卡精品| 91在线精品国自产拍蜜月| 中文字幕亚洲精品专区| av天堂中文字幕网| 日韩三级伦理在线观看| 51国产日韩欧美| 成年女人在线观看亚洲视频 | 三级国产精品欧美在线观看| 麻豆国产97在线/欧美| 久久久久性生活片| 日韩国内少妇激情av| 97人妻精品一区二区三区麻豆| 卡戴珊不雅视频在线播放| 又黄又爽又刺激的免费视频.| 成年女人看的毛片在线观看| 少妇的逼好多水| 在现免费观看毛片| 99久久精品一区二区三区| 2022亚洲国产成人精品| 国产视频首页在线观看| 久久精品夜色国产| 欧美+日韩+精品| 国产不卡一卡二| 2021少妇久久久久久久久久久| 国产探花在线观看一区二区| 男女边摸边吃奶| 久久综合国产亚洲精品| 久久久久久久大尺度免费视频| 老司机影院毛片| 国产欧美另类精品又又久久亚洲欧美| 亚洲最大成人手机在线| 国产激情偷乱视频一区二区| 午夜福利在线在线| 国精品久久久久久国模美| 国产精品蜜桃在线观看| 欧美不卡视频在线免费观看| 国产精品国产三级国产av玫瑰| 麻豆成人午夜福利视频| 久久精品熟女亚洲av麻豆精品 | 菩萨蛮人人尽说江南好唐韦庄| 99久国产av精品国产电影| 97超碰精品成人国产| 欧美日韩精品成人综合77777| 最近手机中文字幕大全| 午夜视频国产福利| 久久国产乱子免费精品| 亚洲婷婷狠狠爱综合网| 午夜免费观看性视频| av国产免费在线观看| 天天躁日日操中文字幕| 亚洲人成网站高清观看| 我的老师免费观看完整版| 日韩强制内射视频| 18禁动态无遮挡网站| 老司机影院毛片| 亚洲内射少妇av| 超碰97精品在线观看| 国产熟女欧美一区二区| 久久久久久久午夜电影| 国产精品一二三区在线看| 精品国产露脸久久av麻豆 | 免费黄色在线免费观看| 国产一区二区三区综合在线观看 | 91久久精品国产一区二区成人| 蜜桃久久精品国产亚洲av| 1000部很黄的大片| av在线天堂中文字幕| 内地一区二区视频在线| 男的添女的下面高潮视频| 综合色丁香网| 国产精品伦人一区二区| 白带黄色成豆腐渣| 亚洲精品乱久久久久久| 三级国产精品片| 色网站视频免费| 久久久久久国产a免费观看| 中文字幕亚洲精品专区| 纵有疾风起免费观看全集完整版 | 久久久久性生活片| 久久鲁丝午夜福利片| 国产国拍精品亚洲av在线观看| 亚洲一区高清亚洲精品| 色网站视频免费| 青青草视频在线视频观看| 九九爱精品视频在线观看| 亚洲电影在线观看av| 成人漫画全彩无遮挡| 久久精品国产鲁丝片午夜精品| 91午夜精品亚洲一区二区三区| 视频中文字幕在线观看| 国产亚洲av片在线观看秒播厂 | 欧美另类一区| 三级国产精品片| 老司机影院成人| 亚洲三级黄色毛片| av播播在线观看一区| 成人无遮挡网站| 中文字幕av在线有码专区| 精品久久久久久久久亚洲| 简卡轻食公司| 国产白丝娇喘喷水9色精品| 亚洲熟妇中文字幕五十中出| 女的被弄到高潮叫床怎么办| 国产有黄有色有爽视频| av一本久久久久| 中文在线观看免费www的网站| 22中文网久久字幕| 欧美日韩视频高清一区二区三区二| 老师上课跳d突然被开到最大视频| 久久99精品国语久久久| 好男人视频免费观看在线| 国产伦理片在线播放av一区| 久久久久久伊人网av| 色视频www国产| 亚洲欧美一区二区三区黑人 | 1000部很黄的大片| 国产一区二区亚洲精品在线观看| 九九爱精品视频在线观看| 久久久午夜欧美精品| 性色avwww在线观看| 高清视频免费观看一区二区 | 日本一本二区三区精品| 亚洲成人中文字幕在线播放| 久久久久网色| 亚洲经典国产精华液单| 日韩强制内射视频| 国产成人精品一,二区| 日韩中字成人| 熟女人妻精品中文字幕| 天天一区二区日本电影三级| 少妇的逼好多水| 成人亚洲精品av一区二区| 免费av不卡在线播放| 免费看美女性在线毛片视频| 九九爱精品视频在线观看| 午夜日本视频在线| 毛片一级片免费看久久久久| 不卡视频在线观看欧美| 国产一区亚洲一区在线观看| 极品教师在线视频| 纵有疾风起免费观看全集完整版 | 精品久久久久久成人av| 神马国产精品三级电影在线观看| 嫩草影院入口| 亚洲精品一二三| 亚洲熟妇中文字幕五十中出| 免费黄网站久久成人精品| 看免费成人av毛片| 亚州av有码| 特大巨黑吊av在线直播| 2022亚洲国产成人精品| 日韩欧美精品免费久久| 日韩国内少妇激情av| 国产精品三级大全| 午夜久久久久精精品| 国产 一区精品| a级毛色黄片| 精品熟女少妇av免费看| 久久韩国三级中文字幕| 欧美变态另类bdsm刘玥| 又大又黄又爽视频免费| 婷婷色av中文字幕| 亚洲精品aⅴ在线观看| 97热精品久久久久久| 人妻一区二区av| 日韩av免费高清视频| 久久精品国产鲁丝片午夜精品| 国产午夜精品久久久久久一区二区三区| 禁无遮挡网站| 伊人久久精品亚洲午夜| av免费在线看不卡| 国产淫片久久久久久久久| 久久久久久久久久成人| 亚洲国产精品专区欧美| 欧美成人一区二区免费高清观看| 免费看不卡的av| 国产亚洲午夜精品一区二区久久 | 国产在视频线精品| 熟妇人妻久久中文字幕3abv| 国产国拍精品亚洲av在线观看| 伊人久久国产一区二区| 精华霜和精华液先用哪个| 熟女人妻精品中文字幕| 一级毛片黄色毛片免费观看视频| 一个人免费在线观看电影| 成人一区二区视频在线观看| 高清日韩中文字幕在线| 啦啦啦啦在线视频资源| 日韩欧美 国产精品| 精品久久久久久久末码| 久久精品国产亚洲av涩爱| 淫秽高清视频在线观看| 秋霞伦理黄片| 美女黄网站色视频| 欧美 日韩 精品 国产| 国产乱人视频| 男插女下体视频免费在线播放| 国产探花极品一区二区| 久久久久久久午夜电影| 99热这里只有是精品50| 热99在线观看视频| 亚洲成人一二三区av| 我的老师免费观看完整版| 插阴视频在线观看视频| 人人妻人人澡人人爽人人夜夜 | 韩国高清视频一区二区三区| 人妻夜夜爽99麻豆av| 毛片女人毛片| 国产精品一区二区在线观看99 | 亚洲一区高清亚洲精品| 麻豆久久精品国产亚洲av| 好男人视频免费观看在线| 亚洲成色77777| 一二三四中文在线观看免费高清| 亚洲精华国产精华液的使用体验| 2018国产大陆天天弄谢| 最近最新中文字幕大全电影3| 久久久久精品久久久久真实原创| 一级av片app| 大话2 男鬼变身卡| 久久久久九九精品影院| 国产精品1区2区在线观看.| 久久精品久久精品一区二区三区| 色5月婷婷丁香| 成人特级av手机在线观看| 联通29元200g的流量卡| 日韩三级伦理在线观看| av在线老鸭窝| 亚洲av.av天堂| 久久精品国产亚洲av天美| 日韩伦理黄色片| 亚洲国产成人一精品久久久| 五月伊人婷婷丁香| 免费看av在线观看网站| 男女边摸边吃奶| 精品一区二区三区人妻视频| 精品久久久久久久末码| 成人亚洲精品一区在线观看 | 欧美97在线视频| 亚洲精品色激情综合| 国产在视频线在精品| 久久久久久久久大av| 精品亚洲乱码少妇综合久久| 久久久久久伊人网av| av黄色大香蕉| 美女国产视频在线观看| 亚洲天堂国产精品一区在线| 七月丁香在线播放| 国产精品麻豆人妻色哟哟久久 | 久久精品国产鲁丝片午夜精品| 久久这里只有精品中国| 亚洲精品国产av蜜桃| eeuss影院久久| 丝瓜视频免费看黄片| 婷婷六月久久综合丁香| 午夜久久久久精精品| 美女黄网站色视频| 亚洲在线自拍视频| 日本猛色少妇xxxxx猛交久久| 亚洲三级黄色毛片| 床上黄色一级片| 80岁老熟妇乱子伦牲交| 成人特级av手机在线观看| 国产成人午夜福利电影在线观看| 18禁在线播放成人免费| 伦精品一区二区三区| 69av精品久久久久久| 亚洲色图av天堂| 一区二区三区乱码不卡18| av.在线天堂| 久久久色成人| 日韩一本色道免费dvd| 非洲黑人性xxxx精品又粗又长| 久久99精品国语久久久| 日本黄色片子视频| 免费观看a级毛片全部| 直男gayav资源| 亚洲自偷自拍三级| 久久这里有精品视频免费| 亚洲欧洲日产国产| 国产午夜精品论理片| 天美传媒精品一区二区| 成年版毛片免费区| 最近手机中文字幕大全| 国产有黄有色有爽视频| 男人爽女人下面视频在线观看| 久久精品人妻少妇| 久久久久性生活片| av免费观看日本| 亚洲美女搞黄在线观看| 免费不卡的大黄色大毛片视频在线观看 | 小蜜桃在线观看免费完整版高清| 国产高清国产精品国产三级 | 亚洲国产最新在线播放| 国产精品久久久久久久电影| 乱系列少妇在线播放| 免费看a级黄色片| 成人亚洲精品av一区二区| 草草在线视频免费看| 午夜久久久久精精品| 少妇的逼好多水| kizo精华| 国产免费视频播放在线视频 | 丰满少妇做爰视频| 免费观看av网站的网址| 最近的中文字幕免费完整| 亚洲成人精品中文字幕电影| 极品教师在线视频| 国产精品久久久久久精品电影|