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

基于聲信號(hào)識(shí)別的水下焊接質(zhì)量檢測(cè)方法研究

2023-11-10 06:32:44紀(jì)曉東程天宇華亮張新松
關(guān)鍵詞:信號(hào)實(shí)驗(yàn)檢測(cè)

紀(jì)曉東, 程天宇, 華亮, 張新松

(南通大學(xué) 電氣工程學(xué)院, 江蘇 南通 226000)

隨著科學(xué)技術(shù)的不斷發(fā)展,水下焊接技術(shù)在現(xiàn)代造船業(yè)、航空制造業(yè)和汽車制造業(yè)等領(lǐng)域均起著關(guān)鍵性作用[1]。對(duì)于應(yīng)用場(chǎng)景較多的船舶檢修而言,在發(fā)生結(jié)構(gòu)性故障時(shí),水下焊接技術(shù)無疑是搶修的重要手段,確保焊接焊縫質(zhì)量顯得至關(guān)重要[2]。焊接質(zhì)量檢測(cè)一般分為焊后檢測(cè)和焊中檢測(cè),其中焊中檢測(cè)是目前的研究熱點(diǎn)。判定焊接質(zhì)量的主要途徑包括光信號(hào)檢測(cè)法、X射線檢測(cè)法、磁粉檢測(cè)法、超聲無損檢測(cè)法、機(jī)器視覺檢測(cè)法和聲信號(hào)檢測(cè)法等。如Creaform公司[3]提出了一種海上風(fēng)電導(dǎo)管架無損檢測(cè)法,可實(shí)現(xiàn)導(dǎo)管架的精準(zhǔn)高效檢測(cè),但該方法僅針對(duì)焊后質(zhì)量檢測(cè),缺乏實(shí)時(shí)性;肖書浩等人[4]提出了基于光電信號(hào)的陶質(zhì)襯墊焊接質(zhì)量在線檢測(cè)方法,解決了焊接質(zhì)量檢測(cè)的實(shí)時(shí)性問題,但該方法僅適用于陸地檢測(cè)。在實(shí)際工程應(yīng)用中,一些焊工依靠多年的焊接經(jīng)驗(yàn),僅通過聲音即可大致判斷焊接質(zhì)量[5]。由此可見,聲信號(hào)檢測(cè)法可以直接通過焊接過程中產(chǎn)生的聲信息進(jìn)行焊接質(zhì)量判定,實(shí)時(shí)性強(qiáng)且效率高,被認(rèn)為是最具潛力的水下焊接質(zhì)量檢測(cè)技術(shù)之一[6]。

如今,通過構(gòu)建基于聲信號(hào)采集的監(jiān)測(cè)系統(tǒng)來實(shí)時(shí)檢測(cè)焊接質(zhì)量的技術(shù)已逐步發(fā)展起來。如劉亮[7]對(duì)鎢極氣體電弧焊過程中產(chǎn)生的聲信號(hào)進(jìn)行采集并開展時(shí)域特性分析,并將其應(yīng)用于熔透焊接質(zhì)量的實(shí)時(shí)控制。但是,在實(shí)際焊接過程中,因燃燒電弧而導(dǎo)致的體積膨脹、過度熔滴對(duì)熔池沖刷會(huì)引起振蕩,焊接時(shí)無疑會(huì)夾雜著巨大的聲音,同時(shí)還存在設(shè)備運(yùn)作聲、熔滴飛濺等雜音以及焊接缺陷、焊接過程不穩(wěn)定等狀況,使得焊接聲信號(hào)極易受環(huán)境干擾[8]。隨著神經(jīng)網(wǎng)絡(luò)技術(shù)的發(fā)展,其在質(zhì)量檢測(cè)中的應(yīng)用日漸成熟。為解決類似水下焊接質(zhì)量檢測(cè)等較為復(fù)雜的問題,許多學(xué)者開始研究基于人工神經(jīng)網(wǎng)絡(luò)的質(zhì)量檢測(cè)技術(shù)[9]。其中,雙權(quán)值神經(jīng)網(wǎng)絡(luò)(double-weight neural network, DWNN)模型的非線性逼近能力強(qiáng),具有學(xué)習(xí)速度快、網(wǎng)絡(luò)結(jié)構(gòu)好等特性[10],適用于處理小樣本非線性數(shù)據(jù)的擬合問題,如高維數(shù)據(jù)非線性擬合等,可針對(duì)性地解決水下焊接質(zhì)量檢測(cè)不確定因素多等問題。

