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

    FSAE賽車彎道氣動特性數值模擬與試驗研究

    2023-11-15 06:40:16賴晨光劉麗華溫世豪翟光濤
    重慶理工大學學報(自然科學) 2023年10期
    關鍵詞:氣動力賽車氣動

    賴晨光,劉麗華,馮 帥,溫世豪,翟光濤

    (重慶理工大學 車輛工程學院,重慶 400054)

    0 引言

    大學生方程式比賽簡稱為FSAE(Formula SAE),于1978年由美國汽車工程學會(SAE)舉辦,是一項由學生自主設計、制造與測試一輛小型單座開輪式賽車的比賽。合理的氣動設計可以為賽車產生下壓力,以提高其在彎道中的速度。因此,氣動設計是FSAE賽車的重要環(huán)節(jié)之一。在動態(tài)賽中,彎道眾多,賽車大部分時間均處于彎道工況。賽車在彎道中氣動力與氣動力矩均會發(fā)生改變,從而影響賽車在極限條件下的操縱穩(wěn)定性。因此,在FSAE賽車氣動設計中僅考慮直道工況較為不妥。

    從20世紀60年代開始[1],空氣動力學套件就被應用于賽車以提高操縱穩(wěn)定性。FSAE賽車與一級方程式賽車均屬于開輪式賽車,因此具有較為相似的氣動布局。FSAE賽車平均時速約60 km/h,但氣動套件對FSAE賽車性能具有很大提升,這一點已經得到了Wordley等的驗證[2]。他們研究了FSAE賽車在直線加速、八字繞環(huán)(skid pad)、高速避障和耐久賽等不同動態(tài)項目中加裝倒置翼片的效果,表明加裝了倒置翼片的FSAE賽車在直線加速所用時間更長,但在其他動態(tài)賽中與沒有安裝倒置翼片相比具有很大優(yōu)勢。

    有學者利用計算流體力學(CFD)方法對賽車在彎道處的空氣動力學特性進行了研究。Douglas等[3]考慮了橫擺角對賽車氣動性能的影響,發(fā)現偏航時翼片端板對下壓力有不良影響,特別是橫擺角較大時。Watanabe等[4]提出了一種新的方法——移動計算域(moving computational domain,MCD)來模擬高速汽車通過發(fā)夾彎的情況。Tsubokura等[5]對車輛動態(tài)轉向時的氣動力進行了研究,指出高速轉彎時產生的不穩(wěn)定氣動力會影響汽車的穩(wěn)定性。Nara等[6]用移動計算域的方法對方程式賽車進行了彎道模擬,模擬了賽車在富士賽道上行駛,并解釋了在滑移角和速度變化上存在氣動力差異的原因。Keogh等[7]提出了分析彎道中車輛空氣動力學方法,數值模擬中需要模擬車輛周圍呈現曲率運動的流場情況,并指出沿彎道徑向的流速會呈梯度變化。還通過數值方法研究了Ahmed模型在彎道中的空氣動力學,觀察到偏航力矩引起的氣動力和不對稱尾跡的變化[8]。Sagmo等[9]對FSAE賽車的彎道空氣動力學進行模擬,研究了彎道工況與橫擺工況相比氣動力變化,發(fā)現穩(wěn)態(tài)彎道工況與橫擺角之間的流場存在差異。Patel等[10]研究了多元倒置翼片的彎道和地面效應時的氣動特性,發(fā)現了渦軌跡變化。Nakashima等[11]對Ahmed模型的彎道氣動特性進行了研究,發(fā)現了氣動阻力增加和氣動力矩的產生。Piechna等[12]通過數值模擬和道路實驗研究了主動式氣動部件在轉彎時的性能。宋世達[13]通過仿真研究了不同工況下的氣動性能差異,探究了各氣動裝置對整車性能的影響。梁德慧等[14]利用場函數方法實現了對FSAE賽車穩(wěn)態(tài)彎道的仿真。

    以往對于彎道空氣動力學研究大多集中在過彎模擬,缺乏對彎道過程中流場分析。近年來,賽車空氣動力學設計日趨復雜,越來越多工況被納入考慮范圍之內,因此,需較為深入的研究賽車在彎道處流場變化細節(jié),以更好地理解彎道工況中氣動力變化。筆者在前人研究的基礎上,通過數值仿真研究FSAE賽車在彎道工況中流場變化導致氣動力變化的原因。重疊網格被證明是一種較好的模擬運動物體的方法,目前應用廣泛[15-19],因此,本研究選擇重疊網格來模擬FSAE賽車勻速過彎工況。本研究可為開輪式賽車的空氣動力學設計提供指導意義。

    1 數值仿真

    1.1 空氣動力學模型

    為保證本研究普適性,建立一簡化FSAE賽車模型。如圖1所示,它具有FSAE賽車大部分特征,并符合FSAE賽事規(guī)則,在該模型中忽略懸掛桿系,并簡化后艙與座艙等區(qū)域。該模型主要由前翼、車身、尾翼、車輪和擴散器組成。模型基本幾何參數如表1所示。

    表1 FSAE模型基本幾何參數

    圖1 FSAE賽車模型

    1.2 直道與彎道工況邊界條件

    為保證仿真一致性,在直道和彎道工況下均采取重疊網格方法。直道工況計算域尺寸與邊界如圖2所示。直道工況計算域寬約模型寬11倍,高約9倍,阻塞比約0.63%。由于擬模擬時間為2 s,計算域長度應大于2 s的運動距離,考慮空氣黏度對邊界的影響,故計算域長度為36 m。計算域兩端設有壓力出口,允許氣流進出,其他邊界為壁面。如圖3所示為直道工況重疊網格邊界,重疊網格由背景域、重疊域和交接面組成,FSAE賽車運動速度為15 m/s。

    圖2 直道工況計算域尺寸與邊界

    圖3 直道工況重疊網格邊界

    彎道工況計算域如圖4所示,寬度和高度約為模型寬度的8倍,高度的7倍,阻塞比約1.2%。計算域內半徑為5 m,計算域外半徑為15 m,因此轉彎半徑為10 m。邊界條件與直道保持一致。圖5為彎道工況下重疊網格邊界,背景區(qū)域為圓形弧線,運動速度為15 m/s,軌跡沿虛線所示半徑為 10 m的曲線運動。

    圖4 彎道工況計算域

    圖5 彎道工況重疊網格邊界

    1.3 車身姿態(tài)

    前輪轉角是一個不可忽視的因素。在忽略車輪側偏角情況下,轉彎時前輪轉角為:

    δ=arcsin(L/R)

    (1)

    式中:L為軸距;R為彎道半徑。算出無側偏角下前輪轉角為8.9°。

    賽車在極限條件下具有較大的輪胎側偏角,導致前輪需要更大轉向角[20]。根據FSAE輪胎的特性和彎道條件[21],側偏角取值為4°。忽略阿克曼角,因此總轉向角為12.9°。默認賽車處于中性轉向,忽略俯仰角和橫擺角,側傾角為2°。因此FSAE賽車在彎道中的姿態(tài)如圖6所示。

    圖6 彎道工況中的姿態(tài)

    1.4 網格

    此模型幾何較為復雜,采用非結構化網格可以更好地貼合表面。利用ANSA軟件進行幾何處理和曲面網格劃分。為平衡計算精度與速度,將曲面劃分不同面網格尺寸,如表2所示。為捕捉翼片邊界層流動細節(jié),翼片表面具有最小的表面網格尺寸。直線工況和彎道工況具有相同的表面網格尺寸以保持一致性。

    表2 面網格尺寸

    由于SSTk-ω湍流模型對壁面要求,需要更多網格來捕捉壁面附近邊界層的流動。因此,壁面y+值約為1。棱柱層參數如表3所示。

    表3 棱柱層參數

    網格數量對計算結果影響較大,有必要驗證網格無關性。下壓力作為本文中主要監(jiān)測數值,其準確性至關重要。網格獨立性驗證如表4所示,復雜流動區(qū)和尾跡區(qū)對結果有重要影響,因此設置2個網格加密區(qū),其一對前翼、擴散器和尾翼加密,其二對賽車尾流加密。從表4可以看出,隨著體網格數增加,下壓力急劇下降,當網格數達到2 000萬后,下壓力略有變化。因此,考慮到計算效率,確定流動復雜區(qū)和尾跡區(qū)網格尺寸分別為8 mm和16 mm,總網格數量控制在2 000萬~3 000萬。

    表4 網格獨立性驗證

    具體體網格數如表5所示。由于運動軌跡不同,直線和彎道背景域存在一定差異。2種工況重疊域網格數均保持在2 000萬左右。

    表5 體網格數

    1.5 湍流模型

    本研究主要應用Reynold-averaged Navier-Stokes(RANS)進行數值模擬,目前得到普遍應用[22-24]。RANS湍流模型包括Realizablek-ε、Standardk-ε、SSTk-ω等。SSTk-ω湍流模型被證明是一種能較好地預測賽車下壓力的模型[25],并被廣泛應用[26-27]。SSTk-ω結合了k-ω模型在壁面附近的準確性和k-ε在自由流中的魯棒性,其考慮了湍流剪應力的傳輸,可以很好地預測逆壓梯度下的流動起始和分離量[28]。

    2 結果與分析

    2.1 氣動力結果

    圖7為總氣動力結果。從圖7可以看出,與直道相比,FSAE賽車在彎道中總氣動力出現較明顯變化。彎道工況下阻力和下壓力有不同程度的增加,氣動阻力為188.2 N,與直道相比增加1.5%,下壓力為399.2 N,增加6.5%。表明,彎道工況中下壓力比阻力更敏感。

    圖7 總氣動力結果

    圖8為各部件的氣動力結果。圖8(a)為下壓力變化,除前翼外,其余部件在彎道時下壓力均有不同幅度地增加,增幅主要來自擴散器和車身,分別增加11%和18.3%,尾翼下壓力增加2.5%,前翼下壓力減少約3%。圖8(b)顯示各部件氣動阻力變化,變化最大的為前翼阻力的增加,與直道相比增加約14.8%。擴散器和車身阻力分別增加約11%和2.6%,尾翼和車輪阻力分別減少約1%和7.4%。下壓力變化對彎道中賽車的性能影響更大。

    圖8 各部件氣動力結果

    2.2 氣動力矩結果

    除了氣動力外,氣動力矩對賽車操縱穩(wěn)定性有較大影響[29]。氣動載荷變化與總氣動力矩如表6所示,%F表示前軸氣動載荷與總氣動載荷之比,Mx和Mz分別為氣動側傾力矩和橫擺力矩。氣動力矩均以圖6所示坐標系為基準,彎道工況下%F為30.8%,與直道相比下降1.1%,彎道工況下Mx和Mz分別為-13.5 N·m和-43.1 N·m。

    表6 彎道中氣動力矩的變化

    為直觀理解氣動力矩對操縱穩(wěn)定性的影響,彎道工況中氣動力矩如圖9所示,%F的降低會減少前軸載荷與前輪抓地力,增加彎道中的轉向不足趨勢。負側傾力矩會減小賽車車身左右兩側的氣動載荷之差。負橫擺力矩會抑制賽車在彎道中的入彎趨勢,進而導致轉向不足。綜合以上分析,彎中產生的氣動力矩會增加賽車的轉向不足,導致賽車在過彎時轉向更為遲鈍,對整車操縱穩(wěn)定性較為不利。

    圖9 彎道工況中氣動力矩的影響

    2.3 流場分析

    根據上述直道和彎道工況中氣動力和力矩值,兩截面的速度場分布對比如圖10所示,速度參考系選取為運動中的賽車。本文中主要關注前翼與尾翼端板附近流場,因此,在縱向上創(chuàng)建了Z=-200 mm和Z=700 mm的截面,如圖10(c)所示。圖10(a)為Z=-200 mm處的速度分布對比,如箭頭所示,前翼端板在彎道外側有擾亂高速氣流現象,同時前輪尾跡存在不同,在彎道外側有更多低速區(qū)。圖10(b)為Z=700 mm的速度分布對比。從圖中可觀察到尾翼附近流場,尾翼后方有一個速度較低的尾流區(qū)域(如圖中虛線所示),彎道工況下尾流區(qū)域與直道相比有明顯的氣流扭轉現象。如箭頭所示,尾翼前端內側端板出現氣流分離現象。

    圖10 速度分布對比

    壓力系數分布如圖11所示。圖11(a)為上表面,從中心軸左右對比來看,彎道中表面壓力呈現較大不對稱分布,如箭頭所示,大部分正壓力積聚在翼片的彎道外部一側(包括前翼和后翼),這會造成一定的側傾力矩,加大左右兩側氣動載荷分布。底板表面壓力分布如圖11(b)所示,與正壓分布相比,負壓分布不對稱趨勢更強。圖中箭頭所示,較多的負壓分布在彎道外側,包括前翼和擴散器,原因其一為外部氣流流速較快,其二為彎道外側車身側傾角導致離地間隙減小,增加了擴散器地面效應。

    圖11 壓力系數分布

    為驗證側傾角引起底板負壓分布變化這一假設,對無側傾角工況進行模擬,結果如圖12所示。側傾角存在對底板負壓分布有重要影響,有側傾角下底板負壓較大,產生了更多下壓力,因此側傾角存在是引起彎道中擴散器和底板下壓力增加的主要原因。

    圖12 有無側傾角底板表面壓力對比

    本研究采用Q準則來識別渦結構以及捕捉渦軌跡。Hunt等提出將Q>0的區(qū)域定義為漩渦[30],這意味著在渦區(qū)域內,流體旋轉渦量大于變形渦量。Q準則是一個基于速度梯度張量的特征向量和特征函數,需要設定閾值來顯示渦旋結構。

    圖13是FSAE賽車主要渦結構變化,Q值選擇為5 000 s-2以篩選出明顯渦結構。如圖13(a)所示,在直道工況中前翼兩側產生一對翼尖渦,命名為FW-TV(front wing-tip vortex),尾翼產生的一對翼尖渦命名為RW-TV(rear wing-tip vortex)。圖13(b)顯示彎道工況下渦的運動軌跡,其渦名稱尾部字母L和R分別代表車身左右側渦,同時也代表彎道內外側的渦。在直線工況下,左右兩側渦的形態(tài)和軌跡均呈對稱分布,且FW-TV和RW-TV具有較長的縱向尾跡。彎道中2對渦軌跡沿一定曲率向下游發(fā)展,車身左右兩側渦對的發(fā)展和軌跡呈不對稱分布,位于彎道內側的FW-TV-L渦與具有轉角的左前輪撞擊而破碎,阻止了其向下游運動,這是導致前翼下壓力損失的主要原因。位于彎道外側的FW-TV-R渦尾跡區(qū)相比于直道更長,這是由于來流的曲率影響使FW-TV-R渦沿車身外側流動,對于彎道工況中的RW-TV渦,則體現出不同趨勢,位于彎道內側的RW-TV-L渦縱向尾渦尾跡區(qū)比彎道外側RW-TV-R的長,這種差異體現出彎道內外側不同的流動特性。

    圖13 整車主要渦結構變化

    3 八字繞環(huán)試驗

    八字繞環(huán)(skid pad)是一種特殊的彎道工況,它可以測試FSAE賽車的操穩(wěn)性。以2021賽季重慶理工大學FSCC賽車為基礎,分別進行模擬仿真和八字繞環(huán)試驗。

    3.1 CFD模擬

    八字繞環(huán)(Skid Pad)是FSAE賽車動態(tài)賽事之一,旨在測試賽車在固定半徑轉彎時在平地上的轉向能力。八字繞環(huán)布局如圖14所示,兩圓心之間的距離為18.25 m,內圓直徑為15.25 m,外圓直徑為21.25 m,賽道寬度為3 m,結合車寬,汽車質心的行駛半徑約8.5 m。

    圖14 八字繞環(huán)布局

    實驗仿真模型如圖15所示,實驗模型采用21賽季重慶理工大學Thimchi 39賽車,實車模型如圖15(a)所示。為驗證八字繞環(huán)實驗,需建立與試驗模型相對應的仿真模型,如圖15(b)所示,該模型保留了實車大部分特征,對翼片和擴散器進行精細化處理,簡化了輪胎、發(fā)動機和減震器,去掉外露的軟管、螺栓等小部件。

    圖15 試驗所用模型

    為保證與八字繞環(huán)試驗條件一致,在數值模擬中,轉向半徑8.5 m,在八字繞環(huán)中環(huán)繞一圈時間為6 s。其余邊界條件保持相同。

    仿真結果如圖16所示,包括氣動力和力矩。圖16(a)為總氣動力,總氣動阻力與下壓力分別為77.3 N和123.9 N,各部件氣動阻力和下壓力分別如圖16(b)和圖16(c)所示,氣動力矩如圖16(d)所示,FSAE賽車在八字繞環(huán)時會產生負氣動力矩。氣動力與氣動力矩變化趨勢與簡化FSAE模型相同,驗證了本研究的普適性。

    圖16 仿真結果

    3.2 試驗過程

    安裝在車身上的線位移傳感器用于測量總下壓力。線位移傳感器的安裝位置如圖17所示,如圖中箭頭所示,線位移傳感器安裝在每個減振器兩端,并伴隨減振器做相同反彈和壓縮運動。在彈簧線性剛度條件下,通過測量減振器變化量來計算整車下壓力。

    圖17 線位移傳感器的安裝位置

    試驗地點在重慶大足汽車試驗場,如圖18所示。試驗場的平整路面減少了慣性力對結果影響,使用樁桶鋪設八字繞環(huán)賽道,繞行一圈時間控制在6 s左右,車手盡量保持勻速繞樁。

    圖18 重慶大足試驗場

    實驗數據處理過程如圖19所示,圖19(a)為數據采集的過程,線位移傳感器將采集的數據轉換為電壓信號并存儲在電子控制單元(ECU)中,用外部電腦提取存儲的數據,對數據進行處理,最后轉換成下壓力值。圖19(b)為線位移傳感器采集的初始數據,使用Motec軟件讀取ECU中數據,包括油門開度、制動踏板開度、方向盤轉角、車速等隨時間變化的數據。本試驗主要關注位移變化如方框所示。如圖19 (c)所示,從初始數據中提取0.5 s波動較小的數據,采樣頻率為500 Hz,因此一個傳感器在此時間段內采集251個數據,隨后對這段數據進行統(tǒng)計和處理,處理結果如圖19(d)所示,FL、FR、RL和RR分別代表前左、前右、后左和后右減振器位移量。

    圖19 數據處理過程

    駕駛員坐入賽車靜止不動時,將傳感器采集到的數據作為靜態(tài)數據,在既定工況駕駛后,對動態(tài)數據中0.5 s內取平均,如表7所示為數據統(tǒng)計,正值(+)表示減震器被拉伸,負值(-)表示減震器被壓縮。

    表7 數據統(tǒng)計

    線位移數據變化如圖20所示。圖20(a)為各個減振器的變化。從圖中可以看到彎道外側減振器被壓縮,彎道內側減振器被拉伸。將兩側的數值加和為前后軸荷的變化,如圖20(b)所示,前軸被壓縮2.68 mm,后軸被壓縮0.8 mm,這是由于存在的氣動下壓力導致車身壓向路面,前后減振器彈簧均被壓縮。

    圖20 線位移變化量可視化

    根據胡克定律[31],減振器默認為線性彈簧,忽略第三彈簧與防傾桿剛度的影響,汽車前后減振器彈簧型號為200磅,根據剛度折算為30 250 N/m。

    DFF=-0.002 68mm×35 025 N/m=-93.9 N

    DFR=-0.000 8 mm×35 025 N/m=-28 N

    DFTotal=DFF+DFR=-121.9 N

    式中:DFF為前軸轉換后的下壓力;DFR為后軸的下壓力。DFF和DFR的值不代表真實值,因為從壓縮量來看,前軸的壓縮量明顯大于后軸,說明有載荷轉移存在,但DFTotal沒有變化。在八字繞環(huán)中,整車氣動下壓力為121.9 N。

    3.3 結果對比

    實驗與仿真的下壓力對比如圖21所示,試驗值為121.9 N,略低于仿真值,誤差約1.6%,造成誤差的原因之一是輪胎變形,實際充氣輪胎會比仿真模型大,另一個原因是散熱器處理,忽略了旋轉冷卻風扇的影響。

    圖21 試驗與仿真下壓力對比

    試驗驗證了重疊網格方法模擬復雜工況的準確性。由于試驗中可測量數據有限,需要更多仿真的數據,因此,CFD仿真結果可用于指導賽車設計和優(yōu)化底盤調整。

    4 結論

    1) FSAE賽車在彎道中氣動力發(fā)生較大變化,其中下壓力增大6.5%,除前翼外,其他部件均有不同程度的增幅。下壓力增加主要來自擴散器和底板。

    2) 彎道中會產生額外氣動力矩,俯仰力矩會使壓力中心向后偏移,橫擺力矩會降低FSAE賽車的轉向靈敏度。由于氣動力矩的存在,賽車會趨于轉向不足,車手了解賽車氣動力和力矩的變化,更好地駕馭賽車。

    3) 彎道中車身表面靜壓分布高度不對稱,尤其是底板處的負壓,這是由于轉彎時存在的側傾角導致擴散器和車身下壓力增大。

    4) 渦結構在彎道中有較大變化,具有轉角的前輪對彎道內側前翼渦運動的阻滯效應是前翼下壓力降低的主要原因,從渦軌跡可以看出,尾翼旋渦的軌跡發(fā)生了一定程度的偏轉,尾翼端板氣流在彎道中發(fā)生氣流分離。

    5) 利用線位移傳感器監(jiān)測和計算FSAE賽車在八字繞環(huán)中的氣動下壓力。結果表明,實測下壓力值與數值模擬相差1.6%,吻合度較好。

    在本研究中,有一些現象值得進一步研究。首先,倒置翼片是賽車中產生下壓力的主要方式之一,在本研究中,彎道下翼片產生的下壓力沒有體現顯著差異,因此需要進一步研究曲率流動下機翼的氣動特性,如升阻性能、臨界攻角變化、壁面附近邊界層流動等。此外,如何減小曲率流動中端板的分離是后續(xù)研究的重點。

    猜你喜歡
    氣動力賽車氣動
    賽車
    幼兒畫刊(2023年11期)2024-01-15 10:00:32
    中寰氣動執(zhí)行機構
    基于NACA0030的波紋狀翼型氣動特性探索
    飛行載荷外部氣動力的二次規(guī)劃等效映射方法
    云上賽車
    基于反饋線性化的RLV氣動控制一體化設計
    測控技術(2018年9期)2018-11-25 07:44:24
    加油吧,賽車!
    側風對拍動翅氣動力的影響
    我家的小賽車
    高速鐵路接觸線覆冰后氣動力特性的風洞試驗研究
    黄色视频,在线免费观看| 欧美精品一区二区免费开放| 亚洲精品成人av观看孕妇| 日韩有码中文字幕| 18禁观看日本| 黄色片一级片一级黄色片| 国产成+人综合+亚洲专区| 亚洲熟妇熟女久久| 午夜久久久在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产精品 国内视频| 国产精品99久久99久久久不卡| 免费观看人在逋| 淫妇啪啪啪对白视频| 99re在线观看精品视频| 午夜精品在线福利| 久久国产精品影院| 国产亚洲精品久久久久5区| 脱女人内裤的视频| 高清视频免费观看一区二区| 一区二区三区激情视频| 国产精品免费视频内射| 国产野战对白在线观看| 久久久精品免费免费高清| 精品国产亚洲在线| 久久久久久久国产电影| 高清视频免费观看一区二区| 母亲3免费完整高清在线观看| 欧美乱妇无乱码| 国产亚洲欧美98| 日韩欧美在线二视频 | 亚洲少妇的诱惑av| 亚洲男人天堂网一区| 黑丝袜美女国产一区| 亚洲国产精品一区二区三区在线| 亚洲人成电影观看| 女人被狂操c到高潮| 高清在线国产一区| av不卡在线播放| 在线观看www视频免费| 日韩欧美国产一区二区入口| 欧美老熟妇乱子伦牲交| 中文字幕另类日韩欧美亚洲嫩草| 69精品国产乱码久久久| 国产人伦9x9x在线观看| 成年版毛片免费区| 亚洲伊人色综图| 老汉色∧v一级毛片| 男女午夜视频在线观看| 天堂动漫精品| 久久精品91无色码中文字幕| 十分钟在线观看高清视频www| 免费在线观看日本一区| 美女国产高潮福利片在线看| 亚洲专区中文字幕在线| videos熟女内射| 久久中文看片网| 99国产综合亚洲精品| 黄片大片在线免费观看| 亚洲成国产人片在线观看| 国产99久久九九免费精品| 又黄又爽又免费观看的视频| 飞空精品影院首页| 99re在线观看精品视频| 国产av精品麻豆| 这个男人来自地球电影免费观看| 首页视频小说图片口味搜索| 国产在线精品亚洲第一网站| 水蜜桃什么品种好| 免费看a级黄色片| 99精品欧美一区二区三区四区| 亚洲七黄色美女视频| 欧美日韩视频精品一区| 久久热在线av| 在线观看66精品国产| 91成人精品电影| 91麻豆精品激情在线观看国产 | 一个人免费在线观看的高清视频| 亚洲熟女精品中文字幕| 一级片'在线观看视频| 成人18禁高潮啪啪吃奶动态图| 9191精品国产免费久久| 身体一侧抽搐| 在线免费观看的www视频| 国产伦人伦偷精品视频| 操美女的视频在线观看| 人人妻人人添人人爽欧美一区卜| 日本黄色视频三级网站网址 | 亚洲欧美一区二区三区黑人| 女人被狂操c到高潮| 国精品久久久久久国模美| 少妇 在线观看| 国产成人系列免费观看| 成人亚洲精品一区在线观看| 中文字幕人妻熟女乱码| 国产男女内射视频| 欧美精品啪啪一区二区三区| 欧美精品高潮呻吟av久久| 99精品欧美一区二区三区四区| 久久久久久久精品吃奶| 免费在线观看黄色视频的| 黄色丝袜av网址大全| 成人永久免费在线观看视频| 久久久国产欧美日韩av| 99精国产麻豆久久婷婷| 中文欧美无线码| 国产亚洲欧美精品永久| 91精品国产国语对白视频| 国产成+人综合+亚洲专区| 久久久久久免费高清国产稀缺| 亚洲成人免费av在线播放| 免费在线观看完整版高清| 亚洲国产精品一区二区三区在线| 欧美成人免费av一区二区三区 | 国产激情欧美一区二区| 久久 成人 亚洲| 免费在线观看视频国产中文字幕亚洲| 国产精品影院久久| 99在线人妻在线中文字幕 | 一本综合久久免费| 999久久久国产精品视频| 91麻豆精品激情在线观看国产 | 母亲3免费完整高清在线观看| 在线观看免费午夜福利视频| 欧美大码av| 亚洲午夜精品一区,二区,三区| 国产在线一区二区三区精| 午夜福利乱码中文字幕| 国产不卡av网站在线观看| 国产免费现黄频在线看| 91国产中文字幕| 一级毛片女人18水好多| 国产成人欧美| 亚洲av熟女| 国产高清国产精品国产三级| 欧美国产精品va在线观看不卡| av欧美777| 十八禁网站免费在线| 男人操女人黄网站| 欧美大码av| 欧美 日韩 精品 国产| 欧美亚洲日本最大视频资源| 深夜精品福利| 久久午夜亚洲精品久久| 欧美日韩国产mv在线观看视频| 动漫黄色视频在线观看| 欧美国产精品va在线观看不卡| 亚洲国产欧美日韩在线播放| 亚洲男人天堂网一区| 19禁男女啪啪无遮挡网站| 久久天堂一区二区三区四区| 日韩欧美三级三区| 国产亚洲精品久久久久5区| 亚洲精品久久成人aⅴ小说| 视频在线观看一区二区三区| 亚洲avbb在线观看| 一二三四在线观看免费中文在| 在线观看一区二区三区激情| 国产有黄有色有爽视频| 日本撒尿小便嘘嘘汇集6| 中文字幕av电影在线播放| 下体分泌物呈黄色| 亚洲中文av在线| 午夜福利在线观看吧| 亚洲avbb在线观看| 亚洲色图av天堂| 搡老岳熟女国产| 久久中文字幕一级| 亚洲一区高清亚洲精品| 免费日韩欧美在线观看| 天天躁日日躁夜夜躁夜夜| 自线自在国产av| 国产不卡av网站在线观看| 狠狠狠狠99中文字幕| 伦理电影免费视频| 日本黄色视频三级网站网址 | av福利片在线| 精品高清国产在线一区| 老熟女久久久| 五月开心婷婷网| 国产单亲对白刺激| 女性被躁到高潮视频| 国产在线观看jvid| 亚洲av第一区精品v没综合| 操美女的视频在线观看| 久久精品亚洲精品国产色婷小说| 一二三四社区在线视频社区8| 两性夫妻黄色片| 夫妻午夜视频| 97人妻天天添夜夜摸| 精品一区二区三卡| 亚洲一码二码三码区别大吗| 国产精品电影一区二区三区 | 美女高潮喷水抽搐中文字幕| 夜夜爽天天搞| 少妇猛男粗大的猛烈进出视频| 免费在线观看影片大全网站| 男女床上黄色一级片免费看| 免费不卡黄色视频| 亚洲va日本ⅴa欧美va伊人久久| 99久久人妻综合| 91麻豆精品激情在线观看国产 | 国产av又大| 正在播放国产对白刺激| 好男人电影高清在线观看| 国产亚洲欧美在线一区二区| 99热网站在线观看| 亚洲专区字幕在线| 嫩草影视91久久| 99在线人妻在线中文字幕 | 亚洲欧美日韩另类电影网站| 人人妻人人爽人人添夜夜欢视频| 一进一出抽搐动态| 国产精品一区二区在线不卡| 午夜两性在线视频| 成人亚洲精品一区在线观看| 久久精品国产99精品国产亚洲性色 | 欧美大码av| 老司机靠b影院| 捣出白浆h1v1| 欧美丝袜亚洲另类 | 麻豆乱淫一区二区| 国产一卡二卡三卡精品| 青草久久国产| 亚洲 欧美一区二区三区| 一进一出抽搐gif免费好疼 | 他把我摸到了高潮在线观看| 久久精品亚洲av国产电影网| 国产在线精品亚洲第一网站| 女人被躁到高潮嗷嗷叫费观| 成人免费观看视频高清| 极品人妻少妇av视频| 搡老熟女国产l中国老女人| 国产一区有黄有色的免费视频| 久久久久久久精品吃奶| 国产人伦9x9x在线观看| 999久久久国产精品视频| 露出奶头的视频| 两人在一起打扑克的视频| 亚洲一区二区三区不卡视频| 女人高潮潮喷娇喘18禁视频| 交换朋友夫妻互换小说| 久久久久精品国产欧美久久久| 99久久国产精品久久久| 岛国在线观看网站| 国产主播在线观看一区二区| 最近最新免费中文字幕在线| 国产人伦9x9x在线观看| 好男人电影高清在线观看| 国内毛片毛片毛片毛片毛片| 亚洲欧洲精品一区二区精品久久久| 国产精品亚洲一级av第二区| 国产欧美日韩精品亚洲av| 欧美日韩一级在线毛片| 色综合婷婷激情| 男女午夜视频在线观看| av超薄肉色丝袜交足视频| 9色porny在线观看| 国产野战对白在线观看| 夜夜爽天天搞| 亚洲在线自拍视频| 亚洲伊人色综图| 精品国产乱码久久久久久男人| 12—13女人毛片做爰片一| 成年版毛片免费区| 宅男免费午夜| 国产极品粉嫩免费观看在线| 黑人操中国人逼视频| 国产男靠女视频免费网站| 香蕉丝袜av| 精品久久蜜臀av无| 在线观看66精品国产| 日韩精品免费视频一区二区三区| a级毛片在线看网站| 色尼玛亚洲综合影院| 日本黄色视频三级网站网址 | av片东京热男人的天堂| 一本综合久久免费| 精品人妻在线不人妻| 最新的欧美精品一区二区| 啦啦啦在线免费观看视频4| 成人国产一区最新在线观看| 777久久人妻少妇嫩草av网站| 韩国精品一区二区三区| 国产精品久久久久久人妻精品电影| 高清在线国产一区| 一进一出抽搐gif免费好疼 | 亚洲av成人av| 亚洲人成伊人成综合网2020| 久久精品亚洲精品国产色婷小说| 国产精品国产av在线观看| 中文字幕人妻熟女乱码| 黑人操中国人逼视频| 亚洲国产毛片av蜜桃av| 亚洲精华国产精华精| 国产97色在线日韩免费| 国产99白浆流出| 国产三级黄色录像| av片东京热男人的天堂| 淫妇啪啪啪对白视频| 黑人巨大精品欧美一区二区蜜桃| 久久中文字幕人妻熟女| 久久国产精品人妻蜜桃| 天天躁夜夜躁狠狠躁躁| 欧美日韩乱码在线| 成人18禁在线播放| 身体一侧抽搐| 大陆偷拍与自拍| 免费久久久久久久精品成人欧美视频| 亚洲熟女毛片儿| 久久精品熟女亚洲av麻豆精品| 欧美精品一区二区免费开放| 国产在视频线精品| 一a级毛片在线观看| 超色免费av| 99久久99久久久精品蜜桃| 国产午夜精品久久久久久| 十分钟在线观看高清视频www| 三上悠亚av全集在线观看| 亚洲七黄色美女视频| 久久人人爽av亚洲精品天堂| 搡老岳熟女国产| 中国美女看黄片| 亚洲色图av天堂| 成年人免费黄色播放视频| 国产单亲对白刺激| 悠悠久久av| 久久99一区二区三区| xxxhd国产人妻xxx| 亚洲第一av免费看| 中文字幕人妻熟女乱码| 亚洲av片天天在线观看| 亚洲熟妇熟女久久| 大片电影免费在线观看免费| 男人舔女人的私密视频| 一级,二级,三级黄色视频| 黄色怎么调成土黄色| 在线观看免费视频网站a站| 女警被强在线播放| 欧美激情久久久久久爽电影 | 国产激情久久老熟女| 建设人人有责人人尽责人人享有的| 人人妻人人澡人人爽人人夜夜| 欧洲精品卡2卡3卡4卡5卡区| av中文乱码字幕在线| 久久精品成人免费网站| 99久久99久久久精品蜜桃| 热99久久久久精品小说推荐| 国产区一区二久久| 精品国产亚洲在线| 99国产精品一区二区蜜桃av | 精品少妇一区二区三区视频日本电影| 90打野战视频偷拍视频| 9色porny在线观看| 波多野结衣av一区二区av| 两性夫妻黄色片| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲七黄色美女视频| 中文字幕高清在线视频| 亚洲精品中文字幕在线视频| 国产主播在线观看一区二区| 亚洲黑人精品在线| 少妇猛男粗大的猛烈进出视频| 国产淫语在线视频| 叶爱在线成人免费视频播放| 操出白浆在线播放| av天堂在线播放| 中国美女看黄片| 精品少妇一区二区三区视频日本电影| 天堂中文最新版在线下载| 女人爽到高潮嗷嗷叫在线视频| 精品久久久久久久毛片微露脸| 国产人伦9x9x在线观看| 国产精品二区激情视频| 精品少妇久久久久久888优播| 亚洲片人在线观看| 午夜成年电影在线免费观看| 国产淫语在线视频| 久久久久久亚洲精品国产蜜桃av| 在线视频色国产色| 可以免费在线观看a视频的电影网站| 国产成人欧美| 嫁个100分男人电影在线观看| 久久精品国产综合久久久| 女性生殖器流出的白浆| 王馨瑶露胸无遮挡在线观看| 精品人妻熟女毛片av久久网站| 人妻久久中文字幕网| 亚洲精品国产区一区二| 免费人成视频x8x8入口观看| 人妻丰满熟妇av一区二区三区 | 巨乳人妻的诱惑在线观看| 欧美乱妇无乱码| 久久精品亚洲熟妇少妇任你| 午夜福利,免费看| 777米奇影视久久| 在线观看66精品国产| 成人免费观看视频高清| 国产精品一区二区免费欧美| 日韩欧美国产一区二区入口| 国产精华一区二区三区| 国产1区2区3区精品| 另类亚洲欧美激情| 女人精品久久久久毛片| 黑人巨大精品欧美一区二区mp4| 99精国产麻豆久久婷婷| 久久久久国产一级毛片高清牌| 少妇粗大呻吟视频| a级毛片在线看网站| 国产成+人综合+亚洲专区| 巨乳人妻的诱惑在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲精品国产色婷婷电影| 美女视频免费永久观看网站| 久久中文字幕一级| 怎么达到女性高潮| 国产欧美日韩一区二区三| 欧美亚洲 丝袜 人妻 在线| 人妻丰满熟妇av一区二区三区 | 亚洲欧美激情在线| 欧美日韩亚洲综合一区二区三区_| 久久久久久久久免费视频了| 欧美 日韩 精品 国产| 欧美精品人与动牲交sv欧美| 黑人操中国人逼视频| 美女高潮喷水抽搐中文字幕| 美女扒开内裤让男人捅视频| 99久久人妻综合| 国产色视频综合| 怎么达到女性高潮| 亚洲片人在线观看| 久久久久视频综合| 亚洲美女黄片视频| 国产精品秋霞免费鲁丝片| 久久久久国产精品人妻aⅴ院 | 亚洲精品美女久久av网站| 亚洲综合色网址| 欧美不卡视频在线免费观看 | 高清av免费在线| 亚洲精品中文字幕在线视频| 亚洲人成电影观看| 精品一区二区三区av网在线观看| 人妻 亚洲 视频| 精品久久蜜臀av无| 黑人巨大精品欧美一区二区蜜桃| 午夜精品在线福利| 国产视频一区二区在线看| 日韩欧美免费精品| 一级作爱视频免费观看| 日韩中文字幕欧美一区二区| 午夜亚洲福利在线播放| 亚洲欧美一区二区三区久久| 久久久久精品人妻al黑| 夜夜爽天天搞| 在线观看舔阴道视频| 亚洲精品美女久久av网站| 国产精品国产高清国产av | 91成人精品电影| 成人黄色视频免费在线看| 精品国产乱子伦一区二区三区| 欧美 日韩 精品 国产| 免费高清在线观看日韩| 免费观看精品视频网站| 天天躁狠狠躁夜夜躁狠狠躁| 人人妻人人澡人人看| 亚洲中文字幕日韩| 这个男人来自地球电影免费观看| 久久久久国产精品人妻aⅴ院 | 天堂√8在线中文| 国产激情久久老熟女| 天天影视国产精品| 视频区欧美日本亚洲| 91老司机精品| 一进一出好大好爽视频| 老司机福利观看| 午夜影院日韩av| 国内毛片毛片毛片毛片毛片| 日日摸夜夜添夜夜添小说| 啦啦啦免费观看视频1| 99riav亚洲国产免费| av免费在线观看网站| 亚洲九九香蕉| 99国产精品99久久久久| 国产免费av片在线观看野外av| 在线观看免费午夜福利视频| 亚洲精品中文字幕在线视频| 在线观看午夜福利视频| 亚洲午夜精品一区,二区,三区| 黄频高清免费视频| 淫妇啪啪啪对白视频| 伊人久久大香线蕉亚洲五| 亚洲精华国产精华精| 久久久久久久久免费视频了| a级片在线免费高清观看视频| 午夜老司机福利片| ponron亚洲| 欧美日韩视频精品一区| 成人特级黄色片久久久久久久| 最近最新中文字幕大全免费视频| 欧美乱码精品一区二区三区| 国产精品影院久久| 免费在线观看影片大全网站| 国产精品av久久久久免费| 美女福利国产在线| 亚洲精品美女久久av网站| 99久久精品国产亚洲精品| ponron亚洲| 99精品欧美一区二区三区四区| 国产成人一区二区三区免费视频网站| 午夜两性在线视频| 久久人妻熟女aⅴ| 高清av免费在线| 精品一区二区三区视频在线观看免费 | 久久精品人人爽人人爽视色| 欧美大码av| 久久精品aⅴ一区二区三区四区| 国产精品久久电影中文字幕 | 99精国产麻豆久久婷婷| 黄色成人免费大全| 一边摸一边抽搐一进一小说 | 亚洲人成电影观看| 欧美激情久久久久久爽电影 | 精品久久久久久,| 女人精品久久久久毛片| 日韩成人在线观看一区二区三区| 精品熟女少妇八av免费久了| 婷婷精品国产亚洲av在线 | 免费黄频网站在线观看国产| 国产精品亚洲av一区麻豆| 亚洲视频免费观看视频| 国产aⅴ精品一区二区三区波| 90打野战视频偷拍视频| 一级毛片精品| 99国产精品一区二区蜜桃av | 中文字幕人妻丝袜一区二区| 亚洲一区二区三区不卡视频| 中文字幕人妻丝袜制服| 亚洲精品在线美女| 午夜精品久久久久久毛片777| 欧美精品高潮呻吟av久久| 国产单亲对白刺激| 国产又色又爽无遮挡免费看| 久久九九热精品免费| 精品国产亚洲在线| 丝袜人妻中文字幕| 99国产极品粉嫩在线观看| 日本wwww免费看| 国产主播在线观看一区二区| 欧美精品av麻豆av| 丝袜美腿诱惑在线| 亚洲欧美一区二区三区久久| 一区二区三区精品91| 国产亚洲精品久久久久久毛片 | 韩国精品一区二区三区| 在线观看免费午夜福利视频| 一区福利在线观看| svipshipincom国产片| 免费高清在线观看日韩| 黑人操中国人逼视频| 一区二区三区精品91| 超碰97精品在线观看| 成人国语在线视频| 一级a爱片免费观看的视频| 50天的宝宝边吃奶边哭怎么回事| 日韩三级视频一区二区三区| 九色亚洲精品在线播放| 青草久久国产| 精品视频人人做人人爽| 在线播放国产精品三级| 日韩欧美免费精品| 色婷婷久久久亚洲欧美| 黄色片一级片一级黄色片| 国产亚洲精品第一综合不卡| 黄色怎么调成土黄色| 国产免费现黄频在线看| 成人精品一区二区免费| 男女之事视频高清在线观看| 三上悠亚av全集在线观看| 日本wwww免费看| 人妻一区二区av| 亚洲情色 制服丝袜| 高清黄色对白视频在线免费看| 777久久人妻少妇嫩草av网站| aaaaa片日本免费| tube8黄色片| 黄色毛片三级朝国网站| 黑丝袜美女国产一区| 欧美日韩黄片免| 99热国产这里只有精品6| 久久精品国产清高在天天线| 黑人欧美特级aaaaaa片| 天天操日日干夜夜撸| 国产精品免费视频内射| 成人国语在线视频| 99re在线观看精品视频| 色综合婷婷激情| 香蕉国产在线看| 乱人伦中国视频| 午夜免费鲁丝| 国产男女超爽视频在线观看| 亚洲精品国产区一区二| 99久久国产精品久久久| 久久精品aⅴ一区二区三区四区| 黑人欧美特级aaaaaa片| 精品福利永久在线观看| 亚洲精品成人av观看孕妇| 亚洲第一欧美日韩一区二区三区| avwww免费| 亚洲一区二区三区不卡视频| av有码第一页| 精品国产亚洲在线| 一区在线观看完整版| 国产av一区二区精品久久| 成人18禁在线播放| 男女之事视频高清在线观看|