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

模糊神經(jīng)網(wǎng)絡(luò)變幅液壓系統(tǒng)故障診斷

2014-04-03 07:32:26馮文潔李萬莉嘉紅霞
關(guān)鍵詞:故障診斷故障系統(tǒng)

馮文潔,李萬莉,嘉紅霞

FENG Wenjie1,LI Wanli2,JIA Hongxia3

1.同濟(jì)大學(xué) 機(jī)械與能源工程學(xué)院,上海 201804

2.上海海事大學(xué) 物流工程學(xué)院,上海 201306

1.School of Mechanical Engineering,Tongji University,Shanghai 201804,China

2.School of Logistics Engineering,Shanghai Maritime University,Shanghai 201306,China

1 引言

變幅液壓系統(tǒng)結(jié)構(gòu)復(fù)雜,由大量的動(dòng)力元件、控制元件和執(zhí)行元件有機(jī)組合而成,各組成元件在一個(gè)工作循環(huán)中通過不同的動(dòng)作完成特定的任務(wù),因此系統(tǒng)內(nèi)部各元件之間耦合性強(qiáng),系統(tǒng)故障及其故障原因具有多樣性、模糊性及偶然性的特點(diǎn)。

模糊神經(jīng)網(wǎng)絡(luò)利用神經(jīng)網(wǎng)絡(luò)的自學(xué)習(xí)功能,優(yōu)化模糊邏輯系統(tǒng)中的模糊規(guī)則、隸屬函數(shù)和模糊決策算法;將神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)結(jié)果轉(zhuǎn)化為模糊邏輯系統(tǒng)的規(guī)則知識(shí),從而更有利于知識(shí)的解釋與利用,非常適用于復(fù)雜對(duì)象的過程監(jiān)測(cè)與故障診斷[1-3],則對(duì)變幅液壓系統(tǒng)故障診斷具有非常大的優(yōu)越性。學(xué)習(xí)樣本的歸納一直是神經(jīng)網(wǎng)絡(luò)應(yīng)用的一大難點(diǎn)[4]。故障樹描述并存貯了系統(tǒng)關(guān)于故障樹頂部事件發(fā)生的全部知識(shí),通過構(gòu)建變幅液壓系統(tǒng)故障樹,由故障樹提取變幅液壓系統(tǒng)模糊神經(jīng)網(wǎng)絡(luò)的訓(xùn)練數(shù)據(jù),既克服了提取訓(xùn)練數(shù)據(jù)的盲目性,又可包含系統(tǒng)的全部故障模式,使神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)數(shù)據(jù)的歸納簡(jiǎn)單易行。

因此,本文以基于故障樹的模糊神經(jīng)網(wǎng)絡(luò)故障診斷方法對(duì)變幅液壓系統(tǒng)進(jìn)行故障診斷。

2 變幅液壓系統(tǒng)工作原理

圖1為變幅液壓系統(tǒng)工作原理圖。變幅換向閥4右位時(shí),變量泵1高壓油,經(jīng)平衡閥5的單向閥,進(jìn)入變幅缸6大腔,小腔回油,缸活塞桿伸長(zhǎng),驅(qū)動(dòng)工作臂逆時(shí)針轉(zhuǎn)動(dòng),減小幅度,提高起升高度。換向閥4左位時(shí),高壓油進(jìn)入小腔,在重力作用下下降,由于平衡閥產(chǎn)生背壓作用,臂平穩(wěn)下降,減少起升高度,加大幅度。換向閥處于中位時(shí),由于平衡閥5具有鎖的作用,使臂在作業(yè)時(shí),不會(huì)產(chǎn)生沉降,保證機(jī)械的穩(wěn)定性。

圖1 變幅液壓系統(tǒng)工作原理圖