為解決水下焊接過程中存在的非線性程度高、參數(shù)耦合性強(qiáng)以及檢測(cè)效率低等問題,筆者以航空母艦用高強(qiáng)度低碳合金鋼——HSLA-115鋼為研究對(duì)象,提出了基于聲信號(hào)識(shí)別的水下焊接質(zhì)量檢測(cè)方法。首先,通過搭建水下焊接實(shí)驗(yàn)平臺(tái),獲取焊接過程中的聲信號(hào),并利用短時(shí)加窗技術(shù)對(duì)聲信號(hào)進(jìn)行特征提取;然后,將所提取的聲信號(hào)特征值作為訓(xùn)練樣本,構(gòu)建基于DWNN的水下焊接質(zhì)量識(shí)別模型,旨在為我國(guó)水下焊接工藝的發(fā)展提供新思路。

1 水下焊接質(zhì)量判定

所搭建的水下焊接實(shí)驗(yàn)平臺(tái)如圖1所示,包含焊接母材、焊條和焊接電源等。其中:焊接母材為航空母艦用高強(qiáng)度低碳合金鋼——HSLA-115鋼,其具有強(qiáng)度高、韌性好和耐腐蝕等特性;焊條采用BROCO藥皮包覆,由石蠟防水層、焊芯和藥皮層組成,其具有可在水下持續(xù)燃燒的性能;焊接電源采用KD-630M型三相逆變式直流焊接電源,其采用IGBT(insulated gate bipolar transistor,絕緣柵雙極型晶體管)功率模塊配合逆變控制技術(shù),具有焊接電流穩(wěn)定的特點(diǎn)。

圖1 水下焊接實(shí)驗(yàn)平臺(tái)Fig.1 Underwater welding experiment platform

在焊接實(shí)驗(yàn)開始前,先將HSLA-115鋼加工成200 mm×100 mm的標(biāo)準(zhǔn)化焊件,再將部分焊件切V字口。在焊接實(shí)驗(yàn)中,采用2種焊接速度(0.6,1.2 m/min)、4種板材厚度(6,8,10,12 mm)和5種焊接電流(110,140,150,160,180 A)分別進(jìn)行水下焊接處理,每一種實(shí)驗(yàn)組合均進(jìn)行6次,共開展240(2×4×5×6=240)組實(shí)驗(yàn)。在焊接實(shí)驗(yàn)結(jié)束后,通過觀察焊縫表面的光滑程度及其熔透狀態(tài)來判定焊接質(zhì)量:若焊縫表層光滑,并與焊件邊緣深度熔合且連接順暢,則判定質(zhì)量合格;若焊縫出現(xiàn)毛糙、多氣孔、未熔合等現(xiàn)象,則判定質(zhì)量不合格。部分水下焊接實(shí)驗(yàn)的測(cè)定結(jié)果如表1所示,焊接鋼板實(shí)物如圖2所示。

表1 部分水下焊接實(shí)驗(yàn)的測(cè)定結(jié)果Table 1 Measurement results of partial underwater welding experiments

圖2 水下焊接鋼板實(shí)物Fig.2 Physical objects of underwater welded steel plate

2 聲信號(hào)的采集與處理

