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

    單面柱聲子晶體板低頻帶隙特性與機(jī)理分析

    2016-11-25 20:39:43張昭蘇開創(chuàng)韓星凱
    計(jì)算機(jī)輔助工程 2016年5期
    關(guān)鍵詞:單胞散射體截止頻率

    張昭++蘇開創(chuàng)+韓星凱

    摘要: 基于有限元法計(jì)算分析散射體材料屬性和散射體形狀對(duì)單面柱板結(jié)構(gòu)局域共振型聲子晶體帶隙特性的影響,并通過計(jì)算單胞位移特征模式解釋聲子晶體帶隙特性變化的物理機(jī)理.結(jié)果表明:散射體的密度和彈性模量以及散射體形狀通過改變聲子晶體單胞局域狀態(tài),對(duì)聲子晶體帶隙特性產(chǎn)生很大影響.總結(jié)聲子晶體板結(jié)構(gòu)低頻帶隙特性變化規(guī)律,為工程減振降噪提供參考.

    關(guān)鍵詞: 聲子晶體; 散射體; 基體; 帶隙特性; 位移特征模式; 密度; 彈性模量; 減振降噪

    中圖分類號(hào): O481.1; TB535 文獻(xiàn)標(biāo)志碼: A

    0 引 言

    聲子晶體是由一種或幾種彈性材料周期排列在基底材料中形成的一類人工復(fù)合材料.[1]周期性復(fù)合材料或結(jié)構(gòu)的聲子晶體可以阻擋某些頻率的彈性波傳播,這種現(xiàn)象稱為聲子晶體帶隙或禁帶.[2]頻率落在帶隙范圍內(nèi)的彈性波理論上在聲子晶體中被禁止傳播,對(duì)應(yīng)到有限尺寸的聲子晶體周期結(jié)構(gòu),帶隙頻率范圍內(nèi)的彈性波在傳播過程中會(huì)產(chǎn)生指數(shù)衰減.[3]研究認(rèn)為,聲子晶體帶隙產(chǎn)生的機(jī)理有2種:布拉格散射型和局域共振型.[4]前者主要是結(jié)構(gòu)的周期性起主導(dǎo)作用,當(dāng)入射彈性波的波長(zhǎng)與結(jié)構(gòu)的特征長(zhǎng)度(晶格常數(shù))相近時(shí),將受到結(jié)構(gòu)強(qiáng)烈的散射;后者則是單個(gè)散射體的共振特性起主導(dǎo)作用.[5]LIU等[6]提出局域共振型聲子晶體理論.根據(jù)局域共振理論,相同特征長(zhǎng)度的局域共振型聲子晶體帶隙頻率要低于布拉格散射型聲子晶體帶隙頻率2個(gè)數(shù)量級(jí)[7],這使得聲子晶體小尺寸控制大波長(zhǎng)成為可能.為實(shí)現(xiàn)低頻減振降噪領(lǐng)域聲子晶體的實(shí)際應(yīng)用,目前出現(xiàn)針對(duì)各種聲子晶體結(jié)構(gòu)[8-11]的研究,其中關(guān)于Lamb波在聲子晶體板結(jié)構(gòu)中傳播[12-14]的研究尤其引人注意,特別是橡膠吸聲層聲子晶體結(jié)構(gòu)廣泛應(yīng)用于潛艇等水下裝備,用于吸收干擾聲波或噪聲[15].YU等[16]提出一種極具應(yīng)用前景的局域共振型單面柱板結(jié)構(gòu)聲子晶體,對(duì)該結(jié)構(gòu)聲子晶體低頻帶隙特性受部分幾何參數(shù)影響進(jìn)行計(jì)算研究和物理解釋.除聲子晶體單胞幾何參數(shù)外,材料參數(shù)和散射體形狀對(duì)聲子晶體帶隙特性也可能有很大影響,對(duì)其進(jìn)行研究將非常有利于低頻帶隙特性聲子晶體的實(shí)際應(yīng)用.

    本文采用有限元法,使用多物理場(chǎng)分析軟件COMSOL建立單面柱板結(jié)構(gòu)聲子晶體單胞模型,計(jì)算分析聲子晶體散射體材料密度、彈性模量和改變射體形狀對(duì)聲子晶體帶隙特性的影響;同時(shí)使用多物理場(chǎng)分析軟件COMSOL計(jì)算聲子晶體單胞位移特征模式,給出聲子晶體帶隙特性變化的物理解釋.這些計(jì)算和分析可為這種極具前景的聲子晶體板結(jié)構(gòu)的實(shí)際應(yīng)用鋪平道路,也為其他類似具有低頻帶隙特性的局域共振型聲子晶體結(jié)構(gòu)的設(shè)計(jì)和應(yīng)用提供參考.

    1 聲子晶體帶隙計(jì)算的有限元法

    利用多物理場(chǎng)分析軟件COMSOL求解彈性波在聲子晶體中的波動(dòng)方程.由于聲子晶體是周期結(jié)構(gòu),依照Bloch定理,計(jì)算可在一個(gè)單胞內(nèi)完成.[17]對(duì)模型劃分網(wǎng)格后,單胞離散形式的特征值方程為

    2 帶隙特性計(jì)算和位移特征模式分析

    首先,將本文計(jì)算結(jié)果與文獻(xiàn)[16]對(duì)比驗(yàn)證本文計(jì)算方法準(zhǔn)確性.圓柱散射體聲子晶體單胞有限元模型見圖1,單胞基體底面(正方形)邊長(zhǎng)a=10 mm,基體厚度e=2 mm,散射體高度h=10 mm,圓柱散射體半徑r=4.5 mm.基體所用材料為橡膠,散射體所用材料為鋼,材料參數(shù)見表1.

    使用多物理場(chǎng)分析軟件COMSOL計(jì)算得到的聲子晶體帶隙特性見圖2.圓柱散射體半徑r=4.5 mm時(shí)計(jì)算得到的聲子晶體第一帶隙數(shù)據(jù)與文獻(xiàn)[16]對(duì)比結(jié)果見表2.由此可以看出:第一帶隙中心頻率、帶寬和截止頻率誤差都很小,驗(yàn)證本文模型和計(jì)算方法準(zhǔn)確可信.

    取散射體半徑r=4 mm,單胞模型其他幾何尺寸不變.單胞的基體材料為橡膠,散射體材料以鋼為基準(zhǔn),改變散射體材料的密度,計(jì)算分析聲子晶體帶隙特性變化規(guī)律,結(jié)果見圖3.

    由圖3可知:散射體密度對(duì)第一帶隙截止頻率影響很小幾乎可以忽略,但對(duì)第一帶隙起始頻率影響很大;當(dāng)散射體密度很小時(shí),幾乎沒有帶隙產(chǎn)生;而散射體密度足夠大時(shí),就會(huì)產(chǎn)生明顯帶隙,且散射體密度大于7 000 kg/m3后,第一帶隙起始頻率下降變得平緩,產(chǎn)生的低頻帶隙帶寬慢慢趨于穩(wěn)定.

    計(jì)算聲子晶體單胞的位移特征模式,進(jìn)而分析散射體材料的密度變化影響聲子晶體帶隙特性的物理機(jī)制.分別計(jì)算散射體密度為1 000,3 000,5 000,11 000,15 000和19 000 kg/m3時(shí)的聲子晶體帶隙(第一帶隙)起始頻率和截止頻率頂點(diǎn)對(duì)應(yīng)波矢點(diǎn)的6組單胞位移特征模式,見表3.

    由此可以發(fā)現(xiàn):隨著散射體密度逐漸增大,起始頻率頂點(diǎn)對(duì)應(yīng)波矢點(diǎn)的單胞振動(dòng)能量集中從基體上轉(zhuǎn)移到散射體上,呈現(xiàn)明顯的局域狀態(tài),而截止頻率頂點(diǎn)的單胞位移特征模式幾乎沒有變化.由此可以解釋上面計(jì)算得到的帶隙特性隨散射體密度變化規(guī)律(見圖3):當(dāng)散射體密度增大時(shí),根據(jù)簡(jiǎn)單彈簧模型公式f=k/m/2π,起始頻率的單胞散射體振子(振動(dòng)能量集中的部分)質(zhì)量逐漸變大而截止頻率的單胞基體振子(此時(shí)單胞振動(dòng)能量總是集中在單胞基體的4個(gè)邊角上)質(zhì)量幾乎不變,彈簧剛度總是散射體和4個(gè)基體角邊相連接的基體部分的剛度且?guī)缀醣3植蛔?,故聲子晶體帶隙起始頻率降低而截止頻率基本不變,低頻帶隙曲線打開變寬.

    仍取散射體半徑r=4 mm,單胞模型其他幾何尺寸不變.單胞的基體材料為橡膠,散射體材料以鋼為基準(zhǔn),改變散射體材料的密度,計(jì)算分析聲子晶體帶隙特性變化規(guī)律,結(jié)果見圖4.

    由圖4可知:改變散射體彈性模量,其他保持不變,第一帶隙起始頻率和截止頻率隨彈性模量的增

    大逐漸變大.當(dāng)散射體彈性模量大過橡膠基體的彈

    性模量1.175×105 Pa后,截止頻率的變大速率明顯更快,產(chǎn)生的低頻帶隙帶寬明顯變大.當(dāng)散射體彈性模量大于1.175×107 Pa后,聲子晶體帶隙趨于穩(wěn)

    定.

    計(jì)算散射體彈性模量為1.175×103,1.175×104,1.175×105,1.175×106,1.175×107,1.175×108 Pa時(shí)的聲子晶體帶隙(第一帶隙)起始頻率和截止頻率頂點(diǎn)對(duì)應(yīng)波矢點(diǎn)的6組單胞位移特征模式,見表4.由此可以看出:隨著散射體彈性模量逐漸增大,起始頻率頂點(diǎn)對(duì)應(yīng)波矢點(diǎn)的單胞位移特征模式單胞振子在散射體上由大變小.彈簧部分是與散射體振子相連接的基體橡膠,其彈性模量不變,根據(jù)簡(jiǎn)單彈簧模型公式f=k/m/2π,m變小而k不變.截止頻率單胞振子由鋼散射體轉(zhuǎn)移到橡膠基體的4個(gè)角上,彈簧部分總是位于與散射體相連接的橡膠基體上,故m變小且變化很大,k不變.這也可解釋前文計(jì)算得到的聲子晶體帶隙特性隨散射體彈性模量的變化規(guī)律(見圖4),當(dāng)散射體彈性模量變大時(shí),起始頻率和截止頻率均變大但截止頻率變化更顯著,低頻帶隙打開變寬.

    除圓柱散射體幾何尺寸大?。▓A柱散射體半徑和高度)會(huì)影響聲子晶體帶隙特性外,散射體的幾何形狀也可能會(huì)對(duì)聲子晶體帶隙特性有重要影響.為此,把散射體幾何形狀由圓柱體改為長(zhǎng)方體,單胞有限元模型見圖5,長(zhǎng)方散射體的高度h=10 mm保持不變,改變長(zhǎng)方散射體的底面(正方形)邊長(zhǎng),計(jì)算得到帶隙特性變化規(guī)律見圖6,且與圓柱散射體聲子晶體帶隙特性作對(duì)比,見圖7.

    由圖6和7可以看到:散射體底面邊長(zhǎng)對(duì)第一帶隙起始頻率影響不大,但對(duì)第一帶隙截止頻率有很大影響,與圓柱散射體情況下圓柱半徑對(duì)第一帶隙的影響一致.

    對(duì)比發(fā)現(xiàn),散射體的高度和體積相同時(shí),2種聲子晶體的第一帶隙起始頻率幾乎相同,但長(zhǎng)方散射體聲子晶體的第一帶隙截止頻率明顯大于圓柱散射體聲

    子晶體,因此長(zhǎng)方散射體聲子晶體比圓柱散射體聲子晶體有更好的低頻帶隙特性.

    聲子晶體帶隙特性,二者起始頻率和截止頻率對(duì)比見表5.由此可知:等高等體積的長(zhǎng)方散射體聲子晶體帶隙特性明顯優(yōu)于圓柱散射體聲子晶體.計(jì)算2種散射體截面的剛度見表6,分別對(duì)比二者起始頻率帶隙頂點(diǎn)對(duì)應(yīng)波矢點(diǎn)的單胞位移特征模式和截止頻率帶隙頂點(diǎn)對(duì)應(yīng)波矢點(diǎn)的單胞位移特征模式見表7.由表7可知,關(guān)于起始頻率帶隙頂點(diǎn)對(duì)應(yīng)波矢點(diǎn)的單胞位移特征模式,二者都是彎曲振動(dòng),但由于長(zhǎng)方散射體截面抗彎剛度略大于圓柱散射體截面抗彎剛度,導(dǎo)致前者帶隙起始頻率略大于后者且相差不大;關(guān)于截止頻率帶隙頂點(diǎn)對(duì)應(yīng)波矢點(diǎn)的單胞位移特征模式,長(zhǎng)方散射體聲子晶體單胞主要是基體的扭轉(zhuǎn)和彎曲振動(dòng),圓柱散射體聲子晶體單胞主要是基體的彎曲振動(dòng),故前者的總剛度明顯大于后者,導(dǎo)致長(zhǎng)方散射體聲子晶體帶隙截止頻率明顯高于圓柱散射體聲子晶體帶隙截止頻率.

    3 聲學(xué)應(yīng)用實(shí)例

    綜合前文研究結(jié)果,在實(shí)際常用材料中選用最優(yōu)材料組合,采用合適的聲子晶體單胞構(gòu)型,得到低頻帶隙性能最優(yōu)的聲子晶體板用于聲學(xué)降噪隔聲.漢語(yǔ)的聲能密度一般都集中在200~700 Hz,實(shí)際常用材料見表8.

    為得到低頻彈性波帶隙,基體材料采用橡膠,分別選用鎢、銅、鋼、鋁作為散射體材料;聲子晶體單胞采用低頻帶隙特性更好的長(zhǎng)方散射體單胞構(gòu)型,散射體底邊邊長(zhǎng)b=9 mm,其他幾何參數(shù)不變.計(jì)算得到4種聲子晶體帶隙特性見表9.對(duì)比發(fā)現(xiàn):鋁作為散射體材料得到的聲子晶體帶隙特性明顯不如另外3種材料;鎢、銅、鋼分別作為散射體材料時(shí),都能符合漢語(yǔ)聲能密度集中在200~700 Hz頻率范圍的降噪隔聲要求.銅和鋼的密度遠(yuǎn)小于鎢,鋼遠(yuǎn)比銅經(jīng)濟(jì)且輕質(zhì),因此,基體材料為橡膠、散射體材料為鋼的長(zhǎng)方散射體單面柱板結(jié)構(gòu)聲子晶體是降噪隔聲最佳方案.

    取8×8聲子晶體有限周期板考察振動(dòng)傳輸特性,驗(yàn)證單面柱板結(jié)構(gòu)聲子晶體設(shè)計(jì)和帶隙特性計(jì)算的有效性.選擇基體材料為橡膠、散射體材料為鋼的長(zhǎng)方散射體聲子晶體板,8×8聲子晶體板結(jié)構(gòu)幾何模型見圖8.使用COMSOL軟件進(jìn)行頻率響應(yīng)分析,從板的一端輸入幅值為單位值的位移頻譜,x軸正方向?yàn)槲灰祁l譜方向,頻率范圍為100~1 500 Hz,得到板另一端即輸出端的位移響應(yīng)平均值,結(jié)果見圖9.

    由圖9可知:帶隙頻率范圍內(nèi)聲子晶體周期板減振效果明顯,且頻響分析減振頻率范圍和前文計(jì)算得出的帶隙范圍幾乎重合,驗(yàn)證本文單面柱板結(jié)構(gòu)聲子晶體設(shè)計(jì)和帶隙計(jì)算的有效性.

    4 結(jié) 論

    (1)單面柱板結(jié)構(gòu)聲子晶體材料參數(shù)中,散射體的密度和彈性模量足夠大時(shí),都會(huì)產(chǎn)生優(yōu)良的低頻帶隙,二者對(duì)聲子晶體帶隙特性有很大影響.

    (2)將單面柱板結(jié)構(gòu)聲子晶體散射體形狀改變?yōu)殚L(zhǎng)方體,發(fā)現(xiàn)等高等體積散射體條件下,長(zhǎng)方散射體聲子晶體比圓柱散射體聲子晶體具有更好的低頻帶隙特性.

    (3)從聲子晶體單胞位移特征模式可以得出,聲子晶體帶隙特性隨散射體材料密度和彈性模量變化規(guī)律可以用簡(jiǎn)單彈簧振動(dòng)模型來解釋,簡(jiǎn)單彈簧模型的振子質(zhì)量發(fā)生變化(彈簧剛度幾乎不變),進(jìn)而導(dǎo)致聲子晶體帶隙特性改變.改變散射體形狀,保持散射體體積和材料常數(shù)不變,結(jié)合聲子晶體帶隙頂點(diǎn)對(duì)應(yīng)波矢點(diǎn)的單胞位移特征模式,可以認(rèn)為單胞振動(dòng)截面剛度的變化導(dǎo)致聲子晶體帶隙的變化.

    (4)計(jì)算對(duì)比得出以鋼為散射體、橡膠為基體的長(zhǎng)方散射體單面柱板結(jié)構(gòu)聲子晶體能很好地應(yīng)用于實(shí)際聲學(xué)低頻降噪隔聲,說明合理選用材料組合和單胞構(gòu)型,聲子晶體可以在工程實(shí)際中發(fā)揮重要作用.

    參考文獻(xiàn):

    [1]

    溫熙森, 溫激鴻, 郁殿龍, 等. 聲子晶體[M]. 北京: 國(guó)防工業(yè)出版社, 2009.

    [2] 郁殿龍. 基于聲子晶體理論的梁板類周期結(jié)構(gòu)振動(dòng)帶隙特性研究[D]. 長(zhǎng)沙: 國(guó)防科學(xué)技術(shù)大學(xué), 2006.

    [3] 劉耀宗, 孟浩, 李黎, 等. 基于遺傳算法的聲子晶體梁振動(dòng)傳輸特性優(yōu)化設(shè)計(jì)[J]. 振動(dòng)與沖擊, 2008, 27(9): 47-50. DOI: 10.3969/j.issn.1000-3835.2008.09.012.

    LIU Y Z, MENG H, LI L, et al. Optical design of phononic crystal beams for attenuating vibration propagation[J]. Journal of Vibration and Shock, 2008, 27(9): 47-50. DOI: 10.3969/j.issn.1000-3835.2008.09.012.

    [4] 舒海生, 劉少剛, 王威遠(yuǎn), 等. 集中質(zhì)量邊界條件下聲子晶體桿的縱向振動(dòng)傳遞特性研究[J]. 振動(dòng)與沖擊, 2012, 31(19): 113-117. DOI: 10.3969/j.issn.1000-3835.2012.19.025.

    SHU H S, LIU S G, WANG W Y, et al. Transmission characteristics of longitudinal vibration of a phononic crystal rod with concentrated mass boundary condition[J]. Journal of Vibration and Shock, 2012, 31(19): 113-117. DOI: 10.3969/j.issn.1000-3835.2012.19.025

    [5] ZHAO H, LIU Y, WANG G, et al. Resonance modes and gap formation in a two-dimensional solid phononic crystal[J]. Physical Review B: Condensed Matter, 2005, 72(1): 012301. DOI: 10.1103/PhysRevB.72.012301.

    [6] LIU Z Y, ZHANG X X, MAO Y W, et al. Locally resonant sonic material[J]. Science, 2000, 289(5485): 1734-1736. DOI: 10.1126/science.289.5485.1734

    [7] WANG P, CHEN T N, YU K P, et al. Lamb wave band gaps in a double-sided phononic plate[J]. Journal of Applied Physics, 2013, 113(5): 053509. DOI: 10.1063/1.4790301.

    [8] LIU Z, CHAN C T, SHENG P. Three-component elastic wave band-gap material[J]. Physical Review B: Condensed Matter, 2002, 65(16): 165116. DOI: 10.1103/PhysRevB.65.165116.

    [9] OUDICH M, LI Y, ASSOUAR B M, et al. A sonic band gap based on the locally resonant phononic plates with stubs[J]. New Journal of Physics, 2010, 12(2): 201-206. DOI: 10.1088/1367-2630/12/8/083049.

    [10] HSU J C, WU T T. Lamb waves in binary locally resonant phononic plates with two-dimensional lattices[J]. Applied Physics Letters, 2007, 90(20): 201904-1-201904-3. DOI: 10.1063/1.2739369.

    [11] HSU J C. Local resonances-induced low-frequency band gaps in two-dimensional phononic crystal slabs with periodic stepped resonators[J]. Journal of Physics D: Applied Physics, 2011, 44(5): 055401. DOI: 10.1088/0022-3727/44/5/055401.

    [12] CHENG Y, LIU X J. Three dimensional multilayered acoustic cloak with homogeneous isotropic materials[J]. Applied Physics A, 2009, 94(1): 25-30. DOI: 10.1007/s00339-008-4882-7.

    [13] HSU J C, WU T T. Lamb waves in binary locally resonant phononic plates with two-dimensional lattices[J]. Applied Physics Letters, 2007, 90(20): 201904. DOI: 10.1063/1.2739369.

    [14] 陳阿麗, 郭鳳丹, 汪越勝. 一維準(zhǔn)周期聲子晶體板中蘭姆波的禁帶特性分析[J]. 功能材料, 2015, 46(19): 19046-19051. DOI: 10.3969/j.issn.1001-9731.2015.19.010.

    CHEN A L, GUO D F, WANG Y S. Study on the band gaps of the Lamb wave in one-dimensional quasi periodic phononic crystal plates[J]. Journal of Functional Materials, 2015, 46(19): 19046-19051. DOI: 10.3969/j.issn.1001-9731.2015.19.010.

    [15] 趙宏剛, 溫激鴻, 楊海濱, 等. 一種含柱形空腔結(jié)構(gòu)橡膠層的吸聲機(jī)理及優(yōu)化[J]. 物理學(xué)報(bào), 2014, 63(13): 134303-1-134303-6. DOI: 10.7498/aps.63.134303.

    ZHAO H G, WEN J H, YANG H B, et al. Acoustic absorption mechanism and optimization of a rubber slab with cylindrical cavities[J]. Acta Physica Sinica, 2014, 63(13): 134303. DOI: 10.7498/aps.63.134303.

    [16] YU K, CHEN T, WANG X. Band gaps in the low-frequency range based on the two-dimensional phononic crystal plates composed of rubber matrix with periodic steel stubs[J]. Physica B: Condensed Matter, 2013, 416: 12-16. DOI: 10.1016/j.physb.2013.02.011.

    [17] 黃佳, 尹進(jìn), 張昭, 等. 圓管型局域共振聲子晶體三維構(gòu)型振動(dòng)帶隙研究[J]. 計(jì)算力學(xué)學(xué)報(bào), 2015, 32(3): 353-358. DOI: 10.7511/jslx201503009.

    HUANG J, YIN J, ZHANG Z, et al. Study on the 3D polarization band gaps of pipe-shaped locally resonant phononic crystals[J]. Chinese Journal of Computational Mechanics, 2015, 32(3): 353-358. DOI: 10.7511/jslx201503009.

    [18] 黃佳, 尹進(jìn), 張昭, 等. 基于聲子晶體的層疊式方柱型減振結(jié)構(gòu)設(shè)計(jì)[J]. 計(jì)算機(jī)輔助工程, 2015, 24(2): 53-58. DOI: 10.13340/j.cae.2015.02.010.

    HUANG J, YIN J, ZHANG Z, et al. Vibration of reduction structure design with stacked square columns based on phononic crystal[J]. Computer Aided Engineering, 2015, 24(2): 53-58. DOI: 10.13340/j.cae.2015.02.010.

    猜你喜歡
    單胞散射體截止頻率
    基于規(guī)范圖像的光電成像系統(tǒng)采樣響應(yīng)研究
    激光與紅外(2023年8期)2023-09-22 09:01:10
    基于超聲Lamb波截止頻率的雙層薄板各層厚度表征
    基于NURBS的點(diǎn)陣材料參數(shù)化建模方法
    一種基于單次散射體定位的TOA/AOA混合定位算法*
    復(fù)合材料周期結(jié)構(gòu)數(shù)學(xué)均勻化方法的一種新型單胞邊界條件
    二維結(jié)構(gòu)中亞波長(zhǎng)缺陷的超聲特征
    低頻射頻識(shí)別系統(tǒng)中的RC放大器電路性能分析與研究
    考慮界面層影響的三維機(jī)織復(fù)合材料單胞模型研究
    高斯波包散射體成像方法
    梯度飽和多孔材料中彈性波的截止頻率
    国产精品嫩草影院av在线观看| 免费不卡的大黄色大毛片视频在线观看| 久久久久久久大尺度免费视频| 在线精品无人区一区二区三 | 亚洲欧美日韩卡通动漫| 九色成人免费人妻av| 永久网站在线| av福利片在线观看| 九九在线视频观看精品| 亚洲av综合色区一区| 日本一二三区视频观看| 美女xxoo啪啪120秒动态图| 国产男女超爽视频在线观看| kizo精华| 亚洲精品第二区| 1000部很黄的大片| 性色av一级| 国产精品久久久久久久电影| 精品久久久精品久久久| 国产一区二区三区综合在线观看 | 特大巨黑吊av在线直播| 欧美日韩视频高清一区二区三区二| 国产免费一级a男人的天堂| 久久久久国产精品人妻一区二区| 久久人人爽av亚洲精品天堂 | 在线 av 中文字幕| 男人和女人高潮做爰伦理| 日韩大片免费观看网站| 在线播放无遮挡| 自拍偷自拍亚洲精品老妇| 国产毛片在线视频| 国产精品国产av在线观看| 最近2019中文字幕mv第一页| 偷拍熟女少妇极品色| 久久精品国产自在天天线| 纵有疾风起免费观看全集完整版| av视频免费观看在线观看| 欧美 日韩 精品 国产| 国产成人a∨麻豆精品| 777米奇影视久久| 久久女婷五月综合色啪小说| videos熟女内射| 国产精品久久久久久久电影| 亚洲性久久影院| 在线精品无人区一区二区三 | 天堂中文最新版在线下载| 亚洲欧美精品自产自拍| 久久国产精品男人的天堂亚洲 | 最近手机中文字幕大全| 国产精品久久久久久av不卡| 啦啦啦视频在线资源免费观看| 亚洲欧洲日产国产| 网址你懂的国产日韩在线| 欧美+日韩+精品| 九草在线视频观看| 十八禁网站网址无遮挡 | 在线免费十八禁| 精品一区二区免费观看| 国内少妇人妻偷人精品xxx网站| 亚洲欧美一区二区三区黑人 | 亚洲美女视频黄频| 亚洲精品日韩在线中文字幕| 亚洲久久久国产精品| 亚洲,欧美,日韩| 国产黄片美女视频| 日韩,欧美,国产一区二区三区| 免费黄色在线免费观看| 国产精品熟女久久久久浪| 观看av在线不卡| 少妇丰满av| 十分钟在线观看高清视频www | 欧美+日韩+精品| 三级经典国产精品| 亚洲av日韩在线播放| 波野结衣二区三区在线| 熟女av电影| 黄色视频在线播放观看不卡| 男男h啪啪无遮挡| 国精品久久久久久国模美| 亚洲一级一片aⅴ在线观看| 国产乱人视频| 1000部很黄的大片| 尾随美女入室| 乱码一卡2卡4卡精品| 日韩一区二区三区影片| av国产久精品久网站免费入址| 久久ye,这里只有精品| 一边亲一边摸免费视频| 免费观看a级毛片全部| 国产熟女欧美一区二区| 岛国毛片在线播放| 高清在线视频一区二区三区| 亚洲精品乱码久久久v下载方式| 成人二区视频| 国产深夜福利视频在线观看| 亚洲人成网站高清观看| 久久这里有精品视频免费| 久久久久国产网址| 三级经典国产精品| 新久久久久国产一级毛片| 国产精品一二三区在线看| 下体分泌物呈黄色| 精品久久久精品久久久| 久久国内精品自在自线图片| 一区二区三区免费毛片| 最近手机中文字幕大全| 亚洲精品成人av观看孕妇| 麻豆成人午夜福利视频| 国产成人精品婷婷| 自拍欧美九色日韩亚洲蝌蚪91 | 22中文网久久字幕| 国产欧美亚洲国产| 成人午夜精彩视频在线观看| 毛片一级片免费看久久久久| 亚洲精品国产av成人精品| 国产精品熟女久久久久浪| 在现免费观看毛片| 丰满少妇做爰视频| 性色avwww在线观看| 一本一本综合久久| 少妇人妻精品综合一区二区| 毛片一级片免费看久久久久| 久久久久网色| 免费黄色在线免费观看| 久久久久久久久久久丰满| 国产日韩欧美亚洲二区| 美女中出高潮动态图| 国产精品人妻久久久久久| 久久精品国产亚洲av天美| 一级黄片播放器| 麻豆精品久久久久久蜜桃| 性色avwww在线观看| 精品少妇久久久久久888优播| 国国产精品蜜臀av免费| 亚洲成色77777| 波野结衣二区三区在线| 亚洲性久久影院| 国产无遮挡羞羞视频在线观看| 欧美成人一区二区免费高清观看| 男人舔奶头视频| 免费黄频网站在线观看国产| 观看免费一级毛片| 水蜜桃什么品种好| 久久97久久精品| 国产精品一区二区在线不卡| 丰满乱子伦码专区| 在线观看三级黄色| 色吧在线观看| 天天躁夜夜躁狠狠久久av| 九草在线视频观看| 国产精品一区二区在线观看99| 人人妻人人澡人人爽人人夜夜| 成人毛片a级毛片在线播放| 91午夜精品亚洲一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 久久影院123| 欧美zozozo另类| 一本久久精品| 18禁在线播放成人免费| 男男h啪啪无遮挡| 在线观看三级黄色| 最新中文字幕久久久久| 亚洲自偷自拍三级| 欧美zozozo另类| 亚洲av免费高清在线观看| 六月丁香七月| 亚洲欧美一区二区三区国产| 男人和女人高潮做爰伦理| 偷拍熟女少妇极品色| 久久久色成人| 一本久久精品| 国产一级毛片在线| 99热这里只有是精品50| 久久毛片免费看一区二区三区| 九草在线视频观看| 国产乱人视频| 亚洲自偷自拍三级| 欧美少妇被猛烈插入视频| 成人特级av手机在线观看| 国产高清有码在线观看视频| 99久国产av精品国产电影| 在线亚洲精品国产二区图片欧美 | 一级毛片黄色毛片免费观看视频| 国产有黄有色有爽视频| 亚洲精品乱码久久久久久按摩| 成人国产麻豆网| 狂野欧美激情性xxxx在线观看| 亚洲婷婷狠狠爱综合网| 中文乱码字字幕精品一区二区三区| 国产爽快片一区二区三区| 国产伦精品一区二区三区四那| 蜜桃久久精品国产亚洲av| 中文字幕亚洲精品专区| 一级毛片aaaaaa免费看小| 国产色婷婷99| 欧美激情极品国产一区二区三区 | 97在线人人人人妻| 国产精品国产av在线观看| 成人午夜精彩视频在线观看| 久久久a久久爽久久v久久| 青春草视频在线免费观看| 国产高清不卡午夜福利| 久久国产乱子免费精品| 成人漫画全彩无遮挡| 美女福利国产在线 | 国产在线一区二区三区精| 日韩欧美一区视频在线观看 | 亚洲欧美日韩另类电影网站 | 天堂俺去俺来也www色官网| 三级国产精品片| 亚洲不卡免费看| av在线老鸭窝| 午夜免费观看性视频| 少妇的逼水好多| 亚洲中文av在线| 99精国产麻豆久久婷婷| 午夜福利在线在线| 又大又黄又爽视频免费| 男女国产视频网站| 欧美性感艳星| 51国产日韩欧美| 丝瓜视频免费看黄片| 亚洲欧美成人精品一区二区| 日本免费在线观看一区| 国产成人精品久久久久久| 亚洲精品久久久久久婷婷小说| 爱豆传媒免费全集在线观看| 在线观看国产h片| 精品少妇久久久久久888优播| 国产精品久久久久成人av| 高清欧美精品videossex| 尤物成人国产欧美一区二区三区| 伊人久久国产一区二区| 精品亚洲乱码少妇综合久久| 国产精品99久久久久久久久| 身体一侧抽搐| 日韩不卡一区二区三区视频在线| 欧美性感艳星| 下体分泌物呈黄色| 久久婷婷青草| 我的老师免费观看完整版| 日产精品乱码卡一卡2卡三| 九九久久精品国产亚洲av麻豆| 免费大片18禁| 免费在线观看成人毛片| h视频一区二区三区| 免费观看在线日韩| 色婷婷av一区二区三区视频| 最近中文字幕高清免费大全6| 少妇熟女欧美另类| 精品久久久噜噜| 国产v大片淫在线免费观看| 国产成人a区在线观看| 成人高潮视频无遮挡免费网站| 婷婷色综合大香蕉| 99久久综合免费| 亚洲av综合色区一区| 亚洲久久久国产精品| 国产在线一区二区三区精| 亚洲久久久国产精品| 成人国产麻豆网| 国产精品久久久久久av不卡| 日韩av在线免费看完整版不卡| 国产欧美日韩一区二区三区在线 | 国产精品国产三级专区第一集| 美女主播在线视频| 大码成人一级视频| 国产成人免费无遮挡视频| 男人和女人高潮做爰伦理| 亚洲精品国产av蜜桃| 夜夜爽夜夜爽视频| 国产探花极品一区二区| h日本视频在线播放| 国产 一区精品| 美女国产视频在线观看| 欧美日韩视频高清一区二区三区二| 成人无遮挡网站| 色婷婷久久久亚洲欧美| 一级av片app| 中文字幕人妻熟人妻熟丝袜美| 赤兔流量卡办理| 亚洲色图av天堂| 久久精品久久久久久噜噜老黄| 久久人人爽人人爽人人片va| 久久精品久久精品一区二区三区| 亚洲va在线va天堂va国产| 青春草国产在线视频| 国产精品不卡视频一区二区| 最近中文字幕2019免费版| 精品人妻视频免费看| 日韩大片免费观看网站| 99九九线精品视频在线观看视频| 亚洲婷婷狠狠爱综合网| 亚洲图色成人| 久久精品国产亚洲av涩爱| 国产免费福利视频在线观看| 国产精品福利在线免费观看| 五月玫瑰六月丁香| 国产精品久久久久久精品电影小说 | av不卡在线播放| 尤物成人国产欧美一区二区三区| 美女脱内裤让男人舔精品视频| 美女国产视频在线观看| 搡老乐熟女国产| 日韩亚洲欧美综合| 熟妇人妻不卡中文字幕| h视频一区二区三区| 日本黄色日本黄色录像| av在线观看视频网站免费| 各种免费的搞黄视频| 午夜免费鲁丝| 亚洲不卡免费看| 大香蕉久久网| 日本wwww免费看| 国产视频首页在线观看| a 毛片基地| 国产精品99久久99久久久不卡 | 久久国内精品自在自线图片| 午夜激情福利司机影院| 国产片特级美女逼逼视频| 久久精品久久久久久久性| 亚洲国产精品999| av播播在线观看一区| 国产综合精华液| 中文字幕制服av| 99国产精品免费福利视频| 欧美亚洲 丝袜 人妻 在线| 久久人人爽av亚洲精品天堂 | 亚洲av中文av极速乱| 色视频在线一区二区三区| 国产精品麻豆人妻色哟哟久久| 亚洲精品亚洲一区二区| 欧美另类一区| 久久国内精品自在自线图片| 日本欧美视频一区| 久久久a久久爽久久v久久| kizo精华| 精品国产露脸久久av麻豆| 天天躁夜夜躁狠狠久久av| 18禁裸乳无遮挡动漫免费视频| 久久99热这里只频精品6学生| 国产片特级美女逼逼视频| 新久久久久国产一级毛片| 熟女电影av网| 日韩中字成人| 久久综合国产亚洲精品| 亚洲av综合色区一区| 久久久欧美国产精品| 国产女主播在线喷水免费视频网站| www.色视频.com| 性色avwww在线观看| 高清毛片免费看| 久久久久久久久久成人| 亚洲欧美精品专区久久| 国产乱人视频| 一级片'在线观看视频| 嫩草影院新地址| 少妇熟女欧美另类| 波野结衣二区三区在线| 视频中文字幕在线观看| 人妻系列 视频| 涩涩av久久男人的天堂| 黄色一级大片看看| 91aial.com中文字幕在线观看| av在线播放精品| 日本黄色片子视频| 一本一本综合久久| tube8黄色片| 国产精品一区二区性色av| 国产亚洲欧美精品永久| 日韩 亚洲 欧美在线| 尾随美女入室| 我的女老师完整版在线观看| 亚洲人与动物交配视频| 黄色日韩在线| 在线观看国产h片| 亚洲国产精品国产精品| av天堂中文字幕网| 精品久久久久久电影网| 深爱激情五月婷婷| 嫩草影院新地址| tube8黄色片| 大又大粗又爽又黄少妇毛片口| av网站免费在线观看视频| 日韩欧美精品免费久久| 极品教师在线视频| 最近最新中文字幕免费大全7| 黄色一级大片看看| 高清日韩中文字幕在线| 国产精品福利在线免费观看| 一本—道久久a久久精品蜜桃钙片| 狂野欧美激情性xxxx在线观看| 久久久久久久久久久丰满| 国产成人精品久久久久久| 少妇的逼好多水| 免费观看无遮挡的男女| 一本一本综合久久| 国产美女午夜福利| 亚洲美女搞黄在线观看| 国产毛片在线视频| 欧美三级亚洲精品| 国产真实伦视频高清在线观看| 高清日韩中文字幕在线| 在线观看免费高清a一片| 亚洲av综合色区一区| h视频一区二区三区| 岛国毛片在线播放| 日韩不卡一区二区三区视频在线| 99久久综合免费| 老熟女久久久| 交换朋友夫妻互换小说| 亚洲av福利一区| 欧美精品亚洲一区二区| kizo精华| 少妇人妻一区二区三区视频| 午夜福利影视在线免费观看| 亚洲精品中文字幕在线视频 | 久久久久国产精品人妻一区二区| 亚洲欧美成人精品一区二区| 这个男人来自地球电影免费观看 | 亚洲第一区二区三区不卡| 大香蕉久久网| 91aial.com中文字幕在线观看| 日日摸夜夜添夜夜爱| 青春草视频在线免费观看| 特大巨黑吊av在线直播| 日韩一区二区视频免费看| 观看av在线不卡| 欧美xxxx黑人xx丫x性爽| 涩涩av久久男人的天堂| 人妻夜夜爽99麻豆av| 色视频在线一区二区三区| 丝袜脚勾引网站| 色5月婷婷丁香| a级毛片免费高清观看在线播放| 毛片女人毛片| 亚洲欧美成人精品一区二区| 久久久久精品久久久久真实原创| 欧美另类一区| 免费看日本二区| 亚洲无线观看免费| 男人和女人高潮做爰伦理| 晚上一个人看的免费电影| 久久国产乱子免费精品| 国内揄拍国产精品人妻在线| 日本wwww免费看| 日韩av免费高清视频| 狂野欧美激情性xxxx在线观看| 一区二区三区四区激情视频| 嫩草影院入口| av线在线观看网站| 女的被弄到高潮叫床怎么办| 成人一区二区视频在线观看| 伦精品一区二区三区| 免费av不卡在线播放| 久久国产亚洲av麻豆专区| 久久久久久久精品精品| 国产在线一区二区三区精| 欧美亚洲 丝袜 人妻 在线| 18禁在线无遮挡免费观看视频| 99久久精品一区二区三区| 久久6这里有精品| 国产一区有黄有色的免费视频| 一区在线观看完整版| 亚洲人与动物交配视频| 亚洲欧美精品自产自拍| 91精品国产国语对白视频| 国产淫语在线视频| 久久久成人免费电影| 男人狂女人下面高潮的视频| 久久精品人妻少妇| 国产精品嫩草影院av在线观看| 国产在线免费精品| 天堂8中文在线网| 久热这里只有精品99| 国产精品国产三级国产av玫瑰| 久久久久久久国产电影| 日日啪夜夜撸| 水蜜桃什么品种好| 男女免费视频国产| 久久午夜福利片| 人妻制服诱惑在线中文字幕| 久久影院123| h视频一区二区三区| 黑人高潮一二区| 精品亚洲乱码少妇综合久久| 欧美日韩精品成人综合77777| 三级国产精品片| av网站免费在线观看视频| 2021少妇久久久久久久久久久| 日本黄色日本黄色录像| 国产视频首页在线观看| 久久这里有精品视频免费| 亚洲av福利一区| 亚洲人成网站高清观看| 少妇丰满av| 嘟嘟电影网在线观看| 国产极品天堂在线| 乱系列少妇在线播放| 成年美女黄网站色视频大全免费 | 成人亚洲精品一区在线观看 | 国产爽快片一区二区三区| 国产午夜精品久久久久久一区二区三区| 亚洲欧美日韩另类电影网站 | 亚洲欧洲日产国产| 久久国产精品男人的天堂亚洲 | 成人漫画全彩无遮挡| 18+在线观看网站| 国产男女内射视频| 国产高清三级在线| 国产一区亚洲一区在线观看| 一级爰片在线观看| 国产成人免费观看mmmm| 97在线人人人人妻| 99久久精品热视频| 在线免费观看不下载黄p国产| 日韩成人伦理影院| 多毛熟女@视频| 国产高清三级在线| 国产欧美另类精品又又久久亚洲欧美| 久久久久久久久久久丰满| 22中文网久久字幕| 国产亚洲91精品色在线| 日韩 亚洲 欧美在线| 成人免费观看视频高清| 国产精品一区二区在线不卡| 小蜜桃在线观看免费完整版高清| 777米奇影视久久| 99热这里只有是精品50| 我的女老师完整版在线观看| 亚洲人与动物交配视频| 亚洲av中文av极速乱| 简卡轻食公司| 一级毛片黄色毛片免费观看视频| 欧美成人一区二区免费高清观看| 黄色日韩在线| 王馨瑶露胸无遮挡在线观看| www.av在线官网国产| 啦啦啦中文免费视频观看日本| 久久鲁丝午夜福利片| 国产av国产精品国产| 久久鲁丝午夜福利片| 国产成人免费无遮挡视频| www.av在线官网国产| 高清不卡的av网站| 国产在线视频一区二区| 一区二区三区四区激情视频| 高清视频免费观看一区二区| 超碰av人人做人人爽久久| 午夜免费观看性视频| 精品一区在线观看国产| 亚洲精品一二三| 成人18禁高潮啪啪吃奶动态图 | 自拍偷自拍亚洲精品老妇| 久久久久性生活片| 亚洲精品自拍成人| a级毛色黄片| 一区在线观看完整版| av线在线观看网站| 国产精品秋霞免费鲁丝片| 99热这里只有精品一区| 日本wwww免费看| 国产成人精品一,二区| 男女边吃奶边做爰视频| 另类亚洲欧美激情| 视频中文字幕在线观看| 伦理电影大哥的女人| 久久99热这里只频精品6学生| 国产成人精品福利久久| 最后的刺客免费高清国语| 身体一侧抽搐| 久久久久国产精品人妻一区二区| 黄片无遮挡物在线观看| 国产亚洲一区二区精品| 日韩伦理黄色片| 国产日韩欧美亚洲二区| 亚洲精品色激情综合| 日韩视频在线欧美| 亚洲精品国产成人久久av| 人妻一区二区av| 午夜福利视频精品| 久久人妻熟女aⅴ| 亚洲人与动物交配视频| 久久国产乱子免费精品| 亚洲美女搞黄在线观看| 亚州av有码| 国产精品国产av在线观看| 亚洲精品国产色婷婷电影| 免费看光身美女| 国产亚洲最大av| 内射极品少妇av片p| 成人无遮挡网站| 大话2 男鬼变身卡| 成年免费大片在线观看| 国产 一区精品| 免费黄网站久久成人精品| 中国三级夫妇交换| 久久久久久人妻| 免费不卡的大黄色大毛片视频在线观看| 99久久精品一区二区三区| 日韩成人av中文字幕在线观看| 三级经典国产精品| 久久久久网色| 亚洲欧美精品专区久久| 99热全是精品| 2022亚洲国产成人精品| 成年av动漫网址| 久久久久视频综合| 少妇的逼水好多| 91久久精品国产一区二区成人| 日韩在线高清观看一区二区三区| 久久99热这里只频精品6学生| 国产av精品麻豆| 丝袜脚勾引网站| 少妇人妻一区二区三区视频|