999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

鈦合金深層裂紋缺陷渦流檢測(cè)仿真與實(shí)驗(yàn)研究?

2021-06-16 10:35:38張程杰胡明慧徐小雄華東理工大學(xué)承壓與安全教育部重點(diǎn)實(shí)驗(yàn)室上海200237
傳感技術(shù)學(xué)報(bào) 2021年3期
關(guān)鍵詞:檢測(cè)

張程杰胡明慧徐小雄(華東理工大學(xué)承壓與安全教育部重點(diǎn)實(shí)驗(yàn)室,上海200237)

為了檢測(cè)、定位、評(píng)價(jià)缺陷,評(píng)估構(gòu)件的結(jié)構(gòu)完整性,無(wú)損檢測(cè)通過不損害被測(cè)對(duì)象的未來(lái)用途和功能的方式,對(duì)原材料和零部件進(jìn)行檢測(cè)[1-2]。渦流檢測(cè)(Eddy Current Testing,ET)基于電磁感應(yīng)原理,獲取導(dǎo)電材料的表面或近表面的質(zhì)量信息。其檢測(cè)具有無(wú)需耦合介質(zhì)、檢測(cè)速度快、對(duì)表面和近表面裂紋缺陷的檢測(cè)靈敏度高等優(yōu)點(diǎn)[3-5],因此在航空航天、核電、鐵路、建筑橋梁等領(lǐng)域都有著廣泛的應(yīng)用。

近年來(lái),渦流檢測(cè)在航空航天無(wú)損檢測(cè)中的應(yīng)用越來(lái)越廣泛,在檢測(cè)飛機(jī)蒙皮和多層結(jié)構(gòu)的隱藏缺陷工作中,渦流檢測(cè)的占比已經(jīng)超過50%[6]。在飛機(jī)結(jié)構(gòu)件材料的渦流檢測(cè)研究方面,曾輝耀[7]等針對(duì)單向碳纖維增強(qiáng)復(fù)合材料分層缺陷,設(shè)計(jì)了一個(gè)共面三矩形線圈,通過有限元仿真研究實(shí)現(xiàn)了分層缺陷的檢測(cè)和定量化評(píng)估。孟杞鳳[8]等針對(duì)碳纖維增強(qiáng)復(fù)合材料深層缺陷,開發(fā)了一種激勵(lì)線圈與檢測(cè)線圈相互垂直的新型渦流探頭,并通過仿真與實(shí)驗(yàn)研究證明了探頭能有效檢測(cè)碳纖維復(fù)合材料埋深層2 mm~4 mm處缺陷。在飛機(jī)關(guān)鍵部件渦流檢測(cè)方面,Michal Janovec[9]等通過諧波渦流激勵(lì)陣列探頭,對(duì)飛機(jī)機(jī)身與機(jī)翼鉚接接頭可能存在的裂紋或隱藏腐蝕進(jìn)行檢測(cè)。Wang[10]等應(yīng)用激勵(lì)線圈和隧道磁阻傳感器(TMR)陣列的方式,在100 Hz和1 kHz頻率激勵(lì)下,實(shí)現(xiàn)了對(duì)飛機(jī)鈦合金蒙皮多層鉚接結(jié)構(gòu)深埋缺陷的檢測(cè)。

鈦合金材料具有密度小、強(qiáng)度高、耐蝕、耐高溫、無(wú)磁、可焊、使用溫度范圍寬等優(yōu)點(diǎn)[11],因此廣泛應(yīng)用于航空發(fā)動(dòng)機(jī)、飛機(jī)機(jī)身、航空緊固件等方面。然而目前常規(guī)渦流檢測(cè)方法無(wú)法檢測(cè)出鈦合金深層裂紋缺陷,因此,本文以鈦合金TC4板件深層微小裂紋缺陷為研究對(duì)象,設(shè)計(jì)了一種基于TMR陣列的差分式雙矩形線圈激勵(lì)渦流檢測(cè)探頭。通過建立三維仿真模型,模擬探頭移動(dòng)過程,得到不同尺寸深層缺陷的鈦合金板件表面的磁場(chǎng)圖像,分析了該探頭對(duì)3 mm~12 mm不同長(zhǎng)度和0.2 mm~2 mm不同寬度尺寸的深層缺陷的檢測(cè)能力。研制實(shí)驗(yàn)探頭,搭建實(shí)驗(yàn)系統(tǒng),通過實(shí)驗(yàn)與仿真結(jié)果進(jìn)行比較驗(yàn)證TMR渦流探頭的檢測(cè)性能,為鈦合金材料的深層缺陷渦流檢測(cè)探頭設(shè)計(jì)提供了參考。

