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

一類存在參數攝動的線性隨機系統的魯棒間歇故障診斷方法

2016-08-11 06:18:29鄢镕易何瀟周東華
自動化學報 2016年7期
關鍵詞:故障診斷故障檢測

鄢镕易  何瀟  周東華,2

一類存在參數攝動的線性隨機系統的魯棒間歇故障診斷方法

鄢镕易1何瀟1周東華1,2

間歇故障(Intermittent faults,IFs)具有隨機性,其檢測要求在本次間歇故障消失之前檢測出間歇故障的發生,在下一次間歇故障發生之前檢測出間歇故障的消失.本文針對一類存在未知時變參數攝動的離散線性隨機動態系統,研究了其魯棒間歇故障檢測與分離問題.基于降維未知輸入觀測器,通過引入滑動時間窗口,本文設計了一組與未知時變攝動解耦的結構化截斷殘差,并提出其存在的一個充分條件.與傳統殘差相比,截斷殘差信號更為顯著地反映了間歇故障的發生和消失.為滿足間歇故障的檢測要求,本文提出兩個假設檢驗分別用于檢測間歇故障的發生時刻和消失時刻,并給出了一個詳細算法.最后,在沿參考軌道運行的衛星模型上對所述方法進行了仿真實驗,結果表明該方法能夠有效檢測出間歇故障的所有發生時刻和消失時刻,并準確實現故障分離.

間歇故障,魯棒故障診斷,時變參數攝動,降維未知輸入觀測器,假設檢驗

引用格式鄢镕易,何瀟,周東華.一類存在參數攝動的線性隨機系統的魯棒間歇故障診斷方法.自動化學報,2016,42(7): 1004-1013

間歇故障(Intermittent faults,IFs)是實際工業系統中一種普遍存在的故障類型[1[5],約90%的數字電路系統崩潰由間歇故障引起;混合電路中,間歇故障發生頻率是持續故障的10~30倍[6].在航空航天系統中,元器件老化、高負荷振動、裝配不良等因素都有可能導致間歇故障NFF(No-fault-found)發生[7].在機械傳動系統中,器件磨損、載荷過重以及閥門、氣缸密閉性不良都會引起間歇故障[8-10].電力電氣系統發生間歇故障的原因包括外部環境污染、電氣接觸點腐蝕和松動等[11-13].間歇故障的發生會降低系統可靠性和安全性,增加維護維修成本[14].據統計,在軍工系統中,間歇故障導致的不必要維護維修、元部件過早更換等問題會直接引起巨額經濟損失,降低戰備完好率[6,15].隨著計算機、電子、通信等技術飛速發展,在數字化裝置廣泛普及的工業背景下,間歇故障診斷對有效避免災難性事故發生,提高系統可靠性、可維修性和保障性,降低生產成本具有十分重要的現實意義[6,15].

與持續故障不同,間歇故障具有一定隨機性,持續時間有限,故障幅值未知,無需外部補償措施失,且通常即可自行消會重復發生[16].在故障初期,間歇故障往往以類似小噪聲擾動形式出現;隨著系統運行時間增加,其持續時間和幅值逐漸增加,呈現出明顯的間歇性;在很多情況下,間歇故障能夠進一步演化為永久性故障,造成系統失效[12].考慮到間歇故障的特點,間歇故障診斷要求在每次故障消失(發生)之前檢測出間歇故障的發生(消失)[2],并能準確定位間歇故障.因此,盡管眾多研究學者對持續故障提出了很多行之有效的方法[12,16-18],卻很難滿足上述間歇故障診斷要求.

現階段,間歇故障診斷研究主要采用定性分析方法[1-2,4,12,15-16,19].文獻[12]通過構建實驗平臺模擬間歇故障對嵌入式系統的影響,結果表明間歇故障在整個使用周期中都可能出現.文獻[16]驗證了Petri網模型用于描述計算機系統接口間歇故障的有效性.現有文獻中,基于定量分析方法的間歇故障診斷理論研究成果十分有限[1-2].文獻[17]基于雙線性奉獻觀測器研究了感應電機傳感器間歇故障檢測問題,結果表明上述殘差需要充分時間衰減才能進行下一次間歇故障檢測.文獻[20]針對一類滿足Bernouli分布且均值、方差已知的執行器間歇故障,在均方穩定框架下研究了其容錯控制問題.但上述方法都只考慮間歇故障發生時刻的檢測,卻沒有檢測消失時刻,不滿足其檢測要求.文獻[2]針對一類線性連續隨機動態系統,在不考慮測量噪聲條件下,研究了間歇故障發生時刻和消失時刻的檢測問題.然而,在實際工業環境中,不僅難以獲得精確的系統解析模型,而且存在大量測量噪聲,因此,考慮存在測量噪聲條件下,研究帶有時變參數攝動的線性隨機系統的間歇故障診斷問題是十分必要的.

針對一類帶有時變參數攝動的線性離散隨機動態系統,本文研究了其魯棒間歇故障檢測與分離問題,其主要創新點包括:1)基于降維未知輸入觀測器,通過引入滑動時間窗口,設計了一組結構化截斷殘差,使其與未知時變攝動解耦且對特定方向的間歇故障敏感,以實現故障定位;2)考慮測量噪聲的影響,分析了新殘差信號的統計特性,并提出兩個假設檢驗分別用于檢測間歇故障的發生時刻和消失時刻;3)針對結構化殘差的存在性問題,本文給出了該問題可解的一個充分條件.