圖1所示實(shí)驗(yàn)平臺(tái)采用水聽器獲取水下焊接過程中的聲信號(hào)。水聽器由中科探海公司提供,其型號(hào)為TH-AC0200,監(jiān)聽頻率范圍為0.5~30 kHz,采樣頻率為2 MHz,帶內(nèi)波動(dòng)小于3 dB,動(dòng)態(tài)范圍時(shí)變?cè)鲆鏋?8 dB,固定增益為18.2 dB。該水聽器具有輸入阻抗高、輸出阻抗低和噪聲低等特點(diǎn)。采用聲影軟件對(duì)水聽器獲取的聲信號(hào)進(jìn)行采集。

2.1 聲信號(hào)預(yù)處理

采用Butterworth(巴特沃斯)濾波器對(duì)所采集的聲信號(hào)進(jìn)行濾波降噪處理。該濾波器的幅值特性曲線如圖3所示,其在衰減斜率、線性相位與加載特性三個(gè)方面均具有特性均衡的優(yōu)點(diǎn)[11],能夠有效減少聲信號(hào)中機(jī)械運(yùn)行噪聲與水環(huán)境中氣泡噪聲的干擾。Butterworth濾波器的幅值平方函數(shù)可表示為:

圖3 Butterworth濾波器的幅值特性曲線Fig.3 Amplitude characteristic curve of Butterworth filter

式中:H(ω)為濾波器的幅值;N為濾波器的階數(shù),本文取N=8;?為輸入信號(hào)的頻率;?c為低通濾波器的截止頻率。

圖4所示為在板材厚度為8 mm、焊接電流為140 A、焊接速度為1.2 m/min條件下所采集的聲信號(hào)濾波處理前后的頻譜圖對(duì)比。由圖可以看出,Butterworth濾波器有效抑制了聲信號(hào)中的原始噪聲信號(hào)。

圖4 濾波處理前后聲信號(hào)的頻譜圖對(duì)比Fig.4 Comparison of spectrogram of acoustic signal before and after filtering

2.2 聲信號(hào)特征提取

參照文獻(xiàn)[12]中所提出的特征提取方法,采用短時(shí)加窗技術(shù)對(duì)濾波處理后的聲信號(hào)進(jìn)行特征提取。本文采用的窗函數(shù)為漢明窗,窗口長(zhǎng)度為1 024,重疊50%進(jìn)行分幀。從連續(xù)的聲信號(hào)中截取多個(gè)“分析幀”,并分別在每一個(gè)分析幀的時(shí)域特征中提取短時(shí)能量En、短時(shí)平均幅度Mn、短時(shí)平均過零率Zn這3個(gè)特征值,結(jié)果分別如圖5至圖7所示。

圖5 聲信號(hào)的短時(shí)能量Fig.5 Short-time energy of acoustic signal

圖7 聲信號(hào)的短時(shí)平均過零率Fig.7 Short-time average zero crossing rate of acoustic signal

3 基于DWNN的數(shù)據(jù)擬合算法

針對(duì)水下焊接聲信號(hào)的時(shí)域特征數(shù)據(jù),本文采用文獻(xiàn)[13]中提出的多維函數(shù)擬合逼近算法(即DWNN)進(jìn)行焊接質(zhì)量檢測(cè)。DWNN是一種新型的神經(jīng)網(wǎng)絡(luò),具有分類能力強(qiáng)、學(xué)習(xí)速度快等優(yōu)點(diǎn),其固定結(jié)構(gòu)如圖8所示。

圖8 結(jié)構(gòu)固定的DWNNFig.8 DWNN with fixed structure

在DWNN的數(shù)學(xué)模型中,對(duì)應(yīng)隱含層神經(jīng)元的任意一個(gè)輸出Yj為:

式中:Aj為隱含層第j個(gè)神經(jīng)元的激勵(lì)函數(shù);n為輸入向量的維數(shù);Wij為第i個(gè)輸入節(jié)點(diǎn)到隱含層第j個(gè)神經(jīng)元的方向權(quán)值;xi為第i個(gè)輸入節(jié)點(diǎn)的輸入值;ηij為第i個(gè)輸入節(jié)點(diǎn)到隱含層第j個(gè)神經(jīng)元的核心權(quán)值。

