• <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在线免费观看| 久久婷婷人人爽人人干人人爱| 99久久中文字幕三级久久日本| 色综合站精品国产| 丝袜美腿在线中文| 亚洲精品久久国产高清桃花| 国产视频内射| 亚洲人成伊人成综合网2020| 九九热线精品视视频播放| 国产亚洲精品久久久com| 国产精品一区二区性色av| 精品久久久久久久人妻蜜臀av| 亚洲精品456在线播放app | 天天躁日日操中文字幕| 97热精品久久久久久| 特大巨黑吊av在线直播| 国产探花在线观看一区二区| 免费看光身美女| 久久久久免费精品人妻一区二区| 中文资源天堂在线| 特级一级黄色大片| 神马国产精品三级电影在线观看| 久久婷婷人人爽人人干人人爱| 精品一区二区免费观看| 中亚洲国语对白在线视频| 国产亚洲精品久久久com| 麻豆精品久久久久久蜜桃| 欧美极品一区二区三区四区| 尤物成人国产欧美一区二区三区| 欧美又色又爽又黄视频| 九色成人免费人妻av| 麻豆国产97在线/欧美| 少妇猛男粗大的猛烈进出视频 | 免费在线观看影片大全网站| 一进一出好大好爽视频| av视频在线观看入口| 99热这里只有是精品在线观看| 十八禁国产超污无遮挡网站| 不卡视频在线观看欧美| 国产三级中文精品| 91久久精品国产一区二区三区| 国产亚洲av嫩草精品影院| 极品教师在线视频| 日韩欧美在线乱码| 成人特级黄色片久久久久久久| www日本黄色视频网| or卡值多少钱| 国产视频内射| 韩国av一区二区三区四区| 国产女主播在线喷水免费视频网站 | 日韩亚洲欧美综合| 乱码一卡2卡4卡精品| 欧美不卡视频在线免费观看| 亚洲成人精品中文字幕电影| 午夜亚洲福利在线播放| 国产视频内射| 亚洲无线观看免费| 亚洲自拍偷在线| 亚洲av日韩精品久久久久久密| 日本一二三区视频观看| 99精品在免费线老司机午夜| 亚洲av中文av极速乱 | 韩国av一区二区三区四区| 日韩一本色道免费dvd| 日日撸夜夜添| 成人无遮挡网站| 欧洲精品卡2卡3卡4卡5卡区| 又紧又爽又黄一区二区| 丰满人妻一区二区三区视频av| 啪啪无遮挡十八禁网站| 天堂影院成人在线观看| 一级黄色大片毛片| 久久久久久久久大av| 色精品久久人妻99蜜桃| 级片在线观看| av.在线天堂| 久久中文看片网| 一个人观看的视频www高清免费观看| 成人二区视频| 欧美不卡视频在线免费观看| 国产精品av视频在线免费观看| 制服丝袜大香蕉在线| 国产中年淑女户外野战色| 联通29元200g的流量卡| 亚洲真实伦在线观看| 午夜爱爱视频在线播放| 少妇人妻精品综合一区二区 | 午夜老司机福利剧场| 内地一区二区视频在线| 久久精品91蜜桃| 亚洲av.av天堂| 在线观看午夜福利视频| 亚洲在线自拍视频| 日韩欧美国产一区二区入口| 观看美女的网站| 99热这里只有是精品在线观看| 韩国av一区二区三区四区| 精品人妻1区二区| 成人特级av手机在线观看| 18禁黄网站禁片午夜丰满| 国产av麻豆久久久久久久| 身体一侧抽搐| 十八禁网站免费在线| 日韩亚洲欧美综合| netflix在线观看网站| 亚洲精品亚洲一区二区| 欧美日韩乱码在线| av黄色大香蕉| 中文字幕av成人在线电影| 99在线人妻在线中文字幕| 国产欧美日韩一区二区精品| 国产亚洲精品综合一区在线观看| 亚洲成人精品中文字幕电影| 久久午夜福利片| 中文字幕人妻熟人妻熟丝袜美| 女同久久另类99精品国产91| 99在线人妻在线中文字幕| 麻豆国产av国片精品| 国产欧美日韩精品一区二区| 欧美又色又爽又黄视频| 国产精品女同一区二区软件 | 国产精品久久久久久精品电影| 中国美女看黄片| 黄色丝袜av网址大全| 国产精品98久久久久久宅男小说| 哪里可以看免费的av片| 日本撒尿小便嘘嘘汇集6| 欧美人与善性xxx| 亚洲成人中文字幕在线播放| 欧美在线一区亚洲| 国产精品一区二区免费欧美| 嫩草影视91久久| 亚洲精品粉嫩美女一区| 成年女人毛片免费观看观看9| 日本熟妇午夜| 99热只有精品国产| 麻豆国产av国片精品| 午夜视频国产福利| 精品人妻视频免费看| 男女之事视频高清在线观看| 日韩欧美国产在线观看| 又粗又爽又猛毛片免费看| 听说在线观看完整版免费高清| 国产中年淑女户外野战色| 91狼人影院| 少妇丰满av| 性色avwww在线观看| 人妻丰满熟妇av一区二区三区| 国产精品综合久久久久久久免费| 听说在线观看完整版免费高清| 久久精品国产99精品国产亚洲性色| 国产伦一二天堂av在线观看| 黄片wwwwww| 午夜福利在线观看吧| 久久久久久久亚洲中文字幕| 日本与韩国留学比较| 91在线精品国自产拍蜜月| 人妻久久中文字幕网| 一级黄色大片毛片| 久久久久久大精品| 精品99又大又爽又粗少妇毛片 | 制服丝袜大香蕉在线| 精品欧美国产一区二区三| 国产在线男女| 国产不卡一卡二| 国产精品久久视频播放| 国产探花在线观看一区二区| 国产av在哪里看| 欧美zozozo另类| 老女人水多毛片| 日本免费一区二区三区高清不卡| 俄罗斯特黄特色一大片| 日韩欧美国产在线观看| 国产男靠女视频免费网站| 女人十人毛片免费观看3o分钟| 高清在线国产一区| 在线观看av片永久免费下载| 欧美黑人巨大hd| 九色国产91popny在线| 三级国产精品欧美在线观看| 欧美成人性av电影在线观看| 啦啦啦韩国在线观看视频| 日韩精品青青久久久久久| 国内毛片毛片毛片毛片毛片| 精品人妻一区二区三区麻豆 | 神马国产精品三级电影在线观看| 夜夜爽天天搞| 天堂av国产一区二区熟女人妻| 日日啪夜夜撸| 欧美黑人巨大hd| 日本一本二区三区精品| 日本色播在线视频| 99在线视频只有这里精品首页| 午夜免费男女啪啪视频观看 | 最近视频中文字幕2019在线8| 成人亚洲精品av一区二区| 日本黄色视频三级网站网址| 国产精品精品国产色婷婷| 亚洲不卡免费看| 亚洲中文字幕一区二区三区有码在线看| 女生性感内裤真人,穿戴方法视频| 国产69精品久久久久777片| 国语自产精品视频在线第100页| 亚洲精品色激情综合| 乱系列少妇在线播放| 不卡视频在线观看欧美| 精品久久久久久久久av| 久久精品国产鲁丝片午夜精品 | 免费在线观看成人毛片| av福利片在线观看| 免费大片18禁| 精品国内亚洲2022精品成人| 一区福利在线观看| 日本撒尿小便嘘嘘汇集6| .国产精品久久| 女生性感内裤真人,穿戴方法视频| 国产单亲对白刺激| 国产女主播在线喷水免费视频网站 | 床上黄色一级片| 美女高潮的动态| 女生性感内裤真人,穿戴方法视频| 免费看日本二区| 夜夜看夜夜爽夜夜摸| 欧美一区二区国产精品久久精品| 别揉我奶头 嗯啊视频| 动漫黄色视频在线观看| 亚洲av第一区精品v没综合| 日本在线视频免费播放| 欧美+亚洲+日韩+国产| 久久亚洲真实| 91在线观看av| 亚洲综合色惰| 99热6这里只有精品| 成年版毛片免费区| 亚洲无线在线观看| 精品久久久久久成人av| 成人国产一区最新在线观看| 午夜免费激情av| 国产伦一二天堂av在线观看| 欧美黑人欧美精品刺激| 欧洲精品卡2卡3卡4卡5卡区| 村上凉子中文字幕在线| 国产精品一区www在线观看 | 日本欧美国产在线视频| 少妇高潮的动态图| 琪琪午夜伦伦电影理论片6080| 中文字幕免费在线视频6| 我的老师免费观看完整版| 国产精品女同一区二区软件 | 十八禁网站免费在线| 99热这里只有是精品在线观看| 欧美日韩综合久久久久久 | 国产高清三级在线| 村上凉子中文字幕在线| 能在线免费观看的黄片| 在线观看免费视频日本深夜| 黄色视频,在线免费观看| 天堂影院成人在线观看| 一a级毛片在线观看| 国产精品美女特级片免费视频播放器| 国产白丝娇喘喷水9色精品| 久久久久久国产a免费观看| 欧美另类亚洲清纯唯美| 日韩精品青青久久久久久| 国产精品久久视频播放| 国产精品99久久久久久久久| 丰满乱子伦码专区| 国产精品一区二区免费欧美| 搡女人真爽免费视频火全软件 | 久久草成人影院| 12—13女人毛片做爰片一| 久久6这里有精品| av专区在线播放| 婷婷六月久久综合丁香| 中出人妻视频一区二区| 日韩欧美国产一区二区入口| 最近最新免费中文字幕在线| 国产伦一二天堂av在线观看| 天堂动漫精品| 精品午夜福利视频在线观看一区| 最好的美女福利视频网| 国产麻豆成人av免费视频| 免费av不卡在线播放| 精品午夜福利在线看| 亚洲五月天丁香| 亚洲av电影不卡..在线观看| 直男gayav资源| 国内精品宾馆在线| 国产 一区 欧美 日韩| 国产亚洲精品av在线| 亚洲乱码一区二区免费版| 久久6这里有精品| 老司机午夜福利在线观看视频| 色综合婷婷激情| 麻豆成人av在线观看| 久久久成人免费电影| 欧美激情在线99| 精品久久久久久久久久久久久| av在线蜜桃| 欧美性感艳星| 婷婷丁香在线五月| 欧美日韩瑟瑟在线播放| 精品久久久久久久久亚洲 | 中文字幕人妻熟人妻熟丝袜美| 国语自产精品视频在线第100页| 他把我摸到了高潮在线观看| 日韩高清综合在线| 亚洲在线观看片| 国产爱豆传媒在线观看| 日韩强制内射视频| 一级av片app| 日韩国内少妇激情av| 国产三级在线视频| 黄色视频,在线免费观看| 老司机福利观看| 不卡视频在线观看欧美| 少妇丰满av| 一区二区三区四区激情视频 | 日本一二三区视频观看| 别揉我奶头 嗯啊视频| 少妇熟女aⅴ在线视频| 日韩大尺度精品在线看网址| 国产高潮美女av| 偷拍熟女少妇极品色| 国产69精品久久久久777片| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久伊人网av| 一级黄片播放器| 人人妻人人看人人澡| 一区福利在线观看| 99在线视频只有这里精品首页| 日日摸夜夜添夜夜添av毛片 | 精品久久国产蜜桃| 高清毛片免费观看视频网站| 成熟少妇高潮喷水视频| 黄色一级大片看看| 少妇被粗大猛烈的视频| 69av精品久久久久久| xxxwww97欧美| 国产探花在线观看一区二区| 国产精品久久久久久亚洲av鲁大| 精品久久国产蜜桃| 国产v大片淫在线免费观看| av在线亚洲专区| 热99在线观看视频| 亚洲欧美清纯卡通| 联通29元200g的流量卡| 波多野结衣高清无吗| 免费大片18禁| 免费看日本二区| 国内精品一区二区在线观看| 美女黄网站色视频| 亚洲人成伊人成综合网2020| 免费在线观看影片大全网站| 亚洲四区av| 变态另类成人亚洲欧美熟女| 久久久久性生活片| 久久精品人妻少妇| 97热精品久久久久久| 亚洲av五月六月丁香网| 亚洲欧美日韩卡通动漫| 看免费成人av毛片| 91在线精品国自产拍蜜月| 日日摸夜夜添夜夜添av毛片 | 波多野结衣巨乳人妻| 国产av不卡久久| 十八禁国产超污无遮挡网站| 五月玫瑰六月丁香| 欧美一级a爱片免费观看看| 免费电影在线观看免费观看| 国产aⅴ精品一区二区三区波| 一区福利在线观看| 欧美另类亚洲清纯唯美| 少妇丰满av| 如何舔出高潮| 国产真实伦视频高清在线观看 | 精品国产三级普通话版| 别揉我奶头 嗯啊视频| 国产主播在线观看一区二区| 亚洲图色成人| 男女做爰动态图高潮gif福利片| 中文字幕免费在线视频6| 97人妻精品一区二区三区麻豆| 久久久国产成人精品二区| 午夜影院日韩av| 久久99热6这里只有精品| av中文乱码字幕在线| 热99re8久久精品国产| 久久精品人妻少妇| 97超级碰碰碰精品色视频在线观看| 波多野结衣高清无吗| 美女 人体艺术 gogo| 中文字幕av成人在线电影| 久久久国产成人免费| 亚洲成人久久性| 真人一进一出gif抽搐免费| 国产极品精品免费视频能看的| 99久久精品国产国产毛片| 亚洲电影在线观看av| 国内精品宾馆在线| 12—13女人毛片做爰片一| av女优亚洲男人天堂| 国产高清视频在线播放一区| 国产美女午夜福利| 又粗又爽又猛毛片免费看| av女优亚洲男人天堂| 99久久久亚洲精品蜜臀av| 国产人妻一区二区三区在| 日本-黄色视频高清免费观看| 俺也久久电影网| 日本黄色视频三级网站网址| 日本一二三区视频观看| 国产精品一区二区性色av| 好男人在线观看高清免费视频| 国产女主播在线喷水免费视频网站 | 日本 欧美在线| 国产高清不卡午夜福利| 午夜精品在线福利| 精品人妻偷拍中文字幕| x7x7x7水蜜桃| 搡老岳熟女国产| 校园春色视频在线观看| 久久久久久伊人网av| av在线老鸭窝| 精品久久久久久久久av| 在线观看舔阴道视频| 日本免费a在线| 国产探花极品一区二区| 久久久久久久精品吃奶| av福利片在线观看| 男插女下体视频免费在线播放| 国内精品美女久久久久久| 麻豆国产av国片精品| 男女啪啪激烈高潮av片| 成年免费大片在线观看| 午夜a级毛片| 精品人妻一区二区三区麻豆 | 亚洲熟妇中文字幕五十中出| 少妇的逼水好多| 狂野欧美激情性xxxx在线观看| 在线观看舔阴道视频| 一夜夜www| 免费无遮挡裸体视频| 午夜福利欧美成人| 亚洲第一电影网av| 久久精品影院6| 国产在线男女| 91在线精品国自产拍蜜月| 看十八女毛片水多多多| 国产精华一区二区三区| 久久热精品热| 国产精品久久久久久久久免| www.色视频.com| 99热这里只有是精品50| 亚洲精品国产成人久久av| 欧美日韩国产亚洲二区| 色播亚洲综合网| 日韩中文字幕欧美一区二区| 亚洲成人久久爱视频| 亚洲精品在线观看二区| 色av中文字幕| 国产成年人精品一区二区| 国产欧美日韩精品一区二区| 国产伦人伦偷精品视频| 国产国拍精品亚洲av在线观看| 伊人久久国产一区二区| 国产精品免费大片| 亚洲性久久影院| 国产成人免费无遮挡视频| 日本免费在线观看一区| 日韩 亚洲 欧美在线| 国产真实伦视频高清在线观看| 狠狠精品人妻久久久久久综合| 丰满人妻一区二区三区视频av| 免费人成在线观看视频色| 噜噜噜噜噜久久久久久91| 国产欧美亚洲国产| 一本久久精品| 这个男人来自地球电影免费观看 | 男人爽女人下面视频在线观看| 午夜日本视频在线| 国产v大片淫在线免费观看| 国产成人91sexporn| 亚洲三级黄色毛片| 99九九线精品视频在线观看视频| 免费av不卡在线播放| 久久久色成人| 日本av免费视频播放| 老熟女久久久| 亚洲内射少妇av| 内射极品少妇av片p| 亚洲婷婷狠狠爱综合网| 国产永久视频网站| 美女视频免费永久观看网站| 97超视频在线观看视频| 国产又色又爽无遮挡免| av卡一久久| 身体一侧抽搐| 婷婷色av中文字幕| 高清日韩中文字幕在线| 国产成人午夜福利电影在线观看| 日日摸夜夜添夜夜爱| 国产成人精品久久久久久| 秋霞在线观看毛片| 成人国产麻豆网| 男女无遮挡免费网站观看| 国产精品精品国产色婷婷| 亚洲成色77777| 国产淫片久久久久久久久| 亚洲欧美日韩无卡精品| 精品亚洲成国产av| 成人18禁高潮啪啪吃奶动态图 | 纯流量卡能插随身wifi吗| 一本—道久久a久久精品蜜桃钙片| 联通29元200g的流量卡| 人人妻人人添人人爽欧美一区卜 | 久久97久久精品| 欧美日韩一区二区视频在线观看视频在线| 久久精品国产亚洲网站| 日本wwww免费看| 99久久精品一区二区三区| 亚洲精品久久午夜乱码| av.在线天堂| 国产亚洲一区二区精品| 亚洲综合色惰| 另类亚洲欧美激情| 精品久久国产蜜桃| 亚洲欧美日韩东京热| 日本一二三区视频观看| 免费高清在线观看视频在线观看| 日日啪夜夜爽| 亚洲精品国产色婷婷电影| 少妇高潮的动态图| 水蜜桃什么品种好| 最新中文字幕久久久久| 免费人成在线观看视频色| 久久久久久久精品精品| 肉色欧美久久久久久久蜜桃| 青春草国产在线视频| 精品国产乱码久久久久久小说| 精品一区二区三卡| 日本av手机在线免费观看| 色网站视频免费| 久久精品国产亚洲网站| 两个人的视频大全免费| 国产片特级美女逼逼视频| 欧美精品一区二区大全| 国产精品嫩草影院av在线观看| 天堂中文最新版在线下载| 91精品伊人久久大香线蕉| 国产精品.久久久| 国产精品三级大全| 亚洲国产精品国产精品| 一区二区av电影网| 九九爱精品视频在线观看| 国产精品秋霞免费鲁丝片| 极品少妇高潮喷水抽搐| 亚洲av男天堂| 又黄又爽又刺激的免费视频.| 精品午夜福利在线看| 亚洲三级黄色毛片| 色哟哟·www| 久久亚洲国产成人精品v| 欧美日韩视频高清一区二区三区二| 一个人看视频在线观看www免费| 啦啦啦中文免费视频观看日本| 少妇高潮的动态图| 精品久久久久久久末码| 国产熟女欧美一区二区| 中文乱码字字幕精品一区二区三区| av网站免费在线观看视频| 水蜜桃什么品种好| 日本与韩国留学比较| 国产v大片淫在线免费观看| 欧美bdsm另类| 国产在线一区二区三区精| 亚洲自偷自拍三级| 秋霞伦理黄片| 成人18禁高潮啪啪吃奶动态图 | 国产精品三级大全| 亚洲性久久影院| 91精品伊人久久大香线蕉| 边亲边吃奶的免费视频| 熟妇人妻不卡中文字幕| 久久ye,这里只有精品| 高清午夜精品一区二区三区| 国产乱来视频区| 一本—道久久a久久精品蜜桃钙片| 一个人看视频在线观看www免费| 在线 av 中文字幕| 看非洲黑人一级黄片| 久热这里只有精品99| 亚洲国产最新在线播放| 国产精品成人在线| 男女无遮挡免费网站观看| 美女福利国产在线 | 好男人视频免费观看在线| 日日啪夜夜撸| 最近最新中文字幕大全电影3| 欧美性感艳星| 大话2 男鬼变身卡| 丝瓜视频免费看黄片| 久久久久人妻精品一区果冻| 永久网站在线| 狠狠精品人妻久久久久久综合| 成人影院久久| 欧美三级亚洲精品| 成年av动漫网址| 成人影院久久| 三级国产精品欧美在线观看| 久久久久久人妻| 国产精品女同一区二区软件| 身体一侧抽搐| 精品亚洲乱码少妇综合久久|