1 問題描述

考慮一類存在時變參數攝動的線性隨機系統

對上述系統和間歇故障,給出如下假設.

假設2.1)同一時刻僅有一個故障方向發生間歇故障;2)每一間歇故障都有已知幅值下界即每個間歇故障的持續時間/間隔時間具有最小值令假設 ττeˉ先驗τe已知. τe

2 魯棒間歇故障診斷方法

本節針對系統(1)所示一類存在未知時變參數攝動的線性離散隨機動態系統的間歇故障檢測與分離問題,提出了一種魯棒故障診斷方法.

2.1魯棒殘差設計

對系統(1)進行如下改寫

系統(4)改寫為如下l組系統模型,其中第s(s∈ lll)組系統為

定理1.對系統(3)所示一類存在未知時變參數攝動的線性離散隨機動態系統,設計l組殘差生成器(6),使其滿足條件1所示要求的一個充分條件是系統(3)滿足:具有穩定不變零點.

證明.通過線性變換,顯然,系統(4)等價于系統(3).對系統(4),為實現間歇故障分離,設計滿足式(7)所示要求的l組殘差,當滿足時,對第組系統(6),根據文獻[16],可以采用如下算式計算式(6)中參數

使得rs(k)對未知時變參數攝動和間歇故障ms(k)解耦.顯然,由于p=n-1,dim[rs(k)]=1.為簡化表示,令根據文獻[16],條件2)滿足時,殘差生成器(6)穩定,其極點能夠在單位圓內任意配置.下文中,僅以情況為例進行分析.由式(6)和式(7),可得降維估計誤差為

引入滑動時間窗口?ks,構造新的標量截斷殘差

綜上所述,若條件1)~3)滿足,基于殘差生成器(6),能夠對系統(3)設計l組滿足條件1要求的結構化魯棒殘差.

圖1 間歇故障與滑動時間窗口的相對位置關系Fig.1 Relative positions between the intermittent fault and the sliding-time window

2.2魯棒殘差統計特性分析

為簡化表示,記

2.3魯棒間歇故障診斷方法

間歇故障診斷要求在本次間歇故障消失之前確定間歇故障的發生時刻,在下一次間歇故障發生之前確定本次間歇故障的消失時刻,并準確定位故障.因此,針對間歇故障的發生和消失,本節提出兩個假設檢驗分別檢測,并給出魯棒間歇故障診斷算法.

2.3.1間歇故障發生時刻的檢測

2.3.2間歇故障消失時刻的檢測

2.3.3魯棒間歇故障分離策略

2.3.4魯棒間歇故障診斷算法

根據上述分析并參考文獻[22]給出的算法,本節給出如下魯棒間歇故障診斷算法.

步驟1.對系統(1)進行線性變換得系統(4).

步驟3.由式(4)整理得到式(5)所示l組系統;驗證每組系統是否滿足若滿足則繼續;若不滿足,對式(4)進行線性變換使式(5)滿足上述條件.

步驟7.根據式(7)計算Gs和Hs.

步驟12.根據式(13)和式(17)檢測間歇故障的發生和消失,并根據分離策略定位故障方向.

3 仿真驗證

