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

    板狀粘接結(jié)構(gòu)中對稱和反對稱縱波與界面的相互作用

    2017-04-11 06:31:34丁俊才何存富
    振動與沖擊 2017年6期
    關(guān)鍵詞:反射系數(shù)橫波縱波

    丁俊才, 吳 斌, 何存富

    (北京工業(yè)大學(xué) 機械工程與應(yīng)用電子技術(shù)學(xué)院,北京 100124)

    板狀粘接結(jié)構(gòu)中對稱和反對稱縱波與界面的相互作用

    丁俊才, 吳 斌, 何存富

    (北京工業(yè)大學(xué) 機械工程與應(yīng)用電子技術(shù)學(xué)院,北京 100124)

    為了激發(fā)粘接結(jié)構(gòu)中的導(dǎo)波或界面波,通常需要將聲波從兩半無限介質(zhì)同相位或反相位地同時入射多層系統(tǒng)。針對此問題,基于矩陣方法,推導(dǎo)了界面處于理想連接的情況下,對稱或反對稱縱波從上下半無限空間入射時,三層板狀粘接結(jié)構(gòu)中縱、橫波的反射與透射系數(shù)表達式。分析了入射角度、粘接層厚度以及基體材料等對聲波反射(或透射)特性的影響。結(jié)果表明,對稱或反對稱縱波垂直入射時不發(fā)生波型轉(zhuǎn)換。粘接結(jié)構(gòu)中聲波的反射(或透射)特性與入射角度、頻率以及粘接層厚度等參數(shù)密切相關(guān)。在相同的粘接層厚度(或頻率)范圍內(nèi),隨著聲波頻率(或粘接層厚度)的增加,諧振頻率曲線向低頻漂移。該方法可作為粘接結(jié)構(gòu)中體波或?qū)Р▊鞑ヌ匦匝芯康闹匾碚摶A(chǔ)。

    粘接結(jié)構(gòu);矩陣方法;對稱和反對稱縱波;反射和透射系數(shù);諧振

    隨著現(xiàn)代飛機高性能發(fā)展的需要和新技術(shù)、新材料的應(yīng)用,粘接結(jié)構(gòu)因具有比強度、比模量高,以及密封、減振等優(yōu)越性能,已成為飛機機體的重要結(jié)構(gòu)形式[1],其中金屬板-板粘接結(jié)構(gòu)是飛機粘接結(jié)構(gòu)的主體[2-4]。由于粘接界面脫粘或粘接強度不足而導(dǎo)致的災(zāi)難性事件時有發(fā)生[5-6],因此,發(fā)展有效的金屬板-板結(jié)構(gòu)粘接質(zhì)量的檢測方法不僅具有重要的學(xué)術(shù)價值,而且具有廣泛的工程應(yīng)用前景[7]。

    關(guān)于多層板狀粘接結(jié)構(gòu)中超聲波傳播特性的研究一直受到重視。在利用超聲體波對粘接結(jié)構(gòu)的研究上,BREKHOVSKIKH等[8]利用阻抗和聲壓連續(xù)的方法推導(dǎo)了多層均勻介質(zhì)中體波的反射和透射系數(shù)表達式。李明軒等[9]建立了共振匹配理論,利用縱橫波綜合判斷方法從實驗上對粘接界面的強度進行了定量評價。基于傳遞矩陣方法,THOMSON[10]利用斜入射聲波的反射和透射系數(shù)從理論上分析了粘接層的特性。ROKHLIN等[11-12]將粘接界面作為粘彈性層,對有限厚度界面層處于兩半無限各向同性介質(zhì)之間的多層粘接結(jié)構(gòu)中的體波傳播特性進行了研究,較好地區(qū)分出了理想連接和脫粘界面。王耀俊[13]利用傳遞矩陣方法研究了具有剛性(即理想)連接和滑移界面的層狀固體介質(zhì)的聲反射和透射特性,推導(dǎo)了縱波或橫波由單側(cè)基體入射時的反射和透射系數(shù)解析表達式。在彈簧模型法的基礎(chǔ)上,QIU等[14-15]采用垂直入射和斜入射的方法研究了弱粘接結(jié)構(gòu)中的體波傳播特性,并利用COMSOL進行了數(shù)值仿真。吳斌等[16]研究了浸水斜入射條件下粘接結(jié)構(gòu)的縱波透射特性,得出了完整的縱波透射系數(shù)表達式。安志武等[17]在忽略粘接層厚度的前提下結(jié)合彈簧模型法利用“右峰”漂移特性定性地分析了鋁-粘接層-鋁弱粘接結(jié)構(gòu)中超聲波的低頻反射特性。

    利用超聲導(dǎo)波對粘接結(jié)構(gòu)的研究方面,BANERJEE[18]從理論上分析了粘接結(jié)構(gòu)中的導(dǎo)波傳播模式,研究發(fā)現(xiàn),垂直激發(fā)時響應(yīng)信號被一階反對稱模態(tài)所控制。LE-CROM等[19]通過數(shù)值模擬得到了鋁-復(fù)合板粘接結(jié)構(gòu)的SH波頻散曲線。CASTAINGS[20]利用有限元和實驗的方法研究了處于不同界面形式的鋁-樹脂玻璃-鋁搭接結(jié)構(gòu)中SH波的傳播。YEW等[21]利用SH波評估了粘接結(jié)構(gòu)粘接層的質(zhì)量,研究表明,粘接結(jié)構(gòu)中SH波二階模態(tài)的截止頻率主要依賴于粘接層的厚度和力學(xué)性能。張海燕等[22]采用全局矩陣法研究了層狀各向異性復(fù)合板中Lamb波的傳播特性,給出了不同層數(shù)的復(fù)合板中Lamb波的相速度頻散曲線以及模態(tài)A0和S0沿板厚方向的應(yīng)力分布。寧偉等[23]基于彈簧模型法研究了三層粘接結(jié)構(gòu)中Lamb波兩個最低階模態(tài)A0和S0所產(chǎn)生的質(zhì)點位移沿厚度方向的分布。

    由以上分析可知,許多學(xué)者采用不同的方法或界面簡化模型結(jié)合超聲體波或?qū)Р▽φ辰咏Y(jié)構(gòu)開展了研究并取得了相應(yīng)的成果。實際上,多數(shù)針對粘接結(jié)構(gòu)的研究往往指的是對于界面特性的研究。相對于體波來講,利用超聲導(dǎo)波或界面波對粘接結(jié)構(gòu)進行無損檢測可獲取更多關(guān)于界面的局部特征信息。若想激發(fā)結(jié)構(gòu)中的導(dǎo)波或界面波,通常需要將聲波從兩半無限介質(zhì)(即粘接結(jié)構(gòu)的上下基體)同相位或反相位地同時入射進多層粘接系統(tǒng)。但在以往的文獻報道中,大多采用的是超聲波由單側(cè)基體入射的方法,未見有關(guān)于對稱(反相位)或反對稱(同相位)模式的超聲波由兩側(cè)基體同時入射粘接結(jié)構(gòu)的研究。與超聲波由單側(cè)基體入射檢測粘接結(jié)構(gòu)不同的是,若采用超聲波由兩側(cè)基體同時入射,則可以提取到更純凈、更單一的導(dǎo)波模態(tài),因此很有必要研究超聲波由兩側(cè)基體同時入射時粘接結(jié)構(gòu)中聲波的反射與透射情況。

    針對上述未研究的問題,本文基于矩陣方法,首先從理論上推導(dǎo)出粘接結(jié)構(gòu)界面處于理想連接的情況下,對稱或反對稱縱波同時從上下半無限固體空間入射時,三層板狀粘接結(jié)構(gòu)中縱、橫波的反射與透射系數(shù)表達式;接著以實際工業(yè)生產(chǎn)中常用的粘接結(jié)構(gòu)為例,利用數(shù)值模擬的方法計算結(jié)構(gòu)中縱、橫波的反射與透射系數(shù)曲線;最后分析聲波入射角度、粘接層厚度等參數(shù)的變化對聲波反射與透射特性的影響。

    1 理論研究

    圖1所示為板狀粘接結(jié)構(gòu)幾何模型,將笛卡爾坐標系的x軸置于粘接結(jié)構(gòu)的下界面(界面2),板厚沿z方向。半無限固體介質(zhì)2和0分別為粘接結(jié)構(gòu)的上、下基體,介質(zhì)1是厚度為h的粘接層,且介質(zhì)0、1和2均為各向同性彈性固體介質(zhì)。

    圖1 體波斜入射條件下粘接結(jié)構(gòu)的聲學(xué)模型Fig.1 Acoustic model of the adhesive structure underthe condition of body waves oblique incidence

    圖1中φi2、φr2和φi0、φr0分別為基體2和0中入射、反射縱波的位移勢幅值;ψi2、ψr2和ψi0、ψr0分別為基體2和0中入射、反射橫波(此處及下文提及的橫波均指SV橫波)的位移勢幅值。θl2和θl0分別為基體2和0中的縱波入射或反射角;θt2和θt0分別為基體2和0中的橫波入射或反射角。若令角頻率同為ω的簡諧平面縱波和橫波同時從上、下基體入射到多層粘接結(jié)構(gòu),這些波沿著x方向應(yīng)具有相同的波矢量分量δ。容易知道,當聲波遇到兩側(cè)聲阻抗不同的結(jié)合界面時,會發(fā)生反射與折射,形成反射縱、橫波與折射縱、橫波。根據(jù)Snell定律,所有的反射波與折射波的波矢量沿著x方向的分量也為δ。若忽略波傳播的非線性效應(yīng),則在同一介質(zhì)中入射波和反射波可以線性疊加。另外需要說明的是,本文的研究暫不考慮聲波在結(jié)構(gòu)中傳播時的衰減。

    假設(shè)聲波的入射平面為x-z平面,那么在基體0中縱波和橫波的總位移勢Φ0和Ψ0可分別寫為[24-26]

    (1)

    基于式(1),可將基體0中的聲場所產(chǎn)生的位移與應(yīng)力分量寫成

    (2)

    式中:ux和σzx、uz和σzz分別為沿著x,z方向的位移和應(yīng)力分量,即切向與法向的位移和應(yīng)力分量;λ0和μ0為介質(zhì)0的拉梅常數(shù)。

    在z=0處,將式(1)和(2)相結(jié)合可以得到基體0中的聲場和界面2上的力學(xué)參量之間的關(guān)系,將其寫成矩陣的形式如下

    (3)

    采用與式(1)~(3)類似的分析,在z=h處可以寫出界面1上的質(zhì)點位移和應(yīng)力與基體2中的聲場之間的關(guān)系

    (4)

    (5)

    式中:矩陣[B2]=[A2]-1為一與基體2材料特性相關(guān)的四階矩陣,其值見附錄2。

    由于粘接結(jié)構(gòu)各界面均為理想連接,因此界面上的質(zhì)點位移和應(yīng)力連續(xù)。在此種條件下,我們可將上下基體中的聲場聯(lián)系起來并寫成式(6)的形式

    (6)

    式中:矩陣[H]=[B2]·[E1]·[A0],這里矩陣[E1]為與粘接層1材料特性相關(guān)的矩陣,其各個元素的值見附錄3。本文著重研究對稱或反對稱縱波入射粘接結(jié)構(gòu),因此式(6)中ψi2與ψi0的值為0。將ψi2=0和ψi0=0代入式(6)并重新整理可以得出式(7)

    (7)

    式中:Hgm(g,m=1,2,3,4)為矩陣[H]的元素。按照式(7)可求出不同模式的平面縱波入射時,所產(chǎn)生的縱、橫波的反射與透射系數(shù)。一般地,大多研究均采用縱波或橫波從上或下半無限空間入射多層結(jié)構(gòu),這里若令式(7)中的φi0或φi2=0,則可以分別得出縱波從上或下半無限介質(zhì)入射時的反射與透射系數(shù)。作為例子,本文假定任意形式的頻率均為ω的兩列簡諧平面縱波同相位(反對稱)或反相位(對稱)地同時入射多層粘接結(jié)構(gòu),研究不同形式聲波與界面的相互作用。

    針對反對稱縱波同時從上下基體入射多層結(jié)構(gòu),此時令φi0=φi2并代入式(7),可以得出式(8)

    (8)

    對式(8)進行繁瑣的運算,可以得到反對稱縱波入射時結(jié)構(gòu)中所產(chǎn)生的縱、橫波的反射與透射系數(shù)表達式

    (9)

    這里,

    F1=H23H31-H21H33+H21(H12H33-H13H32)+

    H22(H31H13-H11H33)+H23(H11H32-H31H12)

    F2=H43H31-H41H33+H41(H12H33-H13H32)+

    H42(H31H13-H11H33)+H43(H11H32-H31H12)

    式中,R表示反射系數(shù),T表示透射系數(shù),且第一個下標表示入射波型,第二個下標表示反射或透射波型。例如Rlt指采用縱波反對稱入射所產(chǎn)生的橫波的反射系數(shù)。

    針對對稱縱波同時從上下基體入射多層結(jié)構(gòu),此時令φi0=-φi2并代入式(7),則可以得出:

    (10)

    同樣地,對式(10)進行計算,可以得到對稱縱波入射時產(chǎn)生的縱、橫波的反射與透射系數(shù)表達式

    (11)

    這里,

    F3=H23H31-H21H33+H21(H13H32-H12H33)+

    H22(H11H33-H31H13)+H23(H31H12-H11H32)

    F4=H43H31-H41H33+H41(H13H32-H12H33)+

    H42(H11H33-H31H13)+H43(H31H12-H11H32)

    式中,R′、T′的含義與前文的R、T相同,均指反射、透射系數(shù),但前者為對稱縱波入射下的反射、透射系數(shù),后者為反對稱縱波入射下的反射、透射系數(shù)。實際上,式(9)和(11)中的反射和透射系數(shù)一般為復(fù)數(shù),表明反射波和透射波相對于入射波存在相位差。

    需要指出,針對對稱或反對稱縱波入射粘接結(jié)構(gòu),由于上、下基體皆有聲波進入,因此在稱呼上并無反射與透射之分??蓪⒒w2中的Ψr2、φr2稱作反射波,基體0中的Ψr0、φr0稱之為透射波;也可以將基體0中的Ψr0、φr0稱作反射波,基體2中的Ψr2、φr2稱為透射波。本文為便于讀者理解以及避免混淆,暫且將基體2中的Ψr2、φr2分別稱之為反射橫、縱波;將基體0中的Ψr0、φr0分別稱之為透射橫、縱波。物理意義上講,盡管式(9)中的Rll和Tll、Rlt和Tlt的表達式的形式有所不同,但它們的數(shù)值計算結(jié)果應(yīng)該相等,式(11)也應(yīng)如此。

    2 數(shù)值算例及結(jié)果分析

    基于上述模型,利用MATLAB軟件編寫相應(yīng)的程序計算對稱或反對稱縱波入射下的縱、橫波的反射、透射系數(shù),研究縱、橫波的反射、透射特性。文中選取常見的鋁、鋼或銅作為粘接結(jié)構(gòu)的基體材料,環(huán)氧樹脂為粘接劑,它們的材料參數(shù)見表1。表1中ρ、Cl和Ct分別為材料的密度、縱波和橫波速度。

    表1 材料特性Tab.1 Material properties

    2.1 入射角度對聲波反射或透射特性的影響

    圖2和3計算了對稱或反對稱縱波以不同角度入射的情況下,鋁-環(huán)氧樹脂-鋁粘接結(jié)構(gòu)中縱、橫波的反射與透射系數(shù)。圖中縱波的入射角度分別取0°、20°、40°和60°,粘接層的厚度為0.3 mm。

    對比圖2(a)和3(a)、2(b)和3(b)可以發(fā)現(xiàn),當入射角度為20°、40°或60°,對稱和反對稱縱波入射下的縱、橫波的反射(或透射)系數(shù)有較為明顯的區(qū)別。表明聲波傾斜入射時,對稱和反對稱縱波與同一粘接結(jié)構(gòu)的相互作用不同。當入射角為0°(垂直入射)時,無論縱波對稱或反對稱入射,結(jié)構(gòu)中橫波的反射(或透射)系數(shù)的幅值始終為0,而縱波的反射(或透射)系數(shù)的幅值始終為1,這是由于對稱或反對稱縱波垂直入射時不發(fā)生波型轉(zhuǎn)換,使得結(jié)構(gòu)中僅有縱波存在??v波傾斜入射會發(fā)生波型轉(zhuǎn)換形成新的縱、橫波,所有縱、橫波的反射(或透射)系數(shù)曲線隨著頻率的增加出現(xiàn)了諧振,這與多層結(jié)構(gòu)中聲傳播理論相符合。在諧振頻率點處,橫波的反射(或透射)系數(shù)呈現(xiàn)極小值,而縱波的反射(或透射)系數(shù)呈現(xiàn)極大值,這說明當聲波的頻率處于諧振點附近時,粘接結(jié)構(gòu)中大多以縱波的形式存在而形成的橫波較少。此外,仔細觀察圖2和3可知,無論縱波對稱或反對稱入射,也不論縱波的入射角度或頻率為多大,粘接結(jié)構(gòu)中縱波或橫波的反射與透射系數(shù)曲線均完全重合,這證實了前文的猜測:盡管式(9)、(11)中縱波或橫波的反射與透射系數(shù)表達式的形式有所不同,但它們的數(shù)值計算結(jié)果相等,同時也間接反映出本文所推導(dǎo)公式的正確性。

    (a) 橫波

    (b) 縱波

    (a) 橫波

    (b) 縱波

    2.2 粘接層厚度對聲波反射或透射特性的影響

    圖4和5計算了鋁-環(huán)氧樹脂-鋁粘接結(jié)構(gòu)在不同粘接層厚度下的縱、橫波的反射與透射系數(shù),其中粘接層的厚度分別為0.01 mm、0.1 mm和1 mm,縱波的入射角取20°。

    (a) 橫波

    (b) 縱波

    (a) 橫波

    由圖可見,即便粘接層的厚度有微小的變化也會對聲波的反射(或透射)特性造成較大的影響。無論縱波對稱或是反對稱入射,橫波的反射(或透射)系數(shù)的極小值均接近于0,極大值皆不超過0.3??v波的反射(或透射)系數(shù)的極大值均接近于1且所有頻率點處的幅值不低于0.9,表明當入射角為20°時不論粘接層的厚度是多少,粘接結(jié)構(gòu)中縱波的能量始終遠大于橫波。此外還可以發(fā)現(xiàn),在相同的頻率范圍(即0~10 MHz)內(nèi),隨著粘接層厚度的增加,縱、橫波的反射(或透射)系數(shù)曲線逐漸出現(xiàn)了諧振現(xiàn)象,并且粘接層的厚度越大,曲線諧振點(見圖4(b))的數(shù)目越多。這是由于隨著粘接層厚度的增加,聲波在粘接層上下界面(見圖1)來回反射、干涉或疊加的次數(shù)在增多,因此產(chǎn)生了更多的諧振頻率點。但是所形成的多階諧振并無明顯的變化規(guī)律,也即沒有明顯的諧振周期。另外,與2.1節(jié)類似,無論縱波對稱還是反對稱入射,也無論粘接層的厚度是多少,縱波或橫波的反射與透射系數(shù)曲線依然完全重合。

    (b) 縱波

    2.3 基體材料對聲波反射特性的影響

    本節(jié)僅計算縱波與橫波的反射系數(shù),透射系數(shù)圖不再繪制。圖6~圖8所示為對稱或反對稱縱波入射時,不同形式粘接結(jié)構(gòu)中的體波反射特性,其中粘接層的厚度固定在0.3 mm,聲波入射角度仍取20°。

    (a) 橫波

    (b) 縱波

    (a) 橫波

    (b) 縱波

    (a) 對稱縱波入射

    (b) 反對稱縱波入射

    由圖可知,若粘接結(jié)構(gòu)的基體為不同材料,縱波或橫波的反射系數(shù)存在比較明顯的區(qū)別,說明基體材料對聲波的反射特性有一定影響。圖6(a)與7(a)中,無論對稱還是反對稱縱波入射,四種不同形式粘接結(jié)構(gòu)中的橫波反射系數(shù)曲線整體走勢大致相同。鋁-環(huán)氧樹脂-鋼、鋁-環(huán)氧樹脂-銅以及銅-環(huán)氧樹脂-鋁三種形式粘接結(jié)構(gòu)中的橫波反射系數(shù)差別不大,但是鋼-環(huán)氧樹脂-鋁中的橫波反射系數(shù)與上述三種粘接結(jié)構(gòu)差別明顯。對于圖6(b)與7(b),在對稱或反對稱縱波入射的情況下,鋁-環(huán)氧樹脂-鋼和鋁-環(huán)氧樹脂-銅、鋼-環(huán)氧樹脂-鋁和銅-環(huán)氧樹脂-鋁粘接結(jié)構(gòu)中的縱波反射系數(shù)差別較小,但是鋁-環(huán)氧樹脂-鋼(或鋁-環(huán)氧樹脂-銅)和鋼-環(huán)氧樹脂-鋁以及銅-環(huán)氧樹脂-鋁中的縱波反射系數(shù)差別非常明顯。鋁-環(huán)氧樹脂-鋼與鋁-環(huán)氧樹脂-銅結(jié)構(gòu)中的縱波反射系數(shù)介于1~2.5之間,遠大于鋼-環(huán)氧樹脂-鋁與銅-環(huán)氧樹脂-鋁中的縱波反射系數(shù)。此外,比較圖6(a)和6(b)以及7(a)和7(b)可以發(fā)現(xiàn),即使粘接結(jié)構(gòu)和入射聲波的類型固定,結(jié)構(gòu)中縱波與橫波的反射系數(shù)之間的差別也異常顯著。

    與圖6和7不同的是,圖8中粘接結(jié)構(gòu)的上下基體均為同種介質(zhì)。分析可知,對于縱波對稱入射,鋁-環(huán)氧樹脂-鋁、鋼-環(huán)氧樹脂-鋼以及銅-環(huán)氧樹脂-銅中縱波反射系數(shù)的差別不大。類似地,對于縱波反對稱入射,除了9.2 MHz處,三種形式粘接結(jié)構(gòu)中縱波反射系數(shù)的差別亦不大。比較圖8(a)和6(b)可以得出,同為對稱縱波入射的情況下,鋁-環(huán)氧樹脂-鋁和鋁-環(huán)氧樹脂-鋼、鋁-環(huán)氧樹脂-銅(或鋼-環(huán)氧樹脂-鋁、銅-環(huán)氧樹脂-鋁)中縱波反射系數(shù)的差別明顯。這表明,若保持粘接結(jié)構(gòu)的粘接層和上(或下)基體的材料不變,僅改變下(或上)基體,結(jié)構(gòu)中聲傳播模式同樣會發(fā)生較大的變化。圖8(b)和7(b)的分析方法與此類似,這里不做過多描述。

    上述現(xiàn)象表明,當其它參數(shù)固定,粘接結(jié)構(gòu)中聲波的傳播特性與粘接層兩側(cè)基體材料的特征參數(shù)密切相關(guān)。不僅基體材料的種類對聲傳播特性的影響較大,基體所處的位置對聲傳播的影響同樣不可忽略。這對于檢測人員利用對稱或反對稱縱波檢測粘接結(jié)構(gòu)具有較好的參考價值。

    2.4 頻率、粘接層厚度與橫波反射特性的三維關(guān)系研究

    如前文所言,對于粘接結(jié)構(gòu)的檢測,入射角度、頻率和粘接層厚度一直是重要的參數(shù)。因此,有必要研究對稱或反對稱縱波以不同角度入射的情況下,頻率、粘接層厚度與聲波反射(或透射)特性的關(guān)系。圖9~14為對稱或反對稱縱波分別以5°、15°或30°入射鋁-環(huán)氧樹脂-鋁粘接結(jié)構(gòu)時,縱波頻率、粘接層厚度與橫波反射系數(shù)之間的關(guān)系。圖中顏色的深淺代表反射系數(shù)幅值的大小,類似曲線的部分代表反射系數(shù)的極小值點,其余紅色部分代表反射系數(shù)的極大值點。分析可以發(fā)現(xiàn),隨著聲波頻率以及粘接層厚度的改變,橫波反射系數(shù)均出現(xiàn)了多階諧振。

    (a)

    (b)

    (a)

    (b)

    (a)

    (b)

    (a)

    (b)

    (a)

    (b)

    (a)

    (b)

    針對對稱縱波入射的情形(圖9~圖11),隨著入射角度的增加,雖然反射系數(shù)的極大值在逐漸增大,但是諧振頻率曲線的總體數(shù)目并未隨角度的變化而改變。除此以外,還可以得到另外一個重要現(xiàn)象:在特定角度(例如5°)下,當聲波的頻率范圍固定不變,隨著粘接層厚度的增大,諧振頻率曲線的數(shù)目在增多,如圖9(b)中水平橢圓圈所示,當粘接層的厚度為0.2 mm,橢圓圈內(nèi)僅存在兩條諧振頻率曲線,當粘接層厚度增大到0.6 mm時,諧振頻率曲線的數(shù)目增加到七條,反過來也是如此??蓪⑦@種現(xiàn)象理解為:在相同的粘接層厚度(或頻率)范圍內(nèi),隨著聲波頻率(或粘接層厚度)的增加,諧振頻率曲線向低頻產(chǎn)生了漂移。反對稱縱波入射下(圖12~圖14)的分析方法與圖9~圖11類似,這里不再贅述。

    3 結(jié) 論

    假定粘接界面為理想連接的情況下,本文基于矩陣方法建立了對稱或反對稱縱波穿過三層板狀粘接結(jié)構(gòu)的聲反射和透射模型,推導(dǎo)了結(jié)構(gòu)中體波的反射和透射系數(shù)表達式。基于上述數(shù)值模擬結(jié)果,得到以下結(jié)論:

    (1) 與體波由單側(cè)基體入射的情形類似,對稱或反對稱縱波垂直入射時,經(jīng)過界面不發(fā)生波型轉(zhuǎn)換,但是傾斜入射波型轉(zhuǎn)換會發(fā)生。不同形式的粘接結(jié)構(gòu)會引起不同的聲反射(或透射)特性。對于特定的粘接結(jié)構(gòu),超聲波在其中的傳播模式主要依賴于入射角度、頻率等參數(shù),即使這些參數(shù)有微小的變化也會對聲波的反射(或透射)特性造成較大的影響。因此,若采用對稱或反對稱縱波檢測粘接結(jié)構(gòu),合理的選取這些參數(shù)尤為重要。

    (2) 對稱或反對稱縱波入射粘接結(jié)構(gòu)時縱、橫波的反射(或透射)系數(shù)曲線會產(chǎn)生多階諧振。在相同的粘接層厚度(或頻率)范圍內(nèi),隨著聲波頻率(或粘接層厚度)的增加,諧振頻率曲線將向低頻漂移。

    本文的理論分析模型對于界面為理想連接且基體和粘接層均為各向同性彈性固體介質(zhì)的粘接結(jié)構(gòu)具有普適性。對于粘接結(jié)構(gòu)的研究,除了體波以外還有不少學(xué)者采用導(dǎo)波(例如Lamb波)或界面波進行檢測,雖然本文主要針對對稱或反對稱縱波入射下的體波反射(或透射)系數(shù)展開討論,但是所提方法同樣可以作為粘接結(jié)構(gòu)中導(dǎo)波傳播特性研究的重要理論基礎(chǔ),比如若令式(9)或(11)中的分母為0,就可以得到結(jié)構(gòu)中沿x方向傳播的導(dǎo)波模式的頻散方程。需要指出的是,鑒于實驗條件的挑戰(zhàn)性,目前尚未對理論計算結(jié)果開展相關(guān)的仿真或?qū)嶒炦^程,但所述方法的有效性與所推公式的正確性已在文中得以證實。

    [ 1 ] SERNA-MORENO M C, LPEZ CELA J J, MARTINEZ VICENTE J L, et al. Adhesively bonded joints as a dissipative energy mechanism under impact loading [J]. Applied Mathematical Modelling, 2015, 39(12): 3496-3505.

    [ 2 ] BENYAHIA F, ALBEDAH A, BOUIADJRA B B. Analysis of the adhesive damage for different patch shapes in bonded composite repair of aircraft structures [J]. Materials & Design, 2014, 54: 18-24.

    [ 3 ] BENYAHIA F, BOUANANI M F, ALBEDAH A, et al. Effect of water absorption on the adhesive damage in bonded composite repair of aircraft structures [J]. Materials & Design, 2014, 57: 435-441.

    [ 4 ] ZHOU Z G, MA B Q, JIANG J T, et al. Application of wavelet filtering and barker-coded pulse compression hybrid method to air-coupled ultrasonic testing [J]. Nondestructive Testing and Evaluation, 2014, 29(4): 297-314.

    [ 5 ] SU W G, ZOU L, MU Z T, et al. Stress analysis of cracked metallic aircraft structure adhesively repaired with composite patch [J]. Lecture Notes in Electrical Engineering, 2014, 297(75): 369-377.

    [ 6 ] PANTELAKIS S, TSERPES K I. Adhesive bonding of composite aircraft structures: challenges and recent developments [J]. Science China-Physics Mechanics & Astronomy, 2014, 57(1): 2-11.

    [ 7 ] 許巍, 陳力, 張錢城, 等. 粘接界面力學(xué)行為及其表征 [J]. 中國科學(xué): 技術(shù)科學(xué), 2012, 42(12): 1361-1376. XU Wei, CHEN Li, ZHANG Qiancheng, et al. Mechanical behavior and its representation of the bonding interface [J]. Scientia Sinica Technologica, 2012, 42(12): 1361-1376.

    [ 8 ] BREKHOVSKIKH L M, GODIN O A. Acoustics of layered media [M]. Berlin: Springer Verlag, 1990.

    [ 9 ] 李明軒, 王小民, 安志武. 粘接界面特性的超聲檢測與評價 [J]. 應(yīng)用聲學(xué), 2013, 32(3): 190-198. LI Mingxuan, WANG Xiaomin, AN Zhiwu. Ultrasonic testing and evaluation of bonding interface properties [J]. Applied Acoustics, 2013, 32(3): 190-198.

    [10] THOMSON W T. Transmission of elastic waves through a stratified solid medium [J]. Journal of Applied Physics, 1950, 21(2): 88-93.

    [11] ROKHLIN S I, MAROM D. Study of adhesive bonds using low-frequency obliquely incident ultrasonic waves [J]. Journal of the Acoustical Society of America, 1986, 80(2): 585-590.

    [12] ROKHLIN S I, WANG Y J. Analysis of boundary conditions for elastic wave interaction with an interface between two solids [J]. Journal of the Acoustical Society of America, 1991, 89(2): 503-515.

    [13] 王耀俊. 具有剛性聯(lián)接界面和滑移界面的層狀固體媒質(zhì)的聲反射 [J]. 聲學(xué)學(xué)報, 1992, 17(2): 81-92. WANG Yaojun. Sound reflection from layered solid medium with rigid and slip interfaces [J]. Acta Acustica, 1992, 17(2): 81-92.

    [14] QIU Z G, WU B, HE C F. Research on adhesive layer depth and frequency in weak interface of bonded structures [J]. Non-Destructive Testing and Condition Monitoring, 2011, 53(5): 302-306.

    [15] QIU Z G, WU B, HE C F. Ultrasonic beam steering using neumann boundary condition in multiplysics [J]. Acta Mechanica Sinica, 2012, 28(1): 146-150.

    [16] 吳斌, 張婧, 邱兆國, 等. 浸水斜入射條件下黏接結(jié)構(gòu)的透射特性研究 [J]. 機械工程學(xué)報, 2012, 49(10): 45-52. WU Bin, ZHANG Jing, QIU Zhaoguo, et al. Research on the transmission characteristics of underwater bonding structure in oblique incidence [J]. Journal of Mechanical Engineering, 2012, 49(10):45-52.

    [17] 安志武, 王小民, 毛捷, 等. 等厚鋁板粘接件超聲反射頻譜的低頻特性 [J]. 應(yīng)用聲學(xué), 2009, 28(3): 190-194. AN Zhiwu, WANG Xiaomin, MAO Jie, et al. Low-frequency characteristic of ultrasound reflection spectrum from adhesively bonded aluminum plates [J]. Applied Acoustics, 2009, 28(3): 190-194.

    [18] BANERJEE S. Theoretical modeling of guided wave propagation in a sandwich plate subjected to transient surface excitations [J]. International Journal of Solids and Structures, 2012, 49(23/24): 3233-3241.

    [19] LE-CROM B, CASTAINGS M. Shear horizontal guided wave modes to infer the shear stiffness of adhesive bond layers [J]. Journal of the Acoustical Society of America, 2010, 127(4): 2220-2230.

    [20] CASTAINGS M. SH ultrasonic guided waves for the evaluation of interfacial adhesion [J]. Ultrasonics, 2014, 54(7): 1760-1775.

    [21] YEW C H, WENG X W. Using ultrasonic SH waves to estimate the quality of adhesive bonds in plate structures [J]. Journal of the Acoustical Society of America, 1985, 77(5): 1813-1823.

    [22] 張海燕, 劉鎮(zhèn)清, 呂東輝. 全局矩陣法及其在層狀各向異性復(fù)合板中Lamb波傳播特性研究中的應(yīng)用 [J]. 復(fù)合材料學(xué)報, 2004, 21(2): 111-116. ZHANG Haiyan, LIU Zhenqing, Lü Donghui. Global matrix and its application to study on Lamb wave propagation in layered anisotropic composite plates [J]. Acta Materiae Compositae Sinica, 2004, 21(2): 111-116.

    [23] 寧偉, 王耀俊. 三層固體復(fù)合媒質(zhì)中Lamb的傳播 [J]. 聲學(xué)學(xué)報, 1996, 21(4): 414-420. NING Wei, WANG Yaojun. Lamb wave propagation in three-layered composites [J]. Acta Acustica, 1996, 21(4): 414-420.

    [24] ACHENBACH J D. Acoustic fields and waves in solids [M]. Beijing: Science Press, 1982.

    [25] WANG Y J. Acoustic wave reflection and transmission on multilayered media [J]. Journal of Nanjing University (Natural Sciences Edition), 1993, 29(1): 49-62.

    [26] ROSE J L. Ultrasonic waves in solid media [M]. Beijing:Science Press, 2004.

    附錄1

    附錄2

    這里ρ0、Ct0和Cl0分別為基體0的密度、橫波波速與縱波波速,ρ2、Ct2和Cl2分別為基體2的密度、橫波波速與縱波波速,ω為角頻率。

    附錄3

    這里,ρ1、Cl1和Ct1為粘接層的密度、縱波速度和橫波速度;θl1和θt1分別為縱波和橫波在粘接層中的傳播角。

    Interaction of the symmetric and anti-symmetric longitudinal waveswith the interfaces in plate-like adhesive structure

    DING Juncai,WU Bin,HE Cunfu

    (College of Mechanical Engineering and Applied Electronics Technology,Beijing University of Technology,Beijing 100124, China)

    To stimulate the guided or interface waves in adhesive structures, in general sound waves should be injected into a multilayer system from two infinite media in the same or opposed phase at the same time. Aiming at this issue, based on the matrix method, the expressions of reflection and transmission coefficients for longitudinal and transverse waves in three-tier plate-like adhesive structures were deduced when on interfaces there is perfect connection and the symmetric or anti-symmetric longitudinal waves are injected from upper and lower infinite spaces. The influences of the incident angle, thickness of adhesive layer and substrate materials etc on the characteristics of acoustic wave reflection (or transmission) were analyzed. The results show that when the symmetric or anti-symmetric longitudinal waves are vertically incident, the conversions of wave modes do not occur. The reflection (or transmission) characteristics of acoustic waves in adhesive structures are closely related to the incident angle, frequency, thickness of adhesive layer and other parameters. In the same range of adhesive layer thickness (or frequency), with the increase of the sound frequency (or adhesive layer thickness), the curves of resonant frequency tend to drift towards the low frequency. The proposed method can be used as an important theoretical basis for the research on transmission characteristics of body or guided waves in adhesive structures.

    adhesive structure; matrix method; symmetric and anti-symmetric longitudinal waves; reflection and transmission coefficients; resonance

    國家自然科學(xué)基金(11132002;11372016;51475012;51235001)

    2015-10-10 修改稿收到日期: 2016-01-24

    丁俊才 男,博士生,1984年生

    吳斌 男,教授,博士生導(dǎo)師,1962年生 E-mail: wb@bjut.edu.cn

    O327

    A

    10.13465/j.cnki.jvs.2017.06.005

    猜你喜歡
    反射系數(shù)橫波縱波
    橫波技術(shù)在工程物探中的應(yīng)用分析
    多道隨機稀疏反射系數(shù)反演
    石油物探(2020年6期)2020-11-25 02:38:46
    黃257井區(qū)疊前縱波方位各向異性裂縫分布預(yù)測
    球面波PP反射系數(shù)的頻變特征研究
    變截面階梯桿中的縱波傳播特性實驗
    揚眉一顧,妖嬈橫波處
    橫波一顧,傲殺人間萬戶侯
    火花(2015年1期)2015-02-27 07:40:24
    沙質(zhì)沉積物反射系數(shù)的寬帶測量方法
    橫波淺層地震在城市勘探中的應(yīng)用
    基于反射系數(shù)的波導(dǎo)結(jié)構(gòu)不連續(xù)位置識別
    欧美日本中文国产一区发布| 在线视频色国产色| 久久精品91蜜桃| 一本大道久久a久久精品| www.熟女人妻精品国产| 真人一进一出gif抽搐免费| av在线天堂中文字幕 | 免费少妇av软件| 亚洲三区欧美一区| 免费在线观看视频国产中文字幕亚洲| 亚洲精品国产色婷婷电影| 操出白浆在线播放| 国产黄色免费在线视频| 久久中文字幕一级| 99在线人妻在线中文字幕| 一区二区日韩欧美中文字幕| www.自偷自拍.com| 久久精品91蜜桃| 老熟妇乱子伦视频在线观看| 99久久99久久久精品蜜桃| xxx96com| 欧洲精品卡2卡3卡4卡5卡区| aaaaa片日本免费| 真人一进一出gif抽搐免费| 一级片免费观看大全| 俄罗斯特黄特色一大片| 国产精品一区二区三区四区久久 | 免费观看人在逋| 亚洲精品美女久久av网站| 男女做爰动态图高潮gif福利片 | 中文字幕高清在线视频| 国产精品 欧美亚洲| 成人国语在线视频| 免费在线观看黄色视频的| www.www免费av| 国产成人精品无人区| 久久精品人人爽人人爽视色| 深夜精品福利| 久久精品国产亚洲av香蕉五月| 免费看a级黄色片| 黄色 视频免费看| 亚洲五月色婷婷综合| 高清欧美精品videossex| 免费在线观看日本一区| 国产片内射在线| 色综合欧美亚洲国产小说| 视频区图区小说| 一区在线观看完整版| 久久久久国内视频| 欧美日韩亚洲综合一区二区三区_| 在线天堂中文资源库| 久久精品国产清高在天天线| av在线播放免费不卡| 亚洲美女黄片视频| 国产精品影院久久| 三上悠亚av全集在线观看| 亚洲 欧美 日韩 在线 免费| 水蜜桃什么品种好| 两性午夜刺激爽爽歪歪视频在线观看 | 精品国产乱码久久久久久男人| 美女高潮到喷水免费观看| 深夜精品福利| 久久中文看片网| 另类亚洲欧美激情| 亚洲片人在线观看| 在线观看一区二区三区激情| 韩国精品一区二区三区| 国产精品综合久久久久久久免费 | 国内久久婷婷六月综合欲色啪| 亚洲一区高清亚洲精品| 老司机午夜十八禁免费视频| 日本一区二区免费在线视频| 桃红色精品国产亚洲av| 亚洲国产欧美一区二区综合| 亚洲aⅴ乱码一区二区在线播放 | 亚洲七黄色美女视频| 欧美日韩亚洲综合一区二区三区_| 淫秽高清视频在线观看| av天堂久久9| av片东京热男人的天堂| 午夜久久久在线观看| 精品无人区乱码1区二区| 午夜精品久久久久久毛片777| 久久九九热精品免费| 男人的好看免费观看在线视频 | 一个人观看的视频www高清免费观看 | 90打野战视频偷拍视频| 成年人免费黄色播放视频| 91精品三级在线观看| 在线免费观看的www视频| 一进一出好大好爽视频| 国产蜜桃级精品一区二区三区| 国产不卡一卡二| 午夜a级毛片| av欧美777| 一区福利在线观看| 国产不卡一卡二| 一区福利在线观看| 久久亚洲精品不卡| 日本 av在线| 亚洲精品久久成人aⅴ小说| 国产精品香港三级国产av潘金莲| 啦啦啦免费观看视频1| 精品一区二区三卡| 日韩精品免费视频一区二区三区| 99re在线观看精品视频| 日韩 欧美 亚洲 中文字幕| 日韩 欧美 亚洲 中文字幕| av欧美777| 国产蜜桃级精品一区二区三区| 少妇被粗大的猛进出69影院| 少妇被粗大的猛进出69影院| 91麻豆精品激情在线观看国产 | 人妻久久中文字幕网| 精品久久蜜臀av无| 亚洲第一欧美日韩一区二区三区| 精品午夜福利视频在线观看一区| 欧美一级毛片孕妇| 人妻久久中文字幕网| 欧洲精品卡2卡3卡4卡5卡区| 久久久久精品国产欧美久久久| 亚洲精品国产色婷婷电影| xxx96com| 日日夜夜操网爽| 国产亚洲精品综合一区在线观看 | 一本大道久久a久久精品| netflix在线观看网站| 一本综合久久免费| 国产精品久久久久成人av| 女性被躁到高潮视频| a级毛片黄视频| 一级a爱视频在线免费观看| 午夜两性在线视频| 久久久久久免费高清国产稀缺| 每晚都被弄得嗷嗷叫到高潮| 日韩免费av在线播放| 99国产极品粉嫩在线观看| 久久国产乱子伦精品免费另类| 国产av在哪里看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品久久成人aⅴ小说| 国产91精品成人一区二区三区| 天堂动漫精品| 亚洲人成77777在线视频| 久久人妻熟女aⅴ| 一个人观看的视频www高清免费观看 | 1024视频免费在线观看| 国产色视频综合| 黄色片一级片一级黄色片| 久久久久久亚洲精品国产蜜桃av| 色综合欧美亚洲国产小说| 正在播放国产对白刺激| 亚洲 欧美一区二区三区| 韩国av一区二区三区四区| 国产精品国产高清国产av| 一边摸一边做爽爽视频免费| 亚洲成人久久性| av天堂在线播放| 国产真人三级小视频在线观看| 男女做爰动态图高潮gif福利片 | 日韩中文字幕欧美一区二区| 欧美日韩亚洲高清精品| 久久久久亚洲av毛片大全| 啦啦啦免费观看视频1| 国产一卡二卡三卡精品| 欧美激情极品国产一区二区三区| 叶爱在线成人免费视频播放| 欧美在线一区亚洲| 嫁个100分男人电影在线观看| 一进一出抽搐动态| cao死你这个sao货| 亚洲一码二码三码区别大吗| 视频区图区小说| 99久久精品国产亚洲精品| 免费久久久久久久精品成人欧美视频| 黑人猛操日本美女一级片| 人妻丰满熟妇av一区二区三区| 韩国av一区二区三区四区| 9热在线视频观看99| 亚洲 欧美 日韩 在线 免费| 国产深夜福利视频在线观看| 国产xxxxx性猛交| 精品乱码久久久久久99久播| 女生性感内裤真人,穿戴方法视频| 99国产综合亚洲精品| 夜夜爽天天搞| 亚洲五月婷婷丁香| 国产国语露脸激情在线看| 欧美另类亚洲清纯唯美| 怎么达到女性高潮| 亚洲专区字幕在线| 亚洲熟女毛片儿| 天天躁夜夜躁狠狠躁躁| 村上凉子中文字幕在线| 亚洲视频免费观看视频| 国产一区二区在线av高清观看| 中文字幕人妻丝袜制服| 日韩高清综合在线| 日本wwww免费看| 法律面前人人平等表现在哪些方面| 一级毛片女人18水好多| 老司机深夜福利视频在线观看| 国产亚洲精品久久久久久毛片| 黑人巨大精品欧美一区二区蜜桃| 国产精品香港三级国产av潘金莲| 老汉色av国产亚洲站长工具| 18禁美女被吸乳视频| 制服诱惑二区| 欧美中文日本在线观看视频| 欧美+亚洲+日韩+国产| av欧美777| 久久久国产精品麻豆| 麻豆av在线久日| 日韩有码中文字幕| 久久伊人香网站| 男男h啪啪无遮挡| 黄色成人免费大全| 日韩欧美国产一区二区入口| 日韩欧美一区二区三区在线观看| videosex国产| 1024视频免费在线观看| 免费在线观看亚洲国产| 久久伊人香网站| 欧美日韩国产mv在线观看视频| 国产亚洲av高清不卡| 欧美日本亚洲视频在线播放| 欧美黑人欧美精品刺激| 男女之事视频高清在线观看| 国产欧美日韩一区二区精品| 久久草成人影院| 亚洲 欧美一区二区三区| 欧美日韩av久久| 嫩草影视91久久| 午夜福利一区二区在线看| 女性生殖器流出的白浆| 色精品久久人妻99蜜桃| 精品一区二区三区四区五区乱码| 多毛熟女@视频| 午夜91福利影院| 久久精品91蜜桃| 精品一区二区三卡| 成人免费观看视频高清| 日韩高清综合在线| 中出人妻视频一区二区| 久久 成人 亚洲| 一边摸一边做爽爽视频免费| 精品久久久久久,| 国产黄a三级三级三级人| 国产精品 国内视频| 99久久精品国产亚洲精品| 欧美久久黑人一区二区| 久久久久国产精品人妻aⅴ院| 午夜视频精品福利| 久久精品国产亚洲av香蕉五月| 99精国产麻豆久久婷婷| 国产欧美日韩一区二区精品| 国产成人啪精品午夜网站| 久久精品国产综合久久久| 超碰97精品在线观看| 男人舔女人下体高潮全视频| 免费高清在线观看日韩| 青草久久国产| 欧美中文综合在线视频| 欧美乱码精品一区二区三区| 老汉色∧v一级毛片| 男女下面进入的视频免费午夜 | 老司机午夜十八禁免费视频| 91老司机精品| 精品福利观看| 亚洲成人精品中文字幕电影 | 大香蕉久久成人网| 欧美日韩亚洲国产一区二区在线观看| 丝袜美足系列| 999久久久精品免费观看国产| 一区福利在线观看| 日本免费a在线| 国产国语露脸激情在线看| 国产日韩一区二区三区精品不卡| 激情在线观看视频在线高清| av天堂在线播放| 亚洲精品国产精品久久久不卡| 嫁个100分男人电影在线观看| 久久久国产一区二区| 国产成人影院久久av| 国产成人啪精品午夜网站| 在线观看免费高清a一片| 亚洲成人免费电影在线观看| 一a级毛片在线观看| 香蕉丝袜av| 精品久久久久久电影网| 91大片在线观看| 精品福利观看| 久久国产精品男人的天堂亚洲| 国产一区二区三区在线臀色熟女 | 欧美不卡视频在线免费观看 | 波多野结衣av一区二区av| 天堂√8在线中文| 男人的好看免费观看在线视频 | 男女高潮啪啪啪动态图| 免费日韩欧美在线观看| 亚洲精品一区av在线观看| 五月开心婷婷网| 成人黄色视频免费在线看| 丰满迷人的少妇在线观看| 欧美日韩av久久| 99re在线观看精品视频| 国产伦人伦偷精品视频| 日本vs欧美在线观看视频| 成人精品一区二区免费| 亚洲成人久久性| 91精品三级在线观看| 黄片大片在线免费观看| 亚洲精品国产色婷婷电影| 欧美精品亚洲一区二区| 精品一区二区三卡| 久久精品国产亚洲av高清一级| 成人国语在线视频| 精品久久久精品久久久| 国产深夜福利视频在线观看| 一二三四社区在线视频社区8| 高清毛片免费观看视频网站 | 老熟妇乱子伦视频在线观看| 啦啦啦免费观看视频1| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品乱码一区二三区的特点 | 国产深夜福利视频在线观看| 19禁男女啪啪无遮挡网站| 久久香蕉激情| 91老司机精品| 国产一区二区三区视频了| 天堂动漫精品| 成人黄色视频免费在线看| 黄色a级毛片大全视频| 色哟哟哟哟哟哟| 国产精品一区二区三区四区久久 | 精品久久久久久久久久免费视频 | 国产主播在线观看一区二区| 人成视频在线观看免费观看| 成人亚洲精品av一区二区 | 成在线人永久免费视频| 精品久久久久久电影网| 在线观看舔阴道视频| 日本免费一区二区三区高清不卡 | 欧美黄色片欧美黄色片| 青草久久国产| 高潮久久久久久久久久久不卡| 欧美日韩黄片免| 成年人免费黄色播放视频| 日本一区二区免费在线视频| 国产精品野战在线观看 | 日韩大码丰满熟妇| 国产高清国产精品国产三级| 黑人巨大精品欧美一区二区mp4| 宅男免费午夜| 两性夫妻黄色片| 午夜免费鲁丝| 日本免费一区二区三区高清不卡 | 99国产精品免费福利视频| 成人手机av| 国产激情欧美一区二区| 好男人电影高清在线观看| 成人亚洲精品一区在线观看| 性少妇av在线| 欧美 亚洲 国产 日韩一| 亚洲国产看品久久| 国产精品 国内视频| 久久影院123| 久久婷婷成人综合色麻豆| 美女高潮到喷水免费观看| x7x7x7水蜜桃| 免费在线观看视频国产中文字幕亚洲| 国产又爽黄色视频| 亚洲国产欧美一区二区综合| 天堂影院成人在线观看| 亚洲精品国产区一区二| 五月开心婷婷网| 夜夜夜夜夜久久久久| 岛国在线观看网站| 日韩欧美一区二区三区在线观看| 91国产中文字幕| 免费一级毛片在线播放高清视频 | 日本免费一区二区三区高清不卡 | 国产精品日韩av在线免费观看 | 精品无人区乱码1区二区| av电影中文网址| 亚洲片人在线观看| 亚洲五月婷婷丁香| 日日爽夜夜爽网站| 最好的美女福利视频网| 99国产精品一区二区蜜桃av| 免费av中文字幕在线| 999久久久精品免费观看国产| 国产午夜精品久久久久久| 亚洲成a人片在线一区二区| 国产亚洲欧美98| 精品久久久精品久久久| 亚洲国产看品久久| 久99久视频精品免费| 老司机亚洲免费影院| 精品国产美女av久久久久小说| 超色免费av| 亚洲精品av麻豆狂野| 久久国产乱子伦精品免费另类| 老鸭窝网址在线观看| 淫妇啪啪啪对白视频| 成人国产一区最新在线观看| 免费女性裸体啪啪无遮挡网站| svipshipincom国产片| 国产成人欧美| 国内毛片毛片毛片毛片毛片| 免费久久久久久久精品成人欧美视频| 一边摸一边做爽爽视频免费| 欧美激情极品国产一区二区三区| 国产97色在线日韩免费| 国产精品98久久久久久宅男小说| 日韩欧美三级三区| aaaaa片日本免费| 国产精品一区二区精品视频观看| 中文字幕人妻丝袜一区二区| x7x7x7水蜜桃| 水蜜桃什么品种好| 免费高清视频大片| 亚洲免费av在线视频| 一级a爱视频在线免费观看| 国产一区二区三区综合在线观看| 操美女的视频在线观看| 777久久人妻少妇嫩草av网站| 91精品三级在线观看| 桃色一区二区三区在线观看| 亚洲专区国产一区二区| 国产亚洲欧美在线一区二区| www国产在线视频色| 精品国产超薄肉色丝袜足j| 免费高清在线观看日韩| 日本三级黄在线观看| 国产精品综合久久久久久久免费 | 最好的美女福利视频网| 中国美女看黄片| 色综合婷婷激情| 男女做爰动态图高潮gif福利片 | 亚洲久久久国产精品| 久久久久久亚洲精品国产蜜桃av| 麻豆成人av在线观看| videosex国产| svipshipincom国产片| 黄色 视频免费看| 国产精品日韩av在线免费观看 | 中亚洲国语对白在线视频| 黄色怎么调成土黄色| 91成人精品电影| 免费看十八禁软件| 亚洲国产欧美网| 亚洲激情在线av| 99久久精品国产亚洲精品| 又黄又粗又硬又大视频| 久久精品国产清高在天天线| 欧美日韩福利视频一区二区| 精品久久久精品久久久| 一区福利在线观看| 99香蕉大伊视频| 亚洲va日本ⅴa欧美va伊人久久| 黄片小视频在线播放| x7x7x7水蜜桃| 久热爱精品视频在线9| 黄色视频不卡| 久久精品国产99精品国产亚洲性色 | 身体一侧抽搐| 免费观看人在逋| xxx96com| 久久久国产欧美日韩av| 这个男人来自地球电影免费观看| 亚洲欧洲精品一区二区精品久久久| 三上悠亚av全集在线观看| 日日干狠狠操夜夜爽| 一级毛片女人18水好多| 亚洲欧美日韩无卡精品| 久久99一区二区三区| 精品久久久久久久久久免费视频 | 日本 av在线| 国产色视频综合| 91在线观看av| 91麻豆精品激情在线观看国产 | 国产精品亚洲一级av第二区| 午夜福利免费观看在线| 国产成人啪精品午夜网站| 91字幕亚洲| 久久久久国内视频| 欧美乱色亚洲激情| 成年人免费黄色播放视频| 多毛熟女@视频| av中文乱码字幕在线| 中文欧美无线码| 日韩国内少妇激情av| 日韩三级视频一区二区三区| 多毛熟女@视频| 91成人精品电影| 午夜精品久久久久久毛片777| 国产精品99久久99久久久不卡| 嫩草影院精品99| 精品久久久久久久久久免费视频 | 99精国产麻豆久久婷婷| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产看品久久| 日本撒尿小便嘘嘘汇集6| 丝袜在线中文字幕| 国产精品久久视频播放| 一本大道久久a久久精品| 国产黄a三级三级三级人| 亚洲男人的天堂狠狠| 国产乱人伦免费视频| 亚洲av日韩精品久久久久久密| 免费久久久久久久精品成人欧美视频| 日韩欧美国产一区二区入口| 欧美成狂野欧美在线观看| 国产精品1区2区在线观看.| 国产亚洲精品久久久久5区| 在线观看免费高清a一片| 亚洲av美国av| 久久伊人香网站| 巨乳人妻的诱惑在线观看| 99在线视频只有这里精品首页| 欧美日韩精品网址| 欧美日韩乱码在线| 男男h啪啪无遮挡| 亚洲熟妇中文字幕五十中出 | 精品久久久久久,| 两人在一起打扑克的视频| 日韩欧美国产一区二区入口| 久99久视频精品免费| 亚洲精品中文字幕在线视频| 又大又爽又粗| 亚洲va日本ⅴa欧美va伊人久久| videosex国产| 脱女人内裤的视频| 国产一卡二卡三卡精品| 久久久国产成人精品二区 | 亚洲精品美女久久av网站| 高清在线国产一区| 国产一区在线观看成人免费| 桃色一区二区三区在线观看| 悠悠久久av| 91国产中文字幕| 欧美精品亚洲一区二区| 亚洲狠狠婷婷综合久久图片| 免费观看精品视频网站| 欧美午夜高清在线| 午夜成年电影在线免费观看| 亚洲精品一区av在线观看| 女警被强在线播放| 久久久久久免费高清国产稀缺| 丝袜人妻中文字幕| 国产乱人伦免费视频| 高清黄色对白视频在线免费看| 一本大道久久a久久精品| 亚洲人成77777在线视频| 亚洲色图综合在线观看| 亚洲,欧美精品.| 国产精品一区二区在线不卡| 精品人妻在线不人妻| 久久香蕉激情| 热re99久久精品国产66热6| 日韩大码丰满熟妇| 侵犯人妻中文字幕一二三四区| 一个人观看的视频www高清免费观看 | 成熟少妇高潮喷水视频| 亚洲第一av免费看| 老司机午夜福利在线观看视频| 久久久国产欧美日韩av| 两人在一起打扑克的视频| 午夜免费成人在线视频| 三级毛片av免费| 日韩精品青青久久久久久| 在线观看免费午夜福利视频| 12—13女人毛片做爰片一| 男人操女人黄网站| 丝袜在线中文字幕| 一区在线观看完整版| 久久婷婷成人综合色麻豆| 国产av又大| 女人被躁到高潮嗷嗷叫费观| 一区在线观看完整版| 免费在线观看影片大全网站| 亚洲色图 男人天堂 中文字幕| 国产视频一区二区在线看| 久久精品91无色码中文字幕| 亚洲欧美精品综合一区二区三区| 母亲3免费完整高清在线观看| 成年版毛片免费区| 亚洲国产欧美一区二区综合| a级毛片黄视频| 人成视频在线观看免费观看| 激情在线观看视频在线高清| 一个人观看的视频www高清免费观看 | 老熟妇乱子伦视频在线观看| 国产成人精品无人区| 国产乱人伦免费视频| a级毛片在线看网站| 亚洲欧美日韩另类电影网站| 99riav亚洲国产免费| 中文亚洲av片在线观看爽| 黄片大片在线免费观看| 999久久久精品免费观看国产| 制服诱惑二区| 精品电影一区二区在线| 精品一区二区三区四区五区乱码| av天堂在线播放| 欧美中文日本在线观看视频| 国产xxxxx性猛交| 多毛熟女@视频| 91麻豆精品激情在线观看国产 | 日韩精品免费视频一区二区三区| 免费久久久久久久精品成人欧美视频| 久久精品aⅴ一区二区三区四区| 国产片内射在线| 亚洲久久久国产精品|