• <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)鍵詞:故障診斷故障方法

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

    (中國計量大學機電工程學院,浙江杭州 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ù)學模型的控制對象,例如參數(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]采用學習矢量量化(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ù)學模型.

    本文的主要貢獻:一是提出了一種適用于定值和隨動單回路閉環(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一区二区三区| 亚洲av成人不卡在线观看播放网| 男人舔女人的私密视频| 99久久精品热视频| 国产私拍福利视频在线观看| 91久久精品国产一区二区成人 | 欧美日韩一级在线毛片| 精品久久久久久成人av| 国产伦在线观看视频一区| 免费av不卡在线播放| 黑人巨大精品欧美一区二区mp4| av中文乱码字幕在线| 他把我摸到了高潮在线观看| 日韩欧美 国产精品| 中出人妻视频一区二区| 亚洲第一电影网av| 国产精品 国内视频| 看黄色毛片网站| 亚洲五月婷婷丁香| 人妻夜夜爽99麻豆av| 亚洲成人久久爱视频| 97人妻精品一区二区三区麻豆| 成人亚洲精品av一区二区| 国产精品九九99| 黄色日韩在线| 熟女少妇亚洲综合色aaa.| 久久中文字幕一级| 久久欧美精品欧美久久欧美| 国产成人av激情在线播放| 色精品久久人妻99蜜桃| 国产欧美日韩一区二区三| 欧美三级亚洲精品| 怎么达到女性高潮| 1024手机看黄色片| 国产精品久久久人人做人人爽| 日韩欧美精品v在线| 免费高清视频大片| 母亲3免费完整高清在线观看| 免费在线观看成人毛片| 小蜜桃在线观看免费完整版高清| 国产乱人视频| 亚洲国产欧洲综合997久久,| 久久久久国内视频| 天堂影院成人在线观看| 又黄又爽又免费观看的视频| 搡老岳熟女国产| 色综合站精品国产| 亚洲人与动物交配视频| 天堂av国产一区二区熟女人妻| 国产亚洲av高清不卡| 亚洲国产精品sss在线观看| 女生性感内裤真人,穿戴方法视频| 欧美成人一区二区免费高清观看 | ponron亚洲| 欧美午夜高清在线| 18美女黄网站色大片免费观看| 狂野欧美激情性xxxx| 99国产精品99久久久久| 日韩精品青青久久久久久| 国产成年人精品一区二区| 日本 欧美在线| 亚洲第一欧美日韩一区二区三区| 少妇熟女aⅴ在线视频| 婷婷精品国产亚洲av在线| 97超视频在线观看视频| 亚洲黑人精品在线| 亚洲国产中文字幕在线视频| av天堂中文字幕网| 日本 av在线| 51午夜福利影视在线观看| 免费在线观看成人毛片| 日韩欧美在线乱码| 亚洲乱码一区二区免费版| 国产三级黄色录像| 午夜激情福利司机影院| 国产精品电影一区二区三区| 三级国产精品欧美在线观看 | 国产欧美日韩精品一区二区| 亚洲成人免费电影在线观看| 日韩欧美三级三区| 夜夜爽天天搞| 高清在线国产一区| 国语自产精品视频在线第100页| 亚洲美女黄片视频| 国产精品九九99| 麻豆国产av国片精品| 精品久久久久久久久久久久久| 少妇熟女aⅴ在线视频| 村上凉子中文字幕在线| 免费观看人在逋| 久久这里只有精品19| 国产aⅴ精品一区二区三区波| 亚洲欧美日韩东京热| 亚洲av成人不卡在线观看播放网| 欧美性猛交╳xxx乱大交人| 熟女人妻精品中文字幕| 99精品在免费线老司机午夜| 999久久久国产精品视频| 久久精品影院6| 国产亚洲av嫩草精品影院| 欧美大码av| 三级男女做爰猛烈吃奶摸视频| 午夜福利欧美成人| 日本 欧美在线| www.精华液| 免费在线观看视频国产中文字幕亚洲| 国产野战对白在线观看| 51午夜福利影视在线观看| 欧美色欧美亚洲另类二区| 国产激情偷乱视频一区二区| 此物有八面人人有两片| 成人特级av手机在线观看| 久久午夜亚洲精品久久| 精品一区二区三区av网在线观看| 亚洲专区中文字幕在线| 俄罗斯特黄特色一大片| 两人在一起打扑克的视频| 国产激情欧美一区二区| 又粗又爽又猛毛片免费看| 国产午夜精品久久久久久| 免费观看的影片在线观看| xxx96com| 嫩草影视91久久| 亚洲18禁久久av| 久久九九热精品免费| 观看免费一级毛片| 精品一区二区三区视频在线观看免费| 国产精品免费一区二区三区在线| 亚洲国产色片| 好看av亚洲va欧美ⅴa在| 精品人妻1区二区| 成人三级做爰电影| 国产黄片美女视频| 美女高潮喷水抽搐中文字幕| 亚洲人成电影免费在线| 搡老熟女国产l中国老女人| 国产精品久久久久久亚洲av鲁大| 日本a在线网址| 亚洲欧美激情综合另类| 老司机深夜福利视频在线观看| 国产三级在线视频| 欧美黑人巨大hd| 蜜桃久久精品国产亚洲av| 国产精品一区二区三区四区免费观看 | 18禁裸乳无遮挡免费网站照片| 少妇的丰满在线观看| 亚洲国产日韩欧美精品在线观看 | 国产私拍福利视频在线观看| 久9热在线精品视频| 亚洲欧美日韩高清在线视频| av天堂在线播放| 听说在线观看完整版免费高清| 国产日本99.免费观看| 欧美日本视频| 国产主播在线观看一区二区| 高清毛片免费观看视频网站| www.精华液| 国产亚洲av高清不卡| 亚洲av美国av| 久久这里只有精品19| 香蕉久久夜色| 婷婷六月久久综合丁香| 午夜福利在线观看吧| 日韩欧美国产在线观看| 91字幕亚洲| 听说在线观看完整版免费高清| 欧美黄色淫秽网站| 亚洲熟女毛片儿| 免费电影在线观看免费观看| 黄色片一级片一级黄色片| 91在线观看av| 久久这里只有精品中国| 99在线视频只有这里精品首页| 脱女人内裤的视频| 无遮挡黄片免费观看| 18禁美女被吸乳视频| 国产激情久久老熟女| 在线a可以看的网站| 色吧在线观看| 在线看三级毛片| 欧美乱色亚洲激情| 精品乱码久久久久久99久播| 国产日本99.免费观看| 国产精品乱码一区二三区的特点| 久久中文看片网| av女优亚洲男人天堂 | 岛国在线免费视频观看| 亚洲avbb在线观看| 三级国产精品欧美在线观看 | 老司机深夜福利视频在线观看| 一级a爱片免费观看的视频| cao死你这个sao货| 欧美国产日韩亚洲一区| 一进一出抽搐gif免费好疼| 久久精品综合一区二区三区| 色综合欧美亚洲国产小说| 亚洲av五月六月丁香网| 亚洲av成人av| 色综合欧美亚洲国产小说| 国产精品99久久99久久久不卡| 国产私拍福利视频在线观看| x7x7x7水蜜桃| 亚洲av五月六月丁香网| 国产精品九九99| av片东京热男人的天堂| 黑人巨大精品欧美一区二区mp4| 亚洲精品色激情综合| 成人性生交大片免费视频hd| 中文资源天堂在线| 国产精品久久久人人做人人爽| 一本一本综合久久| 日韩欧美在线乱码| 国产三级中文精品| 男女之事视频高清在线观看| 级片在线观看| 免费看a级黄色片| 亚洲黑人精品在线| 国产激情欧美一区二区| 女人高潮潮喷娇喘18禁视频| 黄片大片在线免费观看| 看黄色毛片网站| 最新中文字幕久久久久 | 一个人看的www免费观看视频| 99精品欧美一区二区三区四区| 美女 人体艺术 gogo| 亚洲va日本ⅴa欧美va伊人久久| 国产免费av片在线观看野外av| 亚洲五月婷婷丁香| 亚洲午夜理论影院| 18禁国产床啪视频网站| 成人鲁丝片一二三区免费| 日本成人三级电影网站| a级毛片a级免费在线| 一卡2卡三卡四卡精品乱码亚洲| 全区人妻精品视频| 18禁美女被吸乳视频| xxx96com| 亚洲国产看品久久| av女优亚洲男人天堂 | 国产亚洲精品久久久com| 欧美日本视频| 中文字幕人妻丝袜一区二区| 欧美3d第一页| 夜夜看夜夜爽夜夜摸| 亚洲性夜色夜夜综合| 日本三级黄在线观看| 国内少妇人妻偷人精品xxx网站 | 麻豆久久精品国产亚洲av| 香蕉丝袜av| 国产精品一区二区免费欧美| 亚洲国产精品sss在线观看| 国产精华一区二区三区| 国产成+人综合+亚洲专区| 国产免费av片在线观看野外av| 国产精华一区二区三区| 国产精品一区二区三区四区免费观看 | 美女黄网站色视频| 老司机福利观看| 两个人的视频大全免费| а√天堂www在线а√下载| 99久久成人亚洲精品观看| 身体一侧抽搐| 国产精品 国内视频| 窝窝影院91人妻| 美女免费视频网站| 国产精品久久视频播放| 99久久成人亚洲精品观看| 999久久久精品免费观看国产| 欧美日韩综合久久久久久 | 美女扒开内裤让男人捅视频| 99久久成人亚洲精品观看| 国产 一区 欧美 日韩| 一卡2卡三卡四卡精品乱码亚洲| 法律面前人人平等表现在哪些方面| 国产视频内射| 亚洲第一欧美日韩一区二区三区| 欧美一级毛片孕妇| 国产精品 欧美亚洲| 天堂网av新在线| 精品久久久久久久人妻蜜臀av| 99久久99久久久精品蜜桃| 亚洲成av人片在线播放无| 国产免费男女视频| 欧美乱色亚洲激情| 国内精品一区二区在线观看| 亚洲电影在线观看av| 香蕉久久夜色| 久久久久国产精品人妻aⅴ院| 国产三级在线视频| 精品国产乱码久久久久久男人| 很黄的视频免费| 99热这里只有是精品50| 麻豆久久精品国产亚洲av| 久久久久久久久久黄片| 国产麻豆成人av免费视频| 老熟妇乱子伦视频在线观看| 国产精品自产拍在线观看55亚洲| 一级毛片女人18水好多| 国产主播在线观看一区二区| 国产精品亚洲av一区麻豆| 999久久久精品免费观看国产| 亚洲精品中文字幕一二三四区| 99热这里只有是精品50| 草草在线视频免费看| 三级毛片av免费| 女同久久另类99精品国产91| 亚洲欧美一区二区三区黑人| 久久久久久久久中文| 久久精品国产清高在天天线| 91麻豆av在线| 亚洲精品一区av在线观看| 久久性视频一级片| 黄色日韩在线| 成人午夜高清在线视频| 色综合婷婷激情| 亚洲,欧美精品.| 亚洲欧美精品综合久久99| 国产爱豆传媒在线观看| 这个男人来自地球电影免费观看| 757午夜福利合集在线观看| 久久香蕉国产精品| 久久热在线av| 黄色女人牲交| 亚洲欧美精品综合久久99| 亚洲专区中文字幕在线| 久久久色成人| 久久香蕉精品热| 91av网站免费观看| 亚洲 欧美一区二区三区| 久久国产乱子伦精品免费另类| 欧美乱妇无乱码| 欧美精品啪啪一区二区三区| 国产午夜精品久久久久久| 国产高清videossex| 国产精品久久久av美女十八| 午夜激情欧美在线| 黄色 视频免费看| 最近最新免费中文字幕在线| 岛国在线观看网站| 国产伦人伦偷精品视频| 中文字幕av在线有码专区| 99热只有精品国产| 在线播放国产精品三级| 一a级毛片在线观看| 欧美色欧美亚洲另类二区| 身体一侧抽搐| 黄片大片在线免费观看| 又紧又爽又黄一区二区| 母亲3免费完整高清在线观看| 麻豆一二三区av精品| 久久精品国产清高在天天线| 国产97色在线日韩免费| 国产精品一区二区三区四区免费观看 | 18美女黄网站色大片免费观看| 亚洲av电影不卡..在线观看| 曰老女人黄片| 国产高清激情床上av| 午夜福利在线观看吧| svipshipincom国产片| 熟女电影av网| 国产精品,欧美在线| 99在线视频只有这里精品首页| 日日夜夜操网爽| 国产乱人视频| 性色av乱码一区二区三区2| 美女扒开内裤让男人捅视频| 又紧又爽又黄一区二区| 床上黄色一级片| 99久久精品国产亚洲精品| 国产亚洲精品一区二区www| 搡老岳熟女国产| 夜夜爽天天搞| 国产成人精品久久二区二区免费| 日韩三级视频一区二区三区| 一级毛片女人18水好多| 国产精品一区二区三区四区久久| 国产精品一区二区免费欧美| 丰满人妻一区二区三区视频av | 亚洲av成人精品一区久久| 99久久精品热视频| 欧美一级毛片孕妇| 91av网一区二区| 日本撒尿小便嘘嘘汇集6| 亚洲中文字幕一区二区三区有码在线看 | 淫妇啪啪啪对白视频| 国产免费av片在线观看野外av| 精品久久蜜臀av无| 国产精品久久久av美女十八| 女人被狂操c到高潮| 国产三级中文精品| 无人区码免费观看不卡| 一进一出抽搐gif免费好疼| 精品久久久久久久毛片微露脸| 欧美成人免费av一区二区三区| 亚洲专区国产一区二区| 欧美性猛交黑人性爽| 91麻豆av在线| 精品不卡国产一区二区三区| 亚洲精品在线观看二区| 欧美国产日韩亚洲一区| 欧美国产日韩亚洲一区| 中文字幕av在线有码专区| 成人特级黄色片久久久久久久| 久久99热这里只有精品18| 国内精品一区二区在线观看| 国产精品久久久久久亚洲av鲁大| 性色av乱码一区二区三区2| 一本精品99久久精品77| 亚洲专区中文字幕在线| 黄频高清免费视频| 国产成人啪精品午夜网站| 午夜视频精品福利| 亚洲av熟女| 亚洲在线自拍视频| 亚洲真实伦在线观看| 一个人观看的视频www高清免费观看 | 亚洲五月婷婷丁香| 久久精品国产综合久久久| 国产91精品成人一区二区三区| 国产成人一区二区三区免费视频网站| 国产乱人视频| 免费在线观看成人毛片| 神马国产精品三级电影在线观看| 欧美色欧美亚洲另类二区| 人人妻,人人澡人人爽秒播| 无人区码免费观看不卡| av天堂在线播放| 老司机午夜十八禁免费视频| av天堂在线播放| 97超级碰碰碰精品色视频在线观看| 男人舔奶头视频| 亚洲人成伊人成综合网2020| 午夜影院日韩av| 久久这里只有精品19| 亚洲av成人不卡在线观看播放网| 国产亚洲欧美98| 日本免费一区二区三区高清不卡| 久久精品人妻少妇| 天天添夜夜摸| 成人永久免费在线观看视频| 亚洲狠狠婷婷综合久久图片| 91九色精品人成在线观看| 久久中文字幕一级| 久久国产精品人妻蜜桃| 啦啦啦免费观看视频1| 免费在线观看成人毛片| 成人av一区二区三区在线看| 黄片大片在线免费观看| 国产精品国产高清国产av| 俄罗斯特黄特色一大片| 成人一区二区视频在线观看| 18禁国产床啪视频网站| 精品久久久久久成人av| 亚洲成人精品中文字幕电影| 在线观看免费午夜福利视频| 一区二区三区激情视频| 欧美日韩乱码在线| 丁香六月欧美| 国产免费av片在线观看野外av| 99久久精品国产亚洲精品| 国产爱豆传媒在线观看| 麻豆国产av国片精品| 亚洲人成网站在线播放欧美日韩| 欧美日韩乱码在线| 看片在线看免费视频| 香蕉av资源在线| 国产成人啪精品午夜网站| 久久久久久九九精品二区国产| 欧美黑人巨大hd| 搞女人的毛片| 国产成人福利小说| 精品99又大又爽又粗少妇毛片 | 长腿黑丝高跟| 熟女少妇亚洲综合色aaa.| 男人舔奶头视频| 久久久色成人| 亚洲欧美精品综合一区二区三区| 极品教师在线免费播放| 脱女人内裤的视频| 他把我摸到了高潮在线观看| 久久草成人影院| 人人妻,人人澡人人爽秒播| 亚洲成av人片免费观看| 色老头精品视频在线观看| 久久久精品欧美日韩精品| 看免费av毛片| 免费高清视频大片| 少妇的逼水好多| 我要搜黄色片| 网址你懂的国产日韩在线| 日韩高清综合在线| 国产高潮美女av| tocl精华| www.自偷自拍.com| 可以在线观看的亚洲视频| 亚洲电影在线观看av| 在线观看美女被高潮喷水网站 | 午夜福利在线在线| 午夜福利在线观看吧| 搞女人的毛片| 色综合站精品国产| 制服人妻中文乱码| 日韩成人在线观看一区二区三区| 精品一区二区三区视频在线 | 久久精品人妻少妇| 香蕉国产在线看| 亚洲av五月六月丁香网| 日韩欧美免费精品| 99国产精品99久久久久| 亚洲精品美女久久久久99蜜臀| 中文字幕最新亚洲高清| 国产成人系列免费观看| 中文字幕精品亚洲无线码一区| 亚洲中文字幕一区二区三区有码在线看 | 最近在线观看免费完整版| 在线观看免费视频日本深夜| 成人午夜高清在线视频| 国产精品99久久久久久久久| 亚洲中文av在线| 日韩大尺度精品在线看网址| 99视频精品全部免费 在线 | 日韩av在线大香蕉| xxx96com| 日韩欧美三级三区| 国产aⅴ精品一区二区三区波| 婷婷六月久久综合丁香| 亚洲欧美精品综合久久99| 婷婷精品国产亚洲av| 国产免费av片在线观看野外av| 欧美丝袜亚洲另类 | 婷婷精品国产亚洲av在线| bbb黄色大片| 亚洲色图av天堂| 国产黄色小视频在线观看| 制服人妻中文乱码| 国产一区二区在线观看日韩 | 女警被强在线播放| 日韩中文字幕欧美一区二区| 亚洲精品久久国产高清桃花| 老司机在亚洲福利影院| 噜噜噜噜噜久久久久久91| 91在线观看av| 国产精品久久久人人做人人爽| 男女视频在线观看网站免费| 日韩国内少妇激情av| 琪琪午夜伦伦电影理论片6080| 欧美3d第一页| 伊人久久大香线蕉亚洲五| 亚洲五月天丁香| 国产成人精品久久二区二区91| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品999在线| 男插女下体视频免费在线播放| 免费看光身美女| 91av网站免费观看| 韩国av一区二区三区四区| 麻豆国产97在线/欧美| 国产激情久久老熟女| 国产乱人视频| 男插女下体视频免费在线播放| av天堂中文字幕网| 在线观看美女被高潮喷水网站 | 久久久国产欧美日韩av| 午夜激情福利司机影院| 天天一区二区日本电影三级| 一个人免费在线观看电影 | 黄色视频,在线免费观看| 一级毛片高清免费大全| 琪琪午夜伦伦电影理论片6080| 国产激情欧美一区二区| 亚洲精品粉嫩美女一区| 久久人妻av系列| 丰满人妻一区二区三区视频av | 国产成人av教育| 亚洲专区中文字幕在线| 亚洲欧美激情综合另类| 色在线成人网| 无人区码免费观看不卡| 99久久久亚洲精品蜜臀av| 欧美乱妇无乱码| 夜夜看夜夜爽夜夜摸| 男插女下体视频免费在线播放| 日本成人三级电影网站| 欧美xxxx黑人xx丫x性爽| 在线a可以看的网站| a级毛片在线看网站| 欧美日韩综合久久久久久 | 99精品久久久久人妻精品| 午夜视频精品福利| 国产成人aa在线观看| 大型黄色视频在线免费观看| 午夜激情欧美在线| av天堂在线播放| 性欧美人与动物交配| 午夜福利欧美成人| 99热这里只有精品一区 | 亚洲五月天丁香| 日本熟妇午夜| 夜夜爽天天搞| 99视频精品全部免费 在线 | www日本在线高清视频| 狠狠狠狠99中文字幕| 三级国产精品欧美在线观看 | 欧美黑人欧美精品刺激| 国产成+人综合+亚洲专区| 久久九九热精品免费| 午夜影院日韩av| 91老司机精品| 黄频高清免费视频| 麻豆久久精品国产亚洲av| 欧美成人免费av一区二区三区| 99在线人妻在线中文字幕| 9191精品国产免费久久| 午夜福利在线观看吧| 午夜免费成人在线视频|