為驗證上述方法有效性,考慮在橢圓參考軌道運行的某衛星的間歇故障診斷問題[23],其受到未知參數攝動和間歇故障影響的動力學模型為

仿真中,ai(k)服從[-1,1]的均勻分布為獨立零均值高斯白噪聲,其協方差分別為Rw=采用狀態反饋跟蹤控制律Kx=[60,7,0,8,0,0;0,-8,2,3,0,0;0,0,0,0,-14,3]使衛星沿參考軌道運行.

易知,系統(20)滿足定理1所示條件,因此,基于降維未知輸入觀測器,能夠設計3組殘差生成器,使得對未知攝動間歇故障ms(k)解耦,而對間歇故障mj(k)(j 6=s)敏感.根據上述算法,3組標量殘差生成器參數為

設計第三組殘差生成器時,可得

圖2 正常運行時的系統輸出Fig.2 Normal output of the satellite system(20)

圖3 在k=500時發生間歇m3(k)的系統輸出Fig.3 Output of system(20)subject to the IF m3(k)

當沒有故障發生時,衛星系統運動狀態如圖2所示.可以看出,由于系統動力參數存在未知時變攝動以及受到隨機噪聲影響,衛星的位移和速度出現波動.以 FF3方向發生間歇故障為例,本文給出了該衛星在間歇故障m3(k)影響下的運動狀態.如圖3所示,在k=500(即第5秒)時, FF3方向發生最小幅值ρ=0.6、最小持續/間隔時間為0.4s的間歇故障m3(k).間歇故障m3(k)的發生使系統不能按原定軌跡運行,根據圖3,無法確定m3(k)的發生時刻和消失時刻,更不能確定發生故障的執行器通道.

對上述發生間歇故障的系統(20),設計3組殘差生成器(6),利用式(13)和式(17)進行間歇故障診斷.設定滑動時間窗口選擇為顯然,滿足仿真結果如圖4~6所示.由條件1可知m3(k)敏感,與m1(k)解耦m3(k)敏感,與m2(k)解耦m2(k)敏感,與m3(k)解耦.從圖4~6可以看出,當m3(k)發生時,殘差能夠快速超過間歇故障發生時刻的檢測閾值,從而迅速確定本次間歇故障的發生時刻;間歇故障消失之后,殘差能夠在下一次間歇故障發生之前衰減到消失時刻的檢測閾值之下,從而確定本次間歇故障的消失時刻.由于與間歇故障m3(k)解耦,因此,其一直位于閾值范圍內.根據分離策略,我們可以判斷 FF3方向發生間歇故障,其發生時刻、消失時刻及實際檢測值如表1所示,檢測結果如圖7所示.可以看出,本文方法能夠迅速檢測出間歇故障的發生時刻和消失時刻并準確定位故障.

圖4 初始殘差信號r1(k)和新殘差信號r1(k,?k1)Fig.4 Comparing r1(k,?k1)with r1(k)

圖5 初始殘差信號r2(k)和新殘差信號r2(k,?k2)Fig.5 Comparing r2(k,?k2)with r2(k)

圖6 初始殘差信號r3(k)和新殘差信號r3(k,?k3)Fig.6 Comparing r3(k,?k3)with r3(k)

表1 間歇故障發生(消失)時刻及其實際檢測值Table 1 The detection result of m3(k)by using the proposed method

為了進一步說明上述方法的有效性,對于發生相同間歇故障m3(k)的系統(20),基于Kalman濾波器得到系統狀態估計值設計殘差信號為選擇構造如下的殘差評價函數其仿真結果如圖8所示,可以看出,根據此殘差值無法檢測出間歇故障m3(k)的發生時刻和消失時刻.

4 結論

本文針對一類存在未知時變參數攝動的線性離散隨機動態系統的間歇故障診斷問題,提出一種魯棒診斷方法.基于降維未知輸入觀測器,通過引入滑動時間窗口,本文設計了一組對系統未知參數攝動解耦的新的結構化標量殘差,該組殘差對間歇故障發生和消失更為敏感.基于對其統計特性的分析,根據間歇故障與滑動時間窗口相對位置關系,本文提出兩個假設檢驗用于檢測間歇故障的發生時刻和消失時刻.并利用結構化殘差集,準確實現故障定位.通過沿參考軌道運行的某衛星系統的仿真實驗,驗證了本文方法的有效性.