分析可知變幅液壓系統(tǒng)結(jié)構(gòu)復(fù)雜,由動(dòng)力元件、控制元件和執(zhí)行元件有機(jī)組合而成,各組成元件在一個(gè)工作循環(huán)中通過不同的動(dòng)作完成特定的任務(wù),因此系統(tǒng)內(nèi)部各元件之間耦合性強(qiáng),系統(tǒng)故障及其故障原因具有多樣性、模糊性及偶然性的特點(diǎn)?;谝陨瞎收咸攸c(diǎn)分析,本文提出以模糊神經(jīng)網(wǎng)絡(luò)作為變幅液壓系統(tǒng)的故障診斷方法。

3 變幅液壓系統(tǒng)模糊神經(jīng)網(wǎng)絡(luò)故障診斷模型

根據(jù)故障樹知識(shí)以及模糊神經(jīng)網(wǎng)絡(luò)的特點(diǎn),首先,從故障樹知識(shí)中提取模糊神經(jīng)網(wǎng)絡(luò)輸入?yún)?shù);其次,收集一定量的樣本數(shù)據(jù),運(yùn)用統(tǒng)計(jì)參數(shù)法中的3σ準(zhǔn)則法確立模糊化處理的模糊隸屬函數(shù),并對(duì)模糊神經(jīng)網(wǎng)絡(luò)輸入?yún)?shù)進(jìn)行模糊化處理;再次,建立神經(jīng)網(wǎng)絡(luò)推理計(jì)算的權(quán)值和閾值的訓(xùn)練算法;最后,對(duì)模糊化的輸入?yún)?shù)進(jìn)行訓(xùn)練并清晰化訓(xùn)練結(jié)果。由此構(gòu)建變幅液壓系統(tǒng)模糊神經(jīng)網(wǎng)絡(luò)故障診斷模型[5-6],如圖2所示。

3.1 變幅液壓系統(tǒng)故障診斷輸入輸出節(jié)點(diǎn)的確定

在變幅液壓系統(tǒng)故障診斷模糊神經(jīng)網(wǎng)絡(luò)中,網(wǎng)絡(luò)的輸入變量表達(dá)系統(tǒng)的故障征兆或故障模式的特征參數(shù),輸出變量則表達(dá)系統(tǒng)的故障原因或故障模式。為了確定變幅液壓系統(tǒng)故障診斷網(wǎng)絡(luò)的輸入及輸出變量,以變幅無動(dòng)作為例,建立變幅無動(dòng)作的故障樹,如圖3所示。

圖2 變幅液壓系統(tǒng)故障診斷模糊神經(jīng)網(wǎng)絡(luò)

圖3 變幅無動(dòng)作故障樹

故障樹中各基本事件的表述如表1所示。

表1 變幅無動(dòng)作故障樹基本事件

根據(jù)故障樹最小割集的定義及上行法計(jì)算可知,變幅液壓系統(tǒng)故障樹中每一個(gè)基本事件就是一個(gè)最小割集,只要檢測(cè)主泵壓力、變幅大腔壓力、變幅小腔壓力及平衡閥控制油口壓力即可分辨無故障、{X1-X2}、{X3}、{X4-X5}、{X6-X7}及{X8-X10}這六種故障模式,因此這4個(gè)特征參數(shù)即為故障診斷網(wǎng)絡(luò)輸入變量,輸出變量為6個(gè)。

3.2 變幅液壓系統(tǒng)故障診斷輸入變量的模糊化處理

3.2.1 變幅液壓系統(tǒng)故障診斷輸入變量模糊化隸屬度函數(shù)

通過隸屬度函數(shù)[7]將輸入變量轉(zhuǎn)換成對(duì)應(yīng)模糊子集上的隸屬度。模糊層的作用即利用隸屬度函數(shù)將輸入變量進(jìn)行模糊量化預(yù)處理,使之符合模糊神經(jīng)網(wǎng)絡(luò)的輸入要求,模糊化采用的隸屬度函數(shù),如圖4所示[8-9]。

圖4 模糊隸屬度函數(shù)

