999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

空調(diào)室外機系統(tǒng)的實驗?zāi)B(tài)分析與模型更新

2019-09-02 08:39:54王栢村張博鈞陳冠銂梁秀瑋
聲學(xué)技術(shù) 2019年4期
關(guān)鍵詞:模態(tài)有限元實驗

王栢村,張博鈞,陳冠銂,梁秀瑋

空調(diào)室外機系統(tǒng)的實驗?zāi)B(tài)分析與模型更新

王栢村1,張博鈞1,陳冠銂1,梁秀瑋2

(1. 屏東科技大學(xué)機械工程系,臺灣屏東 91201;2.瑞智精密股份有限公司,臺灣桃園 32849)

有限元(Finite Element, FE)模型的正確性和可靠性對于確保結(jié)構(gòu)仿真達到其分析目的相當(dāng)重要。文章介紹了空調(diào)室外機(Outdoor Unit of Conditioner, OUC)有限元模型的建模技術(shù)與模型更新。建構(gòu)OUC的有限元模型必須呈現(xiàn)結(jié)構(gòu)的物理特性,包括幾何形狀、材料性質(zhì)、接觸接口和邊界條件等,并施以實驗進行驗證。一般的方法是對結(jié)構(gòu)進行實驗?zāi)B(tài)分析(Experimental Modal Analysis, EMA)以獲得結(jié)構(gòu)模態(tài)參數(shù)。同時,也可由有限元模型的數(shù)值分析求得理論模態(tài)參數(shù)。模型更新是調(diào)整有限元模型參數(shù),使分析模型和實際結(jié)構(gòu)的結(jié)構(gòu)模態(tài)特性相符。結(jié)果顯示,更新后的有限元模型在結(jié)構(gòu)模態(tài)和頻率域特性上,能充分呈現(xiàn)空調(diào)室外機實際結(jié)構(gòu)的振動特性。

空調(diào)室外機;實驗?zāi)B(tài)分析;有限元分析;模型更新

0 引 言

空調(diào)室外機(Outdoor Unit of Conditioner, OUC)即分離式冷氣機的室外機體,其振動噪聲往往通過底盤、梁架、管路、墻面?zhèn)鬟f至室內(nèi),影響到自己或鄰居,引起身心的不舒適,因此在開發(fā)低振動噪聲室外機的過程中,需要通過計算機輔助工程分析并結(jié)合模型驗證程序,探討實際室外機結(jié)構(gòu)的振動模態(tài)特性,且有助于OUC的振動噪聲診斷;雖然電腦輔助工程(Computer Aided Engineering, CAE)為功能強大的應(yīng)用工具,但仍然需要實驗的驗證與認證。本文為對OUC之分析模型進行模型驗證以確認分析模型的可靠度,透過分析模型觀察空調(diào)室外機內(nèi)部結(jié)構(gòu)的振動特性,找出壓縮機產(chǎn)生異響的可能零件,以輔助空調(diào)室外機的噪聲診斷。

