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

        ?

        船舶尾部艙室噪聲預(yù)報及控制分析

        2017-09-11 12:43:10桂洪斌
        艦船科學(xué)技術(shù) 2017年8期
        關(guān)鍵詞:船舶結(jié)構(gòu)模型

        馮 博,桂洪斌,楊 群

        (哈爾濱工業(yè)大學(xué)(威海) 船舶與海洋工程學(xué)院,山東 威海 264209)

        船舶尾部艙室噪聲預(yù)報及控制分析

        馮 博,桂洪斌,楊 群

        (哈爾濱工業(yè)大學(xué)(威海) 船舶與海洋工程學(xué)院,山東 威海 264209)

        基于VA one軟件平臺對某散貨船尾部試驗?zāi)P瓦M(jìn)行噪聲預(yù)報及控制效果分析。分析中采用的是統(tǒng)計能量分析法(SEA),在驗證計算方法正確的基礎(chǔ)上,首先采用新規(guī)范標(biāo)準(zhǔn)對10個主要艙室的噪聲水平進(jìn)行預(yù)報研究和比對,分析其中艙室噪聲超標(biāo)的可能原因。然后對該分析模型進(jìn)行不同控制技術(shù)的研究分析——吸聲技術(shù)、隔聲技術(shù)和阻尼減振技術(shù),并比較在不同位置敷設(shè)阻尼材料的降噪程度。研究表明:采用吸聲、隔聲和阻尼減振技術(shù)對降低船舶艙室噪聲有顯著效果,在激勵源艙室敷設(shè)阻尼材料,僅對非激勵源艙室降噪效果明顯,且約束阻尼要比自由阻尼結(jié)構(gòu)對噪聲控制效果更有優(yōu)勢。研究結(jié)論可以作為船舶艙室噪聲實際控制的參考。

        船舶艙室噪聲;統(tǒng)計能量分析;高頻噪聲;噪聲控制

        0 引 言

        國際海事組織International Marine Organization(IMO)海上安全委員會(MSC)通過幾次會議的討論,在第91次會議上正式通過了《船上噪聲等級規(guī)則》修訂草案[1],對之前的船舶噪聲規(guī)定決議要求強(qiáng)制執(zhí)行;并且規(guī)定在現(xiàn)有的IMO標(biāo)準(zhǔn)[2]基礎(chǔ)上再降低5 dB。該修訂草案已于2014年7月1日正式生效。IMO發(fā)布的這些相關(guān)強(qiáng)制性規(guī)定對我國的造船企業(yè)產(chǎn)生了嚴(yán)重的影響。

        目前船舶艙室噪聲數(shù)值預(yù)報方法比較常用的是有限元法、邊界元法和統(tǒng)計能量分析法[3]。趙德有[4]等人基于統(tǒng)計能量分析方法對船舶結(jié)構(gòu)進(jìn)行三維聲學(xué)建模,并研究討論阻尼結(jié)構(gòu)對船舶艙室噪聲的影響。吳衛(wèi)國等[5]人運用SEA法對某高速船舶噪聲問題進(jìn)行計算分析,討論舷外水對艙室噪聲傳播的影響,發(fā)現(xiàn)舷外水對高速船舶艙室降噪有10%以內(nèi)的貢獻(xiàn)。胡忠平等[6]為了證明有限元法的準(zhǔn)確性及其實用性,利用有限元法對某船舶的結(jié)構(gòu)噪聲進(jìn)行了分析預(yù)報,將得到的分析結(jié)果與實際測量數(shù)據(jù)進(jìn)行了比較。楊德慶等[7]人在海洋工程的噪聲預(yù)報設(shè)計工作中運用邊界元法,對FPSO上層建筑的艙室噪聲進(jìn)行了數(shù)值模擬和分析預(yù)報。楊群[8]利用統(tǒng)計能量分析法對某艙室模型高頻噪聲進(jìn)行了預(yù)報和控制研究。

        船舶是一種大型復(fù)雜結(jié)構(gòu)物,通常處于高頻激勵作用下,因此本文采用統(tǒng)計能量分析法(SEA)進(jìn)行研究。首先VA one軟件和702所船舶聲學(xué)設(shè)計評估軟件建立相同驗證模型,分別在702所軟件的經(jīng)驗公式平臺、2個軟件的SEA平臺預(yù)報計算,將結(jié)果對比驗證,驗證VA one可以預(yù)報船舶艙室,然后以工信部高技術(shù)船舶科研項目設(shè)計的尾部模型(該模型以某散貨船為母型船)為研究對象進(jìn)行預(yù)報研究,建立相應(yīng)模型并選取統(tǒng)計能量分析法計算所需的參數(shù),如模態(tài)數(shù)/模態(tài)密度、內(nèi)損耗因子、耦合損耗因子以及輸入功率等,添加激勵進(jìn)行預(yù)報分析。最后通過對尾部艙室噪聲采用吸聲、隔聲和阻尼減振技術(shù)研究控制方法的降噪效果,以及不同位置采用控制方法的降噪效果。

        1 SEA法的基本理論及數(shù)值模型驗證

        1.1 理論基礎(chǔ)

        統(tǒng)計能量分析法(Statistical Energy Analysis),也稱為SEA法,適用于解決復(fù)雜結(jié)構(gòu)處于寬頻帶高頻動力下的相應(yīng)響應(yīng)難題。在研究振動問題的過程中,當(dāng)某個振動系統(tǒng)處于穩(wěn)態(tài)情況時其輸入功率Pin必然會與損耗功率Pd相平衡。當(dāng)分析的系統(tǒng)由3個及以上的子系統(tǒng)構(gòu)建而成,如果用Pi,in來表示第i個子系統(tǒng)的輸入功率,那么可以得到對應(yīng)系統(tǒng)的能量平衡基本關(guān)系式[9],即

        根據(jù)互易原理可以得到:

        將得到的式(2)代入式(1),通過整理分析可以得到需要的統(tǒng)計能量分析法基本分析方程:

        式中:

        1.2 數(shù)值模型及驗證

        結(jié)合文獻(xiàn)[10]對船舶SEA模型設(shè)計,本文通過在VA one軟件和702所船舶聲學(xué)設(shè)計評估軟件建立相同驗證模型,分別在702所軟件的經(jīng)驗公式平臺、2個軟件的SEA平臺對比驗證。計算其主機(jī)艙及集控室的噪聲級。該SEA模型的主要結(jié)構(gòu)尺寸以及激勵源設(shè)備參數(shù)如下:

        1)主機(jī)艙(噪聲源艙室)

        長寬高分別為10 m,10 m和10 m,艙壁為10 mm鋼板。

        2)集控室(接受艙室)

        長寬高分別為10 m,5 m和10 m,艙壁為10 mm鋼板。

        模型包括艙壁模型和聲腔模型。選取板壁材料為鋼材,賦予相應(yīng)數(shù)值,并查詢數(shù)據(jù)庫,配置鋼的自損耗因子。計算頻率范圍選為倍頻程63 ~ 8 000 Hz。2個軟件模型分別如圖1和圖2所示。

        主機(jī)艙內(nèi)的噪聲源主要由柴油機(jī)組產(chǎn)生的空氣噪聲以及機(jī)艙壁輻射噪聲所構(gòu)成。根據(jù)文獻(xiàn)[11],柴油機(jī)噪聲總聲功率計算方法及文獻(xiàn)中高速柴油機(jī)噪聲頻譜特性,計算出主機(jī)艙內(nèi)噪聲激勵源在倍頻程的噪聲強(qiáng)度如表1所示,之后將其施加到驗證模型上。

        表 1 主機(jī)艙噪聲源強(qiáng)度Tab. 1 Noise source strength of machinery space

        利用702所軟件和VA one軟件運用SEA法分別進(jìn)行預(yù)報計算,最后在經(jīng)驗公式平臺分析,計算結(jié)果如圖3和圖4所示。將得到的3個預(yù)報結(jié)果進(jìn)行對比驗證,結(jié)果如表2所示。

        由圖表可知,本文所用統(tǒng)計能量分析法和VA one軟件可以較準(zhǔn)確地對所研究結(jié)構(gòu)進(jìn)行計算預(yù)報。

        表 2 不同方式得到的預(yù)報平均值比較Tab. 2 Comparison of predicted mean values obtained by different methods

        2 某散貨船艙室噪聲預(yù)報分析

        2.1 目標(biāo)船舶SEA模型的建立

        根據(jù)目標(biāo)船舶的主要參數(shù)(見表3),按照SEA法的基本原理進(jìn)行聲學(xué)建模,按照相應(yīng)的要求和規(guī)范建立相應(yīng)的目標(biāo)船舶的結(jié)構(gòu)模型和聲腔模型,在運用SEA法建模時,對目標(biāo)船舶的模型中加強(qiáng)筋、支柱、縱骨、桁架等結(jié)構(gòu)進(jìn)行適當(dāng)?shù)暮喕?,外板、甲板平臺和艙壁采用加筋板子系統(tǒng),并把與聲場流體介質(zhì)耦合的結(jié)構(gòu)子系統(tǒng)與聲場相連。

        本文主要對該散貨船的船尾部分進(jìn)行預(yù)報研究,該分段的長度為19.6 m。分別對模型船體與聲場介質(zhì)(空氣)賦予屬性,建立完成的目標(biāo)船舶SEA模型的縮放模型和耦合模型分別如圖5和圖6所示。

        表 3 目標(biāo)船舶的主要參數(shù)Tab. 3 Main parameters of target ship

        2.2 預(yù)報結(jié)果及分析

        確定了SEA參數(shù)以后,選取31.5~8 000 Hz之間的倍頻程作為分析頻帶,通過VA one軟件計算可得各聲腔子系統(tǒng)在各個頻率下的聲壓級水平,平均聲壓級水平云圖如圖7所示,預(yù)報結(jié)果如圖8所示。

        通過比較,發(fā)現(xiàn)機(jī)艙、機(jī)艙工作區(qū)以及廚房等與機(jī)艙相近的艙室超過IMO發(fā)布的規(guī)范,分析原因有以下幾點:

        1)機(jī)艙區(qū)域噪聲值超標(biāo)由于在用統(tǒng)計能量法建模時,沒有考慮船底龍骨、縱骨、外板加強(qiáng)筋等結(jié)構(gòu),使得模型在質(zhì)量和剛度上較有限元模型小了很多,從而造成船體結(jié)構(gòu)振動劇烈,以至增大了聲輻射;

        2)數(shù)值模擬模型并沒有考慮船舶的舾裝結(jié)構(gòu),即艙室壁吸聲等材料;

        3)子系統(tǒng)的內(nèi)損耗因子以及系統(tǒng)之間的耦合損耗因子這2個參數(shù)難以確定,造成了預(yù)報結(jié)果的誤差。

        3 艙室噪聲控制技術(shù)研究

        目前比較有效及常見的船舶艙室噪聲控制技術(shù)有吸聲技術(shù)、隔聲技術(shù)、隔振技術(shù)以及阻尼減振技術(shù)四類。本文選取吸聲技術(shù)、隔聲技術(shù)以及阻尼減振技術(shù)進(jìn)行研究。

        3.1 吸聲技術(shù)的應(yīng)用

        吸聲技術(shù)的原理是使聲波透射到材料內(nèi)部中去,利用材料與空氣的相互摩擦作用,從而使聲能以熱能的形式耗散掉[12]。前文利用VA one軟件建立的船舶模型其10個主要艙室,分別在每個艙室的艙壁上附加一層密度50 kg/m3、彈性模量7.86×106 Pa、泊松比是0.450 2、厚度為0.002 m的氯丁阻燃橡膠CR322[9],將船舶各主要艙室的平均噪聲水平在附加吸聲材料前后的預(yù)報結(jié)果進(jìn)行對比,如表4所示??梢园l(fā)現(xiàn)在船舶艙室的內(nèi)部艙壁附加一定厚度的吸聲材料,對降低船舶艙室噪聲聲壓級有顯著效果。

        表 4 各主要艙室平均預(yù)報結(jié)果 dB(A)Tab. 4 Average prediction results of main cabins before and after the addition of sound absorption materials dB(A)

        3.2 隔聲技術(shù)的應(yīng)用

        隔聲是在聲波的擴(kuò)散途徑中,利用結(jié)構(gòu)把激勵源與接受源隔離起來的措施。通常的做法是把聲源、或者接受源封閉在一個有限的空間內(nèi),從而產(chǎn)生一個封閉的空間,以達(dá)到降低噪聲的目的。本文主要分析討論的是對空氣噪聲的隔絕。所以對于集控室(機(jī)艙工作區(qū))艙壁板選取雙層隔聲結(jié)構(gòu),分別把雙層隔聲材料粘貼在集控室艙壁外側(cè)和內(nèi)側(cè),比較2種不同粘貼位置對集控室的噪聲水平的作用效果進(jìn)行分析,可以得到結(jié)果如圖9所示。

        可以發(fā)現(xiàn),在集控室外側(cè)粘貼隔聲材料比在集控室內(nèi)側(cè)粘貼隔聲材料會產(chǎn)生更好的隔聲效果,因為噪聲源主要是機(jī)艙內(nèi)的空氣噪聲,所以在集控室艙壁外側(cè)粘貼隔聲材料更有利于減少空氣噪聲的影響,從而使集控室噪聲水平得以降低。

        3.3 阻尼減振技術(shù)的應(yīng)用

        阻尼減振技術(shù)的原理即為敷設(shè)某些特定材料或結(jié)構(gòu)來提高結(jié)構(gòu)的損耗因數(shù)(或阻尼系數(shù)),通過敷設(shè)阻尼層的變形來損耗系統(tǒng)結(jié)構(gòu)的能量,即使固體機(jī)械振動能量轉(zhuǎn)換成熱能耗散,進(jìn)而實現(xiàn)降低噪聲水平的目標(biāo)。

        目前基板敷設(shè)阻尼層的類型主要有自由阻尼和約束阻尼2種。自由阻尼是只在基板上敷設(shè)一層阻尼材料,如圖10(a)所示;約束阻尼是在自由阻尼方法的基礎(chǔ)上,在敷設(shè)的阻尼材料外再敷設(shè)一層新的約束材料,如圖10(b)所示。

        本文自由阻尼材料選用VA one材料庫中的聚酯泡沫,其主要的物理參數(shù)是:ρ=30 kg/m3;泊松比v=0.35;彈性模量5.4×105 Pa;阻尼材料損耗因子為0.1;阻尼層厚度0.1 m。約束阻尼材料物理屬性:阻尼層布置在結(jié)構(gòu)的中間位置,約束層分布在阻尼層的上方,厚度一般選用基材船舶鋼板厚度的一半;阻尼層選用聚酯泡沫材料,具體材料參數(shù)與自由阻尼層的相同。

        在機(jī)艙位置同時施加相同的結(jié)構(gòu)噪聲激勵和空氣噪聲激勵,然后先后在機(jī)艙敷設(shè)自由阻尼材料和約束阻尼材料,比較敷設(shè)前后居住艙室內(nèi)噪聲水平的變化情況,其變化曲線如圖11所示。

        由圖可以發(fā)現(xiàn),在機(jī)艙位置粘貼相應(yīng)的阻尼材料無論是自由阻尼材料還是約束阻尼材料,都可以在一定程度上降低居住艙室的噪聲水平,而且約束阻尼結(jié)構(gòu)比自由阻尼結(jié)構(gòu)的控制效果更有優(yōu)勢。

        選用VA one材料庫中阻尼系數(shù)較高的粘彈性阻尼材料橡膠,其基本的物理屬性為:ρ=1 100 kg/m3;泊松比v=0.48;剪切模量7.7×108 Pa;拉伸模量2.3×109 Pa;阻尼材料損耗因子0.1;阻尼層厚度0.1 m。將其以自由阻尼形式分別敷設(shè)在機(jī)艙位置,居住艙室的下甲板和后壁板,然后對這2個艙室的敷設(shè)阻尼材料前后的噪聲水平進(jìn)行比較研究。2個艙室位置的噪聲變化情況分別如圖12和圖13所示。

        分析比較,可以發(fā)現(xiàn)在機(jī)艙位置粘貼阻尼材料對機(jī)艙本身的噪聲水平基本沒有任何的作用效果;在居住艙室敷設(shè)阻尼材料能夠明顯降低接受聲源艙室的聲壓級水平,可以明顯降低居住艙室能量的輸入,從而起到很好的減噪效果。

        4 結(jié) 語

        本文基于統(tǒng)計能量分析法對船舶艙室的噪聲進(jìn)行了預(yù)報和控制方法的分析研究,主要得出如下結(jié)論:

        1)通過3種方式對比驗證,發(fā)現(xiàn)預(yù)報結(jié)果三者之間的誤差在3 dB以內(nèi),間接證明了本文所用VA one軟件可以較準(zhǔn)確地對所研究結(jié)構(gòu)進(jìn)行計算預(yù)報。

        2)通過應(yīng)用吸聲技術(shù),以及對比VA one軟件測得的在附加吸聲材料前后各主要艙室預(yù)報結(jié)果,可以發(fā)現(xiàn)在船舶艙室的內(nèi)部艙壁附加一定厚度的吸聲材料,對降低船舶艙室噪聲聲壓級有顯著效果。

        3)在集控室艙壁外側(cè)粘貼隔聲材料更有利于減少空氣噪聲的影響,從而使集控室噪聲水平得以降低。

        4)在激勵源所在艙室(機(jī)艙附近)敷設(shè)阻尼材料,機(jī)艙等艙室噪聲降噪效果不明顯,但是其對其他艙室有很明顯的降噪效果,且約束阻尼結(jié)構(gòu)比自由阻尼結(jié)構(gòu)對非激勵源艙室的噪聲控制效果更有優(yōu)勢。

        5)在非激勵源艙室敷設(shè)一定厚度的阻尼材料,對非激勵源艙室有很明顯的降噪效果。

        [ 1 ]張輝, 張新波. 《船上噪聲等級規(guī)則》修訂案解析及應(yīng)對措施[J]. 船舶標(biāo)準(zhǔn)化工程師, 2013, 46(01): 28–29. ZHANG Hui, ZHANG Xin-bo. The influence on new ships of code on noise levels on-board ships[J]. Ship Standardization Engineer, 2013, 46(01): 28–29.

        [ 2 ]International Marine Organization (IMO)[J]. Code on noise levels on board ships (IMO A468). 1981: 1.

        [ 3 ]李卓亮, 季振林. 車客渡船艙室噪聲預(yù)報與控制[J]. 中國艦船研究, 2014, 9(5): 93–98, 120. LI Zhuo-liang, JI Zhen-lin. Cabin noise prediction and control of car-passenger ferry[J]. Chinese Journal of Ship Research, 2014, 9(5): 93–98, 120.

        [ 4 ]于大鵬, 趙德有, 汪玉. 船舶聲學(xué)建模和阻尼結(jié)構(gòu)對艙室噪聲影響研究[J]. 船舶力學(xué), 2010, 14(5): 539–548. YU Da-peng, ZHAO De-you, WANG Yu. Influence of damped material and the ship model of acoustic on the ship cabin noise [J]. Journal of Ship Mechanics, 2010, 14(5): 539–548.

        [ 5 ]金學(xué)之, 吳衛(wèi)國. 高速船舷外水對艙室噪聲傳播影響[J]. 船海工程, 2009(6): 46–48. JIN Xue-zhi, WU Wei-guo. Study on the out-shipboard fluid affecting the cabin noise propagation for high speed ship[J]. Ship & Ocean Engineering, 2009(6): 46–48.

        [ 6 ]胡忠平. 船體結(jié)構(gòu)噪聲的有限元計算分析[J]. 船舶, 2008, 19(3): 30–33. HU Zhong-ping. Finite element analysis on hull structural noise[J]. Ship & Boat, 2008, 19(3): 30–33.

        [ 7 ]楊德慶, 金咸定, 周永青. 浮式生產(chǎn)儲油船上層建筑振動噪聲數(shù)值預(yù)報[J]. 上海造船, 2003(2): 15–19. YANG De-qing, JIN Xian-ding, ZHOU Yong-qing, et al. Numerical prediction of vibro-acoustics in superstructure of FPSO[J]. Shanghai Shipbuilding, 2003(2): 15–19.

        [ 8 ]楊群. 某艙室模型高頻噪聲預(yù)報及控制技術(shù)研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2015. YANG Qun. High frequency noise prediction and control technology study of certain cabin model[D]. Harbin: Harbin Institute of Technology, 2015.

        [ 9 ]于大鵬. 應(yīng)用統(tǒng)計能量分析方法預(yù)報船舶艙室噪聲[D]. 大連: 大連理工大學(xué)碩士論文, 2007. YU Da-peng. Study on the prediction of shipyard cabin noise with statistical energy analysis[D]. Dalian: Dalian University of Technology, 2007.

        [10]孫增華, 王丹, 于洋, 等. 船舶艙室噪聲的工程預(yù)報方法[J].船舶工程, 2014, 1. SUN Zeng-hua, WANG Dan, YU Yang, et al. Engineering prediction approach for shipboard cabin noise[J]. Ship Engineering, 2014, 1.

        [11]MAGRAB E B. Environmental noise control[M]. John Wiley & Sons, 1975.

        [12]俞微, 沈榮瀛. 浮筏隔振系統(tǒng)功率流的實驗研究[J]. 振動、測試與診斷, 1996, 16(2): 44–48. YU Wei, SHEN Rong-ying. Experimental study of power flow of a raft vibration-isolating system[J]. Journal of Vibration. Measurement & Diagnosis, 1996, 16(2): 44–48.

        Research on prediction method and control technology of shipboard aft cabin noise

        FENG Bo, GUI Hong-bin, YANG Qun
        (School of Naval Architecture and Ocean Engineering, Harbin Institute of Technology, Weihai 264209, China)

        In this paper, according to the VA one software platform and statistical energy analysis (SEA), a research was studied the noise prediction and control technology of a bulk carrier. By using the statistical energy analysis method in the research process (SEA), based on the verification of correct calculation method. Firstly, according to the new standard, the noise level of 10 main cabin were predicted, compared, and analyzed the reason of the cabin noise may exceed the standard. At the end of this model was studied in three different control methods: the sound absorption, sound insulation and vibration damping technology, and compared the degree of noise reduction in different positions for damping material. Research shows that: the three control methods have significant effect to reduce cabin noise. The conclusion of this paper can as the actual control of the ship cabin noise.

        shipboard cabin noise;statistical energy analysis;high-frequency noise;noise control

        U661.44

        A

        1672 – 7649(2017)08 – 0022 – 06

        10.3404/j.issn.1672 – 7649.2017.08.005

        2016 – 09 – 28;

        2016 – 11 – 04

        馮博(1992 – ),男,碩士研究生,研究方向為船舶艙室噪聲預(yù)報及控制技術(shù)。

        猜你喜歡
        船舶結(jié)構(gòu)模型
        一半模型
        計算流體力學(xué)在船舶操縱運動仿真中的應(yīng)用
        《船舶》2022 年度征訂啟事
        船舶(2021年4期)2021-09-07 17:32:22
        《形而上學(xué)》△卷的結(jié)構(gòu)和位置
        重要模型『一線三等角』
        重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
        船舶!請加速
        論結(jié)構(gòu)
        中華詩詞(2019年7期)2019-11-25 01:43:04
        論《日出》的結(jié)構(gòu)
        3D打印中的模型分割與打包
        国产女人高潮的av毛片| 首页 综合国产 亚洲 丝袜| 无码人妻丰满熟妇区毛片| 国产精品99久久久精品免费观看| 亚洲午夜精品国产一区二区三区| 视频在线观看国产自拍| 精品免费久久久久久久| 成人无码午夜在线观看| 黄色av三级在线免费观看| 三级日韩视频在线观看| 亚洲成av人在线播放无码| 狠狠色婷婷久久一区二区| 久久久久久国产福利网站| 日本一区二区三区区视频| 国产色xx群视频射精| 国产成人免费a在线视频| 日韩精品视频免费福利在线观看 | 国产日产一区二区三区四区五区 | 18无码粉嫩小泬无套在线观看| 啪啪网站免费观看| 精品国产中文久久久免费| 亚洲熟妇av一区| 久久精品国产9久久综合| 免费国产在线精品三区| 日本系列中文字幕99| 国产精品亚洲一区二区在线观看 | 国产亚洲人成a在线v网站| 男人天堂av在线成人av| 久久精品国产在热亚洲不卡| av免费不卡国产观看| 欧美日韩亚洲精品瑜伽裤| 大量老肥熟女老女人自拍| 亚洲中文久久精品字幕| 国产情侣久久久久aⅴ免费| 亚洲性无码一区二区三区| 国产精品亚洲五月天高清| 日韩av在线不卡一二三区| 国产激情视频免费在线观看| 无码精品人妻一区二区三区人妻斩 | 国产亚洲欧美在线| 97超碰中文字幕久久|