• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    靜壓下帶金屬孔板聲學(xué)覆蓋層的耐壓與吸聲特性

    2024-12-31 00:00:00賈薪宇靳國永葉天貴閆燕
    振動(dòng)工程學(xué)報(bào) 2024年7期
    關(guān)鍵詞:靜水壓力

    摘要: 為改善靜水壓力作用下傳統(tǒng)聲學(xué)覆蓋層的吸聲性能,提出了將金屬孔板插入傳統(tǒng)聲學(xué)覆蓋層中的水下聲學(xué)覆蓋層結(jié)構(gòu)。以靜力學(xué)方法研究了該聲學(xué)覆蓋層結(jié)構(gòu)在不同靜水壓力下的形變情況。通過建立靜水壓力作用下的聲學(xué)有限元方程,分析了該聲學(xué)結(jié)構(gòu)在不同靜水壓力下的吸聲效果。與傳統(tǒng)聲學(xué)覆蓋層相比,在靜水壓力為0~6 MPa時(shí),該聲學(xué)結(jié)構(gòu)在中、高頻段獲得了更好的寬頻吸聲效果。此外,還討論孔板厚度、孔板材料及孔板孔隙率對(duì)該聲學(xué)結(jié)構(gòu)吸聲性能的影響。研究表明設(shè)計(jì)含金屬孔板聲學(xué)覆蓋層可有效改善靜水壓力條件下水下聲學(xué)結(jié)構(gòu)的吸聲性能。

    關(guān)鍵詞: 聲學(xué)覆蓋層; 金屬孔板; 靜水壓力; 形變; 吸聲性能

    中圖分類號(hào): TB564;TB535""" 文獻(xiàn)標(biāo)志碼: A""" 文章編號(hào): 1004-4523(2024)07-1211-10

    DOI:10.16385/j.cnki.issn.1004-4523.2024.07.013

    收稿日期: 2022-10-03; 修訂日期: 2022-11-14

    基金項(xiàng)目:"國家自然科學(xué)基金資助項(xiàng)目(52225109,52241101,52271309);黑龍江省自然科學(xué)基金資助項(xiàng)目(YQ2022E104);哈爾濱工程大學(xué)博士科創(chuàng)基金資助項(xiàng)目(3072023GIP0302)。

    引 言

    聲學(xué)覆蓋層是敷設(shè)在水下航行器外表面用以實(shí)現(xiàn)聲隱身性能的關(guān)鍵部件。隨著反潛探測(cè)技術(shù)的不斷發(fā)展[1],為進(jìn)一步提高水下航行器的生存能力,提高和改善聲學(xué)覆蓋層吸聲性能的研究工作得到了蓬勃發(fā)展,涌現(xiàn)了許多新型的結(jié)構(gòu)設(shè)計(jì)[2?6]與研究方法[7?10],在較寬的頻率范圍能夠獲得良好的吸聲效果[11?13]。

    目前,對(duì)聲學(xué)覆蓋層的研究工作多數(shù)是在常壓下進(jìn)行的,但水下航行器是在深水中工作,必然受高靜水壓力的作用。因此,在靜水壓力下聲學(xué)覆蓋層的設(shè)計(jì)與研究得到了許多學(xué)者的關(guān)注。姜聞文等[14]利用有限元方法研究了在不同靜水壓力下,聲學(xué)覆蓋層內(nèi)不同腔體的變形情況及聲學(xué)性能的變化規(guī)律,其研究表明,吸聲頻帶隨靜水壓力增加而向高頻方向移動(dòng),其研究工作指明了空腔有效容積與吸聲曲線頻帶的對(duì)應(yīng)關(guān)系。姚熊亮等[15]基于非均勻波導(dǎo)理論推導(dǎo)出多層結(jié)構(gòu)的吸聲系數(shù),研究了靜水壓力與腔型對(duì)聲學(xué)覆蓋層吸聲系數(shù)的影響,其研究成果進(jìn)一步指出了靜水壓力對(duì)聲學(xué)覆蓋層吸聲性能具有重要影響,在實(shí)際應(yīng)用中不可忽略。Humphrey等[16]利用耐壓水箱實(shí)驗(yàn)方法,模擬了深度達(dá)700 m的深海環(huán)境,同時(shí)測(cè)試了聲學(xué)覆蓋層在靜水壓力為2.8 MPa時(shí)的聲反射和聲透射性能,為實(shí)驗(yàn)?zāi)M研究高靜水壓力下聲學(xué)覆蓋層的聲學(xué)性能提供了思路。張沖等[17]利用移動(dòng)網(wǎng)格建模方法,在結(jié)構(gòu)受壓變形的基礎(chǔ)上直接計(jì)算分析了帶有球腔的聲學(xué)覆蓋層的吸聲效果,并對(duì)比了腔壓對(duì)腔體形變和吸聲性能的影響,該研究成果改進(jìn)了以往先求結(jié)構(gòu)變形再對(duì)結(jié)構(gòu)重新建模的計(jì)算方法,簡(jiǎn)化了建模過程,為后續(xù)工作提供了簡(jiǎn)便方法。Thieury等[18]通過數(shù)值模型研究了靜水壓力對(duì)橡膠聲學(xué)覆蓋層聲反射的影響,并考慮了靜水壓力引起的空腔體積變化,其研究工作指明了在靜水壓力作用下,空腔大小、空腔距離以及彈性基體的流變性對(duì)結(jié)構(gòu)聲學(xué)性能的影響。Gu等[19]采用實(shí)驗(yàn)的方法,對(duì)帶有局域共振單元的聲學(xué)覆蓋層的性能進(jìn)行了分析,證明了帶有局域共振單元的吸聲結(jié)構(gòu)在靜水壓力下可以獲得更好的低頻吸聲效果,為耐壓型聲學(xué)覆蓋層的設(shè)計(jì)提供了借鑒思路。董文凱等[20]在研究靜壓力下水下聲學(xué)覆蓋層吸聲性能的同時(shí)考慮了腔壓對(duì)形變的影響,仿真結(jié)果顯示靜水壓力越大,腔壓對(duì)吸聲體聲學(xué)性能的作用越明顯。Zhang等[21]基于唯象理論研究了靜水壓力下含有亞波長(zhǎng)壓電陣列的半主動(dòng)聲學(xué)覆蓋層的幾何變形和聲學(xué)性能,并探討了相應(yīng)吸聲效果的變化規(guī)律,其研究工作既以提高聲學(xué)覆蓋層吸聲性能為研究目的,又同時(shí)考慮靜水壓力對(duì)吸聲效果的影響。對(duì)于靜水壓力下材料參數(shù)的變化,尚未有成熟的理論與經(jīng)驗(yàn)公式,仍主要以實(shí)驗(yàn)研究為主[22?23]。因此,目前對(duì)不同靜水壓力下聲學(xué)覆蓋層的研究,主要以結(jié)構(gòu)的形變對(duì)其聲學(xué)性能的影響為主。

    基于以上討論,本文以靜水壓力條件下改善聲學(xué)覆蓋層吸聲特性為研究目標(biāo),設(shè)計(jì)了一種帶有金屬孔板的水下聲學(xué)覆蓋層結(jié)構(gòu)。討論在不同靜水壓力下聲學(xué)覆蓋層的形變量,研究該聲學(xué)結(jié)構(gòu)的耐壓性能?;谟邢拊椒ń㈧o水壓力作用下結(jié)構(gòu)的聲學(xué)有限元方程,并通過已有文獻(xiàn)驗(yàn)證所建立模型的有效性和正確性。基于所建立的聲學(xué)模型,開展在不同靜水壓力下所設(shè)計(jì)的聲學(xué)覆蓋層吸聲特性的研究。此外,通過計(jì)算對(duì)比分析了孔板材料、孔板厚度以及孔板孔隙率等因素對(duì)所設(shè)計(jì)聲學(xué)結(jié)構(gòu)吸聲性能的影響規(guī)律。

    1 模型建立與計(jì)算方法

    1.1 模型建立

    本文提出的聲學(xué)復(fù)合結(jié)構(gòu)主要由黏彈性基體層、金屬孔板以及鋼板背襯組成。其中,基體層中分布著沿x和y方向呈周期性排列的球型空腔結(jié)構(gòu),在此結(jié)構(gòu)的基礎(chǔ)上,按一定規(guī)律插入沿x方向呈周期分布、y方向無限大的金屬孔板,具體形式如圖1所示。根據(jù)結(jié)構(gòu)的設(shè)計(jì)特點(diǎn),周期結(jié)構(gòu)的晶格常數(shù)a為30 mm,在黏彈性基體層后連接鋼背襯以模擬水下航行器實(shí)際的應(yīng)用條件,背襯后為空氣。圖2為聲學(xué)覆蓋層單一周期單元結(jié)構(gòu)及其截面示意圖,其中基體層厚度t1 = 50 mm,圓形空腔直徑d =24 mm。鋼板背襯厚度ts = 5 mm。胞元內(nèi)插入的金屬孔板尺寸b × c=18 mm × 30 mm,厚度tb為1 mm,在金屬板上呈均勻分布的圓柱孔孔徑為1 mm,孔隙率σ為16%?;w層前端為半無限大流體域,流體介質(zhì)假設(shè)為海水。入射平面波垂直于x?y平面,沿z正方向從半無限流體域入射。

    1.2 計(jì)算理論

    1.2.1 聲學(xué)有限元方程

    本文利用有限元法對(duì)聲學(xué)結(jié)構(gòu)的吸聲性能進(jìn)行分析,假設(shè)流體域?yàn)槔硐肓黧w,滿足均勻、無黏滯吸收、靜態(tài)流速為零、聲波小振幅等理想條件,聲波的波動(dòng)方程可以寫成如下形式:

    (1)

    式中 p為聲壓;c為流體介質(zhì)中的聲速。

    由于復(fù)合聲學(xué)結(jié)構(gòu)內(nèi)部呈現(xiàn)周期性分布,在建立有限元模型時(shí)只需對(duì)其中一個(gè)周期單元進(jìn)行分析,即可滿足計(jì)算要求。通過一系列數(shù)學(xué)推導(dǎo),基于有限元方法,獲得有限元單元的質(zhì)量矩陣、剛度矩陣、流?固耦合矩陣,通過求和運(yùn)算,集總成所研究周期單元的流體域和結(jié)構(gòu)域整體剛度矩陣、質(zhì)量矩陣以及流?固耦合矩陣,進(jìn)而得到流體域離散方程和結(jié)構(gòu)中的控制方程分別為[24]:

    (2)

    (3)

    式中 Mf和Kf分別為選取周期單元流體域的質(zhì)量矩陣和剛度矩陣;Ms和Ks分別為選取周期單元結(jié)構(gòu)域中的質(zhì)量矩陣和剛度矩陣;R為選取周期單元的流?固耦合矩陣;u為結(jié)構(gòu)域中節(jié)點(diǎn)位移的向量矩陣;p為流體域中節(jié)點(diǎn)的聲壓向量矩陣;Ψf為作用在流體域界面的聲壓梯度;Fs為施加在結(jié)構(gòu)域的激勵(lì);為流體密度。

    根據(jù)公式(2)和(3),周期單元的流?固耦合方程可以表示為:

    (4)

    式中 表示流體域邊界法向上的節(jié)點(diǎn)壓力載荷向量矩陣。Ff表示的是施加在流體邊界的等效載荷矩陣;為角頻率。

    由于聲學(xué)覆蓋層中的圓形空腔與金屬孔板結(jié)構(gòu)呈周期性分布,流體域也沿著結(jié)構(gòu)周期方向無限延伸,且流體介質(zhì)為均勻的理想流體,因此周期單元中結(jié)構(gòu)域中節(jié)點(diǎn)的位移與流體域中節(jié)點(diǎn)的聲壓均滿足Bloch方程,即:

    (5)

    式中 代表分布函數(shù),可表示為某節(jié)點(diǎn)的位移或聲壓;dx和dy分別表示x和y方向上的晶格單元距離。

    反射和透射聲壓可以用級(jí)數(shù)之和的形式表示為:

    (6)

    (7)

    式中 kx,ky,kz分別為x,y,z方向的波數(shù);Rnm和Tnm分別為反射聲波和透射聲波的m, n階系數(shù)。

    通過流?固耦合方程(4)以及周期邊界條件(5),可以求出入射面上的反射聲壓與透射面上的透射聲壓,進(jìn)而求解出通過該吸聲結(jié)構(gòu)的平面波反射系數(shù)R及透射系數(shù)T,表達(dá)式如下:

    (8)

    式中 N為反射邊界或透射邊界的有限元節(jié)點(diǎn)數(shù)量,因此結(jié)構(gòu)的吸聲系數(shù)。實(shí)際上,鋼的特性阻抗與空氣相差極大,導(dǎo)致二者阻抗失配,透射系數(shù)幾乎為零,所以最終計(jì)算出通過該結(jié)構(gòu)的吸聲系數(shù)。

    1.2.2 靜水壓下的有限元方程

    目前,普遍采用二次建模法來研究水下聲學(xué)覆蓋層受靜水壓力作用下的聲學(xué)性能。即先求解結(jié)構(gòu)受靜水壓作用后的變形量,再根據(jù)變形量重新建立規(guī)則的幾何模型及網(wǎng)格劃分,此方法過程繁瑣且后續(xù)的計(jì)算誤差較大。因此,本文利用線性化理論,將結(jié)構(gòu)受靜水壓力作用產(chǎn)生的勢(shì)能直接加載到聲學(xué)有限元方程中,形成考慮靜水壓力作用的流?固耦合有限元方程。

    考慮靜水壓力作用后,結(jié)構(gòu)的振動(dòng)控制方程可以寫成:

    (9)

    式中 K0為與作用于聲波入射面A0上的靜水壓力P0相關(guān)的剛度矩陣;Kr0為與靜水壓力作用下結(jié)構(gòu)域Vs的靜態(tài)預(yù)應(yīng)力Pr0相關(guān)的剛度矩陣,其結(jié)構(gòu)形式如下[25]:

    (10)

    (11)

    式中 Np為面單元的形函數(shù);Nδ為有限單元的位移插值形函數(shù)。結(jié)合流體域的有限元方程式(2),靜水壓力作用下的周期單元的流?固耦合方程可以表示為:

    (12)

    通過求解方程(12),可以得到流體域的節(jié)點(diǎn)聲壓,進(jìn)而根據(jù)前文分析的方法,求解出聲學(xué)結(jié)構(gòu)在靜水壓力條件下的吸聲系數(shù)。

    1.3 數(shù)值計(jì)算模型準(zhǔn)確性驗(yàn)證

    為了驗(yàn)證本文有限元計(jì)算模型的正確性與適用性,利用文本計(jì)算方法建立文獻(xiàn)[26]中的水下聲學(xué)覆蓋層模型,并與文獻(xiàn)中對(duì)聲學(xué)覆蓋層吸聲系數(shù)的計(jì)算結(jié)果進(jìn)行比對(duì)。從圖3的對(duì)比結(jié)果中可以看出,利用本文計(jì)算方法得到的聲學(xué)覆蓋層吸聲系數(shù)計(jì)算結(jié)果與文獻(xiàn)[26]中的數(shù)值仿真結(jié)果基本一致。而改變聲學(xué)覆蓋層內(nèi)部的周期結(jié)構(gòu)形式,相當(dāng)于只改變了結(jié)構(gòu)內(nèi)部不同材料之間的耦合邊界條件,并沒有改變聲?固耦合邊界。在利用有限元計(jì)算時(shí),只是代入方程的條件參數(shù)不同,而未改變方程形式。文獻(xiàn)[26]中的聲學(xué)覆蓋層的結(jié)構(gòu)同樣是在黏彈性基體層具有空腔與金屬類填充物,與文本結(jié)構(gòu)的邊界條件一致。因此,可以驗(yàn)證本文所用的聲學(xué)計(jì)算方法是正確的,并可適用于內(nèi)部結(jié)構(gòu)更為復(fù)雜的聲學(xué)覆蓋層計(jì)算中。此外,在有限元單元網(wǎng)格劃分過程中,網(wǎng)格單元的最小尺寸遵循小于材料中剪切波波長(zhǎng)1/4的設(shè)計(jì)原則[24],以滿足計(jì)算精度的要求。

    為了進(jìn)一步說明本文建立的在靜水壓力作用下的聲學(xué)有限元方程的準(zhǔn)確性,利用本文計(jì)算方法與文獻(xiàn)[20]中聲學(xué)覆蓋層在1.5 MPa靜水壓力下的吸聲系數(shù)作了對(duì)比,如圖4所示。從對(duì)比結(jié)果中可以得出,利用本文計(jì)算方法得到聲學(xué)覆蓋層在1.5 MPa靜水壓力下的吸聲系數(shù)計(jì)算結(jié)果與文獻(xiàn)[20]中解析法求得結(jié)果較吻合,進(jìn)一步說明了本文所建立的數(shù)值求解方法的正確性。

    2 不同靜壓下聲學(xué)覆蓋層的形變與吸聲性能

    2.1 不同靜壓下聲學(xué)覆蓋層的形變分析

    目前,水下航行器聲學(xué)覆蓋層主要以含有空腔結(jié)構(gòu)的橡膠材料為主,在大潛深條件下,水下航行器在幾百米的水下通常要承受幾兆帕的靜水壓力,這將使聲學(xué)覆蓋層產(chǎn)生較大變形,特別是內(nèi)部空腔結(jié)構(gòu)。為了更好地分析靜水壓力對(duì)聲學(xué)覆蓋層的影響,本文通過靜力學(xué)分析方法在一個(gè)周期單元的聲壓入射面施加均布載荷來模擬聲學(xué)覆蓋層在不同靜壓下的形變情況。因x,y方向?yàn)榧僭O(shè)的無限大平面,周期邊界處指定x和y方向位移為零。背襯處邊界為固定約束以模擬實(shí)際水下航行器應(yīng)用條件。表1給出了本文所分析聲學(xué)結(jié)構(gòu)的靜態(tài)材料參數(shù)。圖5給出了不同靜水壓作用下無金屬孔板的傳統(tǒng)聲學(xué)覆蓋層的形變情況。從圖5中可以看出,帶有球腔的傳統(tǒng)聲學(xué)覆蓋層在靜水壓力載荷的作用下,會(huì)產(chǎn)生較大形變。隨著靜水壓的增加,結(jié)構(gòu)的變形越大,空腔占有率(穿孔率)越低。圖6給出了不同靜水壓作用下帶有金屬孔板的聲學(xué)覆蓋層的形變情況。從圖6中可以看出,插入金屬孔板后,聲學(xué)覆蓋層的受壓形變量有所減小,并隨著壓力的增大,形變量減小得更明顯。圖7所示為傳統(tǒng)聲學(xué)覆蓋層和帶有金屬孔板聲學(xué)覆蓋層在6 MPa靜水壓力下的內(nèi)部球腔形變對(duì)比??梢钥闯?,相對(duì)于傳統(tǒng)聲學(xué)覆蓋層,在6 MPa的高靜水壓力下,插入金屬孔板的聲學(xué)覆蓋層中的空腔形變量減少了5%以上。因此可以得出,因金屬的楊氏模量較大,通過在傳統(tǒng)聲學(xué)覆蓋層中插入金屬孔板,相當(dāng)于提高了聲學(xué)覆蓋層整體的等效模量,在一定程度上提高聲學(xué)結(jié)構(gòu)在大潛深、高靜水壓力下的耐壓性能。

    2.2 不同靜壓下聲學(xué)覆蓋層的吸聲性能分析

    為了研究靜水壓力對(duì)聲學(xué)覆蓋層吸聲性能的影響規(guī)律,本文分別計(jì)算了在不同靜水壓力下無金屬孔板與帶金屬孔板兩種聲學(xué)結(jié)構(gòu)的吸聲系數(shù)。在進(jìn)行聲學(xué)計(jì)算時(shí),黏彈性材料參數(shù)選用由實(shí)驗(yàn)擬合數(shù)據(jù)得到的動(dòng)態(tài)力學(xué)參數(shù)[27],此時(shí)作為基體材料橡膠的楊氏模量E與損耗因子η隨頻率f變化為E=8.9×106+2.3×103-10?2和η=0.2+1.2×10?4f-8.2×10?9。此外,計(jì)算過程中選用海水作為流體域,其密度和聲速分別為1070 kg/m3和1540 m/s。

    圖8為不受靜水壓力下帶有金屬孔板的聲學(xué)覆蓋層與無金屬孔板的傳統(tǒng)聲學(xué)覆蓋層的吸聲性能對(duì)比。如圖8所示,與傳統(tǒng)聲學(xué)覆蓋層相比,本文提出的帶有金屬孔板的聲學(xué)覆蓋層的吸聲系數(shù)在中、高頻段增大。

    圖9給出了常壓下無金屬孔板聲學(xué)覆蓋層與帶金屬孔板聲學(xué)覆蓋層在低頻吸聲峰(601 Hz)和中、高頻處(1851 Hz, 6000 Hz)的振動(dòng)能量模態(tài)圖。從圖9可以得出,在低頻吸聲峰頻率處,由于空腔的諧振作用,空腔上方的基體層沿垂直于覆蓋層方向產(chǎn)生結(jié)構(gòu)彎曲振動(dòng);此外,基體層作為彈性層和剛性背襯組成另一個(gè)共振系統(tǒng),其共同作用使結(jié)構(gòu)在較低頻率處產(chǎn)生了較強(qiáng)的吸聲效果,此時(shí)振動(dòng)能量主要集中在背襯與基體。當(dāng)在傳統(tǒng)聲學(xué)覆蓋層中插入金屬孔板后,并沒有改變背襯與基體的組成結(jié)構(gòu),兩種聲學(xué)覆蓋層的空腔結(jié)構(gòu)的穿孔率一致,所以二者在低頻處的吸聲效果基本一致。而在中、高頻段,由于基體內(nèi)傳遞的聲波在空腔邊界處產(chǎn)生散射,散射波之間相長(zhǎng)干涉,導(dǎo)致聲波的波形轉(zhuǎn)換,縱波轉(zhuǎn)換為損耗更大的剪切波,增加了聲能耗散。而相比于傳統(tǒng)聲學(xué)覆蓋層,當(dāng)在聲學(xué)覆蓋層插入金屬孔板后,金屬孔板本身及孔板的空腔結(jié)構(gòu),都增強(qiáng)了對(duì)聲波的散射作用,加大了波形轉(zhuǎn)換強(qiáng)度,進(jìn)而加強(qiáng)了對(duì)能量的耗散。如圖8所示,在中、高頻段,當(dāng)插入金屬孔板后加強(qiáng)了聲學(xué)覆蓋層的吸聲效果。

    圖10和11分別為無金屬孔板聲學(xué)覆蓋層與帶金屬孔板的聲學(xué)覆蓋層在不同靜水壓力下的吸聲系數(shù)對(duì)比圖。從圖中可以得出,當(dāng)靜水壓力增大時(shí),兩種聲學(xué)覆蓋層吸聲曲線的變化趨勢(shì)與常壓時(shí)保持一致。但隨著靜水壓力的增大,兩種聲學(xué)結(jié)構(gòu)的空腔變小,使得聲學(xué)覆蓋層穿孔率變小,故兩者的吸聲曲線總體上均向高頻方向變化,且在中、高頻段,兩種結(jié)構(gòu)的吸聲效果有所增強(qiáng)。另外,兩種聲學(xué)覆蓋層的低頻吸聲峰向高頻移動(dòng)的同時(shí),吸聲系數(shù)有所增大。

    圖12對(duì)比了在不同靜水壓力下無金屬孔板的傳統(tǒng)聲學(xué)覆蓋層與帶金屬孔板的聲學(xué)覆蓋層吸聲系數(shù)隨頻率的變化規(guī)律。從圖10~12中可以看出,兩種聲學(xué)覆蓋層在靜水壓力下的吸聲系數(shù)的變化規(guī)律與常壓下的變化趨勢(shì)一致,在中、高頻段帶有金屬孔板的聲學(xué)覆蓋層的吸聲系數(shù)大于無金屬孔板的聲學(xué)覆蓋層。而靜水壓力相同的條件下,兩種聲學(xué)覆蓋層在低頻吸聲峰的吸聲效果基本一致。

    為了更好地揭示在靜水壓力下,帶金屬孔板和無金屬孔板聲學(xué)覆蓋層在變形后的吸聲機(jī)理,本文分析了6 MPa靜水壓力下兩種聲學(xué)覆蓋層的振動(dòng)能量聚集形式。圖13給出了6 MPa靜水壓力下兩種聲學(xué)覆蓋層在低頻吸聲峰(1200 Hz)和中、高頻處(1851 Hz, 6000 Hz)的振動(dòng)能量模態(tài)圖。從圖13中可以看出,由于在靜水壓力下,兩種聲學(xué)覆蓋層的穿孔率變小,所以二者的低頻吸聲峰向高頻移動(dòng)并略有升高,此時(shí)的振動(dòng)能量主要是由空腔諧振以及背襯與基體產(chǎn)生的結(jié)構(gòu)共振引起的。在聲學(xué)覆蓋層中插入金屬孔板后,對(duì)此影響不大,故兩者在低頻的吸聲效果基本一致。在中、高頻段,空腔結(jié)構(gòu)受壓變形后,增強(qiáng)了對(duì)入射聲波的散射,改變了聲波的傳播方向,加大了聲波的傳播距離,并伴隨著縱?橫波的轉(zhuǎn)換,加劇了結(jié)構(gòu)的彎曲振動(dòng),使空腔諧振作用加強(qiáng),增加了聲波的能量損失。因此,如圖10~12所示,相比常壓條件,當(dāng)聲學(xué)覆蓋層受到靜水壓力作用后,兩種聲學(xué)覆蓋層在中、高頻段的吸聲效果得到提高,且隨著靜水壓力的提高,中、高頻段的吸聲效果提高得更加明顯。此外,通過對(duì)結(jié)構(gòu)的吸聲機(jī)理分析可知,金屬孔板的小空腔結(jié)構(gòu)是所設(shè)計(jì)結(jié)構(gòu)提高聲學(xué)性能的關(guān)鍵。金屬材質(zhì)作為空腔的支撐結(jié)構(gòu),能夠使其在靜水壓力條件下保持空腔結(jié)構(gòu)形式,以提高靜水壓力條件下結(jié)構(gòu)聲學(xué)性能的穩(wěn)定性。

    3 金屬孔板影響覆蓋層吸聲性能的規(guī)律分析

    為了進(jìn)一步改善聲學(xué)覆蓋層的吸聲性能,本文研究了金屬孔板材料、厚度、孔隙率對(duì)所設(shè)計(jì)聲學(xué)覆蓋層的聲學(xué)性能的影響規(guī)律。

    3.1 金屬孔板材料對(duì)覆蓋層吸聲性能的影響

    在保證其他材料參數(shù)不變的情況下,本文研究了不同金屬孔板材料對(duì)所設(shè)計(jì)聲學(xué)結(jié)構(gòu)的吸聲性能的影響規(guī)律。金屬孔板材料參數(shù)如表2所示。為探究靜水壓力對(duì)覆蓋層內(nèi)金屬孔板的影響,本文給出了在6 MPa靜水壓力作用下,聲學(xué)覆蓋層內(nèi)3種不同材質(zhì)的金屬孔板的形變情況,如圖14所示。從圖中可以看出,在6 MPa靜水壓力下,覆蓋層內(nèi)以鉛、鋁、鋼作為金屬孔板材料時(shí),金屬孔板在靜水壓力作用方向分別產(chǎn)生了4.6°,3.1°和1°左右的弧度變形。因金屬材料的楊氏模量較大,具有一定的耐壓、抗變形能力,在靜水壓力的作用下,覆蓋層內(nèi)的金屬孔板的形狀參數(shù)變化較小,且6 MPa是本文研究的最大靜水壓力,在1~4.5 MPa的外壓下,覆蓋層內(nèi)金屬孔板形狀參數(shù)的變化幾乎可以忽略。因此,本節(jié)以及后續(xù)章節(jié)金屬孔板材料參數(shù)對(duì)覆蓋層吸聲性能的影響規(guī)律的研究主要是在常壓下進(jìn)行的。

    圖15給出了常壓下不同金屬孔板材料的聲學(xué)覆蓋層吸聲系數(shù)對(duì)比。結(jié)果表明,當(dāng)金屬孔板材料為鉛時(shí),聲學(xué)覆蓋層的吸聲性能在中、高頻段優(yōu)于其他兩種材料。根據(jù)圖9所示的振動(dòng)能量分布,低頻吸聲峰主要取決于空腔諧振以及背襯與基體產(chǎn)生的結(jié)構(gòu)共振,所以金屬孔板材料的改變對(duì)低頻吸聲峰幾乎沒有影響。而鉛的質(zhì)量密度較大,使基體層的等效阻抗增大,阻抗越大,在中、高頻段對(duì)聲波能量的衰減耗散更為充分。

    3.2 金屬孔板厚度對(duì)覆蓋層吸聲性能的影響

    保證其他材料參數(shù)不變,研究金屬板厚度變化對(duì)所設(shè)計(jì)聲學(xué)結(jié)構(gòu)吸聲性能的影響規(guī)律。圖16為常壓下,1 mm,2 mm,3 mm厚金屬孔板的聲學(xué)覆蓋層吸聲系數(shù)對(duì)比。從圖16中可以看出,聲學(xué)覆蓋層中的金屬孔板厚度越薄,在中、高頻的吸聲效果越好。這是因?yàn)?,金屬孔板較薄時(shí),在中、高頻段,受聲波激勵(lì)產(chǎn)生的振幅更大,振動(dòng)更為劇烈,而產(chǎn)生的彎曲振動(dòng)能量越大,聲能損失越多。所以,薄金屬孔板的聲學(xué)覆蓋層在中、高頻段獲得的吸聲效果較好。

    3.3 金屬孔板孔隙率對(duì)覆蓋層吸聲性能的影響

    保證其他研究參數(shù)不變,研究金屬孔板孔隙率變化對(duì)聲學(xué)覆蓋層吸聲性能的影響規(guī)律。圖17所示為常壓下,金屬孔板孔隙率分別為16.3%和31.4%以及無孔金屬板的聲學(xué)覆蓋層的吸聲性能對(duì)比圖。從圖17中可以看出,隨著孔隙率的增大,能夠有效改善該聲學(xué)結(jié)構(gòu)在中、高頻段的吸聲效果,但較高的孔隙率也加大了工藝加工的難度,阻礙了實(shí)際應(yīng)用。

    4 結(jié) 論

    本文設(shè)計(jì)了一種含有金屬孔板的聲學(xué)覆蓋層結(jié)構(gòu),以靜力學(xué)方法研究了該聲學(xué)結(jié)構(gòu)在不同靜水壓力下的形變情況。通過建立靜水壓力下的聲學(xué)有限元方程,研究了該結(jié)構(gòu)在靜水壓0~6 MPa范圍內(nèi)的吸聲特性。在傳統(tǒng)空腔型聲學(xué)覆蓋層內(nèi)插入金屬孔板后,在保持結(jié)構(gòu)低頻段的吸聲效果的基礎(chǔ)上,進(jìn)一步豐富了聲學(xué)結(jié)構(gòu)內(nèi)的吸聲機(jī)制,改善了該結(jié)構(gòu)在2~10 kHz頻帶范圍的吸聲效果,且在受到靜水壓力產(chǎn)生形變下,該聲學(xué)結(jié)構(gòu)在中、高頻段的吸聲效果依然優(yōu)于傳統(tǒng)聲學(xué)覆蓋層,提高了靜水壓力下結(jié)構(gòu)聲學(xué)性能的穩(wěn)定性。此外,研究表明,調(diào)整孔板厚度、孔板材料及孔板孔隙率可以進(jìn)一步提高該聲學(xué)覆蓋層的吸聲性能。因此,帶金屬孔板聲學(xué)覆蓋層的結(jié)構(gòu)設(shè)計(jì)能夠在靜水壓力條件下改善結(jié)構(gòu)的吸聲性能,兼顧了結(jié)構(gòu)耐壓性與吸聲性能。

    參考文獻(xiàn):

    [1]"""" Zhang T X, Liu S W, He X, et al. Underwater target tracking using forward-looking sonar for autonomous underwater vehicles[J]. Sensors, 2019, 20(1): 1-28.

    [2]"""" 劉國強(qiáng), 樓京俊, 何其偉. 多層材料組合腔型吸聲覆蓋層的吸聲特性研究[J]. 武漢理工大學(xué)學(xué)報(bào):交通科學(xué)與工程版, 2016, 40(5): 850-858.

    Liu Guoqiang, Lou Jingjun, He Qiwei. Absorption characteristics of multi-layered material anechoic coating on composite cavities[J]. Journal of Wuhan University of Technology (Transportation Science amp; Engineering), 2016, 40(5): 850-858.

    [3]"""" 高南沙, 侯宏. 三維局域共振型聲子晶體低頻帶隙特性研究[J]. 材料導(dǎo)報(bào), 2018, 32(1): 322-326.

    Gao Nansha, Hou Hong. Low frequency bandgap characteristics of three-dimensional local resonance phononic crystal[J]. Material Reports, 2018, 32(1): 322-326.

    [4]"""" Zhao H G, Wang Y, Yu D L, et al. A double porosity material for low frequency sound absorption[J]. Composite Structures, 2020, 239: 111978.

    [5]"""" 師康康, 靳國永, 葉天貴, 等. 含空腔的功能梯度聲學(xué)覆蓋層水下吸聲特性[J]. 聲學(xué)學(xué)報(bào), 2021, 46(3): 394-404.

    Shi Kangkang, Jin Guoyong, Ye Tiangui,et al. Underwater sound absorption performance of functionally graded anechoic coating with cavities[J]. Acta Acustica, 2021, 46(3): 394-404.

    [6]"""" 侯九霄, 朱海潮, 廖金龍,等. 水下彈性微穿孔吸聲結(jié)構(gòu)吸聲系數(shù)研究[J]. 聲學(xué)學(xué)報(bào), 2021, 46(1): 135-142.

    Hou Jiuxiao, Zhu Haichao, Liao Jinlong,et al. Sound absorption characteristics of underwater flexible micro-perforated sound absorption structure[J]. Acta Acustica, 2021, 46(1): 135-142.

    [7]"""" 王育人, 繆旭弘, 姜恒,等. 水下吸聲機(jī)理與吸聲材料[J]. 力學(xué)進(jìn)展, 2017, 47(1): 92-121.

    Wang Yuren, Miao Xuhong, Jiang Heng, et al. Review on underwater sound absorption materials and mechanisms[J]. Advances in Mechanics, 2017, 47(1): 92-121.

    [8]"""" 馬驍晨, 李建龍, 宋昊, 等. 聲學(xué)覆蓋層吸聲系數(shù)的多通道逆濾波測(cè)量[J]. 聲學(xué)學(xué)報(bào), 2019, 44(4): 726-734.

    Ma Xiaochen, Li Jianlong, Song Hao, et al. Measurement of the absorption coefficient for acoustical coating samples using multichannel inverse filter[J]. Acta Acustica, 2019, 44(4): 726-734.

    [9]"""" 周帥龍, 陳理添, 劉小俠. 聲學(xué)覆蓋層的優(yōu)化設(shè)計(jì)[J]. 噪聲與振動(dòng)控制, 2021, 41(4): 35-41.

    Zhou Shuailong, Chen Litian, Liu Xiaoxia. Optimal design of anechoic coatings[J]. Noise and Vibration Control, 2021, 41(4): 35-41.

    [10]""" 甄冬, 王梓宇, 焦湘和, 等. 微穿孔板在隨機(jī)-區(qū)間混合不確定性理論下的結(jié)構(gòu)優(yōu)化[J]. 振動(dòng)工程學(xué)報(bào), 2022, 35(1): 228-235.

    Zhen Dong, Wang Ziyu, Jiao Xianghe, et al. Structural optimization of micro-perforated panel using random-interval mixed uncertainty theory[J]. Journal of Vibration Engineering, 2022, 35(1): 228-235.

    [11]""" 寧景鋒, 趙桂平, 穆林, 等. 含有空氣背襯層的分層多孔材料的吸聲性能研究[J]. 振動(dòng)工程學(xué)報(bào), 2014, 27(5): 734-740.

    Ning Jingfeng, Zhao Guiping, Mu Lin, et al. Study on sound absorption properties of multilayer porous material structure backed with an air gap[J]. Journal of Vibration Engineering, 2014, 27(5): 734-740.

    [12]""" Zhang Y N, Cheng L. Ultra-thin and broadband low-frequency underwater acoustic meta-absorber[J]. International Journal of Mechanical Sciences, 2021, 210: 106732.

    [13]""" Wang L B, Ma C Z, Wu J H, et al. Realizing high-efficiency low frequency sound absorption of underwater meta-structures by acoustic siphon effect[J]. Modern Physics Letters B, 2021, 35(19): 2150319.

    [14]""" 姜聞文, 陳光冶, 朱彥. 靜水壓變化下橡膠結(jié)構(gòu)吸聲性能的計(jì)算與分析[J]. 噪聲與振動(dòng)控制, 2006, 26(5): 55-57.

    Jiang Wenwen, Chen Guangye, Zhu Yan. Computation and analysis of sound absorption performance of rubber structures under variable hydraulic pressure[J]. Noise and Vibration Control, 2006, 26(5): 55-57.

    [15]""" 姚熊亮, 劉文賀, 劉慶杰, 等. 水深與腔型對(duì)隔聲去耦瓦吸聲系數(shù)的影響[J]. 哈爾濱工程大學(xué)學(xué)報(bào), 2007, 28(6): 605-610.

    Yao Xiongliang, Liu Wenhe, Liu Qingjie,et al. Influence of depth and cavity shape on absorption coefficients of sound-isolating decoupled tiles[J]. Journal of Harbin Engineering University, 2007, 28(6): 605-610.

    [16]""" Humphrey V F, Robinson S P, Smith J D, et al. Acoustic characterization of panel materials under simulated ocean conditions using a parametric array source[J]. Journal of the Acoustical Society of America, 2008, 124(2): 803-814.

    [17]""" 張沖, 何世平, 易少強(qiáng). 靜壓下球形空腔吸聲覆蓋層的建模與性能分析[J]. 船舶力學(xué), 2017, 21(1): 99-106.

    Zhang Chong, He Shiping, Yi Shaoqiang. Model and absorption performance of anechoic coating embedding sphere cavities[J]. Journal of Ship Mechanics, 2017, 21(1): 99-106.

    [18]""" Thieury M,Tourin A,Dassé J, et al. Effect of hydrostatic pressure on a bubble anechoic metascreen[C]//13th International Congress on Artificial Materials for Novel Wave Phenomena-Metamaterials, Italy, 2019.

    [19]""" Gu Y H, Zhou H B, Bao B, et al. Experimental investigation of underwater locally multi-resonant metamaterials under high hydrostatic pressure for low frequency sound absorption[J]. Applied Acoustics, 2021, 172: 107605.

    [20]""" 董文凱, 陳美霞. 靜壓下考慮腔壓的吸聲覆蓋層吸聲性能分析[J].中國艦船研究,2022, 17(1):132-140.

    Dong Wenkai, Chen Meixia. Sound absorption performance analysis of anechoic coating under hydrostatic pressure considering cavity pressure[J]. Chinese Journal of Ship Research, 2022, 17(1):132-140.

    [21]""" Zhang Z F, Li S D, Wang J X, et al. Low-frequency broadband absorption of semi-active composite anechoic coating with subwavelength piezoelectric arrays in hydrostatic environments[J]. Results in Physics, 2021, 30: 104879.

    [22]""" 黃修長(zhǎng), 朱蓓麗, 胡碰, 等. 靜水壓力下橡膠動(dòng)態(tài)力學(xué)參數(shù)的聲管測(cè)量方法[J].上海交通大學(xué)學(xué)報(bào), 2013, 47(10): 1503-1519.

    Huang Xiuchang, Zhu Beili, Hu Peng, et al. Measurement of dynamic properties of rubber under hydrostatic pressure by water-filled acoustic tube[J]. Journal of Shanghai Jiao Tong University, 2013, 47(10): 1503-1519.

    [23]""" 陶猛. 靜水壓力下高分子材料黏彈性動(dòng)力學(xué)參數(shù)測(cè)量[J]. 振動(dòng)與沖擊, 2015, 35(6): 59-63.

    Tao Meng. Measurement of viscoelastic dynamic parameters of polymer materials under hydrostatic pressure[J]. Journal of Vibration and Shock, 2016, 35(6): 59-63.

    [24]""" 王曼. 水聲吸聲覆蓋層理論與實(shí)驗(yàn)研究[D]. 哈爾濱: 哈爾濱工程大學(xué), 2004.

    Wang Man. Theoretical and experimental study on underwater anechoic coating[D]. Harbin: Harbin Engineering University, 2004.

    [25]""" Yang H B, Zhao H G, Wen J H. Theory and numerical method for the effects of hydrostatic pressure on sound absorption of underwater acoustic coatings with air cavities[J]. Journal of Sound and Vibration, 2022, 533: 116985.

    [26]""" Sharma G S,Skvortsov A,MacGillivray I, et al. Sound absorption by rubber coatings with periodic voids and hard inclusions[J]. Applied Acoustics, 2019, 143: 200-210.

    [27]""" 陶猛, 卓琳凱. 靜水壓力下吸聲覆蓋層的聲學(xué)性能分析[J]. 上海交通大學(xué)學(xué)報(bào), 2011, 45(9): 1340-1350.

    Tao Meng, Zhuo Linkai. Effect of hydrostatic pressure on acoustic performance of sound absorption coating[J]. Journal of Shanghai Jiao Tong University, 2011, 45(9): 1340-1350.

    Characteristics of pressure resistance and sound absorption on anechoic coating with metal perforated plate under hydrostatic pressure

    JIA Xin-yu, JIN Guo-yong, YE Tian-gui, YAN Yan

    (College of Power and Energy Engineering, Harbin Engineering University, Harbin 150001, China)

    Abstract: An underwater acoustic structure with metal perforated plate inserted into the traditional anechoic coating is proposed to improve the pressure resistance and sound absorption performance of the structure. The deformation of the acoustic structure under different hydrostatic pressure is studied by the static method. By establishing the acoustic finite element equation under the action of hydrostatic pressure, the sound absorption effect of the acoustic structure under different hydrostatic pressure is analyzed. Compared with the traditional anechoic coating, under the hydrostatic pressure of 0 to 6 MPa, the acoustic structure achieves better broadband sound absorption in the mid-to-high frequency range. In addition, the effects of the thickness, material and porosity of the metal perforated plate on the sound absorption performance of the acoustic structure are discussed. The research shows that the underwater anechoic coating with metal perforated plate is an effective design to improve the pressure resistance and sound absorption performance of underwater acoustic structure.

    Key words: anechoic coating;metal perforated plate;hydrostatic pressure;deformation;sound absorption performance

    作者簡(jiǎn)介: 賈薪宇(1989―),男,博士研究生。E-mail:jiaxy@hrbeu.edu.cn。

    通訊作者: 靳國永(1980―),男,博士,教授。E-mail:guoyongjin@hrbeu.edu.cn。

    猜你喜歡
    靜水壓力
    植物根壓研究進(jìn)展
    廣西植物(2022年4期)2022-05-13 08:24:53
    二次供水豎向分區(qū)技術(shù)分析
    如何做好救生筏壓力釋放器及相關(guān)部件的連接
    中國船檢(2018年6期)2018-06-22 09:40:22
    不同靜水壓力對(duì)大鼠髁突軟骨細(xì)胞的TRPM7和Fas蛋白表達(dá)的影響
    老空水特殊存在形式的猜想
    關(guān)于小高層住宅建筑給水排水設(shè)計(jì)分析
    水下環(huán)肋圓柱殼的振動(dòng)耦合頻率影響因素分析
    基于數(shù)值-解析法測(cè)量靜壓條件下阻尼材料動(dòng)態(tài)力學(xué)參數(shù)
    靜水壓力下高分子材料黏彈性動(dòng)力學(xué)參數(shù)測(cè)量
    巖土類材料的靜水壓力效應(yīng)分析
    山西建筑(2015年23期)2015-04-05 21:06:46
    国产一区亚洲一区在线观看| 亚洲欧美日韩另类电影网站| 亚洲第一av免费看| 天堂中文最新版在线下载| 精品熟女少妇八av免费久了| 免费在线观看日本一区| 日韩一区二区三区影片| 欧美在线黄色| 国产av精品麻豆| 香蕉国产在线看| av一本久久久久| 2021少妇久久久久久久久久久| 亚洲 国产 在线| www.999成人在线观看| 18在线观看网站| 亚洲国产中文字幕在线视频| 一级a爱视频在线免费观看| tube8黄色片| 亚洲av电影在线进入| 欧美日韩一级在线毛片| 国产成人精品久久二区二区免费| 精品人妻一区二区三区麻豆| 免费女性裸体啪啪无遮挡网站| 少妇人妻 视频| 啦啦啦在线免费观看视频4| 精品国产超薄肉色丝袜足j| 少妇被粗大的猛进出69影院| 亚洲一区二区三区欧美精品| 国产日韩一区二区三区精品不卡| 婷婷成人精品国产| 九草在线视频观看| 精品一区二区三区四区五区乱码 | 飞空精品影院首页| 女人久久www免费人成看片| 亚洲天堂av无毛| 91老司机精品| 999精品在线视频| 深夜精品福利| 日韩伦理黄色片| 久久午夜综合久久蜜桃| 国产一级毛片在线| 成人亚洲精品一区在线观看| 成人国语在线视频| 2021少妇久久久久久久久久久| 午夜福利视频在线观看免费| 国产福利在线免费观看视频| 激情五月婷婷亚洲| 精品国产超薄肉色丝袜足j| 欧美黄色片欧美黄色片| 青草久久国产| 色视频在线一区二区三区| 亚洲三区欧美一区| 久久精品国产亚洲av涩爱| 久久九九热精品免费| 国产精品一区二区免费欧美 | 久久天堂一区二区三区四区| 亚洲色图 男人天堂 中文字幕| 黄片播放在线免费| 亚洲国产欧美日韩在线播放| 波野结衣二区三区在线| 大码成人一级视频| 麻豆乱淫一区二区| 亚洲,欧美,日韩| 亚洲熟女精品中文字幕| 精品卡一卡二卡四卡免费| 黑人巨大精品欧美一区二区蜜桃| 老司机影院毛片| 日本vs欧美在线观看视频| 久久 成人 亚洲| 999久久久国产精品视频| 日韩欧美一区视频在线观看| 免费高清在线观看日韩| 欧美黄色淫秽网站| 天堂8中文在线网| 秋霞在线观看毛片| 亚洲,欧美,日韩| 日本午夜av视频| 国产av国产精品国产| 久久久久久久大尺度免费视频| 十八禁网站网址无遮挡| 亚洲欧美成人综合另类久久久| 亚洲成人手机| 波多野结衣av一区二区av| 免费观看人在逋| 精品久久蜜臀av无| 老司机在亚洲福利影院| 欧美人与善性xxx| 观看av在线不卡| 亚洲,一卡二卡三卡| 男人舔女人的私密视频| 男女免费视频国产| 国产高清videossex| 亚洲熟女毛片儿| 亚洲久久久国产精品| 日本av免费视频播放| 啦啦啦在线免费观看视频4| 黑人巨大精品欧美一区二区蜜桃| 日本a在线网址| 大片电影免费在线观看免费| 性色av一级| 99re6热这里在线精品视频| 婷婷色综合大香蕉| 一本大道久久a久久精品| 日韩视频在线欧美| 久久九九热精品免费| 久久精品成人免费网站| 久久精品久久精品一区二区三区| 国产亚洲精品第一综合不卡| 亚洲国产欧美一区二区综合| 亚洲av电影在线观看一区二区三区| 性色av一级| 欧美激情极品国产一区二区三区| 久久久精品国产亚洲av高清涩受| 中文乱码字字幕精品一区二区三区| 考比视频在线观看| 女警被强在线播放| 欧美国产精品va在线观看不卡| 精品国产一区二区三区四区第35| 欧美日韩视频精品一区| 国精品久久久久久国模美| 精品人妻1区二区| 免费女性裸体啪啪无遮挡网站| 日韩大片免费观看网站| 午夜两性在线视频| 亚洲九九香蕉| 一本—道久久a久久精品蜜桃钙片| www日本在线高清视频| 久久精品aⅴ一区二区三区四区| cao死你这个sao货| 夫妻午夜视频| 国产欧美日韩一区二区三区在线| 亚洲欧美精品自产自拍| 国产亚洲午夜精品一区二区久久| 啦啦啦啦在线视频资源| 免费观看av网站的网址| 青春草亚洲视频在线观看| 国产精品一二三区在线看| 大香蕉久久网| 高清av免费在线| 国产日韩欧美视频二区| 大香蕉久久网| √禁漫天堂资源中文www| 国产日韩欧美视频二区| 香蕉丝袜av| 免费av中文字幕在线| 久久精品久久久久久久性| 啦啦啦 在线观看视频| 操美女的视频在线观看| 97人妻天天添夜夜摸| 黄片小视频在线播放| 操出白浆在线播放| 亚洲欧美精品综合一区二区三区| 视频区欧美日本亚洲| 日韩制服骚丝袜av| 欧美精品亚洲一区二区| 人人妻人人澡人人爽人人夜夜| 可以免费在线观看a视频的电影网站| 黄色视频在线播放观看不卡| 国产三级黄色录像| 两人在一起打扑克的视频| 777久久人妻少妇嫩草av网站| 亚洲av美国av| 日本vs欧美在线观看视频| 9热在线视频观看99| 亚洲,欧美,日韩| 日本vs欧美在线观看视频| 精品福利观看| 久久精品人人爽人人爽视色| 一级片'在线观看视频| 亚洲中文字幕日韩| 国产精品一区二区在线观看99| 久久久久国产一级毛片高清牌| 久久人人爽人人片av| 女人久久www免费人成看片| 91麻豆av在线| 伊人久久大香线蕉亚洲五| 国产成人精品久久久久久| 久久ye,这里只有精品| 男人舔女人的私密视频| 视频在线观看一区二区三区| 丰满少妇做爰视频| 成人午夜精彩视频在线观看| 国产成人一区二区三区免费视频网站 | 国产一级毛片在线| 老鸭窝网址在线观看| 亚洲,一卡二卡三卡| 国产在线观看jvid| 国产极品粉嫩免费观看在线| 久久人人爽人人片av| 亚洲成色77777| 亚洲精品国产一区二区精华液| 国产视频一区二区在线看| 女人久久www免费人成看片| 国产精品免费大片| 亚洲国产精品国产精品| 久久精品亚洲av国产电影网| 中文乱码字字幕精品一区二区三区| 久久精品国产亚洲av涩爱| 国产一区亚洲一区在线观看| 好男人电影高清在线观看| 精品福利永久在线观看| 免费在线观看日本一区| 久久毛片免费看一区二区三区| 精品国产乱码久久久久久小说| 女性被躁到高潮视频| 91精品三级在线观看| 久久久久国产一级毛片高清牌| 欧美日韩黄片免| 精品免费久久久久久久清纯 | 亚洲九九香蕉| 国产99久久九九免费精品| 日韩视频在线欧美| 亚洲人成电影免费在线| 两个人看的免费小视频| 免费在线观看影片大全网站 | 亚洲免费av在线视频| 精品福利观看| 国产伦人伦偷精品视频| 久久99热这里只频精品6学生| 国产成人系列免费观看| 制服诱惑二区| 97人妻天天添夜夜摸| 视频区欧美日本亚洲| www.熟女人妻精品国产| 男人舔女人的私密视频| 免费人妻精品一区二区三区视频| 下体分泌物呈黄色| 一级毛片黄色毛片免费观看视频| 日日摸夜夜添夜夜爱| 午夜福利乱码中文字幕| 爱豆传媒免费全集在线观看| 嫩草影视91久久| 亚洲激情五月婷婷啪啪| 久久精品亚洲熟妇少妇任你| 黑丝袜美女国产一区| 欧美精品av麻豆av| 久久精品熟女亚洲av麻豆精品| 欧美日韩一级在线毛片| 一二三四在线观看免费中文在| 波野结衣二区三区在线| 国产av国产精品国产| 一区二区日韩欧美中文字幕| 精品久久蜜臀av无| 亚洲精品国产av蜜桃| 久久精品成人免费网站| 国产片特级美女逼逼视频| 亚洲中文av在线| 欧美少妇被猛烈插入视频| 亚洲国产精品国产精品| 日本av免费视频播放| a级毛片黄视频| 亚洲成色77777| 两个人看的免费小视频| 美女中出高潮动态图| 国产成人一区二区三区免费视频网站 | 欧美黑人精品巨大| 777米奇影视久久| 999久久久国产精品视频| 久久久国产精品麻豆| 久久午夜综合久久蜜桃| 久久精品aⅴ一区二区三区四区| 国产成人91sexporn| 久久精品久久久久久久性| 大型av网站在线播放| 丝袜美足系列| 大香蕉久久网| 啦啦啦中文免费视频观看日本| 国产精品香港三级国产av潘金莲 | 欧美精品高潮呻吟av久久| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜两性在线视频| 亚洲av国产av综合av卡| 精品福利观看| 国产女主播在线喷水免费视频网站| 久久青草综合色| 国产欧美日韩一区二区三 | 久久人人爽av亚洲精品天堂| 国产精品久久久人人做人人爽| 日本一区二区免费在线视频| 亚洲精品国产区一区二| 免费高清在线观看视频在线观看| 一级a爱视频在线免费观看| 亚洲成色77777| 人人澡人人妻人| 精品一品国产午夜福利视频| 女人精品久久久久毛片| 亚洲精品国产av蜜桃| 久久性视频一级片| 天堂俺去俺来也www色官网| 黄频高清免费视频| 亚洲自偷自拍图片 自拍| 久久精品人人爽人人爽视色| 国产免费视频播放在线视频| 老汉色av国产亚洲站长工具| 男女高潮啪啪啪动态图| 亚洲精品久久成人aⅴ小说| 一本色道久久久久久精品综合| 国产欧美日韩一区二区三区在线| 99精品久久久久人妻精品| 亚洲中文字幕日韩| 欧美日韩成人在线一区二区| 一级a爱视频在线免费观看| 久久精品国产亚洲av高清一级| 9191精品国产免费久久| 国产成人精品在线电影| 大片免费播放器 马上看| 捣出白浆h1v1| 我的亚洲天堂| 亚洲成av片中文字幕在线观看| 国产视频首页在线观看| 午夜免费男女啪啪视频观看| 欧美精品av麻豆av| 国产一卡二卡三卡精品| 久久国产精品人妻蜜桃| 国产不卡av网站在线观看| 一边亲一边摸免费视频| 欧美在线黄色| 一边亲一边摸免费视频| 国产精品麻豆人妻色哟哟久久| 菩萨蛮人人尽说江南好唐韦庄| 欧美亚洲日本最大视频资源| 午夜av观看不卡| 亚洲五月色婷婷综合| 在线看a的网站| 女人爽到高潮嗷嗷叫在线视频| 欧美激情极品国产一区二区三区| 日日爽夜夜爽网站| 亚洲av日韩精品久久久久久密 | 少妇被粗大的猛进出69影院| 日韩av免费高清视频| 无遮挡黄片免费观看| 国产麻豆69| 成年av动漫网址| 最黄视频免费看| 国产精品久久久久久精品古装| 最新的欧美精品一区二区| 曰老女人黄片| 国产成人精品久久二区二区91| 婷婷色综合www| 精品国产一区二区三区久久久樱花| 丝袜美足系列| 亚洲九九香蕉| 欧美国产精品va在线观看不卡| 久久人妻福利社区极品人妻图片 | 亚洲一区二区三区欧美精品| 秋霞在线观看毛片| 青草久久国产| 每晚都被弄得嗷嗷叫到高潮| av福利片在线| 久久人人97超碰香蕉20202| 日韩伦理黄色片| 国产高清视频在线播放一区 | 亚洲 欧美一区二区三区| 人人妻人人爽人人添夜夜欢视频| 高潮久久久久久久久久久不卡| 国产成人一区二区在线| 高潮久久久久久久久久久不卡| 人妻一区二区av| 国产日韩一区二区三区精品不卡| 国产精品二区激情视频| 亚洲 欧美一区二区三区| 日本黄色日本黄色录像| 日日爽夜夜爽网站| 夫妻性生交免费视频一级片| 一本综合久久免费| 最近最新中文字幕大全免费视频 | 日本av免费视频播放| 另类精品久久| 在线观看免费午夜福利视频| 亚洲精品一区蜜桃| 亚洲成国产人片在线观看| 国产亚洲精品第一综合不卡| 亚洲国产欧美在线一区| 在线观看免费视频网站a站| 午夜久久久在线观看| 久久亚洲精品不卡| 国产高清国产精品国产三级| 99国产精品免费福利视频| 我的亚洲天堂| 女人久久www免费人成看片| 少妇猛男粗大的猛烈进出视频| 19禁男女啪啪无遮挡网站| 脱女人内裤的视频| 看十八女毛片水多多多| 亚洲国产最新在线播放| 大话2 男鬼变身卡| 少妇 在线观看| 国产黄色免费在线视频| 热re99久久精品国产66热6| 中文欧美无线码| 国产亚洲一区二区精品| 美女主播在线视频| 亚洲少妇的诱惑av| 大片电影免费在线观看免费| 国产欧美日韩一区二区三 | 99九九在线精品视频| 人体艺术视频欧美日本| 日本欧美国产在线视频| 一级,二级,三级黄色视频| 亚洲欧美日韩另类电影网站| 欧美日韩视频精品一区| 亚洲图色成人| 大码成人一级视频| 国产又色又爽无遮挡免| 校园人妻丝袜中文字幕| 成人免费观看视频高清| 一本综合久久免费| 一级片免费观看大全| 久热爱精品视频在线9| 超色免费av| 亚洲av综合色区一区| av福利片在线| 丁香六月天网| √禁漫天堂资源中文www| 蜜桃在线观看..| 人人澡人人妻人| 欧美激情极品国产一区二区三区| 丰满少妇做爰视频| 国产成人精品久久二区二区91| 黄色视频在线播放观看不卡| 久久九九热精品免费| 纯流量卡能插随身wifi吗| 国产成人av教育| 精品亚洲成国产av| svipshipincom国产片| 亚洲第一青青草原| 91精品国产国语对白视频| 人体艺术视频欧美日本| 欧美日韩精品网址| 国产精品一国产av| 亚洲av成人不卡在线观看播放网 | 中国国产av一级| 国产精品二区激情视频| 精品卡一卡二卡四卡免费| 国产伦人伦偷精品视频| 侵犯人妻中文字幕一二三四区| 两性夫妻黄色片| 久久久久久久大尺度免费视频| 亚洲伊人久久精品综合| avwww免费| 国产一区二区 视频在线| 欧美激情极品国产一区二区三区| 最近手机中文字幕大全| 日韩一本色道免费dvd| 天堂俺去俺来也www色官网| 99香蕉大伊视频| 99热全是精品| 黄色片一级片一级黄色片| 国产日韩欧美亚洲二区| 欧美中文综合在线视频| 在线精品无人区一区二区三| 咕卡用的链子| 国产日韩欧美在线精品| 国产高清videossex| 老汉色av国产亚洲站长工具| 人妻人人澡人人爽人人| 亚洲激情五月婷婷啪啪| 纯流量卡能插随身wifi吗| 黄色片一级片一级黄色片| 国产深夜福利视频在线观看| 中国美女看黄片| 亚洲午夜精品一区,二区,三区| 亚洲人成电影观看| 91字幕亚洲| 久久久久国产一级毛片高清牌| 久久久欧美国产精品| 一区二区三区四区激情视频| 精品国产乱码久久久久久男人| 国产精品一区二区精品视频观看| 美女主播在线视频| 国产日韩欧美在线精品| 9色porny在线观看| 成在线人永久免费视频| 精品福利观看| 欧美日韩福利视频一区二区| 欧美av亚洲av综合av国产av| 丰满迷人的少妇在线观看| 国产精品.久久久| 在现免费观看毛片| 成年人黄色毛片网站| av天堂在线播放| 国产精品三级大全| 男女床上黄色一级片免费看| 国产在视频线精品| 亚洲人成电影观看| 乱人伦中国视频| 亚洲精品国产一区二区精华液| 人人妻人人爽人人添夜夜欢视频| 一区二区三区乱码不卡18| 欧美在线一区亚洲| 少妇粗大呻吟视频| 中文字幕色久视频| a级片在线免费高清观看视频| 久久99热这里只频精品6学生| 丝袜美腿诱惑在线| 又紧又爽又黄一区二区| 女警被强在线播放| 国产主播在线观看一区二区 | 亚洲av日韩精品久久久久久密 | 欧美精品高潮呻吟av久久| 一区二区三区精品91| 久久精品国产亚洲av高清一级| 91精品伊人久久大香线蕉| 成年人免费黄色播放视频| 精品少妇久久久久久888优播| 久久av网站| 黄片播放在线免费| 国产一区二区三区综合在线观看| 精品一区二区三区av网在线观看 | 国产极品粉嫩免费观看在线| 好男人电影高清在线观看| 黄色片一级片一级黄色片| 亚洲专区中文字幕在线| 久久久精品国产亚洲av高清涩受| 嫩草影视91久久| 少妇人妻久久综合中文| 韩国精品一区二区三区| 日本vs欧美在线观看视频| 涩涩av久久男人的天堂| 一本—道久久a久久精品蜜桃钙片| 欧美av亚洲av综合av国产av| 久久久国产一区二区| 欧美黑人精品巨大| 欧美大码av| 91麻豆精品激情在线观看国产 | 少妇的丰满在线观看| 男人添女人高潮全过程视频| 亚洲欧洲日产国产| 熟女av电影| 亚洲色图综合在线观看| 久久久国产精品麻豆| 一级毛片我不卡| 肉色欧美久久久久久久蜜桃| www日本在线高清视频| 午夜福利免费观看在线| 曰老女人黄片| 亚洲欧洲国产日韩| 国产精品国产av在线观看| 精品卡一卡二卡四卡免费| 视频在线观看一区二区三区| 女人高潮潮喷娇喘18禁视频| 搡老乐熟女国产| 欧美日韩亚洲综合一区二区三区_| 国产黄色视频一区二区在线观看| 亚洲av成人精品一二三区| 精品人妻在线不人妻| 午夜日韩欧美国产| 超色免费av| 无遮挡黄片免费观看| 欧美日韩视频高清一区二区三区二| 黄频高清免费视频| 欧美成人精品欧美一级黄| 亚洲精品国产一区二区精华液| 亚洲中文av在线| 久热爱精品视频在线9| 少妇人妻久久综合中文| 欧美黄色片欧美黄色片| 国产精品熟女久久久久浪| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品熟女亚洲av麻豆精品| 国产1区2区3区精品| 久久久久久久久久久久大奶| 国产在线一区二区三区精| 中国国产av一级| 性色av乱码一区二区三区2| 老司机靠b影院| 一边亲一边摸免费视频| 亚洲成色77777| 午夜两性在线视频| 国产一区亚洲一区在线观看| www.av在线官网国产| 欧美精品亚洲一区二区| 99热网站在线观看| 99国产精品免费福利视频| 亚洲国产精品成人久久小说| 日本wwww免费看| 午夜福利,免费看| 在现免费观看毛片| 视频区图区小说| 国产91精品成人一区二区三区 | 99久久综合免费| 18在线观看网站| 搡老岳熟女国产| 真人做人爱边吃奶动态| 欧美亚洲 丝袜 人妻 在线| 国产在线观看jvid| 真人做人爱边吃奶动态| 十八禁人妻一区二区| 欧美精品高潮呻吟av久久| 丁香六月欧美| 亚洲人成电影观看| 精品一区二区三卡| 亚洲av片天天在线观看| 成人手机av| 精品国产一区二区三区四区第35| 热re99久久国产66热| 久久精品亚洲av国产电影网| 美女中出高潮动态图| 99久久99久久久精品蜜桃| 亚洲美女黄色视频免费看| 国产又爽黄色视频| 日日摸夜夜添夜夜爱| 成年美女黄网站色视频大全免费| 久久国产精品影院| 在线观看免费午夜福利视频| 涩涩av久久男人的天堂| 国产精品三级大全| 亚洲成人国产一区在线观看 | 热99久久久久精品小说推荐| 亚洲情色 制服丝袜| 国产一区二区三区综合在线观看| av国产久精品久网站免费入址| 人体艺术视频欧美日本| 国产精品一区二区精品视频观看| 国产色视频综合|