王栢村等[1]運用理論模態(tài)分析與實驗?zāi)B(tài)分析,對空調(diào)室外機系統(tǒng)進行理論與實驗測量分析,結(jié)合模型驗證的程序,比對模態(tài)分析所獲得的模態(tài)參數(shù),確認空調(diào)室外機分析模型與實際結(jié)構(gòu)之間的等效性,最終獲得等效于實際結(jié)構(gòu)的空調(diào)室外機分析模型,借此探討室外機內(nèi)部各次系統(tǒng)的振動模態(tài)特性,有助于對空調(diào)室外機系統(tǒng)的振動噪聲進行診斷。丘曉宏等[2]對除濕機底盤進行理論模態(tài)分析,以除濕機壓縮機的倍頻與底盤的前三個模態(tài)作為除濕機底盤共振頻率與底盤強度的參考,最后以調(diào)整加強條、工字筋、凹槽結(jié)構(gòu)幾何特征方式改善底盤的強度與固有頻率,借此消除壓縮機與底盤共振產(chǎn)生的噪聲問題。SILVEIRA等[3]對往復(fù)式壓縮機進行振動與噪聲測量,并量化和確定噪聲的發(fā)生源,同時也對振動傳遞到系統(tǒng)的方式進行探討,準(zhǔn)確地描述所有的噪聲問題,最后通過新開發(fā)的組合件來變更現(xiàn)有的產(chǎn)品,達到靜音的效果。PARK等[4]通過特征值的拓撲優(yōu)化設(shè)計來降低壓縮機產(chǎn)生的噪聲,主要以工程結(jié)構(gòu)噪聲改善的實際案例為例,由實驗?zāi)B(tài)分析與噪聲測量得知,模態(tài)振型與高噪聲有明顯的關(guān)聯(lián)性,因此通過有限元分析之特征值的拓撲優(yōu)化設(shè)計來改善結(jié)構(gòu)的噪聲。KIM等[5]通過表面振動以及聲音強度的測量,確認壓縮機的異響傳遞方式,將壓縮機本體進行次組合的模型驗證程序。第1層次進行泵組合件的驗證,第2層次對轉(zhuǎn)子成品與泵的組合進行驗證,第3層次加入本體外殼進行驗證,第4層次則是加上本體外殼的零組件,最后進行整機壓縮機的模型驗證,共分為5個案例。為了降低壓縮機整體的噪聲量,進行模型的修整,于過濾瓶新增L型支架以固定內(nèi)管結(jié)構(gòu),并在過濾瓶的本體加上橡膠墊圈,整體模型的修改使得壓縮機的噪聲頻譜下降了20 dB。MAGLIE等[6]介紹了工具機從精密復(fù)雜的有限元模型到簡化表示的過程,在計算時間和資源方面的效率均有顯著提升,并應(yīng)用自動化求解,通過不同的分析驗證,獲得工具機的簡化有限元模型。

本文運用有限元分析與實驗?zāi)B(tài)分析,對OUC進行分析模型更新,主要目標(biāo)如下所示:

(1) 對OUC系統(tǒng)與次系統(tǒng)的分析模型進行更新,以O(shè)UC內(nèi)部的冷凝器與回轉(zhuǎn)式壓縮機為更新對象,同時以實驗?zāi)B(tài)分析進行驗證;確認理論與實際結(jié)構(gòu)之間的等效性。

(2) 將更新完成的冷凝器與回轉(zhuǎn)式壓縮機分析模型加入OUC分析模型,再次進行模型的驗證,確認結(jié)構(gòu)的物理特性,包括幾何形狀、材料性質(zhì)、接觸接口和邊界條件等。以零件、次系統(tǒng)、整機系統(tǒng)逐步驗證的方式,更新OUC的有限元模型,最終獲得等效于實際結(jié)構(gòu)的空調(diào)室外機分析模型。

(3) 未來可借由分析模型探討OUC內(nèi)部各次系統(tǒng)的振動模態(tài)特性,有助于空調(diào)室外機系統(tǒng)的振動噪聲診斷。

1 空調(diào)室外機的實驗?zāi)B(tài)分析

為了確認更新后的OUC模型與實際結(jié)構(gòu)之間的等效性,需要進行實驗驗證。本文主要以實驗?zāi)B(tài)分析為實驗方法,將所獲得的模態(tài)參數(shù)作為校正理論分析模型的基準(zhǔn)。

首先,針對要更新的零件與次系統(tǒng)結(jié)構(gòu)進行實驗?zāi)B(tài)分析,表1為室外機系統(tǒng)與次系統(tǒng)的實驗規(guī)劃總表。實驗?zāi)B(tài)分析的規(guī)劃可分為自由邊界與系統(tǒng)邊界;自由邊界下的實驗?zāi)B(tài)參數(shù),主要作為校正有限元模型的幾何形狀、材料性質(zhì)與接觸條件的基準(zhǔn);另外,系統(tǒng)邊界則是校正次系統(tǒng)結(jié)構(gòu)的邊界參數(shù)。

表1 室外機系統(tǒng)與次系統(tǒng)的實驗規(guī)劃總表

本文主要是對壓縮機與冷凝器進行有限元模型的更新,因此須對實際壓縮機與冷凝器進行實驗?zāi)B(tài)分析,并獲得結(jié)構(gòu)的模態(tài)參數(shù)作為校正理論分析模型的接觸參數(shù);冷凝器在自由邊界下的實驗架設(shè)主要以釣魚線懸吊來模擬自由邊界,以固定沖擊槌移動加速度計的方式進行實驗,布點規(guī)劃則是在冷凝器左右銅管;回轉(zhuǎn)式壓縮機在自由邊界下的實驗架設(shè)上,也是以釣魚線懸吊來模擬自由邊界,運用固定驅(qū)動器沖擊錘,移動傳感器加速規(guī)的實驗方式進行。

