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

    基于ADINA的轎車外部流場計算

    2017-10-16 01:40:15袁小慧歸文強
    西安航空學(xué)院學(xué)報 2017年5期
    關(guān)鍵詞:升力轎車流場

    袁小慧,歸文強

    (西安航空學(xué)院 車輛工程學(xué)院,西安 710077)

    基于ADINA的轎車外部流場計算

    袁小慧,歸文強

    (西安航空學(xué)院 車輛工程學(xué)院,西安 710077)

    對汽車外部流場進行計算流體力學(xué)(CFD)分析已成為現(xiàn)今車身設(shè)計的必要環(huán)節(jié)。計算轎車的外部流場,將計算流體力學(xué)軟件ADINA與CATIA相結(jié)合,利用CATIA軟件獲得轎車三維模型,將其導(dǎo)入 ADINA軟件中,采用RNG K-ε模型,對兩種不同車速40m/s和20m/s情形進行計算機仿真和后期處理,獲得轎車外部流場的壓力分布、流速分布云圖、空氣阻力和空氣升力,進而得到該轎車模型的空氣阻力系數(shù)為0.59,20m/s和40m/s速度載荷下的空氣升力系數(shù)分別為0.013和0.133,依此為汽車車身設(shè)計提供依據(jù)。

    空氣動力學(xué);外部流場;轎車

    汽車(特別是轎車)的空氣動力學(xué)特征直接影響其動力性、燃油經(jīng)濟性、操縱穩(wěn)定性、舒適性和安全性等重要特性[1]。因此,研究轎車的外部流場對于轎車的研發(fā)與評價至關(guān)重要并將逐步完善。目前,主要采用數(shù)值模擬與風(fēng)洞試驗相結(jié)合的方法進行研究。相對后者,前者方便改善性能,節(jié)約研究資金,提高研究效率,省時省工[2]。計算流體力學(xué)在汽車中的應(yīng)用始于20世紀80年代的歐美地區(qū),現(xiàn)已發(fā)展到包括倒車鏡、擾流板、復(fù)雜底板、移動地面、發(fā)動機艙、橫風(fēng)過度特性等多樣化計算模擬。國內(nèi)起步較晚,目前出現(xiàn)的由我國自主開發(fā)的應(yīng)用于三維流場汽車的計算軟件,是胡善龍等提出的數(shù)字樣車(DMU)技術(shù)中前期設(shè)計的一個重要組成部分[3]。多物理場耦合方面的幾大軟件,諸如ANSYS、ADINA、ABAQUS、MSC都可以做到結(jié)構(gòu)、流體和熱的耦合分析。ADINA可以直接流固耦合,但其它三個軟件必須與其他軟件聯(lián)合使用,進行迭代分析。

    1 理論研究

    汽車周圍流場的典型流動特征為三維、粘性、分離和非定常。本文將研究空氣視為“連續(xù)介質(zhì)”的理想氣體,只考慮外力作用下的宏觀運動,因流速與聲速相比足夠小,且馬赫數(shù)Ma=U/α≤0.25,可認為汽車周圍流場不可壓縮。假設(shè)汽車周圍流場溫度不變,具有各向同性,且不考慮粘性,而流場中任一質(zhì)點介質(zhì)具有定常流動特點,可采用K-ε數(shù)值分析方法。

    1.1 計算流體力學(xué)基礎(chǔ)

    由汽車雷諾數(shù)Rei=v∞l/v>104及上述分析,汽車外部流場為定常、等溫和不可壓三維流場,因流場中各物理量在隨時間和空間變化時都符合統(tǒng)計規(guī)律,采用時均法對汽車周圍流場的三維湍流進行數(shù)值模擬,其控制方程為

    (1)

    (2)

    湍流動能方程(即K方程)

    (3)

    湍流耗散方程(即ε方程)

    (4)

    (5)

    (6)

    湍流動能k為湍流中單位流體的平均脈動動能,即

    (7)

    湍流耗散ε為湍流中單位質(zhì)量流體的脈動粘性耗散,即

    (8)

    其中,K表示流體的彈性模量,即體積彈性模量;V表示來流速度;ρ表示流體密度[4]。

    考慮到湍流中渦流因素的影響和雷諾數(shù)效應(yīng)時RNG K-ε模型比標準K-ε計算精度更高,需要的內(nèi)存少,較適用于汽車外部復(fù)雜流場的計算,本文采用RNG K-ε模型,將高雷諾方程與壁函數(shù)結(jié)合起來使用。

    1.2 計算流體力學(xué)軟件ADINA簡介

    CFD方法求解汽車周圍三維湍流流場是采用計算數(shù)學(xué)的算法,將其控制方程離散到一系列網(wǎng)格節(jié)點上并求其離散的數(shù)值解。按其采用的數(shù)值解題方式可以歸納為有限差分法(Finite Different Method,F(xiàn)DM)、有限元法(Finite Element Method,F(xiàn)EM)和有限體積法(Finite Volume Method,BEM)三大類[5]。

    有限體積法計算效率高,近年來發(fā)展迅速,目前廣泛應(yīng)用于CFD領(lǐng)域,大多數(shù)商用CFD軟件,包括ADINA軟件,均采用這種方法。ADINA具有強大的非線性求解能力,能夠解決包括結(jié)構(gòu)、流體、熱等的很多耦合場問題。在針對流體流動模擬仿真時,基本假設(shè)如下:流場計算采用 Navier-Stokes 方程或者歐拉方程;認為流體為不可壓縮或完全可壓縮;分析時可進行穩(wěn)態(tài)或者瞬變分析;流體類型可根據(jù)實際情況選擇層流或者湍流;流體流動時可設(shè)置為熱傳遞或者無熱傳遞條件下的流動;此外,還有質(zhì)量傳遞。

    2 模型與計算

    根據(jù)ADINA計算需要,建立簡化轎車模型,并將其轉(zhuǎn)換成于Parasolid格式,導(dǎo)入ADINA中,定義外部流場,劃分單元網(wǎng)格,生成有限元模型;定義材料及邊界條件,進行求解,獲得外部流場壓力分布、流速分布及空氣阻力系數(shù),并與實際相比較,分析所得結(jié)果。

    本文采用RNGK-ε模型進行計算,除去其本身自帶物理量外,還用到的物理量如下:空氣密度ρ=1.225kg/m3,動力粘性系數(shù)μ=2×10-5Ns/m2,速度載荷40m/s和20m/s。

    2.1 模型的建立

    2.1.1 幾何模型

    查閱資料,得到奧迪A5的CAD圖,確定其原始幾何模型參數(shù)4625mm×1854mm×1372mm,并將上述CAD圖導(dǎo)入三維建模軟件CATIA中,得到豪華轎車奧迪A5的CATIA模型。在實際應(yīng)用中,可以先建立CATIA模型,導(dǎo)入分析軟件,根據(jù)計算結(jié)果對其進行優(yōu)化設(shè)計;或者將現(xiàn)車CAD圖導(dǎo)入CATIA中建立三維模型進行計算,分析其優(yōu)化可能性及主要方向。仿真分析的模型為1:1模型,并對實車模型作了如下簡化。

    忽略車身外部突起物如刮雨器、后視鏡等部分,將車輪與車身做成一體;對車身底部進行平整處理。這些改變對流場沒有太大影響,并且降低了計算的強度和對計算機的要求,提高了計算經(jīng)濟性,符合計算要求。圖1所示為CATIA中建立的1:1模型的截圖。

    圖1 轎車CATIA模型

    2.1.2 物理模型

    汽車周圍的流體即空氣,在一個標準大氣壓、海平面高度下,15℃時的空氣密度ρ=1.225kg/m3,動力粘性系數(shù)μ=1.7894×10-5Ns/m2,運動粘度v=μ/ρ=1.4607×10-5m/s。

    由于本文數(shù)值模擬風(fēng)速為40m/s和20m/s,可以認為汽車繞流不可壓縮。據(jù)本文第1節(jié),將汽車繞流簡化為定常、等溫、不可壓縮的三維流場。考慮到氣流的分離影響,按湍流處理,采用RNGK-ε模型,將外部流場視為15m×5m×5m的長方體除去轎車模型所占空間。

    2.2 外部流場計算

    定義模型控制參數(shù),將得到的三維模型導(dǎo)入分析軟件ADINA中,如圖2所示。

    圖2 奧迪A5模型導(dǎo)入ADINA圖

    設(shè)汽車行駛方向為Y軸正方向,垂直于地面向上的方向為Z軸正方向。建立流場模型,定義約束與載荷,輸入邊界條件。其中,車輛外部上下表面為wall類型,其他邊界條件如表1所示。

    表1 邊界條件的設(shè)定

    定義材料(選擇RNG K-Epsilon,密度1.225,粘度為2×10-5)與單元組類型(網(wǎng)格類型為3-D Fluid,單元組為1,單元節(jié)點數(shù)為4),并根據(jù)流場復(fù)雜程度進行網(wǎng)格劃分,網(wǎng)格數(shù)為70683,如圖3所示。

    圖3 單元網(wǎng)格

    本文主要對多數(shù)轎車正常行駛速度范圍內(nèi)的兩個速度(二倍關(guān)系速度)下的壓力分布和流速分布進行對比分析。進口邊界條件見表1,求解,進行計算流體力學(xué)CFD分析以及后處理,獲得其外部流場流速分布圖和壓力分布圖,見圖4-圖7。

    圖4 40m/s的流場對應(yīng)的壓力分布

    圖5 20m/s的流場對應(yīng)的壓力分布

    圖6 40m/s的流場對應(yīng)的流速分布

    圖7 20m/s的流場對應(yīng)的流速分布

    由圖4和圖5可知,負壓區(qū)(藍色區(qū)域)主要出現(xiàn)在下壁面以下、車頂以及車尾部分;正壓區(qū)主要出現(xiàn)在車頭以及發(fā)動機艙蓋與前擋風(fēng)玻璃交接處。其中,保險杠附近出現(xiàn)正壓最大值,下壁面下邊的最前端出現(xiàn)負壓最大值,這是由于作為迎風(fēng)面的車頭直接受到氣流的沖擊作用而出現(xiàn)正壓的最大值,而汽車下壁面則由于離地面距離小出現(xiàn)負壓的最大值。

    由圖6和圖7可知,汽車前部氣流分布比較均勻,而汽車尾部的氣流較為紊亂,形成明顯的漩渦,流速達最大值,在遠離壁面區(qū)域流速接近速度載荷,仿真結(jié)果見表2。

    表2 仿真結(jié)果

    2.3 結(jié)果分析

    將所得結(jié)果與一般轎車外部流場壓力分布圖進行對比說明,并對比40m/s與20m/s對應(yīng)的壓力分布圖與流速分布圖,分析其中的差異,針對建立的模型提出改進方案。其中,一般轎車外部流場壓力分布圖如圖8所示[6]。

    圖8 轎車外部流場壓力分布圖

    凹向車內(nèi)的區(qū)域即為正壓力分布區(qū)域,車外的區(qū)域即為負壓力分布區(qū)域。由此可見,上述壓力分布計算結(jié)果(圖4與圖5)正確。對比兩種不同速度載荷下的流場壓力分布與流速分布可知:兩者壓力分布與流速分布圖相似,只是對應(yīng)數(shù)值有差異;速度載荷為20m/s的流場壓力最大值203.2N,最小值-211.5N,與速度載荷為40m/s的流場壓力的最大值813.4N,最小值-844.4N相比,約是1/4倍數(shù)關(guān)系,而速度載荷僅是1/2倍數(shù)關(guān)系,可見壓力變化幅度是速度變化幅度的平方關(guān)系;速度載荷為20m/s的流場流速最大值39.69m/s,與速度載荷為40m/s的流場流速最大值78.61m/s相比,約是1/2倍數(shù)關(guān)系,與速度載荷變化相同。對比兩者空氣阻力即1388.1N與348.9N可知,速度載荷為二倍關(guān)系,對應(yīng)空氣阻力即為四倍關(guān)系。

    根據(jù)空氣阻力公式Fy=0.5·ρ·v2·S·CD和空氣升力公式Fz=0.5·ρ·v2·S·CL,可得空氣阻力系數(shù)CD=2Fy/(ρ·v2·S)和空氣升力系數(shù)CL=2Fy/(ρ·v2·S),其中ρ為空氣密度,v為相對速度(車速),迎風(fēng)面積S為2.4m2。根據(jù)仿真結(jié)果,可得該汽車模型在不同速度載荷下的空氣動力學(xué)性能數(shù)據(jù),如表3所示。

    表3 空氣動力學(xué)性能結(jié)果

    本文模型為簡化模型,兩種速度載荷下空氣阻力系數(shù)均為0.59。對比20m/s和40m/s速度載荷下的空氣升力系數(shù)值,可得20m/s速度載荷下的升力系數(shù)僅為40m/s速度載荷下的升力系數(shù)的0.1,因此,車輛在高速行駛時,空氣升力對車輛的操縱穩(wěn)定性影響很大。將空氣阻力系數(shù)和空氣升力系數(shù)作為汽車空氣動力學(xué)性能的評價指標,優(yōu)化汽車模型,可以通過兼顧美學(xué)要求,改變汽車前風(fēng)窗玻璃角度、前圍板結(jié)構(gòu)、前后車燈結(jié)構(gòu)、后車窗玻璃、尾部造型等,達到優(yōu)化車型的目的。

    3 結(jié)論

    本文基于奧迪A5的CAD圖建立CATIA模型,獲得了轎車外部流場在速度載荷分別為40m/s與20m/s下的壓力分布和流速分布云圖,進而得到該轎車模型的空氣阻力系數(shù)為0.59,20m/s速度載荷下的空氣升力系數(shù)為0.013,40m/s速度載荷下的空氣升力系數(shù)為0.133,為獲得理想空氣動力學(xué)特性的車型提供了參考。另外,本文在簡化模型基礎(chǔ)上進行了理論探討,僅考慮了40m/s與20m/s的穩(wěn)態(tài)情況,未對來流大小、方向和湍流度等隨時間的變化情況進行瞬態(tài)模擬,在實際應(yīng)用中還需要進行風(fēng)洞試驗,考慮諸多因素以接近路面的真實情況。

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

    [2] 閆江.汽車外流場數(shù)值模擬研究[C]//中國汽車工程學(xué)會.2013中國汽車工程學(xué)會年會論文集,北京,2013:4.

    [3] 胡善龍,錢榕,牛勝福,等.數(shù)字樣車技術(shù)在現(xiàn)代轎車開發(fā)中的應(yīng)用及發(fā)展前景[J].上海汽車,2009,35(2):32-35.

    [4] 李毓洲.基于ANSYS的汽車空氣動力學(xué)特性分析[J].機械設(shè)計與制造,2010,47(2):113-115.

    [5] 曲震,薛澄岐,韓飛聽.基于CFD方法的轎車車身外圍流場分析[J].電子機械工程,2009, 25(1):33-41.

    [6] 杜廣生.汽車空氣動力學(xué)[M].北京:中國標準出版社,1999:132.

    [責(zé)任編輯、校對:東艷]

    Abstract:Nowadays,it has become an essential part of the body design to analyze the car external flow field using CFD method.This article is mainly to calculate the car external flow field.The computational fluid dynamics software ADINA and CATIA are combined.The car 3-D model,obtained through CATIA software,is imported into ADINA software.Through the RNG K-ε model,and under the velocities of 20m/s and 40m/s, both pressure cloud picture and velocity images are obtained,with pressure gotten.The air resistance coefficient of the car model is 0.59.The air lift coefficients at 20m/s and 40m/s are 0.013 and 0.133 respectively.The results can provide the basis for car body design.

    Keywords:aerodynamics;external flow field;car

    AnalysisofCarExternalFlowFieldBasedonADINA

    YUANXiao-hui,GUIWen-qiang

    (School of Vehicle Engineering,Xi′an Aeronautical University,Xi′an 710077,China)

    U461.1

    A

    1008-9233(2017)05-0077-04

    2017-03-14

    袁小慧(1989-),女,河南方城人,助教,主要從事汽車主動安全技術(shù)研究。

    猜你喜歡
    升力轎車流場
    高速列車車頂–升力翼組合體氣動特性
    大型空冷汽輪發(fā)電機轉(zhuǎn)子三維流場計算
    一汽奔騰轎車有限公司
    汽車文摘(2021年2期)2021-02-01 07:02:36
    無人機升力測試裝置設(shè)計及誤差因素分析
    等著的轎車
    基于自適應(yīng)偽譜法的升力式飛行器火星進入段快速軌跡優(yōu)化
    坐著轎車學(xué)聲調(diào)
    轉(zhuǎn)杯紡排雜區(qū)流場與排雜性能
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場統(tǒng)計分析
    升力式再入飛行器體襟翼姿態(tài)控制方法
    99re在线观看精品视频| 免费在线观看成人毛片| 国产成人av激情在线播放| 国产精品99久久久久久久久| www.自偷自拍.com| 丰满的人妻完整版| 18禁裸乳无遮挡免费网站照片| 亚洲aⅴ乱码一区二区在线播放| 成人性生交大片免费视频hd| 一区二区三区国产精品乱码| 精品国产超薄肉色丝袜足j| 婷婷六月久久综合丁香| 国产精品爽爽va在线观看网站| 麻豆av在线久日| 老司机福利观看| 亚洲成av人片在线播放无| 99久久久亚洲精品蜜臀av| 欧美黄色片欧美黄色片| 麻豆国产av国片精品| 久久午夜综合久久蜜桃| 九九在线视频观看精品| 亚洲精品色激情综合| 国产成人精品久久二区二区91| 夜夜夜夜夜久久久久| 一级毛片高清免费大全| 女人高潮潮喷娇喘18禁视频| 香蕉丝袜av| 一进一出抽搐动态| 国产伦在线观看视频一区| 神马国产精品三级电影在线观看| 亚洲国产欧美人成| 日本撒尿小便嘘嘘汇集6| 一区二区三区国产精品乱码| 亚洲自拍偷在线| 极品教师在线免费播放| 久久草成人影院| 亚洲专区字幕在线| 18禁黄网站禁片免费观看直播| 国产乱人视频| 精品久久久久久久毛片微露脸| 在线观看免费午夜福利视频| 在线观看舔阴道视频| 波多野结衣高清无吗| 国产一区二区激情短视频| 午夜精品久久久久久毛片777| 国产精品一区二区免费欧美| 熟女电影av网| 国产av不卡久久| 真人一进一出gif抽搐免费| 国产又黄又爽又无遮挡在线| 亚洲中文字幕一区二区三区有码在线看 | 欧美一区二区精品小视频在线| 色吧在线观看| 免费看十八禁软件| 高潮久久久久久久久久久不卡| 成人特级黄色片久久久久久久| 国产精品一及| 亚洲欧美日韩高清专用| 中文字幕高清在线视频| 久久精品人妻少妇| 日韩欧美免费精品| 国产精品 欧美亚洲| 欧美日韩瑟瑟在线播放| 免费av毛片视频| 一边摸一边抽搐一进一小说| 日韩大尺度精品在线看网址| 男人舔奶头视频| 欧美性猛交黑人性爽| 99热6这里只有精品| 日日摸夜夜添夜夜添小说| 精品一区二区三区视频在线观看免费| 在线免费观看的www视频| 美女cb高潮喷水在线观看 | 国产成人福利小说| 欧美激情在线99| 又爽又黄无遮挡网站| 亚洲午夜理论影院| 亚洲成人精品中文字幕电影| 岛国视频午夜一区免费看| 一个人观看的视频www高清免费观看 | 亚洲av电影不卡..在线观看| 黑人巨大精品欧美一区二区mp4| 久久天堂一区二区三区四区| 毛片女人毛片| 国产av麻豆久久久久久久| 999久久久国产精品视频| 午夜福利免费观看在线| 99久国产av精品| 国产精品一区二区三区四区免费观看 | 国产真人三级小视频在线观看| 国产伦人伦偷精品视频| 国内久久婷婷六月综合欲色啪| 又紧又爽又黄一区二区| 88av欧美| 国产视频内射| 手机成人av网站| 久久精品夜夜夜夜夜久久蜜豆| 国产精品 国内视频| 高清在线国产一区| 两个人的视频大全免费| 国产精品女同一区二区软件 | 黄片小视频在线播放| 青草久久国产| 日韩大尺度精品在线看网址| 日韩精品青青久久久久久| 中文字幕熟女人妻在线| 婷婷丁香在线五月| 亚洲 欧美 日韩 在线 免费| a级毛片a级免费在线| 两人在一起打扑克的视频| 欧美乱码精品一区二区三区| 亚洲美女视频黄频| 久久久国产成人免费| 中文字幕av在线有码专区| av中文乱码字幕在线| 男女之事视频高清在线观看| 国产aⅴ精品一区二区三区波| 全区人妻精品视频| 90打野战视频偷拍视频| 亚洲av电影不卡..在线观看| 黄色 视频免费看| 国产一区二区在线av高清观看| 黄色片一级片一级黄色片| 亚洲国产日韩欧美精品在线观看 | 岛国在线免费视频观看| 看免费av毛片| 日韩成人在线观看一区二区三区| 国产成人精品无人区| 国产高清videossex| 人妻夜夜爽99麻豆av| 成人特级av手机在线观看| 国产97色在线日韩免费| 麻豆成人av在线观看| 亚洲激情在线av| 一级作爱视频免费观看| 免费在线观看视频国产中文字幕亚洲| 99久久精品热视频| 一进一出抽搐动态| 老司机在亚洲福利影院| 麻豆国产97在线/欧美| 在线看三级毛片| 国产精品一区二区三区四区免费观看 | 丰满人妻一区二区三区视频av | 99久久成人亚洲精品观看| svipshipincom国产片| 搡老熟女国产l中国老女人| 国产毛片a区久久久久| 亚洲 国产 在线| 国产一区在线观看成人免费| 噜噜噜噜噜久久久久久91| 俺也久久电影网| 欧美中文日本在线观看视频| 中文字幕高清在线视频| 国产亚洲av高清不卡| 久久国产乱子伦精品免费另类| 黄片小视频在线播放| 免费人成视频x8x8入口观看| av女优亚洲男人天堂 | 日本与韩国留学比较| 国产av不卡久久| 国产精品 国内视频| 一个人看的www免费观看视频| 亚洲真实伦在线观看| 日本精品一区二区三区蜜桃| 一边摸一边抽搐一进一小说| 国产激情偷乱视频一区二区| 久久亚洲精品不卡| 亚洲人成伊人成综合网2020| 亚洲熟女毛片儿| 欧美日韩中文字幕国产精品一区二区三区| 日本三级黄在线观看| 欧美三级亚洲精品| 老汉色av国产亚洲站长工具| 欧美极品一区二区三区四区| 黄色日韩在线| 午夜免费观看网址| 首页视频小说图片口味搜索| 99视频精品全部免费 在线 | 成人特级黄色片久久久久久久| 午夜福利免费观看在线| 国产aⅴ精品一区二区三区波| 国产三级中文精品| 日韩高清综合在线| 搞女人的毛片| 老汉色av国产亚洲站长工具| 国产精品 欧美亚洲| 日本免费一区二区三区高清不卡| 精品久久久久久久末码| 亚洲专区字幕在线| 九色成人免费人妻av| 欧美一级毛片孕妇| 我的老师免费观看完整版| 亚洲九九香蕉| 免费在线观看日本一区| 少妇人妻一区二区三区视频| 成人性生交大片免费视频hd| 亚洲人成网站在线播放欧美日韩| 狠狠狠狠99中文字幕| 亚洲人成网站高清观看| 国内揄拍国产精品人妻在线| 在线看三级毛片| 在线观看免费视频日本深夜| а√天堂www在线а√下载| 99热这里只有是精品50| 一区福利在线观看| 亚洲成a人片在线一区二区| 又紧又爽又黄一区二区| 最近最新中文字幕大全免费视频| 亚洲中文日韩欧美视频| 亚洲,欧美精品.| 欧美3d第一页| 亚洲专区中文字幕在线| 亚洲成av人片免费观看| 一个人看的www免费观看视频| 久久久久久久午夜电影| 日本一本二区三区精品| 国产不卡一卡二| 亚洲人成电影免费在线| 熟女电影av网| 男女视频在线观看网站免费| 好看av亚洲va欧美ⅴa在| 亚洲,欧美精品.| 国产单亲对白刺激| 日韩精品青青久久久久久| 久久久久国内视频| 亚洲午夜精品一区,二区,三区| 女生性感内裤真人,穿戴方法视频| 特大巨黑吊av在线直播| 欧美又色又爽又黄视频| 国产淫片久久久久久久久 | 亚洲国产中文字幕在线视频| 国产av在哪里看| 国产精华一区二区三区| 一个人看的www免费观看视频| 久久久久久大精品| 亚洲欧美日韩高清在线视频| 久久精品国产清高在天天线| 国产v大片淫在线免费观看| 国产高清视频在线播放一区| 免费观看的影片在线观看| 超碰成人久久| 琪琪午夜伦伦电影理论片6080| 在线播放国产精品三级| 日韩 欧美 亚洲 中文字幕| 啪啪无遮挡十八禁网站| 搡老熟女国产l中国老女人| 日韩欧美三级三区| 国产精品永久免费网站| 美女午夜性视频免费| 亚洲成人久久爱视频| 又紧又爽又黄一区二区| 午夜亚洲福利在线播放| 一级毛片女人18水好多| 久久中文看片网| 成熟少妇高潮喷水视频| 亚洲乱码一区二区免费版| 黄片大片在线免费观看| 在线a可以看的网站| 俄罗斯特黄特色一大片| 一区二区三区国产精品乱码| 哪里可以看免费的av片| 久久久精品大字幕| 18禁国产床啪视频网站| 99久久成人亚洲精品观看| 天天躁狠狠躁夜夜躁狠狠躁| aaaaa片日本免费| 精品国产亚洲在线| 俺也久久电影网| 757午夜福利合集在线观看| 色综合站精品国产| 啦啦啦韩国在线观看视频| 草草在线视频免费看| 99久久国产精品久久久| 久久久久久人人人人人| 嫁个100分男人电影在线观看| 一个人看视频在线观看www免费 | 免费看日本二区| 久久久久久人人人人人| 99在线视频只有这里精品首页| 亚洲午夜精品一区,二区,三区| 国产精品自产拍在线观看55亚洲| 中文字幕av在线有码专区| 五月玫瑰六月丁香| 两个人看的免费小视频| 国产免费av片在线观看野外av| 国产69精品久久久久777片 | 欧美国产日韩亚洲一区| 国产主播在线观看一区二区| 久久久久国产一级毛片高清牌| 特级一级黄色大片| 两个人的视频大全免费| 精品一区二区三区视频在线 | 午夜精品久久久久久毛片777| 国产毛片a区久久久久| 国产蜜桃级精品一区二区三区| 欧美激情在线99| 国产又黄又爽又无遮挡在线| 日本黄色视频三级网站网址| 午夜免费观看网址| 夜夜爽天天搞| 色尼玛亚洲综合影院| 亚洲av五月六月丁香网| 久久久成人免费电影| 在线观看午夜福利视频| 亚洲在线自拍视频| 精品电影一区二区在线| 成在线人永久免费视频| 久久草成人影院| 99视频精品全部免费 在线 | 91在线观看av| 久久亚洲精品不卡| 好男人在线观看高清免费视频| 琪琪午夜伦伦电影理论片6080| 欧美一级毛片孕妇| 成年女人看的毛片在线观看| 噜噜噜噜噜久久久久久91| 黄色丝袜av网址大全| 一进一出抽搐gif免费好疼| 欧美黄色片欧美黄色片| 九色国产91popny在线| 久久99热这里只有精品18| 国产亚洲欧美98| 色老头精品视频在线观看| 精品国产三级普通话版| 亚洲精华国产精华精| 亚洲av成人精品一区久久| 日韩欧美一区二区三区在线观看| 高清在线国产一区| 国产欧美日韩一区二区三| 美女扒开内裤让男人捅视频| 男女下面进入的视频免费午夜| 久久久久性生活片| 午夜免费成人在线视频| 国产亚洲欧美98| 久久久久亚洲av毛片大全| 成人性生交大片免费视频hd| 国产乱人视频| 国内精品久久久久久久电影| av天堂在线播放| 女警被强在线播放| 国产精品影院久久| 日本黄色片子视频| 成人三级黄色视频| 亚洲欧美精品综合久久99| 国产成人系列免费观看| 亚洲av电影不卡..在线观看| 久久这里只有精品中国| 精品久久久久久成人av| 久久婷婷人人爽人人干人人爱| 国产成人啪精品午夜网站| av国产免费在线观看| 99久久国产精品久久久| 亚洲真实伦在线观看| aaaaa片日本免费| 神马国产精品三级电影在线观看| 狂野欧美白嫩少妇大欣赏| 午夜福利高清视频| 亚洲人成伊人成综合网2020| 亚洲美女黄片视频| 特大巨黑吊av在线直播| 久久精品aⅴ一区二区三区四区| 成人亚洲精品av一区二区| 国产精品久久久av美女十八| 亚洲avbb在线观看| 日本 av在线| 国产精品久久电影中文字幕| 黄色女人牲交| 男女那种视频在线观看| АⅤ资源中文在线天堂| 麻豆一二三区av精品| 亚洲av日韩精品久久久久久密| 波多野结衣巨乳人妻| 久久久国产精品麻豆| 少妇的逼水好多| 亚洲中文av在线| 国产欧美日韩精品一区二区| 嫩草影院精品99| 亚洲成av人片在线播放无| 亚洲性夜色夜夜综合| 成人精品一区二区免费| 免费大片18禁| 久久久久国产精品人妻aⅴ院| 成人午夜高清在线视频| 亚洲精品一卡2卡三卡4卡5卡| 久久精品91无色码中文字幕| 国产午夜福利久久久久久| 亚洲av第一区精品v没综合| 国产伦精品一区二区三区四那| 琪琪午夜伦伦电影理论片6080| 小蜜桃在线观看免费完整版高清| 国产精品av久久久久免费| 欧美乱色亚洲激情| 欧美激情在线99| h日本视频在线播放| 精品久久久久久久末码| 欧美日韩黄片免| 别揉我奶头~嗯~啊~动态视频| 每晚都被弄得嗷嗷叫到高潮| 午夜激情欧美在线| 少妇人妻一区二区三区视频| 午夜成年电影在线免费观看| 宅男免费午夜| 成人欧美大片| 99国产综合亚洲精品| 无限看片的www在线观看| 国产高清有码在线观看视频| 精品不卡国产一区二区三区| 精品一区二区三区视频在线 | 视频区欧美日本亚洲| 在线国产一区二区在线| 黄频高清免费视频| 亚洲精品一卡2卡三卡4卡5卡| 十八禁人妻一区二区| 国内精品美女久久久久久| av女优亚洲男人天堂 | 亚洲男人的天堂狠狠| 99久久精品热视频| 怎么达到女性高潮| 舔av片在线| 老司机深夜福利视频在线观看| 久久久久久久精品吃奶| 夜夜看夜夜爽夜夜摸| 久9热在线精品视频| 欧美+亚洲+日韩+国产| 久久久久性生活片| 少妇丰满av| 黄色丝袜av网址大全| 深夜精品福利| 色综合婷婷激情| 91av网站免费观看| 亚洲精品一卡2卡三卡4卡5卡| 嫁个100分男人电影在线观看| 悠悠久久av| 很黄的视频免费| 美女免费视频网站| 国产综合懂色| 视频区欧美日本亚洲| 国产伦一二天堂av在线观看| 亚洲国产欧美一区二区综合| av在线蜜桃| 很黄的视频免费| 搡老岳熟女国产| 黄色丝袜av网址大全| 亚洲精品美女久久久久99蜜臀| 久久久久精品国产欧美久久久| 国产av在哪里看| 成人亚洲精品av一区二区| 在线视频色国产色| 天堂影院成人在线观看| 国产午夜精品久久久久久| 人人妻人人看人人澡| 一个人免费在线观看电影 | 熟女人妻精品中文字幕| 少妇人妻一区二区三区视频| 中文在线观看免费www的网站| 亚洲熟妇熟女久久| 国产伦一二天堂av在线观看| 麻豆av在线久日| 人妻夜夜爽99麻豆av| 老司机午夜十八禁免费视频| 精品国产三级普通话版| 国产探花在线观看一区二区| 少妇裸体淫交视频免费看高清| 国产精品一区二区精品视频观看| 一级作爱视频免费观看| 亚洲国产精品999在线| 日韩欧美在线二视频| 亚洲国产欧洲综合997久久,| 搡老熟女国产l中国老女人| 啪啪无遮挡十八禁网站| 每晚都被弄得嗷嗷叫到高潮| 少妇的丰满在线观看| 国产精品久久久人人做人人爽| 一级毛片女人18水好多| 午夜亚洲福利在线播放| 国产精品久久久久久久电影 | 久久精品国产综合久久久| 日日干狠狠操夜夜爽| 亚洲成人免费电影在线观看| 很黄的视频免费| 国内揄拍国产精品人妻在线| 精品国内亚洲2022精品成人| 日本免费a在线| av福利片在线观看| bbb黄色大片| 美女午夜性视频免费| 真实男女啪啪啪动态图| 舔av片在线| 婷婷精品国产亚洲av| 1024手机看黄色片| 97人妻精品一区二区三区麻豆| АⅤ资源中文在线天堂| 欧美绝顶高潮抽搐喷水| 这个男人来自地球电影免费观看| 日韩精品青青久久久久久| 一卡2卡三卡四卡精品乱码亚洲| 美女高潮喷水抽搐中文字幕| 亚洲国产日韩欧美精品在线观看 | 黄色丝袜av网址大全| 国产午夜精品久久久久久| 免费看光身美女| 国产亚洲精品久久久com| av天堂在线播放| 久久伊人香网站| 成年女人永久免费观看视频| 手机成人av网站| 国产淫片久久久久久久久 | 国产精品一区二区三区四区免费观看 | 久99久视频精品免费| 巨乳人妻的诱惑在线观看| 国产 一区 欧美 日韩| 久久久久久国产a免费观看| 岛国在线观看网站| 亚洲五月婷婷丁香| 久久国产乱子伦精品免费另类| 国产 一区 欧美 日韩| 国产精品自产拍在线观看55亚洲| 女同久久另类99精品国产91| 很黄的视频免费| 久久久成人免费电影| 国产成人精品久久二区二区免费| 嫩草影院入口| 亚洲av成人不卡在线观看播放网| 国产成人啪精品午夜网站| 亚洲精品国产精品久久久不卡| 欧美激情久久久久久爽电影| 国产真实乱freesex| 99re在线观看精品视频| 窝窝影院91人妻| 日韩av在线大香蕉| 免费在线观看影片大全网站| 成熟少妇高潮喷水视频| 亚洲精品国产精品久久久不卡| 国产精品国产高清国产av| 好男人在线观看高清免费视频| 亚洲成人中文字幕在线播放| 国产探花在线观看一区二区| 91久久精品国产一区二区成人 | 啦啦啦免费观看视频1| 高潮久久久久久久久久久不卡| 一夜夜www| 人妻夜夜爽99麻豆av| 精品午夜福利视频在线观看一区| 五月玫瑰六月丁香| 久久久久国产精品人妻aⅴ院| av女优亚洲男人天堂 | 欧美最黄视频在线播放免费| 亚洲av熟女| 亚洲专区中文字幕在线| 国产精品永久免费网站| 嫩草影院精品99| 全区人妻精品视频| 首页视频小说图片口味搜索| 亚洲熟妇中文字幕五十中出| 又黄又爽又免费观看的视频| 女警被强在线播放| 国产精品久久久久久精品电影| 久久久久久久久中文| www.www免费av| 无遮挡黄片免费观看| 国内精品久久久久精免费| 一夜夜www| 一区二区三区高清视频在线| 精品久久蜜臀av无| 女生性感内裤真人,穿戴方法视频| 国产精品一区二区精品视频观看| aaaaa片日本免费| 欧美黄色片欧美黄色片| 少妇熟女aⅴ在线视频| 最近在线观看免费完整版| 男女下面进入的视频免费午夜| 亚洲中文日韩欧美视频| 日韩高清综合在线| 91av网一区二区| 国产激情欧美一区二区| 色精品久久人妻99蜜桃| 亚洲精品在线观看二区| 两个人看的免费小视频| 国产成人精品无人区| 亚洲av成人不卡在线观看播放网| 久久人妻av系列| а√天堂www在线а√下载| 国内少妇人妻偷人精品xxx网站 | 久久久久久国产a免费观看| 久久国产精品影院| 又爽又黄无遮挡网站| 国产一区二区在线av高清观看| h日本视频在线播放| e午夜精品久久久久久久| 五月玫瑰六月丁香| 亚洲色图av天堂| 成年免费大片在线观看| 亚洲精品美女久久av网站| 黄色成人免费大全| 欧美黄色淫秽网站| 国产成人福利小说| 在线观看免费视频日本深夜| 99久国产av精品| 午夜福利免费观看在线| 精品一区二区三区视频在线观看免费| 亚洲av电影不卡..在线观看| or卡值多少钱| av在线蜜桃| 99久久成人亚洲精品观看| 国产黄色小视频在线观看| 好男人电影高清在线观看| 久久热在线av| 欧美在线一区亚洲| 男女午夜视频在线观看| 国内精品久久久久精免费| 色精品久久人妻99蜜桃| 在线视频色国产色| 波多野结衣巨乳人妻| 给我免费播放毛片高清在线观看| 97人妻精品一区二区三区麻豆|