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

    基于CFD的汽車前端結(jié)構(gòu)匹配研究

    2017-11-06 07:24:58孫躍東吳旭陵陸飛龍于卓弘
    汽車工程學(xué)報(bào) 2017年5期
    關(guān)鍵詞:風(fēng)洞試驗(yàn)散熱器整車

    劉 念,孫躍東,吳旭陵,陸飛龍,于卓弘

    (1.上海理工大學(xué) 機(jī)械工程學(xué)院,上海 200093;2.上海汽車集團(tuán)股份有限公司 商用車技術(shù)中心,上海 200438)

    基于CFD的汽車前端結(jié)構(gòu)匹配研究

    劉 念1,孫躍東1,吳旭陵2,陸飛龍2,于卓弘2

    (1.上海理工大學(xué) 機(jī)械工程學(xué)院,上海 200093;2.上海汽車集團(tuán)股份有限公司 商用車技術(shù)中心,上海 200438)

    采用CFD仿真技術(shù)結(jié)合正交試驗(yàn)方法對(duì)主要前端結(jié)構(gòu)進(jìn)行了匹配研究。基于Fluent平臺(tái),搭建了整車的前端結(jié)構(gòu)分析模型,設(shè)計(jì)了9種正交試驗(yàn)方案,系統(tǒng)地研究了前端格柵開口、冷凝器導(dǎo)流板、冷卻風(fēng)扇直徑對(duì)整車前端進(jìn)氣、氣動(dòng)阻力的影響。通過整車風(fēng)洞試驗(yàn)及熱管理環(huán)境艙試驗(yàn)對(duì)最優(yōu)配置進(jìn)行了驗(yàn)證,相對(duì)誤差值約2.28%。改進(jìn)方案的Cd值為0.421 06,比改進(jìn)前降低3%,且動(dòng)力總成冷卻系統(tǒng)(PTC)性能良好。研究結(jié)果對(duì)商用車前端結(jié)構(gòu)的設(shè)計(jì)和匹配具有指導(dǎo)意義。

    前端結(jié)構(gòu)匹配;正交試驗(yàn);CFD仿真;整車試驗(yàn)

    引用格式:

    前端配置包括進(jìn)氣格柵開口、冷卻模塊(冷凝器、散熱器、中冷器)和風(fēng)扇總成等結(jié)構(gòu)部件。目前,國內(nèi)外學(xué)者對(duì)前端結(jié)構(gòu)的研究主要集中在通風(fēng)散熱方面[1-4]。WATANABE等[1]基于一維與三維耦合仿真模擬手段,分析了前端結(jié)構(gòu)對(duì)汽車發(fā)動(dòng)機(jī)艙熱管理的影響。COSTA[2]采用CFD軟件研究了進(jìn)氣格柵開口大小、形狀等幾何因素對(duì)發(fā)動(dòng)機(jī)艙內(nèi)部散熱性能的影響。通用汽車的楊志剛等[3]基于三維仿真軟件STAR-CD,通過對(duì)比搭載傳統(tǒng)和全新前端冷卻模塊的發(fā)動(dòng)機(jī)艙內(nèi)部速度場(chǎng)、溫度場(chǎng),結(jié)合實(shí)車試驗(yàn)的驗(yàn)證,提出了全新的汽車?yán)鋮s模塊概念。哈飛汽車的李玉平[4]通過熱管理環(huán)境艙試驗(yàn)研究了前端進(jìn)氣格柵、散熱器以及冷卻風(fēng)扇單體設(shè)計(jì)對(duì)發(fā)動(dòng)機(jī)散熱性能的影響。眾多研究成果加深了人們對(duì)前端結(jié)構(gòu)的理解,也推動(dòng)了前端結(jié)構(gòu)匹配問題的研究。

    在整車開發(fā)流程中,為縮短開發(fā)周期,降低研發(fā)成本,前端結(jié)構(gòu)的匹配主要利用CFD仿真技術(shù)進(jìn)行研究改進(jìn)??紤]到前端配置的匹配需要兼顧動(dòng)力總成冷卻系統(tǒng)(PTC)性能和整車空氣動(dòng)力學(xué)性能。因此,本研究針對(duì)某企業(yè)的皮卡車型,基于Fluent平臺(tái)搭建有效的仿真模型,借助于正交試驗(yàn)來探究不同前端結(jié)構(gòu)因素對(duì)整車前端進(jìn)氣、氣動(dòng)阻力的影響規(guī)律,以此匹配最優(yōu)的前端結(jié)構(gòu)配置方案。

    1 前端結(jié)構(gòu)仿真分析基礎(chǔ)

    1.1 基本方程組和湍動(dòng)能模型

    基于Fluent軟件的前端進(jìn)氣分析和氣動(dòng)阻力分析,涉及了汽車的內(nèi)外流場(chǎng)耦合的過程,需要滿足下列流動(dòng)與傳熱的基本方程組[5]。

    1.1.1 質(zhì)量守恒方程

    將空氣看作不可壓縮氣體,ρu為常數(shù),取值為1.225,單位為kg/m3。

    1.1.2 動(dòng)量守恒方程

    式中:P為流體壓力;u為液體粘度;S為廣義源項(xiàng);i和j為張量指標(biāo),取值范圍為1,2,3。

    1.1.3 能量守恒方程

    式中:keff為有效熱傳導(dǎo)系數(shù);h為焓;J為擴(kuò)散流量;Sh為熱源項(xiàng)。

    1.1.4 可實(shí)現(xiàn)的湍動(dòng)能方程

    湍動(dòng)能模型選用Realizable k-ε方程,其計(jì)算精度高,適用范圍廣,可以更為精準(zhǔn)地模擬汽車的復(fù)雜流場(chǎng)分布情況[6]。

    湍動(dòng)能k方程:

    湍動(dòng)能耗散率ε方程:

    其中C1、η、S的表達(dá)式如下:

    式(4)~(5)中:ρu為密度,kg/m3;Ui為i方向的速度,m/s;Gk為平均速度梯度產(chǎn)生的湍動(dòng)能,m2/s2;Gb為浮力產(chǎn)生的湍動(dòng)能,m2/s2;Sk、Sε為用戶定義項(xiàng);C1ε、C2ε、C3ε為常數(shù),設(shè)定為1.44、1.9、0.8;為K方程的Prandtl數(shù); 為ε方程的Prandtl數(shù)。

    1.2 離散方法和格式

    由于有限體積法對(duì)網(wǎng)格的精度、數(shù)量要求不高,且積分守恒精準(zhǔn),所以通過該離散方法對(duì)計(jì)算區(qū)域內(nèi)的所有控制體積單元進(jìn)行積分,將原先的連續(xù)性微分方程轉(zhuǎn)換成離散方程組[7]。

    一般情況下,對(duì)于面網(wǎng)格采用三角形網(wǎng)格、體網(wǎng)格采用四面體網(wǎng)格的復(fù)雜流動(dòng)問題,通常需要采用二階精度的離散格式,因此采用二階迎風(fēng)格式進(jìn)行插值離散求解。

    二階迎風(fēng)格式的對(duì)流-擴(kuò)散離散方程為:

    式中各參數(shù)用下式表達(dá):

    其中,流動(dòng)方向?yàn)檎?1;流動(dòng)方向?yàn)樨?fù),α=10。

    2 建立仿真模型

    2.1 網(wǎng)格劃分

    對(duì)含有發(fā)動(dòng)機(jī)艙和底盤件的1∶1整車模型進(jìn)行幾何清理,局部特征網(wǎng)格進(jìn)行細(xì)化,網(wǎng)格尺寸為3~5 mm,整車總體區(qū)域面網(wǎng)格尺寸為10~15 mm,面網(wǎng)格形式為三角形網(wǎng)格,網(wǎng)格總數(shù)為5 464 634個(gè)。軟件窗口界面中shaded-elements狀態(tài)的整車模型,如圖1所示。

    圖1 整車數(shù)模

    為了保證整車內(nèi)、外流場(chǎng)能夠充分延展,將整車置于計(jì)算域中間位置,整車的長(zhǎng)、寬、高為L(zhǎng)、W、H,計(jì)算域的長(zhǎng)度為10 L,寬度為8 W,高度為5.5 H,側(cè)面距離汽車邊界為3.5 W。將質(zhì)量檢查合格的面網(wǎng)格模型直接導(dǎo)入T-Grid軟件,并生成計(jì)算區(qū)域和整車的體網(wǎng)格。體網(wǎng)格形式主要為四面體網(wǎng)格,多孔介質(zhì)區(qū)域體網(wǎng)格形式為六面體網(wǎng)格,計(jì)算網(wǎng)格總數(shù)為42 672 845個(gè),整車附近流場(chǎng)區(qū)域進(jìn)行網(wǎng)格加密處理。截取αY =10截面的計(jì)算域網(wǎng)格,如圖2所示。

    圖2 體網(wǎng)格處理

    2.2 邊界條件設(shè)置

    2.2.1 多孔介質(zhì)

    在Fluent軟件的條件設(shè)置中,冷凝器、中冷器、散熱器的簡(jiǎn)化數(shù)學(xué)模型和多孔介質(zhì)的數(shù)學(xué)模型相同,即通過在動(dòng)量方程中增加粘性阻力項(xiàng)和慣性損失項(xiàng)來表述流體計(jì)算域中換熱器的多孔性材料對(duì)前端進(jìn)氣流體的流動(dòng)阻力[8]。

    將冷卻模塊定義為簡(jiǎn)易的均勻多孔介質(zhì),Darcy動(dòng)量源項(xiàng)簡(jiǎn)化為式(7),多孔介質(zhì)參數(shù)的具體設(shè)置情況見表1。

    式中:Si為多孔介質(zhì)元件的壓降;μ為粘性系數(shù);Pi為慣性阻力系數(shù);ρu為空氣的密度;V為多孔介質(zhì)表面的垂直速度分量;Pv為粘性阻力系數(shù)。

    表1 多孔介質(zhì)參數(shù)設(shè)置

    2.2.2 冷卻風(fēng)扇MRF模型

    針對(duì)冷卻風(fēng)扇選擇MRF模型,建立獨(dú)立的風(fēng)扇旋轉(zhuǎn)區(qū)域,旋轉(zhuǎn)中心坐標(biāo)為(1 242.7,1.646,237.185),旋轉(zhuǎn)軸矢量為(1,0,0),旋轉(zhuǎn)速度為3 000 r/min,網(wǎng)格模型如圖3所示。

    2.2.3 邊界條件

    使用數(shù)值模擬的方法對(duì)整車數(shù)模進(jìn)行前端進(jìn)氣和空氣動(dòng)力學(xué)分析,實(shí)際是在模擬氣體經(jīng)過汽車的流動(dòng)過程,具體的模型和參數(shù)設(shè)置見表2。

    圖3 冷卻風(fēng)扇網(wǎng)格模型

    表2 邊界條件設(shè)置

    工程設(shè)計(jì)上,動(dòng)力總成冷卻系統(tǒng)(PTC)性能的前期開發(fā)是參考汽車前端進(jìn)氣狀態(tài)來設(shè)計(jì)考慮的。按照整車設(shè)計(jì)規(guī)范,前端進(jìn)氣的仿真計(jì)算工況為惡劣工況(4擋、80 km/h、爬坡)的散熱器進(jìn)風(fēng)量,氣動(dòng)阻力的仿真計(jì)算工況為車速120 km/h(國內(nèi)的高速路段行駛速度上限)的氣動(dòng)阻力系數(shù)。

    3 正交試驗(yàn)設(shè)計(jì)

    正交試驗(yàn)方法能夠根據(jù)正交表選取部分試驗(yàn)代替全因子試驗(yàn),從而精簡(jiǎn)試驗(yàn)次數(shù),提高試驗(yàn)效率。

    3.1 因素及水平選取

    整車前端進(jìn)氣與氣動(dòng)阻力的主要前端配置影響因素很多,為了避免造成試驗(yàn)矩陣過大,不適宜全部選取。由于整車開發(fā)過程中,冷卻模塊布置形式已經(jīng)確定為中冷器-冷凝器-散熱器。因此,根據(jù)工程經(jīng)驗(yàn)和工程技術(shù)成本顧慮,不再考察冷卻模塊布置形式、格柵形狀的影響。

    基于上述考慮,從車輛燃油經(jīng)濟(jì)性理論分析的角度出發(fā),提出了影響試驗(yàn)結(jié)果的3個(gè)影響因素:

    (1)格柵開口率A,即格柵的開口面積到散熱器表面積的正投影率,如圖4a所示。

    (2)冷凝器導(dǎo)流板角度B,即冷凝器下端導(dǎo)流板與z軸的夾角,如圖4b所示。

    (3)冷卻風(fēng)扇直徑C,即風(fēng)扇的電機(jī)直徑,如圖4c所示。

    圖4 前端配置因素示意圖

    前端配置的因素和水平確定后,將利用正交表進(jìn)行3因素3水平的L9(33)正交試驗(yàn),其試驗(yàn)因素水平安排見表3。

    3.2 正交試驗(yàn)方案

    由上文設(shè)計(jì)的試驗(yàn)因素及水平表可知,這是一個(gè)3因素3水平的正交試驗(yàn),并且不考慮各因素之間的交互作用,因此可以選擇L9(33)正交試驗(yàn)矩陣表安排試驗(yàn),正交試驗(yàn)表及試驗(yàn)結(jié)果見表4。

    表3 前端配置試驗(yàn)因素水平表

    表4 前端配置的正交試驗(yàn)矩陣

    3.3 仿真計(jì)算結(jié)果

    采用Fluent軟件對(duì)上述前端配置的試驗(yàn)方案進(jìn)行前端進(jìn)氣及氣動(dòng)阻力仿真分析,得到各個(gè)方案的響應(yīng)目標(biāo)值,見表5。優(yōu)化前前端結(jié)構(gòu)因素、水平配置為A3B2C3,仿真分析得出風(fēng)阻系數(shù)Cd為0.419 810 4,散熱器芯體進(jìn)風(fēng)量為57.73 CMM。

    表5 正交試驗(yàn)結(jié)果表

    注:Cd表示整車風(fēng)阻系數(shù);Q表示散熱器芯體進(jìn)風(fēng)量,均為正交試驗(yàn)設(shè)計(jì)的響應(yīng)指標(biāo)。PTC性能設(shè)計(jì)要求:4擋、80 km/h、爬坡工況下的散熱器芯體進(jìn)風(fēng)量不可低于45 CMM。

    3.4 數(shù)據(jù)分析

    3.4.1 極差分析

    為了找出各因素對(duì)各響應(yīng)指標(biāo)的影響規(guī)律,采用極差分析法對(duì)數(shù)據(jù)進(jìn)行分析。極差是一組數(shù)據(jù)中最大值和最小值之差,極差值越大表示該因素對(duì)響應(yīng)指標(biāo)的影響越大,該因素越關(guān)鍵;極差值越小表示該因素對(duì)響應(yīng)指標(biāo)的影響越小,該因素重要度越低。表6和表7為各因素不同水平對(duì)各響應(yīng)指標(biāo)的極差分析結(jié)果。

    表6 響應(yīng)目標(biāo)Cd極差分析

    表7 響應(yīng)目標(biāo)-Q極差分析

    比較不同列的極差大小可以得到選定因素對(duì)響應(yīng)目標(biāo)的影響程度大小。Rj值大表明該因素j對(duì)響應(yīng)目標(biāo)影響大,反之則表明該因素對(duì)響應(yīng)目標(biāo)影響小。根據(jù)極差分析結(jié)果,得出各響應(yīng)指標(biāo)的響應(yīng)均值圖,如圖5和圖6所示。

    圖5 因素水平-響應(yīng)均值Cd分布

    圖6 因素水平-響應(yīng)均值Q分布

    由圖5可知,整車風(fēng)阻系數(shù)Cd隨格柵開口率A增大而增大,隨冷卻風(fēng)扇直徑C增大而減小。由于響應(yīng)指標(biāo)Cd是望小最優(yōu),因此各因素影響強(qiáng)弱順序?yàn)锳>B>C,且最優(yōu)水平分別為14.8%(A1)、+15(B3)、480 mm(C3)。

    由圖6可知,散熱器芯體進(jìn)風(fēng)量Q隨格柵開口率A、冷卻風(fēng)扇直徑C增大而增大。由于響應(yīng)指標(biāo)Q是望大最優(yōu),因此各因素影響強(qiáng)弱順序?yàn)锳>B>C,且最優(yōu)水平分別為24.8%(A3)、-15(B1)、480 mm(C3)。

    兼顧前端進(jìn)氣和氣動(dòng)阻力性能匹配最優(yōu)的前端結(jié)構(gòu)配置,應(yīng)先確保前端進(jìn)氣量滿足設(shè)計(jì)要求,再選擇氣動(dòng)阻力性能最優(yōu)的配置形式,以此節(jié)約整車能耗。綜上所述,方案3(A1B3C3)是最佳的前端結(jié)構(gòu)匹配方案。

    3.4.2 速度云圖分析

    發(fā)動(dòng)機(jī)艙散熱狀況和整車氣動(dòng)阻力深受氣流的流動(dòng)狀態(tài)影響,所以研究發(fā)動(dòng)機(jī)艙內(nèi)的空氣流動(dòng)特性十分必要。CFD仿真分析可以得到前端配置改進(jìn)前后αy =10截面的速度云圖,如圖7所示。

    圖7 速度云圖

    由圖7a可知,改進(jìn)前前端結(jié)構(gòu)的前端進(jìn)氣量大,內(nèi)流阻力較大,氣動(dòng)性能較差。格柵開口率大,外界氣流會(huì)快速流入進(jìn)氣格柵。由于冷凝器下部導(dǎo)流板與z軸平行,氣流經(jīng)過時(shí)沒有經(jīng)過任何導(dǎo)流的作用,會(huì)較快地通過冷卻模塊。由圖7b可知,改進(jìn)后空氣流入前端進(jìn)氣格柵的速度明顯減小,經(jīng)過冷凝器下部導(dǎo)流板時(shí)速度變化梯度平緩,且冷凝器芯體迎風(fēng)面區(qū)域流場(chǎng)速度較優(yōu)化前變低,整個(gè)前端的流場(chǎng)得到較大的改善,內(nèi)部阻力損失有效降低。雖然前端進(jìn)氣量較優(yōu)化前有所降低,但是散熱器芯體表面進(jìn)風(fēng)風(fēng)量滿足設(shè)計(jì)規(guī)范要求,發(fā)動(dòng)機(jī)冷卻性能穩(wěn)定。

    4 試驗(yàn)驗(yàn)證

    4.1 整車風(fēng)洞試驗(yàn)

    整車風(fēng)洞試驗(yàn)設(shè)備主要包括:上海地面交通工具風(fēng)洞中心的五帶移動(dòng)地面系統(tǒng)、六分量測(cè)試天平、實(shí)車模型和計(jì)算機(jī)數(shù)據(jù)采集系統(tǒng)。參照整車風(fēng)洞試驗(yàn)規(guī)范,將整車固定于測(cè)試地面中心,進(jìn)行整車氣動(dòng)阻力風(fēng)洞試驗(yàn),如圖8所示。調(diào)整整車姿態(tài),測(cè)得整車正投影面積為2.940 2 m2。

    圖8 氣動(dòng)阻力風(fēng)洞試驗(yàn)

    設(shè)定整車風(fēng)洞試驗(yàn)考核工況為120 km/h下,風(fēng)速穩(wěn)定狀態(tài)的整車風(fēng)阻系數(shù)。試驗(yàn)全過程,風(fēng)機(jī)工作帶動(dòng)氣流從氣動(dòng)聲學(xué)風(fēng)洞噴口流出并吹向試驗(yàn)車輛,風(fēng)速由0 逐漸增加至120 km/h,待風(fēng)速穩(wěn)定在120 km/h狀態(tài),計(jì)算機(jī)數(shù)據(jù)采集系統(tǒng)會(huì)自動(dòng)采集該工況下的整車風(fēng)阻系數(shù)。數(shù)據(jù)采集完畢,風(fēng)機(jī)停止工作,車輛前端氣流停止流動(dòng)。

    對(duì)比響應(yīng)目標(biāo)Cd的最優(yōu)前端配置方案的仿真分析結(jié)果和整車風(fēng)洞試驗(yàn)結(jié)果,見表8。試驗(yàn)結(jié)果表明,改進(jìn)前后的仿真數(shù)值與試驗(yàn)數(shù)值變化趨勢(shì)相同,仿真誤差平均值約3.36%,誤差較小。CFD氣動(dòng)阻力分析可以使用于前端配置風(fēng)阻優(yōu)化的前期設(shè)計(jì)階段。

    表8 CFD仿真與風(fēng)洞試驗(yàn)結(jié)果

    4.2 整車熱管理環(huán)境艙試驗(yàn)

    整車熱管理環(huán)境艙試驗(yàn)設(shè)備主要包括:熱電偶、數(shù)據(jù)采集系統(tǒng)和溫度傳感器等裝置。試驗(yàn)依據(jù)發(fā)動(dòng)機(jī)冷卻系統(tǒng)的設(shè)計(jì)目標(biāo),對(duì)整車發(fā)動(dòng)機(jī)艙內(nèi)所關(guān)注的部件進(jìn)行布點(diǎn),如圖9所示。

    圖9 發(fā)動(dòng)機(jī)艙內(nèi)布點(diǎn)示意圖

    根據(jù)商用車熱管理試驗(yàn)法規(guī),在某技術(shù)中心的環(huán)境艙內(nèi)對(duì)整車進(jìn)行完整的城市、城郊、4擋50 km/h爬坡7.2%以及4擋80 km/h爬坡7.2%等工況熱管理測(cè)試。由于本研究在熱管理試驗(yàn)中僅關(guān)注發(fā)動(dòng)機(jī)冷卻系統(tǒng)能力是否滿足設(shè)計(jì)目標(biāo),因此重點(diǎn)關(guān)注惡劣考核工況的試驗(yàn)數(shù)據(jù)。

    依據(jù)數(shù)采系統(tǒng)的輸出數(shù)據(jù),整理惡劣考核工況下各個(gè)散熱部件的測(cè)點(diǎn)溫度,見表9。4擋、80 km/h、7.2%爬坡工況下的發(fā)動(dòng)機(jī)出水溫度低于115℃,油底殼機(jī)油溫度均低于135℃,變速器油溫均低于132℃,滿足發(fā)動(dòng)機(jī)冷卻系統(tǒng)設(shè)計(jì)要求。最優(yōu)前端配置方案的發(fā)動(dòng)機(jī)系統(tǒng)冷卻性能良好,方案可行性高。

    表9 熱管理試驗(yàn)結(jié)果

    5 結(jié)論

    本文為了匹配合理的前端結(jié)構(gòu)方案,通過仿真與整車試驗(yàn)相結(jié)合的方法對(duì)某皮卡的前端結(jié)構(gòu)進(jìn)行了匹配優(yōu)化。主要結(jié)論如下:

    (1)利用Fluent平臺(tái)對(duì)氣動(dòng)阻力特性及前端進(jìn)氣進(jìn)行了仿真研究,得出了前端結(jié)構(gòu)因素的影響規(guī)律,格柵開口率對(duì)氣動(dòng)阻力和前端進(jìn)氣影響最大,冷卻風(fēng)扇直徑對(duì)氣動(dòng)阻力影響最弱,冷凝器導(dǎo)流板角度對(duì)前端進(jìn)氣影響最弱。

    (2)仿真數(shù)據(jù)分析表明,最優(yōu)前端結(jié)構(gòu)配置為A1B3C3,風(fēng)阻系數(shù)為0.410 764 9,散熱器表面進(jìn)風(fēng)量為46.92 CMM。

    (3)在整車風(fēng)洞試驗(yàn)中,測(cè)得前端結(jié)構(gòu)改進(jìn)后的整車風(fēng)阻系數(shù)為0.421 06,比改進(jìn)前降低3%,相對(duì)誤差約2.82%,誤差較小。在整車熱管理環(huán)境艙試驗(yàn)中,測(cè)得考核工況下的發(fā)動(dòng)機(jī)出水溫度、油底殼機(jī)油溫度、變速器油溫均滿足設(shè)計(jì)要求,改進(jìn)方案的動(dòng)力總成冷卻系統(tǒng)(PTC)性能良好。

    參考文獻(xiàn)(References):

    [1]WATANABEN,KUBOM,YOMODAN. An 1D-3D Integrating Numerical Simulation for Engine Cooling Problem[C]//SAE Technical Papers,2006-01-1603,2006.

    [2]COSTAE A. CFD Approach on Underhood Thermal Management of Passenger Cars and Trucks [C]//SAE Technical Papers,2003-01-3577,2003.

    [3]YANG Zhigang,BOZEMANJ,SHENF Z,et al. CFRM Concept for Vehicle Thermal System[C]//SAE Technical Papers,2002-01-1207,2002.

    [4]李玉平. HF5整車熱平衡研究 [D]. 長(zhǎng)春:吉林大學(xué),2006.LIYuping. Study on HF5 Vehicle Thermal Balance [D].Changchun:Jilin University,2006. (in Chinese)

    [5]李進(jìn)良. 精通FLUENT6.3流場(chǎng)分析 [M]. 北京:化學(xué)工業(yè)出版社,2009.LI Jinliang. Pro fi cient in FLUENT6.3 Flow Field Analysis[M]. Beijing:ChemicalIndustry Press,2009. (in Chinese)

    [6]王福軍.計(jì)算流體動(dòng)力學(xué)分析[M].北京:清華大學(xué)出版社,2004.WANG Fujun. Computational Fluid Dynamics Analysis [M].Beijing:Tsinghua University Press,2004. (in Chinese)

    [7]張英朝. 汽車空氣動(dòng)力學(xué)數(shù)值模擬技術(shù) [M]. 北京:北京大學(xué)出版社,2011.ZHANG Yingchao. Automotive Aerodynamics Numerical Simulation Technology[M]. Beijing:Peking University Press,2011. (in Chinese)

    [8]沈俊,傅立敏,范士杰. CFD在汽車空氣動(dòng)力學(xué)設(shè)計(jì)中的應(yīng)用[J]. 汽車技術(shù),2000(10):1-4.SHEN Jun,F(xiàn)u Limin,F(xiàn)an Shijie. Application of CFD in Automotive Aerodynamics Design[J]. Automotive Technology,2000(10):1-4. (in Chinese)

    Research on Matching of Automobile Front-end Structures Based on CFD

    LIU Nian1,SUN Yuedong1,WU Xuling2,LU Feilong2,YU Zhuohong2
    (1. College of Mechanical Engineering, University of Shanghai for Science and Technology,Shanghai 200093,China;2. SAIC Motor Commercial Vehicle Technical Center,Shanghai 200438,China)

    The paper investigated the matching of the main front-end structure by using CFD simulation technology and the orthogonal test method.The simulation model of vehicle front-end structure was built based on Fluent software. Nine orthogonal test schemes were designed to study the effects of the front grille opening, the guide plate of condenser and the cooling fan diameter on the front-end intake and aerodynamic resistance of the automobile. The optimal configuration was verified by the wind tunnel and the thermal management environmental chamber tests, and the simulation error was about 2.28%. The drag coefficient of the improved scheme is 0.421 06, which is about 3% lower than before the improvement,and the good performance of the powertrain in cooling (PTC) system is achieved.The paper is instructive for design and matching of commercial front-end structures.

    front-end structure matching; orthogonal test; CFD simulation; vehicle test

    U461.1;U467.1+3

    A

    10.3969/j.issn.2095-1469.2017.05.05

    2017-04-24 改稿日期:2017-05-29

    上海市研究生創(chuàng)新基金(JWCXSL 1022)

    劉念,孫躍東,吳旭陵,等. 基于CFD的汽車前端結(jié)構(gòu)匹配研究 [J]. 汽車工程學(xué)報(bào),2017,7(5):342-349.

    LIU Nian,SUN Yuedong,WU Xuling,et al. Research on Matching of Automobile Front-end Structures Based on CFD [J]. Chinese Journal of Automotive Engineering,2017,7(5):342-349. (in Chinese)

    作者介紹

    責(zé)任作者:孫躍東(1965-),男,江蘇興化人。教授,博士生導(dǎo)師,主要研究方向?yàn)檐囕v動(dòng)力系統(tǒng)及匹配技術(shù)、車輛發(fā)動(dòng)機(jī)性能測(cè)試。

    E-mail:syd@usst.edu.cn

    劉念(1992-),女,安徽六安人。碩士研究生,主要研究方向?yàn)槠嚳諝鈩?dòng)力學(xué)、汽車熱管理技術(shù)。

    Tel:13301898576

    E-mail:1457737023@qq.com

    猜你喜歡
    風(fēng)洞試驗(yàn)散熱器整車
    基于六自由度解耦分析的整車懸置設(shè)計(jì)
    ◆ 散熱器
    散熱器
    ◆ 散熱器
    散熱器
    低風(fēng)壓架空導(dǎo)線的風(fēng)洞試驗(yàn)
    電線電纜(2017年5期)2017-10-18 00:52:03
    滾轉(zhuǎn)機(jī)動(dòng)載荷減緩風(fēng)洞試驗(yàn)
    整車低頻加速噪聲研究及改進(jìn)
    遮擋條件下超高層建筑風(fēng)洞試驗(yàn)研究
    重慶建筑(2014年12期)2014-07-24 14:00:32
    高速鐵路接觸線覆冰后氣動(dòng)力特性的風(fēng)洞試驗(yàn)研究
    99热这里只有精品一区 | 国产高清视频在线播放一区| 国产精品,欧美在线| 偷拍熟女少妇极品色| 亚洲成av人片免费观看| 国内精品久久久久久久电影| 亚洲第一电影网av| 亚洲欧美精品综合久久99| 精品午夜福利视频在线观看一区| 午夜福利成人在线免费观看| 91老司机精品| 99国产精品99久久久久| 99国产精品一区二区三区| 欧美激情久久久久久爽电影| 又紧又爽又黄一区二区| 在线观看66精品国产| 精品国产亚洲在线| 日本黄色片子视频| 亚洲av电影不卡..在线观看| 窝窝影院91人妻| h日本视频在线播放| 男人舔奶头视频| 18禁黄网站禁片免费观看直播| 最近视频中文字幕2019在线8| 日韩高清综合在线| 男女那种视频在线观看| 动漫黄色视频在线观看| 午夜福利欧美成人| 亚洲aⅴ乱码一区二区在线播放| 成年免费大片在线观看| 成人高潮视频无遮挡免费网站| 人人妻,人人澡人人爽秒播| 亚洲精品美女久久av网站| 十八禁人妻一区二区| 91麻豆av在线| 亚洲中文日韩欧美视频| www.自偷自拍.com| 九九热线精品视视频播放| 91在线精品国自产拍蜜月 | 淫妇啪啪啪对白视频| 国产亚洲精品久久久久久毛片| 久久精品夜夜夜夜夜久久蜜豆| 搡老妇女老女人老熟妇| 99久国产av精品| 一本综合久久免费| 欧美在线黄色| 99国产综合亚洲精品| 亚洲欧美激情综合另类| 欧美日韩国产亚洲二区| 好男人电影高清在线观看| 九色成人免费人妻av| 香蕉国产在线看| 国产高清三级在线| 免费在线观看影片大全网站| 国产av不卡久久| 久久久色成人| 欧美最黄视频在线播放免费| 一区二区三区高清视频在线| 国产私拍福利视频在线观看| 两性夫妻黄色片| 久久午夜综合久久蜜桃| 亚洲 欧美 日韩 在线 免费| 久久久久久国产a免费观看| av女优亚洲男人天堂 | 欧美性猛交黑人性爽| 久久久色成人| 国产亚洲精品综合一区在线观看| 18禁黄网站禁片免费观看直播| 久久久久久国产a免费观看| 两性夫妻黄色片| 美女高潮的动态| 欧美极品一区二区三区四区| 中文字幕最新亚洲高清| 很黄的视频免费| e午夜精品久久久久久久| 久久性视频一级片| 欧美黄色淫秽网站| 制服人妻中文乱码| 欧美乱码精品一区二区三区| 欧美日韩综合久久久久久 | 成人三级黄色视频| 日韩三级视频一区二区三区| 国产伦人伦偷精品视频| 国产成人影院久久av| 欧美乱妇无乱码| 可以在线观看的亚洲视频| 精品一区二区三区四区五区乱码| 亚洲成a人片在线一区二区| 国产人伦9x9x在线观看| 午夜日韩欧美国产| av天堂中文字幕网| 神马国产精品三级电影在线观看| 国产主播在线观看一区二区| 国产午夜精品久久久久久| 亚洲av熟女| 成人性生交大片免费视频hd| 99久国产av精品| 怎么达到女性高潮| 老司机在亚洲福利影院| 日韩欧美国产一区二区入口| 99国产综合亚洲精品| 亚洲第一欧美日韩一区二区三区| 黑人巨大精品欧美一区二区mp4| 精品乱码久久久久久99久播| 久9热在线精品视频| 国产伦精品一区二区三区四那| 亚洲精品在线观看二区| 久久中文字幕一级| 国内久久婷婷六月综合欲色啪| 黑人巨大精品欧美一区二区mp4| 美女高潮喷水抽搐中文字幕| 国产日本99.免费观看| 成人永久免费在线观看视频| 国产蜜桃级精品一区二区三区| 高清毛片免费观看视频网站| 一区二区三区国产精品乱码| 少妇丰满av| 成人无遮挡网站| 亚洲avbb在线观看| 在线十欧美十亚洲十日本专区| 亚洲精品一区av在线观看| 91av网站免费观看| 欧美国产日韩亚洲一区| 日韩免费av在线播放| 噜噜噜噜噜久久久久久91| 欧美3d第一页| 在线观看66精品国产| 亚洲国产欧洲综合997久久,| 日本黄色片子视频| 夜夜看夜夜爽夜夜摸| 欧美一区二区国产精品久久精品| 最新在线观看一区二区三区| 在线a可以看的网站| 他把我摸到了高潮在线观看| 亚洲av日韩精品久久久久久密| 欧美成人免费av一区二区三区| 黄色片一级片一级黄色片| 亚洲成人中文字幕在线播放| 哪里可以看免费的av片| 一个人观看的视频www高清免费观看 | 婷婷精品国产亚洲av在线| 久久99热这里只有精品18| 精品国产乱码久久久久久男人| 久久九九热精品免费| 色av中文字幕| 好看av亚洲va欧美ⅴa在| 麻豆久久精品国产亚洲av| 日韩欧美三级三区| 久久精品国产99精品国产亚洲性色| 日本熟妇午夜| 免费高清视频大片| 欧美色欧美亚洲另类二区| 日韩精品中文字幕看吧| 国产成人福利小说| 成年女人毛片免费观看观看9| 国产黄片美女视频| 国产精品日韩av在线免费观看| e午夜精品久久久久久久| 亚洲色图av天堂| 久久久久久久久久黄片| 嫩草影院入口| 欧美3d第一页| 精品国内亚洲2022精品成人| 国产精品1区2区在线观看.| 人人妻,人人澡人人爽秒播| 国产黄a三级三级三级人| 久久久久国内视频| 婷婷丁香在线五月| 99久久99久久久精品蜜桃| 国产午夜福利久久久久久| 国产高清激情床上av| 国产一区二区三区在线臀色熟女| 久久午夜亚洲精品久久| 制服人妻中文乱码| 三级男女做爰猛烈吃奶摸视频| av欧美777| 午夜免费激情av| 亚洲av中文字字幕乱码综合| 国产高清激情床上av| 别揉我奶头~嗯~啊~动态视频| 久久久久免费精品人妻一区二区| 亚洲国产欧洲综合997久久,| www.精华液| 国产高清视频在线播放一区| 午夜福利在线观看吧| 美女被艹到高潮喷水动态| 亚洲国产精品999在线| 久久久久国内视频| www.熟女人妻精品国产| 色综合婷婷激情| 国产亚洲精品久久久com| 国产久久久一区二区三区| 国产亚洲精品久久久久久毛片| 国产蜜桃级精品一区二区三区| 久久天躁狠狠躁夜夜2o2o| 亚洲国产精品久久男人天堂| av天堂中文字幕网| 久久精品aⅴ一区二区三区四区| cao死你这个sao货| 在线播放国产精品三级| 网址你懂的国产日韩在线| 欧美av亚洲av综合av国产av| 淫秽高清视频在线观看| 毛片女人毛片| svipshipincom国产片| 亚洲欧美日韩东京热| 最好的美女福利视频网| 国产精品日韩av在线免费观看| 欧美色视频一区免费| 久久人妻av系列| 黑人巨大精品欧美一区二区mp4| 国产高清激情床上av| 亚洲七黄色美女视频| 免费在线观看亚洲国产| 亚洲 欧美一区二区三区| 国内精品美女久久久久久| 99精品欧美一区二区三区四区| 国内毛片毛片毛片毛片毛片| 欧美黑人欧美精品刺激| 亚洲欧美日韩卡通动漫| 亚洲欧洲精品一区二区精品久久久| 黑人巨大精品欧美一区二区mp4| 深夜精品福利| 成人av在线播放网站| 成人永久免费在线观看视频| 午夜a级毛片| 亚洲五月天丁香| 久久人妻av系列| 老熟妇仑乱视频hdxx| 在线观看美女被高潮喷水网站 | 日韩精品青青久久久久久| 久久中文字幕人妻熟女| 99久久综合精品五月天人人| 国产午夜福利久久久久久| 亚洲精品国产精品久久久不卡| 久久久国产成人免费| 精品一区二区三区av网在线观看| 亚洲成人免费电影在线观看| 一级a爱片免费观看的视频| 欧美zozozo另类| 99久久成人亚洲精品观看| 国产视频一区二区在线看| 国产精品久久久av美女十八| 午夜福利免费观看在线| 国产成人aa在线观看| 脱女人内裤的视频| 亚洲成人久久性| 偷拍熟女少妇极品色| 国产激情久久老熟女| 欧美日韩综合久久久久久 | 天堂√8在线中文| 狂野欧美白嫩少妇大欣赏| 少妇熟女aⅴ在线视频| 久久精品国产清高在天天线| 午夜精品在线福利| 亚洲专区字幕在线| 国产一区二区在线观看日韩 | 99精品欧美一区二区三区四区| 国产美女午夜福利| 男女那种视频在线观看| 成人av一区二区三区在线看| 国产精品电影一区二区三区| 亚洲专区国产一区二区| 香蕉av资源在线| 真人一进一出gif抽搐免费| 天堂av国产一区二区熟女人妻| 午夜免费成人在线视频| 日韩欧美免费精品| 熟女少妇亚洲综合色aaa.| 小说图片视频综合网站| 久久精品影院6| 国产日本99.免费观看| 白带黄色成豆腐渣| 欧美黄色片欧美黄色片| 亚洲人成伊人成综合网2020| 免费看光身美女| 久久99热这里只有精品18| 国产亚洲av嫩草精品影院| 又黄又爽又免费观看的视频| 日本a在线网址| 国产av麻豆久久久久久久| 18禁国产床啪视频网站| 国产成人啪精品午夜网站| 国产av在哪里看| 久久午夜亚洲精品久久| 小说图片视频综合网站| 黄频高清免费视频| 无限看片的www在线观看| 床上黄色一级片| 伦理电影免费视频| 国产一区在线观看成人免费| 精品无人区乱码1区二区| av在线蜜桃| 最近最新免费中文字幕在线| 久久精品国产综合久久久| 99精品欧美一区二区三区四区| 日韩欧美在线二视频| 国产成人系列免费观看| 国产精品亚洲av一区麻豆| ponron亚洲| 给我免费播放毛片高清在线观看| 欧美日韩综合久久久久久 | 91九色精品人成在线观看| 亚洲熟妇中文字幕五十中出| 久久人妻av系列| 男女做爰动态图高潮gif福利片| 最新在线观看一区二区三区| 午夜影院日韩av| 无人区码免费观看不卡| 97碰自拍视频| 黄色片一级片一级黄色片| 国内精品美女久久久久久| 18禁裸乳无遮挡免费网站照片| www日本在线高清视频| 久久中文字幕一级| 国产精品亚洲美女久久久| 午夜福利免费观看在线| 一区二区三区国产精品乱码| 亚洲精品美女久久av网站| 男女那种视频在线观看| 他把我摸到了高潮在线观看| 级片在线观看| 亚洲,欧美精品.| 99久久成人亚洲精品观看| 一区福利在线观看| 日本在线视频免费播放| 美女午夜性视频免费| 欧美丝袜亚洲另类 | 成年免费大片在线观看| 午夜视频精品福利| 亚洲国产色片| 国产高潮美女av| 午夜精品在线福利| 免费在线观看影片大全网站| 香蕉国产在线看| 夜夜看夜夜爽夜夜摸| 免费看十八禁软件| 国产毛片a区久久久久| 观看美女的网站| 91九色精品人成在线观看| 90打野战视频偷拍视频| 99精品在免费线老司机午夜| 亚洲国产欧美人成| 天堂av国产一区二区熟女人妻| 一二三四社区在线视频社区8| 黑人欧美特级aaaaaa片| 免费高清视频大片| 国产av麻豆久久久久久久| 天堂√8在线中文| 天堂网av新在线| 99精品欧美一区二区三区四区| 久久人妻av系列| 免费电影在线观看免费观看| 最近最新中文字幕大全电影3| 亚洲激情在线av| 国产精品永久免费网站| 亚洲片人在线观看| 日韩欧美国产在线观看| 88av欧美| 国内精品久久久久精免费| 亚洲国产欧洲综合997久久,| 国产欧美日韩精品一区二区| 综合色av麻豆| 9191精品国产免费久久| av在线天堂中文字幕| 好看av亚洲va欧美ⅴa在| 久久中文字幕人妻熟女| 欧美日韩国产亚洲二区| 两个人视频免费观看高清| 老司机深夜福利视频在线观看| 俄罗斯特黄特色一大片| 熟妇人妻久久中文字幕3abv| cao死你这个sao货| 麻豆成人av在线观看| 18禁裸乳无遮挡免费网站照片| 国产精品亚洲一级av第二区| 午夜福利视频1000在线观看| 亚洲中文字幕一区二区三区有码在线看 | 国产探花在线观看一区二区| 欧美日韩中文字幕国产精品一区二区三区| 黑人操中国人逼视频| 国产精品爽爽va在线观看网站| 亚洲成av人片在线播放无| 99国产综合亚洲精品| 韩国av一区二区三区四区| 国产精品av久久久久免费| 韩国av一区二区三区四区| 久久中文字幕一级| 色老头精品视频在线观看| 男女视频在线观看网站免费| 久久国产精品影院| 午夜精品在线福利| 变态另类成人亚洲欧美熟女| 国产精品av久久久久免费| 久久久久久久久免费视频了| 日本 av在线| 国产视频内射| 成年女人看的毛片在线观看| a级毛片a级免费在线| 日韩欧美在线二视频| 国产三级中文精品| 久久天躁狠狠躁夜夜2o2o| 欧美日韩国产亚洲二区| www.999成人在线观看| 国产激情久久老熟女| 波多野结衣高清无吗| 天天躁日日操中文字幕| 久久久久久久久中文| 丰满的人妻完整版| 国模一区二区三区四区视频 | 在线永久观看黄色视频| a级毛片在线看网站| 麻豆一二三区av精品| 久久久久久久久久黄片| 精品日产1卡2卡| 好男人电影高清在线观看| 美女cb高潮喷水在线观看 | 搡老妇女老女人老熟妇| 可以在线观看毛片的网站| 黄色女人牲交| 成年女人毛片免费观看观看9| 亚洲一区二区三区色噜噜| 在线观看美女被高潮喷水网站 | 国产精品免费一区二区三区在线| 久久久国产成人精品二区| 老熟妇仑乱视频hdxx| 午夜免费激情av| 老司机在亚洲福利影院| 香蕉丝袜av| 麻豆一二三区av精品| 亚洲精品色激情综合| 日韩中文字幕欧美一区二区| 久久久久久久精品吃奶| 搡老熟女国产l中国老女人| 长腿黑丝高跟| 久久热在线av| 免费av毛片视频| 色av中文字幕| 午夜福利18| 在线国产一区二区在线| 成人一区二区视频在线观看| 亚洲国产精品999在线| 久久久久性生活片| 此物有八面人人有两片| 欧美在线一区亚洲| 一级毛片高清免费大全| 亚洲成人精品中文字幕电影| 给我免费播放毛片高清在线观看| 女同久久另类99精品国产91| 亚洲aⅴ乱码一区二区在线播放| 免费看光身美女| 久久天躁狠狠躁夜夜2o2o| 亚洲精品在线美女| 99久久国产精品久久久| www日本黄色视频网| 午夜成年电影在线免费观看| 婷婷六月久久综合丁香| 国产精品电影一区二区三区| 午夜亚洲福利在线播放| av福利片在线观看| 真人一进一出gif抽搐免费| 神马国产精品三级电影在线观看| 校园春色视频在线观看| 久久久久久久久久黄片| 欧美日本视频| 一个人观看的视频www高清免费观看 | 亚洲精品在线美女| xxxwww97欧美| 丰满的人妻完整版| 黑人欧美特级aaaaaa片| 草草在线视频免费看| 久久精品国产99精品国产亚洲性色| 亚洲av成人一区二区三| 精品久久久久久久人妻蜜臀av| 99精品欧美一区二区三区四区| 国内精品久久久久精免费| 天堂√8在线中文| 国产伦精品一区二区三区四那| 国产亚洲欧美在线一区二区| 免费在线观看日本一区| 久久99热这里只有精品18| 欧美大码av| 国产欧美日韩精品亚洲av| 精品国产乱子伦一区二区三区| 1024香蕉在线观看| 国产精品 欧美亚洲| 欧美在线黄色| 亚洲第一欧美日韩一区二区三区| 精品人妻1区二区| 亚洲人成伊人成综合网2020| 久久欧美精品欧美久久欧美| 午夜精品久久久久久毛片777| 亚洲电影在线观看av| 黄片大片在线免费观看| 变态另类成人亚洲欧美熟女| 亚洲欧美日韩卡通动漫| 又黄又爽又免费观看的视频| 日日干狠狠操夜夜爽| 免费看光身美女| 一夜夜www| 超碰成人久久| 亚洲七黄色美女视频| 日韩国内少妇激情av| а√天堂www在线а√下载| 亚洲最大成人中文| 免费搜索国产男女视频| 亚洲精品一卡2卡三卡4卡5卡| 久久亚洲真实| 成年女人看的毛片在线观看| 成人鲁丝片一二三区免费| 99re在线观看精品视频| 精品一区二区三区av网在线观看| 免费观看人在逋| 国产亚洲精品久久久久久毛片| 亚洲人成伊人成综合网2020| 一进一出好大好爽视频| 久久精品国产综合久久久| 中亚洲国语对白在线视频| ponron亚洲| 熟女少妇亚洲综合色aaa.| 亚洲熟妇熟女久久| 成人高潮视频无遮挡免费网站| 欧美激情在线99| 国产精品99久久99久久久不卡| 床上黄色一级片| 国产精华一区二区三区| 国产极品精品免费视频能看的| a在线观看视频网站| 99久久久亚洲精品蜜臀av| 成人特级av手机在线观看| 亚洲国产欧美网| 国产精品女同一区二区软件 | 亚洲欧美一区二区三区黑人| 人人妻人人看人人澡| 舔av片在线| 曰老女人黄片| 免费人成视频x8x8入口观看| 久久精品91蜜桃| 午夜精品久久久久久毛片777| av中文乱码字幕在线| av片东京热男人的天堂| 精华霜和精华液先用哪个| 亚洲人成伊人成综合网2020| 18禁观看日本| 又爽又黄无遮挡网站| 热99在线观看视频| 九色国产91popny在线| 长腿黑丝高跟| 午夜日韩欧美国产| 午夜亚洲福利在线播放| 欧美日韩精品网址| 中国美女看黄片| 老汉色∧v一级毛片| 无人区码免费观看不卡| 午夜福利18| 亚洲国产中文字幕在线视频| 宅男免费午夜| 99热这里只有是精品50| 亚洲国产欧美人成| 国产伦一二天堂av在线观看| cao死你这个sao货| 黑人巨大精品欧美一区二区mp4| 夜夜爽天天搞| 熟女电影av网| 狂野欧美白嫩少妇大欣赏| 美女大奶头视频| 精品日产1卡2卡| 丰满的人妻完整版| 怎么达到女性高潮| 久久久国产精品麻豆| 免费在线观看视频国产中文字幕亚洲| 一个人观看的视频www高清免费观看 | 国产精品久久久久久久电影 | 婷婷丁香在线五月| 一进一出好大好爽视频| 国产精品电影一区二区三区| 亚洲av成人av| 黑人操中国人逼视频| 欧美在线一区亚洲| 亚洲国产高清在线一区二区三| 国产真实乱freesex| 神马国产精品三级电影在线观看| 国产单亲对白刺激| 一进一出好大好爽视频| 脱女人内裤的视频| 国内毛片毛片毛片毛片毛片| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久视频播放| 国产午夜福利久久久久久| 麻豆国产av国片精品| 午夜影院日韩av| 国产成人啪精品午夜网站| 久久中文看片网| 99精品在免费线老司机午夜| 成在线人永久免费视频| 老熟妇乱子伦视频在线观看| 给我免费播放毛片高清在线观看| 国产三级中文精品| 99久国产av精品| 18禁国产床啪视频网站| 99国产精品99久久久久| 黄频高清免费视频| 在线观看免费视频日本深夜| 啦啦啦观看免费观看视频高清| 99国产综合亚洲精品| 婷婷精品国产亚洲av| 麻豆国产av国片精品| 国产野战对白在线观看| 欧美国产日韩亚洲一区| 99在线人妻在线中文字幕| 一本综合久久免费| 久久天躁狠狠躁夜夜2o2o| 男的添女的下面高潮视频| 极品教师在线视频| АⅤ资源中文在线天堂|