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

    并列行駛的類車體氣動(dòng)特性數(shù)值模擬

    2022-04-12 01:42:51羅建斌米珂李銘森李龍杰韋信
    關(guān)鍵詞:數(shù)值模擬

    羅建斌 米珂 李銘森 李龍杰 韋信

    摘? 要:為了研究并列行駛工況下汽車的橫向間距對(duì)其氣動(dòng)特性的影響,采用后傾角為35°的Ahmed車體模型,結(jié)合定常的Realizable k-ε湍流模型進(jìn)行數(shù)值模擬研究。首先,針對(duì)單個(gè)Ahmed模型,以阻力系數(shù)為評(píng)價(jià)指標(biāo),通過與風(fēng)洞實(shí)驗(yàn)數(shù)據(jù)比較,對(duì)數(shù)值模擬方法進(jìn)行驗(yàn)證;然后,分析3輛車在并列行駛的情況下,其阻力系數(shù)、側(cè)向力系數(shù)和升力系數(shù)隨均勻間距變化的規(guī)律;最后,通過汽車車身表面壓力和尾部流場(chǎng)的分布分析了氣動(dòng)力發(fā)生變化的原因。結(jié)果表明:相比單輛車行駛工況,在各間距下,并列行駛的類車體阻力系數(shù)均有較大的增加,其中橫向間距為0.25W時(shí),第一輛車、第二輛車和第三輛車的阻力系數(shù)分別增加22%、35%和22%;第一輛車和第三輛車的阻力系數(shù)隨不同橫向間距的變化趨勢(shì)類似,而第二輛車的阻力系數(shù)始終大于其他2輛車;對(duì)于側(cè)向力系數(shù),即使未受到橫向側(cè)風(fēng)的影響,由于氣動(dòng)干擾,第一輛車和第三輛車均受到較大的側(cè)向力作用,而第二輛車由于流場(chǎng)結(jié)構(gòu)的對(duì)稱性,不管間距如何變化,均不受到側(cè)向力作用。研究結(jié)果可為智能交通系統(tǒng)規(guī)劃中對(duì)實(shí)際汽車并列行駛的氣動(dòng)特性分析提供工程借鑒和理論參考。

    關(guān)鍵詞:Ahmed模型;并列行駛;汽車空氣動(dòng)力學(xué);數(shù)值模擬

    中圖分類號(hào):U461.1? ? ? ? ? ? DOI:10.16375/j.cnki.cn45-1395/t.2022.02.001

    0? ? 引言

    對(duì)單輛車空氣動(dòng)力學(xué)特性進(jìn)行分析的目的是期望通過降低阻力來改善燃油經(jīng)濟(jì)性。目前,對(duì)單輛車的氣動(dòng)特性研究已較為成熟,對(duì)多車隊(duì)列行駛的研究表明其同樣可以降低車輛阻力,且可以緩解交通擁堵。汽車高速行駛時(shí)受到復(fù)雜氣流的作用會(huì)對(duì)汽車產(chǎn)生氣動(dòng)力和氣動(dòng)力矩(即氣動(dòng)六分力),極大地影響了汽車行駛的安全性和操縱穩(wěn)定性[1]。因此,在關(guān)注阻力變化的同時(shí),對(duì)側(cè)向力和升力的關(guān)注也同等重要。

    不同的橫向間距影響著汽車行駛時(shí)的阻力、側(cè)向力和升力,尤其是側(cè)向穩(wěn)定性,這對(duì)自動(dòng)駕駛、智能交通系統(tǒng)的開發(fā)和汽車主動(dòng)安全技術(shù)的研究至關(guān)重要。Vegendla等[2]通過數(shù)值模擬方法研究了不同排列方式下多個(gè)高速公路卡車的氣動(dòng)力影響,分析了2輛車和3輛車在單車道和多車道中5種不同配置下的氣動(dòng)特性,其中,多車道關(guān)注的是并列行駛的影響。目前國內(nèi)關(guān)于汽車空氣動(dòng)力學(xué)研究大致包括單輛車、會(huì)車、超車和隊(duì)列行駛等方面,其中隊(duì)列行駛[3-6]和并列行駛的情況涉及較少。傅立敏等[7]使用轎車模型研究了兩車在并列行駛時(shí)的湍流特性對(duì)交通安全的影響。王靖宇等[8]使用MIRA轎車模型研究了兩車并列行駛的氣動(dòng)特性,在速度分布、車身表面壓力分布等方面與單輛車外流場(chǎng)進(jìn)行了比較分析。陳洪業(yè)等[9]使用某汽車模型研究了汽車在不同間距下并列行駛的湍流特性。王征等[10]使用MIRA轎車模型研究了兩車并排行駛的氣動(dòng)特性,分析了兩車之間的流場(chǎng)分布和x方向上不同截面的尾流結(jié)構(gòu)。

    上述研究只涉及到2輛車并列行駛的情形,并且研究的間距范圍較窄,不能很好地反映出氣動(dòng)力隨間距變化的規(guī)律。因此,本文采用Ahmed模型,通過數(shù)值模擬研究3輛車在0.25W~4.0W(W為車寬)間距范圍內(nèi)行駛的氣動(dòng)特性的變化。主要探究在不同的橫向間距下,中間車在左右兩車相互作用下產(chǎn)生的變化,這可為汽車自動(dòng)駕駛和智能交通系統(tǒng)開發(fā)中的間距和車速控制情況提供理論基礎(chǔ),還可為多車在超車過程中發(fā)生的氣動(dòng)特性變化提供理論參考依據(jù)。

    1? ? 數(shù)值模擬方法設(shè)置

    1.1? ?模型和計(jì)算域的建立

    汽車實(shí)車模型曲面造型復(fù)雜,且底部結(jié)構(gòu)不規(guī)則,在數(shù)值模擬時(shí)常進(jìn)行簡(jiǎn)化處理。Ahmed模型是由Ahmed等[11]在1984年提出的一種類車體模型,由圓形前部、位于車身后部的可變化斜面及連接前部和后部斜面的長(zhǎng)方體組成,后部的可變化斜面主要用于研究不同傾斜角度時(shí)車體的阻力系數(shù)和尾流分離現(xiàn)象。該模型雖然簡(jiǎn)單,但保留了汽車車輛的基本特征,能比較準(zhǔn)確地反映汽車周圍流場(chǎng)的情況。在本文的數(shù)值模擬中,選用后部?jī)A角為35°的Ahmed模型,其實(shí)體建模形狀如圖1所示,主要結(jié)構(gòu)尺寸如圖2所示。

    在仿真計(jì)算中采用長(zhǎng)方體計(jì)算域,其中,入口距離為Ahmed模型的 2 倍長(zhǎng)度,出口距離為Ahmed模型的6倍長(zhǎng)度,頂部距離為 Ahmed 模型的6倍高度,Ahmed 模型兩側(cè)為3倍車寬,離地間隙為50 mm,如圖3所示,此時(shí)阻塞比為2%。

    1.2? ?網(wǎng)格劃分

    網(wǎng)格的類型和精度直接關(guān)系到計(jì)算的時(shí)間和準(zhǔn)確度,因此,網(wǎng)格劃分對(duì)于仿真模擬計(jì)算極其重要。因?yàn)锳hmed模型的幾何外形相對(duì)來說比較規(guī)整和簡(jiǎn)單,所以采用質(zhì)量較好的六面體結(jié)構(gòu)化網(wǎng)格。但其頭部為鈍型且曲率較大,結(jié)構(gòu)捕捉比較困難,因此,需要在其頭部單獨(dú)創(chuàng)建C-Block,可以更好地捕捉類車體頭部結(jié)構(gòu)。為了使網(wǎng)格在車體周圍更好地細(xì)化,在車體周圍生成O-Block之后,再生成一層C-Block,最后生成的網(wǎng)格如圖4所示。

    1.3? ?網(wǎng)格無關(guān)性驗(yàn)證

    一般而言,網(wǎng)格數(shù)量越多,計(jì)算結(jié)果的精度會(huì)越高,就能比較容易捕捉物體表面的細(xì)節(jié)特征。但是網(wǎng)格數(shù)量越多,需要的計(jì)算時(shí)間和計(jì)算內(nèi)存越多,就要求計(jì)算機(jī)的配置越高。由于計(jì)算時(shí)間和計(jì)算機(jī)硬件資源有限,因此,有必要進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證。

    本文僅對(duì)單輛車網(wǎng)格進(jìn)行獨(dú)立性驗(yàn)證。表1為網(wǎng)格無關(guān)性驗(yàn)證結(jié)果,表1中的5套網(wǎng)格類型都采用相同的網(wǎng)格拓?fù)浣Y(jié)構(gòu),將不同網(wǎng)格數(shù)量的阻力系數(shù)(CD)與文獻(xiàn)[12]的結(jié)果進(jìn)行對(duì)比。由表1可知,計(jì)算結(jié)果與網(wǎng)格數(shù)量存在一定的關(guān)系,網(wǎng)格數(shù)量較少的計(jì)算結(jié)果與網(wǎng)格數(shù)量較多的計(jì)算結(jié)果存在較大差異。Mesh 4與Mesh 5的計(jì)算結(jié)果非常相近,并且與實(shí)驗(yàn)數(shù)據(jù)的誤差最小。由于研究的是3輛車狀況,計(jì)算網(wǎng)格數(shù)量會(huì)隨著車輛間距的增大而增大,考慮到本次計(jì)算機(jī)的配置限度,故采用Mesh 3作為本次研究的網(wǎng)格方案。

    1.4? ?數(shù)值模擬參數(shù)設(shè)置及驗(yàn)證

    基于模型長(zhǎng)度的雷諾數(shù) Re=1.8×106 ,湍流強(qiáng)度為1.8%。湍流模型選用計(jì)算精度較高的高雷諾數(shù)Realizable k-ε模型[13],對(duì)于近壁區(qū)域的低雷諾數(shù)黏性流動(dòng),采用壁面函數(shù)法近似求解。邊界層內(nèi)第一層網(wǎng)格厚度△y =0.7 mm,此時(shí)y+分布在32~90,符合計(jì)算要求。壁面函數(shù)選用非平衡壁面函數(shù)(non-equilibrium wall functions),因?yàn)榉瞧胶獗诿婧瘮?shù)在計(jì)算時(shí)考慮了壁面附近的壓力梯度效應(yīng),當(dāng)流場(chǎng)涉及分離、再附著等情況時(shí)求解精度較高[14]。采用SIMPLEC壓力速度耦合算法,壓力插值為二階,控制方程的對(duì)流項(xiàng)和黏性項(xiàng)均采用二階離散格式。計(jì)算監(jiān)控殘差設(shè)為10-6,求解初始化從入口區(qū)用設(shè)定的初始值。當(dāng)所有監(jiān)控量(阻力、壓力)不隨計(jì)算的迭代次數(shù)而發(fā)生變化時(shí),認(rèn)為計(jì)算收斂,終止迭代。邊界條件設(shè)置如表2所示。

    表3為單輛車的模擬計(jì)算結(jié)果與實(shí)驗(yàn)值的比較。從表3可以看出,阻力系數(shù)的模擬值與實(shí)驗(yàn)值基本吻合,誤差僅為1.04%,而升力系數(shù)模擬值與實(shí)驗(yàn)值存在較大誤差,這是因?yàn)閷?shí)驗(yàn)時(shí)加裝了4個(gè)支撐車體的小圓柱,導(dǎo)致車底流動(dòng)變化。此外,目前的湍流模型并不能準(zhǔn)確地預(yù)測(cè)升力系數(shù),在相關(guān)研究[3]中,也同樣出現(xiàn)了升力系數(shù)預(yù)測(cè)不準(zhǔn)的? ? ? 情況。

    2? ? 結(jié)果與討論

    2.1? ?氣動(dòng)力系數(shù)分析

    在本次仿真計(jì)算中,為了研究橫向間距對(duì)汽車并列行駛的氣動(dòng)特性的影響,確定了0.25W、0.5W、0.75W、1.0W、1.5W、2.0W、2.5W、3.0W、3.5W和4.0W共10個(gè)車間距,其中W為車寬。圖5—圖7分別為阻力系數(shù)(CD)、升力系數(shù)(CL)和側(cè)向力系數(shù)(CS)隨車間距的變化圖,其中并列行駛時(shí),各個(gè)汽車的阻力系數(shù)和升力系數(shù)分別用單輛車的阻力系數(shù)(CD0)和升力系數(shù)(CL0)進(jìn)行標(biāo)準(zhǔn)化,這樣可以更好地比較汽車的氣動(dòng)力在并列行駛中的變化。在無側(cè)風(fēng)的情況下,單輛車側(cè)向力很小,幾乎為0,因此,對(duì)側(cè)向力系數(shù)不再進(jìn)行標(biāo)準(zhǔn)化處理。

    從圖5中可以看出,3輛車的阻力系數(shù)隨車間距的增加逐漸減小,并且都高于單輛車的阻力系數(shù)。第二輛車的阻力系數(shù)始終大于第一輛車和第三輛車。第一輛車和第三輛車的阻力系數(shù)幾乎一樣,說明兩車的流動(dòng)結(jié)構(gòu)和表面壓力的變化很相似。在車間距為0.25W時(shí),第二輛車的阻力系數(shù)比單輛車增加約35%。這是因?yàn)樵摴r下間距很近,車與車之間擠壓的高速氣流與第二輛車的尾流結(jié)構(gòu)相互影響,導(dǎo)致前后的壓差增加,從而阻力增加。隨著間距增加到4.0W,3輛車的阻力系數(shù)都在向單輛車值靠近,這時(shí)車與車之間的氣流干擾作用已經(jīng)很小。從圖6中可以看出,3輛車升力系數(shù)的變化趨勢(shì)一致。其中,第一輛車和第三輛車的升力系數(shù)變化曲線幾乎重疊,并且都大于第二輛車的升力系數(shù)。在車間距為0.25W~0.5W范圍內(nèi),3輛車的升力系數(shù)急劇下降,且都為負(fù)值,說明在這個(gè)間距范圍內(nèi),汽車的附著力增強(qiáng)。由圖7可以看出,第二輛車同時(shí)受到第一輛車和第三輛車的氣流擠壓作用,且由于3輛車的結(jié)構(gòu)一樣,擠壓作用對(duì)稱分布,大小相同,方向相反,正好相互抵消。因此,不管間距如何變化,其側(cè)向力系數(shù)接近于0。在車間距為0.25W~2.0W范圍內(nèi),第一輛車和第三輛車的側(cè)向力系數(shù)變化很大。在0.25W間距時(shí),第一輛車和第三輛車側(cè)向力系數(shù)分別比單輛車增加約25%和36%。因?yàn)殚g距較近,車與車兩側(cè)之間的氣流速度較高且相互干擾效應(yīng)非常強(qiáng),從而導(dǎo)致第一輛車和第三輛車的側(cè)向力都指向第二輛車。第一輛車和第三輛車均受到第二輛車的吸引力作用,這將直接影響汽車的安全性和操縱穩(wěn)定性。

    2.2? ?平均靜壓分析

    為了更好地體現(xiàn)出橫向間距下汽車行駛時(shí)的氣動(dòng)力變化,對(duì)車體表面的靜壓進(jìn)行分析。由圖8和圖9可以看出,車體頭部、尾部和車體左右側(cè)面的表面靜壓分布,靜壓系數(shù)[Cp]的表達(dá)式如下:

    [Cp=P?P00.5ρ0V20] ,? ? ? ? ? ? ? ? ? ? ? ? ? ?(1)

    式中:[P]是壓力,Pa;[P0]是自由流的參考?jí)毫?,相?dāng)于大氣壓力,Pa;[ρ0]是自由流的密度,相當(dāng)于空氣密度,kg/m3;[V0]是入口自由流速度,m/s。

    由圖8可以看出,在車間距為0.25W時(shí),3輛車頭部壓力均比單輛車頭部壓力增加,而尾部壓力相對(duì)于單輛車都有很大程度的降低。對(duì)于第二輛車,頭部壓力的增加和尾部壓力的降低都比其他2輛車多,因此,其阻力系數(shù)最高。當(dāng)車間距為0.75W時(shí),也出現(xiàn)了類似的現(xiàn)象。當(dāng)車間距進(jìn)一步增加到2.0W時(shí),3輛車頭部壓力的變化與單輛車相比幾乎沒有改變,但是3輛車尾部壓力的變化相比于單輛車仍然有所降低。在車間距為4.0W時(shí),由于氣流的相互作用變得很弱,3輛車頭部和尾部的壓力分布同單輛車的壓力分布一樣,因此,阻力系數(shù)也在向單輛車的值靠近。

    圖9給出了并列行駛中3輛車側(cè)向靜壓分布隨車間距的變化。在車間距為0.25W時(shí),第一輛車的右側(cè)面有明顯的壓力增加,左側(cè)面(靠近第二輛車的一面)有明顯的壓力降低,并且負(fù)壓的區(qū)域有所增加。第三輛車的側(cè)面壓力分布正好與第一輛車相反,右側(cè)面(靠近第二輛車的一面)壓力降低,左側(cè)面壓力增加,導(dǎo)致它的側(cè)向力是負(fù)號(hào),其合力方向指向第二輛車。因此,在這個(gè)間距下行駛,第一輛車和第三輛車均受到朝向第二輛車的吸引力作用。這樣會(huì)大大影響汽車的安全性和操縱穩(wěn)定性,輕微的方向偏移就會(huì)導(dǎo)致汽車行駛軌跡的劇烈改變。而對(duì)于第二輛車,不管處于哪種間距下,其左右側(cè)面靜壓壓力均呈現(xiàn)出對(duì)稱分布,因此,所受側(cè)向力均為0。車間距為0.75W時(shí),兩側(cè)面壓力變化的趨勢(shì)與車間距為0.25W時(shí)的變化非常相似,只是壓力分布大小的程度比車間距為0.25W時(shí)有所降低。當(dāng)車間距增加到2.0W和4.0W時(shí),3輛車兩側(cè)面的壓力分布幾乎一樣,這種現(xiàn)象在車間距為4.0W時(shí)更為顯著。

    2.3? ?平均流場(chǎng)分析

    壓力變化主要是因?yàn)榱鲃?dòng)結(jié)構(gòu)改變而造成的。圖10給出了單輛車及3輛車在z/H = 0.2平面上(H為車高)的流線圖和流向速度云圖,其中流向速度(V)用入口自由流速度(V0)進(jìn)行標(biāo)準(zhǔn)化。

    從圖10中可以看出,單輛車尾流結(jié)構(gòu)中有上下2個(gè)對(duì)稱且大小相等的旋渦區(qū)域。在車間距為0.25W時(shí),3輛車的尾流結(jié)構(gòu)均發(fā)生了改變。其中第一輛車的上下2個(gè)旋渦不再對(duì)稱,并且下旋渦變小,兩車之間沖出的高速氣流導(dǎo)致尾流閉合線向下偏移,因此,其側(cè)向力指向第二輛車。第三輛車尾流結(jié)構(gòu)的變化與第一輛車中的變化非常相似,這可以對(duì)應(yīng)阻力系數(shù)(圖5)和側(cè)向力系數(shù)(圖7)的變化。第二輛車中的尾流結(jié)構(gòu)沒有發(fā)生變化,因?yàn)樗艿?輛車之間共同的擠壓作用,只是尾流中的旋渦區(qū)域比單輛車中的旋渦區(qū)域大,因此,其阻力系數(shù)在0.25W間距下變得比較高。這里側(cè)向力增加的原因可以運(yùn)用伯努利定律來分析,在0.25W間距中,車與車之間的氣流速度增加,導(dǎo)致壓力降低,因此,車體兩側(cè)形成較高的壓力差,其側(cè)向力增加。

    在車間距為0.75W時(shí),3輛車尾流結(jié)構(gòu)的變化情況與車間距為0.25W時(shí)的變化情況一樣,只是車與車之間的氣流速度開始降低,因此,其側(cè)向力系數(shù)也開始降低。當(dāng)間距增加到2.0W時(shí),3輛車中的尾流結(jié)構(gòu)變得幾乎與單輛車一樣,第一輛車和第二輛車的上下旋渦只有極其微小的變化,這說明車與車之間的氣流相互作用已經(jīng)很弱。當(dāng)最大車間距為4.0W時(shí),3輛車變得更加獨(dú)立,車與車之間的氣流幾乎沒有相互作用。因此,其阻力系數(shù)和側(cè)向力系數(shù)均向單輛車的值靠近。

    圖11為單輛車及3輛車在x/L= 0.522的平面上的流線圖和壓力云圖。從圖11中可以看出,在單輛車狀況下,車體兩側(cè)有4個(gè)左右相互對(duì)稱的旋渦,上旋渦比下旋渦要小,且車體周圍的壓力分布沒有發(fā)生變化。當(dāng)3輛車并列行駛時(shí),車與車之間的流動(dòng)結(jié)構(gòu)和靜壓分布發(fā)生較大的改變。在間距為0.25W時(shí),第一輛車和第三輛車中的4個(gè)旋渦變成3個(gè)旋渦,其大小和位置也發(fā)生了改變。第二輛車的4個(gè)旋渦依然是左右對(duì)稱分布,但是車體兩側(cè)的一對(duì)上旋渦移動(dòng)到車體頂部,且變得比下面的旋渦大。由于3輛車間距很近,氣流流動(dòng)的空間突然變小,速度增加,車與車之間的相互擠壓作用增強(qiáng),因此,車與車間隙中的旋渦被擠壓到車體頂部。同時(shí),車與車之間存在著較高的負(fù)壓區(qū),第一輛車和第三輛車兩側(cè)形成較高的壓力差,因此,2輛車的側(cè)向力系數(shù)變得比較高。在間距為0.75W時(shí),第一輛車和第三輛車又形成了4個(gè)旋渦,第二輛車中的旋渦變得與0.25W間距中的相反,其旋渦大小有明顯的改變,說明車與車之間的相互擠壓作用開始變?nèi)酰虼藗?cè)向力系數(shù)開始降低。隨著間距進(jìn)一步增加到2.0W和4.0W,3輛車中的4個(gè)旋渦在大小和位置上都變得與單輛車非常相似。無論在哪種間距下,第一輛車和第三輛車中的旋渦變化幾乎是一致的,第二輛車中的旋渦始終是左右對(duì)稱的,這種旋渦變化情況可以反映到圖7中3輛車側(cè)向力系數(shù)的變化。

    3? ? 結(jié)論

    本文在不同的橫向間距下對(duì)3個(gè)后傾角為35°的Ahmed 模型進(jìn)行了并列行駛的數(shù)值模擬研究,結(jié)論如下:

    1)3輛車并列行駛時(shí),第一輛車和第三輛車中壓力分布和流動(dòng)結(jié)構(gòu)的變化基本一致。第二輛車在第一輛車和第三輛車共同的擠壓作用下,左右壓力分布和流動(dòng)結(jié)構(gòu)呈現(xiàn)對(duì)稱分布,因此,其側(cè)向力系數(shù)幾乎不變。

    2)相比于單輛車行駛,并列行駛工況下汽車所受到的阻力系數(shù)均有所增加,正好與隊(duì)列行駛中汽車阻力系數(shù)減小的趨勢(shì)相反。在間距很近的并列行駛中,車與車之間沖出的高速氣流與尾流相互影響很大,導(dǎo)致3輛車中的尾流結(jié)構(gòu)變得復(fù)雜,使得前后壓差增加,阻力系數(shù)變大。

    3)不同編排形式的車輛在道路上行駛時(shí)會(huì)產(chǎn)生不同的流動(dòng)結(jié)構(gòu)。隊(duì)列行駛主要影響的是汽車頭、尾部的流動(dòng)變化,而并列行駛主要影響的是汽車側(cè)部的流動(dòng)變化,導(dǎo)致隊(duì)列行駛和并列行駛的汽車產(chǎn)生不同的氣動(dòng)力變化。

    參考文獻(xiàn)

    [1]? ? ?傅立敏.汽車空氣動(dòng)力學(xué)[M].北京:機(jī)械工業(yè)出版社,1998.

    [2]? ? ?VEGENDLA P,SOFU T,SAHA R,et al.Investigation of aerodynamic influence on truck platooning[C]//SAE Commercial Vehicle Engineering Congress,2015.

    [3]? ? ?傅立敏,吳允柱,賀寶琴.隊(duì)列行駛車輛的空氣動(dòng)力特性[J].吉林大學(xué)學(xué)報(bào)(工學(xué)版),2006,36(6):871-875.

    [4]? ? ?高偉,余偉,吳永興,等.隊(duì)列行駛貨車氣動(dòng)特性仿真及影響因素分析[J].湖北汽車工業(yè)學(xué)院學(xué)報(bào),2020,34(4):33-37.

    [5]? ? ?羅建斌,吳量,苗明達(dá),等.汽車隊(duì)列行駛氣動(dòng)特性數(shù)值模擬[J].廣西科技大學(xué)學(xué)報(bào),2018,29(1):8-13.

    [6]? ? ?張英朝,王中儉,蘇暢,等.摩托車隊(duì)列行駛的空氣動(dòng)力特性[C]//2020中國汽車工程學(xué)會(huì)年會(huì)暨展覽會(huì),2020.

    [7]? ? ?傅立敏,扶原放.轎車并列行駛湍流特性的數(shù)值模擬[J].吉林大學(xué)學(xué)報(bào)(工學(xué)版),2005,35(4):358-362.

    [8]? ? ?王靖宇,劉順航,胡興軍,等.并列行駛兩輛汽車外流場(chǎng)的數(shù)值模擬[J].沈陽理工大學(xué)學(xué)報(bào),2007,26(4):83-87.

    [9]? ? ?陳洪業(yè),王愛玲,鄭智貞,等.汽車并列行駛湍流數(shù)值模擬[J].機(jī)械設(shè)計(jì)與制造,2010(4):196-197.

    [10]? ?王征,劉順航.兩車并排行駛空氣動(dòng)力學(xué)特性的研究[J].科學(xué)之友,2011(13):18-20.

    [11]? ?AHMED S R,RAMM G,F(xiàn)AITIN G.Some salient features of the time-averaged ground vehicle wake[C]//SAE International Congress and Exposition,1984.

    [12]? ?STRACHAN R K,KNOWLES K,LAWSON N J.The vortex structure behind an Ahmed reference model in the presence of a moving ground plane[J]. Experiments in Fluids,2007,42(5):659-669.

    [13]? ?梁建永,梁軍,范士杰,等.轎車外流場(chǎng)CFD分析中常用k-ε湍流模型的對(duì)比[J].汽車工程,2008(10):846-852.

    [14]? ?賴晨光,王擎宇,付中正.基于湍流模型與邊界層網(wǎng)格策略的汽車外流場(chǎng)數(shù)值計(jì)算精度分析[C]//2017汽車空氣動(dòng)力學(xué)分會(huì)學(xué)術(shù)年會(huì),2017.

    Numerical simulation of aerodynamic characteristics of Ahmed model driving side by side

    LUO Jianbin1, MI Ke1, LI Mingsen1, LI Longjie2, WEI Xin3

    (1.School of Mechanical and Automotive Engineering, Guangxi University of Science and Technology, Liuzhou 545616, China; 2. Liuzhou Wuling Automobile Industry Co., Ltd., Liuzhou 545007, China;

    3. Shanghai Automotive Industry Corp General Motors Wuling Automobile Co., Ltd., Liuzhou 545007, China)

    Abstract: The numerical simulation is conducted with steady Realizable k-ε turbulence model to study the effect of lateral spacing on the aerodynamic characteristics of vehicles under the condition of? ? ? driving side by side, using the Ahmed model with a rear slant angle of 35°. Firstly, for a single Ahmed model, the numerical simulation method is verified with the drag coefficient as evaluation index through comparing with the wind tunnel experimental data. Secondly, the variation law of drag? ? ? ? ? ?coefficient, side force coefficient and lift coefficient with uniform spacing is analyzed in the case of three vehicles driving side by side. Finally, the causes for aerodynamic changes are analyzed by the? ?distribution of surface pressure and tail flow field. The results show that, compared with the single car driving condition, the drag coefficient of the Ahmed model driving side by side is greatly increased at each spacing. When the lateral spacing is 0.25W, the drag coefficient of the first car, the second car and the third car increases by 22%, 35% and 22%, respectively. The change trend of the drag coefficient of the first car and the third car with different lateral spacing is similar, while the drag coefficient of the second car is always higher than that of the other two cars. Regarding the side force coefficient, even if it is not affected by the lateral crosswind, due to aerodynamic interference, both the first car and the? second car are subjected to greater side force. However, the second car is not subjected to the side force no matter how the spacing changes because of the symmetry of the flow field structure. The results of the study can provide engineering and theoretical reference for the analysis of aerodynamic? ? ? ? ? ? ? characteristics of actual cars driving side by side in intelligent transportation system planning.

    Key words: Ahmed model; driving side by side; vehicle aerodynamics; numerical simulation

    (責(zé)任編輯:黎? ?婭)

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場(chǎng)信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測(cè)
    鉆孔灌注樁樁底沉渣對(duì)樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測(cè)試、預(yù)處理和加工研究報(bào)告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場(chǎng)數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    好看av亚洲va欧美ⅴa在| 熟女电影av网| 日韩精品免费视频一区二区三区| 51午夜福利影视在线观看| 免费在线观看黄色视频的| 久久国产精品人妻蜜桃| 国产亚洲av高清不卡| avwww免费| 久久精品国产亚洲av高清一级| 久久久久性生活片| 欧美日韩国产亚洲二区| 搡老岳熟女国产| 午夜福利免费观看在线| 97超级碰碰碰精品色视频在线观看| 国内精品久久久久精免费| 又紧又爽又黄一区二区| 国产私拍福利视频在线观看| 在线观看午夜福利视频| 正在播放国产对白刺激| 一边摸一边做爽爽视频免费| 久久天堂一区二区三区四区| 又大又爽又粗| 男女那种视频在线观看| 黄色 视频免费看| 美女黄网站色视频| 亚洲欧美日韩高清专用| 99精品在免费线老司机午夜| 99久久久亚洲精品蜜臀av| 国产精品免费一区二区三区在线| 久久精品国产亚洲av高清一级| 一个人免费在线观看电影 | 女警被强在线播放| 最近最新中文字幕大全电影3| 天天一区二区日本电影三级| 久久精品aⅴ一区二区三区四区| 男女之事视频高清在线观看| 欧美极品一区二区三区四区| 亚洲真实伦在线观看| 欧美日韩中文字幕国产精品一区二区三区| 两性夫妻黄色片| 欧美久久黑人一区二区| 女同久久另类99精品国产91| 亚洲av成人av| 99热6这里只有精品| 精品久久久久久久人妻蜜臀av| 久久婷婷成人综合色麻豆| 美女高潮喷水抽搐中文字幕| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美日韩高清在线视频| 嫩草影视91久久| 99精品欧美一区二区三区四区| 久久精品成人免费网站| 高潮久久久久久久久久久不卡| 亚洲中文字幕日韩| 国产一区二区激情短视频| 精品熟女少妇八av免费久了| 在线国产一区二区在线| 国产一区在线观看成人免费| 麻豆av在线久日| 正在播放国产对白刺激| 久久久久久九九精品二区国产 | netflix在线观看网站| 97碰自拍视频| 丝袜美腿诱惑在线| www国产在线视频色| 超碰成人久久| av国产免费在线观看| а√天堂www在线а√下载| 免费看日本二区| 国产精品乱码一区二三区的特点| 国产不卡一卡二| 国产成人精品久久二区二区91| 91麻豆av在线| 国产精品99久久99久久久不卡| 久久99热这里只有精品18| 午夜福利在线观看吧| 欧洲精品卡2卡3卡4卡5卡区| 变态另类成人亚洲欧美熟女| 久久精品影院6| 99精品久久久久人妻精品| 国产精品电影一区二区三区| 国产野战对白在线观看| 亚洲aⅴ乱码一区二区在线播放 | 露出奶头的视频| 免费在线观看成人毛片| 日韩国内少妇激情av| 99re在线观看精品视频| a级毛片在线看网站| 日韩欧美国产在线观看| 午夜福利高清视频| 欧美日韩乱码在线| 一卡2卡三卡四卡精品乱码亚洲| 日韩精品中文字幕看吧| 黄片小视频在线播放| 两性夫妻黄色片| 亚洲精品av麻豆狂野| 神马国产精品三级电影在线观看 | av免费在线观看网站| 精品国产美女av久久久久小说| 嫁个100分男人电影在线观看| 黑人欧美特级aaaaaa片| 老汉色∧v一级毛片| 黄色视频不卡| 久久久久国产一级毛片高清牌| 又黄又爽又免费观看的视频| 日韩欧美三级三区| 三级毛片av免费| 欧美日韩国产亚洲二区| 久久天堂一区二区三区四区| 午夜老司机福利片| 黑人操中国人逼视频| 老司机在亚洲福利影院| 日本免费一区二区三区高清不卡| 美女午夜性视频免费| 亚洲男人天堂网一区| 欧美性猛交╳xxx乱大交人| 国内精品久久久久久久电影| 国产熟女xx| 久久久久久人人人人人| 女生性感内裤真人,穿戴方法视频| 制服诱惑二区| 国产精品一区二区免费欧美| 一二三四在线观看免费中文在| 久久人人精品亚洲av| 亚洲精品粉嫩美女一区| 久久人妻福利社区极品人妻图片| 亚洲精品久久国产高清桃花| 精品欧美国产一区二区三| 久久久久性生活片| 成人18禁在线播放| 真人做人爱边吃奶动态| 黑人巨大精品欧美一区二区mp4| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲熟妇熟女久久| 丝袜人妻中文字幕| 性欧美人与动物交配| 精品国内亚洲2022精品成人| 国产黄片美女视频| 亚洲国产欧洲综合997久久,| 欧美丝袜亚洲另类 | 在线视频色国产色| 欧美色视频一区免费| 亚洲国产欧美网| 97人妻精品一区二区三区麻豆| 免费搜索国产男女视频| 日韩大尺度精品在线看网址| 成人三级黄色视频| 禁无遮挡网站| 欧美 亚洲 国产 日韩一| 国产免费av片在线观看野外av| 制服诱惑二区| 国产免费男女视频| 亚洲,欧美精品.| 中文字幕精品亚洲无线码一区| 天天一区二区日本电影三级| 国产一区二区在线av高清观看| 国产日本99.免费观看| www日本在线高清视频| 丰满的人妻完整版| 欧美黄色淫秽网站| 在线观看www视频免费| av有码第一页| 成熟少妇高潮喷水视频| 我的老师免费观看完整版| 久久久精品欧美日韩精品| 蜜桃久久精品国产亚洲av| 亚洲 欧美 日韩 在线 免费| 在线观看美女被高潮喷水网站 | 99在线视频只有这里精品首页| 免费在线观看日本一区| 男人舔女人下体高潮全视频| 黄色片一级片一级黄色片| 日本熟妇午夜| 19禁男女啪啪无遮挡网站| 国产成人欧美在线观看| 久久久久久九九精品二区国产 | 不卡一级毛片| 国产成人精品无人区| 两个人看的免费小视频| 国产成人系列免费观看| 国产成人aa在线观看| www日本在线高清视频| 欧美黑人欧美精品刺激| 97超级碰碰碰精品色视频在线观看| 国产精品99久久99久久久不卡| 九色成人免费人妻av| 久久久国产成人精品二区| 日韩精品中文字幕看吧| 男人舔女人下体高潮全视频| or卡值多少钱| 久久这里只有精品19| 婷婷亚洲欧美| 91成年电影在线观看| 国产激情久久老熟女| 免费在线观看影片大全网站| 欧美成人性av电影在线观看| 亚洲精品av麻豆狂野| 美女午夜性视频免费| tocl精华| 在线a可以看的网站| 视频区欧美日本亚洲| 熟女少妇亚洲综合色aaa.| 日本黄大片高清| 一夜夜www| 久久热在线av| 老汉色∧v一级毛片| 麻豆国产97在线/欧美 | 婷婷亚洲欧美| 久久天堂一区二区三区四区| 在线观看午夜福利视频| 国产成人av激情在线播放| 精品电影一区二区在线| videosex国产| 非洲黑人性xxxx精品又粗又长| 欧美性长视频在线观看| 香蕉av资源在线| 亚洲五月婷婷丁香| 亚洲欧美日韩高清在线视频| 97碰自拍视频| 欧美一区二区国产精品久久精品 | 黄色丝袜av网址大全| 手机成人av网站| 欧美黑人精品巨大| 麻豆av在线久日| 精品第一国产精品| 国产99白浆流出| 亚洲精品国产精品久久久不卡| 久久香蕉激情| 看片在线看免费视频| 91字幕亚洲| a级毛片a级免费在线| 亚洲免费av在线视频| 无限看片的www在线观看| 丁香欧美五月| 母亲3免费完整高清在线观看| 亚洲男人的天堂狠狠| 久热爱精品视频在线9| 他把我摸到了高潮在线观看| 国内精品久久久久久久电影| 日韩欧美 国产精品| 一级片免费观看大全| 成人精品一区二区免费| 长腿黑丝高跟| 男女那种视频在线观看| 精品国产美女av久久久久小说| 国产97色在线日韩免费| 亚洲 国产 在线| 免费看十八禁软件| 久久久久国产精品人妻aⅴ院| 久久久久久国产a免费观看| 日本三级黄在线观看| 亚洲真实伦在线观看| 欧美乱妇无乱码| 又黄又爽又免费观看的视频| 少妇裸体淫交视频免费看高清 | 国产91精品成人一区二区三区| 又爽又黄无遮挡网站| 一区二区三区高清视频在线| 精品久久久久久久末码| 亚洲自偷自拍图片 自拍| 免费高清视频大片| 国产精品一区二区免费欧美| 欧美午夜高清在线| 法律面前人人平等表现在哪些方面| 97碰自拍视频| 女警被强在线播放| 久久精品91蜜桃| 色播亚洲综合网| 亚洲国产欧美网| 国产高清视频在线观看网站| 中文字幕最新亚洲高清| 一区福利在线观看| av片东京热男人的天堂| 成人国产一区最新在线观看| 99在线视频只有这里精品首页| 国产精品免费视频内射| 男女视频在线观看网站免费 | 一边摸一边抽搐一进一小说| 又黄又爽又免费观看的视频| 亚洲av第一区精品v没综合| 一级毛片高清免费大全| 日韩三级视频一区二区三区| 免费观看精品视频网站| 亚洲国产欧洲综合997久久,| 国产在线精品亚洲第一网站| 久久99热这里只有精品18| 狂野欧美激情性xxxx| 天天躁夜夜躁狠狠躁躁| 精品熟女少妇八av免费久了| 12—13女人毛片做爰片一| 日韩大码丰满熟妇| av福利片在线观看| 免费搜索国产男女视频| 搞女人的毛片| 一个人免费在线观看电影 | 国产午夜精品论理片| 欧美黑人欧美精品刺激| 老司机福利观看| 在线观看免费视频日本深夜| 亚洲aⅴ乱码一区二区在线播放 | 精品欧美国产一区二区三| av福利片在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲自偷自拍图片 自拍| 一个人免费在线观看的高清视频| 黄色毛片三级朝国网站| av有码第一页| 波多野结衣巨乳人妻| 欧美黑人欧美精品刺激| 精品久久久久久久久久免费视频| 草草在线视频免费看| www国产在线视频色| 午夜免费成人在线视频| aaaaa片日本免费| 成人av一区二区三区在线看| 亚洲色图 男人天堂 中文字幕| 婷婷精品国产亚洲av| 两个人视频免费观看高清| 亚洲成av人片在线播放无| 麻豆国产97在线/欧美 | 美女高潮喷水抽搐中文字幕| 日本一本二区三区精品| 丰满人妻一区二区三区视频av | 久久久久国内视频| avwww免费| 中出人妻视频一区二区| 午夜影院日韩av| 国产麻豆成人av免费视频| 久热爱精品视频在线9| 免费人成视频x8x8入口观看| 欧美精品亚洲一区二区| 69av精品久久久久久| 一本综合久久免费| 黄频高清免费视频| www日本在线高清视频| 日韩欧美免费精品| 不卡av一区二区三区| 我要搜黄色片| 国产精品乱码一区二三区的特点| 国产高清有码在线观看视频 | 琪琪午夜伦伦电影理论片6080| 91老司机精品| 天堂√8在线中文| 国产成人av激情在线播放| 亚洲成人国产一区在线观看| 午夜成年电影在线免费观看| 黄色女人牲交| 久久国产精品人妻蜜桃| 中文字幕人妻丝袜一区二区| 国产成人av激情在线播放| 日韩欧美三级三区| 中文字幕久久专区| 亚洲精品美女久久av网站| 九色国产91popny在线| 一进一出好大好爽视频| 国产精品一区二区免费欧美| 老汉色av国产亚洲站长工具| 可以在线观看毛片的网站| 国产精品一区二区三区四区久久| 亚洲 欧美一区二区三区| 久久精品影院6| 高潮久久久久久久久久久不卡| 中文字幕人妻丝袜一区二区| 可以在线观看毛片的网站| 国产黄a三级三级三级人| 啦啦啦观看免费观看视频高清| 国产成年人精品一区二区| 国产又黄又爽又无遮挡在线| 男女床上黄色一级片免费看| 窝窝影院91人妻| 国产一区二区激情短视频| 日日摸夜夜添夜夜添小说| 女生性感内裤真人,穿戴方法视频| 国产探花在线观看一区二区| 国产精品九九99| 精品国内亚洲2022精品成人| 亚洲无线在线观看| 中文资源天堂在线| 日日干狠狠操夜夜爽| 在线十欧美十亚洲十日本专区| 久久精品国产亚洲av高清一级| 欧美日韩瑟瑟在线播放| 日本 欧美在线| 一级黄色大片毛片| 亚洲成av人片免费观看| 老熟妇仑乱视频hdxx| 久99久视频精品免费| 好男人电影高清在线观看| 欧美日韩一级在线毛片| 三级男女做爰猛烈吃奶摸视频| 香蕉丝袜av| 亚洲av成人精品一区久久| 亚洲av成人不卡在线观看播放网| 岛国在线免费视频观看| 夜夜躁狠狠躁天天躁| avwww免费| 亚洲欧美激情综合另类| 亚洲国产高清在线一区二区三| 国产亚洲精品久久久久5区| 欧美一级a爱片免费观看看 | 国产激情久久老熟女| 国产av在哪里看| 麻豆成人午夜福利视频| 国产乱人伦免费视频| 俄罗斯特黄特色一大片| 成人一区二区视频在线观看| 国产av一区二区精品久久| av片东京热男人的天堂| 老汉色∧v一级毛片| 夜夜爽天天搞| 老司机深夜福利视频在线观看| 性色av乱码一区二区三区2| 日韩精品中文字幕看吧| 亚洲精华国产精华精| 国产视频一区二区在线看| 久久性视频一级片| 欧美+亚洲+日韩+国产| 国产亚洲欧美98| 巨乳人妻的诱惑在线观看| 日韩欧美在线二视频| 91麻豆精品激情在线观看国产| 最近最新中文字幕大全电影3| 天堂动漫精品| 精品高清国产在线一区| 在线观看美女被高潮喷水网站 | 91av网站免费观看| 国产精品一区二区精品视频观看| 成年人黄色毛片网站| 精品电影一区二区在线| 又爽又黄无遮挡网站| 色噜噜av男人的天堂激情| 黄色丝袜av网址大全| 国产高清视频在线观看网站| 国产亚洲精品综合一区在线观看 | 亚洲国产高清在线一区二区三| 99在线视频只有这里精品首页| 国产三级中文精品| 91麻豆av在线| 精品乱码久久久久久99久播| 丁香六月欧美| 波多野结衣高清作品| 一本精品99久久精品77| 久久国产乱子伦精品免费另类| 麻豆成人午夜福利视频| 成人av一区二区三区在线看| 国产亚洲欧美98| 男女下面进入的视频免费午夜| 久久午夜综合久久蜜桃| 黄色成人免费大全| 亚洲七黄色美女视频| 一二三四在线观看免费中文在| 给我免费播放毛片高清在线观看| 欧美+亚洲+日韩+国产| 久久久久性生活片| 国产三级中文精品| 欧美人与性动交α欧美精品济南到| 在线观看66精品国产| 国产黄色小视频在线观看| 老司机深夜福利视频在线观看| 最近最新免费中文字幕在线| 我的老师免费观看完整版| 97超级碰碰碰精品色视频在线观看| 日韩 欧美 亚洲 中文字幕| 制服人妻中文乱码| 欧美性猛交黑人性爽| 亚洲欧美日韩东京热| a在线观看视频网站| ponron亚洲| 首页视频小说图片口味搜索| 欧美午夜高清在线| 亚洲中文日韩欧美视频| 免费人成视频x8x8入口观看| 身体一侧抽搐| 国产成人欧美在线观看| 亚洲,欧美精品.| 中文字幕高清在线视频| 两个人视频免费观看高清| 日韩中文字幕欧美一区二区| 久久久久性生活片| 黄色女人牲交| 日日爽夜夜爽网站| 黄频高清免费视频| 精品国产亚洲在线| 久久久久性生活片| 国产精品久久久人人做人人爽| 日日夜夜操网爽| 国产精品爽爽va在线观看网站| 中文字幕高清在线视频| 琪琪午夜伦伦电影理论片6080| 亚洲av第一区精品v没综合| 两个人看的免费小视频| 长腿黑丝高跟| 欧美中文日本在线观看视频| 亚洲专区中文字幕在线| 午夜精品在线福利| 久久婷婷成人综合色麻豆| 国产av又大| 两个人的视频大全免费| 一本久久中文字幕| 日本黄色视频三级网站网址| 黄频高清免费视频| 中文字幕高清在线视频| 动漫黄色视频在线观看| 午夜精品在线福利| www.自偷自拍.com| 国产69精品久久久久777片 | 男人舔奶头视频| 在线观看日韩欧美| 亚洲成av人片在线播放无| 老汉色av国产亚洲站长工具| 夜夜躁狠狠躁天天躁| 99精品欧美一区二区三区四区| 国产在线精品亚洲第一网站| 精品乱码久久久久久99久播| 波多野结衣巨乳人妻| 夜夜爽天天搞| 欧美+亚洲+日韩+国产| 亚洲国产中文字幕在线视频| 亚洲中文av在线| 欧美激情久久久久久爽电影| 久9热在线精品视频| 亚洲欧美日韩无卡精品| 淫妇啪啪啪对白视频| 久久中文字幕人妻熟女| 午夜两性在线视频| 欧美丝袜亚洲另类 | 成人手机av| 午夜影院日韩av| 九色国产91popny在线| 亚洲九九香蕉| 国产午夜福利久久久久久| 亚洲成a人片在线一区二区| 国产私拍福利视频在线观看| 亚洲午夜理论影院| 免费看十八禁软件| 国产野战对白在线观看| 91老司机精品| 欧美成人性av电影在线观看| 午夜视频精品福利| 国产精品久久电影中文字幕| videosex国产| 女同久久另类99精品国产91| 久久久久性生活片| 99热只有精品国产| 搞女人的毛片| 一级作爱视频免费观看| 男女下面进入的视频免费午夜| 欧美久久黑人一区二区| 19禁男女啪啪无遮挡网站| 在线国产一区二区在线| 男女视频在线观看网站免费 | 亚洲成a人片在线一区二区| 国产精品影院久久| xxxwww97欧美| 黄色女人牲交| 亚洲国产日韩欧美精品在线观看 | 18禁美女被吸乳视频| 不卡一级毛片| 99热这里只有精品一区 | 亚洲aⅴ乱码一区二区在线播放 | 亚洲第一欧美日韩一区二区三区| 99精品欧美一区二区三区四区| 天天一区二区日本电影三级| 真人一进一出gif抽搐免费| av在线播放免费不卡| 中文字幕高清在线视频| 18禁裸乳无遮挡免费网站照片| 午夜免费成人在线视频| 91老司机精品| 日本a在线网址| 悠悠久久av| 熟女少妇亚洲综合色aaa.| 国产一区二区激情短视频| 在线观看66精品国产| 国产69精品久久久久777片 | 麻豆国产av国片精品| 午夜福利欧美成人| 色综合婷婷激情| 久久天躁狠狠躁夜夜2o2o| 国产野战对白在线观看| 免费观看精品视频网站| 国产精品久久久久久人妻精品电影| 伊人久久大香线蕉亚洲五| 中文字幕最新亚洲高清| 香蕉丝袜av| 午夜视频精品福利| 久久久久精品国产欧美久久久| av在线天堂中文字幕| 久久婷婷成人综合色麻豆| 熟女少妇亚洲综合色aaa.| 中文字幕熟女人妻在线| 欧美日韩黄片免| 校园春色视频在线观看| 舔av片在线| 日本成人三级电影网站| 夜夜夜夜夜久久久久| 99在线人妻在线中文字幕| 在线观看66精品国产| 国产97色在线日韩免费| 亚洲精品粉嫩美女一区| 亚洲专区国产一区二区| 夜夜夜夜夜久久久久| 变态另类丝袜制服| 婷婷亚洲欧美| xxxwww97欧美| 日韩三级视频一区二区三区| 黄色成人免费大全| 成年女人毛片免费观看观看9| 黄色视频不卡| 又紧又爽又黄一区二区| 久久久久国产一级毛片高清牌| 欧美性猛交黑人性爽| 日韩中文字幕欧美一区二区| 熟女少妇亚洲综合色aaa.| 亚洲九九香蕉|