• <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)控制方法
    亚洲一区二区三区色噜噜| 最近视频中文字幕2019在线8| 啦啦啦观看免费观看视频高清| 久久天躁狠狠躁夜夜2o2o| 国内少妇人妻偷人精品xxx网站| 亚洲avbb在线观看| 欧美xxxx黑人xx丫x性爽| 人人妻,人人澡人人爽秒播| 免费一级毛片在线播放高清视频| 岛国在线免费视频观看| 久久精品久久久久久噜噜老黄 | 日本撒尿小便嘘嘘汇集6| 精品免费久久久久久久清纯| 久久中文看片网| 亚洲成人精品中文字幕电影| av天堂中文字幕网| 小蜜桃在线观看免费完整版高清| 亚洲天堂国产精品一区在线| 国产69精品久久久久777片| 久久久精品欧美日韩精品| 村上凉子中文字幕在线| 国产精品1区2区在线观看.| 十八禁国产超污无遮挡网站| 亚洲人与动物交配视频| 校园人妻丝袜中文字幕| 精品一区二区免费观看| 欧美三级亚洲精品| 成人高潮视频无遮挡免费网站| 在线播放国产精品三级| 亚洲性夜色夜夜综合| 又紧又爽又黄一区二区| 少妇人妻精品综合一区二区 | 亚洲电影在线观看av| 18禁黄网站禁片午夜丰满| 欧美日韩乱码在线| 午夜影院日韩av| 亚洲性久久影院| 亚洲一区高清亚洲精品| 国内毛片毛片毛片毛片毛片| 国产精品一区二区性色av| 日本一本二区三区精品| 中亚洲国语对白在线视频| 国内精品一区二区在线观看| 免费看av在线观看网站| 欧美激情在线99| 午夜福利欧美成人| 国产精品综合久久久久久久免费| 亚洲专区国产一区二区| 欧美极品一区二区三区四区| 99riav亚洲国产免费| 久久中文看片网| 日本黄大片高清| 波多野结衣巨乳人妻| 男人和女人高潮做爰伦理| 久久久午夜欧美精品| 尾随美女入室| 欧美极品一区二区三区四区| 毛片女人毛片| 超碰av人人做人人爽久久| ponron亚洲| 91在线观看av| 老司机午夜福利在线观看视频| 91久久精品国产一区二区三区| 国产精品98久久久久久宅男小说| 91在线观看av| 99riav亚洲国产免费| 亚洲国产欧洲综合997久久,| 久久精品91蜜桃| 日韩欧美在线二视频| 亚洲精品色激情综合| 看黄色毛片网站| 国产国拍精品亚洲av在线观看| 最近最新中文字幕大全电影3| av在线观看视频网站免费| 国产人妻一区二区三区在| 国产精品1区2区在线观看.| 中亚洲国语对白在线视频| 日韩精品中文字幕看吧| 午夜精品久久久久久毛片777| 国产亚洲精品久久久久久毛片| 麻豆成人午夜福利视频| 亚洲美女黄片视频| 人妻久久中文字幕网| 精品乱码久久久久久99久播| 婷婷精品国产亚洲av| 99久久成人亚洲精品观看| 久久天躁狠狠躁夜夜2o2o| 丰满人妻一区二区三区视频av| 成人特级黄色片久久久久久久| 黄色丝袜av网址大全| 成人性生交大片免费视频hd| a在线观看视频网站| 国产精品日韩av在线免费观看| 亚洲熟妇中文字幕五十中出| 久久精品夜夜夜夜夜久久蜜豆| 国产真实伦视频高清在线观看 | 日本 av在线| 一个人看的www免费观看视频| 中国美白少妇内射xxxbb| 日韩大尺度精品在线看网址| 中文资源天堂在线| 一个人看视频在线观看www免费| 欧美一级a爱片免费观看看| 中文字幕久久专区| 国产欧美日韩精品一区二区| 欧美高清性xxxxhd video| 国产精品永久免费网站| 极品教师在线视频| 久久草成人影院| 直男gayav资源| 男女那种视频在线观看| 精品日产1卡2卡| 久9热在线精品视频| 亚洲最大成人av| 91麻豆精品激情在线观看国产| 成人毛片a级毛片在线播放| 成年版毛片免费区| 99精品在免费线老司机午夜| 成年免费大片在线观看| 国内精品久久久久久久电影| 一区二区三区免费毛片| 国产精品爽爽va在线观看网站| 一级a爱片免费观看的视频| 国产精品98久久久久久宅男小说| 婷婷丁香在线五月| 免费搜索国产男女视频| 亚洲一级一片aⅴ在线观看| av.在线天堂| 熟妇人妻久久中文字幕3abv| 人妻久久中文字幕网| 午夜爱爱视频在线播放| 亚洲精品久久国产高清桃花| 国产成人福利小说| 国产精品无大码| 女同久久另类99精品国产91| 国产色爽女视频免费观看| 麻豆国产av国片精品| 搡老妇女老女人老熟妇| 国产av不卡久久| 亚洲成人中文字幕在线播放| 人妻夜夜爽99麻豆av| 国产精品人妻久久久久久| 亚洲天堂国产精品一区在线| 啦啦啦啦在线视频资源| 男人的好看免费观看在线视频| 男女边吃奶边做爰视频| 精品人妻1区二区| 亚洲av不卡在线观看| 一个人免费在线观看电影| 国产一区二区在线av高清观看| 国产主播在线观看一区二区| 婷婷精品国产亚洲av在线| 色综合色国产| 欧美最黄视频在线播放免费| 直男gayav资源| 在线免费观看不下载黄p国产 | 中文字幕精品亚洲无线码一区| 99九九线精品视频在线观看视频| 热99在线观看视频| 1000部很黄的大片| 亚洲va日本ⅴa欧美va伊人久久| 赤兔流量卡办理| 联通29元200g的流量卡| 国产高清三级在线| 99热这里只有精品一区| 免费观看的影片在线观看| 色哟哟·www| 琪琪午夜伦伦电影理论片6080| 欧美黑人欧美精品刺激| 久久精品影院6| 国产一区二区激情短视频| 最近视频中文字幕2019在线8| 午夜免费激情av| 嫩草影院入口| 一区二区三区四区激情视频 | 天美传媒精品一区二区| 波多野结衣高清作品| 婷婷六月久久综合丁香| 欧美日韩精品成人综合77777| 日韩中文字幕欧美一区二区| 色综合婷婷激情| 亚洲avbb在线观看| 国内精品宾馆在线| 色哟哟·www| 久久婷婷人人爽人人干人人爱| 在线天堂最新版资源| 国产精品99久久久久久久久| 免费观看精品视频网站| 亚洲国产精品久久男人天堂| 欧美成人一区二区免费高清观看| 欧美日本亚洲视频在线播放| 欧美最新免费一区二区三区| 真实男女啪啪啪动态图| 美女cb高潮喷水在线观看| 色综合站精品国产| av在线观看视频网站免费| 国产老妇女一区| 中文字幕av在线有码专区| 久久草成人影院| 国产亚洲欧美98| 日本欧美国产在线视频| 深爱激情五月婷婷| 日本黄大片高清| 中文字幕av在线有码专区| 亚洲av电影不卡..在线观看| 国产视频内射| 美女cb高潮喷水在线观看| 国产亚洲av嫩草精品影院| 精品久久久久久久末码| 午夜爱爱视频在线播放| 久久久国产成人免费| 非洲黑人性xxxx精品又粗又长| 久9热在线精品视频| 久久天躁狠狠躁夜夜2o2o| 国产精品自产拍在线观看55亚洲| 亚洲精品粉嫩美女一区| 久久精品国产清高在天天线| 少妇人妻一区二区三区视频| 亚洲精品亚洲一区二区| 国产精品综合久久久久久久免费| 亚洲无线观看免费| 国产精品野战在线观看| 日韩,欧美,国产一区二区三区 | 久久热精品热| 国产精品三级大全| 好男人在线观看高清免费视频| 国产69精品久久久久777片| 国产av一区在线观看免费| 欧美日韩国产亚洲二区| 亚洲国产高清在线一区二区三| 欧美日本亚洲视频在线播放| 亚洲美女视频黄频| 丰满乱子伦码专区| 久久久久久久精品吃奶| 三级毛片av免费| 精品国内亚洲2022精品成人| 国产亚洲精品久久久久久毛片| 人妻少妇偷人精品九色| 日本一本二区三区精品| 国产精品亚洲一级av第二区| 99久久精品热视频| 露出奶头的视频| 热99re8久久精品国产| 亚洲第一电影网av| 一级黄色大片毛片| 亚洲自拍偷在线| 亚洲精品成人久久久久久| av国产免费在线观看| 亚洲综合色惰| 国产精品一及| 精品午夜福利在线看| 午夜a级毛片| 国产一区二区在线观看日韩| 1000部很黄的大片| 中文字幕av成人在线电影| 成年女人看的毛片在线观看| 国产一级毛片七仙女欲春2| 夜夜看夜夜爽夜夜摸| 国产极品精品免费视频能看的| 九九在线视频观看精品| 国产精品亚洲一级av第二区| 身体一侧抽搐| 国产高潮美女av| 成人av一区二区三区在线看| ponron亚洲| 亚洲真实伦在线观看| 婷婷色综合大香蕉| 在线免费观看的www视频| 亚洲av免费在线观看| 99热网站在线观看| 国产精品,欧美在线| 有码 亚洲区| 精品欧美国产一区二区三| 日本在线视频免费播放| 亚洲成av人片在线播放无| 国产av在哪里看| 露出奶头的视频| 久久久久性生活片| 午夜日韩欧美国产| 婷婷精品国产亚洲av在线| av在线天堂中文字幕| 美女 人体艺术 gogo| 赤兔流量卡办理| 国产三级中文精品| 99在线人妻在线中文字幕| 久久午夜亚洲精品久久| 人妻丰满熟妇av一区二区三区| 日本黄色视频三级网站网址| 久久久久久久久久黄片| 美女cb高潮喷水在线观看| 观看免费一级毛片| 欧美日韩黄片免| 国产精品国产三级国产av玫瑰| 亚洲精品久久国产高清桃花| 蜜桃亚洲精品一区二区三区| 亚洲国产精品合色在线| 在线观看一区二区三区| 免费看美女性在线毛片视频| 无人区码免费观看不卡| 国产伦精品一区二区三区四那| 成人毛片a级毛片在线播放| 变态另类丝袜制服| 国产真实乱freesex| 亚洲中文日韩欧美视频| 97碰自拍视频| 舔av片在线| 精品99又大又爽又粗少妇毛片 | 天美传媒精品一区二区| 欧美xxxx性猛交bbbb| 欧美最黄视频在线播放免费| 国产精品无大码| 午夜免费激情av| 欧美潮喷喷水| 国产三级在线视频| 亚洲aⅴ乱码一区二区在线播放| 两人在一起打扑克的视频| 日韩,欧美,国产一区二区三区 | 男人舔奶头视频| 亚洲无线在线观看| 干丝袜人妻中文字幕| 欧美激情久久久久久爽电影| 九九久久精品国产亚洲av麻豆| 色精品久久人妻99蜜桃| 九九爱精品视频在线观看| 又黄又爽又免费观看的视频| 日日啪夜夜撸| 国产一区二区亚洲精品在线观看| 一进一出抽搐动态| 欧美日韩国产亚洲二区| 久久精品久久久久久噜噜老黄 | 一区二区三区四区激情视频 | 听说在线观看完整版免费高清| 免费看日本二区| 波野结衣二区三区在线| 天堂动漫精品| 很黄的视频免费| 伦理电影大哥的女人| 亚洲综合色惰| av女优亚洲男人天堂| 中文字幕av在线有码专区| 91久久精品电影网| 18禁黄网站禁片免费观看直播| 免费高清视频大片| 美女xxoo啪啪120秒动态图| 午夜福利视频1000在线观看| 久久久久九九精品影院| 亚洲一区高清亚洲精品| 欧美潮喷喷水| 身体一侧抽搐| 波多野结衣高清作品| 国产亚洲精品久久久com| 亚洲天堂国产精品一区在线| 婷婷精品国产亚洲av| 韩国av在线不卡| 中文字幕人妻熟人妻熟丝袜美| 亚洲中文日韩欧美视频| 黄色一级大片看看| 欧美三级亚洲精品| 99久久精品热视频| а√天堂www在线а√下载| 久久久久性生活片| 亚洲va日本ⅴa欧美va伊人久久| 十八禁网站免费在线| 成年女人毛片免费观看观看9| 国产女主播在线喷水免费视频网站 | 18禁裸乳无遮挡免费网站照片| 国产欧美日韩精品一区二区| 中文字幕av成人在线电影| 一区二区三区四区激情视频 | 校园人妻丝袜中文字幕| 精品久久久久久久久久免费视频| 国产男靠女视频免费网站| 午夜老司机福利剧场| 国产亚洲精品综合一区在线观看| 久久精品国产亚洲av香蕉五月| 国产精品美女特级片免费视频播放器| 国产毛片a区久久久久| 久久久久久国产a免费观看| 国内久久婷婷六月综合欲色啪| 神马国产精品三级电影在线观看| videossex国产| 亚洲国产日韩欧美精品在线观看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品成人久久久久久| 国产人妻一区二区三区在| 波野结衣二区三区在线| 日韩欧美免费精品| 亚洲三级黄色毛片| 搡女人真爽免费视频火全软件 | 啦啦啦观看免费观看视频高清| 五月玫瑰六月丁香| 天堂av国产一区二区熟女人妻| 色综合站精品国产| 欧美日本视频| 久久久久久久亚洲中文字幕| 精品国产三级普通话版| 亚洲男人的天堂狠狠| 国产精品av视频在线免费观看| 亚洲男人的天堂狠狠| 久久精品国产亚洲av香蕉五月| 国产av一区在线观看免费| 久久精品国产清高在天天线| 成人美女网站在线观看视频| 久久精品综合一区二区三区| 91狼人影院| 校园人妻丝袜中文字幕| 精品福利观看| 男人的好看免费观看在线视频| 午夜福利在线观看免费完整高清在 | 亚洲综合色惰| 欧美一级a爱片免费观看看| 天堂√8在线中文| 日日夜夜操网爽| 欧美日韩中文字幕国产精品一区二区三区| 中国美白少妇内射xxxbb| 在线观看美女被高潮喷水网站| 天天一区二区日本电影三级| 国产精品久久视频播放| 日本黄大片高清| 可以在线观看的亚洲视频| 嫩草影视91久久| 国产精品嫩草影院av在线观看 | 22中文网久久字幕| 久久久午夜欧美精品| 啪啪无遮挡十八禁网站| 亚洲国产精品久久男人天堂| 成人国产一区最新在线观看| 日韩亚洲欧美综合| 成人国产综合亚洲| 狠狠狠狠99中文字幕| 熟妇人妻久久中文字幕3abv| 久久久久久大精品| 久久久精品大字幕| 亚洲美女视频黄频| 久久精品国产99精品国产亚洲性色| 国产极品精品免费视频能看的| 91久久精品国产一区二区成人| 欧美日本视频| 床上黄色一级片| 国产高清三级在线| 久久这里只有精品中国| 亚洲精品456在线播放app | 国产高清有码在线观看视频| 18+在线观看网站| 日本精品一区二区三区蜜桃| 精品久久久久久成人av| 国产午夜精品久久久久久一区二区三区 | 波多野结衣高清无吗| 一本一本综合久久| 小说图片视频综合网站| 亚洲黑人精品在线| 国产精品99久久久久久久久| a级毛片a级免费在线| 成人一区二区视频在线观看| 97超视频在线观看视频| 看片在线看免费视频| 免费大片18禁| 网址你懂的国产日韩在线| 禁无遮挡网站| 嫩草影视91久久| 91精品国产九色| 国产主播在线观看一区二区| 日韩大尺度精品在线看网址| 黄色视频,在线免费观看| 日韩精品中文字幕看吧| 国产精华一区二区三区| 波多野结衣高清无吗| 精品午夜福利视频在线观看一区| 我的女老师完整版在线观看| 真人做人爱边吃奶动态| 亚洲欧美日韩卡通动漫| 日韩中文字幕欧美一区二区| x7x7x7水蜜桃| 国产精品av视频在线免费观看| 精品久久久久久久久久免费视频| 乱人视频在线观看| 一个人看的www免费观看视频| 亚洲av免费高清在线观看| 波多野结衣高清无吗| 欧美极品一区二区三区四区| 在线免费观看的www视频| 麻豆一二三区av精品| 日本爱情动作片www.在线观看 | 亚洲一级一片aⅴ在线观看| 国产激情偷乱视频一区二区| 亚洲欧美日韩卡通动漫| av.在线天堂| 午夜激情福利司机影院| 国产高清激情床上av| 国产大屁股一区二区在线视频| 久久精品国产亚洲av香蕉五月| 九九爱精品视频在线观看| 国产欧美日韩精品亚洲av| 成人永久免费在线观看视频| 欧美成人一区二区免费高清观看| 日韩一本色道免费dvd| 亚洲精品456在线播放app | 国产伦精品一区二区三区四那| 免费大片18禁| 日韩中文字幕欧美一区二区| 麻豆精品久久久久久蜜桃| 男女下面进入的视频免费午夜| 成人美女网站在线观看视频| 免费大片18禁| 日本黄色视频三级网站网址| 成年免费大片在线观看| 久久中文看片网| 91在线观看av| 最后的刺客免费高清国语| 国产三级中文精品| 国内久久婷婷六月综合欲色啪| 精品久久国产蜜桃| 老熟妇乱子伦视频在线观看| 国产人妻一区二区三区在| 亚洲精品成人久久久久久| 成人美女网站在线观看视频| 午夜亚洲福利在线播放| 999久久久精品免费观看国产| 美女高潮的动态| 久久精品国产清高在天天线| 91在线观看av| 久久久久九九精品影院| 午夜影院日韩av| 老女人水多毛片| 99九九线精品视频在线观看视频| 在线天堂最新版资源| 麻豆精品久久久久久蜜桃| 香蕉av资源在线| 亚洲午夜理论影院| 久久久久国内视频| 人人妻人人澡欧美一区二区| av黄色大香蕉| 精品久久久久久久久av| 亚洲国产精品合色在线| 国产人妻一区二区三区在| 3wmmmm亚洲av在线观看| 啦啦啦韩国在线观看视频| 欧美成人性av电影在线观看| 日韩欧美 国产精品| 一卡2卡三卡四卡精品乱码亚洲| 变态另类丝袜制服| 免费无遮挡裸体视频| 18禁黄网站禁片午夜丰满| 久久这里只有精品中国| 免费av观看视频| www.www免费av| 国产三级在线视频| 国产在线精品亚洲第一网站| 我要搜黄色片| 麻豆久久精品国产亚洲av| 久久草成人影院| 国产av不卡久久| 国产毛片a区久久久久| 人妻少妇偷人精品九色| 男女做爰动态图高潮gif福利片| 久久久久久久久中文| 色精品久久人妻99蜜桃| 日本五十路高清| av在线老鸭窝| 亚洲七黄色美女视频| ponron亚洲| 直男gayav资源| 亚洲性夜色夜夜综合| 欧美三级亚洲精品| 特大巨黑吊av在线直播| 国国产精品蜜臀av免费| 国产老妇女一区| 欧美丝袜亚洲另类 | 真人做人爱边吃奶动态| videossex国产| 亚洲精品一区av在线观看| 老女人水多毛片| 嫩草影院新地址| 国内少妇人妻偷人精品xxx网站| 亚洲图色成人| 国产精品精品国产色婷婷| 国产午夜福利久久久久久| 成熟少妇高潮喷水视频| 综合色av麻豆| 久久精品91蜜桃| 欧美极品一区二区三区四区| 91久久精品国产一区二区三区| 中文亚洲av片在线观看爽| 欧美丝袜亚洲另类 | 亚洲中文日韩欧美视频| 欧美激情在线99| 色哟哟哟哟哟哟| xxxwww97欧美| av国产免费在线观看| 级片在线观看| 国产精品日韩av在线免费观看| 伊人久久精品亚洲午夜| 久99久视频精品免费| 日韩欧美一区二区三区在线观看| 国产伦一二天堂av在线观看| 麻豆成人午夜福利视频| 日韩欧美一区二区三区在线观看| 欧美人与善性xxx| 午夜爱爱视频在线播放| 久久久久久久久大av| 免费看光身美女| 可以在线观看的亚洲视频| 99久久中文字幕三级久久日本| 国产精品乱码一区二三区的特点| 18禁黄网站禁片午夜丰满| 我的女老师完整版在线观看| 看免费成人av毛片| 国产成人av教育| 国产精品人妻久久久影院| 国产毛片a区久久久久| 国内精品久久久久久久电影| 五月玫瑰六月丁香| 亚洲精品一卡2卡三卡4卡5卡| 久久久国产成人精品二区| 国产精品三级大全|