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

        ?

        制冷劑管道參數(shù)對(duì)冰箱振動(dòng)影響特性

        2016-01-15 06:12:20孔祥強(qiáng),陳麗娟,李瑛
        振動(dòng)與沖擊 2015年20期
        關(guān)鍵詞:模態(tài)分析固有頻率冰箱

        第一作者孔祥強(qiáng)男,博士,副教授,1976年生

        郵箱:kxqiangly@126.com

        制冷劑管道參數(shù)對(duì)冰箱振動(dòng)影響特性

        孔祥強(qiáng),陳麗娟,李瑛

        (山東科技大學(xué)機(jī)械電子工程學(xué)院,山東青島266590)

        摘要:為減小冰箱振動(dòng)、降低噪聲,利用Pro/Engineer軟件分別建立冰箱制冷劑管道內(nèi)部氣柱及管道三維實(shí)體模型,用有限元分析軟件ANSYS分析模態(tài),獲得氣柱及管道固有頻率。氣柱固有頻率與壓縮機(jī)激發(fā)頻率相差較大不會(huì)產(chǎn)生共振,而管道固有頻率在共振區(qū)內(nèi)會(huì)產(chǎn)生共振。在制冷劑管道材料不變條件下分析不同管道長度、彎管圓角半徑及固定支撐位置對(duì)冰箱振動(dòng)特性影響。結(jié)果表明,減小制冷劑管道長度會(huì)增大固有頻率并逐漸遠(yuǎn)離共振區(qū),從而避免管道共振;增加制冷劑管道彎管圓角半徑,其固有頻率略有增大,一階固有頻率仍在共振區(qū)內(nèi)導(dǎo)致管道共振;在制冷劑管道不同位置增加固定支撐,能有效提高管道固有頻率,使其遠(yuǎn)離共振區(qū),且存在最佳安裝位置。

        關(guān)鍵詞:冰箱;制冷劑管道;振動(dòng);模態(tài)分析;固有頻率

        基金項(xiàng)目:山東省高等學(xué)校科技計(jì)劃資助項(xiàng)目(J11LD63);山東省泰山學(xué)者建設(shè)工程專項(xiàng)經(jīng)費(fèi)資助項(xiàng)目;青島市黃島區(qū)科技項(xiàng)目(2014-1-40)

        收稿日期:2014-08-01修改稿收到日期:2014-09-30

        中圖分類號(hào):TB535

        文獻(xiàn)標(biāo)志碼:A

        DOI:10.13465/j.cnki.jvs.2015.20.035

        Abstract:In order to reduce the vibration and noise of a refrigerator, geometric models of the refrigerant pipe and the gas column in the pipe were created respectively by using Pro/Engineer software, and natural frequencies of the pipe and the gas column were obtained by the modal analysis based on ANSYS software. It is found that gas column resonance will not appear due to the large difference between compressor excitation frequency and the gas column natural frequency, but pipe resonance will occur when the pipe natural frequency is in the resonance zone. The vibration characteristics of the pipes mode of the same material yet with different length, fillet radius and support location were analyzed. The model analysis results show that, the pipe natural frequency increases continuously and gets far away gradually from the resonance zone with the decrease of refrigerant pipe length. The pipe natural frequency is slightly getting higher with the increase of fillet radius, however the pipe resonance will still occur with the first order natural frequency in the resonance zone. The pipe natural frequency can be improved effectively and get far away from the resonance zone by adding supports in different locations, and a best location exists.

        Influence of refrigerant pipe parameters on vibration characteristics of refrigerator

        KONGXiang-qiang,CHENLi-juan,LIYing(College of Mechanical and Electronic Engineering, Shandong University of Science and Technology, Qingdao 266590, China)

        Key words:refrigerator; refrigerant pipe; vibration; modal analysis; natural frequency

        電冰箱作為廣大家庭不可或缺的電器運(yùn)行時(shí)噪聲大小頗受關(guān)注,并與節(jié)能、環(huán)保同為評(píng)定家電產(chǎn)品質(zhì)量?jī)?yōu)劣的三大要素。因而國標(biāo)[1]對(duì)冰箱、空調(diào)、洗衣機(jī)、微波爐、吸油煙機(jī)及電風(fēng)扇等6類家電噪聲限值進(jìn)行明確規(guī)定。其中,冰箱容積<250 L的噪聲值不超過45 dB,容積>250 L的噪聲值不超過48 dB。

        對(duì)冰箱降噪已有諸多研究。如Gue等[2]通過對(duì)冰箱風(fēng)扇噪聲試驗(yàn)研究及數(shù)值模擬提出低噪聲風(fēng)扇設(shè)計(jì)方案。Seo等[3]利用穿孔板結(jié)構(gòu)降低冰箱風(fēng)扇噪聲。褚志剛等[4]設(shè)計(jì)出冰箱風(fēng)扇振動(dòng)質(zhì)量在線檢測(cè)系統(tǒng)。陳建良等[5]對(duì)冰箱壓縮機(jī)殼體輻射噪聲進(jìn)行數(shù)值分析。仲崇明等[6]對(duì)冰箱往復(fù)式壓縮機(jī)振動(dòng)進(jìn)行有限元數(shù)值模擬,驗(yàn)證其用于振動(dòng)響應(yīng)分析的可行性。郭維等[7]用錘擊脈沖激勵(lì)法對(duì)冰箱壓縮機(jī)殼進(jìn)行試驗(yàn)?zāi)B(tài)分析,獲得整機(jī)殼體振動(dòng)特性。趙科等[8]對(duì)冰箱用動(dòng)磁式直線壓縮機(jī)的動(dòng)態(tài)特性進(jìn)行仿真。

        以上研究主要集中于冰箱風(fēng)扇及壓縮機(jī)降噪,而冰箱內(nèi)部制冷劑管道振動(dòng)對(duì)冰箱噪聲亦有較大影響。管道振動(dòng)主要由內(nèi)部氣流脈動(dòng)引起,壓縮機(jī)為冰箱內(nèi)部制冷劑管道振動(dòng)激發(fā)源。Kim等[9]通過測(cè)量冰箱壓縮機(jī)排氣管試驗(yàn)?zāi)P偷臍饬髅}動(dòng)認(rèn)為,氣流脈動(dòng)與噪聲間存在一定關(guān)系。Lee等[10]試驗(yàn)研究冰箱壓縮機(jī)管道的振動(dòng)特性。劉益才等[11]通過研究冰箱毛細(xì)管出口處噴射噪聲,獲得湍流強(qiáng)度、壓力及流速沿流動(dòng)方向變化及流場(chǎng)分布特性。Han等[12]通過改變蒸發(fā)器入口管形狀及布局,使冰箱輻射噪聲在315~3 150 Hz頻率范圍內(nèi)減少約2~5 dB。李春銀等[13]用有限元分析軟件對(duì)空調(diào)旋葉式壓縮機(jī)排氣閥片進(jìn)行模態(tài)分析,獲得排氣閥片振動(dòng)特性。

        用有限元分析軟件對(duì)冰箱制冷劑管道進(jìn)行模態(tài)分析研究較少。為此,本文基于多自由度系統(tǒng)振動(dòng)微分方程,用有限元軟件ANSYS對(duì)冰箱局部制冷劑管道進(jìn)行模態(tài)分析,計(jì)算管道系統(tǒng)固有頻率,獲得管道結(jié)構(gòu)參數(shù)對(duì)冰箱振動(dòng)的影響特性。

        1制冷劑管道系統(tǒng)數(shù)學(xué)模型

        氣流脈動(dòng)引起的管道振動(dòng)非常復(fù)雜,且存在兩振動(dòng)系統(tǒng),即氣柱及管道機(jī)械振動(dòng)系統(tǒng)[14]。制冷劑管道內(nèi)部充滿氣體稱為氣柱,由于氣體可壓縮膨脹,可將氣柱視為有連續(xù)質(zhì)量的彈簧振動(dòng)系統(tǒng),具有固有頻率及振型,氣柱受壓縮機(jī)周期性吸排氣激發(fā)會(huì)產(chǎn)生振動(dòng);而管道本身也具有固有頻率,只要激發(fā)源作用于管道,則會(huì)產(chǎn)生機(jī)械振動(dòng)響應(yīng)。此兩種均為多自由度振動(dòng)系統(tǒng),其振動(dòng)微分方程為一組互相耦合的常微分方程組,直接求解、分析較困難[15]。為此,本文采用模態(tài)分析方法分析求解。

        1.1管道內(nèi)部氣柱固有頻率計(jì)算數(shù)學(xué)模型

        管道內(nèi)部氣柱固有頻率常采用轉(zhuǎn)移矩陣法計(jì)算,但轉(zhuǎn)移矩陣法建立在一維數(shù)學(xué)模型基礎(chǔ)上,用于計(jì)算管道內(nèi)部氣柱固有頻率較繁瑣且精度不高。本文用模態(tài)分析方法計(jì)算冰箱制冷劑管道內(nèi)部氣柱固有頻率。

        設(shè)制冷劑氣體無粘性、可壓縮;氣體平均密度、壓力、速度均為定值;介質(zhì)中傳播為小振幅聲波。

        氣柱聲學(xué)波動(dòng)方程[16]為

        (1)

        式中:▽2為拉普拉斯算子;c為氣體聲速;P為氣體聲壓;t為時(shí)間。

        利用微分方程等效積分形式并結(jié)合加權(quán)余量法的伽遼金(Galerkin)法,對(duì)管道內(nèi)部氣柱振動(dòng)系統(tǒng)建立有限元方程為

        (2)

        式中:[M],[K]分別為氣柱質(zhì)量、剛度矩陣;{P}為壓力向量;[F]為管道內(nèi)部氣柱所受載荷矩陣。

        設(shè)式(2)中[F]=0,即管道氣柱固有頻率數(shù)學(xué)模型為

        (3)

        1.2管道固有頻率計(jì)算數(shù)學(xué)模型

        冰箱制冷劑管道系統(tǒng)振動(dòng)微分方程為

        (4)

        設(shè)式(4)中[C],[F]均為0,即為管道固有頻率數(shù)學(xué)模型

        (5)

        2壓縮機(jī)激發(fā)頻率

        圖1 冰箱內(nèi)部結(jié)構(gòu)示意圖 Fig.1 Schematic diagram of refrigeration system

        冰箱制冷系統(tǒng)內(nèi)部結(jié)構(gòu)見圖1。采用單級(jí)活塞式制冷壓縮機(jī),工作轉(zhuǎn)速3 000r/min;管式蒸發(fā)器,材料銅鋁合金,總長14.96 m,外徑8 mm,壁厚0.6 mm;翹片管式冷凝器,材料表面鍍鋅鋼管,長8.76 m,外徑9.52 mm,壁厚0.35 mm;制冷劑管道、冰箱殼體全部用螺釘固定;0點(diǎn)位置為排氣管原有固定支撐點(diǎn),1~5點(diǎn)分別為新設(shè)不同固定支撐添加位置。

        冰箱管道振動(dòng)激發(fā)源為壓縮機(jī),激發(fā)頻率算式為

        (6)

        式中:n為壓縮機(jī)轉(zhuǎn)速;i為壓縮機(jī)作用方式,對(duì)單級(jí)制冷壓縮機(jī)i=1;m為簡(jiǎn)諧階次。

        壓縮機(jī)激發(fā)頻率計(jì)算結(jié)果見表1。

        表1 壓縮機(jī)激發(fā)頻率

        3管道內(nèi)部氣柱固有頻率

        壓縮機(jī)周期性吸排氣會(huì)激發(fā)其相連管道振動(dòng),為獲得管道參數(shù)對(duì)冰箱振動(dòng)影響規(guī)律,以圖1壓縮機(jī)排氣口及固定支撐間制冷劑管道為分析對(duì)象。該管道長900 mm,彎管圓角半徑30 mm,管道內(nèi)部制冷劑熱物理參數(shù)設(shè)為:壓力1.8 MPa,溫度50 ℃,制冷劑氣體密度2.544 kg/m3,聲速340 m/s。利用Pro/E建立管道內(nèi)部氣柱三維模型,再導(dǎo)入ANSYS進(jìn)行模態(tài)分析。選ANSYS流體分析模塊,單元類型選三維聲單元,冰箱管道內(nèi)氣柱有限元模型見圖2。

        圖2 管道內(nèi)部氣柱有限元模型 Fig.2 Finite element model of gas column in pipe

        計(jì)算得管道內(nèi)氣柱固有頻率見表2。由表2看出,管道內(nèi)氣柱固有頻率較小,與壓縮機(jī)激發(fā)頻率差距較大。因此氣柱不會(huì)發(fā)生共振及由此產(chǎn)生強(qiáng)烈噪聲。

        表2 管道內(nèi)部氣柱固有頻率

        4管道固有頻率

        冰箱內(nèi)部制冷劑管道為連續(xù)彈性體結(jié)構(gòu),具有無窮多自由度及固有頻率。系統(tǒng)固有頻率作為結(jié)構(gòu)的固有屬性,其大小只與結(jié)構(gòu)質(zhì)量、剛度有關(guān),與外載荷無關(guān)[17]。如果冰箱管道結(jié)構(gòu)設(shè)計(jì)不合理,其固有頻率會(huì)接近壓縮機(jī)激發(fā)頻率基頻或倍頻,使壓縮機(jī)運(yùn)行頻率范圍內(nèi)出現(xiàn)管道系統(tǒng)固有模態(tài),導(dǎo)致管道產(chǎn)生共振。

        仍以壓縮機(jī)排氣口與固定支撐間局部制冷劑管道為分析對(duì)象,管道有限元模型見圖3。計(jì)算所得管道固有頻率見表3。

        表3 管道固有頻率

        圖3 管道有限元模型 Fig.3 Finite element model of pipe

        固有頻率為激發(fā)頻率的0.8~1.2倍(即共振區(qū))時(shí)管道會(huì)發(fā)生共振。管道固有頻率與共振區(qū)關(guān)系見圖4。由圖4看出,一階固有頻率落在共振區(qū)內(nèi),五階靠近共振區(qū),此時(shí)管道將發(fā)生共振,且以較大振幅強(qiáng)烈撞擊冰箱殼體,繼而產(chǎn)生較大噪聲。因此,本文通過改變管道長度、彎管處圓角半徑及添加固定支撐等進(jìn)行優(yōu)化設(shè)計(jì),能避免或防止共振發(fā)生,達(dá)到降噪聲目的。

        圖4 管道固有頻率與共振區(qū) Fig.4 Natural frequency of pipe and resonance region

        4.1管道長度對(duì)固有頻率影響

        管道長度分別取700 mm、800 mm、900 mm、1 000 mm及1 100 mm五組數(shù)據(jù),固有頻率計(jì)算結(jié)果見圖5。由圖5可見,隨管道長度增加固有頻率減??;長度由900 mm增到1 000 mm時(shí),其一階~五階固有頻率均落在共振區(qū);長度為1 100 mm時(shí)固有頻率全部落在共振區(qū);長度減小到800 mm或700 mm時(shí)固有頻率遠(yuǎn)離共振區(qū),且隨長度減小固有頻率增大,因而能有效避免共振發(fā)生。減小管道長度可使固有頻率與壓縮機(jī)激發(fā)頻率錯(cuò)開,既能避免共振又能降低管道振動(dòng)及噪聲。

        圖5 管道長度對(duì)固有頻率影響 Fig.5 Effect of length on pipe natural frequency

        4.2彎管圓角半徑對(duì)管道固有頻率影響

        管道其它結(jié)構(gòu)參數(shù)不變條件下,彎管圓角半徑分別取20 mm、30 mm、40 mm、50 mm及60 mm。計(jì)算所得管道固有頻率與共振區(qū)關(guān)系見圖6。由圖6可知,隨圓角半徑增加管道固有頻率略有增大,但增幅較小。圓角半徑增大到60 mm時(shí)管道一階固有頻率仍落在共振區(qū)。故改變彎管圓角半徑不能達(dá)到減振降噪目的。

        圖6 圓角半徑對(duì)管道固有頻率影響 Fig.6 Effect of round radius on pipe natural frequency

        4.3固定支撐對(duì)管道固有頻率影響

        落在共振區(qū)的管道一階固有頻率對(duì)應(yīng)振型見圖7。由圖7可見,管道發(fā)生共振時(shí)拐彎處振動(dòng)最強(qiáng)烈,振幅最大,故此處產(chǎn)生噪聲較大。增加固定支撐,改變?cè)艿兰s束條件,增強(qiáng)剛度,其固有頻率會(huì)隨之增大,可有效避免管道產(chǎn)生共振,不會(huì)產(chǎn)生強(qiáng)烈噪聲。

        圖7 管道一階固有頻率對(duì)應(yīng)振型 Fig.7 Vibration model at first natural frequency

        分別在圖1的1~5點(diǎn)處即管道拐點(diǎn)及離拐點(diǎn)40 mm、60 mm、80 mm、100 mm五處增加固定支撐,計(jì)算結(jié)果見圖8。由圖8看出,增加固定支撐可使管道固有頻率大幅度增加,遠(yuǎn)離共振區(qū)。值得注意的是,管道固有頻率并不隨固定支撐與拐點(diǎn)間距離增加呈線性變化。本文設(shè)定工況條件下,固定支撐與拐點(diǎn)間距離為60 mm時(shí)管道固有頻率增加最大,高于壓縮機(jī)激發(fā)頻率,制冷劑管道不會(huì)發(fā)生共振,亦能有效降低噪聲。

        圖8 固定支撐對(duì)管道固有頻率影響 Fig.8 Effect of support location on pipe natural frequency

        5結(jié)論

        利用ANSYS模態(tài)分析模塊計(jì)算制冷劑管道內(nèi)氣柱及管道固有頻率,分析長度、彎管圓角半徑及增加固定支撐對(duì)管道固有頻率影響,結(jié)論如下:

        (1)氣柱固有頻率與壓縮機(jī)激發(fā)頻率相差較大,不會(huì)發(fā)生共振,即不會(huì)產(chǎn)生強(qiáng)烈噪聲。

        (2)減小長度能使管道固有頻率遠(yuǎn)離共振區(qū),進(jìn)而降低管道的振動(dòng)、噪聲。

        (3)增加彎管圓角半徑對(duì)管道固有頻率增幅較小,在激發(fā)源作用下仍會(huì)發(fā)生共振,產(chǎn)生較大噪聲。

        (4)在合適位置增加固定支撐,可使管道固有頻率大幅度增加,能有效避免共振及降低噪聲。

        參考文獻(xiàn)

        [1]GB/T 19606-2004,家用和類似用途電器噪聲限值[S].

        [2]Gue F, Cheong C, Kim T. Development of low-noise axial cooling fans in a household refrigerator[J]. Journal of Mechanical Science and Technology, 2011, 25(12): 2995-3004.

        [3]Seo J Y, Kim W J, Won J S. Design and application of a perforated panel system to a household refrigerator for cooling fan noise reduction[J]. Journal of Mechanical Engineering Science, 2012, 226(3): 785-797.

        [4]褚志剛,蔣忠翰,周亞男,等. 冰箱風(fēng)扇質(zhì)量在線檢測(cè)系統(tǒng)設(shè)計(jì)[J]. 振動(dòng)與沖擊, 2012, 31(18): 96-99.

        CHU Zhi-gang, JIANG Zhong-han, ZHOU Ya-nan, et al. Design of an online detection system for vibration quality of refrigerator fans[J]. Journal of Vibration and Shock, 2012, 31(18): 96-99.

        [5]陳建良,金濤,孟曉宏,等. 冰箱壓縮機(jī)殼體噪聲輻射數(shù)值分析[J]. 浙江大學(xué)學(xué)報(bào):工學(xué)版, 2007, 41(5): 794-798.

        CHEN Jiang-liang, JIN Tao, MENG Xiao-hong,et al. Numerical analysis of noise radiation from refrigerator compressor shell[J]. Journal of Zhejiang University:Engineering Science, 2007, 41(5): 794-798.

        [6]仲崇明,萬泉,蔣偉康. 往復(fù)式壓縮機(jī)振動(dòng)的有限元數(shù)值分析與實(shí)驗(yàn)研究[J]. 振動(dòng)與沖擊, 2011, 30(5): 156-160.

        ZHONG Chong-ming, WAN Quan, JIANG Wei-kang.Numerical analysis and tests for vibration response of a reciprocating compressor[J]. Journal of Vibration and Shock, 2011, 30(5): 156-160.

        [7]郭維,劉斌,馮濤. 冰箱壓縮機(jī)機(jī)殼實(shí)驗(yàn)?zāi)B(tài)分析[J]. 噪聲與振動(dòng)控制, 2010 (3): 67-70.

        GUO Wei, LIU Bin, FENG Tao. Experimentalmodal analysis of a refrigeration compressor’s shell[J]. Vibration and Noise Control, 2010(3): 67-70.

        [8]趙科,金濤,童水光,等. 冰箱用動(dòng)磁式直線壓縮機(jī)的動(dòng)態(tài)特性仿真[J]. 浙江大學(xué)學(xué)報(bào):工學(xué)版, 2009, 43(1): 138-142.

        ZHAO Ke, JIN Tao, TONG Shui-guang,et al. Simulation of dynamic characteristics of moving magnet linear compressor for refrigerator[J]. Journal of Zhejiang University: Engineering Science, 2009, 43(1): 138-142.

        [9]Kim H, Cho M G, Park J, et al. Prediction of gas pulsation of an industrial compressor[J]. Journal of Central South University, 2013, 20(10): 2735-2740.

        [10]Lee S H, Ryu S M, Jeong W B. Vibration analysis of compressor piping system with fluid pulsation[J]. Journal of Mechanical Science and Technology, 2010, 26(12): 3903-3909.

        [11]劉益才,曹立宏,楊智輝,等. 冰箱毛細(xì)管出口氣液兩相流理論[J]. 中南大學(xué)學(xué)報(bào):自然科學(xué)版, 2007, 38(3): 450-453.

        LIU Yi-cai, CAO Li-hong, YANG Zhi-hui,et al. Theoretical research on gas-liquid two-phase flow at outlet of refrigerator capillary[J]. Journal of Central South University: Science and Technology, 2007, 38(3): 450-453.

        [12]Han H S, Jeong W B, Kim M S, et al. Reduction of the refrigerant-induced noise from the evaporator-inlet pipe in a refrigerator[J]. International Journal of Refrigerator, 2010, 33(7): 1478-1488.

        [13]李春銀,王樹林. 汽車空調(diào)旋葉式壓縮機(jī)排氣閥片的振動(dòng)特性[J]. 振動(dòng)與沖擊, 2014, 33(8): 186-191.

        LI Chun-yin, WANG Shu-lin.Vibration behavior of a discharge valve for vane compressor of a car air conditioner[J]. Journal of Vibration and Shock, 2014, 33(8): 186-191.

        [14]徐斌,馮全科,余小玲. 壓縮機(jī)復(fù)雜管路壓力脈動(dòng)即管道振動(dòng)研究[J]. 核動(dòng)力工程, 2008, 29(4): 79-83.

        XU Bin, FENG Quan-ke, YU Xiao-ling.Study on pressure pulsation and piping vibration of complex piping of reciprocating compressor[J]. Nuclear Power Engineering, 2008, 29(4): 79-83.

        [15]楊國安. 機(jī)械振動(dòng)基礎(chǔ)[M]. 北京:中國石化出版社, 2012.

        [16]盛美萍,王敏慶,孫進(jìn)才. 噪聲與振動(dòng)控制技術(shù)基礎(chǔ)[M].北京:科學(xué)出版社, 2007.

        [17]舒歌群,高文志,劉月輝. 動(dòng)力機(jī)械振動(dòng)與噪聲[M]. 天津:天津大學(xué)出版社, 2008.

        猜你喜歡
        模態(tài)分析固有頻率冰箱
        冰箱里的那些“食”
        現(xiàn)場(chǎng)測(cè)定大型水輪發(fā)電機(jī)組軸系的固有頻率
        冰箱里的奪命派對(duì)
        冰箱是誰發(fā)明的
        肩上冰箱
        大學(xué)生(2017年9期)2017-09-08 13:51:21
        基于ANSYS workbench六片斜葉圓盤渦輪攪拌器的模態(tài)分析
        基于Ansys的礦用局部通風(fēng)機(jī)葉輪模態(tài)分析
        某調(diào)速型液力偶合器泵輪的模態(tài)分析
        民用飛機(jī)沖壓渦輪機(jī)的動(dòng)剛度分析
        科技視界(2015年25期)2015-09-01 16:34:55
        總溫總壓測(cè)頭模態(tài)振型變化規(guī)律研究
        欧美疯狂做受xxxx高潮小说| 美女被内射很爽的视频网站| 亚洲一区二区三区在线看| 亚洲一区二区三区乱码在线中国| 色欲网天天无码av| 亚洲精品无码久久久久久| 国产精品美女白浆喷水| 国产成人精品视频网站| 亚洲中文字幕一二区精品自拍| 国产剧情一区二区三区在线| 亚洲精品成人网站在线播放| 国产精品美女久久久久 | аⅴ资源天堂资源库在线| 精品无码久久久久久久动漫| 久久99国产亚洲高清观看首页| 91久久国产露脸国语对白| 色熟妇人妻久久中文字幕| 亚洲最大av网站在线观看| 日日鲁鲁鲁夜夜爽爽狠狠视频97| 无码av一区在线观看| 国产91成人自拍视频| 国产精品自线一区二区三区| 亚洲av无码乱码国产精品| 成人性生交大片免费看r| 久九九久视频精品网站| 国产美女高潮流白浆视频| 久久久久亚洲精品无码系列| 国产无遮挡又黄又爽在线观看| 动漫在线无码一区| 国产精品亚洲精品日韩动图| 开心久久婷婷综合中文字幕| 国内精品久久久久久久97牛牛| 一本一道波多野结衣一区| 亚洲无线码一区在线观看| 国产黄片一区二区三区| 伊人久久大香线蕉av不变影院| 日本japanese丰满多毛| 亚洲AV无码成人精品区网页| 真人直播 免费视频| 国产综合第一夜| 男的和女的打扑克的视频|