圖7 間歇故障檢測結果Fig.7 The detection result of m3(k)by using the proposed method

圖8 基于Kalman濾波方法的殘差信號Fig.8 The Kalman filter based residual

References

1 Zhou Dong-Hua,Shi Jian-Tao,He Xiao.Review of intermittent fault diagnosis techniques for dynamic systems.Acta Automatica Sinica,2014,40(2):161-171(周東華,史建濤,何瀟.動態系統間歇故障診斷技術綜述.自動化學報,2014,40(2):161-171)

2 Chen M Y,Xu G B,Yan R Y,Ding S X,Zhou D H.Detecting scalar intermittent faults in linear stochastic dynamic systems.International Journal of Systems Science,2015,46(8):1337-1348

3 Correcher A,Garc′?a E,Morant F,Quiles E,Blasco-Gimenez R.Intermittent failure diagnosis in industrial processes.In:Proceedings of the 2003 IEEE International Symposium on Industrial Electronics.Rio de Janeiro,Brazil:IEEE,2003. 723-728

4 Rashid L,Pattabiraman K,Gopalakrishnan S.Characterizing the impact of intermittent hardware faults on programs. IEEE Transactions on Reliability,2015,64(1):297-310

5 Shivakumar P,Kistler M,Keckler S W,Burger D,Alvisi L.Modeling the impact of device and pipeline scaling on the soft error rate of processor elements.Computer Science Department,University of Texas at Austin,2002.

6 Zhou Dong-Hua,Wei Mu-Heng,Si Xiao-Sheng.A survey on anomaly detection,life prediction and maintenance decision for industrial processes.Acta Automatica Sinica,2013,39(6):711-722(周東華,魏慕恒,司小勝.工業過程異常檢測、壽命預測與維修決策的研究進展.自動化學報,2013,39(6):711-722)

7 Sorensen B A,Kelly G,Sajecki A,Sorensen P W.An analyzer for detecting intermittent faults in electronic devices. In:Proceedings of AUTOTESTCON′94 IEEE Conference on Systems Readiness Technology— “Cost Effective Support into the Next Century”.Anaheim,USA:IEEE,1994. 417-421

8 Yesilyurt I,Gu F S,Ball A D.Gear tooth stiffness reduction measurement using modal analysis and its use in wear fault severity assessment of spur gears.NDT and E International,2003,36(5):357-372

9 Zanardelli W G,Strangas E G,Aviyente S.Identification of intermittent electrical and mechanical faults in permanent-magnet AC drives based on time-frequency analysis.IEEE Transactions on Industry Applications,2007,43(4):971-980

10 Ma Jie,Li Gang,Chen Mo.Nonlinear fault reconstruction based fault prognosis for rotating machinery.Acta Automatica Sinica,2014,40(9):2045-2049(馬潔,李剛,陳默.基于非線性故障重構的旋轉機械故障預測方法.自動化學報,2014,40(9):2045-2049)

11 Hamel A,Gaudreau A,Cote M.Intermittent arcing fault on underground low-voltage cables.IEEE Transactions on Power Delivery,2004,19(4):1862-1868

12 Correcher A,Garc′?a E,Morant F,Quiles E,Rodriguez L. Intermittent failure dynamics characterization.IEEE Transactions on Reliability,2012,61(3):649-658

13 Kim C J.Electromagnetic radiation behavior of low-voltage arcing fault.IEEE Transactions on Power Delivery,2009,24(1):416-423

14 Zhou Dong-Hua,Chen Mao-Yin,Xu Zheng-Guo.The Reliabibility Prediction and Optimal Maintenance Technology. Hefei:Press of University of Science and Technology of China,2013.(周東華,陳茂銀,徐正國.可靠性預測與最優維護技術.合肥:中國科學技術大學出版社,2013.)

15 Xu Gui-Bin.Researches on Fault Diagnosis and Prediction in Dynamic Systems[Master dissertation],Tsinghua University,China,2011.(徐貴斌.動態系統故障診斷及預測研究[碩士學位論文],清華大學,中國,2011.)

16 Krasnobaev V A,Krasnobaev L A.Application of Petri nets for the modeling of detection and location of intermittent faults in computers.Automation and Remote Control,1989,49(9):1198-1204

