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

    波紋芯體夾層板分頻段隔聲性能多參數(shù)優(yōu)化

    2016-10-14 07:54:04李會(huì)榮楊浩森米永振
    噪聲與振動(dòng)控制 2016年3期
    關(guān)鍵詞:參考模型芯體聲功率

    李會(huì)榮,鄭 輝,2,楊浩森,米永振

    波紋芯體夾層板分頻段隔聲性能多參數(shù)優(yōu)化

    李會(huì)榮1,鄭輝1,2,楊浩森1,米永振1

    (1.上海交通大學(xué) 振動(dòng)、沖擊、噪聲研究所,上海 200240;2.高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海 200240)

    應(yīng)用遺傳算法與梯度算法相結(jié)合的混合優(yōu)化算法,在結(jié)構(gòu)質(zhì)量及整體剛度約束下,分兩個(gè)頻段對(duì)波紋芯體夾層板在中置點(diǎn)聲源激勵(lì)下平均透射聲功率進(jìn)行單參數(shù)和多參數(shù)優(yōu)化分析,其中優(yōu)化目標(biāo)函數(shù),即夾層板透射聲功率用譜元法和Rayleigh積分進(jìn)行建模計(jì)算。優(yōu)化結(jié)果表明,單參數(shù)和多參數(shù)優(yōu)化都能明顯改善波紋芯體夾層板在高頻段(701 Hz~1 500 Hz)的隔聲性能,但兩種方法在低頻段(1 Hz~700 Hz)效果均不顯著,證明單純針對(duì)結(jié)構(gòu)參數(shù)優(yōu)化能夠有效地抑制由局部模態(tài)引起的透射聲峰值,但是對(duì)結(jié)構(gòu)整體模態(tài)引起的透射聲峰值影響有限。進(jìn)一步對(duì)比表明,對(duì)具有多個(gè)結(jié)構(gòu)參數(shù)的波紋芯體夾層板,其隔聲性能多參數(shù)聯(lián)合優(yōu)化結(jié)果明顯優(yōu)于單參數(shù)優(yōu)化結(jié)果。

    聲學(xué);夾層板;波紋芯體;聲透射;參數(shù)優(yōu)化;譜元法

    波紋芯體夾層板是上下兩層金屬板之間由波紋狀芯體板連接而成,因其具有剛度-質(zhì)量比(比剛度)和強(qiáng)度-質(zhì)量比(比強(qiáng)度)大的優(yōu)點(diǎn),被廣泛應(yīng)用于對(duì)輕量化要求較高的高速載運(yùn)工具。通常來(lái)說(shuō),高速運(yùn)載工具在行駛時(shí)隨著運(yùn)載工具運(yùn)行速度的提高,往往艙室噪聲會(huì)急劇增大。為保證乘客擁有良好的乘坐舒適性,要求作為壁板的波紋芯體夾層板具有良好的隔聲性能。關(guān)于夾層板隔聲性能的研究很多。在理論研究方面,Lee和Kim采用針對(duì)周期性?shī)A層板的空間諧波建模方法,在將筋板簡(jiǎn)化為剛體的基礎(chǔ)上研究了加筋板的聲透射問(wèn)題[1]。Xin等在研究三明治夾層板中的聲透射問(wèn)題采用波傳導(dǎo)模型,考慮了筋板剪切、彎曲和扭轉(zhuǎn)運(yùn)動(dòng),理論模型相對(duì)完整[2]。孫加平等采用傳遞矩陣法研究了三明治板的材料和結(jié)構(gòu)參數(shù)對(duì)高速列車底板隔聲性能的影響[3]。除了理論研究外,數(shù)值計(jì)算方法在對(duì)夾層板的振動(dòng)-聲學(xué)問(wèn)題的研究中也被廣泛應(yīng)用。黃勇和馬克儉根據(jù)空腹夾層板的構(gòu)造特點(diǎn),提出了基于板-塊體元模型的空腹夾層板有限元分析方法[4]。陳波等也使用有限元方法研究了工程中常見的夾層板結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)[5]。畢波等在考慮夾層板的剪切變形的同時(shí)計(jì)入了芯板的橫向變形,建立了等效單層模型的有限元?jiǎng)恿W(xué)方程[6]。

    夾層板結(jié)構(gòu)的隔聲性能與結(jié)構(gòu)的材料參數(shù)和幾何參數(shù)密切相關(guān),因此很多學(xué)者致力于研究夾層板結(jié)構(gòu)參數(shù)對(duì)其振動(dòng)聲學(xué)性能的影響,特別是針對(duì)夾層板結(jié)構(gòu)參數(shù)的優(yōu)化。Anthony等使用遺傳算法對(duì)一個(gè)二維加筋懸臂梁進(jìn)行了尺寸優(yōu)化,使傳遞到最外端梁的振動(dòng)能量最?。?]。Wang用高階夾層梁理論研究了在滿足機(jī)械剛度與隔聲性能前提下結(jié)構(gòu)的輕量化問(wèn)題[8]。李智等用有限元法對(duì)復(fù)合層板建立了數(shù)值模型,分析了板的結(jié)構(gòu)的材料和幾何參數(shù)對(duì)減振降噪的作用[9]。這些研究中基本都采用了有限元計(jì)算方法。然而,有限元方法建模過(guò)程復(fù)雜,變參數(shù)計(jì)算中結(jié)構(gòu)參數(shù)每次變化都需要重復(fù)建模。同時(shí),為保證較好的高頻計(jì)算精度,要求單元尺寸足夠小,導(dǎo)致高頻域隔聲性能優(yōu)化時(shí)間過(guò)長(zhǎng)。

    相對(duì)而言,譜元法(spectral element method,SEM)則為需要多次循環(huán)計(jì)算的優(yōu)化問(wèn)題提供了一種更好的選擇。該方法于1984年被提出,最初用于流體力學(xué)計(jì)算[10]。之后學(xué)者們又將譜元法應(yīng)用到結(jié)構(gòu)振動(dòng)方面的研究。Chakraborty和Gopalakrishnan用譜元法研究了夾層板中的波傳播問(wèn)題[11]。Lee U采用譜元法研究了兩梁之間接頭的動(dòng)力學(xué)特性,且計(jì)算結(jié)果與實(shí)驗(yàn)符合得較好[12]。吳志靜等建立了譜單元?jiǎng)恿W(xué)剛度矩陣,對(duì)鋼架結(jié)構(gòu)的動(dòng)力學(xué)響應(yīng)進(jìn)行了分析,證明了譜元法相較于傳統(tǒng)有限元法具有計(jì)算精度高的優(yōu)點(diǎn),特別是高頻計(jì)算精度高[13]。Wu等研究了周期性板結(jié)構(gòu)的振動(dòng)響應(yīng),證明了譜元法應(yīng)用于周期性結(jié)構(gòu)振動(dòng)建模計(jì)算的高效性[14]。譜元法的這些優(yōu)勢(shì)使其非常適合用于求解波紋芯體夾層板這類周期結(jié)構(gòu)隔聲優(yōu)化問(wèn)題時(shí)所需的循環(huán)計(jì)算。

    本文以波紋芯體夾層板為研究對(duì)象,采用遺傳算法和梯度算法組合構(gòu)成的優(yōu)化算法對(duì)中置點(diǎn)聲源激勵(lì)下構(gòu)件的頻段平均透射聲功率分別進(jìn)行了單參數(shù)和多參數(shù)優(yōu)化。優(yōu)化以結(jié)構(gòu)質(zhì)量及整體剛度為約束條件、分低頻和高頻兩個(gè)頻段進(jìn)行,目標(biāo)為頻段平均透射聲功率最小。其中優(yōu)化頻段的選取基于波紋芯體夾層板參考模型的振動(dòng)模態(tài)分布,低頻段中所有結(jié)構(gòu)模態(tài)均為整體模態(tài),高頻段則包含部分局部模態(tài);而優(yōu)化目標(biāo)函數(shù),即結(jié)構(gòu)在聲激勵(lì)下的輻射聲功率,先運(yùn)用譜元法計(jì)算得到結(jié)構(gòu)振動(dòng)響應(yīng)、然后通過(guò)Rayleigh積分計(jì)算得到。最后對(duì)兩個(gè)不同頻段的多參數(shù)和單參數(shù)優(yōu)化結(jié)果進(jìn)行對(duì)比分析討論。

    2 透射聲功率計(jì)算模型

    2.1結(jié)構(gòu)及聲透射計(jì)算理論模型

    典型波紋芯體夾層板由上板、波紋芯體板和下板組成,其三維結(jié)構(gòu)如圖1所示。

    圖1 波紋芯體夾層板三維結(jié)構(gòu)

    圖2給出夾層板的二維截面圖。夾層板總長(zhǎng)為L(zhǎng),空腔高度(上下面板間距)為H,上下面面板厚度分別為h1、h3,中間波紋芯體板厚為h2,筋板傾角為α。

    因?yàn)閵A層板平面內(nèi)沿波紋垂直方向的橫向彎曲剛度遠(yuǎn)大于波紋方向,通常假設(shè)夾層板沿波紋垂直方向無(wú)限延伸,使得三維夾層板的隔聲問(wèn)題簡(jiǎn)化為圖2所示的二維結(jié)構(gòu)(夾層梁)的隔聲問(wèn)題,其聲透射計(jì)算模型如圖3所示。

    圖2 波紋芯體夾層板二維尺寸圖

    圖3 波紋芯體夾層梁聲透射示意圖

    在圖3所示的聲透射計(jì)算模型中,波紋芯體夾層梁兩端簡(jiǎn)支并置于半無(wú)限長(zhǎng)聲障板之中,下簡(jiǎn)支梁受到點(diǎn)聲源的激勵(lì)而振動(dòng),其振動(dòng)能量通過(guò)芯梁結(jié)構(gòu)及梁內(nèi)聲腔傳至上梁,引起上梁振動(dòng)并向外輻射聲波。考慮到波紋芯體夾層梁剛度較空氣剛度大很多,可將空腔的聲傳遞路徑忽略[15],僅考慮結(jié)構(gòu)路徑傳聲。

    2.2結(jié)構(gòu)振動(dòng)響應(yīng)的譜元法計(jì)算

    基于圖3所示聲透射模型,采用譜元法(Spectral element method,SEM)計(jì)算結(jié)構(gòu)在聲激勵(lì)下的振動(dòng)響應(yīng)。譜元法結(jié)合了傳統(tǒng)有限元方法和譜方法的優(yōu)點(diǎn),與傳統(tǒng)有限元方法的區(qū)別在于此方法基于動(dòng)態(tài)剛度矩陣求解[16]。李蘭清等關(guān)于加筋雙層板聲透射的研究中通過(guò)有限元方法驗(yàn)證了譜元法計(jì)算夾層板結(jié)構(gòu)的精度[17],該文中給出了譜元法計(jì)算夾層板在入射簡(jiǎn)諧聲波激勵(lì)下振動(dòng)響應(yīng)的詳細(xì)數(shù)學(xué)推導(dǎo),此處不再贅述。圖4給出了譜元法及傳統(tǒng)有限元計(jì)算波紋芯體夾層梁在下梁受到1 Hz~2 000 Hz平面波聲激勵(lì)作用下上梁振動(dòng)能量的結(jié)果比較。

    圖4 上梁振動(dòng)能量譜元法計(jì)算結(jié)果與有限元對(duì)比

    從圖4可以看出,譜元法和有限元法得到的振動(dòng)能量-頻率曲線在整個(gè)計(jì)算頻段內(nèi)吻合很好。更重要的是,譜元法計(jì)算耗時(shí)不超過(guò)有限元法計(jì)算時(shí)間的5%。這對(duì)于下文將要討論的波紋芯體結(jié)構(gòu)的隔聲優(yōu)化具有重要意義。

    2.3聲輻射計(jì)算

    應(yīng)用譜元法計(jì)算得到結(jié)構(gòu)的振動(dòng)響應(yīng)后,結(jié)合Rayleigh積分可以計(jì)算結(jié)構(gòu)向外的輻射聲場(chǎng)[18]。根據(jù)Rayleigh積分,圖3中半徑為R的半圓面上任意一點(diǎn)由結(jié)構(gòu)上梁振動(dòng)所引起的聲壓為

    在整個(gè)半圓面S上對(duì)聲強(qiáng)積分即可得到波紋芯體夾層梁通過(guò)上梁向外輻射的聲功率,即結(jié)構(gòu)的透射聲功率[19]

    3 隔聲性能的優(yōu)化分析

    在進(jìn)行波紋芯體夾層梁的隔聲性能優(yōu)化之前,需要首先確定參考模型,并根據(jù)參考模型的頻響特性選擇優(yōu)化目標(biāo)頻段以及優(yōu)化約束條件。

    3.1優(yōu)化頻段選擇

    參考模型截面如圖2所示,芯體為6周期波紋梁(12筋梁),上、下梁及芯體的材料都相同,結(jié)構(gòu)的材料參數(shù)和幾何參數(shù)如表1所列。對(duì)參考模型進(jìn)行模態(tài)分析,得出參考模型橫向振動(dòng)模態(tài)在1 Hz~700 Hz頻段全部為整體彎曲模態(tài),頻率超過(guò)700 Hz后局部模態(tài)才開始出現(xiàn),701 Hz~1 500 Hz頻段內(nèi)整體模態(tài)占比約為50%。因此,為了研究結(jié)構(gòu)模態(tài)對(duì)隔聲性能優(yōu)化的影響,將優(yōu)化的目標(biāo)頻段分為1 Hz~700 Hz(低頻段)和701 Hz~1 500 Hz(高頻段)。

    表1 參考模型的材料及幾何參數(shù)

    3.2優(yōu)化問(wèn)題的數(shù)學(xué)表述

    3.2.1單參數(shù)優(yōu)化問(wèn)題表述

    相對(duì)于隔聲領(lǐng)域廣泛研究的單層板和雙層板結(jié)構(gòu),波紋芯體夾層梁影響隔聲性能優(yōu)化目標(biāo)函數(shù)的結(jié)構(gòu)參數(shù)較多,若進(jìn)行多參數(shù)優(yōu)化,容易產(chǎn)生計(jì)算量大、優(yōu)化效率低的問(wèn)題,對(duì)于需要考慮時(shí)間成本的情況,單參數(shù)優(yōu)化是較好的選擇,因此首先進(jìn)行單參數(shù)優(yōu)化。

    以入射聲波頻段內(nèi)平均透射聲功率來(lái)評(píng)價(jià)構(gòu)件在該頻段的隔聲性能,將目標(biāo)頻段內(nèi)的平均透射聲功率最小化作為優(yōu)化目標(biāo),分別進(jìn)行兩個(gè)目標(biāo)頻段的優(yōu)化:低頻段(1 Hz~700 Hz)和高頻段(701 Hz~1 500 Hz)。目標(biāo)函數(shù)計(jì)算式為

    其中ωh、ωl分別為頻段的圓周頻率上下限。

    考慮到對(duì)結(jié)構(gòu)有輕量化要求,限定其總質(zhì)量不得超過(guò)參考模型的質(zhì)量,即其截面面積不得超過(guò)參考模型的截面面積。另外,為了保證夾層梁具有一定的整體剛度,其1階固有頻率不可過(guò)低[20],因此要求第1階固有頻率不得低于參考模型的1階固有頻率值。

    波紋芯體夾層梁結(jié)構(gòu)參數(shù)中影響透射聲功率的主要參數(shù)有上下梁厚度h1、h3,芯體厚度h2以及芯體傾角α[21]。在單參數(shù)優(yōu)化中,僅保留一個(gè)設(shè)計(jì)參數(shù)作為優(yōu)化變量,并通過(guò)附加面積不變約束條件建立與其他設(shè)計(jì)變量的關(guān)系,從而將其引入優(yōu)化問(wèn)題中。據(jù)此,設(shè)計(jì)了兩個(gè)單參數(shù)優(yōu)化組,優(yōu)化組一(Optimization I)和優(yōu)化組二(Optimization II),其中優(yōu)化組一考察了上梁厚度的優(yōu)化效果,而優(yōu)化組二是考察了芯體傾角的優(yōu)化效果。單參數(shù)優(yōu)化兩個(gè)優(yōu)化組的優(yōu)化變量取值范圍以及約束條件如表2所列。

    表2 單參數(shù)優(yōu)化變量及約束

    3.2.2多參數(shù)優(yōu)化問(wèn)題表述

    考慮到若同時(shí)優(yōu)化對(duì)目標(biāo)函數(shù)有影響的四個(gè)主要的結(jié)構(gòu)參數(shù),則有可能得到更好的優(yōu)化效果,進(jìn)一步進(jìn)行多參數(shù)的聯(lián)合優(yōu)化。多參數(shù)聯(lián)合優(yōu)化的目標(biāo)函數(shù)和單參數(shù)優(yōu)化相同。優(yōu)化變量及其取值范圍為

    約束條件為截面積約束和1階固有頻率約束

    其中S為夾層梁截面積,S0為參考模型夾層梁截面積;f為11階固有頻率,f0為參考模型的1階固有頻率。

    3.3優(yōu)化方法

    對(duì)上述優(yōu)化問(wèn)題使用遺傳算法和梯度算法的組合優(yōu)化算法進(jìn)行優(yōu)化計(jì)算。遺傳算法因計(jì)算的隨機(jī)性使其具有良好的全局搜索能力,然而后期計(jì)算效率低。而梯度算法則能快速地找到局部最優(yōu)解。因此利用遺傳算法的全局探索能力找到全局優(yōu)化解,然后將此全局優(yōu)化解作為梯度優(yōu)化算法的初始值,在其周圍進(jìn)行局部搜索,尋找更加精確的最優(yōu)解。

    為了保證優(yōu)化解的可靠性,遺傳算法每次優(yōu)化至少60代,每代個(gè)體數(shù)為500,當(dāng)連續(xù)4代最優(yōu)結(jié)果相對(duì)誤差小于0.1%時(shí),即認(rèn)為優(yōu)化收斂。

    3.4優(yōu)化結(jié)果及討論

    3.4.1單參數(shù)優(yōu)化

    低頻段為1 Hz~700 Hz,在此頻段內(nèi)分別根據(jù)優(yōu)化組一和優(yōu)化組二對(duì)夾層梁的透射聲功率進(jìn)行優(yōu)化,優(yōu)化結(jié)果以及優(yōu)化模型與參考模型性能指標(biāo)參數(shù)對(duì)比如表3所示。這里僅給出了優(yōu)化變量的優(yōu)化結(jié)果,其余設(shè)計(jì)變量可通過(guò)約束條件求出,或與參考模型保持一致。

    表3 低頻段單參數(shù)優(yōu)化結(jié)果

    由表3可知,兩個(gè)優(yōu)化組透射聲功率相對(duì)于參考模型均降低了1.1 dB,降幅較小。將優(yōu)化模型與參考模型在優(yōu)化頻段內(nèi)的透射聲功率譜進(jìn)行對(duì)比,圖5給出了兩組優(yōu)化結(jié)果與參考模型的透射聲功率頻譜曲線對(duì)比。從圖中可以看出,兩個(gè)單參數(shù)優(yōu)化組在低頻段內(nèi)對(duì)透射聲功率優(yōu)化效果表現(xiàn)在各峰值有所降低,然而峰值個(gè)數(shù)并沒有減少,優(yōu)化效果不明顯。

    高頻段為701 Hz~1 500 Hz,同樣分別對(duì)兩個(gè)優(yōu)化組進(jìn)行單參數(shù)優(yōu)化,優(yōu)化結(jié)果如表4所列。優(yōu)化前后的頻譜特性曲線如圖6所示。

    表4 高頻段單參數(shù)優(yōu)化結(jié)果

    由表4可以看出在高頻段下,相對(duì)于參考模型,兩個(gè)優(yōu)化組的優(yōu)化模型平均透射聲功率大幅降低,其中優(yōu)化組一降低了31.3 dB;優(yōu)化組二降低了29.9 dB。

    從圖6可以看出,高頻段下兩組單參數(shù)優(yōu)化相對(duì)于參考模型均得到了較好的優(yōu)化效果:頻段內(nèi)平均透射聲功率大幅降低;頻段內(nèi)峰值個(gè)數(shù)也有減少。

    圖5 低頻段下兩組單參數(shù)優(yōu)化結(jié)果

    圖6 高頻段下兩組單參數(shù)優(yōu)化結(jié)果

    3.4.2多參數(shù)優(yōu)化

    多參數(shù)優(yōu)化仍然分別在低頻段(1 Hz~700 Hz)和高頻段(701 Hz~1500 Hz)進(jìn)行,優(yōu)化后模型的結(jié)構(gòu)參數(shù)如表5所示。表6給出了優(yōu)化后模型的1階固有頻率、對(duì)應(yīng)頻段的平均透射聲功率及截面相對(duì)面積(相對(duì)于參考模型)。

    表5 優(yōu)化后模型的結(jié)構(gòu)參數(shù)

    表6 優(yōu)化后模型的性能指標(biāo)

    圖7和圖8分別給出了兩個(gè)頻段優(yōu)化前后模型的透射聲功率曲線。從圖7所示的透射聲功率頻譜曲線比較可以看出,對(duì)于低頻段,優(yōu)化后的結(jié)構(gòu)相對(duì)于參考模型峰值個(gè)數(shù)并沒有減少,只是使一些峰值有所降低。圖8中結(jié)果表明,在高頻段,優(yōu)化后的結(jié)構(gòu)隔聲性能明顯優(yōu)于參考模型,優(yōu)化后模型不僅擁有更低的透射聲功率峰值,而且透射峰數(shù)量相對(duì)于參考模型也大量減少。

    經(jīng)過(guò)多參數(shù)聯(lián)合優(yōu)化,相對(duì)參考模型,低頻段平均透射聲功率降低了2.6 dB,高頻段平均透射聲功率降低了35.6 dB。

    3.4.3單參數(shù)優(yōu)化與多參數(shù)優(yōu)化結(jié)果比較

    表7列出了分別使用單/多參數(shù)優(yōu)化得到的優(yōu)化結(jié)果。可以看出,對(duì)低頻段而言,不管采用單參數(shù)優(yōu)化還是多參數(shù)優(yōu)化方法,透射聲功率降低效果都不明顯;反之,對(duì)高頻段,兩種優(yōu)化方法都得到了很好的優(yōu)化效果。

    由優(yōu)化目標(biāo)頻段選取可知,參考模型在低頻段內(nèi)結(jié)構(gòu)的模態(tài)全部為整體模態(tài),而高頻段中約50%為局部模態(tài)。本文對(duì)波紋芯體夾層梁結(jié)構(gòu)的優(yōu)化是幾何尺寸的優(yōu)化,而不改變模型的拓?fù)錁?gòu)型。單純進(jìn)行尺寸優(yōu)化對(duì)結(jié)構(gòu)的整體模態(tài)影響有限,因而對(duì)低頻段的優(yōu)化效果不佳;幾何尺寸優(yōu)化會(huì)改變局部模態(tài)從而使因局部模態(tài)引起的共振峰值有效降低,故對(duì)于有較多局部模態(tài)的高頻段優(yōu)化效果顯著。

    圖7 低頻段透射聲功率優(yōu)化結(jié)果

    圖8 高頻段透射聲功率優(yōu)化結(jié)果

    表7 單參數(shù)與多參數(shù)優(yōu)化結(jié)果比較

    為進(jìn)一步證明這一結(jié)論,本文選取了更高的頻率段,即1 000 Hz~2 000 Hz,對(duì)參考模型在此進(jìn)行了多參數(shù)優(yōu)化。參考模型的模態(tài)分析結(jié)果顯示,此頻段內(nèi)局部模態(tài)占比為約70%。優(yōu)化結(jié)果如圖9所示。從圖中可以看出,優(yōu)化模型在優(yōu)化頻段內(nèi)相對(duì)于參考模型透射聲功率大幅度降低,透射聲功率峰值數(shù)量減少很多。

    圖9 1 000 Hz~2 000 Hz優(yōu)化結(jié)果

    將單參數(shù)優(yōu)化與多參數(shù)優(yōu)化進(jìn)行同頻段優(yōu)化結(jié)果對(duì)比,發(fā)現(xiàn)多參數(shù)優(yōu)化不但使透射聲功率降低更多,并且使結(jié)構(gòu)剛度增大、結(jié)構(gòu)變輕。因此多參數(shù)優(yōu)化效果明顯優(yōu)于單參數(shù)優(yōu)化。

    4 結(jié)語(yǔ)

    本文針對(duì)波紋芯體夾層梁的隔聲性能分別進(jìn)行了單參數(shù)和多參數(shù)優(yōu)化。優(yōu)化結(jié)果表明,針對(duì)高頻段(700 Hz~1 500 Hz)的隔聲性能優(yōu)化結(jié)果明顯優(yōu)于低頻段(1 Hz~700 Hz),這是由于隨著頻率的增加,結(jié)構(gòu)局部模態(tài)增多,改變夾層梁芯體結(jié)構(gòu)參數(shù)能夠有效地改變局部模態(tài),但是對(duì)整個(gè)結(jié)構(gòu)的整體模態(tài)影響有限。

    通過(guò)對(duì)比多參數(shù)優(yōu)化與單參數(shù)優(yōu)化的結(jié)果說(shuō)明,對(duì)于結(jié)構(gòu)參數(shù)較多的復(fù)雜構(gòu)型,多參數(shù)優(yōu)化能夠更大限度地提高整體結(jié)構(gòu)的隔聲性能。雖然如此,考慮到結(jié)構(gòu)在低頻段隔聲性能通過(guò)單純結(jié)構(gòu)參數(shù)優(yōu)化難以得到很好的效果,若要進(jìn)一步提高結(jié)構(gòu)的隔聲性能,應(yīng)從其他途徑考慮,如阻尼處理、主動(dòng)控制等。

    [1]LEE J H,KIM J.Analysis of sound transmission through periodically stiffened panels by space-harmonic expansion method[J].Journal of Sound and Vibration,2002,251(2):349-366.

    [2]XIN F X,LU T J.Analytical modeling of wave propagationinorthogonallyrib-stiffenedsandwich structures:Sound radiation[J].Computers&Structures,2011,89(5):507-516.

    [3]孫加平,張麗榮,孫海榮,等.高速列車夾芯地板結(jié)構(gòu)隔聲特性研究[J].噪聲與振動(dòng)控制,2014,34(4):39-43.

    [4]黃 勇,馬克儉.基于板-塊體模型的空腹夾層板有限元分析[J].貴州工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2002(31):10-14.

    [5]陳波,黃勇,馬克儉.空腹夾層板有限元分析的兩種模型[J].建筑結(jié)構(gòu),2001:58-60.

    [6]畢波,徐干成.一種夾層板結(jié)構(gòu)的動(dòng)力有限元方法[J].空軍工程大學(xué)學(xué)報(bào)(自然科學(xué)版),2002,3(3):77-80.

    [7]ANTHONY D K,ELLIOTT S J,KEANE A J.Robustness ofoptimaldesignsolutionstoreducevibrationtransmission in a lightweight 2-D structure,part I:Geometric design[J].Journal of Sound and Vibration,2000,229(3):505-528.

    [8]WANG T,LI S,NUTT S R.Optimal design of acoustical sandwich panels with a genetic algorithm[J].Applied Acoustics,2009,70(3):416-425.

    [9]李智,殷祥超,何興華,等.復(fù)合層板減振降噪特性的數(shù)值模擬研究[J].噪聲與振動(dòng)控制,2006,26(6):27-30.

    [10]PATERA A T.A spectral element method for fluid dynamics:laminar flow in a channel expansion[J]. Journal of Computational Physics,1984,54(3):468-488.

    [11]CHAKRABORTY A,GOPALAKRISHNAN S.A spectral finite element model for wave propagation analysis in laminated composite plate[J].Journal of Vibration and Acoustics,2006,128(4):477-488.

    [12]LEE U.Dynamic characterization of the joints in a beam structure by using spectral element method[J].Shock and Vibration,2001,8(6):357-366.

    [13]吳志靜,李鳳明,胡恒山,等.譜元法在求解剛架結(jié)構(gòu)動(dòng)力學(xué)問(wèn)題中的應(yīng)用[J].動(dòng)力學(xué)與控制學(xué)報(bào),2012,10(1):71-75.

    [14]WU Z J,LI F M,WANG Y Z.Study on vibration characteristics in periodic plate structures using the spectral element method[J].Acta Mechanica,2013,224(5):1089-1101.

    [15]CHENG L,LI Y Y,GAO J X.Energy transmission in a mechanically-linked double-wall structure coupled to an acoustic enclosure[J].The Journal of the Acoustical Society ofAmerica,2005,117(5):2742-2751.

    [16]LEE U.Spectral element method in structural dynamics [M].John Wiley&Sons,2009.

    [17]李蘭清,鄭輝.基于譜元法的加筋雙層板聲透射分析[J].噪聲與振動(dòng)控制,2015,35(3):181-185.

    [18]FAHY F J,GARDONIO P.Sound and structural vibration:radiation,transmissionandresponse[M]. Academic press,2007.

    [19]杜功煥,朱哲民.聲學(xué)基礎(chǔ)[M].南京:南京大學(xué)出版社,2001.

    [20]BERRY A,NICOLASJ.Structuralacousticsand vibrationbehaviorofcomplexpanels[J].Applied Acoustics,1994,43(3):185-215.

    [21]SHEN C,XIN F X,LU T J.Theoretical model for sound transmission through finite sandwich structures with corrugatedcore[J].InternationalJournalofNon-Linear Mechanics,2012,47(10):1066-1072.

    Multi-parameter Optimization of the Sound Insulation Performance of a Corrugated Core Sandwich Plate in Two Different Frequency Ranges

    LI Hui-rong1,ZHENGHui1,2,YANG Hao-sen1,MI Yong-zhen1
    (1.Institute of Vibration,Shock and Noise,Shanghai Jiaotong University,Shanghai 200240,China;2.Collaborative Innovation Center forAdvanced Ship and Deep-sea Exploration(CISSE),Shanghai 200240,China)

    A hybrid optimization combining genetic algorithm with gradient-based algorithm is applied to minimize the average radiation sound power of a corrugated-core sandwich plate under point sound source excitation in two specific frequency bands with its structural mass and stiffness constrained.With the sound transmission power of the sandwich plate as the target function of optimization,the spectral element method and Rayleigh's integral are applied to the model establishment and analysis.The results show that the optimization results of both single-variable and multi-variable optimizations can significantly improve the isolation performance of the corrugated core sandwich plate in the high frequency band(701 Hz-1500 Hz).But both optimization methods are less effective in low frequency band(1 Hz-700 Hz).It implies that the corrugated core parameters optimization can effectively suppress the sound radiation peaks caused by local modes,but has very little suppression effect on the peaks caused by global modes.Further comparison of the results obtained through single-variable optimization and multi-variable optimization shows that for the structure of corrugated core sandwich plate with multiple parameters,the multi-variable optimization is much better than the single-variable optimization in soundinsulating performance.

    acoustics;sandwich plate;corrugated core;sound transmission;parameter optimization;spectral element method

    O327 O42

    ADOI編碼:10.3969/j.issn.1006-1335.2016.03.039

    1006-1355(2016)03-0190-07

    2015-11-16

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51275285)

    李會(huì)榮(1990-),女,湖北省漢川市人,碩士研究生,主要研究方向?yàn)樵肼暸c振動(dòng)控制。

    鄭輝,男,教授,博士生導(dǎo)師。E-mail:huizheng@sjtu.edu.cn

    猜你喜歡
    參考模型芯體聲功率
    復(fù)合芯體技術(shù)進(jìn)展概況
    生活用紙(2021年1期)2021-01-09 10:27:30
    核燃料環(huán)形芯體壓燒模具的設(shè)計(jì)
    專題:吸收性衛(wèi)生用品芯體結(jié)構(gòu)及SAP選擇研究
    生活用紙(2019年3期)2019-02-13 04:26:18
    整體道床軌道扣件剛度對(duì)鋼軌聲功率特性的影響
    不含絨毛漿超薄芯體新技術(shù)
    生活用紙(2017年9期)2017-09-19 12:43:13
    基于環(huán)境的軍事信息系統(tǒng)需求參考模型
    語(yǔ)義網(wǎng)絡(luò)P2P參考模型的查詢過(guò)程構(gòu)建
    自由風(fēng)扇聲功率級(jí)測(cè)量方法與測(cè)量不確定度
    一種新的壓縮機(jī)噪聲聲功率的室內(nèi)測(cè)量方法
    基于參考模型的制造業(yè)MES需求建模
    国产三级在线视频| 精品午夜福利在线看| 一个人免费在线观看电影| 精品人妻熟女av久视频| 天堂网av新在线| 黄色一级大片看看| 国内揄拍国产精品人妻在线| 校园人妻丝袜中文字幕| 亚洲一区高清亚洲精品| 日韩在线高清观看一区二区三区| 国产高潮美女av| 高清日韩中文字幕在线| 亚洲精品国产av成人精品| 成年av动漫网址| 欧美97在线视频| 白带黄色成豆腐渣| 亚洲怡红院男人天堂| 国产黄片视频在线免费观看| 亚洲欧美精品综合久久99| 国产亚洲精品av在线| 亚洲精品自拍成人| 一级二级三级毛片免费看| 国产老妇伦熟女老妇高清| 中文在线观看免费www的网站| 看十八女毛片水多多多| 精品国产三级普通话版| 亚洲国产最新在线播放| 亚洲人成网站高清观看| 国产极品精品免费视频能看的| 精品一区二区免费观看| 一本久久精品| 在线播放国产精品三级| 国产老妇女一区| 欧美日韩精品成人综合77777| 国产精品人妻久久久影院| 97超视频在线观看视频| 亚洲欧美成人精品一区二区| 久久精品国产亚洲av涩爱| 熟女电影av网| 国产av码专区亚洲av| 一卡2卡三卡四卡精品乱码亚洲| 嘟嘟电影网在线观看| 在线免费观看不下载黄p国产| 日韩av在线免费看完整版不卡| 成人毛片60女人毛片免费| 黄片wwwwww| 国产成年人精品一区二区| 搡女人真爽免费视频火全软件| 久久99蜜桃精品久久| 久久亚洲精品不卡| 亚洲av成人av| 精品免费久久久久久久清纯| 99久国产av精品国产电影| 在线播放国产精品三级| 美女被艹到高潮喷水动态| 欧美变态另类bdsm刘玥| 看免费成人av毛片| 中文天堂在线官网| 国内揄拍国产精品人妻在线| 亚洲18禁久久av| 91久久精品电影网| 国产一级毛片七仙女欲春2| 国产亚洲5aaaaa淫片| 大话2 男鬼变身卡| 草草在线视频免费看| 我要看日韩黄色一级片| 久久久久久久午夜电影| 久久精品综合一区二区三区| 精品久久国产蜜桃| 男女国产视频网站| 69av精品久久久久久| 青春草亚洲视频在线观看| 免费大片18禁| 日本爱情动作片www.在线观看| 美女xxoo啪啪120秒动态图| 亚洲精品国产成人久久av| 长腿黑丝高跟| 床上黄色一级片| 在线观看美女被高潮喷水网站| 男女边吃奶边做爰视频| 免费搜索国产男女视频| 国产美女午夜福利| 亚洲人成网站在线播| 中文资源天堂在线| 欧美最新免费一区二区三区| 精品人妻偷拍中文字幕| 亚洲国产精品sss在线观看| 热99在线观看视频| 亚洲熟妇中文字幕五十中出| 91在线精品国自产拍蜜月| 国产人妻一区二区三区在| 亚洲精品乱码久久久v下载方式| 波多野结衣巨乳人妻| 免费av不卡在线播放| 黄色一级大片看看| 熟女电影av网| 97人妻精品一区二区三区麻豆| 最近视频中文字幕2019在线8| 91精品伊人久久大香线蕉| 99热6这里只有精品| 成人欧美大片| 白带黄色成豆腐渣| 日韩制服骚丝袜av| 综合色av麻豆| 在线免费观看不下载黄p国产| 直男gayav资源| 麻豆成人av视频| 国产精品综合久久久久久久免费| 久久国产乱子免费精品| 国产一区二区在线观看日韩| 看非洲黑人一级黄片| 亚洲精品成人久久久久久| 国产一区二区在线观看日韩| 国产私拍福利视频在线观看| 欧美zozozo另类| 三级国产精品片| 赤兔流量卡办理| 亚洲电影在线观看av| 久久久久久久久大av| 欧美成人一区二区免费高清观看| 欧美成人精品欧美一级黄| 亚洲最大成人中文| 日韩精品青青久久久久久| av免费在线看不卡| 99国产精品一区二区蜜桃av| 狂野欧美白嫩少妇大欣赏| 少妇的逼水好多| 麻豆久久精品国产亚洲av| 亚洲成人av在线免费| 成人毛片a级毛片在线播放| 菩萨蛮人人尽说江南好唐韦庄 | kizo精华| 不卡视频在线观看欧美| 久久久久久久久久成人| 国产白丝娇喘喷水9色精品| 国产老妇女一区| 亚洲精品,欧美精品| 26uuu在线亚洲综合色| 九九热线精品视视频播放| 亚洲精品国产av成人精品| 亚洲不卡免费看| 不卡视频在线观看欧美| 91aial.com中文字幕在线观看| АⅤ资源中文在线天堂| 欧美激情在线99| 91久久精品国产一区二区成人| 亚洲精品影视一区二区三区av| 日韩高清综合在线| 最新中文字幕久久久久| 黑人高潮一二区| 男女边吃奶边做爰视频| 99久国产av精品国产电影| 国产av一区在线观看免费| 九草在线视频观看| 国产午夜精品久久久久久一区二区三区| 有码 亚洲区| 国产亚洲精品久久久com| 99久久人妻综合| 天天一区二区日本电影三级| 只有这里有精品99| 亚洲国产欧美人成| 亚洲欧美中文字幕日韩二区| 午夜精品一区二区三区免费看| 最近中文字幕2019免费版| 久久久成人免费电影| 成人无遮挡网站| 久久这里有精品视频免费| 亚洲国产精品sss在线观看| 欧美成人一区二区免费高清观看| 尾随美女入室| 国产亚洲91精品色在线| 精品久久久久久成人av| 久久欧美精品欧美久久欧美| 91精品伊人久久大香线蕉| 国产探花在线观看一区二区| 91aial.com中文字幕在线观看| 成人漫画全彩无遮挡| 成年免费大片在线观看| 中国国产av一级| 毛片女人毛片| 亚洲精品一区蜜桃| 国产伦在线观看视频一区| 亚洲真实伦在线观看| 99九九线精品视频在线观看视频| 国产麻豆成人av免费视频| 免费人成在线观看视频色| 1000部很黄的大片| 天堂√8在线中文| 亚洲国产精品成人久久小说| 秋霞伦理黄片| 久久热精品热| 国产精品三级大全| 久久久久网色| 日韩视频在线欧美| 欧美精品一区二区大全| 国产亚洲午夜精品一区二区久久 | 久久久午夜欧美精品| 99九九线精品视频在线观看视频| 99热这里只有是精品在线观看| 亚洲欧美成人精品一区二区| 午夜日本视频在线| 97在线视频观看| 建设人人有责人人尽责人人享有的 | 国产白丝娇喘喷水9色精品| 搡女人真爽免费视频火全软件| 能在线免费观看的黄片| 最近中文字幕2019免费版| 狂野欧美激情性xxxx在线观看| 九九爱精品视频在线观看| 纵有疾风起免费观看全集完整版 | 2021少妇久久久久久久久久久| 色5月婷婷丁香| 七月丁香在线播放| 晚上一个人看的免费电影| 亚洲中文字幕一区二区三区有码在线看| 听说在线观看完整版免费高清| 高清在线视频一区二区三区 | 99在线视频只有这里精品首页| 国产免费男女视频| 91久久精品国产一区二区三区| 99久久成人亚洲精品观看| 日产精品乱码卡一卡2卡三| 欧美xxxx黑人xx丫x性爽| 国产精品人妻久久久久久| 在线观看一区二区三区| 免费黄色在线免费观看| 寂寞人妻少妇视频99o| 国产高清有码在线观看视频| 国产精品不卡视频一区二区| АⅤ资源中文在线天堂| 久久精品国产鲁丝片午夜精品| 2022亚洲国产成人精品| 好男人在线观看高清免费视频| 伊人久久精品亚洲午夜| 我要看日韩黄色一级片| 啦啦啦观看免费观看视频高清| 色噜噜av男人的天堂激情| 人妻少妇偷人精品九色| 高清午夜精品一区二区三区| 99国产精品一区二区蜜桃av| 国产在线男女| 日韩三级伦理在线观看| 精品久久久久久久久久久久久| 2022亚洲国产成人精品| 午夜爱爱视频在线播放| 成人亚洲欧美一区二区av| 欧美性感艳星| 赤兔流量卡办理| 最近最新中文字幕免费大全7| 美女高潮的动态| 三级国产精品欧美在线观看| 一边摸一边抽搐一进一小说| 亚洲性久久影院| 免费电影在线观看免费观看| 午夜福利成人在线免费观看| 精品人妻视频免费看| 国产淫片久久久久久久久| 99久久九九国产精品国产免费| 男女国产视频网站| av专区在线播放| 七月丁香在线播放| 欧美成人精品欧美一级黄| 久久久亚洲精品成人影院| 丝袜美腿在线中文| 特级一级黄色大片| 狂野欧美激情性xxxx在线观看| 国产一区有黄有色的免费视频 | 又粗又硬又长又爽又黄的视频| 九九爱精品视频在线观看| 自拍偷自拍亚洲精品老妇| 国产老妇伦熟女老妇高清| 人人妻人人澡人人爽人人夜夜 | 97热精品久久久久久| 久久久久久久久久成人| 欧美性猛交╳xxx乱大交人| 青春草视频在线免费观看| 久久欧美精品欧美久久欧美| 免费观看性生交大片5| 色综合站精品国产| 91在线精品国自产拍蜜月| 久热久热在线精品观看| 看十八女毛片水多多多| 国产私拍福利视频在线观看| 99在线视频只有这里精品首页| 免费av观看视频| 精华霜和精华液先用哪个| 亚洲综合精品二区| 中文字幕av在线有码专区| 亚洲欧美日韩无卡精品| 亚洲精品影视一区二区三区av| 99热这里只有是精品50| 亚洲精品自拍成人| 亚洲欧美日韩卡通动漫| 久久久久久久亚洲中文字幕| 男女啪啪激烈高潮av片| 岛国在线免费视频观看| 日韩欧美国产在线观看| 国产乱人视频| 成人午夜高清在线视频| av在线老鸭窝| 婷婷色麻豆天堂久久 | 国产伦精品一区二区三区四那| av免费在线看不卡| videos熟女内射| www日本黄色视频网| 三级国产精品片| 22中文网久久字幕| 国产人妻一区二区三区在| 久久精品91蜜桃| 六月丁香七月| 国产精品1区2区在线观看.| 一夜夜www| 国产精品不卡视频一区二区| 成年女人永久免费观看视频| 亚洲av中文字字幕乱码综合| 国产精品日韩av在线免费观看| 久久久精品94久久精品| 久久人人爽人人爽人人片va| 秋霞在线观看毛片| 成年版毛片免费区| 又粗又硬又长又爽又黄的视频| 国产精品一二三区在线看| 女人久久www免费人成看片 | 又爽又黄无遮挡网站| 国产精品99久久久久久久久| 亚洲性久久影院| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 晚上一个人看的免费电影| 少妇高潮的动态图| 国产精品日韩av在线免费观看| 一级av片app| 久久久久久久国产电影| 在线观看66精品国产| 国产精品国产三级国产av玫瑰| 在线免费观看不下载黄p国产| av天堂中文字幕网| 亚洲国产欧美人成| 亚洲国产成人一精品久久久| 国产高潮美女av| 91在线精品国自产拍蜜月| 成人美女网站在线观看视频| 久久国产乱子免费精品| 毛片女人毛片| 久久亚洲国产成人精品v| 久久精品夜夜夜夜夜久久蜜豆| 在线免费十八禁| 免费一级毛片在线播放高清视频| av播播在线观看一区| 午夜免费激情av| av在线亚洲专区| 国产欧美另类精品又又久久亚洲欧美| 成年免费大片在线观看| av.在线天堂| 日本免费a在线| 国产色婷婷99| 男女下面进入的视频免费午夜| 国产精品一区www在线观看| 国语自产精品视频在线第100页| 久久精品影院6| 亚洲av不卡在线观看| 国产精品综合久久久久久久免费| 亚洲欧美日韩卡通动漫| 秋霞伦理黄片| 亚洲怡红院男人天堂| 久久人人爽人人片av| 你懂的网址亚洲精品在线观看 | 国产亚洲av片在线观看秒播厂 | 欧美性猛交黑人性爽| 亚洲丝袜综合中文字幕| 久久久久久久午夜电影| 97超碰精品成人国产| 亚洲精品国产成人久久av| 成人午夜精彩视频在线观看| 中文字幕熟女人妻在线| 麻豆成人av视频| 青青草视频在线视频观看| 欧美高清成人免费视频www| 国产亚洲午夜精品一区二区久久 | 蜜桃久久精品国产亚洲av| 欧美bdsm另类| 精品人妻一区二区三区麻豆| 丰满人妻一区二区三区视频av| 日韩欧美国产在线观看| 日韩欧美三级三区| 国产色婷婷99| 综合色av麻豆| 久久精品国产亚洲av涩爱| 国产黄片视频在线免费观看| 一个人看的www免费观看视频| 亚洲婷婷狠狠爱综合网| 久久精品91蜜桃| 精品欧美国产一区二区三| 国产免费一级a男人的天堂| 亚洲欧美日韩高清专用| 亚洲精品日韩在线中文字幕| 成年免费大片在线观看| 男人舔女人下体高潮全视频| 在线播放国产精品三级| 国产精品一区二区三区四区久久| 欧美变态另类bdsm刘玥| 视频中文字幕在线观看| 大香蕉久久网| 国模一区二区三区四区视频| 麻豆久久精品国产亚洲av| 亚洲成人久久爱视频| 午夜福利在线在线| 老司机福利观看| 亚洲av成人精品一二三区| 精品人妻一区二区三区麻豆| 天堂√8在线中文| 国产午夜精品一二区理论片| 我要搜黄色片| 免费av不卡在线播放| 男插女下体视频免费在线播放| 三级国产精品欧美在线观看| 非洲黑人性xxxx精品又粗又长| 国产av在哪里看| 超碰97精品在线观看| 久久精品国产自在天天线| av线在线观看网站| 中文精品一卡2卡3卡4更新| 免费一级毛片在线播放高清视频| 免费大片18禁| 欧美最新免费一区二区三区| 老司机影院成人| 久久精品人妻少妇| 女人被狂操c到高潮| 日韩精品有码人妻一区| 国产成人a区在线观看| 国产一区二区在线av高清观看| 赤兔流量卡办理| 免费观看人在逋| 九草在线视频观看| 美女cb高潮喷水在线观看| 久久久久久九九精品二区国产| 能在线免费观看的黄片| 亚洲av日韩在线播放| 春色校园在线视频观看| 国产精品不卡视频一区二区| 少妇高潮的动态图| 国产高清国产精品国产三级 | 国产精品国产高清国产av| 我要看日韩黄色一级片| 精品国内亚洲2022精品成人| 国产精品av视频在线免费观看| 超碰97精品在线观看| 久久99蜜桃精品久久| 日韩 亚洲 欧美在线| 精品熟女少妇av免费看| 一区二区三区乱码不卡18| 三级国产精品片| 国产一级毛片七仙女欲春2| 国产精品熟女久久久久浪| 日韩人妻高清精品专区| 晚上一个人看的免费电影| 亚洲国产精品专区欧美| 欧美97在线视频| 成人漫画全彩无遮挡| 色综合色国产| 国产老妇伦熟女老妇高清| 日本五十路高清| 久久久久久久国产电影| 亚洲国产欧美在线一区| 日本三级黄在线观看| 麻豆av噜噜一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 变态另类丝袜制服| 国产91av在线免费观看| 国产午夜精品久久久久久一区二区三区| 国产真实乱freesex| 夜夜看夜夜爽夜夜摸| 亚洲五月天丁香| 天天一区二区日本电影三级| 91久久精品国产一区二区三区| 中文字幕久久专区| 在线免费观看的www视频| 亚洲国产精品国产精品| av免费观看日本| 啦啦啦啦在线视频资源| 男人舔奶头视频| 国产真实乱freesex| 精品熟女少妇av免费看| 97热精品久久久久久| 成年版毛片免费区| 丝袜美腿在线中文| 日韩国内少妇激情av| 日本免费a在线| 亚洲成人av在线免费| 国产成人aa在线观看| www日本黄色视频网| 男女那种视频在线观看| 国产精品精品国产色婷婷| 国产精品一区二区在线观看99 | 国产单亲对白刺激| 99久久中文字幕三级久久日本| 国产在线男女| 毛片一级片免费看久久久久| 日韩欧美 国产精品| 亚洲精品456在线播放app| 亚洲伊人久久精品综合 | 亚洲精品乱久久久久久| 日日摸夜夜添夜夜添av毛片| 久久午夜福利片| 中文天堂在线官网| 久久久久久久国产电影| 变态另类丝袜制服| 岛国毛片在线播放| av女优亚洲男人天堂| 日本午夜av视频| 成人美女网站在线观看视频| 日本免费a在线| 欧美三级亚洲精品| 女人被狂操c到高潮| 美女黄网站色视频| 国产午夜精品论理片| 亚洲aⅴ乱码一区二区在线播放| 伦理电影大哥的女人| 美女内射精品一级片tv| 欧美区成人在线视频| 午夜精品在线福利| 日韩一本色道免费dvd| 99热6这里只有精品| 97人妻精品一区二区三区麻豆| 精品无人区乱码1区二区| av播播在线观看一区| 韩国av在线不卡| 2021天堂中文幕一二区在线观| 精品人妻熟女av久视频| 22中文网久久字幕| 国产成人aa在线观看| 国产乱人视频| 少妇人妻精品综合一区二区| 国产乱人视频| 免费黄色在线免费观看| 神马国产精品三级电影在线观看| 午夜激情福利司机影院| 2022亚洲国产成人精品| 青春草国产在线视频| 亚洲真实伦在线观看| 免费看av在线观看网站| av播播在线观看一区| 啦啦啦啦在线视频资源| 日本色播在线视频| 午夜精品国产一区二区电影 | 蜜桃亚洲精品一区二区三区| 亚洲av免费在线观看| 国产免费男女视频| 国产片特级美女逼逼视频| 少妇被粗大猛烈的视频| 看片在线看免费视频| 精品久久久久久久末码| 久久99精品国语久久久| 在线播放国产精品三级| 18禁在线播放成人免费| 国产成人精品久久久久久| 99热这里只有是精品在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 三级毛片av免费| 搡女人真爽免费视频火全软件| 中文字幕av在线有码专区| 丰满少妇做爰视频| 亚洲成人精品中文字幕电影| 亚洲三级黄色毛片| 毛片一级片免费看久久久久| 国产极品天堂在线| 国产又黄又爽又无遮挡在线| 日本免费a在线| 国产高潮美女av| 久久精品国产亚洲av天美| 噜噜噜噜噜久久久久久91| 夜夜看夜夜爽夜夜摸| 国产精品久久久久久久久免| 国产成人精品婷婷| 亚洲经典国产精华液单| 青春草国产在线视频| av又黄又爽大尺度在线免费看 | 一级二级三级毛片免费看| 午夜福利在线观看免费完整高清在| 久久精品国产99精品国产亚洲性色| 天堂√8在线中文| 国语自产精品视频在线第100页| 99热全是精品| 国产精品人妻久久久影院| 国产精品一区二区三区四区免费观看| 热99re8久久精品国产| 色吧在线观看| 亚洲欧美日韩东京热| 最近2019中文字幕mv第一页| av免费在线看不卡| 免费看光身美女| 欧美xxxx性猛交bbbb| 美女黄网站色视频| www.av在线官网国产| 中文字幕精品亚洲无线码一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费看av在线观看网站| 变态另类丝袜制服| 色网站视频免费| 蜜桃亚洲精品一区二区三区| 97在线视频观看| 你懂的网址亚洲精品在线观看 | 久久精品影院6| av福利片在线观看| 一级av片app| 色综合站精品国产| 免费观看在线日韩| 亚洲久久久久久中文字幕| av在线天堂中文字幕| 纵有疾风起免费观看全集完整版 | 人人妻人人看人人澡| 亚洲精品国产av成人精品| 看非洲黑人一级黄片| 国产亚洲精品av在线| av女优亚洲男人天堂| 国产精品一及|