為避免DWNN模型在優(yōu)化過程中陷入局部最優(yōu)解,引入粒子群優(yōu)化(particle swarm optimization,PSO)算法[14]。PSO算法通過實(shí)時(shí)更新粒子的速度和位置來修正和獲取DWNN的連接權(quán)值和幅值,以幫助DWNN跳出局部最優(yōu)解,從而得到全局最優(yōu)解。在整個(gè)求解過程中,粒子自身會(huì)通過不斷的迭代更新求得一個(gè)最優(yōu)解,此最優(yōu)解即粒子自身的個(gè)體極值pbest;除此解之外,還將從整個(gè)粒子群中找到一個(gè)最優(yōu)解,即全局極值gbest[15]。根據(jù)確定好的2個(gè)最優(yōu)解,粒子通過式(3)和式(4)進(jìn)行速度和位置的實(shí)時(shí)更新:

式中:k為迭代次數(shù);c1、c2為學(xué)習(xí)因子,用于調(diào)節(jié)2個(gè)最優(yōu)解的步長(zhǎng);r1、r2為0~1之間的隨機(jī)數(shù);pbestid,k、gbestid,k分別為當(dāng)前的個(gè)體極值和整個(gè)群落中的全局極值;vid,k、xid,k分別為第i個(gè)粒子經(jīng)過k次迭代后在d維分量上的速度及位置;w為慣性權(quán)重,用于改善PSO算法的收斂性。

將PSO算法加入DWNN后,最終的數(shù)據(jù)擬合公式如下:

式中:xi=(xi1,xi2, …,xis)、vi=(vi1,vi2, …,vis)分別為第i個(gè)粒子在s維向量搜索空間中的位置和速度。

根據(jù)上文,在DWNN模型中引入PSO算法的具體步驟如下。

步驟1:初建一個(gè)基礎(chǔ)的DWNN,設(shè)定初始隱含層神經(jīng)元j=1,最大隱含層神經(jīng)元個(gè)數(shù)為L(zhǎng)max,閾值精度為ε,參照隱含層神經(jīng)元j選取輸入樣本和對(duì)應(yīng)的期望輸出,并將其進(jìn)行歸一化處理[16]。

步驟2:引入PSO算法,首先確定種群大小n,粒子的維數(shù)s,搜索范圍qi∈[qi,min,qi,max],速度vi∈[vi,min,vi,max],參數(shù)c1、c2,初始迭代次數(shù)k=1,最大迭代次數(shù)kmax等。隨后,根據(jù)設(shè)定的參數(shù),產(chǎn)生n個(gè)粒子,具體公式如下:

式中:i∈{1, 2,…,n},d∈{1, 2, …,s};r(x)為0~1之間的隨機(jī)數(shù)。

步驟3:計(jì)算粒子的適應(yīng)度值F,其由粒子的優(yōu)化目標(biāo)決定,用于評(píng)價(jià)粒子的搜索性能,指導(dǎo)粒子群的搜索過程。F可表示為DWNN的實(shí)際輸出與樣本期望輸出的均方誤差:

式中:N2為訓(xùn)練樣本數(shù);Yh為第h個(gè)樣本的網(wǎng)絡(luò)實(shí)際輸出值;Dh,j為第h個(gè)樣本對(duì)應(yīng)的第j個(gè)神經(jīng)元的期望輸出值。

步驟4:設(shè)定粒子適應(yīng)度值為粒子自身最優(yōu)解,即個(gè)體極值pbest,將其所在位置默認(rèn)為最優(yōu)位置X1,并找到粒子群中的全局極值gbest,其對(duì)應(yīng)位置默認(rèn)為最優(yōu)位置X2。

步驟5:首先根據(jù)步驟4所得的2個(gè)最優(yōu)解,更新粒子的位置和速度[17],計(jì)算更新后粒子的適應(yīng)度值F,將其與個(gè)體極值pbest對(duì)比擇優(yōu),選定最大值并保留其位置信息。同理,再將適應(yīng)度值F與全局極值gbest比較,選定最大值并保留其位置信息。