1 探頭結(jié)構(gòu)及其原理

1.1 探頭結(jié)構(gòu)設(shè)計(jì)

電磁感應(yīng)原理表明,當(dāng)導(dǎo)體處于交變磁場(chǎng)中時(shí),內(nèi)部會(huì)感應(yīng)出形成閉合回路的渦流。當(dāng)導(dǎo)體中存在裂紋時(shí),導(dǎo)體材料不連續(xù)導(dǎo)致裂紋處導(dǎo)電性能發(fā)生變化,渦流出現(xiàn)擾動(dòng),從而導(dǎo)致渦流產(chǎn)生的二次感應(yīng)磁場(chǎng)發(fā)生變化[12-13],進(jìn)而可直接由TMR傳感器測(cè)得磁場(chǎng)變化并判斷缺陷是否存在。

圖1 TMR傳感器R-H響應(yīng)曲線[10]

對(duì)于TMR傳感器,在很小的外磁場(chǎng)作用下,其磁隧道結(jié)(MTJs)結(jié)構(gòu)中的鐵磁層磁化方向就會(huì)發(fā)生變化,從而實(shí)現(xiàn)其穿隧電阻的巨大變化[14-15]。TMR傳感器的R-H響應(yīng)曲線如圖1所示,其線性高靈敏度區(qū)域在零外場(chǎng)附近,當(dāng)外加磁場(chǎng)強(qiáng)度高于靈敏度區(qū)域時(shí),R-H曲線線性度下降,TMR傳感器的檢測(cè)靈敏度會(huì)受到極大影響。

因此,本文所設(shè)計(jì)的檢測(cè)探頭由兩個(gè)激勵(lì)線圈和TMR傳感器陣列組成,如圖2所示。其中,激勵(lì)線圈為對(duì)稱布置的差分式雙線圈結(jié)構(gòu),兩個(gè)激勵(lì)線圈結(jié)構(gòu)尺寸、匝數(shù)等參數(shù)均相同,線圈串聯(lián)接入激勵(lì)電路。TMR傳感器陣列為測(cè)量元件,由3個(gè)TMR傳感器線性排列構(gòu)成,其靈敏度為25 mV/V/Oe,工作電壓為5 V,由此可以通過輸出信號(hào)反算得到試件表面磁感應(yīng)強(qiáng)度大小:

式中:Bz為試件表面磁感應(yīng)強(qiáng)度;Sen為傳感器靈敏度;U為輸出信號(hào)幅值;U0為傳感器工作電壓。

圖2 探頭結(jié)構(gòu)

根據(jù)畢奧-薩伐爾定律,線圈在空間某點(diǎn)處產(chǎn)生的磁感應(yīng)強(qiáng)度與電流的大小成正比:

式中:I為線圈電流;er為電流元指向觀測(cè)點(diǎn)的單位矢量;μ0為真空磁導(dǎo)率;r為電流源點(diǎn)與觀測(cè)點(diǎn)之間的距離。在TMR傳感器所在觀測(cè)點(diǎn)處,當(dāng)傳感器處于均勻介質(zhì)中時(shí),兩個(gè)激勵(lì)線圈通以不同方向的電流,產(chǎn)生大小相等,方向相反的磁場(chǎng),從而使得兩線圈產(chǎn)生的磁感應(yīng)強(qiáng)度在傳感器陣列所在的中心線上相互抵消,確保TMR陣列處于高靈敏度的狀態(tài)。當(dāng)缺陷進(jìn)入線圈下方的檢測(cè)區(qū)域時(shí),兩側(cè)線圈產(chǎn)生的磁場(chǎng)不能完全抵消,故而在傳感器所在中心線上表現(xiàn)為磁感應(yīng)強(qiáng)度的變化。

2 缺陷電磁檢測(cè)仿真分析

2.1 深層缺陷三維仿真模型

本文采用COMSOL Multiphysics建立深層缺陷檢測(cè)仿真模型,在AC/DC模塊下選擇磁場(chǎng)為物理場(chǎng),在頻域下進(jìn)行求解。模型網(wǎng)格劃分如圖3所示。

