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

        ?

        阻尼環(huán)車輪減振降噪特性分析

        2021-01-29 03:03:14賈尚帥趙新利張文敏徐涆文袁杰韓健
        機(jī)械 2020年12期
        關(guān)鍵詞:模態(tài)振動(dòng)模型

        賈尚帥,趙新利,張文敏,徐涆文,袁杰,韓健

        阻尼環(huán)車輪減振降噪特性分析

        賈尚帥1,趙新利1,張文敏1,徐涆文*,2,袁杰3,韓健2

        (1.中車唐山機(jī)車車輛有限公司 技術(shù)研究中心,河北 唐山 064000;2.西南交通大學(xué) 牽引動(dòng)力國(guó)家重點(diǎn)實(shí)驗(yàn)室,四川 成都 610031;3.西南交通大學(xué) 機(jī)械工程學(xué)院,四川 成都 610031)

        根據(jù)ISO 3745-2003標(biāo)準(zhǔn),在半消聲室內(nèi)采用力錘敲擊和落球撞擊的激勵(lì)方式,對(duì)兩種安裝阻尼環(huán)的車輪進(jìn)行模態(tài)阻尼比和輻射聲功率的測(cè)試和對(duì)比評(píng)價(jià),基于實(shí)驗(yàn)測(cè)試數(shù)據(jù)建立車輪有限元模型,研究在采用阻尼環(huán)車輪降噪措施以及多降噪措施組合后的輪軌振動(dòng)噪聲降低情況。發(fā)現(xiàn)車輪在安裝阻尼環(huán)后能夠有效提升1000~4500 Hz頻率范圍內(nèi)的車輪模態(tài)阻尼比,從而有效降低高頻的輪軌輻射噪聲,特別在3150 Hz位置采用阻尼環(huán)車輪措施后輪軌噪聲降低了13.3 dBA,但輪軌總輻射噪聲降低0.1 dBA。采用鋼軌TMD降噪措施后輪軌噪聲在630~2500 Hz范圍得到有效衰減,和阻尼環(huán)車輪降噪措施組合后阻尼環(huán)車輪的對(duì)輪軌總噪聲的降噪量增加為0.3 dBA。結(jié)果表明,輪軌輻射噪聲的顯著頻段為630~2000 Hz且由鋼軌占主導(dǎo),阻尼環(huán)車輪主要降低高頻位置的輪軌輻射噪聲,因此組合軌道降噪措施后能夠有效提升阻尼環(huán)車輪輪軌總輻射噪聲的降低效果。

        阻尼環(huán)車輪;減振降噪;模態(tài)阻尼比

        城市軌道交通中,車輪和鋼軌表面的不平順會(huì)導(dǎo)致輪軌振動(dòng)噪聲,進(jìn)而影響車輛運(yùn)行平穩(wěn)性和舒適性[1],且車外噪聲對(duì)人們的干擾主要也是來源于輪軌滾動(dòng)噪聲。車輪作為輪軌滾動(dòng)噪聲的主要輻射源,對(duì)高頻2000 Hz以上的輪軌滾動(dòng)噪聲起著主導(dǎo)作用,這是由于在高頻車輪的模態(tài)阻尼較低導(dǎo)致車輪振動(dòng)劇烈從而輻射出更多噪聲。因此增加車輪的阻尼可有效衰減車輪的振動(dòng)能量。在車輪上加入阻尼環(huán)結(jié)構(gòu)即為一種增加車輪模態(tài)阻尼比的方式。

        為進(jìn)一步認(rèn)識(shí)和了解阻尼環(huán)車輪的減振降噪性能,國(guó)內(nèi)外研究學(xué)者通過理論和實(shí)驗(yàn)對(duì)阻尼環(huán)車輪結(jié)構(gòu)進(jìn)行了詳細(xì)研究。J. F. Brunel通過仿真和試驗(yàn)方法分析了阻尼環(huán)車輪的降噪原理,結(jié)果表明,車輪振動(dòng)能量的衰減主要是由于阻尼環(huán)和車輪之間耦合振動(dòng)摩擦效應(yīng)[2-3]。為研究阻尼環(huán)預(yù)緊力、加裝阻尼環(huán)的數(shù)目以及阻尼環(huán)的直徑對(duì)車輪振動(dòng)聲輻射的影響,劉玉霞在實(shí)驗(yàn)中對(duì)上述阻尼車輪結(jié)構(gòu)進(jìn)行了測(cè)試[4],但其研究的阻尼環(huán)數(shù)目較少,具有一定的局限性。劉謀凱[5]詳細(xì)研究了在加裝五個(gè)阻尼環(huán)后車輪振動(dòng)噪聲的降低情況,并分析了在車輪開孔后其固有頻率的變化情況。目前對(duì)于阻尼環(huán)車輪僅局限在研究其自身的降噪情況,對(duì)阻尼環(huán)車輪對(duì)輪軌滾動(dòng)噪聲的降低情況分析還較為欠缺。

        本文通過實(shí)驗(yàn)方法研究了單雙兩種阻尼環(huán)車輪在懸掛狀態(tài)下的振動(dòng)聲輻射特性,并根據(jù)測(cè)試的車輪模態(tài)阻尼比數(shù)據(jù)和實(shí)際車輪結(jié)構(gòu)尺寸建立了阻尼環(huán)車輪振動(dòng)聲輻射模型,利用模型分析了在采用阻尼環(huán)車輪措施后輪軌滾動(dòng)噪聲的降低情況。

        1 阻尼車輪測(cè)試

        1.1 車輪頻響測(cè)試

        將車輪通過彈性繩進(jìn)行懸掛以模擬車輪的自由狀態(tài),將加速度傳感器布置在車輪踏面、輪輞和腹板位置,用力錘對(duì)車輪的名義滾動(dòng)圓位置進(jìn)行徑向敲擊,從而獲取車輪的頻響數(shù)據(jù)?,F(xiàn)場(chǎng)測(cè)試圖以及測(cè)點(diǎn)布置如圖1所示。

        圖1 阻尼環(huán)車輪及加速度測(cè)點(diǎn)位置

        1.2 車輪聲輻射測(cè)試

        本文在半消聲室內(nèi)對(duì)阻尼環(huán)車輪進(jìn)行聲學(xué)特性測(cè)試,圖2為車輪的現(xiàn)場(chǎng)測(cè)試圖片,車輪同樣采用彈性繩進(jìn)行懸掛,用落球徑向沖擊車輪,鋼球直徑為2.5 cm,為了保證落球的滾動(dòng)平穩(wěn)性,自動(dòng)落球裝置采用“V”型的凹槽結(jié)構(gòu)。振動(dòng)測(cè)點(diǎn)布置同圖1,根據(jù)聲學(xué)測(cè)點(diǎn)布置標(biāo)準(zhǔn)ISO 3745-2003[6],在半球型陣列架上布置20個(gè)傳聲器,通過所有傳聲器的聲壓級(jí),可以獲得車輪整體的聲能量級(jí)[5]。

        圖2 阻尼環(huán)車輪聲學(xué)特性測(cè)試現(xiàn)場(chǎng)圖

        2 車輪模型

        根據(jù)城市軌道列車車輪的實(shí)際結(jié)構(gòu)尺寸建立車輪的有限元模型[7],由于加裝阻尼環(huán)后車輪的振動(dòng)模態(tài)變化較小[5],本文采用無阻尼環(huán)結(jié)構(gòu)的車輪實(shí)體有限元模型,計(jì)算車輪的固有頻率和振型,并導(dǎo)入實(shí)測(cè)模態(tài)阻尼比以模擬阻尼環(huán)的阻尼耗能作用,修正仿真模型,計(jì)算得到車輪表面的振動(dòng)特性。

        3 結(jié)果分析

        車輪的模態(tài)振型主要為面外的軸向振動(dòng)模態(tài)以及面內(nèi)的、均通過節(jié)徑數(shù)表征的徑向振動(dòng)模態(tài)和周向振動(dòng)模態(tài),可以依次表示為(,)、(,)、(,);其中,為節(jié)圓數(shù),為節(jié)徑數(shù),為徑向,為周向。

        如表1所示,基于有限元法計(jì)算和頻率響應(yīng)函數(shù)測(cè)試,給出了車輪軸向、徑向模態(tài)的振型、仿真計(jì)算和實(shí)測(cè)固有頻率,由于未對(duì)車輪的周向模態(tài)進(jìn)行測(cè)試,僅給出了車輪周向模態(tài)的振型、仿真計(jì)算固有頻率。由數(shù)據(jù)可得,車輪有限元模型能夠較好地反映出車輪的固有模態(tài)頻率。仿真結(jié)果同兩種阻尼環(huán)車輪的測(cè)試結(jié)果最大相差值小于5%。結(jié)果表明車輪有限元模型能夠較為準(zhǔn)確地反映阻尼環(huán)車輪的振動(dòng)特性,驗(yàn)證了模型的準(zhǔn)確性。

        圖3~圖5分別為車輪軸向、徑向和周向的模態(tài)振型。

        表1 車輪軸向、徑向、周向模態(tài)頻率和振型

        圖6為加裝阻尼環(huán)后車輪各部分頻響函數(shù)的測(cè)試結(jié)果,可知,加裝阻尼環(huán)后,車輪在踏面和腹板位置的振動(dòng)峰值均減小,同時(shí)加裝兩個(gè)阻尼環(huán)后車輪的振動(dòng)峰值較單一阻尼環(huán)車輪振動(dòng)峰值衰減效果更好,表明增加車輪上阻尼環(huán)數(shù)量能一定程度提升車輪的減振性能。

        由圖7可得,標(biāo)準(zhǔn)車輪模態(tài)阻尼比在全頻段都比較低,裝配了金屬環(huán)后,車輪的模態(tài)阻尼比在1000~4500 Hz各頻段均有些許提高,但在4500 Hz以上提高不明顯;雙側(cè)都裝配了金屬環(huán)后,車輪的模態(tài)阻尼比在單阻尼環(huán)車輪的基礎(chǔ)上又有明顯提高,尤其在1000~2000 Hz提高十分顯著。

        圖3 車輪軸向模態(tài)

        圖4 車輪徑向模態(tài)

        自由狀態(tài)下,金屬環(huán)阻尼車輪在較寬的頻段內(nèi)有較好的降噪效果,具體如表2所示。

        圖5 車輪周向模態(tài)

        圖6 頻響函數(shù)

        表2 金屬環(huán)阻尼車輪降噪效果

        4 輪軌滾動(dòng)噪聲模型

        4.1 軌道模型

        根據(jù)雙塊式無砟軌道的實(shí)際尺寸模型和結(jié)構(gòu)參數(shù)對(duì)其進(jìn)行有限元建模,如圖8所示,相關(guān)結(jié)構(gòu)參數(shù)如表3所示。

        圖7 五種車輪阻尼比測(cè)試結(jié)果比較

        圖8 雙塊式無砟軌道有限元模型

        表3 軌道模型參數(shù)

        模型中鋼軌、軌枕和道床采用實(shí)體單元進(jìn)行離散,扣件系統(tǒng)采用線性的彈簧阻尼單元進(jìn)行模擬。單位外激勵(lì)力作用于鋼軌跨中軌頭正上方,軌道板下部進(jìn)行全約束[7]。

        4.2 輪軌相互作用模型

        4.2.1 接觸濾波

        實(shí)際中,輪軌之間的接觸是一個(gè)橢圓形的接觸斑,當(dāng)粗糙度的波長(zhǎng)尺寸小于或等于接觸斑的尺寸時(shí),輪軌表面的粗糙度就不會(huì)對(duì)輪軌的相互作用產(chǎn)生影響。接觸濾波的函數(shù)為:

        式中:為鋼軌車輪表面粗糙度的相關(guān)系數(shù);為粗糙度波數(shù)量;為接觸圓尺寸,mm;1()為貝塞爾函數(shù)。

        4.2.2 接觸剛度

        鋼軌和車輪之間的接觸剛度為[9-10]:

        式中:為接觸結(jié)構(gòu)的表面常數(shù),通過文獻(xiàn)[11]求得;為鋼軌和車輪的彈性模量,Pa;為泊松比;0為單個(gè)車輪的靜載荷,kN;R為車輪的半徑,mm;R為鋼軌軌頭曲率半徑,mm。

        輪軌力為:

        式中:為車輪軌道聯(lián)合粗糙度,mm;α為車輪柔度,N/m;α為鋼軌柔度,N/m;α為接觸剛度倒數(shù),N/m。

        4.3 輪軌聲輻射模型

        在求得車輪軌道的振動(dòng)響應(yīng)后,基于聲學(xué)邊界元理論提取軌道以及車輪表面網(wǎng)格作為聲學(xué)計(jì)算的邊界元網(wǎng)格模型,從而求得車輪和軌道的輻射聲功率,從而得到輪軌總輻射噪聲。為了保證聲學(xué)分析的準(zhǔn)確性,聲學(xué)邊界元模型的網(wǎng)格尺寸需要小于最短計(jì)算波長(zhǎng)的1/6[12]。

        5 阻尼環(huán)車輪降噪特性分析

        由測(cè)試結(jié)果可得,安裝阻尼環(huán)后,車輪頻響函數(shù)的高頻振動(dòng)峰值得到有效削弱,同時(shí),在整個(gè)頻段內(nèi)車輪聲能量得到有效衰減。結(jié)果表明,阻尼環(huán)能夠有效降低車輪的振動(dòng)和聲輻射。然而實(shí)際中,采取阻尼環(huán)車輪降噪措施后,輪軌滾動(dòng)噪聲的降低情況較為有限,本章對(duì)采取阻尼環(huán)車輪后的輪軌滾動(dòng)噪聲降低情況進(jìn)行分析。

        5.1 阻尼環(huán)車輪降噪分析

        基于建立的輪軌噪聲預(yù)測(cè)模型,本節(jié)分析采用阻尼環(huán)車輪措施后輪軌滾動(dòng)噪聲的降低情況,列車運(yùn)行速度為80 km/h,車輪和鋼軌的表面粗糙度等級(jí)為C。

        由圖9可知,阻尼環(huán)車輪降噪措施主要對(duì)高頻輪軌總輻射噪聲起作用,采用阻尼環(huán)車輪后,在頻率3150 Hz處輪軌總輻射噪聲為83.4 dBA,較采用普通標(biāo)準(zhǔn)車輪的輪軌總輻射噪聲降低13.3 dBA。但在整個(gè)頻段,采用阻尼環(huán)車輪降噪措施后輪軌總輻射噪聲僅降低0.1 dBA。

        圖9 標(biāo)準(zhǔn)車輪與阻尼環(huán)車輪輪軌總輻射噪聲對(duì)比

        由上述分析可得,在車輪上加裝阻尼環(huán)對(duì)輪軌總輻射噪聲的降低能力較為有限。為了分析其主要原因,圖10對(duì)車輪、鋼軌和輪軌總輻射噪聲進(jìn)行了分析,車輪對(duì)于輪軌總輻射噪聲的貢獻(xiàn)主要在2500 Hz以上的高頻輪軌噪聲,而2500 Hz以下較寬的頻帶輪軌總輻射噪聲均由鋼軌的輻射噪聲占主導(dǎo)。因此,車輪噪聲的降低能夠有效降低輪軌高頻的噪聲,但在整個(gè)頻段內(nèi)的噪聲控制能力還較為有限。

        圖10 輪軌各部分噪聲貢獻(xiàn)量

        5.2 多措施組合對(duì)輪軌噪聲降低情況分析

        鋼軌TMD(Tune Mass Damping,調(diào)諧質(zhì)量阻尼器)結(jié)構(gòu)是一種鋼軌的減振降噪措施,主要是通過調(diào)節(jié)其固有頻率同主振動(dòng)結(jié)構(gòu)發(fā)生諧振,從而降低主振動(dòng)結(jié)構(gòu)的振動(dòng),以此來抑制輪軌總輻射噪聲。為了分析在采用鋼軌TMD及車輪阻尼環(huán)兩種降噪措施組合情況下輪軌總輻射輻射噪聲的降低情況,本文在上述的輪軌噪聲預(yù)測(cè)模型中加入鋼軌TMD結(jié)構(gòu)模型,結(jié)構(gòu)的質(zhì)量塊模型由實(shí)體單元模擬,橡膠結(jié)構(gòu)則通過彈簧阻尼單元來模擬,如圖11、表4所示。

        圖11 鋼軌TMD結(jié)構(gòu)

        表4 鋼軌TMD模型參數(shù)

        由圖12可得,加裝鋼軌TMD結(jié)構(gòu)后,輪軌總輻射噪聲降低2.4 dBA,且輪軌總輻射噪聲在630~2500 Hz頻段內(nèi)得到了有效緩解。根據(jù)上述分析可知,此頻段的輪軌總輻射噪聲由鋼軌占主導(dǎo),鋼軌TMD結(jié)構(gòu)能夠有效降低鋼軌在630~2500 Hz范圍內(nèi)的振動(dòng)和噪聲輻射,從而對(duì)輪軌總輻射噪聲的影響較大。

        圖12 采用鋼軌TMD結(jié)構(gòu)措施輪軌總輻射噪聲

        由圖13可得,加裝鋼軌TMD結(jié)構(gòu)后,采用阻尼環(huán)車輪對(duì)輪軌總輻射噪聲的降低效果得到有效提升,從112.8 dBA降低到了112.5 dBA。

        圖13 雙重降噪措施組合后輪軌總噪聲降低情況

        上述研究表明,阻尼環(huán)車輪降噪措施能有效降低輪軌高頻位置的輻射噪聲,但由于輪軌輻射噪聲的主要頻段在400~2500 Hz的中高頻段,而此頻段的輪軌輻射噪聲往往是由軌道占主導(dǎo)地位,在采用鋼軌TMD結(jié)構(gòu)后降低了此頻段的鋼軌輻射噪聲從而使得輪軌在這個(gè)頻段的噪聲也得到有效降低,那么車輪部分的貢獻(xiàn)就開始變得突出,此時(shí)采用車輪阻尼環(huán)降噪措施就能有效提升其降噪效果。

        6 結(jié)論

        本文根據(jù)實(shí)驗(yàn)測(cè)試的阻尼環(huán)車輪模態(tài)阻尼比建立了阻尼環(huán)車輪有限元模型,基于TWINS建立了輪軌噪聲預(yù)測(cè)模型,利用模型分析了在采用阻尼環(huán)車輪后輪軌噪聲的降低情況,得出以下結(jié)論:

        (1)在實(shí)驗(yàn)室測(cè)試中,車輪加裝阻尼環(huán)結(jié)構(gòu)能夠有效提升車輪的阻尼,使得車輪在高頻位置的峰值得到有效衰減;

        (2)對(duì)于輪軌滾動(dòng)噪聲,阻尼環(huán)車輪措施能夠有效降低2500 Hz以上頻段的高頻噪聲。但對(duì)于輪軌總噪聲輻射降低效果較為有限;

        (3)在采用鋼軌TMD結(jié)構(gòu)后阻尼環(huán)車輪的降噪效果得到了一定的提升,表明在中高頻段輪軌振動(dòng)噪聲由軌道占主導(dǎo)地位,車輪對(duì)輪軌總噪聲的影響主要在高頻段。

        [1]張國(guó)平,張志波,吳興文,等. 多邊形激勵(lì)下動(dòng)車組動(dòng)態(tài)響應(yīng)研究[J]. 機(jī)械,2020,5(42):67-74.

        [2]BRUNEL J F. Modelling of squeal noise attenuation ofring damped wheels[J]. Applied Acoustics,2004(65):457-471

        [3]BRUNEL J F,DUFRéNOY,CHARLEY J. Analysis of theattenuation of railway squeal noise by preloaded ringsinserted in wheels[J]. Journal of the Acoustical SocietyofAmerica,2010,127(3):1300-1306.

        [4]劉玉霞. 阻尼環(huán)車輪振動(dòng)聲輻射特性試驗(yàn)研究[D]. 成都:西南交通大學(xué),2016.

        [5]劉謀凱,周信,韓健,等. 迷宮式阻尼環(huán)裝置對(duì)車輪的減振降噪效果[J]. 噪聲與振動(dòng)控制,2018,38(4):221-225.

        [6]ISO 3745-2003,Acoustics-determination of sound powerlevels of noise sources using sound pressure-Precisionmethods for anechoic and hemi-anechoic rooms,International Standard [S]. International Organization forStandardization,Geneva,Switzerland,2012.

        [7]韓健,王瑞乾,王諦,等. 徑向激勵(lì)下輻板型式對(duì)車輪振動(dòng)聲輻射影響仿真和試驗(yàn)研究[J]. 機(jī)械工程學(xué)報(bào),2014,50(10):103-111.

        [8]徐涆文,韓健,劉曉龍,等. 彈性軌枕軌道振動(dòng)特性研究[J].噪聲與振動(dòng)控制,2019,39(2):156-161.

        [9]THOMPSON D J,F(xiàn)ODIMAN P,MAHE′ H. Experimental validation of the TWINS prediction program, Part 2:results[J]. Journal of Sound and Vibration,1996(193):137-147

        [10]D J Thompson. Wheel-rail noise generation, Part I: introduction and interaction model[J]. Journal of Sound and Vibration,1993,161(3):387-400.

        [11]Thompson D. Railway Noise and Vibration[M]. Oxford, UK:Elsevier’s Science & Technology Rights Department,2009:6-10.

        [12]韓健. 地鐵列車-嵌入式軌道系統(tǒng)動(dòng)力學(xué)性能及振動(dòng)噪聲特性研究[D]. 成都:西南交通大學(xué),2018.

        Analysis of Vibration and Noise Reduction Characteristics of Damping Ring Wheels

        JIA Shangshuai1,ZHAO Xinli1,ZHANG Wenmin1,XU Hanwen2,YUAN Jie3,HAN Jian2

        ( 1.Technology Research Center, CRRC TANGSHAN Co., Ltd., Tangshan 064000, China;2.State Key Laboratory of Traction Power, Southwest Jiaotong University,Chengdu 610031, China;3.School of Mechanical Engineering, Southwest Jiaotong University, Chengdu 610031, China )

        In order to be able to have a better understanding of the wheel/rail noise decrease after using themetal damping ring wheel, in this paper, according to ISO 3745-2012 standard, in the semi-dead room, using force hammer percussion and falling ball impact of incentives, the wheels of installation of the two damping rings modal damping ratio and the radiation sound power test and contrast evaluation, wheel finite element model was built based on the finite element method, and use the wheel to wheel/rail rolling noise reduction is analyzed.The structure shows that the modal damping ratio of the wheel is effectively improved after the installation of damping ring. As for wheel-rail rolling noise, metal damping ring wheel mainly reduces the high-frequency wheel-rail noise radiation above 2500 Hz.

        metal damping ring wheel;vibration and noise reduction;modal damping ratio

        U270.1+6

        A

        10.3969/j.issn.1006-0316.2020.12.002

        1006-0316 (2020) 12-0010-07

        2020-07-03

        國(guó)家重點(diǎn)研發(fā)計(jì)劃資助項(xiàng)目(2016YFB1200503-02)

        賈尚帥(1982-),男,河北唐山人,博士,高級(jí)工程師,主要研究方向?yàn)檠邪l(fā)設(shè)計(jì)、振動(dòng)噪聲。*通信作者:徐涆文(1993-),男,四川成都人,博士研究生,主要研究方向?yàn)楦咚倭熊嚋p振降噪、舒適性研究,E-mail:xheaven9@163.com。

        猜你喜歡
        模態(tài)振動(dòng)模型
        一半模型
        振動(dòng)的思考
        重要模型『一線三等角』
        振動(dòng)與頻率
        重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
        中立型Emden-Fowler微分方程的振動(dòng)性
        3D打印中的模型分割與打包
        國(guó)內(nèi)多模態(tài)教學(xué)研究回顧與展望
        基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識(shí)別
        UF6振動(dòng)激發(fā)態(tài)分子的振動(dòng)-振動(dòng)馳豫
        风韵丰满妇啪啪区老老熟女杏吧| 亚洲成av人片在www鸭子| 亚洲综合一区二区三区天美传媒| 亚洲av精二区三区日韩| 日本中文字幕一区二区高清在线 | 久久狼人国产综合精品| 国产精品人妻熟女男人的天堂| 中文字幕国产精品一二三四五区 | 国产白浆精品一区二区三区| 久久久噜噜噜久久熟女| 精品国产一区二区三区av免费| 人妻少妇乱子伦无码视频专区| 午夜性无码专区| 久久精品国产99精品国偷| 少妇勾引视频网站在线观看| av剧情演绎福利对白| 丰满少妇被粗大猛烈进人高清| 国产精品一区二区久久乐下载| 久久久久亚洲av无码专区体验| 免费国产h视频在线观看86| 人妻中文字幕在线一二区| 日本动漫瀑乳h动漫啪啪免费 | 亚洲欧美激情在线一区| 欧美亚洲国产精品久久高清 | 亚洲av无码成人yellow| 国产小视频一区二区三区| 日本女优久久精品观看| 亚洲av中文无码乱人伦下载| 精品性高朝久久久久久久| 蜜桃一区二区三区自拍视频| 国产91色综合久久高清| 少妇人妻中文字幕hd| 国产乱子伦精品免费无码专区 | 亚洲熟女网站| 色偷偷亚洲精品一区二区 | 国产成人综合日韩精品无码| 比比资源先锋影音网| 91最新免费观看在线| 亚洲综合久久久中文字幕| 插入日本少妇一区二区三区| 夜夜揉揉日日人人青青|