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

    側(cè)風(fēng)下高速列車臨界傾覆風(fēng)速研究

    2019-12-18 06:21:26劉榮姚松許嬌娥
    關(guān)鍵詞:轉(zhuǎn)向架車速力矩

    劉榮,姚松,許嬌娥

    側(cè)風(fēng)下高速列車臨界傾覆風(fēng)速研究

    劉榮1, 2, 3,姚松1, 2, 3,許嬌娥1, 2, 3

    (1. 軌道交通安全教育部重點(diǎn)實(shí)驗(yàn)室,湖南 長(zhǎng)沙 410075;2. 軌道交通安全關(guān)鍵技術(shù)國際合作聯(lián)合實(shí)驗(yàn)室,湖南 長(zhǎng)沙 410075;3. 軌道交通列車安全保障技術(shù)國家地方聯(lián)合工程研究中心,湖南 長(zhǎng)沙 410075)

    采用數(shù)值模擬方法計(jì)算橫風(fēng)下高速列車的氣動(dòng)力及力矩系數(shù),利用EN14067的五質(zhì)量模型研究橫風(fēng)下車輛臨界傾覆風(fēng)速曲線及不同參數(shù)對(duì)其傾覆的影響。研究結(jié)果表明:臨界傾覆風(fēng)速隨著車速的增大而減小,隨著風(fēng)向角的增大先減小后增大,最小值在 80°左右時(shí)出現(xiàn),且隨著未平衡橫向加速度增大而減小。五質(zhì)量模型中增加考慮的點(diǎn)頭力矩和搖頭力矩對(duì)臨界傾覆風(fēng)速有一定影響,其中若不考慮點(diǎn)頭力矩,設(shè)置車輛臨界風(fēng)速限制時(shí)偏高,對(duì)于車輛運(yùn)行安全性有不利影響。一系懸掛和二系橫向剛度對(duì)傾覆系數(shù)影響不大。隨著二系垂向剛度增加,前轉(zhuǎn)向架輪對(duì)傾覆系數(shù)減小,后轉(zhuǎn)向架增大。橫向止擋間隙增大前后傾覆系數(shù)均增大。當(dāng)抗側(cè)滾扭桿減小到原值60%以上時(shí)傾覆系數(shù)略有增大,幅度不超過10%。車輛質(zhì)心越偏向車輛前端時(shí),前傾覆系數(shù)增加,后傾覆系數(shù)減小。

    高速客車;橫風(fēng)環(huán)境;五質(zhì)量模型;臨界傾覆風(fēng)速

    隨著列車發(fā)展趨向更高速度且輕量化。橫風(fēng)下的列車氣動(dòng)問題越顯突出,列車空氣動(dòng)力學(xué)性能惡化,影響列車橫向穩(wěn)定性,甚至可能使列車脫軌乃至傾覆[1?3]。世界上也多次出現(xiàn)由于強(qiáng)風(fēng)作用導(dǎo)致列車事故的案例[4]。因此,研究大風(fēng)條件下列車運(yùn)行安全速度限值及影響因素是必需的。國內(nèi)外學(xué)者提出了一些能夠評(píng)估鐵路車輛在傾覆風(fēng)險(xiǎn)方面的運(yùn)行安全性的方法。“臨界傾覆風(fēng)速”是評(píng)價(jià)列車橫風(fēng)安全性的重要指標(biāo)之一。國枝正春[5]提出了作用于車輛上外力的靜態(tài)平衡方程,即國枝正春式。日比野有等[6]在國枝式基礎(chǔ)上修正形成了更加詳細(xì)的詳細(xì)解析式法。Thomas等[7]運(yùn)用多體仿真方法模擬研究了軌道彎曲、軌道不平順、車體橫移等影響因素對(duì)列車側(cè)風(fēng)穩(wěn)定性的影響。Baker等[8?10]采用動(dòng)模型風(fēng)洞試驗(yàn)方法研究了列車橫風(fēng)氣動(dòng)特性,并通過靜態(tài)列車模型研究了擋風(fēng)屏障對(duì)列車抗傾覆性能影響。歐洲標(biāo)準(zhǔn)EN14067中提出了計(jì)算列車特征風(fēng)速曲線的方法:三質(zhì)量方法,五質(zhì)量方法(改進(jìn)的準(zhǔn)靜態(tài)方法)以及依靠時(shí)間域的多體系統(tǒng)動(dòng)力學(xué)仿真方法。上述方法中常用的有詳細(xì)解析式法、三質(zhì)量方法和多體仿真方法[11]。多體動(dòng)力學(xué)仿真方法是目前相對(duì)來說考慮因素最全面,計(jì)算臨界風(fēng)速最為準(zhǔn)確的方法;詳細(xì)解析式法和三質(zhì)量模型方法雖然簡(jiǎn)單易推導(dǎo)、考慮因素較少,但是其計(jì)算簡(jiǎn)便,容易定量把握要素影響;五質(zhì)量模型方法相對(duì)來說考慮因素較為全面,更貼近實(shí)際,且相比多體動(dòng)力學(xué)仿真方法要簡(jiǎn)單,因此本文選取較少使用研究的五質(zhì)量模型方法計(jì)算國內(nèi)某型高速列車臨界傾覆風(fēng)速,以及氣動(dòng)力和車輛參數(shù)對(duì)列車運(yùn)行傾覆性能的影響。

    1 列車氣動(dòng)性能數(shù)值模擬

    1.1 高速列車數(shù)值計(jì)算模型建立

    計(jì)算采用流體力學(xué)軟件STAR-CCM+進(jìn)行數(shù)值模擬。為了與風(fēng)洞試驗(yàn)結(jié)果進(jìn)行對(duì)比,采用一節(jié)半列車。模型列車周圍的局部區(qū)域、尾部的尾流區(qū)域、背風(fēng)側(cè)區(qū)域等均加密網(wǎng)格,設(shè)置4個(gè)網(wǎng)格控制體,如圖1所示。采用RANS方法SST-模型,二階中心差分格式擴(kuò)散項(xiàng),對(duì)流項(xiàng)首先使用一階迎風(fēng)格式進(jìn)行離散,再采用二階迎風(fēng)格式[12]。使用SIMPLE法耦合壓力?速度場(chǎng),用迭代法修正壓力。計(jì)算給定空氣來流速度為60 m/s。

    圖1 體網(wǎng)格模型截面

    1.2 風(fēng)洞試驗(yàn)

    為了確定數(shù)值模擬方法對(duì)高速列車側(cè)風(fēng)安全性分析的適用性及檢驗(yàn)數(shù)值模擬的準(zhǔn)確性,將風(fēng)洞試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果進(jìn)行對(duì)比分析。在德國DNW風(fēng)洞試驗(yàn)中,列車模型簡(jiǎn)化了轉(zhuǎn)向架,完整的幾何模型符合EN14067-6:2010要求。列車模型由一節(jié)頭車加一節(jié)中間車的一半組成,縮尺比例為1:15。模型試驗(yàn)中測(cè)定了升力系數(shù),阻力系數(shù),側(cè)向力系數(shù),俯仰力矩系數(shù),側(cè)偏力矩系數(shù)等氣動(dòng)力系數(shù),計(jì)算公式如下:

    式中:為大氣密度;0和0分別為列車參考面積和長(zhǎng)度;i和M分別為列車所受相應(yīng)氣動(dòng)力和力矩;v為作用在車輛上的合成速度,在CFD模擬中也為相對(duì)于車輛的空氣來流速度。其可以分解為車速v和風(fēng)速v為風(fēng)向角,而側(cè)滑角為車輛合成速度與車輛運(yùn)行方向夾角,速度三角形如圖2所示。

    圖2 速度三角形分解

    將模型試驗(yàn)與數(shù)值模擬的結(jié)果進(jìn)行對(duì)比,表1為側(cè)滑角為10°時(shí)的結(jié)果對(duì)比。數(shù)值模擬與模型試驗(yàn)的氣動(dòng)力及力矩系數(shù)的誤差在10%以內(nèi),在工程誤差允許范圍之內(nèi)。存在誤差的主要原因有:試驗(yàn)?zāi)P偷暮?jiǎn)化(轉(zhuǎn)向架),列車模型與地面相對(duì)運(yùn)動(dòng)的差異(地面效應(yīng)模擬),湍流模型以及數(shù)值計(jì)算采用的一些假設(shè)等。由此可見,數(shù)值模擬方法適用于計(jì)算氣動(dòng)載荷系數(shù)。

    表1 數(shù)值計(jì)算與試驗(yàn)結(jié)果對(duì)比

    1.3 CFD計(jì)算結(jié)果

    通過對(duì)列車的氣動(dòng)性能數(shù)值模擬仿真,得到了橫風(fēng)下列車運(yùn)行時(shí)的氣動(dòng)力及力矩系數(shù)。不同側(cè)滑角(0°~90°)的側(cè)風(fēng)環(huán)境下列車氣動(dòng)力系數(shù)如圖3 所示。

    圖3 列車氣動(dòng)力系數(shù)

    由圖3可以看出,除了點(diǎn)頭力矩變化規(guī)律較復(fù)雜,其他力/力矩系數(shù)變化規(guī)律簡(jiǎn)單:側(cè)向力系數(shù)、傾覆力矩系數(shù)的數(shù)值大小均隨著側(cè)滑角的增大先增大后減小,且最大系數(shù)出現(xiàn)在70°左右。傾覆力矩主要與側(cè)向力大小有關(guān),因此兩者變化規(guī)律相似。升力系數(shù)和搖頭力矩系數(shù)的絕對(duì)值大小隨著側(cè)滑角的增大先增大后減小,其數(shù)值最大出現(xiàn)在風(fēng)向角為40°~50°左右。點(diǎn)頭力矩系數(shù)的變化規(guī)律比較復(fù)雜,其數(shù)值最大出現(xiàn)在風(fēng)向角40°左右。

    2 列車臨界傾覆風(fēng)速計(jì)算方法

    2.1 臨界傾覆風(fēng)速定義

    臨界傾覆風(fēng)速是車輛受到橫風(fēng)后開始發(fā)生側(cè)滾時(shí)的風(fēng)速。假如車體受到橫風(fēng)施加的外力,在最不利情況下,迎風(fēng)側(cè)輪重達(dá)到一定臨界值,繼而若風(fēng)速繼續(xù)增大則將引起車輛傾覆,將達(dá)到該臨界值時(shí)的風(fēng)速定義為臨界傾覆風(fēng)速。EN14067-6標(biāo)準(zhǔn)中將輪軌力相對(duì)于靜態(tài)輪軌力的減少值Δ和靜態(tài)輪軌力之比作為傾覆判定指標(biāo),且在減載率大于0.9時(shí)認(rèn)為車輛運(yùn)行有傾覆風(fēng)險(xiǎn)。

    2.2 五質(zhì)量模型方法

    EN14067-6五質(zhì)量模型方法(改進(jìn)的準(zhǔn)靜態(tài)方法)是一種考慮車輛主要部件運(yùn)動(dòng)學(xué)的靜態(tài)方法,不考慮減振器等引起的動(dòng)態(tài)作用,得到的臨界風(fēng)速一般小于多體動(dòng)力學(xué)軟件方法得到的結(jié)果,偏保守。該模型適用于具有兩個(gè)轉(zhuǎn)向架且不考慮相鄰車輛或轉(zhuǎn)向架的車輛。本文研究的高速列車滿足該方法適用條件。

    該方法采用耦合元件連接集中質(zhì)量,模型(見圖4)共包括5個(gè)質(zhì)量塊,車體(CB)、前轉(zhuǎn)向架(BG1)、后轉(zhuǎn)向架(BG2)和前后轉(zhuǎn)向架輪對(duì)(WS1,WS2),一系和二系懸掛位于這些質(zhì)量塊之間,構(gòu)建為具有一定剛度的彈簧單元。此外若有抗側(cè)滾扭桿也應(yīng)構(gòu)建為位于車體和轉(zhuǎn)向架中間的具有繞軸旋轉(zhuǎn)剛度的扭桿。系統(tǒng)坐標(biāo)如圖所示,方向?yàn)檐囕v運(yùn)行方向,為車輛橫向,為垂向。車體具有5個(gè)自由度,沿和方向的平移和繞所有軸旋轉(zhuǎn)的自由度。轉(zhuǎn)向架具有3個(gè)自由度,沿和方向的平移和繞列車運(yùn)行方向旋轉(zhuǎn)的自由度。輪對(duì)具有0個(gè)自由度。系統(tǒng)整體共有11個(gè)自由度。

    在五質(zhì)量模型建立中,需要計(jì)算不同工況下的車輛平衡位置,平衡位置需要由作用在車輛上的外力確定。在列系統(tǒng)方程組時(shí),需要考慮的力和力矩有:質(zhì)量力、彈性力和氣動(dòng)力。質(zhì)量力的計(jì)算與重力加速度和未被平衡加速度有關(guān)。彈性力的確定與變形量相關(guān),除了一二系彈簧力之外,還應(yīng)考慮車體和轉(zhuǎn)向架之間的抗側(cè)滾扭桿和止擋作用,金屬止檔的剛度應(yīng)固定為1×108N/m,在止擋發(fā)生作用前的剛度應(yīng)設(shè)置在0.5 mm的范圍內(nèi)采用拋物線過渡。對(duì)于作用于車體的氣動(dòng)力。不考慮氣動(dòng)阻力作用,計(jì)算如下:

    根據(jù)圖2中的合成速度圖,合成風(fēng)速可以由風(fēng)速、車速、風(fēng)向與線路夾角計(jì)算得到。合成風(fēng)速和側(cè)滑角計(jì)算公式如下:

    考慮系統(tǒng)中所有質(zhì)量體的11個(gè)自由度,設(shè)定車輛系統(tǒng)各剛體的位置矢量,有:

    計(jì)算各個(gè)體上的力和力矩可以構(gòu)建如下關(guān)于車體和轉(zhuǎn)向架力和力矩的系統(tǒng)平衡方程(輪對(duì)自由度為0):

    五質(zhì)量模型的具體數(shù)學(xué)表達(dá)式推導(dǎo)參見文獻(xiàn)[11]。EN14067-6中規(guī)定應(yīng)檢驗(yàn)五質(zhì)量方法的計(jì)算結(jié)果,應(yīng)做檢驗(yàn),與標(biāo)準(zhǔn)附錄中列出的示例車輛的臨界傾覆風(fēng)速曲線的計(jì)算進(jìn)行對(duì)照,計(jì)算風(fēng)向角為90°時(shí)得到的結(jié)果相對(duì)于示例車輛最大偏差應(yīng)在±0.5 m/s內(nèi)。本文中的計(jì)算方法與標(biāo)準(zhǔn)附錄算例進(jìn)行對(duì)比后,結(jié)果符合要求。

    3 列車臨界傾覆風(fēng)速曲線

    利用五質(zhì)量模型方法對(duì)車輛進(jìn)行傾覆計(jì)算,得到了不同風(fēng)速、不同運(yùn)行速度和不同風(fēng)向角下高速列車的傾覆系數(shù)。為防止車輛傾覆,借助高速列車運(yùn)行安全性標(biāo)準(zhǔn),在EN-14067.6中規(guī)定傾覆系數(shù)D <0.9時(shí)車輛安全,以0.9作為臨界值從而確定臨界傾覆風(fēng)速曲線。

    圖5 臨界傾覆風(fēng)速隨車速變化曲線

    圖6 車速為200 km/h時(shí)臨界傾覆風(fēng)速隨風(fēng)向角變化曲線

    圖5是風(fēng)向角為90°時(shí)不同未被平衡橫向加速度下(a=0,a=0.5 m/s2,a=1 m/s2)列車臨界傾覆風(fēng)速隨車速變化曲線??梢钥闯鲭S著車速增加臨界傾覆風(fēng)速下降,在車速在80~100 km/h之間,CWC值迅速下降,而當(dāng)速度大于100 km/h 時(shí),CWC下降速度減慢。且在未被平衡橫向加速度方向?yàn)檎龝r(shí),隨著其增大,列車臨界運(yùn)行安全速度減小,不同未被平衡橫向加速度下CWC曲線變化規(guī)律相同。圖6是車速為200 km/h時(shí)列車臨界傾覆風(fēng)速隨著風(fēng)向角變化曲線,可以看出CWC曲線最小值出現(xiàn)時(shí)的風(fēng)向角并不是垂直于列車運(yùn)行方向,而是在=80°。原因是由于側(cè)風(fēng)作用下列車所受氣動(dòng)系數(shù)的特性,對(duì)車輛傾覆影響最大的是側(cè)向力和傾覆力矩,而其系數(shù)的最大數(shù)值點(diǎn)不是在90°時(shí)出現(xiàn),當(dāng)列車以200 km/h運(yùn)行時(shí),隨著風(fēng)向角增大,相對(duì)于列車的偏航角也增大,但臨界角度對(duì)應(yīng)的系數(shù) 減小。

    4 不同參數(shù)對(duì)列車傾覆特性影響

    4.1 點(diǎn)頭力矩和搖頭力矩影響

    在運(yùn)用較多的詳細(xì)解析式法和三質(zhì)量模型方法中,考慮造成列車傾覆的氣動(dòng)力和力矩時(shí),只考慮了側(cè)向力、升力和傾覆力矩,由于靜態(tài)模型不考慮列車運(yùn)行,運(yùn)行方向的阻力不被考慮。而在五質(zhì)量模型方法中,不僅考慮了這3個(gè)力和力矩,點(diǎn)頭和搖頭力矩也同時(shí)被考慮。因此對(duì)比分析了點(diǎn)頭力矩和搖頭力矩對(duì)臨界傾覆風(fēng)速曲線的影響。

    圖7是不考慮點(diǎn)頭或搖頭力矩時(shí)臨界風(fēng)速變化對(duì)比??梢钥闯?,不考慮搖頭力矩時(shí),臨界風(fēng)速值變化不大。在車速小于200 km/h時(shí),該條件下臨界風(fēng)速小于正??紤]力矩,車速大于200 km/h時(shí),該條件下臨界風(fēng)速稍大于考慮全力矩時(shí),但相差很小。不考慮點(diǎn)頭力矩時(shí),臨界風(fēng)速值變化較大,在車速小于100 km/h時(shí),該條件下臨界風(fēng)速稍小于考慮全力矩時(shí),車速大于100 km/h時(shí),該條件下臨界風(fēng)速大于正??紤]力矩情況,且相差最大值達(dá)2.2 m/s。當(dāng)點(diǎn)頭力矩和搖頭力矩均不考慮時(shí),臨界傾覆風(fēng)速與正常考慮力矩時(shí)相差最大值達(dá)2.5 m/s。總體來說是否考慮搖頭力矩對(duì)臨界傾覆風(fēng)速曲線影響很小,而車速較高時(shí),點(diǎn)頭力矩對(duì)臨界傾覆風(fēng)速曲線有一定影響。說明詳細(xì)解析式法或三質(zhì)量模型方法中不考慮搖頭力矩的作用有一定合理性。但考慮點(diǎn)頭力矩更合理。不考慮點(diǎn)頭力矩時(shí)臨界傾覆風(fēng)速較高,對(duì)車輛運(yùn)行安全有不利影響。

    4.2 車輛參數(shù)影響

    車輛結(jié)構(gòu)參數(shù)尤其懸掛參數(shù)對(duì)列車動(dòng)力學(xué)性能有重要影響,尤其是在側(cè)風(fēng)下運(yùn)行的列車,車輛產(chǎn)生大幅振動(dòng),合適的車體參數(shù)能很好的防止列車偏移,保證列車運(yùn)行安全性。一系懸掛系統(tǒng)能實(shí)現(xiàn)輪對(duì)的彈性定位,二系懸掛系統(tǒng)能提高轉(zhuǎn)向架平穩(wěn)性且抑制蛇形運(yùn)動(dòng)。因此優(yōu)化車輛動(dòng)力學(xué)性能,需要考慮各車輛結(jié)構(gòu)參數(shù)對(duì)其動(dòng)力學(xué)性能的影響。改變某個(gè)參數(shù)并保持其他參數(shù)不變,分析不同車體參數(shù)對(duì)車輛系統(tǒng)動(dòng)力學(xué)性能的影響。

    (a) 不考慮搖頭力矩;(b) 不考慮點(diǎn)頭力矩;(c) 不考慮搖頭力矩和點(diǎn)頭力矩

    由于有些車輛參數(shù)的變化對(duì)臨界傾覆風(fēng)速的影響不明顯,因此以判別車輛臨界狀態(tài)的傾覆系數(shù)作為目標(biāo)結(jié)果??紤]五質(zhì)量車輛模型中的關(guān)鍵特征參數(shù),選定了7個(gè)車輛參數(shù)進(jìn)行分析,分別是:一系垂向/橫向剛度、二系垂向/橫向剛度、二系橫向止擋間隙、抗側(cè)滾扭桿剛度、及車輛偏心。由于各參數(shù)的量綱差異,為統(tǒng)一比較其對(duì)車輛傾覆系數(shù)的影響程度,使指標(biāo)值的變化具有可對(duì)比性,在分析過程中將車輛各原參數(shù)作為比較基準(zhǔn),以參數(shù)變化率為橫坐標(biāo)分析變化情況。設(shè)定各參數(shù)在原值上增減90%,每隔10%取一個(gè)值,其中車輛偏心值過大不符合實(shí)際情況,因此考慮車輛偏心范圍在±1 m以內(nèi)。由于環(huán)境風(fēng)速和列車速度工況較多,取車速為200 km/h,風(fēng)速取25 m/s,未被平衡橫向加速度為0,風(fēng)向角為90°時(shí),計(jì)算車輛參數(shù)變化對(duì)傾覆系數(shù)的影響。在五質(zhì)量模型中,前后轉(zhuǎn)向架輪對(duì)的傾覆系數(shù)被分開考慮,以較大的傾覆系數(shù)值作為傾覆判斷值。

    圖8分別為不同車輛參數(shù)變化對(duì)車輛傾覆系數(shù)的影響??梢钥闯鲲L(fēng)速為25 m/s時(shí),由于此時(shí)合速度作用下的搖頭力矩對(duì)車輛影響較大,因此原車輛參數(shù)下所計(jì)算的列車前轉(zhuǎn)向架輪對(duì)傾覆系數(shù)比后傾覆系數(shù)大,且約高出0.1左右,也說明以單個(gè)轉(zhuǎn)向架的傾覆系數(shù)作為判別指標(biāo)的合理性。可以看出一系懸掛參數(shù)變化對(duì)傾覆系數(shù)無明顯影響,其中一系橫向剛度在其減小量大于70%時(shí)傾覆系數(shù)略有增大。二系懸掛參數(shù)中二系橫向剛度減小對(duì)傾覆系數(shù)的影響很小,其增大時(shí)后轉(zhuǎn)向架輪對(duì)傾覆系數(shù)稍有減小。二系垂向剛度影響較大,隨著二系垂向剛度的減小,前轉(zhuǎn)向架輪對(duì)傾覆系數(shù)增大,后轉(zhuǎn)向架輪對(duì)傾覆系數(shù)減小,且參數(shù)減小到60%時(shí),變化曲線斜率絕對(duì)值增大。橫向止檔自由間隙變化對(duì)傾覆系數(shù)有一定影響,隨著間隙增大前后轉(zhuǎn)向架輪對(duì)傾覆系數(shù)均增大,當(dāng)間隙大于原值時(shí),后轉(zhuǎn)向架輪對(duì)傾覆系數(shù)變化不大??箓?cè)滾扭桿剛度變化的影響不大,在抗側(cè)滾扭桿剛度減小到原值的60%以上時(shí)傾覆系數(shù)略有增大,增幅不超過10%。另外車輛偏心對(duì)車輛傾覆系數(shù)有較大影響,變化率為正時(shí)表示車輛質(zhì)心偏向車輛前端,隨著車輛越向前端偏心,前轉(zhuǎn)向架輪對(duì)的傾覆系數(shù)增加,后轉(zhuǎn)向架輪對(duì)的傾覆系數(shù)減小。

    (a) 一系橫向剛度;(b) 一系垂向剛度;(c) 二系橫向剛度;(d) 二系垂向剛度; (e) 橫向止檔間隙;(f) 抗側(cè)滾扭桿剛度;(g) 車輛偏心

    5 結(jié)論

    1) 臨界傾覆風(fēng)速隨著車速增加而減小,且隨著橫向加速度的增大而減小,但曲線變化規(guī)律相同。隨著風(fēng)向角的增大(0°~90°)臨界傾覆風(fēng)速先增大后略有減小,曲線最小值出現(xiàn)在風(fēng)向角為80°左右。

    2) 五質(zhì)量模型算法較為全面的考慮了靜態(tài)車輛計(jì)算所需的氣動(dòng)力,增加考慮了點(diǎn)頭和搖頭力矩的影響。若不考慮搖頭力矩,當(dāng)車速小于200 km/h時(shí),臨界風(fēng)速小于考慮全力矩時(shí),車速大于200 km/h時(shí),該條件下臨界風(fēng)速變大。點(diǎn)頭力矩對(duì)臨界傾覆風(fēng)速影響較大,在車速小于100 km/h時(shí),該條件下臨界風(fēng)速稍小于考慮全力矩時(shí),車速大于100 km/h時(shí),該條件下臨界風(fēng)速變大,且相差最大值約達(dá) 2.2 m/s。

    3) 通過五質(zhì)量模型算法計(jì)算了不同車輛參數(shù)對(duì)車輛傾覆系數(shù)的影響。一系剛度和二系橫向剛度對(duì)傾覆系數(shù)影響較小。隨著二系垂向剛度增加,前轉(zhuǎn)向架傾覆系數(shù)減小,后轉(zhuǎn)向架傾覆系數(shù)增大。隨著二系橫向止檔間隙增大前后轉(zhuǎn)向架傾覆系數(shù)均增大。抗側(cè)滾扭桿影響較小,當(dāng)其減小到原值的60%以上時(shí)傾覆系數(shù)略有增大,幅度變化不超過10%。若車輛偏心時(shí)質(zhì)心越偏向車輛前端,前轉(zhuǎn)向架傾覆系數(shù)隨之增加,后轉(zhuǎn)向架傾覆系數(shù)減小。說明在計(jì)算列車臨界傾覆風(fēng)速時(shí),應(yīng)重點(diǎn)分析二系垂向剛度、橫向止擋間隙和車輛偏心等因素的影響。

    [1] 高廣軍, 田紅旗, 姚松, 等. 蘭新線強(qiáng)橫風(fēng)對(duì)車輛傾覆穩(wěn)定性的影響[J]. 鐵道學(xué)報(bào), 2004(4): 36?40. GAO Guangjun, TIAN Hongqi, YAO Song, et al. Effect of strong cross-wind on the stability of trains running on the Lanzhou-Xinjiang railway line[J]. Journal of the China Railway Society, 2004(4): 36?40.

    [2] 劉加利, 于夢(mèng)閣, 張繼業(yè), 等. 基于大渦模擬的高速列車橫風(fēng)運(yùn)行安全性研究[J]. 鐵道學(xué)報(bào), 2011, 33(4): 13?21. LIU Jiali, YU Mengge, ZHANG Weihua, et al. Study on running safety of high-speed train under crosswind by large eddy simulation[J]. Journal of the China Railway Society, 2011, 33(4): 13?21.

    [3] Andersson E, Haggstrom J, Sima M, et al. Assessment of train-overturning risk due to strong cross-winds[J]. Proceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit, 2004, 218(3): 213?223.

    [4] 葛盛昌. 新疆鐵路風(fēng)區(qū)大風(fēng)天氣列車安全運(yùn)行辦法研究[J]. 鐵道運(yùn)輸與經(jīng)濟(jì), 2009, 31(8): 32?34, 84. GE Shengchang. Study on train safety operation method in the Xinjiang railway under cross wind[J]. Railway Transport and Economy, 2009, 31(8): 32?34, 84.

    [5] Kunieda M, Theoretical study on the mechanics of overturn of railway rolling stock[R]. Railway Technical Research Report, 1972(793): 1?15.

    [6] 日比野有, 石田弘明.車両の転覆限界に関する靜解析式法[J].鉄道総研報(bào)告, 2003, 17(4): 39?44. Hibino Tamotsu, Ishida Hiroaki. On the static analysis method for the limit of covering[J]. Report of Railway, 2003, 17(4): 39?44.

    [7] Thomas D, Diedrichs B, Berg M, et al. Dynamics of a high-speed rail vehicle negotiating curves at unsteady crosswind[J]. Journal of Rail and Rapid Transit, 2010(224): 567?579.

    [8] Baker C J. The wind tunnel determination of crosswind forces and moments on a high speed train[M]. Berlin: Springer- Verlag Berlin Heidelberg, 2002.

    [9] Baker C J. The behaviour of road vehicles in unsteady cross winds[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1993(49): 439?448.

    [10] Baker C J. Ground vehicles in high crosswinds parts I: Steady aerodynamic forces[J]. Journal of Fluids and Structures, 1991(5): 69?90.

    [11] EN 14067-6: 2010, Railway applications-aerodynamics- part 6: Requirements and test procedures for cross wind assessment[S].

    [12] 王福軍. 計(jì)算流體動(dòng)力學(xué)分析[M]. 北京: 清華大學(xué)出版社, 2004: 210?215. WANG Fujun. Analysis of computational fluid dynamics [M]. Beijing: Tsinghua University Press, 2004: 210?215.

    Study on critical overturning wind speed of high-speed train under cross wind

    LIU Rong1, 2, 3, YAO Song1, 2, 3, XU Jiaoe1, 2, 3

    (1. Key Laboratory of Traffic Safety on Track, Central South University, Ministry of Education, Changsha 410075, China;2. Joint International Research Laboratory of Key Technology for Rail Traffic Safety, Changsha 410075, China; 3. National & Local Joint Engineering Research Center of Safety Technology for Rail Vehicle, Changsha 410075, China)

    The aerodynamic coefficients of high-speed train under crosswind were calculated by numerical simulation method. The five-mass model of EN14067 was used to study the critical overturning wind speed curve of vehicles under crosswind and the influence of different parameters on its overturning. Results show that: the critical overturning wind speed decreases with the increase of vehicle speed, critical overturning wind speed first decreases and then increases with the increase of wind direction angle, and the minimum value appears at about 80°, and decreases with the increase of unbalanced lateral acceleration. The shaking moment and nodding moment have a certain impact on the critical overturning wind speed, If the nodding moment is not taken into account, the critical wind speed limit of the vehicle will be set higher, which has a negative impact on the safety of the train operation. The primary suspension parameters and the second lateral stiffness have little effect on the overturning coefficient. With the increase of the vertical stiffness of the secondary system, the overturning coefficient of the front bogie wheel set decreases and that of the rear bogie increases. With the increase of lateral stop clearance both overturning coefficients increase. When the anti-roll torsion bar is reduced to more than 60% of the original value, the overturning coefficient increases slightly, and the range is not more than 10%. When the centroid of the vehicle leans towards the front end of the vehicle, the overturning coefficient of the front bogie wheel set increases and the overturning coefficient of the rear bogie wheel set decreases.

    high-speed train; crosswind environment; five-mass model; critical overturning wind speed

    U271.1

    A

    1672 ? 7029(2019)11? 2643 ? 08

    10.19713/j.cnki.43?1423/u.2019.11.001

    2019?02?26

    國家重點(diǎn)研發(fā)計(jì)劃資助項(xiàng)目(2016YFB1200506-03)

    姚松(1975?),男,湖北公安人,教授,博士,從事軌道交通安全研究;E?mail:song_yao@csu.edu.cn

    (編輯 蔣學(xué)東)

    猜你喜歡
    轉(zhuǎn)向架車速力矩
    20t軸重米軌轉(zhuǎn)向架裝用搖枕的研制
    基于SPS模式的轉(zhuǎn)向架軸箱裝配線仿真研究
    基于結(jié)構(gòu)應(yīng)力的轉(zhuǎn)向架構(gòu)架焊縫疲勞強(qiáng)度研究
    2012款奔馳R300車修改最高車速限制
    發(fā)動(dòng)機(jī)阻力矩計(jì)算和起動(dòng)機(jī)介紹
    山東青年(2016年12期)2017-03-02 18:22:48
    小型力矩電機(jī)波動(dòng)力矩的測(cè)量
    彈性負(fù)載力矩下舵偏轉(zhuǎn)角度的測(cè)量方法
    209 P型轉(zhuǎn)向架軸箱彈簧斷裂的原因及預(yù)防
    基于D-最優(yōu)化理論的陀螺儀力矩反饋測(cè)試法
    北京現(xiàn)代途勝車車速表不工作
    色哟哟哟哟哟哟| 国产亚洲精品久久久久久毛片| 一级av片app| 久久久国产成人免费| 亚洲av日韩精品久久久久久密| 欧美在线黄色| av福利片在线观看| 婷婷丁香在线五月| 国产精品三级大全| 国产精品一区二区三区四区久久| 亚洲国产精品999在线| 人人妻人人澡欧美一区二区| 熟女电影av网| 麻豆成人av在线观看| 赤兔流量卡办理| www.熟女人妻精品国产| 欧美xxxx性猛交bbbb| 日本一二三区视频观看| 国产av不卡久久| 日本精品一区二区三区蜜桃| 久久亚洲精品不卡| av天堂在线播放| 国产黄片美女视频| 中文字幕av在线有码专区| 日韩欧美在线二视频| 麻豆av噜噜一区二区三区| 亚洲五月天丁香| 欧美性猛交╳xxx乱大交人| 精品久久国产蜜桃| 国产欧美日韩精品亚洲av| 人人妻人人澡欧美一区二区| 在线观看免费视频日本深夜| 午夜两性在线视频| av专区在线播放| 尤物成人国产欧美一区二区三区| 99热这里只有是精品在线观看 | 老司机午夜十八禁免费视频| 日本一二三区视频观看| 一区二区三区免费毛片| 亚洲av一区综合| 免费高清视频大片| 久久人人爽人人爽人人片va | 欧美性猛交黑人性爽| 精品久久久久久久久av| 一a级毛片在线观看| 成人国产一区最新在线观看| 国产精品国产高清国产av| 看片在线看免费视频| 精品午夜福利视频在线观看一区| 好男人在线观看高清免费视频| 夜夜爽天天搞| 色综合站精品国产| 99热只有精品国产| 亚洲自偷自拍三级| АⅤ资源中文在线天堂| 男女之事视频高清在线观看| 日本精品一区二区三区蜜桃| 少妇被粗大猛烈的视频| 人人妻,人人澡人人爽秒播| 国产一区二区三区视频了| 久久精品国产自在天天线| 淫秽高清视频在线观看| 性插视频无遮挡在线免费观看| 成人av在线播放网站| 午夜久久久久精精品| netflix在线观看网站| 成人特级黄色片久久久久久久| 一进一出抽搐gif免费好疼| 久久精品国产自在天天线| 综合色av麻豆| 亚洲五月婷婷丁香| 欧美日韩亚洲国产一区二区在线观看| 免费在线观看日本一区| 人妻制服诱惑在线中文字幕| 亚洲电影在线观看av| 午夜免费成人在线视频| 亚洲18禁久久av| 欧美xxxx黑人xx丫x性爽| 亚洲成av人片免费观看| 桃色一区二区三区在线观看| 中文字幕人妻熟人妻熟丝袜美| 在线播放无遮挡| 亚洲电影在线观看av| 免费看a级黄色片| 五月玫瑰六月丁香| 亚洲内射少妇av| 免费大片18禁| av视频在线观看入口| 嫩草影院新地址| 麻豆成人午夜福利视频| 精品一区二区三区视频在线| 90打野战视频偷拍视频| 岛国在线免费视频观看| 18禁黄网站禁片午夜丰满| 直男gayav资源| 天堂网av新在线| 特级一级黄色大片| 老熟妇仑乱视频hdxx| 色噜噜av男人的天堂激情| 五月玫瑰六月丁香| 黄片小视频在线播放| 黄色视频,在线免费观看| 国产欧美日韩一区二区精品| 又爽又黄无遮挡网站| 欧美一级a爱片免费观看看| 亚洲专区国产一区二区| 色综合亚洲欧美另类图片| 天堂动漫精品| 久久久久久久午夜电影| 麻豆国产97在线/欧美| 制服丝袜大香蕉在线| 日本a在线网址| 精品熟女少妇八av免费久了| 天堂影院成人在线观看| av专区在线播放| 欧美精品啪啪一区二区三区| 又黄又爽又免费观看的视频| h日本视频在线播放| 大型黄色视频在线免费观看| 九九久久精品国产亚洲av麻豆| 欧美成人性av电影在线观看| 亚洲成av人片免费观看| 可以在线观看的亚洲视频| 日韩欧美在线二视频| 直男gayav资源| 亚洲不卡免费看| 亚洲不卡免费看| 久久精品国产亚洲av天美| 亚洲av电影不卡..在线观看| 成熟少妇高潮喷水视频| 欧美高清性xxxxhd video| 校园春色视频在线观看| 欧美极品一区二区三区四区| 国产成人欧美在线观看| 国产视频内射| 少妇人妻精品综合一区二区 | ponron亚洲| 成人av一区二区三区在线看| 亚洲中文字幕日韩| 夜夜看夜夜爽夜夜摸| 国产视频一区二区在线看| 亚洲人成网站在线播放欧美日韩| 变态另类成人亚洲欧美熟女| .国产精品久久| 亚洲成人免费电影在线观看| 国产精品影院久久| 国产不卡一卡二| 国产精品久久久久久精品电影| 麻豆国产97在线/欧美| 一个人免费在线观看的高清视频| 久久精品综合一区二区三区| 嫩草影院精品99| 日日干狠狠操夜夜爽| 国产亚洲av嫩草精品影院| 色视频www国产| 欧美极品一区二区三区四区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | a在线观看视频网站| 精品免费久久久久久久清纯| 精品免费久久久久久久清纯| 99国产精品一区二区蜜桃av| 搡老妇女老女人老熟妇| 少妇裸体淫交视频免费看高清| 18禁在线播放成人免费| 美女xxoo啪啪120秒动态图 | 欧美中文日本在线观看视频| 麻豆av噜噜一区二区三区| 国产精品久久久久久久久免 | .国产精品久久| 亚洲精品在线美女| 午夜福利18| 精品国产三级普通话版| 久久婷婷人人爽人人干人人爱| 国产91精品成人一区二区三区| 国产精品,欧美在线| 欧美性猛交黑人性爽| 国产黄色小视频在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲人成电影免费在线| 久9热在线精品视频| 狂野欧美白嫩少妇大欣赏| 婷婷丁香在线五月| 99久国产av精品| 99久国产av精品| 高清在线国产一区| 丁香六月欧美| 国产熟女xx| 色5月婷婷丁香| 亚洲无线观看免费| 国产免费一级a男人的天堂| 丁香欧美五月| 亚洲内射少妇av| 一个人免费在线观看的高清视频| 偷拍熟女少妇极品色| 听说在线观看完整版免费高清| 极品教师在线免费播放| 亚洲av第一区精品v没综合| 国模一区二区三区四区视频| 男人和女人高潮做爰伦理| 简卡轻食公司| 国产成人欧美在线观看| 直男gayav资源| 99精品在免费线老司机午夜| 好看av亚洲va欧美ⅴa在| 免费无遮挡裸体视频| 亚洲欧美日韩卡通动漫| 日本免费a在线| 日日干狠狠操夜夜爽| 久久久色成人| 桃色一区二区三区在线观看| 国产 一区 欧美 日韩| 日日摸夜夜添夜夜添av毛片 | 国产精品野战在线观看| 两性午夜刺激爽爽歪歪视频在线观看| av天堂中文字幕网| 一区二区三区高清视频在线| 最好的美女福利视频网| 欧美黄色片欧美黄色片| 国产精品精品国产色婷婷| 亚洲五月婷婷丁香| 一进一出抽搐动态| 亚洲成a人片在线一区二区| 国产成人啪精品午夜网站| 午夜免费男女啪啪视频观看 | 啪啪无遮挡十八禁网站| 国产成人aa在线观看| 在线播放国产精品三级| 免费在线观看日本一区| 久久久久国产精品人妻aⅴ院| 91字幕亚洲| 亚洲综合色惰| 亚洲av电影在线进入| 亚洲精品影视一区二区三区av| 国产麻豆成人av免费视频| 夜夜夜夜夜久久久久| 午夜视频国产福利| 久久久久亚洲av毛片大全| 露出奶头的视频| 亚洲人成电影免费在线| 搡老熟女国产l中国老女人| 欧美日本亚洲视频在线播放| 一级毛片久久久久久久久女| 国产成人欧美在线观看| 日本熟妇午夜| 久久精品91蜜桃| 午夜a级毛片| 丰满乱子伦码专区| 成人无遮挡网站| 久99久视频精品免费| 国产精品电影一区二区三区| 国产精品久久久久久久电影| 国产一区二区在线观看日韩| 国产一区二区在线av高清观看| 日本免费a在线| 亚洲av免费高清在线观看| 国产精品永久免费网站| 脱女人内裤的视频| 亚洲av中文字字幕乱码综合| 999久久久精品免费观看国产| 久久精品影院6| av女优亚洲男人天堂| 免费人成在线观看视频色| 看片在线看免费视频| 欧美激情国产日韩精品一区| .国产精品久久| 国产精品影院久久| 有码 亚洲区| 男人和女人高潮做爰伦理| 国产在线精品亚洲第一网站| 夜夜看夜夜爽夜夜摸| 久久人人爽人人爽人人片va | 亚洲美女搞黄在线观看 | 91在线精品国自产拍蜜月| 亚洲一区高清亚洲精品| netflix在线观看网站| 国产精品一区二区性色av| 免费av毛片视频| www.www免费av| 久久精品国产亚洲av天美| 国产精品亚洲一级av第二区| 国产高清视频在线观看网站| 中亚洲国语对白在线视频| 国产在视频线在精品| 动漫黄色视频在线观看| 88av欧美| 国产亚洲av嫩草精品影院| 九色国产91popny在线| 亚洲第一区二区三区不卡| 亚洲av免费在线观看| 亚洲五月天丁香| 欧美极品一区二区三区四区| 99国产精品一区二区三区| 自拍偷自拍亚洲精品老妇| 亚洲av免费在线观看| 最后的刺客免费高清国语| 亚洲不卡免费看| 一区二区三区高清视频在线| 国产精品三级大全| 天天躁日日操中文字幕| 久久久久精品国产欧美久久久| 亚洲一区二区三区色噜噜| 国产成人福利小说| 国产真实乱freesex| 亚洲久久久久久中文字幕| 人妻制服诱惑在线中文字幕| 亚洲片人在线观看| 色播亚洲综合网| 九色国产91popny在线| 国产精品国产高清国产av| 一进一出好大好爽视频| 好男人在线观看高清免费视频| 91狼人影院| 婷婷六月久久综合丁香| 最近在线观看免费完整版| 可以在线观看毛片的网站| 日韩人妻高清精品专区| 麻豆国产97在线/欧美| 中国美女看黄片| 亚洲avbb在线观看| 黄色一级大片看看| 国产高清视频在线播放一区| 欧美黄色片欧美黄色片| a在线观看视频网站| 久久精品综合一区二区三区| 日本 欧美在线| 99热这里只有是精品在线观看 | 精品久久久久久久久亚洲 | a级毛片a级免费在线| 欧美xxxx黑人xx丫x性爽| 久久久国产成人精品二区| 黄色配什么色好看| 精品人妻视频免费看| 少妇裸体淫交视频免费看高清| 老司机午夜福利在线观看视频| 少妇高潮的动态图| 国产精品av视频在线免费观看| 亚洲av第一区精品v没综合| 老司机午夜福利在线观看视频| 757午夜福利合集在线观看| 一区福利在线观看| 啪啪无遮挡十八禁网站| 国产又黄又爽又无遮挡在线| 国产亚洲精品久久久久久毛片| 在线免费观看的www视频| 俺也久久电影网| 国产在线精品亚洲第一网站| 欧美日韩综合久久久久久 | 亚洲成av人片在线播放无| 亚洲专区中文字幕在线| 九九热线精品视视频播放| 国产亚洲欧美98| 欧美绝顶高潮抽搐喷水| av国产免费在线观看| 亚洲精品影视一区二区三区av| 校园春色视频在线观看| 精品不卡国产一区二区三区| 天堂动漫精品| 一进一出抽搐动态| 一本综合久久免费| 亚洲美女视频黄频| 中国美女看黄片| 国产精品久久电影中文字幕| 好看av亚洲va欧美ⅴa在| 午夜影院日韩av| 一区二区三区四区激情视频 | 毛片一级片免费看久久久久 | 搡老妇女老女人老熟妇| 特大巨黑吊av在线直播| 国产精品人妻久久久久久| 亚洲精品在线观看二区| 免费观看人在逋| 精品久久久久久成人av| 日本成人三级电影网站| 男女床上黄色一级片免费看| 最近中文字幕高清免费大全6 | 欧美性猛交╳xxx乱大交人| xxxwww97欧美| netflix在线观看网站| 亚洲成人精品中文字幕电影| 狠狠狠狠99中文字幕| 蜜桃久久精品国产亚洲av| 欧美日韩综合久久久久久 | 精品不卡国产一区二区三区| 久久99热这里只有精品18| 一级作爱视频免费观看| 在线十欧美十亚洲十日本专区| 亚洲av.av天堂| 久99久视频精品免费| 欧美日韩国产亚洲二区| 此物有八面人人有两片| 超碰av人人做人人爽久久| 成熟少妇高潮喷水视频| 亚洲中文字幕日韩| 淫秽高清视频在线观看| bbb黄色大片| 乱码一卡2卡4卡精品| 国产乱人伦免费视频| 亚洲国产精品999在线| 国产精品久久久久久人妻精品电影| 蜜桃亚洲精品一区二区三区| 欧美日本视频| ponron亚洲| 国产成人福利小说| 欧美黄色淫秽网站| 国产精品久久久久久久电影| 超碰av人人做人人爽久久| 亚洲av成人精品一区久久| av欧美777| 国产视频内射| 亚洲成人久久性| 99久久无色码亚洲精品果冻| 亚洲午夜理论影院| 长腿黑丝高跟| 国产 一区 欧美 日韩| 亚洲欧美日韩无卡精品| 18禁裸乳无遮挡免费网站照片| av专区在线播放| 色吧在线观看| 蜜桃亚洲精品一区二区三区| 91午夜精品亚洲一区二区三区 | 俄罗斯特黄特色一大片| 老司机午夜福利在线观看视频| 一区二区三区高清视频在线| 禁无遮挡网站| 成人美女网站在线观看视频| 色哟哟哟哟哟哟| 给我免费播放毛片高清在线观看| 在线观看66精品国产| 制服丝袜大香蕉在线| 丰满人妻熟妇乱又伦精品不卡| 成人国产一区最新在线观看| 嫩草影院精品99| 精品久久久久久久久久免费视频| 亚洲真实伦在线观看| 1024手机看黄色片| 久久香蕉精品热| 最近最新免费中文字幕在线| 搡老熟女国产l中国老女人| a级毛片a级免费在线| 国产色婷婷99| 热99在线观看视频| 麻豆av噜噜一区二区三区| 久久精品人妻少妇| 日韩亚洲欧美综合| 亚洲人成网站在线播放欧美日韩| 丝袜美腿在线中文| 色5月婷婷丁香| 亚洲av日韩精品久久久久久密| 国产 一区 欧美 日韩| 热99在线观看视频| 欧美成人免费av一区二区三区| av天堂在线播放| av女优亚洲男人天堂| 亚洲电影在线观看av| 亚洲精品影视一区二区三区av| 伊人久久精品亚洲午夜| 亚洲片人在线观看| 免费看日本二区| 自拍偷自拍亚洲精品老妇| 日韩有码中文字幕| 日本黄大片高清| 18禁裸乳无遮挡免费网站照片| 国产白丝娇喘喷水9色精品| 一个人免费在线观看电影| 在线观看舔阴道视频| 亚洲av日韩精品久久久久久密| 亚洲三级黄色毛片| 国产欧美日韩一区二区三| 国产精品久久久久久久久免 | 免费大片18禁| 国产伦精品一区二区三区四那| 美女黄网站色视频| 色av中文字幕| 最新中文字幕久久久久| 午夜精品一区二区三区免费看| 午夜两性在线视频| 免费观看的影片在线观看| 又粗又爽又猛毛片免费看| 日本黄色片子视频| 久久久久久久久大av| 国产伦精品一区二区三区四那| 男女之事视频高清在线观看| 日本黄色视频三级网站网址| 日韩亚洲欧美综合| 亚洲欧美日韩无卡精品| 级片在线观看| 一个人免费在线观看的高清视频| 久久精品国产亚洲av涩爱 | 白带黄色成豆腐渣| 麻豆成人午夜福利视频| 亚洲精品亚洲一区二区| 成熟少妇高潮喷水视频| 十八禁网站免费在线| 国产精品亚洲美女久久久| 国模一区二区三区四区视频| 日韩精品青青久久久久久| 99久久精品国产亚洲精品| 最新中文字幕久久久久| 成人欧美大片| 午夜亚洲福利在线播放| 99视频精品全部免费 在线| 精品欧美国产一区二区三| 精品人妻1区二区| 亚洲欧美精品综合久久99| 我的女老师完整版在线观看| 久久久久久久久中文| 国产伦人伦偷精品视频| 午夜两性在线视频| 一级黄片播放器| 草草在线视频免费看| 亚洲三级黄色毛片| 男女之事视频高清在线观看| 欧美日韩黄片免| 欧美3d第一页| www日本黄色视频网| 老司机福利观看| 少妇丰满av| 亚洲在线观看片| 精品久久国产蜜桃| 在线国产一区二区在线| 亚洲无线在线观看| 三级男女做爰猛烈吃奶摸视频| 简卡轻食公司| 亚洲片人在线观看| 午夜视频国产福利| 国产精品久久视频播放| 国产精品影院久久| 午夜久久久久精精品| 最近视频中文字幕2019在线8| 亚洲av五月六月丁香网| 一本一本综合久久| 国产精品久久久久久亚洲av鲁大| 午夜福利成人在线免费观看| av欧美777| 尤物成人国产欧美一区二区三区| 91麻豆精品激情在线观看国产| 日韩欧美国产在线观看| 他把我摸到了高潮在线观看| 淫妇啪啪啪对白视频| 国产黄a三级三级三级人| 波多野结衣巨乳人妻| 午夜福利免费观看在线| 熟妇人妻久久中文字幕3abv| 99久久无色码亚洲精品果冻| 国产精品1区2区在线观看.| 欧美黑人巨大hd| 成年版毛片免费区| 中国美女看黄片| 亚洲第一区二区三区不卡| 精品熟女少妇八av免费久了| 午夜免费成人在线视频| 欧美成人免费av一区二区三区| 亚洲欧美日韩无卡精品| xxxwww97欧美| 亚洲精品在线观看二区| 无人区码免费观看不卡| 国产精品一区二区性色av| 国产欧美日韩一区二区精品| av福利片在线观看| 日韩av在线大香蕉| 亚洲成av人片在线播放无| 乱人视频在线观看| 免费一级毛片在线播放高清视频| 国产私拍福利视频在线观看| av黄色大香蕉| 一a级毛片在线观看| 一卡2卡三卡四卡精品乱码亚洲| 成人国产一区最新在线观看| 一级作爱视频免费观看| ponron亚洲| 国产真实乱freesex| 国产高潮美女av| 久99久视频精品免费| 日本黄色视频三级网站网址| 18禁裸乳无遮挡免费网站照片| 女同久久另类99精品国产91| 性色avwww在线观看| 国产老妇女一区| 高清毛片免费观看视频网站| 亚洲人成网站在线播放欧美日韩| 欧美黑人巨大hd| 小说图片视频综合网站| 国产精品久久视频播放| 国产老妇女一区| 女人十人毛片免费观看3o分钟| 欧美日韩国产亚洲二区| 免费大片18禁| 好看av亚洲va欧美ⅴa在| 亚洲av.av天堂| 精品人妻一区二区三区麻豆 | eeuss影院久久| 亚洲av不卡在线观看| 黄色女人牲交| 欧美高清性xxxxhd video| 久久这里只有精品中国| 日韩欧美在线二视频| 他把我摸到了高潮在线观看| 一本精品99久久精品77| 国产精品av视频在线免费观看| 亚洲成a人片在线一区二区| 成人国产一区最新在线观看| 国产日本99.免费观看| 无遮挡黄片免费观看| 色综合欧美亚洲国产小说| 欧美日本亚洲视频在线播放| 国产中年淑女户外野战色| 少妇人妻一区二区三区视频| 国产v大片淫在线免费观看| 亚洲av成人不卡在线观看播放网| 能在线免费观看的黄片| 国产男靠女视频免费网站| 国产精品亚洲av一区麻豆| 最好的美女福利视频网| 性欧美人与动物交配| avwww免费|