17 Bennett S M,Patton R J,Daley S,Newton D A.Torque and flux estimation for a rail traction system in the presence of intermittent sensor faults.In:Proceedings of United Kingdom Automatic Control Council International Conference on Control′96.Exeter University,UK:IET,1996.72-77

18 Wang Y,Xu G H,Zhang Q,Liu D,Jiang K S.Rotating speed isolation and its application to rolling element bearing fault diagnosis under large speed variation conditions. Journal of Sound and Vibration,2015,348:381-396

19 Zhou C J,Huang X F,Xiong N X,Qin Y Q,Huang S. A class of general transient faults propagation analysis for networked control systems.IEEE Transactions on Systems,Man,and Cybernetics:Systems,2015,45(4):647-661

20 Gu Zhou,Zhang Jian-Hua,Du Li-Long.Fault tolerant control for a class of time-delay systems with intermittent actuators failure.Control and Decision,2011,26(12):1829-1834(顧洲,張建華,杜黎龍.一類具有間歇性執行器故障的時滯系統的容錯控制.控制與決策,2011,26(12):1829-1834)

21 Edelmayer A,Bokor J,Szigeti F,Keviczky L.Robust detection filter design in the presence of time-varying system perturbations.Automatica,1997,33(3):471-475

22 Kudva P,Viswanadham N,Ramakrishna A.Observers for linear systems with unknown inputs.IEEE Transactions on Automatic Control,1980,25(1):113-115

23 Meskin N,Khorasani K.Fault detection and isolation of discrete-time Markovian jump linear systems with application to a network of multi-agent systems having imperfect communication channels.Automatica,2009,45(9):2032-2040

鄢镕易清華大學自動化系博士研究生.主要研究方向為間歇故障診斷與容錯控制,高速列車故障診斷與容錯控制.

E-mail:yry10@mails.tsinghua.edu.cn

(YAN Rong-YiPh.D.candidate in the Department of Automation,Tsinghua University.His research interest covers fault diagnosis and tolerance control of intermittent faults,fault diagnosis for the information control system of high-speed trains.)

何瀟清華大學自動化系副教授.主要研究方向為網絡化系統的魯棒濾波、故障診斷與容錯控制,無人機(群)智能自主控制中的安全性問題,高速列車信息控制系統的故障診斷.

E-mail:hexiao@tsinghua.edu.cn

(HE XiaoAssociate professor in the Department of Automation,Tsinghua University.His research interest covers robust estimation,fault diagnosis and tolerant control of networked systems,safety problems in intelligent autonomous control of unmanned aerial vehicles,fault diagnosis for the information control system of high-speed trains.)

周東華山東科技大學電氣與自動化工程學院教授,清華大學自動化系教授.主要研究方向為動態系統的故障診斷與容錯控制,故障預測與智能維護技術.本文通信作者.

E-mail:zdh@mail.tsinghua.edu.cn

(ZHOUDong-HuaProfessor at the College of Electrical Engineering and Automation,Shandong University of Science and Technology,and the Department of Automation,Tsinghua University.His research interest covers fault diagnosis and tolerant control,fault prediction and intelligent maintenance. Corresponding author of this paper.)

Robust Diagnosis of Intermittent Faults for Linear Stochastic Systems Subject to Time-varying Perturbations

YAN Rong-Yi1HE Xiao1ZHOU Dong-Hua1,2

Since intermittent faults(IFs)have an intermittency property,the detection of IFs requires:the current appearing time of an IF must be detected before its disappearing time;the current disappearing time of an IF must be detected before the subsequent appearing time.In this paper,the robust detection problem of IFs for a class of linear discrete-time stochastic systems subject to unknown time-varying perturbations is investigated.Based on reducedorder unknown input observers(UIOs),a novel set of structured truncated residuals is designed to detect and isolate IFs by introducing sliding-time windows,and a sufficient condition is proposed for the existence of the residual generators. Compared to traditional residuals,the novel truncated residuals,which get decoupled from time-varying perturbations,are more sensitive to the IFs.Based on the analysis of these novel residuals,two hypothesis tests are proposed to detect all the appearing times and the disappearing times of an IF.In addition,a detailed algorithm is provided to perform the given scheme.Finally,simulation results on a model of a satellite moving in a circular reference orbit are presented to illustrate the effectiveness of the proposed method.

Intermittent faults(IFs),robust fault diagnosis,unknown time-varying perturbations,reduced-order unknown input observer,hypothesis tests