然而,我們的“貪得無厭”總歸要有一個天花板,即使銀行按分鐘、按秒來發(fā)放利息,我們一年后最多能獲得的收益就是這樣一個極限,而這個極限正是e。

最后對OUC系統(tǒng)內(nèi)部管路、壓縮機以及外部鈑金進行系統(tǒng)邊界下的實驗?zāi)B(tài)分析,獲得系統(tǒng)邊界下的模態(tài)參數(shù),并作為完整OUC有限元模型的校正基準(zhǔn),以利修正接觸參數(shù)與邊界參數(shù)。OUC系統(tǒng)的實驗方式是采用固定沖擊錘移動三軸向加速度計進行敲擊實驗,通過敲擊實驗結(jié)果可獲得結(jié)構(gòu)的頻率響應(yīng)函數(shù),再經(jīng)由ME'scopeVES軟件截取實驗?zāi)B(tài)參數(shù),包含固有頻率、模態(tài)振型與模態(tài)阻尼比,完成OUC系統(tǒng)與次系統(tǒng)的實驗?zāi)B(tài)分析。

2 空調(diào)室外機的有限元分析模型更新方法

有限元分析模型的等效性,對于響應(yīng)預(yù)測以及設(shè)計變更顯得相當(dāng)?shù)闹匾?,更關(guān)系到結(jié)構(gòu)仿真能否達到分析目的。本文延續(xù)文獻[1]中的空調(diào)室外機系統(tǒng)的有限元模型,對其進行模型變更,使得OUC系統(tǒng)的分析模型在結(jié)構(gòu)模態(tài)和頻率域特性上能充分呈現(xiàn)空調(diào)室外機實際結(jié)構(gòu)的振動特性。本節(jié)將對有限元分析模型的更新方法、分析設(shè)定、模型更新結(jié)果進行說明。

首先說明空調(diào)室外機模型更新方法。圖1為空調(diào)室外機系統(tǒng)與次系統(tǒng)的模型更新架構(gòu)圖,更新的分析模型主要為冷凝器與回轉(zhuǎn)式壓縮機。在冷凝器模型更新上,改善了文獻[1]的冷凝器分析模型幾何上過于簡化的缺點,更新后的冷凝器模型是以實體幾何構(gòu)建冷凝器的銅管,使其在物理性質(zhì)上更貼近于實際結(jié)構(gòu);在回轉(zhuǎn)式壓縮機的更新上,更新后的壓縮機模型加入了實體幾何結(jié)構(gòu),包含馬達轉(zhuǎn)定子、泵組合件以及過濾瓶等重要組合件,主要改善了壓縮機幾何過于簡化的缺點;對更新后的分析模型均用模型驗證程序進行驗證,以實驗?zāi)B(tài)參數(shù)為基準(zhǔn),修正更新的分析模型。接著將完成驗證后的冷凝器與壓縮機模型加入OUC系統(tǒng)的模型內(nèi),替換須改善的零件與次系統(tǒng)結(jié)構(gòu),最后進行理論模態(tài)分析,求解理論模態(tài)參數(shù),以利于后續(xù)與實驗?zāi)B(tài)參數(shù)做比對驗證,確認更新后的OUC模型是否等效于實際結(jié)構(gòu)。

在OUC的理論模態(tài)分析設(shè)定上,將說明空調(diào)室外機系統(tǒng)有限元模型的建立方式、接觸接口的分析假設(shè)、接觸接口的仿真方式,其設(shè)定方式如下:

(1) OUC的有限元模型均采用二次立體元素結(jié)構(gòu)元素(Solid186、Solid187)構(gòu)建,冷凝器管路采用梁元素(188),外殼鈑金采用殼元素(Shell181)來構(gòu)建有限元素模型,質(zhì)量元素(Mass21)仿真未加入風(fēng)扇。

圖1 空調(diào)室外機系統(tǒng)與次系統(tǒng)的模型更新架構(gòu)圖

(2) OUC的材料參數(shù)設(shè)定值則參照文獻[1]中的材料參數(shù),該篇作者是對OUC各零件獨立進行模型驗證,獲得各零件的材料性質(zhì)。