步驟6:確定最終的最優(yōu)適應(yīng)度值及其對(duì)應(yīng)的粒子位置和速度,將此作為第j個(gè)神經(jīng)元的方向權(quán)值、核心權(quán)值和幅值[18],并將最終的最優(yōu)適應(yīng)度值作為DWNN的計(jì)算精度[19]。

步驟7:判別DWNN是否滿足閾值精度ε或達(dá)到最大隱含層神經(jīng)元個(gè)數(shù),若其中一個(gè)條件滿足,則結(jié)束算法,否則繼續(xù)迭代,使j=j+1,跳至步驟2。

引入PSO的DWNN算法的總體流程如圖9所示。

圖9 引入PSO的DWNN算法流程Fig.9 Flow of DWNN algorithm introducing PSO

4 實(shí)驗(yàn)驗(yàn)證及分析

利用水聽器采集水下焊接實(shí)驗(yàn)過程中產(chǎn)生的聲信號(hào),隨后對(duì)其進(jìn)行濾波降噪處理和特征提取。以特征值En、Mn、Zn作為輸入,以焊接質(zhì)量作為輸出,基于特征樣本搭建圖10所示的DWNN,完成基于聲信號(hào)識(shí)別的水下焊接質(zhì)量檢測(cè)。

圖10 水下焊接質(zhì)量檢測(cè)用DWNN的結(jié)構(gòu)Fig.10 Structure of DWNN for underwater welding quality detection

從240組測(cè)試樣本中隨機(jī)選取50組焊接質(zhì)量合格的樣本及50組焊接質(zhì)量不合格的樣本形成訓(xùn)練樣本集,用于構(gòu)建DWNN模型。同時(shí),選取支持向量機(jī)(support vector machine, SVM)及反向傳播神經(jīng)網(wǎng)絡(luò)(back propagation neural network, BPNN)作為對(duì)照。SVM及BPNN的核心思想也是將特征映射至高維空間中并通過線性回歸完成分類任務(wù)[20]。然后,分別從剩余樣本中隨機(jī)選取40組作為驗(yàn)證樣本集,其中20組焊接質(zhì)量合格的驗(yàn)證結(jié)果用“1”表示,20組焊接質(zhì)量不合格的用“0”表示,結(jié)果如圖11所示。

圖11 基于不同方法的水下焊接質(zhì)量識(shí)別結(jié)果(100組訓(xùn)練樣本)Fig.11 Underwater welding quality recognition results based on different methods (100 groups of training samples)

在實(shí)驗(yàn)中,設(shè)定驗(yàn)證值與真實(shí)值的偏差大于0.5時(shí)為識(shí)別錯(cuò)誤。由圖11可知,對(duì)于40組測(cè)試樣本,DWNN模型均可實(shí)現(xiàn)有效檢測(cè),識(shí)別率達(dá)到了100%;而SVM模型出現(xiàn)了3次錯(cuò)誤識(shí)別,BPNN模型出現(xiàn)了2次錯(cuò)誤識(shí)別。為了進(jìn)一步描述不同方法的識(shí)別偏差,用誤差(真實(shí)值與測(cè)試值之間的差值)來量化識(shí)別模型的檢測(cè)精度,結(jié)果如圖12所示,誤差分析結(jié)果如表2所示。

表2 不同水下焊接質(zhì)量識(shí)別模型的誤差對(duì)比Table 2 Error comparison of different underwater welding quality recognition models

圖12 基于不同方法的水下焊接質(zhì)量識(shí)別誤差(100組訓(xùn)練樣本)Fig.12 Underwater welding quality recognition error based on different methods (100 groups of training samples)