圖4中,狀態(tài)變量域低的轉(zhuǎn)換采用降半梯形分布的隸屬函數(shù),其表達(dá)式如(1)所示。

狀態(tài)變量域正常的轉(zhuǎn)換采用梯形分布的隸屬函數(shù),其表達(dá)式如(2)所示。

狀態(tài)變量域高的轉(zhuǎn)換采用升半梯形分布的隸屬函數(shù),其表達(dá)式如(3)所示。

特征參數(shù)的模糊化處理是將輸入的精確特征參數(shù)值轉(zhuǎn)換成模糊隸屬度值。根據(jù)上述分析,首先將輸入的特征參數(shù)變換到各自所屬的論域范圍即低、正常和高,然后將已經(jīng)變換到論域范圍的特征信號(hào)參數(shù)根據(jù)各個(gè)論域的隸屬度函數(shù)進(jìn)行模糊化處理,使原來精確的輸入量變成0~1之間的模糊隸屬度值。

3.2.2 變幅液壓系統(tǒng)故障診斷隸屬度函數(shù)參數(shù)的計(jì)算

在收集一定量的設(shè)備正常工作時(shí)的樣本數(shù)據(jù)基礎(chǔ)上,運(yùn)用統(tǒng)計(jì)參數(shù)法中的3σ準(zhǔn)則法[10-12]確定隸屬函數(shù)各參數(shù)的值。設(shè)備正常工作時(shí)的各類診斷參數(shù)即模糊神經(jīng)網(wǎng)絡(luò)輸入變量xi服從正態(tài)分布規(guī)律,對(duì)正常工作時(shí)各類診斷參數(shù)xi收集一定量的訓(xùn)練樣本數(shù)據(jù),樣本個(gè)數(shù)設(shè)為N,然后求取各輸入變量的樣本均值,其中xik為模糊神經(jīng)網(wǎng)絡(luò)第i個(gè)輸入變量的第k個(gè)樣本:

設(shè)σ為標(biāo)準(zhǔn)方差,則:

若診斷參數(shù)xi服從正態(tài)分布,那么,在±3σ范圍內(nèi)包含了99.73%的質(zhì)量特性值。按正態(tài)分布理論根據(jù)3σ原理,如果測(cè)試樣本與均值距離在3倍方差以上,則認(rèn)為該測(cè)試樣本是不確定的,故可判為異常,以此可以確定控制界限[11-12]:故障界限:X±3σ;警告界限:X±σ或者 X±2σ 。

將3σ原理應(yīng)用于隸屬函數(shù)各參數(shù)的確定,a為故障下限,d為故障上限,b為警告下限,c為警告上限,可得到模糊神經(jīng)網(wǎng)絡(luò)中第i個(gè)輸入變量的隸屬函數(shù)參數(shù)計(jì)算公式為:

3.3 變幅液壓系統(tǒng)故障診斷神經(jīng)網(wǎng)絡(luò)推理

3.3.1 變幅液壓系統(tǒng)故障診斷神經(jīng)網(wǎng)絡(luò)權(quán)值和閾值的求解

第一層故障樹層可知輸入變量的個(gè)數(shù),設(shè)輸入變量為n,則節(jié)點(diǎn)數(shù)為n,其輸入輸出為:i=1,2,…,n。故障樹層直接將輸入變量傳遞給第二層模糊化層,在基于故障樹的模糊神經(jīng)網(wǎng)絡(luò)中,輸入變量為系統(tǒng)的故障征兆。因此第二層的輸出節(jié)點(diǎn)數(shù)仍然為n,其輸入輸出為:

第三層為神經(jīng)網(wǎng)絡(luò)推理層,神經(jīng)網(wǎng)絡(luò)推理層神經(jīng)元采用的傳遞函數(shù)是tansig。

設(shè)該隱含層節(jié)點(diǎn)為q個(gè),節(jié)點(diǎn)的輸入輸出為:

