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

    金屬薄板缺陷的超聲多途Lamb波檢測

    2016-09-07 05:50:27郭杜斌沈希忠
    計算機測量與控制 2016年1期
    關(guān)鍵詞:延遲時間圓孔薄板

    郭杜斌,沈希忠,王 磊

    (上海應(yīng)用技術(shù)學(xué)院 電氣與電子工程學(xué)院,上?!?01418)

    ?

    金屬薄板缺陷的超聲多途Lamb波檢測

    郭杜斌,沈希忠,王磊

    (上海應(yīng)用技術(shù)學(xué)院 電氣與電子工程學(xué)院,上海201418)

    金屬薄板在動車等現(xiàn)代化設(shè)計中應(yīng)用廣泛,其是否存在缺陷等檢測具有很重要的意義;為此,利用超聲多途提供的額外數(shù)據(jù)來檢測金屬薄板的缺陷,提出相應(yīng)的算法,多途路徑在直接反射的路徑上是可以預(yù)測和識別的,直接反射路徑是相對比較容易識別的;同時利用ANSYS有限元建模不同的薄板缺陷模型,在不同的傳感器位置、缺陷位置、缺陷大小情況下分別進行仿真;通過ANSYS仿真和實驗檢測都得出了良好的結(jié)果,誤差在可控制范圍之內(nèi),證明了超聲多途應(yīng)用與金屬薄板Lamb波檢測的是可行的;仿真和實驗表明超聲多途通過實際和多途虛擬擴展的傳感器數(shù)據(jù)的分析證明了超聲多途能夠?qū)Ρ“迦毕葸M行定位并檢測其大小;提出的超聲多途檢測方法,通過單次測量即可檢測出薄板缺陷的位置和大小,具有良好的檢測性能。

    超聲波;薄板;多途;ANSYS有限元仿真;缺陷;尺寸

    0 引言

    材料缺陷尺寸的檢測在超聲無損評價(NDE)中是個很重要的問題[1-2]。如果在鋁合金配件表面或內(nèi)部存在有缺陷(尤其是裂紋,容易在應(yīng)力作用下導(dǎo)致缺陷進一步增大),在外界載荷作用下材料上分布的集中應(yīng)力將會促使板中缺陷迅速擴大,進而影響鋁合金材料的使用壽命,給社會經(jīng)濟工作帶來安全隱患。超聲波檢測具有對人體無輻射、操作簡單、檢測效率高等優(yōu)點,在鋁合金材料的缺陷檢測中,由于鋁合金材料屬于小顆粒晶體,超聲波在傳播過程中的能量衰減小,對于較厚的鋁合金試件超聲波也能進行有效的檢測,超聲波的上述特性,有利于對檢測結(jié)果的分析,從而促進了超聲波在鋁合金材料探傷領(lǐng)域的推廣應(yīng)用。檢測缺陷尺寸的一種直接方法是從不同的兩個側(cè)面進行B型掃描成像[3]。在金屬薄板的探傷中,Lamb波因其自身在薄板檢測中的優(yōu)勢而被廣泛應(yīng)用,特別是在薄板表面的缺陷探傷[4-6]。

    國外對于Lamb波應(yīng)用于板結(jié)構(gòu)的檢測由來已久,美國的F.A.Firestone等人于20世紀50年代首次將Lamb波用于薄板材料結(jié)構(gòu)的探傷檢測,R.J.Dewhurst[5]等利用激光發(fā)射Lamb波來對金屬薄板的厚度進行檢測,D.A.Hutchins[6]采用激光產(chǎn)生的Lamb波對鋁板和金屬玻璃材料板進行檢測,并研究了Lamb波在鋁和金屬玻璃材料板中傳播波形的頻散特性。J.R.Rose[7]一直在從事研究導(dǎo)波的特性,進行了大量的實驗驗證,并取得了豐碩的研究成果。隨著超聲檢測技術(shù)的不斷發(fā)展,各界學(xué)者對于Lamb的研究也越來越多,并開始引入一些先進的算法對裂紋進行識別,W.Ostachowicz與P.Kudela[8]利用Lamb在鋁合金薄板的傳播特性,提出了一種金屬薄板缺陷定位的算法,仿真和試驗結(jié)果表明該算法能夠較好地對缺陷進行定位。近幾年,基于時間反轉(zhuǎn)的超聲導(dǎo)波[9]以及超聲多途[10]檢測也成為了Lamb波檢測的研究熱點。國內(nèi)的一些研究單位和大學(xué)對Lamb波檢測展開了積極的研究工作。例如南京航空航天大學(xué)[11]在智能材料和結(jié)構(gòu)航空科技重點實驗室開展了對材料結(jié)構(gòu)健康監(jiān)控技術(shù)的研究,并取得了一定的成效;北京工業(yè)大學(xué)[12-13]利用導(dǎo)波在管道、板結(jié)構(gòu)檢測方面開展了大量工作;上海交通大學(xué)[14]、同濟大學(xué)聲學(xué)研究所[15]、中國科學(xué)院聲學(xué)所[16]等也都對板結(jié)構(gòu)的Lamb波檢測進行了大量的理論研究和實驗驗證。

    本文提出了一種排成一行的傳感器傳感器陣列的超聲收發(fā)測量裝置來檢測一個孤立缺陷,不同的是[7],我們對多途路徑進行了分析并結(jié)合了不同的路徑,利用ANSYS有限元仿真對不同模型進行仿真計算,通過實際和多途虛擬擴展的傳感器數(shù)據(jù)的分析證明了超聲多途能夠?qū)Ρ“迦毕葸M行定位并檢測其大小,最終用實驗對其進行了驗證。

    1 超聲多途

    超聲多途在文獻[3]中有過較為詳細的應(yīng)用,該文作者主要將其應(yīng)用于鋁合金塊的缺陷檢測,并成功對缺陷進行定位并檢測出其大小。本文主要利用超聲多途作金屬鋁合金薄板的Lamb波檢測。

    如圖1為一個超聲多途DRP,其超聲信號r(t)建模為

    (1)

    式中,s(t)表示超聲波回波脈沖波,αi和τi分別表示第i次路徑的反射率和延遲時間。

    圖1 標(biāo)準的超聲檢測系統(tǒng)中的直接反射(DRP)場景

    圖2表示多途MP-W路徑的多途場景。MP-W是指Lamb波在薄板中傳播多途路徑中檢測缺陷底部的路徑,它的傳播路徑類似英文字母’W’。

    圖2 多途MP-W的底部反射場景

    假設(shè)缺陷的形狀為類圓形的凸面體,其中心點的坐標(biāo)為fp(xf,p,yf,p,zf,p),第m個發(fā)射傳感器位于St,m(xt,m,yt,m,zt,m),第n個接收傳感器位于Sr,n(xr,n,yr,n,zr,n)。與DRP(直達反射路徑)相對應(yīng)的第n個接收傳感器的信號表示為:

    (2)

    式中,α0,mn(ft)表示缺陷的DRP(直接反射路徑)反射率,τ0,mn(ft)表示信號經(jīng)過傳播一周后的延時時間。

    假設(shè)一塊均勻介質(zhì)的材料中具有恒定的超聲傳播速度vg,與缺陷定點ft相對應(yīng)的延遲時間可以由如下的計算得到

    (3)

    a表示薄板的寬,b表示薄板的長,以薄板的正表面左下角端點為原點建立坐標(biāo)系,Stw,m=(xt,m,yt,m-2a,zt,m),Srw,n=(xr,n,yr,n-2a,zr,n)表示在MP-W檢測路徑中虛擬傳感器的坐標(biāo),詳細見于圖2中。MP-W檢測路徑相對應(yīng)的延遲時間可由如下計算表示

    (4)

    上式中,ft=(xf,t,yf,t,zf,t)與fb=(xf,b,yf,b,zf,b)分別表示缺陷的頂點坐標(biāo)與底點坐標(biāo)。

    當(dāng)缺陷回波與底面強烈的反射波距離足夠遠到可以分離開來[8],利用傳統(tǒng)的技術(shù)(例如,缺陷檢測方法)就能夠很容易地識別出DRP(直達反射路徑)的缺陷回波。

    2 ANSYS有限元仿真

    有限元薄板檢測模型如圖3所示,我們選用被測6061型號鋁板的尺寸為500 mm×100 mm×2 mm。

    圖3 ANSYS有限元仿真模型

    6061鋁板的彈性模量為68.9 GPa,泊松比0.330,密度2.8×103kg/m3;在仿真中我們采用節(jié)點激勵模擬傳感器發(fā)射A0模態(tài)的Lamb波,即從薄板的上邊界激勵。我們用Lamb波的A0模態(tài)作為激勵模態(tài),仿真中的激勵信號[9-10]為

    (5)

    即利用窄帶激勵激發(fā)Lamb波信號,用單一頻率正弦信號加窗函數(shù)(Hanning窗)作為激勵信號,如圖4所示。

    圖4 激勵信號

    A表示信號的幅值A(chǔ)=10 μV;fc表示信號的中心頻率fc=200 kHz;n表示信號的周期數(shù)n=3;H(t)為Heaviside函數(shù),其表示當(dāng)t>0時H(t)=1,t<0時H(t)=0,t=0時H(t)=0.5;我們對ANSYS鋁板模型進行2 mm的網(wǎng)格劃,分割模型單元選擇殼單元,ANSYS的波形采樣頻率為4 MHz。

    2.1相同傳感器位置,相同缺陷位置不同缺陷大小

    傳感器的激勵節(jié)點位于(156,100),接收接收節(jié)點位于(216,100),缺陷的為中心位于(186,50)的圓孔缺陷模型。

    當(dāng)缺陷的直徑為4 mm時,在點(216,100)觀察到的基準模型與缺陷模型的回波如圖5所示。

    圖5 坐標(biāo)(216,100)節(jié)點接收的基準模型與直徑4 mm圓孔缺陷模型的回波

    圖5中,虛線代表4 mm圓孔缺陷模型下的(216,100)節(jié)點觀察回波,實線代表基準模型下(216,100)節(jié)點的觀察回波。由圖中我們可以看到,在有無缺陷模型的情況下首次抵達波在相同接收點的波形完全重合。并且我們可以明顯看出在4 mm圓孔缺陷模型中的DRP(直達反射回波)與MP-W回波。由于在仿真中存在著誤差,且在有缺陷模型中由于缺陷衍射的原因,形成與底面回波的干涉,以致形成波形模態(tài)的轉(zhuǎn)換導(dǎo)致在缺陷模型下底面回波到達觀察點的時間與基準模型下的底面回波形成相位差。圖6表示從仿真模型中提取的缺陷回波并進行Hilbert變換。

    圖6 直徑4 mm圓孔缺陷提取回波與Hilbert變換

    從圖6(b)中我們可以看到提取的缺陷回波經(jīng)過Hilbert變換之后的波形,DRP與MP-W各自缺陷回波信號的延遲時間分別為37.8 μs與100 μs。在ANSYS建模仿真中,我們通過在網(wǎng)格節(jié)點(216,100)測得的首次抵達波的時間為19.33 μs,從而可以得到Lamb波在6061鋁合金薄板中傳播的速度[11-12]為vg=2987.29 m/s。我們將測得的DRP與MP-W的延遲時間分別代入公式(3)與公式(4),可以解得缺陷的頂點坐標(biāo)與底點坐標(biāo)分別為(186,52.17)和(186,46.32)。由此我們可以得到缺陷的大小為5.85 mm。

    同上所述,對直徑8 mm圓孔缺陷模型與12 mm圓孔缺陷模型進行了ANSYS仿真計算,在同樣的激勵點(156,100)進行了仿真,并在相同地(216,100)節(jié)點進行回波觀察,我們得到一下波形,如圖7至圖8所示。

    圖7 坐標(biāo)(216,100)節(jié)點接收的基準模型與直徑8 mm圓孔缺陷模型的回波

    圖8 直徑8 mm圓孔缺陷提取回波與Hilbert變換

    圖9 坐標(biāo)(216,100)節(jié)點接收的基準模型與直徑12 mm圓孔缺陷模型的回波

    圖10 直徑12 mm圓孔缺陷提取回波與Hilbert變換

    由圖8中我們可以分別看出直徑8 mm圓孔缺陷模型的DRP與MP-W的延遲時間分別為37 μs與99.3 μs,圖10中直徑12 mm圓孔缺陷模型的DRP與MP-W的延遲時間分別為35.8 μs與98.5 μs。將兩個模型的延遲時間分為帶入公式(3)與公式(4),我們可以得到8 mm直徑圓孔缺陷模型的頂點坐標(biāo)為(186,53.59),底點坐標(biāo)為(186,45.25);12 mm直徑圓孔缺陷模型的頂點坐標(biāo)為(186,55.74),底點坐標(biāo)為(186,44.03)。

    由此,我們可以得到直徑為8 mm圓孔缺陷模型與12 mm圓孔缺陷模型經(jīng)過ANSYS仿真所測得的缺陷的大小分別為8.34 mm和11.71 mm。表1表示不同直徑缺陷模型所測圓孔大小的相對誤差。

    表1 不同圓孔缺陷模型之間的比較

    通過表1我們可以知道當(dāng)缺陷從4 mm增大到12 mm的過程中,通過多途反射所測得的缺陷的相對誤差逐漸減小,且由圖6(b)、8(b)與10(b)中可知,當(dāng)圓孔缺陷直徑逐漸增大時,DRP與MP-W的缺陷反射回波幅值越來越大,出現(xiàn)時期愈加明顯。因此我們可以得出這樣的結(jié)論,在模型不變的情況下,圓孔缺陷越大多途檢測效果越好,缺陷尺寸愈加精確相近。從上述仿真測得的結(jié)果可知,在相同節(jié)點中觀察到不同模型DRP與MP-W傳播路徑的缺陷回波,DPR是用來檢測缺陷的頂點位置,MP-W是用來檢測缺陷的地點坐標(biāo),通過對觀察到的缺陷回波進行Hilbert變換可以直接讀取不同路徑缺陷回波的延遲時間,通過DRP與MP-W路徑的延遲時間可以較為精確的計算出缺陷的大小,缺陷越大,結(jié)果越精確。所以超聲多途檢測方法只需通過單次測量就可以對缺陷進行定位并測出其大小。仿真結(jié)果最后所造成的缺陷的誤差主要是由于ANSYS建模在圓孔缺陷以及傳感器附近網(wǎng)格劃分不均,以及當(dāng)缺陷越小時,經(jīng)過缺陷的衍射波會對底面波形成更加嚴重的干涉所造成的。

    2.2相同缺陷位置及大小,不同傳感器位置

    當(dāng)缺陷的尺寸不變,發(fā)射傳感器與接收傳感器沿著-y軸方向每次按步長10 mm進行移動,對每次移動后的模型進行檢測,由于薄板的寬度不變?yōu)?00 mm,所以總共移動10次,假設(shè)移動一次,則發(fā)射與接收傳感器的編號變?yōu)镾t1與Rt1,則若移動m次,則變?yōu)镾tm與Rtm,如圖11所示。

    圖11 同種缺陷不同收發(fā)傳感器位置的模型

    對應(yīng)的圖11,假設(shè)圓孔缺陷的直徑為8 mm,中心坐標(biāo)為(186,50),在其他條件不變的情況下即與2.1的仿真條件類似。每種情形下收發(fā)傳感器St1、Rt1至St10、Rt10對應(yīng)的坐標(biāo)分別為:

    收發(fā)傳感器St1、Rt1至St10、Rt10每種情況下的仿真結(jié)果如圖12~圖21所示。

    圖12 收發(fā)傳感器位置分別在(156,90)與(216,90)時的仿真結(jié)果

    圖13 收發(fā)傳感器位置分別在(156,80)與(216,80)時的仿真結(jié)果

    圖12(a)、(b)與(c)分別表示基準模型與缺陷模型的回波信號與在網(wǎng)格觀察節(jié)點(216,90)提取出的缺陷回波以及缺陷回波的Hilbert變換。從圖12(a)中可以看出,直接反射回波DRP的延時時間還是非常容易判斷的,但是檢測缺陷底部的MP-W的反射回波卻被掩蓋了,不能較為明顯的觀察出其時間,主要是由于當(dāng)收發(fā)射傳感器的向-y方向移動1次后,上下邊界的反射回波對MP-W的判斷出現(xiàn)了干擾出現(xiàn)了重疊現(xiàn)象,以及Lamb波在傳播過程中存在不同路徑波相互之間的干涉所造成的。

    從圖13(a)中難以判斷出直接反射回波DRP出現(xiàn)的時間,主要是由于當(dāng)收發(fā)射傳感器的向-y方向移動2次后,邊界面B的回波會與DRP路徑的回波造成重疊,且到達時間差不多一致,所以較難判斷;MP-W路徑的回波已基本被其他回波所掩蓋,更加難以識別。

    圖14 收發(fā)傳感器位置分別在(156,70)與(216,70)時的仿真結(jié)果

    圖15 收發(fā)傳感器位置分別在(156,60)與(216,60)時的仿真結(jié)果

    圖16 收發(fā)傳感器位置分別在(156,50)與(216,50)時的仿真結(jié)果

    圖17 收發(fā)傳感器位置分別在(156,40)與(216,40)時的仿真結(jié)果

    圖18 收發(fā)傳感器位置分別在(156,30)與(216,30)時的仿真結(jié)果

    圖19 收發(fā)傳感器位置分別在(156,20)與(216,20)時的仿真結(jié)果

    圖20 收發(fā)傳感器位置分別在(156,10)與(216,10)時的仿真結(jié)果

    圖21 收發(fā)傳感器位置分別在(156,0)與(216,0)時的仿真結(jié)果

    從圖12~圖21可以看出,圖12與圖20、圖13與圖19、圖14與圖18、圖15與圖17這四組仿真結(jié)果圖兩兩相對都是一樣的,這是因為當(dāng)收發(fā)射傳感器的向-y方向逐漸移動,他們的位置關(guān)系會關(guān)于y=50對稱;圖16中收發(fā)傳感器的位置正好位于y=50這條直線上,且由于缺陷的中心位置也在這條線上,所以接收不到直接反射路徑DRP的回波,由于位于中間位置,上下底面的回波會與MP-W的路徑嚴重重疊,且兩個底面反射回來的波干涉嚴重,會使分析更加困難。

    在圖13至圖19中,難以判斷DRP與MP-W的出現(xiàn)時間,主要是由于收發(fā)射傳感器的向-y方向逐漸移動時,上下底面的反射回波會造成嚴重的干擾,這給回波分析帶來了困難。在圖12與圖20中,兩者的方針效果一樣,且能夠較為清楚的判斷直接反射路徑DRP的出現(xiàn)時間,但不能確定MP-W出現(xiàn)的時間,主要是由于當(dāng)收發(fā)射傳感器的向-y方向移動1次后,上下邊界的反射回波對MP-W的判斷出現(xiàn)了干擾出現(xiàn)了重疊現(xiàn)象,以及Lamb波傳播過程中的不同路徑波的干涉所造成的。

    從圖21中可以容易判斷出DRP與MP-W的延遲時間,由于對稱關(guān)系,圖21與2.1節(jié)模型中圖7、圖8的結(jié)果一致,由于圖21的收發(fā)傳感器位于邊界位置,且是對邊界節(jié)點進行的激勵加載,所以避免了底面反射回波所帶來的干擾,能夠得到較為滿意的結(jié)果。所以當(dāng)缺陷位于中心位置不變時,收發(fā)傳感器處在邊界上會是仿真觀察的最佳效果。

    2.3相同傳感器位置,相同缺陷大小不同缺陷位置

    與2.2節(jié)相類似的,為了從不同角度驗證多途超聲檢測法在薄板上的檢測能力,我們令收發(fā)傳感器的位置保持在St與Rt不變,即坐標(biāo)(156,100)與(216,100)。讓缺陷的大小固定為8 mm直徑的圓孔,使其從缺陷中心為(186,50)的位置按每次步長為10 mm向-y方向移動,總共移動5次,并對每次移動后的模型進行ANSYS建模仿真。其原理圖如下圖22所示。

    圖22 相同傳感器位置不同缺陷位置的模型

    圖22中,對應(yīng)的不同缺陷位置的中心坐標(biāo)分別為:

    F1:(186,40)、F2:(186,30)、F3:(186,20)、

    F4:(186,10)、F5:(186,0)。

    其他的條件相對與2.2節(jié)的仿真條件不變。當(dāng)缺陷位于不同的位置,即F1、F2、F3、F4與F5時,每種情況下ANSYS的仿真結(jié)果如圖23~圖27所示。由于當(dāng)缺陷的中心位于F5:(186,0)時,已不適用與超聲多途檢測方法的應(yīng)用,于是不對此種情形進行建模仿真。

    圖23 缺陷的中心坐標(biāo)為(186,50)時的仿真結(jié)果

    圖24 缺陷的中心坐標(biāo)為(186,40)時的仿真結(jié)果

    從圖23中我們可以清楚地看到在缺陷中心坐標(biāo)為(186,50)時,ANSYS建模仿真的結(jié)果,從圖23我們可以得到直接反射路徑DRP的延遲時間為37 μs,MP-W的時間為99.3 μs,由于之前在2.1節(jié)中我們計算得到Lamb的群速度為vg=2 987.29 m/s,所以將時間帶入公式(3)與公式(4),可以計算得到缺陷的頂點坐標(biāo)與底點坐標(biāo)分別為(186,53.59)與(186,45.25),此時可以計算得到缺陷的直徑為8.34 mm與ANSYS實際模型模擬的缺陷8 mm十分的接近。

    仿真結(jié)果最后所造成的缺陷的誤差主要是由于ANSYS建模在圓孔缺陷以及傳感器附近網(wǎng)格劃分不均,以及當(dāng)缺陷越小時,經(jīng)過缺陷的衍射波會對底面波形成更加嚴重的干涉所造成的。

    從圖24中還是可以較為清楚地看到在缺陷中心坐標(biāo)為(186,40)時,對其進行ANSYS建模仿真的結(jié)果,從圖24(c)中缺陷回波的Hilbert變換可以得到,直接反射路徑DRP的延遲時間為41.8 μs,MP-W路徑的延遲時間為93.3 μs,同樣的將其帶入公式(3)和公式(4)中,可以計算得到缺陷的頂點坐標(biāo)與底點坐標(biāo)分別為(186,45.25)與(186,36.09),此時可以計算得到缺陷的直徑為9.16 mm。

    圖25 缺陷的中心坐標(biāo)為(186,30)時的仿真結(jié)果

    圖26 缺陷的中心坐標(biāo)為(186,20)時的仿真結(jié)果

    圖25表示在缺陷中心坐標(biāo)為(186,30)時,ANSYS建模仿真的結(jié)果。從圖25(c)中缺陷回波的Hilbert變換可以得到,直接反射路徑DRP的延遲時間為46.8 μs,但是對于MP-W的延遲時間判斷出現(xiàn)了干擾,我們經(jīng)過思考最終認為MP-W路徑的延遲時間為90 μs,可以見到在圖25(c)中MP-W周圍出現(xiàn)了3個Hilbert變換峰值,主要是由于隨著缺陷中心的下移,底部反射回波的傳播路徑與缺陷底點較為相近,從而部分的回波造成了干涉,導(dǎo)致MP-W的路徑回波時間更加持久。

    同樣的將其帶入公式(4)與公式(3),可以計算得到缺陷的頂點坐標(biāo)與底點坐標(biāo)分別為(186,36.86)與(186,31.04),此時可以計算得到缺陷的直徑為5.82 mm。

    圖26表示在缺陷中心坐標(biāo)為(186,20)時,ANSYS建模仿真的結(jié)果。從圖26(a)中我們可以看到,直接反射路徑DRP回波的延遲時間已經(jīng)伴隨著缺陷中心的下移而逐漸往后移,并與底面回波的延遲時間較為相近,但此時還是可以區(qū)分出DRP的延時時間,從圖26(c)中缺陷回波的Hilbert變換可以得到,直接反射路徑DRP的延遲時間為50.3 μs。對于MP-W的延遲時間判斷較為困難,因為隨著缺陷中心的下移,MP-W傳播的路程變短,導(dǎo)致底面回波對于MP-W路徑的影響加大,我們經(jīng)過思考最終認為MP-W路徑的延遲時間為90.3 μs,可以見到在圖26(c)中MP-W周圍出現(xiàn)了若干個Hilbert變換峰值,主要是由于隨著缺陷中心的下移,底部反射回波的傳播路徑與缺陷底點較為相近,從而部分的回波造成了干涉,并且兩者到達接收節(jié)點的時間相近,導(dǎo)致MP-W的路徑的延遲時間確定更為困難。

    圖27 缺陷的中心坐標(biāo)為(186,10)時的仿真結(jié)果

    同樣的將其帶入公式(3)與公式(4),可以計算得到缺陷的頂點坐標(biāo)與底點坐標(biāo)分別為(186,31.12)與(186,31.5),此時可以計算得到缺陷的直徑為-0.38 mm。在圖27中,此時缺陷中心的坐標(biāo)為(186,10),與底面十分的接近。從圖27(a)中看出直接反射路徑DRP與底面回波已經(jīng)重疊在一起,不能判斷出其延遲時間。再者缺陷的底點與地面的距離只有6 mm,MP-W的路徑檢測會受到底面多次反射、散射的干擾,或者與底面回波重疊,從而很難對MP-W的延遲時間進行估計。

    從圖23至圖27的仿真結(jié)果可以看出,在缺陷中心位于(186,50)與(186,40)時,對于Lamb波多途檢測路徑DPR與MP-W延遲時間的判定不存在干擾;當(dāng)缺陷中心位于(186,30)與(186,20)時,此時對于路徑DPR延遲時間的估計較為清晰,但是對于MP-W延遲時間的確定卻出現(xiàn)了較為嚴重的干擾;當(dāng)圓孔缺陷中心位于(186,10)時,已經(jīng)不能判斷出DRP與MP-W的延遲時間。

    從表2中我們可以看到,在圖23至圖27的仿真結(jié)果中,ANSYS建模時缺陷的大小為直接8 mm的圓孔,在不同缺陷中心位置,建模仿真并計算得到了不同模型的缺陷大小以及相對誤差。從表2可以看到在缺陷中心位置為(186,50)時,仿真所得的結(jié)果最為相近,隨著缺陷中心的下移,ANSYS仿真得到的缺陷的大小相對誤差越大,乃至在(186,20)與(186,10)時,判斷缺陷的大小出現(xiàn)了嚴重的錯誤。

    表2 相同傳感器位置不同缺陷位置模型之間的比較

    綜上所述,當(dāng)缺陷中心位于(186,50)時,結(jié)果最為準確,仿真效果最好,且此時的多途檢測能力最強,隨著缺陷中心的下移,對于多途路徑DRP與MP-W的延遲時間的判斷越來越難,主要是由于底面回波到達接收節(jié)點的時間與缺陷回波到達接收節(jié)點的時間越來越相近,并且當(dāng)缺陷位于底面附近時,底面回波會與MP-W路徑的回波相成較為嚴重的干涉,導(dǎo)致傳播時間無法預(yù)估。所以當(dāng)缺陷位于板的中心附近,且從薄板的邊界激勵,Lamb波多途檢測效果最佳。且超聲多途檢測方法只需通過單次測量就可以對缺陷進行定位并測出其大小是可行的。仿真結(jié)果最后所造成的缺陷的誤差主要是由于ANSYS建模在圓孔缺陷以及傳感器附近網(wǎng)格劃分不均,以及當(dāng)缺陷越小時,經(jīng)過缺陷的衍射波會對底面波形成更加嚴重的干涉所造成的。

    3 結(jié)論

    利用超聲多途提供的額外數(shù)據(jù)來檢測金屬薄板的缺陷,提出相應(yīng)的算法,多途路徑在直接反射的路徑上是可以預(yù)測和識別的,直接反射路徑是相對比較容易識別的。同時利用ANSYS有限元建模不同的薄板缺陷模型,在不同的傳感器位置、缺陷位置、缺陷大小情況下分別進行仿真。實驗表明超聲多途所提供的額外數(shù)據(jù)能夠用于薄板缺陷的檢測。

    本文主要根據(jù)超聲多途Lamb波檢測法原理,建立了一種單發(fā)射單接收的傳感器檢測模型,研究了其中的關(guān)鍵技術(shù),并應(yīng)用于金屬鋁合金薄板探傷,在應(yīng)用領(lǐng)域和應(yīng)用技術(shù)方面取得了新的進展;通過ANSYS仿真和實驗檢測驗證了超聲多途的可行性;提出的超聲多途檢測方法,通過單次測量即可檢測出薄板缺陷的位置和大小,具有良好的檢測性能。

    [1]Meksen T M, Boudraa B, Drai R, et al. Automatic crack detection and characterization during ultrasonic inspection [J].Journal of Nondestructive Evaluation, 2010,29(3):169-174.

    [2]Lester W.Schmerr Jr, Sung-Jin Song. Ultrasonic Nondestructive Evaluation Systems Models and Measurements [M]. Springer Science Business Media, New York, USA,2007.

    [3]Fukuda M, Imano K, Yamagishi H, et al. Detection of the Second Harmonics of Lamb Waves in Fatigued Magnesium Plates[J]. Acoustical Science and Technology,2011,32(6):271-275.

    [4]Dewhurst R.J, Edwords C, Mckie A.D.W, et al. Estimation of the thickness of thin metal sheet using laser generated ultrasound [J]. Applied Physics Letters,1987,51(14):1066-1068.

    [5]Hutchins D A, Lundgren K, Plamer S B. A Laser Study of Transient Lamb Waves in Thin Materials [J]. Journal of the Acoustical Society of America,1989,85(4):1441-1448.

    [6]Rose J L. Ultrasonic Waves in Solid Media [M]. UK: Cambridge University Press,1999.

    [7]Ostachowicz W, Kudela P, Malinowski P, et al. Damage Localization in Plate-like Structures Based on PZT Sensors [J]. Mechanical Systems and Signal Processing,2009,23:1805-1829.

    [8]Park H W, Sohn H, Law K H, et al. Time Reversal Active Sensing for Health Monitoring of A Composite Plate [J]. Journal of Sound and Vibration, 2007,302(1-2):50-66.

    [9]馮勇明,周麗.基于Lamb波時間反轉(zhuǎn)的復(fù)合材料結(jié)構(gòu)傷檢測[J].中國機械工程,2011,22(20):2476-2480.

    [10]Shen X Z, Pan L. material Flaw Sizing By Ultrasonic Multipath Detection[J].Periodical of Advanced Materials research vol.712-715,pp:1067-1070.

    [11]袁慎芳.結(jié)構(gòu)健康監(jiān)控[M].北京:國防工業(yè)出版社,2007.

    [12]何存富,李隆濤,吳斌.周向超聲導(dǎo)波在薄壁管道中的傳播研究[J].實驗力學(xué),2002,17(4):419-424.

    [13]趙彩萍,王維斌,何存富,等.基于神經(jīng)網(wǎng)絡(luò)模型的超聲導(dǎo)波管道缺陷識別[J].傳感器與微系統(tǒng),2009,28(11):19-21.

    [14]李富才,孟光.窄頻帶Lamb波頻散特性研究[J].實驗力學(xué),2002,17(4):419-424.

    [15]王軍,王寅觀.板中正交靜應(yīng)力與任意方向蘭姆波波速的關(guān)系[J].力學(xué)學(xué)報,2008,40(3):345-354.

    [16]代樹武,張海瀾,粗糙界面板Lamb波的傳播特征[J].聲學(xué)學(xué)報,2003,28(3):193-200.

    [17]閻石,張海鳳,蒙彥宇.Lamb波頻散曲線的數(shù)值計算及試驗驗證[J].華中科技大學(xué)學(xué)報(城市科學(xué)版),2010,27(1):1-4.

    Ultrasonic Mutipath Lamb Wave Detection of Sheet Metal Defects

    Guo Dubin, Shen Xizhong, Wang Lei

    (School of Electrical and Electronic Engineering, Shanghai Institute of Technology, Shanghai201418,China)

    Sheet metal in the moving car is widely used in modern design, such as detection of defects has very important significance. For this purpose, the use of ultrasound of multiplicity provide additional data to detect the defects of sheet metal, and the algorithm of multiplicity is predictable path in the direct reflection of the path and identify, directly reflect the path is relatively easy to identify. Using ANSYS finite element modeling different sheet defect model at the same time, the different sensor location, defect location, defect size cases are simulated respectively. Obtained by ANSYS simulation and experimental tests are good results, the error in the range of controlled, proved that the ultrasonic application of multiplicity and sheet metal Lamb wave detection is feasible; Simulation and experiments show that ultrasonic way through actual and virtual extension of multiplicity of sensor data analysis proved that the ultrasonic multiplicity can for detecting sheet to locate defects and their size. Ultrasound of multiplicity of detection method is put forward, through a single measurement can detect the defect size and location of the plate, has good detection performance.

    ultrasound; sheet; multipath; ANSYS finite element simulation; flaw; size

    2015-07-13;

    2015-09-16。

    上海市科委項目(15ZR1440700)。

    郭杜斌(1990-),男,浙江人,碩士研究生,主要從事超聲檢測方向的研究。

    沈希忠(1968-),男,上海人,主要從事信號處理、超聲檢測方向的研究。

    1671-4598(2016)01-0020-08

    10.16526/j.cnki.11-4762/tp.2016.01.006

    TB114.3

    A

    猜你喜歡
    延遲時間圓孔薄板
    超聲心動圖在可降解房間隔封堵器治療卵圓孔未閉術(shù)中的應(yīng)用價值
    一角點支撐另一對邊固支正交各向異性矩形薄板彎曲的辛疊加解
    二氧化碳對乙烷燃燒著火延遲時間的影響
    煤氣與熱力(2021年3期)2021-06-09 06:16:22
    LTE 系統(tǒng)下行鏈路FDRX 節(jié)能機制研究
    10MN鋁合金薄板拉伸機組的研制
    重型機械(2020年2期)2020-07-24 08:16:12
    基于分層COX模型的跟馳反應(yīng)延遲時間生存分析
    延遲時間對氣輔注射成型氣體穿透行為影響的數(shù)值模擬和實驗研究
    中國塑料(2016年8期)2016-06-27 06:35:02
    鋁薄板高速DP-GMAW焊接性能的研究
    焊接(2016年5期)2016-02-27 13:04:42
    基于PLC的薄板激光焊機控制系統(tǒng)
    焊接(2015年10期)2015-07-18 11:04:46
    基于光電檢測和PLC控制的自動化圓孔倒角設(shè)備的設(shè)計與實現(xiàn)
    91成年电影在线观看| 一进一出抽搐动态| 精品第一国产精品| 老熟妇乱子伦视频在线观看| 亚洲全国av大片| 啦啦啦免费观看视频1| 国产一区二区在线av高清观看| 午夜福利在线观看吧| 很黄的视频免费| 中文字幕久久专区| 成人av一区二区三区在线看| 日本精品一区二区三区蜜桃| 大型av网站在线播放| 国产亚洲欧美98| 国产精品野战在线观看| 老熟妇仑乱视频hdxx| 99热只有精品国产| 亚洲精品久久国产高清桃花| 伦理电影免费视频| 久久久久久免费高清国产稀缺| 他把我摸到了高潮在线观看| 亚洲人成电影观看| 色播在线永久视频| 大型黄色视频在线免费观看| 国产亚洲欧美在线一区二区| 欧美成狂野欧美在线观看| 色尼玛亚洲综合影院| 99国产精品一区二区三区| 嫩草影院精品99| 国产精品爽爽va在线观看网站 | 国产av精品麻豆| 国产欧美日韩一区二区精品| 中文字幕av电影在线播放| 不卡一级毛片| 亚洲中文av在线| 制服诱惑二区| 淫妇啪啪啪对白视频| 女同久久另类99精品国产91| 一进一出好大好爽视频| 91国产中文字幕| 久久影院123| av免费在线观看网站| 亚洲国产看品久久| 国产成人精品久久二区二区91| 欧美色欧美亚洲另类二区 | 亚洲七黄色美女视频| 操美女的视频在线观看| 久久精品亚洲熟妇少妇任你| 午夜福利欧美成人| 一级片免费观看大全| 亚洲avbb在线观看| 亚洲五月色婷婷综合| 男女做爰动态图高潮gif福利片 | 亚洲午夜精品一区,二区,三区| 波多野结衣av一区二区av| 国产av精品麻豆| 久久久久久亚洲精品国产蜜桃av| 欧美日韩瑟瑟在线播放| 在线观看免费视频网站a站| 亚洲av熟女| 人妻久久中文字幕网| 欧美日韩瑟瑟在线播放| 午夜福利一区二区在线看| 夜夜爽天天搞| 亚洲国产精品久久男人天堂| 国产精品亚洲一级av第二区| 美女高潮喷水抽搐中文字幕| 一进一出抽搐动态| 色综合婷婷激情| 好男人电影高清在线观看| 欧美激情极品国产一区二区三区| 人人妻人人澡人人看| 日日摸夜夜添夜夜添小说| 999久久久精品免费观看国产| 欧美黑人精品巨大| 桃红色精品国产亚洲av| 一级毛片高清免费大全| 天天躁狠狠躁夜夜躁狠狠躁| 51午夜福利影视在线观看| www国产在线视频色| 欧美 亚洲 国产 日韩一| 免费少妇av软件| 叶爱在线成人免费视频播放| 大码成人一级视频| 欧美精品啪啪一区二区三区| 久久精品aⅴ一区二区三区四区| 日本 欧美在线| 亚洲男人的天堂狠狠| 国产乱人伦免费视频| 少妇被粗大的猛进出69影院| 久久久国产成人免费| 久久 成人 亚洲| 日韩精品青青久久久久久| 亚洲性夜色夜夜综合| 在线av久久热| 久久香蕉精品热| 欧美黄色片欧美黄色片| 操出白浆在线播放| 国产精品乱码一区二三区的特点 | 制服人妻中文乱码| www.自偷自拍.com| 99国产精品一区二区蜜桃av| 91老司机精品| 男女午夜视频在线观看| 男女之事视频高清在线观看| 久久影院123| 黄色 视频免费看| 亚洲成av片中文字幕在线观看| 午夜精品久久久久久毛片777| 中文字幕av电影在线播放| 婷婷精品国产亚洲av在线| 91av网站免费观看| 国产人伦9x9x在线观看| 国产精品综合久久久久久久免费 | 人人澡人人妻人| 色哟哟哟哟哟哟| 视频区欧美日本亚洲| 夜夜夜夜夜久久久久| 在线十欧美十亚洲十日本专区| 一a级毛片在线观看| 在线观看午夜福利视频| 丝袜美足系列| 日韩精品免费视频一区二区三区| 一边摸一边抽搐一进一小说| 午夜免费鲁丝| 亚洲成人精品中文字幕电影| 久久精品91无色码中文字幕| 日日摸夜夜添夜夜添小说| 男女午夜视频在线观看| 免费看十八禁软件| 国产1区2区3区精品| 欧美日本亚洲视频在线播放| 国产成人影院久久av| 亚洲情色 制服丝袜| 亚洲精品在线美女| 亚洲国产毛片av蜜桃av| 久久精品影院6| 十八禁人妻一区二区| 一级毛片女人18水好多| 国产在线观看jvid| 久久性视频一级片| 香蕉久久夜色| 欧美激情 高清一区二区三区| 99在线视频只有这里精品首页| 国产成人精品久久二区二区免费| 亚洲精品久久成人aⅴ小说| 两性夫妻黄色片| 中文字幕高清在线视频| 两个人免费观看高清视频| 欧美成狂野欧美在线观看| 热re99久久国产66热| 啦啦啦免费观看视频1| 这个男人来自地球电影免费观看| 看免费av毛片| 午夜a级毛片| 制服诱惑二区| 老司机深夜福利视频在线观看| 男女做爰动态图高潮gif福利片 | 欧美亚洲日本最大视频资源| 性色av乱码一区二区三区2| 色尼玛亚洲综合影院| 黄色 视频免费看| 叶爱在线成人免费视频播放| 成人18禁在线播放| 婷婷精品国产亚洲av在线| 国产精品久久久久久人妻精品电影| 国产精品一区二区在线不卡| 午夜影院日韩av| 国产精品久久久久久人妻精品电影| 亚洲在线自拍视频| 国产又色又爽无遮挡免费看| 女警被强在线播放| 国产乱人伦免费视频| 亚洲精品美女久久av网站| 午夜老司机福利片| 日韩精品中文字幕看吧| 国产精品二区激情视频| 校园春色视频在线观看| 伊人久久大香线蕉亚洲五| 欧美日韩亚洲综合一区二区三区_| 少妇的丰满在线观看| 精品不卡国产一区二区三区| 视频区欧美日本亚洲| 中文字幕色久视频| 一区二区三区高清视频在线| 老鸭窝网址在线观看| 欧美久久黑人一区二区| 日本五十路高清| 老鸭窝网址在线观看| 可以在线观看毛片的网站| 免费看a级黄色片| 色精品久久人妻99蜜桃| 国产亚洲精品一区二区www| 精品人妻在线不人妻| 久久香蕉激情| 美女国产高潮福利片在线看| 国产一区二区激情短视频| 成人亚洲精品av一区二区| 不卡一级毛片| 夜夜爽天天搞| 精品久久久精品久久久| 91麻豆精品激情在线观看国产| 热99re8久久精品国产| 9热在线视频观看99| 国产欧美日韩精品亚洲av| 91麻豆av在线| 成人国产一区最新在线观看| 亚洲人成电影观看| 女性被躁到高潮视频| 中文字幕人妻丝袜一区二区| 亚洲av电影不卡..在线观看| 大型黄色视频在线免费观看| 午夜免费激情av| 国产亚洲精品综合一区在线观看 | 国产野战对白在线观看| 一本久久中文字幕| 老司机深夜福利视频在线观看| 免费在线观看日本一区| 美女大奶头视频| 91国产中文字幕| 超碰成人久久| 亚洲中文av在线| 欧美中文综合在线视频| 午夜福利18| 男人的好看免费观看在线视频 | 亚洲在线自拍视频| 国产成人系列免费观看| 一进一出抽搐动态| 国产午夜福利久久久久久| 国内毛片毛片毛片毛片毛片| 人妻丰满熟妇av一区二区三区| 亚洲国产精品sss在线观看| 国产色视频综合| 免费在线观看视频国产中文字幕亚洲| 国产一区二区激情短视频| 国产成人啪精品午夜网站| 国产成人欧美| 久久久久国内视频| 天天一区二区日本电影三级 | 一本大道久久a久久精品| 深夜精品福利| 人妻久久中文字幕网| 日本免费一区二区三区高清不卡 | 女人被躁到高潮嗷嗷叫费观| 麻豆av在线久日| 黑人操中国人逼视频| 欧美午夜高清在线| 校园春色视频在线观看| 亚洲成人精品中文字幕电影| 国产成人av教育| 在线十欧美十亚洲十日本专区| 久久精品影院6| 午夜久久久在线观看| 自线自在国产av| 嫁个100分男人电影在线观看| 好看av亚洲va欧美ⅴa在| 99热只有精品国产| 丰满人妻熟妇乱又伦精品不卡| 日本五十路高清| 国产色视频综合| 日韩精品中文字幕看吧| 久久久精品欧美日韩精品| 人人澡人人妻人| 精品国产一区二区三区四区第35| 夜夜夜夜夜久久久久| 老汉色∧v一级毛片| 久久久久久久精品吃奶| 国产精品亚洲美女久久久| АⅤ资源中文在线天堂| 午夜成年电影在线免费观看| 国产精品久久久av美女十八| 精品一区二区三区四区五区乱码| 身体一侧抽搐| 午夜福利视频1000在线观看 | 嫁个100分男人电影在线观看| 亚洲自偷自拍图片 自拍| 老汉色av国产亚洲站长工具| 亚洲国产精品久久男人天堂| 日韩 欧美 亚洲 中文字幕| 午夜老司机福利片| 每晚都被弄得嗷嗷叫到高潮| 禁无遮挡网站| 国产亚洲精品久久久久久毛片| 国产高清有码在线观看视频 | 免费在线观看亚洲国产| 12—13女人毛片做爰片一| av超薄肉色丝袜交足视频| 精品卡一卡二卡四卡免费| 亚洲久久久国产精品| 午夜成年电影在线免费观看| 免费看a级黄色片| 69精品国产乱码久久久| 欧美日韩黄片免| 最近最新免费中文字幕在线| 亚洲av五月六月丁香网| 51午夜福利影视在线观看| 精品久久久久久久人妻蜜臀av | 精品欧美一区二区三区在线| 午夜福利18| 午夜精品久久久久久毛片777| 777久久人妻少妇嫩草av网站| 精品国产乱码久久久久久男人| 国产1区2区3区精品| 国产91精品成人一区二区三区| 可以免费在线观看a视频的电影网站| 99久久久亚洲精品蜜臀av| 亚洲国产欧美一区二区综合| 在线十欧美十亚洲十日本专区| 久久精品人人爽人人爽视色| 人妻丰满熟妇av一区二区三区| 国产精品久久久人人做人人爽| 欧美色欧美亚洲另类二区 | 在线av久久热| 国产欧美日韩一区二区三区在线| 亚洲色图 男人天堂 中文字幕| 又紧又爽又黄一区二区| 啪啪无遮挡十八禁网站| 女性被躁到高潮视频| 90打野战视频偷拍视频| 国产激情欧美一区二区| 啦啦啦 在线观看视频| 级片在线观看| 色综合亚洲欧美另类图片| 少妇的丰满在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲人成伊人成综合网2020| 老汉色∧v一级毛片| 国产成人影院久久av| 国产精品99久久99久久久不卡| 免费观看人在逋| 热re99久久国产66热| 啦啦啦 在线观看视频| 淫妇啪啪啪对白视频| 91麻豆精品激情在线观看国产| 成人亚洲精品av一区二区| 亚洲国产精品sss在线观看| 日本撒尿小便嘘嘘汇集6| 每晚都被弄得嗷嗷叫到高潮| av福利片在线| 欧美一区二区精品小视频在线| 日韩大码丰满熟妇| 国产精品1区2区在线观看.| 18禁美女被吸乳视频| 国产亚洲精品一区二区www| 久久精品人人爽人人爽视色| 一区二区三区国产精品乱码| 欧美激情极品国产一区二区三区| 亚洲精品av麻豆狂野| 国产一区二区三区综合在线观看| 久久青草综合色| 88av欧美| 欧美最黄视频在线播放免费| 国产一区二区三区综合在线观看| 丁香六月欧美| 女性生殖器流出的白浆| 丝袜美足系列| 久久中文字幕人妻熟女| 国产成+人综合+亚洲专区| 黑丝袜美女国产一区| 99国产精品99久久久久| 免费高清在线观看日韩| 欧美乱码精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 久久精品亚洲精品国产色婷小说| or卡值多少钱| 一本综合久久免费| 亚洲欧美一区二区三区黑人| 免费在线观看视频国产中文字幕亚洲| 久久精品国产亚洲av高清一级| 咕卡用的链子| 免费看美女性在线毛片视频| 欧美乱码精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 岛国在线观看网站| 亚洲精品中文字幕一二三四区| 好看av亚洲va欧美ⅴa在| 日韩视频一区二区在线观看| 国产伦一二天堂av在线观看| 亚洲欧美一区二区三区黑人| 女性被躁到高潮视频| 亚洲av熟女| 欧美精品啪啪一区二区三区| 最近最新中文字幕大全免费视频| 757午夜福利合集在线观看| 少妇裸体淫交视频免费看高清 | 夜夜夜夜夜久久久久| 精品人妻1区二区| 热99re8久久精品国产| 欧美一区二区精品小视频在线| 久99久视频精品免费| 男女下面插进去视频免费观看| 天天一区二区日本电影三级 | 不卡一级毛片| 免费在线观看黄色视频的| 嫩草影视91久久| x7x7x7水蜜桃| 欧美+亚洲+日韩+国产| 日韩有码中文字幕| 乱人伦中国视频| 亚洲人成伊人成综合网2020| 午夜影院日韩av| 亚洲第一欧美日韩一区二区三区| 男女下面进入的视频免费午夜 | 久久青草综合色| 天堂动漫精品| 波多野结衣巨乳人妻| 国产亚洲欧美在线一区二区| av电影中文网址| 制服丝袜大香蕉在线| 日韩国内少妇激情av| 免费高清在线观看日韩| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品一区二区精品视频观看| 18禁裸乳无遮挡免费网站照片 | 日本免费一区二区三区高清不卡 | 日韩精品免费视频一区二区三区| 国产区一区二久久| 每晚都被弄得嗷嗷叫到高潮| 在线观看日韩欧美| 人人妻,人人澡人人爽秒播| 亚洲人成电影观看| 老熟妇乱子伦视频在线观看| 高清毛片免费观看视频网站| 国产av在哪里看| 两性夫妻黄色片| 色哟哟哟哟哟哟| 好男人电影高清在线观看| 99久久国产精品久久久| 村上凉子中文字幕在线| 美女扒开内裤让男人捅视频| 岛国视频午夜一区免费看| 看免费av毛片| 国产精品av久久久久免费| netflix在线观看网站| 一夜夜www| 一边摸一边抽搐一进一出视频| 亚洲精华国产精华精| 黄网站色视频无遮挡免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲 国产 在线| 亚洲欧美激情综合另类| 国产不卡一卡二| 亚洲精品国产色婷婷电影| 欧美黑人欧美精品刺激| 国产成人精品久久二区二区免费| 最近最新中文字幕大全免费视频| 麻豆成人av在线观看| 给我免费播放毛片高清在线观看| 久久欧美精品欧美久久欧美| 精品午夜福利视频在线观看一区| 国产一卡二卡三卡精品| 99精品在免费线老司机午夜| 久久中文字幕一级| 久久国产精品影院| 久久久国产成人免费| 人人妻,人人澡人人爽秒播| 精品久久蜜臀av无| 国产私拍福利视频在线观看| 国产精品电影一区二区三区| 制服丝袜大香蕉在线| 老熟妇乱子伦视频在线观看| 精品欧美国产一区二区三| 999久久久国产精品视频| 欧美中文日本在线观看视频| 国产精品久久久久久精品电影 | 欧美绝顶高潮抽搐喷水| 久久精品91无色码中文字幕| 每晚都被弄得嗷嗷叫到高潮| 性欧美人与动物交配| 日本 av在线| 国产日韩一区二区三区精品不卡| 琪琪午夜伦伦电影理论片6080| 成人永久免费在线观看视频| 一级,二级,三级黄色视频| 身体一侧抽搐| 成人国语在线视频| 欧美午夜高清在线| 热re99久久国产66热| 亚洲精品中文字幕一二三四区| 麻豆成人av在线观看| 久99久视频精品免费| 天天躁狠狠躁夜夜躁狠狠躁| 制服人妻中文乱码| 久久久国产精品麻豆| 欧美 亚洲 国产 日韩一| www.www免费av| 女人被躁到高潮嗷嗷叫费观| 大码成人一级视频| 日韩欧美三级三区| 精品卡一卡二卡四卡免费| 国产精品久久久久久人妻精品电影| 妹子高潮喷水视频| 国产区一区二久久| 欧美老熟妇乱子伦牲交| 成人精品一区二区免费| 老鸭窝网址在线观看| 老汉色∧v一级毛片| 国产97色在线日韩免费| 久久精品国产亚洲av香蕉五月| 久久青草综合色| 免费观看人在逋| 国产视频一区二区在线看| 日本一区二区免费在线视频| 欧美在线黄色| 黄片播放在线免费| 黑丝袜美女国产一区| 国产精品免费视频内射| 女警被强在线播放| 一区二区日韩欧美中文字幕| 怎么达到女性高潮| 亚洲精品粉嫩美女一区| 啦啦啦韩国在线观看视频| 亚洲狠狠婷婷综合久久图片| 欧美不卡视频在线免费观看 | 少妇被粗大的猛进出69影院| 麻豆一二三区av精品| 亚洲av第一区精品v没综合| 日本黄色视频三级网站网址| 操美女的视频在线观看| 亚洲成人国产一区在线观看| 国产亚洲精品av在线| 在线观看午夜福利视频| 日韩欧美国产一区二区入口| 高清在线国产一区| 午夜福利免费观看在线| 久久国产精品影院| 很黄的视频免费| 国产av一区二区精品久久| 99国产精品99久久久久| 老鸭窝网址在线观看| 亚洲专区中文字幕在线| 午夜福利成人在线免费观看| 一二三四社区在线视频社区8| 精品久久久精品久久久| 禁无遮挡网站| 天天躁狠狠躁夜夜躁狠狠躁| 黄色a级毛片大全视频| 免费观看人在逋| 国产午夜福利久久久久久| 在线av久久热| 欧美激情 高清一区二区三区| 搡老熟女国产l中国老女人| 视频在线观看一区二区三区| av天堂在线播放| 在线播放国产精品三级| 国产av精品麻豆| 午夜两性在线视频| 亚洲成国产人片在线观看| 亚洲久久久国产精品| 超碰成人久久| 国产欧美日韩一区二区精品| av在线播放免费不卡| 一级毛片精品| 亚洲精华国产精华精| 欧美+亚洲+日韩+国产| 国产片内射在线| 黑丝袜美女国产一区| 亚洲中文字幕日韩| 国内精品久久久久久久电影| 免费av毛片视频| 亚洲国产日韩欧美精品在线观看 | 久久午夜综合久久蜜桃| 久久中文字幕人妻熟女| 黄网站色视频无遮挡免费观看| 女人被狂操c到高潮| 亚洲av美国av| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美午夜高清在线| 精品久久久精品久久久| 制服诱惑二区| 一a级毛片在线观看| 成人精品一区二区免费| 亚洲精品久久国产高清桃花| 国产三级在线视频| 人妻丰满熟妇av一区二区三区| 久久久久久久精品吃奶| www.自偷自拍.com| 最新美女视频免费是黄的| 免费看美女性在线毛片视频| 久久婷婷成人综合色麻豆| 国产亚洲欧美精品永久| 美女高潮到喷水免费观看| 老熟妇乱子伦视频在线观看| x7x7x7水蜜桃| 在线免费观看的www视频| 亚洲最大成人中文| 国产精品美女特级片免费视频播放器 | 国内精品久久久久久久电影| 国产精品,欧美在线| 国产熟女xx| 国产亚洲av嫩草精品影院| 自线自在国产av| 亚洲天堂国产精品一区在线| 18禁美女被吸乳视频| 18禁观看日本| 成人18禁高潮啪啪吃奶动态图| 免费女性裸体啪啪无遮挡网站| 黄色 视频免费看| 日本在线视频免费播放| 精品国产乱码久久久久久男人| 悠悠久久av| 丰满人妻熟妇乱又伦精品不卡| 99国产综合亚洲精品| 老司机靠b影院| 欧美中文综合在线视频| 丝袜美腿诱惑在线| 欧美亚洲日本最大视频资源| 欧洲精品卡2卡3卡4卡5卡区| 亚洲狠狠婷婷综合久久图片| 色哟哟哟哟哟哟| 日本三级黄在线观看| 大陆偷拍与自拍| 国产成人欧美| 国产精品永久免费网站| 99国产综合亚洲精品|