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

    基于非穩(wěn)態(tài)模擬的SAE 車模氣動(dòng)減阻降噪研究

    2023-03-08 06:33:18楊小龍龔繁龔政黃元康
    關(guān)鍵詞:尾流聲壓級(jí)尾部

    楊小龍,龔繁 ,龔政,黃元康

    (湖南大學(xué) 機(jī)械與運(yùn)載工程學(xué)院,湖南 長沙 410082)

    進(jìn)入21 世紀(jì)以來,汽車行駛速度隨著汽車行業(yè)、高速公路的快速發(fā)展而不斷提高.汽車高速行駛時(shí)的氣動(dòng)阻力對(duì)汽車的燃油經(jīng)濟(jì)性產(chǎn)生重要的影響;噪聲問題也隨著行駛速度的增加而嚴(yán)重影響乘客的乘坐舒適性.降低汽車空氣阻力能有效地節(jié)省燃油,而降低氣動(dòng)噪聲能有效地提升乘客的舒適度.減少汽車能源消耗與降低汽車氣動(dòng)噪聲成為當(dāng)代車輛設(shè)計(jì)中考慮的兩大問題,并且成為近些年汽車行業(yè)研究的熱點(diǎn),國內(nèi)外都對(duì)此開展了一定的研究.

    在氣動(dòng)減阻研究方面,Choi 等人[1-2]對(duì)簡單的汽車模型應(yīng)用不同的減阻裝置,研究減阻裝置對(duì)流場(chǎng)結(jié)構(gòu)的控制效果.Pujals 等人[3]研究了在模型頂面放置圓形渦發(fā)生器時(shí)的減阻能力.Sharma 等人[4]利用CFD 方法研究了渦發(fā)生器(VGS)對(duì)MIRA 模型尾流結(jié)構(gòu)的影響.Eulalie 等人[5]以方背Ahmed 為研究對(duì)象,采用了被動(dòng)和主動(dòng)的控制策略,最大可獲得10%的減阻效果.Bruneau 等人[6]以方背 Ahmed 模型為研究對(duì)象,通過數(shù)值模擬方法研究了被動(dòng)和主動(dòng)控制方法耦合時(shí)的減阻效果,結(jié)果表明結(jié)合被動(dòng)和主動(dòng)控制方法實(shí)現(xiàn)了高達(dá)30%的減阻率.Wieser等人[7]以階梯背轎車為研究對(duì)象,通過風(fēng)洞試驗(yàn)方法研究了渦發(fā)生器安裝位置對(duì)汽車氣動(dòng)特性的影響.黃莎等人[8]采用基于Realizablek-ε的數(shù)值計(jì)算方法探索在列車尾部不同位置、不同速度射流對(duì)列車周圍湍流流場(chǎng)和氣動(dòng)阻力的影響規(guī)律,在不同的車速下找到最佳射流位置,以提高整車的減阻率.張迪分別以Ahmed、GTS 和 MIRA 模型為對(duì)象,研究了非光滑表面對(duì)三種模型的氣動(dòng)減阻效果,并且對(duì)Ahmed 和GTS 模型設(shè)計(jì)了優(yōu)化方案,優(yōu)化后最大的減阻效果分別為4.69%和5.18%[9].楊小龍等人[10]研究了非光滑表面對(duì)汽車氣動(dòng)性能的影響,并簡要闡述了非光滑表面的減阻機(jī)理.許建民等人[11]采用計(jì)算流體動(dòng)力學(xué)方法研究了側(cè)風(fēng)環(huán)境下減阻裝置對(duì)重型貨車氣動(dòng)特性的影響,發(fā)現(xiàn)復(fù)合減阻裝置在所有側(cè)風(fēng)工況下均具有顯著的減阻效果.

    對(duì)于氣動(dòng)噪聲研究,Oettle等人[12]綜合介紹了汽車氣動(dòng)噪聲的產(chǎn)生部位以及不同部位氣動(dòng)噪聲產(chǎn)生的原理.Purohit 等人[13]通過計(jì)算氣動(dòng)聲學(xué)混合方法研究了外激勵(lì)對(duì)柔性結(jié)構(gòu)遠(yuǎn)場(chǎng)氣動(dòng)噪聲的影響,結(jié)果表明有外部激勵(lì)時(shí)產(chǎn)生的遠(yuǎn)場(chǎng)噪聲相較于無外部激勵(lì)時(shí)降低了近40%.Manikandan 等人[14]以一款商用救護(hù)車模型為研究對(duì)象,分析了被動(dòng)控制方法對(duì)救護(hù)車模型氣動(dòng)噪聲的影響.汪碩[15]研究了仿生非光滑表面對(duì)汽車外流場(chǎng)特性和氣動(dòng)特性的影響,分析了凹坑分布規(guī)律、尺寸和在汽車表面的分布位置等對(duì)阻力系數(shù)和氣動(dòng)噪聲的影響.陳鑫等人[16]使用 DrivAer 汽車模型來研究仿生非光滑車外后視鏡罩減阻降噪機(jī)理,結(jié)果表明車外后視鏡上應(yīng)用仿生非光滑結(jié)構(gòu),使整車阻力降低5.9%,側(cè)窗外響度降低 19.4%.林肖輝[17]以某品牌量產(chǎn)車模型為研究對(duì)象,分析了周圍流場(chǎng)的氣動(dòng)特性和車身表面的噪聲源及遠(yuǎn)場(chǎng)監(jiān)測(cè)點(diǎn)處的聲壓級(jí)頻譜圖,并通過對(duì)車輪輪轂和后視鏡的造型優(yōu)化使得整車噪聲得到有效的控制.

    上述文獻(xiàn)表明,相關(guān)學(xué)者分別對(duì)氣動(dòng)減阻以及氣動(dòng)噪聲做了一定的研究.但是,針對(duì)同時(shí)降低氣動(dòng)阻力和氣動(dòng)噪聲的研究較少.被動(dòng)減阻裝置依靠改變流場(chǎng)結(jié)構(gòu)來達(dá)到降低氣動(dòng)阻力的目的,但有可能引入新的擾動(dòng),從而增大噪聲.如何在降低氣動(dòng)阻力的同時(shí)有效降低噪聲是一個(gè)比較新穎的研究方向.此外,氣動(dòng)噪聲研究需要知道流場(chǎng)中非定常脈動(dòng)壓力信息,并不能適用于傳統(tǒng)的穩(wěn)態(tài)模擬,需要采用非穩(wěn)態(tài)模型,這對(duì)模擬也提出了一定的挑戰(zhàn).另外,注意到實(shí)際車型幾何構(gòu)造復(fù)雜,底盤部件眾多,導(dǎo)致計(jì)算量過大,同時(shí)不同車型之間缺乏對(duì)比和參考,因此目前學(xué)術(shù)研究仍主要采用標(biāo)準(zhǔn)化車模.同時(shí),模型本身及其研究數(shù)據(jù),可在空氣動(dòng)力學(xué)數(shù)據(jù)定義、實(shí)驗(yàn)標(biāo)準(zhǔn)、計(jì)算標(biāo)準(zhǔn)等方面成為溝通語言,便于數(shù)據(jù)共享,可大量減少汽車開發(fā)過程中的重復(fù)性工作[18].

    本文以20°后背傾角階梯背式SAE 模型為參考模型,在模型后傾斜面頂部和尾部添加新型渦發(fā)生器(VGS)和溝槽(RTS),基于非穩(wěn)態(tài)模擬的同時(shí)研究了被動(dòng)減阻裝置在汽車氣動(dòng)減阻和氣動(dòng)降噪方面的積極影響.本文的研究內(nèi)容對(duì)于研究汽車協(xié)同減阻降噪的方法具有一定的借鑒意義.

    1 仿真模型的建立

    1.1 SAE模型與計(jì)算域

    為了不同研究者對(duì)比研究的方便,有多種標(biāo)準(zhǔn)化的車模被提出,包括Ahmed[19]、MIRA[20]、SAE[21]、GTS[22]等.本文選擇20°后背傾角階梯背式SAE 模型作為研究對(duì)象,因?yàn)槠渚哂械湫偷钠囂攸c(diǎn),同時(shí)有相應(yīng)的實(shí)驗(yàn)數(shù)據(jù)作為參考[23].圖1 為SAE 參考模型,1∶ 5 尺寸縮放的模型在諸如風(fēng)洞阻塞比、可制造性、可實(shí)現(xiàn)雷諾數(shù)范圍方面提供了折中的效果.SAE 模型是MIRA 模型的改進(jìn)模型,模型頭部是一個(gè)30°傾斜表面且有著較大的前緣半徑,底板靠近尾部位置有一個(gè)用于抵消前傾斜面產(chǎn)生的大俯仰力矩的6°擴(kuò)散器結(jié)構(gòu),模型的尺寸為L(長) ×W(寬) ×H(高)=840 mm × 320 mm × 240 mm.如圖2 所示,計(jì)算域的尺寸為Ω=10L× 7W× 5H,阻塞比等于2.86%.SAE模型頭部倒角面離入口表面的距離為2L,尾部垂直面距離出口表面7L.

    圖1 SAE參考模型Fig.1 SAE reference model

    圖2 SAE參考模型計(jì)算域尺寸Fig.2 Calculated domain dimensions for the SAE reference model

    1.2 計(jì)算模型

    汽車流場(chǎng)具有三維非定常大分離的特點(diǎn),而小尺度的減阻裝置增加了流動(dòng)的復(fù)雜性,為了捕捉非定常流場(chǎng)的渦結(jié)構(gòu),考慮到傳統(tǒng)的雷諾平均模擬(RANS)很難滿足計(jì)算精度要求,本文采用了基于SST 模型的改進(jìn)延遲分離渦模擬(IDDES)方法.IDDES[24-25]湍流模型既能捕捉非定常流場(chǎng)結(jié)構(gòu),又相比DNS、LES 具有計(jì)算量小的優(yōu)點(diǎn),目前在工程模擬方面應(yīng)用比較多,但在汽車模擬方面應(yīng)用還比較少.而近場(chǎng)噪聲源使用計(jì)算氣動(dòng)聲學(xué)方法(Computa?tional Aeroacoustics,CAA)對(duì)氣動(dòng)噪聲直接進(jìn)行計(jì)算,CAA 方法把聲音的傳播也當(dāng)成流體現(xiàn)象,因此在求解 N-S 方程時(shí)能夠同時(shí)計(jì)算聲音的產(chǎn)生和傳播,在計(jì)算流場(chǎng)的同時(shí)對(duì)聲壓脈動(dòng)進(jìn)行求解.

    遠(yuǎn)場(chǎng)噪聲則采用了FW-H 聲類比方法.Light hill聲類比理論控制方程是由N-S 方程推導(dǎo)出的波動(dòng)方程,可表示為:

    FW-H方程如下:

    式中:c為聲速;Tij是Light hill 應(yīng)力張量;H(f)為Heaviside 函數(shù);δ(f)為Dirac delta 函數(shù);Pij為壓應(yīng)力張量;p′和ρ分別為聲壓和密度脈動(dòng);ρ0為平均密度;u為流體流速;v為平面速度;n為指向外的單位法向向量.方程中聲源面被定義為f=0,x表示監(jiān)測(cè)點(diǎn)的位置.

    聲壓級(jí)(SPL)和總聲壓級(jí)(OASPL)定義如下:

    式中:prms為實(shí)測(cè)聲壓;pmean為平均聲壓;pref=2×10-5Pa,為參考聲壓.

    1.3 網(wǎng)格劃分和邊界條件

    本文采用多面體網(wǎng)格進(jìn)行數(shù)值仿真,圖3 展示的是計(jì)算域流向中心截面的網(wǎng)格分布,并對(duì)模型周圍進(jìn)行了網(wǎng)格加密.

    圖3 SAE原模型截面多面體網(wǎng)格Fig.3 Polyhedral mesh of the section of the SAE original model

    流場(chǎng)計(jì)算域入口為U0=40 m/s 的均勻流速,基于SAE 模型長度的雷諾數(shù)為ReL=2.32 × 106.出口邊界為壓力等于大氣壓的壓力出口,SAE 模型和計(jì)算域底面設(shè)置為無滑移壁面,計(jì)算域的兩個(gè)側(cè)面和頂面都設(shè)置為對(duì)稱面.時(shí)間步長為?t=0.000 1 s,有效頻率范圍計(jì)算時(shí)間為1.1 s,時(shí)均時(shí)間約為52個(gè)特征周期.

    1.4 網(wǎng)格無關(guān)性分析與模型驗(yàn)證

    本研究一共采用粗糙網(wǎng)格、中等網(wǎng)格和精細(xì)網(wǎng)格3 種網(wǎng)格分別進(jìn)行了計(jì)算.圖4 展示了尾部斜面、行李艙蓋和尾部垂直面上的時(shí)均壓力系數(shù),并與Daniel Wood 等人[23]的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了比較.圖4(a)、4(b)、4(c)分別對(duì)應(yīng)粗糙網(wǎng)格、中等網(wǎng)格和精細(xì)網(wǎng)格的結(jié)果,圖4(d)為Daniel Wood 等人的實(shí)驗(yàn)測(cè)量結(jié)果.從圖中可知,粗糙網(wǎng)格結(jié)果與中等網(wǎng)格和精細(xì)網(wǎng)格結(jié)果存在很直觀的誤差,而中等網(wǎng)格和精細(xì)網(wǎng)格結(jié)果之間誤差較小,同時(shí)這兩種結(jié)果較好地滿足實(shí)驗(yàn)測(cè)量結(jié)果.為了進(jìn)一步對(duì)尾部斜面的流場(chǎng)進(jìn)行定量分析,圖5給出了中心截面上沿x方向的尾部斜面和行李艙蓋上方的無量綱時(shí)均速度分布,結(jié)果表明粗糙網(wǎng)格沒有準(zhǔn)確預(yù)測(cè)尾部斜面上的時(shí)均速度,而中等網(wǎng)格和精細(xì)網(wǎng)格能較好地與實(shí)驗(yàn)結(jié)果相吻合,同時(shí)這也驗(yàn)證了計(jì)算模型的準(zhǔn)確性.

    圖4 仿真模擬結(jié)果與實(shí)驗(yàn)結(jié)果的時(shí)均壓力系數(shù)對(duì)比Fig.4 Comparison of time-mean pressure coefficients between simulation results and experimental results

    圖5 中心截面上沿x方向的無量綱時(shí)均速度分布Fig.5 Dimensionless mean velocity distribution along the X direction on the central section

    除此之外,阻力系數(shù)的計(jì)算結(jié)果如表1 所示.粗糙網(wǎng)格相對(duì)于中等網(wǎng)格的誤差為3.19%,而精細(xì)網(wǎng)格相對(duì)于中等網(wǎng)格的誤差為0.49%.為了節(jié)約計(jì)算資源,本文后續(xù)的研究使用具有網(wǎng)格數(shù)量較少,但能保證計(jì)算精度的中等網(wǎng)格.

    表1 SAE模型3種網(wǎng)格分辨率的網(wǎng)格數(shù)量及阻力系數(shù)Tab.1 The number of meshes and the drag coefficient of the three mesh resolutions of the SAE model

    2 VGS減阻效果研究

    2.1 VGS結(jié)構(gòu)尺寸

    目前被動(dòng)減阻裝置有渦發(fā)生器、仿生非光滑表面、溝槽和多孔介質(zhì)等,這里設(shè)計(jì)了一種新型的半圓柱型渦發(fā)生器和溝槽.半圓柱型渦發(fā)生器相比圓柱或其他渦發(fā)生器,具有尺寸小、結(jié)構(gòu)簡單、易于布置等特點(diǎn).本文將渦發(fā)生器應(yīng)用在SAE 模型尾部斜面的上沿線附近,結(jié)構(gòu)如圖6 所示.22 個(gè)渦發(fā)生器均勻分布在尾部斜面的橫向方向上,直徑D=5 mm,其中相鄰渦發(fā)生器間隔L=12.5 mm.

    圖6 附加VGS的SAE模型的尺寸參數(shù)Fig.6 Dimensional parameters of the SAE model with VGS

    2.2 尾流結(jié)構(gòu)對(duì)比分析

    為了研究VGS 的減阻效果,首先對(duì)比分析了尾流流線.圖7 為SAE 模型和附加VGS 的模型在z=0截面上的流向時(shí)均速度流線.附加VGS 的模型尾部斜面上的分離渦渦1 消失,因?yàn)閂GS 抑制了流體的分離,使流體經(jīng)過斜面時(shí)始終附著在斜面上,阻礙了分離渦的形成.通過比較尾部反向旋轉(zhuǎn)渦的渦核中心位置,可知VGS對(duì)尾渦的影響相對(duì)較小,僅表現(xiàn)出對(duì)尾渦的擾動(dòng)作用,使得渦3 渦核向下偏移.渦核負(fù)壓區(qū)遠(yuǎn)離尾部以及尾部斜面的分離渦的消失,都會(huì)增加尾部近壁面的壓力,從而使得模型頭部和尾部的壓差阻力減小,達(dá)到一定的減阻效果.

    圖7 時(shí)均速度流線圖(1)Fig.7 Time-mean velocity flow diagrams(1)

    Q 準(zhǔn)則等值面通??梢杂脕盹@示鈍體模型的三維渦系結(jié)構(gòu),圖8 對(duì)比了SAE 模型和附加VGS 的SAE 模型Q 準(zhǔn)則等值面,展示了兩種模型尾流區(qū)的三維渦結(jié)構(gòu).通過對(duì)比圖8 中兩種模型的局部放大圖可以發(fā)現(xiàn),相比于SAE 模型尾部斜面形成的分離渦,附加VGS 的模型產(chǎn)生了一系列抑制尾部斜面分離渦形成的小渦.如圖8(a)所示,當(dāng)流體從模型頂部流到尾部斜面時(shí),剛開始的分離屬于在平面上的二維形式的分離,然后受C 柱拖拽渦的影響變成更加復(fù)雜的三維分離流.在與原模型相對(duì)應(yīng)的位置,VGS后面的小渦使二維分離消失,同時(shí)近壁低速區(qū)域明顯減小,說明VGS 結(jié)構(gòu)增加了下游近壁區(qū)域內(nèi)的動(dòng)能,使近壁區(qū)邊界層下方和邊界層上方速度梯度減小,延遲了流體分離.

    圖8 兩種模型的瞬態(tài)Q準(zhǔn)則(1)Fig.8 Transient Q criterion for the two models(1)

    2.3 阻力系數(shù)對(duì)比分析

    在SAE模型尾部斜面上邊沿上添加的VGS模型與原模型的氣動(dòng)阻力系數(shù)對(duì)比結(jié)果如表2 所示,結(jié)果表明添加VGS 有一定的減阻效果,減阻率約為2.41%.

    表2 SAE模型和附加VGS模型的阻力系數(shù)和減阻率Tab.2 Drag coefficients and drag reduction rate for SAE model and SAE model with VGS

    3 RTS減阻效果研究

    3.1 RTS結(jié)構(gòu)尺寸

    本文在模型尾部添加了如圖9所示的V 型溝槽.圖9 中展示了尾部垂直面添加V 型溝槽的SAE 模型的全局視圖和尾部垂直面的放大視圖,以及溝槽的局部放大視圖以及尺寸參數(shù).其中溝槽間距S=3 mm,深度H=3 mm.

    圖9 尾部垂直面添加V型溝槽的SAE模型Fig.9 SAE model with V-riblets added to the vertical surface of the tail

    3.2 尾流結(jié)構(gòu)對(duì)比分析

    圖10 為SAE 模型和附加RTS 的模型在z=0 截面上的流向時(shí)均速度流線.與VGS 情況類似,圖10(b)中尾部垂直面上的非光滑表面使得近尾流區(qū)的兩渦中心位置向下游移動(dòng),特別是上渦的渦核中心,說明尾部面上的RTS 使得行李艙蓋尾沿處的分離延遲,在靠近尾部面上邊沿處還產(chǎn)生了一個(gè)衍生渦.另外我們從圖中發(fā)現(xiàn)尾部垂直面上的溝槽還影響了尾部斜面上分離渦的形成.

    圖10 時(shí)均速度流線圖(2)Fig.10 Time-mean velocity flow diagrams(2)

    圖11 進(jìn)一步對(duì)溝槽內(nèi)部流動(dòng)情況與光滑表面流動(dòng)情況進(jìn)行了比較,由圖可知,在圖11(b)尾部面溝槽內(nèi)部存在許多小渦,這些小渦阻礙了更外層的高速流體與近壁面直接接觸,能夠降低摩擦阻力.同時(shí),V 型溝槽能夠限制流體的橫向運(yùn)動(dòng),使得流動(dòng)更加穩(wěn)定,相應(yīng)減小了湍流猝發(fā)強(qiáng)度.

    圖11 時(shí)均速度絕對(duì)值流線圖Fig.11 Flow diagram of absolute value of time-mean velocity

    SAE 模型和附加RTS 模型的Q 準(zhǔn)則等值面對(duì)比如圖12 所示,其展示了兩種模型尾流區(qū)的三維渦結(jié)構(gòu).圖12(b)模型尾部斜面上渦流強(qiáng)度相較于原模型有所減小,近壁面低速區(qū)域范圍也相對(duì)減小.再對(duì)比尾部面上的渦結(jié)構(gòu),附加RTS 的模型尾部面上有許多條狀的渦結(jié)構(gòu),溝槽內(nèi)的這些渦結(jié)構(gòu)能夠限制近壁面更外層的高速流體與壁面直接接觸,從而有利于減小摩擦阻力.

    圖12 兩種模型的瞬態(tài)Q準(zhǔn)則(2)Fig.12 Transient Q criterion for the two models(2)

    3.3 阻力系數(shù)對(duì)比分析

    在SAE 模型尾部斜面上邊沿上添加RTS 的模型與原模型的氣動(dòng)阻力系數(shù)對(duì)比結(jié)果如表3 所示,結(jié)果表明添加RTS 有一定的減阻效果,達(dá)到2.76%的減阻率.

    表3 SAE原模型和附加RTS模型的阻力系數(shù)和減阻率Tab.3 Drag coefficients and drag reduction rate for SAE model and SAE model with RTS

    4 VGS和RTS對(duì)氣動(dòng)噪聲的影響研究

    由于車模后傾斜面頂部和尾部添加了VGS 和RTS,尾流結(jié)構(gòu)會(huì)產(chǎn)生變化,從而影響氣流脈動(dòng)頻率及脈動(dòng)壓力,所以會(huì)使得聲源的分布以及強(qiáng)弱發(fā)生變化.減阻裝置在減阻的同時(shí)對(duì)車輛噪聲產(chǎn)生了影響,本節(jié)通過模擬設(shè)置近場(chǎng)噪聲監(jiān)測(cè)點(diǎn)與遠(yuǎn)場(chǎng)麥克風(fēng)對(duì)被動(dòng)減阻產(chǎn)生的噪聲進(jìn)行研究.

    4.1 噪聲監(jiān)測(cè)點(diǎn)的布置

    通過計(jì)算獲得的z=0 截面的脈動(dòng)壓力如圖13所示,圖13(a)、13(b)、13(c)分別為原模型、附加VGS 的模型和附加RTS 的模型的結(jié)果.由圖可知,壓力劇烈變化的位置主要集中在尾流區(qū)域,特別是在尾部斜面上的分離渦位置和近尾流區(qū)域的兩個(gè)旋轉(zhuǎn)渦位置上,因此尾流區(qū)是主要的氣動(dòng)噪聲聲源.尾流區(qū)近場(chǎng)噪聲通過仿真中設(shè)置的監(jiān)測(cè)點(diǎn)獲得,監(jiān)測(cè)點(diǎn)的分布如圖14所示,它們都分布在圖13所示壓力脈動(dòng)劇烈的區(qū)域附近以及近尾流下游方向.圖14(a)展示了z=0截面上11個(gè)監(jiān)測(cè)點(diǎn)的位置,實(shí)際上每個(gè)監(jiān)測(cè)點(diǎn)都有左右兩個(gè)額外的監(jiān)測(cè)點(diǎn),如圖14(b)所示.

    圖13 z=0截面脈動(dòng)壓力分布Fig.13 Pulsating pressure distribution at section z=0

    圖14 流場(chǎng)仿真中監(jiān)測(cè)點(diǎn)的分布Fig.14 Distribution of monitoring points in flow field simulation

    另外,遠(yuǎn)場(chǎng)噪聲也是我們需要關(guān)注的.通過仿真模擬設(shè)置麥克風(fēng)來獲取遠(yuǎn)場(chǎng)輻射噪聲聲壓級(jí),遠(yuǎn)場(chǎng)麥克風(fēng)的位置分布在以SAE 模型為中心的半徑為10 m 的圓或者半圓上,如圖15 所示,分別為縱向、橫向和水平面上的麥克風(fēng)分布示意圖.其中縱、橫向截面上的監(jiān)測(cè)點(diǎn)P1~P7與P8~P14對(duì)應(yīng)于極坐標(biāo)中的0~180°,水平截面上的監(jiān)測(cè)點(diǎn)P15~P27對(duì)應(yīng)于極坐標(biāo)中的0~330°.

    圖15 監(jiān)測(cè)遠(yuǎn)場(chǎng)噪聲的麥克風(fēng)位置分布Fig.15 Microphone position distribution for monitoring far-field noise

    4.2 監(jiān)測(cè)點(diǎn)總聲壓級(jí)對(duì)比分析

    大量渦團(tuán)引起的氣流脈動(dòng)壓力就是汽車氣動(dòng)噪聲產(chǎn)生的根源,相對(duì)穩(wěn)定的渦結(jié)構(gòu)可以減少再附著渦量,降低由于氣流脈動(dòng)引起的氣動(dòng)噪聲.如前所述,添加在SAE 原模型斜面頂部的VGS 產(chǎn)生的小渦群抑制了尾部斜面上分離渦的形成,產(chǎn)生了相對(duì)穩(wěn)定的渦結(jié)構(gòu),而尾部添加的RTS 限制了流體的橫向運(yùn)動(dòng),使得流動(dòng)更加穩(wěn)定,這在一定程度上都能減小湍流猝發(fā)強(qiáng)度,降低氣流的脈動(dòng)頻率,由氣流脈動(dòng)產(chǎn)生的脈動(dòng)壓力也相應(yīng)降低,這將有利于氣動(dòng)噪聲的降低.圖16 為通過計(jì)算聲學(xué)方法所得的SAE 模型和附加VGS的模型以及附加RTS的模型各監(jiān)測(cè)點(diǎn)處的總聲壓級(jí)對(duì)比.從圖中可知,對(duì)于附加VGS 的模型,在1、3、4、6、10 左側(cè)監(jiān)測(cè)點(diǎn)處,1、2、3、4、5、7、10、11中間監(jiān)測(cè)點(diǎn)處,1、3、4、6、7、10 右側(cè)監(jiān)測(cè)點(diǎn)處總聲壓級(jí)均比原模型小,最大差值在點(diǎn)10 中間點(diǎn)處,為9.55 dB.而附加RTS 的模型在3、4、6 左側(cè)監(jiān)測(cè)點(diǎn)處,所有中間監(jiān)測(cè)點(diǎn)處,以及2、3、4、5、6、7、9、10 右側(cè)監(jiān)測(cè)點(diǎn)處的總聲壓級(jí)均比原模型小,最大差值在點(diǎn)7的中間點(diǎn)位置,為5.46 dB.所以SAE 模型添加VGS或RTS 后,均有一定的降噪效果.我們注意到,即使近場(chǎng)噪聲在大部分監(jiān)測(cè)點(diǎn)相較于原模型均有所降低,但同時(shí),也有部分監(jiān)測(cè)點(diǎn)噪聲高于原模型,這是由于VGS與RTS部分監(jiān)測(cè)點(diǎn)位置的渦團(tuán)數(shù)量要大于原模型.

    圖16 SAE模型和附加VGS、RTS的模型尾流區(qū)監(jiān)測(cè)點(diǎn)總聲壓級(jí)對(duì)比Fig.16 Comparison of total sound pressure levels in wake area monitoring points between SAE model and SAE model with VGS and RTS

    4.3 遠(yuǎn)場(chǎng)監(jiān)測(cè)點(diǎn)聲壓指向性分析

    添加VGS與RTS的模型與原模型遠(yuǎn)場(chǎng)輻射噪聲存在許多共同特征.如圖17 所示,在縱向截面以及橫向截面上,總聲壓級(jí)最大的位置出現(xiàn)在模型正上方;在水平截面上,總聲壓級(jí)最大的位置是模型兩側(cè),而頭部和尾部的總聲壓級(jí)相對(duì)較小.在縱向截面上,添加被動(dòng)減阻裝置對(duì)模型上方的噪聲減小有著積極的影響,而其他方向上的總聲壓級(jí)差別基本不大.而在橫向截面上,添加VGS 或RTS,各個(gè)方向的總聲壓級(jí)均小于原模型,且添加RTS 后的降噪效果優(yōu)于添加VGS.在水平方向上,除頭部和尾部差別不大外,其他方向上添加VGS 或RTS 后均有一定降噪效果,尤其是在模型兩側(cè),降噪效果明顯.

    圖17 SAE模型與添加VGS、RTS模型在各監(jiān)測(cè)點(diǎn)的總聲壓級(jí)指向性Fig.17 Total sound pressure level directivity at each monitoring point of SAE model and SAE model with VGS and RTS

    5 結(jié)論

    本文使用計(jì)算流體力學(xué)和計(jì)算氣動(dòng)聲學(xué)數(shù)值仿真分析方法,結(jié)合風(fēng)洞試驗(yàn)數(shù)據(jù)進(jìn)行準(zhǔn)確性和可靠性驗(yàn)證,以20°后背傾角階梯背式SAE模型為研究對(duì)象,分別研究了半圓柱形渦發(fā)生器(VGS)和溝槽(RTS)對(duì)模型外流場(chǎng)以及氣動(dòng)聲學(xué)的影響,并結(jié)合流場(chǎng)參數(shù)的變化簡單闡述了附加裝置的減阻降噪機(jī)理.本文的研究內(nèi)容和研究結(jié)果表明:

    1)相比于SAE 原模型,添加新型半圓柱型渦發(fā)生器及溝槽對(duì)氣動(dòng)減阻有積極的影響,其減阻率分別為2.41%、2.76%.

    2)在噪聲的控制方面,相較于SAE 模型,添加新型半圓柱型渦發(fā)生器及溝槽分別實(shí)現(xiàn)了最大9.55 dB、5.46 dB的降噪效果.

    3)本文對(duì)尾流區(qū)噪聲源研究的意義在于綜合評(píng)估被動(dòng)減阻裝置對(duì)尾流區(qū)氣動(dòng)噪聲源的影響,對(duì)協(xié)同減阻降噪具有一定的借鑒意義.

    猜你喜歡
    尾流聲壓級(jí)尾部
    機(jī)器噪聲平均聲壓級(jí)計(jì)算方法差異性實(shí)證研究
    船舶尾部響應(yīng)特性試驗(yàn)與計(jì)算
    超聲及磁共振診斷骶尾部藏毛竇1例
    一種計(jì)算消聲室聲壓級(jí)的新方法
    全新DXR mkll有源揚(yáng)聲器
    演藝科技(2019年4期)2019-03-30 03:21:46
    飛機(jī)尾流的散射特性與探測(cè)技術(shù)綜述
    錐形流量計(jì)尾流流場(chǎng)分析
    彎式尾部接頭注塑模具設(shè)計(jì)
    中國塑料(2015年4期)2015-10-14 01:09:32
    水面艦船風(fēng)尾流效應(yīng)減弱的模擬研究
    Diodes1.9W D類音頻放大器提供高聲壓級(jí)水平并延長電池壽命
    国产精品免费大片| 午夜久久久在线观看| 久久中文字幕人妻熟女| 一进一出好大好爽视频| 中文字幕人妻丝袜一区二区| 啦啦啦免费观看视频1| 国产日韩欧美在线精品| 亚洲精品久久成人aⅴ小说| 少妇精品久久久久久久| 亚洲精品久久成人aⅴ小说| 国产不卡一卡二| 午夜福利一区二区在线看| 日韩 欧美 亚洲 中文字幕| 一区福利在线观看| 国产av一区二区精品久久| 考比视频在线观看| 99国产精品免费福利视频| 男女下面插进去视频免费观看| 欧美日韩视频精品一区| 国产精品亚洲av一区麻豆| 老汉色av国产亚洲站长工具| 一区二区三区激情视频| 久久国产亚洲av麻豆专区| 777米奇影视久久| 黄色丝袜av网址大全| 香蕉久久夜色| 日韩免费av在线播放| 欧美+亚洲+日韩+国产| 天天躁狠狠躁夜夜躁狠狠躁| 国产不卡av网站在线观看| 国产一卡二卡三卡精品| 欧美日韩亚洲综合一区二区三区_| 亚洲熟女精品中文字幕| 一区二区av电影网| 久久久久视频综合| 最新美女视频免费是黄的| 精品国内亚洲2022精品成人 | 在线观看免费视频网站a站| 国产欧美日韩一区二区三| 日韩欧美一区二区三区在线观看 | 国产无遮挡羞羞视频在线观看| 天天躁夜夜躁狠狠躁躁| 国产一区二区三区在线臀色熟女 | 汤姆久久久久久久影院中文字幕| 欧美久久黑人一区二区| 悠悠久久av| 欧美变态另类bdsm刘玥| 99精品欧美一区二区三区四区| 国产伦理片在线播放av一区| 国产一区二区三区综合在线观看| 满18在线观看网站| 香蕉国产在线看| 午夜福利在线免费观看网站| 日韩 欧美 亚洲 中文字幕| 亚洲国产看品久久| 人人妻人人爽人人添夜夜欢视频| 一二三四社区在线视频社区8| 成人亚洲精品一区在线观看| 欧美精品一区二区大全| 男女下面插进去视频免费观看| 精品国产乱码久久久久久小说| 国产区一区二久久| 少妇 在线观看| 下体分泌物呈黄色| 女同久久另类99精品国产91| 日韩制服丝袜自拍偷拍| 国产日韩欧美在线精品| 久久精品国产a三级三级三级| 久久人人97超碰香蕉20202| 国产成人系列免费观看| 91精品三级在线观看| 日本a在线网址| 中亚洲国语对白在线视频| 成人国语在线视频| 怎么达到女性高潮| 黄色怎么调成土黄色| 色婷婷av一区二区三区视频| 正在播放国产对白刺激| 成人精品一区二区免费| 久久青草综合色| 一边摸一边抽搐一进一出视频| 下体分泌物呈黄色| 乱人伦中国视频| 两个人免费观看高清视频| 亚洲,欧美精品.| 午夜激情av网站| 三级毛片av免费| 99热网站在线观看| 91老司机精品| 亚洲国产欧美网| 精品少妇久久久久久888优播| 最新美女视频免费是黄的| 国产成人精品久久二区二区免费| 国产老妇伦熟女老妇高清| 亚洲欧美精品综合一区二区三区| 人人妻,人人澡人人爽秒播| 亚洲精品一二三| 国产精品久久久久久人妻精品电影 | 99热国产这里只有精品6| 51午夜福利影视在线观看| 午夜福利在线免费观看网站| 欧美成狂野欧美在线观看| 午夜成年电影在线免费观看| 男女床上黄色一级片免费看| 亚洲av美国av| 肉色欧美久久久久久久蜜桃| 99riav亚洲国产免费| www.精华液| 99久久99久久久精品蜜桃| 热re99久久精品国产66热6| 热99国产精品久久久久久7| 欧美精品av麻豆av| 大片电影免费在线观看免费| 日本wwww免费看| 久久ye,这里只有精品| 中文字幕色久视频| 精品国产一区二区久久| 亚洲国产中文字幕在线视频| 亚洲第一青青草原| 午夜福利视频在线观看免费| 天堂动漫精品| 国产欧美日韩一区二区三区在线| 人人澡人人妻人| 下体分泌物呈黄色| www.自偷自拍.com| 国产在线视频一区二区| 亚洲精品自拍成人| 黄色片一级片一级黄色片| 欧美另类亚洲清纯唯美| 丰满迷人的少妇在线观看| 超碰成人久久| 蜜桃国产av成人99| 一个人免费看片子| 99在线人妻在线中文字幕 | 一个人免费在线观看的高清视频| 免费观看av网站的网址| 国产激情久久老熟女| 国产精品亚洲一级av第二区| 啦啦啦免费观看视频1| 精品乱码久久久久久99久播| 色婷婷久久久亚洲欧美| 日韩免费高清中文字幕av| 免费在线观看视频国产中文字幕亚洲| 老汉色∧v一级毛片| 涩涩av久久男人的天堂| 午夜福利,免费看| 成人免费观看视频高清| 久久久久久久精品吃奶| 久久九九热精品免费| 免费av中文字幕在线| 亚洲欧美日韩另类电影网站| 精品亚洲成国产av| 五月开心婷婷网| 美女主播在线视频| 两人在一起打扑克的视频| 美女午夜性视频免费| 国产主播在线观看一区二区| 捣出白浆h1v1| 日韩欧美一区二区三区在线观看 | 精品欧美一区二区三区在线| √禁漫天堂资源中文www| 一级毛片精品| 热re99久久国产66热| 黄网站色视频无遮挡免费观看| 日韩成人在线观看一区二区三区| videosex国产| 99久久国产精品久久久| 黄色怎么调成土黄色| 在线观看www视频免费| 巨乳人妻的诱惑在线观看| 久久精品国产亚洲av高清一级| 亚洲一码二码三码区别大吗| 亚洲国产欧美一区二区综合| 两性夫妻黄色片| 免费一级毛片在线播放高清视频 | 黄色视频在线播放观看不卡| 国产xxxxx性猛交| 岛国毛片在线播放| 夜夜爽天天搞| 免费久久久久久久精品成人欧美视频| 欧美日韩国产mv在线观看视频| 又紧又爽又黄一区二区| 亚洲美女黄片视频| www.999成人在线观看| 精品福利永久在线观看| 操出白浆在线播放| 看免费av毛片| 丝袜在线中文字幕| 亚洲av成人不卡在线观看播放网| 99久久人妻综合| 免费不卡黄色视频| 亚洲色图 男人天堂 中文字幕| 最黄视频免费看| 欧美中文综合在线视频| 一区二区三区国产精品乱码| 国产成人欧美| 精品国产一区二区三区四区第35| 中文字幕高清在线视频| 欧美黑人精品巨大| 国产成人av教育| 国产老妇伦熟女老妇高清| 黄片小视频在线播放| 极品人妻少妇av视频| 岛国在线观看网站| 亚洲久久久国产精品| 久久久久久亚洲精品国产蜜桃av| 亚洲欧美激情在线| 久久中文字幕人妻熟女| 久久人人爽av亚洲精品天堂| 国产精品香港三级国产av潘金莲| 精品国产一区二区久久| 丰满饥渴人妻一区二区三| 人人妻人人爽人人添夜夜欢视频| 亚洲av日韩在线播放| 国产单亲对白刺激| 亚洲精品久久午夜乱码| 满18在线观看网站| 香蕉国产在线看| 日韩欧美免费精品| 男人操女人黄网站| 制服人妻中文乱码| 国产精品1区2区在线观看. | 免费在线观看完整版高清| 成人免费观看视频高清| av福利片在线| 欧美久久黑人一区二区| 侵犯人妻中文字幕一二三四区| 老熟妇仑乱视频hdxx| 日本精品一区二区三区蜜桃| 国产精品成人在线| 免费在线观看日本一区| 欧美黑人精品巨大| 首页视频小说图片口味搜索| 18禁观看日本| 男女边摸边吃奶| 精品国产国语对白av| 国产欧美日韩一区二区三| 亚洲人成伊人成综合网2020| 欧美成人午夜精品| 满18在线观看网站| 亚洲欧美一区二区三区久久| 日本黄色日本黄色录像| 亚洲国产中文字幕在线视频| 色在线成人网| 国产欧美日韩综合在线一区二区| 国产日韩欧美视频二区| 中文字幕另类日韩欧美亚洲嫩草| 久久99热这里只频精品6学生| 国产精品麻豆人妻色哟哟久久| 后天国语完整版免费观看| 不卡av一区二区三区| 欧美黑人欧美精品刺激| 少妇的丰满在线观看| 另类精品久久| 亚洲一区二区三区欧美精品| 欧美 亚洲 国产 日韩一| 精品少妇黑人巨大在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 在线亚洲精品国产二区图片欧美| 一边摸一边抽搐一进一小说 | 午夜福利视频在线观看免费| 别揉我奶头~嗯~啊~动态视频| 两性午夜刺激爽爽歪歪视频在线观看 | 精品高清国产在线一区| 欧美在线黄色| 女人精品久久久久毛片| 每晚都被弄得嗷嗷叫到高潮| 午夜福利在线免费观看网站| 人人妻人人澡人人看| 亚洲精品乱久久久久久| 精品少妇黑人巨大在线播放| av超薄肉色丝袜交足视频| 国产精品秋霞免费鲁丝片| cao死你这个sao货| 日本av免费视频播放| 精品国产乱子伦一区二区三区| 欧美国产精品va在线观看不卡| 亚洲国产成人一精品久久久| 日韩欧美一区视频在线观看| 中国美女看黄片| 日本wwww免费看| 女人爽到高潮嗷嗷叫在线视频| 成年动漫av网址| 日本撒尿小便嘘嘘汇集6| 欧美 亚洲 国产 日韩一| 好男人电影高清在线观看| 国产成人欧美在线观看 | 久久精品91无色码中文字幕| 亚洲国产欧美在线一区| 手机成人av网站| avwww免费| 成人亚洲精品一区在线观看| 大香蕉久久成人网| 国产高清国产精品国产三级| 悠悠久久av| 欧美精品人与动牲交sv欧美| 我要看黄色一级片免费的| 亚洲精品国产精品久久久不卡| 国产高清国产精品国产三级| 国产av一区二区精品久久| 亚洲成国产人片在线观看| 每晚都被弄得嗷嗷叫到高潮| 免费观看a级毛片全部| 欧美精品亚洲一区二区| videosex国产| 国产成人免费无遮挡视频| 80岁老熟妇乱子伦牲交| 免费av中文字幕在线| 少妇精品久久久久久久| 三级毛片av免费| 国产一区二区三区视频了| 亚洲av国产av综合av卡| 精品福利观看| 国产成人精品久久二区二区91| 亚洲成人手机| 国产亚洲欧美在线一区二区| 一级毛片精品| 99精品欧美一区二区三区四区| 国产97色在线日韩免费| 成人特级黄色片久久久久久久 | 久久中文字幕人妻熟女| 两个人免费观看高清视频| 黑丝袜美女国产一区| 操美女的视频在线观看| 亚洲自偷自拍图片 自拍| 两个人免费观看高清视频| 午夜福利在线观看吧| av国产精品久久久久影院| 老熟女久久久| 国产欧美日韩一区二区三区在线| 十八禁网站网址无遮挡| 国产精品秋霞免费鲁丝片| 99国产精品一区二区三区| 精品视频人人做人人爽| 成人亚洲精品一区在线观看| 亚洲精品美女久久av网站| 久久亚洲真实| 精品国内亚洲2022精品成人 | 精品午夜福利视频在线观看一区 | 免费看a级黄色片| 91av网站免费观看| 欧美乱妇无乱码| 精品高清国产在线一区| 日韩欧美一区二区三区在线观看 | 在线天堂中文资源库| 亚洲七黄色美女视频| 一区二区av电影网| 中文字幕av电影在线播放| 亚洲精品中文字幕在线视频| 亚洲熟妇熟女久久| 免费在线观看黄色视频的| 日日爽夜夜爽网站| 免费黄频网站在线观看国产| 国产福利在线免费观看视频| 亚洲精品久久午夜乱码| 精品亚洲成a人片在线观看| 中文字幕另类日韩欧美亚洲嫩草| 两人在一起打扑克的视频| 国产成+人综合+亚洲专区| 波多野结衣av一区二区av| a级毛片在线看网站| 国产精品美女特级片免费视频播放器 | 亚洲精品一二三| 久久这里只有精品19| 少妇 在线观看| 最近最新中文字幕大全免费视频| 99国产精品一区二区三区| 国产一区二区三区在线臀色熟女 | 国产成人啪精品午夜网站| 欧美av亚洲av综合av国产av| 欧美黄色淫秽网站| 在线亚洲精品国产二区图片欧美| 最黄视频免费看| 不卡一级毛片| 黄片大片在线免费观看| 一区二区三区激情视频| 午夜福利在线观看吧| 在线十欧美十亚洲十日本专区| 久久中文字幕人妻熟女| 国产成+人综合+亚洲专区| 国产三级黄色录像| 亚洲 欧美一区二区三区| 日本精品一区二区三区蜜桃| 淫妇啪啪啪对白视频| kizo精华| 黄色丝袜av网址大全| av在线播放免费不卡| 精品国产一区二区久久| 国产在线一区二区三区精| 中文字幕人妻熟女乱码| 久久久久久久大尺度免费视频| 免费看十八禁软件| 777久久人妻少妇嫩草av网站| 久久99一区二区三区| 我要看黄色一级片免费的| 极品教师在线免费播放| 亚洲专区字幕在线| 大型av网站在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品91无色码中文字幕| 免费观看人在逋| 国产高清国产精品国产三级| 亚洲少妇的诱惑av| a级毛片在线看网站| 亚洲欧美一区二区三区久久| 老司机靠b影院| 国产野战对白在线观看| 叶爱在线成人免费视频播放| 亚洲国产中文字幕在线视频| 国产深夜福利视频在线观看| 欧美精品一区二区免费开放| 亚洲精品成人av观看孕妇| 欧美在线一区亚洲| 成人永久免费在线观看视频 | 日韩中文字幕视频在线看片| 亚洲av成人不卡在线观看播放网| 如日韩欧美国产精品一区二区三区| 亚洲av成人一区二区三| 日韩 欧美 亚洲 中文字幕| 国产欧美日韩综合在线一区二区| 亚洲欧美激情在线| 亚洲精品中文字幕一二三四区 | 中文字幕人妻丝袜制服| 午夜福利视频精品| 亚洲成av片中文字幕在线观看| 国产欧美日韩综合在线一区二区| 国产精品98久久久久久宅男小说| 亚洲成人国产一区在线观看| 午夜久久久在线观看| 亚洲自偷自拍图片 自拍| 99精品欧美一区二区三区四区| 亚洲第一欧美日韩一区二区三区 | 最近最新免费中文字幕在线| 欧美 亚洲 国产 日韩一| 成在线人永久免费视频| 啦啦啦中文免费视频观看日本| 国产精品影院久久| 成年人黄色毛片网站| 老司机福利观看| 国产深夜福利视频在线观看| 中文字幕人妻熟女乱码| 飞空精品影院首页| 国产黄频视频在线观看| 久久久久久久大尺度免费视频| 久久中文字幕人妻熟女| av天堂久久9| 91成人精品电影| 中国美女看黄片| 亚洲精华国产精华精| 怎么达到女性高潮| 欧美黑人精品巨大| 亚洲精品乱久久久久久| 亚洲全国av大片| 99riav亚洲国产免费| 国产老妇伦熟女老妇高清| 国产男女超爽视频在线观看| 丁香六月欧美| 深夜精品福利| 男女高潮啪啪啪动态图| 久热这里只有精品99| 两个人看的免费小视频| 高清视频免费观看一区二区| 91九色精品人成在线观看| 国产精品久久久久久人妻精品电影 | 国产一区有黄有色的免费视频| 亚洲欧洲日产国产| 午夜福利,免费看| 精品高清国产在线一区| 精品福利永久在线观看| 悠悠久久av| 一本—道久久a久久精品蜜桃钙片| 久久精品国产亚洲av香蕉五月 | 亚洲成人免费电影在线观看| 国产极品粉嫩免费观看在线| 亚洲精品在线观看二区| 美女扒开内裤让男人捅视频| 欧美黑人欧美精品刺激| 午夜福利影视在线免费观看| 又大又爽又粗| 好男人电影高清在线观看| 久久久久视频综合| 午夜福利,免费看| 中文字幕高清在线视频| cao死你这个sao货| 飞空精品影院首页| 国产97色在线日韩免费| 国产精品国产av在线观看| 啦啦啦视频在线资源免费观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美精品啪啪一区二区三区| 高清黄色对白视频在线免费看| 久久精品亚洲精品国产色婷小说| 午夜激情av网站| 午夜福利免费观看在线| 黄色怎么调成土黄色| 久久狼人影院| 成人国产av品久久久| 久久人人97超碰香蕉20202| 99香蕉大伊视频| 午夜视频精品福利| 成人黄色视频免费在线看| 成人特级黄色片久久久久久久 | 免费在线观看日本一区| 99精品欧美一区二区三区四区| 人人澡人人妻人| 国产欧美日韩精品亚洲av| 亚洲欧美日韩高清在线视频 | netflix在线观看网站| 亚洲七黄色美女视频| 久久久久久免费高清国产稀缺| 无人区码免费观看不卡 | 日本av手机在线免费观看| 免费高清在线观看日韩| 黄色成人免费大全| 别揉我奶头~嗯~啊~动态视频| 日韩成人在线观看一区二区三区| 久久热在线av| 国内毛片毛片毛片毛片毛片| 天堂俺去俺来也www色官网| 精品久久久久久电影网| 考比视频在线观看| 国产1区2区3区精品| 久久久久久久精品吃奶| 蜜桃国产av成人99| 亚洲国产欧美在线一区| 日本精品一区二区三区蜜桃| 无限看片的www在线观看| 黄片大片在线免费观看| 国产日韩欧美亚洲二区| av一本久久久久| 久久久国产精品麻豆| 九色亚洲精品在线播放| 视频区图区小说| 国产真人三级小视频在线观看| 看免费av毛片| 亚洲国产精品一区二区三区在线| 999精品在线视频| 欧美亚洲 丝袜 人妻 在线| 麻豆成人av在线观看| 一级,二级,三级黄色视频| 久久精品国产a三级三级三级| 国产亚洲午夜精品一区二区久久| 美女高潮到喷水免费观看| 亚洲午夜精品一区,二区,三区| 欧美人与性动交α欧美精品济南到| 一个人免费看片子| 一级毛片女人18水好多| 久久久国产精品麻豆| 757午夜福利合集在线观看| av超薄肉色丝袜交足视频| 久久久精品区二区三区| 亚洲人成电影免费在线| 久久久久国内视频| 精品国产国语对白av| 肉色欧美久久久久久久蜜桃| 中文字幕av电影在线播放| 国产午夜精品久久久久久| 精品卡一卡二卡四卡免费| 丝袜美腿诱惑在线| 我要看黄色一级片免费的| 国产福利在线免费观看视频| 黄色成人免费大全| 久久性视频一级片| 老熟妇仑乱视频hdxx| 精品一区二区三区四区五区乱码| 国产成人系列免费观看| 国产亚洲午夜精品一区二区久久| a级片在线免费高清观看视频| 啦啦啦中文免费视频观看日本| 露出奶头的视频| 美国免费a级毛片| av超薄肉色丝袜交足视频| 韩国精品一区二区三区| 欧美中文综合在线视频| 老司机影院毛片| 欧美成人午夜精品| 老司机午夜福利在线观看视频 | 久久这里只有精品19| 国产高清激情床上av| 午夜91福利影院| 日本av免费视频播放| avwww免费| 精品国产国语对白av| 欧美精品一区二区免费开放| 啦啦啦免费观看视频1| 中文字幕高清在线视频| 汤姆久久久久久久影院中文字幕| 无遮挡黄片免费观看| 日本av手机在线免费观看| 欧美在线一区亚洲| 午夜免费成人在线视频| 9色porny在线观看| 天天操日日干夜夜撸| 好男人电影高清在线观看| 欧美精品av麻豆av| 亚洲精品中文字幕一二三四区 | 久久精品国产99精品国产亚洲性色 | 久久久久国内视频| 无遮挡黄片免费观看| 日日摸夜夜添夜夜添小说| 麻豆av在线久日| tube8黄色片| 国产精品电影一区二区三区 | 夜夜夜夜夜久久久久| 在线十欧美十亚洲十日本专区| 日韩视频在线欧美| 老鸭窝网址在线观看| 国产精品98久久久久久宅男小说| 高清毛片免费观看视频网站 | 高清黄色对白视频在线免费看| 久久婷婷成人综合色麻豆| 桃花免费在线播放| 在线永久观看黄色视频| 久久国产亚洲av麻豆专区| av超薄肉色丝袜交足视频|