圖3 模型網(wǎng)格劃分

模型中仿真試件材料為TC4,相對(duì)介電常數(shù)和相對(duì)磁導(dǎo)率為1,電導(dǎo)率為6×105S/m,試件的長(zhǎng)、寬、高分別為120 mm、120 mm和5 mm??諝夂腿毕莸碾妼?dǎo)率設(shè)為0.1,相對(duì)介電常數(shù)和相對(duì)磁導(dǎo)率均設(shè)為1。缺陷設(shè)置于試件背面,其尺寸根據(jù)仿真需求確定。激勵(lì)線圈設(shè)置為差分式雙激勵(lì)線圈,單個(gè)線圈外尺寸為14 mm×84 mm,內(nèi)尺寸為4 mm×74 mm,線圈高度10 mm,線徑0.4 mm,匝數(shù)200匝,兩線圈間距離為2 mm。激勵(lì)電流設(shè)置為0.2 A,兩個(gè)激勵(lì)線圈中的電流方向相反,當(dāng)探頭在無(wú)缺陷區(qū)域時(shí),兩激勵(lì)線圈在中線處產(chǎn)生的感應(yīng)磁場(chǎng)相互抵消,中線處磁感應(yīng)強(qiáng)度為零。線圈結(jié)構(gòu)尺寸參數(shù)與后續(xù)實(shí)驗(yàn)驗(yàn)證所用線圈參數(shù)相同。

2.2 激勵(lì)頻率

根據(jù)標(biāo)準(zhǔn)透入深度公式可得不同激勵(lì)頻率下的標(biāo)準(zhǔn)透入深度理論值,對(duì)于非鐵磁性材料磁導(dǎo)率,有μ≈μ0=4π×10-7H/m,則標(biāo)準(zhǔn)透入深度為[16]

式中:δ為標(biāo)準(zhǔn)透入深度,f為激勵(lì)頻率,μ為材料磁導(dǎo)率,σ為材料電導(dǎo)率??芍?對(duì)于同一待測(cè)試件,渦流標(biāo)準(zhǔn)透入深度只與激勵(lì)頻率有關(guān)。

為研究探頭檢測(cè)深層缺陷的最佳激勵(lì)頻率范圍,通過理論計(jì)算不同激勵(lì)頻率下的標(biāo)準(zhǔn)透入深度與在該激勵(lì)頻率下仿真得到的缺陷引起的磁感應(yīng)強(qiáng)度變化值比較。設(shè)置1 kHz~30 kHz不等的激勵(lì)頻率檢測(cè)距試件表面4 mm深度下缺陷,缺陷的長(zhǎng)、寬、高分別為6 mm、0.2 mm、1 mm,得到不同激勵(lì)頻率下缺陷引起的磁感應(yīng)強(qiáng)度變化Bz分量如表1所示。

表1 不同激勵(lì)頻率下缺陷引起的磁感應(yīng)強(qiáng)度變化

由此可知,距表面4 mm下缺陷引起的磁感應(yīng)強(qiáng)度在10 kHz~30 kHz激勵(lì)頻率均有較大的變化,在20 kHz左右達(dá)到最大值;20 kHz的標(biāo)準(zhǔn)透入深度也接近4 mm。因此,可以通過標(biāo)準(zhǔn)透入深度公式初步確定檢測(cè)不同深度缺陷的最佳激勵(lì)頻率范圍。

2.3 不同長(zhǎng)度缺陷的仿真分析

仿真研究探頭對(duì)不同長(zhǎng)度缺陷檢測(cè)的檢出效果,設(shè)置缺陷位于試件表面下3 mm深度,缺陷寬度和高度分別為0.2 mm、2 mm,長(zhǎng)度分別為3 mm、6 mm、8 mm、10 mm、12 mm,激勵(lì)頻率10 kHz,得到不同長(zhǎng)度缺陷的磁感應(yīng)強(qiáng)度Bz分量分布的仿真結(jié)果如圖4所示。

圖4 不同長(zhǎng)度(3 mm~12 mm)缺陷仿真結(jié)果

由圖4可知,仿真結(jié)果能直接反映缺陷的長(zhǎng)度變化,缺陷長(zhǎng)度越長(zhǎng),則仿真結(jié)果中磁感應(yīng)強(qiáng)度變化的范圍對(duì)應(yīng)增長(zhǎng)。