(3) 接觸分析設(shè)定主要可分為Bonded、No Separation以及Joint-fix。Bonded用于模擬焊接與螺栓鎖固,No Separation用于仿真壓縮機內(nèi)部零件間的滑動接觸,Joint-fix用于模擬OUC鈑金間的螺栓鎖固。

(4) 空調(diào)室外機系統(tǒng)與次系統(tǒng)的分析模型均采用自由分割(Free Mesh)方式進行元素分割,共427 281個節(jié)點、203 234個元素。

(6) 進行理論模態(tài)分析,求得OUC系統(tǒng)結(jié)構(gòu)的理論模態(tài)參數(shù),最后也進行簡諧響應(yīng)分析,給予、、三軸方向為1 N的簡諧激振外力,輸出獲得OUC結(jié)構(gòu)三軸方向的頻率響應(yīng)函數(shù)。

3 模型更新結(jié)果與討論

本文更新了OUC系統(tǒng)內(nèi)部的冷凝器與壓縮機結(jié)構(gòu),并進行了實驗?zāi)B(tài)分析、理論模態(tài)分析以及簡諧響應(yīng)分析,獲得了模態(tài)域的固有頻率、模態(tài)振型以及頻率域的頻率響應(yīng)函數(shù)。將模態(tài)域與頻率域的結(jié)果進行比對驗證,確認有限元分析模型的等效性。表2為空調(diào)室外機系統(tǒng)的外部機殼模態(tài)振型,表3為空調(diào)室外機系統(tǒng)的壓縮機模態(tài)振型,表4為空調(diào)室外機系統(tǒng)的管路系統(tǒng)模態(tài)振型。圖2為OUC模態(tài)域與頻率域方向響應(yīng)結(jié)果比較圖,圖3為模態(tài)域與頻率域方向響應(yīng)結(jié)果比較圖。驗證結(jié)果說明如下:

(1) 外部機殼模態(tài)振型的物理意義,實驗與分析的模態(tài)振型均有良好的對應(yīng),理論分析模型的固有頻率接近實驗的固有頻率。

(2) 由表3中的OUC系統(tǒng)下的壓縮機模態(tài)域結(jié)果,可觀察到OUC理論分析模型能完整仿真實際結(jié)構(gòu)單純壓縮機擺動的Local mode現(xiàn)象,固有頻率的比對也相當(dāng)接近。

(3) 通過表4的結(jié)果可知,理論分析模型能夠完全仿真壓縮機做剛體運動而牽動管路擺動的現(xiàn)象;OUC的Global mode條件下,管路的模態(tài)振型特性也具有良好的對應(yīng)關(guān)系,固有頻率的比對也相當(dāng)接近。

(4) 圖2、3為空調(diào)室外機的模態(tài)振型與頻率響應(yīng)函數(shù)(Frequency Response Function, FRF)比較圖,紅色虛線為簡諧響應(yīng)分析求得的頻率響應(yīng)函數(shù),黑色曲線為實驗的頻率響應(yīng)函數(shù),兩曲線的峰值位置為結(jié)構(gòu)的固有頻率,且所對應(yīng)的模態(tài)振型也具有相同的物理意義。

表2 空調(diào)室外機系統(tǒng)的外部機殼的振動模態(tài)

表3 空調(diào)室外機系統(tǒng)壓縮機的振動模態(tài)

表4 空調(diào)室外機管路系統(tǒng)的振動模態(tài)

4 結(jié)論

本文介紹了空調(diào)室外機(OUC)有限元模型的建模技術(shù)與模型更新方法,主要對冷凝器與回轉(zhuǎn)式壓縮機進行模型變更,同時進行模型驗證程序,以實驗?zāi)B(tài)分析所獲得的模態(tài)參數(shù)做為修正分析模型的基準(zhǔn);最后將更新完成的冷凝器與回轉(zhuǎn)式壓縮機分析模型加入OUC分析模型中,再一次進行模型驗證程序;綜合上述結(jié)果與討論,總結(jié)如下:

(1) 室外機(OUC)內(nèi)部的冷凝器與回轉(zhuǎn)式壓縮機進行模型變更,并通過模型驗證程序,確認兩種分析模型與實際間的等效性,完成結(jié)構(gòu)的物理性質(zhì)確認,包括幾何形狀、材料性質(zhì)、接觸接口和邊界條件等。

