• <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)
    国产精品免费视频内射| 最新在线观看一区二区三区| 大型av网站在线播放| cao死你这个sao货| 欧美成人性av电影在线观看| 久久中文字幕一级| 国产精品精品国产色婷婷| 亚洲,欧美精品.| 亚洲成av片中文字幕在线观看| 在线永久观看黄色视频| 波多野结衣高清作品| 19禁男女啪啪无遮挡网站| 亚洲精品粉嫩美女一区| 日本免费a在线| 国产精品 国内视频| 日韩有码中文字幕| 精品免费久久久久久久清纯| 国产91精品成人一区二区三区| 国产精品98久久久久久宅男小说| 欧美性长视频在线观看| 精品少妇一区二区三区视频日本电影| 麻豆成人av在线观看| 黄网站色视频无遮挡免费观看| 午夜视频精品福利| 免费看a级黄色片| 最好的美女福利视频网| 亚洲片人在线观看| 欧美av亚洲av综合av国产av| 午夜福利一区二区在线看| 满18在线观看网站| 午夜成年电影在线免费观看| 一边摸一边做爽爽视频免费| 波多野结衣高清作品| 给我免费播放毛片高清在线观看| 黄色丝袜av网址大全| 亚洲午夜精品一区,二区,三区| 波多野结衣av一区二区av| 国产三级黄色录像| 国产黄色小视频在线观看| 免费人成视频x8x8入口观看| 给我免费播放毛片高清在线观看| 国产真实乱freesex| 可以在线观看的亚洲视频| 嫩草影院精品99| 国产高清有码在线观看视频 | 久久精品亚洲精品国产色婷小说| 欧美激情久久久久久爽电影| 亚洲国产精品合色在线| 精品国产超薄肉色丝袜足j| 久久精品人妻少妇| 美女免费视频网站| 草草在线视频免费看| 美女高潮喷水抽搐中文字幕| 国产一区二区激情短视频| 亚洲国产精品成人综合色| 老司机午夜福利在线观看视频| 国产熟女xx| av天堂在线播放| 欧美性猛交黑人性爽| 在线观看免费午夜福利视频| 黑人操中国人逼视频| 男人舔奶头视频| 欧美日韩中文字幕国产精品一区二区三区| e午夜精品久久久久久久| 看免费av毛片| 岛国在线观看网站| 亚洲一区二区三区色噜噜| 大型黄色视频在线免费观看| 亚洲,欧美精品.| 中国美女看黄片| 欧洲精品卡2卡3卡4卡5卡区| 婷婷丁香在线五月| 两个人看的免费小视频| 久久久久国产精品人妻aⅴ院| 国产精品综合久久久久久久免费| 久久天堂一区二区三区四区| 嫩草影院精品99| 午夜免费激情av| 亚洲七黄色美女视频| 午夜精品在线福利| 精品日产1卡2卡| 国产伦在线观看视频一区| 黄色成人免费大全| 91成人精品电影| 久久精品国产清高在天天线| 国产成人精品久久二区二区91| 中文字幕人妻熟女乱码| 高清在线国产一区| 国产伦一二天堂av在线观看| 日韩免费av在线播放| 黄片小视频在线播放| 亚洲欧美一区二区三区黑人| 欧美人与性动交α欧美精品济南到| 日韩欧美国产在线观看| 成年版毛片免费区| 午夜福利在线在线| 久久久久久国产a免费观看| 欧美成狂野欧美在线观看| 国产精品 欧美亚洲| 亚洲一区二区三区不卡视频| 色哟哟哟哟哟哟| 日韩欧美一区视频在线观看| 日韩av在线大香蕉| 国产97色在线日韩免费| 国产亚洲精品综合一区在线观看 | 夜夜看夜夜爽夜夜摸| 丁香欧美五月| 国产成人精品久久二区二区91| 人人妻人人澡人人看| 50天的宝宝边吃奶边哭怎么回事| 黑人巨大精品欧美一区二区mp4| 国产主播在线观看一区二区| 俄罗斯特黄特色一大片| 精品一区二区三区四区五区乱码| 国产1区2区3区精品| 国产一级毛片七仙女欲春2 | 国产激情久久老熟女| 日本熟妇午夜| 最好的美女福利视频网| 免费高清视频大片| av超薄肉色丝袜交足视频| 成人国语在线视频| 欧美中文综合在线视频| 丁香六月欧美| 两性午夜刺激爽爽歪歪视频在线观看 | 丝袜美腿诱惑在线| 在线观看免费日韩欧美大片| 在线看三级毛片| 怎么达到女性高潮| АⅤ资源中文在线天堂| 制服诱惑二区| 国产久久久一区二区三区| 欧美最黄视频在线播放免费| 日本免费a在线| 亚洲精品一区av在线观看| 午夜福利欧美成人| 一个人免费在线观看的高清视频| 青草久久国产| 久久这里只有精品19| 啦啦啦韩国在线观看视频| 国产精品乱码一区二三区的特点| 色哟哟哟哟哟哟| 亚洲,欧美精品.| 久久亚洲真实| 精品无人区乱码1区二区| 后天国语完整版免费观看| 色尼玛亚洲综合影院| 在线观看舔阴道视频| 久久国产精品影院| 97碰自拍视频| 亚洲精品色激情综合| 国产成年人精品一区二区| 色播亚洲综合网| 国产成人av教育| 亚洲成a人片在线一区二区| 国产单亲对白刺激| 国内久久婷婷六月综合欲色啪| 欧美性长视频在线观看| 久久欧美精品欧美久久欧美| 亚洲国产精品999在线| 国产精品精品国产色婷婷| 好男人在线观看高清免费视频 | 欧美中文综合在线视频| 久久精品亚洲精品国产色婷小说| 大型黄色视频在线免费观看| 两个人免费观看高清视频| а√天堂www在线а√下载| 国产精华一区二区三区| 国产97色在线日韩免费| 男男h啪啪无遮挡| 国产99白浆流出| 男人舔女人下体高潮全视频| 欧美日韩中文字幕国产精品一区二区三区| 欧美在线一区亚洲| 久久这里只有精品19| 国内精品久久久久精免费| 变态另类成人亚洲欧美熟女| 999久久久国产精品视频| 精品国产美女av久久久久小说| 俄罗斯特黄特色一大片| 香蕉国产在线看| 99热这里只有精品一区 | 国产亚洲精品久久久久久毛片| 日韩成人在线观看一区二区三区| 精品国产美女av久久久久小说| 1024视频免费在线观看| 特大巨黑吊av在线直播 | avwww免费| 免费高清视频大片| 国产精品久久视频播放| 欧美一级毛片孕妇| 午夜日韩欧美国产| 亚洲电影在线观看av| 最好的美女福利视频网| 露出奶头的视频| 国产欧美日韩一区二区精品| 在线播放国产精品三级| 久久久久国产精品人妻aⅴ院| 悠悠久久av| 国产av一区二区精品久久| 久久人妻福利社区极品人妻图片| 男人的好看免费观看在线视频 | 国产成人欧美在线观看| 伦理电影免费视频| 韩国av一区二区三区四区| 亚洲专区中文字幕在线| 久久婷婷人人爽人人干人人爱| 国内毛片毛片毛片毛片毛片| 国产亚洲欧美98| 黑人巨大精品欧美一区二区mp4| 国产一区二区三区在线臀色熟女| 男女那种视频在线观看| 在线观看一区二区三区| 午夜老司机福利片| 欧美三级亚洲精品| 日本 欧美在线| 99国产精品99久久久久| 91麻豆精品激情在线观看国产| 三级毛片av免费| 侵犯人妻中文字幕一二三四区| 精品久久久久久久毛片微露脸| 精品久久久久久久久久久久久 | 国产真实乱freesex| 久久久久国产精品人妻aⅴ院| 国产黄a三级三级三级人| 欧美另类亚洲清纯唯美| 最好的美女福利视频网| 午夜精品在线福利| 精品不卡国产一区二区三区| 欧美zozozo另类| 免费观看精品视频网站| 久久国产精品男人的天堂亚洲| 精品欧美国产一区二区三| 又黄又爽又免费观看的视频| 叶爱在线成人免费视频播放| 国产激情欧美一区二区| 午夜免费观看网址| 精品久久久久久成人av| 久久久国产欧美日韩av| 在线观看免费午夜福利视频| 久热爱精品视频在线9| 淫秽高清视频在线观看| 日日夜夜操网爽| 久久精品夜夜夜夜夜久久蜜豆 | 午夜福利成人在线免费观看| 国产男靠女视频免费网站| 国产成人影院久久av| 香蕉久久夜色| or卡值多少钱| 午夜免费鲁丝| 亚洲,欧美精品.| 成人三级做爰电影| 制服人妻中文乱码| 成人三级黄色视频| 国产日本99.免费观看| 国产99白浆流出| 人成视频在线观看免费观看| 超碰成人久久| 99国产综合亚洲精品| 人人澡人人妻人| 男女下面进入的视频免费午夜 | 国产一区二区三区视频了| 欧美zozozo另类| 国产成人精品久久二区二区免费| АⅤ资源中文在线天堂| 亚洲在线自拍视频| 日韩精品免费视频一区二区三区| 久久久久久亚洲精品国产蜜桃av| 大型黄色视频在线免费观看| 别揉我奶头~嗯~啊~动态视频| 欧美精品亚洲一区二区| 免费女性裸体啪啪无遮挡网站| 日日摸夜夜添夜夜添小说| 香蕉av资源在线| 91麻豆av在线| 男人舔女人的私密视频| 欧美日韩一级在线毛片| 国产亚洲精品久久久久久毛片| 日本a在线网址| 一级毛片精品| 欧美黄色片欧美黄色片| 欧美大码av| 日韩精品免费视频一区二区三区| 国产成人欧美在线观看| 午夜影院日韩av| 男人操女人黄网站| 精品国产亚洲在线| 一区二区三区精品91| 成年版毛片免费区| 色播亚洲综合网| 欧美乱色亚洲激情| 久久久久久大精品| 欧美大码av| 国产又色又爽无遮挡免费看| 国产v大片淫在线免费观看| 日本免费一区二区三区高清不卡| 1024手机看黄色片| 身体一侧抽搐| 午夜福利视频1000在线观看| 亚洲欧美精品综合一区二区三区| 欧美性猛交╳xxx乱大交人| 久久青草综合色| 在线永久观看黄色视频| 国内久久婷婷六月综合欲色啪| 亚洲自拍偷在线| 女警被强在线播放| 亚洲 国产 在线| 亚洲色图 男人天堂 中文字幕| 黄色成人免费大全| 麻豆av在线久日| 亚洲午夜精品一区,二区,三区| 俄罗斯特黄特色一大片| 女人爽到高潮嗷嗷叫在线视频| 99在线人妻在线中文字幕| 久久午夜综合久久蜜桃| 日韩精品免费视频一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 国产又爽黄色视频| 色av中文字幕| 国产不卡一卡二| 国产97色在线日韩免费| 精品久久久久久,| x7x7x7水蜜桃| 欧美又色又爽又黄视频| 制服诱惑二区| 一本综合久久免费| 国产亚洲欧美在线一区二区| 久久久久久久久免费视频了| 99精品欧美一区二区三区四区| 午夜成年电影在线免费观看| 免费看a级黄色片| 一个人观看的视频www高清免费观看 | 亚洲成人精品中文字幕电影| 两人在一起打扑克的视频| 免费在线观看影片大全网站| 成人手机av| 手机成人av网站| 欧美成人性av电影在线观看| 色老头精品视频在线观看| 黄色视频,在线免费观看| 男男h啪啪无遮挡| 国产精品九九99| 精品久久久久久久久久免费视频| 日本精品一区二区三区蜜桃| 久久青草综合色| 亚洲国产欧洲综合997久久, | 亚洲精品粉嫩美女一区| 免费观看精品视频网站| 亚洲精品国产精品久久久不卡| 一级毛片女人18水好多| 波多野结衣高清作品| 观看免费一级毛片| 欧美黑人欧美精品刺激| 精品久久久久久久久久免费视频| 天堂动漫精品| 看免费av毛片| 久久久久久久午夜电影| 男男h啪啪无遮挡| 国产精品日韩av在线免费观看| 亚洲在线自拍视频| 啦啦啦 在线观看视频| 脱女人内裤的视频| 久久久久久久午夜电影| 亚洲国产日韩欧美精品在线观看 | 老司机深夜福利视频在线观看| 欧美日韩亚洲综合一区二区三区_| 亚洲av片天天在线观看| 亚洲av日韩精品久久久久久密| 亚洲精品av麻豆狂野| 国产一区二区在线av高清观看| 我的亚洲天堂| 亚洲国产欧洲综合997久久, | 久久中文字幕一级| 天天一区二区日本电影三级| 国产成人精品无人区| 久久久久国产精品人妻aⅴ院| 两个人视频免费观看高清| 国产成+人综合+亚洲专区| 91av网站免费观看| 又黄又爽又免费观看的视频| a级毛片在线看网站| 国产精品一区二区精品视频观看| 色精品久久人妻99蜜桃| 亚洲国产高清在线一区二区三 | 国产精品乱码一区二三区的特点| 免费搜索国产男女视频| 在线观看免费视频日本深夜| 91国产中文字幕| 国产一级毛片七仙女欲春2 | 国产成人系列免费观看| 亚洲精品国产一区二区精华液| 国内揄拍国产精品人妻在线 | 日本五十路高清| 精品无人区乱码1区二区| 午夜免费成人在线视频| 成人国语在线视频| av欧美777| 一卡2卡三卡四卡精品乱码亚洲| 国产不卡一卡二| 欧洲精品卡2卡3卡4卡5卡区| 国产精品免费一区二区三区在线| 欧美最黄视频在线播放免费| 亚洲av电影在线进入| 在线观看舔阴道视频| 制服丝袜大香蕉在线| 成人永久免费在线观看视频| 人成视频在线观看免费观看| 熟女电影av网| 久99久视频精品免费| 精品电影一区二区在线| 国产视频内射| 久久狼人影院| 久久国产精品人妻蜜桃| 啦啦啦 在线观看视频| 成人精品一区二区免费| 国内毛片毛片毛片毛片毛片| www日本在线高清视频| 母亲3免费完整高清在线观看| 女性被躁到高潮视频| 国产高清有码在线观看视频 | 亚洲熟女毛片儿| 成人三级黄色视频| 日本 欧美在线| 99国产精品99久久久久| 欧美黑人欧美精品刺激| 午夜福利成人在线免费观看| 亚洲美女黄片视频| 人人妻人人澡人人看| 国产一区二区三区视频了| 最近最新免费中文字幕在线| 老熟妇乱子伦视频在线观看| 搡老妇女老女人老熟妇| 韩国av一区二区三区四区| videosex国产| 国产99久久九九免费精品| 亚洲精品久久成人aⅴ小说| 操出白浆在线播放| 国产一区二区三区视频了| 亚洲黑人精品在线| 欧美 亚洲 国产 日韩一| 长腿黑丝高跟| 女人被狂操c到高潮| 国内少妇人妻偷人精品xxx网站 | 在线十欧美十亚洲十日本专区| 哪里可以看免费的av片| svipshipincom国产片| 最好的美女福利视频网| 国产欧美日韩一区二区精品| 日韩欧美国产一区二区入口| 午夜福利视频1000在线观看| 好男人电影高清在线观看| 免费在线观看成人毛片| 巨乳人妻的诱惑在线观看| 精品久久久久久久久久久久久 | 在线看三级毛片| 日韩免费av在线播放| 午夜福利欧美成人| 国产av不卡久久| 国产成人啪精品午夜网站| 久久久久久人人人人人| 精品熟女少妇八av免费久了| 变态另类成人亚洲欧美熟女| 啪啪无遮挡十八禁网站| 一级黄色大片毛片| 亚洲一码二码三码区别大吗| 中文资源天堂在线| 欧美日韩中文字幕国产精品一区二区三区| 看免费av毛片| 亚洲精品色激情综合| 久久国产精品人妻蜜桃| 日韩 欧美 亚洲 中文字幕| 久久香蕉激情| 一级黄色大片毛片| 两个人视频免费观看高清| 免费在线观看日本一区| 国产又爽黄色视频| 91麻豆av在线| 18禁黄网站禁片免费观看直播| 人人妻人人澡人人看| 自线自在国产av| xxxwww97欧美| 精品电影一区二区在线| bbb黄色大片| 日韩免费av在线播放| 免费在线观看黄色视频的| 久久久久久免费高清国产稀缺| 久久九九热精品免费| 在线观看免费视频日本深夜| 99riav亚洲国产免费| 亚洲av五月六月丁香网| 啪啪无遮挡十八禁网站| 国产精品一区二区精品视频观看| 在线观看午夜福利视频| 亚洲人成网站在线播放欧美日韩| 色播亚洲综合网| 国产乱人伦免费视频| 麻豆av在线久日| 国产亚洲精品久久久久久毛片| 久久性视频一级片| 侵犯人妻中文字幕一二三四区| 日韩欧美免费精品| 欧美性长视频在线观看| 99国产精品一区二区蜜桃av| 老熟妇乱子伦视频在线观看| 最近最新中文字幕大全免费视频| 久久国产亚洲av麻豆专区| 国产黄a三级三级三级人| 人人澡人人妻人| 国产精品免费一区二区三区在线| 欧美大码av| 色综合亚洲欧美另类图片| 婷婷丁香在线五月| 国产在线观看jvid| 久久草成人影院| 日韩成人在线观看一区二区三区| 特大巨黑吊av在线直播 | cao死你这个sao货| 久久久久国产精品人妻aⅴ院| 日日干狠狠操夜夜爽| 真人做人爱边吃奶动态| 一二三四社区在线视频社区8| 窝窝影院91人妻| 国产午夜精品久久久久久| 看免费av毛片| 99在线人妻在线中文字幕| 男女那种视频在线观看| 亚洲成av片中文字幕在线观看| 怎么达到女性高潮| 亚洲精品粉嫩美女一区| 久久 成人 亚洲| 亚洲熟妇中文字幕五十中出| 99精品在免费线老司机午夜| 国产精品九九99| 热re99久久国产66热| 欧美性猛交╳xxx乱大交人| 欧美 亚洲 国产 日韩一| 国产成人精品无人区| 免费无遮挡裸体视频| 满18在线观看网站| 国产精品免费视频内射| 亚洲成av人片免费观看| 999精品在线视频| 国产麻豆成人av免费视频| 2021天堂中文幕一二区在线观 | 日韩成人在线观看一区二区三区| www国产在线视频色| 一区福利在线观看| 亚洲avbb在线观看| 亚洲激情在线av| 最新美女视频免费是黄的| 99在线人妻在线中文字幕| 亚洲av成人一区二区三| 国产精品自产拍在线观看55亚洲| 精品国产国语对白av| 亚洲欧美一区二区三区黑人| 一级a爱视频在线免费观看| 俺也久久电影网| 老熟妇仑乱视频hdxx| 最新在线观看一区二区三区| 国内揄拍国产精品人妻在线 | 亚洲 国产 在线| 狂野欧美激情性xxxx| 丰满的人妻完整版| 久久青草综合色| 国产精品av久久久久免费| 亚洲午夜精品一区,二区,三区| 一区二区三区精品91| 97碰自拍视频| 成人三级做爰电影| 国产91精品成人一区二区三区| 国产99久久九九免费精品| 黄色视频,在线免费观看| 少妇熟女aⅴ在线视频| 久久精品91无色码中文字幕| 久久久久久免费高清国产稀缺| 美女 人体艺术 gogo| 亚洲国产毛片av蜜桃av| svipshipincom国产片| 十八禁网站免费在线| 中文字幕av电影在线播放| 香蕉丝袜av| 久久香蕉精品热| 精品一区二区三区四区五区乱码| 亚洲片人在线观看| 夜夜躁狠狠躁天天躁| 免费观看人在逋| 一级黄色大片毛片| 午夜老司机福利片| 久久亚洲精品不卡| 亚洲精品av麻豆狂野| 欧美成人性av电影在线观看| 精品一区二区三区视频在线观看免费| aaaaa片日本免费| 精品国产美女av久久久久小说| 成人精品一区二区免费| 啪啪无遮挡十八禁网站| 老熟妇仑乱视频hdxx| 色av中文字幕| 变态另类丝袜制服| 日韩精品中文字幕看吧| 搡老妇女老女人老熟妇| 国产视频一区二区在线看| 精品欧美国产一区二区三| 免费看a级黄色片| 亚洲男人的天堂狠狠| 深夜精品福利| 99久久99久久久精品蜜桃| 免费人成视频x8x8入口观看| 日韩精品中文字幕看吧| 一级作爱视频免费观看| 亚洲国产毛片av蜜桃av| 黑人欧美特级aaaaaa片| 欧美性猛交╳xxx乱大交人| 精品少妇一区二区三区视频日本电影| 99久久国产精品久久久|