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

    定值和隨動單閉環(huán)系統(tǒng)傳感器故障診斷

    2020-10-12 14:42:32那文波劉甜甜
    控制理論與應(yīng)用 2020年9期
    關(guān)鍵詞:故障診斷故障方法

    那文波,高 宇,李 明,劉甜甜

    (中國計量大學(xué)機電工程學(xué)院,浙江杭州 310000)

    1 引言

    定值和隨動單回路閉環(huán)控制系統(tǒng)在生產(chǎn)中應(yīng)用廣泛,大致占控制系統(tǒng)的90%左右.傳感器處在反饋回路當中,是構(gòu)成閉環(huán)控制的主要單元,由于閉環(huán)系統(tǒng)的實時控制作用,使得其故障診斷變得復(fù)雜[1–2].

    故障診斷技術(shù)可以分為3類方法[3–4].第1類方法是基于解析模型法,主要針對于可以建立精確數(shù)學(xué)模型的控制對象,例如參數(shù)估計法[5];第2類方法是基于定性經(jīng)驗法,通過先驗知識建立起非量化的故障判斷依據(jù),例如專家系統(tǒng)[6]和故障元等效樹法[7];第3類方法是基于數(shù)據(jù)驅(qū)動法,不依賴于控制對象的解析模型,利用實時數(shù)據(jù)在線處理分析故障類型,或者根據(jù)大量已知樣本總結(jié)出故障特征從而進行推理,例如卡爾曼濾波[8]、神經(jīng)網(wǎng)絡(luò)[9].

    目前基于數(shù)據(jù)驅(qū)動的故障診斷方法是研究熱點[10].如文獻[11]根據(jù)小波分析思想構(gòu)造新型機械設(shè)備故障診斷裝置;如文獻[12]采用學(xué)習(xí)矢量量化(learning vector quantization,LVQ)神經(jīng)網(wǎng)絡(luò)作為模式分類器進行故障處理;如文獻[13]利用主成分分析(principal component analysis,PCA)方法以及神經(jīng)網(wǎng)絡(luò)方法,對R¨ossler混沌方程預(yù)測仿真算法.

    針對閉環(huán)系統(tǒng)傳感器故障診斷問題,如文獻[14]提出一種基于子空間識別和卡爾曼濾波技術(shù),采用平均值推導(dǎo)、廣義似然比檢測(generalized x likelihood ratio test,GLRT)等方法處理殘差來進行故障檢測,但不能較好的處理誤差以及檢測慢漂移故障;文獻[15]針對閉環(huán)系統(tǒng)傳感器故障診斷問題,提出一種用于四旋翼無人飛行器加速度計和陀螺儀故障同時發(fā)生的故障檢測與隔離以及故障偏差值估計的非線性診斷方法;文獻[16]針對閉環(huán)系統(tǒng)中的傳感器故障,提出了基于Kullback-Leibler(KL)距離的微小故障在線檢測與估計方法;文獻[17]分析故障情況下系統(tǒng)殘差的一步轉(zhuǎn)移方程遞推表達式,研究閉環(huán)系統(tǒng)故障情況下估計誤差信號和殘差信號的動態(tài)演化規(guī)律;文獻[18]針對牽引電機閉環(huán)控制系統(tǒng),考慮電機運行過程中的諧波干擾情況,基于閉環(huán)控制策略的故障診斷方法對傳感器進行故障診斷,但該方法若對轉(zhuǎn)速變化的復(fù)雜瞬態(tài)過程中發(fā)生的故障進行診斷,則需要更為精確的數(shù)學(xué)模型.

    本文的主要貢獻:一是提出了一種適用于定值和隨動單回路閉環(huán)系統(tǒng)的傳感器故障診斷方法;二是基于定值和隨動單回路閉環(huán)系統(tǒng)正常和故障動態(tài)特性分析和數(shù)據(jù)驅(qū)動方法,建立了故障診斷模型,并給出了利用實時數(shù)據(jù)標定故障診斷模型的流程;三是所提出的方法,一定程度上抑制了隨機噪聲和系統(tǒng)及環(huán)境固有噪聲的影響,適應(yīng)緩慢漂移故障的診斷,診斷時間短實時性強,并方便工程組態(tài)實現(xiàn),通過實驗驗證了方法的可行性和故障診斷的高準確度.

    2 控制系統(tǒng)正常和故障動態(tài)特性分析

    本文以單容對象的定值和隨動單回路閉環(huán)系統(tǒng)為例,如圖1所示.其中:Ksor是傳感器無故障時放大系數(shù);kKsor為乘性故障模塊,k為乘性故障增益;Ksor+a為加性故障模塊,a為加性故障偏差;通過模擬機械開關(guān)實現(xiàn)了乘性故障模塊、無故障模塊和加性故障模塊之間的切換.

    圖1 單回路控制系統(tǒng)方框圖Fig.1 Single loop control system block diagram

    系統(tǒng)無故障狀態(tài)下的傳遞函數(shù):

    其中:G(s)=G1(s)G2(s)KV;定值系統(tǒng)輸入信號R(s)為一恒定值,隨動系統(tǒng)中R(s)為一隨時間變化值,從上述式(1)–(3)可以看出當不同系統(tǒng)發(fā)生不同故障時,系統(tǒng)的動態(tài)響應(yīng)不同,下面介紹利用不同故障時系統(tǒng)的動態(tài)響應(yīng)特征的故障診斷方法.

    3 故障診斷方法

    以傳感器為故障診斷對象,基于動態(tài)數(shù)據(jù)流實時趨勢的特征值提取建立故障診斷模型.解決方案具體分為故障監(jiān)測、故障估計和故障分離.

    3.1 故障監(jiān)測

    數(shù)據(jù)殘差是一個判定標準,可以用于統(tǒng)計系統(tǒng)發(fā)生故障時的邊界值

    其中:e為殘差,A為可觀察數(shù)據(jù)值,A?為期望數(shù)據(jù)值.當系統(tǒng)采樣周期為T,數(shù)據(jù)點總數(shù)為N,選取數(shù)據(jù)窗口大小t(數(shù)據(jù)點n),則有

    采用堆棧結(jié)構(gòu)將新采集的動態(tài)數(shù)據(jù)填充進窗口并剔除最早的數(shù)據(jù),設(shè)共有m個數(shù)據(jù)窗口,則

    數(shù)據(jù)的窗口殘差和由該窗口內(nèi)所有數(shù)據(jù)點殘差相加取絕對值得到,設(shè)Ei為第i個窗口殘差和,那么

    式中:Ai為第i個實際數(shù)據(jù)點,為第i個期望數(shù)據(jù)點.

    假設(shè)正常數(shù)據(jù)殘差和滿足正態(tài)分布,采集l組正常數(shù)據(jù)的最大殘差和,計算l組數(shù)據(jù)最大殘差和平均值

    標準差為

    將2σ作為去除異常數(shù)據(jù)的方法,則故障發(fā)生時殘差和

    故障監(jiān)測方法:即是將窗口殘差和與閾值作比較,若超出了閾值范圍,則判斷系統(tǒng)發(fā)生了故障,反之系統(tǒng)屬于正常運行狀態(tài).

    3.2 故障估計

    根據(jù)故障監(jiān)測可以將有效診斷數(shù)據(jù)范圍隔離在一個數(shù)據(jù)窗口內(nèi),即得到故障窗口數(shù)據(jù),其中變化率最大(i)的點i對應(yīng)的時刻即為故障發(fā)生時刻.

    定值系統(tǒng)加性故障估計方法:當系統(tǒng)的故障類型為加性故障時,此時加性故障偏差a為

    當a的符號為正時,系統(tǒng)故障值為增量,反之為減量.

    定值系統(tǒng)乘性故障估計方法:當系統(tǒng)的故障類型為乘性故障時,此時乘性故障增益k為

    當k >1時,系統(tǒng)故障值為增量;當k <1時,系統(tǒng)故障值為減量.

    式(11)–(12)即為定值系統(tǒng)故障估計模型.

    隨動系統(tǒng)加性故障估計方法:當系統(tǒng)的故障類型為加性故障時,此時加性故障偏差a為

    隨動系統(tǒng)乘性故障估計方法:當系統(tǒng)的故障類型為乘性故障時,此時乘性故障增益k為

    式(13)–(14)即為隨動系統(tǒng)故障估計模型.

    3.3 故障分離

    系統(tǒng)發(fā)生不同強度的同類故障或發(fā)生相同強度的不同類型故障時,系統(tǒng)的初始響應(yīng)斜率均不同,基于此實現(xiàn)故障分離.

    由最小二乘法擬合方程

    式(15)為一元線性回歸方程模型,其中β0和β1為回歸方程系數(shù).

    選取故障數(shù)據(jù)窗故障發(fā)生時刻開始的數(shù)據(jù)變化最為明顯的數(shù)據(jù)集,根據(jù)最小二乘法得到這段時間內(nèi)數(shù)據(jù)動態(tài)趨勢的擬合方程

    對系統(tǒng)的加性故障狀態(tài)和乘性故障狀態(tài)不同強度的模擬動態(tài)數(shù)據(jù)擬合處理,將式(16)中一次項系數(shù)作為提取到的特征值EV:

    其中:P為加性故障特征值,Q為乘性故障特征值.再次采用一元線性回歸法對不同故障類型的一次項系數(shù)進行擬合,以加性故障偏差作為自變量,加性故障一次項系數(shù)作為因變量,得到加性故障特征值擬合方程

    同樣得到乘性故障特征值擬合方程

    在系統(tǒng)發(fā)生故障時,基于故障估計分別得到加性和乘性的故障強度,將加性故障偏差和乘性故障增益代入相應(yīng)擬合方程中得到相應(yīng)的.

    式(16)(18)–(19)即為故障分離的靜態(tài)模型,實際應(yīng)用時要結(jié)合實際系統(tǒng)在線數(shù)據(jù)進行標定.

    較小者視為有效結(jié)果,由此判斷故障類型為加性故障或者乘性故障.

    式(20)即為故障分離的診斷模型.

    3.4 在線標定及故障診斷流程

    實際工程應(yīng)用中傳感器系統(tǒng)在線標定及故障診斷的流程圖如圖2所示.

    4 驗證實例

    驗證基于四容水箱的“復(fù)雜過程系統(tǒng)創(chuàng)新實驗平臺”,該平臺將OPC服務(wù)器連接到下位機PLC,MATLAB/Simulink虛擬控制器能夠通過PID運算進行定值或隨動單容液位對象的單回路控制,結(jié)合組態(tài)Wincc監(jiān)控界面的實時監(jiān)控和MATLAB中Workspace的在線動態(tài)數(shù)據(jù)進行實時數(shù)據(jù)分析和建立診斷模型.通過調(diào)用Simulink函數(shù)編寫程序生成M文件動態(tài)鏈接到MATLAB軟件中實現(xiàn)故障診斷方法的驗證.通過實驗可知一般在線診斷時間為200 ms.

    4.1 在線數(shù)據(jù)采集

    實驗中設(shè)最終期望液位值為10 cm,采樣周期Ts=0.5 s,單次實驗中設(shè)備運行時間設(shè)定為1000 s.定值系統(tǒng)設(shè)定輸入幅值為10 cm的階躍信號,零時刻開始輸入信號;隨動系統(tǒng)設(shè)定輸入斜率為0.01 cm/s的斜坡信號,零時刻開始輸入初始值為零的信號.

    4.1.1 加性故障狀態(tài)

    選取500 s時刻加入大小分別為±2.0 cm范圍內(nèi)不同強度的加性故障.如圖3所示為隨動控制系統(tǒng),以加性故障偏差為?1.0 cm為例,得到加性故障狀態(tài)下的系統(tǒng)響應(yīng)曲線.

    圖2 在線標定及故障診斷流程圖Fig.2 Flow chart of on-line calibration and fault diagnosis

    圖3 隨動控制系統(tǒng)加性故障?1.0 cm液位響應(yīng)曲線Fig.3 Additive fault ?1.0 cm level response curve of servo control system

    4.1.2 乘性故障狀態(tài)

    在其他參數(shù)保持不變的情況下,對于乘性故障,同樣在500 s時加入增益分別為0.8~1.2范圍內(nèi)不同強度的乘性故障,可以采集到乘性故障數(shù)據(jù).

    4.2 實時故障監(jiān)測

    以隨動系統(tǒng)為例,觀察各種類型故障響應(yīng)曲線可知系統(tǒng)從開始發(fā)生故障到系統(tǒng)自身調(diào)節(jié)到設(shè)定值附近的時間為25 s,設(shè)定t=25 s為數(shù)據(jù)窗口大小,則每個窗口數(shù)據(jù)點n為n=t/T=50.

    一次實驗有N=2001個數(shù)據(jù),共有m個數(shù)據(jù)窗口

    根據(jù)加窗殘差和的方法,求得多組正常無故障狀態(tài)下數(shù)據(jù)的最大窗口殘差和,為集合EMax:

    根據(jù)式(8)–(9)得到這10組數(shù)據(jù)的平均值及標準差

    由式(10)確定故障閾值E?=4.6259.

    對故障診斷方法靈敏度做邊界測試,確定實時數(shù)據(jù)的漏檢情況.同時對多組不同故障強度的情況進行故障監(jiān)測,在隨動控制系統(tǒng)中當加性故障偏差在設(shè)定值的0.8%范圍外時,乘性故障增益在設(shè)定值的1.8%范圍外時,故障監(jiān)測均未出現(xiàn)漏檢情況.

    對于定值控制系統(tǒng),這種方法同樣可以達到故障監(jiān)測的目的,加性故障和乘性故障大小均在設(shè)定值的0.8%范圍外時,故障監(jiān)測均未出現(xiàn)漏檢情況.

    4.3 實時故障估計

    以隨動系統(tǒng)為例,當發(fā)生乘性故障1.20時的變化率曲線如圖4所示,可以判斷故障發(fā)生時間約為500 s.

    圖4 隨動控制系統(tǒng)故障狀態(tài)下液位變化率Fig.4 Level change rate of servo control system

    在期望值的5%范圍內(nèi),選取多組加性故障和乘性故障進行實驗,表1和表2給出了故障估計驗證結(jié)果.

    表1 加性故障估計驗證Table 1 Additive fault verification

    表2 乘性故障估計驗證Table 2 Multiplicative fault verification

    計算可得加性故障定值精度2.0%,乘性故障定值精度2.2%.

    對于定值系統(tǒng),同樣在期望值的范圍內(nèi)進行故障估計實驗,得到加性故障定值精度1.8%,乘性故障定值精度2.0%.

    一般的系統(tǒng)在強故障情況下診斷精度高,可以看出在弱到5%故障強度范圍內(nèi)也有很高的診斷精度.

    4.4 實時故障分離

    定值和隨動系統(tǒng)故障分離的區(qū)別在于,故障發(fā)生后系統(tǒng)調(diào)節(jié)過程的特征不相同,需要根據(jù)實際的實驗結(jié)果分別得到定值和隨動系統(tǒng)的擬合曲線,進行故障分離.當確定故障發(fā)生后,根據(jù)最小二乘法計算隨動控制系統(tǒng)中不同類型故障下發(fā)生不同強度故障時的數(shù)據(jù)動態(tài)趨勢擬合方程,如表3、表4、表5和表6所示為32組實驗的擬合方程的一次項擬合系數(shù)以及擬合曲線的RMSE.

    以表3為例,根據(jù)加性故障偏差a<0的一次項擬合系數(shù)得到不同程度故障狀態(tài)下的擬合函數(shù)曲線如圖5所示.同樣對表4、表5和表6中的數(shù)據(jù)分別進行二次擬合,如表7不同故障狀態(tài)下擬合函數(shù).

    根據(jù)表7的擬合結(jié)果進行分析,r2→1,RMSE →0,數(shù)據(jù)較好地擬合了對應(yīng)的故障響應(yīng).如表8和表9所示為故障分離方法的驗證結(jié)果.

    表3 加性故障偏差小于0的一次項擬合系數(shù)Table 3 Fitting coefficient of first term with additive fault deviation less than 0

    表4 加性故障偏差大于0的一次項擬合系數(shù)Table 4 Fitting coefficient of first term with additive fault deviation greater than 0

    表5 乘性故障增益小于1的一次項擬合系數(shù)Table 5 Fitting coefficient of first term with multiplicative fault gain less than 1

    表6 乘性故障增益大于1的一次項擬合系數(shù)Table 6 Fitting coefficient of first term with multiplicative fault gain greater than 1

    表7 二次擬合結(jié)果分析Table 7 Analysis of quadratic fitting results

    圖5 同程度故障狀態(tài)下的二次擬合曲線Fig.5 Quadratic fitting curve under different fault conditions

    表8 加性故障分離驗證結(jié)果Table 8 Validation results of additive fault isolation

    表9 乘性故障分離驗證結(jié)果Table 9 Validation results of multiplicative fault isolation

    可見,當隨動控制系統(tǒng)傳感器發(fā)生的故障大小在偏離期望液位值的2%以內(nèi)時,不能分離.對于定值系統(tǒng),同樣驗證得到故障在偏離期望液位值的1.5%的范圍內(nèi),不能分離.

    綜上所述,基于動態(tài)趨勢的隨動系統(tǒng)傳感器故障診斷的故障檢測強度死區(qū)是1.8%,故障估計精度為2.2%,故障分離死區(qū)是2.0%;對于定值系統(tǒng),故障診斷的故障檢測強度死區(qū)是0.8%,故障估計精度為2.0%,故障分離死區(qū)是1.5%.驗證了方法的有效性和診斷的高精度[17].

    5 結(jié)語

    定值和隨動單閉環(huán)系統(tǒng)傳感器故障診斷,以一階定值和隨動閉環(huán)控制系統(tǒng)中的傳感器為故障診斷對象,基于數(shù)據(jù)驅(qū)動的思想,利用系統(tǒng)正常和故障動態(tài)特性分析,建立了系統(tǒng)故障監(jiān)測、故障估計和故障分離的靜態(tài)模型,并給出了在線應(yīng)用的診斷模型標定方法和診斷流程.通過實驗驗證了方法的有效性和診斷的高精度.

    方法適用于一般的一階閉環(huán)定值或隨動控制系統(tǒng)傳感器的故障診斷;結(jié)合在線仿真數(shù)據(jù)標定故障診斷靜態(tài)模型,不需建立系統(tǒng)模型,避免了系統(tǒng)模型建立不精確的缺陷,同時也避免了環(huán)境干擾和系統(tǒng)其它組成部分特性的影響.考慮系統(tǒng)各組成部分磨損,也可采取周期性標定的方法,提高診斷系統(tǒng)的準確度.

    工程中的噪聲對故障診斷是有影響的[19],本文研究方法是基于數(shù)據(jù)動態(tài)趨勢分析,選擇數(shù)據(jù)窗口的長度覆蓋了故障發(fā)生和隨之的調(diào)節(jié)過程,監(jiān)測閾值的確定采用多組數(shù)據(jù)最大殘差和平均值結(jié)合標準差來擇優(yōu)確定,且各在線診斷模型均采用實際工程系統(tǒng)在線數(shù)據(jù)標定,對隨機噪聲和系統(tǒng)及環(huán)境固有噪聲起到了一定程度的抑制作用,獲得了較高的診斷準確度.如何進一步抑制噪聲影響,進一步提高診斷準確度是下一步重點研究的問題.

    工程對故障診斷的實時性有一定要求[20],造成時間滯后的原因主要有傳感器測量時間、數(shù)據(jù)傳輸時間和故障診斷算法處理時間,本方法在實驗系統(tǒng)實時故障診斷時間為200 ms,基本能滿足一般工程系統(tǒng)的實時性需要.進一步提高故障診斷的實時性,也是下一步研究工作的重點之一.

    猜你喜歡
    故障診斷故障方法
    故障一點通
    可能是方法不對
    奔馳R320車ABS、ESP故障燈異常點亮
    用對方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    捕魚
    故障一點通
    江淮車故障3例
    基于LCD和排列熵的滾動軸承故障診斷
    亚洲婷婷狠狠爱综合网| 看免费成人av毛片| 又大又黄又爽视频免费| 婷婷色综合www| 日韩不卡一区二区三区视频在线| 在线观看一区二区三区激情| 18禁国产床啪视频网站| 一区在线观看完整版| 精品亚洲乱码少妇综合久久| 美国免费a级毛片| 视频区图区小说| 成年人免费黄色播放视频| 一边摸一边抽搐一进一出视频| 99香蕉大伊视频| 久久毛片免费看一区二区三区| 亚洲综合精品二区| 国产一区亚洲一区在线观看| 亚洲av国产av综合av卡| 熟女av电影| 建设人人有责人人尽责人人享有的| 久久热在线av| 亚洲av中文av极速乱| 国产又色又爽无遮挡免| 丝袜脚勾引网站| 99热网站在线观看| av有码第一页| 女人高潮潮喷娇喘18禁视频| 男女国产视频网站| 久久久久久久久久久久大奶| 国产色婷婷99| 亚洲伊人色综图| 国产av国产精品国产| 亚洲国产av影院在线观看| 亚洲av国产av综合av卡| 久久免费观看电影| 免费观看av网站的网址| 日韩,欧美,国产一区二区三区| 一区福利在线观看| 老熟女久久久| 老司机靠b影院| 国产精品99久久99久久久不卡 | 黑人猛操日本美女一级片| 亚洲天堂av无毛| 亚洲一区二区三区欧美精品| 男女下面插进去视频免费观看| 老汉色av国产亚洲站长工具| 成人免费观看视频高清| 深夜精品福利| 亚洲精华国产精华液的使用体验| 国产日韩欧美在线精品| 亚洲美女搞黄在线观看| 久热爱精品视频在线9| 国产1区2区3区精品| 国产精品 欧美亚洲| 可以免费在线观看a视频的电影网站 | 女人久久www免费人成看片| 欧美黑人精品巨大| 国产深夜福利视频在线观看| 久久久久久久大尺度免费视频| 最近手机中文字幕大全| 亚洲国产最新在线播放| 少妇人妻 视频| 男人舔女人的私密视频| 国产一区有黄有色的免费视频| 美女高潮到喷水免费观看| 亚洲五月色婷婷综合| 日韩中文字幕视频在线看片| 99精品久久久久人妻精品| 街头女战士在线观看网站| 午夜福利网站1000一区二区三区| 精品免费久久久久久久清纯 | 国产成人免费无遮挡视频| 精品第一国产精品| 国产麻豆69| 一区在线观看完整版| 日韩av不卡免费在线播放| 久久精品国产亚洲av高清一级| 波野结衣二区三区在线| 亚洲av日韩精品久久久久久密 | 无限看片的www在线观看| 亚洲色图综合在线观看| 一级黄片播放器| www.熟女人妻精品国产| 91国产中文字幕| 成人手机av| 男女边摸边吃奶| 天天躁夜夜躁狠狠久久av| 交换朋友夫妻互换小说| 你懂的网址亚洲精品在线观看| 亚洲精品美女久久av网站| 欧美黑人精品巨大| 秋霞伦理黄片| 亚洲国产毛片av蜜桃av| 999精品在线视频| 婷婷成人精品国产| 丁香六月天网| 亚洲国产成人一精品久久久| 久久精品久久久久久久性| 一级a爱视频在线免费观看| 婷婷色综合www| 一本久久精品| 国产精品熟女久久久久浪| 一本久久精品| 午夜91福利影院| 欧美 亚洲 国产 日韩一| 午夜福利影视在线免费观看| av国产久精品久网站免费入址| 高清不卡的av网站| 久久国产亚洲av麻豆专区| 国产99久久九九免费精品| 9色porny在线观看| 中国国产av一级| 两个人看的免费小视频| 青春草亚洲视频在线观看| 亚洲情色 制服丝袜| 久久97久久精品| 9191精品国产免费久久| 日韩 亚洲 欧美在线| 久久人人97超碰香蕉20202| 日韩一区二区视频免费看| 岛国毛片在线播放| 国产成人精品无人区| 别揉我奶头~嗯~啊~动态视频 | 成人午夜精彩视频在线观看| 19禁男女啪啪无遮挡网站| 亚洲成人av在线免费| 中文字幕另类日韩欧美亚洲嫩草| 最新的欧美精品一区二区| 亚洲成国产人片在线观看| 777米奇影视久久| 久久青草综合色| av不卡在线播放| 亚洲精品成人av观看孕妇| 中文字幕制服av| 久久狼人影院| 亚洲在久久综合| 免费高清在线观看日韩| 久久人妻熟女aⅴ| 久久精品亚洲av国产电影网| av片东京热男人的天堂| 中文字幕最新亚洲高清| 久久 成人 亚洲| 国产精品嫩草影院av在线观看| 久久人人爽人人片av| 人人妻,人人澡人人爽秒播 | 日日爽夜夜爽网站| 一区二区三区乱码不卡18| 97精品久久久久久久久久精品| 久久久久久人人人人人| 亚洲国产看品久久| 看十八女毛片水多多多| 电影成人av| 欧美人与性动交α欧美精品济南到| 色吧在线观看| 欧美另类一区| 啦啦啦在线免费观看视频4| 水蜜桃什么品种好| 中文字幕亚洲精品专区| 青春草视频在线免费观看| 九草在线视频观看| 男的添女的下面高潮视频| bbb黄色大片| √禁漫天堂资源中文www| 久久精品久久精品一区二区三区| 一本大道久久a久久精品| 青春草国产在线视频| 新久久久久国产一级毛片| 亚洲精品一二三| 五月开心婷婷网| 视频在线观看一区二区三区| 操美女的视频在线观看| 久久ye,这里只有精品| 欧美人与性动交α欧美软件| 老汉色∧v一级毛片| 一区福利在线观看| 免费高清在线观看日韩| 街头女战士在线观看网站| 大片电影免费在线观看免费| 国产1区2区3区精品| 少妇人妻精品综合一区二区| 在线观看人妻少妇| 国产精品二区激情视频| 日韩av不卡免费在线播放| 日日撸夜夜添| 亚洲精品国产色婷婷电影| 午夜91福利影院| 超碰97精品在线观看| 久久精品人人爽人人爽视色| 日韩av在线免费看完整版不卡| 亚洲精品国产一区二区精华液| 无限看片的www在线观看| 亚洲国产成人一精品久久久| 久久亚洲国产成人精品v| 国产精品 国内视频| 操美女的视频在线观看| 校园人妻丝袜中文字幕| 黄色视频在线播放观看不卡| 国产亚洲精品第一综合不卡| 日本黄色日本黄色录像| 你懂的网址亚洲精品在线观看| 日韩伦理黄色片| 丰满乱子伦码专区| 我要看黄色一级片免费的| 午夜精品国产一区二区电影| 精品人妻熟女毛片av久久网站| 汤姆久久久久久久影院中文字幕| 精品国产乱码久久久久久小说| 久久ye,这里只有精品| 9热在线视频观看99| 精品亚洲成国产av| 精品酒店卫生间| 男女边摸边吃奶| 一级毛片 在线播放| 精品少妇内射三级| 久久久精品区二区三区| 青春草视频在线免费观看| 免费高清在线观看日韩| 中文字幕人妻熟女乱码| 久久热在线av| 黄色视频在线播放观看不卡| 午夜激情av网站| 国产深夜福利视频在线观看| 在线观看一区二区三区激情| 19禁男女啪啪无遮挡网站| 婷婷色麻豆天堂久久| 成人国语在线视频| 亚洲国产欧美网| 欧美 日韩 精品 国产| 黄色一级大片看看| 免费观看a级毛片全部| 丰满迷人的少妇在线观看| 欧美日韩精品网址| 99久久人妻综合| 亚洲精品美女久久av网站| 18在线观看网站| 在线观看三级黄色| 免费在线观看视频国产中文字幕亚洲 | 午夜免费男女啪啪视频观看| 久久99一区二区三区| 欧美日韩一级在线毛片| 国产日韩一区二区三区精品不卡| 色婷婷久久久亚洲欧美| 亚洲精品中文字幕在线视频| 久久久欧美国产精品| 亚洲国产精品一区二区三区在线| 国产精品秋霞免费鲁丝片| 啦啦啦中文免费视频观看日本| 国产精品99久久99久久久不卡 | 日本wwww免费看| 少妇猛男粗大的猛烈进出视频| 精品国产超薄肉色丝袜足j| 亚洲,欧美精品.| 亚洲精品,欧美精品| 青春草亚洲视频在线观看| 国产精品免费视频内射| h视频一区二区三区| 欧美激情高清一区二区三区 | av在线app专区| 午夜福利在线免费观看网站| 少妇人妻久久综合中文| 国产精品女同一区二区软件| 少妇人妻精品综合一区二区| 国产av一区二区精品久久| 亚洲一区二区三区欧美精品| 日韩av免费高清视频| 又黄又粗又硬又大视频| 乱人伦中国视频| 亚洲男人天堂网一区| 欧美精品一区二区免费开放| 成人影院久久| 天堂俺去俺来也www色官网| 久久天堂一区二区三区四区| 视频在线观看一区二区三区| av一本久久久久| 日韩 亚洲 欧美在线| 成人午夜精彩视频在线观看| 中文天堂在线官网| 2018国产大陆天天弄谢| 夜夜骑夜夜射夜夜干| 亚洲精品一区蜜桃| 精品少妇黑人巨大在线播放| 国产 一区精品| 观看美女的网站| 飞空精品影院首页| 欧美黄色片欧美黄色片| 日韩一本色道免费dvd| 午夜久久久在线观看| 下体分泌物呈黄色| kizo精华| 波多野结衣一区麻豆| 中文字幕av电影在线播放| 久久久久久久国产电影| 国产精品.久久久| 亚洲一卡2卡3卡4卡5卡精品中文| av国产精品久久久久影院| 老司机靠b影院| 伦理电影大哥的女人| 日本wwww免费看| 中文字幕高清在线视频| 久久久久久人妻| 亚洲成人av在线免费| 日韩一卡2卡3卡4卡2021年| 国产成人a∨麻豆精品| 亚洲精品一区蜜桃| 国产有黄有色有爽视频| 在线观看三级黄色| 亚洲欧美日韩另类电影网站| 国产成人午夜福利电影在线观看| 亚洲精品,欧美精品| av国产精品久久久久影院| 韩国精品一区二区三区| 男女国产视频网站| 捣出白浆h1v1| 一区二区av电影网| 午夜91福利影院| 搡老岳熟女国产| 久久天躁狠狠躁夜夜2o2o | 2018国产大陆天天弄谢| 欧美精品亚洲一区二区| 一本一本久久a久久精品综合妖精| 日韩大码丰满熟妇| 国产av国产精品国产| 免费黄网站久久成人精品| 中文字幕另类日韩欧美亚洲嫩草| 夫妻午夜视频| 啦啦啦啦在线视频资源| 岛国毛片在线播放| 日韩大码丰满熟妇| 黑人欧美特级aaaaaa片| 国产精品蜜桃在线观看| 欧美日韩一区二区视频在线观看视频在线| 可以免费在线观看a视频的电影网站 | 18禁动态无遮挡网站| 亚洲人成77777在线视频| 又粗又硬又长又爽又黄的视频| 99热全是精品| 日韩熟女老妇一区二区性免费视频| 免费不卡黄色视频| 在线 av 中文字幕| 国产又色又爽无遮挡免| 99香蕉大伊视频| 最近手机中文字幕大全| kizo精华| 精品国产乱码久久久久久小说| 欧美乱码精品一区二区三区| 一区二区日韩欧美中文字幕| 亚洲色图 男人天堂 中文字幕| 色播在线永久视频| 国产免费又黄又爽又色| 国产黄频视频在线观看| 久久韩国三级中文字幕| 亚洲国产看品久久| 大香蕉久久成人网| 日日爽夜夜爽网站| 波多野结衣av一区二区av| 一级,二级,三级黄色视频| 日日爽夜夜爽网站| 美女午夜性视频免费| 亚洲精品日韩在线中文字幕| 亚洲色图综合在线观看| 国产亚洲精品第一综合不卡| 国产精品 国内视频| 久久97久久精品| 涩涩av久久男人的天堂| 一级片免费观看大全| 观看美女的网站| 国产亚洲精品第一综合不卡| 日韩av不卡免费在线播放| 亚洲色图综合在线观看| 国产又爽黄色视频| 又大又爽又粗| av国产久精品久网站免费入址| 考比视频在线观看| 视频区图区小说| 国产精品一区二区精品视频观看| 久久久久精品人妻al黑| 日韩中文字幕视频在线看片| 日韩中文字幕欧美一区二区 | 性高湖久久久久久久久免费观看| 亚洲成国产人片在线观看| 国产黄色视频一区二区在线观看| 不卡视频在线观看欧美| 18禁裸乳无遮挡动漫免费视频| 精品少妇一区二区三区视频日本电影 | 午夜福利,免费看| 亚洲av男天堂| 最近的中文字幕免费完整| 成人免费观看视频高清| 日韩av免费高清视频| 欧美 日韩 精品 国产| 满18在线观看网站| 高清不卡的av网站| 日韩精品免费视频一区二区三区| 欧美日韩福利视频一区二区| 老司机影院成人| 高清视频免费观看一区二区| 久久久精品国产亚洲av高清涩受| 美女国产高潮福利片在线看| 国产精品二区激情视频| 中国三级夫妇交换| 最近的中文字幕免费完整| 人体艺术视频欧美日本| 国产日韩一区二区三区精品不卡| 99久久人妻综合| 日韩av在线免费看完整版不卡| 男男h啪啪无遮挡| 国产一区有黄有色的免费视频| 亚洲熟女精品中文字幕| 亚洲精品日本国产第一区| 19禁男女啪啪无遮挡网站| 高清av免费在线| 熟妇人妻不卡中文字幕| 又大又爽又粗| 又黄又粗又硬又大视频| 赤兔流量卡办理| 国产成人免费观看mmmm| 免费看av在线观看网站| 天天操日日干夜夜撸| 热re99久久国产66热| 久久亚洲国产成人精品v| 国产成人欧美在线观看 | 各种免费的搞黄视频| 精品国产超薄肉色丝袜足j| 制服丝袜香蕉在线| 丁香六月欧美| 欧美变态另类bdsm刘玥| 国产成人91sexporn| 亚洲欧美激情在线| 看十八女毛片水多多多| 免费av中文字幕在线| 亚洲成人一二三区av| 男女之事视频高清在线观看 | 视频在线观看一区二区三区| 观看美女的网站| 天天影视国产精品| 99香蕉大伊视频| 热99国产精品久久久久久7| 亚洲在久久综合| 色综合欧美亚洲国产小说| 女人被躁到高潮嗷嗷叫费观| 纵有疾风起免费观看全集完整版| 亚洲色图综合在线观看| 美女中出高潮动态图| 成人手机av| 免费日韩欧美在线观看| 久久鲁丝午夜福利片| 一级毛片电影观看| 涩涩av久久男人的天堂| 欧美日韩视频精品一区| 成年动漫av网址| 国产黄频视频在线观看| 最近手机中文字幕大全| 高清视频免费观看一区二区| 狂野欧美激情性xxxx| 韩国高清视频一区二区三区| 亚洲av日韩在线播放| 国精品久久久久久国模美| 精品一区二区三卡| 日韩伦理黄色片| 免费不卡黄色视频| 精品一区二区三区av网在线观看 | 中文乱码字字幕精品一区二区三区| 黄色视频不卡| 男女高潮啪啪啪动态图| 久久久国产欧美日韩av| 如日韩欧美国产精品一区二区三区| 99九九在线精品视频| 精品一品国产午夜福利视频| 成人手机av| 一本久久精品| 亚洲久久久国产精品| 日韩大码丰满熟妇| 国产精品一区二区在线不卡| 性色av一级| 免费看不卡的av| 中文欧美无线码| 国产精品久久久久成人av| 在线精品无人区一区二区三| 日本vs欧美在线观看视频| a 毛片基地| 精品酒店卫生间| 欧美国产精品一级二级三级| 免费黄网站久久成人精品| 亚洲国产精品一区二区三区在线| 亚洲精品久久午夜乱码| av在线app专区| 十八禁人妻一区二区| 超碰成人久久| 精品人妻熟女毛片av久久网站| 国产97色在线日韩免费| 18禁动态无遮挡网站| 亚洲av中文av极速乱| 9热在线视频观看99| 叶爱在线成人免费视频播放| 午夜免费鲁丝| avwww免费| 狂野欧美激情性xxxx| 十八禁高潮呻吟视频| 国产精品无大码| 韩国av在线不卡| 纵有疾风起免费观看全集完整版| 中文字幕av电影在线播放| 亚洲久久久国产精品| 中国国产av一级| 蜜桃在线观看..| 色综合欧美亚洲国产小说| 亚洲精品一区蜜桃| 日韩伦理黄色片| 热re99久久精品国产66热6| 夫妻性生交免费视频一级片| 精品第一国产精品| 欧美av亚洲av综合av国产av | 国产一级毛片在线| 欧美日韩成人在线一区二区| 丝袜脚勾引网站| 一级片免费观看大全| 老司机在亚洲福利影院| 亚洲国产精品一区二区三区在线| 最近手机中文字幕大全| 国产有黄有色有爽视频| 女人爽到高潮嗷嗷叫在线视频| 老司机影院毛片| 午夜福利视频精品| 国产成人精品福利久久| 亚洲色图综合在线观看| 大香蕉久久网| 中文字幕高清在线视频| 性高湖久久久久久久久免费观看| 亚洲,一卡二卡三卡| 国产一区二区三区综合在线观看| 搡老乐熟女国产| 国产不卡av网站在线观看| 一边亲一边摸免费视频| 精品国产一区二区三区四区第35| 亚洲一级一片aⅴ在线观看| 亚洲成人免费av在线播放| 成人国语在线视频| 国产日韩欧美亚洲二区| 亚洲av男天堂| 丁香六月天网| 爱豆传媒免费全集在线观看| 免费黄频网站在线观看国产| 少妇人妻久久综合中文| 日日啪夜夜爽| av片东京热男人的天堂| 国产一区有黄有色的免费视频| 亚洲精品乱久久久久久| 99久久人妻综合| 国产国语露脸激情在线看| 老司机靠b影院| 99热国产这里只有精品6| 亚洲精品中文字幕在线视频| 欧美人与性动交α欧美精品济南到| 晚上一个人看的免费电影| 精品国产乱码久久久久久小说| 久久久久久久大尺度免费视频| 国产国语露脸激情在线看| 一区二区三区激情视频| 两性夫妻黄色片| 晚上一个人看的免费电影| 日韩伦理黄色片| 亚洲美女搞黄在线观看| 99热网站在线观看| 巨乳人妻的诱惑在线观看| 中文精品一卡2卡3卡4更新| 极品少妇高潮喷水抽搐| 亚洲欧美成人综合另类久久久| 国产精品一二三区在线看| 伊人亚洲综合成人网| 亚洲综合精品二区| 国产探花极品一区二区| 免费av中文字幕在线| 一个人免费看片子| 精品国产乱码久久久久久小说| 老司机在亚洲福利影院| 成人亚洲欧美一区二区av| 午夜免费鲁丝| 欧美变态另类bdsm刘玥| 久久久久精品国产欧美久久久 | 午夜av观看不卡| 欧美日韩亚洲高清精品| 日本wwww免费看| 亚洲欧美色中文字幕在线| 一级片'在线观看视频| 少妇被粗大的猛进出69影院| 中国国产av一级| 黄片小视频在线播放| 亚洲 欧美一区二区三区| 久久女婷五月综合色啪小说| 午夜福利影视在线免费观看| 免费观看人在逋| 丝袜脚勾引网站| 欧美亚洲 丝袜 人妻 在线| 欧美成人精品欧美一级黄| 日韩制服骚丝袜av| 巨乳人妻的诱惑在线观看| √禁漫天堂资源中文www| a级片在线免费高清观看视频| 观看美女的网站| 亚洲成色77777| 热99国产精品久久久久久7| 日韩av不卡免费在线播放| 性色av一级| 精品一区二区三区四区五区乱码 | 免费少妇av软件| 免费av中文字幕在线| 人人妻人人爽人人添夜夜欢视频| 毛片一级片免费看久久久久| 国产成人精品久久久久久| 麻豆精品久久久久久蜜桃| bbb黄色大片| 精品一区二区三卡| 国产精品女同一区二区软件| 黑丝袜美女国产一区| 一级a爱视频在线免费观看|