圖2 空調(diào)室外機結(jié)構(gòu)的頻率相應(yīng)曲線和相應(yīng)的結(jié)構(gòu)z方向位移的振動模態(tài)

圖3 空調(diào)室外機結(jié)構(gòu)的頻率相應(yīng)曲線和相應(yīng)的結(jié)構(gòu)y方向位移的振動模態(tài)

(2) 將更新后的冷凝器與回轉(zhuǎn)式壓縮機模型加入到OUC模型內(nèi),再一次進行模型驗證,在模態(tài)域的結(jié)果表明,無論是外部機殼或是內(nèi)部的壓縮機

與管路系統(tǒng),在模態(tài)參數(shù)的比對上均有良好的對應(yīng)性,表示更新后的OUC模型足以顯示實際OUC系統(tǒng)的結(jié)構(gòu)振動特性。

(3) 在OUC的頻率域比對中,實驗與分析的頻率響應(yīng)函數(shù)的峰值尖銳度有差別,推測為理論分析模型的簡諧響應(yīng)是以Mode Superposition進行計算的,阻尼效應(yīng)的輸入是以一個常數(shù)值進行設(shè)定,導(dǎo)致頻率響應(yīng)函數(shù)的振幅有差別,但峰值頻率所對應(yīng)的模態(tài)振型有良好的對應(yīng)關(guān)系,因此確定完成空調(diào)室外機的模型更新。

(4) 本文進行OUC系統(tǒng)各零件的更新技巧、邊界仿真、接觸型式以及分析模型的驗證方法,為空調(diào)機設(shè)計開發(fā)人員建立空調(diào)系統(tǒng)仿真應(yīng)用提供了標(biāo)準(zhǔn)的作業(yè)流程,減少了設(shè)計開發(fā)所需的時間。

未來可進一步探討OUC分析模型應(yīng)用于振動與噪聲的診斷方法,通過分析模型的幫助觀察到實驗儀器無法測量到的內(nèi)部振動特性,借此輔助OUC的振動與噪聲的診斷。

[1] 王栢村, 張博鈞, 陳冠銂, 等. 空調(diào)室外機系統(tǒng)之模擬技術(shù)與實驗驗證[C]//屏東科技大學(xué)暨北京科技大學(xué)第十二屆學(xué)術(shù)交流研討會. 屏東: 屏東科技大學(xué)出版社, 2017: A15.

[2] 丘曉宏, 溫驪. 基于ANSYS的除濕機底盤強度分析及結(jié)構(gòu)優(yōu)化[J]. 家電科技, 2015, 18(5): 73-75.

QIU Xiaohong, WEN Li. Modal analysis for structure of floor-standing air conditioner[J]. China Appliance Technology, 2015, 18(5): 73-75.

[3] SILVEIRA M. Noise and vibration reduction in compressors for commercial applications[C]//International Compressor Engineering Conference, 2004: C065.

[4] PARK J, WANG S. Noise reduction for compressors by modes control using topology optimization of eigenvalue[J]. Journal of Sound and Vibration, 2008, 315(4-5): 836-848.

[5] KIM J D, LEE B C, BAE J O, KIM Y H. Noise reduction of a rotary compressor using structural modification of the accumulator[C]//International Compressor Engineering Conference, 1998: 1265.

[6] MAGLIE P, CARBINI R, WEIKERT S, et al. Efficient mechatronic evaluation of machine tool designs using model reduction[C]//Proceedings of the 12th Mechatronic Forum Biennial International Conference, 2010: 285-292.

Experimental modal analysis and model updating on outdoor unit of conditioner

WANG Bor-tsuen1, ZHANG Bo-jun1, CHEN Guan-zhou1, LIANG Hsiu-wei2

(1. Department of Mechanical Engineering, Pingtung University of Science and Technology, Pingtung 91201, Taiwan, China;2. Rechi Precision Co., LTD, Taoyuan 32849, Taiwan, China)

