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

        ?

        基于附加約束的小型渦輪增壓器轉(zhuǎn)子模態(tài)測試

        2017-06-15 17:02:23賓光富
        中國機械工程 2017年11期
        關(guān)鍵詞:增壓器振型渦輪

        賓光富 周 偉 王 鋼 陳 文

        湖南科技大學(xué)機械設(shè)備健康維護湖南省重點實驗室,湘潭,411201

        ?

        基于附加約束的小型渦輪增壓器轉(zhuǎn)子模態(tài)測試

        賓光富 周 偉 王 鋼 陳 文

        湖南科技大學(xué)機械設(shè)備健康維護湖南省重點實驗室,湘潭,411201

        針對轉(zhuǎn)子模態(tài)參數(shù)難以有效獲取的問題,進行了基于附加約束的小型渦輪增壓器轉(zhuǎn)子模態(tài)測試研究。以多種材料和零部件組成的某型汽油發(fā)動機渦輪增壓器轉(zhuǎn)子結(jié)構(gòu)為例,采用有限元法構(gòu)建了其有限元模型,得到自由邊界條件下轉(zhuǎn)子前三階模態(tài)頻率。然后分別采用臥式懸掛、立式懸掛、海綿墊、硬紙墊四種不同附加約束條件,開展轉(zhuǎn)子的力錘模態(tài)測試,通過穩(wěn)態(tài)圖識別模態(tài)頻率,與有限元仿真結(jié)果進行對比發(fā)現(xiàn),采用懸掛方式較采用直接固定約束方式的識別值要低,而采用立式懸掛約束方式獲得的轉(zhuǎn)子模態(tài)測試綜合效果最佳。

        模態(tài)參數(shù);附加約束;渦輪增壓器轉(zhuǎn)子;錘擊模態(tài)測試;穩(wěn)態(tài)圖

        0 引言

        渦輪增壓器是一種能夠在不改變發(fā)動機體積和質(zhì)量的前提下,提高發(fā)動機功率、降低排放,同時還能降低油耗的裝置。目前渦輪增壓器在車用發(fā)動機、船舶、航空航天等領(lǐng)域的應(yīng)用越來越廣泛。轉(zhuǎn)子系統(tǒng)是渦輪增壓器的核心部件,它一般由不同材質(zhì)的渦輪、葉輪、轉(zhuǎn)軸等零部件構(gòu)成,經(jīng)常在高速、高溫下工作,渦輪端的溫度在600 ℃左右,工作轉(zhuǎn)速在幾萬到幾十萬r/min之間,屬于典型小型輕質(zhì)高速旋轉(zhuǎn)機械[1],所以在設(shè)計上不僅要求渦輪增壓器轉(zhuǎn)子具有高的強度和可靠性,同時要求具備好的動力學(xué)特性,以滿足渦輪增壓器的功能和壽命條件。模態(tài)頻率是渦輪增壓器轉(zhuǎn)子動力學(xué)設(shè)計的重要參數(shù),模態(tài)頻率的精確度和準確性直接影響到渦輪增壓器轉(zhuǎn)子動力學(xué)特性的評定[2]。由于渦輪增壓器轉(zhuǎn)子常要求快速響應(yīng)、高速運轉(zhuǎn),加上轉(zhuǎn)子結(jié)構(gòu)輕巧緊湊,工作在較高溫度與壓力的密閉型腔內(nèi),難以通過直接安裝傳感器有效測量其工作動態(tài)響應(yīng)[3-4],因此,正確的自由模態(tài)測試方法和精確的模態(tài)分析技術(shù)是獲得高質(zhì)量模態(tài)頻率等參數(shù)的基礎(chǔ),也是小型高速柔性轉(zhuǎn)子仿真計算的重要驗證途徑[5]。在模態(tài)試驗中,測試方法、支撐約束方式、傳感器位置以及測試傳感器的附加質(zhì)量等均會影響轉(zhuǎn)子模態(tài)參數(shù)的精確識別[6-7]。任軍等[8]開展了模態(tài)測試中力傳感器附加質(zhì)量辨識及消除方法研究。林賢坤等[9]探討了針對橋梁結(jié)構(gòu)的基于附加質(zhì)量的試驗?zāi)B(tài)振型質(zhì)量歸一化方法。趙俊等[10]研究了結(jié)構(gòu)模態(tài)測試中傳感器位置優(yōu)化方案。在邊界約束條件對模態(tài)參數(shù)影響方面,蔡力鋼等[11]進行了多約束狀態(tài)下主軸有限元建模及模態(tài)分析。程小勇等[12]研究了初始安裝應(yīng)力對管道固有頻率的影響。而支撐約束方式對轉(zhuǎn)子模態(tài)頻率測試的影響分析,尤其是針對小型渦輪增壓器轉(zhuǎn)子結(jié)構(gòu)的研究目前相對較少。

        本文以某型汽油發(fā)動機渦輪增壓器轉(zhuǎn)子為例,開展這類帶約束的小型轉(zhuǎn)子自由模態(tài)測試研究。通過設(shè)置臥式懸掛、立式懸掛、海綿墊、硬紙墊四種不同附加約束條件,分別開展錘擊法模態(tài)測試研究,將測試結(jié)果與仿真計算進行對比,從而找出最佳測試方案,為轉(zhuǎn)子模態(tài)仿真參數(shù)驗證以及實際工況下轉(zhuǎn)子振動特性分析提供參考。

        1 基于有限元的轉(zhuǎn)子自由模態(tài)仿真分析

        渦輪增壓器轉(zhuǎn)子總長為152 mm,質(zhì)量為215.7 g,且渦輪盤的質(zhì)量屬性和外形均不同于壓氣機盤,轉(zhuǎn)子重心靠近渦輪盤端軸承,即轉(zhuǎn)子具有軸向非對稱特點。該轉(zhuǎn)子主要由渦輪端葉輪、壓氣端葉輪、轉(zhuǎn)軸、鎖緊螺母四部分組成。其中,渦輪端葉輪材料為K418,壓氣端葉輪材料為ZL105,轉(zhuǎn)軸和鎖緊螺母材料為42CrMo。利用CAD軟件對轉(zhuǎn)子進行實體建模,將三維實體模型導(dǎo)入到有限元軟件中,由于模態(tài)分析主要關(guān)注轉(zhuǎn)子臨界轉(zhuǎn)速、振型等宏觀特性,故可將轉(zhuǎn)子各零部件間結(jié)合部位進行簡化等效,直接利用ANSYS的Workbench平臺中自動劃分網(wǎng)格功能進行轉(zhuǎn)子網(wǎng)格劃分,結(jié)果如圖1所示。本文只考慮轉(zhuǎn)子自由邊界條件下的模態(tài),故忽略軸承和基礎(chǔ)的影響;同時考慮到壓氣機葉輪前端部螺紋、葉片圓角等容易在網(wǎng)格劃分過程中產(chǎn)生畸變,但對臨界轉(zhuǎn)速影響很小,故忽略其特征[13-14]。

        圖1 渦輪增壓器轉(zhuǎn)子有限元網(wǎng)格Fig.1 FE model of turbocharger rotor

        根據(jù)建立的渦輪增壓器轉(zhuǎn)子有限元模型,采用Modal求解器分析轉(zhuǎn)子前三階自由模態(tài)頻率值(ω1、ω2、ω3分別為974.0 Hz、1443.7 Hz、2511.8 Hz)對應(yīng)的模態(tài)振型結(jié)果如圖2所示。顯然,第一、第三階振型轉(zhuǎn)軸上應(yīng)變能占比大,為典型轉(zhuǎn)子彎曲振型,而第二階振型則主要在壓氣機葉輪上,轉(zhuǎn)軸上幾乎為零,表現(xiàn)為典型的扭轉(zhuǎn)振型。

        (a)ω1 = 974.0 Hz

        (b)ω2 = 1443.7 Hz

        (c)ω3 = 2511.8 Hz圖2 渦輪增壓器轉(zhuǎn)子前三階模態(tài)振型Fig.2 First three modes calculated by FE model

        2 基于附加約束的錘擊模態(tài)試驗與分析

        2.1 模態(tài)試驗方案

        采用丹麥B&K PULSE噪聲、振動多功能分析儀進行渦輪增壓器轉(zhuǎn)子模態(tài)頻率錘擊法試驗,試驗裝置如圖3所示。其中采用的拾振器為ICP加速度傳感器,靈敏度為1.02 mV/(m·s-2),頻響范圍為1~10 kHz;激勵裝置為8206-002型脈沖力錘,選用的鋁合金錘帽頻率帶寬為0~5 kHz;數(shù)據(jù)采集卡為PULSE前端系統(tǒng)LAN-XI,具體型號為3053-B120,分析頻率范圍為0~25.6 kHz;采用Reflex模態(tài)測試分析處理軟件,選取有理分式多項式擬合法進行穩(wěn)態(tài)圖中模態(tài)參數(shù)識別。

        圖3 錘擊法模態(tài)試驗設(shè)備Fig.3 Impact modal testing setup for turbocharger rotor

        由于渦輪增壓器轉(zhuǎn)子的固有頻率值較高,故選用較硬的鋁合金錘帽進行模態(tài)錘擊試驗??紤]到轉(zhuǎn)子結(jié)構(gòu)尺寸小,測試空間有限,若布置多個拾振傳感器則會因過大的傳感器附件質(zhì)量而對測試結(jié)果產(chǎn)生較大偏差,本文采用多點激勵、單點響應(yīng)的錘擊測試方案,將拾振傳感器布置于轉(zhuǎn)軸中間位置,用強力膠黏結(jié)固定,以避免力錘沖擊激勵下傳感器與轉(zhuǎn)子間間隙松動造成能量損失,拾振點和錘擊點的布置位置如圖4所示。在模態(tài)測試試驗過程中,每次盡可能以較快的錘擊速度來提高激振能量,且每個錘擊點使用大致相同大小的力各錘擊5次,然后對測試響應(yīng)信號進行平均化處理,采用指數(shù)窗函數(shù),以減少噪聲、人為等隨機因素對測試結(jié)果造成的誤差,提高測試精度。

        圖4 轉(zhuǎn)子模態(tài)測試測點布置Fig.4 Impact location and sensor placement of rotor

        2.2 基于附加約束的錘擊模態(tài)試驗

        針對臥式懸掛、立式懸掛、海綿墊、硬紙墊四種不同附加約束條件,分別開展錘擊法模態(tài)測試研究,對渦輪增壓器轉(zhuǎn)子模態(tài)頻率參數(shù)進行識別,以便將測試結(jié)果與有限元仿真結(jié)果進行對比分析,找出最佳附加約束條件下的轉(zhuǎn)子自由模態(tài)測試方案。

        自由模態(tài)測試過程中,被測結(jié)構(gòu)物附加約束點位置一般根據(jù)所關(guān)注模態(tài)階次盡可能選擇在模態(tài)節(jié)點處。由于渦輪增壓器轉(zhuǎn)子前三階模態(tài)變形主要發(fā)生在葉輪及轉(zhuǎn)軸上,為便于力錘敲擊測試,結(jié)合轉(zhuǎn)子的實際情況,自由懸掛方式采用彈性繩臥式與立式約束兩種方案進行試驗,如圖5所示。考慮到空間有限,為便于力錘敲擊,臥式懸掛約束點分別位于轉(zhuǎn)子兩端螺母處,而立式懸掛約束點則位于葉輪端螺母處,以使渦輪端位于底部。由于渦輪增壓器轉(zhuǎn)子小、質(zhì)量輕,考慮到懸掛約束方式下力錘敲擊會導(dǎo)致轉(zhuǎn)子沿著錘擊方向擺動,為限制這種擺動以及避免“連擊”現(xiàn)象,采用直接固定附加約束測試方案,將轉(zhuǎn)子分別水平臥式置于海綿墊、硬紙墊兩種較大附加支撐約束上,如圖6所示。

        接下來分別采用多點激勵單點響應(yīng)方法進行模態(tài)測試,分析得到4種附加約束條件下的傳遞函數(shù)的頻響函數(shù),并進行幾何加權(quán)獲得穩(wěn)態(tài)圖,自動識別模態(tài)參數(shù)結(jié)果如圖7所示。圖中“▽”表示頻率和特征向量穩(wěn)定,“*”表示頻率和阻尼穩(wěn)定,“×”表示頻率穩(wěn)定,“◇”表示頻率、特征向量、阻尼都穩(wěn)定。

        (a)水平懸掛

        (b)立式懸掛圖5 基于彈性繩自由懸掛約束的模態(tài)測試Fig.5 Free hanging with elastic rope for impact testing

        (a)海綿墊支撐約束

        (b)硬紙墊支撐約束圖6 基于支撐約束的模態(tài)測試Fig.6 Fixed supporting for impact testing

        2.3 試驗結(jié)果分析與討論

        從渦輪增壓器轉(zhuǎn)子采用自由懸掛約束方式測試結(jié)果來看,圖7a中在1400 Hz附近沒有明顯的峰值,而圖7b中各階波峰則清晰,說明采用彈性繩臥式懸掛會存在部分模態(tài)丟失,而立式懸掛則能較全面地識別出各階模態(tài)。從采用直接固定附加約束測試方案獲得的圖7c、圖7d來看,在1100 Hz左右處的頻響函數(shù)峰值均很明顯,而1400 Hz附近則均沒有明顯峰值,這可能與附加約束剛度值較大有關(guān)。將這四種測試方案所識別的模態(tài)頻率值與仿真值進行對比及對相對誤差率δ進行分析,結(jié)果如表1和表 2所示。

        (a)臥式懸掛約束

        (d)硬紙墊支撐約束圖7 穩(wěn)態(tài)圖與模態(tài)頻率識別結(jié)果Fig.7 Stability diagram and modal frequencies identify

        表1 四種約束下模態(tài)頻率匯總

        表2 測試值與仿真值的相對誤差率δ

        從表1與表2中可以看出采用臥式懸掛、海綿墊和硬紙墊三種附加約束方式測試所得第一階模態(tài)頻率值ω1與有限元計算值相差較大,相對誤差率在10%以上,只有立式懸掛約束所得測試值與仿真值誤差為4.5%,且直接固定約束比懸掛約束的誤差要大,這與懸掛約束會使轉(zhuǎn)子更軟這一因素有關(guān);ω2只有立式懸掛和海綿墊才有識別值,這可能與第二階為扭轉(zhuǎn)模態(tài)振型有關(guān);ω3的誤差有所減小,且直接固定約束誤差比懸掛約束的誤差也要大。因此,從試驗結(jié)果來看:采用懸掛方式測試結(jié)果較直接固定約束方式所測得結(jié)果要更優(yōu),其中又以彈性繩立式懸掛約束方式的測試結(jié)果最佳。這也說明采用立式懸掛對轉(zhuǎn)子約束最小,與自由模態(tài)邊界條件分析結(jié)果最吻合。

        3 結(jié)論

        (1)通過分別采用臥式懸掛、立式懸掛、海綿墊、硬紙墊四種不同附加約束條件進行渦輪增壓器轉(zhuǎn)子模態(tài)試驗測試,測試結(jié)果有效驗證了有限元仿真模型,同時找出了最佳約束的測試方案。

        (2) 采用懸掛法測試識別的轉(zhuǎn)子模態(tài)頻率值較直接固定約束方式識別的轉(zhuǎn)子模態(tài)頻率值要偏小,這與采用不同約束方式下給測試轉(zhuǎn)子帶來的附加約束有關(guān),即懸掛約束相對直接固定約束會使轉(zhuǎn)子變得更軟。而第二階的相對誤差率更高,這可能與轉(zhuǎn)子扭轉(zhuǎn)模態(tài)振型有關(guān)。

        (3) 采用立式懸掛約束所測前三階模態(tài)頻率結(jié)果與自由邊界條件下有限元仿真結(jié)果最匹配。說明采用約束最小的立式懸掛法對轉(zhuǎn)子模態(tài)測試影響最小,與真實結(jié)果最接近。本文測試研究可為這類轉(zhuǎn)子動力學(xué)設(shè)計和分析提供參考??紤]到實際工作中軸承、轉(zhuǎn)速等因素的影響,今后還需結(jié)合實際工況進行工作運行模態(tài)分析。

        [1] KIRK R G. Experimental Evaluation of Hydrodynamic Bearings for a High Speed Turbocharger[J]. Journal of Engineering for Gas Turbines and Power,2014,136(7):072501.

        [2] 李紅偉,于文濤,劉淑琴. 基于ANSYS的磁懸浮撓性轉(zhuǎn)子模態(tài)分析與設(shè)計[J]. 中國機械工程,2014,25(11):1447-1452. LI Hongwei,YU Wentao,LIU Shuqin. Mode Analysis and Design of Maglev Flexible Rotor with ANSYS[J]. China Mechanical Engineering,2014,25(11):1447-1452.

        [3] WANG Longkai, BIN Guangfu, LI Xuejun, et al. Effects of Unbalance Location on Dynamic Characteristics of High-speed Gasoline Engine Turbocharger with Floating Ring Bearings[J]. Chinese Journal of Mechanical Engineering, 2016, 29(2):271-280.

        [4] 應(yīng)廣馳,孟光,龍新華,等. 用自由界面模態(tài)綜合法研究發(fā)動機-排氣管-渦輪增壓器系統(tǒng)的振動[J]. 機械強度,2010,32(1):10-16. YING Guangchi,MENG Guang,LONG Xinhua,et al. Study on Engine Manifold Turbocharger System’s Vibration by Free Interface Modal Synthesis Method[J]. Journal of Mechanical Strength,2010,32(1):10-16.

        [5] 羅君杰,龍新華,孟光. 基于實驗?zāi)B(tài)分析的渦輪增壓器壓縮端振動響應(yīng)分析[J]. 噪聲與振動控制,2010,30(5):10-14. LUO Junjie,LONG Xinhua,MENG Guang. Vibration Response Analysis of Turbocharger Compressor Housing Based on Experimental Modal Analysis[J]. Noise and Vibration Control,2010,30(5):10-14.

        [6] 董小瑞,楊世文. 復(fù)雜結(jié)構(gòu)模態(tài)測試方法研究[J]. 兵工學(xué)報,2008,29(4):474-477. DONG Xiaorui,YANG Shiwen. Research on Modal Measurement Method of Complicated Structure[J]. Acta Armamentarii, 2008,29(4):474-477.

        [7] 羅忠,王宇,孫寧,等. 不同邊界條件下旋轉(zhuǎn)薄壁短圓柱殼的強迫振動響應(yīng)計算[J]. 機械工程學(xué)報,2015,51(9):64-70. LUO Zhong,WANG Yu,SUN Ning,et al. Forced Vibration Response Calculation of Rotating Short Thin Cylindrical Shells for Various Boundary Conditions[J]. Journal of Mechanical Engineering,2015,51(9):64-70.

        [8] 任軍,畢樹生. 模態(tài)測試中力傳感器附加質(zhì)量辨識及消除方法研究[J]. 振動與沖擊,2014,33(14): 108-112. REN Jun,BI Shusheng. Identification and Correction of Force Transducer Mass Effects in Modal Testing[J]. Journal of Vibration and Shock,2014,33(14):108-112.

        [9] 林賢坤,覃柏英,張令彌,等. 基于附加質(zhì)量的試驗?zāi)B(tài)振型質(zhì)量歸一化[J]. 振動、測試與診斷, 2012,32(5):784-790. LIN Xiankun,QIN Boying,ZHANG Lingmi,et al. Way of Getting Mass-normalized Experimental Mode Shapes Based on Mass Changes[J]. Journal of Vibration, Measurement & Diagnosis,2012,32(5):784-790.

        [10] 趙俊,聶振華,馬宏偉. 拱結(jié)構(gòu)模態(tài)測試中傳感器優(yōu)化配置[J]. 振動、測試與診斷,2011,31(2):217-222. ZHAO Jun,NIE Zhenhuai,MA Hongwei. Sensor Placement Optimization for Modal Test of Arch Structure[J]. Journal of Vibration Measurement & Diagnosis,2011,31(2):217-222.

        [11] 蔡力鋼,馬仕明,趙永勝,等. 多約束狀態(tài)下重載機械式主軸有限元建模及模態(tài)分析[J]. 機械工程學(xué)報, 2012,48(3):165-173. CAI Ligang,MA Shiming,ZHAO Yongsheng,et al. Finite Element Modeling and Modal Analysis of Heavy-duty Mechanical Spindle under Multiple Boundary conditions[J]. Journal of Mechanical Engineering,2012,48(3):165-173.

        [12] 程小勇,陳果,劉明華,等. 初始安裝應(yīng)力對管道固有頻率的影響分析及試驗驗證[J]. 中國機械工程,2015,26(4):512-517. CHENG Xiaoyong,CHEN Guo,LIU Minghua,et al. Analysis and Experimental Verification to Effects of Pipe Initial Installation Stress on Pipe’s Natural Frequencies[J]. China Mechanical Engineering,2015,26(4):512-517.

        [13] 李元生,李磊,敖良波,等. 渦輪增壓器轉(zhuǎn)子系統(tǒng)三維動力學(xué)特性分析[J]. 機械強度,2011,33(2):175-178. LI Yuansheng,LI Lei,AO Liangbo,et al. Analysis of the Three-dimensional Rotor Dynamics Characteristic of the Turbocharger[J]. Journal of Mechanical Strength,2011,33(2):175-178.

        [14] MA Hui,LI Hui,NIU Heqiang,et al. Numerical and Experimental Analysis of the First and Second-mode Instability in a Rotor-bearing System[J]. Archive of Applied Mechanics,2014,84(4):519-541.

        (編輯 王艷麗)

        Modal Testing for Small-size Turbocharger Rotor with Additional Constraints

        BIN Guangfu ZHOU Wei WANG Gang CHEN Wen

        Health Maintenance for Mechanical Equipment Key Lab of Hunan Province, Hunan University of Science and Technology, Xiangtan,Hunan,411201

        As modal parameters usually were difficult to obtain for small-size turbocharger rotors, a rotor modal testing with additional constraints was proposed. A gasoline engine turbocharger rotor including multi materials and components was used as the example to demonstrate the proposed approach. A finite element (FE) model for turbocharger rotor with free boundary conditions was built to calculate the modal frequencies for the former three orders. Then impact modal testing was carried out based on four different additional constraints including horizontal hanging, vertical hanging, supporting on spongy cushion, and supporting on hardboard. The modal frequencies were evaluated from the stability diagrams. It is found the recognition values from constraint way of hanging are lower than that of the constraint way of supporting compared to the FE simulation. However, the modal frequencies identified under the vertical hanging constraint is the best for all kinds of constraints.

        modal parameter; additional constraint; turbocharger rotor; hammer impact testing; stability diagram

        2016-07-04

        國家自然科學(xué)基金資助項目(51575176,11672106);湖南省教育廳優(yōu)秀青年項目(15B085);湖南省科技重大專項(2015GK1003)

        TH212;TH213.3

        10.3969/j.issn.1004-132X.2017.11.004

        賓光富,男,1981年生。湖南科技大學(xué)機械設(shè)備健康維護湖南省重點實驗室副教授。主要研究方向為旋轉(zhuǎn)機械動力學(xué)與振動控制。獲發(fā)明專利12項。發(fā)表論文50余篇。E-mail:abin811025@163.com。周 偉,男,1991年生。湖南科技大學(xué)機械設(shè)備健康維護湖南省重點實驗室碩士研究生。王 鋼,男,1971年生。湖南科技大學(xué)機械設(shè)備健康維護湖南省重點實驗室教授。陳 文,男,1990年生。湖南科技大學(xué)機械設(shè)備健康維護湖南省重點實驗室碩士研究生。

        猜你喜歡
        增壓器振型渦輪
        關(guān)于模態(tài)綜合法的注記
        縱向激勵下大跨鋼桁拱橋高階振型效應(yīng)分析
        2014款寶馬525Li渦輪增壓壓力過低
        塔腿加過渡段輸電塔動力特性分析
        博格華納由可變截面渦輪增壓器向電子渦輪增壓器發(fā)展
        結(jié)構(gòu)振型幾何辨識及應(yīng)用研究
        山西建筑(2015年14期)2015-06-05 09:37:07
        渦輪增壓發(fā)動機與雙離合變速器的使用
        Opel公司新型1.0L渦輪增壓直接噴射汽油機
        小型化發(fā)動機用新型羅茨增壓器
        采用兩級渦輪增壓器提高功率密度
        久久免费网站91色网站| 国产男女免费完整视频| 国产成人免费一区二区三区| 国产精品亚洲一区二区杨幂| 亚洲av影片一区二区三区| 成人在线观看av毛片| 亚洲av无码专区亚洲av网站| 人妻少妇av无码一区二区| 亚洲国产欧美久久香综合| 白色白在线观看免费2| 刺激一区仑乱| 无套内谢的新婚少妇国语播放| 亚洲国产欧美另类va在线观看| 国产人妖直男在线视频| 国产玉足榨精视频在线观看| 午夜三级a三级三点| 国产亚洲第一精品| 亚洲综合小综合中文字幕| 麻豆国产精品va在线观看不卡| 亚洲色欲色欲www在线播放| 国产va在线播放| 国产一区二区中文字幕在线观看| 久久久久久欧美精品se一二三四| 久久丫精品国产亚洲av| 制服丝袜中文字幕在线| 亚洲国产区男人本色| 97久久成人国产精品免费| 久久综合久久综合久久| 屁屁影院ccyy备用地址 | 97精品人人妻人人| 亚洲另类自拍丝袜第五页| 精品黄色av一区二区三区| 在线视频色系中文字幕| 久久久久久国产精品mv| 免费国产黄线在线播放| 亚洲av日韩精品一区二区| 四虎国产精品永久在线| 这里有精品可以观看| 国产在线视频网站不卡| 人人妻人人澡人人爽精品日本| 无码手机线免费观看|