式中ωij為第二層第i個(gè)節(jié)點(diǎn)到第三層第 j個(gè)節(jié)點(diǎn)的連接權(quán)值,θj為第 j個(gè)節(jié)點(diǎn)的閾值,f(s)為tansig函數(shù)。

第四層為輸出層,輸出層神經(jīng)元采用的傳遞函數(shù)是logsig。

式中ωjl為第三層第 j個(gè)節(jié)點(diǎn)到第四層第l個(gè)節(jié)點(diǎn)的連接權(quán)值,θl為第l個(gè)節(jié)點(diǎn)的閾值,f(s)為logsig函數(shù)。

3.3.2 變幅液壓系統(tǒng)故障診斷神經(jīng)網(wǎng)絡(luò)權(quán)值和閾值的訓(xùn)練

運(yùn)用Levenberg-Marquardt[13]優(yōu)化算法對(duì)神經(jīng)網(wǎng)絡(luò)的連接權(quán)值和閾值進(jìn)行訓(xùn)練,網(wǎng)絡(luò)誤差指標(biāo)函數(shù)表達(dá)為[14-15]:

式中,Yq表示實(shí)際的神經(jīng)網(wǎng)絡(luò)輸出向量,Yq'表示期望的神經(jīng)網(wǎng)絡(luò)輸出向量,P為神經(jīng)網(wǎng)絡(luò)訓(xùn)練樣本數(shù)目,w為神經(jīng)網(wǎng)絡(luò)權(quán)值和閾值組成的向量,eq(w)為誤差。設(shè)wk表示第k次迭代的權(quán)值和閾值所組成的向量,wk+1=wk+Δw,運(yùn)用Levenberg-Marquardt算法對(duì)網(wǎng)絡(luò)的連接權(quán)值和閾值進(jìn)行優(yōu)化,權(quán)值與閾值組成向量的增量Δw計(jì)算公式為:

式中,J(w)為Jacobian矩陣,I為單位矩陣,μ為用戶定義的學(xué)習(xí)率。

4 變幅液壓系統(tǒng)故障診斷實(shí)驗(yàn)

通過搭建變幅液壓系統(tǒng)故障診斷實(shí)驗(yàn)臺(tái)以實(shí)現(xiàn)變幅液壓系統(tǒng)故障診斷實(shí)驗(yàn)數(shù)據(jù)采集,為基于故障樹的模糊神經(jīng)網(wǎng)絡(luò)各輸入變量隸屬函數(shù)參數(shù)值的計(jì)算及故障診斷提供實(shí)測(cè)數(shù)據(jù),對(duì)基于故障樹的模糊神經(jīng)網(wǎng)絡(luò)故障診斷方法的正確性和有效性進(jìn)行驗(yàn)證。

變幅液壓系統(tǒng)變幅無動(dòng)作實(shí)驗(yàn)信號(hào)采集傳感器包括4個(gè)壓力傳感器,由變幅無動(dòng)作故障樹可知采集的信號(hào)變量包括:PP-主泵壓力;PA-變幅小腔壓力;PB-變幅大腔壓力;PC-平衡閥控制油口壓力,壓力傳感器外形及安裝如圖5所示。

圖5 壓力傳感器及安裝

系統(tǒng)軟件的實(shí)現(xiàn)采用LabVIEW軟件的特色技巧——狀態(tài)機(jī)模式。關(guān)于狀態(tài)機(jī)的一個(gè)極度確切的描述是它是一個(gè)有向圖形,由一組節(jié)點(diǎn)和一組相應(yīng)的轉(zhuǎn)移函數(shù)組成。將流程圖各個(gè)過程轉(zhuǎn)換為狀態(tài)圖,采用狀態(tài)機(jī)結(jié)構(gòu)實(shí)現(xiàn)了上述程序流程圖。數(shù)據(jù)采集程序框圖如圖6所示。數(shù)據(jù)采集軟件圖形界面如圖7所示。

圖6 變幅液壓系統(tǒng)故障診斷數(shù)據(jù)采集程序框圖

