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

    基于有限元—統(tǒng)計(jì)能量數(shù)值混響室法的聲學(xué)包設(shè)計(jì)與優(yōu)化

    2017-08-05 01:36:48楊德慶石嘉欣郁揚(yáng)
    中國(guó)艦船研究 2017年4期
    關(guān)鍵詞:有限元優(yōu)化設(shè)計(jì)

    楊德慶,石嘉欣,郁揚(yáng)

    1高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海200240 2上海交通大學(xué)海洋工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海200240 3上海交通大學(xué)船舶海洋與建筑工程學(xué)院,上海200240

    基于有限元—統(tǒng)計(jì)能量數(shù)值混響室法的聲學(xué)包設(shè)計(jì)與優(yōu)化

    楊德慶1,2,3,石嘉欣1,2,3,郁揚(yáng)1,2,3

    1高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海200240 2上海交通大學(xué)海洋工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海200240 3上海交通大學(xué)船舶海洋與建筑工程學(xué)院,上海200240

    [目的]為解決傳統(tǒng)聲學(xué)包設(shè)計(jì)法耗時(shí)、耗資高的問(wèn)題,[方法]利用有限元—統(tǒng)計(jì)能量(FE-SEA)數(shù)值混響室法在全頻段進(jìn)行聲學(xué)包的設(shè)計(jì)與優(yōu)化。設(shè)計(jì)空氣層聲學(xué)包及無(wú)空氣層聲學(xué)包這2種類型共8種聲學(xué)包,計(jì)算其隔聲量。采用遺傳算法,以達(dá)到最大隔聲量值為目標(biāo)對(duì)每一種聲學(xué)包的材料層厚度進(jìn)行優(yōu)化。[結(jié)果]研究表明,有空氣層存在的聲學(xué)包的隔聲量值較無(wú)空氣層聲學(xué)包的高2~6 dB,且與原始設(shè)計(jì)相比,優(yōu)化后的聲學(xué)包的隔聲性能得到大幅提高。[結(jié)論]所得研究成果對(duì)實(shí)際工程中船舶聲學(xué)包設(shè)計(jì)有重要的指導(dǎo)價(jià)值。

    聲學(xué)包;優(yōu)化;FS-SEA數(shù)值混響室法;隔聲量

    0 引 言

    隨著新能源及綠色船舶設(shè)計(jì)技術(shù)的發(fā)展,噪聲作為影響舒適性的重要因素,越來(lái)越受到人們的重視。國(guó)際海事組織和各國(guó)船級(jí)社也陸續(xù)出臺(tái)了嚴(yán)格的船舶噪聲規(guī)范,降噪成為船舶行業(yè)的技術(shù)熱點(diǎn)。聲學(xué)包作為解決降噪問(wèn)題的新技術(shù)和新產(chǎn)品,受到廣泛關(guān)注。所謂聲學(xué)包,是指將阻尼材料、吸聲材料、隔聲材料以及防火材料等通過(guò)裝飾面板、襯板或封裝面板,制作成一定形狀的層合結(jié)構(gòu),以實(shí)現(xiàn)降低指定頻段范圍內(nèi)一定量級(jí)噪聲及振動(dòng)的聲學(xué)處理元件。汽車行業(yè)是較早開始采用聲學(xué)包技術(shù)來(lái)提升汽車聲學(xué)環(huán)境舒適性的一個(gè)領(lǐng)域[1-4]。本文擬利用數(shù)值模擬手段研究聲學(xué)包的設(shè)計(jì)和優(yōu)化方法,探索聲學(xué)包在船舶領(lǐng)域的應(yīng)用。數(shù)值模擬是聲學(xué)包設(shè)計(jì)的重要手段,可規(guī)避實(shí)驗(yàn)測(cè)試的成本高、周期長(zhǎng)等問(wèn)題,在理論研究及產(chǎn)品開發(fā)階段更能縮短設(shè)計(jì)周期,適于產(chǎn)品多樣性設(shè)計(jì)。范明偉[5]和高處[6]利用統(tǒng)計(jì)能量(SEA)數(shù)值混響室法進(jìn)行了聲學(xué)包設(shè)計(jì),這種方法的缺點(diǎn)在于單板及兩側(cè)聲腔子系統(tǒng)在低頻段(100 Hz以下)不滿足SEA法對(duì)子系統(tǒng)單位帶寬模態(tài)數(shù)大于5的要求,其計(jì)算結(jié)果在低頻段不精確、不可信。在聲學(xué)包的優(yōu)化方面,蘇朝勇[7]利用SEA數(shù)值混響室法對(duì)聲學(xué)包進(jìn)行了優(yōu)化設(shè)計(jì),但在低頻段仍然存在計(jì)算結(jié)果不可信的問(wèn)題。常規(guī)聲學(xué)包主要用于中、高頻的降噪,但在降噪設(shè)計(jì)中聲學(xué)包的低頻性能也是需要考慮的,不能總是忽略不計(jì)。在有些情況下,需要在聲學(xué)包內(nèi)加入阻尼材料來(lái)進(jìn)行低頻減振降噪設(shè)計(jì),這時(shí)就需進(jìn)行低頻段計(jì)算,所以未來(lái)的聲學(xué)包設(shè)計(jì)不應(yīng)僅針對(duì)中、高頻,還需覆蓋全頻段。

    本文將采用適用于全頻段的有限元—統(tǒng)計(jì)能量(FE-SEA)數(shù)值混響室法進(jìn)行聲學(xué)包的設(shè)計(jì)[8],將其作為評(píng)價(jià)聲學(xué)包性能的參數(shù),得到聲學(xué)包的隔聲量。FE-SEA數(shù)值混響室法可以很好地解決低頻域聲學(xué)包材料單位帶寬模態(tài)數(shù)不夠的問(wèn)題,從而得到較精確的聲學(xué)包全頻段的隔聲量。此外,還將利用遺傳算法對(duì)含空氣層與不含空氣層這2種類型的聲學(xué)包進(jìn)行優(yōu)化設(shè)計(jì),提升聲學(xué)包的性能,達(dá)到更佳的降噪效果。同時(shí),分別探討在約束聲學(xué)包整體厚度情況下聲學(xué)包的優(yōu)化效果。所得結(jié)果用于為實(shí)際工程設(shè)計(jì)提供參考。

    1 聲學(xué)性能數(shù)值分析理論與方法

    聲學(xué)包的聲學(xué)性能是指在某一頻段聲學(xué)包可以實(shí)現(xiàn)的降噪指標(biāo)或振動(dòng)衰減量,如隔聲量或插入損失。聲學(xué)包聲學(xué)性能評(píng)定的常規(guī)手段是聲學(xué)實(shí)驗(yàn)測(cè)試,如混響室法和駐波管法。本文采用FE-SEA數(shù)值混響室法代替聲學(xué)實(shí)驗(yàn)測(cè)試,具體原理如下。

    1.1 FE-SEA耦合聲學(xué)分析原理

    FE-SEA耦合分析方法將聲場(chǎng)分為直接場(chǎng)和混響場(chǎng),前者定義為由入射波產(chǎn)生的聲場(chǎng),后者指經(jīng)過(guò)一次以上反射形成的聲場(chǎng)。結(jié)構(gòu)響應(yīng)q可以寫為以下形式[9-10]:

    子系統(tǒng)j輸入到直接場(chǎng)的平均功率為:

    以上式中:ω為圓頻率;Pin,j為子系統(tǒng)j輸入到直接場(chǎng)的平均功率;為外載荷輸入到子系統(tǒng)j的功率;Pout,j為子系統(tǒng)j輸出到混響場(chǎng)的平均功率;Sqq為響應(yīng)q的互譜矩陣,系統(tǒng)有限元模型的動(dòng)力剛度矩陣;Sff為外載荷的互譜矩陣;D(k)dir為有限元與統(tǒng)計(jì)能量模型連接處的動(dòng)力剛度矩陣;Dd,rs為系統(tǒng)有限元模型的動(dòng)力剛度矩陣;Ek,nk分別為子系統(tǒng)k的振動(dòng)能量和模態(tài)密度;ηj為子系統(tǒng)j的損耗因子;ηjk為統(tǒng)計(jì)能量模型的耦合損耗因子;ηd,j為有限元模型計(jì)算的損耗因子;下標(biāo)rs為隨機(jī)子系統(tǒng)j確定性邊界上的自由度。其中式(7)和式(8)是FE-SEA混合分析理論的主要方程。

    1.2 FE-SEA數(shù)值混響室法

    FE-SEA數(shù)值混響室法[8]簡(jiǎn)稱為FS數(shù)值混響室法,是通過(guò)建立構(gòu)件的有限元模型,在構(gòu)件的一側(cè)建立半無(wú)限流場(chǎng)(SIF)并與構(gòu)件聯(lián)接來(lái)模擬聲學(xué)及結(jié)構(gòu)特性。在構(gòu)件另一側(cè),施加混響激勵(lì)以模擬實(shí)際聲學(xué)邊界條件,并在有限元構(gòu)件處施加簡(jiǎn)支邊界條件。試件尺寸可根據(jù)實(shí)際結(jié)構(gòu)設(shè)置,若無(wú)具體要求,在評(píng)價(jià)某材料聲學(xué)性能時(shí)建議選取測(cè)試件大小為350 mm×220 mm的長(zhǎng)方形結(jié)構(gòu),半無(wú)限流場(chǎng)與待測(cè)試件之間的距離建議取為500mm。作為入射聲源的混響聲源的聲壓值,應(yīng)在各測(cè)試頻率下均相同且穩(wěn)定,建議取為1 Pa(等效于94 dB)。求解該FE-SEA耦合模型,通過(guò)入射聲源激勵(lì)和半無(wú)限流場(chǎng)處的聲壓值,可以得到待測(cè)試件的隔聲量。上述算法實(shí)現(xiàn)及模型建立是采用VA One 2013軟件完成的,原理如圖1所示。

    圖1 FS數(shù)值混響室法示意圖Fig.1 Schematic of FS numerical reverberation method

    2 聲學(xué)包隔聲量的FE-SEA數(shù)值混響室法計(jì)算

    2.1 FE-SEA耦合聲學(xué)分析算例

    建立如圖1所示的FS數(shù)值混響室法模型,構(gòu)件為鋁板,尺寸為350 mm×220 mm×1 mm,彈性模量為71.0 GPa,泊松比為0.33,密度為2.81×103kg/m3,阻尼損耗系數(shù)為0.1%。半無(wú)限流場(chǎng)距鋁板500 mm,混響聲源強(qiáng)度的取值為1 Pa(等效于94 dB)。結(jié)構(gòu)約束邊界條件為簡(jiǎn)支,計(jì)算頻率16~125 Hz。低頻段(125 Hz以下)FS數(shù)值混響室法與瑞利李茲理論解法的隔聲量結(jié)果[8]如圖2所示。

    計(jì)算結(jié)果表明,125 Hz以下FS數(shù)值混響室法得到的鋁板的隔聲量與瑞利李茲理論計(jì)算方法所得結(jié)果吻合較好。而有關(guān)125 Hz以上該模型的準(zhǔn)確性,文獻(xiàn)[11]已進(jìn)行了驗(yàn)證,證明其隔聲量計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果吻合較好,F(xiàn)S數(shù)值混響室法可用于全頻段隔聲量的精確計(jì)算。

    圖2 低頻段(125 Hz以下)FS數(shù)值混響室法模型驗(yàn)證Fig.2 The modal validation of FS numerical reverberation method(below 125 Hz)

    2.2 艙室降噪聲學(xué)包設(shè)計(jì)與性能評(píng)估

    采用1.2節(jié)所述的FS數(shù)值混響室模型進(jìn)行聲學(xué)包性能分析。本文設(shè)計(jì)的聲學(xué)包主要用于船舶艙室降噪,選用鋼板作為襯板,尺寸為350 mm×220 mm×7 mm,彈性模量為210.0 GPa,泊松比為0.3,密度為7.8×103kg/m3,損耗系數(shù)為0.1%。半無(wú)限流場(chǎng)距襯板500 mm,混響聲源強(qiáng)度取為1 Pa,聲學(xué)包邊界條件為四邊簡(jiǎn)支約束。鑒于船舶降噪設(shè)計(jì)的主要頻段為中、高頻,因此計(jì)算頻段取為63~8 000 Hz。分別設(shè)計(jì)8種聲學(xué)包(Case 1~Case 8),聲學(xué)包假設(shè)設(shè)計(jì)為7層(Layer 1~Layer 7),其材料參數(shù)和結(jié)構(gòu)參數(shù)分別如表1及表2所示[12]。計(jì)算得到的8種聲學(xué)包的隔聲量曲線及總隔聲量(Transmission Loss,TL)如圖3和圖4所示。

    其中,Case 1和Case 2,Case 3和Case 4,Case 5和Case 6,Case 7和Case 8兩兩為一組,其材料設(shè)置均相同,差別在于后者在層間插入了5 mm厚的空氣層。也即Case 1,Case 3,Case 5和Case 7為未加空氣層的聲學(xué)包;Case 2,Case 4,Case 6和Case 8為加入了空氣層的聲學(xué)包。

    表1 8種聲學(xué)包設(shè)計(jì)方案的材料特性參數(shù)Table 1 Material characteristic parameters of eight sound packages

    表2 8種聲學(xué)包設(shè)計(jì)方案的結(jié)構(gòu)參數(shù)Table 2 Structural characteristic parameters of eight sound packages

    圖3 8種聲學(xué)包的隔聲量曲線Fig.3 The TL curves of eight sound packages

    圖4 8種聲學(xué)包總隔聲量曲線Fig.4 The overall TL values of eight sound packages

    計(jì)算表明:8種聲學(xué)包均可將總隔聲量再提高1.5~9 dB,其中Case 4的隔聲量最大,性能最佳。

    3 聲學(xué)包聲學(xué)性能優(yōu)化設(shè)計(jì)

    以2.2節(jié)中的8種聲學(xué)包為對(duì)象,采用遺傳算法進(jìn)行優(yōu)化。以獲得最大隔聲量為優(yōu)化目標(biāo),約束條件為材料層(即聲學(xué)包)總厚度,聲學(xué)包的材料層厚度為設(shè)計(jì)變量,變量下限為0,上限為整個(gè)聲學(xué)包的厚度。按照約束條件的不同,分2種情況進(jìn)行計(jì)算:優(yōu)化算例1,約束設(shè)定為優(yōu)化后聲學(xué)包總厚度不大于原始設(shè)計(jì)厚度;優(yōu)化算例2,優(yōu)化后聲學(xué)包總厚度不大于150 mm。

    3.1 優(yōu)化算例1

    優(yōu)化時(shí),要求聲學(xué)包的總厚度不大于原始設(shè)計(jì)的厚度。各方案優(yōu)化前、后的參數(shù)對(duì)比如表3所示,優(yōu)化前、后隔聲量曲線的對(duì)比結(jié)果如圖5所示。

    表3和圖5共同表明,在控制聲學(xué)包厚度的基礎(chǔ)上:

    1)Case 1優(yōu)化后,總隔聲量上升了1.55 dB,總厚度減小了0.5mm,總重量減小了0.36×10-6kg/m2;Case 2優(yōu)化后,總隔聲量整體上升了22.56 dB,和初始設(shè)計(jì)相比提高了94.8%,總厚度減小了3.1 mm,重量降低為原始設(shè)計(jì)的34%;Case 3優(yōu)化后,總隔聲量整體上升了1.38 dB,厚度減小了0.41 mm,重量與初始設(shè)計(jì)時(shí)的基本相同;Case 4優(yōu)化后,總隔聲量整體上升了18 dB,和原始設(shè)計(jì)相比提高了71%,厚度減小了0.1 mm,重量降低為原始設(shè)計(jì)的34.8%;Case 5優(yōu)化后,總隔聲量整體上升了1.23 dB,厚度減小了7.32 mm,重量減小了0.4×10-6kg/m2;Case 6優(yōu)化后,總隔聲量整體上升了20.36 dB,和原始設(shè)計(jì)相比提高了86.9%,厚度減小了5.21 mm,重量變?yōu)樵瓉?lái)的34%;Case 7優(yōu)化后,總隔聲量整體上升了1.68 dB,厚度減小了0.75 mm,重量和原始設(shè)計(jì)相比略有增重;Case 8優(yōu)化后,總隔聲量整體上升了19.87 dB,和原始設(shè)計(jì)相比提高了100%,厚度減小了9.87 mm,重量變?yōu)樵荚O(shè)計(jì)的22.5%。

    2)不含空氣層聲學(xué)包(Case1,Case 3,Case 5和Case 7)優(yōu)化后,隔聲量平均上升了1.46 dB,厚度略有減小,重量與原始設(shè)計(jì)相比大致持平,說(shuō)明不含空氣層聲學(xué)包的聲學(xué)性能會(huì)有一定的提升,但提升的潛力較小。

    3)含空氣層的聲學(xué)包(Case2,Case 4,Case 6和Case 8)優(yōu)化后,隔聲量平均上升20 dB以上,和原始設(shè)計(jì)相比平均提升了88%,厚度平均減小了4.57 mm,并且由于空氣層厚度的增加,使得聲學(xué)包重量大幅下降,平均降低為原始設(shè)計(jì)的31.3%。由此可見(jiàn),含空氣層聲學(xué)包的性能提升較大。

    4)從優(yōu)化后的結(jié)果來(lái)看,含空氣層聲學(xué)包(Case2,Case 4,Case 6和Case 8)在優(yōu)化后空氣層的厚度會(huì)大幅增加,由此可見(jiàn)空氣層在提升聲學(xué)包性能方面具有重要作用。

    3.2 優(yōu)化算例2

    要求聲學(xué)包的總厚度不大于150 mm。各方案優(yōu)化前、后參數(shù)的對(duì)比如表4所示,優(yōu)化前、后隔聲量曲線的對(duì)比結(jié)果如圖6所示。

    表3 Case 1~Case 8優(yōu)化前后變量值對(duì)比(算例1)Table 3 Comparison of parameters before and after optimization of Case 1~Case 8 in example 1

    表4 Case 1~Case 8優(yōu)化前后變量值對(duì)比(算例2)Table 4 Comparison of parameters before and after optimization of Case 1~Case 8 in example 2

    圖5 Case 1~Case 8優(yōu)化前后隔聲量曲線對(duì)比(算例1)Fig.5 Comparison of TL curves before and after optimization of Case 1~Case 8 in example 1

    圖6 Case 1~Case 8優(yōu)化前后隔聲量曲線對(duì)比(算例2)Fig 6 Comparison of TL curves before and after optimization of Case 1~Case 8 in example 2

    表4和圖6共同表明,聲學(xué)包厚度控制在150 mm以內(nèi)時(shí):

    1)Case 1優(yōu)化后總隔聲量上升了6.5 dB,厚度增加了37.3 mm,相比于初始設(shè)計(jì)重量增加了23.7%。Case 2優(yōu)化后總隔聲量整體上升了27 dB,厚度增加了37.4 mm,重量降低為原始設(shè)計(jì)的58.8%。Case 3優(yōu)化后總隔聲量上升了8 dB,厚度增加50 mm,重量和原始設(shè)計(jì)相比增加了36.4%。Case 4優(yōu)化后總隔聲量整體上升了28.8 dB,厚度增加48 mm,重量下降為原始設(shè)計(jì)的79.7%。Case 5優(yōu)化后總隔聲量整體上升了7.6 dB,厚度增加59.2 mm,重量和原始設(shè)計(jì)相比增加了47.8%。Case 6優(yōu)化后總隔聲量整體上升了24 dB,厚度增加49.3 mm,重量降低為初始設(shè)計(jì)的58%。Case 7優(yōu)化后總隔聲量整體上升了8.2 dB,厚度增加54.2 mm,相比于原始設(shè)計(jì)重量增加了39.6%。Case 8優(yōu)化后總隔聲量整體上升了26.3 dB,厚度增加60.8 mm,重量降低為原始設(shè)計(jì)的51.9%。

    2)不含空氣層的聲學(xué)包(Case 1,Case 3,Case 5和Case 7)優(yōu)化后,隔聲量平均上升了7.6 dB,最高(Case 7)可提高8.21 dB,厚度平均增加50 mm,重量相對(duì)原始設(shè)計(jì)平均增加了36.9%。這說(shuō)明當(dāng)聲學(xué)包厚度增加時(shí),不含空氣層聲學(xué)包的性能會(huì)進(jìn)一步提升;與3.1節(jié)中數(shù)據(jù)進(jìn)行對(duì)比可以看出,隔聲量雖然上升了約6.2 dB,但厚度增加了50 mm左右,重量增加也超過(guò)原始設(shè)計(jì)的35%以上,說(shuō)明更好的聲學(xué)性能是以犧牲厚度和重量為代價(jià)的。

    3)含空氣層的聲學(xué)包(Case 2,Case 4,Case 6和Case 8)在優(yōu)化后,隔聲量平均上升了26.5 dB以上,最高(Case 4)可達(dá)28.81 dB,厚度平均增加了48.9 mm,但聲學(xué)包的重量平均降低為原始設(shè)計(jì)的62.2%。由此可見(jiàn),放開厚度約束,含空氣層聲學(xué)包的性能還有進(jìn)一步提升的潛力。與3.1節(jié)中數(shù)據(jù)進(jìn)行對(duì)比可以看出,含空氣層聲學(xué)包的隔聲量在上升6.5 dB的同時(shí)厚度增加了48.9 mm以上,重量變?yōu)樵荚O(shè)計(jì)的60%以上,也是控制厚度(3.1節(jié)優(yōu)化算例1)優(yōu)化后重量的兩倍,這說(shuō)明含有空氣層聲學(xué)包聲學(xué)性能的進(jìn)一步提升也需要以犧牲厚度和重量來(lái)達(dá)到。

    不同的是,相比于3.1節(jié)中的優(yōu)化結(jié)果,重量有所上升,但相比于原始設(shè)計(jì),其重量還是有大幅減小,因此當(dāng)主要任務(wù)為提升性能而可以適當(dāng)增加厚度時(shí),這種優(yōu)化約束也是適用的。

    4 結(jié) 論

    本文采用FE-SEA數(shù)值混響室法評(píng)價(jià)各聲學(xué)包設(shè)計(jì)的隔聲量;使用遺傳算法對(duì)聲學(xué)包各層材料厚度進(jìn)行優(yōu)化,在給定約束條件下獲得聲學(xué)包最優(yōu)的聲學(xué)性能。將算例聲學(xué)包分為含空氣層4種設(shè)計(jì)方案和不含空氣層4種設(shè)計(jì)方案進(jìn)行了性能分析及優(yōu)化。研究表明:

    1)當(dāng)控制優(yōu)化后聲學(xué)包的厚度不增加時(shí),不含空氣層聲學(xué)包的隔聲量約可提升1.4 dB,而含空氣層聲學(xué)包的隔聲量可提升20 dB以上,并且重量有大幅度的下降,說(shuō)明含空氣層聲學(xué)包的性能更好。

    2)當(dāng)設(shè)計(jì)條件允許加大聲學(xué)包厚度值時(shí),兩種類型聲學(xué)包的隔聲量均會(huì)進(jìn)一步增加,但其厚度和重量也會(huì)上升,因此使用時(shí)需要綜合考慮隔聲量、厚度、重量以及造價(jià)等一系列因素。

    3)FE-SEA數(shù)值混響室法可應(yīng)用于全頻段聲學(xué)包的優(yōu)化。本文給出的聲學(xué)包計(jì)算分析與優(yōu)化設(shè)計(jì)方法簡(jiǎn)單易行,對(duì)實(shí)際應(yīng)用具有一定的參考價(jià)值。

    [1]王世明.綠色船舶的現(xiàn)狀和發(fā)展前景分析[J].中國(guó)造船,2008,49(增刊 1):48-57.WANG S M.The analysis on present situation and de?velopment prospect for green ship[J].Shipbuilding of China,2008,49(Supp 1):48-57(in Chinese).

    [2]中國(guó)船級(jí)社.船舶及產(chǎn)品噪聲控制與檢測(cè)指南[M].北京:中國(guó)船級(jí)社,2013.China Classification Society.Guidance on noise control and testing of ships and products[M].Beijing:China Classification Society,2013(in Chinese).

    [3]康洪軍,蔣維鋼.聲學(xué)包仿真技術(shù)在高速綜合檢測(cè)列車降噪設(shè)計(jì)上的應(yīng)用與實(shí)踐[J].鐵道車輛,2012,50(4):17-21.KANG H J,JIANG W G.Application and practice of the sound package simulation technology in the noise reduction design of the high speed comprehensive in?spection trains[J].Rolling Stock,2012,50(4):17-21(in Chinese).

    [4]SAHA P.Future trends in the automotive sound pack?age industry[J].The Journal of the Acoustical Society of America,2001,109(5):2368.

    [5]范明偉.高速船艙室噪聲的統(tǒng)計(jì)能量分析與降噪優(yōu)化研究[D].武漢:武漢理工大學(xué),2012.FAN M W.Prediction and optimization of cabin noiseof high-speed ship based on statistical energy analysis method[D].Wuhan:Wuhan University of Technology,2012(in Chinese).

    [6]高處.船舶艙室聲學(xué)優(yōu)化設(shè)計(jì)理論與方法研究[D].上海:上海交通大學(xué),2015.GAO C.Optimization theory and methods of ship cab?ins acoustics design[D].Shanghai:Shanghai Jiao Tong University,2015(in Chinese).

    [7]蘇朝勇.基于VA One的聲結(jié)構(gòu)優(yōu)化系統(tǒng)的二次開發(fā)[D].大連:大連理工大學(xué),2009.SU Z Y.Secondary development of structural vibra?tion-acoustic optimization based on VA One[D].Da?lian:Dalian University of Technology,2009(in Chi?nese).

    [8]石嘉欣,楊德慶,郁揚(yáng).隔聲量計(jì)算的邊界元—有限元—統(tǒng)計(jì)能量法[J].噪聲與振動(dòng)控制,2016,36(6):62-66.SHI J X,YANG D Q,YU Y.BEM/FEM/SEA hybrid method for sound transmission loss calculation of acoustic materials[J].Noise and Vibration Control,2016,36(6):62-66(in Chinese).

    [9]COTONI V,SHORTER P,LANGLEY R.Numerical and experimental validation of a hybrid finite ele?ment-statistical energy analysis method[J].The Jour?nal of the Acoustical Society of America,2007,122(1):259-270.

    [10]LANGLEY R S,CORDIOLI J A.Hybrid determinis?tic-statistical analysis of vibro-acoustic systems with domain couplings on statistical components[J].Jour?nal of Sound and Vibration,2009,321(3/4/5):893-912.

    [11]應(yīng)申舜,盧奐采,姜偉.多孔橡膠材料聲傳遞損失性能分析:模型描述[J].輕工機(jī)械,2013,31(6):1-4.YING S S,LU H C,JIANG W.Sound transmission loss analysis of porous rubber material:descriptions of model[J].Light Industry Machinery,2013,31(6):1-4(in Chinese).

    [12]SGARD F C,ATALLA N,NICOLAS J.A numerical model for the low frequency diffuse field sound trans?mission loss of double-wall sound barriers with elas?tic porous linings[J].The Journal of the Acoustical Society of America,2000,108(6):2865-2872.

    Design and optimization of sound package using FE-SEA numerical reverberation method

    YANG Deqing1,2,3,SHI Jiaxin1,2,3,YU Yang1,2,3
    1 Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration,Shanghai 200240,China 2 State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University,Shanghai 200240,China 3 School of Naval Architecture,Ocean and Civil Engineering,Shanghai Jiao Tong University,Shanghai 200240,China

    To solve the low efficiency problem of the conventional sound package,the FE-SEA numerical reverberation method is applied in the design and optimization of sound packages in the whole frequency domain.Two types of sound package with and without air layers are designed and optimized.The FE-SEA numerical reverberation method is used to calculate sound transmission loss.The thicknesses of the air layers of the sound packages are defined as design variables.The optimization results reveal that the sound package with air layer has better acoustic performance than the one without air layers,and the sound transmission loss can be raised 2-6 dB after optimization.The research results have important guiding values for ship sound package design in actual engineering.

    sound package;optimization;FE-SEA numerical reverberation method;transmission loss

    U661.44

    A

    10.3969/j.issn.1673-3185.2017.04.005

    http://kns.cnki.net/kcms/detail/42.1755.TJ.20170727.1027.028.html期刊網(wǎng)址:www.ship-research.com

    楊德慶,石嘉欣,郁揚(yáng).基于有限元—統(tǒng)計(jì)能量數(shù)值混響室法的聲學(xué)包設(shè)計(jì)與優(yōu)化[J].中國(guó)艦船研究,2017,12(4):26-34.

    YANG D Q,SHI J X,YU Y.Design and optimization of sound package using FE-SEA numerical reverberation method[J].Chinese Journal of Ship Research,2017,12(4):26-34.

    2017-03-10< class="emphasis_bold">網(wǎng)絡(luò)出版時(shí)間:

    時(shí)間:2017-7-27 10:27

    國(guó)家高技術(shù)船舶科研計(jì)劃資助項(xiàng)目(2012-533);國(guó)家自然科學(xué)基金資助項(xiàng)目(51479115)

    楊德慶(通信作者),男,1968年生,博士,教授,博士生導(dǎo)師。研究方向:船舶結(jié)構(gòu)減振降噪理論與方法。E-mail:yangdq@sjtu.edu.cn

    猜你喜歡
    有限元優(yōu)化設(shè)計(jì)
    超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
    民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    瞞天過(guò)?!律O(shè)計(jì)萌到家
    設(shè)計(jì)秀
    海峽姐妹(2017年7期)2017-07-31 19:08:17
    有種設(shè)計(jì)叫而專
    Coco薇(2017年5期)2017-06-05 08:53:16
    磨削淬硬殘余應(yīng)力的有限元分析
    基于SolidWorks的吸嘴支撐臂有限元分析
    箱形孔軋制的有限元模擬
    上海金屬(2013年4期)2013-12-20 07:57:18
    欧美另类亚洲清纯唯美| 黄色成人免费大全| 午夜视频精品福利| 亚洲专区中文字幕在线| 黄片小视频在线播放| 国产伦人伦偷精品视频| 国产成人免费观看mmmm| 91九色精品人成在线观看| 一进一出抽搐gif免费好疼 | 亚洲精品粉嫩美女一区| 叶爱在线成人免费视频播放| 免费观看a级毛片全部| √禁漫天堂资源中文www| 又黄又粗又硬又大视频| 十八禁高潮呻吟视频| 天堂动漫精品| 超色免费av| 国产男女内射视频| 在线十欧美十亚洲十日本专区| 久久精品国产a三级三级三级| 啦啦啦 在线观看视频| 黄频高清免费视频| 女人久久www免费人成看片| 欧美黑人欧美精品刺激| 国产一区二区激情短视频| 欧美日韩中文字幕国产精品一区二区三区 | 成人国产一区最新在线观看| 婷婷精品国产亚洲av在线 | 性少妇av在线| 在线av久久热| 18禁黄网站禁片午夜丰满| 很黄的视频免费| 免费在线观看黄色视频的| 18在线观看网站| 午夜老司机福利片| 亚洲第一欧美日韩一区二区三区| 高清在线国产一区| 正在播放国产对白刺激| 欧美最黄视频在线播放免费 | 成人国产一区最新在线观看| 狠狠婷婷综合久久久久久88av| 精品国产一区二区久久| 91精品国产国语对白视频| 每晚都被弄得嗷嗷叫到高潮| 91成人精品电影| 1024香蕉在线观看| 如日韩欧美国产精品一区二区三区| 成人影院久久| 黑人巨大精品欧美一区二区mp4| 精品视频人人做人人爽| 99国产精品一区二区蜜桃av | 一区二区日韩欧美中文字幕| 高清毛片免费观看视频网站 | 国精品久久久久久国模美| 天堂俺去俺来也www色官网| 国产精品免费一区二区三区在线 | 热99国产精品久久久久久7| 国产亚洲精品第一综合不卡| 亚洲av日韩在线播放| 久久久国产成人精品二区 | 亚洲国产精品sss在线观看 | 精品熟女少妇八av免费久了| 亚洲欧美日韩高清在线视频| 嫩草影视91久久| 成人免费观看视频高清| 中出人妻视频一区二区| 黄色成人免费大全| 久久久水蜜桃国产精品网| 午夜精品在线福利| 黄色毛片三级朝国网站| 自线自在国产av| 久久久国产精品麻豆| 国产成人av教育| 国产1区2区3区精品| 久久久久久久国产电影| 国产精品自产拍在线观看55亚洲 | 久久精品aⅴ一区二区三区四区| 国产精品 国内视频| 欧洲精品卡2卡3卡4卡5卡区| 十八禁网站免费在线| 久久人妻福利社区极品人妻图片| 国产精品亚洲av一区麻豆| 不卡av一区二区三区| 在线观看一区二区三区激情| 狠狠狠狠99中文字幕| 美女扒开内裤让男人捅视频| 又紧又爽又黄一区二区| 五月开心婷婷网| 无人区码免费观看不卡| 在线观看舔阴道视频| 男女午夜视频在线观看| 久久亚洲真实| 中文字幕最新亚洲高清| 国产成人欧美| 亚洲成av片中文字幕在线观看| 如日韩欧美国产精品一区二区三区| 精品国产乱子伦一区二区三区| 激情视频va一区二区三区| 亚洲男人天堂网一区| 欧美成狂野欧美在线观看| 91老司机精品| 国产高清视频在线播放一区| 国内毛片毛片毛片毛片毛片| svipshipincom国产片| 麻豆成人av在线观看| 亚洲成人国产一区在线观看| 成年人黄色毛片网站| 亚洲精品久久午夜乱码| 日韩免费av在线播放| 女性生殖器流出的白浆| 欧美成人免费av一区二区三区 | 极品教师在线免费播放| 久久精品国产99精品国产亚洲性色 | 久久亚洲真实| 中文字幕av电影在线播放| 身体一侧抽搐| 久久久精品免费免费高清| 亚洲精品美女久久av网站| 成年人免费黄色播放视频| 亚洲国产精品一区二区三区在线| 麻豆av在线久日| 十分钟在线观看高清视频www| 久久人人爽av亚洲精品天堂| 国产在视频线精品| 建设人人有责人人尽责人人享有的| 国产野战对白在线观看| 在线免费观看的www视频| 亚洲人成77777在线视频| 中文欧美无线码| 大片电影免费在线观看免费| 中文字幕av电影在线播放| 在线观看66精品国产| 午夜免费成人在线视频| 亚洲,欧美精品.| 少妇裸体淫交视频免费看高清 | 黄色a级毛片大全视频| 免费人成视频x8x8入口观看| 久久这里只有精品19| 国产男女内射视频| 热re99久久国产66热| 女性生殖器流出的白浆| 麻豆av在线久日| 婷婷丁香在线五月| 在线十欧美十亚洲十日本专区| 免费看a级黄色片| 亚洲熟妇中文字幕五十中出 | 中文字幕av电影在线播放| 欧美精品亚洲一区二区| 亚洲精品国产区一区二| 亚洲第一av免费看| 久久人妻熟女aⅴ| 亚洲熟女精品中文字幕| 免费久久久久久久精品成人欧美视频| 精品少妇久久久久久888优播| 欧美不卡视频在线免费观看 | 欧美激情久久久久久爽电影 | 久久久久久久久免费视频了| av不卡在线播放| 日韩免费av在线播放| 国产精品免费一区二区三区在线 | 最近最新中文字幕大全电影3 | 新久久久久国产一级毛片| 大型av网站在线播放| 最近最新中文字幕大全电影3 | 国产欧美日韩综合在线一区二区| 色老头精品视频在线观看| 99re6热这里在线精品视频| 99久久国产精品久久久| 国产精品亚洲一级av第二区| 极品教师在线免费播放| 亚洲欧美日韩高清在线视频| 亚洲第一av免费看| 一进一出抽搐动态| 久热爱精品视频在线9| 久久久久精品人妻al黑| 大片电影免费在线观看免费| 少妇粗大呻吟视频| 日韩欧美免费精品| 国产aⅴ精品一区二区三区波| 午夜精品国产一区二区电影| 欧美精品人与动牲交sv欧美| 日日夜夜操网爽| 精品久久久久久,| 中国美女看黄片| 亚洲精品av麻豆狂野| 午夜福利视频在线观看免费| 久久九九热精品免费| 久久精品国产综合久久久| 少妇被粗大的猛进出69影院| 操美女的视频在线观看| 成人特级黄色片久久久久久久| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精华国产精华精| 国产精品一区二区免费欧美| 午夜免费成人在线视频| 国产精华一区二区三区| 色婷婷久久久亚洲欧美| 亚洲av电影在线进入| 男女高潮啪啪啪动态图| 极品少妇高潮喷水抽搐| 如日韩欧美国产精品一区二区三区| 一区二区三区国产精品乱码| 亚洲全国av大片| 99精国产麻豆久久婷婷| 久久草成人影院| 手机成人av网站| 国产欧美日韩一区二区三| 黄色片一级片一级黄色片| 精品福利观看| 黑丝袜美女国产一区| 夫妻午夜视频| 午夜影院日韩av| 18禁裸乳无遮挡免费网站照片 | 人妻 亚洲 视频| 大型av网站在线播放| 后天国语完整版免费观看| 午夜成年电影在线免费观看| 黄色毛片三级朝国网站| 超碰97精品在线观看| 久久久久精品人妻al黑| 热99re8久久精品国产| 大片电影免费在线观看免费| 国产成+人综合+亚洲专区| 欧美乱色亚洲激情| 婷婷精品国产亚洲av在线 | 成人18禁在线播放| 麻豆乱淫一区二区| 午夜免费观看网址| 91国产中文字幕| 国产在线精品亚洲第一网站| 飞空精品影院首页| 亚洲熟妇熟女久久| 久99久视频精品免费| av国产精品久久久久影院| 久久久久久久国产电影| 国产一区二区三区在线臀色熟女 | 国产精品欧美亚洲77777| 国产亚洲av高清不卡| www.999成人在线观看| 久久香蕉激情| 一区二区三区国产精品乱码| 波多野结衣av一区二区av| 高清毛片免费观看视频网站 | 动漫黄色视频在线观看| 亚洲av欧美aⅴ国产| 黄色毛片三级朝国网站| 十八禁网站免费在线| 美女扒开内裤让男人捅视频| 十八禁高潮呻吟视频| 两性夫妻黄色片| 又黄又爽又免费观看的视频| 国产男女超爽视频在线观看| 啦啦啦视频在线资源免费观看| 另类亚洲欧美激情| 亚洲第一欧美日韩一区二区三区| 极品少妇高潮喷水抽搐| 国产成人一区二区三区免费视频网站| 女人被躁到高潮嗷嗷叫费观| 欧美成狂野欧美在线观看| 激情在线观看视频在线高清 | 国产精品国产高清国产av | 一边摸一边抽搐一进一出视频| 亚洲精品粉嫩美女一区| 亚洲欧美日韩高清在线视频| 欧美午夜高清在线| 色尼玛亚洲综合影院| 色在线成人网| e午夜精品久久久久久久| 999精品在线视频| cao死你这个sao货| 日韩免费高清中文字幕av| 中文字幕人妻丝袜一区二区| 黄色丝袜av网址大全| 亚洲欧美色中文字幕在线| 后天国语完整版免费观看| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区三区在线臀色熟女 | 纯流量卡能插随身wifi吗| 一边摸一边抽搐一进一小说 | 十八禁人妻一区二区| www日本在线高清视频| 露出奶头的视频| 怎么达到女性高潮| 欧美老熟妇乱子伦牲交| 欧美精品av麻豆av| 精品亚洲成a人片在线观看| 一区在线观看完整版| 一级黄色大片毛片| 久久人妻熟女aⅴ| 在线观看免费高清a一片| 国产一卡二卡三卡精品| 久久ye,这里只有精品| 国产不卡一卡二| 最新在线观看一区二区三区| 国产男靠女视频免费网站| 99香蕉大伊视频| 亚洲三区欧美一区| 1024视频免费在线观看| 一边摸一边抽搐一进一出视频| 国内毛片毛片毛片毛片毛片| xxxhd国产人妻xxx| 国产91精品成人一区二区三区| 大码成人一级视频| 日本欧美视频一区| 国产单亲对白刺激| 亚洲欧美色中文字幕在线| 极品人妻少妇av视频| 女人被躁到高潮嗷嗷叫费观| 亚洲色图 男人天堂 中文字幕| 777久久人妻少妇嫩草av网站| 男女之事视频高清在线观看| 男男h啪啪无遮挡| 中文字幕人妻丝袜一区二区| 国产高清videossex| 丝瓜视频免费看黄片| 免费看十八禁软件| 大香蕉久久网| 国产精华一区二区三区| 女警被强在线播放| 老司机在亚洲福利影院| 亚洲视频免费观看视频| 国产精品 国内视频| 久久人妻av系列| 黑人巨大精品欧美一区二区蜜桃| 亚洲av欧美aⅴ国产| 免费一级毛片在线播放高清视频 | 自线自在国产av| 国产三级黄色录像| av网站在线播放免费| 精品福利观看| 久久久水蜜桃国产精品网| 可以免费在线观看a视频的电影网站| 黑人巨大精品欧美一区二区mp4| 国产成人精品在线电影| √禁漫天堂资源中文www| 999久久久国产精品视频| 男女免费视频国产| 免费观看精品视频网站| 亚洲熟妇熟女久久| 中出人妻视频一区二区| 一二三四在线观看免费中文在| 女人久久www免费人成看片| 熟女少妇亚洲综合色aaa.| 动漫黄色视频在线观看| 久久久久久久午夜电影 | 在线观看www视频免费| 日韩大码丰满熟妇| 日韩中文字幕欧美一区二区| 成年人黄色毛片网站| 亚洲中文字幕日韩| 亚洲色图 男人天堂 中文字幕| 精品一区二区三区视频在线观看免费 | 久久人人97超碰香蕉20202| 中文欧美无线码| 国产成+人综合+亚洲专区| 99热只有精品国产| 久久久久久久久免费视频了| 日本黄色日本黄色录像| 黄色怎么调成土黄色| 亚洲久久久国产精品| 丰满迷人的少妇在线观看| 日韩欧美一区视频在线观看| 免费av中文字幕在线| 两个人免费观看高清视频| 国产精品欧美亚洲77777| 国产高清激情床上av| 夫妻午夜视频| 黄片小视频在线播放| 国产精品永久免费网站| 在线永久观看黄色视频| 日韩欧美在线二视频 | 日日爽夜夜爽网站| 国产亚洲欧美在线一区二区| 久久国产亚洲av麻豆专区| 亚洲av日韩精品久久久久久密| 日日夜夜操网爽| 亚洲中文av在线| 黄片播放在线免费| 亚洲熟妇中文字幕五十中出 | 999久久久国产精品视频| 亚洲一卡2卡3卡4卡5卡精品中文| 老熟妇仑乱视频hdxx| 一级作爱视频免费观看| 91字幕亚洲| bbb黄色大片| 亚洲专区中文字幕在线| 精品无人区乱码1区二区| av天堂在线播放| 久久久久久亚洲精品国产蜜桃av| 亚洲国产毛片av蜜桃av| 一级片'在线观看视频| 亚洲第一av免费看| 丰满人妻熟妇乱又伦精品不卡| 国产精品国产av在线观看| 男女高潮啪啪啪动态图| 高清视频免费观看一区二区| 国产成人免费无遮挡视频| 美女高潮到喷水免费观看| 色综合婷婷激情| 久久中文字幕一级| 在线十欧美十亚洲十日本专区| 久久久久久久精品吃奶| 国产精品免费视频内射| 午夜精品在线福利| 老熟女久久久| 久久久国产成人精品二区 | 亚洲国产欧美一区二区综合| 国产视频一区二区在线看| 超色免费av| 老汉色∧v一级毛片| 女警被强在线播放| 在线视频色国产色| 国内久久婷婷六月综合欲色啪| 国产91精品成人一区二区三区| 黄色a级毛片大全视频| 国产午夜精品久久久久久| 看片在线看免费视频| 成在线人永久免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品偷伦视频观看了| 久久婷婷成人综合色麻豆| 免费观看a级毛片全部| 纯流量卡能插随身wifi吗| x7x7x7水蜜桃| 在线观看一区二区三区激情| 国产片内射在线| 国产精品免费一区二区三区在线 | 丰满饥渴人妻一区二区三| 人妻一区二区av| 久久久久久久国产电影| 欧美黑人精品巨大| 精品国产一区二区久久| 激情视频va一区二区三区| 亚洲五月婷婷丁香| 国产色视频综合| 午夜福利视频在线观看免费| 18禁裸乳无遮挡免费网站照片 | 精品少妇久久久久久888优播| 国产不卡av网站在线观看| 国产单亲对白刺激| 欧美精品高潮呻吟av久久| 午夜久久久在线观看| 黄片小视频在线播放| 日日夜夜操网爽| 久久国产精品大桥未久av| 老司机亚洲免费影院| 视频在线观看一区二区三区| 国产欧美日韩一区二区三区在线| 久久精品国产99精品国产亚洲性色 | 中文字幕最新亚洲高清| 久久香蕉国产精品| 欧美性长视频在线观看| 最近最新免费中文字幕在线| 欧美 日韩 精品 国产| 女人精品久久久久毛片| 久久精品国产亚洲av香蕉五月 | 午夜亚洲福利在线播放| 97人妻天天添夜夜摸| 欧美日韩成人在线一区二区| 午夜精品国产一区二区电影| 大码成人一级视频| 免费在线观看完整版高清| 99精品在免费线老司机午夜| 国产亚洲精品久久久久5区| 一级a爱视频在线免费观看| 成人黄色视频免费在线看| 精品人妻1区二区| 国产男女内射视频| 捣出白浆h1v1| 亚洲自偷自拍图片 自拍| 在线观看午夜福利视频| 久久人人97超碰香蕉20202| 亚洲第一欧美日韩一区二区三区| 欧美激情 高清一区二区三区| 国产激情欧美一区二区| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲av高清不卡| 日本一区二区免费在线视频| 欧美精品高潮呻吟av久久| 中文字幕最新亚洲高清| 人人妻人人澡人人爽人人夜夜| 久久久久国产精品人妻aⅴ院 | 涩涩av久久男人的天堂| 午夜久久久在线观看| 老汉色av国产亚洲站长工具| 欧美 亚洲 国产 日韩一| 日韩大码丰满熟妇| 精品人妻1区二区| 国产精品 欧美亚洲| 啦啦啦 在线观看视频| 超碰成人久久| 亚洲欧美色中文字幕在线| 久久久久精品国产欧美久久久| 老汉色av国产亚洲站长工具| 热re99久久国产66热| 国产成人精品在线电影| 亚洲av片天天在线观看| 天天添夜夜摸| 99久久99久久久精品蜜桃| 久久香蕉国产精品| 亚洲精品av麻豆狂野| 亚洲中文av在线| 热re99久久国产66热| 成人亚洲精品一区在线观看| 91精品国产国语对白视频| 久久久水蜜桃国产精品网| 99精品在免费线老司机午夜| 日韩三级视频一区二区三区| 窝窝影院91人妻| 午夜福利一区二区在线看| 国产欧美日韩精品亚洲av| 十分钟在线观看高清视频www| 精品人妻1区二区| 日韩一卡2卡3卡4卡2021年| 动漫黄色视频在线观看| 一边摸一边抽搐一进一小说 | 露出奶头的视频| 欧美日韩国产mv在线观看视频| 国产精品一区二区在线观看99| 国产1区2区3区精品| 午夜福利免费观看在线| 亚洲av片天天在线观看| 黄色视频,在线免费观看| 亚洲国产欧美网| 丝袜美腿诱惑在线| 国产精品久久久久成人av| 国产伦人伦偷精品视频| 在线十欧美十亚洲十日本专区| 久久久久久久精品吃奶| 午夜激情av网站| 最新的欧美精品一区二区| 午夜两性在线视频| 如日韩欧美国产精品一区二区三区| 老汉色∧v一级毛片| 在线观看免费午夜福利视频| 日韩欧美国产一区二区入口| 国产成人免费观看mmmm| 啦啦啦免费观看视频1| 午夜成年电影在线免费观看| 久久精品国产清高在天天线| 精品国产乱子伦一区二区三区| 精品免费久久久久久久清纯 | 国产亚洲欧美98| 欧美激情 高清一区二区三区| 午夜福利影视在线免费观看| 亚洲色图 男人天堂 中文字幕| 日本欧美视频一区| 中文字幕av电影在线播放| 啪啪无遮挡十八禁网站| 韩国精品一区二区三区| 中文字幕最新亚洲高清| 免费观看人在逋| 久久久精品区二区三区| 91麻豆av在线| av有码第一页| 麻豆成人av在线观看| 女人精品久久久久毛片| 中文字幕人妻丝袜一区二区| 又黄又爽又免费观看的视频| 99久久国产精品久久久| 1024香蕉在线观看| av不卡在线播放| 国产av又大| 国产精品香港三级国产av潘金莲| 中文字幕av电影在线播放| 最近最新免费中文字幕在线| 女人高潮潮喷娇喘18禁视频| 成人18禁高潮啪啪吃奶动态图| 国产亚洲精品久久久久5区| 久久婷婷成人综合色麻豆| 欧美日韩亚洲国产一区二区在线观看 | 国产精品成人在线| 啦啦啦免费观看视频1| 国产1区2区3区精品| 一边摸一边抽搐一进一小说 | 久久午夜综合久久蜜桃| 狠狠婷婷综合久久久久久88av| 免费在线观看影片大全网站| 午夜成年电影在线免费观看| 国产精品自产拍在线观看55亚洲 | 99久久精品国产亚洲精品| 丰满的人妻完整版| 91九色精品人成在线观看| 一区二区三区精品91| 天天影视国产精品| 男男h啪啪无遮挡| 99riav亚洲国产免费| 91成人精品电影| av国产精品久久久久影院| 久久天躁狠狠躁夜夜2o2o| 国产深夜福利视频在线观看| 亚洲情色 制服丝袜| 国产精品自产拍在线观看55亚洲 | 黄色视频,在线免费观看| 欧美乱妇无乱码| 国产无遮挡羞羞视频在线观看| 99香蕉大伊视频| 啦啦啦免费观看视频1| 久久热在线av| 久久久久精品人妻al黑| 脱女人内裤的视频| 99国产综合亚洲精品| 国产成+人综合+亚洲专区| 久久亚洲真实| 午夜两性在线视频| 久久精品亚洲av国产电影网| 欧美日韩黄片免| 国产成人一区二区三区免费视频网站| 亚洲国产精品sss在线观看 | 69av精品久久久久久| 超碰97精品在线观看| 日本wwww免费看| 久久狼人影院| 黄色怎么调成土黄色| 狠狠婷婷综合久久久久久88av|