亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        B柱擾流器對風(fēng)振噪聲的影響及其機理 *

        2023-10-12 02:17:06陳小勇郭正鑫王永亮
        汽車工程 2023年9期
        關(guān)鍵詞:風(fēng)振聲壓級空腔

        陳小勇,姚 璐,郭正鑫,劉 浩,羅 挺,王永亮

        (東風(fēng)汽車集團有限公司技術(shù)中心,武漢 430056)

        前言

        在汽車乘坐舒適性指標(biāo)中,風(fēng)噪大一直被眾多消費者所抱怨,特別是電動車的風(fēng)噪,因為沒有發(fā)動機噪聲的掩蓋效應(yīng)。在汽車開窗行駛過程中,車內(nèi)空氣會產(chǎn)生一種低頻但波動劇烈的振動噪聲[1],即風(fēng)振噪聲,該噪聲會嚴(yán)重影響車內(nèi)乘員的駕乘體驗。因此對風(fēng)振噪聲的特性進行研究并對其進行控制具有重要意義。

        風(fēng)振噪聲可分為天窗風(fēng)振和側(cè)窗風(fēng)振兩種,其中天窗風(fēng)振已經(jīng)有較多研究,且得到了較好控制[2-4]。側(cè)窗風(fēng)振的控制相對較難,研究表明單開后側(cè)車窗的風(fēng)振噪聲要顯著大于單開前側(cè)車窗[5-6],車內(nèi)各點的風(fēng)振噪聲幅值和頻率是一致的[7]。側(cè)窗風(fēng)振是由自由剪切層的自激振蕩導(dǎo)致車內(nèi)產(chǎn)生壓力脈動,壓力脈動再與車內(nèi)空腔發(fā)生赫姆霍茲共振而導(dǎo)致的[7]。已有研究總結(jié)了不同車速、不同側(cè)窗開啟工況下的風(fēng)振噪聲變化規(guī)律[5,8-10],為風(fēng)振噪聲的控制提供參考。對于側(cè)窗風(fēng)振的控制,眾多研究都提出了被動控制方案,例如B 柱擾流器、車窗立柱、后視鏡結(jié)構(gòu)等[5,7,11-14]。風(fēng)振研究一般以樣機風(fēng)洞試驗和數(shù)值模擬相結(jié)合的方式進行[15],研究表明瞬態(tài)大渦模擬能較準(zhǔn)確地捕捉風(fēng)振噪聲[7],不同湍流度對側(cè)窗風(fēng)振噪聲會有一定影響[16]。Mendonca 等[17]的研究發(fā)現(xiàn)簡易縮比模型可以預(yù)測全尺寸整車模型的風(fēng)振噪聲。以上研究用仿真和試驗的方法對風(fēng)振噪聲的產(chǎn)生機理、不同工況下風(fēng)振噪聲的特性以及風(fēng)振噪聲的控制方法進行了探索,為理解和控制汽車風(fēng)振噪聲提供了思路。

        在這些研究的基礎(chǔ)上,進一步探索不同形式B柱擾流器對后側(cè)窗風(fēng)振噪聲的影響,且揭示不同擾流器對風(fēng)振的影響機理。首先,用不同車速下的風(fēng)振仿真與試驗結(jié)果對標(biāo)驗證了仿真方法的有效性;然后結(jié)合風(fēng)振仿真中的速度場和壓力場解釋了風(fēng)振噪聲的產(chǎn)生機理,同時提出了風(fēng)振噪聲的降噪策略;接著通過仿真方法探究了4 種B 柱擾流器降噪方案對后側(cè)窗風(fēng)振噪聲的影響,并揭示了不同B 柱擾流器對風(fēng)振的影響機理;最后用試驗驗證了兩種B 柱擾流器方案的有效性。本文研究內(nèi)容對于進一步理解和控制風(fēng)振噪聲具有一定理論意義和工程應(yīng)用價值。

        1 風(fēng)振仿真方法的有效性驗證

        擬采用流場仿真和試驗相結(jié)合的方法探究不同形式B 柱擾流器對風(fēng)振的影響及其機理,因此首先需要將風(fēng)振仿真結(jié)果和試驗結(jié)果進行對比,以驗證風(fēng)振仿真方法的有效性。

        1.1 仿真方法

        將整車外表面和內(nèi)飾的封閉幾何模型導(dǎo)入Star-CCM+流體仿真軟件,幾何模型的左后車窗全開,以模擬產(chǎn)生風(fēng)振的工況。在Star-CCM+中新建計算域以模擬風(fēng)洞,計算域大小為13 倍車長,12 倍車寬,6 倍車高,以保證計算域邊界處的流場不受汽車周圍流場的影響。

        接著在計算域和整車內(nèi)外飾之間生成用于仿真計算的有限元體網(wǎng)格。為保證計算精度,同時提高計算效率,對體網(wǎng)格進行分區(qū)細(xì)化。距車身0.15 m內(nèi)的網(wǎng)格目標(biāo)尺寸為0.01 m,距車身0.15 ~ 1.3 m以內(nèi)的網(wǎng)格目標(biāo)尺寸為0.05 m,距車身1.3 m 以上的網(wǎng)格目標(biāo)尺寸為0.5 m。最終的體網(wǎng)格數(shù)量約為1 875萬,如圖1所示。

        圖1 仿真計算體網(wǎng)格模型

        設(shè)定計算域和車身內(nèi)外表面的邊界條件,計算域入口為速度入口(40/70/100 km/h),出口為壓力出口(標(biāo)準(zhǔn)大氣壓力)。車身內(nèi)外表面和車體下的計算域為不可滑移邊界(模擬固體界面),其他表面為可滑移邊界(模擬無邊界)。

        為提升計算效率,先以穩(wěn)態(tài)物理模型進行計算,收斂后再以瞬態(tài)可壓分離渦模擬(detached-eddy simulation, DES)模型進行計算。在瞬態(tài)計算過程中,監(jiān)測車內(nèi)4 個乘員耳點位置的壓力波動。瞬態(tài)計算時間步長0.001 s,總時長2 s,0.5 s 后壓力波動趨于穩(wěn)定,選取穩(wěn)定后的1~2 s時間段內(nèi)的壓力波動數(shù)據(jù),經(jīng)傅里葉變換后得到乘員耳點的聲壓級-頻率曲線。

        1.2 試驗方法

        風(fēng)振噪聲的峰值頻率很低,在20 Hz 內(nèi),且幅值較大,在100 dB 以上[5]。與外造型風(fēng)噪、加速噪聲和路噪等其他噪聲相比,風(fēng)振噪聲的頻率很低、量級很大,因此道路試驗的風(fēng)振測試結(jié)果基本不受其他噪聲的影響,可以用來與仿真結(jié)果進行對比。

        道路風(fēng)振試驗在直線水泥路上進行,汽車以定速(40/70/100 km/h)行駛,汽車左后車窗全開,與仿真工況保持一致,在車內(nèi)4 個乘員的耳點位置布置傳聲器,采集汽車行進過程中的壓力波動,采集時間為15 s,經(jīng)傅里葉變換轉(zhuǎn)換成聲壓級-頻率曲線。此外,因道路測試時,自然風(fēng)會對風(fēng)振噪聲產(chǎn)生一定影響,所以同一車速下須測量往返兩次數(shù)據(jù)。

        1.3 仿真與試驗結(jié)果對比

        70 km/h 車速下主駕乘員耳點的仿真與試驗結(jié)果曲線的對比如圖2 所示。從圖中可以看出,由于自然風(fēng)的影響,往返的數(shù)據(jù)存在差異,因此需要對往返數(shù)據(jù)取均值。

        圖2 主駕乘員耳點風(fēng)振噪聲(70 km/h)的仿真與試驗曲線

        表1 統(tǒng)計了各工況下各乘員耳點的仿真結(jié)果和往返取均值后的試驗結(jié)果,以及仿真和試驗結(jié)果之間的誤差,其中峰值頻率的誤差在1.7 Hz 以內(nèi),風(fēng)振噪聲峰值的誤差在3.9 dB 以內(nèi)。與所調(diào)研的最新研究[13]的精度(0~2 Hz,4.7~6 dB)相比,本文的峰值頻率精度與其相當(dāng),但峰值大小的精度要更高一些。這表明本文仿真方法的準(zhǔn)確性較好。

        表1 乘員耳點風(fēng)振噪聲的仿真與試驗結(jié)果對比

        2 風(fēng)振噪聲的產(chǎn)生機理與降噪策略

        文獻[8]中利用仿真方法分析了簡易空腔的風(fēng)振產(chǎn)生原理,即自激蕩原理和赫姆霍茲共振原理。接下來根據(jù)這兩個原理,結(jié)合本課題所做的汽車風(fēng)振仿真結(jié)果,對汽車風(fēng)振產(chǎn)生機理進行簡要分析,并由此得到風(fēng)振噪聲的降噪策略。

        2.1 風(fēng)振噪聲的自激蕩原理

        汽車在運動時,其表面因為不可滑移,會在周圍的空氣中形成一個邊界層。邊界層內(nèi)的空氣流速呈梯度變化,接觸車身的邊界和汽車速度一致,接觸空氣的邊界和空氣速度一致。

        在開啟的后車窗處,因為幾何形狀突變,導(dǎo)致空氣流場的邊界層分離,產(chǎn)生自由剪切層,如圖3 中后車窗位置的流場所示。

        圖3 一個風(fēng)振周期內(nèi)的流場變化

        汽車在運動時,車窗內(nèi)外流場存在速度差,這個速度差會導(dǎo)致自由剪切層不穩(wěn)定,促成湍流渦的產(chǎn)生,如圖3 中的紅色箭頭所示。產(chǎn)生的渦沿著車窗從前向后運動,一個周期后,渦撞擊車窗后緣,舊的渦消失,新的渦在車窗前緣生成。

        從0T到3T/8,渦在車窗前半部分運動,車外的空氣受渦的影響侵入車內(nèi),自由剪切層向車內(nèi)偏移,車內(nèi)的氣壓增大,如圖4所示。

        圖4 一個風(fēng)振周期內(nèi)的壓力場變化

        從T/2 到7T/8,渦在車窗后半部分,車內(nèi)的氣壓已經(jīng)增大到一定程度,車內(nèi)的空氣受渦的影響開始撤出車內(nèi),自由剪切層向車外偏移,車內(nèi)的氣壓減小。

        在T時刻,渦撞擊車窗后緣,惡化了空腔前緣剪切層的不穩(wěn)定性,誘發(fā)前緣生成一個新的渦結(jié)構(gòu),形成渦運動反饋機制。

        上述的渦的生成和反饋機制即為汽車風(fēng)振的自激蕩原理,兩個渦生成的時間差即為自激蕩的周期,其倒數(shù)為自激蕩頻率。

        2.2 風(fēng)振噪聲的赫姆霍茲共振原理

        對于汽車內(nèi)部的幾何空腔結(jié)構(gòu),存在一個空腔固有頻率??涨坏墓逃蓄l率由空腔的幾何形狀和空腔內(nèi)的空氣屬性決定。文獻[8]中給出了簡單空腔的固有頻率經(jīng)驗公式,即

        式中:c為聲速;A為空腔開口處面積;hc為開口厚度;V為空腔體積。

        本文中空腔開口處面積A為0.23 m2,開口厚度hc為0.08 m,空腔體積V為3.21 m3。根據(jù)式(1)可以計算得到汽車內(nèi)部空腔的固有頻率f為19.8 Hz。隨著空氣流速的變化,汽車風(fēng)振自激蕩的周期和頻率也會相應(yīng)變化。當(dāng)自激蕩的頻率接近或等于汽車內(nèi)部空腔的固有頻率時,就會激發(fā)赫姆霍茲共振現(xiàn)象。

        根據(jù)表1 可以知道40、70、100 km/h 車速下的自激蕩頻率分別為8.7、14.0、18.2 Hz。由此可以推斷,40 km/h下的自激蕩頻率離汽車內(nèi)部空腔的固有頻率最遠(yuǎn),因此赫姆霍茲共振效果最弱,車內(nèi)氣壓脈動最小,聲壓級峰值最小。而100 km/h 下的自激蕩頻率離汽車內(nèi)部空腔的固有頻率最近,赫姆霍茲共振效果最強,車內(nèi)氣壓脈動最大,聲壓級峰值最大。

        2.3 風(fēng)振噪聲的降噪策略

        由汽車風(fēng)振的自激蕩原理可知,后側(cè)車窗內(nèi)外流場速度差形成的湍流渦是風(fēng)振的誘因。渦的大小和生成頻率決定了風(fēng)振噪聲峰值的頻率和大小。其他條件相同時,生成的渦越大,風(fēng)振峰值越大;渦的生成頻率越大,風(fēng)振峰值頻率越大;此外,根據(jù)赫姆霍茲共振原理,渦的生成頻率越接近汽車內(nèi)部空腔的固有頻率,風(fēng)振峰值也越大,反之亦然。

        根據(jù)以上分析可以得到的風(fēng)振噪聲降噪策略有:

        (1)減小圖3 中渦的強度,使流入車內(nèi)和從車內(nèi)流出的空氣量減少,可以通過在B 柱上設(shè)計擾流器實現(xiàn);

        (2)減小風(fēng)振的自激蕩頻率,削弱赫姆霍茲共振效果,可通過改變車窗的寬度實現(xiàn);

        (3)增大汽車內(nèi)部空腔的固有頻率,削弱赫姆霍茲共振效果,可通過改變車窗面積、車窗厚度和車內(nèi)空腔體積實現(xiàn)。

        在汽車的設(shè)計過程中,一旦到達樣車試制階段,車窗面積、寬度、厚度和車內(nèi)空腔等參數(shù)都已經(jīng)完全固化。此時若有風(fēng)振降噪的需求,便只能采用B 柱擾流器的降噪策略。因此,本文將對B 柱擾流器這一降噪策略進行深入研究,探討其對風(fēng)振噪聲的影響,并進一步分析其機理。

        因為不同車型的結(jié)構(gòu)不一樣,對B 柱擾流器的要求也不一樣,所以研究B 柱擾流器對風(fēng)振的影響及其機理可以為更有效的B 柱擾流器設(shè)計提供指導(dǎo)。

        3 B柱擾流器對風(fēng)振的影響及其機理

        3.1 不同形式B柱擾流器對風(fēng)振的影響

        為探究B 柱擾流器對風(fēng)振的影響及其機理,本文設(shè)計了4 種B 柱擾流器,擾流器的結(jié)構(gòu)如圖5 所示,安裝位置如圖6所示。

        圖5 B柱擾流器的結(jié)構(gòu)

        圖6 B柱擾流器的安裝位置(藍(lán)色零件為B柱)

        對安裝這4 種擾流器的車型進行70 km/h 車速下的風(fēng)振仿真,得到其風(fēng)振的聲壓級-頻率曲線。其中,主駕乘員耳點的風(fēng)振聲壓級-頻率曲線如圖7 所示。從曲線上可獲得聲壓級的峰值和峰值對應(yīng)的頻率,如表2所示。

        表2 4種擾流器方案的風(fēng)振峰值頻率和峰值

        圖7 主駕乘員耳點的風(fēng)振噪聲曲線

        從表2 可以計算得到相對于原結(jié)構(gòu),方案A、B、C、D對峰值和峰值頻率的影響,結(jié)果如表3所示。

        表3 4種擾流器方案對風(fēng)振噪聲的影響

        3.2 影響機理分析

        為進一步探究4種方案對風(fēng)振的影響機理,圖8展示了4 種不同方案的左后車窗處速度矢量圖。根據(jù)此前分析,減小圖3 中渦的強度,使流入車內(nèi)和從車內(nèi)流出的空氣量減少,可以降低風(fēng)振聲壓級峰值。而渦的強度可以用自由剪切層偏離B 柱外表面的最大距離大小(如圖8中距離d,以下簡稱偏移距離)來衡量。偏移距離越大,渦的強度越大;偏移距離越小,渦的強度越小。

        圖8 左后車窗處的速度矢量圖

        圖8中所示速度矢量圖的時刻均為圖3中的T/2時刻。因為從0T到3T/8,自由剪切層向車內(nèi)偏移;從T/2 到7T/8,自由剪切層向車外偏移。所以在T/2時刻,偏移距離最大,易于對比分析。對圖8 中5 種結(jié)構(gòu)的偏移距離進行測量,結(jié)果如表3所示。

        對比表3 中的聲壓級峰值和偏移距離,可以得出擾流器的形式對風(fēng)振的影響規(guī)律:

        (1)擾流器對偏移距離與對聲壓級峰值的影響呈正相關(guān)。即若擾流器增大了偏移距離,則渦的強度變大,自由剪切層振蕩加劇,風(fēng)振聲壓級增大;反之亦然。

        (2)擾流器對峰值的頻率基本無影響。

        (3)擾流器方案A 和B 能減小偏移距離和渦的強度,從而減小風(fēng)振聲壓級峰值;方案C 和D 會增大偏移距離和渦的強度,從而增大風(fēng)振聲壓級峰值。

        由以上結(jié)論可知,要想降低風(fēng)振聲壓級峰值,可以在汽車B 柱上設(shè)計方案A 或B 所示的擾流器結(jié)構(gòu),這兩種結(jié)構(gòu)可以減小偏移距離和渦的強度,從而減小風(fēng)振聲壓級峰值。此外,其他能減小偏移距離的擾流器結(jié)構(gòu),也都能起到降低風(fēng)振聲壓級峰值的效果。

        3.3 方案B的參數(shù)尋優(yōu)

        由表3 可知,方案B 的風(fēng)振降噪效果最好,達到4.1 dB,但降噪后的主駕風(fēng)振噪聲仍有121.2 dB,處于較高水平。為達到更優(yōu)的降噪效果,本文對方案B 的截面高度進行調(diào)節(jié),并仿真分析了新方案的風(fēng)振降噪效果。

        新方案的截面高度在已有方案B 的基礎(chǔ)上進行了加高,分別加高了5 mm(方案B-1)和10 mm(方案B-2),如圖9所示。新方案在B柱上的位置與圖6中的方案B位置保持一致。

        圖9 方案B、B-1和B-2的截面圖

        仿真得到方案B-1和B-2的降噪效果,如表4所示。與表2 中的原結(jié)構(gòu)風(fēng)振峰值相比,方案B-2 的風(fēng)振降噪效果顯著,達到了12.8 dB。

        表4 方案B-1和B-2風(fēng)振峰值頻率和峰值

        圖10 進一步對比了方案B、B-1、B-2 左后車窗處的速度矢量圖。由圖可知,方案B-1 的偏移距離與圖8 中的原方案偏移距離相等,即渦的強度和原方案相近,均為d,導(dǎo)致其降噪效果很小。方案B-2的自由剪切層出現(xiàn)了分層(如圖10 中的紅圈所示),上層氣流可以認(rèn)為在車內(nèi)循環(huán),只有下層會導(dǎo)致車內(nèi)外氣流交換,從而導(dǎo)致車內(nèi)壓力波動。因此,方案B-2 的偏移距離只測量下自由剪切層,即0.5d。該偏移距離顯著小于原方案的偏移距離d,即渦的強度顯著小于原方案,使其降噪效果明顯。以上結(jié)論與3.2節(jié)中的影響機理一致。

        圖10 方案B、B-1、B-2的速度矢量圖

        4 風(fēng)振噪聲降噪方案的試驗驗證

        仿真表明圖5 中方案A 和B 所示的擾流器結(jié)構(gòu)可以減小風(fēng)振聲壓級峰值,為進一步探究其實際改善效果,對這兩種方案進行試驗驗證。試驗的車型、車速與3.1節(jié)中的仿真保持一致,即車速為70 km/h,試驗方法與1.2 節(jié)一致。仿真與試驗的擾流器方案對比如圖11 所示。試驗方案的擾流器用3D 打印制作,仿真與試驗的方案結(jié)構(gòu)基本一致。

        圖11 仿真與試驗的擾流器方案

        試驗得到的主駕乘員耳點風(fēng)振噪聲的聲壓級-頻率曲線與仿真曲線對比如圖12 所示。從圖12 可獲得方案A 和B 對聲壓級峰值的試驗降噪效果和仿真降噪效果,結(jié)果如表5 所示。根據(jù)表中數(shù)據(jù),可以得到以下結(jié)論。

        表5 方案A和B的試驗與仿真降噪效果對比

        圖12 主駕乘員耳點的試驗與仿真風(fēng)振噪聲曲線對比

        (1)定性來看,試驗和仿真的方案都能減小聲壓級峰值,且試驗和仿真結(jié)果都表明方案B 的效果優(yōu)于方案A,因此試驗和仿真具有較好的一致性;

        (2)定量來看,方案A 和方案B 的試驗降噪效果均比仿真降噪效果小1 dB 左右,推測主要原因在于試驗結(jié)果是在道路上測試得到,而仿真結(jié)果是模擬的風(fēng)洞,因此道路測試時其他噪聲(路噪、胎噪等)難免會對風(fēng)噪產(chǎn)生一定影響。

        總而言之,試驗驗證結(jié)果表明,本文給出的方案A 和B 兩種降噪方案降低風(fēng)振聲壓級的效果比較明顯,且試驗的降噪效果與仿真的降噪效果一致性較好。

        5 結(jié)論

        探索不同形式B 柱擾流器對后側(cè)窗風(fēng)振噪聲大小的影響,且揭示了不同擾流器對風(fēng)振的影響機理。首先驗證了流場仿真方法的有效性;然后解釋了風(fēng)振噪聲的產(chǎn)生機理并提出風(fēng)振噪聲的降噪策略;接著探究4 種B 柱擾流器降噪方案對后側(cè)窗風(fēng)振噪聲大小的影響,并揭示其影響機理;最后用試驗驗證了兩種B柱擾流器方案的有效性。主要結(jié)論如下。

        (1)風(fēng)振仿真的峰值和峰值頻率誤差在3.9 dB和1.7 Hz以內(nèi),仿真方法的準(zhǔn)確性較好。

        (2)汽車風(fēng)振由自激蕩原理和赫姆霍茲共振原理導(dǎo)致,設(shè)計B 柱擾流器,改變車窗寬度、面積、厚度和車內(nèi)空腔體積等降噪措施可改善風(fēng)振噪聲。

        (3)B 柱擾流器對偏移距離與對聲壓級峰值的影響呈正相關(guān),而對峰值的頻率基本無影響。方案A 和B 以及其他能減小偏移距離的擾流器結(jié)構(gòu),可以減小湍流渦的強度,從而改善風(fēng)振噪聲。方案B-2的風(fēng)振降噪效果最好,可達12.8 dB。

        (4)道路試驗表明方案A 和B 能實際起到降低風(fēng)振噪聲的效果;且試驗的降噪效果與仿真的降噪效果一致性較好。

        本文研究結(jié)論對于進一步理解和控制風(fēng)振噪聲具有一定理論意義和工程應(yīng)用價值。

        猜你喜歡
        風(fēng)振聲壓級空腔
        機器噪聲平均聲壓級計算方法差異性實證研究
        電動工具(2024年1期)2024-02-29 01:40:24
        基于邊光滑有限元法的二維復(fù)合彈性空腔聲振特性分析
        一種計算消聲室聲壓級的新方法
        全新DXR mkll有源揚聲器
        演藝科技(2019年4期)2019-03-30 03:21:46
        張力對接觸網(wǎng)系統(tǒng)風(fēng)振響應(yīng)影響的風(fēng)洞試驗研究
        高塔設(shè)備風(fēng)振失效原因分析及改善措施
        空腔參數(shù)對重力壩穩(wěn)定的影響分析
        前置污水去油池
        前置污水去油池
        大跨屋蓋結(jié)構(gòu)MTMD風(fēng)振控制最優(yōu)性能研究
        亚洲国产成人久久精品一区| 色欲国产精品一区成人精品| 国产精品亚洲综合天堂夜夜 | 亚洲综合一区二区三区在线观看 | 一区二区精品天堂亚洲av| 蜜臀av在线观看| 中国xxx农村性视频| 亚洲熟妇大图综合色区| 国产精品美女一区二区av| 亚洲va中文字幕无码一二三区| 国产真实夫妇交换视频| 国产精品国产三级国产专播| 日本女u久久精品视频| 大又大又粗又硬又爽少妇毛片| 黑人巨茎大战欧美白妇| 91美女片黄在线观看| 国产自拍成人在线免费视频| av无码一区二区三区| 亚洲av无码一区二区三区在线| 亚洲无码激情视频在线观看| 国产不卡在线观看视频| 中文字幕乱码亚洲精品一区| 久久国产偷| 在线观看国产av一区二区 | 97久久草草超级碰碰碰| 人妻丰满熟妇av无码区hd| 亚洲精品女优中文字幕| 人妻少妇-嫩草影院| 国产无遮挡又黄又爽又色| 久久aⅴ无码av高潮AV喷| 日韩高清不卡一区二区三区| 欧美俄罗斯40老熟妇| 久久亚洲AV无码精品色午夜| 国产精品久久三级精品| 韩国三级在线观看久| 久久精品国产9久久综合| 亚洲中文字幕第二十三页| 亚洲啪啪视频一区二区| 特黄a级毛片免费视频| 99在线无码精品秘 入口九色 | 天天弄天天模|