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

    一種新型的水下低頻共鳴器?

    2022-07-07 07:23:32邱燕萍周浩萌尋天雨郁高坤
    應(yīng)用聲學(xué) 2022年3期
    關(guān)鍵詞:聲阻抗共振頻率聲場

    邱燕萍 周浩萌 尋天雨 郁高坤

    (中國海洋大學(xué) 青島 266100)

    0 引言

    聲音對海洋動物的生存非常重要。近年來由人類活動產(chǎn)生的水下噪聲嚴(yán)重地破壞了海洋環(huán)境平衡,正在威脅著海洋哺乳動物及其他水生生物的生命健康,引起了廣泛的關(guān)注。由人為活動產(chǎn)生的噪聲包括船舶噪聲、聲吶工程噪聲以及水下作業(yè)平臺噪聲等,對海洋哺乳動物的危害主要有生理損傷、生理功能障礙和行為改變[1]。目前人們在呼吁各國聯(lián)合制定與海洋噪聲污染相關(guān)的法律、規(guī)范來保護海洋生物的同時,還從降噪的角度出發(fā),致力于研究削弱噪聲機制以降低污染影響,其中比較成熟的方法是利用氣泡進行水下低頻聲降噪[2?3]。但是氣泡的抗壓性差,并不適用于高壓的水下環(huán)境,因此設(shè)計一種具備耐壓性的水下低頻共鳴器十分必要。

    亥姆霍茲共鳴器[4?5]作為經(jīng)典的聲學(xué)共振結(jié)構(gòu),通常由一個腔體和一個與外界連通的小孔組成,因其結(jié)構(gòu)簡單、聲學(xué)效果顯著而被人們廣泛應(yīng)用于噪聲控制等領(lǐng)域。目前已經(jīng)有不少研究者用它進行改良構(gòu)造成亞波長的完美吸收體[6?9],或者將其按照一定的規(guī)律進行排列組成聲學(xué)超構(gòu)材料[10]來實現(xiàn)對聲波和彈性波的調(diào)控。但是亥姆霍茲共鳴器也存在其固有缺陷[11],比如對頻率的選擇性很強。雖然在共振頻率處,亥姆霍茲共鳴器的消聲效果很好,但工作頻段受品質(zhì)因素影響大,一旦稍微超出范圍,消聲量便會快速下降,同時受到自身結(jié)構(gòu)限制,其低頻消聲效果不理想。特別是在水下環(huán)境中,還要考慮材料特性阻抗、耐壓性等問題,應(yīng)用起來并非易事。

    本文主要基于共振理念提出了在寬開口的剛性腔體內(nèi)嵌入聚氨酯泡沫來構(gòu)造一種新型的水下共鳴器,以實現(xiàn)低頻共振、品質(zhì)因素低的效果。首先是掌握聚氨酯泡沫必要的聲學(xué)參數(shù)。在此基礎(chǔ)上建立鑲嵌聚氨酯泡沫共鳴器的理論模型,根據(jù)聲場形式獲得相應(yīng)的聲阻抗率和共振頻率,并對模型進行簡化檢驗,建立集中參數(shù)系統(tǒng)分析其共振機制,解釋鑲嵌聚氨酯泡沫的共鳴器能夠降低共振頻率、維持較小品質(zhì)因素的原因。最后利用輻射阻抗對理論模型的共振頻率進行補充修正,運用有限元軟件完成建模仿真工作,與理論模型作對比。

    1 聚氨酯泡沫的聲速測定

    聚氨酯泡沫塑料是一種非常受歡迎的吸聲材料,具有質(zhì)量輕、易成型、適用范圍廣等優(yōu)點。要想將聚氨酯泡沫應(yīng)用到聲學(xué)設(shè)計之中,掌握其聲學(xué)特性十分必要。而聚氨酯泡沫塑料泊松比滿足Kerner-Rusch經(jīng)驗公式[12]:

    其中,?=ρavg/ρs是基體材料體積比,ρs和ρavg分別是基體材料和泡沫塑料的密度。根據(jù)文獻[13],假設(shè)E0是聚氨酯基體材料的楊氏模量,那么聚氨酯泡沫楊氏模量Ef的預(yù)測可以采用公式(2):

    聚氨酯泡沫的縱波聲速和橫波聲速為

    已知聚氨酯基體材料的楊氏模量為E0=2433 MPa,密度為ρs=1200 kg·m?3。如果這時知道聚氨酯泡沫的密度,就可以根據(jù)式(1)計算泊松比,再根據(jù)式(2)得到楊氏模量,進而代入式(3)和式(4)獲得相應(yīng)的縱波聲速和橫波聲速,具體的參數(shù)理論估算結(jié)果見表1。

    表1 聚氨酯泡沫參數(shù)理論估算結(jié)果Table 1 Theoretical estimated results for the parameters of polyurethane foam

    因為聚氨酯泡沫的聲速參數(shù)關(guān)系到后續(xù)的理論建模以及將來的實驗研究,所以還需要通過水池實驗測定。在實驗中,可以由兩個水聽器之間有無聚氨酯泡沫板的時間差以及板的厚度來計算聚氨酯泡沫的聲速,具體操作如圖1所示。在實驗過程中,聲源、水聽器的聲中心、聚氨酯泡沫板的幾何中心要在同一水平線上,確保聲波能夠垂直穿過聚氨酯泡沫板。

    圖1 聚氨酯泡沫聲速測量實驗示意圖Fig.1 Schematic diagram of sound velocity measurement for polyurethane foam

    假設(shè)兩個水聽器之間的距離記為l,聚氨酯泡沫板的厚度為d,沒有聚氨酯泡沫板時兩個水聽器接收到聲波的時間差記為t1,則水的聲速有

    當(dāng)聚氨酯泡沫板存在時,兩個水聽器接收到聲波的時間差記為t2,其聲速測量方程可整理為

    其中,d=0.07 m,時間差t1、t2可以通過水聽器信號分別與功放信號做相關(guān)處理后做差獲得。考慮到聲波衍射以及水面散射的影響,兩個水聽器的距離應(yīng)該盡可能接近,所以在實驗過程中l(wèi)主要采用了0.35 m和0.45 m兩個距離。

    為了保證實驗結(jié)果的正確性,選取4 天不同的時間段進行多次測量,具體處理結(jié)果見表2。其中聚氨酯泡沫的縱波聲速可以直接由實驗測量得到。而橫波速度無法直接測量,所以需要先通過式(3)計算出對應(yīng)的楊氏模量,再根據(jù)式(4)獲得。最后對這4 天的測量結(jié)果取平均,聚氨酯泡沫的縱波聲速為970.75 m/s,橫波聲速為495.33 m/s,與理論估算結(jié)果差不多,符合鑲嵌聚氨酯泡沫共鳴器結(jié)構(gòu)設(shè)計的預(yù)期要求。

    表2 聚氨酯泡沫聲速測量處理結(jié)果Table 2 Measured sound velocity of polyurethane foam

    2 建立鑲嵌聚氨酯泡沫的水下共鳴器理論模型

    2.1 理論模型

    求解聲場的理論解法主要有積分方程法和分離變量法[14]。通常認(rèn)為,分離變量法適用于求解有界空間的定解問題,而積分變換法則適用于求解無界空間的定解問題。通常嚴(yán)格理論解只針對規(guī)則目標(biāo),對于現(xiàn)實中比較復(fù)雜的散射體,則需要用到近似解析法和數(shù)值解法。圖2是本文設(shè)計的鑲嵌聚氨酯泡沫共鳴器結(jié)構(gòu),腔壁用鋼制備,中心的圓柱彈性體則為聚氨酯泡沫,其中腔體高度記為h,腔體半徑記為b,鋼壁厚度記為q,聚氨酯泡沫的半徑記為a。在建立理論模型時,假定水為理想流體,聚氨酯泡沫所制成的彈性圓柱體為理想彈性體,腔壁剛性,聲波的傳播絕熱且結(jié)構(gòu)上界面處的聲壓均勻分布。

    圖2 鑲嵌聚氨酯泡沫的水下共鳴器Fig.2 An underwater resonator embedded with polyurethane foam

    已知標(biāo)量勢函數(shù)Φ和矢量勢函數(shù)Ψ可以分別描述固體中的縱波和橫波[15]。因此由氨酯泡沫所制成的彈性體中的聲場滿足關(guān)于標(biāo)量勢和矢量勢的兩個方程:

    彈性體周圍的水中聲場p滿足理想流體中的波動方程:

    其中,c0是水的聲速。因為Laplace算子在柱坐標(biāo)系下的形式為

    而當(dāng)Laplace 算子在柱坐標(biāo)系下作用于一個矢量F時有

    利用Laplace 算子對式(7)、式(8)、式(9)展開,應(yīng)用積分變換法,代入邊界條件即可求解聲場[16?17]。忽略時間因子e?jωt并假定腔體聲場呈軸對稱分布,在z軸方向上對聲場做傅里葉展開,則聚氨酯泡沫彈性棒的縱波、橫波及其周圍水中的聲場可以整理為以下積分形式

    共鳴器聲場模型的邊界條件考慮了聲波在不同媒質(zhì)之間界面上的反射、折射與透射,結(jié)合圖2鑲嵌聚氨酯泡沫共鳴器的結(jié)構(gòu)特點,本文將涉及到的邊界條件整理歸納成表3。根據(jù)表3,只需要將邊界條件逐一代入聲場函數(shù)表達式就可以列出各個系數(shù)之間的方程組,然后進行聯(lián)立求解即可。注意在求解聲場展開系數(shù)的超定方程組時需要用到對角加載奇異值分解方法[18]提高數(shù)值穩(wěn)定性。

    表3 邊界條件Table 3 Boundary conditions

    對于理論模型而言,假定鑲嵌聚氨酯泡沫共鳴器開口處的聲壓是均勻的,則該處質(zhì)點振速分布v(r)的形式可以寫為

    其中,ρ0表示的是水的密度。此時質(zhì)點振速分布與徑向位置有關(guān),對其取平均就可以得到共鳴器開口處的平均質(zhì)點振速:

    如果在邊界條件中設(shè)定開口處的聲壓為1,則聲阻抗率可以寫為

    當(dāng)共鳴器發(fā)生共振時,聲阻抗率的虛部為零,也就是說聲阻抗率虛部為零時所對應(yīng)的頻率值即為共鳴器的共振頻率。

    2.2 聲場模型的簡化檢驗

    為了確保鑲嵌聚氨酯泡沫共鳴器理論模型的正確性,有必要對結(jié)構(gòu)中的聲場模型進行簡化檢驗,所謂簡化就是在建模過程用到的近似處理。如果鑲嵌聚氨酯泡沫共鳴器的理論模型成立,那么它理應(yīng)能夠描述一些在原結(jié)構(gòu)基礎(chǔ)上簡化而來的聲場模式,這里主要考慮了聚氨酯泡沫趨于無以及聚氨酯泡沫趨于填滿整個共鳴腔兩種情況分別用于檢驗水中聲場模式與聚氨酯泡沫彈性體聲場模式。

    首先是水中聲場模式的檢驗。當(dāng)共鳴腔內(nèi)的聚氨酯泡沫半徑趨于無窮小時,模型近似于末端被封閉的有限長直管。而末端被封閉的有限長直管可以直接利用阻抗轉(zhuǎn)移公式得到開口處的聲阻抗率:

    其中,k0=ω/c0。此時將聚氨酯泡沫半徑趨于0 時的理論模型的聲阻抗率與直管結(jié)構(gòu)的聲阻抗率進行對比,若是能夠?qū)?yīng)吻合,則可驗證模型正確。而且為了說明一般性,本文假設(shè)所有的聲速量和密度量均以水為參考進行歸一化,而所有的結(jié)構(gòu)長度量都以腔體半徑作為參考。

    圖3是當(dāng)鑲嵌聚氨酯泡沫共鳴器的腔體半徑b=0.5,聚氨酯泡沫半徑a=0.0001,不同腔體高度的理論模型與對應(yīng)高度的直管結(jié)構(gòu)的聲阻抗率對比結(jié)果,其中星號、圈、四方形、五角星標(biāo)記代表h依次取值為0.3、0.4、0.5、0.6時的理論模型的聲阻抗率,紅色、綠色、玫紅色、黃色實線則為相應(yīng)高度直管的聲阻抗率,二者幾乎對應(yīng)重合。而且通過理論模型計算出來的聲阻抗率的實部均遠遠小于與其相對應(yīng)的虛部,所以可以忽略實部將聲阻抗率看成一個純虛數(shù),這在形式上與式(18)一致。從這兩點來看,水中聲場的建模成立。

    圖3 理論模型與直管模型聲阻抗率對比Fig.3 Comparison of specific acoustic impedance between the theoretical model and a tube model

    然后是聚氨酯泡沫彈性體聲場模式的檢驗。當(dāng)聚氨酯泡沫的半徑趨于腔體半徑時,模型近似于一端固定、一端受到均勻聲壓作用的棒。為了更好地貼合實際棒振動情況,本文從棒的縱振動方程出發(fā),同時考慮徑向運動,得到彈性體頂端處的聲阻抗率為

    其中,kP是聚氨酯泡沫的縱波波數(shù),而Ef=ρavgc2L是聚氨酯泡沫的等效楊氏模量。然后將聚氨酯半徑接近腔體半徑的理論模型的聲阻抗率與棒縱振動的聲阻抗率進行對比,觀察二者是否一致。

    圖4是當(dāng)鑲嵌聚氨酯泡沫共鳴器的腔體半徑b=0.5,聚氨酯泡沫半徑a=0.4999,不同腔體高度的理論模型與對應(yīng)高度的彈性棒縱振動模型的聲阻抗率對比結(jié)果,其中星號、圈、四方形、五角星標(biāo)記代表h依次取值為0.3、0.4、0.5、0.6 時的理論模型的聲阻抗率,紅色、綠色、玫紅色、黃色實線則為相應(yīng)高度彈性棒縱振動的聲阻抗率,二者也幾乎對應(yīng)吻合。而且同樣的,由理論模型計算出來的聲阻抗率實部很小,可以舍去不計,只考慮虛部情況即可,其聲阻抗率在形式上與式(19)一致。說明了聚氨酯泡沫彈性體聲場模式正確。

    圖4 理論模型與棒縱振動模型聲阻抗率對比Fig.4 Comparison of specific acoustic impedance between the theoretical model and the longitudinal vibration model of a rod

    從水中聲場以及聚氨酯泡沫彈性體聲場模型的檢驗結(jié)果來看,鑲嵌聚氨酯泡沫共鳴器的理論模型完全可以描述兩種極端假設(shè)下的聲場情形,且模型的求解結(jié)果與嚴(yán)格的理論值高度吻合。由此可以認(rèn)為,本文所建立的關(guān)于鑲嵌聚氨酯泡沫的共鳴腔內(nèi)聲場的理論模型是正確的。

    2.3 建立鑲嵌聚氨酯泡沫的水下共鳴器的集中參數(shù)模型

    本文在2.2 節(jié)聲場模型的簡化檢驗中討論了直管結(jié)構(gòu)以及棒縱振動的聲阻抗率,證明了理論模型的正確性,但實際上鑲嵌聚氨酯泡沫共鳴器的聲阻抗率應(yīng)該由水和聚氨酯泡沫彈性體共同決定。由于鑲嵌聚氨酯泡沫共鳴器的整體尺寸遠小于本文所關(guān)心頻段的波長,故在描述其聲學(xué)特性時,還能夠采用建立集中參數(shù)等效模型的方法。這時鑲嵌聚氨酯泡沫共鳴器開口處的聲阻抗可以用水中聲場在管口處的聲阻抗與彈性棒縱振動自由端的聲阻抗并聯(lián)的形式來進行等效。

    因為聚氨酯泡沫的存在,開口處的水呈環(huán)形分布,故在管口處水貢獻的聲阻抗為

    聚氨酯泡沫貢獻的聲阻抗為

    由并聯(lián)規(guī)律求出該理論下開口處總的聲阻抗為

    或者用聲阻抗率來表示:

    與簡化檢驗時的做法一樣,繼續(xù)將理論模型和并聯(lián)模型的聲阻抗率進行對比。

    圖5是當(dāng)鑲嵌聚氨酯泡沫共鳴器的腔體半徑b=0.5,腔體高度h=0.5,不同半徑的聚氨酯泡沫的理論模型與并聯(lián)模型聲阻抗率對比結(jié)果,其中星號、圈、四方形標(biāo)記代表a依次取值為0.25、0.35、0.45 時的理論模型的聲阻抗率,紅色、綠色、玫紅色實線則為相應(yīng)并聯(lián)模型的聲阻抗率,二者在數(shù)值以及走勢方面都高度一致,說明了并聯(lián)模型的猜想正確。

    圖5 理論模型與并聯(lián)結(jié)構(gòu)模型聲阻抗率對比Fig.5 Comparison of specific acoustic impedance between the theoretical model and a parallel model

    除此之外建立集中參數(shù)系統(tǒng)也有利于繼續(xù)討論分析內(nèi)嵌聚氨酯泡沫對共鳴器共振機制的影響。所以繼續(xù)對結(jié)構(gòu)開口處的聲阻抗率式(23)進行展開:

    當(dāng)Zst=0 時,得到或者其中n=,1,2,...。由這兩個等式求得兩類共振頻率,分別記為frw和frP,則有

    其中,frw是沒有聚氨酯泡沫時大開口共鳴器的共振頻率,而frP則是純聚氨酯泡沫彈性棒做縱振動時的共振頻率。由于cL

    從另一個角度來看,聚氨酯泡沫的存在可以增大共鳴器的聲容。由式(18)和式(24)可以分別得到無聚氨酯泡沫和有聚氨酯泡沫時共鳴器的聲阻抗率。當(dāng)結(jié)構(gòu)作為低頻散射體存在時,對其做泰勒展開,并只取第一項做近似得到

    進而得到二者的聲容:

    為了比較無聚氨酯泡沫時共鳴器腔體的聲容和有聚氨酯泡沫時共鳴器腔體聲容的大小,將這二者做差?=Cat?Caw并化簡得到

    因為前面在對聚氨酯泡沫做聲速測定的時候,已經(jīng)知道c20ρ0>Ef,所以?>0,嚴(yán)格證明嵌入聚氨酯泡沫確實使得聲容變大了。

    3 模型仿真研究

    前面已經(jīng)建立了鑲嵌聚氨酯泡沫共鳴器的理論模型,并做了大量的理論分析和檢驗,本節(jié)主要應(yīng)用有限元軟件進行建模仿真以實現(xiàn)進一步的驗證,更好地分析鑲嵌聚氨酯泡沫共鳴器的聲學(xué)特性。

    考慮到現(xiàn)實中共鳴器是作為水下散射體存在的,整個水下環(huán)境可以被認(rèn)為是無限大的自由空間,聲場模態(tài)均需要外界入射聲波的激勵才能產(chǎn)生,此時共鳴器中的水和聚氨酯泡沫彈性體會發(fā)生振動,并向外輻射聲能量。輻射抗作為同振質(zhì)量的存在對系統(tǒng)的共振頻率影響很大,所以在描述鑲嵌聚氨酯泡沫共鳴器的散射特性時,需要加上輻射阻抗。假定平面波沿軸向入射,且腔體壁是剛性的,輻射波只能經(jīng)由管口傳出,這與無限大剛性障板上圓形活塞的輻射有一定的類似性,因此可以直接利用活塞的輻射阻抗對所求得的聲阻抗率進行修正。已知無限大剛性障板上活塞的輻射聲阻抗率的實部和虛部[19]分別為

    在這里,r1是活塞的半徑,R1和X1分別是活塞的阻函數(shù)和抗函數(shù),如果令x=2kr1,則它們的具體形式為

    在原來聲阻抗率虛部的基礎(chǔ)上增添一項輻射聲阻抗虛部就能得到鑲嵌聚氨酯泡沫共鳴器的修正聲阻抗率虛部,進而得到修正后的共振頻率。然后將修正后的理論模型與有限元軟件仿真結(jié)果進行對比驗證。

    為了與理論模型對應(yīng)起來,本文在運用有限元軟件建模仿真時同樣也使用了歸一化參數(shù),即所有材料的聲速量和密度量均以水為參考進行歸一化,而所有的結(jié)構(gòu)長度量都以腔體半徑作為參考,具體使用到的參數(shù)見表4。因為聚氨酯泡沫半徑a是本文的設(shè)計變量,故沒有在表4中給出具體的參數(shù)歸一化結(jié)果。

    表4 鑲嵌聚氨酯泡沫共鳴器模型歸一化參數(shù)Table 4 Normalized parameters of the resonator embedded with polyurethane foam in simulations

    為了更加充分地進行模型驗證,根據(jù)表4的參數(shù),再令聚氨酯泡沫半徑a分別為0.0001、0.05、0.1、0.15、0.2、0.25、0.3、0.35、0.4、0.45、0.47 來構(gòu)造不同的鑲嵌聚氨酯泡沫共鳴器,通過有限元軟件建模仿真,計算其共振頻率以及對應(yīng)的聲場分布情況。

    圖6是鑲嵌聚氨酯泡沫共鳴器(a=0.45)歸一化模型在共振頻率0.11578?0.0063008j 下的聲場圖與振型圖,其聲壓和振動位移也分別用最大值進行了歸一化處理,圖例顏色越深代表相對幅值越大。其中0.11578?0.0063008j 的實部表示的是共振頻率的大小,虛部表示結(jié)構(gòu)的聲波輻射能力,當(dāng)虛部很小時可認(rèn)為不具備傳播性。圖6(a)表示鑲嵌聚氨酯泡沫共鳴器在共振頻率下的歸一化聲壓分布,顯然開口處的聲壓呈向外輻射狀態(tài),證明了本文在做理論模型分析時考慮聲輻射并用輻射阻抗對共振頻率進行修正的做法是合理且必要的。同時,高聲壓值主要集中于聚氨酯泡沫彈性棒的頂端,由彈性棒的圓心向四周逐漸減小,說明此時發(fā)生的確實是聚氨酯泡沫彈性棒的軸向共振。而圖6(b)是模型在該共振頻率下的歸一化振動位移狀態(tài),以縱振動模式進行振動,頂端振動幅度最大,證明了聚氨酯泡沫彈性棒確實能夠通過其縱振動降低共鳴器的共振頻率。

    圖6 歸一化的聲壓與總位移(a=0.45) (共振頻率為0.11578?0.0063008j)Fig.6 Normalized pressure distribution and total displacement at the complex frequency 0.11578?0.0063008j for a=0.45

    圖7是鑲嵌聚氨酯泡沫共鳴器擁有不同半徑聚氨酯泡沫時的品質(zhì)因素。品質(zhì)因素越小,共振頻帶越寬,而共鳴器共振頻率的降低往往是以犧牲共振帶寬為代價的。但是從圖7來看,共鳴器嵌入聚氨酯泡沫并沒有使品質(zhì)因素產(chǎn)生大幅度的變化,且數(shù)值基本維持在10 以下。這意味著本文所設(shè)計的鑲嵌聚氨酯泡沫共鳴器在實現(xiàn)低共振頻率的同時依然擁有令人滿意的共振帶寬。

    圖7 鑲嵌聚氨酯泡沫共鳴器的品質(zhì)因素Fig.7 Quality factors of the resonator embedded with polyurethane foam

    圖8則是鑲嵌聚氨酯泡沫共鳴器擁有不同半徑的聚氨酯泡沫時的歸一化共振頻率。帶圈、帶星、帶菱形的3 條曲線分別對應(yīng)了未修正的理論模型、輻射修正后的理論模型以及通過有限元軟件仿真的模型計算得到的共振頻率。當(dāng)聚氨酯泡沫半徑接近0 也就是共鳴器沒有嵌入聚氨酯泡沫的時候,三者計算得到的共振頻率依次大約為0.5、0.323、0.318,與考慮了輻射修正的鑲嵌聚氨酯泡沫共鳴器相比,產(chǎn)生差別的原因是前者沒有考慮輻射阻抗,而后者鋼壁參與了振動。當(dāng)嵌入聚氨酯泡沫以后,沒有經(jīng)過修正的模型獲得的共振頻率始終為0.323,而考慮了輻射修正的理論模型與有限元軟件仿真模型計算的結(jié)果基本一致,且隨著聚氨酯泡沫半徑的增加而呈現(xiàn)降低趨勢,當(dāng)聚氨酯泡沫半徑增加至0.45時,二者的共振頻率分別降低至0.113和0.116,與空共鳴器相比,頻率降低至1/3。但是不可否認(rèn),二者確實還存在著一定的差異,這可以從兩個方面進行解釋。

    圖8 鑲嵌聚氨酯泡沫共鳴器模型與有限元仿真模型的共振頻率對比Fig.8 Resonant frequencies of the resonator embedded with polyurethane foam obtained by the theoretical model and the finite element simulation

    一方面,無限大剛性障板上圓形活塞的輻射阻抗不能完全替代鑲嵌聚氨酯泡沫共鳴器的真實輻射阻抗。本文用來對共振頻率進行修正的輻射阻抗公式是建立在活塞表面振速為均勻分布的假設(shè)下的,而從有限元軟件仿真出來的結(jié)果看,開口處的聲壓并非均勻分布,所以質(zhì)點振速也就不可能滿足均勻分布的條件,導(dǎo)致其實質(zhì)上與活塞輻射沒有嚴(yán)格等價。而開口處振速不均勻最本質(zhì)的原因是聚氨酯泡沫和水各自在開口處的聲阻抗不同。再加上共鳴器的腔體淺,水中的共振基頻又高,因此在做有限元仿真時共鳴器腔體內(nèi)的水可能包含了非平面波模態(tài),而本文的理論模型只考慮了平面波模態(tài)。

    另一方面,鋼壁沒有滿足理論上嚴(yán)格意義的絕對剛要求,在水中會產(chǎn)生輕微振動,這點可以從圖6(b)的振型圖看出來。考慮到現(xiàn)實中人們對水下物體很難做到徹底隔振[20],因此在誤差影響允許的范圍之內(nèi),依然認(rèn)為二者是對應(yīng)一致的,也就是本文建立的鑲嵌聚氨酯泡沫共鳴器的理論模型是有效的。

    4 結(jié)論

    本文設(shè)計了一種鑲嵌聚氨酯泡沫共鳴器,應(yīng)用近似解析法求解腔體內(nèi)聲場建立了理論模型,獲得了鑲嵌聚氨酯泡沫共鳴器的聲阻抗率以及共振頻率的理論表達式。通過與一端封閉的有限長直管結(jié)構(gòu)以及棒縱振動的聲阻抗率對比,證明了理論模型的正確性。在此基礎(chǔ)上,建立集中參數(shù)系統(tǒng),分析鑲嵌聚氨酯泡沫共鳴器共振頻率降低的原因,一是聚氨酯泡沫彈性體的嵌入,使得共振結(jié)構(gòu)的一階共振頻率變成了彈性棒的共振基頻;二是聚氨酯泡沫彈性體的存在增加了共鳴腔的聲容。利用有限元軟件建模仿真,與經(jīng)過輻射修正的理論模型進行對比,兩者的共振頻率相差不大,進一步體現(xiàn)理論模型的正確性,也說明鑲嵌聚氨酯泡沫共鳴器在維持大開口共鳴器品質(zhì)因素小、共振頻帶寬優(yōu)勢的同時實現(xiàn)了共振頻率的降低,為設(shè)計新型水下低頻聲學(xué)共振結(jié)構(gòu)提出了一種新的思路。

    猜你喜歡
    聲阻抗共振頻率聲場
    多激勵下結(jié)構(gòu)反共振頻率的計算和配置
    基于BIM的鐵路車站聲場仿真分析研究
    金屬材料聲發(fā)射信號傳播的聲阻抗特性研究
    基于回波能量的HIFU治療區(qū)聲阻抗測量方法
    探尋360°全聲場發(fā)聲門道
    Multi-functional roles of TaSSI2 involved in Fusarium head blight and powdery mildew resistance and drought tolerance
    三層軟磁鎳超晶格薄膜的共振頻率*
    水管路消聲器聲阻抗測試技術(shù)研究
    懸臂梁共振頻率與其形狀的關(guān)系探討
    機械工程師(2014年6期)2014-11-28 08:16:52
    板結(jié)構(gòu)-聲場耦合分析的FE-LSPIM/FE法
    丁香欧美五月| 757午夜福利合集在线观看| 麻豆一二三区av精品| 精品一区二区三区av网在线观看| 亚洲av电影不卡..在线观看| 国产精品 欧美亚洲| 亚洲欧美精品综合久久99| 老司机午夜十八禁免费视频| 久久性视频一级片| 久久精品成人免费网站| 午夜福利在线观看吧| a级毛片在线看网站| 又大又爽又粗| 麻豆av在线久日| www.www免费av| 桃红色精品国产亚洲av| 日韩高清综合在线| 色精品久久人妻99蜜桃| 麻豆av在线久日| 在线观看午夜福利视频| 狠狠狠狠99中文字幕| 亚洲九九香蕉| 搡老妇女老女人老熟妇| 午夜亚洲福利在线播放| 免费女性裸体啪啪无遮挡网站| 久久精品国产亚洲av高清一级| 国产高清有码在线观看视频 | 熟女电影av网| 麻豆成人av在线观看| 亚洲人成网站在线播放欧美日韩| 欧美不卡视频在线免费观看 | 久久久久国内视频| 人成视频在线观看免费观看| 男女视频在线观看网站免费 | 最近最新免费中文字幕在线| 十八禁人妻一区二区| 久9热在线精品视频| 别揉我奶头~嗯~啊~动态视频| 人人妻人人澡人人看| 国产亚洲欧美98| 午夜福利成人在线免费观看| 亚洲国产精品久久男人天堂| 伦理电影免费视频| 淫妇啪啪啪对白视频| 精品久久久久久久人妻蜜臀av| 一本综合久久免费| 国产精品98久久久久久宅男小说| 国产精品亚洲一级av第二区| 大型黄色视频在线免费观看| 久久久国产欧美日韩av| 日本 av在线| 哪里可以看免费的av片| or卡值多少钱| 大型av网站在线播放| 精品久久久久久久末码| 国产av又大| 在线观看免费午夜福利视频| 中出人妻视频一区二区| 国产久久久一区二区三区| 亚洲天堂国产精品一区在线| 脱女人内裤的视频| 国产伦一二天堂av在线观看| 超碰成人久久| 免费无遮挡裸体视频| 999久久久精品免费观看国产| 一进一出抽搐gif免费好疼| 99久久久亚洲精品蜜臀av| 波多野结衣av一区二区av| 欧美黑人欧美精品刺激| 99热只有精品国产| 国产精品免费一区二区三区在线| 视频区欧美日本亚洲| 免费人成视频x8x8入口观看| 国产av又大| 夜夜夜夜夜久久久久| 国产区一区二久久| 91麻豆av在线| 国产精品 欧美亚洲| 男男h啪啪无遮挡| 777久久人妻少妇嫩草av网站| 变态另类成人亚洲欧美熟女| 午夜视频精品福利| 波多野结衣av一区二区av| 18禁黄网站禁片免费观看直播| 亚洲成人精品中文字幕电影| 免费av毛片视频| 美女高潮到喷水免费观看| 亚洲中文av在线| 欧美黑人巨大hd| 国产亚洲av嫩草精品影院| 久99久视频精品免费| 99久久精品国产亚洲精品| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品国产一区二区精华液| 亚洲国产中文字幕在线视频| 黄色丝袜av网址大全| 国产亚洲精品av在线| 久久国产精品人妻蜜桃| 日日爽夜夜爽网站| 精品久久蜜臀av无| 亚洲国产毛片av蜜桃av| 日韩欧美一区二区三区在线观看| 精品一区二区三区视频在线观看免费| 国内久久婷婷六月综合欲色啪| 日韩大尺度精品在线看网址| 老司机深夜福利视频在线观看| 超碰成人久久| 国产伦一二天堂av在线观看| 国产亚洲精品久久久久久毛片| 男人舔奶头视频| 男人舔女人下体高潮全视频| 国产激情久久老熟女| 老司机福利观看| 夜夜躁狠狠躁天天躁| 国产成+人综合+亚洲专区| 女性被躁到高潮视频| 一级黄色大片毛片| www.自偷自拍.com| 国内揄拍国产精品人妻在线 | 久久久精品欧美日韩精品| 欧美黑人巨大hd| 欧美中文日本在线观看视频| 精品国产乱子伦一区二区三区| 无限看片的www在线观看| 一区二区三区高清视频在线| 精品国产一区二区三区四区第35| АⅤ资源中文在线天堂| 夜夜爽天天搞| 两人在一起打扑克的视频| 成人18禁高潮啪啪吃奶动态图| 又黄又爽又免费观看的视频| 成人精品一区二区免费| 白带黄色成豆腐渣| 波多野结衣高清无吗| 一进一出好大好爽视频| 嫩草影院精品99| 国产精华一区二区三区| 99久久综合精品五月天人人| 两个人看的免费小视频| 免费在线观看视频国产中文字幕亚洲| 国产成人欧美| 成人三级黄色视频| 日韩欧美三级三区| 91成年电影在线观看| 黄片大片在线免费观看| 免费在线观看黄色视频的| 国内久久婷婷六月综合欲色啪| 精品国产国语对白av| 一二三四在线观看免费中文在| av中文乱码字幕在线| 亚洲精品美女久久av网站| 免费搜索国产男女视频| xxxwww97欧美| 看片在线看免费视频| 精品国产亚洲在线| 好看av亚洲va欧美ⅴa在| 非洲黑人性xxxx精品又粗又长| 精品久久久久久成人av| 成年人黄色毛片网站| 欧美日韩亚洲综合一区二区三区_| 亚洲国产精品sss在线观看| 2021天堂中文幕一二区在线观 | 久久久久久人人人人人| 美国免费a级毛片| 亚洲一区二区三区不卡视频| 美女扒开内裤让男人捅视频| 免费看日本二区| 欧美一级a爱片免费观看看 | 亚洲精品一卡2卡三卡4卡5卡| 日韩欧美 国产精品| 精品乱码久久久久久99久播| 国产精品久久久久久人妻精品电影| 欧美成人免费av一区二区三区| 国产亚洲精品av在线| 老司机靠b影院| 99在线人妻在线中文字幕| 国产成人精品久久二区二区91| 欧美日韩福利视频一区二区| 男人操女人黄网站| 91麻豆av在线| 亚洲成av人片免费观看| 国产亚洲精品久久久久5区| 久久精品国产亚洲av香蕉五月| 一二三四在线观看免费中文在| 成人一区二区视频在线观看| 手机成人av网站| 97碰自拍视频| 欧美黑人欧美精品刺激| 丝袜人妻中文字幕| 性色av乱码一区二区三区2| 日韩高清综合在线| 国产精品永久免费网站| 老熟妇乱子伦视频在线观看| www.自偷自拍.com| 一区福利在线观看| 亚洲成av片中文字幕在线观看| 欧美不卡视频在线免费观看 | 国产精品乱码一区二三区的特点| 老熟妇仑乱视频hdxx| 久久久久久久久久黄片| 黄片大片在线免费观看| 天堂影院成人在线观看| 欧美国产精品va在线观看不卡| 成人av一区二区三区在线看| 久久精品国产亚洲av香蕉五月| 男女之事视频高清在线观看| 可以在线观看的亚洲视频| 观看免费一级毛片| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产高清在线一区二区三 | 久久精品成人免费网站| 亚洲男人的天堂狠狠| 欧美乱色亚洲激情| 人人妻,人人澡人人爽秒播| 女生性感内裤真人,穿戴方法视频| 视频在线观看一区二区三区| 欧美国产日韩亚洲一区| 欧美绝顶高潮抽搐喷水| 国产精品久久久久久亚洲av鲁大| 亚洲av电影在线进入| 亚洲成人久久爱视频| 欧美一级a爱片免费观看看 | 亚洲精品中文字幕在线视频| 免费在线观看完整版高清| 亚洲国产精品成人综合色| 久久精品国产99精品国产亚洲性色| 欧美乱码精品一区二区三区| 欧美日本亚洲视频在线播放| 99国产综合亚洲精品| 亚洲人成电影免费在线| 国产在线精品亚洲第一网站| www日本在线高清视频| 午夜福利在线在线| 久久香蕉国产精品| 国产色视频综合| 欧美大码av| 亚洲中文字幕日韩| av欧美777| 亚洲九九香蕉| 一二三四在线观看免费中文在| 大型黄色视频在线免费观看| 亚洲成人国产一区在线观看| 亚洲国产欧美一区二区综合| 亚洲精品国产一区二区精华液| 少妇熟女aⅴ在线视频| 午夜免费成人在线视频| 精品久久久久久久毛片微露脸| 午夜老司机福利片| 一级毛片精品| 国产爱豆传媒在线观看 | 淫秽高清视频在线观看| 99精品久久久久人妻精品| 精品第一国产精品| 欧美三级亚洲精品| 亚洲久久久国产精品| 香蕉国产在线看| 国产精品电影一区二区三区| 麻豆久久精品国产亚洲av| 国产91精品成人一区二区三区| 久久精品成人免费网站| 高清在线国产一区| 国语自产精品视频在线第100页| 亚洲国产看品久久| 久久亚洲真实| 国产人伦9x9x在线观看| 亚洲五月色婷婷综合| 男人舔女人下体高潮全视频| 一级片免费观看大全| www.精华液| 亚洲国产高清在线一区二区三 | 中亚洲国语对白在线视频| √禁漫天堂资源中文www| 国产精品一区二区精品视频观看| 一区二区日韩欧美中文字幕| 久久久久久久精品吃奶| 国产精品 欧美亚洲| 在线观看66精品国产| 国产私拍福利视频在线观看| 久久久久九九精品影院| 夜夜躁狠狠躁天天躁| 国产高清视频在线播放一区| 琪琪午夜伦伦电影理论片6080| 黄色片一级片一级黄色片| 99久久99久久久精品蜜桃| 久久人人精品亚洲av| 亚洲av熟女| 午夜久久久在线观看| 国产精品1区2区在线观看.| 日韩精品免费视频一区二区三区| 久久久久九九精品影院| 国产激情欧美一区二区| 日本免费一区二区三区高清不卡| 黄色视频不卡| 欧美日韩福利视频一区二区| 热re99久久国产66热| 村上凉子中文字幕在线| 狠狠狠狠99中文字幕| 国内久久婷婷六月综合欲色啪| 一边摸一边做爽爽视频免费| 天堂动漫精品| 精品熟女少妇八av免费久了| netflix在线观看网站| 色在线成人网| www.www免费av| 一级毛片高清免费大全| 午夜福利在线在线| 免费看美女性在线毛片视频| 国产激情久久老熟女| 又紧又爽又黄一区二区| 中文字幕人妻丝袜一区二区| 午夜免费激情av| 99国产精品99久久久久| 国产男靠女视频免费网站| 制服丝袜大香蕉在线| 少妇的丰满在线观看| 女性被躁到高潮视频| 精品乱码久久久久久99久播| av福利片在线| 国产亚洲精品av在线| 在线视频色国产色| 色播亚洲综合网| 淫妇啪啪啪对白视频| 久久欧美精品欧美久久欧美| 90打野战视频偷拍视频| 最新在线观看一区二区三区| 日韩精品中文字幕看吧| 日本 av在线| 久久久久国内视频| 国产99久久九九免费精品| 2021天堂中文幕一二区在线观 | 岛国在线观看网站| 久久中文看片网| 欧美人与性动交α欧美精品济南到| 麻豆一二三区av精品| 美女午夜性视频免费| 18美女黄网站色大片免费观看| 啦啦啦观看免费观看视频高清| 日韩精品青青久久久久久| 搡老妇女老女人老熟妇| 丰满人妻熟妇乱又伦精品不卡| 99久久久亚洲精品蜜臀av| 国产精品自产拍在线观看55亚洲| 欧美日韩亚洲国产一区二区在线观看| 在线观看www视频免费| 久久久精品欧美日韩精品| 色av中文字幕| 国产亚洲精品久久久久久毛片| 欧美激情极品国产一区二区三区| 欧美黄色淫秽网站| 国产乱人伦免费视频| 亚洲成人国产一区在线观看| 国内精品久久久久久久电影| 欧美日韩精品网址| 他把我摸到了高潮在线观看| 99riav亚洲国产免费| 亚洲天堂国产精品一区在线| 久久久久久九九精品二区国产 | 亚洲第一av免费看| 欧美黑人巨大hd| 特大巨黑吊av在线直播 | 亚洲欧美一区二区三区黑人| 日本三级黄在线观看| 久久精品91蜜桃| 亚洲精品粉嫩美女一区| 老司机午夜十八禁免费视频| 白带黄色成豆腐渣| 国产91精品成人一区二区三区| 操出白浆在线播放| 精品久久久久久久人妻蜜臀av| 午夜影院日韩av| 母亲3免费完整高清在线观看| 中文字幕最新亚洲高清| 亚洲精品在线美女| 成人亚洲精品一区在线观看| 国产极品粉嫩免费观看在线| 欧美久久黑人一区二区| 一边摸一边做爽爽视频免费| 亚洲熟妇中文字幕五十中出| 97超级碰碰碰精品色视频在线观看| 熟妇人妻久久中文字幕3abv| 一区福利在线观看| 久久香蕉激情| 国产伦在线观看视频一区| 国产真人三级小视频在线观看| 在线观看免费日韩欧美大片| 精品高清国产在线一区| 他把我摸到了高潮在线观看| 91国产中文字幕| 可以在线观看毛片的网站| 欧美+亚洲+日韩+国产| 国产成人精品无人区| 黄色a级毛片大全视频| 好男人电影高清在线观看| 日本a在线网址| 欧美成人免费av一区二区三区| 动漫黄色视频在线观看| 亚洲精品美女久久av网站| 热99re8久久精品国产| 精品国产国语对白av| 欧美av亚洲av综合av国产av| 99国产精品一区二区蜜桃av| 淫妇啪啪啪对白视频| 午夜福利高清视频| 国产成人精品无人区| 中文字幕久久专区| 亚洲国产精品999在线| 丝袜人妻中文字幕| 精品第一国产精品| 久久久久免费精品人妻一区二区 | 亚洲男人天堂网一区| 久久这里只有精品19| 欧美日韩精品网址| 女生性感内裤真人,穿戴方法视频| 久久久国产成人精品二区| 999久久久精品免费观看国产| av在线天堂中文字幕| 久久国产精品人妻蜜桃| 在线视频色国产色| 亚洲人成电影免费在线| 一二三四在线观看免费中文在| 免费观看人在逋| 成熟少妇高潮喷水视频| 久久久久久九九精品二区国产 | 欧美激情 高清一区二区三区| 精品日产1卡2卡| 此物有八面人人有两片| 国产精品 国内视频| 看片在线看免费视频| 美女午夜性视频免费| 免费在线观看黄色视频的| 美女国产高潮福利片在线看| 婷婷丁香在线五月| 久久 成人 亚洲| 久热爱精品视频在线9| 亚洲成人精品中文字幕电影| 午夜福利在线在线| 一个人免费在线观看的高清视频| 国产av又大| 亚洲精品久久国产高清桃花| 国产精品亚洲av一区麻豆| 免费高清在线观看日韩| 午夜亚洲福利在线播放| 中文字幕人成人乱码亚洲影| 国产精品免费一区二区三区在线| 国产精品久久久久久人妻精品电影| 精品高清国产在线一区| 一区二区日韩欧美中文字幕| 欧美黄色淫秽网站| 日韩免费av在线播放| 亚洲精品久久国产高清桃花| 亚洲专区字幕在线| 天天添夜夜摸| 日韩欧美国产在线观看| 日本 av在线| 成人三级做爰电影| 男女午夜视频在线观看| 三级毛片av免费| 又黄又爽又免费观看的视频| 男女之事视频高清在线观看| 久久久久免费精品人妻一区二区 | 欧美日韩一级在线毛片| 老司机深夜福利视频在线观看| 亚洲无线在线观看| 首页视频小说图片口味搜索| 久久欧美精品欧美久久欧美| 人人妻人人澡人人看| ponron亚洲| 天天添夜夜摸| 国产久久久一区二区三区| 亚洲国产精品999在线| 桃色一区二区三区在线观看| 少妇粗大呻吟视频| 久久久久久久久久黄片| 特大巨黑吊av在线直播 | 一本久久中文字幕| 国产伦在线观看视频一区| 国产精品久久久久久精品电影 | 午夜视频精品福利| 又大又爽又粗| 婷婷精品国产亚洲av| 色精品久久人妻99蜜桃| 日韩精品中文字幕看吧| 大型av网站在线播放| 亚洲午夜理论影院| 欧美黑人欧美精品刺激| 精品第一国产精品| 午夜福利欧美成人| 精品欧美一区二区三区在线| 亚洲真实伦在线观看| 中文字幕精品亚洲无线码一区 | 19禁男女啪啪无遮挡网站| 啦啦啦韩国在线观看视频| 久热这里只有精品99| 非洲黑人性xxxx精品又粗又长| 婷婷精品国产亚洲av在线| 久久久精品欧美日韩精品| 精品久久久久久久久久久久久 | 十八禁人妻一区二区| 一本大道久久a久久精品| 99国产综合亚洲精品| 亚洲国产欧美网| 亚洲狠狠婷婷综合久久图片| 啦啦啦韩国在线观看视频| 久久久水蜜桃国产精品网| 很黄的视频免费| 午夜视频精品福利| 日本一区二区免费在线视频| 日本免费一区二区三区高清不卡| 一级作爱视频免费观看| 欧美日韩精品网址| 中亚洲国语对白在线视频| 国产精品亚洲一级av第二区| 99久久无色码亚洲精品果冻| 999久久久精品免费观看国产| 叶爱在线成人免费视频播放| 国产成+人综合+亚洲专区| 久久精品国产亚洲av高清一级| 亚洲av五月六月丁香网| 人成视频在线观看免费观看| 精品久久久久久久毛片微露脸| 最新美女视频免费是黄的| 欧美在线黄色| 成年女人毛片免费观看观看9| 国产一区在线观看成人免费| 国产97色在线日韩免费| 午夜免费成人在线视频| 久久久精品国产亚洲av高清涩受| 亚洲精品一区av在线观看| 国产色视频综合| 国产精品综合久久久久久久免费| 国产极品粉嫩免费观看在线| 午夜福利高清视频| 18禁黄网站禁片午夜丰满| 999久久久国产精品视频| 禁无遮挡网站| 久热爱精品视频在线9| 18禁裸乳无遮挡免费网站照片 | 超碰成人久久| 日韩欧美一区视频在线观看| 欧美色欧美亚洲另类二区| 午夜免费成人在线视频| 欧美人与性动交α欧美精品济南到| 免费在线观看影片大全网站| 午夜精品久久久久久毛片777| 嫩草影视91久久| 亚洲人成伊人成综合网2020| 国产亚洲精品第一综合不卡| 日韩一卡2卡3卡4卡2021年| 欧美日本亚洲视频在线播放| 首页视频小说图片口味搜索| www日本在线高清视频| 亚洲精华国产精华精| 亚洲中文字幕一区二区三区有码在线看 | 日韩一卡2卡3卡4卡2021年| 亚洲自拍偷在线| 婷婷丁香在线五月| 日韩欧美一区视频在线观看| 黄片播放在线免费| 国产成年人精品一区二区| 中文字幕av电影在线播放| 天天添夜夜摸| www日本在线高清视频| 久久久久国内视频| 亚洲国产欧美日韩在线播放| 别揉我奶头~嗯~啊~动态视频| 日韩欧美国产一区二区入口| 亚洲成国产人片在线观看| 亚洲精品久久成人aⅴ小说| 91大片在线观看| www.999成人在线观看| 色播在线永久视频| 久久精品国产亚洲av高清一级| 在线观看66精品国产| 1024手机看黄色片| 18禁观看日本| 亚洲成国产人片在线观看| 亚洲精品久久成人aⅴ小说| 成人18禁高潮啪啪吃奶动态图| 长腿黑丝高跟| 国产成+人综合+亚洲专区| 免费搜索国产男女视频| 精品福利观看| 午夜精品在线福利| 一边摸一边做爽爽视频免费| 好男人在线观看高清免费视频 | 精品欧美国产一区二区三| 国产高清视频在线播放一区| 十分钟在线观看高清视频www| 男女下面进入的视频免费午夜 | 每晚都被弄得嗷嗷叫到高潮| 中文字幕人妻熟女乱码| 男女午夜视频在线观看| √禁漫天堂资源中文www| 亚洲国产欧洲综合997久久, | 俄罗斯特黄特色一大片| 欧美 亚洲 国产 日韩一| 国产成年人精品一区二区| 真人做人爱边吃奶动态| 久久精品国产亚洲av高清一级| 精品少妇一区二区三区视频日本电影| 久久久久九九精品影院| 久久久久久九九精品二区国产 | 久久伊人香网站| 色综合亚洲欧美另类图片| 怎么达到女性高潮| 免费看美女性在线毛片视频| 在线观看66精品国产| 给我免费播放毛片高清在线观看| 日韩欧美国产一区二区入口| 熟女少妇亚洲综合色aaa.| www.999成人在线观看| 9191精品国产免费久久| 99久久无色码亚洲精品果冻| 亚洲成人久久性| 亚洲成国产人片在线观看|