• <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)研究
    一二三四中文在线观看免费高清| 最近中文字幕高清免费大全6| 99精国产麻豆久久婷婷| 欧美精品高潮呻吟av久久| 黄色毛片三级朝国网站| 国产精品三级大全| 亚洲欧洲日产国产| 9热在线视频观看99| av黄色大香蕉| 在线天堂中文资源库| 蜜臀久久99精品久久宅男| 亚洲av.av天堂| 国产免费视频播放在线视频| 韩国av在线不卡| 成年动漫av网址| 日本免费在线观看一区| 国产一区亚洲一区在线观看| 99久国产av精品国产电影| 97在线人人人人妻| 在线天堂最新版资源| 免费在线观看黄色视频的| 国产日韩一区二区三区精品不卡| 欧美3d第一页| 在线天堂中文资源库| 国产 一区精品| 亚洲国产精品一区三区| 亚洲av综合色区一区| 久久久欧美国产精品| 男的添女的下面高潮视频| 黑人高潮一二区| 日韩欧美精品免费久久| 在线观看一区二区三区激情| 天天操日日干夜夜撸| 各种免费的搞黄视频| 99久久综合免费| 亚洲国产成人一精品久久久| 久久女婷五月综合色啪小说| 亚洲国产av新网站| 亚洲av免费高清在线观看| 丁香六月天网| 亚洲av电影在线进入| 看非洲黑人一级黄片| 观看美女的网站| 国产69精品久久久久777片| 久久久欧美国产精品| 亚洲一区二区三区欧美精品| 天堂8中文在线网| 免费看av在线观看网站| 国产一区二区激情短视频 | 又大又黄又爽视频免费| 91精品伊人久久大香线蕉| 中文字幕制服av| 免费看光身美女| 国产熟女午夜一区二区三区| 大片免费播放器 马上看| 天天躁夜夜躁狠狠躁躁| 性高湖久久久久久久久免费观看| 桃花免费在线播放| 亚洲精品日本国产第一区| av在线老鸭窝| 一级片免费观看大全| 男女边吃奶边做爰视频| 亚洲欧美色中文字幕在线| av黄色大香蕉| 亚洲激情五月婷婷啪啪| 亚洲色图综合在线观看| 女人被躁到高潮嗷嗷叫费观| 黄片无遮挡物在线观看| 精品国产一区二区久久| kizo精华| 黄色一级大片看看| 大话2 男鬼变身卡| 日本午夜av视频| 亚洲av福利一区| av国产精品久久久久影院| 美女国产高潮福利片在线看| 亚洲国产成人一精品久久久| 久久精品aⅴ一区二区三区四区 | 啦啦啦啦在线视频资源| 国产精品麻豆人妻色哟哟久久| 亚洲成av片中文字幕在线观看 | 免费看不卡的av| av.在线天堂| 极品人妻少妇av视频| 国产成人91sexporn| 午夜福利视频精品| tube8黄色片| 国产成人精品久久久久久| 777米奇影视久久| 免费日韩欧美在线观看| 免费少妇av软件| 美女大奶头黄色视频| 欧美少妇被猛烈插入视频| 99热6这里只有精品| 亚洲精华国产精华液的使用体验| av免费观看日本| 精品少妇内射三级| 欧美精品高潮呻吟av久久| 天天操日日干夜夜撸| 激情五月婷婷亚洲| 99热这里只有是精品在线观看| 精品少妇黑人巨大在线播放| 精品国产一区二区久久| 亚洲成人手机| 亚洲欧美色中文字幕在线| 成人毛片a级毛片在线播放| 精品久久国产蜜桃| 欧美3d第一页| 亚洲人成77777在线视频| 美女国产视频在线观看| 色婷婷久久久亚洲欧美| 深夜精品福利| 亚洲欧洲国产日韩| 99九九在线精品视频| 91成人精品电影| 亚洲欧美精品自产自拍| 成人漫画全彩无遮挡| 国产精品蜜桃在线观看| 大陆偷拍与自拍| 成人国语在线视频| 在线 av 中文字幕| 内地一区二区视频在线| 久久久精品94久久精品| 在线精品无人区一区二区三| 精品99又大又爽又粗少妇毛片| 99久久中文字幕三级久久日本| 青春草亚洲视频在线观看| 精品久久久精品久久久| 妹子高潮喷水视频| 日韩精品免费视频一区二区三区 | 国产av一区二区精品久久| 另类亚洲欧美激情| 久久精品国产亚洲av涩爱| 人成视频在线观看免费观看| 又大又黄又爽视频免费| 人人妻人人添人人爽欧美一区卜| 黄色视频在线播放观看不卡| 极品少妇高潮喷水抽搐| www.av在线官网国产| 亚洲熟女精品中文字幕| 9色porny在线观看| 高清欧美精品videossex| 午夜福利视频在线观看免费| 欧美亚洲 丝袜 人妻 在线| 久久久久久久精品精品| 男女午夜视频在线观看 | 日韩一区二区三区影片| 国产av码专区亚洲av| 一级a做视频免费观看| 国产无遮挡羞羞视频在线观看| 成人免费观看视频高清| 又粗又硬又长又爽又黄的视频| 亚洲精品成人av观看孕妇| 九草在线视频观看| 一边亲一边摸免费视频| 国产精品蜜桃在线观看| 欧美精品一区二区大全| av卡一久久| 国产日韩欧美在线精品| 国产 一区精品| 美国免费a级毛片| 亚洲丝袜综合中文字幕| 精品国产乱码久久久久久小说| 国产精品久久久久成人av| av线在线观看网站| 亚洲三级黄色毛片| 在线观看三级黄色| 91国产中文字幕| 婷婷色综合大香蕉| 国产成人91sexporn| 国产男女内射视频| av国产久精品久网站免费入址| 亚洲国产成人一精品久久久| 成人午夜精彩视频在线观看| 成年动漫av网址| 精品国产露脸久久av麻豆| 色网站视频免费| 成年人午夜在线观看视频| 韩国av在线不卡| 高清毛片免费看| 黑人猛操日本美女一级片| 国产高清国产精品国产三级| 99久国产av精品国产电影| 久久国产亚洲av麻豆专区| 国产一级毛片在线| 午夜福利,免费看| 久久精品久久久久久噜噜老黄| 少妇猛男粗大的猛烈进出视频| av在线app专区| 色网站视频免费| 大陆偷拍与自拍| 精品福利永久在线观看| 中文字幕人妻熟女乱码| 日日摸夜夜添夜夜爱| 久久韩国三级中文字幕| 午夜免费鲁丝| 大码成人一级视频| 精品少妇久久久久久888优播| 韩国高清视频一区二区三区| 热99久久久久精品小说推荐| 边亲边吃奶的免费视频| 又粗又硬又长又爽又黄的视频| 老司机影院成人| 亚洲欧美色中文字幕在线| 中文字幕另类日韩欧美亚洲嫩草| 久久这里只有精品19| av卡一久久| 99久久人妻综合| 99精国产麻豆久久婷婷| 女人精品久久久久毛片| 国产免费一区二区三区四区乱码| 亚洲图色成人| 老司机影院毛片| 少妇被粗大猛烈的视频| 精品少妇久久久久久888优播| 九九爱精品视频在线观看| 国产高清不卡午夜福利| 成人手机av| 18禁在线无遮挡免费观看视频| 秋霞伦理黄片| 国产国语露脸激情在线看| www日本在线高清视频| 日韩,欧美,国产一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 国产精品国产三级国产专区5o| 99国产综合亚洲精品| 欧美精品av麻豆av| 伊人久久国产一区二区| 日本与韩国留学比较| 欧美人与性动交α欧美软件 | 午夜激情久久久久久久| 亚洲国产欧美日韩在线播放| 日韩一本色道免费dvd| 毛片一级片免费看久久久久| 免费在线观看完整版高清| 色吧在线观看| 亚洲成人一二三区av| 大话2 男鬼变身卡| 午夜老司机福利剧场| 熟女av电影| 成人国语在线视频| 久久精品国产亚洲av天美| 日韩av在线免费看完整版不卡| 国产精品久久久久久av不卡| 亚洲精品美女久久久久99蜜臀 | 乱人伦中国视频| 久久久久人妻精品一区果冻| 亚洲综合精品二区| 亚洲图色成人| 欧美激情 高清一区二区三区| 中文天堂在线官网| 日日撸夜夜添| 天堂8中文在线网| 丰满少妇做爰视频| 黄色配什么色好看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲伊人久久精品综合| 性色avwww在线观看| 亚洲少妇的诱惑av| 成人免费观看视频高清| 亚洲一级一片aⅴ在线观看| av不卡在线播放| 日韩制服丝袜自拍偷拍| 少妇精品久久久久久久| 国产探花极品一区二区| 一级,二级,三级黄色视频| 91久久精品国产一区二区三区| 精品人妻一区二区三区麻豆| 看非洲黑人一级黄片| 亚洲五月色婷婷综合| 肉色欧美久久久久久久蜜桃| 午夜91福利影院| 久久av网站| 欧美精品国产亚洲| 91成人精品电影| 中文字幕人妻熟女乱码| 欧美亚洲日本最大视频资源| 亚洲精品久久久久久婷婷小说| 亚洲综合精品二区| 色网站视频免费| 久久青草综合色| 国产精品人妻久久久久久| 亚洲精品色激情综合| 中文字幕亚洲精品专区| 大香蕉久久成人网| 少妇人妻精品综合一区二区| 日本-黄色视频高清免费观看| 少妇猛男粗大的猛烈进出视频| 宅男免费午夜| 狂野欧美激情性xxxx在线观看| 国产又爽黄色视频| 黑人高潮一二区| 亚洲图色成人| 国产精品秋霞免费鲁丝片| a级毛色黄片| 国产福利在线免费观看视频| 亚洲精品国产av成人精品| 捣出白浆h1v1| 久久综合国产亚洲精品| a级片在线免费高清观看视频| 在线观看三级黄色| 伊人久久国产一区二区| 国产一区二区在线观看av| 欧美3d第一页| 亚洲美女黄色视频免费看| 亚洲欧美中文字幕日韩二区| 纵有疾风起免费观看全集完整版| 亚洲精品日本国产第一区| 插逼视频在线观看| 久久 成人 亚洲| 中国美白少妇内射xxxbb| 免费高清在线观看视频在线观看| 9色porny在线观看| 9热在线视频观看99| 久久久精品94久久精品| 狂野欧美激情性bbbbbb| 久久久久久久久久成人| 一区二区三区精品91| 成人影院久久| 在线亚洲精品国产二区图片欧美| 久久国产精品大桥未久av| 夫妻午夜视频| 久热久热在线精品观看| 精品国产露脸久久av麻豆| www.熟女人妻精品国产 | 亚洲欧美一区二区三区黑人 | 少妇猛男粗大的猛烈进出视频| 国产又爽黄色视频| 国产精品无大码| 久久精品国产综合久久久 | 国产片内射在线| 18禁观看日本| 亚洲av欧美aⅴ国产| 亚洲 欧美一区二区三区| 丝袜喷水一区| 我的女老师完整版在线观看| 午夜福利视频在线观看免费| 国产一区亚洲一区在线观看| 欧美97在线视频| 女人精品久久久久毛片| 丝袜喷水一区| 欧美精品国产亚洲| 成年女人在线观看亚洲视频| 亚洲内射少妇av| 精品一区二区三卡| 一级a做视频免费观看| 青春草视频在线免费观看| 大码成人一级视频| 成年女人在线观看亚洲视频| 一级片'在线观看视频| 免费少妇av软件| 亚洲成国产人片在线观看| 国产日韩欧美亚洲二区| 各种免费的搞黄视频| 老女人水多毛片| 亚洲伊人久久精品综合| 日本91视频免费播放| 七月丁香在线播放| 天天影视国产精品| 有码 亚洲区| 日韩视频在线欧美| 中文精品一卡2卡3卡4更新| 999精品在线视频| 国产精品麻豆人妻色哟哟久久| 久久久久久久国产电影| 久久婷婷青草| 精品国产露脸久久av麻豆| 国产亚洲午夜精品一区二区久久| 美女视频免费永久观看网站| 你懂的网址亚洲精品在线观看| 国产日韩一区二区三区精品不卡| 丁香六月天网| 街头女战士在线观看网站| 性色avwww在线观看| 人妻一区二区av| 色婷婷av一区二区三区视频| 99热6这里只有精品| 香蕉丝袜av| 亚洲情色 制服丝袜| 国产又爽黄色视频| 精品国产乱码久久久久久小说| 欧美日韩av久久| 色网站视频免费| 综合色丁香网| 免费少妇av软件| 欧美精品国产亚洲| 国产精品国产av在线观看| 多毛熟女@视频| 日韩欧美一区视频在线观看| av一本久久久久| 18禁国产床啪视频网站| 看免费av毛片| 久久久久久久久久久免费av| 欧美日韩国产mv在线观看视频| 美女脱内裤让男人舔精品视频| 一本久久精品| 日本黄大片高清| 五月天丁香电影| 国产成人精品无人区| 熟女人妻精品中文字幕| 亚洲精品久久久久久婷婷小说| 国产精品一二三区在线看| 亚洲美女黄色视频免费看| 欧美日韩亚洲高清精品| 久久久久久伊人网av| 欧美精品高潮呻吟av久久| 久久人人爽人人爽人人片va| 国产精品人妻久久久久久| 黄片播放在线免费| 性高湖久久久久久久久免费观看| 亚洲欧美清纯卡通| 亚洲一码二码三码区别大吗| 最近2019中文字幕mv第一页| 国产精品国产三级专区第一集| 免费看av在线观看网站| 日韩不卡一区二区三区视频在线| tube8黄色片| 2018国产大陆天天弄谢| 欧美激情 高清一区二区三区| 人人妻人人添人人爽欧美一区卜| 国产成人精品在线电影| 日本欧美国产在线视频| 亚洲av中文av极速乱| 另类精品久久| 成人国产麻豆网| 丁香六月天网| 欧美国产精品一级二级三级| 在线精品无人区一区二区三| 18禁在线无遮挡免费观看视频| 狠狠婷婷综合久久久久久88av| 美女国产高潮福利片在线看| 欧美日韩一区二区视频在线观看视频在线| 高清不卡的av网站| 一区在线观看完整版| 亚洲精品色激情综合| 久热久热在线精品观看| 大陆偷拍与自拍| 国产成人a∨麻豆精品| 亚洲在久久综合| 九色成人免费人妻av| 国产成人欧美| www日本在线高清视频| 精品亚洲成a人片在线观看| 亚洲欧美精品自产自拍| 亚洲国产欧美在线一区| 狂野欧美激情性bbbbbb| 欧美人与善性xxx| 欧美97在线视频| 26uuu在线亚洲综合色| 国产日韩欧美视频二区| 边亲边吃奶的免费视频| 国产精品久久久久成人av| 波野结衣二区三区在线| 中文字幕av电影在线播放| 欧美人与善性xxx| 亚洲欧洲精品一区二区精品久久久 | 国产成人免费无遮挡视频| 日韩在线高清观看一区二区三区| 国产男女内射视频| 日本欧美视频一区| 一本—道久久a久久精品蜜桃钙片| 亚洲国产成人一精品久久久| 午夜av观看不卡| 亚洲五月色婷婷综合| 国产永久视频网站| 亚洲四区av| 亚洲欧美一区二区三区黑人 | 午夜福利,免费看| 国产成人a∨麻豆精品| 国产免费福利视频在线观看| 国产亚洲一区二区精品| 亚洲伊人久久精品综合| 人妻系列 视频| 青春草国产在线视频| 中文字幕人妻丝袜制服| 好男人视频免费观看在线| 日韩熟女老妇一区二区性免费视频| av国产精品久久久久影院| 亚洲精品美女久久久久99蜜臀 | 大陆偷拍与自拍| 中文字幕av电影在线播放| 久久久久国产网址| 黄色毛片三级朝国网站| www.av在线官网国产| 久久韩国三级中文字幕| 老司机影院毛片| 夜夜骑夜夜射夜夜干| 中文字幕另类日韩欧美亚洲嫩草| 99九九在线精品视频| 亚洲中文av在线| 国产免费视频播放在线视频| 成年人免费黄色播放视频| 日本av手机在线免费观看| 国产成人一区二区在线| 国产亚洲欧美精品永久| 天堂俺去俺来也www色官网| av.在线天堂| 91精品国产国语对白视频| 日韩熟女老妇一区二区性免费视频| 最近中文字幕高清免费大全6| 免费黄网站久久成人精品| 美女大奶头黄色视频| 我要看黄色一级片免费的| 欧美97在线视频| 日韩伦理黄色片| 99热全是精品| 亚洲av日韩在线播放| 卡戴珊不雅视频在线播放| 精品人妻熟女毛片av久久网站| 亚洲精品国产色婷婷电影| 18禁观看日本| 精品国产一区二区久久| 久久久久久久大尺度免费视频| 国产精品一国产av| a级毛色黄片| 男的添女的下面高潮视频| 一二三四在线观看免费中文在 | 亚洲精品美女久久av网站| 亚洲欧美成人综合另类久久久| 女人被躁到高潮嗷嗷叫费观| 91国产中文字幕| 国产爽快片一区二区三区| 亚洲综合色惰| 久久婷婷青草| 黄片播放在线免费| 色94色欧美一区二区| 亚洲欧美清纯卡通| 欧美精品人与动牲交sv欧美| 国产日韩欧美视频二区| 欧美bdsm另类| tube8黄色片| 激情五月婷婷亚洲| 久久鲁丝午夜福利片| 一级毛片我不卡| 伊人亚洲综合成人网| 欧美3d第一页| 极品人妻少妇av视频| 青春草视频在线免费观看| 亚洲av电影在线观看一区二区三区| 国产成人精品福利久久| 国语对白做爰xxxⅹ性视频网站| 亚洲av欧美aⅴ国产| 妹子高潮喷水视频| 精品一区二区免费观看| 亚洲国产看品久久| 亚洲情色 制服丝袜| 波野结衣二区三区在线| a级毛片黄视频| 欧美精品一区二区免费开放| 久久精品国产a三级三级三级| 日韩中文字幕视频在线看片| 考比视频在线观看| 亚洲色图综合在线观看| 你懂的网址亚洲精品在线观看| 欧美xxxx性猛交bbbb| 欧美人与性动交α欧美精品济南到 | 久久久精品94久久精品| 成年人免费黄色播放视频| 国产免费一区二区三区四区乱码| 久久久久久久久久久免费av| 在线观看美女被高潮喷水网站| 国产精品不卡视频一区二区| 久久热在线av| 不卡视频在线观看欧美| 毛片一级片免费看久久久久| 美女脱内裤让男人舔精品视频| 2021少妇久久久久久久久久久| 国产精品久久久久久久电影| 我要看黄色一级片免费的| 久久免费观看电影| 国产精品国产三级国产专区5o| 黄片无遮挡物在线观看| 国产爽快片一区二区三区| 亚洲国产欧美在线一区| 尾随美女入室| 午夜福利视频精品| 欧美xxⅹ黑人| 色视频在线一区二区三区| 亚洲人成77777在线视频| 一本大道久久a久久精品| 男女啪啪激烈高潮av片| 999精品在线视频| 色吧在线观看| 国产不卡av网站在线观看| 看免费av毛片| 国产一区二区三区综合在线观看 | 赤兔流量卡办理| 精品一区在线观看国产| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 五月伊人婷婷丁香| 国产亚洲一区二区精品| 日韩精品有码人妻一区| 极品少妇高潮喷水抽搐| 国产精品国产三级国产av玫瑰| 日韩电影二区| 欧美精品av麻豆av| 中国美白少妇内射xxxbb| 一二三四中文在线观看免费高清| 亚洲av电影在线进入| 色婷婷久久久亚洲欧美| 国产精品国产三级专区第一集| 在线亚洲精品国产二区图片欧美| 国产成人免费观看mmmm| 久久毛片免费看一区二区三区| 春色校园在线视频观看| 国产一区二区在线观看日韩| 热99久久久久精品小说推荐| 国产精品久久久久成人av| 国产极品粉嫩免费观看在线| 久久亚洲国产成人精品v| 国产片特级美女逼逼视频| 日本vs欧美在线观看视频| 国产探花极品一区二区| 久久久久国产精品人妻一区二区| 少妇人妻久久综合中文|