圖7 變幅液壓系統(tǒng)故障診斷數(shù)據(jù)采集程序框圖

5 變幅液壓系統(tǒng)故障診斷實(shí)例分析

(1)變幅液壓系統(tǒng)故障診斷模糊隸屬度函數(shù)的參數(shù)計(jì)算

在實(shí)施中,在線采集檢測(cè)信號(hào)進(jìn)行故障診斷,利用式(4)至(6),可分別計(jì)算基于故障樹的模糊神經(jīng)網(wǎng)絡(luò)四個(gè)輸入變量對(duì)應(yīng)的隸屬度函數(shù)的參數(shù)值。

①主泵壓力 PP,Yˉ=71.176 3 bar,σ=5.597 6,隸屬度函數(shù)參數(shù)值為:{a,b,c,d}={54.383 5,65.578 7,76.773 9,87.969 2}。

②變幅小腔壓力PA,Yˉ=3.907 2 bar,σ=0.009 5,隸屬度函數(shù)參數(shù)值為:{a,b,c,d}={3.878 7,3.897 7,3.916 7,3.935 7}。

③變幅大腔壓力PB,Yˉ=53.712 8 bar,σ=0.025 6,隸屬度函數(shù)參數(shù)值為:{a,b,c,d}={53.63,593.687 1,53.738 4,53.789 7}。

④平衡閥控制油口壓力PC,Yˉ=21.095 5 bar,σ=0.062 3,隸屬度函數(shù)參數(shù)值為:{a,b,c,d}={20.908 7,21.033 2,21.157 7,21.282 3}。

(2)變幅液壓系統(tǒng)故障診斷神經(jīng)網(wǎng)絡(luò)權(quán)值和閾值的優(yōu)化訓(xùn)練

對(duì)每個(gè)輸入變量取模糊子集0.1、0.5、0.9分別表示壓力過低、正常、過高。由此提取神經(jīng)網(wǎng)絡(luò)權(quán)值及閾值的訓(xùn)練樣本,如表2所示。

表2 變幅無動(dòng)作神經(jīng)網(wǎng)絡(luò)訓(xùn)練樣本

運(yùn)用Levenberg-Marquardt算法對(duì)神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練,神經(jīng)網(wǎng)絡(luò)輸入節(jié)點(diǎn)為4個(gè),輸出節(jié)點(diǎn)為6個(gè),根據(jù)Kolmogorov定理可確定隱含層的節(jié)點(diǎn)數(shù)為8個(gè),訓(xùn)練結(jié)果如圖8所示,從圖可以看出,神經(jīng)網(wǎng)絡(luò)經(jīng)過5次訓(xùn)練即已達(dá)到期望目標(biāo)。

圖8 目標(biāo)函數(shù)訓(xùn)練曲線

Levenberg-Marquardt法實(shí)際上是梯度下降法和牛頓法的結(jié)合,它的優(yōu)點(diǎn)在于網(wǎng)絡(luò)權(quán)值數(shù)目較少時(shí)收斂非常迅速。應(yīng)用Levenberg-Marquardt優(yōu)化算法比傳統(tǒng)的BP及其他改進(jìn)算法迭代次數(shù)少,收斂速度快,精確度高。因此Levenberg-Marquardt優(yōu)化算法與其他訓(xùn)練算法相比在BP網(wǎng)絡(luò)學(xué)習(xí)中具有一定優(yōu)越性,神經(jīng)網(wǎng)絡(luò)訓(xùn)練算法比較如表3所示。

表3 神經(jīng)網(wǎng)絡(luò)訓(xùn)練算法比較