磁感應(yīng)強(qiáng)度Bz分量峰值與缺陷長(zhǎng)度的關(guān)系如圖5所示。從圖5中可以看出,在缺陷長(zhǎng)度6 mm以下時(shí),磁感應(yīng)強(qiáng)度Bz分量隨著缺陷長(zhǎng)度的增加而增大;當(dāng)缺陷長(zhǎng)度超過6 mm后,繼續(xù)增加缺陷長(zhǎng)度,磁感應(yīng)強(qiáng)度Bz分量增長(zhǎng)減緩。這是由于線圈繞線寬度為5 mm,繼續(xù)增加缺陷長(zhǎng)度時(shí),缺陷長(zhǎng)度超出線圈寬度,缺陷長(zhǎng)度變化對(duì)線圈產(chǎn)生的渦流的擾動(dòng)效果影響減弱。

圖5 磁感應(yīng)強(qiáng)度Bz幅值和缺陷長(zhǎng)度的關(guān)系

2.4 不同寬度缺陷的仿真分析

仿真研究探頭對(duì)不同寬度缺陷檢測(cè)的檢出效果,設(shè)置缺陷位于試件表面下3 mm深度,根據(jù)上節(jié)內(nèi)容確定缺陷長(zhǎng)度和高度分別為6 mm、2 mm,寬度分別為0.2 mm、0.4 mm、0.8 mm、1.2 mm、1.6 mm、2.0 mm,激勵(lì)頻率10 kHz,得到不同寬度缺陷的磁感應(yīng)強(qiáng)度Bz分量分布的仿真結(jié)果如圖6所示。

從圖6可以看出,不同寬度缺陷的仿真結(jié)果表現(xiàn)為磁感應(yīng)強(qiáng)度Bz分量的強(qiáng)弱變化,難以直接判定缺陷的寬度。這是由于缺陷寬度尺寸較小,而單個(gè)傳感器所掃查的寬度大于缺陷寬度。因此無(wú)法通過Bz分量的變化直接定量判定缺陷的寬度變化。

圖6 不同寬度(0.2 mm~2 mm)缺陷仿真結(jié)果

磁感應(yīng)強(qiáng)度Bz分量峰值與缺陷寬度的關(guān)系如圖7所示。從圖7中可以看出,磁感應(yīng)強(qiáng)度Bz幅值與缺陷寬度成正相關(guān),但缺陷寬度變化對(duì)磁感應(yīng)強(qiáng)度Bz幅值的影響較小,當(dāng)缺陷寬度由0.2 mm增大10倍至2 mm時(shí),磁感應(yīng)強(qiáng)度Bz幅值由0.012 mT增加至0.022 mT,增大約80%。

圖7 磁感應(yīng)強(qiáng)度Bz幅值和缺陷寬度的關(guān)系

3 實(shí)驗(yàn)驗(yàn)證

3.1 渦流檢測(cè)實(shí)驗(yàn)裝置

在仿真研究的基礎(chǔ)上,利用實(shí)驗(yàn)驗(yàn)證所設(shè)計(jì)陣列探頭的實(shí)際檢測(cè)效果,整體實(shí)驗(yàn)流程如圖8(a)所示。渦流檢測(cè)實(shí)驗(yàn)系統(tǒng)主要包括信號(hào)發(fā)生器、功率放大器、線性TMR陣列探頭、被測(cè)試件、示波器、NI/DAQ采集卡和LABVIEW虛擬儀器數(shù)據(jù)采集模塊,如圖8(b)所示。其中,信號(hào)發(fā)生器用于產(chǎn)生頻率幅值可調(diào)的正弦信號(hào),正弦信號(hào)經(jīng)過功率放大器放大后施加于線性TMR陣列探頭的激勵(lì)線圈上,從而使待測(cè)試件中產(chǎn)生渦流。TMR傳感器陣列拾取磁場(chǎng)信號(hào)變化后進(jìn)入調(diào)理電路,采集卡采集信號(hào)進(jìn)入計(jì)算機(jī),由LABVIEW軟件搭建的采集模塊進(jìn)行數(shù)據(jù)顯示、采集和保存工作。示波器用于實(shí)時(shí)監(jiān)測(cè)輸入/輸出信號(hào)的波形變化,確保系統(tǒng)工作正常。

圖8 渦流檢測(cè)實(shí)驗(yàn)系統(tǒng)