由圖12可知,DWNN模型識(shí)別的平均誤差約為0.002 2,最大誤差約為0.011 4,均方誤差僅為0.000 387,顯著優(yōu)于SVM模型和BPNN模型,表明基于DWNN的水下焊接質(zhì)量檢測(cè)方法的識(shí)別精度更高,抗干擾能力更強(qiáng)。具體分析誤差產(chǎn)生原因可知,由于焊縫不平整與噪聲等因素干擾,在對(duì)焊接質(zhì)量合格樣本進(jìn)行質(zhì)量識(shí)別時(shí)會(huì)存在一定誤差。而對(duì)于焊接質(zhì)量不合格的樣本,可實(shí)現(xiàn)精準(zhǔn)、有效識(shí)別,即識(shí)別誤差均為0,說明所構(gòu)建的DWNN模型對(duì)不合格焊件的識(shí)別精度更高。

一般焊件材料如304不銹鋼、T10工具鋼等的價(jià)格低廉,用途廣泛,測(cè)試樣本大小不受限制。但是,本文所研究的HSLA-115鋼的價(jià)格昂貴,多用于航空母艦,這類金屬焊件在一般情況下不具備獲取大量訓(xùn)練樣本的背景,不適合使用監(jiān)督式深度學(xué)習(xí)神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練和識(shí)別。因此,需要在小樣本情況下,對(duì)DWNN模型的有效性和識(shí)別率進(jìn)行測(cè)試。在240組測(cè)試樣本中重新抽取10組焊接質(zhì)量合格的樣本及10組焊接質(zhì)量不合格的樣本形成訓(xùn)練樣本集,并重新構(gòu)建DWNN模型。同時(shí),從剩余樣本中隨機(jī)選取40組作為驗(yàn)證樣本集,同樣用“1”表示焊接質(zhì)量合格,用“0”表示焊接質(zhì)量不合格,識(shí)別結(jié)果和誤差分析結(jié)果分別如圖13和圖14所示。

圖13 基于不同方法的水下焊接質(zhì)量識(shí)別結(jié)果(20組訓(xùn)練樣本)Fig.13 Underwater welding quality recognition results based on different methods (20 groups of training samples)

圖14 基于不同方法的水下焊接質(zhì)量識(shí)別誤差(20組訓(xùn)練樣本)Fig.14 Underwater welding quality recognition error based on different methods (100 groups of training samples)

由圖13可知,僅用20組訓(xùn)練樣本構(gòu)建識(shí)別模型時(shí),SVM模型出現(xiàn)了5次錯(cuò)誤識(shí)別,BPNN模型出現(xiàn)了3次錯(cuò)誤識(shí)別。此外,通過對(duì)比圖12和圖14可知,SVM模型與BPNN模型的識(shí)別誤差顯著增大,而DWNN模型的識(shí)別誤差無明顯變化。

綜合2次實(shí)驗(yàn)的驗(yàn)證結(jié)果可得,當(dāng)訓(xùn)練樣本總數(shù)為100組時(shí),DWNN模型相對(duì)于SVM模型和BPNN模型的識(shí)別精度更高,且識(shí)別過程更穩(wěn)定;當(dāng)訓(xùn)練樣本總數(shù)為20組時(shí),由于訓(xùn)練樣本數(shù)量的減少,SVM模型和BPNN模型的識(shí)別率明顯降低,而DWNN模型依然保持高識(shí)別率,說明DWNN模型不用依靠大量訓(xùn)練樣本就能通過高維空間的特征映射完成特征向量的最優(yōu)覆蓋。

5 結(jié) 論

為解決水下焊接質(zhì)量檢測(cè)中存在的實(shí)時(shí)性差、識(shí)別精度低等問題,基于不同板材厚度、焊接電流和焊接速度條件下水下焊接過程中產(chǎn)生的聲信號(hào)的特征差異,通過提取短時(shí)能量、短時(shí)平均幅度和短時(shí)平均過零率這3個(gè)時(shí)域特征值作為訓(xùn)練樣本,以構(gòu)建DWNN模型,實(shí)現(xiàn)了水下焊接質(zhì)量的實(shí)時(shí)檢測(cè)。實(shí)驗(yàn)結(jié)果表明,不同焊接參數(shù)下水下焊接過程中的聲信號(hào)特征與焊接質(zhì)量之間存在映射關(guān)系,DWNN模型能夠在不依賴大樣本條件下精確表征映射關(guān)系,為基于聲信號(hào)識(shí)別的水下焊接質(zhì)量檢測(cè)提供了有效方法,同時(shí)也為HSLA-115鋼的水下焊接參數(shù)優(yōu)化以及水下焊件專家知識(shí)庫(kù)的構(gòu)建提供了參考依據(jù)。

