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

    管道有源噪聲控制中壁面分布次級(jí)聲源的空間分布優(yōu)化

    2022-11-21 01:11:06羅平展張芳杰徐健李曉東
    應(yīng)用聲學(xué) 2022年5期
    關(guān)鍵詞:聲源聲波高階

    羅平展 張芳杰 徐健 李曉東?

    (1 中國(guó)科學(xué)院噪聲與振動(dòng)重點(diǎn)實(shí)驗(yàn)室(聲學(xué)研究所)北京 100190)

    (2 中國(guó)科學(xué)院大學(xué)北京 100049)

    0 引言

    有源噪聲控制技術(shù)已成功應(yīng)用于管道中的低頻噪聲控制。相比于控制管道中的平面波,高階模式聲波的控制要求更高,是有源控制的難點(diǎn)問(wèn)題[1]。由于控制高階模式聲波需要多個(gè)換能器,其中一個(gè)難點(diǎn)是設(shè)計(jì)次級(jí)聲源和誤差傳聲器等換能器在管道中的空間布放。已有文獻(xiàn)對(duì)控制管道中的高階模式聲波進(jìn)行了一些探索,但是換能器的最優(yōu)空間分布尚未得到闡明[2-3]??紤]到實(shí)際應(yīng)用中諸如風(fēng)洞和地鐵等大型通風(fēng)管道,為了保障暢通,要求換能器分布在管道側(cè)壁[4-5]。因此,本文主要研究次級(jí)聲源在被約束于管道壁面條件下的空間分布優(yōu)化問(wèn)題。

    管道中次級(jí)聲源常被布放于管道截面。由于高階模式主要由管道截面特性所決定,當(dāng)次級(jí)聲源合理分布在截面時(shí),能有效控制各高階模式聲波。理論上使用N個(gè)次級(jí)聲源最多可以控制前N個(gè)模式聲波,即工作頻率上限逼近第N+1階模式聲波的截止頻率[2]。已有的研究工作表明,次級(jí)聲源布放在管道截面時(shí)能夠有效控制工作頻率范圍內(nèi)的高階模式聲波,阻礙其通過(guò)截面?zhèn)鞑ブ凉艿老掠危瑥亩玫絻?yōu)異的全局控制效果[3,6-7]。這種位于截面的布放方式應(yīng)用廣泛,對(duì)于在管道和腔體等出口處抑制噪聲輻射同樣有效[8]。因此,將次級(jí)聲源布放在管道截面上成為有源噪聲控制中的常見(jiàn)方案。

    然而在部分應(yīng)用中,由于實(shí)際條件限制使得次級(jí)聲源被約束在管道壁面。此時(shí)噪聲控制效果相比于次級(jí)聲源位于管道截面時(shí)通常有所劣化,表現(xiàn)出降噪量下降且工作頻率范圍顯著低于理論最優(yōu)值等現(xiàn)象,更加凸顯了優(yōu)化設(shè)計(jì)管道壁面上次級(jí)聲源空間分布的重要性[9]。已有研究工作針對(duì)壁面分布次級(jí)聲源的空間分布優(yōu)化問(wèn)題的研究尚不充分。Stell等[2]指出次級(jí)聲源的布放應(yīng)當(dāng)避開(kāi)高階模式聲波的節(jié)線,以保證次級(jí)聲源能較好地激勵(lì)高階模式聲波。Zander等[10]針對(duì)2個(gè)次級(jí)聲源的簡(jiǎn)單情形討論了軸向間距的選擇。Laugesen[11]考慮到在剛性壁假設(shè)下所有高階模式聲波在管道角落都能取到最大值,建議將次級(jí)聲源布放在管道角落處以提高控制效率。在仿真實(shí)驗(yàn)中他還發(fā)現(xiàn)將次級(jí)聲源沿管道軸向展開(kāi),可以增加系統(tǒng)控制的魯棒性。Zhang等[12]對(duì)比了次級(jí)聲源與誤差傳聲器的徑向分布,發(fā)現(xiàn)當(dāng)次級(jí)聲源和誤差傳聲器的數(shù)量較少時(shí),在相對(duì)的平行壁面上布放它們可以更好地控制高階模式聲波。但是這些研究工作由于缺乏詳細(xì)的理論分析,沒(méi)有得到系統(tǒng)性的完整結(jié)論,難以成為次級(jí)聲源空間分布優(yōu)化設(shè)計(jì)的具體指導(dǎo)準(zhǔn)則。迄今為止,次級(jí)聲源被約束在管道壁面時(shí)對(duì)于控制高階模式聲波的局限性,以及優(yōu)化設(shè)計(jì)次級(jí)聲源空間分布的原理和方法均未得到清晰徹底的解釋。

    針對(duì)管道中高階模式聲波的有源控制問(wèn)題,本文研究了次級(jí)聲源在管道壁面的空間分布優(yōu)化問(wèn)題。本文首先介紹了管道中次級(jí)聲源激勵(lì)高階模式聲波的理論模型;然后顯式推導(dǎo)了次級(jí)聲源空間分布對(duì)控制性能的影響,分析了次級(jí)聲源空間分布的優(yōu)化準(zhǔn)則;最后在最優(yōu)控制策略下,通過(guò)數(shù)值仿真,驗(yàn)證了次級(jí)聲源空間分布優(yōu)化準(zhǔn)則的有效性。

    1 理論模型

    1.1 管道中次級(jí)聲源激勵(lì)高階模式聲波建模

    不失一般性,假設(shè)一無(wú)限長(zhǎng)剛性壁面的矩形管道,其截面為S=Lx×Ly。次級(jí)聲源被約束在管道壁面,如圖1所示。建立直角坐標(biāo)系,約定z軸方向?yàn)楣艿垒S向,xy方向?yàn)楣艿乐芟?。將次?jí)聲源根據(jù)其振速分布建模為W(r′),其中r′=(x′,y′,z′)表示聲源位置。矩形管道中的格林函數(shù)為G(r|r′),其中r=(x,y,z)表示觀測(cè)點(diǎn)位置。為了表述簡(jiǎn)便,在下文的推導(dǎo)中省去時(shí)間項(xiàng)eiωt。

    圖1 壁面分布次級(jí)聲源示意圖Fig.1 The sketch of the boundary-located secondary sources

    管道中次級(jí)聲源所激勵(lì)聲場(chǎng)聲壓根據(jù)格林函數(shù)可以表示成[13]

    其中,ω為激勵(lì)頻率,ρ0為空氣密度,格林函數(shù)G(r|r′)可以表示為

    其中,φn表示管道的高階模式聲波。對(duì)于四周為剛性壁面的矩形管道,φn滿足

    對(duì)應(yīng)的歸一化系數(shù)為

    第n階模式聲波的波數(shù)k在各方向上的分量可以分別表示為

    下文開(kāi)始討論壁面分布次級(jí)聲源的空間分布對(duì)控制管道高階模式聲波的影響。假設(shè)次級(jí)聲源的分布范圍限于0≤z≤Lz范圍內(nèi)的4個(gè)壁面。構(gòu)造次級(jí)聲源分布Wm(r′)滿足

    其中,δ()表示狄拉克函數(shù)的定義與式(4)中相同,即等于第m階模式聲波各方向上的波數(shù)。顯而易見(jiàn),Wm(r′)代表了在管道壁面上空間頻率與第m階模式聲波的頻率相同的4個(gè)面聲源分別是這4個(gè)面聲源的強(qiáng)度。根據(jù)傅里葉變換的原理,壁面上任意有限的次級(jí)聲源的空間分布均能表示成不同Wm(r′)的疊加[14]。因此,通過(guò)研究Wm(r′)代表的次級(jí)聲源分布對(duì)高階模式聲波的激勵(lì)機(jī)制,可以推演得到任意有限分布的次級(jí)聲源對(duì)高階模式聲波的激勵(lì)。

    不失一般性,只考慮向管道下游傳播的高階模式聲波,管道內(nèi)的聲場(chǎng)可以分解成模式聲波疊加的形式[13],即

    次級(jí)聲源Wm(r′)在管道中所激勵(lì)的高階模式聲波,可以通過(guò)將式(6)代入式(1)中得到,對(duì)比式(7),得到Wm(r′)所激勵(lì)的各階模式聲波的系數(shù)為

    式(8)中包含了x、y和z三個(gè)方向的積分,彼此獨(dú)立,分別計(jì)算如下。

    首先,次級(jí)聲源在x方向空間分布對(duì)于激勵(lì)高階模式聲波的貢獻(xiàn)可以寫(xiě)成

    其中,δmn為克羅內(nèi)克函數(shù),它的值僅當(dāng)m=n為1,其他時(shí)候均為0。式(9)中x方向的積分表明,在x方向按分布的次級(jí)聲源,只激勵(lì)mx相關(guān)的高階模式聲波,不會(huì)對(duì)其他模式聲波產(chǎn)生影響。同時(shí),管道中x方向上兩個(gè)相對(duì)的平行墻面上次級(jí)聲源的貢獻(xiàn)分別為X1和X-1。可以令,則對(duì)于my為奇數(shù)的模式聲波,X1和X-1反向,X為0,此時(shí)只激勵(lì)my為偶數(shù)的模式聲波;同理,可以令,此時(shí)只激勵(lì)my為奇數(shù)的模式聲波。這可以解釋為,兩個(gè)沿x方向沿伸的平行墻面等同于y方向的兩點(diǎn)空間采樣,可以控制my相關(guān)的兩個(gè)模態(tài),使得根據(jù)my的奇偶性可以區(qū)分和獨(dú)立控制相關(guān)模式聲波。

    同理可得次級(jí)聲源在y方向的空間分布對(duì)于激勵(lì)高階模式聲波的貢獻(xiàn)為

    式(10)中y方向的積分表明,在y方向按kmy分布的次級(jí)聲源,只激勵(lì)my相關(guān)的高階模式聲波,不會(huì)對(duì)其他模式聲波產(chǎn)生影響。同時(shí),可以通過(guò)選擇amx和a′mx的符號(hào)來(lái)決定Y1和Y-1的值,以選擇性的激勵(lì)mx相關(guān)的奇數(shù)階或偶數(shù)階模式聲波。

    最后,次級(jí)聲源在z方向的空間分布對(duì)于激勵(lì)高階模式聲波的貢獻(xiàn)可以寫(xiě)成

    其中,sinc(x)=sinx/x。式(11)表明壁面分布次級(jí)聲源無(wú)法對(duì)z方向上的模式聲波進(jìn)行嚴(yán)格的獨(dú)立控制,必然會(huì)激勵(lì)各種高階模式聲波,這是壁面分布次級(jí)聲源控制高階模式聲波局限性的體現(xiàn)。但是應(yīng)注意到,式(11)中關(guān)于Lz/2只有一個(gè)線性項(xiàng),其他均為隨Lz的震蕩衰減的sinc函數(shù)。因此,當(dāng)Lz增大時(shí),線性項(xiàng)增大從而顯著大于其他sinc函數(shù),使得在z方向上按照kmz分布的次級(jí)聲源在軸向的尺寸Lz較大時(shí)能顯著激勵(lì)z方向上波數(shù)為kmz的模式聲波。進(jìn)一步考慮到,當(dāng)x>π時(shí),有sinc(x)?1,因此一般要求壁面分布次級(jí)聲源的軸向尺寸Lz>2π/min(knz-kmz),其中min()表示在各種nz和mz的組合中選擇波數(shù)差的最小值。

    將式(9)、式(10)和式(11)代入式(8),得到壁面分布次級(jí)聲源激勵(lì)的高階模式聲波可以寫(xiě)成

    式(12)表明,壁面分布次級(jí)聲源在x方向和y方向,即在管道周向的空間分布對(duì)激勵(lì)高階模式聲波的貢獻(xiàn)會(huì)直接線性相加,再與z方向即管道軸向上空間分布的貢獻(xiàn)相乘,共同激勵(lì)高階模式聲波。當(dāng)次級(jí)聲源在管道壁面上按照Wm(r′)分布時(shí),主要激勵(lì)第m階模式聲波;同理,第m階模式聲波也主要由Wm(r′)分布的次級(jí)聲源激勵(lì)。次級(jí)聲源的空間分布頻率和激勵(lì)高級(jí)模式聲波存在對(duì)應(yīng)關(guān)系,是對(duì)高階模式獨(dú)立控制的理論依據(jù)。為了形象展示壁面分布次級(jí)聲源在各方向上的空間分布對(duì)于控制高階模式聲波的貢獻(xiàn),圖2給出了X、Y和Z對(duì)激勵(lì)第m階(mx,my)模式聲波的貢獻(xiàn)示意圖,其中第m階模式聲波的幅度被歸一化到1以方便比較。圖2表明,壁面分布次級(jí)聲源的空間分布對(duì)其控制性能非常重要。當(dāng)次級(jí)聲源在合理優(yōu)化空間分布后能近似獨(dú)立地控制管道中各高階模式聲波,從而實(shí)現(xiàn)較好的控制效果;但是若未被合理的布放,在激勵(lì)目標(biāo)模式聲波時(shí)也將同時(shí)激勵(lì)其他模式聲波,無(wú)法實(shí)現(xiàn)獨(dú)立激勵(lì)各高階模式聲波,從而喪失對(duì)高階模式的控制效果。

    圖2 壁面分布次級(jí)聲源Wm(r′)激勵(lì)(mx,my)模式聲波的示意圖Fig.2 Illustration of the mechanism of generating the(mx,my)mode by boundary-located secondary sources Wm(r′)

    1.2 點(diǎn)聲源控制策略

    在實(shí)際應(yīng)用中,常以揚(yáng)聲器作為有源控制系統(tǒng)的次級(jí)聲源。由于其尺寸相對(duì)聲波波長(zhǎng)較小,此時(shí)次級(jí)聲源應(yīng)被建模成點(diǎn)聲源。點(diǎn)聲源可以看作聲源連續(xù)分布W(r′)的空間采樣,因此在滿足采樣定律的條件下前文所得關(guān)于次級(jí)聲源空間分布的結(jié)論同樣適用。次級(jí)聲源最優(yōu)驅(qū)動(dòng)強(qiáng)度可以通過(guò)以最小化管道中聲能流為目標(biāo)的全局控制策略得到。本小節(jié)介紹控制策略的推導(dǎo),以便后文對(duì)比次級(jí)聲源空間分布對(duì)控制效果的影響。

    為了表述方便,定義管道中的N個(gè)傳播模式聲波為a=[a1,a2,···,aN]T,M個(gè)次級(jí)點(diǎn)聲源的驅(qū)動(dòng)強(qiáng)度為q=[q1,q2,···,qM]T,對(duì)應(yīng)的空間位置為對(duì)式(8)在上述位置進(jìn)行空間采樣,受激模式聲波與次級(jí)聲源的關(guān)系可以用矩陣的形式來(lái)表示:

    定義N×M維阻抗矩陣Ω(r′),其各元素為

    根據(jù)式(14)的定義,將式(13)重寫(xiě)為

    下面討論以最小化管道中聲能流為目標(biāo)的全局控制策略以及次級(jí)聲源最優(yōu)驅(qū)動(dòng)強(qiáng)度的求解[3]。管道中聲能流可以表示為

    其中,Re()表示取復(fù)數(shù)的實(shí)數(shù)部分,Uz(r)是聲場(chǎng)在z方向上的粒子振速,*表示共軛。注意只有傳播模式聲波能傳遞聲能流。

    定義一個(gè)N×N的對(duì)角矩陣,表示為

    因此,式(16)中的聲能流可以寫(xiě)成

    選用下標(biāo)p和s分別表示與初級(jí)聲源和次級(jí)聲源相關(guān)的參量。初級(jí)聲源的驅(qū)動(dòng)強(qiáng)度、阻抗矩陣和激勵(lì)的模式聲波分別記作qp、Ωp和ap;同理,次級(jí)聲源的驅(qū)動(dòng)強(qiáng)度、阻抗矩陣和激勵(lì)的模式聲波分別記作qs、Ωs和αs。選擇代價(jià)函數(shù)J為管道中初級(jí)聲源和次級(jí)聲源共同激勵(lì)的總聲能流與表示次級(jí)聲源強(qiáng)度的懲罰項(xiàng)之和[2],寫(xiě)成

    其中,β為正實(shí)數(shù),用來(lái)限制次級(jí)聲源的驅(qū)動(dòng)強(qiáng)度大小[15]。取代價(jià)函數(shù)J對(duì)qs的導(dǎo)數(shù)為零,可得最優(yōu)次級(jí)聲源驅(qū)動(dòng)強(qiáng)度qopt為

    2 仿真實(shí)驗(yàn)與分析

    下文通過(guò)數(shù)值仿真對(duì)比壁面分布次級(jí)聲源在前文所述控制策略下的各種空間分布對(duì)高階模式控制效果的影響,以驗(yàn)證本文所提出的理論模型和結(jié)論的有效性。

    仿真實(shí)驗(yàn)中無(wú)限長(zhǎng)矩形管道的截面積設(shè)為4.0 m×3.0 m(Lx×Ly)。管道中初級(jí)聲場(chǎng)由布放在(0.0,0.0,-10.0)m處的點(diǎn)聲源所激勵(lì),仿真中聲波的頻率范圍設(shè)為30~210 Hz。在這個(gè)頻率范圍內(nèi),最多存在18個(gè)傳播模式聲波,對(duì)應(yīng)的截止頻率為fn=c/2[(nx/Lx)2+(ny/Ly)2]1/2。它們的階數(shù)和截止頻率的具體取值如表1所示。

    表1 管道中各階模式聲波的截止頻率Table 1 The higher-order modes in the duct and the corresponding cut-on frequencies

    對(duì)于前18階模式聲波,其階數(shù)滿足nx≤4和ny≤3,即在x方向上不超過(guò)5個(gè)模式聲波,在y方向上不超過(guò)4個(gè)模式聲波。因此根據(jù)采樣定律,在x方向上和y方向上分別需要5個(gè)和4個(gè)次級(jí)聲源。在z方向上對(duì)于可能的多層分布次級(jí)聲源,其層間距和軸向尺寸分別記為ΔLz和Lz。為了方便比較,壁面分布次級(jí)聲源在各方向上都服從均勻分布。仿真中β始終設(shè)為矩陣最大特征值的1/1000以保證矩陣求逆魯棒性。為了與前文所述最小化管道聲能流的控制策略一致,本文采用管道中殘余的聲能流作為控制效果的評(píng)價(jià)標(biāo)準(zhǔn),并將得到的聲能流做歸一化處理,其中0 dB對(duì)應(yīng)1.0×10-12W。

    圖3給出了x和y方向上壁面分布次級(jí)聲源的空間分布示意圖和相應(yīng)的控制效果。其中,空間分布(1)~(4)分別對(duì)應(yīng)前文中的X1、X、Y1和Y,即圖2(a)~圖2(d)。圖3所示的結(jié)果表明,在低頻段高階模式聲波數(shù)目相對(duì)較少時(shí),各種空間分布的次級(jí)聲源均能取得較好的控制效果。這是因?yàn)樵谶@些頻率范圍內(nèi)壁面分布次級(jí)聲源對(duì)所有高階模式聲波均能實(shí)現(xiàn)獨(dú)立控制,從而取得了接近最優(yōu)的控制效果。但是,圖3所示的壁面分布次級(jí)聲源對(duì)高階模式控制效果存在明顯的工作頻率上限。例如,對(duì)于圖3中的次級(jí)聲源空間分布(1),其工作頻率上限為56.7 Hz,對(duì)應(yīng)第3階(0,1)模式聲波的截止頻率。這是因?yàn)楫?dāng)頻率低于56.7 Hz時(shí),管道中僅有第1階(0,0)和第2階(1,0)模式聲波??臻g分布(1)為一個(gè)x方向上的線陣,可以獨(dú)立控制不同nx的模式聲波。當(dāng)頻率高于56.7 Hz時(shí)管道中出現(xiàn)了第3階(0,1)模式聲波,而空間分布(1)所示的次級(jí)聲源無(wú)法獨(dú)立控制(0,0)和(0,1)兩個(gè)nx=0的簡(jiǎn)并模式,導(dǎo)致控制效果顯著下降。而對(duì)于空間分布(2),通過(guò)增加了一個(gè)x方向上的線陣,使得其可以獨(dú)立控制ny相關(guān)的兩個(gè)模式,因此可以同時(shí)控制(0,0)和(0,1)模式聲波,相應(yīng)的工作頻率上限為(0,2)模式聲波的截止頻率113.3 Hz。此時(shí)空間分布(2)對(duì)應(yīng)的次級(jí)聲源無(wú)法同時(shí)獨(dú)立控制(0,0)、(0,1)和(0,2)這3個(gè)nx=0的簡(jiǎn)并模式。同理可以解釋空間分布(3)的工作頻率上限為42.5 Hz是因?yàn)闊o(wú)法獨(dú)立控制(0,0)和(1,0)模式,而空間分布(4)的工作頻率上限為85.0 Hz是因?yàn)闊o(wú)法獨(dú)立控制(0,0)、(1,0)和(2,0)三個(gè)簡(jiǎn)并模式。

    圖3 壁面分布次級(jí)聲源的空間分布與對(duì)應(yīng)的噪聲控制效果Fig.3 The spatial distribution and the corresponding control performances of the boundary-located secondary sources

    圖4給出了管道中單層和多層壁面分布次級(jí)聲源的空間分布示意圖和相應(yīng)的控制效果。其中,空間分布(1)和空間分布(2)分別對(duì)應(yīng)前文的X+Y和(X+Y)·Z,即對(duì)應(yīng)了圖2(e)和圖2(f)。通過(guò)將x方向上和y方向上組合起來(lái)的單層壁面分布次級(jí)聲源,對(duì)高階模式聲波的控制效果更好。此時(shí),其工作頻率的上限提升至141.7 Hz,對(duì)應(yīng)第11階(2,2)模式聲波的截止頻率,因?yàn)?2,2)模式聲波在x和y方向上均有3個(gè)簡(jiǎn)并模式而單層壁面分布次級(jí)聲源無(wú)法同時(shí)控制這些簡(jiǎn)并模式。通過(guò)在z方向增加一層次級(jí)聲源,構(gòu)成軸向間距ΔLz=1 m的雙層分布,相比于單層分布的次級(jí)聲源能夠在高于141.7 Hz以上的頻段取得20 dB以上的降噪量提升。該結(jié)果表明,次級(jí)聲源沿著管道軸向分布時(shí)能更好地控制高階模式聲波,在頻率相對(duì)較高和高階模式聲波數(shù)目較多時(shí)對(duì)于控制效果的改進(jìn)更為明顯。

    圖4 壁面分布次級(jí)聲源的空間分布與對(duì)應(yīng)的噪聲控制效果Fig.4 The spatial distribution and the corresponding control performances of the boundary-located secondary sources

    進(jìn)一步研究壁面分布次級(jí)聲源的數(shù)量和間距對(duì)高階模式聲波控制效果的影響。首先討論次級(jí)聲源在管道周向上數(shù)目的影響。圖5中展示了管道周向上6個(gè)、14個(gè)和22個(gè)壁面分布單層次級(jí)聲源以及對(duì)應(yīng)的控制效果。相比于周向上14個(gè)次級(jí)聲源,周向上僅有6個(gè)次級(jí)聲源時(shí)其工作頻率上限降至113.3 Hz,對(duì)應(yīng)著(0,2)模式聲波的截止頻率。這是因?yàn)榇藭r(shí)次級(jí)聲源在y方向上僅有2個(gè)次級(jí)聲源,根據(jù)空間采樣定律只能獨(dú)立控制ny=0和ny=1兩類(lèi)模式聲波,當(dāng)頻率超過(guò)113.3 Hz后(0,2)模式聲波開(kāi)始在管道傳播時(shí),壁面分布次級(jí)聲源將失去控制效果。但是,次級(jí)聲源在管道周向上的數(shù)目增加至22個(gè)并沒(méi)有提升控制效果。根據(jù)采樣定律這是因?yàn)?2個(gè)次級(jí)聲源屬于過(guò)采樣,不能改善控制效果。

    圖5 壁面分布次級(jí)聲源的周向上次級(jí)聲源數(shù)量對(duì)控制效果的影響Fig.5 The influence of the number of the single-layer boundary-located secondary sources in the circumferential direction on the control performance

    接下來(lái)研究次級(jí)聲源在管道軸向上的數(shù)目和間距的影響。將壁面分布次級(jí)聲源在周向上固定分布為14個(gè),即按照?qǐng)D5(a)中空間分布(1),層間距固定為ΔLz=1 m,軸向上不同層數(shù)的次級(jí)聲源對(duì)高階模式的控制效果如圖6所示。圖6中14×2表示在管道周向上有14個(gè)次級(jí)聲源,在軸向上共有2層,其余表述以此類(lèi)推。圖6的結(jié)果表明,通過(guò)增加壁面分布次級(jí)聲源在管道軸向上的層數(shù),擴(kuò)大次級(jí)聲源在管道軸向上的分布范圍,能進(jìn)一步提高控制效果,尤其在相對(duì)高頻段,高階模式聲波數(shù)目較多時(shí)控制效果提升更為明顯。當(dāng)次級(jí)聲源的層數(shù)超過(guò)6層時(shí),即軸向尺寸Lz超過(guò)5 m,壁面分布次級(jí)聲源對(duì)高階模式聲波的控制效果達(dá)到接近最優(yōu),實(shí)現(xiàn)近60 dB的降噪量。

    圖6 軸向上不同層數(shù)的壁面分布次級(jí)聲源對(duì)高階模式控制效果Fig.6 The control performances achieved by different layers of boundary-located secondary sources

    最后討論壁面分布次級(jí)聲源的軸向間距對(duì)控制效果的影響。本文提出的理論模型表明次級(jí)聲源分布的軸向尺寸較大時(shí)有利于對(duì)控制高階高階模式,因此通過(guò)使用雙層分布次級(jí)聲源并增加層間距也能起到增大軸向尺寸的效果。壁面分布的14×2雙層次級(jí)聲源不同軸向間距對(duì)高階模式聲波的控制效果如圖7所示。圖7的結(jié)果表明,采樣壁面分布雙層次級(jí)聲源并適當(dāng)增大軸向間距同樣能夠獲得較好的控制效果,尤其在相對(duì)中高頻段增大軸向間距可以顯著提高控制效果。但應(yīng)該注意到在軸向間距很大時(shí)可能在部分頻段出現(xiàn)降噪量下降的現(xiàn)象,例如對(duì)于軸向間距ΔLz=5 m相對(duì)于ΔLz=2 m在140 Hz左右對(duì)高階模式的控制效果反而出現(xiàn)了惡化。這主要是軸向間距過(guò)大導(dǎo)致在軸向上不滿足采樣定律,使得控制效果惡化。盡管如此,總體來(lái)看采用壁面分布雙層次級(jí)聲源并適當(dāng)增大軸向?qū)娱g距,能夠以相對(duì)較少的次級(jí)聲源實(shí)現(xiàn)對(duì)高階模式聲波的控制。

    圖7 壁面分布雙層次級(jí)聲源的不同軸向間距對(duì)高階模式聲波的控制效果Fig.7 The control performances achieved by different axial spacings of the boundary-located double-layer secondary sources

    本節(jié)的仿真實(shí)驗(yàn)?zāi)M了理論模型中所述幾種典型的壁面分布次級(jí)聲源的控制性能。仿真結(jié)果顯示,壁面分布次級(jí)聲源控制高階模式聲波的工作頻率范圍明顯受限于次級(jí)聲源的空間分布。其中,僅沿x方向分布的次級(jí)聲源,其工作頻率上限不超過(guò)(0,2)模式聲波的截止頻率;僅沿y方向分布的次級(jí)聲源,其工作頻率不超過(guò)(2,0)模式聲波的截止頻率;同時(shí)沿x和y方向分布的次級(jí)聲源的工作頻率上限則可以提高至(2,2)模式聲波的截止頻率。根據(jù)理論模型推理可得,這是因?yàn)榇渭?jí)聲源沿x和y方向的空間分布對(duì)高階模式聲波的貢獻(xiàn)只能控制mx和my相關(guān)的高階模式,且由于對(duì)于矩形管道在x和y方向只有兩面墻,使得在x和y方向最多控制兩個(gè)不同的高階模式。此時(shí),在更高頻率段次級(jí)聲源對(duì)高階模式的控制性能主要取決于其軸向分布范圍。當(dāng)次級(jí)聲源在管道軸向上范圍較大時(shí),能顯著提高對(duì)高階模式聲波的控制效果。

    在理論模型中通過(guò)將次級(jí)聲源的分布離散化,根據(jù)采樣定律也能得到次級(jí)聲源在各方向上的數(shù)量和間距對(duì)控制效果的影響。壁面分布次級(jí)聲源在管道周向上的數(shù)目及間距需要滿足采樣定律,但更多的次級(jí)聲源會(huì)導(dǎo)致過(guò)采樣,對(duì)提升控制效果沒(méi)有幫助。在管道軸向上,盡管壁面分布次級(jí)聲源在軸向間距過(guò)大時(shí)可能造成控制性能下降,但由于增大了軸向分布范圍因此整體上依然能提高控制性能。

    3 結(jié)論

    針對(duì)管道有源噪聲控制中的壁面分布次級(jí)聲源空間分布優(yōu)化設(shè)計(jì)問(wèn)題,本文提出了壁面分布次級(jí)聲源獨(dú)立控制高階模式聲波的理論模型,顯式推導(dǎo)了次級(jí)聲源的空間分布對(duì)于控制高階模式聲波的貢獻(xiàn),分別指出次級(jí)聲源在管道周向和軸向空間分布的優(yōu)化準(zhǔn)則;然后通過(guò)將次級(jí)聲源空間分布離散化介紹了點(diǎn)聲源控制策略,并根據(jù)采樣定律分析了次級(jí)聲源的數(shù)目和間距對(duì)于控制效果的影響。通過(guò)對(duì)比多種次級(jí)聲源空間分布以及對(duì)應(yīng)的控制性能,其結(jié)果表明為了優(yōu)化壁面分布次級(jí)聲源的空間分布,需要在管道周向所有壁面布放足夠數(shù)量的次級(jí)聲源以滿足采樣定律,同時(shí)在軸向上增大間距來(lái)獲得較大的分布范圍。

    猜你喜歡
    聲源聲波高階
    虛擬聲源定位的等效源近場(chǎng)聲全息算法
    有限圖上高階Yamabe型方程的非平凡解
    高階各向異性Cahn-Hilliard-Navier-Stokes系統(tǒng)的弱解
    滾動(dòng)軸承壽命高階計(jì)算與應(yīng)用
    哈爾濱軸承(2020年1期)2020-11-03 09:16:02
    基于GCC-nearest時(shí)延估計(jì)的室內(nèi)聲源定位
    電子制作(2019年23期)2019-02-23 13:21:12
    愛(ài)的聲波 將愛(ài)留在她身邊
    聲波殺手
    自適應(yīng)BPSK在井下鉆柱聲波傳輸中的應(yīng)用
    運(yùn)用內(nèi)積相關(guān)性結(jié)合迭代相減識(shí)別兩點(diǎn)聲源
    “聲波驅(qū)蚊”靠譜嗎
    桃花免费在线播放| 亚洲一区二区三区欧美精品| 免费看十八禁软件| 亚洲国产av新网站| 国产亚洲欧美精品永久| 桃花免费在线播放| 亚洲精品av麻豆狂野| 国产免费一区二区三区四区乱码| 老司机靠b影院| av欧美777| 精品久久久久久久毛片微露脸 | 中文字幕另类日韩欧美亚洲嫩草| 国产精品成人在线| 国产精品免费大片| 超碰成人久久| 国产成人系列免费观看| 男女午夜视频在线观看| 日韩欧美一区视频在线观看| 90打野战视频偷拍视频| 亚洲精品久久久久久婷婷小说| 黄色一级大片看看| 人体艺术视频欧美日本| 亚洲人成电影观看| 精品一区二区三区四区五区乱码 | 中文乱码字字幕精品一区二区三区| 少妇 在线观看| 中文字幕制服av| 搡老岳熟女国产| 中文字幕人妻丝袜一区二区| 亚洲av日韩在线播放| 国产精品成人在线| 亚洲欧洲日产国产| 我要看黄色一级片免费的| 国产精品国产三级国产专区5o| 乱人伦中国视频| 蜜桃在线观看..| 成人亚洲欧美一区二区av| 亚洲精品日韩在线中文字幕| 久久这里只有精品19| 男女国产视频网站| 欧美人与性动交α欧美软件| 国产欧美日韩一区二区三区在线| 日本欧美视频一区| 亚洲中文av在线| 男女边摸边吃奶| 中文欧美无线码| 午夜激情久久久久久久| 大码成人一级视频| av国产久精品久网站免费入址| 亚洲精品第二区| 极品少妇高潮喷水抽搐| 在线天堂中文资源库| 久久精品国产亚洲av高清一级| 国产精品一国产av| 亚洲伊人久久精品综合| 啦啦啦中文免费视频观看日本| 老司机深夜福利视频在线观看 | 久久久久视频综合| 啦啦啦中文免费视频观看日本| 搡老岳熟女国产| 91字幕亚洲| 婷婷色av中文字幕| 老司机深夜福利视频在线观看 | 免费在线观看影片大全网站 | 久久热在线av| 久久久久视频综合| 日本色播在线视频| www.av在线官网国产| 久久亚洲精品不卡| tube8黄色片| 亚洲精品中文字幕在线视频| 亚洲国产精品成人久久小说| 国产在视频线精品| 亚洲成人手机| 大片免费播放器 马上看| 丝袜美足系列| 欧美亚洲日本最大视频资源| 99国产精品免费福利视频| 99国产精品一区二区蜜桃av | 国产成人系列免费观看| 成年美女黄网站色视频大全免费| 亚洲欧洲日产国产| 一区二区三区激情视频| av有码第一页| 啦啦啦在线观看免费高清www| 中文字幕高清在线视频| 中文字幕制服av| 亚洲国产日韩一区二区| 又大又爽又粗| 菩萨蛮人人尽说江南好唐韦庄| 国产精品 国内视频| 99久久精品国产亚洲精品| 国产一区二区 视频在线| 超碰成人久久| 天天躁日日躁夜夜躁夜夜| 亚洲av片天天在线观看| 婷婷成人精品国产| 久久鲁丝午夜福利片| 日韩av不卡免费在线播放| 波野结衣二区三区在线| 日本a在线网址| 天堂中文最新版在线下载| 欧美精品一区二区免费开放| 美女扒开内裤让男人捅视频| 男人爽女人下面视频在线观看| www.自偷自拍.com| 悠悠久久av| 亚洲成人免费av在线播放| 少妇粗大呻吟视频| 亚洲少妇的诱惑av| 美女高潮到喷水免费观看| 性色av乱码一区二区三区2| 国产极品粉嫩免费观看在线| 欧美久久黑人一区二区| 美女国产高潮福利片在线看| 国产成人系列免费观看| 亚洲精品国产av蜜桃| a级片在线免费高清观看视频| 人体艺术视频欧美日本| 久久中文字幕一级| 欧美日韩福利视频一区二区| 黄色视频在线播放观看不卡| 国产精品 国内视频| 国产成人精品无人区| 观看av在线不卡| 热re99久久国产66热| 久久人人97超碰香蕉20202| 国产主播在线观看一区二区 | 两性夫妻黄色片| 日韩制服骚丝袜av| 国产国语露脸激情在线看| 最新在线观看一区二区三区 | 丝袜美足系列| 国产爽快片一区二区三区| 精品人妻熟女毛片av久久网站| 欧美日韩成人在线一区二区| 丝袜美足系列| 精品国产一区二区三区四区第35| 亚洲美女黄色视频免费看| 国产高清videossex| 人人澡人人妻人| 18禁观看日本| 丝袜美腿诱惑在线| 国产人伦9x9x在线观看| 亚洲欧美精品综合一区二区三区| 1024香蕉在线观看| 欧美日韩成人在线一区二区| 精品高清国产在线一区| 天天躁狠狠躁夜夜躁狠狠躁| 免费看十八禁软件| 日本av手机在线免费观看| 麻豆乱淫一区二区| h视频一区二区三区| 啦啦啦啦在线视频资源| 少妇精品久久久久久久| 国产男女超爽视频在线观看| 久久国产精品人妻蜜桃| 久久久精品区二区三区| 国产有黄有色有爽视频| 中文精品一卡2卡3卡4更新| 狂野欧美激情性bbbbbb| av在线app专区| 亚洲精品美女久久av网站| 一二三四社区在线视频社区8| 免费高清在线观看日韩| 男女之事视频高清在线观看 | 国产亚洲av片在线观看秒播厂| 在线观看人妻少妇| 欧美av亚洲av综合av国产av| 极品少妇高潮喷水抽搐| 国产av精品麻豆| 极品人妻少妇av视频| 成年人午夜在线观看视频| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕人妻丝袜制服| 国产成人91sexporn| 国产av国产精品国产| 亚洲精品在线美女| 美女中出高潮动态图| 男女之事视频高清在线观看 | 一边摸一边做爽爽视频免费| 精品亚洲成国产av| 一级,二级,三级黄色视频| 免费av中文字幕在线| 中文字幕亚洲精品专区| 19禁男女啪啪无遮挡网站| 欧美日韩福利视频一区二区| 午夜91福利影院| 国产精品免费视频内射| 成人免费观看视频高清| 又粗又硬又长又爽又黄的视频| 麻豆乱淫一区二区| 女性生殖器流出的白浆| av一本久久久久| 精品少妇一区二区三区视频日本电影| 老汉色av国产亚洲站长工具| 可以免费在线观看a视频的电影网站| 无限看片的www在线观看| 乱人伦中国视频| 男女下面插进去视频免费观看| 2018国产大陆天天弄谢| 建设人人有责人人尽责人人享有的| 午夜福利乱码中文字幕| 色婷婷av一区二区三区视频| 丝袜脚勾引网站| 啦啦啦 在线观看视频| 操出白浆在线播放| 咕卡用的链子| 久久精品aⅴ一区二区三区四区| 高清黄色对白视频在线免费看| 校园人妻丝袜中文字幕| 三上悠亚av全集在线观看| 韩国精品一区二区三区| 美女脱内裤让男人舔精品视频| 69精品国产乱码久久久| 黄色毛片三级朝国网站| 18禁观看日本| 中文字幕另类日韩欧美亚洲嫩草| 免费久久久久久久精品成人欧美视频| 免费看av在线观看网站| 久久久亚洲精品成人影院| 大片电影免费在线观看免费| 色视频在线一区二区三区| 免费看av在线观看网站| 国产精品麻豆人妻色哟哟久久| 中文字幕最新亚洲高清| 精品久久久久久电影网| 十八禁人妻一区二区| 蜜桃国产av成人99| 国产成人精品久久二区二区免费| 精品福利观看| 精品亚洲乱码少妇综合久久| 日韩熟女老妇一区二区性免费视频| 永久免费av网站大全| 最黄视频免费看| 欧美中文综合在线视频| 一二三四在线观看免费中文在| 精品亚洲乱码少妇综合久久| 一区福利在线观看| 视频区欧美日本亚洲| 亚洲美女黄色视频免费看| 久久久久国产一级毛片高清牌| 中文乱码字字幕精品一区二区三区| 中文字幕制服av| 午夜福利,免费看| 香蕉丝袜av| 亚洲欧美精品自产自拍| 在线观看免费午夜福利视频| 精品一区在线观看国产| 精品卡一卡二卡四卡免费| 国产精品久久久久成人av| 久久久精品区二区三区| 国产精品香港三级国产av潘金莲 | 国产成人精品久久久久久| 欧美亚洲 丝袜 人妻 在线| 久久99精品国语久久久| 自拍欧美九色日韩亚洲蝌蚪91| 在现免费观看毛片| 成年av动漫网址| 一二三四社区在线视频社区8| 欧美激情 高清一区二区三区| 国产成人精品久久二区二区91| 国产精品一国产av| 操美女的视频在线观看| 亚洲国产欧美网| 天天躁夜夜躁狠狠躁躁| 国产成人精品在线电影| 国产欧美日韩综合在线一区二区| 欧美黄色淫秽网站| 99久久人妻综合| 十八禁高潮呻吟视频| 岛国毛片在线播放| 久久九九热精品免费| 亚洲精品久久午夜乱码| 曰老女人黄片| 亚洲欧美色中文字幕在线| 另类亚洲欧美激情| 丝袜美腿诱惑在线| 菩萨蛮人人尽说江南好唐韦庄| 91老司机精品| 欧美日韩视频高清一区二区三区二| 亚洲精品一卡2卡三卡4卡5卡 | 伊人亚洲综合成人网| 欧美日韩福利视频一区二区| 男女国产视频网站| 在现免费观看毛片| 国产精品国产三级国产专区5o| av线在线观看网站| 99国产精品免费福利视频| 男女下面插进去视频免费观看| 国产成人91sexporn| 精品少妇久久久久久888优播| 美女国产高潮福利片在线看| 国产亚洲精品第一综合不卡| 视频区欧美日本亚洲| 你懂的网址亚洲精品在线观看| 在线av久久热| 丰满迷人的少妇在线观看| 国产亚洲一区二区精品| 久久天堂一区二区三区四区| 亚洲七黄色美女视频| 操美女的视频在线观看| 91精品伊人久久大香线蕉| 成人国产一区最新在线观看 | 午夜免费男女啪啪视频观看| 亚洲一码二码三码区别大吗| 亚洲国产欧美一区二区综合| 欧美日韩一级在线毛片| 黑人欧美特级aaaaaa片| 手机成人av网站| 伊人久久大香线蕉亚洲五| 九草在线视频观看| 亚洲美女黄色视频免费看| 久久久久久久精品精品| 韩国高清视频一区二区三区| www.999成人在线观看| 伊人亚洲综合成人网| 亚洲欧美日韩高清在线视频 | 国产精品亚洲av一区麻豆| a级毛片黄视频| 国产精品香港三级国产av潘金莲 | 18在线观看网站| 色婷婷av一区二区三区视频| 一级毛片我不卡| 亚洲 国产 在线| 亚洲av美国av| e午夜精品久久久久久久| 亚洲三区欧美一区| 国产亚洲精品久久久久5区| 国产精品免费大片| 日韩 亚洲 欧美在线| 欧美黄色片欧美黄色片| 亚洲国产欧美一区二区综合| www.熟女人妻精品国产| 又紧又爽又黄一区二区| 日本91视频免费播放| 美女午夜性视频免费| 交换朋友夫妻互换小说| 丁香六月欧美| videos熟女内射| 日本a在线网址| 在线av久久热| 免费在线观看日本一区| 久久久久久人人人人人| 在线天堂中文资源库| 亚洲一卡2卡3卡4卡5卡精品中文| 男女高潮啪啪啪动态图| 九色亚洲精品在线播放| 国产主播在线观看一区二区 | 久久久久国产一级毛片高清牌| 久久精品久久精品一区二区三区| 亚洲国产成人一精品久久久| 在线天堂中文资源库| 777米奇影视久久| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲国产欧美在线一区| 国产1区2区3区精品| 国产亚洲av片在线观看秒播厂| 叶爱在线成人免费视频播放| 波多野结衣一区麻豆| 亚洲欧洲精品一区二区精品久久久| 一级黄片播放器| 老熟女久久久| 久久精品国产亚洲av涩爱| 日本av免费视频播放| 日日爽夜夜爽网站| 在线观看www视频免费| 成人免费观看视频高清| 亚洲欧美成人综合另类久久久| 亚洲色图综合在线观看| 国产亚洲av片在线观看秒播厂| www.999成人在线观看| 最黄视频免费看| 精品久久久精品久久久| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品美女久久久久99蜜臀 | 亚洲中文日韩欧美视频| 国产爽快片一区二区三区| 亚洲九九香蕉| 国产一区二区 视频在线| 久久久久久人人人人人| 国产一区二区激情短视频 | 亚洲专区国产一区二区| 国产在线观看jvid| 女性生殖器流出的白浆| 亚洲av电影在线观看一区二区三区| 中文字幕人妻熟女乱码| 黑人欧美特级aaaaaa片| 天天躁夜夜躁狠狠躁躁| 成人亚洲精品一区在线观看| 国产成人91sexporn| 妹子高潮喷水视频| 色94色欧美一区二区| 少妇人妻 视频| av在线播放精品| 三上悠亚av全集在线观看| 少妇裸体淫交视频免费看高清 | 久久久久国产一级毛片高清牌| 中文字幕精品免费在线观看视频| 丝袜脚勾引网站| av欧美777| 超色免费av| av网站在线播放免费| 亚洲精品国产av成人精品| 日韩 欧美 亚洲 中文字幕| 欧美日韩视频精品一区| 夜夜骑夜夜射夜夜干| 91麻豆精品激情在线观看国产 | 色精品久久人妻99蜜桃| 国精品久久久久久国模美| 精品少妇一区二区三区视频日本电影| 成人国产一区最新在线观看 | 日韩 欧美 亚洲 中文字幕| 永久免费av网站大全| 久久久久精品国产欧美久久久 | 亚洲一区二区三区欧美精品| 精品欧美一区二区三区在线| 亚洲欧美中文字幕日韩二区| 午夜福利乱码中文字幕| 日韩大码丰满熟妇| 国产一区二区 视频在线| 蜜桃国产av成人99| 最近中文字幕2019免费版| 成人18禁高潮啪啪吃奶动态图| av在线老鸭窝| 汤姆久久久久久久影院中文字幕| 国产99久久九九免费精品| 亚洲成人免费av在线播放| 男人操女人黄网站| 日本av手机在线免费观看| 午夜福利一区二区在线看| 亚洲人成电影免费在线| 一本久久精品| 欧美黑人精品巨大| 麻豆国产av国片精品| 男人爽女人下面视频在线观看| 国产91精品成人一区二区三区 | a级毛片在线看网站| www.999成人在线观看| 超色免费av| 啦啦啦中文免费视频观看日本| 亚洲图色成人| 最近最新中文字幕大全免费视频 | 人人澡人人妻人| 欧美 日韩 精品 国产| 一级毛片黄色毛片免费观看视频| 婷婷色综合大香蕉| a 毛片基地| 99re6热这里在线精品视频| 久久精品久久久久久噜噜老黄| 欧美性长视频在线观看| 一本一本久久a久久精品综合妖精| 777米奇影视久久| 一二三四在线观看免费中文在| 黄色一级大片看看| av网站在线播放免费| 国产无遮挡羞羞视频在线观看| videosex国产| 欧美日韩视频高清一区二区三区二| 亚洲精品日本国产第一区| 一本一本久久a久久精品综合妖精| 欧美日韩福利视频一区二区| 婷婷色综合www| 一二三四社区在线视频社区8| 中文精品一卡2卡3卡4更新| 一区二区三区激情视频| 亚洲国产看品久久| 免费在线观看黄色视频的| 王馨瑶露胸无遮挡在线观看| 男人添女人高潮全过程视频| 婷婷色综合大香蕉| 国产精品一二三区在线看| 亚洲九九香蕉| 精品高清国产在线一区| 精品亚洲成a人片在线观看| 女警被强在线播放| 9热在线视频观看99| 最近中文字幕2019免费版| 男女边吃奶边做爰视频| 日韩制服丝袜自拍偷拍| 亚洲av欧美aⅴ国产| √禁漫天堂资源中文www| videos熟女内射| 久久精品亚洲av国产电影网| 99热全是精品| 亚洲九九香蕉| 国产免费视频播放在线视频| av视频免费观看在线观看| 国产精品三级大全| 悠悠久久av| 两性夫妻黄色片| 老司机影院成人| 久久国产精品大桥未久av| 久久综合国产亚洲精品| 精品久久久精品久久久| 搡老岳熟女国产| 波多野结衣一区麻豆| 午夜福利影视在线免费观看| 视频区欧美日本亚洲| 可以免费在线观看a视频的电影网站| 丝瓜视频免费看黄片| 老鸭窝网址在线观看| 色婷婷久久久亚洲欧美| 91精品三级在线观看| 中文精品一卡2卡3卡4更新| 啦啦啦啦在线视频资源| 亚洲色图 男人天堂 中文字幕| 又粗又硬又长又爽又黄的视频| 免费少妇av软件| 少妇人妻 视频| 成人国产一区最新在线观看 | 久久久久久久精品精品| av在线app专区| 国产色视频综合| 又紧又爽又黄一区二区| 19禁男女啪啪无遮挡网站| av一本久久久久| 波多野结衣av一区二区av| 国产97色在线日韩免费| 免费在线观看影片大全网站 | 777米奇影视久久| 欧美精品高潮呻吟av久久| 亚洲av国产av综合av卡| 欧美日韩av久久| 男女之事视频高清在线观看 | 久久国产精品影院| 国产一区二区激情短视频 | 精品熟女少妇八av免费久了| 免费日韩欧美在线观看| 欧美xxⅹ黑人| av国产精品久久久久影院| 十八禁高潮呻吟视频| 最黄视频免费看| 国产精品久久久久久精品电影小说| 国产av精品麻豆| 丰满迷人的少妇在线观看| 亚洲,欧美精品.| 亚洲精品乱久久久久久| 亚洲伊人色综图| 又大又爽又粗| 久热爱精品视频在线9| 丝袜美腿诱惑在线| 夜夜骑夜夜射夜夜干| 亚洲中文av在线| 一区二区三区乱码不卡18| 久久鲁丝午夜福利片| 少妇猛男粗大的猛烈进出视频| 欧美在线黄色| 七月丁香在线播放| 黑人巨大精品欧美一区二区蜜桃| 99精品久久久久人妻精品| 少妇 在线观看| 一级,二级,三级黄色视频| 一本久久精品| 狠狠婷婷综合久久久久久88av| 欧美少妇被猛烈插入视频| 国产欧美日韩一区二区三 | 亚洲中文字幕日韩| 性色av乱码一区二区三区2| 亚洲精品国产色婷婷电影| 人妻 亚洲 视频| 两人在一起打扑克的视频| 麻豆乱淫一区二区| 十八禁网站网址无遮挡| 国产男女内射视频| 侵犯人妻中文字幕一二三四区| 男的添女的下面高潮视频| 夜夜骑夜夜射夜夜干| 免费观看av网站的网址| 大片电影免费在线观看免费| 男人添女人高潮全过程视频| 桃花免费在线播放| 国产日韩欧美在线精品| 黄色片一级片一级黄色片| 国产亚洲av高清不卡| 天天操日日干夜夜撸| 亚洲欧美一区二区三区久久| 久久精品aⅴ一区二区三区四区| 99热网站在线观看| 亚洲中文字幕日韩| 操出白浆在线播放| av在线app专区| 国产国语露脸激情在线看| 黑人巨大精品欧美一区二区蜜桃| 欧美成狂野欧美在线观看| 免费看av在线观看网站| 日韩熟女老妇一区二区性免费视频| 精品少妇一区二区三区视频日本电影| 欧美成狂野欧美在线观看| 国产精品香港三级国产av潘金莲 | 国产精品一二三区在线看| 欧美精品人与动牲交sv欧美| 久久久久国产精品人妻一区二区| 成人国产av品久久久| 男女边摸边吃奶| 亚洲欧美日韩另类电影网站| 18禁黄网站禁片午夜丰满| 国产午夜精品一二区理论片| 久久久久久亚洲精品国产蜜桃av| 18禁国产床啪视频网站| 五月天丁香电影| 国产黄频视频在线观看| 日韩精品免费视频一区二区三区| 久久精品国产综合久久久| 肉色欧美久久久久久久蜜桃| 纯流量卡能插随身wifi吗| 精品一区在线观看国产| 国产亚洲欧美在线一区二区| 国产精品国产三级国产专区5o| 日韩大码丰满熟妇| 亚洲av在线观看美女高潮| 精品第一国产精品| 欧美黑人欧美精品刺激| 国产片特级美女逼逼视频| 99热网站在线观看|