實(shí)驗(yàn)采用TC4板試件進(jìn)行檢測(cè),試件尺寸為200 mm×400 mm×5 mm。在試件同一水平線方向上加工共6個(gè)槽型缺陷,從反面對(duì)缺陷進(jìn)行檢測(cè),探頭掃查方向如圖9所示,從左往右依次為1~6號(hào)缺陷,缺陷尺寸見表2。其中,缺陷1、5、6的寬度和深度均相同,長(zhǎng)度不同;缺陷1和缺陷3的長(zhǎng)度和深度均相同,寬度不同;缺陷1、3和缺陷2、4長(zhǎng)度和寬度各自對(duì)應(yīng)相同,其缺陷距表面深度和加工深度不同。

圖9 探頭掃查方向示意

表2 缺陷尺寸參數(shù)

3.2 實(shí)驗(yàn)結(jié)果及分析

實(shí)驗(yàn)采用激勵(lì)頻率為10 kHz的正弦激勵(lì)信號(hào)。探頭沿缺陷長(zhǎng)度方向依次掃查1~6號(hào)缺陷,采集得到電壓信號(hào)輸出結(jié)果如圖10所示,其中橫軸代表電壓信號(hào)隨探頭掃查過程產(chǎn)生的變化,由于掃查速度由人為控制,存在誤差,因此以時(shí)間單位表示。

圖10 缺陷掃描檢測(cè)結(jié)果

缺陷引起的信號(hào)幅值變化數(shù)據(jù)及其反算得到的試件表面磁感應(yīng)強(qiáng)度Bz分量如表3所示。將實(shí)驗(yàn)中5、1、6號(hào)缺陷反算得到的表面磁感應(yīng)強(qiáng)度Bz值與仿真結(jié)果中不同長(zhǎng)度缺陷的磁感應(yīng)強(qiáng)度Bz值進(jìn)行對(duì)比,如圖11所示。

表3 信號(hào)幅值變化數(shù)據(jù)

圖11 不同長(zhǎng)度缺陷實(shí)驗(yàn)結(jié)果與仿真結(jié)果對(duì)比

結(jié)合表3數(shù)據(jù),對(duì)比5、1、6號(hào)缺陷實(shí)驗(yàn)結(jié)果可知,輸出電壓信號(hào)幅值變化受缺陷長(zhǎng)度變化影響較大,缺陷長(zhǎng)度由3 mm增加至6 mm時(shí),缺陷引起的電壓信號(hào)幅值變化增大約3倍,當(dāng)缺陷長(zhǎng)度由6 mm增加至12 mm時(shí),缺陷引起的信號(hào)幅值變化增大約30%。由圖11可知,反算得到的磁感應(yīng)強(qiáng)度Bz分量變化趨勢(shì)與仿真結(jié)果基本一致。對(duì)比1、3號(hào)缺陷實(shí)驗(yàn)結(jié)果可知,電壓信號(hào)幅值變化受缺陷寬度變化影響較小,缺陷寬度由0.2 mm增加至0.4 mm時(shí),幅值變化增大約8%。對(duì)比1、2號(hào)缺陷實(shí)驗(yàn)結(jié)果可知,電壓信號(hào)幅值變化受缺陷深度變化影響較大,實(shí)驗(yàn)表明該探頭能有效檢出TC4材料深層4 mm下微小缺陷。實(shí)驗(yàn)對(duì)不同長(zhǎng)度、寬度的深層裂紋缺陷的檢測(cè)結(jié)果與上節(jié)得到的仿真分析結(jié)果具有較好的一致性。

4 結(jié)論

本文針對(duì)鈦合金板件深層裂紋缺陷檢測(cè),設(shè)計(jì)了一種基于TMR傳感器陣列的雙矩形線圈差分激勵(lì)式平面渦流探頭,通過有限元仿真,模擬了探頭的移動(dòng)掃查過程,得到含不同尺寸深層缺陷的鈦合金板件表面的磁場(chǎng)圖像,結(jié)合圖像分析預(yù)測(cè)了該探頭對(duì)不同尺寸深層缺陷的檢測(cè)能力。并制作實(shí)驗(yàn)探頭,通過實(shí)驗(yàn)驗(yàn)證了探頭對(duì)鈦合金深層微小裂紋缺陷渦流檢測(cè)的可行性。主要得出了以下結(jié)論:

①探頭對(duì)缺陷長(zhǎng)度尺寸變化的靈敏度有如下規(guī)律:鈦合金板件表面3 mm下長(zhǎng)12 mm、6 mm、3 mm的微裂紋引起的感應(yīng)磁場(chǎng)強(qiáng)度分別為0.01992 mT、0.0152 mT、0.00528 mT;表面4 mm下長(zhǎng)6 mm的微裂紋引起的感應(yīng)磁場(chǎng)強(qiáng)度大約為0.00448 mT。所設(shè)計(jì)的TMR陣列渦流探頭對(duì)缺陷長(zhǎng)度變化敏感。通過傳感器信號(hào)幅值反算得到的試件表面磁感應(yīng)強(qiáng)度數(shù)值與仿真結(jié)果能保持較好的一致性。

②缺陷寬度變化與信號(hào)幅值的變化基本呈線性關(guān)系,當(dāng)缺陷寬度由0.2 mm增加至0.4 mm時(shí),輸出信號(hào)幅值增大約8%,寬度變化對(duì)信號(hào)幅值的影響不顯著。③該探頭能夠檢出鈦合金板件表面3 mm和4 mm下的深層微小裂紋缺陷,初步驗(yàn)證了探頭對(duì)非鐵磁性材料深層裂紋缺陷的檢測(cè)性能。

猜你喜歡
檢測(cè)
QC 檢測(cè)
“不等式”檢測(cè)題
“一元一次不等式”檢測(cè)題
“一元一次不等式組”檢測(cè)題
“幾何圖形”檢測(cè)題
“角”檢測(cè)題
“有理數(shù)的乘除法”檢測(cè)題
“有理數(shù)”檢測(cè)題
“角”檢測(cè)題
“幾何圖形”檢測(cè)題
主站蜘蛛池模板: 亚洲a级毛片| 欧美国产菊爆免费观看 | 四虎免费视频网站| 国产在线一区视频| 久久精品电影| 国产精品成人一区二区不卡| 久久久久无码精品| 麻豆精品国产自产在线| 亚洲精品在线影院| 99视频精品在线观看| 中文成人无码国产亚洲| www.亚洲国产| 国产手机在线小视频免费观看| 丁香五月婷婷激情基地| 国产午夜福利在线小视频| 亚洲欧美日韩中文字幕在线| 九色视频最新网址 | 污视频日本| 国产尤物在线播放| 国产欧美日韩资源在线观看| 午夜无码一区二区三区| 国产精品美乳| 三上悠亚精品二区在线观看| 99精品伊人久久久大香线蕉| 91九色最新地址| 99久久99这里只有免费的精品| 亚洲精品无码高潮喷水A| 自拍亚洲欧美精品| 国产精品va免费视频| 国产成人高清在线精品| 欧美色视频网站| 国产美女无遮挡免费视频| 欧美成人日韩| 九月婷婷亚洲综合在线| 久久精品女人天堂aaa| 女人爽到高潮免费视频大全| 日韩在线永久免费播放| 精品一区二区三区无码视频无码| 国产在线视频自拍| 国产精品久久久久久影院| 精品福利国产| 国产成人精品一区二区| 亚洲高清中文字幕| 日本一区二区不卡视频| 日本道综合一本久久久88| 亚洲成人一区二区三区| 九色综合伊人久久富二代| 毛片在线看网站| 青草精品视频| 天堂亚洲网| 制服丝袜在线视频香蕉| 亚洲国产成人在线| 四虎AV麻豆| 无码粉嫩虎白一线天在线观看| 91福利国产成人精品导航| 亚洲人成高清| 99久久人妻精品免费二区| 青青操国产| 国产欧美视频一区二区三区| 欧美天天干| 亚洲无线观看| 欧美一区中文字幕| 色香蕉影院| 国产成人精品视频一区视频二区| 久久久久久久蜜桃| 免费国产小视频在线观看| 久久成人免费| 亚洲美女一区| 九色国产在线| 欧美日韩一区二区三区四区在线观看| 亚洲人成影院午夜网站| 在线精品亚洲一区二区古装| 波多野结衣视频网站| 中文字幕伦视频| 久久久久青草大香线综合精品 | 国产无码性爱一区二区三区| 亚洲精品综合一二三区在线| 成人欧美日韩| 久久综合亚洲色一区二区三区| 欧美天堂久久| 亚洲第一区精品日韩在线播放| 成人永久免费A∨一级在线播放|