猜你喜歡
信號(hào)實(shí)驗(yàn)檢測(cè)
記一次有趣的實(shí)驗(yàn)
“不等式”檢測(cè)題
“一元一次不等式”檢測(cè)題
“一元一次不等式組”檢測(cè)題
信號(hào)
鴨綠江(2021年35期)2021-04-19 12:24:18
完形填空二則
做個(gè)怪怪長(zhǎng)實(shí)驗(yàn)
基于FPGA的多功能信號(hào)發(fā)生器的設(shè)計(jì)
電子制作(2018年11期)2018-08-04 03:25:42
NO與NO2相互轉(zhuǎn)化實(shí)驗(yàn)的改進(jìn)
實(shí)踐十號(hào)上的19項(xiàng)實(shí)驗(yàn)
太空探索(2016年5期)2016-07-12 15:17:55
主站蜘蛛池模板: 毛片手机在线看| 欧美精品在线看| 久久精品国产91久久综合麻豆自制| 亚洲成aⅴ人在线观看| 中文字幕 91| 成人看片欧美一区二区| 免费无码又爽又黄又刺激网站| 午夜啪啪网| 亚洲色婷婷一区二区| 国产成人无码Av在线播放无广告| 国产二级毛片| 2019年国产精品自拍不卡| 国内熟女少妇一线天| 日韩毛片免费观看| 国产v欧美v日韩v综合精品| 欧美精品一区二区三区中文字幕| 日本黄网在线观看| 亚洲专区一区二区在线观看| 亚洲自拍另类| 黄色网站不卡无码| 91精品综合| 欧美日韩一区二区在线免费观看| 精品国产一区二区三区在线观看 | 国产无遮挡猛进猛出免费软件| 亚洲第一色网站| 动漫精品中文字幕无码| 色网在线视频| 亚洲人网站| 亚洲精品在线影院| 国产三区二区| 中文字幕在线一区二区在线| Aⅴ无码专区在线观看| 国产网站一区二区三区| 亚洲精品图区| 久久综合五月婷婷| 欧美性色综合网| 国产成在线观看免费视频| 国产av色站网站| 综合社区亚洲熟妇p| 久久99精品国产麻豆宅宅| 免费一级毛片不卡在线播放| 精品99在线观看| 色婷婷天天综合在线| 噜噜噜综合亚洲| 亚洲国产成人自拍| 国产人成乱码视频免费观看| 国产精品午夜福利麻豆| 亚洲乱亚洲乱妇24p| swag国产精品| 精品在线免费播放| 日韩一级二级三级| 免费黄色国产视频| 国产一区二区色淫影院| 中文字幕资源站| 日韩精品成人在线| 国内精品免费| 久久精品无码中文字幕| 在线99视频| 国产一区二区丝袜高跟鞋| 免费jizz在线播放| 国产精品成人一区二区不卡 | 国产资源站| 午夜日韩久久影院| 亚洲人免费视频| 午夜国产大片免费观看| 久草视频一区| 日韩一区二区在线电影| 亚洲精品动漫| 亚洲人成影视在线观看| 国产女人18毛片水真多1| 久久青草免费91观看| 成年人免费国产视频| 成人在线第一页| 国产无码性爱一区二区三区| 午夜在线不卡| 亚洲一级无毛片无码在线免费视频| 久久免费观看视频| 成人永久免费A∨一级在线播放| 秘书高跟黑色丝袜国产91在线| 91系列在线观看| аⅴ资源中文在线天堂| 国产97公开成人免费视频|