為了驗(yàn)證故障診斷方法的有效性,將變幅無動(dòng)作下采集的故障數(shù)據(jù){PP,PA,PB,PC}={70.536 8 bar,4.015 6 bar,53.674 2 bar,4.018 9 bar}作為模糊神經(jīng)網(wǎng)絡(luò)的輸入,將輸入數(shù)據(jù)進(jìn)行模糊化處理,求得各輸入變量對(duì)應(yīng)論域的隸屬度為 {0.124 7,0.374 8,0.783 0,0.421 5},經(jīng)過網(wǎng)絡(luò)計(jì)算得到輸出為:{0.001 8,0.001 3,0.001 4,0.002 5,0.998 7,0.001 6},由表2中故障模式描述可知,該輸出與故障模式中的F5對(duì)應(yīng)程度非常高,即X6控制器或手柄損壞或X7先導(dǎo)比例減壓閥損壞,實(shí)際狀況是操作手柄與變幅換向閥斷接導(dǎo)致控制腔壓力未建立,診斷結(jié)果與實(shí)際狀況一致。

6 結(jié)束語

變幅液壓系統(tǒng)是各機(jī)構(gòu)控制回路中最復(fù)雜的系統(tǒng)之一,傳統(tǒng)的故障推理算法難以滿足這類復(fù)雜系統(tǒng)故障診斷的要求。本文提出了一種基于故障樹的模糊神經(jīng)網(wǎng)絡(luò)故障推理機(jī)制,該機(jī)制融合了故障樹和模糊神經(jīng)網(wǎng)絡(luò)兩種方法的優(yōu)點(diǎn),研究了變幅液壓系統(tǒng)模糊神經(jīng)網(wǎng)絡(luò)故障診斷推理機(jī)制的研制過程,以頂事件變幅無動(dòng)作為對(duì)象,通過實(shí)例驗(yàn)證了該故障診斷方法的可行性和有效性。

[1]吳宗彥,韓煜,張建軍,等.基于模糊神經(jīng)網(wǎng)絡(luò)的自動(dòng)生產(chǎn)線故障診斷方法研究[J].中國(guó)機(jī)械工程,2008(9):1072-1075.

[2]Yznaga Blanco I,Sun D,He Y K.Neural-fuzzy technique for inverter faults diagnosis in PMSM DTC system[J].Electric Machines and Control,2008,12(3):132-139.

[3]張敬芬,孟光,趙德有.基于模糊神經(jīng)網(wǎng)絡(luò)的薄板不同指標(biāo)裂紋診斷[J].機(jī)械工程學(xué)報(bào),2006(3):145-149.

[4]陳煒,吳志良.故障診斷專家系統(tǒng)結(jié)合故障樹技術(shù)的結(jié)構(gòu)研究[J].航海技術(shù),2005(6):47-49.

[5]艾芳菊.模糊神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)優(yōu)化研究[D].成都:中國(guó)科學(xué)院研究生院成都計(jì)算機(jī)應(yīng)用研究所,2006.

[6]姜長(zhǎng)元.模糊神經(jīng)網(wǎng)絡(luò)模型及其應(yīng)用研究[D].南京:南京師范大學(xué),2005.

[7]彭廣平,李泳鮮.模糊神經(jīng)網(wǎng)絡(luò)在汽車發(fā)動(dòng)機(jī)電子點(diǎn)火系統(tǒng)故障診斷中的應(yīng)用[J].汽車科技,2011(3):47-50.

[8]孫增國(guó).神經(jīng)網(wǎng)絡(luò)和模糊專家系統(tǒng)在故障診斷中的應(yīng)用[D].遼寧大連:大連理工大學(xué),2004.

[9]張婷.模糊神經(jīng)網(wǎng)絡(luò)在冷水機(jī)組故障診斷中的應(yīng)用研究[D].天津:天津大學(xué),2009.

[10]Pltkeisheim F.The three sigma rule[J].The American Statistician,1994,48(2):88-91.

[11]趙偉.船舶柴油機(jī)監(jiān)測(cè)與故障診斷技術(shù)研究[D].遼寧大連:大連海事大學(xué),2006.

[12]李瓊.信息融合技術(shù)在航空發(fā)動(dòng)機(jī)智能監(jiān)測(cè)中的應(yīng)用研究[D].南京:南京航空航天大學(xué),2011.