10.16383/j.aas.2016.c150756

Yan Rong-Yi,He Xiao,Zhou Dong-Hua.Robust diagnosis of intermittent faults for linear stochastic systems subject to time-varying perturbations.Acta Automatica Sinica,2016,42(7):1004-1013

2015-11-11錄用日期2016-03-20
Manuscript received November 11,2015;accepted March 20,2016
國家自然科學基金(61490701,61290324,61473163,61522309),山東省泰山學者優勢特色學科人才團隊支持計劃(魯政辦字[2015]73),清華大學自主科研項目(025-陳茂銀-Z09)資助
Supported by National Natural Science Foundation of China (61490701,61290324,61473163,61522309),Research Fund for the Taishan Scholar Project of Shandong Province([2015]73),and Tsinghua University Initiative Scientific Research Program (025-CMY-Z09)
本文責任編委鐘麥英
Recommended by Associate Editor ZHONG Mai-Ying
1.清華大學自動化系北京1000842.山東科技大學電氣與自動化學院青島266590
1.Department of Automation,Tsinghua University,Beijing 1000842.College of Electrical Engineering and Automation,Shandong University of Science and Technology,Qingdao 266590

猜你喜歡
故障診斷故障檢測
“不等式”檢測題
“一元一次不等式”檢測題
“一元一次不等式組”檢測題
故障一點通
奔馳R320車ABS、ESP故障燈異常點亮
小波變換在PCB缺陷檢測中的應用
因果圖定性分析法及其在故障診斷中的應用
故障一點通
江淮車故障3例
基于LCD和排列熵的滾動軸承故障診斷
主站蜘蛛池模板: 久久久国产精品无码专区| 国产精品久久久精品三级| 国产人人干| 在线播放真实国产乱子伦| 欧美一级在线看| 亚洲国产精品日韩av专区| 国产精品永久不卡免费视频| 嫩草在线视频| 天天干伊人| 亚洲精品图区| 午夜啪啪网| 欧美专区在线观看| 国产拍在线| 看你懂的巨臀中文字幕一区二区 | 色欲综合久久中文字幕网| 中文字幕一区二区人妻电影| 国产亚洲高清视频| 欧美色视频日本| 亚洲国产综合精品一区| 国产大片喷水在线在线视频| 欧美三级日韩三级| 亚洲精品无码不卡在线播放| 亚洲成年人网| 老色鬼欧美精品| 91原创视频在线| 亚洲精品男人天堂| 亚洲精品自拍区在线观看| 一级一毛片a级毛片| 亚卅精品无码久久毛片乌克兰| 亚洲国产精品国自产拍A| 国产午夜小视频| 成人在线不卡视频| 在线观看欧美精品二区| 欧美亚洲欧美| 日韩福利视频导航| 国产美女免费| 欧洲熟妇精品视频| 久久人与动人物A级毛片| 亚洲精品欧美日本中文字幕| 欧美α片免费观看| 91年精品国产福利线观看久久 | 97青青青国产在线播放| 久久久精品国产亚洲AV日韩| 欧美一区二区福利视频| 野花国产精品入口| 欧美a级在线| 色噜噜在线观看| 亚洲欧美人成电影在线观看| 亚洲人成影院在线观看| 亚洲欧美另类久久久精品播放的| 国产福利免费在线观看| 亚洲AV无码乱码在线观看裸奔| 91成人免费观看| 精品成人免费自拍视频| 人人爱天天做夜夜爽| 精品视频第一页| 91激情视频| 一级毛片不卡片免费观看| AⅤ色综合久久天堂AV色综合 | 欧美亚洲欧美| 国产精品男人的天堂| 丰满少妇αⅴ无码区| 九九久久99精品| 色丁丁毛片在线观看| 亚洲天堂伊人| 欧美成人午夜视频免看| 国产精品999在线| 久久免费观看视频| 啪啪免费视频一区二区| 国产玖玖视频| 亚洲AV无码一二区三区在线播放| 日本高清在线看免费观看| 亚洲综合九九| av一区二区无码在线| 欧美另类图片视频无弹跳第一页| 欧美一区国产| 蜜桃视频一区二区| 国产女人18水真多毛片18精品| 国产毛片高清一级国语 | 中文字幕在线看| 国产a v无码专区亚洲av| 国产国模一区二区三区四区|