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

    聲振耦合聲場(chǎng)分析與聲輻射結(jié)構(gòu)優(yōu)化

    2021-12-13 02:26:04王旭東朱擁勇王德石
    關(guān)鍵詞:優(yōu)化結(jié)構(gòu)

    王旭東,朱擁勇,王德石

    (1.海軍工程大學(xué) 兵器工程學(xué)院,武漢 430033;2.中國(guó)人民解放軍92064部隊(duì),廣東 東莞 523900)

    同時(shí)具有內(nèi)部聲場(chǎng)和外部自由聲場(chǎng)的聲振耦合系統(tǒng)廣泛存在于汽車、飛機(jī)和潛艇等運(yùn)載工具和武器平臺(tái)之中。結(jié)構(gòu)與聲介質(zhì)因振動(dòng)相互作用而產(chǎn)生的中低頻噪聲會(huì)對(duì)人員及設(shè)備產(chǎn)生負(fù)面影響,輻射噪聲影響人員舒適性和造成環(huán)境污染,在軍事上,聲輻射過(guò)大會(huì)影響武器系統(tǒng)的隱身性和自身的聲吶探測(cè)距離[1]。因此在準(zhǔn)確計(jì)算聲輻射的基礎(chǔ)上,通過(guò)適當(dāng)設(shè)計(jì)材料的結(jié)構(gòu)分布以使耦合系統(tǒng)在所關(guān)心的激勵(lì)頻段下具有較低的聲輻射值的研究十分必要。

    作為新興的結(jié)構(gòu)動(dòng)力優(yōu)化方法,不同于尺寸優(yōu)化[2]和形狀優(yōu)化[3]僅能改變已知設(shè)計(jì)部件的尺寸,拓?fù)鋬?yōu)化的目標(biāo)是尋求最優(yōu)的結(jié)構(gòu)材料分布以滿足性能和結(jié)構(gòu)工藝等多方面要求[4]。國(guó)、內(nèi)外諸多研究者在減振降噪問(wèn)題上采用拓?fù)鋬?yōu)化法開(kāi)展了研究。在內(nèi)聲場(chǎng)優(yōu)化問(wèn)題中,Akl等[5]采用有限元法建立內(nèi)聲場(chǎng)聲輻射模型,采用MMA法優(yōu)化柔性板厚度,降低了聲腔內(nèi)的聲輻射,數(shù)值計(jì)算結(jié)果與實(shí)驗(yàn)吻合較好;Dühring等[6]結(jié)合混合有限元法和MMA優(yōu)化法實(shí)現(xiàn)在無(wú)法明確耦合邊界的情況下解決聲-結(jié)構(gòu)耦合優(yōu)化問(wèn)題;隨后Kook[7]用Comsol軟件實(shí)現(xiàn)將混合有限元法擴(kuò)展至BESO優(yōu)化方法中,優(yōu)化結(jié)果與采用傳統(tǒng)有限元法的文獻(xiàn)[8]結(jié)果一致。在外聲場(chǎng)優(yōu)化問(wèn)題的研究中,陳爐云等[9]結(jié)合ESO優(yōu)化法和單向耦合的邊界元法研究了簡(jiǎn)支板在低頻率下外場(chǎng)聲輻射優(yōu)化問(wèn)題;張軍[10]推導(dǎo)了空氣為聲介質(zhì)的有限元-邊界元耦合方法的聲壓級(jí)和輻射聲功率靈敏度公式,并以板厚為設(shè)計(jì)變量對(duì)薄板聲輻射進(jìn)行優(yōu)化;針對(duì)兩相材料板結(jié)構(gòu)的聲輻射拓?fù)鋬?yōu)化問(wèn)題,文獻(xiàn)[11-12]通過(guò)優(yōu)化布局兩種彈性材料的方式降低聲輻射;此后,黃其柏在上述基礎(chǔ)上研究了板結(jié)構(gòu)聲輻射計(jì)算方法和輻射特性,提出約束阻尼復(fù)合板的拓?fù)鋬?yōu)化方法;Zhang等[13]采用狀態(tài)空間中的復(fù)模態(tài)疊加法求解非比例全局阻尼的復(fù)合板聲響應(yīng),伴隨變量法進(jìn)行聲壓靈敏度分析,通過(guò)合理分配阻尼層材料達(dá)到最小化指定點(diǎn)聲壓的目標(biāo)。

    目前,聲振耦合方法的研究大多針對(duì)單一的有限內(nèi)聲場(chǎng)[5-8]或無(wú)限大/半無(wú)限大自由聲場(chǎng)問(wèn)題[9-13],且在處理自由聲場(chǎng)靈敏度分析時(shí)往往忽略了聲介質(zhì)對(duì)結(jié)構(gòu)的反作用,無(wú)法適用于內(nèi)聲場(chǎng)和外部自由聲場(chǎng)均存在的聲輻射問(wèn)題和強(qiáng)耦合問(wèn)題?;谏鲜鰡?wèn)題,本文結(jié)合聲振耦合有限元方程和邊界元法,推導(dǎo)了內(nèi)外聲場(chǎng)聲壓的耦合計(jì)算公式,并提出一種新的強(qiáng)耦合聲功率的靈敏度求解方法,通過(guò)不同設(shè)計(jì)條件下的聲輻射優(yōu)化設(shè)計(jì)結(jié)果驗(yàn)證了聲輻射預(yù)報(bào)的準(zhǔn)確性和靈敏度方法在優(yōu)化過(guò)程中的適用性。

    1 聲輻射問(wèn)題

    1.1 聲學(xué)基本方程

    聲振耦合示意圖如圖1所示,假定流體是可壓縮的,無(wú)黏且無(wú)旋,忽略自由表面重力效應(yīng)的均勻理想聲學(xué)介質(zhì),聲壓滿足如下的波動(dòng)方程:

    (1)

    ?2pf+(ω/c0)2pf=0

    (2)

    在流固交界面上,應(yīng)滿足如下的運(yùn)動(dòng)學(xué)和力學(xué)條件:

    (3)

    usns+ufnf=0

    (4)

    σs|ns=-pf

    (5)

    式中:ρf為流體密度,nf為耦合表面從聲場(chǎng)指向結(jié)構(gòu)的法向量,ns為耦合表面從結(jié)構(gòu)指向聲場(chǎng)的法向量,σs為結(jié)構(gòu)應(yīng)力,us為結(jié)構(gòu)位移向量。

    圖1 聲振耦合示意

    1.2 內(nèi)聲場(chǎng)聲振耦合的有限元法

    結(jié)構(gòu)域的平衡狀態(tài)方程為

    (6)

    式中:ρs為結(jié)構(gòu)質(zhì)量密度,fs為結(jié)構(gòu)所受外激勵(lì)。

    聯(lián)立式(2)、(3)~式(6),并采用加權(quán)余量的伽遼金法可得到聲振耦合的有限元方程[14]:

    (7)

    式中:Ks、Kf分別為結(jié)構(gòu)域和流體域的系統(tǒng)剛度矩陣,Ms、Mf分別為結(jié)構(gòu)域和流體域系統(tǒng)質(zhì)量矩陣,C為系統(tǒng)耦合矩陣。

    為方便后文推導(dǎo)計(jì)算,式(7)可簡(jiǎn)寫為

    (Kfe-ω2Mfe)Xfe=Ffe

    (8)

    1.3 自由場(chǎng)聲振耦合的邊界元法

    穩(wěn)態(tài)諧激勵(lì)下,為求解結(jié)構(gòu)表面上受到介質(zhì)的壓力,需要建立介質(zhì)中任一點(diǎn)的聲壓積分方程。對(duì)聲壓和格林函數(shù)應(yīng)用格林第二公式,由Sommerfeld輻射定律,外域自由聲場(chǎng)通過(guò)Helmholtz積分方程求解:

    (9)

    式中:px為源點(diǎn)聲壓;py為場(chǎng)點(diǎn)聲壓;c為空間角系數(shù),光滑面上取0.5,聲場(chǎng)內(nèi)取1;G(x,y)為自由空間格林函數(shù);Γbe為自由場(chǎng)聲輻射邊界。

    以二維空間為例,對(duì)于二維聲場(chǎng)問(wèn)題,格林函數(shù)G(x,y)和其對(duì)場(chǎng)點(diǎn)處的法向?qū)?shù)的基本階函數(shù)如下:

    (10)

    (11)

    聯(lián)立式(3),式(9)~式(11),并通過(guò)Kronecker函數(shù)合并聲壓項(xiàng)得到邊界積分方程:

    Hpbe-Gvbe=0

    (12)

    式中:H、G分別為采用邊界積分方程得到的頻率相關(guān)的矩陣,pbe、vbe分別為邊界上的聲壓和法向速度向量。

    1.4 聲振耦合方程及輻射聲功率

    當(dāng)聲場(chǎng)由有限大內(nèi)聲場(chǎng)和無(wú)限大自由聲場(chǎng)組成時(shí),式(8)的平衡方程須引入自由場(chǎng)聲介質(zhì)對(duì)結(jié)構(gòu)表面的反作用力項(xiàng)為

    (Kfe-ω2Mfe)Xfe=Ffe+Lpbe

    (13)

    式中L為結(jié)構(gòu)與自由聲場(chǎng)的邊界壓強(qiáng)與邊界法向作用力的耦合矩陣。

    (14)

    由穩(wěn)態(tài)諧波激勵(lì)下的界面位移與界面法向速度的界面連續(xù)性條件可知:

    (15)

    式中:Nu為結(jié)構(gòu)單元的位移形函數(shù),Γe為自由聲場(chǎng)與結(jié)構(gòu)的單元耦合邊界,m為耦合邊界單元總數(shù)。上式求積分并求逆移項(xiàng)后,將結(jié)構(gòu)位移從邊界單元擴(kuò)充至全局響應(yīng)向量Xfe得到:

    (16)

    式中:Mbe為邊界質(zhì)量矩陣,其表達(dá)式即為式(15)左端積分項(xiàng);聯(lián)立式(12)、(13)和式(16),得到修正后含外部自由聲場(chǎng)的聲振耦合方程

    (17)

    直接求解式(17)即可得到內(nèi)外聲場(chǎng)聲壓值以及結(jié)構(gòu)邊界位移。同式(14),采用形函數(shù)Np對(duì)自由場(chǎng)邊界聲速和聲壓進(jìn)行采用統(tǒng)一的插值離散并代入上式,得到聲功率W的數(shù)值計(jì)算公式為[15]

    (18)

    式中Re為取計(jì)算值的實(shí)部。

    式(18)在聲功率靈敏度求解中不便于解耦變量,為方便后文求解單元靈敏度,需要將聲功率公式修改為用結(jié)構(gòu)位移響應(yīng)表達(dá)的方式。將邊界積分方程式(12)代入上式并由復(fù)共軛矩陣乘法規(guī)則得到采用邊界振速表示的結(jié)構(gòu)聲功率的表達(dá)式:

    (19)

    聲功率靈敏度是求聲功率對(duì)結(jié)構(gòu)單元狀態(tài)參數(shù)xi求微分,為方便計(jì)算,應(yīng)使式(19)盡量簡(jiǎn)單。由實(shí)部操作符得到式(19)的簡(jiǎn)單表達(dá)式:

    (20)

    其中

    將式(16)代入式(20),得到輻射聲功率關(guān)于Xfe的關(guān)系式:

    (21)

    其中

    2 單元靈敏度分析與優(yōu)化算法

    拓?fù)鋬?yōu)化問(wèn)題需要對(duì)目標(biāo)函數(shù)相對(duì)設(shè)計(jì)變量的靈敏度進(jìn)行分析,由輻射聲功率的表達(dá)式可知,在外載荷和輻射邊界不發(fā)生變化的情況下,輻射聲功率僅與結(jié)構(gòu)表面法向速度vbe有關(guān),即是結(jié)構(gòu)響應(yīng)Xfe的函數(shù)。Kang等[16-17]利用伴隨變量法(AVM)推導(dǎo)出廣義結(jié)構(gòu)響應(yīng)函數(shù)g(Xfe)的靈敏度分析方法,該方法明確指出靈敏度依賴于結(jié)構(gòu)位移響應(yīng)的幅值,隨后他們團(tuán)隊(duì)采用這種伴隨變量法計(jì)算僅考慮單向耦合的聲壓靈敏度[13],取得了一定的成果。本文將該方法推廣到強(qiáng)耦合的聲功率靈敏度分析中。

    2.1 聲功率單元靈敏度的AVM方法

    在本文中,首先將AVM方法擴(kuò)展到重介質(zhì)的靈敏度分析中。對(duì)于依賴于位移響應(yīng)Xfe的結(jié)構(gòu)響應(yīng)g(Xfe),其目標(biāo)函數(shù)可以通過(guò)引入兩個(gè)拉格朗日乘子的和改寫為

    (22)

    式中,Kd=Kfe-ω2Mfe。上式中存在與Xfe不統(tǒng)一的變量pbe,聯(lián)立式(12)和式(16)可得

    Lpbe=BXfe

    (23)

    (24)

    將式(24)兩邊對(duì)設(shè)計(jì)變量xi求導(dǎo),可得

    (25)

    令伴隨向量滿足以下方程:

    (26)

    式中μ1、μ2滿足μ1=(μ2)H。式(25)可以化簡(jiǎn)為

    (27)

    由聲振耦合方程(17)求解的Xfe是復(fù)向量,將其分解為實(shí)部與虛部之和的形式代入式(21),得到輻射聲功率的表達(dá)式:

    (28)

    式中( )R、( )I分別為取實(shí)部和虛部。式(28)對(duì)Xfe求導(dǎo)得到:

    (29)

    (30)

    將式(29)、(30)代入式(26)可得到伴隨向量的表達(dá)式:

    (31)

    求解式(31)即可得到拉格朗日乘子μ1,進(jìn)而輻射聲功率靈敏度可由式(27)得出

    (32)

    由上述伴隨變量法的推導(dǎo)過(guò)程可以發(fā)現(xiàn),該方法推導(dǎo)過(guò)程較為繁瑣,中間變量較多,容易出錯(cuò)。因此,本文提出更為簡(jiǎn)潔高效的直接求導(dǎo)法求解輻射聲功率靈敏度。

    2.2 聲功率單元靈敏度的直接求導(dǎo)法

    式(21)中的輻射聲功率W直接對(duì)單元變量xi求導(dǎo),并由實(shí)部操作符得到聲功率的簡(jiǎn)單形式:

    (33)

    式(13)和式(23)分別對(duì)xi求導(dǎo)得到:

    (34)

    由式(34)得到:

    (35)

    將式(35)代入式(33)得到輻射聲功率的單元靈敏度公式:

    (36)

    對(duì)比本文的直接求導(dǎo)法和AVM法可以發(fā)現(xiàn),本文方法的推導(dǎo)過(guò)程較為簡(jiǎn)潔高效,中間變量少。最終的靈敏度公式中的各矩陣都在AVM法的最終公式中出現(xiàn),不需要計(jì)算額外的變量,說(shuō)明靈敏度計(jì)算量也并未增加。

    2.3 目標(biāo)函數(shù)與優(yōu)化算法

    體積約束條件下,以輻射聲功率C最小化為目標(biāo)函數(shù),則結(jié)構(gòu)拓?fù)鋬?yōu)化數(shù)學(xué)模型為:

    (37)

    式中:V*為結(jié)構(gòu)有效體積,xi為單元狀態(tài)變量(取xmin或1)。

    在涉及到特征值問(wèn)題的優(yōu)化中,低密度區(qū)域容易產(chǎn)生局部模態(tài)。在使用BESO優(yōu)化方法時(shí),Huang等[18]對(duì)質(zhì)量和剛度采用相同的罰函數(shù)冪率,該模型對(duì)部分優(yōu)化問(wèn)題是有效的。但是在計(jì)算中發(fā)現(xiàn),采用這種方法后有些優(yōu)化問(wèn)題在失效單元區(qū)域仍然會(huì)出現(xiàn)局部模態(tài),這導(dǎo)致優(yōu)化過(guò)程中結(jié)構(gòu)拓?fù)渫蛔?,迭代不穩(wěn)定難以收斂。Pedersen[19]提出線性化單元?jiǎng)偠鹊姆椒?,該方法可以使單元在極低密度時(shí)仍然能保持較高的剛度質(zhì)量比,可有效避免偽本征模,使優(yōu)化結(jié)果更加合理,罰函數(shù)模型如下:

    (38)

    式中:ρ0為結(jié)構(gòu)初始密度,E0為結(jié)構(gòu)初始彈性模量,p為罰因子。

    結(jié)合BESO方法中單元狀態(tài)變量xi僅可取xmin或1的特點(diǎn),同時(shí)為保證罰函數(shù)模型的連續(xù)性,本文提出以下改進(jìn)的罰函數(shù)模型:

    (39)

    拓?fù)鋬?yōu)化過(guò)程中需要采用基于半徑濾波法和單元靈敏度歷史平均方法等來(lái)避免“孤立網(wǎng)格”和加速收斂。具體的耦合優(yōu)化的計(jì)算流程如下:

    1)初始化結(jié)構(gòu)參數(shù)。定義初始耦合域的物理參數(shù)、定義BESO算法的各初始參數(shù);

    2)生成數(shù)值網(wǎng)格、生成單元矩陣和耦合矩陣并組裝生成有限元聲振耦合方程,生成結(jié)構(gòu)外表面邊界元的各矩陣;

    3)根據(jù)步驟2)計(jì)算各矩陣,求解系統(tǒng)響應(yīng)。根據(jù)式(38)計(jì)算單元靈敏度。進(jìn)行靈敏度過(guò)濾和靈敏度歷史平均[20];

    4)判斷失效單元是否與聲學(xué)單元相鄰,如果相鄰,則將失效單元轉(zhuǎn)變?yōu)槁晫W(xué)單元,并更新結(jié)構(gòu)域和流體域;

    5)根據(jù)更新后的各域,識(shí)別結(jié)構(gòu)域和流體域的耦合單元和耦合邊界;

    6)進(jìn)行收斂性判斷,如不滿足收斂條件,則重復(fù)步驟2)~步驟5),直到滿足收斂條件。

    3 數(shù)值算例

    3.1 聲輻射數(shù)值計(jì)算結(jié)果驗(yàn)證

    為驗(yàn)證本文文建立的聲振耦合方程的正確性,設(shè)置驗(yàn)證算例。無(wú)限大水中放置一個(gè)內(nèi)徑為1.0 m外徑為1.2 m的鋁質(zhì)圓環(huán),結(jié)構(gòu)密度為2 700 kg/m3,彈性模量為69 GPa,泊松比為0.3,內(nèi)聲場(chǎng)和外聲場(chǎng)介質(zhì)為水,其密度為1 000 kg/m3,聲速為1 450 m/s。圓環(huán)內(nèi)側(cè)頂點(diǎn)處施加頻率為10 Hz,幅值為1 N,方向?yàn)閥軸方向的激勵(lì)力。

    采用Matlab編程計(jì)算圓環(huán)內(nèi)外的耦合聲場(chǎng),并與無(wú)反射邊界法的計(jì)算計(jì)算結(jié)果進(jìn)行對(duì)比。在無(wú)反射邊界中,取圓環(huán)外側(cè)半徑5 m的同心圓為外側(cè)計(jì)算域,外邊界設(shè)置為無(wú)反射邊界,材料屬性和其他邊界條件與聲振耦合法相同,兩種方式的聲壓計(jì)算結(jié)果如圖2所示。對(duì)比可知,本文的聲振耦合法與無(wú)反射邊界法計(jì)算的內(nèi)、外場(chǎng)聲壓分布基本一致。

    圖2 圓環(huán)內(nèi)外聲場(chǎng)聲壓分布

    圖3給出圓環(huán)外邊界聲壓實(shí)部的分布曲線,為方便比較,設(shè)置起始點(diǎn)為圓環(huán)左側(cè)中點(diǎn)(與圖3中弧度為0處相對(duì)應(yīng)),沿順時(shí)針?lè)较颉煞N數(shù)值方法計(jì)算的整體匹配度較好,但在谷值處誤差相對(duì)較大,原因在于有限元的無(wú)反射邊界對(duì)正入射波吸收性較好,但對(duì)傾斜入射波有一定的反射作用,因此與實(shí)際值有一定誤差。此外可明顯看出無(wú)反射邊界法的曲線在對(duì)稱節(jié)點(diǎn)上的聲壓值略微不對(duì)稱,而邊界條件和外載荷均是對(duì)稱,其邊界聲壓實(shí)部也應(yīng)該是對(duì)稱的。出現(xiàn)這種現(xiàn)象的原因是劃分有限元網(wǎng)格無(wú)法做到完全對(duì)稱,造成一定的誤差。因此,可以認(rèn)為采用聲振耦合方程可獲得良好的聲壓預(yù)報(bào)結(jié)果。

    圖3 結(jié)構(gòu)表面聲壓實(shí)部分布

    3.2 算例1

    圖4為算例1的設(shè)計(jì)模型,設(shè)計(jì)域?yàn)橐煌獠堪霃?.0 m,內(nèi)部半徑2.0 m的半圓環(huán),外圈厚度0.1 m部分為非設(shè)計(jì)域,初始設(shè)計(jì)域?yàn)槿O(shè)計(jì)域的80%(即圖中半徑2.2 m部分),半圓環(huán)兩側(cè)的底部固定約束,內(nèi)部流體的底部邊界為零聲壓邊界。內(nèi)外聲場(chǎng)的聲介質(zhì)均為水,聲介質(zhì)屬性和結(jié)構(gòu)材料屬性同聲輻射數(shù)值計(jì)算結(jié)果。內(nèi)聲場(chǎng)中點(diǎn)Pin點(diǎn)為一幅值大小為100 Pa的聲壓激勵(lì),激勵(lì)頻率為125 Hz。本算例以結(jié)構(gòu)外表面輻射聲功率最小化為目標(biāo)函數(shù),設(shè)置目標(biāo)體積為全結(jié)構(gòu)域的80%不變。計(jì)算域用最大0.03 m尺寸的四邊形四節(jié)點(diǎn)等參數(shù)單元離散,過(guò)濾半徑rmin=0.1 m,罰函數(shù)因子p=3,xmin=0.000 1。

    圖4 拱形梁設(shè)計(jì)模型

    為驗(yàn)證本文所提靈敏度方法的設(shè)計(jì)效果,同條件下采用AVM靈敏度法進(jìn)行對(duì)比。迭代過(guò)程的聲功率級(jí)和靈敏度計(jì)算消耗時(shí)間如圖5所示,其中基準(zhǔn)聲功率為1×10-12W。兩種方法的目標(biāo)函數(shù)在迭代過(guò)程中除了在開(kāi)始有小幅波動(dòng)之外,均很快達(dá)到平穩(wěn)狀態(tài),說(shuō)明兩種靈敏度方法的魯棒性較好,且最終的優(yōu)化目標(biāo)函數(shù)一致。對(duì)比兩種方法在迭代過(guò)程中靈敏度部分消耗的時(shí)間可以發(fā)現(xiàn),兩種靈敏度求解方法在每一個(gè)迭代步中的分析時(shí)間沒(méi)有明顯差異。從式(32)和式(36)看出,兩個(gè)靈敏度計(jì)算公式都包含相似的變量,因此計(jì)算時(shí)間沒(méi)有太大差異。

    圖5 目標(biāo)函數(shù)和靈敏度計(jì)算時(shí)間的迭代曲線

    結(jié)構(gòu)對(duì)外聲輻射情況也體現(xiàn)在外部聲場(chǎng)的聲壓級(jí)分布情況上。圖6為優(yōu)化前、后結(jié)構(gòu)拓?fù)鋱D與內(nèi)部聲場(chǎng)和部分外部聲場(chǎng)的聲壓級(jí)響應(yīng)分布,外部聲場(chǎng)截取10 m×6 m范圍。對(duì)比內(nèi)部聲場(chǎng),優(yōu)化前、后僅分布范圍改變,強(qiáng)度未發(fā)生明顯變化,但從外部聲場(chǎng)的聲壓級(jí)分布可以看出,初始設(shè)計(jì)的外聲場(chǎng)存在大范圍的高聲壓級(jí)區(qū)域。優(yōu)化設(shè)計(jì)模型中,該區(qū)域聲壓級(jí)明顯低于初始設(shè)計(jì),這也印證了優(yōu)化后結(jié)構(gòu)的聲輻射得到有效的降低。

    圖6 優(yōu)化前、后結(jié)構(gòu)拓?fù)浜吐晧杭?jí)分布

    3.3 算例2

    算例1為鋁在水介質(zhì)中的聲輻射優(yōu)化,為研究不同聲介質(zhì)和結(jié)構(gòu)材料的優(yōu)化效果,算例2的設(shè)計(jì)模型同算例1,改變聲介質(zhì)和結(jié)構(gòu)材料。初始設(shè)計(jì)域?yàn)槿O(shè)計(jì)域,半圓環(huán)兩側(cè)的底部固定約束,流體下邊界為剛性邊界。圓環(huán)內(nèi)部和外部輻射聲場(chǎng)的聲介質(zhì)均為空氣,密度為1.21 kg/m3,聲速為343 m/s。結(jié)構(gòu)的材料密度為800 kg/m3,彈性模量為2 GPa,泊松比為0.3。聲壓激勵(lì)源位于圓心處Pin點(diǎn),幅值大小為100 Pa,激勵(lì)頻率為205 Hz。以結(jié)構(gòu)外表面輻射聲功率最小化為目標(biāo)函數(shù),設(shè)置目標(biāo)體積為全結(jié)構(gòu)域的65%。單次迭代的體積縮減率為2%,過(guò)濾半徑rmin=0.1 m,罰函數(shù)因子p=3,xmin=0.000 1。

    為方便對(duì)比外聲場(chǎng)情況,截取外部半徑為6 m的圓環(huán)區(qū)域。圖7為優(yōu)化前、后同色度范圍的聲壓級(jí)分布圖和結(jié)構(gòu)拓?fù)湫螒B(tài)。結(jié)構(gòu)拓?fù)涞倪吔缜逦?,且較為規(guī)整,內(nèi)部聲場(chǎng)和外部聲場(chǎng)的聲壓級(jí)得到大幅的降低。這與圖8的迭代曲線相對(duì)應(yīng),結(jié)構(gòu)拓?fù)錆u變穩(wěn)定,迭代曲線未出現(xiàn)大幅波動(dòng),聲功率級(jí)從初始的104.8 dB降低至62.1 dB,降低幅度達(dá)到40.7%。從優(yōu)化前、后內(nèi)部聲場(chǎng)的分布情況看,由激勵(lì)源聲波、邊界反射聲波以及結(jié)構(gòu)振動(dòng)輻射聲波相互疊加形成環(huán)形的低聲壓帶,由于結(jié)構(gòu)拓?fù)涞母淖儯吐晧簬Х植嘉恢煤蛡€(gè)數(shù)發(fā)生改變。為方便觀察內(nèi)部聲壓帶與外部輻射聲場(chǎng)的關(guān)系,圖7還給出聲壓級(jí)高對(duì)比圖,從圖中圈出部分可以看出,圓環(huán)內(nèi)邊界上的聲壓帶個(gè)數(shù)與外聲場(chǎng)高(或低)聲壓是對(duì)應(yīng)的。

    圖7 優(yōu)化前、后結(jié)構(gòu)拓?fù)浜吐晧杭?jí)分布

    圖8 目標(biāo)函數(shù)和約束體積的優(yōu)化迭代曲線

    3.4 算例3

    圖9為算例3的設(shè)計(jì)模型,結(jié)構(gòu)域?yàn)?.0 m×0.2 m的固支梁,梁頂部0.02 m部分為非設(shè)計(jì)域,初始設(shè)計(jì)域即為全設(shè)計(jì)域,結(jié)構(gòu)密度為2 700 kg/m3,彈性模量為69 GPa,泊松比為0.3,設(shè)置目標(biāo)體積為全結(jié)構(gòu)域的85%。結(jié)構(gòu)域下側(cè)聲介質(zhì)域?yàn)?.0 m×0.2 m的矩形域,聲介質(zhì)為水,梁上側(cè)聲介質(zhì)為空氣,各介質(zhì)屬性與上文相同。聲壓激勵(lì)源為矩形流體域的下邊界,幅值大小為100 Pa,激勵(lì)頻率為1 940 Hz。同樣以聲功率最小化為目標(biāo)函數(shù)。設(shè)計(jì)域離散為200×40個(gè)四邊形單元,體積縮減率設(shè)置為1%,過(guò)濾半徑rmin=0.015 m,罰函數(shù)因子p=3,xmin=0.000 1。

    圖9 固支梁設(shè)計(jì)模型

    截取固支梁上方2.0 m×1.0 m區(qū)域的聲壓級(jí)分布來(lái)比較優(yōu)化結(jié)果,如圖10所示??梢钥闯?,優(yōu)化后的外聲場(chǎng)聲壓得到明顯降低。圖中的等值線上數(shù)值表明,初始設(shè)計(jì)的輻射聲場(chǎng)最大聲壓級(jí)達(dá)到100 dB,且占據(jù)固支梁上方的大部分區(qū)域,優(yōu)化構(gòu)型的最大輻射聲壓級(jí)降低至65 dB,降低幅度約35%。圖11顯示了兩種模型在不同激勵(lì)頻率下的聲功率級(jí)曲線,激勵(lì)頻率位于初始設(shè)計(jì)模型的第2個(gè)響應(yīng)峰值處,優(yōu)化后的結(jié)構(gòu)輻射聲功率峰值頻率降低,避開(kāi)激勵(lì)頻率所在位置,使其位于響應(yīng)谷值附近,聲功率級(jí)從初始的70.6 dB降低至54.7 dB,降低幅度達(dá)22.5%。圖12為優(yōu)化迭代曲線,與前文一致,迭代曲線平緩,結(jié)構(gòu)拓?fù)湫问轿窗l(fā)生突變,說(shuō)明本文算法的穩(wěn)定性較好。

    圖10 優(yōu)化前、后結(jié)構(gòu)拓?fù)浜吐晧杭?jí)分布

    圖11 聲功率頻率響應(yīng)曲線

    圖12 目標(biāo)函數(shù)和約束體積的優(yōu)化迭代曲線

    4 結(jié) 論

    1)建立的含內(nèi)部聲場(chǎng)的強(qiáng)耦合有限元/邊界元模型計(jì)算精度較高,可用于求解耦合系統(tǒng)的聲輻射問(wèn)題。

    2)提出的直接求導(dǎo)法與經(jīng)典的伴隨變量法得到的優(yōu)化結(jié)果一致,且魯棒性較好。

    3)本文的線性化剛度模型和靈敏度分析法適用于不同結(jié)構(gòu)材料與聲介質(zhì)問(wèn)題,因此使用本文的優(yōu)化算法可以得到最優(yōu)拓?fù)洹?/p>

    猜你喜歡
    優(yōu)化結(jié)構(gòu)
    超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
    《形而上學(xué)》△卷的結(jié)構(gòu)和位置
    民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
    論結(jié)構(gòu)
    新型平衡塊結(jié)構(gòu)的應(yīng)用
    模具制造(2019年3期)2019-06-06 02:10:54
    論《日出》的結(jié)構(gòu)
    創(chuàng)新治理結(jié)構(gòu)促進(jìn)中小企業(yè)持續(xù)成長(zhǎng)
    亚洲av成人av| 伦理电影免费视频| 中文字幕最新亚洲高清| 色老头精品视频在线观看| 亚洲色图av天堂| 丝袜在线中文字幕| 国产精品国产av在线观看| 精品国产乱码久久久久久男人| 后天国语完整版免费观看| 正在播放国产对白刺激| 成人特级黄色片久久久久久久| 国产主播在线观看一区二区| 狠狠狠狠99中文字幕| 亚洲免费av在线视频| 啦啦啦在线免费观看视频4| 在线观看www视频免费| 午夜福利一区二区在线看| 国产伦人伦偷精品视频| 熟女少妇亚洲综合色aaa.| 国产熟女xx| ponron亚洲| 亚洲久久久国产精品| 日韩高清综合在线| 国产精品一区二区免费欧美| 身体一侧抽搐| 亚洲中文av在线| 亚洲一卡2卡3卡4卡5卡精品中文| 一级毛片高清免费大全| 涩涩av久久男人的天堂| 久久精品国产亚洲av高清一级| 桃色一区二区三区在线观看| 老司机亚洲免费影院| 久久人人精品亚洲av| 90打野战视频偷拍视频| 国产精品亚洲av一区麻豆| 亚洲国产欧美网| 97碰自拍视频| 国产av精品麻豆| 久久久久九九精品影院| 国产av又大| 亚洲欧美精品综合久久99| 一级毛片高清免费大全| 热99国产精品久久久久久7| 麻豆成人av在线观看| 91成人精品电影| 亚洲av成人不卡在线观看播放网| 亚洲av成人不卡在线观看播放网| 搡老岳熟女国产| 免费观看人在逋| 久久久久久亚洲精品国产蜜桃av| 色婷婷久久久亚洲欧美| 在线观看一区二区三区激情| 日本免费a在线| 欧美黄色淫秽网站| 亚洲自偷自拍图片 自拍| 久久人人精品亚洲av| 九色亚洲精品在线播放| 国产精品免费一区二区三区在线| 69精品国产乱码久久久| 久久 成人 亚洲| 欧美日韩亚洲高清精品| 亚洲性夜色夜夜综合| 精品无人区乱码1区二区| 久久久久久久久中文| 麻豆久久精品国产亚洲av | 精品欧美一区二区三区在线| 色婷婷久久久亚洲欧美| 久久精品国产清高在天天线| 国产色视频综合| 黄片播放在线免费| 日韩欧美国产一区二区入口| 51午夜福利影视在线观看| 性欧美人与动物交配| 日韩免费高清中文字幕av| 亚洲成人精品中文字幕电影 | 91老司机精品| 国产深夜福利视频在线观看| 国产成人精品无人区| 久久久国产一区二区| 亚洲国产精品合色在线| 国产精品免费一区二区三区在线| 美女 人体艺术 gogo| 99久久99久久久精品蜜桃| 三上悠亚av全集在线观看| 女生性感内裤真人,穿戴方法视频| 亚洲人成电影免费在线| 国产1区2区3区精品| 亚洲 国产 在线| 大陆偷拍与自拍| 女人精品久久久久毛片| 精品少妇一区二区三区视频日本电影| 国产成年人精品一区二区 | 久久人人97超碰香蕉20202| 久久久久久久午夜电影 | 妹子高潮喷水视频| 一a级毛片在线观看| av在线天堂中文字幕 | 夜夜看夜夜爽夜夜摸 | av网站在线播放免费| 成年版毛片免费区| 美女高潮喷水抽搐中文字幕| 国产欧美日韩综合在线一区二区| 色综合婷婷激情| 最好的美女福利视频网| 午夜亚洲福利在线播放| 国产高清国产精品国产三级| 亚洲全国av大片| 亚洲性夜色夜夜综合| 性欧美人与动物交配| www日本在线高清视频| 在线观看一区二区三区| 国产欧美日韩一区二区精品| 手机成人av网站| 777久久人妻少妇嫩草av网站| 纯流量卡能插随身wifi吗| 亚洲精品一区av在线观看| 在线视频色国产色| 在线视频色国产色| 精品人妻在线不人妻| avwww免费| 99在线人妻在线中文字幕| 狠狠狠狠99中文字幕| av福利片在线| 高清在线国产一区| 精品卡一卡二卡四卡免费| 亚洲专区国产一区二区| 亚洲精品一区av在线观看| www.999成人在线观看| 久久精品国产亚洲av高清一级| 精品久久久久久电影网| 亚洲avbb在线观看| 欧美乱色亚洲激情| 80岁老熟妇乱子伦牲交| 国产极品粉嫩免费观看在线| 女性被躁到高潮视频| 在线永久观看黄色视频| 国产欧美日韩综合在线一区二区| 曰老女人黄片| 久久香蕉国产精品| 高清av免费在线| 黑人巨大精品欧美一区二区mp4| 亚洲精品美女久久久久99蜜臀| 一进一出抽搐动态| 手机成人av网站| 国产91精品成人一区二区三区| 99精品久久久久人妻精品| 亚洲av片天天在线观看| 一进一出好大好爽视频| 老司机午夜十八禁免费视频| 国产精品98久久久久久宅男小说| 女性生殖器流出的白浆| 自拍欧美九色日韩亚洲蝌蚪91| 一级a爱视频在线免费观看| 人人妻人人爽人人添夜夜欢视频| 老司机午夜十八禁免费视频| 夜夜看夜夜爽夜夜摸 | 在线观看日韩欧美| 亚洲国产精品合色在线| 国产三级黄色录像| 久久久久国产一级毛片高清牌| 亚洲五月色婷婷综合| 久久精品影院6| 亚洲成a人片在线一区二区| 老司机在亚洲福利影院| av有码第一页| 免费av中文字幕在线| 悠悠久久av| 叶爱在线成人免费视频播放| 久久国产精品人妻蜜桃| 久久人人精品亚洲av| 亚洲精华国产精华精| 免费观看精品视频网站| 久久国产精品男人的天堂亚洲| 夜夜爽天天搞| netflix在线观看网站| 国产精品香港三级国产av潘金莲| 久久午夜亚洲精品久久| 热99国产精品久久久久久7| a级片在线免费高清观看视频| 在线观看免费午夜福利视频| 久久久久国产精品人妻aⅴ院| 99精品欧美一区二区三区四区| 精品一区二区三区四区五区乱码| 国产精品1区2区在线观看.| 国产精品一区二区三区四区久久 | 欧美日韩中文字幕国产精品一区二区三区 | 国产精品 国内视频| 久久青草综合色| 亚洲成人精品中文字幕电影 | 久久人人爽av亚洲精品天堂| 一区二区三区激情视频| 天堂中文最新版在线下载| 日韩有码中文字幕| 一本综合久久免费| 欧美黑人欧美精品刺激| 91麻豆av在线| 亚洲 国产 在线| 国产成人精品久久二区二区91| 久久热在线av| 日韩精品免费视频一区二区三区| 久久久国产成人免费| 高清欧美精品videossex| 国产91精品成人一区二区三区| 一进一出抽搐动态| 日本黄色视频三级网站网址| 欧美成人性av电影在线观看| 热re99久久国产66热| 男人舔女人下体高潮全视频| 日本wwww免费看| 亚洲av五月六月丁香网| 午夜福利免费观看在线| 国产精品一区二区免费欧美| 国产不卡一卡二| 国产av又大| 99国产精品一区二区三区| 99精品欧美一区二区三区四区| 性欧美人与动物交配| 亚洲熟女毛片儿| 国产av精品麻豆| 国产精品一区二区在线不卡| 在线观看一区二区三区激情| 日本wwww免费看| 精品少妇一区二区三区视频日本电影| 亚洲av第一区精品v没综合| 最新美女视频免费是黄的| 国产亚洲精品第一综合不卡| 后天国语完整版免费观看| 欧美国产精品va在线观看不卡| 岛国视频午夜一区免费看| 久久久久久亚洲精品国产蜜桃av| 久久久久久人人人人人| 午夜福利在线观看吧| 国产成人精品久久二区二区免费| 美女高潮喷水抽搐中文字幕| 日本三级黄在线观看| 亚洲av成人不卡在线观看播放网| 如日韩欧美国产精品一区二区三区| 热99国产精品久久久久久7| 99国产精品免费福利视频| 成年女人毛片免费观看观看9| 亚洲精品在线美女| 国产精品 国内视频| 精品卡一卡二卡四卡免费| 在线永久观看黄色视频| 成人三级做爰电影| 女人高潮潮喷娇喘18禁视频| av福利片在线| 国产无遮挡羞羞视频在线观看| 欧美日韩乱码在线| 一区二区三区国产精品乱码| 国产精品二区激情视频| 亚洲欧洲精品一区二区精品久久久| 婷婷精品国产亚洲av在线| 久久久久国内视频| 又黄又爽又免费观看的视频| 极品教师在线免费播放| 在线观看免费高清a一片| 午夜免费观看网址| 亚洲男人的天堂狠狠| 日韩欧美国产一区二区入口| 日韩中文字幕欧美一区二区| 一级片免费观看大全| 亚洲午夜理论影院| 国产欧美日韩一区二区三| 欧美日韩乱码在线| 国产极品粉嫩免费观看在线| 中文字幕高清在线视频| 黄色成人免费大全| 亚洲国产看品久久| 性色av乱码一区二区三区2| 久久亚洲真实| 午夜精品久久久久久毛片777| 黄色a级毛片大全视频| 女人高潮潮喷娇喘18禁视频| 视频区图区小说| 天堂俺去俺来也www色官网| 天堂中文最新版在线下载| 丁香六月欧美| 欧美日本中文国产一区发布| 淫妇啪啪啪对白视频| av欧美777| 免费在线观看黄色视频的| 一二三四在线观看免费中文在| 亚洲九九香蕉| 亚洲国产欧美网| 十八禁网站免费在线| 国产精品久久电影中文字幕| 国产xxxxx性猛交| 国产麻豆69| 在线观看www视频免费| av网站免费在线观看视频| 好看av亚洲va欧美ⅴa在| 亚洲av成人av| 精品国产乱子伦一区二区三区| 80岁老熟妇乱子伦牲交| 日韩欧美在线二视频| 国产三级在线视频| 免费在线观看亚洲国产| 好男人电影高清在线观看| 久久性视频一级片| 丰满饥渴人妻一区二区三| 国产成人精品无人区| 女人爽到高潮嗷嗷叫在线视频| 一二三四社区在线视频社区8| 久久伊人香网站| 亚洲自拍偷在线| 午夜老司机福利片| 老司机靠b影院| 久久精品成人免费网站| 久久国产精品影院| 中亚洲国语对白在线视频| 免费日韩欧美在线观看| 亚洲欧美一区二区三区黑人| 男男h啪啪无遮挡| 精品乱码久久久久久99久播| 国产日韩一区二区三区精品不卡| 久久中文看片网| 夫妻午夜视频| 久久精品国产亚洲av香蕉五月| 亚洲一区高清亚洲精品| 乱人伦中国视频| 久久香蕉激情| 午夜视频精品福利| 91国产中文字幕| 午夜激情av网站| 99精国产麻豆久久婷婷| 91在线观看av| 亚洲人成77777在线视频| 亚洲精品av麻豆狂野| 午夜福利,免费看| 丰满人妻熟妇乱又伦精品不卡| 亚洲少妇的诱惑av| bbb黄色大片| 亚洲精品国产区一区二| 午夜精品在线福利| 亚洲国产欧美日韩在线播放| 久9热在线精品视频| 久久精品91蜜桃| 又大又爽又粗| 国产亚洲精品久久久久久毛片| 交换朋友夫妻互换小说| 曰老女人黄片| 亚洲五月天丁香| 欧美人与性动交α欧美软件| 日日干狠狠操夜夜爽| 中文字幕精品免费在线观看视频| 男女之事视频高清在线观看| 色婷婷久久久亚洲欧美| 亚洲黑人精品在线| 手机成人av网站| а√天堂www在线а√下载| 99久久99久久久精品蜜桃| 国产精品九九99| 91老司机精品| 亚洲精品美女久久久久99蜜臀| 黄片播放在线免费| 亚洲精品美女久久av网站| 黑人欧美特级aaaaaa片| 免费高清视频大片| 可以在线观看毛片的网站| 欧美黄色片欧美黄色片| xxx96com| www.999成人在线观看| 精品久久久久久久久久免费视频 | 露出奶头的视频| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品sss在线观看 | 一边摸一边做爽爽视频免费| 超色免费av| 一边摸一边做爽爽视频免费| 久久九九热精品免费| 午夜两性在线视频| 国产色视频综合| 激情视频va一区二区三区| 十八禁人妻一区二区| 亚洲五月色婷婷综合| 亚洲一区二区三区欧美精品| 夜夜看夜夜爽夜夜摸 | 色老头精品视频在线观看| 亚洲av熟女| 又黄又爽又免费观看的视频| 夜夜爽天天搞| 亚洲成人免费电影在线观看| 亚洲色图 男人天堂 中文字幕| 精品国产超薄肉色丝袜足j| 久久九九热精品免费| 亚洲精品国产区一区二| 在线观看免费午夜福利视频| 女人高潮潮喷娇喘18禁视频| 精品午夜福利视频在线观看一区| 午夜a级毛片| 日韩欧美在线二视频| 国产精品国产av在线观看| 国产高清国产精品国产三级| 一二三四社区在线视频社区8| 欧美最黄视频在线播放免费 | 精品第一国产精品| 99精品欧美一区二区三区四区| 黄色怎么调成土黄色| 国产欧美日韩精品亚洲av| 我的亚洲天堂| 欧美日韩中文字幕国产精品一区二区三区 | 久久久久精品国产欧美久久久| 久久天堂一区二区三区四区| 日韩大尺度精品在线看网址 | 黄色片一级片一级黄色片| 国产成人欧美| 18禁美女被吸乳视频| 黄色丝袜av网址大全| 色哟哟哟哟哟哟| 一区在线观看完整版| 日本免费一区二区三区高清不卡 | 日日摸夜夜添夜夜添小说| 女性生殖器流出的白浆| 又黄又粗又硬又大视频| 国产精品 国内视频| 午夜两性在线视频| 欧美日韩福利视频一区二区| 国产av在哪里看| 国产精品秋霞免费鲁丝片| 婷婷精品国产亚洲av在线| 久久久国产成人精品二区 | 天堂影院成人在线观看| 十八禁网站免费在线| 91精品国产国语对白视频| 99精国产麻豆久久婷婷| 激情在线观看视频在线高清| 国产精华一区二区三区| 国产欧美日韩一区二区三区在线| 在线av久久热| 国产高清激情床上av| 级片在线观看| www.999成人在线观看| 窝窝影院91人妻| 日韩视频一区二区在线观看| 欧美日韩瑟瑟在线播放| www.www免费av| 精品一区二区三区视频在线观看免费 | 国产精品一区二区在线不卡| 美女大奶头视频| 日本三级黄在线观看| 亚洲,欧美精品.| 国产激情欧美一区二区| av在线天堂中文字幕 | 久久中文字幕人妻熟女| 国产精品国产高清国产av| 国产黄色免费在线视频| 老司机靠b影院| 极品人妻少妇av视频| 日本欧美视频一区| 久久伊人香网站| 丰满的人妻完整版| 淫妇啪啪啪对白视频| 国产免费现黄频在线看| 热99国产精品久久久久久7| 一级片免费观看大全| 精品国产超薄肉色丝袜足j| 国产成人av教育| 精品久久久久久,| 欧美中文综合在线视频| 亚洲av美国av| 亚洲av片天天在线观看| cao死你这个sao货| 精品一区二区三区视频在线观看免费 | 美女高潮喷水抽搐中文字幕| 三级毛片av免费| 黄色 视频免费看| 十八禁人妻一区二区| 黄色成人免费大全| 人成视频在线观看免费观看| 又黄又粗又硬又大视频| 免费观看人在逋| 色婷婷久久久亚洲欧美| 少妇裸体淫交视频免费看高清 | 真人做人爱边吃奶动态| 国产欧美日韩精品亚洲av| 精品久久久久久久久久免费视频 | 91国产中文字幕| 久久香蕉精品热| 久久天躁狠狠躁夜夜2o2o| 午夜激情av网站| 日韩av在线大香蕉| 超碰成人久久| 中文字幕人妻丝袜制服| 亚洲色图综合在线观看| 美女大奶头视频| 日本wwww免费看| 国产精品免费视频内射| a级毛片在线看网站| 国产精品一区二区精品视频观看| 少妇裸体淫交视频免费看高清 | 看免费av毛片| 91精品三级在线观看| 黄色a级毛片大全视频| 黄片小视频在线播放| 级片在线观看| 1024香蕉在线观看| 成人18禁高潮啪啪吃奶动态图| 欧美午夜高清在线| 级片在线观看| 国产一卡二卡三卡精品| 成人影院久久| 18禁国产床啪视频网站| 最新美女视频免费是黄的| 91麻豆精品激情在线观看国产 | 一进一出好大好爽视频| www.熟女人妻精品国产| 久久精品国产亚洲av香蕉五月| 国产激情久久老熟女| 80岁老熟妇乱子伦牲交| 黄色丝袜av网址大全| 国产亚洲欧美在线一区二区| 老司机在亚洲福利影院| 操出白浆在线播放| 亚洲国产看品久久| 啦啦啦在线免费观看视频4| 日日干狠狠操夜夜爽| 亚洲精品国产区一区二| 久久草成人影院| 人人妻人人爽人人添夜夜欢视频| 我的亚洲天堂| 国产精品日韩av在线免费观看 | 久久精品国产综合久久久| 叶爱在线成人免费视频播放| 国产伦一二天堂av在线观看| 日韩欧美免费精品| 色婷婷久久久亚洲欧美| av网站免费在线观看视频| 又紧又爽又黄一区二区| 一个人观看的视频www高清免费观看 | 国产av一区在线观看免费| 亚洲av日韩精品久久久久久密| 久久精品亚洲熟妇少妇任你| 大型av网站在线播放| 国产精品乱码一区二三区的特点 | 91麻豆av在线| 亚洲精品美女久久av网站| 女性被躁到高潮视频| 69精品国产乱码久久久| 午夜a级毛片| 真人做人爱边吃奶动态| 久久天躁狠狠躁夜夜2o2o| 韩国精品一区二区三区| 亚洲欧洲精品一区二区精品久久久| 一个人免费在线观看的高清视频| 精品少妇一区二区三区视频日本电影| 日本vs欧美在线观看视频| 国产一卡二卡三卡精品| 国产伦一二天堂av在线观看| 香蕉国产在线看| 国产男靠女视频免费网站| 欧美中文日本在线观看视频| 一a级毛片在线观看| 一级毛片精品| 色精品久久人妻99蜜桃| 国产男靠女视频免费网站| 丰满迷人的少妇在线观看| 亚洲精品美女久久av网站| 91在线观看av| 十八禁人妻一区二区| 日韩精品青青久久久久久| 中出人妻视频一区二区| 国产精品国产av在线观看| 人人妻人人添人人爽欧美一区卜| 国产三级黄色录像| 91精品国产国语对白视频| 黑丝袜美女国产一区| 91在线观看av| 大型av网站在线播放| 黑人巨大精品欧美一区二区mp4| 如日韩欧美国产精品一区二区三区| 国产一区二区三区视频了| 国产av又大| 91字幕亚洲| 久久久精品国产亚洲av高清涩受| 女人被躁到高潮嗷嗷叫费观| 中文字幕人妻熟女乱码| 波多野结衣高清无吗| 免费在线观看亚洲国产| 另类亚洲欧美激情| 国产成+人综合+亚洲专区| 麻豆国产av国片精品| av欧美777| 精品国产国语对白av| 中文欧美无线码| 午夜福利欧美成人| 色播在线永久视频| 国产色视频综合| 国产不卡一卡二| 色播在线永久视频| 香蕉国产在线看| 精品一品国产午夜福利视频| 久久精品国产综合久久久| 欧美成狂野欧美在线观看| 美女午夜性视频免费| 看黄色毛片网站| 青草久久国产| av欧美777| 在线观看舔阴道视频| 亚洲色图 男人天堂 中文字幕| 日本 av在线| 精品久久久精品久久久| 日韩高清综合在线| 国产xxxxx性猛交| 在线国产一区二区在线| 青草久久国产| 国产主播在线观看一区二区| 女人爽到高潮嗷嗷叫在线视频| 国产欧美日韩一区二区精品| av视频免费观看在线观看| 欧美人与性动交α欧美软件| 一进一出抽搐动态| 欧美黄色片欧美黄色片| 色婷婷久久久亚洲欧美| 精品久久久久久电影网|