[13]楊柳,陳艷萍.一種新的Levenberg-Marquardt算法的收斂性[J].計(jì)算數(shù)學(xué),2005(1):55-62.

[14]毛明明,柳益君,湯嘉立.基于L-M神經(jīng)網(wǎng)絡(luò)的齒輪故障診斷[J].計(jì)算機(jī)技術(shù)與發(fā)展,2011(1):210-213.

[15]柳益君,古春生,朱廣萍,等.LM神經(jīng)網(wǎng)絡(luò)在環(huán)境質(zhì)量評(píng)估中的應(yīng)用[J].江南大學(xué)學(xué)報(bào):自然科學(xué)版,2010(2):213-216.

猜你喜歡
故障診斷故障系統(tǒng)
Smartflower POP 一體式光伏系統(tǒng)
WJ-700無人機(jī)系統(tǒng)
ZC系列無人機(jī)遙感系統(tǒng)
故障一點(diǎn)通
連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
奔馳R320車ABS、ESP故障燈異常點(diǎn)亮
因果圖定性分析法及其在故障診斷中的應(yīng)用
故障一點(diǎn)通
江淮車故障3例
基于LCD和排列熵的滾動(dòng)軸承故障診斷
主站蜘蛛池模板: 国产AV毛片| 欧美视频在线不卡| 精品人妻一区二区三区蜜桃AⅤ| 久久黄色视频影| 国产麻豆精品久久一二三| 欧美区在线播放| 欧美成人一级| 色老二精品视频在线观看| 国产一级精品毛片基地| 亚洲黄色网站视频| 欧美无专区| 青草视频久久| WWW丫丫国产成人精品| 亚洲永久色| 久久精品亚洲热综合一区二区| 国产大全韩国亚洲一区二区三区| 亚洲最大情网站在线观看| 国模极品一区二区三区| 国产导航在线| 五月婷婷综合色| 久久国产亚洲偷自| Jizz国产色系免费| 亚洲一区无码在线| 天天色综网| 免费毛片全部不收费的| 国模视频一区二区| 亚洲一区二区三区国产精品| 日韩欧美国产精品| 免费观看无遮挡www的小视频| 精品少妇人妻一区二区| 国产女人在线观看| 欧美一区二区三区国产精品| 日韩无码黄色| 国产网站免费观看| 亚洲永久精品ww47国产| 一级成人a做片免费| 亚洲清纯自偷自拍另类专区| 亚洲精品国产精品乱码不卞| 亚洲免费三区| 亚洲精品男人天堂| 亚洲v日韩v欧美在线观看| 国产网站一区二区三区| 一级毛片在线免费视频| 亚洲色大成网站www国产| 久久毛片基地| 国产精品第5页| 国产国产人免费视频成18| 亚洲无码高清一区| 国产在线观看第二页| 97在线观看视频免费| 亚洲中文字幕97久久精品少妇| 欧美激情伊人| 国产精品hd在线播放| 亚洲欧洲日韩综合色天使| 精品一区二区久久久久网站| 狠狠亚洲婷婷综合色香| 97青青青国产在线播放| 亚洲综合国产一区二区三区| 精品久久久无码专区中文字幕| 99尹人香蕉国产免费天天拍| 国产系列在线| 伊人久久影视| 最新亚洲人成无码网站欣赏网 | 亚洲无码久久久久| 四虎成人免费毛片| 91探花在线观看国产最新| 波多野结衣在线se| 欧美性精品| 欧美成人第一页| 国产综合在线观看视频| 香蕉久久国产超碰青草| 欧美精品啪啪一区二区三区| 91精品专区国产盗摄| 99这里只有精品免费视频| 77777亚洲午夜久久多人| A级毛片无码久久精品免费| 黄色网站在线观看无码| a级毛片免费播放| 在线色国产| 久久这里只精品国产99热8| 婷婷六月激情综合一区| 国产亚洲视频在线观看|