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

    基于細(xì)觀混凝土模型的時(shí)間逆轉(zhuǎn)損傷成像方法1)

    2017-08-12 11:57:05趙國(guó)旗仇亞萍駱英馮侃
    力學(xué)學(xué)報(bào) 2017年4期
    關(guān)鍵詞:波場(chǎng)細(xì)觀時(shí)域

    趙國(guó)旗仇亞萍 駱英 馮侃

    (江蘇大學(xué)土木工程與力學(xué)學(xué)院,江蘇鎮(zhèn)江212013)

    生物、工程及交叉力學(xué)

    基于細(xì)觀混凝土模型的時(shí)間逆轉(zhuǎn)損傷成像方法1)

    趙國(guó)旗2)仇亞萍 駱英 馮侃

    (江蘇大學(xué)土木工程與力學(xué)學(xué)院,江蘇鎮(zhèn)江212013)

    提出了一種針對(duì)混凝土結(jié)構(gòu)損傷檢測(cè)的時(shí)間逆轉(zhuǎn)損傷成像方法.以檢測(cè)混凝土結(jié)構(gòu)中與骨料尺寸相近的微小損傷為目的,引入細(xì)觀混凝土隨機(jī)骨料模型,該模型將混凝土結(jié)構(gòu)視為由水泥漿基底、骨料及粘接層組成的三相復(fù)合材料,基于Monte Carlo隨機(jī)樣本原理并結(jié)合真實(shí)試件的骨料級(jí)配曲線建立.在數(shù)值模擬分析中,將生成含損傷的細(xì)觀模型導(dǎo)入有限元分析軟件進(jìn)行超聲波場(chǎng)模擬,同時(shí)采用自適應(yīng)性強(qiáng)的時(shí)間逆轉(zhuǎn)模型(time reversedmodel,TRM)進(jìn)行損傷定位.TRM分為正向檢測(cè)和逆時(shí)成像兩個(gè)部分:正向檢測(cè)過(guò)程得到包含損傷的一系列散射回波信號(hào),從數(shù)值角度進(jìn)行時(shí)間反演并作為逆時(shí)過(guò)程的輸入信號(hào);逆時(shí)成像過(guò)程選用等效彈性參數(shù)模型,幾何尺寸與隨機(jī)骨料模型相同,時(shí)反信號(hào)在相應(yīng)幾何位置同時(shí)加載形成時(shí)反波場(chǎng),時(shí)反波場(chǎng)在損傷位置會(huì)發(fā)生干涉疊加從而導(dǎo)致能量峰值的出現(xiàn),通過(guò)確定干涉峰值時(shí)刻,并獲取該時(shí)刻對(duì)應(yīng)原始波場(chǎng)以及小波變換能量場(chǎng)完成成像.與原始數(shù)據(jù)波場(chǎng)圖相比,小波變換處理成像結(jié)果消除了雜波干擾,成像結(jié)果更加清晰.進(jìn)一步對(duì)等效彈性參數(shù)的取值進(jìn)行討論,并且在骨料尺寸范圍內(nèi)調(diào)整損傷大小,結(jié)果顯示成像結(jié)果匹配度高,對(duì)于非均質(zhì)混凝土結(jié)構(gòu)的損傷檢測(cè)能很好滿足損傷定位需求.由此證明,時(shí)間逆轉(zhuǎn)成像方法對(duì)于具有復(fù)雜結(jié)構(gòu)的混凝土材料的損傷檢測(cè)具有較好的適用性.

    混凝土結(jié)構(gòu),細(xì)觀模型,時(shí)間逆轉(zhuǎn)法,超聲檢測(cè),損傷成像

    引言

    當(dāng)前國(guó)家經(jīng)濟(jì)發(fā)展快速,基礎(chǔ)建設(shè)規(guī)模空前,從樓房、橋梁到機(jī)場(chǎng)、海港,其數(shù)量與日俱增,而混凝土正是這些基礎(chǔ)設(shè)施的重要組成部分[1].混凝土結(jié)構(gòu)內(nèi)部在生產(chǎn)和使用過(guò)程中易產(chǎn)生損傷,損傷的積累將導(dǎo)致結(jié)構(gòu)失效,從而誘發(fā)災(zāi)難性事故.無(wú)損檢測(cè)技術(shù)成為目前預(yù)防災(zāi)難性事故的重要手段.

    先前的損傷檢測(cè)研究將混凝土視為均質(zhì)材料,在有限元軟件模型庫(kù)中采用平均化參數(shù)處理方法;然而當(dāng)結(jié)構(gòu)中損傷的尺寸與骨料相近時(shí),平均化處理會(huì)引起檢測(cè)精度降低.基于細(xì)觀力學(xué)為背景的細(xì)觀混凝土模型已經(jīng)取得重大成果,細(xì)觀模型從細(xì)觀尺度出發(fā)將混凝土結(jié)構(gòu)視為由骨料、基底及粘接層組成的復(fù)合材料.在早期研究中以Wittmann等[2]的二維模型以及Schlangen等[3]的格構(gòu)模型模擬了混凝土不同組分的力學(xué)特性.近年來(lái)混凝土細(xì)觀分析的基礎(chǔ)是生成形狀、尺寸、位置分布都相似于真實(shí)混凝土粗骨料的隨機(jī)骨料結(jié)構(gòu),文獻(xiàn)[4-5]分別對(duì)圓形骨料及橢圓形骨料的生成方法和投放技巧進(jìn)行了報(bào)道;Wang等[6]較早提出隨機(jī)骨料模型并進(jìn)行有限元網(wǎng)格的劃分;馬懷發(fā)等[7]依據(jù)Walaraven公式生成二維二級(jí)配和四級(jí)配混凝土隨機(jī)骨料模型;楊華等[8]對(duì)隨機(jī)骨料模型的混凝土彈性模量的確定進(jìn)行了預(yù)測(cè)性研究.本文研究將細(xì)觀隨機(jī)骨料模型引入無(wú)損檢測(cè)以提高檢測(cè)精度,結(jié)合超聲檢測(cè)技術(shù)完成損傷定位成像.

    近年來(lái),超聲檢測(cè)技術(shù)由于具有能量高、穿透力強(qiáng)、檢測(cè)速度快以及安全方便等優(yōu)點(diǎn),逐漸成為國(guó)內(nèi)外研究機(jī)構(gòu)關(guān)于混凝土無(wú)損檢測(cè)的研究重點(diǎn),并得到廣泛應(yīng)用.超聲波用于混凝土檢測(cè)始于1949年,Leslie和Cheesman[9]首次采用超聲脈沖檢測(cè)混凝土結(jié)構(gòu)缺陷;然而由于實(shí)驗(yàn)局限性,諸多影響因素在當(dāng)時(shí)未得到驗(yàn)證,因此未能用于實(shí)際工程檢測(cè).21世紀(jì)以來(lái),隨著電子技術(shù)的發(fā)展,超聲檢測(cè)技術(shù)進(jìn)入數(shù)字時(shí)代,數(shù)值仿真模擬技術(shù)的突破促使超聲成像技術(shù)得到進(jìn)一步發(fā)展,其中合成孔徑聚焦技術(shù)[10-11]、超聲相控陣技術(shù)[1213]、層析成像技術(shù)[14-15]等都取得了一定成果.時(shí)間逆轉(zhuǎn)技術(shù)由于其自適應(yīng)性特點(diǎn)逐漸由地球科學(xué)領(lǐng)域推廣到無(wú)損檢測(cè)領(lǐng)域(nondestructive testing,NDT).邱雷等[16-18]針對(duì)復(fù)合材料板結(jié)構(gòu)檢測(cè)中存在的信號(hào)信噪比低、頻散以及模式混疊問(wèn)題提出了時(shí)間反轉(zhuǎn)聚焦成像法,有效提高了有用信號(hào)能量.Zhu等[1921]對(duì)Lamb的時(shí)逆聚焦技術(shù)提出了基于頻率--波數(shù)域內(nèi)的快速成像法.近年來(lái),將時(shí)間逆轉(zhuǎn)技術(shù)應(yīng)用于混凝土檢測(cè)的嘗試逐漸引起學(xué)者關(guān)注.時(shí)間逆轉(zhuǎn)模型(time reversedmodeling,TRM)是探測(cè)地球物理學(xué)中一項(xiàng)波源定位技術(shù)[22].Saenger等[2326]將TRM技術(shù)應(yīng)用于多相高散射混凝土復(fù)合材料的檢測(cè),TRM被用于對(duì)聲發(fā)射源進(jìn)行定位和定性,與之前的檢測(cè)技術(shù)相比無(wú)需進(jìn)行首時(shí)識(shí)別或預(yù)知檢測(cè)通道的相關(guān)信息.

    為了彌補(bǔ)先前研究中將混凝土視為均質(zhì)材料而造成的仿真實(shí)驗(yàn)實(shí)用性不足的缺陷,本文依據(jù)真實(shí)試件的骨料級(jí)配曲線建立更具有實(shí)際應(yīng)用價(jià)值的細(xì)觀混凝土模型,結(jié)合有限元軟件進(jìn)行模擬,采用TRM技術(shù)完成損傷檢測(cè).同時(shí)對(duì)比不同成像條件結(jié)果的精確性,提出更適用于實(shí)際混凝土結(jié)構(gòu)檢測(cè)的損傷成像方法.

    1 時(shí)間逆轉(zhuǎn)成像法

    時(shí)間逆轉(zhuǎn)成像法包含正向檢測(cè)和逆時(shí)成像兩個(gè)步驟:正向檢測(cè)過(guò)程完成損傷信號(hào)的采集,采用單點(diǎn)激勵(lì)信號(hào)得到正向波場(chǎng)及損傷散射波場(chǎng),通過(guò)邊界傳感器接收散射回波信號(hào);逆時(shí)成像過(guò)程是從數(shù)值角度將損傷散射信號(hào)在時(shí)域內(nèi)反轉(zhuǎn)并在相同位置加載,依據(jù)波前干涉疊加確定損傷位置.

    1.1 正向檢測(cè)

    確定檢測(cè)區(qū)域Ω,邊界表示為?Ω,以x∈Ω表示檢測(cè)區(qū)域中的位置.檢測(cè)時(shí)域確定為t∈[0,T],T為終止時(shí)間.正向檢測(cè)過(guò)程用來(lái)獲取邊界一系列位移場(chǎng)時(shí)域信號(hào)

    式中,x(k)代表傳感器位置,k為傳感器編號(hào).

    式中,N為邊界傳感器的個(gè)數(shù)[25].若加載點(diǎn)與損傷間的距離為L(zhǎng),傳感器與損傷間的距離為L(zhǎng)k,傳感器接收到損傷散射信號(hào)的時(shí)間為tk=(L+Lk)/vp,其中vp為縱波波速.

    1.2 逆時(shí)成像

    逆時(shí)過(guò)程采用等效彈性參數(shù)模型[2627],模型參數(shù)由正向檢測(cè)過(guò)程實(shí)際測(cè)試得到,包括縱波波速vp,eff、橫波波速vs,eff以及密度ρeff.正向檢測(cè)過(guò)程中邊界?Ω記錄的位移信號(hào)為

    式中,t∈[0,T].將位移信號(hào)在時(shí)域內(nèi)取反

    時(shí)反信號(hào)將作為逆時(shí)成像過(guò)程的激勵(lì)信號(hào).正向檢測(cè)過(guò)程記錄的k組信號(hào)所對(duì)應(yīng)的幾何位置將作為時(shí)反信號(hào)加載的位置,k組信號(hào)同時(shí)加載形成時(shí)反波場(chǎng).信號(hào)中包含的損傷散射信號(hào)擁有同一個(gè)波源即損傷,該損傷又稱二次波源.時(shí)反處理后,傳感器加載的損傷散射信號(hào)到達(dá)損傷位置時(shí)間為

    將tk=(L+Lk)/vp代入式(5)得tRk=T-L/vp,由此可見,時(shí)反過(guò)程中損傷散射信號(hào)回到損傷位置所用的時(shí)間與各個(gè)時(shí)反加載點(diǎn)幾何位置無(wú)關(guān),僅與正向加載點(diǎn)與損傷的幾何位置以及波速有關(guān),則k組損傷信號(hào)將同時(shí)回到波源,形成波前干涉,造成能量疊加.

    1.3 成像條件

    為了更好地呈現(xiàn)時(shí)反波場(chǎng)中的波前干涉疊加現(xiàn)象,引入成像條件進(jìn)行干涉情況表征和損傷位置確定.全時(shí)域全場(chǎng)位移最大值成像[28]

    最大位移值umax(x)為時(shí)反波場(chǎng)中每一點(diǎn)x∈Ω在t∈[0,T]內(nèi)的位移最大值.此外,另提取時(shí)反應(yīng)力場(chǎng),依據(jù)全時(shí)域應(yīng)力最大值成像

    式中,最大應(yīng)力值σmax(x)為時(shí)反波場(chǎng)中每一點(diǎn)x∈Ω在t∈[0,T]內(nèi)的應(yīng)力最大值.由于時(shí)反加載過(guò)程僅選取來(lái)自邊界數(shù)量有限位置上的信號(hào),受雜波干擾會(huì)造成難以識(shí)別損傷信號(hào)形成的匯聚點(diǎn),引入表征物理特性總和的能量場(chǎng)成像法[28]

    式中,在全時(shí)域內(nèi)網(wǎng)格點(diǎn)處讀取應(yīng)力、應(yīng)變,文獻(xiàn)[29]證明了能量法對(duì)于混凝土中聲發(fā)射源的定位最為精確.

    然而針對(duì)本文中混凝土損傷的主動(dòng)檢測(cè)問(wèn)題,由于正向檢測(cè)過(guò)程中提取的損傷散射信號(hào)比直達(dá)波、邊界反射波微弱,加上混凝土結(jié)構(gòu)對(duì)波形衰減嚴(yán)重[30]以及邊界臨近傳感器加載時(shí)會(huì)形成波形混疊等因素,全時(shí)域內(nèi)成像結(jié)果并不理想.若能獲取時(shí)反波場(chǎng)中損傷散射信號(hào)的干涉疊加產(chǎn)生峰值時(shí)刻對(duì)應(yīng)的波場(chǎng),便能顯著提高對(duì)損傷位置的確定精度.因此,本文提出單時(shí)刻波場(chǎng)成像法,為了準(zhǔn)確定位能量聚焦的位置,對(duì)波場(chǎng)中各檢測(cè)點(diǎn)采用小波變換法提取中心頻率對(duì)應(yīng)的能量幅值,并確定對(duì)應(yīng)的最大幅值出現(xiàn)時(shí)刻te,max,同時(shí)獲取該時(shí)刻時(shí)反波場(chǎng)以及小波變換能量場(chǎng)

    成像定位誤差用下式進(jìn)行表征

    式中,?L表示定位位置與真實(shí)損傷之間的距離,Rmax為混凝土結(jié)構(gòu)中骨料的最大粒徑.

    2 細(xì)觀混凝土模型的建立

    從專注層面以及研究方法等側(cè)重點(diǎn)的不同,研究過(guò)程中將混凝土結(jié)構(gòu)分為3個(gè)尺度:宏觀尺度(macro-scale)、細(xì)觀尺度(meso-scale)以及微觀尺度(m icro-scale)[31].當(dāng)檢測(cè)損傷尺寸與骨料尺寸相近時(shí),選取細(xì)觀角度分析更為精確.

    混凝土骨料的分布是隨機(jī)的,依據(jù)Monte Carlo[32]法進(jìn)行隨機(jī)數(shù)的生成.生成的骨料采用經(jīng)典最大密實(shí)度曲線指導(dǎo)尺寸以及各尺寸對(duì)應(yīng)的數(shù)量分布

    其中,D0為篩選直徑,p為經(jīng)過(guò)D0篩選的骨料所占的質(zhì)量百分比,Dmax為骨料最大粒徑.根據(jù)上式,Walraven推出二維截面內(nèi)D<D0粒徑骨料出現(xiàn)的概率[33]計(jì)算公式

    式中,pk為骨料體積占總體積百分比.據(jù)此可得出橫截面上骨料分布情況.本研究采用MATLAB軟件生成具有隨機(jī)骨料尺寸和隨機(jī)骨料分布的混凝土模型.

    2.2 模型生成實(shí)例

    依據(jù)實(shí)際混凝土試件,骨料粒徑為10mm~40mm,水泥富余系數(shù)為1.08,水、水泥、砂和石子之間的比例為0.42:1:1.152:2.449;混凝土試件整體尺寸大小為400mm×500mm×500mm,據(jù)此試件相關(guān)參數(shù)建立混凝土剖面--損傷--傳感器二維截面模型如圖1所示.

    圖1 細(xì)觀混凝土仿真模型Fig.1 Meso-scale concretemodel

    將模型導(dǎo)入有限元軟件進(jìn)行參數(shù)設(shè)定,模型尺寸為500mm×400mm,設(shè)定損傷位置為圓形通孔損傷,中心位置在(225,265)mm處,半徑為r=30mm.假設(shè)傳感器為點(diǎn)源,坐標(biāo)(150,400)mm至(300,400)mm每間隔12.5mm設(shè)置一個(gè)傳感器,編號(hào)為1~13.模型主要材料參數(shù)如表1所示.

    表1 混凝土仿真模型主要參數(shù)Table 1 Main propertiesof the concretemodel

    3 數(shù)值仿真及損傷成像

    圖2 激勵(lì)信號(hào)Fig.2 Excication signal

    3.1 仿真成像

    (2)因?yàn)樗訟E=OF=FB′.所以EO=EB′.因而點(diǎn)O在⊙E上.所以從而∠BAC=∠BED.所以DE∥CA.

    將建立的混凝土隨機(jī)骨料模型導(dǎo)入COMSOL Multiphysics 5.1有限元軟件建立如圖1所示的混凝土剖面--損傷--傳感器模型,選用7號(hào)傳感器作為激勵(lì)源進(jìn)行單點(diǎn)激發(fā)多點(diǎn)接收模式.激勵(lì)信號(hào)為Hanning窗調(diào)制的五周期正弦信號(hào),如圖2(a)所示,以等效力形式加載;圖2(b)為激勵(lì)信號(hào)對(duì)應(yīng)的頻譜.由于縱波波長(zhǎng)最大為27mm,為保證數(shù)值計(jì)算的從而出現(xiàn)匯聚點(diǎn).

    但由于正向檢測(cè)過(guò)程中提取的損傷散射信號(hào)較直達(dá)波、邊界反射波微弱,再加上混凝土結(jié)構(gòu)對(duì)波形衰減嚴(yán)重以及邊界臨近傳感器加載時(shí)信號(hào)會(huì)在邊界及路徑中出現(xiàn)不同形式的混疊,損傷位置之外在邊界以及路徑中也出現(xiàn)了像素亮點(diǎn),致使全時(shí)域內(nèi)成像結(jié)果并不理想.因此,本文提出基于小波變換的單時(shí)刻波場(chǎng)成像法,選取波形干涉出現(xiàn)能量峰值時(shí)刻的全場(chǎng)波形來(lái)進(jìn)行損傷成像,所選基函數(shù)為Gabor小波,該基函數(shù)可以同時(shí)提供時(shí)域和頻域局部化的信息[34].穩(wěn)定性與準(zhǔn)確性,故設(shè)置模型網(wǎng)格單元尺寸不大于3×10-4m,時(shí)間步長(zhǎng)為0.5μs.

    激勵(lì)源進(jìn)行信號(hào)激勵(lì)形成正向波場(chǎng)和損傷散射回波,全部傳感器采集信號(hào).波形中主波包主要由直達(dá)波、損傷散射回波、邊界反射波組成,根據(jù)波速以及模型幾何尺寸關(guān)系,截取僅包含以直達(dá)波和損傷散射回波為主波形的信號(hào)(t∈[0μs,150μs]).

    逆時(shí)成像過(guò)程選用等效彈性參數(shù)模型,幾何尺寸與隨機(jī)骨料模型相同,采用軟件ANSYS 14.0建立,依據(jù)文獻(xiàn)[26-27]進(jìn)行等效參數(shù)的獲取,縱波波速vp,eff=4030.5m/s、橫波波速vs,eff=2586.3m/s以及密度ρeff=2600kg/m3.網(wǎng)格尺寸為2mm,波場(chǎng)傳播時(shí)間與信號(hào)時(shí)間窗值均為T=150μs.時(shí)反處理后的13組信號(hào)在相應(yīng)幾何位置同時(shí)加載形成時(shí)反波場(chǎng).分別得到時(shí)反位移場(chǎng)、時(shí)反應(yīng)力場(chǎng)和時(shí)反能量場(chǎng)信息.

    圖3(a)~圖3(c)分別為依據(jù)式(6)~式(8)的成像結(jié)果,圖中每個(gè)像素點(diǎn)的值為該點(diǎn)上對(duì)應(yīng)全時(shí)域內(nèi)的場(chǎng)量最大值;由于時(shí)反信號(hào)中損傷回波形成的波包在同時(shí)刻匯聚于損傷位置形成波前干涉疊加,

    圖3 全時(shí)域成像Fig.3 Whole time domain imaging

    基于小波變換的峰值時(shí)刻波場(chǎng)成像法的主要步驟如下:(1)首先對(duì)波場(chǎng)中各檢測(cè)點(diǎn)采用小波變換提取加載主頻率f0=150kHz對(duì)應(yīng)的信號(hào)幅值,基于小波變換后的全場(chǎng)數(shù)據(jù)進(jìn)行最大幅值出現(xiàn)時(shí)刻te,max的識(shí)別.如圖4(a)所示,橫坐標(biāo)表示時(shí)間,時(shí)域T=150μs,步長(zhǎng)為0.5μs;縱坐標(biāo)表示每一時(shí)刻對(duì)應(yīng)全場(chǎng)的能量最大值,最大峰值左側(cè)出現(xiàn)的第一個(gè)峰值是由波形在邊界處的混疊造成,全場(chǎng)最大能量出現(xiàn)時(shí)刻te,max=109.2μs;(2)獲取該時(shí)刻對(duì)應(yīng)原始數(shù)據(jù)波場(chǎng)圖.本文中為了提高計(jì)算效率采用間距2mm的網(wǎng)格點(diǎn),成像過(guò)程為提高精度采用插值處理.如圖4(b)所示,圖中顯示了波形在該時(shí)刻的干涉疊加,同時(shí)也出現(xiàn)了雜波干擾,圖像聚焦點(diǎn)對(duì)應(yīng)實(shí)際坐標(biāo)(223.9,289.5)mm;(3)對(duì)峰值時(shí)刻小波變換數(shù)據(jù)進(jìn)行成像,如圖4(c)所示,圖像聚焦點(diǎn)對(duì)應(yīng)實(shí)際坐標(biāo)(224.8,291.0)mm.

    在正向檢測(cè)過(guò)程中損傷回波信號(hào)是由于損傷外表面反射引起,由結(jié)構(gòu)邊界特定位置布置的傳感器采集得到.損傷表面被認(rèn)為是二次波源,時(shí)反過(guò)程中信號(hào)在波源處干涉疊加,實(shí)際損傷反射點(diǎn)坐標(biāo)為(225,295)mm;檢測(cè)結(jié)果與實(shí)際損傷位置相差4.005mm,根據(jù)式(10)進(jìn)行誤差表征,誤差為10.01%,成像結(jié)果準(zhǔn)確地匹配了損傷上表面位置,與原始數(shù)據(jù)波場(chǎng)圖相比,小波變換后的數(shù)據(jù)成像結(jié)果消除了雜波干擾,成像結(jié)果更加清晰.

    圖4 能量場(chǎng)成像圖Fig.4 Image of theenergy fiel

    圖5 不同等效參數(shù)成像結(jié)果Fig.5 Damage imaging of di ff erentequivalentparameters

    3.2 關(guān)于等效參數(shù)取值及成像精度的討論

    原始模型為非均質(zhì)結(jié)構(gòu),與實(shí)際結(jié)構(gòu)相比,均勻化等效彈性參數(shù)的取值可能引起損傷定位誤差.現(xiàn)依據(jù)文獻(xiàn)[23]的方法討論等效波速取值對(duì)損傷定位精度的影響,對(duì)波速等效分別設(shè)置為:(1)vp=vp,eff+100m/s,vs=vs,eff+100m/s;(2)vp=vp,eff-100m/s,vs=vs,eff-100m/s.

    圖5(a)為設(shè)置一的模擬結(jié)果,時(shí)反波場(chǎng)數(shù)據(jù)經(jīng)小波變換后能量峰值出現(xiàn)時(shí)刻為te,max=114.5μs,該時(shí)刻成像聚焦點(diǎn)坐標(biāo)為(222.7,278.9)mm,實(shí)際損傷位置間距為16.1mm,據(jù)式(10)計(jì)算誤差為40.25%;圖5(b)為設(shè)置二的模擬結(jié)果,能量峰值出現(xiàn)時(shí)刻為te,max=111.8μs,該時(shí)刻成像聚焦點(diǎn)坐標(biāo)為(221.7,297.4)mm,實(shí)際損傷位置間距為4.08mm,誤差為10.2%.結(jié)果顯示,損傷定位點(diǎn)并沒(méi)有出現(xiàn)大的偏移,邊界傳感器一定程度上平衡了參數(shù)取值的差異性誤差.

    本文研究的時(shí)間逆轉(zhuǎn)損傷成像方法以檢測(cè)骨料粒徑尺寸范圍內(nèi)的損傷為目的,適用于“不密實(shí)區(qū)”及“空洞”類型缺陷的檢測(cè);此類缺陷一般是由于石子架空現(xiàn)象造成,即損傷尺寸與骨料石子尺寸相近,屬于微小型缺陷.

    本文選用的混凝土結(jié)構(gòu)骨料粒徑為10mm~40mm,故縮小損傷直徑d為40mm和20mm,成像結(jié)果見圖6.圖6(a)焦點(diǎn)坐標(biāo)為(226.1,274.1)mm,對(duì)應(yīng)損傷實(shí)際坐標(biāo)為(225,285)mm,間距為10.95mm,誤差為27.3%.圖6(b)焦點(diǎn)坐標(biāo)為(224.1,260.1)mm,對(duì)應(yīng)損傷實(shí)際坐標(biāo)為(225,275)mm,間距為15.02mm,誤差為37.5%.

    從圖6可看出,損傷成像結(jié)果精度隨損傷尺寸的減小而降低,當(dāng)損傷尺寸小于骨料平均尺寸(27.63mm)時(shí),成像焦點(diǎn)識(shí)別度降低(圖6(b)),當(dāng)損傷尺寸小于骨料最小尺寸時(shí),成像結(jié)果將無(wú)法識(shí)別.然而對(duì)于非均質(zhì)混凝土結(jié)構(gòu)中“空洞”型損傷的檢測(cè),本文的定位結(jié)果已滿足損傷定位需求.由此證明了時(shí)間逆轉(zhuǎn)模型損傷檢測(cè)法及基于小波變換的單時(shí)刻成像方法對(duì)于具有復(fù)雜結(jié)構(gòu)的混凝土材料的損傷檢測(cè)方法依然具有較好的適用性,且操作簡(jiǎn)便、結(jié)果清晰.

    圖6 不同損傷尺寸成像結(jié)果Fig.6 Damage imaging of di ff erentdiameters

    4 結(jié)論

    本文基于細(xì)觀混凝土模型,結(jié)合時(shí)間逆轉(zhuǎn)成像法進(jìn)行損傷檢測(cè);通過(guò)正向檢測(cè)和逆時(shí)成像兩個(gè)過(guò)程完成損傷回波信號(hào)的提取和時(shí)反加載.時(shí)反加載在建立的與細(xì)觀模型相同尺寸的等效模型中進(jìn)行,損傷回波信號(hào)在損傷處發(fā)生干涉疊加;通過(guò)提取時(shí)反過(guò)程波場(chǎng)數(shù)據(jù),并進(jìn)一步結(jié)合小波變換法提取信號(hào)主頻率對(duì)應(yīng)的能量幅值,最終對(duì)基于小波變換后的全場(chǎng)數(shù)據(jù)進(jìn)行最大幅值出現(xiàn)時(shí)刻的識(shí)別并獲取該時(shí)刻對(duì)應(yīng)原始波場(chǎng)以及小波變換能量場(chǎng),完成成像,成像結(jié)果清晰且成功匹配了損傷表面的幾何位置.數(shù)值計(jì)算結(jié)果表明,時(shí)間逆轉(zhuǎn)模型對(duì)復(fù)雜混凝土結(jié)構(gòu)的損傷定位具有強(qiáng)的適用性,基于小波變換的單時(shí)刻成像法能夠消除噪聲干擾提高成像精度.本文采用依據(jù)真實(shí)試件的骨料級(jí)配曲線建立細(xì)觀混凝土模型提高了研究的實(shí)際應(yīng)用價(jià)值,提出的適用于復(fù)雜混凝土結(jié)構(gòu)的損傷定位方法,對(duì)混凝土結(jié)構(gòu)內(nèi)部損傷檢測(cè)的工程實(shí)踐具有較強(qiáng)的指導(dǎo)及應(yīng)用價(jià)值.

    1萬(wàn)征,姚仰平,孟達(dá).復(fù)雜加載下混凝土的彈塑性本構(gòu)模型.力學(xué)學(xué)報(bào),2016,48(5):1159-1171(Wan Zheng,Yao Yangping,Meng Da.An elastoplastic constitutivemodel of concrete under complicated load.Chinese JournalofTheoreticaland Applied Mechanics,2016,48(5):1159-1171(in Chinese))

    2 Wittmann FH,Roelfstra PE.Simulation and analysis of composite structures.Materials Science and Engineering,1984,68(2):239-248

    3 Schlangen E,Van M ier JGM.Simple latticemodel for numerical simulation of fracture of concretematerials and structures.Materialsand Structures,1992,25(9):534-542

    4 W riggers P,Moftah SO.Mesoscalemodels for concrete:homogenization and damage behaviour.Finite Elements in Analysis and Design,2006,42(7):623-636

    5 Leite JPB,Slow ik V,Apel J.Computationalmodel ofmesoscopic structure of concrete for simulation of fracture processes.Comput-ersand Structures,2007,85(17/18):1293-1303

    6 Wang ZM,Kwan AKH,Chan HC.Mesoscopic study of concrete I:generation of random aggregate structure and finit elementmesh.Computers&Structures,1999,70(5):533-544

    7馬懷發(fā),陳厚群,黎保琨.混凝土試件細(xì)觀結(jié)構(gòu)的數(shù)值模擬.水利學(xué)報(bào),2004,35(10):27-35(Ma Huaifa,Chen Houqun,Li Baokun.Meso-structure numerical simulation of concrete specimens.JournalofHydraulic Engineering,2004,35(10):27-35(in Chinese))

    8楊華,李宗利,惠弘毅.基于隨機(jī)骨料模型的混凝土彈性模量預(yù)測(cè)研究.長(zhǎng)江科學(xué)院院報(bào),2016,33(2):100-105(Yang Hua,Li Zongli,HuiHongyi.Prediction ofelasticmodulusof concretebased on random aggregatemodel.JournalofYangtze River Scientifi Research Institute,2016,33(2):100-105(in Chinese))

    9 Leslie JR,CheesmanWJ.An ultrasonicmethod of studying deterioration and cracking in concrete structures.Journalofthe American Concrete Institute,1949,21(1):17-36

    10 Clayton D,Barker A,A lbright A,et al.Improved synthetic aperture focusing technique resultsof thick concrete specimens through frequency banding//Chimenti DE,Bond LJ.AIP Conference Proceedings.AIPPublishing,2016,1706(1):020012

    11 FendtKT,Mooshofer H,Rupitsch SJ,etal.Ultrasonic defect characterization in heavy rotor forgingsbymeansof the synthetic aperture focusing technique and optimizationmethods.IEEE Transactions on Ultrasonics,Ferroelectrics,and Frequency Control,2016,63(6):874-885

    12 BeniwalS,Ghosh D,GanguliA.Ultrasonic imaging of concreteusing scattered elastic wave modes.NDT&E International,2016,82(1):26-35

    13 Liu Z,Sun K,Song G,etal.Damage localization in alum inum plate w ith compact rectangularphased piezoelectric transducerarray.MechanicalSystemsand SignalProcessing,2016,70(1):625-636

    14 Vassholz M,Koberstein-Schwarz B,Ruhlandt A,etal.New X-ray tomography method based on the 3d Radon transform compatible w ith anisotropic sources.Physical Review Letters,2016,116(8):088101-1

    15李飛,余西龍,林鑫.基于TDLAS的層析成像技術(shù)TDLAT.力學(xué)學(xué)報(bào),2014,46(1):55-59(LiFei,Yu Xilong,Lin Xin.Parallel-beam tomography based on TDLAS.Chinese Journal ofTheoretical and Applied Mechanics,2014,46(1):55-59(in Chinese))

    16邱雷,袁慎芳,蘇永振.基于Shannon復(fù)數(shù)小波和時(shí)間反轉(zhuǎn)聚焦的復(fù)合材料結(jié)構(gòu)多源沖擊成像定位方法.航空學(xué)報(bào),2010,31(12):2417-2424(Qiu Lei,Yuan Shenfang.Multiple impactsource imaging and localization on composite structure based on shannon complex waveletand time reversal focusing.Acta Aeronautica EtAstronautica Sinica,2010,31(12):2417-2424(in Chinese))

    17王強(qiáng),袁慎芳.主動(dòng)Lamb波結(jié)構(gòu)健康監(jiān)測(cè)中信號(hào)增強(qiáng)與損傷成像方法.航空學(xué)報(bào),2008,29(4):1061-1067(Wang Qiang,Yuan Shenfang.Amplifying signal and imaging damagemethod for active lamb-wave structure health monitoring.Acta Aeronautica Et Astronautica Sinica,2008,29(4):1061-1067(in Chinese))

    18王強(qiáng),袁慎芳.無(wú)參考主動(dòng)Lamb波結(jié)構(gòu)損傷時(shí)反成像監(jiān)測(cè)方法.航空學(xué)報(bào),2010,31(1):178-183(Wang Qiang,Yuan Shenfang.No baseline time reversal imagingmethod for active lamb wave structuraldamagemonitoring.Acta Aeronautica Et Astronautica Sinica,2010,31(1):178-183(in Chinese))

    19 Zhu R,Huang GL,Yuan FG.Fastdamage imaging using the timereversal technique in the frequency-wavenumber domain.Smart Materialsand Structures,2013,22(7):075028

    20 He J,Harb MS,Yuan FG.Non-contact ultrasonic technique for Lambwave characterization in composite plates.Ultrasonics,2016,64(1):162-169

    21 Chen C,Yuan FG.Impact source identificatio in finit isotropic platesusing a time-reversalmethod:theoreticalstudy.SmartMaterialsand Structures,2010,19(10):105028

    22 Fink M,Cassereau D,Derode A,et al.Time-reversed acoustics.Reportson Progress in Physics,2000,63(12):1933

    23 Saenger EH.Time reverse characterization of sources in heterogeneousmedia.NDT&E International,2011,44(8):751-759

    24 Kocur GK,Saenger EH.Deconvolution of acoustic em issions for source localization using time reversemodeling.Journal ofSound and Vibration,2017,387(1):66-78

    25 Saenger EH,Kocur GK,Jud R,et al.Application of time reverse modeling on ultrasonic non-destructive testing of concrete.Applied MathematicalModelling,2011,35(2):807-816

    26 SaengerEH.Numericalmethods to determinee ff ectiveelastic properties.International Journal ofEngineering Science,2008,46(6):598-605

    27 Saenger EH,Shapiro SA.E ff ective velocities in fracturedmedia:a numerical study using the rotated staggered finite-d ff erence grid.Geophysical Prospecting,2002,50(2):183-194

    28 Steiner B,Saenger EH,Schmalholz SM.Time reversemodeling of low-frequency m icrotremors:application to hydrocarbon reservoir localization.GeophysicalResearch Letters,2008,35(3):L03307

    29 Kocur GK,Saenger EH,Grosse CU,et al.Time reversemodeling of acoustic em issions in a reinforced concrete beam.Ultrasonics,2016,65(1):96-104

    30蔡偉,陳文.復(fù)雜介質(zhì)中任意階頻率依賴耗散聲波的分?jǐn)?shù)階導(dǎo)數(shù)模型.力學(xué)學(xué)報(bào),2016,48(6):1265-1280(CaiWei,ChenWen.Fractuionalderivativemodeling of frequency-dependentdissipative mechanism forwave propagation in complexmedia.Chinese JournalofTheoreticaland Applied Mechanics,2016,48(6):1265-1280(in Chinese))

    31唐春安,朱萬(wàn)成.混凝土損傷與斷裂一數(shù)值試驗(yàn).北京:科學(xué)出版社,2003(Tang Chunan,ZhuWancheng.A Numerical Experiment of Concrete Damage and Fracture.Beijing:Science Press,2003(in Chinese))

    32 Doucet A,Freitas DN,Gordon N.An introduction to sequential MonteCarlomethods.Statistics forEngineering&Information Science,2001,18(1):3-14

    33 Walraven JC,Reinhardt HW.Theory and experiments on themechanical behaviour of cracks in plain and reinforced concrete subjected to shear loading.HERON,1981,26(1A):1-68

    34 Liu Y,Li Z,ZhangW.Crack detection of fibr reinforced composite beams based on continuouswavelet transform.Nondestructive Testing&Evaluation,2010,25(25):25-44

    DAMAGEDETECTIONOFMESO-SCALECONCRETE STRUCTURESBASED ON TIME REVERSALMETHOD1)

    Zhao Guoqi2)Qiu Yaping Luo Ying Feng Kan
    (Faculty ofCivil Engineering and MechanicsofJiangsu University,Zhenjiang 212013,Jiangsu,China)

    A specifi time reversal imaging method is proposed in this article to detect defects in concrete structures.In order to detect the damage that the scale is the same as the aggregate,ameso-scale concretemodel is introduced in thisarticle.As the concrete isa compositematerial composed by cement,aggregate,waterand concrete admixtures,the Monte Carlo random modeland the aggregate grading curve of real concrete samples are introduced for designing this finit elementmodel.Then,thedamagedmodelwasanalyzed by employing aself-adapted time reversedmodel to achieve the ultrasonicwave fiel simulation.This imagingmethod contains two steps:the firs is the forward detection.A series of reflecte echo signalsw ith damage information are obtained in this section.These

    signals can be reversed in Matlab to serve as the incident signals in the next Time Reversal process;the second step is to image the damage location via interfering thewave-frontsactuated by di ff erent transducers to illustrate the peaks ofwaveform amplitudes.By determ ining the interfering wave peak time,obtaining the originalwave fiel of thatmoment,we can form ing the wavelet transform energy field and then complete the damage imaging of the concretemodelw ith defect.In the time reversalprocess,we introduced theequivalentelastic parametersas the same geometric dimension of the originalmesosacleconcretemodel to locate thedamagee ff ectively.Comparedw ith theoriginalwave field theenergy fiel modifie by wavelet transform can lower thee ff ectof theenvironmentalnoise.Finally,we discussed theequivalentelastic parameters and the damage sizes to verify the robustnessof thismethod which isapplicable inmonitoring and evaluating the damage in concrete structures.

    concrete structure,meso-scalemodel,time reversalmethod,ultrasonic detection,damage imaging

    TU317+.8,O347.4+1

    A

    10.6052/0459-1879-17-007

    2017-01-04收稿,2017-03-23錄用,2017-03-24網(wǎng)絡(luò)版發(fā)表.

    1)國(guó)家自然科學(xué)重點(diǎn)國(guó)際合作項(xiàng)目(11520101001)、江蘇大學(xué)高級(jí)人才啟動(dòng)基金(5501480007)資助.

    2)趙國(guó)旗,副教授,主要研究方向:結(jié)構(gòu)健康監(jiān)測(cè)與檢測(cè)新技術(shù)研究及應(yīng)用.E-mail:flagzhao@163.co

    趙國(guó)旗,仇亞萍,駱英,馮侃.基于細(xì)觀混凝土模型的時(shí)間逆轉(zhuǎn)損傷成像方法.力學(xué)學(xué)報(bào),2017,49(4):953-960

    Zhao Guoqi,Qiu Yaping,Luo Ying,Feng Kan.Damage detection ofmeso-scale concrete structures based on time reversalmethod.Chinese JournalofTheoreticaland Applied Mechanics,2017,49(4):953-960

    猜你喜歡
    波場(chǎng)細(xì)觀時(shí)域
    基于細(xì)觀結(jié)構(gòu)的原狀黃土動(dòng)彈性模量和阻尼比試驗(yàn)研究
    地震研究(2021年1期)2021-04-13 01:05:24
    基于時(shí)域信號(hào)的三電平逆變器復(fù)合故障診斷
    彈性波波場(chǎng)分離方法對(duì)比及其在逆時(shí)偏移成像中的應(yīng)用
    交錯(cuò)網(wǎng)格與旋轉(zhuǎn)交錯(cuò)網(wǎng)格對(duì)VTI介質(zhì)波場(chǎng)分離的影響分析
    基于Hilbert變換的全波場(chǎng)分離逆時(shí)偏移成像
    基于極大似然準(zhǔn)則與滾動(dòng)時(shí)域估計(jì)的自適應(yīng)UKF算法
    基于時(shí)域逆濾波的寬帶脈沖聲生成技術(shù)
    基于時(shí)域波形特征的輸電線雷擊識(shí)別
    旋轉(zhuǎn)交錯(cuò)網(wǎng)格VTI介質(zhì)波場(chǎng)模擬與波場(chǎng)分解
    基于四叉樹網(wǎng)格加密技術(shù)的混凝土細(xì)觀模型
    男女之事视频高清在线观看 | 日韩 亚洲 欧美在线| 亚洲专区中文字幕在线| 亚洲国产精品999| 久久人妻熟女aⅴ| 女警被强在线播放| av国产精品久久久久影院| 亚洲国产最新在线播放| 成人18禁高潮啪啪吃奶动态图| 男女之事视频高清在线观看 | 国产一卡二卡三卡精品| 欧美日韩亚洲高清精品| 国产精品一二三区在线看| 50天的宝宝边吃奶边哭怎么回事| 国产成人av教育| 免费看不卡的av| 精品亚洲乱码少妇综合久久| 午夜视频精品福利| 亚洲人成网站在线观看播放| 男女边摸边吃奶| 高清av免费在线| 老司机深夜福利视频在线观看 | 女人久久www免费人成看片| 中文欧美无线码| 久久久久精品国产欧美久久久 | 亚洲精品成人av观看孕妇| 在线观看一区二区三区激情| 狠狠婷婷综合久久久久久88av| 午夜老司机福利片| 亚洲精品av麻豆狂野| 天堂俺去俺来也www色官网| 久久毛片免费看一区二区三区| 永久免费av网站大全| 一本色道久久久久久精品综合| 热re99久久国产66热| 999久久久国产精品视频| 一二三四社区在线视频社区8| 超色免费av| 欧美国产精品一级二级三级| 久久久精品区二区三区| 亚洲精品美女久久久久99蜜臀 | 精品亚洲成a人片在线观看| 午夜福利免费观看在线| 国产成人av教育| 91精品国产国语对白视频| 晚上一个人看的免费电影| 黄片小视频在线播放| 1024视频免费在线观看| 久久久久精品国产欧美久久久 | 国产高清视频在线播放一区 | 欧美日韩一级在线毛片| 99久久综合免费| 伊人久久大香线蕉亚洲五| 19禁男女啪啪无遮挡网站| 亚洲熟女精品中文字幕| 久久 成人 亚洲| 激情视频va一区二区三区| 国产精品国产av在线观看| 天天躁夜夜躁狠狠久久av| 日韩av免费高清视频| 亚洲一区中文字幕在线| 国产又色又爽无遮挡免| 久久精品亚洲av国产电影网| 亚洲精品av麻豆狂野| 三上悠亚av全集在线观看| www.av在线官网国产| 最黄视频免费看| 久久精品久久精品一区二区三区| 丝瓜视频免费看黄片| 国产三级黄色录像| 亚洲免费av在线视频| 一本综合久久免费| 国产精品麻豆人妻色哟哟久久| 成年美女黄网站色视频大全免费| 亚洲国产欧美在线一区| 一级黄片播放器| 免费看不卡的av| 嫁个100分男人电影在线观看 | 色网站视频免费| 两人在一起打扑克的视频| 国语对白做爰xxxⅹ性视频网站| 女人高潮潮喷娇喘18禁视频| 成人亚洲欧美一区二区av| 国产成人91sexporn| 午夜日韩欧美国产| 国产精品久久久av美女十八| 男女下面插进去视频免费观看| 日韩一卡2卡3卡4卡2021年| 精品一区在线观看国产| 久久久久精品国产欧美久久久 | 亚洲图色成人| 精品久久久精品久久久| 大片免费播放器 马上看| 国产成人a∨麻豆精品| 久久国产亚洲av麻豆专区| 精品卡一卡二卡四卡免费| 在线 av 中文字幕| 亚洲情色 制服丝袜| svipshipincom国产片| 国产亚洲欧美精品永久| 可以免费在线观看a视频的电影网站| 丝袜美足系列| 一区福利在线观看| 欧美日韩视频精品一区| av电影中文网址| 欧美精品高潮呻吟av久久| 亚洲一码二码三码区别大吗| 一级毛片电影观看| 天堂俺去俺来也www色官网| 亚洲精品美女久久av网站| av福利片在线| 欧美 亚洲 国产 日韩一| 美女中出高潮动态图| 老鸭窝网址在线观看| 国产日韩一区二区三区精品不卡| 老司机靠b影院| 国产精品九九99| 久久久精品国产亚洲av高清涩受| 丝袜脚勾引网站| 午夜免费男女啪啪视频观看| 亚洲欧美精品自产自拍| 中国美女看黄片| 亚洲中文av在线| av有码第一页| 丝袜人妻中文字幕| 日日摸夜夜添夜夜爱| 亚洲精品一区蜜桃| 电影成人av| 欧美久久黑人一区二区| 国产日韩欧美亚洲二区| 又大又爽又粗| 王馨瑶露胸无遮挡在线观看| 免费看不卡的av| 国产熟女午夜一区二区三区| www日本在线高清视频| 人妻 亚洲 视频| 大香蕉久久成人网| 婷婷色综合www| 成年动漫av网址| 18在线观看网站| 校园人妻丝袜中文字幕| 国产亚洲av片在线观看秒播厂| 男女下面插进去视频免费观看| 国产人伦9x9x在线观看| 成年av动漫网址| 欧美激情高清一区二区三区| 国产精品一区二区免费欧美 | 午夜免费男女啪啪视频观看| 日本vs欧美在线观看视频| 久久久久网色| 欧美黑人欧美精品刺激| 少妇人妻久久综合中文| 亚洲欧美激情在线| 美女脱内裤让男人舔精品视频| 精品久久蜜臀av无| 欧美精品亚洲一区二区| 妹子高潮喷水视频| 天天躁夜夜躁狠狠躁躁| 欧美精品亚洲一区二区| 亚洲精品第二区| 久久久久久久大尺度免费视频| 精品人妻1区二区| 日韩制服丝袜自拍偷拍| 亚洲av日韩精品久久久久久密 | 亚洲欧美激情在线| 免费观看人在逋| 国产精品欧美亚洲77777| 免费看av在线观看网站| 最近最新中文字幕大全免费视频 | 日韩熟女老妇一区二区性免费视频| 十分钟在线观看高清视频www| 777久久人妻少妇嫩草av网站| 人妻 亚洲 视频| 性色av乱码一区二区三区2| 午夜av观看不卡| 国产成人欧美在线观看 | 热99国产精品久久久久久7| 丰满饥渴人妻一区二区三| 天堂俺去俺来也www色官网| 大码成人一级视频| 飞空精品影院首页| 色精品久久人妻99蜜桃| 黄色毛片三级朝国网站| 午夜福利乱码中文字幕| 亚洲国产av影院在线观看| 成年人午夜在线观看视频| 亚洲男人天堂网一区| 亚洲av男天堂| 老汉色∧v一级毛片| 99精国产麻豆久久婷婷| 青青草视频在线视频观看| 亚洲欧美一区二区三区国产| 韩国精品一区二区三区| 男人添女人高潮全过程视频| 三上悠亚av全集在线观看| 亚洲国产精品国产精品| 18禁国产床啪视频网站| 天天操日日干夜夜撸| 丝袜人妻中文字幕| 亚洲一区二区三区欧美精品| 久久精品久久久久久噜噜老黄| 大码成人一级视频| 国产高清视频在线播放一区 | av福利片在线| 国产老妇伦熟女老妇高清| 免费看不卡的av| 99香蕉大伊视频| 国产精品一区二区在线不卡| 午夜福利视频精品| 精品国产国语对白av| 2021少妇久久久久久久久久久| 久久国产精品影院| 精品国产超薄肉色丝袜足j| 狠狠精品人妻久久久久久综合| 日本91视频免费播放| 欧美变态另类bdsm刘玥| 尾随美女入室| 久久99精品国语久久久| 欧美性长视频在线观看| 亚洲中文av在线| 夫妻午夜视频| 国产精品久久久av美女十八| 国产成人av教育| 国产黄色视频一区二区在线观看| 国产伦人伦偷精品视频| 精品高清国产在线一区| 欧美人与善性xxx| 可以免费在线观看a视频的电影网站| 亚洲精品在线美女| 日韩av不卡免费在线播放| 91精品伊人久久大香线蕉| 亚洲人成网站在线观看播放| 国产精品人妻久久久影院| 老汉色∧v一级毛片| 亚洲国产精品一区二区三区在线| 热re99久久国产66热| kizo精华| 午夜久久久在线观看| 久久人人97超碰香蕉20202| 亚洲激情五月婷婷啪啪| av国产久精品久网站免费入址| 欧美老熟妇乱子伦牲交| 国产精品免费大片| 少妇精品久久久久久久| 亚洲av成人精品一二三区| 少妇 在线观看| 99精国产麻豆久久婷婷| 国产成人免费无遮挡视频| 考比视频在线观看| 亚洲综合色网址| 91字幕亚洲| 最新的欧美精品一区二区| 婷婷色综合www| 老汉色av国产亚洲站长工具| 久久影院123| 午夜影院在线不卡| 七月丁香在线播放| 嫩草影视91久久| 欧美亚洲 丝袜 人妻 在线| av电影中文网址| 久久精品国产亚洲av涩爱| 国产精品免费大片| 午夜福利视频在线观看免费| 国产国语露脸激情在线看| 亚洲精品第二区| 狠狠婷婷综合久久久久久88av| 看免费av毛片| kizo精华| 一级毛片电影观看| 一区福利在线观看| xxx大片免费视频| 欧美+亚洲+日韩+国产| 亚洲精品久久久久久婷婷小说| 精品久久久精品久久久| 午夜福利在线免费观看网站| 亚洲,一卡二卡三卡| 亚洲人成网站在线观看播放| www.精华液| 极品少妇高潮喷水抽搐| 黄色片一级片一级黄色片| 国产精品国产三级国产专区5o| 日本a在线网址| 99久久综合免费| 亚洲成av片中文字幕在线观看| 国产激情久久老熟女| 国产野战对白在线观看| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品久久成人aⅴ小说| 国产成人一区二区三区免费视频网站 | 国产精品麻豆人妻色哟哟久久| 又大又黄又爽视频免费| 在线亚洲精品国产二区图片欧美| 两性夫妻黄色片| 精品福利观看| 国产成人精品久久二区二区91| 99九九在线精品视频| 在线天堂中文资源库| av天堂久久9| 日韩 欧美 亚洲 中文字幕| 午夜老司机福利片| 好男人电影高清在线观看| 制服诱惑二区| 啦啦啦啦在线视频资源| 国产成人影院久久av| 手机成人av网站| 香蕉丝袜av| 91精品国产国语对白视频| 欧美人与善性xxx| 亚洲欧洲精品一区二区精品久久久| 日韩视频在线欧美| 操美女的视频在线观看| 日韩中文字幕欧美一区二区 | 亚洲国产欧美在线一区| 日日爽夜夜爽网站| 女人被躁到高潮嗷嗷叫费观| 国产伦理片在线播放av一区| 国产麻豆69| 久久精品久久久久久久性| 亚洲黑人精品在线| 97人妻天天添夜夜摸| 男女免费视频国产| netflix在线观看网站| 九草在线视频观看| 亚洲欧美日韩高清在线视频 | av在线老鸭窝| 日韩大片免费观看网站| 精品人妻1区二区| 18禁黄网站禁片午夜丰满| 老汉色∧v一级毛片| 久久99精品国语久久久| 激情视频va一区二区三区| 这个男人来自地球电影免费观看| 国产精品免费视频内射| 久久久久精品国产欧美久久久 | 69精品国产乱码久久久| 大片电影免费在线观看免费| 国产片特级美女逼逼视频| 国产免费又黄又爽又色| 日韩av不卡免费在线播放| 国产成人av教育| 精品国产国语对白av| 新久久久久国产一级毛片| av有码第一页| 日韩中文字幕视频在线看片| av电影中文网址| 精品视频人人做人人爽| 日韩一本色道免费dvd| 久久久国产一区二区| 大片免费播放器 马上看| 在线观看一区二区三区激情| 久久99精品国语久久久| 天天躁日日躁夜夜躁夜夜| 18禁国产床啪视频网站| bbb黄色大片| 亚洲欧美日韩另类电影网站| 人人妻人人爽人人添夜夜欢视频| 精品久久久久久电影网| 黄片播放在线免费| 999精品在线视频| 在线天堂中文资源库| 一本色道久久久久久精品综合| 女性被躁到高潮视频| 国产免费又黄又爽又色| 午夜福利免费观看在线| 2018国产大陆天天弄谢| 亚洲欧美精品综合一区二区三区| 一个人免费看片子| 嫁个100分男人电影在线观看 | 欧美黑人精品巨大| 成年人黄色毛片网站| 精品少妇一区二区三区视频日本电影| 视频在线观看一区二区三区| 又紧又爽又黄一区二区| 久久精品亚洲av国产电影网| 悠悠久久av| 国产熟女午夜一区二区三区| 自线自在国产av| 亚洲av成人精品一二三区| 久久人人97超碰香蕉20202| 黄色片一级片一级黄色片| 好男人电影高清在线观看| 精品国产乱码久久久久久小说| 久久久久久久国产电影| 色婷婷av一区二区三区视频| tube8黄色片| 中文字幕色久视频| 男人添女人高潮全过程视频| www.熟女人妻精品国产| 欧美老熟妇乱子伦牲交| 久久精品久久精品一区二区三区| 成人亚洲精品一区在线观看| 一本大道久久a久久精品| 国产免费福利视频在线观看| 在现免费观看毛片| 视频在线观看一区二区三区| 亚洲精品国产区一区二| 少妇人妻 视频| 精品人妻一区二区三区麻豆| 色综合欧美亚洲国产小说| 一区福利在线观看| 亚洲国产精品一区三区| 精品人妻熟女毛片av久久网站| 久久综合国产亚洲精品| 岛国毛片在线播放| 人妻人人澡人人爽人人| cao死你这个sao货| 久久人人爽人人片av| 超色免费av| 国产亚洲一区二区精品| 精品熟女少妇八av免费久了| 男女午夜视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 国产亚洲精品久久久久5区| 日本欧美国产在线视频| 电影成人av| 亚洲欧美激情在线| 欧美日韩福利视频一区二区| 精品人妻1区二区| 成年女人毛片免费观看观看9 | 国产成人免费观看mmmm| 国产淫语在线视频| 国产人伦9x9x在线观看| 啦啦啦 在线观看视频| 真人做人爱边吃奶动态| 久久人妻熟女aⅴ| 9热在线视频观看99| 日韩一本色道免费dvd| 韩国高清视频一区二区三区| 国产国语露脸激情在线看| 精品亚洲成a人片在线观看| 婷婷丁香在线五月| 亚洲视频免费观看视频| 日日爽夜夜爽网站| 你懂的网址亚洲精品在线观看| 老司机影院成人| 一级,二级,三级黄色视频| 久久久精品区二区三区| 久久久精品国产亚洲av高清涩受| 巨乳人妻的诱惑在线观看| 黑人猛操日本美女一级片| 1024香蕉在线观看| 亚洲,一卡二卡三卡| 欧美日韩综合久久久久久| 2021少妇久久久久久久久久久| 伊人亚洲综合成人网| 久久久精品区二区三区| 夫妻性生交免费视频一级片| 国产精品免费大片| 深夜精品福利| 精品少妇内射三级| 亚洲一码二码三码区别大吗| 亚洲成人免费电影在线观看 | 丝袜喷水一区| 国产精品欧美亚洲77777| 男女床上黄色一级片免费看| 9191精品国产免费久久| 国产精品久久久久久人妻精品电影 | 一级片'在线观看视频| 午夜福利一区二区在线看| 9色porny在线观看| 国产高清不卡午夜福利| 男人爽女人下面视频在线观看| 国产黄色免费在线视频| 女人被躁到高潮嗷嗷叫费观| 国产成人欧美在线观看 | 免费一级毛片在线播放高清视频 | 岛国毛片在线播放| av不卡在线播放| 亚洲av国产av综合av卡| 亚洲欧美精品综合一区二区三区| 国产福利在线免费观看视频| 9191精品国产免费久久| 欧美日本中文国产一区发布| 日韩av免费高清视频| av天堂在线播放| 一区二区三区精品91| 中文乱码字字幕精品一区二区三区| 国产成人精品在线电影| 脱女人内裤的视频| 久久国产精品大桥未久av| 一边摸一边抽搐一进一出视频| 91国产中文字幕| 成人国产一区最新在线观看 | 国产一卡二卡三卡精品| 大片免费播放器 马上看| 亚洲一码二码三码区别大吗| 亚洲精品国产区一区二| 成人国语在线视频| 国产又色又爽无遮挡免| 亚洲中文日韩欧美视频| 亚洲黑人精品在线| 午夜视频精品福利| 啦啦啦在线观看免费高清www| 亚洲第一青青草原| 国产成人精品久久二区二区免费| 精品福利永久在线观看| 午夜视频精品福利| 成人三级做爰电影| 午夜视频精品福利| av又黄又爽大尺度在线免费看| a级片在线免费高清观看视频| 美女福利国产在线| 男的添女的下面高潮视频| 国产无遮挡羞羞视频在线观看| 国产麻豆69| 久久精品aⅴ一区二区三区四区| 免费在线观看日本一区| 国产亚洲精品久久久久5区| 国产亚洲av高清不卡| 少妇人妻 视频| 99久久人妻综合| 中文字幕人妻熟女乱码| 热re99久久精品国产66热6| 国产精品国产三级专区第一集| 黄色视频在线播放观看不卡| 久久人妻福利社区极品人妻图片 | 国产欧美日韩一区二区三 | 久久午夜综合久久蜜桃| 成在线人永久免费视频| 国产视频一区二区在线看| 久久这里只有精品19| 人妻 亚洲 视频| 免费一级毛片在线播放高清视频 | 色网站视频免费| 一二三四在线观看免费中文在| 亚洲人成电影免费在线| av欧美777| 丁香六月欧美| 欧美在线一区亚洲| 18禁观看日本| 国产精品欧美亚洲77777| 性高湖久久久久久久久免费观看| 777久久人妻少妇嫩草av网站| 日韩伦理黄色片| 美女午夜性视频免费| 免费在线观看日本一区| 永久免费av网站大全| 亚洲国产成人一精品久久久| av网站在线播放免费| 成年人免费黄色播放视频| 亚洲精品日本国产第一区| 麻豆国产av国片精品| 精品一区二区三卡| 久久 成人 亚洲| 美女大奶头黄色视频| 999久久久国产精品视频| 午夜久久久在线观看| 日韩熟女老妇一区二区性免费视频| 777久久人妻少妇嫩草av网站| 国产国语露脸激情在线看| 人妻 亚洲 视频| 色播在线永久视频| 亚洲人成电影观看| 一个人免费看片子| 亚洲精品久久午夜乱码| 国产女主播在线喷水免费视频网站| 中文字幕亚洲精品专区| 国产男女内射视频| 午夜福利免费观看在线| 99热全是精品| 亚洲成国产人片在线观看| 久久亚洲精品不卡| av天堂久久9| 国产一区二区三区综合在线观看| 国产av国产精品国产| 国产亚洲欧美在线一区二区| av国产精品久久久久影院| 尾随美女入室| 久久综合国产亚洲精品| 老汉色av国产亚洲站长工具| 精品少妇久久久久久888优播| 在线精品无人区一区二区三| 国产av国产精品国产| 黄色视频不卡| 91精品国产国语对白视频| 人人妻人人澡人人爽人人夜夜| 午夜日韩欧美国产| 亚洲精品自拍成人| 午夜福利视频在线观看免费| 一级毛片女人18水好多 | 亚洲中文字幕日韩| xxx大片免费视频| 亚洲中文av在线| 久久人人爽人人片av| 国产精品亚洲av一区麻豆| 老汉色av国产亚洲站长工具| 国产精品99久久99久久久不卡| 性少妇av在线| 在线精品无人区一区二区三| 手机成人av网站| 性少妇av在线| 亚洲成av片中文字幕在线观看| a级毛片黄视频| 国产精品秋霞免费鲁丝片| 黄片小视频在线播放| 亚洲精品日本国产第一区| 亚洲av成人精品一二三区| 熟女av电影| 2021少妇久久久久久久久久久| 日韩伦理黄色片| a级片在线免费高清观看视频| 99国产精品一区二区三区| 99精品久久久久人妻精品| 另类亚洲欧美激情| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产av新网站| 宅男免费午夜| 久久久久网色| 菩萨蛮人人尽说江南好唐韦庄| 亚洲久久久国产精品| 下体分泌物呈黄色| 桃花免费在线播放| 欧美激情极品国产一区二区三区| 国产视频一区二区在线看| 色综合欧美亚洲国产小说|