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

    基于CFD的變道超車過(guò)程氣動(dòng)特性分析

    2020-08-17 12:23:12唐洪濤龍世桀楊峻程
    汽車工程學(xué)報(bào) 2020年4期
    關(guān)鍵詞:極值側(cè)向小車

    唐洪濤,龍世桀,楊峻程,王 旭

    (天津科技大學(xué) 機(jī)械工程學(xué)院,天津 300222)

    汽車高速行駛的操作穩(wěn)定性、安全性和舒適性在很大程度上取決于汽車的氣動(dòng)特性[1]。汽車在路面上行駛的工況越來(lái)越復(fù)雜,尤其是汽車在高速超車時(shí),車身周圍流場(chǎng)不斷急劇變化,這不僅嚴(yán)重影響汽車的氣動(dòng)特性,而且會(huì)使車輛發(fā)生偏移,影響車輛行駛時(shí)的穩(wěn)定性和安全性[2]。

    2000年,AZIM等應(yīng)用實(shí)際汽車1/60的模型基于煙霧可視化方法,對(duì)超車過(guò)程的湍流特性進(jìn)行了試驗(yàn)[3],結(jié)果表明,干擾車周圍的流動(dòng)結(jié)構(gòu)根據(jù)雷諾數(shù)的改變而發(fā)生顯著變化,但是模型比例與真實(shí)汽車相差較大,而且可視化方法提供的定量信息很少。2005年,NOGER 等放大了模型尺寸,采用動(dòng)態(tài)試驗(yàn)對(duì)兩個(gè)1/5尺寸的模型進(jìn)行了超車過(guò)程的瞬態(tài)模擬[4],研究表明,縱向間距、橫向間距和相對(duì)速度都會(huì)對(duì)車輛的氣動(dòng)力和力矩產(chǎn)生影響,但是模型較為簡(jiǎn)單,精確程度低。隨后在2008年,CORIN等對(duì)超車過(guò)程進(jìn)行了數(shù)值模擬[5],發(fā)現(xiàn)了相對(duì)速度會(huì)對(duì)超車過(guò)程中車輛周圍的壓力場(chǎng)產(chǎn)生影響,從而使車輛的受力狀態(tài)發(fā)生改變。2013年,BRUZELIUS等基于RANS方法進(jìn)行了數(shù)值模擬[6],分析了超車過(guò)程引起的阻力、側(cè)向力和橫擺力矩系數(shù)的變化規(guī)律。2014年,HOWELL等通過(guò)試驗(yàn)的方法研究了小型汽車超越大型卡車時(shí)的空氣動(dòng)力學(xué)特性[7],得到了試驗(yàn)數(shù)據(jù),并討論了試驗(yàn)過(guò)程的局限性。ALHOMOUD等應(yīng)用SolidWorks軟件設(shè)計(jì)了小型汽車和大型卡車模型并進(jìn)行了模擬仿真[8],探索超車過(guò)程中4個(gè)不同位置上阻力系數(shù)對(duì)兩車的影響,但是該模擬并不是連續(xù)的超車過(guò)程,與真實(shí)的超車過(guò)程有差異。上述研究主要針對(duì)直道超車,對(duì)實(shí)際生活中的直道超車現(xiàn)象具有一定的指導(dǎo)意義。但是,實(shí)際生活中更多情況是變道超車,很少有人將變道超車過(guò)程中兩車的流場(chǎng)變化以及流場(chǎng)中渦的變化與車身受力相結(jié)合進(jìn)行分析。

    綜上所述,本文基于前人的研究,采用動(dòng)網(wǎng)格技術(shù)[9]和SIMPLE算法[10]對(duì)相同車型不同速度、相同車型不同間距、不同車型的變道超車過(guò)程進(jìn)行了數(shù)值模擬,得到了相應(yīng)的壓力云圖、速度矢量圖、流線圖以及側(cè)向力等結(jié)果,并對(duì)側(cè)向力數(shù)據(jù)進(jìn)行分析處理,得出變道超車過(guò)程不同位置處兩車周圍流場(chǎng)變化、車身周圍束縛渦與流場(chǎng)中自由渦的變化過(guò)程,從而解釋車身受力狀態(tài)發(fā)生改變的原因,為汽車行駛提供理論建議。

    1 基本控制方程

    本文選用的汽車行駛速度不超過(guò)100 km/h,馬赫數(shù)約為0.081 7,小于0.3,因此可以認(rèn)為此時(shí)氣體是不可壓縮的[11]。

    連續(xù)性方程:

    x方向上的動(dòng)量守恒方程:

    y方向上的動(dòng)量守恒方程:

    z方向上的動(dòng)量守恒方程:

    式(2)~(4)中:ρ為密度;t為時(shí)間;u、v、w分別表示速度在x、y、z方向上的速度分量;Fx、Fy、Fz分別表示體積力在x、y、z方向上的分量。

    標(biāo)準(zhǔn)k-ε模型的湍動(dòng)能及耗散率運(yùn)輸方程為:

    式(5)~(6)中:μi為湍流粘度;Gk是由平均速度梯度引起的湍動(dòng)能產(chǎn)生項(xiàng);Gb是由浮力影響引起的湍動(dòng)能產(chǎn)生項(xiàng);YM為可壓縮湍流脈動(dòng)膨脹對(duì)總的耗散率的影響;C1ε、C1ε、C3ε為經(jīng)驗(yàn)常數(shù);σk、σε分別為湍動(dòng)能k和湍動(dòng)耗散率ε對(duì)應(yīng)的普朗特?cái)?shù),F(xiàn)LUENT中默認(rèn)值為σk=1.0、σε=1.3;Sk、Sε表示用戶自定義源項(xiàng)。

    2 數(shù)值模擬

    2.1 汽車模型的選取和計(jì)算域的設(shè)定

    汽車模型選用某款家用轎車,其基本尺寸為:長(zhǎng)4 866 mm,寬1 832 mm,高1 464 mm。按長(zhǎng)寬的比例簡(jiǎn)化為平面矩形,矩形的長(zhǎng)用L表示,矩形的寬用W表示,車頭與車尾的棱角作光滑處理,兩車從右向左行駛,主超車命名為CarA,被超車命名為CarB,兩車初始位置距離計(jì)算域右邊界3L,兩車齊頭并進(jìn)。為保證變道超車時(shí)兩車有足夠的行駛距離,且不受兩邊壁面的影響,計(jì)算域尺寸設(shè)計(jì)為:長(zhǎng)30L,寬15W。如圖1 所示。

    圖1 汽車模型和計(jì)算域

    2.2 網(wǎng)格劃分

    模擬采用三角形非結(jié)構(gòu)網(wǎng)格,模擬二維變道超車過(guò)程,小車車身線的網(wǎng)格尺寸設(shè)為0.1,計(jì)算域網(wǎng)格尺寸設(shè)為0.6。這樣的劃分方法可以使小車附近的區(qū)域網(wǎng)格相對(duì)較為密集,而遠(yuǎn)離小車的區(qū)域網(wǎng)格相對(duì)較為稀疏,從而可以提高小車周圍計(jì)算域的精度,更好地完成模擬。如圖2所示。

    圖2 網(wǎng)格模型

    2.3 邊界條件和初始條件的設(shè)定

    計(jì)算域的右邊界設(shè)置為壓力進(jìn)口,左邊界設(shè)置為壓力出口,兩車車身設(shè)置為wall,求解器選用非定常求解器,計(jì)算采用k-epsilon湍流模型,定義流體材料中選擇ideal-gas(理想氣體)。動(dòng)網(wǎng)格參數(shù)的設(shè)置選擇Smoothing(網(wǎng)格光順)和Remeshing(網(wǎng)格重劃分),Smoothing中的彈性因子設(shè)置為0.05[12]。動(dòng)網(wǎng)格中需要對(duì)局部網(wǎng)格進(jìn)行重新劃分,重新劃分網(wǎng)格時(shí)需要判定網(wǎng)格畸形率和網(wǎng)格尺寸是否符合要求,因此,在Remeshing中選用Must Improve Skewness(改善畸形率)和Size Function(尺寸函數(shù)),從而使局部網(wǎng)格完成更好的重新劃分,提高網(wǎng)格質(zhì)量[12]。為了避免動(dòng)網(wǎng)格計(jì)算時(shí)出現(xiàn)負(fù)體積的問(wèn)題,網(wǎng)格的最小體積應(yīng)該大于小車的最快速度與計(jì)算步長(zhǎng)的乘積。

    2.4 變道超車方案設(shè)計(jì)

    記CarA的中心與CarB的中心沿x軸正方向的距離為X,沿y軸正方向的距離為Y,用X/L和Y/W表示兩車的位置。變道超車方案見(jiàn)表1。

    表1 變道超車方案

    3 模擬結(jié)果分析

    3.1 變道超車壓力云圖分析

    選擇Case2的壓力云圖進(jìn)行分析。選取超車過(guò)程的5個(gè)位置,分別是X/L=0.5、Y/W=1.5,X/L=1、Y/W=1.25,X/L=1.5、Y/W=1,X/L=2、Y/W=0.75,X/L=3.5、Y/W=0。單車行駛時(shí),在不考慮側(cè)風(fēng)的影響下,小車所受的側(cè)向力和旋轉(zhuǎn)力矩幾乎為0,但復(fù)數(shù)車輛變道超車過(guò)程中,兩車會(huì)產(chǎn)生相互影響,圖例單位為Pa。

    圖3 變道超車壓力云圖

    如圖3a所示,CarA在x方向上超越CarB半個(gè)車身長(zhǎng)時(shí),CarA前方的正壓區(qū)與CarB前方的正壓區(qū)相互作用,產(chǎn)生不均衡的側(cè)向力。兩車之間的區(qū)域形成負(fù)壓區(qū),產(chǎn)生相互吸引的趨勢(shì)。CarA兩側(cè)的氣壓在短時(shí)間內(nèi)發(fā)生劇烈變化,局部受力不均衡,車頭內(nèi)側(cè)壓力高于外側(cè)壓力,車尾內(nèi)側(cè)壓力低于外側(cè)壓力,產(chǎn)生逆時(shí)針的趨勢(shì),車身發(fā)生振動(dòng)。由于局部受力不均衡,CarB車身內(nèi)側(cè)壓力低于外側(cè)壓力,產(chǎn)生逆時(shí)針旋轉(zhuǎn)的趨勢(shì)。此時(shí),CarB的車頭有很大幾率撞擊CarA的車身,屬于變道超車過(guò)程中的危險(xiǎn)位置。

    如圖3b所示,CarA在x方向上完成對(duì)CarB的超越時(shí),CarA內(nèi)側(cè)氣壓受到CarB車頭處高壓的影響,相較外側(cè)呈現(xiàn)高壓狀態(tài),兩側(cè)產(chǎn)生壓力差,CarA車身受到內(nèi)側(cè)指向外側(cè)的側(cè)向力,車尾處內(nèi)側(cè)區(qū)域呈現(xiàn)負(fù)壓區(qū),受到外側(cè)指向內(nèi)側(cè)的側(cè)向力,車身產(chǎn)生逆時(shí)針旋轉(zhuǎn)趨勢(shì),車身振動(dòng)加劇。CarB內(nèi)側(cè)氣壓受到CarA尾部氣壓的影響相較外側(cè)呈現(xiàn)高壓,CarB車身受到內(nèi)側(cè)指向外側(cè)的側(cè)向力,車頭處受到外側(cè)指向內(nèi)側(cè)的側(cè)向力,車身同樣產(chǎn)生逆時(shí)針旋轉(zhuǎn)趨勢(shì),此時(shí)兩車容易出現(xiàn)剮蹭現(xiàn)象,因此,兩車應(yīng)該在橫向上保持安全距離。

    如圖3c所示,CarA在x方向上超出CarB半個(gè)車身長(zhǎng)時(shí),CarA車身內(nèi)側(cè)仍然受到CarB車頭處高壓的影響,兩側(cè)產(chǎn)生壓力差,產(chǎn)生側(cè)向力,隨著兩車x方向間距的增大,側(cè)向力逐漸減小。同理,CarB的內(nèi)側(cè)氣壓受到CarA尾部氣流的影響減小,產(chǎn)生的側(cè)向力也隨之減小。

    如圖3d和圖3e所示,CarA超越CarB變道至被超車車道的過(guò)程中,隨著兩車x方向間距的增大,兩車所受到的側(cè)向力逐漸減小。CarA的尾流對(duì)CarB的內(nèi)側(cè)影響逐漸減小。

    3.2 速度矢量圖分析

    在變道超車過(guò)程中,運(yùn)動(dòng)的小車周圍形成不均衡流動(dòng)甚至于偏置的流場(chǎng),兩個(gè)流場(chǎng)區(qū)域相互干涉,不斷變化,從而引起小車氣動(dòng)力發(fā)生變化。圖例單位m/s。

    圖4 變道超車速度矢量圖

    圖4b為圖4a的局部放大圖,CarA在x方向上超越CarB半個(gè)車身長(zhǎng)時(shí),兩車之間的區(qū)域氣流相互干涉,該區(qū)域氣流明顯加速,形成低壓區(qū)。CarA車身內(nèi)側(cè)的氣流與CarB車頭內(nèi)側(cè)的氣流方向相反,相互作用產(chǎn)生排斥力,CarA車頭處受到內(nèi)側(cè)指向外側(cè)的力,車尾處受到低壓區(qū)的影響,產(chǎn)生外側(cè)指向內(nèi)側(cè)的力,產(chǎn)生逆時(shí)針旋轉(zhuǎn)的趨勢(shì),CarA車身發(fā)生振動(dòng),對(duì)駕駛員的操作有很大的影響。CarB車頭處受到由外側(cè)指向內(nèi)側(cè)的力,產(chǎn)生逆時(shí)針旋轉(zhuǎn)的趨勢(shì)。

    圖4d為圖4c的局部放大圖,CarA在x方向上完成對(duì)CarB的超越,在X/L=0.5到X/L=1的過(guò)程中,CarA所受到的側(cè)向力逐漸增大至最大,CarA內(nèi)側(cè)氣流干擾加劇,車身振動(dòng)加劇。CarB內(nèi)側(cè)氣流受到CarA車尾處氣流影響,氣流方向、大小都發(fā)生了變化,產(chǎn)生渦,由于渦的渦量和位置發(fā)生變化,導(dǎo)致CarB車身受力變化明顯,對(duì)CarB的行駛穩(wěn)定性造成嚴(yán)重的影響。

    如圖4e所示,CarA在x方向上超出CarB半個(gè)車身長(zhǎng)時(shí),隨著兩車x方向距離的增大,兩車的相互干涉減小,CarA受到的側(cè)向力逐漸減小。CarA尾部氣流與CarB車頭處氣流相互干涉,氣流加速,在該區(qū)域形成低壓區(qū)。CarB受到的側(cè)向力,同時(shí)隨著兩車x方向的距離增大而減小。

    如圖4f和圖4g所示,CarA超越CarB變道至被超車車道的過(guò)程中,隨著兩車x方向距離的增大,兩車流場(chǎng)的相互干涉減小,兩車基本恢復(fù)單車行駛狀態(tài)。

    3.3 流線圖分析

    對(duì)Case2的流線圖進(jìn)行瞬態(tài)流線圖采集,時(shí)間步長(zhǎng)為0.000 1 s,每100個(gè)時(shí)間步長(zhǎng)采集一張圖片。通過(guò)對(duì)比圖5所示X/L=0.1到X/L=1.5的瞬態(tài)流線圖分析發(fā)現(xiàn),從X/L=0.6開(kāi)始在兩車壁面的內(nèi)外兩側(cè)分別出現(xiàn)一對(duì)明顯束縛渦,隨著兩車相對(duì)位置的變化,束縛渦的渦量發(fā)生變化,使兩車周圍的流場(chǎng)比較紊亂;在CarA的尾流中存在自由渦,這個(gè)渦產(chǎn)生、變大、運(yùn)動(dòng)、消失的過(guò)程,引起兩車之間流場(chǎng)相互干涉。

    圖5 變道超車流線圖

    在X/L=0.1到X/L=0.5的過(guò)程中,兩車內(nèi)側(cè)流線相互作用,氣流加速,此處形成低壓區(qū),兩車產(chǎn)生相互吸引的趨勢(shì)。隨著兩車y方向距離減小,兩車相互吸引的趨勢(shì)逐漸增強(qiáng)。在此過(guò)程中,兩車相互吸引的趨勢(shì)占主導(dǎo)地位。兩車所受到的側(cè)向力逐漸增大,在X/L=0.5時(shí),CarB所受到的側(cè)向力增至最大。

    在車身周圍循環(huán)流動(dòng)通常發(fā)生在車身表面的旋渦稱為束縛渦;在車身上形成的循環(huán)流動(dòng)向外側(cè)擴(kuò)散,并最終消失在大氣中的旋渦稱為自由渦。

    隨著兩車縱向間距的增大,在X/L=0.6時(shí),存在4個(gè)比較明顯的束縛渦,分別命名為1、2、3、4渦,其中1渦、2渦在貼近CarA壁面的位置,1渦為逆時(shí)針?lè)较蛐D(zhuǎn),渦尺寸較小,2渦為順時(shí)針?lè)较蛐D(zhuǎn),渦尺寸較大,3渦、4渦在貼近CarB壁面的位置,3渦為逆時(shí)針?lè)较蛐D(zhuǎn),渦尺寸較大,4渦為順時(shí)針?lè)较蛐D(zhuǎn),渦尺寸較小。1渦、4渦在CarA、CarB的外側(cè)受到的影響較小,發(fā)生的變化也較小。

    束縛渦(2渦)的產(chǎn)生、變化消耗流場(chǎng)中的能量,引起周圍壓強(qiáng)的改變,進(jìn)而影響CarA的受力狀態(tài),同時(shí)兩車縱向間距增大,低壓區(qū)范圍減小,相互吸引的趨勢(shì)減弱,因此,CarA受到由內(nèi)側(cè)指向外側(cè)的側(cè)向力增大。

    在X/L=0.7時(shí),2渦渦尺寸增大,3渦渦尺寸增大,5渦是從由2渦脫離車身表面進(jìn)入大氣內(nèi)形成的自由渦,5渦呈順時(shí)針?lè)较蛐D(zhuǎn),渦尺寸較小。

    自由渦的產(chǎn)生消耗流場(chǎng)中的能量,產(chǎn)生自由渦的位置,壓強(qiáng)會(huì)發(fā)生急劇變化,同時(shí)兩車相互吸引的趨勢(shì)減弱,CarB受到由外側(cè)指向內(nèi)側(cè)的側(cè)向力減小。

    在X/L=0.8時(shí),2渦向內(nèi)側(cè)運(yùn)動(dòng),3渦渦尺寸減小,5渦向后運(yùn)動(dòng),逐漸脫離CarA的尾流流場(chǎng),渦尺寸增大,呈逆時(shí)針?lè)较蛐D(zhuǎn)。

    在X/L=0.9時(shí),2渦向后運(yùn)動(dòng),3渦渦尺寸增大,5渦持續(xù)向后運(yùn)動(dòng),渦尺寸持續(xù)減小,隨著5渦從CarA的尾流場(chǎng)中脫離出來(lái),在CarA的尾流場(chǎng)中產(chǎn)生6渦,6渦渦尺寸較小,呈順時(shí)針?lè)较蛐D(zhuǎn),與5渦旋轉(zhuǎn)方向相反。

    在X/L=1.0時(shí),2渦渦尺寸持續(xù)增大,3渦向前運(yùn)動(dòng),5渦持續(xù)向后運(yùn)動(dòng)從捕捉平面消失,6渦向后運(yùn)動(dòng),渦量增大,呈逆時(shí)針?lè)较蛐D(zhuǎn),6渦展現(xiàn)出5渦的特征。

    束縛渦渦尺寸、位置的變化過(guò)程,引起流場(chǎng)中壓強(qiáng)發(fā)生改變,束縛渦(2渦)渦尺寸持續(xù)增加,位置不斷發(fā)生變化,對(duì)CarA受力狀態(tài)的影響增大,因此,CarA受到的側(cè)向力不斷增加,在X/L=1.0時(shí),側(cè)向力增大至最大。

    自由渦產(chǎn)生、變大、運(yùn)動(dòng)、消失的過(guò)程消耗流場(chǎng)的能量,隨著5渦、6渦渦尺寸增大和位置向后運(yùn)動(dòng),對(duì)CarB受力狀態(tài)的影響逐漸增大,在圖6b側(cè)向力變化曲線中的表現(xiàn)為CarB所受側(cè)向力的斜率增大。

    圖6 不同超車速度側(cè)向力變化曲線

    在X/L=1.1時(shí),2渦渦尺寸減小,3渦渦尺寸減小,6渦與5渦的運(yùn)動(dòng)狀態(tài)相同,持續(xù)向后運(yùn)動(dòng)從捕捉平面消失,CarA尾流場(chǎng)中7渦逐漸產(chǎn)生,7渦渦尺寸較小,呈順時(shí)針?lè)较蛐D(zhuǎn)。在X/L=1.2時(shí),2渦渦尺寸增大,3渦向后運(yùn)動(dòng),7渦向后運(yùn)動(dòng),渦尺寸明顯增大,呈順時(shí)針?lè)较蛐D(zhuǎn),7渦呈現(xiàn)出與5渦、6渦相同的特征。在X/L=1.3時(shí),3渦向前運(yùn)動(dòng),渦尺寸增大,7渦持續(xù)向后運(yùn)動(dòng),逐漸從CarA的尾流中脫離出來(lái),形成一個(gè)獨(dú)立的呈逆時(shí)針?lè)较蛐D(zhuǎn)的自由渦。在X/L=1.4時(shí),3渦消失,7渦持續(xù)向后運(yùn)動(dòng),渦尺寸減小,呈逆時(shí)針?lè)较蛐D(zhuǎn)。在X/L=1.5時(shí),7渦持續(xù)向后運(yùn)動(dòng)從捕捉平面消失。隨著兩車x方向間距增大,CarA的尾流場(chǎng)中不再產(chǎn)生自由渦。

    束縛渦(2渦)在X/L=1.1時(shí)渦尺寸開(kāi)始減弱,CarA受到的側(cè)向力逐漸減弱。隨著兩車縱向間距的增大,兩車之間的相互干涉減小,側(cè)向力減小,兩車行駛狀態(tài)逐漸趨于平穩(wěn)。

    自由渦(7渦)在X/L=1.1時(shí)產(chǎn)生,X/L=1.2時(shí)變大,X/L=1.3時(shí)向后運(yùn)動(dòng)的過(guò)程中,對(duì)CarB的影響持續(xù)增大。隨著自由渦向后運(yùn)動(dòng)至車尾處后消失,自由渦對(duì)CarB受力狀態(tài)的影響開(kāi)始減弱。在圖6b側(cè)向力變化圖中的表現(xiàn)為:在X/L=1.3附近側(cè)向力趨勢(shì)發(fā)生改變。

    通過(guò)以上對(duì)兩車之間流場(chǎng)流線的瞬態(tài)仿真分析,在變道超車過(guò)程中,高速氣流在兩車內(nèi)外兩側(cè)壁面位置和兩車之間區(qū)域變化較為復(fù)雜,隨著兩車相對(duì)位置的變化,兩對(duì)束縛渦的渦尺寸大小持續(xù)變化,旋轉(zhuǎn)方向相反。變道超車過(guò)程中,CarA尾流場(chǎng)中出現(xiàn)自由渦,自由渦產(chǎn)生、變大、運(yùn)動(dòng)、消失的過(guò)程,加劇了兩車周圍流場(chǎng)的相互干涉,因此,研究束縛渦渦量和位置的變化、自由渦從產(chǎn)生到消失的過(guò)程對(duì)兩車受力狀態(tài)的影響具有重要意義。

    3.4 不同車速數(shù)據(jù)分析

    對(duì)比分析Case1~3在X/L=0.5、Y/W=1.5位置的壓力云圖可知,隨著兩車的速度加快,兩車車頭處正壓值依次增加,相互作用依次加強(qiáng),兩車之間的負(fù)壓區(qū)范圍和最大負(fù)壓值依次增大。因此,在負(fù)壓區(qū)的影響下,主超車和被超車都受到側(cè)向力和旋轉(zhuǎn)力矩的作用,隨著兩車速度的增加,側(cè)向力和旋轉(zhuǎn)力矩依次增加。

    3.5 不同車距數(shù)據(jù)分析

    對(duì)比分析Case4、 Case2和Case5在X/L=0.5位置的壓力云圖可以得出,隨著兩車的間距增大,兩車之間的負(fù)壓區(qū)范圍和最大負(fù)壓值依次減小。隨著變道超車的間距變大,主超車和被超車所受到的側(cè)向力和旋轉(zhuǎn)力矩依次降低。

    圖7 不同超車速度壓力云圖

    圖8 不同超車間距壓力云圖

    3.6 側(cè)向力數(shù)據(jù)分析

    3.6.1 不同車速的變道超車數(shù)據(jù)分析

    本文中側(cè)向力的方向是垂直于小車行駛方向的力。在單車行駛狀態(tài),不考慮側(cè)風(fēng)影響的條件下,所受到的側(cè)向力基本為0。在變道超車過(guò)程中,小車所受到的側(cè)向力有明顯的變化,在超車過(guò)程中相對(duì)位置不同,所受到的側(cè)向力也不同。當(dāng)側(cè)向力發(fā)生劇烈變化時(shí),小車的行駛穩(wěn)定性會(huì)受到嚴(yán)重的影響。

    主超車CarA在Case1、Case2、Case3變道車的超車過(guò)程中所受到的側(cè)向力如圖6a所示,主超車CarA所受到的側(cè)向力先增大再減小后趨于穩(wěn)定,在X/L=1時(shí)側(cè)向力達(dá)到最大。 Case2與Case1、Case3對(duì)比發(fā)現(xiàn):控制兩車間距不變,兩車速度越快,側(cè)向力的變化幅度越大。主超車CarA在Case1的極值比Case2的極值增大了56%, Case1的極值比Case3的極值增大了156%,Case2的極值比Case3的極值增大了64%。

    被超車CarB在Case1、Case2、Case3變道車超車過(guò)程中所受到的側(cè)向力如圖 6b所示,被超車CarB所受到的側(cè)向力先增大再減小后趨于穩(wěn)定,在X/L=0.5時(shí)側(cè)向力達(dá)到最大;與主超車CarA相似,隨著車速的增加,側(cè)向力的變化幅度增大。Case2與Case1、Case3對(duì)比發(fā)現(xiàn):被超車CarB在Case1的極值比Case2的極值增大了68%,Case1的極值比Case3的極值增大了199%,Case2的極值比Case3的極值增大了78%。根據(jù)數(shù)據(jù)可知,改變車速變道超車時(shí),被超車受到的影響大于主超車。

    3.6.2 不同間距的變道超車數(shù)據(jù)分析

    主超車CarA在 Case2、Case4、Case5變道車超車過(guò)程中所受到的側(cè)向力如圖9a所示,控制兩車車速相同,改變兩車間距對(duì)側(cè)向力的影響。隨著兩車間距的增加,主超車受到的側(cè)向力變化幅度降低。主超車CarA在Case4的極值比Case2的極值增大了42%,Case4的極值比Case5的極值增大了74%,Case2的極值比Case5的極值增大了23%。

    被超車CarB在Case2、Case4、Case5變道車超車過(guò)程中所受到的側(cè)向力如圖9b所示,與主超車相同,隨著兩車間距的增大,被超車受到的側(cè)向力減小。被超車CarB在Case4的極值比Case2的極值增大了18%,Case4的極值比Case5的極值增大了41%,Case2的極值比Case5的極值增大了18%。根據(jù)數(shù)據(jù)可知,改變間距變道超車時(shí),主超車受到的影響大于被超車。

    圖9 不同超車間距側(cè)向力變化曲線

    3.7 不同車型變道超車的分析

    3.7.1 小車超越大車

    以上是對(duì)相似車型變道超車的研究,下面對(duì)不同車型變道超車進(jìn)行探索,對(duì)小車超越大車的壓力圖和側(cè)向力的變化進(jìn)行分析。該方案命名為Case6,小車命名CarA,大車命名CarB。選取X/L=0.5和X/L=1的位置進(jìn)行分析。

    如圖10a所示:不同車型變道超車相較于相同車型變道超車,兩車車頭處正壓區(qū)的相互作用明顯增大,兩車之間的負(fù)壓區(qū)范圍和最大負(fù)壓值有所增大,說(shuō)明兩車相互吸引的趨勢(shì)增強(qiáng),側(cè)向力增大。CarB的尺寸和重量較CarA相差較大,側(cè)向力對(duì)CarA的影響更為明顯,逆時(shí)針旋轉(zhuǎn)趨勢(shì)增強(qiáng),CarA的操作穩(wěn)定性受到嚴(yán)重影響,車身出現(xiàn)振動(dòng),有很大可能與CarB相撞。

    如圖10b所示:不同車型變道超車時(shí),兩車正壓區(qū)作用范圍明顯增強(qiáng),CarA內(nèi)側(cè)壓力升高,兩側(cè)壓力差增大,CarA受到的側(cè)向力增加至最大,逆時(shí)針旋轉(zhuǎn)趨勢(shì)增強(qiáng),CarA的操作穩(wěn)定性受到影響,車身的振動(dòng)現(xiàn)象加劇,有一定幾率被CarB“吸引”過(guò)來(lái),發(fā)生車禍。

    圖10 不同車型變道超車壓力云圖(小車超越大車)

    對(duì)比不同車型變道超車的側(cè)向力發(fā)現(xiàn):主超車CarA在變道超車的過(guò)程中,所受到側(cè)向力的變化幅度明顯增加。主超車CarA在Case6的極值比Case2的極值增大了60%。被超車CarB受到的側(cè)向力也有所增加,被超車CarB在Case6的極值比Case2的極值增大了27%,主超車CarA所受到側(cè)向力的變化幅度更大,對(duì)小車的操作穩(wěn)定性影響更大。

    圖11 不同車型變道超車側(cè)向力變化曲線(小車超越大車)

    3.7.2 大車超越小車

    大車變道超越小車的情況,同樣對(duì)壓力圖和側(cè)向力的變化進(jìn)行分析,該方案命名為Case7,大車命名為CarC,小車命名為CarD。同樣選取X/L=0.5和X/L=1的位置進(jìn)行分析。

    如圖12a所示:不同車型變道超車時(shí),由于CarC的尺寸和質(zhì)量較大,車頭處的正壓區(qū)范圍和正壓值較相同車型變道超車明顯增大,兩車之間負(fù)壓區(qū)范圍和負(fù)壓最大值有所增加。此時(shí)CarD所受的側(cè)向力達(dá)到最大值,對(duì)CarD的操作穩(wěn)定性有很大的影響,有很高機(jī)率與CarC相撞。因此,應(yīng)控制兩車保持安全的距離。

    如圖12b所示:不同車型變道超車相較于相同車型變道超車,兩車的車頭處正壓區(qū)的相互作用明顯增強(qiáng),CarC的內(nèi)側(cè)壓力增大,導(dǎo)致兩側(cè)壓差增大,此時(shí)CarC受到的側(cè)向力也增加至最大,車身振動(dòng)加劇,對(duì)CarC的操作穩(wěn)定性造成影響。因此,CarD應(yīng)增大與CarC的間距,降低行車速度,避免發(fā)生事故。

    圖12 不同車型變道超車壓力云圖(大車超越小車)

    當(dāng)大車超越小車時(shí),兩車變道超車過(guò)程中,各個(gè)位置所受到的側(cè)向力明顯增大。主超車CarC在Case7的極值比Case2的極值增大了94%,被超車CarD在Case7的極值比Case2的極值增大了47%,主超車CarC所受到的側(cè)向力變化幅度增大,對(duì)大車的操作穩(wěn)定性影響增大。由于變道超車過(guò)程中車型的不同,兩車所受的側(cè)向力變化幅度增大。

    圖13 不同車型變道超車側(cè)向力變化曲線(大車超越小車)

    4 結(jié)論

    (1)變道超車過(guò)程中壓力、速度、側(cè)向力隨著兩車的相對(duì)位置的變化而變化,兩車所受側(cè)向力都會(huì)先增大后減小,然后趨于穩(wěn)定,主超車在X/L=1時(shí),側(cè)向力達(dá)到最大值,被超車在X/L=0.5時(shí),側(cè)向力達(dá)到最大。

    (2)變道超車過(guò)程中,汽車的行駛穩(wěn)定性與側(cè)向力的大小有密切的聯(lián)系,在X/L=0.5到X/L=1的位置的過(guò)程中,側(cè)向力較大,此時(shí)屬于變道超車的危險(xiǎn)位置。

    (3)變道超車過(guò)程相較于直道超車過(guò)程流場(chǎng)變化更為復(fù)雜,變道超車過(guò)程中車身兩側(cè)的束縛渦渦量大小持續(xù)變化,旋轉(zhuǎn)方向相反;流場(chǎng)中自由渦經(jīng)歷產(chǎn)生、變大、發(fā)展、消失的過(guò)程,消耗流場(chǎng)中的能量,加劇兩車周圍的流線發(fā)生偏置,從而影響兩車的受力狀態(tài)。

    (4)通過(guò)方案對(duì)比發(fā)現(xiàn):變道超車過(guò)程中,兩車車速增大,側(cè)向力增大,被超車受到的影響大于主超車;兩車間距減小,側(cè)向力增大,主超車受到的影響大于被超車。

    (5)不同車型變道超車時(shí),側(cè)向力的變化會(huì)有所不同,小車超越大車時(shí),小車受到的影響遠(yuǎn)大于大車,嚴(yán)重影響駕駛員的操作穩(wěn)定性,為了避免發(fā)生交通事故,應(yīng)降低超車速度,增加超車間距。

    猜你喜歡
    極值側(cè)向小車
    極值點(diǎn)帶你去“漂移”
    軍航無(wú)人機(jī)與民航航班側(cè)向碰撞風(fēng)險(xiǎn)評(píng)估
    快樂(lè)語(yǔ)文(2020年36期)2021-01-14 01:10:32
    極值點(diǎn)偏移攔路,三法可取
    自制小車來(lái)比賽
    一類“極值點(diǎn)偏移”問(wèn)題的解法與反思
    劉老師想開(kāi)小車
    文苑(2018年22期)2018-11-19 02:54:18
    兩輪自平衡小車的設(shè)計(jì)與實(shí)現(xiàn)
    電子制作(2018年8期)2018-06-26 06:43:02
    彎月薄鏡的側(cè)向支撐
    側(cè)向風(fēng)場(chǎng)中無(wú)人機(jī)的飛行研究
    91成年电影在线观看| 视频区图区小说| 久久久国产成人精品二区 | 日韩成人在线观看一区二区三区| 深夜精品福利| 午夜精品在线福利| 黄片播放在线免费| 国产亚洲精品久久久久久毛片| 免费看a级黄色片| 久久久水蜜桃国产精品网| 亚洲国产欧美一区二区综合| 美女扒开内裤让男人捅视频| 欧美激情久久久久久爽电影 | 久久精品影院6| 一级片免费观看大全| 最新在线观看一区二区三区| 91在线观看av| 波多野结衣高清无吗| av视频免费观看在线观看| 黄频高清免费视频| 女人精品久久久久毛片| 91九色精品人成在线观看| 一区二区三区激情视频| 1024香蕉在线观看| 久久99一区二区三区| 精品一区二区三区四区五区乱码| 看黄色毛片网站| 三上悠亚av全集在线观看| 老鸭窝网址在线观看| 很黄的视频免费| 夜夜躁狠狠躁天天躁| 国产精品98久久久久久宅男小说| 久久狼人影院| 日韩欧美免费精品| 无人区码免费观看不卡| 国产aⅴ精品一区二区三区波| 欧美性长视频在线观看| 最好的美女福利视频网| 男女之事视频高清在线观看| 视频区图区小说| 最近最新中文字幕大全电影3 | e午夜精品久久久久久久| 国产成人av激情在线播放| av免费在线观看网站| 一级a爱片免费观看的视频| 日韩视频一区二区在线观看| 三级毛片av免费| 亚洲精品久久午夜乱码| 三上悠亚av全集在线观看| 深夜精品福利| 午夜精品在线福利| 50天的宝宝边吃奶边哭怎么回事| 久久精品91蜜桃| 在线观看免费视频网站a站| 黄频高清免费视频| √禁漫天堂资源中文www| 国产精品久久久久成人av| 91字幕亚洲| 国产一区二区在线av高清观看| 国产精品偷伦视频观看了| 伦理电影免费视频| 动漫黄色视频在线观看| 国产精品九九99| 欧美乱妇无乱码| svipshipincom国产片| 不卡av一区二区三区| 国产三级在线视频| 精品国产乱码久久久久久男人| 国产成人av教育| 亚洲色图综合在线观看| 自线自在国产av| 视频区欧美日本亚洲| 国产av又大| 精品第一国产精品| 每晚都被弄得嗷嗷叫到高潮| 国产免费现黄频在线看| 精品久久久久久电影网| 国产亚洲精品一区二区www| 午夜精品国产一区二区电影| 亚洲成av片中文字幕在线观看| 人成视频在线观看免费观看| 少妇的丰满在线观看| 一本综合久久免费| 天天躁狠狠躁夜夜躁狠狠躁| 99精国产麻豆久久婷婷| 十分钟在线观看高清视频www| 久久国产亚洲av麻豆专区| 久久久久久免费高清国产稀缺| 精品乱码久久久久久99久播| 99香蕉大伊视频| 亚洲情色 制服丝袜| 侵犯人妻中文字幕一二三四区| 亚洲国产看品久久| 免费看十八禁软件| 人妻丰满熟妇av一区二区三区| 久热爱精品视频在线9| 久热这里只有精品99| av网站免费在线观看视频| 精品久久久久久成人av| 少妇 在线观看| 日韩精品免费视频一区二区三区| 色精品久久人妻99蜜桃| 色尼玛亚洲综合影院| 亚洲欧美精品综合一区二区三区| 亚洲视频免费观看视频| 中文欧美无线码| 国产精品一区二区免费欧美| 欧美一区二区精品小视频在线| 五月开心婷婷网| 99久久国产精品久久久| 女人爽到高潮嗷嗷叫在线视频| 精品人妻1区二区| 亚洲性夜色夜夜综合| 午夜老司机福利片| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产毛片av蜜桃av| 亚洲aⅴ乱码一区二区在线播放 | videosex国产| 熟女少妇亚洲综合色aaa.| 9热在线视频观看99| 在线观看一区二区三区| 又大又爽又粗| 脱女人内裤的视频| 最好的美女福利视频网| 69av精品久久久久久| 亚洲国产中文字幕在线视频| 亚洲欧美激情在线| 十八禁网站免费在线| 如日韩欧美国产精品一区二区三区| 黄色a级毛片大全视频| 久久久久久免费高清国产稀缺| 91精品三级在线观看| 国产色视频综合| 两个人看的免费小视频| www.自偷自拍.com| 亚洲欧美一区二区三区黑人| 老司机午夜福利在线观看视频| 高清欧美精品videossex| 精品久久久久久,| 91av网站免费观看| 国产人伦9x9x在线观看| 日韩精品青青久久久久久| 久久久精品欧美日韩精品| 亚洲一区高清亚洲精品| 一级黄色大片毛片| 久久精品人人爽人人爽视色| 欧美中文综合在线视频| 一进一出抽搐gif免费好疼 | 丰满饥渴人妻一区二区三| 国产欧美日韩一区二区三| 久久精品亚洲av国产电影网| 亚洲自偷自拍图片 自拍| 国产在线观看jvid| 成年人黄色毛片网站| 在线观看免费高清a一片| 亚洲美女黄片视频| 日韩高清综合在线| 欧美乱码精品一区二区三区| 国产在线观看jvid| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美乱色亚洲激情| 国产精品免费一区二区三区在线| 狂野欧美激情性xxxx| 亚洲精品成人av观看孕妇| 老汉色av国产亚洲站长工具| e午夜精品久久久久久久| 国产精品国产高清国产av| 在线观看免费视频网站a站| tocl精华| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品国产精品久久久不卡| 很黄的视频免费| 久久久久国产精品人妻aⅴ院| 午夜福利一区二区在线看| 大香蕉久久成人网| 香蕉久久夜色| 国产精品二区激情视频| e午夜精品久久久久久久| 国产乱人伦免费视频| 欧美黑人欧美精品刺激| 99久久精品国产亚洲精品| 国产免费男女视频| 日本精品一区二区三区蜜桃| 99久久精品国产亚洲精品| 色婷婷久久久亚洲欧美| 久久精品aⅴ一区二区三区四区| 99riav亚洲国产免费| 精品一区二区三区av网在线观看| 99久久人妻综合| 18禁美女被吸乳视频| 亚洲精品国产一区二区精华液| 亚洲欧美日韩另类电影网站| 19禁男女啪啪无遮挡网站| 国产区一区二久久| 激情视频va一区二区三区| 91精品三级在线观看| 欧美黄色淫秽网站| 久热爱精品视频在线9| 国产精品九九99| 可以免费在线观看a视频的电影网站| 9热在线视频观看99| 一二三四社区在线视频社区8| 久久香蕉激情| 老熟妇乱子伦视频在线观看| 亚洲精品粉嫩美女一区| 校园春色视频在线观看| 欧美成人午夜精品| 午夜影院日韩av| 亚洲九九香蕉| 黄色视频,在线免费观看| 无遮挡黄片免费观看| 亚洲精品粉嫩美女一区| 日韩大码丰满熟妇| 免费av中文字幕在线| 亚洲午夜精品一区,二区,三区| 免费在线观看影片大全网站| 国产亚洲精品第一综合不卡| 国产一区二区三区综合在线观看| 亚洲黑人精品在线| 精品久久久久久电影网| 看片在线看免费视频| 少妇粗大呻吟视频| 18禁国产床啪视频网站| 交换朋友夫妻互换小说| 亚洲欧美精品综合久久99| 黑人操中国人逼视频| 亚洲精品中文字幕在线视频| 他把我摸到了高潮在线观看| 久久久久久久久久久久大奶| 精品欧美一区二区三区在线| 在线观看舔阴道视频| 久久久久亚洲av毛片大全| 99在线人妻在线中文字幕| tocl精华| 国产成人av激情在线播放| 国产精品影院久久| 国产精华一区二区三区| 亚洲自偷自拍图片 自拍| 亚洲成人国产一区在线观看| 欧美中文日本在线观看视频| 亚洲精品国产区一区二| 日本免费a在线| 丰满迷人的少妇在线观看| 91字幕亚洲| av网站免费在线观看视频| 国产1区2区3区精品| 天天影视国产精品| 久99久视频精品免费| 色老头精品视频在线观看| 老汉色av国产亚洲站长工具| 搡老乐熟女国产| 美女大奶头视频| 国产成人影院久久av| 欧美 亚洲 国产 日韩一| 一二三四社区在线视频社区8| 亚洲精品国产一区二区精华液| 国产精品综合久久久久久久免费 | 91麻豆精品激情在线观看国产 | 男人舔女人的私密视频| 亚洲狠狠婷婷综合久久图片| 久久午夜亚洲精品久久| 精品国产美女av久久久久小说| 国产伦人伦偷精品视频| 国产区一区二久久| 欧美乱色亚洲激情| 国产又爽黄色视频| 亚洲欧美一区二区三区黑人| 亚洲第一欧美日韩一区二区三区| 又黄又粗又硬又大视频| 村上凉子中文字幕在线| 香蕉丝袜av| 国产成人欧美在线观看| 一本大道久久a久久精品| 国产野战对白在线观看| 黑丝袜美女国产一区| 日本撒尿小便嘘嘘汇集6| 国产伦一二天堂av在线观看| 最好的美女福利视频网| 大型av网站在线播放| 99久久综合精品五月天人人| av免费在线观看网站| a级片在线免费高清观看视频| 欧美精品一区二区免费开放| 99精品在免费线老司机午夜| 精品电影一区二区在线| 久久久久久久午夜电影 | 激情在线观看视频在线高清| 国产不卡一卡二| 18禁裸乳无遮挡免费网站照片 | 日韩欧美国产一区二区入口| 一二三四社区在线视频社区8| 免费在线观看完整版高清| 男人的好看免费观看在线视频 | 欧美精品啪啪一区二区三区| 亚洲精品久久成人aⅴ小说| 国产精品久久视频播放| 亚洲欧美日韩高清在线视频| 妹子高潮喷水视频| 亚洲精品国产精品久久久不卡| 亚洲成a人片在线一区二区| 亚洲欧美一区二区三区久久| 精品久久久久久,| 欧美乱妇无乱码| 亚洲国产精品合色在线| 久久中文字幕人妻熟女| 美女福利国产在线| 欧美激情极品国产一区二区三区| 50天的宝宝边吃奶边哭怎么回事| a级片在线免费高清观看视频| 亚洲第一欧美日韩一区二区三区| 国产av在哪里看| 久久精品亚洲av国产电影网| 国产精品国产高清国产av| 淫秽高清视频在线观看| 怎么达到女性高潮| av视频免费观看在线观看| 狂野欧美激情性xxxx| 日本免费一区二区三区高清不卡 | 亚洲男人的天堂狠狠| 韩国精品一区二区三区| 一级作爱视频免费观看| 久久香蕉精品热| 亚洲七黄色美女视频| 丁香六月欧美| 亚洲第一青青草原| 级片在线观看| 国产一区在线观看成人免费| 精品一区二区三区四区五区乱码| 久久香蕉精品热| 男女下面插进去视频免费观看| 麻豆一二三区av精品| 国产一卡二卡三卡精品| 亚洲第一欧美日韩一区二区三区| 国产亚洲精品第一综合不卡| 国产精品美女特级片免费视频播放器 | 一本大道久久a久久精品| 天堂√8在线中文| 久久久久精品国产欧美久久久| 亚洲国产毛片av蜜桃av| 国产成人欧美在线观看| 亚洲色图av天堂| 亚洲成人国产一区在线观看| 1024视频免费在线观看| 精品午夜福利视频在线观看一区| 水蜜桃什么品种好| 新久久久久国产一级毛片| 一级a爱片免费观看的视频| 日本精品一区二区三区蜜桃| 国产亚洲欧美在线一区二区| 久久性视频一级片| 一级毛片女人18水好多| 三上悠亚av全集在线观看| 91av网站免费观看| 高清在线国产一区| 国产欧美日韩一区二区三区在线| 精品久久久久久久毛片微露脸| 18禁裸乳无遮挡免费网站照片 | av超薄肉色丝袜交足视频| 久99久视频精品免费| 黄片播放在线免费| 久久久久久人人人人人| 亚洲久久久国产精品| 麻豆国产av国片精品| 日韩免费高清中文字幕av| 少妇被粗大的猛进出69影院| 嫁个100分男人电影在线观看| 日日夜夜操网爽| 一区福利在线观看| 一区在线观看完整版| 咕卡用的链子| 在线永久观看黄色视频| 69精品国产乱码久久久| 亚洲中文日韩欧美视频| 一级a爱片免费观看的视频| 超色免费av| 亚洲一码二码三码区别大吗| 亚洲欧洲精品一区二区精品久久久| 久久精品亚洲熟妇少妇任你| 窝窝影院91人妻| 757午夜福利合集在线观看| 不卡一级毛片| 亚洲第一av免费看| 亚洲 国产 在线| videosex国产| 老鸭窝网址在线观看| 国产亚洲欧美98| 可以在线观看毛片的网站| 悠悠久久av| 国产精品1区2区在线观看.| 国产野战对白在线观看| 国产精品 国内视频| 最新美女视频免费是黄的| a级毛片黄视频| a级片在线免费高清观看视频| 欧美日韩视频精品一区| 级片在线观看| 精品久久久精品久久久| 日本免费一区二区三区高清不卡 | 黄色视频,在线免费观看| 精品熟女少妇八av免费久了| 国产高清国产精品国产三级| 国产一区在线观看成人免费| 在线av久久热| 91字幕亚洲| 母亲3免费完整高清在线观看| 欧美一级毛片孕妇| 午夜免费成人在线视频| 99riav亚洲国产免费| 日韩国内少妇激情av| 久久人妻福利社区极品人妻图片| 国产激情欧美一区二区| 亚洲精品av麻豆狂野| 久久国产精品男人的天堂亚洲| 桃色一区二区三区在线观看| 神马国产精品三级电影在线观看 | 国产97色在线日韩免费| 国产一区二区三区视频了| 欧美在线一区亚洲| 亚洲一区二区三区不卡视频| 天天添夜夜摸| 久久精品国产综合久久久| 欧美激情极品国产一区二区三区| 夜夜躁狠狠躁天天躁| 1024香蕉在线观看| 正在播放国产对白刺激| 亚洲国产精品sss在线观看 | 成年版毛片免费区| 国产伦人伦偷精品视频| av电影中文网址| 琪琪午夜伦伦电影理论片6080| 亚洲国产中文字幕在线视频| 777久久人妻少妇嫩草av网站| 日韩精品青青久久久久久| 超碰成人久久| 身体一侧抽搐| 国产xxxxx性猛交| 视频区欧美日本亚洲| 性少妇av在线| 欧美日韩亚洲高清精品| 欧美日韩精品网址| 精品免费久久久久久久清纯| 又黄又粗又硬又大视频| 国产午夜精品久久久久久| 国产成人免费无遮挡视频| 91精品三级在线观看| 最近最新中文字幕大全电影3 | a在线观看视频网站| 亚洲熟女毛片儿| 久久性视频一级片| 国产一区二区三区视频了| 成人黄色视频免费在线看| 在线观看免费视频日本深夜| 日韩大码丰满熟妇| 一区二区三区激情视频| 午夜精品久久久久久毛片777| 欧美日韩黄片免| 香蕉国产在线看| 亚洲国产欧美日韩在线播放| 欧美人与性动交α欧美精品济南到| 欧美亚洲日本最大视频资源| 亚洲色图 男人天堂 中文字幕| 丝袜在线中文字幕| 国产成人一区二区三区免费视频网站| 久久久久亚洲av毛片大全| 精品国产一区二区久久| 日本五十路高清| 国产成人欧美在线观看| 亚洲国产看品久久| 美女扒开内裤让男人捅视频| 国产蜜桃级精品一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 18禁国产床啪视频网站| 怎么达到女性高潮| 午夜福利免费观看在线| 国产aⅴ精品一区二区三区波| 国产91精品成人一区二区三区| 黄片播放在线免费| 五月开心婷婷网| 精品卡一卡二卡四卡免费| 嫩草影院精品99| 色综合婷婷激情| 80岁老熟妇乱子伦牲交| 成人亚洲精品av一区二区 | 色综合婷婷激情| 国产国语露脸激情在线看| 色婷婷久久久亚洲欧美| 51午夜福利影视在线观看| 少妇的丰满在线观看| 国产三级在线视频| 久久精品国产亚洲av高清一级| 一级毛片精品| 久久天堂一区二区三区四区| 久久久久久久久中文| 俄罗斯特黄特色一大片| 女性被躁到高潮视频| 国产一区二区三区在线臀色熟女 | 国产高清视频在线播放一区| 桃色一区二区三区在线观看| 丝袜在线中文字幕| 亚洲人成伊人成综合网2020| 纯流量卡能插随身wifi吗| 精品久久蜜臀av无| 国产精品一区二区精品视频观看| 午夜免费激情av| 好男人电影高清在线观看| 极品教师在线免费播放| 99在线视频只有这里精品首页| 美女大奶头视频| 成人国产一区最新在线观看| 在线观看一区二区三区激情| 琪琪午夜伦伦电影理论片6080| 中文字幕人妻熟女乱码| netflix在线观看网站| 免费不卡黄色视频| 欧美激情久久久久久爽电影 | 日韩人妻精品一区2区三区| 国产在线精品亚洲第一网站| ponron亚洲| 精品国产乱码久久久久久男人| 国产在线观看jvid| 淫秽高清视频在线观看| 午夜老司机福利片| 国产欧美日韩一区二区三区在线| 精品国产亚洲在线| 国产91精品成人一区二区三区| 亚洲成a人片在线一区二区| 高清欧美精品videossex| 巨乳人妻的诱惑在线观看| 黄片大片在线免费观看| 国产成人欧美在线观看| 一区二区三区激情视频| 国产免费现黄频在线看| 亚洲成国产人片在线观看| 国产亚洲精品久久久久久毛片| 亚洲性夜色夜夜综合| 久久国产亚洲av麻豆专区| 丁香欧美五月| 日本vs欧美在线观看视频| 超碰成人久久| 十八禁人妻一区二区| 韩国av一区二区三区四区| 久久天躁狠狠躁夜夜2o2o| 日韩免费av在线播放| 在线观看日韩欧美| 国产又爽黄色视频| 国产欧美日韩一区二区精品| 亚洲国产看品久久| 国产99久久九九免费精品| 亚洲中文av在线| 免费日韩欧美在线观看| 中文欧美无线码| 九色亚洲精品在线播放| 一进一出抽搐动态| 日韩欧美一区视频在线观看| 国产成人啪精品午夜网站| 久久香蕉国产精品| 一区二区三区国产精品乱码| 天天添夜夜摸| 成人精品一区二区免费| 久久久久国内视频| 无遮挡黄片免费观看| 久久久精品国产亚洲av高清涩受| 精品人妻在线不人妻| 91在线观看av| 露出奶头的视频| 色老头精品视频在线观看| 国产极品粉嫩免费观看在线| 好男人电影高清在线观看| 免费av中文字幕在线| 国产精品自产拍在线观看55亚洲| 老熟妇乱子伦视频在线观看| 黄片大片在线免费观看| 久久国产乱子伦精品免费另类| 亚洲一区二区三区欧美精品| 国产视频一区二区在线看| 两性夫妻黄色片| 亚洲中文字幕日韩| 日韩有码中文字幕| 一区福利在线观看| 欧美大码av| 亚洲精品中文字幕一二三四区| 国产视频一区二区在线看| 交换朋友夫妻互换小说| 嫩草影视91久久| av片东京热男人的天堂| 十八禁网站免费在线| 极品人妻少妇av视频| 嫩草影视91久久| 国产欧美日韩一区二区精品| 免费女性裸体啪啪无遮挡网站| 中文字幕色久视频| 久久香蕉精品热| 日韩欧美三级三区| 国产熟女午夜一区二区三区| 日韩中文字幕欧美一区二区| 长腿黑丝高跟| 香蕉国产在线看| 欧美不卡视频在线免费观看 | 法律面前人人平等表现在哪些方面| 天天躁夜夜躁狠狠躁躁| 三级毛片av免费| 亚洲熟女毛片儿| 国产精品国产av在线观看| 在线播放国产精品三级| 日韩欧美国产一区二区入口| 国产精华一区二区三区| 亚洲中文字幕日韩| 高清黄色对白视频在线免费看| 精品卡一卡二卡四卡免费| 国产蜜桃级精品一区二区三区| 欧美黄色淫秽网站| 国产又爽黄色视频| 日本vs欧美在线观看视频| 夜夜看夜夜爽夜夜摸 | 夜夜夜夜夜久久久久|