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

    1種未知輸入觀測器的推廣設計方法

    2017-11-09 09:07:46龍一夫
    航空發(fā)動機 2017年5期
    關鍵詞:故障診斷發(fā)動機測量

    龍一夫,王 曦

    (1.北京航空航天大學能源與動力工程學院,2.先進航空發(fā)動機協(xié)同創(chuàng)新中心:北京100191)

    1種未知輸入觀測器的推廣設計方法

    龍一夫1,2,王 曦1,2

    (1.北京航空航天大學能源與動力工程學院,2.先進航空發(fā)動機協(xié)同創(chuàng)新中心:北京100191)

    為使未知輸入觀測器的應用范圍更廣,提出1種解決在未知輸入觀測器(UIO)設計過程中對象測量輸出線性相關問題的方法。通過對測量輸出的線性相關性分析,測量輸出可以被合理分組,使得各組內包含部分互相線性無關的測量輸出并且保證這些分組涵蓋了所有需要診斷的測量輸出傳感器;以這種方式來滿足設計UIO的存在條件,并在此基礎上針對不同分組來設計UIO以進行關于傳感器的故障診斷算法設計。針對1個工程實例,使用UIO進行傳感器的故障診斷設計驗證了該方法的可行性。

    未知輸入觀測器;故障診斷;線性相關;傳感器

    0 引言

    現代控制系統(tǒng)越來越復雜,對于其可靠性和安全性的要求也就越來越被重視;如果能在破壞性的故障發(fā)生前發(fā)現并采取措施,就可以避免災難性的損失。因而,控制系統(tǒng)的故障診斷十分必要[1-2]。

    被監(jiān)控系統(tǒng)的數學模型是基于模型故障診斷技術的基礎;因此,一般模型對真實系統(tǒng)的描述越精確,故障診斷的能力和性能也就越好。但是對于航空發(fā)動機而言,建模的誤差與干擾是不可避免的,因而故障診斷系統(tǒng)的魯棒性是非常重要的,即處于模型與真實系統(tǒng)之間存在差別或受到干擾的情況下,故障診斷系統(tǒng)生成的殘差仍然能夠保持對故障敏感而對那些干擾不敏感[3-5]。

    本文介紹的基于未知輸入觀測器(UIO)的魯棒傳感器故障檢測方法,是1種基于模型的FDI方法[6-7],與傳統(tǒng)的卡爾曼濾波器不同[8-10],由于該方法能實現在保證其生成的殘差對故障敏感的同時保持對未知干擾不敏感這一特點,可以很好地用于航空發(fā)動機傳感器的故障診斷中。

    1 UIO的基本原理

    UIO本質上是針對不確定系統(tǒng)設計的系統(tǒng)觀測器,其中系統(tǒng)的不確定性可以歸納為動態(tài)方程中的加性干擾項,這類不確定的系統(tǒng)可以描述為

    可以發(fā)現,這種描述形式下的系統(tǒng)與嚴格正則的線性系統(tǒng)相比,在系統(tǒng)的狀態(tài)方程中多出了1項未知輸入項Ed(t),其中E為未知輸入分布矩陣,d(t)為未知干擾向量[1,7]。

    UIO定義為針對式(1)所描述的系統(tǒng),如果其觀測器的狀態(tài)估計誤差無論是否存在未知輸入都漸近地趨于零,那么稱該觀測器為未知輸入觀測器。

    本文中討論的UIO為全階觀測器且具有如下的常規(guī)結構

    也就是說如果知道式(1)中的所有參數矩陣,就可以設計出 1 個如式(2)所示的 UIO;其中 F、T、K、H為設計參數矩陣,用以滿足未知輸入干擾和狀態(tài)估計誤差解耦以及其他設計要求。

    式(2)所示的系統(tǒng)為式(1)所示系統(tǒng)的UIO的充分必要條件[1]為:rank(CE)=rank(E),(C,A1)矩陣對可檢測。其中

    在滿足以上2個條件的情況下,可以計算出相應的UIO設計參數矩陣,限于篇幅,具體過程這里未給出,可以參考文獻[1]。

    由以上內容可知,UIO設計的1個關鍵之處在于未知輸入干擾項Ed(t)的獲取。對于航空發(fā)動機而言,未知輸入干擾很多,包括建模誤差、線性化誤差、飛行條件變化干擾等[11-15]。一般很難直接獲得其未知輸入干擾信息,故針對式(1)的系統(tǒng)可以假定未知干擾向量d(t)是緩慢時變矢量,對于航空發(fā)動機來說這個假定是合理的;那么系統(tǒng)可以改寫為

    其中d1(t)=Ed(t)。那么如果知道系統(tǒng)的輸入輸出,就可以針對式(4)所示的系統(tǒng)設計1個狀態(tài)觀測器來獲得d1(t)的估計d^1(t),以此獲得有關分布矩陣E的一些信息。

    式(4)系統(tǒng)中的能觀矩陣為

    等式右邊第2個矩陣為滿秩陣,故可以知道

    當且僅當rank(W0)=2n時,式(4)所描述的系統(tǒng)是能觀的,這就要求。從這里可知要求系統(tǒng)擁有n個(狀態(tài)維數)互相線性無關的測量輸出,并且(C,A)矩陣對是可觀測的。(C,A)可觀測這個條件對于一般的航空發(fā)動機來說都是可以滿足的;所以主要限制在于要求系統(tǒng)擁有n個互相線性無關的測量輸出[1,7]。

    2 UIO的推廣方法及實例

    本文主要提出了1種當系統(tǒng)的獨立測量輸出小于系統(tǒng)的狀態(tài)維數時也能觀測未知干擾項的方法。

    系統(tǒng)的獨立測量輸出個數小于系統(tǒng)的狀態(tài)維數,一般有以下2種情況,假設A矩陣為n行n列(n為系統(tǒng)狀態(tài)的維數),C矩陣為m行n列(m即為發(fā)動機測量輸出的個數):

    (1)rank(C)=m,m<n

    (2)rank(C)=p<m,m<n

    針對第1種情況下的系統(tǒng),觀測器存在的最好辦法是將發(fā)動機的線性模型降階使其滿足獨立的測量量個數不小于系統(tǒng)狀態(tài)的維數;但是如果系統(tǒng)的階次不能再降低這種情況出現時,則可以考慮引入幾個與現有測量量線性無關的虛擬輸出量,使得系統(tǒng)滿足rank(C)=n。例如,從發(fā)動機的狀態(tài)量中取與測量輸出線性無關的幾個量使得增廣后的虛擬系統(tǒng)滿足rank(C)=n。值得一提的是,虛擬輸出量的引入僅僅起到設計增廣系統(tǒng)觀測器的橋梁作用,并不會增加實際發(fā)動機輸出,而只是在設計階段引入一些輔助量,這些量的獲取可以從發(fā)動機模型的計算結果中直接得到。

    如果是第2種情況下的系統(tǒng)則稍微復雜,因為此時系統(tǒng)的測量量之間不是互相獨立的,需要對系統(tǒng)做些處理。一般方法是對這些測量量進行分組,對于航空發(fā)動機來說,互相線性相關的測量輸出不會太多,一般分成2組就基本上可保證每組中的測量量互相線性無關;同時使這2組測量量包含發(fā)動機所有的測量輸出。

    下面以某型渦扇航空發(fā)動機的全階線性模型為例,說明如何解決在其UIO設計時出現的測量輸出線性相關問題。

    發(fā)動機線性模型如下

    根據上述發(fā)動機線性模型來設計發(fā)動機的未知輸入觀測器。其中發(fā)動機4個輸入為TBV(高壓壓氣機放氣活門角度)、VBV(增壓級可調放氣活門面積)、VSV(高壓壓氣機可調靜子導葉角)和Wf(燃油流量),9個測 量 輸 出 分 別 為 :T3、P25、P3、P13、T42、N2、N1、T5、T25;12個 狀 態(tài) 分 別 為 T4、P4、T48、P48、T17、P17、T25、P25、T7、P7、N2和N1;其中N2和N1表示高低壓轉子轉速,3截面為高壓壓氣機出口截面,13截面為風扇外涵出口截面,17截面為外涵噴管出口截面,25截面為高壓壓氣機進口截面,42截面為高壓渦輪出口截面,48截面為低壓渦輪進口截面,5截面為低壓渦輪出口截面,7截面為內涵噴管出口截面。

    發(fā)動機的輸入數由控制調節(jié)計劃決定,共有4個控制輸入。輸出傳感器的選取包括控制計劃的需要以及狀態(tài)監(jiān)視的要求。而發(fā)動機的狀態(tài)數與發(fā)動機的階數相對應,即是模型中考慮到的微分方程的個數。具體控制計劃內容不在本文研究范圍內。

    按照式(4)對發(fā)動機線性模型增廣得到線性模型A、B、C、D;為了獲取系統(tǒng)的未知干擾信息,需要對此增廣系統(tǒng)設計觀測器。如前所述,此系統(tǒng)能觀的條件是rank(C)=n;但實際從上述發(fā)動機線性模型中可以得到rank(C)=8,而發(fā)動機的狀態(tài)共有12個,測量輸出共有9個,而C矩陣的每行對應1個發(fā)動機的測量輸出,說明發(fā)動機的測量輸出中并不是相互線性無關的,且符合前文所說的第2種情況。于是可以將這些輸出分成相互線性無關的2組:第1組包含按照順序的后8個測量輸出,第2組包含前8個測量輸出。再分別從發(fā)動機的狀態(tài)中找出T17、T48、T7和T4這4個與輸出線性無關的量作為輔助。這樣C矩陣就被增廣成2個12×12并且滿秩的方陣,滿足觀測器的存在條件。由于實際發(fā)動機中并不包含這4個輔助信號的傳感器,故在實際使用中這4個信號的值可以由發(fā)動機模型計算得到。

    由此設計出來的觀測器就可以觀測未知輸入干擾項d^1(t)用于設計系統(tǒng)的UIO。

    在發(fā)動機飛行條件以及油門桿指令變化情況如圖1所示。利用未知輸入干擾項估計出來的幾個重要傳感器輸出值與發(fā)動機非線性模型傳感器輸出值的對比分別如圖2~6所示。從圖中可見,在不同工況下發(fā)動機模型輸出與觀測器的估計輸出都十分接近,可見觀測器對于未知輸入干擾的觀測效果較好。

    利用UIO實現的發(fā)動機控制系統(tǒng)傳感器的故障診斷系統(tǒng)包括最優(yōu)未知輸入觀測器組、故障檢測模塊以及故障隔離和重構模塊,如圖7所示。

    對于單回路傳感器故障,首先,根據發(fā)動機標稱穩(wěn)態(tài)點的線性狀態(tài)空間模型設計出針對每個傳感器的最優(yōu)未知輸入觀測器組,將控制系統(tǒng)給出的發(fā)動機控制量和傳感器測量值作為未知輸入觀測器組的輸入,每個未知輸入觀測器都已除去1路傳感器外其它的傳感器測量值作為輸入。進而,由觀測器組獲得相應傳感器狀態(tài)估計的殘差加權平方和(Weighted Sum-Squared Residual,WSSR)作為故障指示信號并分析,由故障檢測模塊得到故障信息[11-13]。

    帶有未知干擾的發(fā)動機線性系統(tǒng)為

    針對第i個傳感器設計的最優(yōu)UIO為

    式中:x^i(t)為針對第i個傳感器設計的觀測器的狀態(tài)量最優(yōu)估計值;yik為測量值yk除第i行值的子集;y^i為yi的估計值;Ci為C除去第i行的子集。不失一般性,在這里假設系統(tǒng)均為嚴格正則系統(tǒng),E為系統(tǒng)對應的干擾分布矩陣。

    由于本文的測量輸出被分成2組,所以對應設計的UIO也同樣有2組,各8個。將發(fā)動機的油門桿輸入信號設置(如圖8所示),仍然保持在海平面標準狀態(tài)下的高度0 km和馬赫數0。以P3(高壓壓氣機出口壓力)傳感器在35 s時刻發(fā)生硬故障為例給出一些診斷結果。P3傳感器輸出隨時間的響應如圖9所示。

    圖9 中在15 s時刻P3與油門桿輸入信號的階躍相對應,進入穩(wěn)態(tài)之后在35 s時刻斷路故障被注入。

    去除T3傳感器信號后設計出的第1組各UIO的殘差響應如圖10所示。

    除去T25傳感器信號后構成的8個相互獨立的輸出量各自的最優(yōu)UIO殘差響應如圖11所示。

    從圖11中可見,當Ps3傳感器故障時,只有針對Ps3傳感器設計的最優(yōu)UIO的估計殘差保持在很小的范圍內不變,而其他UIO的估計殘差都會在Ps3傳感器發(fā)生故障時瞬時變大并快速保持為較大的穩(wěn)態(tài)值,正是通過這種特性,可以將故障進行檢測與隔離;值得一提的是,在15 s時刻的發(fā)動機油門桿指令階躍引起的殘差波動相對于故障時引起的殘差變化來說是很小的,所以可以通過合理的設置故障診斷的閾值來避免發(fā)動機工況的變化被誤認為是故障。

    同樣還以圖8作為發(fā)動機的油門桿指令輸入、圖9所示的P3傳感器發(fā)生斷路故障為例,給出UIO設計時不分組各UIO的殘差響應情況作為對比,如圖12所示。

    從圖12中可見,在不進行傳感器分組時,同樣飛行條件和故障注入下,針對T5以及T25傳感器設計的UIO殘差相性相比分組時的殘差響應要差,并且不分組時針對T25傳感器設計的UIO殘差響應在沒有注入故障時就已經很大,而這時其他UIO估計殘差均十分小,這樣會引起診斷系統(tǒng)的誤報警。

    3 總結

    從上述結果可知,可以通過分組的方式來避免由于發(fā)動機測量輸出線性相關導致的UIO無法設計的情況。此外,通過與不分組情況的UIO殘差進行的1組對比可知,分組能夠更加有效地實現故障診斷,降低誤報警率;通過這種方法分組后設計的UIO能實現對故障敏感且同時保證了對發(fā)動機小工況范圍內變化的不敏感,能夠很好地實現故障的診斷及隔離,達到診斷的目的。同時也推廣了UIO方法在故障診斷中的應用范圍,不再局限于要求應用該方法時需要保證線性無關的測量輸出數不小于系統(tǒng)狀態(tài)的維數,有一定的工程應用價值。

    [1]Chen Jie,Patton R J.Robust model-based fault diagnosis for dynamic systems[M].Netherlands:Kluwer Academic Publishers,1999:17-52.

    [2]Chen C T.Linear system theory and design[M].New York Oxford:Oxford university press,1999:12-31.

    [3]Ballin M G.A high fidelity real-time simulation of a small turboshaft engine[R].NASA-TM-100991,1988.

    [4]Caliskan F,Hajiyev C M.Aircraft sensor fault diagnosis based on Kalman filter innovation sequence[C]//Proceedings of the 37th IEEE Conference,1998:42-67.

    [5]Kobayashi T,Simon D L.Application of a bank of kalman filters for aircraft engine fault diagnostics[R].ASME 2003-GT-38550.

    [6]Al-Hamdan Q Z,Ebaid M.Modeling and simulation of a gas turbine engine for power generation[J].ASME Journal of Engineering for Gas Turbines and Power,2006,128(2):302-311.

    [7]覃道亮,何皚,孔祥興,等.基于UIO的航空發(fā)動機控制系統(tǒng)傳感器故障診斷[J].航空動力學報,2011,26(6):1396-1404.TAN Daoliang,HE Ai,KONG Xiangxing,et al.UIO-based sensor fault diagnosis for aero-engine control systems[J].Journal of Aerospace Power,2011,26(6):1396-1404.(in Chinese)

    [8]張鵬.基于卡爾曼濾波的航空發(fā)動機故障診斷技術研究 [D].南京:南京航空航天大學,2009.ZHANG Peng.Aero engine fault diagnosis based on kalman filter[D].Nanjing:Nanjing University of Aeronautics and Astronautics,2009.(in Chinese)

    [9]方昌德.發(fā)動機狀態(tài)監(jiān)視和故障診斷系統(tǒng)的發(fā)展 [J].國際航空,2005(6):66-68.FANG Changde.Condition monitoring and fault diagnosis system of engine[J].International Aviation,2005(6):66-68.(in Chinese)

    [10]Luppold R H,Gallops G W,Kerr L J,et al.Estimating in-flight engine performance variations using kalman filter concepts[R].AIAA-89-2584.

    [11]陳毅.航空發(fā)動機控制系統(tǒng)傳感器故障診斷研究 [D].南京:南京航空航天大學,2007.CHEN Yi.Research on sensor fault diagnosis of aero engine control system[D].Nanjing:Nanjing University of Aeronautics and Astronautics,2007.(in Chinese)

    [12]Garg S.NASA Glenn research in controls and diagnostics for intelligent aerospace propulsion systems[R].NASA-TM-2005-214036.

    [13]Howlett J J.UH-60A Black hawk engineering simulation program:volume I-mathematical model[R].NASA-CR-166309,1981.

    [14]Szuch J R.Advancements in real-time engine simulation technology[R].NASA-TM-82825,1980.

    [15]Rama K.Yedavalli,Praveen Shankar,Majid Siddiqi.Modeling,Diagnostics and prognostics ofa two-spoolturbofan engine[R].AIAA-2005-4344.

    A Promoted Design Approach for Unknown Input Observer

    LONG Yi-fu1,2,WANG Xi1,2
    (1.School of Energy and Power Engineering,Beijing University of Aeronautics and Astronautics,2.Collaborative Innovation Center for Advanced Aero-Engine,Beijing 100191,China)

    This paper presents a promoted Unknown Input Observer(UIO)design approach to solve linear correlation problem while designing UIO.It is a generalized method in UIO designing.Through the linear correlation analysis of measurement output,the measured output can be reasonable grouping so that each group contains partially linearly independent output and ensure that these groups cover all the measurement output sensors that need to be diagnosed.Based on those conditions above,the existence condition of those UIOs are met and they were designed for different groups,those UIOs were used to carry out on the sensor fault diagnosis algorithm.Ultimately,a case study in aero engine has employed for confirming the feasibility of this new approach.It has been certificated that the new approach has extensively good engineering application result.

    UIO;fault diagnosis;linear correlation;sensor

    V233.7+1

    A

    10.13477/j.cnki.aeroengine.2017.05.002

    2017-03-27 基金項目:國家重大基礎研究項目資助

    龍一夫(1992),男,在讀碩士研究生,研究方向為航空發(fā)動機數字控制;E-mail:357773591@qq.com。

    龍一夫,王曦.1種未知輸入觀測器的推廣設計方法[J].航空發(fā)動機,2017,43(5):7-13.LONG Yifu,WANG Xi.A promoted design approach for unknown input observer[J].Aeroengine,2017,43(5):7-13.

    (編輯:李華文)

    猜你喜歡
    故障診斷發(fā)動機測量
    把握四個“三” 測量變簡單
    發(fā)動機空中起動包線擴展試飛組織與實施
    滑動摩擦力的測量和計算
    滑動摩擦力的測量與計算
    測量
    因果圖定性分析法及其在故障診斷中的應用
    新一代MTU2000發(fā)動機系列
    基于LCD和排列熵的滾動軸承故障診斷
    基于WPD-HHT的滾動軸承故障診斷
    機械與電子(2014年1期)2014-02-28 02:07:31
    高速泵的故障診斷
    河南科技(2014年3期)2014-02-27 14:05:48
    91在线精品国自产拍蜜月| 国产精品不卡视频一区二区| 内地一区二区视频在线| 亚洲国产欧美日韩在线播放| 国产精品99久久99久久久不卡 | 91久久精品国产一区二区三区| 精品久久蜜臀av无| 一区二区av电影网| 五月伊人婷婷丁香| 国产精品无大码| 夜夜看夜夜爽夜夜摸| 大陆偷拍与自拍| 性高湖久久久久久久久免费观看| 久久av网站| 久久毛片免费看一区二区三区| 天天影视国产精品| 夜夜爽夜夜爽视频| 久久久国产欧美日韩av| 又黄又爽又刺激的免费视频.| 国产色爽女视频免费观看| 亚洲少妇的诱惑av| 午夜福利视频在线观看免费| av在线app专区| 精品亚洲成国产av| av在线老鸭窝| 黄色欧美视频在线观看| h视频一区二区三区| 久久久精品免费免费高清| 国产午夜精品一二区理论片| h视频一区二区三区| 亚洲精品久久成人aⅴ小说 | 欧美xxxx性猛交bbbb| 91久久精品国产一区二区三区| 国产精品偷伦视频观看了| 蜜桃在线观看..| 秋霞在线观看毛片| 国产片内射在线| 久久免费观看电影| 国产色婷婷99| 久久久久精品性色| 国产国拍精品亚洲av在线观看| 母亲3免费完整高清在线观看 | 国产日韩欧美在线精品| 男人爽女人下面视频在线观看| 97精品久久久久久久久久精品| 欧美丝袜亚洲另类| 免费黄网站久久成人精品| 丝袜喷水一区| 在线观看免费视频网站a站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 我的老师免费观看完整版| 汤姆久久久久久久影院中文字幕| 国产 一区精品| 成人国产麻豆网| 中文精品一卡2卡3卡4更新| 一级a做视频免费观看| av在线播放精品| 春色校园在线视频观看| 日韩亚洲欧美综合| a级毛色黄片| 亚洲精品一区蜜桃| 国产免费又黄又爽又色| 久久久精品免费免费高清| 尾随美女入室| 香蕉精品网在线| 五月天丁香电影| 国产乱来视频区| 一个人免费看片子| 91精品一卡2卡3卡4卡| 日本黄色日本黄色录像| 久久国产精品大桥未久av| 纵有疾风起免费观看全集完整版| 亚洲激情五月婷婷啪啪| 国产一区二区三区av在线| 亚洲国产av影院在线观看| 国产亚洲午夜精品一区二区久久| 国产亚洲最大av| 国精品久久久久久国模美| 国精品久久久久久国模美| 国模一区二区三区四区视频| 午夜91福利影院| 精品人妻熟女av久视频| 麻豆精品久久久久久蜜桃| 特大巨黑吊av在线直播| 天天影视国产精品| 91精品三级在线观看| 欧美最新免费一区二区三区| 边亲边吃奶的免费视频| 99热这里只有精品一区| av播播在线观看一区| 成人二区视频| 交换朋友夫妻互换小说| av卡一久久| 欧美最新免费一区二区三区| 久久久久久久精品精品| av播播在线观看一区| 中文欧美无线码| 日日撸夜夜添| 中文字幕最新亚洲高清| 卡戴珊不雅视频在线播放| 99九九在线精品视频| 亚洲国产精品成人久久小说| 久久99热6这里只有精品| 97在线人人人人妻| 久热这里只有精品99| 欧美老熟妇乱子伦牲交| 国产成人精品福利久久| 欧美日韩一区二区视频在线观看视频在线| 男女啪啪激烈高潮av片| 最新中文字幕久久久久| 久久久久精品久久久久真实原创| 女人精品久久久久毛片| 成人影院久久| av国产精品久久久久影院| 下体分泌物呈黄色| 狠狠精品人妻久久久久久综合| 亚洲精品久久成人aⅴ小说 | 日本黄大片高清| 免费人妻精品一区二区三区视频| 3wmmmm亚洲av在线观看| 亚洲av福利一区| 精品少妇久久久久久888优播| 免费观看的影片在线观看| 五月天丁香电影| 如何舔出高潮| 国产色婷婷99| 久久久久国产网址| 狠狠婷婷综合久久久久久88av| 蜜臀久久99精品久久宅男| 日韩精品有码人妻一区| av又黄又爽大尺度在线免费看| 免费大片18禁| 日韩欧美精品免费久久| 狂野欧美激情性bbbbbb| 国产亚洲av片在线观看秒播厂| 午夜免费观看性视频| 免费人成在线观看视频色| 夜夜爽夜夜爽视频| 日韩,欧美,国产一区二区三区| 天天操日日干夜夜撸| 国产熟女午夜一区二区三区 | 极品少妇高潮喷水抽搐| 欧美变态另类bdsm刘玥| 久久久久久久久久久免费av| 精品人妻偷拍中文字幕| 只有这里有精品99| 777米奇影视久久| 啦啦啦视频在线资源免费观看| 在线免费观看不下载黄p国产| 欧美成人午夜免费资源| 18+在线观看网站| 99视频精品全部免费 在线| 国产午夜精品久久久久久一区二区三区| 极品人妻少妇av视频| 亚洲国产精品一区三区| 插阴视频在线观看视频| 久久精品久久久久久久性| 一个人免费看片子| 欧美日本中文国产一区发布| 全区人妻精品视频| 欧美日韩在线观看h| 免费久久久久久久精品成人欧美视频 | 最近手机中文字幕大全| 欧美丝袜亚洲另类| 久久久精品免费免费高清| 亚洲国产精品一区二区三区在线| 精品国产一区二区久久| 纯流量卡能插随身wifi吗| 久久国产精品男人的天堂亚洲 | 亚洲性久久影院| 日韩av在线免费看完整版不卡| 日韩三级伦理在线观看| 高清不卡的av网站| 青春草国产在线视频| videos熟女内射| 99国产精品免费福利视频| 男女啪啪激烈高潮av片| 一本色道久久久久久精品综合| 国产高清有码在线观看视频| 国产精品久久久久久精品电影小说| 亚洲图色成人| 免费少妇av软件| 久久国产精品大桥未久av| 日产精品乱码卡一卡2卡三| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品av麻豆狂野| 最近的中文字幕免费完整| 十八禁高潮呻吟视频| 十八禁网站网址无遮挡| 免费高清在线观看视频在线观看| 大码成人一级视频| 好男人视频免费观看在线| 免费看av在线观看网站| 国产色婷婷99| 九九在线视频观看精品| 制服诱惑二区| 国产片内射在线| 国产精品久久久久久精品电影小说| 新久久久久国产一级毛片| 国产精品一区二区三区四区免费观看| 亚洲精品久久久久久婷婷小说| 国产精品人妻久久久影院| 卡戴珊不雅视频在线播放| 在线观看www视频免费| 国产淫语在线视频| 免费观看在线日韩| 国产伦精品一区二区三区视频9| 啦啦啦中文免费视频观看日本| 韩国高清视频一区二区三区| 日韩欧美一区视频在线观看| 99视频精品全部免费 在线| 人人妻人人澡人人看| 最近中文字幕高清免费大全6| 国产精品一区www在线观看| 少妇精品久久久久久久| 夫妻性生交免费视频一级片| 国产黄频视频在线观看| 国产精品欧美亚洲77777| 亚洲精品日本国产第一区| 国产精品久久久久久久久免| 少妇被粗大的猛进出69影院 | 午夜福利视频在线观看免费| .国产精品久久| 菩萨蛮人人尽说江南好唐韦庄| 国产日韩一区二区三区精品不卡 | 99久久精品一区二区三区| 男女免费视频国产| 中国三级夫妇交换| 亚洲高清免费不卡视频| 蜜桃国产av成人99| 91精品三级在线观看| 天天躁夜夜躁狠狠久久av| 日韩三级伦理在线观看| 精品久久久精品久久久| 亚洲美女搞黄在线观看| 国产探花极品一区二区| 18在线观看网站| 男女高潮啪啪啪动态图| 简卡轻食公司| 国产高清三级在线| 久久国产精品大桥未久av| 成人国产麻豆网| 在线观看免费高清a一片| 韩国高清视频一区二区三区| 国产av一区二区精品久久| 美女大奶头黄色视频| 一级a做视频免费观看| 国产精品99久久99久久久不卡 | 少妇高潮的动态图| 国产精品成人在线| 久久久久人妻精品一区果冻| 午夜激情av网站| 精品视频人人做人人爽| 99精国产麻豆久久婷婷| 少妇人妻久久综合中文| 美女福利国产在线| 三上悠亚av全集在线观看| videosex国产| 精品一区二区免费观看| 久久 成人 亚洲| 97在线人人人人妻| 精品国产露脸久久av麻豆| 一本—道久久a久久精品蜜桃钙片| 中文字幕精品免费在线观看视频 | 欧美97在线视频| av网站免费在线观看视频| 一边亲一边摸免费视频| 国产 精品1| 欧美老熟妇乱子伦牲交| 一区二区三区精品91| 黑人欧美特级aaaaaa片| 成人二区视频| 亚洲国产日韩一区二区| 亚洲国产av新网站| videossex国产| 精品亚洲成a人片在线观看| 欧美日韩一区二区视频在线观看视频在线| 日韩中文字幕视频在线看片| av不卡在线播放| 久久这里有精品视频免费| 亚洲精品自拍成人| 99九九线精品视频在线观看视频| 亚洲精品456在线播放app| 亚洲,欧美,日韩| 久久久久精品久久久久真实原创| 好男人视频免费观看在线| 少妇精品久久久久久久| 国产精品 国内视频| 看非洲黑人一级黄片| 欧美日韩视频高清一区二区三区二| 亚洲av不卡在线观看| 欧美最新免费一区二区三区| 成人国语在线视频| 日韩亚洲欧美综合| 人妻 亚洲 视频| 亚洲一区二区三区欧美精品| 天堂8中文在线网| 婷婷色av中文字幕| 欧美激情国产日韩精品一区| 视频区图区小说| 国产精品久久久久成人av| 亚洲第一av免费看| 久久久久久久久久久久大奶| 啦啦啦中文免费视频观看日本| 国产在线一区二区三区精| 日韩熟女老妇一区二区性免费视频| 少妇人妻 视频| 亚洲综合精品二区| tube8黄色片| 午夜日本视频在线| 精品少妇久久久久久888优播| 亚洲不卡免费看| 欧美xxⅹ黑人| 成人影院久久| 亚洲av综合色区一区| 亚洲中文av在线| 日日撸夜夜添| 国产老妇伦熟女老妇高清| 97超碰精品成人国产| 看非洲黑人一级黄片| 国语对白做爰xxxⅹ性视频网站| 精品一区二区三卡| 汤姆久久久久久久影院中文字幕| xxxhd国产人妻xxx| 亚洲不卡免费看| 永久网站在线| 大片电影免费在线观看免费| 亚洲人与动物交配视频| 精品少妇黑人巨大在线播放| 日韩精品有码人妻一区| 色婷婷久久久亚洲欧美| 亚洲欧美成人精品一区二区| 免费高清在线观看视频在线观看| 久久精品久久精品一区二区三区| 在线亚洲精品国产二区图片欧美 | 久久青草综合色| 极品少妇高潮喷水抽搐| 亚洲激情五月婷婷啪啪| 久久久久久久精品精品| 久久久国产欧美日韩av| 日韩一本色道免费dvd| 天堂8中文在线网| 精品亚洲成国产av| 久久久午夜欧美精品| 91在线精品国自产拍蜜月| 中文字幕人妻熟人妻熟丝袜美| 伊人亚洲综合成人网| 久久青草综合色| 国产高清三级在线| 婷婷色麻豆天堂久久| 国模一区二区三区四区视频| 秋霞伦理黄片| 欧美bdsm另类| 国产精品嫩草影院av在线观看| 18+在线观看网站| 3wmmmm亚洲av在线观看| 嫩草影院入口| 两个人免费观看高清视频| 久久精品久久久久久噜噜老黄| 一本大道久久a久久精品| 欧美老熟妇乱子伦牲交| 少妇丰满av| 亚洲欧美中文字幕日韩二区| 99国产精品免费福利视频| 我的女老师完整版在线观看| 午夜激情av网站| 亚洲精品美女久久av网站| 国产一区亚洲一区在线观看| 久久人人爽人人片av| 成人免费观看视频高清| 爱豆传媒免费全集在线观看| 免费黄色在线免费观看| 岛国毛片在线播放| 2018国产大陆天天弄谢| 天堂俺去俺来也www色官网| 美女脱内裤让男人舔精品视频| 一级片'在线观看视频| 国产精品久久久久久精品古装| 免费高清在线观看日韩| 国产精品欧美亚洲77777| 18禁动态无遮挡网站| 国产免费一级a男人的天堂| 久久99热这里只频精品6学生| 精品亚洲成国产av| 婷婷色麻豆天堂久久| 国精品久久久久久国模美| 日韩一区二区视频免费看| 又大又黄又爽视频免费| 国产老妇伦熟女老妇高清| 免费看光身美女| 又粗又硬又长又爽又黄的视频| 国产精品久久久久久久久免| 亚洲第一区二区三区不卡| 亚洲国产精品一区三区| 中文字幕人妻熟人妻熟丝袜美| 日产精品乱码卡一卡2卡三| 丰满少妇做爰视频| 蜜桃在线观看..| 国产精品国产三级国产av玫瑰| 国产精品一二三区在线看| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产精品一区三区| 国产日韩欧美亚洲二区| 婷婷色麻豆天堂久久| 王馨瑶露胸无遮挡在线观看| 久久毛片免费看一区二区三区| 精品国产乱码久久久久久小说| 性色avwww在线观看| 在线精品无人区一区二区三| 我的女老师完整版在线观看| 妹子高潮喷水视频| 大话2 男鬼变身卡| 国产日韩一区二区三区精品不卡 | 成人亚洲精品一区在线观看| 97在线视频观看| 自线自在国产av| 一区二区三区精品91| 欧美日韩国产mv在线观看视频| 日韩成人伦理影院| 亚洲国产欧美在线一区| 国产69精品久久久久777片| 日本-黄色视频高清免费观看| 青青草视频在线视频观看| 97在线人人人人妻| 亚洲精品,欧美精品| 美女cb高潮喷水在线观看| 亚洲精品美女久久av网站| 国产亚洲午夜精品一区二区久久| 久久久久久人妻| 国产精品偷伦视频观看了| 欧美 亚洲 国产 日韩一| 亚洲性久久影院| 免费久久久久久久精品成人欧美视频 | 中文天堂在线官网| 啦啦啦视频在线资源免费观看| 天天躁夜夜躁狠狠久久av| 卡戴珊不雅视频在线播放| 久久狼人影院| 69精品国产乱码久久久| 成人手机av| 一区在线观看完整版| 欧美日韩av久久| 美女脱内裤让男人舔精品视频| 校园人妻丝袜中文字幕| 久久精品熟女亚洲av麻豆精品| 国产国语露脸激情在线看| 青春草国产在线视频| 在线看a的网站| 亚洲成人手机| 少妇的逼水好多| 色吧在线观看| 欧美 日韩 精品 国产| freevideosex欧美| 国产老妇伦熟女老妇高清| 熟女电影av网| 日本爱情动作片www.在线观看| 国产又色又爽无遮挡免| 日韩三级伦理在线观看| 丝瓜视频免费看黄片| 亚洲色图 男人天堂 中文字幕 | 久久久久久久精品精品| 青春草视频在线免费观看| 桃花免费在线播放| 亚洲人成网站在线播| 黄片无遮挡物在线观看| 亚洲精品第二区| 交换朋友夫妻互换小说| 99久久精品国产国产毛片| 久久女婷五月综合色啪小说| 国产综合精华液| 国产成人午夜福利电影在线观看| 少妇的逼好多水| 久久久a久久爽久久v久久| 日韩制服骚丝袜av| 亚洲国产欧美日韩在线播放| 只有这里有精品99| 九色亚洲精品在线播放| 99精国产麻豆久久婷婷| 在线观看免费高清a一片| 国产在线视频一区二区| 中国国产av一级| 青春草视频在线免费观看| 久久久久久久久大av| 中文字幕制服av| av在线观看视频网站免费| 亚洲精品第二区| 22中文网久久字幕| 色婷婷久久久亚洲欧美| 中文字幕人妻丝袜制服| 人妻人人澡人人爽人人| 国产欧美亚洲国产| 视频中文字幕在线观看| 亚洲无线观看免费| 亚洲精品国产色婷婷电影| 国产淫语在线视频| 久久人人爽人人爽人人片va| 国产亚洲av片在线观看秒播厂| 美女国产视频在线观看| 亚洲人成网站在线播| 日韩中字成人| 嘟嘟电影网在线观看| 久久99精品国语久久久| 亚洲精品第二区| 久久久亚洲精品成人影院| 丝袜脚勾引网站| 美女大奶头黄色视频| 国产精品久久久久久精品电影小说| 天堂中文最新版在线下载| 国产亚洲精品久久久com| 国产免费现黄频在线看| 啦啦啦中文免费视频观看日本| 人人妻人人爽人人添夜夜欢视频| 狂野欧美激情性xxxx在线观看| 午夜免费男女啪啪视频观看| 精品久久久久久久久av| 亚洲三级黄色毛片| 在线观看三级黄色| 少妇精品久久久久久久| 一本—道久久a久久精品蜜桃钙片| 高清在线视频一区二区三区| 亚洲欧洲国产日韩| 自线自在国产av| 少妇人妻 视频| 乱人伦中国视频| 亚洲国产精品专区欧美| 亚洲美女搞黄在线观看| 搡女人真爽免费视频火全软件| 免费黄网站久久成人精品| 精品国产一区二区久久| 特大巨黑吊av在线直播| 制服诱惑二区| 亚洲精品中文字幕在线视频| 你懂的网址亚洲精品在线观看| 欧美日韩视频高清一区二区三区二| 日日摸夜夜添夜夜添av毛片| 国产精品人妻久久久影院| 亚洲精品456在线播放app| 午夜av观看不卡| 香蕉精品网在线| 少妇人妻 视频| 日日啪夜夜爽| 国产成人91sexporn| 丁香六月天网| 色吧在线观看| tube8黄色片| 一边亲一边摸免费视频| 少妇被粗大的猛进出69影院 | 亚洲欧美成人综合另类久久久| 国产深夜福利视频在线观看| 少妇人妻 视频| 新久久久久国产一级毛片| 亚洲一级一片aⅴ在线观看| 国产色爽女视频免费观看| 97超视频在线观看视频| 国产免费又黄又爽又色| 一边亲一边摸免费视频| 午夜免费观看性视频| 日韩强制内射视频| 欧美人与善性xxx| 日本av手机在线免费观看| 亚洲精品美女久久av网站| 久久97久久精品| 性高湖久久久久久久久免费观看| 午夜福利视频精品| 国产成人精品在线电影| 男男h啪啪无遮挡| 亚洲av二区三区四区| 精品少妇内射三级| 大陆偷拍与自拍| 热99国产精品久久久久久7| 欧美精品一区二区大全| 国产老妇伦熟女老妇高清| 22中文网久久字幕| 国产日韩欧美视频二区| 最新的欧美精品一区二区| 三级国产精品片| 欧美亚洲 丝袜 人妻 在线| 久久久精品区二区三区| 一区在线观看完整版| 亚洲精品国产av蜜桃| 美女国产视频在线观看| 亚洲欧美色中文字幕在线| 2018国产大陆天天弄谢| 日韩欧美一区视频在线观看| 亚洲av电影在线观看一区二区三区| 99re6热这里在线精品视频| 黄色毛片三级朝国网站| 人体艺术视频欧美日本| av女优亚洲男人天堂| 欧美日韩成人在线一区二区| 超碰97精品在线观看| 草草在线视频免费看| 久久久久久久亚洲中文字幕| 日韩中文字幕视频在线看片| 日韩一区二区视频免费看| 9色porny在线观看| 久久久国产精品麻豆| 高清欧美精品videossex| 91午夜精品亚洲一区二区三区| 国产亚洲一区二区精品| 69精品国产乱码久久久| 97在线人人人人妻| 又粗又硬又长又爽又黄的视频| 精品少妇黑人巨大在线播放| 免费观看无遮挡的男女| 人体艺术视频欧美日本| 国产精品久久久久久精品电影小说| 婷婷色麻豆天堂久久| 纯流量卡能插随身wifi吗| 最近最新中文字幕免费大全7| 日韩一区二区视频免费看| 欧美日韩视频高清一区二区三区二| 亚洲中文av在线| 99久国产av精品国产电影| 母亲3免费完整高清在线观看 |