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

    風(fēng)能轉(zhuǎn)換系統(tǒng)隨機建模與H∞容錯控制

    2015-02-23 10:47:22史運濤侯彥嬌孫德輝李正熙
    電機與控制學(xué)報 2015年3期
    關(guān)鍵詞:發(fā)電機風(fēng)速故障

    史運濤,侯彥嬌,孫德輝,李正熙

    (北方工業(yè)大學(xué)現(xiàn)場總線技術(shù)與自動化北京市重點實驗室,北京 100144)

    風(fēng)能轉(zhuǎn)換系統(tǒng)隨機建模與H∞容錯控制

    史運濤,侯彥嬌,孫德輝,李正熙

    (北方工業(yè)大學(xué)現(xiàn)場總線技術(shù)與自動化北京市重點實驗室,北京 100144)

    針對帶有傳感器、執(zhí)行器故障的風(fēng)能轉(zhuǎn)換系統(tǒng)(WECS)的建模與容錯控制問題,提出了基于隨機分段仿射(PWA)模型的建模方法和H∞容錯控制方法。建立了風(fēng)電機組多工作區(qū)域的隨機PWA正常模型;針對WECS的槳距角傳感器故障、功率測量傳感器故障、槳距角執(zhí)行器故障、電機電磁轉(zhuǎn)矩執(zhí)行器故障,建立風(fēng)能轉(zhuǎn)換系統(tǒng)的隨機PWA故障模型;基于WECS系統(tǒng)的正常和故障的隨機PWA模型,利用非線性耗散系統(tǒng)的穩(wěn)定性理論與H∞魯棒控制方法,提出并證明了WECS系統(tǒng)的H∞容錯控制器存在性定理,給出控制器求解方法。仿真結(jié)果表明本文所提出的基于隨機PWA模型的風(fēng)能轉(zhuǎn)換系統(tǒng)的H∞容錯控制方法對于存在的傳感器、執(zhí)行器增益故障具有較好的容錯控制性能。

    容錯控制;隨機PWA;風(fēng)能轉(zhuǎn)換系統(tǒng);H∞控制

    0 引言

    風(fēng)能作為重要的清潔能源,得到世界各國的重視。風(fēng)能轉(zhuǎn)換系統(tǒng)是由風(fēng)輪葉片系統(tǒng)、齒輪箱主傳動鏈系統(tǒng)、發(fā)電機系統(tǒng)以及其控制系統(tǒng)組成,是典型的機、電、液一體的復(fù)雜的非線性動力學(xué)系統(tǒng),在隨機、間歇性的風(fēng)力作用下表現(xiàn)出隨機、切換的特性。對于這樣的系統(tǒng)實現(xiàn)建模與控制是一個比較困難的問題;同時,風(fēng)力發(fā)電機組的控制系統(tǒng)中存在傳感器、執(zhí)行器等常見故障[1-2],控制系統(tǒng)的故障嚴(yán)重的影響風(fēng)力發(fā)電機組的性能與安全。因此,風(fēng)能轉(zhuǎn)換系統(tǒng)的容錯控制技術(shù)的研究對于風(fēng)電機組的安全、高效運行具有非常重要的意義。

    近年來,國內(nèi)外已經(jīng)開始了對于風(fēng)力發(fā)電系統(tǒng)的容錯控制研究[2-9]:在文獻[6]中,Sloth等人利用線性變參數(shù)(LPV)控制器設(shè)計方法設(shè)計了風(fēng)電機組的主動和被動容錯控制器;在文獻[7]中,作者提出了雙饋式感應(yīng)發(fā)電機的基于模型的故障檢測和控制回路重組的主動容錯方法;文獻[8]給出了直驅(qū)式風(fēng)力發(fā)電機變流器的容錯控制方法;在文獻[9]中,作者設(shè)計了一種針對變槳距控制系統(tǒng)的容錯開關(guān)磁阻電機。目前對于風(fēng)力發(fā)電機組的建模與容錯控制的研究的主要缺點是沒有深入描述風(fēng)力發(fā)電機組的隨機、非線性切換的動力學(xué)特性。一些研究者僅考慮非線性特性,利用滑模變結(jié)構(gòu)、反饋線性化等技術(shù)手段設(shè)計非線性控制器,但是系統(tǒng)建模中沒有考慮其隨機與切換特性。就作者所知,還沒有在隨機、非線性動力學(xué)建??蚣芟拢槍︼L(fēng)能轉(zhuǎn)換系統(tǒng)的傳感器與執(zhí)行器故障的容錯控制器設(shè)計方法的研究。

    分段仿射系統(tǒng)是分析和設(shè)計非線性系統(tǒng)的強有力工具,很多非線性系統(tǒng)可以用PWA系統(tǒng)以任意精度進行逼近,PWA系統(tǒng)還與幾類典型的混雜系統(tǒng)(如混合邏輯動態(tài)系統(tǒng)、線性混合自動機等)等價;隨機的PWA模型則是在PWA模型的基礎(chǔ)上考慮了系統(tǒng)受到隨機信號的干擾作用下的情況,是對于PWA模型的一個特別有意義的擴展,從而大大的增強其隨機建模的能力。目前在保證隨機PWA系統(tǒng)穩(wěn)定和H∞控制方面已經(jīng)取得一系列的成果[11-15]。Cuzzola在文獻[15]中討論了離散PWA系統(tǒng)H∞控制器設(shè)計方法,但結(jié)論只適用于沒有仿射項的系統(tǒng);隨后Feng和Xu分別在文獻[10]、[16]中重新給出結(jié)論適合于具有仿射項的系統(tǒng)的離散PWA的H∞控制器的設(shè)計方法。

    本文利用隨機PWA模型,將非高斯的風(fēng)力信號(符合威布爾分布)轉(zhuǎn)化為湍流信號與一個平均風(fēng)速的疊加,從而把風(fēng)速信號進行分解:低頻部分(平均風(fēng)速)轉(zhuǎn)化為PWA模型的位移常數(shù)項(displacement),高頻部分(湍流部分)表示為普通高斯信號經(jīng)過(與該點平均風(fēng)速相關(guān))二階成型濾波器,據(jù)此建立了風(fēng)電機組多工作區(qū)域的隨機PWA模型,針對WECS的傳感器故障(槳距角傳感器故障、功率測量傳感器故障)和執(zhí)行器故障(槳距角執(zhí)行器故障、電機電磁轉(zhuǎn)矩執(zhí)行器故障)建立風(fēng)能轉(zhuǎn)換系統(tǒng)的隨機PWA故障模型;在隨機PWA模型框架下,推導(dǎo)并證明了H∞容錯控制器存在條件和設(shè)計方法。本文由以下幾部分組成:第1部分簡單描述了WECS的隨機、非線性動態(tài)模型和切換控制策略; WECS的隨機PWA故障建模方法在第2部分中做了介紹;針對WECS的傳感器故障、執(zhí)行器故障的H∞容錯控制方法在第3部分中展開討論;仿真結(jié)果和相關(guān)分析在第4部分給出;最后給出結(jié)論部分。

    1 WECS的動態(tài)模型和控制策略

    WECS(圖1)的輸入為風(fēng)速v(t),槳距角設(shè)定值βref(t)和電磁轉(zhuǎn)矩設(shè)定值Tgref(t),系統(tǒng)的輸出為發(fā)電機功率Pg(t)和高速軸轉(zhuǎn)速ωg(t)。

    圖1 WECS的結(jié)構(gòu)圖Fig.1The structure of WECS

    1.1 風(fēng)模型

    風(fēng)由兩部分疊加產(chǎn)生[13],如式(1)所示,vm(t)是低頻部分(描述長期,低頻變化量),代表平均風(fēng)速;vs(t)是湍流部分(對應(yīng)于快速,高頻變化量)。

    1.2 空氣動力學(xué)模型

    穿過整個風(fēng)輪旋轉(zhuǎn)平面面積的風(fēng)的有效功率表示為:

    式中:A為風(fēng)輪旋轉(zhuǎn)平面面積,Pw為風(fēng)的有效功率,v為風(fēng)輪采集到的有效風(fēng)速,ρ為空氣密度,被考慮為常量?;陲L(fēng)的有效功率,風(fēng)輪捕獲功率可根據(jù)功率系數(shù)Cp(λ,β)得出,功率系數(shù)的大小依賴于葉尖速比λ和槳距角β。風(fēng)輪捕獲功率為:

    式中:Pa為轉(zhuǎn)子捕獲功率。

    葉尖速比定義為葉尖速度和風(fēng)輪有效風(fēng)速的比值:

    式中:ωr是低速軸轉(zhuǎn)速。

    1.3 傳動鏈模型

    傳動鏈?zhǔn)莻鬟f風(fēng)能的主要裝置,由低速軸、高速軸、齒輪以及彈性裝置組成。傳動鏈動態(tài)方程如下:

    式中:Jr和Jg分別為低速軸和高速軸的轉(zhuǎn)動慣量,Ks是傳動鏈的彈性裝置的剛性系數(shù),Ds為傳動鏈的彈性裝置的阻尼系數(shù),Ng為齒輪比,δ為柔性傳動鏈的扭轉(zhuǎn)度

    1.4 槳距系統(tǒng)模型

    變槳距控制系統(tǒng)可以由一個二階模型來表示[14],其狀態(tài)方程為:

    1.5 發(fā)電機和變流器模型

    就整體而言,風(fēng)力發(fā)電機與整流器模型可以用一階模型代替,即電磁子系統(tǒng):

    式中:τg為時間常數(shù)。時發(fā)功率可描述為:

    綜上,由各個部分的模型可以得到整個風(fēng)電機組的完整模型:

    1.6 控制策略

    風(fēng)電機組的控制策略如圖2所示。

    圖2 風(fēng)電機組的控制策略Fig.2The control strategy for wind turbine

    風(fēng)電系統(tǒng)在不同的風(fēng)速時對功率和轉(zhuǎn)速均有不同的控制要求,一般來說在風(fēng)速達不到系統(tǒng)切入風(fēng)速時,風(fēng)機不啟動;在風(fēng)速未達到額定風(fēng)速前,風(fēng)機需最大效率的捕獲風(fēng)能;當(dāng)風(fēng)速超過額定風(fēng)速則需要限制功率上升并保持系統(tǒng)穩(wěn)定。

    如圖3在A—B段即切入風(fēng)速前,系統(tǒng)處于停機狀態(tài);

    圖3 β、Tg和Pg隨風(fēng)速的變化趨勢Fig.3The tendency of the β、Tgand Pgfollowing the wind speed

    在B-C段又稱部分負荷段,控制的主要目標(biāo)在于功率的最大捕獲。通過調(diào)節(jié)槳距角和電磁轉(zhuǎn)矩使得功率系數(shù)Cp(λ,β)獲得最大值,當(dāng)風(fēng)速繼續(xù)增大時,應(yīng)注意控制高速軸轉(zhuǎn)速ωg的值不要超過電機的閾值轉(zhuǎn)速。

    在C—D段又稱滿負荷段,此時的控制目標(biāo)是保持時發(fā)功率Pg維持在額定值Pgnom,同時限制高速軸轉(zhuǎn)速ωg維持在[ωgnomωgmax]之間。變槳控制的主要目的則變?yōu)橄拗骑L(fēng)能的捕獲即限制功率,避免增加系統(tǒng)的機械負荷。

    在D—E段接近切出風(fēng)速后,改變槳距角的值以盡量減少風(fēng)能的捕獲。同時系統(tǒng)與電網(wǎng)斷開,抱閘停機。

    2 WECS在隨機PWA框架下的故障建模

    綜上所述,風(fēng)電系統(tǒng)主要工作在四個工作區(qū)間,本節(jié)首先介紹隨機PWA建模方法,然后給出了風(fēng)電系統(tǒng)的隨機PWA建模過程:

    2.1 隨機PWA模型框架

    離散時間隨機PWA系統(tǒng)的狀態(tài)空間方程如下所示:

    式中:xk∈ RRn為系統(tǒng)狀態(tài),uk∈ RRm為控制輸入,zk∈R Rm是一個可建模的性能輸出,所有可能的向量的集合? RRn+m既可以是n+m維的實空間 RRn+m,也可以是包含原點的多面體,{χi}si=1是的一個多面體劃分,ai∈ RRn是常數(shù)向量。

    把每一個χi當(dāng)作一個單元,假定每一個單元都是多面體形式的,由矩陣定義為如下形式:

    Sj:={所有的i,?x,u和x∈,xTu [T]T∈

    2.2 WECS的隨機PWA模型

    某一靜態(tài)工作點,可以利用泰勒級數(shù)展開的方法實現(xiàn)非線性動力學(xué)特性的線性化展開的方法,將χi}。用I={1,…,s}代表所有χ下標(biāo)i可能滿足的集合,同時J={1,…,t}代表所有下標(biāo)j可能滿足的集合;同時必須滿足傳動鏈模型進行線性化,線性化后的傳動鏈動態(tài)方程如下所示:

    由此,可以得到一個線性化的WECS整體動態(tài)狀態(tài)空間模型:

    隨機的風(fēng)力載荷提供驅(qū)動風(fēng)機的能量,可以看為平均風(fēng)速vm和隨機部分vs的疊加,平均風(fēng)速可以是一個常量。根據(jù)文獻[15-16],高頻湍流部分vs可以看作是白噪聲經(jīng)過二階濾波器之后的“點風(fēng)”,二階濾波器模擬了由風(fēng)輪葉片掃過的盤狀區(qū)域的影響。在頻域中,vs的功率譜密度Svs可以寫為Svs(f,vm)=Sp(f)Sf(f,vm),式中:Sp(f)是“點風(fēng)”的頻譜,Sf(f,vm)代表基于平均風(fēng)速的濾波器。然后風(fēng)速模型的湍流部分表示成白噪聲過程近似為一個線性二階傳遞函數(shù):

    式中:ω1=vs,e∈N(0,1),a1、a2、a3為基于平均風(fēng)速的參數(shù)。

    則式(12)中的模型可以轉(zhuǎn)化成如下形式:

    根據(jù)圖4,選擇3個工作點,風(fēng)機本身的參數(shù)如轉(zhuǎn)動慣量等可直接得到,各個工作點的的值通常由一個有效的風(fēng)力估計器根據(jù)不同的風(fēng)速求取[17-18]。

    圖4 風(fēng)速和相對應(yīng)的工作點Fig.4Wind speed vs.the corresponding working point

    其中,風(fēng)力發(fā)電機組的模型參數(shù)為:低速軸轉(zhuǎn)動慣量Jr=90 000 kg·m2,高速軸轉(zhuǎn)動慣量Jg=10 kg·m2;傳動鏈的彈性裝置的剛性系數(shù)Ks=8×106N·m/rad,傳動鏈的彈性裝置的阻尼系數(shù)Ds=8×104kg·m2/ (rad·s);時間常數(shù)τβ=0.15 s,τg=0.1 s;時發(fā)電功率額定值Pg,nom=225 kw;低速軸轉(zhuǎn)速額定值ωr,nom=4.29 rad/s,高速軸轉(zhuǎn)速額定值ωg,nom=105.534 rad/s,低速軸轉(zhuǎn)速最低值ωr,min=3.5 rad/s,高速軸轉(zhuǎn)速最低值ωg,min=86.1 rad/s。其中,風(fēng)力發(fā)電機組不同工作點的線性模型參數(shù)如表1所示。

    表1 不同工作點的線性模型參數(shù)Table 1Parameters of linearized model in different working points

    2.3 WECS的執(zhí)行器故障和傳感器故障模型

    在本文的工作中,同時考慮執(zhí)行器故障和傳感器故障。用uj表示第j個執(zhí)行器或傳感器,uFj表示第j個失效的執(zhí)行器或第j個失效的傳感器。將執(zhí)行器或者傳感器中的一種增益丟失表示為:

    式中:αj表示第j個執(zhí)行器或傳感器的失效比,αMj是對應(yīng)的最大失效比,αj=0表示對應(yīng)的執(zhí)行器或傳感器沒有失效,0<αj<1表示對應(yīng)的執(zhí)行器或傳感器部分失效;αj=1表示對應(yīng)的執(zhí)行器或傳感器完全失效。定義α=diag{α1,α2,…,αm},則uF= Γu,其中Γ=(1-α)。文中考慮的故障如表2所示。

    表2 文中考慮的故障Table 2The considered faults

    其中Γa和Γs均為二階對角矩陣。

    帶有執(zhí)行器增益丟失Γa的WECS系統(tǒng)的隨機PWA模型如下所示:

    同理,帶有傳感器增益丟失Γs的WECS系統(tǒng)的隨機的PWA模型如下所示:

    3 風(fēng)能轉(zhuǎn)換系統(tǒng)的H∞容錯控制

    3.1 執(zhí)行器故障容錯控制

    根據(jù)執(zhí)行器故障的PWA模型(16),將輸入干擾信號wk擴展為:

    本文考慮的H∞框架是基于l2有限時域定義的,因此,對干擾輸入的擴展是合理的。

    由此文獻[5]中提到的控制方法可以直接應(yīng)用到擴展系統(tǒng)(19)。在標(biāo)準(zhǔn)的假設(shè)下:

    通過這種方式,控制器可以考慮到如下置換ai=D,其中D=[0I]。

    通過將控制律(22)應(yīng)用到PWA系統(tǒng)(19),可以得到閉環(huán)PWA執(zhí)行器故障系統(tǒng):

    則系統(tǒng)(23)滿足(25)式所示的耗散不等式,同時,系統(tǒng)(23)是PWQ穩(wěn)定的。

    證明回憶初始條件x0=0,則式(25)有如下的關(guān)系式,?N≥0

    明顯的,如果條件(26)成立,那么不等式(29)也成立。另外,認為Cij≥0,所以由式(26)的元素(1,1),可以得到

    此式表明系統(tǒng)(23)是PWQ穩(wěn)定的。

    以下針對(16)式所描述的WECS系統(tǒng)中存在執(zhí)行器增益故障,尋找一種形如(22)式的狀態(tài)反饋H∞控制律實現(xiàn)系統(tǒng)被動容錯控制,主要結(jié)果可歸納為如下的定理:

    定理1對于PWA系統(tǒng)(16),存在一個形如(22)的狀態(tài)反饋控制律能夠保證系統(tǒng)PWQ李雅普諾夫穩(wěn)定并滿足如下供應(yīng)率:可以得到狀態(tài)反饋矩陣K1j(j∈J)為:

    證明通過引用舒爾引理,可以將式(26)改寫為:

    明顯看出,對于Qi>0,?i∈I,控制矩陣可以重寫為式(33)的形式。

    3.2 傳感器故障容錯控制

    定理2對于PWA系統(tǒng)(17),存在一個形如(22)的狀態(tài)反饋控制律能夠保證系統(tǒng)PWQ李雅普諾夫穩(wěn)定并滿足供應(yīng)率為

    4 仿真結(jié)果與分析

    本文在平均風(fēng)速為7.5、10、16 m/s3個工作點建立了WECS系統(tǒng)的隨機PWA模型,并設(shè)計了WECS系統(tǒng)正常情況下的H∞狀態(tài)反饋控制器;建立了WECS系統(tǒng)槳距執(zhí)行器及發(fā)電機執(zhí)行器的PWA故障模型(故障因數(shù)Γa=0.6),根據(jù)定理1設(shè)計了執(zhí)行器故障下的H∞容錯控制器;建立了WECS系統(tǒng)電機轉(zhuǎn)速傳感器及電機功率傳感器故障的隨機PWA模型(故障因數(shù)Γs=1.2),根據(jù)定理2設(shè)計了傳感器故障下的H∞容錯控制器;以下給出了WECS系統(tǒng)在三個工作點的隨機PWA模型,并給出了根據(jù)定理1和定理2計算的狀態(tài)反饋系數(shù)矩陣;分別在4.1、4.2、4.3節(jié)中給出了WECS系統(tǒng)正常、執(zhí)行器故障和傳感器故障三種情況下的仿真結(jié)果。

    根據(jù)第2部分的建模過程,可以得到WECS的如下隨機PWA模型:

    如果0<E1xk≤8,有

    狀態(tài)空間分割矩陣:E1=[0 0 0 0 0 1 0]。

    根據(jù)本文所推導(dǎo)的定理1和定理2可以成功求解WECS系統(tǒng)在正常、傳感器故障、執(zhí)行器故障情況下的分段線性狀態(tài)反饋矩陣K1,K2,K3,保證WECS系統(tǒng)滿足PWQ李亞普諾夫穩(wěn)定。三種情況下求取的反饋矩陣如下:

    WECS系統(tǒng)正常時反饋增益矩陣為

    4.1 正常WECS系統(tǒng)的H∞控制仿真結(jié)果

    在仿真中,利用式(9)所描述的WECS系統(tǒng)的非線性模型作為實際的被控對象。圖5給出了無故障的WECS系統(tǒng)在本文所設(shè)計的H∞控制器作用下的動態(tài)響應(yīng)結(jié)果。圖5(a)給出了用于測試的風(fēng)速信號,其平均風(fēng)速vm分別取值7.5、10、16m/s,在100、200、300s進行切換;圖5(b)給出了WECS系統(tǒng)在H∞控制器作用下的葉尖速比曲線,圖5(c)則給出了最大功率點跟蹤算法(MPPT)計算出的最優(yōu)發(fā)電功率設(shè)定值的曲線;圖5(d)和圖5(e)給出了發(fā)電機轉(zhuǎn)速和實時功率的響應(yīng)曲線,從圖中可以看出實時功率較好的跟蹤了功率設(shè)定值,發(fā)電機轉(zhuǎn)速控制性能良好;圖5(f)和圖5(g)給出了控制器輸出的槳距角和發(fā)電機轉(zhuǎn)矩的曲線;圖5的仿真結(jié)果表明:本文所設(shè)計的WECS系統(tǒng)基于隨機PWA模型的H∞控制器取得較好的控制效果。

    圖5 風(fēng)電機組基于隨機PWA模型的H∞可靠容錯控制仿真結(jié)果Fig.5Simulation results of H∞r(nóng)eliable control for normal WECS based on stochastic PWA model

    4.2 帶有執(zhí)行器故障的WECS模型的H∞控制仿真結(jié)果

    圖6 風(fēng)電機組存在執(zhí)行器故障但無容錯控制的仿真結(jié)果Fig.6Simulation results of actuators fault of WECS without fault tolerant control

    圖6和圖7給出了風(fēng)電系統(tǒng)存在槳距系統(tǒng)執(zhí)行器故障增益損耗Γβ及發(fā)電機執(zhí)行器故障增益ΓTg狀態(tài)下(本文中令增益損耗Γβ=ΓTg=0.6)的仿真結(jié)果;圖6給出了風(fēng)力發(fā)電系統(tǒng)在存在兩個執(zhí)行器故障,無容錯控制器作用下的響應(yīng)曲線。圖7給出了根據(jù)定理1設(shè)計的執(zhí)行器故障容錯控制器作用下的控制效果。圖6(a),(b),(c)曲線與圖5(a),(b),(c)相比較,曲線變化很小,說明MPPT算法在傳感器故障下也能很好的工作。另一方面,相對于圖5的(d)和(e)曲線,圖6中的(d)發(fā)電機功率Pg和(e)發(fā)電機速度ωg的性能惡化,并開始展現(xiàn)出失去穩(wěn)定的趨勢。圖6(f)和6(g)中的控制變量顯示出了大范圍的振蕩。由此可以看出,對于風(fēng)能轉(zhuǎn)換系統(tǒng)正常模型下的H∞控制器無法處理風(fēng)能轉(zhuǎn)換系統(tǒng)存在執(zhí)行器增益因子損失的故障。和圖6(d),(e)進行比較,圖7 (d),(e)給出了執(zhí)行器故障容錯控制器作用下的WECS系統(tǒng)的發(fā)電機轉(zhuǎn)速和發(fā)電機的功率的響應(yīng)曲線。從圖7中可以看出,存在容錯控制器比沒有容錯設(shè)計的正??刂破鞯男阅茱@著提高了,發(fā)電機功率Pg和發(fā)電機速度ωg保持穩(wěn)定的趨勢。

    圖7 風(fēng)電機組存在執(zhí)行器故障下的H∞可靠容錯控制仿真結(jié)果Fig.7Simulation results of H∞fault tolerant control for actuators fault of WECS

    4.3 帶有傳感器故障的WECS模型的H∞控制仿真結(jié)果

    圖8 風(fēng)電機組存在傳感器故障但無容錯控制的仿真結(jié)果Fig.8Simulation results of sensors fault of WECS without fault tolerant control

    圖9 風(fēng)電機組存在傳感器故障下的H∞可靠容錯控制仿真結(jié)果Fig.9Simulation results of H∞fault tolerant control for sensors fault of WECS

    圖7和圖8給出了風(fēng)電系統(tǒng)在存在電機速度傳感器故障增益Γω及發(fā)電機功率傳感器故障增益ΓP狀態(tài)下(本文中令增益故障因數(shù)Γω=ΓP=1.2)的仿真結(jié)果;圖8給出了風(fēng)力發(fā)電系統(tǒng)在存在兩個傳感器故障,無容錯控制器作用下的響應(yīng)曲線。圖9給出了根據(jù)定理2設(shè)計的傳感器故障容錯控制器作用下的控制效果。在圖8(d),(e)中發(fā)電機功率和發(fā)電機轉(zhuǎn)速的性能已經(jīng)惡化,而且完全失去穩(wěn)定,同時在圖8 (f),(g)中的控制變量也顯示出大范圍的振蕩;而圖9(f),(g)表明在容錯控制器的作用下,風(fēng)能轉(zhuǎn)換系統(tǒng)的輸出獲得了較好的性能,發(fā)電機功率很好的跟蹤了MPPT算法的輸出,發(fā)電機的轉(zhuǎn)速響應(yīng)速度快并且無任何靜態(tài)誤差,發(fā)電功率也能夠跟蹤功率的最優(yōu)設(shè)定。

    5 結(jié)論

    本文針對具有傳感器、執(zhí)行器故障的風(fēng)能轉(zhuǎn)換系統(tǒng)的建模與容錯控制問題,提出了一種隨機PWA的建模方法解決了風(fēng)能轉(zhuǎn)換系統(tǒng)(WECS)正常與故障情況的建模問題;利用非線性耗散系統(tǒng)的穩(wěn)定性理論與H∞容錯控制方法,以定理形式給出了WECS系統(tǒng)在傳感器與執(zhí)行器故障情況下H∞容錯控制器設(shè)計方法。實際的仿真結(jié)果表明本文所提出的隨機PWA建模與H∞容錯控制方法能夠很好的解決隨機風(fēng)力載荷作用下的風(fēng)能轉(zhuǎn)換系統(tǒng)的建模與容錯控制問題。

    [1]ODGAARD P F,STOUSTRUP J,KINNAERT M.Fault-tolerant control of wind turbines:a benchmark model[J].IEEE Transactions on Control Systems Technology,2013,21(4):1168-1182.

    [2]SLOTH C,ESBENSEN T,STOUSTRUP J.Robust and fault-tolerant linear parameter-varying control of wind turbines[J].Mechatronics,2011,21(4):645-659.

    [3]ROTONDO D,NEJJARI F,PUIG V,et al.Fault tolerant control of the wind turbine benchmark using virtual sensors/actuators[C]// Proceedings of the 8th IFAC Symposium on Fault Detection,Supervision and Safety of Technical Processes,August 29-31,2012,Mexico City,Mexico.2012:114-119.

    [4]ZHANG H B,DANG C H,ZHANG J.H∞control of piecewise-linear systems under unreliable communication links[J].Circuits Syst Signal Process,2012,31:1297-1318.

    [5]BURKART R,MARGELLOS K,LYGEROS J.Nonlinear control of wind turbines:an approach based on switched linear systems and feedback linearization[C]//50th IEEE Conference on Decision and Control and European Control Conference,Dec 12-15,2011,Orlando,USA.2011:5485-5490.

    [6]SLOTH C,ESBENSEN T,STOUSTRUP J.Active and passive fault-tolerant LPV control of wind turbines[C]//2010 American Control Conference,June 30-July 2,2010,Baltimore,USA.2010:4640-4646.

    [7]ROTHENHAGEN K,F(xiàn)UCH F W.Doubly fed induction generator model-based sensor fault detection and control loop reconfiguration[J].IEEE Transactions on Industrial Electronics,2009,56(10): 4229-4238.

    [8]PARKER M A,LI Ran,,F(xiàn)INNEY S J.Distributed control of a faulttolerant modular multilevel inverter for direct-drive wind turbine grid interfacing[J].IEEE Transactions on Industrial Electronics, 2013,60(2):509-522.

    [9]Ruba M,Szabo L,Jurca F.Fault tolerant switched reluctance machine for wind turbine blade pitch control[C]//2009 International Conference on Clean Electrical Power,June 9-11,2009,Capri,Italy.2009:721-726.

    [10]CUZZOLA F A,MORARI M.A generalized approach for analysis and control of discrete-time piecewise affine and hybrid systems[J].Hybrid Systems:Computation and Control-Lecture Notes in Computer Sciences 2034,2001:189-203.

    [11]JOHANSSON M,RANTZER A.Computation of piecewise quadratic Lyapunov functions for hybrid systems[J].IEEE Transactions on Automatic Control,1998,43(4):555-559.

    [12]RANTZER A,JOHANSSON M.Piecewise linear quadratic optimal control[J].IEEE Transactions on Automatic Control,2000,45 (4):629-637.

    [13]HASSIBI A,BOYD S.Quadratic stabilization and control of piecewise-linear systems[C]//Proceedings of the American Control Conference,June 21-26,1998,Philadelphia USA.1998:3659 -3664.

    [14]PETTERSSON S,LENNARTSON B.LMI approach for stability and robustness of hybrid systems[C]//Proceedings of the American Control Conference,June 4-6,1997,Albuquerque,USA.1997:1714-1718.

    [15]SLUPPHAUG O,F(xiàn)OSS B A.Constrained quadratic stabilization of discrete time uncertain nonlinear multi-model systems using piecewise affine state feedback[J].International Journal of Control,1999,72:686-701.

    [16]MIGNONE D,F(xiàn)ERRARI-TRECATE G,MORARI M.Stability and stabilization of piecewise affine and hybrid systems:an LMI approach[C]//Proceedings of the 39th IEEE Conference on Decision and Control,December 12-15,2000,Sydney,Australia.2000: 504-510.

    [17]年珩,周義杰,李嘉文.基于開繞組結(jié)構(gòu)的永磁風(fēng)力發(fā)電機控制策略[J].電機與控制學(xué)報,2013,17(4):79-85.

    NIAN Heng,ZHOU Yijie,LI Jiawen.Control strategy of permanent magnet wind generator based on open winding configuration[J].Electric Machines and Control,2013,17(4):79-85.

    [18]包廣清,鄭文鵬,毛開富.基于橫磁通發(fā)電機的永磁直驅(qū)風(fēng)力發(fā)電系統(tǒng)[J].電機與控制學(xué)報,2012,16(11):39-44.

    BAO Guangqing,ZHENG Wenpeng,MAO Kaifu.Direct-driven wind turbine system with multiphase transverse flux permanent magnet generator[J].Electric Machines and Control,2012,16 (11):39-44.

    [19]GANG F.Stability analysis of piecewise discrete-time linear systems[J].IEEE Transactions on Automatic Control,2002,47(7): 1108-1112.

    [20]ESBENSEN T,SLOTH C.Fault diagnosis and fault-tolerant control of wind turbine[D].Aalborg:Aalborg University,2009.

    [21]H?JSTRUP J.Velocity spectra in the unstable planetary boundary layer[J].Journal of the Atmospheric Sciences,1982,39(10):2239 -2248.

    [22]Munteanu I,Cutululis N-A,Bratcu A L,et al.Optimal control of wind energy systems[M].London:Springer,2008.

    [23]?STERGAARD K Z,BRATH P,STOUSTRUP J.Estimation of effective wind speed[J].Journal of Physics:Conference Series,2007,75(1):1-9.

    (編輯:張詩閣)

    Stochastic modelling and H∞fault tolerance control of WECS

    SHI Yun-tao,HOU Yan-jiao,SUN De-hui,LI Zheng-xi
    (Key Laboratory of Beijing for Field-bus Technology&Automation,North China University of Technology,Beijing 100144,China)

    Modelling and H∞fault tolerant control problem of wind energy conversion system(WECS)was studied with sensor and actuator faults.The normal stochastic piecewise affine models for WECS with multiple working points were established.Also the stochastic piece wise affine(PWA)models were proposed for WECS with generator speed sensor,power sensor gain loss faults and pitch angle actuator,generator actuator gain loss faults.A reliable piecewise linear quadratic state feedback regulator was designed for WECS with the actuator and sensor faults.A sufficient condition for existence of passive fault tolerant controller was derived and formulated as feasibility of a set of linear matrix inequalities.Simulation results of control systems confirm that they are capable of controlling the wind turbine exposed to multiple simultaneous sensors faults or actuators faults.Consequently,the reliability of wind turbines is improved.

    fault tolerant control;stochastic PWA;WECS;H∞r(nóng)obust control

    10.15938/j.emc.2015.03.016

    TM 614

    A

    1007-449X(2015)03-0100-11

    2014-03-03

    國家自然科學(xué)基金面上項目(61473002);北京市自然科學(xué)基金(4132021)

    史運濤(1975—),男,博士,副教授,研究方向為風(fēng)力發(fā)電機組的容錯控制;

    侯彥嬌(1987—),女,碩士研究生,研究方向為風(fēng)力發(fā)電機組的容錯控制;

    孫德輝(1962—),男,博士,教授,研究方向為現(xiàn)場總線自動化;

    李正熙(1955—),男,博士,教授,研究方向為電力電子與電氣傳動、新能源。

    史運濤

    猜你喜歡
    發(fā)電機風(fēng)速故障
    基于Kmeans-VMD-LSTM的短期風(fēng)速預(yù)測
    基于最優(yōu)TS評分和頻率匹配的江蘇近海風(fēng)速訂正
    海洋通報(2020年5期)2021-01-14 09:26:54
    故障一點通
    大型發(fā)電機勵磁用旋轉(zhuǎn)變換器的開發(fā)和應(yīng)用
    隨身攜帶的小發(fā)電機
    軍事文摘(2016年16期)2016-09-13 06:15:49
    奔馳R320車ABS、ESP故障燈異常點亮
    基于GARCH的短時風(fēng)速預(yù)測方法
    故障一點通
    考慮風(fēng)速分布與日非平穩(wěn)性的風(fēng)速數(shù)據(jù)預(yù)處理方法研究
    江淮車故障3例
    高清不卡的av网站| 中文资源天堂在线| 久久久久久久久久久久大奶| 在线观看三级黄色| 在线观看免费高清a一片| 国产精品一二三区在线看| 美女中出高潮动态图| 成人亚洲欧美一区二区av| 国产一区二区三区av在线| a级毛片免费高清观看在线播放| av不卡在线播放| 亚洲精品一区蜜桃| 你懂的网址亚洲精品在线观看| av免费观看日本| 成人18禁高潮啪啪吃奶动态图 | 国产综合精华液| 午夜激情福利司机影院| 人人妻人人看人人澡| 天天操日日干夜夜撸| 日韩大片免费观看网站| 在现免费观看毛片| 在线观看一区二区三区激情| 亚洲欧美一区二区三区黑人 | 赤兔流量卡办理| 肉色欧美久久久久久久蜜桃| 国产一区亚洲一区在线观看| 久久综合国产亚洲精品| 免费不卡的大黄色大毛片视频在线观看| 国产探花极品一区二区| 赤兔流量卡办理| 亚洲第一av免费看| av专区在线播放| 精品少妇黑人巨大在线播放| .国产精品久久| 中文天堂在线官网| 国产真实伦视频高清在线观看| 菩萨蛮人人尽说江南好唐韦庄| 高清av免费在线| 9色porny在线观看| 欧美日韩国产mv在线观看视频| 亚洲精品第二区| 丝瓜视频免费看黄片| 国产色婷婷99| 免费大片黄手机在线观看| 国内少妇人妻偷人精品xxx网站| 少妇人妻一区二区三区视频| 国产精品人妻久久久影院| 最后的刺客免费高清国语| 九九爱精品视频在线观看| 欧美激情极品国产一区二区三区 | 九草在线视频观看| 精华霜和精华液先用哪个| 美女xxoo啪啪120秒动态图| 建设人人有责人人尽责人人享有的| av视频免费观看在线观看| 亚洲av欧美aⅴ国产| 欧美精品国产亚洲| 国语对白做爰xxxⅹ性视频网站| 午夜精品国产一区二区电影| 美女脱内裤让男人舔精品视频| 少妇人妻 视频| 99九九线精品视频在线观看视频| 久久6这里有精品| 亚洲av福利一区| 天天躁夜夜躁狠狠久久av| 免费观看a级毛片全部| 99热网站在线观看| 又大又黄又爽视频免费| 3wmmmm亚洲av在线观看| 最新中文字幕久久久久| 爱豆传媒免费全集在线观看| 18禁在线播放成人免费| 国模一区二区三区四区视频| 国产乱来视频区| 极品少妇高潮喷水抽搐| 在线播放无遮挡| 99精国产麻豆久久婷婷| av在线观看视频网站免费| 久久久久视频综合| 欧美日韩视频高清一区二区三区二| 国模一区二区三区四区视频| 啦啦啦啦在线视频资源| 国产免费又黄又爽又色| 国产在线视频一区二区| 最新中文字幕久久久久| 男人爽女人下面视频在线观看| 国产日韩欧美在线精品| 国产黄色免费在线视频| 久久精品久久精品一区二区三区| 欧美成人精品欧美一级黄| 欧美日本中文国产一区发布| 亚洲国产精品专区欧美| 能在线免费看毛片的网站| 少妇精品久久久久久久| 一区二区三区精品91| 天美传媒精品一区二区| 亚洲av在线观看美女高潮| 欧美激情国产日韩精品一区| 亚洲久久久国产精品| 偷拍熟女少妇极品色| 我要看日韩黄色一级片| 视频中文字幕在线观看| 人妻人人澡人人爽人人| 国产老妇伦熟女老妇高清| 亚洲国产欧美日韩在线播放 | 亚洲精品日韩在线中文字幕| 日韩制服骚丝袜av| 日韩强制内射视频| 久久久午夜欧美精品| 亚洲熟女精品中文字幕| 日韩强制内射视频| 少妇猛男粗大的猛烈进出视频| 狂野欧美白嫩少妇大欣赏| 青春草视频在线免费观看| 久久国产精品大桥未久av | 一级毛片久久久久久久久女| 中文欧美无线码| 22中文网久久字幕| 自拍欧美九色日韩亚洲蝌蚪91 | 老女人水多毛片| 两个人免费观看高清视频 | 亚洲四区av| 亚洲精品乱码久久久v下载方式| 国产黄片美女视频| 最新的欧美精品一区二区| 97超视频在线观看视频| 99九九线精品视频在线观看视频| 男女国产视频网站| 亚洲四区av| 一级av片app| 自线自在国产av| av有码第一页| 免费播放大片免费观看视频在线观看| 9色porny在线观看| 精品一区二区三区视频在线| 五月天丁香电影| 极品少妇高潮喷水抽搐| 80岁老熟妇乱子伦牲交| 岛国毛片在线播放| 久久久久久久久久久久大奶| 亚洲国产精品专区欧美| 日韩成人av中文字幕在线观看| 亚洲国产色片| 国产欧美日韩一区二区三区在线 | 亚洲无线观看免费| 交换朋友夫妻互换小说| 亚洲美女视频黄频| 91在线精品国自产拍蜜月| 亚洲人与动物交配视频| 日韩一区二区视频免费看| 人妻夜夜爽99麻豆av| 欧美bdsm另类| 一级二级三级毛片免费看| 日韩免费高清中文字幕av| 麻豆精品久久久久久蜜桃| 人妻 亚洲 视频| 久久国产精品大桥未久av | 亚洲欧美一区二区三区国产| 国产高清三级在线| 久久久久久久久久久免费av| 亚洲美女黄色视频免费看| 少妇的逼好多水| 精品卡一卡二卡四卡免费| 嫩草影院入口| 自拍欧美九色日韩亚洲蝌蚪91 | 精品亚洲成a人片在线观看| 人妻一区二区av| 国产极品天堂在线| 亚洲精品国产成人久久av| 夜夜骑夜夜射夜夜干| 中文在线观看免费www的网站| 青春草视频在线免费观看| 久久女婷五月综合色啪小说| 精华霜和精华液先用哪个| 国产精品免费大片| 国产精品久久久久成人av| 三上悠亚av全集在线观看 | 国产高清国产精品国产三级| 精品少妇久久久久久888优播| 久久久久久久久大av| 精品卡一卡二卡四卡免费| 在线观看人妻少妇| 你懂的网址亚洲精品在线观看| 亚洲国产av新网站| 99久久精品国产国产毛片| 午夜精品国产一区二区电影| 亚洲在久久综合| 两个人的视频大全免费| 女人精品久久久久毛片| 国精品久久久久久国模美| 婷婷色av中文字幕| 成年av动漫网址| 不卡视频在线观看欧美| 亚洲国产av新网站| 亚洲国产精品一区三区| 99热6这里只有精品| 国产毛片在线视频| 国产 一区精品| 国产黄片视频在线免费观看| 亚洲欧美一区二区三区国产| 欧美老熟妇乱子伦牲交| 久久久亚洲精品成人影院| 亚洲av福利一区| 女性被躁到高潮视频| 久久国产精品男人的天堂亚洲 | 最近2019中文字幕mv第一页| 成年av动漫网址| 九九在线视频观看精品| 五月玫瑰六月丁香| 国产精品一区www在线观看| 中国美白少妇内射xxxbb| 国产精品福利在线免费观看| 国产高清国产精品国产三级| 99热6这里只有精品| 国产在线视频一区二区| 亚洲美女黄色视频免费看| 久久国产亚洲av麻豆专区| 99热这里只有精品一区| 国产男女超爽视频在线观看| 一区二区av电影网| 中文精品一卡2卡3卡4更新| 亚洲人与动物交配视频| 久久毛片免费看一区二区三区| 国产免费一级a男人的天堂| 久久久精品94久久精品| 少妇熟女欧美另类| kizo精华| 亚洲精品日韩av片在线观看| 久久综合国产亚洲精品| 黄色视频在线播放观看不卡| 久久精品国产亚洲av天美| 老司机影院毛片| 国产欧美另类精品又又久久亚洲欧美| 精品卡一卡二卡四卡免费| 亚洲av在线观看美女高潮| 成人亚洲精品一区在线观看| 欧美日韩视频精品一区| 久久6这里有精品| 欧美xxxx性猛交bbbb| 欧美一级a爱片免费观看看| 寂寞人妻少妇视频99o| 一本—道久久a久久精品蜜桃钙片| 亚洲,一卡二卡三卡| 国产精品伦人一区二区| 内射极品少妇av片p| freevideosex欧美| 伦理电影免费视频| 国产69精品久久久久777片| 亚洲精品一二三| 久久久久久久久久久久大奶| 免费大片18禁| 免费看av在线观看网站| 国产精品一区二区在线观看99| 黑人高潮一二区| 最近的中文字幕免费完整| 久久精品国产亚洲av涩爱| 国产一区二区在线观看日韩| 一区二区三区四区激情视频| 欧美日韩亚洲高清精品| 国产伦理片在线播放av一区| 亚洲精品日韩av片在线观看| 国产成人a∨麻豆精品| 免费少妇av软件| 精品人妻偷拍中文字幕| 青春草亚洲视频在线观看| 久久久久久久大尺度免费视频| 免费看不卡的av| 3wmmmm亚洲av在线观看| 18禁在线无遮挡免费观看视频| 日日啪夜夜爽| 免费观看av网站的网址| 久久久国产一区二区| 99热全是精品| 国产黄片美女视频| 亚洲精品色激情综合| 久久久久精品性色| 日韩一区二区视频免费看| 日韩人妻高清精品专区| 精品午夜福利在线看| 秋霞在线观看毛片| 中文资源天堂在线| 亚洲欧美精品专区久久| 中文字幕精品免费在线观看视频 | 26uuu在线亚洲综合色| av国产久精品久网站免费入址| 精品亚洲成a人片在线观看| 在线看a的网站| 国产精品免费大片| 亚洲不卡免费看| 亚洲一区二区三区欧美精品| 午夜福利视频精品| videos熟女内射| 三级国产精品欧美在线观看| 久久久国产一区二区| 欧美成人精品欧美一级黄| 日本黄大片高清| 黄色怎么调成土黄色| 久久久a久久爽久久v久久| 看免费成人av毛片| 高清av免费在线| 亚洲精品中文字幕在线视频 | 午夜久久久在线观看| 超碰97精品在线观看| 欧美少妇被猛烈插入视频| 国产精品熟女久久久久浪| 久久久久久久久久成人| 哪个播放器可以免费观看大片| 我要看黄色一级片免费的| 99re6热这里在线精品视频| 久久久国产一区二区| 午夜精品国产一区二区电影| 边亲边吃奶的免费视频| 在线观看免费高清a一片| 久久精品熟女亚洲av麻豆精品| 3wmmmm亚洲av在线观看| 免费观看的影片在线观看| 人人妻人人看人人澡| 国产一区二区在线观看av| 插阴视频在线观看视频| 欧美高清成人免费视频www| 国产精品无大码| 午夜91福利影院| 成人美女网站在线观看视频| 自拍偷自拍亚洲精品老妇| 成人18禁高潮啪啪吃奶动态图 | 毛片一级片免费看久久久久| 欧美精品高潮呻吟av久久| 男人舔奶头视频| 在线观看免费高清a一片| 久久久久久久亚洲中文字幕| 中文欧美无线码| 91精品国产国语对白视频| 在线亚洲精品国产二区图片欧美 | 久久久久久久大尺度免费视频| 伦理电影大哥的女人| a级毛片免费高清观看在线播放| 亚洲国产精品国产精品| 这个男人来自地球电影免费观看 | 欧美三级亚洲精品| 日韩中字成人| 人人妻人人添人人爽欧美一区卜| 插阴视频在线观看视频| 丝袜在线中文字幕| 欧美高清成人免费视频www| 我的老师免费观看完整版| 在线观看免费日韩欧美大片 | 日韩av在线免费看完整版不卡| 亚洲国产av新网站| 高清av免费在线| 国产老妇伦熟女老妇高清| 高清在线视频一区二区三区| 精品少妇黑人巨大在线播放| 欧美精品一区二区大全| 国产一区二区在线观看av| 国产亚洲最大av| 中文字幕精品免费在线观看视频 | 国产熟女欧美一区二区| 国产一区二区三区av在线| 午夜激情久久久久久久| 亚洲色图综合在线观看| 最近的中文字幕免费完整| 日韩视频在线欧美| av有码第一页| 国产免费福利视频在线观看| 日韩制服骚丝袜av| 啦啦啦中文免费视频观看日本| 国产精品一二三区在线看| 美女中出高潮动态图| 国产成人精品无人区| 久久久久精品久久久久真实原创| 一级毛片电影观看| 中文乱码字字幕精品一区二区三区| 久久久久久人妻| 亚洲在久久综合| 免费高清在线观看视频在线观看| 一级a做视频免费观看| 大片免费播放器 马上看| 丝袜在线中文字幕| 久久毛片免费看一区二区三区| 天堂中文最新版在线下载| 观看av在线不卡| 好男人视频免费观看在线| 国产精品蜜桃在线观看| a级片在线免费高清观看视频| 久久久久久久久久久久大奶| xxx大片免费视频| 欧美日韩国产mv在线观看视频| 女的被弄到高潮叫床怎么办| 中国美白少妇内射xxxbb| 亚洲一区二区三区欧美精品| 久久久国产一区二区| videossex国产| 内地一区二区视频在线| 黄色配什么色好看| 一级,二级,三级黄色视频| 午夜激情久久久久久久| 日韩免费高清中文字幕av| 九九爱精品视频在线观看| 精品少妇黑人巨大在线播放| 男女啪啪激烈高潮av片| √禁漫天堂资源中文www| av女优亚洲男人天堂| 蜜桃久久精品国产亚洲av| 久久鲁丝午夜福利片| 国产 一区精品| 校园人妻丝袜中文字幕| 欧美3d第一页| 人妻一区二区av| 国内少妇人妻偷人精品xxx网站| 高清av免费在线| 熟女人妻精品中文字幕| 日韩一区二区视频免费看| 交换朋友夫妻互换小说| av一本久久久久| 国产乱来视频区| 内射极品少妇av片p| 日本午夜av视频| 免费观看性生交大片5| 韩国高清视频一区二区三区| 午夜福利,免费看| 日韩大片免费观看网站| 99热全是精品| 国产成人精品无人区| 啦啦啦视频在线资源免费观看| 黄色欧美视频在线观看| 国产一区二区在线观看日韩| 亚洲av国产av综合av卡| 国产一区亚洲一区在线观看| 久久99热6这里只有精品| 九草在线视频观看| 国产伦精品一区二区三区视频9| 亚洲国产最新在线播放| 国产成人精品一,二区| 永久免费av网站大全| 男女国产视频网站| 曰老女人黄片| 尾随美女入室| 人人妻人人澡人人看| 精品卡一卡二卡四卡免费| 成年女人在线观看亚洲视频| 自拍偷自拍亚洲精品老妇| 亚洲精品国产成人久久av| 亚州av有码| 亚洲欧美日韩另类电影网站| 女性生殖器流出的白浆| 国产日韩一区二区三区精品不卡 | 久久精品久久久久久久性| 久久久久久久久久久久大奶| 一本久久精品| 国产精品一区二区三区四区免费观看| 久久影院123| 另类亚洲欧美激情| 亚洲av福利一区| 久久久久久久久久久丰满| 哪个播放器可以免费观看大片| 99热6这里只有精品| 欧美精品人与动牲交sv欧美| 免费观看性生交大片5| 亚洲久久久国产精品| 一级毛片久久久久久久久女| 久久av网站| 男女国产视频网站| kizo精华| 国产国拍精品亚洲av在线观看| 青春草视频在线免费观看| 免费不卡的大黄色大毛片视频在线观看| 国产91av在线免费观看| 国产精品福利在线免费观看| 久久久久国产精品人妻一区二区| 日韩欧美 国产精品| 老司机亚洲免费影院| 你懂的网址亚洲精品在线观看| 亚洲综合精品二区| 欧美xxⅹ黑人| av不卡在线播放| .国产精品久久| 午夜免费鲁丝| 国产欧美亚洲国产| 国精品久久久久久国模美| 亚洲欧美日韩东京热| 一区二区三区四区激情视频| 精品酒店卫生间| 卡戴珊不雅视频在线播放| 亚洲在久久综合| 亚洲情色 制服丝袜| 国产免费又黄又爽又色| 国产成人精品婷婷| 久久精品国产鲁丝片午夜精品| 99热网站在线观看| 观看免费一级毛片| 久久人人爽人人爽人人片va| 国产色婷婷99| 午夜激情久久久久久久| 国产一区有黄有色的免费视频| 久久久久精品久久久久真实原创| 香蕉精品网在线| 看非洲黑人一级黄片| 丝袜在线中文字幕| 街头女战士在线观看网站| 超碰97精品在线观看| 自拍偷自拍亚洲精品老妇| 在线观看www视频免费| 亚洲美女黄色视频免费看| 超碰97精品在线观看| 观看av在线不卡| 久久久久久久大尺度免费视频| 精华霜和精华液先用哪个| 晚上一个人看的免费电影| 精品国产乱码久久久久久小说| 久久久久久久国产电影| av福利片在线| 亚洲人成网站在线播| 精品一区二区免费观看| 丁香六月天网| 五月玫瑰六月丁香| 少妇的逼水好多| 美女福利国产在线| 国内少妇人妻偷人精品xxx网站| 女的被弄到高潮叫床怎么办| 自拍欧美九色日韩亚洲蝌蚪91 | 夫妻午夜视频| 少妇高潮的动态图| 我要看黄色一级片免费的| 三级国产精品欧美在线观看| 国产精品一区www在线观看| a级毛色黄片| 日本黄大片高清| 亚洲国产精品成人久久小说| 国产永久视频网站| 人妻少妇偷人精品九色| 男女边摸边吃奶| 少妇熟女欧美另类| 亚洲av国产av综合av卡| 最近的中文字幕免费完整| 日产精品乱码卡一卡2卡三| 久久精品夜色国产| 2022亚洲国产成人精品| 国产精品伦人一区二区| 久久精品久久久久久噜噜老黄| 国产深夜福利视频在线观看| 女人精品久久久久毛片| 久久久国产欧美日韩av| 成人国产av品久久久| 黄色配什么色好看| 成人无遮挡网站| 制服丝袜香蕉在线| 人妻少妇偷人精品九色| 亚洲精品自拍成人| 人妻系列 视频| 国产精品国产三级专区第一集| 日韩一区二区三区影片| 只有这里有精品99| 亚洲av综合色区一区| 少妇被粗大猛烈的视频| 亚洲精品国产av蜜桃| 国产一区二区在线观看日韩| 亚洲图色成人| 色吧在线观看| 精品久久久噜噜| 国产亚洲最大av| 蜜桃在线观看..| 中国国产av一级| 欧美激情极品国产一区二区三区 | av网站免费在线观看视频| 国产精品熟女久久久久浪| 亚洲欧美一区二区三区黑人 | 十八禁高潮呻吟视频 | 高清在线视频一区二区三区| 精品久久久久久电影网| 日本91视频免费播放| 久久久国产一区二区| 99久久精品一区二区三区| 97在线视频观看| 久久久欧美国产精品| 久久ye,这里只有精品| 日日摸夜夜添夜夜添av毛片| 精品少妇久久久久久888优播| 三级经典国产精品| 伦精品一区二区三区| 黄色毛片三级朝国网站 | 热re99久久国产66热| 色网站视频免费| 麻豆成人午夜福利视频| 这个男人来自地球电影免费观看 | av在线播放精品| 国产黄片美女视频| 最近手机中文字幕大全| 永久免费av网站大全| 国产日韩欧美视频二区| 国产极品天堂在线| 女的被弄到高潮叫床怎么办| 久久久久久久久大av| 欧美精品高潮呻吟av久久| 国产精品欧美亚洲77777| 国产成人一区二区在线| 91精品伊人久久大香线蕉| av福利片在线观看| 久久久久网色| av在线观看视频网站免费| 中文字幕免费在线视频6| 各种免费的搞黄视频| 国产精品国产av在线观看| 人人妻人人添人人爽欧美一区卜| 亚洲成色77777| 麻豆成人午夜福利视频| 国产精品成人在线| av黄色大香蕉| 有码 亚洲区| 免费看光身美女| 高清欧美精品videossex| 国产黄色免费在线视频| 国产欧美日韩一区二区三区在线 | 一级毛片 在线播放| 日韩一区二区视频免费看| 久久久久久久精品精品|