The correctness and reliability of finite element (FE) model is crucial to ensure structural simulation to fulfill its purpose. This work introduces the FE modeling techniques and model updating for outdoor unit of conditioner (OUC). The FE model of OUC must reveal the reality of structure, including geometry, material, interface and boundary conditions, and to be validated by experiments. The general approach is to perform experimental modal analysis (EMA) on the real structure and obtain structural modal parameters. At the meantime, the structural modal parameters can be numerically obtained from FE model. Model updating is to adjust the FE model such that both FE model and real structure are equivalent based on the agreement of structural modal characteristics. Results show the updated FE model of OUC can reveal good agreement with real structure in modal and frequency domain properties.

outdoor unit of conditioner(OUC); experimental modal analysis(EMA); finite element analysis; model updating

TB657.2

A

1000-3630(2019)-04-0433-07

10.16300/j.cnki.1000-3630.2019.04.013

2018-10-16;

2019-02-11

瑞智精密股份有限公司產(chǎn)學(xué)合作案(10511193)

王栢村(1960-), 男, 臺灣屏東人, 博士, 研究方向為結(jié)構(gòu)振動。

王栢村,E-mail: wangbt@mail.npust.edu.tw

猜你喜歡
模態(tài)有限元實驗
記一次有趣的實驗
做個怪怪長實驗
NO與NO2相互轉(zhuǎn)化實驗的改進
實踐十號上的19項實驗
太空探索(2016年5期)2016-07-12 15:17:55
國內(nèi)多模態(tài)教學(xué)研究回顧與展望
基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
磨削淬硬殘余應(yīng)力的有限元分析
由單個模態(tài)構(gòu)造對稱簡支梁的抗彎剛度
計算物理(2014年2期)2014-03-11 17:01:39
基于SolidWorks的吸嘴支撐臂有限元分析
箱形孔軋制的有限元模擬
上海金屬(2013年4期)2013-12-20 07:57:18
主站蜘蛛池模板: 亚洲第七页| 中国国产一级毛片| 91午夜福利在线观看| 国产精品一区不卡| 国产成人亚洲综合A∨在线播放| 人人看人人鲁狠狠高清| 在线色国产| h网站在线播放| 欧美三级不卡在线观看视频| 国产成人禁片在线观看| 99久久精品无码专区免费| 91久久国产热精品免费| 国产精品页| 精品无码日韩国产不卡av| 国产人妖视频一区在线观看| 综1合AV在线播放| 国产亚洲精品精品精品| 亚洲无码免费黄色网址| 中文字幕久久亚洲一区| 97视频精品全国免费观看| 久久情精品国产品免费| 国产资源免费观看| www欧美在线观看| 黄片在线永久| 亚洲国产欧美国产综合久久| 69综合网| 久久这里只有精品8| 亚洲最大情网站在线观看| 精品久久国产综合精麻豆| 亚洲精品综合一二三区在线| 自拍中文字幕| 国模极品一区二区三区| 在线播放国产一区| 中国国产A一级毛片| 午夜无码一区二区三区| 日韩精品一区二区深田咏美| 国产青青草视频| 天天躁夜夜躁狠狠躁图片| 国产精品爽爽va在线无码观看| 亚洲区视频在线观看| 久久久精品国产亚洲AV日韩| 色悠久久久久久久综合网伊人| 高清亚洲欧美在线看| 欧美国产视频| 国内黄色精品| 亚洲第一色视频| 中文字幕永久视频| 色综合综合网| 最近最新中文字幕在线第一页| 91黄视频在线观看| 韩日无码在线不卡| 国产十八禁在线观看免费| 亚洲成a人片| 日本免费一区视频| 午夜不卡福利| 国产精品毛片一区视频播| 国产欧美日韩va另类在线播放| 亚洲一级毛片在线观播放| 久青草国产高清在线视频| 性激烈欧美三级在线播放| 欧美福利在线播放| 国产超薄肉色丝袜网站| 亚洲人成网站18禁动漫无码| 色哟哟国产精品| 亚洲aaa视频| 久久精品亚洲专区| 制服丝袜亚洲| 最新国产麻豆aⅴ精品无| 国精品91人妻无码一区二区三区| 黄色免费在线网址| 亚洲第一成人在线| 色综合五月婷婷| 免费人成黄页在线观看国产| 国产jizz| 99草精品视频| 精品無碼一區在線觀看 | 成人午夜免费观看| 国产国产人成免费视频77777| 免费午夜无码18禁无码影院| 国产高颜值露脸在线观看| 91系列在线观看| 呦视频在线一区二区三区|